I’m implementing a D3Q19 model, with bounce back at walls and wrap around in the direction of flow using body force forcing.
My problem is, that after running a couple(about 10000) simulations, i have found that the highest reynold’s number at which the velocity profile converged is very small- about 300(cylindrical pipe flow). At higher reynolds no., and a lot of lower ones as well, the parabola collapsed into random fluctuations.
What (if any) among these: tao,body force, lattice dimensions, reynolds no., max velocity or perhaps some combination, is the parameter which can definitively predict numerical collapse?(from my data i’m certain it is not Re, as i was led to believe)
Is any research available on the subject?