OneEqKsgs.H

Go to the documentation of this file.

#ifndef ONEEQKSGS_H
#define ONEEQKSGS_H
 
#include <string>
#include "amr-wind/turbulence/TurbModelBase.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace amr_wind::turbulence {

template <typename Transport>
class OneEqKsgs : public TurbModelBase<Transport>
{
public:
    explicit OneEqKsgs(CFDSim& sim);
 
    ~OneEqKsgs() override;
 
protected:
    Field& m_vel;
 
    Field& m_turb_lscale;
    Field& m_shear_prod;
    Field& m_buoy_prod;
    Field& m_dissip;
    Field& m_rho;
 
    Field* m_tke{nullptr};

    amrex::Real m_Ce{0.1_rt};
    amrex::Real m_Ceps{0.93_rt};
};

template <typename Transport>
class OneEqKsgsM84 : public OneEqKsgs<Transport>
{
public:
    static std::string identifier()
    {
        return "OneEqKsgsM84-" + Transport::identifier();
    }
 
    explicit OneEqKsgsM84(CFDSim& sim);
 
    ~OneEqKsgsM84() override;
 
    std::string model_name() const override { return "OneEqKsgsM84"; }

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

    void post_advance_work() override;

    void update_alphaeff(Field& alphaeff) override;

    void update_scalar_diff(Field& deff, const std::string& name) override;

    void parse_model_coeffs() override;

    TurbulenceModel::CoeffsDictType model_coeffs() const override;
 
private:
    Field& m_temperature;

    amrex::Vector<amrex::Real> m_gravity{0.0_rt, 0.0_rt, -9.81_rt};

    bool m_hybrid_rl{false};
    Field* m_sdr;
};

template <typename Transport>
class OneEqKsgsS94 : public OneEqKsgs<Transport>
{
public:
    static std::string identifier()
    {
        return "OneEqKsgsS94-" + Transport::identifier();
    }
 
    explicit OneEqKsgsS94(CFDSim& sim);
 
    ~OneEqKsgsS94() override;
 
    std::string model_name() const override { return "OneEqKsgsS94"; }

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

    void post_advance_work() override {}

    void update_scalar_diff(Field& deff, const std::string& name) override;

    void parse_model_coeffs() override;

    TurbulenceModel::CoeffsDictType model_coeffs() const override;
};
 
} // namespace amr_wind::turbulence
 
#endif /* ONEEQKSGS_H */