Benchmark problems for transcranial ultrasound simulation: Intercomparison of compressional wave models.

Computational models of acoustic wave propagation are frequently used in transcranial ultrasound therapy, for example, to calculate the intracranial pressure field or to calculate phase delays to correct for skull distortions. To allow intercomparison between the different modeling tools and techniques used by the community, an international working group was convened to formulate a set of numerical benchmarks. Here, these benchmarks are presented, along with intercomparison results.

On-the-Fly Calculation of Time-Averaged Acoustic Intensity in Time-Domain Ultrasound Simulations Using a k-Space Pseudospectral Method.

This article presents a method to calculate the average acoustic intensity during ultrasound simulation using a new approach that exploits compression of intermediate results. One of the applications of high-intensity focused ultrasound (HIFU) simulations is the calculation of the thermal dose, which indicates the amount of tissue destroyed using a state-of-the-art k-space pseudospectral method. The thermal simulation is preceded by the calculation of the average intensity within the acoustic simulation.

Nonlinear ultrasound simulation in an axisymmetric coordinate system using a k-space pseudospectral method.

A full-wave model for nonlinear ultrasound propagation through a heterogeneous and absorbing medium in an axisymmetric coordinate system is developed. The model equations are solved using a nonstandard or k-space pseudospectral time domain method. Spatial gradients in the axial direction are calculated using the Fourier collocation spectral method, and spatial gradients in the radial direction are calculated using discrete trigonometric transforms.

Broadband All-Optical Plane-Wave Ultrasound Imaging System Based on a Fabry-Perot Scanner.

A broadband all-optical plane-wave ultrasound imaging system for high-resolution 3-D imaging of biological tissues is presented. The system is based on a planar Fabry-Perot (FP) scanner for ultrasound detection and the photoacoustic generation of ultrasound in a carbon-nanotube-polydimethylsiloxane (CNT-PDMS) composite film. The FP sensor head was coated with the CNT-PDMS film which acts as an ultrasound transmitting layer for pulse-echo imaging.

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.

Representing arbitrary acoustic source and sensor distributions in Fourier collocation methods.

Accurately representing acoustic source distributions is an important part of ultrasound simulation. This is challenging for grid-based collocation methods when such distributions do not coincide with the grid points, for instance when the source is a curved, two-dimensional surface embedded in a three-dimensional domain. Typically, grid points close to the source surface are defined as source points, but this can result in "staircasing" and substantial errors in the resulting acoustic fields.

The Effect of Tissue Physiological Variability on Transurethral Ultrasound Therapy of the Prostate.

Therapeutic ultrasound is an investigational modality which could potentially be used for minimally invasive treatment of prostate cancer. Computational simulations were used to study the effect of natural physiological variations in tissue parameters on the efficacy of therapeutic ultrasound treatment in the prostate. The simulations were conducted on a clinical ultrasound therapy system using patient computed tomography (CT) data.

Full Modeling of High-Intensity Focused Ultrasound and Thermal Heating in the Kidney Using Realistic Patient Models.

Objective: High-intensity focused ultrasound (HIFU) therapy can be used for noninvasive treatment of kidney (renal) cancer, but the clinical outcomes have been variable. In this study, the efficacy of renal HIFU therapy was studied using a nonlinear acoustic and thermal simulations in three patients.

Methods: The acoustic simulations were conducted with and without refraction in order to investigate its effect on the shape, size, and pressure distribution at the focus.

Transurethral ultrasound therapy of the prostate in the presence of calcifications: A simulation study.

Purpose: Transurethral ultrasound therapy is an investigational treatment modality which could potentially be used for the localized treatment of prostate cancer. One of the limiting factors of this therapy is prostatic calcifications. These attenuate and reflect ultrasound and thus reduce the efficacy of the heating.

Rapid calculation of acoustic fields from arbitrary continuous-wave sources.

A Green's function solution is derived for calculating the acoustic field generated by phased array transducers of arbitrary shape when driven by a single frequency continuous wave excitation with spatially varying amplitude and phase. The solution is based on the Green's function for the homogeneous wave equation expressed in the spatial frequency domain or k-space. The temporal convolution integral is solved analytically, and the remaining integrals are expressed in the form of the spatial Fourier transform.

Full Modeling of High-Intensity Focused Ultrasound and Thermal Heating in the Kidney Using Realistic Patient Models.

Objective: High-intensity focused ultrasound (HIFU) therapy can be used for noninvasive treatment of kidney (renal) cancer, but the clinical outcomes have been variable. In this study, the efficacy of renal HIFU therapy was studied using nonlinear acoustic and thermal simulations in three patients.

Methods: The acoustic simulations were conducted with and without refraction in order to investigate its effect on the shape, size, and pressure distribution at the focus.

Nonlinear 3-D simulation of high-intensity focused ultrasound therapy in the Kidney.

Kidney cancer is a severe disease which can be treated non-invasively using high-intensity focused ultrasound (HIFU) therapy. However, tissue in front of the transducer and the deep location of kidney can cause significant losses to the efficiency of the treatment. The effect of attenuation, refraction and reflection due to different tissue types on HIFU therapy of the kidney was studied using a nonlinear ultrasound simulation model.

Control of broadband optically generated ultrasound pulses using binary amplitude holograms.

In this work, the use of binary amplitude holography is investigated as a mechanism to focus broadband acoustic pulses generated by high peak-power pulsed lasers. Two algorithms are described for the calculation of the binary holograms; one using ray-tracing, and one using an optimization based on direct binary search. It is shown using numerical simulations that when a binary amplitude hologram is excited by a train of laser pulses at its design frequency, the acoustic field can be focused at a pre-determined distribution of points, including single and multiple focal points, and line and square foci.

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.