Numerical modelling of an acoustically-driven bubble collapse near a solid boundary

被引:17
作者
Boyd, Bradley [1 ]
Becker, Sid [1 ]
机构
[1] Univ Canterbury, Dept Mech Engn, Christchurch, New Zealand
关键词
WENO; ultrasound-induced collapse; compressible multiphase flow; acoustic cavitation; LOW-FREQUENCY SONOPHORESIS; SHOCK-INDUCED COLLAPSE; CAVITATION BUBBLES; COMPRESSIBLE FLOWS; FLUID METHOD; TRANSIENT CAVITIES; MULTIPHASE FLOWS; RIGID BOUNDARY; WENO SCHEMES; FINAL STAGE;
D O I
10.1088/1873-7005/aad58a
中图分类号
O3 [力学];
学科分类号
08 ; 0801 ;
摘要
The acoustically-driven collapse of a bubble results from the oscillating position of an ultrasound transducer face in a liquid medium. Existing fully compressible models of bubble collapse have been applied to represent the Rayleigh collapse and the shock-induced collapse. These applications may not adequately represent the conditions associated with the acoustically-driven collapse. The current work presents a fully compressible model that is the first to capture the collapse of a bubble set in a liquid medium subjected to an ultrasound transducer. The oscillating transducer face is represented by an immersed moving reflective boundary. The flow is simulated using a conservative interface capturing method, which includes the use of a high-order WENO reconstruction, a maximum-principle-satisfying and positivity-preserving limiter, and the HLLC approximate Riemann flux. The numerical method is verified by quantitative comparison to the benchmark shock-bubble problem. A simulation is conducted of an acoustically-driven collapse of an initially spherical bubble near a solid boundary (wall). The radius of the bubble at the beginning of the collapse, R-max = 40.56 mm, is given by the Rayleigh-Plesset growth from a R-0 = 10 mu m bubble. The distance from the centre of the bubble to the wall (standoff distance) is S = 1.1R(max). The acoustic field resulted from a transducer face oscillating with a frequency of 30 kHz and a displacement amplitude of 0.4174 mu m, with a subsequent acoustic pressure amplitude of p(A) approximate to 120 kPa. The collapse shape is found to be qualitatively consistent with previous work and a jet is observed after 4.7 mu s. The resulting maximum pressure at the wall is about 6.7 MPa. We compare the wall pressures of the acoustically-driven collapse to the Rayleigh collapse. Compared to the acoustically-driven collapse model, the Rayleigh collapse can easily overestimate or underestimate the resulting maximum wall pressure depending on the choice of the initial pressure distribution.
引用
收藏
页数:32
相关论文
共 73 条
[1]   Computations of compressible multifluids [J].
Abgrall, R ;
Karni, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 2001, 169 (02) :594-623
[2]   Numerical analysis of Rayleigh-Plesset equation for cavitating water jets [J].
Alehossein, H. ;
Qin, Z. .
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, 2007, 72 (07) :780-807
[3]   A five-equation model for the simulation of interfaces between compressible fluids [J].
Allaire, G ;
Clerc, S ;
Kokh, S .
JOURNAL OF COMPUTATIONAL PHYSICS, 2002, 181 (02) :577-616
[4]  
[Anonymous], 1998, Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws, DOI 10.1007/BFb0096355
[5]  
[Anonymous], THESIS
[6]   THE FORMATION OF TOROIDAL BUBBLES UPON THE COLLAPSE OF TRANSIENT CAVITIES [J].
BEST, JP .
JOURNAL OF FLUID MECHANICS, 1993, 251 :79-107
[7]   TRANSIENT CAVITIES NEAR BOUNDARIES .1. RIGID BOUNDARY [J].
BLAKE, JR ;
TAIB, BB ;
DOHERTY, G .
JOURNAL OF FLUID MECHANICS, 1986, 170 :479-497
[8]  
Brennen C.E., 1995, CAVITATION BUBBLE DY, V44
[9]   The final stage of the collapse of a cavitation bubble close to a rigid boundary [J].
Brujan, EA ;
Keen, GS ;
Vogel, A ;
Blake, JR .
PHYSICS OF FLUIDS, 2002, 14 (01) :85-92
[10]  
Chahine G.L., 2014, ADV EXPT NUMERICAL T, P123