ConvectingTaylorVortex.H

Go to the documentation of this file.

#ifndef CONVECTINGTAYLORVORTEX_H
#define CONVECTINGTAYLORVORTEX_H
 
#include "src/core/Physics.H"
#include "src/core/Field.H"
#include "src/CFDSim.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::ctv {
 
namespace {
struct UExact
{
    AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real operator()(
        amrex::Real /*u0*/,
        amrex::Real /*v0*/,
        amrex::Real /*omega*/,
        amrex::Real /*x*/,
        amrex::Real /*y*/,
        amrex::Real /*t*/) const;
    const int m_comp{0};
};
 
struct VExact
{
    AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real operator()(
        amrex::Real /*u0*/,
        amrex::Real /*v0*/,
        amrex::Real /*omega*/,
        amrex::Real /*x*/,
        amrex::Real /*y*/,
        amrex::Real /*t*/) const;
    const int m_comp{1};
};
 
struct WExact
{
    AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real operator()(
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/) const;
    const int m_comp{2};
};
 
struct GpxExact
{
    AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real operator()(
        amrex::Real /*u0*/,
        amrex::Real /*unused*/,
        amrex::Real /*omega*/,
        amrex::Real /*x*/,
        amrex::Real /*unused*/,
        amrex::Real /*t*/) const;
    const int m_comp{0};
};
 
struct GpyExact
{
    AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real operator()(
        amrex::Real /*unused*/,
        amrex::Real /*v0*/,
        amrex::Real /*omega*/,
        amrex::Real /*unused*/,
        amrex::Real /*y*/,
        amrex::Real /*t*/) const;
    const int m_comp{1};
};
 
struct GpzExact
{
    AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real operator()(
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/,
        amrex::Real /*unused*/) const;
    const int m_comp{2};
};
 
} // namespace

class ConvectingTaylorVortex : public Physics::Register<ConvectingTaylorVortex>
{
public:
    static std::string identifier() { return "ConvectingTaylorVortex"; }
 
    explicit ConvectingTaylorVortex(const CFDSim& sim);
 
    ~ConvectingTaylorVortex() override = default;
 
    void initialize_fields(int level, const amrex::Geometry& geom) override;
 
    template <typename T>
    amrex::Real compute_error(const Field& /*field*/);
 
    void post_init_actions() override;
 
    void post_regrid_actions() override {}
 
    void pre_advance_work() override {}
 
    void post_advance_work() override;
 
    [[nodiscard]] amrex::Real get_u0() const { return m_u0; }
 
    [[nodiscard]] amrex::Real get_v0() const { return m_v0; }
 
    [[nodiscard]] amrex::Real get_omega() const { return m_omega; }
 
private:
    const kynema_sgf::SimTime& m_time;
    const kynema_sgf::CFDSim& m_sim;
    const FieldRepo& m_repo;
    const amrex::AmrCore& m_mesh;
 
    Field& m_velocity;
    Field& m_gradp;
    Field& m_density;
 
    void output_error();

    amrex::Real m_rho{1.0_rt};

    amrex::Real m_u0{1.0_rt};

    amrex::Real m_v0{1.0_rt};

    amrex::Real m_omega;

    const int m_w{18};

    std::string m_output_fname{"ctv.log"};
 
    bool m_activate_pressure{false};
    bool m_mesh_mapping{false};
};
} // namespace kynema_sgf::ctv
 
#endif /* ConvectingTaylorVortex_H */