A comparison of methods for focusing the field of a HIFU array transducer through human ribs

Aberration correction in abdominal histotripsy

In transabdominal histotripsy, ultrasound pulses are focused on the body to noninvasively destroy soft tissues via cavitation. However, the ability to focus is limited by phase aberration, or decorrelation of the ultrasound pulses due to spatial variation in the speed of sound throughout heterogeneous tissue. Phase aberration shifts, broadens, and weakens the focus, thereby reducing the safety and efficacy of histotripsy therapy. This paper reviews and discusses aberration effects in histotripsy and in related therapeutic ultrasound techniques (e.g., high intensity focused ultrasound), with an emphasis on aberration by soft tissues. Methods for aberration correction are reviewed and can be classified into two groups: model-based methods, which use segmented images of the tissue as input to an acoustic propagation model to predict and compensate phase differences, and signal-based methods, which use a receive-capable therapy array to detect phase differences by sensing acoustic signals backpropagating from the focus. The relative advantages and disadvantages of both groups of methods are discussed. Importantly, model-based methods can correct focal shift, while signal-based methods can restore substantial focal pressure, suggesting that both methods should be combined in a 2-step approach. Aberration correction will be critical to improving histotripsy treatments and expanding the histotripsy treatment envelope to enable non-invasive, non-thermal histotripsy therapy for more patients.

Soft Tissue Aberration Correction for Histotripsy Using Acoustic Emissions From Cavitation Cloud Nucleation and Collapse

OBJECTIVE

Phase aberration from soft tissue limits the efficacy of histotripsy, a therapeutic ultrasound technique based on acoustic cavitation. Previous work has shown that the acoustic emissions from cavitation can serve as "point sources" for aberration correction (AC). This study compared the efficacy of soft tissue AC for histotripsy using acoustic cavitation emissions (ACE) from bubble cloud nucleation and collapse.

METHODS

A 750-kHz, receive-capable histotripsy array was pulsed to generate cavitation in ex vivo porcine liver through an intervening abdominal wall. Received ACE signals were used to determine the arrival time differences to the focus and compute corrective delays. Corrections from single pulses and from the median of multiple pulses were tested.

DISCUSSION

On average, ACE AC obtained 96% ± 3% of the pressure amplitude obtained by hydrophone-based correction (compared with 71% ± 5% without AC). Both nucleation- and collapse-based corrections obtained >96% of the hydrophone-corrected pressure when using medians of ≥10 pulses. When using single-pulse corrections, nucleation obtained a range of 49%-99% of the hydrophone-corrected pressure, while collapse obtained 95%-99%.

CONCLUSION

The results suggest that (i) ACE AC can recover nearly all pressure amplitude lost owing to soft tissue aberration and that (ii) the collapse signal permits robust AC using a small number of pulses.

Tissue specific considerations in implementing high intensity focussed ultrasound under magnetic resonance imaging guidance

High-intensity focused ultrasound can ablate a target permanently, leaving tissues through which it passes thermally unaffected. When delivered under magnetic resonance (MR) imaging guidance, the change in tissue relaxivity on heating is used to monitor the temperatures achieved. Different tissue types in the pre-focal beam path result in energy loss defined by their individual attenuation coefficients. Furthermore, at interfaces with different acoustic impedances the beam will be both reflected and refracted, changing the position of the focus. For complex interfaces this effect is exacerbated. Moreover, blood vessels proximal to the focal region can dissipate heat, altering the expected region of damage. In the target volume, the temperature distribution depends on the thermal conductivity (or diffusivity) of the tissue and its heat capacity. These are different for vascular tissues, water and fat containing tissues and bone. Therefore, documenting the characteristics of the pre-focal and target tissues is critical for effective delivery of HIFU. MR imaging provides excellent anatomic detail and characterization of soft tissue components. It is an ideal modality for real-time planning and monitoring of HIFU ablation, and provides non-invasive temperature maps. Clinical applications involve soft-tissue (abdomino-pelvic applications) or bone (brain applications) pre-focally and at the target (soft-tissue tumors and bone metastases respectively). This article addresses the technical difficulties of delivering HIFU effectively when vascular tissues, densely cellular tissues, fat or bone are traversed pre-focally, and the clinical applications that target these tissues. The strengths and limitations of MR techniques used for monitoring ablation in these tissues are also discussed.

Development of a Numerical Model of High-Intensity Focused Ultrasound Treatment in Mobile and Elastic Organs: Application to a Beating Heart

High-intensity focused ultrasound (HIFU) is a promising method used to treat cardiac arrhythmias, as it can induce lesions at a distance throughout myocardium thickness. Numerical modeling is commonly used for ultrasound probe development and optimization of HIFU treatment strategies. This study was aimed at describing a numerical method to simulate HIFU thermal ablation in elastic and mobile heart models. The ultrasound pressure field is computed on a 3-D orthonormal grid using the Rayleigh integral method, and the attenuation is calculated step by step between cells. The temperature distribution is obtained by resolution of the bioheat transfer equation on a 3-D non-orthogonally structured curvilinear grid using the finite-volume method. The simulation method is applied on two regions of the heart (atrioventricular node and ventricular apex) to compare the thermal effects of HIFU ablation depending on deformation, motion type and amplitude. The atrioventricular node requires longer sonication than the ventricular apex to reach the same lesion volume. Motion considerably influences treatment duration, lesion shape and distribution in cardiac HIFU treatment. These results emphasize the importance of considering local motion and deformation in numerical studies to define efficient and accurate treatment strategies.

Complementary Acoustic Metamaterial for Penetrating Aberration Layers

Impedance-matched acoustic materials were developed to improve ultrasound penetration through the aberration layer. The traditional ultrasound layer matching material is called a couplant, which can only enhance ultrasound transmission to soft biological media such as the cartilage and muscle but cannot penetrate hard media such as the bone. Here, we propose a phase-modulated complementary acoustic metamaterial based on the principle of impedance matching, which enables ultrasound to penetrate the bone, and use the equivalent parameter technology of acoustic metamaterials for parameter design. Ultrasonic layer adjustment is performed through 3D printing and corrects bone aberrations. Several configurations were investigated through numerical simulations and experiments in non-reflecting tanks. Specifically, the bone matching layer can be optimally designed for a specific bone thickness and a specific operating frequency of the ultrasound probe, thereby amplifying the ultrasound to penetrate the matching layer and bone. The experimental and simulation results show that the proposed acoustic metamaterial can improve the transmission efficiency of ultrasound through the aberration layer.

"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.

A fast full-wave solver for calculating ultrasound propagation in the body

Therapeutic ultrasound is a promising non-invasive method for inducing various beneficial biological effects in the human body. In cancer treatment applications, high-power ultrasound is focused at a target tissue volume to ablate the malignant tumour. The success of the procedure depends on the ability to accurately focus ultrasound and destroy the target tissue volume through coagulative necrosis whilst preserving the surrounding healthy tissue. Patient-specific treatment planning strategies are therefore being developed to increase the efficacy of such therapies, while reducing any damage to healthy tissue. These strategies require to use high-performance computing methods to solve ultrasound wave propagation in the body quickly and accurately. For realistic clinical scenarios, all numerical methods which employ volumetric meshes require several hours or days to solve the full-wave propagation on a computer cluster. The boundary element method (BEM) is an efficient approach for modelling the wave field because only the boundaries of the hard and soft tissue regions require discretisation. This paper presents a multiple-domain BEM formulation with a novel preconditioner for solving the Helmholtz transmission problem (HTP). This new formulation is efficient at high-frequencies and where high-contrast materials are present. Numerical experiments are performed to solve the HTP in multiple domains comprising: (i) human ribs, an idealised abdominal fat layer and liver tissue, (ii) a human kidney with a perinephric fat layer, exposed to the acoustic field generated by a high-intensity focused ultrasound (HIFU) array transducer. The time required to solve the equations associated with these problems on a single workstation is of the order of minutes. These results demonstrate the great potential of this new BEM formulation for accurately and quickly solving ultrasound wave propagation problems in large anatomical domains which is essential for developing treatment planning strategies.

The Optimization of Transcostal Phased Array Refocusing Using the Semidefinite Relaxation Method

Tumors in organs partially obscured by the rib cage represent a challenge for high-intensity focused ultrasound (HIFU) therapy. The ribs distort the HIFU beams in a manner that reduces the focusing gain at the target, which could result in treatment-limiting collateral damage. In fact, skin burns are a common complication during the ablation of hepatic tumors. This problem can be addressed by employing optimal refocusing algorithms that are designed to achieve a specified focusing gain at the target while controlling the exposure to the ribs in the path of the HIFU beam. However, previously proposed optimal refocusing algorithms did not allow for the controlled transmission through the ribs. In this article, we introduce a new approach for refocusing that can more efficiently steer power toward the target while limiting the power deposition on the ribs. The approach utilizes the semidefinite relaxation (SDR) technique to approximate the original (nonconvex) optimization problem. An important advantage of the SDR-based method over previously proposed optimization methods is the control of the side lobes in the focal plane. The method also allows for specifying an acceptable level of exposure to the ribs. Simulation results using a 1-MHz spherical concave phased array focused on an inhomogeneous medium are presented to demonstrate the performance of the SDR refocusing approach. A finite-difference time-domain propagation model was used to model the propagation in the inhomogeneous tissues, including the ribs. Temperature simulations based on the inhomogeneous transient bioheat transfer equation (tBHTE) demonstrate the significance of the improvements in the focusing gain when using the limited power deposition (LPD) method. The results also demonstrate that the LPD method yields well-behaved array excitation vectors, realizable by currently existing drivers.

Focused Ultrasound Hyperthermia for Targeted Drug Release from Thermosensitive Liposomes: Results from a Phase I Trial

Purpose To demonstrate the feasibility and safety of using focused ultrasound planning models to determine the treatment parameters needed to deliver volumetric mild hyperthermia for targeted drug delivery without real-time thermometry. Materials and Methods This study was part of the Targeted Doxorubicin, or TARDOX, phase I prospective trial of focused ultrasound-mediated, hyperthermia-triggered drug delivery to solid liver tumors ( ClinicalTrials.gov identifier NCT02181075). Ten participants (age range, 49-68 years; average age, 60 years; four women) were treated from March 2015 to March 2017 by using a clinically approved focused ultrasound system to release doxorubicin from lyso-thermosensitive liposomes. Ultrasonic heating of target tumors (treated volume: 11-73 cm³ [mean ± standard deviation, 50 cm³ ± 26]) was monitored in six participants by using a minimally invasive temperature sensor; four participants were treated without real-time thermometry. For all participants, CT images were used with a patient-specific hyperthermia model to define focused ultrasound treatment plans. Feasibility was assessed by comparing model-prescribed focused ultrasound powers to those implemented for treatment. Safety was assessed by evaluating MR images and biopsy specimens for evidence of thermal ablation and monitoring adverse events. Results The mean difference between predicted and implemented treatment powers was -0.1 W ± 17.7 (n = 10). No evidence of focused ultrasound-related adverse effects, including thermal ablation, was found. Conclusion In this 10-participant study, the authors confirmed the feasibility of using focused ultrasound-mediated hyperthermia planning models to define treatment parameters that safely enabled targeted, noninvasive drug delivery to liver tumors while monitored with B-mode guidance and without real-time thermometry. Published under a CC BY 4.0 license. Online supplemental material is available for this article. See also the editorial by Dickey and Levi-Polyachenko in this issue.

Validation of hybrid angular spectrum acoustic and thermal modelling in phantoms

In focused ultrasound (FUS) thermal ablation of diseased tissue, acoustic beam and thermal simulations enable treatment planning and optimization. In this study, a treatment-planning methodology that uses the hybrid angular spectrum (HAS) method and the Pennes' bioheat equation (PBHE) is experimentally validated in homogeneous tissue-mimicking phantoms. Simulated three-dimensional temperature profiles are compared to volumetric MR thermometry imaging (MRTI) of FUS sonications in the phantoms, whose acoustic and thermal properties are independently measured. Additionally, Monte Carlo (MC) uncertainty analysis is performed to quantify the effect of tissue property uncertainties on simulation results. The mean error between simulated and experimental spatiotemporal peak temperature rise was +0.33°C (+6.9%). Despite this error, the experimental temperature rise fell within the expected uncertainty of the simulation, as determined by the MC analysis. The average errors of the simulated transverse and longitudinal full width half maximum (FWHM) of the profiles were -1.9% and 7.5%, respectively. A linear regression and local sensitivity analysis revealed that simulated temperature amplitude is more sensitive to uncertainties in simulation inputs than in the profile width and shape. Acoustic power, acoustic attenuation and thermal conductivity had the greatest impact on peak temperature rise uncertainty; thermal conductivity and volumetric heat capacity had the greatest impact on FWHM uncertainty. This study validates that using the HAS and PBHE method can adequately predict temperature profiles from single sonications in homogeneous media. Further, it informs the need to accurately measure or predict patient-specific properties for improved treatment planning of ablative FUS surgeries.

Self-adaptive ultrasonic beam amplifiers: application to transcostal shock wave therapy

Ultrasound shock wave therapy is increasingly used for non-invasive surgery. It requires the focusing of very high pressure amplitude in precisely controlled focal spots. In transcostal therapy of the heart or the liver, the high impedance mismatch between the bones and surrounding tissues gives rise to strong aberrations and attenuation of the therapeutic wavefront, with potential risks of injury at the tissue-bone interface. An adaptive propagation of the ultrasonic beam through the intercostal spaces would be required. Several solutions have been developed so far, but they require a prior knowledge of the patient's anatomy or an invasive calibration process, not applicable in clinic. Here, we develop a non-invasive adaptive focusing method for ultrasound therapy through the ribcage using a time reversal cavity (TRC) acting as an ultrasonic beam amplifier. This method is based on ribcage imaging through the TRC and a projection orthogonally to the strongest identified reflectors. The focal pressure of our device was improved by up to 30% using such self-adaptive processing, without degrading the focal spots size and shape. This improvement allowed lesion formation in an Ultracal^® phantom through a ribcage without invasive calibration of the device. This adaptive method could be particularly interesting to improve the efficiency and the safety of pulsed cavitational therapy of the heart or the liver.

Automated geometric optimization for robotic HIFU treatment of liver tumors

BACKGROUND

High intensity focused ultrasound (HIFU) represents a non-invasive method for the destruction of cancerous tissue within the body. Heating of targeted tissue by focused ultrasound transducers results in the creation of ellipsoidal lesions at the target site, the locations of which can have a significant impact on treatment outcomes. Towards this end, this work describes a method for the optimization of lesion positions within arbitrary tumors, with specific anatomical constraints.

MATERIALS & METHODS

A force-based optimization framework was extended to the case of arbitrary tumor position and constrained orientation. Analysis of the approximate reachable treatment volume for the specific case of treatment of liver tumors was performed based on four transducer configurations and constraint conditions derived. Evaluation was completed utilizing simplified spherical and ellipsoidal tumor models and randomly generated tumor volumes. The total volume treated, lesion overlap and healthy tissue ablated was evaluated. Two evaluation scenarios were defined and optimized treatment plans assessed.

RESULTS

The optimization framework resulted in improvements of up to 10% in tumor volume treated, and reductions of up to 20% in healthy tissue ablated as compared to the standard lesion rastering approach. Generation of optimized plans proved feasible for both sub- and intercostally located tumors.

CONCLUSIONS

This work describes an optimized method for the planning of lesion positions during HIFU treatment of liver tumors. The approach allows the determination of optimal lesion locations and orientations, and can be applied to arbitrary tumor shapes and sizes.

Numerical and Experimental Study of Mechanisms Involved in Boiling Histotripsy

The aim of boiling histotripsy is to mechanically fractionate tissue as an alternative to thermal ablation for therapeutic applications. In general, the shape of a lesion produced by boiling histotripsy is tadpole like, consisting of a head and a tail. Although many studies have demonstrated the efficacy of boiling histotripsy for fractionating solid tumors, the exact mechanisms underpinning this phenomenon are not yet well understood, particularly the interaction of a boiling vapor bubble with incoming incident shockwaves. To investigate the mechanisms involved in boiling histotripsy, a high-speed camera with a passive cavitation detection system was used to observe the dynamics of bubbles produced in optically transparent tissue-mimicking gel phantoms exposed to the field of a 2.0-MHz high-intensity focused ultrasound (HIFU) transducer. We observed that boiling bubbles were generated in a localized heated region and cavitation clouds were subsequently induced ahead of the expanding bubble. This process was repeated with HIFU pulses and eventually resulted in a tadpole-shaped lesion. A simplified numerical model describing the scattering of the incident ultrasound wave by a vapor bubble was developed to help interpret the experimental observations. Together with the numerical results, these observations suggest that the overall size of a lesion induced by boiling histotripsy is dependent on the sizes of (i) the heated region at the HIFU focus and (ii) the backscattered acoustic field by the original vapor bubble.

Thermal ablation of a confluent lesion in the porcine kidney with a clinically available MR-HIFU system

The incidence of small renal masses (SRMs) sized  <4 cm has increased over the decades (as co-findings/or due to introduction of cross sectional imaging). Currently, partial nephrectomy (PN) or watchful waiting is advised in these patients. Ultimately, 80-90% of these SRMs require surgical treatment and PN is associated with a 15% complication rate. In this aging population, with possible comorbidities and poor health condition, both PN and watchful waiting are non-ideal treatment options. This resulted in an increased need for early, non-invasive treatment strategies such as MR-guided high intensity focused ultrasound (MR-HIFU). (i) To investigate the feasibility of creating a confluent lesion in the kidney using respiratory-gated MR-HIFU under clinical conditions in a pre-clinical study and (ii) to evaluate the reproducibility of the MR-HIFU ablation strategy. Healthy pigs (n  =  10) under general anesthesia were positioned on a clinical MR-HIFU system with integrated cooling. A honeycomb pattern of seven overlapping ablation cells (4  ×  4  ×  10 mm³, 450 W, <30 s) was ablated successively in the cortex of the porcine kidney. Both MR thermometry and acoustic energy delivery were respiratory gated using a pencil beam navigator on the contralateral kidney. The non-perfused volume (NPV) was visualized after the last sonication by contrast-enhanced (CE) T ₁-weighted MR (T ₁ w) imaging. Cell viability staining was performed to visualize the extent of necrosis.

RESULTS

a median NPV of 0.62 ml was observed on CE-T ₁ w images (IQR 0.58-1.57 ml, range 0.33-2.75 ml). Cell viability staining showed a median damaged volume of 0.59 ml (IQR 0.24-1.35 ml, range 0-4.1 ml). Overlooking of the false rib, shivering of the pig, and too large depth combined with a large heat-sink effect resulted in insufficient heating in 4 cases. The NPV and necrosed volume were confluent in all cases in which an ablated volume could be observed. Our results demonstrated the feasibility of creating a confluent volume of ablated kidney cortical tissue in vivo with MR-HIFU on a clinically available system using respiratory gating and near-field cooling and showed its reproducibility.

Extending the Utility of the Parabolic Approximation in Medical Ultrasound Using Wide-Angle Diffraction Modeling

Wide-angle parabolic models are commonly used in geophysics and underwater acoustics but have seen little application in medical ultrasound. Here, a wide-angle model for continuous-wave high-intensity ultrasound beams is derived, which approximates the diffraction process more accurately than the commonly used Khokhlov-Zabolotskaya-Kuznetsov (KZK) equation without increasing implementation complexity or computing time. A method for preventing the high spatial frequencies often present in source boundary conditions from corrupting the solution is presented. Simulations of shallowly focused axisymmetric beams using both the wide-angle and standard parabolic models are compared to assess the accuracy with which they model diffraction effects. The wide-angle model proposed here offers improved focusing accuracy and less error throughout the computational domain than the standard parabolic model, offering a facile method for extending the utility of existing KZK codes.

Metamaterial bricks and quantization of meta-surfaces

Controlling acoustic fields is crucial in diverse applications such as loudspeaker design, ultrasound imaging and therapy or acoustic particle manipulation. The current approaches use fixed lenses or expensive phased arrays. Here, using a process of analogue-to-digital conversion and wavelet decomposition, we develop the notion of quantal meta-surfaces. The quanta here are small, pre-manufactured three-dimensional units-which we call metamaterial bricks-each encoding a specific phase delay. These bricks can be assembled into meta-surfaces to generate any diffraction-limited acoustic field. We apply this methodology to show experimental examples of acoustic focusing, steering and, after stacking single meta-surfaces into layers, the more complex field of an acoustic tractor beam. We demonstrate experimentally single-sided air-borne acoustic levitation using meta-layers at various bit-rates: from a 4-bit uniform to 3-bit non-uniform quantization in phase. This powerful methodology dramatically simplifies the design of acoustic devices and provides a key-step towards realizing spatial sound modulators.

Rapid full-wave phase aberration correction method for transcranial high-intensity focused ultrasound therapies

Background

Non-invasive high-intensity focused ultrasound (HIFU) can be used to treat a variety of disorders, including those in the brain. However, the differences in acoustic properties between the skull and the surrounding soft tissue cause aberrations in the path of the ultrasonic beam, hindering or preventing treatment.

Methods

We present a method for correcting these aberrations that is fast, full-wave, and allows for corrections at multiple treatment locations. The method is simulation-based: an acoustic model is built based on high-resolution CT scans, and simulations are performed using the hybrid angular spectrum (HAS) method to determine the phases needed for correction.

Results

Computation of corrections for clinically applicable resolutions can be achieved in approximately 15 min. Experimental results with a plastic model designed to mimic the aberrations caused by the skull show that the method can recover 95 % of the peak pressure obtained using hydrophone-based time-reversal methods. Testing using an ex vivo human skull flap resulted in recovering up to 70 % of the peak pressure at the focus and 61 % when steering (representing, respectively, a 1.52- and 1.19-fold increase in the peak pressure over the uncorrected case). Additionally, combining the phase correction method with rapid HAS simulations allows evaluation of such treatment metrics as the effect of misregistration on resulting pressure levels.

Conclusions

The method presented here is able to rapidly compute phases required to improve ultrasound focusing through the skull at multiple treatment locations. Combining phase correction with rapid simulation techniques allows for evaluation of various treatment metrics such as the effect of steering on pressure levels. Since the method computes 3D pressure patterns, it may also be suitable for predicting off-focus hot spots during treatments—a primary concern for transcranial HIFU. Additionally, the plastic-skull method presented here may be a useful tool in evaluating the effectiveness of phase correction methods.

A fast boundary element method for the scattering analysis of high-intensity focused ultrasound

High-intensity focused ultrasound (HIFU) techniques are promising modalities for the non-invasive treatment of cancer. For HIFU therapies of, e.g., liver cancer, one of the main challenges is the accurate focusing of the acoustic field inside a ribcage. Computational methods can play an important role in the patient-specific planning of these transcostal HIFU treatments. This requires the accurate modeling of acoustic scattering at ribcages. The use of a boundary element method (BEM) is an effective approach for this purpose because only the boundaries of the ribs have to be discretized instead of the standard approach to model the entire volume around the ribcage. This paper combines fast algorithms that improve the efficiency of BEM specifically for the high-frequency range necessary for transcostal HIFU applications. That is, a Galerkin discretized Burton-Miller formulation is used in combination with preconditioning and matrix compression techniques. In particular, quick convergence is achieved with the operator preconditioner that has been designed with on-surface radiation conditions for the high-frequency approximation of the Neumann-to-Dirichlet map. Realistic computations of acoustic scattering at 1 MHz on a human ribcage model demonstrate the effectiveness of this dedicated BEM algorithm for HIFU scattering analysis.

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.

Ultrasound-guided therapeutic focused ultrasound: current status and future directions

This paper reviews ultrasound imaging methods for the guidance of therapeutic focused ultrasound (USgFUS), with emphasis on real-time preclinical methods. Guidance is interpreted in the broadest sense to include pretreatment planning, siting of the FUS focus, real-time monitoring of FUS-tissue interactions, and real-time control of exposure and damage assessment. The paper begins with an overview and brief historical background of the early methods used for monitoring FUS-tissue interactions. Current imaging methods are described, and discussed in terms of sensitivity and specificity of the localisation of the FUS effects in both therapeutic and sub-therapeutic modes. Thermal and non-thermal effects are considered. These include cavitation-enhanced heating, tissue water boiling and cavitation. Where appropriate, USgFUS methods are compared with similar methods implemented using other guidance modalities, e.g. magnetic resonance imaging. Conclusions are drawn regarding the clinical potential of the various guidance methods, and the feasibility and current status of real-time implementation.