An empirical mean-field model of symmetry-breaking in a turbulent wake.

Improved turbulence modeling remains a major open problem in mathematical physics. Turbulence is notoriously challenging, in part due to its multiscale nature and the fact that large-scale coherent structures cannot be disentangled from small-scale fluctuations. This closure problem is emblematic of a greater challenge in complex systems, where coarse-graining and statistical mechanics descriptions break down. This work demonstrates an alternative data-driven modeling approach to learn nonlinear models of the coherent structures, approximating turbulent fluctuations as state-dependent stochastic forcing. We demonstrate this approach on a high-Reynolds number turbulent wake experiment, showing that our model reproduces empirical power spectra and probability distributions. The model is interpretable, providing insights into the physical mechanisms underlying the symmetry-breaking behavior in the wake. This work suggests a path toward low-dimensional models of globally unstable turbulent flows from experimental measurements, with broad implications for other multiscale systems.

Dimensionally consistent learning with Buckingham Pi.

In the absence of governing equations, dimensional analysis is a robust technique for extracting insights and finding symmetries in physical systems. Given measurement variables and parameters, the Buckingham Pi theorem provides a procedure for finding a set of dimensionless groups that spans the solution space, although this set is not unique. We propose an automated approach using the symmetric and self-similar structure of available measurement data to discover the dimensionless groups that best collapse these data to a lower dimensional space according to an optimal fit. We develop three data-driven techniques that use the Buckingham Pi theorem as a constraint: (1) a constrained optimization problem with a non-parametric input-output fitting function, (2) a deep learning algorithm (BuckiNet) that projects the input parameter space to a lower dimension in the first layer and (3) a technique based on sparse identification of nonlinear dynamics to discover dimensionless equations whose coefficients parameterize the dynamics. We explore the accuracy, robustness and computational complexity of these methods and show that they successfully identify dimensionless groups in three example problems: a bead on a rotating hoop, a laminar boundary layer and Rayleigh-Bénard convection.

Learning dominant physical processes with data-driven balance models.

Throughout the history of science, physics-based modeling has relied on judiciously approximating observed dynamics as a balance between a few dominant processes. However, this traditional approach is mathematically cumbersome and only applies in asymptotic regimes where there is a strict separation of scales in the physics. Here, we automate and generalize this approach to non-asymptotic regimes by introducing the idea of an equation space, in which different local balances appear as distinct subspace clusters. Unsupervised learning can then automatically identify regions where groups of terms may be neglected. We show that our data-driven balance models successfully delineate dominant balance physics in a much richer class of systems. In particular, this approach uncovers key mechanistic models in turbulence, combustion, nonlinear optics, geophysical fluids, and neuroscience.

Population annealing simulations of a binary hard-sphere mixture.

Population annealing is a sequential Monte Carlo scheme well suited to simulating equilibrium states of systems with rough free energy landscapes. Here we use population annealing to study a binary mixture of hard spheres. Population annealing is a parallel version of simulated annealing with an extra resampling step that ensures that a population of replicas of the system represents the equilibrium ensemble at every packing fraction in an annealing schedule. The algorithm and its equilibration properties are described, and results are presented for a glass-forming fluid composed of a 50/50 mixture of hard spheres with diameter ratio of 1.4:1. For this system, we obtain precise results for the equation of state in the glassy regime up to packing fractions φ≈0.60 and study deviations from the Boublik-Mansoori-Carnahan-Starling-Leland equation of state. For higher packing fractions, the algorithm falls out of equilibrium and a free volume fit predicts jamming at packing fraction φ≈0.667. We conclude that population annealing is an effective tool for studying equilibrium glassy fluids and the jamming transition.