CanteraThermo.h

Go to the documentation of this file.

 
#ifndef RIZER_CANTERA_THERMO_H
#define RIZER_CANTERA_THERMO_H
 
#include "thermo/ThermoAccessor.h"
 
#include "cantera/thermo/PlasmaPhase.h"
 
#include <memory>
#include <string>
#include <vector>
 
namespace rizer
{

    class CanteraThermo : public ThermoAccessor
    {
    public:
        CanteraThermo(const std::string &mech, const std::string &phase);

        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;

        StateProperties stateProperties(double r, double t, double rho, double Tg,
                                        double Te,
                                        const double *Y) const override;

        std::size_t nSpecies() const { return m_nsp; }

        std::size_t electronIndex() const { return m_ie; }
 
    private:
        void setState(double rho, double Tg, double Te, const double *Y) const;

        void heavyEnergyAndCv(double Tg, const double *Y, double &internal_energy,
                              double &cv) const;

        void heavyCpAndCv(double &cp, double &cv) const;

        double soundSpeedAtSetState(double p_heavy, double p_electron,
                                    double Tg) const;

        double heavyEnergyAtSetState() const;
 
        std::shared_ptr<Cantera::ThermoPhase> m_phase;
        Cantera::ThermoPhase *m_thermo = nullptr;
        Cantera::PlasmaPhase *m_plasma = nullptr;
        std::size_t m_nsp = 0;
        std::size_t m_ie = 0;               
        double m_W_electron = 0.0;          
        mutable std::vector<double> m_Yc;   
        mutable std::vector<double> m_work; 
    };
 
} // namespace rizer
 
#endif // RIZER_CANTERA_THERMO_H