Une technique variationnelle numérique qui approchera des solutions de problèmes aux limites en remplaçant des espaces d'essai de dimension infinie par des espaces polynomiaux par morceaux de dimension finie définis sur un maillage, conduisant à des systèmes algébriques creux résolus par solveurs linéaires ou non linéaires.

Principe

Discrétiser la forme faible : choisir un maillage et un degré polynomial, assembler les contributions locales des éléments en une matrice de raideur (et de masse) globale via des intégrales variationnelles, et résoudre le système algébrique obtenu ; la convergence découle des propriétés d'approximation et de stabilité des espaces discrets.

Démonstration

Pour résoudre −Δu=f avec Dirichlet homogène sur un domaine polygonal : trianguler le domaine, choisir des fonctions de base P1 (linéaires par morceaux) à support local, assembler la matrice de raideur A_{ij}=∫Ω ∇ϕ_i·∇ϕ_j et le vecteur charge b_i=∫Ω fϕ_i, puis résoudre Au=b pour les valeurs nodales de u_h.

Mauvaise application

Employer des éléments non conformes sans correction, omettre un raffinement là où apparaissent des singularités, ou ignorer les conditions de stabilité pour problèmes mixtes (inf-sup) produit des approximations inexactes ou non convergentes.

Conséquence

Conduit à des approximations pratiques et convergentes avec estimations d'erreur calculables, stratégies d'adaptation, et solveurs efficaces ; permet de traiter des géométries complexes et des simulations à l'échelle industrielle lorsqu'il est correctement mis en œuvre.

Inversion

Abandonner la base polynomiale locale et la structure en maillage (par ex. recourir à des bases spectrales globales ou méthodes sans maillage) modifie les propriétés d'approximation, le conditionnement et les compromis d'implémentation et peut être préférable ou non selon la régularité et la classe de problème.

Limite

S'applique aux problèmes admettant une formulation variationnelle sur des domaines où l'on peut construire des maillages ; exclut les approches purement sans maillage, certains opérateurs non locaux sans localisation adéquate et les cas où le choix d'élément doit être adapté (coefficients singuliers, domaines irréguliers).

Tension sémantique

Souvent comparée aux méthodes aux différences finies, volumes finis et spectrales : la FEM privilégie la consistance variationnelle, le support local et l'adaptivité de maillage, tandis que les alternatives insistent sur la simplicité, la conservation ou la convergence exponentielle pour solutions lisses.

Synthèse

La méthode des éléments finis transforme une EDP variationnelle en un système algébrique de dimension finie par approximation polynomiale locale sur un maillage, équilibrant ordre d'approximation, résolution de maillage et stabilité pour obtenir des solutions numériques fiables et un contrôle d'erreur.