Formulation de Schvab-Zeldovitch

La formulation de Schvab-Zeldovitch (en anglais Shvab–Zeldovich formulation) a été développée par V. A. Schvab (1948) et Iakov Zeldovitch (1949) pour décrire simplement la dynamique d'un mélange gazeux réactif sous l'hypothèse d'un nombre de Prandtl et d'un nombre de Lewis unités. Elle montre d'une manière générale l'analogie formelle entre l'enthalpie totale et les fractions élémentaires dans le milieu. Elle est employée dans le cas plus spécifique de la combustion où Amable Liñán l'a généralisée pour un nombre de Lewis différent de l'unité.

Équations de conservation pour un milieu hors équilibre chimique

Conservation des espèces et de l'énergie

Les équation de conservation de la masse et de l'énergie pour un milieu réactif comportant $N_e$ espèces s'écriventModèle:,Modèle:,Modèle:, :

$\frac{\partial (\rho c_i)}{\partial t}+\overrightarrow{\mathrm{\nabla }}\cdot [\rho c_i\overrightarrow{V}+\rho {\overrightarrow{V}}_{D_i}]={\omega }_i,i=1,N_e$

$\frac{\partial (\rho E)}{\partial t}+\overrightarrow{\mathrm{\nabla }}\cdot (\rho E\overrightarrow{V})=-p\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{V}-\overrightarrow{\mathrm{\nabla }}\cdot \overrightarrow{q}+\frac{1}{2\mu }||\mathrm{\Sigma }|{|}_F^2$

où

• $\rho$ est la masse volumique ;
• $\overrightarrow{V}$ la vitesse ;
• $E=e+\frac{||\overrightarrow{V}|{|}^2}{2}$ est l'énergie totale par unité de masse et $e=\rho \sum _{i=1}^{N_e}{c}_i{e}_i$ l'énergie interne par unité de masse ;
• $\overrightarrow{q}=-\lambda \overrightarrow{\mathrm{\nabla }}T+\rho \sum _{i=1}^{N_e}{c}_i{h}_i{\overrightarrow{V}}_{D_i}$ est le flux de chaleur par unité de surface où on néglige l'effet Dufour ;
• $\lambda$ est la conductivité thermique ;
• $h_i=e_i+\frac{p_i}{\rho }$ l'enthalpie par unité de masse de l'espèce Modèle:Mvar ;
• $c_i$ est la fraction massique :
• $p=\sum _{i=1}^{N_e}{p}_i$ est la pression et $p_i$ la pression partielle de l'espèce Modèle:Mvar ;
• ${\overrightarrow{V}}_{D_i}=-{𝒟}_i\overrightarrow{\mathrm{\nabla }}{c}_i$ est une approximation de la vitesse de diffusion par gradient de concentration, l'effet des gradients de température (effet Soret) et de pression sont négligés ;
• ${\omega }_i$ est la masse formée par unité de volume et par unité de temps ;
• $||\mathrm{\Sigma }|{|}_F$ est la norme de Frobenius du tenseur des contraintes visqueuses ;
• $\mu$ est la viscosité dynamique.

Le rang du système de conservation des espèces est $N_e-1$ puisque par définition $\sum _{i=1}^{N_e}{c}_i=1$, $\sum _{i=1}^{N_e}\rho {\overrightarrow{V}}_{D_i}=0$ et qu'il ny a pas de création de masse $\sum _{i=1}^{N_e}{\omega }_i=0$. La somme de ces équations donne la conservation de masse globale :

$\frac{\partial \rho }{\partial t}+\overrightarrow{\mathrm{\nabla }}\cdot (\rho \overrightarrow{V})=0$

Conservation des éléments

On peut également écrire une équation de conservation pour chacun des $N_A$ éléments présents dans le milieu, dont les fractions massiques sont :

${\tilde{c}}_k={ℳ}_k\sum _{i=1}^N{n}_{ki}\frac{c_i}{{ℳ}_i},k=1,N_A$

où $n_{ki}$ est le nombre d'éléments Modèle:Mvar dans l'espèce Modèle:Mvar et ${ℳ}_k$ la masse molaire de l'atome Modèle:Mvar.

Dans l'hypothèse où tous les coefficients de diffusion sont égaux $\mathrm{\forall }i{𝒟}_i=𝒟$ (c'est le cas si l'on suppose tous les coefficients binaire égaux, hypothèse assez mal vérifiée) et compte tenu de ${\omega }_k=0$ on obtient un système de conservation pour les fractions élémentaires en sommant les équations sur les espèces prémultipliées par ${ℳ}_k\sum _{i=1}^{N_e}\frac{n_{ki}}{{ℳ}_i}$ :

$\frac{\partial (\rho {\tilde{c}}_k)}{\partial t}+\overrightarrow{\mathrm{\nabla }}\cdot (\rho {\tilde{c}}_k\overrightarrow{V})=\overrightarrow{\mathrm{\nabla }}\cdot (\rho 𝒟\overrightarrow{\mathrm{\nabla }}{\tilde{c}}_k),k=1,N_A$

ou, en tenant compte de la conservation de masse globale :

$\rho \frac{\partial {\tilde{c}}_k}{\partial t}+\rho \overrightarrow{V}\cdot \overrightarrow{\mathrm{\nabla }}{\tilde{c}}_k=\overrightarrow{\mathrm{\nabla }}\cdot (\rho 𝒟\overrightarrow{\mathrm{\nabla }}{\tilde{c}}_k),k=1,N_A$

Méthode de Shvab-Zeldovich

L'équation de conservation de l'énergie peut être réécrite sous la forme suivanteModèle:, :

$\frac{\partial h_T}{\partial t}+\rho \overrightarrow{V}\cdot \overrightarrow{\mathrm{\nabla }}{h}_T=\overrightarrow{\mathrm{\nabla }}\cdot [\frac{\mu }{𝒫𝓇}\overrightarrow{\mathrm{\nabla }}{h}_T+\mu (1-\frac{1}{𝒫𝓇})\overrightarrow{\mathrm{\nabla }}\left(\frac{V^2}{2}\right)+(1-\frac{1}{ℒℯ})\sum _{i=1}^{N_e}\rho 𝒟h_i\overrightarrow{\mathrm{\nabla }}{c}_i]$

où

• $𝒫𝓇=\frac{\mu C_p}{\lambda }$ est le nombre de Prandtl ;
• $ℒℯ=\frac{\lambda }{\rho C_p𝒟}$ est le nombre de Lewis ;
• $h_T=\sum _{i=1}^{N_e}{c}_i{h}_i+\frac{V^2}{2}$ est l'enthalpie totale.

Ces relations sont liées par :

$h_i=h_i^0+\underset{T_0}{\overset{T}{\int }}{C}_{p}d{T}^{\prime },C_p=\sum _{i=1}^{N_e}{c}_i{C}_{p_i}$

On voit que l'hypothèse $ℒℯ=1$ fait disparaître la diffusion de l'équation. Si de plus on suppose $𝒫𝓇=1$ on obtient $\mu =\rho 𝒟$ et l'équation de conservation de l'énergie se simplifie considérablement :

$\frac{\partial h_T}{\partial t}+\rho \overrightarrow{V}\cdot \overrightarrow{\mathrm{\nabla }}{h}_T=\overrightarrow{\mathrm{\nabla }}\cdot (\rho 𝒟\overrightarrow{\mathrm{\nabla }}{h}_T)$

Si on rapproche cette équation de celle de conservation des fractions élémentaires on voit qu'elles sont identiques : leurs solutions sont donc linéairement dépendantes. Toutefois la connaissance des fractions élémentaires ne permet pas en général de remonter à celles des fraction des espèces. Toutefois, dans le cas où il n'existe qu'une seule réaction, il est possible de définir un terme de production indépendant de l'espèce. Cette approximation est souvent utilisée dans les problèmes de combustion dans lesquels la diffusion est un phénomène prépondérant et les effets visqueux négligés.

Méthode de Shvab-Zeldovich-Liñán

Les nombres de Lewis sont plus ou moins proches de l'unité : Modèle:Unité pour $CO$, 0,99 pour $C{H}_4$, 1,1 pour l'air mais 0,41 pour $H_2$ (valeurs à Modèle:Unité). La méthode de Shvab-Zeldovich a été modifiée par Amable Liñán en 1991 pour corriger des effets du nombre de LewisModèle:,.

Liñán considère le cas de la combustion de combustible et d'oxydant, décrite par une flamme infiniment mince (limite de Burke–Schumann correspondante à une chimie infiniment rapide), dans l'approximation de faible nombre de Mach. La réaction est à une seule étape, les fractions massiques de l'oxydant et du combustible sont $c_O$ et $c_C$, respectivement. La température est adimensionnée par $\frac{Q}{C_p}$ où $C_p$ est la chaleur massique et $Q$ l'énergie par unité de masse issue de la combustion complète d'un mélange stœchiométrique. Ce milieu supposé ouvert (pression donnée) est décrit par les équations suivantes en négligeant la gravité :

$\begin{array}{cc}\rho \frac{\partial c_C}{\partial t}+\rho \overrightarrow{V}\cdot \overrightarrow{\mathrm{\nabla }}{c}_C& =\frac{1}{{ℒℯ}_C}\overrightarrow{\mathrm{\nabla }}\cdot (\rho {𝒟}_T\mathrm{\nabla }{c}_C)-\omega \\ \rho \frac{\partial c_O}{\partial t}+\rho \overrightarrow{V}\cdot \overrightarrow{\mathrm{\nabla }}{c}_O& =\frac{1}{{ℒℯ}_O}\overrightarrow{\mathrm{\nabla }}\cdot (\rho {𝒟}_T\mathrm{\nabla }{c}_O)-S\omega \\ \rho \frac{\partial T}{\partial t}+\rho \overrightarrow{V}\cdot \overrightarrow{\mathrm{\nabla }}T& =\mathrm{\nabla }\cdot (\rho {𝒟}_T\overrightarrow{\mathrm{\nabla }}T)+\omega \end{array}$

${ℒℯ}_C$ et ${ℒℯ}_O$ sont les nombres de Lewis. $S$ est le rapport de mélange. Le taux de production chimique adimensionné est :

$\omega =𝒟𝒶c_C{c}_O{e}^{-\frac{E_a}{ℛT}}$

où $𝒟𝒶$ est le nombre de Damköhler et $E_a$ l'énergie d'activation.

Dans l'hypothèse de la flamme mince combustible et oxydant ne coexistent pas $c_C{c}_O=0$ presque partout. $\omega$ est une distribution de Dirac.

Pour résoudre ce problème on introduit les fonctions suivantes :

$\begin{array}{cc}Z=\frac{S{c}_C-c_O+1}{S+1},& \tilde{Z}=\frac{\tilde{S}{c}_C-c_O+1}{\tilde{S}+1},\\ H=\frac{T-T_0}{T_s-T_0}+c_C+c_O-1,& \tilde{H}=\frac{T-T_0}{T_s-T_0}+\frac{c_O}{L{e}_O}+\frac{c_C-1}{L{e}_C}\end{array}$

où $\tilde{S}=SL{e}_O/L{e}_C$, $T_0$ est la température du flux de carburant et $T_s$ est la température de flamme adiabatique, toutes deux réduites par la température du flux d'oxydant. L'introduction de ces fonctions réduit les équations de conservation à :

$\begin{array}{cc}\rho \frac{\partial Z}{\partial t}+\rho 𝐯\cdot \mathrm{\nabla }Z& =\frac{1}{L{e}_m}\mathrm{\nabla }\cdot (\rho {𝒟}_T\mathrm{\nabla }\tilde{Z})\\ \rho \frac{\partial H}{\partial t}+\rho 𝐯\cdot \mathrm{\nabla }H& =\mathrm{\nabla }\cdot (\rho {𝒟}_T\mathrm{\nabla }\tilde{H})\end{array}$

où $L{e}_m=L{e}_O(S+1)/(\tilde{S}+1)$ est le nombre de Lewis moyen (ou effectif). La relation entre $Z$ et $\tilde{Z}$ et entre $H$ et $\tilde{H}$ peut être déduite des quantités décrivant la flamme :

• à la surface stœchiométrique (la surface de la flamme), $c_C$ et $c_O$ sont nuls, ce qui conduit à $Z=Z_s=1/(S+1)$, $\tilde{Z}={\tilde{Z}}_s=1/(\tilde{S}+1)$, $H=H_s=(T_f-T_0)/(T_s-T_0)-1$ et $\tilde{H}={\tilde{H}}_s=(T_f-T_0)/(T_s-T_0)-1/L{e}_F$, où $T_f$ est la température de la flamme (réduite par la température du flux d'oxydant) qui n'est, en général, pas égale à $T_s$ sauf si $L{e}_F=L{e}_O=1$.
• Côté combustible, puisque $c_C-1=c_O=T-T_0=0$, on a $Z-1=\tilde{Z}-1=H=\tilde{H}=0$.
• De même, côté oxydant, puisque $Y_F=Y_O-1=T-1=0$, on a $Z=\tilde{Z}=H-(1-T_0)/(T_s-T_0)=\tilde{H}-(1-T_0)/(T_s-T_0)-1/L{e}_O+1/L{e}_C=0$.

La condition d'équilibre définit $\tilde{Z}$ :

$\begin{array}{cc}\tilde{Z}<{\tilde{Z}}_s& c_C=0,c_O=1-\frac{\tilde{Z}}{{\tilde{Z}}_s}=1-\frac{Z}{Z_s},\\ \tilde{Z}>{\tilde{Z}}_s& c_O=0,c_C=\frac{\tilde{Z}-{\tilde{Z}}_s}{1-{\tilde{Z}}_s}=\frac{Z-Z_s}{1-Z_s}\end{array}$

Les relations ci-dessus définissent la fonction par morceaux $Z(\tilde{Z})$ :

$Z=\{\begin{array}{ccc}\tilde{Z}/L{e}_m& \text{si}& \tilde{Z}<{\tilde{Z}}_s\\ Z_s+Le(\tilde{Z}-{\tilde{Z}}_s)/L{e}_m& \text{si}& \tilde{Z}>{\tilde{Z}}_s\end{array}$

où $L{e}_m={\tilde{Z}}_s/Z_s=(S+1)/(S/L{e}_F+1)$ est un nombre de Lewis moyen.

Cela conduit à une équation non linéaire pour $\tilde{Z}$. Étant donné que $H-\tilde{H}$ n'est fonction que de $Y_F$ et $Y_O$, les expressions ci-dessus peuvent être utilisées pour définir la fonction $H(\tilde{Z},\tilde{H})$ :

$H=\tilde{H}+\{\begin{array}{ccc}(1/L{e}_F-1)-(1/L{e}_O-1)(1-\tilde{Z}/{\tilde{Z}}_s)& \text{si}& \tilde{Z}<{\tilde{Z}}_s\\ (1/L{e}_F-1)(1-\tilde{Z})/(1-{\tilde{Z}}_s)& \text{si}& \tilde{Z}>{\tilde{Z}}_s\end{array}$

Le problème peut être résolu avec des conditions aux limites appropriées pour $\tilde{H}$.

On montre que $\tilde{Z}$ et $\tilde{H}$ sont des fonctions continues à dérivées continues à la traversée de la flamme alors que $Z$ et $H$ présentent des sauts du gradient.