Matthew R. Myers

Influence of incidence angle on sound generation by airfoils interacting with high-frequency gusts

A theoretical model is developed for the sound generated when a convected vortical or entropic gust encounters an airfoil at non-zero angle of attack. The theory is based on a linearization of the Euler equations about the steady subsonic flow past the airfoil. High-frequency gusts, whose wavelengths are short compared to the airfoil chord, but long compared to the displacement of the mean-flow stagnation point from the leading edge, are considered. The analysis utilizes singular-perturbation techniques and involves four asymptotic regions. Local regions, which scale on the gust wavelength, are present at the airfoil leading and trailing edges. Behind the airfoil a ‘transition’ region, which is similar to the transition zone between illuminated and shadow zones in optical problems, is present. In the outer region, far away from the airfoil edges and wake, the solution has a geometric-acoustics form. The primary sound generation is found to be concentrated in the local leading-edge region. The trailing edge plays a secondary role as a scatterer of the sound generated in the leading-edge region. Parametric calculations are presented which illustrate that moderate levels of airfoil steady loading can significantly affect the sound field produced by airfoil-gust interactions.

HIFU procedures at moderate intensities : effect of large blood vessels

A three-dimensional computational model is presented for studying the efficacy of high-intensity focused ultrasound (HIFU) procedures targeted near large blood vessels. The analysis applies to procedures performed at intensities below the threshold for cavitation, boiling and highly nonlinear propagation, but high enough to increase tissue temperature a few degrees per second. The model is based upon the linearized KZK equation and the bioheat equation in tissue. In the blood vessel the momentum and energy equations are satisfied. The model is first validated in a tissue phantom, to verify the absence of bubble formation and nonlinear effects. Temperature rise and lesion-volume calculations are then shown for different beam locations and orientations relative to a large vessel. Both single and multiple ablations are considered. Results show that when the vessel is located within about a beam width (few mm) of the ultrasound beam, significant reduction in lesion volume is observed due to blood flow. However, for gaps larger than a beam width, blood flow has no major effect on the lesion formation. Under the clinically representative conditions considered, the lesion volume is reduced about 40% (relative to the no-flow case) when the beam is parallel to the blood vessel, compared to about 20% for a perpendicular orientation. Procedures involving multiple ablation sites are affected less by blood flow than single ablations. The model also suggests that optimally focused transducers can generate lesions that are significantly larger (>2 times) than the ones produced by highly focused beams.

Multilaboratory Particle Image Velocimetry Analysis of the FDA Benchmark Nozzle Model to Support Validation of Computational Fluid Dynamics Simulations

This study is part of a FDA-sponsored project to evaluate the use and limitations of computational fluid dynamics (CFD) in assessing blood flow parameters related to medical device safety. In an interlaboratory study, fluid velocities and pressures were measured in a nozzle model to provide experimental validation for a companion round-robin CFD study. The simple benchmark nozzle model, which mimicked the flow fields in several medical devices, consisted of a gradual flow constriction, a narrow throat region, and a sudden expansion region where a fluid jet exited the center of the nozzle with recirculation zones near the model walls. Measurements of mean velocity and turbulent flow quantities were made in the benchmark device at three independent laboratories using particle image velocimetry (PIV). Flow measurements were performed over a range of nozzle throat Reynolds numbers (Re(throat)) from 500 to 6500, covering the laminar, transitional, and turbulent flow regimes. A standard operating procedure was developed for performing experiments under controlled temperature and flow conditions and for minimizing systematic errors during PIV image acquisition and processing. For laminar (Re(throat)=500) and turbulent flow conditions (Re(throat)≥3500), the velocities measured by the three laboratories were similar with an interlaboratory uncertainty of ∼10% at most of the locations. However, for the transitional flow case (Re(throat)=2000), the uncertainty in the size and the velocity of the jet at the nozzle exit increased to ∼60% and was very sensitive to the flow conditions. An error analysis showed that by minimizing the variability in the experimental parameters such as flow rate and fluid viscosity to less than 5% and by matching the inlet turbulence level between the laboratories, the uncertainties in the velocities of the transitional flow case could be reduced to ∼15%. The experimental procedure and flow results from this interlaboratory study (available at http://fdacfd.nci.nih.gov) will be useful in validating CFD simulations of the benchmark nozzle model and in performing PIV studies on other medical device models.

In vitro modeling of spinal anesthesia. A digital video image processing technique and its application to catheter characterization.

BackgroundMaldistribution of intrathecal local anesthetic has recently been implicated as a contributor to neurotoxic injury. In vitro modeling can be used to understand the distribution of anesthetic agents within the subarachnoid space. We describe an in vitro modeling technique that uses digital video image processing and its application to catheter injection of local anesthetic. MethodsA clear plastic model of the subarachnoid space, including a simulated spinal cord and cauda equina, was filled with lactated Ringers solution. Phthalocyanine blue dye of known concentration was injected into the model through small-bore (28-G) and large-bore (18-G) catheters. Injections were performed at a variety of controlled rates and sacral catheter positions, and the propagation of dye throughout the model was recorded on videotape, digitized by computer, and converted to a two-dimensional image of dye concentration. A subset of data was compared with results obtained from spectrophotometric analysis. ResultsThere was a strong correlation (r = 0.98) between data obtained with analysis by digital video image processing and those obtained spectrophotometrically. Catheter size, catheter angle, and injection rate significantly influenced the distribution and peak concentration of simulated anesthetic. No major differences in distribution or peak concentration were observed with the two types of 28-G catheters. ConclusionsThe digital video image processing technique can be used to quantify anesthetic distribution rapidly within a model of the subarachnoid space without disturbing the distribution. The current results demonstrate a strong dependence of anesthetic distribution on catheter angle, catheter size, and injection rate. Comparisons between 28-G catheters suggest that the difference in reported incidence of cauda equlna syndrome associated with different 28-G catheters cannot be explained on the basis of differences in anesthetic distribution

Thermal effects generated by high-intensity focused ultrasound beams at normal incidence to a bone surface

Experiments and computations were performed to study factors affecting thermal safety when high-intensity focused ultrasound (HIFU) beams are normally incident (i.e., beam axis normal to the interface) upon a bone/soft-tissue interface. In particular, the temperature rise and thermal dose were determined as a function of separation between the beam focus and the interface. Under conditions representative of clinical HIFU procedures, it was found that the thermal dose at the bone surface can exceed the threshold for necrosis even when the beam focus is more than 4 cm from the bone. Experiments showed that reflection of the HIFU beam from the bone back into the transducer introduced temperature fluctuations of as much as +/-15% and may be an important consideration for safety analyses at sufficiently high acoustic power. The applicability of linear propagation models in predicting thermal dose near the interface was also addressed. Linear models, while underpredicting thermal dose at the focus, provided a conservative (slight overprediction) estimate of thermal dose at the bone surface. Finally, temperature rise due to absorption of shear waves generated by the HIFU beam in the bone was computed. Modeling shear-wave propagation in the thermal analysis showed that the predicted temperature rise off axis was as much as 30% higher when absorption of shear waves is included, indicating that enhanced heating due to shear-wave absorption is potentially important, even for normally incident HIFU beams.

Characterization of high intensity focused ultrasound transducers using acoustic streaming

A new approach for characterizing high intensity focused ultrasound (HIFU) transducers is presented. The technique is based upon the acoustic streaming field generated by absorption of the HIFU beam in a liquid medium. The streaming field is quantified using digital particle image velocimetry, and a numerical algorithm is employed to compute the acoustic intensity field giving rise to the observed streaming field. The method as presented here is applicable to moderate intensity regimes, above the intensities which may be damaging to conventional hydrophones, but below the levels where nonlinear propagation effects are appreciable. Intensity fields and acoustic powers predicted using the streaming method were found to agree within 10% with measurements obtained using hydrophones and radiation force balances. Besides acoustic intensity fields, the streaming technique may be used to determine other important HIFU parameters, such as beam tilt angle or absorption of the propagation medium.

A subarray approach to matched‐field processing

A simplified approach is presented for matched‐field processing of large vertical arrays in typical ocean environments. The large array is divided into subarrays such that a plane‐wave decomposition of the signal field is a good approximation over each subarray. A simple rule, based on the sound‐speed profile, is presented for subarray sizing. Plane‐wave beamforming is then performed to provide multiple beams from each subarray. For long‐range situations, only those beams need to be formed that are steered to receive rays that are not strongly bottom interacting. When elements are spaced at one‐half wavelength, the typical result of considering only these beams is a significant reduction in the size of the problem with far fewer beams required than elements. This first step is very robust and does not require a detailed knowledge of the signal field, or an assumption of signal coherence across the full array aperture. Beam outputs can then be combined to achieve full matched‐field processing. The relative...

Mechanical Bioeffects of Pulsed High Intensity Focused Ultrasound on a Simple Neural Model

PURPOSE To study how pressure pulses affect nerves through mechanisms that are neither thermal nor cavitational, and investigate how the effects are related to cumulative radiation-force impulse (CRFI). Applications include traumatic brain injury and acoustic neuromodulation. METHODS A simple neural model consisting of the giant axon of a live earthworm was exposed to trains of pressure pulses produced by an 825 kHz focused ultrasound transducer. The peak negative pressure of the pulses and duty cycle of the pulse train were controlled so that neither cavitation nor significant temperature rise occurred. The amplitude and conduction velocity of action-potentials triggered in the worm were measured as the magnitude of the pulses and number of pulses in the pulse trains were varied. RESULTS The functionality of the axons decreased when sufficient pulse energy was applied. The level of CRFI at which the observed effects occur is consistent with the lower levels of injury observed in this study relative to blast tubes. The relevant CRFI values are also comparable to CRFI values in other studies showing measureable changes in action-potential amplitudes and velocities. Plotting the measured action-potential amplitudes and conduction velocities from different experiments with widely varying exposure regimens against the single parameter of CRFI yielded values that agreed within 21% in terms of amplitude and 5% in velocity. A predictive model based on the assumption that the temporal rate of decay of action-potential amplitude and velocity is linearly proportional the radiation force experienced by the axon predicted the experimental amplitudes and conduction velocities to within about 20% agreement. CONCLUSIONS The functionality of axons decreased due to noncavitational mechanical effects. The radiation force, possibly by inducing changes in ion-channel permeability, appears to be a possible mechanism for explaining the observed degradation. The CRFI is also a promising parameter for quantifying neural bioeffects during exposure to pressure waves, and for predicting axon functionality.PURPOSE To study how pressure pulses affect nerves through mechanisms that are neither thermal nor cavitational, and investigate how the effects are related to cumulative radiation-force impulse (CRFI). Applications include traumatic brain injury and acoustic neuromodulation. METHODS A simple neural model consisting of the giant axon of a live earthworm was exposed to trains of pressure pulses produced by an 825 kHz focused ultrasound transducer. The peak negative pressure of the pulses and duty cycle of the pulse train were controlled so that neither cavitation nor significant temperature rise occurred. The amplitude and conduction velocity of action-potentials triggered in the worm were measured as the magnitude of the pulses and number of pulses in the pulse trains were varied. RESULTS The functionality of the axons decreased when sufficient pulse energy was applied. The level of CRFI at which the observed effects occur is consistent with the lower levels of injury observed in this study relative to blast tubes. The relevant CRFI values are also comparable to CRFI values in other studies showing measureable changes in action-potential amplitudes and velocities. Plotting the measured action-potential amplitudes and conduction velocities from different experiments with widely varying exposure regimens against the single parameter of CRFI yielded values that agreed within 21% in terms of amplitude and 5% in velocity. A predictive model based on the assumption that the temporal rate of decay of action-potential amplitude and velocity is linearly proportional the radiation force experienced by the axon predicted the experimental amplitudes and conduction velocities to within about 20% agreement. CONCLUSIONS The functionality of axons decreased due to noncavitational mechanical effects. The radiation force, possibly by inducing changes in ion-channel permeability, appears to be a possible mechanism for explaining the observed degradation. The CRFI is also a promising parameter for quantifying neural bioeffects during exposure to pressure waves, and for predicting axon functionality.

Beam localization in HIFU temperature measurements using thermocouples, with application to cooling by large blood vessels

Experimental studies of thermal effects in high-intensity focused ultrasound (HIFU) procedures are often performed with the aid of fine wire thermocouples positioned within tissue phantoms. Thermocouple measurements are subject to several types of error which must be accounted for before reliable inferences can be made on the basis of the measurements. Thermocouple artifact due to viscous heating is one source of error. A second is the uncertainty regarding the position of the beam relative to the target location or the thermocouple junction, due to the error in positioning the beam at the junction. This paper presents a method for determining the location of the beam relative to a fixed pair of thermocouples. The localization technique reduces the uncertainty introduced by positioning errors associated with very narrow HIFU beams. The technique is presented in the context of an investigation into the effect of blood flow through large vessels on the efficacy of HIFU procedures targeted near the vessel. Application of the beam localization method allowed conclusions regarding the effects of blood flow to be drawn from previously inconclusive (because of localization uncertainties) data. Comparison of the position-adjusted transient temperature profiles for flow rates of 0 and 400ml/min showed that blood flow can reduce temperature elevations by more than 10%, when the HIFU focus is within a 2mm distance from the vessel wall. At acoustic power levels of 17.3 and 24.8W there is a 20- to 70-fold decrease in thermal dose due to the convective cooling effect of blood flow, implying a shrinkage in lesion size. The beam-localization technique also revealed the level of thermocouple artifact as a function of sonication time, providing investigators with an indication of the quality of thermocouple data for a given exposure time. The maximum artifact was found to be double the measured temperature rise, during initial few seconds of sonication.

Gaussian representation of high-intensity focused ultrasound beams

A method for fast numerical simulation of high-intensity focused ultrasound beams is derived. The method is based on the frequency-domain representation of the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation, and assumes for each harmonic a Gaussian transverse pressure distribution at all distances from the transducer face. The beamwidths of the harmonics are constrained to vary inversely with the square root of the harmonic number, and as such this method may be viewed as an extension of a quasilinear approximation. The technique is capable of determining pressure or intensity fields of moderately nonlinear high-intensity focused ultrasound beams in water or biological tissue, usually requiring less than a minute of computer time on a modern workstation. Moreover, this method is particularly well suited to high-gain simulations since, unlike traditional finite-difference methods, it is not subject to resolution limitations in the transverse direction. Results are shown to be in reasonable agreement with numerical solutions of the full KZK equation in both tissue and water for moderately nonlinear beams.