Effects of Stone Size on the Comminution Process and Efficiency in Shock Wave Lithotripsy

Improving burst wave lithotripsy effectiveness for small stones and fragments by increasing frequency: theoretical modeling and ex vivo study

INTRODUCTION AND OBJECTIVE

In clinical trial NCT03873259, a 2.6-mm lower pole stone was treated transcutaneously and ex vivo with 390-kHz burst wave lithotripsy (BWL) for 40 minutes and failed to break. The stone was subsequently fragmented with 650-kHz BWL after a 4-minute exposure. This study investigated how to fragment small stones and why varying BWL frequency may more effectively fragment stones to dust.

METHODS

A linear elastic model was used to calculate the stress created inside stones from shock wave lithotripsy (SWL) and different BWL frequencies mimicking the stone's size, shape, lamellar structure, and composition. To test model predictions about the impact of BWL frequency, matched pairs of stones (1-5 mm) were treated at 1) 390 kHz, 2) 830 kHz, and 3) 390 kHz followed by 830 kHz. The mass of fragments greater than 1 and 2 mm was measured over 10 minutes of exposure.

RESULTS

The linear elastic model predicts that the maximum principal stress inside a stone increases to more than 5.5 times the pressure applied by the ultrasound wave as frequency is increased, regardless of composition tested. The threshold frequency for stress amplification is proportionate to the wave speed divided by the stone diameter. Thus, smaller stones may be likely to fragment at higher frequency, but not lower frequency below a limit. Unlike with SWL, this amplification in BWL occurs consistently with spherical and irregularly shaped stones. In water tank experiments, stones smaller than the threshold size broke fastest at high frequency (p=0.0003), whereas larger stones broke equally well to sub-millimeter dust at high, low, or mixed frequency.

CONCLUSIONS

For small stones and fragments, increasing frequency of BWL may produce amplified stress in the stone causing the stone to break. Using the strategies outlined here, stones of all sizes may be turned to dust efficiently with BWL.

Multiphysics Analysis of Ultrasonic Shock Wave Lithotripsy and Side Effects on Surrounding Tissues

Background:

Today, the most common method for kidney stone therapy is extracorporeal shock wave lithotripsy. Current research is a numerical simulation of kidney stone fragmentation via ultrasonic shock waves. Most numerical studies in lithotripsy have been carried out using the elasticity or energy method and neglected the dissipation phenomenon. In the current study, it is solved by not only the linear acoustics equation, but also the Westervelt acoustics equation which nonlinearity and dissipation are involved.

Objective:

This study is to compare two methods for simulation of shock wave lithotripsy, clarifying the effect of shock wave profiles and stones’ material, and investigating side effects on surrounding tissues

Material and Methods:

Computational study is done using COMSOL Multiphysics, commercial software based on the finite element method. Nonlinear governing equations of acoustics, elasticity and bioheat-transfer are coupled and solved.

Results:

A decrease in the rise time of shock wave leads to increase the produced acoustic pressure and enlarge focus region. The shock wave damages kidney tissues in both linear and nonlinear simulation but the damage due to high temperature is very negligible compared to the High Intensity Focused Ultrasound (HIFU).

Conclusion:

Disaffiliation of wave nonlinearity causes a high incompatibility with reality. Stone’s material is an important factor, affecting the fragmentation

Variations of stress field and stone fracture produced at different lateral locations in a shockwave lithotripter field

During clinical procedures, the lithotripter shock wave (LSW) that is incident on the stone and resultant stress field is often asymmetric due to the respiratory motion of the patient. The variations of the LSW-stone interaction and associated fracture pattern were investigated by photoelastic imaging, phantom experiments, and three-dimensional fluid-solid interaction modeling at different lateral locations in a lithotripter field. In contrast to a T-shaped fracture pattern often observed in the posterior region of the disk-shaped stone under symmetric loading, the fracture pattern gradually transitioned to a tilted L-shape under asymmetric loading conditions. Moreover, the model simulations revealed the generation of surface acoustic waves (SAWs), i.e., a leaky Rayleigh wave on the anterior boundary and Scholte wave on the posterior boundary of the stone. The propagation of SAWs on the stone boundary is accompanied by a progressive transition of the LSW reflection pattern from regular to von Neumann and to weak von Neumann reflection near the glancing incidence and, concomitantly, the development and growth of a Mach stem, swirling around the stone boundary. The maximum tensile stress and stress integral were produced by SAWs on the stone boundary under asymmetric loading conditions, which drove the initiation and extension of surface cracks into the bulk of the stone that is confirmed by micro-computed tomography analysis.

Towards an understanding of the chemo-mechanical influences on kidney stone failure via the material point method

This paper explores the use of the meshfree computational mechanics method, the Material Point Method (MPM), to model the composition and damage of typical renal calculi, or kidney stones. Kidney stones are difficult entities to model due to their complex structure and failure behavior. Better understanding of how these stones behave when they are broken apart is a vital piece of knowledge to medical professionals whose aim is to remove these stone by breaking them within a patient’s body. While the properties of individual stones are varied, the common elements and proportions are used to generate synthetic stones that are then placed in a digital experiment to observe their failure patterns. First a more traditional engineering model of a Brazil test is used to create a tensile fracture within the center of these stones to observe the effect of stone consistency on failure behavior. Next a novel application of MPM is applied which relies on an ultrasonic wave being carried by surrounding fluid to model the ultrasonic treatment of stones commonly used by medical practitioners. This numerical modeling of Extracorporeal Shock Wave Lithotripsy (ESWL) reveals how these different stones failure in a more real-world situation and could be used to guide further research in this field for safer and more effective treatments.

An investigation of elastic waves producing stone fracture in burst wave lithotripsy

Burst wave lithotripsy is a method to noninvasively fragment urinary stones by short pulses of focused ultrasound. In this study, physical mechanisms of stone fracture during burst wave lithotripsy were investigated. Photoelasticity imaging was used to visualize elastic wave propagation in model stones and compare results to numerical calculations. Epoxy and glass stone models were made into rectangular, cylindrical, or irregular geometries and exposed in a degassed water bath to focused ultrasound bursts at different frequencies. A high-speed camera was used to record images of the stone during exposure through a circular polariscope backlit by a monochromatic flash source. Imaging showed the development of periodic stresses in the stone body with a pattern dependent on frequency. These patterns were identified as guided wave modes in cylinders and plates, which formed standing waves upon reflection from the distal surfaces of the stone model, producing specific locations of stress concentration in the models. Measured phase velocities compared favorably to numerically calculated modes dependent on frequency and material. Artificial stones exposed to bursts produced cracks at positions anticipated by this mechanism. These results support guided wave generation and reflection as a mechanism of stone fracture in burst wave lithotripsy.

Shock-Induced Damage and Dynamic Fracture in Cylindrical Bodies Submerged in Liquid

Understanding the response of solid materials to shock loading is important for mitigating shock-induced damages and failures, as well as advancing the beneficial use of shock waves for material modifications. In this paper, we consider a representative brittle material, BegoStone, in the form of cylindrical bodies and submerged in water. We present a computational study on the causal relationship between the prescribed shock load and the resulting elastic waves and damage in the solid material. A recently developed three-dimensional computational framework, FIVER, is employed, which couples a finite volume compressible fluid solver with a finite element structural dynamics solver through the construction and solution of local, one-dimensional fluid-solid Riemann problems. The material damage and fracture are modeled and simulated using a continuum damage mechanics model and an element erosion method. The computational model is validated in the context of shock wave lithotripsy and the results are compared with experimental data. We first show that after calibrating the growth rate of microscopic damage and the threshold for macroscopic fracture, the computational framework is capable of capturing the location and shape of the shock-induced fracture observed in a laboratory experiment. Next, we introduce a new phenomenological model of shock waveform, and present a numerical parametric study on the effects of a single shock load, in which the shock waveform, magnitude, and the size of the target material are varied. In particular, we vary the waveform gradually from one that features non-monotonic decay with a tensile phase to one that exhibits monotonic decay without a tensile phase. The result suggests that when the length of the shock pulse is comparable to that of the target material, the former waveform may induce much more significant damage than the latter one, even if the two share the same magnitude, duration, and acoustic energy.

An in vitro evaluation of laser settings and location in the efficiency of the popcorn effect

To examine different locations and laser settings' effects on the efficiency of the "popcorn" method of laser lithotripsy, which consists of placing the laser in a group of small stones and firing continuously to break them into smaller particles. Pre-fragmented BegoStones were created between 2 and 4 mm to mimic typical popcorning conditions. A 0.5 g collection of fragments was placed into 3D-printed models (a spherical calyx and ellipsoid pelvis model) and a 200-µm laser fiber was positioned above the stones. The laser was fired for 2 min with irrigation, with 5 trials at each setting: 0.2 J/50 Hz, 0.5 J/20 Hz, 0.5 J/40 Hz, 1 J/20 Hz, 0.2 J/80 Hz, 0.5 J/80 Hz. After drying, fragmentation efficiency was determined by calculating the mass of stones reduced to sub-2 mm particles. Statistical analysis was performed with ANOVA and Student's t test. The trials within the calyx model were significantly more efficient compared to the pelvis (0.19 vs 0.15 g, p = 0.01). When comparing laser settings, there was a difference between groups by one-way ANOVA [F(5,54) = 8.503, p = 5.47 × 10^-6]. Post hoc tests showed a power setting of 0.5 J/80 Hz was significantly more efficient than low-power settings 0.2 J/50 Hz and 0.5 J/20 Hz (p < 0.05). Additionally, 0.2 J/50 Hz was significantly less efficient than 0.5 J/40 Hz, 1 J/20 Hz, and 0.2 J/80 Hz. Popcorning is most efficient in smaller spaces; we recommend displacement of stones into a calyx before popcorning. No difference was seen between high-power settings, although 0.5 J/40 Hz and 0.5 J/80 Hz performed best, suggesting that moderate energy popcorning methods with at least 0.5 J per pulse are most efficient.

An experimentally-calibrated damage mechanics model for stone fracture in shock wave lithotripsy

A damage model suggested by the Tuler-Butcher concept of dynamic accumulation of microscopic defects is obtained from experimental data on microcrack formation in synthetic kidney stones. Experimental data on appearance of microcracks is extracted from micro-computed tomography images of BegoStone simulants obtained after subjecting the stone to successive pulses produced by an electromagnetic shock-wave lithotripter source. Image processing of the data is used to infer statistical distributions of crack length and width in representative transversal cross-sections of a cylindrical stone. A high-resolution finite volume computational model, capable of accurately modeling internal reflections due to local changes in material properties produced by material damage is used to simulate the accumulation of damage due to successive shocks. Comparison of statistical distributions of microcrack formation in computation and experiment allows calibration of the damage model. The model is subsequently used to compute fracture of a different aspect-ratio cylindrical stone predicting concurrent formation of two main fracture areas as observed experimentally.