PrescribedPressureThermo.h

Go to the documentation of this file.

 
#ifndef RIZER_PRESCRIBED_PRESSURE_THERMO_H
#define RIZER_PRESCRIBED_PRESSURE_THERMO_H
 
#include "thermo/GammaLawThermo.h"
 
#include <functional>
 
namespace rizer
{

    class PrescribedPressureThermo : public ThermoAccessor
    {
    public:
        using PressureFn = std::function<double(double r, double t)>;

        PrescribedPressureThermo(PressureFn pressure, GammaLawThermo backing);
 
        void pressures(double r, double t, double rho, double Tg, double Te,
                       const double *Y, double &p_heavy,
                       double &p_electron) const override;
        double soundSpeed(double r, double t, double rho, double Tg, double Te,
                          const double *Y) const override;
        double heavyEnergy(double Tg, const double *Y) const override;
        double electronCv() const override;
        double densityFromPressure(double pressure, double Tg, double Te,
                                   const double *Y) const override;
        double heavyTemperature(double rho, const double *Y,
                                double heavy_internal_energy_target,
                                double Tg_seed) const override;
 
    private:
        PressureFn m_pressure;
        GammaLawThermo m_backing;
    };
 
} // namespace rizer
 
#endif // RIZER_PRESCRIBED_PRESSURE_THERMO_H