wmapplication_resistive_diffusion.h

Go to the documentation of this file.

#pragma once
 
// WarpM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
 
class WmApplication_Resistive_Diffusion : public WmApplication
{
public:
    WmApplication_Resistive_Diffusion();
 
    ~WmApplication_Resistive_Diffusion() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _inputVariables;
    }
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _auxVariables;
    }
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _outputVariables;
    }
 
    // redefinition of fluxes from parent
    real numerical_flux(const real* q_l,
                        const real* q_r,
                        const real* aux_l,
                        const real* aux_r,
                        const solverVariables_t* pFV,
                        real* numericalFlux) const override;
    // real internal_flux(const real * q, const real * aux, const
    // elementGeometry_t * pEG, real *internalFlux[8]) const;
    real internal_flux(const real* q,
                       const real* aux,
                       const solverVariables_t* pSV,
                       std::vector<std::vector<real>>& internalFlux) const override;
 
    // specific calls
    real numerical_flux_ldg(const real* q_l,
                            const real* q_r,
                            const real* aux_l,
                            const real* aux_r,
                            const solverVariables_t* pFV,
                            real* numericalFlux) const;
 
    real numerical_flux_ip(const real* q_l,
                           const real* q_r,
                           const real* aux_l,
                           const real* aux_r,
                           const solverVariables_t* pFV,
                           real* numericalFlux) const;
 
protected:
    real _gas_gamma;
    real _skin_depth_norm;
    real _nu_p_norm;
    real _omega_p_norm;
    real _Ai;
    real _Zi;
    real _lnlam;
 
    real _constant_resistivity_eta;
    bool _const_resisitivity;
 
    // std::string _imhd_numerical_flux_type; // this is the type of numerical flux you
    // want to use (rusanov, roe, hll, hlld, etc.)
 
    std::vector<int> _inputVariables;
    std::vector<int> _auxVariables;
    std::vector<int> _outputVariables;
 
    std::string _resistivity_type;
    std::string _gradient_method;
 
    real _min_density;
    real _min_pressure;
    real _max_eta;
    real _min_eta;
 
    real _n_vac;
    real _eta_vac;
 
    real _cfl_diff;
    
    std::string _coordinate_system;
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
    // flux tensor is used to calculate left and right fluxes for numerical flux
    // mutable because this will continuously be changed upon call, but functions are
    // labled const.
    // this is just a helper vector so this behavior should be justified.
    mutable std::vector<std::vector<real>> _flux_tensor;
 
private:
    WmApplication_Resistive_Diffusion&
    operator=(const WmApplication_Resistive_Diffusion& var);
    WmApplication_Resistive_Diffusion(const WmApplication_Resistive_Diffusion& var);
};
 
} // namespace rmhd
} // namespace apps
} // namespace wxm