artificial_dissipation.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "apps/wmapplication.h"
#include "apps/mhd_two_temperature/app_setup_helper.h"
 
// std includes
#include <memory>
 
namespace wxm
{
namespace apps
{
namespace mhd_two_temperature
{
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 q_indexes_;
    }
 
    const std::vector<int>& getAuxiliaryVariableIndexes(int flag) const override
    {
        return grad_q_indexes_;
    }
 
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return q_indexes_;
    }
 
    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<wxm::apps::mhd_two_temperature::AppSetupHelper> setup_helper_;
 
    real gas_gamma_;
 
    real Ai_;
 
    real Zi_;
 
    real min_density_floor_;
 
    real min_pressure_floor_;
 
    real coeff_dens_;
 
    real coeff_mom_;
 
    real coeff_ion_heat_;
 
    real coeff_electron_heat_;
 
    real diff_cap_;
 
    bool use_v_anom_;
 
    real cfl_diff_;
 
    real dx_cut_;
 
    bool is_cyl_;
 
    int ir_;
    int ith_;
    int iz_;
 
    std::string gradient_method_;
 
    std::vector<int> q_indexes_;
    std::vector<int> grad_q_indexes_;
 
    mutable std::vector<std::vector<real>> flux_tensor_;
};
 
} // namespace limiters
} // namespace mhd_two_temperature
} // namespace apps
} // namespace wxm