WallFunction.H

Go to the documentation of this file.

#ifndef WALLFUNCTION_H
#define WALLFUNCTION_H
 
#include "amr-wind/CFDSim.H"
#include "amr-wind/core/FieldBCOps.H"
#include "amr-wind/utilities/FieldPlaneAveragingFine.H"
#include "amr-wind/boundary_conditions/wall_models/LogLaw.H"
#include "amr-wind/boundary_conditions/wall_models/MOSD.H"
 
namespace amr_wind {

class WallFunction
{
public:
    explicit WallFunction(CFDSim& sim);
 
    amrex::Real utau() const { return m_log_law.utau_mean; }
    LogLaw log_law() const { return m_log_law; }
    MOSD mosd() const { return m_mosd; }

    void update_umean();
    void update_utau_mean();
    void update_time();
 
    ~WallFunction() = default;
 
private:
    const CFDSim& m_sim;
 
    const amrex::AmrCore& m_mesh;

    LogLaw m_log_law;
    int m_direction{2}; 
    VelPlaneAveragingFine m_pa_vel;
 
    MOSD m_mosd;
};

class VelWallFunc : public FieldBCIface
{
public:
    VelWallFunc(Field& velocity, WallFunction& wall_func);
 
    void operator()(Field& velocity, const FieldState rho_state) override;
 
    static void wall_model(
        Field& velocity, const FieldState rho_state, const amrex::Real utau);
 
    template <typename ShearStressSimple>
    static void wall_model(
        Field& velocity,
        const FieldState rho_state,
        const ShearStressSimple& tau);
 
private:
    WallFunction& m_wall_func;
    std::string m_wall_shear_stress_type{"constant"};
};
} // namespace amr_wind
 
#endif /* WALLFUNCTION_H */