briowu_shocktube.h

Go to the documentation of this file.

#ifndef wmapplication_function_mhd_briowu_shocktube_h
#define wmapplication_function_mhd_briowu_shocktube_h
 
// WARPXM includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace functions
{
namespace mhd
{
 
class BrioWuShockTube : public WmApplication
{
public:
    BrioWuShockTube();
 
    ~BrioWuShockTube() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _output_fluid;
    }
 
    void evaluate_function(const real* q,
                           const real* aux,
                           const solverVariables_t* pSV,
                           real* result) const override;
 
protected:
    std::vector<int> _output_fluid;
 
    real _rho_l;
    real _rho_r;
    real _vx_l;
    real _vx_r;
    real _vy_l;
    real _vy_r;
    real _vz_l;
    real _vz_r;
    real _Bx;
    real _By_l;
    real _By_r;
    real _Bz_l;
    real _Bz_r;
    real _P_l;
    real _P_r;
    real _transitionWidth;
 
    real _gas_gamma;
    // real _skin_depth_param;
 
    std::vector<real> _direction;
    std::vector<real> _center;
 
private:
    BrioWuShockTube& operator=(const BrioWuShockTube& var);
    BrioWuShockTube(const BrioWuShockTube& var);
};
} // namespace mhd
} // namespace functions
} // namespace apps
} // namespace wxm
 
#endif