olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM > Struct Template Reference

#include <collisionLES.h>

Collaboration diagram for olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >:

Public Types
using	MomentaF = typename MOMENTA::template type<DESCRIPTOR>
using	CollisionO = typename COLLISION::template type<DESCRIPTOR, MOMENTA, EQUILIBRIUM>

Public Member Functions
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
V	computeEffectiveOmega (CELL &cell, PARAMETERS &parameters) any_platform
template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>
CellStatistic< V >	apply (CELL &cell, PARAMETERS &parameters) any_platform

Detailed Description

template<typename COLLISION, typename DESCRIPTOR, typename MOMENTA, typename EQUILIBRIUM>
struct olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 81 of file collisionLES.h.

Member Typedef Documentation

CollisionO

template<typename COLLISION , typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::CollisionO = typename COLLISION::template type<DESCRIPTOR, MOMENTA, EQUILIBRIUM>

Definition at line 83 of file collisionLES.h.

MomentaF

template<typename COLLISION , typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

using olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::MomentaF = typename MOMENTA::template type<DESCRIPTOR>

Definition at line 82 of file collisionLES.h.

Member Function Documentation

apply()

template<typename COLLISION , typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>

CellStatistic< V > olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 107 of file collisionLES.h.

                                                                          {
    parameters.template set<descriptors::OMEGA>(
      computeEffectiveOmega(cell, parameters));
    const auto iT = parameters.template get<descriptors::LATTICE_TIME>();
    V piNeqNormSqr { };
    MomentaF().computePiNeqNormSqr(cell, piNeqNormSqr);
    V avShear = cell.template getField<descriptors::AV_SHEAR>();
    avShear = (avShear*iT + util::sqrt(piNeqNormSqr)) / (iT+1);
    cell.template setField<descriptors::AV_SHEAR>(avShear);
    return CollisionO().apply(cell, parameters);
  }

References olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::computeEffectiveOmega(), and olb::util::sqrt().

Here is the call graph for this function:

computeEffectiveOmega()

template<typename COLLISION , typename DESCRIPTOR , typename MOMENTA , typename EQUILIBRIUM >

template<typename CELL , typename PARAMETERS , typename V = typename CELL::value_t>

V olb::collision::detail::ShearSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::computeEffectiveOmega ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 86 of file collisionLES.h.

                                                                           {
    V piNeqNormSqr { };
    MomentaF().computePiNeqNormSqr(cell, piNeqNormSqr);
    const V rho = MomentaF().computeRho(cell);
    const auto iT = parameters.template get<descriptors::LATTICE_TIME>();
    const V omega = parameters.template get<descriptors::OMEGA>();
    const V smagorinsky = parameters.template get<collision::LES::Smagorinsky>();
    V avShear = cell.template getField<descriptors::AV_SHEAR>();
    V piNeqNorm = util::sqrt(piNeqNormSqr);
    V preFactor = smagorinsky*smagorinsky
                * descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
                * 2 * util::sqrt(2);
    avShear = (avShear*iT + piNeqNorm) / (iT+1);
    V tauMol = V{1} / omega;
    V piNeqNormSISM = piNeqNorm - avShear;
    V tauTurb = V{0.5} * (util::sqrt(tauMol*tauMol+(preFactor*piNeqNormSISM/rho))-tauMol);
    V tauEff = tauMol + tauTurb;
    return  V{1} / tauEff;
  }

References olb::util::sqrt().

Here is the call graph for this function:

Here is the caller graph for this function:

---

The documentation for this struct was generated from the following file:

• src/dynamics/collisionLES.h