VortexPatchScalarVel.H

Go to the documentation of this file.

#ifndef VortexPatchScalarVel_H
#define VortexPatchScalarVel_H
 
#include "amr-wind/core/Physics.H"
#include "amr-wind/core/Field.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;


 
namespace amr_wind {
 
class VortexPatchScalarVel : public Physics::Register<VortexPatchScalarVel>
{
    static_assert(
        AMREX_SPACEDIM == 3,
        "VortexPatchScalarVel requires 3 dimensional mesh");
 
public:
    static std::string identifier() { return "VortexPatchScalarVel"; }
 
    explicit VortexPatchScalarVel(CFDSim& sim);
 
    ~VortexPatchScalarVel() override = default;

 
    void initialize_fields(int level, const amrex::Geometry& geom) override;
 
    void post_init_actions() override {}
 
    void post_regrid_actions() override {}
 
    void pre_advance_work() override;
 
    void post_advance_work() override;
 
private:
    const CFDSim& m_sim;
 
    Field& m_velocity;
    Field& m_levelset;
    Field& m_density;

    amrex::Vector<amrex::Real> m_loc{0.35_rt, 0.35_rt, 0.35_rt};

    amrex::Real m_radius{0.15_rt};

    amrex::Real m_sfactor{1.0_rt};

    amrex::Real m_TT{3.0_rt};
};
 
} // namespace amr_wind
 
#endif /* VortexPatchScalarVel_H */