Skip to content

Commit

Permalink
Fix pyro mech naming issue
Browse files Browse the repository at this point in the history
  • Loading branch information
MTCam committed May 5, 2022
1 parent d9c3867 commit c62ae04
Showing 1 changed file with 16 additions and 15 deletions.
31 changes: 16 additions & 15 deletions combozzle.py
Original file line number Diff line number Diff line change
Expand Up @@ -742,10 +742,10 @@ def vol_max(x):
else:
from mirgecom.thermochemistry import get_pyrometheus_wrapper_class
from mirgecom.mechanisms.uiuc import Thermochemistry
pyro_mech = get_pyrometheus_wrapper_class(Thermochemistry)(actx.np)
nspecies = pyro_mech.num_species
species_names = pyro_mech.species_names
eos = PyrometheusMixture(pyro_mech, temperature_guess=temperature_seed)
pyro_mechanism = get_pyrometheus_wrapper_class(Thermochemistry)(actx.np)
nspecies = pyro_mechanism.num_species
species_names = pyro_mechanism.species_names
eos = PyrometheusMixture(pyro_mechanism, temperature_guess=temperature_seed)
init_y = [0.06372925, 0.21806609, 0., 0., 0., 0., 0.71820466]

# {{{ Initialize simple transport model
Expand Down Expand Up @@ -1032,7 +1032,7 @@ def my_post_step(step, t, dt, state):

return make_obj_array([cv, fluid_state.temperature]), dt

from mirgecom.inviscid import inviscid_flux_rusanov
from mirgecom.inviscid import inviscid_facial_flux_rusanov

def dummy_pre_step(step, t, dt, state):
if logmgr:
Expand All @@ -1059,15 +1059,17 @@ def cfd_rhs(t, state):
temperature_seed=tseed)
if inviscid_only:
fluid_rhs = \
euler_operator(discr, state=fluid_state, time=t,
boundaries=boundaries, gas_model=gas_model,
inviscid_numerical_flux_func=inviscid_flux_rusanov,
quadrature_tag=quadrature_tag)
euler_operator(
discr, state=fluid_state, time=t,
boundaries=boundaries, gas_model=gas_model,
inviscid_numerical_flux_func=inviscid_facial_flux_rusanov,
quadrature_tag=quadrature_tag)
else:
fluid_rhs = \
ns_operator(discr, state=fluid_state, time=t, boundaries=boundaries,
gas_model=gas_model, quadrature_tag=quadrature_tag,
inviscid_numerical_flux_func=inviscid_flux_rusanov)
ns_operator(
discr, state=fluid_state, time=t, boundaries=boundaries,
gas_model=gas_model, quadrature_tag=quadrature_tag,
inviscid_numerical_flux_func=inviscid_facial_flux_rusanov)

if inert_only:
chem_rhs = 0*fluid_rhs
Expand All @@ -1077,9 +1079,8 @@ def cfd_rhs(t, state):
if av_on:
alpha_f = compute_av_alpha_field(fluid_state)
av_rhs = av_laplacian_operator(discr, fluid_state=fluid_state,
boundaries=boundaries,
boundary_kwargs={"time": t,
"gas_model": gas_model},
boundaries=boundaries, time=t,
gas_model=gas_model,
alpha=alpha_f, s0=s0_sc, kappa=kappa_sc)
else:
av_rhs = 0*fluid_rhs
Expand Down

0 comments on commit c62ae04

Please sign in to comment.