Articulo de referencia

ecuación G

En la combustión , la ecuación G es un escalar. GRAMO ( incógnita , t ) {\displaystyle G(\mathbf {x} ,t)} Ecuación de campo que describe la posición instantánea de la llama, int...

En la combustión , la ecuación G es un escalar.GRAMO(incógnita,t){\displaystyle G(\mathbf {x} ,t)}Ecuación de campo que describe la posición instantánea de la llama, introducida por Forman A. Williams en 1985 [ 1 ] [ 2 ] en el estudio de la combustión turbulenta premezclada. La ecuación se deriva a partir del método de conjunto de nivel . La ecuación fue estudiada por primera vez por George H. Markstein , en una forma restrictiva para la velocidad de combustión y no como un conjunto de nivel de un campo. [ 3 ] [ 4 ] [ 5 ] La ecuación G se lee [ 6 ]

ρ(GRAMOt+vGRAMO)=metro˙|GRAMO|{\displaystyle \rho \left({\frac {\partial G}{\partial t}}+\mathbf {v} \cdot \nabla G\right)={\dot {m}}|\nabla G|}

dóndeρ{\displaystyle \rho }es la densidad del flujo,v{\displaystyle v}es la velocidad del flujo ymetro˙=metro˙(incógnita,t){\displaystyle {\dot {m}}={\dot {m}}(\mathbf {x} ,t)}es el flujo de masa normal que ingresa a cualquier conjunto de niveles particularGRAMO(incógnita,t)={\displaystyle G(\mathbf {x} ,t)=}constante.

Descripción matemática

La ecuación G se lee como [ 7 ] [ 8 ]

GRAMOt+vGRAMO=ST|GRAMO|{\displaystyle {\frac {\partial G}{\partial t}}+\mathbf {v} \cdot \nabla G=S_{T}|\nabla G|}

dónde

  • v{\displaystyle \mathbf {v} }es el campo de velocidad del flujo ,
  • ST=metro˙/ρ{\displaystyle S_{T}={\dot {m}}/\rho _{u}}es la velocidad de combustión local con respecto al gas no quemado con densidadρ{\displaystyle \rho _{u}}.

La ubicación de la llama viene dada porGRAMO(incógnita,t)=GRAMOo{\displaystyle G(\mathbf {x} ,t)=G_{o}}que puede definirse arbitrariamente de tal manera queGRAMO(incógnita,t)>GRAMOo{\displaystyle G(\mathbf {x} ,t)>G_{o}}es la región de gas quemado yGRAMO(incógnita,t)<GRAMOo{\displaystyle G(\mathbf {x} ,t)<G_{o}}es la región de gas no quemado. El vector normal a la llama, apuntando hacia el gas quemado, esnorte=GRAMO/|GRAMO|{\displaystyle \mathbf {n} =\nabla G/|\nabla G|}.

La ecuación G tiene la forma de la ecuación de Hamilton-Jacobi , una ecuación de mecánica analítica utilizada para modelar la dinámica de partículas como propagación de ondas. [ 9 ] [ 10 ]

Velocidad de combustión local

Según la teoría de Matalon-Matkowsky-Clavin-Joulin , la velocidad de combustión de la llama estirada , para una curvatura y deformación pequeñas, viene dada por

ST=SL+METROdoδL(SLvnorte)norteMETROtδLtvt{\displaystyle S_{T}=S_{L}+{\mathcal {M}}_{c}\delta _{L}(S_{L}-\mathbf {v} \cdot \mathbf {n} )\nabla \cdot \mathbf {n} -{\mathcal {M}}_{t}\delta _{L}\nabla _{t}\cdot \mathbf {v} _ {t}}

dónde

  • SL{\displaystyle S_{L}}es la velocidad de combustión de la llama no estirada con respecto al gas no quemado
  • METROdo{\displaystyle {\mathcal {M}}_{c}}yMETROt{\displaystyle {\mathcal {M}}_{t}}son los dos números de Markstein , asociados con la curvatura y la deformación tangencial;tvt=nortenorte:v(vnorte)norte{\displaystyle \nabla _ {t}\cdot \mathbf {v} _ {t}=-\mathbf {n} \otimes \mathbf {n} :\nabla \mathbf {v} -(\mathbf {v} \cdot \mathbf {n} )\nabla \cdot \mathbf {n} } es la divergencia superficial de la velocidad tangencialvt=(Inortenorte)v{\displaystyle \mathbf {v} _ {t}=(\mathbf {I} -\mathbf {n} \otimes \mathbf {n} )\mathbf {v} }
  • δL{\displaystyle \delta _{L}}es el espesor de la llama laminar.

Un ejemplo sencillo: quemador de ranura

quemador de ranura

La ecuación G tiene una expresión exacta para un quemador de ranura simple. Consideremos un quemador de ranura planar bidimensional de ancho de ranurab{\displaystyle b}La mezcla de reactivos premezclada se alimenta a través de la ranura desde la parte inferior a una velocidad constante.v=(0,U){\displaystyle \mathbf {v} =(0,U)}donde la coordenada(incógnita,y){\displaystyle (x,y)}se elige de tal manera queincógnita=0{\displaystyle x=0}se encuentra en el centro de la ranura yy=0{\displaystyle y=0}se encuentra en la ubicación de la boca de la ranura. Cuando se enciende la mezcla, se desarrolla una llama premezclada desde la boca de la ranura hasta cierta altura.y=L{\displaystyle y=L}en forma de cuña bidimensional con un ángulo de cuñaα{\displaystyle \alpha }Para simplificar, supongamos que...ST=SL{\displaystyle S_{T}=S_{L}}, lo cual es una buena aproximación excepto cerca de la esquina de la cuña donde los efectos de curvatura se vuelven importantes. En el caso estacionario, la ecuación G se reduce a

UGRAMOy=SL(GRAMOincógnita)2+(GRAMOy)2{\displaystyle U{\frac {\partial G}{\partial y}}=S_{L}{\sqrt {\left({\frac {\partial G}{\partial x}}\right)^{2}+\left({\frac {\partial G}{\partial y}}\right)^{2}}}}

Si una separación de la formaGRAMO(incógnita,y)=y+F(incógnita){\displaystyle G(x,y)=y+f(x)}se introduce, entonces la ecuación se convierte en

U=SL1+(Fincógnita)2,Fincógnita=U2SL2SL{\displaystyle U=S_{L}{\sqrt {1+\left({\frac {\partial f}{\partial x}}\right)^{2}}},\quad \Rightarrow \quad {\frac {\partial f}{\partial x}}={\frac {\sqrt {U^{2}-S_{L}^{2}}}{S_{L}}}}

que al integrarse da

F(incógnita)=(U2SL2)1/2SL|incógnita|+do,GRAMO(incógnita,y)=(U2SL2)1/2SL|incógnita|+y+do{\displaystyle f(x)={\frac {\left(U^{2}-S_{L}^{2}\right)^{1/2}}{S_{L}}}|x|+C,\quad \Rightarrow \quad G(x,y)={\frac {\left(U^{2}-S_{L}^{2}\right)^{1/2}}{S_{L}}}|x|+y+C}

Sin pérdida de generalidad, elija la ubicación de la llama enGRAMO(incógnita,y)=GRAMOo=0{\displaystyle G(x,y)=G_{o}=0}Dado que la llama está unida a la boca de la ranura|incógnita|=b/2, y=0{\displaystyle |x|=b/2,\ y=0}, la condición de contorno esGRAMO(b/2,0)=0{\displaystyle G(b/2,0)=0}, que puede utilizarse para evaluar la constantedo{\displaystyle C}. Por lo tanto, el campo escalar es

GRAMO(incógnita,y)=(U2SL2)1/2SL(|incógnita|b2)+y{\displaystyle G(x,y)={\frac {\left(U^{2}-S_{L}^{2}\right)^{1/2}}{S_{L}}}\left(|x|-{\frac {b}{2}}\right)+y}

En la punta de la llama, tenemosincógnita=0, y=L, GRAMO=0{\displaystyle x=0,\ y=L,\ G=0}, que nos permiten determinar la altura de la llama

L=b(U2SL2)1/22SL{\displaystyle L={\frac {b\left(U^{2}-S_{L}^{2}\right)^{1/2}}{2S_{L}}}}

y el ángulo de la llamaα{\displaystyle \alpha },

broncearseα=b/2L=ST(U2SL2)1/2{\displaystyle \tan \alpha ={\frac {b/2}{L}}={\frac {S_{T}}{\left(U^{2}-S_{L}^{2}\right)^{1/2}}}}

Utilizando la identidad trigonométricabroncearse2α=pecado2α/(1pecado2α){\displaystyle \tan ^{2}\alpha =\sin ^{2}\alpha /\left(1-\sin ^{2}\alpha \right)}, tenemos

pecadoα=SLU.{\displaystyle \sin \alpha ={\frac {S_{L}}{U}}.}

De hecho, la fórmula anterior se utiliza a menudo para determinar la velocidad de combustión planar.SL{\displaystyle S_{L}}, midiendo el ángulo de la cuña.

Referencias

  1. Williams, FA (1985). Combustión turbulenta. En Las matemáticas de la combustión (pp. 97-131). Sociedad de Matemáticas Industriales y Aplicadas.
  2. Kerstein, Alan R., William T. Ashurst y Forman A. Williams. "Ecuación de campo para la propagación de la interfaz en un campo de flujo homogéneo no estacionario". Physical Review A 37.7 (1988): 2728.
  3. GH Markstein. (1951). Interacción de pulsaciones de flujo y propagación de llama. Journal of the Aeronautical Sciences, 18(6), 428-429.
  4. Markstein, GH (Ed.). (2014). Propagación de llamas no estacionarias: AGARDograph (Vol. 75). Elsevier.
  5. Markstein, GH, & Squire, W. (1955). Sobre la estabilidad de un frente de llama plano en flujo oscilante. The Journal of the Acoustical Society of America, 27(3), 416-424.
  6. Rajamanickam, P., & Daou, J. (2024). Teoría hidrodinámica de llamas premezcladas bajo la ley de Darcy. Física de fluidos, 36(12).
  7. Peters, Norbert. Combustión turbulenta. Cambridge University Press, 2000.
  8. Williams, Forman A. "Teoría de la combustión." (1985).
  9. Xin, J.; Yu, Y.; Ronney, P. Métodos lagrangianos, de teoría de juegos y de EDP para promediar ecuaciones G en combustión turbulenta: existencia y más allá . arXiv:2401.14575 [math.AP], 2024.
  10. G. Pagnini, RAD Akkermans, N. Buchmann, A. Mentrelli, Enfoques de la ecuación de reacción-difusión y la ecuación G reconciliados en el modelado de la combustión premezclada turbulenta , XXXVIII Reunión de la Sección Italiana del Instituto de Combustión, BCAM – Centro Vasco de Matemáticas Aplicadas (Bilbao, España); Universidad Técnica de Braunschweig (Alemania); Universidad de Bolonia (Italia).