AMDNoTherm.H

Go to the documentation of this file.

#ifndef AMDNOTHERM_H
#define AMDNOTHERM_H
 
#include <string>
#include "amr-wind/incflo_enums.H"
#include "amr-wind/turbulence/TurbModelBase.H"
#include "amr-wind/core/FieldRepo.H"
#include "amr-wind/fvm/stencils.H"
#include "AMReX_Config.H"
#include "AMReX_AmrCore.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace amr_wind::turbulence {
template <typename Transport>
class AMDNoTherm : public TurbModelBase<Transport>
{
public:
    static std::string identifier()
    {
        return "AMDNoTherm-" + Transport::identifier();
    }
 
    explicit AMDNoTherm(CFDSim& sim);

    std::string model_name() const override { return "AMDNoTherm"; }

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

    TurbulenceModel::CoeffsDictType model_coeffs() const override;

    void parse_model_coeffs() override;

    void post_advance_work() override {};
 
private:
    amrex::Real m_C{0.333333333333333_rt};
    const Field& m_vel;
    const Field& m_rho;
};
 
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real amd_base_muvel(
    int i,
    int j,
    int k,
    const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> dx, // Grid spacing
    amrex::Real C,                                         // Poincare const
    amrex::Array4<amrex::Real const> const& gradVel) noexcept
{
 
    amrex::Real num_shear = 0;
    amrex::Real denom = 0;
    for (int ii = 0; ii < AMREX_SPACEDIM; ++ii) {
        for (int jj = 0; jj < AMREX_SPACEDIM; ++jj) {
            denom = denom + gradVel(i, j, k, ii * AMREX_SPACEDIM + jj) *
                                gradVel(i, j, k, ii * AMREX_SPACEDIM + jj);
            amrex::Real sij =
                0.5_rt * (gradVel(i, j, k, ii * AMREX_SPACEDIM + jj) +
                          gradVel(i, j, k, jj * AMREX_SPACEDIM + ii));
            for (int kk = 0; kk < AMREX_SPACEDIM; ++kk) {
                amrex::Real dkui = gradVel(i, j, k, ii * AMREX_SPACEDIM + kk);
                amrex::Real dkuj = gradVel(i, j, k, jj * AMREX_SPACEDIM + kk);
                num_shear = num_shear + dkui * dkuj * dx[kk] * dx[kk] * sij;
            }
        }
    }
    denom = amrex::max<amrex::Real>(
        std::numeric_limits<amrex::Real>::epsilon() * 1.0e1_rt, denom);
    num_shear = -C * num_shear;
 
    return amrex::max<amrex::Real>(
        std::numeric_limits<amrex::Real>::epsilon() * 1.0e1_rt,
        num_shear / denom);
}
 
} // namespace amr_wind::turbulence
 
#endif /* AMDNoTherm_H */