Quantum-Inspired Simulation of 2D Turbulent Rayleigh-B\'enard Convection

TL;DR

This paper demonstrates that Matrix Product State (MPS) methods can efficiently simulate 2D turbulent Rayleigh-Bénard convection at high Rayleigh numbers, showing favorable scaling and significant reduction in computational complexity.

Contribution

It introduces the application of MPS to buoyancy-driven flows with thermal coupling, revealing scalable simulation capabilities for high-Rayleigh-number turbulence.

Findings

MPS bond dimension grows without saturation up to Ra=10^{11} in a priori analysis.

Dynamical simulations at fixed bond dimension recover statistical observables with less complexity.

Achieved 1.8% error in Nusselt number at Ra=10^{10} with nearly 9-fold reduction in degrees of freedom.

Abstract

Turbulent thermal convection governs heat transport in systems ranging from stellar interiors to industrial heat exchangers. Two-dimensional Rayleigh-B\'enard convection serves as a paradigm for these flows, reproducing key features such as thin boundary layers, large-scale circulation, and sustained plume dynamics. While Matrix Product State (MPS) methods have demonstrated significant compression of isothermal turbulent fields, their application to buoyancy-driven flows with active thermal coupling has remained unexplored. We apply MPS to two-dimensional Rayleigh-B\'enard convection with dynamical simulations up to $\mathrm{Ra} = 10^{10}$. An a priori decomposition of DNS snapshots up to $\mathrm{Ra} = 10^{11}$ shows that the bond dimension $\chi$ required to represent the flow fields grows without saturation, in contrast to the plateauing of $\chi$ reported for velocity fields in isothermal 2D turbulence. Crucially, however, dynamical simulations solving the governing equations directly in the compressed MPS format at fixed $\chi$ show that the $\chi$ required to recover statistical observables, such as the Nusselt number, scales significantly more favorably with $\mathrm{Ra}$ than the a priori complexity suggests. At $\mathrm{Ra} = 10^{10}$, a relative error of $1.8\%$ in the mean Nusselt number is achieved with a nearly 9-fold reduction in degrees of freedom, using a $\chi$ comparable to that required at $\mathrm{Ra} = 10^{9}$. Spectral analysis confirms the progressive recovery of spatial and temporal scales with increasing $\chi$. These findings establish MPS as a scalable tool for simulating thermally driven turbulence, suggesting the method may remain viable for investigations of the ultimate regime at substantially higher $\mathrm{Ra}$.