8#include <compile_time_options.h>
14 namespace EquationOfStateLibrary
33 this->add_parameter(
"gamma", gamma_,
"The ratio of specific heats");
42 "gas constant R", R_,
"The specific gas constant R");
44 cv_ = R_ / (gamma_ - 1.);
48 "covolume b", b_,
"The maximum compressibility constant");
51 this->add_parameter(
"reference specific internal energy",
53 "The reference specific internal energy");
57 "reference pressure", pinf_,
"The reference pressure p infinity");
60 ParameterAcceptor::parse_parameters_call_back.connect([
this] {
73 double pressure(
double rho,
double e)
const final
75 return (gamma_ - 1.) * rho * (e - q_) / (1. - b_ * rho) -
88 const auto numerator = (p + gamma_ * pinf_) * (1. - b_ * rho);
89 const auto denominator = rho * (gamma_ - 1.);
90 return q_ + numerator / denominator;
101 return (e - q_ - pinf_ * (1. / rho - b_)) / cv_;
113 const auto covolume = 1. - b_ * rho;
115 (rho * (e - q_) - pinf_ * covolume) / (covolume * covolume * rho);
116 radicand *= gamma_ * (gamma_ - 1.);
117 return std::sqrt(radicand);
virtual double specific_internal_energy(double rho, double p) const =0
virtual double pressure(double rho, double e) const =0
virtual double speed_of_sound(double rho, double e) const =0
virtual double temperature(double rho, double e) const =0
double interpolation_pinfty_
double speed_of_sound(double rho, double e) const final
NobleAbelStiffenedGas(const std::string &subsection)
double temperature(double rho, double e) const final
double pressure(double rho, double e) const final
double specific_internal_energy(double rho, double p) const final