| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
Département de Sismologie (URA 195) Institut de Physique du Globe de Paris, 4, Place Jussieu, 75252—Paris Cedex 05, France
Abstract
We present the spectral element method to simulate elastic-wave propagation in realistic geological structures involving complieated free-surface topography and material interfaces for two- and three-dimensional geometries. The spectral element method introduced here is a high-order variational method for the spatial approximation of elastic-wave equations. The mass matrix is diagonal by construction in this method, which drastically reduces the computational cost and allows an efficient parallel implementation. Absorbing boundary conditions are introduced in variational form to simulate unbounded physical domains. The time discretization is based on an energy-momentum conserving scheme that can be put into a classical explicit-implicit predictor/multi-corrector format. Long-term energy conservation and stability properties are illustrated as well as the efficiency of the absorbing conditions. The associated Courant condition behaves as 
tC < O (nel–1/ndN–2), with nel the number of elements, nd the spatial dimension, and N the polynomial order. In practice, a spatial sampling of approximately 5 points per wavelength is found to be very accurate when working with a polynomial degree of N = 8. The accuracy of the method is shown by comparing the spectral element solution to analytical solutions of the classical two-dimensional (2D) problems of Lamb and Garvin. The flexibility of the method is then illustrated by studying more realistic 2D models involving realistic geometries and complex free-boundary conditions. Very accurate modeling of Rayleigh-wave propagation, surface diffraction, and Rayleigh-to-body-wave mode conversion associated with the free-surface curvature are obtained at low computational cost. The method is shown to provide an efficient tool to study the diffraction of elastic waves by three-dimensional (3D) surface topographies and the associated local effects on strong ground motion. Complex amplification patterns, both in space and time, are shown to occur even for a gentle hill topography. Extension to a heterogeneous hill structure is considered. The efficient implementation on parallel distributed memory architectures will allow to perform real-time visualization and interactive physical investigations of 3D amplification phenomena for seismic risk assessment.
This article has been cited by other articles:
|
|
 |
|
  S.-J. Lee, H.-W. Chen, Q. Liu, D. Komatitsch, B.-S. Huang, and J. Tromp Three-Dimensional Simulations of Seismic-Wave Propagation in the Taipei Basin with Realistic Topography Based upon the Spectral-Element Method Bulletin of the Seismological Society of America, February 1, 2008; 98(1): 253 - 264. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
T. Ichimura, M. Hori, and H. Kuwamoto Earthquake Motion Simulation with Multiscale Finite-Element Analysis on Hybrid Grid Bulletin of the Seismological Society of America, August 1, 2007; 97(4): 1133 - 1143. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
Z. Ge and X. Chen Wave Propagation in Irregularly Layered Elastic Models: A Boundary Element Approach with a Global Reflection/Transmission Matrix Propagator Bulletin of the Seismological Society of America, June 1, 2007; 97(3): 1025 - 1031. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
T. Yang and Y. Shen Frequency-Dependent Crustal Correction for Finite-Frequency Seismic Tomography Bulletin of the Seismological Society of America, December 1, 2006; 96(6): 2441 - 2448. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Park and E. Kausel Response of Layered Half-Space Obtained Directly in the Time Domain, Part I: SH Sources Bulletin of the Seismological Society of America, October 1, 2006; 96(5): 1795 - 1809. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
F. J. Sanchez-Sesma and U. Iturraran-Viveros The Classic Garvin's Problem Revisited Bulletin of the Seismological Society of America, August 1, 2006; 96(4A): 1344 - 1351. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
L. Zhao, T. H. Jordan, K. B. Olsen, and P. Chen Frechet Kernels for Imaging Regional Earth Structure Based on Three-Dimensional Reference Models Bulletin of the Seismological Society of America, December 1, 2005; 95(6): 2066 - 2080. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
P. Dangla, J.-F. Semblat, H. Xiao, and N. Delepine A Simple and Efficient Regularization Method for 3D BEM: Application to Frequency-Domain Elastodynamics Bulletin of the Seismological Society of America, October 1, 2005; 95(5): 1916 - 1927. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J.-P. Ampuero and F. A. Dahlen Ambiguity of the Moment Tensor Bulletin of the Seismological Society of America, April 1, 2005; 95(2): 390 - 400. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Bielak Reply to "Comment on 'Domain Reduction Method for Three-Dimensional Earthquake Modeling in Localized Regions, Part I: Theory,' by J. Bielak, K. Loukakis, Y. Hisada, and C. Yoshimura, and 'Part II: Verification and Applications,' by C. Yoshimura, J. Bielak, Y. Hisada, and A. Fernandez," by E. Faccioli, M. Vanini, R. Paolucci, and M. Stupazzini Bulletin of the Seismological Society of America, April 1, 2005; 95(2): 770 - 773. [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.
|
|
|
|
|
|
Spectral-Element Moment Tensor Inversions for Earthquakes in Southern California Bulletin of the Seismological Society of America, October 1, 2004; 94(5): 1748 - 1761.
|
|
|
|
|
|
Hybrid Modeling of Elastic P-SV Wave Motion: A Combined Finite-Element and Staggered-Grid Finite-Difference Approach Bulletin of the Seismological Society of America, August 1, 2004; 94(4): 1557 - 1563.
|
|
|
|
|
|
D. Komatitsch, Q. Liu, J. Tromp, P. Suss, C. Stidham, and J. H. Shaw Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method Bulletin of the Seismological Society of America, February 1, 2004; 94(1): 187 - 206. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
A Comparison of Synthetic Seismograms for 2D Structures: Semianalytical versus Numerical Bulletin of the Seismological Society of America, December 1, 2003; 93(6): 2752 - 2757.
|
|
|
|
|
|
A Nearly Analytic Discrete Method for Acoustic and Elastic Wave Equations in Anisotropic Media Bulletin of the Seismological Society of America, April 1, 2003; 93(2): 882 - 890.
|
|
|
|
|
|
Scattering Attenuation of 2D Elastic Waves: Theory and Numerical Modeling Using a Wavelet-Based Method Bulletin of the Seismological Society of America, April 1, 2003; 93(2): 922 - 938.
|
|
|
|
 |
|
 D. Komatitsch, J. Ritsema, and J. Tromp The Spectral-Element Method, Beowulf Computing, and Global Seismology Science, November 29, 2002; 298(5599): 1737 - 1742. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
R. Paolucci Numerical evaluation of the effect of cross-coupling of different components of ground motion in site response analyses Bulletin of the Seismological Society of America, August 1, 1999; 89(4): 877 - 887. [Abstract] [PDF]
|
|
|
|
|
|
E. Faccioli and A. Quarteroni Comment on "The spectral element method: An efficient tool to simulate the seismic response of 2D and 3D geological structures," by D. Komatitsch and J.-P. Vilotte Bulletin of the Seismological Society of America, February 1, 1999; 89(1): 331 - 331. [PDF]
|
|
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
