Skip to content

Commit d4f382a

Browse files
hyeoksu-leehyeoksu-lee
committed
fix swapped vapor gas properties
1 parent 2902000 commit d4f382a

File tree

20 files changed

+544
-789
lines changed

20 files changed

+544
-789
lines changed

‎examples/2D_lagrange_bubblescreen/case.py‎

Lines changed: 24 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -26,8 +26,8 @@
2626
R_uni = 8314 # Universal gas constant - J/kmol/K
2727
MW_g = 28.0 # Molar weight of the gas - kg/kmol
2828
MW_v = 18.0 # Molar weight of the vapor - kg/kmol
29-
gamma_g = 1.4 # Specific heat ratio of the gas
30-
gamma_v = 1.333 # Specific heat ratio of the vapor
29+
gam_g = 1.4 # Specific heat ratio of the gas
30+
gam_v = 1.333 # Specific heat ratio of the vapor
3131
pv = 2350 # Vapor pressure of the host - Pa
3232
cp_g = 1.0e3 # Specific heat of the gas - J/kg/K
3333
cp_v = 2.1e3 # Specific heat of the vapor - J/kg/K
@@ -128,6 +128,8 @@
128128
# Lagrangian Bubbles
129129
"bubbles_lagrange": "T",
130130
"bubble_model": 2, # Keller-Miksis model
131+
"thermal": 3,
132+
"polytropic": "F",
131133
"lag_params%nBubs_glb": 1194, # Number of bubbles
132134
"lag_params%solver_approach": 2,
133135
"lag_params%cluster_type": 2,
@@ -139,33 +141,34 @@
139141
"lag_params%valmaxvoid": 0.9,
140142
"lag_params%write_bubbles": "F",
141143
"lag_params%write_bubbles_stats": "F",
142-
"lag_params%c0": c0,
143-
"lag_params%rho0": rho0,
144-
"lag_params%T0": T0,
145-
"lag_params%x0": x0,
146-
"lag_params%diffcoefvap": diffVapor,
147-
"lag_params%Thost": T_host,
148-
"lag_params%charwidth": z_virtual / x0,
144+
# Bubble parameters
145+
"bub_pp%R0ref": 1.0,
146+
"bub_pp%p0ref": 1.0,
147+
"bub_pp%rho0ref": 1.0,
148+
"bub_pp%T0ref": 1.0,
149+
"bub_pp%ss": sigBubble / (rho0 * x0 * c0 * c0),
150+
"bub_pp%pv": pv / p0,
151+
"bub_pp%vd": diffVapor / (x0 * c0),
152+
"bub_pp%mu_l": mu_host / (rho0 * x0 * c0),
153+
"bub_pp%gam_v": gam_v,
154+
"bub_pp%gam_g": gam_g,
155+
"bub_pp%M_v": MW_v,
156+
"bub_pp%M_g": MW_g,
157+
"bub_pp%k_v": k_v * (T0 / (x0 * rho0 * c0 * c0 * c0)),
158+
"bub_pp%k_g": k_g * (T0 / (x0 * rho0 * c0 * c0 * c0)),
159+
"bub_pp%cp_v": cp_v * (T0 / (c0 * c0)),
160+
"bub_pp%cp_g": cp_g * (T0 / (c0 * c0)),
161+
"bub_pp%R_v": (R_uni / MW_v) * (T0 / (c0 * c0)),
162+
"bub_pp%R_g": (R_uni / MW_g) * (T0 / (c0 * c0)),
149163
# Fluids Physical Parameters
150164
# Host medium
151165
"fluid_pp(1)%gamma": 1.0 / (gamma_host - 1.0),
152166
"fluid_pp(1)%pi_inf": gamma_host * (pi_inf_host / p0) / (gamma_host - 1.0),
153167
"fluid_pp(1)%Re(1)": 1.0 / (mu_host / (rho0 * c0 * x0)),
154-
"fluid_pp(1)%mul0": mu_host,
155-
"fluid_pp(1)%ss": sigBubble,
156-
"fluid_pp(1)%pv": pv,
157-
"fluid_pp(1)%gamma_v": gamma_v,
158-
"fluid_pp(1)%M_v": MW_v,
159-
"fluid_pp(1)%k_v": k_v,
160-
"fluid_pp(1)%cp_v": cp_v,
161168
# Bubble gas state
162-
"fluid_pp(2)%gamma": 1.0 / (gamma_g - 1.0),
169+
"fluid_pp(2)%gamma": 1.0 / (gam_g - 1.0),
163170
"fluid_pp(2)%pi_inf": 0.0e00,
164171
"fluid_pp(2)%Re(1)": 1.0 / (mu_g / (rho0 * c0 * x0)),
165-
"fluid_pp(2)%gamma_v": gamma_g,
166-
"fluid_pp(2)%M_v": MW_g,
167-
"fluid_pp(2)%k_v": k_g,
168-
"fluid_pp(2)%cp_v": cp_g,
169172
}
170173
)
171174
)

‎examples/3D_lagrange_bubblescreen/case.py‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -162,16 +162,16 @@
162162
"bub_pp%pv": pv / p0,
163163
"bub_pp%vd": diffVapor / (x0 * c0),
164164
"bub_pp%mu_l": mu_host / (rho0 * x0 * c0),
165-
"bub_pp%gam_v": gam_g,
166-
"bub_pp%gam_g": gam_v,
167-
"bub_pp%M_v": MW_g,
168-
"bub_pp%M_g": MW_v,
169-
"bub_pp%k_v": k_g * (T0 / (x0 * rho0 * c0 * c0 * c0)),
170-
"bub_pp%k_g": k_v * (T0 / (x0 * rho0 * c0 * c0 * c0)),
171-
"bub_pp%cp_v": cp_g * (T0 / (c0 * c0)),
172-
"bub_pp%cp_g": cp_v * (T0 / (c0 * c0)),
173-
"bub_pp%R_v": (R_uni / MW_g) * (T0 / (c0 * c0)),
174-
"bub_pp%R_g": (R_uni / MW_v) * (T0 / (c0 * c0)),
165+
"bub_pp%gam_v": gam_v,
166+
"bub_pp%gam_g": gam_g,
167+
"bub_pp%M_v": MW_v,
168+
"bub_pp%M_g": MW_g,
169+
"bub_pp%k_v": k_v * (T0 / (x0 * rho0 * c0 * c0 * c0)),
170+
"bub_pp%k_g": k_g * (T0 / (x0 * rho0 * c0 * c0 * c0)),
171+
"bub_pp%cp_v": cp_v * (T0 / (c0 * c0)),
172+
"bub_pp%cp_g": cp_g * (T0 / (c0 * c0)),
173+
"bub_pp%R_v": (R_uni / MW_v) * (T0 / (c0 * c0)),
174+
"bub_pp%R_g": (R_uni / MW_g) * (T0 / (c0 * c0)),
175175
# Fluids Physical Parameters
176176
# Host medium
177177
"fluid_pp(1)%gamma": 1.0 / (gamma_host - 1.0),

‎examples/3D_lagrange_shbubcollapse/case.py‎

Lines changed: 10 additions & 10 deletions
Original file line numberDiff line numberDiff line change
@@ -164,16 +164,16 @@
164164
"bub_pp%pv": pv / p0,
165165
"bub_pp%vd": diffVapor / (x0 * c0),
166166
"bub_pp%mu_l": mu_host / (rho0 * x0 * c0),
167-
"bub_pp%gam_v": gam_g,
168-
"bub_pp%gam_g": gam_v,
169-
"bub_pp%M_v": MW_g,
170-
"bub_pp%M_g": MW_v,
171-
"bub_pp%k_v": k_g * (T0 / (x0 * rho0 * c0 * c0 * c0)),
172-
"bub_pp%k_g": k_v * (T0 / (x0 * rho0 * c0 * c0 * c0)),
173-
"bub_pp%cp_v": cp_g * (T0 / (c0 * c0)),
174-
"bub_pp%cp_g": cp_v * (T0 / (c0 * c0)),
175-
"bub_pp%R_v": (R_uni / MW_g) * (T0 / (c0 * c0)),
176-
"bub_pp%R_g": (R_uni / MW_v) * (T0 / (c0 * c0)),
167+
"bub_pp%gam_v": gam_v,
168+
"bub_pp%gam_g": gam_g,
169+
"bub_pp%M_v": MW_v,
170+
"bub_pp%M_g": MW_g,
171+
"bub_pp%k_v": k_v * (T0 / (x0 * rho0 * c0 * c0 * c0)),
172+
"bub_pp%k_g": k_g * (T0 / (x0 * rho0 * c0 * c0 * c0)),
173+
"bub_pp%cp_v": cp_v * (T0 / (c0 * c0)),
174+
"bub_pp%cp_g": cp_g * (T0 / (c0 * c0)),
175+
"bub_pp%R_v": (R_uni / MW_v) * (T0 / (c0 * c0)),
176+
"bub_pp%R_g": (R_uni / MW_g) * (T0 / (c0 * c0)),
177177
# Fluids Physical Parameters
178178
# Host medium
179179
"fluid_pp(1)%gamma": 1.0 / (gamma_host - 1.0),

‎tests/016C1B8B/golden-metadata.txt‎

Lines changed: 54 additions & 85 deletions
Some generated files are not rendered by default. Learn more about customizing how changed files appear on GitHub.

‎tests/016C1B8B/golden.txt‎

Lines changed: 9 additions & 9 deletions
Some generated files are not rendered by default. Learn more about customizing how changed files appear on GitHub.

0 commit comments

Comments
 (0)