期刊
GEOPHYSICAL JOURNAL INTERNATIONAL
卷 179, 期 1, 页码 459-474出版社
OXFORD UNIV PRESS
DOI: 10.1111/j.1365-246X.2009.04305.x
关键词
Numerical solutions; Computational seismology; Wave propagation
资金
- China Scholarship Council
- National '863' Program of China [2007AA06Z218]
P>We derive explicit and new implicit staggered-grid finite-difference (FD) formulas for derivatives of first order with any order of accuracy by a plane wave theory and Taylor's series expansion. Furthermore, we arrive at a practical algorithm such that the tridiagonal matrix equations are formed by the implicit FD formulas derived from the fractional expansion of derivatives. Our results demonstrate that the accuracy of a (2N + 2)th-order implicit formula is nearly equivalent to or greater than that of a (4N)th-order explicit formula. The new implicit method only involves solving tridiagonal matrix equations. We also demonstrate that a (2N + 2)th-order implicit formulation requires nearly the same amount of memory and computation as those of a (2N + 4)th-order explicit formulation but attains the accuracy achieved by a (4N)th-order explicit formulation when additional cost of visiting arrays is not considered. Our analysis of efficiency and numerical modelling results for elastic wave propagation demonstrates that a high-order explicit staggered-grid method can be replaced by an implicit staggered-grid method of some order, which will increase the accuracy but not the computational cost.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据