8494 articles  [version française]

inria-00435786, version 2

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

Marco Picasso () 1

N° RR-7115 (2009)

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.

  • 1:  GAMMA (INRIA Rocquencourt)
  • INRIA
  • Domain : Computer Science/Modeling and Simulation
  • Keywords : Wave equation – Adaptive finite elements – Anisotropic a posteriori error estimates
  • Internal note : RR-7115
  • Available versions :  v1 (2009-11-24) v2 (2010-01-06)
 
  • inria-00435786, version 2
  • oai:hal.inria.fr:inria-00435786
  • From: 
  • Submitted on: Wednesday, 6 January 2010 15:00:46
  • Updated on: Tuesday, 2 February 2010 11:30:53