carrier image

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

Casarotti, Emanuele, Marco Stupazzini, Shiann Jong Lee, Dimitri Komatitsch, Antonio Piersanti, and Jeroen Tromp

Proceedings, 16th International Meshing Roundtable, Springer-Verlag, pp.579-598, October 14-17 2007

IMR
PROCEEDINGS

16th International Meshing Roundtable
Seattle, Washington, U.S.A.
October 14-17, 2007

Seismological Laboratory, California Institute of Technology. 1200 E. California
Blvd., MS 252-21, Pasadena, California, 91125, USA.
[email]emanuele@gps.caltech.edu
Department of Earth- and Environmental Sciences, Geophysics section,
Ludwig-Maximilians-Universit®at. Theresienstrasse 41, 80333, Munich, Germany.
[email]stupa@geophysik.uni-muenchen.de
Institute of Earth Science, Academia Sinica. No.128, Section2, Academia Road,
Nankang, Taipei, Taiwan 115.
sjlee@earth.sinica.edu.tw
Laboratoire de Modlisation et díImagerie en Gosciences UMR 5212, Universit de
Pau et des Pays de líAdour. Btiment IPRA - Avenue de líUniversit, BP 1155,64013, Pau, France.
dimitri.komatitsch@univ-pau.fr
Istituto Nazionale di Geofisica e Vulcanologia. Via di Vigna Murata 605, 00143,Rome, Italy.
[email]piersanti@ingv.it
Seismological Laboratory, California Institute of Technology. 1200 E. California
Blvd., MS 252-21, Pasadena, California, 91125, USA. jtromp@gps.caltech.edu

Abstract
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, cubit.sandia.gov), 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 solidsolid 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).

Download Full Paper (PDF Format)


Contact author(s) or publisher for availability and copyright information on above referenced article