Plasma0DReactor.h

Go to the documentation of this file.

 
#ifndef RIZER_PLASMA_0D_REACTOR_H
#define RIZER_PLASMA_0D_REACTOR_H
 
#ifndef NOMINMAX
#define NOMINMAX
#endif
#include "cantera/base/Solution.h"
#include "ReactorRHS.h"
#include "transport/CollisionModel.h"
 
#include <memory>
#include <string>
#include <vector>
 
namespace rizer {
 
class Plasma0DReactor {
public:
    struct Config {
        std::string mech, phase;      
        bool reacting = true;
        std::vector<CollisionModel::SpeciesSpec> specs;
        double Te_min = 300.0;
        double Te_max = 1.0e5;
        std::size_t Te_n = 400;
        bool spitzer = true;
        double gap = 0.0;   
        std::vector<double> species_h0k;
    };

    explicit Plasma0DReactor(const Config& cfg);
 
    std::size_t nSpecies() const { return m_rhs->nSpecies(); }
    std::vector<std::string> speciesNames() const;

    double conductivity(double Tg, double Te, const double* Y, double rho) const {
        return m_rhs->conductivity(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 ReactorRHS::Diagnostics& diagnostics() const { return m_diag; }
 
private:
    std::shared_ptr<Cantera::Solution> m_sol;
    std::unique_ptr<ReactorRHS> m_rhs;
    std::size_t m_ie;
    double m_gap;
    mutable ReactorRHS::Diagnostics m_diag;
};
 
} // namespace rizer
 
#endif