Numerical simulations of a stochastic dynamics leading to cascades and loss of regularity: applications to fluid turbulence and generation of fractional Gaussian fields

Geoffrey Beck,Charles-Edouard Bréhier,Laurent Chevillard,Ricardo Grande,Wandrille Ruffenach
2024-03-08
Abstract:Motivated by the modeling of the spatial structure of the velocity field of three-dimensional turbulent flows, and the phenomenology of cascade phenomena, a linear dynamics has been recently proposed able to generate high velocity gradients from a smooth-in-space forcing term. It is based on a linear Partial Differential Equation (PDE) stirred by an additive random forcing term which is delta-correlated in time. The underlying proposed deterministic mechanism corresponds to a transport in Fourier space which aims at transferring energy injected at large scales towards small scales. The key role of the random forcing is to realize these transfers in a statistically homogeneous way. Whereas at finite times and positive viscosity the solutions are smooth, a loss of regularity is observed for the statistically stationary state in the inviscid limit. We here present novel simulations, based on finite volume methods in the Fourier domain and a splitting method in time, which are more accurate than the pseudo-spectral simulations. We show that the novel algorithm is able to reproduce accurately the expected local and statistical structure of the predicted solutions. We conduct numerical simulations in one, two and three spatial dimensions, and we display the solutions both in physical and Fourier spaces. We additionally display key statistical quantities such as second-order structure functions and power spectral densities at various viscosities.
Fluid Dynamics,Numerical Analysis,Computational Physics
What problem does this paper attempt to address?
This paper mainly discusses a stochastic dynamical model for simulating energy cascade and instability phenomena in three-dimensional turbulent fluids. This linear partial differential equation (PDE) model is driven by a time-correlated time-additive stochastic forcing term, which generates high velocity gradients in space, transferring energy from smooth initial conditions to small scales. Although solutions are smooth for finite time and positive viscosity, instability phenomena are observed in the limit of zero viscosity, resulting in statistical equilibrium states. The paper presents more accurate numerical simulations using finite volume methods and time-splitting techniques, which reproduce the local and statistical structures of the predicted solutions, including solutions in one, two, and three-dimensional spaces as well as solutions in Fourier space. The study also highlights that the model is able to simulate second-order statistical behaviors, such as viscosity-independent mean square velocity, power spectral density with power-law decay, and power-law behavior of second-order structure functions, which are fundamental characteristics of turbulence.