Thèse de Jérôme Giordano soutenue en 2004 - iusti
Thèse de Jérôme Giordano soutenue en 2004 - iusti
Thèse de Jérôme Giordano soutenue en 2004 - iusti
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2.1. Modèles physiques <strong>de</strong> la dynamique <strong>de</strong>s flui<strong>de</strong>s compressibles<br />
2.1.2 Expression <strong>de</strong>s termes sources chimiques<br />
L’étu<strong>de</strong> d’ecoulem<strong>en</strong>ts à haute <strong>en</strong>thalpie génératrice impose <strong>de</strong> ne pas préjuger <strong>de</strong> l’état<br />
thermo-chimique du gaz. Ainsi les équations <strong>de</strong> relaxation chimique et vibrationnelle du gaz<br />
doiv<strong>en</strong>t être modélisées afin <strong>de</strong> r<strong>en</strong>dre compte d’un év<strong>en</strong>tuel état à l’équilibre, figé et hors d’équilibre<br />
[81]. Considérons une réaction chimique écrite sous sa forme générique (réactions chimiques<br />
<strong>de</strong> l’air à haute température, combustion, etc.) :<br />
∑Nsc<br />
ν k ′ A k ⇐⇒<br />
k=1<br />
∑Nsc<br />
k=1<br />
ν ′′<br />
k A k (2.6)<br />
où ν<br />
k ′ et ν′′<br />
k<br />
représ<strong>en</strong>t<strong>en</strong>t, respectivem<strong>en</strong>t, les coeffici<strong>en</strong>ts stoechiométriques directs et inverses<br />
<strong>de</strong> chaque produits A k ; Nsc étant le nombre total d’espèces qui peuv<strong>en</strong>t réagir. Le taux <strong>de</strong><br />
production molaire <strong>de</strong> l’espèce k se calcule <strong>de</strong> la manière suivante :<br />
dn k<br />
dt<br />
= (ν ′′<br />
k − ν′ k )(K d<br />
∏<br />
l<br />
n ν′ l<br />
l<br />
− K i<br />
∏<br />
l<br />
n ν′′ l<br />
l<br />
) (2.7)<br />
où k d et K i sont les vitesses <strong>de</strong> la réaction directe et inverse et n k = ρ k /M k représ<strong>en</strong>te le nombre<br />
<strong>de</strong> môles <strong>de</strong> l’espèce k par unité <strong>de</strong> volume. Ainsi, pour obt<strong>en</strong>ir le taux <strong>de</strong> production massique<br />
total, on additionne le terme dρ k<br />
dt<br />
sur les N r réactions que conti<strong>en</strong>t le modèle réactif choisi, soit :<br />
∑N r<br />
( ) dnk<br />
ω k = M k<br />
dt<br />
r<br />
r=1<br />
(2.8)<br />
Les vitesses <strong>de</strong> réactions directe K d,r et inverse K i,r repos<strong>en</strong>t sur <strong>de</strong>s lois dép<strong>en</strong>dant ess<strong>en</strong>tiellem<strong>en</strong>t<br />
<strong>de</strong> la température <strong>de</strong> la translation T du mélange. La vitesse directe s’écrit :<br />
K d,r = A r T Br exp( −θ r<br />
T ) (2.9)<br />
où les constantes A r , B r et θ r sont respectivem<strong>en</strong>t le facteur pré-expon<strong>en</strong>tiel, l’exposant préexpon<strong>en</strong>tiel<br />
et la température d’activation <strong>de</strong> la réaction numéro r. Ces constantes <strong>de</strong> réaction<br />
sont généralem<strong>en</strong>t déduites par un ajustem<strong>en</strong>t à <strong>de</strong>s mesures expérim<strong>en</strong>tales [81]. Les vitesses<br />
<strong>de</strong> réaction inverse K b,r sont, suivant le modèle, soit données <strong>de</strong> la même manière, soit liées à la<br />
vitesse directe au moy<strong>en</strong> la constante d’équilibre K eq,r par la relation :<br />
2.1.3 Modèle réactionnel<br />
K i,r = K d,r<br />
K eq,r<br />
(2.10)<br />
Une réaction chimique est généralem<strong>en</strong>t composée <strong>de</strong> nombreuses réactions élém<strong>en</strong>taires faisant<br />
interv<strong>en</strong>ir <strong>de</strong>s intermédiaires réactionnels. Lorsque l’on veut simuler cette réaction, il faut<br />
<strong>en</strong> fait le faire pour l’<strong>en</strong>semble <strong>de</strong>s réactions élém<strong>en</strong>taires [76]. Mais le plus souv<strong>en</strong>t on peut<br />
13