L. Biferale, F. Bonaccorso, M. Buzzicotti, P. Clark di Leoni
We present TURB-Rot, a new open database of 3d and 2d snapshots of turbulent velocity fields, obtained by Direct Numerical Simulations (DNS) of the original Navier-Stokes equations in the presence of rotation. The aim is to provide the community interested in data-assimilation and/or computer vision with a new testing-ground made of roughly 300K complex images and fields. TURB-Rot data are characterized by multi-scales strongly non-Gaussian features and rough, non-differentiable, fields over almost two decades of scales. In addition, coming from fully resolved numerical simulations of the original partial differential equations, they offer the possibility to apply a wide range of approaches, from equation-free to physics-based models. TURB-Rot data are reachable at smart-turb.roma2.infn.it
L Biferale, F Bonaccorso, M Buzzicotti, PC Di Leoni, K Gustavsson
To find the path that minimizes the time to navigate between two given points in a fluid flow is known as Zermelo's problem. Here, we investigate it by using a Reinforcement Learning (RL) approach for the case of a vessel which has a slip velocity with fixed intensity, Vs , but variable direction and navigating in a 2D turbulent sea. We show that an Actor-Critic RL algorithm is able to find quasi-optimal solutions for both time-independent and chaotically evolving flow configurations. For the frozen case, we also compared the results with strategies obtained analytically from continuous Optimal Navigation (ON) protocols. We show that for our application, ON solutions are unstable for the typical duration of the navigation process, and are therefore not useful in practice. On the other hand, RL solutions are much more robust with respect to small changes in the initial conditions and to external noise, even when V s is much smaller than the maximum flow velocity. Furthermore, we show how the RL approach is able to take advantage of the flow properties in order to reach the target, especially when the steering speed is small
L Biferale, F Bonaccorso, M Buzzicotti, KP Iyer
A class of spectral subgrid models based on a self-similar and reversible closure is studied with the aim to minimize the impact of subgrid scales on the inertial range of fully developed turbulence. In this manner, we improve the scale extension where anomalous exponents are measured by roughly 1 order of magnitude when compared to direct numerical simulations or to other popular subgrid closures at the same resolution. We find a first indication that intermittency for high-order moments is not captured by many of the popular phenomenological models developed so far
M Buzzicotti, L Biferale, F Toschi
By means of high-resolution numerical simulations, we compare the statistical properties of homogeneous and isotropic turbulence to those of the Navier-Stokes equation where small-scale vortex filaments are strongly depleted, thanks to a non-linear extra viscosity acting preferentially on high vorticity regions. We show that the presence of such smart small-scale drag can strongly reduce intermittency and non-Gaussian fluctuations. Our results pave the way towards a deeper understanding on the fundamental role of degrees of freedom in turbulence as well as on the impact of (pseudo)coherent structures on the statistical small-scale properties. Our work can be seen as a first attempt to develop smart-Lagrangian forcing/drag mechanisms to control turbulence
P. Clark Di Leoni, A. Mazzino, L. Biferale
Nudging is an important data assimilation technique where partial field measurements are used to control the evolution of a dynamical system and/or to reconstruct the entire phase-space configuration of the supplied flow. Here, we apply it to the toughest problem in fluid dynamics: three dimensional homogeneous and isotropic turbulence. By doing numerical experiments we perform a systematic assessment of how well the technique reconstructs large- and small-scales features of the flow with respect to the quantity and the quality/type of data supplied to it. The types of data used are: (i) field values on a fixed number of spatial locations (Eulerian nudging), (ii) Fourier coefficients of the fields on a fixed range of wavenumbers (Fourier nudging), or (iii) field values along a set of moving probes inside the flow (Lagrangian nudging). We present state-of-the-art quantitative measurements of the scale-by-scale {\it transition to synchronization} and a detailed discussion of the probability distribution function of the reconstruction error, by comparing the nudged field and the {\it truth} point-by-point. Furthermore, we show that for more complex flow configurations, like the case of anisotropic rotating turbulence, the presence of cyclonic and anticyclonic structures leads to unexpectedly better performances of the algorithm. We discuss potential further applications of nudging to a series of applied flow configurations, including the problem of field-reconstruction in thermal Rayleigh-Bénard convection and in magnetohydrodynamics (MHD), and to the determination of optimal parametrisation for small-scale turbulent modeling. Our study fixes the standard requirements for future applications of nudging to complex turbulent flows
G. Margazoglou, L. Biferale, R. Grauer, K. Jansen, D. Mesterházy, T. Rosenow, and R. Tripiccione
We introduce a variant of the Hybrid Monte Carlo (HMC) algorithm to address large-deviation statistics in stochastic hydrodynamics. Based on the path-integral approach to stochastic (partial) differential equations, our HMC algorithm samples space-time histories of the dynamical degrees of freedom under the influence of random noise. First, we validate and benchmark the HMC algorithm by reproducing multiscale properties of the one-dimensional Burgers equation driven by Gaussian and white-in-time noise. Second, we show how to implement an importance sampling protocol to significantly enhance, by orders of magnitudes, the probability to sample extreme and rare events, making it possible to estimate moments of field variables of extremely high order (up to 30 and more). By employing reweighting techniques, we map the biased configurations back to the original probability measure in order to probe their statistical importance. Finally, we show that by biasing the system towards very intense negative gradients, the HMC algorithm is able to explore the statistical fluctuations around instanton configurations. Our results will also be interesting and relevant in lattice gauge theory since they provide unique insights into reweighting techniques
L Biferale, D Khomenko, V L'vov, A Pomyalov, I Procaccia, G Sahoo
Three-dimensional anisotropic turbulence in classical fluids tends towards isotropy and homogeneity with decreasing scales, allowing—eventually—the abstract model of homogeneous and isotropic turbulence to be relevant. We show here that the opposite is true for superfluid 4He turbulence in three-dimensional counterflow channel geometry. This flow becomes less isotropic upon decreasing scales, becoming eventually quasi-two-dimensional. The physical reason for this unusual phenomenon is elucidated and supported by theory and simulations
L Biferale, K Gustavsson, R Scatamacchia
We present numerical and theoretical results concerning the properties of turbulent flows with strong multi-scale helical injection. We perform direct numerical simulations of the Navier–Stokes equations under a random helical stirring with power-law spectrum and with different intensities of energy and helicity injections. We show that there exists three different regimes where the forward energy and helicity inertial transfers are: (i) both leading with respect to the external injections, (ii) energy transfer is leading and helicity transfer is sub-leading and (iii) both are sub-leading and helicity is maximal at all scales. As a result, the cases (ii)–(iii) give flows with Kolmogorov-like inertial energy cascade and tuneable helicity transfers/contents. We further explore regime (iii) by studying its effect on the kinetics of point-like isotropic helicoids, particles whose dynamics is isotropic but breaks parity invariance. We investigate small-scale fractal clustering and preferential sampling of intense helical flow structures. Depending on their structural parameters, the isotropic helicoids either preferentially sample co-chiral or anti-chiral flow structures. We explain these findings in limiting cases in terms of what is known for spherical particles of different densities and degrees of inertia. Furthermore, we present theoretical and numerical results for a stochastic model where dynamical properties can be calculated using analytical perturbation theory. Our study shows that a suitable tuning of the stirring mechanism can strongly modify the small-scale turbulent helical properties and demonstrates that isotropic helicoids are the simplest particles able to preferentially sense helical properties in turbulence
A Alexakis, L Biferale
Turbulent flows are characterized by the non-linear cascades of energy and other inviscid invariants across a huge range of scales, from where they are injected to where they are dissipated. Recently, new experimental, numerical and theoretical works have revealed that many turbulent configurations deviate from the ideal three and two dimensional homogeneous and isotropic cases characterized by the presence of a strictly direct and inverse energy cascade, respectively. New phenomena appear that alter the global and local transfer properties. In this review, we provide a critical summary of historical and recent works from a unified point of view and we present a classification of all known transfer mechanisms. Beside the classical cases of direct and inverse energy cascades, the different scenarios include: split cascades for which an invariant flows both to small and large scales simultaneously, multiple/dual cascades of different quantities, bi-directional cascades where direct and inverse transfers of the same invariant coexist in the same scale-range and finally equilibrium states where no cascades are present, including the case when a large scale condensate is formed. We classify all possible transitions from one scenario to another as the control parameters are changed and we analyse when and why different configurations are observed. Our discussion is based on a set of paradigmatic applications: helical turbulence, rotating and/or stratified flows, magnetohydrodynamics (MHD) turbulence, and passive/active scalars where the transfer properties are altered as one changes the embedding dimensions, the thickness of the domain or other relevant control parameters, as, e.g., the Reynolds, Rossby, Froude, Peclet, or Alfven numbers. We briefly discuss the presence of anomalous scaling laws in 3D hydrodynamics and in other configurations, in connection with the intermittent nature of the energy dissipation in configuration space. A quick overview is also provided concerning the importance of cascades in other applications such as bounded flows, quantum fluids, relativistic and compressible turbulence, and active matter, together with a discussion of the implications for turbulent modelling. Finally, we present a series of open problems and challenges that future work needs to address
S Colabrese, K Gustavsson, A Celani, L Biferale
We performed a numerical study to train smart inertial particles to target specific flow regions with high vorticity through the use of reinforcement learning algorithms. The particles are able to actively change their size to modify their inertia and density. In short, using local measurements of the flow vorticity, the smart particle explores the interplay between its choices of size and its dynamical behavior in the flow environment. This allows it to accumulate experience and learn approximately optimal strategies of how to modulate its size in order to reach the target high-vorticity regions. We consider flows with different complexities: a two-dimensional stationary Taylor-Green-like configuration, a two-dimensional time-dependent flow, and finally a three-dimensional flow given by the stationary Arnold-Beltrami-Childress (ABC) helical flow. We show that smart particles are able to learn how to reach extremely intense vortical structures in all the tackled cases
P Clark Di Leoni, A Mazzino, L Biferale
Inferring physical parameters of turbulent flows by assimilation of data measurements is an open challenge with key applications in meteorology, climate modeling and astrophysics. Up to now, spectral nudging was applied for empirical data-assimilation as a mean to improve deterministic and statistical predictability in the presence of a restricted set of field measurements only. Here, we explore under which conditions a nudging protocol can be used for two novel objectives: to unravel the value of the physical flow parameters and to reconstruct large-scale turbulent properties starting from a sparse set of information in space and in time. First, we apply nudging to quantitatively infer the unknown rotation rate and the shear mechanism for turbulent flows. Second, we show that a suitable spectral nudging is able to reconstruct the energy containing scales in rotating turbulence by using a blind set-up, ie without any input about the external forcing mechanisms acting on the flow. Finally, we discuss the broad potentialities of nudging to other key applications for physics-informed data-assimilation in environmental or applied flow configurations
M Buzzicotti, H Aluie, L Biferale, and M Linkmann
It is known that rapidly rotating turbulent flows are characterized by the emergence of simultaneous upscale and downscale energy transfer. Indeed, both numerics and experiments show the formation of large-scale anisotropic vortices together with the development of small-scale dissipative structures. However the organization of interactions leading to this complex dynamics remains unclear. Two different mechanisms are known to be able to transfer energy upscale in a turbulent flow. The first is characterized by two-dimensional interactions among triads lying on the two-dimensional, three-component (2D3C)/slow manifold, namely on the Fourier plane perpendicular to the rotation axis. The second mechanism is three-dimensional and consists of interactions between triads with the same sign of helicity (homochiral). Here, we present a detailed numerical study of rotating flows using a suite of high-Reynolds-number direct numerical simulations (DNS) within different parameter regimes to analyze both upscale and downscale cascade ranges. We find that the upscale cascade at wave numbers close to the forcing scale is generated by increasingly dominant homochiral interactions which couple the three-dimensional bulk and the 2D3C plane. This coupling produces an accumulation of energy in the 2D3C plane, which then transfers energy to smaller wave numbers thanks to the two-dimensional mechanism. In the forward cascade range, we find that the energy transfer is dominated by heterochiral triads and is dominated primarily by interaction within the fast manifold where kz≠0. We further analyze the energy transfer in different regions in the real-space domain. In particular, we distinguish high-strain from high-vorticity regions and we uncover that while the mean transfer is produced inside regions of strain, the rare but extreme events of energy transfer occur primarily inside the large-scale column vortices
M Buzzicotti, M Linkmann, H Aluie, L Biferale, J Brasseur and C Meneveau
The effects of different filtering strategies on the statistical properties of the resolved-to-subfilter scale (SFS) energy transfer are analysed in forced homogeneous and isotropic turbulence. We carry out a-priori analyses of the statistical characteristics of SFS energy transfer by filtering data obtained from direct numerical simulations with up to 20483 grid points as a function of the filter cutoff scale. In order to quantify the dependence of extreme events and anomalous scaling on the filter, we compare a sharp Fourier Galerkin projector, a Gaussian filter and a novel class of Galerkin projectors with non-sharp spectral filter profiles. Of interest is the importance of Galilean invariance and we confirm that local SFS energy transfer displays intermittency scaling in both skewness and flatness as a function of the cutoff scale. Furthermore, we quantify the robustness of scaling as a function of the filtering type
K P Iyer, F Bonaccorso, L Biferale, F Toschi
We use high-resolution direct numerical simulations to study the anisotropic contents of a turbulent, statistically homogeneous flow with random transitions among multiple energy containing states. We decompose the velocity correlation functions on different sectors of the three-dimensional group of rotations, SO(3), using a high-precision quadrature. Scaling properties of anisotropic components of longitudinal and transverse velocity fluctuations are accurately measured at changing Reynolds numbers. We show that independently of the anisotropic content of the energy containing eddies, small-scale turbulent fluctuations recover isotropy and universality faster than previously reported in experimental and numerical studies. The discrepancies are ascribed to the presence of highly anisotropic contributions that have either been neglected or measured with less accuracy in the foregoing works. Furthermore, the anomalous anisotropic scaling exponents are devoid of any sign of saturation with increasing order. Our study paves the way to systematically assess persistence of anisotropy in high-Reynolds-number flows
L Biferale, AA Mailybaev, G Parisi
We discuss a theoretical framework to define an optimal subgrid closure for shell models of turbulence. The closure is based on the ansatz that consecutive shell multipliers are short-range correlated, following the third hypothesis of Kolmogorov formulated for similar quantities for the original three-dimensional Navier-Stokes turbulence. We also propose a series of systematic approximations to the optimal model by assuming different degrees of correlations across scales among amplitudes and phases of consecutive multipliers. We show numerically that such low-order closures work well, reproducing all known properties of the large-scale dynamics including anomalous scaling. We found small but systematic discrepancies only for a range of scales close to the subgrid threshold, which do not tend to disappear by increasing the order of the approximation. We speculate that the lack of convergence might be due to a structural instability, at least for the evolution of very fast degrees of freedom at small scales. Connections with similar problems for large eddy simulations of the three-dimensional Navier-Stokes equations are also discussed
G Sahoo, A Alexakis, L Biferale
Inviscid invariants of flow equations are crucial in determining the direction of the turbulent energy cascade. In this work we investigate a variant of the three-dimensional Navier-Stokes equations that shares exactly the same ideal invariants (energy and helicity) and the same symmetries (under rotations, reflections, and scale transforms) as the original equations. It is demonstrated that the examined system displays a change in the direction of the energy cascade when varying the value of a free parameter which controls the relative weights of the triadic interactions between different helical Fourier modes. The transition from a forward to inverse cascade is shown to occur at a critical point in a discontinuous manner with diverging fluctuations close to criticality. Our work thus supports the observation that purely isotropic and three-dimensional flow configurations can support inverse energy transfer when interactions are altered and that inside all turbulent flows there is a competition among forward and backward transfer mechanisms which might lead to multiple energy-containing turbulent states
S Colabrese, K Gustavsson, A Celani, L Biferale
Smart active particles can acquire some limited knowledge of the fluid environment from simple mechanical cues and exert a control on their preferred steering direction. Their goal is to learn the best way to navigate by exploiting the underlying flow whenever possible. As an example, we focus our attention on smart gravitactic swimmers. These are active particles whose task is to reach the highest altitude within some time horizon, given the constraints enforced by fluid mechanics. By means of numerical experiments, we show that swimmers indeed learn nearly optimal strategies just by experience. A reinforcement learning algorithm allows particles to learn effective strategies even in difficult situations when, in the absence of control, they would end up being trapped by flow structures. These strategies are highly nontrivial and cannot be easily guessed in advance. This Letter illustrates the potential of reinforcement learning algorithms to model adaptive behavior in complex flows and paves the way towards the engineering of smart microswimmers that solve difficult navigation problems
L Biferale, F Bonaccorso, I M Mazzitelli, M AT van Hinsberg, A S Lanotte, S Musacchio, P Perlekar, F Toschi
By using direct numerical simulations (DNS) at unprecedented resolution, we study turbulence under rotation in the presence of simultaneous direct and inverse cascades. The accumulation of energy at large scale leads to the formation of vertical coherent regions with high vorticity oriented along the rotation axis. By seeding the flow with millions of inertial particles, we quantify—for the first time—the effects of those coherent vertical structures on the preferential concentration of light and heavy particles. Furthermore, we quantitatively show that extreme fluctuations, leading to deviations from a normal-distributed statistics, result from the entangled interaction of the vertical structures with the turbulent background. Finally, we present the first-ever measurement of the relative importance between Stokes drag, Coriolis force, and centripetal force along the trajectories of inertial particles. We discover that vortical coherent structures lead to unexpected diffusion properties for heavy and light particles in the directions parallel and perpendicular to the rotation axis