Estimating and localizing the algebraic and total numerical errors using flux reconstructions

Estimating and localizing the algebraic and total numerical errors using flux reconstructions This paper presents a methodology for computing upper and lower bounds for both the algebraic and total errors in the context of the conforming finite element discretization of the Poisson model problem and an arbitrary iterative algebraic solver. The derived bounds do not contain any unspecified constants and allow estimating the local distribution of both errors over the computational domain. Combining these bounds, we also obtain guaranteed upper and lower bounds on the discretization error. This allows to propose novel mathematically justified stopping criteria for iterative algebraic solvers ensuring that the algebraic error will lie below the discretization one. Our upper algebraic and total error bounds are based on locally reconstructed fluxes in $${\mathbf {H}}(\mathrm{div},\varOmega )$$ H ( div , Ω ) , whereas the lower algebraic and total error bounds rely on locally constructed $$H^1_0(\varOmega )$$ H 0 1 ( Ω ) -liftings of the algebraic and total residuals. We prove global and local efficiency of the upper bound on the total error and its robustness with respect to the approximation polynomial degree. Relationships to the previously published estimates on the algebraic error are discussed. Theoretical results are illustrated on numerical experiments for higher-order finite element approximations and the preconditioned conjugate gradient method. They in particular witness that the proposed methodology yields a tight estimate on the local distribution of the algebraic and total errors over the computational domain and illustrate the associated cost. http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Numerische Mathematik Springer Journals

Estimating and localizing the algebraic and total numerical errors using flux reconstructions

Loading next page...
 
/lp/springer_journal/estimating-and-localizing-the-algebraic-and-total-numerical-errors-PdxewpYBfQ
Publisher
Springer Berlin Heidelberg
Copyright
Copyright © 2017 by Springer-Verlag GmbH Deutschland
Subject
Mathematics; Numerical Analysis; Mathematics, general; Theoretical, Mathematical and Computational Physics; Mathematical Methods in Physics; Numerical and Computational Physics, Simulation; Mathematical and Computational Engineering
ISSN
0029-599X
eISSN
0945-3245
D.O.I.
10.1007/s00211-017-0915-5
Publisher site
See Article on Publisher Site

Abstract

This paper presents a methodology for computing upper and lower bounds for both the algebraic and total errors in the context of the conforming finite element discretization of the Poisson model problem and an arbitrary iterative algebraic solver. The derived bounds do not contain any unspecified constants and allow estimating the local distribution of both errors over the computational domain. Combining these bounds, we also obtain guaranteed upper and lower bounds on the discretization error. This allows to propose novel mathematically justified stopping criteria for iterative algebraic solvers ensuring that the algebraic error will lie below the discretization one. Our upper algebraic and total error bounds are based on locally reconstructed fluxes in $${\mathbf {H}}(\mathrm{div},\varOmega )$$ H ( div , Ω ) , whereas the lower algebraic and total error bounds rely on locally constructed $$H^1_0(\varOmega )$$ H 0 1 ( Ω ) -liftings of the algebraic and total residuals. We prove global and local efficiency of the upper bound on the total error and its robustness with respect to the approximation polynomial degree. Relationships to the previously published estimates on the algebraic error are discussed. Theoretical results are illustrated on numerical experiments for higher-order finite element approximations and the preconditioned conjugate gradient method. They in particular witness that the proposed methodology yields a tight estimate on the local distribution of the algebraic and total errors over the computational domain and illustrate the associated cost.

Journal

Numerische MathematikSpringer Journals

Published: Sep 6, 2017

References

You’re reading a free preview. Subscribe to read the entire article.


DeepDyve is your
personal research library

It’s your single place to instantly
discover and read the research
that matters to you.

Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.

All for just $49/month

Explore the DeepDyve Library

Search

Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly

Organize

Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.

Access

Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.

Your journals are on DeepDyve

Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.

All the latest content is available, no embargo periods.

See the journals in your area

DeepDyve

Freelancer

DeepDyve

Pro

Price

FREE

$49/month
$360/year

Save searches from
Google Scholar,
PubMed

Create lists to
organize your research

Export lists, citations

Read DeepDyve articles

Abstract access only

Unlimited access to over
18 million full-text articles

Print

20 pages / month

PDF Discount

20% off