Skip to content

Commit

Permalink
small fixes for CI
Browse files Browse the repository at this point in the history
  • Loading branch information
@mrodrig6
mrodrig6 committed Jan 31, 2025
1 parent f2c9c61 commit 2984369
Show file tree
Hide file tree
Showing 2 changed files with 93 additions and 267 deletions.
182 changes: 93 additions & 89 deletions examples/1D_copper/case.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2,98 +2,102 @@
import math
import json

#Numerical setup
c_l = 3910
Nx = 256
cfl = 0.1
leng = 1.
dx = leng/(Nx+1)
mydt = cfl*dx/c_l
Tend = 100.E-6
Nt = int(Tend/mydt)
#mydt = Tend/(1.*Nt)
vel1 = 0.0
vel2 = 0.0
# Numerical setup
c_l = 3910
Nx = 256
cfl = 0.1
leng = 1.0
dx = leng / (Nx + 1)
mydt = cfl * dx / c_l
Tend = 100.0e-6
Nt = int(Tend / mydt)
# mydt = Tend/(1.*Nt)
vel1 = 0.0
vel2 = 0.0
theta_0 = 298.0

#Initial condition
theta_0 = 298 #K
P_0 = 1.000E5 #Pa
compression_ratio = 1.1 #rho/rho_0 in the shocked region
rho_0 = 8924.0 #kg/m^3
vel0 = 0.0
s = 1.51
# Initial condition
theta_0 = 298 # K
P_0 = 1.000e5 # Pa
compression_ratio = 1.1 # rho/rho_0 in the shocked region
rho_0 = 8924.0 # kg/m^3
vel0 = 0.0
s = 1.51
tilde_rho = compression_ratio

xi = 1.0-1.0/tilde_rho
ps = P_0 + c_l*c_l*rho_0*xi/((1.0-s*xi)*(1.0-s*xi))
vel = vel0 + math.sqrt((ps-P_0)*xi/rho_0)
xi = 1.0 - 1.0 / tilde_rho
ps = P_0 + c_l * c_l * rho_0 * xi / ((1.0 - s * xi) * (1.0 - s * xi))
vel = vel0 + math.sqrt((ps - P_0) * xi / rho_0)

#phi = math.exp(gamma_suc*(1-1/tilde_rho))
# phi = math.exp(gamma_suc*(1-1/tilde_rho))
# Configuring case dictionary
print(json.dumps({
# Logistics
"run_time_info": "T",
# Computational Domain Parameters
"x_domain%beg": 0.E+00,
"x_domain%end": 1.E+00,
"m": Nx,
"n": 0,
"p": 0,
"dt": mydt,
"t_step_start": 0,
"t_step_stop": int(Nt),
"t_step_save": int(math.ceil(Nt/100.)),
# Simulation Algorithm Parameters
"num_patches": 2,
"model_eqns": 5,
"alt_soundspeed": "F",
"num_fluids": 1,
"MGEoS_model": 1,
"mpp_lim": "F",
"mixture_err": "F",
"time_stepper": 3,
"weno_order": 5,
"weno_eps": 1.E-16,
"weno_Re_flux": "F",
"weno_avg": "T",
"mapped_weno": "T",
"null_weights": "T",
"mp_weno": "T",
"riemann_solver": 2,
"wave_speeds": 1,
"avg_state": 2,
"bc_x%beg": -3,
"bc_x%end": -3,
# Formatted Database Files Structure Parameters
"format": 1,
"precision": 2,
"prim_vars_wrt":"T",
"parallel_io":"T",
# Patch 1: Equilibrium state
"patch_icpp(1)%geometry": 1,
"patch_icpp(1)%x_centroid": 0.5,
"patch_icpp(1)%length_x": leng,
"patch_icpp(1)%vel(1)": vel0,
"patch_icpp(1)%pres": P_0,
"patch_icpp(1)%alpha_rho(1)": rho_0,
"patch_icpp(1)%alpha(1)": 1.0,
# Patch 2: Shocked state
"patch_icpp(2)%geometry": 1,
"patch_icpp(2)%x_centroid": 0.0625,
"patch_icpp(2)%length_x": 0.125,
"patch_icpp(2)%alter_patch(1)": "T",
"patch_icpp(2)%vel(1)": vel,
"patch_icpp(2)%pres": ps,
"patch_icpp(2)%alpha_rho(1)": tilde_rho*rho_0,
"patch_icpp(2)%alpha(1)": 1.0,
# Fluids Physical Parameters
"fluid_pp(1)%gamma": 1.96E0, # Gruneisen constant
"fluid_pp(1)%pi_inf": 1.0E5, # p0
"fluid_pp(1)%qv": 0.0, # e0
"fluid_pp(1)%qvp": 1.0, # Gruneisen exponent
"fluid_pp(1)%mg_a": 3910.0, # c0
"fluid_pp(1)%mg_b": 1.51, # s
"fluid_pp(1)%rho0": 8924.0, # refernce density
"fluid_pp(1)%cv": 385.0, # specific heat capacity
}))
print(
json.dumps(
{
# Logistics
"run_time_info": "T",
# Computational Domain Parameters
"x_domain%beg": 0.0e00,
"x_domain%end": 1.0e00,
"m": Nx,
"n": 0,
"p": 0,
"dt": mydt,
"t_step_start": 0,
"t_step_stop": int(Nt),
"t_step_save": int(math.ceil(Nt / 100.0)),
# Simulation Algorithm Parameters
"num_patches": 2,
"model_eqns": 5,
"alt_soundspeed": "F",
"num_fluids": 1,
"MGEoS_model": 1,
"mpp_lim": "F",
"mixture_err": "F",
"time_stepper": 3,
"weno_order": 5,
"weno_eps": 1.0e-16,
"weno_Re_flux": "F",
"weno_avg": "T",
"mapped_weno": "T",
"null_weights": "T",
"mp_weno": "T",
"riemann_solver": 2,
"wave_speeds": 1,
"avg_state": 2,
"bc_x%beg": -3,
"bc_x%end": -3,
# Formatted Database Files Structure Parameters
"format": 1,
"precision": 2,
"prim_vars_wrt": "T",
"parallel_io": "T",
# Patch 1: Equilibrium state
"patch_icpp(1)%geometry": 1,
"patch_icpp(1)%x_centroid": 0.5,
"patch_icpp(1)%length_x": leng,
"patch_icpp(1)%vel(1)": vel0,
"patch_icpp(1)%pres": P_0,
"patch_icpp(1)%alpha_rho(1)": rho_0,
"patch_icpp(1)%alpha(1)": 1.0,
# Patch 2: Shocked state
"patch_icpp(2)%geometry": 1,
"patch_icpp(2)%x_centroid": 0.0625,
"patch_icpp(2)%length_x": 0.125,
"patch_icpp(2)%alter_patch(1)": "T",
"patch_icpp(2)%vel(1)": vel,
"patch_icpp(2)%pres": ps,
"patch_icpp(2)%alpha_rho(1)": tilde_rho * rho_0,
"patch_icpp(2)%alpha(1)": 1.0,
# Fluids Physical Parameters
"fluid_pp(1)%gamma": 1.96e0, # Gruneisen constant
"fluid_pp(1)%pi_inf": 1.0e5, # p0
"fluid_pp(1)%qv": 0.0, # e0
"fluid_pp(1)%qvp": 1.0, # Gruneisen exponent
"fluid_pp(1)%mg_a": 3910.0, # c0
"fluid_pp(1)%mg_b": 1.51, # s
"fluid_pp(1)%rho0": 8924.0, # reference density
"fluid_pp(1)%cv": 385.0, # specific heat capacity
}
)
)
178 changes: 0 additions & 178 deletions examples/2D_lungwave/case.py
View file

This file was deleted.

0 comments on commit 2984369

Please sign in to comment.