bennett_equilibrium.h

Go to the documentation of this file.

#pragma once
 
// WARPXM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace functions
{
namespace mhd
{
class BennettEquilibrium : public WmApplication
{
public:
    BennettEquilibrium();
 
    ~BennettEquilibrium() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _mhd_fluid;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
    void bc_q(const real* q_in,
              const real* aux_in,
              const real* aux_out,
              const solverVariables_t* pFV,
              real* q_out) const override;
 
protected:
    std::vector<int> _mhd_fluid;
 
    std::string _coordinate_system;
    int _ir, _ith, _iz; // r, theta, and z indices
 
    real _gas_gamma; // gamma
    real _Ai;        // mass
    real _Zi;        // charge
    real _T;         // uniform temperature
    real _a;         // pinch radius
    real _I;         // pinch current
    real _wct;       // (omega_c tau)
 
    real _pert;   // perturbation fraction on pressure
    real _lambda; // wavelength of perturbation in z-coordinate
 
private:
    BennettEquilibrium& operator=(const BennettEquilibrium& var);
    BennettEquilibrium(const BennettEquilibrium& var);
};
 
} // namespace mhd
} // namespace functions
} // namespace apps
} // namespace wxm