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

#include <collisionLES.h>

Collaboration diagram for olb::collision::detail::ConSmagorinskyEffectiveOmega< 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::ConSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >

Definition at line 292 of file collisionLES.h.

Member Typedef Documentation

CollisionO

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

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

Definition at line 294 of file collisionLES.h.

MomentaF

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

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

Definition at line 293 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::ConSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::apply ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 351 of file collisionLES.h.

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

References olb::collision::detail::ConSmagorinskyEffectiveOmega< 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::ConSmagorinskyEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::computeEffectiveOmega ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 297 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 H[util::TensorVal<DESCRIPTOR >::n];
    V conSmagoR[DESCRIPTOR::q];
    V S[util::TensorVal<DESCRIPTOR >::n];
    V tauMol = V{1} / omega;
    V cs2 = V{1} / descriptors::invCs2<V,DESCRIPTOR>();
    V piNeqNorm = util::sqrt(2*piNeqNormSqr);
    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 rate 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 preFactor = smagorinsky*smagorinsky
                * descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
                * 2 * util::sqrt(2);
    //consistent Samagorinsky additional R term
    //for (int q=0; q < DESCRIPTOR::q; ++q) {
    {
      unsigned q = 0;
      V t = descriptors::t<V,DESCRIPTOR>(q);
      //Hermite-Polynom H = c*c-cs^2*kronDelta
      H[0] = descriptors::c<DESCRIPTOR>(q,0)*descriptors::c<DESCRIPTOR>(q,0)-cs2;
      H[1] = descriptors::c<DESCRIPTOR>(q,0)*descriptors::c<DESCRIPTOR>(q,1);
      H[2] = descriptors::c<DESCRIPTOR>(q,1)*descriptors::c<DESCRIPTOR>(q,1)-cs2;
      if constexpr (util::TensorVal<DESCRIPTOR>::n == 6) {
        H[2] = descriptors::c<DESCRIPTOR>(q,0)*descriptors::c<DESCRIPTOR>(q,2);
        H[3] = descriptors::c<DESCRIPTOR>(q,1)*descriptors::c<DESCRIPTOR>(q,1)-cs2;
        H[4] = descriptors::c<DESCRIPTOR>(q,1)*descriptors::c<DESCRIPTOR>(q,2);
        H[5] = descriptors::c<DESCRIPTOR>(q,2)*descriptors::c<DESCRIPTOR>(q,2)-cs2;
      }
      //contraction or scalar product H*S
      V contractHS = H[0]*S[0] + 2.0*H[1]*S[1] + H[2]*S[2];
      if constexpr (util::TensorVal<DESCRIPTOR>::n == 6) {
        contractHS += H[2]*S[2] + H[3]*S[3] + 2.0*H[4]*S[4] + H[5]*S[5];
      }
      //additional term
      conSmagoR[q] = t*preFactor*SNorm*contractHS;
    }
    return conSmagoR[0];
  }

References olb::descriptors::c(), olb::descriptors::invCs2(), olb::util::sqrt(), and olb::descriptors::t().

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