Méthode
des différences
finies
MDF
Plan du cours
principes généraux
Différences finies 1D
Différences finies 2D
Etude de quelques cas réels
principes généraux
Introduction
L’évolution actuelle de la technologie amène l’ingénieur à réaliser des projets de plus en
plus complexes, coûteux, et soumis à des contraintes de sécurité de plus en plus sévères.
Pour dominer ces projets, l’ingénieur a besoin des modèles qui lui permettent de simuler le
comportement de systèmes physiques complexes.
Des problèmes, considérés insolubles avec les méthodes classiques, sont maintenant
facilement accessibles avec les méthodes numériques.
Une de ces méthodes est la Méthode des Différences Finies (MDF).
Introduction
En analyse numérique, la méthode des différences finies est une technique courante de
recherche de solutions approchées d'équations aux dérivées partielles qui consiste à résoudre
un système de relations (schéma numérique) liant les valeurs des fonctions inconnues en
certains points suffisamment proches les uns des autres.
la plus ancienne méthodes d’approximation pour la résolution des équations différentielles
(ordinaires et partielles). Elle est introduite en 18 siècle (en 1768) par Euler.
C’ est la plus simple des nombreuses méthode d'approximation et la méthode la plus
efficace pour les problèmes avec des géométries de régions simples.
Elle se base sur la solution directe de l’EDP (pas d’ équations d'intégrale).
Introduction
Objectif : transformer une équation « continue » valable sur un domaine
continu en un système à N équations pour N inconnues associées
à un domaine discret appelé maillage
Méthode : écrire sous forme discrète (i-1, i, i+1 …) tous les termes de
dérivées présents dans l’équation d’équilibre appliquée en i
ainsi que dans les C.L.
Introduction
Intérêts des méthodes de différences finies
La méthode des différences finies est une méthode de résolutions d'EDP couramment
pratiquée, car elle est la plus facile d'accès, puisqu'elle repose sur deux notions bien connus:
la discrétisation des opérateurs de dérivation/différentiation (assez intuitive) d'une part,
la convergence du schéma numérique ainsi obtenu d'autre part.
Les avantages de cette méthode sont:
simplicité de mise en œuvre
efficacité
possibilité de construire des approximations d'ordre élevé
analyse locale(simple) de la précision et de la convergence
les autres méthodes (Eléments Finis et Volume Finis) peuvent souvent s'interpréter comme
des schémas différences finies dans le cas de maillage régulier.
Intérêts des méthodes de différences finies
mais elle possède un certain nombre de limitations
domaine de calcul simple (maillage régulier)
transformation géométrique possible, mais plus complexe
traitement des conditions aux limites
Méthodologie
On couvre d'abord le domaine physique avec une grille régulière (maillage structuré, grille
structurée).
f(x)
2D 1D
Méthodologie
L'équation différentielle est approximée en chaque nœud du maillage en remplaçant les
dérivées partielles par des différences divisées.
f(x)
Le résultat est une expression qui contient des valeurs ponctuelles de la fonction en un
nombre fini de points discrets.
Classification des équations aux dérivées partielles
Il est très recommandable d'examiner un problème d’après ces points suivants:
– L’ordre du problème
– Problème linéaire ou non-linéaire
– Le nombre de variables indépendantes
– Type du problème (elliptique, parabolique, hyperbolique)
Classification des équations aux dérivées partielles
Classification des équations aux dérivées partielles
Classification des équations aux dérivées partielles
EDP linéaire du second ordre à 2 variables indépendantes
Classification des équations aux dérivées partielles
EDP linéaire du second ordre à 2 variables indépendantes
Classification des équations aux dérivées partielles
EDP elliptique
Classification des équations aux dérivées partielles
EDP parabolique
Classification des équations aux dérivées partielles
EDP hyperbolique
Principe de la MDF
La méthode consiste à remplacer les dérivées partielles par des déférences divisées ou
combinaisons de valeurs ponctuelles de la fonction en un nombre fini de points discrets ou
nœuds du maillage.
se déduit directement de la définition de la dérivée.
f(x)
f(x) f(xi+1)
f(xi-1) f(xi)
f\(xi)
Dx Dx x
xi-1 xi xi+1
Principe de la MDF
f ( x dx) f ( x)
x f
dx
f ( x dx) f ( x)
x f
dx
f ( x dx) f ( x)
x f
dx
Example
Résoudre l’ EDO du 2ème ordre suivant
avec les conditions initiales
y(0) = 5
y(10) = -1.5, et Dx = 2
En utilisant les équations des différences divisées:
Example
Substituer dans le ED original pour d2y / dx2 et dy/dx
avec Dx =2:
0 2 4 6 8 10
x
Dx =2
Example
Example
Ou sous forme matriciel:
Example
1-
-1 -0.5 0 0.5 1
x
Dx =0.5
Example
2- 1
Example
3-
Example
4-
Example
5- Solution
L‘erreur de troncature et l‘ordre de l’approximation
Approximation au moyen des développements en série de Taylor
La méthode des déférences finies consiste à approximer les dérivées des équations de la
physique au moyen des développements de Taylor et se déduit directement de la définition de
la dérivée. Elle est due aux travaux de plusieurs mathématiciens du 18éme siècle (Euler,
Taylor, Leibniz...).
Approximation au moyen des développements en série de Taylor
Approximation au moyen des développements en série de Taylor
Approximation au moyen des développements en série de Taylor
Approximation au moyen des développements en série de Taylor
Approximation au moyen des développements en série de Taylor
Approximation au moyen des développements en série de Taylor
Exercice
démontrer que la dérivée première de f en (xi) de l’ordre 2 , 3 et 4 est donnée par:
