HPC Modeling of Coupled Elastic-Acoustic Wave Propagation in Biological Media: Numerical Validation
Abstract
Accurate numerical models of sound propagation through biological media are an important tool for many applications, from medical physics to studying the auditory system of humans or other animals. We model high-frequency elastic and acoustic wave propagation through the head anatomy of a common bottlenose dolphin (Tursiops truncatus), by means of the open-source software packages \textsc{specfem3d}, based on the spectral-element method, and \textsc{{\it k}-wave}, based on the pseudospectral method. To achieve sufficiently high performance, we ported the latter solver to C++, for Multi-GPU CUDA support via the slab-decomposition of three-dimensional Fast Fourier Transform approach. Because the two schemes differ fundamentally in how the propagation medium is discretized and internal (in particular fluid-solid) interfaces are treated, similarity between modeled signals across methods is a legitimate measure of model accuracy. Plane waves, depending on time like four-cycle sinusoidal bursts of varying central frequency (20--100 kHz), are numerically propagated through a computed-tomography-based anatomy model using both solvers. The sound is ``recorded'' in front of the rostrum and at the right inner-ear locations for comparison. We successfully cross-validate both solvers on high-fidelity High-Performance Computing (HPC) clusters. We find that the stability of results from both solvers grows as the corresponding spatial resolution is refined. At their highest resolutions, the two methods show excellent agreement, with normalized correlation exceeding 0.99 across the entire 20--100 kHz frequency range. Together, our results provide a validated HPC-simulation framework for wave propagation in biological media, with broader implications, e.g., for biosonar research, auditory biomechanics, and medical ultrasound.
Related
- LES-derived turbulence tank optical propagation with numerical adaptive-optics modeling: quantifying OAM modal crosstalk and mitigation
- Light pulse propagation through a dielectric slab. II. Two roads to the same field: numerical cross-validation of local (FDTD-ADE) and non-local (retarded integral-equation) solutions
- Numerical Investigation of Acoustic Streaming-Induced Shear Stress in 3D-Printed Bone-Mimicking Scaffolds
- Numerical modelling of the poto-poto by coupling of finite elements and boundary elements for acoustic characterization
- Boundary-induced medium mapping enables air-equivalent acoustic propagation and perfect absorption in water
BibTeX
@misc{ali2026modeling,
title = {{HPC Modeling of Coupled Elastic-Acoustic Wave Propagation in Biological Media: Numerical Validation}},
author = {Fawad Ali and Carlos García and Lapo Boschi},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2609.03644},
eprint = {2609.03644},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2609.03644},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - HPC Modeling of Coupled Elastic-Acoustic Wave Propagation in Biological Media: Numerical Validation
AU - Fawad Ali
AU - Carlos García
AU - Lapo Boschi
PY - 2026
DA - 2026/09/03
PB - arXiv
DO - 10.48550/arxiv.2609.03644
UR - https://arxiv.org/abs/2609.03644
AB - Accurate numerical models of sound propagation through biological media are an important tool for many applications, from medical physics to studying the auditory system of humans or other animals. We model high-frequency elastic and acoustic wave propagation through the head anatomy of a common bottlenose dolphin (Tursiops truncatus), by means of the open-source software packages \textsc{specfem3d}, based on the spectral-element method, and \textsc{{\it k}-wave}, based on the pseudospectral method. To achieve sufficiently high performance, we ported the latter solver to C++, for Multi-GPU CUDA support via the slab-decomposition of three-dimensional Fast Fourier Transform approach. Because the two schemes differ fundamentally in how the propagation medium is discretized and internal (in particular fluid-solid) interfaces are treated, similarity between modeled signals across methods is a legitimate measure of model accuracy. Plane waves, depending on time like four-cycle sinusoidal bursts of varying central frequency (20--100 kHz), are numerically propagated through a computed-tomography-based anatomy model using both solvers. The sound is ``recorded'' in front of the rostrum and at the right inner-ear locations for comparison. We successfully cross-validate both solvers on high-fidelity High-Performance Computing (HPC) clusters. We find that the stability of results from both solvers grows as the corresponding spatial resolution is refined. At their highest resolutions, the two methods show excellent agreement, with normalized correlation exceeding 0.99 across the entire 20--100 kHz frequency range. Together, our results provide a validated HPC-simulation framework for wave propagation in biological media, with broader implications, e.g., for biosonar research, auditory biomechanics, and medical ultrasound.
ER -