Commun. Comput. Phys., 7 (2010), pp. 631-638.


On the Use of Adjoint-Based Sensitivity Estimates to Control Local Mesh Refinement

Joseph D. Frazier 1, Peter K. Jimack 2*, Robert M. Kirby 1

1 School of Computing, University of Utah, Salt Lake City, UT, USA.
2 School of Computing, University of Leeds, Leeds, LS2 9JT, UK.

Received 31 March 2009; Accepted (in revised version) 4 September 2009
Available online 9 October 2009
doi:10.4208/cicp.2009.09.060

Abstract

The goal of efficient and robust error control, through local mesh adaptation in the computational solution of partial differential equations, is predicated on the ability to identify in an a posteriori way those localized regions whose refinement will lead to the most significant reductions in the error. The development of a posteriori error estimation schemes and of a refinement infrastructure both facilitate this goal, however they are incomplete in the sense that they do not provide an answer as to where the maximal impact of refinement may be gained or what type of refinement - elemental partitioning (h-refinement) or polynomial enrichment (p-refinement) - will best lead to that gain. In essence, one also requires knowledge of the sensitivity of the error to both the location and the type of refinement. In this communication we propose the use of adjoint-based sensitivity analysis to discriminate both where and how to refine. We present both an adjoint-based and an algebraic perspective on defining and using sensitivities, and then demonstrate through several one-dimensional model problem experiments the feasibility and benefits of our approach.

AMS subject classifications: 65M60, 65M70, 65N50

Notice: Undefined variable: pac in /var/www/html/issue/abstract/readabs.php on line 164
Key words: Adaptivity, error estimation, hp-finite elements, sensitivity, discrete adjoint.

*Corresponding author.
Email: jfrazier@cs.utah.edu (J. D. Frazier), p.k.jimack@leeds.ac.uk (P. K. Jimack), kirby@sci.utah.edu (R. M. Kirby)
 

The Global Science Journal