mhd_two_temperature.h

Go to the documentation of this file.

#pragma once
 
#include "warpxm/warpxm_config.h"
 
#include <limits>
 
namespace wxm
{
namespace apps
{
namespace mhd_two_temperature
{
 
real getRhoFloored(const real* q,
                   const real rho_min = std::numeric_limits<real>::epsilon());
 
real getIonNumberDensity(const real* q,
                         const real Ai,
                         const real rho_min = std::numeric_limits<real>::epsilon());
 
real getIonPressureFloored(const real* q,
                           const real gas_gamma,
                           const real rho_min = std::numeric_limits<real>::epsilon(),
                           const real press_min = std::numeric_limits<real>::epsilon());
 
real getElectronPressureFloored(
    const real* q,
    const real gas_gamma,
    const real rho_min = std::numeric_limits<real>::epsilon(),
    const real press_min = std::numeric_limits<real>::epsilon());
 
real getFluidPressure(const real press_i, const real press_e);
 
void getGradientElectronPressure(const real* grad_q, const real gas_gamma, real* grad_pe);
 
void getGradientIonPressure(const real* q,
                            const real* grad_q,
                            const real gas_gamma,
                            real* grad_press_i,
                            const real rho_min = std::numeric_limits<real>::epsilon());
 
real getElectronTemperatureFloored(const real* q,
                                   const real gas_gamma,
                                   const real Ai,
                                   const real Zi,
                                   const real rho_min,
                                   const real press_min);
 
real getIonTemperatureFloored(const real* q,
                              const real gas_gamma,
                              const real Ai,
                              const real rho_min,
                              const real press_min);
 
void getGradientIonTemperature(const real* q,
                               const real* grad_q,
                               const real gas_gamma,
                               const real Ai,
                               real* grad_Ti,
                               const real rho_min,
                               const real press_min);
 
void getGradientElectronTemperature(const real* q,
                                    const real* grad_q,
                                    const real gas_gamma,
                                    const real Ai,
                                    const real Zi,
                                    real* grad_Te,
                                    const real rho_min,
                                    const real press_min);
 
real current_density_cart(const real skin_depth_norm, const real* grad_B[3], real* J);
 
real current_density_cyl(const real skin_depth_norm,
                         const real r,
                         const real* B,
                         const real* grad_B[3],
                         const int ir,
                         const int ith,
                         const int iz,
                         real* J);
 
} // namespace mhd_two_temperature
} // namespace apps
} // namespace wxm