forked from dshaver14/wirehexmesher
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathbase.rea
180 lines (180 loc) · 6.37 KB
/
base.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
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
****** PARAMETERS *****
2.610000 NEKTON VERSION
3 DIMENSIONAL RUN
103 PARAMETERS FOLLOW
1.000000 p01 DENSITY
-100.0000 p02 VISCOS
0.0000000E+00
0.0000000E+00
0.0000000E+00
0.0000000E+00
1.000000 p07 RHOCP
1.000000 p08 CONDUCT
0.0000000E+00
0.0000000E+00 p10 FINTIME
1.000 p11 NSTEPS
-5.0000000E-03 p12 DT
0.0000000E+00 p13 IOCOMM
1.0000000E+00 p14 IOTIME
0.0000000E+00 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-10 p24 TOLREL
1.0000000E-10 p25 TOLABS
0.2500000 p26 COURANT/NTAU
2.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
0.0000000E+00 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.0000000E+00
0.0000000E+00
0.0000000E+00
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
0.0000000E+00 p93 Numbr of prev pressure solns saved
0.0000000E+00 p94 start projecting vel. after p94 step
0.0000000E+00 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.0000001E-02 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 T 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
8.00000 8.00000 -4.000000 -4.00000 XFAC,YFAC,XZERO,YZERO
**MESH DATA** 1st line is X of corner 1,2,3,4. 2nd line is Y.
1 3 1 NEL,NDIM,NELV
ELEMENT 1 [ 1a] GROUP 0
0.0000000E+00 1.0000000E+00 1.0000000E+00 0.0000000E+00
0.0000000E+00 0.0000000E+00 1.0000000E+00 1.0000000E+00
0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
0.0000000E+00 1.0000000E+00 1.0000000E+00 0.0000000E+00
0.0000000E+00 0.0000000E+00 1.0000000E+00 1.0000000E+00
1.0000000E+00 1.0000000E+00 1.0000000E+00 1.0000000E+00
***** CURVED SIDE DATA *****
0 Curved sides follow IEDGE,IEL,CURVE(I),I=1,5, CCURVE
***** BOUNDARY CONDITIONS *****
***** FLUID BOUNDARY CONDITIONS *****
1 1 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 2 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 3 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 4 0.000000 0.000000 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 4 0.000000 0.000000 0.0000000E+00 0.0000000E+00 0.0000000E+00
1 4 0.000000 0.000000 0.0000000E+00 0.0000000E+00 0.0000000E+00
***** NO THERMAL BOUNDARY CONDITIONS *****
0 PRESOLVE/RESTART OPTIONS *****
7 INITIAL CONDITIONS *****
C Default
C Default
C Default
C Default
C Default
C Default
C Default
***** DRIVE FORCE DATA ***** BODY FORCE, FLOW, Q
4 Lines of Drive force data follow
C
C
C
C
***** Variable Property Data ***** Overrrides Parameter data.
1 Lines follow.
0 PACKETS OF DATA FOLLOW
***** HISTORY AND INTEGRAL DATA *****
0 POINTS. Hcode, I,J,H,IEL
***** OUTPUT FIELD SPECIFICATION *****
6 SPECIFICATIONS FOLLOW
T COORDINATES
F VELOCITY
F PRESSURE
T TEMPERATURE
F TEMPERATURE GRADIENT
0 PASSIVE SCALARS
***** OBJECT SPECIFICATION *****
0 Surface Objects
0 Volume Objects
0 Edge Objects
0 Point Objects