PrescribedPressureThermo.cpp

Go to the documentation of this file.

#include "thermo/PrescribedPressureThermo.h"
 
#include <cmath>
#include <stdexcept>
#include <string>
#include <utility>
 
namespace rizer
{
 
    // ------------------------------------------------------------------------------ //
    //   PrescribedPressureThermo - constructor                                       //
    // ------------------------------------------------------------------------------ //
 
    PrescribedPressureThermo::PrescribedPressureThermo(PressureFn pressure,
                                                       GammaLawThermo backing)
        : m_pressure(std::move(pressure)), m_backing(std::move(backing))
    {
        if (!m_pressure)
        {
            throw std::invalid_argument(
                "PrescribedPressureThermo: pressure function is empty.");
        }
    }
 
    // ------------------------------------------------------------------------------ //
    //   PrescribedPressureThermo - public methods                                    //
    // ------------------------------------------------------------------------------ //
 
    void PrescribedPressureThermo::pressures(
        double r,
        double t,
        double /*rho*/,
        double /*Tg*/,
        double /*Te*/,
        const double * /*Y*/,
        double &p_heavy,
        double &p_electron) const
    {
        p_heavy = m_pressure(r, t);
        // The callback is arbitrary caller code and its result is the pressure
        // every flux at this face is built from, so a non-finite one is rejected
        // where it enters rather than where it finally shows.
        if (!std::isfinite(p_heavy))
        {
            throw std::runtime_error(
                "PrescribedPressureThermo::pressures: the prescribed pressure at "
                "r = " +
                std::to_string(r) + " m, t = " + std::to_string(t) + " s is " +
                std::to_string(p_heavy) + " Pa, which is not finite.");
        }
        p_electron = 0.0;
    }
 
    double PrescribedPressureThermo::soundSpeed(
        double r,
        double t,
        double rho,
        double Tg,
        double Te,
        const double *Y) const
    {
        return m_backing.soundSpeed(r, t, rho, Tg, Te, Y);
    }
 
    double PrescribedPressureThermo::heavyEnergy(double Tg, const double *Y) const
    {
        return m_backing.heavyEnergy(Tg, Y);
    }
 
    double PrescribedPressureThermo::electronCv() const
    {
        return m_backing.electronCv();
    }
 
    double PrescribedPressureThermo::densityFromPressure(
        double pressure,
        double Tg,
        double Te,
        const double *Y) const
    {
        return m_backing.densityFromPressure(pressure, Tg, Te, Y);
    }
 
    double PrescribedPressureThermo::heavyTemperature(
        double rho,
        const double *Y,
        double heavy_internal_energy_target,
        double Tg_seed) const
    {
        return m_backing.heavyTemperature(rho, Y, heavy_internal_energy_target, Tg_seed);
    }
 
} // namespace rizer