Combustion Turbuulente PDF

C OMBUSTION TURBULENTE DE SPRAYS Combustion
MASTER 2
Auteur: Julien Reveillon de Rouen Universite
de Rouen - [email protected] p.1/55 J. Reveillon - Universite
Introduction - Gnralits
Comme nous lavons voqu dans lintroduction gnrale, une description prcise de la dispersion et de lvaporation du spray associe des modles ddis au micro-mlange de la vapeur est ncessaire an de modliser correctement les phnomnes de combustion au sein de la chambre. En fonction des caractristiques du systme tudi, la combustion peut avoir lieu soit aprs lvaporation complte du combustible, soit au sein du spray en cours dvaporation.
Dans le premier cas, mme si lutilisation des modles de combustion dvelopps dans le cadre dcoulements purement gazeux est envisageable, la topologie de la fraction de mlange issue de lvaporation peut tre trs diffrente dune injection gazeuse classique et de nombreuses zones partiellement prmlanges peuvent apparatre. Si des ammes se propagent dans le milieu alors que le spray nest pas encore totalement vapor, des interactions trs complexes apparaissent entre le spray et son vaporation, la turbulence et les phnomnes de combustion.
dispersion evaporation coalescence breakup wall interactions
SPRAY
models ?
mixing dissipation
TURBULENCE
COMBUSTION
height of ignition flame length combustion regimes extinction flashback
le spray modie localement les transferts de chaleur et de
quantit de mouvement ainsi que la fraction de mlange et le
taux de dissipation. Tous ces facteurs affectent la structure
de la amme mais aussi les caractristiques de la turbulence.
Richesse du mlange
Si lon tudie la combustion des coulements purement gazeux, la dnition de la richesse locale du mlange (et donc de la fraction de mlange) ne pose pas de difcult quel que soit le rgime de combustion prmlang ou partiellement prmlang. En fonction de la topologie du champ de fraction de mlange, il est possible de prvoir avec une assez grande prcision le rgime local de combustion et le taux de raction en fonction de la quantit de combustible inject dans la chambre.
Richesse du mlange
Dans le cadre de la combustion diphasique, connatre le ux massique de combustible inject dans la chambre nest pas sufsant pour la dtermination des richesses locales de combustion.
En effet, aprs sa pulvrisation, le combustible liquide se dplace au sein de loxydant et peut traverser des zones dont les proprits de temprature et de pression vont inuer sur (1) le taux dvaporation des gouttes (2) la densit locale doxydant. Ce mme phnomne se rencontre lorsque les gouttelettes sapprochent dun front de amme.
Richesse du mlange
Considrons une goutte isole dans un volume de contrle V contenant de loxydant gazeux. Le diamtre de cette goutte dtermine la richesse liquide initiale du problme par rapport la quantit doxydant au sein du volume V . Loxydant entourant la goutte subit une lvation de temprature dun facteur 5 ce qui correspond un ratio classique: temprature gaz brl / temprature gaz frais.
Richesse du mlange
c
5 4 3
acements2
1
Temps
G /L u /b
time
v
Richesse du mlange
Deux paramtres entrent en jeu dans cette analyse: les temps caractristiques dvaporation v de la goutte et dlvation de la temprature ou temps de chauffage c . Le domaine V ntant pas ferm, la pression demeure gale la pression atmosphrique et une augmentation de la temprature du mlange doxydant et de combustible vapor provoquera son expansion. De ce fait, la quantit doxydant au sein du volume de contrle diminuera.
Pendant ce temps, la goutte svapore dans ce volume V . Il savre donc que la richesse gazeuse nale dpend du paramtre Dv qui est le rapport entre le temps dvaporation v et le temps de chauffage c : de Rouen - [email protected] p.9/55 J. Reveillon - Universite Dv = v /c .
Richesse du mlange
Un temps dvaporation trs bref (Dv trs petit) provoque la disparition de la gouttelette avant la diminution de la masse doxydant et la richesse gazeuse est trs proche de la richesse liquide initiale. En revanche, si Dv est grand, la perte de masse doxydant implique une forte augmentation de la richesse gazeuse locale. Le ratio maximal thorique de richesse peut donc atteindre le ratio de temprature gaz brls/gaz frais. Dun point de vu pratique le paramtre Dv reste petit ce qui correspond nanmoins la zone de forte variation du rapport de richesse.
Richesse du mlange
u /b
g replacements
v c
2
Temps
1 0 10 Dv 20
G /L
Flammes sprays laminaires

La seconde tape ltude de la combustion diphasique est lanalyse de la propagation et de la structure des ammes prmlanges diphasiques dans un coulement compressible.
En gnral, deux hypothses sont adoptes: soient toutes les gouttes svaporent avant datteindre le front de amme, soit chaque goutte brle avec une amme de diffusion autour delle. Par la suite, de nombreuses tudes ont t consacres aux cas monodisperss ou polydisperss mais trs souvent avec une prsence initiale de vapeur de combustible et des gouttes qui modient faiblement la richesse locale.

Cependant, dans ces tudes analytiques, la dilatation des gaz nest pas prise en compte or nous avons pu voir que cette dilatation peut provoquer de fortes modications des proprits locales du mlange combustible/oxydant et donc de la structure et de la vitesse de propagation des ammes . Grce la DNS, la propagation des ammes diphasiques monodimensionelles au sein dun coulement compressible peut tre tudie. Une premire conguration simple a t mise en uvre: des gouttes monodisperses sont rgulirement positionnes dans notre domaine de calcul au sein dun coulement initialement au repos.

Gaseous fuel 1. Temperature
computational domain
Gaseous fuel Droplets 2.
Temperature
Droplets 3.
Temperature Gaseous fuel
Index de amme
An de pouvoir diffrencier le dgagement de chaleur du une combustion non-prmlange par rapport une combustion prmlange, Takeno a propos dutiliser un index de amme bas sur le produit scalaire des gradients de fraction de mlange de combustible et doxydant. Lindex de amme de Takeno, GF O , scrit:
G F O = Y F Y O
Il est positif pour une amme prmlange et ngatif pour une amme non-prmlange.
Index de amme
A partir de cette dnition de base, un index de amme normalis p peut tre dni tel quil soit unitaire pour une combustion purement prmlange et nulle pour une combustion non-prmlange. Soit
1 p = 2 GF O 1+ |GF O |
Bien sur, la dtermination du rgime de combustion na lieu quaux points o le taux de raction est non-nul.

1.2 B 1
0.06 0.08
0.8 Wemax
0.6
Wvmax
0.04
0.4
0.02
0.2 A 0 325 C335 E 345 Time 355 365

0 325 335 345 Time 355 365
Rgime puls priodique

74
0.08
72 70 68 66
X(Flam) X(Evap)
0.06 B
62 60 58 56 54 52 50 325
Dire
Wvmax
64
0.04
C 0.02
ction de p ropa gatio
A E D 0 0 0.2 0.4 0.6 Wemax 0.8 1 1.2
335
345 Time
355
365
Ces pulsations ne proviennent donc pas de lvaporation goutte goutte de la phase liquide du domaine mais dune dynamique plus gnrale de la amme.
Lintensit du taux de raction pour les deux pics est trs proche de lunit ce qui correspond au taux de raction maximal de la amme stchiomtrique gazeuse.
front de flamme front devaporation
Premier front de flamme associe au front devaporation
Un second front de flamme suit le premier
Diagramme de combustion
Diagramme de Chiu : nuage de gouttes plong au sein dun gaz chaud au repos. Les rgimes de combustion ont ainsi t classis en fonction dun nombre de groupe dnot G. Ce nombre sans dimension peut tre vu comme le rapport entre la vitesse caractristique dvaporation et la vitesse de diffusion de la temprature au sein du nuage.
G 5N 2/3 /S dpend du nombre de groupe G, le nombre total N de gouttes dans le nuage et le paramtre de sparation S . Ce dernier est en fait le rapport S = s /rf entre la distance moyenne intergoutte s et le rayon rf
de la amme de diffusion qui entourerait une unique goutte ayant les proprits moyennes des gouttes du de Rouen - [email protected] p.21/55 J. Reveillon - Universite nuage.
De ce fait, lorsque le nombre de sparation dcrot, il existe un point o la topologie de la amme passe dun rgime de combustion individuel un rgime de groupe.
En xant S et en faisant varier N , deux modes principaux pour la combustion de spray peuvent tre identis. Le premier mode est quali dexterne, le second mode dinterne.
Dans le premier cas, si G >> 1, les gouttes sont trop proches les unes des autres pour permettre la diffusion de temprature au sein du nuage et seule une couche externe de liquide est vapore. La amme reste donc une distance xe du bord du nuage. Dans le second cas, si G << 1, les gouttes sont sufsamment loignes les unes des autres pour que les gaz chauds pntrent dans le nuage. Alors un processus individuel dvaporation et de combustion a lieu autour de chaque goutte.
External sheath combustion External group combustion External flame
nonvaporizing droplets G>>1 G<1 main flame vaporizing droplets G>1 G<<1
Individual droplet burning Internal group combustion single droplet combustion
Dans un deuxime temps, Chang puis Borghi ont ajout cette analyse la possibilit de prendre en compte les proprits des ammes. Ils ont donc introduit lpaisseur de amme f et son temps caractristique f quils ont associs au temps dvaporation moyen v des gouttes.
Daprs le diagramme quils ont pu driver, lorsque v << f , les gouttes svaporent trs rapidement et une combustion de type prmlange a lieu quelles que soient la distance intergouttes ou lpaisseur de amme.
Si le temps dvaporation est plus important, pour f > s , les gouttes pntrent dans les zones de raction/diffusion car la amme est plus paisse que lespace inter-goutte s . Cette situation mne rapidement un paississement de la amme.
Aprs le passage dun premier front partiellement prmlang, des zones secondaires de raction peuvent demeurer. La topologie de cette zone secondaire dpendra du paramtre S . Pour les petites valeur de S , les gouttes brlent individuellement ou alors sont amalgames en petits groupes (ou clusters) entours par une amme. Il sagit du rgime de combustion de groupe. Lorsque la valeur du paramtre S augmente, Borghi et al. ont aussi discern lapparition dun rgime de percolation et dun rgime de poches
ln(tv/tf) S > 0.73 Pockets Thick flame
0.73>S>0.41 Percolation S < 0.41 Groups
0 Prevaporized flame 0 ln(ds/df)
prevaporization zone reaction zone
Fuel rich premixed reaction zone
diffusion reaction zone
Fuel lean premixed reaction zone
Hot Products + Oxidizer
flame zone
(left) Prevaporized spray ame (middle) Prevaporized spray ame with sparse droplets (right) Thick ame
Thick reaction zone
(left) Pocket (middle) Percolation (right) Group
A la fois Chiu et al et Borghi et al ont dni ces structures de amme avec deux restrictions majeures: lcoulement est initialement au repos la notion de richesse fuel/oxydant na pas t abord Pourtant, dans un systme de combustion de spray cette richesse est un point clef capable dinuencer les structures de amme, la stabilit de combustion ainsi que lensemble des phnomnes physiques prenant place au sein de la chambre.
A la fois Chiu et al et Borghi et al ont dni ces structures de amme avec deux restrictions majeures: lcoulement est initialement au repos la notion de richesse fuel/oxydant na pas t abord Pourtant, dans un systme de combustion de spray cette richesse est un point clef capable dinuencer les structures de amme, la stabilit de combustion ainsi que lensemble des phnomnes physiques prenant place au sein de la chambre.
Flamme spray "prmlange"
Flamme spray "non-prmlange"

1 1.5 0.9
Premixed
1
0.8 0.7 0.6 0.5
0.5 0.4
Ox.
0.3 0.2
Diffusion Fuel*15
0 0.5 1 1.5
0.1 0
Distance
Flamme double
YF*15 and YO
Reaction rate

0.14 0.13 0.12 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3
Reacting rate, diffusion regime
0.11 0.1 0.09 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0 0 1 2 3 4
Fuel
Oxidizer
Fuel
0.2 0.1 0
Distance
Flamme de diffusion centrale
YF*15, YO

Fuel
0.2
prem.
1.8
prem.
0.1
prem.
1.6 1.4 1.2
Reaction rate
diff.
-0.1
1 0.8
-0.2
diff.
0.6
Oxidizer
-0.3 0.4 0.2 0
-0.4 0 1 2 3
Distance
Flamme de prmlange centrale
Y1*15, Y2

Short evaporation delay No combustion Interspace Medium evaporation delay No combustion Interspace
Flammability limits
Flammability limits
Closed external combustion

ER=0.87
Closed external combustion 1

ER=0.87 ER=1 ER=1.13 ER=0.87 ER=1 ER=1.13 ER=1.48
1
ER=1 ER=1.13 ER=1.48
Group combustion (vapor pocket)
Group combustion
ER=1.48
Thick flame limit
ER=1.85
ER=1.85
ER=1.85
Hybrid combustion
Hybrid combustion Central diffusion flame

ER=2.65 ER=2.65
ER=2.65 d=248
Central premixed flame

d=41 d=25
Open external combustion Short evaporation delay limit Equivalence ratio Equivalence ratio Open external combustion

Evaporation delay
Interspace Thick flame limit

Group combustion
Hybrid combustion Cent. Diff.
Cent. Prem
Equivalence ratio
prevaporized flame
Interspace
Closed external combustion 1
it
Fl
am
ab
ili
ty
lim
Group combustion
Hybrid combustion
Open external flame
Equivalence ratio
Gnralits prmlange/diffusion
YO 1 lean side P<1, Z<Zst P=1, Z=Zst R Fig. A Fig. B YO Species through the premixed flame, lean case P=cte
RMAX
P>1, Z>Zst PS rich side RMIN YF isoP
isop p=0.66 PS 0 1 YF p=0.83 p=1
Taux de combustion prmlange

Dnissons Wp (x), la quantit de combustible brlant en rgime prmlang.
p F dy Wp (x) =
Ly
F dy
Ly
Grce ce paramtre, il est possible de voir, en effectuant les statstiques selon la direction perpendiculaire lcoulement, que la amme est tout dabord domine par le prmlange sa base.
Taux de combustion prmlange

1
0.8
0.8
0.6 PBR
0.6 PBR
0.4
0.4
0.2
0.2
ements
0 0 10
PSfrag replacements
g x/f
20 30 40 50
10
20
g x/f
30
40
50
Gaseous ame
Spray ame
Dcomposition terme source

Les variables de progrs ont t utilises dans le cadre des modlisations de type amelette. Les ammes prmlange peuvent tre dcrites par la variable c(x, t) qui est gale 1 dans les gaz brls et 0 dans les gaz frais. Des modles comme le BML (cf cours comustion turbulente) utilise la variable de progrs dans le cas de rgimes de type amelettes.

Un rsultat classique de la combustion gazeuse scrit
ci + u ci = wci |ci | t
avec la vitesse de lisosurface

1 w ci = |ci | +
N
m ci 2 (Dci ) + ci + d
k=1 N
Aik Z,ck
2
k=1
Aik
ck Z
2 ci Z 2

Le taux de consommation de combustible en rgime prmlang scrit:
p Eq F = (wc |c| (Dc))(YF,0 YF (Z ))
De mme dans le cas de ammes de diffusion, il est possible dapplique la relation

d dYF d Z 2 = F (Z, Z ) dZ
La contribution locale complte de la combustion prmlange et non prmlange peut tre modlise par la relation associant les deux taux de raction issues de lanalyse de type ammelettes.
p d F = p F (c, Z ) + (1 p ) F (Z, Z )
A partir de ce point nimporte quel modle ddi la combustion gazeuse prmlange et non-prmlange peut tre envisag. Le modle de type ammelette est un exemple.

0
0.0
0.2
1
0.4
0.6
0.8
3
1.0
ements
4 0 10
PSfrag replacements
g x/f
20
1.2
30
40
50
10
20
g x/f
30
40
50
Gaseous ame Spray ame Comparaisons modle/DNS pour le taux de combustion prmlange.

0.0
0.1
0.2
0.2
0.4
0.3
0.6
0.4
ements
0.5 0 10
PSfrag replacements
g x/f
20 30 40 50
0.8
10
20
g x/f
30
40
50
Gaseous ame Spray ame Comparaison modle/DNS pour le taux de combustion non-prmlange.

Combustion Turbuulente PDF

Transféré par

Droits d'auteur :

Formats disponibles

Combustion Turbuulente PDF

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Combustion Turbuulente PDF

Transféré par

Droits d'auteur :

Formats disponibles

C OMBUSTION TURBULENTE DE SPRAYS Combustion

de Rouen - [email protected] p.1/55 J. Reveillon - Universite

de Rouen - [email protected] p.2/55 J. Reveillon - Universite

de Rouen - [email protected] p.3/55 J. Reveillon - Universite

height of ignition flame length combustion regimes extinction flashback

le spray modie localement les transferts de chaleur et de

quantit de mouvement ainsi que la fraction de mlange et le

taux de dissipation. Tous ces facteurs affectent la structure

de la amme mais aussi les caractristiques de la turbulence.

de Rouen - [email protected] p.4/55 J. Reveillon - Universite

de Rouen - [email protected] p.5/55 J. Reveillon - Universite

de Rouen - [email protected] p.6/55 J. Reveillon - Universite

de Rouen - [email protected] p.7/55 J. Reveillon - Universite

de Rouen - [email protected] p.8/55 J. Reveillon - Universite

de Rouen - [email protected] p.10/55 J. Reveillon - Universite

de Rouen - [email protected] p.11/55 J. Reveillon - Universite

Flammes sprays laminaires

de Rouen - [email protected] p.12/55 J. Reveillon - Universite

Flammes sprays laminaires

de Rouen - [email protected] p.13/55 J. Reveillon - Universite

Flammes sprays laminaires

Gaseous fuel Droplets 2.

Temperature Gaseous fuel

de Rouen - [email protected] p.14/55 J. Reveillon - Universite

de Rouen - [email protected] p.15/55 J. Reveillon - Universite

de Rouen - [email protected] p.16/55 J. Reveillon - Universite

Flammes sprays laminaires

0.2 A 0 325 C335 E 345 Time 355 365

Rgime puls priodique

de Rouen - [email protected] p.17/55 J. Reveillon - Universite

Flammes sprays laminaires

A E D 0 0 0.2 0.4 0.6 Wemax 0.8 1 1.2

de Rouen - [email protected] p.18/55 J. Reveillon - Universite

Flammes sprays laminaires

de Rouen - [email protected] p.19/55 J. Reveillon - Universite

Flammes sprays laminaires

front de flamme front devaporation

Premier front de flamme associe au front devaporation

Un second front de flamme suit le premier

de Rouen - [email protected] p.20/55 J. Reveillon - Universite

de Rouen - [email protected] p.22/55 J. Reveillon - Universite

de Rouen - [email protected] p.23/55 J. Reveillon - Universite

Individual droplet burning Internal group combustion single droplet combustion

de Rouen - [email protected] p.24/55 J. Reveillon - Universite

de Rouen - [email protected] p.25/55 J. Reveillon - Universite

de Rouen - [email protected] p.26/55 J. Reveillon - Universite

de Rouen - [email protected] p.27/55 J. Reveillon - Universite

0.73>S>0.41 Percolation S < 0.41 Groups

0 Prevaporized flame 0 ln(ds/df)

de Rouen - [email protected] p.28/55 J. Reveillon - Universite

Fuel rich premixed reaction zone

diffusion reaction zone

Fuel lean premixed reaction zone

Hot Products + Oxidizer

Thick reaction zone

de Rouen - [email protected] p.29/55 J. Reveillon - Universite

(left) Pocket (middle) Percolation (right) Group

de Rouen - [email protected] p.30/55 J. Reveillon - Universite

de Rouen - [email protected] p.31/55 J. Reveillon - Universite