Théorie de Frank-Kamenetskii

Avant David Frank-Kamenetskii, son directeur de thèse Nikolaï Semionov a proposé une théorie de l'explosion thermique avec un modèle plus simple dans lequel il suppose une fonction linéaire pour le processus de conduction thermique en lieu et place de l'opérateur Laplacien. L'équation de Semionov en géométrie plane s'écrit :

${\displaystyle {\frac {d\theta }{d\tau }}=e^{\theta }-{\frac {\theta }{\delta }},\quad \theta (0)=0}$

dans laquelle le terme exponentiel ${\displaystyle e^{\theta }}$ tend à augmenter la température réduite ${\displaystyle \theta }$ tandis que le terme linéaire ${\displaystyle -\theta /\delta }$ tend à la diminuer. L'importance relative entre les deux termes est déterminée par le nombre de Damköhler ${\displaystyle \delta }$. La solution numérique de l'équation ci-dessus pour différentes valeurs de ${\displaystyle \delta }$ est illustrée dans la figure.

Lorsque ${\displaystyle 0<\delta <e^{-1}}$, le terme linéaire finit par dominer et le système atteint un état stationnaire lorsque ${\displaystyle \tau \rightarrow \infty }$. À l'état stationnaire (${\displaystyle d\theta /d\tau =0}$), l'équilibre est donné par l'équation :

${\displaystyle \delta e^{\theta }=\theta ,\quad \Rightarrow \quad \theta =-W(-\delta )}$

où ${\displaystyle W}$ représente la fonction W de Lambert. À partir des propriétés de cette fonction il est facile de voir que la température à l'état stationnaire fournie par l'équation ci-dessus n'existe que lorsque ${\displaystyle \delta \leq \delta _{c}=1/e}$, où ${\displaystyle \delta _{c}}$ est appelé paramètre de Frank-Kamenetskii en tant que point critique où le système bifurque de l'existence d'un état stationnaire à un état explosif pour les grandes valeurs du temps.

Pour ${\displaystyle \delta _{c}<\delta <\infty }$, le système explose puisque le terme exponentiel domine lorsque la durée augmente. En raison du forçage exponentiel ${\displaystyle \theta \rightarrow \infty }$ pour une valeur finie de ${\displaystyle \tau }$. Ce temps est interprété comme le temps d'allumage ou le temps d'induction du système. Lorsque ${\displaystyle \delta \gg \delta _{c}}$ le terme de conduction thermique ${\displaystyle -\theta /\delta }$ peut être négligé, auquel cas le problème admet une solution explicite :

${\displaystyle {\frac {d\theta }{d\tau }}=e^{\theta }\quad \Rightarrow \quad \theta =\ln \left({\frac {1}{1-\tau }}\right)}$

À l'instant ${\displaystyle \tau =1}$, le système explose. Cet instant est également appelé « période d'induction adiabatique » puisque le terme de conduction thermique ${\displaystyle -\theta /\delta }$ est négligé.

Dans la condition quasi-critique, c'est-à-dire lorsque ${\displaystyle \delta -\delta _{c}\rightarrow 0^{+}}$, le système met très longtemps à exploser. L'analyse de cette limite a été réalisée pour la première fois par Frank-Kamenetskii^[12] bien que des asymptotiques appropriées n'aient été montrées que plus tard par D. R. Kassoy et Amable Liñán^[13] y compris la consommation de réactifs car celle-ci n'est pas négligeable lorsque ${\displaystyle \tau \sim \beta \gamma }$.

Une analyse simplifiée ignorant la consommation de réactifs est présentée ici. Définissons un petit paramètre ${\displaystyle \epsilon \rightarrow 0^{+}}$ tel que ${\displaystyle \delta =\delta _{c}(1+\epsilon )}$. Dans ce cas, l'évolution temporelle de ${\displaystyle \theta }$ est la suivante : elle croît d'abord jusqu'à la valeur de température à l'état stationnaire correspondant à ${\displaystyle \delta =\delta _{c}}$, qui est donnée par ${\displaystyle \theta =-W(-\delta _{c})=1}$ à des instants d'ordre ${\displaystyle \tau \sim O(1)}$, puis il reste très proche de cette valeur à l'état stationnaire pendant un temps long avant de finalement exploser à un temps long. La quantité qui nous intéresse est l'estimation à long terme de l'explosion. Pour trouver l'estimation on introduit les transformations ${\displaystyle \zeta ={\sqrt {\epsilon }}\tau }$ et ${\displaystyle \theta =1+{\sqrt {\epsilon }}\psi (\zeta )-\epsilon }$ utilisables dans la région où ${\displaystyle \theta }$ reste proche de ${\displaystyle 1}$ dans l'équation de référence. En ne retenant que les termes d'ordre un on obtient :

${\displaystyle \delta _{c}{\frac {d\psi }{d\zeta }}=1+{\frac {\psi ^{2}}{2}},\quad \psi (0)\rightarrow -(1-\theta )/{\sqrt {\epsilon }}=-\infty }$

où la condition limite est dérivée en faisant correspondre avec la région initiale dans laquelle ${\displaystyle \tau ,\theta \sim 1}$. La solution de cette équation est donnée par :

${\displaystyle \psi =-{\sqrt {2}}\cot \left({\frac {\zeta }{{\sqrt {2}}\delta _{c}}}\right)}$

qui montre que ${\displaystyle \psi \rightarrow \infty }$ lorsque ${\displaystyle \zeta ={\sqrt {2}}\pi \delta _{c}}$. En écrivant cette condition en termes de ${\displaystyle \tau }$, le temps d'explosion dans la condition quasi-critique est :

${\displaystyle \tau ={\sqrt {\frac {2\pi ^{2}\delta _{c}^{3}}{\delta -\delta _{c}}}}}$

L'instant d'explosion ${\displaystyle \tau \rightarrow \infty }$ quand ${\displaystyle \delta -\delta _{c}\rightarrow 0^{+}}$.

Le seul paramètre caractérisant l'explosion est le nombre de Damköhler ${\displaystyle \delta }$. Lorsque ${\displaystyle \delta }$ est très élevé, le temps de conduction est plus long que le temps de réaction chimique et le système explose quand il atteint une haute température. D'autre part, lorsque ${\displaystyle \delta }$ est très faible, le temps de conduction thermique est beaucoup plus rapide que le temps de réaction chimique, de sorte que toute la chaleur produite par la réaction chimique est immédiatement conduite vers la paroi, il n'y a donc pas d'explosion, on passe à un état presque stationnaire. Amable Liñán a qualifié ce phénomène de mode de réaction lente. À un nombre de Damköhler critique ${\displaystyle \delta _{c}}$, le système passe du mode de réaction lente au mode explosif. Au lieu de résoudre le problème complet pour trouver ce ${\displaystyle \delta _{c}}$, Frank-Kamenetskii a résolu le problème de l'état stationnaire pour différents nombres de Damköhler jusqu'à la valeur critique au-delà de laquelle aucune solution stationnaire n'existe. Le problème à résoudre est décrit par l'équation de Liouville – Bratu – Gelfand^[14],[15] :

${\displaystyle \nabla ^{2}\theta +\delta e^{\theta }=0}$

avec des conditions initiale ${\displaystyle \theta (1)=0}$ et la condition de symétrie ${\displaystyle {\frac {d\theta }{d\eta }}_{\eta =0}=0}$.

Dans le cas plan, l'équation s'écrit :

${\displaystyle {\frac {d^{2}\theta }{d\eta ^{2}}}=-\delta e^{\theta }}$

On introduit les transformations ${\displaystyle \Theta =\theta _{m}-\theta }$ et ${\displaystyle \xi ^{2}=\delta e^{\theta _{m}}\eta ^{2}}$, où ${\displaystyle \theta _{m}}$ est la température maximale qui se produit à ${\displaystyle \eta =0}$ en raison de la symétrie, alors :

${\displaystyle {\frac {d^{2}\Theta }{d\xi ^{2}}}=e^{-\Theta },\quad \Theta (0)=0,\quad {\frac {d\Theta }{d\xi }}_{\xi =0}=0}$

En intégrant une fois et en utilisant la condition de symétrie l'équation devient :

${\displaystyle {\frac {d\Theta }{d\xi }}={\sqrt {2(1-e^{-\Theta })}}}$

En intégrant à nouveau :

${\displaystyle e^{(\theta _{m}-\theta )/2}=\cosh \left(\eta {\sqrt {\frac {\delta e^{\theta _{m}}}{2}}}\right)}$

L'équation ci-dessus est la solution exacte, mais la température maximale de ${\displaystyle \theta _{m}}$ est inconnue. En utilisant la condition limite au bord ${\displaystyle \theta (1)=0}$ la température maximale est obtenue à partir de l'expression implicite :

${\displaystyle e^{\theta _{m}/2}=\cosh {\sqrt {\frac {\delta e^{\theta _{m}}}{2}}}\quad {\text{ou}}\quad \delta =2e^{-\theta _{m}}\left(\operatorname {arcosh} e^{\theta _{m}/2}\right)^{\!2}}$

La valeur critique ${\displaystyle \delta _{c}}$ est obtenue en trouvant le maximum de l'équation ${\displaystyle d\delta /d\theta _{m}=0}$ à ${\displaystyle \delta _{c}}$ (voir figure) :

${\displaystyle {\frac {d\delta }{d\theta _{m}}}=0\quad \Rightarrow \quad e^{\theta _{m,c}/2}-{\sqrt {e^{\theta _{m,c}}-1}}\operatorname {arcosh} e^{\theta _{m,c}/2}=0}$

${\displaystyle \theta _{m,c}=1,1868\quad \Rightarrow \quad \delta _{c}=0,8785}$

Le système n'a pas d'état stationnaire (ou explose) pour ${\displaystyle \delta >\delta _{c}=0.8785}$.

Pour ${\displaystyle \delta <\delta _{c}=0.8785}$ le système est dans un état quasi-stationnaire avec une réaction très lente.

Paul L. Chambré a fourni une solution exacte du problème en 1952^[16]. H. Lemke a également fourni une solution sous une forme quelque peu différente en 1913^[17].

L'équation s'écrit ici :

${\displaystyle {\frac {1}{\eta }}{\frac {d}{d\eta }}\left(\eta {\frac {d\theta }{d\eta }}\right)=-\delta e^{\theta }}$

Si les transformations ${\displaystyle \omega =\eta d\theta /d\eta }$ et ${\displaystyle \chi =e^{\theta }\eta ^{2}}$ sont introduites il vient :

${\displaystyle {\frac {d\omega }{d\chi }}=-{\frac {\delta }{2+\omega }},\quad \omega (0)=0}$

La solution générale est ${\displaystyle \omega ^{2}+4\omega +C=-2\delta \chi }$. ${\displaystyle C=0}$ d'après la condition de symétrie. En réécrivant dans la variable d'origine, l'équation est :

${\displaystyle \eta ^{2}\left({\frac {d\theta }{d\eta }}\right)^{2}+4\eta {\frac {d\theta }{d\eta }}=-2\delta \eta ^{2}e^{\theta }}$

Mais l'équation d'origine multipliée par ${\displaystyle 2\eta ^{2}}$ est :

${\displaystyle 2\eta ^{2}{\frac {d^{2}\theta }{d\eta ^{2}}}+2\eta {\frac {d\theta }{d\eta }}=-2\delta \eta ^{2}e^{\theta }}$

En soustrayant les deux dernières équations l'une de l'autre, on obtient :

${\displaystyle {\frac {d^{2}\theta }{d\eta ^{2}}}-{\frac {1}{\eta }}{\frac {d\theta }{d\eta }}-{\frac {1}{2}}\left({\frac {d\theta }{d\eta }}\right)^{2}=0}$

On introduit la fonction ${\displaystyle g(\eta )=d\theta /d\eta }$ qui vérifie donc :

${\displaystyle {\frac {dg}{d\eta }}-{\frac {g}{\eta }}-{\frac {g^{2}}{2}}=0}$

Il s'agit d'une équation différentielle de Bernoulli d'ordre ${\displaystyle 2}$ dont la solution est :

${\displaystyle g(\eta )={\frac {d\theta }{d\eta }}=-{\frac {4\eta }{B'+\eta ^{2}}}}$

En intégrant celle-ci on obtient ${\displaystyle \theta =A-2\ln(B\eta ^{2}+1)}$ où ${\displaystyle B=1/B'}$. Il reste une condition pour les deux constantes ${\displaystyle A,\ B}$. Mais celles-ci sont liées comme l'on peut le constater en remplaçant la solution ci-dessus dans l'équation de départ. On trouve ainsi ${\displaystyle A=\ln(8B/\delta )}$. Par conséquent, la solution est :

${\displaystyle \theta =\ln {\frac {8B/\delta }{(B\eta ^{2}+1)^{2}}}}$

Si nous utilisons maintenant l'autre condition limite au bord du domaine ${\displaystyle \theta (1)=0}$, nous obtenons une équation pour ${\displaystyle B}$ telle que ${\displaystyle \delta (B+1)^{2}-8B=0}$. La valeur maximale de ${\displaystyle \delta }$ pour laquelle une solution est possible est lorsque ${\displaystyle B=1}$, donc le paramètre critique de Frank-Kamentskii est ${\displaystyle \delta _{c}=2}$. Le système n'a pas d'état stationnaire (ou explose) pour ${\displaystyle \delta >\delta _{c}=2}$ et pour ${\displaystyle \delta <\delta _{c}=2}$, le système passe à un état stationnaire avec une réaction très lente. La température maximale ${\displaystyle \theta _{m}=\theta (\eta =0)}$ vaut :

${\displaystyle \theta _{m}=\ln {\frac {8B}{\delta }}\quad {\text{ou}}\quad \delta =8Be^{-\theta _{m}}}$

Pour chaque valeur de ${\displaystyle \delta }$, nous avons deux valeurs de ${\displaystyle \theta _{m}}$ puisque ${\displaystyle B}$ est à valeurs multiples. La température critique maximale est ${\displaystyle \theta _{m,c}=\ln 4}$.

Pour la géométrie sphérique, il n'existe pas de solution explicite connue pour l'équation :

${\displaystyle {\frac {1}{\eta ^{2}}}{\frac {d}{d\eta }}\left(\eta ^{2}{\frac {d\theta }{d\eta }}\right)=-\delta e^{\theta }}$

Si les transformations ${\displaystyle \Theta =\theta _{m}-\theta }$ et ${\displaystyle \xi ^{2}=\delta e^{\theta _{m}}\eta ^{2}}$, où ${\displaystyle \theta _{m}}$ est la température maximale qui se produit à ${\displaystyle \eta =0}$ en raison de la symétrie, sont introduites il vient :

${\displaystyle {\frac {1}{\xi ^{2}}}{\frac {d}{d\xi }}\left(\xi ^{2}{\frac {d\Theta }{d\xi }}\right)=e^{-\Theta },\quad \Theta (0)=0,\quad {\frac {d\Theta }{d\xi }}_{\xi =0}=0}$

Cette équation est l'équation de Emden-Chandrasekhar^[18]. Contrairement au cas plan et cylindrique, le cas sphérique a une infinité de solutions pour ${\displaystyle \delta <\delta _{c}}$, oscillant autour du point ${\displaystyle \delta =2}$^[19] au lieu de seulement deux solutions, ce qui a été montré par Israel Gelfand^[20]. La branche basse correspond au comportement explosif.

À partir de la résolution numérique, on constate que le paramètre critique de Frank-Kamenetskii est ${\displaystyle \delta _{c}=3.3220}$. Le système n'a pas d'état stationnaire (ou explose) pour ${\displaystyle \delta >\delta _{c}=3.3220}$. Pour ${\displaystyle \delta <\delta _{c}=3.3220}$ le système passe à un état stationnaire avec une réaction très lente. La température maximale ${\displaystyle \theta _{m}}$ se produit à ${\displaystyle \eta =0}$ et la température critique maximale est ${\displaystyle \theta _{m,c}=1,6079}$.

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Frank-Kamenetskii theory » (voir la liste des auteurs).

• Portail de la physique