artificial_dissipation.h

Go to the documentation of this file.

#pragma once
 
// wxm includes
#include "apps/wmapplication.h"
#include "apps/5-moment/app_setup_helper.h"
 
// STL includes
#include <memory>
 
namespace wxm
{
namespace apps
{
namespace five_moment
{
namespace limiters
{
 
class ArtificialDissipation : public WmApplication
{
public:
    ArtificialDissipation();
 
    ~ArtificialDissipation() override;
 
    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;
    }
 
    // redefinition of fluxes from parent
 
    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;
 
    // specific calls
 
    real numerical_flux_ldg(const real* q_l,
                            const real* q_r,
                            const real* aux_l,
                            const real* aux_r,
                            const solverVariables_t* pFV,
                            real* numericalFlux) const;
 
    real numerical_flux_ip(const real* q_l,
                           const real* q_r,
                           const real* aux_l,
                           const real* aux_r,
                           const solverVariables_t* pFV,
                           real* numericalFlux) const;
 
protected:
    std::unique_ptr<AppSetupHelper> _setup_helper;
 
    real _min_density;  
    real _min_pressure; 
    real _gas_gamma;    
    real _mass;
 
    // coefficients on artificial momentum and heat diffusivities
    real _coeff_dens;
    real _coeff_mom;
    real _coeff_heat;
    real _diff_cap;
 
    real _cfl_diff;
 
    bool _is_cyl;
    bool _use_v_anom;
 
    int _ir;
    int _ith;
    int _iz;
 
    real _dx_cut;
 
    // gradient calculation method ("ip" or "ldg")
    std::string _gradient_method;
    std::vector<real> _params;
 
    std::vector<int> _input_variables;
    std::vector<int> _aux_variables;
    std::vector<int> _output_variables;
 
    // flux tensor is used to calculate left and right fluxes for numerical flux
    // mutable because this will continuously be changed upon call, but functions are
    // labled const.
    // this is just a helper vector so this behavior should be justified.
    mutable std::vector<std::vector<real>> _flux_tensor;
};
 
} // namespace limiters
} // namespace five_moment
} // namespace apps
} // namespace wxm