Skip to content

Commit

Permalink
remove last use of PyroWrapper.wts (#1054)
Browse files Browse the repository at this point in the history
  • Loading branch information
majosm authored Aug 16, 2024
1 parent 159f8e1 commit 0d680df
Showing 1 changed file with 1 addition and 3 deletions.
4 changes: 1 addition & 3 deletions mirgecom/thermochemistry.py
Original file line number Diff line number Diff line change
Expand Up @@ -85,7 +85,6 @@ class PyroWrapper(pyro_class):
# This only affects chemistry-related evaluations and does not interfere
# with the actual fluid state.
def get_concentrations(self, rho, mass_fractions):
# concs = self.inv_molecular_weights * rho * mass_fractions
concs = self.inv_molecular_weights * rho * mass_fractions
# ensure non-negative concentrations
zero = self._pyro_zeros_like(concs[0])
Expand Down Expand Up @@ -177,8 +176,7 @@ def get_heat_release_rate(self, state):

heat_rls = state.cv.mass*0.0
for i in range(self.num_species):
# heat_rls = heat_rls - h_a[i]*w_dot[i]/(self.molecular_weights[i])
heat_rls = heat_rls - h_a[i]*w_dot[i]/(self.wts[i])
heat_rls = heat_rls - h_a[i]*w_dot[i]/(self.molecular_weights[i])

return heat_rls*self.gas_constant*state.temperature

Expand Down

0 comments on commit 0d680df

Please sign in to comment.