Quick Search:    advanced search

GSW Home   GeoRef Home   My GSW Alerts   Contact GSW   About GSW   Journals List   Help

This Article Figures Only Full Text Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in Web of Science Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Web of Science (3) Citing Articles via Google Scholar Google Scholar Articles by Zeng, Y. Q. Articles by Liu, Q. H. Search for Related Content GeoRef GeoRef Citation

A Multidomain PSTD Method for 3D Elastic Wave Equations

Department of Electrical and Computer Engineering
Duke University
Durham, North Carolina 27708

A 3D multidomain pseudospectral time-domain method is developed for elastic wave equations. The method is based on the spectral derivative operator approximated by Chebyshev or Lagrange polynomials. Unlike the Fourier method that assumes periodic boundary conditions, the Chebyshev pseudospectral method allows for the incorporation of various boundary conditions (such as the free surface boundary condition) into the numerical scheme. In this multidomain scheme, the computational domain is decomposed into a set of subdomains conformal to the problem geometry. Each curved subdomain is then mapped onto a cube in the curvilinear coordinates so that a tensor-product Chebyshev grid can be utilized without the staircasing error. An unsplit-field perfectly matched layer is developed as the absorbing boundary condition. Numerical examples show that this scheme is efficient for simulating elastic waves phenomena in the presence of complex objects. The method is found to be significantly more efficient than the finite-difference time-domain method in terms of memory and run-time requirements.