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

#include <collisionLES.h>

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

Definition at line 92 of file collisionLES.h.

Member Typedef Documentation

CollisionO

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

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

Definition at line 94 of file collisionLES.h.

MomentaF

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

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

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

inline

Definition at line 123 of file collisionLES.h.

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

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

inline

Molecular realaxation time

Turbulent realaxation time

Effective realaxation time

Definition at line 97 of file collisionLES.h.

                                                                           {
    V piNeqNormSqr { };
    MomentaF().computePiNeqNormSqr(cell, piNeqNormSqr);
    V u[DESCRIPTOR::d];
    MomentaF().computeU(cell, u);
    //const V rho = MomentaF().computeRho(cell);
    const V rho = cell.template getField<descriptors::RHO>();
    const auto nablaRho = cell.template getField<descriptors::NABLARHO>();
    V nablaRhoSqr = util::normSqr<V,DESCRIPTOR::d>(nablaRho);
    V nablaRhoU = nablaRho[0]*u[0]+nablaRho[1]*u[1];
    const V omega = parameters.template get<descriptors::OMEGA>();
    const V smagorinsky = parameters.template get<collision::LES::SMAGORINSKY>();
    V piNeqNorm = util::sqrt(piNeqNormSqr);
    V preFactor = smagorinsky*smagorinsky
                * descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
                * 2 * util::sqrt(2);
    V tauMol = V{1} / omega;
    V tauTurb = V{0.5} * (util::sqrt(tauMol*tauMol + preFactor / rho * piNeqNorm) - tauMol);
    V tauEff = tauMol + tauTurb;
    return  V{1} / tauEff;
  }

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