View in EDS

Iterative Born solver for the acoustic Helmholtz equation with heterogeneous sound speed and density.

Saved in:
Bibliographic Details
Title: Iterative Born solver for the acoustic Helmholtz equation with heterogeneous sound speed and density.
Authors: Stanziola, Antonio1 (AUTHOR) stanziola.antonio@gmail.com, Arridge, Simon R.2 (AUTHOR), Treeby, Bradley E.1 (AUTHOR), Cox, Benjamin T.1 (AUTHOR) b.cox@ucl.ac.uk
Source: Journal of the Acoustical Society of America. Feb2026, Vol. 159 Issue 2, p1457-1470. 14p.
Subjects: Helmholtz equation, Born approximation, Iterative methods (Mathematics), Density, Ultrasonic imaging, Inhomogeneous materials, Speed of sound, Fast Fourier transforms
Abstract: Efficient numerical solution of the acoustic Helmholtz equation in heterogeneous media remains challenging, particularly for large-scale problems with spatially varying density—a limitation that restricts applications in biomedical acoustics and seismic imaging. A fast iterative solver that extends the convergent Born series [Osnabrugge, Leedumrongwatthanakun, and Vellekoop, J. Comput. Phys. 322, 113–124 (2016)] method to handle arbitrary variations in sound speed, density, and absorption simultaneously is presented. This approach reformulates the Helmholtz equation as a first-order system and applies the universal split-preconditioner from Vettenburg and Vellekoop [arXiv:2207.14222v2 (2022)], yielding a matrix-free algorithm that leverages Fast Fourier Transforms for computational efficiency. Unlike existing Born series methods, this solver accommodates heterogeneous density without requiring expensive matrix decompositions or preprocessing steps, making it suitable for large-scale three-dimensional problems with minimal memory overhead. The method provides forward and adjoint solutions, enabling its application for inverse problems. Accuracy is validated through comparison against an analytical solution and the solver's practical utility is demonstrated through transcranial ultrasound simulations. The solver achieves convergence for strong scattering scenarios, offering a computationally efficient alternative to time-domain methods and matrix-based Helmholtz solvers for applications ranging from medical ultrasound treatment planning to seismic exploration. [ABSTRACT FROM AUTHOR]
Copyright of Journal of the Acoustical Society of America is the property of American Institute of Physics and its content may not be copied or emailed to multiple sites without the copyright holder's express written permission. Additionally, content may not be used with any artificial intelligence tools or machine learning technologies. However, users may print, download, or email articles for individual use. This abstract may be abridged. No warranty is given about the accuracy of the copy. Users should refer to the original published version of the material for the full abstract. (Copyright applies to all Abstracts.)
Database: Engineering Source
Description
Abstract:Efficient numerical solution of the acoustic Helmholtz equation in heterogeneous media remains challenging, particularly for large-scale problems with spatially varying density—a limitation that restricts applications in biomedical acoustics and seismic imaging. A fast iterative solver that extends the convergent Born series [Osnabrugge, Leedumrongwatthanakun, and Vellekoop, J. Comput. Phys. 322, 113–124 (2016)] method to handle arbitrary variations in sound speed, density, and absorption simultaneously is presented. This approach reformulates the Helmholtz equation as a first-order system and applies the universal split-preconditioner from Vettenburg and Vellekoop [arXiv:2207.14222v2 (2022)], yielding a matrix-free algorithm that leverages Fast Fourier Transforms for computational efficiency. Unlike existing Born series methods, this solver accommodates heterogeneous density without requiring expensive matrix decompositions or preprocessing steps, making it suitable for large-scale three-dimensional problems with minimal memory overhead. The method provides forward and adjoint solutions, enabling its application for inverse problems. Accuracy is validated through comparison against an analytical solution and the solver's practical utility is demonstrated through transcranial ultrasound simulations. The solver achieves convergence for strong scattering scenarios, offering a computationally efficient alternative to time-domain methods and matrix-based Helmholtz solvers for applications ranging from medical ultrasound treatment planning to seismic exploration. [ABSTRACT FROM AUTHOR]
ISSN:00014966
DOI:10.1121/10.0042259