Application of the difference Gaussian rules to solution of hyperbolic problems

被引:23
作者
Asvadurov, S
Druskin, V
Knizhnerman, L
机构
[1] Schlumberger Doll Res Ctr, Ridgefield, CT 06877 USA
[2] Cent Geophys Expedit, Moscow 123298, Russia
关键词
finite differences; optimal grids; exponential convergence; hyperbolic problems; linear elasticity; wave propagation;
D O I
10.1006/jcph.1999.6410
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Two of the authors earlier suggested a method of calculating special grid steps for three point finite-difference schemes which yielded exponential superconvergence of the Neumann-to-Dirichlet map. We apply this approach to solve the two-dimensional time-domain wave problem and the 2.5-D elasticity system in cylindrical coordinates. Our numerical experiments exhibit exponential convergence at prescribed points, with the cost per grid node close to that of the standard second order finite-difference scheme. The scheme demonstrates high accuracy with slightly more than two grid points per wavelength. The reduction of the grid size by one order compared to the standard scheme with the equidistant grids is observed. (C) 2000 Academic Press.
引用
收藏
页码:116 / 135
页数:20
相关论文
共 4 条
[1]  
DRUSKIN V, IN PRESS SIAM J NUME
[2]  
DRUSKIN V, 1997, EMB0029722 SCHLUMB D
[3]  
Fornberg B., 1996, A Practical Guide to Pseudospectral Methods
[4]  
PARLETT B. N., 1980, The Symmetric Eigenvalue Problem