BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Luca Biferale

Luca Biferale appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynnlin.CDcond-mat.stat-mechcond-mat.softMachine Learningeess.SYnlin.CGphysics.geo-phSystems and ControlArtificial IntelligenceBiological Physicscond-mat.mes-hallcond-mat.otherhep-latmath-phmath.APmath.MPMultiagent Systemsphysics.ao-ph

Catalog footprint

What is connected

45works
19topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2026arXiv

Multi-Scale Data Assimilation in Turbulent Models

We explore the potential of Data-Assimilation (DA) within the multi-scale framework of a shell model of turbulence, with a focus on the Ensemble Kalman Filter (EnKF). The central objective is to understand how measuring mesoscales (i.e., inertial-range scales) enhances the prediction of both large-scale and small-scale intermittent variables, by systematically varying observation frequency and the set of measured scales. We demonstrate that measurements conducted at frequencies that exceed those of the observed scales enable full synchronization of larger scales, provided that at least two adjacent mesoscale are measured. In addition, we benchmark the EnKF against two other DA methods, namely Nudging and Ensemble 4D-Var. EnKF is clearly superior to the former, and comparable with the latter but achieving the result with a lower computational complexity. Moreover, our results underscore the need for a tailored, scale-aware inflation technique to stabilize the assimilation process, preventing filter divergence and ensuring robust convergence.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A minimal phase-coupling model for intermittency in turbulent systems

Turbulent systems exhibit a remarkable multi-scale complexity, in which spatial structures induce scale-dependent statistics with strong departures from Gaussianity. In Fourier space, this is reflected by pronounced phase synchronization. A quantitative relation between real-space structure, statistics, and phase synchronization is currently missing. Here, we address this problem in the framework of a minimal phase-coupling model, which enables a detailed investigation by means of dynamical systems theory and multi-scale high-resolution simulations. We identify the spectral power-law steepness, which controls the phase coupling, as the control parameter for tuning the non-Gaussian properties of the system. Whereas both very steep and very shallow spectra exhibit close-to-Gaussian statistics, the strongest departures are observed for intermediate slopes comparable to the ones in hydrodynamic and Burgers turbulence. We show that the non-Gaussian regime of the model coincides with a collapse of the dynamical system to a lower-dimensional attractor and the emergence of phase synchronization, thereby establishing a dynamical-systems perspective on turbulent intermittency.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mesoscale Modelling of the Tolman Length in Multi-component Systems

In this paper we analyze the curvature corrections to the surface tension in the context of the Shan-Chen (SC) multi-component Lattice Boltzmann method (LBM). We demonstrate that the same techniques recently applied in the context of the Shan-Chen multi-phase model can be applied to multi-component mixtures. We implement, as a new application, the calculation of the surface of tension radius $R_s$ through the minimization of the generalized surface tension $σ[R]$. In turn we are able to estimate the Tolman length, i.e. the first order coefficient of the curvature expansion of the surface tension $σ(R)$, as well as the higher order corrections, i.e. the curvature- and the Gaussian-rigidity coefficients. The SC multi-component model allows to model both fully-symmetric as well as asymmetric interactions among the components. By performing an extensive set of simulations we present a first example of tunable Tolman length in the mesoscopic model, being zero for symmetric interactions and different from zero otherwise. This result paves the way for controlling such interface properties which are paramount in presence of thermal fluctuations. All reported results can be independently reproduced through the "idea.deploy" framework available at https://github.com/lullimat/idea.deploy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Non-equilibrium Ensembles for the three-dimensional Navier-Stokes Equations

At the molecular level fluid motions are, by first principles, described by time reversible laws. On the other hand, the coarse grained macroscopic evolution is suitably described by the Navier-Stokes equations, which are inherently irreversible, due to the dissipation term. Here, a reversible version of three-dimensional Navier-Stokes is studied, by introducing a fluctuating viscosity constructed in such a way that enstrophy is conserved, along the lines of the paradigm of microcanonical versus canonical treatment in equilibrium statistical mechanics. Through systematic simulations we attack two important questions: (a) What are the conditions that must be satisfied in order to have a statistical equivalence between the two non-equilibrium ensembles? (b) What is the empirical distribution of the fluctuating viscosity observed by changing the Reynolds number and the number of modes used in the discretization of the evolution equation? The latter point is important also to establish regularity conditions for the reversible equations. We find that the probability to observe negative values of the fluctuating viscosity becomes very quickly extremely small when increasing the effective Reynolds number of the flow in the fully resolved hydro dynamical regime, at difference from what was observed previously.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Reinforcement learning for pursuit and evasion of microswimmers at low Reynolds number

We consider a model of two competing microswimming agents engaged in a pursue-evasion task within a low-Reynolds-number environment. Agents can only perform simple maneuvers and sense hydrodynamic disturbances, which provide ambiguous (partial) information about the opponent's position and motion. We frame the problem as a zero-sum game: The pursuer has to capture the evader in the shortest time, while the evader aims at deferring capture as long as possible. We show that the agents, trained via adversarial reinforcement learning, are able to overcome partial observability by discovering increasingly complex sequences of moves and countermoves that outperform known heuristic strategies and exploit the hydrodynamic environment.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

$λ$-Navier-Stokes turbulence

We investigate numerically the model proposed in Sahoo et al [Phys. Rev. Lett. 118, 164501, (2017)] where a parameter $λ$ is introduced in the Navier-Stokes equations such that the weight of homochiral to heterochiral interactions is varied while preserving all original scaling symmetries and inviscid invariants. Decreasing the value of $λ$ leads to a change in the direction of the energy cascade at a critical value $λ_c \sim 0.3$. In this work, we perform numerical simulations at varying $λ$ in the forward energy cascade range and at changing the Reynolds number $\mathrm{Re}$. We show that for a fixed injection rate, as $λ\to λ_c$, the kinetic energy diverges with a scaling law $\mathcal{E} \propto (λ-λ_c)^{-2/3}$. The energy spectrum is shown to display a larger bottleneck as $λ$ is decreased. The forward heterochiral flux and the inverse homochiral flux both increase in amplitude as $λ_c$ is approached while keeping their difference fixed and equal to the injection rate. As a result, very close to $λ_c$ a stationary state is reached where the two opposite fluxes are of much higher amplitude than the mean flux and large fluctuations are observed. Furthermore, we show that intermittency as $λ_c$ is approached is reduced. The possibility of obtaining a statistical description of regular Navier-Stokes turbulence as an expansion around this newly found critical point is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A Mesoscale Perspective on the Tolman Length

We demonstrate that the multi-phase Shan-Chen lattice Boltzmann method (LBM) yields a curvature dependent surface tension $σ$ as computed from three-dimensional hydrostatic droplets/bubbles simulations. Such curvature dependence is routinely characterized, at first order, by the so-called {\it Tolman length} $δ$. LBM allows to precisely compute $σ$ at the surface of tension $R_s$ and determine the Tolman length from the coefficient of the first order correction. The corresponding values of $δ$ display universality for different equations of state, following a power-law scaling near the critical temperature. The Tolman length has been studied so far mainly via computationally demanding molecular dynamics (MD) simulations or by means of density functional theory (DFT) approaches playing a pivotal role in extending Classical Nucleation Theory. The present results open a new hydrodynamic-compliant mesoscale arena, in which the fundamental role of the Tolman length, alongside real-world applications to cavitation phenomena, can be effectively tackled. All the results can be independently reproduced through the "idea.deploy" framework.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Acceleration of tracer and light particles in compressible homogeneous isotropic turbulence

The accelerations of tracer and light particles in compressible homogeneous isotropic turbulence (CHIT) is investigated by using data from direct numerical simulations (DNS) up to turbulent Mach number $M_t =1$. For tracer particles, the flatness factor of acceleration components, $F_a$, increases gradually for $M_t \in [0.3, 1]$. On the contrary, $F_a$ for light particles develops a maximum around $M_t \sim 0.6$. The PDF of longitudinal acceleration of tracers is increasingly skewed towards the negative value as $M_t$ increases. By contrast, for light particles, the skewness factor of longitudinal acceleration, $S_a$, firstly becomes more negative with the increase of $M_t$, and then goes back to $0$ when $M_t$ is larger than $0.6$. Similarly, differences among tracers and light particles appear also in the zero-crossing time of acceleration correlation. It is argued that all these phenomenons are intimately linked to the flow structures in compression regions, e.g. close to shocklets.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Cascade or not cascade? Energy transfer and elastic effects in active nematics

We numerically study the multi-scale properties of a 2d active gel to address the energy transfer mechanism. We find that activity is able to excite long-ranged distortions of the nematic pattern giving rise to spontaneous laminar flows and to a chaotic regime by further increasing the rate of active energy injection. By means of a scale-to-scale spectral analysis we find that the gel is basically driven by the local balancing between active injection and viscous dissipation, without any signal of non-linear hydrodynamical transfer and turbulent cascades. Furthermore, elasticity may qualitatively play an important role by transferring energy from small to larger scales through nemato-hydrodynamic interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Optimal control of point-to-point navigation in turbulent time-dependent flows using Reinforcement Learning

We present theoretical and numerical results concerning the problem to find the path that minimizes the time to navigate between two given points in a complex fluid under realistic navigation constraints. We contrast deterministic Optimal Navigation (ON) control with stochastic policies obtained by Reinforcement Learning (RL) algorithms. We show that Actor-Critic RL algorithms are able to find quasi-optimal solutions in the presence of either time-independent or chaotically evolving flow configurations. For our application, ON solutions develop unstable behavior within the typical duration of the navigation process, and are therefore not useful in practice. We first explore navigation of turbulent flow using a constant propulsion speed. Based on a discretized phase-space, the propulsion direction is adjusted with the aim to minimize the time spent to reach the target. Further, we explore a case where additional control is obtained by allowing the engine to power off. Exploiting advection of the underlying flow, allows the target to be reached with less energy consumption. In this case, we optimize a linear combination between the total navigation time and the total time the engine is switched off. Our approach can be generalized to other setups, for example, navigation under imperfect environmental forecast or with different models for the moving vessel.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Lattice Boltzmann simulations of droplet breakup in confined and time-dependent flows

We study droplet dynamics and breakup in generic time-dependent flows via a multicomponent Lattice Boltzmann algorithm, with emphasis on flow start up conditions. We first study droplet breakup in a confined oscillatory shear flow via two different protocols. In one set up, we start from an initially spherical droplet and turn on the flow abruptly ("shock method"); in the other protocol, we start from an initially spherical droplet as well, but we progressively increase the amplitude of the flow, by allowing the droplet to relax to the steady state for each increase in amplitude, before increasing the flow amplitude again ("relaxation method"). The two protocols are shown to produce substantially different breakup scenarios. The mismatch between these two protocols is also studied for variations in the flow topology, the degree of confinement and the inertia of the fluid. All results point to the fact that under extreme conditions of confinement the relaxation protocols can drive the droplets into metastable states, which break only for very intense flow amplitudes, but their stability is prone to external perturbations, such as an oscillatory driving force.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Collective olfactory search in a turbulent environment

Finding the distant source of an odor dispersed by a turbulent flow is a vital task for many organisms, either for foraging or for mating purposes. At the level of individual search, animals like moths have developed effective strategies to solve this very difficult navigation problem based on the noisy detection of odor concentration and wind velocity alone. When many individuals concurrently perform the same olfactory search task, without any centralized control, sharing information about the decisions made by the members of the group can potentially increase the performance. But how much of this information is actually valuable and exploitable for the collective task ? Here we show that, in a model of a swarm of agents inspired by moth behavior, there is an optimal way to blend the private information about odor and wind detections with the publicly available information about other agents' heading direction. At optimality, the time required for the first agent to reach the source is essentially the shortest flight time from the departure point to the target. Conversely, agents who discard public information are several fold slower and groups that do not put enough weight on private information perform even worse. Our results then suggest an efficient multi-agent olfactory search algorithm that could prove useful in robotics, for instance in the identification of sources of harmful volatile compounds.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Multi-scale control of active emulsion dynamics

We numerically study the energy transfer in a multi-component $2d$ film, made of an active polar gel and a passive isotropic fluid in presence of surfactant favoring emulsification. We show that by confining the active behavior into the localized component, the typical scale where chemical energy is transformed in mechanical energy can be substantially controlled. Quantitative analysis of kinetic energy spectra and fluxes shows the presence of a multi-scale dynamics due to the existence of a flux induced by the active stress only, without the presence of a turbulent cascade. An increase in the intensity of active doping induces drag reduction due to the competition of elastic and dissipative stresses against active forces. Furthermore we show that a non-homogeneous activity pattern induces localized response, including a modulation of the slip length, opening the way toward the control of active flows by external doping.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The statistical properties of turbulence in the presence of a smart small-scale control

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Zermelo's problem: Optimal point-to-point navigation in 2D turbulent flows using Reinforcement Learning

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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Non-universal behaviour of helical two-dimensional three-component turbulence

The dynamics of two-dimensional three-component (2D3C) flows is relevant to describe the long-time evolution of strongly rotating flows and/or of conducting fluids with a strong mean magnetic field. We show that in the presence of a strong helical forcing, the out-of-plane component ceases to behave as a passive advected quantity and develops a nontrivial dynamics which deeply changes its large-scale properties. We show that a small-scale helicity injection correlates the input on the 2D component with the one on the out-of-plane component. As a result, the third component develops a non-trivial energy transfer. The latter is mediated by homochiral triads, confirming the strong 3D nature of the leading dynamical interactions. In conclusion, we show that the out-of-plane component in a 2D3C flow enjoys strong non-universal properties as a function of the degree of mirror symmetry of the small-scale forcing.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Lagrangian Statistics for Navier-Stokes Turbulence under Fourier-mode reduction: Fractal and Homogeneous Decimations

We study small-scale and high-frequency turbulent fluctuations in three-dimensional flows under Fourier-mode reduction. The Navier-Stokes equations are evolved on a restricted set of modes, obtained as a projection on a fractal or homogeneous Fourier set. We find a strong sensitivity (reduction) of the high-frequency variability of the Lagrangian velocity fluctuations on the degree of mode decimation, similarly to what is already reported for Eulerian statistics. This is quantified by a tendency towards a quasi-Gaussian statistics, i.e., to a reduction of intermittency, at all scales and frequencies. This can be attributed to a strong depletion of vortex filaments and of the vortex stretching mechanism. Nevertheless, we found that Eulerian and Lagrangian ensembles are still connected by a dimensional bridge-relation which is independent of the degree of Fourier-mode decimation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coherent structures and extreme events in rotating multiphase turbulent flows

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 forces 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extended self-similarity in moment-generating-functions in wall-bounded turbulence at high Reynolds number

In wall-bounded turbulence, the moment generating functions (MGFs) of the streamwise velocity fluctuations $\left<\exp(qu_z^+)\right>$ develop power-law scaling as a function of the wall normal distance $z/δ$. Here $u$ is the streamwise velocity fluctuation, $+$ indicates normalization in wall units (averaged friction velocity), $z$ is the distance from the wall, $q$ is an independent variable and $δ$ is the boundary layer thickness. Previous work has shown that this power-law scaling exists in the log-region {\small $3Re_τ^{0.5}\lesssim z^+$, $z\lesssim 0.15δ$}, where $Re_τ$ is the friction velocity-based Reynolds numbers. Here we present empirical evidence that this self-similar scaling can be extended, including bulk and viscosity-affected regions $30<z^+$, $z<δ$, provided the data are interpreted with the Extended-Self-Similarity (ESS), i.e. self-scaling of the MGFs as a function of one reference value, $q_o$. ESS also improves the scaling properties, leading to more precise measurements of the scaling exponents. The analysis is based on hot-wire measurements from boundary layers at $Re_τ$ ranging from $2700$ to $13000$ from the Melbourne High-Reynolds-Number-Turbulent-Boundary-Layer-Wind-Tunnel. Furthermore, we investigate the scalings of the filtered, large-scale velocity fluctuations $u^L_z$ and of the remaining small-scale component, $u^S_z=u_z-u^L_z$. The scaling of $u^L_z$ falls within the conventionally defined log region and depends on a scale that is proportional to {\small $l^+\sim Re_τ^{1/2}$}; the scaling of $u^{S}_z$ extends over a much wider range from $z^+\approx 30$ to $z\approx 0.5δ$. Last, we present a theoretical construction of two multiplicative processes for $u^L_z$ and $u^S_z$ that reproduce the empirical findings concerning the scalings properties as functions of $z^+$ and in the ESS sense.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Intermittency in fractal Fourier hydrodynamics: Lessons from the Burgers Equation

We present theoretical and numerical results for the one-dimensional stochastically forced Burgers equation decimated on a fractal Fourier set of dimension $D$. We investigate the robustness of the energy transfer mechanism and of the small-scale statistical fluctuations by changing $D$. We find that a very small percentage of mode-reduction ($D \lesssim 1$) is enough to destroy most of the characteristics of the original non-decimated equation. In particular, we observe a suppression of intermittent fluctuations for $D <1$ and a quasi-singular transition from the fully intermittent ($D=1$) to the non-intermittent case for $D \lesssim 1$. Our results indicate that the existence of strong localized structures (shocks) in the one-dimensional Burgers equation is the result of highly entangled correlations amongst all Fourier modes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On the Vortex Dynamics in Fractal Fourier Turbulence

Incompressible, homogeneous and isotropic turbulence is studied by solving the Navier-Stokes equations on a reduced set of Fourier modes, belonging to a fractal set of dimension $D$. By tuning the fractal dimension parameter, we study the dynamical effects of Fourier decimation on the vortex stretching mechanism and on the statistics of the velocity and the velocity gradient tensor. In particular, we show that as we move from $D=3$ to $D \sim 2.8$, the statistics gradually turns into a purely Gaussian one. This result suggests that even a mild fractal mode reduction strongly depletes the stretching properties of the non-linear term of the Navier-Stokes equations and suppresses anomalous fluctuations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Phase and precession evolution in the Burgers equation

We present a phenomenological study of the phase dynamics of the one-dimensional stochastically forced Burgers equation, and of the same equation under a Fourier mode reduction on a fractal set. We study the connection between coherent structures in real space and the evolution of triads in Fourier space. Concerning the one-dimensional case, we find that triad phases show alignments and synchronisations that favour energy fluxes towards small scales --a direct cascade. In addition, strongly dissipative real-space structures are associated with entangled correlations amongst the phase precession frequencies and the amplitude evolution of Fourier triads. As a result, triad precession frequencies show a non-Gaussian distribution with multiple peaks and fat tails, and there is a significant correlation between triad precession frequencies and amplitude growth. Links with dynamical systems approach are briefly discussed, such as the role of unstable critical points in state space. On the other hand, by reducing the fractal dimension $D$ of the underlying Fourier set, we observe: i) a tendency toward a more Gaussian statistics, ii) a loss of alignment of triad phases leading to a depletion of the energy flux, and iii) the simultaneous reduction of the correlation between the growth of Fourier mode amplitudes and the precession frequencies of triad phases.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Preferential sampling of helicity by isotropic helicoids

We present a theoretical and numerical study on the motion of isotropic helicoids in complex flows. These are particles whose motion is invariant under rotations but not under mirror reflections of the particle. This is the simplest, yet unexplored, extension of the much studied case of small spherical particles. We show that heavy isotropic helicoids, due to the coupling between translational and rotational degrees of freedom, preferentially sample different helical regions in laminar or chaotic advecting flows. This opens the way to control and engineer particles able to track complex flow structures with potential applications to microfluidics and turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Role of helicity for large- and small-scale turbulent fluctuations

The effect of the helicity on the dynamics of the turbulent flows is investigated. The aim is to disentangle the role of helicity in fixing the direction, the intensity and the fluctuations of the energy transfer across the inertial range of scales. We introduce an external parameter, $α$, that controls the mismatch between the number of positive and negative helically polarized Fourier modes. We present the first set of direct numerical simulations of Navier-Stokes equations from the fully symmetrical case, $α=0$, to the fully asymmetrical case, $α=1$, when only helical modes of one sign survive. We found a singular dependency of the direction of the energy cascade on $α$, measuring a positive forward flux as soon as only a few modes with different helical polarities are present. On the other hand, small-scales fluctuations are sensitive only to the degree of mode-reduction, leading to a vanishing intermittency already for values of $α\sim 0.1$ and independently of the degree of mirror symmetry-breaking. Our findings suggest that intermittency is the result of a global mode-coupling in Fourier space.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rotating turbulence under "precession-like" perturbation

The effects of changing the orientation of the rotation axis on homogeneous turbulence is considered. We perform direct numerical simulations on a periodic box of $1024^3$ grid points, where the orientation of the rotation axis is changed (a) at a fixed time instant (b) regularly at time intervals commensurate with the rotation time scale. The former is characterized by a dominant inverse energy cascade whereas in the latter, the inverse cascade is stymied due to the recurrent changes in the rotation axis resulting in a strong forward energy transfer and large scale structures that resemble those of isotropic turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Turbulence on a Fractal Fourier set

A novel investigation of the nature of intermittency in incompressible, homogeneous and isotropic turbulence is performed by a numerical study of the Navier-Stokes equations constrained on a fractal Fourier set. The robustness of the energy transfer and of the vortex stretching mechanisms is tested by changing the fractal dimension, D, from the original three dimensional case to a strongly decimated system with D=2.5, where only about $3\%$ of the Fourier modes interact. This is a unique methodology to probe the statistical properties of the turbulent energy cascade, without breaking any of the original symmetries of the equations. While the direct energy cascade persists, deviations from the Kolmogorov scaling are observed in the kinetic energy spectra. A model in terms of a correction with a linear dependency on the co-dimension of the fractal set, $E(k) \sim k^{-5/3 + 3-D}$, explains the results. At small scales, the intermittency of the vorticity field is observed to be quasi-singular as a function of the fractal mode reduction, leading to an almost Gaussian statistics already at $D \sim 2.98$. These effects must be connected to a genuine modification in the triad-to-triad nonlinear energy transfer mechanism.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Clustering of vertically constrained passive particles in homogeneous, isotropic turbulence

We analyze the dynamics of small particles vertically confined, by means of a linear restoring force, to move within a horizontal fluid slab in a three-dimensional (3D) homogeneous isotropic turbulent velocity field. The model that we introduce and study is possibly the simplest description for the dynamics of small aquatic organisms that, due to swimming, active regulation of their buoyancy, or any other mechanism, maintain themselves in a shallow horizontal layer below the free surface of oceans or lakes. By varying the strength of the restoring force, we are able to control the thickness of the fluid slab in which the particles can move. This allows us to analyze the statistical features of the system over a wide range of conditions going from a fully 3D incompressible flow (corresponding to the case of no confinement) to the extremely confined case corresponding to a two-dimensional slice. The background 3D turbulent velocity field is evolved by means of fully resolved direct numerical simulations. Whenever some level of vertical confinement is present, the particle trajectories deviate from that of fluid tracers and the particles experience an effectively compressible velocity field. Here, we have quantified the compressibility, the preferential concentration of the particles, and the correlation dimension by changing the strength of the restoring force. The main result is that there exists a particular value of the force constant, corresponding to a mean slab depth approximately equal to a few times the Kolmogorov length scale, that maximizes the clustering of the particles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Deformation statistics of sub-Kolmogorov-scale ellipsoidal neutrally buoyant drops in isotropic turbulence

Small droplets in turbulent flows can undergo highly variable deformations and orientational dynamics. For neutrally buoyant droplets smaller than the Kolmogorov scale, the dominant effects from the surrounding turbulent flow arise through Lagrangian time histories of the velocity gradient tensor. Here we study the evolution of representative droplets using a model that includes rotation and stretching effects from the surrounding fluid, and restoration effects from surface tension including a constant droplet volume constraint, while assuming that the droplets maintain an ellipsoidal shape. The model is combined with Lagrangian time histories of the velocity gradient tensor extracted from DNS of turbulence to obtain simulated droplet evolutions. These are used to characterize the size, shape and orientation statistics of small droplets in turbulence. A critical capillary number, $Ca_c$ is identified associated with unbounded growth of one or two of the droplet's semi-axes. Exploiting analogies with dynamics of polymers in turbulence, the $Ca_c$ number can be predicted based on the large deviation theory for the largest Finite Time Lyapunov exponent. Also, for sub-critical $Ca$ the theory enables predictions of the slope of the power-law tails of droplet size distributions in turbulence. For cases when the viscosities of droplet and outer fluid differ in a way that enables vorticity to decorrelate the shape from the straining directions, the large deviation formalism based on the stretching properties of the velocity gradient tensor loses validity and its predictions fail. Even considering the limitations of the assumed ellipsoidal droplet shape, the results highlight the complex coupling between droplet deformation, orientation and the local fluid velocity gradient tensor to be expected when small viscous drops interact with turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Disentangling the triadic interactions in Navier-Stokes equations

We study the role of helicity in the dynamics of energy transfer in a modified version of the Navier-Stokes equations with explicit breaking of the mirror symmetry. We select different set of triads participating in the dynamics on the basis of their helicity content. In particular, we remove the negative helically polarized Fourier modes at all wavenumbers except for those falling on a localized shell of wavenumber, $|{\bf k}| \sim k_m$. Changing $k_m$ to be above or below the forcing scale, $k_f$, we are able to assess the energy transfer of triads belonging to different interaction classes. We observe that when the negative helical modes are present only at wavenumber smaller than the forced wavenumbers, an inverse energy cascade develops with an accumulation of energy on a stationary helical condensate. Vice versa, when negative helical modes are present only at wavenumber larger than the forced wavenumbers, a transition from backward to forward energy transfer is observed in the regime when the minority modes become energetic enough.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse energy cascade in nonlocal helical shellmodels of turbulence

Following the exact decomposition in eigenstates of helicity for the Navier-Stokes equations in Fourier space [F. Waleffe, Phys. Fluids A 4, 350 (1992)] we introduce a modified version of helical shell models for turbulence with non-local triadic interactions. By using both an analytical argument and numerical simulation, we show that there exists a class of models, with a specific helical structure, that exhibits a statistically stable inverse energy cascade, in close analogy with that predicted for the Navier-Stokes equations restricted to the same helical interactions. We further support the idea that turbulent energy transfer is the result of a strong entanglement among triads possessing different transfer properties.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the evolution of particle-puffs in turbulence

We study the evolution of turbulent puffs by means of high-resolution numerical simulations. Puffs are bunches of passive particles released from an initially spherical distribution at regular time intervals of the order of the Kolmogorov time. The instantaneous shapes of particle puffs, in particular their asphericity and prolateness, are characterized by measuring the gyration tensor. Analysis has been performed by following, up to one large scale eddy-turn-over time, more than $10^4$ different puffs, each made of 2000 particle tracers, emitted from different places in a homogeneous and isotropic turbulent fluid with Taylor-scale Reynolds number $Re_λ\sim 300$. We also analyze the probability of hitting a given target placed downstream with respect to the local wind at the time of emission, presenting data for three different cases: (i) without any reconstruction of the shape, i.e. considering the bunch of point tracers, and approximating the particle-puff as a (ii) sphere or as an (iii) ellipsoid. The results show a strong dependence on the fluctuations of the instantaneous wind at the moment of the emission and appear to be robust with respect to the approximations (i)-(iii).

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the Peterlin approximation for turbulent flows of polymer solutions

We study the impact of the Peterlin approximation on the statistics of the end-to-end separation of poly- mers in a turbulent flow. The FENE and FENE-P models are numerically integrated along a large number of Lagrangian trajectories resulting from a direct numerical simulation of three-dimensional homogeneous isotropic turbulence. Although the FENE-P model yields results in qualitative agreement with those of the FENE model, quantitative differences emerge. The steady-state probability of large extensions is overesti- mated by the FENE-P model. The alignment of polymers with the eigenvectors of the rate-of-strain tensor and with the direction of vorticity is weaker when the Peterlin approximation is used. At large Weissenberg numbers, both the correlation times of the extension and of the orientation of polymers are underestimated by the FENE-P model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The importance of chemical potential in the determination of water slip in nanochannels

We investigate the slip properties of water confined in graphite-like nano-channels by non-equilibrium molecular dynamics simulations, with the aim of identifying and analyze separately the influence of different physical quantities on the slip length. In a system under confinement but connected to a reservoir of fluid, the chemical potential is the natural control parameter: we show that two nanochannels characterized by the same macroscopic contact angle -- but a different microscopic surface potential -- do not exhibit the same slip length unless the chemical potential of water in the two channels is matched. Some methodological issues related to the preparation of samples for the comparative analysis in confined geometries are also discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new assessment of the second order moment of Lagrangian velocity increments in turbulence

The behavior of the second-order Lagrangian structure functions on state-of-the-art numerical data both in two and three dimensions is studied. On the basis of a phenomenological connection between Eulerian space-fluctuations and the Lagrangian time-fluctuations, it is possible to rephrase the Kolmogorov $4/5$-law into a relation predicting the linear (in time) scaling for the second order Lagrangian structure function. When such a function is directly observed on current experimental or numerical data, it does not clearly display a scaling regime. A parameterization of the Lagrangian structure functions based on Batchelor model is introduced and tested on data for $3d$ turbulence, and for $2d$ turbulence in the inverse cascade regime. Such parameterization supports the idea, previously suggested, that both Eulerian and Lagrangian data are consistent with a linear scaling plus finite-Reynolds number effects affecting the small- and large-time scales. When large-time saturation effects are properly accounted for, compensated plots show a detectable plateau already at the available Reynolds number. Furthermore, this parameterization allows us to make quantitative predictions on the Reynolds number value for which Lagrangian structure functions are expected to display a scaling region. Finally, we show that this is also sufficient to predict the anomalous dependency of the normalized root mean squared acceleration as a function of the Reynolds number, without fitting parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Regularization of the slip length divergence in water nanoflows by inhomogeneities at the Angstrom scale

We performed non-equilibrium Molecular Dynamics simulations of water flow in nano-channels with the aim of discriminating {\it static} from {\it dynamic} contributions of the solid surface to the slip length of the molecular flow. We show that the regularization of the slip length divergence at high shear rates, formerly attributed to the wall dynamics, is controlled instead by its static properties. Surprisingly, we find that atomic displacements at the Angstrom scale are sufficient to remove the divergence of the slip length and realize the no-slip condition. Since surface thermal fluctuations at room temperature are enough to generate these displacements, we argue that the no-slip condition for water can be achieved also for ideal surfaces, which do not present any surface roughness.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lagrangian single particle turbulent statistics through the Hilbert-Huang Transform

The Hilbert-Huang transform is applied to analyze single particle Lagrangian velocity data from numerical simulations of hydrodynamic turbulence. The velocity trajectory is described in terms of a set of intrinsic mode functions, C_{i}(t), and of their instantaneous frequency, ω_{i}(t). On the basis of this decomposition we define the ω-conditioned statistical moments of the C_{i} modes, named q-order Hilbert Spectra (HS). We show that such new quantities have enhanced scaling properties as compared to traditional Fourier transform- or correlation-based (Structure Functions) statistical indicators, thus providing better insights into the turbulent energy transfer process. We present a clear empirical evidence that the energy-like quantity, i.e. the second-order HS, displays a linear scaling in time in the inertial range, as expected from dimensional analysis and never observed before. We also measure high order moment scaling exponents in a direct way, without resorting the Extended Self Similarity (ESS) procedure. This leads to a new estimate of the Lagrangian structure functions exponents which are consistent with the multifractal prediction in the Lagrangian frame as proposed in [Biferale et al., Phys. Rev. Lett. vol. 93, 064502 (2004)].

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lattice Monte Carlo methods for systems far from equilibrium

We present a new numerical Monte Carlo approach to determine the scaling behavior of lattice field theories far from equilibrium. The presented methods are generally applicable to systems where classical-statistical fluctuations dominate the dynamics. As an example, these methods are applied to the random-force-driven one-dimensional Burgers' equation - a model for hydrodynamic turbulence. For a self-similar forcing acting on all scales the system is driven to a nonequilibrium steady state characterized by a Kolmogorov energy spectrum. We extract correlation functions of single- and multi-point quantities and determine their scaling spectrum displaying anomalous scaling for high-order moments. Varying the external forcing we are able to tune the system continuously from equilibrium, where the fluctuations are short-range correlated, to the case where the system is strongly driven in the infrared. In the latter case the nonequilibrium scaling of small-scale fluctuations are shown to be universal.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On the Global Regularity of a Helical-decimated Version of the 3D Navier-Stokes Equations

We study the global regularity, for all time and all initial data in $H^{1/2}$, of a recently introduced decimated version of the incompressible 3D Navier-Stokes (dNS) equations. The model is based on a projection of the dynamical evolution of Navier-Stokes (NS) equations into the subspace where helicity (the $L^2-$scalar product of velocity and vorticity) is sign-definite. The presence of a second (beside energy) sign-definite inviscid conserved quadratic quantity, which is equivalent to the $H^{1/2}-$Sobolev norm, allows us to demonstrate global existence and uniqueness, of space-periodic solutions, together with continuity with respect to the initial conditions, for this decimated 3D model. This is achieved thanks to the establishment of two new estimates, for this 3D model, which show that the $H^{1/2}$ and the time average of the square of the $H^{3/2}$ norms of the velocity field remain finite. Such two additional bounds are known, in the spirit of the work of H. Fujita and T. Kato \cite{kato1,kato2}, to be sufficient for showing well-posedness for the 3D NS equations. Furthermore, they are directly linked to the helicity evolution for the dNS model, and therefore with a clear physical meaning and consequences.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stick-Slip Sliding of Water Drops on Chemically Heterogeneous Surfaces

We present a comprehensive study of water drops sliding down chemically heterogeneous surfaces formed by a periodic pattern of alternating hydrophobic and hydrophilic stripes. Drops are found to undergo a stick-slip motion whose average speed is an order of magnitude smaller than that measured on a homogeneous surface having the same static contact angle. This motion is the result of the periodic deformations of the drop interface when crossing the stripes. Numerical simulations confirm this view and are used to elucidate the principles underlying the experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Breakup of small aggregates driven by turbulent hydrodynamic stress

Breakup of small solid aggregates in homogeneous and isotropic turbulence is studied theoretically and by using Direct Numerical Simulations at high Reynolds number, Re_λ \simeq 400. We show that turbulent fluctuations of the hydrodynamic stress along the aggregate trajectory play a key role in determining the aggregate mass distribution function. Differences between turbulent and laminar flows are discussed. A novel definition of the fragmentation rate is proposed in terms of the typical frequency at which the hydrodynamic stress becomes sufficiently high to cause breakup along each Lagrangian path. We also define an Eulerian proxy of the real fragmentation rate, based on the joint statistics of the stress and its time derivative, which should be easier to measure in any experimental set-up. Both our Eulerian and Lagrangian formulations define a clear procedure for the computation of the mass distribution function due to fragmentation. Contrary, previous estimates based only on single point statistics of the hydrodynamic stress exhibit some deficiencies. These are discussed by investigating the evolution of an ensemble of aggregates undergoing breakup and aggregation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Droplet size distribution in homogeneous isotropic turbulence

We study the physics of droplet breakup in a statistically stationary homogeneous and isotropic turbulent flow by means of high resolution numerical investigations based on the multicomponent lattice Boltzmann method. We verified the validity of the criterion proposed by Hinze (1955) for droplet breakup and we measured the full probability distribution function (pdf) of droplets radii at different Reynolds numbers and for different volume fraction. By means of a Lagrangian tracking we could follow individual droplets along their trajectories, define a local Weber number based on the velocity gradients and study its cross-correlation with droplet deformation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Inverse energy cascade in three-dimensional isotropic turbulence

In turbulent flows kinetic energy is spread by nonlinear interactions over a broad range of scales. Energy transfer may proceed either toward small scales or in the reverse direction. The latter case is peculiar of two-dimensional (2D) flows. Interestingly, a reversal of the energy flux is observed also in three-dimensional (3D) geophysical flows under rotation and/or confined in thin layers. The question is whether this phenomenon is enforced solely by external anisotropic mechanisms or it is intimately embedded in the Navier-Stokes (NS) equations. Here we show that an inverse energy cascade occurs also in 3D isotropic flow. The flow is obtained from a suitable surgery of the NS equations, keeping only triadic interactions between sign-defined helical modes, preserving homogeneity and isotropy and breaking reflection invariance. Our findings highlight the role played by helicity in the energy transfer process and show that both 2D and 3D properties naturally coexist in all flows in nature.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Multi-Time Multi-Scale Correlation Functions in Hydrodynamic Turbulence

High Reynolds numbers Navier-Stokes equations are believed to break self-similarity concerning both spatial and temporal properties: correlation functions of different orders exhibit distinct decorrelation times and anomalous spatial scaling properties. Here, we present a systematic attempt to measure multi-time and multi-scale correlations functions, by using high Reynolds numbers numerical simulations of fully homogeneous and isotropic turbulent flow. The main idea is to set-up an ensemble of probing stations riding the flow, i.e. measuring correlations in a reference frame centered on the trajectory of distinct fluid particles (the quasi-Lagrangian reference frame introduced by Belinicher & L'vov, Sov. Phys. JETP 66, 303 (1987)). In this way we reduce the large-scale sweeping and measure the non-trivial temporal dynamics governing the turbulent energy transfer from large to small scales. We present evidences of the existence of dynamic multiscaling: multi-time correlation functions are characterized by an infinite set of characteristic times.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Droplet breakup in homogeneous and isotropic turbulence

This fluid dynamics video shows the breakup of a droplet in a stationary homogeneous and isotropic turbulent flow. We consider droplets with the same density of the transporting fluid. The droplets and the fluid are numerically modelled by means of a multicompo- nent Lattice-Boltzmann method. The turbulent fluid is maintained through a large scale stirring force and the radius of stable droplets, for the parameters in our simulation, is larger than the Kolmogorov scale. Events of droplet deformation, break-up and aggregation are clearly visible from the movie. With the present database droplet evo- lution can be studied from both an Eulerian and Lagrangian point of view. The Kolmogorov-Hinze criteria for droplets break-up can be tested also by means of simulations with different viscosity contrast between the two components.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Lattice Boltzmann Methods for thermal flows: continuum limit and applications to compressible Rayleigh-Taylor systems

We compute the continuum thermo-hydrodynamical limit of a new formulation of lattice kinetic equations for thermal compressible flows, recently proposed in [Sbragaglia et al., J. Fluid Mech. 628 299 (2009)]. We show that the hydrodynamical manifold is given by the correct compressible Fourier- Navier-Stokes equations for a perfect fluid. We validate the numerical algorithm by means of exact results for transition to convection in Rayleigh-Bénard compressible systems and against direct comparison with finite-difference schemes. The method is stable and reliable up to temperature jumps between top and bottom walls of the order of 50% the averaged bulk temperature. We use this method to study Rayleigh-Taylor instability for compressible stratified flows and we determine the growth of the mixing layer at changing Atwood numbers up to At ~ 0.4. We highlight the role played by the adiabatic gradient in stopping the mixing layer growth in presence of high stratification and we quantify the asymmetric growth rate for spikes and bubbles for two dimensional Rayleigh- Taylor systems with resolution up to Lx \times Lz = 1664 \times 4400 and with Rayleigh numbers up to Ra ~ 2 \times 10^10.

0 citations0 reviewsNo code linkedOpen access