LES-ADE coupling for multiple reactions. More...

#include <navierStokesAdvectionDiffusionCoupling.h>

Collaboration diagram for olb::LESReactionCoupling< T, numComp >:

Classes
struct	LATTICE_REACTION_COEFF

struct	OMEGA_NSE

struct	OMEGAS_ADE

struct	REACTION_ORDER

struct	SCHMIDT

struct	SMAGORINSKY_PREFACTOR

struct	STOCH_COEFF

Public Types
using	parameters = meta::list<LATTICE_REACTION_COEFF, STOCH_COEFF, REACTION_ORDER, SMAGORINSKY_PREFACTOR, SCHMIDT, OMEGA_NSE, OMEGAS_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, int numComp>
struct olb::LESReactionCoupling< T, numComp >

LES-ADE coupling for multiple reactions.

Definition at line 233 of file navierStokesAdvectionDiffusionCoupling.h.

Member Typedef Documentation

◆ parameters

template<typename T , int numComp>

using olb::LESReactionCoupling< T, numComp >::parameters = meta::list<LATTICE_REACTION_COEFF, STOCH_COEFF, REACTION_ORDER, SMAGORINSKY_PREFACTOR, SCHMIDT, OMEGA_NSE, OMEGAS_ADE>

Definition at line 244 of file navierStokesAdvectionDiffusionCoupling.h.

Member Function Documentation

◆ apply()

template<typename T , int numComp>

template<typename CELLS , typename PARAMETERS >

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

inline

Molecular realaxation time

Turbulent realaxation time

Definition at line 247 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);
    cells.template get<names::Concentration1>().template setField<descriptors::VELOCITY>(u);
    cells.template get<names::Concentration2>().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>();
    auto omegasAD = parameters.template get<OMEGAS_ADE>();
    T tau_mol_AD0 = T{1} / omegasAD[0];
    T tau_mol_AD1 = T{1} / omegasAD[1];
    T tau_mol_AD2 = T{1} / omegasAD[2];
    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
    auto Sc = parameters.template get<SCHMIDT>();
    T tauTurbADPrefactor0 = descriptors::invCs2<T,DESCRIPTOR_ADE>() / descriptors::invCs2<T,DESCRIPTOR>() / Sc[0];
    T tauTurbADPrefactor1 = descriptors::invCs2<T,DESCRIPTOR_ADE>() / descriptors::invCs2<T,DESCRIPTOR>() / Sc[1];
    T tauTurbADPrefactor2 = descriptors::invCs2<T,DESCRIPTOR_ADE>() / descriptors::invCs2<T,DESCRIPTOR>() / Sc[2];
    T tau_turb_AD0 = tau_turb_NS * tauTurbADPrefactor0;
    T tau_turb_AD1 = tau_turb_NS * tauTurbADPrefactor1;
    T tau_turb_AD2 = tau_turb_NS * tauTurbADPrefactor2;
    cells.template get<names::Concentration0>().template setField<descriptors::OMEGA>(T{1} / (tau_mol_AD0 + tau_turb_AD0));
    cells.template get<names::Concentration1>().template setField<descriptors::OMEGA>(T{1} / (tau_mol_AD1 + tau_turb_AD1));
    cells.template get<names::Concentration2>().template setField<descriptors::OMEGA>(T{1} / (tau_mol_AD2 + tau_turb_AD2));
 
    auto stochCoeff = parameters.template get<STOCH_COEFF>();
    auto reactionOrder = parameters.template get<REACTION_ORDER>();
    T source = parameters.template get<LATTICE_REACTION_COEFF>();
    {
      T conc = cells.template get<names::Concentration0>().computeRho();
      source *= util::pow(conc, reactionOrder[0]);
    }
    {
      T conc = cells.template get<names::Concentration1>().computeRho();
      source *= util::pow(conc, reactionOrder[1]);
    }
    {
      T conc = cells.template get<names::Concentration2>().computeRho();
      source *= util::pow(conc, reactionOrder[2]);
    }
    {
      cells.template get<names::Concentration0>().template setField<descriptors::SOURCE>(stochCoeff[0]*source);
      cells.template get<names::Concentration1>().template setField<descriptors::SOURCE>(stochCoeff[1]*source);
      cells.template get<names::Concentration2>().template setField<descriptors::SOURCE>(stochCoeff[2]*source);
    }
  }