ABLWallFunction.H

Go to the documentation of this file.

#ifndef ABLWALLFUNCTION_H
#define ABLWALLFUNCTION_H
 
#include "amr-wind/CFDSim.H"
#include "amr-wind/utilities/FieldPlaneAveraging.H"
#include "amr-wind/utilities/FieldPlaneAveragingFine.H"
#include "amr-wind/core/FieldBCOps.H"
#include "amr-wind/wind_energy/MOData.H"
#include "amr-wind/utilities/constants.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace amr_wind {

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

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

    void init_log_law_height();

    void update_umean(
        const VelPlaneAveragingFine& vpa, const FieldPlaneAveragingFine& tpa);
 
    void update_tflux(const 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::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, const FieldState rho_state) override;
 
    template <typename ShearStress>
    void wall_model(
        Field& velocity, const 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, const FieldState rho_state) override;
 
    template <typename HeatFlux>
    void wall_model(
        Field& temperature, const 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 amr_wind
 
#endif /* ABLWALLFUNCTION_H */