Numerical Study Of Shock-bubble Interaction

Numerical simulations of the interaction between shock waves and bubbles/bubbles groups were carried out in this paper.The interaction was simulated by the front tracking approach and anisotropic mesh adaptive technology respectively,at the meantime,we investigate the development of shock wave system and the evolution process of bubble surface during the interaction.The front tracking approach,which is one of the most important method to achieve sharp resolution of interface or contact discontinuity between distinct materials,is currently used to study the effects of incident Mach number on interaction of shock with an array of side-by-side bubbles after the capture effect of shock wave and interface was verified.Totally 6 numerical scenarios were considered and discussed in the present bubble array simulations with different Atwood number and shock Mach number.The shock behaviors and geometry deformation was presented in details by using numerical schlieren imaging technique.The mechanism of interactions between the wave structures reflected from the bubbles and the way in which these interactions impact on the deformation of these bubbles were revealed extensively.Through the description of vorticity generation and transport process,the internal mechanism of bubble deformation is also revealed.In order to perform a quantitative assessment on the shocked flow field,the volumetric compression effect and the evolution of the circulation were also considered,the results show that the strength of the incident shock has little influence on the approaching time to quasi-steady state and there is no further vorticity deposition over the bubble interface or in the shocked flow field at the later time of evolution.To improve the accuracy of the numerical algorithm in capturing the details of the flow field and reduce the computational scale effectively,the anisotropic mesh adaptive method is developed and used to simulate the interaction between shock waves and bubbles.Firstly,two classical problems of interaction between shock wave and boundary layer and liquid jet are used to verify the identification and refinement of shock wave and interface.After that,the interaction between shock wave and He bubble was calculated.It demonstrates that the adaptive mesh can well capture the main information of the flow field,adjust the mesh over time,and perform adaptive refinement at the location where the flow field changes quickly.By comparing the mesh size and results of the adaptive mesh with the fixed mesh,the advantages of adaptive generation of anisotropic meshes in reducing mesh size and improving solution accuracy are demonstrated.