blob: 4af00bc8f3a8281432283acd549d9bd29d9842fe (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
|
****************************************************************
* *
* Procedure : ThmComponent.c2m *
* Purpose : Thermo-hydraulics component *
* Author : A. Hebert *
* *
* CALL : *
* Fmap Thm := ThmComponent Fmap Thm *
* :: <<iter>> <<densB>> >>CONV<< ; *
* *
****************************************************************
PARAMETER Fmap Thm ::
::: LINKED_LIST Fmap Thm ; ;
MODULE THM: GREP: END: ;
INTEGER iter ;
:: >>iter<< ;
DOUBLE DdensB ;
:: >>DdensB<< ;
REAL densB := DdensB D_TO_R ;
REAL erra1 ;
LOGICAL CONV ;
IF iter 1 = THEN
Thm Fmap := THM: Fmap ::
EDIT 2
TIME 0.0 (* s *)
CONV 50 50 50 0.5 1.0E-3
CRITFL 1.80E+6 (* W/m^2 *)
ASSMB 264 25
CWSECT 3.86 (* m^2 *) 67398.0 (* m^3/hr *)
INLET 15500000.0 (* Pa *) 560.95 (* K *)
RADIUS 0.4095E-2 0.418E-2 0.475E-2 0.6025E-2 (* m *)
SET-PARAM 'C-BORE' <<densB>>
RODMESH 5 10
CONV 50 50 50 1.0 1.0E-4
RELAX 0.5
;
ELSE
Thm Fmap := THM: Thm Fmap ;
ENDIF ;
*--
* Convergence test
*--
GREP: Thm :: GETVAL 'ERROR-T-FUEL' 1 >>erra1<< ;
ECHO "error on fuel temperature at iteration" iter "=" erra1 ;
EVALUATE CONV := erra1 1.0E-1 < iter 50 > + ;
ECHO "at iter=" iter "erra1=" erra1 " CONV=" CONV ;
ECHO "-------------------------------" ;
:: <<CONV>> ;
END: ;
|
