République Algérienne Démocratique et
Populaire
Université Echahid Hamma Lakhdar – El Oued
Faculté des sciences et de latechnologie
2éme Master Energétique
Menu projet
CODE MATLAB
Présenté par : PROF :
S.MEGDOUD
HENKA .Ali
AZEB ABDE ALLAH .YOUSEF
OUASSA .MOUNIER
Année universitaire: 2019/2020
MATLAB
Serie de TD N : 01
Exercice N :01
Editor
clear all
clc
TCe=100; % la temperature de entrie du fluide chaud[k]
TCs=45; % la temperature de sortie du fluide chaud[k]
TFe=15; % la temperature de entrie du fluide froid[k]
TFs=30; % la temperature de sortie du fluide froid[k]
L=10; % longuer
x=0:0.1:L;
%A contre courant
'delta Temeperature moyne contre courant [c]'
deltaTLMcc=((TCs-TFe)-(TCe-TFs))/log((TCs-TFe)/(TCe-TFs))
%A co courant
'delta Temeperature moyne co courant[c]'
deltaTLMcp=((TCs-TFs)-(TCe-TFe))/log((TCs-TFs)/(TCe-TFe))
%plot A contre courant
yTC=((TCe-TCs)/(0-L))*x+TCe;
yTF=((TFs-TFe)/(0-L))*x+TFs;
figure(1);
plot(x,yTC,x,yTF,'linewidth',2)
ylabel('Temeperature(c°)');
xlabel('L(m)');
title('A contre courant');
legend('Tc','Tf');
grid on
%plot A co courant
zTC=((TCe-TCs)/(0-L))*x+TCe;
zTF=((TFe-TFs)/(0-L))*x+TFe;
figure(2);
plot(x,zTC,x,zTF,'linewidth',2);
ylabel('Temeperature(c°)');
xlabel('L(m)');
title('A co courant');
legend('Tc','Tf');
grid on
1
commed window
delta Temeperature moyne contre courant [c]
deltaTLMcc =
47.2089
delta Temeperature moyne co courant[c]
deltaTLMcp =
40.3551
2
Exercice N :02
Editor
clear all
clc
TCe=110; % la temperature de entrie du fluide chaud[k]
TCs=30; % la temperature de sortie du fluide chaud[k]
TFe=12; % la temperature de entrie du fluide froid[k]
DebitC=5000/3600; % debit du fluide chaud[g/s]
DebitF=12000/3600; % debit du fluide froid[g/s]
CPc=2100; % les caracteristique thermo physiques du fluide chaud[J/kg*k]
CPf=4180; % les caracteristique thermo physiques du fluide froid[J/kg*k]
r=0.5; % rayonement
L=2; % longuer
x=0:0.001:L;
'surfaces d echange [m2]'
S=2*pi*L*r
'Temeperature TFs[c]'
Q=DebitC*CPc*(TCe-TCs)
TFs=(Q/(DebitF*CPf))+TFe
%A contre courant
'la puissance de l ehangeur contre courant'
'delta Temeperature moyne contre courant[c]'
deltaTLMcc=((TCs-TFe)-(TCe-TFs))/log((TCs-TFe)/(TCe-TFs))
Ucc=Q/(S*deltaTLMcc)
%A co courant
'delta Temeperature moyne co courant[c]'
deltaTLMcp=((TCs-TFs)-(TCe-TFe))/log((TCs-TFs)/(TCe-TFe))
Ucp=Q/(S*deltaTLMcp)
%plot A contre courant
yTC=((TCe-TCs)/(0-L))*x+TCe;
yTF=((TFs-TFe)/(0-L))*x+TFs;
figure(1);
plot(x,yTC,x,yTF,'linewidth',2)
ylabel('Temeperature(c°)');
xlabel('L(m)');
title('A contre courant');
legend('Tc','Tf');
grid on
%plot A co courant
zTC=((TCe-TCs)/(0-L))*x+TCe;
zTF=((TFe-TFs)/(0-L))*x+TFe;
figure(2);
plot(x,zTC,x,zTF,'linewidth',2)
ylabel('Temeperature(c°)');
xlabel('L(m)');
title('A co courant');
legend('Tc','Tf');
grid on
3
commed window
surfaces d echange [m2]
S=
6.2832
Temeperature TFs[c]
Q=
2.3333e+05
TFs =
28.7464
la puissance de l ehangeur contre courant
delta Temeperature moyne contre courant[c]
deltaTLMcc =
41.9675
Ucc =
884.8783
delta Temeperature moyne co courant[c]
deltaTLMcp =
22.1949
Ucp =
1.6732e+03
4
5
Exercice N :03
Editor
clear all
clc
TCe=180+273; % la temperature de entrie du fluide chaud[k]
TFe=40+273; % la temperature de sortie du fluide chaud[k]
TFs=60+273; % la temperature de sortie du fluide froid[k]
DebitC=1000*1000/3600; % debit du fluide chaud[g/s]
DebitF=2000*1000/3600; % debit du fluide froid[g/s]
CPc=4315; % les caracteristique thermo physiques du fluide chaud[J/kg*k]
CPf=4315; % les caracteristique thermo physiques du fluide froid[J/kg*k]
Re=10000; % nobre de raynolds
D=0.02; % diametre intereur[m]
K=450; % le coeffcient d echange global[w/m2*k]
u=0.00019; % la viscosité dynamique[kg/m*s]
ron=920; % la masse volumique[kg/m3]
'la puissance echangee[J]'
Q=DebitF*CPf*(TFs-TFe)
'la temperature de sortie du fluide chaud[k]'
TCs=TCe-(Q/(DebitC*CPc))
'delta Temeperature moyne contre courant[k]'
deltaTLMcc=((TCs-TFe)-(TCe-TFs))/log((TCs-TFe)/(TCe-TFs))
'la suface d echenge neeessaire[m2]'
S=Q/(K*deltaTLMcc)
'la vitesse dans les tubes[m/s]'
V=(Re*u)/(ron*D)
'la section totale des tubes[m2]'
s=pi*D*1
6
commed window
la puissance echangee[J]
Q=
4.7944e+07
la temperature de sortie du fluide chaud[k]
TCs =
413
delta Temeperature moyne contre courant[k]
deltaTLMcc =
109.6963
la suface d echenge neeessaire[m2]
S=
971.2562
la vitesse dans les tubes[m/s]
V=
0.1033
la section totale des tubes[m2]
s=
0.0628
