Wayne Kreider

Characterization of a multi-element clinical HIFU system using acoustic holography and nonlinear modeling

High-intensity focused ultrasound (HIFU) is a treatment modality that relies on the delivery of acoustic energy to remote tissue sites to induce thermal and/or mechanical tissue ablation. To ensure the safety and efficacy of this medical technology, standard approaches are needed for accurately characterizing the acoustic pressures generated by clinical ultrasound sources under operating conditions. Characterization of HIFU fields is complicated by nonlinear wave propagation and the complexity of phased-array transducers. Previous work has described aspects of an approach that combines measurements and modeling, and here we demonstrate this approach for a clinical phased-array transducer. First, low amplitude hydrophone measurements were performed in water over a scan plane between the array and the focus. Second, these measurements were used to holographically reconstruct the surface vibrations of the transducer and to set a boundary condition for a 3-D acoustic propagation model. Finally, nonlinear simulations of the acoustic field were carried out over a range of source power levels. Simulation results were compared with pressure waveforms measured directly by hydrophone at both low and high power levels, demonstrating that details of the acoustic field, including shock formation, are quantitatively predicted.

Fragmentation of Urinary Calculi In Vitro by Burst Wave Lithotripsy

PURPOSE We developed a new method of lithotripsy that uses short, broadly focused bursts of ultrasound rather than shock waves to fragment stones. We investigated the characteristics of stone comminution by burst wave lithotripsy in vitro. MATERIALS AND METHODS Artificial and natural stones (mean ± SD size 8.2 ± 3.0 mm, range 5 to 15) were treated with ultrasound bursts using a focused transducer in a water bath. Stones were exposed to bursts with focal pressure amplitude of 6.5 MPa or less at a 200 Hz burst repetition rate until completely fragmented. Ultrasound frequencies of 170, 285 and 800 kHz were applied using 3 transducers, respectively. Time to fragmentation for each stone type was recorded and fragment size distribution was measured by sieving. RESULTS Stones exposed to ultrasound bursts were fragmented at focal pressure amplitudes of 2.8 MPa or greater at 170 kHz. Fractures appeared along the stone surface, resulting in fragments that separated at the surface nearest to the transducer until the stone was disintegrated. All natural and artificial stones were fragmented at the highest focal pressure of 6.5 MPa with a mean treatment duration of 36 seconds for uric acid stones to 14.7 minutes for cystine stones. At a frequency of 170 kHz the largest artificial stone fragments were less than 4 mm. Exposure at 285 and 800 kHz produced only fragments less than 2 mm and less than 1 mm, respectively. CONCLUSIONS Stone comminution with burst wave lithotripsy is feasible as a potential noninvasive treatment method for nephrolithiasis. Adjusting the fundamental ultrasound frequency allows for stone fragment size to be controlled.

Passive Cavitation Detection during Pulsed HIFU Exposures of Ex Vivo Tissues and In Vivo Mouse Pancreatic Tumors

Pulsed high-intensity focused ultrasound (pHIFU) has been shown to enhance vascular permeability, disrupt tumor barriers and enhance drug penetration into tumor tissue through acoustic cavitation. Monitoring of cavitation activity during pHIFU treatments and knowing the ultrasound pressure levels sufficient to reliably induce cavitation in a given tissue are therefore very important. Here, three metrics of cavitation activity induced by pHIFU and evaluated by confocal passive cavitation detection were introduced: cavitation probability, cavitation persistence and the level of the broadband acoustic emissions. These metrics were used to characterize cavitation activity in several ex vivo tissue types (bovine tongue and liver and porcine adipose tissue and kidney) and gel phantoms (polyacrylamide and agarose) at varying peak-rare factional focal pressures (1-12 MPa) during the following pHIFU protocol: frequency 1.1 MHz, pulse duration 1 ms and pulse repetition frequency 1 Hz. To evaluate the relevance of the measurements in ex vivo tissue, cavitation metrics were also investigated and compared in the ex vivo and in vivo murine pancreatic tumors that develop spontaneously in transgenic KrasLSL.G12 D/+; p53 R172 H/+; PdxCretg/+ (KPC) mice and closely re-capitulate human disease in their morphology. The cavitation threshold, defined at 50% cavitation probability, was found to vary broadly among the investigated tissues (within 2.5-10 MPa), depending mostly on the water-lipid ratio that characterizes the tissue composition. Cavitation persistence and the intensity of broadband emissions depended both on tissue structure and lipid concentration. Both the cavitation threshold and broadband noise emission level were similar between ex vivo and in vivo pancreatic tumor tissue. The largest difference between in vivo and ex vivo settings was found in the pattern of cavitation occurrence throughout pHIFU exposure: it was sporadic in vivo, but it decreased rapidly and stopped over the first few pulses ex vivo. Cavitation activity depended on the interplay between the destruction and circulation of cavitation nuclei, which are not only used up by HIFU treatment but also replenished or carried away by circulation in vivo. These findings are important for treatment planning and optimization in pHIFU-induced drug delivery, in particular for pancreatic tumors.

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.

Acoustic holography as a metrological tool for characterizing medical ultrasound sources and fields

Acoustic holography is a powerful technique for characterizing ultrasound sources and the fields they radiate, with the ability to quantify source vibrations and reduce the number of required measurements. These capabilities are increasingly appealing for meeting measurement standards in medical ultrasound; however, associated uncertainties have not been investigated systematically. Here errors associated with holographic representations of a linear, continuous-wave ultrasound field are studied. To facilitate the analysis, error metrics are defined explicitly, and a detailed description of a holography formulation based on the Rayleigh integral is provided. Errors are evaluated both for simulations of a typical therapeutic ultrasound source and for physical experiments with three different ultrasound sources. Simulated experiments explore sampling errors introduced by the use of a finite number of measurements, geometric uncertainties in the actual positions of acquired measurements, and uncertainties in the properties of the propagation medium. Results demonstrate the theoretical feasibility of keeping errors less than about 1%. Typical errors in physical experiments were somewhat larger, on the order of a few percent; comparison with simulations provides specific guidelines for improving the experimental implementation to reduce these errors. Overall, results suggest that holography can be implemented successfully as a metrological tool with small, quantifiable errors.

Design of HIFU Transducers for Generating Specified Nonlinear Ultrasound Fields

Various clinical applications of high-intensity focused ultrasound have different requirements for the pressure levels and degree of nonlinear waveform distortion at the focus. The goal of this paper is to determine transducer design parameters that produce either a specified shock amplitude in the focal waveform or specified peak pressures while still maintaining quasi-linear conditions at the focus. Multiparametric nonlinear modeling based on the Khokhlov–Zabolotskaya–Kuznetsov (KZK) equation with an equivalent source boundary condition was employed. Peak pressures, shock amplitudes at the focus, and corresponding source outputs were determined for different transducer geometries and levels of nonlinear distortion. The results are presented in terms of the parameters of an equivalent single-element spherically shaped transducer. The accuracy of the method and its applicability to cases of strongly focused transducers were validated by comparing the KZK modeling data with measurements and nonlinear full diffraction simulations for a single-element source and arrays with 7 and 256 elements. The results provide look-up data for evaluating nonlinear distortions at the focus of existing therapeutic systems as well as for guiding the design of new transducers that generate specified nonlinear fields.

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.

The dynamics of histotripsy bubbles

Histotripsy describes treatments in which high‐amplitude acoustic pulses are used to excite bubbles and erode tissue. Though tissue erosion can be directly attributed to bubble activity, the genesis and dynamics of bubbles remain unclear. Histotripsy lesions that show no signs of thermal coagulative damage have been generated with two different acoustic protocols: relatively long acoustic pulses that produce local boiling within milliseconds and relatively short pulses that are higher in amplitude but likely do not produce boiling. While these two approaches are often distinguished as ‘boiling’ versus ‘cavitation’, such labels can obscure similarities. In both cases, a bubble undergoes large changes in radius and vapor is transported into and out of the bubble as it oscillates. Moreover, observations from both approaches suggest that bubbles grow to a size at which they cease to collapse violently. In order to better understand the dynamics of histotripsy bubbles, a single‐bubble model has been developed ...

Cavitation-induced damage of soft materials by focused ultrasound bursts: A fracture-based bubble dynamics model.

A generalized Rayleigh-Plesset-type bubble dynamics model with a damage mechanism is developed for cavitation and damage of soft materials by focused ultrasound bursts. This study is linked to recent experimental observations in tissue-mimicking polyacrylamide and agar gel phantoms subjected to bursts of a kind being considered specifically for lithotripsy. These show bubble activation at multiple sites during the initial pulses. More cavities appear continuously through the course of the observations, similar to what is deduced in pig kidney tissues in shock-wave lithotripsy. Two different material models are used to represent the distinct properties of the two gel materials. The polyacrylamide gel is represented with a neo-Hookean elastic model and damaged based upon a maximum-strain criterion; the agar gel is represented with a strain-hardening Fung model and damaged according to the strain-energy-based Griffiths fracture criterion. Estimates based upon independently determined elasticity and viscosity of the two gel materials suggest that bubble confinement should be sufficient to prevent damage in the gels, and presumably injury in some tissues. Damage accumulation is therefore proposed to occur via a material fatigue, which is shown to be consistent with observed delays in widespread cavitation activity.

Beamwidth measurement of individual lithotripter shock waves

New lithotripters with narrower foci and higher peak pressures than the original Dornier HM3 electrohydraulic lithotripter have proven to be less effective and less safe. Hence, accurate measurements of the focal characteristics of lithotripter shock waves are important. The current technique for measuring beamwidth requires a collection of single-point measurements over multiple shock waves, thereby introducing error as a result of any shock-to-shock variability. This work reports on the construction of a hydrophone array sensor and on array measurements of individual lithotripter shock waves. Beamwidths for an electrohydraulic lithotripter with a broad-focus HM3-style reflector and a narrow-focus modified reflector were measured using both new and worn electrodes as well as two different electrical charging potentials. The array measured the waveform, beamwidth, and focal location of individual shock waves. The HM3-style reflector produced repeatable focal waveforms and beam profiles at an 18 kV charging potential with new and worn electrodes. Corresponding measurements suggest a narrower beamwidth than reported previously from averaged point measurements acquired under the same conditions. In addition, a lack of consistency in the measured beam profiles at 23 kV underscores the value of measuring individual shock waves.