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 <compile_time_options.h>
 
#include "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
{
  namespace EquationOfStateLibrary
  {
    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 EquationOfStateLibrary
} /* namespace ryujin */