Aeraulique

Télécharger au format pdf ou txt
Télécharger au format pdf ou txt
Vous êtes sur la page 1sur 53

ARAULIQUE

Mai 2000

Pierre NEVEU

1re PARTIE :

ECOULEMENT EN CONDUITE

1 EQUATION DE BERNOULLI GENERALISEE. ................................................................................................ 4


1.1 coulement isotherme, pas d'changeur.......................................................................................... 5
1.2 Prsence dun changeur ou temprature non-uniforme ................................................................. 7
2 FLUIDE PARFAIT , APPLICATION A LA MESURE DE DEBIT ........................................................................... 8
2.1 Tube de Pitot. ................................................................................................................................ 8
2.2 Phnomne de Venturi ................................................................................................................. 11
3 ECOULEMENT EN CONDUITE DES FLUIDES REELS (NEWTONIENS), PERTES DE CHARGES ............................ 13
3.1 Rgime tabli............................................................................................................................... 13
3.2 Notion de pertes de charge........................................................................................................... 14
3.3 Pertes de charge rgulires ......................................................................................................... 14
3.4 Conduites non cylindriques. ......................................................................................................... 18
3.5 Pertes de charge singulires. ....................................................................................................... 19
QUELQUES OUVRAGES UTILES. ............................................................................................................... 20

L'hydraulique s'intresse aux relations relatives au mouvement des fluides incompressibles (quivolumiques), dans une veine d'coulement, que ce soit une conduite (circulaire ou non), un changeur de chaleur
ou une pompe. Toutefois, les relations reprises ici s'appliqueront dans la plupart des cas l'araulique,
puisqu'alors, l'hypothse, bien qu'rrone, de fluide incompressible, est en gnral accepte pour le calcul des
rseaux de traitement d'air.
Les relations de base sont celles issues de la mcanique des fluides, qui, partir des relations de conservation
locales (masse, impulsion (ou quantit de mouvement), et nergie) permettent d'exprimer sous forme d'quations
au drives partielles les relations entre pressions, vitesses et tempratures (cf. poly. mcanique des fluides). Ce
sont les clbres relations :

d'EULER pour l'coulement d'un fluide parfait (non visqueux)

de NAVIER-STOCKES pour l'coulement des fluides newtowniens.


La rsolution de ces quations reste trs difficile, et n'admet des solutions analytiques que pour quelques cas
d'cole trs simples. Pendant longtemps, seules les quipes de recherche les ont utilises, pour produire
notamment diverses abaques ou corrlations donnant directement les caractristiques macroscopiques des
coulements (pertes de charge, coefficient d'change) directement utilisables par l'ingnieur. Soulignons
cependant un regain d'intrt de cette formulation locale due au dveloppement la fois du matriel informatique
et des mthodes de rsolution numrique. En effet, le dimensionnement optimal d'un composant quelconque
sige d'un coulement, ncessite la prise en compte du couplage entre les phnomnes thermiques et
hydrauliques dont seules ces quations locales permettent la simulation. Les efforts de recherche entrepris tout
d'abord en mtorologie, puis pour l'optimisation des chambres de combustion permettent aujourd'hui d'avoir
notre disposition des mthodes numriques de rsolution, la fois fiables et puissantes, pour le dimensionnement
d'objets plus simples que sont les changeurs et les conduites. Outre les outils "maison" dvelopps par les
manufacturiers d'changeurs, plusieurs logiciels sont actuellement commercialiss pour l'aide au
dimensionnement d'changeurs (par exemple CETUC, CANUT ) et plus spcifiquement pour la caractrisation
des coulements (par exemple FIDAP, FLUINT).
Si on se place maintenant au niveau d'un avant-projet, il est clair que les quations locales ne sont d'aucun
secours, puisque leur rsolution numrique impose de dfinir, et donc de connaitre, le domaine d'intgration
(conditions aux limites), c'est dire les dimensions gomtriques du composant rel que l'on veut simuler, et qui
sont inconnues a priori. Il apparat donc ici deux niveaux d'tude dans la conduite d'un projet:

tout d'abord, le pr-dimensionnement de toute installation doit permettre d'estimer les dimensions
gomtriques de tous les composants. On utilise ici des quations macroscopiques, en ne s'intressant qu'aux
grandeurs d'entre/sortie de chaque composant (dbit, pression, temprature).

ce pr-dimensionnement tant effectu, les quations locales (ou plutt les logiciels permettant leur
rsolution) permettent alors d'affiner et d'optimiser chaque composant.
C'est au premier niveau d'tude que l'on s'intresse ici. D'ailleurs, nombreux sont les bureaux d'tude qui s'y
arrtent, le deuxime niveau restant (malheureusement!) rserv aux diffrents fabricants des composants
intervenant dans le rseau (pompes, changeurs, vannes, etc.)
Nous rappelons ici les principales quations utiliss en hydraulique (quation de Bernoulli) et pour l'tude des
turbo-machines.

quation de Bernoulli gnralise.


La forme gnrale d'un rseau peut prsenter des changements de section (et donc de vitesse d'coulement),
d'altitude et de direction. Il comprend en outre au moins un organe moteur (pompe,ventilateur), et
ventuellement un ou plusieurs changeurs (figure 1).
En appliquant les deux principes de la thermodynamique sur , on obtient:
dE t & &

V 2 VS2
& ( h E h S ) + E
= Q+W+m
+ g ( z E z S )

dt
2

&
dS
Q
= +m
& (s E sS ) + P(S)
dt T
En rgime permanent, [1] conduit :

[1]

V 2 VS2
& +W
& +m
& ( h E h S ) + E
+ g ( z E z S )
0 = Q

&
0 = Q + m
& (s E s S ) + P(S)

[2]

D'autre part, pour un corps pur, les variations lmentaires d'enthalpie et d'entropie peuvent tre mise sous la
forme (cf. poly. thermodynamique):

v
dh = c p dT + v T dP
T P

ds = c dT v dP
p

T T P

[3]

Dans le cas qui nous intresse ici (fluide incompressible), on a:


v
= Cste, c'est dire v = 1 = Cste , d'o = 0 . Les quations [3] se simplifient alors en :
T P
dh = c p dT + v dP

ds = c dT
p

[4]

que l'on peut introduire dans [2] aprs intgration de TE TS et de PE PS:

V 2 VS2
& +W
& +m
& c p TE / S + v PE / S + E
+ g ( z E z S )
0 = Q
2

&
T
Q
E
0 = + m
& c p Ln
+ P(S)

TS
T

[5]

.
Q

.
W

ZS

.
m
machine

changeur

ZE

Figure 1
coulement isotherme, pas d'changeur
Dans ce cas, TE = TS, et [5] conduit :

V 2 VS2
& +W
& +m
& v PE / S + E
+ g ( z E z S )
0 = Q
2

Q
& = T P(S)

[6]

Il faut souligner ici que l'coulement isotherme implique un change de chaleur avec l'extrieur (Cf. 2nd
principe). En fait, l'exprience montre qu'un fluide circulant dans une boucle ferme s'chauffe mme en
l'absence d'un dispositif de chauffage. Ce sont les frottements internes (T. P(S)) qui impliquent cet chauffement.
La temprature se stabilise une temprature Tfluide > Tambiant permettant l'vacuation de cette chaleur cre (non
compense) vers l'extrieur. On regroupe en gnral le 1er et le 2nd principe en une seule quation en faisant
apparatre la masse volumique du fluide ( = 1/v),:
p VS2 VE2
&
&
m
+
+ g (z S z E ) = W
T. P(S)

&
W

[7]

: Puissance mcanique change avec lexterieur

T. P(S)

: Puissance dissipe au sein du fluide (frottements)

& peut tre soit positif (ventilateur, pompe), soit ngatif (turbine). Par contre, T.
Dans le cas le plus gnral, W
P(S) est toujours positif . La forme habituelle de l'quation de BERNOULLI gnralise s'crit en introduisant le
& :
travail des forces de frottement Tf (toujours ngatif), et en divisant par le dbit m
p S p E VS2 VE2

+
+ g (z S z E ) = w T. i s = w + Tf

[8]

L'quation [8] exprime que la variation d'enthalpie totale d'une masse unitaire de fluide (terme de gauche) est
gale au travail mcanique w chang avec le milieu extrieur diminu du travail des forces de frottement. Deux
formes sont plus particulirement utilises en hydraulique, la premire exprime en pression (Pa), la seconde
exprime en hauteur manomtrique (m).
La premire forme s'obtient en multipliant [8] par :
pS p E +

VS2 VE2
+ g (z S z E ) = W + Tf
2

[9]

soit, en regroupant les termes relatifs l'entre et la sortie du systme:


V2
V2
pS + S + g z S p E + E + g z E = Pm Pf
2
2


On nomme :

V2
Pt = p +
+ g z

P* = p + gz
V2
Pd =
2
V2
Pa = p +
2
Pm
Pf = Tf

[10]

: Pression totale
: Pression motrice
: Pression dynamique
: Pression darrt
: Diffrence de pression totale induite par la machine
: Perte de charge exprime en Pascal

La deuxime forme s'obtient en divisant [8] par = .g (poids spcifique),


p S p E VS2 VE2
w T
+
+ (z S z E ) = + f

2g
g g

[11]

qui amne:
p S VS2
p

V2
+
+ z S E + E + z E = H m H f

g
2
g
2

On nomme :
p V2

H = +
+ z
2g

[12]

: Charge hydraulique (ou simplement charge)

H* = + z

: Hauteur piezomtrique

Hm (ou Hm)
Hf

: Hauteur manomtrique de la machine ( > 0 ou < 0 )


: Perte de charge ( > 0 ), exprime en mtre

Remarques:
a) Dans le cas o il n'y a pas de pertes de charge (Hf ou Pf = 0), on retrouve l'quation de Bernoulli classique,
valide dans le cas des fluides parfaits (non visqueux).
b) Dans le cas o il n'y a ni pertes de charge, ni change d'nergie mcanique, on retrouve l'quation de
l'hydrostatique (fluide au repos).

c) l'quation [8] peut tre dduite des quations de conservation locales (en particulier de Navier-Stockes).
Prsence dun changeur ou temprature non-uniforme
Dans ce cas, les quations prcdentes ne sont plus valables. En toute rigueur, il faudrait ici utiliser
ncessairement les quations locales. En fait, on ne le fait jamais pour un avant projet. On simplifie le problme
en dcoupant le rseau en tronons isothermes et lorsque que ce n'est pas possible (changeurs), en dcouplant
les phnomnes thermiques et hydrauliques. Par exemple, pour l'changeur du rseau reprsent figure 1, on
remplacerait la transformation relle 12 par:
1) un chauffage isobare 12':
& +m
0 = Q
& c p (T1 T2' )
12 '

&
Q 12'
T

& c p Ln( 1 ) + P(S)12'


0 = T + m
T
2'

& =m
Q
& h12'
12'
&

Q
12'
& s12'
P(S)12' = T + m

2) une dtente isotherme 2'2 du fluide incompressible de masse volumique moyenne : m =

& +m
& + p
& ( h 2 h 2' ) = Q
0 = Q
2' 2
2' 2
m

& = T P(S)
Q
2' 2
2' 2

[13]
1 + 2'
2
[14]

En intgrant [13 et 14] dans le bilan global suivant:

V 2 V22
& +m
& (h 1 h 2 ) + 1
+ g (z 1 z 2 )
0 = Q
2

&
Q
0 = + m
& (s1 s 2 ) + P(S)

[15]

on obtient:

V 2 V22
& +Q
& +m
& (h 1 h 2' + h 2' h 2 ) + 1
+ g ( z E z S )
0 = Q
12 '
2' 2
2

&
&
0 = Q 12' + Q 2' 2 + m
& (s1 s 2' ) + P(S)12' + P(S) 2' 2

[16]

soit finalement

p 2' p 2 V12 V22


& +m
&

+
+ g ( z 1 z 2 )
0 = Q
22'

&
0 = Q 22' + P(S)
22'

[17]

& = 0 et dans laquelle la masse volumique prise en compte est


On retrouve un expression similaire [8] avec W
une masse volumique moyenne : Les trois formes de l'quation de Bernoulli gnralise sont similaires celles
trouves prcdemment (eqs [10, 11 et 12]):
En nergie :
p1 V12
p

V2

+
+ g z S 2 + 2 + g z 2 = T. i S 22 ' = Tf

2
2
m

En pression :


V2
V2
pS + m S + m g z S p E + m E + m g z E = Pf
2
2


En hauteur :

pS VS2
p

V2
+
+ z S E + E + z E = H f

m 2g
m 2g

La figure 2 prsente la variation des diffrents termes de l'quation [2] le long d'une canalisation. Pour un fluide
parfait (non visqueux), le terme pf disparat et l'on retrouve l'quation classique de Bernoulli. Cette dernire
permet de mettre en vidence deux phnomnes utiliss pratiquement pour la mesure de dbit.

tuyre de Venturi,

Tube de Pitot.
L'tude de ces deux dispositifs sera expos daans une premire partie, afin de mettre en vidence la diffrence
entre pression statique, pression dynamique et pression totale.
Les autres paragraphes concernent le dimensionnement des rseaux. Concernant des fluides rels (visqueux), ce
dimensionnement ncessite le calcul des pertes de charge Pf. En effet, les conditions d'entre et de sortie tant
en gnral fixes par le cahier des charges (cf. fig.1), le calcul des pertes de charge permet de dterminer la
pression totale de la machine installer dans le rseau (cf equations [2],[3]).

Sortie
pf

Entre

pm

Pression atmosphrique

pression totale

pression statique
pression dynamique

Figure 2 : volution des diffrents termes de [10] lors de la traverse du circuit


Fluide parfait , application la mesure de dbit
Tube de Pitot.
Principe.
On considre le montage reprsent sur la figure 3, constitu de deux tubes TA et TB.
B1
ZA2

A1
ZA1
TA

TB

u
A2

B2

Figure 3 : mesures de la pression statique (TA) et de la pression darrt (TB).


En appliquant lquation classique de Bernoulli aux deux tubes TA et TB, et en remarquant que :
u A 2 = u A1 = u B1 = 0 , que p A1 = p B1 = p o (pression atmosphrique), et que p A 2 = p B2 = p (pression dans la
veine fluide dite pression statique) il vient :

pour le tube A : z A2 +

p
p
= z A1 + o
g
g

(18)

pour le tube B : z B2 +

u 2B
p
p
+ 2 = z B1 + o
g 2g
g

(19)

En posant z A 2 = z B2 = 0 (rfrence dans laxe de la conduite), on obtient une mesure de la pression statique p
par le tube TA et une mesure de la pression darrt pa par le tube TB:
p = p o + g z A1

u2
= p o + g z B1
pa = p +
2

Dautre part, en soustrayant les deux quations (18) et (19), on obtient :

p A u B2
pB
z B2 z A2 + 2 2 + 2 = z B1 z A1
g g 2g

Or, z B2 z A 2 = 0 et p B2 p A2 = 0 , do finalement :
z B1 z A1 =

u 2B1

(20)

2g

ce qui permet dexprimer la pression dynamique en fonction de z B1 et z A1 :


u2
= g z B1 z A1
pd =
2

(21)

La relation (21) permet de remonter au dbit volumique. Dans le cas o la vitesse dcoulement u est constante
dans toute la section S de la conduite, on a :

& = u n ds = u S = S 2g z z
V
B1
A1

(22)

Dispositif pratique.
En pratique, on modifie le montage de la figure 3 en utilisant des manomtres tube en U, et ceci pour deux
raisons :

le reprage de z nest pas facile par rapport un repre fictif (axe de la conduite),

la hauteur des colonnes de fluide impliques par un tel montage seraient trop leves : dans le cas de leau
par exemple, une pression de 2 bars conduit :
z=

p po
105
=
= 10,2 m !!
g
1000 9,81

On utilise donc un manomtre tube en U, rempli dun fluide auxiliaire dont la masse volumique est plus leve
que celle du fluide circulant dans la conduite (figure 4). En gnral, le mercure (Hg = 13600 kg/m3) est utilis
sur les conduites de liquide, et leau sur des conduites de gaz. Mais tout dpend en fait de la valeur des pressions
statiques p et darrt pa et de la compatibilit des deux fluides.
On peut crire lquation de Bernoulli pour les deux tubes TA et TB en distinguant le fluide circulant dans la
conduite (indice fl) et le fluide auxiliaire (indice aux) du manomtre :

A2

B2

u
TA

TB

B0
A0
zB

zA

A1
B1

Figure 4 : Mesures des pressions statiques et darrt, dispositif avec tubes en U.

p A1

p
= zA1 +
z A 2 +
fl. g
fl g

pour le tube A :
p
p0
A1
z +
= zA 0 +
A1 aux . g
aux g

(23)

p B1
u 2B
p
z B2 +
+ 2 = z B1 +
fl g 2 g
fl g

pour le tube B :
p B1
p0

z B1 + g = z B0 + g
aux .
aux

(24)

En combinant les quations de ces deux systmes, on obtient :


z A = z A 0 z A1 =

z B = z B0 z B1 =

p A1 p o
aux . g

p B1 p 0
aux g

(p +
p+

fl

))

g z A 2 z A1 p o
aux . g

fl u 2B2
2

(25)

+ fl g z B2 z B1 p 0

(26)

aux g

que lon peut crire plus simplement :


z A =

p po
+ fl z A 2 z A1
aux. g aux.

(27)

u2
p + fl B2 p o

z B =
+ fl z B2 z B1
aux g
aux

(28)

En rgle gnrale, les termes fl aux ( z 2 z1 ) restent ngligeables devant les z. Le dispositif de la figure 4
permet donc de mesurer la pression statique (tube A) et la pression darrt (tube B) dans la veine fluide :
p aux g z A et p a aux g z B . Notons que ces deux dernires relations donnent des pressions relatives

et ngligent la colonne de fluide remplissant partiellement le tube z A2 z A1 et z B2 z B1 .


Les relations (27) et (28) permettent dobtenir immdiatement une expression de la pression dynamique :

10

p d = p a p = aux . g ( z B z A ) fl g z A1 z B1

(29)

L encore, le dernier terme est gnralement ngligeable. Lexpression de la vitesse u devient alors :
u=

2 pd

2 aux g ( z B z A )
fl
fl

(30)

La figure 5 montre un dispositif appel tube de Pitot. On pourra vrifier que cet appareil mesure la pression
dynamique, et donc la vitesse u. Il est ainsi utilis pour mesurer des dbits volumiques. Les expressions de la
pression dynamique et de la vitesse peuvent tre dduites de la mme manire que prcdemment. On trouve :


p d = p a p = aux . g z 1 fl
aux.

u=

(31)

2 pd
= 2 g z aux 1
fl
fl

(32)

Figure 5 : Tube de Pitot


Phnomne de Venturi
coulement dun fluide parfait dans un convergent/divergent.
On sintresse ici lcoulement dun fluide incompressible parfait au travers dun dispositif comprenant un
convergent et un divergent (figure 6). En appliquant lquation de Bernoulli entre les sections S1 et S2, on
obtient:
p1 +

u12
u2
+ z1 = p 2 + 2 + z 2
2
2

(33)

Or, z1 = z2 do :
p 2 = p1 +

u12 u 22
2

(34)

& = u1 S1 = u 2 S 2 , il vient :
En remarquant que V
p 2 = p1 +

S2
u12
(1 12 )
2
S2

(35)

Lquation (35) illustre le fait que, la section de passage diminuant, la vitesse du fluide va augmenter : il y a
conversion de pression statique en pression dynamique, cest dire conversion dnergie potentielle en nergie
cintique (On se rappellera que la pression est homogne une nergie par unit de volume). En effet, S2 tant
infrieur S1, (35) montre que la pression p2 sera infrieure la pression p1 (voir figure 6). Entre les sections S2
et S3, la mme analyse conduit une relation analogue :
p3 = p2 +

u 22
S2
u 23 S 23
(1 32 ) = p 2 +
( 2 1)
2
2
S2
S2

(36)

11

qui montre que dans le divergent, on assiste la conversion inverse : il y a augmentation de la pression
statique et diminution de la pression dynamique.
Pression totale pt (P 0)
pd =

u 2
2

u1

u2

S2
S1

Col

Convergent

Divergent

Figure 6 : ensemble convergent-divergent. volution de la pression statique p et de la pression dynamique pd, en


supposant la pression totale constante (fluide parfait).
Applications
Trompe eau, jecteurs.
Les quations (35) et (36) montrent que la pression au col est infrieure la pression dentre du convergent et
de sortie du divergent. Un tel dispositif peut donc produire une dpression. Cette proprit est utilise dans les
trombes eau et les jecteurs pour pomper un autre fluide.
Eau
(fluide moteur )

Col
P < Po

Po

Figure 7 : Principe dune trompe eau.


Dans une trompe eau, la pression de sortie est gale la pression atmosphrique (figure 6). En consquence, la
pression au col sera infrieure cette dernire. Les trompes eau sont ainsi utilises pour vidanger des
rcipients, ou pour raliser un pr-vide (Pmin 50/100 mbars) qui facilitera le travail ultrieur dune pompe
vide.
Tuyre de Venturi.
Lquation (35) permet dexprimer la vitesse u1 (et donc de remonter au dbit volumique), en fonction de la
diffrence de pression statique entre lentre et le col :
u1 =

2 ( p 2 p1 )

(37)

S2
1 12
S2

En munissant lentre du convergent et le col dun manomtre tube en U (figure 8), on obtient une mesure
& = u S . En effet, on a
directe du p = p1- p2, et par suite une mesure indirecte de la vitesse u1 et du dbit V
1
1
(raisonnement analogue lobtention de (31))

12



p = p1 p 2 = aux. g z 1 fl
aux .

(38)

En combinant (37) et (38), il vient :


2 aux . g z 1 fl
aux .

2
=
aux. 1 g z
u1 =
2
2
S

S1
fl

2 1
fl 12 1
S2
S2

(39)

S2
En posant K = 2 12 1 , qui est une caractristique de lensemble convergent-divergent, (39) se rduit :
S2

u1 = K aux . 1 g z
fl

(40)

Un tel dispositif, reprsent figure 8 est appel tuyre de Venturi ou simplement Venturi.

u1

Figure 8 : Venturi
coulement en conduite des fluides rels (newtoniens), pertes de charges
Rgime tabli
Tout ce qui va suivre nest valable quen rgime tabli, cest dire lorsque le profil de vitesse ne dpend plus de
la distance lentre de la conduite. Cette notion est trs importante, puisqu chaque fois que le fluide traverse
un composant particulier (vanne, changeur, rtrcissement, largissement ....), il y a perturbation de ce profil et
le rtablissement de celui-ci, qui ncessite une longueur de lordre dune cinquantaine de diamtres, va
engendrer une perte de charge supplmentaire, que lon appelle perte de charge singulire. Par opposition, la
perte de charge due la viscosit du fluide est appele perte de charge rgulire ou linaire (on verra quelle est
directement proportionnelle la longueur de conduite).
Cette longueur dtablissement correspond en fait au dveloppement de la couche limite lintrieur de la
conduite. Deux cas sont ici envisager (figure 9) :

Re < Rec, le rgime est laminaire . A lentre de la conduite la couche limite laminaire commence se
dvelopper, puis partir dune distance L 50 D, envahit toute la conduite. On retrouve le profil de vitesse
parabolique.

Re < Rec, le rgime est turbulent. La vitesse du fluide est telle que la couche limite ne peut pas envahir toute
la section de la conduite. De plus, il y a dveloppement d'une couche limite turbulente en T, puis
tablissement de celle-ci lorsque son paisseur atteint une valeur proche de 0,2D. La rgion centrale
prsente alors un profil presque plat (dissipation visqueuse ngligeable).

13

Re < Rec

ue

u = u(r) pour r <


u = ue pour r >

u = u(r)

sous couche laminaire


couche limite turbulente

Re < Rec
ue

u = u(r) pour r <


u cste pour r >

Figure 9 : Le rgime est tabli lorsque le profil de vitesse ne varie plus selon laxe de la conduite. Les couches
limites se sont alors compltement dveloppes.
Notion de pertes de charge.
L'application de l'quations (1) un tronon de conduite ne comportant pas de machine se simplifie en:
p p E u S2 u E2

w f = T. i s = S
+
+ g ( z S z E )
2

(41)

Rappelons que le terme wf correspond au travail des forces de frottement qui s'oppose au mouvement du fluide
(<0). Ce travail est totalement dgrad en chaleur, qui est vacue par les parois de la conduite. En consquence,
la diffrence d'nergie totale entre l'entre et la sortie (2me membre de (41)) est obligatoirement ngative. Il en
est donc de mme pour les diffrences de pression totale et de charge, au travers de (2) et (3):

u 2 u 2E
p tS p tE = p S p E + S
+ g (z S z E ) < 0
2

(42)

p p E u 2S u 2E

HS HE = S
+
+ (z S z E ) < 0
2g

(43)

On appelle perte de charge loppos de cette dernire diffrence : H = H E H S , qui est par consquent
toujours positif, et qui correspond bien une perte (une diminution) de la charge dans le sens de lcoulement.
Par extension, on appelle galement perte de charge son quivalent en pression g H = p t (equ 42), qui
correspond une chute de pression totale :
Pertes de charge rgulires .
On peut montrer par lanalyse dimensionnelle que cette chute de pression peut scrire sous la forme:
p t ( Re, r ) u 2
=

L
D
2

(44)

p t L
: perte de charge linaire (Pa/m),

: coefficient de pertes de charge,


Re : nombre de Reynolds,

: masse volumique du fluide (kg/m3),


r :rugosit rduite de la conduite.
u
: vitesse moyenne (m/s).
La vitesse u prendre en compte est parfois appele vitesse moyenne de dbit. Elle correspond au rapport du
& par la section de passage S : u = V
& S.
dbit volumique V

avec

14

La longueur caractristique entrant dans le calcul du nombre de Reynolds est le diamtre D pour une conduite
cylindrique.
La rugosit relative r est le rapport entre la rugosit absolue de la conduite et son diamtre : r = D . La
rugosit absolue reprsente la hauteur des asprits prsentes dans toutes les conduites industrielles. Le tableau
I donne quelques valeurs de qui dpendent essentiellement du matriau en contact avec le fluide.
Le coefficient de perte de charge peut tre dduit, pour les problmes simples, de la rsolution des quations de
Navier-Stockes, mais est maintenant donn directement sous forme de corrlations.
Remarque : Notons ici que certains ouvrage (notamment anglo-saxons), utilisent non pas le coefficient de pertes
de charge, mais le coefficient de frottement (Friction factor) qui relie la contrainte tangentielle de frottement
la paroi o la pression dynamique :
o = C f ( Re, r )

u 2
2

(45)

Contrainte de frottement et perte de charge tant deux aspects dun mmes problme, on peut exprimer lun en
fonction de lautre. On montre que :
o =

D p

4 L

(46)

et on dduit de (44), (45) et (46) que :


= 4 Cf

(47)

p
u2
= 2 C f ( Re, r )
L
D

(48)

Tableau I : Rugosit relative de quelques matriaux


Matriau
Qualit
rugosit absolue e(mm)
Acier Inox
lisse
0,03
PVC
Aluminium
Acier galvanis joint spiral
moyennement lisse
0,09
Acier galvanis joint longitudinal
0,15
Fibre de verre
moyennement rugueux
0,9
Flexible
rugueux
3
Avant dutiliser une corrlation, on sassurera donc du coefficient considr dans louvrage utilis, dautant plus
que les appellations anglo-saxonnes sont trs voisines (Loss friction factor pour , et friction factor pour Cf)
et que ces deux coefficients sans dimension sont souvent dsigns par la mme lettre f. Il est donc
indispensable de reprer lexpression des pertes de charge donne souvent la premire page et de la comparer
aux expressions (44) et (48) pour pouvoir conclure.
Rgime laminaire.
Les quations de Navier-Stockes fournissent immdiatement les lments caractristiques de lcoulement
savoir :

un profil de vitesse parabolique,


4
& = R p ).

un dbit volumique proportionnel la perte de charge (formule de Poiseuille: V


8 L
En faisant apparatre la vitesse u et le diamtre D dans l'quation de Poiseuille, on obtient :
& = u S = u D =
V
4
2

( )

D2
8

p

L

do,
p
32 64 u 2
= u 2 =

L
UD 2D
D

(49)

15

soit :
p 64 u 2
=

L
Re 2D

(50)

qui montre, en comparant avec (44), que = 64 Re


Rgime turbulent.
Remarques prliminaires
Pour les coulement en conduite, le rgime laminaire ne se maintient pas (sauf prcautions particulires) des
nombres de Reynolds suprieurs 2000/2300. Ceci correspond par exemple pour une conduite de diamtre D
=10 cm, une vitesse moyenne dcoulement de 2 cm/s pour leau ( = 1000 kg/m3, =10-3 Pl) et de 26 cm/s
pour de lair ( = 1,29 kg/m3, =1,7 10-5 Pl). Une telle vitesse ne correspond pas la pratique industrielle o les
vitesses sont de lordre de 1 2 m/s pour les liquides et de 5 10 m/s pour les gaz.
On peut galement noter que si le rgime laminaire se maintenait grande vitesse, les coulements naturels dans
les fleuves et les rivires auraient des vitesses normes. Par exemple, la rsolution des quations de NavierStockes pour un coulement laminaire en nappe appliquer au Rhne Lyon, o la pente moyenne du lit est de 6
10-4 , conduirait des vitesses suprieures 1000 m/s ds que la hauteur deau dpasserait 1 m ! , ce qui est
absurde. En ralit, la vitesse ne dpasse pas 2 m/s. Cette exemple simplissime illustre le caractre fortement
dissipatif de la turbulence, notamment dans la couche limite dont ltablissement ralentit normment
lcoulement.
Coefficients de pertes de charge dans les conduites industrielles.
Plusieurs travaux tant thoriques (modlisation de la couche limite) quexprimentaux (Blasius, Karman,
Prandtl, Nikuradz en sont les auteurs les plus connus) ont permis dtablir des abaques donnant en fonction du
nombre de Reynolds. On trouvera ci-aprs (figure 10) labaque le plus utilis, labaque de Moody, qui est
paramtr en rugosit relative r. Colbrook a propos une corrlation pour calculer le coefficient de perte de
charge, qui correspond sensiblement aux donnes de labaques de Moody pour lcoulement turbulent :
2.51

1
= 2 log 10
+ r
3
.71

Re

(51)

Cette expression tant implicite en , on procde par approximations successives, en partant dune valeur initiale
0. Citons galement la corrlation de Blasius, valable uniquement pour les tubes lisses, mais qui peut tre
utilise pour calculer une valeur initiale 0:
= 0,3164 Re -1/4, pour 2000< Re < 105 (Blasius)

(52)

Par exemple, pour Re = 106 , et r =0,01, on pose 0 = 0,3164 Re 1 / 4 =0,01


On calcule alors 1 avec (51) :
2.51

= 2 log10
+ r
1
Re o 3.71

1= 0,0379

puis 2 daprs 1 , et ainsi de suite jusqu ce que i+1 i.


Ici, on trouve :
2.51
1

= 2 log 10
+ r
2
Re 1 3.71

2= 0,0379

Cette dernire valeur correspond sensiblement celle donne par labaque de Moody :
Moody= 0,038.

16

figure 10 : abaque de Moody

17

Conduites non cylindriques.


Dans ce cas (conduites ovales, carrs, rectangulaires), les quations (44) (51) restent valides condition de
remplacer le diamtre D de la conduite circulaire, par le diamtre hydraulique Dh de la conduite considre,
dfinie comme tant gale 4 fois le rapport de la section du conduit par le primtre mouill P:
Dh =

4 S
P

Pour une conduite cylindrique, on retrouve bien videmment le diamtre de la conduite :


D 2
P = D

Dh = D
,
4
Pour une conduite rectangulaire (l x h), on aura par exemple :
S=

2lh
l+h
Ce concept de diamtre hydraulique peut tre interprt de la manire suivante : une conduite non circulaire de
& S , la mme
diamtre hydraulique Dh et de section de passage S provoquera, pour une vitesse moyenne u = V
perte de charge linaire quune conduite circulaire de diamtre D=Dh pour la mme vitesse u. Puisque les pertes
de charge sont exprimes en fonction de la pression dynamique (u2/2) et donc de la vitesse u, ce concept permet
de dterminer le comportement hydraulique de toutes conduites partir des donnes et des corrlations
dveloppes pour les conduites circulaires.
S = lh ,

P = 2 ( l + h)

Dh =

& ) , diamtre quivalent Deq.


Remarque : Abaques p = f ( V
Certains ouvrages et documentations techniques prsentent des abaques et corrlations permettant de dterminer
les pertes de charge en fonction non pas de la vitesse (ou de la pression dynamique), mais directement en
fonction du dbit volumique. Comme prcdemment, ces abaques ont t tablis pour des conduites circulaires.
Mais dans ce cas, le concept de diamtre hydraulique ne sapplique pas puisque pour une mme vitesse moyenne
& sera diffrent.
u, le dbit volumique V
Pour pouvoir utiliser les donnes relatives aux conduits cylindriques, le concept de diamtre quivalent Deq a t
& , la mme
dfini: il reprsente le diamtre dune conduite circulaire qui provoquerait, mme dbit volumique V
perte de charge p. La dtermination de cette grandeur est exprimentale. Ce diamtre quivalent est donn en
gnral par le constructeur de conduit, sous forme de corrlations ou dabaque. Le problme est ici que
lexpression de ce diamtre dpend de louvrage considr. Par exemple, si on reprend lexpression des pertes de
charge (44):
p
U2
= ( Re, r )
L
2D
En remplaant D par le diamtre hydraulique Dh = 2 a b ( a + b) , dun conduit rectangulaire a x b , on obtient :
p
( a + b) u 2
= ( Re, r )
L
4ab
& = u S = u ab
et, en faisant apparatre le dbit volumique : V
p
& 2 ( a + b)
= ( Re, r ) V
L
4a 3 b 3
Daprs la dfinition du diamtre quivalent, il faut chercher le diamtre de la conduite circulaire, qui, mme
dbit volumique, prsenterait la mme perte de charge. Daprs (44) et dune manire analogue ci-dessus, on
trouve :
p
&2 8
= ( Re, r ) V
,
L
2 D5

montrant que Deq doit vrifier : 8 2 D eq 5 = ( a + b) 4a 3 b 3


soit :

18

D eq

a 3b 3

= 1,265
a + b

0, 2

(52)

(32) correspond donc lapplication directe des dfinitions du diamtre hydraulique et du diamtre quivalent.
On retrouve un expression proche de celle-ci recommande par lASHRAE :
0.625
ab)
(
D eq = 1,30
( a + b) 0.250

alors que le note EUROVENT prconise, avec la condition a/b

D eq

1+ n
2 n
a
1 +

= 2b

3
a

4:

1
n 5

, avec n =

1
1,05 Log(Re) 0,45

qui prend en compte linfluence du nombre de Reynolds.


Ceci montre quavant dutiliser une abaque donnant la perte de charge en fonction du dbit volumique, il faut
rechercher dans le mmes ouvrage lexpression de Deq prconise.
En conclusion, on se souviendra de la signification de ces deux diamtres caractristiques pour les conduites non
circulaires :

Dh
: mme vitesse moyenne, mme perte de charge ,

Deq
: mme dbit volumique, mme perte de charge .

Pertes de charge singulires.


Lors de lcoulement dun fluide dans une conduite, les pertes de charge rgulires correspondent une
dissipation dnergie sous forme de chaleur lorsque le rgime est tabli, cest dire lorsque le profil de vitesse
dans la veine fluide reste inchang le long de la conduite. Par opposition, les pertes de charge singulires
apparaissent lors de la traverse de composants qui modifie le profil de vitesse (en grandeur et/ou en direction)
dans la veine fluide, et donc va engendrer une modification de la pression dynamique (les anglo-saxons appellent
ces pertes de charge singulires dynamic losses ). La figure 11 donnent quelques exemples de dformations
dun profil. Sans entrer dans les dtails, on voit que cette dformation a pour consquence daccrotre les
gradients de vitesse au sein du fluide, et donc de provoquer une dissipation dnergie supplmentaire.
Les pertes de charge singulires apparatront donc chaque changement de direction (coude), de vitesse
(largissement, rtrcissement), chaque sparation ou jonction, et enfin chaque entre/sortie.

Figure 11 : les pertes de charge singulires sont dues la dformation du profil de vitesse.
On exprime les pertes de charge singulires par une expression analogue (44) :
U2
2
avec : coefficient de perte de charge singulire (sans dimensions).
p =

19

(53)

Le facteur est donne soit par le constructeur de llment considr, soit par des abaques ou corrlations que
lon peut trouver dans des ouvrages spcialises (cf. bibliographie en fin de polycopi)
Pour un avant-projet sommaire, on pourra utiliser en premire approximation les valeurs donnes ci-dessous :

pour un coude 90, arrondis section circulaire, carr ou rectangulaire grand cot perpendiculaire au
rayon de courbure

0,5
1

0,75
0,45

1
0,30

1,5
0,20

2
0,20

pour un coude 90, arrondis section rectangulaire grand cot suivant le rayon de courbure
R/d

0,5
1,5

0,75
0,7

1
0,45

1,5
0,3

2
0,3

avec arrondi extrieur = 1

pour un coude 90, brusque,

sans arrondi extrieur = 1,5

largissement brusque, en ligne droite


S1

S2

1 rapport la vitesse dans S1 : 1 = 1 S1 S

sortie de conduit 1=1 (S2 = )

rtrcissement brusque, en ligne droite


2 rapport la vitesse dans S2
S
S
2 =0,5
entre de conduit 2=0,5
1

Quelques ouvrages utiles.

Guide de lAICVF , Araulique : Principes de laraulique appliqus au gnie climatique, PYC Edition,
Paris, 1991.
RECNAGEL-SPRENGER, Manuel pratique du gnie climatique, PYC Edition, Paris, 1995.
ASHRAE Handbook: 1989 Fundamentals.
IDELCHICK Handbook of hydraulic resistance

20

21

2me PARTIE

VENTILATEURS

1 DESCRIPTION D' UNE MACHINE CENTRIFUGE ......................................................................................................23


1.1 La roue
23
1.2 Le distributeur
24
1.3 Le diffuseur 24
1.4 La volute
24
1.5 Le divergent 24
2 VENTILATEUR CENTRIFUGE IDEAL....................................................................................................................24
2.1 triangle des vitesses 24
2.2 quation de continuit25
2.3 quation de Bernoulli 25
2.4 Caractristiques de fonctionnement
26
3 VENTILATEUR AXIAL IDEAL. ............................................................................................................................27
3.1 Triangle des vitesses 27
3.2 Caractristiques de fonctionnement
28
4 CARACTERISTIQUE DE LA MACHINE IDEALE, COEFFICIENT DE RATEAU ..............................................................28
5 MACHINES REELLES, RENDEMENT....................................................................................................................29
6 UTILISATION DE LA SIMILITUDE .......................................................................................................................30
7 SELECTION D' UN VENTILATEUR. .......................................................................................................................30
8 REGLAGE DES VENTILATEURS ..........................................................................................................................31
ANNEXE : ANALYSE DIMENSIONNELLE ...................................................................................................33
9 PRINCIPE, THEOREME DEVASHY-BUCKINGHAM (OU THEOREME DES )..............................................................33
10 EXEMPLES ...................................................................................................................................................33
10.1
Rsistance lavancement dun bateau.
33
10.2
Pertes de charges dans une conduite.
35

22

Les ventilateurs constituent le moteur du rseau et doivent donc tre adapts le plus prcisment possible au rseau, de
faon satisfaire le cahier des charges (dbit correct) tout en gardant un comportement nergtique optimum. Deux
types sont principalement utilis en araulique :

les ventilateurs axiaux (ou hlico des),

les ventilateurs centrifuge.


Ces deux types font partie des turbomachines. Leurs noms proviennent du mouvement de lair durant la compression
(figure 1). Pour les ventilateurs axiaux, lair a un mouvement hlico dale, cest dire que la distance entre la veine
fluide et laxe de la roue reste constante. Dans un ventilateur centrifuge en revanche, lair sloigne de cet axe au fur et
mesure de la traverse du ventilateur.
axial

centrifuge

r=cste

Figure 1 : ventilateurs centrifuges et axiaux.


Il rsulte de cet loignement de laxe une augmentation de la vitesse de lair, et donc de la pression dynamique qui va
sajouter la pression statique cre par les pales. Ce type de ventilateur pourra donc crer une diffrence de pression
totale entre lamont et laval suprieure celle que pourra crer un ventilateur axial. En revanche, il apparat clairement
qu diamtre de roue gale, la section dentre du ventilateur centrifuge est plus faible, et par suite sa capacit de dbit
volumique sera infrieure celle du ventilateur axial.
description d'une machine centrifuge
La figure 2 reprsente les diffrents lments constituant une machine centrifuge.

figure 2 : reprsentation schmatique d'une pompe centrifuge


La roue
Elle comporte des canalisations formes par le disque, les flasques de recouvrement et les aubages.
23

Le distributeur
Le rle du distributeur, situ en amont de la roue (non reprsent sur la figure) est de permettre une orientation correcte
des filets fluides leur entre dans la roue. Il n'existe en fait que dans les pompes plusieurs tages de compression. Il
permet alors d'effectuer la liaison entre l'tage amont (au niveau de la volute) et l'tage aval (au niveau des ou es).
Le diffuseur
Avec ou sans aubages, cet organe permet de transformer en partie l'nergie cintique (pression dynamique) en nergie
potentielle (pression statique) et d'orienter correctement le fluide la sortie de la roue pour son entre dans la volute ou
dans le distributeur de l'tage suivant.
La volute
Son rle consiste essentiellement transformer le mouvement de rotation du fluide en mouvement de translation.
Le divergent
Comme le diffuseur, il permet de transformer une fraction de l'nergie cintique du fluide en nergie potentielle, ceci
afin de limiter les pertes de charge dans les canalisations de transport du fluide.
Ventilateur centrifuge idal.
triangle des vitesses
Le vecteur vitesse du fluide en un point quelconque de la roue peut tre dcompos en trois composantes (figure 2a) : ca
:vitesse axiale, cd : vitesse diamtrale (ou radiale) et ct : vitesse tangentielle.
En introduisant la vitesse mridienne cm =cd +ca, (figure 2b), on dfini le plan mridien (figure 2c), dans lequel la
vitesse absolue c peut tre dcompose selon deux mouvements:

un mouvement d'entranement, qui correspond au mouvement du rotor, et donc du plan mridien: c e = r ,


r tant la distance l'axe et la vitesse angulaire du rotor,

un mouvement relatif, qui correspond au mouvement du fluide par rapport au rotor.


En un point quelconque, on a alors : c = ce + cr. Le triangle des vitesses est la reprsentation des compositions de ces
diffrents vecteurs vitesse dans le plan mridien. On le trace en gnral l'entre et la sortie du rotor (figure 4).

ca

c
cd

cm

ct

ct
(a) : C = Cd +Ca +Ct

(b) : C = Cm +Ct

cm

cr

ce
C = Cr +Ce

Plan mridien

ct

(c)

Figure 3 : dfinitions des vitesses et du plan mridien


A l'entre de la roue, afin de rduire les chocs, la vitesse relative du fluide doit tre tangente l'aubage ( = ' ). Dans
le cas contraire, il y a choc l'entre de la roue ce qui diminue le rendement. A la sortie, le fluide tant guid par
l'aubage, la condition ' = est toujours vrifie.

24

c2
cr2

cm2
2

cr1

cm1

ce2
ct2

c1

1
1 c
t1

'1

ce1

Figure 4 : Triangle des vitesses en entre (indice 1) et en sortie (indice 2) de la roue.


quation de continuit
L'quation de continuit (conservation de la masse) peut s'crire ici selon la conservation des dbits volumiques puisque
le fluide est considr comme incompressible. Ce dbit est gal au produit de la vitesse mridienne par la section de
passage normale la direction de l'coulement relatif. En premire approximation, cette section de passage est gale
2 r.e, o e reprsente la largeur de la section normale l'coulement mridien. On peut alors crire :
& = 2 r e c = 2 r e c = 2 r e c
V
m
1
1
m1
2
2
m2

[1]

Notons qu' la sortie, l'coulement est toujours perpendiculaire l'axe. On a donc cm2 = cd2.
Pour un dbit volumique donn et un rayon de sortie donn, [1] montre que cm2 ne dpend pas de l'angle de sortie 2. De
mme, la vitesse d'entranement ce2 ne dpend que du rayon r2 et de la vitesse de rotation . Le triangle des vitesses en
sortie (fig 4) montre alors que la vitesse absolue c est d'autant plus grande que l'angle 2 est faible.
Ainsi un rotor aubage en avant (2 </2) conduit une vitesse absolue de sortie suprieure celle qui est procure par
un rotor aubage en arrire (2 > /2). Les pertes de charge tant proportionnelles aux vitesses d'coulement, les
aubages en avant sont utiliss dans les ventilateur de faible puissance, les aubages en arrire tant rservs aux
puissances plus leves (figure 5).
2>
/2

2=
/2

cd2
cr2

cr2 = cd2

c2
2

2<
/2

c2

ce2

c2

cr2

cd2
2

ce2

ce2

Figure 5 : influence de 2 sur la vitesse de sortie c2.


quation de Bernoulli
On s'intresse ici, au ventilateur idal : il n'y a donc ni choc l'entre ('1 = 1) ni frottement (fluide incompressible
idal : hf = 0). L'application de l'quation de Bernoulli applique entre l'entre (1) et la sortie (2) s'crit alors:
p 2 c 22
p

c2
+
+ z 2 1 + 1 + z 1 = H m

2g
2g

[2]

En mouvement relatif, le thorme de Bernoulli crit le long de la trajectoire relative amne:


p 2 c 2r 2
p
c 2 c 2e1
c2
+
+ z 2 1 + r 1 + z1 = e 2

2g
2g
2g

[3]

Or, d'aprs le triangle des vitesses (figure 4), on a


c 2r = c 2m + (c e c t ) 2 = c 2 c t 2 + (c e c t ) 2 = c 2e + c 2 2c e c t
En combinant [2], [3] et [4] on obtient finalement:

25

[4]

H m =

c e 2 c t 2 c e1 c t1
g

[5]

Frquemment, l'entre du fluide dans le rotor se fait sans prgiration. La composante tangentielle de la vitesse d'entre
est alors nulle. La hauteur manomtrique du ventilateur devient alors:
H m =

c e2 c t 2
g

[6]

Dans le cas o le fluide entre dans le rotor avec un mouvement de prgiration (obtenu avec un distributeur, appele
ventelles pour un ventilateur), [5] devient:
H m =

ce2
g

r
c t 2 1 c t1
r2

[7]

Caractristiques de fonctionnement
Il est intressant de connatre Hm non pas en fonction de la vitesse, mais en fonction du dbit volumique
&
V
D'aprs l'quation de continuit [1], on a: c m =
.
2 r e
c
c
D'autre part, le triangle des vitesses amne : c t 2 = c e 2 + m 2 et c t 1 = m1
tg 2
tg1
En introduisant ces deux expressions dans [7], on obtient:
H m =

&
c e22
ce2
ce2
V
+

g 2 r2 e 2 g tg 2 2 r2 e1 g tg 2

[8]

En remarquant que la vitesse d'entrainement ce2 est lie la vitesse angulaire par: ce2 = r2., il vient:
r2
Cotg 2
Cotg1
& = k 2 + k V
&
V

H m = 2 2 +
1
2
g
2
e
g
2
e
g

2
1

[9]

o k1 et k2 ne dpendent que des dimensions gomtriques de la roue. k1 est toujours positif, alors que le signe de k2
dpend des angles 2, et 1. Dans le cas d'un fonctionnement sans prgiration ( 1 = /2, cotg1 = 0), on a :

2 < /2 cotg 2 > 0

k2 > 0

2 = /2 cotg 2 = 0

k2 = 0

2 > /2 cotg 2 < 0

k2 <0
& en fonction du dbit volumique
A chaque cas correspond un type de variation de hm et donc de la puissance utile W
u
&
V:
& = p V
& = g H V
&
W
u
m
m

[10]

La figure 6 reprsente schmatiquement ces volutions selon l'angle 2.


h

&
W
u

< /2

< /2

= /2

= /2

k1.2

> /2

> /2
&
V

&
V

Figure 6 : courbes caractristiques des ventilateurs centrifuges idaux.


En ce qui concerne l'angle d'incidence l'entre 1 (prgiration), sa valeur reste proche de /2. Par suite, l'influence de
ce paramtre est moins sensible, mais permet une certaine plage de rglage lorsque le ventilateur est muni de directrices
orientables (ventelles). Par exemple, pour une roue raction (2 >/2, k2 <0):
26

si 1 < /2
cotg 1 > 0. Le coefficient
k'2 devient suprieur en valeur absolue au coefficient k2 correspondant un fonctionnement sans prgiration.
& ) devient suprieure au cas de non-prgiration.
Ainsi, la pente de la droite h( V
si 1 < /2
cotg 1 > 0. On assiste au
phnomne inverse.

/2
>
/2
<

2 < /2

1 > /2
1 < /2

2 = /2

>
/
< 2
/2
1

2 > /2

&
V
Figure 7 : influence de la prgiration sur les caractristiques d'un ventilateur centrifuge idal
1

Ventilateur axial idal.

La figure 8 reprsente schmatiquement un ventilateur axial. La trajectoire tant hlico de, on remarquera que les
vitesses d'entranement l'entre et la sortie sont identiques.
ce2 (= ce1)
ct

ce1

ca = cm

figure 8 : ventilateur axial


Triangle des vitesses
1

cm =ca

c1

1
2
cr1

c2

cr2

ct2 -ct1

Figure 9 : triangles des vitesses pour un ventilateur axial.


On adopte ici une reprsentation condense des deux triangles des vitesses (entre et sortie) puisque la vitesse
d'entranement unique peut servir de base commune aux deux triangles. D'autre part, si les sections de passage l'entre
27

1 et la sortie 2 sont gales, la composante axiale de la vitesse ca, gale la composante mridienne cm, a une valeur
unique, donne par l'quation de continuit:
c a1 =

&
&
V
V
=
= c a2 = ca = c m
1 2

[11]

On obtient alors les deux triangles des vitesses reprsents figure 9.


Caractristiques de fonctionnement
On suit le mme raisonnement que pour les ventilateurs axiaux:

p
p
c2
c2
Bernoulli : 2 + 2 + z 2 1 + 1 + z 1 = H m

2g
2g

[12]

c 2 c 2e1
p
p
c2
c2
=0
Bernoulli en mouvement relatif : 2 + r 2 + z 2 1 + r1 + z 1 = e 2
2g

2g
2g

[13]

On peut vrifier que [4] est toujours valable et en combinant cette dernire quation avec [12] et [13], on obtient:
H m =

ce
( c t 2 c t1 )
g

[14]

Le triangle des vitesses (figure 9) permet de relier ct2 - ct1 aux deux angles 1 et 2.
c t 2 = c e + c a cot g 2

c t1 = c a cot g1

c t 2 c t1 = c e + c a (cot g 2 cot g 1 )

ce qui amne, en utilisant [11]:


H m =

&
c 2e c e V
+
(cot g 2 cot g 1 )
g
g

[15]

En faisant apparatre la vitesse angulaire (ce = .r), on obtient une relation analogue [9]:
r2
r

& = k' 2 + k' V


&
H m = 2 +
(cot g 2 cot g 1 ) V
1
2
g
g

Les caractristiques obtenues sont donc analogues celles des ventilateurs centrifuges. La figure 10 donnent
schmatiquement ces caractristiques tablies pour un angle 2 >0 et plusieurs valeurs de 1.

k'1.2
1 > /2

1 < /2

1 = /2

&
V

Figure 10 : caractristiques d'un ventilateur axial idal (2 > /2)


2

Caractristique de la machine idale, coefficient de Rateau

On s'intresse ici une machine idale (pas de frottement ni de choc) traitant un fluide incompressible. D'aprs ce qui
prcde, on a:
2
&
H = k1 + k2 . . V

[16]

28

On en dduit la puissance mcanique transmise au fluide :


& =V
& P = V
& g H = g (k 2 V
& + k V
& 2)
W
1
2

[17]

On peut adimensionnaliser l'quation [16] en divisant chaque membre par le rapport ( 2 r22 ) g , homogne une
longueur. On obtient alors :
&
V
gH
2 2 = 1 + a (1 , 1 , r2 )
3
r2
r2

[18]

En remarquant que la vitesse de l'aubage en sortie de roue u2 = .r2, et en introduisant les groupements adimensionnels:
g H g H
= 2 2= 2
r2 u 2
& V
&
V
=
=
3
2
r2 u 2 r2
[20] prend la forme de l'quation d'une droite de coefficient directeur a, et d'ordonne l'origine 1:
=1+a.

[19]

On peut galement adimensionnaliser l'quation [18] donnant la puissance transmise au fluide en faisant apparatre les
trois groupements , a et :
2 2

& = g r2 V
& + a ( , , r ) V
& 2 = u 2 V
& + a ( , , r ) u 2 V
& 2
W

2
1
2
2
1 2 2
g

r
2

soit,
2

&
&
& = (u 3 r 2 ) V + cot g 2 r2 ( u 3 r 2 ) V
W
2
2
2
2
2
2
e2
u 2 r2 2
u 2 r2

Finalement, on obtient la forme adimensionnelle suivante :

&

W
= + a 2
=
3
2
(
u
r
)

2
2

[20]

On peut remarquer que le coefficient [a] caractrise les dimensions gomtriques de la roue. Les trois autres coefficients
sont nomms:

: coefficient manomtrique de Rateau,

: coefficient de dbit de Rateau,

: coefficient de puissance de Rateau.


[19] et [20] montrent que toutes machines semblables (ou homologue, c'est dire prsentant le mme coefficient [a])
auront les mmes caractristiques adimensionnelles = f() et = f().
3

Machines relles, rendement.

Les quations [16 19] prcdentes ne sont valables que pour une machine idale, c'est dire sans frottement, ni choc,
& ) et W
& = f (V
& ) s'en
ni fuite entre aspiration et refoulement. Pour les machines relles, les caractristiques H = f (V
trouvent modifies. La figure 11 prsente titre d'exemple l'volution de la hauteur manomtrique en fonction du dbit
volumique pour une machine centrifuge aubages en arrire. la hauteur thorique H dont la variation est
linaire(quation [16] et figure 6), il faut soustraire:

Les pertes de charge rgulires (dues uniquement aux frottements) qui varient sensiblement avec le carr du dbit,

Les pertes de charge singulires (essentiellement les chocs) qui passent par un minimum pour un dbit particulier,
& .
correspondant une entre dans la roue sans choc V
SC
On obtient alors une caractristique d'allure parabolique. La mme figure reprsente galement l'volution de la
& lorsque l'on prend en compte de la mme faon ces diffrentes pertes d'nergie.
puissance reue par le fluide W
f

29

&
W
f

pertes dues aux chocs


pertes dues aux frottements

&
V

&
V

Figure 11 : influences des pertes nergtiques sur les caractristiques relles d'une roue aubages courbs vers
l'arrire.
Ces pertes nergtiques conduisent dfinir le rendement global de la machine:
g =
&
avec W
a

&
&
W
gH V
f
=
&
&
W
W
a
a

[21]

: puissance absorbe sur l'arbre.

& , W
& , et g.
La figure 7 prsente les courbes de fonctionnement d'un ventilateur centrifuge : H, W
f
a

h
&
W
a

g
&
W
u
&
V
Figure 12 : Caractristiques relles d'un ventilateur centrifuge
4

Utilisation de la similitude

Bien que les caractristiques des ventilateurs relles sont diffrentes de celles issues de l'tude de la machine idale, les
rgles de similitude restent valables. Les quations [19] et [20] sont simplement remplaces par des lois plus complexes
= f(), f = f(),
Si l'on connat les caractristiques relles d'une machine dans certaines conditions de fonctionnement (par exemple une
vitesse de rotation) :
&), W
& = f (V
&), W
& = f (V
& ) et = f (V
&),
H = f (V
g
f
a
on peut tout d'abord en dduire les caractristiques adimensionnelles:
= f(), f = f(), a = f() et g = f().
Un fois ces caractristiques connues, on peut calculer les caractristiques relles pour une autre vitesse de rotation, ou
bien en dduire les caractristiques relles d'une autre machine semblable (c'est dire qui prsente le mme coefficient
[a]).
5

Slection d'un ventilateur.

Les ventilateurs sont caractriss par leurs courbes dbit-pression (ou dbit-hauteur manomtrique), qui sont diffrentes
suivant le type de ventilateur (figure 13) :

ventilateur centrifuge : la caractristique part dun dbit nul, augmente jusqu une valeur maximum, puis diminue
lorsque le dbit continue augmenter.
30

ventilateur axiaux : on a thoriquement le mme type de caractristique, mais faible dbit ( gauche du maximum),
il existe une zone dinstabilit dans laquelle le ventilateur dcroche : langle dincidence de lcoulement par
rapport aux pales et trop lev et la roue emprisonne un certain volume dair, auquel elle ne communique quun
mouvement de giration. Aux dbits encore plus faibles, il peut enfin stablir un mouvement centrifuge : entre prs
de laxe o la vitesse est la moins leve, et refoulement loin de laxe.
Dans les deux cas, la zone dutilisation est situe droite du maximum, afin dviter le phnomne de pompage.

Pt ou ht

zone
dutilisation

axial

centrifuge
&
V

Figure 13 : Caractristiques des ventilateurs axiaux et centrifuges, domaine dutilisation.


De tels caractristiques, fournies par le constructeur, sont issues en gnral de tests effectus sur le ventilateur associ
une longueur droite de tube permettant dtablir le rgime (cf. 1re partie), et appele longueur effective Le. Si dans le
montage rel cette longueur ne peut pas tre mise en place (piquages ou coudes placs une distance infrieure Le), il
faut tenir compte des effets dun tel branchement qui vont se traduire par une perte de charge singulire supplmentaire
exprime par la relation classique:
p =

U 20
2

: coefficient de perte du la connexion ventilateur/rseau,


Uo
: vitesse de sortie, calcule daprs la section de sortie pour les ventilateurs centrifuges (outlet area
sur la figure 20), daprs la section de la roue pour les ventilateurs axiaux (blast area sur la figure 20).

avec

Comme pour les autres pertes de charge singulires, est donn pour plusieurs assemblages spcifiques (voir
assemblages 7, annexe A). Sa valeur dpend de Le. Si Le est inconnue, on peut lestimer par :
si Uo 13 m/s, L e =

U o A o 0, 5
350

si Uo > 13 m/s, L e =

U o A o 0, 5
4500

relations empiriques issues de lexprience.


Pour les conditions dentre, le mme problme apparat mais les phnomnes sont plus complexes et il nest pas
possible dexprimer cet effet par une relation linaire de la pression dynamique. On se rappellera que des mauvaises
conditions dentre peuvent avoir des consquences dramatiques sur les performances dun ventilateur. Les conditions
idales rechercher doivent permettent une entre dair uniforme et axial, sans rotation ni turbulence.
6

Rglage des ventilateurs

Bien que les caractristiques relles des ventilateurs diffrent grandement de celles issues de la machine idale,
l'influence des angles 1 (prgiration) et 2 (angle de sortie des aubages) , ainsi que de la vitesse de rotation (equ 16]
reste valable qualitativement. Le rglage des ventilateurs peut donc s'effectuer:

par modification de la vitesse de rotation (figure 14a)

par modification de l'angle 1 par utilisation de ventelles l'aspiration (figure 14b).

par laminage (figure 15).


Ce dernier moyen consiste introduire une perte de charge supplmentaire sur le circuit de branchement du ventilateur
(figure 15). On remarque que dans le cas d'un ventilateur angle 2 lev, le laminage conduit une variation non
31

ngligeable de la hauteur manomtrique, ce qui grve considrablement le rendement global de l'installation. Ce mode
de rglage est donc mieux adapt au ventilateur faible b2, et en particulier aux ventilateurs aubages lgrement en
avant.

>/2
N3

N2

</2

N1

=/2

&
V

&
V

(a)

(b)

Figure 15 : rglage par variation devitesse (N1>N2>N3), rglage par modification de l'incidence l'entre.

s faible
h*

h*

s lev

&
V

&
V

Figure 16 : rglage par laminage.

32

Annexe : analyse dimensionnelle


7

Principe, thorme deVashy-Buckingham (ou thorme des )

Lanalyse dimensionnelle est fonde sur le principe fondamental trs simple selon lequel une relation entre grandeurs
physiques doit tre dimensionnellement homogne, cest dire indpendante du systme dunits de mesure choisi.
Ce principe de base est traduit par le thorme deVashy-Buckingham (ou thorme des ) que lon peut noncer ainsi :
Toute fonction G de p variables indpendantes x1, x2, ,xp, mesure par q units fondamentales, et traduisant la
variation dun phnomne en fonction de p causes indpendantes, o p > q sexprime ncessairement partir dune
relation fonctionnelle de la forme :
g

G = x1g1 x g2 2 L x q q F(q +1 , q +2 ,L, p )


les variables x1, x2, ,xq tant choisies dimensionnellement indpendantes. On dmontre en outre que les fonctions
i, sont des groupements adimensionnels des variables x1, x2, ,xp.
Dans la pratique, on choisit pour x1, x2, ,xp les paramtres que lon considre comme essentiels pour le problme
considr et que lon veut voir figurer explicitement dans lexpression de G.
Remarque :
En mcanique, on a habituellement 3 units fondamentales :

la longueur L,

la masse M,

le temps T,
do q = 3 est lexpression prcdente devient :
G ( x1 , x 2 ,L, x p ) = x1g1 x g22 x 3g3 F( 4 , 5 ,L, p ) .
Par contre, ds que des phnomnes thermiques interviennent, il convient dintroduire la temprature comme
quatrime grandeur fondamentale, do q = 4 est lexpression de G devient :
G (x1 , x 2 ,L, x p ) = x1g1 x g2 2 x g3 3 x g44 F(5 , 6 ,L, p )
8

Exemples

Rsistance lavancement dun bateau.


Supposons que G reprsente la rsistance lavancement dun bateau. On peut penser que ce phnomne dpend de la
taille l et de la vitesse u du bateau, de la masse volumique , de la viscosit et de la tension superficielle de leau
ainsi que de lacclration de la pesanteur g.
Choix des variables principales
La rsistance lavancement est une force, dont la dimension est [LMT-2]. On a donc q=3 (M,L,T). Il est assez
logique de considrer quune grandeur caractristique de sa taille (longueur l), sa vitesse et la masse volumique de leau
sont choisir comme variables principales, de prfrence la tension superficielle, la viscosit ou lacclration de
la pesanteur. Ces deux derniers paramtres sont importants mais, si la viscosit pourrait la rigueur tre choisie, g ne
constituerait pas avec L et u un systme de variables dimensionnellement indpendantes, puisque lon a :
[g]=LT-2 =

(LT 1 ) 2 [u]2
=
L
[l]

Dtermination de la loi G(,u,l,,g,)


On a q=3 et p=6 . On cherche donc une relation du type :
G (, u, l, , g, ) = g1 u g 2 l g3 F(4 , 5 , 6 )
On dresse alors le tableau suivant o sont regroupes les dimensions des diffrentes variables, ainsi que les
exposants ei quil faudra leur affecter pour vrifier lquation aux dimensions de lexpression:
G = e1 u e2 l e3 e4 g e5 e6
La premire ligne correspond la grandeur recherche (R en N), les trois suivantes aux variables principales , u et
l. On termine par les trois variables restantes , g, .
33

-2

-3

e1

-1

e2

e3

-1

-1

e4

-2

e5

-2

e6

Daprs ce tableau, on a
1 = -3e1 + e2 + e3 - e4 + e5
1 = e1 + e4 + e6
-2 = - e2 - e4 - 2e5 - 2e6
On exprime alors les exposants des variables principales en fonctions des autres. Il vient:
e1 = 1- e4 - e6
e2 = 2 - e4 - 2e5 - 2e6
e3 = 1 + 3(1- e4 - e6 ) - (2 - e4 - 2e5 - 2e6) + e4 - e5 = 2 - e4 + e5 - e6
R sexprime alors par :
R = 1e4 e6 u 2e4 2 e5 2 e6 l 2e4 +e5 e6 e4 g e5 e6
On regroupe alors les termes en e4, e5 et e6 pour faire apparatre les groupements adimensionnels :
e

e
4 gl 5

2
R = u 2 l 2
2

u l u u l

e6

On voit alors apparatre trois nombres adimensionnels :


ul
,qui traduit linfluence de la viscosit, donc du frottement de leau sur la
le nombre de Reynolds : Re =

coque,
u
, qui traduit linfluence de la pesanteur, donc du sillage (systme de vagues

le nombre de Froude : Fr =
gl

produit par le bateau),


u2 l

le nombre de Weber : Wb =
, qui traduit linfluence des forces de tension superficielles.

On peut donc crire :

R =

u2
S F(Re, Fr, Wb)
2

en posant S = l2 et en introduisant le coefficient pour faire apparatre lnergie cintique.


Ainsi nous voyons que, sans nous donner la forme prcise de la relation, lanalyse dimensionnelle nous montre
comment interviennent les paramtres. Ltude exprimentale de la rsistance lavancement dun bateau se ramne
alors ltude de la fonction F, que lon appelle le coefficient de rsistance lavancement, en fonction de trois
paramtres seulement.

34

Pertes de charges dans une conduite.


Les variables prendre en compte sont ici :

le diamtre D et la longueur L de la conduite,

la vitesse u, la masse volumique et la viscosit du fluide,

lacclration de la pesanteur g.
On a donc p = 5
Choix des variables principales
Une perte de charge peut tre considre soit comme une perte dnergie (dimension [ML2T-2]), soit comme une chute
de pression (dimension [ML-1 T-2]). Notons que son expression en hauteur (dimension [L]) ne peut pas tre prise ici,
puisquune hauteur est toujours relatif un fluide donn (cf. lexpression m CE ou m.CA), cest dire un poids
volumique donn (dimension [ML-2 T-2]). On a donc q=3 (M,L,T).
Comme prcdemment, on peut choisir , u ,et une dimension caractristique de la conduite (D).
Lexpression recherche, par exemple en pression, est donc de la forme :
P(, u, D, L, , g) = g1 u g 2 D g 3 F( 4 , 5 , 6 )
Dtermination de la loi P(,u,D,l,,g)
L

e1

e2

e3

e4

e5

e6

Dtermination de la loi P L (,u,D,,)

P/L

e1

e2

e3

e4

e5

35

36

3me PARTIE

VENTILATION DES LOCAUX

1 OBJECTIFS DU CONDITIONNEMENT D'AIR ..........................................................................................................38


2 DIFFUSION DE L'AIR .........................................................................................................................................40
2.1 Les jets d'air ...........................................................................................................................................40
2.2 Le soufflage de l'air.................................................................................................................................42
2.3 Les diffrents types de bouche .................................................................................................................43
3 VENTILATION PAR DEPLACEMENT ....................................................................................................................44
4 DEBIT DE SOUFFLAGE ......................................................................................................................................45
4.1 Dbit d'air neuf .......................................................................................................................................45
4.2 Dbit de soufflage ...................................................................................................................................45

37

Objectifs du conditionnement d'air

Un systme de conditionnement d'air maintient les conditions de temprature et d'humidit d'un local. Dans la plupart
des cas, de nombreuses solutions technologiques peuvent rpondre l'objectif souhait. Le rle de l'ingnieur est donc
d'analyser ces diffrentes possibilits, et de comparer leur comportement par rapport un objectif d'ambiance. En effet,
il faut se rappeler que le dimensionnement d'un systme se fait en gnral sur la base de 2 points de fonctionnement,
correspondant aux situations les plus svres en t et en hiver. L'objectif d'ambiance consiste dfinir une zone dans
laquelle les conditions de temprature et d'humidit sont considres comme acceptables par la majorit des occupants.
Cet objectif dpend donc de la finalit du local (locaux d'habitation, bureaux, ateliers ...) et est issu d'tudes la fois
physiologiques et psychologique, puisque la notion de confort reste trs subjective. Elle dpend en particulier des
conditions climatiques extrieures: dans le secteur rsidentiel par exemple, il est bien connu qu'une ambiance 18/19C
sera ressentie comme confortable en hiver, mais peut apparatre frache et mme froide si la temprature extrieure est
suprieure 35C. La figure 1 prsente la zone de confort dfinie par l'ASHRAE, pour les deux modes hiver et t.

Figure 1 : zone de confort, d'aprs ASHRAE


Outre les conditions de temprature et d'humidit, la qualit de l'air comprend galement un critre de puret. Trois
mthodes d'essai de filtres coexistent afin de couvrir toutes les conditions d'utilisation. Il en rsulte trois dsignations de
filtres selon leur efficacit:

efficacit gravimtrique,

efficacit opacimtrique,

trs haute efficacit.


Les deux premires concernent les installations de traitement d'air "communes" (rsidentiel, tertiaire), la troisime
s'adapte aux locaux spcialiss (salles blanches). Le tableau 1 prsente les zones de recouvrement de ces diffrentes
catgories, et leurs classes EU correspondantes dfinies par EUROVENT. La lgislation impose pour les locaux
affects au travail une filtration de l'air recycl sur des filtres d'efficacit opacimtrique de 50% minimum. En fait, les
organismes professionnels recommandent plutt une efficacit opacimtrique de 85%, pour lutter contre l'accumulation
d'arocontaminants, plus sre pour la captation des bactries, et d'un cot d'exploitation peine plus lev.
Les filtres de classe 1 4 (efficacit gravimtrique) places en amont du ventilateur de soufflage (appel quelquefois
prfiltres) permet de protger les lments du traitement d'air (ventilateur, batteries, conduites de distribution et
bouches, appareils de contrle) en vitant leur encrassement.
Enfin, le traitement des odeurs apparat maintenant comme un autre critre important. Il peut tre soit trait par des
filtres spcifiques (par exemple, charbon actif), soit "bipass". Dans ce dernier cas, on ne peut pas parler de vritable
38

traitement d'odeurs puisque cela consiste ne pas recycler l'air provenant des locaux pollutions spcifiques (cuisines,
salle fumeur,...) et donc augmenter globalement le taux de renouvellement d'air.

Tableau 1: Classification des filtres


La conception d'une installation comprend donc plusieurs postes, tous interactifs. Il y a donc lieu de dfinir, ds
l'initiation du projet, les diffrents objectifs de confort que l'on s'impose, sans oublier les contraintes requises par la
lgislation. Il est difficile de gnraliser une mthode de dimensionnement puisque chaque btiment est un cas
particulier, et que plusieurs solutions peuvent rpondre au mme problme. Toutefois, on pourra suivre la procdure
suivante:
1. Calcul des charges hydriques et thermiques. En gnral, c'est le bureau d'architecture qui effectue le calcul des
charges sensibles externes. Le calcul des charges sensibles et hydriques internes doit tre effectu par le climaticien,
et ncessite de dfinir un scnario d'occupation.
2. Choix du systme : tout air, centralis ou localis, air/eau, tout eau.
3. Calcul du dbit d'air neuf, impos en gnral par la lgislation, et du dbit de soufflage, qui dpend du systme
choisi,
4. Dessin du rseau araulique et, le cas chant hydraulique et dimensionnement du (ou des) rseau(x),
5. Dfinition des quipements de chauffage et de rfrigration.
Nous nous intressons ici plus particulirement aux rseaux arauliques, et donc aux systmes tout air, centraliss ou
localiss.
Dans les systmes centraliss, l'air souffl doit rpondre trois objectifs principaux:

Renouvellement d'air suffisant,

Maintien des conditions de temprature,

Maintien des conditions d'humidit.


ces trois objectifs, se greffe une quatrime contrainte concernant la diffusion de l'air dans les locaux. C'est en grande
partie le sujet de ce polycopi.

39

Dans les systmes localiss, l'air souffl doit rpondre principalement aux contraintes du renouvellement d'air, et au
contrle de l'humidit. Le contrle de la temprature est alors effectue indpendamment dans chaque zone, soit par
ventiloconvecteurs, soit par climatiseurs locaux ("windows" ou "split").
Un systme de conditionnement d'air peut par ailleurs associer les deux systmes: l'installation centralise effectue un
prtraitement de l'air souffl, que le conditionnement local va ajuster en fonction de la demande.
10

Diffusion de l'air

Toute installation de ventilation comportant un rseau araulique de distribution d'air implique la mise en uvre d'un
systme efficace de diffusion de cet air, susceptible d'assurer une rpartition uniforme d'air frais sans provoquer une
sensation de courant d'air pour les occupants. Cette rpartition est influence principalement par l'emplacement, le type
et la vitesse d'air des dispositifs de soufflage ou d'aspiration.
10.1

Les jets d'air

10.1.1
Jet d'air isotherme
La figure 3, qui reprsente schmatiquement la forme d'un jet d'air isotherme (air souffl la mme temprature que
l'air ambiant), permet de distinguer quatre zones en fonction de la vitesse axiale de l'air.

Zone 1 : zone initiale. c'est une zone assez courte dont la longueur se situe entre 4 et 6 diamtres (ou largeurs de
fentes). La vitesse reste pratiquement constante dans le dard, qui la forme d'un cne de rvolution.

Zone 2 : zone de transition. Cette zone s'tend sur environ 8 diamtres et la vitesse axiale diminue
proportionnellement au carr de la distance l'ouverture.

Zone 3 : zone d'panouissement. La longueur de cette zone peut atteindre entre 25 et 100 diamtres, et dpend de
la forme et de la section de la bouche de soufflage, de la vitesse initiale et de l'espace dans lequel le jet d'air
s'panouit. L'coulement dans cette zone est turbulent, avec formation de petits tourbillons. La vitesse dcrot de
faon inversement proportionnelle la distance de la bouche.

Zone 4 : zone terminale. La vitesse du jet y dcrot rapidement (sur une longueur de quelques diamtres) jusqu'
une valeur de 0,25 m/s.

Figure 2 : Les diffrentes zones d'un jet d'air isotherme


Caractrisation des jets d'air
a)
Induction, vitesse rsiduelle
C'est dans les zones 2 4 qu'apparat le phnomne d'induction: le jet d'air principal met en mouvement l'air du local
(figure 3), provoquant son mlange avec l'air souffl, et conduit au mouvement de cet air de mlange une vitesse
rsiduelle vr dans la zone d'occupation. Cette vitesse rsiduelle est proportionnelle en premire approximation la
vitesse terminale vt: vr = R vt, R tant compris entre 0,35 et 0,5 pour les bouches usuelles.
Le tableau 2 prsente les vitesses d'air rsiduelles prconises pour diverses applications.
& et d'air
On dfinit le taux d'induction i d'une bouche par le rapport entre les dbits volumiques d'air de mlange V
m
& , que l'on peut galement exprimer en fonction du dbit d'air induit V
& :
souffl V
s

&
& +V
&
V
V
i
i= m = s
&
&
Vs
Vs
Tableau 2: vitesses rsiduelles recommandes
Locaux d'hbergement
40

vr = 0,12 m/s

Hpitaux
Locaux d'enseignement
Locaux de runions et de spectacle
Bureaux et locaux assimils
Locaux commerciaux
Atelier
Locaux sportifs
Grands magasins
Btiments du transport
Locaux industriels

vr = 0,17 m/s

vr = 0,25 m/s

Figure 3 : Induction d'air secondaire, porte et chute d'un jet d'air


b)
Porte
C'est la distance horizontale parcourue par un jet d'air sortant d'une bouche de soufflage (figure 3). Cette distance est
mesure depuis la bouche de soufflage jusqu' l'extrmit de la zone 4, c'est dire lorsque le jet atteint la vitesse
terminale vt. Les constructeurs dfinissent la porte en choisissant une vitesse de rfrence, souvent gale 0,2 m/s.
c)
Chute
C'est la distance verticale maximale entre l'axe du jet d'air et le point le plus bas relative une vitesse donne (vitesse
terminale vt) (figure 3).
La temprature de soufflage a une grande influence sur la chute, puisqu'un air plus froid que l'ambiance tendance
tomber alors qu'un air plus chaud a tendance monter. Les dispositifs de soufflage comprennent en gnral des
dflecteurs qui permettent d'orienter le jet suivant la temprature de soufflage. La figure 4 montre par exemple les
trajectoires thoriques suivies par un jet d'air selon la temprature de l'air souffl et l'angle de soufflage du jet.
10.1.2

Influence des parois

a)
Adhrence : l'effet Couanda
Nous avons vu qu'un jet d'air a tendance mettre en mouvement l'air dans lequel il s'panouit. Lorsqu'un air est souffl
prs d'une paroi, il va donc aspirer une partie de l'air situ entre le jet et la paroi (figure 5). Par suite, il apparat une
chute de pression statique prs de la paroi, qui va tendre plaquer le jet le long de cette paroi. Dans le cas d'un soufflage
horizontal, le jet d'air a donc tendance adhrer contre le plafond.
L'effet Couanda est trs utile en climatisation puisque dans ce cas, l'air souffl est une temprature infrieure l'air
ambiant et le jet a donc tendance s'incurver vers le bas (cf. ci-dessus). L'effet Couanda permet au jet d'air de pntrer
dans le local sans s'incurver trop rapidement.
On recommande en gnral une distance infrieure 0,3 m du plafond pour profiter de cet effet. Plus l'air souffl est
froid, plus la bouche de soufflage doit tre proche du plafond.
b)
Point de sparation
Malgr l'adhrence du jet au plafond d l'effet Couanda, le jet finit par ne plus adhrer au plafond et s'incurve pour
tomber dans la zone d'occupation (figure 5).

c)
Obstacle
La prsence d'un obstacle perpendiculaire au jet d'air (poutre, luminaire,...) peut perturber la diffusion de l'air et
supprimer l'adhrence (figure 6). On dfinit une hauteur critique de l'obstacle au-del de laquelle le jet d'air va chuter.
Cette hauteur critique peut tre estime par la corrlation suivante:

41

h c = d ( 0,08 0,005 T )
avec

d
T

: distance entre l'obstacle et la bouche de sortie,


: diffrence de temprature entre air souffl et air ambiant.

60

45

chauffage
refroidissement

30

x
Figure 4 : trajectoire d'un jet d'air en fonction de la temprature de soufflage et de l'angle initiale

dpression

Point de sparation

Figure 5 : Effet Couanda

Figure 6: Influence des obstacles


10.2
Le soufflage de l'air
Le soufflage de l'air doit permettre une induction suffisante pour avoir une temprature de l'ambiance quasi-constante
dans la zone d'occupation, dfinie par une hauteur de 1,8 m, et dont la surface au sol est dlimite par une distance de
0,15 m par rapport aux murs verticaux (cf. figure 7a). Trois implantations de bouches de soufflage sont possibles:

Soufflage mural en partie haute (figure 7a),

Soufflage en partie centrale de plafond (figure 7b),

Soufflage au sol ou en allge de fentre (figure 7c).


Le choix de la bouche doit conduire l'galit de sa porte P avec la dimension L du local dans le sens de l'coulement
du jet. Lorsque plusieurs bouches sont utilises, il y a lieu de vrifier que tout le local est bien trait (absences de zones
mortes) et que les bouches ne se perturbent pas mutuellement. Pour deux bouches en vis vis, on aura donc P = L/2.
42

1,80

vr

2,40 3,00

vt

(a)

0,15
(c)

(b)

vt

vt

vr

vr

Figure 7 : Les diffrents modes de soufflage


10.3

Les diffrents types de bouche

10.3.1
Les grilles
Une grille est une bouche d'air comprenant de multiples passages pour l'air (figure 8). Les plus simples comportent des
ailettes fixes, mais les types les plus utiliss comportent des ailettes rglables, verticales ou horizontales permettant une
certaine souplesse pour le rglage de la dflexion du jet. D'autres modles comportent la fois des ailettes verticales et
horizontales (double dflexion).

(a)

(b)

Figure 8 : grilles ailettes fixes verticales (a) et horizontales (b)


Les grilles doivent tre disposes de faon ce que le jet n'entre pas dans la zone d'occupation. Elles sont donc places
gnralement en partie haute de mur en soufflage horizontal, afin de profiter au mieux de l'effet Couanda, ou disposes
au sol ou en allge de fentre pour un soufflage vertical.
Certaines de ces grilles peuvent tre utilises indiffremment en soufflage ou en reprise. D'autres n'ont qu'une utilisation
: soufflage ou reprise.
Les dbits traits par les grilles sont trs varis suivant les modles, et la gamme couverte s'tend d'une centaine
quelques milliers de m3/h.
Souvent, des systmes de rglage de dbit sont associs aux grilles de transfert (figure 9).

43

(a)

(b)

Figure 9 : Rglage du dbit, par "pelle " (a), par clapets (b)
10.3.2
Les diffuseurs plafonniers
La structure des diffuseurs plafonniers permet l'aspiration d'air ambiant dans le diffuseur, qui facilitent le mlange avec
l'air primaire souffl, et conduit des taux d'induction importants (jusqu' 30). On peut ainsi accrotre la diffrence
entre la temprature de l'air souffl et de l'air ambiant, puisqu'il y a mlange avant le soufflage proprement dit.
10.3.3
Les bouches (ou diffuseurs) linaires
Elles ne se distinguent des grilles que par leurs dimensions gomtriques : leurs rapports longueur sur largeur tant
suprieurs 10. dbit et surface gaux, le traitement d'une zone est donc plus homogne et la forme du jet (trs plat)
conduit des taux d'induction plus importants.
10.3.4
Plafonds perfors
Il s'agit de composants raliss en tles perfores ou en matriaux poreux. Par leur grande section de passage, ils
permettent des dbits importants sous de faibles vitesses.
10.3.5
Classification des diffrentes bouches
Le tableau suivant donne quelques caractristiques moyennes en refroidissement suivant le dispositif de soufflage et
permet une prslection de matriel.
Tableau 3: Caractristiques des dispositifs de soufflage.
dispositif

grilles

bouches linaires

Diffuseurs
plafonniers
16 90

11 21
14 36
& (m3/h par m2 de
V
plancher)
Charge sensible
57
98
330
maximum (W/m2)
Taux de brassage*
7
12
30
(Volume/h)
8
8
11
T** (K)
vr (m/s)
0,15 0,175
0,1 0,175
0,1 0,25
* : hauteur de local de 3 m, ** entre temprature de soufflage et temprature ambiante.
11

Plafonds perfors
18 180
670
60
11
0,150

Ventilation par dplacement

Rcemment, il est apparu un nouveau concept de ventilation : la ventilation par dplacement, qui consiste introduire
de l'air basse vitesse dans la partie infrieure du local une temprature plus basse que celle de l'air ambiant.
On utilise ici le fait bien connu que l'air chaud s'lve naturellement dans une atmosphre plus froide (phnomne de
convection naturelle). L'air froid ainsi introduit se rpand sur toute la surface du plancher. Au contact des sources
thermiques (homme, machine), cet air se rchauffe et un courant ascendant se cre spontanment (figure 10). La
stratification de l'air qui se met en place conduit une concentration en polluant dans la zone d'occupation plus faible
que dans le cas d'une ventilation par diffusion. En effet, le principe mme de la ventilation par diffusion, qui implique le
mlange d'air souffl et d'air ambiant conduit une homognisation de la qualit de l'air dans tout le local, alors qu'ici,
les polluants sont entrans vers le haut du local, o se trouve la bouche de reprise.
La mise au point de tels systmes reste complexe, puisque, comme tout phnomne de convection naturelle, l'quilibre
qui se met en place est prcaire, et peut tre facilement perturb. Dans l'industrie, o la position des machines est
dfinie une fois pour toutes, la modlisation du local permet une mise au point optimum de ce concept. En revanche,
son application au tertiaire ou au rsidentiel reste difficile, puisque les diffrentes sources de chaleur, qui constituent le
vritable moteur du systme, sont susceptibles d'voluer la fois dans le temps et dans l'espace.
Soulignons que ce procd permet des conomies d'nergies substantielles: la stratification de l'air conduit un gradient
de temprature quasiment constant sur toute la hauteur de la pice. La temprature maintenir ne concernant que la
zone d'occupation (h <1,8 m), on peut se permettre des tempratures d'air repris plus leves, et diminuer ainsi la charge
frigorifique de l'installation. Ce dernier point est particulirement vrai pour des locaux industriels comprenant des
machines fortes dperditions thermiques.

44

(a)

(b)

Figure 10: Ventilation par diffusion (a) et par dplacement (b)


12

Dbit de soufflage

Il y a lieu de distinguer ici le dbit d'air neuf, fix par la rglementation, et le dbit de soufflage, qui doit permettre
d'extraire les charges thermiques et hydriques du local pour assurer, soit les conditions d'hygine (hpitaux, cuisines,..)
et/ou de confort (habitat, tertiaire), soit la qualit du produit fabriqu ou conditionn (micro-lectronique, agroalimentaire, pharmaceutique,...).
12.1
Dbit d'air neuf
Les deux tableaux suivants regroupent les dbits d'air neuf minimaux pour les locaux pollution non spcifique, en
m3/h et par occupant (tableau 4), et pour les locaux pollution spcifique, en m3/h (tableau 5). Toute rglementation
tant amene voluer, il faut bien entendu surveiller les dcrets relatifs aux rglements sanitaires et thermiques.
Soulignons ici, que les recommandations prconises par l'ASHRAE (ASHRAE Standard) ont fortement volu l'anne
dernire et il ne serait pas tonnant que le mme processus arrive en Europe au cours des prochaines annes.
12.2
Dbit de soufflage
La dtermination du dbit de soufflage ncessite le calcul des charges thermique et hydriques du local ou de l'ensemble
de plusieurs locaux, et le choix du systme de climatisation.
Dans le cas d'une installation centralise, le dbit de soufflage sera dtermin par rapport la droite de soufflage, et
requiert donc la connaissance de la temprature de soufflage, qui dpend elle-mme du type de bouche mis en place.
Dans le cas d'une climatisation localise (ventilo-convecteur, climatiseur), le dbit de soufflage peut se rduire au dbit
d'air neuf. Toutefois, on utilise en gnral une centrale de pr-conditionnement, en particulier lorsque l'humidit de l'air
extrieur est importante.
En France, il faut souligner que le dbit maximal d'air neuf ne doit pas excder 20% en zone H1 et H2, 30% en zone H3
du dbit total trait. En cas de dpassement de ces limites suprieures, il y a lieu de prvoir un systme de rcupration
thermique sur l'air extrait.
Enfin, la ventilation doit pouvoir tre arrte en cas de non-occupation ou de non-pollution des locaux. Ceci revient
souvent prvoir le comportement de l'installation centralis lorsqu' une ou plusieurs zones ne sont plus ventiles.

Tableau 4 : Locaux pollution non spcifique. Dbit minimal d'air neuf en m3/h et par occupant (air 1,2 kg/m3)
Destination des locaux

avec
interdiction de
fumer

Locaux d'enseignement:
Classes, salles d'tude, laboratoires ( l'exclusion de ceux pollution
spcifique:

Maternelles, primaires et secondaires du premier cycle

secondaires du 2me cycle et universitaires

Ateliers
Locaux d'hbergement:

Chambres collectives (plus de trois personnes*, dortoirs, cellules,


salles de repos

Bureaux et locaux assimils (tels que locaux d'accueil,


bibliothques, bureaux de poste, banques)
Locaux de runions:
45

sans
interdiction de
fumer

15
18
18

non fumeur
non fumeur
non fumeur

18

25

18

25

Salles de runions, de spectacles, de culte, clubs, foyers


18
30
Locaux de ventes

boutiques, supermarchs
22
non fumeur
Locaux de restauration:

Cafs, bars, restaurants, cantines, salles manger


22
30
Locaux usage sportif
par sportif:
22

Piscines
25

Autres
30
18
Par spectateur

: pour les chambres de moins de trois personnes, le dbit minimal prvoir est de 30 m3/h par local.

Tableau 5 : Locaux pollution spcifique. Dbit minimal d'air neuf en m3/h.


Pices usage individuel

Salle de bains ou de douche, cabinet d'aisances


15 par local
Pices usage collectif
30

cabinet d'aisances isol


45

Salle de bains ou de douches isoles


60

Salle de bains ou de douches commune avec un cabinet d'aisance


30+15N*

Bains, douche et cabinet d'aisance groups


10+15N*

Lavabos groups
5 par m2 de local**

Salle de lavage, schage et repassage du linge


Cuisines collectives
15/repas

Office relais
25/repas

Moins de 150 repas servis simultanment


20/repas, min: 3750 m3/h

De 151 500 repas servis simultanment


15/repas, min: 10000 m3/h

De 501 1500 repas servis simultanment


10/repas, min: 22500 m3/h

Plus de 1500 repas servis simultanment


* : N est le nombre d'quipement dans le local
**: Compte tenu des contraintes techniques, les dbits retenus seront de prfrence arrondis au multiple suprieur de
15

46

4me partie

Rseaux arauliques, dimensionnement

1. CHOIX DE LA VITESSE DECOULEMENT.


12
2. CALCUL DES PERTES DE CHARGE. 13
3. CALCUL DES RESEAUX
15
3.1.
Les diffrentes mthodes. 15
3.2. La mthode du regain en pression statique.
3.3. Choix du ventilateur 17

47

16

Le rseau araulique doit permettre dassurer les diffrents flux dair, issus de ltude des charges thermiques du
btiment et du taux de renouvellement dair. Une fois les dbits connus, il faut alors dterminer la totalit du
rseau araulique, cest dire ces diffrents lments (ventilateur, filtre, bouches de soufflage et daspiration )
et la manire de les connecter (formes et longueurs des gaines et des coudes, drivations).
Lquation gnrale utilise pour le calcul des rseaux est lquation de Bernoulli gnralise, applique une
conduite :
p t = p tA p tB = p A p B +

U 2A U 2B

+ g(z A z B )
2
2

(1)

Pour les rseaux arauliques, cette quation peut tre simplifie en introduisant les pressions relatives
p r = p p o , avec po: pression atmosphrique du lieu considr. On obtient:
p t = p rA p rB +

U 2A U 2B

+ g(z A z B ) + p o (z A ) p o (z B )
2
2

(2)

avec po(zi)
: pression atmosphrique laltitude zi.
Or, p o (z A ) p o (z B ) = air g(z A z B ) , et en considrant que la masse volumique de lair circulant dans la
conduite est peu diffrente de celle de lair extrieur, (2) se rduit :
p t = p rA p rB +

U 2A U 2B

2
2

(3)

Nota : Pour allger lcriture, les pressions relatives seront notes p dans tout ce qui suit.
Plusieurs mthodes de dimensionnement sont utilises qui diffrent dune part sur la mthode de calcul des
pertes de charge et dautre part sur la manire de dvelopper le rseau. L'objectif est d'obtenir un rseau
quilibr, c'est dire un rseau dont toutes les branches en parallle induisent la mme perte de charge, lorsque
les dbits souhaits circulent dans les diffrentes branches. Par exemple, le rseau reprsent figure 1 sera
& , V
& et V
& issues du cahier des charges on a:
quilibr si pour les diffrents dbits V
E
F
G
PAE = PAF = PAG
La branche racine (AC) tant commune aux trois dbits volumiques, on remarquera que cette conditions est
quivalentes :
PCE = PCF = PCG
A

.
V+ .
E V
.
F +V

C
.
VE

.
V+ .
F V

.
VG

F
.
VF

Figure 1 : exemple de rseau araulique trois branches drivs L'quilibrage impose l'galit des pertes de
charge entre AE, AF et AG
1.1.

Choix de la vitesse dcoulement.

Ce choix rsulte dun compromis entre le prix de linstallation et celui de son exploitation. A ces deux
contraintes principales, sajoutent plusieurs contraintes secondaires qui permettent dencadrer les valeurs
possibles :

48

grandes sections limites par des contraintes architecturales, et disolation thermique,


grandes vitesses (c'est dire petites sections) limites par le bruit engendr.
La figure 2 donnent, en fonction du dbit volumique, les vitesses recommandes par lASHRAE pour les
installations de conditionnement dair , ainsi que les perte de charge linaires quelles crent suivant le diamtre
de la conduite. (Cet abaque tant relatif au dbit volumique, le diamtre prendre en compte dans le cas dune
conduite non circulaire est le diamtre quivalent Deq).
Dans le cas o lacoustique est une contrainte primordiale (salle de spectacle, auditorium, chambre dhtel,
etc.), on ne dpassera en aucun cas 6 m/s.

1.2.

Calcul des pertes de charge.

On a vu dans la premire partie que les pertes de charges pouvait se mettre sous la forme:

U2
L
p = (Re, r )
+ i

Dh
2
i

L
i

(4)

: longueur total du tronon considr,


: coefficient de perte de charge singulire de laccident i.

Contrairement ce qui a t expos dans la partie 1, certains ouvrages ne caractrisent pas les pertes de charges
singulires par un coefficient , mais assimilent les accidents une longueur de canalisation de mme diamtre,
prsentant la mme perte de charge. Cette longueur de canalisation est dite longueur quivalente laccessoire
considr. Cette mthode de calcul simplifie grandement le calcul du rseau puisque chaque branche peut tre
caractrise par une simple longueur fictive, mais conduit des erreurs plus importantes que celles associes la
considration dun coefficient de perte de charge singulire .
La raison en est que le coefficient de forme dpend peu de la viscosit du fluide, et donc du nombre de
Reynolds, alors que la caractrisation de laccident par une longueur quivalente L:
p = (Re, r )

U2
L
2D

(5)

avec L : longueur quivalente de laccident,


entrane cette dpendance de Re par lintermdiaire du coefficient .
En consquence, lutilisation de cette mthode doit tre vite, ou du moins rserve aux avant-projets
sommaires.
Remarques :
Longueur totale L
La longueur totale L inclue non seulement les longueurs droites de canalisation, mais aussi les longueurs des
diffrents accessoires (coudes, largissement, rtrcissement, branchement ), que lon dtermine en prenant en
compte laxe de la conduite. Par exemple, un coude arrondi 90, de rayon de courbure R, aura deux
contributions dans le calcul des pertes de charge :

une contribution rgulire caractrise par une longueur l = R 2 ,

une contribution singulire caractrise par un coefficient .


On veillera ne pas confondre la longueur l, laquelle correspond une perte de charge rgulire (viscosit), avec
la longueur quivalente Leq dfinie ci-dessus, qui caractrise lensemble des pertes de charge rgulires et
singulires.
coefficient
Les tables ou corrlations donnant le coefficient prcisent toujours avec quelle vitesse il faut calculer la
pression dynamique. Dans le cas dune branche section variable, ou comportant une drivation, il faudra donc
rapporter ce coefficient la vitesse prise en compte dans lquation (3), qui est en gnral la vitesse dans le
tronon le plus en amont. Par exemple, la figure 4 montre une branche comportant une drivation en A et un
coude en B. Cette branche comporte donc trois tronons I,II et III. Soit A et B les coefficients relatifs
respectivement la drivation A et au coude B. On a, pour les pertes de charge singulires dues ces deux
accessoires: p A = A

UI2
2

p B = B

U II 2
2

49

Figure 2 : vitesses dair et pertes de charge linaires admissibles en fonction du dbit volumique (daprs
ASHRAE).

50

&I
V

II

& II
V

III

&
V
III

Figure 3 : branche de rseau comportant une drivation A, et un coude B.


Pour pouvoir se ramener la forme (4), il faut exprimer A en fonction de UII. On a:
2

U
U 2
U 2
p A = A I = B I II
2
2
U II
ce qui permet de dduire la valeur de *B , coefficient de perte de charge de la drivation rapport UII:
2

U
*A = A I
U II
et dexprimer la perte de charge singulire due la drivation sous la forme :
U 2
p A = *A II
2
La perte de charge totale du tronon (AB) devient alors analogue (4).:

U2
L
p AB = (Re, r ) II + B + *A II
D II
2

1.3.

Calcul des rseaux

Le dimensionnement dun rseau araulique ncessite tout dabord la connaissance des diffrents flux dair
mettre en uvre, puis les caractristiques et lemplacement des diffrents accessoires (bouches dair, filtres,
coude, etc... ).
On peut alors dessiner le rseau araulique, en utilisant une reprsentation unifilaire en perspective
axonomtrique. Ce schma doit comporter tous les accidents (coudes, branchements) auxquels on attribuera
un repre. Seront indiqus le long de chaque branche le dbit qui la traverse, et, sur chaque partie rectiligne, la
longueur sparant deux accidents.
L encore, de nombreuses mthodes existent dont les plus connues sont :

le dimensionnement perte de charge linaire constante,

le dimensionnement vitesse constante,

le dimensionnement regain de pression statique,

loptimisation technico-conomique.
1.3.1. Les diffrentes mthodes.
Le dimensionnement perte de charge linaire constante est particulirement bien adapt lorsque les pertes de
charge singulires sont caractrises par des longueurs quivalentes aux accidents. Elle consiste dterminer le
diamtre des diffrents tronons de faon avoir une perte de charge linaire constante dans tout le rseau. Le
choix de cette perte de charge peut tre effectue partir de la figure 2 mais dpend des contraintes technicoconomiques: lorsque le cot de lnergie (cot de fonctionnement) est lev et le cot des quipements (cot
dinvestissement) faible, un faible taux de perte de charge linaire sera plus conomique. Inversement, si le cot
de lnergie est faible et le cot des quipements importants , on pourra choisir un taux de pertes de charge plus
important.
Le dimensionnement vitesse constante est galement possible. La procdure consiste dimensionner les
diffrentes sections de passage de faon avoir une vitesse constante dans toute linstallation.
Dans les deux cas, il est ncessaire, aprs un premier dimensionnement initial, deffectuer le calcul des pertes de
charge totale dans chaque branche du rseau et de modifier certains diamtres de faon quilibrer le rseau,
cest dire dobtenir les dbits dair voulus dans chaque branche. En gnral, on garde le critre de dpart
(vitesse ou perte de charge constante) dans la branche la plus rsistante, et on adapte les pertes de charge des

51

branches secondaires de faon assurer la contrainte de dbits, soit en modifiant les diamtres, soit en intgrant
des organes de rglage induisant des pertes de charge supplmentaires.
Lavantage de ces deux mthodes est quelles sont simples et rapides, mais leur inconvnient rside dans le fait
quelle ne correspondent aucun optimum puisque lquilibrage du rseau reste artificiel : on provoque une
dissipation supplmentaire dnergie dans les branches qui conduisent un dbit trop important. Cest pour cette
raison quon leur prfre aujourdhui la mthode dite de regain en pression statique qui ralise loptimisation
nergtique du rseau et qui est expos au paragraphe suivant.
Dautres mthodes ont t rcemment dveloppes. Elle consiste en une procdure doptimisation qui minimise
une fonction objectif. Cette fonction objectif prend en compte non seulement le cot nergtique du
fonctionnement de linstallation, mais intgre galement le cot dinvestissement initial, les taux dinflation et
les taux dintrt. On peut citer par exemple la mthode dveloppe par Tsal et Coll., appele T-Method, qui
permet de remplacer les diffrentes branches du rseau par un seul tronon quivalent, caractris par un
diamtre fictif. Le calcul des pertes de charge et donc la slection du ventilateur est alors possible. Cette mthode
permet ainsi doptimiser lensemble de linstallation, ventilateur compris, dun point de vue non seulement
nergtique mais aussi conomique. Elle reste malheureusement relativement complexe puisquelle fait appel
plusieurs comptences : calcul araulique, calcul conomique et analyse numrique. Les tudiants intresss
pourront se rfrer directement la publication de Tsal sous la rfrence suivante :
Tsal R.J.,H.F.Behls, and R. Mangel; T-method for duct design. Part 1 : Optimization Theory; Part 2 : Calculation
procedure and economic analysis. ASHRAE Transaction 94(2), pp90-111.
1.3.2. La mthode du regain en pression statique.
Lobjectif de cette mthode est dobtenir une pression statique quasi-constante chaque nud du rseau (les
nuds dun rseau sont les points o apparaissent des drivations ou des jonctions), pour dfinir le mieux
possible le dbit des drivations successives. Ceci est effectu en adaptant les diamtres des tronons qui relient
les diffrents nuds de faon ce que la perte de charge du tronon soit gale la diffrence de pression
dynamique en amont des deux nuds considrs (figure 4).
(3) montre que pour avoir galit des pressions statiques en A et B, il faut vrifier :
U 2A U 2B

p t =

(6)
2
2

La dtermination de la perte de charge pt ncessitant la connaissance de UB, la recherche du diamtre vrifiant


(6) se fait par itrations.

Pt
ps
PdI

pl

pt

PdII

&
V
II

&
V
I
A

&
V
1

&
V
2

Figure 4 : Regain en pression statique. Le diamtre du tronon reliant les drivations 1 et 2 est tel que la chute de
pression dynamique compense la perte de charge entre les deux sections amont A et B des deux drivations. ps
reprsente la perte de charge singulire due la culotte de drivation, pl correspond la perte de charge linaire
entre A et B.
Il faut noter que (6) ne peut tre vrifie qui si la vitesse diminue aprs chaque drivation. Ceci peut impliquer
une augmentation du diamtre aval. Si cette solution nest pas envisageable, on garde le mme diamtre, et on
transforme lquation (3) pour calculer la pression statique au nud B :

52

p rA p rB = p t

U 2A U 2B
+
2
2

(7)

Pour commencer la procdure de dimensionnement, on slectionne une vitesse maximum pour la branche racine.
On pourra l encore se reporter la figure 2, et on retrouve le mme critre conomique que dans le cas dun
dimensionnement perte de charge constante mais exprim cette fois ci en vitesse: lorsque le cot de lnergie
est lev et le cot des quipements faible, une faible vitesse sera plus conomique. Inversement, si le cot de
lnergie est faible et le cot des quipements importants, on pourra choisir une vitesse plus importante. Pour
toutes les autres branches, exceptes les branches terminales, la vitesse est dtermine par itration laide de
lquation (6). Pour les branches terminales, la pression statique est impose par les caractristiques des bouches
dair. La connaissance des pressions statiques en amont de ces bouches, permet, par lutilisation de (6), de
dterminer le diamtres des branches terminales. Dans le cas dune sortie direct (conduit donnant directement
dans le local), la pression statique terminale est bien videmment 0 (pression relative).
Lintrt de cette mthode est quelle prend en compte ds le dpart lquilibrage du rseau: les dbits dans les
diffrentes branches seront donc bien dfinis. Dautre part, on voit que le calcul du rseau conduit
intrinsquement lutilisation optimale de la pression totale dlivre par le ventilateur.
Son inconvnient est quelle ncessite des calculs itratifs sur chaque branche du rseau qui deviennent vite
fastidieux lorsque le rseau est complexe. Lutilisation de linformatique devient alors vite recommande.
1.3.3.

Choix du ventilateur

&),
Une fois le rseau dfini, on peut dterminer la caractristique du rseau (appele aussi ouverture) p = f (V
& dbit volumique dans la branche racine (ou dbit principal):
avec V
&2
p = K V
Cette caractristique a une allure parabolique, avec toutefois une partie linaire pour les petits dbits,
correspondant au rgime laminaire. Cette zone ntant jamais atteinte en fonctionnement, on peut lignorer. Pour
trouver le point de fonctionnement de linstallation, il suffit de reporter cette caractristique sur les
caractristiques des ventilateurs fournies par les constructeurs. Le point de fonctionnement sera dfini par
lintersection des deux courbes (figure 5).
Pt ou ht

&
V
Figure 5 : Point de fonctionnement
& , on peut calculer la puissance consomme par le ventilateur :
Connaissant le point de fonctionnement pt, V
&
& = p t V
W

53

Vous aimerez peut-être aussi