electric_field.h

Go to the documentation of this file.

#pragma once
 
// General includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
 
class ElectricField : public WmApplication
{
public:
    ElectricField();
 
    ~ElectricField() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _fluid_idcs;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _gradient_fluid_idcs;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _electric_field_idcs;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    std::vector<int> _fluid_idcs;          // conserved variable vector to output to
    std::vector<int> _gradient_fluid_idcs; // conserved variable vector to output to
    std::vector<int> _electric_field_idcs; // conserved variable vector to output to
    real _gas_gamma;
    real _skin_depth_norm;
    real _nu_p_norm;
    real _omega_p_norm;
    real _Ai;
    real _Zi;
    real _lnlam;
 
    std::string _resistivity_type;
    real _constant_resistivity_eta;
    bool _const_resisitivity;
    real _min_density;
    real _min_pressure;
    real _max_eta;
    real _min_eta;
    real _n_vac;
    real _eta_vac;
 
    // if cylindrical curlB has to be calculated differently
    std::string _coordinate_system;
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
private:
    ElectricField& operator=(const ElectricField& var);
    ElectricField(const ElectricField& var);
};
 
} // namespace rmhd
} // namespace apps
} // namespace wxm