EOSModel.H

Go to the documentation of this file.

#ifndef EOSMODEL_H
#define EOSMODEL_H
 
#include "AMReX_REAL.H"
#include "amr-wind/utilities/constants.H"
 
namespace amr_wind::eos {
 
struct GammaLaw
{
    using eos_type = GammaLaw;
    static std::string identifier() { return "GammaLaw"; }
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real p_rth(
        const amrex::Real rho,
        const amrex::Real theta,
        const amrex::Real qv = 0.0) const
    {
        const auto rhotheta = rho * theta;
        return m_p0 * std::pow(
                          m_air_gas_constant * rhotheta *
                              (1.0 + (m_water_vapor_gas_constant /
                                      m_air_gas_constant) *
                                         qv) *
                              m_ip0,
                          m_gamma);
    }
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real p_rth( {…}
 
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real dp_constanttheta(
        const amrex::Real rho,
        const amrex::Real theta,
        const amrex::Real qv = 0.0) const
    {
        return m_gamma * m_p0 *
               std::pow(
                   (m_air_gas_constant * theta *
                    (1.0 +
                     m_water_vapor_gas_constant / m_air_gas_constant * qv) *
                    m_ip0),
                   m_gamma) *
               std::pow(rho, m_gamma - 1.0);
    }
    AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real dp_constanttheta( {…}
 
    template <class... Args>
    AMREX_GPU_HOST_DEVICE explicit GammaLaw(const amrex::Real p0 = 1e-5)
        : m_p0{p0}, m_ip0{1.0 / p0}
    {}
    AMREX_GPU_HOST_DEVICE explicit GammaLaw(const amrex::Real p0 = 1e-5) {…}
 
    const amrex::Real m_gamma{constants::HEAT_CAPACITY_RATIO};
    const amrex::Real m_p0{1e5};
    const amrex::Real m_ip0{1e-5};
    const amrex::Real m_air_gas_constant{
        constants::UNIVERSAL_GAS_CONSTANT / constants::MOLAR_MASS_AIR};
    const amrex::Real m_air_gas_constant {…};
    const amrex::Real m_water_vapor_gas_constant{
        constants::UNIVERSAL_GAS_CONSTANT / constants::MOLAR_MASS_WATER_VAPOR};
    const amrex::Real m_water_vapor_gas_constant {…};
};
struct GammaLaw {…};
} // namespace amr_wind::eos
namespace amr_wind::eos {…}
 
#endif /* EOSMODEL_H */