Integración de Ecuaciones Diferenciales rígidas de valor de contorno en problemas de combustión con cinética de reacción detallada
Modelos detallados de la llama de difusión plana, laminar de jets opuestos requieren la solución de las ecuaciones de conservación de cantidad de movimiento, energía y especies químicas en conjunto con una cinética de reacción química stiff. El problema tiene soluciones autosimilares y se puede reso...
| Autores: | , |
|---|---|
| Tipo de recurso: | artículo |
| Fecha de publicación: | 1988 |
| País: | España |
| Institución: | Universitat Politècnica de Catalunya (UPC) |
| Repositorio: | UPCommons. Portal del coneixement obert de la UPC |
| Idioma: | español |
| OAI Identifier: | oai:upcommons.upc.edu:2099/7450 |
| Acceso en línea: | https://hdl.handle.net/2099/7450 |
| Access Level: | acceso abierto |
| Palabra clave: | Numerical Methods Equacions diferencials Àrees temàtiques de la UPC::Matemàtiques i estadística::Equacions diferencials i integrals |
| Sumario: | Modelos detallados de la llama de difusión plana, laminar de jets opuestos requieren la solución de las ecuaciones de conservación de cantidad de movimiento, energía y especies químicas en conjunto con una cinética de reacción química stiff. El problema tiene soluciones autosimilares y se puede resolver mediante integración numérica de un conjunto de ecuaciones diferenciales ordinarias stiff de segundo orden, con condiciones de contorno, cada una de las cuales tiene un punto de torno de primer orden debido a la ecuación de conservación de cantidad de movimiento. El uso de discretización en diferencias finitas standard (en el dominio espacial) y la expansión en serie de Taylor alrededor de la iteración anterior (en el espacio temporal) del término cinético de las reacciones químicas lleva a una ecuación de matrices en la cual la matriz de coeficientes es del tipo tridiagonal en bloques y muy grande de tamaño. También representa una matriz de banda y se puede obtener una solución al problema mediante una descomposición LU. Inestabilidades numéricas debidas a la malla finita de espaciamiento irregular y a la difusión de errores numéricos hacia los bordes forzaron el uso de una malla regular con un gran número de puntos. Esto restringió el estudio a la solución de problemas con esquemas cinéticos relativamente simples, debido a la alta necesidad de memoria de computador. Este dilema fue resuelto por el desarrollo de un esquema de diferenciación modificado, que asume que la solución puntual del problema está dada por una suma de exponenciales. A Usando esta nueva técnica fue posible hallar la solución de una llama con 150 reacciones envolviendo las ecuaciones de ccnservación de 70 especies químicas. En este artículo formulamos el problema y describimos resultados de los varios esquemas numéricos mediante los que se intentó resolver el mismo. El foco principal, sin embargo, está u dirigido hacia la modificación del esquema de diferencias finitas, acoplado con la expansión en serie de Taylor de los términos de reacción química. |
|---|