@@ -31,7 +31,7 @@ k_n = 0.02556
3131gamma_gas = 1.4
3232
3333#reference bubble size
34- R0ref = 10 .E-06
34+ R0ref = 1 .E-06
3535
3636pa = 0.1 * 1.E+06 / 101325.
3737
@@ -44,7 +44,8 @@ Ca = 1.
4444# We = rho0*(uu**2.)*R0ref/ss
4545We = 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
5051vf0 = XXVF
@@ -56,8 +57,8 @@ t0 = x0/c0
5657nbubbles = 1
5758myr0 = R0ref
5859
59- cfl = 0.1
60- Nx = 400
60+ cfl = XXCFL
61+ Nx = XXNX
6162Ldomain = 20.E-03
6263L = Ldomain / x0
6364dx = L / float (Nx )
@@ -67,16 +68,14 @@ Tpulse = Lpulse/cact
6768Tfinal = 0.25 * 10. * Tpulse * c0 / x0
6869Nt = 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 .
7676Nout = int (math .ceil (Nt / Nfiles ))
7777Nt = int (Nout * Nfiles )
7878
79-
8079print ('Web' , We )
8180print ('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