Efficient monolithic projection method with staggered time discretization for natural convection problems

For a more efficient algorithm, we introduce staggered time discretization to improve the previous method (Pan et al., 2017), fully decoupled monolithic projection method with one Poisson equation (FDMPM-1P), to solve time-dependent natural convection problems. The momentum and energy equations are discretized using the Crank–Nicolson scheme at the staggered time grids, in which temperature and pressure fields are evaluated at half-integer time levels (n+ [Formula presented] ), while the velocity fields are evaluated at integer time levels (n+1). Numerical simulations of two-dimensional (2D) Rayleigh–Bénard convection show that the proposed method is more computationally efficient and stable than FDMPM-1P, while preserving the second-order spatial and temporal accuracy. Further, the proposed method provides accurate predictions of 2D Rayleigh–Bénard convection under different thermal boundary conditions for a Rayleigh number up to 10¹⁰, three-dimensional turbulent Rayleigh–Bénard convection in the range of 1×10⁵–2×10⁷ in horizontal periodic domain, and three-dimensional turbulent Rayleigh–Bénard convection in the range of 1×10⁶–1×10⁷ in bounded domain. Finally, we theoretically confirmed that the proposed method is second-order in time and it is more stable than FDMPM-1P for small Ra.