Motores Termicos
Motores Termicos
Motores Termicos
TOMO I
JOHN R. AGUDELO
UNIVERSIDAD DE ANTIOQUIA
MOTORES TÉRMICOS
TOMO I
Principios básicos de Motores de Combustión Interna Alternativos
Universidad de Antioquia
Facultad de Ingeniería
Departamento de Mecánica
Medellín, 2000
A Diana y a Isabel
Agradecimientos
Quiero expresar mis agradecimientos sinceros a D. Manuel Muñoz Torralbo, hasta 1998
Catedrático de la unidad docente de Motores Térmicos de la Universidad Politécnica de Madrid,
y a su colega y director de mi Tesis Doctoral D. Jesús Casanova Kindelan por su aportación en
mi formación y en la discusión de parámetros claves para la ejecución de este proyecto.
Reconocimientos a mi hermano Andrés Felipe, quien actualmente realiza sus estudio de
Ingeniería Mecánica en la Universidad Nacional de Colombia sede Medellín, por su colaboración
en la transcripción de este trabajo y en la edición de las figuras.
A la Universidad de Antioquia por brindarme el espacio para transmitir los conocimientos
adquiridos durante mi formación y por ser promotora en la difusión del texto.
Finalmente a todos mis amigos y colegas que de una u otra forma aportaron para que el
sueño se convirtiera en realidad.
i
INTRODUCCIÓN
Este primer tomo ha sido concebido de manera integral, como material de apoyo y consulta para
estudiantes de ingeniería que se encuentren en su proceso de formación con el tema de los
motores de combustión interna alternativos.
El libro no es de mi completa autoría. La mayor parte es una recopilación de diferentes
obras de reconocidos autores en el campo de los motores de combustión interna alternativos (en
adelante MCIA). Principalmente he seguido el libro de John B. Heywood titulado “Internal
Combustion Engine Fundamentals”, el cual fue publicado por McGraw – Hill en Nueva York el
año de 1988 y que hasta el momento no ha sido traducido al español. He incluido muchos datos
nuevos, fruto de mis investigaciones y de trabajos recientes de otros investigadores.
En el capítulo 1 se hace una breve historia de los motores de combustión en general, con
el fin de introducir al estudiante en el mundo de los motores térmicos. Para recorrer los capítulos
posteriores del libro, el estudiante necesitará los fundamentos estudiados previamente en
Mecánica de Fluidos, Termodinámica y Transferencia de masa y de calor. Además se revisarán a
lo largo del estudio de los motores algunos criterios de Diseños Térmicos, Diseño Mecánico,
Cinemática y Dinámica, ya que todas estas áreas del conocimiento se aplican a profundidad en
las Máquinas Térmicas y muy especialmente en los motores de combustión interna alternativos
que son el pilar de este primer tomo.
En los capítulos 2 a 4, se establecen los principios geométricos y de funcionamiento de
los MCIA, se hace una revisión de la termoquímica de las mezclas aire – combustible y se
establece la forma como se trabajará con las propiedades de los fluidos de trabajo. Se ha
decidido adoptar para este último caso las técnicas computacionales basadas en las Tablas de
JANAF, sin embargo se mencionan las cartas termodinámicas.
En el capítulo 5 se estudian los ciclos de trabajo teóricos y se comparan, hacia el final del
capítulo, con los ciclos de trabajo reales, haciendo especial énfasis en las irreversibilidades
propias del proceso de combustión. Los cálculos termodinámicos de los ciclos de trabajo teóricos
nos permiten conocer el máximo trabajo por unidad de combustible que en situaciones ideales se
pudiera alcanzar. Sin embargo dichos cálculos no nos permiten definir la máxima potencia para
un motor determinado, pues depende del número de revoluciones por minuto permisibles, de las
características de los materiales empleados, de ciertas propiedades de los combustibles y
lubricantes usados, etc.; pero sí permiten definir un máximo valor del rendimiento en la
conversión energética, a la cual debe tender dicho motor. En estos cálculos generalmente se toma
el poder calorífico másico inferior del combustible, determinado a presión constante, como el
ii
máximo trabajo obtenible de él. Esto solo sería posible, si la totalidad del proceso, incluyendo la
combustión, se realizara reversiblemente. No es así, ya que la realidad impone sus
irreversibilidades en procesos como: estrangulación, rozamiento, turbulencia, transmisión de
calor a las paredes, entre otros. Además los fluidos que realizan trabajo deberían regresar
reversiblemente a las condiciones que tenían antes de iniciar el ciclo, lo que tampoco ocurre en
los procesos reales.
En los capítulos siguientes se revisan los principios de transferencia de calor y
refrigeración, lubricación y desgaste, tipos de combustibles y sus sistemas de suministro al
interior del cilindro y, de manera muy detallada se dedican los tres capítulos finales al estudio de
la combustión y la formación y determinación experimental de las emisiones contaminantes en
los dos tipos de motores objeto de este primer tomo: los Motores de Encendido Provocado (en
adelante MEP) y los Motores de Encendido por Compresión (en adelante MEC).
A lo largo del libro se citan las referencias bibliográficas consultadas y en la mayoría de
los casos, se ha agregado al final, un apartado con las referencias bibliográficas recomendadas
para ampliar información.
iii
CONTENIDO
Capítulo 1
Breve historia de los motores de combustión
Introducción................................................................................................................................. 1
Las máquinas de fuego ................................................................................................................ 2
Vacío a partir de la condensación de vapor................................................................................. 3
Motores de aire caliente............................................................................................................... 8
Del Calórico a la termodinámica ................................................................................................. 9
Evolución histórica de los motores de combustión interna alternativos ................................... 10
El motor de Lenoir................................................................................................................. 12
Principios de Beau de Rochas............................................................................................... 13
El motor de Brayton .............................................................................................................. 14
El motor Otto......................................................................................................................... 14
El motor de dos tiempos ........................................................................................................ 16
El motor Diesel...................................................................................................................... 16
El motor Wankel .................................................................................................................... 19
Literatura Recomendada............................................................................................................ 22
Referencias ................................................................................................................................ 22
Capítulo 2
Tipos de motores y su funcionamiento
La máquina de fluido................................................................................................................. 23
El motor térmico........................................................................................................................ 24
Motores de combustión externa e interna.................................................................................. 25
Clasificación de los MCIA ........................................................................................................ 25
Ciclos de operación de los motores ........................................................................................... 26
Componentes del motor............................................................................................................. 29
Operación del MEP ................................................................................................................... 31
Operación del MEC ................................................................................................................... 34
Motores de carga estratificada................................................................................................... 38
Literatura recomendada ............................................................................................................. 40
Referencias ................................................................................................................................ 40
Capítulo 3
Parámetros de operación y diseño de motores
Características importantes de los motores ............................... ¡Error! Marcador no definido.
Propiedades geométricas de los motores alternativos ............................................................... 42
Par y Potencia efectivos............................................................................................................. 44
Trabajo indicado por ciclo......................................................................................................... 45
iv
Rendimiento Mecánico.............................................................................................................. 47
Presión media efectiva............................................................................................................... 48
Consumo especifico de combustible y Rendimiento................................................................. 49
Relaciones aire - combustible y combustible - aire................................................................... 51
Rendimiento Volumétrico ......................................................................................................... 51
Peso y volumen específicos del motor ...................................................................................... 52
Factores de corrección para la potencia y el rendimiento volumétrico ..................................... 53
Emisiones específicas e índice de emisiones............................................................................. 54
Relaciones entre los parámetros de funcionamiento ................................................................. 55
Diseño de motores y datos de funcionamiento.......................................................................... 57
Referencias ................................................................................................................................ 62
Capítulo 4
Termoquímica de las mezclas aire - combustible y propiedades de los fluidos
de trabajo
Caracterización de llamas.......................................................................................................... 63
Revisión del modelo de gas ideal .............................................................................................. 65
Propiedades termodinámicas de los gases ideales ............................................................... 66
Mezclas de gases ideales ....................................................................................................... 68
Composición del aire y los combustibles .................................................................................. 69
Aire ........................................................................................................................................ 69
Combustibles ......................................................................................................................... 71
Combustión estequiométrica ..................................................................................................... 73
Balance energético de primera ley en un MCIA ....................................................................... 75
Calores de reacción............................................................................................................... 75
Entalpía de formación ........................................................................................................... 76
Poder calorífico..................................................................................................................... 77
Propiedades de los fluidos de trabajo ........................................................................................ 78
Cálculo de las propiedades ........................................................................................................ 80
Mezclas sin quemar ................................................................................................................... 81
Mezclas quemadas..................................................................................................................... 87
Referencias ................................................................................................................................ 94
Capítulo 5
Ciclos de trabajo teóricos y reales
Ciclos Teóricos .......................................................................................................................... 95
Ciclos con aportación de calor a volumen constante ........................................................... 97
Ciclos con aportación de calor a presión constante ............................................................. 97
Relaciones Termodinámicas para los diferentes procesos en motores.................................... 100
Relaciones para calcular ciclos de trabajo ideales ............................................................ 101
Ciclos Reales ........................................................................................................................... 102
v
Literatura recomendada ........................................................................................................... 107
Referencias .............................................................................................................................. 107
Capítulo 6
Proceso de renovación de la carga
Introducción............................................................................................................................. 108
Procesos de admisión y escape en motores de cuatro tiempos................................................ 108
Rendimiento volumétrico ........................................................................................................ 109
Efectos cuasi – estacionarios .............................................................................................. 110
Combinación de efectos dinámicos y cuasi – estacionarios ............................................... 113
Efectos dinámicos................................................................................................................ 116
Turboalimentación................................................................................................................... 118
Ecuaciones básicas.............................................................................................................. 120
Efecto de la altitud............................................................................................................... 124
Pérdidas mecánicas............................................................................................................. 126
El Compresor....................................................................................................................... 126
La Turbina ........................................................................................................................... 130
Literatura recomendada ........................................................................................................... 140
Referencias .............................................................................................................................. 140
Capítulo 7
Pérdidas de calor y refrigeración
Introducción............................................................................................................................. 141
Condiciones de funcionamiento a carga parcial ...................................................................... 144
Importancia de la transferencia de calor.................................................................................. 144
Cálculos de transferencia de calor en MCIA........................................................................... 145
Tasa de Transferencia de Calor en el cilindro ......................................................................... 147
Procedimiento de Cálculo de las pérdidas de calor Convectivo en el cilindro .................. 150
Cálculos de Refrigeración en el motor .................................................................................... 150
Procedimiento de cálculo del calor cedido al refrigerante................................................. 152
Balance de Energía en el motor............................................................................................... 152
Referencias .............................................................................................................................. 164
Capítulo 8
Pérdidas Mecánicas y Lubricación
Introducción............................................................................................................................. 165
Formas de medir las pérdidas de fricción en los MCIA .......................................................... 168
Pérdidas de fricción en el motor .............................................................................................. 170
Sistema de lubricación............................................................................................................. 172
vi
Requerimientos del lubricante............................................................................................. 173
Referencias .............................................................................................................................. 174
Capítulo 9
Capítulo 10
Capítulo 11
Capítulo 12
Capítulo 13
vii
CAPÍTULO 1
INTRODUCCIÓN
Desde siempre la primera energía mecánica utilizada por el hombre ha tenido el origen en su
propia energía muscular, incrementada a través de mecanismos simples como la palanca y la
polea. Otro paso importante en la obtención de la energía necesaria para realizar las labores más
extenuantes consistió en la domesticación de los animales para la sustitución de la energía
humana por la animal, en general más potente. A lo largo de los años la mejora progresiva de
diversos aparejos permitió ir incrementando la capacidad de tiro de los animales. La importancia
del aprovechamiento de la energía de los animales en la sustitución de tareas de los hombres se
comprende fácilmente teniendo en cuenta que aún hoy en día se sigue utilizando.
Un paso importante ha sido el aprovechamiento de las energías potenciales de la naturaleza
como la hidráulica y la eólica, que eran las más fácilmente detectables por el hombre. Así, las
primeras ruedas hidráulicas, precursoras de las actuales turbinas, constituyeron las primeras
actividades de producción en serie en las que se utilizó maquinaria.
El uso de la energía hidráulica pasó de los molinos de grano a las demás actividades
industriales, hasta tal punto que la revolución industrial se inició en un momento en el que las
máquinas hidráulicas, sobre todo motores, estaban muy extendidas. De hecho, como lo veremos
luego, los primeros motores de vapor se usaron en muchas ocasiones para bombear agua y mover
ruedas hidráulicas. Respecto a la energía eólica aparte de su uso convencional para mover
barcos, comenzó a utilizarse más tarde que la hidráulica con fines industriales.
Tanto la energía hidráulica como la eólica tienen el inconveniente de ser fuentes de energía de
carácter estacional cuya utilización está limitada a la presencia de saltos de agua y a la relativa
constancia de la velocidad y dirección del viento, respectivamente. La aparición de las primeras
máquinas de vapor potenció y aceleró el gran desarrollo tecnológico que dio lugar a la revolución
industrial a principios del siglo XVIII.
1
Adaptado de la literatura recomendada que aparece al final del capítulo
Capítulo 1. Breve Historia de los motores
Figura 1.1. Eolípila o primera turbina de reacción pura que se conoce en la historia. Siglo I a.C.
2
Capítulo 1. Breve Historia de los motores
residual (generación de vacío), junto con la presión atmosférica, actuando sobre la cara opuesta
del pistón empujaban a éste hacia abajo realizando la carrera de trabajo.
3
Capítulo 1. Breve Historia de los motores
provocaba un vacío que elevaba un volumen de agua que quedaba atrapado por encima de una
válvula antirretorno. En ese momento se readmitía vapor a la cámara y la presión del mismo
forzaba dicho volumen de agua hacia arriba pasando por una segunda válvula antirretorno. Hay
que señalar que las válvulas de admisión de vapor a la cámara y de agua de refrigeración para
provocar la condensación eran manuales. Además, la altura a la que podía elevarse el agua era
muy dependiente de la presión del vapor en la caldera. Se cree que el motor de Savery alcanzaba
los cinco ciclos por minuto y que su rendimiento térmico era muy inferior al 0.5%.
A pesar de las limitaciones que tenía, el motor de Savery tuvo un breve éxito comercial hasta
la aparición en 1712 de la primera máquina de vapor construida por Newcomen. Uno de los
inventos que ha influido de forma más decisiva en la historia de la humanidad es, sin duda, la
máquina de vapor, desarrollada por el herrero de profesión, Thomas Newcomen a principios del
siglo XVIII. El papel trascendental que ha desarrollado este motor térmico lo indica el hecho de
que fue el único utilizado por el hombre durante aproximadamente 150 años.
Una de las más importantes mejoras de la máquina de vapor construida por Newcomen,
ayudado por John Calley, fue la idea de provocar la condensación del vapor inyectando agua
pulverizada en el interior del cilindro. Se cuenta que una fuga accidental de agua hacia el interior
del cilindro a través de un poro produjo una condensación tan rápida del vapor que el pistón
rompió la cadena que la unía a un balancín, así como el fondo del cilindro y la parte superior de
la caldera. Este accidente providencial condujo al descubrimiento de la condensación por
mezcla.
4
Capítulo 1. Breve Historia de los motores
En la máquina de vapor de Newcomen (Figura 1.4), el vapor, generado en una caldera situada
bajo el cilindro de potencia, se introducía en el cilindro, elevando el pistón hasta su máxima
altura vertical. En este momento, se provocaba la condensación del vapor mediante una
inyección de agua fría en el interior del cilindro, de modo que por la acción de la presión
atmosférica aplicada en la otra cara del pistón (de ahí el nombre de motor atmosférico), éste
realizaba la carrera descendente de trabajo.
Entre las innovaciones introducidas en la máquina de vapor de Newcomen cabe destacar las
siguientes:
• Empleo del balancín para transmitir el movimiento lineal del pistón
• Válvulas accionadas automáticamente desde el balancín
• Válvula para expulsar el aire del interior del cilindro
• Empleo de segmentos (anillos) de cuero, sobre los que se mantenía un cierto nivel de
agua, para asegurar el sellado entre el cilindro y el pistón. Este dispositivo constituye el
primer cierre hidráulico de la historia de la tecnología.
Las mejoras introducidas por Smeaton alrededor de 1768 aumentaron el valor del
rendimiento de la máquina de vapor hasta un 1%. En esta época, su funcionamiento era ya muy
5
Capítulo 1. Breve Historia de los motores
fiable, por lo que, sólo en Inglaterra se construyeron aproximadamente 200 unidades. Fueron, sin
embargo, las mejoras introducidas por el escocés James Watt, el cual, trabajando en la
Universidad de Glasgow, intervino en la reparación de un modelo reducido de máquina de vapor.
El interés que despertó en él este motor provocó que dedicase el resto de su vida a su mejora y
construcción.
Quizá la mejora más sustancial que introdujo Watt en la máquina de Newcomen fue la
invención del condensador separado (Figura 1.5). Al hacer esto se conseguía eliminar el proceso
repetitivo de recalentar el cilindro en cada ciclo. Se puede decir que esta mejora elevó por sí
misma el rendimiento al doble del logrado por Smeaton y a cuatro veces el original de
Newcomen. La invención del condensador separado proporcionó a Watt una patente en 1769,
cuya duración fue extendida en 1775 hasta el año 1800. El título de la misma fue “Método para
reducir el consumo de vapor y combustible de las máquinas de fuego”, fechas entre las cuales en
unión con Mathew Boulton lograron poner en operación cerca de 500 motores.
6
Capítulo 1. Breve Historia de los motores
• La primera definición de una unidad precisa de potencia, que fue llamada “horsepower”,
definida como la potencia necesaria para levantar 33.000 libras a un pie de altura durante
un minuto.
El uso de la máquina de vapor se popularizó tanto que llegó a ser empleada en la navegación
hacia 1807 por R. Fulton y mas tarde, en 1825 G. Stephenson la usó para accionar una
locomotora. En la fotografía 1.1 se observa una máquina de vapor empleada en el banco de
España para estampación de moneda. Esta máquina fue donada a la Escuela Técnica Superior de
Ingenieros Industriales de Madrid y se encuentra actualmente en el salón principal del edificio
como pieza de museo.
La aparición de la turbina de vapor a finales del siglo XIX desplazó a la máquina de vapor del
campo de las grandes potencias, aunque ésta siguió utilizándose en aplicaciones con grados de
carga y regímenes variables. Sin embargo, fue la aparición de los motores de combustión interna
alternativos el hecho que constituyó el principio del fin de la máquina de vapor, que desapareció
prácticamente del mercado a mediados del siglo actual.
Fotografía 1.1. Máquina de vapor ubicada como pieza de museo en la ETSII de Madrid
7
Capítulo 1. Breve Historia de los motores
Figura 1.6. Esquema del motor Stirling y funcionamiento básico del motor
8
Capítulo 1. Breve Historia de los motores
en la cámara B donde se calienta de nuevo al entrar en contacto con las paredes superiores del
cilindro expuestas al calor procedente de un horno. Gracias a este calentamiento el aire se
expande en la cámara B y mueve el pistón con el regenerador hacia abajo, sin embargo, sigue
pasando aire de la cámara A hacia la B, ya que el pistón de potencia sigue su carrera ascendente,
hasta que ambos pistones se juntan en la posición superior del pistón de potencia. Con todo el
aire contenido ahora en la cámara B, se produce la expansión del mismo, empujando hacia abajo
ambos pistones hacia la posición inferior. En este momento, el pistón con el regenerador es
elevado por la inercia del motor, mientras el pistón de potencia permanece quieto. El
regenerador, que se ha refrigerado al proporcionar calor al aire que ha pasado a su través, absorbe
calor del aire contenido en la cámara B que es forzado a pasar hacia abajo a través de él. De este
modo el regenerador se vuelve a calentar y el aire se enfría al entrar en la cámara A, iniciándose
un nuevo ciclo.
Los motores de aire caliente se dejaron de fabricar a principios del siglo actual, eclipsados por
los motores de combustión interna alternativos. Sin embargo, los avances en metalurgia y en la
tecnología de transmisión de calor han revitalizado el interés en el motor de Stirling, debido
principalmente a sus bajas emisiones contaminantes y bajo ruido. Así, no se puede descartar que,
en el futuro, este motor pueda suponer una cierta competencia para los actuales.
9
Capítulo 1. Breve Historia de los motores
Sin embargo, fue hasta muy entrado el siglo XIV cuando se logró rebatir completamente la
teoría del Calórico. Se puede decir que la ciencia que explica los procesos que ocurren en los
motores térmicos, la termodinámica, se encuentra en los pensamientos del ingeniero militar
francés Nicholas Leonard Sadi Carnot, publicados en 1824 con el título de “Reflexiones sobre la
potencia motriz del fuego y sobre las máquinas adecuadas a recibir esta potencia”. De esta
manera, y hacia finales de 1860, con la aportación de grandes genios de la historia, se logra
plantear la termodinámica como cuerpo de doctrina.
En este motor, el aire se inducía al cilindro bombeándolo con la mano, haciendo que el pistón
subiera un poco. Luego se inyectaba combustible líquido y se producía el encendido debido al
calentamiento de la parte inferior del cilindro ubicado en un horno. La combustión y expansión
de los gases forzaban el pistón hacia arriba y realizaba trabajo, no solamente para elevar el
pesado pistón, sino también para bombear agua a una elevación superior.
La parte superior del cilindro estaba rodeado de una chaqueta de agua para refrigerar los
gases de combustión y producir así un vacío parcial, el cual, combinado con el efecto de la
gravedad, hacían que el pistón se impulsara hacia abajo y elevara el agua desde el pozo.
10
Capítulo 1. Breve Historia de los motores
32 años más tarde, en 1826, el inglés Samuel Brown era acreditado como el inventor del
motor de vacío, usaba un gran tanque lleno de aire dentro del cual inyectaba combustible
gaseoso. El encendido de la mezcla producía una llama que forzaba el aire hacia fuera del
tanque. La válvula de escape del tanque se cerraba cuando la presión disminuía hasta la
atmosférica, y en ese momento inyectaba un chorro de agua dentro de los gases calientes
provocando un enfriamiento rápido de los gases y obteniendo un vacío parcial en el tanque. Una
cara del pistón estaba sometida a la presión atmosférica y la otra al interior del tanque, de tal
forma que la diferencia de presiones unidas al peso de éste, proveían la fuerza para mover el
pistón y realizar trabajo.
El 13 de mayo de 1854 los Padres Eugenio Barsanti y Felix Matteucci patentaban en
Inglaterra su método de “obtención de potencia motriz por explosión de gases”. En la figura 1.8
muestra el pistón y la biela descendiendo debido al efecto del enfriamiento de gases, al peso del
pistón y a la presión atmosférica hasta que la cara superior del pistón descubre el puerto superior
en la parte inferior izquierda del cilindro. Esto obliga a pasar los gases de combustión hacia la
parte superior del pistón. La parte superior del cilindro está abierto a la atmósfera, de tal forma
que los productos de combustión se mezclan con el aire inducido dentro del cilindro durante la
carrera de bajada del pistón. La aplicación de la patente requiere del uso de dos cilindros, así que
la bajada del segundo pistón inicia la subida del primero y el aire es inducido a través del puerto
empleado para escape. Se inyecta hidrógeno a través de una válvula ubicada en la parte inferior
derecha del cilindro. El movimiento del mecanismo de esta válvula produce una chispa selectiva
la cuál enciende los gases. El pistón y la biela eran libres durante la carrera de subida, pero
realizaban trabajo en la carrera descendente.
Figura 1.8. Vista parcial del motor de pistón libre de Barsanti y Matteucci
11
Capítulo 1. Breve Historia de los motores
El motor de Lenoir
En 1860 el francés Étienne Lenoir desarrolló el motor de gas sin compresión (figura 1.9)
fuertemente inspirado en las máquinas de vapor de la época. Se usaba una mezcla aire y gas de
hulla que se quemaba a presión atmosférica. La mezcla aire – gas se admitía durante parte de la
carrera de admisión a través de una válvula de corredera, cuya apertura era controlada por un
sistema de varillas movido desde el cigüeñal. Cuando el pistón había realizado aproximadamente
la mitad de la carrera, las válvulas se cerraban y, poco después, saltaba una chispa, originada por
una bujía, lo que inflamaba la mezcla elevando la presión y transmitiendo potencia al pistón
durante la segunda mitad de la carrera. En la carrera de regreso se descargaban los gases del
cilindro. La eficiencia de este motor era baja debido a lo corta que era la carrera de expansión y a
la ausencia de compresión previa de la mezcla. Pese a esto, se construyeron cerca de 500
motores de este tipo en Francia e Inglaterra.
En 1866 los alemanes Otto y Langen construyen una versión mejorada del motor de Barsanti
y Matteucci, incluyendo además las mejoras del motor de Lenoir como la válvula de corredera y
el encendido provocado de la mezcla. El motor Otto-Langen, como se llamó después, empleaba
una volante de inercia para subir el pistón desde la posición inferior e inducir carga fresca. Este
motor aprovechaba el incremento de presión por la combustión de la mezcla al principio de la
carrera de subida de un pistón libre para acelerar dicho pistón, que estaba solidario a una
cremallera. La subida rápida del pistón libre, que tenía una inercia considerable, producía la
formación de vacío en el interior del cilindro. Posteriormente, con la cremallera del pistón
engranada al cigüeñal, la presión atmosférica empujaba al pistón en su carrera descendente, que
era la que proporcionaba trabajo al pistón. Por otra parte la válvula de corredera controlaba los
12
Capítulo 1. Breve Historia de los motores
Las primeras dos condiciones van encaminadas a reducir al mínimo las pérdidas de calor,
permitiendo así que toda la energía disponible en los gases de escape se transforme en trabajo.
La tercera condición prevé la expansión de los gases hasta el máximo, de tal manera que se
obtenga el máximo trabajo de la expansión. La cuarta condición reconoce el hecho de que entre
más alta sea la presión inicial, mayor es la presión a través de la carrera para una relación de
expansión dada y hace posible mayores relaciones de expansión, resultado ambas en más trabajo.
Beau de Rochas, también indicó el método de operación deseado para el motor de combustión
interna a pistón:
Han transcurrido cerca de 140 años y los principios siguen siendo prácticamente los mismos,
únicamente se han variado los tiempos de aperturas y cierres de válvulas con el fin de
incrementar las velocidades y las potencias.
13
Capítulo 1. Breve Historia de los motores
El motor de Brayton
En 1873, el estadounidense Brayton, desarrolló un motor que tenía las características de proceso
de combustión a presión constante y expansión completa. Se usaba un cilindro para comprimir
aire o la mezcla aire – combustible; el otro cilindro se usaba como un cilindro de trabajo y era lo
suficientemente grande para obtener la expansión completa hasta presión atmosférica. El
compresor descargaba la mezcla a una cámara, y la mezcla fluía luego de ésta al motor, siendo
encendida y quemada a presión constante a medida que entraba al motor (ver figura 1.10).
Figura 1.10. Esquema del motor de Brayton y representación en el P – V del proceso ideal
Se mantenía una llama para el encendido de la mezcla mediante un baypass desde la cámara ,
y como medida de seguridad, para evitar la retrollama a la cámara se ponía una malla.
El motor de Brayton, aunque construido relativamente en gran número, no representó una
competencia considerable con el motor de pistón libre de Otto – Langen debido principalmente a
su bajo rendimiento por las elevadas pérdidas de calor y de fricción. Posteriormente éste se
abandonó completamente a favor del proceso de Beau de Rochas cuando el motor de Otto
penetró en el mercado estadounidense. Teóricamente el ciclo de Otto y del Brayton deberían
tener el mismo rendimiento térmico. Aunque el proceso de Brayton fue desechado para su
aplicación en motores de pistones, fue sin embargo usado como proceso para las turbinas de gas y
por esta razón a veces se dice que éstas trabajan según el ciclo de Brayton.
El motor Otto
A pesar de todos los diseños de motores de combustión interna alternativos descritos hasta ahora,
la máquina de vapor (e incluso los coches de caballos) era aun el motor más utilizado. Sin
embargo, esta situación cambió rotundamente en el año 1876, cuando el alemán Nikolaus
14
Capítulo 1. Breve Historia de los motores
August Otto (ver fotografía 1.2) construyó el primer motor de cuatro tiempos de la historia (ver
figura 1.11) basado en los principios de Beau de Rochas. Este motor ingresaría al banco de
patentes con el número 365,701.
A pesar de que la patente de Otto fue invalidada en Alemania porque se basaba en los
principios de Beau de Rochas, nadie acusó a éste de haber copiado, sino que simplemente se
pensó que a Otto se le había ocurrido el ciclo de cuatro tiempos 14 años después que a Beau de
Rochas. Las repercusiones de este litigio llegan hasta nuestros días, ya que el ciclo de cuatro
tiempos se conoce en Francia como ciclo de Beau de Rochas, mientras que, prácticamente en el
resto del mundo, se le llama ciclo Otto.
15
Capítulo 1. Breve Historia de los motores
Los motores de dos tiempos de compresión en carcaza (o en cárter) se usan más comúnmente
debido a su bajo costo y pequeño tamaño especialmente en motores de fuera de borda y en
motocicletas. Los motores de dos tiempos con soplador de barrido independiente (figura 1.12)
se emplean en motores grandes para trabajo pesado. En la actualidad los primeros tienden a
desaparecer del mercado por su excesiva descarga de emisiones contaminantes a la atmósfera.
El motor Diesel
En 1892 el alemán Rudolf Diesel, (ver fotografía 1.3) basado en que el trabajo de compresión
rápida en un medio incrementa su energía interna y su temperatura, propuso comprimir sólo aire
hasta que alcanzara una temperatura lo suficientemente elevada para que lograra encender el
16
Capítulo 1. Breve Historia de los motores
Así pues, Diesel recibiría el 16 de julio de 1895 en los Estados Unidos la patente número
542.846 por este diseño, usando carbón pulverizado. Esta patente igualmente incluía el sistema
de inyección de combustible líquido en el cual una bomba suministraba la cantidad de
combustible deseado mediante una válvula de aguja y aire a presión arrastraba el combustible
dentro del cilindro.
El motor diesel admite y comprime hasta un instante antes de iniciarse la inyección del
combustible, únicamente aire, lo que le permitió tener relaciones de compresión mucho más
elevadas que las de todos los motores de la época (los cuales presentaban problemas de
detonación de la premezcla aire - combustible al interior de la cámara de combustión, debido a
las características propias del combustible que empleaban) y, por tanto, se pudo obtener un
rendimiento más elevado y un menor consumo de combustible. El primer motor experimental de
Diesel se construyó en 1893, pero nunca llegó a funcionar. Después de superar algunos
problemas, el primer motor operativo se construyó en 1897 (Fotografía 1.4). Este motor
proporcionaba 17,8 CV a 154 rpm con un rendimiento del 26,2% valor muy elevado para aquella
época.
17
Capítulo 1. Breve Historia de los motores
En la fotografía 1.5 se puede observar a Rudolf Diesel con sus ayudantes, en su taller de
Augsburg (Alemania) preparando sus experimentos.
18
Capítulo 1. Breve Historia de los motores
En la fotografía 1.6 se observa la primera planta de generación eléctrica a gran escala con
motores Diesel instalada en Kiev (Rusia) en el año de 1904.
Fotografía 1.6. Planta de generación eléctrica empleando motores Diesel. Kiev (Rusia) 1904
Se puede afirmar que después de la aparición del motor Diesel se completó el desarrollo
básico de los motores de combustión interna alternativos con respecto a su concepción mecánica.
El diseño de estos motores durante el siglo actual ha estado marcado por la mejora de
prestaciones (potencia y rendimiento – consumo, sobre todo), conseguida en gran parte por el
importante desarrollo de los combustibles, sobre todo en el período que siguió a la 2da. Guerra
Mundial con el descubrimiento de sustancias antidetonantes, que condujo a un mayor
conocimiento de la forma como la naturaleza del combustible afectaba la combustión.
Otro factor más recientes que ha influido en el diseño de los motores de combustión interna
alternativos es la disminución de las emisiones contaminantes químicas y acústicas.
El motor Wankel
Hacia 1956 el ingeniero alemán Federico Wankel inventaba un nuevo motor de combustión
interna cuyo principio no era alternativo sino rotativo. Es un motor de encendido provocado
mecánicamente distinto a los anteriormente mencionados. Como se muestra en la figura 1.13,
19
Capítulo 1. Breve Historia de los motores
este motor consta de un cárter provisto de una cámara de forma especial dentro de la cual gira un
émbolo rotativo que tiene forma de triángulo curvilíneo; el émbolo está dotado de un dentado
interior que engrana con un piñón solidario al cigüeñal del motor, sobre el cual gira el émbolo
excéntricamente. La gran ventaja del motor Wankel, como se le suele llamar en memoria a su
inventor, consiste en que sus piezas no suben ni bajan, sino que giran. Pesa menos, es más
compacto y tiene menor número de partes móviles que el motor convencional. El eje de salida
gira al triple de velocidad que el rotor triangular. El rotor mantiene sus tres vértices en continuo
contacto con la superficie interior lo que hace difícil el sellado de las cámaras de combustión.
Entre otras desventajas se encuentran las elevadas fugas y las grandes pérdidas por transferencia
de calor. En la figura 1.14 se aprecia el funcionamiento del motor Wankel.
Actualmente este tipo de motores lo llevan los japoneses Mazda en su serie RX. A nivel
experimental, esta misma casa ha sacado el Mazda Miata de hidrógeno con motor Wankel,
conocido como el MX – 5.
20
Capítulo 1. Breve Historia de los motores
El motor Wankel, sin duda, ha sido un gran marginado en la historia reciente de los motores.
Los principales fabricantes de motores en el mundo, han hecho grandes esfuerzos en optimizar
los motores convencionales de pistones y han dejado de lado este prometedor prototipo.
Investigando más en mejorar sus desventajas, especialmente las fugas y el sellado, podría llegar a
obtenerse una gran potencia específica (relación potencia / peso o potencia / volumen) propia de
la concepción original del ingeniero F. Wankel.
21
Capítulo 1. Breve Historia de los motores
LITERATURA RECOMENDADA
1. Burón, J.M., “Proyecto docente para optar al cargo de profesor titular en la Universidad
Politécnica de Madrid”, Escuela Técnica Superior de Ingenieros Industriales (ETSII),
Madrid, 1996
2. Obert F. Edward, “Motores de Combustión Interna”, Cecsa, Mexico, 1997
3. Heywood J.B., “Internal Combustion Engine Fundamentals”, McGraw-Hill, New York, 1989
4. Lichty, L.C. “Combustion Engine Processes”, International student edition, McGraw-Hill,
1967
REFERENCIAS
[1] Salanki, P.A., Wallace, J.S., (1996 ), “Evaluation of the hydrogen fueled rotary engine for
hybrid vehicle applications”, SAE paper No. 960204.
22
CAPÍTULO 2
1
Los números entre corchetes indican la referencia bibliográfica en orden de aparición al final del capítulo.
Capítulo 2. Tipos de motores y su funcionamiento
como una máquina hidráulica. En las máquinas térmicas, por el contrario, evolucionan fluidos
compresibles que tienen una compresibilidad no despreciable. La compresibilidad juega un papel
muy importante en el intercambio energético que tiene lugar entre el fluido y el eje de la
máquina, ya que, la variación del volumen específico es el mecanismo que permite la
transformación de energía química en mecánica y, por tanto, su posterior aparición en el eje de la
máquina.
Las máquinas de fluido se pueden clasificar además en otros dos grupos característicos:
Máquinas de desplazamiento positivo o volumétricas y Turbomáquinas. En las primeras existe
una cierta cantidad bien definida de fluido que atraviesa la máquina en cada instante. En las
turbomáquinas, por el contrario, el volumen o la masa desplazada no está materializada por un
contorno definido, sino que el flujo es continuo. En éstas últimas la transferencia de energía del
fluido al eje se basa en el teorema del momento cinético, que conduce a la ecuación de Euler ó
ecuación básica de las turbomáquinas.
Las máquinas de desplazamiento positivo las podemos a su vez dividir en máquinas
alternativas (de cilindro y pistón) y máquinas rotativas.
24
Capítulo 2. Tipos de motores y su funcionamiento
25
Capítulo 2. Tipos de motores y su funcionamiento
Todas las clasificaciones anteriores son importantes e ilustran el amplio espectro de diseños
disponibles de motores. A lo largo de éste texto nos guiaremos fundamentalmente desde el punto
de vista del método de encendido MEC o MEP. En importancia seguirá, el combustible
empleado, el método de preparación de la mezcla, el diseño de la cámara de combustión, el
método de control de la carga, algunos detalles del proceso de combustión, y finalmente las
características de operación y emisiones de los motores.
26
Capítulo 2. Tipos de motores y su funcionamiento
2. Carrera de compresión, ocurre cuando las dos válvulas están cerradas y la mezcla al interior
del cilindro es comprimida hasta una pequeña fracción de su volumen inicial. Hacia el final
de la carrera de compresión, se inicia la combustión y la presión en el cilindro incrementa
más rápidamente.
3. Carrera de potencia, o carrera de expansión, la cual inicia con el pistón en el PMS y termina
en el PMI en la medida que los gases a elevada temperatura y presión empujan el pistón hacia
abajo forzando a girar la manivela. Durante ésta carrera se logra obtener cerca de cinco veces
más trabajo que el realizado por el pistón durante la compresión. En la medida en que el
pistón se aproxima al PMI la válvula de escape se abre para iniciar el proceso de escape y cae
la presión en el cilindro hasta un valor cercano a la presión en el escape.
4. Carrera de escape, donde los gases quemados salen del cilindro: primero, porque la presión
en el cilindro puede ser sustancialmente más alta que la presión en el escape y luego porque
son barridos por el propio movimiento del pistón hacia el PMS. A medida que el pistón se
aproxima al PMS abre la válvula de admisión. Justo después del PMS se cierra la válvula de
escape y el ciclo inicia de nuevo.
Los motores de cuatro tiempos requieren, para cada cilindro, dos revoluciones del cigüeñal
por cada carrera de potencia. Para obtener potencias más elevadas a partir de un tamaño de
27
Capítulo 2. Tipos de motores y su funcionamiento
motor dado, se desarrolló el motor con ciclo de trabajo de dos tiempos. Los dos tiempos son
aplicables tanto a MEP como a MEC.
En la Figura 2.3 se muestra el ciclo de operación de dos tiempos. Los puertos en la camisa
del cilindro son abiertos y cerrados por el movimiento del pistón controlando así los flujos de
admisión y escape a medida que el pistón se acerca al PMI. Los dos tiempos son:
1. Carrera de compresión, la cual inicia con el cierre de los puertos de admisión y escape y
entonces comprime el contenido del cilindro mientras se arrastra carga fresca al interior del
cárter. A medida que el pistón se acerca al PMS se inicia la combustión.
2. Carrera de potencia o expansión, similar a los cuatro tiempos hasta que el pistón se aproxima
al PMI, cuando primero se descubre el puerto de escape y luego el de admisión. Cuando los
puertos de admisión de descubren la carga fresca que había sido comprimida en el cárter fluye
hacia el cilindro. El pistón y los puertos generalmente tienen formas especiales para cambiar
bruscamente la dirección del flujo que entra permitiendo así un barrido efectivo de los gases
residuales.
Figura 2.3 Ciclo de operación de dos tiempos. Barrido por cárter [2]
Para lograr un ciclo completo de potencia del motor se requiere solamente una revolución del
cigüeñal. Sin embargo, existen grandes dificultades de llenar completamente el volumen
desplazado con carga fresca, y parte de ésta fluye directamente hacia fuera del cilindro durante el
proceso de barrido siendo ésta la principal razón por la que los MEP de dos tiempos están en gran
desventaja debido a que la carga fresca contiene combustible y aire.
28
Capítulo 2. Tipos de motores y su funcionamiento
29
Capítulo 2. Tipos de motores y su funcionamiento
Figura 2.4 Dibujo en corte de un motor mostrando sus principales componentes. [2]
La culata sella los cilindros y se fabrica generalmente en hierro fundido o en aluminio. Esta
debe ser lo suficientemente fuerte y rígida para distribuir la fuerza de los gases que actúan de la
manera más uniformemente posible a través del bloque. La culata contiene las bujías si es un
MEP o los inyectores si es un MEC y en motores que tienen válvulas sobrecabeza suele llevar
partes del mecanismo de las válvulas.
30
Capítulo 2. Tipos de motores y su funcionamiento
31
Capítulo 2. Tipos de motores y su funcionamiento
potencia requerida (a cualquier velocidad del motor) está por debajo de la máxima, la cual se
obtiene con la mariposa totalmente abierta. El colector de admisión generalmente se calienta
para facilitar la evaporación rápida del combustible líquido y obtener así una distribución
uniforme de combustible entre todos los cilindros.
A partir de 1998 se exige en Colombia por ley, el uso de sistemas de inyección electrónico en
lugar de la carburación con fines de mejorar las emisiones de los vehículos. Los sistemas de
inyección de combustible pueden ser monopunto cuando se trata de un solo inyector ubicado en
el colector de admisión o multipunto cuando se usan tantos inyectores como cilindros ubicados
en los puertos de admisión. En el sistema de inyección controlado electrónicamente el flujo de
aire se mide directamente; las válvulas de inyección son actuadas dos veces por cada revolución
del árbol de levas por pulsos de inyección cuya duración viene determinada por la unidad de
control electrónica para suministrar la cantidad deseada de combustible por ciclo y por cilindro.
La secuencia de eventos que tienen lugar al interior del cilindro se ilustra en la Figura 2.5
Figura 2.5 Secuencia de eventos en un MEP de 4 tiempos en función del ángulo de giro del cigüeñal [2]
RCA: Retraso al cierre de la válvula de admisión; AAA: Avance de Apertura de la válvula de Admisión;
AAE: Avance de Apertura de la válvula de Escape; RCE: Retraso al cierre de la válvula de escape
32
Capítulo 2. Tipos de motores y su funcionamiento
Como se observa en la Figura 2.5 es más cómodo utilizar como variable independiente el
ángulo de giro del cigüeñal en lugar del tiempo debido a que la mayoría de los procesos en el
motor ocurren a intervalos de ángulo de giro del cigüeñal casi constante sobre un amplio rango de
velocidades del motor. Para mantener flujos de mezcla elevados a altas velocidades (y por lo
tanto altas potencias) la válvula de admisión que abre antes del PMS y cierra bastante después
del PMI, produciéndose de ésta manera una mezcla entre la carga fresca y los gases quemados
residuales que permanecieron del ciclo anterior.
Una vez cierra la válvula de admisión el contenido del cilindro es comprimido por encima de
la presión y temperatura atmosférica a medida que el volumen del cilindro es reducido. Se
genera alguna transferencia de calor hacia el pistón, la culata y las paredes del cilindro, pero el
efecto sobre las propiedades del gas sin quemar aún es despreciable.
Entre 10 y 40 grados de cigüeñal antes del PMS, la chispa genera una descarga eléctrica entre
los dos electrodos iniciándose el proceso de combustión. Un interruptor rotativo o distribuidor
arrastrado por el árbol de levas interrumpe la corriente desde la batería hacia el circuito primario
de la bobina. El cableado secundario de la bobina de encendido conectado a la bujía produce un
voltaje elevado a través de los electrodos a medida que se colapsa el campo magnético.
Actualmente el distribuidor es electrónico. La descarga eléctrica hace que se desarrolle un frente
de llama turbulento que se propaga a través de la mezcla aire, combustible y gases residuales en
el cilindro y se extingue en las paredes de la cámara de combustión. La duración de éste proceso
varía según el diseño y la operación del motor pero generalmente está entre 40 y 60 grados de
ángulo de giro del cigüeñal como se observa en la Figura 2.5. A medida que se quema la mezcla
en el frente de llama, incrementa la presión en el cilindro por encima del nivel que lo haría si solo
se comprimiera (línea punteada). Esta última curva también llamada curva de presión en cilindro
a motor arrastrado, se obtiene arrastrando el motor desde el exterior o sin combustión. En la
realidad los procesos de admisión y compresión del motor con combustión y del motor arrastrado
no son exactamente las mismas debido a la presencia de gases quemados del ciclo anterior. Es
importante notar que debido a las diferencias en los patrones de flujo y a la composición de la
mezcla entre cilindros, y dentro de cada cilindro ciclo a ciclo, cada proceso de combustión difiere
en algo. Este fenómeno se conoce como dispersión cíclica. Como consecuencia de esto, la forma
de la curva de presión versus ángulo de cigüeñal en cada cilindro y ciclo a ciclo no es
exactamente la misma.
Existe un momento óptimo para el salto de la chispa el cual, para una masa dada de aire +
combustible en el interior del cilindro genera el par máximo. Un avance muy alejado (más
retirado del PMS) o muy retardado (muy cercano al PMS) del óptimo daría siempre menor
33
Capítulo 2. Tipos de motores y su funcionamiento
34
Capítulo 2. Tipos de motores y su funcionamiento
motores de pequeño tamaño los dos tiempos son bastante competitivos con los cuatro tiempos, en
gran parte debido, a que en el ciclo diesel únicamente se producen pérdidas de aire durante el
proceso de barrido. La operación de un MEC típico de aspiración natural de cuatro tiempos se
ilustra en la Figura 2.6.
La relación de compresión de los motores diesel es mucho más alta que la de los MEP típicos,
estando entre el rango de 12 a 24 dependiendo del tipo de MEC y de si el motor es de aspiración
natural o turboalimentado. La sincronización de válvulas usada es similar a los MEP. El aire a
una presión cercana a la atmosférica es inducido durante la carrera de admisión y luego
comprimido a una presión cercana a los 4000 kPa (40 bar) y a una temperatura cercana a los
800º K durante la carrera de compresión. Unos 20 grados antes del PMS comienza la inyección
de combustible al cilindro. Un perfil típico de la tasa de inyección de combustible se puede ver
en la Figura 2.6 b. El chorro de combustible líquido se atomiza en gotas y penetra en el aire.
Las gotas van reduciendo de tamaño a medida que el combustible líquido se evapora; el vapor de
combustible se mezcla luego con el aire. La temperatura y presión del aire están por encima del
punto de encendido del combustible. Por lo tanto después de un corto período de retraso, ocurre
un encendido espontáneo (autoignición) de aquellas regiones dentro del cilindro en las que se han
logrado las proporciones adecuadas aire – combustible iniciándose y propagándose el proceso de
combustión. La presión en el cilindro sube por encima del nivel de presión sin combustión
(Figura 2.6 c). La llama se esparce rápidamente a través de aquellas porciones de aire –
combustible que ya se han mezclado. A medida que se da el proceso de expansión continua la
mezcla entre el aire, el combustible y los gases quemados, acompañados por más combustión
(Figura 2.6 d). A plena carga, la masa de combustible inyectado es aproximadamente un 5% de
la masa de aire en el cilindro. La cantidad de combustible que puede ser quemado eficientemente
está limitado por los niveles de humo negro en el escape. El proceso de escape es similar al del
MEP de cuatro tiempos. AL final de la carrera de escape el ciclo inicia de nuevo.
En el MEC de dos tiempos la compresión, inyección de combustible, combustión y expansión
son similares a su equivalente de cuatro tiempos; son las presiones en la admisión y en el escape
las que difieren.
El sistema de inyección de combustible convencional de un motor diesel consiste de una
bomba de inyección, tuberías de suministro de combustible, y los inyectores. Existen diferentes
tipos de bombas de inyección y de inyectores.
35
Capítulo 2. Tipos de motores y su funcionamiento
Figura 2.6 Secuencia de eventos en un MEC típico de cuatro tiempos de aspiración natural [2]
Donde RCA y AAE son las mismas definiciones de la Figura 2.5. AII: Angulo de Inicio de la Inyección;
AIF: Angulo de Fin de la Inyección; AIC: Angulo de Inicio de la Combustión y AFC: Angulo de Fin de la
Combustión
En una bomba de inyección corriente (como el diseño de bomba en línea que se muestra en la
Figura 2.7) un conjunto de émbolos movidos por levas, uno por cada cilindro, operan en
cilindros muy ajustados. Al inicio de la carrera del émbolo el puerto de entrada está cerrado y el
combustible atrapado por encima de éste es forzado a través de una válvula de cheque hacia la
36
Capítulo 2. Tipos de motores y su funcionamiento
línea de inyección. La tobera de inyección tiene generalmente varios agujeros a través de los
cuales el chorro de combustible fluye al cilindro (ver Figura 2.8). Una válvula normalmente con
carga de un resorte mantiene cerrados esos agujeros hasta que la presión en la línea de inyección
actuando sobre parte de la superficie de la válvula vence la fuerza del resorte y abre la válvula.
La inyección inicia ligeramente después de que la presión en la línea empieza a aumentar. Así
pues, El inicio de la inyección está controlado por la puesta en fase entre el árbol de levas de la
bomba y el cigüeñal. La inyección finaliza cuando el puerto de admisión de la bomba es
descubierto por una ranura helicoidal labrada en el émbolo de la bomba, debido a que se libera la
elevada presión ejercida por encima de éste. La cantidad de combustible inyectado, (la cual
controla la carga) está determinada por el diseño de la leva de la bomba de inyección y por la
posición de la ranura helicoidal. Así pues, para un diseño de leva dado, se puede variar la carga
rotando el émbolo.
Los MEC pequeños, suelen llevar bombas de inyección centrífugas. Estas tienen únicamente
un émbolo el cual mide y distribuye el combustible hacia todas las toberas de inyección. En el
Capítulo 9 se detalla más este tipo de bombas. La unidad contiene una bomba de combustible de
baja presión a la izquierda y una bomba de inyección de alta presión a la derecha, un regulador de
37
Capítulo 2. Tipos de motores y su funcionamiento
Figura 2.8 Detalles del sistema de inyección de combustible de una bomba en línea de un MEC [2]
38
Capítulo 2. Tipos de motores y su funcionamiento
Estos motores a menudo son llamados motores de carga estratificada debido a la necesidad de
producir en el proceso de mezclado entre el chorro de combustible y el aire en el cilindro una
mezcla estratificada, con una composición fácilmente encendible por la bujía. Debido a que éstos
motores no requieren como los MEP combustibles con elevado octanaje ni combustibles con una
elevada calidad de encendido como los MEC, son por lo tanto apropiados para muchos tipos de
combustibles y operan con un amplio rango de combustibles líquidos.
En la Figura 2.9 se ilustran los principios de motores de carga estratificada
multicombustibles. Generalmente la cámara de combustión viene labrada en forma de copa en la
corona del pistón y requiere un alto grado de torbellino creado durante la admisión y aumentado
durante la compresión para alcanzar un mezclado rápido entre el aire y el combustible.
39
Capítulo 2. Tipos de motores y su funcionamiento
LITERATURA RECOMENDADA
REFERENCIAS
[1] Muñoz Torralbo, M. y Payri Gonzáles, F. (1978), “Turbomáquinas Térmicas”. Sección de
publicaciones de la Escuela Técnica Superior de Ingenieros Industriales de la Universidad
Politécnica de Madrid, Madrid.
[2] Heywood, J. B. (1988), “Internal Combustion Engine Fundamentals”. McGraw Hill, New
York.
40
CAPÍTULO 3
Este tomo, entendido como guía para estudiantes de pregrado de ingeniería, está relacionado
fundamentalmente con el funcionamiento, rendimiento y emisiones características de los
motores; la omisión de los otros factores listados anteriormente no reduce, de ninguna manera, su
gran importancia.
El funcionamiento del motor está definido más precisamente por:
1. La máxima potencia (o el par máximo) disponible en cada velocidad dentro del rango de
operación del motor.
2. El rango de velocidades y potencia sobre los que la operación del motor es satisfactoria
Figura 3.1 Geometría del cilindro, pistón, biela, y cigüeñal donde B = diámetro, L = carrera, l =
longitud de la biela, a = manivela, θ = ángulo del cigüeñal. [1]
42
Capítulo 3. Parámetros de operación y de diseño de motores
(
s = aCos θ + l 2 − a 2 Sen 2θ ) 1/ 2
(3.5)
El ángulo θ, definido como se muestra en la Figura 3.1, es llamado ángulo del cigüeñal. La
ecuación (3.4) con las definiciones anteriores puede ser reacomodada:
V
Vc
1
2
⎡
⎣
(
= 1 + (rc − 1)⎢ R + 1 − Cosθ − R 2 − Sen θ
2
)
1/ 2 ⎤
⎥
⎦
(3.6)
43
Capítulo 3. Parámetros de operación y de diseño de motores
La Figura 3.2 muestra como Sp varía con cada carrera para R = 3.5.
La resistencia al flujo de gas dentro del motor o los esfuerzos debido a la inercia de las partes en
movimiento limita la máxima velocidad media del pistón dentro del rango de 8 a 15 m/s. Los
motores de automóviles operan al extremo más alto de este rango; el extremo más bajo es típico
de los grandes motores marinos diesel.
Figura 3.2 Velocidad instantánea del pistón/velocidad media del pistón como función del ángulo del
cigüeñal para R = 3.5. [1]
M = Fb (3.12)
44
Capítulo 3. Parámetros de operación y de diseño de motores
La potencia P entregada por el motor y absorbida por el dinamómetro es el producto del par y
la velocidad angular:
P = 2πnM (3.13a)
( )
P(kW ) = 2πn rev T ( N .m ) × 10 −3
s
(3.13b)
Note que el par es la medida de la capacidad de un motor para realizar trabajo; la potencia es
la velocidad (o tasa) a la que se realiza el trabajo. La potencia de un motor medida como se
describió atrás es llamada potencia efectiva o al freno Pe.
1
El término indicado se usa porque tales diagramas p-V suelen ser generados directamente con un dispositivo
llamado un indicador de motores.
45
Capítulo 3. Parámetros de operación y de diseño de motores
Figura 3.4 Ejemplos de diagramas p-V para (a) motor de dos tiempos, (b) motor de cuatro tiempos; (c)
MEP de cuatro tiempos, carreras de admisión y escape (lazo de bombeo) a carga parcial. AE: Apertura
escape; AA: Apertura Admisión; AVA: Apertura válvula admisión; RCA: Retraso cierre admisión [1]
En motores de dos tiempos (Figura 3.4a), la aplicación de la ecuación (3.14) es directa. Con
la adición de carreras de admisión y escape para los cuatro tiempos, se hace necesario aclarar el
término de indicado en el ciclo:
Trabajo bruto indicado por ciclo Wc,ig. Trabajo entregado al pistón en las carreras de
compresión y expansión únicamente.
Trabajo neto indicado por ciclo Wc,in. Trabajo entregado al pistón en todo el ciclo incluyendo
el lazo de bombeo.
En las Figuras 3.4b y c, Wc,ig es (área A + C) y Wc,in es (área A + C) – (área B + C) lo que
equivale a (área A – B), donde cada una de estas áreas está considerada como una cantidad
positiva. El área A + C es el trabajo transferido entre el pistón y los gases del cilindro durante las
carreras de admisión y escape y se llama trabajo de bombeo Wp. La transferencia del trabajo de
bombeo será hacia los gases del cilindro si la presión durante la carrera de admisión es menor que
la presión durante la carrera de escape. Esta es la situación con motores de aspiración natural.
La transferencia del trabajo de bombeo será desde los gases del cilindro hacia el pistón si la
presión de la carrera de escape es más baja que la de la admisión, que es normalmente el caso con
los motores turboalimentados.
La potencia por cilindro está relacionada con el trabajo indicado por ciclo por:
Wc ,i n
Pi = (3.15)
i
donde i es el número de revoluciones del cigüeñal por cada carrera de potencia del cilindro. Para
cuatro tiempos, i es 2; para dos tiempos, i es 1. Esta potencia es la potencia indicada; es decir, la
tasa de transferencia de trabajo del gas dentro del cilindro hacia el pistón. Es un poco diferente
46
Capítulo 3. Parámetros de operación y de diseño de motores
de la potencia efectiva debido a la potencia absorbida en la fricción del motor, los accesorios
activos del motor, y (en el caso de la potencia bruta indicada) la potencia de bombeo.
En las cantidades indicadas de los motores de cuatro tiempos, tal como el trabajo o potencia
por ciclo, la definición usada para “indicado” (i.e. bruta o neta) se expresará siempre
explícitamente. El rendimiento bruto indicado, la definición más comúnmente usada, seguirá
siendo escogida a lo largo de este texto por las siguientes razones: las cantidades indicadas son
usadas primordialmente para identificar el impacto de los procesos de compresión, combustión y
expansión en el funcionamiento del motor, etc. El rendimiento bruto indicado es, por
consiguiente, la definición más apropiada. Ésta representa la cantidad de trabajo útil disponible
en el eje y el trabajo requerido para vencer todas las pérdidas del motor. Además las normas
técnicas para ensayo de motores definen procedimientos para medir la potencia al freno y la
potencia de fricción (el ensayo de potencia de fricción es una buena aproximación a la potencia
total de pérdidas en el motor). La suma de la potencia al freno y la de fricción proporciona una
forma alternativa de estimar la potencia indicada; el valor obtenido es una buena aproximación a
la potencia bruta indicada.
Los términos al freno (o efectiva) e indicado, son usados para describir otros parámetros tales
como la presión media efectiva (pme), el consumo específico de combustible (gef), y las
emisiones específicas en un modo similar al usado en el trabajo y potencia por ciclo.
RENDIMIENTO MECÁNICO
Hemos visto que parte del trabajo o potencia indicados por ciclo son usados para expulsar los
gases de escape e inducir carga fresca. Una porción adicional se usa para vencer la fricción de
los cojinetes, pistones y otros componentes mecánicos del motor y para impulsar sus accesorios.
Todos éstos requerimientos de potencia se agrupan y son llamados potencia de fricción Pf . De
éste modo:
Pig = Pe + Pf (3.16)
47
Capítulo 3. Parámetros de operación y de diseño de motores
presión del gas en el pistón y los anillos son más bajas en la prueba a motor arrastrado que con
combustión y las temperaturas del aceite en la pared del cilindro son también más bajas bajo
motor arrastrado.
La relación entre la potencia efectiva y la potencia indicada se llama rendimiento mecánico
ηm :
Pe Pf
ηm = = 1− (3.17)
Pig Pig
Puesto que la potencia de fricción incluye la potencia requerida para bombear el gas hacia
dentro y fuera del motor, el rendimiento mecánico depende de la posición de la mariposa así
como del diseño y velocidad del motor. Para motores de automoción, algunos valores típicos del
rendimiento mecánico a plena carga son cerca del 90% para regímenes en torno a 1800 – 2400
rpm, disminuyendo al 75% cuando el motor alcanza la máxima velocidad nominal. A medida
que se rebaja la carga en el motor, disminuye el rendimiento mecánico, llegando a valer cero
eventualmente en operación al ralentí.
P ⋅i
pme = (3.18a)
Vd ⋅ n
Para unidades SI,
P(kW ) ⋅ i × 10 −3
pme(kPa ) = (3.18b)
( )(
Vd dm 3 n rev
s
)
48
Capítulo 3. Parámetros de operación y de diseño de motores
La presión media efectiva puede ser expresada también en términos del par usando la
ecuación (3.13):
6.28 ⋅ i ⋅ M ( N .m )
pme(kPa ) = (3.19)
(
V d dm 3 )
La presión media efectiva máxima al freno ya está prefijada para un motor bien diseñado y es
prácticamente constante en un amplio rango de tamaños de motores. Así, una manera de medir la
capacidad del diseñador de un motor consiste en comparar la pme. También, cuando se realicen
cálculos de diseño en motores, se pueden estimar el desplazamiento del motor requerido para
proveer un par o potencia dados, a una velocidad especificada asumiendo valores apropiados para
la pme.
Para MEP de aspiración natural, los valores máximos de la pme están en el rango de 850 a
1050 kPa en la velocidad del motor en la que se obtiene el máximo par (cerca de 3000 rpm). En
la máxima potencia nominal, los valores de la pme son 10 a 15% más bajos. Para MEP
turboalimentados la máxima pme está en el rango de los 900 a 1400 kPa. Para MEC de
aspiración natural y cuatro tiempos, la máxima pme está en el rango de los 700 a 900 kPa, con la
pme en la máxima potencia nominal de cerca de 700 kPa. Los máximos valores de la pme en los
MEC de cuatro tiempos turboalimentados están generalmente en el rango de 1000 a 1200 kPa;
para los motores turboalimentados posenfriados ésta puede llegar a 1400 kPa. En la máxima
potencia nominal, la pme está alrededor de los 850 a 950 kPa. Los MEC de dos tiempos tienen
un funcionamiento similar a los motores de cuatro tiempos. Los motores grandes de baja
velocidad de dos tiempos pueden alcanzar valores de pme de cerca de 1600 kPa.
m f ⎛⎜ g ⎞⎟
⎛
gef ⎜ g ⎞ = ⎝ h⎠ (3.21)
⎟
⎝ kW .h ⎠ P(kW )
49
Capítulo 3. Parámetros de operación y de diseño de motores
Lo más deseable es alcanzar bajos valores del consumo específico de combustible. Para
MEP los mejores valores del consumo específico de combustible están cerca de los 270 g/kW-h.
Para motores de encendido por compresión, los mejores valores son más bajos y en los motores
más grandes pueden estar por debajo de los 200 g/kW-h.
El consumo específico de combustible tiene unidades, sin embargo existe un parámetro
adimensional que relaciona el rendimiento deseado del motor (trabajo o potencia por ciclo) con
la entrada necesaria (flujo de combustible), lo que daría un valor más representativo. La relación
entre el trabajo producido y la cantidad de energía del combustible suministrado que puede ser
liberada en el proceso de combustión se usa comúnmente para este propósito. Es una medida de
la eficiencia o rendimiento del motor. La energía del combustible suministrado que puede ser
liberada por combustión está dada por el producto entre el flujo másico de combustible y su
poder calorífico másico inferior. El poder calorífico de un combustible, HC , define su contenido
de energía. Es determinado en un procedimiento de ensayo estandarizado en el que una masa
conocida de combustible se quema completamente con aire, y la energía térmica liberada por el
proceso de combustión es absorbida por un calorímetro a medida que los productos de
combustión se enfrían a su temperatura original.
Esta medida del rendimiento de un motor, que en adelante se llamará rendimiento de
conversión de combustible ηf, 2 está dado por:
ηf =
Wc
=
(
P ⋅i
n
) =
P (3.22)
m f H C ⎛ m f ⋅ i ⎞ m f H C
⎜ n ⎟H C
⎝ ⎠
2
El término rendimiento de conversión de combustible se emplea porque es más descriptivo, y distingue mejor este
rendimiento de otras definiciones de rendimiento de un motor. Note que hay muchas definiciones diferentes del
poder calorífico. En éste texto, se usará el poder calorífico inferior para el cálculo del rendimiento.
50
Capítulo 3. Parámetros de operación y de diseño de motores
3600
ηf = (3.23b)
gef ⎛⎜ g ⎞⎟ H ⎛⎜ MJ ⎞⎟
⎝ kW .h ⎠ C⎝ kg ⎠
Los poderes caloríficos típicos para los combustibles comerciales usados en los motores están
en el rango de 42 a 44 MJ/kg. Así, el consumo específico de combustible es inversamente
proporcional al rendimiento de conversión de combustible. Note que la energía del combustible
suministrado al motor por ciclo no es completamente liberada como energía térmica en el proceso
de combustión porque el proceso real de combustión es incompleto. Cuando hay suficiente aire
presente en el cilindro para oxidar el combustible completamente, la mayoría (cerca del 96%) de
esta energía del combustible suministrado es transferida como energía térmica al fluido de
trabajo. Cuando no hay suficiente aire presente para oxidar el combustible completamente, la
falta de oxígeno evita liberar completamente la energía del combustible suministrado.
a y el flujo másico de
En el ensayo de motores, normalmente se miden el flujo másico de aire m
f . La relación entre estos flujos es útil al momento de definir las condiciones
combustible m
operativas del motor:
51
Capítulo 3. Parámetros de operación y de diseño de motores
rendimiento volumétrico ηv, como parámetro que indica la efectividad del proceso de llenado del
cilindro en el motor. Éste se usa solo con motores de cuatro tiempos los cuales tienen un proceso
distinto de inducción. Se define como la relación entre el flujo de aire que entra en condiciones
reales de operación del motor y el flujo de aire que entraría teóricamente a llenar el cilindro a las
condiciones del sitio:
2m a
ηv = (3.26a)
ρ a,iVd n
donde ρa,i es la densidad del aire admitido a las condiciones locales. Una definición alternativa
equivalente para el rendimiento volumétrico es:
ma
ηv = (3.26b)
ρ a,iVd
La densidad del aire en la admisión puede ser tomada como la del aire a las condiciones
ambiente del lugar (en tal caso, ηv mide la capacidad de bombeo de todo el sistema de admisión)
– es lo más usual - o puede ser tomada como la densidad del aire en el múltiple de admisión (en
tal caso, ηv mide la capacidad de bombeo del puerto de admisión y la válvula solamente) –
menos usual -. Los valores máximos típicos de ηv para motores de aspiración natural están en el
rango de 80 a 90%. El rendimiento volumétrico para los MEC es un poco más alto que para
motores los MEP como se verá en el capítulo 6.
Estos parámetros indican la efectividad con que el diseñador ha usado los materiales del
motor y sus componentes.
52
Capítulo 3. Parámetros de operación y de diseño de motores
(γ +1) / γ ⎫ 1/ 2
AE p o ⎧⎪ 2γ ⎡⎛ p ⎞
2/γ
⎛ p ⎞ ⎤⎪
m = ⎨ ⎢ ⎜ ⎟ − ⎜ ⎟ ⎥⎬ (3.29)
R g To ⎪γ − 1 ⎢⎣⎜⎝ p o ⎟⎠ ⎜p ⎟
⎝ ⎠ ⎥⎦ ⎪⎭
⎩ o
Al derivar esta ecuación, se ha asumido que el fluido es un gas ideal con Rg como constante
(
del gas y que la relación entre calores específicos C p / C v = γ ) es constante; po y To son la
presión y la temperatura total (estática + dinámica) aguas arriba de la restricción, p es la presión
en la restricción y AE es el área efectiva de flujo. Si en el motor p/po se asume como constante
para condiciones de plena carga o mariposa totalmente abierta, entonces para un sistema de
a varía según:
admisión y un motor dados, la razón de flujo másico de aire seco m
po
m a ∝ (3.30)
To
Pi ,s = C F Pi ,m (3.31)
53
Capítulo 3. Parámetros de operación y de diseño de motores
1/ 2
p s ,d⎛ Tm ⎞
CF = ⎜ ⎟
p m − pv ,m ⎜⎝ Ts ⎟⎠
(3.32)
La potencia efectiva se corrige usando la ecuación (3.32) para corregir así la potencia
indicada asumiendo que la potencia de fricción no cambia. Así:
Pe ,s = C F Pi ,m − Pf ,m (3.33)
η v,s ⎛ Ts ⎞
1/ 2
C F' = =⎜ ⎟ (3.34)
η v,m ⎜⎝ Tm ⎟⎠
54
Capítulo 3. Parámetros de operación y de diseño de motores
multiplicada por 106 o por 102, respectivamente). Las emisiones específicas son la relación del
flujo másico de contaminante por unidad de potencia de salida:
m NOx
sNO x = (3.35a)
P
m CO
sCO = (3.35b)
P
m HC
sHC = (3.35c)
P
m Part
sPart = (3.35d)
P
Al igual que en casos anteriores, se pueden definir las emisiones específicas efectivas o al
freno y las indicadas. Las unidades de uso común son, g / kW − h y g / hp − h .
Alternativamente a las emisiones específicas, se suelen usar los índices de emisiones de la
siguiente manera:
m NOx ( g / s )
EI NOx = (3.36)
m f (kg / s )
55
Capítulo 3. Parámetros de operación y de diseño de motores
η f ma nH C (F / A)
P= (3.37)
i
η f η v nVd H C ρ a ,i (F / A)
P= (3.38)
2
Para el par M:
η f η v V d H C ρ a ,i ( F / A)
M= (3.39)
4π
pme = η f η v H C ρ a ,i (F / A) (3.40)
La potencia por unidad de área del pistón, comúnmente llamada potencia específica, es una
medida del éxito del diseñador del motor en el uso del área disponible del pistón para cualquier
tamaño del cilindro. De la ecuación (3.38), la potencia específica es
P η f η v nLH C ρ a ,i (F / A )
= (3.41)
Ap 2
La velocidad lineal media del pistón puede ser introducida con la ecuación (3.9) para dar:
P η f η v S p H C ρ a,i (F / A)
= (3.42)
Ap 4
56
Capítulo 3. Parámetros de operación y de diseño de motores
Velocidad media del pistón. Mide el éxito comparativo en el manejo de cargas debido a la
inercia de las partes, resistencia al flujo de aire, y / o fricción del motor.
Presión media efectiva al freno. En motores de aspiración natural la pme no es limitada por
esfuerzos mecánicos. Ésta entonces refleja el producto del rendimiento volumétrico
(habilidad para inducir aire), la relación combustible / aire (efectividad de la utilización del
aire en la combustión), y el rendimiento de conversión de combustible. En motores
sobrealimentados, la pme indica el grado de éxito en el manejo de altas presiones de gas y
carga térmica.
Potencia por unidad de área del pistón. Mide la efectividad con que se usa el área del pistón, sin
tomar en cuenta el tamaño del cilindro.
Peso específico. Indica la economía relativa con la que se usan los materiales.
Volumen específico. Indica el rendimiento en la utilización del espacio en el motor
2. En todas las velocidades en las que el motor será usado a plena carga:
Presión media efectiva al freno. Mide la capacidad de obtener un alto flujo de aire y usarlo
efectivamente en todo el rango de operación del motor.
57
Capítulo 3. Parámetros de operación y de diseño de motores
Tabla 3.1
Datos típicos para el diseño y funcionamiento de los MEP [1]
Máximo Nominal
Ciclo de Relación de Diámetro, Carrera/ Velocidad, pme, Potencia Relación gef al freno
operación compresión m diámetro rev/min atm específica peso/Potencia aprox.
kW/dm3 kg/kW g/kW.h
MEP
Pequeños (e.g., 2S,4S 6-11 0.05- 1.3.0.9 4500- 4-10 20-60 5.5-3.5 350
motocicletas) 0.085 7500
Carros de 4S 8-10 0.07-0.1 1.1-0.9 4500- 7-10 20-50 4-2 270
pasajeros 6500
Camiones 4S 7-9 0.09-0.13 1.3.0.7 3600- 6.5-7 25-30 6.5-3.5 300
5000
Motores grandes 2S,4S 8-12 0.23.0.45 1.1-1.4 300-900 6.8-12 3-7 23-35 200
a gas
Motores Wankel 4S ≈9 0.57 dm3 por cámara 6000- 9.5- 35-45 1.6-0.9 300
8000 10.5
Tabla 3.2
Datos típicos para el diseño y funcionamiento de los MEC [1]
Máximo Nominal
Ciclo de Relación de Diámetro, Carrera/ Velocidad, pme, Potencia Relación gef al freno
operación compresión m diámetro rev/min atm específica peso/Potencia aprox.
kW/dm3 kg/kW g/kW.h
MEC
Carros de pasajeros 4S 17-23 0.075-0.1 1.3.0.9 4000- 5-7.5 18-22 5-3.5 250
5000
Camiones 4S 16-22 0.1-0.15 1.3-0.8 2100- 6-9 15-22 7-4 210
(Aspiración Natur.)
4000
Camiones 4S 14-20 0.1-0.15 1.3-0.8 2100- 13.18 18-26 7-3.5 200
(Turboalimentados.)
4000
Locomotoras, 4S,2S 13.18 0.15-0.4 1.1-1.3 425-1800 7-23 5-20 6-18 190
industriales,
marinos
Motores grandes, 2S 10-12 0.4-1 1.3.3 110-400 9-17 3.8 13.50 180
marinos y estáticos
58
Capítulo 3. Parámetros de operación y de diseño de motores
Los datos típicos de funcionamiento para los MEP y MEC en el rango de tamaños de
producción normal están resumidos en las Tablas 3.1 y 3.2 respectivamente. El ciclo de cuatro
tiempos predomina excepto en los rangos extremos de tamaño de motor. Los motores más
grandes son turboalimentados o supercargados. La máxima velocidad nominal del motor
disminuye a medida que el tamaño del motor aumenta, manteniendo la máxima velocidad lineal
media del pistón en el rango de 8 a 15 m/s. La máxima presión media efectiva al freno para
motores turboalimentados y supercargados es más alta que para los motores de aspiración natural.
Debido a que la máxima relación combustible - aire para MEP es más alta que para MEC, sus
niveles de pme de aspiración natural máximos son más altos. A medida que el tamaño del motor
aumenta, el consumo específico de combustible al freno disminuye y el rendimiento de
conversión de combustible aumenta, debido a la poca importancia de las pérdidas de calor y
fricción. Para los MEC más grandes, se pueden obtener rendimientos de conversión de
combustible efectivos cercanos del 50% y rendimientos de conversión de combustible indicado
de cerca de 55%.
Para un motor de cuatro cilindros, el volumen desplazado, diámetro y carrera se relacionan por
π
Vd = 4 × B2L
4
Asumiendo B = L ; esto da B = L = 86 mm
El régimen de giro a máxima potencia se puede estimar de un valor apropiado de la velocidad
lineal media del pistón, 15 m/s.
59
Capítulo 3. Parámetros de operación y de diseño de motores
La potencia máxima al freno se puede estimar de un valor típico de la pme a potencia máxima,
800 kPa, usando la ecuación (3.18):
N emax (kW ) =
( )
pme(kPa )V dm 3 nmax (rev / s ) 800 × 2 × 87
= = 70kW
i × 10 3 2 × 10 3
Ejercicio 3.2
De un motor diesel de automoción de 4 cilindros, aspiración natural y 4 tiempos funcionando en
Medellín (Condiciones ambiente promedio: 850 mbar y 20 ºC), se conocen los siguientes datos:
Se pide: 1. Carrera (L) y diámetro (B) del pistón; 2. Presión media efectiva; y 3. Potencia
efectiva; 4. Par efectivo.
Solución
1. Carrera y Diámetro
πB 2 πB 2
Vd = ⋅L⋅Z = 1,1B ⋅ 4 = 1,1πB 3 = 1906Cm 3
4 4
De donde B = 8,2 Cm y por tanto L = 9,02 Cm
2. Presión media efectiva:
pme(kPa ) = η v ⋅ ρ a ,i ⋅ F ⋅ η f ⋅ H c
1000 N
85kN ⋅
p amb kN kg
ρ a ,i = = =1 3
J N −m m
287 Tamb m 2 ⋅ 287 ⋅ (273 + 20) K
kg ⋅ K kg ⋅ K
1
gef =
η f ⋅ Hc
60
Capítulo 3. Parámetros de operación y de diseño de motores
1 1
ηf = = = 0,2873
g kJ 290 1h kJ
290 * 43200 * * 43200
kW − h kg 1000 3600 s kg
pme(kPa ) = η v ⋅ ρ a ,i ⋅ F ⋅ η f ⋅ H c = 0,78 ⋅ 1 ⋅
1
⋅ 0,2873 ⋅ 43200 = 537,8kPa ♦
18
3. Potencia efectiva:
pme(kPa ) ⋅ V d (m 3 ) ⋅ n(rps ) 537,8 ⋅ 1906 / 10 6 ⋅ 4000 / 60
N e (kW ) = = = 34kW ♦
i 2
4. Par efectivo:
Es importante en este ejercicio observar la variación de la densidad sobre las variables mecánicas
del motor, si por ejemplo tomáramos los datos ambientales promedio de Barranquilla, donde la
densidad media del aire está en torno a 1,2 kg/m3, tendríamos:
Ejercicio 3.3
El consumo específico de combustible de un MEP de cuatro tiempos en condiciones de plena
carga es 330 g/kW h. Sabiendo que la cilindrada total es 1100 cm3, el rendimiento mecánico es
0,81 y el poder calorífico de la gasolina es 44200 kJ/kg, calcular el rendimiento indicado, la
potencia efectiva y el par motor a 5500 rpm.
61
Capítulo 3. Parámetros de operación y de diseño de motores
ηe η 0,2468
ηm = ⇒ ηi = e = = 0,3 ♦
ηi ηm 0,81
2. Potencia efectiva
η v ⋅ V d ⋅ ρ a ,i ⋅ n 0,82 ⋅ 1100 ⋅ 10 −6 ⋅ 5000 kg
m a = = = 0,03758
2 2 ⋅ 60 s
1 kg
m f = ⋅ 0,03758 = 0,002087
18 s
kg kJ
Pe = η e m f Hc = 0,2468 ⋅ 0,002087 ⋅ 44200 = 22,8kW ♦
s kg
3. Par motor
kN − m
22,8
Pe s
Me = = = 43,5 N ♦
ϖ 5000 ⋅ 2π
60
REFERENCIAS
[1] Heywood, J. B. (1988), “Internal Combustion Engine Fundamentals”. McGraw Hill, New
York.
62
CAPÍTULO 4
64
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
pV = mRg T
Ru
pV = m T
Mg
(4.1)
pV = nRu T
p
= Rg T
ρg
65
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
⎛ ∂h ⎞ dh
Cp = ⎜ ⎟ = (4.5)
⎝ ∂T ⎠ p dT
h = u + pv
(4.6)
dh du ⎧ dv dp ⎫
= + ⎨p +v ⎬
dT dT ⎩ dT dT ⎭
pv = R g T
d ( pv )
= Rg (4.7)
dT
dv dp
p +v = Rg
dT dT
66
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
C p − Cv = R g (4.8)
Con frecuencia se realiza una hipótesis restrictiva adicional que consiste en tomar los calores
específicos como constantes, lo cual no es una parte necesaria en las relaciones de los gases
ideales.
En general, la energía interna y la entalpía de una gas ideal a una temperatura T relativa a su
energía interna y entalpía en alguna temperatura de referencia To, están dadas por:
u = u o + ∫ Cv (T )dT
T
(4.10)
To
h = ho + ∫ C p (T )dT
T
y (4.11)
To
La entropía en T, v, y p, relativa a la entropía en algún estado de referencia To, vo, po, se puede
obtener de las relaciones:
Cv dv C p dp
ds = dT + R g = dT − R g (4.12)
T v T p
T Cv v
s = so + ∫ dT + R g ln (4.13)
To T vo
T Cp p
y s = so + ∫ dT − R g ln (4.14)
To T po
67
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Las propiedades u, h, y s pueden ser evaluadas en una base másica o molar. En una base
másica, Cv, Cp, y Rg tendrían como unidades J/kg-K; en una base molar u, h, y s se reemplazan
~, h~ , y ~
por u s . Cv y Cp se reemplazan por
~ ~
Cv y C p y sus unidades al igual que las de Ru
serían J/kmol-K.
Presión parcial pi. Es la presión que cada componente ejercería si únicamente él ocupara el
volumen de la mezcla a la temperatura de la misma.
Partes por volumen Vi/V. Es la fracción del volumen de mezcla total que cada componente
ocuparía si fuera separado de la mezcla, a la temperatura y presión de ésta.
Fracción másica xi. Es la masa de cada componente mi, dividida por la masa total de mezcla
m.
Fracción molar ~
xi . Es el número de moles de cada componente ni, divido por el número
total de moles de la mezcla n.
pi Vi M ~
= = xi = xi (4.15)
p V Mi
Las propiedades termodinámicas de una mezcla de gases ideales se pueden calcular mediante
las siguientes relaciones:
1
M= ∑
n i
ni M i = ∑ ~
xi M i (4.16)
i
68
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
u = ∑ xi u i
i
h = ∑ xi hi (4.17)
i
s = ∑ xi s i
i
Aire
En los motores usualmente se quema el combustible con aire. El aire seco es una mezcla que
tiene una composición representativa por volumen de 20.95% de oxígeno, 78.09% de nitrógeno,
0.93% de argón, y algunas trazas de CO2, neón, helio, metano y otros gases. En la Tabla 4.1 se
muestran los principales constituyentes del aire seco.
Tabla 4.1
Principales Constituyentes del Aire Seco
Gas ppm por volumen Peso molecular Fracción molar Relación molar
O2 209500 31,998 0.2095 1
N2 780900 28,012 0.7905 3,773
Ar 9300 38,948
CO2 300 44,009
Aire 1000000 28,962 1 4,773
⎡ kg ⎤ 3,483 ⋅ p[kPa]
ρ aire ⎢ =
⎣ m 3 ⎥⎦
(4.18)
T [K ]
69
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Así pues, el valor para la densidad del aire “seco” a 1 atm (1013 mbar) y 25C (298K) es de
1,184kg/m3. Para Medellín, que en condiciones normales la presión atmosférica es de 852 mbar
(85,2 kPa) y para una temperatura media de 20C (293K), la densidad del aire seco es de 1,013
kg/m3.
En la realidad, y muy especialmente en Colombia, al ser un país tropical, el aire no es seco.
Normalmente contiene vapor de agua cuya cantidad depende de la temperatura y del grado de
saturación. La proporción en masa del contenido de vapor de agua en el aire puede variar entre 1
y 4 o 5 % en condiciones extremas. A continuación definiremos una serie de parámetros
necesarios para estudiar el aire “real”:
Humedad relativa: es la relación entre la presión parcial del vapor de agua presente en la
mezcla con el de la presión de saturación del agua a la temperatura de la mezcla.
pv
φ= (4.19)
pg
Temperatura de rocío (Drew point): Es aquella temperatura a la cuál empiezan a aparecer las
primeras gotas de condensado en una mezcla de gases (Tg).
Temperatura de bulbo seco (Dry - bulb Temp.): Es la temperatura del aire medida con un
termómetro corriente cuyo bulbo se expone al ambiente
Temperatura de bulbo húmedo (Wet – bulb Temp.): Esta se mide con un termómetro cuyo
bulbo cubierto con un algodón se humedece en agua. Es más baja que la de bulbo seco debido a
la evaporación del agua del algodón. Generalmente se asume que ésta es aproximadamente igual
a la temperatura de saturación adiabática.
Humedad absoluta o específica: Es la relación entre la masa de agua y la masa de aire se
expresa en kgagua/kgaire.
Normalmente cuando se trabaja con el aire “real”, el cual contiene vapor de agua, se emplea
la carta sicrométrica. Es importante resaltar que en la mayoría de los textos europeos y
estadounidenses se suele despreciar el efecto del contenido de vapor de agua en el aire, algunos
autores como Heywood (1988) [1] sugieren una corrección por contenido de vapor de agua
inferior a 0,03; sin embargo, en los principales asentamientos urbanos de nuestro país, se suelen
tener durante todo el año humedades relativas bastante elevadas (en Medellín por ejemplo,
siempre supera el 65%). Esto significa que se debe tener mayor consideración en todos los
cálculos termodinámicos que incluyan el aire como fluido de trabajo.
70
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Combustibles
La mayoría de los combustibles comúnmente usados en los MCIA (gasolina y acpm) son mezclas
de muchos hidrocarburos que se obtienen mediante la refinación del crudo de petróleo. Estos
combustibles son predominantemente carbono e hidrógeno (generalmente el 86% de C y el 14%
de H2 en masa), aunque algunos ACPM llegan a tener hasta un 1% en masa de azufre. Otros
combustibles de interés son los alcoholes (que contienen radicales oxigenados), combustibles
gaseosos (GNC y GLP), compuestos simples de hidrocarburos (metano, propano, iso-octano,
hidrógeno) que comúnmente se emplean en la investigación de motores. Las propiedades de los
combustibles más comunes en los MCIA se resumen en las Tablas 4.2 y 4.3.
Tabla 4.2
Entalpía de formación standard y pesos moleculares [1]
~
Especie Fórmula M (g/mol) Estado 1 Δh fo
MJ/kmol
Oxígeno O2 32 Gas 0
Nitrógeno N2 28,01 Gas 0
Carbón C 12,011 Sólido 0
Monóxido de Carbono CO 28,01 Gas -110,5
Dióxido de Carbono CO2 44,01 Gas -393,5
Hidrógeno H2 2,016 Gas 0
Gas -241,8
Agua H2O 18,02
Líquida -285,8
Metano CH4 16,04 Gas -74,9
Propano C3H8 44,1 Gas -103,8
Gas -224,1
Iso-octano C8H18 114,23
Líquido -259,28
Cetano C16H34 226,44 Líquido -454,5
Gas -201,2
Alcohol Metílico CH3OH 32,04
Líquido -238,6
Gas -234,6
Alcohol Etílico C2H5OH 46,07
Líquido -277
1
A 25C (298.15K) y 1atm (1013 mbar)
71
Tabla 4.3
Propiedades de algunos combustibles [1]
Calor específico Poder Octanaje
Calorífico
Gravedad Poder Poder
Calor de Inferior de
específica Vapor Calorífico Calorífico
Combustible Fórmula (estado) M vaporiz. Líquido la mezcla (A/F)est (F/A)est
(densidad, Cp Superior Inferior RON MON
kJ/kg kJ/kg-K estequiom
kg/m3) kJ/kg-K MJ/kg MJ/kg
étrica
MJ/kg
Combustibles corrientes
Gasolina CnH1,87n(l) ~110 0,72-0,78 350 2,4 ~1,7 47,3 44 2,83 14,6 0.0685 91-99 82-89
ACPM liviano CnH1,8n(l) ~170 0,78-0,84 270 2,2 ~1,7 46,1 43,2 2,79 14,5 0.069 --- ---
ACPM pesado CnH1,7n(l) ~200 0,82-0,88 230 1,9 ~1,7 45,5 42,8 2,85 14,4 0.0697 --- ---
Gas Natural CnH3,8nN0,1n(g) ~18 (~0,79) --- --- ~2 50 45 2,9 14,5 0.069 --- ---
Hidrocarburos puros
Metano CH4(g) 16,04 (0,72) 509 0,63 2,2 55,5 50 2,72 17,23 0.058 120 120
Propano C3H8(g) 44,1 0,51(2) 426 2,5 1,6 50,4 46,4 2,75 15,67 0.0638 112 97
Iso-octano C8H18(l) 114,23 0,692 308 2,1 1,63 47,8 44,3 2,75 15,13 0.0661 100 100
Cetano C16H34(l) 226,44 0,773 358 --- 1,6 47,3 44 2,78 14,82 0.0675 --- ---
Benceno C6H6(l) 78,11 0,879 433 1,72 1,1 41,9 40,2 2,82 13,27 0,0753 --- 115
Tolueno C7H8(l) 92,14 0,867 412 1,68 1,1 42,5 40,6 2,79 13,5 0,0741 120 109
Alcoholes
Metanol CH4O(l) 32,04 0,792 1103 2,6 1,72 22,7 20 2,68 6,47 0,155 106 92
Etanol C2H6O(l) 46,07 0,785 840 2,5 1,93 29,7 26,9 2,69 9 0,111 107 89
Otros combustibles
Carbono C(s) 12,01 ~2 --- --- --- 33,8 33,8 2,7 11,51 0,0869 --- ---
Monóx. de Carb CO(g) 28,01 (1,25) --- --- 1,05 10,1 10,1 2,91 2,467 0,405 --- ---
Hidrógeno H2 2,015 (0,09) --- --- 1,44 142 120 3,4 34,3 0,0292 --- ---
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
COMBUSTIÓN ESTEQUIOMÉTRICA
En esta sección se repasan algunas relaciones entre la composición de los reactivos y de los
productos. En una combustión completa ideal, los únicos productos de la combustión son CO2 y
agua. Por ejemplo, a continuación se muestra la ecuación química para la combustión completa
de una mol de propano:
⎛ b⎞ ⎛ b⎞
C a H b + ⎜ a + ⎟(O2 + 3,773N 2 ) = aCO2 + H 2 O + 3,773⎜ a + ⎟ N 2
b
(4.21)
⎝ 4⎠ 2 ⎝ 4⎠
−1
⎛ A⎞ 12 + y
Fest = ⎜ ⎟ = (4.22)
⎝ F ⎠ est 34,56(4 + y )
73
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
En algunas situaciones se puede llegar a quemar más o menos aire del requerido
estequiométricamente así, por ejemplo, para una combustión pobre en combustible, es decir, con
exceso de aire suele aparecer el oxígeno en los productos de la combustión. Por ejemplo, para la
combustión del iso-octano con un 25% de exceso de aire, es decir, 1,25 veces más del aire
requerido estequiométricamente se tendría:
Con menos aire del requerido estequiométricamente, es decir, una combustión rica en
combustible, no habría suficiente oxígeno para oxidar completamente el C y H del combustible a
CO2 y H2O. Bajo esta condición los productos serían una mezcla de CO2, H2O, CO, H2 y N2. En
este caso la composición de los productos no se puede determinar con un balance elemental
únicamente. Aquí es necesario realizar algunas hipótesis adicionales sobre la composición
química de los productos.
Debido entonces, a que la composición de los productos de la combustión es diferente para
mezclas pobres y ricas, y que además la relación combustible/aire (o F/A o más comúnmente
llamada dosado) depende de la composición del combustible, se considera que es un parámetro
más informativo el dosado relativo, que se define como la relación entre el dosado de
funcionamiento del motor y el dosado estequiométrico:
⎛F⎞
⎜ ⎟
⎝ A ⎠ abs Fabs
FR = = (4.24)
⎛F⎞ Fest
⎜ ⎟
⎝ A ⎠ est
Aquí se puede observar que el inverso del dosado relativo (FR) es lo que se denomina Lambda
(λ) o relación aire/combustible (A/F) relativa a la estequiométrica:
⎛ A⎞
⎜ ⎟
−1 ⎝ F ⎠ abs
λ = FR = (4.25)
⎛ A⎞
⎜ ⎟
⎝ F ⎠ est
74
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
QR − P − WR − P = U P − U R (4.26)
El agua en los productos de combustión puede estar en fase líquida o gaseosa. La energía
interna (o entalpía) de los productos en el proceso a volumen constante (o a presión constante)
descrito anteriormente, dependerá de la cantidad de agua y de su estado (líquido o gaseoso). En
la Figura 4.1 se muestran los casos extremos para cuando el agua en los productos de combustión
sale en estado gaseoso y en estado líquido. La diferencia de energía interna entre estas dos
curvas será:
(ΔU )V ,T ,H Oliq
2
− (ΔU )V ,T , H 2Ovap = m H 2O u fgH 2O (4.27)
donde mH2O será la masa de agua en los productos y ufgH2O es la energía interna de vaporización
del agua a la temperatura y presión de los productos.
75
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Figura 4.1 Energía interna en función de la Temperatura. Efecto del agua en los productos [1]
Entalpía de formación
Se define como entalpía de formación de un compuesto químico al incremento de entalpía
asociado con la reacción de formación de una mole de un compuesto dado a partir de sus
elementos, estando cada sustancia en su estado termodinámico estándar a la temperatura dada.
Comúnmente el estado estándar es de 1atm y 25ºC (298,15 K). Se asume que la entalpía de
formación para los elementos en este estado es 0. En la Tabla 4.2 se muestran algunos datos de
entalpía de formación para diferentes especies químicas.
Para una reacción de combustión dada, la entalpía de los productos en el estado standard está
dada por:
~
H Po = ∑ ni Δh fo,i (4.28)
productos
~
H Ro = ∑ ni Δh fo,i (4.29)
Re activos
76
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
H C p = −(ΔH ) p , 298 K
(4.30)
H CV = −(ΔU )V , 298 K
Combustión completa significa que todo el carbono se convierte en CO2 y todo el hidrógeno
se convierte en H2O y cualquier presencia de azufre se convertiría en SO2. El poder calorífico
usualmente se expresa en J/kg de combustible. No es necesario especificar cuanto oxidante se
mezcló con el combustible , ni tampoco si éste es aire u oxígeno puro.
En adelante emplearemos el término poder calorífico superior cuando el agua en los
productos de escape se condensa hasta el estado líquido (ha cedido su calor de vaporización). Y
el término poder calorífico inferior cuando el agua en los productos de escape sale en fase
gaseosa. Los dos poderes caloríficos a presión constante estarían relacionados mediante la
siguiente ecuación:
⎛ m H 2O ⎞
H C sup p = H C inf p + ⎜ ⎟h (4.31)
⎜ m f ⎟ fgH 2O
⎝ ⎠
(
donde m H 2O m f ) es la relación entre la masa de agua producida y la masa de combustible
quemado. Una expresión similar, la energía interna de vaporización del agua se aplica para los
poderes caloríficos superior e inferior a volumen constante. El poder calorífico a presión
constante es comúnmente el más usado, incluso se llega a omitir el término “a presión constante”.
77
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
aA + bB = cC + dD
ΔH R − P = [cH C + dH D − aH A − bH B ]T
⎛ ∂ΔH R − P ⎞
⎜ ⎟ = ΔC p R − P
⎝ ∂T ⎠ p
T T
∫298 ∂ΔH R − P = ∫
298
ΔC p R − P ∂T
T
ΔH T = ΔH 298 + ∫ ΔC p R − P dT
298
donde : ΔC p R − P = cC pC + dC pD − aC pA − bC pB (4.32)
78
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Tabla 4.4
Fluidos de trabajo según proceso y tipo de motor [1]
Proceso MEP MEC
Aire, combustible (líquido y vapor), Aire, gases recirculados (en algunos
Admisión gases recirculados (en algunos casos) casos para controlar los NOx) y gases
y gases residuales residuales
Aire, combustible (vapor), gases Aire, gases recirculados (en algunos
Compresión recirculados (en algunos casos) y casos para controlar los NOx) y gases
gases residuales residuales
Productos de combustión (mezclas de Productos de combustión (mezclas de
Expansión N2, H2O, CO2, CO, H2, O2, NO, OH, N2, H2O, CO2, CO, H2, O2, NO, OH,
O, H. . .) O, H. . .)
Productos de combustión
Productos de combustión
Escape (principalmente N2, H2O, CO2 y O2
(principalmente N2, H2O, CO2 y O2)
si FR < 1 ó CO y H2 si FR > 1)
Tabla 4.5
Modelos para el cálculo de las propiedades termodinámicas [1]
Mezcla sin quemar Mezcla quemada
1. Un solo gas ideal para todo el ciclo con Cv y Cp constantes
2. Gas ideal; Cv,u constante Gas ideal; Cv,b constante
3. Mezcla de gases ideales, Cv,i (T) Mezcla de gases ideales, Cv,i (T)
4. Mezcla de gases ideales, Cv,i (T) Aproximación al equilibrio químico
Mezcla de gases ideales que reaccionan en
5. Mezcla de gases ideales, Cv,i (T)
equilibrio químico
Nota: El subíndice i indica la especie i en la mezcla de gases, u indica mezcla sin quemar y b indica
mezcla quemada.
79
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
La primera categoría únicamente se usa con propósitos ilustrativos debido a que los calores
específicos para la mezcla quemada y sin quemar son muy diferentes. La segunda categoría es
una aproximación analítica simple de poca precisión debido a que los calores específicos varían
fuertemente con la temperatura.
Los valores de las propiedades termodinámicas de las mezclas quemadas y sin quemar
relevantes para los cálculos en motores se pueden obtener de cartas, tablas y relaciones
algebraicas obtenidas mediante regresiones de datos tabulados.
CÁLCULO DE LAS PROPIEDADES
Las cartas termodinámicas son un método de presentar las propiedades termodinámicas de
mezclas de gases quemadas y sin quemar para cálculos en MCIA. Existen dos grupos de cartas
de uso común: las desarrolladas por Hottel y otros [5] y las desarrolladas por Newhall y
Starkman [6]. Ambas utilizan unidades del sistema Británico. Las cartas de propiedades
termodinámicas se usan poco actualmente; ya que están siendo reemplazadas por modelos de
computador para el cálculo de las propiedades de los fluidos de trabajo en motores. Sin embargo,
éstas son útiles para propósitos ilustrativos y representan un método fácil y preciso cuando se
requiere un número limitado de cálculos.
Cuando se quieren realizar cálculos bastante precisos para determinar la composición y las
propiedades de los gases sin quemar durante todo el ciclo de operación del motor se pueden
emplear rutinas para computador. Ésta varían considerablemente en rango de aplicación y en
precisión. Los modelos más completos se basan en ajustes polinomiales de curvas de datos
termodinámicos para cada especie en la mezcla con las siguientes hipótesis: (1) la mezcla sin
quemar no varía de composición y (2) que la mezcla quemada está en equilibrio químico. En éste
texto usaremos las tablas de datos termodinámicos de JANAF.
Para cada especie i en su estado estándar y a la temperatura T (K) se cumple:
c~p ,i
= ai1 + ai 2T + ai 3T 2 + ai 4T 3 + ai 5T 4 (4.33)
Ru
~
hi a a a a a
= ai1 + i 2 T + i 3 T 2 + i 4 T 3 + i 5 T 4 + i 6 (4.34)
Ru T 2 3 4 5 T
80
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
~
si a a a
= ai1 ln T + ai 2T + i 3 T 2 + i 4 T 3 + i 5 T 4 + ai 7 (4.35)
Ru 2 3 4
Los valores de los coeficientes aij para CO2, H2O, CO, H2, O2, N2, OH, NO, O y H del
programa de la NASA se dan en la Tabla 4.6. Se tienen en dicha tabla dos rangos de
temperaturas. El rango de 300 a 1000 K es apropiados para el cálculo de las propiedades de las
mezclas sin quemar. El rango de 1000 a 5000 K es apropiado para el cálculo de las propiedades
de la mezcla quemada.
MEZCLAS SIN QUEMAR
Se han ajustado funciones polinómicas para varios combustibles (en estado gaseoso), una de ellas
responde a la siguiente ecuación:
Af 5
c~p , f = A f 1 + A f 2 t + A f 3t 2 + A f 4 t 3 + 2 (4.36)
t
~ t2 t3 t 4 Af 5
h f = Af 1 + Af 2 + Af 3 + Af 4 − + Af 6 + Af 8 (4.37)
2 3 4 t
donde t = T(K)/1000. Af6 es la constante para el dato de cero entalpía para C, H2, O2 y N2 a
298.15 K. Para un dato de 0 K, Af8 se suma a Af6. Para hidrocarburos puros, los coeficientes Afi
se encontraron ajustando las ecuaciones (4.36) y (4.37) a los datos de Rossini y otros [1]. Los
valores para los combustibles puros importantes se dan en la Tabla 4.7. Las unidades para
~
c~p , f son cal/gmol-K, y para h f son kcal/gmol.
Así pues, se pueden ahora obtener las propiedades termodinámicas de la mezcla no quemada
conociendo: las moles de cada especie por mol de O2 y ni, determinada de la Tabla 4.8; la masa
de mezcla por mol de O2 y mR-P, determinada de la Tabla 4.9 y, empleando las ecuaciones (4.38
a 4.40).
81
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Tabla 4.6
Coeficientes para las propiedades termodinámicas de las especies. Fuente Heywood [1]
Especie Rango de T (K) ai1 ai2 ai3 ai4 ai5 ai6 ai7
1000-5000 0,44608(+1) 0,30982(-2) -0,12393(-5) 0,22741(-9) -0,15526(-13) -0,48961(+5) -0,98636(0)
CO2
300-100 0,24008(+1) 0,87351(-2) -0,66071(-5) 0,20022(-8) 0,63274(-15) -0,48378(+5) 0,96951(+1)
Tabla 4.7
Coeficientes de los polinomios (ecuaciones (4.36) y (4.37)) para la entalpía del combustible y el calor específico [1]
Combus
Fórmula M (A/F)est (F/A)est Af1 Af2 Af3 Af4 Af5 Af6 Af8
tible
Metano CH4 16,04 17,23 0,0580 -0,29149 26,327 -10,610 1,5656 0,16573 -18,331 4,3000
Propano C3H8 44,10 15,67 0,0638 -1,4867 74,339 -39,065 8,0543 0,01219 -27,313 8,852
Hexano C6H14 86,18 15,24 0,0656 -20,777 210,48 -164,125 52,832 0,56635 -39,836 15,611
Iso-octano C8H18 114,2 15,14 0,0661 -0,55313 181,62 -97,787 20,402 -0,03095 -60,751 20,232
Metanol CH3OH 32,04 6,47 0,1546 -2,7059 44,168 -27,501 7,2193 0,20299 -48,288 5,3375
Etanol C2H5OH 46,07 9,00 0,1111 6,990 39,741 -11,926 0 0 -60,214 7,6135
C8,26H15,5 114,8 14,64 0,0683 -24,078 256,63 -201,68 64,750 0,5808 -27,562 17,792
Gasolina
C7,76H13,1 106,4 14,37 0,0696 -22,501 227,99 -177,26 56,048 0,4845 -17,578 15,235
ACPM C10,8H18,7 148,6 14,4 0,0694 -9,1063 246,97 -143,74 32,329 0,0518 -50,128 23,514
Nota:
~ ~
Las unidades de A fi al igual que las de h f están en kcal/gmol y las de C p , f están en cal/gmol-K con t = T(K)/1000
( )
A f 6 da el dato de la entalpía a 298,15K; A f 6 + A f 8 da el dato de la entalpía a 0 K.
83
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
1
c p ,u =
mR−P
∑ ni c~p,i (4.38)
i
1 ~
hu =
mR−P
∑ ni hi (4.39)
i
1 ⎧⎪ ⎡~ o ⎛ ni ⎞⎤ ⎫⎪
su = ⎨∑ i ⎢ si − Ru ln⎜⎜
n ⎟⎟⎥ − nu Ru ln p ⎬ (4.40)
mR− P ⎪⎩ ⎣ ⎝ nu ⎠⎦ ⎪⎭
donde p se da en atmósferas.
Las Figuras 4.2 y 4.3, obtenidas con las relaciones anteriores, muestran como Cp,u y
γ u (= c p ,u cv,u ) varían con la temperatura, el dosado relativo (FR), y la fracción de gases
quemados (xb), para una mezcla gasolina – aire.
Tabla 4.8
Composición de la mezcla sin quemar [1]
ni, moles / mol O2 reactivo
Especie
FR ≤ 1 FR > 1
Combustible 4(1 − xb )(1 + 2ε )FR / M f
O2 1 − xb FR 1 − xb
N2 ψ ψ
CO2 xb εFR xb (εFR − c )
H2O 2 xb (1 − ε )FR xb [2(1 − εFR ) + c ]
CO 0 xb c
H2 0 xb [2(FR − 1) − c ]
Sum nu nu
84
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Tabla 4.9
Factores para relacionar las propiedades en bases molar y másica [1]
Cantidad por mol de O2 en la Ecuación para la mezcla
Ecuación general
mezcla Iso-octano (C8H18 ) – aire
nb = (1 − ε )FR + 1 + ψ , FR ≤ 1 → nb = 0.36 FR + 4.773
nb = (2 − ε )FR + ψ ,
Moles de mezcla quemada nb
FR > 1 → nb = 1.36 FR + 3.773
⎧ FR ≤ 1 → nu = 0.08 FR + 4.773 + 0.28 xb FR
⎨
⎩ FR > 1 → nu = 0.08 FR + 4.773 + xb (1.28 FR − 1)
Moles de mezcla sin quemar nu
⎡ 4(1 − 2ε )FR ⎤
nu = (1 − xb )⎢ + 1 + ψ ⎥ + xb nb
⎢⎣ Mf ⎥⎦
85
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
86
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
MEZCLAS QUEMADAS
La aproximación más precisa para el cálculo de las propiedades y composición de la mezcla
quemada consiste usar un programa de equilibrio termodinámico a temperaturas por encima de
los 1700 K y considerar que la composición de la mezcla no varía por debajo de 1700 K. Las
propiedades de cada especie a altas y bajas temperaturas se dan mediante funciones polinómicas
como las ecuaciones (4.33) a (4.35) y sus coeficientes en la Tabla 4.6.
El método más empleado consiste en calcular las propiedades y/o composición de los
productos suponiendo que no existe disociación y luego emplear métodos iterativos o
correcciones que permitan tener en cuenta los efectos de la disociación.
Olikara y Borman [7], han desarrollado un programa de computador para calcular las
propiedades de equilibrio de los productos de combustión. Éste ha sido diseñado específicamente
para el uso en aplicaciones de los MCIA. Se especifican en éste la composición del combustible
(CnHmOlNk), el dosado relativo, y la presión y temperatura de los productos. Las especies
incluidas en la mezcla de productos son: CO2, H2O, CO, H2, O2, N2, Ar, NO, OH, O, H y N. Las
ecuaciones de balance de los elementos y las constantes de equilibrio para siete reacciones no
redundantes dan como resultado un total de 11 ecuaciones requeridas para solucionar la
concentración de estas especies. Las constantes de equilibrio son curvas ajustadas a partir de los
datos de las tablas de JANAF. La estimación inicial de las fracciones molares para comenzar el
proceso iterativo es la composición sin disociar. Una vez se determina la composición de la
mezcla, se calculan las propiedades termodinámicas y sus derivadas con respecto a la
temperatura, presión y dosado relativo. Se ha encontrado que este conjunto limitado de especies
es lo suficientemente preciso para el cálculo de los gases quemados de motores, y es mucho más
rápido que el extenso programa de equilibrio de la NASA.
Se han desarrollado numerosas técnicas para estimar las propiedades termodinámicas de los
gases quemados a altas temperaturas en aplicaciones de motores. Una aproximación
comúnmente usada es la desarrollada por Krieger y Borman [8]. La energía interna y la constante
de los gases de los productos de combustión sin disociar aparecen como polinomios en función
de la temperatura del gas. El segundo paso consistió en limitar el rango de T y p a valores
encontrados en MCI. A continuación las desviaciones entre los datos de las propiedades de
equilibrio termodinámico publicados por Newhall y Starkman y los valores sin disociar
calculados, se ajustaron mediante una función exponencial de T, p y FR. Para FR ≤ 1 , resulta un
solo conjunto de ecuaciones. Para FR ≥ 1 , se desarrollaron varios conjuntos de ecuaciones, cada
uno aplicado a un valor específico del dosado relativo (ver la referencia [8]). En general, el
87
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
ajuste para la energía interna está en el rango de 2 1/2 % por encima de la presión y temperatura
de interés, y el error en gran parte del rango es menor al 1%. Para muchas aplicaciones, las
ecuaciones sin disociar para las propiedades termodinámicas son lo suficientemente precisas.
Ejercicio 4.1
Un motor diesel de inyección directa turboalimentado, de cuatro tiempos con seis cilindros en
línea, tiene un diámetro de cilindro de 112 mm y una carrera de 130 mm, y se ha ensayado en
banco de prueba bajo condiciones de funcionamiento estacionarias (ver mediciones en tabla
anexa). Las condiciones ambientales de la celda de ensayos son presión 852 mbar y temperatura
20 ºC.
El combustible que se ha utilizado es un ACPM cuya fórmula química es C10.8H18.7, con una
densidad de 0,82 kg/dm3. Suponiendo que todo el combustible se quema estequiométricamente,
se pide:
1. Dosado estequiométrico
2. Poder calorífico másico inferior en [kJ/kgcombustible]
3. Gasto másico de combustible en [g/s]
4. Masa de combustible inyectada en [g/(ciclo*cilindro)]
5. Gasto másico de aire en [kg/s]
6. Exceso o defecto de aire en [%]
7. Potencia efectiva del motor en [kW]
8. Presión media efectiva del motor en [bar]
9. Consumo específico de combustible en [g/kW-h]
10. Rendimiento efectivo del motor
Datos:
MAGNITUD UNIDAD MEDICIÓN
Régimen de giro del motor min-1 2000
Fuerza en el freno N 700
Tiempo consumo 1 litro (1dm3) de combustible s 93,28
3
Flujo volumétrico de aire m /min 12
Brazo del freno electromagnético mm 955
ENTALPÍAS DE FORMACIÓN[Δhf en MJ/kmol]
N2 y O2 0
CO2 -393,52
H2Ogas -241,83
H2Oliq -285,84
C10.8H18.7 -104,55
88
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Solución
C10.8 H 18.7 + 15,475(O2 + 3,773N 2 ) → 10.8CO2 + 9,35H 2O + 46,4 N 2
1) Dosado estequiométrico
mamb 12 *10,8 + 1*18,7
f stq = = = 0,0693 ♦
maire 15,475 * 4,773 * 28,962
H c = − ΔH = (− 104,55 + 6511,12 )
MJ 1kmolcomb 1000kJ
* *
kmolcomb 148,3kg 1MJ
kJ
H c = 43200 ♦
kg comb
3) Gasto másico de combustible
1dm 3 kg 1000 g
m f = * 0,82 3 * = 8,79 g ♦
93,28seg dm 1kg s
4) Masa de combustible inyectada por ciclo y por cilindro
g rev 1kg
8,79 * 2 *
m f ⋅ i s cic 1000 g kg
mf = = = 8,79 × 10 −5 ♦
n ⋅ Z 2000 rev * 1min * 6cilind . cilind ⋅ cic
min 60 s
5) Gasto másico de aire
⎛ p T ⎞ kg m 3 1min
m a = ρ a' ⋅ Va = ⎜⎜ ρ st * m * st ⎟⎟ 3 *12 *
⎝ p st Tm ⎠m min 60 s
♦
⎛ 850 (273 + 15,6 ) ⎞ 12 12,12
= ⎜⎜1,293 * * ⎟⎟ * = ⇒ m a = 0,2136 kg
⎝ 1013 (273 + 20 ) ⎠ 60 60 s
6) Exceso de aire
m f 8,79
Freal = = = 0,041
m a 213,6
89
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Freal 0,041
FR = = = 0,59
Fstq 0,0693
1
λ = FR−1 = = 1,69 ⇒ 69% ♦
0,59
7) Potencia efectiva
Ne = M e ⋅ n
M e = F ⋅ b = 700 N * 0,955m = 668,5 N ⋅ m
rev 2πrad 1min 1kN
N e = 668,5 N ⋅ m * 2000 * * * = 140kW ♦
min 1rev 60 s 1000 N
8) Presión media efectiva
2π ⋅ i ⋅ M e ( N .m )
pme = = 1093kPa[10,9bar ] ♦
(
Vd dm 3 )
π
donde Vd = D p2 * S * Z = 7,685dm 3
4
9) Consumo específico de combustible
m f 8,79 g * 3600s g
gef = = = 226 ♦
N e seg *1h *140kW kW ⋅ h
10) Rendimiento efectivo
1 kW ⋅ h *1kJ * 3600 s *1000 g * kg
ηe = = = 0,368 ♦
gef ⋅ H c 226 g *1kW * seg *1h *1kg * 43200kJ
Ejercicio 4.2
Un MEP de cuatro tiempos y 2000 cm3 de cilindrada funciona con etanol, realizándose la
combustión de forma estequiométrica. En las condiciones nominales de funcionamiento a 4300
rpm, la masa de aire aspirado es de 60 g/s siendo la presión y temperatura ambiente 0,9 bar y 3C,
respectivamente. Se supone que el nitrógeno del aire no reacciona en la combustión y que la
composición del aire es de 79% nitrógeno y 21% oxígeno (en volumen). Se pide:
1. Rendimiento volumétrico
2. Dosado de funcionamiento
3. Potencia efectiva nominal del motor
4. Par motor
90
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
Datos: Poder calorífico del etanol 26600 kJ/kg. Masa atómicas: Oxígeno 16, Nitrógeno 14,
Hidrógeno 1, Carbono 12. Rendimiento indicado 44%. Rendimiento mecánico 89%. Densidad
del aire a 25 ºC y 1013 mbar , 1,18 kg/m3.
Solución
1) Rendimiento volumétrico
poVo p1V1 pV T m poT1 m
= → V1 = o o ⋅ 1 → = ⋅
To T1 To p1 ρ1 To p1 ρ o
ρ oTo p1 1,18 * 298 * 0,9
ρ1 = = ⇒ ρ1 = 1,13 kg 3
poT1 276 *1,013 m
m a 0,06 0,06 kg / s
ηv = = = = 0,74 ♦
m T V ⋅ n ⋅ ρ
1
× −6 4300 1 m / s * kg / m 3
3
T 2000 10 * * * 1,13
2 60 2
2) Dosado de funcionamiento
C 2 H 6O + (O2 + 3,76 N 2 ) → CO2 + H 2O
C 2 H 6O + 3(O2 + 3,76 N 2 ) → 2CO2 + 3H 2O + 11,28 N 2
Proporciones estequiométricas
46 g cvomb / mol → 411,84 g aire / mol
m f 411,84 ♦
Fstq = = = 8,95
m a 46
3) Potencia efectiva
1 kg kJ
Pe = m f ⋅ H c ⋅ η e = F ⋅ m a ⋅ H c ⋅ η e = ⋅ 0,06 ⋅ 26600 ⋅ η e
8,95 s kg
donde: η e = η i ⋅ η m = 0,44 ⋅ 0,89 = 0,3916
1
Ne = ⋅ 0,066 ⋅ 26600 ⋅ 0,3916 ≈ 70kW ♦
8,95
4) Par efectivo
M=
Ne
=
70000
(W ⋅ s ) = 155 N ⋅ m ♦
ω 2π
n
60
Ejercicio 4.3
De un motor diesel turboalimentado de inyección directa, con seis cilindros en línea y 7,685 litros
de cilindrada, se tienen los siguientes datos de funcionamiento:
91
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
n = 2004rpm
kg g
ma = 0,00165 = 1,65
ciclo ⋅ cil ciclo ⋅ cil
g
m f = 0,0424618
ciclo ⋅ cil
FR = 0,35
Encontrar la composición de los gases quemados cuando su temperatura es inferior a 1700K.
Solución
Suponiendo que la composición química del ACPM corresponde a C10.8H18.7 (de Tabla 4.7),
tenemos:
4
ε=
4+ y
H 18,7
donde : y = = = 1,7315
C 10,8
⇒ ε = 0,6979
moles ⋅ de ⋅ N 2
Además, ψ = 3,773
moles ⋅ de ⋅ O2
De la Tabla 4.8 se tiene:
mol. ⋅ especie
Especie FR ≤ 1 mol. ⋅ O2 ⋅ real
CO2 (0,6979)(0,35) 0,24426
H2O 2(1 − 0,6979 )(0,35) 0,21147
O2 1 – 0,35 0,65
N2 3,773 3,773
∑ moles = n b
(1 − 0,6979)0,35 + 1 + 3,773 4,878735
92
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
moles ⋅ O2 ⋅ reactante
moles ⋅ O2 ⋅ reactante = 0,010828
ciclo ⋅ cil
2. Ahora podemos conocer el número de moles de cada especie en el cilindro (quemados)
multiplicando el resultado anterior por la tercera columna de la tabla anterior:
Especie FR ≤ 1 ηi moles de la especie
CO2 (0,24426)(0,010828) 0,002645
H2O (0,21147 )(0,010828) 0,0022898
O2 (0,65)(0,010828) 0,0070382
N2 (3,773)(0,10828) 0,040854
∑ moles = n b
(4,878735)(0,10828) 0,0528269
3. Obtenemos ahora la fracción de moles de la especie con respecto a las moles totales: ♦
Especie FR ≤ 1 Fracción (%)
CO2 0,002645 / 0,052827 5
H2O 0,0022898 / 0,05827 4,33
O2 0,0070382 / 0,052827 13,32
N2 0,040854 / 0,052827 77,35
93
Capítulo 4. Termoquímica de las mezclas aire – combustible y propiedades de los fluidos de trabajo
REFERENCIAS
[1] Heywood, J.B. (1988), “Internal Combustion Engine Fundamentals”. McGraw Hill, New York.
[2] Spalding, D. Brian. (1979), “Combustión y transferencia de masa”. CECSA, Mexico.
[3] Glassman I. (1977), “Combustion”. Academic Press, New York
[4] Mölenkamp, H. (1976), “Zur Genarigkeit Der Brenngezetzrechnung eines Dieselmotors mit
Nichtunterteiltem Brennraum” (“Sobre la exactitud del cálculo de la ley de combustión de un motor de
diesel con cámara de combustión no dividida”), MTZ Vol. 37 7/8
[5] Hottel, H.C., Williams, G.C., and Satterfield, C.N. (1949), “Thermodinamic charts for combustion
processes”, John Wiley
[6] Newhall, H.K., and Starkman, E.S. (1964), “Thermodinamic properties of octane and air for engine
performance calculations” in Digital calculations of engine cycles, Progress in technology, vol. TP-7, pp.
38-48, SAE.
[7] Olikara, C and Borman, G. L. (1975), “A Computer Program for Calculating Propperties of
Equilibrium Combustion Products with some applications to IC Engines”. SAE paper No. 750468
[8] Krieger, R. B., and Borman, G. L. (1966), “The Computation of Apparent Heat Release for
Internal Combustion Engines” in Proc. Diesel Gas Power, ASME paper 66-WA/DGP-4.
94
CAPÍTULO 5
Un ciclo como el anterior se conoce con el nombre de “Ciclo de Aire ideal”; y según los
procesos del ciclo, se habla entonces, del ciclo ideal de aire OTTO, ciclo ideal de aire DIESEL, y
otros como el de Lenoir, Brayton y Atkinson. El análisis detallado de estos ciclos ideales, nos
indicará cómo es el comportamiento entre el rendimiento del ciclo dependiendo de características
como la relación de compresión y de expansión, además, la forma como varía el rendimiento con
la cantidad de calor aportado, o si permanece constante en este caso, según el ciclo, y también su
relación con el método empleado para la aportación y rechazo de calor, etc.
Esta aproximación no permite cuantificar de manera precisa las variables del ciclo, debido a
que en la realidad ocurren una serie de transformaciones físicas y químicas en el fluido de trabajo
cuya actuación conjunta es difícil de determinarse mediante métodos analíticos meramente. No
Capítulo 5. Ciclos de trabajo teóricos y reales
obstante es posible acercar el análisis en grado considerable a la situación real mediante cálculos
más refinados. Esto implicaría tener en consideración elementos como:
96
Capítulo 5. Ciclos de trabajo teóricos y reales
Combinando todos estos posibles procesos podemos clasificar los ciclos en dos categorías:
97
Capítulo 5. Ciclos de trabajo teóricos y reales
Para comparar el rendimiento de los ciclos (Figura 5.6) se puede seguir un orden similar a los
expuestos anteriormente, teniendo en cuenta que:
• El proceso de aportación de calor a volumen constante, al permitir una mayor expansión
es más eficiente que el isobárico en términos generales.
• El enfriamiento a presión constante también es más eficiente que el isométrico, pues
permite igualmente una mayor expansión.
• El ciclo que permita una mayor compresión también será más eficiente, pues también
aumenta la expansión.
98
Capítulo 5. Ciclos de trabajo teóricos y reales
99
Capítulo 5. Ciclos de trabajo teóricos y reales
(el cual, desde que el rendimiento de la combustión sea 1, equivale al rendimiento térmico
indicado de conversión de combustible ηt ,i )
La presión media indicada (pmi):
Wc ,i m f H Cinf η f ,i
pmi = = (5.2)
Vd Vd
Wc,i, el trabajo indicado por ciclo, es la suma de los trabajos de las carreras de compresión y
expansión:
Wc ,i = WC + WE (5.3)
100
Capítulo 5. Ciclos de trabajo teóricos y reales
Carrera de compresión:
v1
= rc (5.4)
v2
Mientras el proceso sea adiabático y reversible:
s 2 = s1 (5.5)
El trabajo de compresión es:
WC = U 1 − U 2 = m(u1 − u 2 ) (5.6)
Proceso de combustión:
Para el ciclo a volumen constante:
v3 = v2 u3 − u 2 = 0 (5.7a, b)
Para el ciclo a presión constante:
p3 = p 2 h3 − h2 = 0 (5.7c, d)
Para el ciclo a presión limitada:
v3 a = v 2 p3b = p3a (5.7e, f)
Carrera de expansión:
Para el ciclo a volumen constante:
v4
= rc s 4 = s3 (5.8a, b)
v3
y el trabajo de expansión es:
WE = U 3 − U 4 = m(u3 − u 4 ) (5.9)
Para el ciclo a presión constante:
v4
p3 = p 2 = rc s 4 = s3 (5.10a, b, c)
v2
101
Capítulo 5. Ciclos de trabajo teóricos y reales
CICLOS REALES
Los ciclos de trabajo ideales, como se dijo antes, difieren bastante de la realidad, las principales
razones según cada proceso dentro del motor son:
Compresión
El proceso de compresión real es prácticamente igual cuando se analiza teóricamente, es decir,
cuando se analiza el proceso adiabático y reversible. Se dice que es adiabático porque la
diferencia de temperaturas entre el gas y la pared es pequeña (teóricamente). Y se dice que es
reversible porque las velocidades a las que se realiza el proceso son relativamente bajas y el
proceso se puede considerar mecánicamente reversible.
102
Capítulo 5. Ciclos de trabajo teóricos y reales
Combustión
La combustión real no es igual a la ideal por diversas razones. En el diagrama p – V se pueden
diferenciar fácilmente las pérdidas más relevantes en el proceso de combustión real respecto al
teórico:
1. Pérdidas de tiempo.
Tienen lugar debido a que en la realidad el proceso no es instantáneo, es necesario dar un
cierto tiempo para que se produzca la combustión. Esto se refleja en una pérdida de área en el
diagrama.
2. Pérdidas de calor.
En la práctica hay que refrigerar los motores por varias razones:
• Porque los materiales no resisten temperaturas tan elevadas
• Para mantener la temperatura del aceite en un valor tal que éste no se degrade
La refrigeración implica una disminución del rendimiento del motor, ya que parte del calor
que podría transformarse en trabajo útil se pierde a través de las paredes del volumen de
control.
103
Capítulo 5. Ciclos de trabajo teóricos y reales
Ejercicio 5.1
Un ciclo de aire equivalente de presión limitada, tiene en el instante inicial de la compresión una
presión p1 = 0,9 bar y una temperatura T1 = 41 ºC, su relación de compresión volumétrica es rc =
14, la presión máxima de combustión es p3 = 70 bar y el calor total aportado al ciclo es q = 2000
kJ/kg.
Calcular el rendimiento del ciclo, el trabajo específico y la pmi.
Datos:
p
3 3A kJ
q 2−3 A = 2000
70 kg
kJ
Cp =1
kg ⋅ K
2 4 J
Ra = 287
kg ⋅ K
1
0,9 γ aire = 1,4
rc = 14
V
V2 = V3 V1 = V4
104
Capítulo 5. Ciclos de trabajo teóricos y reales
punto 1
p1 = 0,9bar
T1 = 41C = 314 K
Rg T1 287 J 314 K 1bar ⋅ m 2 N ⋅ m m3
V1 = = * * * = 1
p1 kg ⋅ K 0,9bar 105 N J kg
V1 + V2 V1 1 m3
rc = 14 = ⇒ V2 = = = 0,077 = V3
V2 rc − 1 13 kg
punto 2
3
V2 = 0,077 m
kg
pV γ = cte ⇒ p1V1γ = p2V2γ
γ
⎛V ⎞
1.4
⎛ 1 ⎞
p2 = p1 ⎜⎜ 1 ⎟⎟ = 0,9bar ⎜ ⎟ = 32,6bar
⎝ V2 ⎠ ⎝ 0,077 ⎠
pV 32,6bar * 0,077 m 3 ⋅ kg ⋅ K ⋅ 105 N
T2 = 2 2 = = 874,5K
Rg 287 N ⋅ m ⋅ kg ⋅ 1bar ⋅ m 2
punto 3
3
V3 = 0,077 m
kg
p3 = 70bar
p3V3 70bar * 0,077 *10 5
T3 = = = 1878K
kg 287
pV pV pV p
ó también se puede obtener mediante: = cte ⇒ 2 2 = 3 3 ⇒ T3 = 3 T2
T T2 T3 p2
punto 3A
p3 A = p3 = 70bar
Para determinar T3A es necesario hacer un balance de energía:
q3 − 3 A = C p (T3 A − T3 ) (2)
porque es a presión constante
105
Capítulo 5. Ciclos de trabajo teóricos y reales
qtotal − q2 − 3
(1) = (2): C p (T3 A − T3 ) = qtotal − q2 − 3 ⇒ T3 A = T3 +
Cp
qtotal − Cv (T3 − T2 )
T3 A = T3 + porque es a volumen constante
Cp
Cp
donde γ= y Rg = C p − Cv (3)
Cv
Rg
(3) / Cv : = γ −1 donde
Cv
Rg 0,287 kJ / kg ⋅ K
Cv = ⇒ Cv = = 0,7175 kJ
γ −1 1,4 − 1 kg ⋅ K
− 0,7175(1878 K − 874,5 K )
kJ kJ
2000
kg kg ⋅ K
T3 A = 1878K + = 3158K
kJ
1
kg ⋅ K
punto 4
3
p4V4γ = p3 AV3γA donde V4 = V1 = 1m
kg
γ
⎛V ⎞
1.4
⎛ 0,129 ⎞
p4 = p3 A * ⎜⎜ 3 A ⎟⎟ = 70bar ⎜ ⎟ = 3,98bar
⎝V4 ⎠ ⎝ 1 ⎠
pV 3,98bar * 1m 3 ⋅ kg ⋅ K 10 5 N
T4 = 4 4 = * = 1387 K
Rg kg * 287 N ⋅ m 1bar ⋅ m 2
q1−4 = Cv (T4 − T1 )
q1−4 = 0,7175
kJ
(1387 − 314)K = 769,8 kJ
kg ⋅ K kg
106
Capítulo 5. Ciclos de trabajo teóricos y reales
LITERATURA RECOMENDADA
107
CAPÍTULO 6
INTRODUCCIÓN
Los esfuerzos en las investigaciones actuales en los MCIA van dirigidos hacia mejorar
principalmente el proceso de combustión y el proceso de llenado y vaciado del cilindro o
renovación de la carga. Esto es comprensible debido a la importancia de la masa de aire en la
obtención de la potencia del motor. El propósito de las carreras de admisión y escape consiste en
realizar el barrido de la carga de tal manera que se remuevan los gases al final de la carrera de
potencia y se admita carga fresca para el siguiente ciclo. La calidad del proceso de renovación de
la carga en motores para un motor de cuatro tiempos está representado por el rendimiento
volumétrico el cual depende entre otras cosas del diseño del colector de admisión, puerto de
admisión, válvula de admisión, tipo de combustible, así como de las condiciones de operación del
motor. Una técnica común para mejorar el llenado del cilindro consiste en emplear la
turboalimentación, ésta incrementa el flujo de aire y por lo tanto la densidad de potencia.
En este capítulo estudiaremos básicamente los principios para el cálculo del proceso de
renovación de la carga en motores de cuatro tiempos.
suele abrir entre 40 y 60 grados antes del PMI y suele cerrar entre 15 y 30 grados después del
PMS. La válvula de admisión abre entre 10 y 20 grados antes del PMS. Ambas válvulas
permanecen abiertas durante un período llamado cruce de válvulas. La ventaja de este cruce de
válvulas consiste en mejorar el rendimiento volumétrico especialmente a elevadas velocidades
del motor. La válvula de admisión permanece abierta hasta unos 50 a 70 grados después del
PMI.
En el sistema de admisión de los MEC no hay ni carburador o sistema electrónico de
inyección de combustible ni mariposa. Los MEC suelen llevar turboalimentador con más
frecuencia que los MEP. En los MEC turboalimentados cuando abre la válvula de escape, los
gases quemados fluyen hacia una turbina la cual, mueve un compresor que a su vez se encarga de
comprimir el aire que entra al cilindro.
RENDIMIENTO VOLUMÉTRICO
El rendimiento volumétrico se usa como una medida de la capacidad de llenado del cilindro
durante el proceso de renovación de la carga en un motor de cuatro tiempos. Se define como la
relación entre el flujo másico de aire que entra al cilindro durante el funcionamiento del motor y
el flujo másico de aire que sería desplazado en condiciones teóricas por el motor a las
condiciones de entrada:
ηv =
(masa ⋅ de ⋅ aire )real =
m a ⋅ i
(6.1)
(masa ⋅ de ⋅ aire )teórica ρ a ,i ⋅ Vd ⋅ n
109
Capítulo 6. Proceso de renovación de la carga
p i = p a ,i + p f , i + p H 2 O , i
donde:
pi p f ,i p H 2O ,i
= 1+ +
pa ,i pa ,i pa ,i
y
pa ,i 1
=
pi ⎡ mf m H 2O ⎤
⎢ ⎥
⎢1 + Mf M H 2O ⎥
+
⎢ ma ma ⎥
⎢ Ma M a ⎥⎥
⎣⎢ ⎦
para llegar finalmente a:
pa ,i 1
= (6.2)
pi ⎡ M a m H 2O ⋅ M a ⎤
⎢1 + Fa + ⎥
⎢⎣ M f ma ⋅ M H 2O ⎥⎦
110
Capítulo 6. Proceso de renovación de la carga
pa ,i
En la Figura 6.1 se muestra la variación de en función del dosado relativo.
pi
Figura 6.1 Efecto del tipo de combustible en la presión parcial del aire en la admisión [1]
⎢⎣ a a
( )
⎡m h + 1− x m h + x m h ⎤ = Q +⎛⎜m
e f f , L e f f ,V ⎥⎦ A ⎝
aha + m
⎞
f hf , L ⎟
⎠B
(6.3)
111
Capítulo 6. Proceso de renovación de la carga
112
Capítulo 6. Proceso de renovación de la carga
En la Figura 6.2, para el caso en que pe > pi, es decir a la derecha, mejora el rendimiento
volumétrico al aumentar la relación de compresión, ya que cuando esto sucede, el volumen final
de compresión disminuye, con lo cual hay menos volumen de gases quemados que puedan
retroceder hacia la admisión.
En la misma figura se puede observar, para el caso en que pi > pe, (región a la izquierda),
como mejora el rendimiento volumétrico al disminuir la relación de compresión, ya que cuando
esto sucede, el volumen final de compresión es mayor, con lo cual mejora el barrido de la mezcla.
1. Pérdidas de fricción
Durante la carrera de admisión la presión en el cilindro cae por debajo de la presión atmosférica
en una cantidad que depende del cuadrado de la velocidad del fluido debido a la fricción en cada
parte que compone el sistema de admisión. Esta caída total de presión es la suma de las pérdidas
de presión en cada componente del sistema de admisión: filtro de aire, carburador y mariposa,
colector, puerto, y válvula de admisión. Cada pérdida es un pequeño porcentaje siendo los
principales contribuyentes la válvula y el puerto, como resultado la presión en el cilindro durante
el proceso de admisión cuando el pistón se está moviendo en un punto cercano a su máxima
velocidad puede llegar a ser un 10 a un 20% inferior a la atmosférica. Para cada componente en el
sistema de admisión (y en el escape) la ecuación de Bernoulli da:
Δp = ξ ⋅ ρ ⋅ v 2
j j j
v j A j = S p Ap
113
Capítulo 6. Proceso de renovación de la carga
Donde Aj y Ap son el área de flujo mínima del componente y el área del pistón
respectivamente. Por lo tanto, la caída de presión total cuasi – estacionaria debida a la fricción
es:
2
⎛ Ap ⎞
patm − pc = ∑ Δp j = ∑ ξ j ρv = ρS ∑ ξ j ⎜
2 2 ⎟ (6.4)
j p ⎜A ⎟
⎝ j ⎠
Esta ecuación indica la importancia de un área de flujo grande para reducir las pérdidas de
fricción y la dependencia de éstas con la velocidad del motor.
La Figura 6.3 muestra un ejemplo de la caída de presión debido a la fricción a través del filtro
de aire, carburador, mariposa y colector de un sistema de admisión de un automóvil con motor de
cuatro cilindros standard.
En el sistema de escape se producen pérdidas de presión equivalentes especialmente en el
puerto y en el colector donde los niveles de presión promedios son más altos que la presión
atmosférica. A altas velocidades y altas cargas el colector de escape opera a presiones
substancialmente por encima de la atmosférica.
Figura 6.3 Caída de presión en el sistema de admisión de un MEP de cuatro tiempos determinada
bajo condiciones de flujo estacionario [1]
114
Capítulo 6. Proceso de renovación de la carga
2. Efecto de empuje
La presión en el colector de admisión varía de un proceso de admisión al siguiente para un mismo
cilindro debido a la variación de la velocidad del pistón, a la variación en la apertura de la
válvula, y a los efectos de flujo no estacionario que resultan de esas variaciones geométricas. La
masa de aire inducida dentro del cilindro, y por lo tanto el rendimiento volumétrico, está casi
completamente determinado por el nivel de presión en el puerto de admisión durante el corto
período anterior al cierre de la válvula de admisión. A altas velocidades del motor, la inercia del
gas en el sistema de admisión medida que la válvula de admisión está cerrando incrementa la
presión en el puerto continuando así el proceso de carga a medida que el pistón se desliza hacia el
PMI e inicia su carrera de compresión. Este efecto se vuelve progresivamente más elevado a
medida que incrementa la velocidad del motor. La válvula se cierra entre 40 a 60 grados después
del PMI, en parte para tomar ventaja de este efecto de empuje.
115
Capítulo 6. Proceso de renovación de la carga
cilindro. Si se ajustan adecuadamente los tiempos de esas ondas, las ondas de presión positivas
pueden llegar a incrementar la presión en la válvula de admisión hacia el final de la carrera de
admisión por encima de la presión de admisión nominal. Esto incrementaría la masa de aire
admitida. Cuando sucede este fenómeno se dice que el sistema de admisión está sintonizado.
En la Figura 6.4 se muestran algunos ejemplos de las variaciones de presión en los sistemas
de admisión y escape de un MEP de cuatro cilindros a plena carga. La complejidad del
fenómeno que está ocurriendo es evidente. La amplitud de las fluctuaciones de la presión
aumenta sustancialmente con el incremento de la velocidad del motor. La frecuencia primaria en
la admisión y escape corresponde a la frecuencia de los procesos de admisión y escape de cada
cilindro individual.
Figura 6.4 Presión instantánea en los colectores de admisión y escape de un MEP 4 tiempos y 4
cilindros, a plena carga. p1 colector de admisión a 150 mm desde el cilindro 1; p2 colector de escape
a 200 mm desde el cilindro 1; p3 colector de escape a 700 mm desde el cilindro 1. IO : válvula de
admisión abierta, EO: válvula de escape abierta [1]
Efectos dinámicos
El efecto del flujo sobre el rendimiento volumétrico depende de la velocidad de la mezcla fresca
en el colector, puerto y válvula de admisión. Debido a que las dimensión del sistema y válvula
de admisión se escalan con el diámetro del cilindro, entonces las velocidades en el sistema de
admisión se escalarán con la velocidad del pistón. Por lo tanto los rendimientos volumétricos
116
Capítulo 6. Proceso de renovación de la carga
como una función de la velocidad para diferentes motores deberían ser comparables a la misma
velocidad lineal media del pistón. La Figura 6.5 muestra unas curvas típicas de rendimiento
volumétrico en función de la velocidad media del pistón para un MEC de cuatro cilindros de
inyección indirecta y para un MEP de seis cilindros. Los rendimientos volumétricos de los MEP
generalmente son bajos que los valores de los MEC debido a las pérdidas de flujo en el
carburador y la mariposa, calentamiento del colector de admisión, presencia de vapor de
combustible, y facción de gases residuales más elevada. El doble pico en la curva de los MEC
muestra el efecto de la sintonización del sistema de admisión.
Figura 6.5 Rendimiento volumétrico en función de la velocidad lineal media del pistón [1]
La forma de las curvas del rendimiento volumétrico se puede explicar con la ayuda de la
Figura 6.6. Esta muestra en forma esquemática como varía con la velocidad cada uno de los
fenómenos antes discutidos que afectan el rendimiento volumétrico. Los efectos que no
dependen de la velocidad (tales como la presión del vapor de combustible) son los que bajan la
curva del 100% (curva A). El calentamiento de la carga en el colector de admisión y en el
cilindro hacen que la curva A se baje hasta la curva B. Esta tiene mayor efecto a bajas
velocidades del motor debido a los tiempos prolongados de residencia del gas. Las pérdidas de
fricción del flujo aumentan con el cuadrado de la velocidad del motor y hacen bajar la curva B
hasta la curva C. A elevadas velocidades del motor, el flujo que entra al motor hacia el final de la
carrera de admisión se frena. Una vez que esto ocurre, incrementos adicionales en la velocidad
del motor no aumentan significativamente el flujo de tal manera que el rendimiento volumétrico
disminuye agudamente (de la curva C a la D). El efecto de empuje a altas velocidades hace que
117
Capítulo 6. Proceso de renovación de la carga
Figura 6.6 Efecto sobre el rendimiento volumétrico de diferentes fenómenos que afectan el flujo de aire
en función de la velocidad del motor [1]
TURBOALIMENTACIÓN
La potencia máxima que puede suministrar un motor está limitada por la cantidad de
combustible que se pueda quemar eficientemente al interior del cilindro. Éste está limitado por la
cantidad de aire que se introduce en cada cilindro cada ciclo. Si el aire inducido logra ser
comprimido para que alcance una densidad más alta que la del ambiente, antes de entrar al
cilindro, las potencia máxima de un motor de dimensiones fijas, puede llegar a ser incrementada.
Este es el objeto primordial de la turboalimentación, las ecuaciones (3.38) a (3.40) muestran
cómo la potencia, el par y la presión media efectiva son directamente proporcionales a la
densidad del aire en la admisión [2]. En la turboalimentación, el turboalimentador – un
compresor y una turbina en un solo eje – se usa para incrementar la densidad del aire en la
admisión. La energía disponible en la corriente de gases de escape se emplea para mover una
turbina, la cual a su vez mueve un compresor que incrementa la densidad del flujo que entra en
118
Capítulo 6. Proceso de renovación de la carga
cada cilindro. La Figura 6.7 muestra algunos esquemas típicos de diferentes sistemas de
incremento de la densidad del aire en la admisión. Los sistemas más comunes son la
sobrealimentación mecánica (Figura 6.7a) y la turboalimentación (Figura 6.7b). El esquema de
la Figura 6.7c se usa en grandes motores marinos. La turboalimentación de doble etapa (Figura
6.7d) es una forma de alcanzar presiones muy altas en la admisión (de 4 a 7 atmósferas) con el
fin de incrementar la presión media efectiva. Existen algunos motores, especialmente de uso
militar, llamados hiperbáricos y que consisten de un sistema como el mostrado en la Figura 6.7d
pero que incluye además poscombustión. El esquema mostrado en la Figura 6.7f es de uso
común en los motores actuales. La refrigeración del aire con un intercambiador de calor después
de la compresión y antes de la entrada a los cilindros se usa para incrementar la densidad del aire.
119
Capítulo 6. Proceso de renovación de la carga
Ecuaciones básicas
La ecuación de conservación de la energía para flujo estacionario, aplicada un volumen de
control alrededor de un componente cualquiera se puede expresar de la siguiente manera:
⎡⎛ ⎤
⎢⎜ C 2 ⎞⎟ ⎛
⎜ C 2 ⎞⎟ ⎥
Q − W = m h + − h+ (6.5)
⎢⎜ 2 ⎟ ⎜ 2 ⎟ ⎥
⎣⎝ ⎠ sale ⎝ ⎠ entra ⎦
donde Q es la tasa de transferencia de calor hacia dentro del volumen de control, W es la tasa de
transferencia de trabajo en el eje que sale del volumen de control, m es el flujo másico, h es la
entalpía específica, y C2/2 es la energía cinética específica.
La entalpía total o de parada, ho, se puede definir como:
C2
ho = h + (6.6)
2
Para un gas ideal, con calores específicos constantes, de la ecuación (6.6) se obtiene la
temperatura total o de parada:
C2
To = T + (6.7)
2c p
También se define la presión total o de parada como: la presión que se obtiene si el gas es
dirigido isentrópicamente:
γ (γ −1)
⎛T ⎞
po = p⎜ o ⎟ (6.8)
⎝T ⎠
120
Capítulo 6. Proceso de renovación de la carga
potencia ⋅ reversible
ηC = (6.10)
potencia ⋅ real
La Figura 6.8 muestra los estados finales del gas que pasa a través de un compresor en un
diagrama h-s. Se muestran ambas, las líneas de presión constante estática (p1, p2) y de parada
(p01, p02). El rendimiento isentrópico total – a – total (de condiciones de parada reversibles a
condiciones de parada reales) es, de la ecuación (6.10):
h02 s − h01
ηCTT = (6.11)
h02 − h01
la cual, debido a que cp es prácticamente constante para el aire, o mezcla aire – combustible, se
convierte en:
T02 s − T01
ηCTT = (6.12)
T02 − T01
(γ −1) γ
⎛p ⎞
T02 s = T01 ⎜⎜ 02 ⎟⎟
⎝ p01 ⎠
ηCTT =
( p02 p01 )
(γ −1) γ
−1
(6.13)
(T02 T01 ) − 1
121
Capítulo 6. Proceso de renovación de la carga
T − T01 ( p 2 p 01 )
(γ −1) γ
−1
η CTS = 2s = (6.14)
T02 − T01 (T02 T01 ) − 1
(γ −1) γ
m i c p ,i T01 ⎡⎛ p 02 ⎞ ⎤
− W C = m i c p ,i (T02 − T01 ) = ⎢⎜⎜ ⎟⎟ − 1⎥ (6.15)
η CTT ⎢⎝ p01 ⎠ ⎥⎦
⎣
Donde el subíndice i indica las propiedades de la mezcla de admisión. Si ηCTS se usa para
definir el funcionamiento del compresor, entonces p2 reemplaza a p02 en la ecuación (6.15). La
122
Capítulo 6. Proceso de renovación de la carga
ecuación (6.15) se emplea para calcular la potencia termodinámica requerida para mover el
compresor. También habrán pérdidas mecánicas en el compresor. De este modo la potencia
requerida para mover el aparato será:
W
− WC , D = − C (6.16)
ηm
donde ηm es el rendimiento mecánico del compresor.
La Figura 6.9 muestra los estados del gas a la entrada y salida de una turbina en un diagrama
h-s. El estado 03 es el estado de parada en la entrada; los estados 4 y 04 son respectivamente los
estados estático y de parada de la salida. Los estados 4s y 04s definen los estados estático y de
parada isentrópicos a la salida de la turbina. El rendimiento isentrópico de la turbina se define
como:
potencia ⋅ real
ηT = (6.17)
potencia ⋅ reversible
h03 − h04
ηTTT = (6.18)
h03 − h04 s
Si el gases de escape se modela como un gas ideal con calores específicos constantes,
entonces la ecuación (6.18) se pude escribir así:
Obsérvese que para el gas de escape en el rango de temperaturas de interés, cp puede variar
significativamente con la temperatura.
Dado que la energía cinética en la salida de una turbina usualmente consumida, es más real un
rendimiento isentrópico total – a – estático de la turbina, donde la potencia isentrópica es la que
se obtiene entre las condiciones de parada de la entrada y la presión estática de la salida:
123
Capítulo 6. Proceso de renovación de la carga
La potencia entregada por la turbina está dada por (ecuaciones (6.9) a (6.19)):
⎡ ⎛ p ⎞ (γ e −1) γ e ⎤
WT = m e (h03 − h04 ) = m e c p ,e (T03 − T04 ) = m e c p ,eηTTT T03 ⎢1 − ⎜⎜ 04 ⎟⎟ ⎥ (6.21)
⎢⎣ ⎝ p03 ⎠ ⎥⎦
donde el subíndice e indica las propiedades del gas de escape. Si el rendimiento total – a –
estático de la turbina (ηTTS ) se usa en la relación para WT , entonces p4 reemplaza a p04 en la
ecuación (6.21).
Efecto de la altitud
A medida que incrementa la altura sobre el nivel del mar, la densidad del aire disminuye. A
manera de ejemplo, aplicando la ley de los gases ideales al aire en diferentes ciudades de
Colombia tendríamos:
p amb
ρ aire =
RaireTamb
124
Capítulo 6. Proceso de renovación de la carga
Ciudad Altura s.n.m. pamb (mbar) Tamb (C) Densidad (kg / m3) Potencia (kW) 1
Barranquilla 0 1013 30 1,1648 69,7 (93,3 hp)
Medellín 1500 853 22 1,0075 60,3 (80,7 hp)
Bogotá 2600 752 18 0,90042 54 (72 hp)
1
Calculada con la ecuación (3.38), suponiendo como único dato variable en la ecuación el valor de la densidad del
aire (ηV = 0,8; ηf = 0,3; Vd = 2000cc; nm = 5000rpm; Fabs = 1 / 14,7; Hc = 44000 kJ / kg).
125
Capítulo 6. Proceso de renovación de la carga
Pérdidas mecánicas
En un turboalimentador o turbocompresor, la turbina está mecánicamente unida al compresor,
por lo tanto a velocidad constante, se cumple que,
Siendo ηm el rendimiento mecánico del grupo turbocompresor. Las pérdidas mecánicas son
básicamente producidas por los rodamientos. El rendimiento mecánico generalmente se combina
con el de turbina debido a que estas pérdidas son difíciles de determinar. Sin embargo,
Winterbone et al, [3] han propuesto calcular el par de pérdidas de fricción por rodamiento
mediante la siguiente correlación:
⎡ ⎛ n (rpm ) − 30000 ⎞ ⎤
M fric ( N ⋅ m ) = ⎢0,85⎜ TC
60
⎟ + 0,25⎥ ⋅ 1000,325 ⋅ (6.23)
⎣ ⎝ 30000 ⎠ ⎦ 2π ⋅ nTC (rpm )
El Compresor
La mayoría de MCIA empleados en automoción que llevan turbocompresor, suelen llevar
compresor centrífugo de una sola etapa. Se caracterizan por manejar elevados flujos másicos a
relaciones de compresión relativamente bajas (máximo de 3.5). Su mejor funcionamiento se
consigue a elevadas velocidades angulares por lo que se suelen acoplar en un solo eje con la
turbina. Consisten de una carcaza estacionaria, un rotor con álabes, un difusor estacionario (con
o sin aspas), y una voluta (colector) que lleva al motor el aire que sale del difusor (Figura 6.10)
126
Capítulo 6. Proceso de renovación de la carga
En la Figura 6.11 se indica en un diagrama h-s la manera como cada componente contribuye
al incremento de presión a través del compresor.
127
Capítulo 6. Proceso de renovación de la carga
128
Capítulo 6. Proceso de renovación de la carga
La Figura 6.13 muestra las curvas características de un compresor centrífugo real. En la práctica,
las variables corregidas (llamadas pseudoadimensionales), velocidad y flujo másico se definen de
la siguiente manera:
Tref
ncor = nTC
To, in
(6.24)
To, in ⎛ p ref ⎞
m cor = m ⎜ ⎟
Tref ⎜⎝ po, in ⎟⎠
129
Capítulo 6. Proceso de renovación de la carga
La Turbina
La turbina es movida por la energía disponible en los gases de escape. En la Figura 6.14 se
observa la energía disponible ideal. Ésta está compuesta por el trabajo producido por el soplado
de los gases de escape al momento de la apertura de la válvula hasta la presión atmosférica (área
abc) y (para el motor de cuatro tiempos) por el trabajo ascendente realizado por el pistón en el
desplazamiento de los gases remanentes en el cilindro después del soplado (área cdef).
Figura 6.14 Diagrama p-V ciclo a volumen constante mostrando la energía disponible en el escape. [1]
Los MCIA son por naturaleza aparatos de flujo pulsatorio no estacionario. Las turbinas
pueden ser diseñadas para aceptar este tipo de flujo, pero éstas operan más eficientemente bajo
condiciones de flujo estacionario. Los tipos de turbinas empleadas en automoción generalmente
son de flujo radial. Éstas son aparentemente iguales al compresor centrífugo, la diferencia radica
en que en la turbina el flujo entra radialmente mientras que en el compresor el flujo sale
radialmente. En la Figura 6.15 se muestra una turbina de flujo radial de una sola entrada. Ésta
consiste de una carcaza o caracol, un conjunto de pequeñas toberas (a menudo se omiten en
turbinas pequeñas) y el rotor de turbina.
130
Capítulo 6. Proceso de renovación de la carga
Figura 6.16 (a) diagrama h-s para una turbina real. [1]
131
Capítulo 6. Proceso de renovación de la carga
132
Capítulo 6. Proceso de renovación de la carga
Ejercicio 6.1
De un motor diesel de inyección directa turboalimentado, 6 cilindros en línea, se han obtenido en
banco de ensayos, en régimen estacionario, los datos que aparecen en la tabla adjunta. El
combustible empleado es un ACPM con una densidad de 0.82 kg/dm3. Las condiciones
atmosféricas del lugar son pamb = 933 mbar y Tcelda = 25C. Calcular la potencia del compresor,
la potencia de turbina y el rendimiento volumétrico.
Solución
Para el compresor tenemos:
1
2
m 1a = m 2 a
133
Capítulo 6. Proceso de renovación de la carga
ρ1V1 = ρ 2V2
m a = ρ a1 ' ⋅Va1
T = 25C (298K )
W C = m a1 .(h2 − h1 ) → 1
T2 = 127C (400 K )
Con estas temperaturas obtengo, de tabla de propiedades del aire [1], los valores de las entalpías:
kg kJ
WC = 0,218 .(560 − 460) = 21,8kW ♦
s kg
Para la turbina tenemos:
4
3
T3 = 600 C T4 = 535 C
p3 = 859 mbar p4 = 36,7 mbar
m e = m a + m f
kg dm 3 kg
m f = 0,82 .0,0107204 = 8,7907.10 − 3
dm 3 s s
kg
m e = 0,218 + 0,0087907 = 0,2268
s
134
Capítulo 6. Proceso de renovación de la carga
Método 1
kJ
Suponiendo ηTs = 0,7 , C pe = 1,3 y γ e = 1,33
kg º K
⎧ ⎟⎫
⎛ 0 , 33 ⎞
⎜
⎪ ⎡ (51,48 + 933) ⎤ ⎜⎝ 1,33 ⎟⎠ ⎪
WT = 0,2268 × 0,7 × 1,3 × 873⎨1 − ⎢ ⎥ ⎬ = 24,8kW ♦
⎪ ⎣ (859,6 + 933) ⎦ ⎪
⎩ ⎭
Este resultado implica que el grupo turbocompresor tiene un rendimiento mecánico de:
W C 22
ηm = = = 0,88
WT 24'8
Método 2
Partiendo de la composición del combustible y suponiendo combustión estequiométrica:
[ ]
1 C10,8 H18,7 + 15,475(O2 + 3,773 N 2 ) → 10,8CO2 + 9,35H 2 O + 58,186 N 2
kgComb
M f = 10,8(12) + 18,7(1) = 148,3 .1kMol = 148,3kgComb
kMol
kgAire
M a = (15,475 × 4,773)kMol × 28,962 = 2139,2kgAire
kMol
148,3 A
Fstq = = 0,06932 ⇒ = 14,42
2139,2 F
Pero el motor está funcionando a unas condiciones diferentes de las estequiométricas, por lo
tanto es necesario conocer el dosado de funcionamiento o dosado absoluto:
m f 8,7907.10 −3
Fabs = = = 0,0403
m a 0,218
135
Capítulo 6. Proceso de renovación de la carga
0,0403
FR = = 0,5817
0,06932
1
FR −1 = λ = = 1,72 ⇒ 72%
0,5817
Con esto, puedo determinar una ecuación de reacción más próxima a la realidad:
[ ]
1 C10,8 H18,7 + 1,72 × 15,475(O2 + 3,773N 2 ) → 10,8CO2 + 9,35H 2 O + 100,4 N 2 + 11,14O2
Para el CO2:
~ kJ
C p co = 8,3143 .{2,4008 + 0,87351 × 10 − 2 (873) −
2 kMol º K
0,66071 × 10 − 5 (873) 2 + 0,20022 × 10 − 8 (873) 3 + 0,63274 × 10 −15 (873) 4 }
~ kJ
C PCO = 52,576
2 kMol º K
H2O:
~
C PH = 8,3143 × {4,0701 − 0,11084 × 10 −2 (873) + 0,41521 × 10 −5 (873) 2 −
2O
136
Capítulo 6. Proceso de renovación de la carga
N2:
~ ~ kJ
C PN = 8,3143 × {....} ⇒ C PN = 31,92
2 2 kMol º K
O2:
~ kJ
C PO = 8,3143 × {.....} ⇒ C PO = 34,21
2 2 kMol º K
∑ X i ⋅ C Pi
~ ~ ~ ~ n
C PPr od = Donde X i = Fracción molar del componente i = i
i
n
~
~ kJ C PPr od
C PPr od = 43,43 Donde C PPr od =
kMol º K M Pr od
M Pr od =
1
[10,8 × 44 + 9,35 × 18 + 100,4 × 28 + 11,14 × 32]
131,69
kg
M Pr od = 28,94
kMol
kJ kMol kJ
C PPr od = 43,43 × = 1,5
kMol º K 28,94 Kg kg º K
Para determinar γ:
8,3134
C p − C v = R g ⇒ Cv = C p − R g Donde Rg =
28,94
137
Capítulo 6. Proceso de renovación de la carga
Cp 1,5
γ = = = 1,237
Cv 1,2127
873 − 808
ηT = = 0,69
⎛ 0 , 237
⎞
⎜ ⎛ 984,5 ⎞ 1, 237 ⎟
873⎜1 − ⎜ ⎟ ⎟
⎜ ⎝ 1792,6 ⎠ ⎟
⎝ ⎠
⎡ 0 , 237
⎤
⎢ ⎛ 984,5 ⎞ 1, 237
⎥ = 22,3kW
WT = 0,2268 × 1,5 × 0,69 × 873 1 − ⎜ ⎟ ♦
⎢ ⎝ 1792,6 ⎠ ⎥
⎢⎣ ⎥⎦
i.m a p2 2,0126 × 10 5 Pa kg
ηv = , Donde ρ 2 = = = 1,753
ρ a, i .η m .Z .Vd R g T2
287,3
J
× 400 K m3
kg º K
kg
2 × 0,218
ηv = s = 0,97
kg rev 1min −3 3
1,753 × 2000 × × 7,685 × 10 m
m3 min 60seg
Este valor del rendimiento volumétrico mide la capacidad de llenado de todo el sistema de
admisión. Obsérvese que se ha tomado el valor de densidad después del compresor y no el valor
de la densidad del sitio, para efectos de cálculo. Para este último caso, el rendimiento daría
mayor que 1, con lo cual hablaremos de un coeficiente de llenado, propio de motores fuertemente
turboalimentados.
Ejercicio 6.2
De un motor diesel de automoción de aspiración natural funcionando en Medellín (850 mbar
y 20 ºC), se conocen los siguientes datos:
138
Capítulo 6. Proceso de renovación de la carga
Solución
Sin turbo:
η v ⋅ ρ a , i ⋅ Vd ⋅ n 0,85 × 1 × 2 × 10 −3 × 4200
m a = = = 0,0595kg / s
i 2 × 60
m f
Fa = FR ⋅ FS = ⇒ m f = m a ⋅ FR ⋅ FS
m a
2,77389 × 10 −3 kg / s
gef = = 0,2 kg = 200 g / kWh
50kW *1h / 3600 s kWh
Con turbo:
η v = 0,9
gef = 0,2 × 0,9 = 0,180kg / kWh
FR = 0,68
p2 1,6 × 10 5 kg
ρ2 = = = 1,553
R g T2 287 × 359 m3
139
Capítulo 6. Proceso de renovación de la carga
m f
Fa = FR ⋅ FS = ⇒ m f = m a ⋅ FR ⋅ FS
m a
4,4309 × 10 −3 kg / s
Pe = = 88kW ♦
kg 1h
0,180 ×
kW − h 3600s
LITERATURA RECOMENDADA
REFERENCIAS
[1] Heywood, J.B., (1988), “Internal Combustion Engine Fundamentals”, McGraw-Hill, New
York.
[2] Watson, N. And Janota, M.S., (1982), “Turbocharging the Internal Combustion Engine”,
Macmillan Press, London.
[3] Winterbone, D.E., Horlock, J.H., (1986), “The Thermodynamics and Gas Dynamics of
Internal Combustion Engines”, Clarendon Press, Oxford.
140
CAPÍTULO 7
INTRODUCCIÓN
El calor cedido a las paredes del cilindro y arrastrado por el agua de refrigeración puede ser
determinado realizando algunas hipótesis simplificadoras. Éste incluye [1]:
No incluye el calor arrastrado por el aceite lubricante, ni las pérdidas por radiación y
convección de las paredes externas del cilindro, pipas, colectores, etc. Tampoco incluye una
pequeña proporción de calor que absorbe el aire en el momento de entrar al cilindro.
Es casi imposible distinguir la proporción de calor transportada por el aceite lubricante (que
de otra manera hubiera pasado a las camisas del cilindro), de aquel que ha sido generado por la
fricción en los cojinetes, etc.
Sin embargo es posible evaluar con una aproximación regular la pérdida de calor por
radiación y convección de las paredes externas y esto puede alcanzar una proporción bastante
grande si las pruebas son llevadas a cabo bajo una fuerte corriente de aire, o con temperaturas del
medio refrigerante bastante altas.
eficiencia que corresponde a la eficiencia de la expansión sola (es decir, excluyendo el trabajo
negativo hecho durante la compresión) lo cual en un motor con una relación de compresión de
5:1 tendría una recuperación de cerca de 40%. El 60% restante del calor, así recuperado, sería en
cualquier caso, rechazado al escape después de la expansión.
Calor generado por la fricción del pistón. Este tema, aunque es importante, es difícil de evaluar
debido a la dificultad de separar el calor generado por la fricción de aquel que ha entrado al
142
Capítulo 7. Pérdidas de calor y refrigeración
pistón proveniente de los gases de combustión y ha sido transferido a las paredes del cilindro.
También variaría, en una gran extensión, según el diseño del pistón, el número de anillos, la
viscosidad del lubricante y otros factores. No todo el calor que entra al pistón proveniente de
cualquier fuente, encuentra necesariamente camino hacia las paredes del cilindro, puesto que una
proporción sustancial será llevada fuera por circulación tanto de aceite como del aire que está
dentro del cárter del motor.
Pruebas separadas realizadas por medio de motores de ensayo, bajo condiciones lo más
aproximadas posible a las condiciones normales de rodamiento, muestran que el calor generado
por la fricción del pistón usualmente está en rangos entre el 1% y el 1.5% del poder calorífico
inferior total del combustible; la mayor parte de este calor por fricción encontrará su camino a
través de las paredes del cilindro.
En los MEC, la distribución del calor varía mucho respecto al tamaño del cilindro y a la
forma de la cámara de combustión, y es por lo tanto imposible hacer generalizaciones razonables
como en el caso de los MEP.
El rechazo de calor al agua de enfriamiento a plena carga y a velocidad máxima, para los
motores de cámara de turbulencia de tamaño pequeño, es alto, llegando hasta 1.25 veces el calor
equivalente de la potencia efectiva. Esto se compara con cerca de 0.8 veces que es de un motor
de gasolina típico. Este calor tan alto crea dificultades cuando se quiere cambiar el motor de
gasolina de un vehículo por un diesel. Si se necesita un poder similar máximo, el radiador puede
ser inadecuado pero afortunadamente en muchos casos se adopta una potencia menor y el mismo
radiador sirve.
Sin embargo, con cargas parciales la posición se invierte (excepto a la velocidad máxima).
Esto se debe al efecto de un ciclo de temperatura más bajo que predomina sobre una densidad de
carga alta y un intenso movimiento de aire.
Como se mencionó anteriormente al discutir el MEP, notamos que la pérdida de calor durante
el período de combustión, era calor que de otra manera habría podido convertirse en trabajo
indicado, a una eficiencia de cerca del 40% para una relación de compresión de 5:1 y
correspondientes mayores para relaciones más altas. En los MEC, especialmente en los del tipo
de turbulencia por compresión, las pérdidas de calor tienden a ser mayores debido a un aumento
de densidad y un movimiento de aire intenso. Además la eficiencia de conversión de calor en
trabajo sería más alta, de aquí que la importancia de la pérdidas de calor durante la fase de
combustión sea mayor. Argumentos similares se aplican a la fase de expansión, aunque la
eficiencia de la utilización del calor ahorrado sería menor.
143
Capítulo 7. Pérdidas de calor y refrigeración
144
Capítulo 7. Pérdidas de calor y refrigeración
alcancen valores de hasta 10MW / m2 durante el período de combustión. Sin embargo, en otras
partes del ciclo, el flujo de calor llega a ser prácticamente cero. El flujo varía sustancialmente
con la ubicación; aquellas regiones de la cámara que están en contacto con el movimiento rápido
de los gases quemados a elevada temperatura son las que están sometidas a los flujos más altos.
En estas regiones, se deben mantener las temperaturas lo suficientemente bajas para lograr evitar
la fatiga térmica (inferiores a 400C para fundiciones grises y de 300C para aleaciones de
aluminio). Las paredes del cilindro del lado del gas, deberán mantenerse a una temperatura
inferior a 180C para evitar el deterioro de la película de aceite de lubricación. La bujía y las
válvulas, especialmente la de escape, tienen que refrigerarse para evitar problemas de detonación
del gas final (end gas) y de pre – encendido.
La transferencia de calor afecta el funcionamiento, el rendimiento y las emisiones del motor.
Para una masa de combustible dada dentro del cilindro, se cumple que a mayor transferencia de
calor hacia las paredes de la cámara, serán más bajas las presiones y temperaturas medias de los
gases de combustión, lo que reduce el trabajo por ciclo transmitido al pistón.
Figura 7.1 Esquema del proceso de transferencia de calor del cilindro al refrigerante en un MCIA. qCV,
qR y qCN son las tasas de transferencia de calor por convección, radiación y conducción, respectivamente.
T es temperatura y los subíndices g, p y r indican gas, pared y refrigerante, respectivamente. [2]
145
Capítulo 7. Pérdidas de calor y refrigeración
q = qCond =
k
(Tpared ,g − Tpared ,r ) (7.2)
e
q = q CV = hr (T pared ,r − Tr ) (7.3)
donde hr es el coeficiente de película del lado del refrigerante que lo determinaremos en este libro
con el método propuesto por Taylor y Toong (ecuación (7.10)). [2]
Reuniendo términos para hallar el flujo de calor y la diferencia de temperaturas se llega a:
q =
1
Tg − Tr( ) (7.4)
1 e ⎛ 1 Ag ⎞
+ +⎜ ⋅ ⎟
hg k w ⎜⎝ hr Ar ⎟⎠
(T pg − Tpr = ) kw
1
⎛ k w Ag ⎞
(
Tg − Tr ) (7.5)
+ 1 + ⎜⎜ ⋅ ⎟⎟
e ⋅ hg ⎝ e ⋅ hr Ar ⎠
⎛ k Ag ⎞
1 + ⎜⎜ w ⋅ ⎟⎟
⋅
Tpg = Tr + ⎝ e h r A r ⎠
(Tg − Tr ) (7.6)
kw ⎛ k Ag ⎞
+ 1 + ⎜⎜ w ⋅ ⎟⎟
e ⋅ hg ⎝ e ⋅ hr Ar ⎠
146
Capítulo 7. Pérdidas de calor y refrigeración
En la Tabla 7.1 se puede observar el efecto del incremento del gasto másico, del tamaño y de
la velocidad del motor sobre los tres parámetros de las ecuaciones (7.4) a (7.6).
El coeficiente de película (hg) se obtiene a partir de la relación entre los números de Nusselt,
Reynolds y Prandtl: Nu = a Re b Pr c . El método propuesto por G. Woschni [3] se empleará en
este libro para el cálculo del hg.
En términos de los parámetros característicos ésta expresión se puede expresar como:
c
⎛ Cpμ ⎞
b
hg L ⎛ ρvL ⎞ ⎜ ⎟
= a⎜⎜ ⎟⎟
k pared ⎝ μ ⎠ ⎜k ⎟ (7.7)
⎝ pared ⎠
ó lo que es lo mismo:
Tabla 7.1
Efecto del incremento del gasto másico, tamaño y velocidad del motor en la tasa de transferencia
de calor
Variable q Tpg - Tpr Tpg Debido a:
Gasto másico ↑ ↑ ↑ Re, hg, F
Tamaño motor ↓ ↑ ↑ e, hg
Velocidad motor ↑ ↑ ↑ Re
Debido a que la solución simultánea de las ecuaciones (7.4) a (7.7) supone una gran
dificultad, especialmente debido a la vairación espacial y temporal del coeficiente de película, se
suelen emplear técnicas de cálculo simplificadoras que emplean correlaciones semiempíricas.
147
Capítulo 7. Pérdidas de calor y refrigeración
Según Watson y Janota [4] el que mayor influencia tiene de los tres sobre el calor transferido
total es el primero de ellos (del gas a las paredes). Por esa razón la mayoría de los modelos de
transferencia de calor en los cilindros de los MCIA se concentran en la zona gas - paredes. Así
pues la transferencia de calor dependerá del gradiente de temperaturas en la capa límite en las
diferentes superficies que rodean el sistema (el gas). Algunos autores incluyen en sus estudios la
transmisión de calor por radiación, atribuyéndole distintos grados de importancia dependiendo
del tipo de motor. La contribución de la radiación en el calor total transmitido es difícil de
determinar debido a la dificultad de medir la temperatura instantánea de la llama y el flujo de
calor. Las diferentes estimaciones encontradas en la literatura le atribuyen desde un 0 hasta un
45% del calor total transmitido [5 a 7].
La expresión que se emplea para calcular el calor transmitido es la siguiente:
[ ] (
Q = hg ⋅ Ap ⋅ (Tg − Tp ) + Acul ⋅ (Tg − Tcul ) + Acamisa ⋅ (Tg − Tcamisa ) + εσ Tg4 − Tpared
4
,g ) (7.9)
[
Q = hg ⋅ Ap ⋅ (Tg − Tp ) + Acul ⋅ (Tg − Tcul ) + Acamisa ⋅ (Tg − Tcamisa )] (7.10)
Para el cálculo del coeficiente medio de transferencia de calor (hg), Woschni utiliza la relación
entre los números de Nusselt (htl/k)y de Reynolds ( ρCl / μ ) de la siguiente manera:
148
Capítulo 7. Pérdidas de calor y refrigeración
b
ht l ⎛ ρCl ⎞
= a ⋅ ⎜⎜ ⎟⎟ (7.11)
k ⎝ μ ⎠
donde:
l = longitud característica (diámetro del pistón)
k = conductividad térmica del gas
ρ = densidad del gas
μ = viscosidad dinámica del gas
C = velocidad característica del gas (proporcional a la velocidad lineal media del
pistón)
a y b son coeficientes que se ajustan experimentalmente según el tipo de motor y
sus condiciones de funcionamiento.
C = k4·Cm
Durante la combustión, Woschni establece que el movimiento del gas en el cilindro debería
ser una función del incremento de presión por encima de la presión a motor arrastrado. Así pues
durante las fases de combustión y expansión el término de la velocidad característica queda:
⎛ V ⋅T ⎞
C = 2.28 ⋅ C m + 0.00324 ⋅ ⎜⎜ D RCA ⎟⎟ ⋅ ( p − p arrastrado ) (7.12)
⎝ p RCA ⋅ V RCA ⎠
donde:
Cm = velocidad lineal media del pistón (m/s)
VD = cilindrada del motor (m3)
149
Capítulo 7. Pérdidas de calor y refrigeración
TRCA, PRCA y VRCA = temperatura (K), presión (kPa) y volumen (m3) en el momento en que
cierra la válvula de admisión y
(p – parrastrado) = presión en el cilindro menos presión a motor arrastrado
b
⎛ W ⎞ ⎛m⎞
h g ⎜ 2 ⎟ = 3,26 ⋅ D p (m ) ⋅ p (kPa ) ⋅ T (K )
b −1 (0.75−1.62b )
⋅ C⎜ ⎟
b
(7.13)
⎝m K ⎠ ⎝s⎠
donde:
b = 0.8
Dp = diámetro del pistón (m)
P = presión en el cilindro (kPa)
T = temperatura en el cilindro (K)
150
Capítulo 7. Pérdidas de calor y refrigeración
kg kg
hm = 10.4 ⋅ Re 0.75 = 10.4 (ρ ⋅ Cm )0.75 D p−0.25 (7.16)
Dp μ 0.75
g
Los valores de las propiedades de transporte para la ecuación de Taylor y Toong se obtienen
de la Figura 7.2.
Figura 7.2 Correlación global de transferencia de calor en el motor: número de Nusselt del lado del gas
- número de Reynolds para diferentes tipos de motores. Tg,a temperatura media del gas, μg viscosidad
absoluta del gas y kg conductividad térmica del gas. Datos obtenidos para una Tamb de 27C [2].
151
Capítulo 7. Pérdidas de calor y refrigeración
Tabla 7.2
Pérdidas en el motor en función del poder calorífico inferior del combustible
Porcentaje respecto al poder calorífico inferior del combustible
Motor Ne Qref Qaceite He,inc Mhe
MEP 25 - 28 17 - 26 3 - 10 2-5 34 - 45
MEC 34 - 38 16 - 35 2-6 1-2 22 - 35
152
Capítulo 7. Pérdidas de calor y refrigeración
Figura 7.3 mfHc = potencia suministrada; Pi = potencia indicada; Pe = potencia efectiva; Ptf = potencia
de fricción total; Ppf = potencia de fricción en el pistón; He = entalpía en el escape; He,s,a = flujo de
entalpía que se va a la atmósfera; He,ic = flujo entálpico asociado a la combustión incompleta; Ec,k = flujo
de energía cinética en el escape; Qpared= calor a las paredes; Qr,e = calor cedido al refrigerante en los
puertos de escape; Qe,r = flujo de calor radiante desde el escape
153
Capítulo 7. Pérdidas de calor y refrigeración
Consideremos el caso de un MEC de aspiración natural, usado para automoción, con una
relación de compresión de 15. El rendimiento indicado es 45%, y 25% de la energía del
combustible es arrastrada por el agua de refrigeración. De ese 25%, cerca del 2% es debido a la
fricción. De el 23% restante, cerca del 8% son pérdidas de calor durante la combustión, 6% son
pérdidas de calor durante la expansión, y 9% pérdidas de calor durante el escape. Del 8% de las
pérdidas durante la combustión, casi un 4% podría ser convertida en trabajo útil en el pistón. Del
6% de las pérdidas de calor durante la expansión, cerca del 2% podría haber sido utilizada. Así,
del 25% perdido en el refrigerante, sólo cerca de un 6% podría recuperarse en forma de trabajo
útil, lo que permitiría incrementar el rendimiento indicado de 45% a 51%.
Para un MEP, la conversión a trabajo útil es aún más baja, debido a las bajas relaciones de
compresión.
Ejercicio 7.1
La tasa instantánea de transferencia de calor ( Q ) de los gases del cilindro a las paredes de la
cámara de combustión en un MEP puede ser estimada aproximadamente de la ecuación:
Q = hg A(Tg − T p )
Nu = c(Re) m (Pr) n
La velocidad y longitud características usadas en esta relación son la velocidad media del
pistón y el diámetro del cilindro respectivamente. Asumiendo valores apropiados para la
geometría del motor y condiciones de funcionamiento de mariposa totalmente abierta, con la
temperatura media de pared a 400 K, a una velocidad del motor de 2500 rpm, y usando la presión
del cilindro vs. ángulo del cigüeñal que se anexa, calcular:
154
Capítulo 7. Pérdidas de calor y refrigeración
1. La temperatura media del gas en el cilindro para: θ = -180º, -90º, 0º, 20º, 40º, 90º, 150º.
2. El coeficiente instantáneo de transferencia de calor hg y la tasa de transferencia de calor Q
del gas a las paredes de la cámara de combustión de un cilindro en estos ángulos de cigüeñal.
Grafique estos resultados en función de θ.
3. Estime la fracción de energía del combustible que es transferida a las paredes del cilindro
durante la compresión y la expansión.
Datos:
Asuma para el gas que la viscosidad μ = 7 ×10 −5 kg/m s , la conductividad térmica
k = 1.5 × 10 −1 J / m ⋅ s ⋅ K , el peso molecular = 28 y el número Prandtl es 0.8. Asuma que la
cámara de combustión tiene forma de disco con Dp = 102 mm, L = 88 mm y rc = 9.
Solución
Se leen de la figura los datos de P (θ) y T (θ) para los ángulos que se Piden. Para calcular el
número de Reynolds, se calcula la densidad ρ (θ) como si fuera un gas ideal. Se consideran
constantes dentro del ciclo las propiedades de transporte suministradas por el problema. Con el
Re y el Nu se determina el coeficiente de película (hg) y finalmente se determina el flujo de calor
155
Capítulo 7. Pérdidas de calor y refrigeración
(Q ) convectivo a las paredes del cilindro. Se debe tener en cuenta que el área de transferencia de
calor es la suma de: pistón, culata y camisa. La relación biela - manivela para un motor Otto de
estas características suele estar entre 3 y 4. Haremos el ejercicio para un solo ángulo y a
continuación se generalizará la solución, presentando los resultados que se piden en una tabla.
Usando la curva de presión vs. ángulo de giro para un valor de –180º tenemos que:
Datos:
ρ ⋅ Cm ⋅ φ p
Nu = 0.4 Re 0.75 Pr 0.4 Re =
μg
rev 1m 1min
Cm = 2 Ln = 2 ⋅ 88mm ⋅ 2500 ⋅ 3 ⋅ = 7,3m / s
min 10 mm 60 seg
pM
pv = nRu T ; Tomando gas dentro del cilindro como gas ideal, tenemos que: ρ=
Ru T
kg
de donde ρ = 0.842
m3
kg m
0842 3
⋅ 7,3 ⋅ 0.102m
Re = m s = 8996,2
−5 kg
7 ⋅ 10
ms
W
337,95 ⋅ 1,5 ⋅ 10 −1
Nu ⋅ kg mK ≈ 497 W
hg = =
φp 0,102m m2 K
AT = AP + AC +
πDS
2
(R + 1 − cosθ − )
R 2 − sen 2 θ = 2,4 A p +
πDS
2
(4,5 − cosθ − 3,5 2 − sen 2 θ )
156
Capítulo 7. Pérdidas de calor y refrigeración
π
AP = (0.102) 2 = 8.171× 10 −3 m 2
4
AT = 2.4(8.171×10 −3 ) +
π (0.102)(0.088)
2
(3.5 +1− cos(−180) − )
3S 2 − sen 2 (−180) = 47,81⋅10 −3 m 2
Repitiendo el procedimiento para los otros valores del ángulo de giro y asumiendo que μ, K y
Pr no cambian con la temperatura ni con la presión, se tiene:
θ T p ρ hg AT
Re Nu Q (W)
grados K kPa kg/m3 W/m2ºK m *10-3
2
157
Capítulo 7. Pérdidas de calor y refrigeración
Q hg
100000 2000
2
1600
60000
1200
40000
800
20000
400
0
-20000 0
-200 -100 0 100 200
Ángulo de giro de cigüeñal [º]
Tasa de calor convectivo y coeficiente medio de transferencia de calor en función del ángulo de
giro de cigüeñal, para los resultados del ejercicio 7.1
Ejercicio 7.2
Se tiene un motor diesel turboalimentado con intercooler, su cilindrada son Vd = 9 lt, Este motor
es capaz de proporcionar 235 kW a 2200 rpm, Dp = 116 mm, L = 142 mm. El motor funciona con
un dosado relativo de 0,7, la presión a la salida del compresor es de 1,7 bar, el rendimiento
isentrópico del compresor es 0,75, el rendimiento del intercooler es 0,85, la temperatura media
del refrigerante es de 95 ºC, si Tamb = 20 ºC y Pamb = 1 bar, calcule el flujo de calor convectivo
que tiene que arrastrar el refrigerante.
El motor tiene 6 cilindros en línea y el intercooler es del tipo aire-aire.
Solución
Para un FR = 0’7 se hace de la figura de Taylor y Toong:
T = 600 K
μ g = 3,0 × 10 −5 Pa ⋅ S
W
K g = 4,5 × 10 −2
mK
(
Q C = z ⋅ Ac ⋅ hg Tg' − Tref ) (1)
Donde: Z = Número de cilindros
158
Capítulo 7. Pérdidas de calor y refrigeración
Ac = Area convectiva
hg = Coeficiente de película convectivo
Tg ' = Temperatura media del gas corregida
hg se calcula con la ecuación de Taylor y Toong
kg
hg = 10'4 * * Re 0.75 (2)
Dp
Por lo tanto para poder determinar hg se necesita antes conocer el número de Reynolds:
ρ 2' ⋅ C m ⋅ D p
Re = (3)
μg
P2 '
ρ 2' = (4)
Rg T2 '
293 K 1
1 bar 4
ηC = 0,75 C T
2 1,7 bar
3
ηi = 0,85
293 K Ta
1 bar
2’
Motor
159
Capítulo 7. Pérdidas de calor y refrigeración
En la ecuación (5) es necesario determinar T2 , para ello se utiliza la ecuación del rendimiento
isentrópico del compresor:
h2 S − h1 T2 S − T1 T2 S − T1
ηc = = despejando: T2 = T1 + (6)
h2 − h1 T2 − T1 ηc
γ −1 γ −1
T2 S ⎛ P2 ⎞ γ ⎛P ⎞ γ
=⎜ ⎟ así pues: T2 S = T1 ⋅ ⎜⎜ 2 ⎟⎟ donde:
T1 ⎜⎝ P1 ⎟⎠ ⎝ P1 ⎠
341 − 293
T2 = 293 + = 357 K (8)
0'75
Llevando este valor (8) a la ecuación (5)
1,7 × 10 5 Pa kg
ρ 2' = = 1,955 3 (10)
J m
287 ⋅ 303K
kgK
Para conocer entonces el número de Reynolds se necesita C m (velocidad lineal media del
pistón):
C m = 2 Ln = 2(0,142) ⋅ 2200 / 60 = 10,41m / s
160
Capítulo 7. Pérdidas de calor y refrigeración
kg m
1,955 3
⋅ 10,41 ⋅ 0,116m
Re = m s = 78693
−5 N 1kg 9,81m
3 × 10 ⋅s⋅ ⋅ 2
m2 9,81N s
4,5 × 10 −2
⋅ (78693) = 18960 2
W
hg = 10,4 ⋅
0 , 75
0,116 m K
⎡ kg ⎤ Q c
m H 2O ⎢ ⎥ =
⎣ s ⎦ C PH 2O ⋅ ΔT
Donde ΔT es el salto de temperaturas en el intercambiador de calor (radiador) y el Cp del agua se
puede asumir con un valor de 4,18 kJ / kg-K
Nota
El problema se podía solucionar igualmente determinando el número de Nusselt a partir del
Reynolds, leyendo este dato de la gráfica de Taylor y Toong.
Ejercicio 7.3
Calcular el flujo de calor convectivo que tiene que evacuar el sistema de refrigeración de un
motor diesel de 6 cilindros en línea turboalimentado con intercooler. Se conocen los siguientes
161
Capítulo 7. Pérdidas de calor y refrigeración
datos: diámetro pistón 112 mm, carrera 130 mm. Rendimiento isentrópico del compresor 0.7,
rendimiento intercooler aire – aire 0.8, presión a la salida del compresor 1,5 bares, dosado de
funcionamiento del motor 0.6, temperatura media del refrigerante 95 ºC. Condiciones ambiente
20 ºC y 850 mbar. El motor desarrolla 150 kW a 2000 rpm.
Solución
T1 =293 ºK
1 4
2 3
2’
W N ⋅s Tg = 550º K
kg = 4,4 × 10 − 2 ; μ g = 3,5 × 10 −5 ;
m⋅K m2 Tr = 368º K
Qr = hm (Tg − Tr )Ap
π
Ap = (0.112)2 = 9.852 × 10 −3 m 2
4
kg
hm = 10.4 Re 0.75
Dp
Cm D p ρ
Re =
μg
162
Capítulo 7. Pérdidas de calor y refrigeración
γ −1
⎡ ⎤
W c = m a (h2 R − h1 ) = m a
(h2 − h1 ) = m a ρ p (T − T1 ) =
m a C pT1 ⎢⎛ p2 ⎞ γ
⎜⎜ ⎟⎟ − 1⎥
ηc ηc 2
ηc ⎢ ⎝ p1 ⎠ ⎥
⎣⎢ ⎥⎦
p1V1γ = p 2V2γ
1−γ
γ γ γ γ
p1 ⎛ V2 ⎞ ⎛ nRT2 p1 ⎞ ⎛ p ⎞ ⎛T ⎞ T ⎛p ⎞ γ
= ⎜⎜ ⎟⎟ = ⎜⎜ ⋅ ⎟⎟ = ⎜⎜ 1 ⎟⎟ ⋅ ⎜⎜ 2 ⎟⎟ ⇒ 2 = ⎜⎜ 1 ⎟⎟
p 2 ⎝ V1 ⎠ ⎝ p 2 nRT1 ⎠ ⎝ p 2 ⎠ ⎝ T1 ⎠ T1 ⎝ p 2 ⎠
1−γ 1−1.33
⎛p ⎞ γ ⎛ 0,85bar ⎞ 1.33
T2 = T1 ⎜⎜ 1 ⎟⎟ = 293º K ⎜ ⎟ = 337,34º K
⎝ p2 ⎠ ⎝ 1,5bar ⎠
q real T2 − T2'
η int = =
q máx imo T2 − Tamb
1,5 273
ρ = 1,293 ⋅ ⋅ = 1,731kg / m 3
1,013 301,87
Suponiendo que en el radiador logremos tener un salto térmico de 10K, podríamos determinar el
caudal de agua necesaria para arrastrar este calor:
Q r 143kW
m H 2O = = = 3,421kg / s
C p ⋅ ΔT kJ
4,18 ⋅ 10 K
kg − K
163
Capítulo 7. Pérdidas de calor y refrigeración
Para las condiciones de funcionamiento que da este problema, se tiene un gasto másico de
combustible de 8,79 g/s de combustible, así pues:
Q ref = m agua ⋅ C p ⋅ ΔT
75,952 kg kg l
m agua = = 1,817 = 1,817 ♦
41,8 s s s
REFERENCIAS
[1] Ricardo, H., (1953) “The High Speed Internal Combustion Engine”, Blackie and sons,
London
[2] Heywood, J.B., (1988), “Internal Combustion Engine Fundamentals”, McGraw-Hill, New
York.
[3] G. Woschni, “A Universally Applicable Equation for the Instantaneous Heat Transfer
Coefficient in the Internal Combustion Engine”, SAE paper No. 670931 (1967)
[4] N. Watson and M.S. Janota, Turbocharging the internal combustion engine, The Macmillan
Press Ltd. London (1982)
[5] D.N. Assanis and J.B. Heywood, “Development and use of a Computer Simulation of the
Tubocompounded Diesel system for Engine Performance and Component Heat Transfer
Studies”, SAE paper No. 860329 (1986)
[6] W.J.D. Annand, “Heat Transfer in the Cylinders of Reciprocating Internal Combustion
Engines”, Proc. I. Mech. E, Vol. 177 No. 36 (1963)
[7] P. Flynn, M. Mizusawa, O.A. Uyehara and P.S. Myers, “An Experimental Determination of
the Instantaneous Potential Radiant Heat Transfer Within an Operating Diesel Engine”, SAE
paper No. 720022 (1972).
164
CAPÍTULO 8
INTRODUCCIÓN
Hasta ahora hemos considerado el rendimiento del motor en términos de la proporción del calor
total convertido en trabajo útil en el cilindro, pero ésta no es de ninguna manera la historia
completa; en lo que estamos realmente interesados es en la proporción de trabajo útil obtenible en
el del cigüeñal y en realidad no ganamos nada consumiendo por fricción interna o por bombeo de
aire lo que ganamos con algunos mejoramientos en el rendimiento térmico indicado, ya que
generalmente al tratar de aumentar ésta, aumentan los anteriores.
Tomada en su extensión la fricción interna mecánica de un motor es una función de la
máxima presión, puesto que esto, en gran parte, determina tanto el área de las partes deslizantes
como el peso de las partes móviles; también determina la fricción de los anillos del pistón, o por
lo menos del anillo superior del pistón contra las paredes del cilindro, pues para funcionar, es
esencial que este anillo tenga detrás de sí la presión del gas.
Así permaneciendo los demás factores iguales, mientras más baja sea la relación de la presión
máxima a la presión media efectiva, mayor será el rendimiento mecánica del motor. Además,
una gran cantidad de las pérdidas se debe a la viscosidad del aceite y otra parte sustancial, al
trabajo del bombeo de aire para llenar y vaciar el cilindro, ambos son independientes de la
presión de trabajo, por lo tanto se sigue también que, mientras más alta sea la presión efectiva
media, mayor será el rendimiento mecánico.
Desde el punto de vista de el rendimiento mecánico, necesitamos entonces la máxima pme
posible, combinada con la menor relación posible entre la presión máxima y la pme.
Desafortunadamente los recursos a los cuales podríamos acudir para mejorar el rendimiento
térmico indicado (tales como usar una relación de compresión muy alta o una prolongación del
tiempo de expansión) tienden a aumentar demasiado la relación entre la presión media y la
máxima, y una proporción sustancial de nuestras ganancias se perdería en la fricción entre el
cilindro y el cigüeñal.
De nuevo, a medida que aumentamos la velocidad de rotación del motor, ganamos en
rendimiento indicado debido a la reducción de las pérdidas de calor, pero igualmente, a medida
que aumentamos la velocidad aumentamos también la fricción resultante de las fuerzas dinámicas
Capítulo 8. Pérdidas mecánicas y lubricación
debidas a la inercia del pistón, etc. Las fuerzas dinámicas aumentan con el cuadrado de la
velocidad de rotación como se demostró en el Capítulo 6, como también la resistencia viscosa del
lubricante.
No todo el trabajo transferido por los gases al interior del cilindro - Trabajo indicado- se
convierte en trabajo disponible en el eje. Una fracción de este se transforma en fricción y a éste
se le denomina trabajo de fricción.
El trabajo o potencia de fricción es una fracción considerable del trabajo indicado, llegando a
ser cerca del 10% a plena carga y del 100% al ralentí, punto de funcionamiento en el cual toda la
energía generada en el cilindro se emplea en vencer las pérdidas de fricción. Gran parte de las
pérdidas de fricción se disipan en forma de calor y son arrastradas por el refrigerante y/o por el
aceite.
El trabajo de fricción, definido como la diferencia entre el trabajo durante la compresión y
expansión y el trabajo en el eje del motor se consume de la siguiente forma:
• En arrastrar mezcla fresca durante la carrera de admisión y en sacar los gases quemados
desde el cilindro durante la carrera de escape. Este usualmente se llama trabajo de
bombeo (durante la renovación de la carga).
• Para vencer el movimiento relativo de las partes móviles del motor. Incluye fricción entre
segmentos (anillos) y camisa, cigüeñal, rodamientos del árbol de levas, mecanismo de las
válvulas, piñones, poleas o bandas.
• Para mover los accesorios del motor. Ventilador, bomba de agua, bomba de aceite,
bomba de combustible, alternador, y en algunos casos aire acondicionado.
Trabajo de bombeo: Es el trabajo neto por ciclo hecho por el pistón durante la renovación de la
carga (carreras de admisión y escape). Wb se define únicamente para motores de 4 tiempos.
Trabajo de fricción por rozamiento: Es el trabajo por ciclo disipado para sobreponerse a la
fricción debida al movimiento relativo entre componentes adyacentes del motor.
Trabajo de accesorios o auxiliares: Es el trabajo por ciclo requerido para conducir los accesorios
del motor, como ventilador, bombas, alternador, compresor del aire acondicionado, etc.
Normalmente se incluyen en los cálculos solamente los más esenciales
166
Capítulo 8. Pérdidas mecánicas y lubricación
Wt = Wb + Wr + Wa (8.1)
Figura 8.1 Comparación entre las distintas formas de fricción del motor: pmf a diferentes cargas y
velocidades para un MEP y para un MEC [1]
La mejor forma de medir la fricción en un motor es por vía termodinámica, es decir, restando
de la potencia indicada (medida mediante un captador de presión en la cámara de combustión),
la potencia efectiva (o al freno). Este método sin embargo es difícil de usar en motores
multicilíndricos, debido a la variación entre cilindros de la potencia indicada y debido a la
dificultad de obtener datos de presión en cámara suficientemente precisos, ya que el ruido
inducido por el funcionamiento de los demás cilindros afecta las mediciones. Por esta razón, la
potencia de fricción usualmente se mide a motor arrastrado (sin combustión). Si bien es cierto
que las pérdidas son diferentes que cuando hay combustión, se suelen dar buenas aproximaciones
167
Capítulo 8. Pérdidas mecánicas y lubricación
NEXT I
Entre los puntos rca y aae (ciclo termodinámico cerrado) se almacena en una matriz de datos
los valores del trabajo indicado Wi dados como el producto pdV. Como recomendación de
algunos autores, se suma un 5.5% para acabar de cerrar el ciclo, ya que el rca y el aae no
coinciden con el PMI [2].
168
Capítulo 8. Pérdidas mecánicas y lubricación
Teniendo ahora la potencia indicada sólo basta restar la potencia efectiva para obtener
entonces la potencia de fricción global del motor.
2. Motor arrastrado. consisten en arrastrar el motor con un medio externo (otro MCIA o un
motor eléctrico), bajo condiciones de operación lo más parecidas posibles a cuando hay
combustión. Es necesario elevar primero la temperatura del aceite y del refrigerante del
motor. Otro método consiste en llevar el motor a las condiciones de operación normales y a
partir de allí retirar, rápidamente y por unos segundos, el sistema de encendido de aquel
cilindro donde se encuentre el captador de presión (la bujía en un MEP o la inyección de
combustible en un MEC). El resto de cilindros arrastran al cilindro en cuestión. Es importante
tener en cuenta que las pérdidas de fricción calculadas de esta manera incluyen las pérdidas
de bombeo.
169
Capítulo 8. Pérdidas mecánicas y lubricación
Figura 8.2 Método de las líneas de Williams para determinar la pmf [1]
170
Capítulo 8. Pérdidas mecánicas y lubricación
Donde pmf es la presión media de fricción (bar), pmax es la presión máxima de combustión
(bar) y Cm es la velocidad lineal media del pistón (m/s).
Winterbone y Tennant [4] han desarrollado su propia correlación basados en la ecuación de
Chen y Flynn con el fin de utilizarla durante la simulación del régimen transitorio del motor. Su
ecuación es la siguiente:
Ledger [5], Winterbone [6] y Watson [7] emplearon la misma estructura de la ecuación (8.3)
en sus modelos transitorios modificando únicamente los valores de los coeficientes de acuerdo a
sus resultados experimentales. Agudelo [8], obtuvo su propia correlación para un motor IVECO
diesel, turboalimentado, de 6 cilindros en línea, inyección directa y una cilindrada de 7,685 litros:
Donde pmf es la presión media de fricción (kPa), pmax es la presión máxima de combustión
(kPa) y nm es el régimen de giro del motor (rpm). En la Figura 8.3 se comparan los resultados
experimentales y los calculados por las diferentes correlaciones mencionadas de la presión media
de fricción. En la Figura 8.4 se muestra la correspondencia entre los datos experimentales y los
calculados con la correlación (8.5).
300
250
200
pmf (kPa)
150
100
50
0
800 1000 1200 1400 1600 1800 2000 2200
Régimen de giro (min-1)
Figura 8.3 Presión media de fricción experimental y calculada con diferentes correlaciones
171
Capítulo 8. Pérdidas mecánicas y lubricación
300
250
150
100
50
0
0 50 100 150 200 250 300
Calculados con Ecuación 2.35 (kPa)
Figura 8.4 Comparación de la pmf experimental y calculada con la correlación (8.5) [x]
Para MEP, Heywood [1] propone en su libro la siguiente correlación:
⎛ n(rpm ) ⎞
2
⎛ n ⎞
pmf total (bar ) = 0.97 + 0.15⎜ ⎟ + 0.05⎜ ⎟ (8.6)
⎝ 1000 ⎠ ⎝ 1000 ⎠
SISTEMA DE LUBRICACIÓN
El lubricante y el sistema de lubricación (Figura 8.5) desempeñan las siguientes funciones
principales:
• Reducir las pérdidas de fricción y asegurar el máximo rendimiento mecánico del motor
• Proteger el motor contra el desgaste
• Contribuir a la refrigeración del pistón y de aquellas partes por las cuales se disipa el
trabajo de fricción
• Remover impurezas de las zonas lubricadas
• Mantener las fugas de gas y de aceite (especialmente en la región de los anillos) a un
nivel mínimo aceptable.
172
Capítulo 8. Pérdidas mecánicas y lubricación
173
Capítulo 8. Pérdidas mecánicas y lubricación
REFERENCIAS
[1] Heywood, J.B., (1988), “Internal Combustion Engines Fundamental”, McGraw-Hill, New
York
[2] J.H. Horlock and D.E. Winterbone, (1986), “The Thermodynamics and Gas Dynamics of
Internal Combustion Engines”, Clarendon Press, Vol. 2, Oxford
1
SAE indica Society of Automotive Engineering
174
Capítulo 8. Pérdidas mecánicas y lubricación
[3] S.K. Chen and P.F. Flynn, “Development of a Single Cylinder Compression Ignition
Research Engine”, SAE paper No. 650733 (1965)
[4] D.E. Winterbone, D.W.H. Tennant, “The Variation of Friction and Combustion Rates During
Diesel Engine Transients”, SAE paper No. 810339 (1981)
[5] J.D. Ledger, N.D. Walmsley, "Computer Simulation of a Turbocharged Diesel Engine
Operating Under Transient Load Conditions", SAE paper No. 710177 (1971)
[6] D.E. Winterbone, C. Thiruarooran, P.E. Wellstead, "A Wholly Dynamic Model of a
Turbocharged Diesel Engine for Transfer Function Evaluation", SAE paper No. 770124 (1977)
[7] Watson N, Marzouk M., "A non-linear Digital Simulation of Turbocharged Diesel Engines
under Transient Conditions" SAE paper No. 770123 (1977).
175
CAPÍTULO 9
COMBUSTIBLES
Clasificación
La gran mayoría de los combustibles empleados en automoción son hidrocarburos. Se
caracterizan porque únicamente presentan carbono e hidrógeno en su estructura molecular. Se
clasifican principalmente de acuerdo con la presencia o no de enlaces múltiples o de anillos,
dando origen a las clases principales que se resumen en la Tabla 9.1.
Una clasificación pertinente respecto a su posible utilización en vehículos automotores puede
ser, según su estado en la naturaleza: sólido, líquido o gaseoso.
COMBUSTIBLES SÓLIDOS
Básicamente el carbono mineral. En los MCIA se han realizado especialmente en Estados
Unidos algunas investigaciones del quemado del “slurry” en motores. Éste es una mezcla pastosa
de carbón finamente pulverizado con agua (para mayor información, se recomienda buscar
aquellos números del Transactions of The ASME, Journal of Engineering for Gas Turbine and
Power dedicados a motores de combustión interna alternativos).
COMBUSTIBLES LÍQUIDOS
La mayoría proceden del petróleo. De éstos los más frecuentes son las gasolinas (combustibles
tradicionales para los MEP) y el ACPM (aceite combustible para motores tradicionalmente
empleado para los MEC). Los MEC de gran tamaño empleados en nuestra industria para
generación eléctrica, especialmente en horas pico, o para otras demandas de gran potencia,
pueden quemar combustibles llamados fuel oils, que son aceites de elevada viscosidad (se
recomienda consultar la norma ASTM D396 para combustóleos).
En la Tabla 9.2 se muestran las principales especificaciones de las gasolinas comerciales
(extra y corriente), según la norma Colombiana ICONTEC No. 1380. Igualmente, en la Tabla
9.3 aparecen las especificaciones según la norma Colombiana ICONTEC No. 1871 para la
gasolina de aviación tipo 100/130.
Capítulo 9. Combustibles y sistemas de suministro de combustible
Tabla 9.1
Clasificación de los hidrocarburos [1]
Particularidad de las
Clase Fórmula general Ejemplos
moléculas
Grupo Alquil
Parafínicos
Todos son llamados
compuestos n-normales e
Todos los enlaces simples. iso, respectivamente. CH4,
Alcanos CnH2n+2
No presentan anillos metano; C2H6, etano; C3H8,
propano; C8H18, n-octano e
iso-octano.
Cicloparafínicos o naftenos
C3H6, ciclo-propano (un
Todos los enlaces son anillo de 3 átomos de C);
Cíclicos simples. Anillos alifáticos, CnH2n C4H8, ciclo-butano (anillo
saturados o no. de 4 átomos de C); C5H10,
ciclo-pentano.
Olefinas
Tienen uno o más enlaces
C2H4, etileno; C3H6,
Alquenos dobles, por lo tanto son CnH2n
propileno; C4H8, butileno;
insaturados
Acetilenos
Uno o más enlaces triples.
El primer miembro es
Alquinos CnH2n-2 C2H2, acetileno.
acetileno,
H −C ≡ C − H .
Grupo de compuestos Aromáticos
Al menos un anillo de seis
carbonos con enlaces
dobles alternos. La
C7H8, tolueno; C8H10,
Aromáticos estructura de los CnH2n-6
xileno.
hidrocarburos aromáticos
es el anillo del benceno
(C6H6)
Grupo de los Alcoholes
En éstos compuestos
CH3OH, metanol; C2H5OH,
Alcoholes mohídricos orgánicos, un grupo CnH2n+1OH
etanol.
hidroxil (- OH)
177
Capítulo 9. Combustibles y sistemas de suministro de combustible
Tabla 9.2
Especificaciones de las gasolinas motor [2]
Especificaciones Norma ICONTEC No. 1380 Extra corriente
Número octano (mínimo) 92 80
Punto inicial de ebullición (ºC) 30 30
30 min 30 min
Porcentaje destilado a 100 ºC
70 máx 70 máx
Porcentaje destilado a 190 ºC (mínimo) 90 90
Punto final de ebullición (ºC) 225 225
Corrosión a 50 ºC (3 hrs, máximo) 1 1
Presión de vapor residual (máx, en kPa) 78.42 78.42
Contenido de “goma” en mg/100 cc máximo 5 5
Color Natural Natural
Azufre ( porcentaje en peso, máx) 0.15 0.15
Tabla 9.3
Especificaciones de la gasolina para aviación. Grado 100/130 [3]
Propiedad – Norma ICONTEC No. 1871 Valor
Valor antidetonante mínimo, número octano (mezcla pobre) 100
Valor antidetonante mínimo, número octano (mezcla rica) 130
Color Verde
Tetraetilplomo máximo en cc/dm3 1.0
Propiedades de destilación, temperatura máxima en ºC
¾ Para evaporar el 10% 75
¾ Para evaporar el 50% 105
¾ Para evaporar el 90% 135
Temperatura mínima para evaporar el 40% en ºC 75
Punto final de ebullición (máx en ºC) 170
Suma de temperatura para evaporar el 10%, 50%, (min en ºC) 135
Porcentaje de destilado recobrado (mínimo) 97
Porcentaje del residuo de destilación (máx) 1.5
Porcentaje de pérdida por destilación (máximo) 1.5
Calor neto de combustión en MJ/kg, mínimo 43.5
Presión de vapor en kPa 38 a 48
Corrosión a la lámina de cobre, 2 h a 100 ºC (en bomba), máx No. 1
Goma potencial (5 h de envejecimiento), máx en mg/100cc 6
Precipitado de plomo visible, mg/100cc máx 3
Azufre, porcentaje en masa, máx 0.05
Punto de congelación, máx (en ºC) -58
Interacción agua – combustible, el cambio de volumen no debe exceder del máx, en
2
cc
178
Capítulo 9. Combustibles y sistemas de suministro de combustible
COMBUSTIBLES GASEOSOS
Los tipos de combustibles gaseosos más comúnmente quemados en motores son el Gas Natural
Comprimido – GNC - y el GLP o Gas Licuado de Petróleo, que para automoción es generalmente
una mezcla de butano y propano en diferentes proporciones.
El Gas Natural debido a sus dificultades de licuefacción a temperatura ambiente se suministra
a lo motores comprimido (en torno a los 230 bar) lo que hace difícil su almacenamiento (los
depósitos deben cumplir normas muy exigentes – ASME 1 sec. 8, div. 1). La elevada presión de
almacenamiento permite almacenar suficiente combustible en el depósito para brindar autonomía
al vehículo. Puede haber recarga lenta o rápida, según el tipo de estación de llenado (la rápida
para un autobús tarda un mínimo de 20 minutos). Su componente fundamental es el metano (en
torno al 90-96%). Es un combustible ecológico hasta cierto punto ya que en su reacción de
combustión suele emitir tanto o más NOx que los otros combustibles. Su ventaja fundamental es
la homogeneidad de la mezcla por lo que produce muy pocos hidrocarburos sin quemar HC, sin
embargo, la emisión de metano al ambiente es mucho más severa desde el punto vista del efecto
invernadero que el CO2 (del orden de 30 veces más dañino).
Los vehículos que usan gas pueden ser de dedicación exclusiva, que tienen la ventaja de ser
diseñados con la relación de compresión adecuada para incrementar potencia, debido a las buenas
propiedades del GNC a la autodetonación. Estos sistemas permiten también emplear sistemas de
inyección electrónicos exclusivos para Gas Natural con una única centralita de control
electrónico lo que mejora enormemente el proceso de combustión en el cilindro. Otro tipo de
vehiculos que usan el GNC son los de funcionamiento dual (mediante el accionamiento de un
botón, una electroválvula deja pasar gas o gasolina, según convenga), contrario al anterior,
presentan la desventaja de pérdida de potencia y par cuando funcionan con Gas Natural debido a
que el motor inicialmente ha sido concebido (relación de compresión) para funcionar con
Gasolina Motor. Generalmente se necesita un avanzador electrónico de chispa (curva dual), para
lograr que su combustión sea lo más parecida posible a la de la Gasolina.
GLP. Su principal ventaja es que a presión y temperatura ambiente es líquido, lo cual facilita
enormemente su almacenamiento (bajas presiones). Su componente fundamental es el propano
aunque también contiene butano en forma de mezcla. En algunos países tienen un uso más
generalizado que el Gas Natural en automoción, debido principalmente a su similitud en el
proceso de combustión con la Gasolina. En cuanto a emisiones contaminantes es muy similar al
gas natural.
1
ASME indica American Society of Mechanical Engineering.
179
Capítulo 9. Combustibles y sistemas de suministro de combustible
Otros combustibles
Las políticas energéticas y ambientales actuales están propiciando la investigación de nuevos
combustibles alternativos tales como los alcoholes (metanol y etanol), los biocombustibles o
combustibles de origen vegetal (EMG , éster metílico de girasol; EMC, éster metílico de colza,
etc.).
Son de uso frecuente en países donde hay excedentes agrarios y se suelen mezclar
directamente con las gasolinas o con el ACPM sin hacer cambios en el motor.
Se habla del hidrógeno como el combustible del futuro ya que además de encontrarse en
abundancia en la naturaleza se puede emplear en pilas de combustible las cuales mediante un
proceso electroquímico (con muy bajas o casi sin emisiones contaminantes) transforman
directamente la energía química en fuerza electromotriz.
180
Capítulo 9. Combustibles y sistemas de suministro de combustible
2. Volatilidad
La volatilidad es una medida del porcentaje en volumen de combustible que se evapora en
función de la temperatura y se mide con la curva de destilación.
Los porcentajes más comunes en la curva de destilación son el 10%, 50%, y 90% de
combustible.
Si la volatilidad, cuando se ha destilado el 10% del combustible, es muy alta, hay problemas
de arranque en frío porque, cuando el motor opere a bajas temperaturas no se evapora el
suficiente combustible como para que se mezcle correctamente con el aire.
Si el punto del 10% es bajo, se pueden producir problemas por la formación de burbujas de
vapor.
181
Capítulo 9. Combustibles y sistemas de suministro de combustible
de combustible y aceite empleados en la operación del ensayo anterior. Existen varios métodos
para determinar el IO de un combustible. Los dos métodos más comúnmente usados son :
El código ASTM exige determinar el IO (MON y RON) en un motor CFR (Cooperative Fuel
Research). Este motor es monocilíndrico con relación de compresión variable (culata móvil).
El método motor es más exigente porque las condiciones de ensayo son más severas.
Algunos países como Alemania exigen la determinación del IO únicamente mediante el MON.
182
Capítulo 9. Combustibles y sistemas de suministro de combustible
Carburadores
Un carburador es un aparato que se emplea para controlar el flujo de combustible que entra al
sistema de admisión del motor, distribuyendo el combustible a través de la corriente de aire. El
aire fluye a través de una tobera convergente – divergente o venturi. La diferencia de presión
entre la entrada del carburador y la garganta (o estrangulamiento - la cual depende del flujo de
aire-) se emplea para medir la cantidad apropiada de flujo de combustible para ese flujo de aire.
El combustible entra a la corriente de aire a través del tubo de descarga de combustible (o puertos
en el cuerpo del carburador) y se atomiza para ser arrastrado a través de una válvula de mariposa
hacia el colector de admisión.
El combustible se empieza a evaporar dentro del carburador y continua en el colector a
medida que las gotas se mueven en la corriente de aire y a medida que fluye en forma de película
líquida sobre la válvula de mariposa y sobre las paredes del colector.
183
Capítulo 9. Combustibles y sistemas de suministro de combustible
184
Capítulo 9. Combustibles y sistemas de suministro de combustible
m f = C Do Ao 2 ρ f Δp f (9.4)
Donde CDo y Ao son el coeficiente de descarga y el área del orificio, respectivamente. Δpf es
la diferencia de presión a través del orificio, y el área del orificio se asume que es mucho menor
que el área del venturi.
Los coeficientes de descarga representan los efectos de todas las posibles desviaciones del
flujo ideal isentrópico unidimensional.
185
Capítulo 9. Combustibles y sistemas de suministro de combustible
Figura 9.2 Sistema de inyección electrónico multipunto. Medición de aire: sistema L-Jetronic de Bosch
[1]
Los sistemas tradicionales de inyección multipunto, tienen aún desventajas con la formación de
películas de combustible en los alrededores de las válvulas y del puerto de entrada, esto se discute
en el capítulo 12. Existe actualmente una tecnología más reciente en la cual los inyectores se
ubican en la parte inferior del puerto de admisión y no encima como comúnmente se hace. Esta
disposición permite que el combustible entre directamente al cilindro al momento de la apertura
de la válvula de admisión (Figura 9.4).
186
Capítulo 9. Combustibles y sistemas de suministro de combustible
Figura 9.4 Sistema de inyección multipunto con inyector montado en la parte inferior del puerto [5]
Sistemas de inyección monopunto
Consisten de uno o dos inyectores (ver Figura 9.5), controlados electrónicamente, que
suministran el combustible justo encima del cuerpo de la válvula de mariposa. Son más
económicos que los sistemas de inyección multipunto, pero tienen prácticamente las desventajas
asociadas al carburador desde el punto de vista de transporte de combustible
187
Capítulo 9. Combustibles y sistemas de suministro de combustible
Figura 9.6 Sistema de inyección con bomba de inyector único – unit injector –
188
Capítulo 9. Combustibles y sistemas de suministro de combustible
Para MEC que generan menos de 30 kW por cilindro (automoción de trabajo liviano,
automóviles) se emplean las bombas de inyección rotativas con presiones de inyección de hasta
75.000 kPa (750 bar). Para motores de tamaño medio (empleado en transporte pesado, camiones
y autobuses) que generan entre 40 y 100 kW por cilindro, se emplean las bombas de inyección en
línea (ver Figura 9.7). Éstas alcanzan presiones de suministro similares a las de inyector único.
En los mecanismos de las Figuras 9.6 y 9.7 el combustible es suministrado al inyector a
través de un orificio comunicado con el colector de combustible, el cual entra a ocupar la cavidad
justo por encima del pistón (plunger). Mediante un sistema de seguidor – leva que va entre el
pistón y el árbol de la bomba de inyección se empuja al pistón provocando el incremento de
presión del combustible. Cuando se vence la válvula de cheque ubicada a la salida de la bomba,
fluye el combustible hacia el inyector a través de las pipas o tubos de inyección. La cantidad de
combustible inyectada es función de la posición de la cremallera (rack) girando un engranaje
solidario al pistón, el cuál libera combustible por un canal labrado sobre él en forma helicoidal.
189
Capítulo 9. Combustibles y sistemas de suministro de combustible
LITERATURA RECOMENDADA
REFERENCIAS BIBLIOGRÁFICAS
[1] Heywood, J.B., (1988), “Internal Combustion Engine Fundamentals”, McGraw-Hill, New
York
[2] Norma Técnica Colombiana ICONTEC 1380, “Gasolina para Motores de Combustión
Interna”, 1992
[3] Norma Técnica Colombiana ICONTEC 1871, “Gasolina de aviación”, 1992
[4] M. Lapuerta, F. Losilla, O. Armas, (1997), “Incertidumbres Asociadas ala Utilización del
Poder Calorífico del Combustible en el Diagnóstico de la Combustión en Motores de Combustión
Interna”, Actas del XII Congreso Nacional de Ingeniería Mecánica, Vol. 2. Bilbao, España.
[5] Catálogo: Engine Videoscope AVL 513D. Visualisation of combustion and injection
precesses. 1998
[6] Garret, T. K. (1994), “Automotive Fuels and Fuel Systems”, Volume 1: Diesel, Pentech
Press, London
190
CAPÍTULO 10
Figura 10.1 Presión en el cilindro, fracción de masa quemada, y fracción de volumen quemado para
cinco ciclos consecutivos en un MEP en función del ángulo de cigüeñal. Salto de la chispa 30° APMS,
mariposa totalmente abierta, 1044 r.p.m., dosado = 0.98 [1]
Combustión Centrada
En una curva de presión - ángulo de cigüeñal (Figura 10.2) se puede observar como el par motor
(a un régimen y a unas condiciones en la admisión dadas) varía en función del instante en que
salte la chispa. Si se avanza la chispa (alejarla del PMS), el pistón transmite mayor cantidad de
192
Capítulo 10. Combustión en MEP
trabajo a los gases, lo que implica una presión pico muy elevada y muy próxima al PMS. Si por
el contrario la chispa se retrasa (acercarla al PMS), la presión pico se da muy avanzada la carrera
de expansión, reduciendo el trabajo transmitido por los gases al pistón y su magnitud será muy
baja respecto al caso anterior. El instante óptimo para que se produzca la chispa es aquel que
da el máximo par al freno (MBT), el cual ocurre en una posición intermedia a las dos
tendencias anteriores.
Figura 10.2 (a) Presión en el cilindro vs ángulo de cigüeñal para diferentes instantes de salto de la
chispa. (b) Efecto del avance de la chispa en el par al freno a relación Aire/Combustible y velocidad
constante, mariposa totalmente abierta. [1]
193
Capítulo 10. Combustión en MEP
Combustión anormal
El objetivo de esta sección es determinar el estado termodinámico de los gases en el interior del
cilindro, cuya distribución de presiones, temperaturas y densidades está variando debido a que la
combustión se produce en un proceso de propagación de una llama turbulenta confinada en un
espacio de geometría variable con el tiempo.
Durante la combustión, la presión en el cilindro aumenta gracias a la liberación de la energía
química de la mezcla aire/combustible. Cuando un elemento de mezcla aire/combustible arde, su
densidad se reduce aproximadamente la cuarta parte. Este elemento de gas ahora quemado,
inicialmente se expande y comprime de un lado la mezcla sin quemar que está delante del frente
de llama empujándola hacia las paredes de la cámara de combustión, y de otro lado a los
elementos de gas que ya han quemado comprimiéndolos hacia la bujía. Así, aquellos elementos
de la mezcla sin quemar que arden en instantes diferentes tienen temperaturas y presiones
diferentes justo antes de la combustión y por lo tanto se darán también estados termodinámicos
diferentes después de ésta. Esto implica que no uniformidad ni en composición ni en estado
1 Knock: nombre onomatopéyico en inglés asignado al ruido que caracteriza a la autoignición espontánea de aquella
porción de gas que es comprimida tanto por el frente de llama de turbulento como por el movimiento del pistón.
194
Capítulo 10. Combustión en MEP
Figura 10.3 Esquema del cilindro durante la combustión: gases sin quemar (U) a la izquierda del frente
de llama, gases quemados a la derecha. A representa el núcleo adiabático de gases quemados., BL
representa la capa límite en los gases quemados, W es la tasa de trabajo transmitido al pistón, Q es la
tasa de transferencia de calor a las paredes de la cámara de combustión [1]
Se toma como sistema termodinámico el cilindro, una vez que ha iniciado la combustión
(Figura 10.3). Asumiendo condiciones de operación normales y de momento cero fugas, la
presión en cualquier posición angular del cigüeñal a través del cilindro es aproximadamente
uniforme. Las condiciones en los gases quemados y sin quemar están determinadas por la
conservación de la masa:
VTot Vb Vu
= +
m m m
VTot m m
= υb ⋅ b + υu u
m m m
Donde m = mb + mu :
= υb ⋅ xb + (1 − xb ) ⋅ υu
VTot
m
x
b 1
VTot
= ∫ υb ⋅ dx + ∫ υu ⋅ dx (10.1)
m 0 xb
195
Capítulo 10. Combustión en MEP
ΔU = Q + W
U0 - U = Q + W
U = U0 - Q - W
Ub + Uu = U0 - Q - W
ub ⋅ mb + uu ⋅ mu = U 0 − W − Q
De donde:
U0 − W − Q
ub ⋅ xb + uu ⋅ (1 − xb ) =
m
U0 − W − Q
xb 1
∫u
0
b ⋅ dx + ∫ uu ⋅ dx =
xb
m
(10.2)
Se pueden obtener resultados ilustrativos suponiendo que los gases quemados y sin quemar con
sus calores específicos respectivos, son gases ideales diferentes; i.e.,
pυb = Rb Tb ub = Cv ,b Tb + h f ,b (10.3)
pυu = Ru Tu uu = Cv ,u Tu + h f ,u (10.4)
Asumiendo que :
x 1
1 b 1
xb ∫0 b 1 − xb x∫b u
Tb = ⋅ T dx Tu = ⋅ T dx
b x 1
pV RT RT
= p ∫ b b dx + p ∫ u u dx
m 0
p xb
p
196
Capítulo 10. Combustión en MEP
= Rb Tb xb + Ru Tu (1 − xb )
pV
(10.5)
m
pV R T R T
Tb = − u u + u u
mRb xb Rb xb Rb
Ru pV − mRu Tu
Tb = ⋅ Tu + (10.6)
Rb mRb xb
de (10.2) tenemos
U0 − W − Q
xb
∫ (C ) ( )
1
T + h f ,b ⋅ dx + ∫ Cv ,u Tu + h f ,u ⋅ dx =
v ,b b
0 xb
m
Integrando,
U0 − W − Q
( ) (
xb Cv ,b Tb + h f ,b + (1 − xb ) Cv ,u Tu + h f ,u = ) m
Cv 1
Reemplazamos Tb de (10.6) y además sabiendo que Cp - Cv = R, y que por lo tanto = :
R γ −1
⎡ ⎛R pV − mRu Tu ⎞ ⎤
U0 − W − Q
m
= xb ⎢Cv ,b ⎜ u Tu +
⎢⎣ ⎝ Rb mRb xb ⎠
(
⎟ + h f ,b ⎥ + Cv ,u Tu + h f ,u − xb Cv ,u Tu + h f ,u
⎥⎦
) ( )
Cv ,b ⎛ pV − mRu Tu ⎞ Cv ,b
U0 − W − Q
m
= xb Ru Tu
Rb ⎝
+⎜
m
⎟⋅
⎠ Rb
(
+ xb h f ,b + Cv ,u Tu + h f ,u − xb Cv ,u Tu + h f ,u ) ( )
Ru Tu ⎛ pV − mRu Tu ⎞
U0 − W − Q
m
= xb +⎜ (
⎟ + xb h f ,b + Cv ,u Tu + h f ,u − xb Cv ,u Tu + h f ,u
γ b − 1 ⎜⎝ m(γ b − 1) ⎟⎠
) ( )
U0 − W − Q ⎛R T ⎞ ⎛ pV − mRu Tu ⎞
m
(
= xb ⋅ ⎜ u u + h f ,b − Cv ,u Tu + h f ,u ⎟ + ⎜⎜
⎝γb −1
) ⎟ + Cv ,u Tu + h f ,u
⎠ ⎝ m(γ b − 1) ⎟⎠
( )
197
Capítulo 10. Combustión en MEP
U 0 − W − Q pV − mRu Tu
m
−
m(γ b − 1)
− Cv ,u Tu + h f ,u ( ) (−)
xb = ⋅
(−)
Ru Tu
(
+ h − Cv ,u Tu + h f ,u
γ b − 1 f ,b
)
pV − mRu Tu Q + W − U0
+ + Cv ,u Tu + h f ,u
m(γ b − 1) m
xb =
Ru Tu
h f ,u − h f ,b + Cv ,u Tu −
γb −1
xb =
pV − mRu Tu + (γ b − 1)( Q + W − U 0 ) + m(γ b − 1) Cv ,u Tu + h f ,u ( )
[ ( )
m (γ b − 1) h f ,u − h f ,b + (γ b − 1)Cv ,u Tu − Ru Tu ]
Reemplazando U0 = m[Cv,uT0 + hf,u] y p0V0 = mRuT0 se tendría :
xb =
(
pV − p0V0 + mRu T0 − mRu Tu + (γ b − 1)(Q + W ) − m(γ b − 1) Cv ,u T0 + h f ,u + m(γ b − 1) Cv ,u Tu + h f ,u ) ( )
[ ( ) (
m (γ b − 1) h f ,u − h f ,b + Tu γ b Cv ,u − Cv ,u − R )] u
xb =
( )
pV − p0V0 − mRu Tu − T0 + (γ b − 1)(Q + W ) + m(γ b − 1) Cv ,u Tu − Cv ,u T0 ( )
[ ( ) (
m (γ b − 1) h f ,u − h f ,b + Tu γ b Cv ,u − C p ,u )]
xb =
( )
pV − p0V0 − mRu Tu − T0 + (γ b − 1)(Q + W ) + m(γ b − 1) Cv ,u Tu − Cv ,u T0 ( )
[ ( )
m (γ b − 1) h f ,u − h f ,b + Tu Cv ,u (γ b − γ u ) ]
xb =
( )
pV − p0V0 + (γ b − 1)(Q + W ) − mRu Tu − T0 + mCv ,u (γ b − 1) Tu − T0 ( )
[ ( )
m (γ b − 1) h f ,u − h f ,b + (γ b − γ u )Cv ,u Tu ]
xb =
(
pV − p0V0 + (γ b − 1)(Q + W ) − m Tu − T0 γ b Cv ,u − Cv ,u − Ru )( )
[ ( )
m (γ b − 1) h f ,u − h f ,b + (γ b − γ u )Cv ,u Tu ]
198
Capítulo 10. Combustión en MEP
xb =
pV − p0V0 + (γ b − 1)(Q + W ) + (γ b − γ u )mCv ,u Tu − T0( )
[ ( )
m (γ b − 1) h f ,u − h f ,b + (γ b − γ u )Cv ,u Tu ] (10.7)
Una aproximación más simplificada consiste en suponer que los gases sin quemar son
comprimidos isentrópicamente, se puede decir entonces que :
(γ u −1)
Tu ⎛ p ⎞ γu
=⎜ ⎟ (10.8)
T0 ⎝ p0 ⎠
Esta ecuación, junto con las ecns. (10.7) y (10.8) permiten determinar xb y Tb a partir de las
propiedades termodinámicas de los gases quemados y sin quemar (Capítulo 4), y de los valores
de p, V, m y Q. Si por el contrario se conoce xb, se puede entonces determinar p. La fracción de
masa quemada y la presión en el cilindro están directamente relacionadas.
Se puede pensar en llegar a una ecuación para xb en función de la relación de densidades
entre la parte quemada y la parte sin quemar y de la relación entre el volumen total y el volumen
de los quemados (yb) de la siguiente manera :
mTot = mu + mb V = Vu + Vb yb = Vb/V
mb 1 1
xb = = =
mu + mb mu + mb mu
1+
mb mb
−1 −1
⎡ ρ ⎛V ⎞ ⎤ ⎡ ρ ⎛ V − Vb ⎞ ⎤
xb = ⎢1 + u ⎜ u ⎟ ⎥ = ⎢1 + u ⎜ ⎟⎥
⎣ ρb ⎝ Vb ⎠ ⎦ ⎣ ρb ⎝ Vb ⎠ ⎦
−1
⎡ ρ ⎛ 1 ⎞⎤
xb = ⎢1 + u ⎜ − 1⎟ ⎥ (10.9)
⎣ ρb ⎝ yb ⎠⎦
199
Capítulo 10. Combustión en MEP
Mediante la ecuación (10.8) se puede determinar como hemos visto la temperatura media de
los gases quemados, pero en la realidad la distribución de las temperaturas no es uniforme. La
mezcla que arde al iniciar la combustión es además comprimida a medida que el resto de la carga
se quema. Por otra parte, la mezcla que arde al final de la combustión ha sido previamente
comprimida por el frente de llama y por el movimiento del pistón y por lo tanto su estado final es
diferente. Existe un gradiente de temperaturas a través de los gases quemados, estando aquellos
que ardieron al inicio del proceso a una temperatura más alta. Esto se puede resolver a través de
dos modelos:
• modelo de mezcla total, donde se asume que cada elemento de mezcla que arde se mezcla
instantáneamente con el resto de los gases quemados y por tanto se tiene una temperatura
uniforme ( Tb ).
• modelo de no-mezclado, donde se asume que no existe mezcla entre elementos de gas que
arden en momentos diferentes, y que cada elemento de gas es por lo tanto comprimido (y
eventualmente expandido) isentrópicamente después de la combustión. Así:
γ b −1
Tb ( x , b , xb ) ⎡ p( x b ) ⎤ γb
=⎢ ⎥
Tb ( x , b ) ⎢⎣ p( x b ) ⎥⎦
,
h f ,u − h f ,b + C p ,u Tu ( xb, )
Tb ( x ,
b )= C p ,b
es la temperatura resultante de una combustión isentálpica de los gases sin quemar a Tu(x’b),
p(x’b).
En la Figura 10.4 se pude ver una distribución de temperaturas calculadas con este modelo.
La temperatura media es más cercana a la mas baja de las curvas. Se puede observar un
gradiente aproximadamente constante de unos 400K, siendo más alta la temperatura en el gas que
arde al comienzo de la combustión.
200
Capítulo 10. Combustión en MEP
Figura 10.4 p, Tb y xb durante la combustión. Los subíndices u y b indican sin quemar y quemado,
respectivamente, los subíndices ‘e’ y ‘l’ significan al comienzo y hacia el final, respectivamente [1]
201
Capítulo 10. Combustión en MEP
Figura 10.5 Temperatura de los gases de escape en función de la presión en el cilindro. Medida a través
de ventanas de cuarzo en la culata del motor, usando técnicas espectroscópicas. Temperaturas medidas
cerca de la bujía tienen valores más altos. Las líneas a trazos muestras el comportamiento isentrópico[1]
Resultados más precisos de las propiedades termodinámicas del gas se pueden obtener de la
siguiente manera:
= υb xb + υu (1 − xb )
V
(10.10)
m
U0 − W − Q
ub ⋅ xb + uu ⋅ (1 − xb ) = (10.11)
m
donde,
x x
1 b 1 b
xb ∫0 b 1 − xb ∫0 u
υb = υ dx y υu = υ dx
Igualmente para ub y uu. Para un dosado, un combustible y una fracción de gases quemados
determinados:
hu = hu(Tu) hb = hb(Tb,p) (10.12)
~ ⎞
⎛ R ⎛ R~ ⎞
pυ u = ⎜ ⎟T pυb = ⎜ ⎟T (10.13)
⎝ Mu ⎠ u ⎝ Mb ⎠ b
202
Capítulo 10. Combustión en MEP
Asumiendo que la transferencia de calor de los gases sin quemar es despreciable, entonces se
puede considerar una compresión isentrópica. Tu está especificada por algún estado inicial de los
gases sin quemar por P0, V0, Mu y la masa m. Entonces:
Tu = Tu ( p, s)
∂ Tu ∂T
d Tu = dp + u ds
∂p ∂s
⎛ ∂h ⎞
υ−⎜ ⎟
⎛ ∂Tu ⎞ ⎝ ∂p ⎠ T
⎜ ⎟ = (10.15)
⎝ ∂p ⎠ s ⎛ ∂h ⎞
⎜ ⎟
⎝ ∂T ⎠ p
⎛ ∂h ⎞
como para los gases sin quemar h = h(T) (sólo función de la temperatura) entonces ⎜ ⎟ = 0
⎝ ∂p ⎠ T
⎛ ∂T ⎞ dh ⎞
y por tanto ⎜ u ⎟ se convierte en una diferencial ordinaria siendo ⎟ = C p , reemplazando
⎝ ∂p ⎠ s dT ⎠ p
tendríamos:
υ
dTu = dp
C p (T )
se cumple que,
~ ⎞
⎛ R
pυ = ⎜ ⎟T
⎝ Mu ⎠ u
con lo cual;
~ ⎞
⎛ R
d Tu dp
∫ Cp (T) Tu
= ∫⎜ ⎟
⎝ Mu ⎠ p
Teniendo entonces datos de la curva de presión en cámara de combustión, se puede determinar
Tu y luego despejando de (10.10) y (10.11) a xb e igualando ambos tendríamos
⎛V ⎞ ⎛U ⎞
⎜ ⎟ − υu ⎜ ⎟ − uu
⎝ m⎠ ⎝ m⎠
= (10.16)
υb − υu ub − uu
203
Capítulo 10. Combustión en MEP
Esta ecuación se puede resolver por una técnica iterativa apropiada para determinar Tb
La presión en el cilindro cambia con el ángulo del cigüeñal como resultado de:
• variación del volumen debido al movimiento del pistón;
• combustión;
• transferencia de calor a las paredes de la cámara de combustión;
• flujos de gases que entran y salen de las fisuras que quedan dentro de la cámara y
• fugas a través de los segmentos del pistón.
La Figura 10.6 muestra una cadena de adquisición de datos de presión de la que se puede obtener
después mediante la ecn. (10.17) una relación entre la presión y el volumen instantáneo del
cilindro [2].
( )
⎡ ⎛ Rc − 1⎞ ⎤ Vc = volumen muerto
V (θ ) = Vc ⋅ ⎢1 + ⎜ ⎟ ⋅ R + 1 − cosθ − R 2 − sen θ 2 ⎥ (10.17)
⎣ ⎝ 2 ⎠
Rc = relación de compresión
⎦
R = relación Biela-Manivela
θ = ángulo de cigüeñal
Figura 10.6 Montaje experimental para análisis de datos de la presión dinámica en cámara de
combustión [2]
204
Capítulo 10. Combustión en MEP
Igualmente en la referencia [2], se ha hecho un análisis de los parámetros claves que se deben
considerar cuando se trata de analizar datos de presión en cámara de combustión, algunos de
ellos son:
• determinación adecuada de la presión de referencia
• precisión del codificador angular
• determinación correcta del PMS
• variaciones en la señal del captador de presión piezoeléctrico debido a los gradientes térmicos,
etc
En la curva logp-logV de la Figura 10.7 el proceso de compresión es una línea recta cuya
pendiente (exponente politrópico) es 1,3. Se puede identificar el ángulo de inicio de la
combustión por el momento en que se separan la línea recta y la curva.
El ángulo de final de la combustión se puede ubicar con el mismo criterio anterior. El
proceso de expansión es una línea recta cuya pendiente (n) es 1,33. La diferencia entre los
exponentes politrópicos durante la compresión y la expansión se debe únicamente a las pérdidas
de calor por refrigeración. Para combustibles convencionales n = 1,33 ± 0,05.
Figura 10.7 (a) Diagrama presión - volumen ; (b) logp - log(V/Vmax). 1500r.p.m., MBT timinig, pme =
513 kPa. F = 0,8, Rc = 8,72 [1]
205
Capítulo 10. Combustión en MEP
Una técnica para estimar la fracción de masa quemada xb a partir de los datos de presión en el
cilindro, es la desarrollada por Rassweiler y Withrow, según la cuál en cualquier intervalo de
ángulo de cigüeñal Δθ, el cambio en la presión real Δp es la suma de el aumento presión debido
a la combustión Δpc y el cambio en la presión debido al cambio en el volumen Δpv.
Las presiones y los volúmenes al inicio y al final del intervalo Δθ, en ausencia de
combustión, están relacionadas por:
piVi n = p jV jn
piVi n
Δpv = p j − pi = − pi
V jn
⎡⎛ V ⎞ n ⎤
Δpv = pi ⎢⎜⎜ i ⎟⎟ − 1⎥ (10.19)
⎢⎝ V j ⎠ ⎥
⎣ ⎦
mb ( i ) ∑ Δp c
=
0
N (10.20)
mb ( total )
∑ Δp
0
c
206
Capítulo 10. Combustión en MEP
Otra vía para analizar los datos de presión es aplicando la primera ley de la termodinámica
aplicada a un sistema abierto (Figura 10.8). La principal ventaja de este método es que relaciona
el aumento real de la presión con la cantidad de energía química contenida en el combustible.
Dependiendo de lo fino que se quiera ser, se puede trabajar con un modelo de una sola zona, de
dos o incluso de tres zonas si se consideran los residuales.
Figura 10.8 Límites del sistema abierto para el análisis de liberación de calor en la cámara de
combustión. [1]
U = m·u(T) donde m es la masa dentro de los límites del sistema y T es la temperatura media
de la carga dentro del sistema
dU = m·du(T) + u(T)·dm
207
Capítulo 10. Combustión en MEP
Reemplazando en (10.21):
δQCH = mCvdT - udmcr + δQHT + pdV + h’dmcr
pdV Vdp
mdT = + − Tdm
R R
⎛ Cv ⎞ ⎛ Cv ⎞
δQCH = ⎜ ⎟Vdp + ⎜ + 1⎟ pdV + (h ′ − u + Cv T )dmcr + δQHT (10.23)
⎝ R⎠ ⎝ R ⎠
⎛ Cv ⎞ ⎛ Cv ⎞
δQCH − δQHT − (h ′ − u + Cv T )dmcr = ⎜ ⎟Vdp + ⎜ + 1⎟ pdV
⎝ R⎠ ⎝ R ⎠
Cuando a la ecuación del ‘calor liberado neto’ se le normaliza para que su máximo valor sea
la unidad y se representa gráficamente contra el ángulo de cigüeñal, se obtiene la curva de
fracción de masa quemada (que en términos más precisos es la fracción de energía liberada).
208
Capítulo 10. Combustión en MEP
Un modelo para considerar el flujo de energía de los volúmenes de las fisuras y grietas es el
siguiente:
⎡T′ 1 ⎛ γ − 1 ⎞ ⎤ dp
(h ′ − u + C T )dm = Vcr ⎢ +
T
+ ln⎜ ⎟⎥
⎢⎣ Tw Tw ( γ − 1) bTw ⎝ γ ′ − 1⎠ ⎥⎦ dθ
v cr
Un ejemplo del uso de esta ecuación para analizar los datos experimentales de presión contra
ángulo de cigüeñal para un MEP convencional se muestra en la Figura 10.9
Figura 10.9 Resultados del análisis del calor liberado. Se pueden observar los efectos del calor
transmitido, las fisuras, y la ineficacia de la combustión
209
Capítulo 10. Combustión en MEP
Ángulo de quemado rápido Δθb. Intervalo de grados de cigüeñal requerido para quemar el
grueso de la carga. Se define como el intervalo entre el final de la etapa de desarrollo de la llama
(Usualmente 10% de la fracción de masa quemada o energía liberada) y el final del proceso de
propagación de la llama (normalmente el 90% de la fracción de masa quemada o energía
liberada).
Ángulo de quemado total Δθ0. Es la duración total del proceso de combustión. Es la suma de
los dos anteriores.
210
Capítulo 10. Combustión en MEP
Figura 10.10 Definición del ángulo de desarrollo de la llama, y del ángulo de quemado rápido, en una
curva de fracción de masa quemada contra ángulo de cigüeñal
Una manera usual de representar la curva de fracción de masa quemada contra el ángulo de
cigüeñal es la función de Wiebe:
θ es el ángulo de cigüeñal
⎡ ⎛ θ − θ0 ⎞ m +1
⎤ θ0 es el ángulo de inicio de la combustión
xb = 1 − exp ⎢− a ⎜ ⎟ ⎥ (10.25)
⎣ ⎝ Δθ ⎠ ⎦ Δθ es la duración de la combustión total
a y m son parámetros ajustables, a = 5 y m = 2
211
Capítulo 10. Combustión en MEP
intensidad depende de las características del combustible y de las características anti-knock del
motor.
• La detonación, establece que el frente de llama se acelera a la velocidad del sonido (se genera
una onda de choque) y consume todo el gas terminal a una tasa mucho más elevada de lo que
ocurriría con un frente de llama que se propagase a velocidad normal.
• El oído humano
• pruebas ópticas
• detectores de ionización y
• captadores de intensidad de knock (Figura 10.11)
Figura 10.11 Presión en el cilindro contra ángulo de cigüeñal con (a) combustión normal, (b) knock
suave, y c) knock fuerte. 4000 r.p.m., mariposa totalmente abierta, 381 cm3. Motor monocilíndrico [1]
212
Capítulo 10. Combustión en MEP
213
Capítulo 10. Combustión en MEP
REFERENCIAS
[1] Heywood, J.B. “Internal Combustion Engine Fundamentals”, McGraw-Hill, New York,
1988, 930p.
[3] Woschni, G. “A Universally Applicable Equation for the Instantaneous heat Transfer
Coefficient in the Internal Combustion Engine”, SAE paper No. 670931. 1967.
214
CAPÍTULO 11
INTRODUCCIÓN
Tal y como lo expresaban N. Watson y M.S. Janota [1] en el año 1982: “El proceso de
combustión es el aspecto más importante de cualquier motor de combustión interna.
Desafortunadamente es también el más complejo y el menos entendido”. Con el paso de los años
se ha podido profundizar más en el estudio del proceso de combustión en motores, gracias al
avance alcanzado en tecnologías de visualización e indicación en el interior del cilindro.
El proceso de combustión en un motor de encendido por compresión, además de ser no
estacionario, heterogéneo y multidimensional ocurre en muy poco tiempo, del orden de
milisegundos, en un medio confinado y con geometría variable, lo cuál por lo demás confirma lo
expuesto en el párrafo anterior.
1) Tiempo de retraso (a-b). Es el tiempo (o ángulo) que transcurre desde el momento en que
se inyecta el combustible dentro de la cámara de combustión, hasta el momento en que inicia
la combustión. Este último punto es muy difícil de determinar con exactitud y por lo tanto
existen diversos criterios en la literatura [3].
Tiempo de retraso
Existen dos vías para calcular el tiempo de retraso: la primera consiste en emplear correlaciones
semiempíricas como la de Hardenberg y Hase [4] ampliamente aceptada en el campo de los
motores. Esta expresión y otras como la de Wolfer [2] igualmente reconocida, emplean un
arreglo del tipo Arrehnius, en el cual la presión y la temperatura durante el tiempo de retraso
quedan dentro de la parte exponencial.
216
Capítulo 11. Combustión en MEC
217
Capítulo 11. Combustión en MEC
Figura 11.2 Curva típica de liberación de calor aparente de un motor diesel de inyección directa
identificando las 3 fases de la combustión [2].
En la Tabla 11.1 se hace una comparación entre los resultados calculados con la ecuación
(11.1) y los resultados obtenidos experimentalmente en un motor IVECO diesel de inyección
directa turboalimentado [2].
⎡ (0 . 36 + 0 . 22 C m )⎤ ⎡ ⎛ 1 1 ⎞ ⎛ 21 . 2 ⎞
0 . 63 ⎤
t ret = ⎢ ⎥ exp ⎢ E A ⎜ ~ − 17190 ⎟ ⎜⎜ p − 12 . 4 ⎟⎟ ⎥ (11.1)
⎣ 0 . 006 n m ⎦ ⎢⎣ ⎝ RT ⎠⎝ ⎠ ⎥⎦
donde:
Cm = velocidad lineal media del pistón (m/s)
~
R = constante universal de los gases ideales (8.3143 J/mol·K)
618840 es la energía de activación aparente (J/mol)
EA =
CN + 25
CN = índice de cetano del combustible
nm = régimen de giro del motor (rpm)
218
Capítulo 11. Combustión en MEC
m(τ )
M (τ ) =
mt
[
= 1 − exp − k 2τ (k1 +1) ] (11.2)
[
M (τ ) = (k + 1)k 2τ k1 exp − k 2τ ( K 1 +1) ] (11.3)
Donde:
(θ − θ i )
τ= (11.4)
Δθ
219
Capítulo 11. Combustión en MEC
ecuación (11.2) lo ponen como función del tiempo de retraso (tr), el régimen de giro del motor
(nm), la presión (pc) y la temperatura (Tc) en el inicio de la compresión:
a b
⎡ t r ⎤ ⎡ p ⎤ ⎡T ⎤ ⎡ n ⎤
K1 = k10 ⎢ 0 ⎥ ·⎢ c ⎥·⎢ c 0 ⎥·⎢ m 0 ⎥ (11.5)
⎣ t r ⎦ ⎣ pc 0 ⎦ ⎣ Tc ⎦ ⎣ nm ⎦
donde el subíndice 0 representa las condiciones de referencia basada en una curva conocida y los
exponentes a y b son constantes empíricas.
La Figura 11.2 mostraba una curva de liberación de calor aparente típica de un motor diesel
de inyección directa, en ella se observan dos picos claramente definidos y que corresponden a las
fases de combustión premezclada y combustión por difusión turbulenta. Ni la función de Wiebe
ni la de Woschni y Anisits pueden representar estas dos fases de combustión con precisión. Por
esta razón en 1980 Watson, Pilley y Marzouk [6] propusieron una nueva correlación de la
combustión que comprende dos expresiones algebraicas, una para la fase de premezcla y otra
para la fase de difusión:
m t = m p + m d (11.6)
mp
β = (11.7)
mt
donde mt = m p + md (11.8)
[
M p (τ ) = 1 − 1 − τ ]
c p1 c p 2
(Acumulado) (11.9)
220
Capítulo 11. Combustión en MEC
(c −1)
[
M p (τ ) = C p1 ·C p 2 ·τ p1 · 1 − τ p1
c
](
c p 2 −1)
(Tasa) (11.10)
[
M d (τ ) = 1 − exp − cd 1 ·τ c d 2 ] (Acumulado) (11.11)
[
M d (τ ) = cd 1 ·cd 2 ·τ (c d 2 −1) ·exp − cd 1 ·τ c d 2 ] (Tasa) (11.12)
Reuniendo términos:
donde:
Fr b
β =1− a (11.14)
t rc
Fr = dosado relativo
tr = tiempo de retraso (ms)
Cp1 = 2.0 + 1.25 · 10-8(tr · nm)2.4
Cp2 = 5000
−k 2
Cd1 = k1 · Fr
C d 2 = k3 ·Cd 1k 4
nm = régimen de giro del motor en rpm
k1, k2, k3 y k4 son coeficientes que permiten ser ajustados para diferentes condiciones de
operación incluyendo el régimen transitorio y diferentes tipos de motores diesel.
A mediados de la década de los ochenta Miyamoto et al [7] (1985), reportaron una correlación
similar basada en dos funciones de Wiebe con la siguiente forma:
221
Capítulo 11. Combustión en MEC
⎡ (M +1)
⎛θ ⎞ p ⎤
Mp
⎛θ ⎞
dQ
dθ
= 6.9
Qp
θp
( )
M p +1 ⎜ ⎟
⎜θ p ⎟
exp⎢− 6.9⎜ ⎟
⎢ ⎜θ p ⎟
⎥+
⎥
⎝ ⎠ ⎣ ⎝ ⎠ ⎦ (11.15)
⎡ (M +1)
⎛θ ⎞ d ⎤
Md
Qd ⎛θ ⎞
6.9 (M d + 1)⎜⎜ ⎟⎟ exp⎢− 6.9⎜⎜ ⎟⎟ ⎥
θd ⎝θd ⎠ ⎢⎣ ⎝θd ⎠ ⎥⎦
Donde los subíndices p y d se refieren a las fases de combustión premezclada y por difusión
respectivamente. Qp y Qd son las cantidades de calor liberado, θp y θd son las duraciones
de la liberación de calor, y Mp y Md son los factores de forma. Estos seis parámetros son
ajustables al igual que en la correlación anterior de acuerdo con las condiciones de operación y el
tipo de motor.
En el presente trabajo se ha adoptado la ley de Watson para calcular la tasa de quemado de
combustible ya que es ampliamente usada en el campo de la simulación de la combustión en
motores diesel, además requiere ajustar un parámetro (coeficiente) menos que la correlación
propuesta por Miyamoto et al. [7].
En la Tabla 11.2 se han ilustrado los valores de los coeficientes de la ley de Watson obtenidos
por diferentes investigadores. En este trabajo se han ajustado las constantes durante la fase de
calibración del modelo para el motor IVECO 8360.46.417 usado en los ensayos [8].
Mediante la ley de Watson se obtiene, como se demostró antes, la tasa de combustible
f ) . Este dato es necesario en el balance termodinámico que se realiza dentro del
quemado ( m
cilindro para obtener la tasa de liberación de calor aparente. La manera de relacionarlo es a través
del poder calorífico másico inferior del combustible mediante la ecuación (11.16):
222
Capítulo 11. Combustión en MEC
Tabla 11.2
Valores de las constantes de la ley de Watson para diferentes tipos de motores diesel [8]
Tipo de Motor
Constantes
Empíricas Cummins
A* B* C* IVECO[5] Ford [6]
[11]
A 0.926 0.95 0.81 1.05 0.926 1.0569
B 0.37 0.41 0.28 0.36 0.37 0.489
C 0.26 0.28 0.51 0.25 0.26 0.55
k1 14.2 16.67 7.54 3.05 14.2 79.705
k2 0.644 0.6 0.65 0.644 0.644 1.484
k3 0.79 1.2 0.93 1 1.05 0.581
k4 0.25 0.13 0.22 0.25 0.25 0.739
Motor A: Seis cilindros en línea, turboalimentado, cuatro tiempos, inyección directa, con
puerto de generación de swirl y cámara de combustión tipo “deep bowl”
Motor B: V8, turboalimentado con intercooler, diesel de inyección directa, cuatro tiempos.
Motor C: Turboalimentado con intercooler pero de mayor potencia que los motores A y B.
Los anteriores modelo llevan a predecir el proceso de combustión al interior del cilindro. Los
parámetros más relevantes que se pueden llegar a obtener mediante esta simulación son:
En la Figura 11.3 se puede observar un chorro de combustible inyectado al interior del cilindro
de un MEC, visualizado mediante una ventana de cuarzo con luz estroboscópica y con captura
digital de imágenes mediante el equipo engine videoscope AVL 513D [10]. La imagen de la
derecha muestra un análisis de temperatura en un plano transversal del cilindro. Obsérvese que
para las condiciones de funcionamiento del motor en las que fueron filmadas estas imágenes, se
alcanzan en el interior del cilindro temperaturas máximas (zonas rojas) del orden de 2275 K.
223
Capítulo 11. Combustión en MEC
224
Capítulo 11. Combustión en MEC
INICIO
RCA a SI SI
AAI ALMACENAR VALORES
? PASO ANTERIOR
QL = 0
CÁLCULO QL
i < AAE NO
?
FIN
225
Capítulo 11. Combustión en MEC
226
Capítulo 11. Combustión en MEC
REFERENCIAS
[1] Watson, N., Janota, M.S. (1982), “Turbocharging the internal combustion engine”, The
Macmillan Press Ltd. London
[2] Heywood, J.B., (1988) “Internal combustion engine fundamentals”, McGraw-Hill, New York
[3] Agudelo, J.R., Casanova, J., (1998), “Combustion Analysis Based on Low Cost Personal
Computer and Related Instruments”, Proceedings of the international simposium of
Autotest’98, España
[4] Hardenberg, H.O. and Hase F.W., (1979), “An Empirical Formula for Computing the
Pressure Rise Delay of a Fuel from its Cetane Number and from the Relevant Parameters of
Direct-Injection Diesel Engines”, SAE paper No. 790493
[5] Woschni, G., Anisits, F., (1974), “Experimental Investigation and Mathematical Presentation
of Rate of Heat Release in Diesel Engines Dependent upon Engine Operating Conditions”,
SAE paper No 740086
[6] Watson, N., Pilley, A.D., and Marzouk, M., (1980) “A Combustion Correlation for Diesel
Engines Simulation”, SAE paper No 800029
[7] N.Miyamoto, Takemi Chikahisa, Tadashi Murayama And Robert Sawyer, “Description and
Analysis of Diesel Engine Rate of Combustion and Performance Using Wiebe’s
Functions”, SAE paper No 850107 (1985).
[8] Agudelo J.R., “Análisis Teórico – Experimental de Motores de Encendido por Compresión
Turboalimentados en Régimen Transitorio”, Tesis Doctoral, Universidad Politécnica de
Madrid, España. 1998.
[9] G. Woschni, “A Universally Applicable Equation for the Instantaneous Heat Transfer
Coefficient in the Internal Combustion Engine”, SAE paper No 670931 (1967).
[10] Catálogo: Engine Videoscope AVL 513D. Visualisation of combustion and injection
precesses. 1998
227
CAPÍTULO 12
EMISIONES CONTAMINANTES
INTRODUCCIÓN
A lo largo de la histortia de la humanidad se ha destacado la importancia del aire y sobre todo de los
efectos que tiene respirar aire contaminado, pues a causa de este hecho han habido numerosas
muertes. En 1943 en Inglaterra murieron alrededor de 4000 personas por causa del smog [3] (smog
proviene de la fusión de las palabras inglesas smoke = humo y fog = niebla [2] ). Estos son los
efectos más funestos que sabemos se pueden esperar de la contaminación del aire, pero hay otros
efectos que nos tocan de cerca como el que narra Eilliam H. Crouse [2]: “...Si has esperimentado un
caso desfavorable de smog, nunca lo olvidarás. No puedes ver a lo lejos, tienes difcultad para
respirar. Y el aire que respiras no se siente bien. Tus ojos y tu garganta se irritan. Los alimentos y
las plantas se ven seriamente dañados. La pintura de tu casa se ve afectada...”.
Bastantes enfermadades respiratorias han sido causadas por la contaminación del aire en los
seres humanos y en los animales, pero no podemos olvidar otras implicaciones negativas de este
hacho como son el envenenamiento de lagos en Estados Unidos de América, la muerte de los
bosques en Alemania y el deterioro de los monumentos y las edificaciones en ciudades como
Londres y Madrid.
Los primeros vestigios de control en las emisiones contaminantes tienen origen en Inglaterra a
principios de siglo XIII cuando la corona decidió cobrar impuestos a quienes utilizaran carbón como
medio de calefacción debido a la gran concentración de humo que éste producía [3]. A principios de
1910 en Alemania se empieza a hablar muy tímidamente del control de las emisiones provenientes de
los motores de los vehículos [3], pero es solo hasta principios de la década de los sesenta en la ciudad
de los Angeles (California) cuando se empiezan a regular por vía legal. La razón por la cual se inició
allí todo el movimiento que restringía las emisiones contaminantes se debió a que la carencia casi
total de industria y los pocos calefactores existentes en esta zona, llevaron a concluir que eran los
automóviles los causantes de un 60% del total de la emisión de contaminantes [14]. Desde entonces
todos los países industrializados han adoptado sus métodos para restringir las emisiones provenientes
de los vehículos, siendo cada vez más estrechos los límites de permisividad, lo que ha llevado a los
fabricantes de motores a enfrentar grandes retos tecnológicos que les permitan cada vez ajustarse más
a la legislación anticontaminante con el fin de poder sobrevivir en el mercado.
Capítulo 12. Emisiones Contaminantes
Así pues, se ha logrado incluir el estudio del nivel de emisiones como uno de los parámetros
fundamentales en la evaluación del desempeño de los motores [1].
En el presente Capítulo se empieza haciendo una descripción de las fuentes de emisiones
contaminantes en el vehículo y los principales parámetros que influyen para que éstas se produzcan,
se hace también una descripción de como se forman los diversos contaminantes en cada tipo de
motor y finalmente se establece una comparación por tipo de combustible en la aportación de
emisiones contaminantes. Luego se describen las características que deben poseer la muestra a ser
analizada proveniente de los gases de escape de los MCIA y la diferentes técnicas de muestreo
utilizadas en la medición de emisiones contaminantes provenientes de los motores, más adelante se
hace una revisión de los principios analíticos empleados para su determinación y finalmente se toman
uno a uno y se describen cuales son los métodos más ampliamente utilizados en el ensayo de
motores.
229
Capítulo 12. Emisiones Contaminantes
los MEC, dado que la principal fuente de contaminación de los motores proviene de los gases de
escape productos de la combustión. Los MEC emiten partículas carbonáceas que se pueden
considerar virtualmente ausentes en los motores que queman gasolina (MEP).
El dióxido de carbono (CO2) es un producto normal de la combustión de cualquier combustible
que contenga carbón, no es tóxico y solo se considera contaminante si supera la concentración
normal en la naturaleza (330ppm), ya que desplaza al aire al ser 1,5 veces más denso que él,
haciéndole irrespirable [16].
b. Gases procedentes del cárter del motor, que fluyen al exterior por el respiradero de éste. Estan
compuestos fundamentalmente por HC, aunque, dependiendo del estado del motor, pueden
contener también productos procedentes de la combustión. Son responsables de otro 25% del
total de HC emitidos por el vehículo. Esta fuente de emisiones puede eliminarse fácilmente
ediante el recirculado de estos vapores hacia la admisión (blow-by) [2].
Las cantidades de contaminantes varían con el tipo de motor, su diseño, su geometría interna, sus
condiciones de funcionamiento, el tipo de combustible y sus aditivos, el tipo de aceite lubricante, la
forma y disposición de la cámara de combustión, el sistema de suministro de aire, el sistema de
suministro de combustible, entre otros [8] (en esta referencia se analizan las mejoras realizadas
230
Capítulo 12. Emisiones Contaminantes
sobre el diseño de los MEC hasta el año 1991), además existen otros parámetros que influyen en las
emisiones del automóvil. Los órdenes de magnitud de las cantidades emitidas están dados en la
Tabla 12.1 [6], otros contaminantes adicionales son el bromuro de plomo y el óxido de azufre, los
cuales varían directamente con el contenido de plomo y de azufre respectivamente de los
combustibles del motor, y los aldehídos, los cuales resultan del uso de alcoholes como combustible.
Tabla 12.1
Orden promedio de magnitud de contaminantes emitidos [6].
Cantidad relativa emitida
Contaminante Concentración (ppm)
(g/kg de combustible)
NOx 500 a 1.000 20
CO 10.000 a 20.000 200
HC (ppm en base C) 3.000 25
Partículas --- 2a5
La variable más importante que gobierna la emisión de contaminantes en los MCIA, independiente
de que sea MEP o MEC es la relación aire - combustible o su inversa combustible - aire (dosado), la
cual afecta directamente las emisiones relativas de CO, NOx, aldehídos y HC (Figura 12.1 [3]).
Nótese que para un dosado estequiométrico se tendrían las bajas emisiones de CO y HC, pero las
emisiones de NOx estarían cercanas a su máximo, esto se debe entre otras cosas a la formación de
NOx a altas temperaturas. Para estudiar detalladamente los procesos de formación de estos
contaminantes durante el proceso de combustión que tiene lugar en los MCIA, se recomienda mirar
en [3] Capítulo 10; en [6] Capítulo 11, en [9] Capítulo 18 y en [14] numeral 4.
Las emisiones provenientes del vehículo no dependen únicamente del tipo de motor y de su potencia
como podríamos pensar hasta el momento, también dependen de otros factores como son: pendiente
de la carretera, altitud, humedad, temperatura del ambiente, límites de velocidad, desgaste del
vehículo, condiciones de tráfico, entre muchas otras. Estas condiciones indican que las emisiones
reales en "tráfico", relacionadas con la distancia viajada por el vehículo, no siempre coinciden con
las medidas reguladoras tomadas en condiciones del ciclo de conducción claramente definido.
231
Capítulo 12. Emisiones Contaminantes
Hidrocarburos (HC)
Los gases de escape de los MEP contienen entre 1,000 a 3,000 ppmC, lo que corresponde
aproximadamente de un 1 a 2.5% del combustible alimentado al motor [6]. En la Figura 12.1 se veía
la dependencia de la cantidad de HC sin quemar con la relación aire - combustible. Para mezcla ricas
y mezcla pobres, aumenta dicha cantidad, debido a que el proceso de combustión no se desarrolla en
buenas condiciones y se puede dar el apagado de la llama. Para mezcla ligeramente pobres, donde
existe exceso de oxígeno y la temperatura todavía es elevada (es decir, que la llama no tiene
problemas para progresar), la aparición de HC en el escape es mínima, puesto que éstos se oxidan,
aunque solo sea parcialmente (formación de CO), en el colector y tubo de escape.
232
Capítulo 12. Emisiones Contaminantes
Las emisiones de HC en un MEP son producidas al cesar las reacciones de combustión en las
proximidades de las paredes (efecto pared o wall quenching), por cortocircuito de la carga fresca y
por combustión incompleta (flame quenching) [9] de la carga, si bien esta última causa es la menos
importante en las condiciones más frecuentes de operación.
Los HC sin quemar que normalmente aparecen en los gases de escape [6] se ven en la Tabla 12.2
Tabla 12.2
Ejemplo de distribución de diferentes HC en los gases de escape [6]
Carbono, porcentaje de HC totales
Parafinas Olefinas Acetilénicos Aromáticos
33 27 8 32
De todos ellos están considerados como mayores contaminantes el benzol y los aromáticos
polinucleares (PAH) por ser cancerígenos.
El efecto de enfriamiento de pared [10], está regido por el balance térmico de la mezcla en las
proximidades de la superficie de la cámara de combustión. En esta zona, si las pérdidas de calor por
233
Capítulo 12. Emisiones Contaminantes
conducción y radiación en la pared son mayores que el calor proveniente de los gases quemados
adyacentes, los radicales activos difundidos en la mezcla no alcanzan la temperatura suficiente para
continuar la combustión. Las variables que más influyen en el espesor del wall quenching son la
temperatura, la presión, y el dosado de la mezcla, así como el material y la temperatura de la pared de
la cámara de combustión. Cuanto más alta sea la presión y la temperatura de la mezcla y su dosado
más próximo al estequiométrico, menor será este espesor.
Las emisiones de HC en el escape son menores que las producidas en los cilindros cuando la
temperatura de escape supera los 600 C y existe oxígeno disponible [2]. En estas condiciones parte
de los HC sin quemar se oxidan en el escape produciendo CO.
Las emisiones de HC sin quemar junto con los óxidos de nitrógeno son los responsables de la
formación del "smog fotoquímico".
El óxido nítrico (NO) y el dióxido de nitrógeno (NO2), usualmente están agrupados como NOx,
siendo el NO bastante predominante. La fuente principal de NO es el nitrógeno molecular del aire
usado como comburente en el motor. Los combustibles diesel y gasolina en sí mismos contienen
muy poco nitrógeno, esto hace que su contribución a la formación de NO sea poco significativa [3].
Las emisiones de NOx son máximas cuando el motor funciona a par máximo, es decir, cuando la
presión media efectiva (pme) del motor es máxima. Así pues, las formación de NOx está muy
afectada por las presiones y temperaturas en las carreras de expansión y escape, los óxidos de
nitrógeno se disociasen en oxígeno y nitrógeno desplazándose el equilibrio químico hacia estos
compuestos. Esto no ocurre en la práctica y la razón puede residir en la complejidad del mecanismo
de formación [14]. De lo anterior se puede deducir que los óxidos de nitrógeno aparecen en los
MEP, pero tienen mayor influencia en los MEC dado que este tipo de motor funciona a más altas
presiones.
234
Capítulo 12. Emisiones Contaminantes
Partículas 1
En los MEP las partículas tienen tres distintos orígenes [3]: el plomo de los combustibles,
adicionado para incrementar el número de octano y evitar así el fenómeno de las detonaciones
espontáneas (knocking): el azufre propio del combustible y el humo negro (hollín), dependiendo
directamente del control que se tenga en el contenido de plomo y de azufre en el combustible y del
buen ajuste del motor, se pueden controlar dichas emisiones.
A diferencia de los MEP, en este tipo de motor el combustible es inyectado dentro del cilindro justo
antes de que se inicie la combustión. La distribución heterogénea del combustible dentro de la
cámara de combustión durante el proceso de combustión da como resultado la distribución
heterogénea de temperaturas y de composición de los gases quemados. Para entender mejor el
proceso de formación de los contaminantes provenientes de la combustión de los MEC, conviene
decir que comprende dos fases: una fase en la que una fracción del combustible es quemada durante
las premezclas inmediatamente después del tiempo de retraso, seguido por una fase de llama de
difusión. En la fase de llama premezclada, la composición de la mezcla varía bastante de la
estequiométrica. En la fase de llama de difusión, la mezcla se acerca a la estequiométrica.
En los MEC, como en los MEP, la máxima temperatura alcanzada gobierna la formación de NO. La
fracción de mezcla quemada en la región de premezcla tiene una gran influencia, porque cuando ésta
es comprimida en la siguiente fase, aumenta considerablemente su presión y su temperatura,
incrementándose así, la formación de NO. Esos gases se expanden luego en la carrera de potencia y
se mezclan con aire o con gases más fríos, lo cual congela la concentración de NO formada. La
presencia de aire frío en la cámara de combustión es propia de los MEC y explica por qué la
disminución de NO ocurre más rápido en este tipo de motores que en los MEP.
1
"Partícula": En el documento ISO 8178 es definida como efluente del escape atrapado en un papel de filtro después
de dilución a una temperatura del filtro no mayor de 52 C [5, 18].
235
Capítulo 12. Emisiones Contaminantes
pueden causar emisiones de CO. Esto podría ocurrir a bajas cargas y a máximas cargas a altas
velocidades [3].
En los MEC, la emisión de HC sin quemar es de importancia secundaria, debido a que funcionan con
dosados relativos inferiores a la unidad. Sin embargo, se debe destacar que la mayoría de los
hidrocarburos pesados son absorbidos por las partículas carbonáceas (humo) en forma de Fracción
Orgánica Soluble (SOF) [12], como se indicará más adelante.
Partículas
Si durante el proceso de combustión existen zonas con dosados muy ricos, al alcanzarse elevadas
presiones y temperaturas con falta de oxígeno, puede suceder que la cadena del hidrocarburo
comience a fracturarse y deshidrogenizarse, pudiendo quedar fácilmente convertido en carbón. Este
fenómeno suele darse, bien en la zona interior del chorro de combustible que proviene del inyector,
donde el dosado es muy rico y, por tanto, hay falta de oxígeno, o bien en el combustible depositado
en las paredes. Afortunadamente, la mayor parte de estas partículas se oxidan posteriormente para
formar dióxidos de carbono, con lo que no aparecen en el escape [14]. Estas partículas generadas
durante la combustión, adsorben varias especies orgánicas (SOF) que contienen HC, derivados
oxigenados (ketonas, esteres, aldehídos, latonas, éteres, ácidos orgánicos) e hidrocarburos
aromáticos policíclicos acompañados por sus derivados de nitrógeno y oxígeno. Adicionado a éstos,
existe un poco de derivados inorgánicos (SO2, NO2, sulfatos) [3].
Además de las tres variables consideradas como fundamentales en el control de las emisiones de
partículas de los MEC que son: 1) el combustible, 2) el proceso de combustión, y 3) procesos de
postratamiento, se debe incluir una cuarta variable: el consumo de aceite lubricante [19]. Aunque
quizá el parámetros más determinante, por encima de las mejoras tecnológicas antes mencionadas, es
la pauta de conducción.
Los diseñadores de los MEC, buscan cada vez elevar más la presión media efectiva (pme) del
motor, con lo que se producen mayores distribuciones de temperatura en la zona del anillo del pistón,
lo que puede propiciar el desplazamiento del aceite lubricante hacia la cámara de combustión con las
indeseables consecuencias que trae con sigo, tales como disminución de la potencia efectiva, altas
tasas de consumo de aceite lubricante, mayor formación de partículas, mayor número de etapas para
reparación y revisión del motor y posibles daños en la cámara de combustión. Una sugerencia
acertada es disminuir la tasa de consumo de aceite lubricante hasta el punto óptimo, es decir, hasta el
punto en que no perjudique la vida del motor ni reduzca una de las principales ventajas de un MEC
que es la de su largo período de duración sin abrir el motor [19].
236
Capítulo 12. Emisiones Contaminantes
Esencialmente podemos decir que las emisiones de humo crecen con la cantidad de combustible
inyectada después del tiempo de retraso, ó dicho en otras palabras, dada una cantidad de combustible
inyectado fija, crecen con todos aquellos parámetros que tiendan a disminuir el tiempo de retraso,
como pueden ser: aumento de la relación de compresión, aumento del número de cetano (ver
Capítulo 9) del combustible, disminución del avance a inyección, etc.
A manera de resumen, en la Tabla 12.3, se presenta la importancia relativa de cada uno de los
diversos contaminantes emitidos por los MEP y los MEC.
Tabla 12.3
Importancia relativa de los distintos contaminantes con cada tipo de motor [14]
Tipo de motor
Contaminante
MEP MEC MEC Sobrealimentado
NOx Importante Muy importante Muy importante
CO Muy importante Poco importante Poco importante
HC Muy importante Poco importante Poco importante
ALDEHÍDOS Poco importante Importante Importante
Plomo Importante Nula Nula
Óxidos de azufre Poco importante Importante Importante
HUMOS Poco importante* Importante Muy importante
* Excepto para dosados anormalmente ricos.
Establecer una comparación entre diferentes contaminantes de vehículos que utilizan diferentes
combustibles puede llevar a falsas conclusiones si no se define claramente la base de comparación.
El gas natural ha sido objeto de creciente interés como una alternativa de bajas emisiones para
motores de automoción convencional, especialmente para autobuses de servicio urbano y camiones.
Las tendencias actuales consisten en convertir los antiguos motores MEP o MEC a Gas Natural
Comprimido (GNC) operando algunos cambios sobre diversos subsistemas del motor [4, 20], los
resultados obtenidos de esta conversión pueden ser adversos si no se hace una selección
adecuada de los nuevos aditamentos que requiere el GNC. Particularmente en la ciudad de
Houston, el Southwest Research Institute, realizó un estudio de conversión de motores a GNC,
obteniendo incrementos en los niveles de emisiones de NOx, CO, además el consumo específico de
combustible efectivo (gef) empeoró en un 25% [4]. Cuando se considere la conversión de vehículos
a otro combustible para el cual el motor ha sido diseñado, el equipo a comprar para la conversión
debe cumplir con los requerimientos de niveles de emisión legislados.
237
Capítulo 12. Emisiones Contaminantes
Si además de lo anterior, se logra instalar un buen sistema de control para mantener el nivel de
emisiones bajo, seguramente se podrá comprobar que el GNC es una alternativa al igual que el GLP
muy llamativa desde el punto de vista de disminución en las emisiones contaminantes. El sistema de
control debe cumplir con los siguientes requerimientos [13]: que permita mantener el dosado
estequiométrico cuando se opera con catalizador de tres vías, que permita controlar el sistema de
recirculación de gases de escape (EGR), control sobre el tiempo de ignición, detección y control de
las explosiones (knock), control sobre la potencia y control sobre la velocidad al ralentí y la
mariposa.
Los combustibles gaseosos, particularmente el GN y el GLP son excelentes combustibles en ciclos de
combustión de quemado pobre. La combustión completa puede obtenerse con mezclas aire -
combustible muy pobres. Con un factor λ = 1,4 o 1,5 (i.e., un exceso de aire del 40% al 50%) se
logran minimizar las emisiones de óxidos de nitrógeno, estando aún a rendimientos elevados del
motor [20]. Los fabricantes de motores están encaminando sus esfuerzos en el desarrollo de motores
de mezcla pobre.
Las emisiones de los automóviles están reguladas en todos los países industrializados. Desde
principios de la década de los 60 se ha venido incrementando la lista de contaminantes que son
regulados, y a su vez se ha ido cerrando el límite de emisiones que son permitidas. A principios de
los 70, los valores de concentración se convirtieron en valores de masa de contaminantes. Estos
valores fueron expresados como masas emitidas por unidad de distancia viajada, basada en un ciclo
de recorrido normalizado que corresponde a una secuencia de velocidades en función del tiempo,
simulando las condiciones de marcha del vehículo en un banco dinamométrico [3, 4, 12, 15, y 19].
Últimamente se empiezan a plantear alternativas de simulación en bancos de prueba dinámicos que
eliminan en gran número las faltas de repetitividad debidas a la participación directa del hombre en la
realización de los ensayos convencionales [15]. El método de muestreo a volumen constante (CVS)2
se utiliza para acular las emisiones correspondientes al ciclo de conducción normalizado. Los
mismos ciclos son usados también para medir el consumo específico de combustible (gef) de los
vehículos. Desafortunadamente, no se ha adoptado un procedimiento aplicable universalmente, por
lo que existen diversos tipos de ciclos de prueba, entre los cuales no existe ninguna correlación para
convertir de uno a otro. Los se resumen en la Tabla 2.4.
2
Para ampliar infromación sobre este método de recolección de muestras, ver SAE Recomended
Practice, (1978), "Constant volume sampler system for exhaust emission measurement", SAE
J1094a, 24p.
238
Capítulo 12. Emisiones Contaminantes
LA MUESTRA
En cualquier análisis, la mezcla debe ser uniforme y representativa del producto a analizar. Los
valores de los contaminantes recogidos en los ciclos de prueba mencionados anteriormente, deben
coincidir con el valor real de emisiones en el tubo de escape del vehículo.
La toma de muestras en ensayo de motores a través de la técnica CVS, se ha normalizado desde
1982 [3]. En este método, los gases de escape son diluidos con aire del ambiente filtrado para
mantener flujo total constante (gases de escape + aire) bajo todas las condiciones de prueba del
ciclo, usando flujómetros de gas seco se puede medir el flujo másico de aire de dilución, por lo tanto,
restando de la cantidad anterior, se puede conocer el volumen real de gases de escape [5]. Con esto
se pretende simular las condiciones reales de dilución de los gases de escape en el aire de la
atmósfera. El flujo másico total puede ser determinado de dos maneras:
239
Capítulo 12. Emisiones Contaminantes
motor, la ISO 3 definió el "Sistema de muestreo total el cual extrae una fracción extremadamente
pequeña del escape el motor y luego la diluye y la filtra, esto representa una aproximación adecuada
[5].
Tabla 12.4
Parámetros comparativos de los ciclos de conducción [3]
Japan Japan
FTP(b) FTP
Parámetro Unid. ECE(a) Calif. 10 12
72 75
modes modes
Velocidad media km/h 18.7 35.6 31.5 34.1 17.7 30.6
Velocidad media
km/h 27.1 41.7 38.3 41.6 24.1 39.1
(sin contar el ralentí)
Aceleración media m/s2 0.75 0.68 0.60 0.67 0.54 0.64
Deceleración media m/s2 0.75 0.66 0.70 0.71 0.65 0.60
Tiempo medio (de 1 secuencia
s 45 117 66 70 50 94
desde el arranque)
Número promedio (de cambios de
acel./decl. y vice-versa en una 1 3 5.6 6 2 5
secuencia)
Ralentí %Tiempo 30.8 14.6 17.8 18.0 26.7 21.7
Aceleraciones %Tiempo 18.5 31.4 33.5 33.1 24.4 34.2
Velocidad constante %Tiempo 32.3 21.9 20.1 20.4 23.7 13.3
Deceleraciones %Tiempo 18.5 32.1 28.6 28.5 25.2 308
(a) Economic Commission for Europe (b) Federal Test Procedure
Precauciones
A parte de las posibles interferencias analíticas que pueden aparecer durante el muestreo, se deben
tener una serie de precauciones para obtener resultados de las medidas reproducibles y confiables [5,
3
ISO indica International Standars Organization
240
Capítulo 12. Emisiones Contaminantes
• Estado del vehículo a ser ensayado, lo cual indica preacondicionar el vehículo antes de probarlo
• Estado térmico del motor claramente definido antes de iniciar el procedimiento normalizado
• Purgar el sistema de escape y las líneas de atrapado de partículas
• Revisar las posibles absorciones o adsorciones de ciertos productos en las paredes de las tuberías
conectoras, las cuales pueden causar deposiciones que alteran fuertemente las mediciones
• Revisar que los tubos de transferencia plásticos no sean propicios para que haya difusión de los
gases a ser analizados
• Para escapes con diámetros hidráulicos substancialmente grandes y que se sepa que van a tener
partículas, se recomienda hacer muestreo isocinético
• Es importante asegurar que el número de Reynolds en el túnel de dilución sea lo suficientemente
alto para garantizar flujo turbulento y que la distribución de las especies a ser medidas sea
uniforme a través de toda la sección del túnel [5]
• Cuando se va a analizar partículas, a menudo ocurre que filtros tomados del mismo "baño" tienen
características de filtración diferentes, lo cual puede alterar significativamente los resultados, por
lo que se recomienda cortar los filtros de la misma pieza de material [12]
• El aire de dilución debe ser analizado bajo las mismas condiciones en cada secuencia de medida
• Los instrumentos de medidas deben también ser considerados ya que los resultados obtenidos
pueden variar de acuerdo al tipo de instrumento, fabricante y sus grados de calibración y
mantenimiento
• Para el Detector de Ionización de Llama FID, es importante especificar el material del detector,
sus parámetros de operación (flujo másico de gas, presión, temperatura, voltaje) y su geometría
interna[3]
Posibles interferencias
En las diferentes técnicas analíticas empleadas, otros compuestos presentes en el escape pueden
generar una respuesta en el rango de medida correspondiente al compuesto que está siendo analizado.
Para NO analizado por Quimioluminiscencia (QL), otros compuestos de nitrógeno además del
NO2 (NH3, HCN, aminas, etc.) Pueden ser convertidos a NO, y moléculas inertes (O2, N2, CO2, y
H2O, etc.) pueden descargar por colisión con el NO2 especies excitadas sin emisión de luz para así
producir falsas lecturas. El agua muestra la interferencia más seria [3]. Su remoción por
condensación elimina su efecto de descarga y mantiene valores correctos para el NO, el cual es
insoluble en agua. Disminuyendo la presión, aumentando la temperatura por encima del punto de
241
Capítulo 12. Emisiones Contaminantes
rocío, o removiendo el agua por difusión a través de una membrana semipermeable se supera este
porcentaje. Sin embargo, un incremento excesivo de la temperatura en la líneas de la muestra puede
causar una reducción de NO2 a NO por los hidrocarburos que están presentes, los cuales son
catalizados por los metales usados en las líneas, alternando así la relación NO / NO2.
TÉCNICAS DE NUESTREO
• Métodos de absorción: en los que el gas se hace burbujear a través de una solución absorbente
que retiene el contaminante por disolución o por reacción química, dando lugar a la formación de
compuestos que posteriormente pueden ser determinados por diversos métodos.
• Trampas criogénicas: en los que la muestra, generalmente muy volátil, se pasa a través de unas
columnas capilares inmersas en argón líquido, para que al disminuir la temperatura se condensen
y permanezcan inalterados para más tarde ser analizados en el laboratorio. Es necesario eliminar
previamente las sustancias o partículas como el agua que, al congelarse, obturaría los conductos
de aspiración.
• Bolsas de Teflón 4 , Tedlar 5 y Mylar 6 : Las cuales tienen capacidades desde 10 hasta 100 litros de
gas, pero son delicadas, difíciles de limpiar, y fácilmente se pueden perder materiales por
permeabilidad. Tedlar usualmente es la más recomendada [3].
• Frascos de vidrio: Los cuales son de fácil limpieza y de confiabilidad en el almacenamiento de
la muestra, pero son frágiles y tienen volúmenes muy limitados [16]
• Recipientes metálicos: los cuales son usados una vez se les ha efectuado un vacío. Son fáciles
de limpiar y se puede utilizar la misma muestra en ellos contenidos por varios análisis [16]
4 Politetrafluoroetileno
5 Pilivinil fluoride
6 Politereftalano
242
Capítulo 12. Emisiones Contaminantes
Métodos espectroscópicos
La absorción espectroscópica está basada en la medición de los efectos creados por la interacción de
la radiación electromagnética con los estados de energía interna de las moléculas [18]. Una molécula
con un momento dipolar eléctrico o magnético, absorbe energía (quanta) específica en la radiación
incidente, dependiendo de la longitud de onda, mientras ésta sufre transiciones de energía desde el
estado básico hasta los niveles más altos de energía permitidos. Esta excitación de los niveles más
altos de energía puede afectar los electrones de la molécula, las vibraciones interatómicas, o las
rotaciones de la molécula alrededor de los diferentes ejes. La determinación espectroscópica de
contaminantes atmosféricos se basa en las propiedades físicas de las moléculas y no sobre su
actividad química [17].
Cada gas muestra longitudes de onda características por las cuales absorben radiación
electromagnética que excita sus moléculas. Esto se refleja por un espectro de absorción del gas.
Para la mayor parte de los contaminantes atmosféricos hay bandas (longitudes de onda) específicas
de absorción en el infrarrojo, el visible o el UV cercano, basadas en la ley de Beer – Lambert [18].
Esta ley expresa la relación entre la intensidad de la radiación I a través de un recipiente de longitud
L transparente a la radiación y que contiene una sustancia absorbente en concentración c:
I = e − kcL
Donde k es el coeficiente de absorción de la sustancia a ser analizada.
Para tener una idea del grado de precisión de estos instrumentos en la determinación de
contaminantes, veamos la Tabla 12.5 [18]:
Absorción en el infrarrojo
La región infrarroja abarca las regiones del espectro comprendidas entre los números de onda 7 de
12800 a 10 cm-1 aproximadamente, lo que corresponde a las longitudes de onda 0.78 a 1000 μm .
Tanto desde el punto de vista de las aplicaciones como de los instrumentos es conveniente subdividir
la región infrarroja del espectro en tres porciones denominadas infrarrojo cercano, medio y lejana
[18].
7 Número de onda ( σ ) se define como el número de ondas por centímetro, y es otra forma de describir la radiación
electromagnética. Cuando la longitud de onda en el vacío se expresa en centímetros, el número de onda es igual a
1 / λ [18].
243
Capítulo 12. Emisiones Contaminantes
Tabla 12.5
Grado de precisión de los equipos espectroscópicos en la determinación de contaminantes
Concentración Concentración Error relativo
Contaminante
conocida en ppm medida en ppm %
CO 50 49.1 1.8
Metiletilcetona 100 98.3 1.7
Alcohol metílico 100 99.0 1.0
Óxido de etileno 50 49.9 0.2
Cloroformo 100 99.5 0.5
Las moléculas formadas por al menos dos átomos diferentes (CO, NO, CO2, NO2, H2O, etc.)
excepto los gases que están constituidos por átomos idénticos (H2, O2, Cl2, N2, etc.), absorben
radiación infrarroja convirtiendo la energía lumínica recibida en energía vibracional - rotacional de
sus moléculas, lo cual puede ser detectado en forma de calor.
Para análisis de gases de combustión, especialmente de los MCIA, se utilizan analizadores sin
dispersión espectral, los cuales son llamados infrarrojos no dispersivos NDIR. En lugar de utilizar un
monocromador para separar la luz, los instrumentos NDIR utilizan la absorción total sobre la
longitud de onda dada [6].
En los NDIR de filtro positivo, un rayo pasa a través de una celda de absorción que contiene el
gas a ser analizado y el otro cruza una celda que contiene un gas de referencia no – absorbente (como
el nitrógeno) y los dos rayos pasan luego a una celda de compensación y finalmente al detector,
provocando un desplazamiento del diafragma de sensor capacitivo como se observa en la Figura 12.3
[3]. El NDIR de filtro negativo no tiene celda de referencia, por lo que una sola celda contiene el gas
a ser analizado, el cual es atravesado por ambos rayos. En la salida, en el primer rayo pasa a través
de una celda de compensación y el segundo rayo pasa a través de una celda llena con el gas a ser
analizado y luego ambos llegan al sensor diferencial.
Los analizadores NDIR se utilizan para medir muchos gases, incluyendo CO, CO2, CH4, NO,
SO2, e hidrocarburos. En la contaminación de los MCIA son utilizados principalmente para
determinar CO y CO2. El CO absorbe la radiación infrarroja de 4,6 μm y para su determinación se
utiliza el analizador NDIR del tipo de filtro positivo [16]. Especial cuidado se debe tener al emplear
este tipo de analizador en la determinación de HC, dado que no tiene la sensibilidad suficiente para
distinguir el tipo de hidrocarburo que está en la celda, por lo que los resultados se dan en términos de
244
Capítulo 12. Emisiones Contaminantes
hexano 8 Lla recomendación más adecuada es utilizar un detector de ionización de llama (FID) para
determinar estos últimos.
Figura 12.3. Ilustración esquemática de un sensor infrarrojo. 1 Cámara receptora con volúmenes de
compensación V1 y V2. 2. Sensor de flujo. 3 Celda de medición. 4. Disco ranurado con motor. 5.
Emisor infrarrojo [3]
NO + O3 ⇒ NO2 * +O2
NO2 * ⇒ NO2 + hv
8 Tomado de las notas de clase de la asignatura Métodos Analíticos para el Control de la Contaminación Ambiental.
245
Capítulo 12. Emisiones Contaminantes
donde la corriente de aire u Oxígeno es expuesta a una radiación ultravioleta proveniente de una
lámpara. La reacción es:
3O2 + hv → 2O3
Figura 12.4. Ilustración de un analizador de gases de combustión por ionización por llama.[3]
246
Capítulo 12. Emisiones Contaminantes
Figura 12.5. Ilustración del quemador de una analizador de gases de combustión por ionización por
llama [3]
Métodos electroquímicos
E = E 0 + 2,3(RT / nF )LogA
247
Capítulo 12. Emisiones Contaminantes
vías en la técnica potenciométrica. Esta es aplicada en el sensor de oxígeno λ usado para mantener
el dosado cercano al estequiométrico aguas arriba del catalizador, aprovechando la propiedad de
ciertos óxidos refractarios tales como el zirconio para dar una conductividad iónica a altas
temperaturas (> 300C) debido a la movilidad de los iones O--.
Métodos Cromatográficos
• El tiempo de elución, el cual es ajustado para un sistema y un compuesto dado: esto ayuda a
caracterizar el compuesto.
• El área de pico, la cual sirve para determinar el componente de la mezcla.
• El ancho del pico, cuya estreches caracteriza la calidad de la separación del componente en la
mezcla.
Para analizar los contaminantes del aire usando substratos contenidos en la columna, se utilizan
dos métodos principalmente:
Cromatografía de gases (CG): en la que la muestra se volatiliza y se inyecta en la cabeza de un
columna cromatográfica. La elución se produce por el flujo de una fase móvil de un gas inerte, y a
diferencia de la mayoría de los tipos de cromatografía, la fase móvil no interacciona con las
moléculas del analito; su única función es transportar el analito a través de la columna. Actualmente
se está incrementando el uso de columnas capilares de diámetros muy pequeños de 0.1 a 0.5 mm y
248
Capítulo 12. Emisiones Contaminantes
muy largas, hasta varias docenas de metros de longitud 9 las cuales son mucho más efectivas,
presentan mayor permeabilidad y mejor selectividad y permiten análisis mucho más rápidos.
Cromatografía líquida de alta eficacia: Se lleva a cabo a bajas temperaturas pero requiere altas
presiones (hasta 50-60 MPa) para transportar la mezcla a ser separada y los líquidos, permitiendo
una desorción selectiva de los componentes adsorbidos. El solvente de desorción puede mantener
una composición constante o tener un gradiente de composición como una función del tiempo para
mejorar las separaciones.
Uno de los componentes esenciales del GC y el HPLC es el detector, el cual se encarga de indicar y
medir la cantidad de analito saliente de la columna cromatográfica. Existen diversos tipos de
detectores de los cuales los más usados son los siguientes:
Detector de ionización de llama (FID) antes visto, es un detector sensible a la masa, más que un
sistema sensible a la concentración. En consecuencia, este detector tiene la ventaja de que los
cambios en el caudal de la fase móvil tienen poco efecto sobre la respuesta del detector [11].
Catarómetro, (o detector de la conductividad térmica TCD), el cual hace una comparación
continua del componente de mezcla. El flujo de calor es generado por un resistor eléctrico [3]. Este
detector tiene poca sensibilidad.
Detector de captura electrónica (ECD), el cual es altamente selectivo siendo muy sensible a las
moléculas que contienen grupos funcionales electronegativos tales como halógenos, peróxidos y
grupos nitro, en cambio no es sensible a grupos funcionales como aminas, alcoholes e hidrocarburos.
Un electrón del emisor provoca la ionización del gas portador y la producción de una ráfaga de
electrones, resultando una corriente constante entre un par de electrodos. Sin embargo, la corriente
disminuye en presencia de moléculas orgánicas que tiendan a capturar los electrones [11].
Detector termoiónico, también llamado NPD (nitrogen Phosphorus detectors), es altamente
selectivo. Este detector es altamente sensitivo para medir nitrógeno y derivados fosforosos.
Detector fotométrico de llama (FPD), es altamente selectivo y filtra la luz emitida por la
combustión de los componentes de una llama de hidrógeno y mide la intensidad de la luz emitida.
Este detector es altamente sensitivo para uso en medidas de compuestos de azufre [11].
Detector de fotoionización (PID), en este detector, el eluyente de la columna se irradia con un
haz intenso de radiación ultravioleta que provoca la ionización de las moléculas al aplicar un
potencial a través de una celda que contiene los iones producidos, se origina una corriente de iones
que es amplificada y registrada [11].
9 Tomado de las notas del curso Métodos Analíticos en el Control de Contaminación Ambiental
249
Capítulo 12. Emisiones Contaminantes
Compuestos de Carbono
La técnica NDIR descrita antes, se utiliza prácticamente en todas partes para determinar CO2 y CO.
La precisión de la medición es esencialmente gobernada por una calibración adecuada. Debido a que
este método es sensible a la presión, este parámetro deberá mantenerse fijo durante el análisis y la
calibración.
Históricamente el NO2 fue primero analizado por el método Griess - Saltzmann descrito en la
referencia [16], en la cual se explican otras técnicas muy usuales como la del ácido fenol disulfónico.
También se han propuesto métodos electroquímicos. La espectrometría Ultravioleta (NDUV)
presenta buena aplicabilidad para determinar NO2, pero requiere remover previamente el vapor de
agua.
La técnica de la Quimioluminiscencia descrita antes es la más usada universalmente [6]. Se
encuentran instrumentos disponibles a operar bajo presión atmosférica y aún presiones reducidas.
Amoniaco y aminas
• Método Nessler, se basa en la formación del amónico con el Yoduro de Mercurio, por lo que
presenta problemas por la toxicidad del Hg. Se encuentra en la norma BOE 11-5-76.
• Método Alcalimétrico, el cual permite determinar sustancias en concentraciones altas, y el
• Método del Fenolato Sódico, en el cual el reactivo absorbente es fenol sódico mas
hipoclorito sódico, forman un compuesto de color azul que absorbe a 626 um.
El amoniaco puede ser medido también por Quimioluminiscencia en forma de NO, después de una
oxidación térmica. Otro método recolecta el amoniaco en un filtro y desarrolla anillos migratorios
que son detectados por una solución de o - ftaldicarboxilaldeido. El amoniaco puede también ser
medido por cromatografía iónica. Este es el método recomendado por la EPA después de atraparla en
ácido sulfúrico diluido.
Las aminas han sido atrapadas en los gases de escape por absorción en una solución ácida o en
una lana de vidrio impregnada con ácido oxálico. Las soluciones obtenidas son analizadas
espectrónicamente por ideometría o cromatrográficamente por GC con un detector termoiónico
250
Capítulo 12. Emisiones Contaminantes
calibrado en nitrógeno. Sin embargo, estos métodos no so lo suficientemente sensitivos para detectar
los contenidos tan bajos de aminas presentes en el escape. En el método de la EPA el atrapado se
realiza en ácido sulfúrico diluido. La solución obtenida es metida por GC con detección especifica de
nitrógeno [3].
Compuestos orgánicos
Al ser estos de características tan variadas, se pueden utilizar varios métodos para su determinación.
En el análisis de contaminantes de los MEP, estos pueden ser recolectados en bolsas plásticas no
reactivas o en botellas de acero inoxidable cuando son hidrocarburos livianos no condensables y su
análisis debe hacerse rápidamente [16].
Generalmente se utiliza un GC con detector FID, filamento incandescente y captura electrónica.
Cada columna y detector deben ser calibrados para el gas a determinar. La concentración total de
hidrocarburos puede ser determinada con un único detector FID, la cual generalmente permite la
separación de hidrocarburos con más de seis carbones. Sin embargo, la elección de la trampa
absorbente es crucial. El carbón activo se considera de uso universal para lograr este efecto.
Los PAH son absorbidos por las partículas de humo de los MEC o están presentes en la fase gaseosa
de la fracción absorbida en las partículas, especialmente para PAH con cuatro anillos tales como el
fluoranteno y pireno, de los cuales el 90% permanece en las partículas diesel. La mayoría de los PAH
con tres anillos permanecen en la fase gaseosa para los MEP y el 50% de ellos en los MEC, y más de
un 99% de los PAH con cinco anillos o más son retenidos en las partículas [3]. Los PAH recolectados
en el túnel de dilución durante la filtración de las partículas se hallan en forma de SOF. El tipo de
filtro puede afectar los resultados analíticos de los PAH por lo que se recomiendan filtros fabricados
en teflón - vidrio como el “Pallflex” [12]. Para el análisis de los PAH presentes en las emisiones del
escape de los MCIA, se usa como solvente orgánico el diclorometano (DCM) porque permite además
el análisis de la mutagenicidad.
Análisis de partículas
Para determinar la contribución de los compuestos orgánicos a las partículas diesel, las técnicas
universalmente empleadas para la determinación de partículas son 1) la extracción para solventes que
mide la fracción orgánica soluble (SOF) y 2) el método del horno al vacío que mide la fracción de
volátiles totales (TVF) contenidos en las partículas. En la actualidad se han desarrollado otras
técnicas como la DFI - GC (análisis por GC de inyección directa de filtro) para medir la fracción de
orgánicos volátiles (VOF) absorbidos por las partículas [12].
251
Capítulo 12. Emisiones Contaminantes
Las cantidades VOF, SOF, TVF no son exactamente lo mismo, pero están estrechamente
relacionadas. Específicamente, la VOF es el material orgánico que responde a un detector FID y
tiene un punto de ebullición ligeramente inferior, pero muy cercano a los 600C. La SOF es cualquier
material (independiente del punto de ebullición) que es soluble en cloro - metileno en unas
condiciones de ensayo especificadas. La TVF incluye todo el material filtrado que puede ser
volatilizado a 225C y una presión absoluta de 5 Torr [12].
Conclusiones
Muchas preguntas quedan sin responder hasta el momento en el campo de las emisiones
contaminantes sin embargo, se han venido presentando alternativas para su disminución en la
automoción, entre ellas se pueden destacar dos tendencias con proyección para el futuro y sobre las
cuales todavía queda por investigar: 1) el vehículo eléctrico el cual para obtener la energía mecánica
necesaria para su movimiento no lo hace por medio del motor de combustión interna alternativo sino
a través de una batería y un generador; el punto débil de esta tecnología radica justo en la batería,
debido a que presenta en la actualidad problemas de autonomía del vehículo y tiempo de recarga muy
largo. En España han entrado al mercado los primeros coches eléctricos de fabricación con
tecnología propias a partir del primero de Abril de 1.996, representando competencia sólo al grupo de
los vehículos de transporte de cargas liviano ya que por las razones antes mencionadas, esta
tecnología aún no compite con los vehículos tipo turismos o deportivos. El uso del coche eléctrico se
limita a la ciudad, por lo que también se ha planteado la posibilidad de utilizar los coches híbridos,
que funcionan como vehículo eléctrico en la ciudad y como vehículo con uso de combustible para
realizar largos recorridos, su desventaja radica en el alto costo de inversión y 2) el uso de las
llamadas pilas de combustible o celdas de combustible se saltan toda la termodinámica, obteniéndose
energía mecánica sin liberación de calor, lo cual parece muy prometedor.
Las limitaciones cada vez más restringidas de las emisiones contaminantes producidas por los
MCIA, han llevado a los fabricantes de motores a desarrollar nuevas técnicas y diseños para mejorar
su producto. En general las ganancias más importantes se han logrado a Través de cambios en el
diseño de la cámara de combustión, mejoras en el sistema de combustión, utilización de dispositivos
de pos - tratamiento para los gases de escape, disminución de la tasa de consumo de aceite lubricante
y tratamientos tanto al aire como a los gases de escape que puedan ser inyectados en el colector de
admisión.
Establecer una comparación entre los diferentes contaminantes producidos por los dos tipos de
motores analizados en este documento requiere de mucho tacto, se puede concluir que cada
contaminante tiene mayor o menor importancia en función de cuatro parámetros principales que son:
252
Capítulo 12. Emisiones Contaminantes
El monitoreo de las emisiones de los motores diesel se concentra en determinar los niveles de
emisión de material particulado. Se emplean dos métodos estandarizados:
Método de filtrado
Una cantidad especificada de los gases de emisión se hacen pasar a través de un filtro. El nivel de
decoloración del filtro provee un índice de la cantidad de hollín contenido en la muestra.
Figura 12.6. Medidor de opacidad. 1. Filtro de papel, 2. Inyector de gas, 3. Fotómetro, 4. Transporte de
pape, 5. Medida de volumen, 6. Electroválvula para paso de aire fresco, 7. Bomba [3]
253
Capítulo 12. Emisiones Contaminantes
El equipo extrae una cantidad especificada de los productos de combustión y la hace pasar a
través de un filtro de papel en forma de banda. Para obtener resultados comparables en diferentes
pruebas se debe gravar el volumen de gas procesado en cada prueba. El equipo convierte la
información a una forma estandarizada. El sistema también monitorea y compensa por otros factores
como presión, temperatura y volumen muerto entre el tubo de muestreo y el filtro. Un fotómetro
reflectivo se usa para evaluar la opacidad del filtro de papel. El resultado es indicado como un
numero BSZ o como concentración másica (mg/m3).
Este método esta basado en la atenuación en la intensidad de una haz de luz el cual es transmitido a
través de la muestra. Durante la prueba, una bomba diverge una porción del los gases de
combustión hacia la cámara de pruebas. Este procedimiento previene que la presión y sus
fluctuaciones deterioren los resultados obtenidos.
Dentro de la cámara de prueba, un haz de luz es transmitido a través de los productos de
combustión. Un equipos fotoeléctrico monitorea la reducción en la intensidad del haz de luz, lo cual
es posteriormente indicado como un porcentaje de opacidad o como un coeficiente de absorción. La
longitud de la cámara de prueba debe ser perfectamente definida y la ventana de la cámara de prueba
debe ser equipada con protección térmica contra los depósitos de hollín para asegurar resultados
precisos y consistentes. Los resultados asumen la forma de un proceso continuo en el tiempo. El
equipo de medición automáticamente determina el valor máximo y calcula el promedio de varias
pruebas.
REFERENCIAS
[1] Beshouri, G.M. (1991), “On the Development of modern analysis techniques for single
cylinder testing of large - bore engines”, Transaction of the ASME, journal of engineering for
gas turbines and power, 113, 390 – 398
[2] Crouse, W.H., (1971), “Automotive Emission Control”, McGraw-Hill, New York
[3] Degobert, P. (1995), “Automobile and Pollution”, Technip, Paris
[4] Fritz, S.G., et al, (1992), “Emission from Heavy Duty trucks converted to CNG”,
Transaction of the ASME, journal of engineering for gas turbines and power, 114, 561 – 567
[5] Graze, R.R., (1994), “Development and test of a fractional sampling system for diesel engine
particulate measurement”, Transaction of the ASME, journal of engineering for gas turbines
power, 116, 763 – 765
[6] J.B. Heywood, (1988), “Internal Combustion Engine Fundamentals”, McGraw-Hill, New
York
254
Capítulo 12. Emisiones Contaminantes
[7] Hollemans, B., et al. (1993), “Emission Comparison of LPG / Gasoline / diesel in passenger
cars”, TNO road vehiculo research institute.
[8] Kahir, M.K., (1992), “Progress in diesel engine emission control” Transaction of the ASME,
journal of engineering for gas turbines and power, 114, 568 – 576
[9] Muñoz M., (1989), “Motores de Combustión Interna Alternativos”, Sección de publicaciones
de la ETSII, Madrid.
[10] Notas del autor en su formación doctoral, (1994), “seminario Procesos de combustión I”,
Madrid.
[11] Narros, G.A., (1995), “Métodos de separación”, 49 pp
[12] Newkirk, M.S., (1992), “Application of methods for determining total organic contribution
to diesel particulates”, Transaction of the ASME, journal of engineering for gas turbines
and power, 114, 578 – 589
[13] Nobel, A.D., et al., (1991), “Control System for a low emission natural gas engine for urban
vehicles”, SAE paper No. 910255
[14] Payri, F., et al., (1979), “Soluciones a la contaminación de motores térmicos alternativos de
automoción”, Monografía sobre medio ambiente, serie tecnología y recursos naturales, Vol
1.
[15] Quissek, F., Bandel, W., et al., (1992), “Transient development of SI engines on a dynamic
engine test bed”, AVL – List Gmbh, Graz – Austria.
[16] Sánchez, M.S., et al., (1986), “Determinación de gases”, Revista: medida de la
contaminación atmosférica.
[17] Sánchez, M.S., et al., (1986), “Técnicas instrumentales”, Revista: medida de la
contaminación atmosférica.
[18] Skoog, D.A. y West, D.M., (1984), “Análisis Instrumental”, Interamericana, México
[19] Springer, K.J., (1989), “Diesel lube oils – fourth dimension of diesel particulate control”,
Transaction of the ASME, journal of Engineering for gas turbines and power, 111, 355 –
360
[20] TNO Systems for Gas Engines, Road Vehicules Research institute department of internal
combustion engines, Catalog.
255