1
|
Xu P, Wu H, Shen G. Characterization of weakly nonlinear effects in relationship to transducer parameters in focused ultrasound therapy. Med Phys 2024. [PMID: 38935266 DOI: 10.1002/mp.17270] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2023] [Revised: 06/04/2024] [Accepted: 06/13/2024] [Indexed: 06/28/2024] Open
Abstract
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.
Collapse
Affiliation(s)
- Peng Xu
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, People's Republic of China
| | - Hao Wu
- Shanghai Shende Green Medical Era Healthcare Technology Co., Ltd, Shanghai, People's Republic of China
| | - Guofeng Shen
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, People's Republic of China
| |
Collapse
|
2
|
Kagami S, Kanagawa T. Weakly nonlinear focused ultrasound in viscoelastic media containing multiple bubbles. ULTRASONICS SONOCHEMISTRY 2023; 97:106455. [PMID: 37271029 DOI: 10.1016/j.ultsonch.2023.106455] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/19/2023] [Revised: 05/12/2023] [Accepted: 05/23/2023] [Indexed: 06/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Shunsuke Kagami
- Department of Engineering Mechanics and Energy, Degree Program of Systems and Information Engineering, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8573, Japan.
| | - Tetsuya Kanagawa
- Department of Engineering Mechanics and Energy, Degree Program of Systems and Information Engineering, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8573, Japan.
| |
Collapse
|
3
|
Kagami S, Kanagawa T. Weakly nonlinear propagation of focused ultrasound in bubbly liquids with a thermal effect: Derivation of two cases of Khokolov-Zabolotskaya-Kuznetsoz equations. ULTRASONICS SONOCHEMISTRY 2022; 88:105911. [PMID: 35810619 PMCID: PMC9696949 DOI: 10.1016/j.ultsonch.2022.105911] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/11/2021] [Revised: 01/04/2022] [Accepted: 01/05/2022] [Indexed: 05/06/2023]
Abstract
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.
Collapse
Affiliation(s)
- Shunsuke Kagami
- Department of Engineering Mechanics and Energy, Graduate School of Systems and Information Engineering, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8573, Japan
| | - Tetsuya Kanagawa
- Department of Engineering Mechanics and Energy, Faculty of Engineering, Information and Systems, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8573, Japan.
| |
Collapse
|
4
|
Yuldashev PV, Karzova MM, Kreider W, Rosnitskiy PB, Sapozhnikov OA, Khokhlova VA. "HIFU Beam:" A Simulator for Predicting Axially Symmetric Nonlinear Acoustic Fields Generated by Focused Transducers in a Layered Medium. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2021; 68:2837-2852. [PMID: 33877971 PMCID: PMC8486313 DOI: 10.1109/tuffc.2021.3074611] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
"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.
Collapse
|
5
|
Martin E, Jaros J, Treeby BE. Experimental Validation of k-Wave: Nonlinear Wave Propagation in Layered, Absorbing Fluid Media. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2020; 67:81-91. [PMID: 31535990 DOI: 10.1109/tuffc.2019.2941795] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
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.
Collapse
|
6
|
Numerical Simulations of the Nonlinear Interaction of a Bubble Cloud and a High Intensity Focused Ultrasound Field. ACOUSTICS 2019. [DOI: 10.3390/acoustics1040049] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
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.
Collapse
|
7
|
Gu J, Jing Y. Simulation of the Second-Harmonic Ultrasound Field in Heterogeneous Soft Tissue Using a Mixed-Domain Method. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2019; 66:669-675. [PMID: 30640608 PMCID: PMC6492553 DOI: 10.1109/tuffc.2019.2892753] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/30/2023]
Abstract
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.
Collapse
|
8
|
Ted Christopher P, Parker KJ. The nonlinear ultrasound needle pulse. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2018; 144:861. [PMID: 30180703 DOI: 10.1121/1.5050519] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/08/2018] [Accepted: 07/26/2018] [Indexed: 06/08/2023]
Abstract
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.
Collapse
Affiliation(s)
| | - Kevin J Parker
- Departments of Electrical & Computer and of Biomedical Engineering, University of Rochester, Rochester, New York 14627, USA
| |
Collapse
|
9
|
Haddadi S, Ahmadian MT. Numerical and Experimental Evaluation of High-Intensity Focused Ultrasound-Induced Lesions in Liver Tissue Ex Vivo. JOURNAL OF ULTRASOUND IN MEDICINE : OFFICIAL JOURNAL OF THE AMERICAN INSTITUTE OF ULTRASOUND IN MEDICINE 2018; 37:1481-1491. [PMID: 29193279 DOI: 10.1002/jum.14491] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/01/2017] [Revised: 08/13/2017] [Accepted: 09/03/2017] [Indexed: 06/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Samaneh Haddadi
- School of Mechanical Engineering, Sharif University of Technology, Tehran, Iran
| | - Mohammad Taghi Ahmadian
- Center of Excellence in Design, Robotics, and Automation, School of Mechanical Engineering, Sharif University of Technology, Tehran, Iran
| |
Collapse
|
10
|
Qiao S, Jackson E, Coussios CC, Cleveland RO. Simulation of nonlinear propagation of biomedical ultrasound using pzflex and the Khokhlov-Zabolotskaya-Kuznetsov Texas code. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2016; 140:2039. [PMID: 27914432 PMCID: PMC5849034 DOI: 10.1121/1.4962555] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/05/2016] [Revised: 08/12/2016] [Accepted: 08/26/2016] [Indexed: 05/22/2023]
Abstract
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.
Collapse
Affiliation(s)
- Shan Qiao
- Institute of Biomedical Engineering, Department of Engineering, University of Oxford, Old Road Campus Research Building, Oxford, OX3 7DQ, United Kingdom
| | - Edward Jackson
- Institute of Biomedical Engineering, Department of Engineering, University of Oxford, Old Road Campus Research Building, Oxford, OX3 7DQ, United Kingdom
| | - Constantin C Coussios
- Institute of Biomedical Engineering, Department of Engineering, University of Oxford, Old Road Campus Research Building, Oxford, OX3 7DQ, United Kingdom
| | - Robin O Cleveland
- Institute of Biomedical Engineering, Department of Engineering, University of Oxford, Old Road Campus Research Building, Oxford, OX3 7DQ, United Kingdom
| |
Collapse
|
11
|
Du Y, Fan R, Li Y, Chen S, Jensen JA. A framework for simulating ultrasound imaging based on first order nonlinear pressure-velocity relations. ULTRASONICS 2016; 69:152-165. [PMID: 27107165 DOI: 10.1016/j.ultras.2016.03.015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/25/2015] [Revised: 03/27/2016] [Accepted: 03/27/2016] [Indexed: 06/05/2023]
Abstract
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%.
Collapse
Affiliation(s)
- Yigang Du
- National-Regional Key Technology Engineering Laboratory for Medical Ultrasound, Guangdong Key Laboratory for Biomedical Measurements and Ultrasound Imaging, Department of Biomedical Engineering, School of Medicine, Shenzhen University, CN-518060 Shenzhen, China; Shenzhen Mindray Bio-Medical Electronics Co. Ltd., CN-518057 Shenzhen, China; College of Information Engineering, Shenzhen University, CN-518060 Shenzhen, China.
| | - Rui Fan
- Shenzhen Mindray Bio-Medical Electronics Co. Ltd., CN-518057 Shenzhen, China
| | - Yong Li
- Shenzhen Mindray Bio-Medical Electronics Co. Ltd., CN-518057 Shenzhen, China
| | - Siping Chen
- National-Regional Key Technology Engineering Laboratory for Medical Ultrasound, Guangdong Key Laboratory for Biomedical Measurements and Ultrasound Imaging, Department of Biomedical Engineering, School of Medicine, Shenzhen University, CN-518060 Shenzhen, China
| | - Jørgen Arendt Jensen
- Center for Fast Ultrasound Imaging, Department of Electrical Engineering, Technical University of Denmark, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
12
|
Gu J, Jing Y. Modeling of wave propagation for medical ultrasound: a review. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2015; 62:1979-1993. [PMID: 26559627 DOI: 10.1109/tuffc.2015.007034] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
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.
Collapse
|
13
|
Solovchuk MA, Hwang SC, Chang H, Thiriet M, Sheu TWH. Temperature elevation by HIFU inex vivoporcine muscle: MRI measurement and simulation study. Med Phys 2014; 41:052903. [DOI: 10.1118/1.4870965] [Citation(s) in RCA: 29] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023] Open
|
14
|
Jing Y, Tao M, Cannata J. An improved wave-vector frequency-domain method for nonlinear wave modeling. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2014; 61:515-524. [PMID: 24569255 DOI: 10.1109/tuffc.2014.2935] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
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.
Collapse
|
15
|
Solovchuk M, Sheu TWH, Thiriet M. Simulation of nonlinear Westervelt equation for the investigation of acoustic streaming and nonlinear propagation effects. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2013; 134:3931-3942. [PMID: 24180802 DOI: 10.1121/1.4821201] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Maxim Solovchuk
- Center of Advanced Study in Theoretical Sciences (CASTS), National Taiwan University, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan
| | | | | |
Collapse
|
16
|
Hasani MH, Gharibzadeh S, Farjami Y, Tavakkoli J. Unmitigated numerical solution to the diffraction term in the parabolic nonlinear ultrasound wave equation. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2013; 134:1775-1790. [PMID: 23967912 DOI: 10.1121/1.4774278] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
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.
Collapse
Affiliation(s)
- Mojtaba H Hasani
- Department of Biomedical Engineering, Amirkabir University of Technology, Tehran, Iran
| | | | | | | |
Collapse
|
17
|
Albin N, Bruno OP, Cheung TY, Cleveland RO. Fourier continuation methods for high-fidelity simulation of nonlinear acoustic beams. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2012; 132:2371-2387. [PMID: 23039433 DOI: 10.1121/1.4742722] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/01/2023]
Abstract
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.
Collapse
Affiliation(s)
- Nathan Albin
- Mathematics, Kansas State University, 138 Cardwell Hall, Manhattan, Kansas 66506, USA
| | | | | | | |
Collapse
|
18
|
Jing Y, Wang T, Clement GT. A k-space method for moderately nonlinear wave propagation. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2012; 59:1664-73. [PMID: 22899114 PMCID: PMC3777432 DOI: 10.1109/tuffc.2012.2372] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
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.
Collapse
Affiliation(s)
- Yun Jing
- Department of Mechanical and Aerospace Engineering, North Carolina State University, Raleigh, NC, USA.
| | | | | |
Collapse
|
19
|
Treeby BE, Jaros J, Rendell AP, Cox BT. Modeling nonlinear ultrasound propagation in heterogeneous media with power law absorption using a k-space pseudospectral method. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2012; 131:4324-36. [PMID: 22712907 DOI: 10.1121/1.4712021] [Citation(s) in RCA: 177] [Impact Index Per Article: 14.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/07/2023]
Abstract
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.
Collapse
Affiliation(s)
- Bradley E Treeby
- Research School of Engineering, College of Engineering and Computer Science, The Australian National University, Canberra ACT 0200, Australia.
| | | | | | | |
Collapse
|
20
|
Dagrau F, Rénier M, Marchiano R, Coulouvrat F. Acoustic shock wave propagation in a heterogeneous medium: a numerical simulation beyond the parabolic approximation. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2011; 130:20-32. [PMID: 21786874 DOI: 10.1121/1.3583549] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/31/2023]
Abstract
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.
Collapse
Affiliation(s)
- Franck Dagrau
- Université Pierre et Marie Curie - Paris 6, Institut Jean le Rond d'Alembert (UMR UPMC/CNRS 7190), 4 place Jussieu, 75252 Paris Cedex 05, France
| | | | | | | |
Collapse
|
21
|
Yuldashev PV, Khokhlova VA. SIMULATION OF THREE-DIMENSIONAL NONLINEAR FIELDS OF ULTRASOUND THERAPEUTIC ARRAYS. ACOUSTICAL PHYSICS 2011; 57:334-343. [PMID: 21804751 PMCID: PMC3145364 DOI: 10.1134/s1063771011030213] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/23/2023]
Abstract
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.
Collapse
|
22
|
Averiyanov M, Blanc-Benon P, Cleveland RO, Khokhlova V. Nonlinear and diffraction effects in propagation of N-waves in randomly inhomogeneous moving media. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2011; 129:1760-72. [PMID: 21476633 DOI: 10.1121/1.3557034] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/08/2023]
Abstract
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.
Collapse
|
23
|
Shepherd MR, Gee KL, Wochner MS. Short-range shock formation and coalescence in numerical simulation of broadband noise propagation. THE JOURNAL OF THE ACOUSTICAL SOCIETY OF AMERICA 2009; 126:2886-2893. [PMID: 20000901 DOI: 10.1121/1.3243466] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
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.
Collapse
Affiliation(s)
- Micah R Shepherd
- Applied Research Laboratory, The Pennsylvania State University, P.O. Box 30, State College, Pennsylvania 16804, USA.
| | | | | |
Collapse
|
24
|
Näsholm SP, Hansen R, Måsøy SE, Johansen TF, Angelsen BAJ. Transmit beams adapted to reverberation noise suppression using dual-frequency SURF imaging. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2009; 56:2124-2133. [PMID: 19942500 DOI: 10.1109/tuffc.2009.1295] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/28/2023]
Abstract
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.
Collapse
Affiliation(s)
- Sven Peter Näsholm
- Dept. of Circulation & Imaging, Norwegian Univ. of Sci. & Technol., Trondheim, Norway.
| | | | | | | | | |
Collapse
|
25
|
Computational aspects in high intensity ultrasonic surgery planning. Comput Med Imaging Graph 2009; 34:69-78. [PMID: 19740625 DOI: 10.1016/j.compmedimag.2009.08.001] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2009] [Revised: 07/24/2009] [Accepted: 08/07/2009] [Indexed: 01/04/2023]
Abstract
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.
Collapse
|
26
|
Marchiano R, Coulouvrat F, Ganjehi L, Thomas JL. Numerical investigation of the properties of nonlinear acoustical vortices through weakly heterogeneous media. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2008; 77:016605. [PMID: 18351949 DOI: 10.1103/physreve.77.016605] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2007] [Indexed: 05/26/2023]
Abstract
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.
Collapse
Affiliation(s)
- Régis Marchiano
- Institut Jean le Rond d'Alembert (UMR CNRS 7190), Université Pierre et Marie Curie-Paris 6 4, place Jussieu 75252, Paris, Cedex 05, France.
| | | | | | | |
Collapse
|