Multiscale Turbulence Synthesis: Validation in 2D Hydrodynamics

TL;DR

This paper introduces MuScaTS, a multiscale turbulence synthesis method that generates realistic 2D turbulent flow fields efficiently, bridging the gap between simple models and costly numerical simulations.

Contribution

The paper presents a novel multiscale approach for synthesizing turbulence that accurately reproduces statistical properties and non-Gaussian structures in 2D hydrodynamics with reduced computational cost.

Findings

Reproduces power spectra, increments, and structure functions up to one-third of the turnover time.

Uses scattering transform statistics to characterize non-Gaussian features.

Scales logarithmically with resolution, offering computational efficiency.

Abstract

Numerical simulations can follow the evolution of fluid motions through the intricacies of developed turbulence. However, they are rather costly to run, especially in 3D. In the past two decades, generative models have emerged which produce synthetic random flows at a computational cost equivalent to no more than a few time-steps of a simulation. These simplified models qualitatively bear some characteristics of turbulent flows in specific contexts (incompressible 3D hydrodynamics or magnetohydrodynamics), but generally struggle with the synthesis of coherent structures. We aim at generating random fields (e.g. velocity, density, magnetic fields, etc.) with realistic physical properties for a large variety of governing partial differential equations and at a small cost relative to time-resolved simulations. We propose a set of approximations applied to given sets of partial differential equations, and test the validity of our method in the simplest framework: 2D decaying incompressible hydrodynamical turbulence. We compare results of 2D decaying simulations with snapshots of our synthetic turbulence. We assess quantitatively the difference first with standard statistical tools: power spectra, increments and structure functions. These indicators can be reproduced by our method during up to about a third of the turnover time scale. We also consider recently developed scattering transforms statistics, able to efficiently characterise non-Gaussian structures. This reveals more significant discrepancy, which can however be bridged by bootstrapping. Finally, the number of Fourier transforms necessary for one synthesis scales logarithmically in the resolution, compared to linearly for time-resolved simulations. We have designed a multiscale turbulence synthesis (MuScaTS) method to efficiently short-circuit costly numerical simulations to produce realistic instantaneous fields.