Scalable Fully Implicit Runge-Kutta Solvers for Time-Dependent Rayleigh-Bénard Convection
Abstract
We consider the efficient numerical solution of time-dependent Rayleigh-Bénard convection governed by the incompressible Boussinesq equations. The coupled velocity-pressure-temperature system presents significant computational challenges arising from incompressibility, buoyancy coupling, convection, and the simultaneous resolution of multiple physical time scales. We present a systematic study of fully implicit Runge-Kutta methods for time-dependent Rayleigh-Bénard convection, together with scalable finite element solvers for the resulting stage-coupled systems. We employ high-order, L-stable Radau IIA methods for temporal integration and investigate Taylor-Hood and Alfeld-Sorokina finite elements in two dimensions, along with Taylor-Hood and enriched macroelement discretizations in three dimensions. The macroelement approaches offer strong local inf-sup stability and patch structures well suited to Vanka-type relaxation. In two dimensions, monolithic multigrid achieves nearly stage-independent outer Krylov iteration counts and maintains robust convergence under mesh refinement. On the same mesh with incremental pressure correction, selected Radau configurations reduce elapsed solve time by 18-32% and yield smaller sampled transient differences relative to a numerical reference solution. Three-dimensional experiments demonstrate favorable parallel scalability with Taylor-Hood elements and substantially improved discrete mass conservation with the enriched macroelement discretization. These results identify regimes in which coupled multistage fully implicit integration is computationally competitive when equipped with appropriate multigrid preconditioners.
Related
- Space-Time Parallel Direct Numerical Simulation of Rayleigh-Bénard Convection
- A second-order linear decoupled scheme with a relaxed zero-energy variable for Rayleigh-Bénard convection
- GPU-Accelerated Nested Schur-Complement Solvers for Atmospheric Flow Simulations Using Horizontally-Explicit, Vertically-Implicit Schemes and Discontinuous Galerkin Methods
- Fast and scalable multi-GPU sliced FFT solvers for asymptotically compatible discretizations of 3D nonlocal diffusion problems
- Energy stability of explicit--implicit--null time discretizations for nonlinear degenerate parabolic equations with state-dependent coefficients
BibTeX
@misc{ali2026scalable,
title = {{Scalable Fully Implicit Runge-Kutta Solvers for Time-Dependent Rayleigh-Bénard Convection}},
author = {Ahsan Ali and Robert C. Kirby},
howpublished = {arXiv},
year = {2026},
month = nov,
doi = {10.48550/arxiv.2611.03361},
eprint = {2611.03361},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2611.03361},
}
RIS
TY - UNPB
TI - Scalable Fully Implicit Runge-Kutta Solvers for Time-Dependent Rayleigh-Bénard Convection
AU - Ahsan Ali
AU - Robert C. Kirby
PY - 2026
DA - 2026/11/02
PB - arXiv
DO - 10.48550/arxiv.2611.03361
UR - https://arxiv.org/abs/2611.03361
AB - We consider the efficient numerical solution of time-dependent Rayleigh-Bénard convection governed by the incompressible Boussinesq equations. The coupled velocity-pressure-temperature system presents significant computational challenges arising from incompressibility, buoyancy coupling, convection, and the simultaneous resolution of multiple physical time scales. We present a systematic study of fully implicit Runge-Kutta methods for time-dependent Rayleigh-Bénard convection, together with scalable finite element solvers for the resulting stage-coupled systems. We employ high-order, L-stable Radau IIA methods for temporal integration and investigate Taylor-Hood and Alfeld-Sorokina finite elements in two dimensions, along with Taylor-Hood and enriched macroelement discretizations in three dimensions. The macroelement approaches offer strong local inf-sup stability and patch structures well suited to Vanka-type relaxation. In two dimensions, monolithic multigrid achieves nearly stage-independent outer Krylov iteration counts and maintains robust convergence under mesh refinement. On the same mesh with incremental pressure correction, selected Radau configurations reduce elapsed solve time by 18-32% and yield smaller sampled transient differences relative to a numerical reference solution. Three-dimensional experiments demonstrate favorable parallel scalability with Taylor-Hood elements and substantially improved discrete mass conservation with the enriched macroelement discretization. These results identify regimes in which coupled multistage fully implicit integration is computationally competitive when equipped with appropriate multigrid preconditioners.
ER -