ThermalPlasmaColumn1D.h

Go to the documentation of this file.

 
#ifndef RIZER_THERMAL_PLASMA_COLUMN_1D_H
#define RIZER_THERMAL_PLASMA_COLUMN_1D_H
 
#ifndef NOMINMAX
#define NOMINMAX
#endif
#include "cantera/oneD/Domain1D.h"
#include "PropertyTable.h"
 
#include <string>
 
namespace rizer {
 
class ThermalPlasmaColumn1D : public Cantera::Domain1D {
public:
    ThermalPlasmaColumn1D(double R, std::size_t npoints,
                          double electric_field, double T_wall,
                          PropertyTable sigma, PropertyTable kappa,
                          PropertyTable p_rad,
                          double T_center_guess,
                          double rho_cp,
                          PropertyTable init_profile = PropertyTable());
 
    std::string domainType() const override { return "thermal-plasma-column"; }
 
    // Residual of the cylindrical heat equation. span-based signature per
    // Cantera main (>3.2). xg/rg are GLOBAL spans; the local part is
    // xg.subspan(loc(), size()).

    void eval(std::size_t jg, Cantera::span<const double> xg,
              Cantera::span<double> rg, Cantera::span<int> maskg,
              double rdt) override;
 
    std::string componentName(std::size_t n) const override { return "T"; }
    std::size_t componentIndex(const std::string& name,
                               bool checkAlias = true) const override;
    bool hasComponent(const std::string& name,
                      bool checkAlias = true) const override;

    void getValues(const std::string& component,
                   Cantera::span<double> values) const override;

    double initialValue(std::size_t n, std::size_t j) override;

    void resetBadValues(Cantera::span<double> xg) override;

    void setElectricField(double E) { m_E = E; }
    double electricField() const { return m_E; }
    double radius() const { return m_R; }

    double sigmaAt(double T) const { return m_sigma.eval(T); }
    double kappaAt(double T) const { return m_kappa.eval(T); }
    double pRadAt(double T) const { return m_prad.eval(T); }
 
private:
    double m_R;        
    double m_E;        
    double m_Twall;    
    double m_Tcenter;  
    double m_rho_cp;   
    PropertyTable m_sigma;  
    PropertyTable m_kappa;  
    PropertyTable m_prad;   
    PropertyTable m_init;   
};
 
} // namespace rizer
 
#endif