Georgy N. Sankin
Zheng Fang1†
Pei Zhong- 1Thomas Lord Department of Mechanical Engineering and Materials Science, Duke University, Durham, NC, United States
- 2Department of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC, United States
- 3Graduate Program in Acoustics, The Pennsylvania State University, University Park, PA, United States
Background and objective: All clinical shock wave lithotripters produce an axisymmetric acoustic field without accounting for the anatomic features of the kidney or respiratory motion of the patient. This work presents a steerable and adjustable focusing electrohydraulic (SAFE) shock wave generator design with variable beam size and shape.
Materials and methods: 90 electrohydraulic transducers are mounted concentrically on a spherical basin with adjustable connection to individual transducers. Each transducer consists of 45 3D-printed titanium microelectrodes embedded in epoxy with a tip diameter of 0.3 mm. All the transducers are arranged in 5 concentric rings and sub-divided into 6 sectors.
Results: By changing the connections of individual transducers, the focused pressure field produced by the transducer array can be either axisymmetric with a -6 dB focal width of 14.8 mm in diameter, or non-axisymmetric with a long axis of 22.7 mm and a short axis of 15.1 mm. The elongated beam produces a peak positive pressure of 33.7 ± 4.1 MPa and comminution efficiency of 42.2 ± 3.5%, compared to 36.2 ± 0.7 MPa and 28.6 ± 6.1% for axisymmetric beam after 150 pulses at 20 kV.
Conclusions: We have demonstrated that the SAFE shock wave generator can produce an elongated non-axisymmetric pressure field with higher stone comminution efficiency. The SAFE shock wave generator may provide a flexible and versatile design to achieve accurate, stable, and safe lithotripsy for kidney stone treatment.
Introduction
Kidney stone disease (or urolithiasis) affects 5-15% of the population worldwide (1). In the United States, shock wave lithotripsy (SWL) is still used as the first line therapy for about 50% of the stone patients in spite of the technical advances and growing popularity of ureteroscopy in the past decade (2). During clinical SWL, however, elevated dose delivered at fast rate is often preferred, in association with the high pressure and pulse energy generated by the contemporary lithotripters, which increases the risk of renal injury, i.e., hemorrhage and perinephric hematomas (3–6). It is therefore highly desirable to reduce the total number of high energy shock waves delivered to the patient to minimize the risk of tissue injury while maintaining successful stone comminution in SWL.
Two primary factors have been identified that contribute to the elevated shock wave exposure during SWL: 1) patient’s respiratory motion and 2) residual fragment dispersion (7, 8). Using B–mode ultrasound, Sorensen et al. have shown that kidney stones may move by 10-40 mm and ureteral stones by 7-10 mm in patients under sedation during SWL, leading to about 40% of the lithotripter-generated shock waves (LSWs) missing the target stones (7). In addition, following the initial breakup, residual fragments may spread out inside the renal collecting system to an area exceeding the fragmentation zone of the lithotripter, rendering the ensuing shock waves less effective (8).
The respiration trajectory of the kidney is roughly aligned in parallel to the spine along the head-to-toe direction of the patient. In contrast, the spreading of residual fragments will depend on the anatomic features of the kidney (i.e., renal calyces and pelvis), or the orientation of the ureter where the stone resides. Matching the pressure or energy distribution with these anatomic features or respiratory motion patterns in stone patients during SWL are difficult using contemporary clinical lithotripters because they can only generate axisymmetric pressure fields along the LSW propagation direction. On the other hand, it is conceivable that an optimized distribution of the shock wave energy to maximally match with the stone or fragments trajectory during SWL will improve treatment efficiency while minimizing tissue injury.
We have previously used a foam mask placed above the acoustic lens of an electromagnetic shock wave source to transform its pressure distribution in the focal plane from an axisymmetric to an elongated non-axisymmetric field (9, 10). We have shown in vitro that by matching the pattern of respiratory stone-movement in patient, 30% improvement in targeting accuracy can be achieved (9). The mask design, however, cannot be used easily for flexible re-configuration of the acoustic field to track a moving stone target.
In this study, we describe the development of a prototype steerable and adjustable focusing electrohydraulic (SAFE) shock wave generator for flexible lithotripter field control utilizing a re-configurable multi-spark transducer array. We present the design principle, preliminary tests of acoustic field characterization and in vitro stone comminution, demonstrating several unique advantages of this new shock wave source that may provide flexible beam forming and stone targeting capability during clinical SWL procedures.
Materials and methods
Design concept
In the SAFE shock wave generator, 90 individual electrohydraulic transducers are mounted on a spherical carrier made of dielectric material. Each transducer consists of 45 3D-printed titanium electrodes embedded in epoxy resin and each polished to yield a tip diameter of 0.3 mm following previous studies (11, 12). The transducers are arranged in 5 concentric rings divided into 6 sectors (Figure 1).
Figure 1 (A) CAD model of the 45-pin transducer; (B) enlarged photo of the 45-pin transducer on the SAFE; (C) illustration of wave propagated from the SAFE and focus at kidney stones, with two side sections in red shaded areas disconnected, leading to a focal zone elongation in the side-section direction.
In linear acoustics, the beam width is inversely proportional to the aperture diameter of the acoustic source (13). Therefore, for a given focal length, by reducing the active section dimension of the shock wave generator, for example along a 45-degree angle with respect to the y-axis in Figure 2, we can effectively increase the corresponding beam width in a particular direction. This design strategy is used to transform an axisymmetric pressure field into a non-axisymmetric, elongated focal zone (Figure 2).
Figure 2 Schematic view of the transducers distribution on the spherical carrier and resulting pressure contour in the focal plane in (A) and (C) Case A, and (B) and (D) Case B.
The numerical simulations throughout this pilot study were carried out using an open-source toolbox k-Wave (14) in time-domain under linear wave propagation assumption. While the peak pressure may be underestimated, the linear wave propagation model is anticipated to qualitatively capturing the focal pressure distribution for the purpose of comparison between the axisymmetric vs. elongated pressure fields. In the simulation, the voxel size was 1.875 mm, whereas the time step size was 0.156 microseconds. The excitation waveform was taken from Eq. 8 in (15). More realistic waveform generated by micro-spark at individual pin tip will be used in future simulations.
Experiments were conducted using two different beam configurations, designed for employing the same total number of activated pins (Np) yet with different transducer distributions in the SAFE shock wave generator. To produce a non-axisymmetric or elongated beam (Case A), 12 transducers on the two opposite sectors (12x2 = 24) were disconnected (Figure 1C). In contrast, to produce an axisymmetric beam (Case B), all 4 transducers in the third ring for each of the six sectors (6x4 = 24) were disconnected (see Figure 2). Hence, the total number of transducers (Nt) where Nt = 90 - 24 = 66 were activated in both case A and case B. If we assume that the total input electric energy [Etotal = CU2/2 = 0.5×3μF×(15 kV)2 = 0.34 kJ] was evenly distributed to Np activated pins, the electric energy delivered to each activated pin (Ep) can then be calculated by:
where Np = Nt × 45 = 2970, C = 3 µF is the capacitance, and U is the charging voltage.
Acoustic field characterization
The SAFE shock wave generator, installed at the bottom of a cylindrical acrylic water tank (Ø254×H305 mm), was triggered by a digital delay pulse generator (BNC Model 555, Berkley Nucleonics). A metal mesh with grid size 12.5×12.5 mm was mounted above the surface of the titanium electrodes and grounded. The tank was filled with electrolyte (1% sodium chloride in water) to ensure synchronization of the spark discharges from all transducer tips. A fiber optic probe hydrophone (FOPH 500, RP Acoustics, Leutenbach, Germany) mounted on a computer-controlled 3D translational stage (VXM-2 step motors with BiSlide-M02 lead screw, Velmex, Bloomfield, NY) was used for scanning and pressure measurements in the shock wave focal plane.
The pressure measurements were performed at U = 15 kV for beam configuration comparison. In addition, pressure was measured at the beam focus under U = 20 kV, which was used for stone comminution experiments. The pressure measurements were repeated at least three times at each point.
Stone fragmentation assessment
Stone fragmentation experiments were performed using cylindrical soft BegoStone (5:2 powder to water mixing ratio) of Ø6×H6 mm in size and 0.33 g in weight (16). Stone phantoms were soaked in water for more than 30 minutes before the fragmentation test. As shown in Figure 3, a polyurethane rubber stone holder (Ø48mm outer diameter, 30 mm height) in elliptical shape (long axis: 24 mm, and short axis: 12 mm) was used to approximate the anatomic geometry of the renal pelvis while allowing residual stone fragments to be dispersed laterally during SWL.
Figure 3 (A) Photo of the stone holder with the illustration of stone phantom location; (B) diagram of the experimental setup for stone fragmentation test.
A water circulation system was constructed to facilitate removal of bubble remnants accumulated underneath the stone holder during SWL (17). The pre-soaked stone sample was placed inside the polyurethane holder and aligned with the focus of the SAFE shock wave generator. After treatment with 150 pulses produced at 20 kV, stone fragments were collected and air dried overnight. Afterwards, fragments were filtered through 2 mm and 2.8 mm grid sieves (W.S. Tyler, Mentor, OH), and weighted to calculate the stone comminution efficiency based on the percentage of the residual fragments over the original stone weight. Six stones were treated for each case. Data were post-processed in Excel (Microsoft, Redmond, WA) and presented in bar chart with mean, standard deviation, and p-value.
Results and discussion
Figure 4A shows the peak pressure (p+) distribution along two orthogonal directions in the focal plane (z = 0 mm). The -6 dB focal width, estimated by the full width at half maximum using a Gaussian curve fitting for p+, was found to be 22.7 mm (along the side-section direction) by 15.1 mm (in the orthogonal direction) for the elongated beam. In comparison, the -6 dB focal width of the axisymmetric focal zone was about 14.8 mm in both directions. Moreover, the pressure distributions along the z-axis were found to be comparable between the two configurations (Figure 4B). The slope of pressure change post-focally (i.e., z ≥ 0) is steeper than its counterpart pre-focally (i.e., z ≤ 0). These results are summarized in Table 1. In comparison, based on linear wave model simulation, the -6 dB focal width is 18.5 mm along the side-section direction (Case A, non-axisymmetric, Figure 2) by 11.5 mm. In contrast, the -6 dB focal width of the axisymmetric focal zone (Case B, axisymmetric, Figure 2) is 13.0 mm along both directions. In general, the trends in the focal width change of the shock wave generator between experimental measurements and model simulations are similar, with an average discrepancy about 20%. These discrepancies are likely to be reduced when the nonlinear wave propagation is included in the future modeling work.
Figure 4 Peak positive pressure (p+) and peak negative pressure (p-) distribution comparison on (A) focal plane (x- and y-axis), and (B) along wave propagation path (z-axis) between Case A and Case B.
Table 1 Characteristics comparison of the SAFE in different cases.
As the input voltage increases from 15 kV to 20 kV (with the corresponding total input electric energy varying from 338 J to 600 J), the value of p+ increases from 21.1 ± 1.1 MPa to 33.7 ± 4.1 MPa for the elongated beam. Similarly, p+ increases from 21.8 ± 0.8 MPa to 36.2 ± 0.7 MPa for the axisymmetric beam (Figure 5). At 20 kV, the rise time of the shock wavefront is in the range of 0.2 μs to 0.3 μs (Table 2). These acoustic field parameters of the SAFE shock wave generator are comparable to the corresponding values in an HM3 lithotripter, except the longer rise time of the shock front (see Table 2). In comparison to the HM3, the SAFE shock wave generator has the unique advantage of transforming the axisymmetric acoustic field in the lithotripter focal plane to an elongated (oval shape) pressure field that can better match with the anatomical features of the kidney and/or the trajectory of respiratory motion of the patients during SWL. More importantly, the SAFE shock wave generator has the potential of flexible control of the lithotripter focal beam size, shape, and orientations.
Table 2 Comparison of the acoustic fields and stone comminution (SC) efficiency produced by the SAFE and Dornier HM3.
Figure 5 Pressure waveforms measured at focal point when input voltage U = 15 kV and 20 kV in (A) and (B) Case A, and (C) and (D) Case B.
A total of 12 stones were treated. However, one outlier in each group was detected and removed, resulting in a final sample size of n = 5 (Figure 6A). For the elongated beam, stone fragmentation rate for fragments less than 2.0 mm was 42.2 ± 3.5% compared to 28.6 ± 6.1% for the axisymmetric beam. For stone fragments smaller than 2.8 mm, stone comminution for the elongated beam (Case A: 80% ± 9%) is greater than that of the axisymmetric beam (Case B: 47.7% ± 5.1%). Two-tail t-test of the data from the two groups show p-values of 0.0043 (<2 mm) and 0.0003 for (<2.8 mm), respectively, indicating statistically significant difference in stone fragmentation produced by the two configurations of the SAFE shock wave generator. Graphically, Figures 6B, C display the stone fragments after SWL treatment. Stones treated by the elongated beam have more fragments and smaller sizes compared to those produced by the axisymmetric beam. These results suggest that higher stone comminution efficiency may be produced by adjusting the beam size and shape to better match with the target stone/fragments trajectory during clinical SWL procedures.
Figure 6 (A) Stone comminution (n = 5) of fragments smaller than 2 mm (left) and 2.8 mm (right) after 150 pulses of shock wave treatment by the SAFE in case A (blue) and case B (red). Indicated p-values smaller than 0.05 for each group of results. Stone fragments after 150 pulses of shock wave treatment by the SAFE in both (B) Case A and (C) Case B.
Conclusion
In this study, we present the prototype design and evaluation of a steerable and adjustable focusing electrohydraulic (SAFE) shock wave generator. Such a design provides beam-forming flexibility in SWL, which allows us to better match the acoustic field of the lithotripter with anatomic features and spreading of residual fragments in stone patients. Improved stone comminution efficiency has been demonstrated in an elliptical stone holder. Future work is warranted in electronic control of beam forming and steering, evaluation of safety and treatment efficiency in animal models. Pressure measurements and model simulation along the direction of the shock wave propagation will also be conducted.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material. Further inquiries can be directed to the corresponding author.
Author contributions
GS designed and constructed the multi-spark shock wave generator, and together with ZF carried out the pressure measurement and stone comminution tests. YJ performed numerical model simulation of the axisymmetric and non-axisymmetric pressure field with technical assistance of JG. PZ supervised the development of the study plan, data analysis and preparation of the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
We want to thank Patrick McGuire for 3D printing titanium parts, Eric Stach for help with SolidWorks design, and Jun Ji for extracting the coordinates of the pins. This work was partially supported by National Institute of Health (NIH) through Grants 5R01DK052985-25 and 1R01EB025205.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
SWL, shockwave lithotripsy; LSW, lithotripter shock wave; FW, focal width; SC, stone comminution; CAD, computer-aided design.
References
1. Moe OW. Kidney stones: Pathophysiology and medical management. Lancet (2006) 367(9507):333–44. doi: 10.1016/S0140-6736(06)68071-9
2. Pearle MS, Calhoun EA, Curhan GC. Urologic diseases in America project: Urolithiasis. J Urol (2005) 173(3):848–57. doi: 10.1097/01.ju.0000152082.14384.d7
3. Delius M, Enders G, Xuan ZR, Liebich HG, Brendel W. Biological effects of shock-waves - kidney damage by shock-waves in dogs - dose dependence. Ultrasound Med Biol (1988) 14(2):117–22. doi: 10.1016/0301-5629(88)90178-0
4. Delius M, Jordan M, Liebich HG, Brendel W. Biological effects of shock waves: Effect of shock waves on the liver and gallbladder wall of dogs–administration rate dependence. Ultrasound Med Biol (1990) 16(5):459–66. doi: 10.1016/0301-5629(90)90168-C
5. Willis LR, Evan AP, Connors BA, Blomgren P, Fineberg NS, Lingeman JE. Relationship between kidney size, renal injury, and renal impairment induced by shock wave lithotripsy. J Am Soc Nephrol (1999) 10(8):1753–62. doi: 10.1681/ASN.V1081753
6. Willis LR, Evan AP, Connors BA, Shao Y, Blomgren PM, Pratt JH, et al. Shockwave lithotripsy: Dose-related effects on renal structure, hemodynamics, and tubular function. J Endourol (2005) 19(1):90–101. doi: 10.1089/end.2005.19.90
7. Sorensen MD, Bailey MR, Shah AR, Hsi RS, Paun M, Harper JD. Quantitative assessment of shockwave lithotripsy accuracy and the effect of respiratory motion. J Endourol (2012) 26(8):1070–4. doi: 10.1089/end.2012.0042
8. Qin J, Simmons WN, Sankin G, Zhong P. Effect of lithotripter focal width on stone comminution in shock wave lithotripsy. J Acoustical Soc America (2010) 127(4):2635–45. doi: 10.1121/1.3308409
9. Yang TW, Lautz J, Sankin G, Zhong P. Respiration gating beam steering to precisely target a movable stone in shock wave lithotripsy, in: The 5th Annual IEEE International Conference on Cyber Technology in Automation, Control and Intelligent Systems. (Shenyang, China: IEEE) (2015). pp. 1174–9.
10. Lautz J. Non-axisymmetric and steerable acoustic field for enhanced stone comminution in shock wave lithotripsy. PhD Dissertation. Durham NC: Duke University (2014).
11. Teslenko VS, Zhukov AI, Mitrofanov VV, Drozhzhin AP. Generation and focusing of shock-acoustic waves in a liquid by a multicenter electric discharge. Tech Phys (1999) 44(4):476–7. doi: 10.1134/1.1259328
12. Sankin GN, Drozhzhin AP, Lomanovich KA, Teslenko VS. A multisite electric-discharge diaphragm generator of shock waves in a liquid. Instrum Exp Techniques (2004) 47(4):525–8. doi: 10.1023/B:INET.0000038402.96944.eb
13. Pierce AD. Acoustics: An introduction to its physical principles and applications. Phys Today (1981) 34(12):56. doi: 10.1063/1.2914388
14. Treeby BE, Cox BT. K-wave: MATLAB toolbox for the simulation and reconstruction of photoacoustic wave fields. J Biomed Optics (2010) 15(2):1–12, 12. doi: 10.1117/1.3360308
15. Zhou Y, Zhong P. The effect of reflector geometry on the acoustic field and bubble dynamics produced by an electrohydraulic shock wave lithotripter. J Acoust Soc Am (2006) 119(6):3625–36. doi: 10.1121/1.2195074
16. Esch E, Simmons WN, Sankin G, Cocks HF, Preminger GM, Zhong P. A simple method for fabricating artificial kidney stones of different physical properties. Urol Res (2010) 38(4):315–9. doi: 10.1007/s00240-010-0298-x
17. Lautz J, Sankin G, Zhong P. Turbulent water coupling in shock wave lithotripsy. Phys Med Biol (2013) 58(3):735–48. doi: 10.1088/0031-9155/58/3/735
18. Zhong P. Shock wave lithotripsy. In: Delale CF, editor. Bubble dynamics &Shock waves. Berin Heidelberg: Springer-Verlag (2013). p. 291–338.
19. Coleman AJ, Saunders JE. A review of the physical properties and biological effects of the high amplitude acoustic fields used in extracorporeal lithotripsy. Ultrasonics (1993) 31(2):75–89. doi: 10.1016/0041-624X(93)90037-Z
20. Eisenmenger W, Du XX, Tang C, Zhao S, Wang Y, Rong F, et al. The first clinical results of "wide-focus and low-pressure" ESWL. Ultrasound Med Biol (2002) 28(6):769–74. doi: 10.1016/S0301-5629(02)00507-0
21. Cleveland RO, McAteer JA. The physics of shock wave lithotripsy. In: Smith’s textbook on endourology. Berlin Heidelberg New York: Blackwell Publishing Ltd (2007).
Keywords: shock wave lithotripsy, electrohydraulic lithotriptor, pressure field characteristics, stone comminution, beam steering
Citation: Sankin GN, Fang Z, Gu J, Jing Y and Zhong P (2023) A multi-spark electrohydraulic shock wave generator with adjustable pressure field distribution and beam steering capability. Front. Urol. 3:1057723. doi: 10.3389/fruro.2023.1057723
Received: 29 September 2022; Accepted: 27 February 2023;
Published: 14 March 2023.
Edited by:
Dirk Lange, University of British Columbia, CanadaReviewed by:
Alberto Pérez-Lanzac, Alfonso X el Sabio University, SpainYufeng Zhou, Chongqing Medical University, China
Dong Zhang, Nanjing University, China
Mahesh R. Desai, Muljibhai Patel Urological Hospital, India
Copyright © 2023 Sankin, Fang, Gu, Jing and Zhong. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Pei Zhong, cHpob25nQGR1a2UuZWR1
†These authors have contributed equally to this work