Forecasting small-scale dynamics of fluid turbulence using deep neural networks.

Turbulence in fluid flows is characterized by a wide range of interacting scales. Since the scale range increases as some power of the flow Reynolds number, a faithful simulation of the entire scale range is prohibitively expensive at high Reynolds numbers. The most expensive aspect concerns the small-scale motions; thus, major emphasis is placed on understanding and modeling them, taking advantage of their putative universality. In this work, using physics-informed deep learning methods, we present a modeling framework to capture and predict the small-scale dynamics of turbulence, via the velocity gradient tensor. The model is based on obtaining functional closures for the pressure Hessian and viscous Laplacian contributions as functions of velocity gradient tensor. This task is accomplished using deep neural networks that are consistent with physical constraints and explicitly incorporate Reynolds number dependence to account for small-scale intermittency. We then utilize a massive direct numerical simulation database, spanning two orders of magnitude in the large-scale Reynolds number, for training and validation. The model learns from low to moderate Reynolds numbers and successfully predicts velocity gradient statistics at both seen and higher (unseen) Reynolds numbers. The success of our present approach demonstrates the viability of deep learning over traditional modeling approaches in capturing and predicting small-scale features of turbulence.

Saturation and Multifractality of Lagrangian and Eulerian Scaling Exponents in Three-Dimensional Turbulence.

Inertial-range scaling exponents for both Lagrangian and Eulerian structure functions are obtained from direct numerical simulations of isotropic turbulence in triply periodic domains at Taylor-scale Reynolds number up to 1300. We reaffirm that transverse Eulerian scaling exponents saturate at ≈2.1 for moment orders p≥10, significantly differing from the longitudinal exponents (which are predicted to saturate at ≈7.3 for p≥30 from a recent theory). The Lagrangian scaling exponents likewise saturate at ≈2 for p≥8. The saturation of Lagrangian exponents and transverse Eulerian exponents is related by the same multifractal spectrum by utilizing the well-known frozen hypothesis to relate spatial and temporal scales. Furthermore, this spectrum is different from the known spectra for Eulerian longitudinal exponents, suggesting that Lagrangian intermittency is characterized solely by transverse Eulerian intermittency. We discuss possible implications of this outlook when extending multifractal predictions to the dissipation range, especially for Lagrangian acceleration.

Vorticity-Strain Rate Dynamics and the Smallest Scales of Turbulence.

Building upon the intrinsic properties of Navier-Stokes dynamics, namely the prevalence of intense vortical structures and the interrelationship between vorticity and strain rate, we propose a simple framework to quantify the extreme events and the smallest scales of turbulence. We demonstrate that our approach is in excellent agreement with the best available data from direct numerical simulations of isotropic turbulence, with Taylor-scale Reynolds numbers up to 1300. We additionally highlight a shortcoming of prevailing intermittency models due to their disconnection from the observed correlation between vorticity and strain. Our work accentuates the importance of this correlation as a crucial step in developing an accurate understanding of intermittency in turbulence.

Scaling of Acceleration Statistics in High Reynolds Number Turbulence.

The scaling of acceleration statistics in turbulence is examined by combining data from the literature with new data from well-resolved direct numerical simulations of isotropic turbulence, significantly extending the Reynolds number range. The acceleration variance at higher Reynolds numbers departs from previous predictions based on multifractal models, which characterize Lagrangian intermittency as an extension of Eulerian intermittency. The disagreement is even more prominent for higher-order moments of the acceleration. Instead, starting from a known exact relation, we relate the scaling of acceleration variance to that of Eulerian fourth-order velocity gradient and velocity increment statistics. This prediction is in excellent agreement with the variance data. Our Letter highlights the need for models that consider Lagrangian intermittency independent of the Eulerian counterpart.

Generation of intense dissipation in high Reynolds number turbulence.

Intense fluctuations of energy dissipation rate in turbulent flows result from the self-amplification of strain rate via a quadratic nonlinearity, with contributions from vorticity (via the vortex stretching mechanism) and pressure-Hessian-which are analysed here using direct numerical simulations of isotropic turbulence on up to [Formula: see text] grid points, and Taylor-scale Reynolds numbers in the range 140-1300. We extract the statistics involved in amplification of strain and condition them on the magnitude of strain. We find that strain is self-amplified by the quadratic nonlinearity, and depleted via vortex stretching, whereas pressure-Hessian acts to redistribute strain fluctuations towards the mean-field and hence depletes intense strain. Analysing the intense fluctuations of strain in terms of its eigenvalues reveals that the net amplification is solely produced by the third eigenvalue, resulting in strong compressive action. By contrast, the self-amplification acts to deplete the other two eigenvalues, whereas vortex stretching acts to amplify them, with both effects cancelling each other almost perfectly. The effect of the pressure-Hessian for each eigenvalue is qualitatively similar to that of vortex stretching, but significantly weaker in magnitude. Our results conform with the familiar notion that intense strain is organized in sheet-like structures, which are in the vicinity of, but never overlap with tube-like regions of intense vorticity due to fundamental differences in their amplifying mechanisms. This article is part of the theme issue 'Scaling the turbulence edifice (part 1)'.

Turbulence is an Ineffective Mixer when Schmidt Numbers Are Large.

We solve the advection-diffusion equation for a stochastically stationary passive scalar θ, in conjunction with forced 3D Navier-Stokes equations, using direct numerical simulations in periodic domains of various sizes, the largest being 8192^{3}. The Taylor-scale Reynolds number varies in the range 140-650 and the Schmidt number Sc≡ν/D in the range 1-512, where ν is the kinematic viscosity of the fluid and D is the molecular diffusivity of θ. Our results show that turbulence becomes an ineffective mixer when Sc is large. First, the mean scalar dissipation rate ⟨χ⟩=2D⟨|∇θ|^{2}⟩, when suitably nondimensionalized, decreases as 1/logSc. Second, 1D cuts through the scalar field indicate increasing density of sharp fronts on larger scales, oscillating with large excursions leading to reduced mixing, and additionally suggesting weakening of scalar variance flux across the scales. The scaling exponents of the scalar structure functions in the inertial-convective range appear to saturate with respect to the moment order and the saturation exponent approaches unity as Sc increases, qualitatively consistent with 1D cuts of the scalar.

Small-Scale Isotropy and Ramp-Cliff Structures in Scalar Turbulence.

Passive scalars advected by three-dimensional Navier-Stokes turbulence exhibit a fundamental anomaly in odd-order moments because of the characteristic ramp-cliff structures, violating small-scale isotropy. We use data from direct numerical simulations with grid resolution of up to 8192^{3} at high Péclet numbers to understand this anomaly as the scalar diffusivity, D, diminishes, or as the Schmidt number, Sc=ν/D, increases; here ν is the kinematic viscosity of the fluid. The microscale Reynolds number varies from 140 to 650 and Sc varies from 1 to 512. A simple model for the ramp-cliff structures is developed and shown to characterize the scalar derivative statistics very well. It accurately captures how the small-scale isotropy is restored in the large-Sc limit, and additionally suggests a possible correction to the Batchelor length scale as the relevant smallest scale in the scalar field.

Twisting vortex lines regularize Navier-Stokes turbulence.

Fluid flows are intrinsically characterized via the topology and dynamics of underlying vortex lines. Turbulence in common fluids like water and air, mathematically described by the incompressible Navier-Stokes equations (INSE), engenders spontaneous self-stretching and twisting of vortex lines, generating a complex hierarchy of structures. While the INSE are routinely used to describe turbulence, their regularity remains unproven; the implicit assumption being that the self-stretching is ultimately regularized by viscosity, preventing any singularities. Here, we uncover an inviscid regularizing mechanism stemming from self-stretching itself, by analyzing the flow topology as perceived by an observer aligned with the vorticity vector undergoing amplification. While, initially, vorticity amplification occurs via increasing twisting of vortex lines, a regularizing anti-twist spontaneously emerges to prevent unbounded growth. By isolating a vortex, we additionally demonstrate the genericity of this self-regularizing anti-twist. Our work, directly linking dynamics of vortices to turbulence statistics, reveals how the Navier-Stokes dynamics avoids the development of singularities even without the aid of viscosity.

Self-attenuation of extreme events in Navier-Stokes turbulence.

Turbulent fluid flows are ubiquitous in nature and technology, and are mathematically described by the incompressible Navier-Stokes equations. A hallmark of turbulence is spontaneous generation of intense whirls, resulting from amplification of the fluid rotation-rate (vorticity) by its deformation-rate (strain). This interaction, encoded in the non-linearity of Navier-Stokes equations, is non-local, i.e., depends on the entire state of the flow, constituting a serious hindrance in turbulence theory and even establishing regularity of the equations. Here, we unveil a novel aspect of this interaction, by separating strain into local and non-local contributions utilizing the Biot-Savart integral of vorticity in a sphere of radius R. Analyzing highly-resolved numerical turbulent solutions to Navier-Stokes equations, we find that when vorticity becomes very large, the local strain over small R surprisingly counteracts further amplification. This uncovered self-attenuation mechanism is further shown to be connected to local Beltramization of the flow, and could provide a direction in establishing the regularity of Navier-Stokes equations.