A compressible Lagrangian framework for the simulation of the underwater implosion of large air bubbles

A fully Lagrangian compressible numerical framework for the simulation of underwater implosion of a large air bubble is presented. Both air and water are considered compressible and the equations for the Lagrangian shock hydrodynamics are stabilized via a variationally consistent multiscale method. A nodally perfect matched definition of the interface is used and then the kinetic variables, pressure and density, are duplicated at the interface level. An adaptive mesh generation procedure, which respects the interface connectivities, is applied to provide enough refinement at the interface level. This framework is verified by several benchmarks which evaluate the behavior of the numerical scheme for severe compression and expansion cases. This model is then used to simulate the underwater implosion of a large cylindrical bubble, with a size in the order of cm. We observe that the conditions within the bubble are nearly uniform until the converging pressure wave is strong enough to create very large pressures near the center of the bubble. These bubble dynamics occur on very small spatial (0.3 mm), and time (0.1 ms) scales. During the final stage of the collapse Rayleigh–Taylor instabilities appear at the interface and then disappear when the rebounce starts. At the end of the rebounce phase the bubble radius reaches 50% of its initial value and the bubble recover its circular shape. It is when the second collapse starts, with higher mode shape instabilities excited at the bubble interface, that leads to the rupture of the bubble. Several graphs are presented and the pressure pulse detected in the water is compared by experiment.

[1]  Eugenio Oñate,et al.  Derivation of stabilized equations for numerical solution of advective-diffusive transport and fluid flow problems , 1998 .

[2]  William J. Rider,et al.  A conservative nodal variational multiscale method for Lagrangian shock hydrodynamics , 2010 .

[3]  Eugenio Oñate,et al.  On the analysis of heterogeneous fluids with jumps in the viscosity using a discontinuous pressure field , 2010 .

[4]  R. Codina Finite element approximation of the hyperbolic wave equation in mixed form , 2008 .

[5]  W. F. Noh Errors for calculations of strong shocks using an artificial viscosity and artificial heat flux , 1985 .

[6]  Andrea Prosperetti,et al.  The equation of bubble dynamics in a compressible liquid , 1987 .

[7]  Christopher Earls Brennen Cavitation and Bubble Dynamics: Cavitation Bubble Collapse , 2013 .

[8]  Guglielmo Scovazzi,et al.  A geometrically-conservative, synchronized, flux-corrected remap for arbitrary Lagrangian-Eulerian computations with nodal finite elements , 2011, J. Comput. Phys..

[9]  E. Oñate,et al.  The particle finite element method. An overview , 2004 .

[10]  Joseph B. Keller,et al.  Damping of Underwater Explosion Bubble Oscillations , 1956 .

[11]  Riccardo Rossi,et al.  Migration of a generic multi-physics framework to HPC environments , 2013 .

[12]  T. Belytschko,et al.  An enriched finite element method and level sets for axisymmetric two‐phase flow with surface tension , 2003 .

[13]  R. Codina,et al.  Improving Eulerian two‐phase flow finite element approximation with discontinuous gradient pressure shape functions , 2005 .

[14]  Xiaoliang Wan,et al.  Comput. Methods Appl. Mech. Engrg. , 2010 .

[15]  Eugenio Oñate,et al.  Multi-fluid flows with the Particle Finite Element Method , 2009 .

[16]  Kenneth E. Jansen,et al.  Hydrodynamic simulation of air bubble implosion using a level set approach , 2006, J. Comput. Phys..

[17]  Detlef Lohse,et al.  Analysis of Rayleigh–Plesset dynamics for sonoluminescing bubbles , 1998, Journal of Fluid Mechanics.

[18]  S. Turner Underwater implosion of glass spheres. , 2007, The Journal of the Acoustical Society of America.

[19]  Jonathan Richard Shewchuk,et al.  Delaunay refinement algorithms for triangular mesh generation , 2002, Comput. Geom..

[20]  M. Schoenberg,et al.  Acoustic signatures from deep water implosions of spherical cavities , 1976 .

[21]  Eugenio Oñate,et al.  Possibilities of Finite Calculus in Computational Mechanics , 2001 .

[22]  Veselin Dobrev,et al.  Curvilinear finite elements for Lagrangian hydrodynamics , 2011 .

[23]  Thomas J. R. Hughes,et al.  A comparative study of different sets of variables for solving compressible and incompressible flows , 1998 .

[24]  Holland,et al.  Decay of large underwater bubble oscillations , 2000, The Journal of the Acoustical Society of America.

[25]  L. Rayleigh VIII. On the pressure developed in a liquid during the collapse of a spherical cavity , 1917 .

[26]  Eugenio Oñate,et al.  An Object-oriented Environment for Developing Finite Element Codes for Multi-disciplinary Applications , 2010 .

[27]  Charbel Farhat,et al.  A higher-order generalized ghost fluid method for the poor for the three-dimensional two-phase flow computation of underwater implosions , 2008, J. Comput. Phys..

[28]  T. Hughes,et al.  The variational multiscale method—a paradigm for computational mechanics , 1998 .

[29]  A. J. Barlow,et al.  A compatible finite element multi‐material ALE hydrodynamics algorithm , 2008 .

[30]  L Howarth Similarity and Dimensional Methods in Mechanics , 1960 .

[31]  Thomas J. R. Hughes,et al.  Stabilized shock hydrodynamics: I. A Lagrangian method , 2007 .

[32]  Ted Taylor Los Alamos National Laboratory , 2005 .

[33]  T. Hughes Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods , 1995 .

[34]  Eugenio Oñate,et al.  The particle finite element method: a powerful tool to solve incompressible flows with free‐surfaces and breaking waves , 2004 .

[35]  Bogoyavlenskiy Single-bubble sonoluminescence: shape stability analysis of collapse dynamics in a semianalytical approach , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[36]  A. Prosperetti,et al.  Bubble Dynamics and Cavitation , 1977 .

[37]  R. D. Richtmyer,et al.  A Method for the Numerical Calculation of Hydrodynamic Shocks , 1950 .

[38]  Andrea Prosperetti,et al.  Bubble dynamics in a compressible liquid. Part 1. First-order theory , 1986, Journal of Fluid Mechanics.

[39]  Guglielmo Scovazzi,et al.  Lagrangian shock hydrodynamics on tetrahedral meshes: A stable and accurate variational multiscale approach , 2012, J. Comput. Phys..

[40]  M. Shashkov,et al.  Elimination of Artificial Grid Distortion and Hourglass-Type Motions by Means of Lagrangian Subzonal Masses and Pressures , 1998 .

[41]  Guglielmo Scovazzi,et al.  Algebraic Flux Correction and Geometric Conservation in ALE Computations , 2012 .

[42]  Mikhail J. Shashkov,et al.  A Compatible Lagrangian Hydrodynamics Algorithm for Unstructured Grids , 2003 .

[43]  Robert Hickling,et al.  Collapse and rebound of a spherical bubble in water , 1964 .

[44]  William J. Rider,et al.  Revisiting Wall Heating , 2000 .

[45]  Guglielmo Scovazzi,et al.  Stabilized shock hydrodynamics: II. Design and physical interpretation of the SUPG operator for Lagrangian computations☆ , 2007 .