Skip to content

Re: Smagorinsky/Turbulence – Modification Code

#2729
Marc
Participant

Hey Pedro,
if you want to use a Dynamics pointer you have to adapt the arguments in the function prepareLattice.
The easiest way to run the model is to change the line:
BGKdynamics<T, DESCRIPTOR> bulkDynamics( converter.getOmega(), instances::getBulkMomenta<T, DESCRIPTOR>() );
to
SmagorinskyBGKdynamics<T, DESCRIPTOR>bulkDynamics( converter.getOmega(), instances::getBulkMomenta<T, DESCRIPTOR>(), 0.15);
Good luck with your turbulence Simulation!
Best,
Marc