linear_acoustics.h

Go to the documentation of this file.

#ifndef WXM_APPS_LINEAR_ACOUSTICS_H
#define WXM_APPS_LINEAR_ACOUSTICS_H
 
// Wm includes
#include "apps/wmapplication.h"
 
namespace wxm
{
namespace apps
{
 
//  1d linear advection
class linear_acoustics_t : public WmApplication
{
public:
    linear_acoustics_t();
 
    ~linear_acoustics_t() override;
 
    void setup(const WxCryptSet& wxc) override;
 
    const std::vector<int>& getInputVariableIndexes(int flag) const override
    {
        return _variables;
    }
    const std::vector<int>& getOutputVariableIndexes(int flag) const override
    {
        return _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;
 
protected:
    real _cs;   // sound speed
    real _rho0; // density of fluid
    real _u0;   // average fluid speed
 
    real _K0; // bulk modulus of compressibility -> calculated from
              // _rho0 and _cs
 
    std::vector<int> _variables;
 
private:
    linear_acoustics_t& operator=(const linear_acoustics_t& var) = delete;
    linear_acoustics_t(const linear_acoustics_t& var) = delete;
};
} // namespace apps
} // namespace wxm
 
#endif // WXM_APPS_LINEAR_ACOUSTICS_H