density_diffusion_flux.h

Go to the documentation of this file.

#pragma once
 
// Wm includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
namespace rmhd
{
 
class DensityDiffusionFlux : public WmApplication
{
public:
    DensityDiffusionFlux();
 
    ~DensityDiffusionFlux() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    enum GradientMethod
    {
        LDG
    };
 
    enum DDiffType
    {
        BASIC,
        NIMROD
    };
 
    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 _input_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;
 
    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;
 
protected:
    GradientMethod _gradient_method; 
    DDiffType _ddiff_type; 
    std::vector<int>
        _input_variables; 
    std::vector<int>
        _aux_variables; 
    real _ddiff;       
    real _ddiff_limit; 
    real _cfl_diff;    
    mutable std::vector<std::vector<real>> _flux_tensor;
 
private:
    DensityDiffusionFlux& operator=(const DensityDiffusionFlux& var) = delete;
    DensityDiffusionFlux(const DensityDiffusionFlux& var) = delete;
};
} // namespace rmhd
} // namespace apps
} // namespace wxm