Observations of the collapses and rebounds of millimeter-sized lithotripsy bubbles

Bubbles excited by lithotripter shock waves undergo a prolonged growth followed by an inertial collapse and rebounds. In addition to the relevance for clinical lithotripsy treatments, such bubbles can be used to study the mechanics of inertial collapses. In particular, both phase change and diffusion among vapor and noncondensable gas molecules inside the bubble are known to alter the collapse dynamics of individual bubbles. Accordingly, the role of heat and mass transport during inertial collapses is explored by experimentally observing the collapses and rebounds of lithotripsy bubbles for water temperatures ranging from 20 to 60 °C and dissolved gas concentrations from 10 to 85% of saturation. Bubble responses were characterized through high-speed photography and acoustic measurements that identified the timing of individual bubble collapses. Maximum bubble diameters before and after collapse were estimated and the corresponding ratio of volumes was used to estimate the fraction of energy retained by the bubble through collapse. The rebounds demonstrated statistically significant dependencies on both dissolved gas concentration and temperature. In many observations, liquid jets indicating asymmetric bubble collapses were visible. Bubble rebounds were sensitive to these asymmetries primarily for water conditions corresponding to the most dissipative collapses.     

Cavitation bubble behavior and bubble–shock wave interaction near a gelatin surface as a study of in vivo bubble dynamics

The collapse of a single cavitation bubble near a gelatin surface, and the interaction of an air bubble attached to a gelatin surface with a shock wave, were investigated. These events permitted the study of the behavior of in vivo cavitation bubbles and the subsequent tissue damage mechanism during intraocular surgery, intracorporeal and extracorporeal shock wave lithotripsy. Results were obtained with high-speed framing photography. The cavitation bubbles near the gelatin surface did not produce significant liquid jets directed at the surface, and tended to migrate away from it. The period of the motion of a cavitation bubble near the gelatin surface was longer than that of twice the Rayleigh's collapse time for a wide range of relative distance, L/R_max, excepting for very small L/R_max values (L was the stand-off distance between the gelatin surface and the laser focus position, and R_max was the maximum bubble radius). The interaction of an air bubble with a shock wave yielded a liquid jet inside the bubble, penetrating into the gelatin surface. The liquid jet had the potential to damage the gelatin. The results predicted that cavitation-bubble-induced tissue damage was closely related to the oscillatory bubble motion, the subsequent mechanical tissue displacement, and the liquid jet penetration generated by the interaction of the remaining gas bubbles with subsequent shock waves. The characteristic bubble motion and liquid jet formation depended on the tissue's mechanical properties, resulting in different damage mechanisms from those observed on hard materials.     

Numerical investigation of shock-induced bubble collapse dynamics and fluid–solid interactions during shock-wave lithotripsy

In this paper we investigate the bubble collapse dynamics under shock-induced loading near soft and rigid bio-materials, during shock wave lithotripsy. A novel numerical framework was developed, that employs a Diffuse Interface Method (DIM) accounting for the interaction across fluid–solid-gas interfaces. For the resolution of the extended variety of length scales, due to the dynamic and fine interfacial structures, an Adaptive Mesh Refinement (AMR) framework for unstructured grids was incorporated. This multi-material multi-scale approach aims to reduce the numerical diffusion and preserve sharp interfaces. The presented numerical framework is validated for cases of bubble dynamics, under high and low ambient pressure ratios, shock-induced collapses, and wave transmission problems across a fluid–solid interface, against theoretical and numerical results. Three different configurations of shock-induced collapse applications near a kidney stone and soft tissue have been simulated for different stand-off distances and bubble attachment configurations. The obtained results reveal the detailed collapse dynamics, jet formation, solid deformation, rebound, primary and secondary shock wave emissions, and secondary collapse that govern the near-solid collapse and penetration mechanisms. Significant correlations of the problem configuration to the overall collapse mechanisms were found, stemming from the contact angle/attachment of the bubble and from the properties of solid material. In general, bubbles with their center closer to the kidney stone surface produce more violent collapses. For the soft tissue, the bubble movement prior to the collapse is of great importance as new structures can emerge which can trap the liquid jet into induced crevices. Finally, the tissue penetration is examined for these cases and a novel tension-driven tissue injury mechanism is elucidated, emanating from the complex interaction of the bubble/tissue interaction during the secondary collapse phase of an entrapped bubble in an induced crevice with the liquid jet.     

Study on shock wave-induced cavitation bubbles dissolution process

This study investigated dissolution processes of cavitation bubbles generated during in vivo shock wave(SW)-induced treatments. Both active cavitation detection(ACD) and the B-mode imaging technique were applied to measure the dissolution procedure of bi Spheres contrast agent bubbles by in vitro experiments. Besides, the simulation of SW-induced cavitation bubbles dissolution behaviors detected by the B-mode imaging system during in vivo SW treatments, including extracorporeal shock wave lithotripsy(ESWL) and extracorporeal shock wave therapy(ESWT), were carried out based on calculating the integrated scattering cross-section of dissolving gas bubbles with employing gas bubble dissolution equations and Gaussian bubble size distribution. The results showed that(i) B-mode imaging technology is an effective tool to monitor the temporal evolution of cavitation bubbles dissolution procedures after the SW pulses ceased, which is important for evaluation and controlling the cavitation activity generated during subsequent SW treatments within a treatment period;(ii) the characteristics of the bubbles, such as the bubble size distribution and gas diffusion, can be estimated by simulating the experimental data properly.     

A dual passive cavitation detector for localized detection of lithotripsy-induced cavitation in vitro

A passive cavitation detector (PCD) identifies cavitation events by sensing acoustic emissions generated by the collapse of bubbles. In this work, a dual passive cavitation detector (dual PCD), consisting of a pair of orthogonal confocal receivers, is described for use in shock wave lithotripsy. Cavitation events are detected by both receivers and can be localized to within 5 mm by the nature of the small intersecting volume of the focal areas of the two receivers in association with a coincidence detection algorithm. A calibration technique, based on the impulse response of the transducer, was employed to estimate radiated pressures at collapse near the bubble. Results are presented for the in vitro cavitation fields of both a clinical and a research electrohydraulic lithotripter. The measured lifetime of the primary growth-and-collapse of the cavitation bubbles increased from 180 to 420 microseconds as the power setting was increased from 12 to 24 kV. The measured lifetime compared well with calculations based on the Gilmore-Akulichev formulation for bubble dynamics. The radiated acoustic pressure 10 mm from the collapsing cavitation bubble was measured to vary from 4 to 16 MPa with increasing power setting; although the trends agreed with calculations, the predicted values were four times larger than measured values. The axial length of the cavitation field correlated well with the 6-dB region of the acoustic field. However, the width of the cavitation field (10 mm) was significantly narrower than the acoustic field (25 mm) as bubbles appeared to be drawn to the acoustic axis during the collapse. The dual PCD also detected signals from "rebounds," secondary and tertiary growth-and-collapse cycles. The measured rebound time did not agree with calculations from the single-bubble model. The rebounds could be fitted to a Rayleigh collapse model by considering the entire bubble cloud as an effective single bubble. The results from the dual PCD agreed well with images from high-speed photography. The results indicate that single-bubble theory is sufficient to model lithotripsy cavitation dynamics up to time of the main collapse, but that upon collapse bubble cloud dynamics becomes important.     

Experimental investigation of the collapse of laser-generated cavitation bubbles near a solid boundary

The oscillation of a laser-generated single cavitation bubble near a solid boundary is investigated by a fiber-optic diagnostic technique based on optical beam deflection (OBD). The maximum bubble radii and collapse time for each oscillation cycle are determined from a sequence of bubble oscillations. Furthermore, by combining the revised Rayleigh theory, the prolongation factor κ at different dimensionless parameter γ (γ=L/R_max, where R_max is the maximum bubble radius and L is the distance of a cavity inception point from a boundary) is obtained. In addition, the prolongation factor of the collapse time versus laser energy is also derived, which are valuable in the fields of hydraulic cavitation, laser lithotripsy and laser ophthalmology.     

Modified shock waves for extracorporeal shock wave lithotripsy: a simulation based on the Gilmore formulation

Extracorporeal shock wave lithotripsy (SWL) is a reliable therapy for the treatment of urolithiasis. Nevertheless, improvements to enhance stone fragmentation and reduce tissue damage are still needed. During SWL, cavitation is one of the most important stone fragmentation mechanisms. Bubbles with a diameter between about 7 and 55 μm have been reported to expand and collapse after shock wave passage, forming liquid microjets at velocities of up to 400 m/s that contribute to the pulverization of renal calculi. Several authors have reported that the fragmentation efficiency may be improved by using tandem shock waves. Tandem SWL is based on the fact that the collapse of a bubble can be intensified if a second shock wave arrives tenths or even a few hundredths of microseconds before its collapse. The object of this study is to determine if tandem pulses consisting of a conventional shock wave (estimated rise time between 1 and 20 ns), followed by a slower second pressure profile (0.8 μs rise time), have advantages over conventional tandem SWL. The Gilmore equation was used to simulate the influence of the modified pressure field on the dynamics of a single bubble immersed in water and compare the results with the behavior of the same bubble subjected to tandem shock waves. The influence of the delay between pulses on the dynamics of the collapsing bubble was also studied for both conventional and modified tandem waves. For a bubble of 0.07 mm, our results indicate that the modified pressure profile enhances cavitation compared to conventional tandem waves at a wide range of delays (10-280 μs). According to this, the proposed pressure profile could be more efficient for SWL than conventional tandem shock waves. Similar results were obtained for a ten times smaller bubble.     

A reduced-order, single-bubble cavitation model with applications to therapeutic ultrasound

Cavitation often occurs in therapeutic applications of medical ultrasound such as shock-wave lithotripsy (SWL) and high-intensity focused ultrasound (HIFU). Because cavitation bubbles can affect an intended treatment, it is important to understand the dynamics of bubbles in this context. The relevant context includes very high acoustic pressures and frequencies as well as elevated temperatures. Relative to much of the prior research on cavitation and bubble dynamics, such conditions are unique. To address the relevant physics, a reduced-order model of a single, spherical bubble is proposed that incorporates phase change at the liquid-gas interface as well as heat and mass transport in both phases. Based on the energy lost during the inertial collapse and rebound of a millimeter-sized bubble, experimental observations were used to tune and test model predictions. In addition, benchmarks from the published literature were used to assess various aspects of model performance. Benchmark comparisons demonstrate that the model captures the basic physics of phase change and diffusive transport, while it is quantitatively sensitive to specific model assumptions and implementation details. Given its performance and numerical stability, the model can be used to explore bubble behaviors across a broad parameter space relevant to therapeutic ultrasound.     

Effects of surface tension on the dynamics of a single micro bubble near a rigid wall in an ultrasonic field

Acoustic cavitation is a very important hydrodynamic phenomenon, and is often implicated in a myriad of industrial, medical, and daily living applications. In these applications, the effect mechanism of liquid surface tension on improving the efficiency of acoustic cavitation is a crucial concern for researchers. In this study, the effects of liquid surface tension on the dynamics of an ultrasonic driven bubble near a rigid wall, which could be the main mechanism of efficiency improvement in the applications of acoustic cavitation, were investigated at the microscale level. A synchronous high-speed microscopic imaging method was used to clearly record the temporary evolution of single acoustic cavitation bubble in the liquids with different surface tension. Meanwhile, the bubble dynamic characteristics, such as the position and time of bubble collapse, the size and stability of the bubbles, the speed of bubble boundaries and the micro-jets, were analyzed and compared. In the case of the single bubbles near a rigid wall, it was found that low surface tension reduces the stability of the bubbles in the liquid medium. Meanwhile, the bubbles collapse earlier and farther from the rigid wall in the liquids with lower surface tension. In addition, the surface tension has no significant influence on the speed of the first micro-jet, but it can substantially increase the speed of second and the third micro-jets after the first collapse of the bubble. These effects of liquid surface tension on the bubble dynamics can explain the mechanism of surfactants in numerous fields of acoustic cavitation for facilitating its optimization and application.     

Oscillation characteristics of a laser-induced cavitation bubble in water at different temperatures

Comprehensive numerical and experimental analyses of the effect of temperature on cavitation oscillations are performed. In the experimental study, the oscillation of a laser-generated single cavitation bubble near a rigid boundary is obtained using a fiber-optic diagnostic technique based on optical beam detection (OBD). The maximum and minimum bubble radii as well as the oscillation times for each oscillation cycle are determined according to the characteristic signals. And cavitation bubble tests are performed using water at different temperatures, and its temperature ranges from freezing point (0 °C) to near boiling. Furthermore, a modified Rayleigh-Plesset equation is derived for calculating the temporal development of the bubble radius at different temperatures. Both the experimental and the numerical results show that the maximum bubble radius and bubble lifetime both increase as temperature increases. The mechanism behind it has also been discussed.     

枝晶生长和气泡形成的数值模拟

本文建立了二维的格子玻尔兹曼方法-元胞自动机(lattice Boltzmann method-cellular automaton, LBM-CA)耦合模型, 对凝固过程中枝晶生长和气泡形成进行模拟研究. 本模型采用CA方法模拟枝晶的生长, 根据界面溶质平衡法计算枝晶生长的驱动力. 采用基于Shan-Chen多相流的LBM模拟气泡在液相中的生长和运动. 在LBM-CA的耦合模型中包含了固-液-气三相之间的相互作用. 应用Laplace定理和模拟气-液-固三相之间的润湿现象对模型进行了验证. 应用所建立的LBM-CA耦合模型模拟研究了气-液相互作用系数对单气泡生长的影响. 发现单气泡的生长速度和平衡半径随气-液相互作用系数的增大而增大. 定向凝固过程中枝晶和气泡生长的模拟结果再现了枝晶的择优生长、 气泡的优先形核位置、气泡的长大、合并、在枝晶间受挤变形以及在液相通道中的运动等物理现象, 与实验结果符合良好. 此外, 初始气体含量越高, 凝固结束时气泡的体积分数也相对较高. 本模型的模拟结果可以揭示在凝固过程中气泡形核、 生长和运动演化以及与枝晶生长相互作用的物理机理.     

Optical and acoustic monitoring of bubble cloud dynamics at a tissue-fluid interface in ultrasound tissue erosion

Short, high-intensity ultrasound pulses have the ability to achieve localized, clearly demarcated erosion in soft tissue at a tissue-fluid interface. The primary mechanism for ultrasound tissue erosion is believed to be acoustic cavitation. To monitor the cavitating bubble cloud generated at a tissue-fluid interface, an optical attenuation method was used to record the intensity loss of transmitted light through bubbles. Optical attenuation was only detected when a bubble cloud was seen using high speed imaging. The light attenuation signals correlated well with a temporally changing acoustic backscatter which is an excellent indicator for tissue erosion. This correlation provides additional evidence that the cavitating bubble cloud is essential for ultrasound tissue erosion. The bubble cloud collapse cycle and bubble dissolution time were studied using the optical attenuation signals. The collapse cycle of the bubble cloud generated by a high intensity ultrasound pulse of 4-14 micros was approximately 40-300 micros depending on the acoustic parameters. The dissolution time of the residual bubbles was tens of ms long. This study of bubble dynamics may provide further insight into previous ultrasound tissue erosion results.     

纳米尺度下气泡核化生长的分子动力学研究

采用分子动力学方法模拟纳米尺度下液体在固体壁面上发生核化沸腾的过程,主要研究壁面浸润性对气泡初始核化过程和气泡生长速率的影响以及固-液界面效应在液体核化沸腾的能量传递过程中所起到的作用.研究结果发现：壁面浸润性越强,气泡在固壁处越容易核化.该结果与经典核化理论中“疏水壁面易于产生气泡”的现象产生了明显的区别.其根本原因是在纳米尺度下,固-液界面热阻效应不能被忽略.一方面,在相同的壁温下,通过增强固-液相互作用,可以显著降低界面热阻,使得热量传递效率提高,导致靠近壁面处的流体温度升高,气泡核化等待时间缩短,有利于液体沸腾核化.另一方面,气泡的生长速率随着壁面浸润性的增强而明显升高.当气泡体积生长到一定程度时,会在壁面处形成气膜,从而导致壁面传热性能恶化.因此,通过壁面的热流密度呈现出先增大后减小的规律.     

垂直低温两相流管底液氮汽泡上升速度的实验研究

本文对垂直低温两相流管底部液氮汽泡的上升速度运用高速摄像机进行了可视化实验研究,对所采集的实时图像进行处理.在Cole和Mendelson汽泡上升速度经验公式的基础上,通过对液氮汽泡的上升速度的分析,提出了圆管管路底部液氮汽泡的上升速度的拟合公式,分析研究其变化规律,为进一步探求圆管内液氮弹状汽泡的生成机理提供依据.     

Physical mechanisms of the therapeutic effect of ultrasound (a review)

Therapeutic ultrasound is an emerging field with many medical applications. High intensity focused ultrasound (HIFU) provides the ability to localize the deposition of acoustic energy within the body, which can cause tissue necrosis and hemostasis. Similarly, shock waves from a lithotripter penetrate the body to comminute kidney stones, and transcutaneous ultrasound enhances the transport of chemotherapy agents. New medical applications have required advances in transducer design and advances in numerical and experimental studies of the interaction of sound with biological tissues and fluids. The primary physical mechanism in HIFU is the conversion of acoustic energy into heat, which is often enhanced by nonlinear acoustic propagation and nonlinear scattering from bubbles. Other mechanical effects from ultrasound appear to stimulate an immune response, and bubble dynamics play an important role in lithotripsy and ultrasound-enhanced drug delivery. A dramatic shift to understand and exploit these nonlinear and mechanical mechanisms has occurred over the last few years. Specific challenges remain, such as treatment protocol planning and real-time treatment monitoring. An improved understanding of the physical mechanisms is essential to meet these challenges and to further advance therapeutic ultrasound.     

鼓泡床中超声驻波的模拟及其对气泡的调制机理

采用计算流体动力学（CFD）的方法，数值生成了鼓泡床中一对声换能器以16kHz高频振动引发的超声场。数值计算是基于包括粘性影响的可压缩流体基本守恒方程，并耦合了水的状态方程。模拟结果表明，在本研究所用的几何布置和换能器与时间相关的速度入口边界条件下，反应器中形成了一个稳定的驻波声场；由于波的非线性以及水的粘性，压力波节点呈现出轻微的时间漂移性。模拟结果与前人的实验结果定性吻合。在模拟的声压分布的基础上，分析了驻波声场调制气泡的机理。如比较熟知，气泡在驻波声场作用下或者向压力波节点运动或者向压力波腹点运动，取决于气泡尺寸与共振尺寸的关系。     

Laser-induced microjet: wavelength and pulse duration effects on bubble and jet generation for drug injection

The expansion of the laser-induced bubble is the main mechanism in the developed microjet injector. In this study, Nd:YAG and Er:YAG lasers are used as triggers of the bubble formation. The impact of the laser parameters on the bubble dynamics is studied and the performance of the injector is evaluated. We found that the main cause of the differences in the bubble behavior comes from the pulse duration and wavelength. For Nd:YAG laser, the pulse duration is very short relative to the bubble lifetime making the behavior of the bubble close to that of the cavitation bubble, while in Er:YAG case, the high absorption in the water and long pulse duration change the initial behavior of the bubble making it close to a vapor bubble. The contraction and subsequent rebound are typical for cavitation bubbles in both cases. The results show that the laser-induced microjet injector generates velocity which is sufficient for the drug delivery for both laser beams of different pulse duration. We estimate the typical velocity within 30–80 m/s range and the breakup length to be larger than 1 mm suitable for trans-dermal drug injection.     

A comprehensive numerical analysis of heat and mass transfer phenomenons during cavitation sono-process

The present study treats the effects of mass transport, heat transfer and chemical reactions heat on the bubble dynamics by spanning a range of ambient bubble radii. The thermodynamic behavior of the acoustic bubble was shown for three wave frequencies, 355, 515 and 1000 kHz. The used acoustic amplitude ranges from 1 to 3 atm. It has been demonstrated that the ambient bubble radius, R₀, of the maximal response (i.e., maximal bubble temperature and pressure, T_max and P_max) is shifted toward lower values if the acoustic amplitude (at fixed frequency) or the ultrasonic frequency (at fixed amplitude) are increased. The range of the ambient bubble radius narrows as the ultrasonic frequency increases. Heat exchange at the bubble interface was found to be the most important mechanism within the bubble internal energy balance for acoustic amplitudes lower than 2.5 and 3 atm for ultrasonic frequencies of 355 and 515 kHz, respectively. For acoustic amplitudes greater or equal to 2.5 and 3 atm, corresponding to 355 and 515 kHz, respectively, mass transport mechanism (i.e., evaporation and condensation of water vapor) becomes dominant compared to the other mechanisms. At 1000 kHz, the mechanism of heat transfer persists to be dominant for all the used acoustic amplitudes (from 1 to 3 atm). Practically, all the above observations were maintained for bubbles at and around the optimum bubble radius, whereas no significant impact of the three energetic mechanisms was observed for bubbles of too lower and too higher values of R₀ (limits of the investigated ranges of R₀).     

超声振动珩磨作用下空化泡动力学及影响参数

为了合理利用超声振动珩磨作用下的空化效应,以磨削区单个空化泡为研究对象,考虑珩磨头合成扰动速度和珩磨压力的作用建立了磨削区空化泡的动力学模型。数值模拟了空化泡初始半径,珩磨压力,液体静压力和超声声压幅值对磨削区空化效应的影响。研究表明考虑超声振动珩磨作用时,空化泡膨胀的幅值会受到抑制,其溃灭时间也会缩短,而且较容易出现稳态空化。珩磨压力和液体静压力对磨削区空化主要起抑制作用,超声波声压幅值在一定范围内能够促进磨削区空化效果的提升。本文的研究为进一步理解超声振动珩磨的空化机理提供了理论支持。