TerrainDrag.H

Go to the documentation of this file.

#ifndef TerrainDrag_H
#define TerrainDrag_H
 
#include "amr-wind/core/Physics.H"
#include "amr-wind/core/Field.H"
#include "amr-wind/CFDSim.H"
 
namespace amr_wind::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();
 
    std::string wave_velocity_field_name() const
    {
        return m_wave_velocity_name;
    }
    std::string wave_velocity_field_name() const {…}
 
    std::string wave_negative_elevation_name() const
    {
        return m_wave_negative_elevation_name;
    }
    std::string wave_negative_elevation_name() const {…}
 
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"};
 
    Field& m_terrainz0;
    Field& m_terrain_height;
 
    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"};
};
class TerrainDrag : public Physics::Register<TerrainDrag> {…};
} // namespace amr_wind::terraindrag
 
#endif