Merge pull request #91 from LLNL/bugfix/StiffenedGasEOS2

error in stiffened gas bulk mod made it way too large
LLNL · Oct 4, 2021 · 25bb5d7 · 25bb5d7
2 parents e47116d + d92cd1d
commit 25bb5d7
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/Material/StiffenedGas.cc b/src/Material/StiffenedGas.cc
@@ -230,19 +230,19 @@ StiffenedGas<Dimension>::
 bulkModulus(const Scalar massDensity,
             const Scalar specificThermalEnergy) const {
   CHECK(valid());
-  return massDensity*max(0.0, mGamma*mGamma1*(specificThermalEnergy+mP0/massDensity));
+  return massDensity*max(0.0, mGamma*mGamma1*(specificThermalEnergy-mP0/massDensity));
 }
 
 //------------------------------------------------------------------------------
-// Calculate an entropy. This should be double checked
+// Calculate an entropy.
 //------------------------------------------------------------------------------
 template<typename Dimension>
 typename Dimension::Scalar
 StiffenedGas<Dimension>::
 entropy(const Scalar massDensity,
         const Scalar specificThermalEnergy) const {
   CHECK(valid());
-  return (mGamma*mP0+this->pressure(massDensity, specificThermalEnergy))*safeInvVar(pow(massDensity, mGamma));
+  return (mP0+this->pressure(massDensity, specificThermalEnergy))*safeInvVar(pow(massDensity, mGamma));
 }
 
 //------------------------------------------------------------------------------