Spark-generated bubble collapse near or inside a circular aperture and the ensuing vortex ring and droplet formation

This paper aims to study the oscillation of a sparkgenerated submerged bubble located near or inside a circular aperture made in a flat plate using high-speed visualization technique. In the case of a bubble oscillating near an aperture the initial free surface of the water was set at the bottom surface of the plate. The effects of aperture size and bubblefree surface distance on the bubble behavior as well as on the ensuing droplet dynamics are investigated. It was found that the direction of the bubble reentrant jet was towards the aperture or away from it respectively when the normalized aperture size was smaller or greater than a certain critical value. In addition, a toroidal vortex ring was observed to form, which rotated inwards as it moved away from the aperture. It was also found that if the bubble was incepted at a distance sufficiently away from a supercritical size aperture a single droplet could be produced. In the case of a bub- ble initiated in the middle of a circular aperture submerged just beneath the water free surface, the bubble was found to take the shape of an ellipsoid during its expansion. Then a reentrant jet was initiated and pierced the bubble from its top side.

Level set method for numerical simulation of a cavitation bubble,its growth, collapse and rebound near a rigid wall

A level set method of non-uniform grids is used to simulate the whole evolution of a cavitation bubble, including its growth, collapse and rebound near a rigid wall. Single-phase Navier-Stokes equation in the liquid region is solved by MAC projection algorithm combined with second-order ENO scheme for the advection terms. The moving inter-face is captured by the level set function, and the interface velocity is resolved by ＂one-side＂ velocity extension from the liquid region to the bubble region, complementing the second-order weighted least squares method across the interface and projection inside bubble. The use of non-uniform grid overcomes the difficulty caused by the large computational domain and very small bubble size. The computation is very stable without suffering from large flow-field gradients, and the results are in good agreements with other studies. The bubble interface kinematics, dynamics and its effect on the wall are highlighted, which shows that the code can effectively capture the ＂shock wave＂-like pressure and velocity at jet impact, toroidal bubble, and complicated pressure structure with peak, plateau and valley in the later stage of bubble oscillating.

Numerical study of spike characteristics due to the motions of a non-spherical rebounding bubble

The boundary integral method（BIM） is used to simulate the 3-D gas bubble, generated within the two bubble pulsation periods in proximity to a free surface in an inviscid, incompressible and irrotational flow. The present method is well validated by comparing the calculated shapes of the bubble and the free surface with both the experimental results and the numerical ones obtained by the Axisymmetric BIM code. The expansion, the collapse of the gas bubble and the further evolution of the rebounding non-spherical bubble are simulated. The various variation patterns of the free surface spike and the bubble centroid for different standoff distances, the buoyancy parameters and the strength parameters are obtained to reveal the nonlinear interaction between the bubble and the free surface. The amplitude of the second maximum bubble volume and the four typical patterns of the bubble jet and the free surface spike are examined in the context of the standoff distance. The large buoyancy is used to elevate the spray dome rather than the free surface spike.