SciELO - Scientific Electronic Library Online

 
vol.27 issue2Experimental determination of bagasse stiffness coefficientGrade-two resistance screw shear connector behaviour for a 28 MPa concrete section author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Ingeniería e Investigación

Print version ISSN 0120-5609

Ing. Investig. vol.27 no.2 Bogotá May/Aug. 2007

 

Modelo predictivo de combustión de mezclas de gas natural en motores de encendido por chispa (MECH)

A predictive model of natural gas mixture combustion in internal combustion engines

Henry Espinoza,1 Luís Patino,2 Yordy González3 y Lezama Irving4


1 Dr. Ingenierio mecánico. Centro de Termofluidodinámica y Mantenimiento (CTYM), Universidad de Oriente, Puerto La Cruz, Venezuela. hespinoza@cantv.net
2 M. Sc., Ingeniero mecánico. Centro de Termofluidodinámica y Mantenimiento (CTYM), Universidad de Oriente, Puerto La Cruz, Venezuela. luipaca1@doctor.upv.es
3 Especialista, Ingeniero mecánico. Centro de Termofluidodinámica y Mantenimiento (CTYM), Universidad de Oriente, Puerto La Cruz, Venezuela. yordygonzalez@cantv.net
4 Ingeniero mecánico. Centro de Termofluidodinámica y Mantenimiento (CTYM), Universidad de Oriente, Puerto La Cruz, Venezuela.


RESUMEN

En el presente trabajo se muestra el desarrollo de un modelo predictivo de la combustión de mezclas de gas natural para motores de encendido por chispa. El modelo se basa en la solución de dos zonas: una de mezcla combustible sin quemar y otra de productos de combustión. En cada una de ellas se resuelve para cada ángulo de giro del cigüeñal las ecuaciones de energía y conservación de la materia. Las ecuaciones diferenciales no lineales de la energía para cada una de las fases consideradas: compresión, combustión y expansión, se resuelven por el método de Runge-Kutta de cuarto orden. El modelo, además de permitir estudiar la combustión de diferentes composiciones de los componentes del gas natural, facilita evaluar la combustión en presencia de aire seco y húmedo. Se muestran los resultados de la validación con datos experimentales demostrando precisión y exactitud en los resultados que arroja el programa. Se presentan resultados de presión en el cilindro, temperatura en las zonas de mezcla sin quemar y quemada, fracción de masa quemada y calor de reacción de la combustión para el motor modelado usando una mezcla de gas natural.

Palabras clave: Modelo predictivo, combustión, gas natural, motores de encendido por chispa.


ABSTRACT

This study shows the development of a predictive natural gas mixture combustion model for conventional combustion (ignition) engines. The model was based on resolving two areas; one having unburned combustion mixture and another having combustion products. Energy and matter conservation equations were solved for each crankshaft turn angle for each area. Nonlinear differential equations for each phase’s energy (considering compression, combustion and expansion) were solved by applying the fourth-order Runge-Kutta method. The model also enabled studying different natural gas components’ composition and evaluating combustion in the presence of dry and humid air. Validation results are shown with experimental data, demonstrating the software’s precision and accuracy in the results so produced. The results showed cylinder pressure, unburned and burned mixture temperature, burned mass fraction and combustion reaction heat for the engine being modelled using a natural gas mixture.

Keywords: predictive model, combustion, natural gas, internal combustion engine.


Recibido: diciembre 17 de 2006
Aceptado: junio 28 de 2007

Introducción

Debido a las grandes reservas mundiales de gas natural se está realizando cada vez más investigación y desarrollo para el uso de este combustible alterno en motores de encendido por chispa y de compresión. En países como Venezuela, con reservas importantes de gas natural, existen programas nacionales destinados al uso masivo de este combustible con el nombre genérico de gas natural para vehículos GNV (PDVSA, 1997). A pesar que el uso de este combustible no es nuevo, aún se requieren mayores estudios para conocer mejor el comportamiento del motor con él, en cuanto a la combustión y emisiones para diferentes mezclas, desgaste, confiabilidad, lubricación y diseño de las partes del motor, así como en los sistemas de alimentación. En el presente trabajo se desarrolla un modelo predictivo que permite estudiar el comportamiento de las variables características de la combustión del GNV para MECH, como presión y gradiente de presión en el cilindro, presión máxima, temperatura en el cilindro, velocidad de combustión, fracción de calor liberado y productos de la combustión en función de los parámetros operativos, composición del gas natural y parámetros de diseño.

Descripción del modelo

Balance de energía en la cámara de combustión

Para obtener las variables presión, temperatura, volumen, masa y flujo de calor, se efectuaron balances en los tres procesos que se presentan durante la combustión.

Etapa de compresión: aplicando la primera ley de la termodinámica en función del ángulo del cigüeñal y la ecuación de estado para la mezcla de gases ideales, se obtuvo un sistema de ecuaciones que contiene las variables de interés en esta etapa, presión y temperatura. Las ecuaciones de balance de energía en términos de presión y temperatura para la mezcla combustible sin quemar son:

Donde, dp/dα es el gradiente de presión dentro del cilindro, Rm la constante de de los gases ideales para la mezcla, V el volumen interior dentro del cilindro para el ángulo α, Cvm el calor específico a volumen constante de la mezcla, dQm/dα la de transferencia de calor instantánea, p la presión en el cilindro, dV/dα la variación del volumen con el ángulo de giro de cigüeñal, y Tm la temperatura de la mezcla y α el ángulo del cigüeñal. El subíndice m denota la mezcla sin quemar.

La variación del volumen con respecto al ángulo de giro del cigüeñal se calcula por la ecuación de cinemática biela-manivela y el flujo de calor del gas a la pared se obtuvo teniendo en cuenta la transferencia hacia las paredes del cilindro, el pistón y la culata de la cámara de combustión, calculándose el coeficiente convectivo de transferencia de calor por la ecuación de Woschni (Günter, et al., 2002).

Donde D es el diámetro del pistón [m]; p la presión instantánea del gas [New/m2]; la velocidad media del pistón [m/seg]; po la presión instantánea con motor arrastrado o correspondiente al ciclo del motor [New/m2] Vs el volumen recorrido por el cilindro [m3]; T1, p1, V1 la temperatura, presión y volumen de referencia o condiciones iniciales [K, New/m2, m3].

Las constantes C1 y C2 varían con la fase de trabajo de la siguiente manera:

- C1 = 2,28 para compresión, combustión y expansión [m/segK].

- C2 = 0 para compresión.

- C2 = 3,24x10-3 para combustión y expansión [m/segK].

Este coeficiente convectivo de transferencia de calor se asume constante para todas las superficies del cilindro. A medida que el proceso de compresión se avanza, las variables son recalculadas usando la expresión general:

Donde M es cualquier variable y δMn el incremento de la variable durante un paso de tiempo; se calcula utilizando el algoritmo de Runge-Kutta de cuarto orden.

Período de espera (Δα): se asumió que el núcleo de gases quemados en el período de espera es de forma semiesférica. Para calcularlo se recurrió alas siguientes fórmulas:

Donde N son las rpm, rt el radio del frente de llama, (Vt) la velocidad de combustión turbulenta, Vt la velocidad de llama turbulenta y ff un factor de proporcionalidad entre las velocidades de combustión laminar y turbulenta.

Para calcular la rapidez de la llama laminar asociada con el núcleo de inicio de la combustión, se usó la expresión de Ryan T. W. y Lestz (Sierens, 1998).

La velocidad de llama turbulenta depende del factor ff. Este se obtiene por ensayo y error, comparando los resultados P-α obtenidos con la aplicación del presente modelo con resultados experimentales para el motor simulado. El valor de ff que mejor ajuste los resultados es el que se toma como parámetro de simulación.

Una vez que se conoce el Δα por la ecuación (5), se continuó la compresión durante este período de tiempo. En este momento la combustión se asume que ha comenzado y se actualizan las variables p, T, m y V para cada zona por las siguientes ecuaciones (la temperatura inicial de los productos, Tp, se obtiene por la siguiente ecuación).

Esta fórmula se logra por la aplicación de la ecuación de la energía para un sistema suponiendo calor específico constante, que establece que la energía de los productos es igual a la inicial de la mezcla más la aportada por la combustión (tercer término de las ecuaciones 9 y 10). Donde fi es la relación de la mezcla aire-combustible de toda la masa dentro del cilindro. Considerando la formación del núcleo de llama se calcula el volumen de la mezcla y de los productos con el fin de estimar el radio del núcleo. El radio del núcleo de llama se obtiene por la ecuación (11).

Etapa de propagación y expansión con dos zonas. El efecto de la propagación de la llama y los cambios en las variables termodinámicas se halló usando el método numérico de integración Runge-Kutta de cuarto orden. Las ecuaciones para estimar la presión y la temperatura de la mezcla y de los productos en esta etapa son obtenidas del balance de energía (ecuaciones 1 y 2):

Donde p, T y V son la presión, temperatura y volumen dentro del cilindro, U la energía interna, R la constante de los gases, la masa, Cp el calor específico a presión constante, Cv el calor específico a volumen constante y α el ángulo de giro del cigüeñal; m y p denotan mezcla combustible sin quemar y productos de la combustión, respectivamente.

El término dmp/dα se determina por la ecuación de continuidad,

donde ρm es la densidad de la mezcla sin quemar [Kg/m3], Vt la velocidad de llama turbulenta [m/s], As el área frontal de la llama entre la zona quemada y sin quemar [m2] y ω la velocidad angular [rad/seg]. El área del frente de llama se calcula suponiéndola esférica y determinando para cada instante su posición mediante la velocidad turbulenta de la llama.

La transferencia de calor para cada zona se calcula usando la ecuación (16).

Las variables p, Tp y Tm son actualizadas utilizando la ecuación (4) y el algoritmo Runge-Kutta de cuarto orden. Se asume que la combustión se ha completado cuando el volumen de la mezcla sin quemar es igual a cero.

Expansión con los productos de la combustión. Para esta etapa se emplearon las ecuaciones (1) y (2) con el fin de modelar el proceso y calcular p y T. Las ecuaciones se transforman en:

El procedimiento para resolver estas ecuaciones es el mismo utilizado en la parte (a), con la diferencia de que el fluido de trabajo es el producto en fase gaseosa (Moreno et al., 1994; Horlock y Winterbone, 1986).

Balance de masa en la cámara de combustión

La Figura 1 muestra las corrientes de entrada y salida de la cámara de combustión en el motor ECH.

Partiendo del balance de la ecuación (19), en el que aparecen 14 incógnitas (na, nf y los ni de los productos) fue necesario encontrar para su resolución otras ecuaciones, razón por la cual se aplicó la conservación de la masa a cada elemento con el fin de proporcionar cinco ecuaciones y establecer un sistema de ecuaciones, dado por:

Balance de carbono:

donde “c” son los átomos de carbono y “b” los de hidrógeno, y el cual se realiza para cada una de las fracciones de hidrocarburo presentes en el combustible (c = 1, 2,3 ,4 y 5 y b = 4, 6, 8, 10 y 12, desde metano hasta pentano).

Balance de hidrógeno:

Balance de oxígeno:

Balance de nitrógeno:

Balance de argón:

Además de estas cinco ecuaciones se pueden añadir ciertas hipótesis dependiendo de las condiciones en que tiene lugar la reacción. La considerada en este trabajo fue la de una reacción de combustión en equilibrio con disociación, para la cual se necesitan siete ecuaciones de disociación. La condición de concentración se utiliza como ecuación auxiliar:

Ecuaciones auxiliares para resolver el balance de masa

Mediante el análisis CHON de doce especies se plantean las reacciones de disociación de los reactivos y productos en la cámara de combustión (Araque, 1997).

Con el fin de estimar la constante de equilibrio en términos de las fracciones molares se aplicó la ecuación (33) a cada una de las reacciones anteriores y se hicieron las transformaciones correspondientes.

Donde, a y b son los coeficientes de los reactantes de la ecuación y c y d los estequiométricos de los productos. Dados los valores de p, T, φ, nf y r se estiman los yi. Inicialmente se evalúan las constantes de equilibrio Kp en función de la temperatura y la presión.

Métodos numéricos empleados en el modelo

El método numérico usado para hallar las variables de funcionamiento fue el Runge-Kutta de cuarto orden para sistemas diferenciales de primer orden, en donde el incremento del ángulo de rotación δα es el parámetro que controla la exactitud de la solución y el tiempo de cálculo. Para hallar las fracciones molares de los productos nos valimos del método de Newton-Raphson para una sola variable no lineal y para sistemas de varias variables no lineales con una tolerancia de 0,000001 (Burden y Faires, 1985).

Desarrollo del software y validación del modelo

Teniendo como base los conceptos teóricos y las ecuaciones expresadas se elaboró un software, denominado Duspemec, en lenguaje Visual Basic versión 4.0, para realizar los cálculos con el modelo (Lezama, 1999). En vista de que sólo se pudo obtener la curva experimental de combustión del propano con aire en base húmeda, se utilizaron estos datos como patrón para validar los resultados del programa. Actualmente se llevan a cabo los estudios experimentales para la validación con mezclas de gas natural.

Los datos experimentales se obtuvieron en motor Toyota Starlet 2E, de cuatro cilindros en línea, con 1.296 cc de cilindrada. Para establecer los resultados de presión en función del ángulo de giro del cigüeñal se utilizó un sensor de presión tipo 6053 Bsp, instalado en la cámara de combustión, que está conectado por medio de un cable a un amplificador de carga tipo 5037B1, a la unidad de control y adquisición de datos HP3852A. Esta unidad es a su vez conectada a una computadora que tiene acoplada una tarjeta procesadora de lenguaje HP8232, la cual actúa como interfase entre la unidad de control y adquisición de datos y la computadora. El banco de prueba cuenta con un sistema medidor del ángulo de giro del cigüeñal llamado CAM, que es un codificador óptico que va acoplado a la unidad de control y de adquisición de datos para sincronizar los giros del cigüeñal con las presiones dentro del cilindro.

En la Figura 3 se representan los valores de presión del sistema obtenidos por el software (línea continua) y los experimentales en función del ángulo de giro del cigüeñal (puntos), con un factor de llama (ff) de 4. Se observa una clara similitud en ambas curvas.

Para demostrar cuantitativamente la reproducibilidad de los resultados del software se determinó el modelo de regresión lineal simple que mejor ajusta los datos de presión calculados en función de los medidos experimentalmente, obteniéndose también el análisis de varianza del modelo. La curva de ajuste sigue un modelo lineal, pcalc = [0,146855 + 0,9915 pexp] x105 N/m2, con un coeficiente de correlación del 99,94% y un nivel de confianza del 95%.

Para estimar el grado de significancia de las pendientes de la curva experimental y calculada (P vs α) se aplicó la prueba de hipótesis relativa, que utiliza la estadística de prueba “t de Student” (Montgomery, 1991). Se determinó que no hay diferencias significativas entre los datos computados y los experimentales.

Con la finalidad de obtener un modelo de ajuste para la curva de presión frente al ángulo del cigüeñal, se aplicó análisis de regresión lineal múltiple, encontrando la ecuación polinomial siguiente, con un coeficiente de correlación del 98,2% y un 95% de confianza. Las unidades de α y p son en grados y N/m2, respectivamente.

Predicción de los parámetros de operación

Se corrió el programa con los datos del GNV dados en la Tabla 1, lográndose los siguientes resultados.

La Figura 4 indica la evolución de la presión en función del ángulo de giro del cigüeñal con un factor de llama de 5,35, y la presión máxima alcanzada fue de 46,873x105 N/m2.

La temperatura de los productos con respecto al ángulo del cigüeñal se muestra en la Figura 5, siendo sus máximos de 2.678,41 °K y 1.260,08 °K, respectivamente.

La fracción de masa quemada dentro del cilindro en función del ángulo de giro del cigüeñal se señala en la Figura 6. El comportamiento observado representa fielmente el encontrado en los motores de encendido por chispa, según las referencias bibliográficas (Günter et al., 2002).

Predicción del calor de reacción para el balance de energía. En la Figura 7 se representa la evolución del calor de reacción (ΔHr) a diferentes temperaturas de los productos.

Esta curva tiene un comportamiento parecido al de la fracción de masa quemada (Figura. 6) debido a que la cantidad de reactivos es la variable que determina el calor que libera la reacción en la transformación de la mezcla de reactivos a productos durante el proceso de combustión (Payri y Tinaut, 1997).

Conclusiones

El modelo desarrollado predice los parámetros de combustión, como son: presión, temperatura, fracción de mezcla sin quemar y quemada, durante la combustión de un motor de combustión interna encendido por chispa usando mezclas de gases como combustible, por lo cual puede ser destinado a estudiar el funcionamiento del motores que utilicen gas natural como combustible.

El modelo requiere de un ajuste inicial con datos experimentales del motor a simular, en particular para obtener el factor de llama más adecuado y establecer el cálculo de la velocidad de combustión turbulenta.

Agradecimiento

El presente estudio fue realizado gracias al apoyo parcial del Conicit, a través de cofinanciamiento del Laboratorio de Usos Múltiples de Ingeniería Mecánica, N° (F-96000645).

Bibliografía

Araque, J., Modelo de Combustión para un Sistema CHON de 18 Especies., Memorias del Segundo Congreso Nacional de Ingeniería Mecánica, Merida, Universidad de los Andes, 1997.        [ Links ]

Burden R. y Faires J., Análisis Numérico., Editorial Iberoamericana S. A., México, 1985.        [ Links ]

Günter P., Schwarz C., Stiesch G. and Otto F., Simulation of Combustion and Pollutant Formation for Engine- Development., Editorial Springer, España, 2002.        [ Links ]

Horlock, J. and Winterbone D., The Thermodynamics and Gas Dynamics of Internal Combustion Engines., Clarendon Press, Volumen II, Oxford, 1986.        [ Links ]

Lezama, I., Desarrollo de un software para la predicción de las emisiones en un motor de encendido por chispa que utiliza GNV., Tesis presentada en la Universidad de Oriente, Puerto la Cruz, Venezuela, para optar al grado de Ingeniero Químico, 1999.        [ Links ]

Montgomery, D., Diseño y Análisis de Experimentos., Grupo Editorial Iberoamericano, México, 1991.        [ Links ]

Moreno J., Espinoza H., y Baduy F., Modelo de Combustión de Dos Zonas para Motor Encendido por Chispa., Memorias del V Congreso Latinoamericano de Transferencia de Calor y Materia, Caracas, Universidad Simón Bolívar, 1994, pp. IVB-2.2.        [ Links ]

Payri, F., y Tinaut F., Modelo de Cálculo de Emisiones Contaminantes para Motores de Combustión Interna Alternativos., Memorias del VIII Congreso Nacional de Ingeniería Mecánica, Anales de Ingeniería Mecánica, Barcelona, Vol. 2, 1997, pp. 193-199.        [ Links ]

PDVSA., GNV: Gas Natural para Vehículos., Boletín Informativo a la Gerencia General de Mercadeo Nacional, Caracas, 1997.        [ Links ]

Sierens, R., Combustion Model for Rotary Engines., Publicaciones del Congreso Modelado en Motores de Combustión Interna, Universidad Politécnica de Valencia, España, 1998.        [ Links ]

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License