PlasmaState.h

Go to the documentation of this file.

 
#ifndef RIZER_PLASMA_STATE_H
#define RIZER_PLASMA_STATE_H
 
#include "units.h"
 
#include "cantera/thermo/PlasmaPhase.h"
 
namespace rizer
{

    constexpr double TemperatureFloor = 200.0;

    constexpr double TemperatureCeiling = 1.0e5;

    constexpr double GasConstant = units::R_kmol;
    constexpr double Boltzmann = units::k_b;

    inline double monatomicElectronCv(double electron_molecular_weight)
    {
        return 1.5 * GasConstant / electron_molecular_weight;
    }

    struct TwoTemperatureState
    {
        double pressure = 0.0;          
        double electron_pressure = 0.0; 
        double heavy_pressure = 0.0;    
        double electron_fraction = 0.0; 
        double mean_temperature = 0.0;  
        double electron_density = 0.0;  
    };

    inline TwoTemperatureState readTwoTemperatureState(const Cantera::PlasmaPhase &phase)
    {
        const std::size_t electron = phase.electronSpeciesIndex();
        TwoTemperatureState state;
        state.pressure = phase.pressure();
        state.electron_pressure = phase.electronPressure();
        state.heavy_pressure = state.pressure - state.electron_pressure;
        state.electron_fraction = phase.moleFraction(electron);
        state.mean_temperature = phase.meanTemperature();
        // concentration() is in kmol/m^3 and Avogadro's number is per kmol, so
        // the product is the number density the collision model wants.
        state.electron_density = phase.concentration(electron) * Cantera::Avogadro;
        return state;
    }
 
} // namespace rizer
 
#endif // RIZER_PLASMA_STATE_H