Modélisation de l'évaporation de gouttes multi-composants

ANNEXE A : PROPRIETES PHYSIQUESB =∞∞ 2∞ 3( bpv00+ bpv01P+ bpv02P+ bpv03P)∞∞ 2∞ 3+ ( bpv10+ bpv11P+ bpv12P+ bpv13P) T∞∞ 2∞ 3+ ( b) 2pv20+ bpv21P+ bpv22P+ bpv23PT(A.53)Si la variable de distribution I est la température normale d’ébullition, elle s’écrit :P( I ) P [ ( s RT )( T T ( I ))],= ∞ exp τ −(A.54)sat jj fgb,javec T ( I )bIj=∞ln( P )(A.55),1−s Rfget le coefficient correcteur τ = τ + τ∞ + τ P (A.56)jab P c∞ 2Température normale d’ébullition :Température critique :T= T T I(A.57)nb ba+Tbb= T T I(A.58)c ca+cbEnthalpie de vaporisation :⎛ Tca+ TcbI− T ⎞l = ( + )⎜⎟vlvbalvbbI(A.59)⎝ Tc− Tnb⎠0.38Densité molaire du liquide :c( c + c T ) + ( c + c T ) I + ( c c T ) 2= (A.60)l lA1 lA2lB1lB2lC1+lC 2ILa densité molaire du liquide c l est une corrélation du second ordre en fonction de la variablede distribution I car il faut une très bonne précision sur cette propriété. En effet, elle estutilisée pour calculer la décroissance de la taille de la goutte et une erreur peut entraîner unedifférence de 10% sur sa durée de vie en utilisant une corrélation du premier ordre.Chaleur spécifique du liquide :CpViscosité dynamique du liquide :log( Cp + Cp T ) + ( Cp Cp T )I= (A.61)l lA1 lA2lB1+lB2⎛ 12( μ ) ( ) ⎜⎟ l= Visb1+ Visb2I+ Visb3I−2⎝ T Vist01+ Vist02I + Vist03I ⎠1⎞(A.62)239