pn_bc_gradients.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
namespace bc
{
/*
 * Sets gradient BC for plasma-neutral model for acceleration region inlet in axisymmetric
 * zpinch
 * Current implementation only works for single-fluid MHD. Plasma-neutral model has
 * not been implemented.
 *
 */
class pn_bc_gradients : public WmApplication
{
public:
    pn_bc_gradients();
 
    ~pn_bc_gradients() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _var_idcs;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return _aux_var_idcs;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _var_idcs;
    }
 
    // redefinition of bc
    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:
    int _bc_type;
    real _gas_gamma;
    real _min_density;
 
    std::vector<int> _var_idcs;
    std::vector<int> _aux_var_idcs;
 
private:
    pn_bc_gradients& operator=(const pn_bc_gradients& var);
    pn_bc_gradients(const pn_bc_gradients& var);
};
 
} // namespace bc
} // namespace rmhd
} // namespace apps
} // namespace wxm