GammaLawThermo.h

Go to the documentation of this file.

 
#ifndef RIZER_GAMMA_LAW_THERMO_H
#define RIZER_GAMMA_LAW_THERMO_H
 
#include "thermo/ThermoAccessor.h"
 
#include <vector>
 
namespace rizer
{

    class GammaLawThermo : public ThermoAccessor
    {
    public:
        GammaLawThermo(double gamma, std::vector<double> molecular_weights,
                       int electron_index);
 
        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:
        double heavyMolesPerMass(const double *Y) const;
        double electronMolesPerMass(const double *Y) const;
 
        double m_gamma; 
        std::vector<double> m_W;
        int m_electron_index;
    };
 
} // namespace rizer
 
#endif // RIZER_GAMMA_LAW_THERMO_H