rmhd.h

Go to the documentation of this file.

#pragma once
 
#include "apps/mhd/imhd/imhd.h"
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
 
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);
 
real get_eta(const real* q,
             const real J[3],
             const bool const_resistivity,
             const real constant_resistivity_eta,
             const std::string resistivity_type,
             const real max_eta,
             const real min_eta,
             const real n_vac,
             const real eta_vac,
             const real gas_gamma,
             const real ompt,
             const real nupt,
             const real lnlam = 10.0,
             const real Ai = 1.0,
             const real Zi = 1.0,
             const real Ae = 1.0 / 1836.0,
             const real rho_min = std::numeric_limits<real>::epsilon(),
             const real press_min = std::numeric_limits<real>::epsilon());
 
real get_resistivity(const real* q,
                     const real gas_gamma,
                     const real lnlam = 10,
                     const real Ai = 1.0,
                     const real Zi = 1.0,
                     const real Ae = 1.0 / 1836.0,
                     const real rho_min = std::numeric_limits<real>::epsilon(),
                     const real press_min = std::numeric_limits<real>::epsilon());
 
real get_resistivity_chodura(const real* q,
                             const real J[3],
                             const real gas_gamma,
                             const real ompt,
                             const real nupt,
                             const real lnlam = 10,
                             const real Ai = 1.0,
                             const real Zi = 1.0,
                             const real Ae = 1.0 / 1836.0,
                             const real rho_min = std::numeric_limits<real>::epsilon(),
                             const real press_min = std::numeric_limits<real>::epsilon());
 
real rmhd_internal_diffusion_flux(
    const real nu_p_norm,
    const real skin_depth_norm,
    const solverVariables_t* pSV,
    const real* q,
    const real* J,
    const real eta,
    std::vector<std::vector<real>>& internal_diffusion_flux);
 
void noslipWallCondition_MHD(const real gas_gamma,
                             const solverVariables_t* pSV,
                             const real* QC,
                             const real vx_w,
                             const real vy_w,
                             const real vz_w,
                             const real bx_w,
                             const real by_w,
                             const real bz_w,
                             real* QC_w,
                             const real rho_min = std::numeric_limits<real>::epsilon(),
                             const real press_min = std::numeric_limits<real>::epsilon());
 
real set_time_step(const solverVariables_t* pSV,
                   const real j_coefficient,
                   const real diffusion_coefficient,
                   const real eta,
                   const real cfl_diff = 0.5);
 
namespace bc
{
// helper functions for boundary condition current profiles
real ramp(const real t, const real per, const real Imax);
real staggered_ramp(const real t, const real per, const real Imax);
real half_period_sine_wave(const real t, const real t_qtr, const real Imax);
 
} // namespace bc
 
} // namespace rmhd
} // namespace apps
} // namespace wxm