Source: http://docplayer.es/1427628-Simulacion-numerica-de-la-distribucion-de-temperaturas-y-flujo-de-aire-en-el-horno-de-una-cocina-de-uso-domestico.html
Timestamp: 2018-04-24 01:41:21+00:00

Document:
Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico - PDF
Download "Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico"
Blanca Torregrosa Fuentes
1 Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico
3 UNIVERSIDAD POLITÉCNICA SALESIANA SEDE CUENCA UNIDAD DE POSGRADOS MAESTRÍA EN MÉTODOS NUMÉRICOS PARA DISEÑO EN INGENIERÍA Tema: Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico Tesis previa a la obtención del Grado de: Magíster en Métodos Numéricos para Diseño en Ingeniería Autor: Ing. Olena Neira Director: MSc. Ing. Paúl Álvarez Cuenca, 1 de junio del 2011
4 Datos de catalogación bibliográfica OLENA NEIRA RODAS Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico. Universidad Politécnica Salesiana Cuenca-Ecuador, 2011 DEPERTAMENTO DE POSGRADO Formato 170x240 Páginas 104 Breve reseña del autor e información de contacto Olena Neira Rodas Ingeniero Mecánico Egresada de la Maestria en Métodos Numéricos para Diseño en Ingeniería Unidad de Posgrados Universidad Politécnica Salesiana Dirigido por: Paúl Alvarez Lloret Ingeniero Mecánico Master en Métodos Numéricos Docente de la Universidad Politécnica Salesiana Todos los derechos reservados. Queda prohibida, salvo excepción prevista en la Ley, cualquier forma de reproducción, distribución, comunicación pública y transformación de esta obra para fines comerciales, sin contar con autorización de los titulares de propiedad intelectual. La infracción de los derechos mencionados puede ser constitutiva de delito contra la propiedad intelectual. Se permite la libre difusión de este texto con fines académicos o investigativos por cualquier medio, con la debida notificación a los autores. DERECHOS RESERVADOS 2011 Universidad Politécnica Salesiana. Cuenca-Ecuador-Sudamérica OLENA NEIRA RODAS Simulación numérica de la distribución de temperaturas y flujo de aire en el horno de una cocina de uso doméstico. IMPRESO EN ECUADOR - PRINTED IN ECUADOR
5 ÍNDICE GENERAL CAPÍTULO I. INTRODUCCIÓN Las técnicas numéricas en Ingeniería Antecedentes Justificación CAPÍTULO II. PLANTEAMIENTO DEL PROBLEMA Problemas existentes en los hornos de cocinas de uso doméstico Necesidad de solución numérica Descripción del sistema analizado Identificación del dominio Selección de software CAPÍTULO III. PREPROCESO Introducción Creación de una geometría tridimensional Discretización espacial Desarrollo del modelo matemático Introducción Ecuaciones de conservación Modelo de turbulencia Las ecuaciones de conservación de especies Modelo de combustión Modelo de radiación Condiciones iniciales y condiciones de contorno CAPÍTULO IV. RESOLUCIÓN DE SISTEMAS DE ECUACIONES DE GOBIERNO Método de solución Método de las diferencias finitas Método de los elementos finitos Método de los volúmenes finitos CAPÍTULO V. POST-PROCESO Análisis de resultados de simulaciones Validación de resultados... 92
6 CAPÍTULO VI. CONCLUSIONES Y RECOMENDACIONES Conclusiones Recomendaciones BIBLIOGRAFÍA ÍNDICE DE FIGURAS Figura 1. Esquema de funcionamiento del horno de la cocina de uso doméstico Figura 2. Quemador del horno de la cocina de uso doméstico Figura 3 a. Cocina completa Figura 3 b. Volumen ocupado por el fluido en diferentes zonas de la cocina Figura 3 c. Quemador Figura 3 d. Calientaplatos Figura 3 e. Horno Figura 4. División del dominio en subdominios Figura 5. Elementos usados para mallar geometrías bidimensionales Figura 6. Elementos usados para mallar geometrías tridimensionales Figura 7. Ejemplo de mallado body-fitted Figura 8. Ejemplo de mallado en coordenadas cartesianas Figura 9. Mallas con nodos centrados en los elementos y centrados en los vértices.. 26 Figura 10 a,b,c. Quemador mallado Figura 11 a, b, c. Calientaplatos mallado Figura 12 a,b. Horno mallado Figura 13. Flujo másico que entra y sale de un elemento fluido Figura 14. Esfuerzos superficiales sobre las caras de un elemento fluido Figura 15. (a), (b). Esfuerzos normales y tangenciales actuando sobre un elemento fluido Figura 16. Vectores de calor de un elemento fluido Figura 17. Discretización de variables Figura 18 a), b), c): Campo de velocidad de flujo de mezcla butano-aire primario dentro del quemador Figura 19. El rango de velocidad del flujo dentro del quemador Figura 20. La forma de las llamas. a). Resultado de la simulación. b). La imagen real Figura 21. Distribución de la temperatura en la bandeja Figura 22. Distribución de la temperatura: a) y b) en el interior del calientaplatos en los planos paralelos a la bandeja c) en la bandeja, obtenida experimentalmente... 89
7 Figura 23. El rango de velocidad y temperatura en calientaplatos Figura 24. El campo de velocidad en el calientaplados: a) en los orificios de salida; b) y c) en los planos transversal y longitudinal del calientaplatos en los planos paralelos a la bandeja Figura 25. Distribución de la temperatura en el interior del horno en los planos paralelos a la bandeja Figura 26. Flujo de gases de combustión en el horno Figura 27. Rango de temperatura en el horno Figura 28. Fuentes de error en los cálculos CFD Figura 29. La posición de termocuplas utilizada para medir la temperatura dentro del horno ÍNDICE DE TABLAS Y GRÁFICAS Tabla 1. Parámetros seleccionados para el mallado del quemador Tabla 2. Parámetros seleccionados para el mallado del calientaplatos Tabla 3. Parámetros seleccionados para el mallado del horno Tabla 4. Datos informativos de la malla del quemador Tabla 5. Datos informativos de la malla del calientaplatos Tabla 6. Datos informativos de la malla del horno Tabla 7. Ecuaciones de conservación Tabla 8. Ecuaciones promediadas de Navier-Stokes, RANS Tabla 9. Coeficientes constantes. Modelo k-ε Tabla 10. Comparación de los modelos de radiación Tabla 11. Resumen de los modelos empleados en las simulaciones Tabla 12. Tipos clásicos de condiciones frontera, que se presentan en los problemas de fenómenos de transporte Tabla 13. Los valores de la temperatura el los puntos indicados en la figura 28, obtenida experimentalmente y mediante la simulación Gráfica 1. La velocidad de salida de la mezcla butano-aire primario en los centros de los agujeros del quemador en la parte frontal y posterior Gráfica 2. La velocidad de salida de los gases de combustión en los centros de los agujeros de la bandeja en la parte izquierda y derecha Gráfica 3 a) y b). Comparación de la temperatura de los gases de combustión medida en los puntos indicados en la figura del horno y obtenida mediante la simulación... 97
9 CAPÍTULO I INTRODUCCIÓN 1.1. Las técnicas numéricas en Ingeniería [8] Las técnicas numéricas en Ingeniería han experimentado un gran desarrollo en las últimas décadas. De esta tendencia no se ha apartado una rama tan característica de la Ingeniería Mecánica como es la Mecánica de Fluidos. La dinámica de fluidos computacional (o CFD, acrónimo de las palabras inglesas Computational Fluid Dynamics ) consiste en el análisis del movimiento de los fluidos mediante simulaciones con ordenadores. Su objetivo es la búsqueda de una solución aproximada de las ecuaciones que gobiernan el movimiento de los fluidos, discretizando o dividiendo el dominio de cálculo en pequeños elementos y resolviendo allí dichas ecuaciones. Los métodos numéricos aplicados a la Mecánica de Fluidos resultan una herramienta muy útil para el diseño y análisis de las distintas situaciones prácticas en las que se utilizan los fluidos. Hoy en día, con las enormes posibilidades de los ordenadores y los desarrollos de los métodos numéricos, la CFD se está convirtiendo en una herramienta muy práctica y eficiente para el análisis de fluidos. La CFD puede ser muy útil en un amplio espectro de industrias y representa una poderosa ayuda para el ingeniero de diseño, de producción e incluso de mantenimiento. Se está utilizando en sectores tan variados como la industria química, aeroespacial, metalúrgica, eléctrica, etc. Es también ampliamente utilizada en el estudio para predecir el movimiento de aire y el confort (sistemas de ventilación), así como para analizar fenómenos termofluídicos en el medio ambiente, como dispersión de contaminantes en la atmósfera o de vertidos en los ríos, etc. 9
10 Básicamente son tres los tipos de problemas que se pueden resolver: Simulación de un equipo existente, con el fin de evaluar cambios en su operación o en su diseño (ahorros de energía, mejoras en la calidad del producto, mejoras medioambientales), o para diagnosticar problemas operacionales. Mejoras en el diseño de equipos, pudiéndose evaluar un amplio abanico de opciones de diseño, incrementándose así las posibilidades para la innovación tecnológica. Simulación de procesos, con o sin transiciones de fase, interacciones sólidofluido, fluido-medio poroso. El software existente de Dinámica de Fluidos Computacional puede ser utilizado para la simulación numérica de los procesos de flujo de fluidos, mezclado, combustión, transporte de calor y masa, etc. Las principales aplicaciones de la CFD en equipos y sistemas son: Diagnósticos con predicción de flujos de calor, campos de velocidades y temperaturas en conductos de gases, tuberías de líquidos, cambiadores de calor, etc. Diagnósticos con predicción de flujos de calor, campos de velocidades y temperaturas en cámaras de combustión, calderas, hornos, etc. Análisis de flujos en ciclones, separadores, etc. Diseño de componentes de plantas de tratamiento de aguas residuales y aguas limpias. Diseño de quemadores en general, quemadores de bajo NOx, válvulas y otros equipos auxiliares Antecedentes Hasta finales de los años 60 los ordenadores no alcanzaron velocidades de cálculo suficientes como para resolver casos sencillos. Hasta entonces, las técnicas experimentales constituían prácticamente la única herramienta de análisis y diseño de cualquier problema de Mecánica de Fluidos. En la actualidad, los ensayos experimentales siguen siendo necesarios para la comprobación de las prestaciones de diseños complejos, pero los continuos avances en la 10
11 capacidad de los ordenadores y en los algoritmos permiten una reducción importante en el número de ensayos necesarios. Así, por ejemplo, el diseño típico de un modelo de ala de avión, exige en la actualidad 3 o 4 ensayos en túnel aerodinámico, en lugar de los 10 o 15 que eran necesarios anteriormente. En realidad, la palabra que mejor definiría hoy la relación entre ambas herramientas podría ser la de complementariedad. A lo largo de los últimos veinte años las técnicas de CFD han evolucionado mejorando los programas comerciales e introduciéndose en las distintas áreas de la ingeniería hasta obtener un merecido lugar dentro de las necesidades reales de la industria. Dichos programas se vienen usando de manera creciente de manera paralela a la mejora de los sistemas informáticos que les sirven de soporte. Desde su inicial concepción, orientada a la industria aeroespacial, la técnicas numéricas se han extendido a un número creciente de aplicaciones dentro de un amplio espectro de industrias, desde las más clásicas como la automovilística o electrónica, hasta las nuevas aplicaciones en industrias alimentarias y biomédicas. Sin embargo, aún no se emplean las técnicas numéricas como auténticas herramientas de diseño. En realidad, se obtendrán los mayores beneficios cuando se llegue a un nivel de utilización en el día a día del diseño industrial, es decir, cuando su difusión sea tal que puedan llegar a ser utilizadas por ingenieros sin demasiada especialización en estas técnicas y no estar restringida su utilización a los expertos en la materia Justificación El uso de prototipos para el diseño de productos es mucho más costoso que la construcción de modelos computacionales. Por otro lado, el control que se tiene sobre los experimentos virtuales es mucho mayor que el que se tiene sobre los experimentos reales. Además, la cantidad de información que se obtiene de las simulaciones permite entender con mayor detalle los procesos que se desarrollan durante el funcionamiento de un producto industrial. El conocimiento de los procesos térmicos que se desarrollan en los hornos de las cocinas se vuelve indispensable para la toma de decisiones inteligentes sobre ciertos parámetros que regulan su funcionamiento. El uso de prototipos y de experimentos, aunque siempre es deseable en alguna medida, no responde a las necesidades modernas de eficiencia y precisión. Por esta razón, Los procesos térmicos están caracterizados por la distribución de temperaturas y los correspondientes flujos de calor que se producen en el transcurso del tiempo. La determinación de dicha distribución exige la resolución de un problema de valores iniciales y de contorno (PVIC), gobernado por ecuaciones en derivadas parciales. Estas ecuaciones provienen del balance de energía en el medio estudiado. En la formulación del modelo matemático se pueden considerar diversos efectos de transmisión del calor: conducción, convección, radiación. 11
12 La resolución del PVIC se realizará mediante el uso de métodos numéricos, esto es, se transformará el modelo matemático en un modelo numérico mediante un proceso de discretización. La discretización en el espacio se realizará mediante el método de los volúmenes finitos y la discretización en el tiempo mediante diferencias finitas. Vale indicar que en este proyecto se recurrirá a códigos comerciales que tienen ya implementados los modelos numéricos descritos. Las ventajas que proporciona el análisis con técnicas CFD se pueden resumir en las siguientes: Reducción sustancial de tiempos y de costes en los nuevos diseños. Posibilidad de analizar sistemas y condiciones muy difíciles de simular experimentalmente: velocidades supersónicas, temperaturas extremas y elementos en movimiento relativo. Capacidad de estudiar sistemas bajo condiciones peligrosas o más allá de sus condiciones límite de funcionamiento, por ejemplo, accidentes con sustancias tóxicas. Nivel de detalle prácticamente ilimitado. Los métodos experimentales son tanto más caros cuanto mayor es el número de puntos de medida, mientras que los programas CFD pueden generar un gran volumen de resultados sin coste añadido, y con posibilidad hacer estudios paramétricos. Un valor añadido es poder poner en el producto la etiqueta de Diseñado con ayuda del ordenador, y la facilidad para generar gráficos fácilmente interpretables, que estimulan la compra del producto (lo cual, por otro lado, constituye un riesgo). Aunque se razone que se abaratan los costes respecto a las técnicas experimentales, las técnicas CFD no son gratuitas. En primer lugar, son necesarias máquinas de gran capacidad de cálculo (los usuarios de técnicas CFD trabajan habitualmente con los ordenadores más potentes que existen), y un software con precio todavía no accesible al gran público. En segundo lugar, se necesita personal cualificado que sea capaz de dominar los programas y, sobre todo, analizar adecuadamente los resultados. Los desarrollos en el campo de las técnicas numéricas dedicadas al estudio de flujos se están acercando cada vez más a los de otras herramientas de CAE como las de análisis de esfuerzos en sólidos y estructuras. Su retraso se debe a la gran complejidad de las ecuaciones y, sobre todo, a la dificultad de modelar adecuadamente ciertos fenómenos como la turbulencia, los flujos multifásicos y la combustión. 12
13 Uno de los mayores inconvenientes de las técnicas CFD consiste en que no siempre es posible llegar a obtener resultados suficientemente precisos y siempre está presente la posibilidad de cometer graves errores en cuestiones básicas. Esto proviene de: Simplificar el fenómeno que será estudiado para que el hardware y software sean capaces de abordarlo. El resultado será tanto más preciso cuanto más adecuadas hayan sido las hipótesis y las simplificaciones realizadas. La existencia de insuficientes e incompletos modelos para la simular el efecto de la turbulencia, los flujos multifásicos o la combustión, entre otros. La tendencia humana de aceptar sin mayores reservas todo lo que se ha obtenido utilizando un ordenador, sobre todo cuando se presentan los resultados en forma atractiva. 13
14 CAPÍTULO II PLANTEAMIENTO DEL PROBLEMA 2.1. Problemas existentes en los hornos de cocinas de uso doméstico En las cocinas de uso doméstico en la actualidad existe el problema de irregularidad de distribución de temperaturas dentro de los hornos, lo que afecta la calidad de los productos preparados en ellos, sobre todo productos de pastelería. Este problema en algunos hornos se soluciona con el uso de un ventilador, pero esta solución se refleja en los costos de las cocinas, aumento considerable del tiempo de cocción y consumo de gas. Las causas de este comportamiento son desconocidas y los fabricantes de las cocinas intentan solucionar estos problemas mediante cambios en la geometría del horno y del quemador o cambiando los parámetros de funcionamiento como por ejemplo el caudal del gas y del aire primario. Pero este procedimiento es largo y muy costoso y hasta el momento no se ha logrado solucionar el problema Necesidad de solución numérica El diseño de sistemas de combustión y de control de contaminantes tradicionalmente se ha basado en métodos empíricos: experiencia y experimentación práctica. Los modelos físicos y matemáticos que se han usado han sido de un rango de aplicación limitado y de formulación básica. El tema ha sido tratado fundamentalmente desde el punto de vista tecnológico, donde la cuestión principal es el procedimiento técnico o cómo hacerlo. Paralelamente, ha habido un lento desarrollo de la ciencia de la combustión, y algún progreso en contestar los por qué más que los cómo. Recién en los últimos años se ha producido un progreso significativo con el desarrollo de descripciones matemáticas más completas como modelos de fenómenos de combustión. Estos modelos ya permiten 14
15 el análisis detallado del flujo, la transferencia de calor y la distribución espacial de las velocidades de reacción dentro de un sistema, aunque aún no son lo suficientemente versátiles, ni han sido validados experimentalmente a un grado suficiente, como para permitir a un ingeniero especialista en combustión que los utilice para optimizar un combustor o planta de quemado sin recurrir a datos experimentales, fórmulas empíricas o experiencia previa. Sin embargo, estos modelos resultan muy útiles en la etapa de diseño, pues permiten analizar rápida y económicamente la influencia de los cambios en la geometría, las condiciones de entrada, etc, en forma paramétrica. Los métodos tradicionales utilizados para un diseño de hornos aumenta considerablemente el tiempo necesario para mejorar su calidad y los costos de los mismos. La aplicación de métodos numéricos para el diseño permitirá soluciones optimas en tiempos mas cortos y con costos mas bajos. A diferencia de los métodos experimentales CFD proporciona datos de presión velocidad, densidad, temperatura y otros valores de interés en todos los puntos del dominio de flujo simulado. Todo ello hace que el modelado matemático del funcionamiento del horno resulte una herramienta muy interesante en el diseño o rediseño de procesos Descripción del sistema analizado. Para el análisis se escogió la cocina a gas de uso doméstico de marca Indurama modelo Avant 32. Construcción y funcionamiento del horno de la cocina. El horno de la cocina de uso doméstico consta de dos compartimientos (figura 1): la parte inferior es la cámara de combustión, llamada calientaplatos, donde se realiza el proceso de combustión, y la parte superior, llamada horno, donde se realiza la preparación de los alimentos. Las dos partes están separadas por una bandeja con orificios por los cuales los gases calientes pasan del calientaplatos al horno. 15
16 Figura 1. Esquema de funcionamiento del horno de la cocina de uso doméstico. Dentro del calientaplatos está ubicado el quemador atmosférico, llamado quemador Bunsen. El quemador (figura 2) se compone del inyector de gas 1, en el cual en función de la presión de la línea de combustible y de la configuración del agujero de salida, induce el combustible y aire primario en un vénturi (formado por el cono 2, cuello 3 y difusor 4), seguido de la cabeza 5 donde la mezcla reactiva es distribuida a los pequeños agujeros de salida, en los cuales las llamas son estabilizadas. El proceso de inducción de aire primario acontece debido a la acción combinada de efectos difusivos y advectivos ejercidos por el combustible que sale del inyector a presión en el aire en torno a él. Estos efectos fluido-dinámicos dependen de las condiciones de presión, densidad del combustible y geometría del inyector y del vénturi. La estrangulación del fluido a su paso por el vénturi genera una caída de presión axial, induciendo aire al interior del quemador. Un segundo efecto de estrangulación en la garganta 3 del vénturi es la homogeneización espacial de la estequiometría de la mezcla reactiva en función de los valores de velocidad y turbulencia alcanzados. Después de la garganta del vénturi, el fluido comienza a recuperar presión estática debido a la expansión del área de flujo, hasta un valor superior al de la presión atmosférica en un punto localizado antes de los pequeños agujeros, lo anterior permite la salida de la mezcla reactiva por los agujeros en los cuales se forman las llamas. La mezcla aire / gas que circula por el tubo debe tener una velocidad tal que impida que la llama que se forma en el extremo superior después de la ignición inicial, penetre en el interior del quemador. La reacción de combustión no finaliza en la región del cono, existe un manto luminoso exterior al cono donde se completa la combustión del gas residual por medio de la difusión del 16
17 aire secundario dentro de ese manto. La tuerca 7 sirve para regular la cantidad de aire primario inducido. llamas 1 aire gas aire 7 Figura 2: Quemador del horno de la cocina de uso doméstico. 1-inyector (ciclor), 2- cono, 3- cuello, 4-difusor, 5-cabeza con los agujeros. (2,3,4 forman el mezclador o vénturi) 7-tuerca reguladora del caudal de aire primario. Las llamas, en este caso, son parcialmente premezcladas porque el gas entra a la cámara de combustión mezclado con una parte de aire, necesario para la combustión y otra parte, entra al calientaplatos por los agujeros que se encuentran en la parte posterior e inferior. Para iniciar este trabajo el fabricante proporcionó las dimensiones de la cocina, el consumo de gas que depende del tipo del inyector utilizado y la conductividad térmica del material de las paredes y de la bandeja. Datos iniciales: Flujo volumétrico del gas-butano: Qg = m3 /h La conductividad de las paredes y de la bandeja: k = 0,075 W/m*K 17
18 Cantidades por determinar Mediante la solución numérica se pretende determinar las temperaturas y direcciones de flujo de gases calientes dentro del horno. Luego esta información puede ser utilizada para modificar las dimensiones, formas, parámetros de funcionamiento y buscar las soluciones óptimas Identificación del dominio. Para simplificar el problema es necesario aislar el horno del sistema completo (cocina completa), pero esto debe hacerse de tal forma que podamos definir perfectamente las condiciones de contorno. Para la simulación del horno en condiciones de contorno es necesario saber la temperatura, velocidad y dirección de gas caliente que sale del calientaplatos, y de esta información no disponemos, por eso no podemos aislar el horno del calientaplatos. Para el análisis del calientaplatos tenemos que saber la velocidad, dirección y composición de la mezcla gas-aire que sale del quemador, la velocidad y la dirección no sabemos, por eso tampoco podemos aislar el horno y calientaplatos del quemador. Para el análisis del quemador necesitamos el caudal del gas, y la dirección de entrada del gas al quemador y esta información es proporcionada por el fabricante. Aislamos nuestro sistema quemador-calientaplatos-horno del resto de partes de la cocina. También prescindimos de las paredes, porque no nos interesa la distribución de la temperatura dentro de ellas. Considerando lo anterior nuestro sistema tiene que empezar desde la entrada del gas al quemador y terminar en la salida de gas de combustión del horno. En la simulación del desplazamiento de la mezcla gas-aire dentro del quemador y del desplazamiento de gases de combustión dentro del horno el dominio es el fluido; en la simulación del proceso de combustión el dominio es el fluido reactante. Número de dimensiones del modelo A pesar de que los problemas uni y bidimensionales son mucho más fáciles de abordar, en nuestro caso estos modelos se separarían mucho de la realidad, debido a que la irregularidad de temperatura se da en tres direcciones espaciales. Por las mismas razones no se puede usar simetría axial. Modelos físicos incluidos en el análisis En el análisis de distribución de temperaturas y de direcciones de flujo de gases de combustión se incluirá el modelo de turbulencia, de combustión y de radiación. 18
19 La hipótesis Puesto que en el análisis de las condiciones de operación de un horno algunos detalles pueden ser obviados y a fin de obtener un modelo matemático cuya resolución numérica sea abordable con tiempos de cálculo que no pasen de unas pocas horas, es necesario hacer varias hipótesis simplificadoras. 1). Las temperaturas de gases y de las paredes del horno son estacionarias. 2). El fluido se considera incompresible (el número de Mach<1). 3). La variación de la densidad del fluido con la presión puede ser ignorada (en simulaciones de combustión a presión atmosférica, la variación de la densidad se da por parte de la temperatura y no por la presión). 4). No se considerará la distribución de calor en las paredes. 5). La composición de la mezcla gas-aire es homogénea. 6). La combustión es estequiométrica. Debido a su complejidad, el problema debe dividirse en 3 etapas diferenciadas: 2.5. Selección de software Actualmente, si bien se puede encontrar una gran variedad de códigos para la resolución numérica de problemas de dinámica de fluidos, la realidad es que este mercado se en- 19
20 cuentra monopolizado por el grupo norteamericano Fluent Inc. cuyo principal software es FLUENT. La competencia más importante de este grupo dentro del ámbito de los CFD es la del grupo ANSYS con el software ANSYS-CFX. En muchas ocasiones no se trata de dar los mejores resultados en una determinada aplicación para dominar un mercado sino de tener una buena estrategia de marketing y es quizás en esto donde sobresale este grupo: en primer lugar se encuentran presentes en todas las ramas en las que puede ser necesaria aplicación del CFD. En segundo lugar, va creando la necesidad del uso de dicha herramienta en el proceso de formación del ingeniero al ofrecer versiones educacionales en las universidades. Finalmente, y quizás es la gran ventaja que tienen sobre los pequeños grupos, ofrece un enorme asesoramiento técnico a sus clientes: asesora en la estrategia que se debe seguir para alcanzar los objetivos, dan una formación para el manejo del software. Para la realización del presente trabajo se optó por el programa Ansys-CFX v.11, debido a que para el uso de FLUENT se requiere de software adicional para la creación de geometría y mallado. Otra razón por la que Ansys-CFX v.11 fue seleccionado para el trabajo es porque es uno de los numerosos softwares de ANSYS, los cuales son compatibles entre sí. Con ANSYS-CFX se puede realizar el estudio de la interacción fluidos-estructuras. Además, este software ofrece todas las herramientas necesarias para un estudio: creación de la geometría, mallado, resolución y post-procesado. De esta manera, el usuario no tiene por qué enfrentarse a los problemas que ofrecen programas menos conocidos como puede ser el no reconocimiento de la geometría o el mallado a la hora de realizar el estudio en un CFD determinado o la imposibilidad de acoplar resultados del flujo y de la estructura. El programa proporciona rápidamente los resultados exactos del CFD. La versatilidad y capacidad de adecuación a las necesidades del usuario mediante la posibilidad de programación de funciones que permite la puesta en práctica de los nuevos modelos. Aun así, el código CFX pone a su disposición multitud de opciones entre las cuales están modelos para diferentes procesos físicos y químicos que le dan una enorme versatilidad. De esta manera, no sólo se puede realizar simulaciones de flujos laminares o turbulentos (incluye transición), newtonianos o no newtonianos, compresibles o incompresibles, monofásicos o multifásicos, sino también procesos de transferencia de calor por radiación, conducción y por supuesto por convección, así como procesos de fundición y con reacciones químicas, como combustión de gases, líquidos y combustibles sólidos. La disposición del solver y del post-processing interactivos del ANSYS -CFX hacen fácil detener un cálculo, examinar los resultados con el post-procesador, cambiar cualquier ajuste, y después continúar con el cálculo dentro de un solo uso. Los ficheros del caso y de datos pueden también ser leídos en el post-procesador de ANSYS para el análisis adicional con las herramientas avanzadas del post-processing y comparar resultados de diversos casos de lado a lado. 20
21 CAPÍTULO III PREPROCESO 3.1. Introducción En la etapa de preproceso se define la geometría del caso de estudio. El volumen de estudio se discretiza dividiéndolo en celdas de modo que el conjunto forma una malla (mesh). Se establecen también las condiciones de contorno, las condiciones iniciales y el modelo matemático a emplear Creación de una geometría tridimensional Para la simulación, a partir de los planos proporcionados por el fabricante se creó una geometría tridimensional mediante el programa Inventor v.11. La geometría mantiene las características relevantes para capturar el flujo pasando por alto detalles que a nivel de fabricación son esenciales pero desde el punto de vista de los procesos físicos que ocurren se consideran de nula importancia (por ejemplo tornillos, perillas, etc.). La geometría de la cocina (figura 3 a,b) está dividida en 3 dominios: quemador (figura 3 c), calientaplatos (figura 3 d), y horno (figura 3 e), considerando sólo el volumen ocupado por el fluido. La geometría simplificada carece de áreas pequeñas, salvo las de entrada y salida para evitar un número muy grande de los elementos en el proceso de mallado. Para el análisis se escogió la cocina a gas de uso doméstico con las siguientes dimensiones: Cámara de combustión: ancho 660 mm, profundidad mm, alto 220 mm, 7 orificios de entrada de aire secundario de area de mm2 cada uno y 2 ranuras laterales de 11,712.6 mm2 de area; 21
22 Horno: ancho 660 mm, profundidad mm, alto 370 mm, 2 orificios de 725,097 mm2 de area y 2 de 881,097 mm2; la camara de combustión y el horno tienen 16 orificios a los lados de 432 mm2 de área cada uno y un orificio central de mm2 de área. Quemador: 417 mm longitud de la cabeza, 187 mm. desde el inyector de gas hasta la cabeza, 116 agujeros de diámetro 2.38 mm, la distancia entre los agujeros es de 6 mm, el diámetro de la cabeza es de 22.2 mm; 660 mm 4xøx3212,388 mm mm 196,35 mm2 432 mm2 370 mm 10 mm 220 mm Figura 3 a. Cocina completa. Figura 3 b. Volumen ocupado por el fluido en diferentes zonas de la cocina. Figura 3 c. Quemador. Figura 3 d. Calientaplatos 22
23 Figura 3 e. Horno 3.3. Discretización espacial [15], [13] El objetivo de un código de CFD es la resolución de las ecuaciones de Navier- Stokes para un flujo con unas características determinadas y para una geometría concreta. El código de CFD resuelve dichas ecuaciones en una serie de puntos de la geometría en cuestión mediante métodos numéricos que transforman las ecuaciones diferenciales en ecuaciones algebraicas, proceso que se denomina discretización espacial. Por tanto, el primer paso que se da en la resolución de las ecuaciones es la selección de los puntos en los que se llevará a cabo la discretización, operación que se denomina generación del mallado. Para generar la malla, el dominio se divide en pequeños subdominios llamados elementos (figura 4). Control Volume* Figura 4. División del dominio en subdominios. 23
24 La geometría de éstos depende del número de dimensiones sobre las que estemos trabajando y de las características del problema: Problemas unidimensionales: segmentos de la línea que forma el dominio. Problemas bidimensionales: triángulos y cuadriláteros (figura 5). Los cuadriláteros se ajustan mejor a dominios cuadriláteros, mientras que los triángulos se suelen utilizar para dominios con una geometría más compleja. De hecho, cualquier dominio puede ser discretizado exclusivamente con triángulos, y no así con cuadriláteros. Problemas tridimensionales: tetraedros, piramides, hexaedros y prismas (figura 6). Los tetraedros son los elementos más simples y se ajustan bien a cualquier geometría, pero para superficies cilíndricas los hexaedros proporcionan mejor ajuste. Triangle Quadrilateral Figura 5. Elementos usados para mallar geometrías bidimensionales. Figura 6. Elementos usados para mallar geometrías tridimensionales. El mallado generado puede ser de dos tipos según los ejes coordenados que utilice: Mallado body-fitted, (figura 7), donde los puntos del mallado coinciden con la frontera física. Mallado en coordenadas cartesianas, (figura 8), donde las líneas del mallado son paralelas a los ejes cartesianos, independientemente de la forma de la frontera física. 24
25 Figura 7. Ejemplo de mallado body-fitted. Figura 8. Ejemplo de mallado en coordenadas cartesianas. La principal ventaja del mallado body-fitted es que permite soluciones muy exactas del flujo cerca de la frontera en el caso de geometrías muy complicadas a pesar de que requiere técnicas de construcción muy complejas. Por el contrario, el mallado en coordenadas cartesianas, si bien se construye muy fácilmente y tiene una mayor simplicidad en la resolución del flujo, es aplicable únicamente a geometrías muy sencillas ya que no permite un buen tratamiento de las condiciones de frontera. Por tanto, debido a las importantes limitaciones del mallado en coordenadas cartesianas, es el mallado body-fitted el más usado dentro del mundo de la industria, donde abundan geometrías y condiciones de contorno muy complejas. Existen muchos métodos para la discretización de las ecuaciones de Navier- Stokes. Las diferentes técnicas de discretización requieren siempre uno de los siguientes tipos de mallado: Mallado estructurado: Cada punto del mallado se identifica de manera exclusiva con los índices i, j, k del sistema de coordenadas cartesianas. Las celdas serán cuadriláteros en el caso de 2D y hexaedros en el caso de 3D. Mallado no-estructurado: Ni las celdas ni los nodos del mallado siguen un orden particular, es decir, las celdas vecinas no pueden identificarse por sus índices cartesianos (por ejemplo, la celda 6 estaría junto a la celda 128). Las celdas serán triángulos y cuadriláteros en el caso de 2D y tetraedros, prismas, hexaedros y pirámides en el caso de 3D. El tipo de discretización espacial a emplear, depende del tipo de discretización de las ecuaciones empleado, así como de la estructura interna de datos empleados para resolver el flujo. El método de las diferencias finitas, que discretiza las ecuaciones de Navier-Stokes en forma diferencial, requiere una malla estructurada de puntos en los que se guardan las variables de flujo. Por otro lado, en el método de los volúmenes finitos es necesario definir puntos de control en cada volumen generado. Precisamente es, en esos vértices donde se guardan las variables que, dependiendo del método que se 25
26 emplee pueden ser definidos en el centro del elemento (cel l-centred ) o en los vértices (cel l vertex ) de las celdas. (figura 9). Figura 9. Mallas con nodos centrados en los elementos y centrados en los vértices. Mallas estructuradas La principal ventaja de los mallados estructurados reside en la ordenación de los elementos en memoria, ya que de esta forma, el acceso a las celdas vecinas a una dada resulta muy rápido y fácil, sin más que sumar o restar un número al valor del índice correspondiente. Las mallas estructuradas pueden representarse en un sistema cartesiano o curvilíneo. En el primer caso, las líneas que configuran las celdas son siempre paralelas al sistema de ejes coordenados; por el contrario, en los sistemas curvilíneos, el sistema de coordenadas es deformado para adaptarse a la geometría del objeto de estudio. Por otro lado, también pueden clasificarse en mallas ortogonales y mallas no ortogonales. Las mallas ortogonales son aquellas en las que todas las líneas que la configuran se cortan entre sí con un ángulo de 90 o. Algunos flujos pueden considerarse como axisimétricos; en estos casos, las ecuaciones de flujo pueden ser expresadas en coordenadas polares (r, θ), en vez de en coordenadas cartesianas (x, y). Otra metodología, relacionada con los mallados estructurados, es la llamada técnica de la quimera. La idea básica es generar en primer lugar, mallas separadas alrededor de cada una de las entidades geométricas en el dominio. Después, las mallas se combinan de forma que se produzcan solapes de elementos. Lo más importante es la transferencia precisa de las cantidades entre las diferentes mallas y la región donde se produce el solape. La ventaja de esta metodología de mallado es que se pueden generar mallas para problemas particulares sin tener que generar mallas particulares, independientes unas de otras. Por otro lado, el problema de la técnica de la quimera es que la conservación de las propiedades de las ecuaciones de gobierno del fluido no son satisfechas exactamente en la zona donde se produzca solape de mallas. 26
27 Mallas no estructuradas El segundo tipo de mallas, son los mallados no estructurados. Éstos ofrecen gran flexibilidad en el tratamiento de geometrías complejas. La principal ventaja de los mallados no estructurados reside en que los triángulos (2-D) o los tetraedros (3-D), se pueden generar automáticamente, independientemente de la complejidad del dominio. En la práctica, es necesario determinar unos parámetros adecuadamente para obtener una buena calidad de malla. El tiempo requerido por una computadora para generar un mallado no estructurado es mucho menor que el que requiere para uno estructurado. Otra ventaja de este tipo de metodología es que la solución obtenida depende del refinamiento que se realice de la malla. Una desventaja de este tipo de mallado es el espacio que la malla ocupa en una computadora, es decir, requiere más memoria que los mallados estructurados. El mallado estructurado ofrece su mayor ventaja en el menor coste computacional que presenta dado que, a la posición de memoria en que se almacena la información correspondiente a un nodo, le sigue la información correspondiente al nodo adyacente, lo que permite una mayor rapidez en la realización del cálculo. No obstante este tipo de mallado no es muy sencillo adaptarlo a geometrías complejas. Sin embargo existe una solución que permite en geometrías complejas aprovechar la gran ventaja del mallado estructurado y que se denomina mallado en multibloques. Este mallado consiste en dividir la geometría en diferentes bloques, cada uno de ellos con una finura de mallado según las exigencias; así, los bloques mas próximos a la frontera tendrán mallados muy finos mientras que los bloques más alejados tendrán mallados menos finos. Sin embargo presenta la desventaja de requerir un código más complejo ya que será necesario pasar la información de los nodos de un bloque a los nodos de otro bloque diferente. A su vez, los mallados estructurados multibloques pueden ser conformes o no-conformes. Los mallados no-conformes son aquellos en los que los nodos de un bloque no tienen por qué coincidir con los nodos de otro bloque adyacente, por lo que el número de líneas de los bloques no tienen por qué ser iguales añadiendo mayor flexibilidad al mallado en multibloques. No obstante también será requerido un mayor tiempo en la elaboración del mallado. Respecto al mallado no estructurado, su ventaja principal es la posibilidad de la realización de un mallado automático por lo que es muy recomendable para casos de geometrías muy complejas. Además, el tiempo necesario para construir este tipo de mallados es considerablemente inferior al necesario para construir un mallado estructurado en multibloques. No obstante, dada la forma de almacenar la información en memoria, se requieren computadoras con mejores características que las que se podrían utilizar en un mallado estructurado o bien un mayor tiempo de cálculo. En el mundo de la industria, la elección de un mallado adecuado es fundamental a la hora de obtener una solución lo más exacta posible y sobre todo, en términos de ahorro de tiempo, ya que no sólo importa el tiempo empleado en la construcción del mallado 27
28 sino también el tiempo necesario para realizar el cálculo que, como hemos visto, depende de forma importante del tipo de mallado seleccionado. Actualmente, la mayor parte de los códigos emplean bien mallados no-estructurados o bien mallados híbridos que utilizan elementos rectangulares en 2D y hexaédricos en 3D y suponen una buena solución de compromiso entre exactitud, rapidez y requerimientos para la computadora. En el presente proyecto, los mallados empleados en todas partes de la cocina han sido mallados no-estructurados con elementos tetraédricos y, en las zonas de forma cilíndrica como en caso del quemador, con elementos piramidales y prismáticos. Se podría haber elegido también un mallado estructurado de un solo bloque pero la solución ofrecería poca exactitud en determinadas zonas del dominio como pueden ser las entradas / salidas de aire. Otra posibilidad sería un mallado estructurado multibloques, pero la geometría es relativamente compleja por lo que el mallado se realiza más rápidamente mediante el mallado no-estructurado. Por último, las características de las computadoras permitían desarrollar el cálculo en un tiempo dentro de lo habitual a pesar de que el mallado no- estructurado requiere tiempos de computación mayores que para un mallado estructurado. La complejidad de la física involucrada junto al tamaño del dominio definen a grandes rasgos el tamaño del problema y la potencia de cálculo necesaria. La densidad de nodos o elementos puede cambiar de unas regiones a otras debiendo acumular un mayor número de ellos en las zonas donde se esperan fuertes variaciones de alguna variable. En nuestro caso mallamos más intensamente la entrada de mezcla al quemador, salidas de la mezcla del quemador, entradas de aire secundario, fugas de aire, entradas de gases calientes al horno y salida de gases del horno (figuras 10,11 y 12), previamente asignando a estas zonas nombres, para la identificación de estas regiones en el momento de asignación de las condiciones de contorno. Por tratarse de una simulación en tres dimensiones, los elementos usados son tetraedros, pirámides y prismas. Cada celda de la malla será un punto de cálculo. Cuanto más pequeñas sean las celdas, de más valores por unidad de volumen dispondremos, con ello el resultado de la simulación se ajustará más a la realidad. Por otro lado, más ecuaciones a resolver tendrá el simulador CFX, con lo que mayor será el tiempo de cálculo al necesitar mayor consumo de memoria del ordenador. Así pues, se trata de generar una malla óptima, de forma que se ajuste a nuestra estructura y que pueda ser calculada por nuestro ordenador. Las herramientas que permiten un control directo sobre la distribución del tamaño de las celdas, pudiendo especificar los siguientes parámetros: 28 Expansión Factor-tasa de crecimiento de la malla Máximum Edge Length-valor máximo permitido de la malla Mínimum Edge Length-valor mínimo permitido de la malla
29 Los parámetros escogidos de todas las mallas aparecen en las tablas 1, 2 y 3 y los datos informativos de las mallas están en las tablas 4, 5 y 6. Tabla 1. Parámetros seleccionados para el mallado del quemador Tabla 2. Parámetros seleccionados para el mallado del calientaplatos 29
30 Tabla 3. Parámetros seleccionados para el mallado del horno. 10 a 10 b 10 c Figura 10 a,b,c. Quemador mallado 30
31 Tabla 4. Datos informativos de la malla del quemador 11 a 11 b 11 c Figura 11 a, b, c. Calientaplatos mallado 31
32 Tabla 5. Datos informativos de la malla del calientaplatos. a b Figura 12 a,b. Horno mallado Tabla 6. Datos informativos de la malla del horno. 32
33 3.4. Desarrollo del modelo matemático Introducción. Resolver el problema de la combustión de gas butano en un horno de la cocina de uso doméstico, implica la simulación de todos los procesos que tienen lugar en el hogar: 1. Modelo de turbulencia 2. Modelo de combustión 3. Modelo de radiación A continuación se detallan estos modelos, particularizando su aplicación al caso de estudio de este trabajo de tesis. El sistema de ecuaciones que describe el comportamiento del fluido es distinta en cada etapa. En la etapa 1 tenemos el flujo turbulento, isotérmico de la mezcla gas-aire; en esta etapa nos interesa la dirección y velocidad de salida del quemador de la mezcla que importaríamos como condiciones iniciales para la etapa 2. En la etapa 2 tenemos fluido turbulento con reacción química. Al resolver el sistema de ecuaciones correspondiente a este caso obtenemos el campo de temperatura y de velocidad del flujo de gases de combustión en las entradas al horno y también la distribución de temperatura en la bandeja que vamos a utilizar como condiciones iniciales en la etapa 3. En la etapa 3 tenemos el flujo de gases de combustión dentro del horno; en esta etapa determinaremos la distribución de calor dentro del horno y las direcciones de flujo de gases Ecuaciones de conservación [16] [ 8 ] El ingeniero francés Claude-Louis Navier (1822) y el matemático inglés George Stokes ( ) escribieron las ecuaciones básicas que describen el movimiento de un fluido, a las cuales se les conoce como ecuaciones de Navier-Stokes. Estas ecuaciones expresan en el lenguaje del medio continuo las tres leyes de conservación de la masa, ecuación de conservación del momentum y la ecuación de conservación de la energía. 33
34 La ecuación de continuidad (conservación de masa) se basa en la ley de conservación de la masa. Aplicando al concepto de movimiento de un fluido, la tasa de variación de la masa en un volumen de control es equivalente a la diferencia de la masa que entra y sale a través de sus fronteras. La tasa de variación de masa en el volumen fluido es: d! " x" y"z dt ( 3.1 ) El cambio del flujo másico que cruza la cara del elemento fluido que viene dada por el producto de la densidad, el área y la velocidad normal a la cara. En la figura 7 se puede analizar que el cambio neto del flujo másico que pasa a través de las fronteras del elemento fluido se escribe de la siguiente forma: Figura 13. Flujo másico que entra y sale de un elemento fluido. 34
35 ( 3.2 ) La tasa de variación de la masa que se encuentra dentro del elemento se iguala a la tasa neta del flujo másico que cruza las fronteras del mismo. Acomodando todos los términos del resultado del balance másico a la izquierda de la ecuación nos queda: En forma contraída: ( 3.3a ) ( 3.3b ) Asumiendo que el fluido es incompresible, la densidad r es constante y la ecuación se convierte en la siguiente: ( 3.4 ) la cual también puede escribir de forma reducido en notación vectorial: (3.5 ) 35
36 Ecuación de Conservación del momentum Se deriva de la aplicación del concepto de la segunda ley de Newton a un fluido en movimiento. La ecuación del momentum se expresa en términos de la presión y los esfuerzos debido a la viscosidad actuando sobre una partícula fluida. La tasa de variación del momentum en una partícula fluida es la diferencia de las fuerzas totales debido a los esfuerzos de la superficie y las fuerzas volumétricas que actúan sobre ella. La tasa de variación del momentum en cada dirección por unidad de volumen de u volumen fluido viene dado de la siguiente forma en cada dirección coordenada: Sobre un volumen fluido se distinguen dos tipos de fuerzas: Fuerzas de superficie: 36 - Fuerza de presión - Fuerzas viscosas Fuerzas volumétricas: - Fuerza gravitatoria - Fuerza centrífuga - Fuerza electromagnética - etc. ( 3.6 ) Es muy común colocar la contribución de las fuerzas superficiales en términos separados en la ecuación de momentum e incluir también los efectos de las fuerzas volumétricas como términos fuentes. El estado de esfuerzos sobre un volumen fluido esta definido en términos de la presión y las nueve componentes del tensor de esfuerzos mostrados en la figura 2. La presión, que es un esfuerzo normal, se denota como p. Los esfuerzos se denotan como τ. La notación en subíndices τ ij sirve para denotar la dirección del esfuerzo. El sufijo i y el j en τ ij indican que la componente del esfuerzo actúa en la dirección j en una superficie normal en la dirección i.
37 Figura 14. Esfuerzos superficiales sobre las caras de un elemento fluido. Primero se considera el componente x de las fuerzas debido a la presión p y los esfuerzos τ xx, τ yx y τ zx. La magnitud de la fuerza resultante de un esfuerzo de superficie es el producto del esfuerzo y del área. Las fuerzas que actúan en sentido positivo del ee un eje coordenado tienen signo positivo y si la fuerza va en sentido opuesto adquiere signo negativo. La fuerza neta en la dirección x es la suma de las componentes que actúan en esa dirección sobre el volumen fluido. La fuerza neta en las caras paralelas al plano yz en dirección x está dada por: La fuerza neta en las caras paralelas al plano xz en dirección x esta dada por: ( 3.7 ) ( 3.8 ) 37
38 y finalmente la fuerza neta en las caras paralelas al plano xy en dirección x esta dada por: ( 3.9 ) La fuerza total por unidad de volumen sobre el fluido debido a estos esfuerzos de superficie es igual a la suma de ( 3.7 ), ( 3.8 ), ( 3.9 ) y dividida por el volumen! x! y!z ( 3.10 ) Pero a esta ecuación se le deben incluir los efectos de las fuerzas volumétricas como termino fuente SMX de el componente x del momentum por unidad de volumen. El componente x de la ecuación del momentum es el cambio del momentum en el volumen fluido igualado a las fuerzas totales en la dirección x debido a los esfuerzos de superficie más la tasa de variación del momentum en el componente x debido a las fuentes: ( 3.11 ) De la misma forma se verifica que el componente y de la ecuación del momentum está dado por: 38 ( 3.12 ) De la misma forma se verifica que el componente z de la ecuación del momentum está dado por:
39 ( 3.13 ) Tensor de esfuerzos viscosos Los esfuerzos viscosos, originados de la interacción entre el fluido y la superficie del elemento, se representan mediante el tensor de tensiones!. En coordenadas cartesianas se escribe de la siguiente manera: ( 3.14 ) Las componentes τ xx, τ yy y τ zz,representan las tensiones normales y el resto de componentes que forman el tensor,!, representan los esfuerzos cortantes que actúan sobre el dominio fluido. En la figura 15 se representan todos los esfuerzos para un cubo infinitesimal. Figura 15 (a), (b). Esfuerzos normales y tangenciales actuando sobre un elemento fluido. 39
40 La mayoría de los fluidos presentes en problemas comunes, se pueden considerar como newtonianos. De esta forma, los componentes del tensor de esfuerzos viscosos quedan: que se escriben en forma contraída como: ( 3.15 ) ( 3.16 ) ( 3.17 ) en donde λ representa el segundo coeficiente de viscosidad, y µ el coeficiente de viscosidad dinámica, es la función delta Kronecker! ij = 1 si i! j y! ij = 0 si i = j! ij 40
41 Introduciendo las ecuaciones ( 3.15 ) en las ecuaciones ( 3.11 ), ( 3.12), y ( 3.13 ) se obtiene: ( 3.18 ) ( 3.19 ) ( 3.20 ) Las ecuaciones de Navier-Stokes se pueden escribir de la forma más útil para el desarrollo del método de volúmenes finitos: ( 3.21 ) 41
42 ( 3.22 ) ( 3.23 ) La ecuación de conservación de la energía se basa en la primera ley de la termodinámica, es decir, el cambio de energía de una partícula fluida es la diferencia entre el trabajo realizado sobre la partícula y el calor suministrado. La ecuación de la energía describe el transporte de la energía a través de un fluido y sus efectos. La tasa de incremento de la energía de una partícula por unidad de volumen esta dada por: ( 3.24 ) Trabajo realizado por las fuerzas superficiales La tasa de trabajo neto realizado sobre el volumen fluido por una fuerza de superficie es igual al producto de la fuerza y el componente de la velocidad en la dirección de la fuerza. Las fuerzas dadas por las ecuaciones ( 3.7 ), ( 3.8 ), ( 3.9 ) que actúan en la dirección x. El trabajo realizado por estas fuerzas esta dado por: 42
43 ( 3.25 ) La tasa de trabajo neto realizado por estas fuerzas superficiales actuando en la dirección x esta dado por: ( 3.26 ) Las componentes de los esfuerzos de superficie en la dirección y y z también realizan trabajo sobre la partícula y se puede comprobar que siguiendo el proceso anterior se llega a las siguientes ecuaciones respectivamente: ( 3.27 ) ( 3.28 ) La tasa total de trabajo realizado por unidad de volumen sobre la partícula fluida por todas las fuerzas superficiales esta dada por la suma de ( 3.26 ), ( 3.27 ), ( 3.28 ) dividido 43
44 por el volumen! x! y!z. Los términos que contienen la presión se pueden agrupar y escribirse en forma vectorial: ( 3.29 ) Esto conlleva la tasa total de trabajo realizado sobre la partícula por los esfuerzos superficiales. ( 3.30 ) Flujo de energía debido a la conducción de calor El vector de flujo de calor tiene tres componentes q x, q y y q z : Figura 16. Vectores de calor de un elemento fluido. La tasa de transferencia neta de calor hacia un volumen fluido debido al flujo de calor en la dirección x esta dada por la diferencia entre la tasa de calor impuesta sobre la cara oeste y la tasa de calor perdida en la cara este: 44
45 ( 3.31 ) De manera similar la tasa neta de transferencia de calor hacia el fluido debido al flujo de calor en la dirección y y z y ( 3.32 a, b ) La tasa total de calor adicionado al volumen fluido por unidad de volumen debido al flujo de calor a través de las fronteras es la suma de las ecuaciones ( 3.31 ), ( 3.32 a, b ) divididas por el volumen: ( 3.33 ) La ley de Fourier de la conducción de calor relaciona el flujo de calor con el gradiente local de temperatura de la siguiente forma: lo cual también se puede escribir en forma vectorial: ( 3.34 ) ( 3.34 a ) Combinando las ecuaciones ( 3.33 ) y ( 3.34 a) se puede escribir la forma final de la tasa de adición de calor sobre una partícula fluida debido a la conducción de calor: 45
46 ( 3.35 ) Se ha denotado la energía especificada en un fluido como E. La energía de un fluido esta definida como la suma de la energía (térmica) interna i, energía cinética: y energía potencial gravitacional. ( 3.36 ) Para la energía también se puede definir un término fuerte SE por unidad de volumen por unidad de tiempo. La ecuación de conservación de la energía en un volumen fluido establece la igualdad entre la tasa del cambio de la energía en el volumen y la suma de trabajo neto realizado sobre el volumen, la tasa neta de adición de calor y la tasa de variación de energía debido a las fuentes. La ecuación de energía se escribe entonces de la siguiente forma: ( 3.37 ) Es común extraer el cambio de la energía cinética mecánica para obtener una ecuación para la energía interna i o para la temperatura T. La parte de la ecuación de la energía que se atribuye a la energía cinética puede encontrarse multiplicando el componente x de la ecuación de momentum ( 3.11 ) por la componente u de la velocidad, la ecuación del momentum en y ( 3.12 ) por v y la ecuación de momentum en z ( 3.13 ) por w, luego se suman y se acoplan para obtener la ecuación: 46
47 ( 3.38 ) Sustrayendo la ecuación ( 3.37 ) de la ecuación ( 3.38 ) y definiendo un nuevo término fuerte como se llega a la ecuación de la energía interna: ( 3.39 ) Si utilizamos el modelo newtoniano para tensiones viscosas en la ecuación de la energía ( 3.39 ) se obtiene después de un reordenamiento: ( 3.40 ) Todos los efectos debidos a las tensiones viscosas en esta ecuación de la energía interna se describen por la disipación de Φ función que, después de un proceso algebraico considerable, se puede demostrar que es igual a: ( 3.41 ) 47
48 Para un fluido incompresible se tiene que i=ct donde c es el calor especifico a volumen constante y divu=0, esto convierte la ecuación anterior en una ecuación para la temperatura en la siguiente forma: ( 3.42 ) El movimiento de un fluido en tres dimensiones es descrito por un sistema de cinco ecuaciones en derivadas parciales (tabla 7): una de conservación de masa ( 3.3 a ), tres ecuaciones de cantidad de movimiento ( 3.21 ), ( 3.22 ), ( 3.23 ) y una ecuación de energía ( 3.40 ). Entre las incógnitas son siete variables: ρ, p, i y T, u, v, w. Para resolver el sistema la completamos con dos ecuaciones de estado. Podemos describir el estado de una sustancia en equilibrio termodinámico a través de dos variables de estado. Si usamos p y T como variables de estado las ecuaciones de estado para la presión p y para la energía son: Para gases perfectos: y y ( 3.43 ) ( 3.44 ) Los líquidos y los gases que fluyen a bajas velocidades se comportan como fluidos incompresibles. Sin variaciones de densidad no existe un vínculo entre la ecuación de la energía y la conservación de la masa y el momento.y el flujo a menudo se puede resolver teniendo en cuenta la conservación de masa y ecuaciones de momento solamente. La ecuación de la energía sólo hay que resolver junto con los demás si el problema consiste en la transferencia de calor. 48
49 Tabla 7. Ecuaciones de conservación. De las ecuaciones anteriores es claro que hay significativas semejanzas entre varias ecuaciones. Si introducimos una variable general φ, la forma conservativa de todas las ecuaciones de flujo de fluidos, incluyendo las ecuaciones para cantidades escalares tales como la temperatura, concentración de contaminante, entre otras, pueden ser escritas de la siguiente forma: ( 3.45 ) La conservación de una variable de flujo general φ, por ejemplo un componente de la velocidad o la entalpía dentro del volumen de control finito, puede ser expresada como un balance entre varios procesos tendientes a incrementarlo o disminuirlo: La Ecuación ( 3.45 ) es usada como el punto de partida para los procedimientos computacionales en el método de volúmenes finitos. Ajustando φ igual a 1, u, v, w y i ( ó T ó ho) y seleccionando los valores apropiados para el coeficiente de difusión Γ y los térmi- 49
50 nos fuente, obtenemos formas especiales de la tabla 7 para cada una de las ecuaciones parciales: masa, momento, y conservación de energía Modelos de turbulencia [ 15 ] El número de Reynolds Re de un flujo da una medida de la importancia relativa de las fuerzas de inercia, asociadas con los efectos convectivos, y las fuerzas viscosas. En experimentos con fluidos se observa que para valores inferiores a un número de Reynolds denominado crítico ( Re<2000 ), el flujo es intrínsecamente estable y las capas de fluido adyacentes se deslizan unas sobre otras de forma ordenada. El régimen del flujo se denomina laminar. Si el flujo tiene un valor del número de Reynolds por encima del denominado crítico (Re > 4000), se manifiestan en éste unas perturbaciones que dan lugar a un cambio radical del carácter del flujo. El movimiento se vuelve intrínsecamente no estacionario, incluso con condiciones de contorno constantes. Este régimen se denomina flujo turbulento. La turbulencia se define como el estado de movimiento de un fluido en el que las distintas variables relevantes (presión, velocidad, etc.) fluctúan de una forma desordenada. Se trata de un estado no estacionario desde el punto de vista macroscópico en el que las distintas variables adoptan valores dependientes tanto de la posición como del tiempo y estos valores varían de una forma aleatoria y desordenada. El flujo turbulento contiene zonas con vórtices de diferentes tamaños. Los vórtices de grandes escalas extraen energía del flujo y la transfieren a los vórtices de menores escalas, donde la energía es disipada en forma de calor debido a la viscosidad. La descripción del movimiento de las partículas fluidas debido al efecto de la turbulencia resulta altamente complejo y constituye un problema aún sin solución desde el punto de vista de los métodos numéricos. Se han propuesto varias formas de resolver el problema utilizando distintas aproximaciones. A continuación se exponen los métodos conocidos como simulación numérica directa, simulación de los grandes vórtices y promediado temporal de las ecuaciones de Navier-Stokes. Simulación directa de las ecuaciones ( Direct Numerical Simulation, DNS). Este método (cuyas iniciales provienen de la denominación inglesa Direct Numerical Simulation) consiste, en realidad, en no utilizar ningún modelo para la turbulencia, sino realizar discretizaciones temporales y espaciales que sean capaces de simular el flujo en un determinado problema. 50
51 El método DNS resuelve las ecuaciones completas de Navier-Stokes sin ningun modelo para el movimiento turbulento, ya que todas las escalas de la turbulencia son determinadas explícitamente. Este modelo tiene como limitación el gran tiempo de cómputo necesario para su resolución, lo que restringe su uso a estudios muy simples de carácter académico. Simulación de grandes vórtices ( Large Eddy Simulation, LES). La técnica LES resuelve las ecuaciones de Navier-Stokes después de aplicar un filtrado espacial, lo que es similar a realizar un promedio en el tiempo. Dependiendo de la anchura del filtro, que es generalmente una función del tamaño de la rejilla numérica empleada, LES puede capturar las fluctaciones más importantes en gran escala de las cantidades del flujo (explícitamente), dejando solamente las escalas relativamente pequeñas para modelarse empíricamente ( modelado utilizando un modelo de subrejilla). Cuanto más fino es el tamaño de la rejilla, mejor es la resolución de las escalas de turbulencia. Ecuaciones promediadas de Navier-Stokes, RANS La aproximación tipo RANS es conocida en ingles como Reynolds averaget Navier- Stokes, lo cual se traduce ecuaciones de Navier-Stokes promediadas temporalmente, han sido muy estudiados y resultan bastante útiles en la mayoría de los problemas prácticos resueltos mediante técnicas numéricas. El punto de partida es muy sencillo. Se trata de obtener una descomposición de las variables en su valor medio y su valor fluctuante. Por ejemplo, para la velocidad, la descomposición sería: ( 3.46 ) donde el componente medio de la velocidad se obtiene haciendo la integral de la velocidad instantánea: ( 3.47 ) 51
52 donde T es el tiempo promediado de la simulación, usualmente grande comparado con escala de tiempo de fluctuación turbulenta, u es el componente fluctuante. Sustituyendo ( 3.46 ) en la ecuación de momento, se obtienen las ecuaciones promediadas de Navier-Stokes, ( 3.48 ) donde ρ, p, y v representan la densidad, la presión y la viscosidad cinemática del fluido respectivamente. Las ecuaciones promediadas de Navier-Stokes presentan el término desconocido: el cual es la correlación entre las componentes de fluctuación de la velocidad y se le conoce como el término de los esfuerzos de Reynolds. La existencia de los esfuerzos de Reynolds implica que se deben hacer suposiciones sobre su estructura para dar solución al problema de turbulencia. Los modelos RANS han sido desarrollados en base al concepto de que la escala de velocidad y la escala de longitud son suficientes para describir el efecto de la turbulencia de flujo. En instancias iniciales el modelo de turbulencia utilizado era el modeo de viscosidad turbulenta el cual realiza una estimación de las escalas de velocidad y longitud del flujo por medio de los valores locales. Esto se lleva a cabo relacionando la viscosidad turbulenta con el gradiente de velocidad de flujo. Pero este simple modelo isentrópico no simula bien los flujos de mayor complejidad. Los modelos clásicos de dos ecuaciones como k-e son actualmente los más utilizados. Estos se basan en la presunción de que existe una analogía entre la acción de los esfuerzos viscosos y los esfuerzos de Reynolds en el flujo. Ambos esfuerzos aparecen en el lado derecho de la ecuación del momentum. En la ley de la viscosidad de Newton el esfuerzo viscoso es proporcional a la tasa de deformación de los elementos fluidos. Para un fluido incompresible se verifica 52
53 ( 3.49 ) Con el fin de simplificar la notación se utilizan los subíndices i y j. La convención de esta notación es que i o j=1, corresponde a la dirección x, i o j=2 a la dirección y, y la dirección z a i o j=3. Experimentalmente se observa que la turbulencia decae a menos que exista un esfuerzo cortante en el flujo incompresible isotermo. Por otra parte, se encuentra que el esfuerzo turbulento aumenta conforme aumenta la tasa de la deformación. Boussinesq en 1877 propuso que el esfuerzo de Reynolds se puede asociar a las tasas de deformación. Usando la nomenclatura de subíndices se escribe: ( 3.50 ) El transporte turbulento de calor, masa y otras propiedades escalares se estiman de igual forma. La ecuación anterior muestra que el momentum de transporte turbulento es proporcional a los gradientes de velocidad ( gradientes del momentum por unidad de masa). Análogamente el transporte turbulento del escalar se toma proporcional a los gradientes de la cantidad transportada, en notación de subíndices da como resultado: ( 3.51 ) donde Γt es la difusividad turbulenta. Cuando se mezclan los remolinos, se espera de alguna forma que los valores de la difusividad turbulenta Γt sean cercanos al valor de la viscosidad turbulenta μt, por esta razón se incluye el número de Prandtl/Schmidt: 53
54 como: ( 3.52 a ) ( 3.52 b ) ( 3.52 c ) ( 3.52 d ) ( 3.52 e ) Tabla 8. Ecuaciones promediadas de Navier-Stokes, RANS Al hacer el promediado en el tiempo se introducen 6 nuevas incógnitas (los esfuerzos de Reynolds: 3 normales y 3 de cizalladura). Similarmente, en las ecuaciones de transporte escalar aparecen nuevos términos extra (en círculos en las ecuaciones anteriores), conocidos como Fluxes de Reynolds o Fluxes Turbulentos. Ahora tenemos siete ecuaciones (una de continuidad, tres de momento, una de energía y dos de estado) y trece incógnitas: tres para la velocidad media, una para la presión media, una para densidad media, una para la energía media, una para la temperatura media y seis para el tensor de tensiones de Reynolds (3.53 ), ya que por ser este último simétrico solo seis de sus nueve componentes son distintos. El tensor de los esfuerzos de Reynolds: ( 3.53 ) 54
55 donde Entonces, la tarea principal del modelado de la turbulencia es encontrar procedimientos computacionales que permitan encontrar con suficiente exactitud y generalidad para los ingenieros los Esfuerzos de Reynolds y Fluxes de Reynolds. Para propósitos de ingeniería no es necesario resolver los detalles de las fluctuaciones turbulentas sino solo los efectos de la turbulencia sobre el flujo promedio, para ello tenemos los siguientes modelos clásicos: Modelo de una ecuación: modelo k. Modelo de dos ecuaciones: modelos k-ε, k-ω2 o q-ω, modelo RNG. Para este trabajo se seleccionó el modelo k-ε por su robustez y porque en comparación con otros modelos del tipo RANS, es el que menos tiempo de cálculo por iteración requiere. Modelo de turbulencia RANS k-ε Este modelo de turbulencia se centra básicamente en el mecanismo que afecta directamente la energía turbulenta. Algunos términos preliminares se deben tener en cuenta. La energía cinética instantánea k(t) de un flujo turbulento es la suma de la energía cinética promedio ( 3.54 ) y la energía cinética turbulenta por unidad de masa se define a partir de la velocidad fluctuante: ( 3.55 ) 55
56 ( 3.56 ) En las demostraciones realizadas a continuación se requiere utilizar la tasa de deformación y los esfuerzos turbulentos. Para facilitar la subsecuencia de la cálculos es muy común escribir los componentes de la tasa de deformación e ij en forma de tensor: ( 3.57 ) Descomponiendo la tasa de deformación de un elemento fluido en una componente promedio y fluctuante, resulta la siguiente matriz: 56
57 Ecuación para la energía cinética promedio K del flujo La ecuación para la energía cinética promedio puede obtenerse multiplicando el componente x de la ecuación del momentum promediadas temporalmente (4.11) por U, la componente y de la ecuación (4.12) por V y la componente z de la ecuación (4.13) por W. Luego de sumar los resultados y agruparlos y aplicando el álgebra se puede llegar a que la ecuación promediada de la energía cinética de flujo es de la siguiente forma: ( 3.58 ) En otros términos el significado de cada uno de los términos anteriormente escritos es el siguiente: Los términos de transporte (III), (IV) y (V) se caracterizan por la aparición del operador div. El efecto de los esfuerzos viscosos sobre K ha sido divididos en dos partes, el termino (IV), el transporte de K por los esfuerzos viscosos y el término (VI) es la disipación viscosa de la energía cinética K. Los dos términos que contienen los esfuerzos de Reynolds se deben a los efectos turbulentos: el término (V) es el transporte turbulento de K debido a los esfuerzos de Reynolds y el término (VII) es la producción turbulenta de K debido al trabajo de deformación debido a los esfuerzos de Reynolds. Para números de Reynolds altos los términos turbulentos (V) y (VII) son mucho mayores que las contribuciones viscosas (IV) y (VI). Ecuación para la energía cinética turbulenta k Multiplicando cada una de las ecuaciones de Navier-Stokes sin los esfuerzos de Reynolds por su componente de fluctuación de la velocidad (el componente x se multiplica 57
58 por u etc.) y sumando los resultados, seguido de una repetición de este proceso con las ecuaciones de Reynolds, restando las dos ecuaciones resultantes y re ordenando se llega a la siguiente ecuación para la energía cinética turbulenta k: ( 3.59 ) En otros términos el significado de cada uno de los términos anteriormente escrito es el siguiente: Las ecuaciones ( 3.58 ) y ( 3.59 ) parecen muy similares en muchos aspectos, a pesar que los términos cuantitativos del lado derecho de la ecuación muestran que el cambio de la energía cinética turbulenta viene determinado por las interacciones turbulentas. El término (VII) en ambas ecuaciones es igual en magnitud, pero opuesto en signo. El término (VII) proporciona una contribución positiva en la ecuación k y representa un término de producción. En la ecuación K representa una destrucción de la energía cinética. Esto expresa matemáticamente la transferencia de la energía cinética promedio a energía cinética turbulenta. El término (VI) es la disipación debido a la viscosidad 58 ( 3.60 ) y proporciona una contribución negativa a la ecuación ( 3.59 ) debido a la aparición de la suma de cuadrados de la tasa de deformación fluctuante eij. La disipación de la energía cinética turbulenta es debida al trabajo realizado por los vórtices de pequeñas escalas sobre los esfuerzos viscosos. La tasa de disipación por unidad de masa, cuyas dimensiones son m 2 /s 3 es de vital importancia para el estudio de dinámica turbulenta y se denota como:
59 ( 3.61 ) Éste es el término de destrucción de la energía cinética turbulenta, el cual es de un orden de magnitud similar al término de producción. En contraste, cuando se tiene un número de Reynolds alto, el término de transporte viscoso (IV) en la ecuación ( 3.59 ) es siempre mucho mas pequeño comparado con el término de disipación turbulenta (VI). Es posible desarrollar una ecuación de transporte similar para todas las otras variables turbulentas incluyendo la tasa de la disipación viscosa ε. La ecuación ε contiene términos desconocidos y que son difícilmente medibles. La ecuación estándar para el modelo k-ε, tiene dos ecuaciones modelo, una para k y otra para ε. Se usa k para definir el valor de la escala de la velocidad, v, y ε para el valor de la escala de longitud, l respectivamente, representando la mayor escala de la turbulencia, de la siguiente forma: La viscosidad turbulenta se escribe de la siguiente forma: ( 3.62 ) donde Cµ es una constante a dimensional. El modelo estándar k-ε, utiliza las siguientes ecuaciones de transporte: ( 3.63 ) ( 3.64 ) 59
60 En otros términos el significado de las ecuaciones es el siguiente: La constantes son obtenidas de condiciones experimentales, usualmente los valores de estas constantes tabl. son: Tabla 9. Coeficientes constantes. Modelo k-ε. A partir del modelo, los esfuerzos de Reynolds se calculan como: Donde δ ij es el tensor delta Kronecker: ( 3.66 ) " ij = 0 si i! j Las ecuaciones para k y ε se emplean para determinar la viscosidad turbulenta (fórmula 3.62), a partir de la cual se calculan los esfuerzos de Reynolds (fórmula 3.66), que posteriormente se reemplazan en las ecuaciones promediadas de Navier-Stokes, que junto con la ecuación de continuidad forman el sistema de ecuaciones que compone el método. Las principales ventajas e inconvenientes de este modelo son las siguientes: Ventajas: " = 1 si i = ij Sólo se necesita fijar las condiciones iniciales y de contorno. Resultados satisfactorios para una gran cantidad de flujos. j 60
61 Es el modelo turbulento más ampliamente utilizado en la mayoría de flujos en aplicaciones industriales. Se dispone de leyes de pared desarrolladas como condiciones de contorno para este tipo de modelos. Inconvenientes: Implementación más compleja que los modelos algebraicos debido a la introducción de dos ecuaciones diferenciales adicionales. Pobres resultados en casos como: flujos no confinados, flujos con grandes gradientes longitudinales, flujos turbulentos completamente desarrollados en conductos no circulares Las ecuaciones de conservación de especies [ 4 ], [1] Las ecuaciones de Navier-Stokes para los flujos con reacción química se diferencian de las ecuaciones para los flujos habituales en tres aspectos principales: El fluido reactante es una mezcla no isotérmica de múltiples especies (hidrocarburos, oxígeno, agua, dióxido de carbono, etc.) que deben ser seguidos de forma individual. Los datos termodinámicos también son más complejos que en la aerodinámica clásica, porque la capacidad calorífica del gas reaccionado cambia significativamente con la temperatura y composición. La velocidad de reacción química requiere de modelación específica. Puesto que el fluido es una mezcla de gases, los coeficientes de transporte (difusividad calorífica, difusión de especies, viscosidad, etc.) requieren una atención especial. La combustión consiste en las reacciones múltiples de especies químicas, y las ecuaciones de Navier-Stokes tienen que incluir algunos términos adicionales que consideran las propiedades de los componentes que constituyen la mezcla, que determinan las propiedades del fluido. En primer lugar las mezclas se caracterizan por la composición fraccional [1]. Sea una mezcla de i gases de los cuales, en un volumen de control V haya ni moles de cada uno. Definimos la fracción molar: 61
62 ( 3.67 ) Si las masas moleculares de los gases son Wi, definimos la fracción de masa ( 3.68 ) donde m i es la masa de especie i y m es la cantidad de masa total en el volumen de control. De las expresiones anteriores: Luego, como por definición ( 3.69 ) resulta: ( 3.70 ) 62 ( 3.71 ) Las expresiones anteriores nos permiten calcular la composición fraccional en peso (y i ) y en volumen (x i ) de la mezcla, así como el número de moles n y la masa molecular de la mezcla W. La presión parcial se relaciona con la densidad parcial de cada componente de la mezcla por:
63 ( 3.72 ) siendo luego, como y resulta: ( 3.73 ) Donde, T es la temperatura, R=8.314 J/(moleK) es la constante de gas perfecto Las variables primitivas de un flujo tridimensional reactante y compresible son: La densidad ρ = m / V. El campo de velocidades, en tres dimensiones ui Una variable de la energía (presión, entalpía o temperatura T). Las fracciones de masa Yk para N especies reactantes. La ecuación ( 3.73 ) es remplazada a la ecuación de conservación de masa y ésta se resuelve para cada componente del fluido. ( 3.74 ) donde el término fuente SI representa una fuente / sumidero, debido a las reacciones químicas. 63
64 En general, las reacciones químicas se pueden describir en términos de las reacciones elementales K relacionadas con componentes Nc que se puede escribir como: ( 3.75 ) donde vki es el coeficiente estequiométrico para el componente I de la reacción elemental k y Nc es el número de componentes. La tasa de producción/consumo, SI, para el componente I puede ser calculado como la suma de la tasa de progreso para todas las reacciones elementales en los que participa componente I: ( 3.76 ) donde Rk es la velocidad de reacción elemental de progreso para la reacción, que en ANSYS CFX se puede calcular mediante el modelo de Eddy Dissipation y/o el modelo Finite Rate Chemistry Modelo de combustión Simular el proceso de combustión demanda la inclusión de un mecanismo cinético que represente razonablemente el proceso que se quiere abordar numéricamente, a su vez, los mecanismos cinéticos se dividen en dos grandes grupos: Mecanismos Globales y Mecanismos Detallados. Un mecanismo cinético global, consigue reproducir informaciones generales del proceso de combustión, mientras que un mecanismo cinético detallado consigue reproducir particularidades del proceso de combustión. Esta diferencia es función del número de especies y reacciones químicas contempladas en los dos tipos de mecanismos. Una mayor cantidad de especies químicas contempladas en un mecanismo cinético va a permitir colectar una mayor cantidad de información que, dependiendo del proceso en análisis, pueden ser de gran impacto. No siempre es posible realizar una simulación en CFD usando un mecanismo cinético detallado de grandes proporciones. El problema radica en el número de especies químicas que han de tenerse en cuenta en la simulación. Por ejemplo, querer simular un proceso de combustión en CFD tridimensional, no estacionario, turbulento es imposible con los recursos computacionales actuales. Por esto es común en simulaciones CFD de combustión, usar mecanismos cinéticos globales o 64
65 reducidos, que sean posibles de implementar en los recursos computacionales actuales. Son diversos los métodos, con enfoque numérico, para abordar problemas de combustión. Diferentes programas CFD tienen inplementadas varias metodologías para abordar la solución numérica de un proceso de combustión. Estas metodologías envuelven el transporte de alguna variable y a partir de allí, surgen los diferentes modelos: a). Eddy Dissipation Model, adecuado para las reacciones rápidas; b) Finite Rate Chemistry Model, adecuado para reacciones más lentas Individuales y las reacciones de varios pasos c) Combinación de los métodos a y b; d) Laminar Flamelet Model, propicio para las llamas de difusión (la reacción sólo se produce en la superficie de la llama), e) NOx y Soot modelos. En este trabajo se utilizó el modelo de combustión Eddy Dissipation, usando el mecanismo de reacción de un paso (single one-step), que postula que la combustión se produce mediante la reacción irreversible: Combustible+Oxidante=Productos Donde los productos que se forman son CO 2 y H 2 O. Para el butano C 4 H 10 la reacción de combustión completa estequiométrica: C 4 H (O N 2 ) = 4CO 2 + 5H 2 O N 2 ( 3.77 ) El modelo de Eddy Dissipation se basa en el concepto de que la reacción química es rápida en relación con los procesos de transporte en el flujo. Cuando se mezclan los reactivos a nivel molecular, instantáneamente se forman productos. El modelo asume que la velocidad de reacción puede estar relacionada directamente con el tiempo necesario para la mezcla de reactivos a nivel molecular. En flujos turbulentos, la velocidad es proporcional al tiempo de mezcla y se define por la energía cinética turbulenta, k, y la disipación, ε. En el modelo de Eddy Disipation, el ritmo de avance de la reacción elemental k, está determinado por la menor de las dos expresiones siguientes: Limitador de reactivos: ( 3.78 ) 65
66 donde [I] es la concentración molar del componente I y sólo se incluyen los componentes del reactivo. Limitador de productos: ( 3.79 ) donde P iteraciones sobre todos los productos de los componentes en la reacción primaria k. El limitador de productos se desactiva cuando el modelo B coeficiente se establece en un valor negativo. Por tanto solo el paso y los esquemas multi-paso de la reacción, que está desactivada por defecto, (se establece en -1), pero puede ser activada por establecer explícitamente el coeficiente B modelo para un valor positivo (aunque esto no se recomienda para los sistemas de reacción de varios pasos). Al coeficiente A se asigna el valor 4. El modelo Eddy Disipation se utiliza para la combustión premezclada pero asignando el valor 0.5 a la constante B se activa la ecuación de limitador de productos y el modelo puede ser usado para la combustión parcialmente premezclada. En este caso es necesario seleccionar los valores iniciales diferentes de cero para los productos. El limitador de productos no se recomienda para las reacciones de varios pasos. [ 13a ] Modelo de radiación [4] La radiación es la energía emitida por la materia en la forma de ondas electromagnéticas (o fotones), como resultado de los cambios en la configuraciones electrónicas de los átomos o moléculas. A diferencia de la conducción y la convección, la transferencia de energía por radiación no requiere la presencia de un medio interpuesto. De hecho la transferencia de energía por radiación es la más rápida (a la velocidad de la luz) y no sufre la atenuación en el vacío. De esta manera llega a la Tierra la energía del Sol. En nuestro estudio es de interés la radiación térmica, que es la forma de radiación emitida por los cuerpos debido a su temperatura. Es diferente de las otras formas de radiación, como los rayos X, los rayos gamma, las microondas, las ondas de radio y de televisión, que no están relacionadas con la temperatura. Todos los cuerpos a una temperatura arriba del cero absoluto emiten radiación térmica. La radiación es muy importante en las cámaras de combustión, donde las temperaturas son elevadas porque según la ley de Stefan-Boltzmann, la cantidad de energía emitida por radiación depende de la diferencia de las temperaturas absolutas elevada a la potencia cuatro ( pag.28). 66 ( 3.80 )
67 donde A-área de la superficie, s - constante de Boltzmann, T s - temperatura absoluta! = 5.67 "10 #8 W / m 2 $ K 4 En el horno el medio está constituido por una mezcla de gases, algunos de ellos participativos, a alta temperatura. Los gases absorben y emiten radiación térmica dependiente de la temperatura y de la longitud de onda. Realizando un balance de la energía radiactiva a lo largo de una dirección s, el cambio de intensidad radiativa, se determina sumando las contribuciones correspondientes a la emisión, absorción, dispersión entrante en la dirección s y dispersión saliente en la dirección s, de forma que se obtiene una ecuación para la intensidad de la radiación. ( 3.81 ) El cambio de absorción y dispersión la intensidad = emisión - dispersión entrante + saliente en la radiactiva en la dirección s dirección s Donde: ka - coeficiente de absorción espectral, ks - coeficiente de dispersión, ε - emisividad espectral; I b λ - es la función de Planck, describe la intensidad de irradiación de cuerpo negro; Iλ - intensidad de la radiación espectral que depende de la posición ( r ) y dirección ( s ); f - función dispersión de fase; Ω - ángulo sólido. Por la ley de radiación de Kirchhoff, la emisividad espectral ε(λ) y el coeficiente de absorción espectral ka(λ) tomarán el mismo valor. Este valor se puede estimar aplicando un modelo de gases grises, descomponiendo el gas real, en uno claro y varios grises. Pero, para simplificar el problema y evaluar las propiedades del gas, se ha optado por un modelo de gas gris, con un valor de absortividad constante ( ε=1) para todo el espectro de frecuencias. De acuerdo con la ley de Planck, la intensidad de cuerpo negro I b, se obtiene integrando en todo el espectro la intensidad en cada longitud de onda λ (Ecs. y ). y ( 3.82 ) ( 3.83 ) 67
68 En cuanto a la dispersión, ésta se asume isótropa en todo el medio simplificándose así el último término de la ecuación de transporte de la intensidad de la radiación (Ec. 3.66) Es ecuación integro-diferencial de primer orden para Iλ en dirección s, su resolución es es muy costosa computacionalmente. Para resolverla existen muchos modelos de radiación que proporcionan las soluciones aproximadas. [13] Entre los modelos de radiación más usados están: Modelo Rosseland (Modelo de Aproximación de Difusión). Modelo P1 ( también conocido como el modelo Gibb s o Modelo de armónicos esféricos). Modelo de Transferencia Discreta (Modelo Shah). Modelo Monte Carlo. La elección del modelo de radiación térmica afectará no sólo la calidad de la solución, sino también el tiempo de cómputo requerido. La selección adecuada debe hacerse a partir de consideraciones físicas. [12]. La aproximación al medio ópticamente delgado u ópticamente grueso es la más utilizada en la teoría de transmisión por radiación. En el primer caso la probabilidad de absorción de fotones en el volumen estudiado es muy pequeña y, en el segundo caso, muy grande. Esto significa que en el medio ópticamente delgado los fotones de radiación propia fácilmente abandonan el volumen y los fotones de la radiación exterior lo atraviesan sin resistencias. En el caso de aproximación al medio ópticamente grueso, al revés, los fotones de radiación propia y de la radiación exterior son absorbidos con el propio medio. Es importante resaltar que en el proceso de resolución de los problemas de radiación con frecuencia los medios son aproximados a uno de los dos casos. Aproximación al medio ópticamente delgado. k! k! = e " K!S - espesor óptico K! ( 3.84 ) -coeficiente de absorción espectral; S-espesor del medio 68
69 k!!1 En el medio ópticamente delgado, esto permite simplificar considerablemente las ecuaciones de la transferencia por radiación. En aproximación al medio ópticamente grueso aproximación de Rosseland o de difusión. k!! 1 esta aproximación se llama Muchos de los flujos de fluidos de interés práctico se encuentran en situaciones donde el fluido y / o las paredes están calientes. En tales situaciones los efectos de la transferencia de calor radiante pueden llegar a ser importantes. Un ambiente típico donde la radiación juega un papel significativo es un horno o una cámara de combustión. Un extremo es la situación en la que el medio es transparente a la radiación en longitudes de onda en la que se produce la mayoría de la transferencia de calor. En este caso, la radiación sólo afecta al medio por el calentamiento o enfriamiento de las superficies del dominio, sin la transferencia de energía radiante directamente al medio. Para este caso se debe utilizar sólo el modelo Monte Carlo. El modelo de Transferencia Discreta también se ha utilizado en este caso, pero con éxito limitado. En el extremo opuesto está la situación en la que el medio es ópticamente denso y la radiación interactúa con el medio en el interior del dominio, así como en las superficies. Si el medio es ópticamente grueso, la energía radiante está bien dispersa o se absorbe y es reemitida en todas las direcciones con una escala de longitud pequeña en comparación con el tamaño del dominio. Esta situación se conoce como el límite de la difusión, ya que la intensidad de radiación es independiente de la dirección. En este caso, los modelos Rosseland y P1 son una alternativa mas atractiva que los modelos la Transferencia Discreta y Monte Carlo debido a su simplicidad. Para los casos generales, que van desde la óptica delgada (transparente) y ópticamente grueso (difusión), las regiones, como la combustión, los modelos de Transferencia Discreta y Monte Carlo representan la solución de la ecuación de transferencia radiactiva con mayor precisión. [13]. En los problemas de difusión o medio ópticamente grueso, todos los modelos producen casi los mismos resultados. Entonces la mejor alternativa son los modelos Rosseland y P1. A medida que disminuye el espesor óptico y se acerca a 1, el modelo P1 se convierte en la alternativa menos costosa. Y, por último, en el medio ópticamente delgado y transparente sólo los modelos Monte Carlo y Transferencia Discreta debe ser utilizados. Para los modelos con medios gris, el campo de radiación se espera que sea razonablemente homogéneo en todas partes (al menos a nivel local), y se requiere de alta resolución espacial, en este caso el método Transferencia Discreta es mucho más eficiente y proporciona resultados muy precisos si se utiliza suficiente resolución angular. Un problema importante con la transferencia discreta es la falta de información de error. 69
70 1. Modelo Rosseland 2.Modelo P1 70
71 3. Modelo Monte Carlo 71

References: RESOLUCIÓN 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución 
 resolución