arrow
Return

Fictitious domain method with mixed finite elements for elastodynamics

delete2007-01-01
delete5
delete
OA
AI
É
Éliane Bécache *
J
Jerónimo Rodríguez
C
Chrysoula Tsogka
DOI:10.1137/060655821delete
deleteOriginal
deleteOriginal request for help
deleteShare
deleteSave
Abstract

Abstract

En 中文
We consider in this paper the wave scattering problem by an object with Neumann boundary conditions in an anisotropic elastic body. To obtain an efficient numerical method (permitting the use of regular grids) we follow a fictitious domain approach coupled with a first order velocity stress formulation for elastodynamics. We first observe that the method does not always converge when the Q(1)(div) - Q(0) finite element is used. In particular, the method converges for some scattering object geometries but not for others. Note that the convergence of the Q(1)(div) - Q(0) finite element method was shown in [E. Becache, P. Joly, and C. Tsogka, SIAM J. Numer. Anal., 39 (2002), pp. 2109 - 2132] for the elastodynamic problem in the absence of a scattering object (i.e., without the coupling of the mixed finite elements with the fictitious domain method). Therefore we propose here a modi. cation of the Q(1)(div) - Q(0) element following the approach in [E. Becache, J. Rodriguez, and C. Tsogka, On the convergence of the fictitious domain method for wave equation problems, Technical report 5802, INRIA, 2006], where the simpler acoustic case was considered. To study the numerical properties of the new element we carry out a dispersion analysis. Several numerical simulations as well as a numerical convergence analysis show that the proposed method provides a good approximate solution.
Keywords:
mixed finite elements
fictitious domain method
elastic waves
convergence

Journal

SIAM Journal on Scientific Computing cover
SIAM Journal on Scientific Computing
IF:
2.6
Papers:
5.1K
Citations:
1.8W

Organization

No organization information available