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

Enrico Calzavarini

Enrico Calzavarini 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.softphysics.ao-phBiological Physicsphysics.comp-phphysics.geo-ph

Catalog footprint

What is connected

26works
7topics
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

26 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' 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
preprint2016arXiv

A Lagrangian model of copepod dynamics: Clustering by escape jumps in turbulence

Planktonic copepods are small crustaceans that have the ability to swim by quick powerful jumps. Such an aptness is used to escape from high shear regions, which may be caused either by flow per- turbations, produced by a large predator (i.e., fish larvae), or by the inherent highly turbulent dynamics of the ocean. Through a combined experimental and numerical study, we investigate the impact of jumping behaviour on the small-scale patchiness of copepods in a turbulent environment. Recorded velocity tracks of copepods displaying escape response jumps in still water are here used to define and tune a Lagrangian Copepod (LC) model. The model is further employed to simulate the behaviour of thousands of copepods in a fully developed hydrodynamic turbulent flow obtained by direct numerical simulation of the Navier-Stokes equations. First, we show that the LC velocity statistics is in qualitative agreement with available experimental observations of copepods in tur- bulence. Second, we quantify the clustering of LC, via the fractal dimension $D_2$. We show that $D_2$ can be as low as ~ 2.3 and that it critically depends on the shear-rate sensitivity of the proposed LC model, in particular it exhibits a minimum in a narrow range of shear-rate values. We further investigate the effect of jump intensity, jump orientation and geometrical aspect ratio of the copepods on the small-scale spatial distribution. At last, possible ecological implications of the observed clustering on encounter rates and mating success are discussed

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite Volume $vs.$ Streaming-based Lattice Boltzmann algorithm for fluid-dynamics simulations: a one-to-one accuracy and performance study

A new finite volume (FV) discretisation method for the Lattice Boltzmann (LB) equation which combines high accuracy with limited computational cost is presented. In order to assess the performance of the FV method we carry out a systematic comparison, focused on accuracy and computational performances, with the standard $streaming$ (ST) Lattice Boltzmann equation algorithm. To our knowledge such a systematic comparison has never been previously reported. In particular we aim at clarifying whether and in which conditions the proposed algorithm, and more generally any FV algorithm, can be taken as the method of choice in fluid-dynamics LB simulations. For this reason the comparative analysis is further extended to the case of realistic flows, in particular thermally driven flows in turbulent conditions. We report the first successful simulation of high-Rayleigh number convective flow performed by a Lattice Boltzmann FV based algorithm with wall grid refinement.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Microbubbles and microparticles are not faithful tracers of turbulent acceleration

We report on the Lagrangian statistics of acceleration of small (sub-Kolmogorov) bubbles and tracer particles with Stokes number St << 1 in turbulent flow. At decreasing Reynolds number, the bubble accelerations show deviations from that of tracer particles, i.e. they deviate from the Heisenberg-Yaglom prediction and show a quicker decorrelation despite their small size and minute St. Using direct numerical simulations, we show that these effects arise due the drift of these particles through the turbulent flow. We theoretically predict this gravity-driven effect for developed isotropic turbulence, with the ratio of Stokes to Froude number or equivalently the particle drift-velocity governing the enhancement of acceleration variance and the reductions in correlation time and intermittency. Our predictions are in good agreement with experimental and numerical results. The present findings are relevant to a range of scenarios encompassing tiny bubbles and droplets that drift through the turbulent oceans and the atmosphere. They also question the common usage of micro-bubbles and micro-droplets as tracers in turbulence research.

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

Axially-homogeneous Rayleigh-Benard convection in a cylindrical cell

Previous numerical studies have shown that the "ultimate regime of thermal convection" can be attained in a Rayleigh-Benard cell when the kinetic and thermal boundary layers are eliminated by replacing the walls with periodic boundary conditions (homogeneous Rayleigh-Benard convection). Then, the heat transfer scales like Nu ~ Ra^{1/2} and turbulence intensity as Re ~ Ra^{1/2}, where the Rayleigh number Ra indicates the strength of the driving force. However, experiments never operate in unbounded domains and it is important to understand how confinement might alter the approach to this ultimate regime. Here we consider homogeneous Rayleigh-Benard convection in a laterally confined geometry - a small aspect-ratio vertical cylindrical cell - and show evidence of the ultimate regime as Ra is increased: In spite of the confinement and the resulting kinetic boundary layers, we still find Nu ~ Re ~ Ra^{1/2}. The system supports exact solutions composed of modes of exponentially growing vertical velocity and temperature fields, with Ra as the critical parameter determining the properties of these modes. Counterintuitively, in the low Ra regime, or for very narrow cylinders, the numerical simulations are susceptible to these solutions which can dominate the dynamics and lead to very high and unsteady heat transfer. As Ra is increased, interaction between modes stabilizes the system, evidenced by the increasing homogeneity and reduced fluctuations in the r.m.s. velocity and temperature fields. We also test that physical results become independent of the periodicity length of the cylinder, a purely numerical parameter, as the aspect ratio is increased.

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

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' 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
preprint2008arXiv

Acceleration statistics of finite-sized particles in turbulent flow: the role of Faxen forces

The dynamics of particles in turbulence when the particle-size is larger than the dissipative scale of the carrier flow is studied. Recent experiments have highlighted signatures of particles finiteness on their statistical properties, namely a decrease of their acceleration variance, an increase of correlation times -at increasing the particles size- and an independence of the probability density function of the acceleration once normalized to their variance. These effects are not captured by point particle models. By means of a detailed comparison between numerical simulations and experimental data, we show that a more accurate model is obtained once Faxen corrections are included.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Dimensionality and morphology of particle and bubble clusters in turbulent flow

We conduct numerical experiments to investigate the spatial clustering of particles and bubbles in simulations of homogeneous and isotropic turbulence. Varying the Stokes parameter and the densities, striking differences in the clustering of the particles can be observed. To quantify these visual findings we use the Kaplan--Yorke dimension. This local scaling analysis shows a dimension of approximately 1.4 for the light bubble distribution, whereas the distribution of very heavy particles shows a dimension of approximately 2.4. However, clearly separate parameter combinations yield the same dimensions. To overcome this degeneracy and to further develop the understanding of clustering, we perform a morphological (geometrical and topological) analysis of the particle distribution. For such an analysis, Minkowski functionals have been successfully employed in cosmology, in order to quantify the global geometry and topology of the large-scale distribution of galaxies. In the context of dispersed multiphase flow, these Minkowski functionals -- being morphological order parameters -- allow us to discern the filamentary structure of the light particle distribution from the wall-like distribution of heavy particles around empty interconnected tunnels.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Flow organization in non-Oberbeck-Boussinesq Rayleigh-Benard convection in water

Non-Oberbeck-Boussinesq (NOB) effects on the flow organization in two-dimensional Rayleigh-Benard turbulence are numerically analyzed. The working fluid is water. We focus on the temperature profiles, the center temperature, the Nusselt number, and on the analysis of the velocity field. Several velocity amplitudes (or Reynolds numbers) and several kinetic profiles are introduced and studied; these together describe the various features of the rather complex flow organization. The results are presented both as functions of the Rayleigh number Ra (with Ra up to 10^8) for fixed temperature difference (Delta) between top and bottom plates and as functions of Delta ("non-Oberbeck-Boussinesqness") for fixed Ra with Delta up to 60 K. All results are consistent with the available experimental NOB data for the center temperature Tc and the Nusselt number ratio Nu_{NOB}/Nu_{OB} (the label OB meaning that the Oberbeck-Boussinesq conditions are valid). Beyond Ra ~ 10^6 the flow consists of a large diagonal center convection roll and two smaller rolls in the upper and lower corners. In the NOB case the center convection roll is still characterized by only one velocity scale.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Microbubbly drag reduction in Taylor-Couette flow in the wavy vortex regime

We investigate the effect of microbubbles on Taylor-Couette flow by means of direct numerical simulations. We employ an Eulerian-Lagrangian approach with a gas-fluid coupling based on the point-force approximation. Added mass, drag, lift, and gravity are taken into account in the modeling of the motion of the individual bubble. We find that very dilute suspensions of small non-deformable bubbles (volume void fraction below 1%, zero Weber number and bubble Reynolds number <10) induce a robust statistically steady drag reduction (up to 20%) in the so called wavy vortex flow regime (Re = 600-2500). The Reynolds number dependence of the normalized torque (the so-called Torque Reduction Ratio (TRR) which corresponds to the drag reduction) is consistent with a recent series of experimental measurements performed by Murai et al. (J. Phys. 14, 143 (2005)). Our analysis suggests that the physical mechanism for the torque reduction in this regime is due to the local axial forcing, induced by rising bubbles, that is able to break the highly dissipative Taylor wavy vortices in the system. We finally show that the lift force acting on the bubble is crucial in this process. When neglecting it, the bubbles preferentially accumulate near the inner cylinder and the bulk flow is less efficiently modified.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Evidences of Bolgiano scaling in 3D Rayleigh-Benard convection

We present new results from high-resolution high-statistics direct numerical simulations of a tri-dimensional convective cell. We test the fundamental physical picture of the presence of both a Bolgiano-like and a Kolmogorov-like regime. We find that the dimensional predictions for these two distinct regimes (characterized respectively by an active and passive role of the temperature field) are consistent with our measurements.

0 citations0 reviewsNo code linkedOpen access