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"
 
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, 0.0, -9.81};
 
    bool m_tempflux{true};
    amrex::Real m_surf_temp_rate{0.0};
    amrex::Real m_surf_temp_rate_tstart{0.0};
    amrex::Real m_surf_temp_init{300.0};
 
    bool m_inflow_outflow{false};
    amrex::Real m_wf_vmag{0.0};
    amrex::Array<amrex::Real, 2> m_wf_vel{0.0, 0.0};
    amrex::Real m_wf_theta{300.0};
};
 
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"};
};
 
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"};
};
 
} // namespace amr_wind
 
#endif /* ABLWALLFUNCTION_H */