Iterative Born Solver for the Acoustic Helmholtz Equation with Heterogeneous Sound Speed and Density
Fuente:
arXiv
Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866916855576264704 |
|---|---|
| author | Stanziola, Antonio Arridge, Simon R. Treeby, Bradley E. Cox, Benjamin T. |
| author_facet | Stanziola, Antonio Arridge, Simon R. Treeby, Bradley E. Cox, Benjamin T. |
| contents | 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. We present a fast iterative solver that extends the Convergent Born Series method to handle arbitrary variations in sound speed, density, and absorption simultaneously. Our approach reformulates the Helmholtz equation as a first-order system and applies Vettenburg and Vellekoop's universal split-preconditioner, yielding a matrix-free algorithm that leverages Fast Fourier Transforms for computational efficiency. Unlike existing Born series methods, our solver accommodates heterogeneous density without requiring expensive matrix decompositions or pre-processing steps, making it suitable for large-scale 3D problems with minimal memory overhead. The method provides both forward and adjoint solutions, enabling its application for inverse problems. We validate accuracy through comparison against an analytical solution and demonstrate the solver's practical utility 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2507_16087 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Iterative Born Solver for the Acoustic Helmholtz Equation with Heterogeneous Sound Speed and Density Stanziola, Antonio Arridge, Simon R. Treeby, Bradley E. Cox, Benjamin T. Computational Physics 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. We present a fast iterative solver that extends the Convergent Born Series method to handle arbitrary variations in sound speed, density, and absorption simultaneously. Our approach reformulates the Helmholtz equation as a first-order system and applies Vettenburg and Vellekoop's universal split-preconditioner, yielding a matrix-free algorithm that leverages Fast Fourier Transforms for computational efficiency. Unlike existing Born series methods, our solver accommodates heterogeneous density without requiring expensive matrix decompositions or pre-processing steps, making it suitable for large-scale 3D problems with minimal memory overhead. The method provides both forward and adjoint solutions, enabling its application for inverse problems. We validate accuracy through comparison against an analytical solution and demonstrate the solver's practical utility 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. |
| title | Iterative Born Solver for the Acoustic Helmholtz Equation with Heterogeneous Sound Speed and Density |
| topic | Computational Physics |
| url | https://arxiv.org/abs/2507.16087 |