RadialDiffusion.h

Go to the documentation of this file.

 
#ifndef RIZER_RADIAL_DIFFUSION_H
#define RIZER_RADIAL_DIFFUSION_H
 
#include "numeric/PropertyTable.h"
 
#include "cantera/base/ctexceptions.h"
 
#include <cstddef>
#include <vector>
 
namespace rizer
{
 
    class PlasmaChannel1D;

    struct RadialDiffusionCoefficients
    {
        const PropertyTable *kappa = nullptr;
        const PropertyTable *kappa_e = nullptr;
        double species_diffusivity = 0.0;
        double ambient_temperature = 0.0;
        double wall_ghost_density = 0.0;
        const std::vector<double> *Y_amb = nullptr;
        const std::vector<double> *molecular_weights = nullptr;
        const std::vector<double> *species_cp = nullptr;
        std::size_t electron_index = 0;
        std::size_t n_species = 0;
    };

    void addRadialDiffusion(const PlasmaChannel1D &domain,
                            Cantera::span<const double> x,
                            const RadialDiffusionCoefficients &c, std::size_t j,
                            double lhsTg, double lhsTe, double &dTg, double &dTe,
                            double *dYdt, double *scratch);

    struct WallEnergyOutflow
    {
        double conduction = 0.0;
        double diffusion = 0.0;
        std::vector<double> species_mass;
    };

    WallEnergyOutflow wallEnergyOutflow(const PlasmaChannel1D &domain,
                                        Cantera::span<const double> x,
                                        const RadialDiffusionCoefficients &c,
                                        const double *species_enthalpy);
 
} // namespace rizer
 
#endif // RIZER_RADIAL_DIFFUSION_H