ABLWallFunction.H

Go to the documentation of this file.

#ifndef ABLWALLFUNCTION_H
#define ABLWALLFUNCTION_H
 
#include "src/CFDSim.H"
#include "src/utilities/FieldPlaneAveraging.H"
#include "src/core/FieldBCOps.H"
#include "src/wind_energy/MOData.H"
#include "src/utilities/constants.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf {

class ABLWallFunction
{
public:
    explicit ABLWallFunction(const CFDSim& sim);
 
    ~ABLWallFunction() = default;
 
    MOData& mo() { return m_mo; }
    [[nodiscard]] const MOData& mo() const { return m_mo; }

    [[nodiscard]] amrex::Real utau() const { return m_mo.utau; }

    void init_log_law_height(int max_level);

    void
    update_umean(const VelPlaneAveraging& vpa, const FieldPlaneAveraging& tpa);
 
    void update_tflux(amrex::Real tflux);
 
private:
    const CFDSim& m_sim;
 
    const amrex::AmrCore& m_mesh;

    MOData m_mo;
 
    int m_direction{2};             
    bool m_use_fch{true};           
    amrex::Real m_wall_pos{0.0_rt}; 
    amrex::Vector<amrex::Real> m_gravity{0.0_rt, 0.0_rt, -9.81_rt};

    std::string m_surf_temp_timetable;
    amrex::Vector<amrex::Real> m_surf_temp_time;
    amrex::Vector<amrex::Real> m_surf_temp_value;
 
    bool m_tempflux{true};
    bool m_temp_table{false};
    amrex::Real m_surf_temp_rate{0.0_rt};
    amrex::Real m_surf_temp_rate_tstart{0.0_rt};
    amrex::Real m_surf_temp_init{300.0_rt};
 
    bool m_inflow_outflow{false};
    amrex::Real m_wf_vmag{0.0_rt};
    amrex::Array<amrex::Real, 2> m_wf_vel{0.0_rt, 0.0_rt};
    amrex::Real m_wf_theta{300.0_rt};
};

class ABLVelWallFunc : public FieldBCIface
{
public:
    ABLVelWallFunc(Field& velocity, const ABLWallFunction& wall_func);
 
    void operator()(Field& velocity, FieldState rho_state) override;
 
    template <typename ShearStress>
    void
    wall_model(Field& velocity, FieldState rho_state, const ShearStress& tau);
 
private:
    const ABLWallFunction& m_wall_func;
    std::string m_wall_shear_stress_type{"moeng"};
    std::string m_wall_het_model{"none"};
    amrex::Real m_monin_obukhov_length{constants::LARGE_NUM};
};
 
class ABLTempWallFunc : public FieldBCIface
{
public:
    ABLTempWallFunc(Field& temperature, const ABLWallFunction& wall_fuc);
 
    void operator()(Field& temperature, FieldState rho_state) override;
 
    template <typename HeatFlux>
    void
    wall_model(Field& temperature, FieldState rho_state, const HeatFlux& tau);
 
private:
    const ABLWallFunction& m_wall_func;
    std::string m_wall_shear_stress_type{"moeng"};
    std::string m_wall_het_model{"none"};
    amrex::Real m_monin_obukhov_length{constants::LARGE_NUM};
};
 
} // namespace kynema_sgf
 
#endif /* ABLWALLFUNCTION_H */