Time and space adaptivity of the wave equation discretized in time by a second-order scheme

Marco Picasso, Olga Gorynina
2019
Journal paper

Abstract

The aim of this paper is to obtain a posteriori error bounds of optimal order in time and space for the linear second-order wave equation discretized by the Newmark scheme in time and the finite element method in space. An error estimate is derived in the L-infinity-in-time/energy-in-space norm. Numerical experiments are reported for several test cases and confirm equivalence of the proposed estimator and the true error.

Official source

https://infoscience.epfl.ch/record/271813?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts (8)

Related publications (2)

Related MOOCs (4)

We propose a cheaper version of a posteriori error estimator from Gorynina et al. (Namer. Anal. (2017)) for the linear second-order wave equation discretized by the Newmark scheme in time and by the finite element method in space. The new estimator preserves all the properties of the previous one (reliability, optimality on smooth solutions and quasi-uniform meshes) but no longer requires an extra computation of the Laplacian of the discrete solution on each time step.

EDP Sciences2019

We developed a space and time adaptive method to simulate electroosmosis and mass transport of a sample concentration within a network of microchannels. The space adaptive criteria is based on an error estimator derived using anisotropic interpolation estimates and a post-processing procedure. In order to improve the accuracy of the numerical solution and to reduce even further the computational cost of the numerical simulation, a time adaptive procedure is combined with the one in space. To do so, a time error estimator is derived for a first model problem, the linear heat equation discretized in time with the Crank-Nicolson scheme. The main difficulty is then to obtain an optimal second order error estimator. Applying standard energy techniques with a continuous, piecewise linear approximation in time fail in recovering the optimal order. To restore the appropriate rate of convergence, a continuous piecewise quadratic polynomial function in time is needed. For this purpose, two different quadratic functions are introduced and two different time error estimators are then derived. It turns out that the second error estimator is more efficient than the first one when considering our adaptive algorithm. Thus, using the second quadratic polynomial, an upper bound for the error is derived for a second model problem, the time-dependent convection-diffusion problem discretized in time with the Crank-Nicolson scheme. The corresponding space and time error estimators are finally used for the numerical simulation of mass transport of a sample concentration within a complex network of microchannels driven by an electroosmotic flow and/or by a pressure-driven flow. Numerical results presented show the efficiency and the robustness of this approach.

EPFL2009