Stability of 3D thermal LBM model
- This topic has 2 replies, 2 voices, and was last updated 2 years, 10 months ago by .
-
Posts
-
Hello,
I am working on thermal LBM and during the solution of the benchmark problem of natural convection in a differentially heated cavity I found out that the 3D model is not as stable as the 2D model. I use the D3Q19 scheme with the two-relaxation rate approximation for flow whereas the advection/diffusion equation is solved by the conventional finite difference technique. As expected, the instabilities are observed with the Rayleigh number higher than 10^9. However, the 2D model with the D2Q9 scheme works relatively fine with Ra=10^9 and Ra=10^10 whereas the 3D model is unstable with the same values of relaxation rates and grid size along two coordinates. My first thought was to extend the number of speeds to D3Q27. Unfortunately, it did not work. Did someone meet this problem? Or can someone suggest to me the direction of this problem solution?Hi Alex,
3D simulations being easier to be instable sounds resonable. You have more degrees of freedom, more calculations that could potentially intodruce numerical errors. Also, do not forget that turbulence in 2D is a completely different phenomenon than in 3D!
Above Ra>=1e7 these simulations transition to turbulent natural convection and eddies will arise. Are you using a turbulence model for your simulations or are you wanting to resolve all structures? If no turbulence model is applied, check your resolution. It should be high enough to resolve the boundary layer that develops close to the heated and cooled walls. You can find some more info in e.g. https://doi.org/10.1016/j.camwa.2018.08.018.
Best,
Max
- You must be logged in to reply to this topic.