AD coupling with Boussinesq bouancy for Smagorinsky-LES. More...

#include <navierStokesAdvectionDiffusionCoupling.h>

Collaboration diagram for olb::SmagorinskyBoussinesqCoupling:

Classes
struct	FORCE_PREFACTOR

struct	OMEGA_ADE

struct	OMEGA_NSE

struct	PR_TURB

struct	SMAGORINSKY_PREFACTOR

struct	T0

Public Types
using	parameters = meta::list<FORCE_PREFACTOR,SMAGORINSKY_PREFACTOR,PR_TURB,T0,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

AD coupling with Boussinesq bouancy for Smagorinsky-LES.

Definition at line 143 of file navierStokesAdvectionDiffusionCoupling.h.

Member Typedef Documentation

◆ parameters

using olb::SmagorinskyBoussinesqCoupling::parameters = meta::list<FORCE_PREFACTOR,SMAGORINSKY_PREFACTOR,PR_TURB,T0,OMEGA_NSE,OMEGA_ADE>

Definition at line 153 of file navierStokesAdvectionDiffusionCoupling.h.

Member Function Documentation

◆ apply()

template<typename CELLS , typename PARAMETERS >

void olb::SmagorinskyBoussinesqCoupling::apply	(	CELLS &	cells,
		PARAMETERS &	parameters )

inline

Molecular realaxation time

Turbulent realaxation time

Effective realaxation time

Definition at line 156 of file navierStokesAdvectionDiffusionCoupling.h.

  {
    using V = typename CELLS::template value_t<names::NavierStokes>::value_t;
    using DESCRIPTOR = typename CELLS::template value_t<names::NavierStokes>::descriptor_t;
    using DESCRIPTOR_ADE = typename CELLS::template value_t<names::Temperature>::descriptor_t;
 
    auto& cellNSE = cells.template get<names::NavierStokes>();
    auto& cellADE = cells.template get<names::Temperature>();
 
    // computation of the bousinessq force
    auto force = cellNSE.template getFieldPointer<descriptors::FORCE>();
    auto forcePrefactor = parameters.template get<FORCE_PREFACTOR>();
    V temperatureDifference = cellADE.computeRho() - parameters.template get<T0>();
 
    for (unsigned iD=0; iD < DESCRIPTOR::d; ++iD) {
      force[iD] = forcePrefactor[iD] * temperatureDifference;
    }
 
    // Velocity coupling
    auto u = cellADE.template getField<descriptors::VELOCITY>();
    // tau coupling
    auto tauNS = cellNSE.template getFieldPointer<descriptors::TAU_EFF>();
    auto tauAD = cellADE.template getFieldPointer<descriptors::TAU_EFF>();
 
    V rho, pi[util::TensorVal<DESCRIPTOR>::n] { };
    cellNSE.computeAllMomenta(rho, u.data(), pi);
    cellADE.template setField<descriptors::VELOCITY>(u);
    V 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];
    }
    V PiNeqNorm    = util::sqrt(PiNeqNormSqr);
    V tau_mol_NS = V{1} / parameters.template get<OMEGA_NSE>();
    V tau_mol_AD = V{1} / parameters.template get<OMEGA_ADE>();
    V smagoPrefactor = parameters.template get<SMAGORINSKY_PREFACTOR>();
    V tau_turb_NS = V{0.5}*(util::sqrt(tau_mol_NS*tau_mol_NS + smagoPrefactor/rho*PiNeqNorm) - tau_mol_NS);
    tauNS[0] = tau_mol_NS+tau_turb_NS;
 
    V prTurb = parameters.template get<PR_TURB>();
    V tauTurbADPrefactor = descriptors::invCs2<V,DESCRIPTOR_ADE>() / descriptors::invCs2<V,DESCRIPTOR>() / prTurb;
    V tau_turb_AD = tau_turb_NS * tauTurbADPrefactor;
    tauAD[0] = tau_mol_AD+tau_turb_AD;
  }