Flamme de Burke-Schumann

En combustion la flamme de Burke-Schumann est une de flamme de diffusion à l'embouchure des deux cylindres concentriques injectant séparément l'oxydant et le combustible. Son nom vient de S. P. Burke et T. E. W. SchumannModèle:, qui ont prédit la hauteur et la forme de la flamme en utilisant l'hypothèse de chimie infiniment rapide (limite de Burke–Schumann) en 1928 lors du premier symposium sur la combustion.

Description mathématique

Considérons un conduit cylindrique d'axe $z$ et de rayon $a$ à travers lequel le combustible est alimenté par le bas (voir figure). Son embouchure est située à $z=0$. L'oxydant est alimenté selon le même axe dans le tube concentrique de rayon $b>a$. Soit $Y_{Fo}$ la fraction massique dans le tube de carburant (qui peut être dilué dans un gaz neutre) et $Y_{Oo}$ celle de l'oxygène dans le conduit extérieur (qui peut être alimenté en air). Le mélange du carburant et de l'oxygène se produit dans la région $z>0$. Les hypothèses suivantes sont faites dans l'analyse :

• La vitesse moyenne est parallèle à l'axe $z$ .
• Le flux de masse $\rho v$ dans la direction axiale est constant.
• La diffusion axiale est négligeable par rapport à la diffusion transversale (radiale).
• La réaction est infiniment rapide (limite de Burke–Schumann), par conséquent la flamme apparaît comme une nappe de réaction constituant une discontinuité des diverses quantités.
• Les effets de la gravité son négligés.

La réaction chimique est irréversible et en une seule étape, décrite par une loi d'Arrhenius :$\mathrm{C}\mathrm{o}\mathrm{m}\mathrm{b}\mathrm{u}\mathrm{s}\mathrm{t}\mathrm{i}\mathrm{b}\mathrm{l}\mathrm{e}+s{\mathrm{O}}_2\to (1+s)\mathrm{P}\mathrm{r}\mathrm{o}\mathrm{d}\mathrm{u}\mathrm{i}\mathrm{t}\mathrm{s}+\mathrm{C}\mathrm{h}\mathrm{a}\mathrm{l}\mathrm{e}\mathrm{u}\mathrm{r}$

où $s$ est la masse d'oxygène nécessaire pour brûler une unité de masse de carburant.

$\omega$ est la masse de carburant brûlé par unité de volume par unité de temps et on introduit les fractions massiques normalisées de carburant et d'oxydant ainsi que le paramètre de stœchiométrie :

$y_F=\frac{Y_F}{Y_{Fo}},y_O=\frac{Y_O}{Y_{Oo}},S=\frac{s{Y}_{Fo}}{Y_{Oo}}$

Les équations régissant la fraction massique de carburant et d'oxydant s'écrivent :

$\begin{array}{c}\frac{\rho D_T}{r}\frac{\partial }{\partial r}\left(r\frac{\partial y_F}{\partial r}\right)-\rho v\frac{\partial y_F}{\partial z}=\frac{\omega }{Y_{Fo}}\\ \frac{\rho D_T}{r}\frac{\partial }{\partial r}\left(r\frac{\partial y_O}{\partial r}\right)-\rho v\frac{\partial y_O}{\partial z}=S\frac{\omega }{Y_{Fo}}\\ \end{array}$

où le nombre de Lewis des deux espèces est supposé égal à un et $\rho D_T$ est supposé constant, où $D_T$ est la diffusivité thermique. Les conditions aux limites du problème sont

$\begin{array}{cccc}\text{en}z=0:& 0<r<a,& y_F=1,& y_O=0\\ \text{en}z=0:& a<r<b,& y_F=0,& y_O=1\\ \text{en}r=b:& 0<z<\mathrm{\infty },& \frac{\partial y_F}{\partial r}=0,& \frac{\partial y_O}{\partial r}=0\end{array}$

L'équation peut être combinée linéairement pour éliminer le terme de réaction non linéaire $\omega /Y_{Fo}$ et introduire la fraction de mélange :

$Z=\frac{S{y}_F-y_O+1}{S+1}$

Cette quantité prend la valeur unité dans le flux de combustible et de zéro dans le flux d'oxydant et c'est un champ scalaire qui n'est pas affecté par la réaction. L'équation satisfaite par $Z$ est :

$\frac{1}{r}\frac{\partial }{\partial r}\left(r\frac{\partial Z}{\partial r}\right)-\frac{\rho v}{\rho D_T}\frac{\partial Z}{\partial z}=0$

Si les nombres de Lewis du combustible et du comburant ne sont pas égaux à un, alors l'équation satisfaite par $Z$ est non linéaire, donnée par la formulation de Shvab–Zeldovich–Liñán.

On introduit les coordonnées suivante :

$\xi =\frac{r}{b},\eta =\frac{\rho D_T}{\rho v}\frac{z}{b^2},c=\frac{a}{b}$

Avec celles-ci le système devient :

$\frac{1}{\xi }\frac{\partial }{\partial \xi }\left(\xi \frac{\partial Z}{\partial \xi }\right)-\frac{\partial Z}{\partial \eta }=0$

Les conditions aux limites correspondantes deviennent :

$\begin{array}{ccc}\text{en}\eta =0:& 0<\xi <c& Z=1\\ \text{en}\eta =0:& c<\xi <1& Z=0\\ \text{en}\xi =1:& 0<\eta <\mathrm{\infty }& \frac{\partial Z}{\partial \xi }=0\end{array}$

L'équation peut être résolue par séparation des variables :

$Z(\xi ,\eta )=c^2+2c\sum _{n=1}^{\mathrm{\infty }}\frac{1}{{\lambda }_n}\frac{J_1(c{\lambda }_n)}{J_0^2({\lambda }_n)}{J}_0({\lambda }_n\xi )e^{-{\lambda }_n^2\eta }$

où $J_0$ et $J_1$ sont la fonction de Bessel de première espèce et ${\lambda }_n$ est la racine n-ième de $J_1(\lambda )=0$.

Le problème en géométrie plane peut être résolu de la même façonModèle:,.

Forme et hauteur de la flamme

Position de la flamme pour une faible alimentation en air (1) ou une faible alimentation en combustible (2) et diagrammes $Y(Z)$ et $T(Z)$

Dans la limite de Burke–Schumann la flamme est considérée comme une fine couche de réaction en dehors de laquelle le combustible et l'oxygène ne peuvent pas coexister : $y_F{y}_O=0$ presque partout. La couche de réaction elle-même est constituée par la surface stœchiométrique $S{y}_F=y_O$, soit :

$Z=Z_s\equiv \frac{1}{S+1}$

où $Z_s$ est la fraction de mélange stœchiométrique.

La couche de réaction sépare la région du combustible et de l'oxydant. La structure interne de la couche réactionnelle est décrite par l'équation de Liñán. Du côté combustible de la nappe de réaction ($Z>Z_s$) on a :

$y_F=\frac{Z-Z_s}{1-Z_s},y_O=0$

et du côté oxydant $Z<Z_s$ :

$y_F=0,y_O=1-\frac{Z}{Z_s}$

Pour des valeurs données de $Z_s$ et de $c$ la forme de la flamme est donnée par la condition $Z(\xi ,\eta )=Z_s$, c'est-à-dire :

$Z_s=c^2+2c\sum _{n=1}^{\mathrm{\infty }}\frac{1}{{\lambda }_n}\frac{J_1(c{\lambda }_n)}{J_0^2({\lambda }_n)}{J}_0({\lambda }_n\xi )e^{-{\lambda }_n^2\eta }.$

Lorsque $Z_s\to 0$ ($S\to \mathrm{\infty }$), la flamme s'étend de l'embouchure du tube intérieur et se fixe au tube extérieur à une certaine hauteur (cas sous-ventilé) et lorsque $Z_s\to 1$ ($S\to 0$), la flamme part de l'embouchure du tube intérieur et se joint à l'axe à une certaine hauteur éloignée de l'embouchure (cas sur-ventilé). En général, la hauteur de la flamme est obtenue en résolvant $\eta$ dans l'équation ci-dessus après avoir défini $\xi =1$ pour le cas sous-ventilé et $\xi =0$ pour le cas sur-ventilé.

Étant donné que les hauteurs de flamme sont généralement assez grandes la hauteur de la flamme peut être estimée en ne conservant que le premier terme de la série du développement des fonctions de Bessel. Cette approximation prédit les hauteurs de flamme pour les deux cas comme suit :

$\begin{array}{cc}\text{cas sous-ventilé}& \eta =\frac{1}{{\lambda }_1^2}\ln \left[\frac{2c{J}_1(c{\lambda }_1)}{(Z_s-c^2){\lambda }_1{J}_0({\lambda }_1)}\right]\\ \text{cas sur-ventilé}& \eta =\frac{1}{{\lambda }_1^2}\ln \left[\frac{2c{J}_1(c{\lambda }_1)}{(Z_s-c^2){\lambda }_1{J}_0^2({\lambda }_1)}\right]\end{array}$

où ${\lambda }_1=3,8317.$