Research Papers: Multiphase Flows

Euler–Lagrange Simulations of Bubble Cloud Dynamics Near a Wall

[+] Author and Article Information
Jingsen Ma

Dynaflow, Inc.,
10621-J Iron Bridge Road,
Jessup, MD 20794
e-mail: jingsen@dynaflow-inc.com

Chao-Tsung Hsiao, Georges L. Chahine

Dynaflow, Inc.,
10621-J Iron Bridge Road,
Jessup, MD 20794

1Corresponding author.

Contributed by the Fluids Engineering Division of ASME for publication in the JOURNAL OF FLUIDS ENGINEERING. Manuscript received March 12, 2014; final manuscript received October 13, 2014; published online December 3, 2014. Assoc. Editor: Mark R. Duignan.

J. Fluids Eng 137(4), 041301 (Apr 01, 2015) (10 pages) Paper No: FE-14-1128; doi: 10.1115/1.4028853 History: Received March 12, 2014; Revised October 13, 2014; Online December 03, 2014

We present in this paper a two-way coupled Eulerian–Lagrangian model to study the dynamics of clouds of microbubbles subjected to pressure variations and the resulting pressures on a nearby rigid wall. The model simulates the two-phase medium as a continuum and solves the Navier–Stokes equations using Eulerian grids with a time and space varying density. The microbubbles are modeled as interacting singularities representing moving and oscillating spherical bubbles, following a modified Rayleigh–Plesset–Keller–Herring equation and are tracked in a Lagrangian fashion. A two-way coupling between the Euler and Lagrange components is realized through the local mixture density determined by the bubbles' volume change and motion. Using this numerical framework, simulations involving a large number of bubbles were conducted under driving pressures at different frequencies. The results show that the frequency of the driving pressure is critical in determining the overall dynamics: either a collective strongly coupled cluster behavior or nonsynchronized weaker multiple bubble oscillations. The former creates extremely high pressures with peak values orders of magnitudes higher than that of the excitation pressure. This occurs when the driving frequency matches the natural frequency of the bubble cloud. The initial distance between the bubble cloud and the wall also affects significantly the resulting pressures. A bubble cloud collapsing very close to the wall exhibits a cascading collapse, with the bubbles farthest from the wall collapsing first and the nearest ones collapsing last, thus the energy accumulates and this results in very high pressure peaks at the wall. At farther cloud distances from the wall, the bubble cloud collapses quasi-spherically with the cloud center collapsing last.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Brujan, E., Ikeda, T., Yoshinaka, K., and Matsumoto, Y., 2011, “The Final Stage of the Collapse of a Cloud of Bubbles Close to a Rigid Boundary,” Ultrason. Sonochem., 18(1), pp. 59–64. [CrossRef] [PubMed]
Brujan, E., Ikeda, T., and Matsumoto, Y., 2012, “Shock Wave Emission From a Cloud of Bubbles,” Soft Matter, 8(21), pp. 5777–5783. [CrossRef]
Kalumuck, K. M., and Chahine, G. L., 2000, “The Use of Cavitating Jets to Oxidize Organic Compounds in Water,” ASME J. Fluids Eng., 122(3), pp. l465–1470. [CrossRef]
Chahine, G. L., 2014, “Modeling of Cavitation Dynamics and Interaction With Material,” Advanced Experimental and Numerical Techniques for Cavitation Erosion Prediction, Vol. 106, K. H.Kim, G. L.Chahine, J. P.Franc, and A.Karimi, eds., Springer, Dordrecht, The Netherlands, pp. 123–162 [CrossRef].
Plesset, M. S., and Prosperetti, A., 1977, “Bubble Dynamics and Cavitation,” Ann. Rev. Fluid Mech., 9(1), pp. 145–185. [CrossRef]
Blake, J. R., and Gibson, D. C., 1987, “Bubble Dynamics and Cavitation,” Ann. Rev. Fluid Mech., 19(1), pp. 99–123. [CrossRef]
Van Wijngaarden, L., 1964, “On the Collective Collapse of a Large Number of Gas Bubbles in Water,” Proceedings of 11th International Congress of Applied Mechanics, Springer, Berlin, Germany, Aug. 30–Sept. 5, pp. 854–861.
Morch, K. A., 1981, “Cavity Cluster Dynamics and Cavitation Erosion,” ASME Proceeding of the Cavitation and Polyphase Flow Forum, Boulder, CO, June 22–24, pp. 1–10.
Chahine, G. L., 1983, “Cloud Cavitation: Theory,” Proceedings of 14th Symposium on Naval Hydrodynamics, National Academy Press, Ann Arbor, MI, pp. 165–194.
Chahine, G. L., and Liu, H.-L., 1985, “A Singular Perturbation Theory of the Growth of a Bubble Cluster in a Superheated Liquid,” J. Fluid Mech., 156, pp. 257–279. [CrossRef]
Chahine, G. L., and Duraiswami, R., 1992, “Dynamical Interaction in a Multi-Bubble Cloud,” ASME J. Fluids Eng., 114(4), pp. 680–686. [CrossRef]
Esmaeeli, A., and Tryggvason, G., 1998, “Direct Numerical Simulations of Bubbly Flows. Part 1. Low Reynolds Number Array,” J. Fluid Mech., 377, pp. 313–345. [CrossRef]
Lu, J., Biswas, S., and Tryggvason, G., 2006, “A DNS Study of Laminar Bubbly Flows in a Vertical Channel,” Int. J. Multiphase Flow, 32(6), pp. 643–660. [CrossRef]
Seo, J. H., Lele, S. K., and Tryggvason, G., 2010, “Investigation and Modeling of Bubble–Bubble Interaction Effect in Homogeneous Bubbly Flows,” Phys. Fluids, 22(6), p. 063302. [CrossRef]
Crowe, C. T., Troutt, T. R., and Chung, J. N., 1996, “Numerical Models for Two-Phase Turbulent Flows,” Ann. Rev. Fluid Mech., 28, pp. 11–43. [CrossRef]
Balachandar, S., and Eaton, J. K., 2010, “Turbulent Dispersed Multiphase Flow,” Ann. Rev. Fluid Mech., 42(1), pp. 111–133. [CrossRef]
Apte, S. V., Mahesh, K., and Lundgren, T., 2008, “Accounting for Finite-Size Effect in Disperse Two-Phase Flow,” Int. J. Multiphase Flow, 34(3), pp. 260–271. [CrossRef]
Spelt, P. D. M., and Biesheivel, A., 1997, “On the Motion of Gas Bubbles in Homogeneous Isotropic Turbulence,” J. Fluid Mech., 336, pp. 221–244. [CrossRef]
Druzhinin, O. A., and Elghobashi, S. E., 1998, “Direct Numerical Simulations of Bubble-Laden Turbulent Flows Using the Two-Fluid Formulation,” Phys. Fluids, 10(3), pp. 685–697. [CrossRef]
Bunner, B., and Tryggvason, G., 2003, “Effect of Bubble Deformation on the Stability and Properties of Bubbly Flows,” J. Fluid Mech., 495, pp. 77–118. [CrossRef]
Brennen, C. E., 1995, Cavitation and Bubble Dynamics, Oxford University, New York.
Prosperetti, A., Sundaresan, S., Pannala, S., and Zhang, D. Z., 2007, “Segregated Methods for Two-Fluid Models,” Computational Methods for Multiphase Flow, Cambridge University, New York, pp. 320–385.
Ma, J., Oberai, Drew, D. A., Lahey, R. T., Jr., and Moraga, F., 2010, “Two-Fluid Modeling of Bubbly Flows Around Surface Ships Using a Phenomenological Subgrid Air Entrainment Model,” Comput. Fluids, 39(1), pp. 77–86. [CrossRef]
Ma, J., Oberai, A., Hyman, M., Drew, D. A., and Lahey, R. T., Jr., 2011, “Two-Fluid Modeling of Bubbly Flows Around Surface Ships Using a Phenomenological Subgrid Air Entrainment Model,” Comput. Fluids, 52, pp. 50–57. [CrossRef]
Ma, J., Oberai, A., Drew, D. A., and Lahey, R. T., Jr., 2012, “A Two-Way Coupled Polydispersed Two-Fluid Model for the Simulation of Air Entrainment Beneath a Plunging Liquid Jet,” ASME J. Fluids Eng., 134(10), p. 101304. [CrossRef]
Xiang, M., Cheung, S., Yeoh, G., Zhang, W., and Tu, J., 2011, “On the Numerical Study of Bubbly Flow Created by Ventilated Cavity in Vertical Pipe,” Int. J. Multiphase Flow, 37(7), pp. 756–768. [CrossRef]
Raju, R., Singh, S., Hsiao, C.-T., and Chahine, G. L., 2011, “Study of Pressure Wave Propagation in a Two-Phase Bubbly Mixture,” ASME J. Fluids Eng., 133(12), p. 121302. [CrossRef]
Gilmore, F. R., 1952, “The Growth and Collapse of a Spherical Bubble in a Viscous Compressible Liquid,” California Institute of Technology, Hydraulic Laboratory, Report No. 26-4.
Wardlaw, A., Jr., and Luton, J. A., 2000, “Fluid–Structure Interaction for Close-in Explosions,” Shock Vib. J., 7, pp. 265–275. [CrossRef]
Wardlaw, A., Jr., Luton, J. A., Renzi, J. J., and Kiddy, K., 2003, “Fluid–Structure Coupling Methodology for Undersea Weapons,” Fluid Structure Interaction II, WIT Press, Southampton, UK, pp. 251–263.
Ma, J., Singh, S., Hsiao, C.-T., Choi, J.-K., and Chahine, G. L., 2012, “Spherical Bubble Dynamics in a Bubbly Medium,” International Conference on Numerical Methods in Multiscale Flows, Philadelphia, PA, June 12–14.
Ma, J., Hsiao, C.-T., and Chahine, G. L., 2012, “Bubble Dynamics in Bubbly Medium,” 65th Annual Meeting of the APS Division of Fluid Dynamics, San Diego, CA, Nov. 18–20.
Ma, J., Chahine, G. L., and Hsiao, C.-T., “Spherical Bubble Dynamics in a Bubbly Medium Using an Euler–Lagrange Model,” Chem. Eng. Sci., (unpublished).
Hsiao, C.-T., and Chahine, G. L., 2012, “Effect of a Propeller and Gas Diffusion on Bubble Nuclei Distribution in a Liquid,” J. Hydrodyn., Ser. B, 24(6), pp. 809–822. [CrossRef]
Hsiao, C.-T., and Chahine, G. L., 2004, “Prediction of Tip Vortex Cavitation Inception Using Coupled Spherical and Nonspherical Bubble Models and Navier–Stokes Computations,” J. Mar. Sci. Technol., 8(3), pp. 99–108. [CrossRef]
Hsiao, C.-T., and Chahine, G., 2008, “Numerical Study of Cavitation Inception Due to Vortex/Vortex Interaction in a Ducted Propulsor,” J. Ship. Res., 52, pp. 114–123.
Wu, X., Choi, J.-K., Hsiao, C.-T., and Chahine, G. L., 2010, “Bubble Augmented Waterjet Propulsion: Numerical and Experimental Studies,” Proceedings of 28th Symposium on Naval Hydrodynamics, Pasadena, CA, Sept. 12–17, pp. 1–17.
Wu, X., Choi, J.-K., Singh, S., Hsiao, C.-T., and Chahine, G. L., 2012, “Experimental and Numerical Investigation of Bubble Augmented Waterjet Propulsion,” J. Hydrodyn., 24(5), pp. 635–647. [CrossRef]
Singh, S., Choi, J.-K., and Chahine, G. L., 2012, “Optimum Configuration of an Expanding-Contracting-Nozzle for Thrust Enhancement by Bubble Injection,” ASME J. Fluids Eng., 134(1), p. 0113011. [CrossRef]
Hsiao, C. T., Wu, X., Ma, J., and Chahine, G. L., 2013, “Numerical and Experimental Study of Bubble Entrainment Due to a Horizontal Plunging Jet,” Int. Shipbuild. Prog., 60, pp. 435–469.
van Wijngaarden, L., 1968, “On the Equations of Motion for Mixtures of Liquid and Gas Bubbles,” J. Fluid Mech., 33(3), pp. 465–474. [CrossRef]
van Wijngaarden, L., 1972, “One-Dimensional Flow of Liquids Containing Small Gas Bubbles,” Ann. Rev. Fluid Mech., 4, pp. 369–396. [CrossRef]
Commander, K. W., and Prosperetti, A., 1989, “Linear Pressure Waves in Bubbly Liquids: Comparison Between Theory and Experiments,” J. Acoust. Soc. Am., 85(2), pp. 732–746. [CrossRef]
Brennen, C. E., 2005, Fundamentals of Multiphase Flows, Cambridge University, New York.
Awad, M. M., and Muzzchka, Y. S., 2008, “Effective Property Models for Homogeneous Two Phase Flow,” Exp. Therm. Fluid Sci., 33(1), pp. 106–113. [CrossRef]
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, pp. 345–355. [CrossRef]
Darmana, D., Deen, N. G., and Kuipers, J. A. M., 2006, “Parallelization of an Euler–Lagrange Model Using Mixed Domain Decomposition and a Mirror Domain Technique: Application to Disperse Gas–Liquid Two-Phase Flow,” J. Comput. Phys., 220(1), pp. 216–248. [CrossRef]
Shams, E., Finn, J., and Apte, S. V., 2011, “A Numerical Scheme for Euler–Lagrange Simulation of Bubbly Flows in Complex Systems,” Int. J. Num. Methods Fluids, 67(12), pp. 1865–1898. [CrossRef]
Chahine, G. L., 2009, “Numerical Simulation of Bubble Flow Interactions,” J. Hydrodyn., 21(3), pp. 316–332. [CrossRef]
Hsiao, C.-T., Chahine, G. L., and Liu, H.-L., 2003, “Scaling Effects on Prediction of Cavitation Inception in a Line Vortex Flow,” ASME J. Fluids Eng., 125(1), pp. 53–60. [CrossRef]
Johnson, V. E., and Hsieh, T., 1966, “The Influence of the Trajectories of Gas Nuclei on Cavitation Inception,” 6th Symposium on Naval Hydrodynamics, Washington, DC, Sept. 28–Oct. 4, pp. 163–179.
Saffman, P., 1965, “The Lift on a Small Sphere in a Slow Shear Flow,” J. Fluid Mech., 22(2), pp. 385–400. [CrossRef]
Haberman, W. L., and Morton, R. K., 1953, “An Experimental Investigation of the Drag and Shape of Air Bubbles Rising in Various Liquids,” DTMB Report No. 802.
Chorin, A. J., 1967, “A Numerical Method for Solving Incompressible Viscous Flow Problems,” J. Comput. Phys., 2(1), pp. 12–26. [CrossRef]
Roe, P. L., 1981, “Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes,” J. Comput. Phys., 43(2), pp. 357–372. [CrossRef]
Van Leer, B., Thomas, J. L., Roe, P. L., and Newsome, R. W., 1987, “A Comparison of Numerical Flux Formulas for the Euler and Navier–Stokes Equation,” AIAA Paper No. 87-1104-CP. [CrossRef]
Chahine, G. L., Hsiao, C.-T., and Raju, R., 2014, “Scaling of Cavitation Bubble Cloud Dynamics on Propellers,” Advanced Experimental and Numerical Techniques for Cavitation Erosion Prediction, Vol. 106, K. H.Kim, G. L.Chahine, J. P.Franc, and A.Karimi, eds., Springer, Dordrecht, The Netherlands, pp. 345–372.


Grahic Jump Location
Fig. 1

Illustration of the void fraction computation [33]

Grahic Jump Location
Fig. 2

Schematic of a bubble cloud with an equivalent radius A(t) and distributed bubbles of radius rb(t) near a rigid wall. Ox is the axis of symmetry and r is the radial distance from the axis of symmetry. The distance between the center of the cloud and the wall is X(t). AT(t) and AB(t) indicate the distances between the center of the cloud and the top and bottom of the cloud, respectively. The cloud is subjected to an imposed pressure function at its edge, P(t).

Grahic Jump Location
Fig. 3

Time sequences of bubble locations and pressures in a bubble cloud excited by a pressure wave as listed in Table 1: (a) case 1, (b) case 2, and (c) case 3. The bubbles are represented by spheres colored by the gas pressure inside each of them. Lengths are scaled by bubble radii in the top row corresponding to 50 μm. Velocity vectors are represented by arrows. Only velocity vectors larger than 2 m/s extend beyond the bubble representation and can be seen in the snapshots. The bottom row shows only the velocity vectors on the bubbles in each case at t = T.

Grahic Jump Location
Fig. 4

Variations of the total volume of all bubbles in the cloud for the three cases described in Table 1

Grahic Jump Location
Fig. 5

Variations of the outer cloud size for the three cases in Table 1. AB and AT are the distances from the center of the cloud to the top and bottom of the cloud as defined in Fig. 2.

Grahic Jump Location
Fig. 6

Void fraction distributions within the cloud at different times during the first cycle for cases 2 and 3, i.e., without and with a wall, α0 = 5%, ω = 10 kHz. (Note that due to the selected Gaussian spreading radius of A0 (see Eq. (9)) a nonzero value of the void fraction exists here, beyond the cloud radius, i.e., approximately between −2A0 and 2 A0.)

Grahic Jump Location
Fig. 7

Isosurface of 5% void fraction at t = 0.9 T for case 3 (presence of a wall). The arrows denote the translational velocities of the bubbles located at different locations in the cloud.

Grahic Jump Location
Fig. 8

Comparison of the size variation of bubbles at different cloud locations for the three different test cases described in Table 1 (a) ω = 100 kHz, no wall; (b) ω = 10 kHz, no wall; and (c) ω = 10 kHz, with wall

Grahic Jump Location
Fig. 9

Comparison of the pressures at different cloud locations for the three test cases described in Table 1: (a) ω = 100 kHz, no wall; (b) ω = 10 kHz, no wall; and (c) ω = 10 kHz, with wall

Grahic Jump Location
Fig. 10

Size variation of bubbles on the axis of symmetry at seven radial locations in the cloud for case 3 of α0 = 5%, ω = 10 kHz, with wall. The numbers 1–7 on each line indicate the order in which each bubble starts to collapse. Negative values of x correspond to bubble locations in the lower half of the cloud closer to the wall.

Grahic Jump Location
Fig. 11

Variations of the axial (along x in Fig. 2) component of bubble translation velocity for bubble on the axis of symmetry at seven radial locations in the cloud for case 3 of α0 = 5%, ω = 10 kHz, with wall

Grahic Jump Location
Fig. 12

Time evolution of the pressure at different radial locations on the wall for case 3 of α0 = 5%, ω = 10 kHz, with wall




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In