Modeling the propagation of nonlinear three-dimensional acoustic beams in inhomogeneous media

Characterization of weakly nonlinear effects in relationship to transducer parameters in focused ultrasound therapy

BACKGROUND

Focused ultrasound therapy has been widely used for the treatment of various diseases, employing different types of transducers. The focused ultrasound pressure fields inevitably exhibit nonlinear effects, which can influence the ablation region. However, the nonlinear effects exhibit noticeable variations across different applications. The characterization of the nonlinear pressure fields of ultrasound is important for the effective implementation of focused ultrasound therapy.

PURPOSE

The traditional angular spectrum method (ASM) was extended to accurately and efficiently simulate the propagation of weakly nonlinear ultrasound in heterogeneous mediums of clinical model. The nonlinear effects were further analyzed in relationship to the transducer parameters that are different in various applications.

METHODS

The pressure fields were simulated using the extended ASM, incorporating calculations for phase aberration in the frequency domain and magnitude compensation in the spatial domain to account for heterogeneous acoustic impedance mismatch. Validation was performed by comparison to k-Wave simulation results using two simplified clinical models, an abdominal soft tissue and a transcranial skull model. The nonlinear effects were then analyzed in relation to the transducer parameters of f-number and effective source area based on the same acoustic output power. The analysis of nonlinear effects was conducted under both homogeneous medium and the clinical models.

RESULTS

The simulation results demonstrated a maximum error of 3.93% in the calculated harmonic pressure of the abdominal model, and a maximum error of 4.89% within the transcranial model when comparing the extended ASM simulation results to those obtained from k-Wave simulations. The characterization of the nonlinear effects reveals a strong correlation with the transducer parameters. Specifically, the results indicate that the nonlinear effects intensify with an increase in the effective source area and f-number, under the same acoustic output power of the transducer. However, the clinical model also showed an influence on the nonlinear effects in relation to the f-number.

CONCLUSION

The extended ASM was demonstrated as an accurate and efficient simulation tool, and the simulation results provide a reference for evaluating the intensity of nonlinear effects in various transducer designs.

Weakly nonlinear focused ultrasound in viscoelastic media containing multiple bubbles

To facilitate practical medical applications such as cancer treatment utilizing focused ultrasound and bubbles, a mathematical model that can describe the soft viscoelasticity of human body, the nonlinear propagation of focused ultrasound, and the nonlinear oscillations of multiple bubbles is theoretically derived and numerically solved. The Zener viscoelastic model and Keller-Miksis bubble equation, which have been used for analyses of single or few bubbles in viscoelastic liquid, are used to model the liquid containing multiple bubbles. From the theoretical analysis based on the perturbation expansion with the multiple-scales method, the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation, which has been used as a mathematical model of weakly nonlinear propagation in single phase liquid, is extended to viscoelastic liquid containing multiple bubbles. The results show that liquid elasticity decreases the magnitudes of the nonlinearity, dissipation, and dispersion of ultrasound and increases the phase velocity of the ultrasound and linear natural frequency of the bubble oscillation. From the numerical calculation of resultant KZK equation, the spatial distribution of the liquid pressure fluctuation for the focused ultrasound is obtained for cases in which the liquid is water or liver tissue. In addition, frequency analysis is carried out using the fast Fourier transform, and the generation of higher harmonic components is compared for water and liver tissue. The elasticity supresses the generation of higher harmonic components and promotes the remnant of the fundamental frequency components. This indicates that the elasticity of liquid suppresses shock wave formation in practical applications.

Weakly nonlinear propagation of focused ultrasound in bubbly liquids with a thermal effect: Derivation of two cases of Khokolov-Zabolotskaya-Kuznetsoz equations

A physico-mathematical model composed of a single equation that consistently describes nonlinear focused ultrasound, bubble oscillations, and temperature fluctuations is theoretically proposed for microbubble-enhanced medical applications. The Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation that has been widely used as a simplified model for nonlinear propagation of focused ultrasound in pure liquid is extended to that in liquid containing many spherical microbubbles, by applying the method of multiple scales to the volumetric averaged basic equations for bubbly liquids. As a result, for two-dimensional and three-dimensional cases, KZK equations composed of the linear combination of nonlinear, dissipation, dispersion, and focusing terms are derived. Especially, the dissipation term depends on three factors, i.e., interfacial liquid viscosity, liquid compressibility, and thermal conductivity of gas inside bubbles; the thermal conduction is evaluated by using four types of temperature gradient models. Finally, we numerically solve the derived KZK equation and show a moderate temperature rise appropriate to medical applications.

"HIFU Beam:" A Simulator for Predicting Axially Symmetric Nonlinear Acoustic Fields Generated by Focused Transducers in a Layered Medium

"HIFU beam" is a freely available software tool that comprises a MATLAB toolbox combined with a user-friendly interface and binary executable compiled from FORTRAN source code (HIFU beam. (2021). Available: http://limu.msu.ru/node/3555?language=en). It is designed for simulating high-intensity focused ultrasound (HIFU) fields generated by single-element transducers and annular arrays with propagation in flat-layered media that mimic biological tissues. Numerical models incorporated in the simulator include evolution-type equations, either the Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation or one-way Westervelt equation, for radially symmetric ultrasound beams in homogeneous and layered media with thermoviscous or power-law acoustic absorption. The software uses shock-capturing methods that allow for simulating strongly nonlinear acoustic fields with high-amplitude shocks. In this article, a general description of the software is given along with three representative simulation cases of ultrasound transducers and focusing conditions typical for therapeutic applications. The examples illustrate major nonlinear wave effects in HIFU fields including shock formation. Two examples simulate propagation in water, involving a single-element source (1-MHz frequency, 100-mm diameter, 90-mm radius of curvature) and a 16-element annular array (3-MHz frequency, 48-mm diameter, and 35-mm radius of curvature). The third example mimics the scenario of a HIFU treatment in a "water-muscle-kidney" layered medium using a source typical for abdominal HIFU applications (1.2-MHz frequency, 120-mm diameter, and radius of curvature). Linear, quasi-linear, and shock-wave exposure protocols are considered. It is intended that "HIFU beam" can be useful in teaching nonlinear acoustics; designing and characterizing high-power transducers; and developing exposure protocols for a wide range of therapeutic applications such as shock-based HIFU, boiling histotripsy, drug delivery, immunotherapy, and others.

Experimental Validation of k-Wave: Nonlinear Wave Propagation in Layered, Absorbing Fluid Media

Models of ultrasound propagation in biologically relevant media have applications in planning and verification of ultrasound therapies and computational dosimetry. To be effective, the models must be able to accurately predict both the spatial distribution and amplitude of the acoustic pressure. This requires that the models are validated in absolute terms, which for arbitrarily heterogeneous media should be performed by comparison with measurements of the acoustic field. In this article, simulations performed using the open-source k-Wave acoustics toolbox, with a measurement-based source definition, were quantitatively validated against measurements of acoustic pressure in water and layered absorbing fluid media. In water, the measured and simulated spatial-peak pressures agreed to within 3% under linear conditions and 7% under nonlinear conditions. After propagation through a planar or wedge-shaped glycerol-filled phantom, the difference in spatial-peak pressure was 8.5% and 10.7%, respectively. These differences are within or close to the expected uncertainty of the acoustic pressure measurement. The -6 dB width and length of the focus agreed to within 4% in all cases, and the focal positions were within 0.7 mm for the planar phantom and 1.2 mm for the wedge-shaped phantom. These results demonstrate that when the acoustic medium properties and geometry are well known, accurate quantitative predictions of the acoustic field can be made using k-Wave.

Numerical Simulations of the Nonlinear Interaction of a Bubble Cloud and a High Intensity Focused Ultrasound Field

We studied the effects of a small bubble cloud located at the pre-focal area of a high-intensity focused ultrasound field. Our objective is to show that bubbles can modify the bioeffects of an ultrasound treatment in muscle tissue. We model a three-dimensional ultrasound field in an idealized configuration of real operating conditions. Simulations are performed using a combined method based on the Khokhlov-Zabolotskaya-Kuznetsov equation, describing the ultrasound propagation, and a Rayleigh-Plesset equation, modeling the bubble oscillations. The nonlinear interaction of the ultrasound field and the bubble oscillations is considered. Results with and without bubbles for different void fractions of the cloud and different acoustic powers are compared. The cloud induces scattering, nonlinear distortion, and shielding of ultrasound, which increase the mechanical index in the pre-focal zone, shift the location, reduce the size, and modify the shape of the volume of tissue of high mechanical index values, and lower the pressure at the intended focus considerably. Although some hypothesis and parameters used in the models do not fit the real HIFU situations, the simulation results suggest that the effects caused by a bubble cloud located in the pre-focal area should be considered and monitored to ensure the safety of high-intensity focused ultrasound treatments.

Simulation of the Second-Harmonic Ultrasound Field in Heterogeneous Soft Tissue Using a Mixed-Domain Method

A mixed-domain method (MDM) dubbed frequency-specific MDM (FSMDM) is introduced for the simulation of the second-harmonic ultrasound field in weakly heterogeneous media. The governing equation for the second harmonics is derived based on the quasi-linear theory. The speed of sound, nonlinear coefficient, and attenuation coefficient are all spatially varying functions in the equation. The fundamental frequency pressure field is first solved by the FSMDM and it is subsequently used as the source term for the second-harmonics equation. This equation can be again solved by the FSMDM to rapidly obtain the second-harmonic pressure field. Five 2-D cases, including one with a realistic human tissue map, are studied to systematically verify the proposed method. Results from the previously developed transient MDM are used as the benchmark solutions. Comparisons show that the two methods give similar results for all cases. More importantly, the FSMDM has a crucial advantage over the transient MDM in that it can be two orders of magnitude faster.

The nonlinear ultrasound needle pulse

Recent work has established an analytical formulation of broadband fields which extend in the axial direction and converge to a narrow concentrated line. Those unique (needle) fields have their origins in an angular spectrum configuration in which the forward propagating wavenumber of the field ( kz ) is constant across any z plane for all of the propagated frequencies. A 3 MHz-based, finite amplitude distorted simulation of such a field is considered here in a water path scenario relevant to medical imaging. That nonlinear simulation had its focal features compared to those of a comparable Gaussian beam. The results suggest that the unique convergence of the needle pulse to a narrow but extended axial line in linear propagation is also inherited by higher harmonics in nonlinear propagation. Furthermore, the linear needle field's relatively short duration focal pulses, and the asymptotic declines of its radial profiles, also hold for the associated higher harmonics. Comparisons with the Gaussian field highlight some unique and potentially productive features of needle fields.

Numerical and Experimental Evaluation of High-Intensity Focused Ultrasound-Induced Lesions in Liver Tissue Ex Vivo

OBJECTIVES

Recent advances in the field of acoustics and piezoelectric and ultrasound transducers have led to new approaches to the diagnosis and treatment of certain diseases. One method of treatment with ultrasonic waves is high-intensity focused ultrasound (HIFU) treatment, which is a thermal therapeutic method used to treat malignant tumors. Although a variety of treatment-planning strategies using ultrasonic waves have been investigated, little clinical success has been achieved. Computational modeling is a powerful tool for predicting device performance.

METHODS

The heating induced by a concave transducer with operating powers of 85 and 135 W was studied, and the experimental results presented in this article verify its applicability. Numerical simulations of the nonlinear acoustic field were performed by using the Westervelt and Khokhlov-Zabolotskaya-Kuznetsov equations. Heat transfer was measured for the 2 operational powers, and the results were compared with ex vivo experimental results. In addition, thermal dose contours for both the simulation and experimental results were calculated to investigate the ablated area.

RESULTS

Good agreement was found between the experimental and numerical results. The results show that the average temperature deviations calculated at the focal point were 12.8% and 4.3% for transducer powers of 85 and 135 W, respectively.

CONCLUSIONS

This study provides guidance to HIFU practitioners in determining lesion size and identifying nonlinear effects that should be considered in HIFU procedures.

Simulation of nonlinear propagation of biomedical ultrasound using pzflex and the Khokhlov-Zabolotskaya-Kuznetsov Texas code

Nonlinear acoustics plays an important role in both diagnostic and therapeutic applications of biomedical ultrasound and a number of research and commercial software packages are available. In this manuscript, predictions of two solvers available in a commercial software package, pzflex, one using the finite-element-method (FEM) and the other a pseudo-spectral method, spectralflex, are compared with measurements and the Khokhlov-Zabolotskaya-Kuznetsov (KZK) Texas code (a finite-difference time-domain algorithm). The pzflex methods solve the continuity equation, momentum equation and equation of state where they account for nonlinearity to second order whereas the KZK code solves a nonlinear wave equation with a paraxial approximation for diffraction. Measurements of the field from a single element 3.3 MHz focused transducer were compared with the simulations and there was good agreement for the fundamental frequency and the harmonics; however the FEM pzflex solver incurred a high computational cost to achieve equivalent accuracy. In addition, pzflex results exhibited non-physical oscillations in the spatial distribution of harmonics when the amplitudes were relatively low. It was found that spectralflex was able to accurately capture the nonlinear fields at reasonable computational cost. These results emphasize the need to benchmark nonlinear simulations before using codes as predictive tools.

A framework for simulating ultrasound imaging based on first order nonlinear pressure-velocity relations

An ultrasound imaging framework modeled with the first order nonlinear pressure-velocity relations (NPVR) and implemented by a half-time staggered solution and pseudospectral method is presented in this paper. The framework is capable of simulating linear and nonlinear ultrasound propagation and reflections in a heterogeneous medium with different sound speeds and densities. It can be initialized with arbitrary focus, excitation and apodization for multiple individual channels in both 2D and 3D spatial fields. The simulated channel data can be generated using this framework, and ultrasound image can be obtained by beamforming the simulated channel data. Various results simulated by different algorithms are illustrated for comparisons. The root mean square (RMS) errors for each compared pulses are calculated. The linear propagation is validated by an angular spectrum approach (ASA) with a RMS error of 3% at the focal point for a 2D field, and Field II with RMS errors of 0.8% and 1.5% at the electronic and the elevation focuses for 3D fields, respectively. The accuracy for the NPVR based nonlinear propagation is investigated by comparing with the Abersim simulation for pulsed fields and with the nonlinear ASA for monochromatic fields. The RMS errors of the nonlinear pulses calculated by the NPVR and Abersim are respectively 2.4%, 7.4%, 17.6% and 36.6% corresponding to initial pressure amplitudes of 50kPa, 200kPa, 500kPa and 1MPa at the transducer. By increasing the sampling frequency for the strong nonlinearity, the RMS error for 1MPa initial pressure amplitude is reduced from 36.6% to 27.3%.

Modeling of wave propagation for medical ultrasound: a review

Numerical modeling of medical ultrasound has advanced tremendously in the past two decades. This opens up a great number of opportunities for medical ultrasound and associated technologies. Numerous new governing equations and algorithms have emerged and been applied to studying various medical ultrasound applications, including ultrasound imaging, photo-acoustic imaging, and therapeutic ultrasound. In addition, thanks to the rapid development of computers, modeling acoustic wave propagation in three-dimensional, large-scale domains has become a reality. This article will provide an indepth literature and technical review of recent progress on numerical modeling of medical ultrasound. Future challenges will also be discussed.

An improved wave-vector frequency-domain method for nonlinear wave modeling

In this paper, a recently developed wave-vector frequency-domain method for nonlinear wave modeling is improved and verified by numerical simulations and underwater experiments. Higher order numeric schemes are proposed that significantly increase the modeling accuracy, thereby allowing for a larger step size and shorter computation time. The improved algorithms replace the left-point Riemann sum in the original algorithm by the trapezoidal or Simpson's integration. Plane waves and a phased array were first studied to numerically validate the model. It is shown that the left-point Riemann sum, trapezoidal, and Simpson's integration have first-, second-, and third-order global accuracy, respectively. A highly focused therapeutic transducer was then used for experimental verifications. Short high-intensity pulses were generated. 2-D scans were conducted at a prefocal plane, which were later used as the input to the numerical model to predict the acoustic field at other planes. Good agreement is observed between simulations and experiments.

Simulation of nonlinear Westervelt equation for the investigation of acoustic streaming and nonlinear propagation effects

This study investigates the influence of blood flow on temperature distribution during high-intensity focused ultrasound (HIFU) ablation of liver tumors. A three-dimensional acoustic-thermal-hydrodynamic coupling model is developed to compute the temperature field in the hepatic cancerous region. The model is based on the nonlinear Westervelt equation, bioheat equations for the perfused tissue and blood flow domains. The nonlinear Navier-Stokes equations are employed to describe the flow in large blood vessels. The effect of acoustic streaming is also taken into account in the present HIFU simulation study. A simulation of the Westervelt equation requires a prohibitively large amount of computer resources. Therefore a sixth-order accurate acoustic scheme in three-point stencil was developed for effectively solving the nonlinear wave equation. Results show that focused ultrasound beam with the peak intensity 2470 W/cm(2) can induce acoustic streaming velocities up to 75 cm/s in the vessel with a diameter of 3 mm. The predicted temperature difference for the cases considered with and without acoustic streaming effect is 13.5 °C or 81% on the blood vessel wall for the vein. Tumor necrosis was studied in a region close to major vessels. The theoretical feasibility to safely necrotize the tumors close to major hepatic arteries and veins was shown.

Unmitigated numerical solution to the diffraction term in the parabolic nonlinear ultrasound wave equation

Various numerical algorithms have been developed to solve the Khokhlov-Kuznetsov-Zabolotskaya (KZK) parabolic nonlinear wave equation. In this work, a generalized time-domain numerical algorithm is proposed to solve the diffraction term of the KZK equation. This algorithm solves the transverse Laplacian operator of the KZK equation in three-dimensional (3D) Cartesian coordinates using a finite-difference method based on the five-point implicit backward finite difference and the five-point Crank-Nicolson finite difference discretization techniques. This leads to a more uniform discretization of the Laplacian operator which in turn results in fewer calculation gridding nodes without compromising accuracy in the diffraction term. In addition, a new empirical algorithm based on the LU decomposition technique is proposed to solve the system of linear equations obtained from this discretization. The proposed empirical algorithm improves the calculation speed and memory usage, while the order of computational complexity remains linear in calculation of the diffraction term in the KZK equation. For evaluating the accuracy of the proposed algorithm, two previously published algorithms are used as comparison references: the conventional 2D Texas code and its generalization for 3D geometries. The results show that the accuracy/efficiency performance of the proposed algorithm is comparable with the established time-domain methods.

Fourier continuation methods for high-fidelity simulation of nonlinear acoustic beams

On the basis of recently developed Fourier continuation (FC) methods and associated efficient parallelization techniques, this text introduces numerical algorithms that, due to very low dispersive errors, can accurately and efficiently solve the types of nonlinear partial differential equation (PDE) models of nonlinear acoustics in hundred-wavelength domains as arise in the simulation of focused medical ultrasound. As demonstrated in the examples presented in this text, the FC approach can be used to produce solutions to nonlinear acoustics PDEs models with significantly reduced discretization requirements over those associated with finite-difference, finite-element and finite-volume methods, especially in cases involving waves that travel distances that are orders of magnitude longer than their respective wavelengths. In these examples, the FC methodology is shown to lead to improvements in computing times by factors of hundreds and even thousands over those required by the standard approaches. A variety of one-and two-dimensional examples presented in this text demonstrate the power and capabilities of the proposed methodology, including an example containing a number of scattering centers and nonlinear multiple-scattering events.

A k-space method for moderately nonlinear wave propagation

A k-space method for moderately nonlinear wave propagation in absorptive media is presented. The Westervelt equation is first transferred into k-space via Fourier transformation, and is solved by a modified wave-vector time-domain scheme. The present approach is not limited to forward propagation or parabolic approximation. One- and two-dimensional problems are investigated to verify the method by comparing results to analytic solutions and finite-difference time-domain (FDTD) method. It is found that to obtain accurate results in homogeneous media, the grid size can be as little as two points per wavelength, and for a moderately nonlinear problem, the Courant-Friedrichs-Lewy number can be as large as 0.4. Through comparisons with the conventional FDTD method, the k-space method for nonlinear wave propagation is shown here to be computationally more efficient and accurate. The k-space method is then employed to study three-dimensional nonlinear wave propagation through the skull, which shows that a relatively accurate focusing can be achieved in the brain at a high frequency by sending a low frequency from the transducer. Finally, implementations of the k-space method using a single graphics processing unit shows that it required about one-seventh the computation time of a single-core CPU calculation.

Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method

The simulation of nonlinear ultrasound propagation through tissue realistic media has a wide range of practical applications. However, this is a computationally difficult problem due to the large size of the computational domain compared to the acoustic wavelength. Here, the k-space pseudospectral method is used to reduce the number of grid points required per wavelength for accurate simulations. The model is based on coupled first-order acoustic equations valid for nonlinear wave propagation in heterogeneous media with power law absorption. These are derived from the equations of fluid mechanics and include a pressure-density relation that incorporates the effects of nonlinearity, power law absorption, and medium heterogeneities. The additional terms accounting for convective nonlinearity and power law absorption are expressed as spatial gradients making them efficient to numerically encode. The governing equations are then discretized using a k-space pseudospectral technique in which the spatial gradients are computed using the Fourier-collocation method. This increases the accuracy of the gradient calculation and thus relaxes the requirement for dense computational grids compared to conventional finite difference methods. The accuracy and utility of the developed model is demonstrated via several numerical experiments, including the 3D simulation of the beam pattern from a clinical ultrasound probe.

Acoustic shock wave propagation in a heterogeneous medium: a numerical simulation beyond the parabolic approximation

Numerical simulation of nonlinear acoustics and shock waves in a weakly heterogeneous and lossless medium is considered. The wave equation is formulated so as to separate homogeneous diffraction, heterogeneous effects, and nonlinearities. A numerical method called heterogeneous one-way approximation for resolution of diffraction (HOWARD) is developed, that solves the homogeneous part of the equation in the spectral domain (both in time and space) through a one-way approximation neglecting backscattering. A second-order parabolic approximation is performed but only on the small, heterogeneous part. So the resulting equation is more precise than the usual standard or wide-angle parabolic approximation. It has the same dispersion equation as the exact wave equation for all forward propagating waves, including evanescent waves. Finally, nonlinear terms are treated through an analytical, shock-fitting method. Several validation tests are performed through comparisons with analytical solutions in the linear case and outputs of the standard or wide-angle parabolic approximation in the nonlinear case. Numerical convergence tests and physical analysis are finally performed in the fully heterogeneous and nonlinear case of shock wave focusing through an acoustical lens.

SIMULATION OF THREE-DIMENSIONAL NONLINEAR FIELDS OF ULTRASOUND THERAPEUTIC ARRAYS

A novel numerical model was developed to simulate three-dimensional nonlinear fields generated by high intensity focused ultrasound (HIFU) arrays. The model is based on the solution to the Westervelt equation; the developed algorithm makes it possible to model nonlinear pressure fields of periodic waves in the presence of shock fronts localized near the focus. The role of nonlinear effects in a focused beam of a two-dimensional array was investigated in a numerical experiment in water. The array consisting of 256 elements and intensity range on the array elements of up to 10 W/cm(2) was considered. The results of simulations have shown that for characteristic intensity outputs of modern HIFU arrays, nonlinear effects play an important role and shock fronts develop in the pressure waveforms at the focus.

Nonlinear and diffraction effects in propagation of N-waves in randomly inhomogeneous moving media

Finite amplitude acoustic wave propagation through atmospheric turbulence is modeled using a Khokhlov-Zabolotskaya-Kuznetsov (KZK)-type equation. The equation accounts for the combined effects of nonlinearity, diffraction, absorption, and vectorial inhomogeneities of the medium. A numerical algorithm is developed which uses a shock capturing scheme to reduce the number of temporal grid points. The inhomogeneous medium is modeled using random Fourier modes technique. Propagation of N-waves through the medium produces regions of focusing and defocusing that is consistent with geometrical ray theory. However, differences up to ten wavelengths are observed in the locations of fist foci. Nonlinear effects are shown to enhance local focusing, increase the maximum peak pressure (up to 60%), and decrease the shock rise time (about 30 times). Although the peak pressure increases and the rise time decreases in focal regions, statistical analysis across the entire wavefront at a distance 120 wavelengths from the source indicates that turbulence: decreases the mean time-of-flight by 15% of a pulse duration, decreases the mean peak pressure by 6%, and increases the mean rise time by almost 100%. The peak pressure and the arrival time are primarily governed by large scale inhomogeneities, while the rise time is also sensitive to small scales.

Short-range shock formation and coalescence in numerical simulation of broadband noise propagation

The number of jet and rocket noise studies has increased in recent years as researchers have sought to better understand aeroacoustic source and radiation characteristics. Although jet and rocket noise is finite-amplitude in nature, little is known about the existence of shock formation and coalescence close to the source. A numerical experiment is performed to propagate finite-amplitude noise and determine the extent of the nonlinearity over short distances with spherical spreading. The noise is filtered to have a haystack shape in the frequency domain, as is typical of such sources. The effect of the nonlinearity is compared in both the temporal and frequency domains as a function of distance. Additionally, the number of zero-crossings and overall sound pressure level is compared at several distances. The results indicate that the center frequency plays a particularly important role in the amount of coalescence and spectral redistribution that occurs. The general applicability of these results to actual near-field finite-amplitude jet and rocket noise experiments is also presented.

Transmit beams adapted to reverberation noise suppression using dual-frequency SURF imaging

A method that uses dual-frequency pulse complexes of widely separated frequency bands to suppress noise caused by multiple scattering or multiple reflections in medical ultrasound imaging is presented. The excitation pulse complexes are transmitted to generate a second order ultrasound field (SURF) imaging synthetic transmit beam. This beam has reduced amplitude near the transducer, which illustrates the multiple scattering suppression ability of the imaging method. Field simulations solving a nonlinear wave equation are used to calculate SURF imaging beams, which are compared with beams for pulse inversion (PI) and fundamental imaging. In addition, a combined SURF and PI beam generation is described and compared with the beams mentioned above. A quality ratio, relating the energy within the near-field to that within the imaging region, is defined and used to score the multiple scattering and multiple reflection suppression abilities when imaging with the different beams. The realized combined SURF-PI beam scores highest, followed by SURF, PI (that score equally well), and the fundamental. The amplitude in the imaging region and therefore also the SNR is highest for the fundamental followed by SURF, PI, and SURF-PI. The work hence indicates that when substituting PI for SURF, one may trade increased SNR into use of increased imaging frequencies without loss of multiple scattering and multiple reflection noise suppression.

Computational aspects in high intensity ultrasonic surgery planning

Therapeutic ultrasound treatment planning is discussed and computational aspects regarding it are reviewed. Nonlinear ultrasound simulations were solved with a combined frequency domain Rayleigh and KZK model. Ultrasonic simulations were combined with thermal simulations and were used to compute heating of muscle tissue in vivo for four different focused ultrasound transducers. The simulations were compared with measurements and good agreement was found for large F-number transducers. However, at F# 1.9 the simulated rate of temperature rise was approximately a factor of 2 higher than the measured ones. The power levels used with the F# 1 transducer were too low to show any nonlinearity. The simulations were used to investigate the importance of nonlinarities generated in the coupling water, and also the importance of including skin in the simulations. Ignoring either of these in the model would lead to larger errors. Most notably, the nonlinearities generated in the water can enhance the focal temperature by more than 100%. The simulations also demonstrated that pulsed high power sonications may provide an opportunity to significantly (up to a factor of 3) reduce the treatment time. In conclusion, nonlinear propagation can play an important role in shaping the energy distribution during a focused ultrasound treatment and it should not be ignored in planning. However, the current simulation methods are accurate only with relatively large F-numbers and better models need to be developed for sharply focused transducers.

Numerical investigation of the properties of nonlinear acoustical vortices through weakly heterogeneous media

Acoustical vortices (AV) are the acoustical equivalent of optical vortices (OV) that are a key feature in the discipline of singular optics. For linear waves, OV and AV possess the same properties. But as nonlinearities are different in optics and acoustics, the nonlinear behavior of these structures has to be different. In this paper, a numerical investigation of the three-dimensional (3D) nonlinear propagation of acoustical vortices through homogeneous or heterogeneous media is reported. First, an original numerical method is described and compared to existing ones. Then, it is used to study the dynamics of AV in a nonlinear regime. The nonlinear properties of acoustical vortices in a homogeneous medium are investigated. It is shown that shock waves can be produced during propagation, leading to an interesting spatiotemporal wave field with an azimuthal shock. The dynamics of the topological charge, intrinsic property of AV or OV, is studied in the nonlinear regime through different focusing lenses. Inversion of the topological charge is observed if the AV propagates through a 1D focusing medium (cylindrical lens), while the charge remains constant if the medium is 2D (spherical lens). These last results already observed in linear optics are generalized here to the nonlinear behavior through the investigation of harmonics which show the same behavior as the fundamental with respect to inversion.