Skip to content

Commit 1f764b5

Browse files
sbryngelsonsbryngelson
committed
Merge branch 'qbmm' of github.com:ComputationalFlowPhysics/MFC_private into qbmm
2 parents bc2bdd1 + 27fde1b commit 1f764b5

File tree

22 files changed

+2550
-239
lines changed

22 files changed

+2550
-239
lines changed

‎.gitignore‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -21,12 +21,14 @@ tests/*/silo*
2121
tests/*/restart*
2222
tests/*/*.out
2323
tests/*/binary
24+
example_cases/*batch/*/
2425
example_cases/*/D/*
2526
example_cases/*/p*
2627
example_cases/*/D_*
2728
example_cases/*/*.inf
2829
example_cases/*/*.inp
2930
example_cases/*/*.dat
31+
example_cases/*/*.o*
3032
example_cases/*/silo*
3133
example_cases/*/restart*
3234
example_cases/*/*.out

‎example_cases/1D_qbmm/input_qbmm.py‎

Lines changed: 16 additions & 17 deletions
Original file line numberDiff line numberDiff line change
@@ -133,25 +133,25 @@
133133
# 't_step_stop' : 4, \
134134
't_step_save' : 8000, \
135135
# 't_step_save' : 1, \
136-
# ==========================================================
136+
# ==========================================================
137137
\
138138
# Simulation Algorithm Parameters ==========================
139139
'num_patches' : 2, \
140140
'model_eqns' : 2, \
141141
'alt_soundspeed' : 'F', \
142142
'num_fluids' : 1, \
143-
'adv_alphan' : 'T', \
144-
'mpp_lim' : 'F', \
145-
'mixture_err' : 'F', \
146-
'time_stepper' : 3, \
143+
'adv_alphan' : 'T', \
144+
'mpp_lim' : 'F', \
145+
'mixture_err' : 'F', \
146+
'time_stepper' : 3, \
147147
'weno_vars' : 2, \
148148
'weno_order' : 5, \
149149
'weno_eps' : 1.E-16, \
150150
'char_decomp' : 'F', \
151151
'mapped_weno' : 'T', \
152152
'null_weights' : 'F', \
153153
'mp_weno' : 'T', \
154-
'riemann_solver' : 2, \
154+
'riemann_solver' : 2, \
155155
'wave_speeds' : 1, \
156156
'avg_state' : 2, \
157157
'commute_err' : 'F', \
@@ -164,13 +164,13 @@
164164
'format' : 1, \
165165
'precision' : 2, \
166166
'prim_vars_wrt' :'T', \
167-
'parallel_io' :'F', \
168-
'fd_order' : 1, \
167+
'parallel_io' :'F', \
168+
'fd_order' : 1, \
169169
#'schlieren_wrt' :'T', \
170-
'probe_wrt' :'T', \
171-
'num_probes' : 1, \
172-
'probe(1)%x' : 0., \
173-
# ==========================================================
170+
'probe_wrt' :'T', \
171+
'num_probes' : 1, \
172+
'probe(1)%x' : 0., \
173+
# ==========================================================
174174

175175
# Patch 1 _ Background =====================================
176176
'patch_icpp(1)%geometry' : 1, \
@@ -229,14 +229,13 @@
229229
'bubbles' : 'T', \
230230
'bubble_model' : 2, \
231231
'polytropic' : 'T', \
232-
# 'polydisperse' : 'T', \
233-
# 'R0_type' : 2, \
234-
'polydisperse' : 'F', \
235-
# 'poly_sigma' : 0.3, \
232+
'polydisperse' : 'T', \
233+
'R0_type' : 1, \
234+
'poly_sigma' : 0.3, \
236235
'thermal' : 3, \
237236
'R0ref' : myr0, \
238237
# 'nb' : 3, \
239-
'nb' : 1, \
238+
'nb' : 3, \
240239
'Ca' : Ca, \
241240
# 'Web' : We, \
242241
# 'Re_inv' : Re_inv, \

‎example_cases/1D_qbmm_batch/dummy/input_qbmm.dummy‎

Lines changed: 34 additions & 36 deletions
Original file line numberDiff line numberDiff line change
@@ -31,7 +31,7 @@ k_n = 0.02556
3131
gamma_gas = 1.4
3232

3333
#reference bubble size
34-
R0ref = 10.E-06
34+
R0ref = 1.E-06
3535

3636
pa = 0.1 * 1.E+06 / 101325.
3737

@@ -44,7 +44,8 @@ Ca = 1.
4444
# We = rho0*(uu**2.)*R0ref/ss
4545
We = p0*R0ref/ss
4646
#Inv. bubble Reynolds number
47-
Re_inv = mul0/(rho0*uu*R0ref)
47+
# Re_inv = mul0/(rho0*uu*R0ref)
48+
Re_inv = 1./float(XXRE)
4849

4950
#IC setup
5051
vf0 = XXVF
@@ -56,8 +57,8 @@ t0 = x0/c0
5657
nbubbles = 1
5758
myr0 = R0ref
5859

59-
cfl = 0.1
60-
Nx = 400
60+
cfl = XXCFL
61+
Nx = XXNX
6162
Ldomain = 20.E-03
6263
L = Ldomain/x0
6364
dx = L/float(Nx)
@@ -67,16 +68,14 @@ Tpulse = Lpulse/cact
6768
Tfinal = 0.25*10.*Tpulse*c0/x0
6869
Nt = int(Tfinal/dt)
6970

70-
dt = dt * 0.1
71+
# dt = dt * 0.1
7172
# print('dt: ',dt)
7273

73-
Nt = 10000
74-
75-
Nfiles = 2.
74+
# Nt = 10000
75+
Nfiles = XXNFILES.
7676
Nout = int(math.ceil(Nt/Nfiles))
7777
Nt = int(Nout*Nfiles)
7878

79-
8079
print('Web', We)
8180
print('Re_inv', Re_inv)
8281

@@ -118,7 +117,7 @@ case_dict =
118117
'case_dir' : '\'.\'', \
119118
'run_time_info' : 'F', \
120119
'nodes' : 1, \
121-
'ppn' : 12, \
120+
'ppn' : XXNP, \
122121
'queue' : 'normal', \
123122
'walltime' : '24:00:00', \
124123
'mail_list' : '', \
@@ -132,29 +131,29 @@ case_dict =
132131
'm' : Nx, \
133132
'n' : 0, \
134133
'p' : 0, \
135-
'dt' : dt, \
134+
'dt' : dt, \
136135
't_step_start' : 0, \
137136
't_step_stop' : Nt, \
138-
't_step_save' : Nout, \
139-
# ==========================================================
137+
't_step_save' : Nout, \
138+
# ==========================================================
140139
\
141140
# Simulation Algorithm Parameters ==========================
142141
'num_patches' : 2, \
143142
'model_eqns' : 2, \
144143
'alt_soundspeed' : 'F', \
145144
'num_fluids' : 1, \
146-
'adv_alphan' : 'T', \
147-
'mpp_lim' : 'F', \
148-
'mixture_err' : 'F', \
149-
'time_stepper' : 3, \
145+
'adv_alphan' : 'T', \
146+
'mpp_lim' : 'F', \
147+
'mixture_err' : 'F', \
148+
'time_stepper' : 3, \
150149
'weno_vars' : 2, \
151150
'weno_order' : 5, \
152151
'weno_eps' : 1.E-16, \
153152
'char_decomp' : 'F', \
154153
'mapped_weno' : 'T', \
155154
'null_weights' : 'F', \
156155
'mp_weno' : 'T', \
157-
'riemann_solver' : 2, \
156+
'riemann_solver' : 2, \
158157
'wave_speeds' : 1, \
159158
'avg_state' : 2, \
160159
'commute_err' : 'F', \
@@ -167,13 +166,13 @@ case_dict =
167166
'format' : 1, \
168167
'precision' : 2, \
169168
'prim_vars_wrt' :'T', \
170-
'parallel_io' :'F', \
171-
'fd_order' : 1, \
169+
'parallel_io' :'T', \
170+
'fd_order' : 1, \
172171
#'schlieren_wrt' :'T', \
173-
'probe_wrt' :'T', \
174-
'num_probes' : 1, \
175-
'probe(1)%x' : 0., \
176-
# ==========================================================
172+
'probe_wrt' :'T', \
173+
'num_probes' : 1, \
174+
'probe(1)%x' : 0., \
175+
# ==========================================================
177176

178177
# Patch 1 _ Background =====================================
179178
'patch_icpp(1)%geometry' : 1, \
@@ -230,24 +229,23 @@ case_dict =
230229

231230
# Bubbles ==================================================
232231
'bubbles' : 'T', \
233-
'bubble_model' : 2, \
232+
'bubble_model' : 3, \
234233
'polytropic' : 'T', \
235-
# 'polydisperse' : 'T', \
236-
# 'R0_type' : 2, \
237-
'polydisperse' : 'F', \
238-
# 'poly_sigma' : 0.3, \
234+
'polydisperse' : 'T', \
235+
'R0_type' : XXR0TYPE, \
236+
'polydisperse' : 'XXPOLY', \
237+
'poly_sigma' : XXSIG0, \
239238
'thermal' : 3, \
240239
'R0ref' : myr0, \
241-
# 'nb' : 3, \
242-
'nb' : 1, \
240+
'nb' : XXNB, \
243241
'Ca' : Ca, \
244-
# 'Web' : We, \
245-
# 'Re_inv' : Re_inv, \
242+
'Web' : We, \
243+
'Re_inv' : Re_inv, \
246244
'qbmm' : 'T', \
247245
'nnode' : 4, \
248246
'dist_type' : 2, \
249-
'sigR' : 0.1, \
250-
'sigV' : 0.1, \
247+
'sigR' : XXSIGR, \
248+
'sigV' : XXSIGV, \
251249
'rhoRV' : 0.0, \
252250
# ==========================================================
253251

@@ -258,7 +256,7 @@ case_dict =
258256
'Mono(1)%npulse' : 1, \
259257
'Mono(1)%dir' : 1., \
260258
'Mono(1)%pulse' : 1, \
261-
'Mono(1)%mag' : 1*pa, \
259+
'Mono(1)%mag' : XXPFAC*pa, \
262260
'Mono(1)%length' : (1./(300000.))*cact/x0, \
263261
# ==========================================================
264262
}

0 commit comments

Comments
 (0)