Kosovic.H

Go to the documentation of this file.

#ifndef KOSOVIC_H
#define KOSOVIC_H
 
#include <string>
#include "src/turbulence/TurbModelBase.H"
#include "src/core/FieldRepo.H"
#include "src/utilities/constants.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::turbulence {

template <typename Transport>
class Kosovic : public TurbModelBase<Transport>
{
public:
    static std::string identifier()
    {
        return "Kosovic-" + Transport::identifier();
    }
 
    explicit Kosovic(CFDSim& sim);

    [[nodiscard]] std::string model_name() const override { return "Kosovic"; }

    void update_turbulent_viscosity(
        FieldState fstate, DiffusionType /*unused*/) override;

    void post_advance_work() override {}

    void update_alphaeff(Field& alphaeff) override;

    void parse_model_coeffs() override;

    [[nodiscard]] TurbulenceModel::CoeffsDictType model_coeffs() const override;
 
private:
    // Ref: Mirocha et. al "Implementation of a Nonlinear Subfilter Turbulence
    //  Stress Model for Large-Eddy Simulation in the Advanced Research WRF
    //  Model"
    // , MWR 2012.
    amrex::Real m_Cb{0.36_rt};
    amrex::Real m_Cs{0.135_rt};
    amrex::Real m_C1{2.1_rt};
    amrex::Real m_C2{2.1_rt};
    amrex::Real m_Sk{0.5_rt};
    bool m_surfaceRANS{false};
    amrex::Real m_surfaceFactor{0.0_rt};
    amrex::Real m_switchLoc{24.0_rt};
    amrex::Real m_surfaceRANSExp{2.0_rt};
    amrex::Real m_LESTurnOff{constants::LARGE_NUM};
    amrex::Real m_muCoeff{3.0_rt};
    bool m_writeTerms{false};
    const Field& m_vel;
    const Field& m_rho;
    const Field& m_theta;
    Field& m_Nij;
    Field& m_divNij;
    std::string m_wall_het_model{"none"};
    amrex::Real m_monin_obukhov_length{constants::LARGE_NUM};
    amrex::Real m_kappa{0.41_rt};
    amrex::Real m_gamma_m{5.0_rt};
    amrex::Real m_beta_m{16.0_rt};
    amrex::Real m_surface_roughness_z0{0.1_rt};
    amrex::Vector<amrex::Real> m_gravity{0.0_rt, 0.0_rt, -9.81_rt};
};
 
} // namespace kynema_sgf::turbulence
 
#endif /* SMAGORINSKY_H */