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

#include <collisionLES.h>

Collaboration diagram for olb::collision::detail::ConStrainSmagorinskyEffectiveOmega< 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<concepts::Cell CELL, concepts::Parameters PARAMETERS, typename V = typename CELL::value_t>
V	computeEffectiveOmega (CELL &cell, PARAMETERS &parameters) any_platform
template<concepts::Cell CELL, concepts::Parameters 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::ConStrainSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 247 of file collisionLES.h.

Member Typedef Documentation

CollisionO

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

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

Definition at line 249 of file collisionLES.h.

MomentaF

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

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

Definition at line 248 of file collisionLES.h.

Member Function Documentation

apply()

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

template<concepts::Cell CELL, concepts::Parameters PARAMETERS, typename V = typename CELL::value_t>

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

inline

Definition at line 283 of file collisionLES.h.

                                                                         {
    parameters.template set<descriptors::OMEGA>(
      computeEffectiveOmega(cell, parameters));
    return CollisionO().apply(cell, parameters);
  }

References olb::collision::detail::ConStrainSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::computeEffectiveOmega().

Here is the call graph for this function:

computeEffectiveOmega()

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

template<concepts::Cell CELL, concepts::Parameters PARAMETERS, typename V = typename CELL::value_t>

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

inline

Turbulent realaxation time

Effective realaxation time

Definition at line 252 of file collisionLES.h.

                                                                           {
    V pi[util::TensorVal<DESCRIPTOR>::n] { };
    MomentaF().computeStress(cell, pi);
    V piNeqNormSqr { };
    MomentaF().computePiNeqNormSqr(cell, piNeqNormSqr);
    const V rho = MomentaF().computeRho(cell);
    const V omega = parameters.template get<descriptors::OMEGA>();
    const V smagorinsky = parameters.template get<collision::LES::SMAGORINSKY>();
    V piNeqNorm = util::sqrt(2*piNeqNormSqr);
    V S[util::TensorVal<DESCRIPTOR>::n];
    V cs2 = V{1} / descriptors::invCs2<V,DESCRIPTOR>();
    V tauMol = V{1} / omega;
    //Strain Tensor
    V Phi = (-0.5*(-rho*tauMol*cs2+util::sqrt(rho*rho*tauMol*tauMol*cs2*cs2+V{2}*(smagorinsky*smagorinsky)*rho*piNeqNorm))/(smagorinsky*smagorinsky*rho*piNeqNorm));
    for (int n=0; n < util::TensorVal<DESCRIPTOR>::n; ++n) {
      S[n] = Phi*pi[n];
    }
    //Strain Tensor Norm
    V SNormSqr = S[0]*S[0] + 2.0*S[1]*S[1] + S[2]*S[2];
    if constexpr (util::TensorVal<DESCRIPTOR>::n == 6) {
      SNormSqr += S[2]*S[2] + S[3]*S[3] + 2.0*S[4]*S[4] + S[5]*S[5];
    }
    V SNorm = util::sqrt(2*SNormSqr);
    V tauTurb = smagorinsky*smagorinsky*SNorm/cs2;
    V tauEff = tauMol + tauTurb;
    return V{1} / tauEff;
  }

References olb::descriptors::invCs2(), and 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