electric_field.h

Go to the documentation of this file.

#pragma once
 
#include "apps/wmapplication.h"
#include "apps/mhd_two_temperature/app_setup_helper.h"
 
// std includes
#include <memory>
#include <vector>
#include <string>
 
namespace wxm
{
namespace apps
{
namespace mhd_two_temperature
{
namespace resistive_mhd
{
class ElectricField : public WmApplication
{
public:
    ElectricField();
 
    ~ElectricField() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return mhd_fluid_indices_;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return mhd_gradient_fluid_indices_;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return electric_field_indices_;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    std::unique_ptr<wxm::apps::mhd_two_temperature::AppSetupHelper> setup_helper_;
 
    std::vector<int> mhd_fluid_indices_;
 
    std::vector<int> mhd_gradient_fluid_indices_;
 
    std::vector<int> electric_field_indices_;
 
    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_floor_;
    real min_pressure_floor_;
    real max_eta_;
    real min_eta_;
    real n_vac_;
    real eta_vac_;
 
    std::string coordinate_system_;
 
    int radial_index_;
    int ir_, ith_, iz_;
 
private:
    ElectricField& operator=(const ElectricField& var);
    ElectricField(const ElectricField& var);
};
 
} // namespace resistive_mhd
} // namespace mhd_two_temperature
} // namespace apps
} // namespace wxm