Asymptotique des grandes énergies d'activation

L'asymptotique des grandes énergies d'activation est un développement asymptotique utilisée dans le domaine de la combustion. Il utilise le fait que la vitesse de réaction est extrêmement sensible aux changements de température en raison de la grande énergie d'activation de la réaction chimique.

Histoire

La technique a été introduite par Iakov Zeldovitch, David Frank-Kamenetskii et leurs collègues dans les années 30, dans leur étude sur les flammes prémélangées et les déflagrations (théorie de Frank-Kamenetskii). Elle a été utilisée dans l'ensemble de la communauté scientifique au début des années 70, grâce aux travaux de Williams B. Bush, Francis E. Fendell, Forman Williams, Amable LiñánModèle:, et John Frederick ClarkeModèle:,. Par la suite elle a été largement utilisée pour expliquer des problèmes plus complexes de combustion.

Présentation de la méthode

Dans les processus de combustion, la vitesse de réaction $\omega$ dépend de la température $T$ sous la forme suivante (loi d'Arrhenius),

$\omega (T)\propto {\mathrm{e}}^{-E_{\mathrm{a}}/RT},$

où $E_{\mathrm{a}}$ est l'énergie d'activation, et $R$ est la constante universelle des gaz parfaits. En général, la condition $E_{\mathrm{a}}/R\gg T_b$ est satisfaite, où $T_{\mathrm{b}}$ est la température des gaz brûlés. Cette condition constitue la base de l'asymptotique de l'énergie d'activation. En désignant $T_{\mathrm{u}}$ comme température des gaz imbrûlés, on peut définir le nombre de Zeldovitch et le paramètre de dégagement de chaleur comme suit :

$\beta =\frac{E_{\mathrm{a}}}{R{T}_{\mathrm{b}}}\frac{T_{\mathrm{b}}-T_{\mathrm{u}}}{T_{\mathrm{b}}},q=\frac{T_{\mathrm{b}}-T_{\mathrm{u}}}{T_{\mathrm{u}}}.$

De plus, si nous définissons une température adimensionnelle

$\theta =\frac{T-T_{\mathrm{u}}}{T_{\mathrm{b}}-T_{\mathrm{u}}},0\le \theta \le 1$

telle que $\theta$ s'approche de zéro dans la région des gaz imbrûlés et s'approche de l'unité dans la région des gaz brûlés, alors le rapport entre la vitesse de réaction à n'importe quelle température et la vitesse de réaction à la température des gaz brûlés est donné par :

$\frac{\omega (T)}{\omega (T_{\mathrm{b}})}\propto \frac{{\mathrm{e}}^{-E_{\mathrm{a}}/RT}}{{\mathrm{e}}^{-E_{\mathrm{a}}/R{T}_{\mathrm{b}}}}=\exp [-\beta (1-\theta )\frac{1+q}{1+q\theta }].$

Pour $\beta \to \mathrm{\infty }$ (grande énergie d'activation) avec $q\sim O(1)$, la vitesse de réaction est exponentiellement petite, en $O(e^{-\beta })$ et négligeable partout, sauf lorsque $\beta (1-\theta )\sim O(1)$. En d'autres termes, la vitesse de réaction est négligeable presque partout, sauf dans une petite région très proche de la température des gaz brûlés où $1-\theta \sim O(1/\beta )$. On identifie par ces critères deux régions différentes :

• la zone convective-diffusive externe,
• la couche réactive-diffusive interne.

où, dans la zone convective-diffusive, le terme de réaction sera négligé et dans la fine couche réactive-diffusive, les termes convectifs peuvent être négligés. Les solutions dans ces deux régions sont raccordées par les pentes en utilisant la méthode des développements asymptotiques raccordés.

On peut mettre en évidence une structure plus complexe à partir des longueurs caractéristiques présentes : convection, diffusion et réaction.