CUBIT and Seismic Wave Propagation Based Upon the Spectral-Element Method: An Advanced Unstructured Mesher for Complex 3D Geological Media

Unstructured hexahedral mesh generation is a critical part of the modeling process in the Spectral-Element Method (SEM). We present some examples of seismic wave propagation in complex geological models, automatically meshed on a parallel machine based upon CUBIT (Sandia Laboratory), an advanced 3D unstructured hexahedral mesh generator that offers new opportunities for seismologist to design, assess, and improve the quality of a mesh in terms of both geometrical and numerical accuracy. The main goal is to provide useful tools for understanding seismic phenomena due to surface topography and subsurface structures such as low wave-speed sedimentary basins. Our examples cover several typical geophysical problems: 1) “layer-cake” volumes with high-resolution topography and complex solid-solid interfaces (such as the Campi Flegrei Caldera Area in Italy), and 2) models with an embedded sedimentary basin (such as the Taipei basin in Taiwan or the Grenoble Valley in France).

[1]  Roland Martin,et al.  WAVE PROPAGATION IN 2-D ELASTIC MEDIA USING A SPECTRAL ELEMENT METHOD WITH TRIANGLES AND QUADRANGLES , 2001 .

[2]  R. Courant,et al.  Über die partiellen Differenzengleichungen der mathematischen Physik , 1928 .

[3]  D. Komatitsch,et al.  Simulations of Ground Motion in the Los Angeles Basin Based upon the Spectral-Element Method , 2004 .

[4]  Guust Nolet,et al.  Seismic earth : array analysis of broadband seismograms , 2005 .

[5]  M. Bouchon,et al.  Boundary Integral Equations and Boundary Elements Methods in Elastodynamics , 2007 .

[6]  D. Komatitsch,et al.  Simulation of anisotropic wave propagation based upon a spectral element method , 2000 .

[7]  Guido Russo,et al.  A 3D velocity model for earthquake location in Campi Flegrei area: application to the 1982-84 uplift event , 2006 .

[8]  Ru-Shan Wu,et al.  Advances in wave propagation in heterogeneous earth , 2007 .

[9]  M. Dumbser,et al.  An arbitrary high-order Discontinuous Galerkin method for elastic waves on unstructured meshes — III. Viscoelastic attenuation , 2007 .

[10]  J. Bielak,et al.  Domain Reduction Method for Three-Dimensional Earthquake Modeling in Localized Regions, Part II: Verification and Applications , 2001 .

[11]  Steven E. Benzley,et al.  Adaptive sweeping techniques , 2009, Engineering with Computers.

[12]  P. Moczo,et al.  The finite-difference time-domain method for modeling of seismic wave propagation , 2007 .

[13]  C. di Prisco,et al.  Nonlinear SEM numerical analyses of dry dense sand specimens under rapid and dynamic loading , 2007 .

[14]  O. Ghattas,et al.  Parallel Octree-Based Finite Element Method for Large-Scale Earthquake Ground Motion Simulation , 2005 .

[15]  Aldo Zollo,et al.  Geophysical Exploration of the Campi Flegrei (Southern Italy) Caldera' Interiors: Data, Methods and Results , 2006 .

[16]  Matthew L. Staten,et al.  Unconstrained Paving and Plastering: Progress Update , 2006, IMR.

[17]  P. Paolucci,et al.  The “Cubed Sphere” , 1996 .