Abstract

This paper performs numerical investigations on the interaction of shock wave with an ellipsoidal bubble in liquid medium. The governing equations, including the conservative Euler equations and the non-conservative transport equation of the liquid volume fraction, are discretized based on the finite volume method. A tangent of hyperbola for interface capturing (THINC) interface reconstruction scheme is employed for the phasic densities and the liquid volume fraction to maintain the interface sharpness. The major-axis (z-axis) of the bubble is parallel to the incident planar shock wave. Different collapse behaviors are observed for two ellipsoidal geometries, named as the disk-like bubble and rugby-like one. Different collapse patterns of an initial ellipsoidal bubble are presented and the manifestation is that the transverse jets pierce the bubble differently, classified as along the centerline, off-centerline along the circumferential direction or along the meridian line if the aspect ratio is varied. For the disk-like bubble, it presents the strongest collapsing process under certain eccentricity, characterized by the highest water hammer pressure under the same incident shock strength. The second sheeting jet is an important factor that leads to the collapse of the remaining bubbles pierced by the first transverse jet.

Highlights

  • Bubble collapse occurs in a wide variety of physical scenarios, such as cavitation in hydraulic machines,[1] drug delivery and ultrasound lithotripsy in medical therapy,[2] and sonoluminescence.[3]

  • The collapse of an ellipsoidal bubble under 1 GPa shock wave in liquid medium is investigated from numerical aspect

  • The governing equations consist of the conservative Euler equations and the non-conservative liquid volume fraction transport equation

Read more

Summary

INTRODUCTION

Bubble collapse occurs in a wide variety of physical scenarios, such as cavitation in hydraulic machines,[1] drug delivery and ultrasound lithotripsy in medical therapy,[2] and sonoluminescence.[3]. Ding et al.[26] performed numerical simulations of bubble collapse under weak (≤30 MPa) and strong shocks (from 500 to 2000 MPa), but the resolution was not enough to provide detailed flow structures during the bubble collapse. Hawker et al.[22] used the front tracking method in their numerical simulations and studied the bubble collapse under strong shock waves, the wave structures as well as the interface evolution were sufficiently presented in their results. Et al.[28] explored more sophisticated equation of state (EOS) in their numerical studies for the problem of shock bubble interaction and estimated the temperature, density and pressure of different EOS during the jet impact.

Problem description and domain setup
Governing equations
Numerical method
CASE I
First stage
Second stage
Third stage
Parameter studies
CASE II
CONCLUSIONS
Full Text
Published version (Free)

Talk to us

Join us for a 30 min session where you can share your feedback and ask us any queries you have

Schedule a call