Résolution de l’Equation de la Chaleur par la Méthode des Différences Finies.
Encadré par : Prof. Hamid Elouardi Préparé par : Elboutaybi Sara Bourras Ismail
Rapport du Mini-Projet de l’Aanalyse Numérique II .
Année Universitaire Universitaire:: 1431/1432 2010/2011
Génie Informatique. Informatique.
Table des matières Remerciement......................... Remerciement... ............................................. ............................................. ............................................ ...................................... ................ 3 I.
Introduction ..................... ............................................ ............................................. ............................................ ............................................. ....................... 4
II.
Méthode des différences finies ....................... ............................................. ............................................. ...................................... ............... 4
a. Approximation des dérivées par la formule de Taylor ........................................................ ..............................................................5 ......5 b.Maillage .................... ........................................... ............................................. ............................................ ............................................ ............................................. ....................... 4 c.Schéma Numérique : ............................................................... ............................................................................................................................. ..............................................................6 d.Condition aux limites et conditions initiales : ..................................................................................5 ..................................................................................5 III.
Application.......................... Application.... ............................................. .............................................. ............................................. ......................................... ................... 7
a.Etude Numérique ..................................................................... ......................................................................................................................... .................................................... b.SchémaExplicite............................................................ ............................................................................................................................... ........................................................................... ........ c. Schéma Implicite .................................................................. ................................................................................................................................. ............................................................... d.Schéma de CranckNicolson ................................................................................................................. ................................................................................................................. IV.
Programmation Programma tion ................... ......................................... ............................................ ............................................. ......................................... ..................10
a.Script pour le Schéma Explicite ........................................................... .................................................................................................... ......................................... b.Script pour le Schéma Implicite .................................................................................................... .................................................................................................... c.Script pour le Schéma de CranckNicolson ..................................................................... ..................................................................................... ................ V.
Cas de Calcul .................... ........................................... ............................................. ............................................ ............................................ ......................20
a.Solution exacte ............................................................... .................................................................................................................... ..................................................... b.Solution approchée par le schéma explicite ................................................................... ........................................................................ ..... c.Solution approchée approc hée par le schéma implicite ................................................................... ........................................................................ ..... d.Solution approchée par le schéma de CranckNicolson ....................................................... e.Conclusion sur la stabilité et la convergence convergence de ces méthodes .........................................
VI.
Conclusion ................... .......................................... ............................................. ............................................ ............................................. .......................... ...32
VII.
Références ................... .......................................... ............................................. ............................................ ............................................. .......................... ...33
Remarque : Nous avons choisi l’outil matlab pour la programmation et la visualisation des solutions.
2
Table des matières Remerciement......................... Remerciement... ............................................. ............................................. ............................................ ...................................... ................ 3 I.
Introduction ..................... ............................................ ............................................. ............................................ ............................................. ....................... 4
II.
Méthode des différences finies ....................... ............................................. ............................................. ...................................... ............... 4
a. Approximation des dérivées par la formule de Taylor ........................................................ ..............................................................5 ......5 b.Maillage .................... ........................................... ............................................. ............................................ ............................................ ............................................. ....................... 4 c.Schéma Numérique : ............................................................... ............................................................................................................................. ..............................................................6 d.Condition aux limites et conditions initiales : ..................................................................................5 ..................................................................................5 III.
Application.......................... Application.... ............................................. .............................................. ............................................. ......................................... ................... 7
a.Etude Numérique ..................................................................... ......................................................................................................................... .................................................... b.SchémaExplicite............................................................ ............................................................................................................................... ........................................................................... ........ c. Schéma Implicite .................................................................. ................................................................................................................................. ............................................................... d.Schéma de CranckNicolson ................................................................................................................. ................................................................................................................. IV.
Programmation Programma tion ................... ......................................... ............................................ ............................................. ......................................... ..................10
a.Script pour le Schéma Explicite ........................................................... .................................................................................................... ......................................... b.Script pour le Schéma Implicite .................................................................................................... .................................................................................................... c.Script pour le Schéma de CranckNicolson ..................................................................... ..................................................................................... ................ V.
Cas de Calcul .................... ........................................... ............................................. ............................................ ............................................ ......................20
a.Solution exacte ............................................................... .................................................................................................................... ..................................................... b.Solution approchée par le schéma explicite ................................................................... ........................................................................ ..... c.Solution approchée approc hée par le schéma implicite ................................................................... ........................................................................ ..... d.Solution approchée par le schéma de CranckNicolson ....................................................... e.Conclusion sur la stabilité et la convergence convergence de ces méthodes .........................................
VI.
Conclusion ................... .......................................... ............................................. ............................................ ............................................. .......................... ...32
VII.
Références ................... .......................................... ............................................. ............................................ ............................................. .......................... ...33
Remarque : Nous avons choisi l’outil matlab pour la programmation et la visualisation des solutions.
2
Remerciements vant d’entamer ce rapport, nous profitons de l’occasion pour A remercier chaleureusement notre cher professeur pour avoir crée cette occasion (le mini-projet) et nous permettre ainsi de voir nos acquis purement théoriques rencontrer le monde réel par le biais de la programmation. programmation.
Ce travail est le fruit de vos efforts et de votre générosité qui nous a surveillés partout et presque chaque semaine via l’adresse électronique.
I.
Introduction :
L’analyse numérique est une discipline des mathématiques. Elle s’intéresse tant aux fondements théoriques qu’à la mise en pratique des méthodes permettant de résoudre, par
des calculs purement numériques, des problèmes d’analyse mathématique.
algorithmes permettant algorithmes permettant de résoudre résoudre les problèmes de mathématiques continues (distinguées des mathématiques discrètes). discrètes) . Cela Plus formellement, l’analyse numérique est l’étude des
3
signifie qu’elle s’occupe principalement de répondre numériquement à des questions à variable réelle ou complexe comme l’algèbre linéaire numérique sur les champs réels ou complexes, la recherche de solution numérique d’ équations différentielles et d’autres problèmes liés survenant dans les sciences physiques et l’ingénierie.
Dans le domaine de l'analyse numérique, on peut être amené à rechercher la solution d'une équation aux dérivées partielles. Parmi les méthodes de résolutions couramment pratiquées, la méthode des différences finies est la plus facile d'accès, puisqu'elle repose sur deux notions : la discrétisation des opérateurs de dérivation/différentiation (assez intuitive) d'une part, et la convergence du schéma numérique ainsi obtenu d'autre part. En mathématiques et en physique théorique, l'équation de la chaleur est une équation aux dérivées partielles parabolique, introduite initialement en 1811 par Fourier pour décrire le phénomène physique de conduction thermique.
II.
Méthode des differences finies :
Parmi les méthodes de résolution, la méthode des différences finies, qui repose sur deux notions : la discrétisation des opérateurs de dérivation/différentiation par différences finies d'une part, et la convergence du schéma numérique ainsi obtenu d'autre part. En effet Un problème aux dérivées partielles nécessite la donnée de : D’un domaine D’une équation aux dérivées partielles
De conditions aux limites
De conditions initiales
a. Approximation des dérivées par la formule de Taylor :
Grâce aux formules de Taylor, on définit la discrétisation des opérateurs différentiels (dérivées premières, secondes, etc. partielles ou non). La formulation de Taylor-Young est préférable dans son utilisation simple, la formulation de Taylor avec reste intégral de Laplace permet de mesurer les erreurs
Et :
4
Ou toutes les applications convergent vers 0 avec h. Alors :
Et en sommant les développements pour x-h et x+h l'on obtient:
On obtient respectivement des approximations de 1 er ordre et 2nd ordre en h.
Alors on a :
. . On peut aussi montrer de la même façon que : Et
b. Maillage :
Un maillage est un ensemble de points du domaine de définition sur lequel on va appliquer la méthode des différences finies. Pour une application définie sur un segment de , on ajoutera en général les deux extrémités du segment ; pour un maillage en dimension supérieure, on sera amené à choisir, éventuellement, des points des contours du domaine de définition. On appelle le pas du maillage la distance entre deux points successifs du maillage voisins. En dimension 1, cela se simplifie en différence des abscisses. Ce pas n'est pas nécessairement constant, il peut même être judicieux de ne pas le fixer comme tel. Le pas (global) de l'approximation peut être défini comme le plus grand pas du maillage. Ainsi, si ce pas global tend vers 0, cela veut dire que la répartition des points du maillage dans l'intervalle choisi tend à se faire sur tout le domaine d'étude par densité.
5
Exemple :
Un intervalle de validité [0,1] on utilisera (N + 1) points, par exemple {0, h, 2h,...,
N*h=1} pour un pas constante
c. schema numérique :
Écrire un schéma numérique de résolution de l'équation différentielle initiale signifie :
substituer les formulations des dérivées/différentielles obtenues par approximation aux opérateurs eux-mêmes sur tous les points du maillage. réorganiser les équations pour faire apparaître un schéma explicite (ex : les valeurs à la date t+1 données en fonction des valeurs des dates 0 à t) ou implicite (une équation lie les valeurs passées, présentes et futures sans qu'on arrive à exprimer ces dernières seules).
Dans un cadre de modélisation classique d'opérateurs linéaires dans des équations différentielles linéaires, on aboutit à un système d'équations linéaires de dimension égale au nombre de nœuds du maillage (en fait un peu moin s, du fait des données initiales, par exemple).
Résoudre le schéma numérique signifie simplement trouver les valeurs discrètes de la fonction en chaque nœud. Un système issu d'une équation linéaire peut souvent être algébriquement simple à résoudre. Pour simplifier, on peut dire que les schémas explicites engendrent des systèmes d'équation à matrice triangulaire ou trigonalisables, ce qui n'est pas le cas des schémas implicites. d. Condition aux limites et conditions initiales :
Un problème aux limites est une équation aux dérivées partielles munie de conditions aux limites sur la totalité de la frontière du domaine. Un problème de Cauchy est une équation aux dérivées partielles ou, pour la variable de temps, les conditions « au bord« sont des conditions initiales (et pas finales). On dit que le problème A(u) = f est bien posé si pour toute donnée f ; il admet une solution unique u, et si cette solution u dépend continument de la donnée f conditions nécessaire pour faire du calcul numérique. Les équations de type (I1) sont représentatives de problème de type potentiel qui apparaissent dans des études de régime permanent en électricité (électrostatique ou magnétostatique), mécanique (déformation d’un solide, écoulement) et thermique (répartition
des températures), les conditions aux limites associées sont de type :
6
(a)
Dirichlet :
(b)
Neumann :
(c) Mixte
III.
:
u(s)=u0
Application :
Voici l’équation à résoudre
(E ) :
Avec
données.
{
Il s’agit là des conditions de Dirichlet ;
Voici l’organigramme correspondant à notre résolution :
Début
Discrétisation du domaine de la définition de l’équation.
Conditions initiales
7
Expression de la solution par chaque schéma des MDF. Valorisation de l’erreur.
Itération Max ?
Fin
On note :
Avec
. ; et
Maintenant on passe à établir schéma par schéma. a. Schéma Explicite :
En se basant sur ce qui precède, exprimons les 2 dériveés présentes dans (E ) au point (ih, (j-1) k) :
8
b. Schéma Implicite :
En se basant sur ce qui precède, exprimons la première dérivée à (i h, j k) et la deuxième à (i h, (j+1) k) :
{ c. SchémadeCrank Nicolson :
Le schéma de crank Nicolson est un cas spécifique de la Ө-Méthode où
En effet pour ce schéma Ө = .
.
En se basant sur ce qui precède, exprimons la première dérivée à (i h, j k) ; Pour la deuxième cela diffère :
× (approximation à (i h, j k)) + (1- ) × (approximation à (i h, (j+1) k)).
Après le calcul:
Avec :
IV.
et
1
programmation
Et avec les programme suivant on resoud le probleme numeriquement : 1.1 Programme du schéma Explicite :
9
10
11
1.2 Programme du schéma Implicite :
Donc Posant On a
On a alors
En
12
13
En posant :
[] [
(1 2a )
0
a
a
a
0
a
(1 2a )
(1 2a )
0
a
a
[
a
0
1
B=
[
a
0
0
0
] [] [ ]
0 0
0
1
(1 2a )
] [
2
] 1
1
0
0 1
1
2
]
= -I + aT
On utilise alors la méthode indirecte de résolution des systèmes linéaire : JACOBI, GAUSS SEIDEL, RELAXATION) A×x=b
Schéma implicite à l’aide de Jacobi
14
15
1.3 Programme du On a
*
+
* ( ) ( )+ Et
16
17
(1 a)
[
a 2
0
a
a
2
2
0
a 2
(1 a)
]
(1 a )
[
a
2
0
a
a
2
2
0
a
2
(1 a)
]
18
a
a
0
(1 a)
a
a
a
a
2
2
2
2
(1 a)
0
(1 a)
[ ] [ 0
2
a
2
2
a 2
0
(1 a)
]
On utilise alors la méthode indirecte de résolution des systèmes linéaire : JACOBI, GAUSS SEIDEL, RELAXATION)
.
19
V. Cas de calcul : a. Solution Exacte :
{
Notre équation : (E )
Avec maintenant pour application :
20
{
Pour ce faire, on utilise le cour des équations aux dérivées partielles de 1 ère année :
Faire de sorte de simplifier l’expression de la forme solution en l’exprimant
comme é tan t l e produit de deux fonctions in dé pendantes :
Nous posons :
u(x, t) v(x).w(x), nous obtenons alors
v = k
(E)
v.w’=
‘’.w
k une constante.
Ceci car les expressions dans chaque membre de droite et de gauche ne peuvent être en tout temps égales que si elles sont égales à un constante. D’où :
Dé term iner un e solution v(x) qui sati sfasse les condi ti ons aux f r ontiè res.
L’équation (1) a pour solution :
Ou bien
v(x)=ax+b, si k=0.
√ √ √
v(x)=C1 .
+ C2.
v(x)=A.cos(
si k>0.
)+B.
. k<0
Pour ne pas tomber sur des solutions triviales, on suppose que k<0. On a selon les conditions initiales : U(4,t)=v(4).w(t)=0,
D’où A= 0, et (B=0 ou
u(0,t)=v(0).w(t)=0,
.
),
21
Pour ne pas tomber sur une solution triviale on prend le cas A=0 et Ceci dit :
√ =n.
vn(x)=
D’où
.
, n=1,2,3………
Dé ter miner la f or me soluti on w(t) en consé quence des imposé es sur v(x).
Avec la contrainte k=-(
W n ( t)= C n
)², et avec l’équation (2), on trouve :
n=1,2,3,4……………
Exprimer la solution u(x, t) sous forme de série de Fourrier pour qu’elle
satisfasse àpr é sent àla conditi on initi ale sur t.
On a alors la famille des solutions :
n=1,2,3…..
La solution générale u g(x, t) est la combinaison linéaire de toutes ces u n(x, t) :
∑ ∑ ∫ Fixons les coefficients C n :
u(x,0)=
Avec :
Formules d’Euler pour les fonctions p=2
4 périodique :
=>
avec w* est la fonction de prolongement impair fictif entre -4< x <0, de sorte que
constitue une fonction paire.
C-à-d :
22
Après un calcul avec des formules trigonométriques et vérification des conditions, on arrive à la solution exacte :
***********************************
Pour des raisons de facilité d’utilisation, on a regroupé toutes les scripts des méthodes numériques déjà citées dans un seul programme, permettant à l’utilisateur de choisir la méthode qu’il veut.
Voici le M-file utilisé :programmefinale.m On prend :
On trouve : N =19
, pour visualiser la solution du cas de calcul proposé.
et a= 0,4053.
Dont le code est :
23
24
25
b. Solution approchée par le schéma explicite : Choix : 1 dans le programme prog.m Analyse des figures ci-dessous :
Nous pouvons constater que la solution approchée suit en gros l’allure de la solution exacte mais reste, en détails, un peu loin de lasolution exacte ; la courbe d’erreur nouspermet de voir cela : en effet l’erreur est nulle sur des points en particuliers ce quipermet à lacourbe de lasolution approchée de suivre celui dela solution exacte.
26
Fig. courbe solution exacte et approchée et Erreur. A l’aide du mesh :
c. Solution approchée par le schéma implicite : Choix : 2 dans le programme prog.m Analyse des figures ci-dessous :
La meme chose que la solution approchée par le schéma explicite, la solution approchée par le schéma implicite suit l’allure de la solution exacte.
27
Comme avant avec mesh :
Fig. Courbe de la solution approchée et exacte et Courbe d’Erreur.
Fig. Courbe de la solution approchée schéma implicite par la méthode de Jacobi .
28
L’utilisation de Jacobi montre une courbe d’erreur plus importante, la courbe s’écarte de la courbe de la solution exacte. Alors la double boucle se montre
simple et efficace : d. Solution approchée par la méthode de Crank Nicolson : Choix : 3 dans le programme prog.m Analyse des figures ci-dessous : A l’encontre du schéma implicite et explicite, la méthode de CrankNicolson se dévie
clairement de la solution exacte. En appliquant CrankNicolson par la méthode de Jacobi, la solution ma lgré qu’elle soit loin de la solution exacte, elle reste meilleure que CrankNicolson seul.
A l’aide de mesh :
Fig. Courbe de la solution approchée et exacte et Courbe d’Erreur.
29
Fig. Courbe de la solution approchée Crank Nicolson par la méthode de Jacobi.
e. Conclusion sur la convergence et la stabilité de ces méthodes :
Convergence :
E l’erreur de convergence s’exprime comme suivant :
‖‖
.
Pour le schéma explicite :
l’erreur présente un maximum après tend vers 0 d’une façon remarquable . Pour le schéma implicite :
30
l’erreur présente un maximum après tend vers 0 d’une façon remarquable mais
lente par rapport au schéma explicite.
Pour le schéma de Crank Nicolson :
L’erreur est nulle jusqu’à un certain rang,après il augmente d’une façon
importante.
Stabilité :
Là on puise directement du cours de notre cher prof. L’étude de la stabilité se résume dans l’évaluation du facteur λ tel que :
.=d. .
Pour le schéma explicite : Il faut avoir :
λ
31
On a:λ
λ
.
Schéma stable donc.
Pour le schéma implicite : Il faut avoir :
Schéma stable donc.
et on a :
. λ
Pour le schéma de CrankNicolson : On a
, la méthode est inconditionnellement stable.
Synthèse : Le schéma explicite se montre fort avec la rapidité de convergence de l’erreur vers 0.
Tandis que le schéma de CrankNicolson semble clairement diverger de la solution exacte. Le schéma implicite reste acceptable par rapport au schéma de Crank.
VI.
Conclusion : Ce
projet
illustre
bien
l’importance
des
méthodes
numériques
pour
la
résolution des problèmes mathématique, leur variétés, et permet de constater le peu de différences concernant les solutions proposées par chaque méthode, d’où leur efficacité. D’autre part, il nous a été très utile de travailler sur ce projet, sachant que d’une part on a pu mieux concevoir l’idée de devoir résoudre un tel problème mathématique, et d’un autre côté se familiariser davantage avec un outil important pour un élève in génieur qu’est le Matlab, et ainsi reconnaitre son
utilité.
32