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

被引:18
作者
Kamran, K. [1 ]
Rossi, R. [1 ,2 ]
Onate, E. [1 ,2 ]
Idelsohn, S. R. [1 ,3 ]
机构
[1] CIMNE, Barcelona 08034, Spain
[2] Univ Politecn Cataluna, Barcelona, Spain
[3] ICREA, Barcelona, Spain
基金
欧洲研究理事会;
关键词
Lagrangian shock hydrodynamics; Variational multiscale stabilization; Two phase flow; PFEM; Bubble implosion; FINITE-ELEMENT-METHOD; 2-PHASE FLOW; SHOCK HYDRODYNAMICS; DYNAMICS; APPROXIMATION; VISCOSITY; COLLAPSE; SETS;
D O I
10.1016/j.cma.2012.11.018
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
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. (C) 2012 Elsevier B.V. All rights reserved.
引用
收藏
页码:210 / 225
页数:16
相关论文
共 46 条
[1]  
[Anonymous], 2003, SELCUK J APPL MATH
[2]  
[Anonymous], 1993, SIMILARITY DIMENSION, DOI DOI 10.1016/C2013-0-08173-X
[4]   Single-bubble sonoluminescence: Shape stability analysis of collapse dynamics in a semianalytical approach [J].
Bogoyavlenskiy, VA .
PHYSICAL REVIEW E, 2000, 62 (02) :2158-2167
[5]  
Brennan CE., 1995, Cavitation and Bubble Dynamics
[6]   Elimination of artificial grid distortion and hourglass-type motions by means of Lagrangian subzonal masses and pressures [J].
Caramana, EJ ;
Shashkov, MJ .
JOURNAL OF COMPUTATIONAL PHYSICS, 1998, 142 (02) :521-561
[7]   An enriched finite element method and level sets for axisymmetric two-phase flow with surface tension [J].
Chessa, J ;
Belytschko, T .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2003, 58 (13) :2041-2064
[8]   Finite element approximation of the hyperbolic wave equation in mixed form [J].
Codina, Ramon .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2008, 197 (13-16) :1305-1322
[9]   Improving Eulerian two-phase flow finite element approximation with discontinuous gradient pressure shape functions [J].
Coppola-Owen, AH ;
Codina, R .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, 2005, 49 (12) :1287-1304
[10]  
Dadvand P., 2012, COMPUT FLUIDS