Propagating fronts of convection rolls in Rayleigh-Bénard convection
Abstract
We investigate the propagation of counter-rotating convection rolls in Rayleigh-Bénard convection initiated locally in a quiescent fluid layer under supercritical conditions. The velocity of the front separating quiescent fluid from the forming convection rolls, and the wavenumber of the convection rolls remaining behind the front, are explored. We numerically investigate fronts of forming convection rolls over five orders of magnitude of the reduced Rayleigh number, ε, in 2D and 3D domains, for a broad range of boundary conditions, and for different front initiation approaches. In all cases, the front velocity increases as ε1/2 with increasing ε for ε≲ 1 in agreement with predictions using the amplitude equation. The amplitude equation description of the front velocity remains accurate for ε≲ 10 except when the Prandtl number is large which yields a velocity that is faster than predicted for a fluid layer far from threshold. The wavenumber of the convection rolls increases linearly with ε in agreement with the wavenumber that maximizes the growth rate of perturbations in the linear regime. Farther from onset, the wavenumber growth transitions to a reduced scaling of ε1/4 in agreement with predictions using the Swift-Hohenberg equation in the large ε limit. The scalings describing the wavenumber variation with ε are independent of the domain geometry, boundary conditions, and front initiation method. However, the front-selected wavenumber at criticality does not equal the critical wavenumber of the bulk instability, in general, and depends significantly upon these details. We compare our results with experimental measurements where possible.
Versions
Related
- Space-Time Parallel Direct Numerical Simulation of Rayleigh-Bénard Convection
- A Lagrangian analysis of coherent structures in 2D annular Rayleigh-Bénard convection
- Does the Boussinesq approximation correctly describe statistically stationary states of Rayleigh-Bénard convection?
- A second-order linear decoupled scheme with a relaxed zero-energy variable for Rayleigh-Bénard convection
- A method for measuring and analyzing four-dimensional flow fields in viscous Rayleigh-Bénard convection
BibTeX
@misc{mukherjee2026propagating,
title = {{Propagating fronts of convection rolls in Rayleigh-Bénard convection}},
author = {Saikat Mukherjee and Mark Paul},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2608.21728},
eprint = {2608.21728},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2608.21728},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Propagating fronts of convection rolls in Rayleigh-Bénard convection
AU - Saikat Mukherjee
AU - Mark Paul
PY - 2026
DA - 2026/10/08
PB - arXiv
DO - 10.48550/arxiv.2608.21728
UR - https://arxiv.org/abs/2608.21728
AB - We investigate the propagation of counter-rotating convection rolls in Rayleigh-Bénard convection initiated locally in a quiescent fluid layer under supercritical conditions. The velocity of the front separating quiescent fluid from the forming convection rolls, and the wavenumber of the convection rolls remaining behind the front, are explored. We numerically investigate fronts of forming convection rolls over five orders of magnitude of the reduced Rayleigh number, ε, in 2D and 3D domains, for a broad range of boundary conditions, and for different front initiation approaches. In all cases, the front velocity increases as ε^1/2 with increasing ε for ε≲ 1 in agreement with predictions using the amplitude equation. The amplitude equation description of the front velocity remains accurate for ε≲ 10 except when the Prandtl number is large which yields a velocity that is faster than predicted for a fluid layer far from threshold. The wavenumber of the convection rolls increases linearly with ε in agreement with the wavenumber that maximizes the growth rate of perturbations in the linear regime. Farther from onset, the wavenumber growth transitions to a reduced scaling of ε^1/4 in agreement with predictions using the Swift-Hohenberg equation in the large ε limit. The scalings describing the wavenumber variation with ε are independent of the domain geometry, boundary conditions, and front initiation method. However, the front-selected wavenumber at criticality does not equal the critical wavenumber of the bulk instability, in general, and depends significantly upon these details. We compare our results with experimental measurements where possible.
ER -