Distribution de la température à
l’intérieur d’une plaque infini
Présenté par : Encadré par :
CHAKKOUR Ouassima BOUARDI Abdelmjid
EL KAITONI Souhail
2015/2016
Sommaire :
I. Introduction ……………………………………….03
II. Caractéristiques thermiques d’Aluminium pur…04
1. La conductivité ………………………………04
2. La masse volumique …………………………04
3. La chaleur thermique ……………………….04
4. La diffusivité thermique ……………………04
2
I. Introduction :
La Détermination de la distribution de température à l’intérieur d’une plaque infini
en aluminium et d’une épaisseur 4 e.
On a utilisé pour ce projet l’aluminium pur :
L'aluminium est un métal mou, léger, mais résistant avec un aspect argent-gris mat,
dû à une couche mince d'oxydation qui se forme rapidement quand on l'expose à l'air
et qui empêche la corrosion de progresser. À la différence de la plupart des métaux, il
est utilisable même s'il est oxydé en surface.
L'aluminium a une densité environ trois fois plus faible que celle de l'acier ou
du cuivre ; il est malléable, ductile et facilement usiné et moulé. Il possède une
excellente résistance à la corrosion et une grande longévité. Il est également non
magnétique et ne provoque pas d'étincelles. C'est le deuxième métal le plus malléable
et le sixième le plus ductile. En solution, l'aluminium se trouve le plus généralement
sous la forme d'ions Al3+.
3
II- Caractéristiques thermiques d’Aluminium pur:
1. La conductivité :
: La conductivité thermique d’aluminium pur.
λ= 237 watt. .
2. La masse volumique :
ρ= 2.7*10^3 Kg. .
3. La chaleur thermique :
Cp=888 j. .
4. La diffusivité thermique :
a= =98,8*10^ (-6) .
L’épaisseur de la plaque :
4 e=0,4m
4
L’équation de transfert de chaleur :
= (a)
Les conditions aux limites (Neumann) :
(b)
x=0 (c)
= 0 (d)
s On effectue le changement de variable suivant :
=
L’équation de la chaleur devient :
5
= (a)
Donc les conditions aux limites deviennent :
(x,t)= (b)
x=0 (c)
0 (d)
Pour la résolution on a choisi la méthode de la transformée de Laplace de
l’équation de chaleur on obtient:
= Ɵ → - Ɵ =0 et q²=
Ɵ(x,p)=Acosh(qx)+B sinh(qx).
D’après les conditions aux limites on aura :
(qx)+ Bqcosh (qx)=0
X=0 donc Bq=0 alors B=0
Ɵ (x,p)=A cosh (qx)
Déterminant A :
On a x=0
= Φ0/p
Φ0/p
Donc A= Φ0/
6
Ɵ(x,p)= (Φ0 /λpq) *(cosh(qx)/sinh(qL)) avec q=
-Première méthode :
La température T(x,t) peut s’en déduire en appliquant la formule de
« stehfest » pour trouver la transformée de Laplace inverse
T(x,t)= )
Un nombre de terme N=10 est suffisant pour obtenir une précision
satisfaisant :
T(x,t)=
Donc :
-La température d’une plaque infinie en Aluminium
Résultat 1 :
Le transfert de chaleur en x=0,1cm de l’épaisseur
7
Résultat 2 :
Le transfert de chaleur en x=0,25 cm de l’épaisseur
Résultat 3 :
Le transfert de chaleur en x= 0,3cm de l’épaisseur
8
On a utilisé le programme de MATLAB selon l’Algorithme
suivant :
clear ;
format long e ;
v(1)=1./12 ;
v(2)=-385./12;
v(3)=1279;
v(4)=-46871./3;
v(5)=505465./6;
v(6)=-473915./2;
v(7)=1127735./3;
v(8)=-1020215/3;
v(9)=328125./2;
v(10)=-65625./2;
Tfini=10000000;
T0=40;
Ti=20;
DT=Ti-T0;
landa=237;
phi0=0.034;
a=9,89.e-6;
e=0.4;
x=0.1;
x=0.25;
x=0.3;
for it=1:201,
tt=it*Tfini/200;
y(it)=tt;
ln2_on_t=log(2.0)/tt;
sum=0.0;
for n=1:10,
p=n*ln2_on_t;
q=sqrt(p/a);
teta=(phi0*cosh(q*x))/(p*landa*q*sinh(q*e));
sum=sum+v(n)*teta ;
end
T(it)=Ti+sum*ln2_on_t;
end
plot(y,T,'b-')
ylabel('T')
xlabel('t')
9
-Deuxième méthode :
T=Ti+2*Φ0 *
clear;
format long e;
% Teméprature intiale
Ti=20;
% Température imposée aux deux extrémités du domaine
T0=40;
% Diffusivité thermique
a=9,89.e-6;
% Longueur du domaine (épaisseur du mur)
e=0.4;
%Durée de calcul
Tc=10000000;
%Nombre de points temporels
m=101;
%Pas de temps
Dt=Tc/(m-1);
%Nombre de points d'espace
n=21;
%Pas d'espace
Dx=e/(n-1);
Td=T0-Ti;
% Nombre de termes de la série tronquée
NS=101;
% Boucle sur le temps
%Boucle sur l'espace
x=0.1;
x=0.25;
x=0.3;
phi=0.30;
lambda=237;
for i=1:m+1,
t=1.e-6+(i-1)*Tc/m;
TT=Ti;
% Boucle sur les termes de la série tronquées
for k=1:NS,
u1=((2*k+1)*e-x)/(2*sqrt(a*t));
u2=((2*k+1)*e+x)/(2*sqrt(a*t));
c1=exp(-u1^2)/sqrt(pi)-u1*(1-erf(u1));
c2=exp(-u2^2)/sqrt(pi)-u2*(1-erf(u2));
TT=TT+(2*phi*sqrt(a*t)/lambda)*(c1+c2);
end
T(i)=TT;
y(i)=t;
end
%Y=T(:,i);
%plot(X,Y');
%pause(0.1)
plot(y,T,'r-');
xlabel('Temps t (s)');
ylabel('T(t)');
10
Résultat 1 :
Le transfert de chaleur en x=0,1cm de l’épaisseur
Résultat 2 :
Le transfert de chaleur en x=0,25cm de l’épaisseur
11
Résultat 3 :
Le transfert de chaleur en x=0,3cm de l’épaisseur
12