basic_physics.hpp

Go to the documentation of this file.

#ifndef BASIC_PHYSICS_HPP
#define BASIC_PHYSICS_HPP
#include "space_settings.hpp"
#include "constants.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));
}

/* Returns thermal velocity
 * @param[in]   mass     species mass
 * @param[in]   temp_ev  reference temperature
 * return thermal velocity
 */
KOKKOS_INLINE_FUNCTION double thermal_velocity(double mass, double temp_ev){
    return sqrt((temp_ev*EV_2_J)/mass);
}

/* Returns velocity normalized by thermal velocity
 * @param[in]   mass     species mass
 * @param[in]   temp_ev  reference temperature
 * @param[in]   v        input velocity
 * return normalized velocity
 */
KOKKOS_INLINE_FUNCTION double normalize_to_vth(double mass, double temp_ev, double v){
    return v*sqrt(mass/(temp_ev*EV_2_J));
}

/* Returns normalized parallel velocity
 * @param[in]   c_m      species charge over mass
 * @param[in]   mass     species mass
 * @param[in]   B        magnetic field magnitude
 * @param[in]   temp_ev  reference temperature
 * @param[in]   rho      gyroradius
 * return normalized parallel velocity
 */
KOKKOS_INLINE_FUNCTION double normalized_v_para(double c_m, double mass, double B, double temp_ev, double rho){
    return c_m*rho*B*sqrt(mass/(temp_ev*EV_2_J));
}

/* Returns normalized sqrt(mu)
 * @param[in]   B        magnetic field magnitude
 * @param[in]   temp_ev  particle/reference temperature
 * @param[out]  mu       mu
 * return normalized sqrt(mu)
*/
KOKKOS_INLINE_FUNCTION double normalized_sqrt_mu(double B, double temp_ev, double mu){
    return sqrt(mu*2.0*B/(temp_ev*EV_2_J));
}

// 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