12 namespace EquationOfStateLibrary
26 this->add_parameter(
"gamma", gamma_,
"The ratio of specific heats");
29 this->add_parameter(
"vdw a", a_,
"The vdw a constant");
33 "covolume b", b_,
"The maximum compressibility constant");
42 "gas constant R", R_,
"The specific gas constant R");
44 cv_ = R_ / (gamma_ - 1.);
47 ParameterAcceptor::parse_parameters_call_back.connect(
57 double pressure(
double rho,
double e)
const final
59 const auto intermolecular = a_ * rho * rho;
60 const auto numerator = rho * e + intermolecular;
61 const auto covolume = 1. - b_ * rho;
62 return (gamma_ - 1.) * numerator / covolume - intermolecular;
74 const auto intermolecular = a_ * rho * rho;
75 const auto covolume = 1. - b_ * rho;
76 const auto numerator = (p + intermolecular) * covolume;
77 const auto denominator = rho * (gamma_ - 1.);
78 return numerator / denominator - a_ * rho;
89 return (e + a_ * rho) / cv_;
101 const auto covolume = 1. - b_ * rho;
102 const auto numerator = gamma_ * (gamma_ - 1.) * (e + a_ * rho);
103 return std::sqrt(numerator / (covolume * covolume) - 2. * a_ * rho);
VanDerWaals(const std::string &subsection)
double temperature(double rho, double e) const final
double pressure(double rho, double e) const final
double speed_of_sound(double rho, double e) const final
double specific_internal_energy(double rho, double p) const final