q_h_artificial_dissipation.h

Go to the documentation of this file.

#pragma once
 
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace five_moment
{
namespace limiters
{    
class q_h_artificial_dissipation : public WmApplication
{
public:
    q_h_artificial_dissipation();
 
    ~q_h_artificial_dissipation() override = default;
 
    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 _output_variables;
    }
 
    real numerical_flux(const real* q_l,
                        const real* q_r,
                        const real* aux_l,
                        const real* aux_r,
                        const solverVariables_t* pFV,
                        real* numericalFlux) const override;
 
    real internal_flux(const real* q,
                       const real* aux,
                       const solverVariables_t* pSV,
                       std::vector<std::vector<real>>& internalFlux) const override;
private:
    real _gamma;
    // normalized species mass
    real _mass;
    real _min_density;
    real _min_pressure;
    real _multiplier;
 
    // whether eta penalty should depend on delta x or not (inversely proportional?)
    bool _eta_dx;
 
    std::vector<int> _input_variables;
    std::vector<int> _aux_variables;
    std::vector<int> _output_variables;
};
} // namespace limiters    
} // namespace five_moment
} // namespace apps
} // namespace wxm