equation_of_state.h

Go to the documentation of this file.

//
// SPDX-License-Identifier: Apache-2.0 WITH LLVM-exception
// Copyright (C) 2023 - 2026 by the ryujin authors
//
 
#pragma once
 
#include <ryujin/base/compile_time_options.h>
#include <ryujin/base/convenience_macros.h>
 
#include <deal.II/base/exceptions.h>
#include <deal.II/base/parameter_acceptor.h>
#include <deal.II/base/tensor.h>
 
#include <string>
 
namespace ryujin
{
  class EquationOfState : public dealii::ParameterAcceptor
  {
  public:
    EquationOfState(const std::string &name, const std::string &subsection)
        : ParameterAcceptor(subsection + "/" + name)
        , name_(name)
    {
      /*
       * Derived EOS can override the covolume constant b.
       */
      covolume_constant_ = 0.;
 
      /*
       * Derived EOS can override the reference pressure that is used in
       * the interpolatory NASG approach.
       */
      interpolation_pinfty_ = 0.;
 
      /*
       * Derived EOS can override the reference specific internal energy
       * q that is used in the interpolatory NASG approach.
       */
      interpolation_q_ = 0.;
    }
 
    virtual double pressure(double rho, double e) const = 0;
 
    virtual double specific_internal_energy(double rho, double p) const = 0;
 
    virtual double specific_entropy(double /*rho*/, double /*e*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
      __builtin_trap();
      return 0;
    }
 
    virtual double cold_curve_bound(double /*rho*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
      __builtin_trap();
      return 0;
    }
 
    virtual double temperature(double /*rho*/, double /*e*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
      __builtin_trap();
      return 0;
    }
 
    virtual double speed_of_sound(double /*rho*/, double /*e*/) const
    {
      AssertThrow(false, dealii::ExcNotImplemented());
      __builtin_trap();
      return 0;
    }
 
    ACCESSOR_READ_ONLY(covolume_constant)
 
    
    ACCESSOR_READ_ONLY(interpolation_pinfty)
 
    
    ACCESSOR_READ_ONLY(interpolation_q)
 
    
    ACCESSOR_READ_ONLY(name)
 
  protected:
    double covolume_constant_;
 
    double interpolation_pinfty_;
    double interpolation_q_;
 
  private:
    const std::string name_;
  };
} // namespace ryujin