LES-ADE coupling with Schmidt number stabilization. More...

#include <navierStokesAdvectionDiffusionCoupling.h>

Collaboration diagram for olb::LESADECoupling< T >:

Classes
struct	OMEGA_ADE

struct	OMEGA_NSE

struct	SCHMIDT

struct	SMAGORINSKY_PREFACTOR

Public Types
using	parameters = meta::list<SMAGORINSKY_PREFACTOR, SCHMIDT, OMEGA_NSE, OMEGA_ADE>

Public Member Functions
template<typename CELLS , typename PARAMETERS >
void	apply (CELLS &cells, PARAMETERS &parameters) any_platform

Static Public Attributes
static constexpr OperatorScope	scope = OperatorScope::PerCellWithParameters

Detailed Description

template<typename T>
struct olb::LESADECoupling< T >

LES-ADE coupling with Schmidt number stabilization.

Definition at line 311 of file navierStokesAdvectionDiffusionCoupling.h.

Member Typedef Documentation

◆ parameters

template<typename T >

using olb::LESADECoupling< T >::parameters = meta::list<SMAGORINSKY_PREFACTOR, SCHMIDT, OMEGA_NSE, OMEGA_ADE>

Definition at line 319 of file navierStokesAdvectionDiffusionCoupling.h.

Member Function Documentation

◆ apply()

template<typename T >

template<typename CELLS , typename PARAMETERS >

void olb::LESADECoupling< T >::apply	(	CELLS &	cells,
		PARAMETERS &	parameters )

inline

Molecular realaxation time

Turbulent realaxation time

Definition at line 322 of file navierStokesAdvectionDiffusionCoupling.h.

   {
    using DESCRIPTOR = typename CELLS::template value_t<names::NavierStokes>::descriptor_t;
    using DESCRIPTOR_ADE = typename CELLS::template value_t<names::Concentration0>::descriptor_t;
    // Velocity coupling
    auto u = cells.template get<names::Concentration0>().template getField<descriptors::VELOCITY>();
    T rho, pi[util::TensorVal<DESCRIPTOR>::n] { };
    cells.template get<names::NavierStokes>().computeAllMomenta(rho, u.data(), pi);
    cells.template get<names::Concentration0>().template setField<descriptors::VELOCITY>(u);
    // Stress tensor
    T PiNeqNormSqr = pi[0]*pi[0] + 2.0*pi[1]*pi[1] + pi[2]*pi[2];
    if constexpr (util::TensorVal<DESCRIPTOR>::n == 6) {
      PiNeqNormSqr += pi[2]*pi[2] + pi[3]*pi[3] + 2*pi[4]*pi[4] +pi[5]*pi[5];
    }
    T PiNeqNorm    = util::sqrt(PiNeqNormSqr);
    T tau_mol_NS = T{1} / parameters.template get<OMEGA_NSE>();
    T tau_mol_AD = T{1} / parameters.template get<OMEGA_ADE>();
    T smagoPrefactor = parameters.template get<SMAGORINSKY_PREFACTOR>();
    T tau_turb_NS = T{0.5}*(util::sqrt(tau_mol_NS*tau_mol_NS + smagoPrefactor/rho*PiNeqNorm) - tau_mol_NS);
    // Schmidt number stabilization
    T tauTurbADPrefactor = descriptors::invCs2<T,DESCRIPTOR_ADE>() / descriptors::invCs2<T,DESCRIPTOR>() / parameters.template get<SCHMIDT>();
    T tau_turb_AD = tau_turb_NS * tauTurbADPrefactor;
    cells.template get<names::Concentration0>().template setField<descriptors::OMEGA>(T{1} / (tau_mol_AD + tau_turb_AD));
  }