equation_of_state_polytropic_gas.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2023 - 2025 by the ryujin authors
//
 
#pragma once
 
#include <compile_time_options.h>
 
#include "equation_of_state.h"
 
namespace ryujin
{
  namespace EquationOfStateLibrary
  {
    class PolytropicGas : public EquationOfState
    {
    public:
      using EquationOfState::pressure;
      using EquationOfState::specific_internal_energy;
      using EquationOfState::speed_of_sound;
      using EquationOfState::temperature;
 
      PolytropicGas(const std::string &subsection)
          : EquationOfState("polytropic gas", subsection)
      {
        gamma_ = 7. / 5.;
        this->add_parameter("gamma", gamma_, "The ratio of specific heats");
 
        /*
         * R is the specific gas constant with units [J / (Kg K)]. More details
         * can be found at:
         * https://en.wikipedia.org/wiki/Gas_constant#Specific_gas_constant
         */
        R_ = 287.052874;
        this->add_parameter(
            "gas constant R", R_, "The specific gas constant R");
 
        cv_ = R_ / (gamma_ - 1.);
 
        s0_ = 0.;
        this->add_parameter("reference specific entropy",
                            s0_,
                            "The reference specific entropy");
 
        /* Update the EOS interpolation parameters on parameter read in: */
        ParameterAcceptor::parse_parameters_call_back.connect(
            [this] { cv_ = R_ / (gamma_ - 1.); });
      }
 
      double pressure(double rho, double e) const final
      {
        return (gamma_ - 1.) * rho * e;
      }
 
      double specific_internal_energy(double rho, double p) const final
      {
        return p / (rho * (gamma_ - 1.));
      }
 
      double temperature(double /*rho*/, double e) const final
      {
        return e / cv_;
      }
 
      double cold_curve_bound(double /*rho*/) const final
      {
        return 0.;
      }
 
      double specific_entropy(double rho, double e) const final
      {
        const auto gm1 = gamma_ - 1.;
        const auto p = gm1 * rho * e;
        return cv_ * (std::log(p) - gamma_ * std::log(gm1 * cv_ * rho)) + s0_;
      }
 
      double speed_of_sound(double /*rho*/, double e) const final
      {
        return std::sqrt(gamma_ * (gamma_ - 1.) * e);
      }
 
    private:
      double gamma_;
      double R_;
      double cv_;
      double s0_;
    };
  } // namespace EquationOfStateLibrary
} // namespace ryujin