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

#include <collisionLES.h>

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

Definition at line 445 of file collisionLES.h.

Member Typedef Documentation

CollisionO

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

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

Definition at line 447 of file collisionLES.h.

MomentaF

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

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

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

inline

Definition at line 474 of file collisionLES.h.

                                                                          {
    V rho, u[DESCRIPTOR::d], fNeq[DESCRIPTOR::q] { };
    MomentaF().computeRhoU(cell, rho, u);
    lbm<DESCRIPTOR>::computeFneq(cell, fNeq, rho, u); // TODO: Use EQUILIBRIUM instead of 2nd order
    const V omega = parameters.template get<descriptors::OMEGA>();
    const V smagorinsky = parameters.template get<collision::LES::SMAGORINSKY>();
    const V preFactor = smagorinsky*smagorinsky
                      * 3*descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
                      * 2*util::sqrt(2);
    FieldD<V,DESCRIPTOR,typename COLLISION::OMEGA> omegaEff;
    const V tauMol = V{1} / omega;
    for (unsigned iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      V tauTurb = V{0.5}*(util::sqrt(tauMol*tauMol + preFactor/rho * util::fabs(fNeq[iPop])) - tauMol);
      omegaEff[iPop] = V{1} / (tauMol + tauTurb);
    }
    parameters.template set<typename COLLISION::OMEGA>(omegaEff);
    return CollisionO().apply(cell, parameters);
  }

References olb::lbm< DESCRIPTOR >::computeFneq(), olb::util::fabs(), olb::descriptors::invCs2(), and olb::util::sqrt().

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::KrauseEffectiveOmega< COLLISION, DESCRIPTOR, MOMENTA, EQUILIBRIUM >::computeEffectiveOmega ( CELL & cell, PARAMETERS & parameters )

inline

Definition at line 450 of file collisionLES.h.

                                                                           {
    V rho, u[DESCRIPTOR::d], fNeq[DESCRIPTOR::q] { };
    MomentaF().computeRhoU(cell, rho, u);
    lbm<DESCRIPTOR>::computeFneq(cell, fNeq, rho, u);
    const V omega = parameters.template get<descriptors::OMEGA>();
    const V smagorinsky = parameters.template get<collision::LES::SMAGORINSKY>();
    const V preFactor = smagorinsky*smagorinsky
                      * 3*descriptors::invCs2<V,DESCRIPTOR>()*descriptors::invCs2<V,DESCRIPTOR>()
                      * 2*util::sqrt(2);
    FieldD<V,DESCRIPTOR,typename COLLISION::OMEGA> omegaEff;
    const V tauMol = V{1} / omega;
    for (unsigned iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      V tauTurb = V{0.5}*(util::sqrt(tauMol*tauMol + preFactor/rho * util::fabs(fNeq[iPop])) - tauMol);
      omegaEff[iPop] = V{1} / (tauMol + tauTurb);
    }
    V avgOmegaEff = 0.;
    for (int iPop=0; iPop < DESCRIPTOR::q; ++iPop) {
      avgOmegaEff += omegaEff[iPop];
    }
    avgOmegaEff /= DESCRIPTOR::q;
    return avgOmegaEff;
  }

References olb::lbm< DESCRIPTOR >::computeFneq(), olb::util::fabs(), olb::descriptors::invCs2(), and olb::util::sqrt().

Here is the call graph for this function:

---

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

• src/dynamics/collisionLES.h