Source author record

L. Biferale

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

26 published item(s)

preprint2020arXiv

TURB-Rot. A large database of 3d and 2d snapshots from turbulent rotating flows

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 http://smart-turb.roma2.infn.it

preprint2019arXiv

Synchronization to big-data: nudging the Navier-Stokes equations for data assimilation of turbulent flows

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.

preprint2016arXiv

Fluctuating Multicomponent Lattice Boltzmann Model

Current implementations of fluctuating lattice Boltzmann equations (FLBE) describe single component fluids. In this paper, a model based on the continuum kinetic Boltzmann equation for describing multicomponent fluids is extended to incorporate the effects of thermal fluctuations. The thus obtained fluctuating Boltzmann equation is first linearized to apply the theory of linear fluctuations, and expressions for the noise covariances are determined by invoking the fluctuation-dissipation theorem (FDT) directly at the kinetic level. Crucial for our analysis is the projection of the Boltzmann equation onto the ortho-normal Hermite basis. By integrating in space and time the fluctuating Boltzmann equation with a discrete number of velocities, the FLBE is obtained for both ideal and non-ideal multicomponent fluids. Numerical simulations are specialized to the case where mean-field interactions are introduced on the lattice, indicating a proper thermalization of the system.

preprint2016arXiv

Inertial-particle accelerations in turbulence: a Lagrangian closure

The distribution of particle accelerations in turbulence is intermittent, with non-Gaussian tails that are quite different for light and heavy particles. In this article we analyse a closure scheme for the acceleration fluctuations of light and heavy inertial particles in turbulence, formulated in terms of Lagrangian correlation functions of fluid tracers. We compute the variance and the flatness of inertial particle accelerations and we discuss their dependency on the Stokes number. The closure incorporates effects induced by the Lagrangian correlations along the trajectories of fluid tracers, and its predictions agree well with results of direct numerical simulations of inertial particles in turbulence, provided that the effects induced by the inertial preferential sampling of heavy/light particles outside/inside vortices are negligible. In particular, the scheme predicts the correct functional behaviour of the acceleration variance, as a function of Stokes, as well as the presence of a minimum/maximum for the flatness of the acceleration of heavy/light particles, in good qualitative agreement with numerical data. We also show that the closure works well when applied to the Lagrangian evolution of particles using a stochastic surrogate for the underlying Eulerian velocity field. Our results support the conclusion that there exist important contributions to the statistics of the acceleration of inertial particles independent of the preferential sampling. For heavy particles we observe deviations between the predictions of the closure scheme and direct numerical simulations, at Stokes numbers of order unity. For light particles the deviation occurs for larger Stokes numbers.

preprint2016arXiv

Sliding drops across alternating hydrophobic and hydrophilic stripes

We perform a joint numerical and experimental study to sistematically characterize the motion of drops sliding over a periodic array of alternating hydrophobic and hydrophilic stripes with large wettability contrast, and typical width of hundreds of $μ\textrm{m}$. The fraction of the hydrophobic stripes has been varied from about 20% to 80%. The effects of the heterogeneous patterning can be described by a renormalized value of the critical Bond number, i.e. the critical dimensionless force needed to depin the drop before it starts to move. Close to the critical Bond number we observe a jerkily motion characterized by an evident stick-slip dynamics. As a result, dissipation is strongly localized in time, and the mean velocity of the drops can easily decrease by an order of magnitude compared to the sliding on homogeneous surface. Lattice Boltzmann (LB) numerical simulations are crucial for disclosing to what extent the sliding dynamics can be deduced from the computed balance of capillary, viscous and body forces at varying the Bond number, the surface composition and the liquid viscosity. Away from the critical Bond number, we characterize both experimentally and numerically the dissipation inside the droplet by studying the relation between the average velocity and the applied volume forces.

preprint2015arXiv

Evolution of a double-front Rayleigh-Taylor system using a GPU-based high resolution thermal Lattice-Boltzmann model

We study the turbulent evolution originated from a system subjected to a Rayleigh-Taylor instability with a double density at high resolution in a 2 dimensional geometry using a highly optimized thermal Lattice Boltzmann code for GPUs. The novelty of our investigation stems from the initial condition, given by the superposition of three layers with three different densities, leading to the development of two Rayleigh-Taylor fronts that expand upward and downward and collide in the middle of the cell. By using high resolution numerical data we highlight the effects induced by the collision of the two turbulent fronts in the long time asymptotic regime. We also provide details on the optimized Lattice-Boltzmann code that we have run on a cluster of GPUs

preprint2015arXiv

Homogeneous and Isotropic Turbulence: a short survey on recent developments

We present a detailed review of some of the most recent developments on Eulerian and Lagrangian turbulence in homogeneous and isotropic statistics. In particular, we review phenomenological and numerical results concerning the issue of universality with respect to the large scale forcing and the viscous dissipative physics. We discuss the state-of-the-art of numerical versus experimental comparisons and we discuss the dicotomy between phenomenology based on coherent structures or on statistical approaches. A detailed discussion of finite Reynolds effects is also presented.

preprint2015arXiv

Pair separation of magnetic elements in the quiet Sun

The dynamic properties of the quiet Sun photosphere can be investigated by analyzing the pair dispersion of small-scale magnetic fields (i.e., magnetic elements). By using $25$ hr-long Hinode magnetograms at high spatial resolution ($0".3$), we tracked $68,490$ magnetic element pairs within a supergranular cell near the disk center. The computed pair separation spectrum, calculated on the whole set of particle pairs independently of their initial separation, points out what is known as a super-diffusive regime with spectral index $γ=1.55\pm0.05$, in agreement with the most recent literature, but extended to unprecedented spatial and temporal scales (from granular to supergranular). Furthermore, for the first time, we investigated here the spectrum of the mean square displacement of pairs of magnetic elements, depending on their initial separation $r_0$. We found that there is a typical initial distance above (below) which the pair separation is faster (slower) than the average. A possible physical interpretation of such a typical spatial scale is also provided.

preprint2015arXiv

Statistical analysis of the velocity and scalar fields in reacting turbulent wall-jets

The concept of local isotropy in a chemically reacting turbulent wall-jet flow is addressed using direct numerical simulation (DNS) data. Different DNS databases with isothermal and exothermic reactions are examined. The chemical reaction and heat release effects on the turbulent velocity, passive scalar and reactive species fields are studied using their probability density functions (PDF) and higher order moments for velocities and scalar fields, as well as their gradients. With the aid of the anisotropy invariant maps for the Reynolds stress tensor the heat release effects on the anisotropy level at different wall-normal locations are evaluated and found to be most accentuated in the near-wall region. It is observed that the small-scale anisotropies are persistent both in the near-wall region and inside the jet flame. Two exothermic cases with different Damkohler number are examined and the comparison revealed that the Damkohler number effects are most dominant in the near-wall region, where the wall cooling effects are influential. In addition, with the aid of PDFs conditioned on the mixture fraction, the significance of the reactive scalar characteristics in the reaction zone is illustrated. We argue that the combined effects of strong intermittency and strong persistency of anisotropy at the small scales in the entire domain can affect mixing and ultimately the combustion characteristics of the reacting flow.

preprint2014arXiv

Natural convection with mixed insulating and conducting boundary conditions: low and high Rayleigh numbers regimes

We investigate the stability and dynamics of natural convection in two dimensions, subject to inhomogeneous boundary conditions. In particular, we consider a Rayleigh-Bènard (RB) cell, where the horizontal top boundary contains a periodic sequence of alternating thermal insulating and conducting patches, and we study the effects of the heterogeneous pattern on the global heat exchange, both at low and high Rayleigh numbers. At low Rayleigh numbers, we determine numerically the transition from a regime characterized by the presence of small convective cells localized at the inhomogeneous boundary to the onset of bulk convective rolls spanning the entire domain. Such a transition is also controlled analytically in the limit when the boundary pattern length is small compared with the cell vertical size. At higher Rayleigh number, we use numerical simulations based on a lattice Boltzmann method to assess the impact of boundary inhomogeneities on the fully turbulent regime up to $Ra \sim 10^{10}$.

preprint2013arXiv

Hydrodynamics of air entrainment by moving contact lines

We study the dynamics of the interface between two immiscible fluids in contact with a chemically homogeneous moving solid plate. We consider the generic case of two fluids with any viscosity ratio and of a plate moving in either directions (pulled or pushed in the bath). The problem is studied by a combination of two models, namely an extension to finite viscosity ratio of the lubrication theory and a Lattice Boltzmann method. Both methods allow to resolve, in different ways, the viscous singularity at the triple contact between the two fluids and the wall. We find a good agreement between the two models particularly for small capillary numbers. When the solid plate moves fast enough, the entrainment of one fluid into the other one can occur. The extension of the lubrication model to the case of a non-zero air viscosity, as developed here, allows us to study the dependence of the critical capillary number for air entrainment on the other parameters in the problem (contact angle and viscosity ratio).

preprint2012arXiv

Extreme events in the dispersions of two neighboring particles under the influence of fluid turbulence

We present a numerical study of two-particle dispersion from point-sources in 3D incompressible Homogeneous and Isotropic turbulence, at Reynolds number Re \simeq 300. Tracer particles are emitted in bunches from localized sources smaller than the Kolmogorov scale. We report the first quantitative evidence, supported by an unprecedented statistics, of the deviations of relative dispersion from Richardson's picture. Deviations are due to extreme events of pairs separating much faster than average, and of pairs remaining close for long times. The two classes of events are the fingerprint of complete different physics, the former being dominated by inertial subrange and large-scale fluctuations, while the latter by the dissipation subrange. A comparison of relative separation in surrogate white-in-time velocity field, with correct viscous-, inertial- and integral-scale properties allows us to assess the importance of temporal correlations along tracer trajectories.

preprint2011arXiv

Convection in multiphase flows using Lattice Boltzmann methods

We present high resolution numerical simulations of convection in multiphase flows (boiling) using a novel algorithm based on a Lattice Boltzmann method. We first validate the thermodynamical and kinematical properties of the algorithm. Then, we perform a series of 3d numerical simulations at changing the mean properties in the phase diagram and compare convection with and without phase coexistence at $Ra \sim 10^7$. We show that in presence of nucleating bubbles non-Oberbeck Boussinesq effects develops, mean temperature profile becomes asymmetric, heat-transfer and heat-transfer fluctuations are enhanced. We also show that small-scale properties of velocity and temperature fields are strongly affected by the presence of buoyant bubble leading to high non-Gaussian profiles in the bulk.

preprint2011arXiv

Second order closure for stratified convection: bulk region and overshooting

The parameterization of small-scale turbulent fluctuations in convective systems and in the presence of strong stratification is a key issue for many applied problems in oceanography, atmospheric science and planetology. In the presence of stratification, one needs to cope with bulk turbulent fluctuations and with inversion regions, where temperature, density -or both- develop highly non-linear mean profiles due to the interactions between the turbulent boundary layer and the unmixed -stable- flow above/below it. We present a second order closure able to cope simultaneously with both bulk and boundary layer regions, and we test it against high-resolution state-of-the-art 2D numerical simulations in a convective and stratified belt for values of the Rayleigh number, up to Ra = 10^9. Data are taken from a Rayleigh-Taylor system confined by the existence of an adiabatic gradient.

preprint2010arXiv

High resolution numerical study of Rayleigh-Taylor turbulence using a thermal lattice Boltzmann scheme

We present results of a high resolution numerical study of two dimensional (2d) Rayleigh-Taylor turbulence using a recently proposed thermal lattice Boltzmann method (LBT). The goal of our study is both methodological and physical. We assess merits and limitations concerning small- and large-scale resolution/accuracy of the adopted integration scheme. We discuss quantitatively the requirements needed to keep the method stable and precise enough to simulate stratified and unstratified flows driven by thermal active fluctuations at high Rayleigh and high Reynolds numbers. We present data with spatial resolution up to 4096 x 10000 grid points and Rayleigh number up to Ra ~ 10^11 . The statistical quality of the data allows us to investigate velocity and temperature fluctuations, scale-by-scale, over roughly four decades. We present a detailed quantitative analysis of scaling laws in the viscous, inertial and integral range, supporting the existence of a Bolgiano-like inertial scaling, as expected in 2d systems. We also discuss the presence of small/large intermittent deviation to the scaling of velocity/temperature fluctuations and the Rayleigh dependency of gradients flatness.

preprint2010arXiv

Numerical simulations of compressible Rayleigh-Taylor turbulence in stratified fluids

We present results from numerical simulations of Rayleigh-Taylor turbulence, performed using a recently proposed lattice Boltzmann method able to describe consistently a thermal compressible flow subject to an external forcing. The method allowed us to study the system both in the nearly-Boussinesq and strongly compressible regimes. Moreover, we show that when the stratification is important, the presence of the adiabatic gradient causes the arrest of the mixing process.

preprint2009arXiv

Fully developed turbulence and the multifractal conjecture

We review the Parisi-Frisch MultiFractal formalism for Navier--Stokes turbulence with particular emphasis on the issue of statistical fluctuations of the dissipative scale. We do it for both Eulerian and Lagrangian Turbulence. We also show new results concerning the application of the formalism to the case of Shell Models for turbulence. The latter case will allow us to discuss the issue of Reynolds number dependence and the role played by vorticity and vortex filaments in real turbulent flows.

preprint2009arXiv

Hydrokinetic simulations of nanoscopic precursor films in rough channels

We report on simulations of capillary filling of high-wetting fluids in nano-channels with and without obstacles. We use atomistic (molecular dynamics) and hydrokinetic (lattice-Boltzmann) approaches which point out clear evidence of the formation of thin precursor films, moving ahead of the main capillary front. The dynamics of the precursor films is found to obey a square-root law as the main capillary front, z^2(t) ~ t, although with a larger prefactor, which we find to take the same value for the different geometries (2D-3D) under inspection. The two methods show a quantitative agreement which indicates that the formation and propagation of thin precursors can be handled at a mesoscopic/hydrokinetic level. This can be considered as a validation of the Lattice-Boltzmann (LB) method and opens the possibility of using hydrokinetic methods to explore space-time scales and complex geometries of direct experimental relevance. Then, LB approach is used to study the fluid behaviour in a nano-channel when the precursor film encounters a square obstacle. A complete parametric analysis is performed which suggests that thin-film precursors may have an important influence on the efficiency of nanochannel-coating strategies.

preprint2009arXiv

Lattice Boltzmann method with self-consistent thermo-hydrodynamic equilibria

Lattice kinetic equations incorporating the effects of external/internal force fields via a shift of the local fields in the local equilibria, are placed within the framework of continuum kinetic theory. The mathematical treatment reveals that, in order to be consistent with the correct thermo-hydrodynamical description, temperature must also be shifted, besides momentum. New perspectives for the formulation of thermo-hydrodynamic lattice kinetic models of non-ideal fluids are then envisaged. It is also shown that on the lattice, the definition of the macroscopic temperature requires the inclusion of new terms directly related to discrete effects. The theoretical treatment is tested against a controlled case with a non ideal equation of state.

preprint2009arXiv

Turbulent pair dispersion of inertial particles

The relative dispersion of pairs of inertial particles in incompressible, homogeneous, and isotropic turbulence is studied by means of direct numerical simulations at two values of the Taylor-scale Reynolds number $Re_λ \sim 200$ and 400. The evolution of both heavy and light particle pairs is analysed at varying the particle Stokes number and the fluid-to-particle density ratio. For heavy particles, it is found that turbulent dispersion is schematically governed by two temporal regimes. The first is dominated by the presence, at large Stokes numbers, of small-scale caustics in the particle velocity statistics, and it lasts until heavy particle velocities have relaxed towards the underlying flow velocities. At such large scales, a second regime starts where heavy particles separate as tracers particles would do. As a consequence, at increasing inertia, a larger transient stage is observed, and the Richardson diffusion of simple tracers is recovered only at large times and large scales. These features also arise from a statistical closure of the equation of motion for heavy particle separation that is proposed, and which is supported by the numerical results. In the case of light particles with high density ratios, strong small-scale clustering leads to a considerable fraction of pairs that do not separate at all, although the mean separation increases with time. This effect strongly alters the shape of the probability density function of light particle separations.

preprint2008arXiv

Universal intermittent properties of particle trajectories in highly turbulent flows

We present a collection of eight data sets, from state-of-the-art experiments and numerical simulations on turbulent velocity statistics along particle trajectories obtained in different flows with Reynolds numbers in the range $R_λ\in [120:740]$. Lagrangian structure functions from all data sets are found to collapse onto each other on a wide range of time lags, revealing a universal statistics, and calling for a unified theoretical description. Parisi-Frisch Multifractal theory, suitable extended to the dissipative scales and to the Lagrangian domain, is found to capture intermittency of velocity statistics over the whole three decades of temporal scales here investigated.

preprint2005arXiv

Mesoscopic two-phase model for describing apparent slip in micro-channel flows

The phenomenon of apparent slip in micro-channel flows is analyzed by means of a two-phase mesoscopic lattice Boltzmann model including non-ideal fluid-fluid and fluid-wall interactins. The weakly-inhomogeneous limit of this model is solved analytically. The present mesoscopic approach permits to access much larger scales than molecular dynamics, and comparable with those attained by continuum methods. However, at variance with the continuum approach, the existence of a gas layer near the wall does not need to be postulated a priori, but emerges naturally from the underlying non-ideal mesoscopic dynamics. It is therefore argued that a mesoscopic Lattice Boltzmann approach with non-ideal fluid-fluid and fluid-wall interactions might achieve an optimal compromise between physical realism and computational efficiency for the study of channel micro-flows.

preprint2005arXiv

Multi-particle dispersion in fully developed turbulence

The statistical geometry of dispersing Lagrangian clusters of four particles (tetrahedra) is studied by means of high-resolution direct numerical simulations of three-dimensional homogeneous isotropic turbulence. We give the first evidence of a self-similar regime of shape dynamics characterized by almost two-dimensional, strongly elongated geometries. The analysis of four-point velocity-difference statistics and orientation shows that inertial-range eddies typically generate a straining field with a strong extensional component aligned with the elongation direction and weak extensional/compressional components in the orthogonal plane.

preprint1995arXiv

Generalized scaling in fully developed turbulence

In this paper we report numerical and experimental results on the scaling properties of the velocity turbulent fields in several flows. The limits of a new form of scaling, named Extended Self Similarity(ESS), are discussed. We show that, when a mean shear is absent, the self scaling exponents are universal and they do not depend on the specific flow (3D homogeneous turbulence, thermal convection , MHD). In contrast, ESS is not observed when a strong shear is present. We propose a generalized version of self scaling which extends down to the smallest resolvable scales even in cases where ESS is not present. This new scaling is checked in several laboratory and numerical experiment. A possible theoretical interpretation is also proposed. A synthetic turbulent signal having most of the properties of a real one has been generated.

preprint1994arXiv

Transition to Chaos in a Shell Model of Turbulence

We study a shell model for the energy cascade in three dimensional turbulence at varying the coefficients of the non-linear terms in such a way that the fundamental symmetries of Navier-Stokes are conserved. When a control parameter $ε$ related to the strength of backward energy transfer is enough small, the dynamical system has a stable fixed point corresponding to the Kolmogorov scaling. This point becomes unstable at $ε=0.3843...$ where a stable limit cycle appears via a Hopf bifurcation. By using the bi-orthogonal decomposition, the transition to chaos is shown to follow the Ruelle-Takens scenario. For $ε> 0.3953..$ the dynamical evolution is intermittent with a positive Lyapunov exponent. In this regime, there exists a strange attractor which remains close to the Kolmogorov (now unstable) fixed point, and a local scaling invariance which can be described via a intermittent one-dimensional map.