basic_physics.hpp

Go to the documentation of this file.

#ifndef BASIC_PHYSICS_HPP
#define BASIC_PHYSICS_HPP
#include "space_settings.hpp"

// Kinetic energy
// E = (1/2) m v^2
KOKKOS_INLINE_FUNCTION double kinetic_energy(double mass, double v){
    return 0.5*mass*v*v;
}

// Return equilibrium distribution function
// f = (n / sqrt(T^3)) * exp(-E/T)
KOKKOS_INLINE_FUNCTION double maxwellian_dist(double den, double temp, double energy){
    const double EXP_LIM = 64.0;
    return den/(temp*sqrt(temp))*exp(-min(EXP_LIM,energy/temp));
}

// gyro_radius
KOKKOS_INLINE_FUNCTION double gyro_radius(double B, double mu, double c2_2m){
    return sqrt(mu/(B*c2_2m));
}

// Convert from {rho, mu} to {energy, pitch}
KOKKOS_INLINE_FUNCTION void rho_mu_to_en_pitch(double B, double c_m, double c2_2m, double mass, double rho, double mu, double& en, double& pitch){
    en = mu*B + c2_2m*(rho*rho*B*B);
    double v_para = c_m*rho*B;
    double v_perp2 = 2.0*mu*B/mass;
    pitch=v_para/sqrt(v_para*v_para + v_perp2);
}

// Convert from {energy, pitch} to {rho, mu}
KOKKOS_INLINE_FUNCTION void en_pitch_to_rho_mu(double B, double c2_2m, double en, double pitch, double& rho, double& mu){
    rho = pitch*sqrt(en/c2_2m)/B;
    mu = max(0.0,(1.0-pitch*pitch)*en/B);
}

#endif