Numerical simulation of underwater explosion bubble with a refined interface treatment

With the intermediate flow states predicted by local two phase Riemann problem, the modified ghost fluid method (MGFM) and its variant (rGFM) have been widely employed to resolve the interface condition in the simulation of compressible multi-medium flows. In this work, the drawback of the construction procedure of local two phase Riemann problem in rGFM was investigated in detail, and a refined version of the construction procedure was specially developed to make the simulation of underwater explosion bubbles more accurate and robust. Beside the refined rGFM, the fast and accurate particle level set method was also adopted to achieve a more effective and computationally efficient capture of the evolving multi-medium interfaces during the simulation. To demonstrate the improvement brought by current refinement, several typical numerical examples of underwater explosion bubbles were performed with original rGFM and refined rGFM, respectively. The results indicate that, when compared with original rGFM, numerical oscillations were effectively removed with the proposed refinement. Accordingly, with present refined treatment of interface condition, a more accurate and robust simulation of underwater explosion bubbles was accomplished in this work.

[1]  Nikolaus A. Adams,et al.  A conservative interface method for compressible flows , 2006, J. Comput. Phys..

[2]  T. Saito,et al.  Experimental and numerical studies of underwater shock wave attenuation , 2003 .

[3]  Boo Cheong Khoo,et al.  A Real Ghost Fluid Method for the Simulation of Multimedium Compressible Flow , 2006, SIAM J. Sci. Comput..

[4]  H. S. Udaykumar,et al.  Ghost Fluid Method for Strong Shock Interactions Part 1: Fluid-Fluid Interfaces , 2009 .

[5]  S. Osher,et al.  A Non-oscillatory Eulerian Approach to Interfaces in Multimaterial Flows (the Ghost Fluid Method) , 1999 .

[6]  R. Fedkiw,et al.  USING THE PARTICLE LEVEL SET METHOD AND A SECOND ORDER ACCURATE PRESSURE BOUNDARY CONDITION FOR FREE SURFACE FLOWS , 2003 .

[7]  S. Osher,et al.  Spatially adaptive techniques for level set methods and incompressible flow , 2006 .

[8]  Ronald Fedkiw,et al.  An Isobaric Fix for the Overheating Problem in Multimaterial Compressible Flows , 1999 .

[9]  Cho-Chung Liang,et al.  Shock responses of a surface ship subjected to noncontact underwater explosions , 2006 .

[10]  G. Tryggvason,et al.  A front-tracking method for viscous, incompressible, multi-fluid flows , 1992 .

[11]  H. Udaykumar,et al.  A sharp interface method for high-speed multi-material flows: strong shocks and arbitrary materialpairs , 2011 .

[12]  J. Blake,et al.  Growth and collapse of a vapour cavity near a free surface , 1981, Journal of Fluid Mechanics.

[13]  Ronald Fedkiw,et al.  Regular Article: The Ghost Fluid Method for Deflagration and Detonation Discontinuities , 1999 .

[14]  Takahiro Matsuno,et al.  Numerical investigations of shock–bubble interactions in mercury , 2008 .

[15]  姚熊亮,et al.  Interaction of underwater explosion bubble with complex elastic-plastic structure , 2008 .

[16]  Khoon Seng Yeo,et al.  The simulation of compressible multi-medium flow. I. A new methodology with test applications to 1D gas–gas and gas–water cases , 2001 .

[17]  Direct Numerical Simulations of Interaction of Strong Shock Waves with Nonspherical Gas Bubbles near Glass Boundaries in Mercury , 2009 .

[18]  Frank Losasso,et al.  Simulating water and smoke with an octree data structure , 2004, SIGGRAPH 2004.

[19]  Boo Cheong Khoo,et al.  Ghost fluid method for strong shock impacting on material interface , 2003 .

[20]  Chen Wei-zhong,et al.  Dynamics and measurement of cavitation bubble , 2006 .

[21]  R. Fedkiw,et al.  A numerical method for two-phase flow consisting of separate compressible and incompressible regions , 2000 .

[22]  Yin Lu Young,et al.  Multiphase modeling of dynamic fluid–structure interaction during close‐in explosion , 2008 .

[23]  J. Sethian.  Evolution, implementation, and application of level set and fast marching methods for advancing fronts , 2001 .

[24]  Eric Johnsen,et al.  Implementation of WENO schemes in compressible multicomponent flow problems , 2005, J. Comput. Phys..

[25]  R. Fedkiw,et al.  Coupling an Eulerian fluid calculation to a Lagrangian solid calculation with the ghost fluid method , 2002 .

[26]  Boo Cheong Khoo,et al.  The simulation of compressible multi-medium flow: II. Applications to 2D underwater shock refraction , 2001 .

[27]  Danping Peng,et al.  Weighted ENO Schemes for Hamilton-Jacobi Equations , 1999, SIAM J. Sci. Comput..

[28]  Boo Cheong Khoo,et al.  The ghost fluid method for compressible gas-water simulation , 2005 .

[29]  Ronald Fedkiw,et al.  Level set methods and dynamic implicit surfaces , 2002, Applied mathematical sciences.

[30]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[31]  C. W. Hirt,et al.  Volume of fluid (VOF) method for the dynamics of free boundaries , 1981 .

[32]  Tiegang Liu,et al.  The simulation of cavitating flows induced by underwater shock and free surface interaction , 2007 .

[33]  Evert Klaseboer,et al.  Experimental and numerical investigation of the dynamics of an underwater explosion bubble near a resilient/rigid structure , 2005, Journal of Fluid Mechanics.

[34]  James A. Sethian,et al.  Level Set Methods and Fast Marching Methods: Evolving Interfaces in Computational Geometry, Fluid , 2012 .

[35]  P. Helluy,et al.  Comparison and validation of compressible flow simulations of laser-induced cavitation bubbles , 2009 .

[36]  Frank Losasso,et al.  A fast and accurate semi-Lagrangian particle level set method , 2005 .

[37]  Ian M. Mitchell,et al.  A hybrid particle level set method for improved interface capturing , 2002 .

[38]  J. Sethian,et al.  LEVEL SET METHODS FOR FLUID INTERFACES , 2003 .

[39]  B. Khoo,et al.  Isentropic one-fluid modelling of unsteady cavitating flow , 2004 .

[40]  J. A. Sethian,et al.  Fast Marching Methods , 1999, SIAM Rev..