Access the full text.

Sign up today, get DeepDyve free for 14 days.

Physics
, Volume 2020 (2005) – May 28, 2020

/lp/arxiv-cornell-university/bubbly-and-buoyant-particle-laden-turbulent-flows-9Z7zedFpGn

- ISSN
- 1947-5454
- eISSN
- ARCH-3341
- DOI
- 10.1146/annurev-conmatphys-031119-050637
- Publisher site
- See Article on Publisher Site

All rights reserved Fluid turbulence is commonly associated with stronger drag, greater heat transfer, and more eﬃcient mixing than in laminar ﬂows. In many natural and industrial settings, turbulent liquid ﬂows contain suspensions of dispersed bubbles and light particles. Recently, much attention has been devoted to understanding the be- havior and underlying physics of such ﬂows by use of both experiments and high- resolution direct numerical simulations. This review summarizes our present understanding of various phenomenological aspects of bubbly and buoyant par- ticleâĂŞladen turbulent ﬂows. We begin by discussing diﬀerent dynamical regimes, including those of crossing trajectories and wake-induced oscillations of rising particles, and regimes in which bubbles and particles preferentially accumulate near walls or within vortical structures. We then address how cer- tain paradigmatic turbulent ﬂows, such as homogeneous isotropic turbulence, channel ﬂow, TaylorâĂŞCouette turbulence, and thermally driven turbulence, are modiﬁed by the presence of these dispersed bubbles and buoyant particles. We end with a list of summary points and future research questions. arXiv:2005.13774v1 [physics.flu-dyn] 28 May 2020 1. INTRODUCTION AND OBJECTIVES OF REVIEW Turbulent multiphase ﬂows are widely prevalent in nature and industry. Typical examples are pollutants dispersed in the atmosphere, air bubbles and plankton in the oceans, sediment laden river ﬂows, and catalytic particles and bubble columns in process technology. In all of these examples, the particles of the dispersed phase have a diﬀerent mass density from the carrier phase. When the dispersed particle is lighter than the carrier ﬂuid element, there can be major consequences on the kinematics and dynamics of both phases, often triggering a multitude of ﬂow modiﬁcations. Many of these ﬂow modiﬁcations can be viewed as originating from the inherent “buoyancy” of the bubbles and light particles, which renders them capable of adding energy over a range of scales while ascending through the turbulent ﬂow. These have decisive roles in many natural phenomena we observe around us, for example in the distributions of buoyant zooplankton undergoing diel vertical migrations (e.g. Calzavarini et al. 2018; Sengupta et al. 2017), in bubble induced mixing of the upper oceans (e.g. Thorpe and Hall 1987), and in engineering applications of drag reduction, heat transfer and mixing (e.g. AlmÃľras et al. 2015; Ceccio 2010; Gvozdić et al. 2018). The subject of dispersed particles in turbulent ﬂows has been studied extensively from the perspective of small (inertial) particles. For reviews on experimental and numerical techniques for the particle laden ﬂows, readers are referred to Crowe et al. (1996); Elghobashi (1994, 2019); Maxey (2017); Prosperetti (2015); Toschi and Bodenschatz (2009); Voth and Soldati (2017). Over the past couple of decades, many insights has been gained on particle dynamics and ﬂow modulations at various scales of turbulence. The vast majority of these have been on neutrally buoyant and heavy particles. As we will discuss in this review, for a variety of reasons, the heavy and neutrally buoyant particle laden ﬂows are experimentally, numerically, and theoretically more amenable as compared to their bubble laden counterparts (Lohse 2018). Balachandar and Eaton (2010) presented the most recent review discussing important aspects of turbulent bubbly ﬂows. Yet, these reviews do not address the subject from the viewpoint of how the “buoyancy” of gas bubbles and light particles can manifest in various forms in turbulent multiphase settings. The aim of the present review is to examine the subject from this perspective, for bubbles and buoyant particles in turbulence. Therefore, all ﬂows considered will be liquid, and we will henceforth use the term buoyant particle, generically, to refer to both gas bubbles and light particles. Vapor bubbles, which have rich underlying physics of their own (Prosperetti 2017), will not be discussed here. The review is organized into two parts. In the ﬁrst part, we will begin by addressing the dynamics of small bubbles and buoyant particles in a regime where particle inertia is dominant, followed by a regime where buoyancy and inertia are competing eﬀects, leading to clustering and reduced rise velocities. We then explore the unsteady wake induced dynamics of ﬁnite sized buoyant particles and air bubbles, and later, their so-called lift and lateral migration tendencies. In the second part of the review, we discuss the main consequences of the particle dynamics on ﬂow modulations including drag reduction and heat transfer enhancement pertaining to diﬀerent turbulent ﬂow systems. The review ends with a summary and outlook towards open issues for future research. 1.1. Dimensionless groups To cover the vastness of issues encountered in buoyant particle laden turbulence requires an exhaustive compilation of the relevant dimensionless parameters; for a review on the general subject of bubbles see Lohse (2018). In the present review, we restrict our attention to a selected subset of these. We divide the control parameter space into particle parameters and ﬂow parameters. For the bubble or buoyant particle, the following dimensionless groups are important: the density ratio of particle (or bubble) to liquid = , the size ratio, d =l, p l p where d is the particle diameter and l a characteristic length scale of the ﬂow. Combining and , one can obtain the Stokes number St = , where is the response time of the particle and a characteristic time p l p l scale of the ﬂow. We note that the problem is suﬃciently deﬁned once two among the three parameters [, , St] are speciﬁed. Further, the eﬀect of buoyancy (or gravity) can be included by using two additional control parameters: the Froude number Fr and the Galileo number Ga. In the most general form1 Fr = U =(gd ) ug ug p gives the ratio of inertia to buoyancy, where U is a relevant velocity scale and g the gravitational acceleration. Ga gd ( 1)= compares buoyancy eﬀects to viscous eﬀects (Veldhuis et al. 2004), where is the 1Note that many diﬀerent, but equivalent, deﬁnitions of Froude number are available. 2 Mathai et al. Bubbles and Buoyant Particles in Turbulence kinematic viscosity of the carrier liquid. Note that the above deﬁnition for Ga, which is alternately referred to as the Archimedes number Ar is equivalent (Ern et al. 2012). Next, when the buoyant particle under consideration is a gas bubble, its deformability can have an inﬂuence on the particle-turbulence interactions. The Eötvös number Eo = gd = (alternately known as the Bond number Bo) gives the ratio of buoyancy force to capillary force. Here is the liquid density, and is the surface tension of the gas-liquid interface. Furthermore, when the buoyant particle’s rotation is important, its mass moment of inertia ratio I I /I becomes a relevant control p l parameter, where I is the particle moment of inertia and I the moment of inertia of the displaced liquid. For p l the turbulence, we will refer to the Reynolds number (with varying deﬁnitions based on the type of the turbulent ﬂow) and the Rayleigh number Ra (for thermal turbulence). With these parameters speciﬁed, we look at a variety of responses of the particles, as well as the ﬂow modulation they induce. Lastly, the volume fraction of the dispersed phase is a crucial control parameter for the system. For particle laden systems, it is common to specify a low volume fraction threshold for transitions to two-way and four-way coupling regimes (Elghobashi 1994). For bubble and light particle laden turbulence, such a criterion is rarely used. As is often the case, even one buoyant bubble can signiﬁcantly modify the turbulence around it, independently of the volume fraction. Thus, to point out a threshold of bubble volume fraction below which they can be considered as passive is not appropriate. In this review we restrict our discussions to situations of low and moderate volume fractions ( < 5%), where complex issues of coalescence and break up (Jha and Govardhan 2015) are not dominant. In addition to the above mentioned control parameters, the particle Reynolds number Re v d = and p T p a buoyancy parameter R v =u , which are estimable output parameters, will also be discussed. Here v v T T is the measured mean particle rise velocity, and u is the rms velocity of liquid ﬂuctuations of the single phase turbulence. Broadly speaking, the Lagrangian dynamics of the particles and the various ﬂow modiﬁcations they induce will be the output quantities of importance. Our intention here is to present these under one umbrella of familiar terminology. Further details of the above mentioned dimensionless groups will be provided in the relevant sections. 2. PARTICLE MOTION IN A FLUID FLOW The ﬂuid motion in an incompressible multiphase ﬂow is governed by the Navier-Stokes (NS) equations @U rp + (Ur)U = + r U + f ; 1. @t where U is the velocity at the location of the particle, p is the pressure, is the liquid density, is the kinematic viscosity of the liquid, and f is a back reaction force per unit mass on the ﬂuid. Based on the ﬂow setting and the modeling approach, f can be either exactly computed (Elghobashi 2019) or modeled (Mazzitelli et al. 2003a,b) or ignored altogether (Toschi and Bodenschatz 2009). Appropriate boundary conditions need to be applied on the boundaries of the ﬂuid domain. For the particle, when employing fully resolved direct numerical simulations (DNS), the boundary condition is either no-slip or free slip or a combination of the two. We begin by considering the equation of motion of a small buoyant spherical particle advected in a ﬂuid ﬂow with velocity U(X(t); t), and in the presence of gravity. Assuming the particle’s spatial dimension is point-like, one can use a form of equations in the spirit of the celebrated Maxey-Riley equation (Maxey and Riley 1983): DU V X = V + F + F + F + F ; 2. p p p p l M B L D Dt where X is the position of the particle (or bubble),V = its volume, and and are the liquid and particle p p l p mass densities, respectively. In obtaining the above equation, the coupling between translation and rotation of the particle has been assumed negligible. The forces contributing on the right-hand-side besides those due to the accelerated ﬂow (which includes the pressure gradient term) are the added mass F , drag F , buoyancy F M D B and a so-called (shear-induced) lift F . Generally, they are modelled as www.annualreviews.org 3 Buoyancy (a) (b) Added mass + Pressure Vortex Drag Shear induced lift Figure 1 Schematic representations of a buoyant spherical particle in a turbulent ﬂow. (a) Schematic of the forces (represented as point forces) felt by a small buoyant particle (or gas bubble) rising past a vortex (a simpliﬁed representation of turbulent eddy). v is the instantaneous particle velocity vector, and the forces due to buoyancy F , drag F , so-called lift F , and the net eﬀect of B D L added mass and pressure F , are shown using colored arrows giving their relative directions. Here F is obtained by MP MP rearranging the added mass and ﬂuid force terms in eq. (2), and accounts for the pressure force on the buoyant particle. For details of the variables shown, we refer to Eq. 3-6. (b) Schematic representation of the buoyant particle rising through a turbulent ﬂow composed of vortices (or eddies) of various length and time scales. Here d is the particle diameter, is the dissipative length scale, and L is the integral length scale of the ﬂow (Pope 2000; Toschi and Bodenschatz 2009). is an intermediate length scale, commonly known as the Taylor micro-scale. DU F = C V X ; 3. M M l p p Dt F = V ( ) g e ; 4. B p p l y F = C V (X U) (r U); 5. L L l p p _ _ F = C jX Uj(X U); 6. D D l p p where C is the added mass coeﬃcient, g is the gravitational acceleration directed along e ^ , the unit-vector in the M y vertical direction. The shear induced lift depends on the alignment between the vorticity vector (r U) and the relative velocity of the particle (X U), and C and C the constants of proportionality for the lift and drag, p L D respectively (life and drag coeﬃcents). It is important to ﬁrst appreciate the directions of these force vectors for a rising buoyant particle. A simpliﬁed picture is drawn in Figure 1(a), through a cartoon of the buoyant particle rising along the downward side of a “vortex”. The approximate directions of the eﬀective forces acting on the particle are also sketched: the drag (purple arrow) is oriented to oppose the direction of the particle’s velocity, the buoyancy (green arrow) is vertical, the shear-induced lift force (red arrow) acts perpendicular to the plane containing the particle velocity and the vorticity vector, and is opposed by the centrifugal force (blue arrow) that is directed towards the eye of the vortex, as the particle is lighter than the surrounding liquid. Note that the Faxen ´ forces (accounting for ﬂow non-uniformity at the particle scale) and the Basset history force have been omitted for simplicity (Auton et al. 1988; Rensen et al. 2001). The relative importance of these terms for buoyant rising particles is still to be resolved (Calzavarini et al. 2009, 2012; Homann and Bec 2010). When the condition of particle Reynolds number Re 1 is met, the drag coeﬃcient reduces to C = p D 24=Re , which implies the linear drag relation F =3 d (X U). While the basic form of this equation p D p p is founded on a uniﬁed treatment of particles, drops and bubbles, it further assumes that for Re 1, the particle locally sees a Stokes ﬂow (Leal 1980; Maxey and Riley 1983) despite the unsteadiness and turbulence in the 4 Mathai et al. Bubbles and Buoyant Particles in Turbulence carrier ﬂow. The expression for the drag used above assumes a contaminated air-liquid interface for the bubble. For a clean bubble interface, the prefactor of F is slightly modiﬁed (Mougin and Magnaudet 2001), although retaining its functional form. The assumption of a contaminated interface is indeed reasonable in most natural and industrial ﬂows, since the carrier liquid is almost never ultrapure. Assuming potential ﬂow in the outer regions, we can use C = 1=2. Similarly, by considering the momentum ﬂux far from the particle, one can obtain the Auton lift C = 1=2 (Auton et al. 1988; C. R. Hunt et al. 1994), which applies to small spherical bubbles or particles in a shear ﬂow. For the purpose of simpliﬁcation, we will now consider the turbulence to be homogeneous and isotropic. A state of homogeneous isotropic turbulence (HIT) is fully determined by knowing the kinematic viscosity , an outer length scale, and a time scale. For a bubble or light particle rising through HIT (see Figure 1(b)), it sees the largest and smallest length scales, which are the energy-injection scale L and the Kolmogorov (or dissipative) scale , respectively. Additionally, represents an intermediate scale known as the Taylor micro-scale, upon which the Taylor Reynolds number Re u = is based. Since the particle size is comparable to the dissipative scale, it is appropriate to non-dimensionalize the equation of particle motion using the Kolmogorov units of length () and time ( ). This leads to Du 1 1 x = + (u x _ ) + e ^ + (u x _ ) !; 7. p p y p Dt St Fr 3 where the lower bold case variables x , u, and ! denote the new dimensionless vectors for particle position, ﬂow velocity and vorticity, respectively. Here 3=(1 + 2) is an eﬀective density ratio that takes the ﬂuid added mass into account. The Stokes number St and the Froude number Fr are deﬁned generically 12 ( 1)g to be valid for light (1 < 3), heavy (0 < 1) and neutrally buoyant ( = 1) particles. Here a = is the acceleration at the Kolmogorov scale. 3. BUOYANT PARTICLE DYNAMICS 3.1. Regimes of dynamics A wealth of intriguing phenomena have been reported for buoyant particle suspensions in turbulent ﬂows. Figure 2 outlines a simpliﬁed regime diagram for a buoyant spherical particle (or bubble) when the particle parameters St, Fr, and Re are varied, with , , and background turbulent ﬂow unchanged. When extremely small, a bubble (or light particle) can serve as a passive ﬂow visualization tool (see FT regime in Figure 2), while in other scenarios, they have been used to selectively sample intense vortical structures present in ﬂows (Douady et al. 1991; La Porta et al. 2001). In the case of ﬁnite particle inertia (St O(1) and Fr 1), it is widely known that bubbles cluster in high vorticity regions (see IC regime in Figure 2). Wang and Maxey (1993), Spelt and Biesheuvel (1997), and Mazzitelli et al. (2003a) observed these eﬀects in their simulations of homogeneous isotropic turbulence laden with microbubbles. The phenomenon was later extensively explored by others (Calzavarini et al. 2008a,b), thus providing a uniﬁed understanding of the clustering behavior of light, heavy and neutrally buoyant particles. The methods adopted to quantify clustering have been diverse, with Calzavarini et al. (2008b) using the Kaplan- Yorke dimension and Minkowski functionals, while others have used Voronoï tessellations (Fiabane et al. 2012; Monchaux et al. 2010; Obligado et al. 2014; Tagawa et al. 2012) or simply the relative particle concentration in high vorticity ﬂow regions (Mazzitelli et al. 2003b), but all to the same end. An assessment of the actual forces that bubbles are subject to in turbulent environments was performed by Volk et al. (2008a,b), who experimentally investigated the acceleration dynamics of small bubbles (d 75 m) in a relatively intense turbulent ﬂow (Re = 850) generated by a von Kármán ﬂow apparatus. The high intensity of turbulence ensured that the 75 m air bubble had a St= 1:85, with a Froude number Fr 10, which meant that the role of buoyancy was negligible. These bubbles showed an acceleration variance four times that of the ﬂuid, which can be presumed to be a combined eﬀect of inertial forces and preferential accumulation in high vorticity regions of the ﬂow (Calzavarini et al. 2009). For reviews on inertial particle dynamics in the absence of buoyancy, readers are directed to Toschi and Bodenschatz (2009); Voth and Soldati (2017). It is worthwhile to take a step back to appreciate the various regime transitions seen by an air bubble ( = 10 ) when its size is increased, while keeping all other parameters ﬁxed. The * symbols in Figure 2 www.annualreviews.org 5 FT - Fluid tracer-like IC - Inertial Clustering RR - Reduced Rise-velocity CT - Crossing trajectories CT RR IC Buoyancy Turbulence Figure 2 A simpliﬁed regime diagram for a buoyant particle (or air bubble) in a turbulent ﬂow. The dynamics of the particle can be summarized as a function of of three dimensionless parameters, namely the Stokes number, St - a measure of particle inertia, the Froude number, Fr - the ratio between turbulence strength and gravity, and the particle Reynolds number, Re (Calzavarini et al. 2009; Mathai et al. 2016a). The acronyms FT, IC, RR, and CT denote the diﬀerent dynamical regimes, which have been expanded at the top of the ﬁgure. The horizontally hatched (magenta) region is the regime where the eﬀect of crossing trajectories can be seen to various extents depending on the value of St. The symbols (connected by arrows) help to picturize a gradual regime transition seen by an air bubble (density ratio = 10 ) upon increasing its size at constant turbulence level, and assuming that the bubble Reynolds number Re 1. Note that the eﬀects of changing and size ratio () are not represented in this diagram. show this transition if the bubble Re is small. However, in practicality for a laboratory scale turbulent water ﬂow (say Re O(100) and L 100 mm) the actual regime transitions can be more complex. For instance, a microbubble with d < 10 m can still be considered a good tracer of the turbulent ﬂow, while upon increasing its diameter d to around100 m, the eﬀects of buoyancy begin to play a role in the dynamics. When the bubble is a few hundred microns in size, in addition to the buoyancy, the bubble inertia becomes important, and upon increasing the diameter further (d > 500m), one can expect noticeable non-Stokesian and ﬁnite Re p p contributions. In the following sections, we will discuss these regimes of dynamics in more detail. 3.2. Crossing trajectories The dynamics of a particle advected in HIT in a regime where Fr < 1 will be addressed in this section. ( 1)g The importance of buoyancy (Fr < 1) naturally implies that the particle experiences a mean vertical drift through the turbulent ﬂow, and hence the name “crossing trajectories” (see Figure 1(b)). 3.2.1. Non-inertial particles with buoyancy. We begin with bubbles and buoyant particles that are in the non- inertial limit, i.e. St 1. Since a very small St naturally implies a tiny particle dimension, such bubbles are commonly used as tracers in turbulence experiments. Recently, Mathai et al. (2016a) conducted a combined experimental and numerical study on the dynamics of such small bubbles and particles in the non-inertial (St 1) 6 Mathai et al. Bubbles and Buoyant Particles in Turbulence Turbulent wake Unsteady wake driven Stokesian dynamics FT Non-inertial Inertial limit. An interesting consequence of buoyancy is that even small St bubbles are subject to intense accelerations. For Fr 1 and St 1, Eq. (1) is dictated by the balance between just the drag and buoyancy terms. This yields Du St an expression x ' + @ u, where @ u are the gradients of the turbulent ﬂow velocity at the particle’s p y y Dt Fr location. The buoyancy parameter can be exactly computed as R = u St=(u Fr). The acceleration variance th 2 (i component) of the buoyant particle a can be expressed as an enhancement over the ﬂuid acceleration variance a : p St D E ' 1 + ; 8. 2 Fr where = 2=(15a ) for the horizontal component, and = 1=(15a ) for the vertical component, with a x 0 y 0 0 the so-called Heisenberg-Yaglom constant (Voth et al. 2002). These relations follow exactly from the assumption that the turbulence is statistically isotropic (Pope 2000). In other words, the increase in acceleration variance is a direct consequence of the vertical drift of the bubble through the turbulent eddies. While the acceleration variance increases, its decorrelation time compared to the ﬂuid is suppressed (see Figure 3(a)-(b)), as the drifting particle spends comparatively less time within the turbulent eddies. Similarly, the kurtosis of acceleration is diminished, a consequence of the fact that the spatial velocity gradients @ u in turbulence are less intermittent than the ﬂuid acceleration (for further details, see Mathai et al. (2016a)). Finally, we note that in this non-inertial limit, the behavior of buoyant particles is expected to be nearly identical to that of heavy particles, for ﬁxed St/Fr (Csanady 1963; Maxey 1987; Parishani et al. 2015). Of course, the eﬀect is more pronounced for bubbles as their density contrast is generally larger as compared to that of the heavy particles in liquid ﬂows. Thus, a tiny bubble or droplet is not necessarily a good tracer of turbulent acceleration. In reality, the situation of ﬁnite g=a is common for bubbles that drift through the turbulent oceans (g=a 100 1000), and for droplets settling through clouds (g=a 10 100) (Bodenschatz et al. 2010). On the practical side, these also point to the key condition St/Fr 1 that must be met (in addition to St 1) for the usage of small bubbles (or droplets) for ﬂow visualization and particle tracking in turbulent ﬂows. 3.2.2. Inertia and buoyancy. Next, we consider the cases when Fr 1 and St 1, i.e. a regime where both buoyancy and inertia are important, and in competition. Employing an Euler-Lagrangian point particle (PP) approach of a form similar to eq (7), Mazzitelli et al. (2003a) studied the eﬀect of increasing St at ﬁxed Fr for bubbles rising in isotropic turbulence (Re = 62). Snapshots of the bubble distribution obtained from their simulations are shown in Figure 3(c)-(d). At St = 0.33, the particles (bubbles) are nearly uniformly distributed, but at St = 1 clustering is visibly ampliﬁed. Three-dimensional Voronoï analysis was used to quantify the clustering of bubbles in homogeneous isotropic turbulence using data sets from numerics in the point particle limit (without gravity and lift) and an experimental data (Tagawa et al. 2012). In the PDFs of normalized Voronoï volumes (V=V) shown in ﬁgure 3(e), the solid black curve - a -distribution - is representative of non-clustering (randomly distributed) particles. For bubbles, they observed that the probability of ﬁnding both small and large Voronoï volumes V=V is generally higher. The two regions of small and large volumes can be used to identify clusters and voids. A high probability for low values of V=V is signatory of intense clustering. As shown in the ﬁgure, when St increases, the probability of ﬁnding clusters (and voids) increases, and reaches a maximum value at St 1.6, suggesting that the strongest clustering for bubbles occurs in the St range 1-2. A direct quantiﬁcation of the clustering in the presence of buoyancy is shown in Figure 3(f), where the number fraction of bubbles located in a vorticity dominated region is given by the green symbols. Since the buoyancy and lift terms were included by Mazzitelli et al. (2003b), the relative degree of clustering is typically less than what is seen in minimalistic simulations (Calzavarini et al. 2008a) where these terms were neglected. 3.2.3. Non-Stokesian bubbles and particles. For the more commonly encountered situation of air bubbles in water ﬂows, a Stokes number of order one typically almost never satisﬁes the condition of Re 1 (Magnaudet and Eames 2000; Mathai et al. 2018a). Hence a modiﬁed consideration of the drag, added mass and lift forces is essential to predict the trajectory of ﬁnite sized and ﬁnite Reynolds number bubbles. As a model problem, Sridhar and Katz (1995) experimentally studied the entrainment of such air bubbles (d 500 800 m) by a vortex ring. They reported that the drag coeﬃcient was comparable to the Schiller-Naumann parameterization www.annualreviews.org 7 25 (a) (c) (d) St = 0.33 St = 1 1.2 (b) (e) (f) Experiment St = 0.04 0.8 Numerics St = 0.1 0.4 St = 0.6 St = 1.6 St = 4.1 0 10 20 0.1 0.4 0.7 1.0 St Figure 3 Accelerations and clustering behavior of small ( < 1) buoyant particles and bubbles in isotropic turbulence obtained using Euler-Lagrangian DNS with the inclusion of buoyancy and lift. (a) Normalized acceleration variance of buoyant non-inertial (St 1) particles in turbulence versus St/Fr at Re = 75. The acceleration variance here is normalized by the acceleration variance of the ﬂuid. (b) Normalized acceleration decorrelation time (time to reach an autocorrelation value of 0.5) of the particles for the same cases as in (a). The normalization here is with the Kolmogorov times scale of the turbulence. Solid and hollow symbols in (a) and (b) correspond to vertical (gravity) and horizontal components, respectively. The solid and dashed curves in (a) and (b) are theoretical predictions. Figures (a,b) adapted from Mathai et al. (2016a). (c) and (d) Projections of bubble distribution in isotropic turbulence (Re = 62) for two Stokes numbers St = 0.33 and St = 1, respectively. Note that the particle concentration is identical for (c) and (d). Bubble clustering is pronounced in the St = 1 case (d). We estimated the ratio St/Fr ( v =u ) to be 1.0 and 3.0, respectively, in (c) and (d), which implies R < 1. Figures T v adapted from Mazzitelli et al. (2003a). (e) Clustering of light particles in turbulence, quantiﬁed using probability density functions (PDFs) of normalized Voronoï volumes for diﬀerent values of St. The solid curve is a -distribution, which is representative for non-clustering (randomly distributed) particles. Green data points are based on DNS using an Euler-Lagrangian (point particle) model, at Re = 180. Black data points are experimental data for microbubbles (St = 0.04) in turbulence at Re = 162. The results suggest that the maximum clustering for bubbles occurs at St around 1-2. This ﬁgure was adapted from Tagawa et al. (2012). (f) Number of bubbles located in a particular zone of ﬂuid n , normalized by the volume fraction of that zone v , versus St. Here again Re = 62, and St/Fr was kept ﬁxed at 1. The symbol colors refer to zones: eddy (green), shear (orange), streaming (magenta) and convergence (blue). This ﬁgure is adapted from Mazzitelli et al. (2003b). (see Figure 4(a)), while the lift forces did not agree with the existing theoretical or numerical models (Maxey 1987; Maxey and Riley 1983; Tio et al. 1993). More recently, Aliseda and Lasheras (2011) experimentally investigated the behavior of small spherical bubbles (100 – 1000 m) immersed in a homogeneous isotropic turbulent water ﬂow (see Figure 4(b)). Within the turbulent ﬂow the concentration ﬁeld of the bubbles was altered, with preferential accumulation at the smaller scales and reduced rise velocities as compared to the value in stationary liquid. These can be interpreted as occurring due to two phenomena. Firstly, the pressure ﬂuctuations drive the inertial bubbles to the cores of the vortices. Secondly, the lift forces cause the bubbles to be preferentially transported toward downﬂow regions, where, in combination with an increased relative velocity (increased viscous drag), they are further slowed down (see cartoon in Figure 1(a)) (Mazzitelli et al. 2003b). The conditions of small size (d = 1) and nearly spherical shape is typically satisﬁed only for sub- millimeter air bubbles (Aliseda and Lasheras 2011; Sridhar and Katz 1995). However, the vast majority of bubble laden turbulence operates under conditions where the bubbles are of ﬁnite size, free to move, and, most 8 Mathai et al. Bubbles and Buoyant Particles in Turbulence 4 (a) (b) Experiment 0.12 3 Sridhar & Katz (1995) Schiller-Naumann 0.08 0.04 0 20 40 60 80 100 0 100 200 300 400 500 Figure 4 Drag coeﬃcient and rise velocity modiﬁcations for sub-millimetric air bubbles rising in turbulent ﬂows. (a) Comparison of the drag coeﬃcient for 500 800 m diameter air bubbles (in a vortical ﬂow), along with the predictions of the Schiller-Naumann drag model. Data reproduced from Sridhar and Katz (1995). (b) Direct measurements of rise velocity v for small bubbles versus bubble diameter d , in a grid generated turbulent ﬂow at Re 404 (as roughly estimated by us). For comparison the following are given. v : the terminal velocity estimated theoretically using Stokes drag theor C = 24=Re , v : the terminal velocity measured in triply distilled water, and v : the same quantity measured for tap D p dist tap water. Data adapted from Aliseda and Lasheras (2011). importantly, deformable (Deane and Stokes 2002; Duineveld 1995). These call for a more detailed consideration of the coupled interaction between the bubble topology and its wake-induced dynamics arising from the ﬁnite Re and ﬁnite Weber number eﬀects (Ryskin and Leal 1984). A fully resolved treatment of the interaction between isotropic turbulence and large solid spherical particles has been performed in a variety of conﬁgurations (Chouippe and Uhlmann 2015; Homann and Bec 2010; Naso and Prosperetti 2010). In comparison, DNS of turbulent bubbly ﬂows are challenging, owing to the interface deformations and internal circulations, along with the need to resolve a wide range of length and times scales inherent to the turbulent ﬂow. Loisy and Naso (2017) studied an isolated deformable bubble freely rising in an otherwise isotropic turbulent ﬂow using DNS (see Figure 5(a)). The buoyancy parameter (R 2 [0:63; 2:17]) and bubble Reynolds number (Re 2 [17, 62]) were both moderate, v p and hence in quiescent liquid the bubble rises along straight vertical paths. However, with decreasing R the trajectories become erratic and increasingly deviate from vertical paths (see Figure 5(b)-(d)). This was accompanied by a reduction in rise velocities. With regard to the statistics of bubble velocity and acceleration, the probability density functions (PDFs) was nearly Gaussian for the velocity and showed stretched tails for accelerations. Lastly, the bubble showed a preference for increased residence in vorticity-dominated regions, 2 3 here, computed using the discriminant D = 27R + 4Q , where Q and R are the second and third invariants of the velocity gradient tensor, respectively (Naso and Pumir 2005). This was revealed through conditional sampling of the average discriminant proﬁle and discriminant ﬁeld around the bubble (see Figure 5(e) and Figure 5(f)-(h), respectively). The analyses of Loisy and Naso (2017) (given above) show that the dynamics of a moderate Re bubble in turbulence (rise velocity, PDF shapes and clustering) are, at least, qualitatively captured by the point-particle model (Calzavarini et al. 2009), despite their ﬁnite Re , ﬁnite size, and deformed shape. However, this is not reﬂective of the full picture, since certain aspects of the bubble statistics are markedly diﬀerent from the PP predictions. The PDF of the longitudinal acceleration, i.e. the component of bubble acceleration directed along its instantaneous velocity, was found to be negatively skewed, a feature not captured by the PP model even with the inclusion of a back reaction force. Whether the origin of this lies truly in the time-irreversibility (Xu et al. 2014) of turbulence (as postulated by Loisy and Naso (2017)), or in the asymmetry of ﬂow induced forcing on the bubble, remains to be unravelled. www.annualreviews.org 9 Bubble rise velocity (a) (b) (c) (d) -1 0 1 2 3 -2 0 2 4 -4 -2 0 2 4 6 0.5 (e) (f) (g) (h) 0.6 1 1 0.4 0.4 0 0 0 0.3 0.2 -1 0.2 -1 -1 0.1 0 1 2 3 4 5 -2 0 1 2 0 1 2 0 1 2 Figure 5 Direct numerical simulations of a slightly deformed bubble ( 10, Re 2 [17; 62]) rising in a turbulent ﬂow at Re = 30. (a) Snapshot of the (normalized) vertical velocity ﬁeld of the turbulent ﬂow on a plane passing through the bubble center. (b)-(d) Sample bubble trajectories for decreasing buoyancy parameter R v =u . R = 2:17 in (b), 1.11 in (c) and 0.63 in v v (d). The bubble trajectories are colored by their instantaneous vertical velocity normalized by their terminal velocity. (e) Radially averaged proﬁle of discriminant D around the bubble for the three diﬀerent R cases. The abscissa =d is the v p (normalized) radial distance from the centroid of the bubble. The sign of the discriminant D (ordinate) can be used to highlight the vorticity-dominated regions. (f)-(h) Average discriminant ﬁeld around the bubble for decreasing R (left to right). Color scale centered around the mean value D 3 to show that values above D are the vorticity dominated regions. The discriminant ﬁeld clearly demonstrates that clustering in high vorticity regions is increased as the bubble rise velocity decreases. Data and ﬁgures adapted from Loisy and Naso (2017). 3.3. Wake-driven dynamics and path-instabilities With increasing buoyancy over inertia, the bubble or buoyant particle’s Reynolds number Re can reach a few hundreds. Two important changes come into eﬀect in this regime. Firstly, the mean drag coeﬃcient loses its Reynolds number dependence, and C becomes weakly dependent on Re . In addition, such bubbles and particles D p also experience ﬂuctuating components of forces that originate from the instability of their wakes. Although the mean forces on the particle can still be approximated, the instantaneous drag and lift can no longer be described using simpliﬁed coeﬃcients. Furthermore, owing to the lightness of the particle, this regime paves way for a strong coupling between the unsteady wake dynamics and the particle motion, often resulting in vigorous path instabilities (Brücker 1999; Ern et al. 2012; Mathai et al. 2016b, 2017, 2018b; Mougin and Magnaudet 2006). As reported in Risso (2017), there is, today, compelling evidence that the wakes and dynamics of isolated buoyant particles are remarkably robust to turbulent perturbations (Ford and Loth 1998). Therefore, the forces acting on an isolated buoyant particle in a ﬂow can still provide a basis for understanding dispersed two phase ﬂows. 3.3.1. Wake-driven bubbles. Building on the original formulations proposed by Kirchhoﬀ (Galper and Miloh 1995; Lamb 1993) for a rigid particle in a arbitrary irrotational ﬂow, Mougin and Magnaudet (2001) extended the case to a situation where the NS equations governing the liquid ﬂow (Auguste and Magnaudet 2018; Mougin and Magnaudet 2002) are considered in conjunction with Newton’s laws for the translational and rotational 10 Mathai et al. Bubbles and Buoyant Particles in Turbulence (b) (a) (d) Diffusive regime Ballistic regime 10 mm Experiment Prediction (c) (long-time MSD) Figure 6 The horizontal component of velocity and dispersion for 1.8 mm diameter air bubbles (R 10) rising in active-grid-generated turbulence at Re = 110. Here the bubble size ratio 5. (a) A sample image from one of the cameras. (b) Lagrangian temporal autocorrelation function of the horizontal component of the bubble velocity. (c) Kurtosis of the horizontal velocity-increment PDFs as a function of time lag . (d) Horizontal component of the mean squared displacement (MSD) for the bubble as a function of time lag . The well-known ballistic-diﬀusive behavior of ﬂuid tracers is shown for comparison: dashed gray line (/ ), and solid gray line (/ 2T ). The dotted blue line on the right shows the prediction for the reduced long-time dispersion for the bubble, obtained using a crossing-trajectory hypothesis. Data and ﬁgures adapted from Mathai et al. (2018a). dynamics of a “ﬁxed-in-shape” non-spherical bubble or buoyant particle. The Kelvin-Kirchhoﬀ equations, which disregard small-scale deformability eﬀects, demonstrated that wake instability and anisotropic added-mass eﬀects of oblate spheroids (Cano-Lozano et al. 2016; Mougin and Magnaudet 2006) are indeed suﬃcient to explain the experimentally observed path instabilities (Wu and Gharib 2002) of millimetric bubbles. The method, now at our ﬁngertips, is yet to be used widely for turbulent bubble laden ﬂows, and holds the potential to yield useful insights at a reduced computational cost, since the condition of a continuously deforming bubble interface is relaxed. Experimentally, the dynamics of isolated millimetric air bubbles in turbulent ﬂows was studied by Mathai et al. (2018a). In this work, the authors generated a suspension of 1.8 mm diameter air bubbles in an active grid generated homogeneous isotropic turbulent water ﬂow with Re 2 [110; 300]. The bubbles were nearly spherical (see Figure 6(a)), and their volume fraction in the experiments was low ( 5 10 ). At a low level of turbulence (Re = 110), the Lagrangian temporal autocorrelation of the bubble velocity and its kurtosis both showed periodicity (see Figure 6(b)-(c)), clearly indicative of wake shedding at a frequency f 0:1v =d . The viv p p eﬀects of these on the spreading of the bubbles was analyzed using the mean-squared-displacement (MSD), which was then compared to the well-known regimes of Taylor-dispersion for the ﬂuid in turbulence (see Figure 6(d)). At short times, the bubble MSD grows ballistically (/ ), whereas at a larger time scale set by the wake-shedding frequency f , it approaches the diﬀusive regime where the MSD / . Note that for the bubbles, the ballistic viv 02 2 regime lies well above the u prediction of ﬂuid tracers, while the diﬀusive regime MSD for the bubbles lies well below the 2u T of ﬂuid tracers. Here T is the Lagrangian integral time scale which sets the ballistic- L L diﬀusive transition time for ﬂuid tracers in turbulence. To conclude, with high Reynolds number millimetric air bubbles in turbulence, we can appreciate an elegant merger of two classical phenomena, namely the wake-induced velocity ﬂuctuations of the bubbles (at short times), and the reduced dispersion (at longer times) originating from the crossing trajectories eﬀect (Calzavarini et al. 2018; Mathai et al. 2016a, 2018a; Mazzitelli and Lohse 2004). The motion of an even larger bubble (d 9 mm) in turbulence was considered by Ravelet et al. (2011). At these sizes in a turbulent water ﬂow, the bubble shows considerable deformability (We ' 11:6), and its Reynolds number Re 2800. The researchers tracked the bubble motion and orientation in intricate detail using three-dimensional shape recognition, yielding statistics of bubble translation, rotation and deformation in the turbulence (see Figure 7(a)-(d)). The bubble dynamics was found to be governed by three fairly independent mechanisms. The average bubble shape is imposed by the mean motion of the bubble relative to the liquid. Further to this, wake instability of the bubble generates periodic oscillations in its velocity and orientation. Lastly, www.annualreviews.org 11 -135 -135 (a) (b) (e) 100 1 1 1 1 1 -145 -145 -50 -100 -155 -150 -155 -200 -10 -5 0 5 10 -10 -5 0 5 10 -20 0 0 X (mm) Y (mm) -20 (c) (d) (f) -153 -148 -143 -10 -5 X (mm) 10 0 -10 -10 Y (mm) -10 -5 0 5 10 X (mm) Figure 7 Kinematics of a deformable air bubble (d 9:3 mm, Re = 2800) rising in a (rotating) grid-generated turbulent ﬂow, based p p on Ravelet et al. (2011). The buoyancy parameter R lies in the range [5, 10]. (a)-(d) Illustration of image processing steps on the bubble images, with equivalent ellipse (red), center of mass (red circle), major axis (blue) and minor (green) axis of the best-ﬁt ellipsoid. (e) Standard deviation of the liquid velocity measured in one vertical plane and eight horizontal planes, superimposed with a trajectory of the bubble (black). Note that the bubble here is trapped in a region of highest vorticity. (f) Normalized power spectra of horizontal velocity v (black), orientation (red), and major axis a (blue). The spectra are normalized by the frequency f = 15 Hz, as done by Ravelet et al. (2011). the turbulence adds to the random deformations, which under rare circumstances can even lead to bubble breakup. The bubbles were observed to be trapped inside a vortex at the core of the ﬂow (see Figure 7(e)). The temporal spectra of horizontal velocity, orientation and semi-axis lengths (Figure 7(f)) reveal peaks at around 8.5 Hz. This is consistent with a wake instability and yields a Strouhal number St 0:27. 3.3.2. Wake-driven rigid buoyant particles. In the following discussion we provide an overview of studies on large buoyant particles of ﬁxed shape advected in turbulent ﬂows. Initial headway was made using ﬁnite sized neutrally buoyant spherical particles (Bellani and Variano 2012; Zimmermann et al. 2011a,b). From a modeling perspective, the point-particle approach with the so-called Faxen ´ corrections (Faxén and Hilding 1922; Maxey and Riley 1983) for particle size was adopted widely to study neutrally buoyant spherical particles (Calzavarini et al. 2009). Still, the question of ﬁnite slip velocities, which occur in most practical situations, raises important concerns about the validity of this treatment (Prosperetti 2015). By performing experiments using marginally buoyant ( 1:05) rigid ﬁnite sized spheres in turbulence, Mathai et al. (2015) showed that even a slight density mismatch is suﬃcient to cause the dynamics of buoyant particles to deviate signiﬁcantly from the Faxen ´ model predictions. The deviations increase dramatically upon reducing the density ratio, until for very buoyant particles with 2:90, i.e. comparable to the of a bubble, vigorous path oscillations outweigh the turbulence induced motions. Figure 8(a)-(c) demonstrate this eﬀect of increasing buoyancy (Ga from 30 to 3000), whereupon the turbulence induced chaotic dynamics (in (a)) are replaced by wake induced oscillations (in (b) and (c)) reminiscent of Lissajous orbits (Govardhan and Williamson 2005). Three mechanisms contribute to the increasing path-oscillations witnessed here. Firstly, a Ga increase makes the wake-induced forces stronger. Next, an accompanying increase in the buoyancy parameter R v =u occurs. When R increases, the rising v T v spherical particle crosses the turbulent eddies at increasingly higher speeds, thus having little time to respond to the turbulent ﬂuctuations. A third, not so obvious, inﬂuence was revealed experimentally in a recent study by Mathai et al. (2018b), which hints that the observed path-instabilities are augmented by the particle’s rotational motions 12 Mathai et al. Bubbles and Buoyant Particles in Turbulence Z (mm) Y (mm) Z (mm) Z (mm) PSD Z (mm) (a) (b) (c) Ga = 3000 Ga = 30 Ga = 900 (d) (e) (f) (d) (d) (d) (d) (d) (d) (d) (d) (d) (d) (d) Frenet–Serret -1000 -500 0 500 1000 Coordinates Figure 8 Kinematics and dynamics of rigid buoyant spheres ( 10) in a turbulent ﬂow with Re 300. (a)-(c) Three dimensional trajectories of buoyant spheres for increasing Galileo number, Ga. In obtaining these trajectories, the mean rising motion of the spheres were counteracted by matching the sphere rise velocity with the mean downward ﬂow speed. Note that the Ga increase (in (a)-(c)) coincides with an increase in the buoyancy parameter R , which was measured to be around 1.7 in (a), 13.6 in (b) and 45.6 in (c). (d) The Frenet-Serret coordinates: TNB, used to decompose the instantaneous motion of the sphere in turbulence. (e) PDFs of the centripetal acceleration a for two spheres with identical properties ( = 0:89; 100, Ga = 6000), except their mass moments of inertia I = 1.0 (black circles) and 0.6 (red circles). Here a is directed along the N vector in the Frenet-Serret formulae. (e) Reducing I triggers a dramatic increase in the translational and rotational accelerations (see (e) and (f)). Drawings and ﬁgures are adapted from Mathai et al. (2015, 2018b). as well. To explain this, we revisit the Kelvin-Kirchhoﬀ equations expressing linear and angular momentum conservation for a buoyant spherical particle: 1 dV F p Q ( + + B ) + V = + ( 1)g^ e ; 9. 1 p p j 2 dt V l p 1 d p Q ( I + B ) = ; 10. 10 dt V d p p where V is the sphere velocity vector, is the sphere angular velocity vector, g is the acceleration due to p p gravity, I I =I is the moment of inertia ratio, with I the sphere moment of inertia, and I the moment of p l p l inertia of the liquid volume displaced by the spherical particle. F and T represent the “quasi-static” ﬂuid Q Q force and torque vectors, respectively, which result from the existence of vorticity in the ﬂow. These terms can be straightforwardly obtained by integrating the local stress and moment over the sphere surface. Note that = is the dimensionless Stokes boundary layer which develops in a time . The prefactors B = 18 and v 1 B = 2 are analytically obtained from unsteady viscous contributions (Auguste and Magnaudet 2018; Zhang and Stone 1998). Equations (9) & (10) help appreciate the strong coupling that could exist between translation and rotation for a buoyant spherical particle (I < 1). Assuming the time available for the Stokes layer to www.annualreviews.org 13 PDF PDF Prolate Oblate 1.0 0.75 1.5 2.0 0.50 4.0 0.25 Aspect Ratio Figure 9 Trajectories of buoyant ellipsoids of diﬀerent aspect ratios rising in quiescent liquid. The trajectories seen are remarkably robust to background turbulent ﬂuctuations. Here, = 0:8, and Ga = 300 calculated based on the diameter 1=3 d = (6V =) . is used to distinguish oblate ( < 1) and prolate ( > 1) ellipsoids from a spherical particle ( = 1). pe p The oblate ellipsoids provides an approximation of the mean deformed shapes of millimetric air bubbles in water (Ravelet et al. 2011; Zenit and Magnaudet 2008). Figures based on Will et al. (2019) (unpublished). develop scales with the vortex shedding time scale, Mathai et al. (2018b) estimated that up to a moderate Ga, the role of I ought to be insigniﬁcant. Whereas, upon increasing Ga further, I becomes increasingly dominant in Eq. (10). 3.3.3. Rotation induced accelerations. To analyse the role of rotation further, it is advantageous to adopt a Lagrangian frame of reference that is oriented with respect to the sphere’s instantaneous motion (see Figure 8(d)). The mutually orthogonal Frenet-Serret coordinates are the tangent T = x_ =jx_ j, normal N = B T, and p p binormal B = (x _ x )=kx _ x _ k vectors, which are directed along the particle velocity, curvature, and p p p p a direction perpendicular to the trajectory plane, respectively. For neutrally buoyant spheres in turbulence, Zimmermann et al. (2011a) had originally shown the existence of an alignment between translation and rotation. Figure 8(e) and inset show PDFs of linear (a ) and angular accelerations ( ), respectively, of two buoyant N N spheres ( 0:89) that diﬀer solely in their rotational inertia ratios I = 1.0 and 0.6. For the lower I case, a strong coupling between translation and rotation ensues, which reﬂects strongly in the particle’s linear and angular acceleration PDFs (Mathai et al. 2018b). The same qualitative eﬀects were reproduced for buoyant cylindrical particles having diﬀerent I inertia (Mathai et al. 2017). New experiments are being extended to the realm of buoyant ellipsoidal particles (oblate to prolate) in turbulence (Will et al. 2019). Interesting new regimes are being revealed (see Figure 9) due to the coupling between particle buoyancy, particle shape, and turbulence. 3.4. Lift and shear-induced lateral migrations It is well known that bubbles or buoyant particles are inclined to a lift force when subjected to a mean shear in a ﬂow. For spherical particles rising in a simple shear ﬂow, the nature of these lift forces are by now well understood (Loth 2008a; Loth and Dorgan 2009). The case of bubbles, however, is quite diﬀerent due to additional complexities arising from deformability, internal circulations, and surface contamination (Legendre and Magnaudet 1997, 1998; Loth 2008b; Magnaudet and Eames 2000; Takemura et al. 2002; Van Nierop et al. 14 Mathai et al. Bubbles and Buoyant Particles in Turbulence Spherical bubbles Deformable bubbles (a) (b) Eo = 3 Eo = 0.9 (c) (d) Eo = 3 Eo = 3 -1 0 1 -1 0 1 Figure 10 Direct numerical simulations of rising bubbles in vertical channel ﬂows. (a) and (b) Distribution of spherical (Eo = 3) and deformable bubbles (Eo = 0.9), respectively, in an upward ﬂow turbulent channel at a friction Reynolds number Re = 127. The ratio of bubble diameter to channel halfwidth d =h = 0:3, and the mean bubble Reynolds number Re 136. Spherical p p bubbles migrate toward the walls, but deformable bubbles distribute themselves nearly uniformly in the bulk of the ﬂow. This is further elucidated by the tracks of the spherical (c) and deformable bubbles (d). Figures are adapted from Dabiri et al. (2013). 2007). The focus in most situations is to assess the lateral forces which induce bubble migration towards or away from the wall, and this in turn can be expressed as a function of bubble properties (deformability, size), and ﬂow properties (co-ﬂowing channel, counter-ﬂowing channel, or background turbulence level). Bubble laden wall-layers are a common observation, and have been extensively explored experimentally (Kitagawa et al. 2004; So et al. 2002; van Sint Annaland et al. 2006; Zhang et al. 2006). In the presence of turbulence, the two most common ﬂow conﬁgurations are upward (Nakoryakov et al. 1981) and downward (Kashinsky and Randin 1999) turbulent channel ﬂows. Although in single phase ﬂow the two are identical, they diﬀer greatly for two-phase situations, since the bubble buoyancy can now orient diﬀerently with respect to the mean shear near the channel walls. Drew (1993) developed the earliest model to unveil the mechanisms involved. Using an asymptotic analysis, the author was able to qualitatively reproduce the general trends for the velocity proﬁles and void fraction distributions. In comparison, DNS provide the ideal setting where the governing NS equations are solved numerically for both phases in such a way that all the length and time scales are fully resolved. Lu et al. (2006) used DNS with front tracking to examine bubbly ﬂows in a vertical channel. The results showed that for nearly spherical bubbles, the lift-induced lateral migration resulted in two regions: a nearly uniform velocity bulk ﬂow region where the weight of the liquid/bubble mixture balances the imposed pressure gradient, and a wall layer that is free of bubbles for downﬂow and wall-rich for upﬂow (Tryggvason and Lu 2015). The latter situation with bubble clustering near the walls is shown in Figure 10(a). While this result can be explained www.annualreviews.org 15 4 (a) (b) (d) 3.5 5 mm 2.5 1.5 0 1 2 3 4 5 (c) (e) b = 0.5 b = 1.0 b = 2.8 Experiment b = 4.8 -2 Almeras et al. (2017) -2 -4 -4 DNS -6 du Cluseau et al. (2019) -6 -8 -2 0 2 -1 0 1 2 10 10 10 Figure 11 Flow modiﬁcation in turbulent bubbly upﬂows. (a) Experimental images of 2-4 mm diameter deformable air bubbles in a vertical ﬂow at Re = 262. Snapshot adapted from the experiments reported in Alméras et al. (2017). (b) Volume-rendering of deformable bubbles in a turbulent ﬂow, obtained from DNS with front tracking method. Figure adapted from du Cluzeau et al. (2019). The ﬂow structures are colored by isovalues of the so-called criterion (Jeong and Hussain 1995). (c) Probability density functions (vertically shifted for better visibility) of the axial velocity in the bulk of the ﬂow from experimental (black) and DNS (blue) cases shown in (a) and (b), respectively. The PDF shapes are qualitatively similar in experiment and DNS. (d) Normalized vertical velocity ﬂuctuations of the liquid phase for diﬀerent values of the bubblance parameter b. (e) Normalized spectra of the bubble induced liquid ﬂuctuations. Distinct ranges of scalings can be observed, i.e. 5=3 for f=f < 1 and3 for f=f > 1. f is a characteristic cut-oﬀ frequency, which can be calculated using the bubble c c c and ﬂow properties. Figures (d) and (e) adapted from Alméras et al. (2017). in the same spirit as the shear-induced lift of rigid spherical particles, strikingly, when the bubble is deformable the eﬀect is reversed (see Figure 10(b)), and deformable bubbles distribute themselves in the bulk of the channel. The negligible lift for the deformable bubble can be attributed to the pliant nature of its interface, which prevents the buildup of a non-uniform surface pressure distribution. Thus, it is bubble deformability (Tomiyama et al. 2002), and not size, which causes the sign change of the lift force in turbulent upﬂow channels and pipes. The ﬂow modiﬁcations that the bubbles bring about in vertical channel ﬂows will be reserved for a later section. 4. TURBULENCE MODULATION BY BUBBLES In the preceding sections we reviewed a variety of dynamical regimes for bubbles and buoyant particles in turbulent ﬂows. However, it is the collective behavior of these particle which often contributes to sizable eﬀects in most engineering applications of particle-laden turbulence. Next, we will discuss how some of the afore-described mechanisms play a crucial role in triggering the diﬀerent kinds of ﬂow modulations occurring in low to moderate volume fraction ( < 5%) suspensions of bubble laden turbulence. 4.1. Bubbles rising within homogeneous turbulence 16 Mathai et al. Bubbles and Buoyant Particles in Turbulence PDF 4.1.1. Homogeneous isotropic turbulence. Homogeneous and isotropic turbulence (HIT) represents one of the most elementary forms of turbulence imaginable (Batchelor 1953). Notwithstanding how unrealistic this ﬂow might seem from a practical viewpoint, studies of HIT have led to great breakthroughs in our understanding of real-world turbulent ﬂows (Pope 2000). The assumption of statistical homogeneity and isotropy (among others) has been central to many successful theories of turbulent ﬂows. It may be noted that all single phase ﬂows, at high enough Reynolds numbers, will behave as HIT in the universal range and far from boundaries. For single phase turbulence, the 5=3 scaling of the energy spectrum (in the inertial range) is well known (Pope 2000), wherein the energy ﬂux ﬂowing down to smaller length scales is nearly constant up to the dissipative scales. In contrast, for a swarm of high Re bubbles rising within an otherwise quiescent liquid, bubble induced turbulence (BIT) leads to a3 scaling for the energy spectrum. This result, originally observed in the milestone work of Lance and Bataille (1991), has by now become well-established through detailed experiments, direct numerical simulations, and even large eddy simulations (Martínez-Mercado et al. 2007; Riboux et al. 2013; Roghair et al. 2011). A further simpliﬁed approach was undertaken by Mazzitelli and Lohse (2009), using the PP equation of motion with Du an imposed back-reaction on the ﬂow. This cumulative back reaction force f = g V (x x ), R p f p Dt acting at the point x in the ﬂow, did not generate the kind of liquid agitation commonly seen in BIT. Although the two-way coupling proved non-ideal for BIT, it demonstrated that the key ingredient crucial for the scaling (missing in the PP approach) was a model for the unsteady bubble wakes. For a detailed review of the liquid agitation induced by bubbles swarms (BIT), readers are referred to Risso (2017). Both single phase turbulence and bubbly swarms have been studied separately. But the situation where bubbles are injected into an already turbulent background ﬂow, despite its relevance in industrial applications, has only recently begun to be understood. What determines the nature of liquid ﬂuctuations and energy spectrum of such bubble laden turbulent ﬂows? To allow comparisons across diﬀerent levels of turbulence and bubble volume fractions, Lance and Bataille (1991) and building on that Rensen et al. (2005) introduced the so-called “bubblance” parameter b, which compares the intensity of BIT to the intensity of incident turbulence. This ratio of kinetic energies can be written as b = V =u , where V is the mean rise velocity of the bubbles2. Recently, r r Alméras et al. (2017) conducted extensive experiments in such turbulent bubbly ﬂows for a wide range of b. Figure 11(a) shows a snapshot of 2-4 mm diameter bubbles rising in the bulk region of upward channel ﬂow. The liquid velocity ﬂuctuations in bubble laden turbulence was measured using a phase-sensitive hot-ﬁlm anemometry (Alméras et al. 2017), and shown to be positively skewed (see Figure 11(c)), as is also the case for BIT. Similarly skewed velocity PDF have been seen in DNS of bubbles in turbulent upﬂows (du Cluzeau et al. 2019) using the front tracking method (see Figure 11(b) and (c)). In the experiments, turbulence attenuation was observed at low values of b (< 0:25), while the liquid velocity ﬂuctuations were augmented at larger b (see Figure 11(d)). It is important to stress that the liquid agitation produced by high Reynolds number bubbles is anisotropic (Risso 2017), the eﬀects of which are noticeable also in turbulent ﬂows. For details about the anisotropy of liquid velocity ﬂuctuations in bubble laden turbulence, the reader is referred to Alméras et al. (2019). With the magnitude and distributions of liquid agitation revealed, the natural next question to ask is: How is the energy spectrum modiﬁed by the presence of bubbles in an incident turbulent ﬂow? Figure 11(e) shows the spectrum of liquid velocity ﬂuctuations for b in the range [0, 1.3]. With the addition of bubbles, the higher frequencies of the inertial sub-range of single phase turbulence are substituted by a3 scaling of BIT. However, the -5/3 scaling appears to be preserved for the lower frequencies for all values of b tested. The characteristic cut-oﬀ frequency f imposed by the bubble swarm may then be calculated as V = , where = d =C where c r c c p d0 C is the drag coeﬃcient of an isolated bubble in still ﬂuid (Riboux et al. 2010). The above results only partially d0 resolve the complexity of the problem. One may note that the spectrum modiﬁcation by 2-4 mm diameter bubbles (reported above) is at variance with some of the observations in Risso (2017). The reason for this lies in the diﬀerences in the operating conditions, i.e. the size and Reynolds number of the bubbles. With a change in bubble diameter, the cut-oﬀ length scale is diﬀerent. This leaves room for much variability in the frequency range of BIT. These issues can be resolved through careful studies where bubble size and bubble Reynolds number are controlled independently. Such studies will however be challenging. 2Note that a prefactor of 1/2, based on the C of a spherical bubble (Van Wijngaarden 1998), was used in the deﬁnitions of b in Rensen et al. (2005) www.annualreviews.org 17 (a) (b) Figure 12 Flow modiﬁcation in turbulent bubbly upward channel ﬂows for 3% volume fraction of bubbles with a size ratio d =h = 0:3, where h is the channel half width. (a) Flow rate as a function of bubble deformablity, quantiﬁed here using Eötvös number Eo (or equivalently the Bond number Bo), for two diﬀerent values of friction Reynolds number Re . The friction Reynolds number eﬀect is negligible; the phenomenon is governed primarily by the bubble deformability. Bubble Reynolds number Re lies in the range [120, 200]. (b) Velocity proﬁle near the wall for the Eo = 3 (Re = 127) case in (a), as compared to the single phase case. Figures adapted from Dabiri et al. (2013). 4.1.2. Homogeneous shear turbulence. Homogeneous shear turbulence can be considered one of the simplest turbulent ﬂows where the ﬂow relaxes the condition of statistical isotropy, but maintains homogeneity over all spatial scales (Champagne et al. 1970; Pumir 1996). Explorations of bubble-laden homogeneous shear turbulence have been largely on the simulations side (Gualtieri et al. 2015). As a general rule, bubbles enhance the dissipation rate of turbulent kinetic energy (Rosti et al. 2019); however, its production rate can be either enhanced or diminished depending on the ﬂow parameters (Kawamura and Nakatani 2006). While the TKE production rate increases with the turbulent Reynolds number, it was found to decrease with the shear Reynolds number. In addition, the eﬀect of bubble deformability was to enhance the TKE production rates. Although simpliﬁed models have been proposed to explain these, the absence of direct experimental evidence leaves much debate on the validity of these mechanisms in causing turbulence modulation. 4.2. Bubble laden channel ﬂows 4.2.1. Horizontal channels and boundary layers. The injection of bubbles, small or large, in turbulent boundary layers and horizontal channel ﬂows has long been known to modulate the ﬂow, with several experimental explorations in the past few decades (Gutierrez-Torres et al. 2008; Madavan et al. 1984, 1985; Sanders et al. 2006). Recent numerical works include those due to Xu et al. (2002) and Ferrante and Elghobashi (2004) employing Eulerian-Lagrangian models, and those due to Lu et al. (2005) using DNS to reveal how bubbles can modify the near-wall vortical structures causing friction drag reduction even with relatively low volume fractions. Reviews on drag reduction by bubbles in turbulent boundary layers can be found in Murai (2014), Ceccio (2010), and more recently in Rawat et al. (2019). 4.2.2. Vertical channels. Turbulence modulation in bubbly vertical channel ﬂows has historically been the subject of numerous investigations, both experimental (Serizawa et al. 1975) and analytical (Antal et al. 1991). As earlier introduced in section 3.4, with high resolution DNS becoming increasingly feasible, simulations employing the front tracking method (Unverdi and Trygvasson 1992) have led to signiﬁcant breakthroughs in our understanding of these problems (Dabiri and Bhuvankar 2016; du Cluzeau et al. 2019; Lu and Tryggvason 2006, 2008; Tryggvason et al. 2011). Here we will highlight the main eﬀects caused by bubbles on the underlying ﬂow ﬁeld. For a bubble laden channel ﬂow, the relevant dimensionless parameters are the channel Reynolds number Re , the bubble Galileo number Ga, and the bubble Eötvös number Eo. The bubble Reynolds number, which is an output parameter Re , can be expressed as a function of Ga, Eo, and, to a weaker extent, the volume fraction . For upward channels ﬂows, the most prominent eﬀect is a reduction in the net ﬂow rate. The ﬂow rate reduction is 18 Mathai et al. Bubbles and Buoyant Particles in Turbulence (a) (b) (a) WSS Figure 13 Bubbly drag reduction (DR) in turbulent Taylor-Couette ﬂow. (a) Direct numerical simulations of Taylor-Couette (TC) turbulence at inner cylinder Reynolds number Re = 900 for single-phase ﬂow (top) and two-phase ﬂow laden with point-like bubbles (bottom). Wall shear stress (colorbar) is reduced in the two phase ﬂow with small bubbles, leading to about 20% DR. Figures adapted from Sugiyama et al. (2008). (b) DR% as function of time at a high Re = 2 10 . Nearly all DR is lost after injection of a surfactant (Triton X-100), which reduces the two phase ﬂow from a suspension of large deformable bubbles (upper left inset) to microbubbles (upper right inset). Data and ﬁgures from Verschoof et al. (2016). most extreme for spherical bubbles, since they occupy the near wall regions. Highly deformable bubbles, owing to their nearly zero (or slightly negative) lift force, remain in the bulk of the ﬂow, thus having negligible eﬀects on the volumetric ﬂow rate (see Figure 12(a)). The physical mechanism behind the ﬂow rate reduction is a sudden rise in the near-wall viscous dissipation when the spherical bubbles enter the viscous sublayer. As dictated by their relative distributions, the near-wall liquid velocity ﬂuctuations are enhanced for spherical bubbles, but the same occurs in the bulk for deformable bubbles. In addition the turbulent velocity ﬂuctuations in the bubble rich regions are enhanced, a result which directly follows from the knowledge of BIT (Risso 2017). The velocity proﬁle near the wall is retarded in the presence of bubbles (see Figure 12(b)), which is consistent with the observed ﬂow rate reductions. In the less commonly explored downﬂow conﬁguration with bubbles (Lu et al. 2006), a ﬂow rate reduction is accompanied by a suppression of turbulence in the near wall bubble-free layer, and turbulence augmentation in the bubble-rich bulk ﬂow (see also section 3.4). Thus, while in many engineering applications it might seem beneﬁcial to reduce the bubble size and maximize the interfacial area, such eﬀorts should not be made at the expense of a ﬂow rate reduction due to viscous losses of the bubbles accumulating near the walls. Lastly, we note that the above discussions pertain to clean bubbles. In reality, the void fraction proﬁles in turbulent bubbly upﬂow experiments (Serizawa et al. 1975) are possibly not as sharply peaked as seen in the simulations (DNS). When the gas-liquid interfaces are contaminated (with dirt or surfactants), a slight reduction of the lift force and lateral migration tendency can be expected (Clift et al. 1978; Lu et al. 2017; Takagi and Matsumoto 2011). Thus, even at low volume fractions of the dispersed phase, bubble dynamics can show fundamental diﬀerences to that of particles due to the presence of an internal circulation. On the other hand, dense bubbly ﬂows ( 5%) need to be treated as a markedly diﬀerent subject; for a review of experimental work on bubbles in vertical pipe/channel ﬂows, see Guet and Ooms (2006). Reviews discussing ﬂow regimes, operating parameters and design parameters of industrial bubble columns can be found in Besagni and Inzoli (2016). www.annualreviews.org 19 Two phase Single phase Deformable bubbles Microbubbles 4.3. Bubbly Taylor-Couette turbulence Taylor-Couette (TC) ﬂow, the ﬂow between two coaxial co- or counter-rotating cylinders, is one of the paradigmatic systems of ﬂuid physics (Grossmann et al. 2016). Unlike the turbulent ﬂows discussed above, the TC ﬂow is a closed system with exactly derivable balances between driving and dissipation. Owing to these particular beneﬁts, turbulent TC ﬂow experiments have been widely used by researchers as a model system to study bubble induced ﬂow modulations, bubble-vortex interactions and bubbly drag reduction (Chouippe et al. 2014; Fokoua et al. 2015; Murai et al. 2005, 2008; van den Berg et al. 2005; van Gils et al. 2013; Verschoof et al. 2016). A remarkable eﬀect of introducing bubbles in turbulent TC ﬂow is that with only a small percentage ( O(1%)) of bubbles, signiﬁcant turbulence modulation can be achieved. This usually manifests in major drag reduction, with various mechanisms contributing to it, see also section 9 of Lohse (2018). 4.3.1. Drag reduction in the buoyancy dominant regime. We ﬁrst discuss the studies on drag reduction by microbubbles in turbulent TC ﬂow. To this end, Sugiyama et al. (2008) and Spandan et al. (2016) conducted DNS employing an Euler-Lagrangian (PP) two-way coupled approach in a regime where the inner cylinder Reynolds number Re ! r (r r )= was in the range [600, 8000]. Here ! , r , and r are the inner cylinder ic ic ic oc ic ic ic oc angular velocity, inner cylinder radius, and outer cylinder radius, respectively. The drag reduction was deﬁned as hC i hC i f f s tp DR(%) = 100; 11. hC i 2 2 where C = ((1 r =r ) =)G=Re is the friction factor, G = is the dimensionless torque, with f ic oc 2 ic 2` the torque that is necessary to keep the inner cylinder of length ` rotating at constant angular velocity (Zhu et al. 2016). The subscripts ‘s’ and ‘tp’ denote single phase and two phase, respectively. The simulations of Sugiyama et al. (2008) and Spandan et al. (2016) showed that the buoyant motions of microbubbles can disrupt the coherent vortices (“Taylor rolls"), resulting in a reduction of drag (up to 20 %) on the inner cylinder surface (see Figure 13(a)). Spandan et al. (2016) varied the Froude number Fr = ! r =g, representing the ratio of ic ic ic centripetal force strength over buoyancy, in the range [0.16, 2.56]. Keeping the values of Re , and d ﬁxed, ic p drag reduction was found to be signiﬁcant at low Froude number (Fr . 1) and negligible at high Fr & 1. In this ic ic regime, the drag reduction becomes smaller with increasing Reynolds number. For a more generic quantiﬁcation of DR in this low Reynolds number regime, one might incorporate the bubble Stokes number and volume fraction as well, in which case the DR% can likely be expressed as a function of St/Fr and (assuming the bubbles are non-inertial). The hope is to obtain an overarching dimensionless parameter that can explain the degree of drag reduction in the so-called wavy vortex regime (Andereck et al. 1986; Fardin et al. 2014; Marcus 1984) of TC turbulence. 4.3.2. Drag reduction in highly turbulent regime. In the large Re regime, the stable coherent structure of the ic vortices in TC turbulence gets lost. Consequently, the eﬀect of the bubbles on the friction drag diminishes, as mentioned above. However, bubbles can still be used to reduce drag in the highly turbulent TC regime. In this high Reynolds number regime, however, the deformability of the bubble (i.e., large Weber number We > 1) is crucial for DR (van den Berg et al. 2005; van Gils et al. 2013). In this large Reynolds number regime, drag reduction increases with increasing Re . A direct experimental demonstration of the eﬀects of bubble deformability on ic turbulent TC drag reduction can be found in Verschoof et al. (2016). These authors dynamically changed the drag by adding a minute amount of surfactant (Triton X-100) to a highly turbulent TC ﬂow (Re up to 210 ) laden ic with deformable bubbles (see Figure 13(b)). In the original state with only a 4% volume fraction of deformable bubbles, the DR% was over 40% (left half of the ﬁgure). The addition of surfactant initiated a remarkable turn of events (breakup, coalescence prevention, etc) that caused the large deformable bubbles to be fully substituted by tiny microbbubles (right half of the ﬁgure), thereby reducing the DR to just 4%, which corresponds to the trivial eﬀect of the bubbles on the density and viscosity of the liquid. Recent work by Spandan et al. (2018) have used DNS to investigate the physical mechanisms of drag reduction in the turbulent regime (up to Re = 2 10 ). They connected the increase in drag reduction to a decrease in ic the dissipation in the wake of highly deformed bubbles near the inner cylinder. This touches the familiar territory of polymer drag reduction (Benzi and Ching 2018; Procaccia et al. 2008; White and Mungal 2008), and indicates 20 Mathai et al. Bubbles and Buoyant Particles in Turbulence interesting similarities in DR mechanisms where elastic properties of the dispersed phase are being exploited. Yet another important issue when studying bubbly drag reduction in TC ﬂow is the eﬀect of centripetal force on the bubble distribution in the ﬂow. If the TC ﬂow is in a fully laminar state (and Fr 1), all the bubbles ic should be pushed against the inner wall due to the radial pressure gradient induced by the centrifugally driven ﬂow. However, when the system is a highly turbulent state, the bubbles experience liquid velocity ﬂuctuations and pressure ﬂuctuations, which are enough to diﬀuse them towards the bulk region of the ﬂow (van Gils et al. 2013). The resulting bubble distribution in the gap between the cylinders will depend on the competing eﬀects of the turbulent pressure ﬂuctuations induced acceleration a and the centripetal acceleration a (r). This can be pf c deﬁned as a so-called3 centripetal Froude number a u =d pf p Fr (r) = = ; 12. a (r) U =r with U the mean azimuthal liquid velocity and r the radial position in the TC setup under consideration. van Gils 5 6 et al. (2013) estimated that Fr 1:6 at Re 5 10 and Fr 3:4 at Re = 1 10 . The lower Re implies c ic c ic ic a lower Fr and thus the eﬀective centripetal force on the bubbles is higher and, hence, the bubble accumulation is stronger near the inner cylinder wall at the expense of a lower concentration in the bulk. This reasoning is also consistent with their direct experimental observations. Thus to summarize, drag reduction in bubbly Taylor-Couette turbulence is a function of several parameters. While at moderate Re the buoyancy-induced drift of the microbubbles are suﬃcient, in the highly turbulent ic regime, buoyancy, deformability and centripetal eﬀects are all crucial ingredients to DR. In light of the close analogy that exists between TC ﬂows and pipe ﬂows (Eckhardt et al. 2007), the results obtained in turbulent TC ﬂows are of value to drag reduction research. Yet, whether and how the principles of turbulent two-phase TC ﬂow can be extended to pipelines and naval applications needs to be astutely examined. 4.4. Bubbles in turbulent convection In heat transfer systems, the motion of bubbles is known to be able to eﬃciently induce mixing of warm and cold parcels of liquid. For many industrial applications, injecting bubbles in the ﬂow can lead to a 100 times enhancement in the heat transfer coeﬃcient as compared to its single-phase counterpart (Deckwer 1980). Therefore, the eﬀect of bubbles and light particles on heat transfer has been subject of several experimental and numerical investigations. One approach for enhancing heat transport is to create vapor bubbles (Prosperetti 2017) or biphasic (Wang et al. 2019) and thermally expandable particles (Alards et al. 2019). Some of these have resulted in impressive heat transport enhancements as compared to the single phase case (Lakkaraju et al. 2013; Wang et al. 2019; Zhong et al. 2009). However, in the present review, we restrict our attention to gas phase bubbles. Early studies in forced convection systems with bubble injection (Sato et al. 1981a,b; Sekoguchi et al. 1980) showed that the bubbles modify the temperature proﬁle of the system, and the gas volume fraction close to the heated wall is important for heat transfer enhancement, i.e. higher void fractions close to the heated wall lead to an enhanced heat transfer. Using DNS with front tracking, Dabiri and Tryggvason (2015) recently studied the eﬀect of bubbles on the heat transfer rate in a ﬂow between two parallel walls under a constant heat ﬂux condition. They found that the bubbles stir up the viscous layer and reduce the size of the conduction region near the wall, resulting in an improved heat transfer eﬃciency, i.e. a 3% volume fraction of bubbles can increase the Nusselt number by 60%. Studies on the eﬀects of bubble injection on heat transfer in natural convection systems were mostly conducted with injection of micro-bubbles (Kitagawa et al. 2008, 2009) and sub-millimeter-bubbles (Kitagawa and Murai 2013) close to the heated wall. Kitagawa and Murai (2013) investigate the eﬀects of microbubble injection on natural convection heat transfer from a vertical heated plate in water. They found microbubble injection signiﬁcantly increases the heat transfer coeﬃcient in both the laminar and transition regimes. The enhancement 3Note that this Froude number should not be mistaken with the more widespread deﬁnition of Froude number where gravity force appears in the denominator. Here, the body force in the denominator is centrifugal. www.annualreviews.org 21 (a) (b) Thermistors Single phase 0.33 (c) b) ~ 2.5 mm bubbles Capillaries for air bubble injection Figure 14 (a) A sketch of a turbulent vertical convection setup (experiment) with bubble injection. Rectangular bubbly column heated from one sidewall and cooled from the other side, where the height H = 600 mm, and length between heated and cooled side walls L = 230 mm. Bubbles of diameter of about 2.5 mm were injected into the system through 180 capillaries (inner diameter 0.2 mm) placed at the bottom of the apparatus. (b) Nusselt number Nu in bubble laden vertical convection versus Rayleigh number Ra for various gas volume fractions . Blue circles shown for comparison correspond to the single phase turbulent vertical convection case ( = 0%). The size of the symbols corresponds to the typical error-bar in the data. (c) Heat transfer enhancement due to bubble injection Nu=Nu versus Ra . Here, Nu is the Nusselt number of the single-phase case. Figure 0 H 0 adapted from Gvozdić et al. (2018). ratio of the heat transfer coeﬃcient due to bubble injection is 1.6-2.0 in the laminar regime and 1.5-2.0 in the transition regime. The physical reason for heat transfer enhancement in the laminar regime is due to eﬀective mixing, whereas the physical reason for the heat transfer enhancement in the transition region is because microbubble injection accelerates the transition to turbulence (Kitagawa and Murai 2013). Deen and Kuipers (2013) studied wall-to-liquid heat transfer in dispersed gas liquid two-phase ﬂow using DNS, and found that a few high Reynolds number bubbles rising in quiescent liquid can considerably increase the local heat transfer between the liquid and a hot wall. Recently, Gvozdić et al. (2018) studied the eﬀect of deformable bubbles (with diameter of 2-3 mm) on heat transfer in a vertical natural convection setup, which was heated from one side and cooled from the other side (as shown in Figure 14(a)). The air bubbles were injected into the system using 180 capillaries (inner diameter 0.21 mm) uniformly distributed over the bottom of the nature convection system. The gas volume fraction varied 9 10 from 0.5% to 5%, and the Rayleigh number ranged from 4:0 10 to utmost 3:6 10 . Here, Rayleigh number g (T T )H Q =A c c is deﬁned as Ra = , and the Nusselt number as Nu = , where is the thermal c K (T T )=L h c expansion coeﬃcient, T and T are the mean temperatures of the hot and cold walls, respectively, L is the length h c of the setup, A is the surface of the sidewall, the thermal diﬀusivity, K the thermal conductivity of water, and Q is the measured power supplied to the heaters. For the entire range of and Ra , adding bubbles dramatically c H increased the heat transport eﬃciency, as the Nusselt number is about an order of magnitude higher as compared to single-phase ﬂow case (Figure 14(b)). In order to more clearly quantify the heat transport enhancement due to bubble injection, Figure 14(c) shows the ratio of the Nusselt number for two phase bubbly ﬂow, Nu, to that of the single-phase case, Nu , as a function of Ra at diﬀerent . It is shown that heat transfer was enhanced 0 H up to 20 times thanks to the bubble injection, and that the heat transfer enhancement increased with increasing the gas volume fraction and decreasing Rayleigh number Ra . Note that the decreasing trend of Nu=Nu with H 0 Ra occurred because the single-phase Nusselt number Nu increased with Ra whereas the two phase Nusselt H 0 H number did not change with Ra . It was found that the Nusselt number Nu was nearly independent of Ra and H H 0:45 depended solely on the gas volume fraction with the scaling of Nu / , which is suggestive of a diﬀusive transport mechanism as found in the case of the mixing of a passive tracer in a homogeneous bubbly ﬂow for a 22 Mathai et al. Bubbles and Buoyant Particles in Turbulence low gas volume fraction (Alméras et al. (2019); AlmÃľras et al. (2015)). Thus, bubble-induced mixing dominates the eﬃciency of the heat transfer in the moderate Ra bubbly nature convection systems. SUMMARY POINTS 1. The past decade has witnessed tremendous progress in our understanding of buoyant particle and bubble laden turbulent ﬂows. The addition of buoyant particles to turbulent ﬂows can modify key aspects of single phase turbulence, such as spectra or drag. This oﬀers opportunity to employ bubbles or light particles to tailor turbulent ﬂows to our beneﬁt. 2. As in many areas of modern ﬂuid dynamics, fully resolved direct numerical simulations oﬀer great potential, capable of explaining many intricate phenomena of two-phase turbulence. At the same time, a reduced treatment employing the Euler-Lagrangian approach prove suﬃcient in a remarkably large number of situations. The basic point particle formulation have been extended to include rotation through the Kelvin-Kirchhoﬀ equations, which have signiﬁcant predictive capabilities for buoyant particle and bubble dynamics in ﬂows. 3. Deformable bubble dynamics in turbulence is found to be governed by three fairly independent mech- anisms, which are as follows: The average bubble shape is mainly controlled by the relative motion between the bubble and the surrounding ﬂuid; the bubble velocity and orientation are a result of its own wake instability; and the eﬀect of turbulence reﬂects through random deformations of the bubble interface, which under extreme situations can cause bubble breakup. 4. Buoyancy brings about a multitude of modiﬁcations to particle dynamics in turbulence. For small bubbles and particles the crossing trajectories eﬀect leads to augmented particle accelerations, while in ﬁnite sized and ﬁnite particle Reynolds number cases, the wake induced accelerations add to the turbulent forcing. When the particle Reynolds number is increased further, the buoyant particle’s rotation further aids in the development of vigorous accelerated motions. 5. Adding bubbles to turbulent ﬂow is not synonymous with drag reduction. While bubbly drag reduction is possible in horizontal channel ﬂows, boundary layers and Taylor-Couette ﬂows, in vertical channel ﬂows (both upﬂows and downﬂows) the eﬀective drag is enhanced. The reason is that in vertical channels the bubbles increase the energy dissipation rate, while in Taylor-Couette turbulence and other ﬂows they result in suppression of the dissipation. 6. Air bubbles added to (open) turbulent convection systems dramatically enhance the heat transfer, thanks to their induced liquid agitation and mixing. FUTURE DIRECTIONS 1. Rigid buoyant anisotropic particles, including ellipsoidal, chiral and vaned particles, can add signiﬁcant amounts of energy to turbulent ﬂows. Their coupled translational-rotational dynamics is crucial to the liquid agitation. 2. Varying rotational inertia and/or center of mass location of buoyant particles presents exciting opportu- nities for turbulence modulation. 3. Whether the collective wake instabilities of rising bubbles and buoyant particles persist in intense turbulent environments or not is an open question. 4. The issue of energy spectra in bubble laden turbulence is only partially resolved. The eﬀects of the bubble size as compared to the Kolmogorov scale is not clear. Similarly, the spectrum modiﬁcation by the wakes of low to moderate Reynolds number bubbles remains to be elucidated. Furthermore, the behavior of bubble laden turbulent ﬂows in the limit of very large Reynolds number turbulence needs to be studied. 5. Tumbling buoyant particles can be engineered for turbulent downﬂow channels to aﬀect the near wall turbulence, with a potential for heat transfer enhancements. Bubbles cannot be used to achieve this. www.annualreviews.org 23 6. The response of bubbles to homogeneous shear turbulence remains to be experimentally explored. This requires the design of dedicated experimental setups. The relative alignment between buoyancy and shear, in combination with deformability is expected to induce symmetry breaking, and rich variability in bubble dynamics can be expected. DISCLOSURE STATEMENT The authors are not aware of any aﬃliations, memberships, funding, or ﬁnancial holdings that might be perceived as aﬀecting the objectivity of this review. ACKNOWLEDGMENTS We thank all coworkers for their contributions and for the many stimulating discussions over the years. The authors acknowledge support from the Natural Science Foundation of China under grant nos 11988102, 91852202, 11861131005 and 11672156, the Max Planck Center Twente for Complex Fluid Dynamics, NWO, and European Research Council (ERC) via an Advanced Grant for ﬁnancial support over the years. LITERATURE CITED Alards, K. M., Kunnen, R. P., Clercx, H. J., and Toschi, F. (2019). Statistical properties of thermally expandable particles in soft rayleigh-benard convection. arXiv preprint arXiv:1907.00049. Aliseda, A. and Lasheras, J. (2011). Preferential concentration and rise velocity reduction of bubbles immersed in a homogeneous and isotropic turbulent ﬂow. Phys. Fluids, 23(9):093301. Alméras, E., Mathai, V., Lohse, D., and Sun, C. (2017). Experimental investigation of the turbulence induced by a bubble swarm rising within incident turbulence. J. Fluid Mech., 825:1091–1112. Alméras, E., Mathai, V., Sun, C., and Lohse, D. (2019). Mixing induced by a bubble swarm rising through incident turbulence. Int. J. Multiph. Flow, 114:316–322. AlmÃľras, E., Risso, F., Roig, V., Cazin, S., Plais, C., and Augier, F. (2015). Mixing by bubble-induced turbulence. J. Fluid Mech., 776:458âĂŞ474. Andereck, C. D., Liu, S. S., and Swinney, H. L. (1986). Flow regimes in a circular Couette system with independently rotating cylinders. J. Fluid Mech., 164:155. Antal, S., Lahey Jr, R., and Flaherty, J. (1991). Analysis of phase distribution in fully developed laminar bubbly two-phase ﬂow. Int. J. Multiph. Flow, 17(5):635–652. Auguste, F. and Magnaudet, J. (2018). Path oscillations and enhanced drag of light rising spheres. J. Fluid Mech., 841:228–266. Auton, T., Hunt, J., and Prud’Homme, M. (1988). The force exerted on a body in inviscid unsteady non-uniform rotational ﬂow. J. Fluid Mech., 197:241–257. Balachandar, S. and Eaton, J. K. (2010). Turbulent dispersed multiphase ﬂow. Annu. Rev. Fluid Mech., 42:111–133. Batchelor, G. K. (1953). The theory of homogeneous turbulence. Cambridge university press. Bellani, G. and Variano, E. A. (2012). Slip velocity of large neutrally buoyant particles in turbulent ﬂows. New J. Phys., 14(12):125009. Benzi, R. and Ching, E. S. (2018). Polymers in ﬂuid ﬂows. Annu. Rev. Condens. Matter Phys., 9:163–181. Besagni, G. and Inzoli, F. (2016). Comprehensive experimental investigation of counter-current bubble column hydrodynamics: Holdup, ﬂow regime transition, bubble size distributions and local ﬂow properties. Chem. Eng. Sci., 146:259–290. Bodenschatz, E., Malinowski, S. P., Shaw, R. A., and Stratmann, F. (2010). Can we understand clouds without turbulence? Science, 327(5968):970–971. Brücker, C. (1999). Structure and dynamics of the wake of bubbles and its relevance for bubble interaction. Phys. Fluids, 11(7):1781–1796. C. R. Hunt, J., Perkins, R., and Fung, J. (1994). Review of the problem of modelling dispersed two-phase ﬂows. Multiphase Science and Technology, 8:595–643. Calzavarini, E., Cencini, M., Lohse, D., and Toschi, F. (2008a). Quantifying turbulence-induced segregation of inertial particles. Phys. Rev. Lett., 101(8):084504. 24 Mathai et al. Bubbles and Buoyant Particles in Turbulence Calzavarini, E., Huang, Y. X., Schmitt, F. G., and Wang, L. P. (2018). Propelled microprobes in turbulence. Phys. Rev. Fluids, 3(5):054604. Calzavarini, E., Kerscher, M., Lohse, D., and Toschi, F. (2008b). Dimensionality and morphology of particle and bubble clusters in turbulent ﬂow. J. Fluid Mech., 607:13–24. Calzavarini, E., Volk, R., Bourgoin, M., Leveque, E., Pinton, J. F., and Toschi, F. (2009). Acceleration statistics of ﬁnite-sized particles in turbulent ﬂow: the role of faxén forces. J. Fluid Mech., 630:179–189. Calzavarini, E., Volk, R., Leveque, E., Pinton, J. F., and Toschi, F. (2012). Impact of trailing wake drag on the statistical properties and dynamics of ﬁnite-sized particle in turbulence. Physica D: Nonlinear Phenomena, 241D(3):237–244. Cano-Lozano, J. C., Martinez-Bazan, C., Magnaudet, J., and Tchoufag, J. (2016). Paths and wakes of deformable nearly spheroidal rising bubbles close to the transition to path instability. Phys. Rev. Fluids, 1(5):053604. Ceccio, S. L. (2010). Friction drag reduction of external ﬂows with bubble and gas injection. Annu. Rev. Fluid Mech., 42:183–203. Champagne, F., Harris, V., and Corrsin, S. (1970). Experiments on nearly homogeneous turbulent shear ﬂow. J. Fluid Mech., 41(1):81–139. Chouippe, A., Climent, E., Legendre, D., and Gabillet, C. (2014). Numerical simulation of bubble dispersion in turbulent Taylor-Couette ﬂow. Phys. Fluids, 26(4):043304. Chouippe, A. and Uhlmann, M. (2015). Forcing homogeneous turbulence in direct numerical simulation of particulate ﬂow with interface resolution and gravity. Phys. Fluids, 27(12):123301. Clift, R., Grace, J., and Weber, M. (1978). Bubbles, Drops, and Particles. Academic Press. Crowe, C., Troutt, T., and Chung, J. (1996). Numerical models for two-phase turbulent ﬂows. Annu. Rev. Fluid Mech., 28(1):11–43. Csanady, G. (1963). Turbulent diﬀusion of heavy particles in the atmosphere. J. Atmos. Sci., 20(3):201–208. Dabiri, S. and Bhuvankar, P. (2016). Scaling law for bubbles rising near vertical walls. Phys. Fluids, 28(6):062101. Dabiri, S., Lu, J., and Tryggvason, G. (2013). Transition between regimes of a vertical channel bubbly upﬂow due to bubble deformability. Phys. Fluids, 25(10):102110. Dabiri, S. and Tryggvason, G. (2015). Heat transfer in turbulent bubbly ﬂow in vertical channels. Chem. Eng. Sci., 122:106 – Deane, G. B. and Stokes, M. D. (2002). Scale dependence of bubble creation mechanisms in breaking waves. Nature, 418(6900):839–844. Deckwer, W.-D. (1980). On the mechanism of heat transfer in bubble column reactors. Chem. Eng. Sci., 35(6):1341–1346. Deen, N. G. and Kuipers, J. (2013). Direct numerical simulation of wall-to liquid heat transfer in dispersed gas–liquid two-phase ﬂow using a volume of ﬂuid approach. Chem. Eng. Sci., 102:268–282. Douady, S., Couder, Y., and Brachet, M. (1991). Direct observation of the intermittency of intense vorticity ﬁlaments in turbulence. Phys. Rev. Lett., 67(8):983. Drew, D. A. (1993). Analytical modeling of multiphase ﬂow. Particulate two-phase ﬂow. du Cluzeau, A., Bois, G., and Toutant, A. (2019). Analysis and modelling of Reynolds stresses in turbulent bubbly up-ﬂows from direct numerical simulations. J. Fluid Mech., 866:132–168. Duineveld, P. (1995). The rise velocity and shape of bubbles in pure water at high reynolds number. J. Fluid Mech., 292:325–332. Eckhardt, B., Grossmann, S., and Lohse, D. (2007). Fluxes and energy dissipation in thermal convection and shear ﬂows. Europhys. Lett., 24001:78. Elghobashi, S. (1994). On predicting particle-laden turbulent ﬂows. Appl. Sci. Res., 52(4):309–329. Elghobashi, S. (2019). Direct numerical simulation of turbulent ﬂows laden with droplets or bubbles. Annu. Rev. Fluid Mech., 51:217–244. Ern, P., Risso, F., Fabre, D., and Magnaudet, J. (2012). Wake-induced oscillatory paths of bodies freely rising or falling in ﬂuids. Annu. Rev. Fluid Mech., 44:97–121. Fardin, M., Perge, C., and Taberlet, N. (2014). âĂĲThe hydrogen atom of ﬂuid dynamicsâĂİ–Introduction to the Taylor–Couette ﬂow for soft matter scientists. Soft Matter, 10(20):3523–3535. Faxén and Hilding (1922). Der widerstand gegen die bewegung einer starren kugel in einer zÃďhen ﬂÃĳssigkeit. Ann. Phys., 373(10):89. Ferrante, A. and Elghobashi, S. (2004). On the physical mechanisms of drag reduction in a spatially developing turbulent boundary layer laden with microbubbles. J. Fluid Mech., 503:345–355. Fiabane, L., Zimmermann, R., Volk, R., Pinton, J.-F., and Bourgoin, M. (2012). Clustering of ﬁnite-size particles in turbulence. Phys. Rev. E, 86(3):035301. Fokoua, G. N., Gabillet, C., Aubert, A., and Colin, C. (2015). Eﬀect of bubble arrangement on the viscous torque in bubbly Taylor-Couette ﬂow. Phys. Fluids, 27(3):034105. Ford, B. and Loth, E. (1998). Forces on ellipsoidal bubbles in a turbulent shear layer. Phys. Fluids, 10(1):178–188. www.annualreviews.org 25 Galper, A. and Miloh, T. (1995). Dynamic equations of motion for a rigid or deformable body in an arbitrary non-uniform potential ﬂow ﬁeld. J. Fluid Mech., 295:91–120. Govardhan, R. N. and Williamson, C. H. K. (2005). Vortex-induced vibrations of a sphere. J. Fluid Mech., 531:11–47. Grossmann, S., Lohse, D., and Sun, C. (2016). High Reynolds number Taylor-Couette turbulence. Annu. Rev. Fluid Mech., 48(1):53. Gualtieri, P., Battista, F., and Casciola, C. (2015). Transport of micro-bubbles in turbulent shear ﬂows. In J. Phys. Conf. Ser., volume 656, page 012018. IOP Publishing. Guet, S. and Ooms, G. (2006). Fluid mechanical aspects of the gas-lift technique. Annu. Rev. Fluid Mech., 38:225–249. Gutierrez-Torres, C., Hassan, Y., and Jimenez-Bernal, J. A. (2008). Turbulence structure modiﬁcation and drag reduction by microbubble injections in a boundary layer channel ﬂow. J. Fluids Eng., 130(11):111304. Gvozdić, B., Alméras, E., Mathai, V., Zhu, X., van Gils, D. P., Verzicco, R., Huisman, S. G., Sun, C., and Lohse, D. (2018). Experimental investigation of heat transport in homogeneous bubbly ﬂow. J. Fluid Mech., 845:226–244. Homann, H. and Bec, J. (2010). Finite-size eﬀects in the dynamics of neutrally buoyant particles in turbulent ﬂow. J. Fluid Mech., 651:81–91. Jeong, J. and Hussain, F. (1995). On the identiﬁcation of a vortex. J. Fluid Mech., 285:69–94. Jha, N. K. and Govardhan, R. (2015). Interaction of a vortex ring with a single bubble: bubble and vorticity dynamics. J. Fluid Mech., 773:460–497. Kashinsky, O. and Randin, V. (1999). Downward bubbly gas–liquid ﬂow in a vertical pipe. Int. J. Multiph. Flow, 25(1):109–138. Kawamura, T. and Nakatani, T. (2006). Direct numerical simulation of homogeneous turbulent shear ﬂow containing bubbles. In ASME 2006 2nd Joint US-European Fluids Engineering Summer Meeting Collocated With the 14th International Conference on Nuclear Engineering, pages 1571–1578. American Society of Mechanical Engineers. Kitagawa, A., Kosuge, K., Uchida, K., and Hagiwara, Y. (2008). Heat transfer enhancement for laminar natural convection along a vertical plate due to sub-millimeter-bubble injection. Exp. Fluids, 45(3):473–484. Kitagawa, A. and Murai, Y. (2013). Natural convection heat transfer from a vertical heated plate in water with microbubble injection. Chem. Eng. Sci., 99:215–224. Kitagawa, A., Sugiyama, K., and Murai, Y. (2004). Experimental detection of bubble–bubble interactions in a wall-sliding bubble swarm. Int. J. Multiph. Flow, 30(10):1213–1234. Kitagawa, A., Uchida, K., and Hagiwara, Y. (2009). Eﬀects of bubble size on heat transfer enhancement by sub-millimeter bubbles for laminar natural convection along a vertical plate. Int. J. Heat Fluid Flow, 30(4):778–788. La Porta, A., Voth, G. A., Crawford, A. M., Alexander, J., and Bodenschatz, E. (2001). Fluid particle accelerations in fully developed turbulence. Nature, 409(6823):1017–1019. Lakkaraju, R., Stevens, R. J., Oresta, P., Verzicco, R., Lohse, D., and Prosperetti, A. (2013). Heat transport in bubbling turbulent convection. Proc. Natl. Acad. Sci., 110(23):9237–9242. Lamb, H. (1993). Hydrodynamics. Cambridge university press. Lance, M. and Bataille, J. (1991). Turbulence in the liquid phase of a uniform bubbly air–water ﬂow. J. Fluid Mech., 222:95–118. Leal, L. (1980). Particle motions in a viscous ﬂuid. Annu. Rev. Fluid Mech., 12(1):435–476. Legendre, D. and Magnaudet, J. (1997). A note on the lift force on a spherical bubble or drop in a low-Reynolds-number shear ﬂow. Phys. Fluids, 9(11):3572–3574. Legendre, D. and Magnaudet, J. (1998). The lift force on a spherical bubble in a viscous linear shear ﬂow. J. Fluid Mech., 368:81?126. Lohse, D. (2018). Bubble puzzles: From fundamentals to applications. Phys. Rev. Fluids, 3(11):110504. Loisy, A. and Naso, A. (2017). Interaction between a large buoyant bubble and turbulence. Phys. Rev. Fluids, 2(1):014606. Loth, E. (2008a). Lift of a spherical particle subject to vorticity and/or spin. AIAA J., 46(4):801–809. Loth, E. (2008b). Quasi-steady shape and drag of deformable bubbles and drops. Int. J. Multiph. Flow, 34(6):523–546. Loth, E. and Dorgan, A. J. (2009). An equation of motion for particles of ﬁnite Reynolds number and size. Environ. Fluid Mech., 9(2):187–206. Lu, J., Biswas, S., and Tryggvason, G. (2006). A dns study of laminar bubbly ﬂows in a vertical channel. Int. J. Multiph. Flow, 32(6):643–660. Lu, J., Fernández, A., and Tryggvason, G. (2005). The eﬀect of bubbles on the wall drag in a turbulent channel ﬂow. Phys. Fluids, 17(9):095102. Lu, J., Muradoglu, M., and Tryggvason, G. (2017). Eﬀect of insoluble surfactant on turbulent bubbly ﬂows in vertical channels. Int. J. Multiph. Flow, 95:135–143. Lu, J. and Tryggvason, G. (2006). Numerical study of turbulent bubbly downﬂows in a vertical channel. Phys. Fluids, 18(10):103302. Lu, J. and Tryggvason, G. (2008). Eﬀect of bubble deformability in turbulent bubbly upﬂow in a vertical channel. Phys. Fluids, 20(4):040701. 26 Mathai et al. Bubbles and Buoyant Particles in Turbulence Madavan, N., Deutsch, S., and Merkle, C. (1984). Reduction of turbulent skin friction by microbubbles. Phys. Fluids, 27(2):356–363. Madavan, N., Deutsch, S., and Merkle, C. (1985). Measurements of local skin friction in a microbubble-modiﬁed turbulent boundary layer. J. Fluid Mech., 156:237–256. Magnaudet, J. and Eames, I. (2000). The motion of high-Reynolds-number bubbles in inhomogeneous ﬂows. Annu. Rev. Fluid Mech., 32(1):659–708. Marcus, P. S. (1984). Simulation of Taylor-Couette ﬂow. part 2. numerical results for wavy-vortex ﬂow with one travelling wave. J. Fluid Mech., 146:65–113. Martínez-Mercado, J., Palacios-Morales, C. A., and Zenit, R. (2007). Measurement of pseudoturbulence intensity in monodis- persed bubbly liquids for 10< re< 500. Phys. Fluids, 19(10):103302. Mathai, V., Calzavarini, E., Brons, J., Sun, C., and Lohse, D. (2016a). Microbubbles and microparticles are not faithful tracers of turbulent acceleration. Phys. Rev. Lett., 117(2):024501. Mathai, V., Huisman, S. G., Sun, C., Lohse, D., and Bourgoin, M. (2018a). Dispersion of air bubbles in isotropic turbulence. Phys. Rev. Lett., 121(5):054501. Mathai, V., Neut, M. W. M., van der Poel, E. P., and Sun, C. (2016b). Translational and rotational dynamics of a large buoyant sphere in turbulence. Exp. Fluids, 57(4):1–10. Mathai, V., Prakash, V. N., Brons, J., Sun, C., and Lohse, D. (2015). Wake-driven dynamics of ﬁnite-sized buoyant spheres in turbulence. Phys. Rev. Lett., 115(12):124501. Mathai, V., Zhu, X., Sun, C., and Lohse, D. (2017). Mass and moment of inertia govern the transition in the dynamics and wakes of freely rising and falling cylinders. Phys. Rev. Lett., 119(5):054501. Mathai, V., Zhu, X., Sun, C., and Lohse, D. (2018b). Flutter to tumble transition of buoyant spheres triggered by rotational inertia changes. Nat. Commun., 9(1):1792. Maxey, M. (1987). The gravitational settling of aerosol particles in homogeneous turbulence and random ﬂow ﬁelds. J. Fluid Mech., 174:441–465. Maxey, M. (2017). Simulation methods for particulate ﬂows and concentrated suspensions. Annu. Rev. Fluid Mech., 49:171–193. Maxey, M. R. and Riley, J. J. (1983). Equation of motion for a small rigid sphere in a nonuniform ﬂow. Phys. Fluids, 26(4):883–889. Mazzitelli, I. M. and Lohse, D. (2004). Lagrangian statistics for ﬂuid particles and bubbles in turbulence. New J. Phys., 6(1):203. Mazzitelli, I. M. and Lohse, D. (2009). Evolution of energy in ﬂow driven by rising bubbles. Phys. Rev. E, 79:066317. Mazzitelli, I. M., Lohse, D., and Toschi, F. (2003a). The eﬀect of microbubbles on developed turbulence. Phys. Fluids, 15(1):L5–L8. Mazzitelli, I. M., Lohse, D., and Toschi, F. (2003b). On the relevance of the lift force in bubbly turbulence. J. Fluid Mech., 488:283–313. Monchaux, R., Bourgoin, M., and Cartellier, A. (2010). Preferential concentration of heavy particles: a voronoï analysis. Phys. Fluids, 22(10):103304. Mougin, G. and Magnaudet, J. (2001). Path instability of a rising bubble. Phys. Rev. Lett., 88(1):014502. Mougin, G. and Magnaudet, J. (2002). The generalized kirchhoﬀ equations and their application to the interaction between a rigid body and an arbitrary time-dependent viscous ﬂow. Int. J. Multiph. Flow, 28(11):1837–1851. Mougin, G. and Magnaudet, J. (2006). Wake-induced forces and torques on a zigzagging/spiralling bubble. J. Fluid Mech., 567(00):185–194. Murai, Y. (2014). Frictional drag reduction by bubble injection. Exp. Fluids, 55(1773):1–10. Murai, Y., Oiwa, H., and Takeda, Y. (2005). Bubble behavior in a vertical Taylor-Couette ﬂow. In J. Phys. Conf. Ser., volume 14, page 143. IOP Publishing. Murai, Y., Oiwa, H., and Takeda, Y. (2008). Frictional drag reduction in bubbly Taylor-Couette ﬂow. Phys. Fluids, 20(3):034101. Nakoryakov, V., Kashinsky, O., Burdukov, A., and Odnoral, V. (1981). Local characteristics of upward gas-liquid ﬂows. Int. J. Multiph. Flow, 7(1):63–81. Naso, A. and Prosperetti, A. (2010). The interaction between a solid particle and a turbulent ﬂow. New J. Phys., 12(3):033040. Naso, A. and Pumir, A. (2005). Scale dependence of the coarse-grained velocity derivative tensor structure in turbulence. Phys. Rev. E, 72(5):056318. Obligado, M., Teitelbaum, T., Cartellier, A., Mininni, P., and Bourgoin, M. (2014). Preferential concentration of heavy particles in turbulence. J. Turb., 15(5):293–310. Parishani, H., Ayala, O., Rosa, B., Wang, L.-P., and Grabowski, W. (2015). Eﬀects of gravity on the acceleration and pair statistics of inertial particles in homogeneous isotropic turbulence. Phys. Fluids, 27(3):033304. Pope, S. B. (2000). Turbulent Flow. Cambridge University Press, Cambridge. Procaccia, I., L’vov, V. S., and Benzi, R. (2008). Colloquium: Theory of drag reduction by polymers in wall-bounded turbulence. www.annualreviews.org 27 Rev. Mod. Phys., 80:225–247. Prosperetti, A. (2015). Life and death by boundary conditions. J. Fluid Mech., 768:1–4. Prosperetti, A. (2017). Vapor bubbles. Annu. Rev. Fluid Mech., 49:221–248. Pumir, A. (1996). Turbulence in homogeneous shear ﬂows. Phys. Fluids, 8(11):3112–3127. Ravelet, F., Colin, C., and Risso, F. (2011). On the dynamics and breakup of a bubble rising in a turbulent ﬂow. Phys. Fluids, 23(10):103301. Rawat, S., Chouippe, A., Zamansky, R., Legendre, D., and Climent, E. (2019). Drag modulation in turbulent boundary layers subject to diﬀerent bubble injection strategies. Comput. Fluids, 178:73–87. Rensen, J., Bosman, D., Magnaudet, J., Ohl, C.-D., Prosperetti, A., Tögel, R., Versluis, M., and Lohse, D. (2001). Spiraling bubbles: How acoustic and hydrodynamic forces compete. Phys. Rev. Lett., 86(21):4819. Rensen, J., Luther, S., and Lohse, D. (2005). The eﬀect of bubbles on developed turbulence. J. Fluid Mech., 538:153–187. Riboux, G., Legendre, D., and Risso, F. (2013). A model of bubble-induced turbulence based on large-scale wake interactions. J. Fluid Mech., 719:362–387. Riboux, G., Risso, F., and Legendre, D. (2010). Experimental characterization of the agitation generated by bubbles rising at high Reynolds number. J. Fluid Mech., 643:509–539. Risso, F. (2017). Agitation, mixing, and transfers induced by bubbles. Annu. Rev. Fluid Mech., 50:25–48. Roghair, I., Mercado, J. M., Annaland, M. V. S., Kuipers, H., Sun, C., and Lohse, D. (2011). Energy spectra and bubble velocity distributions in pseudo-turbulence: Numerical simulations vs. experiments. Int. J. Multiph. Flow, 37(9):1093–1098. Rosti, M. E., Ge, Z., Suresh, S., Dodd, M., and Brandt, L. (2019). Droplets and bubbles in homogeneous shear turbulence. arXiv:1902.05259. Ryskin, G. and Leal, L. G. (1984). Numerical solution of free-boundary problems in ﬂuid mechanics. part 2. buoyancy-driven motion of a gas bubble through a quiescent liquid. J. Fluid Mech., 148:19–35. Sanders, W. C., Winkel, E. S., Dowling, D. R., Perlin, M., and Ceccio, S. L. (2006). Bubble friction drag reduction in a high-Reynolds-number ﬂat-plate turbulent boundary layer. J. Fluid Mech., 552:353–380. Sato, Y., Sadatomi, M., and Sekoguchi, K. (1981a). Momentum and heat transfer in two-phase bubble ﬂow i. theory. Int. J. Multiph. Flow, 7(2):167 – 177. Sato, Y., Sadatomi, M., and Sekoguchi, K. (1981b). Momentum and heat transfer in two-phase bubble ﬂow ii. a comparison between experimental data and theoretical calculations. Int. J. Multiph. Flow, 7(2):179 – 190. Sekoguchi, K., Nakazatomi, M., Sato, Y., and Tanaka, O. (1980). Forced convective heat transfer in vertical air-water bubble ﬂow. Bulletin of JSME, 23(184):1625–1631. Sengupta, A., Carrara, F., and Stocker, R. (2017). Phytoplankton can actively diversify their migration strategy in response to turbulent cues. Nature, 543(7646):555. Serizawa, A., Kataoka, I., and Michiyoshi, I. (1975). Turbulence structure of air-water bubbly ﬂowâĂŤii. local properties. Int. J. Multiph. Flow, 2(3):235–246. So, S., Morikita, H., Takagi, S., and Matsumoto, Y. (2002). Laser doppler velocimetry measurement of turbulent bubbly channel ﬂow. Exp. Fluids, 33(1):135–142. Spandan, V., Ostilla-Mónico, R., Verzicco, R., and Lohse, D. (2016). Drag reduction in numerical two-phase Taylor-Couette turbulence using an euler–lagrange approach. J. Fluid Mech., 798:411–435. Spandan, V., Verzicco, R., and Lohse, D. (2018). Physical mechanisms governing drag reduction in turbulent Taylor-Couette ﬂow with ﬁnite-size deformable bubbles. J. Fluid Mech., 849. Spelt, P. and Biesheuvel, A. (1997). On the motion of gas bubbles in homogeneous isotropic turbulence. J. Fluid Mech., 336:221–244. Sridhar, G. and Katz, J. (1995). Drag and lift forces on microscopic bubbles entrained by a vortex. Phys. Fluids, 7(2):389–399. Sugiyama, K., Calzavarini, E., and Lohse, D. (2008). Microbubbly drag reduction in Taylor-Couette ﬂow in the wavy vortex regime. J. Fluid Mech., 608:21–41. Tagawa, Y., Mercado, J. M., Prakash, V. N., Calzavarini, E., Sun, C., and Lohse, D. (2012). Three-dimensional lagrangian voronoï analysis for clustering of particles and bubbles in turbulence. J. Fluid Mech., 693:201–215. Takagi, S. and Matsumoto, Y. (2011). Surfactant eﬀects on bubble motion and bubbly ﬂows. Annu. Rev. Fluid Mech., 43:615–636. Takemura, F., Takagi, S., Magnaudet, J., and Matsumoto, Y. (2002). Drag and lift forces on a bubble rising near a vertical wall in a viscous liquid. J. Fluid Mech., 461:277?300. Thorpe, S. A. and Hall, A. J. (1987). Bubble clouds and temperature anomalies in the upper ocean. Nature, 328(6125):48–51. Tio, K.-K., Liñán, A., Lasheras, J. C., and Gañán-Calvo, A. M. (1993). On the dynamics of buoyant and heavy particles in a periodic Stuart vortex ﬂow. J. Fluid Mech., 254:671–699. Tomiyama, A., Tamai, H., Zun, I., and Hosokawa, S. (2002). Transverse migration of single bubbles in simple shear ﬂows. Chem. Eng. Sci., 57(11):1849–1858. 28 Mathai et al. Bubbles and Buoyant Particles in Turbulence Toschi, F. and Bodenschatz, E. (2009). Lagrangian properties of particles in turbulence. Annu. Rev. Fluid Mech., 41:375–404. Tryggvason, G. and Lu, J. (2015). Generalization of the kelvin–kirchhoﬀ equations for the motion of a body through a ﬂuid. Mech. Eng. Rev., 2(2):1–14. Tryggvason, G., Scardovelli, R., and Zaleski, S. (2011). Direct numerical simulations of gas–liquid multiphase ﬂows. Cambridge University Press. Unverdi, S. and Trygvasson, G. (1992). A front-tracking method for the computation of multiphase ﬂow. J. Comput. Phys., 100:25–37. van den Berg, T. H., Luther, S., Lathrop, D. P., and Lohse, D. (2005). Drag reduction in bubbly Taylor-Couette turbulence. Phys. Rev. Lett., 94:044501. van Gils, D. P., Narezo Guzman, D., Sun, C., and Lohse, D. (2013). The importance of bubble deformability for strong drag reduction in bubbly turbulent Taylor-Couette ﬂow. J. Fluid Mech., 722:317–347. Van Nierop, E. A., Luther, S., Bluemink, J. J., Magnaudet, J., Prosperetti, A., and Lohse, D. (2007). Drag and lift forces on bubbles in a rotating ﬂow. J. Fluid Mech., 571:439–454. van Sint Annaland, M., Dijkhuizen, W., Deen, N., and Kuipers, J. (2006). Numerical simulation of behavior of gas bubbles using a 3-d front-tracking method. AIChE J., 52(1):99–110. Van Wijngaarden, L. (1998). On pseudo turbulence. Theor. Comput. Fluid Dyna., 10(1-4):449–458. Veldhuis, C., Biesheuvel, A., van Wijngaarden, L., and Lohse, D. (2004). Motion and wake structure of spherical particles. Nonlinearity, 18(1):C1. Verschoof, R. A., Van Der Veen, R. C., Sun, C., and Lohse, D. (2016). Bubble drag reduction requires large bubbles. Phys. Rev. Lett., 117(10):104502. Volk, R., Calzavarini, E., Verhille, G., Lohse, D., Mordant, N., Pinton, J.-F., and Toschi, F. (2008a). Acceleration of heavy and light particles in turbulence: comparison between experiments and direct numerical simulations. Physica D: Nonlinear Phenomena, 237(14):2084–2089. Volk, R., Mordant, N., Verhille, G., and Pinton, J. F. (2008b). Laser doppler measurement of inertial particle and bubble accelerations in turbulence. Europhys. Lett., 81(3):34002. Voth, G. A., la Porta, A., Crawford, A. M., Alexander, J., and Bodenschatz, E. (2002). Measurement of particle accelerations in fully developed turbulence. J. Fluid Mech., 469:121–160. Voth, G. A. and Soldati, A. (2017). Anisotropic particles in turbulence. Annu. Rev. Fluid Mech., 49:249–276. Wang, L.-P. and Maxey, M. R. (1993). The motion of microbubbles in a forced isotropic and homogeneous turbulence. Appl. Sci. Res., 51(1-2):291–296. Wang, Z., Mathai, V., and Sun, C. (2019). Self-sustained biphasic catalytic turbulence. Nat. Commun., (in Press). White, C. M. and Mungal, M. G. (2008). Mechanics and prediction of turbulent drag reduction with polymer additives. Annu. Rev. Fluid Mech., 40(1):235–256. Will, J., Mathai, V., Krug, D., Huisman, S., Lohse, D., and Sun, C. (2019). Large ellipsoids rising in quiescent ﬂuids and turbulent ﬂows. unpublished. Wu, M. and Gharib, M. (2002). Experimental studies on the shape and path of small air bubbles rising in clean water. Phys. Fluids, 14(7):L49–L52. Xu, H., Pumir, A., Falkovich, G., Bodenschatz, E., Shats, M., Xia, H., Francois, N., and Boﬀetta, G. (2014). Flight–crash events in turbulence. Proc. Natl. Acad. Sci., 111(21):7558–7563. Xu, J., Maxey, M. R., and Karniadakis, G. E. (2002). Numerical simulation of turbulent drag reduction using micro-bubbles. J. Fluid Mech., 468:271–281. Zenit, R. and Magnaudet, J. (2008). Path instability of rising spheroidal air bubbles: a shape-controlled process. Phys. Fluids, 20(6):061702. Zhang, D., Deen, N., and Kuipers, J. (2006). Numerical simulation of the dynamic ﬂow behavior in a bubble column: A study of closures for turbulence and interface forces. Chem. Eng. Sci., 61(23):7593 – 7608. Zhang, W. and Stone, H. (1998). Oscillatory motions of circular disks and nearly spherical particles in viscous ﬂows. J. Fluid Mech., 367:329–358. Zhong, J.-Q., Funfschilling, D., and Ahlers, G. (2009). Enhanced heat transport by turbulent two-phase rayleigh-bénard convection. Phys. Rev. Lett., 102(12):124501. Zhu, X., Ostilla-Monico, R., Verzicco, R., and Lohse, D. (2016). Direct numerical simulation of Taylor-Couette ﬂow with grooved walls: torque scaling and ﬂow structure. J. Fluid Mech., 794:746–774. Zimmermann, R., Gasteuil, Y., Bourgoin, M., Volk, R., Pumir, A., and Pinton, J.-F. (2011a). Rotational intermittency and turbulence induced lift experienced by large particles in a turbulent ﬂow. Phys. Rev. Lett., 106(15):154501. Zimmermann, R., Gasteuil, Y., Bourgoin, M., Volk, R., Pumir, A., Pinton, J.-F., et al. (2011b). Tracking the dynamics of translation and absolute orientation of a sphere in a turbulent ﬂow. Rev. Sci. Inst., 82(3):033906. www.annualreviews.org 29

Physics – arXiv (Cornell University)

**Published: ** May 28, 2020

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.