THESE

PRESENTEE PAR

CHANDRASA

POUR OBTENIR LE TITRE DE DOCTEUR

DE L'INSTITUT NATIONAL POLYTECHNIQUE DE GRENOBLE

(ARRETE MINISTERIEL DU 30 MARS 1992)

SPECIALITE: GENIE DES PROCEDES

ETUDE EXPERIMENTALE DE L'ENCRASSEMENT

PARTICULAIRE EN PHASE GAZEUSE
SUR LES TUBES DES ECHANGEURS DE CHALEUR

DATE DE SOUTENANCE 20 DECEMBRE 1994

DEVANT LE JURY COMPOSE DE :

PRESIDENT: M. FEIDT
RAPPORTEURS: M. FEIDT
C. JALLUT
EXAMINATEURS: M.RENAUD
J.M. GRILLOT
R. VIDIL

Thèse préparée au
Laboratoire de Thermohydraulique des procédés Industriels
Service de Thermohydraulique pour les applications Industrielles
CEAICENG

dans le cadre des activités du Greth

(Groupement pour la Recherche sur les Echangeurs Thermiques)
 à mon épouse, Rla Rayant'
co-auteur Virtuel de ce travail

à mes familles
à mes amis

1
 AVANT-PROPOS

Cette étude s'mscrlt dans le cadre du programme JOULE menée par le Gre th. La CCE
(CommissIOn des communautés Européennes) dans le cadre du programme JOULE finance
un projet sur l'encrassement des échangeurs de chaleur mené par des laboratoires, dont en
France: le Greth, en Ecosse: le NEL et par d'autres partenaires mdustnels (Français. Belge,
Hollandais et Grec). Ce projet a pour objectif de concevOir des méthodes d'évaluatIOn et de
prédictIOn de l'encrassement des échangeurs de chaleur par des fumées. Il s'appUie sur "des
modèles analytiques et sur des moyens de mesure.

2
 REMERCIEMENTS

Je remercœ A1Jvf. J.M DELHA YE, chef du SerVice de Thermohydraulique pOlir les
ApplicatIOns Indll~trIelles - STI, R. V/DIL, Directeur du Groupement pOlir la Recherche ~71r

les Echangeurs Thermiques - Greth, et P. BANDELIER, Chef du Groupe de

Thermohydraulique pour les ApplicatIOns Industrielles, pour les moyens considérables qu'Ils
011t bien voult, mettre à ma d/~posltlOn.

Je remercie également Mlvf. 1. HUYGHE, ancien chef du STI, et F. LAURO anCien chef du
Laboratoire de Thermohydraulique pour les Procédés Industriels - LTPI, pour m'avoIr
accueilli au sem de leurs groupes re~pectifs.

J'exprime toute ma reconllaIs~ance à M M FEIDT. professeur à l'u1l1verslté de Nancy l, pour

avOIr accepté la présldellce dll jllry de thèse et d'étre un rapporteur de ce travail.
Je tiens aussI à remercier M. C. JALLUT. professeur à l'umver~/té de Claude Bernard à
Lyon, qlll a accepté d'étre UII rapporteur de ce travail lors du JUry de thèse. Je le remercie
tOIi/ parllcullèrenumt pour sa par/lclpatlOn acllve e/ COTlSlnlcllve.

Mes remerciements s'adressent également à A1lvf. M RENA lID, professeur à l'lNPG, et J.M
GRILLOT. maÎtre de cOllférence à l'INPG, qui ont as~7ITé, avec beaucoup d'mtérét, la
directIOn SCientifique de cette thèse.
Je liens à remerCier M R. VIDIL pour sa parllcipatioll au JUry de thèse et pour l'mtérét qU'II
a manifesté pour cette étude.

J'exprime ma gratitude à toutes les personnes côtoyées - Ingémeurs, thésards, staglalres,

secrétaIres, documentalistes et techniCiens - qui m'ont supporté et avec les quels j'aI eu des
relallons les plus amICales.

Enfin, je tœns à remercier Badan Pengkajian & Penerapan Telm%gl (BPPT), Lahoralorlum
TermodmamIl(Q, Motor dan slstem Propulsl (LTMP), et PUSPIPTEK Illdonesla, qUi ont
rendu Vlahle mon séjour en France.

3
 RESUME

Une étude expérimentale de l'encrassement partieulaire sur la paroI des tubes d'un échangeur
thermique a été effectuée en réalisant le dépôt d'encrassement à l'aide de partleules de sulfate
de sodium de tOllle uniforme. Un dispositif expérimental équipé d'un générateur de particules
a été utilisé pour exammer le comportement du dépôt pour différentes condlllons
expérimentales.
Deux études ont été réalisées. La première concerne le dépôt sur la paroI d'un luhe IIl11que.
L'effet de différents paramètres (taille de particules, vitesse et température du fllllde) sur la
cinétique de dépôt a été analysé.
Au cours de la seconde, le comportement du dépôt dans un fOlsceau de tubes à OIlelles a été
étudié. La forme du dépôt observée est /'image des conditions locales de /'écolllemelll aulour
des tubes du faisceau. La résistance d'encrassement est plus Importante pOlir les peilles
partieules et lorsque le débit de gaz est plus faible.

4
 SOMMAIRE

Chapitre 1 INTRODUCTION . ........ ....... . . 13

Chapitre II ETUDE BffiLIOGRAPIflQUE .... ... .... . 15

2 1 GENERALITES SUR L'ENCRASSEMENT . 15

2.1 1 Effet sur la performance thermique ...... .. 16

2 12 Effet sur le dimensionnement d'un échangeur 18
2 13 Modélisation de l'encrassement particulaire. . ... . 19

2.2 MECANISMES DE L'ENCRASSEMENT PARTICULAIRE .... 21

2.2 1 TRANSPORT ALAPAROI ...... ........ . .. .. ... .. 21

2.2.1. 1 Diffusion convective....... ...... . . ....... . 22

2.2 1 2 Forces extérieures (phorétiques) .... .. . 28
2.2 1.3 Diffusion thermophorétique .. 32
2 2.1 4 Inertie de particules.. . . .... ..... ... .. . . 34
2.2.1.5 Impaction inertielle sur un cylindre... .. . 36
2.2 1 6 Efficacité de collecte par impaction. .... . 38
2.2.1 7 Impaction inertielle un faisceau de tubes .. . 43
22 1.8 Interception ..... .. .. .... ..... ..... ....... . 47
2 2 1. 9 Sédimentation. ... ........... ....... . ...... . 48

2.22 ADHESION ET REENTRAINEMENT .. ... . 48

2 2.2 1 Probabilité d'adhésion.... .. ... . . . ...... . .49

2.22.2 Réentraînement.. ............... ............ . 50

2 3 BffiUOGRAPHIE SUR L'ENCRASSEMENT DES TUBES AILETES. 51

2.3 1 GAZ ENCRASSANT ARTIFICIEL ..... ......... ..... . .... 51

2.3.1 1 Etude de El Shobokshy et El Zayat [28]. .... . .. 51

2 3.1.2. Etude de Bott et al [29] ..... ..... .. ........ . 53
2.3 1 3. Etude de Reay et al [30]..... .......... ... ..... . . .. 54

2.3 .2 GAZ ENCRASSANT REEL .. .. . ...... .. . .. ... .... . ..... 56

23 2.1. Etude de Reay et al [30] ... .. .......... . .. .. 56

2.3.22. Etude de Lokshin et al [31], [32] ...... .. . . 57
2.3 23. Etude de Lokshin et al. [33]. .. ... ... . . . .59

5
 2 3 2.4. Etude de Silvestrini [34] ......... ........ . . .. . .. .. 60
2 3 2.5. Etude de Casper et al. [35] .. ..... ............ .. .. 61
2.326. Etude de Roberts et Kubasco [36] ... ....... .. 63

2.3 3 CONCLUSION .. . 64

Chapitre III PARTIE EXPERIMENTALE ............ . .. 65

3 1 PRODUCTION DES PARTICULES ................ .. ............... . 65

3. 1 1 Générateur d'aérosol .......................... . .................. ..... . .. 66

3.1.2 Neutralisation des charges électriques ................ ......... .... . 67

3 2 DISPOSITIF EXPERIMENTAL. ........... ... ...... ........... .... . 67

3 2 1 Instrumentation.. ..... ..... . .... ..... .... .. .. 69

3.2.2 Sections d'essais .. ..................... .. ...... . 69
3 2.2 1 TUBE SONDE.......... .. ........ .. ........ .. . ..69
3 2.2 2 FAISCEAU DE TUBES AILETES .. .. ... .. . 70

33 GRANULOMETRIE DES PARTICULES .................... ........ . .. .. 72

3.3 1 Mesure de la granulométrie.. ............................. .. ...... . ...... . . . .. . 73

3.32 Mise en oeuvre de l'appareil..... ........... .. . ........ . ..... ... .74
3.3 3 Résultats des mesures ................ ....... ..... ... .... . .. .. . 75

Chapitre IV - ETUDE DU DEPOT SUR UN TUBE ISOLE . .. ...... ...... .. . ....... 78

4.1 TRANSFERT CONVECTIF SUR UN TUBE ISOLE..... .... .. .. .. . .. . 78

4.l.1 Effet du Blocage. ................................... ............. ...... . ... .. .. ... ... .79
4.l.2 Effet de la Rugosité ... ............... ..... .... . .......... ..... ... 80

42 ESSAIS DE DEPOT SUR le TUBE SONDE..... ............... ..... ............ .. 81

4.2.1 Description des essais ............................ .......... ....... ..... .. . .... . ......... 81
4.2.2 Résultats expérimentaux........ .......................... . . .. .. ....... .. ........... 82
42.2 1 Essais de la série-A.. ......... . .. ...... .... .. .. ..... . ... .... ........ 82
4.2.2.2 Essais de la série B.... ............. ......... ..... .. ... .. .... ... ..... . ... 85
4 2.2.2 Essais de la série C .... .... ..... . ..... ...... . . ..... .. .. .. .... .. ... ...... 86

4.3 INTERPRETATION DES RESULTATS ........ .. .:...... . .............. ......... .... .....88

43.1 Dépot par impaction inertielle ...... .......... ...... ... ... ... ...... .. '.. ... ., 89
4.3 .2 Dépôt par interception et sédimentation .................... ... . ... . .. ... . .......... 91
4.3 3 Dépot par diffusion et thermophorese..... .. ....... .. ............ .... ............... 93
4 3.4 Effet de l'isolation du dépôt.... ........... ...... ....... .... . . ..... . ..... .. ... .. ...... 95
 4 3 4 1 Augmentation de la vitesse de dépôt. .. .. . .96
4 3.4 2 Température de surface du dépôt....... . 97
4.3 4.3 Résultats de calcul .. ......................... . 98

4 4 CONCLUSION . .... ... ... 100

Chapitre V ETUDE DU DEPOT SUR UN FAISCEAU DE TUBES AILETES 102

5 1 ECOULEMENT DANS UN FAISCEAU DE TUBES . 102

5.1 1 Approche numérique de l'écoulement . .. ........ 105

5 1 2 Distribution de température ....................... 106

5 2 FAISCEAU DE TUBES AILETES EN CONDITIONS ENCRASSANTES 108

5 2 1 Perfonnances thenniques en conditions propres ... 109

5 3 ESSAIS DE DEPOT EN TRANSFERT DE CHALEUR.... 109

5.3 1 Résultats expérimentaux......... .. ............. ......... .... .. 111

5 4 ESSAIS DE DEPOT EN ISOTHERME. ..... ..... .. .... . ... . 121

5.5 ESTIMATION DU DEPOT ............. ... ........... ...... ..... 124

551 Dépôt par diffusion ........ ... ........ ......... .,. ... . 125
5.5.2 Dépôt par impaction inertielle........ .... . .. ....... 127
, A • •

5.5 3 Depot par mterceptlon ......... ........... ... ....... . ... 129
554 Dépôt par sédimentation...... ....................... .... ... 130

5 6 COMPARAISON DES RESULTATS. ........ . ...... ... ...... .. . .. 131

5.7 CONCLUSION.. .. .. ... .. ............... .. ...... . ........... ....... .... . ... . 132

CONCLUSION GENERALE . .... ............................... . ...... . . ... .. . 133

BIBLIOGRAPIiIE... .... .... ..... .. ..... ................. ...... ... .. ... ...... ... ..... ... ...... . . .. . 134

Annexe - A Détermination du coefficient d'échange thermique côté gaz 140

Annexe - B Incertitudes sur le coefficient d'échange de chaleur . ..... .... 143

Annexe - C Incertitude sur la resistance d'encrassement....... ..... .. . .. ...... 146

Annexe - D Récapitulatif des essais realises sur le tube sonde ........ .... . 147

7
Annexe - E Récapitulatif des essais realises sur le faisceau de tubes 148

Annexe - F Proprietes thermophysiques de l'air sec . 156

Annexe - G Proprietes physiques associées avec des aerosols 158

8
 NOMENCLATURE

Al surface intérieure de la paroi

A2 surface extérieure de la paroi
Ad surface d'échange en conditions encrassant es
Am surface moyenne de la paroi
Ao surface d'échange en conditons propres
AproJ surface apparente du collecteur
Cp chaleur massique à la pression constante
CT vitesse de thermophorèse
CO facteur de correction pour le coefficient de traînée
Cp concentration moyenne en particules
Cpe concentration en particules en écoulement
Cps concentration en particules à la paroi
Cs concentration de soluté dans la solution
CCU" facteur de correction Cunningham
d diamètre de la buse dans l'impacteur
da diamètre aéraulique moyen de la particule
de diamètre du tube
dp diamètre géometrique de la particule
Db coefficient de diffusion Brownienne
ed épaisseur de dépôt
f coefficient de trainée de la particule
F811g facteur d'augmentation phorètique
Fu force de trainée de particule
g gravité
gJit force extérieure exercée sur une particule
h l ,h 2 coefficient d'échange du fluide 1 et du fluide 2
H coefficient global d'échange de chaleur
Hd coefficient d'échange en conditions encrassantes
Ho coefficient d'échange en conditions propres
hd coefficient de transfert de masse réel
hm coefficient de transfert de masse pour la captation parfaite
kq facteur du blocage
Kf Coefficient de transfert thermique local adimensionné (fig 4 1)
L longueur caractéristique de l'écoulement

9
m masse de la particule sphérique
md masse de particules déposées par unité de surface
m: masse asymptotique de dépôt
mp masse de la particule
Mg masse moléculaire du gaz
Ml masse de particules injectée
N mt paramètre adimensionnel d'interception
p pression normale à la paroi
Pu puissance calorifique
q flux de chaleur surfacique
qr flux de chaleur produit par unité de volume de fluide
qs flux de chaleur transféré à la paroi
rc flux massique de particules produit par une réaction chimique
R rayon du cylindre (expression 276)
R constante universelle des gaz parfaits
Rd résistance thermique de dépôt
R~ résistance asymptotique d'encrassement
Rp résistance thermique de la paroi
s distance d'arrêt de la particule
S probabilité d'adhésion (expression 2.99)
S pas transversal adimensionné des tubes dans un faisceau
t temps
te constante de temps (expression 2 9)
te durée d'essai corrigée (paragraphe 4.2.2)
tr durée d'essai réelle
T température
u composante de la vitesse le long de la paroi
U_ vitesse moyenne de l'écoulement à l'infini amont
v composante de la vitesse normale à la paroi
vf vitesse locale du fluide
v0 vitesse initiale de la particule
v Pt vitesse relative d'une particule soumise à une force extérieure
V vitesse moyenne de l'écoulement principal
x coordonnée le long de la paroi
y coordonnée normale à la paroi

Nombres Adimensionnels
Kn nombre de Knudsen

10
Le nombre de Lewis
Nu nombre de Nusselt pour le transfert de chaleur
NUe nombre de Nusselt du transfert de chaleur autour d'un cylindre
Pern nombre de Peclet basé sur la vitesse phorètique
Pr nombre de Prandtl
Re nombre de Reynolds
Re p nombre de Reynolds basé sur le diamètre de la particule
Re pl nombre de Reynolds local de la particule
Sc nombre de Schmidt
Sh nombre de Sherwood
Stk nombre de Stokes
Stm nombre de Stanton pour le transfert de masse avec la thermophorèse
Stmo nombre de Stanton pour le transfert de masse sans effet de la thermophorèse

Lettres Grecques

diffusivité thermique du gaz

coefficient de diffusion thermophorètique
diffusivité thennophorétique de la particule
rapport diffusivité thennophorètique des particules et celle du gaz porteur
épaisseur de la couche limite massique
coefficient de diffusion turbulente
paramètre de l'impaction inertielle
dissipation visqueuse
efficacité de collecte
efficacité de collecte pour un mécanisme de dépot
conductivité thennique du gaz
conductivité thermique du dépôt
libre parcours moyen moléculaire du gaz
conductivité thermique de la particule
viscosité dynamique du fluide
viscosité cinématique du fluide
masse volumique
masse volumique du dépôt
masse volumique du gaz
masse volumique des particules
densité de soluté dans la solution
écart type géométrique de la distribution granulométrique

Il
 temps caractéristique de la particule
contrainte de cisaillement à la paroi
densité de flux massique convective
densité de flux massique de particules déposées
densité de flux massique par diffusion convective
densité de flux massique par diffusion
densité de flux massique par des forces phorètiques
densité de flux massique de réentraînement
fraction massique en particules
facteur de correction pour le nombre de Stokes effectif
facteur de correction pour les particules non-Stokesiennes
coefficient de résistance au réentraînement du dépôt

Indices

cr critique
d dépôt
dif diffusion
e écoulement
eff effectif
f fluide
g gaz
pour une classe de taille de particule
imp impaction inertielle
int interception
1 coordonnée adimensionnée
o sans effet de la thermophorèse
p particule
s surface ou paroi
séd sédimentation
sp point de stagnation
T thermophorèse
00 infini amont

Exposant

* paramètres adimensionnels

12
 chap 1 IntroductIOn

Chapitre 1

INTRODUCTION

Les échangeurs thermiques sont présents dans un grand nombre de procédés industriels Ils
permettent d'optimiser les dépenses énergétiques de l'ensemble d'un processus Le transfert de
chaleur se fait généralement par l'intermédiaire d'une surface d'échange En phase gazeuse, sa
géométrie peut être simple, comme dans le cas des faisceaux de tubes lisses utilisés dans les
chaudières, ou plus complexe comme dans le cas des faisceaux de tubes ailetés équipant les
échangeurs compacts utilisés dans les applications aéronautiques

L'encrassement est le résultat de l'accumulation d'un dépôt sur la surface d'échange d'un
échangeur de chaleur. Ce dépôt produit une résistance thermique supplémentaire sur la surface
d'échange. Il peut y avoir différentes conséquences telles que l'augmentation du coût due à
une augmentation de la surface d'échange à installer, la perte d'énergie due à la diminution du
transfert de chaleur et à l'augmentation de la perte de pression, l'augmentation des coûts de
maintenance et de réparation.

L'encrassement sur la surface des échangeurs de chaleur est un problème très étudié depuis
quelques années De nombreuses études théoriques et expérimentales ont été faites par
différents auteurs pour comprendre les mécanismes de dépôt sur une paroi de géométrie simple
ou complexe.
Lors du dimensionnement d'un échangeur, on utilise des valeurs empiriques de la résistance
d'encrassement pour déterminer la surface d'échange nécessaire. Mais, cette approche
empirique peut conduire à un mauvais dimensionnement de l'échangeur, à cause de son
imprécision Une meilleure compréhension des processus physiques de l'encrassement permet
d'obtenir à une meilleure estimation de la valeur de la résistance d'encrassement pour le
dimensionnement de l'échangeur

Les travaux présentés sont consacrés à l'étude des mécanismes du dépôt particulaire en phase
gazeuse sur la surface des tubes d'un échangeur de chaleur L'approche envisagée consiste à
étudier le dépôt sur un tube sonde et le dépôt sur un échangeur

13
 chap 1 IntroductIOn

Une sonde d'encrassement est généralement un tube unique servant à mesurer la masse de
dépôt dans des conditions contrôlées L'utilisation de cet appareil permet de vérifier
précisèment les différents modèles de dépôt
Le tube sonde utilisé est un tube cylindrique hydrauliquement lisse placé perpendiculairement à
un écoulement de gaz chargé en particules Nous nous proposons d'étudier le comportement de
dépôt sur cette géométrie simple, puis d'essayer de comparer le comportement vis à vis de
l'encrassement d'un tube lisse isolé par rapport à un faisceau de tubes à ailettes disposés en
ligne avec comme paramètres : la taille de particules, la vitesse d'écoulement et les conditions
limites thermiques

On propose dans le deuxième chapitre une étude bibliographique sur les différentes
mécanismes de dépôt Le troisième chapitre présente le dispositif expérimental (boucle Gazpar)
utilisé lors de cette étude Le choix et la gamme des paramètres étudiés sont précisés.
Le quatrième et le cinquième chapitres sont respectivement consacrés à l'étude du dépôt sur un
tube isolé (tube sonde) et sur un faisceau de tubes ailetés Les résultats expérimentaux sont
comparés aux prévisions des modèles Le sixième chapitre donne la conclusion générale de
cette étude.

14
 Chap Il Etude Bibliographique

Chapitre Il

ETUDE BIBLIOGRAPHIQUE

Dans ce chapitre, on décrit les résultats des études sur l'encrassement particulaire en phase
gazeuse d'un échangeur de chaleur tubulaire La première partie présente quelques généralités
sur l'encrassement des échangeurs de chaleur; la seconde partie, les mécanismes physiques
gouvernant le transport à la paroi de particules en suspension dans un fluide, ainsI qu'une
discussion sur l'adhésion et le réentraînement Enfin la dernière partie présente quelques
résultats concernant l'encrassement particulaire des échangeurs à tubes ailetés

2.1 GENERALITES SUR L'ENCRASSEMENT

L'encrassement est généralement défini comme l'accumulation d'un dépôt de matériaux ou de

substances non désirées sur la surface d'un échangeur de chaleur Cette accumulation de dépôt
produit une réduction des performances thermohydrauliques de l'échangeur telles que
• une diminution de la valeur du coefficient d'échange de chaleur par l'augmentation de la
résistance thermique de l'interface.
• une augmentation de la perte de pression par la réduction de la section de passage du fluide
et par la formation de rugosités.
Il existe plusieurs types d'encrassement [1] .
1 Encrassement particulaire - accumulation des fines particules qui sont en suspension dans
un fluide sur la surface d'un échangeur de chaleur
2 Entartrage - formation d'un solide cristallin à partir d'une solution liquide sur la surface de
l'échangeur de chaleur.
3. Encrassement par une réaction chimique - formation de dépôt du produit d'une réaction
chimique.
4 Encrassement par la corrosion - accumulation du produit de corrosion sur la surface de
l'échangeur de chaleur. La corrosion est le résultat d'une réaction chimique ou
électrochimique entre la surface et le fluide.
5 Encrassement biologique - adhésion des micro-organismes créant un film sur la surface de
l'échangeur de chaleur.

15
 Chap./I Etude BiblIOgraphique

En général, l'encrassement sur la surface d'un échangeur de chaleur se produit par la

combinaison de ces différents processus.

Ce travail est consacré à l'étude de l'encrassement particulaire en phase gazeuse sur la surface
des tubes d'échangeurs de chaleur Ce genre d'encrassement se produit couramment sur les
échangeurs de chaleur utilisés dans les systèmes de récupération sur fumées industrielles, sur
les échangeurs de climatisation, etc. La plupart des écoulements de flUides gazeux industriels
transportent des particules en suspension dont la taille varie d'une fraction à quelques dizaines
de micromètre Lorsque ces gaz sont véhiculés dans les échangeurs de chaleur, le dépôt de ces
particules cause l'encrassement particulaire.
Les processus principaux gouvernant les mécanismes de l'encrassement particulaire en phase
..
gazeuse sont le transport de particules à la paroi, l'adhésion des particules sur la paroi et le
réentraînement du dépôt Si les mécanismes de transport à la paroi sont bien connus, les deux
autres mécanismes sont encore mal compris

2.1.1 EFFET SUR LA PERFORMANCE THERMIQUE

La présence d'un dépôt sur la surface d'échange d'un échangeur de chaleur se manifeste entre
autres par une augmentation de la résistance thermique de la surface au transfert de chaleur

l
fluide 1 fluide 2

1
,,
, 1
1 \
,
\

Rd' 'Rd
1 Rp 2

Figure 2.1 Résistance thermique de part et d'autre d'une paroi

On peut exprimer le coefficient global surfacique de transfert thermique H rapporté à la surface

extérieure de la paroi (A 2 ) sur la figure 2. 1 par:
1 12 A 1 A A R
2
-=-+--+-R +R d2 + -2 dl (2 1)
H h2 AI hl Am P AI

16
 Chap.II Etude Bibliographique

2.1.2 EFFET SUR LE DIMENSIONNEMENT D'UN ECHANGEUR

Lors du dimensionnement d'un échangeur de chaleur, les valeurs correctes des résistances
d'encrassement doivent être choisies par le projeteur. Celles-ci sont généralement fournies soit
par l'utilisateur de l'échangeur qui est guidé par sa propre expérience, soit par des valeurs
empiriques extraites par exemple des tables TEMA [2] Leur utilisation entraîne évidemment
un sur-dimensionnement de l'échangeur par rapport à celui effectué pour des conditIons de
fonctionnement propre. Pour un bon fonctionnement, il est nécessaire de ne pas sous-estimer
ou sur-estimer les valeurs de la résistance d'encrassement. La sous-estimation conduit à un
faible échange thermique lorsque la surface est encrassée. La sur-estimation conduit à un
surcoût inutil et quelque fois à une modification du régime de fonctionnement pour compenser
la surface d'échange trop importante

Figure 2.3 Excès de surface d'échange en fonction Ho et Rd [2]

En effet, pour un échangeur de puissance Pu, l'équation de base est.

Pu = Ho Ao DTLM (2.2)

Dans les conditions encrassantes, l'expression devient

(23)

où Ao représente la surface d'échange nécessaire pour un échangeur propre, Ad la surface

nécessaire lorsque la résistance d'encrassement est incorporée dans le calcul, DTLM la

18
 Chap Il Etude BiblIOgraphique

différence logarithmique moyenne de température et Hd le coefficient d'échange en conditions

encrassantes :

(24)

En combinant ces deux expressions, on obtient :

Ad
--I=H
A 0
Rd (2 5)
o

La figure 2.3 représente cette équation où on peut remarquer que l'influence de

l'encrassement est d'autant plus importante que le coefficient d'échange thermique global en
conditions propres est élevé.

2.1.3 MODELISATION DE L'ENCRASSEMENT PARTICULAIRE

La plupart des modèles traduisant le comportement asymptotique de l'encrassement -des

échangeurs de chaleur industriels est basé sur l'observation de Kem et Seaton [3] qui
proposent un modèle simple pour prendre en compte l'encrassement asymptotique

pente

t 3t
Figure 2.4 Evolution asymptotique de la masse du dépôt au cours du temps [3]

Pour des conditions bien définies, la densité de flux massique des particules qui s'accumulent à
la paroi est égale à la différence entre la densité de flux massique de particules déposées et la
densité de flux massique de réentraînement, soit :
dm d _
dt"-<Pd -<Pr (2.6)

La densité de flux massique de dépôt est reliée au mécanisme de transport à la paroi par

19
 Chap.JJ Etude BIblIOgraphIque

(2 7)

hm représente le coefficient de transfert de masse, Cpe et C ps sont respectivement la

concentration en particules dans l'écoulement principal et à la paroi
Cependant, le mécanisme de réentraÎnement dépend de l'équilibre entre les forces de
cisaillement à la paroi et les forces d'adhésion des particules à la paroi, soit
rnd
=C't _
ln (2 8)
""r s 'If d

où 't. est la contrainte de cisaillement à la paroi, 'If d le coefficient de résistance au

réentraÎnement du dépôt, m d la masse de particules déposées par unité de surface, et C une
constante.

En intégrant (2 6) avec la condition initiale t = 0, md = 0 et en supposant que le flux massique

de particules déposées reste constant, on obtient .

md = m; [1 - exp (- t: )] (2 9)

te est la constante de temps, qui est égale à

t e = rn d (2 10)
<Pr
et rn •d est la masse asymptotique de dépôt donnée par

md*=<Pd·te (2 Il)

La vitesse initiale de dépôt est proportionnelle à

dm, = m~ (2 12)
dt t=o t
c

L'évolution de la masse du dépôt m d peut être également exprimée en fonction de la résistance

d'encrassement. Si on suppose la masse volumique et la conductibilité thermique du dépôt
constantes, on peut alors exprimer la résistance thermique du dépôt par

(2 13)

et l'évolution de la résistance d'encrassement devient .

Rd = R; [ 1- exp (- t:)] (2.14)

où R~ la résistance d'encrassement asymptotique donnée par·

(2 15)

20
 Chap.JJ Etude BiblIOgraphique

De nombreux progrès dans l'étude des mécanismes de dépôt caractéristiques de l'encrassement

particulaire ont été réalisés. Plusieurs modèles sont disponibles dans la littérature Cependant, il
est encore difficile d'estimer la valeur exacte de la résistance d'encrassement pour une condition
donnée car de nombreux paramètres peuvent intervenir simultanément Ces paramètres sont les
conditions de l'écoulement, le type de matériau de la surface d'échange, les propriétés
physiques du dépôt, etc. En général, les effets de ces paramètres sont mal connus

2.2 MECANISMES DE L'ENCRASSEMENT PARTICULAIRE

L'encrassement particulaire de la surface d'échange d'un échangeur de chaleur est un

phénomène complexe. D'une manière générale, l'apparition et le développement de
l'encrassement particulaire peut se décomposer en cinq phases.
(1) l'initiation,
(2) le mécanisme physique de transport des particules à la paroi,
(3) l'adhésion des particules sur la paroi,
(4) le réentraînement des particules déjà déposées,
(5) le vieillissement du dépôt.
Ces mécanismes sont fonction de différents paramètres tels que la température du gaz et celle
de la paroi, la vitesse du gaz, l'intensité de la turbulence, l'état de la surface. les forces
s'excerçant sur les particules, ainsi que les caractéristiques des particules (forme. composition,
concentration, taille). Parmi les cinq phases, les phases (2), (3), et (4) seront développees dans
les paragraphes suivants.

2.2.1 TRANSPORT A LA PAROI

Parmi ces cinq mécanismes conduisant au dépôt, le transport à la paroi est le mécamsme le plus
étudié et donc le plus connu Des modèles expérimentaux ou théoriques sont disponibles dans
la littérature Le transport à la paroi peut se produire par transfert de masse convectif, par
transfert sous l'action des forces extérieures telles que la gravité, la thermophorèse,
l'électrophorèse, la diffusiophorèse, et par effet de l'inertie des particules Ces mécanismes
peuvent coexister dans les échangeurs de chaleur industriels

Pour les particules de taille inférieure ou égale au micromètre, le transfert est généralement
contrôlé par la diffusion [4] La thermophorèse s'ajoute lorsqu'il existe un gradient de
température entre le gaz et la paroi La figure 2.5 représente l'efficacité de collecte d'un tube
avec ou sans échange de chaleur pour différentes tailles de particules Sur cette figure, on peut
remarquer que l'effet de la thermophorèse est important pour l'encrassement des échangeurs
thermiques

21
 C'hap.lI Etude BiblIOgraphique

Figure 2 5 Efficacité de collecte des particules d'après Rosner et al [4]

Pour les particules de tailles supérieures à la dizaine de micromètres, c'est l'inertie qui est le
paramètre contrôlant le dépôt Cependant, pour les classes de tailles intermédiaires, les
différents phénomènes coexistent. De plus, lorsque l'écoulement est turbulent la fluctuation
turbulente de vitesse peut projeter les particules vers la paroi
Dans les paragraphes suivants, l'accent est mis sur les mécanismes de la diffusion convective,
de la thermophorèse, de l'impaction inertielle, de l'interception et de la sédimentation

2.2.1.1 DIFFUSION CONVECTIVE

Le transfert de masse par diffusion dans un écoulement de fluide est généralement appelé la
diffusion convective Ce mécanisme est fondamental dans les systèmes de nettoyage des gaz
comme la filtration, par exemple. Ce mécanisme est généralement important pour des
particules de taille inférieure à 1 Jlm. L'intensité du mouvement de diffusion augmente lorsque
la taille des particules décroît.
On décrit dans ce paragraphe le principe fondamental du processus de diffusion convective
ainsi que les équations décrivant le transfert en présence de ce mécanisme On présente
également le concept d'analogie entre les transferts de chaleur et de masse

Généralités.

Lorsqu'un fluide chargé de particules s'écoule autour d'un obstacle avec une différence de
concentration entre l'écoulement principal et la paroi, il peut se produire un transfert de

22
 Chap.l/ Elude BIbliographIque

particules par convection La densité de flux massique convectif à la paroi (expression 2 7) est
donnée par
(2 16)

Le coefficient de transfert de masse hm dépend des propriétés physiques du fluide et de la

particule, de la géométrie de surface et de la configuration de l'écoulement autour de l'obstacle
Une couche limite de concentration se développe entre l'écoulement principal et la paroi au
sein de laquelle le transfert de masse se produit par une combinaison entre la convection et la
diffusion Cependant, près de la paroi le transfert est contrôlé par la diffusion et la densité de
flux massique des particules à la paroi peut être exprimé par .

=-Db ac
p 1
<Pd ay y=o (2 17)

où Db représente le coefficient de diffusion Brownienne et y la coordonnée normale à la paroi

En utilisant les deux expressions ci-dessus, on peut exprimer le coefficient de transfert de
masse sous la forme .

(2 18)

On peut remarquer que le coefficient de transfert de masse convectif dépend du gradient de

concentration à la paroi et ce dernier dépend des conditions aux limites appliquées

D'un autre côté, le coefficient de transfert de masse, hm peut également être exprimé sous la

10·

s,o. NRTICLU
150011 ,"It

10"' 10" 10° 10' 101

d. ("ml

Figure 26 Nombre de Schmidt en fonction de diamètre de particules de Si0 2 [S]

23
 Chap II Etude BIblIOgraphIque

forme adimensionnelle par intermédiaire du nombre de Sherwood Sh défini par

Sh =hm L (2.19)
Db
qui est généralement fonction du nombre de Reynolds du gaz, du nombre de Schmidt et de la
géométrie de la surface (L est une longueur caractéristique). Le nombre de Schmidt est donné
par :

(220)

où v la viscosité cinématique. Il joue pour le transfert de masse le même rôle que le nombre de
Prandtl pour le transfert de chaleur.
Dans le cas d'un aérosol de fines particules, le nombre de Schmidt est fortement fonction de la
taille de particules. La figure 2.6 montre la variation du nombre de Schmidt pour la diffusion
de particules de Si0 2 dans l'air [5].

Equations décrivant la diffusion convective.

Dans ce paragraphe, on décrit les mécanismes physiques déterminant le comportement de la

couche limite et leur influence sur le transport convectif en développant les équations décrivant
les couches limites hydraulique, thermique et massique
Considèrons les couches limite hydraulique, thermique et massique qui se développent
simultanément à la paroi comme le montre la figure 27 ci-dessous (cité dans [6]) .

CA... "_

cS
1

Figure 2.7 Développement de la couche limite hydraulique, thermique et de concentration

pour une paroi arbitraire [6]

24
 Chap.ll Etude Bibliographique

Le fluide peut être considéré comme un mélange de deux substance A et B (les fines particules
dans un gaz, par exemple). Sur cette figure, le choix des épaisseurs relatives des couches
limites (Ô I > ôe > ô) est arbitraire Ainsi, si on considére l'écoulement établi bidimensionnel
d'un fluide incompressible de vitesse amont uniforme à l'infini, les équations du mouvement
dans la couche limite pour un volume de contrôle fixe s'écrivent (figure 2 7)
dU
-+-=0
av (221)
dx dy

(222)

(223)

où u et v représentent les composantes de la vitesse respectivement suivant les directions x et

y, v la viscosité cinématique du fluide et p la pression statique Le bilan d'énergie se traduit
quant à lui par l'expression suivante.

(224)

avec 1.. la conductivité thermique du fluide, qr le flux de chaleur produit par unité de volume de
fluide, et cJ> la dissipation visqueuse donnée par.

{( dy, + dx )2 + 2[(dX)2+ (dy )2] -"3 (dX + dy )2 }

dU dV dU dV 2 dU dV
<1> Eql
(225)

Le bilan de particules supposées soumises au mouvement de diffusion Brownienne s'exprime

comme suit:

(226)

où rc représente le flux massique de particules produit par unité de volume par une réaction
chimique éventuelle

Les équations développées ci-dessus décrivent les mécanismes physiques conditionnant les
couches limites hydraulique, thermique et massique pour un écoulement établi, bidimensionnel
Les difficultés de résolution de ces équations viennent du fait qu'on ne peut considérer les
paramètres caractérisant le gaz (1.., p et ~) comme des constantes Cependant, un grand nombre
de modélisations s'en tiennent aux cas le plus simple du fluide parfait incompressible Les
phénomènes de transfert de chaleur et de masse sont alors décrits à l'aide de valeurs moyennes.

25
 Chap Il Etude BiblIOgraphique

Dans certains cas ,simples, on peut considérer que le fluide a des propriétés physiques
constantes. En l'absence de réaction chimique et de source de chaleur et en considérant
l'épaisseur de la couche limite très faible, on peut alors appliquer les hypothèses suivantes

1) la composante de vitesse le long de la paroi est beaucoup plus importante que celle normale
à paroi (u»v), et les gradients de vitesse normaux à la paroi sont également plus importants
que ceux le long de la parOi «()y
.du dU av
» dx' ()y» dx) ·
dV

2) les flux de chaleur par conduction et de matière par diffusion dans la direction normale à la
paroi sont plus importants .
dT
-»-
dT dC ac p
- -p » - -
()y dX' ()y dx
3) On peut également supposer que le transfert de masse à la paroi a un effect négligeable sur
la couche limite hydraulique (la vitesse des molécules du fluide reste nulle à la paroi)
Ainsi, comme dans la plupart des cas la concentration moyenne en particules est très f~ible
(Cp == 0), on peut également considérer que les propriétés thermophysiques du mélange gaz-
particules sont égales à celles du fluide porteur
En tenant compte de ces hypothèses, les équations de conservation dans la couche limite
deviennent :
dU av
-+-=0 (227)
dX ()y

du du 1 dp d2 U
u-+v-=---+v- (2.28)
ax dy P ax ()y 2

(229)

dC p dC p d 2C p
u--+v--=D b - - (230)
dX dy ()y 2

avec ex la diffusivité thermique du fluide (ex = ~). Dans la plupart des cas, le terme ~(dU)2
pC p cp dy
dans l'équation (2.29) est négligeable par rapport aux autres termes
Quand on examine les équations 2 28, 2.29 et 2.30, on peut noter qu'il y a une similitude entre
elles. Si on suppose que les termes du gradient de pression dans l'équation 2 28 et de la
dissipation visqueuse dans l'équation 2.29 peuvent être -négligés, on peut remarquer que les
trois équations sont de mêmes formes. On peut noter que ces équations décrivent le
comportement de la couche limite en convection forcée à faible vitesse que l'on peut rencontrer
dans diverses applications industrielles.

26
 Chap Il Etude BIblIOgraphIque

Analogie entre les transferts de chaleur et de masse.

D'une autre façon, on peut également exprimer les équations décrivant les couches limites sous
une forme adimensionnée en utilisant les variables suivantes .
x· : ~ et y. : 1. où L est une longueur caractéristique de la paroi
L L
• u • v
U : - et v : - où U _ est la vitesse moyenne débitante à l'infini amont
U_ U_

T. = T-T. (231)
Te -T.

(232)

•- p (233)
p = (pU_ 2 )

Après substitution de ces variables, les équations de conservation dans la couche limite
deviennent

(234)

• au· • au· ap· 1 a 2u·

u -+v
axo
-
ay.= -axe
- +Re
-- -
ay.2
(235)

aveclesconditionsauxlimites: u·(x·,O)=O ,v·(x·,O)=O et u·(x·,oo) = 1

o o
o aT • aT· 1 a 2T
u -+v --=-- (236)
axo ay· Re Pr ayo2

avec les conditions aux limites' T· (x· ,0) = °, T· (x· , 00) = 1

• ac· • ac· 1 a 2c·
u --+v --=-- (237)
axe ay· Re.Sc ayo2

avec les conditions aux limites' C· (x· , 0) = ° , C· (x0,00) = 1

où
Re: U_L ,Pr:~ Sc:~
va' Db
On peut remarquer sur les équations 236 et 237 que pour des conditions aux limites
semblables les transferts convectifs de chaleur et de masse sont analogues Cette analogie n'est
valable que si les hyphothèses d'analogie thermique et massique sont respectées [7].
Les hypothèses d'analogie thermiques sont :

27
 Chap.Jl Etude BiblIOgraphique

- la composition du fluide est unifonne dans l'espace (pas de réaction chimique, ni de transfert
de masse).
- les conditions aux limites thenniques sont simples (par exemple température de paroi Ts
constante).
- les propriétés thennophysiques du fluide sont supposées constantes
- toutes les sources de dissipation d'énergie, incluant dissipation visqueuse et réaction
chimique, sont négligeables.
Par ailleurs, les hypothèses d'analogie massique sont
- les particules sont très diluées (fraction massique de particule très faible)
- la concentration en particules le long de la paroi est supposée constante.
- les forces extérieures (phorètiques) sont négligeables.
- Il n'y a pas de réaction chimique homogène
Si ces conditions d'analogie sont respectées, les relations décrivant le transfert de chaleur et de
masse sont analogues. Dans ce cas, le coefficients de transfert de masse adimensionnés
(nombre de Sherwood Sh) et le coefficient de transfert de chaleur adimensionné (nombre de
Nusselt NUb) ont la même fonction

2.2.1.2 FORCES EXTERIEURES (PHORETIQUES)

L'un des facteurs conduisant à la rupture de l'analogie entre les transferts de chaleur et de
masse est la présence des forces extérieures (phorétiques) dans le système d'aérosols telles que
la thennophorèse, l'électrophorèse, la diffusiophorèse.
On présente dans ce paragraphe le rôle de ces forces dans le processus de transfert de fines
particules par le mécanisme de diffusion convective Nous nous rappelons qu'en général le
transport par convection peut être décomposé en deux parties : la convection due au
mouvement du fluide transportant les particules, et la convection due aux forces extérieures
s'excerçant sur les particules [7].

Considérons maintenant le cas d'un écoulement de gaz chargé de particules de fraction

massique CJ)p (ou de concentration massique Cp). La densité de flux massique locale convectée
est donnée par .
(2.38)

où p et V sont respectivement la masse volumique et la ~tesse du gaz porteur

Si les particules sont entraînées par l'une des forces d'extérieures avec une vitesse relative
quasi-pennanente v lit ' le flux convectif associé peut être exprimé sous la fonne :
(239)

28
 Chap./J Etude BiblIOgraphique

Ainsi, la présence du gradient de concentration dans le fluide entraîne une diffusion vers la
zone de faible concentration et le flux massique des particules par cette diffusion est donné par
(240)

Ce terme est généralement appelé le flux de diffusion de Fick pour la diffusion des molécules
du fluide Dans le cas des aérosols, ce terme est appelé le flux de diffusion Brownienne Si
l'écoulement est turbulent, c'est le flux de diffusion turbulente
On peut maintenant exprimer le flux total par diffusion convective par la relation
(241)

Dans le cas des aérosols, si on considére que sur chaque particule s'exerçe une force extérieure
gJ:Ù' on peut alors exprimer la vitesse résultante de la particule par

_ mp gph
v ph - f (242)

où f est le coefficient de trainée de la particule (voir paragraphe 2 2 1 4)

Electrophorèse.

Si gJ:Ù est due à la présence d'un champ de forces électrostatiques dans le milieu, la migration
des particules est appelée électrophorèse [7]. Ces forces proviennent soit de la charge
électrique atmosphèrique, soit du contact entre les particules L'effet de la charge
électrostatique est généralement important pour les particules de petite taille Dans le cas des
aérosols, l'effet électrostatique devient plus important lorsque la taille des particules est
inférieure à 0.1 J.1m [8].

DifTusiophorèse.

La diffusiophorèse permet d'enlever les particules du fluide porteur lorsqu'un deuxième

composant diffuse à travers celui-ci sous l'action d'un gradient de concentration, les particules
microniques suivront le composant qui diffuse [8]. La diffusiophorèse contribuera au dépôt des
particules dans les récupérateurs de chaleur dont la paroi est à une température inférieure à la
température de rosée de l'un des constituants du fluide gazeux porteur.
La diffusiophorèse existe dans le cas de condensation sur la paroi. Les molécules de vapeur
venant se condenser sur la paroi entraînent avec elles des particules solides Lorsque la vapeur
se condense sur la surface à une température inférieure à la température de rosée, la vapeur
crée une force sur les particules qui sont en suspension dans le gaz et entraîne celles-ci vers la
surface [2].

29
 Chap Il Etude BiblIOgraphique

Thermophorèse.

Lorsque le mouvement des particules est dû à un gradient de température, le mécanisme est

appelé la thermophorèse Les particules se déplacent de la région de haute température vers la
région de faible température La thermophorèse est importante dans la plupart des systèmes
d'aérosols et dans les systèmes de récupération sur les fumées industrielles Dans ces cas, les
autres forces phorétiques sont généralement négligeables par rapport à la force de
thermophorèse En général, ce mécanisme contrôle le transfert de particules de taille de l'ordre
o 1 J..lm, et il reste important jusqu'à 5 J.lm [7]. La vitesse de thermophorèse est proportionnelle
au gradient de température, et elle peut être exprimée sous la forme .

(243)

où nT représente le coefficient de diffusion thermophorètique, et Db est la diffusivité

Brownienne Les valeurs du produit nTDb dépendent du mouvement des particules dans le
fluide.

Dans le régime moléculaire où la particule est de dimension négligeable par rapport au libre
parcours moyen des molécules du gaz, la distribution des vitesses des molécules gazeuses n'est
pas modifiée par la présence de la particule Dans ce cas, en supposant que les chocs
moléculaires sur la particule sont élastiques, le bilan de quantité de mouvement cédée à la
particule conduit à l'expression suivante (Waldman, cité dans [9])
(244)

avec v f représente la viscosité cinématique du gaz.

Pour les particules de dimension non négligeable par rapport au libre parcours moyen des
molécules du gaz (régime continu et intermédiaire), la présence de la particule modifie dans
son voisinage le mouvement des molécules gazeuses et la diffusion de la chaleur A partir des
équations de Navier-,Stokes et de Fourier, Epstein (cité dans [9] ) a établi

IlT Db = 1.5 v, C1..~:1...) (245)

où À f et Àp représentent respectivement la conductivité thermique du gaz et de la particule.

Selon cet auteur, cette relation n'est valable que pour des particules de faibles conductivités
À
thermiques telles que À; < 10 .
Beaucoup d'auteurs l'améliorent en précisant les variations de température au voisinage de la
particule Derjaguin et al. [10] donnent un modèle vérifié expérimentalement:

30
 Chap.lI Etude BIblIOgraphIque

Âf C Kn
~+ t·
p
(246)
Â
1+2f+2C t Kn
p

avec Kn le nombre de Knudsen et Ct = 2.17 (coefficient d'ajustement empirique)

Par ailleurs, pour les particules de Si0 2 en suspension dans l'air Tabot et al. (cité dans Rosner
et Park [II)) proposent l'expression suivante '

ex.TLe _
2Cs(~f + ctKn)[1+ Kn(A + A e-
p
l
,
2
A3 Kn
/ )]

(247)
p;:-- (1+3CmKn{I+2~: +2C,Kn)
avec At=12 ,A 2 =0.41 ,A 3 =088 ,et cs =117 Ct =2.18 cm =114 et Le est le
nombre de Lewis donné par '
Le: Db (248)
ex.

où Cl la diffusivité thermique du gaz (Cl = ~ )

pC p
Rappelons qu'en physique des aérosols, on définit trois domaines de dimensions caractérisés
par la valeur du nombre de Knudsen Kn Le nombre de Knudsen est le rapport du libre
parcours moyen des molécules de gaz Ig au rayon de la particule rp
1
Kn : ...!. (2 49)
rp

Le libre parcours moyen moléculaire du gaz donné par [12] .

Il ~"M. (2.50)
Ig = -O.-4"-99-p 8 R T

où Mg la masse moléculaire du gaz considéré, R la constante universelle des gaz parfaits et T

la température absolue. Pour l'air, on peut utiliser une relation simplifiée donnée par.

Ig = 2.26 10-5 T (251)

P
où p est la pression absolue en (pa)
Lorsque Kn »1, on se trouve dans le domaine moléculaire et la particule peut être considérée
comme une molécule. Le mouvement des aérosols est alors décrit en utilisant les lois de la
théorie cinétique des gaz Pour Kn -t 0 (jusqu'à Kn < 0.2), le milieu qui entoure la particule
peut être considéré comme continu et on utilisera alors les lois de la mécanique des fluides. La
région intermédiaire, soit 0.2 < Kn < 10, correspond au domaine de transition

31
 Chap.lI Etude BIblIOgraphIque

2.2.1.3 DIFFUSION THERMOPDORETIQUE

L'un des facteurs provoquant la rupture de l'analogie entre les transferts de chaleur et de masse
est la présence des forces phorètiques dans le système des aérosols Les forces phorètiques
(thermophorèse, électrophorèse, diffusiophorèse) vont affecter le profil de concentration et
augmenter le flux massique de particules à la paroi. Rosner [7] a mis au point une approche qui
prend en compte le fait que la thermophorèse n'est pas négligeable En utilisant le nombre de
Stanton pour le transfert de masse St m, on peut exprimer la densité de flux massique à la paroi
par diffusion (expression 2 40) par la relation:

= - P Db aro
p
q>dtf dy 1y=o -= P UmStm(rope -rops) (252)

où U m la vitesse moyenne du gaz et (rope -rops)le potentiel de transfert de masse

Lorsque l'on prend en compte la thermophorèse, le flux massique à la paroi par diffusion
devient ·

Dans ces conditions, la présence de thermophorèse affecte le profil de concentration ce qui

modifie la diffusion à la paroi Rosner [7] a qualifié la thermophorèse comme un facteur
augmentant la diffusion à la paroi et recherche une expression de ce facteur donnant
l'accroissement de la diffusion dû à la thermophorèse .
(254)

avec Fsug le facteur d'augmentation thermophorètique

L'auteur a estimé que le facteur Fsug est fonction du nombre de Peclet basé sur la vitesse
phorètique v Pl' l'épaisseur de la couche limite massique 0m,o et le coefficient de diffusion
Brownienne et qu'il s'exprime sous la forme:

F = Pem (255)
aug -l-exp(-Pe m)

où Pem le nombre de Peclet de la particule donnée par .

Pe m ~ v~:."'. ]
[ (256)

En utilisant l'expression de l'épaisseur de la couche limite massique 0m,odonnée par.

L
o =-
Sh
01,0
(257)
o

on peut alors exprimer Pem sous la forme .

32
 Chap.II Etude BiblIOgraphique

(258)

avec v th est la vitesse de thermophorèse (expression 2 43) donnée par:

V th = - aTD b grad (ln T) (259)

Si on utilise la loi de Fourier (q =-Â. grad T), on peut ecrire

v th = aTD b Igrad
T
TI = aTD b q
Â.T (260)
s

Ainsi, en utilisant l'expression du coefficient adimensionnel de transfert de chaleur surfacique

(nombre de Nusselt) donnée par

NUh = (T q-T) (261)

Â. - s
L
on peut exprimer la vitesse de thermophorèse sous la forme·

v th =a TD b NUh T_ -Ts (262)

L Ts

où q est la densité de flux de chaleur transmis sur une surface considérée et (T_ - Ts ) la
différence de température entre l'écoulement principal et la paroi
Le nombre de Peclet de la particule (expression 2 56) peut être alors exprimé par

Pe =a NU h T_ -Ts (263)
m TSh
,os
T

D'un autre côté, en ~ppliquant l'analogie entre les transferts de chaleur et de masse dans le cas
du transfert convectif autour d'un cylindre, on peut alors établir la relation suivante
Sc
P em =a T Le( -
)2/3 T_ - T s
(264)
Pr Ts

Le terme atLe représente le rapport entre la diffusivité thermophorètique des particules et la

diffusivité thermique du gaz porteur. Pour les particules de Si0 2 en suspension dans un
écoulement d'air à la température 1500 K, Rosner et Femandez de la Mora [13] ont proposé
pour la variation de a tLe en fonction du diamètre des particules les courbes de la figure 2 8
Rosner et al. [4] préconisent également que les valeurs de atLe données sur la figure 28
peuvent être utilisées pour des particules de Na 2S0 4 et de C dans l'air
Enfin, Gokoglu et Rosner (cité dans [7]) ont amélioré l'expression du facteur d'augmentation
par thermophorèse et montré que dans le cas où le nombre de Schmidt est très élevé la valeur
du Faug doit être réduite par un facteur de correction donné par :
T -T
Fe = exp( -aTLe - S) (2.65)
Ts

33
 Chap Il Etude BIblIOgraphIque

100r---~--~----~--'---~--~
080
060 p: 0 lolm
040
GTLe
020

010
008
006
004 SI02 PARTICLES
I500K AIR

002

o OI
OL.4------l
o~.3----l.10-:".2--1.L0•.,....1---L::-----lIO~1- - . . J
102
(~m)

Figure 2 8 Les valeur du rapport entre la diffusivité thermophorètique et

la diffusivité thermique de particules Si0 2 dans l'air [13]

En conclusion, nous avons évoqué dans ce paragraphe les aspects du transfert de masse des
fines particules par, diffusion-convective Nous avons également analysé les effets de la
présence des forces phorètiques sur le phénomène de diffusion Celles-ci peuvent augmenter la
vitesse de transport de fines particules
Pour le cas du transfert de masse des substances nonréactives sans présence de forces
phorètiques, on peut directement appliquer l'analogie entre les transferts de masse et de
chaleur Dans le cas de la présence des forces phorètiques, on peut cependant établir des
facteurs de corrections afin de tenir compte l'effet de ces forces sur le transfert de masse

2.2.1.4 INERTIE DE PARTICULES

Ce paragraphe est consacré au mécanisme de dépôt par impaction inertielle des particules sur
la paroi d'un obstacle immergé dans un écoulement de fluide L'impaction inertielle est
généralement importante pour les particules de taille supérieure à lllm. Ce phénomène
contrôle les mécanismes de dépôt dans les systèmes de filtration, de séparation [14]
Lorsqu'un écoulement rencontre un obstacle, les particules, du fait leur inertie, peuvent quitter
les lignes de courant et elles pourront heurter l'obstacle et s'y fixer Dans la plupart des cas,
c'est l'accélération de l'écoulement qui conduit à l'impaction inertielle des particules.

34
 Chap.JJ Etude BiblIOgraphique

En général, l'interaction entre les particules et l'obstacle est caractérisée par un paramètre sans
dimension appelé le nombre de Stokes, qui est égal au rapport entre les forces d'inertie et les
forces de traînée auxquelles est soumise la particule L'impaction inertielle ne se produit que
lorsque le nombre de Stokes de la particule dépasse une valeur critique qui dépend de la forme
de l'obstacle et de la configuration de l'écoulement autour de ce dernier

Pour une particule sphérique en mouvement avec une vitesse constante dans un fluide
immobile ou dans un écoulement permanent à faible vitesse (Re<< 1), on peut exprimer la force
de frottement à laquelle est soumise la particule (force de traînée) par la loi de Stokes [14]
FT = 37t Jl d p Vp (2 66)
où V p est la vitesse de la particule par rapport au fluide, d p le diamètre de la particule
sphérique Pour des particules soumises à une accélération dans un écoulement de faible
vitesse, Friedlander [14] a proposé pour estimer la force de traînée d'utiliser l'expression de
Millikan-Cunningham
37t Jl d p
Ftr = Vp (267)
C cun

où C cun est le facteur de correction de Cunningham donné par

CClIII =[I+Kn(1.2+0.41 e -0 88 / Kn )] (268)

Si on considère maintenant une particule sphérique de masse m, projetée dans un écoulement

stationnaire avec une vitesse initiale Vo, l'équation de quantité de mouvement de la panicule
dans les conditions d'application de la loi de Stokes, s'écrit:
dVp
m - = -37t Jl d p Vp (269)
dt
En supposant la particule sphérique, sa masse est donnée par ' m = 7tPpd~ /6
on obtient par intégration :
t

Vp = V0 e 't, (270)

où 't p représente le temps caractéristique de la particule donnée par :

_ pp d p 2
't --- (271)
p 18Jl

Lorsque t -+ 00 , la particule parcourt une distance finie s (distance d'arrêt de la particule)

donnée par :
2
s = pp d p Vo (2.72)
18Jl

35
 Chap Il Etude Bibliographique

On peut ainsi exprimer le nombre de Stokes par le rapport entre la distance d'arrêt de la
particule et une longueur caractéristique de l'écoulement L (dimension de l'obstacle ou
épaisseur de la couche visqueuse dans le cas de l'écoulement dans une conduite)

ppd p 2 Vo
Stk = (273)
18JlL

2.2.1.5 IMPACTION INERTIELLE SUR UN CYLINDRE

L'écoulement autour d'un cylindre est un modèle de base pour les systèmes de filtration et
pour l'étude du dépôt sur des tubes et sur d'autres obstacles immergés dans l'écoulement d'un
gaz. Nous avons vu au paragraphe précédent que l'impaction inertielle sur un cylindre dépend
fortement de la configuration de l'écoulement autour de ce dernier est fonction de nombreux
paramètres tels que . la vitesse de l'écoulement incident, la viscosité et la densité du fluide, la
taille et la masse volumique des particules, la dimension et la géométrie de l'obstacle. Ces
paramètres sont généralement exprimés sous formes adimensionnelles par le nombre de Stokes
et le nombre de Reynolds de la particule et l'impaction inertielle sur un obstacle ne peut se
produire que lorsque le nombre de Stokes dépasse une valeur critique.

-U.. --

lignes d'écoulement

trajectoire des partlcules

Figure 2.9 Impaction inertielle de particules sphériques sur un cylindre [14]

Considérons maintenant le cas de l'écoulement d'un gaz chargé de particules au voisinage d'un
obstacle. A proximité de la paroi l'écoulement est accé~éré, et si l'inertie des particules est
importante la vitesse des particules sera différente que celle du fluide. Dans ces conditions,
l'équation de mouvement de la particule s'écrit:
dV
mT!p =-f(V p - Vf ) (2.74)

36
 Chap.ll Etude BIblIOgraphIque

où Vf est la vitesse locale du fluide et f le coefficient de traînée. Dans le cas où la vitesse de

l'écoulement est faible, on peut supposer que les particules suivent la loi de Stokes et le
coefficient de traînée est donné par
(275)

Pour un écoulement potentiel, on peut considérer que la composante de la vitesse suivant la

direction x est donnée par [14] :
(276)

où b est une constante qui dépend de la forme de l'obstacle et Xl ==~, on peut alors exprimer
R
l'équation de mouvement de la particule par :
m du p
---+u -u=O (277)
f dt p

ou

(278)

on peut encore écrire :

d2 x dx
t p R -2_l + R_l + bU_(x l + 1) = 0 (279)
dt dt

Ainsi, en posant . Xll = Xl + 1 et 9 = t ~- , l'équation de quantité du mouvement des

particules Stokesiennes devient:

d 2x I_I +_1_1
Stk-- dx + bx = 0
(280)
d9 2 d9 Il

La solution de cette équation est donnée par :

X -
Il -
A 1eÂ.,8 + A 2 eÂ.2 8 (281)

où AI et A 2 sont les constantes d'intégration, et )..1 et)..2 représentent les racines de l'équation
caractéristique :

(2.82)
Les valeurs de )..1 et)..2 sont données par :

)..12 = __I_[I±(1-4bStk)05] (2.83)

. 2Stk

La valeur de Stk = _1_ représente la valeur limite au dessous de laquelle l'impaction s'annule.
4b
D'autre part, la valeur du coefficient b de l'équation (2.76) dépend de la géométrie de
l'obstacle. D'après Friedlander [14], pour l'écoulement potentiel autour d'un cylindre le

37
 Chap.ll Etude B,bltograph,que

coefficient b est égal à 2, et le nombre de Stokes critique, Stk cr = 0 125 (Pour un écoulement
potentiel autour d'une sphère, le coefficient b est égal à 3).
Cependant, dans la plupart des cas, l'écoulement est visqueux et dans la région proche du point
de stagnation la composante de vitesse ur est proportionnelle à (XI + 1)2 Dans cette région,
les particules ne se comportent pas comme des particules Stokesiennes Néanmoins, les
approches développées ci-dessus peuvent donner une bonne estimation lorsque la taille des
particules est supérieure ou égale à l'épaisseur de la couche limite qui se développe sur la paroi
de l'obstacle [14]

2.2.1.6 EFFICACITE DE COLLECTE PAR IMPACTION

De nombreux auteurs ont essayé de déterminer les trajectoires des particules dans un
écoulement autour d'obstacles de différentes formes. A partir de ces trajectoires on calcule le
flux de dépôt à la paroi par impaction inertielle. En général, ce flux massique de dépôt- est
exprimé en fonction d'une efficacité de collecte définie comme étant la fraction des particules
ayant traversée la surface apparente amont de l'obstacle et qui sont réellement capturées
De nombreuse études théoriques et expérimentales concernant l'impaction inertielle des
particules sont disponibles dans la littérature Le cas le plus étudié est celui de l'écoulement
perpendiculaire à un cylindre isolé car il sert de modèle de base dans de nombreuses
applications industrielles (filtration, encrassement particulaire des tubes d'échangeurs de
chaleur, etc).
Brun et al. [15] ont fait le calcul numérique des efficacités de collecte pour l'impaction
inertielle des particules sur un cylindre dans l'hypothèse d'un écoulement potentiel En
supposant que toutes les particules heurtant la paroi y restent, on représente sur la figure 2 10
leur résultat

1.0 r-:-----rr.~--r'-""'-""T'""r'f_=_:;a~~~"'"
o
C
0
100
6 1.000
08 010.000
V 50.000

06
7]
0.4

02

OL-~~~~--L-~~~~--~~~~~~~
3
10. 1 100 10' 10 2 10
Stk

Figure 2 10 Efficacité de collecte par impaction inertielle d'après Brun et al. [15]

38
 Chap JI Etude Bibliographique

Le nombre de Stokes calculé est basé sur le rayon du cylindre (dimension caractéristique du
Re 2
collecteur) et le paramètre <Il est défini par <Il = St~

D'un autre côté, Dorch et al. (cité dans[4]) ont fait le même calcul pour la sphère La figure
2 Il montre leur résultat. Sur ces résultats on peut remarquer que les évolutions de l'efficacité
de collecte pour les deux types de particules (Stokesiennes et non-Stokesiennes) sont
qualitativement semblables. Ceci conduit à définir un paramètre de similitude permettant de
réduire la dépendance de l'efficacité de collecte à une seule variable

10

0.8 Sphère

06 Re p
'T'J 0 0
04 c 16
0 128
6- 512
9 1024
02 1> 4096
Cl 8192

0 2 5
10.2 10. 1 10 10
5tk

Figure 2 Il Efficacité de collecte pour un sphère (cité dans[4])

Israel et Rosner [5] ont proposé d'utiliser un nombre de Stokes effectif

Stk eff = 'V(Re p )' Stk (284)

où 'V(Re p ) représente un facteur de correction donné par :

08

06
lJI 04
02

o~~~~~~~~~~~~~~~~~~
10. 2 101 10 4
Re p

Figure 2.12. Facteur de correction du nombre de Stokes d'après [5]

39
 Chap.ll Etude BiblIOgraphique

24 Repo dR
",(Re ) = - J e (285)
p Re p 0 Co(Re).Re

Ce facteur de correction ",(Re p ) modifie le nombre de Stokes des particules non-

Stokesiennes On peut trouver l'évolution du terme \jf(Re p ) en fonction du nombre de
Reynolds de la particule sur la figure 2.12.
Par intégration de l'équation (2.85), Wessel et Righi [16] proposent l'expression suivante
3(~ Re p -tan-I(~ Re p »
",(Re p ) = 3/2 (2 86)
c Re p

avec c = 0 158
Par la suite, Israel et Rosner [5] en utilisant le nombre de Stokes effectif ont reporté les
résultats de Brun et al sur la courbe unique de la figure suivante .

100 101 10 2 10 3
Stk eff : o/~ Re,,}.Sfk

Figure 2.13 Efficacité de collecte théorique par impaction d'après [5]

et ils ont proposé (pour Stkeff > 0.14) l'ajustement suivant

11 = (1 + 1. 25(Stko:tr - 1/8)-1 - 1,4.10-2(Stk ctr -1/8)-2 + 5, 08 .1O-s(Stk etr - 1/8) -3) -1 (2 87)

Comme dans le cas précèdent, toutes les particules heurtant la paroi sont supposées y adhérer.

Wessel et Righi [16], en s'appuyant sur le schéma d'Israel et Rosner [5], ont également calculé
les trajectoires des particules en suspension dans un écoulement perpendiculaire à un cylindre
en utilisant la méthode des différences finies. A partir de l'hypothèse d'un écoulement potentiel,
ils ont établi les équations du mouvement des particules le long des lignes de courant depuis un
point en amont loin du cylindre jusqu'au point de stagnation La comparaison entre leurs
résultats et celui de Brun et al est représentée sur la figure 2.14 .

40
 Chap./I Etude BIbliographique

v
.'

~~~+H~~·iR4 . ~.c~:+H~~~~~~~~
~~~~~~~~++~Hflffffi~-r~HHH---TR~ECTORIESMODfl
~~~~~~~H+~~~-H~i-----.~~U~~A~~~'I~~", .•
1111111111111 1. 1 IIIIIIIIIIII~

• 10 20 10 '00 200 500 1000

nombre de Stokes. Stk

Figure 2.14 Efficacité de collecte par impaction pour un cylindre circulaire [16].

Sur cette figure, on peut constater qu'en général la comparaison est bonne Les différences ne
sont sensibles que pour des valeurs élevés du paramètres cI».

Enfin. Wessel et Righi [16] ont également proposé des corrélations générales pour l'efficacité
de collecte par impaction sur un cylindre. Ils proposent pour 0.125 < Stk eff < 05 la relation

" = ~1In(8Stkeff) + ~2 (Stk eff - 0.125) + ~3 (Stk eff - 0.125)2 (288)

avec: ~I =0 01978749. fil =O. 5136545, ~3 =-0.0482858

et pour Stk eff > 0 5 :
" = (1 + ~l (Stk eff - 0.125)-1 + ~2 (Stk eff - 0.125)-2 + fi3 (Stk eff - 0.125)-3 )-1 (289)

avec: ~I =1.54424,' ~2 =-0 538013, ~3 =0.2020116

L'hypothèse de l'écoulement potentiel utilisée dans le calcul de Wessel et Righi [16] limite la
validité de ces calculs à l'écoulement autour d'un cylindre de dimension assez grande placé dans
un écoulement à forte vitesse. Pour les cylindres de dimension plus petite tels que ceux utilisés
dans les systèmes de filtration, l'écoulement est visqueux. Selon les auteurs le calcul précédent
devient alors moins précis.

Par ailleurs, Ilias et Douglas [17] ont fait les calculs numériques pour l'impaction inertielle des
particules sur un cylindre placé perpendiculairement à un écoulement pour 30<Re<40000.
Pour obtenir les trajectoires de particules, ils ont établi les équations de quantité de mouvement
des particules en utilisant la répartition des vitesses obtenue après résolution des équations de
Navier-Stokes. Leurs résultats de calcul sont comparés avec les points expérimentaux de
Wong et Johnstone (cité dans [17]), et ceux-ci sont illustrés sur la figure suivante.

41
 Chap.JJ Etude BiblIOgraphique

101~----------------------~
o réel
1)

écoulement

potentiel

,
,,'
o ,, "
o
QO~-'~--~~~----------
-
,
visqueux
0.05 0 50 5.00
Stk/2
Figure 2 15 Comparaison des efficacités de collecte calculées avec les résultats
expérimentaux de Wong et Johnstone (cité dans [17])

On peut constater sur cette figure que lorsque l'on compare les résultats expérimentaux avec
les modèles de l'écoulement visqueux l'efficacité de collecte calculée est sous-estimée, tandis
que l'hypothèse de l'écoulement potentiel donne une légèrement surestimation Dans la
pratique, le champ d'écoulement est situé entre les deux extrêmes. La résolution complète des
équations de Navier-Stokes est alors nécessaire pour obtenir l'écoulement réel.

Par la suite, en s'appuyant sur la résolution complète des équations de Navier-Stokes, Ilias et
Douglas [17] ont calculé les variations des efficacités de collecte par impaction inertielle sur un
cylindre, et la figure (2.16) illustre les résultats du calcul pour différents nombres de Reynolds
On peut remarquer que ces résultats de calcul sont en bon accord avec les expériences de
Wong et Johnstone.
Enfin, les auteurs ont également proposé la fonnulation suivante
(valable pour 30 < Re < 4.104 et 007 < q, < 5.0)

Tl = (<pJ +~).(A2.<pJ +(A 3 + A 4 .lnRe).<p2 + As.lnRe +A6)-l (2.90)

<p
Stk 4 2
avec <p = - , et Al = 1. 622 10 ,A2 = 1. 031 , AJ = 1 14 , A4 = 4 044 10-
2
As = 1.479 10-2 , A6 = 02013
Cette relation pennet d'obtenir les points avec un écart maximum de 6 4 % et un écart type de
o 08 %. Les auteurs préconisent également que cette relation peut être utilisée pour estimer
des efficacités de collecte par impaction pour l'ensemble des collecteurs dans les systèmes de
filtration.

42
 Chap Il Etude BIbliographIque

Stk/2

Figure 2.16 Efficacités de collecte de l'impaction inertielle pour différentes

valeurs du nombre de Reynolds [17]

2.2.1.7 IMPACTION INERTIELLE SUR UN FAISCEAU DE TUBES

L'impaction inertielle des particules sur un obstacle cylindrique isolé dans un écoulement est
maintenant bien connue. Les résultats des études couvrent une importante gamme de
paramètres et les corrélations nécessaires permettant d'estimer les paramètres caractéristiques
sont disponibles. Cependant, dans le cas de l'impaction inertielle sur un faisceau de tubes, on
connait peu les effets de la disposition des tubes dans un faisceau sur l'impaction inertielle
relative à un tube placé dans une rangée.

Plusieurs auteurs ont essayé d'étudier ce sujet, par exemple Choudary et Gentry [18] et Ingham
et al. [19], mais leurs résultats sont limités pour le cas d'une rangée de tubes ayant un pas
transversal adimensionné S de 5.35 seulement (S représente le pas transversal normalisé avec le
rayon du tube - figure 2.17).
Konstandopoulos et al. [20] ont étudié l'impaction inertielle sur une rangée de tubes lisses
placé perpendiculairement à un écoulement du fluide. Ils ont analysé les effets du nombre de
tube dans une rangée et leurs dispositions sur l'impaction inertielle de particules sur un tube
central dans la rangée.

43
 Chap.// Etude BIbliographique

Figure 2 17 Une rangée des tubes pour S = 4

Ils ont montré que, comme proposé par Israel et Rosner [5], l'impaction inertielle sur un tube
placé dans une rangée de tubes peut être caractérisée en utilisant le nombre de Stokes effectif
défini par :

Stk eff == 2.(dU) (291)

2 dx sp

ou. (dU)
- est 1e grad'lent d e Vitesse
. ' de stagnation.
au pomt .
dx sp
Les valeurs de ce terme adimensionné (normalisé par U ... /R) pour un tube central placé dans
une rangée de tubes en fonction du nombre de tubes Ncyl et leurs disposition dans une rangée
(cité dans [20]) sont données sur la figure 2.18.

Sur cette figure, le terme [dU] atteint une valeur maximale pour la valeur de S égale à 3
dx sp
environ et cette valeur augmente avec le nombre de cylindres (tubes) dans une rangée.
Dans le cas d'un tube isolé dans un écoulement potentiel, la valeur de [dU] est égale à 2
dx sp
(expression 2 76). Cependant, pour une rangée de tubes infinies ayant un pas transversal
adimensionné S égale à 4, la valeur de [dU] pour le tube central est d'environ 2.36
dx sp

44
 Chap 11 Etude BiblIOgraphique

J.OO r----,,......--,.---r---...-__ro-..,...-..,---..,

2.00

1.00 ' - - - - - I L - - - - ' _ - - L _ - - L _......._ ........_ .............-I

2.00 4.00 S.OO 8.00 10.00

pas transversal S
Figure 2 18 Gradient de vitesse au point de stagnation en fonction des valeurs du pas .
transversal S, pour différent nombre de cylindres N cyl [20]

De plus, les auteurs ont montré que l'efficacité de captation sur un tube dans le faisceau évolue
comme le montre la figure 2.19. A partir de ces résultats, il est donc possible de calculer
l'efficacité de collecte par impaction inertielle sur un tube placé dans un faisceau de tubes en
utilisant les valeurs de [dU] données sur la figure 2 18 et en utilisant les corrélations pour
dx op
l'efficacité de collecte par impaction inertielle (expression 2 87 à 2.90)

N..,: 1-00
S : 2.1-00
1] 0.50

Figure 2.19 Efficacité de collecte d'un tube central dans une rangée de tubes [20]

45
 Chap.JJ Etude BibliographIque

Parallèlement, Fuhs [21] a réalisé une étude expérimentale du dépôt par impaction inertielle sur
un faisceau de tubes lisses. Dans ses expériences, l'auteur utilise un faisceau de tubes composé
de 5 rangées de 5 tubes Le faisceau de tubes est placé dans une section d'essai et exposé à un
écoulement d'air chargé de particules d'ammonium fluoresceine de diamètre moyen 4 82 J.1m
Chaque campagne d'essais dure 10 minutes environ et la vitesse moyenne du gaz à l'entrée du
faisceau est de 1320 cmls ce qui correspond à un nombre de Stokes effectif de 0 32 En
mesurant la quantité de particules déposées sur chaque tube dans le faisceau, il obtient les
résultats représentés sur la figure 2.20.

Les résultats expérimentaux obtenus montrent que la quantité de particules collectées dépend
de la disposition des tubes dans le faisceau. La quantité de particules collectées sur la première
rangée est plus importante Elle diminue pour la deuxième rangée, puis augmente légèrement
sur les trois dernières rangées. L'auteur a remarqué que d'une manière générale la répartition
du dépôt dans le faisceau de tubes en arrangement aligné est la conséquence de la structure de
l'écoulement à travers celui-ci. Dans ces conditions, on peut constater que l'impact ion inenielle
se produit essentiellement sur la première rangée de tubes.

c 80
S 70 111 tube 1er
lro
'5 ml tube2e
~50
8.40 IIl1n tube 3e
G»
'0 30 ~ tUbe4e
G»
l5 20
E • tube Se
0 la
z
a
1er 2e 3e 4e Se
rangée de tubes

Figure 2.20 Quantité des particules déposées sur les tubes dans un faisceau de tubes.
en configuration aligné (cité dans [21 D.

En conclusion, l'impaction inertielle des particules sur un obstacle cylindrique, lorsqu'on peut
appliquer l'hypothèse de l'écoulement potentiel, est désormais bien connue et plusieurs modèles
classiques existent dans la littérature. De nombreux auteurs ont calculé les trajectoires de
particules dans un écoulement potentiel autour des obstacles de différentes formes. A partir de
la détermination des trajectoires, il est possible de déterminer la valeur des efficacités de
collecte par impaction.

46
 Chap.1I Etude BiblIOgraphique

Israel et Rosner [5] ont calculé les efficacités de collecte des particules sur un cylindre et une
sphère isolée dans le cas d'un écoulement potentiel en introduisant la notion de nombre de
Stokes effectif pour rendre compte du comportement des particules non-Stokesiennes.
Puis, en utilisant le schéma de calcul proposé par Israel et Rosner, Wessel et Righi [16] ont
également proposé des corrélations pour l'efficacité de collecte d'un cylindre dans l'écoulement
potentiel
Par ailleurs, I1ias et Douglas [17] ont calculé numériquement l'efficacité de collecte par
impaction inertielle sur un cylindre en utilisant les modèles de l'écoulement potentiel, de
l'écoulement de transition et de l'écoulement visqueux pour 30<Re<40,OOO.

Pour ce qui concerne l'impaction inertielle sur les tubes dans un faisceau, Konstandopoulos et
al. [20] ont montré comment l'impaction sur un tube de la première rangée dans un faisceau
était modifiée par la présence des tubes voisins Ils proposent de corriger le nombre de Stokes
par le gradient de vitesse au point de stagnation.

2.2.1.8 INTERCEPTION

Lorsque des aérosols s'écoulent au voisinage d'un obstacle cylindrique et en négligeant l'effet
des forces extérieures, nous avons w que la captation des particules à la paroi peut être
effectuée par impaction inertielle, diffusion et par interception

0-
---:::=-- - .--- -.l-~-
..=iQ.~~Q::::- __+--_

Figure 2.21 Mécanismes de dépôt sur des cylindres perpendiculaires à un écoulement [12]

Sur la figure 2.21, les particules (1) malgré leur inertie peuvent être déviées de la ligne de
courant de l'écoulement et heurtent la paroi. Pour les particules non déviées des lignes de
courant (particules 2), si elles passent à une distance inférieure au rayon de particule, elles
peuvent être interceptées. Selon Hesketh [12], ce mécanisme intervient pour un angle compris
entre 60 0 et 800 à partir du point de stagnation amont.
Pour estimer l'efficacité de collecte par interception, Ranz et Wong (cité dans [22]) ont
proposé l'expression suivante:

47
 Chap.JJ Etude BIblIOgraphIque

1
ll=l+N lIlt - - - - (2.92)
1+ N lOt

où N mt représente le paramètre d'interception .

d
N =_P (293)
mt D
c

avec d p est le diamètre de particule et De le diamètre du collecteur

2.2.1.9 SEDIMENTATION

-
Un autre mécanisme affectant le dépôt est la sédimentation qui est important lorsque les
particules sont denses et lorsque la vitesse est faible Ce mécanisme de dépôt est généralement
important dans des systèmes horizontaux et relativement statiques [12]
Sous l'action des forces de gravité g, l'équation de mouvement des particules Stokesiennes
s'écrit :
dv
m- +31t Il dpv =m g (2.94)
dt
En admettant les particules sphériques, on peut également écrire :
dv +18 Il v_ g =O (295)
dt ppd p

Si on prend comme condition intiale v = 0 à t = 0, on obtient une solution de la forme

v = g 't p (1- e-t/'t p
) (2 96)

Dans ce cas, la vitesse limite de sédimentation est donnée par

(297)

2.2.2 ADHESION ET REENTRAINEMENT

Contrairement aux mécanismes de transport, le phénomène d'adhésion des particules sur la

paroi est mal connu Les particules transportées à la paroi peuvent y être capturées ou peuvent
ensuite être réentraînées Il y a donc une probabilité pour que les particules restent à la paroi
(sans être réentraînées) La raison pour laquelle certaines particules restent et les autres sont
ré entraînées est encore mal connue. Les forces agissant sur les particules lorsqu'elles
s'approchent de la paroi jouent certainement un rôle important [8], [23], [24]. Il en est de
même pour les propriétés physiques des particules (densité, élasticité, état et condition de

48
 Chap Il Etude B,bllographlque

surface). La nature de la paroi, le type de matériau, la rugosité jouent également un rôle

important.
Lorsque des particules ont été capturées par la paroi, les particules suivantes seront collectées
par la couche de particules déjà déposées. Le mécanisme d'adhésion doit donc être analysé en
considérant deux étapes. l'interaction entre particule-paroi et l'interaction entre particules
Néanmoins. pour la plupart des situations d'encrassement particulaire. l'épaisseur du dépôt est
généralement de l'ordre de la centaine de micromètres. n est alors raisonnable de considérer
que l'adhésion des particules dépend principalement de l'interaction entre particules [8], [24]
Cependant, l'interaction entre particule-paroi ne peut être négligée lors de la phase initiale de
dépôt.

D'autre part, le processus d'adhésion dépend des caractéristiques des forces surfaciques Les
particules peuvent se maintenir fixées à la paroi grâce aux forces d'adhésion dues à l'attraction
moléculaire (force de Van der Waals), aux forces électrostatiques et capillaires [2]. [21) La
force capillaire explique l'accroissement de l'adhérence des particules d'aérosols en fonction de
l'humidité relative de l'air. En général au-dessus d'un taux de 70 % il se forme sur la paroi un
film liquide adsorbé extrêmement mince qui s'interpose entre la paroi et les particules déposées
[2] Une liaison chimique entre les particules et la couche du dépôt peut également se produire
pendant le processus d'adhésion. Dans la pratique, l'intéraction entre ces forces existe et le
phénomène devient complexe [8], [24], [25]. De plus, les modèles analytiques traitant ce sujet
sont généralement limités à des situations idéales.

Malgré la compléXité du problème, le phénomène d'adhésion est généralement traité

statistiquement en utilisant le terme de probabilité d'adhésion qui est la probabihté qu'une
particule heurtant la paroi y soit capturée. Cependant, la valeur exacte de la probabilité
d'adhésion est difficile à déterminer. L'estimation de la valeur de cette probabilité est
généralement faite à partir des connaisances de l'ordre de grandeur des forces surfaciques [8],
[24], [25]. Dans le cas où-les forces d'interactions sont supérieures aux forces aérodynamiques,
la valeur de la probabilité d'adhésion sera proche de l'unité Ceci se produit en général lorsque
la taille de particules et la vitesse sont faibles. En revanche, la valeur de la probabilité
d'adhésion est inférieure à l'unité lorsque la force aérodynamique s'excerçant sur la particule est
importante.

2.2.2.1 PROBABILITE D'ADHESION

La plupart des modèles de dépôt particulaire sont basés sur l'hypothèse d'une parfaite
captation Cette hypothèse conduit à une condition limite de concentration nulle à la paroi

49
 Chap Il Etude BiblIOgraphique

[24]. En négligeant l'effet de gravité, on peut exprimer la densité du flux de transport convectif
à la paroi (expression 2.16) par'
<Pd = hm (Cpe -C ps )
Selon Epstein [24], une approche de la probabilité d'adhésion S est celle qui pose Cps = 0 dans
l'expression ci-dessus, faisant temporairement l'hypothèse d'une adhésion parfaite et d'écrire
(298)

hm représente le coefficient de transfert de masse pour la captation parfaite et hd le coefficient

de transfert de masse réel S correspond à la probabilité qu'à une particule se dirigeant vers une
paroi d'adhérer à celle-ci En fait S corrige l'hypothèse précédente (Cps = 0) en tenant compte
de la probabilité de non adhésion des particules On peut alors exprimer la probabilité
d'adhésion par .
S = hd = hd (réel) (299)
hm hd (S = 1)

On peut donc penser que la valeur de la probabilité d'adhésion dépend des caractéristiques de
l'interaction entre les particules et la paroi qui dépend de la taille, la vitesse et les propriétés
physiques des particules, la géométrie et la température de la paroi Mais, la dépendance de la
probabilité d'adhésion avec ces facteurs est encore mal connue.

2.2.2.2 REEN:TRAINEMENT

Comme le phénomène d'adhésion, le phénomène de réentraînement est peu étudié En général,

on considére que le réentraînement se produit lorsque la force aérodynamique est supérieure
aux forces d'interactions entre les particules et la couche du dépôt. Cependant, on peut noter
que plusieurs auteurs ont essayé d'analyser ce phénomène du point de we macroscopique et!ou
. .
nucroscoplque.

Pour Kem et Seaton [3], la dénsité de flux massique de réentraînement est proportionnelle à la
contrainte de cisaillement et à l'épaisseur du dépôt e d :

(2.100)

Taborek et al. [26], [26a] ont modifié cette relation dans le cas du dépôt de particules en phase
liquide. Selon eux, la densité de flux massique de réentraînement doit être proportionnelle à la
contrainte de cisaillement et à la masse du dépôt, et inversement proportionnelle à la résistance
du dépôt au réentraînement :
m
<P =Ct _d (2.101)
r S 'Vd

50
 Chap.JJ Etude BiblIOgraphique

OÙ 't s est la contrainte de cisaillement à la paroi, 'JI d le coefficient de résistance au

réentraÎnement du dépôt, md la masse de particules déposées par unité de surface, et C une
constante empirique Comme la contrainte de cisaillement correspond au facteur de frottement,
on peut donc constater que le flux de réentraînement croît avec à la vitesse du fluide

Cleaver et Yates [27], [27a],[27b] ont montré que la sous couche visqueuse n'est pas tout à
fait laminaire, mais qu'elle est continuellement désorganisée par les remous turbulents
aléatoires qui naissent à la paroi. Ces remous peuvent produire des forces de portance qui sont
suffisament importantes pour ramener les particules déposées dans l'écoulement principal Ces
auteurs ont également constaté que la particule sera détachée de la paroi si la valeur de la
contrainte de cisaillement dépasse une valeur critique, soit :
". d 4/3> ft (2 102)
"s· p -1-'

où d p est le diamètre de particule et J3 une constante qui correspond aux forces d'adhésion
entre les particules et la paroi (ou la couche du dépôt)

2.3 BIBLIOGRAPHIE SUR L'ENCRASSEMENT DES TUBES AILETES

Le nombre de travaux traitant de l'encrassement particulaire en phase gazeuse et portant

spécifiquement sur des tubes ailetés est très réduit En général, il existe deux types d'approche
pour appréhender ces phénomènes d'encrassement. Quelques auteurs ont utilisé l'approche
simulée avec des gaz encrassants artificiels pour permettre de cerner un type de phénomène en
contrôlant de façon satisfaisante les conditions aux limites D'autres auteurs ont utilisé une
approche réaliste avec des fumées réelles (diesel, chaudière, ... ) mais la gamme de variation des
paramètres est en général plus faible Nous allons dans ce paragraphe présenter les résultats
expérimentaux de différentes études

2.3.1 GAZ ENCRASSANT ARTIFICIEL

2.3.1.1. Etude de El Shobokshy et El layat [28]

Les auteurs présentent les résultats d'une étude expérimentale du dépôt de particules sur un
tube aileté isolé, avec variation de la granulométrie de l'aérosol (relativement monodispersé) et
de la température de paroi Les particules sont produites dans un générateur d'aérosols à partir
d'une solution d'uranine La gamme de variation des paramètres est la suivante ·

51
 Chap Il Etude Bibliographique

durée des essais 90mn

vitesse de l'air 3 mis
température de paroi 25 à 65 oC
diamètre des particules 2 à 10 ~m
concentration moyenne 40 mg/m3
diamètre de tube 18mm
hauteur d'ailette 8mm
écartement d'ailette 2mm

,oo~----------------------~

80

60

40

20 oA. 298 Il •• 8 .312 le

• C • 327 Il •• 0 .338 Il

00 2
diamètre de particule d (~)
P

Figure 2.22 Efficacité de collecte pour différentes températures de surface [28]

Le tube, muni d'un réchauffeur électrique, est équipé de deux thermocouples. La concentration
massique des particules et la masse collectée sont déterminées par tluorimétrie. On définit une
efficacité de collecte qui est le rapport de la masse retenu à la quantité totale passant à travers
la surface apparente du tube. On constate (figure 2.22) que celle-ci augmente avec la taille des
particules et diminue avec la température de surface en deçà de 8 ~m d'autant plus que les

52
 Chap.II Etude Blbltographlque

particules sont petites : la thermophorèse s'oppose au dépôt, le gradient de température étant

opposé à celui de concentration

2.3.1.2. Etude de Bott et al. [29]

Cette étude porte sur l'encrassement d'un faisceau quinconcé de quatre rangées de tubes murus
d'ailettes spirales enduites d'un agent collant. Les particules utilisées (du carbonate de calcium)
sont redispersées La gamme de variation des paramètres est la suivante

durée des essais 15 à 150 h

vitesse de l'air 2.4 à 58 mis
température d'air 20 à 50°C
température d'eau 5 à 80°C
diamètre des particules 3 - 30 ~m
concentration moyenne 0.12 à 1 7 glm3
diamètre de tube 254mm
hauteur d'ailette 16mm
écartement d'ailette 2mm

0·0&
0·05

003

0005

0·002

1000 1500 2000 3000 4000 5000

Re
Figure 2.23 Evolution du facteur de Chilton-Colbum j
et du facteur de frottement f [29]

53
 Chap.1I Elude Blbltographlque

Pour différent états d'encrassement, on détermine en condition propre les évolutions du facteur
de Chilton-Colbum j = St h . Pr2/3 et du facteur de frottement f en fonction du nombre de
Reynolds du fluide (figure 2 23) On obtient des droites en coordonnées log-log dont les
pentes sont conservées.

Le suivi au cours du temps montre que j diminue faiblement tandis que f présente une
croissance asymptotique (figure 2 24) Ceci s'explique par un dépôt important conduisant au
bouchage presque complet des faces avant et arrière du faisceau A l'intérieur du faisceau,
l'échange thermique n'est que peu modifié.

005r--r--'---r-~--'---r--;--ÏÎ--r-~---r--r--'
001.

.0 .,a> .~•• _ . ';wu ....... _ • .... -

oIT

001
0001
0001
0005
0- 001.
0·001

0002 •
0-0015
~ -
.... ...__ ........... -- .._.-.
•
1
...........
~...... .
.......... : ........ -.4\-.....

Figure 2.24 Evolution de j et f au cours du temps d'après [29]

2.3.1.3. Etude de Reay et al. [30]

Ce travail porte sur l'encrassement, par une poudre de silice redispersée, d'un tube coudé muni
d'ailettes De l'air chaud ou de l'azote circule à l'intérieur du tube. Il étudie l'influence de la
vitesse des gaz, de la charge électrostatique des particules, des températures et du sens du flux
de chaleur. La gamme de variation des paramètres est la suivante:

54
 Chap JI Etude BiblIOgraphique

durée des essais 8h

vItesse de l'air 20 mis
température gaz encrasst. 20°C
température aIr chaud 70/160°C
température azote, LN 2 -50/-lO o C
dIamètre de partIcule 50%< Illm
concentration 0.5 g/m3
diamètre de tube ISmm
hauteur d'ailette 17.5 mm
écartement d'aIlette 2mm

On étudie l'influence la charge électrostatique en réglant le fonctionnement du neutraliseur

(source de radiations ionisantes) On enregistre la décroissance du coefficient global de
transfert thermique au cours du temps (figure 225) et on détermine la masse déposée par unité
de temps.

-2 -1
7 HG w. m .1 -------------------,
-- ._- -- r-- __ _
T
.---.---.-- --- --.
---._--- - - -e _ _ _ ,..

,
J

°O~--------IOO~---------2·~----------l~~---------'~OO--------~SOO

temps mn

Figure 2.25 Coefficient global de transfert thermique au cours du temps [30]

Le flux net de dépôt augmente avec la charge électrostatique, il diminue avec la vitesse de
passage entre les tubes. L'effet de la thermophorèse n'est pas mis en évidence.
En début d'essai, il se forme une couche uniforme d'épaisseur 0.5 mm, ensuite le dépôt
s'épaissit à l'avant du tube et sur l'extrémité de l'ailette dans le sillage de recirculation (figure
2 26), enfin on observe un pontage entre les ailettes.

55
 Chap Il Etude Btbllographique

3 mm
--
~I
l
--
(J
~ '} ---
\

coupe transversale
--
coupe longitudinale

Figure 2.26 Localisation du dépôt d'après [30]

2.3.2 GAZ ENCRASSANT REEL

2.3.2.1. Etude de Reay et al. [30]

Le même tube aileté que celui de l'étude précédente est monté sur la conduite d'échappement
d'un moteur Diesel. La gamme de variation des paramètres est la suivante'

durée de l'essai 64 h
vitesse des fumées 20 mis
température de fumée 360°C
température d'air 35 à 240 oC
diamètre de particule 90%<1~m
concentration -08 glm3

L'observation du dépôt en fin d'essai met en évidence le pontage presque complet des
premières ailettes du côté froid du tube d'entrée. Ce dépôt préférentiel peut s'expliquer par la
formation d'un film collant, la température du point de rosée acide étant de 130°C. Le dépôt
est léger, pelucheux et tombe facilement. Le reste du tube est recouvert d'un dépôt noir
d'épaisseur inférieure à 0 1 mm, plus tenace bien que cassant ce peut être un film de condensat
formé en début d'essai lorsque le tube était froid.

56
 Chap II Etude BIblIOgraphIque

2.3.2.2. Etude de Lokshin et al. [31), [32]

Les auteurs étudient le comportement d'un faisceau de tubes à ailettes spirales soumis à un
encrassement par des cendres provenant d'une chaudière à charbon Ils font varier la vitesse du
gaz, l'arrangement, l'écartement des tubes ainsi que les paramètres d'ailetage La gamme de
variation des paramètres est la suivante :

durée des essais 05 à 15 h

vitesse de l'air 2 à 13 mis
diamètre de particule 46%< 30 J'm
concentration 36 glm3
diamètre de tube 32 mm
hauteur d'ailette 10mm
écartement d'ailette 10 mm

Les mesures de la résistance d'encrassement (figure 2.27) montrent que celle-ci est plus
importante pour les arrangements en ligne. L'influence de la vitesse est plus marquée pour les
tubes en ligne.

Pour un faisceau en quinconce Rf augmente avec le pas longitudinal car la contraction des
écoulement dirigés en diagonale et balayant l'arrière des tubes tend à disparaitre les vitesses
diagonales diminuent et favorisent le dépôt.
Pour un faisceau en ligne cet effet est peu marqué. Rf augmente avec le pas transversal
Le temps de stabilisation de la couche de dépôt est proportionnel à la vitesse de l'écoulement
Les observations du dépôt (figure 2.28) montrent que'
- celui-ci augmente de la première à la cinquème rangée de tubes,
- dans le cas des tubes en ligne, il y a une zone de stagnation entre les tubes voisins d'une
rangée longitudinale, donc un dépôt à l'avant et à l'arrière. Pour les tubes en quinconce, la
région arrière est balayée par des écoulements transverses de gaz. Ceci explique la différence
de la quantité déposée selon l'arrangement,
- pour une disposition en ligne, l'épaisseur du dépôt augmente à ravant avec la vitesse jusqu'à
2 à 4 mis puis disparait au delà de 6 mis; à l'arrière des tubes, il augmente jusqu'à 6 mis On
n'observe pas cette dépendance avec la vitesse dans le cas de tubes en quinconce.

57
 Chap Il Etude BIbliographique

-1
loi .m2.~ l 1"
quinconce
1.~ 1
1

T
ne de l'angée
1

,
~
6-1
x-z

,
1
.-1
0- •
• a-S
\ ~-I
1
.-10
1
\ l~ 0

, ~ ~ .. . -'- ~

.~ -.. --
O~ __ ~ __ ~ ______- L_ _ ~~ , ~ li-.6... ~ a......:
~ ~

•
1 J 1 3 -1 1 1 1 1 1 n -1
m.s m.s

Figure 2.27 Résistance d'encrassement en fonction de la vitesse et de l'arrangement

du faisceau d'après [31].

Figure 2.28 Localisation du dépôt d'après [31]

58
 Chap /1 Etude BIbliographIque

L'intensité de l'encrassement est proportionnelle au facteur d'ailetage qui est le rapport de la

surface ailetée totale à celle du tube nu. Ceci montre le rôle important que jouent les
1

caractéristiques de l'ailette (hauteur, espacement, inclinaison)

2.3.2.3. Etude de Lakshin et al. [33]

Ce travail concerne l'étude des caractéristiques thermiques et aérodynamiques d'un

économiseur de chaudière à charbon équipé de tubes aHetés en arrangement quinconcé. La
gamme de variation des paramètres est la suivante :

vitesse des fumées 42 à 52 mis

températu re de fumée 345 à 560°C
température d'eau 220 à 300°C
diamètre de tube 32 mm
hauteur d'ailette 10mm
écartement d'ailette 10mm

La résistance d'encrassement varie en fonction de la vitesse des gaz (figure 2 29) L'examen du
dépôt après 2000 heures de fonctionnement montre que celui-ci, léger et friable, se forme
surtout sur les bords de fuite L'épaisseur augmente dans la profondeur du faisceau Les tubes
situées à l'avant sont presque couvert de cendres, ceux à l'arrière sont plus propres et le dépôt
ne s'étend pas au delà des limites des ailettes. Pour une vitesse de 5 mis on note les formes
moyennes typiques des dépôts sur tube lisse et aileté (figure 2 30).

-1 2
W •m . K

o.OOi
~
D,oru " ~
.......
~ kc ~
0.009
I.IOJ
~
r-a- --
fI "1 '.' 'II 1/1 S,IJ s,z mis

Figure 2.29 Résistance d'encrassement en fonction de la vitesse [33]

59
 Chap JI Etude Blbltographlque

Il ,
lei
- -

Il 1 alleté li) llsse b)

Figure 2.30 Formes typiques du dépôt d'après [33]

2.3.2.4. Etude de Silvestrlni [34]

Cet auteur détermine l'effet de la configuration (hauteur, densité et type d'ailette) sur le
développement de l'encrassement par des fumées de diesel. Pour chaque essai on utilise deux
maquettes constituées de quatre tubes à ailettes lisses ou dentelées ; une est refroidie. l'autre
non. La gamme de variation des paramètres est la suivante .

durée des essais 40 à 70 jours

température de fumée 260 à 400°C
température du pomt de rosée acide 117°C
température de paroi 104 et 246°C
hauteur d'ailette 16 et 19 mm
nombre d'ailettes par mètre 118 - 197

On ne constate aucun dépôt sur les tubes non refroidis. Les tubes à forte densité d'ailette
s'encrassent plus . le bouchage complet de l'espace inter-ailette apparait pour une densité de
197 , mais non pour une densité de 118 ailettes par mètre. Pour les tubes à ailettes lisses,
l'épaisseur du dépôt est la même sur le tube ou sur les ailettes en début d'encrassement , puis
celle-ci atteint un maximum sur les ailettes tandis que l'espace inter-ailettes se remplit de suie à
partir des tubes (figure 2 31).

60
 Chap Il Etude BiblIOgraphique

1 II III

Figure 2 31 Croissance du dépôt d'après [34]

Pour les ailettes dentelées on observe le même mécanisme, avec en plus un pontage localisé
aux extrèmités. La forme des courbes de l'encrassement au cours du temps est caractéristique
(figure 2.32). Pour une même densité et une même hauteur, les tubes à ailettes lisses ou
dentelées se comportent de la même façon.

l:~
,1 1

•
épaisseur mm • face avale
• 1
face amont ~
1 --l ......,

L
1

•
V
/
, /
,
.k:::::
.,
l'
v

Il
1
,\
1
01
1
0'
---r
l temps Jo urs
1

Figure 2.32 Evolution de l'encrassement au cours du temps d'après [34]

2.3.2.5. Etude de Casper et al. [35]

Cet article expose les résultats de l'évolution en fonction du temps de l'épaisseur du dépôt et
des performances thermiques pour un échangeur de chaleur à tubes ailetés refroidi par une
circulation d'eau en arrangement quinconcé sur des fumées de diesel La gamme de variation
des paramètres est la suivante

61
 Chap.// Etude BIblIOgraphIque

vitesse des fumées 3 à 9 mis

température de fumée 188 à 304°C
température d'eau 19 à 52°C
hauteur d'ailette 19 mm
nombre d'ailettes par mètre 118

L'épaisseur du dépôt à l'avant et à l'arrière des tubes est déterminée par photographie, elle
augmente rapidement puis se stabilise (figure 2.33). Sur la face avant on observe une substance
épaisse, résistante au lavage (résidus d'hydrocarbures) Sur les faces arrières se dépose une
couche de cendres légères et friables. Après une dizaine de jours de fonctionnement, l'espace
inter-tubes se réduit à 2 mm mais sans qu'apparaisse de pontage
Le coefficient d'échange de chaleur côté fumées chute rapidement puis prend une valeur
asymptotique après 200 heures environ (figure 2.34).

épaisseur mm

20~------r-------~------~------~

, $

10

10 1$ 20
temps Jours

Figure 2.33 Evolution de l'épaisseur du dépôt au cours du temps d'après [35]

•
~
•

•
•
100 ~.

--- .. _..1

temps h

Figure 2.34 Evolution du transfert thermique côté externe au cours du temps [35].

62
 Chap.ll Etude BiblIOgraphique

2.3.2.6. Etude de Roberts et Kubasco [36]

Les auteurs s'intéressent à l'influence de la température de paroi et de la vitesse des fumées sur
l'encrassement de tubes ailetés de récupérateurs sur turbines à gaz La gamme de variation des
paramètres est la suivante .

durée des essais 30 à 240 h

vitesse des fumées 6 à 12 mis
température des fumées 455°C
température de paroi 21 à 150°C
diamètre de particule 70%<1.4 !lm
concentration o 005 à 0 8 glm3
diamètre de tube 15.9 mm
nombre d'ailettes par mètre 71 à 350

On ne note pas de différence entre les rangées quinconcées et celles en ligne On observe des
pontages pour des densités de 138 et 280 ailettes par mètre mais pas pour des den sites de 71
ou 91.

chute d'efficacité Y.

o
• " H JO 40 St .. 70 10 .. 100
temps h

Figure 235 Chute de l'efficacité de l'échangeur d'après [36]

Le dépôt présente en général une double couche : la première, composée d'hydrocarbures, est
dure et adhérente.; la seconde, qui est sèche et qui s'écaille, est formée de particules de suie de
taille inférieure au micromètre On note une diminution asymptotique de l'efficacité de l'échange
thermique (figure 2.35).

63
 Chap.Jl Etude BIbliographique

Les auteurs étudient également l'auto nettoyage de tubes préalablement encrassés que l'on
place dans un four en atmosphère oxydante ou neutre avec balayage par un gaz chaud La
température à partir de laquelle-cet effet apparait dépend de la vitesse des gaz

2.3.3 CONCLUSION

Cette étude bibliographique montre que les travaux expérimentaux portant sur l'encrassement
des tubes ailetés d'échangeurs thermiques sont peu nombreux mais variés En effet, les
géométries étudiées, les conditions thermiques et hydrauliques, les spectres granulométriques
des aérosols encrassants, les dispositif expérimentaux et les moyens de mesure sont très divers
Malgré cela, il tend à se dégager certains recoupements en ce qui concerne . l'effet de la
thennophorèse qui est mis en évidence, la localisation du dépôt et les conditions de dépôt dans
le cas des fumées d'échappement de diesel
On mesure l'importance du chemin qu'il reste à parcourir pour améliorer la compréhension des
phénomènes et par la même tenter de résoudre le problème de l'encrassement des échangeurs
de chaleur.
Par la mise en oeuvre d'un dispositif expérimental pennettant le contrôle précis de tous les
paramètres, nous nous proposons d"étudier le comportement en conditions encrassantes,
d'abord d'une géométrie simple un tube isolé en écoulement perpendiculaire, puis d'une
géométrie complexe une maquette d'échangeur à tubes ailetés alignés

64
 Chap III Partie Expérimentale

Chapitre III
PARTIE EXPERIMENTALE

Dans ce chapitre, on décrit le dispositif expérimental permettant de déterminer les

performances d'une maquette d'un échangeur de chaleur en conditions encrassantes. On décrit
en détail les sections d'essais, les moyens de mesure mis en place pour réaliser l'étude
expérimentale Le dispositif expérimental utilisé (boucle Gazpar) comprend une maquette d'un
échangeur thermique constitué d'un tube unique ou d'un faisceau de tubes placé dans une veine
d'air parcourue par un aérosol. Les particules encrassantes injectées dans le circuit d'air sont
générées par pulvérisation ultrasonore dans un générateur de gouttelettes de solution saline.
Celles-ci sont ensuite introduites dans le circuit d'air chaud où elles sont séchées avant de
traverser la section d'essai.
En suivant l'évolution dans le temps des performances thermiques de l'échangeur, on peut
déterminer l'évolution de la résistance d'encrassement pour une classe de taille de particules et
pour une vitesse données La mesure de la masse du dépôt recueillie en fin d'essai permet
d'obtenir une intégration du phénomène sur la durée de l'essai Le suivi de la résistance
d'encrassement donne des informations sur la cinétique d'encrassement.

3.1 PRODUCTION DES PARTICULES

Les particules encrassant es utilisées pour l'étude sont produites par pulvérisation ultrasonore
dans un générateur d'aérosols qui permet d'obtenir à partir d'une solution saline un débit de
gouttelettes de taille presque uniforme Le séchage de ces gouttelettes permet ensuite d'obtenir
un résidu sec de particules proche de la monodispersion Boulaud [37] a souligné que ce
système de production de particules est le seul capable de produire un aérosol stable de
granulomètrie bien contrôlée et de forte concentration
La dimension d'une particule résiduelle d p dépend de la dimension de la gouttelette d'origine
dg et de la concentration Cs du soluté dans la solution (en pourcentage massique), si Ps est la
densité du soluté ces grandeurs sont reliés par ' .
d 3P Ps -- d3g Cs (3 1)

On constate donc qu'il est possible par des conditions de fonctionnement identiques de
modifier le spectre des particules résiduelles en faisant varier Cs Mais comme le débit de

65
 Chap III PartIe Expérimentale

gouttelettes et dg restent constants, le choix d'un faible d p conduit à une diminution très rapide
du débit massique du résidu sec

3.1.1 GENERATEUR D'AEROSOL

Boulaud et Bourbigot [38] on réalisé une étude préliminaire d'un pulvérisateur ultrasonore
dont le principe est donnée figure 3 1. La vibration d'une céramique piézoélectrique immergée
dans un liquide produit à l'interface liquide-gaz une fontaine ultrasonore, l'atomisation du
liquide au sommet de cette fontaine a pour conséquence la formation d'un brouillard Le débit
de gouttelettes d'aérosol entraînées par le débit Qd d'air de dispersion est injecté dans
l'écoulement principal. L'évaporation du solvant est réalisée à la température ambiante La
granulométrie et la concentration des résidus secs sont mesurées par impacteur cascade
(paragraphe 3.3.1)

Le générateur d'aérosols utilisé est équipé d'une seule céramique dont les paramètres de
fonctionnement sont les suivants :
- source acoustique : fréquence et puissance délivrées,
- solution à générer: pression de vapeur saturante, tension superficielle, viscosité,
concentration en soluté,
- hauteur de liquide dans le pôt de pulvérisation qui agit sur la focalisation de la source, sur la
surface et sur l'entrainement des gouttelettes,
- air de dispersion: point d'injection, débit, pression, et température.

L'optimisation de l'appareiIJage a été réalisée par lM. Grinot [1] Il en ressort que la quantité
d'aérosol produite est sensiblement proportionnelle à la puissance appliquée à la céramique La
fréquence optimale est déterminée selon la forme de la fontaine et l'intensité d'atomisation
Lorsque la concentration en soluté augmente, la pression de vapeur saturante diminue, la
viscosité et la tension interfaciale augmentent , en conséquence le débit de gouttelettes
produites diminue La production augmente lorsqu'on diminue la distance entre l'injection et la
surface du liquide, avec une limite car il faut éviter la dispersion des grosses gouttelettes. Il
existe un maximum de production pour un certain débit de dispersion , au delà, la turbulence
précipite les gouttelettes sur la paroi.
Ce générateur fonctionne de façon optimale pour une fréquence d'excitation de 775 kHz, une
puissance de 200 W, une distance entre l'injection d'air et la surface de la solution de 2 cm, une
hauteur de liquide de 9 cm au dessus de la céramique et un débit d'air de dispersion de 30
Nl/min. Toutefois, pour des essais de longue durée (supérieure à 24 h) la puissance doit être
réduite jusqu'à 70 % de la puissance maximum Dans ces conditions, le débit massique de
liquide pulvérisé et entrainé est en moyenne de 0.370 I/h pour une concentration en soluté de

66
 Chap III PartIe Expérimentale

10 gIl et de 0.240 I/h à 75 gIl Enfin la céramique utilisée fixe la taille des gouttelettes obtenues
à environ 8 ~m Pour Cs = 20% et Ps = 2, la dimension de la particule résiduelle est donc de
37 ~m

_ _ _ _ _ _ Débit d'aérosol
oar

l1 ~~~~~~~:on Distance entre la surface du

o 1// liquide et l'In)ecbon du flUide

d/ { ' nature./' solvant

h st ......... soluté
1-----1---- Solutron • température

~~=::q.. ___ ~~-:""1 PUissance P

• concentration soluté Cs

...........;........., Fréquence tH

Figure 3 1 Schéma du pulvérisateur ultrasonore

3.1.2 NEUTRALISATION DES CHARGES ELECTRIQUES

L'inconvénient majeur du générateur ultrasonore est de produire des gouttelettes qui sont
électriquement chargées Pour éviter des effets électrostatiques indésirables, il est nécessaire de
neutraliser le brouillard formé Liu et Pui [39] ont montré que l'on pouvait atteindre ce résultat
en exposant l'aérosol à une source de radiations ionisantes Pour ce faire, les gouttelettes de
l'aérosol traversant une chambre de neutralisation équipée de 12 sources ORIS de 4
microcuries d'Américium 241 placées à la sortie du générateur et réalisée selon les
spécifications de Charuau [40] du Centre d'Etudes Nucléaires de Saclay
La source de radiations (Am 241 ) produit des ions bi-polaires le brouillard de gouttelettes est
amené à un état électrique proche de celui de Boltzmann pour lequel la charge totale est nulle
et le nombre de particules chargées positivement égal au nombre de particules chargées
négativement par classe granulométrique

3.2 DISPOSITIF EXPERIMENTAL

Le dispositif expérimental schématisé dans la figure 3 2 se compose de la section d'essai, d'un

circuit d'air chaud et d'un circuit d'eau de refroidissement. L'air sec et propre est fourni par des
compresseurs suivis d'un assécheur d'air et d'un filtre micronique Le contrôle du débit d'air se

67
 Chap. III Partie Expémllentale

fait par un débimètre massique et par une vanne électropneumatique Un réchauffeur et un

régulateur permettent de contrôler la température de l'air à l'entrée de la section d'essais
L'ensemble de dispositif expérimental est isolé thermiquement

CIIICVIT 0' EAU

elreulatel.ll"
alero HP acqulslllo~

CIRCUIT D'AIR
rol~tre

cUblh,Hre
r ric:baulteur
•

SECT 1CIl D' E:sSA 1

neul"llsal Ion

asdcheur d'air

pul"trlsaleur
•
dtb\lM~re

Figure 3 2 Schéma général de l'installation

Nous avons choisi comme sel le sulfate de sodium (Na 2 S0 4 de masse volumique de 2680
kg/m 3) du fait de sa bonne solubilité dans l'eau, de sa bonne stabilité à l'état cristanin et de sa
faible corrosivité Un générateur de particules produit des gouttelettes qui sont entraînées par
de l'air de dispersion Le débit massique de ce dernier est contrôlé à l'aide d'un régulateur
L'aérosol est neutralisé par une source ionisante avant d'être introduit et séché dans le circuit
d'air principal. L'air chargé de fines particules est homogénéisé dans une conduite (de diamètre
de 90 mm et d'une longueur de 1250 mm) calorifugée. Il traverse ensuite un convergent-
divergent avant d'entrer dans la section d'essai La granulométrie de l'aérosol est mesurée par
un impacteur cascade Andersen qui détermine la distribution en taille des particules pour des
diamètres aérodynamiques compris entre 0 35 p.m et 10 ~m.

Le débit massique d'air est maintenu constant pendant un essai , la pression de l'air dans le
circuit est la pression atmosphérique. Le circuit interne d'eau de refroidissement de l'échangeur
fonctionne en boucle fermée L'ajustement du débit est réalisé par un circulateur et par des
vannes de réglage manuel, un rotamètre permet de le mesurer. Un réchauffeur et un régulateur

68
 Chap. III Partie Expérimentale

régulateur permettent d'obtenir une température d'eau contrôlée En sortie de la section

d'essais, le refroidissement est assuré par l'intermédiaire d'un échangeur L'ensemble de
l'échangeur d'essais et des divergents d'entrée et de sortie d'air est placé dans une enceinte
remplie de vermiculite pour assurer l'isolation thermique.

3.2.1 INSTRUMENTATION

Le débit d'air du circuit principal est mesuré par un débitmètre massique Hastings et contrôlé
par intermédiaire d'une vanne électropneumatique Adareg grâce à un régulateur Eurotherm Le
débit d'air de dispersion est contrôlé par un régulateur de débit massique Tylan capable de
fournir un débit compris entre 0 et 1.8 kg/h. Les deux débitmètres, étalonnés par les
constructeurs, ont tous une précision annoncée inférieure à ±1% de la mesure

Le chauffage de l'air est assuré par un réchauffeur électrique de 2 kW contrôlé par un

régulateur Eurotherm. Le suivi des performances thermiques de l'échangeur au cours du temps
est assuré par des mesures de températures. Nous avons pris soin d'utiliser des thermocouples
fer-constantan de 0 5 mm de diamètre plongeant de quelques centimètres dans le fluide. Leur
erreur absolue est de l'ordre du dixième de degré. L'entrée et la sortie de l'échangeur côté air
sont équipées chacune de quatre thermocouples placés au niveau des divergents
Pour éviter les erreurs dues à l'accumulation de particules sur les thermocouples, lorsque l'on
travaille avec des grosses particules, il est nécessaire de surveiller l'évolution des températures.
Sur le circuit d'eau sont implantés un thermocouple simple à l'entrée de la section d'essais et un
thermocouple différentiel entre l'entrée et la sortie.
L'instrumentation est complétée par un système d'acquisition : multiplexeur HP-3421 et micro-
ordinateur HP-9807 qui permet l'enregistrement de toutes les données (moyennes sur 10
acquisitions toutes les 10 minutes) et leur traitement ultérieur.

3.2.2 SECTIONS D'ESSAIS

3.2.2.1 TUBE SONDE

Le tube sonde est un tube lisse cylindrique de 20 mm de diamètre et de 40 mm de longueur

placé dans une section d'essai rectangulaire dont les dimensions sont données sur la figure 3 3
La section de passage à l'entrée de la section d'essais est de 40 mm x 80 mm, ce qui donne un
facteur de blocage de 0.25 (on définit le facteur de blocage comme un rapport entre le
diamètre du tube et la hauteur de la section de passage du fluide).

69
 Chap. III PartIe Expérimentale

aérosols

Tube

Figure 3 3 Tube sonde placé dans la section d'essai

Tube Sonde

Diamètre interne 18 mm
Diamètre externe 20 mm
Longueur 40 mm
Maténau: ACier Inoxydable "NS22S"

masse volumique 7900 kglm3

chaleur spécifique 0.51 JlkgK
conductlVlté thermique 16 W/mK
diffuslVlté thermique 0.04 m2/s

Tableau 3.1. Caractéristiques du tube.

Le tube est refroidi par une circulation d'eau d'un débit de 700 lIh, ce qui pennet de maintenir
une température relativement constante sur ta paroi. Les caractéristiques géométriques de ce
tube sont données dans le tableau 3.1.

3.2.2.2 FAISCEAU DE TUBES AILETES

La figure 3 4 donne le schéma de la maquette de l'échangeur de chaleur utilisé, elle est

constituée de 10 tubes munis d'ailettes circulaires planes d'épaisseur constante Les ailettes sont
usinées dans la masse du tube pour éliminer le problème de la résistance de contact.

70
 Chap III PartIe expérimentale

L'arrangement des tubes dans le faisceau est en ligne Les tubes sont reliés en série et refroidis
par une circulation d'eau interne de débit de 300 Vh Le coefficient d'échange interne sera donc
élevé Les conditions réelles d'écoulement existant dans l'échangeur sont simulées grâce à la
présence de demis tubes refroidis de géométries identiques à celles des tubes centraux Les
parties supérieure et inférieure de l'enveloppe contenant les tubes sont démontables, ce qui
facilite l'accès au dépôt. Les caractéristiques géométriques de cet échangeur sont données dans
le tableau 3 2.

10 rangées de tubes

Figure 3 4 Faisceau de tubes ailetés

Faisceau de tubes ailetés

Tubes : diamètre interne 18 mm

diamètre externe 20 mm
pas longitudinal 40 mm
pas transversal 40mm
matériau ACIer Inoxydable
Ailettes : "NS22S"

diamètre 36 mm
épaisseur 1 mm
hauteur 8mm
nombre d'ailette lm tube 100
espacement d'ailette 10mm
surface d'échange totale 167007 nun2
(coté air)
matériau Acier Inoxydable
"NS22S"

Faisceau. encombrement: L x H x 1 40x80x400

secuon minimale de passage 1472 mm2

Tableau 3 2. Caractéristiques du faisceau de tubes ailetés

71
 Chap III PartIe Expérimentale

3.3 GRANULOMETRIE DES PARTICULES

La distribution granulométrique correspond à la variation avec le diamètre de la probabilité que

possède une particule prise au hasard dans une population d'avoir ce diamètre La figure 3.5 en
montre les différentes formes de présentation. Le spectre de fréquence des aérosols se
rapproche en général de la loi logarithmico-normale que l'on représente en portant la loi de
probabilité normale en fonction du logarithme du diamètre des particules.
Le diagramme cumulé est une représentation très commode puisque les méthodes
expérimentales de mesure donnent souvent par rapport à la population totale, la proportion en
nombre. en surface ou en masse des particules de taille supérieure ou inférieure à un diamètre
L'ordonnée est donc proportionnelle à la fréquence cumulée pour un diamètre donné et
exprimée en pourcentage de la population totale

Pour un aérosol distribué suivant une loi log-normale. le diagramme cumulé est une droite- Le
diamètre médian correspond à une fréquence cumulée de 50% et la dispersion de l'aérosol est
caractérisée par son écart type géométrique déterminable graphiquement Plus ce nombre sans
dimension est voisin de 1, moins l'aérosol est dispersé quelle que soit la grandeur de référence

"cumulés
échelle 1 1 1
log probahtés
99 1
9S
95
~.t
,--, 1
90 .....·V
sa Il '/ fréquences
70 /'
./
60 v
50
40 /1'1
30
/' 1 1
20
10
5
/' , 1 6. • diamètre à 84.1 "
diamètre à 50" -

:a 1
1 1 courbes de GAUSS
0.5 échelle
Q2 '01 " 10 loganthmiQue
0' diamètre en microns

distnbutJon loganthmico-nonnale cumulée

Figure 3.5 Présentation des résultats en granulométrie [40]

72
 Chap III PartIe ExpérImentale

nombre, surface, masse des particules. Un autre avantage de la représentation tient dans le fait
que si les extrèmités de la courbe cumulée s'éloignent plus ou moins de la droite logarithmico-
normale, les déterminations se font sur la partie rectiligne sans provoquer d'erreurs sensibles

3.3.1 MESURE DE LA GRANULOMETRIE

On présente dans ce paragraphe d'une part le moyen de mesure mis en oeuvre pour détenniner
la granulométrie de l'aérosol et d'autre part les résultats moyens de mesure de taille de
particules. Pour obtenir des ordres de grandeur de la taille des particules, on utilise un
impacteur cascade (Andersen Mark II) facile de mise en oeuvre et d'interprétation Les bases
de la technique de l'impacteur cascade comme moyen de mesure granulométrique peuvent être
trouvées dans une revue bibliographique de Fuchs [41]

Les particules y sont classées selon leur diamètre aéraulique da (tel que l03d a2 =p p d p 2 )
L'appareil est constitué d'une série de buses d'accélération et de plateaux de collecte : la taille,
la vitesse et la géométrie du jet étant contrôlées, l'inertie des particules au-dessus d'une certaine
taille devient supérieure à la force de traînée' ces particules quittent alors les lignes de courant
d'écoulement, défléchies pour aller heurter les supports de collecte où elles se déposent

• ••
•• • •
~"F'TT~~i~.w-r7~'7'7"'7"7"'7"""""
........

Figure 3 6 Principe de l'impaction en cascade

L'appareil comprend en général huit plateaux et un filtre total et permet d'obtenir, par pesée
des supports de collecte, le diagramme cumulé et éventuellement la concentration massique Le
diamètre des buses est compris entre 2.54 mm pour le premier étage et 0254 mm pour le

73
 Chap. III Partie Expérimentale

dernier. L'efficacité de collecte d'un étage pour des particules d'une taille donnée et un débit
fixé dépend du nombre de Stokes en l'absence de réentraÎnement ou de rebond. Ce nombre est
égale à (expression 2.73) .
d 2V
Stk = Pp p 0
(3 2)
181ld

où d est le diamètre de la buse.

La courbe idéale de séparation d'un étage devrait présenter un seuil net pour la valeur de Stk
considérée, ce qui n'est pas le cas dans la réalité (figure 3.7) En conséquence, on caractérise
habituellement un étage donné par le diamètre de coupure d so ou Stk~~2 qui correspond à une
efficacité de 50 % . L'appareil est en général étalonné par le constructeur avec des aérosols
monodispersés

~R 10~--~----~----~--~~--~

oa

01

02

o 10 Stk

Figure 3.7 Efficacité de collecte d'un étage d'impacteur

3.3.2 MISE EN OEUVRE DE L'APPAREIL

Le prélèvement se fait en sortie de la conduite d'homogénéisation au niveau d'un convergent,

ce qui permet de positionner parfaitement la buse dont le diamètre a été calculé pour respecter
les conditions d'un prélèvement isocinétique L'appareil doit être traversé par un débit
volumique de 28.32 Vmin, il faut donc qu'il travaille à température et pression connues et fixes
Pour que la mesure ait un sens, il faut d'autre part, travailler à une température supérieure à la
température de rosée du fluide chaud soit 20.3°C pour une humidité absolue de 15 glkg
Pour éviter les réentraÎnements dans l'impacteur, la masse maximale déposée par plateau ne
doit pas dépasser 10 mg ce qui permet de fixer le temps du prélèvement Le temps de

74
 Chap. III PartIe Expérimentale

prélèvement varie selon la concentration de la solution pulvérisée • il est limité par la masse
déposée et selon Boulaud [37] les particules devraient se déposer sur 2 à 3 plateaux après
tarage Ils sont pesés après échantillonnage avec une précision de 0 01 mg

Pour respecter ces contraintes sans trop compliquer le système. nous avons choisi de placer
l'ensemble de mesure. c'est à dire l'impacteur et sa pompe dans une enceinte thermostatée à
50°C. Le débit massique qui traverse l'appareil correspond au débit volumique de
fonctionnement à cette température et fixe le diamètre de la buse de prélèvement Pour
minimiser les risques de dépôt dans la tuyauterie de liaison avant l'impacteur, nous avons placé
l'impacteur dans une enceinte à la verticale du point de prélèvement Le refroidissement dans
cette tuyauterie s'effectue naturellement. en contrôlant simplement en cours de prélèvement la
température à l'entrée de la pompe d'extraction qui est toujours restée voisine de 50°C

3.3.3 RESULTATS DES MESURES

En s'appuyant sur les résultats de Soole [42], Boulaud [37a] a mis au point un logiciel de calcul
du diamètre médian massique de la population et de son écart type Nous avons utilisé ce
logiciel qui prend en compte les caractéristiques réelles de coupure des différents étages pour
recalculer par itérations successives la loi log-normale de répartition.
La distribution de la taille de particules est déterminée, par un prélèvement d'une partie
d'aérosol à l'entrée de la section d'essais. en utilisant l'impacteur cascade d'Andersen A partir
de données recueillies sous forme de fractions finies sur les différents étages de l'impacteur, on
peut calculer le diamètre médian massique et la déviation standard géométrique de la loi log-
normale.

Pour les deux concentrations en sel retenues (10 g/I et 75 g/I). nous avons obtenu les résultats
suivants:

Concentration (g/l) da (J.1.m) crg

10 2.21 172
221 1.65
2.20 168

75 394 153
4.05 148
3.97 148

75
 Chap. III Partle expérimentale

40
da moyen = 2.20 ",m
35
.... éc. type = 1.68
"'30
-~ 25
r:r
lli
~
20
c 15
-e
.2
CJ
10
- 5

0
~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
d d d d d d 0
diamètre aéraulique, da ( ~ )

Figure 3.8 Distribution en masse des particules de diamètre aéraulique moyen de

da = 2.20 J.1m

4O.-------------------------------------------------~
da moyen • 3.97 ",m
35
.... éc. type = 1A8

-
"'30
!r:r 25

)20
oc 15
=u
,g 10
5
o +-~~-+--~--+-~---+--~--+

~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~
d d d d d 0 d
diamètre aéraulique, da ( ~ )

Figure 3.9 Distribution en masse des particules de diamètre aéraulique nioyen de

da = 3.97 J.1m

76
 Chap. III Partie Expénmentale

Dans l'impacteur, l'essentiel du dépôt est recueilli sur deux plateaux de plus en plus proche de
la sortie quand la concentration baisse L'aérosol obtenu n'est donc pas monodispersé
Toutefois, si l'on admet comme critère de "monodispersion pratique" (puhs et Sutugin [43])
log Og ~O 2 (Og ~ 1.22), l'écart reste peu important

Nous dénommerons dans la suite les deux tailles de particules obtenues par leur diamètre
aéraulique moyen soient 2.20 /lm et 3.97 /lm . Le détail de la distribution en masse des
particules de ces deux aérosols est représenté sur les figure 3.8 et 3 9
Nous n'avons pas effectué d'autre mesures de contrôle, même après le changement de la
céramique du pot de pulvérisation du fait de la bonne reproductibilité des résultats
expérimentaux

Il est bon de remarquer en conclusion, que l'originalité du dispositif expérimental réside dans la
mise en oeuvre d'un appareil de production d'aérosols faiblement dispersés mais reproductibles
iJ devrait permettre J'étude des mécanismes élémentaires de dépôt

77
 Chap IV Etude du dépôt sur un tube Isolé

CHAPITRE IV - ETUDE DU DEPOT SUR UN TUBE ISOLE

Ce chapitre est consacré à une étude expérimentale de l'accumulation du dépôt des particules
sur la paroi d'un tube d'échangeur de chaleur. En effectuant des essais en conditions
encrassantes, on se propose d'examiner les différents mécanismes de dépôt sur la paroi d'un
tube lisse initialement propre placé perpendiculairement à un écoulement d'aérosols Le
transport des particules à la paroi va dépendre à la fois de la mécanique des fluides et de la
dynamique des aérosols.
Après avoir évoqué de façon générale le problème du transfert convectif sur un tube isolé,
nous présentons dans la deuxième partie le déroulement d'un essai caractéristique et décrivons
ensuite les résultats expérimentaux et leurs interprétations. Nous présentons enfin la
comparaison à la modélisation.

4.1 TRANSFERT CONVECTIF SUR UN TUBE ISOLE

Rappelons que lorsqu'un tube d'échangeur de chaleur est exposé à un écoulement d'aérosols,
une partie des particules va être capturée à la paroi du tube par différents mécanismes . la
diffusion. la thennophorèse, l'impaction inertielle, l'interception et la sédimentation Dans le cas
de la diffusion convective, la densité de flux massique de dépôt peut être exprimée par
(expression 2.7) :

<Pd = hm (Cpe - Cps)

hm représente le coefficient de transfert de masse (ou la vitesse de dépôt), Cpe et C ps sont
respectivement la concentration en particules dans l'écoulement principal et à la paroi.
Dans le cas des aérosols, il est classique de considérer que la diffusion de fines particules se
comporte comme la diffusion de grosses molécules de fluide. Ceci implique que le coefficient
de transport des particules est équivalent au coefficient de transfert de masse qui peut être
déterminé à partir des corrélations empiriques ou des équations théoriques pour le transfert de
chaleur.

78
 Chap IV Etude du dépôt sur un tube Isolé

Dans le cas d'un tube isolé dans un écoulement, par exemple, si on peut supposer que l'analogie
entre les transferts de chaleur et de masse est valable, on peut généralement exprimer le
coefficient global de transfert de masse par l'intermédiaire du nombre de Sherwood donné par

(4 1)

où C et a sont des constantes

4.1.1 EFFET DU BLOCAGE

Dans la pratique, les tubes d'échangeur sont généralement placés dans une conduite où l'effet
du blocage (paragraphe 3.2 2) peut être assez important. Le blocage affecte le profil de vitesse
du fluide autour de la paroi du cylindre, ce qui modifie le transfert à la paroi
Cependant, on constate que les études traitant ce problème sont peu nombreuses (nous n'en
avons répertorié qu'une seule [44] dans la littérature). Zukauskas et al. [44] ont étudié
expérimentalement l'effet du blocage sur le transfert thermique local autour de la paroi d'un
cylindre La figure 4. 1 montre leurs résultats.

KJ
1.1

U
I.J

1.1
Qg

OJ
as 0 JO 60 go 120 ISO ,.
position angulaire
Figure 4.1 Coefficient de transfert thermique local, pour différentes valeurs
du facteur de blocage, en régime subcritique [44].

On remarqu~ sur ces résultats que le blocage se traduit par l'augmentation du transfert de
chaleur sur la face amont du cylindre Selon ces auteurs, cette augmentation est principalement
due à la diminution d'épaisseur de la couche limite du fait de la croissance du gradient de
pression par accélération de l'écoulement. En ce qui concerne l'effet du blocage sur le transfert
thermique moyen, ces auteurs proposent les corrélations suivantes:

79
 Chap.JV Etude du dépôt sur un tube Isolé

Nu =0.408 (1+1.18k q 3) Reo so PrO.37 +0.058 (1+ 1. 18k q 3)l40 Re0 70 Pr0 37 (42)

qui est valable lorsque le facteur de blocage kq < 0.6, et :

Nu =0.408 (1+118k q 3) Reo so Pr0 37 +0.042 (1+1.18k q 3)140 Re0 70 Pr0 37 (43)

lorsque le facteur de blocage kq > 0.6

Selon ces auteurs, ces corrélations sont valables pour 100 < Re < 2.10 s , 07 < Pr < 80 et
0< kq < 0.9.
Dans notre cas, le tube est placé dans une conduite avec un facteur de blocage de 0.25
(paragraphe 3 2.2) Dans ce cas, si on suppose que l'analogie entre les transferts de masse et de
chaleur est valable, on peut exprimer le coefficient global de transfert de masse par
l'intermédiaire du nombre de Sherwood sous la forme:

Sh = 0.416 Reo so SC0 37 +0.060 Re0 70 SC0 37 (44)

ou par l'intermédiaire du nombre de Stanton pour le transfert de masse sous la forme

Stm= 0.416 Re-()so SC-()63 +0.060 Re-()30 SC-()63 (45)

Cependant, lorsqu'il existe un gradient de température à la paroi, les propriétés physiques du

gaz et de la particule changent et le transfert convectif se modifie Si le gradient de température
est suffisament important, la thermophorèse peut devenir le mécanisme de dépôt prédominant
et celle-ci est capable d'augmenter considérablement le transfert de masse à la paroi
(paragraphe 2.2.1.3).

4.1.2 EFFET DE LA RUGOSITE

Le paragraphe précédent évoque le cas du transfert convectif sur une paroi initialement propre.
En général, dès que la quantité déposée est importante la paroi devient rugueuse et le transfert
est modifié. On peut remarquer qu'au cours de la formation du dépôt de nombreux paramètres
évoluent (géométrie de la paroi, rugosité de la surface d'échange, section de passage). L'état de
surface peut varier en fonction du temps et en fonction de la position sur la paroi du tube.
L'état de rugosité de la surface peut avoir un effet important sur le transfert à la paroi La
croissance de la rugosité de la surface peut augmenter le transfert de masse d'un facteur qui
peut atteindre jusqu'à 10 [45]. Mais, il faut noter que la croissance de rugosité se manifeste
également par la croissance de la contrainte de cisaillement à la paroi, ce qui peut faciliter le
réentraînement.

Différents résultats expérimentaux issus de la littérature ont mis en évidence l'importance de la

rugosité de la surface. Les auteurs ([45], [46], [47], [48], [49] et [49a]) ont montré que la

80
 Chap IV Etude du dépdt sur un tube Isolé

rugosité de la surface augmente considérablement le transfert à la paroi, en particulier pour les

particules de taille moyenne inférieure à 5 J.lm ..
Selon Gudmundsson [47] et EI-Shobokshy et al. [48], l'accroissement du transfert à la paroi
est dû à la diminution d'épaisseur de la sous couche visqueuse et éventuellement à
l'augmentation de l'intensité de la turbulence de l'écoulement principal Davies [49] et [49a]
pense plutôt que la diminution de la distance d'arrêt de la particule et l'interception des
particules par les éléments rugueux sont responsables de ('augmentation de la vitesse de dépôt.
Par ailleurs, Shimada et al. [46] montrent que le transport à la paroi augmente avec la rugosité
de la surface, la taille des particules et l'intensité de la turbulence de l'écoulement incident

4.2 ESSAIS DE DEPOT SUR LE TUBE SONDE

Pour examiner le comportement du dépôt sur une paroi de géométrie simple, nous allons
effectuer des essais d'encrassement sur un tube isolé dans un écoulement d'aérosols On se
propose d'étudier l'influence de différents paramètres sur la cinétique de dépôt Le tube (tube
sonde) utilisé est un tube cylindrique hydrauliquement lisse de 20 mm de diamètre et de 40 mm
de longueur placé dans une section d'essai rectangulaire (figure 3.3, page 70).
La section d'essai est traversée par un flux d'aérosols de particules de Sulfate de Sodium
(Na 2S0 4 ) de granulométrie connue et bien contrôlée. Celui-ci est produit à partir des
gouttelettes d'une solution de sulfate de sodium générées par pulvérisation ultrasonore dans un
générateur de particules. En mesurant la masse de particules déposée sur la paroi du tube à des
temps différents, il est possible d'étudier la cinétique de dépôt pour les mêmes conditions
encrassantes.

4.2.1 DESCRIPTION DES ESSAIS

Les essais sont effectués en utilisant le dispositif expérimental décrit au paragraphe 3 2

Chaque essai consiste à encrasser le tube pendant une durée déterminée et pour des paramètres
(taille des particules, débit, température) définis. On utilise des particules de diamètre
aéraulique moyen de 2.20 J.lm et de 3.97 J.lm dont la distribution en masse est donnée sur les
figures 3 8 et 3.9.
Pendant l'essai, la température et le débit d'air à l'entrée de la section d'essai ainsi que la
température de l'eau de refroidissement du tube sont maintenus constants. Les essais sont
effectués à la pression atmosphérique. Chaque essai fournit une mesure de la masse du dépôt.

Les essais sont effectuées en suivant la procédure suivante :

- Avant d'effectuer un essai, le tube est monté dans la section d'essai en conditions propres.
Chaque remontage de la section d'essai est accompagné d'un contrôle de l'absence de fuites.

81
 Chap IV Etude du dépôt sur un tube Isolé

- Après une nuit de mise en température, le dispositif d'injection des particules est mis en route
- Le débit de gouttelettes est contrôlé en suivant en fonction du temps le niveau du réservoir
d'alimentation du pot de pulvérisation.
- A la fin de l'essai, le tube est démonté et le dépôt photographié Le dépôt est ensuite prélevé
soigneusement sur un papier taré à l'aide d'un pinceau et enfin on le pèse.
Pour une série d'essais, nous répétons les essais dans les mêmes conditions encrassantes en
faisant varier leur durée.

4.2.2 RESULTATS EXPERIMENTAUX

Nous avons réalisé pour différentes conditions plusieurs essais dont les résultats peuvent être
regroupés en trois séries. Les résultats expérimentaux sont présentés dans les paragraphes
suivants avec :
da : diamètre aéraulique moyen des particules (!lm)
Ue vitesse débitante à l'entrée de la section d'essai (mis)
Re nombre de Reynolds basé sur le diamètre extérieur du tube
Qa débit massique du gaz traversant la section d'essai (kg/h)
T.:e température moyenne du gaz à l'entrée de la section d'essai (OC)
Tsm : température moyenne du gaz dans la section d'essai (OC)
~ température moyenne de l'eau de refroidissement du tube eC)
tr durée d'essai réelle (h)
te durée d'essai corrigée (h)
m, : masse de particules injectées dans le circuit (g)
Cp: concentration moyenne en particules pour un essai (ppm)
Cpn : concentration moyenne en particules pour une série d'essais (ppm)
md : masse de dépôt par unité de surface (glml).

4.2.2.1 ESSAIS DE LA SERIE-A

Pour cette série d'essais, nous avons utilisé les particules de diamètre moyen de 2.20 !lm. Nous
effectuons les essais en maintenant constants le débit massique du gaz à 16 kg/h, la
température d'entrée du gaz à 90 oC, et la température d'eau de refroidissement du tube à
environ 40 oC. En faisant varier les durées d'essais, nous avons réalisé six essais dont les
résultats sont donnés sur le tableau 4.1.
Afin de vérifier la stabilité des conditions expérimentales, il faut observer la variation des
paramètres au cours du temps. La figure 4.2 montre la variation des paramètres mesurés lors
d'un essai réalisé en continu pendant 24.7 heures (essai AS). On peut remarquer que pendant

82
 Chap.IV Etude du dépôt sur un tube Isolé

l'essai le débit et la température d'air à l'entrée de la section d'essai, ainsi que la température
d'eau de refroidissement du tube sont pratiquement constants.

Tube· A : da:K 2.20 J,lm, Ue .. 2 mIs (Re .. 2360) C~ .. 124~m

essai Qa Tgm Ts t.r ml Cp t.c md
kg!h oC oC h g ppm h Q/m 2

Al 23.44 91.54 42.39 4.87 14.51 124.54 5.21 3.86

A2 23.47 87.47 42.29 7.53 22.96 126.10 8.25 5.37
Al 23.49 90.64 42.69 9.91 29.72 . 127.72 10.67 13.53
AA 23.44 88.24 42.74 15.38 43.66 119.58 15.68 JO.64
AS 23.84 88.14 42.99 24.67 67.82 115.32 24.36 54.11
A6 23.54 90.14 43.76 33.47 100.42 127.45 36.07 71.62
A7 23.54 90.27 43.70 64.68 133.24 87.50 47.85 91.28

Tableau 4.1 Résultats expérimentaux des essais de la série A.

100~--------------------------------------------------------.
90 .- .. _............ -_ ..•.. _-- .. - ........ _... -.- ... .
80 T Essai: AS
70 Tce (OC)
60
Qa (kg/h)
50
40 ..•• -..•. - .•••• - • • • - • - . -.. -1- • - - ••••. - .
30
20
10
.-- .-_._ .. _- --_ .. ---- ---- -.. -*-- ._-----
o +-----------+-----------+-----------+-----------+---------~
o 5 10 15 20 25
durée d'essai ( h )

Figure 4.2 Evolutions des paramètres expérimentaux.

Pour vérifier la reproductibilité des mesures de la masse déposée, nous avons répété l'essai Al
de courte durée et observé que les résultats obtenus en ce qui concerne la masse de dépôt
recueillie à la fin des essais ne différaient que de 5 % (La même démarche a été effectuée pour
l'essai B4 avec des résultats égaux à 3 % près). Nous sommes donc en mesure de penser que la
reproductibilité des essais est bonne.

A partir des résultats donnés par le tableau (4.1), on se propose d'étudier la cinétique de dépôt.

83
 Chap IV Etude du dépôt sur un tube Isolé

Pour prendre en compte de petites variations inévitables de la concentration en particules et

pour pouvoir comparer les différents résultats, nous avons choisi de corriger la durée réelle de
chaque essai En effet, si on considère que le coefficient de transfert est constant, en faisant
l'hypothèse de la captation parfaite, on peut exprimer la masse déposée par·
md =hm Cp t (46)

Pour avoir la même masse déposée, on définit une durée corrigée te

C
t =t _P- (47)
e r C
pm

tr représente la durée d'essai réelle, Cp la concentration réellement mesurée et C(nI la

concentration moyenne pour une série de n essais·
1 n
Cpm =- 2, Cp,n (4 8)
n 1=1

A partir des résultats du tableau 4 1, et avec une concentration moyenne C (ni de 124 ppm; on
peut à présent tracer l'évolution de la masse de dépôt en fonction te (figure 43)

140

- 120

--
1:
m
,Z
100

80
3e

lit

8. Tube - A
(;iJ
'CD
'0 da = 2.20 ~m
Tgm =90°C
8E 40 T5 = 43 oC
Ue = 2 mIs
20 Cpm = 124 ppm

0
0 5 10 15 20 25 30 35 40 45 50
t.e (h)

Figure 4.3 Evolution de la masse de dépôt pour des essais de la série A

La cinétique de croissance de la masse de dépôt se caractérise d'abord par une période à vitesse
de croissance sensiblement constante. On observe ensuite un décrochement correspondant à
une vitesse de croissance plus élevée qui ensuite s'atténue. Enfin, pour les durées les plus
élevées l'évolution semble devenir asymptotique.

84
 Chap.IV Etude du dépôt sur un tube Isolé

4.2.2.2 ESSAIS DE LA SERIE B

Pour la deuxième série d'essais, nous avons utilisé les particules de diamètre moyen de 2.20 /lm
et choisi de travailler avec un débit de gaz de 25 kglh, une température d'entrée du gaz de
ISO°C et une température de l'eau de refroidissement du tube de 40 oC Les conditions
expérimentales et les résultats de mesures sont donnés sur le tableau 4 2 Comme pour les
essais précédents, en observant la variation des paramètres expérimentaux, on peut constater
que la stabilité des conditions expérimentales est bOMe

Tube - 8 : da • 2.20 J,lm, Ue .2.8 mIs (Re. 2340), Cpm • 147 ppm
essai Qa Tgm Ts t.r ml Cp I.e md
kg/h OC OC h 9 ppm h g/m2

81 25.04 179.99 41.65 3.03 13.55 178.48 3.72 8.16

B2 25.04 181.61 38.93 7.35 30.80 167.45 8.45 24.07
' 83 25.01 181.37 38.81 14.68 46.55 126.67 12.77 35.21
B4 23.48' 182.52 38.06 26.52 65.77 105.70 18.04 63.66
as 25.03 182.39 38.95 24.80 99.72 160.63 27.35 81.13
B6 25.06 183.29 40.06 48.10 155.12 128.78 42.55 109.14
87 24.96 182.07 39.30 120.80 478.37 158.60 131.22 251.66

Tableau 4.2 Résultats expérimentaux des essais de la série B

400 ~-------------------------------------------------------------------------------------------------------------------------------------,

-'è 350

-Q
-
Gt
'Gt
300

250
lit 200 Tube - B
8.
'Gt
'0 150
da =2.20 IJm
Tgm = 182 oC
ftE 100 1er Ts = 40 oC
Ue = 2.8 mis
50 Cpm = 147 ppm
o ~~----------+----------------~---------------~---------------~---------------~-------------r---------------~
o 20 80 100 120 140
t.e (h)

Figure 4.4 Evolution de la masse de dépôt pour des essais de la série B

85
 <i= sens de l'éooulem.ent.
FÎg\l1'C' 4.5 Ph<ltograpl1îc du d~1 i la lin d'un essai de la série B .

L'évolution Cn fO.1'I.clion <Ill temps de 1 nW$Il déposée est donnéil sur la figurê 4.4. POUT ces
!isultats. nous avons corri.gé III durée d'I)$:W en IlÛsonmlIIt .5W" une (:~mcll'ntrntion moyenne en
panicule, de 147 l'pm De façon gé~rale, on retrouve la même umd:rnce en ce quiconceme
l'allu e de la courbe de révolution qu'avec les essaîs cie la série pnk6dente. La caurol': montre
que la cin1!liqué de dé;pôt se caractérise par trois résimes diffi!renlS~

ous avons ég/llem.eRl réalisé 'Pour quelques essai$ des photogfllphies du dépôt La pholo de la
figure 4.5 monlre le tube encr=é Ii la fin d'un essai d'une durée de 48 lIelites. On observe une
'variation importante de l'épai$SCur du dépôt, certainement li~ aux varia!ions du coèIDcilmt de
Itansfert sur la paroi du tube.
Au voisinage du polnt de stagnation amon!., on obse,",,'e que le dépôt est plus. impOl'\allt, Nou s
avons w sur la :figure 4.1 que les coefficients d'échange prennent une. valeur importante sur
cette zone. L'lIcvumulation du dépôt est donc: as5Q. lill'gement mvorisée. D'un aulre côté, on
remarque que sur la zoné dupainl de décollem.ent le dépôt est plus faible.

4.2.2. 2 ESSAIS DE LA SERIE C

Cette fois, nous avons effectué des es~is 1l'VeC les particules de diam~tre moyen de 3.97 IJ.m
Le débit du gaz trn"crsallt tll. section d'essai, la ttmpéralure d'enl~t;ie du gu et la (·empérarure

86
 Chap IV Etude du dépôt sur un tube Isolé

moyenne de l'eau de refroidissement du tube sont respectivement maintenus constants à 25

kg/h, ISO°C et 40 oc.
La figure 4.6 montre l'évolution de la masse déposée en fonction de la durée d'essai corrigée
Les résultats expérimentaux montrent la même tendance qu'avec les essais précédents en ce qui
concerne l'allure de la courbe expérimentale. On remarque cependant que la période de
transition est plus marquée et apparait plus tôt. Dans ce cas, l'explication peut venir du fait que
la rugosité de la surface encrassée est probablement plus forte pour le dépôt formé dans ces
conditions expérimentales

Tube - C : da • 3.97 ~m u. 1:1 2.8 mIs (Re = 2340) Cpm. 945 ppm
essai Qa Tgm Ts t.r ml Cp t.c md
kg!h OC OC h g ppm h g/m 2

Cl 25.04 182.70 40.31 5.17 112.38 868.7 4.89 91.28

C2 25.01 183.91 40.63 6.16 85.33 654.7 3.71 65.49
C3 25.02 184.00 41.19 9.35 157.08 671.6 6.83 161.10
C4 23.45 183.32 43.73 6.DO 173.24 1477.0 7.54 194.96
CS 23.50 - 183.40 43.25 13.68 344.37 1063.0 14.98 311.15
C6 24.99 183A9 39.81 17.10 ,".31 1044.7 19.42 329.25

Tableau 4.3 Résultats expérimentaux des essais de la série C

400

- 350 3e

-.:
~ 300
œ
- 250
Tube - C
8,200 da = 3.971Jm
'CD
'0 150 Tgm = 183 oC

~ 100
Ts = 42 oC
1er
Ue = 2.8 mIs
E Cpm=945ppm
50

0
0 2 4 6 8 10 12 14 16 18 20
t.c (h)

Figure 4.6 Evolution de la masse de dépôt des essais de la série C.

87
 ~ .sens de l'écoulement
Figure 4.7 POOlographJe du dépôt Il la fin d'un e.sIlai de la série C

La figure 4.1 montre la photographie de la surface encrassée ~ 1<1 fin (I.e "e;s$ru C6 de durée S
heures. On remarque tà aussi que l'.épai$$cur du dépôt wrie $UiVlll'l! la posi tion autoùr de la
par-oi. La masse déposée sur la face amant du tube e~t importarlle.

En oonclusion, les résultats expérimenlaux montrent que d'une l'lIiInÎère générale la répartition
de La masse déposée sur la paroi du tube corrcspolild IIUX variations des coefficients d'échange
autour de ce dernier. Les courbes montrant la dnélique de croi.ssance du dépôt Ibnt apparaitrc
Irois r~gimes différent Le premier Il. une vitesse de ç;ois:sance constante matS kible. Le
deuxième e~ carn,C1éri$é pli/" une vitesse de croissan.ce initiale plus. élevée,qui ensLlite s'atténLle..
Pendant le troisième régime. l'évolution semble p~ésenltr une CatBC,tÎ:të asymptoti.que.

4.3 INTERPRETATI ON o:es RESULTATS

Au coun; de l'accumulation du dépol. dt nombreux. p.anmel:res peuvent évoluer (modification

dllns le t.emps de l'état de la SI.J rface, des caractmsliques de ]'êcoulemeN, de la température dt
surfaoe). L'evolutioll dans te lemp$ de ce,s paramè-tres dépend f'orterMllt des conditions àIlx
Ilmitesel d.c$ ~acl';ri~tiq\les du dépôt en formatiot'l. n est done difficile de modéliœ.r pour

88
 Chap IV Etude du dépôt sur un tube Isolé

toute période de croissance la cinétique de dépôt à cause des modifications parfois très rapides
de ces paramètres
Dans le cadre d'une approche globale, on se propose dans les paragraphes suivants d'estimer la
masse des particules déposée pour les trois régimes de la cinétique de dépôt proposée

Dans un premier temps, nous estimerons pour les conditions correspondantes à nos essais la
masse déposée lors du premier régime de dépôt (phase à croissance constante) et comparerons
les résultats obtenus aux résultats expérimentaux. On s'intéresse au transfert de particules sur la
paroi d'un tube intialement propre et hydrauliquement lisse .
. En plus des hypothèses concernant les transferts de chaleur et de masse déjà formulées
auparavant (voir chapitre II), on introduit les hypothèses suivantes.
• Le gaz s'écoulant autour du tube est considéré comme un milieu continu (à l'échelle du
diamètre du tube) et également supposé gaz parfait. L'écoulement est supposé stationnaire
• Le tube est considéré comme un cylindre isolé dans un écoulement avec la convection forcée
dominant la convection naturelle
• On suppose que la turbulence de l'écoulement principal est faible. De même, la dissipation
visqueuse est négligeable.
• Pour les particules submicroniques, on suppose que l'analogie entre les transferts de chaleur
et de masse est valable
On considère que le dépôt est contrôlé par impaction inertielle, diffusion, thermophorèse,
interception et sédimentation.

4.3.1 DEPOT PAR IMPACTION INERTIELLE

Dans un premier temps, nous allons vérifier si les conditions d'un régime de dépôt contrôlé par
l'inertie des particules sont réunies Nous avons w au paragraphe 2.2 1.5 que le comportement
des particules dans un écoulement dévié par un obstacle dépendait de la valeur du critère de
Stokes, rapport des forces d'inertie et de traînée auxquelles sont soumises les particules.
Nous allons donc calculer pour les conditions correspondantes à nos essais le nombre de
Stokes effectif en utilisant l'expression (2.84)
Stk elf = 'l'(Re p ). Stk

Dans ce cas, le nombre de Stokes est donné par:

ppd/U_
Stk = (4.9)
9J.Ld c

89
 Chap IV Etude du dép6t sur un tube Isolé

où U .. est la vitesse initiale de la particule supposée égale à U e (vitesse moyenne du gaz à

l'entrée de la section d'essai), Pp et d p sont respectivement la masse volumique et le diamètre
de la particule, et de le diamètre du tube.

Nous retenons les données des figures 3.8 et 3 9 pour la distribution en masse de la taille des
particules considérées. Nous calculons ainsi pour chaque classe de taille de particules le
diamètre de la particule d p à l'aide de la relation suivante [12]

(4 10)

où Pp est la masse volumique de la particule, Pa la masse volumique de référence

(1000 kg / m3 ), et Ccun , Ca sont respectivement les facteurs de correction de Cunningham
(expression 2 67) pour les particules de diamètre d p et d.
En ce qui concerne Pp' nous retenons la valeur de 26799 kg/m3 pour la masse volumique
moyenne de particules de diamètre aéraulique moyen de 2 20 ",m et la valeur de 2625 5 kg/m 3
pour les particules de 3 97 flm (résultats de mesures effectuées grâce à un pycnomètre
automatique - Micromeritics Accupyc 1330).

Pour tenir compte la présence des particules non Stokesiennes (particules qui ne suivent pas la
loi de Stokes), on introduit le facteur de correction 'V(Re p ) que l'on calcule à l'aide de
l'expression (2 86)
3(J; Re p 1/3 - tan -1 ( J; Re p 1/3»
'V(Re p ) = 3/2
C Re p

avec c = 0 158 et Re p le nombre de Reynolds de la particules donné par

dpU_
Re : : - - (4 Il)
p V

où v la viscosité cinématique du gaz.

Enfin, nous effectuons le calcul en supposant que les propriétés therrnophysiques du gaz
peuvent être assimilées à celles de l'air sec (annexe-F) et calculons ces propriétés physiques à la
pression atmosphèrique et à la température de film : Tf::.!.. (Tgm + Ts )
2
où Tgm et Ts respectivement sont la température moyenne du gaz et celle de la surface
En supposant que la température moyenne de la surface du tube est constante et égale à la
température moyenne de l'eau de refroidissement du tube, on obtient les résultats donnés sur la
figure 4.8

90
 Chap.JV Etude du dépôt sur un tube Isolé

0.12

0.1

0.08

006
--=
:=
( .)
Q)

-
Q)
~
C I)

0.04

0.02

0
~ :2 ~ ~ ~ ~
lt') Il)
~
,.....
Il) Il)

d 0 d
Il) Ici <> cO 0:

diamètre aéraulique, da ( 1JI'n)

Figure 4.8 Nombre de Stokes pour des essais des séries A, B et C

Ces résultats montrent que les valeurs calculées sont inférieures à la valeur du nombre de
Stokes critique égal à 0.125 pour un cylindre isolé dans un écoulement potentiel.

De plus, dans notre cas le tube est placé dans une section d'essai avec un facteur de blocage
assez important (0.25). Comme noté au paragraphe 4.1.1, le blocage affecte le profil de
vitesses autour de la paroi du tube Konstandopoulos et al. [20] ont montré que le blocage va
augmenter la valeur du gradient de vitesse sur la face amont du tube (paragraphe 2 2 1 7). La
valeur du nombre de Stokes critique est donc diminuée.
Mais, comme le calcul relatif à notre cas n'est pas disponible, on peut souligner que le nombre
maximum de Stokes dans les conditions les plus défavorables (celles de la série C) est
suffisamment inférieur à 0.125
On peut donc en conclure que dans nos conditions expérimentales, l'impaction inertielle est
négligeable.

4.3.2 DEPOT PAR INTERCEPTION ET SEDIMENTATION

Nous allons à présent estimer la masse des particules déposées par les mécanismes
d'interception et de sédimentation.
Pour l'interception, nous calculons pour chaque classe de taille de particules l'efficacité de
collecte par interception à l'aide de l'expression (2.92) .

91
 Chap IV Etude du dépôt sur un tube Isolé

1
Tl",. = l+N",. - - - -
1+ Nin.
_ dp
où N "'. est le paramètre interception donnée par : N ml - -
de
Puis, nous estimons pour les N classes de taille de particules l'efficacité totale par '
N
Tlml = LTlml.1
1=1
Û)P.I (4 12)

Nous calculons ensuite la vitesse d'interception à l'aide de la relation suivante'

Tlml = Vml (AmI

U m AproJ
J= VIOl
U
.!
2
(4.13)
m

Aproj la surface apparente en amont du tube (= de le)

7t
Dans ce cas, on suppose que la surface nécessaire pour l'interception est égale à AmI =-de le
2
où 1e la longueur du tube

série da Qa Tgm TI Um Vint Vséd

(Jm) (kg/h) (oC) ~C) (mIs) (mIs) (mIs)

A 2.2 23,5 90 43 2.1 3.30E-tl4 2A1E-tl4

B 2.2 2S 182 40 2.8 4.39E-Q4 2.25E-Q4

C 3.97 2S 183 41 2.8 6.92E-Q4 5.32E-Q4

Tableau 4.4 Dépôt par interception et sédimentation

En ce qui concerne la sédimentation, nous estimons pour chaque classe de taille de particule la
vitesse de sédimentation à l'aide de l'expression (2.97) et nous calculons la vitesse moyenne de
sédimentation par
N
Vséd = ~ v
~ séd .1
û)
p.1 (4.14)
1=1

pour les conditions moyennes correspondantes à nos essais, le calcul des vitesses d'interception
et de sédimentation conduit aux résultats donnés dans le tableau 4.4

92
 Chap IV Etude du dépôt sur un tube Isolé

4.3.3 DEPOT PAR DIFFUSION ET THERMOPHORESE

On se propose dans ce paragraphe d'estimer la masse des particules déposées par diffusion sur
la paroi d'un tube isolé dans un écoulement d'aérosols On estime le coefficient de transfert de
masse en appliquant l'analogie entre les transferts de chaleur et de masse.
En supposant que chaque classe de taille de particules a un comportement indépendant, on
calcule le coefficient moyen de transfert de masse (vitesse de dépôt) pour les N classes de taille
de particules par :
N
hm =L hm,. (Op,. (4.15)
.=1
hm,! représente la vitesse de dépôt pour une classe de taille de particules donnée, et (Op,! la
fraction massique de la classe de taille de particules. Nous retenons pour le calcul, comme dans
le cas précédent, les données sur les figures 3 8 et 3 9 pour la distribution en masse de la taille
de particules.
Nous calculons hm,! par l'intérmediaire du nombre de Stanton pour le transfert de masse
h
St m,. =2!:!..
U (4 16)
m

avec U m la vitesse moyenne du gaz dans la section d'essai

En appliquant l'analogie entre les transferts de chaleur et de masse, nous déterminons St m,! à
l'aide de l'expression (4.5) :

Stm = 0 416 Re-{)so SC-{)63 +0.060 Re-{)30 SC-{)63

Dans ce cas, le nombre de Reynolds Re est donné par: Re = de U m

V
où V la viscosité cinématique du gaz évaluée à la température moyenne du fluide et à la
pression atmosphèrique.
Le nombre de Schmidt Sc est donné par: Sc =~
Db
où Db le coefficient de diffusion Brownienne que l'on calcule pour chaque classe de taille de
particule en utilisant la loi de Stokes-Einstein donnée par:

D = kaT (417)
b 31tlldp
23
où T est la température absolue et kB la constante de Boltzmann (1. 3 8 10- J / K)

Du fait du gradient de température, le transfert de masse à la paroi est augmenté par

thermophorèse Dans ce cas, la thermophorèse va augmenter la vitesse de diffusion par un
facteur Faug (paragraphe 2.2.1.3) qui est fonction de différents paramètres

93
 Chap IV Etude du dépôt sur un tube Isolé

Tgm -Ts
(aTLe, Pr, Sc, ).
Ts
En déterminant le nombre de Peclet pour le transfert de masse par l'expression (263), on
calcule Faug à l'aide de l'expression (2.55). Dans ce cas, le terme atLe (rapport entre la
diffusivité thermophorètique des particules et la diffusivité thermiqe du gaz porteur) est calculé
par l'expression (2 47) Ce terme est principalement fonction de la conductivité thermique du
gaz, celle de la particule et du nombre de Knudsen que l'on calcule par l'intermédiaire des
expressions (2 49) et (2 51). On estime la conductivité thermique de la particule en utilisant
l'expression suivante [11]:
070
À = 3 Ox 10-3
(
....!.-) cali s (4 18)
p 300 cmK

qui est valable pour les particules de silice (Si0 2 ) .

En admettant que la température de surface du tube est égale à la température moyenne de l'eau de
refroidissement du tube, nous avons calculé pour les conditions moyennes de nos essais le facteur
Faug , le nombre de Stanton Stm,t ainsi que la vitesse moyenne de diffusion (tableau 44) On
remarque que pour les conditions des séries B et C les valeurs du facteur Faug ont le même ordre
de grandeur, alors que pour la série A il est moins important. Ceci provient essentiellement des
valeurs de la différence de température entre le gaz et la paroi (TgJ1l - Ts ) .

série do Qo Tgm Ts Um Re Pr Sim Foug Slm,t hm

~ (kg/h) (OC) rc> (mIs) (mIs)

A 2.2 23.5 90 43 , 2.1 2161 0.7 8.38E-G7 326.3 2.73E-G4 1.19E-03

B 2.2 25 182 40 2.8 2319 0.7 7.15E-07 1217.8 8.71E-G4 5.52E-03

C 3.97 183 41 2.8 2280 0.7 1.04E-03 4.38E-03

25
'. 7.24E-G7
'"
1431.8

Tableau 4 .5 Dépôt par diffusion avec la thermophorèse

Quand on compare ces résultats aux résultats du tableau 4.4, on peut constater que dans nos
conditions expérimentales la diffusion avec thermophorèse est le mécanisme de dépôt
prédominant.

94
 Chap IV Etude du dépôt sur un tube Isolé

4.3.4 EFFET DE L'ISOLATION DU DEPOT

Après le premier régime, la courbe de la cinétique de dépôt est caractérisée par un

décrochement qui correspond à une augmentation de la vitesse de dépôt Puis, due à l'effet de
l'isolation du dépôt, la vitesse de dépôt s'atténue
On se propose dans ce paragraphe d'estimer l'évolution de la vitesse de dépôt en prenant en
compte de l'attenuation de la thermophorèse due à l'effet de l'isolation du dépôt
Afin de tenir compte l'effet de l'isolation du dépôt, J.M Grillot [50] a introduit la notion d'un
facteur d'atténuation thermophorètique Fa donné par .
_ Tgm -Td Ts
Fa - ~-----"-- (4.19)
Tgrn Tgrn - Ts

TglU représente la température du fluide et Td la température de surface de dépôt

Ainsi, en considérant le transfert de chaleur monodimensionnel sur la paroi du tube isolé par le
dépôt d'épaisseur uniforme et en admettant que le coefficient de transfert de chaleur h reste
constant, l'auteur a proposé la relation suivante .
1
Fa =---=---------
Tgm h
(420)
1+------md
Ts PdÂ,d

h
ou, en posant b =T.grns -,\-,
T Pdl\.d
on exprime Fa par ·

F = 1 (421)
a l+b md

Pd, Â,d sont respectivement la masse volumique et la conductivité thermique du dépôt.

Sur le modèle de base de l'accumulation du dépôt (expression 2.6, page 19), si on considére
que la densité de flux massique de dépôt est constante et celle de réentnlÏnement est
proportionnelle à la masse du dépôt, on peut exprimer :
drnd
--=a-Cm d (422)
dt
En introduissant le facteur Fa, l'expression devient:
dm a
- -d= -Crnd (423)
dt l+b m d

a, C sont des constantes

Lorsque la durée d'essai n'est pas trop longue, il est possible de considérer que le
réentraînement est négligeable (c'est à dire C = 0). Dans ce cas, par intégration de l'expression

95
 Chap IV Etude du dépôt sur un tube Isolé

(423) avec les conditions initiales t =0, md =0, on peut alors exprimer le modèle de la
cinétique de dépôt par :

(424)

Dans notre cas, comme les points expérimentaux du deuxième régime de dépôt sont décalés
(ils ne passent pas par l'origine), on peut exprimer la formule (4 24) en rajoutant un coefficient

(425)

(c x n'a pas de signification physique, mais il permet de tenir compte du décalage des points
dans le repère temporel)
Lorsque la durée d'essai devient plus importante, les effets du réentrainement se manifestent et
le coefficient C de l'équation (423) ne peut plus être négligé L'intégration de l'équation
différentielle donne

Alln lm-aIl +A21n Im-a21 =-bCt (4.26)

I-all l- a 21
La courbe représentative possède une asymptote d'équation m = a)
où a) est la racine positive de l'équation a - Cm - bCm 2 = 0
Toutefois, il n'a pas été possible d'ajuster les points expérimentaux pour obtenir des valeurs
significatives des paramètres b et C
Pour tenir compte du réentraînement, les points expérimentaux des deuxième et troisième
régimes ont été ajustés par
(427)

résultat de l'intégration de (4 22)

4.3.4.1 AUGMENTATION DE LA VITESSE DE DEPOT

Nous avons vu que la transition entre le premier régime de dépôt et le deuxième se caractérise
par une augmentation de la vitesse de dépôt due à la rugosité de la surface On peut exprimer
cette croissance par un facteur Fcr donné par :

(428)

Vd1 (t = 0) représente la vitesse de dépôt du premier régime à l'instant t = 0 et Vd2 (ttnJ) est la
vitesse de dépôt du deuxième régime à l'instant t = ttnJ En posant m d = 0 dans l'expression
(4.25), on obtient·

96
 Chap IV Etude du dépôt sur un tube Isolé

bc; -2c x
t m, =-.::.....-_... (429)
2a
La vitesse de dépôt Vd2 (t uu ) est proportionnelle à la pente de la tangente à la courbe donnée
par :

(430)

4.3.4.2 ESTIMATION DE LA TEMPERATURE DE SURFACE DU DEPOT

On se propose à présent d'estimer la température de surface du dépôt Td en déterminant à

partir des résultats expérimentaux le facteur d'atténuation Fa.
En effet, on peut définir pour chaque point expérimental un facteur d'atténuation expérimental
Faxp par le rapport de la pente de la tangente à la courbe cinétique pour t =te à la pente de la
tangente à la courbe pour t = tin,

F = II-bcxl (431)
a.xp ~I + 2abt c

On cherche ensuite la valeur de la température de surface de dépôt telle que le facteur

d'atténuation calculé Fa eal soit égal à Faxp :

F - Vd(T=Td)
(4.32)
a.eal - Vd (T = Ts )

Vd (T = Td ) est la vitesse de dépôt calculée en se basant sur la température de surface de dépôt

et Vd (T = Ts ) la vitesse de dépôt calculée pour la température de surface égale à la
température moyenne de la paroi du tube.
Le calcul de Vd se fait en utilisant la loi de dépôt sur surface lisse multipliée par le facteur
d'augmentation dû à la rugosité Fer et en prenant en compte l'augmentation de diamètre du
tube. L'estimation de la conductibilité thermique du dépôt est ensuite possible puisque

(4.33)

soit:

(434)

e représente l'épaisseur moyenne du dépôt et r le rayon du tube.

97
 Chap IV Etude du dépôt sur un tube Isolé

4.3.4.3 RESULTATS DE CALCUL

Série A.

Pour le premier régime de la cinétique de dépôt (figure 4.3), nous avons déterminé la pente de
la tangente à l'origine d:d 1t=O = 0 677 g 1m lh

Avec une concentration moyenne en particules de 0 122 g 1 m3 , nous avons calculé pour ces
conditions expérimentales une vitesse de dépôt de 1.547 10-3 mis
D'un autre côté, le calcul du modèle de dépôt par diffusion avec la thermophorèse (tableau 4 5)
conduit à une vitesse de dépôt de 1 189 1O-3 ml s (soit 0.77 fois de la valeur mesurée)

Si on considère maintenant le deuxième régime de dépôt, en ajustant aux points expérimentaux

l'expression (4.25) on obtient les résultats suivants a = 6.77 b = 0019 c x = -35 32
A partir de ces résultats, on peut calculer tuu (expression 4.29) ce qui donne t ml = 6 951 h._
et par l'expression (4 30) . dm d 1t = 4 067 g 1 m 2 h
dt ..
Enfin, en calculant le facteur d'augmentation Fa (expression 428), on peut constater que du
fait de la rugosité de la surface encrassée, la vitesse de dépôt est augmentée par un facteur égal
à 601.

Nous avons également calculé, en appliquant le principe de calcul décrit précédemment, la

température de surface et la conductibilité thermique du dépôt En supposant que dans ces
conditions expérimentales la masse volumique du dépôt est égale à 100 kg 1 m3 • on obtient les
résultats donnés dans le tableau 4 6

Tube - A : da • 2.20 um. U• • 2 ml' (Re • 2360). Cpm • 124 ppm

essai Qa Tgm T. t.e md Td lam.d
kg/h oC oC h g/m 2 oC W/mK
.
A2 23A7 87A7 42.29 8.26 5.37 45.13 0.210
Al 23A9 90.64 42.69 10.67 13.53 48A2 0.212
A4 23M 88.24 42.74 15.68 30.64 53.30 0.213
AS 23.84 88.14 42.99 24.36 5':11 58.51 0.216

lam.d : conducflbllHé thennlaue du dép6t

Tableau 4.6 Température et conductibilité thermique du dépôt (A)

98
 Chap IV Etude du dépôt sur un tube Isolé

Les coefficients AJ, A4 et A5 de la relation (4.27) proposée pour exprimer une cinétique de
dépôt tenant compte du réentrainement, valent respectivement 125.3, 0.03086 et 6 796 (figure
4.3)

Série B.

Comme dans le cas précédent, par ajustement des points expérimentaux du premier régime de
dépôt (figure 4 4), on obtient· d:d 1t=O = 2.742 g / m2 h

Avec une concentration moyenne en particules de 0.114 g / m J, on obtient une vitesse de dépôt
égale à 6 681 1O-J mIs; et l'application du modèle de dépôt par diffusion avec la
thermophorèse donne la valeur de 5.518 lü-Jm 1s (soit 0.83 fois de la valeur mesurée).

En ce qui concerne le deuxième régi~e, en ajustant les point expérimentaux par l'expression
(4.25), on obtient· a =15.55 b =0044 c x =-34.69
Si on calcule ttru par l'expression 4 29, on obtient t ml = 3.933 h, et en calculant la pente de la
tangente à ttnJ' on estime· dm d 1t = 6.155 g 1 m 2 h
dt ..
On peut donc constater que le facteur d'augmentation dû à la rugosité est égal à 2 25

De la même façon, en supposant que dans ces conditions expérimentales la masse volumique
du dépôt est égale à 100 kg 1 mJ, on calcule la température de surface et la conductibilité
thermique du dépôt et on obtient les résultats donnés dans le tableau 4 7

= =
Tube - B : da 2.20 um, Ue 2.8 mIs, Cpm = 147 ppm
essal Qg Tgm TI t.e md Td lam.d
kg/h oC oC h g/m2 oC W/mK

B2 25.04 181.61 38.93 8.45 24.07 73.19 0.068

B3 25.01 181.37 38.81 12.77 35.21 89.24 0.068
B4 23A8 182.52 38.06 18.04 63.66 101.34 0.069
as 25.Q3 182.39 38.95 27.36 81.13 114.17 0.070

lam.d : conductibilité fhennlque du dépôt

Tableau 4.7 Température et conductibilité thermique du dépôt (B)

Les coefficients AJ, A4 et A5 de la relation (4.27) proposée pour exprimer une cinétique de
dépôt tenant compte du réentrainement, valent respectivement 359.5, 0009167 et 0.4323
(figure 4.4)

99
 Chap IV Etude du dépôt sur un tube Isolé

Les coefficients A3 , A4 et As de la relation (4.27) proposée pour exprimer une cinétique de

dépôt tenant compte du réentrainement, valent respectivement 359 5, 0 009167 et 0 4323
(figure 4.4)

Série C.

Comme dans le cas précédent. si on détermine par ajustement des points expérimentaux la
pente de la tangente à l'origine (figure 4.6), on obtient: d:d 1l=O = 18 3 g / m 2 h Puis, avec la

concentration moyenne de 0 731 g / m3 , on obtient une vitesse de dépôt de 6.954 10-3 m / s.

Cependant, le calcul de dépôt par le modèle de diffusion avec la thermophorèse donne la valeur
de 4.378 1O-3 m / s, soit 063 fois la valeur expérimentale Cet écart, plus important que pour
les séries précédentes, est probablement dû au manque des points expérimentaux

Lorsque l'on considère le deuxième régime de dépôt, en ajustant les points expérimentaux 2, 3,
4 (figure 4.6), on obtient une évolution linéaire On n'a pas pu observer l'effet de l'atténuation.
Ceci est également dû au manque des points expérimentaux, car les points 5 et 6 appartiennent
au troisième régime de dépôt
Les coefficients A3 , A4 et As de la relation (4 27) proposée pour exprimer une cinétique de
dépôt tenant compte du réentrainement, valent respectivement 3692, 0 139 et 2 532 (figure
4.6)

4.4 CONCLUSION

Nous avons étudié expérimentalement l'accumulation de particules sur la paroi d'un tube
d'échangeur de chaleur placé perpendiculairement à un écoulement d'aérosols Il s'agit du dépôt
de particules de sulfate de sodium sur la paroi d'un tube cylindrique en acier inoxydable
initialement propre.
Les résultats expérimentaux montrent que d'une manière générale la cinétique de croissance du
dépôt est caractérisée par trois régimes différents Au cours du premier, la vitesse de
croissance est pratiquement constante mais faible Le deuxième régime commence avec une
vitesse de croissance plus élevée qui ensuite s'atténue. Pendant le troisième, l'évolution semble
présenter une caractère asymptotique

Pour les essais des séries A et B, effectués avec les mêmes particules malS une vitesse
légèrement différente et surtout des conditions thermiques différentes, l'interprétation proposée
rend bien compte de la modification du transfert de masse due à l'augmentation de la rugosité

100
 Chop IV Etude du dépôt sur un tube Isolé

de la surface mais aussi de l'effet de l'atténuation de la thennophorèse par l'isolation thermique

crèée par le dépôt
En faisant l'hypothèse d'un dépôt unifonne de masse volumique connue, il a été possible de
calculer pour chaque point expérimental du deuxième régime de dépôt la température de
surface et la conductibilité thennique du dépôt. Les valeurs obtenues pour la conductibilité
thermique sont réalistes

Lorsque la durée des essais est plus importante, le phénomène de réentraînement apparait et
devient prepondérant Au cours de l'essai C réalisé avec des particules de taille supérieure et
une concentration massique de l'aérosol plus importante, le manque d'essais de courte durée
n'a pas permis de mettre en évidence l'effet d'atténuation de la thennophorèse

101
 Chap.5 Etude du dépôt sur unfOlsceau de tubes OI/etés

CHAPITRE V

ETUDE DU DEPOT SUR UN FAISCEAU DE TUBES AILETES

Le présent chapitre est consacré à l'étude expérimentale du comportement d'une maquette

d'échangeur de chaleur à tubes ailetés vis à vis de l'encrassement particulaire en phase gazeuse
Elle fait suite à l'étude menée par J.M Grillot [1] sur l'encrassement d'un échangeur air/air à
plaques corruguées (même nature de particules) Nous proposons dans cette étude d'analyser
les effets des différents paramètres sur le comportement de l'échangeur en conditions
encrassantes
Après avoir procédé à la description générale de l'écoulement dans un faisceau de tubes à
ailettes, nous présentons le déroulement d'un essai caractéristique et décrivons les résultats
expérimentaux obtenus et leurs interprétations Nous présentons enfin la comparaison des
résultats expérimentaux à la modélisation.

5.1 ECOULEMENT DANS UN FAISCEAU DE TUBES

Avant d'évoquer l'écoulement sur les tubes munis d'ailettes, nous présentons l'organisation
générale de l'écoulement dans un faisceau de tubes lisses
D'une manière générale, l'écoulement à travers un faisceau de tubes est complexe et dépend de
la géométrie du faisceau, de la position du tube considéré et des conditions
thermohydrauliques Les variations locales du coefficient de transfert de chaleur et de masse
sont importantes et dépendent fortement de la position des points de décollement et des points
d'impact de l'écoulement.

Dans un faisceau de tubes lisses (figure 5 1), on observe un écoulement secondaire dans
l'espace intertube qui se sépare de l'écoulement principal On peut également distinguer un
point de décollement et un point d'impact de l'écoulement sur le deuxième tube. Selon
Zukauskas et Kami [51], l'emplacement de ces points varie avec la géométrie et les conditions
de l'écoulement. Ces points correspondent à des
- variations brusques de pression à la paroi du tube,
- inversions du sens des vitesses pariétales,

102
 Chap 5 Etude du dépôt sur un faIsceau de tubes a/Jetés

- maxima du transfert thermique ou massique local pour les points d'impact de l'écoulement,
- minima du transfert thermique ou massique local pour les points de décollement

Figure 5 1 Organisation générale de l'écoulement dans un faisceau aligné [51]

'c-...oc. .

Figure 5 2 Modèle de l'écoulement sur une ailette [52], [53]

103
 Chap 5 Etude du dépôt sur un JO/sceau de tubes a/Jetés

Dans le faisceau de tubes à ailettes, l'écoulement est plus complexe, très varié et fonction de
différents paramètres : thermique, hydraulique et géométrique Les conditions de l'écoulement
au niveau de la premiere rangée de tubes sont différentes de celles régnant autour des rangées
suivantes De façon générale, l'écoulement est caractérisé par l'effet de séparation de la couche
limite et par l'effet des intéractions entre les sillages La premiere rangée de tubes se comporte
comme une grille de turbulence permettant d'augmenter le transfert pour les rangées suivantes

UQG _ ____

Figure 5.3 Tourbillon en fer à cheval [52], [54], [55]

Owen [52] et Ireland et al. [53] ont montré que l'on peut séparer l'écoulement en plusieurs
zones (figure 5.2) :
a Zone de détachement de l'écoulement après contact au sommet de l'ailette
b Zone de recirculation de l'écoulement sur la partie amont de l'ailette, causée par
les tourbillons spirales.
c Zone de formation d'un tourbillon en forme de fer à cheval (horse shoe-vortex)
Ce tourbillon crée une paire de tourbillons se dévellopant à l'aval du cylindre
d. Deux points de séparation de la couche limite sur la paroi de tube.
e Zone de recirculation de l'écoulement à l'arrière du tube
On présente sur la figure 5 3 la forme et l'influence du tourbillon en fer à cheval Ce tourbillon
se crée à cause du blocage occasionné par le tube dans le canal inter-ailettes, s'enroule autour
du tube, et finalement se propage dans la suite du faisceau. Ce tourbillon est responsable des
mélanges énergiques intervenant entre le fluide près de la paroi et le fluide du reste du canal
inter-ailette; il favorise ainsi les échanges de chaleur.

104
5.1.1 APPROCHE :NUMERIQUE DE L'ECOULEMENT

Nous ~tons d;ans té pmarapho les rtsultats d'lin calWl num.~riquc décrivant les
canet&btiques locales de 1'6cou1cmc:nt dans notre maquette d'~geur t!lcrmiquc (voir
figure 3.4, page 71) réalisé par P. Mercier el al. (S'J. Ces auteurs Oltt uillisé le logiciel TRIO-
VP d.éveIoppé par le CEA qui pemII!t de résoudre \es ~&tions de N&vier-Stokes couplées
avec Nquation de l'é:nctgie et un modèle de tutbutcnce k-e en utilisant la m6thode des volumes
finis.

~ . -- ....
, '''''': ..... ~......... - -"-'40"'.1.;'11 1 t

VICTIliU!l V:I .QI

.~~~ ~ -- --- - ~ ~---

Figure SA ProOI de vitesse d~ le domaine du fluide pour uo fa~ de tubes sans all~L1C

lOS
 ,
.. t • .. a

-- ~

-- -

Figu re S. S Profil de vitesse dans le domaine du fluide :pour un i.~~u de tub~ 8\'eC ailette

La figure SA manlI@ le profil de viteue d~ le domaine du l'Iù~de pour un ftlsœau. cie It+bes
~ ailettes.. Sur ces r&lltats. 0.11 rtn1&1'tfue .i~méllt une zone de recirculation â j'&rrière du
tube et Un poInt d'limpact devtcoulement sur la face amont du d.eu)(jème mbe.
Sur la figure S.S. on manlTe l~ profil de vitesse dms le dor.illc du fluide pour un f'aisœau de
tubes: ave<: ~Iettes. On mnarque qu le pMnomène observé ptécédemenl se répèle dllls la
panic ai1.~ On observe la fotm;lihOn des siUa,ges chins la. zone l!Oituêe derrière ailette. On
retrouve également le poinl d'impact de l'éooulementsur·les llileltc

5.1.2 DISTRIBUTION DE TEMPERATURE

Nws présentons i présent les résulllll$ de mesure, de Il distributÎon de lemperalure sur la

surfate d'1Û!e:ttc réal~ par Le 00 tU DI, [57]. Les aul,eun. ont. rail des mesures locales de
tèl:TlJ)érature: sur la sumoe d'aileUe plane cin:tllaifle en utilisant la technique de thermographie

106
infruouge. Un faiseeau de tubes ,est exposé ,6 un koulem.cnt d'air à la température ambiante el
Il. III pression umospMrique. Le tt/be de mesure cMpllçablc dans la scçtion d'e5AI est chautré
par une cltculadoo d'eau chaude. La r6Wlall de mesure so,,! ~rêse!ltês sous fol'm.c
adimensiOMelle, àllaidc du pamnètre ydi!fini comme suit:

T.. repre:serue 1. température de l'ailette. T_ l:a tempérlture de l'air IIJJlbiMt ci 1- !il

lem~ure de l'eau chaude.

LI
..
'
1

•
•
•

•
,

, •
, •
•

.Figuree 5.6 Dittributr.on de ternpèrarur1i sur des ailettes [57]

La figure 5.6 montre les rul,!Jtm, demewl'e pouf les ailettes disposées sur Ics tubes au 1er
rong. Sême ran.s et 9ème rang en. fàisceau de tlIbes a~gnés. Les essais ont ét.é effectués pour un
nombre de Reynolds de 6800 environ (oombfe de Reynolds basé sur le diamètre du tube), Sur
ces résultats, on rem~l;Ie clairement qJ,1e 1. distribution de t.empérature sur une ailette n'est pas
urûfOl'l'l'lC.. Da!lS celte figure. les zones froides (en bleu) c.ormpon4ent à on bon IIGhan,gè de
dhaltllr et les toIleS chaudes (en rouge),6 lin mauvais édianse de chaleur.
Enfin, lesauleurs onl soilligné les pointllllllivantll ;

1.07
 Chap 5 Etude du dépôt sur un faisceau de tubes alletés

- La région à fort échange de chaleur en amont de l'ailette est le résultat de la contraction de

la section de passage du fluide, et est due au développement du tourbillon en "fer à cheval" qui
favorise un fort mélange dans le fluide.
- La région à faible échange à l'arrière de l'ailette est due à la faible vitesse d'écoulement.

5.2 FAISCEAU DE TUBES AILETES EN CONDITIONS ENCRASSANTES

Nous proposons dans les paragraphes suivants d'étudier les effets du dépôt de particules sur la
performance d'un échangeur thermique. Dans ce but, nous allons effectuer des essais en
conditions encrassantes sur l'échangeur et examiner l'influence de différents paramètres sur le
comportement de notre maquette expérimentale.

Deux types d'essais seront effectués : des essais en régime de transfert de chaleur et en
isotherme. En régime de transfert de chaleur nous allons nous intéresser à l'évolution en
fonction du temps de la résistance d'encrassement et de la masse de dépôt. En isotherme, nous
ne nous intéressons qu'à l'évolution de la masse de dépôt. Les caractéristiques de l'échangeur
utilisé sont données sur la figure 3.4 (page 71).
Chaque campagne d'essais consiste à encrasser le faisceau de tubes pendant une durée
déterminée et pour des paramètres (débit, température) définis. Les essais sont effectués à la
pression atmosphèrique. En ce qui concerne les aérosols, nous utilisons comme dans le cas de
dépôt sur le tube sonde les particules de sulfate de sodium de diamètre aéraulique moyen de
2.20 J.lm et de 3 97 J.lm dont la distribution en taille est donnée sur les figures 3.8 et 3.9 (page
76).

Les campagnes d'essais sont effectuées en suivant la même procédure que celle du cas de dépôt
sur le tube sonde (paragraphe 4.2.1). Pendant la campagne d'essais, toutes les grandeurs
(températures et débits) sont suivies automatiquement et pour le suivi de l'encrassement on
acquiert 10 séries de mesures toutes les dix minutes et les moyennes sont stockées pour les
calculs ultérieurs.
A la fin de chaque campagne d'essais, après démontage de la section d'essai le dépôt est
photographié et soigneusement déposé sur un papier taré. Puis, on pèse séparément le dépôt
sur les demi-tubes supérieurs et inférieurs de l'enveloppe, ainsi que le dépôt sur les tubes
centraux.
De plus, avant de l'utiliser pour les campagnes d'essais en conditions encrassantes, nous avons
effectué une série d'essai en conditions propres pour déterminer la loi d'échange thermique de
cet échangeur.

108
 Chap.5 Etude du dépôt sur un faisceau de tubes aI/etés

5.2.1 PERFORMANCES THERMIQUES EN CONDITIONS PROPRES

Nous avons réalisé une série d'essais avec de "l'air propre" pour déterminer la loi d'échange.
Les mesures ont été effectuées en faisant varier le débit massique du gaz traversant la section
d'essai pour des valeurs comprises entre 6 kglh et 41 kg/h ce qui correspond à 1100 < Re <
7200 (nombre de Reynolds basé sur le diamètre du tube). Dans ce cas, la température d'entrée
du gaz ainsi que la température de l'eau de refroidissement des tubes sont maintenues
constantes à 185 oC et 40 oC respectivement.
A partir des valeurs enregistrées et en calculant le bilan thermique dans l'échangeur étudié (voir
annexe A), nous déterminons le coefficient d'échange coté gaz pour le faisceau de tubes dont
les résultats sont représentés sur la figure 5.7.
\
-'
70
2 60
,
~~
~50 ~
~
~40 ~
~
-
.: 30
<0
()2Q ~
..-

z=- 10 I~
o
o Hm 2CXD 3(0) 7CXD 8CXD
Nombre de Reynolds, Re

Figure 5.7 Coefficient de transfert thermique côté gaz - faisceau de tubes ailetés

Par ajustement aux points expérimentaux, on représente ces résultats sous la forme
Nu =0.024 Reo 896 PrO 333 (5.1)

Nous retenons ce résultat pour la suite de l'étude.

5.3 ESSAIS DE DEPOT EN TRANSFERT DE CHALEUR

Nous avons réalisé quatre séries d'essais en régime de transfert de chaleur Nous avons choisi
de travailler à la température constante de 185 oC pour le gaz à l'entrée de la section d'essai et
à la température constante de 40 oC pour l'eau de refroidissement des tubes dans le faisceau.
En évaluant le bilan thermique dans la section d'essai, il est possible d'étudier l'évolution de la
résistance d'encrassement pour différentes conditions expérimentales.

109
 Chap 5 Etude du dép6t sur un fOisceau de tubes Ol/etés

Pendant les campagnes d'essais, nous avons observé que le débit massique de matière sèche
injectée dépend fortement de la concentration en sel de la solution pulvérisée puisque le débit
de gouttelettes reste sensiblement constant. C'est pourquoi nous avons choisi de représenter
l'évolution des paramètres non pas en fonction du temps mais en fonction de la masse Mi de
particules injectées dans la section d'essai.
Soit :
Ml =Qg.Cl·t (5 .2)
avec Qg débit de gouttelettes injectées,
CI concentration en sel de la solution pulvérisée, et
t temps écoulé.
Pour la comparaison des essais avec des particules de même taille cette présentation permet
d'éliminer l'influence de petite variation de Qg.
Les résultats expérimentaux de chaque série d'essais sont représentés sur un tableau précisant :
da : diamètre aéraulique moyen (Ilm)
Qa : débit massique de gaz chargé en particules (kglh)
Tce : température moyenne de gaz à l'entrée de la section d'essai (OC)
Tgm . température moyenne de gaz dans la section d'essai eC)
Ts . température moyenne de l'eau de refroidissement des tubes eC)
Ue . vitesse débitante moyenne à l'entrée de la section d'essai (mis)
Um : vitesse débitante moyenne dans le faisceau (mis)
dur durée d'essai (h)
md masse de dépôt total recueillie (g)
Cp concentration moyenne en particules (ppm)
Re nombre de Reynolds donné par : Re = dt U m
V
où dt le diamètre extérieure du tube et v est évalué à la température de film

Résistance d'encrassement.

Nous calculons la résistance d'encrassement à partir des données de mesure de température et

de débit en utilisant la relation suivante
R __1__ 1
(5.3)
d - H(t) H(t = 0)

H(t) représente le coefficient global d'échange à l'instant t et H(t=O) celui à l'instant initial.
Le coefficient d'échange global H ramené à la surface externe des tubes ailetés est calculé par:

H= Qc (5.4)
AtDTLM

110
 Chap.5 Etude du dépôt sur unla/sceau de tubes a//etés

où Qe représente la puissance calorifique échangée, AI la surface d'échange extérieure des

tubes, et DTLM l'écart logarithmique moyen de température donné par .
DTLM= (Tes -Tfs)-(Tee -Tfe ) (55)
ln (Tes - Tf.)
(Tee - Tfe )

avec Tee et Tes respectivement la température d'air à l'entrée et à la sortie de l'échangeur, Tfe et
Tfs la température de l'eau à l'entrée et à la sortie de l'échangeur
Nous utilisons pour le calcul de H la puissance calorifique Qe cédée par l'air
(5 6)

où me représente le débit massique du gaz, Cp la chaleur massique à la pression constante

(calculée à température moyenne), et .1Tc l'écart de température du gaz aux bornes de la
section d'essai. En effet, la mesure du débit massique du gaz est plus précise (débitmètre
massique thermique avec régulation) que le débit de l'eau de refoidissement des tubes
(rotamètre) et de plus la différence de température côté gaz est plus forte, donc plus précise.
En ce qui concerne l'incertitude expérimentale, on peut estimer l'incertitude relative sur la
valeur de la résistance d'encrassement en fonction des différentes variables (voir annexe-C).
A titre d'exemple, pour un essai, pour une incertitude relative de 1% sur le débit d'air et une
incertitude de 1°C sur la température, on peut obtenir une incertitude relative de 5.1 % sur la
valeur de Rd

5.3.1 RESULTATS EXPERIMENTAUX

SERIE - A.

Pour cette première série d'essais, nous avons utilisé les particules de diamètre aéraulique
moyen de 2.20 !lm et choisi un débit du gaz traversant la section d'essai égal à 25 kg/h Nous
avons réalisé six essais dont les résultats figurent dans le tableau en annexe E .
La figure 5.8 montre l'évolution des paramètres mesurés lors d'une campagne d'essais réalisé en
continu pendant 156 heures Sur ces résultats, on observe une bonne stabilité des conditions
expérimentales

L'évolution de la résistance d'encrassement est donnée sur la figure 5.9. On peut remarquer que
la reproductibilité des valeurs expérimentales est bonne Nous notons cependant qu'aucun essai
n'a été assez long pour faire apparaitre une valeur limite de l'encrassement On ne peut pas
mettre en évidence un palier pour lequel le réentraÎnement tendrait à égaler le dépôt, bien que
la masse de particules injectées dépasse 516 g.

111
 Chap.5 Etude du dépôt sur un/a/sceau de tubes a/Jetés

2oor---------------------------------------------------~
1~~------------------~----~--------~------------
160
140
t
Tee (OC)

120
100

~ _~.-_~.~~;~_~~~~.=~~~~.-~.~~~~~~;.;r;.~:=~~~~~,~~~~.-~
O+-----~-------+------~------r-----~------~------+-----~

o 20 40 60 100 120 140 160

durée d'essai (h)

Figure 5 8 Evolutions des paramètres expérimentaux

0.04 -r----------------------------------------------------------,
Rf = 0.0299 ( 1 - exp ( -0.0032 mi ) )
0.035

- 0.03+-----------------------------------~--------

~ 0.025
~
~ 0.02 série - A:
;; 0.015 da =2.20 IJm
0.01 Tgm = 139 oC
Ts =40 oC
0.005 Um =3.1 mIs
O~----_+------+_----_+------r_----~------r_----~----~

o 100 200 300 400 500 600 700 800

masse Injectée ( g )

Figure 5 9 Evolutions de la résistance d'encrassement

Nous avons également réalisé pour quelques essais la photographie de dépôt à la fin des
campagnes d'essais La figure 5. 10 montre le résultat de la photographie réalisée à la fin d'un
essai de durée 218.2 heures. On observe que le dépôt sur la premiere rangée de tubes est plus
important que le dépôt sur les rangées suivantes Le dépôt se localise préférentiellement dans la
zone des points de stagnation amont des tubes et des ailettes. A l'arrière des tubes, on observe
que le dépôt est plus faible. Il semble que ceci soit dû à l'écoulement secondaire de faible
vitesse.

112
 <i= sens de l'écoulement

Flgure 5.10 POOtopphie du dépôt pour lin essais de la série A

( da .. 2.20 J11R, QI. - 14 kgln, Tgm - I4l"C,. Ts; . <II oc. dur - 2 18 b)

A l'intérieur du faisceau. la localisation du dépôt est saRS doule .Ia CQns6quenoe de-
l'orpniJation gtnérale de l'écoulement dans le msceau delube$ en arrIl\gement allgn6. DaI.'\$
OC cas, lM Pitt culos se déposent préférenlieUemenl am la zone de$ po[nts d~mp:act des Wbc$ et
dalilettes.

r!g\lte S. II Schématisation du dépôt clans lin. faisceau de ruiles [58 ]

Quand on compare œs résuTlllls.awc le.!! condilions loea.les de r«Oulemcnl. dw un filisceau dé

tubes (~ SA el S.S), (ln. m!lAffl\le que la position maxlm.1Im ct minimum ~r de
dépôt correspond à. t. position. R\l)Cimllm ft minimum des coeffi.eil!lJts d'échange !!Ur la swfjn:e
des tubes aileli$ dans le rusceau. Ceci rend bim compte de III l'onne gi!n~e d:e d6p6t
sehêmatisée sur III. f®lre S. Il .
 Chap.5 Etude du dépdt SUI' unfaisceau de tubes alletés

A la fin des campagnes d'essais, nous avons pesé séparément le dépôt recueillie sur les demis
tubes supérieurs et inférieurs de l'enveloppe ainsi que le dépôt sur les tubes centraux.
Sur les résultats de mesure (voir annexe-E), on remarque que la masse déposée sur les demis
tubes inférieurs de l'enveloppe est généralement plus importante que celle sur les demis tubes
supérieurs. Ceci est probablement dû à l'effet de sédimentation directe (pour les particules de
plus grande taille) et indirecte (lors du réentraînement, des agglomérats de particules sont
arrachés et peuvent chuter) mais également une chute du dépôt lors du démontage de la
section d'essai (bien que celui-ci soit soigné, on ne peut éviter quelques chocs mécaniques).
Cependant, nous allons dans ce qui suit nous intéresser à la masse déposée sur la totalité de la
surface des tubes du faisceau

25
Tgm = 139 oc
-20 Ts =40 oc
-
Q

G»
'~ 15
da = 2.20 IJm
Qo = 14 kg/h
0
~
'CD (série-A)
"la
G»
la
"E 5

0
a 100 200 300 400 500 600 700 800
masse Injectée ( 9 )

Figure 5.12 Evolution de la masse de dépôt des essais de la série A

La figure 5. 12 montre l'évolution de la masse de dépôt recueillie à la fin des essais. On

remarque que la masse déposée augmente pratiquement linéairement avec la masse des
particules injectées. En déterminant la pente de cette droite, on peut estimer l'ordre de
grandeur du coefficient de transfert de masse (ou la vitesse de dépôt).
Dans ce cas, avec une concentration moyenne de 0.125 g/m3 et un débit massique des
particules de 2 3 g/h, on estime que la vitesse de dépôt pour ces conditions expérimentales est
de l'ordre de 0.8 10-3 mis.

114
 Chap.5 Etude du dépôt su,. unfOlsceau de tubes al/elés

SERIE· B.

Pour la deuxième série d'essais, nous avons utilisé les particules de même taille, mais nous
augmentons le débit à 25 kg/h Nous avons également vérifié la stabilité des paramètres
expérimentaux au cours du temps et constaté que celle-ci est pratiquement bonne

0.02 Rf = 0.0174 ( 1 - exp (-0.0039 mi)

~ 0.015 r------~=::;;:~;;;;;Z~~===
~ série - B:
li =2.20 ~m
- 0.01

0.005
da
Tgm = 145°C
Ts =42 oC
Um = 5.6 mis

O~---r----~--~--~----~---+----+----+----~--~
o 100 200 300 400 500 600 700 800 900 1(0)

masse Injectée ( g )

Figure 5 13 Evolution de la résistance d'encrassement (série B)

0.03

0.025

-
~
~
0.02
série B
~ 0.015 Qa = 25 kg/h

-
Q2 0.01

0.005

0
0 100 200 300 400 500 600 700 800 900
masse Injectée ( g )

Figure 5.14 Evolution de la résistance d'encrassement (série A et B)

115
 ~ sens de réooulement

Figure S.l 5 Pholo8fàphle du dépôt pour un essais de la série B

( da = 2.20 f1m, Qa = 2S kgIb, TSin = 14.9 ·C, Ts = 39 ·C, dur == 166 h )

ro r------------------------------------------------,

O~~~----~----~--~----+_--~----~----~--~
o 100 300 600 700

Figure 5. 16 Evolution de III muse de d.èpôt des ~is A et B

Sur la 6gurc 5. 14, on relTlarque que la rési~tmoe d'encfMSemenl est RlOillS. importante lorsque
le débit ~lptus irnpol'tant. Ceci est SlIns doute lié li l'effet de l'isolallon du depôt (6gure 5.1S).
Cependant. comme pour les essais de la ,séne précédente. 011 ob erve Qu'aucun essai n'a été
assl!Z long pour Ilùre appaiïÙtre Llne valeur limite de !'ellCl1I8Sement bie:n que la. rru!$SC de
parti.culés injectées déplisse 700 g.
En ce qui COnCerne la répaniliDntfu dépôt dllllS le faisceau (figure 5, 15), on observe ]a mSme
tendance que cette des essais de Ill. série prêcëdente. On remarq"e ct:p'enrlant une légère
<;rOimnoe enee qui concerne le dépôt au niveau des impacts des tubes ct dt aileltcs.

116
 Chap.5 Etude du dépôt sur unfDlsceau de tubes DI/etés

La figure 5.16 représente l'évolution de la masse de dépôt. On remarque que par rapport aux
essais précédents la vitesse de dépôt est plus importante. En déterminant la pente de la droite
et en utilisant les valeurs de 0,101 g/m3 et de 3,28 g/h respectivement pour la concentration
moyenne en particules et pour le débit massique des particules injectées, on peut estimer que la
vitesse de dépôt est de l'ordre de 2.22 10-3 mIs.

SERIE· C.

Pour cette série d'essais, nous avons utilisé les particules de diamètre 3.97 !lm Sur la figure
5 18, on remarque que plus les particules sont grandes moins la résistance d'encrassement est
importante L'explication peut venir de l'effet de l'isolation du dépôt La figure 5 19 montre que
pour ces conditions expérimentales les particules se déposent préfèrentiellement au niveau des
impacts des tubes et des ailettes On peut penser que l'effet de l'isolation par le dépôt est moins
important.

L'évolution de la masse déposée est donnée sur la figure 5 20 Pour ces conditions
(concentration moyenne en particules 0.746 glm3 et débit massique des particules 24 36 g/h),
on estime que la vitesse de dépôt est de l'ordre de 3.41 10-3 m / s.

O.OOr-------------------------------------------~
0.018 Rf =0.0153 ( 1 - exp (-0.0028 mi) )

0.016 t-------------------+-----
- 0.014
~ 0.012 série - C:
~ 0.01 da = 3.971Jm
- 0.008 Tgm =144 oC
i52 0.006 Ts = 40 oC
0.004 Um =5.6 mIs
0.002
O~-~----~---+----+----~---~--_+----~--~
o 100 200 300 400 500 600 700 800 900
masse Injectée ( g )

Figure 5.17 Evolution de la résistance d'encrassement (série C)

117
 Q.02
0;018 ~étle a
0!016 do .. 2.20 J.ITI
i' 0,0 14
~. 11012
" om cio .. 3-9i !JIll
! o.oos
<2 Q.(X).<,
sérleC

0,.00:1
UOO2
o
o 100 300 400 500 (IX) 100 800
mosse tnlKtée' g )

Figure 5. 1S Evolullon de 14 résistance d'encrassement (séri e B el C)

<=::::: sens de l'·écoulement

Figure 5. 19 PhOlo.gJ1Iphie du dépôt pOil r un essai de lil série C
( da. '" 3,97 i!lm, Q.a. = 25 4fh. Tgm = 148 ·C. Ts = 39 ·C. dur ~ 2.7 .3 h }

Il.8
 Chap 5 Etude du dépôt sur un faIsceau de tubes QI/etés

50~------------------------------------------------------~
45
-40
Ct da = 3.971Jm
-CD 35 série C
'! 30
8. 25
'CD
"20
1
o
15
E 10
5
Oe=~--~-----+----~------~----+-----~-----+----~----~

o 100 200 300 400 500 600 700 800 900

masse Injectée ( g )

Figure 5 20 Evolution de la masse de dépôt des essais B et C

SERIE - D.

Comme pour les essais de la série précédente, nous avons effectué les essais avec les particules
de 3.97 '.un, mais un débit massique du gaz à 43 kg/h. En ce qui concerne la résistance
d'encrassement (figure 5 22), on remarque qu'elle est d'autant plus faible que le débit est
important. Dans ces conditions expérimentales, le dépôt est très localisé (figure 5.23)

La figure 5.24 montre l'évolution de la masse déposée. Par rapport aux essais de la série
précédente, on remarque que la vitesse de dépôt est plus importante dans ces conditions
expérimentales. Avec une concentration en particule de 0 407 glm3 et un débit de 23.2 gIh, on
estime la vitesse de dépôt est de l'ordre de 7.36 10-3 mIs

119
 Chap 5 Etude du dépôt sur un faisceau de tubes aile tés

0.012 . . . . - - - - - - - - - - - - - - - - - - - - - - - - - - - - ,
Rf =0.0089 ( 1 - exp ( - 0.0044 mi ) }
0.01

~ 0.008
~
série - 0:
1i
-
~
0.006

0.004
da = 3.971Jm
Tgm = 153 oC
Ts =42 oC
Um =9.8 mIs
0.002

O~-------+--------+--------+--------r-------~------~
o 100 200 300 400 600
masse Injectée ( g )

Figure 5 21 Evolution de la résistance d'encrassement (série 0)

0.03

0.025

-
~
~
0.02

~ 0.015 série C
- 0.01
Qa = 25kg/h
~

0.005
!'!:== __ r---sérle 0
Qa =43 kg/h
0
0 100 200 300 400 500 600 700 800 900
masse Injectée ( g )

Figure 5.22 Evolution de la résistance d'encrassement (série Cet D)

120
 *"" sen~ de l'é:cQulement

Figure 5,23 Photographie du dépôt. pour un essai d,é la série D

( da = 3.97 p.m, Qa = 43 kgfh, Tgm = 155"C, Ts = 40 oc. dur = lU Il )

OOr-------------------------------------~
45
.- ; 4)
- 35

i~
1 20

I:~5
o -
o 100 300 «Xl 50J 6CO 100
1'IlOS$e.1rI}ed" (SI )

Figure 5,24 Evolution de la. masse d dépôl des essais C ct D

5.4 ESSAIS DE DEPOT EN ISOTHERME

Pour examiner l'inJ]Il.énce de laihermophorèse nouS elfeciuons ([es essais en ool'ldition~

isothermes Cl Comparons les résultats obtenus avec les résuliats des essais en ul!Ilsfert de
chaleur. Nous effectuons des essais à la. tempèrature de 80 .oC en utilisant la même prooédufC
d'C$SlIique edle pOUf dC$ essais pr6oéd.C!'I1s.

121
 Chap.5 Etude du dépôt sur un/alsceau de tubes alJetés

Nous avons choisi de travailler à même vitesse moyenne dans le faisceau que lors des essais
correspondants en transfert de chaleur ce qui nous a conduit à modifier le débit massique de
gaz.

SERIE· ISOA.

Nous avons réalisé une série d'essais en utilisant des particules 2.20 f.1m et nous avons choisi le
débit massique du gaz égal à 16 kg/h qui correspond à la vitesse moyenne du gaz de 3 1 mis
(identique à celle des essais de la série A).
Sur la figure 5.25, on remarque que par rapport aux essais de la série A les masses recueillies
sont nettement plus faibles. Dans ces conditions, on observe que la vitesse de dépôt est deux
fois inférieure à celle en régime transfert de chaleur (0.48 10-3 m / s en isotherme et
0.80 10-3 m / s pour la série A)

da = 2.20 ~m
Um = 3.1m/s

série - A

100 200 300 400 500 600 700 800 900

masse Injectée ( g )

Figure 5.25 Evolution de la masse de dépôt des essais séries: ISOA & A

SERIE· ISOB.

Pour cette série d'essais, nous avons utilisé des particules de même taille que celles de la série
précédente, mais le débit massique du gaz est fixé à 30 kg/h ce qui correspond à la vitesse
moyenne du gaz de 5.7 mis (identique à celle des essais de la série B).
Les résultats obtenus (figure 5.26) montrent la même tendance en ce qui concerne les masses
de dépôt recueillies. On estime que dans ces conditions la vitesse de dépôt de l'ordre de
1.53 10-3 m / s (elle est en transfert de chaleur de l'ordre de 2.22 10-3 m / s).

122
 Chap 5 Etude du dépôt sur unfatsceau de tubes attelés

OOr-------------------------------------------~
45 da = 2.20 IJm
-40
tJ) Um =5.6 mIs
- 35
,1 30 série - B
8.
25
'G)
"20
8E 15
10
5
...
O~~~----_r----~----~--~----~----+_----~--~
o 100 200 300 400 500 600 700 800 900
masse Injectée ( 9 )

Figure 5 26 Evolution de la masse de dépôt des essais séries ISOB & B

SERIE·ISOC.

Nous avons effectué des essais avec les particules de 3.97 Ilm et avec un débit de 30 kg/h. Les
résultats montrent que par rapport aux essais correspondants en transfert de chaleur (figure
5.27), on n'observe pratiquement aucune différence significative en ce qui concerne la masse
déposée notamment pour les essais de courte durée. Il n'y a donc dans ce cas pas d'effet dû à la
thermophorèse

OOr-------------------------------------------------~
45 da = 3.971Jm série - C
-40
tJ) Um = 5.6 mIs
-35
·1 30
.G)
8. 25
"20
; 15
o
E 10
5
O~----~--~-----+----~----~----+---~~---+----~
o 100 200 300 400 500 600 700 800 900
masse Injectée ( 9 )

Figure 5.27 Evolution de la masse de dépôt des essais séries' ISOC & C

123
 Chap.5 Etude du dépôt sur un jatsceau de tubes attetés

SERIE - ISOO.

Pour cette série d'essais, nous avons utilisé les mêmes particules que celles de la série
précédente, mais le débit est plus fort (51 kglh). Les résultats obtenus comparés avec les
résultats des essais de la série D en régime de transfert de chaleur sont donnés sur la figure
5.28 . En ce qui concerne la vitesse de dépôt, on obtient pour ces conditions expérimentales la
valeur de 6.70 10-3 ml s (7.36 10-3 ml s en transfert de chaleur)

OOr-------------------------------------------------~
45 da = 3.971Jm
-40
Q Um =9.8 mis
-35 série - 0
-1 30
8.
-G)
25
'020
~ 15
E 10
5
O~----~--~-----+----~----~----+---~~--~----~
o 100 200 300 500 700 800 900
masse Injectée ( g )

Figure 5.28 Evolution de la masse de dépôt des essais séries· ISOD & D

En conclusion, d'une manière générale la répartition du dépôt sur les tubes considérés est la
conséquence de la configuration de l'écoulement dans le faisceau. Les particules ont tendance à
se déposer sur la zone où les coefficients d'échange sont importants Sur la première rangée de
tubes, la masse déposée est plus importante que celle sur les rangées suivantes

Les résultats expérimentaux montrent que pour les conditions limites thermiques données la
résistance d'encrassement dépend principalement de la taille des particules et du débit de gaz
traversant la section d'essai. Les plus petites particules produisent une résistance thermique
plus importante. Lorsque le débit est plus faible, la résistance d'encrassement est plus
importante.
La masse déposée sur la surface des tubes augmente avec la masse des particules injectées. Pour
les essais de courte durée, l'évolution de la masse de dépôt est pratiquement linéaire. Pour les
essais réalisés avec les plus petites particules, l'effet de la thermophorèse est mis en évidence.

124
 Chap.5 Etude du dépôt sur unfOlsceau de tubes OI/etés

5.5 ESTIMATION DU DEPOT

Nous avons remarqué sur les résultats expérimentaux que la répartition du dépôt sur les tubes
à ailettes est l'image de la configuration de l'écoulement dans le faisceau de tubes L'épaisseur
du dépôt varie suivant la position sur la surface des tubes et des ailettes et dépend des
conditions limites appliquées
D'une manière générale, le dépôt est par nature transitoire Au cours de l'accumulation du
dépôt, les coefficients d'échange évoluent et dépendent des caractéristiques du dépôt en
formation. Dans ces conditions, la modélisation du dépôt et de son évolution est
particulièrement complexes.

Dans le cadre d'une approche globale, on se propose dans les paragraphes suivants d'estimer la
masse des particules déposée sur la paroi en appliquant les différents modèles de dépôt
existants. Dans ce cas, on s'intéresse aux particules déposées par diffusion avec et sans
thermophorèse, impaction inertielle, interception et sédimentation

DEPOT PAR DIFFUSION

Nous avons vu au chapitre IV que dans le cas du dépôt sur un tube isolé l'analogie entre les
transferts de chaleur et de masse par diffusion est vérifiée. Nous essaierons de vérifier de
manière qualitative si le concept d'analogie est applicable dans le cas d'une configuration de
tubes dans un faisceau.

Sur les photographies des figures 5.10, 5.15, 5.19 et 5.23, nous avons remarqué que d'une
manière générale le dépôt est maximum sur la zone de points d'impact de l'écoulement et
minimum sur la zone de points de décollement
En comparant avec les résultats sur les figures 5 4 à 5 6, on peut remarquer que la position des
maximum et minimum d'épaisseur du dépôt correspond respectivement à la position des zones
échanges de chaleur maximum et minimum. Nous sommes donc en mesure de penser que
l'analogie entre les transferts de chaleur et de masse est vérifiée

En considérant que la surface des tubes est initialement propre et lisse hydrauliquement, on se
propose alors d'estimer la masse des particules déposée par diffusion En plus des hypothèses
concernant les transferts de chaleur et de masse, on admet également que les coefficients
d'échange restent constants pendant cette phase initiale de dépôt
Nous retenons pour cette approche les données sur les figures 3.8 et 3 9 (page 76) pour la
distribution en taille des particules considérées.

125
 Chap.5 Etude du dépôt sur un faisceau de tubes Olletés

Si on suppose que chaque classe de taille des particules est indépendante, on peut alors estimer
le coefficient moyen de transfert de masse par :
N
hm =~
,=1
hm., (OP.' (5.7)

hm.1 représente le coefficient de transfert pour une classe de taille de particule, et (OP.I la
fraction massique de la classe de taille de particule donnée.
Nous calculons h m •1 par l'intermédiaire du nombre de Stanton pour le transfert de masse .
h
St m., =~
U (5.8)
m

U m représente la vitesse moyenne du gaz traversant la section de passage minimal dans le

faisceau En appliquant l'analogie entre les transferts de chaleur et de masse, nous déterminons
Stm,i à l'aide de l'expression (5 1)

Stm = 0 024 Re~'04 Sc~"' (5.9)

Dans ce cas. le nombre de Reynolds est donné par

Re= U m dt (5 10)
v
dt le diamètre de tube et v la viSCOSité cinématique évaluée à la température de film

Pour tenir compte de la thermophorèsc. nous calculons le nombre de Stanton pour le transfert
de masse thermophorètique par l'expression (2.54) Dans ce cas, comme dans le cas du dépôt
sur tube isolé, nous estimons le facteur d'augmentation thermophorétique à l'aide des modèles
donnés au paragraphe 22.1.3.
Cependant, le calcul nécessite la connaissance de la température de la surface des tubes ailetés
que l'on conn ait pas à priori. Nos conditions expérimentales ne nous permettent pas de les
mesurer. Nous avons seulement mesuré lors de l'expérience la température moyenne de l'eau de
refroidissement des tubes
Mais. compte tenu de la dimension du tube (20 mm de diamètre et 40 mm de longeur) et
comme le débit d'eau de refroidissement de tube est suffisament important (300 I/h), on peut
alors supposer que la température moyenne de la surface des tubes est identique à celle de l'eau
de refroidissement.

D'un autre côté, nous avons vu au paragraphe 5 1.2 que la température locale sur la surface de
l'ailette varie suivant la direction radiale et azimutale et elle change d'une rangée à une autre
dans le faisceau Cependant, en vue d'une approche globale, on peut déterminer la température
moyenne de la surface d'échange T. à l'aide de relation

(5 11)

126
 Chap 5 Etude du dépôt sur un jmsceau de tubes m/etés

T8I représente la température moyenne de la surface d'ailette, Tba la température moyenne de la

base de l'ailette supposée égale à la température d'eau de refroidissement des tubes, AlOI et A8I
représentent respectivement la surface d'échange totale et la surface ailette seule.
Nous estimons TOI en utilisant le terme d'efficacité d'ailette donné par'

(5 12)

Les valeurs de T\OI ont été déterminées lors des essais en conditions propres (annexe A)

Résultats du calcul:

Nous avons calculé pour des conditions moyennes correspondantes à nos essais les coefficients
moyens de transfert de masse (ou vitesse de dépôt) par diffusion dont les résultats sont donnés
sur le tableau 5 1

DIFFUSION
essâis da Qa Tgm Ts Um Re - Pr Stm,t hm
série: CUm) (kg/h) (OC) (OC) tln/s) (mis)

A 2.2 14 139.2 42.0 3.1 3140 0.70 8.5E-G4 2.6E-G3

1 2.2 25 144.6 42.0 5.6 5468 0.70 8.0E-G4 4.5E-G3
C 3.91 25 144.0 40.1 5.6 5549 0.70 6.5E-G4 3.6E:œ
D 3.91 43 152.6 41.9 t .. 9439 0.70 6.oE-G4 5.9E-G3

lsoA 2.2 16 80.. 80.1 3..1 2959 0.70 6.1E.Q6 1.89E-DS '
IsoB 2.2 30 81.6 81.2 5.7 -, S420 0.70 2.7E.Q6 1.57E-DS
lsoC 3.91 30 83.0 81.7 5.7 5415 0.70 7..3E.Q6 4.16E-DS
IsoD 3.91 51 80.0 79..3 93 9299 0.70 3.2E.Q6 3.09E-DS

Tableau 5.1 Résultats de calcul du dépôt par diffusion.

DEPOT PAR IMPACTION INERTIELLE

Nous allons à présent estimer la masse des particules déposée par l'impaction inertielle.
Compte tenu de la structure globale de l'écoulement dans un faisceau de tubes en arrangement
aligné, il semblerait que si les conditions sont convenables l'impaction inertielle des particules
ne se produit que sur les tubes situés en première rangée.

127
 Chap.5 Etude du dépôt sur unfOlsceau de tubes OI/etés

En considérant que l'impaction inertielle ne se produit que sur la première rangée, nous
estimons l'efficacité de collecte T\unp à l'aide de l'expression (2 88) .
T\unp = ~1In(8Stkeff) + ~2 (Stk eff - 0.125) + ~3 (Stk eff - 0 125)2
avec: ~1 = 001978749, 132 = 0.5136545, ~3 = -00482858

Stk eff représente le nombre de Stokes effectif pour l'impaction inertielle sur une rangée de
tubes donnée par expression (2 91)
Dans ce cas, [dU] représente le gradient de vitesse adimensionné au point de stagnation
dx sp
amont qui sert à tenir compte de l'effet de la disposition des tubes voisins sur l'impaction
inertielle sur un tube situé dans une rangée. Nous estimons la valeur de ce terme en utilisant les
données obtenus par Konstandopoulos et al. [20] représentées sur la figure 2 18.
Compte tenu du fait que la section d'essai utilisée simule un faisceau de tubes infinis, et en
prenant la valeur de S égale à 4, nous retenons alors la valeur de 2 35 pour le terme de [dU] .
dx sp

En ce qui concerne t p (temps caractéristique de la particule), nous calculons celui-ci pour

chaque classe de taille des particules considérées à l'aide de l'expression (2 71)
Comme dans le cas précédent, nous retenons également les données sur les figures 3.8 et 3.9
pour la distribution de la taille des particules. Nous retenons également la valeur de 2679.9
kg/m3 pour la masse volumique des particules de diamètre moyen de 2 20 /lm et la valeur de
2525.5 kg/m3 pour celle des particules de diamètre moyen de 3 97 /lm

Cependant, pour tenir compte de la présence des particules qui ne suivent pas la loi de Stokes
(particules non Stokesiennes), nous calculons t p avec le facteur de correction de 'l'(Re p )
déterminée par l'expression (2 86)
Nous calculons la vitesse de dépôt par impaction inertielle V unp en utilisant la relation suivante:

Vunp (Aunp
"'hmp -
- U -A -) (5 13)
m proJ

AproJ la surface apparente des tubes et Aunp la surface nécessaire pour l'impaction inertielle.

Résultats du calcul:

Nous avons calculé séparément l'efficacité de collecte par impaction inertielle sur la partie du
tube et sur la partie des ailettes Dans le premier cas, le nombre de Stokes est basé sur le
diamètre du tube Sur le second, il est basé sur le diamètre de l'ailette. Les résultats de calculs

128
 Chap 5 Etude du dépôt sur un fmsceau de tubes a/Jetés

pour des conditions correspondantes à nos essais sont donnés dans le tableau 5.2 Dans ce cas,
l'indice (1) est pour le tube, (2) pour l'ailette

IMPACnON INElll'lEUE
essais da QG Tgm Ts Um etc. VdJmpl etc. Vd.imp2
série: (Jm) (kg/h) rC) rC) tin/$) tube tube ailette ailette
(Ok) (mIs) (%) (m/S)
A 2.2 14 139.2 Q,O 3.1
B 2.2 25 14A.6 Q,O 6.6
C 3.97 25 144.D .,.1 6.6 0.D2 I.6E44
D 3.97 oC3 152.6 41.9 ,.1 0.24 2.3E~ 0.01 7.8E44

IsoA 2.2 16 10.1 10.1 3.1

Isoi 2.2 30 lU lU 5.7
!soC 3.97 30 l1D 11.7 1.1 0.01 6.GE44
IsoD 3.97 51 10.0 79.3 9.1 0.19 UE~ '0.002 2.1E44
Imp1 : Impactlon inertie. . . . tube efo : efftcacl" d. collecte
.
",,','1'::, .. ",'11',11
,
l

. ~""
' ,
1mp2 : ImDactlon inertie.......... "

Tableau 5 2 Résultats de calcul du dépôt par impaction inertielle

DEPOT PAR INTERCEPTION

Comme dans le cas du dépôt sur le tube isolé, nous estimons pour chaque classes de tailles de
particules l'efficacité d'interception à l'aide de l'expression (2 92). Nous calculons ainsi
l'efficacité totale par:
N

l1mt = L
1=1
l1mt.1 CJ)P.I (5 14)

Nous déterminons enfin la vitesse de dépôt par interception Vmt par la relation suivante .

(5 15)

A mt est la surface nécessaire pour l'interception

Résultats de calcul:

Comme dans le cas précédent, nous avons calculé séparément l'efficacité de collecte par
interception sur la partie du tube et sur la partie de l'ailette, et les résultats de calcul sont
donnés dans le tableau 5.3

129
 Chap.5 Etude du dépôt sur unJOIsceau de tubes OIletés

INTERCEPTION
essais da Qa Tgm Ts Um .'c. Vd.lntl efc. Vd.lnt2
tube tube ailette ailette
,
sene:
. Ci.a'n) (kg/h) rC) (OC) Cm/s) (Ok) (mIs) (Ok) (mIs)
A 2.2 14 139.2 42.0 3.1 0.015374 üE-04 0.009 2.7E-04
B 2.2 25 - 144.6 42.0 5.6 0.015374 1.6E-04 0.009 4.8E-04
C 3.97 25 144.G 40.1 5.6 0.G25152 lAEoOO 0.014 7.8E-04
D 3.97 43 152.6 41.9 9A 0.G25152 2.5EoOO 0.014 lAEoOO

IsoA 2.2 16 80.8 80.1 3.1 0.015374 .&.8E-04 0.009 2.6E-04

IsoB 2.2 30 81.6 81.2 1.7 0.015374 I.8E-04 0.009 4.9E-04
!soC 3.97 30 83.0 81.7 1.7 0.G25152 lAEoOO 0.014 8.0E-04
Iso 0 3.97 51 80.0 793 9.7 0.G25152 2AEoOO 0.014 lAEoOO
Inn : IntercepHon lUI tube etc : eftIcacltt de coUecte
1nt2 : IntarcepHon sur ailette

Tableau 5.3 Résultats de calcul du dépôt par interception.

DEPOT PAR SEDIMENTATION

Pour la sédimentation, nous estimons la vitesse de dépôt par sédimentation Vsed à l'aide de
l'expression (2.97) .
vséd = gt p
Les résultats de calcul de la vitesse de sédimentation pour les conditions moyennes
correspondantes à nos essais sont donnés dans le tableau 5.4

SEDIMENTAnON
essais da Tgm Ts Um .fc. Vd.séd
série: rC)
A 2.2 14 139.2 42.0 ' 3.1 0.0236 7.4E-D4
8 r 2.2 2S 144.6 42.0 ' 6.6 D.02351 1.3E-D3
C 2.2 2S 144.0 40'.1 6.6 0.D5494 3.1E-D3
D 3.97 43 162.6 41.9 9.1 0.D5446 6.3E-D3

IsoA 2.2 16 80.1 80.1 3.1 0.D2371 7.4E-D4

Iso8 2.2 30 81.6 81.2 6.7 0.D2367 1.4E-D3
IsoC 3.97 30 83.0 81.7 6.7 0.D5538 3.2E-D3
IsoD 3.97 51 80.0 79.3 9.7 0.05567 6.4E-D3

Tableau 5.4 Résultats de calcul du dépôt par sédimentation.

130
 Chap 5 Etude du dépôt sur unfmsceau de tubes mletés

5.6 COMPARAISON DES RESULTATS

Nous avons estimé dans les paragraphes précédents la masse des particules déposées sur les
tubes en appliquant les différents mécanismes de dépôt Les résultats des calculs montrent que
pour nos conditions expérimentales l'interception et la sédimentation sont pratiquement
négligeables. Pour les conditions des essais des séries A,B et C, on peut constater que la
diffusion et la thermophorèse sont les mécanismes de dépôt prédominants
Pour les essais de la série B particules 3.97 ~m, débit 25 kglh), on remarque déja l'apparition
du régime de dépôt controlé par l'impaction inertielle. Ce mécanisme de dépôt est plus
important dans les essais de la série D.
En régime isotherme, on remarque également que pour les essais ISO-D, l'impaction inertielle
est le mécanisme de dépôt prédominant Dans les conditions ISO-A, ISO-B et ISO-C, ce sont
l'interception et la sédimentation qui sont les plus importants

essais da Um Re Vd.xp Vd.th

série: (

A 2.2 14 139.2 ~.o 3.1 31., a.OE-04 3.6E-03

B 2.2 25 144.6 ~.o 5.6 S46I 2.2E-03 6.0E-03
C 3.97 25 144.0 ".1 5.6 &&9 3AE-03 7.2E-03
D 3.97 43 152.6 n9 9.1 MW 7AE-03 3.4E-02

IsoA 2.2 16 10.1 10.1 3.1 2969 ~-G4 1.0E-03

1so8 2.2 30 lU lU U 5GO 1.5E-03 1.6E-03
IsoC 3.97 30 13.Q 11.7 U 5&li 3.5E-03 3.4E-03
IsoD 3.97 SI 10.0 79.3 9~ 9299 6.7E-03 2.3E-02
Vd.xp : vtt.... de ~p6t exp6rlmentale
Vd.th : vit. . . de

Tableau 5 5 Comparaison des vitesses de dépôt en régime de transfert de chaleur

Le tableau 5.5 donne les valeurs des vitesses de dépôt calculées et celles déterminées
expérimentalement. En comparant ces résultats, on peut remarquer qu'ils ont globalement le
même ordre de grandeur.
On remarque cependant que pour les conditions de la série D l'écart observé est assez
important Cela provient du mécanisme de rebond des particules qui accompagne l'impaction
inertielle.

131
 Chap 5 Etude du dépôt sur unfmsceau de tubes mletés

5.7 CONCLUSION

Nous avons examiné le comportement du dépôt sur un faisceau de tubes à ailettes en réalisant
des essais en régime de transfert de chaleur et en isotherme. La répartition du dépôt sur les
tubes correspond à l'image de la configuration de l'écoulement dans le faisceau Les particules
ont tendance à se déposer sur la zone où les coefficients d'échange sont importants

Pour les essais effectués avec un débit d'aérosols plus faible, la résistance d'encrassement est
plus importante De même la résistance thermique du dépôt est plus importante lorsque l'on
utilise les plus petites particules Pour les essais effectués avec le débit le plus important, le
dépôt est localisé dans les zones d'impact sur les tubes et les ailettes

En appliquant les différents modèles de dépôt, il a été possible d'estimer la masse des particules
déposées sur la paroi des tubes Compte tenu des nombreuses hypothèses et des incertitudes
expérimentales non négligeables, on peut constater que l'approche effectuée permet d'obtenir
des résultats acceptables

132
CONCLUSION GENERALE

L'étude de l'encrassement particulaire sur la surface des tubes d'échangeur de chaleur a été
effectuée en simulant le dépôt d'encrassement à l'aide de particules de taille uniforme Ce choix
a permis de mettre en évidence les différents mécanismes de dépôt.

Deux sections d'essais ont été utilisées.

La première est constituée par un tube simple disposé perpendiculairement à l'écoulement
d'aérosols Le bon contrôle des conditions expérimentales parfaitement reproductibles a permis
par des essais répétés de durée variable d'obtenir l'évolution de la masse déposée en fonction
du temps. Trois régimes de dépôt ont été mis en évidence Au cours du premier, la vitesse de
croissance est pratiquement constante mais faible. Le deuxième régime commence avec une
vitesse de croissance plus élevée qui ensuite s'atténue Pendant le troisième, l'évolution semble
présenter un caractère asymptotique.
L'interprétation de cette cinétique a été réalisée en considérant les mécanismes élémentaires de
dépôt par sédimentation, interception, diffusion convective avec thermophorèse et impaction
inertielle. Il a été vérifié que, pour nos conditions expérimentales, l'impaction inertielle
n'intervient pas La diffusion convective avec thermophorèse a été estimée en s'appuyant sur le
modèle proposé par ROSNER, elle contribue pour la part la plus importante au mécanisme de
dépôt pendant le premier régime et devient prépondérante pour le régime suivant.
L'augmentation de la vitesse de croissance au début de ce régime a été attribuée aux rugosités
créées par le dépôt initial sur la surface du tube lisse Nous avons ensuite démontré que
l'isolation thermique créée par le dépôt sur la surface du tube contribuait à l'atténuation de la
thermophorèse et donc de la vitesse de croissance Moyennant l' hypothèse simplificatrice d'un
dépôt uniforme de masse volumique connue, le calcul de l'évolution de la température de
surface du dépôt et de sa conductivité thermique apparente a conduit à des valeurs réalistes.
L'étude du dépôt sur la géométrie simple du tube isolé a donc permis de mettre en évidence les
mécanismes élémentaires contrôlant le dépôt.

La deuxième section d'essais est constituée d'un faisceau de tubes à ailettes Nous avons étudié
le comportement du dépôt sur cette géométrie plus complexe en élargissant la gamme des
conditions expérimentales pour réaliser des essais avec impaction inertielle. De plus chacun des
essais a aussi été réalisé en condition isotherme
Les essais réalisés avec transfert de chaleur ont permis d'obtenir, à partir de mesures globales
permettant le calcul du coefficient de transfert thermique global, l'évolution de la résistance

133
d'encrassement de la maquette d'échangeur Il faut la aussi souligner la bonne reproductibilité
des essais de durée différente Toutefois les courbes donnant la cinétique de croissance du
dépôt en fonction du temps sont essentiellement linéaires alors que celles donnant l'évolution
de la résistance d'encrassement présentent un caractère asymptotique
Lorsque l'impaction inertielle est négligeable, la répartition du dépôt sur la paroi correspond à
l'image de la configuration de l'écoulement dans un faisceau de tubes Les particules tendent à
se déposer dans la zone où les coefficients d'échange sont importants
La résistance d'encrassement est plus importante pour des essais effectués avec les plus petites
particules. Lorsque le débit du gaz encrassant est important, l'impaction inertielle n'est plus
négligeable, le dépôt est alors très localisé, particulièrement aux niveaux des impacts de
l'écoulement sur les tubes et les ailettes. Dans ces conditions, la résistance d'encrassement est
moins importante
Les vitesses de croissance obtenues à partir des courbes de cinétique ont pu être interprétées
en les comparant aux vitesses de dépôt de chaque mécanisme élémentaire La concordance est
convenable aussi bien en régime de transfert de chaleur qu'en régime isotherme sauf lorsqu'il y
a impaction inertielle car alors le phénomène de rebond des particules sur les tubes et ailettes
ne peut pas être modélisé.

Sur la base des résultats obtenus, il est possible de proposer quelques axes de travail futur
-en ce qui concerne le dépôt sur un cylindre isolé:
-l'étude expérimentale du dépôt de particules monodispersées sur des surfaces de
rugosité initiale connue,
-la mesure et la prise en compte sur notre dispositif expérimental de la forme réelle du
dépôt.
-en ce qui concerne le dépôt sur maquette d'échangeur réel, l'influence de la disposition des
tubes et du nombre d'ailettes par tube.

134
BIBLIOGRAPHIE

[ 1] GRaLOT J M, "Etude du dépôt de particules en phase gazeuse dans des canaux

d'échangeur thermique à plaques", thèse Doct-INPG, 1989.

[2] DUFF AU C, GRaLOT lM, NAV ARRO JM, VIDIL R, "L'encrassement dans les
échangeurs: Description, Prévention, Remèdes", GRETh- Grenoble, 1985

[3] KERN D Q, SEATON R.E : A Theoretical Analysis of Thermal Surface Fouling . Brit.
Chem. Eng, Vol. 4, No.5, pp 258- 262, 1969.

[4] ROSNER D.E, GOKOGLU S.A, ISRAEL R, "Rational Engineering Correlation of

Diffusional and inertial Particle Deposition behaviour in non isothermal forced Convection
environments", Fouling of Heat Exchanger surfaces, Proceedings of the Engineering
Foundation Conference, White Haven, Pennsylvania BRYERS RW, 1982

[5] ISRAEL R, ROSNER D.E, "Use of a Generalized Stokes Number to Determine the
Aerodynamic Capture Efficiency of Non- Stokesian Particles from a Compressible Gas Flow",
Aerosol Sci & Technol. 2, 45-51, 1983.

[6] INCROPERA F, DE WITT, "Fundamentals of heat and mass transfer", 3rd Ed, J Wiley
and Sons, New York, 1985.

[7] ROSNER, D.E, "Transport Processes in Chemically Reacting Flow Systems", Stoneham .
Butterworths, p.47, 1986.

[8] MARNER VI.J, "Progress in Gas Side Fouling of heat Transfer Surfaces", Appl Mech.
Rev, Vo1.43,p 35-36,1990

[9] TALBOT L, "Thermophoresis - A reView American Institute of Aeronautics and

Astronautics", The Twelfth international Symposium on rarefied Gas dynamics Paper 195.
Charlottesville Inc, p. 467-488, 1980.

[10] DERJAGUIN B V, RABINOVICH l, STOROZHILOVA H.I and SHCHERBINA G.I,

"Measurement of the coefficient of thermal slip of gases and the thermophoresis velocity of
large size aeorosol particles", J colloid Interface Science, Vol. 57, p.451-461, 1976

135
[Il] ROSNER D.E and PARK HM, "Thennophoretically Augmented Mass, Momentum and
Energy Transfer Rates in High Particle Mass- Loaded Laminar Forced Convection Systems",
Chem.Engineering Science, Vol. 43, p 2689 - 2704, 1988.

[12] HESKETH HE, "Fine Particles in Gaseous media", second edition, Chelsea, Michigan,
Lewis Publishers Inc, 1986.

[13] ROSNER D.E, Fernandez de la MORA,"Small Particles transport across Turbulent non
Isothermal Boundary Layers", ASME Journal Engineering for Power, Vol 104, 1982

[14] FRIEDLANDER SK, "Smoke, Dust and Haze Fundamentals of aerosol behaviour",
edition J WILEY and sons, New York, 1977

[15] BRUN R J, LEWIS W, PERKINS P.J and SERAFINI J S, "Impingement of cloud

droplets on a cylinder and procedure for measuring liquid water content and droplet sizes in
supercold clouds by rotating multicylinder method, National advisory committee for
aeronautic, Report No. 1215, p 141-183, 1955.

[16] WESSEL R A et RIGHI J, "Generalized Correlations for Inertial Impaction of Particles

on Circular Cylinder ", Aerosol Science and Technology, No. 9, p 29-68, 1988

[17] ILlAS S, DOUGLAS P.L, "Inertial Impaction of Aerosol Partic1es on Cylinder at

Intermediate and High Reynolds Numbers", Chemical Engineering Science, Vol 44, No. 1,
1989.

[18] CHOUDARY K.R, GENTRY J W, "A Model for Particle Collection with Potential Flow
between Two ParaUei Cyliders", The Canadian Journal of Chemical Engineering, Vol 55, p.
403-407, 1977.

[19] INGRAM D.B, HILDYARD M.L, HEGGS P.J, "The Particle Collection Efficiency ofa
Cascade of Cyliders", The Canadian Journal of Chemical Engineering, Vol 67, p 545-553,
1989.

[20] KONST ANDOPOULOS, LABOWSKY, ROSNER, "Inertial Deposition ofPartic1es from

Potential Flows Past Cylinder Arrays", J Aerosol Science, H TCRE # 187, 1993

[21] FUHS S.E, "Studies of inertial deposition of particles onto heat exchanger elements",
Thesis, California Institute ofTechnology, 1988

136
[22] GLEN NF and HOW ARTH J.H, "Modeling refuse incineration fouling ", 2nd United
Kingdom Conf on Heat Transfer, Institution ofMech. Eng London, No 1, pAOI-420, 1988

[23] EPSTEIN N, " Fouling . Technical aspects", Fouling of heat transfer equipment
SOMERSCALES EFC., KNUDSEN JG, Hemisphere Wash DC, 1981

[24] EPSTEIN N, " Particulate Fouling ofHeat Transfer Surfaces' Mechanisms and Models,
L .F Melo et al.Eds, Fouling Science and Technologie, p 143-164, Kluwer, 1988.

[25] MARNER W J, SUIT OR J W, "Fouling with Convective Heat Transfer", Handbook

of Single Phase Convective Heat Transfer KAKAC, SHAH, AUNG, ch 21, 1987

[26] TABOREK J, AOKI T, RITTER R.B, PALEN J W and KNUDSEN J G, "Fouling The
Major Unresolved Problem in Heat Transfer, Chemical Engineering Progress, Vol 68, No.2, P
59-67, 1972.

[26a] TABOREK J, AOKI T, RITTER RB, PALEN JW and KNUDSEN J.G, "Predictive
methods for fouling behavior, Chemical Engineering Progress, Vol. 68, No 7, p. 69-78, 1972.

[27] CLEAVER J W and YATES B, "Mechanism of detachment of colloidal particles from a

flat substrate in a turbulent flow,Journal of Colloid and Interface Science, No 44, p 464-474,
1973.

[27a] CLEAVER J.W and YATES B, "The effect of reentrainment on particle deposition,
Chemical Engineering Science, No 31, p 147-151, 1976.

[27b] CLEAVER J. W and Y ATES B, "A Sub layer model for the deposition of particles from
a turbulent flow, Chemical Engineering Science, No 30, p.983 - 992, 1975

[28] EL SHOBOKSHY M.S, EL ZAYAT R.E, "Deposition of airborne Particles Onto the
Surfaces of air-cooled finned tubes", Canadian Journal of Chemical Engineering, Vol. 59,
1981.

[29] BOTT T R, BEMROSE C R, YEO B.K, "Particulate fouling on the gas side of finned
tube heat exchangers", Journal ofHeat Transfer, Vol 105, 1983

137
[30] REAY D, ALLEN R W, BULL K, CURSLEY C.J, "Fouling of a finned tube heat
exchanger by micron and submicron size particles" AERE HARWELL G 2946 Engineering
Science Division m... 84/2649, 1984

[31] LOKSHIN V.A, FOMINA V.N, PAVLOV N V, PORTYANKO AA, SKORNY

AKOVA L.G, "Experimental investigation of the fouling oftranversely finned bundles of tubes
in a dust laden air flow", Thermal Engineering, Vol 27, No. 6, 1980

[32] PORTYANKO A.A, LOKSHIN V.A, FOMINA V.N, "The fouling oftransversely swept
bundles of tubes with welded spiral tinning", Thermal Engineering, Vol 27, No 12, 1980.

[33] LOKSHIN V.A, ZOZUL y A N V, FOMINA V.N, PORTYANKO A.A, "Tests of an

economiser manufactured of tubes with welded strip-type fins in a boiler operating on
Ekibastuz coal", Thermal Engineering. Vol. 26, No 3, 1979

[34] SILVESTRINI R, "Heat Exchanger Fouling and Corrosion", Chemical Engineering

Progress, Vol. 75, No 12, 1979

[35] CASPER L.A, SEMLER T T. BOGUE J.L, W AITE T H, "Fouling of a finned tube diesel
flue gas heat recuperator", Final repon DOE, 1982.

[36] ROBERTS, KUBASCO, ·Combined cycle steam generator gas side fouling evaluation",
Final report for NS SC department of the Navy PUCN N0024 77C 4386, 1979

[37] BOULAUD D, "Production d'aerosol dans les expériences de projet PITEAS Analyse
préliminaire", CEAlDPT/SPINILPA Fontenay-aux-Roses, Raport Interne SPT, No 260, Sept
1981 .

[37a] BOULAUD D, "Data reduction by EDS method Application to impactor or cyclone",

CEAlDPT/SPINILPA Fontenay- aux-Roses, Rapport interne SPIN, No. 84/9, Mai, 1984

[38] BOULAUD D, BOURBIGOT Y, "Etude des Caractéristiques d'un Pulvérisateur

Acoustique", Projet PITEAS, Rapport interne CEAlDPT/SPINILPA, No 343, 1983

[39] LIU B.Y H, PUI D.Y H, "Electrical neutralisation of aerosols", J. Aerosol Sci, no.5,
p 465-472, 1974.

138
[40] CHARUAU J, "Etude du dépôt de particules dans les conduites", Rapport CEA R 5158,
1982.

[41] FUCHS N.A, "Aerosol Impactors (a review), Fundamentals of Aerosol Science", SHAW
DT WILEY J and Sons, New York, 1978.

[42] SOOLE B.W, "Conceming the calibrations constants of cascade impactors with special
reference to the casella MK2", J Aerosol Sci, No.2, p 1-14, 1971.

[43] FUCHS N.A, SUTUGIN A.G, "Generation and use of monodisperse aerosols dans
Aerosol Science", ed CN DAVIES, Academic Press - London, Chap 1, 1966.

[44] ZUKAUSKAS A, ZruGZDA J, "Heat Transfer of a Cylinder in a crossflow",

G.F.HEWITT, U S.A, 1985.

[45] WOOD N.B, "A Simple method for the calculation of turbulent deposition to smooth and
rough surfaces", 1. Aerosol Sci, Vol 12, No 3, p.275-290, 1981

[46] SHIMADA M, OKUYAMA K, KOUSAKA Y and OHSHIMA K, "Turbulent and

Brownian Diffusive Deposition of Aerosol ParticIes onto a Rough Wall", J. Chem. Eng Japan,
Vol 20, p.57-64, 1987 .

[47] GUDMUNDSSON J.S, "Particulate Fouling in Fouling of Heat Transfer Equipment",

Somerscales EFC & Knudsen JG, eds., pp 357-387, Hemisphere, Washington D.C, 1981

[48] EL-SHOBOKSHY M.S and ISMAIL LA, "Deposition of Aerosol Particles from
Turbulent Flow onto Rough Pipe Wall", Atmos. Environ, Vol 14, No. 3, P 297-304, 1980.

[49] DAVIES J T, "A New Theory of Aerosol Deposition from Turbulent fluids", Chem Eng.
Sci, Vol. 38, p 135-139, 1983.

[49a] DAVIES J T, "A New Theory of the Deposition of colloidal particles from Turbulent
fluids", Annals New York Academy Science, p.313-325, 1983.

[50] GRILLOT J.M, "Fouling of cylindrical probe in crossflow ID a diesel exchaust

enviroment", Seminar 23 Proceedings, Fouling mechanismes Theoritical and pratical aspects,
Grenoble - Editions Europeennes Thermiques et Industrie, p 189, 1992.

139
[51] ZUKAUSKAS A and KARNI J, "High Performance single phase heat exchangers",
London: Hemisphere, Chapitre 14, Heat transfer and flow over finned tubes, p 291-318, 1989

[52] OWEN F.K, "Heat transfer from plain and finned cylinders in crossflow", J Inst. Heating
Vent Eng., Vol 35, No. 1, p.213 - 226, 1967.

[53] IRELAND P T et JONES T. V, "Detailed measurements ofheat transfer on and around

a pedestal in fully developed passage flow", 8 th. Int Heat Transfer Conf San Francisco,
p.975-980, 1986.

[54] GOLDSTEIN R et KARNI J, "The effect ofwall boundary layer on local masse transfert
from a cylinder in crossflow", 1. Heat Transfer, Vol 106, p.260-267, 1984.

[55] SHUZ G et KOTTKE V, "Enhancement ofheat transfert by horse vortice in fin tube heat
exchanger", Eurotherme Seminar No 9 - Heat in Single phase Flows, Bochum, p.9612, 1989.

[56] MERCIER P and VILLAND M, "Generic studies for industrial heat exchanger fouling",
Eurotherm Seminar No.18, Hamburg, 1991

[57] LE GALL R, BOUGRIOU C, SOLECKI J C and LAURENT M, "Local Study of

Temperature Distribution On A Circular Plane Fin", IChMe Symposium Series, No 129, 1993.

[58] JOURNOT P, "Etude de l'encrassement particulaire en phase gazeuse de tubes ailetés

d'échangeur de chaleur", DEA de Génie des Procédés - INPG, CENG - STT - Greth, 1990.

[59] GNIELINSKI, ZUKAUSKAS et SKRINSKA, "Banks of plain and finned tubes", Single
Phase Convective Heat Transfer, Hemisphere Publising Corporation, 1983

[60] CRITTENDEN B.D, DOWNEY 1. L, KOLACZKOWSKI S T, "Acquisition and

interpretation of oil refinery plant data for fouling studies", International Conference on fouling
in process plant, St. Catherine College, Oxford, 1988.

[61] CETIAT, "Propriétés physiques de quelques fluides utilisés dans les échangeurs de
chaleur", Paris - Lavoisier Tec. Doc, p.90, 1982.

140
ANNEXE-A

DETERMINATION DU COEFFICIENT D'ECHANGE THERMIQUE COTE GAZ

Pour un faisceau de tubes en conditions propres, on peut relier le coefficient d'échange

thermique côté gaz avec le coefficient global de transfert de chaleur H par la relation

(Al)

où hl et he représentent les coefficients d'échange des fluides interne et externe, Al et Ae sont

respectivement les surfaces d'échange intérieure et extérieure de la paroi, de et d r sont
respectivement le diamètre extérieure de l'ailette et le diamètre extérieure du tube, 1 est la
hauteur d'ailette, kt la conductivité thermique de tube et T'lg l'efficacité globale de la surface
des tubes ailetés

Pour le faisceau de tubes étudié, compte tenu de la géométrie des tubes (petit diamètre et faible
longueur) et compte tenu du débit de circulation de l'eau de refroidissement des tubes qui est
important (300 l/h), on peut supposer que:

Aussi, comme:

ln(d e / d r ) « 1
2xlk t T'lgheAe

On peut alors calculer le coefficient d'échange thermique côté gaz he par

(A2)

où T'lg représente l'efficacité globale de la surface des tubes ailetés donnée par:

(A3)

avec Am. et T'lm. sont respectivement la surface d'échange de l'ailette seule et l'efficacité d'ailette.
Pour une ailette circulaire plane, on peut calculer l'efficacité d'ailette à l'aide de l'expression
proposée par Gnielinski et al. [59]) .

141
 (A4)

avec le =l+~ où :
2
e m l'épaisseur d'ailette, Te le rayon d'ailette et Tb le rayon de tube, et

(A5)

avec km la conductivité thermique d'ailette.

Nous calculons l'efficacité d'ailette des tubes étudiés en utilisant l'expression (A4) et nous
déterminons le coefficient d'échange thennique côté gaz he par un calcul itératif puisque on
calcule he qui est fonction de Tl g lui même fonction de he'

Par ailleurs, nous calculons le coefficient d'échange global H à l'aide de l'expression (54) :

où Qe représente la puissance calorifique échangée, At la surface d'échange totale des tubes

ailetes et DTLM l'écart logarithmique moyen de température donné par:

avec Tee et Tes respectivement la température de l'air chaud à rentrée et à la sortie de

l'échangeur, Tfe et Tfs la température de l'eau à l'entrée et à la sortie de l'échangeur
Ainsi, nous calculons Qe en utilisant la puissance calorifique cédée par air :

avec me représente le débit massique du gaz, Cp la chaleur massique calculée pour une
température moyenne, et ÔTe l'écart de température du gaz aux bornes de la section d'essai
En effet, la mesure du débit massique du gaz est plus précise (débitmètre massique thermique
avec régulation) que le débit de l'eau de refroidissement des tubes (rotamètre)

142
Nous avons réalisé une série d'essais avec de "l'air propre" dont les mesures ont été effectuées
en faisant varier le débit massique du gaz traversant la section d'essai pour des valeurs entre 6
kglh et 41 kglh ce qui correspond à 1100 < Re < 7200. Dans ce cas, nous maintenons
constantes la température d'entrée du gaz ainsi que la température de l'eau de refrodissement
des tubes à 185 oC et 40 oC respectivement.
Les principaux résultats de ces campagnes d'essais sont donnés sur le tableau Al.

Essais Propres réalisés à Tgaz. 185 oC et Ts • 40 oC

Qa Um (mis) Re Pr aff.ai Hg h.ext Nuh.ext

(kg!h) ~/m2K) CW/m2K)

41.0 9,09 7203.2 0.70 0.748 64.94 95.39 62.10

37.8 .A3 "99.8 0.70 0.761 41.56 88.42 57.62
35.1 7.79 6158A 0.70 0.775 57.69 80.73 52.34
32.2 7.13 5665.7 0.70 0.789 54A3 74.55 ~A4
29.2 6A3 5174.8 0.70 0.802 51.oa 68.48 44.83
26.5 5.81 <1681.5 0.70 0.817 47.83 62.86 41.06
23.6 5.13 G07.1 0.70 0.831 44AI 57.13 37.70
20.1 4.48 3703.8 0.70 0.&47 «J.53 SO.88 33.66
17.8 3.83 3198.9 0.70 0.863 36A6 44.64 29.61
14.9 3.18 2683.2 0.71 0.881 32.07 38.27 .25.53
12.0 2.53 , 2171.2 0.71 0.899 27.24 31.61 21.20
9.0' 1.87 1646A 0.71 0.920 24.89 16.84 16.84
6.0 1.22 ' 1110.9 0.7T 0.943 16.53 11.30 11.3

Tableau Al Résultats d'essais en conditions propres

70

...:::J
e
tIJ
,~ "
50
CD
'C Nuh= 0.024 .Re 1\0.896 .P 1\0.333 ~
'CD .."
1< 40 ,-.-
-
CD
-e 30 \ ~
io"""""
~
cO
,
(,)
s::. 20
I~
:::J
Z
10

o
o lCXXl 3CXXl 7CXXl
nombre de Reynolds

Figure A-I Coefficient d'échange thermique côté gaz du faisceau de tubes étudié.

143
ANNEXE- B

INCERTITUDE SUR LE COEFFICIENT D'ECHANGE DE CHALEUR

Nous avons calculé le coefficient d'échange global H par:

avec Qc la puissance de chaleur échangée, Al la surface d'échange totale de tubes ailetes

DTLM est l'écart logarithmique moyen de température donné par.

où Tee et Tes sont respectivement la température de l'air chaud à l'entrée et à la sortie de

l'échangeur, Tfe et Tfs la température de l'eau à l'entrée et à la sortie de l'échangeur

On peut déterminer l'incertitude relative sur H en fonction de Qc' Tee' Tes, Tfe et Tfs en
considérant Al comme un paramètre :

(B.I)

Puis, en posant [60]

(B.2)

(B.3)

y=----
w (B.4)
(Tee -Tfe )

(B'S)

144
On obtient, après calcul :

BR = BQc +_1
R Qc W ln X
{1~-IIBTfi +Il-~BT.,
Sh0(j ,e +12:....-1IBT
ln X ce
+Il-~IBT
ln X cs
} (B 6)

D'autre part, nous avons calculé Qe à l'aide de l'expression:

où me représente le débit massique du gaz et Cp la chaleur massique à pression constante

En considérant Cp comme un paramètre, on peut exprimer l'incertitude sur Qe par

(B 7)

ou:

(B 8)

ou encore :

(B.9)

Mesures des températures:

Nous avons pris soin pour la mesure des températures d'utiliser des thermocouples chromel-
alumel de 0.5 mm de diamètres plongeant de quelques centimètres dans le fluide Le
calorifugeage de la section d'essais est soigné et les points de mesures sont multipliés.
Les thermocouples ont été étalonnés individuellement dans les conditions de mesure en 25
points entre 0 et 250 oC. De ce fait, l'erreur absolue est certainement inférieure à 0.1 oC.

Mesure des débits :

Les deux débitmètres utilisés, étalonnés par les constructeurs, ont tous une précision annoncée
inférieure à 1 % de la mesure. La mise en série a montré que ce résultat n'était que
partiellement vérifié [1] Après avoir recalculé les courbes d'étalonnage, on peut prendre

145
comme erreur global 1 % de la mesure pour les mesures absolues et 0 2 % pour les mesures
relatives.

Application numérique :

Pour les conditions d'une campagne d'essais donnant :

Tee = 188 4 oC, Tes = 95.8 oC, Tfc = 379°C et Tfs = 39.3 oC
et pour une incertitude relative de 0.2 % sur le débit d'air et une incertitude de 0 1 oC sur la
température, on peut obtenir une incertitude relative sur la valeur de H de 0 64 %

Par ailleurs, nous avons calculé l'incertitude relative sur H pour tous les résultats réalisés lors
de la campagne d'essais en conditions propres, en prenant la valeur de 0 2 % pour l'incertitude
relative sur le débit d'air et la valeur de 0 1 oC pour l'incertitude de température, et nous
obtenons que la valeur maximale pour l'incertitude relative sur H est de 0.82 %

146
ANNEXE-C

INCERTITUDE SUR LA RESISTANCE D'ENCRASSEMENT

Nous avons calculé la résistance d'encrassement à l'aide de l'expression (5.3)

où H t représente le coefficient global d'échange à l'instant t et Ho celui à l'instant initial

Nous calculons l'incertitude sur Rf par :

liR
d
= aH
dR d ôH + dR d IliH
t aH 0
(C 1)
t 0

ou encore:

(C.2)

Application numérique:

Pour les conditions d'une campagne d'essais (essais A6) donnant les conditions intiales .
Tee = 182.6 oC, Tes = 77.6 oC, Tfe = 383°C, Tfs = 39 oC et Ho = 31.6 W/m 2K
et pour une incertitude relative de 0.2 % sur le débit d'air et une incertitude de 0.1 oC sur la
température, on obtient ÔH o = 0.68 %.
Ho
Ainsi, pour leurs conditions de fin d'essais:
Tee = 1872 oC, Tes = 113 5 oC, Tfe = 43.6 oC, Tfs = 442°C et H t = 176 W/m2K
on obtient ÔH t = 0.04 %.
Ht
Enfin, pour ces conditions expérimentales on obtient l'incertitude relative sur Rd est inférieure
à2.38 %.

Nous avons également calculé l'incertitude relative sur Rf pour toutes les campagnes d'essais
réalisées, en prenant la valeur de 0.2 % pour l'incertitude relative sur le débit d'air et la valeur
de 0 1 oC pour l'incertitude de température, et nous trouvons que la valeur de ÔR d ne dépasse
Rd
pas de 8.63 %

147
ANNEXE- D

RECAPITULATIF DES ESSAIS REALISES SUR LE TUBE SONDE

=
Tube - A : da 2.20 Jm, U. :II 2 ml' (Re = 2360), Cpm • 124 ppm
essai Qa Tgm T. t.r ml Cp t.e md
kg!h oC oC h 9 ppm h 91m 2

Al 23.A4 91.54 42.39 4.87 14.51 124.54 5.21 ' 3.86

A2 23A7 87A7 42.29 7.53 22.96 126.10 8.25 5.37
A3 23A9 90.64 42.69 9.91 29.72 127.72 10.67 13.53
AA 23.A4 88.24 42.74 15.38 43.66 119.58 15.68 30.64
A5 23.84 88.14 42.99 24.67 67.82 115.32 24.36 54.11
A6 23.54 90.14 43.76 33A7 100.42 127,45 36.07 71.62
A7 23.54 90.27 43.70 64.68 ' 133.24 87.50 47.85 91.28

=
Tube - B : da 0= 2.20 um Ue 2.8 mIs (Re = 2340) Cpm. 147 ppm
essai Qa Tgm TI t.r ml Cp t.e md
kg!h oC oC h 9 ppm h g/m 2
.
Bl 25.04 179.99 41.65 3.03 13.55 178.48 3.72 8.16
B2 25.04 181.61 38.93 7.35 30.80 167,45 8.45 24.07
B3 25.01 181.37 38.81 14.68 46.55 126.67 12.77 35.21
B4 23.48 182.52 38.06 26.52 65.77 105.70 18.04 63.66
as 25.03 182.39 38.95 24.80 99.72 160.63 . 27.35 81.13
B6 25.06 183.29 40.06 48.10 155.12 128.78 42.55 109.14
B7 24.96 182.07 39.30 120.80 478.37 158.60 131.22 251.66

=
Tube - C : da 3.97 ~, Ue = 2.8 mIl (Re = 2340) Cpm =945 ppm
essai Qa Tgm TI t.r ml Cp I.e md
kg/h oC oC h 9 ppm h 91m 2

Cl 25.04 182.70 40.31 5.17 112.38 868.7 4.89 91.28

C2 25.01 183.91 40.63 6.15 85.33 554.7 3.71 65A9
C3 25.02 184.00 41.19 9.35 157.08 671.6 6.83 161.10
C4 23,45 183.32 43.73 5.00 173.24 1477.0 7.54 194.96
CS 23.50 . 183.40 43.25 13.58 344.37 1053.0 14.98 311.15
C6 24.99 183A9 39.81 17.10 446.31 1044.7 19.42 329.25

148
 ANNEXE-E

RECAPITULATIF DES ESSAIS REAUSES SUR LE FAISCEAU DE TUBES

En régime de Transfert de Chaleur

série - A: da = 2.20 ~m, Um =3.12 mIs, Re == 3140

essais Qa(kg/h) TcerC) Tgm(OC) Ts (OC) dur (h) ml (g) ~m) md (g)

Al 14.10 184.35 136.89 44.52 52.52 102.40 131 .60 3.1310

A2 14.15 185.44 138.19 39.52 79.95 195.88 173.18 . 4.7841
A3 14.11 187.12 ,141.17 43.19 143.95 362.10 178.10 9.3806
A4 14.14 185.20 140.36 '40.49 156.15 363.83 164.51 9.3000
A5 14.17 182.15 136.53 43.87 ' 189.63 413.20 154.50 10.1096
A6 14.20 186.28 142.03 40.51 218.12 " 516.94 167.31 13.2541

=
série - B: da 2.20 ~m, Um =5.6 mIs, Re =5467 .
essais Qa(kg/h) TcerC) Tgm(OC) Ts rC) dur (h) ml (g) ~m) md (g)

81 25.03 185.58 142.21 44.03 23.48 74.68 127.07 2.6904

82 25.02 17872 143.24 42.70 39.13 117.40 119.87 4.5802
83 25.02 178.48 143.24 42.29 41.18 138.37 134.30 5.2948
B4 25.04 186.55 145.17 42.70 96.48 258.57 107.07 10.2497
es 25.00 186.78 148.97 38.04 166.60 681.39 163.59 23.1203

série - C: da .. 3.97 ~m, Um =5.6 mIs,Re .. 5549

essais Qa(kg/h) TcerC) Tgm(OC) TsrC) dur (h) ml (g) ~m) md (g)

Cl 24.97 185.58 141.69 44.04 4.80 101.52 846.72 4.4390

C2 24.95 185.23 141.20 39.11 5.70 152.62 1070.42 8.5953
C3 25.07 186.68 142.56 42.09 8.20 168.59 819.85 6.2551
C4 25.01 187.52 143.72 38.22 8.08 240.67 1189.71 15.2844
C5 25.10 188.61 146.95 39.06 15.72 397.32 1010.81 25.4365
C6 25.12 189.04 147.99 38.01 27.33 616.36 901.99 36.2640
.

=
série - 0: da 3.97 J,m, Um = 9.8 mIs, Re =9439
essais Qa(kg/h) TcerC) Tgm(OC) Ts (OC) dur (h) ml (g) (2pm) md (g)

Dl 43.07 192.16 150.23 43.20 4.40 91 .74 484.13 8.0397

02 43.05 192.88 152.76 43.32 9.52 201.26 491.35 13.1351
03 43.05 195.46 154.91 39.27 11.07 305.45 640.66 24.8277

Qa: débit massique d'aérosol (kg/h)

Tee: température d'entrée du gaz (OC) md : masse du dépôt (g)
Tgm : température du gaz dans la section d 'essai (OC) da : diamètre aéraulique moyen (IJm)
Ts : température de la paroi du tube (OC) Um : vitesse moyen dans le faisceau (mIs)
dur : durée d'essais (h) Re : nombre, Reynolds

149
En ISOTHERME
-

iso - A: da = 2.2O~, Um = 3.12 mIs, Re =2959

essais Qa(kg/h) Tce(DC) Tgm(OC) Ts (OC) dur (h) ml (g) (ppm) md (g)

1 16.39 81.86 81.00 79.95 69.52 132.08 116.05 2.3007

2 16.39 81.54 80.80 79.86 88.87 215.07 147.76 3.8041
3 16.38 82.04 81.35 80.57 120.90 288.95 145.91 4.9683
4 16.36 80.76 79.97 80.18 168.18 477.63 173.50 8.1997

Iso - B : da=2.20 ~m, Um = 5.6 m/s, Re = 5420

essais Qa(kg/h) Tce(OC) Tgm(OC) Ts (DC) dur (h) ml (g) (ppm) md (g)

1 30.10 81.18 81.07 80.93 10.56 42.45 133.55 0.5568

2 30.13 82.42 81.55 266.27 28.68 96.37 111.60 2.2599
3 30.11 81.36 81.23 81.38 31.49 120.61 127.22 2.4192
4 30.11 82.71 82.58 -294.72 48.43 170.00 116.55 5.3903
5 30.11 81.89 81.48 80.07 66.09 235.94 118.59 7.5306
6 29.62 82.21 80.98 77.47 101.40 312.31 102.79 8.7454

lso'-C: da= 3.97 ~m, Um = 5.6 m/s, Re= 5415

essais Qa(kg/h) Tce(OC) Tgm(OC) Ts (OC) dur (h) ml (g) (ppm) md (g)

1 30.13 83.40 83.92 85.31 6.01 101.94 562.54 6.5412

2 30.14 82.71 82.49 80.88 14.50 274.05 627.53 16.9396
3 30.14 82.74 82.86 82.29 22.97 425.46 614.93 24.4907
4 30.12 83.59 83.57 82.50 26.98 439.15 540.24 20.8098

Iso - 0: da=3.97~m, Um=9.8 m/s, Re = 9299

-
essais Qa(kg/h) Tce(OC) Tgm(OC) Ts (oC) dur (h) ml (g) (ppm) md(g)

1 51.43 80.03 80.05 81.11 6.20 88.12 276.62 6.3622

2 51.24 80.98 80.38 79.69 8.77 143.60 319.35 10.5455
3 51.72 79.70 79.51 79.71 24.33 424.71 339.53 26.7163

Qa : débit massique d'aérosol (kg/h)

Tee: température d'entrée du gaz (OC) md : masse du dépôt (g)
Tgm: température du gaz dans la section d'essai caC) da : diamètre aéraulique moyen (~m)
Ts : température de la paroi du tube (OC) Um : vitesse moyen dans le faisceau (mIs)
dur; durée d'essais (h) Re : nombre Reynolds

150
 MASSES COLLECTEES SUR LE FAISCEAU DE TUBES

En régime de Transfert de chaleur

série A: da" 2.20 ~m, Qa. 14 kg/h, Tgm. 140 °CL Ts. 42 oC

essai dur (h) mp (gIs) Cp (ppm) mi (g) md (g) md(g) md (g) -md.tot
Haut Central 8as (g)

Al 52.52 5.2E-Q4 131.6 102.4 0.715 1.460 0.956 3.131

A2 79.95 6.8E-Q4 173.2 195.9 1.060 2.150 1.574 4.784
A3 143.95 7.0E-Q4 178.1 362.1 2.332 4.128 2.920 9.381
A4 156.15 6.5E-Q4 164.5 363.8 1.810 4.642 2.848 9.300
AS 189.63 6.1E-Q4 1S..5 413.2 1.422 5.644 3.044 10.110
A6 218.12 6.6E-04 167.3 516.9 2.238 6.550 4.466 13.254

série 8 : da • 2.20 um, Qa • 25 kath. Tam • 144 oC, Ts .42 oc

81 23.48 8.8E-Q4 127.1 74.7 0.635 1.275 0.780 2.690

82 39.13 8.3E-04 119.9 117.4 0.931 1.410 2.240 4.580
83 41.18 9.3E-04 134.3 138.4 1.130 1.695 2.470 5.295
B4 96.48 7.4E-04 107.' 258.6 1.970 5.175 3.105 10.250
85 166.60 1.1E-03 163.6 681.4 4.780 10.465 7.875 23.120

série C: da" 3.97 ~m, Qa. 25 ka"'. Tom. 144 oC, Ts .. 40 oC

Cl 4.80 5.9E-03 846.7 101.5 0.935 1.860 1.645 4.439

C2 5.70 7.4E-03 1070A 152.6 2.388 2.923 3.284 8.595
C3 8.20 5.7E-03 819.9 168.6 1.530 2.134 2.590 6.255
C4 8.08 8.3E-03 1189.7 240.7 3.128 4.841 7.316 15.284
CS 15.72 7.0E-03 1010.8 397.3 5.245 9.720 10.472 25.437
C6 27.33 6.3E-03 902.0 616.4 8.925 13.494 13.845 36.264

série 0: da. 3.97 ~m, Qa .. 43 kg",. Tom. 152 oC, Ts .. 42 oC

Dl 4.40 5.8E-03 484.' 91.7 1.580 4.030 2.430 8.040

02 9.52 5.9E-03 491 A 201.3 2.845 4.685 5.605 13.135
03 11.07 7.7E-03 640.7 305.4 6.272 8.664 9.892 24.828

dur : durée d'essais

,
mp : débit massique des résidus seches
Cp : concentration en particules
md : masse collectée sur les tubes Haut, Cental, 8as du faisceau
md.tot : masse collectée totale

151
 MASSES COLLECTEES SUR LE FAISCEAU DE TUBES

En ISOTHERME

iso-A: da = 2.20 ~m, Qa. 16 kg/h, T,. 80 oC

essai dur (h) mp Cg/s) Cp (ppm) mi (g) md (g) md (g) md (g) md.tot
Haut Central Bas (g)

1 69.52 5.3E-Q4 116.0 132.1 0.165 1.560 0.575 2.301

2 88.87 6.7E-04 147.8 - 215.1 0.900 1.812 1.092 3.804
3 120.90 6.6E-04 145.9 289.0 1.044 1.608 2.316 4.968
4 168.18 7.9E-04 173.5 477.6 1.380 3.352 3.468 8.200

iso-B: da. 2.20 ~m, Qa a 30 kg/h, T. 81°C

-
1 10.56 - 1.1E-03 133.5 42.5 0.061 0.251 0.246 0.557
2 28.68 9.3E-04 111.6 96.4 0.375 0.925 0.960 2.260
3 31.49 L1E-03 127.2 120.6 OA35 0.905 1.080 2.419
4 48.43 9.7E-04 116.6 170.0 0.355 2.640 2.395 5.390
5 66.09 9.9E-04 118.6 235.9 1.245 3.420 2.865 7.531

iso-C: da = 3.97 ~m. Qa. 30 kg/h, T,. 82 oC

1 6.01 4.7E-03 562.5 101.9 0.995 2.405 3.141 6.541

2 14.50 5.3E-03 627.5 274.1 1.880 7.820 7.240 16.940
3 22.97 5.1E-03 614.9 425.5 0.740 7.850 15.900 24.491
4 26.98 4.5E-03 540.2 439.1 1.145 6.990 12.675 20.810

Iso-O: da,. 3.97 ~m. Qa. 51 kg/h, T. 80 oC

1 6.20 4.0E-03 276.6 88.1 1.483 3.221 1.658 6.362

2 8.77 4.5E-03 319.3 143.6 2.450 5.356 2.740 10.546
3 24.33 4.9E-03 339.5 424.7 4.088 14.469 8.160 26.716

dur : durée- d: essais

mp : débit massique des résidus seches
Cp : concentration en particules
md : masse collectée sur les tubes Haut, Cental, Bas du faisceau
md.tot : masse collectée totale

152
Essai A6 en tranat1 . de dWeur : Tsm "142 oc. Ts '" 4l"C
da - 2.2 ~m.. Qa-' 14.2 kgIb. Cp - L67..3 ppm. dur"' 218 h.

<=: Sens de l'écouklment

Essai ISOA4 en bolhcrme , T - gO OC

da " '.2.2j1m, QI - 16.4

153
lissai B4 en ItaIlsfert de chaleur : Tgm = 145 DC. Ts = 43 oC
da. " 2.2 J.Lm~ Qa " 25 kg/b. Cp = 107 ppm, dUT = 96_S h_

{= Sens de l'éoo ulemeot l

Essai ISOBS en Isotherme a T '"' g 1 oC

da - 2_2 Ilm. QI! ~ 30 \cglh, Cp - 118.6 ppm, dur = 80 h.

$= Sens de l'écoulement

154
Essai CS enlransfel't de chaleur : 'rgm " 147 · C, Ts " 40 · C
da .. 3.97 ~m, Qa " 25 ~ Cp - 10llppm, dur - 15.1 h.

~ Se.ns de l'éooulemer"lt

Essai ISOCJ en Isotherme il T = 8.2 OC

da ~ 3.97

~ Sans de l'écoulement

15S
 • 0')en tmufm de ~ : Tarn" 154 oc. TI - 4,f oc
da - 3.97 )lm. QI - 43 k&'h. Cp" 640 ppn. dur'" lU h"

~ Sens de J'écoulement

EsAi 18003 en Isotherme à T .. 80 oC

da - 3.97 "Ill, Qa - SI kglh, Cp ... l40 ppm, dur - a4.3 h.

<= Sens dé l'écoulement

IS6
ANNEXE-F

PROPRIETES THERMOPHYSIQUES DE L'AIR SEC [61]

1. Masse volumique

La masse volumique de l'air est donnée en fonction de la pression et de la température par la

formule:

p =348.29( ~) (~) (F.1)

avec p est la masse volumique en kg/m3, P la pression absolue en bar, T la température

absolue en K et z un facteur de compressibilité fonction de P et T donné par

z =1-1O-7 p(T -350)2 pour 180 ST S 350K

z=l pour 350 < T S 1300K

2. Chaleur massique à pression constante

La chaleur massique de l'air à pression constante en (kJlkg.K) est donnée pour T variant de
180 K à 450 K par les formules suivantes:

Pour P = 1 bar'
C p =AT3 +BT2 +CT+D (F.2)
avec.
A = -0.3798 10-8 B = 0.4178 10-5 C = -0.1456 10-2 D = 1.1688

Pour P = 5 bar.
3 2
Cp = AT +BT +CT+D (F.3)

avec:
A = -0.645 10-8 B = 0.7102 10-5 C = -0.2536 10-2 D = 1. 307

Pour P = 10 bar
2
Cp =AT 3 +BT +CT+D (F.4)
avec:
A = -0.97106 10-8 B = 0.1075 10-4 C = -0.3912 10-2 D = 1.487

157
3. Conductivité thermique

La conductivité thermique de l'air en (kW/m.K) pour T variant de 180 K à 600 K est donnée
par les formules suivantes

PourP= 1 bar
(F.S)
avec
A =-3.38110-8 B = 9.814 10-5 C =-1.308 10-4

Pour P = 10 bar
(F 6)

avec'
A =-2.956 10-8 B =9.349 10-5 C =1. 395 10-3

4. Viscosité dynamique

La viscosité dynamique de l'air en (pa s) pour T variant de 180 K à 600 K est donnée par .

Pour P = 1 bar:
(F.7)
avec'
A =-2.858 IO-Il B =6.605 10-8 C =1.20 1O...{i

Pour P = 10 bar:
(F.8)

avec:
A =-2.738 10-11 B =6.474 10-8 C =1.620 1O...{i

158
ANNEXE-G

PROPRIETES PHYSIQUES ASSOCIEES AVEC DES AEROSOLS

Nous avons w au paragraphe 2.2.1.1 que le transport par diffusion thennophorètique est
fonction du coefficient adimensionnel aTLe (rapport entre la diffusivité thennophorètique des
particules et la diffusivité thermique du gaz porteur) donné par :

(G.l)

Nous avons utilisé la fonnule de Talbot et al. (cité dans [11]) pour estimer le tenne ntLe qui
dépend de la taille des particules, du rapport des conductivités thermiques du gaz et de la
particule ainsi des interactions gaz-particule (cf expression 247) .

(G.2)

avec AI =1.2 ,A2 =041 ,A3 =0.88 ,et cs =1.17 Ct =2.18 cm =l.14, Kn est le nombre
de Knudsen, kg la conductivité thennique du gaz porteur, et kp la conductivité thermique de la
particule.

Pour un aérosol de forte conductivité thennique (cité dans [9]) :

k
NaCI dans l'air - g = 0 04
kp
k
l'uranine dans l'air ~ =0 065
kp
k
Pour des particules de faible conductivité thennique (goutte d'huile) : ~ = 0 1
kp
On peut calculer le coefficient de a,Le pour différentes valeurs de kp et de Kn à l'aide de
l'expression 4.11 dont les résultats sont donnés sur la figure G-1

159
 0.5
0.45 1.gaz = 180 oC, kg = 0.032 Wjm.K

04 Kn = 1
0.35
0.3

-
....CD 0.25
0
0.2
Kn = 0.1

015
0.1 Kn = 0.01
005
0
0 0.5 1.5 2 2.5
kp (Wjm.K)

Figure G-l Variation de ŒtLe en fonction du kp'

D'autre part, pour les particules de silice (Si0 2 )dans l'air, Rosner et Park [11] proposent
l'expression suivante:

070
k = 3.0 X 10-3 --.!... cali s (G.3)
p ( 300 ) cmK

A titre d'exemple, pour la température moyenne de l'ordre de 180 oC, on a .

kp = 1938 W /m.K

Par ailleurs, Rosner et a/.[4] ont également souligné que l'on peut utiliser l'expression (G 3)
pour estimer la conductivité thermique des particules de Sulfate de sodium et de Calcium dans
l'air.

160
 ŒO
CENG/SERVICES TECHNIQUES/REPROGRAPHIE/1995