title: On the use of Gegenbauer reconstructions for shock wave propagation modeling.
abstract: In therapeutic ultrasound, the presence of shock waves can be significant due to the use of high intensity beams, as well as due to shock formation during inertial cavitation. Although modeling of such strongly nonlinear waves can be carried out using spectral methods, such calculations are typically considered impractical, since accurate calculations often require hundreds or even thousands of harmonics to be considered, leading to prohibitive computational times. Instead, time-domain algorithms which generally utilize Godunov-type finite-difference schemes are commonly used. Although these time domain methods can accurately model steep shock wave fronts, unlike spectral methods they are inherently unsuitable for modeling realistic tissue dispersion relations. Motivated by the need for a more general model, the use of Gegenbauer reconstructions as a postprocess tool to resolve the band-limitations of the spectral methods are investigated. The present work focuses on eliminating the Gibbs phenomenon when representing a steep wave front using a limited number of harmonics. Both plane wave and axisymmetric 2D transducer problems will be presented to characterize the proposed method.
authors: Department of Radiology, Harvard Medical School, Brigham and Women's Hospital, Boston, Massachusetts 02115, USA. jingy@bwh.harvard.edu
authors: Yun
authors: Y
authors: Jing
authors: Greg T
authors: GT
authors: Clement
mesh_terms: D000465:Algorithms; D003198:Computer Simulation; D004867:Equipment Design; D005583:Fourier Analysis; D019220:High-Energy Shock Waves; D008962:Models, Theoretical; D009038:Motion; D009716:Numerical Analysis, Computer-Assisted; D011312:Pressure; D013997:Time Factors; D014159:Transducers; D014465:Ultrasonics / Q000295:instrumentation / Q000379:methods*