Open AccessForward and adjoint simulations of seismic wave propagation on fully unstructured hexahedral meshes
Author(s) -
Peter Daniel,
Komatitsch Dimitri,
Luo Yang,
Martin Roland,
Le Goff Nicolas,
Casarotti Emanuele,
Le Loher Pieyre,
Magi Federica,
Liu Qinya,
Blitz Céline,
NissenMeyer Tarje,
Basini Piero,
Tromp Jeroen
Publication year - 2011
Publication title -
geophysical journal international
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.302
H-Index - 168
eISSN - 1365-246X
pISSN - 0956-540X
DOI - 10.1111/j.1365-246x.2011.05044.x
Subject(s) - polygon mesh , hexahedron , finite element method , domain decomposition methods , mesh generation , computational science , computer science , wave propagation , algorithm , physics , structural engineering , engineering , computer graphics (images) , quantum mechanics
SUMMARY We present forward and adjoint spectral‐element simulations of coupled acoustic and (an)elastic seismic wave propagation on fully unstructured hexahedral meshes. Simulations benefit from recent advances in hexahedral meshing, load balancing and software optimization. Meshing may be accomplished using a mesh generation tool kit such as CUBIT, and load balancing is facilitated by graph partitioning based on the SCOTCH library. Coupling between fluid and solid regions is incorporated in a straightforward fashion using domain decomposition. Topography, bathymetry and Moho undulations may be readily included in the mesh, and physical dispersion and attenuation associated with anelasticity are accounted for using a series of standard linear solids. Finite‐frequency Fréchet derivatives are calculated using adjoint methods in both fluid and solid domains. The software is benchmarked for a layercake model. We present various examples of fully unstructured meshes, snapshots of wavefields and finite‐frequency kernels generated by Version 2.0 ‘Sesame’ of our widely used open source spectral‐element package SPECFEM3D.