NOM PRENOM MATRICULE
LAMARI Omar 202031062952
HOCEINE Chaima 202031088115
COMPTE RENDU DU TP 3 : RESOLUTION DES EQUATIONS AUX
DERIVEES PARTIELLES
I. Le but du TP
L’objectif de ce TP est de résoudre les équations aux dérivées partielles du 2éme ordre linéaire
en appliquant les méthodes numériques étudiées au cours. Les méthodes de résolution à
programmées dans ce TP sont la méthode des différences finies.
II. Différences finies pour les équations elliptiques
II.1. Conditions aux limites Dirichlet
On applique sur une plaque les températures suivantes : coté supérieur100°C, côté droit 50°C,
côté gauche 75°C, coté inférieur 0°C.
La formule utilisée dans la méthode des différences finis pour les EDPs elliptiques :
T i+1 , j +T i−1 , j +T i , j +1+ T i , j−1
T i,j=
4
La matrice de température avec itr=46
Tpr =
75.0000 100.0000 100.0000 100.0000 50.0000
75.0000 78.5714 76.1161 69.6429 50.0000
75.0000 63.1696 56.2500 52.4554 50.0000
75.0000 42.8571 33.2589 33.9286 50.0000
0 0 0 0 0
Remarque
Le nombre de maillage n influence la précision de la solution et la géométrie obtenue.
II.2. Conditions aux limites Neumann ex : un coté isolé
T i+1 ,0 +T i−1 ,0 +2 T i , 1
T i ,0=
4
La matrice de température si le bord inférieur est isolé
Tpr =
75.0000 100.0000 100.0000 100.0000 50.0000
75.0000 83.4109 82.6286 74.2614 50.0000
75.0000 76.0151 72.8420 64.4172 50.0000
75.0000 72.8074 68.3073 60.5652 50.0000
75.0000 71.9073 67.0145 59.5362 50.0000
Ex : deux cotés isolés
Script
clc
clear all
x=10;
y=10;
n=5;
T=zeros(n,n);
T(end,:)=150;
T(:,1)=75;
tol=1e-5;
iter=0;
nr=inf;
while nr>tol
Tpr=T
for j=2:n-1
T(1,j)=0.25*(T(1,j+1)+T(1,j-1)+2*T(2,j));
T(1,n)=0.25*(2*T(1,n-1)+2*T(2,n));
for i=2:n-1
T(i,n)=0.25*(T(i+1,n)+T(i-1,n)+2*T(i,n-1));
T(i,j)=0.25*(Tpr(i+1,j)+Tpr(i-1,j)+Tpr(i,j+1)+Tpr(i,j-1));
end
end
iter=iter+1
nr=norm(Tpr-T)
end
xvec=linspace(0,x,n);
yvec=linspace(0,y,n);
[yy,xx]=meshgrid(yvec,xvec);
mesh(xx,yy,T);
xlabel('x');
ylabel('y');
zlabel('temperature');
Tpr =
75.0000 90.4412 102.5735 110.0183 112.5000
75.0000 92.0956 104.9172 112.5000 114.9816
75.0000 98.0239 112.5000 120.0827 122.4264
75.0000 112.5000 126.9761 132.9044 134.5588
75.0000 150.0000 150.0000 150.0000 150.0000
Remarque
Si le nombre de points de maillage est élevée, donc on pourra obtenir une solution plus précise,
« cela influe sur la représentation de la géométrie obtenue. Ainsi qu’on pourra observer une
meilleure convergence de la solution.
Dans le cas de bord isolé inférieur : on remarque que les températures dans les différents
points changent par rapport à la première méthode (bord inférieur =0)
III. Différences finies pour les équations paraboliques
Utiliser la méthode des différences finies pour résoudre la distribution de température à t=3s
d'une tige longue et mince d'une longueur de 10 cm et les valeurs suivantes : Ax=1cm et
Δt=0,1s. A t=0, la température de la tige est nulle et les conditions aux limites sont fixées à tout
moment à T(0)=100°C et T(10)=75°C. Notez que la tige est en aluminium. Par conséquent,
k=0,835 cm2/s
La formule utilisée dans la méthode des différences finis pour les EDPs paraboliques :
i =T i + λ ( T i+1 −2T i +T i−1 ) Tels que : λ=k ∆ t /(∆ x)
T l+1 l l l l 2
Remarque
La première figure représente la distribution de température dans la tige de t=0 à t=3s avec un
pas de 0,1s.
La deuxième figure représente la distribution de température dans la tige à t =0s et à t=3.
On remarque que la distribution de température s’équilibre au fil du temps.
IV. Conclusion
La précision de la température est directement proportionnelle à l’augmentation du nombre de
points de maillage (n) dans la méthode des différences finies.
