THE SPECTRAL ELEMENT METHOD (SEM):FORMULATION AND IMPLEMENTATION OF THE METHOD FOR ENGINEERING SEISMOLOGY PROBLEMS
Various numerical methods have been proposed and used to investigate wave propagation in realistic earth media.Recently an innovative numerical method, known as the Spectral Element Method (SEM), has been developed and used in connection with wave propagation problems in 3D elastic media (Komatitsch and Tromp, 1999). The SEM combines the flexibility of a finite element method with the accuracy of a spectral method and it is not cumbersome in dealing with non-flat free surface and spatially variable anelastic attenuation. The SEM is a highly accurate numerical method that has its origins in computational fluid dynamics. One uses a weak formulation of the equations of motion, which are solved on a mesh of hexahedral elements that is adapted to the free surface and to the main internal discontinuities of the model. The wavefield on the elements is discretized using high-degree Lagrange interpolants, and integration over an element is accomplished based upon the Gauss-Lobatto-Legendre integration rule. This combination of discretization and integration results in a diagonal mass matrix, which greatly simplifies the algorithm. The most important property of the SEM is that the mass matrix is exactly diagonal by construction,which drastically simplifies the impemenation and reduces the computational cost because one can use an explicit time integration scheme without having to invert a linear system. Furthermore, it allows an efficient parallel implementation. We sketch the formulation of the SEM in matrix form that is familiar to earthquake engineers. We demonstrate the efficiency and effectiveness of the method by implementing the proposed formulation to study various simple/canonical problems.
spectral element method elastodynamics elastic wave propagation perfectly matched layer engineering seismology
Kristel C.Meza-Fajardo Apostolos S.Papageorgiou
Professor,Depto.de Ingenieria Civil,Universidad Nacional Autónoma de Honduras,Tegucigalpa,Honduras Professor,Department of Civil Engineering,University of Patras,GR-26500 Patras,Greece
国际会议
14th World Conference on Earthquake Engineering(第十四届国际地震工程会议)
北京
英文
2008-10-12（万方平台首次上网日期，不代表论文的发表时间）