An adaptive finite element method for the wave equation based on anisotropic a posteriori error estimates in the L2(H1) norm

Abstract : An adaptive finite element algorithm is presented for the wave equation in two space dimensions. The goal of the adaptive algorithm is to control the error in the same norm as for parabolic problems, namely the L2(0,T;H1(\Omega)) norm, where T denotes the final time and Omega the computational domain. The mesh aspect ratio can be large whenever needed, thus allowing a given level of accuracy to be reached with fewer vertices than with classical isotropic meshes. The refinement and coarsening criteria are based on anisotropic, a posteriori error estimates and on an elliptic reconstruction. A numerical study of the effectivity index on non-adapted meshes confirms the sharpness of the error estimator. Numerical results on adapted meshes indicate that the error indicator slightly underestimates the true error. We conjecture that the missing information corresponds to the interpolation error between successive meshes. It is observed that the error indicator becomes sharp again when considering the damped wave equation with a large damping coefficient, thus when the parabolic character of the PDE becomes predominant.
Document type :
Reports
Complete list of metadatas

https://hal.inria.fr/inria-00435786
Contributor : Marco Picasso <>
Submitted on : Wednesday, January 6, 2010 - 3:00:46 PM
Last modification on : Friday, May 25, 2018 - 12:02:03 PM
Long-term archiving on : Thursday, September 23, 2010 - 11:43:35 AM

File

RR-7115.pdf
Files produced by the author(s)

Identifiers

  • HAL Id : inria-00435786, version 2

Collections

Citation

Marco Picasso. An adaptive finite element method for the wave equation based on anisotropic a posteriori error estimates in the L2(H1) norm. [Research Report] RR-7115, INRIA. 2009, pp.22. ⟨inria-00435786v2⟩

Share

Metrics

Record views

252

Files downloads

221