We develop a set of numerical approximations to discretize the non-isothermal incompressible hydrodynamic model for Rayleigh-Bénard convection (RBC), which ensures the negative energy dissipation rate with respect to adiabatic boundary conditions. Using the Crank-Nicolson (CN) method, the second-order backward difference method combined with the pressure-correction method, we propose two kinds of decoupled, linear, both second-order in time energy-dissipation-rate-preserving semi-discrete projection numerical algorithms. Meanwhile, the second-order fully discrete numerical algorithms are obtained by the use of a finite difference method on staggered grids in space. These numerical approximations are proved to preserve the property of the energy dissipation rate at the fully discrete level. Moreover, the numerical approximations are also unconditionally energy stable. The fast Fourier algorithm is applied to numerical implementation to enhance experimental efficiency. Convergence rate tests are conducted to verify the accuracy of the algorithms. Our simulations showcase that both hydrodynamic and thermal effects in resolving RBC within the non-isothermal hydrodynamic model of incompressible viscous fluid flow. In general, our structure-preserving projection schemes and implementation methods accurately and efficiently simulate RBC in nature.
| [1] |
Bodenschatz E, Pesch W, Ahlers G. Recent developments in Rayleigh-Bénard convection. Annu. Rev. Fluid Mech., 2000, 32(1): 709-778
|
| [2] |
Cheng Q, Liu C, Shen J. A new Lagrange multiplier approach for gradient flows. Comput. Methods Appl. Mech. Eng., 2020, 367 Article ID: 113070
|
| [3] |
Favier B, Purseed J, Duchemin L. Rayleigh-Bénard convection with a melting boundary. J. Fluid Mech., 2019, 858: 437-473
|
| [4] |
Gong YZ, Zhao J, Wang Q. Second order fully discrete energy stable methods on staggered grids for hydrodynamic phase field models of binary viscous fluids. SISC, 2018, 40(2): B528-B553
|
| [5] |
Gong YZ, Zhao J, Wang Q. Arbitrarily high-order linear energy stable schemes for gradient flow models. J. Comput. Phys., 2020, 419 Article ID: 109610
|
| [6] |
Gong YZ, Zhao J, Wang Q. Arbitrarily high-order unconditionally energy stable SAV schemes for gradient flow models. Comput. Phys. Commun., 2020, 249 Article ID: 107033
|
| [7] |
Guermond JL, Minev P, Shen J. An overview of projection methods for incompressible flows. Comput. Methods Appl. Mech. Eng., 2006, 195(44): 6011-6045
|
| [8] |
Julien K, Knobloch E, Rubio AM, Vasil GM. Heat transport in Low-Rossby-Number Rayleigh-Bénard convection. Phys. Rev. Lett., 2012, 109(25): Article ID: 254503
|
| [9] |
Li J, Zhao J, Wang Q. Energy and entropy preserving numerical approximations of thermodynamically consistent crystal growth models. J. Comput. Phys., 2019, 328: 202-220
|
| [10] |
Shen, J., Tang, T., Wang, L.: Spectral methods: algorithms, analysis and applications. In: Springer Series in Computational Mathematics, Springer, Berlin-Heidelberg, 41 (2011)
|
| [11] |
Shen J, Xu J, Yang J. The scalar auxiliary variable (SAV) approach for gradient flows. J. Comput. Phys., 2018, 353: 407-416
|
| [12] |
Shen J, Xu J, Yang J. A new class of efficient and robust energy stable schemes for gradient flows. SIAM Rev., 2019, 61(3): 474-506
|
| [13] |
Shen J, Zheng N. Efficient and unconditional energy stable schemes for the micropolar Navier-Stokes equations. CSIAM Trans. Appl. Math., 2022, 3(1): 57-81
|
| [14] |
Sun, S.W., Jing, X.B., Wang, Q.: Error estimates of energy stable numerical schemes for Allen-Cahn equations with nonlocal constraints. J. Sci. Comput. 79, 593–623 (2019)
|
| [15] |
Sun, S.W., Li, J., Zhao, J., Wang, Q.: Structure-preserving numerical approximations to a non-isothermal hydrodynamic model of binary fluid flows. J. Sci. Comput. 83, 3 (2020)
|
| [16] |
Sun SW, Wang Q. Entropy-Production-Rate-Preserving algorithms for a hydrodynamical model of binary fluids. J. Sci. Comput., 2024, 101: 53
|
| [17] |
Watanabe T. Flow pattern and heat transfer rate in Rayleigh-Bénard convection. Phys. Fluids, 2004, 16(4): 972-978
|
| [18] |
Wen B, Goluskin D, Doering CR. Steady Rayleigh-Bénard convection between no-slip boundaries. J. Fluid Mech., 2022, 933: R4
|
| [19] |
Xie X, Zhao W, Lin P. A simple and efficient curved boundary scheme of the lattice Boltzmann method for Robin boundary conditions of convection-diffusion equations. Appl. Math. Lett., 2021, 122 Article ID: 107536
|
| [20] |
Yue P, Feng J, Liu C, Shen J. A diffuse-interface method for simulating two-phase flows of complex fluids. J. Fluid Mech., 2004, 515: 293-317
|
Funding
National Natural Science Foundation of China(12101387)
International scientific and technological cooperation projects of Henan Province, China
Henan International Joint Laboratory of Optical Information Transmission and Application
RIGHTS & PERMISSIONS
Shanghai University
PDF
