REPUBLIQUE ALGERIENNE DEMOCRATIQUE ET POPULAIRE

UNIVERSITE KASDI MERBAH OUARGLA

Faculté des hydrocarbures et des

énergies renouvelables et des

science de la terre et de l’univers

Mémoire de Projet de Fin d’Étude pour l’Obtention du Diplôme de

MASTER PROFESSIONNEL

Département : Forage et MCP

Domaine : Génie pétrolier
Spécialité : Mécanique des chantiers pétroliers
Présenté par :
KHADEM MILOUD
BOUHAMDA SALAH
KABAR DERBEL

Thème

Résolution du problème de diminution de performance

thermique d'une turbine à gaz M5332B par un cycle
de régénération

Devant le jury composé de :

Mr. NECIB HICHEM M.A.C.C [Link] PRÉSIDENT

Mr. BOUKSSANI MOHAMED M.A.C.C [Link] EXAMINATEUR
Dr. ABIDI-SAAD AISSA M.C.B [Link] ENCADREUR

2016-2017
 Annexes

Annexe 1:

Méthodologie de calcul sous forme de programme (MATLAB):

Ce programme de calcule permet de facilement de faire le dimensionnement de l'échangeur de

chaleur et de déterminer son efficacité.

Programme de KERN

clc; clear all; close all;

%{ ************************ coté tube ************************** %}
bool=false;
us;25:250
L=4.8768;
do=0.0254;
di=0.02118;
p=0.03175;
b=0.97536;
m=0.01;
i=0;
%{ Calcule de Nombre de Tube %}
a=3066042.627./(us.*56.58)
n1=a./(pi.*do.*L);
n=round(n1);
%{ Calcule de la Vitess miassique %}
w1=2*(log(di)/log(exp(1)));
z1=exp(w1);
st=n.*pi.*z1/4
gt=32.067./st;
%{ Calcule de Nombre Reynolds %}
ret=di.*gt*3600./(0.034157.*3.6);
%{ Calcule de cofficient de film interne %}
w2=0.8.*(log(ret)/log(exp(1)));
z2=exp(w2);
jh=0.027.*z2;
w3=(1./3).*(log(0.6973)/log(exp(1)));
z3=exp(w3);
hi=jh.*0.047268.*z3./di;
%{ Calcule de cofficient de film interne corrégie %}
hi0=hi.*di/do;
%{ ************************ coté calander ************************** %}

%{ ****************** Calcule de daimétre de calandre(de) ********** %}

while bool==false
ac=n1.*0.001064;
dc=sqrt(4.*ac./pi);
%{ Calcule de la Vitess miassique (get_) %}
set_=dc.*(p-do).*b./p;
get_=30.21./set_;
%{ Calcule de Nombre Reynolds (rect) %}
w4=2*(log(p)/log(exp(1)));
z4=exp(w4);
de=(4*z4/(pi*do))-do;
rect=de.*get_.*3600./(0.03406.*3.6);
%{ Calcule de cofficient de film externe %}
w5=0.55.*(log(rect)/log(exp(1)));
z5=exp(w5);
jhc=0.36.*z5;
w6=(1./3).*(log(0.7764)/log(exp(1)));
 Annexes

z6=exp(w6);
h0=jhc.*0.040327.*z6./de;
%{ Calcule de cofficient de transfert sale %}
w7=(-1).*(log(1./hi0+0.0001+0.00024+1./h0)/log(exp(1)));
z7=exp(w7);
us1=z7;
y=abs((us-us1)./us)
%{ ************************ coté teste de condition ******************** %}
if y>m
%{ Calcule de Nombre de Tube %}
a=3066042.627./(us.*56.58);
n1=a./(pi.*do.*L);
n=round(n1);
%{ Calcule de la Vitess miassique %}
w1=2.*(log(di)/log(exp(1)));
z1=exp(w1);
st=n1.*pi.*z1./4;
gt=32.067./st;
%{ Calcule de Nombre Reynolds %}
ret=di*gt.*3600./(0.034157.*3.6);
%{ Calcule de cofficient de film interne %}
w2=0.8.*(log(ret)/log(exp(1)));
z2=exp(w2);
jh=0.027.*z2;
w3=(1./3).*(log(0.6973)/log(exp(1)));
z3=exp(z3);
hi=jh.*0.047268.*z3./di;
%{ Calcule de cofficient de film interne corrége %}
hi0=hi.*di./do;

%{ ************************ coté calander ************************ %}

%{ Calcule de daimétre de calander(dc) %}

ac=n1.*0.001064;
dc=sqrt(4.*ac./(pi));
%{ Calcule de la Vitess miassique (get_) %}
set_=dc.*(p-di).*b./p;
get_=30.21./set_;
%{ Calcule de Nombre Reynolds (rect) %}
rect=de.*get_.*3600./(0.03406.*3.6);
%{ Calcule de cofficient de film externe %}
w4=0.55.*(log(rect)/log(exp(1)));
z4=exp(w4);
jhc=0.36.*z4;
w5=(1/3).*(log(0.7764)/log(exp(1)));
z5=exp(w5);
h0=jhc.*0.040327.*z5./de;
bool=false;
else
us=us1;
bool=true;
end
i=i+1;
end
%{ **************** Calcule de perte de charge coté tube ******************
%}
w8=(-0.32).*(log(ret)/log(exp(1)));
z8=exp(w8);
ft=0.028+0.25.*z8
dp=(gt.*3600).^2./(1.271E+15.*0.4455E-03);
dpt=((ft.*2.438./di)+2)*dp*1E-04;
%{ **************** Calcule de perte de charge coté calandre **************
%}
 Annexes

dpc=((get_.*3600).^2).*(3.*dc.*0.18);
dpct=dpc./((1.271e+15).*(de).*(0.00492))*1e-4;
up=h0.*hi0./(h0+hi0);
%{ ********************* coté affichage de resultat ********************** %}
disp('Le cofficient de transfert sal est :');
disp(us);
disp('Le Nombre de Tube est :');
disp(n);
disp('Le diametre de la calandre est :');
disp(dc);
disp('Le Nombre Reynolds coté tube est :');
disp(ret);
disp('Le Nombre Reynolds coté calandre est:');
disp(rect);
disp('La Vitess miassique coté tube :');
disp(gt);
disp('La Vitess miassique coté calandre est:');
disp(get_);
disp('Cofficient de correction est :');
disp(ft);
disp('Perte de charge coté tube est :');
disp(dpt);
disp('Perte de charge coté calandre est :');
disp(dpct);
disp('Cofficient d"echange propre ');
disp(up);
disp('i=');
disp(i);
 Annexes

Annexe 2:

Propriété thermodynamique de l'air (d'après [Link] et [Link] transferts thermiques, Mécanique

des fluides,1989.
 Bibliographie

Bibliographie

[1] ABIDI-SAAD AISSA, KHOUILDI HICHEM, Optimisation des conditions de

marche d'une turbine à gaz type MARCK de l'unité MPP2 de HASSI R'MEL,
mémoire de fin d'étude d’ingénieur d’état, Université Ouargla, juin 2006.
[2] DJAHNINE SAIFEDDINE, MEDGUEDEM SMAIL,L'amélioration de la
rentabilité de la turbine à gaz MS5002C avec un régénérateur de chaleur étude
technico-économique, Université Boumerdess,2014.
[3] DONALD Q. KERN , process heat transfer –Mc Graw Hill Company 1950.
[4] HITASHI, manuel de fonctionnement de la turbine M5332B.
[5] HITASHI, manuel de maintenance de la turbine M5332B.
[6] IFP, manuel d' évaluation économique des projets, édition 2002.
[7] JAN KARSINSKI, Turbomachines, théorie générale, tom (1-1994).
[8] Mm. MEFTAH, cours thermodynamiques 3éme année IMP, FHC ex INH
Boumerdess .
[9] [Link] Abdallah, Amélioration des performances d’une turbine a gaz utilisée
dans l'industrie des hydrocarbures par l'injection de vapeur d'eau a l'amont de la
chambre de combustion, FHC ex INH Université Boumerdess.
[10] M. REGGIO ET [Link], théorie des turbomachines, école polytechnique,
monlréal.
[11] M. SILINI, A. MENIDJEL, Echangeurs de chaleur, formation industrie,2006.
[12] [Link], Raffinage et génie chimique, tome1,2éme édition, Edition TECHNIP
1972.
[13] .[Link], Raffinage et génie chimique, tome2,2éme édition, Edition TECHNIP
1972.
[14] SONATRACH, Documentation interne du complexe Hassi R'Mel.
[15] VAN WYLEN, Thermodynamique appliquée, 2éme édition.
[16] [Link].
[17] [Link].
[18] [Link]
Chapitre I Description générale des turbines

I.1. Introduction
Notre étude a été faite sur des turbines au niveau d'une station de pompage à Hassi R'mel, la turbine
étudiée a un cycle simple, de modèle M5332B et a usage industriel, fournie par HITASHI, elle se
trouve au niveau des modules 2, 3 et 4 pour la compression des gaz moyens pression.
I.2. Historique[2]
Dans l’histoire de la turbine à gaz, on peut distinguer trois périodes:
La première période, celle des précurseurs, est très ancienne puisqu’il est classique de la faire
remonter à Héron d’Alexandrie avec son Éolipile, simple sujet de curiosité ou d’amusement.
Viennent ensuite les premiers dépôts de brevets. Pour les turbomoteurs, en 1791, l’Anglais John
Barber brevète un appareil hybride puisque cette turbine à gaz comportait encore un compresseur.
La deuxième, celle des premières réalisations, commence à la, fin du XIX e siècle et peut être
considérée comme achevée en 1951. Entre 1872 et 1900 environ, les premiers turbomoteurs sont
effectivement construits mais ne peuvent atteindre leur autonomie par suite de l’insuffisance des
rendements de compression et de détente. Par contre, entre 1901 et 1906, les recherches des Français
Armengaud et Le Male aboutissent au premier turbomoteur autonome avec un rendement global à
3%. Entre 1935 et 1945, de nombreuses réalisations apparaissent, notamment dans le domaine
aéronautique où les turbines à gaz bénéficient des actives recherches menées au cours de la dernière
guerre mondiale. Le premier vol d’un avion équipé d’un turboréacteur a lieu en Allemagne, fin août
1939 (moteur HE S 3 monté sur avion Heinkel 178 V1), précédant, en mai 1941, une réalisation
voisine en Grande-Bretagne (moteur e Whittle W 1X monté sur avion Gloster E.28). Enfin, 1951
voit deux premières mondiales avec des turbines à gaz de la firme française Turboméca. Le 18 avril,
c’est l’hélicoptère SO 1120 Arriel 3 qui effectue un premier vol propulsé par un turbomoteur,
l’ARTOUSTE. Le 6 novembre, c’est le premier vol d’un turbo- réacteur à double flux, l’ASPIN,
monté sur le Fouga Gémeaux IV.
La troisième, la période industrielle, commence en 1939. C’est, en effet, au cours des cinquante
dernières années que ces machines se sont développées de façon tout à fait spectaculaire. On peut
citer M. Sedille qui, dès 1948, pressentait avec raison cette évolution : « Il est hors de doute que,
dans les années à venir, un effort considérable permettra de multiplier dans toutes les branches
d’utilisation les installations turbomotrices à gaz ». Actuellement, la turbine à gaz fait partie de notre
environne- ment courant : l’aviation commerciale et militaire utilise quasi exclusivement des

2
Chapitre I Description générale des turbines

machines de ce type pour propulser ses aéronefs. Pour les applications industrielles, la turbine à gaz
est maintenant le concurrent direct des moteurs diesels, et cette évolution est loin d’être terminé.
I.3. Les types des turbines
I.1.3. Turbine à vapeur [18]
Une turbine est un dispositif rotatif destiné à utiliser l'énergie cinétique d'un fluide liquide comme
l'eau ou gazeux (vapeur, air, gaz de combustion), pour faire tourner un arbre supportant les aubes de
la turbine (Fig.I.1).
L'énergie du fluide, caractérisée par sa vitesse et son enthalpie. est partiellement convertie en énergie
mécanique pour entraîner un alternateur, un compresseur, une pompe ou tout autre récepteur
mécanique rotatif. L'ensemble est alors respectivement appelé turbo- alternateur, turbocompresseur,
turbopompe, etc.

Fig. I.1. Turbine à vapeur[10]

I.3.1.1. Principe général de fonctionnement

La turbine à vapeur est un moteur à combustion externe, fonctionnant selon le cycle
thermodynamique dit de Clausius-Rankine. Ce cycle se distingue par le changement d’état affectant
le fluide moteur qui est en général de la vapeur d'eau. Il comprend au moins les étapes suivantes :

• L’eau liquide est mise en pression par une pompe et envoyée vers la Chaudière .

3
Chapitre I Description générale des turbines

• L’eau est chauffée, vaporisée et surchauffée .

• La vapeur est envoyée vers la turbine, où elle se détend en fournissant de l’énergie mécanique .
• La vapeur détendue est condensée au contact d'une source froide sous vide partiel.

Le fluide utilisé est donc le même que celui de la machine à vapeur à pistons, mais la turbine en
constitue une évolution exploitant les principaux avantages des turbomachines à savoir :

Puissance massique et puissance volumique élevées .

Rendement amélioré par la multiplication des étages de détente.

La turbine à vapeur est l’aboutissement d’un type de machines thermiques introduit par les machines
à vapeur à piston. Les contraintes inhérentes à leur conception restreignent généralement leur usage
à l'industrie. Dans ce cas on obtient de l'électricité bon marché car l'énergie thermique n'est pas
« gaspillée » dans un condenseur. On appelle ces turbines, « turbines à contrepression ».

I.3.2. Turbine à gaz [2]

La turbine à gaz est un moteur à combustion interne. De tous les points de vue, elle peut être
considérée comme un système autosuffisant, elle aspire et comprime l'air atmosphérique dans son
compresseur axial, augmente la puissance énergétique de l'air comprimé dans sa chambre de
combustion, en utilisant la chaleur dégagée par la combustion, et convertie cette puissance thermique
en énergie mécanique utile pendant le processus de détente qui a lieu dans la section turbine.
L'énergie mécanique qui en résulte est transmise par l'intermédiaire d'un accouplement à une
machine réceptrice qui peut-être un compresseur, un alternateur

I.3.2.1. Principe de fonctionnement de la turbine à gaz

L'air atmosphérique est aspiré dans le compresseur a flux axial ou il passe à travers de multiples
étages de compression et ensuite se déversent dans les chambres de combustion, ou’ il est chauffé à
la température requise dans la chambre de combustion. Le gaz à haute température résultant est
ensuite détendu à travers la turbine de premier étage (HP) qui entraîne le compresseur à flux axial, et
est détendu ensuite à travers la turbine basse pression (BP), les produits de combustion sont ensuite
évacues dans l'atmosphère.

4
Chapitre I Description générale des turbines

I.4. Description générale de la turbine M5332B [4]

La turbine M5332B est une turbine à gaz à deux arbres, fonctionne suivant le cycle
thermodynamique simple de Brayton, elle se compose essentiellement d’un compresseur axial, de 12
chambres de combustion et de deux roues turbines : une roue haute pression et une autre de
puissance à basse pression. Comme il est illustré dans La Fig.I.2. ces composants sont complétés par
les systèmes d'aspiration, d’échappement, de lancement et de contrôle.
La turbine à gaz possède quatre principaux paliers :
• Le palier n° l et n° 2 supportent le compresseur et le rotor de la turbine à haute pression.
• Le palier n° 3 et n° 4 supportent le rotor de la turbine à basse pression.

ASPIRATION

Fig. I.2. Turbine a gaz M5332B [4]

A. Section de compression
Dans le compresseur, l'air est comprimé par une série d'aubes du rotor qui donnent la force
nécessaire, pour comprimer l'air à chaque étage du compresseur et les aubes du stator guident l'air,
pour le faire pénétrer dans les étages successifs du rotor.
Les aubes du rotor sont insérées dans des rainures et maintenues dans une position axiale par
l'empilage et le bouclage au bout des rainures. Les disques et le demi-arbre sont assemblés pour
maintenir la conicité, ils sont maintenus par des tirants(Fig.I.3).

5
Chapitre I Description générale des turbines

Fig. I.3. Le compresseur axial de la turbine [5]

1
Le rôle essentiel du compresseur axial est :
Assurer l'alimentation des chambres de combustion avec l'air comprimé, pour l'opération de
combustion.

Assurer un débit et une pression aussi élevé pour avoir une grande puissance utile.

Assurer l'air utilisé pour le refroidissement des pièces exposées aux fortes contraintes
thermiques.
Dans la section compresseur la partie stator (corps du compresseur) est composée de quatre éléments
principaux qui sont
A.1. Le corps d’admission
Le Corps coté aspiration du compresseur se trouve à la partie avant. Sa fonction est de diriger l'air de
façon uniforme vers le compresseur et il porte le premier palier du stator(Fig.I.4).

Fig. I.4. Les corps d’admission [2]

6
Chapitre I Description générale des turbines

A.2. Le corps avant du compresseur

Le corps constituant la partie avant du compresseur contient les quatre premiers étages du stator, il
transmet également les charges de structure qui viennent du corps adjacent vers le support avant.
A.3. Le corps arrière du compresseur
Corps constituant la partie arrière du compresseur contient les derniers étages du stator. Les orifices
d'extraction prévus dans ce corps permettent de prélever l'air au niveau du dixième étage du
Compresseur(Fig.I.5). Cet air est employé pour le refroidissement, le contrôle des pulsations au
démarrage et à l'arrêt et assure également des fonctions d’étanchéité.

Fig. I.5. Les ailettes fixes et mobiles du compresseur axial [2]

B. Section combustion
La section combustion de la turbine à gaz M5332B comporte l'enveloppe de combustion qui est
composée de douze corps de combustion à l’extérieure, douze ensembles chapeau et chemise de
combustion, douze ensemble de pièces de transition et douze injecteurs de combustible, deux
bougies d'allumage, deux détecteurs de flamme, douze tubes à foyer et divers garnitures(Fig.I.6).
L'enveloppe de combustion est un élément soudé entourant la partie arrière du corps de refoulement
du compresseur et recevant l'air de refoulement du compresseur à flux axial. Le combustible est
envoyé dans chaque chemise des chambres de combustion par un injecteur de combustible monté
dans le couvercle de cette dernière et pénétrant dans la chemise. Lorsque l'allumage se produit dans

7
Chapitre I Description générale des turbines

les tubes foyer, ils vont allumer le mélange air combustible des autres chambres

Fig I.6. La distribution de l'air dans la chambre de combustion[2]

B.1. Enveloppe de combustion

L'enveloppe de combustion soutient les douze corps de combustion qui referment les douze pièces
de transition, c'est une enceinte soudée qui reçoit l'air de refoulement du compresseur à flux
axial(Fig.I.7). La plaque arrière de l’enveloppe de combustion est boulonnée à la bride verticale du
corps de la turbine, la plaque avant est boulonnée à la bride arrière du corps de refoulement.

Fig. I.7. Enveloppe de combustion[2]

C. Section de la turbine
La section de la turbine se compose de deux étages, haute pression et basse pression :
Etage haute pression : Il est composé d’une directrice, de 36 aubes fixes et de 18 secteurs

8
Chapitre I Description générale des turbines

(2 aubes par secteur), et de la roue turbine haute pression (fig.I.8.), de 80 aubes, qui entraîne le
compresseur axial et les accessoires entraîne par l'arbre de commande, sa vitesse est toujours
constante.

Fig. I.8. La roue turbine haute pression[5]

Etage basse pression : Il est composé :
A - d'une directrice de 36 aubes variables, placé entre l'étage haut pression et basse pression de la
turbine. Ces aubes sont régulées de manière à repartir l'énergie disponible entre les deux turbines.
B - de la roue turbine basse pression (fig.I.9.) qui comporte 80 aubes reliées à leur sommet pour
éviter le fluage, et qui entraîne le compresseur centrifuge (la machine couplée), elle est
mécaniquement indépendante du roue haut pression, ce qui leur permet de tourner à des vitesses
différentes (peut être réglée de 50% à 105% de la vitesse nominale).

Fig. I.9. La roue turbine basse pression[5]

9
Chapitre I Description générale des turbines

D. Section d’échappement
Le gaz d'échappement provenant de la turbine est déchargé dans la bâche d'échappement. Le châssis
d'échappement est refroidi par l'air ambiant qui pénètre à l'intérieure
I.5. Caractéristiques et performances de la turbine M5332B [1]
Une turbine à gaz utilise de l'air atmosphérique, donc, ses performances sont considérablement
influencées par tous les facteurs qui ont un effet sur le débit massique de l’air refoulé au
compresseur, Ces facteurs sont :
- La température ambiante.
- La pression atmosphérique.
A cet égard, les caractéristiques et les performances de la turbine M5332B, qui est l'objet de notre
étude, donnés dans le tableau ci-dessous sont mesurés et calculés dans les conditions atmosphériques
de chantier :
- Température moyen : 40°C, et humidité moyen: 67%
- Elévation de 746m au dessus de la mer, pression atmosphérique de 0.948 kg/cm².
Puissance de la turbine 18500 KW
Puissance requise pour le compresseur centrifuge 14500 KW
Débit d'air 95.4 kg/s
Débit d'échappement 96.2 kg/s
Combustible Gaz naturel
Débit de combustible 1.55 kg/s
Pouvoir calorifique du combustible 10660 Kcal/kg
Compresseur d’air Taux de compression 8.2
Température de 329.4 °C
refoulement
Température de combustion 927°C
La perte de pression dans la chambre de combustion 0.4 kg/ cm²

Température d’échappement 491°C

Tableau I.1. Caractéristique de la turbine à gaz M5332B[1]

10
Chapitre II Étude
tude thermodynamique et amélioration des performances

Partie I : étude thermodynamique

Dans cette partie on suppose que le cycle d'une
d'une turbine à gaz est très souple de telle sorte que ses
performances puissent être améliorés en ajoutant des composantes supplémentaires à un cycle
simple.
II.1. Cycles thermodynamique[2]
thermodynamique
Il existe un grand nombre de cycles thermodynamiques. Nous considérerons ici le cycle Brayton
,qui présente la turbine à gaz.
II.1.1. Cycle simple idéal (ou cycle Brayton)
Le cycle de Brayton pour la turbine à gaz est un cycle à pression constante, c’est -à - dire que la
combustion et l’ échappement se font sous pression constante, pour un cycle simple de turbine à
gaz, la puissance est disponible de façon continue, tandis que dans un moteur à mouvement
alternatif, la puissance n’est pas disponible que pendent la course de détente.
Les différentes étapes du cycle
cle sont indiquées sur la [Link].1,au
[Link] point (1), l’air entre dans le
compresseur, pour sortir à haute pression au point(2), il est ensuite mélangé au combustible dans
la chambre de combustion, les gaz chauds continuent à s’écouler vers le point (3) et dans la
turbine ces gaz sont détendus jusqu'à la pression atmosphérique au point (4).1a turbine produit la
puissance mécanique pourr entrainer le compresseur ou la charge.

Fig. II.1. Cycles d’une turbine à gaz

Dans le cycle visualisé dans la figure ci-dessus,

ci dessus, les 4 points correspondent aux mêmes décrits
auparavant. (1- 2) : compression isentropique (2-3)
(2 : combustion isobare (3-
(3 4) : détente
isentropique (4-1) échappement isochore.

11
Chapitre II Étude thermodynamique et amélioration des performances

D'après le premier principe de la thermodynamique on a :

+ =∆ +∆ +∆ (II -1)
En négligeant les variations d’énergie cinétique et potentielle, et en Considérant que la
transformation dans le compresseur et la turbine est adiabatique => = 0 et ; l'air est assimilé à
un gaz parfait.
Donc: =∆ (II -2)
A. Le compresseur:
On applique l’équation du premier principe entre l’entrée et la sortie en obtient :
→ =ℎ −ℎ
donc ː
= → = − (Kj/Kg) (II -3)
: le travail massique ou bien la quantité d’énergie nécessaire pour comprimer 1kg d’air

Suivant l’isentropique (1-2) en introduisant la relation isentropique =

Le taux de compression : = ; =

La puissance absorbée par la Compression est :

Ẇ → = ṁ" × → = ṁ" ℎ − ℎ (II -4)
B. La chambre de combustion
Dans la chambre de combustion, il n'y a pas de travail , donc la chaleur Q2-3, fournie à la
chambre de combustion de (2) à (3) est exprimée par :
→$ = ℎ$ − ℎ = $ − (II -5)
Ẇ →$ = ṁ% $ −
C. La turbine
Le taux de détente :
& (II -6)
= D'ou : ' = '& et ' = '$
$
En raisonnant de la même manière que pour le compresseur (détente adiabatique, avec un gaz
parfait), nous obtenons :
= = ℎ& − ℎ$ Kj/Kg (II -7)
( $→&

La puissance fournie par l a détente :

Ẇ$→& = ṁ% × = ṁ% ℎ& − ℎ$ (II -8)
$→&

ṁ% = ṁ) + ṁ" ; (ṁ% ) : débit du gaz ; (ṁ) ) : débit du fuel

12
Chapitre II Étude thermodynamique et amélioration des performances

Le travail utile :
Wu = Wt – Wc = *% (T3 – T4) – *" (T2 – T1) (II -9)

La puissance utile :
(II-10)
Ẇu = Ẇt – Ẇc
Le rendement thermique est donné par la relation suivante
+,-./0 34 – 36 7 8 2 – 9
– 7 : –
ηth = = = (II-11)
1 2 1 2 7 8 2 –

Si On considère : C*% = C*" = C* = cst

=
< 9> ?
2 – 9 – –
η th = = 1 − =
= or 9
= 2
(II-12)
2 – < 2> ?
=

Donc le rendement thermique :

η th=1- = 1- (II-13)
@A
B
@

Le rendement augmente avec le taux de compression

II.1.2. Cycle réel d'une turbine à gaz
Les cycles réels diffèrent par l'incorporation des irréversibilités dans les transformations réelles .
La compression et la détente sont pratiquement adiabatiques et irréversibles , ce qui augmente la
puissance du compresseur et réduit celle de la turbine. Il en résulte une augmentation de la
température des sorties du compresseur et de la turbine. Ainsi, la puissance engendrée par la
turbine sera inférieure à celle du cycle idéal, alors que celle requise par le compresseur sera
supérieure.
A. Phase de compression
L’air est comprimé, dans le compresseur, de la pression P jusqu'à la pression de sortie PD .
la compression est accompagné d’une élévation de température de T1 à T D .
1-2': le processus réel est accompagne, des pertes, qui se traduisent par une augmentation
d’entropie
∆ S = S D -S Donc pour la compression de P à PD ,l 'augmentation de la température à travers le
compresseur est plus grande pour le cycle réel, et par conséquent le travail et la puissance
absorbée sont plus importants que pour le cycle idéal.
On définie le travail et la puissance par les expressions suivantes :
= *" T D − (II-14)
1→2′
Ẇ = ṁ" *" T D − (II-15)
′
1→2

13
Chapitre II Étude
tude thermodynamique et amélioration des performances

T D : Température (réelle) à l'état 2'.

2
On appelle le rendement isentropique de compression, le rapport entre
en lee travail changé dans une
transformation isentropique et le travail réel.
IJKLKMN MOPQIJRSMTUP
η isc = <1 (II-16)
IJKLKMN JPPN

B. Phase de combustion
• Bilan de débit massique :
ṁ% = ṁV + ṁ" ; (ṁ% ) : débit du gaz ; (ṁ) ) : débit du fuel
• Bilan énergétique:

Le bilan énergétique s'écrit :

ṁK hD + ṁ) ηcc PCI = ṁ) + ṁK h′3 (II-17)
η cc : rendement de la chambre
mbre de combustion
PCI : Pouvoir calorifique inferieur
erieur du carburant
En développant cette équation et en divisant par le débit d'air, nous aurons :
Où: CpK T D + η]] f PCI = 1 + f C*% T$D (II-18)

f=
ṁ_
; f : rapport de débit massique(fuel/air).
ṁ`

T$D : Température à l'état 3''.

C*% : Chaleur
Chaleur spécifique à pression constante des gaz.
gaz
C*" : Chaleur
Chaleur spécifique à pression constante de l'air.
On définie la perte de pression (∆')
( ) dans la chambre de combustion par l'expression
b
> 2
′
∆' = b

; (∆' c 0,05)

L'élévation de la température due à la combustion est :

C*g b + η
cc f PCI
T$D = (II-19)
1+f C*h

La quantité de chaleur fournie à la chambre de combustion est :

Q b →$b = C*% (T$D −T D ) (II-20)

14
Chapitre II Étude thermodynamique et amélioration des performances

La puissance calorifique due à la combustion :

Qk = mk m Cpm T$D − T D = mk K 1 + f Cpm T$D − T D (II-21)

2′ →3′
C. Phase turbine
A l'entrée de la turbine : P$D = 1 − ∆' × P D
La détente des gaz brulés dans la turbine, de la pression P$D jusqu'à la pression de sortie P&D est
accompagnée d'une diminution de température de T$D à T&D .
Le rendement isentropique de la turbine est :
3 b b ′ > ′
IJKLKMN JPPN 2 9 (II-22)
η isc = = 2 9 = D’où :T&D = T$D − η isd T3 − T4
IJKLKMN MOPQIJRSMTUP 32 9 2> 9

T&D : température réelle à la sortie de la turbine

Dans le cas d'une turbine, le processus de détente est adiabatique, la puissance reçue par
Le fluide W$b > &b est inférieure à la valeur absolue à la puissance W$>& qu'aurait théoriquement
développé une transformation isentropique.
Avec l 'hypothèse des chaleurs spécifiques constantes entre l'entrée et la sortie turbine, on peut
écrire que le travail et la puissance réelle fournie par la détente est :
W$b > &b = hD& – hD$ = Cpm (T&D – T$D )

Wk$b > &b = ṁ` + ṁ_ Cpm T&D – T$D (II-23)

D. Le rendement de l’installation
• Le travail net : est défini comme l'expression suivant:
WQPI = W$b > &b − W > b (II-24)
• puissance nette :
La puissance nette, est définie comme la différence entre la puissance réelle fournie par la
turbine, est celle réellement absorbée par le compresseur.
WQPI = Wk$b > &b − Wk > b (II-25)

• Rendement thermique :
Le rendement thermique est défini comme étant le rapport entre la puissance nette fournie par la
turbine et la quantité de chaleur fournie par la combustion.
3uv4
s(t = (II-25)
w b 2b

15
Chapitre II Étude
tude thermodynamique et amélioration des performances

Partie II : Amélioration
A de performances
Dans cette deuxième partie on propose un procédés afin d’améliorer les performances de la
turbine. Le procédé proposé est récupération des calories d’échappement.
II.2.. Cycle de régénération[1]
régénération
Le cycle de régénération consiste à récupérer la chaleur latente
te des gaz d'échappement et à
l’utiliser pour augmenter la température de l'air fourni
fou i par le compresseur avant qu'il soit envoyé
envoy
dans la chambre de combustion.
combustion

Fig. II.2. Cycle de régénération (turbine à gaz à deux arbres)[1]

arbres

Le transfert de chaleur des gaz d'échappement à l'air refoulé par le compresseur a lieu dans un
échangeur (régénérateur). Dans ce système, la chaleur à fournir
f par le combustible dans la
chambre de combustion
ion afin d'atteindre le point 3,
3 Q2’-3, est réduite à l’aire située sous le tronçon
d'isobare 2'-3 dans le diagramme T-s
T de la figure II.2 ci-dessus,, alors que la chaleur économisée
est représentée par l’aire située
situé sous le tronçon d'isobare 2 - 2', la chaleur fournie par le
combustible dans le cycle simple, Q2-3, étant représentée
ée par l’aire située sous l’isobare 2-3.
Le rendement thermique du système est défini par :
wnet
η th = (II-26)
Qh
Dans
ans le cycle de régénération ([Link].2),
(fig. la chaleur fournie par le combustible est seulement Q2’-3
tandis que dans le cas d'un cycle simple elle est égale à Q2.3. Ainsi, le rendement
ndement du cycle avec
régénération (ηr) sera relié à celui du cycle simple (η
( s) par :
Q2−3
ηr = η s (II-27)
Q2'−3

16
Chapitre II Étude thermodynamique et amélioration des performances

Q2−3
Le terme est toujours supérieur à 1, il représente le gain du cycle de régénération par
Q2'−3
rapport au cycle simple, en termes d'économie d'énergie thermique. Ce gain est d'autant plus
grand que la chaleur Q2’-3 fournie par le combustible est plus petite; ceci correspond à la
différence entre la chaleur Q2-3 du cycle simple et la chaleur échangée à l'intérieur du
régénérateur Q2-2’.
L'avantage de la régénération augmente proportionnellement au potentiel représenté par la
différence de température (T4-T2). Pour cette raison, la régénération n’est pas pratiquée sur les
machines à taux de compression élevé (> 14 à 15), ne disposant pas d’un potentiel suffisant ;
elle est même impossible sur les moteurs à réaction d'avion, où T4 est inférieure à T2.

II.3. Application de la régénération sur la turbine M5332B[1]

La turbine M5332B fournie par HITASHI est prédisposée pour fonctionner en cycle de
régénération ; il suffit simplement d’ajouter le régénérateur, les conduites d’air et de gaz et une
cloison qui servira à orienter l’air vers le régénérateur. La turbine est déjà dotée de brides de
raccordement et présente un potentiel de régénération assez intéressant (T4 – T2 ≈ 160°C)

Fig. II.3. Schéma représenter l'orientation de l'air vers le régénérateur[1]

La figure ci-dessus montre le cloison qui empêche l'air d'entrer directement dans la zone de
combustion et le dirige vers la sortie, qui est raccordée à la tuyauterie d'admission du
régénérateur par une bride. Après avoir traversé le régénérateur, l'air entre à nouveau dans la
chambre de combustion à une température plus élevée par une autre bride qui relie la sortie du
régénérateur à la chambre de combustion.
Le régénérateur est un échangeur à tubes et calandre, les gaz d'échappement y traversent des
faisceaux tubulaires, tandis que d'air circule dans l’espace aménagé entre la calandre et ces
derniers, l’échange thermique entre les deux fluides se fait à contre-courant([Link].4).

17
Chapitre II Étude thermodynamique et amélioration des performances

Fig. II.4. Schéma fonctionnel du régénérateur[1]

II.4. Description des échangeurs de chaleur faisceau et calandre

L'échangeur est constitué par un faisceau de tubes montés sur deux plaques tubulaires et portant
un certain nombre de chicanes. A chaque extrémité, sont fixées les boites de distribution qui
assurent la circulation du fluide à l'intérieur du faisceau en une seule ou plusieurs passes. Le
faisceau est logé dans une calandre, munie de tubulures d'entrée et de sortie pour le second fluide
qui circule à l'extérieur des tubes suivant le chemin imposé par les chicanes([Link].5).
Tous les éléments entrant dans la construction de ces échangeurs ont fait l'objet d'une
normalisation publiée par la T.E.M.A (Standards of Tubular Exchanger Manufacturers
Association) qui spécifie les caractéristiques mécaniques et thermiques correspondant aux
diverses conditions de fonctionnement. Cette normalisation est complétée par d'autres codes tels
que ceux de I’A.P.I. et de l'A S M E [12].

Fig. II.5. Echangeur à tubes et calandre[7]

18
Chapitre II Étude thermodynamique et amélioration des performances
Étude

Les tubes sont fixés à chacune de leurs

leur extrémités dans deux plaques tubulaires. La perforation
des trous dans ces plaques est normalisée; elle s'effectue selon une disposit
sition soit au pas carré,
soit au pas triangulaire comme montre la figure ci dessous.
dessous Le pas triangulaire permet de placer
environ 10% de tubes de plus que le pas carré sur une plaque
plaque tubulaire de diamètre donné, mais
m
en contre partie, la disposition des tubes rend
rend impossible leur nettoyage extérieur par insertion de
grattoirs ou de racloirs à travers le faisceau. Pour ces appareils, il faut avoir recours au nettoyage
chimique et réserver leur emploi pour des fluides propres.

Pas triangulaire Pas carré

Fig. II.6. La disposition des tubes du faisceau [12]

Le faisceau de tubes porte également des chicanes transversales qui ont pour but d'allonger le
chemin du fluide circulant dans la calandre et d'améliorer ainsi
ainsi le transfert à l'extérieur des tubes.
Ces chicanes assurent en outre la rigidité du faisceau et sont solidaires de la plaque tubulaire fixe
au moyen de tirants et entretoises.
entretoises Elles sont généralement constituées par un disque ayant un
diamètre légèrement inférieur à celui
celui de la calandre et comportant un segment libre[12].
li

19
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

III.1. Introduction
Dans ce chapitre, nous allons présenter les équations générales qui régissent les phénomènes de
transfert de chaleur et de pertes de charge dans les échangeurs à tubes et calandre, en vue de les
appliquer dans le calcul de dimensionnement du régénérateur.
III.2. Étude du transfert de chaleur dans les échangeurs à tubes et calandre
III.2.1. Équations fondamentales[12]
Quel que soit le type d'échangeur utilisé, si l'on ne prend en considération que les conditions
d'entrée et de sortie des deux fluides, il est possible d'établir le bilan thermique global de
l'appareil en écrivant que la quantité de chaleur Q perdue par le fluide chaud est égale à celle
prise par le fluide froid, si l'on néglige les pertes thermiques :
Q = m(h2 − h1 ) = M (H 1 − H 2 ) (III-1)
Les lettres majuscules sont réservées pour le fluide chaud, les minuscules pour le fluide froid,
alors que les indices 1 et 2 correspondent respectivement aux conditions d'entrée et de sortie. M
et m désignent les débits massiques des fluides; H et h, les enthalpies des fluides en fonction de
leurs températures T et t.
Par ailleurs, on peut appliquer l'équation de Fourier à l'ensemble de l'appareil
∆Tm
Q = AF = UAF∆Tm (III-2)
∑R
avec,
A : surface totale d'échange offerte par l'appareil;
U : coefficient de transfert global;
∆Tm : différence de température moyenne entre les deux fluides.
F : facteur de configuration qui est fonction de la disposition des fluides l’un par rapport à
l’autre (contre courant, parallèle, etc…) et du nombre de passes ; la valeur de F est lue sur un
abaque avec les températures d’entrée et de sortie des fluides comme paramètres.

III.2.2. Coefficient de transfert global U [12]

III.2.2.1. Expression générale

Pour un élément de tube de longueur dL, le flux de chaleur correspondant dQ traversera cinq
résistances([Link].1) :
R0 = 1/h0 : résistance dans le fluide à l'extérieur du tube, égale à l'inverse du coefficient
du film externe.
Rs0 : résistance due au film d'encrassement déposé sur l'extérieur du tube

20
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

Rt : résistance due à la paroi métallique du tube (cette résistance peut être négligée dans
les calculs)
Rsi : résistance due au film d'encrassement déposé à l'intérieur du tube;
Ri = 1/hi : résistance dans le fluide à l'intérieur du tube, égale à l'inverse du coefficient
du film interne.

Fig. III.1 Résistances thermiques[12]

Comme ces résistances s'expriment en h m2 °C/kcal, il est nécessaire de rapporter tous ces termes
à la même surface. On convient de choisir, en référence, la surface extérieure du tube, ce qui
amène à corriger les résistances intérieures Rsi et Ri en les multipliant par le rapport d0/di des
diamètres extérieur et intérieur du tube. Ces résistances corrigées s'écriront:
Rsi0 = Rsi (d0/di) (III-3)
hi0 = hi (di/d0) (III-4)
Dans ces conditions, la résistance globale au transfert, égale à la somme des quatre résistances
précédentes :

1 1 1
Rs = + R s 0 + R si 0 + = (III-5)
h0 hi 0 U s
Us est le coefficient de transfert sale relatif à l'élément de tube dL. Lorsque cet élément est neuf,
les résistances Rso et Rsi sont nulles.
La détermination des coefficients de transfert locaux Us nécessite l'estimation des coefficients de
film hi et h0 ainsi que des résistances d'encrassement Rsi et Rs0.

21
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

III.2.2.2. Calcul du coefficient du film interne hi

Le calcul du coefficient de film en convection forcée à l'intérieur d'un tube pour les gaz se fait
par l'expression [12]
1
0.14
λ  Cpµ 3  µ 
hi = j h     (III-6)
d i  λ   µ t 
1
 d 3 (III-7)
Dans le cas d’un régime laminaire (Re<2100) : jh = 1.86 Re i 
 L
Tandis que pour un régime turbulent (Re>10000) : jh = 0.027(Re )
0.8
(III-8)
diG
le nombre de Reynolds étant défini ici par: Re = (III-9)
µ
Où G = uρ représente la vitesse massique. µ est la viscosité dynamique et ρ la masse
volumique

III.2.2.3. Calcul du coefficient du film externe h0

Quels que soit la géométrie de l'enceinte autour du tube, tous les auteurs et expérimentateurs sont
d'accord pour présenter l'expression du coefficient de film externe sous une forme analogue à
celle de l'équation (III-6) en régime turbulent. La difficulté réside dans le fait qu'il faut
déterminer un débit massique moyen.
L'écoulement dans la calandre s'effectue selon la disposition des chicanes et l'on peut distinguer
grossièrement deux types de courant :
• un courant principal qui évolue, tantôt perpendiculairement au faisceau entre les chicanes,
tantôt parallèlement au faisceau à l'aplomb des chicanes;
• de nombreux courants secondaires de fuite, dus aux tolérances de construction entre
chicanes et calandre, et entre tubes et chicanes.

Fig. III.2. Ecoulement du fluide à travers la calandre[12]

22
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

Cependant, compte tenu de l'effort de normalisation réalisé dans la construction de ces appareils,
il a été possible d'établir des corrélations convenables permettant d'utiliser une forme d'équation
analogue à l’équation (III-6). Kern définit un seul débit moyen pour l'ensemble de l'appareil, il
utilise à cet effet un diamètre équivalent qu’il définit par les équations (III.10 et 11) selon la
disposition des tubes :
4P 2
Pas carré : De = − d0 (III-10)
πd 0

3.464 P 2 (III-11)
Pas triangulaire : De = − d0
πd 0
Où : P est le pas des tubes, et d0 , le diamètre extérieur des tubes.
Pour des chicanes normalisées dont la hauteur du segment libre représente 25% du diamètre
intérieur de la calandre, Kern calcule la vitesse massique transversale Gct, qu'il utilise en
association avec le diamètre équivalent De dans la formule suivante pour la détermination du
coefficient du film externe, formule valable uniquement en régime turbulent (Re> 2 100) :
1
0.14
λ  Cpµ 3  µ 
h0 = j h     (III-12)
De  λ   µ t 
j h = 0.36(Re )
0.55
Avec : (III-13)
De Gct
Où le nombre de Reynolds modifié, Re =
µ (III-14)
•

Gct = m ṁ étant le débit massique (III-15)

act
Dc
act = ( P − d 0 ) B , aire de passage. (III-16)
P
B , l'espace entre les chicanes et Dc le diamètre intérieur de la calandre.

III.2.2.4. Les températures caloriques des deux fluides

Les équations précédentes permettent le calcul des coefficients de film locaux h0 et hi0 pour un
élément d'échange dL. Il ressort de ces équations que les coefficients de film sont fonction des
propriétés physiques des fluides, donc de leurs températures T et t. Il est donc nécessaire de
définir pour les deux fluides des températures moyennes Tc et tc, appelées généralement
températures caloriques [12], qui servirons de référence pour le calcul des propriétés physiques
qui interviennent dans les expressions des coefficients de film h0 et hi0.

23
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

La différence des températures (Tc-tc), calculée en admettant une variation linéaire de U en

fonction de la température, est toujours sensiblement très voisine de la différence de température
moyenne logarithmique (DTLM), calculée en supposant U constant. Par conséquent, la DTLM
permet le calcul direct du facteur Fc, puis des températures caloriques par les expressions
suivantes :
t c = t1 + Fc (t 2 − t1 )
(III-17)
Tc = T2 + Fc (T1 − T2 ) (III-18)
DTLM − ∆t fr
Où Fc = , avec : ∆tch = T1 − t2 et ∆t fr = T2 − t1 (III-19)
∆t ch − ∆t fr

III.2.2.5. Les résistances d'encrassement

Après un certain temps d'utilisation, il se produit des dépôts à l'intérieur et à l'extérieur des tubes,
dépôts qui introduisent des résistances Rsi et Rs0 au transfert de chaleur, de sorte que le
coefficient de transfert global U diminue. Il est donc indispensable de calculer la surface
d'échange de l'échangeur pour la valeur minimale admissible du coefficient de transfert quand
l'échangeur est sale. On a l'habitude de considérer qu'un appareil tubulaire doit pouvoir
fonctionner pendant une année sans nettoyage. Les valeurs des résistances Rsi et Rs0 des films
après un an de service ont été établies empiriquement, grâce à l'expérience industrielle des
utilisateurs, en fonction de la nature des fluides utilisés. On trouvera dans la T.E.M.A. une liste
de valeurs Rs correspondant aux divers produits [12].
III.3. Etude de la perte de charge dans les échangeurs à tubes et calandre
Les deux fluides qui traversent l'échangeur s'écoulent sous l'effet d'un potentiel mesuré par la
différence de pression ∆P entre l'entrée et la sortie et que l'on appelle la perte de charge. La
valeur de cette perte de charge est fonction de la vitesse de circulation des fluides et de leurs
caractéristiques physiques : masse volumique et viscosité, ainsi que de la géométrie de
l'échangeur.
D'après Kern, la perte de charge totale due à la circulation du fluide à l'intérieur des tubes est:
nc × nt × Gt ²  f t × l 
∆Pt = × + 2 avec : (III-20)
1.271× 10 × d  d i × φt
15

nc : nombre de passes côté calandre
nt : nombre de passes côté tube
Gt : vitesse massique côté tube
d ; densité du fluide côté tube
di : diamètre intérieur des tubes
ft : coefficient de frottement dans le tube

24
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

l : longueur des tubes

0.14
µ
Φt =   ≈ 1 : terme correctif
 µt 
En ce qui concerne le coefficient de frottement, on a :
32
Dans le cas d’un régime laminaire (Re<2100) : ft = (III-21)
Re
f t = 0.00280 + 0.250(Re )
−0.32
Tandis que pour un régime turbulent (Re>2100) : (III-22)

Pour la perte de charge à l'extérieur des tubes, Kern simplifie le problème en ne prenant en
considération que les échangeurs pour lesquels la hauteur du segment libre de la chicane est égale
au quart du diamètre intérieur de la calandre et emploie le diamètre équivalent De pour établir la
valeur du nombre de Reynolds :
De × Gct
Re ct = (III-22)
µ'
La perte de charge est calculée par l'équation de Fanning modifiée :
nc × f ct × Gct ² × ( N c + 1) × Dc
∆Pc = (III-23)
1.271× 1015 × d ' × De × φc
Le coefficient de frottement fct est donné en fonction de Rect et la distance parcourue par le
fluide est égale à (Nc + 1) Dc .
Avec :
Gct : la vitesse massique côté calandre
Nc : nombre de chicanes
Dc : diamètre intérieur de la calandre
d' : densité du fluide côté calandre

25
Chapitre III Théorie sur le calcul des échangeurs de chaleurs

Le calcul de dimensionnement de l’échangeur se fait selon la méthode de Kern qui se résume

dans l’organigramme suivant :

Fig III.3. Organigramme de la méthode de Kern[1]

26
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Dans ce chapitre, nous allons en premier lieu déterminer toutes les conditions thermodynamiques
aux entrées et sorties du régénérateur (l’échangeur), ensuite nous procéderons au
dimensionnement du régénérateur en appliquant la méthode de Kern.
IV.1. Calcul du rendement thermique actuel (cycle simple)
Le rendement thermique est donné par la relation suivante :
wnet
η th =
Qh

w& net = 18500 kW

T3
Q& h = Q& 2−3 = m& 3 C p (T3 − T2 )
T2

le débit de mélange fuel/air (ṁ3) qui en résulte se dirige ver la chambre de combustion égale la
somme de 95% du débit d'air aspiré à l’entrée du compresseur et le débit fuel gaz injecté :
m& 3 = 92.183 kg/s , T = 927°C, T = 329.4°C
3 2

Cp (gaz de cc, T3 = 927°C) = 0.2943 kcal/kg = 1.2302 kj/kg

Cp (air, T2=329.4°C) = 0.2523 kcal/kg.°C = 1.0546 kj/kg
T3 C p (T3 ) + C p (T2 ) 1.2302 + 1.0546
Cp = = = 1.1424 kj / kg
T2 2 2
Q& 2 −3 = 92.183 × 1.1424 × (927 − 329.4) = 62933.172 kW

18500
η th = = 29.39%
62933.172

IV.2. Les paramètres thermodynamiques de la turbine M5332B

Selon la fiche technique de

la turbine M5332B, mentionné
dans le tableau I.1
T1 = 40°C, P1 = 0.948 kg/cm²
m& 1 = 95.4 kg/s Fig. V.1 Schéma de la turbine à gaz à
deux arbres
T2 = 329.4°C, τ = 8.2
T3 = 927°C
T4 = 491°C, P4 = 0.948 kg/cm²
m& 4 = 96.2kg/s

m& comb = 93.2 kg/min

Fig IV.1. Les paramètres thermodynamiques da la TAG[1]

27
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Détermination de la température de l’air à la sortie du régénérateur T2'

Sur les turbines M5332B, le rendement thermique (ηth), qui se situe autour de 29% pour le cycle
simple, pourrait atteindre facilement de 32% à 34% dans le cas des cycles avec régénération, il
peut même aller jusqu’ à 36% d’après le constructeur. Cette dernière valeur fixerait les limites
de faisabilité de la régénération sur ce type de turbine. Autrement dit ; au-delà de cette limite, les
dimensions de l’échangeur seraient telles qu’on ne pourrait pas l’adapter à l’installation.
Adoptons donc la valeur ηth = 36% pour le rendement thermique de l’installation avec
régénération. La connaissance de cette valeur du rendement va nous permettre de déterminer la
température de l’air à la sortie du régénérateur. On a en effet,
wnet
ηth =
Qh
Expression dans laquelle :
w& net = 18500kw (voir tableau I.1)
T3
Q& h = Q& 2'−3 = m& 3 C p ' (T3 − T2 ' )
T2

Avec :
m& 3 = 92.183kg/s et T3 = 927°C d’après les données du problème.
La chaleur spécifique moyenne se calcule à partir des chaleurs spécifiques de l’air à T2 et des gaz
de combustion à T3 ; ces dernières étant obtenues par le logiciel "DESIGN II FOR
WINDOWS" [12] on a donc :
Cp (T2' ) ≈ Cp (T2) = Cp (air, T2=329.4°C) = 0.2523 kcal/kg°C =1.0547kj/kg
Cp (gaz de combustion, T3 = 927°C) = 0.2943 kcal/kg = 1.2302kj/kg
T3 C p (329.4°C ) + C p (927°C ) 1.0547 + 1.2302
Cp = = = 1.1424kj / kg
T2' 2 2
L’expression du rendement devient :
w& net
η th = T3
& 3 ×Cp
m × (T3 − T2' )
T2'

D’où l’on tire :

w& net
T2' = T3 − T3
= 439.022°C
η th × m
& 3 ×Cp
T '
2

T2' = 440°C

28
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Détermination de la température des gaz à la sortie du régénérateur (T4' )

Pour déterminer cette température, il suffit à présent de considérer le bilan thermique du
régénérateur supposé bien isolé de l’extérieur au moyen d’un calorifugeage qui empêche toute
déperdition de chaleur. Le bilan s’écrit dans ces conditions :

m& 4 C p
T4'

T4
(T '
&)
4 − T4 = m2 C p
T2'

T2
(T
2
'
− T2 )
et la température T4’ s’obtient alors comme suit
T2'
m& 2 C p
T4' = T4 −
T2
T4'
(T
2
'
− T2 ) (V-1)
m& 4 C p
T4

Avec
Cp (air, T2=329.4°C) = 0.2523 Kcal/kg.°C
Cp (air, T2'=440°C) = 0.2575 Kcal/kg.°C
On aura :
T2' 0.2523 + 0.2575
Cp = = 0.2549 kcal / kg °C
T2 2
Pour le calcul de la chaleur spécifique moyenne des gaz, procédons de la manière suivante :
T4'
posons dans une première étape C p = Cp (gaz d’échappement, T4) = 0.2709kcal/kg
T4

T4'
Pour cette valeur,
leur, l’équation (V-1)
(V donne T4' = 392.96°C ; ensuite on recalcule le C p avec
T4

cette valeur de T4’ et on poursuit ainsi jusqu’à ce que le processus converge. On aboutit
finalement à la valeur :
T4' = 391.95°C
Toutes
outes les conditions thermodynamiques des fluides aux entrées et aux sorties de notre
régénérateur (échangeur) sont connues, elles sont portées sur la représentation schématique
suivante :

Fig IV.2. Les paramétres

paramétres thermodynamiques d'un échangeur
29
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

IV.3. Dimensionement de l’échangeur

IV.3.1. Données de procedé requises pour le dimensionnement de l’échangeur
Les données suivantes doivent être déterminé ou calculé pour nous permettre de dimensionner
l’échangeur.
• Les températures d’entrée et sortie des deux fluides (T2,T2’,T4,T4’)
• Débits côté calandre et côté tubes
• Pertes de charge calculées pour chaque circuit
• Facteur d’encrassementde chaque fluide
• Températures et pressions de calcul
• Nature de métal des tubes
• Nature des fluides ou composition (présence d’hydrogène, par exemple)
• La masse volumique des fluides aux conditions d’entrées et de sorties

IV.3.2. Disposition des fluides côté tube ou côté calandre

Il n’existe pas de règle très précises pour ce choix qui résulte souvent d’un compromis entre les
advantages et les inconvinients de placer un fluide déterminé dans des conditions de
fonctionnement définies, côté tube ou côté calandre.
- Côté tube :
• L’eau (principalement l’eau de mer)
• La vapeur d’eau
• Les gaz chargés en hydrogène
• Les fluides hauts pression
• Les fluides les plus chauds
• Les produits corrosifs, peu visqueux, encrassant.
- Côté calandre :
• Les fluides basse pression
• Les liquides très visqueux
• Les produits peu encrassant et non polymérisant

IV.3.3 Calcul de la différence de température moyenne logarithmique entre les deux

fluides
Dans le cas d’un échange de chaleur en contre courant, la différence de température moyenne
logarithmique s’exprime par

30
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

∆t1 − ∆t 2
∆Tm = DTLM =
 ∆t 
Ln 1 
 ∆t 2 
où, ∆t1 = T4 − T2' et ∆t 2 = T4' − T2
On a par conséquent :

∆Tm =
(491 − 440 ) − (391.95 − 329.4) = 56.58°C
 491 − 440 
Ln 
 391.95 − 329.4 
Soit ;
∆Tm = 56.58°C
IV.3.4. Détermination des températures caloriques
Nous rappelons ici que les températures caloriques seront utilisées dans la procédure de calcul du
régénérateur comme température de référence pour l’évaluation des propriétés physiques des
fluides.
A - Température calorique côté calandre (T2c)
Cette température se calcule à l’aide des équations suivantes:
T2c = T2 + Fc (T2' – T2)
DTLM − ∆T fr
où Fc =
∆Tch − ∆T fr

avec : ∆T fr = T4' − T2 et ∆Tch = T4 − T2'

On a pour notre cas :

∆T fr = 391.95 − 329.4 = 62.55°C ⇒ ∆T fr = 62.55°C

et
∆Tch = 491 − 440 = 51°C ⇒ ∆Tch = 51°C
ce qui donne
56.58 − 62.55
Fc = = 0.5169
51 − 62.55
et par suite :
T2c = 329.4+ 0.5169 (440 – 329.4) = 386.57 °C
Soit :

T2c = 386.57 °C

31
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

B – Température calorique côté tube (T4c)

De la même manière, on peut écrire du côté tube à l’aide de l’équation suivante :
T4c = T4' + Fc (T4 – T4') avec toujours Fc = 0.5169
On obtient :
T4c = 391.95 + 0.5169 (491 – 391.95) = 443.15 °C
Soit :
T4c = 443.15 °C

IV.3.5. Propriétés physiques des deux fluides évaluées aux températures caloriques
D'après le logiciel "DESIGN II FOR WINDOWS", on trouve en introduisant les valeurs des
températures caloriques ci-dessus et la composition de l’air et des gaz d’échappement, les pro
priétés physiques contenues dans le tableau ci-dessous :

Côté froid (calandre) Côté chaud (faisceaux tubulaire)

-air- -gaz d'échappement-
'
Cp = 0.2545 kcal/kg Cp= 0.2683 kcal/kg
ρ' = 4.0169 m3/kg ρ = 0.445467 m3/kg
µ' = 0.03406 cP µ = 0.034157 cP
λ' = 0.040327 kcal/mh°C λ = 0.047268 kcal/mh°C

Tableau IV.1 : Propriétés physiques des fluides de l’échangeur

IV.3.6. Calcul de l'échangeur

Les calculs de l'échangeur seront conduits selon la méthode de Kern exposée dans l’ouvrage de
Pierre Wuithier intitulé "Raffinage et génie chimique". La procédure se résume dans
l’organigramme présenter dans le chapitre III.
Dans ce qui suit, nous allons présenter uniquement les détails des calculs qui concernent la
boucle I et qui correspondent à la solution retenue pour la boucle II, solution qui, dans les
conditions thermodynamiques imposées, permet de minimiser la perte de charge dans le
régénérateur afin de ne pas perturber la puissance envoyée au compresseur centrifuge.

IV.3.6.1. Les caractéristiques de l'échangeur à prévoir

Type d'échangeur : à tubes et calandre.
Type de circulation : contre courant.
Longueur des tubes : l = 16 ft = 4.8768m (valeur maximale de la longueur, pour une surface
d’échange donnée, qui maximise la section de passage du gaz).

32
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Tubes: type -14 BWG, diameters d0 = 1' =0.0254 m , di = 0.02118m.

Disposition de l'emplacement des tubes : pas carré, p =1.25 in= 0.03175 m (pour
maximiser la surface de passage de l'air côté calandre, et par conséquence minimiser les
pertes de charge correspondantes Nature des tubes : acier faiblement allié (allié : à cause de
haut niveau de température, et faiblement allié pour que le prix ne soit pas trop élevé).
Nature de la calandre : acier faiblement allié 15Mo3, et d'une plaque tubulaire en acier allié.
Espacement des chicanes B = 0.97536m, (cette valeur est suffisamment importante pour
contribuer à minimiser la perte de charge côté calandre).
On installe trois régénérateurs identiques en parallèle à deux passe côté tube (nt=1) et une
seule passe également côté calandre (nc=1), à cause des quantités considérables d'air refoulé
(96.2kg/s). Ainsi donc, le flux de chaleur échangé entre les fluides et les débits de ces
derniers seront partagés entre les trois régénérateurs et on aura pour chacun d’entre eux :
- Le flux de chaleur par régénérateur :
1
Q& = m& 2 C p
3
T2'

T2
(T
2
'
)
− T2 = 3 066 042.627 kcal / h

Soit,
Q& = 3066042.627kcal / h
- Le débit d'air dans chaque régénérateur est :
m& 2 90.63
& '2 =
m = = 30.21kg / s
3 3
Soit,
& '2 = 30.21kg / s
m
- Le débit de gaz d'échappement dans chaque régénérateur est :
m4 96.2
& '4 =
m = = 32.067kg / s
3 3
& '4 = 32.067kg / s
Soit, m

IV.3.6.2. Estimation de l'intervalle des valeurs du coefficient de transfert sale (Us)

On trouve dans la référence [12] les coefficients de transfert sale de plusieurs pairs de fluides, et
le cas le plus proche De notre cas, correspond à la vapeur côté tube et à l’air côté calandre ; la
référence sus-mentionnée donne pour ce cas : Us dans l’intervalle [ 25.250]kcal / hm ²°C .
Dans ce qui suit, nous allons effectuer une 1ére itération avec Us = 250 kcal/hm²°C et une
seconde itération avec Us = 25 kcal/hm²°C. Ces deux itérations vont nous permettre de voir si le
choix de l’intervalle ci-dessus est judicieux ou non.

33
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

1ère itération avec Us = 25 kcal/hm²°C

A. Calcul du nombre de tubes
Calculons d’abord la surface d'échange :
Q 3066042.627
A= = = 2167.6m 2
U s × ∆Tm 25 × 56.58
Soit,
A = 2167 .6m 2
La surface d'échange correspondant à un seul tube sera :
s = π × d 0 × l = π × 0.0254 × 4.8768 = 0.3892m 2
Soit,
s = 0.3892 m 2
Le nombre de tubes avec une seule passe côté calandre (nc = 1) est :
A 1083.79
Nt = = = 5570
s × nc 0.3892 × 1
Soit,
N t = 5570tubes
B. calcul du coefficient du film interne corrigé (hi0)
La section de passage des gaz d'échappement avec une seule passe côté tube (nt = 1) est donnée
par :
Nt d i2 3712 0.021182
st = ×π × = ×π × = 1.9627m 2
nt 4 1 4
Soit,
st = 1.9627 m 2
A présent, calculons la vitesse massique (Gt)
m& 4' 32.067
Gt = = = 16.3403kg / m² s
st 1.9627
Soit,
Gt = 16.3403kg / m² s
Ensuite, calculons le nombre de Reynolds (Ret)
d i × Gt 0.02118 × 16.3403 × 3600
Re t = = = 1.0132 × 10 4
µ 0.034157 × 3.6
Soit,
Re t = 1.0132 × 10 4

34
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Avec Ret > 104, l'écoulement des gaz côté tube est turbulent et on a :
jh = 0.027 × (Re t ) 0.8 = 0.027 × (1.0132 × 10 4 ) 0.8 = 43.2444
Le coefficient du film interne (hi) se calcule à présent à l’aide de la relation :

λ 1
µ 0.14
hi = jh × × (Pr) 3 × ( )
di µt
où,
µ 0.14
( ) ≈1
µt
le nombre de Prandtl :
Cp × µ 0.2683 × 0.034157 × 3.6
Pr = = = 0.6973
λ 0.047268
On a par calcul :
1
0.047268
hi = 43.2444 × × (0.6973) 3 = 85.5811kcal / hm² s
0.02118
Soit,
hi = 85.5811kcal / hm²°C
Finalement, le coefficient du film interne corrigé (hi0) est :
di 0.02118
hi 0 = hi × = 85.5811× = 71.3625kcal / hm²°C
d0 0.0254
Soit,
hi0 = 71.3625kcal/hm²°C

C. Calcul du coefficient du film externe (h0)

La section de passage de l'air dans côté calandre (sct) est donnée par :

× ( p − d 0 )× B
Dc
sct =
p
avec : B = 0.97536m, p = 0.03175 m, d0 = 0.0254 m
Pour calculer cette section, il faut d’abord déterminer le diamètre de la calandre (Dc) :
4 × Ac
Dc =
π
ou : Ac = N t × f a et f a = 1.65d 02
Les calculs pour sct se font donc comme suit

35
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

f a = 1.65 × (0.0254 )² = 0.00106m²

Ac = 5570 × 0.00106 = 5.9265m²
4 × 5.9265
Dc = = 2.7470m
π
2.7470
sct = × (0.03175 − 0.0254) × 0.97536 = 0.5359m²
0.03175
Calculons ensuite la vitesse massique côté calandre (Gct) :
m& 2' 30.21
Gct = = = 56.3769kg / m² s
sct 0.5359
Soit,
Gct = 56.3769kg / m² s
Le nombre de Reynolds (Rect) est donné par:
De × Gct
Re ct =
µ'
Pour un pas carré , le diamètre équivalent se calcule comme suit :
4 × p² 4 × 0.03175²
De = − d0 = − 0.0254 = 0.02515m
π × d0 π × 0.0254
Ceci donne un nombre de Reynolds :
0.02515× 56.3769 × 3600
Re ct = = 4.1599 ×10 4
0.03406 × 3.6
Cette valeur correspond à un régime d'écoulement turbulent, on a donc d’après l’équation
suivante:
jh = 0.36 × (Re ct ) 0.55 = 0.36 × (4.1599 ×10 4 ) 0.55 = 124.9671
Soit,
jh = 124.9671
Le coefficient du film interne (h0) sera donc :

λ' µ ' 0.14

1
h0 = jh × × (Pr ) × ( ' ) ' 3

De µt

µ ' 0.14
Avec ( ) ≈1
µt'
Le nombre de Prandtl :
C p' × µ ' 0.2545 × 0.03406 × 3.6
Pr ' = = = 0.7764
λ '
0.040327
On a donc :

36
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

1
0.040327
h0 = 124.9671× × (0.7764) 3 = 184.3027 kcal / hm²°C
0.02515
Soit,
h0 = 184.3027kca/hm²°C

D. Calcul du coefficient de transfert sale (Us)

1 1 1
= + Rs 0 + Rsi 0 +
U s h0 hi 0
Avec, h0 = 184.3027kca/hm²°C, hi0 = 71.3625kcal/hm²°C
L'air est un fluide propre, la résistance d'encrassement côté air (Rs0) est très faible, Rs0 = 10-4
hm²°C/kcal [12].
Les gaz d'échappement ont une composition qui très proche de celle de l'air, leur résistance
d'encrassement est aussi très faible mais elle est légèrement supérieure à celle de l'air, à cause des
produits de combustion du gaz naturel qu'il contient, Rsi = 2 ×10-4 hm²°C/kcal [12].
La résistance d'encrassement corrigée (Rsi0) se calcule comme suit :
d0 0.0254
Rsi 0 = Rsi × = 0.0002 × = 2.3985 × 10 − 4 hm²°C / kcal
di 0.02118
Soit,
Rsi 0 = 2..3985 × 10 −4 hm²°C / kcal

1 1 1
= + 0.0001 + 0.00023985 + = 0.0198hm²°C / kcal
U s 184.3027 71.3625
soit,
1
Us = = 50.5591kcal / hm²°C
0.0198
Ce résultat indique que le coefficient global Us doit être en réalité supérieur à la valeur de 25
kcal/hm²s et c’est pour savoir s’il doit être supérieur à 50 kcal/hm²°C que nous effectuons une
seconde itération.

2eme itération Us = 250 kcal/hm²°C

A. Calcul du nombre de tube

Calculons d’abord la surface d'échange :
Q 3066042.627
A= = = 216.75m2
U s × ∆Tm 250 × 56.58
Soit,
A = 216.75m 2

37
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Le nombre de tubes avec une seule passe côté calandre (nc = 1) avec s = 0.3892, d'après la
première itération est :
A 216.75
Nt = = = 557
s × nc 0.3892 ×1
Soit,
N t = 557tubes
B. Calcul du coefficient du film interne corrigé (hi0)
La section de passage des gaz d'échappement avec une seule passe côté tube (nt = 1) est donnée
par :
Nt d 2 557 0.021182
st = ×π × i = ×π × = 0.1962m 2
nt 4 1 4
Soit,
st = 0.1962 m 2
A présent, calculons la vitesse massique (Gt) :
m& 4' 32.067
Gt = = = 163.4035kg / m² s
st 0.1962
Soit,
Gt = 163.4035kg / m² s
Ensuite, calculons le nombre de Reynolds (Ret) :
d i × Gt 0.02118 × 163.4035 × 3600
Re t = = = 1.01×10 5
µ 0.034157 × 3.6
Soit,
Re t = 1.01 × 10 5
Avec Ret > 104 , l'écoulement des gaz côté tube est turbulent et on a :
jh = 0.027 × (Re t ) 0.8 = 0.027 × (1.01× 10 5 ) 0.8 = 272.8537
Le coefficient du film interne (hi) se calcule à présent à l’aide de la relation :

λ 1
µ 0.14
hi = jh × × (Pr) × (
3
)
di µt
Où,
µ 0.14
( ) ≈1 et Pr = 0.6973 , d'après la première itération
µt
On a par calcul :

38
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

1
0.047268
hi = 272.8537 × × (0.6973) 3 = 539.98kcal / hm²°C
0.02118
Soit,
hi = 539.98kcal / hm²°C
Finalement, le coefficient du film interne corrigé (hi0) est obtenu en faisant appel à l’équation
(IV.4) :
di 0.02118
hi 0 = hi × = 373.04 × = 450.2668kcal / hm²°C
d0 0.0254
Soit,
hi0 = 450.2668 kcal/hm²°C

C. Calcul du coefficient du film externe (h0)

La section de passage de l'air dans côté calandre (sct) est donnée par :

× ( p − d 0 )× B
Dc
sct =
p
avec : B = 0.97536 m, p = 0.03175 m, d0 = 0.0254 m
Pour calculer cette section, il faut d’abord déterminer le diamètre de la calandre (Dc) :

4 × Ac
Dc =
π
ou : Ac = N t × f a et f a = 0.001064m² d'après la première itération
Les calculs pour sct se font donc comme suit
Ac = 557 × 0.001064 = 0.5927m²
4 × 0.5.927
Dc = = 0.8687m
π
0.8687
sct = × (0.03175 − 0.0254) × 0.97536 = 0.1695m²
0.03175
Calculons ensuite la vitesse massique côté calandre (Gct) :
m& 2' 30.21
Gct = = = 178.2796kg / m² s
sct 0.1695
Soit,
Gct = 187.2796kg / m² s
Le nombre de Reynolds (Rect) est donné par:
De × Gct
Re ct =
µ'

39
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Pour un pas triangulaire, et De = 0.02515m d'après la première itération, Le nombre de Reynolds

se calcule comme suit :
0.02515×187.2796 × 3600
Re ct = = 1.31×105
0.03406 × 3.6
Cette valeur correspond à un régime d'écoulement turbulent, on a donc d’après l’équation
(IV.13):
jh = 0.36 × (Re ct ) 0.55 = 0.36 × (1.31 × 10 5 ) 0.55 = 235.3942
Soit,
jh = 235.3942
Le coefficient du film interne (h0) sera donc :

λ' 1
µ ' 0.14
h0 = jh × × (Pr ' ) 3 × ( )
De µ t'

µ ' 0.14
avec : ( ) ≈1 et Pr ' = 0.7764 , d'après la première itération
µt'
On a donc :
1
0.040327
h0 = 235.3942 × × (0.7764) 3 = 347.1617 kcal / hm²°C
0.02515
Soit,
h0 =347.1617kca/hm²°C

D. Calcul du coefficient de transfert sale (Us)

1 1 1
= + Rs 0 + Rsi 0 +
U s h0 hi 0
Avec, h0 = 347.1617kca/hm²°C, hi0 =450.2668 kcal/hm²°C,
Rsi = 2.3985 ×10-4 hm²°C/kcal, Rs0 = 10-4 hm²°C/kcal.
1 1 1
= + 0.0001 + 0.00023985 + = 0.0054hm²°C / kcal
U s 347.1617 450.2668
Soit,
1
Us = = 183.7760kcal / hm²°C
0.0054
Ce résultat signifie que le coefficient global Us doit être en réalité inférieur à la valeur de 250
kcal/hm²s, Ce qui signifie que: 25 kcal/hm²°C < Us < 250kcal/hm²°C

40
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

IV.3.7. Estimation de la valeur exacte du coefficient de transfert sale (Us) et caractéristiques

de l’échangeur optimisé
On estime laa valeur exacte du coefficient de transfert sale (Us) et les caractéristiques de
l’échangeur optimisé par itération comme indiqué dans l’organigramme ci-dessus,
ci ce calcul
itératif a été exécuté à l’aide d’une petite procédure en langage MATLAB..
Les résultats obtenus sont :
Coefficient de transfert sale : Us =135.37 kcal/hm²°C
Nombre de tubes : Nt = 1031
Diamètre de la calandre : Dc =1.1821
= m
Nombre de Reynolds côté tube : Ret =5.474 × 10 4
Nombre de Reynolds côté calandre : Rect = 9.666 × 10 4
Vitesse massique côté tube : Gt = 88.2791Kg/m²s
Vitesse massique côté calandre : Gct = 131Kg/m²s
Coefficient de transfert propre : UP = 140.9 kcal/hm²°C
La surface d’échange : A = 401.40 m2

Fig IV.3. Caractéristiques de l’échangeur optimisé

41
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

IV.3.8. Calcul de la perte de charge dans le régénérateur

A. Côté tube (∆Pt)
nc × nt × Gt ²  f ×l 
∆Pt = × t + 2
1.271× 10 × d  d i × φt
15

0.445467
Où : nc = 1, nt = 1, Gt = 88.2791Kg/m²s, d= 3
kg / m 3 ,
10
0.14
µ 
di = 0.02118 m, Φt =   ≈ 1, l = 4.8768 m.
 µt 
ft = 0.00280 + 0.250 Ret-0.32
ft = 0.00280 + 0.250 (5.474 × 10 4 )-0.32 = 0.0356
1×1× (88.279 × 3600 )²  0.0356 × 4.8768 
∆Pt = −3
× + 2 = 1.088 × 10 − 4 kg / cm²
1.271× 10 × 0.445467 × 10
15
 0 . 02118 × 1 
Soit, ∆Pt = 1.088 10-4kg/cm²
B. Côté calandre (∆Pc)
nc × f ct × Gct ² × ( N c + 1) × Dc
∆Pc =
1.271× 1015 × d ' × De × φc
ρ ' 4.0169
Où : nc = 1, Gct = 131 kg/m²s, d 10 3 = 10 3 ,
'=

0.14
 µ' 
De = 0.02515 m, Φc = =  ' 
 ≈ 1, Dc =1.1821 m,
 µt 
Or Rect = 9.666 × 10 4 , d'après le graphe de l’annexe 1 fct = 0.18
l 4.8768
= = 5 ⇒ l = 3B ⇒ N c = 4 (On a 4 chicanes)
B 0.97536
1× 0.18 × (131× 3600)² × (4 + 1) ×1.1821
∆Pc = −3
= 9.0348 ×10 −5 kg / cm²
1.271×10 × 4.0169 ×10 × 0.02515×1
15

Soit,
∆Pc = 9.0348 10-5 kg/cm²
Les pertes de charges calculées sont inférieures a celle des pertes de charges admissibles
(0,5 kg/cm²)[12].

Le régénérateur (échangeur) dimensionné a les caractéristiques suivantes :

42
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Fig. IV.4 Schéma simplifié du régénérateur

Type d’échangeur à tubes et calandre

Type de circulation Contre-courant
Nombre de tubes 1031
Longueur des tubes 4.8768 m
Diamètre intérieur des tubes 0.02118 m
Diamètre extérieur des tubes 0.0254 m
Pas entre les tubes (carré) 0.03175 m
Nature des tubes Acier faiblement allié
Nombre de passes côté tube Une seule
Diamètre de la calandre 1.1821 m
Nature de la calandre Acier faiblement allié
Nombre de passes côté calandre Une seule
Espace entre les chicanes 0.97536 m
Nombre de chicanes 4 chicanes

Tableau IV.2. Caractéristiques du régénérateur

IV.4. Calcul des performances de la turbine après régénération

Notons d’abord pour les calculs qui vont suivre que la puissance nette de l’installation reste
quasiment la même après régénération, étant donné que les pertes de charge dues à la circulation
des fluides dans le régénérateur ne représentent qu’environ 0.2% de cette puissance nette, ce qui
n’affecte pas la valeur du rendement de manière significative. On prendra donc la valeur de

43
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

18 500 kW pour la puissance nette dans le calcul du rendement de l’installation aussi bien dans le
cas simple que dans celui avec régénération. Sur cette base, nous allons donc évaluer le
rendement de l’installation actuelle, puis celui de l’installation avec régénération en vue de les
comparer et d’en tirer des conclusions quant au gain qu’on pourrait réaliser grâce à la
régénération.
IV.4.1. Calcul du rendement thermique de la turbine avec régénération
T3
Q& 2 '−3 = m& 3 C p (T3 − T2 ' )
T2'

m& 3 = 92.183 kg/s , T = 927°C, T ’ = 440°C

3 2
Cp (gaz de cc, T3 = 927°C) = 0.2943 kcal/kg = 1.2302 kj/kg
Cp (air, T2'=440°C) = 0.2575 kcal/kg°C = 1.0764 kj/kg
T3 C p (T3 ) + C p (T2 ' ) 1.2302 + 1.0764
Cp = = = 1.1533 kj / kg
T2' 2 2
Q& 2 ' −3 = 92.183 × 1.1533 × (927 − 440) = 51774.114 kW

wnet 18500
η th = = = 35.73 ⇒ η th = 35.73 %
Qh 51774.114
Soit,
ηth = 35.73 %

IV.4.2. Evaluation de l’économie d’énergie réalisée

( )
'
T2
Q&2−2' = m
& 2 Cp T2' −T2
T2

Cp
T2'
=
C p (T2 ) + C p T2' ( )
T2 2
Cp (air, T2=329.4°C) = 0.2523 Kcal/kg.°C = 1.0546kj/kg
Cp (air, T2'=440°C) = 0.2575 Kcal/kg.°C = 1.0764kj/kg
T2' 1.0546 + 1.0764
Cp = = 1.0656kcal / kg °C
T2 2
Q& 2−2' = 90.63 ×1.0656 × (440 − 329.4) = 10680.299kW
Soit,
Q2-2’ = 10 680.299 kW

Le rendement de la chambre de combustion est :

T3
m& 3 C p (T3 − T2 ) 92.183 × 1.1424 × (927 − 329.4) 62933.172
η cc = = = = 0.90944
T2

d&c × PCI 1.553 × 10660 × 4.18 69199.8164

44
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

Soit,
η cc = 90.94%
Le débit de gaz naturel économisé sera :
Q& 2−2 ' 10680.299
d&ce = = = 0.26356 kg / s
PCI ×η cc 10660 × 4.18 × 0.9094
Donc on pourrait économiser quotidiennement 22771.584 kg de GN, soit 8311628.16 kg par
an. Au niveau du complexe de Hassi R’Mel, il y a toujours trois turbines de ce type en service,
l’économie annuelle du complexe serait donc d’environ 25 millions de kg qui correspond à 30
millions de Nm3 de GN si toutes ces turbines étaient dotées de régénérateurs.
IV.5. Etude économique d'un échangeur de chaleur
Après que nous avons identifié la quantité de fuel gaz récupérée qui est due a la quantité de
chaleur des gaz d’échappement récupérée. Il faut calculer le coût de cette nouvelle installation et
la durée de vie de fonctionnement de cette installation (pour récupérer ce coût).
IV.5.1. Calcul de quantité de chaleur récupérée en BTU
On à :
1MBTU = 10 BTU ; 1 = 1.055 Kj
On a Q =117.85 (Kj/Kg)
.
Q = = 111.71 BTU = 111.71 × 10−6
.
IV.5.2. La valeur de la quantité de récupérée en dollar ( $)
1 = 10.80 $ [18]
Q =111.71 × 10 × 10.80 = 1206 × 10 $
IV.5.3. Calcul de débit de fuel récupéré
d# = 0.26356 (Kj/s)
$# % .' (
Débit volumique de fuel : q = = = 0.3295 )( /s
& .

q = 0.3295×(3600×24×355.4) = 10117811.52 )( /An

Avec un débit volumique de 0.3295 )( /s on a récupéré une quantité de chaleur de
111.71 × 10 MTBU
IV.5.4. Le prix de fuel récupéré dans une année
0.3295)( /s 1206× 10 $
10117811.52 )( /An x

45
Chapitre IV Calcul des rendements de la turbine et Dimensionnement du régénérateur

× ' × *+
X= = 37032 $
.(',
IV.5.5. Le coût de l’échangeur de chaleur
D’après nos recherches dans les sites internet spécialisés dans la fabrication mécanique on a
trouvé une approximation de coût de l’usinage d’un mètre au carré de l’acier allié de l'ordre 160
$[16].
On a utilisé le logiciel Solid-Works pour construire(déssiner) notre échangeur de chaleur avec les
dimensions qu' on a calculés. Ce logiciel peut calculer la surface, la longueur et l’épaisseur de
chaque élément
Surface des tubes :401.38 )'
Surfaces des chicanes :1.4 )'
Surface de la calendre : 115.48 )'
Surface de porte tube : 0.57 )'
La surface totale usinée est : 518.83 )'
Le coût de l’échangeur de chaleur : 518.83 × 160 = 83012.8 $
Si on divise le coût de l’échangeur sur le gain obtenu dans le fuel récupéré on trouve la duré de
fonctionnement (D. ) de l’installation pour récupéré son argent.
IV.5.6. Coefficient d'amortissement
/û0 1 2 3456 ( '.
D. = = = 2.24
2 6475 8 ( ('

D. = 2 ans et 85 jour

46
 Chapitre I
Description générale des
turbine
 Conclusion générale

Conclusion générale
Cette étude nous a permet d'enrichir nos connaissances théorique et pratique dans le domaine
d'engineering et de design des équipements stratégiques des unités de traitement de gaz, ainsi que
sur le principe de fonctionnement des équipements énergique tels que les turbines et les
échangeurs de chaleur .

Vu l'importance de ces équipements durant l'opération de traitement de gaz, la SONATRACH

prévoit un grand investissement pour assurer leur continuité de marche afin de produire un
production stable et durable, dans cet objectif que cette étude a été entamée afin d'aboutir à des
résultats fiables et réalisables.

Les résultats trouvés dans notre étude montrent que :

- Le rendement thermique actuel des turbines à gaz opérants en cycle simple est de l'ordre
de 29.39 % avec une consommation importante en fuel gaz.
- Le rendement thermique de la turbine avec le cycle régénéré est de l'ordre de 36%.
- Le gain de rendement de 6 points à pour équivalence en gain de consommation en fuel
gaz de 30 million Nm3 par année.
- On a besoin de 2 ans et 85 jours pour récupérer le prix d'un échangeur.

A la lumière des résultats trouvés, on peut suggérer la recommandation suivante:

De prévoir l'emplacement d'une batterie d'échangeurs composée de 3 unités d'échange avec

les caractéristiques suivantes :

a. Longueur des tubes 4.8768 m

b. Diamètre extérieur des tubes de 0.0254 m.

c. Disposition en pas carré.

d. La surface d'échange de chaque unité est de 401.4 m2.

e. Nombre de tubes dans chaque unité est de 1031 tubes .

f. Nuance des tubes est acier faiblement allié.

47
 Chapitre I
Description générale de la
turbine M5332B
 Introduction générale

Introduction générale
Les besoins en puissance des installations industrielles pour l’entraînement des différents
équipements, s’imposent sur tous les plans et surtout dans le domaine pétrolier.
Les machines d’entraînements doivent remplir certaines obligations sur tout lorsqu’il s’agit des
équipements puissants, tels qu’on rencontre aux stations de compression et de réfrigération, et
dans les centrales thermiques et les stations nucléaires, ….
Pour des raisons économiques et de puissance nécessaires, seules les turbines ont été utilisées
dans les grandes installations. En les comparant aux autres moteurs, la turbine à gaz présente une
double particularité :
 La turbine à gaz est une machine à écoulement continu, donc ne comporte pas de variation
périodiques de l’état de fluide, dans toute section donnée.
 La turbine à gaz (au sens large du terme), est le siège de l’ensemble de transformation
constituant le cycle thermodynamique décrit par le fluide, donc, la turbine à gaz est assimilée aux
moteurs à combustion interne mais les transformations s’effectuent dans des enceintes séparées et
juxtaposées dans l’espace.
A cause des exigences de croissance de la demande en énergie et son prix de revient élevé , et à
des problèmes qui touchent à l’environnement et qui sont exacerbés par les taux de gaz de
combustion dégagés dans l’atmosphère. Nous intéressons dans le cadre de ce projet à l’étude de
la faisabilité d’une régénération sur les turbines biarbres de Hassi R’mel afin :
 D’augmenter le rendement thermique pour satisfaire l'exigence de la demande en énergie.
 De réduire la consommation en gaz combustible sur les turbines de gaz utilisées pour
comprimer les gaz de vente dans les stations de compressions du gaz naturel au niveau du
complexe gazier de Hassi R’Mel
 De baisser les émissions gazeuses chaudes vers l’atmosphère pour protéger l'environnement.
Pour réaliser ces objectifs nous commencerons dans le premier chapitre par une description
détaillée de la turbine à gaz étudiée, tout en montrant les différents organes constituants cette
turbine. Dans un second chapitre, une théorie des turbines à gaz et des échangeurs de chaleur a
été présentée. Dans le troisième chapitre on a donné une théorie sur le calcul des échangeurs de
chaleur. Dans le dernier chapitre nous avons présenté le détail de calcul des rendements de la
turbine étudiée et le dimensionnement du régénérateur. Puis on a terminé par une conclusion.

1
 La nomenclature

La nomenclature

Symboles Unités Significations

T K Température ambiante.
T K Température de l'air au refoulement du compresseur axial.
T K Température de l'air à la sortie du régénérateur côté calandre.
T K Température de gaz brûlé à la sortie de la chambre de combustion
T K Température de gaz brûlé à l'échappement.
T K Température de gaz brûlé à la sortie du régénérateur côté tube.
P Kg/cm2 Pression ambiante.
P Kg/cm2 Pression de l'air au refoulement du compresseur axial.
P Kg/cm2 Pression de gaz brûlé à la sortie de la chambre de combustion.
P Kg/cm2 Pression de gaz brûlé à l'échappement.
∆P Kg/cm2 Perte de charge dans la chambre de combustion .
∆P Kg/cm2 Perte de charge côté tube.
∆P Kg/cm2 Perte de charge côté calandre.
τ / Taux de compression dans le compresseur axial.
η / Rendement thermique de la turbine.
Q KW Chaleur absorbée par le système.
Ẇ KW Puissance nette fournie par le système.
d Kg/s Débit de la combustible.
PCI Kcal/kg Pouvoir calorifique inférieur du combustible.
ṁ Kg/s Débit d'air à l'aspiration du compresseur axial.
ṁ Kg/s Débit d'air de combustion.
ṁ Kg/s Débit de combustible.
ṁ Kg/s Débit de gaz dans la chambre de combustion .
ṁ Kg/s Débit des gaz d'échappement.

̇ Kg/s Débit d'air de combustion dans chaque régénérateur.

̇ Kg/s Débit des gaz d'échappement dans chaque régénérateur.

C J/kg.k Chaleur spécifique côté tube.

C J/kg.k Chaleur spécifique côté calandre.

S m2 Surface de passage de l'air (côté calandre) dans le régénérateur.

 La nomenclature

D m Diamètre de la calandre.

G Kg/m2.s Vitesse massique côté calandre.

Re / Nombre de Reynolds côté calandre.

Pr / Nombre de Prandtl côté calandre.

D m Diamètre équivalent.

U Kcal/hm2.°C Coefficient de transfert global sale.

R hm2.°C/Kcal Résistance d'encrassement à l'extérieur des tubes.

R hm2.°C/Kcal Résistance d'encrassement à l'intérieur des tubes.

R hm2.°C/Kcal Résistance d'encrassement corrigé à l'intérieur des tubes.

Φ / Facteur correctif côté tube.

f / Coefficient de friction côté tube.
Φ / Facteur correctif côté calandre.
f / Coefficient de friction côté calandre.
N / Nombre de chicane dans la calandre.
HP / Turbine haut pression.
GG / Générateur de gaz.
BP / Turbine basse pression.
BTU / Brittish Thermal Unit (unité anglo-saxonne d'énergie)
 En premier, nous remercions le Tout Puissant ALLAH,
notre créateur qui nous a donné la force d’accomplir ce
travail.
Nous tenons à remercier vivement, le chargé du suivi de
ce travail, Monsieur Dr ABIDI-SAAD AISSA , pour ses
encouragements, ses conseils précieux et sa disponibilité.
Nous remercions vivement le jury qui à accepter de
juger ce modeste travail, M. BOUKSSANI MOHAMED et
M. NECIB HICHEM.
Nous présentons nos chaleureux remerciements aux
enseignants du département Forage et MCP : BOUKSSANI
MOHAMED, CHEBOUT MOHAMED SEDIK pour leurs aides et
orientations durant notre formation en particulier.
Sans oublier nos amis BENHAMIDA MOHAMED AMINE
et BEKKOUCHE MOHAMED EL AMINE pour leurs efforts
pour nous aider.
Et à tous ceux qui ont participé de prés ou de loin à la
réalisation de ce travail.
 Résumé

‫ﻣﻠﺧص‬

‫ﺗزاﯾد اﺳﺗﮭﻼك اﻟطﺎﻗﺔ ﻓﻲ ﻣﺟﺎل اﻵﻻت اﻟﺗورﺑﯾﻧﯾﺔ ﯾﺣﺗﺎج إﻟﻰ اﺳﺗﺛﻣﺎرات ﻛﺑﯾرة وھﺎﻣﺔ وذﻟك ﻣن ﺧﻼل اﻟﺑﺣث ﻋن طرق‬
. ‫ﻓﻌﺎﻟﺔ وﻣﻣﻛﻧﺔ ﻟﺗﺣﺳﯾن ﻣردود ھذه اﻷﺧﯾرة‬

‫وﻓﻲ ھذا اﻟﻌﻣل اﻗﺗرﺣﻧﺎ ﺗﻐﯾﯾر دورة اﻟﺗورﺑﯾﻧﺎت اﻟﻐﺎزﯾﺔ ﻣن دورة ﺑﺳﯾطﺔ إﻟﻰ دورة ﺗﺟدﯾد ﻋن طرﯾق إﺿﺎﻓﺔ ﻣﺣول‬
‫ ﻣﻊ‬،‫ واﻟﺣد ﻣن اﻹﻧﺑﻌﺎﺛﺎت اﻟﻐﺎزﯾﺔ اﻟﺳﺎﺧﻧﺔ ﻓﻲ اﻟﺟو واﻟﺗﻘﻠﯾل ﻣن اﺳﺗﮭﻼك اﻟوﻗود‬، ‫ﺣراري ﻟرﻓﻊ ﻣﺳﺗوى أداء اﻟﺗﯾرﺑﯾﻧﺎت‬
. ‫ﻣﻌرﻓﺔ وﺗﺟدﯾد أﺑﻌﺎد ھذا اﻟﻣﺣول واﻟطرﯾﻘﺔ اﻟﻣﺳﺗﻌﻣﻠﺔ ﻓﻲ ذﻟك‬

Résumé

La demande en énergie dans le domaine de turbo machines exige investissement considérable et

important, qui peut être réduit par la recherche d'une source de chaleur disponible dans le procès
utilisé.

Dans ce contexte notre travail consiste à modifier le cycle simple d'une turbine à gaz par un cycle
régénéré dont le but d'augmenter son rendement, de réduire les émissions gazeuses chaudes dans
l'atmosphère et d'optimiser la consommation en combustible.

Abstract

The demand in energy in the domain of turbo machinery requires a considerable and important
investment the can be reduced by the research of a source of available heat in the used process.

In this context our work consists in modifying the simple cycle of a gas turbine by a regenerate cycle
of which the object to increase his efficiency, to reduce the hot gas radiant in the atmosphere and to
optimize the consumption in fuel.
 Sommaire

Sommaire

Introduction générale .................................................................................................................. 1

Chapitre I : Description générale des turbines

I.1. Introduction................................................................................................................................... 2

I.1. Historique ...................................................................................................................................... 2

[Link] types des turbines .................................................................................................................... 3

I.3.1. Turbine à vapeur.................................................................................................................... 3

I.3.1.1. Principe général de fonctionnement .................................................................................... 3

I.3.2. Turbine à gaz ....................................................................................................................... 4

I.3.2.1. Principe de fonctionnement de la turbineà gaz ................................................................... 4

I.4. Description générale de turbine M5332B .................................................................................... 5

I.5. Caractéristiques et performances de la turbine M5332B ........................................................... 10

Chapitre II : Étude thermodynamique et amélioration des performances

II.1. Cycles thermodynamique ........................................................................................................... 11

II.1.1. Cycle simple idéal (ou cycle de Brayton) ........................................................................... 11

II.1.2. Cycle réel d'une turbine à gaz ............................................................................................. 13

II.2. Cycle de régénération ................................................................................................................ 16

II.3. Application de la régénération sur la turbine M5332B ........................................................... 17

II.4. Description des échangeurs de chaleur faisceau et calandre..................................................... 18

Chapitre III : Théorie sur le calcul des échangeurs de chaleurs

III.1. Introduction .............................................................................................................................. 20

III.2. Étude du transfert de chaleur dans les échangeurs à tubes et calandre ................................... 20

ІІI.2.1. Équations fondamentales .................................................................................................. 20

III.2.2. Coefficient de transfert global U ....................................................................................... 20

III.2.2.1. Expression générale ..................................................................................................... 20

III.2.2.2. Calcul du coefficient du film interne hi ...................................................................... 22

 Sommaire

ІІI.2.2.3. Calcul du coefficient du film externe h0 ..................................................................... 22

III.2.2.4. Les températures caloriques des deux fluides ........................................................... 23

III.2.2.5. Les résistances d'encrassement ................................................................................. 24

III.3. Etude de la perte de charge dans les échangeurs à tubes et calandre ..................................... 24

Chapitre IV : Calcul des rendements de la turbine et dimensionnement du

régénérateur

IV.1. Calcul du rendement thermique actuel (cycle simple) ............................................................. 27

IV.2. Les paramètres thermodynamiques de la turbineM5332B ...................................................... 27

IV.3. Dimensionement de l’échangeur .............................................................................................. 30

IV.3.1. Données de procedé requises pour le dimensionnement de l’échangeur ......................... 30

IV.3.2. Disposition des fluides côté tube ou côté calandre ........................................................... 30

IV.3.3. Calcul de la différence de température moyenne logarithmique entre les deux fluides .... 30

IV.3.4. Détermination des températures caloriques ..................................................................... 31

IV.3.5. Propriétés physiques des deux fluides évaluées aux températures caloriques .................. 32

IV.3.6. Calcul de l'échangeur ........................................................................................................ 32

IV.3.6.1. Les caractéristiques de l'échangeur à prévoir ................................................................. 32

IV.3.6.2. Estimation de l'intervalle des valeurs du coefficient de transfert sale (Us) ................... 33

IV.3.7. Estimation de la valeur exacte du coefficient de transfert sale (Us) et caractéristiques de

l’échangeur optimisé ..................................................................................................................... 41

IV.3.8. Calcul de la perte de charge dans le régénérateur ............................................................. 42

IV.4. Calcul des performances de la turbine après régénération ...................................................... 43

IV.4 .1 Calcul du rendement thermique de la turbine avec régénération ...................................... 44

IV.4.2. Evaluation de l’économie d’énergie réalisée ..................................................................... 44

IV.5. Etude économique d'nu échangeur de chaleur ...................................................................... 45

IV.5.1. Calcul de quantité de chaleur récupérée en BTU ............................................................ 45

IV.5.2. La valeur de la quantité de récupérée en dollar ( $) .......................................................... 45

IV.5.3. . Calcul débit de fuel récupéré ........................................................................................... 45

 Sommaire

IV.5.4. Le prix de fuel récupéré dans une année ........................................................................... 45

IV.5.5. Le cout de l’échangeur de chaleur ..................................................................................... 46

IV.5.6. Coefficient d amortissement.............................................................................................. 46

Conclusion générale …………………………………………………………...…………………47

Bibliographie
Annexes