|
36 | 36 | k_v = 0.02 # Thermal conductivity of the vapor - W/m/K |
37 | 37 | diffVapor = 2.5e-5 # Diffusivity coefficient of the vapor - m2/s |
38 | 38 | sigBubble = 0.069 # Surface tension of the bubble - N/m |
| 39 | +mu_g = 1.48e-5 |
39 | 40 |
|
40 | 41 | # Acoustic source properties |
41 | 42 | patm = 101325. # Atmospheric pressure - Pa |
|
88 | 89 | # Simulation Algorithm Parameters ========================== |
89 | 90 | 'model_eqns' : 2, |
90 | 91 | 'time_stepper' : 3, |
91 | | - 'num_fluids' : 1, |
| 92 | + 'num_fluids' : 2, |
92 | 93 | 'num_patches' : 1, |
93 | 94 | 'viscous' : 'T', |
94 | 95 | 'mpp_lim' : 'F', |
|
143 | 144 | 'patch_icpp(1)%vel(3)' : 0., |
144 | 145 | 'patch_icpp(1)%pres' : patm/p0, |
145 | 146 | 'patch_icpp(1)%alpha_rho(1)' : rho_host/rho0, |
| 147 | + 'patch_icpp(1)%alpha_rho(2)' : 0., |
146 | 148 | 'patch_icpp(1)%alpha(1)' : 1., |
| 149 | + 'patch_icpp(1)%alpha(2)' : 0., |
147 | 150 | # ========================================================== |
148 | 151 |
|
149 | 152 | # Lagrangian Bubbles =========================== |
|
160 | 163 | 'lag_params%valmaxvoid' : 0.9, |
161 | 164 | 'lag_params%write_bubbles' : 'F', |
162 | 165 | 'lag_params%write_bubbles_stats' : 'F', |
163 | | - 'lag_params%csonhost' : c_host/c0, |
164 | | - 'lag_params%vischost' : mu_host/(rho0*x0*c0), |
165 | | - 'lag_params%Thost' : T_host/T0, |
166 | | - 'lag_params%gammagas' : gamma_g, |
167 | | - 'lag_params%gammavapor' : gamma_v, |
168 | | - 'lag_params%pvap' : pv/p0, |
169 | | - 'lag_params%cpgas' : cp_g*(T0/(c0*c0)), |
170 | | - 'lag_params%cpvapor' : cp_v*(T0/(c0*c0)), |
171 | | - 'lag_params%kgas' : k_g*(T0/(x0*rho0*c0*c0*c0)), |
172 | | - 'lag_params%kvapor' : k_v*(T0/(x0*rho0*c0*c0*c0)), |
173 | | - 'lag_params%Rgas' : (R_uni/MW_g)*(T0/(c0*c0)), |
174 | | - 'lag_params%Rvapor' : (R_uni/MW_v)*(T0/(c0*c0)), |
175 | | - 'lag_params%diffcoefvap' : diffVapor/(x0*c0), |
176 | | - 'lag_params%sigmabubble' : sigBubble/(rho0*x0*c0*c0), |
| 166 | + 'lag_params%c0' : c0, |
| 167 | + 'lag_params%rho0' : rho0, |
| 168 | + 'lag_params%T0' : T0, |
| 169 | + 'lag_params%x0' : x0, |
| 170 | + 'lag_params%diffcoefvap' : diffVapor, |
| 171 | + 'lag_params%Thost' : T_host, |
177 | 172 | # ========================================================== |
178 | 173 |
|
179 | 174 | # Fluids Physical Parameters =============================== |
| 175 | + # Host medium |
180 | 176 | 'fluid_pp(1)%gamma' : 1.0/(gamma_host-1.0), |
181 | 177 | 'fluid_pp(1)%pi_inf' : gamma_host*(pi_inf_host/p0)/(gamma_host-1.0), |
182 | 178 | 'fluid_pp(1)%Re(1)' : 1.0/(mu_host/(rho0*c0*x0)), |
| 179 | + 'fluid_pp(1)%mul0' : mu_host, |
| 180 | + 'fluid_pp(1)%ss' : sigBubble, |
| 181 | + 'fluid_pp(1)%pv' : pv, |
| 182 | + 'fluid_pp(1)%gamma_v' : gamma_v, |
| 183 | + 'fluid_pp(1)%M_v' : MW_v, |
| 184 | + 'fluid_pp(1)%k_v' : k_v, |
| 185 | + 'fluid_pp(1)%cp_v' : cp_v, |
| 186 | + |
| 187 | + # Bubble gas state |
| 188 | + 'fluid_pp(2)%gamma' : 1./(gamma_g-1.), |
| 189 | + 'fluid_pp(2)%pi_inf' : 0.0E+00, |
| 190 | + 'fluid_pp(2)%Re(1)' : 1.0/(mu_g/(rho0*c0*x0)), |
| 191 | + 'fluid_pp(2)%gamma_v' : gamma_g, |
| 192 | + 'fluid_pp(2)%M_v' : MW_g, |
| 193 | + 'fluid_pp(2)%k_v' : k_g, |
| 194 | + 'fluid_pp(2)%cp_v' : cp_g, |
183 | 195 | # ========================================================== |
184 | 196 | })) |
185 | 197 |
|
|
0 commit comments