PlasmaChannelSolver.h

Go to the documentation of this file.

 
#ifndef RIZER_PLASMA_CHANNEL_SOLVER_H
#define RIZER_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;             

        double rtol = 1.0e-7;   
        double atol_T = 1.0e-3; 
        double atol_Y = 1.0e-9; 

        bool independent_nodes = false;

        std::size_t n_threads = 1;

        std::vector<double> Et_t, Et_v;

        double gap = 0.0;

        std::function<double(double, double)> plasma_voltage;
    };
 
    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<double> rho;
        std::vector<double> p;
        std::vector<double> u_r;
        std::vector<std::string> species;

        double boundary_mass_out = 0.0;
        std::vector<double> boundary_species_mass_out;
        double boundary_total_energy_out = 0.0;
        double boundary_electron_energy_out = 0.0;
        double pressure_work_electron = 0.0;

        double boundary_conduction_energy_out = 0.0;
        double boundary_diffusion_energy_out = 0.0;

        std::vector<double> solver_steps;
        std::vector<double> solver_rhs_evals;
        std::vector<double> solver_jac_evals;
        std::vector<double> solver_err_test_fails;
        std::vector<double> solver_nonlin_conv_fails;
        std::vector<double> solver_last_order;
        std::vector<double> solver_lin_rhs_evals;
        std::vector<double> solver_lin_setups;

        std::vector<double> Rp;       
        std::vector<double> Rp_state; 
        std::vector<double> Vp;       
        std::vector<double> I;        
        std::vector<double> E;        
 
        std::size_t nframes = 0; 
        std::size_t npts = 0;    
        std::size_t nsp = 0;     
    };

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