Qu'en est-il de cette simple estimation d'erreur pour la PDE linéaire?


10

Soit un domaine Lipschitz borné polygonalement convexe dans , soit .R 2 f L 2 ( Ω )ΩR2fL2(Ω)

Alors la solution du problème de Dirichlet dans , sur a une solution unique dans et est bien posée, c'est-à-dire que pour un certain constant, nous avons .Ω trace u = 0 Ω H 2 C u H 2C f L 2Δu=fΩtraceu=0ΩH2CuH2CfL2

Pour une approximation d'éléments finis uh , disons, avec des éléments nodaux sur une grille uniforme, nous avons l'estimation d'erreur

uuhH1ChuH2

Il semble (peut-être que je me trompe) que les gens n'utilisent généralement pas l'estimation d'erreur évidente

uuhH1ChfL2

que nous pouvons obtenir en combinant les deux inégalités ci-dessus. Au lieu de cela, des estimateurs d'erreur a posteriori sont développés sous diverses formes. La seule objection que je puisse imaginer contre l'équation ci-dessus est que la constante C pourrait en pratique être trop pessimiste ou ne pas être estimée de manière fiable.

Réponses:


8

La raison pour laquelle les gens préfèrent utiliser la première estimation, à mon avis, est que la première découle naturellement de l'orthogonalité de Galerkin du FEM, de la propriété d'approximation d'interpolation et, surtout, de la coercivité de la forme bilinéaire (pour le problème de valeur aux limites de l'équation de Poisson , elle est équivalente à l'inégalité de Poincaré / Friedrichs pour les fonctions ): u - u h 2 H 1 ( Ω )H01

uuhH1(Ω)2c1(uuh)L2(Ω)2(uuh)L2(Ω)2=Ω(uuh)(uuh)=Ω(uuh)(uIu)(uuh)L2(Ω)(uIu)L2(Ω)(uuh)L2(Ω)(uIu)L2(Ω)c2huH2(Ω)
où dépend de la constante dans l'inégalité de Poincaré / Friedrichs pour les fonctions , est l'interpolation de dans le fini l'espace des éléments etc1H01Iuuc2 dépend des angles minimaux du maillage.

Alors que l'estimation de régularité elliptique est uniquement au niveau PDE, cela n'a rien à voir avec l'approximation, plus l'argument ci-dessus est valable même lorsque est une distribution.uH2(Ω)cfL2(Ω)fH1

Passons maintenant à la raison pour laquelle les estimations d'erreur a posteriori sont largement utilisées, principalement parce que:

  • Elle est calculable, il n'y a pas de constante générique dans l'expression des estimations.

  • L'estimateur a sa forme locale, qui pourrait être l'indicateur d'erreur local utilisé dans la procédure de raffinage adaptatif du maillage. Par conséquent, le problème des singularités ou des géométries vraiment «mauvaises» pourrait être résolu.

Les deux estimations de type a priori que vous avez énumérées sont valides, elles nous fournissent les informations des ordres de convergence, mais aucune d'entre elles ne peut être un indicateur d'erreur local juste pour un triangle / tétraèdre, car aucune d'entre elles n'est calculable en raison de la constante , ils ne sont pas non plus définis localement.

EDIT: Pour plus d'une vue générale du FEM pour les PDE elliptiques, je recommande fortement de lire le chapitre 0 du livre de Brenner et Scott: The Mathematical Theory of Finite Element Methods , qui ne comprend que 20 pages et couvre brièvement presque tous les aspects des méthodes des éléments finis , de la formulation Galerkin de la PDE, à la motivation pour laquelle nous aimerions utiliser FEM adaptatif pour résoudre un problème. J'espère que cela vous aidera davantage.


1

Votre estimation est trop pessimiste sur deux fronts. Vous avez déjà identifié la première ( inclut désormais non seulement la constante d'interpolation mais aussi la constante de stabilité). La seconde est que l'estimation d'erreur lit vraiment Notez que le côté droit a la seminorme , pas la norme. Bien sûr, vous pouvez limiter le rhs par la norme complète, mais vous perdez à nouveau de cette façon.e L 2CCH 2

eL2Ch|u|H2.
H2
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.