I wish all of you are fine these days, I am running a 3D homogeneous case for a square cylinder, when I run the case at high Reynolds 2000 it keeps crashing, I tried
"SpectralVanishingViscosity" it is working but it is filtering my field, what about using a different time integration method
"BDFImplicitOrder2" instead of using
"IMEXOrder2", I think the two of them are 2nd order but using Implicit instead of Implicit explicit will not crash, but also I will keep CFL under 1 to be sure it is fine, even with this method i can go for higher CFL, my Reynolds is 2000 turbulent flow.