PlasmaReactor0D.h

Go to the documentation of this file.

 
#ifndef RIZER_PLASMA_REACTOR_0D_H
#define RIZER_PLASMA_REACTOR_0D_H
 
#ifndef NOMINMAX
#define NOMINMAX
#endif
#include "cantera/base/Solution.h"
#include "PlasmaCell.h"
#include "transport/CollisionModel.h"
 
#include "cantera/thermo/PlasmaPhase.h"
 
#include <memory>
#include <string>
#include <vector>
 
namespace rizer
{
 
    struct PlasmaReactor0DConfig
    {
        std::string mech, phase; 
        bool reacting = true;
        std::vector<CollisionModel::SpeciesSpec> specs;
        double Te_min = 200.0;
        double Te_max = 1.0e5;
        std::size_t Te_n = 400;
        double gap = 0.0;
        std::vector<double> species_h0k;
    };
 
    class PlasmaReactor0D
    {
    public:
        explicit PlasmaReactor0D(const PlasmaReactor0DConfig &cfg);

        std::size_t nSpecies() const { return m_rhs->nSpecies(); }

        std::vector<std::string> speciesNames() const;

        double electricalConductivity(double Tg, double Te, const double *Y,
                                      double rho) const
        {
            return m_rhs->electricalConductivity(Tg, Te, Y, rho);
        }

        double electronDensity(double Tg, double Te, const double *Y, double rho) const
        {
            return m_rhs->electronDensity(Tg, Te, Y, rho);
        }

        void rhs(double Tg, double Te, double V, const double *Y, double E, double mass,
                 double p_ext, double polytropic_index, double *dydt) const;

        const PlasmaCell::Diagnostics &diagnostics() const { return m_diag; }
 
    private:
        Cantera::PlasmaPhase *m_plasma = nullptr;

        std::shared_ptr<Cantera::Solution> m_sol;
        std::unique_ptr<PlasmaCell> m_rhs;
        std::size_t m_ie;
        double m_gap;
        mutable PlasmaCell::Diagnostics m_diag;
    };
 
} // namespace rizer
 
#endif