PlasmaColumnSolver.h

Go to the documentation of this file.

// PlasmaColumnSolver: assembles [Empty1D, ThermalPlasmaColumn1D, Empty1D] into
// a Cantera Sim1D, drives the steady solve with adaptive grid refinement, and
// returns the converged radial profile as plain numeric vectors.
//
// The three-domain layout  [Empty1D | ThermalPlasmaColumn1D | Empty1D]  is the
// standard Cantera pattern for a BVP that does not need inlet/outlet domains.
// The Empty1D terminators satisfy the Sim1D API without contributing any
// residual rows.
//
// No Cantera C++ type leaves this translation unit; the pybind11 bindings in
// bindings.cpp convert the output vectors to numpy arrays.  This keeps all ABI
// coupling with libcantera inside one .so, away from the stock `cantera`
// Python package.
 
#ifndef RIZER_PLASMA_COLUMN_SOLVER_H
#define RIZER_PLASMA_COLUMN_SOLVER_H
 
#include "PropertyTable.h"
 
#include <vector>
#include <cstddef>
 
namespace rizer {

struct ColumnResult {
    std::vector<double> r;        
    std::vector<double> T;        
    std::vector<double> sigma;    
    std::vector<double> kappa;    
    double current = 0.0;         
    double electric_field = 0.0;  
    std::size_t n_points = 0;     
};

struct ColumnOptions {
    std::size_t n_points = 41;     
    double T_center_guess = 10000; 
    double rho_cp = 1.0e3;          
    int loglevel = 0;              
    bool refine_grid = true;       
    double refine_ratio = 4.0;     
    double refine_slope = 0.05;    
    double refine_curve = 0.10;    
    double refine_prune = 0.0;     

    std::vector<double> init_r;    
    std::vector<double> init_T;    
};

ColumnResult solveColumn(double R, double electric_field, double T_wall,
                          const PropertyTable& sigma, const PropertyTable& kappa,
                          const PropertyTable& p_rad, const ColumnOptions& opts);
 
} // namespace rizer
 
#endif