banner image
No data available.
Please log in to see this content.
You have no subscription access to this content.
No metrics data to plot.
The attempt to load metrics for this article has failed.
The attempt to plot a graph for these metrics has failed.
Axisymmetric bubble collapse in a quiescent liquid pool. I. Theory and numerical simulations
Rent this article for


Image of FIG. 1.
FIG. 1.

(a) Numerical shapes of an axisymmetric bubble formed from a vertical underwater nozzle for . Note that, at a length scale of the order of , the shape of the bubble apparently forms two cones with a semiangle of close to the pinch-off point. (b) Detail of the bubble profiles—displaced in so that the minimum is located at —close to the pinch-off region. Note that the profiles become more symmetric with respect to the plane as pinch-off is approached. The values of in curves 1, 2, and 3 are, respectively, , , and . (c) Numerical bubble shapes represented in (b) as a function of the dimensionless stretched coordinates and . Note that the bubble shapes become more and more locally slender as pinch-off is approached. (d) In solid lines, full potential flow numerical simulations take into account the nozzle-bubble interaction and were previously represented in Figs. 1(b) and 1(c). In dashed lines, simplified potential flow numerical simulations in which an isolated bubble breaks symmetrically around the plane for a value of the Bond number, . Observe that there are no appreciable differences between the two types of simulations in the region near the minimum radius.

Image of FIG. 2.
FIG. 2.

Velocity at the minimum radius for a bubble in water, and two different gases (air, and , ). The details of the inviscid potential flow numerical simulations, represented in solid lines, will be provided in Sec. III. Note that, in agreement with Ref. 17 the radial velocity for is . Equation (6) is represented for different initial values of with dashed lines and differ from the numerics no matter how small the initial value of is.

Image of FIG. 3.
FIG. 3.

Values of the local Weber number and of the local gas and liquid Reynolds numbers ( and , respectively) for the inviscid potential flow numerical simulations depicted in Fig. 2. Note that the evolution of the local Weber number is very similar to the experimental one reported in Fig. 46 of Ref. 17.

Image of FIG. 4.
FIG. 4.

(a) Dimensionless axial gas velocity profiles for different values of the local gas Reynolds number, ; (b) continuous line: dimensionless gas pressure gradient defined in Eq. (18) and given by Eq. (22); dotted line: , which proves to be a good approximation to the real solution.

Image of FIG. 5.
FIG. 5.

Sketch of the geometry used for the symmetric type of simulations.

Image of FIG. 6.
FIG. 6.

In solid lines, the velocity at the minimum radius in the collapse of bubbles of two different gases in water [(a) air, and (b) , ]. The numerical computations have been performed using the inviscid symmetric code and the Bernoulli equation (35). Dashed lines: the theoretical result obtained integrating the two-dimensional Rayleigh equations (25) and (26) in the limits Re, with initial conditions for , , and the values of the numerical simulations when (a) , (b) . In (a) and (b) . The results of integrating Eqs. (25) and (26) in the limits Re, , and have also been included in each figure. While no appreciable differences between and are observed in the case of air, gas density plays a key role in the good agreement between theory and numerics in the case of .

Image of FIG. 7.
FIG. 7.

Velocity at the minimum radius of an air bubble that collapses within liquids of different viscosities. The computations are performed using the symmetric code and the Bernoulli equation (40).

Image of FIG. 8.
FIG. 8.

The experimental time evolution of the bubble minimum radius is represented using dashed lines for two different values of the liquid viscosity (adapted from the experiments in Ref. 17, where it is reported that the exponents of the power law for the cases of the liquids of 4.2 and 21 cp are, respectively, 0.6 and 0.67). Solid lines: the resulting time evolution calculated through the integration of the two-dimensional Rayleigh-like equations (25) and (26).

Image of FIG. 9.
FIG. 9.

(a) Comparison between the numerical computations depicted in Figs. 6 and 7 (solid lines) and the theoretical results obtained integrating the two-dimensional Rayleigh equations (25) and (26) using as initial conditions for , , and the values of the numerical simulations when oscillates between 120 and . Similarly, the comparison between numerics and theory in the case of is given in (b). In this case, the initial conditions for , , and are the values of the numerical simulation when . The instant at which , which is the condition for satellite formation, is indicated for each viscosity using a vertical arrow: (a) , (b) .

Image of FIG. 10.
FIG. 10.

(a) Satellite formation process in the case of air bubbles collapsing in a liquid with a viscosity of . In (b), the gas is and the material properties of the liquid are those of water.

Image of FIG. 11.
FIG. 11.

(a) Sketch of the gas flow during the initial stages of bubble collapse, where it is expected that flow separation does not occur. (b) During the latest instants of bubble pinch-off, gas flow may separate and the stagnant gas pressure is not recovered from the exit of the tube to the main bubble.

Image of FIG. 12.
FIG. 12.

Comparison between the results in Fig. 6 (, solid line) and those obtained integrating systems (B3) and (B4) for the same values of the initial conditions (dashed lines) and .


Generic image for table
Table I.

Physical properties of the different gases and liquid (water) considered.


Article metrics loading...


Full text loading...

This is a required field
Please enter a valid email address
752b84549af89a08dbdd7fdb8b9568b5 journal.articlezxybnytfddd
Scitation: Axisymmetric bubble collapse in a quiescent liquid pool. I. Theory and numerical simulations