~fluidity-core/fluidity/exorcised

« back to all changes in this revision

Viewing changes to multiphase/tests/BL_test_10elements_1d_compressible/input.dat

  • Committer: saunde01
  • Date: 2011-03-23 13:15:31 UTC
  • Revision ID: svn-v4:5bf5533e-7014-46e3-b1bb-cce4b9d03719:trunk:3310
As discussed in the Dev meeting and in emails, this commit renames the multiphase directory to legacy_reservoir_prototype, and adds an obvious warning to the top of the README file to warn unsuspecting users of the experimental nature of the code. Iterations on the exact wording very welcome

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
###
2
 
###              BL test COMPRESSIBLE
3
 
### First part - before first *** exit *** statement
4
 
### is for scalar variables that will set up the right 
5
 
### size of the arrays
6
 
###
7
 
# Problem type: options -2, -1, 0, and 1
8
 
problem               1
9
 
# Number of phases
10
 
nphase                2
11
 
# total number of elements
12
 
totele               10
13
 
# Dimensionality
14
 
ndim                  1
15
 
# nlev --> controls u_snloc
16
 
nlev                  3
17
 
# Number of velocity nodes
18
 
u_nloc                6
19
 
# Number of spatial nodes associated w velocity 
20
 
xu_nloc               2
21
 
# number of nodes associated with CV
22
 
cv_nloc               3
23
 
# Number of spatial nodes
24
 
x_nloc                3
25
 
# Number of pressure nodes
26
 
p_nloc                3
27
 
# Number of surface nodes associated with CV
28
 
cv_snloc              1
29
 
# Number of surface nodes associated w velocity
30
 
u_snloc              -1
31
 
# Number of surface pressure nodes
32
 
p_snloc               1
33
 
# Number of surface spatial nodes
34
 
x_snloc               1
35
 
# Total number of surface elements
36
 
stotel                2
37
 
# Number of coeffs for polynomial representation for EOS
38
 
ncoef                10
39
 
# Number of coeffs for polynomial representation of permeability/absorption
40
 
nuabs_coefs           1
41
 
# Option for element type: Velocity mesh
42
 
u_ele_type            2
43
 
# Option for element type: Pressure mesh
44
 
p_ele_type            2
45
 
# Option for element type: Material mesh
46
 
mat_ele_type          1
47
 
# Option for element type: CV mesh
48
 
cv_ele_type           2
49
 
# Option for element type: CV surface mesh
50
 
cv_sele_type          1
51
 
# Option for element type: Velocity surface mesh
52
 
u_sele_type           1
53
 
# Total number of time dumps
54
 
ntime               100
55
 
# Total number of non-linear iterations
56
 
nits                  3
57
 
# Option for modified CMC color algorithm(/=0), if <0==>ndpset=cv_nonods 
58
 
ndpset                0
59
 
# Time-step size
60
 
dt                    1.0e-3
61
 
# Reference pressure (not in use)
62
 
patmos                0.
63
 
# Initial pressure (not in use)
64
 
p_ini                 0.
65
 
# Initial temperature (or any scalar field, not in use)
66
 
t_ini                 0.
67
 
# Conservative(1.)/non-conservative(0.) flag for CV discretisation
68
 
t_beta                0.
69
 
# Conservative(1.)/non-conservative(0.) flag for Volfra discretisation
70
 
v_beta                1.
71
 
# Time-stepping discretisation parameter
72
 
t_theta               0.
73
 
# Time-stepping discretisation parameter
74
 
v_theta               1.
75
 
# Time-stepping discretisation parameter (not in use)
76
 
u_theta               1.
77
 
# Disopt: discretisation option in space and time (see cv-adv-dif.f90)
78
 
# = 8: Finite elements in space (Theta = specified;  Limiting: Downwind+)
79
 
t_disopt              1
80
 
# Disopt: discretisation option in space and time (see cv-adv-dif.f90)
81
 
# = 0: 1st order in space (Theta=specified;  Limiting: universal)
82
 
u_disopt              1
83
 
# Disopt: discretisation option in space and time (see cv-adv-dif.f90)
84
 
# = 0: 1st order in space (Theta=specified;  Limiting: universal)
85
 
v_disopt              8
86
 
#  t_dg_vel_int_opt: interface scalar field calculation option between elements
87
 
t_dg_vel_int_opt      0
88
 
#  u_dg_vel_int_opt: interface velocity calculation option between elements
89
 
u_dg_vel_int_opt      4
90
 
#  v_dg_vel_int_opt: interface velocitycalculation option between elements
91
 
v_dg_vel_int_opt      4
92
 
#  w_dg_vel_int_opt: interface velocitycalculation option between elements
93
 
w_dg_vel_int_opt      0
94
 
# Lump_eqns: Lump multiphase flow equations
95
 
lump_eqns             F
96
 
# Domain length
97
 
domain_length         1.
98
 
###
99
 
### end of scalars variables necessary for the arrays read after
100
 
###
101
 
exit    10
102
 
 
103
 
                                                                              
104
 
# BC variables: 
105
 
# If value < -1000, then a pre-defined function may be used.
106
 
 
107
 
# wic_vol_bc( stotel * nphase )
108
 
wic_vol_bc    -1001 input_wic_vol_bc_BL1_fcn
109
 
 
110
 
# wic_d_bc( stotel * nphase )
111
 
wic_d_bc    -1001 input_wic_d_bc_BL1_fcn
112
 
 
113
 
# wic_u_bc( stotel * nphase )
114
 
wic_u_bc    -1001 input_wic_u_bc_BL1_fcn
115
 
 
116
 
# wic_p_bc( stotel * nphase )
117
 
wic_p_bc    -1001 input_wic_p_bc_BL1_fcn
118
 
 
119
 
# wic_t_bc( stotel * nphase )
120
 
wic_t_bc  0   
121
 
 
122
 
# suf_vol_bc( stotel * cv_snloc * nphase )
123
 
suf_vol_bc    -1001.  input_suf_vol_bc_BL1_fcn
124
 
 
125
 
# suf_d_bc( stotel * cv_snloc * nphase )
126
 
suf_d_bc    -1001.  input_suf_d_bc_BL1_fcn
127
 
 
128
 
# suf_cpd_bc( stotel * cv_snloc * nphase )
129
 
suf_cpd_bc    0.
130
 
 
131
 
# suf_t_bc( stotel * cv_snloc * nphase )
132
 
suf_t_bc   -1001 input_suf_t_bc_BL1_fcn
133
 
 
134
 
# suf_p_bc ( stotel * p_snloc * nphase )
135
 
suf_p_bc       0.
136
 
 
137
 
# suf_u_bc( stotel * u_snloc * nphase )
138
 
suf_u_bc      -1001.  input_suf_u_bc_BL1_fcn
139
 
 
140
 
# suf_v_bc( stotel * u_snloc * nphase )
141
 
suf_v_bc       0.
142
 
  
143
 
# suf_w_bc( stotel * u_snloc * nphase )
144
 
suf_w_bc       0.
145
 
  
146
 
# suf_one_bc( stotel * cv_snloc * nphase )
147
 
suf_one_bc    0.
148
 
 
149
 
# suf_u_bc_rob1( stotel * u_snloc * nphase )
150
 
suf_u_bc_rob1  0.
151
 
 
152
 
# suf_u_bc_rob2( stotel * u_snloc * nphase )
153
 
suf_u_bc_rob2  0.
154
 
 
155
 
# suf_v_bc_rob1( stotel * u_snloc * nphase )
156
 
suf_v_bc_rob1  0.
157
 
 
158
 
# suf_v_bc_rob2( stotel * u_snloc * nphase )
159
 
suf_v_bc_rob2  0.
160
 
 
161
 
# suf_w_bc_rob1( stotel * u_snloc * nphase )
162
 
suf_w_bc_rob1  0.
163
 
 
164
 
# suf_w_bc_rob2( stotel * u_snloc * nphase )
165
 
suf_w_bc_rob2  0.
166
 
 
167
 
# suf_t_bc_rob1( stotel * cv_snloc * nphase )
168
 
suf_t_bc_rob1  0.
169
 
  
170
 
# suf_t_bc_rob2( stotel * cv_snloc * nphase )
171
 
suf_t_bc_rob1  0.
172
 
  
173
 
# x( x_nonods )
174
 
x    0.
175
 
 
176
 
# y( x_nonods )
177
 
y    0.
178
 
 
179
 
# z( x_nonods )
180
 
z    0.
181
 
 
182
 
# xu( xu_nonods )
183
 
xu    0.
184
 
 
185
 
# yu( xu_nonods )
186
 
yu    0.
187
 
 
188
 
# zu( xu_nonods )
189
 
zu    0.
190
 
 
191
 
# nu( u_nonods * nphase )
192
 
nu    1.
193
 
 
194
 
# nv( u_nonods * nphase )
195
 
nv    0.
196
 
 
197
 
# nw( u_nonods * nphase )
198
 
nw    0.
199
 
 
200
 
# ug( u_nonods * nphase )
201
 
ug    0.
202
 
 
203
 
# vg( u_nonods * nphase )
204
 
vg    0.
205
 
 
206
 
# wg( u_nonods * nphase )
207
 
wg    0.
208
 
 
209
 
# uabs_option( nphase )
210
 
uabs_option  3
211
 
 
212
 
# uabs_coefs( nphase, nuabs_coefs )
213
 
uabs_coefs   1.
214
 
 
215
 
# u_abs_stab ( mat_nonods, ndim * nphase, ndim * nphase )
216
 
u_abs_stab    0.
217
 
 
218
 
# u_absorb ( mat_nonods, ndim * nphase, ndim * nphase )
219
 
u_absorb   0.
220
 
 
221
 
# t_absorb ( cv_pha_nonods, nphase, nphase )
222
 
t_absorb    0.
223
 
 
224
 
# v_absorb ( cv_pha_nonods, nphase, nphase )
225
 
v_absorb  0.
226
 
 
227
 
# u_source( u_pha_nonods )
228
 
u_source  0.
229
 
 
230
 
# t_source( cv_pha_nonods )
231
 
t_source   0.
232
 
 
233
 
# v_source( cv_pha_nonods )
234
 
v_source   0.
235
 
 
236
 
#  udiffusion( mat_nonods, ndim, ndim, nphase )
237
 
udiffusion   0.
238
 
 
239
 
#  tdiffusion( mat_nonods, ndim, ndim, nphase )
240
 
tdiffusion   0.
241
 
 
242
 
# u ( u_pha_nonods )
243
 
u    0.
244
 
 
245
 
# v ( u_pha_nonods )
246
 
v    0.
247
 
 
248
 
# w ( u_pha_nonods )
249
 
w     0.
250
 
 
251
 
# den ( cv_pha_nonods )
252
 
den    1.
253
 
 
254
 
# satura ( cv_pha_nonods )
255
 
satura   -1001.  input_satura_BL1_func
256
 
 
257
 
# volfra ( cv_pha_nonods )
258
 
volfra    0.
259
 
 
260
 
# t ( cv_pha_nonods )
261
 
t    0.
262
 
 
263
 
# cv_one( cv_pha_nonods )
264
 
cv_one    0.
265
 
 
266
 
# p ( cv_nonods )
267
 
p     0.
268
 
 
269
 
# cv_p (cv_nonods )
270
 
cv_p    0.
271
 
 
272
 
# volfra_pore ( totele )
273
 
volfra_pore   .5
274
 
 
275
 
# perm ( totele, ndim, ndim )
276
 
perm      -1001.   input_perm_BL1_func
277
 
 
278
 
# eos_option( nphase )
279
 
eos_option   2
280
 
 
281
 
# cp_option( nphase )
282
 
cp_option   0
283
 
 
284
 
# eos_coefs( nphase, ncoef )
285
 
eos_coefs   -1001.  input_eos_coef_BL1_func
286
 
 
287
 
# cp_coefs( nphase, ncp_coefs )
288
 
cp_coefs   1.
289
 
 
290
 
# end of the input file
291
 
exit  10
292
 
 
293