Publication Details

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

KLEPÁRNÍK Petr, ZEMČÍK Pavel, TREEBY Bradley E. and JAROŠ Jiří. On-the-Fly Calculation of Time-Averaged Acoustic Intensity in Time-Domain Ultrasound Simulations Using a k-Space Pseudospectral Method. IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 69, no. 10, 2022, pp. 2917-2929. ISSN 1525-8955.
Czech title
Průběžný výpočet průměrné akustické intenzity v čase v ultrazvukových simulacích v časové oblasti pomocí k-space pseudospektrální metody
Type
journal article
Language
english
Authors
Klepárník Petr, Ing., Ph.D. (DCGM FIT BUT)
Zemčík Pavel, prof. Dr. Ing. (DCGM FIT BUT)
Treeby Bradley E. (UCL)
Jaroš Jiří, doc. Ing., Ph.D. (DCSY FIT BUT)
Keywords

Average acoustic intensity, compression, high intensity focused ultrasound, k-Wave toolbox, ultrasound simulation

Abstract

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. Due to the time staggering between the particle velocity and the acoustic pressure used in such simulations, the average intensity calculation is typically executed offline after the acoustic simulation consuming both disk space and time (the data can spread over terabytes). Our new approach calculates the average intensity during the acoustic simulation using the output coefficients of a new compression method which enables resolving the time staggering on-the-fly with huge disk space savings. To reduce RAM requirements, the article also presents a new 40-bit method for encoding compression complex coefficients. Experimental numerical simulations with the proposed method have shown that disk space requirements are up to 99% lower. The simulation speed was not significantly affected by the approach and the compression error did not affect the prediction accuracy of the thermal dose. From the standpoint of supercomputers, the new approach is significantly more economical. Saving computing resources increases the chances of real use of acoustic simulations in practice. The method can be applied to signals of a similar character, e.g., for electromagnetic radio waves.

Published
2022
Pages
2917-2929
Journal
IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control, vol. 69, no. 10, ISSN 1525-8955
Publisher
Institute of Electrical and Electronics Engineers
DOI
UT WoS
000861442800021
EID Scopus
BibTeX
@ARTICLE{FITPUB12807,
   author = "Petr Klep\'{a}rn\'{i}k and Pavel Zem\v{c}\'{i}k and E. Bradley Treeby and Ji\v{r}\'{i} Jaro\v{s}",
   title = "On-the-Fly Calculation of Time-Averaged Acoustic Intensity in Time-Domain Ultrasound Simulations Using a k-Space Pseudospectral Method",
   pages = "2917--2929",
   journal = "IEEE Transactions on Ultrasonics, Ferroelectrics, and Frequency Control",
   volume = 69,
   number = 10,
   year = 2022,
   ISSN = "1525-8955",
   doi = "10.1109/TUFFC.2022.3199173",
   language = "english",
   url = "https://www.fit.vut.cz/research/publication/12807"
}
Back to top