Méthode de Galerkine

Cet article est une ébauche concernant l’analyse.

Vous pouvez partager vos connaissances en l’améliorant (comment ?) selon les recommandations des projets correspondants.

En mathématiques, dans le domaine de l'analyse numérique, les méthodes de Galerkine sont une classe de méthodes permettant de transformer un problème continu (par exemple une équation différentielle) en un problème discret. Cette approche est attribuée aux ingénieurs russes Ivan Boubnov (1911) et Boris Galerkine (1913).

Approximation de fonctions[modifier | modifier le code]

Cette méthode est couramment utilisée dans la méthode des éléments finis.

On part de la formulation faible du problème. La solution appartient à un espace fonctionnel satisfaisant des propriétés de régularité bien définies. La méthode de Galerkine consiste à utiliser un maillage du domaine d'étude, et considérer la restriction de la fonction solution sur chacune des mailles.

D'un point de vue plus formel, on écrit la formulation faible sous la forme :

Trouver

u\in V

telle que

\forall v\in V,a(u,v)=L(v)

où $a$ est une forme bilinéaire, et $L$ une forme linéaire.

L'ensemble $V$ étant généralement de dimension infinie, on construit un espace $V_{h}\subset V$ avec ${\rm {dim}}V_{h}<+\infty$ , et on réécrit le problème de la façon suivante :

Trouver

u_{h}\in V_{h}

telle que

\forall v_{h}\in V_{h},a(u_{h},v_{h})=L(v_{h})

Typiquement, l'espace $V_{h}$ considéré est l'ensemble des fonctions continues telles que la restriction de la fonction sur une maille soit un polynôme.

Propriétés[modifier | modifier le code]

Orthogonalité de l'erreur[modifier | modifier le code]

L'une des propriétés notables des méthodes de Galerkine se trouvent dans le fait que l'erreur commise sur la solution $e_{h}=u-u_{h}$ est orthogonale aux sous-espaces d'approximation. En effet, les propriétés de la forme bilinéaire a donnent :

\forall v_{h}\in V_{h}\ ,\ a(e_{h},v_{h})=a(u-u_{h},v_{h})=a(u,v_{h})-a(u_{h},v_{h})=L(v_{h})-L(v_{h})=0

.

Forme matricielle du problème[modifier | modifier le code]

Du fait que l'espace d'approximation utilisé $V_{h}$ est de dimension finie $n<+\infty$ , on peut décomposer la solution du problème de Galerkine sur une base de fonctions $\left\lbrace e_{i}\right\rbrace _{1\leqslant i\leqslant n}$ de $V_{h}$ :

u_{h}=\sum _{j=1}^{n}u_{j}e_{j}

Ainsi, en écrivant le problème en choisissant l'une des fonctions de base $v_{h}=e_{i}$ , il vient :

a(u_{h},e_{i})=\sum _{j=1}^{n}u_{j}a(e_{j},e_{i})=L(e_{i})\ ,\ \forall i\in [\![1,n]\!]

On obtient ainsi un système d'équations linéaires de la forme $Au=l$ , en notant

A=\left(a(e_{j},e_{i})\right)_{1\leqslant i,j\leqslant n}

,

u=\left(u_{i}\right)_{1\leqslant i\leqslant n}

,

l=\left(L(e_{i})\right)_{1\leqslant i\leqslant n}

Systèmes symétriques et positifs[modifier | modifier le code]

Il apparait que si la forme bilinéaire $a$ est symétrique, la matrice $A$ est également symétrique. De même, $A$ est une matrice positive (définie positive) si $a$ l'est également.

Résultats sur la solution obtenue[modifier | modifier le code]

Existence et unicité

Dans le cas où $a$ est symétrique, on peut montrer que la solution du problème existe et est unique si on a :

continuité de $a$ sur $V$

\forall (u,v)\in V\times V,\,|a(u,v)|\leqslant C\|u\|\|v\|

coercivité de $a$ sur $V_{h}$

\forall u\in V,\,a(u,u)\geqslant c\|u\|^{2}

Il suffit alors d'appliquer le théorème de Lax-Milgram pour obtenir le résultat voulu.

Le caractère bien posé du problème écrit sur $V_{h}$ en découle naturellement.

Qualité de l'approximation

En utilisant les mêmes propriétés de $a$ , ainsi que l'orthogonalité de l'erreur, on obtient l'inégalité pour tout $v_{h}\in V_{h}$ :

c\|u-u_{h}\|^{2}\leqslant a(u-u_{h},u-u_{h})\leqslant a(u-u_{h},u-v_{h})\leqslant C\|u-u_{h}\|\|u-v_{h}\|

.

En divisant par $\|u-u_{h}\|$ et en passant à la borne inférieure à droite, on obtient le lemme de Céa :

\|u-u_{h}\|\leqslant {\frac {C}{c}}\inf _{v_{h}\in V_{h}}\|u-v_{h}\|

Ainsi, à la constante $C / c$ près, la solution obtenue par la méthode de Galerkine est une des meilleures qu'on puisse obtenir par approximation sur $V_{h}$ .

Voir aussi[modifier | modifier le code]

Référence[modifier | modifier le code]

(en) Philippe Ciarlet, The Finite Element Method for Elliptic Problems, avril 2002, 530 p. (ISBN 978-0-89871-514-9, lire en ligne)
Jean-Christophe Cuillière, Introduction à la méthode des éléments finis - 2e édition, Paris, Dunod, 2016, 288 p. (ISBN 978-2-10-074262-2), chap. 6.5 (« Hypothèse de Galerkine »), p. 104-112

Portail de l'analyse