Energy-Conserving Local Time Stepping Based on High-Order Finite Elements for Seismic Wave Propagation Across a Fluid-Solid Interface
Ronan Madec; Dimitri Komatitsch; and Julien Diaz;

doi:10.3970/cmes.2009.049.163
Source CMES: Computer Modeling in Engineering & Sciences, Vol. 49, No. 2, pp. 163-190, 2009
Download Full length paper in PDF format. Size = 3,296,284 bytes
Keywords Fluid-solid coupling, Time substepping, Seismic wave propagation, Spectral-element method.
Abstract When studying seismic wave propagation in fluid-solid models based on a numerical technique in the time domain with an explicit time scheme it is often of interest to resort to time substepping because the stability condition in the solid part of the medium can be more stringent than in the fluid. In such a case, one should enforce the conservation of energy along the fluid-solid interface in the time matching algorithm in order to ensure the accuracy and the stability of the time scheme. This is often not done in the available literature and approximate techniques that do not enforce the conservation of energy are used instead. We introduce such an energy-conserving local time stepping method, in which we need to solve a linear system along the fluid-solid interface. We validate it based on numerical experiments performed using high-order finite elements. This scheme can be used in any other numerical method with a diagonal mass matrix.
PDF download PDF