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 Citation Map 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 HighWire Citing Articles via Web of Science (2) Citing Articles via Google Scholar Google Scholar Articles by Kang, T.-S. Articles by Baag, C.-E. Search for Related Content GeoRef GeoRef Citation

Finite-Difference Seismic Simulation Combining Discontinuous Grids with Locally Variable Timesteps

School of Earth and Environmental Sciences
Seoul National University
San 56-1, Sillim-dong, Gwanak-gu
Seoul 151-747, Korea
tskang{at}seismic.snu.ac.kr, baagce{at}snu.ac.kr

A locally variable timestep scheme that matches with discontinuous grids in the finite-difference method is developed for the efficient simulation of seismic-wave propagation. The first-order velocity-stress formulations are used to obtain the spatial derivatives using finite-difference operators on a staggered grid. In the case of a media interface with high velocity contrast, the computational domain consists of two regions with different grid spacings. Each region roughly covers the medium of the lower or higher wave propagation velocity. There is a small overlap of the two regions, called the transitional zone, within the higher velocity medium. A grid three times coarser in the high-velocity region compared with the grid in the low-velocity region is used to avoid spatial oversampling. Temporal steps corresponding to the spatial sampling ratio between both regions are determined based on local stability criteria. The wave field in the margin of the region with the smaller timestep is linearly interpolated in time using the values calculated in the region with the larger one within the transitional zone. Since the temporal interpolation is a 1D operation performed separately from the spatial interpolation strategy employed to connect two regions with different grid spacings, the proposed scheme is not restricted to 2D or 3D problems with a specific order of accuracy of the spatial finite-difference approximation. The use of the locally variable timestep scheme with discontinuous grids results in remarkable savings of computation time and reductions in memory requirements, with the efficiency depending on the simulation model.

This article has been cited by other articles:

				Z.-L. Pei, L.-Y. Fu, G.-X. Yu, and L.-X. Zhang A Wavelet-Optimized Adaptive Grid Method for Finite-Difference Simulation of Wave Propagation Bulletin of the Seismological Society of America, February 1, 2009; 99(1): 302 - 313. [Abstract] [Full Text] [PDF]

				An Efficient Finite-Difference Method for Simulating 3D Seismic Response of Localized Basin Structures Bulletin of the Seismological Society of America, October 1, 2004; 94(5): 1690 - 1705.