An efficient parallel algorithm for solving viscoelastic wave equations

The high computational complexity of time-fractional viscoelastic wave equations limits their practical applications in seismic exploration and medical imaging. To address this challenge, this paper proposes a parallelized algorithm to address the practical application demands of viscoelastic wave e...

Celý popis

Uložené v:
Podrobná bibliografia
Vydané v:Mathematics and computers in simulation Ročník 241; s. 71 - 96
Hlavní autori: Li, Yaomeng, Guo, Xu
Médium: Journal Article
Jazyk:English
Vydavateľské údaje: Elsevier B.V 01.03.2026
Predmet:
ISSN:0378-4754
On-line prístup:Získať plný text
Tagy: Pridať tag
Žiadne tagy, Buďte prvý, kto otaguje tento záznam!
Popis
Shrnutí:The high computational complexity of time-fractional viscoelastic wave equations limits their practical applications in seismic exploration and medical imaging. To address this challenge, this paper proposes a parallelized algorithm to address the practical application demands of viscoelastic wave equations. The method leverages the rapid decay property of cubic B-spline local wavelets in their dual space and designs perfectly matched boundary conditions (PMBCs) to achieve closure of local spline coefficients. This approach enables accurate global spline reconstruction with only localized communication between adjacent patches. Furthermore, by integrating the distributed-parallel local spline simulator (DPLS) with a short-memory operator splitting (SMOS) scheme, we develop an efficient solver for viscoelastic wave equations. Numerical examples in one-dimensional (1D), two-dimensional (2D), and three-dimensional (3D) viscoelastic wave propagation scenarios validate the convergence, accuracy, and computational efficiency of the proposed method. •We propose a distributed-parallel spline simulator (DPLS) for solving the viscoelastic wave equation.•Stability conditions for SMOS-DPLS are validated by comparing with 1D analytic solutions.•SMOS-DPLS achieves fourth-order accuracy with perfectly matched boundary conditions.•SMOS-DPLS shows a speedup ratio over 50 on a 2D homogeneous medium with 15003 grid points.•The numerical solutions are implemented for 3D large-scale and variable-Q overthrust models.
ISSN:0378-4754
DOI:10.1016/j.matcom.2025.10.007