Simulation of the effects of cavitation and anatomy in the shock path of model lithotripters

A Source Term Approach for Generation of One-way Acoustic Waves in the Euler and Navier-Stokes equations

We derive a volumetric source term for the Euler and Navier-Stokes equations that mimics the generation of unidirectional acoustic waves from an arbitrary smooth surface in three-dimensional space. The model is constructed as a linear combination of monopole and dipole sources in the mass, momentum, and energy equations. The singular source distribution on the surface is regularized on a computational grid by convolution with a smeared Dirac delta function. The source is implemented in the Euler equations using a Cartesian-grid finite-volume WENO scheme, and validated by comparing with analytical solution for unidirectional planar and spherical acoustic waves. Using the scheme, we emulate a spherical piezoelectric transducer and a multi-array transducer to simulate focused ultrasound fields in water. The simulated ultrasound fields show favorable agreement with previous experiments.

Multiphase fluid-solid coupled analysis of shock-bubble-stone interaction in shockwave lithotripsy

A novel multiphase fluid-solid-coupled computational framework is applied to investigate the interaction of a kidney stone immersed in liquid with a lithotripsy shock wave (LSW) and a gas bubble near the stone. The main objective is to elucidate the effects of a bubble in the shock path to the elastic and fracture behaviors of the stone. The computational framework couples a finite volume 2-phase computational fluid dynamics solver with a finite element computational solid dynamics solver. The surface of the stone is represented as a dynamic embedded boundary in the computational fluid dynamics solver. The evolution of the bubble surface is captured by solving the level set equation. The interface conditions at the surfaces of the stone and the bubble are enforced through the construction and solution of local fluid-solid and 2-fluid Riemann problems. This computational framework is first verified for 3 example problems including a 1D multimaterial Riemann problem, a 3D shock-stone interaction problem, and a 3D shock-bubble interaction problem. Next, a series of shock-bubble-stone-coupled simulations are presented. This study suggests that the dynamic response of a bubble to LSW varies dramatically depending on its initial size. Bubbles with an initial radius smaller than a threshold collapse within 1 μs after the passage of LSW, whereas larger bubbles do not. For a typical LSW generated by an electrohydraulic lithotripter (p_max  = 35.0MPa, p_min  =- 10.1MPa), this threshold is approximately 0.12mm. Moreover, this study suggests that a noncollapsing bubble imposes a negative effect on stone fracture as it shields part of the LSW from the stone. On the other hand, a collapsing bubble may promote fracture on the proximal surface of the stone, yet hinder fracture from stone interior.

Review on Lithotripsy and Cavitation in Urinary Stone Therapy

Cavitation is the sudden formation of vapor bubbles or voids in liquid media and occurs after rapid changes in pressure as a consequence of mechanical forces. It is mostly an undesirable phenomenon. Although the elimination of cavitation is a major topic in the study of fluid dynamics, its destructive nature could be exploited for therapeutic applications. Ultrasonic and hydrodynamic sources are two main origins for generating cavitation. The purpose of this review is to give the reader a general idea about the formation of cavitation phenomenon and existing biomedical applications of ultrasonic and hydrodynamic cavitation. Because of the high number of the studies on ultrasound cavitation in the literature, the main focus of this review is placed on the lithotripsy techniques, which have been widely used for the treatment of urinary stones. Accordingly, cavitation phenomenon and its basic concepts are presented in Section II. The significance of the ultrasound cavitation in the urinary stone treatment is discussed in Section III in detail and hydrodynamic cavitation as an important alternative for the ultrasound cavitation is included in Section IV. Finally, side effects of using both ultrasound and hydrodynamic cavitation in biomedical applications are presented in Section V.

Experimentally validated multiphysics computational model of focusing and shock wave formation in an electromagnetic lithotripter

A multiphysics computational model of the focusing of an acoustic pulse and subsequent shock wave formation that occurs during extracorporeal shock wave lithotripsy is presented. In the electromagnetic lithotripter modeled in this work the focusing is achieved via a polystyrene acoustic lens. The transition of the acoustic pulse through the solid lens is modeled by the linear elasticity equations and the subsequent shock wave formation in water is modeled by the Euler equations with a Tait equation of state. Both sets of equations are solved simultaneously in subsets of a single computational domain within the BEARCLAW framework which uses a finite-volume Riemann solver approach. This model is first validated against experimental measurements with a standard (or original) lens design. The model is then used to successfully predict the effects of a lens modification in the form of an annular ring cut. A second model which includes a kidney stone simulant in the domain is also presented. Within the stone the linear elasticity equations incorporate a simple damage model.

Prediction of far-field acoustic emissions from cavitation clouds during shock wave lithotripsy for development of a clinical device

This study presents the key simulation and decision stage of a multi-disciplinary project to develop a hospital device for monitoring the effectiveness of kidney stone fragmentation by shock wave lithotripsy (SWL). The device analyses, in real time, the pressure fields detected by sensors placed on the patient's torso, fields generated by the interaction of the incident shock wave, cavitation, kidney stone and soft tissue. Earlier free-Lagrange simulations of those interactions were restricted (by limited computational resources) to computational domains within a few centimetres of the stone. Later studies estimated the far-field pressures generated when those interactions involved only single bubbles. This study extends the free-Lagrange method to quantify the bubble–bubble interaction as a function of their separation. This, in turn, allowed identification of the validity of using a model of non-interacting bubbles to obtain estimations of the far-field pressures from 1000 bubbles distributed within the focus of the SWL field. Up to this point in the multi-disciplinary project, the design of the clinical device had been led by the simulations. This study records the decision point when the project's direction had to be led by far more costly clinical trials instead of the relatively inexpensive simulations.

Turbulent water coupling in shock wave lithotripsy

Previous studies have demonstrated that stone comminution decreases with increased pulse repetition frequency as a result of bubble proliferation in the cavitation field of a shock wave lithotripter (Pishchalnikov et al 2011 J. Acoust. Soc. Am. 130 EL87-93). If cavitation nuclei remain in the propagation path of successive lithotripter pulses, especially in the acoustic coupling cushion of the shock wave source, they will consume part of the incident wave energy, leading to reduced tensile pressure in the focal region and thus lower stone comminution efficiency. We introduce a method to remove cavitation nuclei from the coupling cushion between successive shock exposures using a jet of degassed water. As a result, pre-focal bubble nuclei lifetime quantified by B-mode ultrasound imaging was reduced from 7 to 0.3 s by a jet with an exit velocity of 62 cm s(-1). Stone fragmentation (percent mass <2 mm) after 250 shocks delivered at 1 Hz was enhanced from 22 ± 6% to 33 ± 5% (p = 0.007) in water without interposing tissue mimicking materials. Stone fragmentation after 500 shocks delivered at 2 Hz was increased from 18 ± 6% to 28 ± 8% (p = 0.04) with an interposing tissue phantom of 8 cm thick. These results demonstrate the critical influence of cavitation bubbles in the coupling cushion on stone comminution and suggest a potential strategy to improve the efficacy of contemporary shock wave lithotripters.

Size and location of defects at the coupling interface affect lithotripter performance

Study Type--Therapy (case series) Level of Evidence 4. What's known on the subject? and What does the study add? In shock wave lithotripsy air pockets tend to get caught between the therapy head of the lithotripter and the skin of the patient. Defects at the coupling interface hinder the transmission of shock wave energy into the body, reducing the effectiveness of treatment. This in vitro study shows that ineffective coupling not only blocks the transmission of acoustic pulses but also alters the properties of shock waves involved in the mechanisms of stone breakage, with the effect dependent on the size and location of defects at the coupling interface.

OBJECTIVE

• To determine how the size and location of coupling defects caught between the therapy head of a lithotripter and the skin of a surrogate patient (i.e. the acoustic window of a test chamber) affect the features of shock waves responsible for stone breakage.

MATERIALS AND METHODS

• Model defects were placed in the coupling gel between the therapy head of a Dornier Compact-S electromagnetic lithotripter (Dornier MedTech, Kennesaw, GA, USA) and the Mylar (biaxially oriented polyethylene terephthalate) (DuPont Teijin Films, Chester, VA, USA) window of a water-filled coupling test system. • A fibre-optic probe hydrophone was used to measure acoustic pressures and map the lateral dimensions of the focal zone of the lithotripter. • The effect of coupling conditions on stone breakage was assessed using gypsum model stones.

RESULTS

• Stone breakage decreased in proportion to the area of the coupling defect; a centrally located defect blocking only 18% of the transmission area reduced stone breakage by an average of almost 30%. • The effect on stone breakage was greater for defects located on-axis and decreased as the defect was moved laterally; an 18% defect located near the periphery of the coupling window (2.0 cm off-axis) reduced stone breakage by only ~15% compared to when coupling was completely unobstructed. • Defects centred within the coupling window acted to narrow the focal width of the lithotripter; an 8.2% defect reduced the focal width ~30% compared to no obstruction (4.4 mm vs 6.5 mm). • Coupling defects located slightly off centre disrupted the symmetry of the acoustic field; an 18% defect positioned 1.0 cm off-axis shifted the focus of maximum positive pressure ~1.0 mm laterally. • Defects on and off-axis imposed a significant reduction in the energy density of shock waves across the focal zone.

CONCLUSIONS

• In addition to blocking the transmission of shock-wave energy, coupling defects also disrupt the properties of shock waves that play a role in stone breakage, including the focal width of the lithotripter and the symmetry of the acoustic field • The effect is dependent on the size and location of defects, with defects near the centre of the coupling window having the greatest effect. • These data emphasize the importance of eliminating air pockets from the coupling interface, particularly defects located near the centre of the coupling window.