Showing 1–50 of 73 results for author: Rosti, M E

2.5D turbulence in shear-thinning jets

Authors: Christian Amor, Giovanni Soligo, Andrea Mazzino, Marco Edoardo Rosti

Abstract: The dimensional transition in turbulent jets of a shear-thinning fluid is studied via direct numerical simulations. Our findings reveal that under vertical confinement, the flow exhibits a unique mixed-dimensional (or 2.5D) state, where large-scale two-dimensional and small-scale three-dimensional structures coexist. This transition from three-dimensional turbulence near the inlet to two-dimension… ▽ More The dimensional transition in turbulent jets of a shear-thinning fluid is studied via direct numerical simulations. Our findings reveal that under vertical confinement, the flow exhibits a unique mixed-dimensional (or 2.5D) state, where large-scale two-dimensional and small-scale three-dimensional structures coexist. This transition from three-dimensional turbulence near the inlet to two-dimensional dynamics downstream is dictated by the level of confinement: weak confinement guarantees turbulence to remain three-dimensional, whereas strong confinement forces the transition to two-dimensions; the mixed-dimensional state is observed for moderate confinement and it emerges as soon as flow scales are larger than the vertical length. In this scenario, we observed that the mixed-dimensional state is an overall more energetic state and it shows a multi-cascade process, where the direct cascade of energy at small scales and the direct cascade of enstrophy at large scales coexist. The results provide insights into the complex dynamics of confined turbulent flows, relevant in both natural and industrial settings. △ Less

Submitted 6 July, 2024; v1 submitted 1 July, 2024; originally announced July 2024.

From yield stress to elastic instabilities: Tuning the extensional behavior of elastoviscoplastic fluid

Authors: Mohamed S. Abdelgawad, Simon J. Haward, Amy Q. Shen, Marco E. Rosti

Abstract: In this study, we delve into the intricacies of elastoviscoplastic (EVP) fluids, particularly focusing on how polymer additives influence their extensional behavior. Our findings reveal that polymer additives significantly alter the extensional properties of the EVP fluids, such as relaxation time and extensional stresses, while having negligible impact on the shear rheology. Interestingly, the mo… ▽ More In this study, we delve into the intricacies of elastoviscoplastic (EVP) fluids, particularly focusing on how polymer additives influence their extensional behavior. Our findings reveal that polymer additives significantly alter the extensional properties of the EVP fluids, such as relaxation time and extensional stresses, while having negligible impact on the shear rheology. Interestingly, the modified fluids exhibit a transition from yield stress-like behavior to viscoelastic-like behavior under high extensional rates, ultimately leading to destabilization under extreme deformation. This research enhances the fundamental understanding of EVP fluids and highlights potential advancements in applications, especially in precision-demanding fields like 3D printing. △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: 20 pages, 15 figures

Finite-size inertial spherical particles in turbulence

Authors: Alessandro Chiarini, Marco Edoardo Rosti

Abstract: We investigate by direct numerical simulations the fluid-solid interaction of non-dilute suspensions of spherical particles moving in triperiodic turbulence, at the relatively large Reynolds number of $Re_λ\approx 400$. The solid-to-fluid density ratio is varied between $1.3$ and $100$, the particle diameter $D$ ranges between $16 \le D/η\le 123$ ($η$ is the Kolmogorov scale), and the volume fract… ▽ More We investigate by direct numerical simulations the fluid-solid interaction of non-dilute suspensions of spherical particles moving in triperiodic turbulence, at the relatively large Reynolds number of $Re_λ\approx 400$. The solid-to-fluid density ratio is varied between $1.3$ and $100$, the particle diameter $D$ ranges between $16 \le D/η\le 123$ ($η$ is the Kolmogorov scale), and the volume fraction of the suspension is $0.079$. Turbulence is sustained using the Arnold-Beltrami-Childress cellular-flow forcing. The influence of the solid phase on the largest and energetic scales of the flow changes with the size and density of the particles. Light and large particles modulate all scales in a isotropic way, while heavier and smaller particles modulate the largest scales of the flow towards an anisotropic state. Smaller scales are isotropic and homogeneous for all cases. The mechanism driving the energy transfer across scales changes with the size and the density of the particles. For large and light particles the energy transfer is only marginally influenced by the fluid-solid interaction. For small and heavy particles, instead, the classical energy cascade is subdominant at all scales, and the energy transfer is essentially driven by the fluid-solid coupling. The influence of the solid phase on the flow intermittency is also discussed. Besides, the collective motion of the particles and their preferential location in relation with properties of the carrier flow are analysed. The solid phase exhibits moderate clustering; for large particles the level of clustering decreases with their density, while for small particles it is maximum for intermediate values. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Journal ref: J. Fluid Mech. 988 (2024) A17

Laminar planar jets of elastoviscoplastic fluids

Authors: Giovanni Soligo, Marco Edoardo Rosti

Abstract: We perform numerical simulations of planar jets of elastoviscoplastic (EVP) fluid (\citet{saramito2007new} model) at low Reynolds number. Three different configurations are considered: $(i)$ EVP jet in EVP ambient fluid, $(ii)$ EVP jet in Newtonian ambient fluid (miscible), and $(iii)$ EVP jet in Newtonian ambient fluid (immiscible). We investigate the effect of the Bingham number, \textit{i.e.} o… ▽ More We perform numerical simulations of planar jets of elastoviscoplastic (EVP) fluid (\citet{saramito2007new} model) at low Reynolds number. Three different configurations are considered: $(i)$ EVP jet in EVP ambient fluid, $(ii)$ EVP jet in Newtonian ambient fluid (miscible), and $(iii)$ EVP jet in Newtonian ambient fluid (immiscible). We investigate the effect of the Bingham number, \textit{i.e.} of the dimensionless yield stress of the EVP fluid, on the jet dynamics, and find a good agreement with the scaling for laminar, Newtonian jets for the centerline velocity $u_c/u_0\propto (x/h)^{-1/3}$ and for the jet thickness $δ_{mom}/h\propto (x/h)^{2/3}$ at small Bingham number. This is lost once substantial regions of the fluid become unyielded, where we find that the spreading rate of the jet and the decay rate of the centerline velocity increase with the Bingham number, due to the regions of unyielded fluid inducing a blockage effect on the jet. The most striking difference among the three configurations we considered is the extent and position of the regions of unyielded fluid: large portions of ambient and jet fluids (in particular, away from the inlet) are unyielded for the EVP jet in EVP ambient fluid, whereas for the other two configurations, the regions of unyielded fluid are limited to the jet, as expected. We derive a power law scaling for the centerline yield variable and confirm it with the results from our numerical simulations. The yield variable determines the transition from the viscoelastic Oldroyd-B fluid (yielded) to the viscoelastic Kelvin-Voigt material (unyielded). △ Less

Submitted 5 April, 2024; originally announced April 2024.

The dynamics of fibers dispersed in viscoelastic turbulent flows

Authors: M. S. Aswathy, Marco E Rosti

Abstract: This study explores the dynamics of finite-size fibers suspended freely in a viscoelastic turbulent flow. For a fiber suspended in Newtonian flows, two different flapping regimes were identified previously by Rosti et al (2018). Here we explore, how the fiber dynamics is modified by the elasticity of the carrier fluid by performing Direct Numerical Simulations of a two-way coupled fiber-fluid syst… ▽ More This study explores the dynamics of finite-size fibers suspended freely in a viscoelastic turbulent flow. For a fiber suspended in Newtonian flows, two different flapping regimes were identified previously by Rosti et al (2018). Here we explore, how the fiber dynamics is modified by the elasticity of the carrier fluid by performing Direct Numerical Simulations of a two-way coupled fiber-fluid system in a parametric space spanning different Deborah numbers, fiber bending stiffness and the linear density difference between fiber and fluid. We examine how these parameters influence various fiber characteristics such as the frequency of flapping, curvature, and alignment with the fluid strain and polymer stretching directions. Results reveal that the neutrally-bouyant fibers, depending on their flexibility, oscillate with large and small time scales transpiring from the flow, but the smaller time-scales are suppressed as the polymer elasticity increases. Polymer stretching is uncommunicative to denser-than-fluid fibers, which flap with large time scales from the flow when flexible and with their natural frequency when rigid. Thus, the characteristic elastic time scale has a subdominant effect when the fibers are neutrally-bouyant, while its effect is absent when the fibers become more inertial. Additionally, we see that the inertial fibers have larger curvatures and are less responsive to the polymer presence, whereas the neutrally-bouyant fibers show quantitative changes. Also, the neutrally-bouyant fibers show a higher alignment with the polymer stretching directions compared to the denser ones. In a nutshell, the polymers exert a larger influence on neutrally-bouyant fibers compared to the denser ones. The study comprehensively addresses the interplay between polymer elasticity and the fiber structural properties in determining its response behaviour in an elasto-inertial turbulent flow. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Dynamics and fluid-structure interaction in turbulent flows within and above flexible canopies

Authors: Giulio Foggi Rota, Alessandro Monti, Stefano Olivieri, Marco Edoardo Rosti

Abstract: Flexible canopy flows are often encountered in natural scenarios, e.g., when crops sway in the wind or when submerged kelp forests are agitated by marine currents. Here, we provide a detailed characterisation of the turbulent flow developed above and between the flexible filaments of a fully submerged dense canopy and we describe their dynamical response to the turbulent forcing. We investigate a… ▽ More Flexible canopy flows are often encountered in natural scenarios, e.g., when crops sway in the wind or when submerged kelp forests are agitated by marine currents. Here, we provide a detailed characterisation of the turbulent flow developed above and between the flexible filaments of a fully submerged dense canopy and we describe their dynamical response to the turbulent forcing. We investigate a wide range of flexibilities, encompassing the case in which the filaments are completely rigid and standing upright as well as that where they are fully compliant to the flow and deflected in the streamwise direction. We are thus able to isolate the effect of the canopy flexibility on the drag and on the inner-outer flow interactions, as well as the two flapping regimes of the filaments already identified for a single fiber. Furthermore, we offer a detailed description of the Reynolds stresses throughout the wall-normal direction resorting to the Lumley triangle formalism, and we show the multi-layer nature of turbulence inside and above the canopy. The relevance of our investigation is thus twofold: the fundamental physical understanding developed here paves the way towards the investigation of more complex and realistic scenarios, while we also provide a thorough characterisation of the turbulent state that can prove useful in the development of accurate turbulence models for RANS and LES. △ Less

Submitted 13 May, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 33 pages, 23 figures

Forced and natural dynamics of a clamped flexible fiber in wall turbulence

Authors: Giulio Foggi Rota, Morie Koseki, Riya Agrawal, Stefano Olivieri, Marco Edoardo Rosti

Abstract: We characterize the dynamical behavior of a clamped flexible fiber immersed in wall turbulence over a wide range of natural frequencies by means of direct numerical simulations. Only two flapping states are possible: one where the fiber oscillates at the characteristic frequency of the largest turbulent eddies and another where the natural structural response dominates. The former is obtained in t… ▽ More We characterize the dynamical behavior of a clamped flexible fiber immersed in wall turbulence over a wide range of natural frequencies by means of direct numerical simulations. Only two flapping states are possible: one where the fiber oscillates at the characteristic frequency of the largest turbulent eddies and another where the natural structural response dominates. The former is obtained in the more flexible cases, while the latter in the more rigid ones. We observe here for the first time that in the turbulence dominated regime the fiber always sways at a frequency proportional to the largest scale of the flow, regardless of its structural parameters. The hindrance of the clamp to the wall prevents the synchronization of the fiber with turbulent eddies of comparable size. △ Less

Submitted 14 November, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 6+4 pages, 6+3 figures

Journal ref: Physical Review Fluids 9, L012601 (2024)

Anisotropic mean flow enhancement and anomalous transport of finite-size spherical particles in turbulent flows

Authors: Alessandro Chiarini, Ianto Cannon, Marco Edoardo Rosti

Abstract: We investigate the influence of dispersed solid spherical particles on the largest scales of the turbulent Arnold-Beltrami-Childress (ABC) flow. The ABC flow is an ideal instance of a complex flow: it does not have solid boundaries, but possesses an inhomogeneous and three-dimensional mean shear. By tuning the parameters of the suspension, we show that particles deviate towards quasi-straight traj… ▽ More We investigate the influence of dispersed solid spherical particles on the largest scales of the turbulent Arnold-Beltrami-Childress (ABC) flow. The ABC flow is an ideal instance of a complex flow: it does not have solid boundaries, but possesses an inhomogeneous and three-dimensional mean shear. By tuning the parameters of the suspension, we show that particles deviate towards quasi-straight trajectories and exhibit anomalous transport. In doing this, they enhance the mean-flow energy and modulate the largest scales of the flow towards an anisotropic and quasi two-dimensional state. △ Less

Submitted 27 December, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Spheres and fibres in turbulent flows at various Reynolds numbers

Authors: Ianto Cannon, Stefano Olivieri, Marco E. Rosti

Abstract: We perform fully coupled numerical simulations using immersed boundary methods of finite-size spheres and fibres suspended in a turbulent flow for a range of Taylor Reynolds numbers $12.8<Re_λ<442$ and solid mass fractions $0\leq M\leq1$. Both spheres and fibres reduce the turbulence intensity with respect to the single-phase flow at all Reynolds numbers, with fibres causing a more significant red… ▽ More We perform fully coupled numerical simulations using immersed boundary methods of finite-size spheres and fibres suspended in a turbulent flow for a range of Taylor Reynolds numbers $12.8<Re_λ<442$ and solid mass fractions $0\leq M\leq1$. Both spheres and fibres reduce the turbulence intensity with respect to the single-phase flow at all Reynolds numbers, with fibres causing a more significant reduction than the spheres. The particles' effect on the anomalous dissipation tends to vanish as $Re \to \infty$. A scale-by-scale analysis shows that both particle shapes provide a "spectral shortcut" to the flow, but the shortcut extends further into the dissipative range in the case of fibres. Multifractal spectra of the near-particle dissipation show that spheres enhance dissipation in two-dimensional sheets, and fibres enhance the dissipation in structures with a dimension greater than one and less than two. In addition, we show that spheres suppress vortical flow structures, whereas fibres produce structures which completely overcome the turbulent vortex stretching behaviour in their vicinity. △ Less

Submitted 15 May, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: 18 pages, 12 figures

MSC Class: 76F05 ACM Class: J.2

A unified view of elastic and elasto-inertial turbulence in channel flows

Authors: Giulio Foggi Rota, Christian Amor, Soledad Le Clainche, Marco Edoardo Rosti

Abstract: The chaotic flow of elastic fluids at low Reynolds number (Re) is typically distinguished into elasto-inertial and elastic turbulence (EIT/ET). However, the clear separation among these two turbulent regimes in parallel flows with a gradual Re decrease remains elusive. Here, spanning Re over four orders of magnitude, we investigate the statistical and structural nature of the flows computed by ful… ▽ More The chaotic flow of elastic fluids at low Reynolds number (Re) is typically distinguished into elasto-inertial and elastic turbulence (EIT/ET). However, the clear separation among these two turbulent regimes in parallel flows with a gradual Re decrease remains elusive. Here, spanning Re over four orders of magnitude, we investigate the statistical and structural nature of the flows computed by fully resolved 3D numerical simulations, revealing that as inertia becomes sub-dominant all flows share similar dynamical properties. Supported by the results, we thus propose a unified view of fully developed ET and EIT in parallel flows. △ Less

Submitted 3 July, 2024; v1 submitted 8 October, 2023; originally announced October 2023.

Comments: 5+6 pages, 3+4 figures

Bridging Polymeric Turbulence at different Reynolds numbers: From Multiscaling to Multifractality

Authors: Rahul K. Singh, Marco E. Rosti

Abstract: The addition of polymers modifies a flow in a non-trivial way that depends on fluid inertia (given by the Reynolds number Re) and polymer elasticity (quantified by the Deborah number De). Using direct numerical simulations, we show that polymeric flows exhibit a Re and De dependent multiscaling energy spectrum. The different scaling regimes are tied to various dominant contributions -- fluid, poly… ▽ More The addition of polymers modifies a flow in a non-trivial way that depends on fluid inertia (given by the Reynolds number Re) and polymer elasticity (quantified by the Deborah number De). Using direct numerical simulations, we show that polymeric flows exhibit a Re and De dependent multiscaling energy spectrum. The different scaling regimes are tied to various dominant contributions -- fluid, polymer, and dissipation -- to the total energy flux across the scales. At small scales, energy is dissipated away by both polymers and the fluid. Fluid energy dissipation, in particular, is shown to be more intermittent in the presence of polymers, especially at small Re. The more intermittent, singular nature of energy dissipation is revealed clearly by the multifractal spectrum. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Intermittency in the not-so-smooth elastic turbulence

Authors: Rahul K. Singh, Prasad Perlekar, Dhrubaditya Mitra, Marco E. Rosti

Abstract: Elastic turbulence is the chaotic fluid motion resulting from elastic instabilities due to the addition of polymers in small concentrations at very small Reynolds ($\mbox{Re}$) numbers. Our direct numerical simulations show that elastic turbulence, though a low $\mbox{Re}$ phenomenon, has more in common with classical, Newtonian turbulence than previously thought. In particular, we find power-law… ▽ More Elastic turbulence is the chaotic fluid motion resulting from elastic instabilities due to the addition of polymers in small concentrations at very small Reynolds ($\mbox{Re}$) numbers. Our direct numerical simulations show that elastic turbulence, though a low $\mbox{Re}$ phenomenon, has more in common with classical, Newtonian turbulence than previously thought. In particular, we find power-law spectra for kinetic energy $E(k) \sim k^{-4}$ and polymeric energy $E_{\rm p}(k) \sim k^{-3/2}$, independent of the Deborah ($\mbox{De}$) number. This is further supported by calculation of scale-by-scale energy budget which shows a balance between the viscous term and the polymeric term in the momentum equation. In real space, as expected, the velocity field is smooth, i.e., the velocity difference across a length scale $r$, $δu \sim r$ but, crucially, with a non-trivial sub-leading contribution $r^{3/2}$ which we extract by using the second difference of velocity. The structure functions of second difference of velocity up to order $6$ show clear evidence of intermittency/multifractality. We provide additional evidence in support of this intermittent nature by calculating moments of rate of dissipation of kinetic energy averaged over a ball of radius $r$, $\varepsilon_{r}$, from which we compute the multifractal spectrum. △ Less

Submitted 23 April, 2024; v1 submitted 14 August, 2023; originally announced August 2023.

Morphology of clean and surfactant-laden droplets in homogeneous isotropic turbulence

Authors: Ianto Cannon, Giovanni Soligo, Marco E. Rosti

Abstract: We perform direct numerical simulations of surfactant-laden droplets in homogeneous isotropic turbulence with Taylor Reynolds number of 180. The droplets are modelled using the volume of fluid method, and the soluble surfactant is transported using an advection-diffusion equation. Effects of surfactant on the droplet and local flow statistics are well approximated using a lower, averaged value of… ▽ More We perform direct numerical simulations of surfactant-laden droplets in homogeneous isotropic turbulence with Taylor Reynolds number of 180. The droplets are modelled using the volume of fluid method, and the soluble surfactant is transported using an advection-diffusion equation. Effects of surfactant on the droplet and local flow statistics are well approximated using a lower, averaged value of surface tension, allowing us to extend the framework developed by Hinze (1955) and Kolmogorov (1949) for surfactant-free bubbles. The Kolmogorov-Hinze scale $d_H$ is indeed found to be a pivotal length scale in the droplets' dynamics, separating the coalescence-dominated (droplets smaller than $d_H$) and the breakage-dominated (droplets larger than $d_H$) regimes in the droplet size distribution. We find that droplets smaller than $d_H$ have compact, regular, spheroid-like shapes, whereas droplets larger than $d_H$ have long, convoluted, filamentous shapes with a diameter equal to $d_H$. This results in two different scaling laws for the interfacial area of the droplet. The normalised area, $A/d_H^2$, of droplets smaller than $d_H$ is proportional to $d^2$, while the area of droplets larger than $d_H$ is proportional to $d^3$, where $d$ is the droplet size. We further characterise the large filamentous droplets by computing the number of handles (loops of the dispersed phase extending into the carrier phase) and voids (regions of the carrier fluid entirely enclosed by the dispersed phase) on each droplet. The number of handles per unit length of filament scales inversely with surface tension. The number of voids is proportional to the droplet size and independent of surface tension. Handles are indeed an unstable feature of the interface and are destroyed by the restoring effect of surface tension, whereas voids can move freely in the interior of the droplets, unaffected by surface tension. △ Less

Submitted 14 April, 2024; v1 submitted 28 July, 2023; originally announced July 2023.

Comments: 33 pages, 13 figures

MSC Class: 76T06

Journal ref: J. Fluid Mech. 987 (2024) A31

On-off pumping for drag reduction in a turbulent channel flow

Authors: Giulio Foggi Rota, Alessandro Monti, Marco E. Rosti, Maurizio Quadrio

Abstract: We show that the energy required by a turbulent flow to displace a given amount of fluid through a straight duct in a given time interval can be reduced by modulating in time the pumping power. The control strategy is hybrid: it is passive, as it requires neither a control system nor control energy, but it manipulates how pumping energy is delivered to the system (as in active techniques) to incre… ▽ More We show that the energy required by a turbulent flow to displace a given amount of fluid through a straight duct in a given time interval can be reduced by modulating in time the pumping power. The control strategy is hybrid: it is passive, as it requires neither a control system nor control energy, but it manipulates how pumping energy is delivered to the system (as in active techniques) to increase the pumping efficiency. Our control employs a temporally periodic pumping pattern, where a short and intense acceleration (in which the pumping system is on) followed by a longer deceleration (in which the pumping system is off) makes the flow alternately visit a quasi-laminar and a turbulent state. The computational study is for a plane channel flow, and employs direct numerical simulations, which present specific computational challenges, for example the highly varying instantaneous value of the Reynolds number, and the importance of discretisation effects. Particular care is devoted to a meaningful definition of drag reduction in the present context. The ability of the forcing to yield significant savings is demonstrated. Since only a small portion of the parameter space is investigated, the best performance of the control technique remains to be assessed. △ Less

Submitted 29 April, 2023; originally announced May 2023.

Comments: Submitted to J. Fluid Mech

Journal ref: Journal of Fluid Mechanics 966, A12 (2023)

Collective dynamics of dense hairy surfaces in turbulent flow

Authors: Alessandro Monti, Stefano Olivieri, Marco E. Rosti

Abstract: Flexible filamentous beds interacting with a turbulent flow represent a fundamental setting for many environmental phenomena, e.g., aquatic canopies in marine current. Exploiting direct numerical simulations at high Reynolds number where the canopy stems are modelled individually, we provide evidence on the essential features of the honami/monami collective motion experienced by hairy surfaces ove… ▽ More Flexible filamentous beds interacting with a turbulent flow represent a fundamental setting for many environmental phenomena, e.g., aquatic canopies in marine current. Exploiting direct numerical simulations at high Reynolds number where the canopy stems are modelled individually, we provide evidence on the essential features of the honami/monami collective motion experienced by hairy surfaces over a range of different flexibilities, i.e., Cauchy number. Our findings clearly confirm that the collective motion is essentially driven by fluid flow turbulence, with the canopy having in this respect a fully-passive behavior. Instead, some features pertaining to the structural response turn out to manifest in the motion of the individual canopy elements when focusing, in particular, on the spanwise oscillation and/or on sufficiently small Cauchy numbers. △ Less

Submitted 27 March, 2023; v1 submitted 17 March, 2023; originally announced March 2023.

Enhanced axial migration of a deformable capsule in pulsatile channel flows

Authors: Naoki Takeishi, Marco Edoardo Rosti

Abstract: We present numerical analysis of the lateral movement of a deformable spherical capsule in a pulsatile channel flow, with a Newtonian fluid in almost inertialess condition and at a small confinement ratio $a_0/R = 0.4$, where $R$ and $a$ are the channel and capsule radius. We find that the speed of the axial migration of the capsule can be accelerated by the flow pulsation at a specific frequency.… ▽ More We present numerical analysis of the lateral movement of a deformable spherical capsule in a pulsatile channel flow, with a Newtonian fluid in almost inertialess condition and at a small confinement ratio $a_0/R = 0.4$, where $R$ and $a$ are the channel and capsule radius. We find that the speed of the axial migration of the capsule can be accelerated by the flow pulsation at a specific frequency. The migration speed increases with the oscillatory amplitude, while the most effective frequency remains basically unchanged and independent of the amplitude. Our numerical results form a fundamental basis for further studies on cellular flow mechanics, since pulsatile flows are physiologically relevant in human circulation, potentially affecting the dynamics of deformable particles and red blood cells (RBCs), and can also be potentially exploited in cell focusing techniques. △ Less

Submitted 5 May, 2023; v1 submitted 10 January, 2023; originally announced January 2023.

Saving Energy in Turbulent Flows with Unsteady Pumping

Authors: Giulio Foggi Rota, Alessandro Monti, Marco E. Rosti, Maurizio Quadrio

Abstract: Viscous dissipation causes significant energy losses in fluid flows; in ducts, laminar flows provide the minimum resistance to the motion, whereas turbulent currents substantially increase the friction at the wall and the energy requirement for pumping. Great effort is currently being devoted to find new ways of reducing the energy losses induced by turbulence. Here we propose a simple and novel d… ▽ More Viscous dissipation causes significant energy losses in fluid flows; in ducts, laminar flows provide the minimum resistance to the motion, whereas turbulent currents substantially increase the friction at the wall and the energy requirement for pumping. Great effort is currently being devoted to find new ways of reducing the energy losses induced by turbulence. Here we propose a simple and novel drag-reduction technique which achieves substantial energy savings in internal flows, even after the energy cost of the control is accounted for. Our approach consists in driving the flow with a pulsed pumping, unlike the common practice of a constant pumping. We alternate "pump on" phases where the flow accelerates, and "pump off" phases where the flow decays freely. The flow cyclically enters a quasi-laminar state during the acceleration, and transitions to a more classic turbulent state during the deceleration. Our numerical results demonstrate that important net energy savings can be achieved by simply modulating the power injection in the system over time. The physical understanding we have achieved can help the industry in reducing the waste of energy, making it cheaper and preserving the environment by reducing harmful emissions. △ Less

Submitted 19 January, 2023; v1 submitted 17 December, 2022; originally announced December 2022.

Comments: 12 pages, 3 figures

Journal ref: Scientific Reports 13, 1299 (2023)

Scaling and intermittency in turbulent flows of elastoviscoplastic fluids

Authors: Mohamed S. Abdelgawad, Ianto Cannon, Marco E. Rosti

Abstract: Non-Newtonian fluids have a viscosity that varies with applied stress. Elastoviscoplastic fluids, the elastic, viscous and plastic properties of which are interconnected in a non-trivial way, belong to this category. We have performed numerical simulations to investigate turbulence in elastoviscoplastic fluids at very high Reynolds-number values, as found in landslides and lava flows, focusing on… ▽ More Non-Newtonian fluids have a viscosity that varies with applied stress. Elastoviscoplastic fluids, the elastic, viscous and plastic properties of which are interconnected in a non-trivial way, belong to this category. We have performed numerical simulations to investigate turbulence in elastoviscoplastic fluids at very high Reynolds-number values, as found in landslides and lava flows, focusing on the effect of plasticity. We find that the range of active scales in the energy spectrum reduces when increasing the fluid plasticity; when plastic effects dominate, a new scaling range emerges between the inertial range and the dissipative scales. An extended self-similarity analysis of the structure functions reveals that intermittency is present and grows with the fluid plasticity. The enhanced intermittency is caused by the non-Newtonian dissipation rate, which also exhibits an intermittent behaviour. These findings have relevance to catastrophic events in natural flows, such as landslides and lava flows, where the enhanced intermittency results in stronger extreme events, which are thus more destructive and difficult to predict. △ Less

Submitted 4 March, 2023; v1 submitted 14 December, 2022; originally announced December 2022.

Comments: 12 pages and 10 figures, including supplemental material

MSC Class: 76A05; 76F05 ACM Class: J.2

Transient inverse energy cascade in free surface turbulence

Authors: G. Boffetta, A. Mazzino, S. Musacchio, M. E. Rosti

Abstract: We study the statistics of free-surface turbulence at large Reynolds numbers produced by direct numerical simulations in a fluid layer at different thickness with fixed characteristic forcing scale. We observe the production of a transient inverse cascade, with a duration which depends on the thickness of the layer, followed by a transition to three-dimensional turbulence initially produced close… ▽ More We study the statistics of free-surface turbulence at large Reynolds numbers produced by direct numerical simulations in a fluid layer at different thickness with fixed characteristic forcing scale. We observe the production of a transient inverse cascade, with a duration which depends on the thickness of the layer, followed by a transition to three-dimensional turbulence initially produced close to the bottom, no-slip boundary. By switching off the forcing, we study the decaying turbulent regime and we find that it cannot be described by an exponential law. Our results show that boundary conditions play a fundamental role in the nature of turbulence produced in thin layers and give limits on the conditions to produce a two-dimensional phenomenology. △ Less

Submitted 8 December, 2022; originally announced December 2022.

Comments: 9 pages, 8 figures

The role of wetting on the flow of two immiscible fluids in porous media

Authors: Armin Shahmardi, Salar Zamani Salimi, Outi Tammisola, Luca Brandt, Marco Edoardo Rosti

Abstract: We study the role of the capillary number, $Ca$ and of the surface wettability on the dynamics of the interface between an invading and a defending phase in a porous medium by means of numerical simulations. We employ a hybrid phase field-immersed boundary approach to successfully model the contact line dynamics over the solid objects. Using a phase-field method which naturally incorporates dynami… ▽ More We study the role of the capillary number, $Ca$ and of the surface wettability on the dynamics of the interface between an invading and a defending phase in a porous medium by means of numerical simulations. We employ a hybrid phase field-immersed boundary approach to successfully model the contact line dynamics over the solid objects. Using a phase-field method which naturally incorporates dynamic wetting we eliminate the need for empirical contact line models to address contact line singularity. We map the two dominant modes governing the motion of the interface, namely, capillary fingering, and stable penetration, in the ($Ca$-$θ$) plane, with $θ$ the static contact angle prescribed at the solid pores. Capillary fingering dominates at lower values of $Ca$ and pores hydrophobic to the invading phase, while a stable penetration is observed on hydrophillic surfaces. We present new measurements and analyses, including curvature probability density functions (pdf) and average curvature. We also show that the pressure needed for the invading phase to advance at constant flow rate decreases with the capillary number, and increases with the contact angle at the capillary numbers considered. The latter is due to a significant increase of the length of the interface in the case of capillary fingering. Finally, we show that it is possible to identify the different interfacial modes by measuring the penetration length and velocity during the medium filling. △ Less

Submitted 30 May, 2024; v1 submitted 22 October, 2022; originally announced October 2022.

Comments: 42 pages, 14 figures

The effect of particle anisotropy on the modulation of turbulent flows

Authors: Stefano Olivieri, Ianto Cannon, Marco E. Rosti

Abstract: We investigate the modulation of turbulence caused by the presence of finite-size dispersed particles. Bluff (isotropic) spheres vs slender (anisotropic) fibers are considered to understand the influence of the object shape on altering the carrier flow. While at a fixed mass fraction - but different Stokes number - both objects provide a similar bulk effect characterized by a large-scale energy de… ▽ More We investigate the modulation of turbulence caused by the presence of finite-size dispersed particles. Bluff (isotropic) spheres vs slender (anisotropic) fibers are considered to understand the influence of the object shape on altering the carrier flow. While at a fixed mass fraction - but different Stokes number - both objects provide a similar bulk effect characterized by a large-scale energy depletion, a scale-by-scale analysis of the energy transfer reveals that the alteration of the whole spectrum is intrinsically different. For bluff objects, the classical energy cascade is shrinked in its extension but unaltered in the energy content and its typical features, while for slender ones we find an alternative energy flux which is essentially mediated by the fluid-solid coupling. △ Less

Submitted 4 October, 2022; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 11 pages, 6 figures

MSC Class: 76F05 ACM Class: J.2

Viscoelasticity of suspension of red blood cells under oscillatory shear flow

Authors: Naoki Takeishi, Marco Edoardo Rosti, Naoto Yokoyama, Luca Brandt

Abstract: We present a numerical analysis of the rheology of a suspension of red blood cells (RBCs) for different volume fractions in a wall-bounded, effectively inertialess, small amplitude oscillatory shear (SAOS) flow for a wide range of applied frequencies. The RBCs are modeled as biconcave capsules, whose membrane is an isotropic and hyperelastic material following the Skalak constitutive law. The freq… ▽ More We present a numerical analysis of the rheology of a suspension of red blood cells (RBCs) for different volume fractions in a wall-bounded, effectively inertialess, small amplitude oscillatory shear (SAOS) flow for a wide range of applied frequencies. The RBCs are modeled as biconcave capsules, whose membrane is an isotropic and hyperelastic material following the Skalak constitutive law. The frequency-dependent viscoelasticity in the bulk suspension is quantified by the complex viscosity, defined by the amplitude of the particle shear stress and the phase difference between the stress and shear. SAOS flow basically impedes the deformation of individual RBCs as well as the magnitude of fluid-membrane interactions, resulting in a lower specific viscosity and first and second normal stress differences than in steady shear flow. Although it is known that the RBC deformation alone is sufficient to give rise to shear-thinning, our results show that the complex viscosity weakly depends on the frequency-modulated deformations or orientations of individual RBCs, but rather depends on combinations of the frequency-dependent amplitude and phase difference. The effect of the viscosity ratio between the cytoplasm and plasma and of the capillary number are also assessed. △ Less

Submitted 17 March, 2024; v1 submitted 19 July, 2022; originally announced July 2022.

Dense bidisperse suspensions under non-homogenous shear

Authors: Alessandro Monti, Marco Edoardo Rosti

Abstract: We study the rheological behaviour of bidisperse suspensions in three dimensions under a non-uniform shear flow, made by the superimposition of a linear shear and a sinusoidal disturbance. Our results show that i) only a streamwise disturbance in the shear-plane alters the suspension dynamics by substantially reducing the relative viscosity, ii) with the amplitude of the disturbance determining a… ▽ More We study the rheological behaviour of bidisperse suspensions in three dimensions under a non-uniform shear flow, made by the superimposition of a linear shear and a sinusoidal disturbance. Our results show that i) only a streamwise disturbance in the shear-plane alters the suspension dynamics by substantially reducing the relative viscosity, ii) with the amplitude of the disturbance determining a threshold value for the effect to kick-in and its wavenumber controlling the amount of reduction and which of the two phases is affected. We show that, iii) the rheological changes are caused by the effective separation of the two phases, with the large or small particles layering in separate regions. We provide a physical explanation of the phase separation process and of the conditions necessary to trigger it. We test the results in the whole flow curve, and we show that the mechanism remains substantially unaltered, with the only difference being the nature of the interactions between particles modified by the phase separation. △ Less

Submitted 8 October, 2023; v1 submitted 16 June, 2022; originally announced June 2022.

Journal ref: Sci Rep 13, 14310 (2023)

Shear-thickening of dense bidispersed suspensions

Authors: Alessandro Monti, Marco Edoardo Rosti

Abstract: We study the rheological behaviour of a dense bidispersed suspension varying the relative size of the two dispersed phases. The main outcome of our analysis is that an enhanced flowability (reduced relative viscosity) of the suspension can be achieved by increasing the dispersion ratio of the phases. We explain the observed result by showing that the presence of large particles increases the packi… ▽ More We study the rheological behaviour of a dense bidispersed suspension varying the relative size of the two dispersed phases. The main outcome of our analysis is that an enhanced flowability (reduced relative viscosity) of the suspension can be achieved by increasing the dispersion ratio of the phases. We explain the observed result by showing that the presence of large particles increases the packing efficiency of the suspension, leading to a reduction of the contribution of the contacts on the overall viscosity of the suspension in the shear-thickening regime, i.e. where the contacts are the dominating component. △ Less

Submitted 19 February, 2023; v1 submitted 13 June, 2022; originally announced June 2022.

Non-Newtonian turbulent jets at low-Reynolds number

Authors: Giovanni Soligo, Marco Edoardo Rosti

Abstract: We perform direct numerical simulations of planar jets of non-Newtonian fluids at low Reynolds number, in typical laminar conditions for a Newtonian fluid. We select three different non-Newtonian fluid models characterized by shear-thinning (Carreau), viscoelasticity (Oldroyd-B) and shear-thinning and viscoelasticity together (Giesekus), and perform a thorough analysis of the resulting flow statis… ▽ More We perform direct numerical simulations of planar jets of non-Newtonian fluids at low Reynolds number, in typical laminar conditions for a Newtonian fluid. We select three different non-Newtonian fluid models characterized by shear-thinning (Carreau), viscoelasticity (Oldroyd-B) and shear-thinning and viscoelasticity together (Giesekus), and perform a thorough analysis of the resulting flow statistics. We observe that, as the Weissenberg number is increased, the jet transitions from a laminar flow at low Weissenberg number, to a turbulent flow at high Weissenberg number. We show that the different non-Newtonian features and their combination give rise to rather different flowing regimes, originating from the competition of viscous, elastic and inertial effects. We observe that both viscoelasticity and shear-thinning can develop the instability and the consequent transition to a turbulent flowing regime; however, the purely viscoelastic Oldroyd-B fluid exhibits the onset of disordered fluid motions at a lower Weissenberg number than what observed for the purely shear-thinning Carreau fluid. When the two effects are both present, an intermediate condition is found, suggesting that, in this case, the shear-thinning feature is acting against the fluid elasticity. Despite the qualitative differences observed in the flowing regime, the bulk statistics, namely the centerline velocity and jet thickness, follow almost the same power-law scalings obtained for laminar and turbulent Newtonian planar jets. △ Less

Submitted 8 June, 2023; v1 submitted 11 June, 2022; originally announced June 2022.

On the solidity parameter in canopy flows

Authors: Alessandro Monti, Shane Nicholas, Mohammad Omidyeganeh, Alfredo Pinelli, Marco E. Rosti

Abstract: We have performed high-fidelity simulations of turbulent open-channel flows over submerged rigid canopies made of cylindrical filaments of fixed length $l=0.25H$ ($H$ being the domain depth) mounted on the wall with an angle of inclination $θ$. The inclination is the free parameter that sets the density of the canopy by varying its frontal area. The density of the canopy, based on the solidity par… ▽ More We have performed high-fidelity simulations of turbulent open-channel flows over submerged rigid canopies made of cylindrical filaments of fixed length $l=0.25H$ ($H$ being the domain depth) mounted on the wall with an angle of inclination $θ$. The inclination is the free parameter that sets the density of the canopy by varying its frontal area. The density of the canopy, based on the solidity parameter $λ$, is a widely accepted criterion defining the ongoing canopy flow regime, with low values ($λ\ll 0.1$) indicating the sparse regime, and higher values ($λ> 0.1$) the dense regime. All the numerical predictions have been obtained considering the same nominal bulk Reynolds number (i.e. $Re_b=U_b H/ν= 6000$). We consider nine configurations of canopies, with $θ$ varying symmetrically around $0°$ in the range $θ\in [\pm 78.5°]$, where positive angles define canopies inclined in the flow direction (with the grain) and $θ=0°$ corresponds to the wall-normally mounted canopy. The study compares canopies with identical solidity obtained inclining the filaments in opposite angles and assesses the efficacy of the solidity as a representative parameter. It is found that when the canopy is inclined, the actual flow regime differs substantially from the one of a straight canopy that shares the same solidity indicating that criteria solely based on this parameter are not robust. Finally, a new phenomenological model describing the interaction between the coherent structures populating the canopy region and the outer flow is given. △ Less

Submitted 28 June, 2022; v1 submitted 16 May, 2022; originally announced May 2022.

FluTAS: A GPU-accelerated finite difference code for multiphase flows

Authors: Marco Crialesi-Esposito, Nicolo Scapin, Andreas D. Demou, Marco Edoardo Rosti, Pedro Costa, Filippo Spiga, Luca Brandt

Abstract: We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited… ▽ More We present the Fluid Transport Accelerated Solver, FluTAS, a scalable GPU code for multiphase flows with thermal effects. The code solves the incompressible Navier-Stokes equation for two-fluid systems, with a direct FFT-based Poisson solver for the pressure equation. The interface between the two fluids is represented with the Volume of Fluid (VoF) method, which is mass conserving and well suited for complex flows thanks to its capacity of handling topological changes. The energy equation is explicitly solved and coupled with the momentum equation through the Boussinesq approximation. The code is conceived in a modular fashion so that different numerical methods can be used independently, the existing routines can be modified, and new ones can be included in a straightforward and sustainable manner. FluTAS is written in modern Fortran and parallelized using hybrid MPI/OpenMP in the CPU-only version and accelerated with OpenACC directives in the GPU implementation. We present different benchmarks to validate the code, and two large-scale simulations of fundamental interest in turbulent multiphase flows: isothermal emulsions in HIT and two-layer Rayleigh-Bénard convection. FluTAS is distributed through a MIT license and arises from a collaborative effort of several scientists, aiming to become a flexible tool to study complex multiphase flows. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Transport and evaporation of virus-containing droplets exhaled by men and women in typical cough events

Authors: Stefano Olivieri, Mattia Cavaiola, Andrea Mazzino, Marco Edoardo Rosti

Abstract: The spreading of the virus-containing droplets exhaled during respiratory events, e.g., cough, is an issue of paramount importance for the prevention of many infections such as COVID-19. According to the scientific literature, remarkable differences can be ascribed to several parameters that govern such complex and multiphysical problem. Among these, a particular influence appears associated with… ▽ More The spreading of the virus-containing droplets exhaled during respiratory events, e.g., cough, is an issue of paramount importance for the prevention of many infections such as COVID-19. According to the scientific literature, remarkable differences can be ascribed to several parameters that govern such complex and multiphysical problem. Among these, a particular influence appears associated with the different airflows typical of male and female subjects. Focusing on a typical cough event, we investigate this aspect by means of highly-resolved direct numerical simulations of the turbulent airflow in combination with a comprehensive Lagrangian particle tracking model for the droplet motion and evaporation. We observe and quantify major differences between the case of male and female subjects, both in terms of the droplet final reach and evaporation time. Our results can be associated with the different characteristics in the released airflow and thus confirm the influence of the subject gender (or other physical properties providing different exhalation profiles) on both short-range and long-range airborne transmission. △ Less

Submitted 10 January, 2022; originally announced January 2022.

An Eulerian-based immersed boundary method for particle suspensions with implicit lubrication model

Authors: Naoki Hori, Marco Edoardo Rosti, Shu Takagi

Abstract: We describe an immersed boundary method in which the fluid-structure coupling is achieved in an Eulerian framework. The method is an improved extension of the immersed boundary method originally developed by Kajishima et al. [1], which accounts for the inertia of the fictitious fluid inside the particle volume and is thus able to reproduce the behaviour of particles both in the case of neutrally-b… ▽ More We describe an immersed boundary method in which the fluid-structure coupling is achieved in an Eulerian framework. The method is an improved extension of the immersed boundary method originally developed by Kajishima et al. [1], which accounts for the inertia of the fictitious fluid inside the particle volume and is thus able to reproduce the behaviour of particles both in the case of neutrally-buoyant objects and in the presence of density difference between the particles and the fluid. The method is capable to handle the presence of multiple suspended objects, i.e., a suspension, by including a soft-sphere normal collision model, while the lubrication correction typically added to similar immersed boundary methods in order to capture the sub-grid unresolved lubrication force is here treated implicitly, i.e., naturally obtained without any explicit expression, thus no additional computation is required. We show that our methodology can successfully reproduce the rheology of a particle suspension in a shear flow up to a dense regime (with a maximum particle volume fraction around 46%) without any additional correction force. The applicability of this methodology is also tested in a turbulent pressure-driven duct flow at high Reynolds number in the presence of non-negligible inertia and non-uniform shear-rate, showing good agreement with experimental measurements. △ Less

Submitted 10 January, 2022; originally announced January 2022.

On the fully coupled dynamics of flexible fibres dispersed in modulated turbulence

Authors: Stefano Olivieri, Andrea Mazzino, Marco E. Rosti

Abstract: The present work investigates the mechanical behaviour of finite-size, elastic and inertial fibres freely moving in a homogeneous and isotropic turbulent flow at moderate Reynolds number. Fully-resolved, direct numerical simulations, based on a finite difference discretisation and the immersed boundary method, are performed to mutually couple the dynamics of fibres and fluid turbulence, allowing t… ▽ More The present work investigates the mechanical behaviour of finite-size, elastic and inertial fibres freely moving in a homogeneous and isotropic turbulent flow at moderate Reynolds number. Fully-resolved, direct numerical simulations, based on a finite difference discretisation and the immersed boundary method, are performed to mutually couple the dynamics of fibres and fluid turbulence, allowing to account for the backreaction of the dispersed phase to the carrier flow. An extensive parametric study is carried out over the characteristic properties of the suspension, i.e., fibre's linear density (from neutrally-buoyant to heavy fibres), length (from short fibres comparable with the dissipative scale to long fibres comparable with the integral scale) and bending stiffness (from highly flexible to almost rigid fibres), as well as the concentration (from dilute to non-dilute suspensions). Results reveal the existence of a robust turbulence modulation mechanism that is primarily controlled by the mass fraction of the suspension (with only a minor influence of the fibre's bending stiffness), which is characterised in detail by means of a scale-by-scale analysis in Fourier space. Despite such alteration with respect to the single-phase case due to the non-negligible backreaction, fibres experience only two possible flapping states (previously identified in the very dilute condition) while being transported and deformed by the flow. Additionally, we show that the maximum curvature obeys to different scaling laws that can be derived from the fibre dynamical equation. Finally, we explore the preferential concentration and alignment of fibres within the flow, highlighting the peculiar role of inertia and elasticity. △ Less

Submitted 19 May, 2022; v1 submitted 23 November, 2021; originally announced November 2021.

Large is different: non-monotonic behaviour of elastic range scaling in polymeric turbulence at large Reynolds and Deborah numbers

Authors: Marco E Rosti, Prasad Perlekar, Dhrubaditya Mitra

Abstract: We use direct numerical simulations to study homogeneous, and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers. We find that for small wavenumbers $k$, the kinetic energy spectrum shows Kolmogorov--like behavior which crosses over at a larger $k$ to a novel, elastic scaling regime, $E(k) \sim k^{-ξ}$, with $ξ\approx 2.3$. We study the contribution of the p… ▽ More We use direct numerical simulations to study homogeneous, and isotropic turbulent flows of dilute polymer solutions at high Reynolds and Deborah numbers. We find that for small wavenumbers $k$, the kinetic energy spectrum shows Kolmogorov--like behavior which crosses over at a larger $k$ to a novel, elastic scaling regime, $E(k) \sim k^{-ξ}$, with $ξ\approx 2.3$. We study the contribution of the polymers to the flux of kinetic energy through scales, and find that it can be decomposed into two parts: one increase in effective viscous dissipation, and a purely elastic contribution that dominates over the nonlinear flux in the range of $k$ over which the elastic scaling is observed. The multiscale balance between the two fluxes determines the crossover wavenumber which depends non-monotically on the Deborah number. Consistently, structure functions also show two scaling ranges, with intermittency present in both of them in equal measure. △ Less

Submitted 12 February, 2023; v1 submitted 22 November, 2021; originally announced November 2021.

Modulation of homogeneous and isotropic turbulence in emulsions

Authors: Marco Crialesi-Esposito, Marco Edoardo Rosti, Sergio Chibbaro, Luca Brandt

Abstract: We present a numerical study of emulsions in homogeneous and isotropic turbulence at $Re_λ=137$. The problem is addressed via Direct Numerical Simulations (DNS), where the Volume of Fluid (VOF) is used to represent the complex features of the liquid-liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turb… ▽ More We present a numerical study of emulsions in homogeneous and isotropic turbulence at $Re_λ=137$. The problem is addressed via Direct Numerical Simulations (DNS), where the Volume of Fluid (VOF) is used to represent the complex features of the liquid-liquid interface. We consider a mixture of two iso-density fluids, where fluid properties are varied with the goal of understanding their role in turbulence modulation, in particular the volume fraction ($0.03<α<0.5$), viscosity ratio ($0.01<μ_d/μ_c<100$) and large scale Weber number ($10.6<We_\mathcal{L}<106.5$). The analysis, performed by studying integral quantities and spectral scale-by-scale analysis, reveals that energy is consistently transported from large to small scales by the interface, and no inverse cascade is observed. Furthermore, the total surface is found to be directly proportional to the amount of energy transported, while viscosity and surface tension alter the dynamic that regulates energy transport. We also observe the $-10/3$ and $-3/2$ scaling on droplet size distributions, suggesting that the dimensional arguments which led to their derivation are verified in HIT conditions. △ Less

Submitted 19 August, 2021; originally announced August 2021.

A fast and efficient tool to study the rheology of dense suspensions

Authors: Alessandro Monti, Vikram Rathee, Amy Q. Shen, Marco E. Rosti

Abstract: A cutting-edge software is presented to tackle the Newton-Euler equations governing the dynamics of granular flows and dense suspensions in Newtonian fluids. In particular, we propose an implementation of a fixed-radius near neighbours search based on an efficient counting sort algorithm with an improved symmetric search. The adopted search method drastically reduces the computational cost and all… ▽ More A cutting-edge software is presented to tackle the Newton-Euler equations governing the dynamics of granular flows and dense suspensions in Newtonian fluids. In particular, we propose an implementation of a fixed-radius near neighbours search based on an efficient counting sort algorithm with an improved symmetric search. The adopted search method drastically reduces the computational cost and allows an efficient parallelisation even on a single node through the multi-threading paradigm. Emphasis is also given to the memory efficiency of the code since the history of the contacts among particles has to be traced to model the frictional contributions, when dealing with granular flows of rheological interest that consider non-smooth interacting particles. An effective procedure based on an estimate of the maximum number of the smallest particles surrounding the largest one (given the radii distribution) and a sort applied only to the surrounding particles only is implemented, allowing us to effectively tackle the rheology of non-monodispersed particles with high size-ratio in large domains. Finally, we present validations and verification of the numerical procedure, by comparing with previous simulations and experiments, and present new software capabilities. △ Less

Submitted 29 September, 2021; v1 submitted 2 August, 2021; originally announced August 2021.

Journal ref: Physics of Fluids 33, 103314 (2021)

Shear-thinning and shear-thickening emulsions in shear flows

Authors: Marco Edoardo Rosti, Shu Takagi

Abstract: We study the rheology of a two-fluid emulsion in semi-concentrated conditions; the solute is Newtonian while the solvent an inelastic power law fluid. The problem at hand is tackled by means of direct numerical simulations using the volume of fluid method. The analysis is performed for different volume fractions and viscosity ratios under the assumption of negligible inertia and zero buoyancy forc… ▽ More We study the rheology of a two-fluid emulsion in semi-concentrated conditions; the solute is Newtonian while the solvent an inelastic power law fluid. The problem at hand is tackled by means of direct numerical simulations using the volume of fluid method. The analysis is performed for different volume fractions and viscosity ratios under the assumption of negligible inertia and zero buoyancy force. Several carrier fluids are considered encompassing both the shear-thinning and thickening behaviors. We show that the effective viscosity of the system increases for shear-thickening fluids and decreases for the shear-thinning ones for all the viscosity ratios considered. The changes in the emulsion viscosity are mainly due to modifications of the coalescence in the system obtained by changing the carrier fluid property: indeed, local large and low shear rates are found in the regions between two interacting droplets for shear-thickening and thinning fluids, respectively, resulting in increased and reduced local viscosity which ultimately affects the drainage time of the system. This process is independent of the nominal viscosity ratio of the two fluids and we show that it can not be understood by considering only the mean shear rate and viscosity of the two fluids across the domain, but the full spectrum of shear rate must be taken into account. △ Less

Submitted 1 August, 2021; originally announced August 2021.

Puff turbulence in the limit of strong buoyancy

Authors: Andrea Mazzino, Marco Edoardo Rosti

Abstract: We provide a numerical validation of a recently proposed phenomenological theory to characterize the space-time statistical properties of a turbulent puff, both in terms of bulk properties, such as the mean velocity, temperature and size, and scaling laws for velocity and temperature differences both in the viscous and in the inertial range of scales. In particular, apart from the more classical s… ▽ More We provide a numerical validation of a recently proposed phenomenological theory to characterize the space-time statistical properties of a turbulent puff, both in terms of bulk properties, such as the mean velocity, temperature and size, and scaling laws for velocity and temperature differences both in the viscous and in the inertial range of scales. In particular, apart from the more classical shear-dominated puff turbulence, our main focus is on the recently discovered new regime where turbulent fluctuations are dominated by buoyancy. The theory is based on an adiabaticity hypothesis which assumes that small-scale turbulent fluctuations rapidly relax to the slower large-scale dynamics, leading to a generalization of the classical Kolmogorov and Kolmogorov-Obukhov-Corrsin theories for a turbulent puff hosting a scalar field. We validate our theory by means of massive direct numerical simulations finding excellent agreement. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Unraveling the secrets of turbulence in a fluid puff

Authors: Andrea Mazzino, Marco E Rosti

Abstract: Turbulent puffs are ubiquitous in everyday life phenomena. Understanding their dynamics is important in a variety of situations ranging from industrial processes to pure and applied science. In all these fields, a deep knowledge of the statistical structure of temperature and velocity space/time fluctuations is of paramount importance to construct models of chemical reaction (in chemistry), of con… ▽ More Turbulent puffs are ubiquitous in everyday life phenomena. Understanding their dynamics is important in a variety of situations ranging from industrial processes to pure and applied science. In all these fields, a deep knowledge of the statistical structure of temperature and velocity space/time fluctuations is of paramount importance to construct models of chemical reaction (in chemistry), of condensation of virus-containing droplets (in virology and/or biophysics), and optimal mixing strategies in industrial applications. As a matter of fact, results of turbulence in a puff are confined to bulk properties (i.e. average puff velocity and typical decay/growth time) and dates back to the second half of the 20th century. There is thus a huge gap to fill to pass from bulk properties to two-point statistical observables. Here we fill this gap exploiting theory and numerics in concert to predict and validate the space/time scaling behaviors of both velocity and temperature structure functions including intermittency corrections. Excellent agreement between theory and simulations is found. Our results are expected to have profound impact to develop evaporation models for virus-containing droplets carried by a turbulent puff, with benefits to the comprehension of the airborne route of virus contagion. △ Less

Submitted 21 June, 2021; originally announced June 2021.

Journal ref: Phys. Rev. Lett. 127, 094501 (2021)

Universal flapping states of elastic fibers in modulated turbulence

Authors: Stefano Olivieri, Andrea Mazzino, Marco E. Rosti

Abstract: We study the fully-coupled dynamics between a fully-developed turbulent flow and an ensemble of immersed flexible fibers. We vary the concentration of the suspension, the mechanical properties and the length of the fibers in a vast parametric range. For all configurations, the fiber dynamics falls in only two possible dynamical states: (i) the fiber manifests its natural response to the flow forci… ▽ More We study the fully-coupled dynamics between a fully-developed turbulent flow and an ensemble of immersed flexible fibers. We vary the concentration of the suspension, the mechanical properties and the length of the fibers in a vast parametric range. For all configurations, the fiber dynamics falls in only two possible dynamical states: (i) the fiber manifests its natural response to the flow forcing or (ii) its motion fully synchronizes to the hydrodynamic timescales of the turbulent flow. This scenario holds for both a dilute condition, where the carrier flow is not affected by the fluid-structure interaction, as well as in the case where the flow is substantially altered by the presence of the immersed objects. Such backreaction effect can be macroscopically modelled in terms of the mass fraction of the suspension. Our results can be readily extended to any elastic objects interacting with fluid turbulence. △ Less

Submitted 21 August, 2021; v1 submitted 15 June, 2021; originally announced June 2021.

Journal ref: Physics of Fluids 33.7, 071704 (2021)

The effect of droplet coalescence on drag in turbulent channel flows

Authors: Ianto Cannon, Daulet Izbassarov, Outi Tammisola, Luca Brandt, Marco E. Rosti

Abstract: We study the effect of droplet coalescence on turbulent wall-bounded flows, by means of direct numerical simulations. In particular, the volume-of-fluid and front-tracking methods are used to simulate turbulent channel flows containing coalescing and non-coalescing droplets, respectively. We find that coalescing droplets have a negligible effect on the drag, whereas the non-coalescing ones steadil… ▽ More We study the effect of droplet coalescence on turbulent wall-bounded flows, by means of direct numerical simulations. In particular, the volume-of-fluid and front-tracking methods are used to simulate turbulent channel flows containing coalescing and non-coalescing droplets, respectively. We find that coalescing droplets have a negligible effect on the drag, whereas the non-coalescing ones steadily increase drag as the volume fraction of the dispersed phase increases: indeed, at 10\% volume fraction, the non-coalescing droplets show a 30\% increase in drag, whereas the coalescing droplets show less than 4\% increase. We explain this by looking at the wall-normal location of droplets in the channel and show that non-coalescing droplets enter the viscous sublayer, generating an interfacial shear stress which reduces the budget for viscous stress in the channel. On the other hand, coalescing droplets migrate towards the bulk of the channel forming large aggregates, which hardly affect the viscous shear stress while damping the Reynolds shear stress. We prove this by relating the mean viscous shear stress integrated in the wall-normal direction to the centreline velocity. △ Less

Submitted 19 July, 2021; v1 submitted 3 June, 2021; originally announced June 2021.

Comments: 11 pages, 7 figures

MSC Class: 76T10

Journal ref: Physics of Fluids 33, 085112 (2021)

Near-wall turbulence alteration with the Transpiration-Resistance Model

Authors: Seyed Morteza Habibi Khorasani, Uǧis Lācis, Simon Pasche, Marco Edoardo Rosti, Shervin Bagheri

Abstract: A set of boundary conditions called the Transpiration-Resistance Model (TRM) are investigated in altering near-wall turbulence. The TRM has been previously proposed by \citet{Lacis2020} as a means of representing the net effect of surface micro-textures on their overlying bulk flows. It encompasses conventional Navier-slip boundary conditions relating the streamwise and spanwise velocities to thei… ▽ More A set of boundary conditions called the Transpiration-Resistance Model (TRM) are investigated in altering near-wall turbulence. The TRM has been previously proposed by \citet{Lacis2020} as a means of representing the net effect of surface micro-textures on their overlying bulk flows. It encompasses conventional Navier-slip boundary conditions relating the streamwise and spanwise velocities to their respective shears through the slip lengths $\ell_x$ and $\ell_z$. In addition, it features a transpiration condition accounting for the changes induced in the wall-normal velocity by expressing it in terms of variations of the wall-parallel velocity shears through the transpiration lengths $m_x$ and $m_z$. Greater levels of drag increase occur when more transpiration takes place at the boundary plane, with turbulent transpiration being predominately coupled to the spanwise shear component for canonical near-wall turbulence. The TRM can reproduce the effect of a homogeneous and structured roughness up to ${k^+}\,{\approx18}$. In this transitionally rough flow regime, the transpiration lengths of the TRM must be empirically determined. The \emph{transpiration factor} is defined as the product between the slip and transpiration lengths, i.e. $(m\ell)_{x,z}$. This factor contains the compound effect of the wall-parallel velocity occurring at the boundary plane and increased permeability, both of which lead to the transport of momentum in the wall-normal direction. A linear relation between the transpiration factor and the roughness function is observed for regularly textured surfaces in the transitionally rough regime of turbulence. The relations obtained between the transpiration factor and the roughness function show that such effective flow quantities can be suitable measures for characterizing rough surfaces in this flow regime. △ Less

Submitted 23 October, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Journal ref: Journal of Fluid Mechanics , Volume 942 , 10 July 2022 , A45

Finite-size evaporating droplets in weakly compressible homogeneous shear turbulence

Authors: Nicolò Scapin, Federico Dalla Barba, Giandomenico Lupo, Marco Edoardo Rosti, Christophe Duwig, Luca Brandt

Abstract: We perform interface-resolved simulations of finite-size evaporating droplets in weakly-compressible homogeneous shear turbulence (HST). The study is conducted by varying three dimensionless physical parameters: the initial gas temperature over the critical temperature $T_{g,0}/T_c$, the initial droplet diameter over the Kolmogorov scale $d_0/η$ and the surface tension, i.e. the shear-based Weber… ▽ More We perform interface-resolved simulations of finite-size evaporating droplets in weakly-compressible homogeneous shear turbulence (HST). The study is conducted by varying three dimensionless physical parameters: the initial gas temperature over the critical temperature $T_{g,0}/T_c$, the initial droplet diameter over the Kolmogorov scale $d_0/η$ and the surface tension, i.e. the shear-based Weber number, $We_{\mathcal{S}}$. For the smallest $We_{\mathcal{S}}$, we first discuss the impact on the evaporation rate of the three thermodynamic models employed to evaluate the gas thermophysical properties: a constant property model and two variable-properties approaches where either the gas density or all the gas properties are allowed to vary. Taking this last approach as reference, the model assuming constant gas properties and evaluated with the "1/3" rule, is shown to predict the evaporation rate better than the model where the only variable property is the gas density. Moreover, we observe that the well-known Frössling/Ranz-Marshall correlation underpredicts the Sherwood number at low temperatures, $T_{g,0}/T_c=0.75$. Next, we show that the ratio between the actual evaporation rate in turbulence and the one computed in stagnant condition is always much higher than one for weakly deformable droplets: it decreases with $T_{g,0}/T_c$ without approaching unity at the highest $T_{g,0}/T_c$ considered. This suggests an evaporation enhancement due to turbulence also in conditions typical of combustion applications. Finally, we examine the overall evaporation rate and the local interfacial mass-flux at higher $We_{\mathcal{S}}$, showing a positive correlation between evaporation rate and interfacial curvature, especially at the lowest $T_{g,0}/T_c$. △ Less

Submitted 9 February, 2022; v1 submitted 20 April, 2021; originally announced April 2021.

A fully Eulerian hybrid Immersed Boundary-Phase Field Model for contact line dynamics on complex geometries

Authors: Armin Shahmardi, Marco Edoardo Rosti, Outi Tammisola, Luca Brandt

Abstract: We present a fully Eulerian hybrid immersed-boundary/phase-field model to simulate wetting and contact line motion over any arbitrary geometry. The solid wall is described with a volume-penalisation ghost-cell immersed boundary whereas the interface between the two fluids by a diffuse-interface method. The contact line motion on the complex wall is prescribed via slip velocity in the momentum equa… ▽ More We present a fully Eulerian hybrid immersed-boundary/phase-field model to simulate wetting and contact line motion over any arbitrary geometry. The solid wall is described with a volume-penalisation ghost-cell immersed boundary whereas the interface between the two fluids by a diffuse-interface method. The contact line motion on the complex wall is prescribed via slip velocity in the momentum equation and static/dynamic contact angle condition for the order parameter of the Cahn-Hilliard model. This combination requires accurate computations of the normal and tangential gradients of the scalar order parameter and of the components of the velocity. However, the present algorithm requires the computation of averaging weights and other geometrical variables as a preprocessing step. Several validation tests are reported in the manuscript, together with 2D simulations of a droplet spreading over a sinusoidal wall with different contact angles and slip length and a spherical droplet spreading over a sphere, showing that the proposed algorithm is capable to deal with the three-phase contact line motion over any complex wall. The Eulerian feature of the algorithm facilitates the implementation and provides a straight-forward and potentially highly scalable parallelisation. The employed parallelisation of the underlying Navier-Stokes solver can be efficiently used for the multiphase part as well. The procedure proposed here can be directly employed to impose any types of boundary conditions (Neumann, Dirichlet and mixed) for any field variable evolving over a complex geometry, modelled with an immersed-boundary approach (for instance, modelling deformable biological membranes, red blood cells, solidification, evaporation and boiling, to name a few) △ Less

Submitted 25 May, 2021; v1 submitted 1 April, 2021; originally announced April 2021.

Comments: 52 pages, 19 figures, Submitted to Journal of Computational Physics

Direct numerical simulation of flapping flags in grid-induced turbulence

Authors: Stefano Olivieri, Francesco Viola, Andrea Mazzino, Marco E. Rosti

Abstract: A fully-resolved direct-numerical-simulation (DNS) approach for investigating flexible bodies forced by a turbulent incoming flow is designed to study the flapping motion of a flexible flag at moderate Reynolds number. The incoming turbulent flow is generated by placing a passive grid at the inlet of the numerical domain and the turbulence level of the flow impacting the flag can be controlled by… ▽ More A fully-resolved direct-numerical-simulation (DNS) approach for investigating flexible bodies forced by a turbulent incoming flow is designed to study the flapping motion of a flexible flag at moderate Reynolds number. The incoming turbulent flow is generated by placing a passive grid at the inlet of the numerical domain and the turbulence level of the flow impacting the flag can be controlled by changing its downstream distance from the grid. The computational framework is based on the immersed boundary method for dealing with arbitrary geometries and implemented using a graphics-processing-unit (GPU) accelerated parallelisation to increase the computational efficiency. The grid-induced turbulent flow is first characterised by means of the comparison with well-known results for decaying turbulence and a scale-by-scale analysis. Then, the flag-in-the-wind problem is revisited by exploring the effect of the turbulence intensity on self-sustained flapping. Whilst the latter is still manifesting under strong fluctuations, the main features of the oscillation (including its amplitude and frequency) are altered by turbulence, whose fingerprint can also be qualitatively detected by spectral analysis. Besides their relevance for advancing the fundamental understanding of fluid-structure interaction in turbulence, these findings have potential impact for related applications, e.g., aeroelastic energy harvesting or flow control techniques. △ Less

Submitted 20 July, 2021; v1 submitted 15 March, 2021; originally announced March 2021.

Fiber Tracking Velocimetry for two-point statistics of turbulence

Authors: Stefano Brizzolara, Marco Edoardo Rosti, Stefano Olivieri, Luca Brandt, Markus Holzner, Andrea Mazzino

Abstract: We propose and validate a novel experimental technique to measure two-point statistics of turbulent flows. It consists in spreading rigid fibers in the flow and tracking their position and orientation in time and therefore been named ``Fiber Tracking Velocimetry'' (FTV). By choosing different fiber lengths, i.e. within the inertial or dissipative range of scales, the statistics of turbulence fluct… ▽ More We propose and validate a novel experimental technique to measure two-point statistics of turbulent flows. It consists in spreading rigid fibers in the flow and tracking their position and orientation in time and therefore been named ``Fiber Tracking Velocimetry'' (FTV). By choosing different fiber lengths, i.e. within the inertial or dissipative range of scales, the statistics of turbulence fluctuations at the selected lengthscale can be probed accurately by simply measuring the fiber velocity at its two ends, and projecting it along the transverse-to-the-fiber direction. By means of fully-resolved direct numerical simulations and experiments, we show that these fiber-based transverse velocity increments are statistically equivalent to the (unperturbed) flow transverse velocity increments. Moreover, we show that the turbulent energy dissipation rate can be accurately measured exploiting sufficiently short fibers. The technique has been tested against standard Particle Tracking Velocimetry (PTV) of flow tracers with excellent agreement. Our technique overcomes the well-known problem of PTV to probe two-point statistics reliably because of the fast relative diffusion in turbulence that prevents the mutual distance between particles to remain constant at the lengthscale of interest. This problem, making it difficult to obtain converged statistics for a fixed separation distance, is even more dramatic for natural flows in open domains. A prominent example are oceanic currents, where drifters (i.e.~the tracer-particle counterpart used in field measurements) disperse quickly, but at the same time their number has to be limited to save costs. Inspired by our laboratory experiments, we propose pairs of connected drifters as a viable option to solve the issue. △ Less

Submitted 13 July, 2021; v1 submitted 29 October, 2020; originally announced October 2020.

Comments: 15 pages, 15 figures

Journal ref: Phys. Rev. X 11, 031060 (2021)

Suspensions of deformable particles in Poiseuille flows at finite inertia

Authors: Luigi Filippo Chiara, Marco Edoardo Rosti, Francesco Picano, Luca Brandt

Abstract: We analyze a suspension of deformable particles in a pressure-driven flow. The suspension is composed of neutrally buoyant initially spherical particles and a Newtonian carrier fluid, and the flow is solved by means of direct numerical simulations, using a fully Eulerian method based on a one-continuum formulation. The solid phase is modeled with an incompressible viscous hyperelastic constitutive… ▽ More We analyze a suspension of deformable particles in a pressure-driven flow. The suspension is composed of neutrally buoyant initially spherical particles and a Newtonian carrier fluid, and the flow is solved by means of direct numerical simulations, using a fully Eulerian method based on a one-continuum formulation. The solid phase is modeled with an incompressible viscous hyperelastic constitutive relation, and the flow is characterized by three main dimensionless parameters, namely the solid volume fraction, the Reynolds and capillary numbers. The dependency of the effective viscosity on these three quantities is investigated to study the inertial effects on a suspension of deformable particles. It can be observed that the suspension has a shear-thinning behavior, and the reduction in effective viscosity for high shear rates is emphasized in denser configurations. The separate analysis of the Reynolds and capillary numbers reveal that the effective viscosity depends more on the capillary than on the Reynolds number. In addition, our simulations exhibit a consistent tendency for deformable particles to move towards the center of the channel, where the shear rate is low. This phenomenon is particularly marked for very dilute suspensions, where a whole region near the wall is empty of particles. Furthermore, when the volume fraction is increased this near-wall region is gradually occupied, because of higher mutual particle interactions. Deformability also plays an important role in the process. Indeed, at high capillary numbers, particles are more sensitive to shear rate variations and can modify their shape more easily to accommodate a greater number of particles in the central region of the channel. △ Less

Submitted 20 October, 2020; originally announced October 2020.

The impact of porous walls on the rheology of suspensions

Authors: Marco E. Rosti, Parisa Mirbod, Luca Brandt

Abstract: We study the effect of isotropic porous walls on a plane Couette flow laden with spherical and rigid particles. We perform a parametric study varying the volume fraction between $0$ and $30\%$, the porosity between $0.3$ and $0.9$ and the non-dimensional permeability between $0$ and $7.9 \times 10^{-3}$ We find that the porous walls induce a progressive decrease in the suspension effective viscosi… ▽ More We study the effect of isotropic porous walls on a plane Couette flow laden with spherical and rigid particles. We perform a parametric study varying the volume fraction between $0$ and $30\%$, the porosity between $0.3$ and $0.9$ and the non-dimensional permeability between $0$ and $7.9 \times 10^{-3}$ We find that the porous walls induce a progressive decrease in the suspension effective viscosity as the wall permeability increases. This behavior is explained by the weakening of the wall-blocking effect and by the appearance of a slip velocity at the interface of the porous medium, which reduces the shear rate in the channel. Therefore, particle rotation and the consequent velocity fluctuations in the two phases are dampened, leading to reduced particle interactions and particle stresses. Based on our numerical evidence, we provide a closed set of equations for the suspension viscosity, which can be used to estimate the suspension rheology in the presence of porous walls. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Fluid dynamics of COVID-19 airborne infection suggests urgent data for a scientific design of social distancing

Authors: M. E. Rosti, S. Olivieri, M. Cavaiola, A. Seminara, A. Mazzino

Abstract: The COVID-19 pandemic is largely caused by airborne transmission, a phenomenon that rapidly gained the attention of the scientific community. Social distancing is of paramount importance to limit the spread of the disease, but to design social distancing rules on a scientific basis the process of dispersal of virus-containing respiratory droplets must be understood. Here, we demonstrate that avail… ▽ More The COVID-19 pandemic is largely caused by airborne transmission, a phenomenon that rapidly gained the attention of the scientific community. Social distancing is of paramount importance to limit the spread of the disease, but to design social distancing rules on a scientific basis the process of dispersal of virus-containing respiratory droplets must be understood. Here, we demonstrate that available knowledge is largely inadequate to make predictions on the reach of infectious droplets emitted during a cough and on their infectious potential. We follow the position and evaporation of thousands of respiratory droplets by massive state-of-the-art numerical simulations of the airflow caused by a typical cough. We find that different initial distributions of droplet size taken from literature and different ambient relative humidity lead to opposite conclusions: (1) most vs none of the viral content settles in the first 1-2 m; (2) viruses are carried entirely on dry nuclei vs on liquid droplets; (3) small droplets travel less than 2.5m vs more than 7.5m. We point to two key issues that need to be addressed urgently in order to provide a scientific foundation to social distancing rules: (I1) a careful characterisation of the initial distribution of droplet sizes; (I2) the infectious potential of viruses carried on dry nuclei vs liquid droplets. △ Less

Submitted 16 December, 2020; v1 submitted 2 September, 2020; originally announced September 2020.

Turbulence dictates the fate of virus-containing droplets in violent expiratory events

Authors: Marco E. Rosti, Mattia Cavaiola, Stefano Olivieri, Agnese Seminara, Andrea Mazzino

Abstract: Violent expiratory events, such as coughing and sneezing, are highly nontrivial examples of a two-phase mixture of liquid droplets dispersed into an unsteady turbulent airflow. Understanding the physical mechanisms determining the dispersion and evaporation process of respiratory droplets has recently become a priority given the global emergency caused by the SARS-CoV-2 infection. By means of high… ▽ More Violent expiratory events, such as coughing and sneezing, are highly nontrivial examples of a two-phase mixture of liquid droplets dispersed into an unsteady turbulent airflow. Understanding the physical mechanisms determining the dispersion and evaporation process of respiratory droplets has recently become a priority given the global emergency caused by the SARS-CoV-2 infection. By means of high-resolution direct numerical simulations (DNS) of the expiratory airflow and a comprehensive Lagrangian model for the droplet dynamics, we identify the key role of turbulence on the fate of exhaled droplets. Due to the considerable spread in the initial droplet size, we show that the droplet evaporation time is controlled by the combined effect of turbulence and droplet inertia. This mechanism is clearly highlighted when comparing the DNS results with those obtained using coarse-grained descriptions that are employed in the majority of the current state-of-the-art investigations, resulting in errors up to $100\%$ when the turbulent fluctuations are filtered or completely averaged out. △ Less

Submitted 23 December, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

Journal ref: Phys. Rev. Research 3, 013091 (2021)

Low Reynolds number turbulent flows over elastic walls

Authors: Marco E Rosti, Luca Brandt

Abstract: We study the laminar and turbulent channel flow over a viscous hyper-elastic wall and show that it is possible to sustain an unsteady chaotic turbulent-like flow at any Reynolds number by properly choosing the wall elastic modulus. We propose a physical explanation for this effect by evaluating the shear stress and the turbulent kinetic energy budget in the fluid and elastic layer. We vary the bul… ▽ More We study the laminar and turbulent channel flow over a viscous hyper-elastic wall and show that it is possible to sustain an unsteady chaotic turbulent-like flow at any Reynolds number by properly choosing the wall elastic modulus. We propose a physical explanation for this effect by evaluating the shear stress and the turbulent kinetic energy budget in the fluid and elastic layer. We vary the bulk Reynolds number from 2800 to 10 and identify two distinct mechanisms for turbulence production. At moderate and high Reynolds numbers, turbulent fluctuations activate the wall oscillations, which, in turn, amplify the turbulent Reynolds stresses in the fluid. At very low Reynolds number, the only production term is due to the energy input from the elastic wall, which increases with the wall elasticity. This mechanism may be exploited to passively enhance mixing in microfluidic devices. △ Less

Submitted 4 August, 2020; originally announced August 2020.

Dispersed fibers change the classical energy budget of turbulence via nonlocal transfer

Authors: S. Olivieri, L. Brandt, M. E. Rosti, A. Mazzino

Abstract: The back-reaction of dispersed rigid fibers to turbulence is analyzed by means of a state-of-the-art fully-coupled immersed boundary method. The following universal scenario is identified: turbulence at large scales looses a consistent part of its kinetic energy (via a Darcy friction term), which partially re-appears at small scales where a new range of energy-containing scales does emerge. Large-… ▽ More The back-reaction of dispersed rigid fibers to turbulence is analyzed by means of a state-of-the-art fully-coupled immersed boundary method. The following universal scenario is identified: turbulence at large scales looses a consistent part of its kinetic energy (via a Darcy friction term), which partially re-appears at small scales where a new range of energy-containing scales does emerge. Large-scale mixing is thus depleted in favor of a new mixing mechanism arising at the smallest scales. Anchored fibers cause the same back-reaction to turbulence as moving fibers of large inertia. Our results thus provide a link between two apparently separated realms: the one of porous media and the one of suspension dynamics. △ Less

Submitted 8 July, 2020; originally announced July 2020.

Journal ref: Phys. Rev. Lett. 125, 114501 (2020)

An interface capturing method for liquid-gas flows at low-Mach number

Authors: Federico Dalla Barba, Nicolò Scapin, Marco E. Rosti, Francesco Picano, Luca Brandt

Abstract: Multiphase, compressible and viscous flows are of crucial importance in a wide range of scientific and engineering problems. Despite the large effort paid in the last decades to develop accurate and efficient numerical techniques to address this kind of problems, current models need to be further improved to address realistic applications. In this context, we propose a numerical approach to the si… ▽ More Multiphase, compressible and viscous flows are of crucial importance in a wide range of scientific and engineering problems. Despite the large effort paid in the last decades to develop accurate and efficient numerical techniques to address this kind of problems, current models need to be further improved to address realistic applications. In this context, we propose a numerical approach to the simulation of multiphase, viscous flows where a compressible and an incompressible phase interact in the low-Mach number regime. In this frame, acoustics is neglected but large density variations of the compressible phase can be accounted for as well as heat transfer, convection and diffusion processes. The problem is addressed in a fully Eulerian framework exploiting a low-Mach number asymptotic expansion of the Navier-Stokes equations. A Volume of Fluid approach (VOF) is used to capture the liquid-gas interface, built on top of a massive parallel solver, second order accurate both in time and space. The second-order-pressure term is treated implicitly and the resulting pressure equation is solved with the eigenexpansion method employing a robust and novel formulation. We provide a detailed and complete description of the theoretical approach together with information about the numerical technique and implementation details. Results of benchmarking tests are provided for five different test cases. △ Less

Submitted 24 May, 2020; originally announced May 2020.