-
Notifications
You must be signed in to change notification settings - Fork 5
/
Copy pathtop.rea
128 lines (128 loc) · 4.51 KB
/
top.rea
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
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
****** PARAMETERS *****
2.610000 NEKTON VERSION
3 DIMENSIONAL RUN
103 PARAMETERS FOLLOW
1.000000 p01 DENSITY
-1400.000 p02 VISCOS
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00 p07 RHOCP
-4000.000 p08 CONDUCT
0.0000000E+00
0.0000000E+00 p10 FINTIME
4000.000 p11 NSTEPS
-5.0000000E-03 p12 DT
0.0000000E+00 p13 IOCOMM
0.0000000E+00 p14 IOTIME
100.0000 p15 IOSTEP
0.0000000E+00 p16 PSSOLVER
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
1.0000000E-06 p21 DIVERGENCE
1.0000000E-08 p22 HELMHOLTZ
0.0000000E+00 p23 NPSCAL
1.0000000E-04 p24 TOLREL
1.0000000E-05 p25 TOLABS
1.000000 p26 COURANT/NTAU
3.000000 p27 TORDER
0.0000000E+00 p28 TORDER: mesh velocity (0: p28=p27)
0.0000000E+00 p29 magnetic visc if > 0, = -1/Rm if < 0
0.0000000E+00 p30 > 0 ==> properties set in uservp()
0.0000000E+00 p31 NPERT: #perturbation modes
0.0000000E+00 p32 #BCs in re2 file, if > 0
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00 p42 0=gmres/1=pcg
0.0000000E+00 p43 0=semg/1=schwarz
0.0000000E+00 p44 0=E-based/1=A-based prec.
0.0000000E+00 p45 Relaxation factor for DTFS
0.0000000E+00 p46 if >0, dont call SETICS
0.0000000E+00 p47 vnu: mesh matieral prop
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00 p52 IOHIS
0.0000000E+00
0.0000000E+00 p54 1,2,3-->fixed flow rate dir=x,y,z
0.0000000E+00 p55 vol.flow rate(p54>0) or Ubar(p54<0)
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00 p59 !=0 --> use std axhelm for all elem
0.0000000E+00 p60 !=0--> init velocity to small nonzero
0.0000000E+00
0.0000000E+00 p62 >0 --> force byte_swap for output
0.0000000E+00 p63 =8 --> force 8-byte output
0.0000000E+00 p64 =1 --> perturbation restart
1.000000 p65 #iofile(eg,0 or 64); <0 --> sep. dirs
6.000000 p66 write fmt:ONLY postx uses rea value
6.000000 p67 read fmt: same modes as p66
0.0000000E+00 p68 iastep: freq for avg_all
0.0000000E+00
0.0000000E+00
0.580000 p71 = delta_99
1.000000E+00 p72 timestep for particle data dump
12000 p73 Total # of particles at any instant
0.0000000E+00 p74 verbose Helmholtz
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00 p84 !=0 -->sets initial timestep if p12>0
0.0000000E+00 p85 dt ratio if p84 !=0, for timesteps>0
0.0000000E+00 p86 reserved
0.0000000E+00
0.0000000E+00
0.0000000E+00 p89 reserved
0.0000000E+00
0.0000000E+00
0.0000000E+00
20.0000 p93 Numbr of prev pressure solns saved
5.000000 p94 start projecting vel. after p94 step
5.000000 p95 start projecting pr after p95 step
0.0000000E+00
0.0000000E+00
0.0000000E+00
3.000000 p99 dealiasing:if <0 disable
0.0000000E+00 p100 reserved
0.0000000E+00 p101 No. of additional filter modes
0.0000000E+00
5.0000000E-03 p103 weight of stabilizing filter (.01)
4 Lines of passive scalar data follows2 CONDUCT; 2RHOCP
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000 1.00000
1.00000 1.00000 1.00000 1.00000
13 LOGICAL SWITCHES FOLLOW
T IFFLOW
F IFHEAT
T IFTRAN
T F F F F F F F F F F IFNAV & IFADVC (convection in P.S. fields)
F F T T T T T T T T T T IFTMSH (IF mesh for this field is T mesh)
F IFAXIS
F IFSTRS
F IFSPLIT
F IFMGRID
F IFMODEL
F IFKEPS
F IFMVBD
F IFCHAR
25.0000 25.0000 -5.00000 -15.0000 XFAC,YFAC,XZERO,YZERO
**MESH DATA** 6 lines are X,Y,Z;X,Y,Z. Columns corners 1-4;5-8