scholarly journals Non-spherical bubble dynamics in a compressible liquid. Part 2. Acoustic standing wave

2011 ◽  
Vol 679 ◽  
pp. 559-581 ◽  
Author(s):  
Q. X. WANG ◽  
J. R. BLAKE

This paper investigates the behaviour of a non-spherical cavitation bubble in an acoustic standing wave. The study has important applications to sonochemistry and in understanding features of therapeutic ultrasound in the megahertz range, extending our understanding of bubble behaviour in the highly nonlinear regime where jet and toroidal bubble formation may be important. The theory developed herein represents a further development of the material presented in Part 1 of this paper (Wang & Blake, J. Fluid Mech. vol. 659, 2010, pp. 191–224) to a standing wave, including repeated topological changes from a singly to a multiply connected bubble. The fluid mechanics is assumed to be compressible potential flow. Matched asymptotic expansions for an inner and outer flow are performed to second order in terms of a small parameter, the bubble-wall Mach number, leading to weakly compressible flow formulation of the problem. The method allows the development of a computational model for non-spherical bubbles by using a modified boundary-integral method. The computations show that the bubble remains approximately of a spherical shape when the acoustic pressure is small or is initiated at the node or antinode of the acoustic pressure field. When initiated between the node and antinode at higher acoustic pressures, the bubble loses its spherical shape at the end of the collapse phase after only a few oscillations. A high-speed liquid bubble jet forms and is directed towards the node, impacting the opposite bubble surface and penetrating through the bubble to form a toroidal bubble. The bubble first rebounds in a toroidal form but re-combines to a singly connected bubble, expanding continuously and gradually returning to a near spherical shape. These processes are repeated in the next oscillation.

1988 ◽  
Vol 110 (3) ◽  
pp. 408-413 ◽  
Author(s):  
L. J. Ghosn

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.


2015 ◽  
Vol 5 (5) ◽  
pp. 20150048 ◽  
Author(s):  
Qianxi Wang ◽  
Wenke Liu ◽  
A. M. Zhang ◽  
Yi Sui

A bubble initiated near a rigid boundary may be almost in contact with the boundary because of its expansion and migration to the boundary, where a thin layer of water forms between the bubble and the boundary thereafter. This phenomenon is modelled using the weakly compressible theory coupled with the boundary integral method. The wall effects are modelled using the imaging method. The numerical instabilities caused by the near contact of the bubble surface with the boundary are handled by removing a thin layer of water between them and joining the bubble surface with its image to the boundary. Our computations correlate well with experiments for both the first and second cycles of oscillation. The time history of the energy of a bubble system follows a step function, reducing rapidly and significantly because of emission of shock waves at inception of a bubble and at the end of collapse but remaining approximately constant for the rest of the time. The bubble starts being in near contact with the boundary during the first cycle of oscillation when the dimensionless stand-off distance γ = s / R m < 1, where s is the distance of the initial bubble centre from the boundary and R m is the maximum bubble radius. This leads to (i) the direct impact of a high-speed liquid jet on the boundary once it penetrates through the bubble, (ii) the direct contact of the bubble at high temperature and high pressure with the boundary, and (iii) the direct impingement of shock waves on the boundary once emitted. These phenomena have clear potential to damage the boundary, which are believed to be part of the mechanisms of cavitation damage.


2007 ◽  
Vol 570 ◽  
pp. 407-429 ◽  
Author(s):  
M. LEE ◽  
E. KLASEBOER ◽  
B. C. KHOO

The formation of a toroidal bubble towards the end of the bubble collapse stage in the neighbourhood of a solid boundary has been successfully studied using the boundary integral method. The further evolution (rebound) of the toroidal bubble is considered with the loss of system energy taken into account. The energy loss is incorporated into a mathematical model by a discontinuous jump in the potential energy at the minimum volume during the short collapse–rebound period accompanying wave emission. This implementation is first tested with the spherically oscillating bubble system using the theoretical Rayleigh–Plesset equation. Excellent agreement with experimental data for the bubble radius evolution up to three oscillation periods is obtained. Secondly, the incorporation of energy loss is tested with the motion of an oscillating bubble system in the neighbourhood of a rigid boundary, in an axisymmetric geometry, using a boundary integral method. Example calculations are presented to demonstrate the possibility of capturing the peculiar entity of a counterjet, which has been reported only in recent experimental studies.


2000 ◽  
Vol 404 ◽  
pp. 151-177 ◽  
Author(s):  
YONGGANG ZHU ◽  
HASAN N. OĞUZ ◽  
ANDREA PROSPERETTI

The process by which a liquid jet falling into a liquid pool entrains air is studied experimentally and theoretically. It is shown that, provided the nozzle from which the jet issues is properly contoured, an undisturbed jet does not entrap air even at relatively high Reynolds numbers. When surface disturbances are generated on the jet by a rapid increase of the liquid flow rate, on the other hand, large air cavities are formed. Their collapse under the action of gravity causes the entrapment of bubbles in the liquid. This sequence of events is recorded with a CCD and a high-speed camera. A boundary-integral method is used to simulate the process numerically with results in good agreement with the observations. An unexpected finding is that the role of the jet is not simply that of conveying the disturbance to the pool surface. Rather, both the observed energy budget and the simulations imply the presence of a mechanism by which part of the jet energy is used in creating the cavity. A hypothesis on the nature of this mechanism is presented.


2016 ◽  
Vol 42 ◽  
pp. 1660157
Author(s):  
SHUAI LI ◽  
SHI-PING WANG ◽  
A-MAN ZHANG

The dynamics of a toroidal bubble splitting near a rigid wall in an inviscid incompressible fluid is studied in this paper. The boundary integral method is adopted to simulate the bubble motion. After the jet impact, the vortex ring model is used to handle the discontinued potential of the toroidal bubble. When the toroidal bubble is splitting, topology changes are made tear the bubble apart. Then, the vortex ring model is extended to multiple vortex rings to simulate the interaction between two toroidal bubbles. A typical case is discussed in this study. Besides, the velocity fields and pressure contours surrounding the bubble are used to illustrate the numerical results. An annular high pressure region is generated at the splitting location, and the maximum pressure may be much higher than the jet impact. More splits may happen after the first split.


1999 ◽  
Vol 380 ◽  
pp. 339-361 ◽  
Author(s):  
R. P. TONG ◽  
W. P. SCHIFFERS ◽  
S. J. SHAW ◽  
J. R. BLAKE ◽  
D. C. EMMONY

Vapour cavities in liquid flows have long been associated with cavitation damage to nearby solid surfaces and it is thought that the final stage of collapse, when a high- speed liquid jet threads the cavity, plays a vital role in this process. The present study investigates this aspect of the motion of laser-generated cavities in a quiescent liquid when the distance (or stand-off) of the point of inception from a rigid boundary is between 0.8 and 1.2 times the maximum radius of the cavity. Numerical simulations using a boundary integral method with an incompressible liquid impact model provide a framework for the interpretation of the experimental results. It is observed that, within the given interval of the stand-off parameter, the peak pressures measured on the boundary at the first collapse of a cavity attain a local minimum, while at the same time there is an increase in the duration of the pressure pulse. This contrasts with a monotonic increase in the peak pressures as the stand-off is reduced, when the cavity inception point is outside the stated interval. This phenomenon is shown to be due to a splash effect which follows the impact of the liquid jet. Three cases are chosen to typify the splash interaction with the free surface of the collapsing cavity: (i) surface reconnection around the liquid jet; (ii) splash impact at the base of the liquid jet; (iii) thin film splash. Hydrodynamic pressures generated following splash impact are found to be much greater than those produced by the jet impact. The combination of splash impact and the emission of shock waves, together with the subsequent re-expansion, drives the flow around the toroidal cavity producing a distinctive double pressure peak.


1998 ◽  
Vol 369 ◽  
pp. 253-272 ◽  
Author(s):  
WILLIAM W. SCHULTZ ◽  
JEAN-MARC VANDEN-BROECK ◽  
LEI JIANG ◽  
MARC PERLIN

We calculate spatially and temporally periodic standing waves using a spectral boundary integral method combined with Newton iteration. When surface tension is neglected, the non-monotonic behaviour of global wave properties agrees with previous computations by Mercer & Roberts (1992). New accurate results near the limiting form of gravity waves are obtained by using a non-uniform node distribution. It is shown that the crest angle is smaller than 90° at the largest calculated crest curvature. When a small amount of surface tension is included, the crest form is changed significantly. It is necessary to include surface tension to numerically reproduce the steep standing waves in Taylor's (1953) experiments. Faraday-wave experiments in a large-aspect-ratio rectangular container agree with our computations. This is the first time such high-amplitude, periodic waves appear to have been observed in laboratory conditions. Ripple formation and temporal symmetry breaking in the experiments are discussed.


1993 ◽  
Vol 254 ◽  
pp. 437-466 ◽  
Author(s):  
J. M. Boulton-Stone ◽  
J. R. Blake

When a small air bubble bursts from an equilibrium position at an air/water interface, a complex motion ensues resulting in the production of a high-speed liquid jet. This free-surface motion following the burst is modelled numerically using a boundary integral method. Jet formation and liquid entrainment rates from jet breakup into drops are calculated and compared with existing experimental evidence. In order to investigate viscous effects, a boundary layer is included in the calculations by employing a time-stepping technique which allows the boundary mesh to remain orthogonal to the surface. This allows an approximation of the vorticity development in the region of boundary-layer separation during jet formation. Calculated values of pressure and energy dissipation rates in the fluid indicate a violent motion, particularly for smaller bubbles. This has important implications for the biological industry where animal cells in bioreactors have been found to be killed by the presence of small bubbles.


1993 ◽  
Vol 251 ◽  
pp. 79-107 ◽  
Author(s):  
J. P. Best

A spectacular feature of transient cavity collapse in the neighbourhood of a rigid boundary is the formation of a high-speed liquid jet that threads the bubble and ultimately impacts upon the side of the bubble nearest to the boundary. The bubble then evolves into some toroidal form, the flow domain being doubly connected. In this work, the motion of the toroidal bubble is computed by connecting the jet tip to the side of the bubble upon which it impacts. This connection is via a cut introduced into the flow domain and across which the potential is discontinuous, the value of this discontinuity being equal to the circulation in the flow. A boundary integral algorithm is developed to account for this geometry and some example computations are presented. Consideration of the pressure field in the fluid has implications for possible damage mechanisms to structures due to nearby cavity collapse.


Sign in / Sign up

Export Citation Format

Share Document