Journal
SIAM JOURNAL ON SCIENTIFIC COMPUTING
Volume 26, Issue 2, Pages 636-652Publisher
SIAM PUBLICATIONS
DOI: 10.1137/S1064827503427121
Keywords
bounds; PDEs; advection; convection; diffusion; reaction; exact solutions; certificates; functional outputs
Categories
Ask authors/readers for more resources
We present a cost-effective method for computing quantitative upper and lower bounds on linear functional outputs of exact weak solutions to the advection-diffusion-reaction equation and demonstrate a simple adaptive strategy by which such outputs can be computed to a prescribed precision. The bounds are computed from independent local subproblems resulting from a standard finite element approximation of the problem. At the heart of the method lies a local dual problem by which we transform an infinite dimensional minimization problem into a finite dimensional feasibility problem. The bounds hold for all levels of refinement on polygonal domains with piecewise polynomial forcing, and the bound gap converges at twice the rate of the H-1-norm of the error in the finite element solution. The bounds are valid for any choice of model problem parameters for which an equilibrating approximate solution can be found, but they become increasingly pessimistic as the parameters tend toward the advection-dominated or reaction-dominated limits.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available