▪ 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 three-dimensional flow around a spherical bubble moving steadily in a viscous linear shear flow is studied numerically by solving the full Navier-Stokes equations. The bubble surface is assumed to be clean so that the outer flow obeys a zero-shearstress condition and does not induce any rotation of the bubble. The main goal of the present study is to provide a complete description of the lift force experienced by the bubble and of the mechanisms responsible for this force over a wide range of Reynolds number (0.1 Re 500, Re being based on the bubble diameter) and shear rate (0 Sr 1, Sr being the ratio between the velocity difference across the bubble and the relative velocity). For that purpose the structure of the flow field, the influence of the Reynolds number on the streamwise vorticity field and the distribution of the tangential velocities at the surface of the bubble are first studied in detail. It is shown that the latter distribution which plays a central role in the production of the lift force is dramatically dependent on viscous effects. The numerical results concerning the lift coefficient reveal very different behaviours at low and high Reynolds numbers. These two asymptotic regimes shed light on the respective roles played by the vorticity produced at the bubble surface and by that contained in the undisturbed flow. At low Reynolds number it is found that the lift coefficient depends strongly on both the Reynolds number and the shear rate. In contrast, for moderate to high Reynolds numbers these dependences are found to be very weak. The numerical values obtained for the lift coefficient agree very well with available asymptotic results in the low-and high-Reynolds-number limits. The range of validity of these asymptotic solutions is specified by varying the characteristic parameters of the problem and examining the corresponding evolution of the lift coefficient. The numerical results are also used for obtaining empirical correlations useful for practical calculations at finite Reynolds number. The transient behaviour of the lift force is then examined. It is found that, starting from the undisturbed flow, the value of the lift force at short time differs from its steady value, even when the Reynolds number is high, because the vorticity field needs a finite time to reach its steady distribution. This finding is confirmed by an analytical derivation of the initial value of the lift coefficient in an inviscid shear flow. Finally, a specific investigation of the evolution of the lift and drag coefficients with the shear rate at high Reynolds number is carried out. It is found that when the shear rate becomes large, i.e. Sr l O(1), a small but consistent decrease of the lift coefficient occurs while a very significant increase of the drag coefficient, essentially produced by the modifications of the pressure distribution, is observed. Some of the foregoing results are used to show that the well-known equality between the added mass coefficient and the lift coefficient holds only in the limit of w...

This work reports the first part of a series of numerical simulations carried out in order to improve knowledge of the forces acting on a sphere embedded in accelerated flows at finite Reynolds number, Re. Among these forces added mass and history effects are particularly important in order to determine accurately particle and bubble trajectories in real flows. To compute these hydrodynamic forces and more generally to study spatially or temporally accelerated flows around a sphere, the full Navier–Stokes equations expressed in velocity–pressure variables are solved by using a finite-volume approach. Computations are carried out over the range 0.1 ≤ Re ≤ 300 for flows around both a rigid sphere and an inviscid spherical bubble, and a systematic comparison of the flows around these two kinds of bodies is presented.Steady uniform flow is first considered in order to test the accuracy of the simulations and to serve as a reference case for comparing with accelerated situations. Axisymmetric straining flow which constitutes the simplest spatially accelerated flow in which a sphere can be embedded is then studied. It is shown that owing to the viscous boundary condition on the body as well as to vorticity transport properties, the presence of the strain modifies deeply the distribution of vorticity around the sphere. This modification has spectacular consequences in the case of a rigid sphere because it influences strongly the conditions under which separation occurs as well as the characteristics of the separated region. Another very original feature of the axisymmetric straining flow lies in the vortex-stretching mechanism existing in this situation. In a converging flow this mechanism acts to reduce vorticity in the wake of the sphere. In contrast when the flow is divergent, vorticity produced at the surface of the sphere tends to grow indefinitely as it is transported downstream. It is shown that in the case where such a diverging flow extends to infinity a Kelvin–Helmholtz instability may occur in the wake.Computations of the hydrodynamic force show that the effects of the strain increase rapidly with the Reynolds number. At high Reynolds numbers the total drag is dramatically modified and the evaluation of the pressure contribution shows that the sphere undergoes an added mass force whose coefficient remains the same as in inviscid flow or in creeping flow, i.e. CM = ½, whatever the Reynolds number. Changes found in vorticity distribution around the rigid sphere also affect the viscous drag, which is markedly increased (resp. decreased) in converging (resp. diverging) flows at high Reynolds numbers.

We model the problem of path instability of a rising bubble by considering the bubble as a spheroidal body of fixed shape, and we solve numerically the coupled fluid-body problem. Numerical results show that this model exhibits path instability for large enough values of the control parameters. The corresponding characteristics of the zigzag and spiral paths are in good agreement with experimental observations. Analysis of the vorticity field behind the bubble reveals that a wake instability leading to a double threaded wake is the primary cause of the path instability.

Leaves falling in air and bubbles rising in water provide daily examples of nonstraight paths associated with the buoyancy-driven motion of a body in a fluid. Such paths are relevant to a large variety of applicative fields such as mechanical engineering, aerodynamics, meteorology, and the biomechanics of plants and insect flight. Although the problem has attracted attention for ages, it is only recently that the tremendous progress in the development of experimental and computational techniques and the emergence of new theoretical concepts have led to a better understanding of the underlying physical mechanisms. This review attempts to bring together the main recent experimental, computational, and theoretical advances obtained on this fascinating subject. To this end it describes the first steps of the transition in the wake of a fixed body and its connection with the onset and development of the path instability of moving bodies. Then it analyzes the kinematics and dynamics of various types of bodies along typical nonstraight paths and how the corresponding information can be used to build low-dimensional predictive models.

International audienceWe consider the generic problem of wake instabilities past fixed axisymmetric bodies, and focus on the extreme cases of a sphere and a flat disk. Numerical results reveal that the wakes of these two bodies evolve differently as the Reynolds number is increased. Especially, two new vortex shedding modes are identified behind a disk. To interpret these results, we introduce a model based on the theory of mode interactions in presence of O(2) symmetry. This model, which was initially developed for the Taylor-Couette system, allows us to explain the structural differences observed in the evolution of the two types of wakes and to accurately predict the evolution of the lift force

This paper reports on an experimental study of the motion of freely rising axisym- metric rigid bodies in a low-viscosity fluid. We consider flat cylinders with height h smaller than the diameter d and density ρb close to the density ρf of the fluid. We have investigated the role of the Reynolds number based on the mean rise velocity um in the range 80 ≤ Re = umd/ν ≤ 330 and that of the aspect ratio in the range 1.5 ≤ χ = d/h ≤ 20. Beyond a critical Reynolds number, Rec, which depends on the aspect ratio, both the body velocity and the orientation start to oscillate periodically. The body motion is observed to be essentially two-dimensional. Its description is particularly simple in the coordinate system rotating with the body and having its origin fixed in the laboratory; the axial velocity is then found to be constant whereas the rotation and the lateral velocity are described well by two harmonic functions of time having the same angular frequency, ω. In parallel, direct numerical simulations of the flow around fixed bodies were carried out. They allowed us to determine (i) the threshold, Recf1(χ), of the primary regular bifurcation that causes the breaking of the axial symmetry of the wake as well as (ii) the threshold, Recf2(χ), and frequency, ωf, of the secondary Hopf bifurcation leading to wake oscillations. As χ increases, i.e. the body becomes thinner, the critical Reynolds numbers, Recf1 and Recf2, decrease. Introducing a Reynolds number Re* based on the velocity in the recirculating wake makes it possible to obtain thresholds $\hbox{\it Re}^*_{cf1}$ and $\hbox{\it Re}^*_{cf2}$ that are independent of χ. Comparison with fixed bodies allowed us to clarify the role of the body shape. The oscillations of thick moving bodies (χ < 6) are essentially triggered by the wake instability observed for a fixed body: Rec(χ) is equal to Recf1(χ) and ω is close to ωf. However, in the range 6 ≤ χ ≤ 10 the flow corrections induced by the translation and rotation of freely moving bodies are found to be able to delay the onset of wake oscillations, causing Rec to increase strongly with χ. An analysis of the evolution of the parameters characterizing the motion in the rotating frame reveals that the constant axial velocity scales with the gravitational velocity based on the body thickness, $\sqrt{((\rho_f-\rho_b)/\rho_f)\,gh}$, while the relevant length and velocity scales for the oscillations are the body diameter d and the gravitational velocity based on d, $\sqrt{((\rho_f-\rho_b)/\rho_f)\,gd}$, respectively. Using this scaling, the dimensionless amplitudes and frequency of the body's oscillations are found to depend only on the modified Reynolds number, Re*; they no longer depend on the body shape.

This paper reports the results of a numerical investigation of the transient evolution of the flow around a spherical bubble rising in a liquid contaminated by a weakly soluble surfactant. For that purpose the full Navier-Stokes equations are solved together with the bulk and interfacial surfactant concentration equations, using values of the physical-chemical constants of a typical surfactant characterized by a simple surface kinetics. The whole system is strongly coupled by nonlinear boundary conditions linking the diffusion flux and the interfacial shear stress to the interfacial surfactant concentration and its gradient. The influence of surfactant characteristics is studied by varying arbitrarily some physical-chemical parameters. In all cases, starting from the flow around a clean bubble, the results describe the temporal evolution of the relevant scalar and dynamic interfacial quantities as well as the changes in the flow structure and the increase of the drag coefficient. Since surface diffusion is extremely weak compared to advection, part of the bubble (and in certain cases all the interface) tends to become stagnant. This results in a dramatic increase of the drag which in several cases reaches the value corresponding to a rigid sphere. The present results confirm the validity of the well-known stagnant-cap model for describing the flow around a bubble contaminated by slightly soluble surfactants. They also show that a simple relation between the cap angle and the bulk concentration cannot generally be obtained because diffusion from the bulk plays a significant role.