FICTITIOUS DOMAINS, MIXED FINITE ELEMENTS AND PERFECTLY MATCHED LAYERS FOR 2-D ELASTIC WAVE PROPAGATION

We design a new and efficient numerical method for the modelization of elastic wave propagation in domains with complex topographies. The main characteristic is the use of the fictitious domain method for taking into account the boundary condition on the topography: the elastodynamic problem is extended in a domain with simple geometry, which permits us to use a regular mesh. The free boundary condition is enforced introducing a Lagrange multiplier, defined on the boundary and discretized with a nonuniform boundary mesh. This leads us to consider the first-order velocity-stress formulation of the equations and particular mixed finite elements. These elements have three main nonstandard properties: they take into account the symmetry of the stress tensor, they are compatible with mass lumping techniques and lead to explicit time discretization schemes, and they can be coupled with the Perfectly Matched Layer technique for the modeling of unbounded domains. Our method permits us to model wave propagation in complex media such as anisotropic, heterogeneous media with complex topographies, as it will be illustrated by several numerical experiments.

[1]  R. Glowinski,et al.  Error analysis of a fictitious domain method applied to a Dirichlet problem , 1995 .

[2]  P. Raviart,et al.  A mixed finite element method for 2-nd order elliptic problems , 1977 .

[3]  Chrysoula Tsogka,et al.  Application of the PML absorbing layer model to the linear elastodynamic problem in anisotropic hete , 1998 .

[4]  P. Grisvard Singularities in Boundary Value Problems , 1992 .

[5]  Eliane Bécache,et al.  Éléments finis mixtes et condensation de masse en élastodynamique linéaire. (I) Construction , 1997 .

[6]  Michel Fortin,et al.  Mixed and Hybrid Finite Element Methods , 2011, Springer Series in Computational Mathematics.

[7]  Patrick Joly,et al.  An Analysis of New Mixed Finite Elements for the Approximation of Wave Propagation Problems , 2000, SIAM J. Numer. Anal..

[8]  R. Glowinski,et al.  A fictitious domain method for Dirichlet problem and applications , 1994 .

[9]  R. Stenberg A family of mixed finite elements for the elasticity problem , 1988 .

[10]  J. Nédélec A new family of mixed finite elements in ℝ3 , 1986 .

[11]  J. Douglas,et al.  PEERS: A new mixed finite element for plane elasticity , 1984 .

[12]  Patrick Joly,et al.  Fictitious Domain Method for Unsteady Problems , 1997 .

[13]  M. E. Morley A family of mixed finite elements for linear elasticity , 1989 .

[14]  Jean-Pierre Berenger,et al.  A perfectly matched layer for the absorption of electromagnetic waves , 1994 .

[15]  Eliane Bécache,et al.  A space-time variational formulation for the boundary integral equation in a 2D elastic crack problem , 1994 .

[16]  M. A. Dablain,et al.  The application of high-order differencing to the scalar wave equation , 1986 .