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 -