Data-driven discovery of the equations of turbulent convection

TL;DR

This paper compares SINDy and SPIDER algorithms for discovering governing equations of turbulent convection from DNS data, highlighting SPIDER's advantages in recovering equations, boundary conditions, and constraints across different flow regimes.

Contribution

The study demonstrates that SPIDER outperforms pySINDy in recovering complex PDEs, boundary conditions, and constraints in turbulent convection simulations.

Findings

SPIDER requires fewer library terms and recovers equations more easily.

pySINDy can recover equations at lower Rayleigh numbers but struggles at higher ones.

Flow properties like correlation time inform optimal subdomain sizes for analysis.

Abstract

We compare the efficiency and ease-of-use of the Sparse Identification of Nonlinear Dynamics (SINDy) algorithm and Sparse Physics-Informed Discovery of Empirical Relations (SPIDER) framework in recovering the relevant governing equations and boundary conditions from data generated by direct numerical simulations (DNS) of turbulent convective flows. In the former case, a weak-form implementation pySINDy is used. Time-dependent data for two- (2D) and three-dimensional (3D) DNS simulation of Rayleigh-Benard convection and convective plane Couette flow is generated using the Dedalus PDE framework for spectrally solving differential equations. Using pySINDy we are able to recover the governing equations of 2D models of Rayleigh-Benard convection at Rayleigh numbers, R, from laminar, through transitional to moderately turbulent flow conditions, albeit with increasing difficulty with larger Rayleigh number, especially in recovery of the diffusive terms (with coefficient magnitude proportional to 1/R^0.5). SPIDER requires a much smaller library of terms and we are able to recover more easily the governing equations for a wider range of R in 2D and 3D convection and plane flow models and go on to recover constraints (the incompressibility condition) and boundary conditions, demonstrating the benefits and capabilities of SPIDER to go beyond pySINDy for these fluid problems governed by second-order PDEs. [We] demonstrat[e] the potential of machine-learning methods to validate numerical solvers and solutions for such flow problems. We also find that properties of the flow, specifically the correlation time and spatial scales, should inform the initial selection of spatiotemporal subdomain sizes [abbreviated]