Double MRT Thermal Lattice Boltzmann Method for Simulating Natural Convection of Low Prandtl Number Fluids

Zheng Li,Mo Yang,Yuwen Zhang
DOI: https://doi.org/10.1108/HFF-04-2015-0135
2016-01-19
Abstract:The purposes of this paper are testing an efficiency algorithm based on LBM and using it to analyze two-dimensional natural convection with low Prandtl number. Steady state or oscillatory results are obtained using double multiple-relaxation-time thermal lattice Boltzmann method. The velocity and temperature fields are solved using D2Q9 and D2Q5 models, respectively. With different Rayleigh number, the tested natural convection can either achieve to steady state or oscillatory. With fixed Rayleigh number, lower Prandtl number leads to a weaker convection effect, longer oscillation period and higher oscillation amplitude for the cases reaching oscillatory solutions. At fixed Prandtl number, higher Rayleigh number leads to a more notable convection effect and longer oscillation period. Double multiple-relaxation-time thermal lattice Boltzmann method is applied to simulate the low Prandtl number fluid natural convection. Rayleigh number and Prandtl number effects are also investigated when the natural convection results oscillate.
Computational Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to use the Double Multiple - Relaxation - Time Thermal Lattice Boltzmann Method (DMRT - LBM) to simulate the natural convection phenomenon of fluids with low Prandtl number (\( \text{Pr} \)). Specifically, the research objectives include: 1. **Verify the effectiveness of the algorithm**: Test the efficient algorithm based on LBM and use it to analyze the natural convection of two - dimensional low - Prandtl - number fluids. 2. **Numerical simulation and analysis**: Through the DMRT - LBM model, solve the steady - state or oscillatory results under different Reynolds numbers (\( \text{Ra} \)) and Prandtl numbers. 3. **Explore the physical mechanism**: Study the influence of a lower Prandtl number on natural convection under a fixed Reynolds number; and the influence of a higher Reynolds number on natural convection under a fixed Prandtl number. ### Main research content - **Model selection**: Use the D2Q9 model to solve the velocity field and the D2Q5 model to solve the temperature field. - **Boundary conditions**: Set the left boundary of the square cavity as a constant high temperature, the right boundary as a constant low temperature, and the upper and lower boundaries as adiabatic and no - slip boundary conditions. - **Parameter range**: Four different combinations of Reynolds numbers and Prandtl numbers were studied, covering Reynolds numbers from \( 10^4 \) to \( 10^5 \) and Prandtl numbers from 0.005 to 0.01. ### Research findings - **Steady - state and oscillation**: For some parameter combinations, natural convection can reach a steady - state or oscillatory state. - **Influence of Prandtl number**: In the case of a fixed Reynolds number, a lower Prandtl number leads to a weaker convection effect, a longer oscillation period, and a higher oscillation amplitude. - **Influence of Reynolds number**: In the case of a fixed Prandtl number, a higher Reynolds number leads to a more significant convection effect and a longer oscillation period. ### Conclusion This study verifies the effectiveness of DMRT - LBM in simulating the natural convection of low - Prandtl - number fluids and reveals the specific influences of Prandtl number and Reynolds number on the natural convection behavior. These findings are helpful for understanding the nonlinear behavior in complex fluid dynamics and provide theoretical support for industrial applications.