hallmhd.h

Go to the documentation of this file.

#pragma once
 
// wxm includes
#include "apps/mhd/imhd/imhd.h"
 
// stl includes
#include <vector>
#include <limits>
 
namespace wxm
{
namespace apps
{
namespace hallmhd
{
 
real hallmhd_internal_flux(const real gas_gamma,
                           const real skin_depth_norm,
                           const real Ai,
                           const real Zi,
                           const real theta,
                           const solverVariables_t* pSV,
                           const real* q,
                           const real* aux,
                           const real* J,
                           std::vector<std::vector<real>>& internal_flux,
                           const real Ae = 0.0,
                           const real rho_min = std::numeric_limits<real>::epsilon(),
                           const real press_min = std::numeric_limits<real>::epsilon(),
                           const real drift_pe_factor = 1.0,
                           const real jxb_factor = 1.0,
                           const real grad_pe_factor = 1.0);
 
real whistler_wave_speed(const real skin_depth_norm,
                         const real Ai,
                         const real Zi,
                         const solverVariables_t* pSV,
                         const real* q,
                         const real Ae = 0.0,
                         const real rho_min = std::numeric_limits<real>::epsilon());
 
real hallmhd_set_time_step(const solverVariables_t* pSV, const real V_m, const real V_w);
real gyro_dt_limit(const real skin_depth_norm,
                   const real Ai,
                   const real Zi,
                   const solverVariables_t* pSV,
                   const real* q);
 
void calc_wave_speeds(const real gas_gamma,
                      const real skin_depth_norm,
                      const real Ai,
                      const real Zi,
                      const solverVariables_t* pSV,
                      const real* q,
                      const real* J,
                      real& v_plasma,
                      real& v_whistler,
                      real& v_fast_magnetosonic,
                      real& v_hall,
                      const real Ae = 0.0,
                      const real rho_min = std::numeric_limits<real>::epsilon(),
                      const real press_min = std::numeric_limits<real>::epsilon());
 
} // namespace hallmhd
} // namespace apps
} // namespace wxm