Dear Community,
When checking SmagorinskyForcedBGKdynamics source code, I seem to have found the problem. In collisionLES.h file, struct SmagorinskyEffectiveOmega, the function computeEffectiveOmega to compute the turbulent realaxation time, is the preFactor term miscalculated and sqrt(2) redundant?
I refer to the “Application of a lattice Boltzmann method combined with a Smagorinsky turbulence model to spatially resolved heat flux inside a refrigerated vehicle “Section 2.2 of the paper.
https://doi.org/https://doi.org/10.1016/j.camwa.2018.08.018
I hope you can check this part, thank you!