TransientPlasmaChannelSolver.h

Go to the documentation of this file.

// TransientPlasmaChannelSolver: assembles [Empty1D, PlasmaChannel1D, Empty1D]
// into a Sim1D and advances it with implicit Backward-Euler, recording the
// radial profiles at a chosen cadence. Returns flat numeric arrays.
 
#ifndef RIZER_TRANSIENT_PLASMA_CHANNEL_SOLVER_H
#define RIZER_TRANSIENT_PLASMA_CHANNEL_SOLVER_H
 
#include "PlasmaChannel1D.h"
 
#include <string>
#include <vector>
#include <cstddef>
 
namespace rizer {
 
struct ChannelOptions {
    PlasmaChannelConfig cfg;        
    double dt = 1.0e-9;             
    std::size_t n_steps = 1000;     
    std::size_t record_every = 1;   
    int loglevel = 0;                
    std::string integrator = "newton";
};
 
struct ChannelHistory {
    std::vector<double> t;          
    std::vector<double> r;          
    std::vector<double> Tg;         
    std::vector<double> Te;         
    std::vector<double> ne;         
    std::vector<double> Y;          
    std::vector<std::string> species;  
    std::size_t nframes = 0;  
    std::size_t npts = 0;     
    std::size_t nsp = 0;      
};

ChannelHistory solveChannelTransient(const ChannelOptions& opts);
 
} // namespace rizer
 
#endif