Researcher profile

Luca Biferale

Luca Biferale contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.