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 - Funding Information:
The authors thank Luis Rivera, Philip and Rachel Aber-crombie, Roland Martin, Tom Sterling, Emmanuel Chaljub, Yann Capdeville, Jan Lindheim, Ewing Lusk, Hans-Peter Bunge and Paul F. Fischer for fruitful discussions and comments. This material is based in part upon work supported by the National Science Foundation under Grant No. 0003716. This is Caltech GPS contribution No. 8823.
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 -