Reply To: Smagorinsky method in AD lattice

[FBukreev]

Hello,

you can couple turbulent NS simulation with AD using special coupler and AD dynamics.
For AD lattice you can use the following descriptor and dynamics:

using ADDESCRIPTOR = D3Q7;
using MOMENTA = momenta::AdvectionDiffusionBulkTuple;
using ADBulkDynamics = dynamics::Tuple,AdvectionDiffusionExternalVelocityCollision>;

the coupler is declared as following:

SuperLatticeCoupling coupling(
LESADECoupling{},
names::NavierStokes{}, sLatticeNS,
names::Concentration0{}, sLatticeAD);
coupling.setParameter::SMAGORINSKY_PREFACTOR>(0.15);
coupling.setParameter::SCHMIDT>(0.05); // turbulent Schmidt number for stabilization
coupling.setParameter::OMEGA_NSE>(converter.getLatticeRelaxationFrequency());
coupling.setParameter::OMEGA_ADE>(converter.getLatticeRelaxationFrequencyFromDiffusivity(diffusivity));

For the NS lattice you can apply standard SmagorinskyBGKdynamics.

Best wishes
Fedor