INT. COMM. HEAT MASS TRANSFER 0735-1933/91 $3.00+ .

00
VoL 18, pp. 731-740, 1991 Printed intheUnitedStates
ePergamon Press plc

COMPORTEMENT EN REGIME TRANSITOZRE DE DZVRR8 TYPES D'ECHANGEUR8

DE CHALEUR ; MODELIBATION ET CONSEQUENCES

M.Henrion , M.Feidt
G.E.S.P.E - L.E.M.T.A , U.N 1 et I.N.P.L
2 AV. de la for~t de Haye 54504 Vandoeuvre les Nancy Cedex

(Communicated by J. Gosse)

ABSTRACT
This paper contributes to develop simple analysis of
transient behaviour of inverse cycle machines , to
say heat pumps , refrigerating machines . We report
here of an analytical procedure to caracterise
transient behaviour of heat exchangers . Single phase
heat transfer extension of analytical results are
proposed for tubular and "plates" heat exchangers ,
on the basis of recent published works . Examination
and comparison with experimental data on a tubular
heat exchanger conclude to the present limitations of
these models .

Introduction
L'etude du comportement en regime transitoire d'echangeur
peut se faire par integration de syst~mes diff~rentiels selon
des methodes numeriques discretes , auquel cas les calculs sont,
longs et couteux [1] . Depuis de nombreuses annees [2]
toutefois des propositions visant & la mise en place de modeles
plus simples d'emploi et plus rapides d'execution ont v u l e
jour afin de r~pondre aux besoins de r~gulation et
d'automatisation des processus thermiques . J.Fontaine et F.
Lannoy [3] sont parmi les premiers a avoir propose un module
empirique pour les ~changeurs compact h plaques sur la base

731
732 M. Henrion and M. Feidt Vol. 18, No. 5

d'experiences hydrodynamiques en e c o u l e m e n t turbulent

monophasique . Y.A. Gogus et O.E. Ataer [4] ont p a r c o n t r e
examine s u r un m o d e l e plus simple que le p r e c e d e n t , le
comportement en t r a n s i t o i r e d'une batterie ailettee , composant
essentiel des climatisations . Leurs hypotheses sont tres
restrictives ; en p a r t i c u l i e r :
- variations lineaires des t e m p e r a t u r e s
- m a s s e des frontieres negligeables
Le m o d u l e fait intervenir une co,stanCe de t e m p s et un
temps de retard .
Cette m~me approche se r e t r o u v e d a n s un a r t i c l e recent de
D. A z i l i n o n , P. P i e r s o n et J. Padet [5] . Mais contrairement
l'article precedent , ces a u t e u r s considerent & la fois les
echelons de t e m p e r a t u r e et les e c h e l o n s de d e b i t quelque soit
le r e g i m e d'ecoulement des fluides ( laminaire ou t u r b u l e n t ) ,
ou la c o n f i g u r a t i o n de l ' e c h a n g e u r bitube considere ( co ou
contre courant ) . L'originalite de l ' a r t i c l e provient
essentiellement de l ' e x t e n s i o n de la m e t h o d e aux echangeurs a
faisceaux de t u b e s , moyennant des c r i t e r e s d'equivalence
La d e m a r c h e tres judicieuse suivie dans [5] est r e p r i s e
ici m a i s en m o d i f i a n t l'une des c o n d i t i o n s d'equivalence
Alors que [5] s u p p o s e la c o n s e r v a t i o n du c o e f f i c i e n t d'echange
s u r un d e s fluides , nous p r e c o n i s o n s l'hypothese d'egalite des
coefficients d'echange globaux des e c h a n g e u r s De ce fait ,
l'echangeur reel et le b i t u b e fictif equivalent ont les m ~ m e s
caracteristiques et e c h a n g e n t donc la m ~ m e puissance.

Condi~o~$ d,6auivalence hun 6ehanaeur b i t u b e

Ii s ' a g i t de d e t e r m i n e r les 5 i n c o n n u e s d i m e n s i o n n e l l e s de
l ' e c h a n g e u r e q u i v a l e n t ( 4 r a y o n s , la l o n g u e u r ) a l ' a i d e de 5
e q u a t i o n s t r a d u i s a n t la c o n s e r v a t i o n d e s c a p a c i t e s t h e r m i q u e s
des d e u x e n v e l o p p e s , c e l l e d e s f l u i d e s c h a u d et f r o i d et ici
la c o n s e r v a t i o n d u e o e f f i o i e n t d,6ohange global

Echanqeur a Dlaaues

Les s e u l e s p r e c a u t i o n s a p r e n d r e c o n c e r n e n t le c h o i x d e s
c o r r e l a t i o n s de t r a n s f e r t , a i n s i q u e la d i f f e r e n t i a t i o n e n t r e
Vol. 18, No. 5 COMPORTEMENT TRANSITOIRE D'ECHANGEURS DE CHALEUR 733

les p l a q u e s i n t e r m e d i a i r e s et les p l a q u e s d ' e x t r ~ m i t ~ s

e q u i v a l e n t e s r e s p e c t i v e m e n t au tube i n t ~ r i e u r et & la c a l a n d r e
du b i t u b e ~ q u i v a l e n t .

Echanaeur & f a i s c e a u de tubes

Dans ce cas , le r a i s o n n e m e n t est a n a l o g u e a celui

p r ~ s e n t ~ p a r D. A z i l i n o n , P. P i e r s o n et J. Padet [5] .
Toutefois l'~quivalence du c o e f f i c i e n t de t r a n s f e r t global se
traduit par l'~quation suivante :

1
- - = A L -1"4 + B L -I
KS

A et B , f o n c t i o n s c o n n u e s des c o e f f i c i e n t s CI,C2,C3 [6] .

La r e s o l u t i o n n u m ~ r i q u e de c e t t e e q u a t i o n non l i n e a i r e fournit
L la l o n g u e u r du b i t u b e e q u i v a l e n t ; d ' o % sont e n s u i t e d ~ d u i t s
les quatre rayons s i m p l e m e n t [6]

Dimensionnement du b i t u b e e u u i v a l e n t

L ' e g a l i t ~ des c a p a c i t e s c a l o r i f i q u e s des q u a t r e c o m p o s a n t s

de l ' e c h a n g e u r c o m p a c t implique q u a t r e e q u a t i o n s par type
d ' e c h a n g e u r entre 5 inconnues ( r I , r 2 , r 3 , r 4 , L )

Echangeur bitube plaques a f a i s c e a u de tubes

~*L*(r32-r22) = (Npiyl)/2*ipi*L~i*Efc ou= ~ * L t * ( R 3 2 - N t * R 2 2 )

~*L*(rl 2 ) = (Npi-l)/2*ipi*Lpi*Eft ou= ~ * L t * ( N t * R I 2)

~ * L * ( r 2 2 - r l 2) = Npi*ipi*Lpi*Epi ou= ~ * L t * N t ( R 2 2 - R I 2)

~*L*(r42-r32) = 2*ipe*Lpe*Epe ou = ~ * L t * ( R 4 2 - R 3 2 )

La c i n q u i e m e e q u a t i o n e x p r i m e l ' e g a l i t e de la r e s i s t a n c e
thermique globale

1 1 (r2-rl) 1
+ +
KS h c (2~r2L) Ametal(r2+rl)/2*2~L ht (2~rlL)

La r e s i s t a n c e de c o n d u c t i o n c o r r e s p o n d & l ' a p p r o x i m a t i o n
des faibles e p a i s s e u r s . Ce qui est g ~ n ~ r a l e m e n t r e c e v a b l e et
c o n d u i t ~ une s i m p l i f i c a t i o n des c a l c u l s .
Le s y s t e m e non l i n e a i r e r e s o l u p a r a p p r o x i m a t i o n s s u c c e s s i v e s
avec c o m m e c r i t ~ r e de c o n v e r g e n c e l ' ~ g a l i t ~ des c o e f f i c i e n t s
d'echange globaux echangeur compact - echangeur bitube
equivalent .
734 M. Hem-ion and M. Feidt Vol. 18, No. 5

180
,~' 170
160
150
140
150
120
110
1O0
go
J3 80
~ 7o
~ ~o I
m 5o
~ ,o +
~ '°
~o I
I
o , , | ~ , , , , , , , ,
0 2~0 40 60 80 100 120 140 160
KS~IObGI #chCr~eur a fOi$¢ilOu [kW/K]

Fig.l - Ddtermination de la l o n g u e u r du b i t u b e equivalent

0.14

012

J~ Ol

OOg

~L 00s " extlr feur

g o05
0.04

o 002
eur tube i n ~ ; ~ , u r '
o01
10 30 50 70 gO ] 10 130 150
LO.~Ueu" Oitut~e ~ i v a l e m t [m]

Fig.2 - Dimensionnement du b i t u b e d q ~ a i v a l e n t en fonction de sa

longueur

Les r e s u l t a t s p r e c e d e n t s p e u v e n t ~tre t r a d u i t s sous forme

d ' a b a q u e , t r a d u i s a n t la m d t h o d e p r ~ c e d e n t e .

Les f i g u r e s 1 et 2 r e n d e n t c o m p t e sur un e c h a n g e u r
f a i s c e a u de t u b e s de l ' i n f l u e n c e du c o e f f i c i e n t de t r a n s f e r t
g l o b a l sur les c a r a c t d r i s t i q u e s d i m e n s i o n n e l l e s de l ' e c h a n g e u r
dquivalent

La figure 3 r e p r d s e n t e 1 ' d v o l u t i o n de la c o n s t a n t e de
t e m p s c a l c u l d e p o u r d i v e r s e s v a l e u r s des d e b i t s en p a r t i c u l i e r
d ' e a u en t u b e s ( le sens de v a r i a t i o n est p l u s c o h d r e n t d a n s ce
casque d a n s celui p r o p o s ~ p a r D . A z i l i n o n ) . I1 r e s s o r t des
c a l c u l s u n e a u g m e n t a t i o n s e n s i b l e de la c o n s t a n t e de temps par
VoL 18, No. 5 COMPORTEMENT TRANS1TOIRE D ' E C H A N G E U R S DE C H A L E U R 735

r a p p o r t & la m e t h o d e precedente , ce q u i v a d a n s le s e n s d e
l'experience .

14

13

12
t'

7 -i "-- ~'~'i'"-~.~=.l-__ , 115 m31h

6 ~ ~ + . . . . + . . . . +----_ "'''~I" '"v'''"

5 [ I i I +I - - I- - - +
40 60 80 100
Debit eou colQndre [ m 3 / h ]

Fig.3 - Comparaison des constantes de temps

= Azilinon et al [5] ; --+--Presents calculs

D6termination exp6rimentale de la constante de temps

P~esentation de l'echangeur etudie

L ' e c h a n g e u r , d e m a r q u e P i c k e r , d u t y p e A. I I 0 . G . 1 9 . C ,
e s t u n e c h a n g e u r c o m p o s e d ' u n e c a l a n d r e e n a l u m i n i u m d e 125 m m
d e d i a m e t r e e x t ~ r i e u r , II0 m m d e d i a m e t r e i n t e r i e u r et d ' u n
f a i s c e a u d e ll0 t u b e s e n c u i v r e d e 8 m m d e d i a m e t r e e x t e r i e u r ,
7 m m d e d i a m e t r e i n t e r i e u r et 600 m m d e l o n g . Le f a i s c e a u e s t
& pas triangulaire , l'espace entre tubes est de 2 ~m , ce qui
donne une surface d'echange interne tubes de 1.45 m ~ .
9 c h i c a n e s e s p a c e e s d e 60 m m m o d i f i e n t la d i r e c t i o n d e p a s s a g e
du fluide calandre .
736 M. Henrion and M. Feidt Vol. 18, No. 5

ExPloitation d'un essai en t e m p e r a t u r e et p u i s s a n c e

ESSAI2H Qft= 10. lm3/h Qfc=3.95m3/h

7O

60 ~ ~

" 50

• $oct|e Clllllrl~ll I n tllllllOlrlllture

30

• sortie t ~ en temperature

~t 0 T--'----r----~--'----40
30
wr-----JI trct f ~ .i. retard f[uidt c a & w t d r e eft t e t t t r o t u r e

Fig.4 - Evolution des temperatures des fluides au cours du

temps ( en s )

ESSAI2H Qft=iO. lm3/h Qfc=3.95m3/h

2OO
~90 I~

140
~30 t lut~
~ 120
o
~.. I 10 . . . . . . .

I O0 -- - w - - . - - . , . , p ~
90 "tp
80
7O
~ tUl~m on puis~e
} 5o

c
~ 3o
5 20 ]j
3. 10 rmlm ~ nltll~l fLui6l tubes an l~il.IIl~e
o a i ' i " i i ' i
10 20 30 40
trtp
TempsenseconOes
Fig.5 - Evolution de la puissance thermlque ~cbangee
VoL 18, No. 5 COMPORTEMENT TRANSITOIRE D'ECHANGEURS DE CHALEUR 737

La f i g u r e 4 p r 6 s e n t e l ' @ v o l u t i o n des t e m p @ r a t u r e s de
l ' @ c h a n g e u r s o u m i s ~ un "~chelon" en e n t r @ e de c a l a n d r e ; On en
d @ d u i t les d i v e r s e s c o n s t a n t e s de temps et t e m p s de retard.

La f i g u r e 5 p r @ s e n t e l ' @ v o l u t i o n des p u i s s a n c e s @ c h a n g ~ e s
d a n s le m @ m e cas de figure ; on en d @ d u i t de n o u v e l l e s
c o n s t a n t e s de t e m p s et t e m p s de r e t a r d .

D@Douillement des essais

Le p l a n d ' e x p @ r i e n c e s a d o p t @ p o u r la c a m p a g n e de m e s u r e a
p e r m i s d ' @ t u d i e r le c o m p o r t e m e n t de l ' @ c h a n g e u r d a n s un m a x i m u m
de cas de figures que l'on p e u t r e n c o n t r e r d a n s l ' i n d u s t r i e :

E c o u l e m e n t & co et c o n t r e c o u r a n t
D @ b i t fluide c a l a n d r e m a x i et mini
D@bit t u b e s m a x i et m i n i
E c h e l o n de t e m p @ r a t u r e sur le fluide chaud de forte et
faible a m p l i t u d e

Ceci a a m e n @ & faire 16 essais sur cet e c h a n g e u r .

L ' e s s a i 2H p r @ s e n t e c i - a v a n t est un essai a @ c o u l e m e n t de
fluides c o n t r e - c o u r a n t , d 4 ~ i t c a l a n d r e m a x i m u m 3,95 m 3 / h ,
d @ b i t t u b e s m a x i m u m i0,i0 m ~ / h en i m p o s a n t un @ c h e l o n de forte
a m p l i t u d e de t e m p @ r a t u r e sur l ' e n t r ~ e du f l u i d e c h a u d en
calandre .
L ' e n r e g i s t r e m e n t g r a p h i q u e d o n n e l ' @ v o l u t i o n des t e m p @ r a t u r e s
d ' e n t r @ e et de s o r t i e des fluides c a l a n d r e et t u b e s (Fig.4) ;
c o n n a i s s a n t les d @ b i t s m a s s i q u e s des f l u i d e s , il est p o s s i b l e
de c a l c u l e r l ' @ v o l u t i o n de la p u i s s a n c e i n s t a n t a n @ e @ c h a n g @ e
p a r c h a q u e fluide en f o n c t i o n du t e m p s (Fig.5)
L ' e x p l o i t a t i o n g r a p h i q u e des r e s u l t a t s o b t e n u s p e r m e t de
d @ t e r m i n e r p a r la sous t a n g e n t e au p o i n t d ' i n f l e x i o n de la
c o u r b e @ t u d i ~ e la c o n s t a n t e de temps ~ qui p e r m e t de
c a r a c t @ r i s e r la c o u r b e par une @ q u a t i o n de t y p e :

t-tr
T(t) = T~ + ( T O - T ~ ) exp - ( - - ) en t e m p e r a t u r e
T

t-tr
p(t) = p . + ( P0- P-) exp- ( - ) en p u i s s a n c e

Le temps de r e t a r d tr est d o n n @ p a r le t e m p s c o m p r i s e n t r e le
d @ b u t du signal d ' e n t r @ e et l ' i n t e r s e c t i o n de la sous t a n g e n t e
a v e c la t e m p @ r a t u r e ou la p u i s s a n c e i n i t i a l e

Les r @ s u l t a t s de l ' e s s a i d @ c r i t s o n t les s u i v a n t s :

P u i s s a n c e m o y e n n e @ c h a n g @ e en r ~ g i m e p e r m a n e n t : 102.26 kW
E c a r t de t e m p e r a t u r e m o y e n l o g a r i t h m i q u e : 37.32 °C
Coefficient d'~change global exp@rimental : 2.81 k W / ° C
T s o r t i e f l u i d e c a l a n d r e en t e m p @ r a t u r e : 5.36 s
s o r t i e f l u i d e c a l a n d r e en p u i s s a n c e : 3.62 s
738 M. Hera'ion and M. Feidt Vol. 18, No. 5

T sortie fluide tubes en temperature 6.91 s

T sortie fluide tubes en puissance 6.71 s
T e m p s de retard calandre en t e m p e r a t u r e 2.62 s
Temps de retard tubes en temperature 0.83 s
Temps de retard calandre en p u i s s a n c e 3.93 s
T e m p s de retard tubes en puissance 1.14 s

R e s u l t a t s du m o d u l e
R e y n o l d s c a l c u l ~ tubes 4257
Reynolds calcul~ calandre 9043
C o e f f i c i e n t d ' ~ c h a n g e global c a l c u l e 3.91 kW/°C
C o n s t a n t e de temps u n i q u e e c h a n g e u r c a l c u l e e : 1.27 s

Conclusion

L ' e x p l o r a t i o n p r o p o s e e tant t h e o r i q u e q u ' e x p e r i m e n t a l e n'a

pas p o u r but de f o u r n i r des r e s u l t a t s p r e c i s , m a i s de d ~ g a g e r
les t e n d a n c e s quant au f o n c t i o n n e m e n t en r ~ g i m e t r a n s i t o i r e
d ' e c h a n g e u r s d i v e r s dans des c o n d i t i o n s i n d u s t r i e l l e s . En ce
sens , elle p r o l o n g e et r e c o n s i d ~ r e les t r a v a u x p r o p o s e s dans
la b i b l i o g r a p h i e .

ADDroche theoriuue

A i n s i , les m o d e l e s sont a p p l i c a b l e s n o n s e u l e m e n t aux

e c h a n g e u r s a f a i s c e a u t u b u l a i r e s , m a i s aussi aux e c h a n g e u r s
p l a q u e s , aussi b i e n q u ' a u x b a t t e r i e s a i l e t ~ e s
L ' ~ c h a n g e u r reel est t o u j o u r s r e m p l a c e p a r un e c h a n g e u r
fictif b i t u b e e q u i v a l e n t c o n f o r m e m e n t & D . A z i l i n o n [5] , mais
ces e c h a n g e u r s ont les m ~ m e s c a r a c t e r i s t i q u e s t h e r m i q u e s et
e c h a n g e n t aussi la m ~ m e p u i s s a n c e

Approche experimentale

Le p l a n d ' e x p e r i e n c e s u c c i n t r e a l i s e sur un e c h a n g e u r
f a i s c e a u de tubes en e c o u l e m e n t s m o n o p h a s i q u e s t u r b u l e n t s
soumis a un e c h e l o n de t e m p e r a t u r e m o n t r e que m a n i f e s t e m e n t ,
il n'y a pas une c o n s t a n t e de temps u n i q u e p e r m e t t a n t de
c a r a c t e r i s e r un e c h a n g e u r , m a i s q u a t r e :
- C o n s t a n t e de temps t e m p e r a t u r e s o r t i e c a l a n d r e
- C o n s t a n t e de temps t e m p e r a t u r e s o r t i e tubes
- C o n s t a n t e de t e m p s p u i s s a n c e ~ c h a n g e e c a l a n d r e
- C o n s t a n t e de temps p u i s s a n c e e c h a n g e e tubes

En fait la c o n s t a n t e de temps t e m p e r a t u r e s o r t i e tube est

identique ,dans le cas ~ t u d i e , a la c o n s t a n t e de temps
puissance echangee cote tubes .
De m ~ m e , le temps de r e t a r d n ' e s t p a s u n i q u e p o u r tout
l ' e c h a n g e u r . S e l o n les e x p e r i e n c e s , il p e u t r e p r e s e n t e r
j u s q u e 70% de la c o n s t a n t e de temps en p u i s s a n c e c S t e c a l a n d r e
, alors que c o t e tubes , il ne r e p r e s e n t e au m a x i m u m que 25% .
VoL 18, No. 5 COMPORTEMENT TRANSITOIRE D'ECHANOEURS DE CHAI.F~UR 739

A u c u n e c o r r e l a t i o n n ' e s t ~ v i d e n t e p o u r le calcul des t e m p s

de r e t a r d , seuls les essais p e u v e n t en f o u r n i r la v a l e u r
actuellement .

E n r ~ s u m ~ , la m o d ~ l i s a t i o n d ' u n e c h a n g e u r en r e g i m e
t r a n s i t o i r e d o i t i m p e r a t i v e m e n t t e n i r compte du t y p e
d ' u t i l i s a t i o n e n v i s a g ~ : l ' o b j e c t i f est il de r ~ g u l e r une
t e m p e r a t u r e de s o r t i e ( o b j e c t i f i n t e n s i f ) ou d ' ~ v a c u e r une
c e r t a i n e p u i s s a n c e ( o b j e c t i f e x t e n s i f ) et ce en r e p o n s e a
q u e l l e s o l l i c i t a t i o n ( ~ c h e l o n de t e m p e r a t u r e , de d e b i t .... )?

Nomenclature

Efc : Epeisse~r du f i l m fLu|de chsud ~changeur h plaques

Eft : Epe|sseur du f i l m ftu~de froJd 6changeur i, plaques

Epe : EpaJsseur plaques d ' e x t r k m i t t ~changeur ~ plaques

Ep| : Epe|sseur plaques |nterm~iJaire ~r~hangeur & plaques

hc : Coefficient de convect|on cot~ catsndre b|tube

ht : Coeff|cient de corwect|on cot6 tube bitube

KS : Coefficient de t r a n s f e r t gLobaL

L : Longueur du b|tube ¢~Jivotent

tpe : Largeur plaques d'extr(mit~ kohangeur ~ plaques

Lpe : Longueur plaques d ' e x t r ( m i t 6 (~changeur ~ plaques

tp| : Largeur plaques i n t e r m ~ i a t r e ~;changeur ~ plaques

Lpi : Longueun plaques interm~liaire 6changeur ~ plaques

Lt : Longueur de tubes de'tW~,changeur J, faisceau de tubes

Npi : N~e de plaques Interm~l|aire ~Lcha~eur ~ plaques

Ht : H~nbre de tubes 6changeur ~ faJsceau de tubes

P(t) : Pu|ssance |nstantann~e ~,chang~,e en fonot|on du temps

PO : Puissance |nstantanrt6e 6chang~e | n | t i a t e

Pet : Pu|ssance instantannle Ichang6e f i n a l e

Qft : O~q)it ftuide tubes

Qfc : DLrqo|t f t u i d e catandre

r1 : Rayon |nt~rieur tube int6r|eur bitube ~lu|vatent

R1 : Rayon i n t l r i e u r tubes ~;changeur ~ fa|sceaux

r2 : RaYon ext6r|eur tube |nt6rieur bitube ~lUivstent

R2 : Rayon ext6rieur tubes (,changeur ~ faisceaux

R] : Rayon |nt6r|eur catandre {,cheng~Jr ~ faJscesux

r3 : Rayon tnt~r|eur tube ext~rieur bitube ~qu|vatent

R4 : Rayon ext~rieur catandre 6changeur & faisceaux

740 M. Henrion and M. Feidt Vol. 18, No. 5

r4 : Rayon e x t 6 r f e u r tube e x t ~ r i e u r bitube ~:luivatent

tr : Temm de r e t a r d

T(t) : Temp6rature f t u i d e f o r m t i o n du temlps

TO : Temperature f t u i d e i n i t i a t e

To= : Teml~rature f t u i d e f i n a t e

rt~ : Constante de temps

m6tat : C o e f f i c i e n t de conduction paroi tube

[ 1 ] J.WARD MacARTHUR "Comportement de la pompe a chaleur en

regime transitoire : recherche th~orique " , Revue
Internationale du froid , Vol 7 , N"2 , Mars 1984

[ 2 ] M.MENAHEM " Sur l'~tude des modules mathematiques

d'~ohangeurs de chaleur " , R.G.T , N ° I00 , Page 407 ,
Avril 1970

[ 3 ] J.FONTAINE , F.LANNOY " Determination exp~rimentale d'un

module d'~coulement dans un echangeur & plaques " , R.G.T
, N ° 164-165 , Page 573 , Aout-Septembre 1975

[ 4 ] Y.A.GOGUS , O.ATAER , " Effect of fins on transient

behavior of cross-flow air-liquid heat exchangers " ,
I.I.F - I.I.R , commissions BI,B2,EI,E2 , page 219 , 1988
(2) , P U R D U E U . S . A

[ 5 ] D.AZILINON , P.PIERSON , J.PADET " Constante de temps des

echangeurs thermiques " , R.G.T , N ° 338 , page 64 ,
Fevrier 1990

[ 6 ] M.HENRION " Contribution & l'~tude theorique et

experimentale du comportement en r~gime transitoire
d'~changeurs de chaleur " , M~moire Ingenieur C.N.A.M ,
Septembre 1990 , NANCY