HydroStep.h

Go to the documentation of this file.

 
#ifndef RIZER_HYDRO_STEP_H
#define RIZER_HYDRO_STEP_H
 
#include "models/nrp/BoundaryConditions.h"
#include "thermo/ThermoAccessor.h"
 
#include <cstddef>
#include <vector>
 
namespace rizer
{

    struct HydroStepConfig
    {
        std::vector<double> r_node;
        std::vector<double> r_face;
        std::vector<double> V_cell;
        bool cylindrical_source = true;
 
        OuterBc outer_bc_mode = OuterBc::AmbientGhost;
        double p_amb = -1.0;
        double T_amb = -1.0;
        std::vector<double> Y_amb; 

        double dt_total = 0.0;
        std::size_t n_subcycles = 0;
        bool refuse_above_cfl = false;
        double cfl = 0.9;
        double t0 = 0.0;

        int electron_index = -1;

        const ThermoAccessor *thermo = nullptr;
        GhostStateFn ghost_state;
    };

    struct HydroStepState
    {
        std::vector<double> rho;   
        std::vector<double> u_r;   
        std::vector<double> Tg;    
        std::vector<double> Te;    
        std::vector<double> Y;     
        std::size_t n_species = 0; 
    };

    struct HydroStepResult
    {
        std::size_t n_subcycles_used = 0; 
        double dt_cfl = 0.0;              
        double rho_min = 0.0;
        double p_min = 0.0; 

        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;
    };

    HydroStepResult hydroStep(const HydroStepConfig &cfg, HydroStepState &state);
 
} // namespace rizer
 
#endif // RIZER_HYDRO_STEP_H