TerrainDrag.H

Go to the documentation of this file.

#ifndef TERRAINDRAG_H
#define TERRAINDRAG_H
 
#include "src/core/Physics.H"
#include "src/core/Field.H"
#include "src/CFDSim.H"
 
namespace kynema_sgf::terraindrag {

 
class TerrainDrag : public Physics::Register<TerrainDrag>
{
public:
    static std::string identifier() { return "TerrainDrag"; }
 
    explicit TerrainDrag(CFDSim& sim);
 
    ~TerrainDrag() override = default;
 
    void initialize_fields(int level, const amrex::Geometry& geom) override;
 
    void pre_init_actions() override {}
 
    void post_init_actions() override;
 
    void post_regrid_actions() override;
 
    void pre_advance_work() override;
 
    void post_advance_work() override {}

 
    void convert_waves_to_terrain_fields();
 
    [[nodiscard]] std::string wave_velocity_field_name() const
    {
        return m_wave_velocity_name;
    }
 
    [[nodiscard]] std::string wave_negative_elevation_name() const
    {
        return m_wave_negative_elevation_name;
    }
 
private:
    CFDSim& m_sim;
    const FieldRepo& m_repo;
    const amrex::AmrCore& m_mesh;

    IntField& m_terrain_blank;
    IntField& m_terrain_drag;

    std::string m_terrain_file{"terrain.amrwind"};

    std::string m_roughness_file{"terrain.roughness"};

    amrex::Real m_uniform_z0{0.1_rt};
    Field& m_terrainz0;
    Field& m_terrain_height;
    Field& m_terrain_damping;

    amrex::Real m_damp_east_slope{0.0_rt};
    amrex::Real m_damp_east_full{0.0_rt};
    amrex::Real m_horizontal_tau{20.0_rt};
    amrex::Real m_horizontal_abl_height{1e15_rt};
    amrex::Real m_horizontal_slope_end{1e30_rt};
    amrex::Real m_vertical_slope{0.0_rt};
    amrex::Real m_vertical_full{1e30_rt};
    amrex::Real m_damp_west_slope{0.0_rt};
    amrex::Real m_damp_west_full{0.0_rt};
    amrex::Real m_damp_south_slope{0.0_rt};
    amrex::Real m_damp_south_full{0.0_rt};
    amrex::Real m_damp_north_slope{0.0_rt};
    amrex::Real m_damp_north_full{0.0_rt};

    bool m_terrain_is_waves{false};
    Field* m_wave_volume_fraction{nullptr};
    Field* m_wave_negative_elevation{nullptr};
    const std::string m_wave_volume_fraction_name{"ow_vof"};
    const std::string m_wave_negative_elevation_name{"ow_levelset"};
    const std::string m_wave_velocity_name{"ow_velocity"};
};
} // namespace kynema_sgf::terraindrag
 
#endif