A Pseudospectral Chebychev method for the 2D Wave Equation with Domain Stretching and Absorbing Boun

In this paper we develop a method for the simulation of wave propagation on artificially bounded domains. The acoustic wave equation is solved at all points away from the boundaries by a pseudospectral Chebychev method. Absorption at the boundaries is obtained by applying one-way wave equations at the boundaries, without the use of damping layers. The theoretical reflection coefficient for the method is compared to theoretical estimates of reflection coefficients for a Fourier model of the problem. These estimates are confirmed by numerical results. Modification of the method by a transformation of the grid to allow for better resolution at the center of the grid reduces the maximum eigenvalues of the differential operator. Consequently, for stability the maximum timestep isO(1/N) as compared toO(1/N2) for the standard Chebychev method. Therefore, the Chebychev method can be implemented with efficiency comparable to that of the Fourier method. Moreover, numerical results presented demonstrate the superior performance of the new method.

[1]  Robert L. Higdon,et al.  Numerical absorbing boundary conditions for the wave equation , 1987 .

[2]  D. Kosloff,et al.  Solution of the equations of dynamic elasticity by a Chebychev spectral method , 1990 .

[3]  Moshe Reshef,et al.  A nonreflecting boundary condition for discrete acoustic and elastic wave equations , 1985 .

[4]  Lloyd N. Trefethen,et al.  Lax-stability of fully discrete spectral methods via stability regions and pseudo-eigenvalues , 1990 .

[5]  Robert L. Higdon,et al.  Absorbing boundary conditions for di erence approxima-tions to the multidimensional wave equation , 1986 .

[6]  L. Trefethen,et al.  An Instability Phenomenon in Spectral Methods , 1987 .

[7]  T. Driscoll,et al.  Pseudospectra of the Wave Operator with an Absorbing Boundary , 1993 .

[8]  Ekkehart Tessmer,et al.  3-D Elastic Modeling With Surface Topography By a Chebychev Spectral Method , 1994 .

[9]  David A. Kopriva.  Spectral solution of acoustic boundary-value problems , 1992 .

[10]  Dan Kosloff,et al.  A modified Chebyshev pseudospectral method with an O(N –1 ) time step restriction , 1993 .

[11]  Rosemary A. Renaut,et al.  Abstract of papers to appear in future issuesAbsorbing boundary conditions, difference operators, and stability , 1992 .

[12]  K. W. Morton,et al.  Numerical methods for fluid dynamics , 1987 .

[13]  T. A. Zang,et al.  Spectral methods for fluid dynamics , 1987 .

[14]  A. Majda,et al.  Radiation boundary conditions for acoustic and elastic wave calculations , 1979 .

[15]  John L. Volakis,et al.  Absorbing boundary conditions , 1995 .

[16]  Edip Baysal,et al.  Forward modeling by a Fourier method , 1982 .

[17]  R. Kosloff,et al.  Absorbing boundaries for wave propagation problems , 1986 .

[18]  L. Trefethen,et al.  THE EIGENVALUES OF SECOND-ORDER SPECTRAL DIFFERENTIATION MATRICES* , 1988 .

[19]  L. Trefethen,et al.  Wide-angle one-way wave equations. , 1988, The Journal of the Acoustical Society of America.

[20]  J. Miller.  Numerical Analysis , 1966, Nature.

[21]  R. Higdon.  Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation , 1986 .

[22]  Eli Turkel,et al.  Mappings and accuracy for Chebyshev pseudo-spectral approximations , 1992 .

[23]  J. Carcione.  Time-dependent boundary conditions for the 2-D linear anisotropic-viscoelastic wave equation , 1994 .

[24]  N. A. Haskell.  The Dispersion of Surface Waves on Multilayered Media , 1953 .