×

Large-scale simulation of elastic wave propagation in heterogeneous media on parallel computers. (English) Zbl 0961.74056

Summary: This paper reports on the development of a parallel numerical methodology for simulating large-scale earthquake-induced ground motion in highly heterogeneous basins. We target large sedimentary basins with contrasts in wavelengths of over an order of magnitude. Regular grid methods prove intractable for such problems. We overcome the problem of multiple physical scales by using unstructured finite elements on locally-resolved Delaunay triangulations derived from octree-based grids. The extremely large mesh sizes require special mesh generation techniques. Despite the method’s multiresolution capability, large problem sizes necessitate the use of distributed memory parallel supercomputers to solve the elastic wave propagation problem. We develop a system that helps write efficient portable unstructured mesh solvers for distributed memory parallel supercomputers. The numerical methodology and software system have heen used to simulate the seismic response of the San Fernando Valley in Southern California to an aftershock of the 1994 Northridge Earthquake. We report on parallel performance on the Cray T3D for several models of the basin ranging in size from 35 000 to 77 million tetrahedra. The results indicate that, despite the highly irregular structure of the problem, excellent performance and scalability are achieved.

MSC:

74S05 Finite element methods applied to problems in solid mechanics
74J10 Bulk waves in solid mechanics
86A15 Seismology (including tsunami modeling), earthquakes
65Y05 Parallel numerical computation

Software:

Triangle
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Aki, K., Local site effect on ground motion, (Lawrence Von Thun, J., Earthquake Engineering and Soil Dynamics. II: Recent Advances in Ground-Motion Evaluation (1988), ASCE), 103-155
[2] Bao, H.; Bielak, J.; Ghattas, O.; Kallivokas, L. F.; O’Hallaron, D. R.; Shewchuk, J. R.; Xu, J., Earthquake ground motion modeling on parallel computers, (Proc. Supercomputing ’96 (November 1996))
[3] Bardet, J. P.; Davis, C., Engineering observations on ground motions at the Van Norman Complex after the 1994 Northridge Earthquake, Bull. Seism. Soc. Am., 86, 1B, S333-S349 (1996)
[4] Bielak, J.; Christiane, P., On the effective seismic input for nonlinear soil-structure interaction systems, Earthq. Engrg. Struct. Dyn., 12, 107-119 (1984)
[5] Bowyer, A., Computing Dirichlet tessellations, Comput. J., 24, 2, 162-166 (1981)
[6] Cremonini, M. G.; Christiano, P.; Bielak, J., Implementation of effective seismic input for soil-structure interaction systems, Earthq. Engrg. Struct. Dyn., 16, 615-625 (1988)
[7] Available from
[8] Frankel, A.; Vidale, J. E., A three-dimensional simulation of seismic waves in the Santa Clara Valley, California from a Loma Prieta aftershock, Bull. Seism. Soc. Am., 82, 2045-2074 (1992)
[9] Graves, R. W., Modeling three-dimensional site response effects in the Marina District Basin, San Francisco, California, Bull. Seism. Soc. Am., 83, 1042-1063 (1993)
[10] Hartzell, S.; Leeds, A.; Frankel, A.; Michael, J., Site response for urban Los Angeles using aftershocks of the Northridge earthquake, Bull. Seism. Soc. Am., 86, 1B, S168-S192 (1996)
[11] Kawase, H.; Aki, K., A study on the response of a soft basin for incident S, P and Rayleigh waves with special reference to the long duration observed in Mexico City, Bull. Seism. Soc. Am., 79, 1361-1382 (1989)
[12] Liu, H.-L.; Heaton, T., Array analysis of the ground velocities and accelerations from the 1971 San Fernando, California, earthquake, Bull. Seism. Soc. Am., 74, 1951-1968 (1996)
[13] Magistrale, H.; McLaughlin, K. L.; Day, S. M., A geology-based 3-D velocity model of the Los Angeles Basin sediments, Bull. Seism. Soc. Am., 86, 1161-1166 (1996)
[14] McLaughlin, K. L.; Day, S. M., 3D elastic finite difference seismic wave simulations, Comput. Phys. (Nov/Dec 1994)
[15] Miller, G. L.; Teng, S.-H.; Thurston, W.; Vavasis, S. A., Automatic mesh partitioning, (George, A.; Gilbert, J.; Liu, J., Graph Theory and Sparse Matrix Computation. Graph Theory and Sparse Matrix Computation, The IMA Volumes in Mathematics and its Application, Vol. 56 (1993), Springer-Verlag), 57-84 · Zbl 0803.68083
[16] Olsen, K. B.; Archuleta, R. J., Three-dimensional simulation of earthquakes on the Los Angeles Fault System, Bull. Seism. Soc. Am., 86, 575-596 (1996)
[17] Olsen, K. B.; Archuleta, R. J.; Matarese, J. R., Magnitude 7.75 earthquake on the San Andreas fault: Three-dimensional ground motion in Los Angeles, Science, 270, 5242, 1628-1632 (1995)
[18] Sánchez-Sesma, F. J.; Luzón, F., Seismic response of three-dimensional valleys for incident P, S, and Rayleigh waves, Bull. Seism. Soc. Am., 85, 269-284 (1995)
[19] Spudich, P.; Lida, M., The seismic coda, side effects and scattering in alluvial basins studied using aftershocks of the 1986 North Palm Springs, California earthquake as source arrays, Bull. Seism. Soc. Am., 83, 1721-1724 (1993)
[20] Shewchuk, J. R., Robust adaptive floating-point geometric predicates, (Proc. Twelfth Annual Symposium on Computational Geometry (May 1996), Association for Computing Machinery), 141-150
[21] Shewchuk, J. R., Triangle: Engineering a 2D quality mesh generator and Delaunay triangulator, (First Workshop on Applied Computational Geometry (May 1996), Association for Computing Machinery), 124-133
[22] Shewchuk, J. R.; Ghattas, O., A compiler for parallel finite element methods with domain-decomposed unstructured meshes, (Keyes, D. E.; Xu, J., Domain Decomposition Methods in Scientific and Engineering Computing. Domain Decomposition Methods in Scientific and Engineering Computing, Contemporary Mathematics, Vol. 180 (1994), American Mathematical Society), 445-450 · Zbl 0817.65079
[23] Thio, H. K.; Kanamori, H., Source complexity of the 1994 Northridge Earthquake and its relation to aftershock mechanisms, Bull. Seism. Soc. Am., 86, 1B, S84-S92 (1996)
[24] Vidale, J. E.; Helmberger, D. V., Elastic finite-difference modeling of the 1971 San Fernando, California, earthquake, Bull. Seism. Soc. Am., 78, 122-141 (1988)
[25] Watson, D. F., Computing the \(n\)-dimensional Delaunay tessellation with application to Voronoï polytopes, Comput. J., 24, 2, 167-172 (1981)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.