High-performance modeling acoustic and elastic waves using the parallel Dichotomy Algorithm

JOURNAL OF COMPUTATIONAL PHYSICS

卷:230

Article

Fatyanov, Alexey G.¹  Terekhov, Andrew V.^1,2,3 

[1] Inst Computat Math & Math Geophys, Novosibirsk 630090, Russia

[2] Budker Inst Nucl Phys, Novosibirsk 630090, Russia

[3] Novosibirsk State Univ, Novosibirsk 630090, Russia

关键词: Acoustic waves;    Elastic waves;    Tridiagonal matrix algorithm (TDMA);    Parallel Thomas algorithm;    Parallel Dichotomy Algorithm;    Laguerre transform;   

DOI  :  10.1016/j.jcp.2010.11.046

来源: Elsevier

【 摘 要 】

A high-performance parallel algorithm is proposed for modeling the propagation of acoustic and elastic waves in inhomogeneous media. An initial boundary-value problem is replaced by a series of boundary-value problems for a constant elliptic operator and different right-hand sides via the integral Laguerre transform. It is proposed to solve difference equations by the conjugate gradient method for acoustic equations and by the GMRES(k) method for modeling elastic waves. A preconditioning operator was the Laplace operator that is inverted using the variable separation method. The novelty of the proposed algorithm is using the Dichotomy Algorithm 1261, which was designed for solving a series of tridiagonal systems of linear equations, in the context of the preconditioning operator inversion. Via considering analytical solutions, it is shown that modeling wave processes for long instants of time requires high-resolution meshes. The proposed parallel fine-mesh algorithm enabled to solve real application seismic problems in acceptable time and with high accuracy. By solving model problems, it is demonstrated that the considered parallel algorithm possesses high performance and efficiency over a wide range of the number of processors (from 2 to 8192). (C) 2010 Elsevier Inc. All rights reserved.

【 授权许可】

Free