Running Large-Scale Ultrasound Simulations on 256 Salomon Nodes

Local domain decomposition, k-Wave toolbox, GPU

Ultrasound simulation is a critical component of model-based treatment planning for ultrasound therapy. However, the domain sizes are typically thousands of wavelengths in size, leading to large-scale numerical models with 10s of billions of unknowns. This paper presents a novel local Fourier basis domain decomposition for full-wave ultrasound propagation simulations with a custom bell function which ensures that the numerical error is below 0.1% while enabling almost ideal strong scaling.

Realistic benchmarks with 6144 Haswell cores across 256 Salomon nodes show significant improvement over native 3D FFT decomposition.

• pdf poster.pdf 1801 KB