▪ Abstract Predicting the motion of bubbles in dispersed flows is a key problem in fluid mechanics that has a bearing on a wide range of applications from oceanography to chemical engineering. In this review we synthesize the recent progress made in describing bubble motion in inhomogeneous flow. A trident approach consisting of experimental, analytical, and numerical work has given a clearer description of the hydrodynamic forces experienced by isolated bubbles moving either in inviscid flows or in slightly viscous laminar flows. A significant part of the paper is devoted to a discussion of drag, added-mass force, and shear-induced lift experienced by spheroidal bubbles moving in inertially dominated, time-dependent, rotational, nonuniform flows. The important influence of surfactants and shape distortion on bubble motion in a quiescent liquid is highlighted. Examples of bubble motion in inhomogeneous flows combining several of the effects mentioned above are discussed.

The Motion of High-Reynolds-Number Bubbles in Inhomogeneous Flows

We describe the results of 3D numerical simulations of oxygen shell burning and hydrogen core burning in a 23 M☉ stellar model. A detailed comparison is made to stellar mixing-length theory (MLT) for the shell-burning model. Simulations in 2D are significantly different from 3D, in terms of both flow morphology and velocity amplitude. Convective mixing regions are better predicted using a dynamic boundary condition based on the bulk Richardson number than by purely local, static criteria like Schwarzschild or Ledoux. MLT gives a good description of the velocity scale and temperature gradient for shell convection; however, there are other important effects that it does not capture, mostly related to the dynamical motion of the boundaries between convective and nonconvective regions. There is asymmetry between upflows and downflows, so the net kinetic energy flux is not zero. The motion of convective boundaries is a source of gravity waves; this is a necessary consequence of the deceleration of convective plumes. Convective overshooting is best described as an elastic response by the convective boundary, rather than ballistic penetration of the stable layers by turbulent eddies. The convective boundaries are rife with internal and interfacial wave motions, and a variety of instabilities arise that induce mixing through a process best described as turbulent entrainment. We find that the rate at which material entrainment proceeds at the boundaries is consistent with analogous laboratory experiments and simulation and observation of terrestrial atmospheric mixing. In particular, the normalized entrainment rate E = uE/σH is well described by a power-law dependence on the bulk Richardson number RiB = ΔbL/σ for the conditions studied, 20 RiB 420. We find E = ARi, with best-fit values log A = 0.027 ± 0.38 and n = 1.05 ± 0.21. We discuss the applicability of these results to stellar evolution calculations.

https://iopscience.iop.org/article/10.1086/520318/pdf

Turbulent convection in stellar interiors. I. Hydrodynamic simulation

The results of a laboratory experiment designed to study turbulent entrainment at sheared density interfaces are described. A stratified shear layer, across which a velocity difference ΔU and buoyancy difference Δb is imposed, separates a lighter upper turbulent layer of depth D from a quiescent, deep lower layer which is either homogeneous (two-layer case) or linearly stratified with a buoyancy frequency N (linearly stratified case). In the parameter ranges investigated the flow is mainly determined by two parameters: the bulk Richardson number RiB = ΔbD/ΔU2 and the frequency ratio fN = ND=ΔU. When RiB > 1.5, there is a growing significance of buoyancy effects upon the entrainment process; it is observed that interfacial instabilities locally mix heavy and light fluid layers, and thus facilitate the less energetic mixed-layer turbulent eddies in scouring the interface and lifting partially mixed fluid. The nature of the instability is dependent on RiB, or a related parameter, the local gradient Richardson number Rig = N2L/ (∂u/∂z)2, where NL is the local buoyancy frequency, u is the local streamwise velocity and z is the vertical coordinate. The transition from the Kelvin–Helmholtz (K-H) instability dominated regime to a second shear instability, namely growing Hölmböe waves, occurs through a transitional regime 3.2 < RiB < 5.8. The K-H activity completely subsided beyond RiB ∼ 5 or Rig ∼ 1. The transition period 3.2 < RiB < 5 was characterized by the presence of both K-H billows and wave-like features, interacting with each other while breaking and causing intense mixing. The flux Richardson number Rif or the mixing efficiency peaked during this transition period, with a maximum of Rif ∼ 0.4 at RiB ∼ 5 or Rig ∼ 1. The interface at 5 < RiB < 5.8 was dominated by ‘asymmetric’ interfacial waves, which gradually transitioned to (symmetric) Hölmböe waves at RiB > 5:8. Laser-induced fluorescence measurements of both the interfacial buoyancy flux and the entrainment rate showed a large disparity (as large as 50%) between the two-layer and the linearly stratified cases in the range 1.5 < RiB < 5. In particular, the buoyancy flux (and the entrainment rate) was higher when internal waves were not permitted to propagate into the deep layer, in which case more energy was available for interfacial mixing. When the lower layer was linearly stratified, the internal waves appeared to be excited by an ‘interfacial swelling’ phenomenon, characterized by the recurrence of groups or packets of K-H billows, their degeneration into turbulence and subsequent mixing, interfacial thickening and scouring of the thickened interface by turbulent eddies. Estimation of the turbulent kinetic energy (TKE) budget in the interfacial zone for the two-layer case based on the parameter α, where α = (−B + ε)/P, indicated an approximate balance (α ∼ 1) between the shear production P, buoyancy flux B and the dissipation rate ε, except in the range RiB < 5 where K-H driven mixing was active.

Entrainment and mixing in stratified shear flows

Recent developments in the physics and modeling of interfacial layers between regions with different turbulent intensities are reviewed. The flow dynamics across these layers governs exchanges of mass, momentum, energy, and scalars (e.g., temperature), which determine the growth, spreading, mixing, and reaction rates in many flows of engineering and natural interest. Results from several analytical and linearized models are reviewed. Particular attention is given to the case of turbulent/nonturbulent interfaces that exist at the edges of jets, wakes, mixing layers, and boundary layers. The geometry, dynamics, and scaling of these interfaces are reviewed, and future lines of research are suggested. The dynamics of passive and active scalars is also discussed, including the effects of stratification, turbulence level, and internal forcing. Finally, the modeling challenges for one-point closures and subgrid-scale models are briefly mentioned.

https://www.annualreviews.org/doi/pdf/10.1146/annurev-fluid-010313-141357

Interfacial Layers Between Regions of Different Turbulence Intensity

The results of a laboratory experiment designed to study turbulent entrainment at sheared density interfaces are described. A stratified shear layer, across which a velocity difference ΔU and buoyancy difference Δb is imposed, separates a lighter upper turbulent layer of depth D from a quiescent, deep lower layer which is either homogeneous (two-layer case) or linearly stratified with a buoyancy frequency N (linearly stratified case). In the parameter ranges investigated the flow is mainly determined by two parameters: the bulk Richardson number RiB = ΔbD/ΔU2 and the frequency ratio fN = ND=ΔU.When RiB > 1.5, there is a growing significance of buoyancy effects upon the entrainment process; it is observed that interfacial instabilities locally mix heavy and light fluid layers, and thus facilitate the less energetic mixed-layer turbulent eddies in scouring the interface and lifting partially mixed fluid. The nature of the instability is dependent on RiB, or a related parameter, the local gradient Richardson number Rig = N2L/
 (∂u/∂z)2, where NL is the local buoyancy frequency, u is the local streamwise velocity and z is the vertical coordinate. The transition from the Kelvin–Helmholtz (K-H) instability dominated regime to a second shear instability, namely growing Holmboe waves, occurs through a transitional regime 3.2 5:8.Laser-induced fluorescence measurements of both the interfacial buoyancy flux and the entrainment rate showed a large disparity (as large as 50%) between the two-layer and the linearly stratified cases in the range 1.5 < RiB < 5. In particular, the buoyancy flux (and the entrainment rate) was higher when internal waves were not permitted to propagate into the deep layer, in which case more energy was available for interfacial mixing. When the lower layer was linearly stratified, the internal waves appeared to be excited by an ‘interfacial swelling’ phenomenon, characterized by the recurrence of groups or packets of K-H billows, their degeneration into turbulence and subsequent mixing, interfacial thickening and scouring of the thickened interface by turbulent eddies.Estimation of the turbulent kinetic energy (TKE) budget in the interfacial zone for the two-layer case based on the parameter α, where α = (−B + e)/P, indicated an approximate balance (α ∼ 1) between the shear production P, buoyancy flux B and the dissipation rate e, except in the range RiB < 5 where K-H driven mixing was active.

A body moves at uniform speed in an unbounded inviscid fluid. Initially, the body is infinitely far upstream of an infinite plane of marked fluid; later, the body moves through and distorts the plane and, finally, the body is infinitely far downstream of the marked plane. Darwin (1953) suggested that the volume between the initial and final positions of the surface of marked fluid (the drift volume) is equal to the volume of fluid associated with the ‘added-mass’ of the body. We re-examine Darwin's (1953) concept of drift and, as an illustration, we study flow around a sphere. Two lengthscales are introduced: ρmax, the radius of a circular plane of marked particles; and x0, the initial separation of the sphere and plane. Numerical solutions and asymptotic expansions are derived for the horizontal Lagrangian displacement of fluid elements. These calculations show that depending on its initial position, the Lagrangian displacement of a fluid element can be either positive – a Lagrangian drift – or negative – a Lagrangian reflux. By contrast, previous investigators have found only a positive horizontal Lagrangian displacement, because they only considered the case of infinite x0. For finite x0, the volume between the initial and final positions of the plane of marked fluid is defined to be the ‘partial drift volume’, which is calculated using a combination of the numerical solutions and the asymptotic expansions. Our analysis shows that in the limit corresponding to Darwin's study, namely that both x0 and ρmax become infinite, the partial drift volume is not well-defined: the ordering of the limit processes is important. This explains the difficulties Darwin and others noted in trying to prove his proposition as a mathematical theorem and indicates practical, as well as theoretical, criteria that must be satisfied for Darwin's result to hold. We generalize our results for a sphere by re-considering the general expressions for Lagrangian displacement and partial drift volume. It is shown that there are two contributions to the partial drift volume. The first contribution arises from a reflux of fluid and is related to the momentum of the flow; this part is spread over a large area. It is well-known that evaluating the momentum of an unbounded fluid is problematic since the integrals do not converge; it is this first term which prevented Darwin from proving his proposition as a theorem. The second contribution to the partial drift volume is related to the kinetic energy of the flow caused by the body: this part is Darwin's concept of drift and is localized near the centreline. Expressions for partial drift volume are generalized for flow around arbitrary-shaped two- and three-dimensional bodies. The partial drift volume is shown to depend on the solid angles the body subtends with the initial and final positions of the plane of marked fluid. This result explains why the proof of Darwin's proposition depends on the ratio ρmax/x0. An example of drift due to a sphere travelling at the centre of a square channel is used to illustrate the differences between drift in bounded and unbounded flows.

/pdf/drift-partial-drift-and-darwin-s-proposition-4rd3qrzhvu.pdf

Drift, partial drift and Darwin's proposition

This paper presents a theoretical model of turbulence and mixing at a shear-free stable density interface. In one case (single-sided stirring) the interface separates a layer of fluid of density ρ in turbulent motion, with r.m.s. velocity uH and lengthscale LH, from a non-turbulent layer with density ρ+Δρ, while in the second case (double-sided stirring) the lower layer is also in turbulent motion. In both cases, the external Richardson number Ri=ΔbLH/
 u2H (where Δb is the buoyancy jump across the interface) is assumed to be large. Based on the hypotheses that the effect of the interface on the turbulence is as if it were suddenly imposed (which is equivalent to generating irrotational motions) and that linear waves are generated in the interface, the techniques of rapid distortion theory are used to analyse the linear aspects of the distortion of turbulence and of the interfacial motions. New physical concepts are introduced to account for the nonlinear aspects.To describe the spectra and variations of the r.m.s. fluctuations of velocity and displacements, a statistically steady linear model is used for frequencies above a critical frequency ωr/μc, where ωr(=Δb/2uH) is the maximum resonant frequency and μc<1. As in other nonlinear systems, observations below this critical frequency show the existence of long waves on the interface that can grow, break and cause mixing between the two fluid layers. A nonlinear model is constructed based on the fact that these breaking waves have steep slopes (which determines the form of the displacement spectrum) and on the physical argument that the energy of the vertical motions of these dissipative nonlinear waves should be comparable to that of the forced linear waves, which leads to an approximately constant value for the parameter μc. The model predictions of the vertical r.m.s. interfacial velocity, the interfacial wave amplitude and the velocity spectra agree closely with new and published experimental results.An exact unsteady inviscid linear analysis is used to derive the growth rate of the full spectrum, which asymptotically leads to the growth of resonant waves and to the energy transfer from the turbulent region to the wave motion of the stratified layer. Mean energy flux into the stratified layer, averaged over a typical wave cycle, is used to estimate the boundary entrainment velocity for the single-sided stirring case and the flux entrainment velocity for the double-sided stirring case, by making the assumption that the ratio of buoyancy flux to dissipation rate in forced stratified layers is constant with Ri and has the same value as in other stratified turbulent flows. The calculations are in good agreement with laboratory measurements conducted in mixing boxes and in wind tunnels. The contribution of Kelvin–Helmholtz instabilities induced by the velocity of turbulent eddies parallel to the interface is estimated to be insignificant compared to that of internal waves excited by turbulence.

Turbulence, waves and mixing at shear-free density interfaces. Part 1. A theoretical model

A laboratory experimental study was performed to investigate turbulence, waves and mixing at a sharp density interface (with a jump in buoyancy Δb), subjected to shear-free turbulence induced by oscillating grids with typical velocity and length scales of uH and LH, respectively. The cases where turbulence is present on one side (single-sided stirring) or on both sides (double-sided stirring) of the interface were considered. Extensive flow visualization studies and quantitative measurements were performed on the motion field and mixing characteristics at the interface. It was found that, rather than any one mechanism controlling the mixing process, different mechanisms (namely engulfment, generation of waves and their breaking, eddy impingement and Kelvin–Helmholtz billows) play dominant roles over different ranges of the bulk Richardson number Ri(Ri
 =ΔbLH/u2H). For the Ri range where wave generation is significant, certain hypotheses and predictions of the companion paper by Fernando & Hunt (1997) were tested in detail, by flow visualization studies of the qualitative properties of interfacial motions and quantitative measurements of the r.m.s. fluctuations of interfacial velocity and displacement, the local gradient Richardson number within the stratified layer, the frequency spectra and the related fractal properties of the interface. The results are consistent with the hypothesis that, at high values of Ri(>35), the density interface consists of linear internal waves driven by turbulence at high frequencies and breaking waves with sharp horizontal gradients of density at low frequencies.

Turbulence, waves and mixing at shear-free density interfaces. Part 2. Laboratory experiments

We present experimental and theoretical investigations into the problems of wind flow and scalar dispersion within and around a group of buildings.

Atmospheric Flow through Groups of Buildings and Dispersion from Localised Sources

We develop a theoretical analysis of the displacement of inviscid fluid particles and material surfaces caused by the unsteady flow around a solid body that is moving away from a wall. The body starts at position hs from the wall, and the material surface is initially parallel to the wall and at distance hL from it. A volume of fluid Df+ is displaced away from the wall and a volume Df- towards the wall. Df+ and Df- are found to be sensitive to the ratio hL/hs. The results of our specific calculations for a sphere can be extended in general to other shapes of bodies. When the sphere moves perpendicular to the wall the fluid displacement and drift volume Df+ are calculated numerically by computing the flow around the sphere. These numerical results are compared with analytical expressions calculated by approximating the flow around the sphere as a dipole moving away from the wall. The two methods agree well because displacement is an integrated effect of the fluid flow and the largest contribution to displacement is produced when the sphere is more than two radii away from the wall, i.e. when the dipole approximation adequately describes the flow. Analytic expressions for fluid displacement are used to calculate Df+ when the sphere moves at an acute angle α away from the wall. In general the presence of the wall reduces the volume displaced forward and this effect is still significant when the sphere starts 100 radii from the wall. A sphere travelling perpendicular to the wall, α = 0, displaces forward a volume Df+(0) = 4πa3hL/33/2hS when the marked surface starts downstream, or behind the sphere, and displaces a volume Df+(0) [similar] 2πa3/3 forward when it is marked upstream or in front of the body. A sphere travelling at an acute angle away from the wall displaces a volume Df+(α) [similar] Df+(0) cos α forward when the surface starts downstream of the sphere. When the marked surface is initially upstream of the sphere, there are two separate regions displaced forward and a simple cosine dependence on α is not found. These results can all be generalized to calculate material surfaces when the sphere moves at variable speed, displacements no longer being expressed in terms of time, but in relation to the distance travelled by the sphere.

/pdf/displacement-of-inviscid-fluid-by-a-sphere-moving-away-from-29grxkqzwb.pdf

J. C. R. Hunt

Papers

Drift, partial drift and Darwin's proposition

Turbulence, waves and mixing at shear-free density interfaces. Part 1. A theoretical model

Turbulence, waves and mixing at shear-free density interfaces. Part 2. Laboratory experiments

Atmospheric Flow through Groups of Buildings and Dispersion from Localised Sources

Displacement of inviscid fluid by a sphere moving away from a wall