Mon Devoir Modjo
Mon Devoir Modjo
Mon Devoir Modjo
RESOLUTION MANUELLE ET
NUMERIQUE D’UN PBROBLEME DE
DIFFUSION DE CHALEUR AVEC SOURCE
INTRODUCTION :........................................3
I. PRESENTATION DU PROBLEME :........4
II......................RESOLUTION MANUELLE :
4
III. RESOLUTION NUMERIQUE :...............5
IV.GRAPHE DES TEMPERATURES ET
INTERPRETATIONS :..................................8
CONCLUSION :..........................................10
p. 2
INTRODUCTION :
Ce travail consiste à étudier le transfert de chaleur dans un système.
La conduction est un mode de transfert qui résulte d’un transfert d’énergie
cinétique d’une molécule à une autre molécule adjacente. Pour cela, la chaleur
se déplace toujours du point le plus chaud vers le point le moins chaud. Il sera
donc question pour nous de résoudre un problème de conduction thermique :
manuellement, puis numériquement (par programmation en un langage tel que
Matlab, Fortran, tec.) et enfin comparer les solutions (courbes) obtenues par les
deux méthodes.
p. 3
I. PRESENTATION DU PROBLEME :
Nous avons une barre cylindrique ayant les températures T A et T B (T A >T B) à
chacune de ses extrémitées. La chaleur va donc se propager de A vers B.
Pour notre travail, nous prendrons une source de chaleur Q=20 ISU
p. 4
Pour x = 0, on a : T(0) = T(A) = 310 ainsi B = 310
1
Pour x = 1m, on a T(1) = T(B) = 300 ainsi Q(1)2 +A*1 + 310 = 300 ; ainsi A =
2
−(20+Q)
2
Ainsi la solution analytique ou la solution théorique générale est donnée par :
1 −(20+Q)
T(x) = 2 Q x 2 + x + 310
2
{
T 0−2 T 1 +T 2=Q ∆ x 2
2
T 1−2T 2 +T 3 =Q ∆ x
T 3−2T 4 +T 5=Q ∆ x 2
.
.
T i−1−2T i+ T i +1=Q ∆ x 2
. 1
.
.
2
T n−1−2T n +T n +1=Q ∆ x
[ ][ ]
L’écriture sous forme matricielle est :
2
T1 Q ∆ x −T 0
T2 Q ∆ x2
[ ]
T3 Q ∆ x2
−2 1 0
. .
1 −2 1 1 ⋯ 0
0 1 −2 1 . .
. = .
⋮ ⋱ ⋮
. .
0 ⋯ −2 1
. .
1 −2
. .
T n−1 Q ∆ x2
Tn Q ∆ x 2−T n+1
L = 1; %longueur en mètre
p. 6
Te=310; %temperature en entrée
Ts=300; %température en sortie
I = 1; %intensité du courant
n =interface-1;
pas = L/n;
M = zeros(n-1,n-1);
T= zeros(n-1, 1);
C = zeros(n-1, 1);
c = pas*pas*q;
C = c*ones(n-1, 1) + C;
C(n-1,1)= c -Ts;
C(1,1)= c -Te;
for i = 1:n-1,
M(i,i) = -2;
end
M;
for i=1:n-2,
M(i, i+1) =1 ;
end
M;
for i=1:n-2,
M(i+1, i) =1 ;
end
M;
T = inv(M)*C;
p. 7
Tb = zeros(1, n-1);
for i = 1:n-1;
Tb(1, i) = T(i,1);
end
Tb;
Tbb = [ Te, Tb, Ts];
x1=0:pas:1;
x=0:pas2:1;
figure(1)
plot(x1,Tbb,'r');
xlabel('L''espace')
ylabel('Les Temperatures')
title('Courbe Numerique et Analytique')
hold on
plot(x, f, 'b')
title('Courbe Analytique')
p. 8
Courbes pour nombre de points d’interpolation numérique n=2000, pas
numérique = 0.0005 et pour un pas théorique = 0.2
p. 9
Courbes pour nombre de points d’interpolation numérique n=3500, pas
numérique = 0.00028 et pour un pas théorique = 0.125
p. 10
CONCLUSION :
En définitive, il était question pour nous de résoudre manuellement et
numériquement la formulation mathématique (équation différentielle) traduisant
un problème de transfert de chaleur par conduction, à travers une barre soumise
à une source de chaleur. Au vu des résultats obtenus, nous constatons que plus le
nombre de points du maillage est grand (c’est-à-dire que le pas du maillage est
petit), plus on minimise l’erreur ; et par conséquent la solution numérique se
rapproche de plus en plus à la solution réelle.
p. 11