Smagorinsky.H

Go to the documentation of this file.

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

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

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

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

    void post_advance_work() override {}

    void parse_model_coeffs() override;

    [[nodiscard]] TurbulenceModel::CoeffsDictType model_coeffs() const override;
 
private:
    amrex::Real m_Cs{0.135_rt};
 
    const Field& m_vel;
    const Field& m_rho;
};
 
} // namespace amr_wind::turbulence
 
#endif /* SMAGORINSKY_H */