An equilibrated residual method with a computable error approximation for a singularly perturbed reaction-diffusion problem on anisotropic finite element meshes

Sergey Grosman

ESAIM: Mathematical Modelling and Numerical Analysis (2006)

  • Volume: 40, Issue: 2, page 239-267
  • ISSN: 0764-583X

Abstract

top
Singularly perturbed reaction-diffusion problems exhibit in general solutions with anisotropic features, e.g. strong boundary and/or interior layers. This anisotropy is reflected in a discretization by using meshes with anisotropic elements. The quality of the numerical solution rests on the robustness of the a posteriori error estimator with respect to both, the perturbation parameters of the problem and the anisotropy of the mesh. The equilibrated residual method has been shown to provide one of the most reliable error estimates for the reaction-diffusion problem. Its modification suggested by Ainsworth and Babuška has been proved to be robust for the case of singular perturbation. In the present work we investigate the modified method on anisotropic meshes. The method in the form of Ainsworth and Babuška is shown here to fail on anisotropic meshes. We suggest a new modification based on the stretching ratios of the mesh elements. The resulting error estimator is equivalent to the equilibrated residual method in the case of isotropic meshes and is proved to be robust on anisotropic meshes as well. Among others, the equilibrated residual method involves the solution of an infinite dimensional local problem on each element. In practical computations an approximate solution to this local problem was successfully computed. Nevertheless, up to now no rigorous analysis has been done showing the appropriateness of any computable approximation. This demands special attention since an improper approximate solution to the local problem can be fatal for the robustness of the whole method. In the present work we provide one of the desired approximations. We prove that the method is not affected by the approximate solution of the local problem.

How to cite

top

Grosman, Sergey. "An equilibrated residual method with a computable error approximation for a singularly perturbed reaction-diffusion problem on anisotropic finite element meshes." ESAIM: Mathematical Modelling and Numerical Analysis 40.2 (2006): 239-267. <http://eudml.org/doc/249759>.

@article{Grosman2006,
abstract = { Singularly perturbed reaction-diffusion problems exhibit in general solutions with anisotropic features, e.g. strong boundary and/or interior layers. This anisotropy is reflected in a discretization by using meshes with anisotropic elements. The quality of the numerical solution rests on the robustness of the a posteriori error estimator with respect to both, the perturbation parameters of the problem and the anisotropy of the mesh. The equilibrated residual method has been shown to provide one of the most reliable error estimates for the reaction-diffusion problem. Its modification suggested by Ainsworth and Babuška has been proved to be robust for the case of singular perturbation. In the present work we investigate the modified method on anisotropic meshes. The method in the form of Ainsworth and Babuška is shown here to fail on anisotropic meshes. We suggest a new modification based on the stretching ratios of the mesh elements. The resulting error estimator is equivalent to the equilibrated residual method in the case of isotropic meshes and is proved to be robust on anisotropic meshes as well. Among others, the equilibrated residual method involves the solution of an infinite dimensional local problem on each element. In practical computations an approximate solution to this local problem was successfully computed. Nevertheless, up to now no rigorous analysis has been done showing the appropriateness of any computable approximation. This demands special attention since an improper approximate solution to the local problem can be fatal for the robustness of the whole method. In the present work we provide one of the desired approximations. We prove that the method is not affected by the approximate solution of the local problem. },
author = {Grosman, Sergey},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis},
keywords = {A posteriori error estimation; singular perturbations; reaction-diffusion problem; robustness; anisotropic solution; stretched elements.; a posteriori error estimation; stretched elements; mesh adaptivity; finite element meshes; interior layer; error bounds; numerical examples},
language = {eng},
month = {6},
number = {2},
pages = {239-267},
publisher = {EDP Sciences},
title = {An equilibrated residual method with a computable error approximation for a singularly perturbed reaction-diffusion problem on anisotropic finite element meshes},
url = {http://eudml.org/doc/249759},
volume = {40},
year = {2006},
}

TY - JOUR
AU - Grosman, Sergey
TI - An equilibrated residual method with a computable error approximation for a singularly perturbed reaction-diffusion problem on anisotropic finite element meshes
JO - ESAIM: Mathematical Modelling and Numerical Analysis
DA - 2006/6//
PB - EDP Sciences
VL - 40
IS - 2
SP - 239
EP - 267
AB - Singularly perturbed reaction-diffusion problems exhibit in general solutions with anisotropic features, e.g. strong boundary and/or interior layers. This anisotropy is reflected in a discretization by using meshes with anisotropic elements. The quality of the numerical solution rests on the robustness of the a posteriori error estimator with respect to both, the perturbation parameters of the problem and the anisotropy of the mesh. The equilibrated residual method has been shown to provide one of the most reliable error estimates for the reaction-diffusion problem. Its modification suggested by Ainsworth and Babuška has been proved to be robust for the case of singular perturbation. In the present work we investigate the modified method on anisotropic meshes. The method in the form of Ainsworth and Babuška is shown here to fail on anisotropic meshes. We suggest a new modification based on the stretching ratios of the mesh elements. The resulting error estimator is equivalent to the equilibrated residual method in the case of isotropic meshes and is proved to be robust on anisotropic meshes as well. Among others, the equilibrated residual method involves the solution of an infinite dimensional local problem on each element. In practical computations an approximate solution to this local problem was successfully computed. Nevertheless, up to now no rigorous analysis has been done showing the appropriateness of any computable approximation. This demands special attention since an improper approximate solution to the local problem can be fatal for the robustness of the whole method. In the present work we provide one of the desired approximations. We prove that the method is not affected by the approximate solution of the local problem.
LA - eng
KW - A posteriori error estimation; singular perturbations; reaction-diffusion problem; robustness; anisotropic solution; stretched elements.; a posteriori error estimation; stretched elements; mesh adaptivity; finite element meshes; interior layer; error bounds; numerical examples
UR - http://eudml.org/doc/249759
ER -

References

top
  1. M. Ainsworth and I. Babuška, Reliable and robust a posteriori error estimation for singularly perturbed reaction-diffusion problems. SIAM J. Numer. Anal.36 (1999) 331–353 (electronic). See also Corrigendum at .  Zbl0948.65114URIhttp://www.maths.strath.ac.uk/~aas98107/papers.html
  2. M. Ainsworth and J.T. Oden, A unified approach to a posteriori error estimation using element residual methods. Numer. Math.65 (1993) 23–50.  Zbl0797.65080
  3. M. Ainsworth and J.T. Oden, A Posteriori Error Estimation in Finite Element Analysis. Wiley (2000).  Zbl1008.65076
  4. T. Apel, Anisotropic interpolation error estimates for isoparametric quadrilateral finite elements. Computing60 (1998) 157–174.  Zbl0897.65003
  5. T. Apel, Treatment of boundary layers with anisotropic finite elements. Z. Angew. Math. Mech. (1998).  
  6. T. Apel, Anisotropic finite elements: local estimates and applications. B.G. Teubner, Stuttgart (1999).  
  7. T. Apel, S. Grosman, P.K. Jimack and A. Meyer, A new methodology for anisotropic mesh refinement based upon error gradients. Appl. Numer. Math.50 (2004) 329–341.  Zbl1050.65122
  8. T. Apel and G. Lube, Anisotropic mesh refinement for a singularly perturbed reaction diffusion model problem. Appl. Numer. Math.26 (1998) 415–433.  Zbl0933.65136
  9. I. Babuška and W. Rheinboldt, A posteriori error estimates for the finite element method. Int. J. Numer. Meth. Eng.12 (1978) 1597–1615.  Zbl0396.65068
  10. R. Bank and A. Weiser, Some a posteriori error estimators for elliptic partial differential equations. Math. Comp.44 (1985) 283–301.  Zbl0569.65079
  11. H. Bufler and E. Stein, Zur Plattenberechnung mittels finiter Elemente. Ingenier Archiv39 (1970) 248–260.  
  12. P.G. Ciarlet, The finite element method for elliptic problems. North-Holland Publishing Co., Amsterdam. Studies in Mathematics and its Applications, Vol. 4, (1978).  Zbl0383.65058
  13. M. Dobrowolski, S. Gräf and C. Pflaum, On a posteriori error estimators in the infinte element method on anisotropic meshes. Electron. Trans. Numer. Anal. 8 (1999) 36–45.  Zbl0934.65122
  14. S. Grosman, The robustness of the hierarchical a posteriori error estimator for reaction-diffusion equation on anisotropic meshes. SFB393-Preprint 2, Technische Universität Chemnitz, SFB 393 (Germany), (2004).  
  15. R. Hagen, S. Roch, and B. Silbermann, C*-algebras and numerical analysis. Marcel Dekker Inc., New York (2001).  Zbl0964.65055
  16. H. Han and R.B. Kellogg, Differentiability properties of solutions of the equation - ϵ 2 δ u + r u = f ( x , y ) in a square. SIAM J. Math. Anal.21 (1990) 394–408.  Zbl0732.35020
  17. G. Kunert, A posteriori error estimation for anisotropic tetrahedral and triangular finite element meshes. Logos Verlag, Berlin, 1999. Also PhD thesis, TU Chemnitz, .  Zbl0919.65066URIhttp://archiv.tu-chemnitz.de/pub/1999/0012/index.html
  18. G. Kunert, An a posteriori residual error estimator for the finite element method on anisotropic tetrahedral meshes. Numer. Math.86 (2000) 471–490.  Zbl0965.65125
  19. G. Kunert, A local problem error estimator for anisotropic tetrahedral finite element meshes. SIAM J. Numer. Anal.39 (2001) 668–689.  Zbl1004.65112
  20. G. Kunert, Robust a posteriori error estimation for a singularly perturbed reaction-diffusion equation on anisotropic tetrahedral meshes. Adv. Comput. Math.15 (2001) 237–259.  Zbl1049.65121
  21. G. Kunert, Robust local problem error estimation for a singularly perturbed problem on anisotropic finite element meshes. ESAIM: M2AN35 (2001) 1079–1109.  Zbl1041.65072
  22. G. Kunert and R. Verfürth, Edge residuals dominate a posteriori error estimates for linear finite element methods on anisotropic triangular and tetrahedral meshes. Numer. Math.86 (2000) 283–303.  Zbl0964.65120
  23. P. Ladevèze and D. Leguillon, Error estimate procedure in the finite element method and applications. SIAM J. Numer. Anal.20 (1983) 485–509.  Zbl0582.65078
  24. K.G. Siebert, An a posteriori error estimator for anisotropic refinement. Numer. Math.73 (1996) 373–398.  Zbl0873.65098
  25. R. Verfürth, A review of a posteriori error estimation and adaptive mesh-refinement techniques. Wiley-Teubner Series Advances in Numerical Mathematics. Chichester: John Wiley & Sons. Stuttgart: B.G. Teubner (1996).  
  26. M. Vogelius and I. Babuška, On a dimensional reduction method. I. The optimal selection of basis functions. Math. Comp.37 (1981) 31–46.  Zbl0495.65049

NotesEmbed ?

top

You must be logged in to post comments.

To embed these notes on your page include the following JavaScript code on your page where you want the notes to appear.

Only the controls for the widget will be shown in your chosen language. Notes will be shown in their authored language.

Tells the widget how many notes to show per page. You can cycle through additional notes using the next and previous controls.

    
                

Note: Best practice suggests putting the JavaScript code just before the closing </body> tag.