ThermoAccessor.h

Go to the documentation of this file.

 
#ifndef RIZER_THERMO_ACCESSOR_H
#define RIZER_THERMO_ACCESSOR_H
 
namespace rizer
{

    constexpr double ElectronGamma = 5.0 / 3.0;

    struct StateProperties
    {
        double p_heavy = 0.0;      
        double p_electron = 0.0;   
        double sound_speed = 0.0;  
        double heavy_energy = 0.0; 
    };

    class ThermoAccessor
    {
    public:
        virtual ~ThermoAccessor() = default;

        virtual void pressures(double r, double t, double rho, double Tg, double Te,
                               const double *Y, double &p_heavy,
                               double &p_electron) const = 0;

        virtual double soundSpeed(double r, double t, double rho, double Tg,
                                  double Te, const double *Y) const = 0;

        virtual double heavyEnergy(double Tg, const double *Y) const = 0;

        virtual double electronCv() const = 0;

        virtual double densityFromPressure(double pressure, double Tg, double Te,
                                           const double *Y) const = 0;

        virtual double heavyTemperature(double rho, const double *Y,
                                        double heavy_internal_energy_target,
                                        double Tg_seed) const = 0;

        virtual StateProperties stateProperties(double r, double t, double rho,
                                                double Tg, double Te,
                                                const double *Y) const
        {
            StateProperties properties;
            pressures(r, t, rho, Tg, Te, Y, properties.p_heavy,
                      properties.p_electron);
            properties.sound_speed = soundSpeed(r, t, rho, Tg, Te, Y);
            properties.heavy_energy = heavyEnergy(Tg, Y);
            return properties;
        }
    };
 
} // namespace rizer
 
#endif // RIZER_THERMO_ACCESSOR_H