Rev Gin Therm (1997) 36, 289-301

@ Elsevier, Paris

Identification du coefficient de transfert

thermique sur la surface externe d’un khangeur
de chaleur ;3Ltubes et ailettes planes
en fonction de I’humiditk de I’air,
en I’absence de condensation
Riad Benelmir*, Mohamed Khalfi, Michel Feidt
LEMTA, CNRS URA 875, Universite Henri-Poincart? (Nancy I), 2, av de la For.3 de Haye,
545 16 Vandozuvre-/es-Nancy cedex, France

(Recu le 16 septembre 1996 ; accept& le 4 fevrier 1997)

Abridged English version at the end of the text

Summary - Identification of the heat transfer coefficient over the external area of a finned tubes heat exchanger
with respect to the moisture content of the air without condensation. The following study, which is rather oriented
towards experimentation, shows the influence of the humidity content of air on heat transfer. This first article concerns
heat transfer between the external fluid (moist air) and the internal fluid (water containing glycol, whose thermal behavior
inside circular tubes is well-known) in a heat exchanger of the same type as those used in automotive air conditioning
(horizontal copper tubes and plane aluminium fins), in the absence of condensation. The most difficult part of this
experimental work is the measurement and control of the air humidity, since one has to make sure that the measurement
incertainties are not significant compared to the precision of the calculation of the heat transfer coefficient. The conclusion
is that, for this type of exchanger, the heat transfer coefficient decreases with air humidity in the absence of condensation
(dry wall). Some correlations have been developed with respect to the relative air humidity. An analog experimental
investigation, but this time carried out in the presence of condensation (partially or completely wetted wall), is about to
be completed; the obtained results will be communicated later on.

R&sum6 - La pre’sente Ltude, qui est plutBt tourne’e vers l’exp&imentation, met en avidence /‘influence de I’humidite de
/‘air sur le transfert thermique. Ce premier article concerne le transfert de chaleur entre le Puide externe (de /‘air humide)
et le Puide interne (de I’eau glycole’e dont le comportement thermique d Yinte’rieur d’un tube lisse est bien connu) duns un
dchangeur de chaleur du me^me type que ceux utilisk en climatisation automobile (tubes horizontaux en cuivre et ailettes
planes en aluminium), en /‘absence de condensation. La difficultk de ce travail expe’rimental rkide duns la mesure et le
contr6le de I’humiditk, avec le souci d’e’viter que /‘incertitude de la mesure ne devienne importante devant la prtkision
du calcul du coefficient de transfert thermique. En conclusion, il ressort que - du moins pour ce type d’ichangeur - le
coefficient de transfert thermique de’croit avec I’humidite’ en /‘absence de condensation (paroi s&he). Des corre’lations ont
ttL itablies en fonction de l’humidite’ relative de /‘air. La me^me investigation exp&imentale, mais menie cette fois-ci en
prkence de condensation (paroi partiellement ou compktement mouillc’e), est en tours de re’alisation ; /es r&w/tats en
seront communique’s ulte’rieurement.

Nomenclature j facteur de Colburn

lf longueur equivalente d’ailette . m
coefficient numerique M masse molaire . [Link]-’
Fl ezt surface exterieure totale. . . .. m2 n nombre de mesures
Atub surface exterieure des tubes nus . . m2 nombre d’ailettes
b coefficient numerique Nf
NT nombre de rangees de tubes
CP chaleur massique .. .. . . . [Link]-l.K-’
m Nu nombre de Nusselt
Dh diametre hydraulique . . .
P puissance thermique . . . W
De diametre externe des tubes . . . .
h coefficient d’echange convectif. . W.mp2 .K’ Pr nombre de Prandtl
P tot pression totale. . Pa
Re nombre de Reynolds
* Correspondance et tires a part R cd resistance thermique de conduction K.W-’

289
 R Benelmir, M Khalfi, M Feidt

Rth rgsistance thermique K.W--’ paramktres gbom&riques sur le coefficient de

St nombre de Stanton transfert. Parmi les parametres Btudiks on retrouve
T temperature....................... K en premier lieu l’espacement entre les ailettes. Les
t temp&ature....................... “C ou K mkthodes de calcul sont g&&ralement baskes sur la
u conductance globale surfacique W.mV2.KP1 mkthode de Wilson (1915) ou une de ses variantes
V 7n(LI vitesse maximale de l’&oulement de (Shah, 1990 ; Khartabil et Christensen, 1992).
I’air............................... m.s -*
x’t pas des ailettes m Avant de prksenter les correlations reportkes
Y titre molaire dans la littkrature, on dbfinit deux rapports gko-
mktriques souvent rencontrks qui sont le rapport r
Symboles grecs de la surface externe totale de transfert thermique,
AT &art de temperature K Aezt, sur la surface totale des tubes en absence
ATl, &art logarithmique moyen de d’ailettes, A tube, et le rapport s du pas d’ailettes,
temp&ature....................... Xf, sur l’kpaisseur d’ailette, 6.f :
conductivit6 thermique W.m-‘.K-’
efficacit6 de la surface ailettee 7’ = A,Lt/At,,~i, ; s = X,f /6j (1)
Bpaisseur ailettes. m
efficacite d’ailette Les figures 1 et 2 montrent la variation de r et du
viscosit6 dynamique . kg.m-‘.s-’ diamhtre hydraulique Dh en fonction du pas d’ai-
masse volumique kg.m-” lette (l’kpaisseur d’ailette Sf, est con&ante) pour
taux d’humiditt? notre 6changeur expkrimental, d&rit ult&ieure-
humidit relative merit.
Indices
a air set
ah air humide
0 contact 14
cd conduction
c c6t6 externe r 12
en encrassement
.f ailette
c6t.B interne
PI pas longitudinal
1: vapeur d’eau 6
IL’ eau 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

Fig 1. T en fonction de S.

1I INTRODUCTION
Fig 1. r with

5,6E-3 ~j
respect to s.

5,OE-3 j ~
L’utilisation des Bchangeurs de chaleur compacts
4.5E-3
est t&s Gpandue dans les industries de climatisa-
tion et de conditionnement de l’air. La complexit E 4,OE-3
de l’koulement de l’air sur la surface externe de 3$x-3
ces derniers exclut, pour le moment, le d&eloppe- ii
ment d’une expression analytique du coefficient de 3,OE-3

transfert de chaleur. 11 est done nkessaire de faire 2,5E-3

appel & l’exp&imentation dans le but de construire 2,OE-3
des modkles empiriques utilisables. 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22
Le present travail concerne l’identification du
S
coefficient de transfert de chaleur d’une batterie
froide (Bchangeur compact destink au refroidisse- Fig 2. Dh en fonction de s.
ment de l’air), en l’absence de condensation, et met Fig 2. Dh, with respect to s.
en kvidence l’influence de l’humiditb de l’air sur les
performances de l’kchangeur.
Turaga, Lin et Fazio (1966) ont Btudik l’influence
du pas d’ailettes sur le coefficient de transfert. 11s
trouvent que le facteur de Colburn dkcroit avec
le pas d’ailettes et proposent des corrblations en
2I BIBLIOGRAPHIE presence et en absence de condensation :
j (saris cond) = 0, 05r-“.24 (Reoh)-“” (2)
Durant ces vingt dernikres annees, plusieurs Etudes avec : 300 < Reu,l < 1500
ont BtB me&es sur ce sujet, mais la majorit& d’entre
elles se limite B montrer l’influence de certains 20 < r < 50

290
 Identification du coefficient de transfert thermique

j (avec cond) = 0,04 r-o,23 s1,15 (Re~h)-“~~ (3) si Repl est important, le coefficient de transfer-t
croit avec le pas d’ailette
avec : 300 < ReDh < 1500
- pour Xf/De > 0,33, le pas d’ailettes n’a aucune
20 < r < 50 influence sur le coefficient de transfer-t, quelle que
10 < s < 22 soit la valeur du nombre de Reynolds.
Rich (1973) a Btudie l’influence de l’espacement Idem, Jacobi et Goldschmidt (1990) comparent
des ailettes sur le coefficient de transfer% a partir les coefficients de transfert dune paroi s&he
des resultats experimentaux sur huit Bchangeurs a et dune paroi humide en donnant des relations
quatre rangees de tubes. 11conclut que le coefficient correspondant a chacun des cas :
de transfer-t est independant du pas d’ailettes et j (paroi humide) = 0,104 Regk423 (8)
correle le facteur de Colburn j au nombre de
Reynolds, qui depend du pas longitudinal des tubes j (paroi s&he) = 0,145 Re,k4X4 (9)
(espacement des tubes dans le sens de l’ecoulement
de l’air) : avec :
400 < ReDh < 4000
j = 0,195 (R~J$“~~ (4)
avec : Kayansayan (1994) montre, en etudiant dix
1 < v,,, < 9 m.s-’ configurations differentes d’echangeurs a quatre
rangees de tubes, que le coefficient de transfert
0 < Nf < 21 croit avec le pas d’ailettes. I1 donne une relation
qui tient compte du diametre hydraulique et du
MC Quiston (1978) montre que le coefficient de paramhtre T :
transfer-t croit avec le pas d’ailette. Le nombre de
Reynolds depend du diametre exterieur des tubes :
j = 0,0014 + 0,2618 (ReD,)-0’4 F”.15 (5) avec
avec :
500 < Rem, < 30000
100 < ReDe < 4000
11,2 < T < 23,5
El Mahdy et Briggs (1979) avancent que le
coefficient de transfer% decroit avec le pas d’ailette La figure 3 presente la correlation de Kayan-
et donnent une relation ou le nombre de Reynolds sayan pour differentes valeurs de r.
depend du diametre hydraulique : Les relations cites precedemment, faisant inter-
venir un nombre de Reynolds dependant du
j = 0,ll RekL3g (6) diametre hydraulique, sont representees sur la
figure 4. L’observation de ces correlations montre
Canton et Fornasieri (1983) affirment qu’il bien une certaine disparite et justifie que l’identi-
n’existe pas de parametre geometrique influant fication du coefficient de transfert sur la surface
plus particulierement sur le coefficient de transfer%, externe de ce type d’bchangeur reste un sujet d’ac-
a l’exception du nombre de rangees de tubes. 11s tualite. Certains auteurs ont essay6 de montrer
presentent une relation qui prend en consideration l’influence dune surface &he et humide sur le
le nombre de Reynolds (dependant du pas longitudi- coefficient de transfer-t, mais apparemment aucune
nal des tubes) et le nombre de rangees de tubes qui etude n’a et6 me&e sur l’influence de l’humidite
se revele avoir une legere influence sur le coefficient relative sur le coefficient de transfer-t. Ceci sera le
de transfert : theme du present travail.
j = 0,202 Repe’37 NFal (7)

Chen et Ren (1988) ont estime que l’influence

du pas d’ailettes sur le coefficient de transfer-t
n’est pas aussi importante, et que les methodes 3I MtTHODE D’IDENTIFICATION
utilisees dans les travaux precedents, generalement
derivees de la methode de Wilson (1915), ne sont
pas assez precises et appropriees pour mettre en Notre procedure d’identification du coefficient
evidence cet effet. Pour differents pas d’ailettes, ils d’echange est basee sur la methode de Khartabil
utiliserent une technique bake sur la sublimation et Christensen (1992) qui, comme la major-it6 des
du naphtalene qui consiste a estimer le coefficient autres methodes, part de l’equation g&&ale expri-
de transfer% thermique a partir du coefficient de mant la resistance thermique globale de l’echan-
transfert de masse, en utilisant l’analogie entre ce geur :
dernier et le transfer% de chaleur ; ils conclurent
que : Rth = Rth,z + Rth,ct + Rth,cd + Rth,en + Rth .e 01)
- pour Xf/De < 0,33 : soit :
si Repl est faible, l’influence du pas d’ailettes est
Rth = &h,z + CR + (rl ;A)e (12)
insignifiante ;

291
 R Benelmir, M Khalfi, M Feidt

0,009
i
0,008

0,007

0,006

0,005
400 600 800 1000 1200 1400 1600 1800 2000
Re (Dh)
Fig 3. Corre’lation de Kayansayan.
Fig 3. Correlation of Kayansayan.

i AY====al I ,
' 0,008

0,006

0,004

400 600 800 1000 1200 1400 1600 1800 2000

Re(Dh)
Fig 4. j en foncrion de &oh.
Fig 4. j with respect to ReDh.

avec : exterieur des tubes, selon :

CR = &[Link] + &h,cd + [Link] (13) N’ILD~ = b Re& Prli3 (14)

La resistance de contact ailette-tube dans un A partir de la definition du nombre de Nusselt, on

Bchangeur compact a ailettes serties a Bte estimee deduit une expression pour le coefficient d’echange
par L’HGte et Bardon (1989) a 4,6.10e4 K.m’.W-l. convectif :
La resistance totale de conduction est Cvaluee a h = bRe$, PI-‘/~X
e
partir de la configuration de l’echangeur et permet De
d’estimer cette valeur a 4.10~~ K.W-l. La resistance
d’encrassement est toujours difficile a evaluer car En remplacant h dans l’equation (12) par l’ex-
elle depend des fluides utilises et de l’etat de la pression ci-dessus, on obtient :
surface d’echange. Nous supposons que l’air est
1
propre (sans particules) et que l’eau glycolee (fluide Rth = &h,i + CR + (16)
interne a l’echangeur) est non encrassante. qbRe& PI-~/~%
>e
Le transfer-t thermique convectif du cot6 interne
(ecoulement de l’eau dans des tubes droits et lisses) En posant
est tout a fait maitrise. Comme la paroi interne
des tubes est depourvue d’ailettes, l’efficacite de la Rth,z + CR = Ro (17)
surface interne v%est unitaire.
Pour le transfert thermique tote externe (air),
(18)
le nombre de Reynolds dependra du diametre et C = (rl b)e (19)

292
 Identification du coefficient de transfert thermique

l’equation (16) prend la forme suivante : L’equation (25) peut Btre reformulee sous la
for-me suivante :
f&h = n Rth ln ReDe n ln ReDe
C Re;= Wo + RQ
c
1
Re& Wo 1 c 1 (Re&Wo)*
C’est une equation a trois inconnues : C (ou l/C), -- - Ro = 0 (28)
n hReDe c n h ReDe
& et a. La determination du coefficient b est lice a
celle de C. L’inconnue Ro est plut8t une contrainte, c Re& WO c Re& WO
1 1
comme on le verra plus tard.
Nous utiliserons la methode des moindres carres, La resolution du systeme d’equations (261, (27)
qui consiste a minimiser, par rapport aux inconnues et (28) est conduite par iteration et est assujettie
l/C, Ro et a, sur l’ensemble des n points de mesure aux conditions suivantes :
releves pour differents nombres de Reynolds, la O<a<l; C>O; R,>O (29)
somme S definie ci-dessous :

S= 2 l&h - ( CRe;, w. + R”)]z c21)

1 L
4 n TRAVAIL EXPiRIMENTAL
En procedant aux derivations partielles sui-
vantes :
Les donnees experimentales doivent se prdter a la
a&S &!3 methode d’identification d&rite plus haut. L’objectif
z=o; T&=0 (22)
experimental est de collecter, pour chaque valeur
d’humidite relative de l’air, des mesures pour
on obtient, respectivement : differentes valeurs du nombre de Reynolds tote
externe, done du debit d’air.
n
1"
--- R th 1
c Re”,, WO C c 1 (Re&Wo)*
1
4.1 BANC EXPERIMENTAL
- ROE2 =0 (23)
Re”,, WO Le bane experimental est compose de quatre
1
elements. Le premier est une centrale de traitement
d’air humide qui permet de reguler un debit
2 Rth - i 2 1 - n Ro = 0 (24)
1 1 (Re‘$, WO) constant d’air en temperature et en hygrometrie.
Elle est constituee de deux parties :

cn & hReDe
-- 1 n In ReDe - une partie inferieure qui renferme le groupe

1
c Re& Wo C 1 (Re”,, Wo)* frigorifique et le demineraliseur
l’eau d’alimentation
destine a traiter
de la chaudiere generatrice de
n h ReDe vapeur ;
-RoCRe&,Wo=o (25) - une partie superieure constituee par une cuve
1 etanche qui renferme l’evaporateur du groupe frigo-
rifique, neuf resistances chauffantes partitionndes
Le systeme d’equations (23) a (25) est non- en trois groupes et une rampe d’injection de vapeur ;
lineaire, mais on remarque que ces trois equations
sont lineaires par rapport a (l/C) et RD. Par - un ventilateur entrain& par un moteur Blec-
consequent, la resolution des equations (23) et (24) trique g&e par un convertisseur de frequence qui
selon la methode de Cramer donne : permet un reglage fin du debit.
Le second element est une centrale de traitement
d’eau glycolee qui permet de reguler en temperature
un debit d’eau constant. Elle se compose de deux
(26) appareils mantes en cascade de faGon a assurer une
regulation dteau par melange, a savoir :
- un refroidisseur a evaporateur noye pouvant
reguler une reserve d’eau de 225 L a une tempera-
tures de - 5 “C ;
Ro = - un rechauffeur / thermoregulateur muni dune
vanne de melange et de deux groupes de resistances
~(Re~Wo)~(~)~$(Re~fW~)2~R’h chauffantes de 9 kW permettant d’obtenir des
temperatures allant jusqu’a 80 “C pour un debit
d’eau maximum de 3 m3 .h-’ .
($&)2hn$ (Re~~Wo)’ Le troisieme composant est une veine de mesure
constituee par un caisson cylindrique divise en trois
(27) parties, dont deux elements fixes situ& de part et

293
 R Benelmir. M Khalfi, M Feidt

d’autre de la partie mobile qui renferme l’echangeur TABLEAU I / TABLE I

a Studier. Caractiristiques de I’ichangeur
Le dernier Clement est une centrale d’acquisi- Charecteristics of the exchanger
tion qui permet une lecture directe des donnees Type de tubes tubes lisses en cuivre
experimentales, ainsi que leur enregistrement sur Type d’ailettes ailettes planes en
disquette afin de pouvoir les exploiter ulterieure- aluminium
ment sur un ordinateur. Longueur de l’bchangeur 26,2.10e2 m
L’instrumentation est essentiellement composee Largeur de l’echangeur 8.10F2 m
de debitmetres, de capteurs d’humidite, de capteurs Hauteur de l’echangeur 22,8.10-’ m
de temperature, dun manometre et dune balance Nombre de circuits 9
pour la mesure de la masse de condensat (voir Nombre de rangees 5
figure 5).
Nombre de tubes 45
Diametre externe des tubes 0,97.10e2 m
Diametre interne des tubes 0,87.10e2 m
Longueur des tubes 27,5.10-” m
CIRCULATION D’AIR Pas longitudinal des tubes 1,6.10-” m
VEINE DE MESURE Pas transversal des tubes 2,58.10F2 m
Nombre d’ailettes 164
Pas d’ailettes 1,6.10-” m
Epaisseur d’ailette 0,2.1OF” m
Rapport ~(&+/&6~,~) 16
Rapport s (Xf /6f ) 8
Diametre hydraulique Dh
1 : ventllateur 6 vltesse variable
!- 1,98.10-” m
2 dWtmMn B ‘Wet voftex"
3 I capteun capacttifs d’humldite
et sondes platines Le systeme de regulation de la temperature
4 : thenocouples d’entree de l’eau glycolee engendre une leg&e
5 : balance
6 : d6bifmBtm dlectmmaan&ioua fluctuation de type sinusoidale et de periode
#environ trois a quatre minutes, ce qui nous
a incites a fixer l’intervalle d’enregistrement des
Fig 5. SchCma du bane exp&imental. donnees a deux minutes. Le regime etabli est atteint
Fig 5. Experimental apparatus. au bout de 30 a 60 minutes, suivant la consigne en
humidite. La duree des essais peut aller jusqu’a
deux heures en fonction de la stabilite des mesures.
Les caracteristiques principales de l’echangeur
sont presentees dans le tableau I.

4.3. D~POUILLEMENT DES DONN~ES

4.2. PROCEDURE EXPERIMENTALE
Un logiciel a ete concu sous Excel pour calculer tous
les parametres thermophysiques (annexe B) lies
Pour chaque valeur de consigne de l’humidite aux ecoulements des deux fluides et pour Btablir
relative de l’air a l’entree de l’echangeur, on fait les bilans de masse et d’energie. Les parametres
varier le debit d’air en fixant les parametres calcules sont donnes ci-dessous.
suivants : Pour l’air :
- temperature s&he de l’air a l’entree de l’echan- - temperatures s&he et de resee moyennes ;
geur ;
- pression partielle et de saturation de la vapeur
- temperature de l’eau a l’entree de l’echangeur ;
d’eau presente dans l’air ;
- debit d’eau a l’entree de l’echangeur.
- humidite relative, taux d’humidite et degre de
Les parametres experimentaux sont donnes en saturation de l’air ;
annexe A. Les don&es experimentales enregistrees - enthalpie de l’air humide ;
sont les suivantes :
- temperature s&he de l’air en entree/sortie de - masse volumique, viscosite dynamique,
conductivite thermique et chaleur massique de l’air
l’echangeur ;
set et humide ;
- temperature de resee de l’air en entree/sortie
- debit massique de l’air set et humide ;
de l’echangeur ;
- nombre de Prandtl et nombre de Reynolds ;
- temperature de l’eau en entree/sortie de
l’echangeur ; - masse theorique de condensat ;
- debit d’air ; - puissance thermique transferee.
- debit d’eau ; Pour l’eau glycolee :
- masse eventuelle de condensat. - temperature moyenne ;

294
 Identification du coefficient de transfert thermique

- masse volumique, viscosite dynamique, TABLEAU Ill / TABLE III

conductivite thermique et chaleur massique, debit Corr6lations Expirimentales
massique ; Experimental Correlations
- nombre de Prandtl et nombre de Reynolds ;
Humiditk relative CorAlations
- puissance thermique transferee. 0,304
20 % j = 0,1446 x (Rep,)-
Pour l’echangeur : j = 0,0766 X (ReDh)-0'304
- difference de temperature logarithmique 35 % j = 0,4021 x (RepI)-
moyenne ; j = 0,1644 x (ReDh)po’J2”
- nombre d’unites de transfer-t (NTU) et effica-
40 % j = 0,9015 x (Rep,)-0X524
cite. j = 0,301i’ x (ReDh)-0’524
Pour un essai donne, les points experimentaux
50 % j = 1,5423 x (Repl)-0.“g2
presentant le minimum d&art de puissance entre
j = 0, 447 x (ReDh)-0'5g2
l’eau et l’air sont selectionnes, suivant la procedure
presentee en annexe C, et sont transfer-es vers
un fichier de synthese (sous Excel) qui permet temperature) sont explicit&es en annexe E. L&art
de calculer dune part le coefficient de transfert entre nos correlations (fonctions de l’humidite) et
de chaleur experimental (annexe D), et d’autre les correlations de la litterature est d’autant plus
part les coefficients de correlation a et b, 21 prononce que l’humidite est faible.
partir de la methode presentee precedemment. Les La figure 7 permet de comparer nos correlations
limites de variation du debit d’air, dont la gamme a celles de la litterature.
comprend largement celle utilisee en climatisation
automobile, sont retranscrites en fonction des
differents nombres de Reynolds dans le tableau II.

TABLEAU II / TABLE II 6 n CONCLUSION

Valeurs extrhnes de Re et V,,,
Extreme values for Re et V,,,
Debit d’air humide (m3.he1) ) 1000 1 400 Cette etude nous a permis d’etablir que les condi-
tions hygrometriques de I’air exercent une influence
ReDh 1 1100 1 450 ( non negligeable sur le coefficient de transfert. En
1 Rep1 I8800 I3500 I l’absence de condensation, plus l’humidite de l’air
I5500 I2200 1 est importante, moins le transfert thermique est
bon. Cette conclusion devrait avoir un impact impor-
tant sur le dimensionnement d’echangeurs destines
a fonctionner dans des regions a forte humidite.
Les erreurs systematiques sur les mesures sont
de f0,2 “C pour les capteurs de temperatures et
de k 2 % de l’echelle de mesure pour les debits. ANNEXES
L’incertitude sur le facteur de Colburn calcule
experimentalement est de 10 % (Moffat, 1988).
A. Paramdres expkrimentaux

Les parametres correspondants aux quatre essais

men& sont les suivants :
5I RbULTATS Essai no1 :
T a,e =30xt1°C
T 20,e = 2ozk 1 “C
Des correlations ont ete developpees en fonction de = 3f0,l m3.hp1
l’humidite relative (tableau III) avec un nombre de $” =2OItl%
azT = 1000 - 900 - 800 - 700 - 600 - 500 - 400
Reynolds rapport6 au pas longitudinal des tubes
(13) m”.h-’
et au diametre hydraulique, afin de pouvoir les
comparer a celles de la litterature. Essais no 2, 3 et 4 :
Ces correlations sont representees sur la figure 6. T a,e =4Oil”C
Le facteur de Colburn decroit avec l’humidite T w,e = 25 f 1 “C
relative de l’air. Cet effet est accent& pour les Qu, = 3 iO.l m3.h-’
nombres de Reynolds importants. L’influence de
l’humidite sur le coefficient d’echange ainsi que avec pour l’essai no 2 :
la coherence des mesures experimentales (relation =35&l%
entre le coefficient d’echange, l’efficacite d’ailette, azr = 1000 - 900 - 800 - 700 - 600 - 430
la resistance thermique de l’air et l’ecart moyen de (i- 3) m3.h-’

295
 R Benelmir, M Khalfi, M Feidt

1,4E-02
1,3E-O2
1,2E-O2
l,lE-02

i l,OE-02
9,OE-03
6,OE-03

7,OE-03
I
6,OE-03 -i /
400 600 600 1000 1200

Re (Dh)
Fig 6. j en fonction de ReDh.
Fig 6. j with respect to ReDh.

0,014 lyan:r=l6
0,013 ,
0,012
0,011
0,010
j 0,009
0,006
0,007
0,006
0,005
0,004
400 600 600 1000 1200 1400 1600
Re (Dh)
Fig 7. Comparaison avec /es corrilations d’autres auteurs
Fig 7. Comparaison with other correlations.

avec pour l’essai no 3 : oti Cp, et Cp, sont les chaleurs massiques de l’air
=40&l% set et de la vapeur d’eau.
;.+ = ;:"$ ,:o,o-; 800- 700 - 600-500-400 Vzkcosite dynamique (Reid, 1988) :

avec pour l’essai no 4 :

pcLh(Pa.s)= yyiFy + .'i",;, (3b)
=50&l% 2' a a 1'
$ air = 1000 - 900 - 800 - 700 - 600 - 500 - 400 avec :
(9~3) m3.he1

6. ParamGtres thermophysiques de I’air

humide
Masse volumique :
(1 + W)Ptot et
Pah (kg.m-3) = (lb)
461,51 (w + 0,622) Ta
avec : w=11411
ma
Chaleur massique : oii pa, pLvet Y,, Y, sont respectivement les viscositks
dynamiques et les titres molaires de l’air set et de
Cpah([Link]-l.K-I) = Cp, +wCp, (2b) la vapeur d’eau.

296
 Identification du coefficient de transfert thermique

42
40
p 38
- 36
3 34
$ 32
E 30
'g 28
E 26
2 24

2 8 14 20 26 32 38 44 50 56 62 68 74 80
Temps (mn)
Fig 8. ivolution des tempe’ratures.
Fig 8. Temperatures profiles.

Conductivit6 thermique de l’air humide (Reid, L’ecart relatif entre ces deux puissances s’ex-
1988) : prime par :

X,h(W.m-‘.K-‘) = r,x’ Y, xa P - peau

(4b) Ep = 100 azrp
Y,, + CYa + YaSdY, eau
avec :
et est represent6 sur la figure 9. Les points
C=
selectionnes doivent satisfaire la condition :

0,25{1+
[p, (~)“~75=&]o~5}‘~ Ep<5% (3c)

Pour les 40 points de mesure, on a :

d= 95 % des points de mesure avec Ep < 10 %

0,25{1+
[(~)(gj”~75E&]o~5}2~ 70 % des points de mesure avec Ep
17 % des points de mesure avec Ep
< 5 %
-c 1 %.
Le tableau IV presente les resultats de l’analyse
;I,, = 0,76 dm statistique pour cet essai.
;r, = 1,5T,, et A, = 1,5Tea

et T,, , T,, et A,, A, sont respectivement les tempera- D. Cal4 du coefficient expkrimental
tures d’ebullition et les conductivites thermiques de de transfert de chaleur
l’air et de la vapeur d’eau.
Le calcul de la puissance transferee par l’eau P,,
est plus fiable que celui de la puissance transferee
C. Wection des points expkimentaux

Afln d’illustrer le mode de selection des points

experimentaux, nous presentons dans cette annexe
la procedure correspondant a un essai particulier,
dune duree de 80 minutes (ce qui correspond a
40 points de mesure).
La figure 8 represente l’evolution des tempe-
ratures. Les fluctuations des temperatures dent&e/
sortie d’eau ont et6 expliquees precedemment. Le
systeme de commande (Blectrovannes et rampe d’in-
jection de vapeur) et de regulation (proportion-
nelle) de l’humidite de l’air entraine les fluctuations
observees au niveau des temperatures de resee. 0 20 Tem$(mn) 60 80
La selection des points est basee principalement
sur la precision des bilans thermiques. La puissance
transferee par l’air est calculee comme suit : Fig 9. ReprCsentation de l’erreur relative dans /es bilans
thermiques.
Pa =&Ah (lc) Fig 9. Relative error on heat transfer balances.

297
 R Benelmir. M Khalfi. M Feidt

avec l’efficacite de la surface ailettee 71 et de l’ailette

E donnees par les relations suivantes :

(8d)

et ou le parametre m est defini par :

Clod)

Comme l’efficacite de la surface ailettee depend

du coefficient d’echange convectif, le calcul de h, et
de 17se fait par iteration.

par l’air, en raison dune dispersion plus importante

quant a la mesure de l’humidite de I’air. On estime E. CohCrence des mesures
la resistance thermique globale Rt,, a partir de la expkrimentales
relation suivante (en supposant que l’echangeur est
adiabatique) :
Rth = AZ,
~ (Id)
Pu,
otl : En general :
Pw = ti,,. Cp, m AT,, (2d)
Nu=nRebPr”” avec Ocb-cl (2e)
La resistance thermique du tote externe (tote
air) peut etre calculee comme suit : d’oh, le facteur de Colburn devient :
R th, = Rth - RttL, ~ c R (3d) j = (Lfj‘ehpl (3e)
La resistance thermique du cot6 interne est
definie par : La relation (3e) fait apparaitre clairement que
le facteur de Colburn decroit avec le nombre de
R (4d) Reynolds. Par ailleurs, l’expression (le) permet de
constater que pour une nombre de Reynolds don&,
et, par definition du nombre de Nusselt : le facteur de Colburn j varie dans le m6me sens
que le nombre de Nusselt et done que le coefficient
(5(i) convectif h,.
Les figures 10 et 11 representent l’evolution du
Ce nombre de Nusselt est calcule avec la correla-
coefficient d’echange experimental et de l’efficacite
tion de Dittus-Belter (1930) pour un ecoulement
de la surface ailettee qui d’ailleurs en depend
en regime turbulent dans un tube lisse en mode de
implicitement par le biais de l’efficacite d’ailette.
chauffage :
11 en resulte que pour une humidite relative plus
Nu, = 0.023 (Re”.8 PT’.~)~ (6d) importante l’efficacite de la surface ailettee est
augmentee alors que le coefficient d’echange est
sous les conditions suivantes : moindre.
0,7 < Pr < 160 Le produit de l’efficacite d’ailette par le coeffi-
Re > 10000 cient d’echange decroit avec l’humidite (fig 12), ce
L/D > 60 qui correspond a une augmentation de la resistance
thermique de l’air. Cette consequence est aussi vi-
On calcule ainsi la resistance thermique interne, sible a travers la figure 13 qui montre que l’aug-
RttL.%, puis externe, Rth,e, et on en deduit le mentation de cette resistance se traduit par une
coefficient de transfer% de chaleur h, : augmentation de l’ecart de temperature moyen, la
puissance transferee &ant sensiblement constante.

298
 Identification du coefficient de transfert thermique

0,73
0,71
0,69
3 0,67
‘3 0,65
g 0,63
w 0,61
0,59
0,57
0,55
400 600 800
Re(Dh)
Fig 10. Efficacit6 q en fonction de Reoh.
Fig 10. Efficiency 71 with respect to Reoh.

r: 70
60
50
400 600 800
Re(Dh)
Fig 11. Coefficient d’tkhange en fonction de ReDh.
Fig 11. Transfer coefficient with respect to Reoh.

70
65

60
55

50
45

40
400 600 800
Re(Dh)
Fig 12. Produit h 77 en fonction de ReDh.
Fig 12. Product h 7 with respect to ReDp,.
 R Benelmir, M Khalfi, M Feidt

400 600 800 1000 1200

Re(Dh)
Fig 13. ATl, et P, en fonction de Reoh.
Fig 13. ATi, and Pa with respect to Reoh.

REFERENCES L’H6te P, Bardon JP (1989) Mesure des resistances de

contact ailette tubes dans un echangeur compact,
Application au controle non destructif. Rev Ce’n Therm
334/l 0
Canton P, Fornasieri E (1983) About the thermal transmit-
tance on the air-side of finned coils (Sulla trasmittanza MC Quiston FC (1978) Correlation of heat, mass and
termica lato aria nelle batterie alettate). Condiziona- momentum transport coefficients for plate fin-tube
mento dell’aria, Riscaldamento, Refrigerazione 9, 933- heat transfer surface with staggered tubes. ASHRAE
946 Trans 2487, 294-309
Chen ZQ, Ren JX (1988) Effect of fin spacing on the heat Moffat RJ (1988) Describing the Uncertainties in Experi-
transfer and pressure drop of a two-row plate fin and mental Results. Exp Therm fluid Sci, 3-l 7
tube heat exchanger. Int 1 Refrig 11, 356-360
Reid RC, Prausnitz JM (1988) The properties of gas and
Dittus FW, Bcelter LMK (1930) Univ Calif, Berkeley, Pub/ liquid. MC Craw Hill
Eng 2, 443
Rich DR (1973) The effect of fin spacing on the heat
El Mahdy AH, Briggs RC (1979) Finned tube heat
transfer and friction performance of multi-rows, smooth
exchangers: correlation of dry surface heat transfer
plate fin-and-tube heat exchangers. In : Proceedings of
data. ASHRAE Trans 85, Part 2, 262-273
the ASHRAE Annual Meeting, Louisville KY, 137-l 45
ldem SA, Jacobi AM, Coldschmidt VW (1990) Heat transfer
Shah RK (1990) Assessement of modified Wilson plot
characterization of a finned-tube heat exchanger (with
techniques for obtening heat exchanger design data.
and without condensation). j Heat Transfer 112,64-70
Heat Transfer Eng 5, 5 l-56
Kayansayan N (1994) Heat transfer characterization of
plate fin- tube heat exchanger. Rev Int Froid 17, 49-57 Turaga M, Lin S, Fazio PP (1966) Performance of direct
expansion plate finned tube coil for air cooling and
Khartabil HF, Christensen RN (1992) An improved sheme
dehumidification. Int / Refrig 1 1, 78-86
for determining heat transfer correlation from heat
exchanger regression models with three unknowns. Wilson EE (1915) A basis for rational Design of heat
Exp Therm Fluid Sci 5, 808-819 transfer apparatus. ASME Trans 2, 47-70

ABRIDGED ENGLISH VERSION

Identification of the heat transfer coefficient over the external area of a finned-tubes
heat exchanger with respect to the moisture content of the air without condensation

The use of compact heat exchangers is wide In the course of the last 20 years several studies
spread in the field of heating, ventilation and air have been carried out on this topic, the major part
conditioning. The complex nature of the air flow over of them being concerned with the impact of key
the external area of this latter would not actually geometry paramaters, as the fin spacing, on the heat
allow the development of an analytical expression transfer coeficient. Calculation procedures were
for the heat transfer coefficient. Therefore, we have usually based on the original Wilson method (1915)
to carry out experimental work in order to build or on alternative methods (Shah, 1990; Khartabil,
empirical correlations. 1992).

300
 Identification du coefficient de transfert thermique

The present study deals with the identification of measurement of the inlet temperatures and humidity
the heat transfer coefficient for an air cooler with (for air), mass flow rate, mass (for the condensed
respect to the moisture content of the air in the water uapor) as displayed in figure 5. The main
absence of condensation. difficulty in this experimental work concerns the
In the first part of this paper we compare the measurement and the control of the air humidity,
various correlations and show, as displayed in since one has to make sure that the measurement
figure 4, that divergences still exist regarding the incertainties are not significant when compared to
impact of the fin spacing on the heat transfer the precision of the calculation of the heat transfer
coefficient. Kayansayan’s relation for the Colburn coefficient.
factor j is plotted with respect to Reynolds number Correlations with respect to the relative humidity
Rem, in figure 3. In this approach, the ratio r have been set for Reynolds numbers based on the
between the total area on which heat transfer occurs hydraulic diameter and the original spacing of the
and the total area inside the tube was considered tubes in order to make easier the comparison with
as a parameter; a significant decrease of j with r correlations by other authors (table III), reported in
was observed. Moreover, we note the lack of studies figure 6. We can see that the Colburn factor decreases
about the impact of air humidity on the heat transfer with respect to the relative humidity of the air and
coefficient. that this trend is more significant for high Reynolds
Our identification procedure of the heat transfer numbers. Figure 7 illustrates that the difference
coefficient is based on Khartabil’s method (1992) between our correlations, which are a function of the
which originates, as is the case for most of the other humidity, and those from the literature is higher at
methods, from the general relation expressing the low humidity values.
overall thermal resistance of the heat exchanger This work has shown that the moisture content in
The experimental apparatus consists of a heat the air has a significant impact on the heat transfer
exchanger of the same type as those used in process. In the absence of condensation, the heat
automotive air conditioning (horizontal copper tubes transfer coefficient decreases with the humidity of
and plane aluminium fins) where the external fluid the air. This conclusion is likely to have an important
is moist air (with glycol) - which exhibits well- impact on the design of heat exchangers for humid
known thermal properties inside circular tubes -, areas (ouersizing or improvement of the heat transfer
equipped with all the means for control and/or areas).