Hafsa - Mongach GA TP Thermo

TP n°1 :
Prise en main de Pdetool de Matlab
2-
TP N° 2 : Les ailettes
I. Répartition des températures le long des ailettes triangulaire et rectangulaires.
Considérons deux ailettes de formes distinctes, une triangulaire et l’autre rectangulaire. Les deux sont faites du
même matériau de conductivité thermique (lamda) k=80,2 W/m.K et ayant chacune 20 mm d'épaisseur et 80
mm de longueur. L’installation des deux ailettes sur un mur refroidi à l’air a pour but d’augmenter la chaleur
évacuée par ce mur. La température du mur est Tb = 115°C et le coefficient de convection entre la surface du
solide et l’air (T=15°C) est h = 100 W/m2 .K
1. En recourant au logiciel Matlab , on se propose d'étudier seulement la moitié de l’ailette.

Justifier ce choix :
Car l’axe de symétrie est un axe adiabatique, alors il va nous faciliter l’etude de l’ailette par étudier juste la
moitié.
2. Ailette de section triangulaire

En supposant un transfert de chaleur en deux dimensions. Déterminer, par simulation numérique, les
températures à différentes positions 0  x  L, y  0. Tracer sur le graphe de la figure 2 la répartition des
températures le long d'ailette.
X 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07

T(x, y=0) 115 107 100 94 88 82 76 70
3. Ailette de section rectangulaire

En maintenant les mêmes hypothèses que dans l’ailette triangulaire, déterminer par simulation puis tracer
sur même graphe (figure 2) avec un motif différent la répartition des températures le long de l'ailette
rectangulaire. Comparer le comportement des deux ailettes.
X 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07

T(x, y=0) 115 107 100 94 90 86 83 81
II. Confrontation entre les résultats numériques et les résultats analytiques.
Nous proposons ici d'étudier par deux méthodes le comportement thermique d'ailette de section
rectangulaire. Dans cette étude, nous supposons que le flux à la surface extrême x=L est négligeable et égale à
zéro.
1. Dans ses conditions, montrer que la solution analytique du problème s'écrit :
Avec : m=h.P/(.S)
On donne : S=1*e et P=2*(e+1)
A partir de la solution analytique (Eq. 1), calculer les températures T(x).
2. Tracer sur le même graphe ci-dessous (T(x) de la solution analytique tracée auparavant sur
le graphe) la répartition des températures le long d'ailette de la solution numérique.
Comparer les résultats des deux approches et conclure.
3. Déterminer par simulation la densité de flux (0,y) à la base de l’ailette.
Pour afficher les densités des flux, aller dans le menu Plot et sélectionner Paramètres. Dans
la boite de dialogue plot sélectionner apparait ; à la place de la température sélectionner
heat flux. Pour plus de précision, il est préférable de faire un zoom au voisinage de x=0. Pour
cela, sélectionnez le menu Options puis l'item Axes Limits.... Entrez les bornes du domaine
du zoom, par exemple [0, 0.001] en x.
Conclusion :
On peut conclure d’après le tableau ci-dessus que l’effet des fluides influence sur l’efficacité
de l’ailette et laisser l’ailette d’avoir une même efficacité par rapport a n’importe quel effet
de fluide .
TP3 Conduction bidimensionnelle :
1- Approche analytique :
a- Les conditions aux limites des domaines y=0 et y=1 :
On suppose que le flux est nul suivant y=0 et y=1 pour constater un problème de mur (Méthode
Dirichlet)
(Car on a épaisseur > L)
b- Les équations différentielles de la chaleur des deux milieux 1 et 2 :
- Milieu R1 :
𝑇₁(𝑥) = −𝑄𝑣 2𝜆₁ 𝑥 2 + 𝑎𝑥 + 𝑏 = −46.15𝑥 2 + 202.7𝑥 + 500
𝜑₁(𝑥) = Q𝑣 𝑥 − aλ₁ = 60x − 131.75
- Milieu R2 :
𝑇₂(𝑥) = 𝐶₁𝑥 + 𝐶₂ = −241.21𝑥 + 1416.36
𝜑₂ = −C₁λ₂ = 48.24 W ⁄ cm²
x 0 1 2 3 4 5 6
T(x) 500.02 717,61 834,53 858,49 794,71 637,30 391,55
𝜑 (x)
x 0 1 2 3 4 5 6
-131.75 -71.74 -11.76 48.25 48.04 48.24 48.24
C - détermination de T(x=e1) et T(x=e₁+e₂) :
On a e₁= 4 cm T(x=e₁) = -46.15x²+202.7x+500
=-46.15*3² + 202.7*3 +500 = 572˚C
Et pour e₁+e₂= 6 cm
T(x=e₁+e₂) =-46.15x²+202.7x+500
=54.8 ˚C
D -détermination de la position et la valeur de T max :
-la position :
On a T(x) = -46.15x²+202.7x+500 Et (𝑥) 𝑑𝑥 = 0
Donc 0 = 2*(-46.15) x + 202.7 92 .3 x = 202.7
Donc : X = 2.20
CYCLE PAD :
TP 1 : Machine frigorifique à deux étages :
Le logiciel CYCLEPAD permet de simuler des cycles thermodynamiques ainsi que des transformations
particulières rencontrées dans l'industrie. On peut à l'aide de ce logiciel réaliser assez rapidement des
calculs concernant des installations industrielles pouvant être assez complexes.
I. Diagramme :
1. Flow fraction=0.18
2. Q dot=33.74( de l evaporateur HTR2)
3. COPr=3.73( cycle)
4. Tools-sensitivity-lier la pression
5.
COP (70%) COP (80%) COP (85%) COP (90%) COP (100%)
R-12 2,31 2,67 2,84 3,02 3,38
R-22 2,47 2,86 3,06 3,25 3,64
R-134a 2,55 2,95 3,15 3,34 3,73
Azote 1,73 2,07 2,22 2,36 2,66

Ammoniac 2,73 3,17 3,40 3,62 4,06
Méthane 1,89 2,22 1,89 2,52 2,84
TP N°2 Le cycle de Rankine à régénération idéal

Soit une centrale thermique fonctionnant selon le cycle à régénération idéal avec un seul réchauffeur
à mélange (voir la figure). La vapeur d'eau est admise dans la turbine à 15 MPa et à 600°C, et elle est
condensée dans le condenseur à 10 kPa. Une partie de la vapeur est soutirée à 1,2 MPa pour être
détournée vers le réchauffeur à mélange. L'état du fluide à l'entrée des pompes est un état saturé
avec x=0 (liquide saturé). La compression et la détente sont isentropiques. Les phases de
condensation, de chauffage et de mixage ont lieu à pression constante.
Déterminez :
- La fraction de vapeur soutirée de la turbine : 0.2281
- Le rendement thermique du cycle : 46%
Avec ŋ=- ẇ /Qm Ẇ=1283 Qm=2767
- La température d'entrée de l'eau dans la chaudière : 190 C°
Etudier l’effet de la pression du soutirage sur les paramètres suivants :

- le rendement thermique du cycle : 46.4%
- le travail net : 1283 KW
- La puissance thermique consommée : 1300 KW
- La fraction de vapeur soutirée de la turbine. 0.77
- La température d'entrée de l'eau dans la chaudière. 190

Déduire la pression optimale de soutirage, c'est-à-dire la pression de soutirage qui maximise le
rendement thermique du cycle.
Conclusion:
On peut deduire que le travail net et la puissance thermique diminue avec
l’augmentation de la pression par contre la fraction de vapeur augmente avec
L’augmentation de la pression.
TP N°3 Le cycle à resurchauffe et à régénération
idéal
Le cycle a resurchauffe et a régénération idéal
La Fraction de vapeur soutirée de la turbine du second étage est : 0.1314

Le travail produit dans l’étage à haute pression de la turbine : 428.9 KW
Le travail produit dans l’étage à basse pression de la turbine est : 471.7 kw

La chaleur fournie a la chaudière et a la resurchauffe ( Q-dot in) 2.923 kw
La chaleur évacuée : (Q dot out) -1.491 kW

Le rendement thermique : (état thermal) 48.98 %
etage 1
Les valeurs ‘’ from’’ et ‘’ to’’ doivent etre de 500 et 1500 car S8<S11<S15
La pression max est presque 2512 KPa
etage 2
Conclusion :
Avant l’optimisation des deux pressions on a un rendement qui est faible et pour améliorer notre
machine on besoin d’améliorer le rendement de la machine d’abord on fait l’optimisation des deux
pressions ce qui nous permettant d’augmenter le rendement thermique de 48.98 % jusqu'à 49.20% .

Hafsa - Mongach GA TP Thermo

Transféré par

Droits d'auteur :

Formats disponibles

Hafsa - Mongach GA TP Thermo

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Hafsa - Mongach GA TP Thermo

Transféré par

Droits d'auteur :

Formats disponibles

TP n°1 :

Prise en main de Pdetool de Matlab

1. En recourant au logiciel Matlab , on se propose d'étudier seulement la moitié de l’ailette.

2. Ailette de section triangulaire

X 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07

3. Ailette de section rectangulaire

X 0 0.01 0.02 0.03 0.04 0.05 0.06 0.07

1. Dans ses conditions, montrer que la solution analytique du problème s'écrit :

On donne : S=1*e et P=2*(e+1)

A partir de la solution analytique (Eq. 1), calculer les températures T(x).

a- Les conditions aux limites des domaines y=0 et y=1 :

(Car on a épaisseur > L)

b- Les équations différentielles de la chaleur des deux milieux 1 et 2 :

𝑇₁(𝑥) = −𝑄𝑣 2𝜆₁ 𝑥 2 + 𝑎𝑥 + 𝑏 = −46.15𝑥 2 + 202.7𝑥 + 500

𝜑₁(𝑥) = Q𝑣 𝑥 − aλ₁ = 60x − 131.75

𝑇₂(𝑥) = 𝐶₁𝑥 + 𝐶₂ = −241.21𝑥 + 1416.36

𝜑₂ = −C₁λ₂ = 48.24 W ⁄ cm²

C - détermination de T(x=e1) et T(x=e₁+e₂) :

On a e₁= 4 cm T(x=e₁) = -46.15x²+202.7x+500

=-46.15*3² + 202.7*3 +500 = 572˚C

D -détermination de la position et la valeur de T max :

Donc 0 = 2*(-46.15) x + 202.7 92 .3 x = 202.7

Azote 1,73 2,07 2,22 2,36 2,66

TP N°2 Le cycle de Rankine à régénération idéal

- Le rendement thermique du cycle : 46%

Avec ŋ=- ẇ /Qm Ẇ=1283 Qm=2767

- La température d'entrée de l'eau dans la chaudière : 190 C°

Etudier l’effet de la pression du soutirage sur les paramètres suivants :

- le travail net : 1283 KW

- La puissance thermique consommée : 1300 KW

- La fraction de vapeur soutirée de la turbine. 0.77

- La température d'entrée de l'eau dans la chaudière. 190

La Fraction de vapeur soutirée de la turbine du second étage est : 0.1314

Le travail produit dans l’étage à basse pression de la turbine est : 471.7 kw

La chaleur évacuée : (Q dot out) -1.491 kW

La pression max est presque 2512 KPa

Vous aimerez peut-être aussi

On donne : S=1e et P=2(e+1)

=-46.153² + 202.73 +500 = 572˚C