CollisionModel.h

Go to the documentation of this file.

 
#ifndef RIZER_COLLISION_MODEL_H
#define RIZER_COLLISION_MODEL_H
 
#include "numeric/PropertyTable.h"
 
#include <cstddef>
#include <vector>
 
namespace rizer {
 
class CollisionModel {
public:
    struct SpeciesSpec {
        std::vector<double> energy_J;   
        std::vector<double> sigma_m2;   
        double radius = -1.0;           
        int Z = 0;                      
        bool is_electron = false;       
        double molar_mass = 0.0;        
    };
 
    CollisionModel() = default;
    CollisionModel(const std::vector<SpeciesSpec>& specs, double Te_min, double Te_max,
                   std::size_t Te_n);
 
    bool empty() const { return m_species.empty(); }

    void evaluate(double Te, double n_e, const double* n_k, std::size_t nsp,
                  double& sigma, double& nu_E, double& nu_eH, double& nu_eI,
                  double& nu_eH_mass_weighted) const;

    static double electronElectronCollisionFrequency(double n_e, double Te);

    static double meanCrossSection(const std::vector<double>& energy_J,
                                    const std::vector<double>& sigma_m2, double Te,
                                    double x_max = 20.0, int N = 1000);

    static std::vector<double> meanCrossSectionTable(
        const std::vector<double>& energy_J, const std::vector<double>& sigma_m2,
        const std::vector<double>& Te, double x_max = 20.0, int N = 1000);
 
private:
    enum class SpeciesCollisionKind { Skip, Tabulated, HardSphere, Ion };
    struct SpeciesCollisionConfig {
        SpeciesCollisionKind kind = SpeciesCollisionKind::Skip;
        PropertyTable Qbar;          
        double cst = 0.0;            
        int Z = 0;
        double two_me_over_mh = 0.0; 
        double inv_mh = 0.0;         
    };
    std::vector<SpeciesCollisionConfig> m_species;
};
 
} // namespace rizer
 
#endif