TY - GEN

T1 - Modeling of Seismic Wave Propagation at the Scale of the Earth on a Large Beowulf

AU - Komatitsch, Dimitri

AU - Tromp, Jeroen

N1 - Publisher Copyright:
© 2001 ACM.

PY - 2001/11/10

Y1 - 2001/11/10

N2 - We use a parallel spectral-element method to simulate the propagation of seismic waves generated by earthquakes in the entire 3-D Earth. The method is implemented using MPI on a large PC cluster (Beowulf) with 151 processors and 76 Gb of RAM. It is based upon a weak formulation of the equations of motion and combines the flexibility of a finite-element method with the accuracy of a pseudospectral method. The finite-element mesh honors all discontinuities in the Earth velocity model. To maintain a relatively constant number of grid points per seismic wavelength, the size of the elements is increased with depth in a conforming fashion, thus retaining a diagonal mass matrix. The effects of attenuation and anisotropy are incorporated. We benchmark spectral-element synthetic seismograms against a normalmode reference solution for a spherically symmetric Earth velocity model. The two methods are in excellent agreement for all waves with periods greater than 20 seconds.

AB - We use a parallel spectral-element method to simulate the propagation of seismic waves generated by earthquakes in the entire 3-D Earth. The method is implemented using MPI on a large PC cluster (Beowulf) with 151 processors and 76 Gb of RAM. It is based upon a weak formulation of the equations of motion and combines the flexibility of a finite-element method with the accuracy of a pseudospectral method. The finite-element mesh honors all discontinuities in the Earth velocity model. To maintain a relatively constant number of grid points per seismic wavelength, the size of the elements is increased with depth in a conforming fashion, thus retaining a diagonal mass matrix. The effects of attenuation and anisotropy are incorporated. We benchmark spectral-element synthetic seismograms against a normalmode reference solution for a spherically symmetric Earth velocity model. The two methods are in excellent agreement for all waves with periods greater than 20 seconds.

UR - http://www.scopus.com/inward/record.url?scp=85129592964&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85129592964&partnerID=8YFLogxK

U2 - 10.1145/582034.582076

DO - 10.1145/582034.582076

M3 - Conference contribution

AN - SCOPUS:85129592964

T3 - Proceedings of the International Conference on Supercomputing

SP - 42

BT - Proceedings of the 2001 ACM/IEEE Conference on Supercomputing, SC 2001

PB - Association for Computing Machinery

T2 - 2001 ACM/IEEE Conference on Supercomputing, SC 2001

Y2 - 10 November 2001 through 16 November 2001

ER -