Accueil🇫🇷Chercher

Méthode de Galerkine

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

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 telle que

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 avec , et on réécrit le problème de la façon suivante :

Trouver telle que

Typiquement, l'espace 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

Orthogonalité de l'erreur

L'erreur est montrée ici orthogonale à l'espace d'approximation

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

.

Forme matricielle du problème

Du fait que l'espace d'approximation utilisé est de dimension finie , on peut décomposer la solution du problème de Galerkine sur une base de fonctions de :

Ainsi, en écrivant le problème en choisissant l'une des fonctions de base , il vient :

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

, ,

Systèmes symétriques et positifs

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

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
  • coercivité de a sur

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 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 :

.

En divisant par et en passant à la borne inférieure à droite, on obtient le lemme de Céa :

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 .

Voir aussi

Référence

  • (en) Philippe Ciarlet, The Finite Element Method for Elliptic Problems, , 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, , 288 p. (ISBN 978-2-10-074262-2), chap. 6.5 (« Hypothèse de Galerkine »), p. 104-112
Cet article est issu de wikipedia. Text licence: CC BY-SA 4.0, Des conditions supplémentaires peuvent s’appliquer aux fichiers multimédias.