euler_vortex.h

#ifndef __EULER_VORTEX_H__
#define __EULER_VORTEX_H__

#include "dg/dg_base.hpp"
#include "parameters/all_parameters.h"
#include "physics/euler.h"
#include "tests.h"

namespace PHiLiP {
namespace Tests {

// Masatsuka2018 section 6.3 & section 7.13.3

template <int dim, typename real>
class EulerVortexFunction : public dealii::Function<dim,real>
{
public:

    EulerVortexFunction (
        const Physics::Euler<dim, dim+2, real> euler_physics,
        const dealii::Point<dim> initial_vortex_center,
        const real vortex_strength,
        const real vortex_stddev_decay);

    const Physics::Euler<dim, dim+2, real> euler_physics; 
    const real vortex_characteristic_length; 
    const dealii::Point<dim> initial_vortex_center; 
    const real vortex_strength; 
    const real vortex_stddev_decay; 


    real value (const dealii::Point<dim> &point, const unsigned int istate = 0) const;

private:
    dealii::Point<dim> advected_location(const dealii::Point<dim> old_location) const;

};

template <int dim, int nstate>
class EulerVortex: public TestsBase
{
public:
    EulerVortex () = delete;

    explicit EulerVortex(const Parameters::AllParameters *const parameters_input);


    int run_test () const;

};


} // Tests namespace
} // PHiLiP namespace
#endif