#include <collisionLES.h>

Collaboration diagram for olb::collision::detail::LocalVanDriestSmagorinskyEffectiveOmega< 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

Static Public Attributes
static constexpr bool	is_vectorizable = false

Detailed Description

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

Definition at line 165 of file collisionLES.h.

Member Typedef Documentation

◆ CollisionO

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

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

Definition at line 167 of file collisionLES.h.

◆ MomentaF

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

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

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

inline

Definition at line 199 of file collisionLES.h.

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

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

inline

Molecular realaxation time

Turbulent realaxation time

Effective realaxation time

Definition at line 172 of file collisionLES.h.

                                                                           {
    V piNeqNormSqr { };
    MomentaF().computePiNeqNormSqr(cell, piNeqNormSqr);
    const V rho = MomentaF().computeRho(cell);
    const V omega = parameters.template get<descriptors::OMEGA>();
    V tauEff {};
    V nuT = cell.template getField<descriptors::VISCOSITY>();
    if ( nuT != V(0) ) {
      tauEff = nuT*descriptors::invCs2<V,DESCRIPTOR>() + V(0.5);
    } else {
      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);
      tauEff = tauMol + tauTurb;
    }
    cell.template setField<descriptors::OMEGA>(V{1} / tauEff);
    return  V{1} / tauEff;
  }