EkmanSpiral.H

Go to the documentation of this file.

#ifndef EKMANSPIRAL_H
#define EKMANSPIRAL_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 {

class EkmanSpiral : public Physics::Register<EkmanSpiral>
{
public:
    static std::string identifier() { return "EkmanSpiral"; }
 
    explicit EkmanSpiral(const CFDSim& sim);
 
    ~EkmanSpiral() 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;
 
private:
    const kynema_sgf::SimTime& m_time;
    const FieldRepo& m_repo;
    const amrex::AmrCore& m_mesh;
    Field& m_velocity;
    Field& m_density;
 
    void output_error();

    amrex::Real m_rho{1.0_rt};
    amrex::Real m_DE{0.0_rt};
    amrex::Real m_vel{15.0_rt};

    const int m_w{18};

    const std::string m_output_fname{"ekman_spiral.log"};
};
} // namespace kynema_sgf
 
#endif /* EkmanSpiral_H */