Fictitious Domains, Mixed Finite Elements and Perfectly Matched Layers for 2D Elastic Wave Propagation
Résumé
We design a new and efficient numerical method for the modelization of elastic wave propagation in domains with complex topographies. The first specificity is the use of the fictitious domain method to take into account the boundary condition on the topography: the elastodynamic problem is extended in a domain with simple geometry, which permits the use of regular meshes. The free boundary condition is enforced introducing a Lagrange multiplier, defined on the boundary and discretized with a non uniform 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 non-standard properties: they take into account the symmetry of the stress tensor, they are compatible with mass lumping techniques and lead to explicit time discretisation schemes, and they can be coupled with the Perfectly Matched Layer technique for the modeling of unbounded domains. Our method permits to model wave propagation in complex media such as anisotropic, heterogeneous media with complex topographi- es or/and with cracks, as it will be illustratred by several numerical experiments.