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 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.);
      }
 
      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 speed_of_sound(double /*rho*/, double e) const final
      {
        return std::sqrt(gamma_ * (gamma_ - 1.) * e);
      }
 
    private:
      double gamma_;
      double R_;
      double cv_;
    };
  } // namespace EquationOfStateLibrary
} // namespace ryujin