BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Enrico Calzavarini

Enrico Calzavarini contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dyncond-mat.softnlin.CDphysics.ao-phBiological Physicsphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
6topics
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

17 published item(s)

preprint2025arXiv

Freezing and ice aging dynamics in saline water under natural convection

Understanding the coupled dynamics of liquid-solid phase change and fluid flows is crucial in a wide range of geophysical and industrial applications. When freezing occurs in saline water, the newly formed ice is mushy, with a porous structure that traps the brine within the ice. In this work, which combines experiments and theoretical analyses, we investigate the long-term evolution of saline ice, comprehensively accounting for the coupled dynamics of multiscale fluid flow, heat and mass transfer, and phase change. We show that in a closed convective system the rapid formation of a mushy ice layer is followed by desalination (i.e, the expulsion of salt from the ice) processes that might lead to a slow asymptotic decrease of the ice thickness. Desalination of mushy ice reduces its porosity, which alters the dynamic thermal equilibrium and ice thickness by weakening buoyancy-driven convection within the mushy layer. In turn, changes in brine convection and ice thickness affect the further desalination of the ice. The long-term dynamics of the system can be accurately predicted by a one-dimensional model based on appropriate parameterizations of global heat and mass transfer properties. Furthermore, within the same theoretical model we explore the ice dynamics across a broader parameter space. Our findings advance the understanding of the coupled phase-change physics of saline solutions in the presence of convective fluid flows and provide a basis for explaining and predicting real-world phenomena such as the aging of sea ice.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Dynamics of finite-size spheroids in turbulent flow: the roles of flow structures and particle boundary layers

We study the translational and rotational dynamics of neutrally-buoyant finite-size spheroids in hydrodynamic turbulence by means of fully resolved numerical simulations. We examine axisymmetric shapes, from oblate to prolate, and the particle volume dependences. We show that the accelerations and rotations experienced by non-spherical inertial-scale particles result from volume filtered fluid forces and torques, similar to spherical particles. However, the particle orientations carry signatures of preferential alignments with the surrounding flow structures, which is reflected in distinct axial and lateral fluctuations for accelerations and rotation rates. The randomization of orientations does not occur even for particles with volume equivalent diameter size in the inertial range, here up to 60 $η$ at $Re_λ=120$. Additionally, we demonstrate that the role of fluid boundary layers around the particles cannot be neglected to reach a quantitative understanding of particle statistical dynamics, as they affect the intensities of angular velocities, and the relative importance of tumbling with respect to spinning rotations. This study brings to the fore the importance of inertial-scale flow structures in homogeneous and isotropic turbulence and their impacts on the transport of neutrally-buoyant bodies with size in the inertial range.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Residence time of inertial particles in 3D thermal convection: implications for magma reservoirs

The dynamic behavior of crystals in convecting fluids determines how magma bodies solidify. In particular, it is often important to estimate how long crystals stay in suspension in the host liquid before being deposited at its bottom (or top, for light particles). We perform a systematic 3D numerical study of particle-laden Rayleigh-Benard convection, and derive a robust model for the particle residence time. For Rayleigh numbers higher than 10^7, inertial particles&#39; trajectories exhibit a monotonic transition from fluid tracer-like to free-fall dynamics, the control parameter being the ratio between particle Stokes velocity and the fluid velocity. The average settling rate is proportional to the particle Stokes velocity in both the end-member regimes, but the distribution of the residence times differs markedly from one to the other. For lower Rayleigh numbers (<10^7), an interaction between large-scale circulation and particle motion emerges, increasing the settling rates on average. Nevertheless, the mean residence time does not exceed the terminal time, i.e. the settling time from a quiescent fluid, by a factor larger than four. An exception are simulations with only a slightly super-critical Rayleigh number (~ 10^4), for which stationary convection develops and some particles become trapped indefinitely. 2D simulations of the same problem overestimate the flow-particle interaction - and hence the residence time - for both high and low Rayleigh numbers, which stresses the importance of using 3D geometries for simulating particle-laden flows. We outline how our model can be used to explain depth changes of crystal size distribution in sedimentary layers of magmatic intrusions that are thought to have formed via settling of a crystal cargo, and discuss how the micro-structural observations of solidified intrusions can be used to infer the past convective velocity of magma.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Surface light modulation by sea ice and phytoplankton survival in a convective flow model

Plankton dynamics depend in a complex manner on a variety of physical phenomena, according to both experimental and numerical data. In particular, experimental field studies have highlighted the relation between phytoplankton survival and turbulent upwelling and downwelling from thermal convection. Recent numerical works have also shown the importance of accounting for advective transport by persistent structures in simulation models. In nutrient-rich polar marine environments phytoplankton blooms are critically limited by light availability under ice-covered waters. Such heterogeneity of the light intensity distribution, in association with a large-scale coherent fluid flow, can give rise to nontrivial growth dynamics. In this work we extend a previous advection-reaction-diffusion model of phytoplankton light-limited vertical dynamics in the presence of convective transport. Specifically, we consider horizontally heterogeneous light conditions through the use of two regions with different production regimes, modelling the absence (presence) of light under (in between) obstacles. Such a model is intended as an idealized representation of nonuniformly ice-covered polar waters. By means of numerical simulations, we find that the main role of advective transport is to hinder phytoplankton growth, but also that such effect depends on the positions of the obstacles with respect to the upwelling and downwelling flow regions. Furthermore, we show that the sinking speed due to the density difference between phytoplankton organisms and water, while small, plays an important role, which depends on how it adds to the flow. These results indicate that advective transport can have a crucial impact on the survival conditions of sinking phytoplankton species in polar environments.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Three-dimensional turbulence effects on plankton dynamics behind an obstacle

We study a predator-prey model of plankton dynamics in the two and three-dimensional wakes of turbulent flows behind a cylinder, focusing on the impact of the three-dimensional character of the carrying velocity field on population variance spectra and spatial distributions. By means of direct numerical simulations, we find that the qualitative behavior of the biological dynamics is mostly independent of the space dimensionality, which suggests that only the relation between the typical flow and biological timescales is crucial to observe persistent blooms. Similarly, in both cases, we find that the spectral properties of the planktonic populations are essentially indistinguishable from those of an inert tracer. The main difference arising from the comparison of the two and three-dimensional configurations concerns the local spatial distribution of plankton density fields. Indeed, the three-dimensional turbulent dynamics tend to destroy the localized coherent structures characterizing the two-dimensional flow, in which the planktonic species are mostly concentrated, thus reducing the phytoplankton biomass in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Anisotropic particles in two-dimensional convective turbulence

The orientational dynamics of inertialess anisotropic particles transported by two-dimensional convective turbulent flows display a coexistence of regular and chaotic features. We numerically demonstrate that very elongated particles (rods) align preferentially with the direction of the fluid flow, i.e., horizontally close to the isothermal walls and dominantly vertically in the bulk. This behaviour is due to the the presence of a persistent large scale circulation flow structure, which induces strong shear at wall boundaries and in up/down-welling regions. The near-wall horizontal alignment of rods persists at increasing the Rayleigh number, while the vertical orientation in the bulk is progressively weakened by the corresponding increase of turbulence intensity. Furthermore, we show that very elongated particles are nearly orthogonal to the orientation of the temperature gradient, an alignment independent of the system dimensionality and which becomes exact only in the limit of infinite Prandtl number. Tumbling rates are extremely vigorous adjacent to the walls, where particles roughly perform Jeffery orbits. This implies that the root-mean-square near-wall tumbling rates for spheres are much stronger than for rods, up to $\mathcal{O}(10)$ times at $Ra\simeq 10^9$. In the turbulent bulk the situation reverses and rods tumble slightly faster than isotropic particles, in agreement with earlier observations in two-dimensional turbulence.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluctuations and correlations of reactive scalars near chemical equilibrium in incompressible turbulence

The statistical properties of species undergoing chemical reactions in a turbulent environment are studied. We focus on the case of reversible multi-component reactions of second and higher orders, in a condition close to chemical equilibrium sustained by random large-scale reactant sources, while the turbulent flow is highly developed. In such a state a competition exists between the chemical reaction that tends to dump reactant concentration fluctuations and enhance their correlation intensity and the turbulent mixing that on the contrary increases fluctuations and remove relative correlations. We show that a unique control parameter, the Damkhöler number ($Da_θ$) that can be constructed from the scalar Taylor micro-scale, the reactant diffusivity and the reaction rate characterises the functional dependence of fluctuations and correlations in a variety of conditions, i.e., at changing the reaction order, the Reynolds and the Schmidt numbers. The larger is such a Damkhöler number the more depleted are the scalar fluctuations as compared to the fluctuations of a passive scalar field in the same conditions, and vice-versa the more intense are the correlations. A saturation in this behaviour is observed beyond $Da_θ\simeq \mathcal{O}(10)$. We provide an analytical prediction for this phenomenon which is in excellent agreement with direct numerical simulation results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rotation of anisotropic particles in Rayleigh-Bénard turbulence

Inertialess anisotropic particles in a Rayleigh-Bénard turbulent flow show maximal tumbling rates for weakly oblate shapes, in contrast with the universal behaviour observed in developed turbulence where the mean tumbling rate monotonically decreases with the particle aspect ratio. This is due to the concurrent effect of turbulent fluctuations and of a mean shear flow whose intensity, we show, is determined by the kinetic boundary layers. In Rayleigh-Bénard turbulence prolate particles align preferentially with the fluid velocity, while oblate ones orient with the temperature gradient. This analysis elucidates the link between particle angular dynamics and small-scale properties of convective turbulence and has implications for the wider class of sheared turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Rotational dynamics of bottom-heavy rods in turbulence from experiments and numerical simulations

We successfully perform the three-dimensional tracking in a turbulent fluid flow of small asymmetrical particles that are neutrally-buoyant and bottom-heavy, i.e., they have a non-homogeneous mass distribution along their symmetry axis. We experimentally show how a tiny mass inhomogeneity can affect the particle orientation along the preferred vertical direction and modify its tumbling rate. The experiment is complemented by a series of simulations based on realistic Navier-Stokes turbulence and on a point-like particle model that is capable to explore the full range of parameter space characterized by the gravitational torque stability number and by the particle aspect ratio. We propose a theoretical perturbative prediction valid in the high bottom-heaviness regime that agrees well with the observed preferential orientation and tumbling rate of the particles. We also show that the heavy-tail shape of the probability distribution function of the tumbling rate is weakly affected by the bottom-heaviness of the particles.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Modelling sea ice and melt ponds evolution: sensitivity to microscale heat transfer mechanisms

We present a mathematical model describing the evolution of sea ice and meltwater during summer. The system is described by two coupled partial differential equations for the ice thickness $h$ and pond depth $w$ fields. We test the sensitivity of the model to variations of parameters controlling fluid-dynamic processes at the pond level, namely the variation of turbulent heat flux with pond depth and the lateral melting of ice enclosing a pond. We observe that different heat flux scalings determine different rates of total surface ablations, while the system is relatively robust in terms of probability distributions of pond surface areas. Finally, we study pond morphology in terms of fractal dimensions, showing that the role of lateral melting is minor, whereas there is evidence of an impact from the initial sea ice topography.

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
preprint2012arXiv

How gravity and size affect the acceleration statistics of bubbles in turbulence

We report results from the first systematic Lagrangian experimental investigation in the regime of very light (air bubbles in water) and large particles (D/η >> 1) in turbulence. Using a traversing camera setup and particle tracking, we study the Lagrangian acceleration statistics of ~3 mm diameter (D) bubbles in a water tunnel with nearly homogeneous and isotropic turbulence generated by an active-grid. The Reynolds number (Re_λ) is varied from 145 to 230, resulting in size ratios, D/η in the range of 7.3--12.5, where η is the Kolmogorov length scale. The experiments reveal that gravity increases the acceleration variance and reduces the intermittency of the PDF in the vertical direction. Once the gravity offset is subtracted, the variances of both the horizontal and vertical acceleration components are about 5\pm2 times larger than the one measured in the same flow for fluid tracers. Moreover, for these light particles, the experimental acceleration PDF shows a substantial reduction in intermittency at growing size ratios, in contrast to neutrally buoyant or heavy particles. All these results are closely matched by numerical simulations of finite-size bubbles with the Faxén corrections.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Rotation Rate of Rods in Turbulent Fluid Flow

The rotational dynamics of anisotropic particles advected in a turbulent fluid flow are important in many industrial and natural setting. Particle rotations are controlled by small scale properties of turbulence that are nearly universal, and so provide a rich system where experiments can be directly compared with theory and simulations. Here we report the first three-dimensional experimental measurements of the orientation dynamics of rod-like particles as they are advected in a turbulent fluid flow. We also present numerical simulations that show good agreement with the experiments and allow extension to a wide range of particle shapes. Anisotropic tracer particles preferentially sample the flow since their orientations become correlated with the velocity gradient tensor. The rotation rate is heavily influenced by this preferential alignment, and the alignment depends strongly on particle shape.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Three-dimensional Lagrangian Voronoi analysis for clustering of particles and bubbles in turbulence

Three-dimensional Voronoi analysis is used to quantify the clustering of inertial particles in homogeneous isotropic turbulence using data from numerics and experiments. We study the clustering behavior at different density ratios and particle response times (i.e. Stokes numbers St). The Probability Density Functions (PDFs) of the Voronoi cell volumes of light and heavy particles show a different behavior from that of randomly distributed particles -i.e. fluid tracers-implying that clustering is present. The standard deviation of the PDF normalized by that of randomly distributed particles is used to quantify the clustering. Light particles show maximum clustering for St around 1-2. The results are consistent with previous investigations employing other approaches to quantify the clustering. We also present the joint PDFs of enstrophy and Voronoi volumes and their Lagrangian autocorrelations. The small Voronoi volumes of light particles correspond to regions of higher enstrophy than those of heavy particles, indicating that light particles cluster in higher vorticity regions. The Lagrangian temporal autocorrelation function of Voronoi volumes shows that the clustering of light particles lasts much longer than that of heavy or neutrally buoyant particles. Due to inertial effects, the Lagrangian autocorrelation time-scale of clustered light particles is even longer than that of the enstrophy of the flow itself.

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&#39;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

Dynamics of inertial particles in a turbulent von Karman flow

We study the dynamics of neutrally buoyant particles with diameters varying in the range [1, 45] in Kolmogorov scale units (eta) and Reynolds numbers based on Taylor scale between 580 and 1050. One component of the particles&#39; velocity is measured using extended Laser Doppler Velocimetry at the center of a Von-Karman (VK) flow, acceleration is derived by differentiation. We find that although the particle acceleration variance decreases with increasing their diameter with scaling close to (D/eta)^-2/3, in agreement with previous observations, the characteristic time of acceleration autocorrelation increases much strongly than previously reported, and linearly with D/eta. A new analysis of the probability density functions of the acceleration shows smaller wings for larger particles; the flatness indeed decreases as also expected from the behavior of Eulerian pressure increments in the VK flow. We contrast our measurements with former observations in wind-tunnel turbulent flows and numerical simulations, and discuss if the observed differences arise from inherent properties of the VK flow.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of trailing wake drag on the statistical properties and dynamics of finite-sized particle in turbulence

We study by means of an Eulerian-Lagrangian model the statistical properties of velocity and acceleration of a neutrally-buoyant finite-sized particle in a turbulent flow statistically homogeneous and isotropic. The particle equation of motion, beside added mass and steady Stokes drag, keeps into account the unsteady Stokes drag force - known as Basset-Boussinesq history force - and the non-Stokesian drag based on Schiller-Naumann parametrization, together with the finite-size Faxen corrections. We focus on the case of flow at low Taylor-Reynolds number, Re_lambda ~ 31, for which fully resolved numerical data which can be taken as a reference are available (Homann & Bec 651 81-91 J. Fluid Mech. (2010)). Remarkably, we show that while drag forces have always minor effects on the acceleration statistics, their role is important on the velocity behavior. We propose also that the scaling relations for the particle velocity variance as a function of its size, which have been first detected in fully resolved simulations, does not originate from inertial-scale properties of the background turbulent flow but it is likely to arise from the non-Stokesian component of the drag produced by the wake behind the particle. Furthermore, by means of comparison with fully resolved simulations, we show that the Faxen correction to the added mass has a dominant role in the particle acceleration statistics even for particle with size in the inertial range.

0 citations0 reviewsNo code linkedOpen access