Adaptive finite elements for the time-dependent wave equation

Abstract : This thesis focuses on the a posteriori error analysis for the linear second-order wave equation discretized by the second order Newmark scheme in time and the finite element method in space. We adopt the particular choice for the parameters in the Newmark scheme, namely β = 1/4, γ = 1/2, since it provides a conservative method with respect to the energy norm. We derive a posteriori error estimates of optimal order in time and space for the fully discrete wave equation. The error is measured in a physically natural norm: H1 in space, Linf in time. Numerical experiments demonstrate that our error estimators are of optimal order in space and time. The resulting estimator in time is referred to as the 3-point estimator since it contains the discrete solution at 3 points in time. The 3-point time error estimator contains the Laplacian of the discrete solution which should be computed via auxiliary finite element problems at each time step. We propose an alternative time error estimator that avoids these additional computations. The resulting estimator is referred to as the 5-point estimator since it contains the fourth order finite differences in time and thus involves the discrete solution at 5 points in time at each time step. We prove that our time estimators are of optimal order at least on sufficiently smooth solutions, quasi-uniform meshes in space and uniform meshes in time. The most interesting finding of this analysis is the crucial importance of the way in which the initial conditions are discretized: a straightforward discretization, such as the nodal interpolation, may ruin the error estimators while providing quite acceptable numerical solution. We also extend the a posteriori error analysis to the general second order Newmark scheme (γ = 1/2) and present numerical comparasion between the general 3-point time error estimator and the staggered grid error estimator proposed by Georgoulis et al. In addition, using obtained a posteriori error bounds, we implement an efficient adaptive algorithm in space and time. We conclude with numerical experiments that show that the manner of interpolation of the numerical solution from one mesh to another plays an important role for optimal behavior of the time error estimator and thus of the whole adaptive algorithm.
Document type :
Liste complète des métadonnées

Cited literature [3 references]  Display  Hide  Download
Contributor : Abes Star <>
Submitted on : Monday, September 24, 2018 - 10:39:39 AM
Last modification on : Tuesday, September 25, 2018 - 1:12:13 AM
Document(s) archivé(s) le : Tuesday, December 25, 2018 - 1:07:22 PM


Version validated by the jury (STAR)


  • HAL Id : tel-01879629, version 1



Olga Gorynina. Adaptive finite elements for the time-dependent wave equation. Functional Analysis [math.FA]. Université Bourgogne Franche-Comté, 2018. English. 〈NNT : 2018UBFCD024〉. 〈tel-01879629〉



Record views


Files downloads