Post-processed Galerkin approximation of improved order for wave equations
For researchers in numerical analysis of wave equations, this provides a method to improve temporal accuracy and enable better a-posteriori error control, though it is an incremental extension of existing finite element techniques.
This paper introduces a post-processing technique for variational space-time approximations of wave equations that lifts continuous finite element solutions to continuously differentiable ones, increasing the temporal convergence order by one. The method achieves optimal-order error estimates and superconvergence at discrete time nodes, confirmed by numerical experiments.
We introduce and analyze a post-processing for a family of variational space-time approximations to wave problems. The discretization in space and time is based on continuous finite element methods. The post-processing lifts the fully discrete approximations in time from continuous to continuously differentiable ones. Further, it increases the order of convergence of the discretization in time which can be be exploited nicely, for instance, for a-posteriori error control. The convergence behavior is shown by proving error estimates of optimal order in various norms. A bound of superconvergence at the discrete times nodes is included. To show the error estimates, a special approach is developed. Firstly, error estimates for the time derivative of the post-processed solution are proved. Then, in a second step these results are used to establish the desired error estimates for the post-processed solution itself. The need for this approach comes through the structure of the wave equation providing only stability estimates that preclude us from using absorption arguments for the control of certain error quantities. A further key ingredient of this work is the construction of a new time-interpolate of the exact solution that is needed in an essential way for deriving the error estimates. Finally, a conservation of energy property is shown for the post-processed solution which is a key feature for approximation schemes to wave equations. The error estimates given in this work are confirmed by numerical experiments.