cyl_source.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "apps/wmapplication.h"
 
// STL includes
#include <memory>
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
class CylSource : public WmApplication
{
public:
    CylSource();
 
    ~CylSource();
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _input_variables;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _aux_variables;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _input_variables;
    }
 
    real source(const real* q,
                const real* aux,
                const elementGeometry_t* pEG,
                real* source) const override;
 
protected:
    real _gas_gamma;
    real _skin_depth_norm;
    real _nu_p_norm;
    real _omega_p_norm;
    real _Ai;
    real _Zi;
    real _lnlam;
    int _radial_index;
    int _ir, _ith, _iz; // r, theta, and z indices
 
    real _constant_resistivity_eta;
    bool _const_resisitivity;
    real _min_density_floor;
    real _min_pressure_floor;
    real _max_eta;
    real _min_eta;
    std::string _resistivity_type;
 
    real _n_vac;
    real _eta_vac;
 
    std::vector<int> _input_variables;
    std::vector<int> _aux_variables;
 
    // pointer to subapplications and boolean
    std::unique_ptr<WmApplication> _floor_app;
    bool _has_floor_app;
 
private:
    CylSource& operator=(const CylSource& var) = delete;
    CylSource(const CylSource& var) = delete;
};
} // namespace rmhd
} // namespace apps
} // namespace wxm