Modelado de la combustin en motores

Disel: revisin del estado del arte

Modelling of combustion in Diesel engines:
a review of the state of the art
Gabriel Fernando Garca Snchez; Jorge Luis Chacn Velasco*; Arlex Chaves Guerrero
Grupo de Investigacin en Energa y Medio Ambiente GIEMA, CICT-Energas,
Universidad Industrial de Santander (UIS),
Cra. 27 Calle 9, Bucaramanga, Colombia
* jchacon@uis.edu.co
Fecha Recepcin: 02 de noviembre de 2012
Fecha Aceptacin: 07 de junio de 2013

Resumen
El estudio del proceso de combustin en motores Disel ha sido, durante aos, un tema de gran inters
debido principalmente a la creciente necesidad de mejorar el rendimiento de los motores, reducir el
consumo de combustible y reducir la emisin de contaminantes. Los estudios experimentales presentan
la desventaja de requerir grandes recursos econmicos y tiempo en su realizacin, razn por la cual se
utilizan modelos computacionales que permiten estudiar el proceso de combustin de una forma ms
econmica y prctica. En este artculo se presenta una revisin de varios de los modelos de combustin
en motores Disel ms relevantes, desarrollados a nivel nacional e internacional, con el fin de dar al
lector una visin de la evolucin y el estado actual de este tipo de modelos.
Palabras clave: modelado de la combustin, motores diesel, modelos termodinmicos, modelos multidimensionales.

Abstract
The study of the combustion process in Diesel engines has been for years a subject of great interest
mainly due to the increasing need to improve engine performance, reduce fuel consumption and reduce
the pollutant emissions. Experimental studies have the disadvantage of require enormous efforts,
money and time, that is the reason why computer models are used to study the combustion process in
a more economical and practical way. This article presents a review of the most relevant Diesel engine
combustion models developed nationally and internationally, to give the reader an insight of the evolution
and current status of such models.
Keywords: combustion modeling, diesel engines, thermodynamics models, multidimensional models.

Introduccin

avances en lo que respecta a la reduccin de

emisiones desde los aos 80 [2], ellos siguen
siendo los responsables de la mayor parte de la
contaminacin en los principales centros urbanos,
incrementando as los problemas de salud de la
poblacin y el efecto invernadero [3]; por esta
razn actualmente hay grandes exigencias
con respecto al desarrollo de motores Disel
cada vez ms eficientes y limpios, lo cual se ve

Los Motores de Combustin Interna (MCI),

particularmente los motores Disel, son
ampliamente utilizados como fuerza motriz
en aplicaciones industriales como: transporte,
agricultura y generacin de energa [1], lo que
los hace indispensables para la sociedad actual.
Sin embargo, pese a que se han visto grandes
rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).
41

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

reflejado en las legislaciones mundiales sobre el

tema [1,4]. Con base en lo anterior, los estudios
computacionales sobre desempeo de los motores
Disel cobran gran importancia en nuestros das,
ya que permiten realizar, mejoras significativas
en el desempeo de los nuevos motores; tanto
es as que se puede afirmar que sin modelos de
simulacin que cuantificaran los efectos de los
procesos fundamentales en los sistemas del motor,
los avances logrados en los motores actuales
habran sido imposibles [4,5].
El proceso de combustin es el aspecto ms
importante de un motor de combustin interna [6],
por lo que su modelado es de gran relevancia al
simular el comportamiento global de un motor
Disel. Los modelos de la combustin se suelen
clasificar de acuerdo a la representacin espacial
del proceso, en este punto hay pequeas
diferencias en las categoras de clasificacin
utilizadas dependiendo del autor citado. De acuerdo
a Heywood [4] hay dos grandes grupos en los cuales
se pueden clasificar los modelos de combustin en
motores Disel: los modelos flui-dinmicos y los
modelos termodinmicos, a su vez estos ltimos
se pueden clasificar como cero-dimensionales,
fenomenolgicos o cuasi-dimensionales. Fenollosa
[2] propone una clasificacin de los mtodos de
clculo para la descripcin de la combustin en
dos grupos principales: mtodos de descripcin
detallada y mtodos de descripcin global.
Bracco [7], clasifica los mtodos en tres grupos:
cero-dimensionales,
cuasi-dimensionales
o
fenomenolgicos y multidimensionales. Esta ltima
forma de clasificacin ser la seguida por este
artculo, cabe destacar que las fronteras entre los
diferentes tipos de modelos en ocasiones son un
poco artificiales y algunos de ellos podran caber
dentro de varias categoras.
En los ltimos aos, debido al gran inters que ha
despertado el uso de biocombustibles, han surgido
varios modelos de combustin que buscan evaluar
el desempeo de motores Disel alimentados con
nuevos combustibles. Este es un tema de especial
inters para la regin nororiental colombiana,
la cual cuenta con una capacidad importante de
produccin de residuos de biomasa, como ha sido
descrito por Escalante et al. [8]. Razn por la cual
en este artculo, adems de presentar la evolucin
de los modelos de combustin Disel en el mundo,
se har nfasis en los modelos de combustin
recientes que buscan estudiar la factibilidad del
uso de biocombustibles.
Este artculo iniciar con la revisin de los modelos

multidimensionales, cero-dimensionales y cuasidimensionales de la combustin, comenzando

en cada caso con la presentacin de los trabajos
desarrollados en el mbito internacional y
finalizando con los trabajos desarrollados en
el mbito nacional; y finalizar presentando
una comparacin de los modelos desarrollados
recientemente con el fin de estudiar el desempeo
de biocombustibles en motores Disel.
Modelos multidimensionales
Los modelos multidimensionales o CFD
(Computational fluid dynamics) son modelos que
tienen el potencial de describir individualmente
cada uno de los fenmenos fsicos y qumicos
que tienen lugar dentro del proceso inyeccincombustin. Estos se basan en la resolucin
numrica de las ecuaciones diferenciales que
gobiernan los procesos que ocurren dentro de la
cmara de combustin [2].
El desarrollo de este tipo de modelos se remonta
a finales de los aos 70 [9-11], ellos suelen
utilizar dos marcos de referencia principales:
la descripcin euleriana continua [12,13] y la
descripcin lagrangiana discreta [14-16]; este
ltimo suele ser el ms usado debido a que
tiene ciertas ventajas desde el punto de vista de
solucin numrica [2]. A partir de la metodologa
de clculo multidimensional han nacido varios
cdigos CFD comerciales ampliamente utilizados,
como SPEED [17,18], FLUENT [19] y KIVA
[20,21], siendo este ltimo el ms usado a nivel
mundial [22]; y recientemente, se ha desarrollado
el software de cdigo abierto OpenFOAM
[23,24], el cual ha despertado gran inters en la
comunidad cientfica [22]. Recientemente los
modelos multidimensionales han sido utilizados
principalmente para estudiar los efectos de
la inyeccin y del flujo de combustible en el
desempeo y generacin de emisiones del motor,
como ejemplo de esto se pueden citar los trabajos
de Zheng et al. [25] en 2005, quienes utilizaron un
cdigo CFD para estudiar los efectos del ngulo de
inyeccin, la temperatura inicial y la composicin
del combustible en el ngulo de ingnicin, la presin
en cmara y las emisiones de NOx; Jayashankara
y Ganesan [26], quienes en 2010 presentaron
un modelo para el estudio de los efectos de la
presin de entrada y el ngulo de inyeccin sobre
el desempeo del motor y sus emisiones de NOx
y material particulado; y Shi y Reitz [27] que en el
mismo ao utilizaron un modelo para estudiar las
estrategias de inyeccin ptimas en un motor de

42

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

encendido por compresin, alimentado con Disel,

gasolina y mezcla E10, y operando a condiciones
de mediana y alta carga.
En Colombia se han presentado varios estudios
con modelos de este tipo en la Universidad de
Antioquia. En 2008, Gutirrez [28] estudi la
incidencia de la naturaleza del combustible y las
condiciones locales de presin y temperatura
en el proceso de inyeccin Disel, as como
las variables que ms afectan dicho proceso,
por medio de un modelo numrico CFD de la
atomizacin y evaporacin del chorro Disel. En
2009, Agudelo y Agudelo [29] desarrollaron un
modelo para simular los principales subprocesos
que ocurren en un chorro Disel usando un cdigo
CFD de libre acceso, dicho modelo se us para
evaluar el efecto del tipo de combustible, la presin
de inyeccin y la presin del gas ambiente en la
penetracin de la punta del chorro, el dimetro
medio de Sauter (SMD) y la masa de combustible
evaporada. En 2010, Delgado [30] present un
estudio comparativo del proceso de combustin
en trminos de los parmetros caractersticos que
lo describen, utilizando para ello el cdigo CFD
OpenFOAM.
Este tipo de modelos presenta como ventaja su
capacidad de describir cada uno de los fenmenos
que tienen lugar en el proceso inyeccincombustin Disel, teniendo en cuenta tanto la
geometra instantnea del cilindro como la variacin
espacial del campo de flujo. No obstante, tienen
la desventaja de un mayor consumo de tiempo
computacional en comparacin con los modelos
termodinmicos [22,31-33] y su sensibilidad a los
submodelos utilizados en los clculos, algunos de
los cuales an no estn bien definidos [2,31].

Disel de inyeccin directa, en donde se presenta

una fuerte liberacin de calor en la combustin
premezclada, se descubri que una sola ecuacin
de Wiebe no era adecuada para modelar la
liberacin de calor caracterstica del proceso, por
lo que se empez a utilizar una combinacin de
dos funciones de Wiebe (o funciones similares a
ella), esto se puede apreciar en los trabajos de
Dyechenko et al. [37,38], Ghoel [39,40], Watson
et al. [41], Miyamoto et al. [42] y Witt et al. [43].
En motores Disel modernos, en los que se han
incorporado nuevas tecnologas y avanzados
sistemas de combustin para disminuir las
emisiones y el consumo de combustible, una o dos
ecuaciones de Wiebe ya no resultan apropiadas
para simular el calor liberado en la combustin, por
lo que se han aplicado tres, cuatro o inclusive ms
como se puede observar en diferentes trabajos
[44-46]. Una descripcin detallada del proceso
de desarrollo y las aplicaciones de las diferentes
formas de la ecuacin de Wiebe se presenta en el
trabajo de Ghojel [36].
Otros modelos del tipo cero-dimensional describen
el problema teniendo en cuenta la fsico-qumica
de los fenmenos ocurridos durante la combustin
de forma simplificada, ellos en ocasiones
tambin son clasificados dentro del grupo de los
fenomenolgicos [2]; entre los modelos de este
ltimo tipo se encuentran los trabajos realizados
por: Way [47], en el que se presentan dos mtodos
para el clculo de la composicin y propiedades
termodinmicas del gas dentro de la cmara
durante la combustin; Hardenberg y Hase [48],
en el cual se obtuvo una ecuacin emprica
para calcular el tiempo de retardo en la ignicin,
permitiendo evaluar los efectos del combustible,
los parmetros del motor y las condiciones del
motor sobre dicho tiempo de retardo; Gardner y
Henein [49], en el que se desarroll un modelo
matemtico con el que se analizaron los efectos
de la relacin de compresin en los parmetros de
desempeo del motor; Harris y Pearce [50], donde
se realiz un modelo que expresa el par, la potencia
y el consumo de combustible como funcin de la
velocidad del motor a cero par y la disminucin de
la velocidad del motor por debajo de la velocidad a
cero par; y ms recientemente el trabajo de Arrgle
et al. [31], el cual identifica y cuantifica la influencia
de los parmetros de inyeccin y las condiciones
de operacin en la combustin Disel, para
esto se define un nuevo parmetro denominado
Tiempo Aparente de Combustin (ACT, de sus

Modelos cero-dimensionales
Estos son modelos basados principalmente en la
ley de conservacin de la energa, en los que se
asume homogeneidad espacial de las variables
implicadas en el proceso [34].
Muchos modelos cero-dimensionales se basan en
la descripcin matemtica de la forma de la ley de
liberacin de calor sin relacin directa con la fsica
que controla el proceso de inyeccin-combustin
[2]. El ms conocido de estos mtodos es la funcin
de Wiebe [35] que se ha utilizado desde los aos
70 y se ha vuelto tan comn que investigadores
en varias partes del mundo han dejado de citar la
fuente de la cual se origin [36]. Dicha funcin ha
sido adaptada a lo largo de los aos para su uso en
motores de diferentes caractersticas. En motores

43

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

siglas en ingls), que permite relacionar la ley

de liberacin de calor con la ley de inyeccin en
el proceso de inyeccin-combustin Disel; y el
trabajo de Chmela et al. [33], el cual presenta una
estrategia genrica para el desarrollo de modelos
de combustin cero-dimensionales usando las
ecuaciones bsicas para la tasa de reaccin de
Arrhenius y Magnussen.
Recientemente, se han realizado varios modelos
cero-dimensionales que buscan predecir el
desempeo de motores Disel alimentados
con diferentes tipos de biocombustibles. En
2006,Ramadhas et al. [1] desarrollaron un modelo
que analiza las caractersticas de rendimiento
de los motores Disel alimentados por biodisel
de semilla de caucho y sus mezclas. En 2010,
Gogoi y Baruah [51] desarrollaron un modelo que
predice el desempeo de un motor Disel mono
cilndrico de 4 tiempos alimentado con mezclas de
biodisel de aceite de Karanja; Bueno et al. [52]
obtuvieron un modelo por medio de un anlisis
exergtico con el fin de analizar el impacto que
ocasiona el uso de ster etlico de aceite de soya
en un motor Disel turbo-cargado, Colao et al.
[53] simularon el perfil de temperatura en un pistn
de un motor Disel bajo varias condiciones de
carga, Colao et al. [3] tambin presentaron un
modelo que simula la presin en cmara de un
motor Disel operando bajo varias condiciones y
usando diferentes mezclas de biodisel. En 2011,
Payri et al. [54] presentaron un modelo para el
estudio del desempeo del motor que cubre varias
de las falencias que otros modelos presentan
debido a las aproximaciones asumidas durante
su desarrollo, como la consideracin de la cmara
de combustin como un sistema cerrado, la no
consideracin de los cambios en la composicin
del gas al calcular los calores especficos de
ste y la consideracin de una ley de liberacin
de calor poco realista. Y en 2012, Patil y Akarte
[55] estudiaron el desempeo de un motor Disel
alimentado con metil ester de aceite de palma.
En la Tabla 1 se presenta una descripcin ms
detallada de estos ltimos modelos.
Estos modelos son relativamente sencillos teniendo
como ventaja un bajo tiempo computacional frente
a los modelos ms complicados, como los CFD
[22,31,33,56], sin embargo presentan la desventaja
de no poder describir en detalle los subprocesos que
ocurren durante el proceso global de combustin
Disel [33].
En Colombia se han desarrollado modelos cerodimensionales en la Universidad de Antioquia

y la Universidad Industrial de Santander. En

la Universidad de Antioquia, en el ao 2000,
Agudelo et al. [57] presentaron una metodologa
para la simulacin del proceso de combustin
en un motor de encendido por compresin
de aspiracin natural o turboalimentado;
ms adelante en el 2008 Agudelo et al. [58]
desarrollaron un modelo de diagnstico
exergtico de una zona y dos especies que se
utiliza para caracterizar la operacin de un motor
Disel. Por su parte en la Universidad Industrial
de Santander se desarrollaron trabajos tericos
por parte de Barbosa y Salcedo [59] y Rodrguez
[60], quienes desarrollaron herramientas de
simulacin de ciclo de un motor de cuatro
tiempos con fines acadmicos, y Chacn [61],
quien utiliz el mtodo de Levenberg-Marquardt
(LM) para ajustar los parmetros de la ecuacin
de Wiebe.
Modelos cuasi-dimensionales o fenomenolgicos
Los modelos cuasi-dimensionales son modelos
que presentan un planteamiento cero-dimensional
bsico, pero incluyen algn aspecto geomtrico
caracterstico del proceso (como la geometra
del chorro) o detalles espaciales adicionales
para alguno de los fenmenos involucrados en
la combustin (atomizacin, evaporacin, mezcla
de combustibles, etc.) [34]. Ellos suelen dividir la
cmara de combustin en varias zonas, haciendo
necesaria la resolucin de las ecuaciones de
conservacin de masa y energa para cada una
de stas.
Los modelos cuasi-dimensionales permiten
calcular la formacin/destruccin de emisiones y
tienen un bajo tiempo de clculo computacional
frente a los modelos multidimensionales. Como
inconveniente presentan la necesidad de ajustar
algunos de sus coeficientes con resultados
experimentales, debido al gran nmero de
simplificaciones que requieren [34].
Uno de los primeros trabajos bajo esta metodologa
fue el realizado por Probert [62], quien examin
la evaporacin del chorro de combustible usando
una funcin de distribucin de gotas; este trabajo
fue posteriormente ampliado por Tanaswa [63],
quien relacion explcitamente los procesos
de evaporacin y combustin por medio de un
coeficiente de combustin. Otros de los trabajos
pioneros de modelos fenomenolgicos son los de
Austen y Lyn [64] y Lyn [65] quienes plantearon
un modelo que relaciona cuantitativamente la

44

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

tasa de inyeccin de combustible con la tasa

de calor liberado; Cook [66,67] Nagao et al.
[68] y Shipinsky et al. [69,70] que describieron
el proceso de combustin de manera similar a
Austen y Lyn pero dividiendo la tasa de inyeccin
en zonas que se pueden evaporar, autoencender
y quemar de forma independiente.
Estos primeros trabajos estn basados en la
hiptesis de que la evaporacin es el fenmeno
que controla la combustin, posteriormente
aparecieron los trabajos de Whitehouse y Way
[71,72], en los cuales se desarroll un mtodo
para el clculo de la tasa de liberacin de calor
que utiliza ecuaciones simples basadas en un
modelo de una zona y que tiene en cuenta
los efectos de la tasa de mezcla entre el aire
y el combustible; Whitehouse y Sareen [73] en
el que se present un modelo de dos zonas
basado en el trabajo de Whitehouse y Way;
Grigg y Syed [74] en donde se consideraron
los efectos, sobre la liberacin de calor, de la
tasa de incorporacin del aire en el chorro de
combustible, la tasa de mezclado turbulento
de aire y combustible dentro del chorro, y la
cintica qumica de la combustin; y Khan et
al. [75] en el que se incluyeron los efectos de
choque de la pared y swirl. Otros investigadores
trabajaron en estrategias basadas en la teora
de chorro estacionario de Abramovich [76],
como Adler y Lyn [77] quienes la adaptaron al
estudio de chorros transitorios pero sin validarla
para condiciones de motor. Posteriormente, Rife
y Heywood [78] lograron adaptar el trabajo de
Alder y Lyn a las condiciones de motor.
Chiu et al. [79] usaron ecuaciones empricas
para describir la penetracin del chorro, su
ngulo de apertura y su trayectoria, con el fin de
describir la naturaleza transitoria de los chorros
Disel. Hiroyasu et al. [80] utilizan un enfoque
similar al de Chiu y otros para describir el chorro
Disel. Dent et al. [81,82] utilizan un tratamiento
del chorro similar al de Hiroyasu y otros, que
posteriormente es ampliado por el trabajo de
Kyriakides et al. [83]. Simultneamente a este
ltimo, Tinaut [84] trabaj en un modelo multizona en el que la masa inyectada se divide en
paquetes en los que se considera la atomizacin,
evaporacin y englobamiento del aire. Ms
adelante Lapuerta [85] y Desantes et al. [86]
incluyen efectos de interaccin del chorro con
el vrtice y clculos de la trayectoria de ste al
trabajo de Tinaut.
A partir de los 90 la mayora de los modelos del

tipo cuasi-dimensional han sido modelos multizona basados en el tratamiento mixto del chorro
inyectado [2]. Ejemplos de esta aproximacin se
pueden observar en los trabajos de Rakopoulos
y Hountalas [56], Bazari [87], Cui et al. [88],
Kouremenos et al. [89], Ottikkutti et al. [90],
Shenghua et al. [91], Torkzadeh et al. [92] y Bi
et al. [93]. Algunos de los modelos multi-zona
desarrollados en los ltimos aos son los de
Rakopoulos et al. [32] en 2007, cuyo fin es el
de estudiar el desempeo de un motor Disel
de inyeccin directa (ID) que usa aceite de
semilla de algodn y su Disel derivado como
combustible; Rakopoulos et al. [94] en 2008, el
cual es utilizado para examinar los mecanismos
de formacin de las emisiones en la combustin
al usar mezcla etanol-Disel; Komninos et
al. [95] en 2010, el cual se desarroll con el
fin de investigar la emisin de hidrocarburos
quemados y la formacin de CO en motores
HCCI alimentados con etanol puro e iso-octano
puro; y Komninos y Kosmadakis [96] en 2011,
cuyo objetivo es estimar el flujo de calor por las
paredes del cilindro de un motor HCCI.
En los ltimos aos tambin se han desarrollado
varios modelos de dos zonas, la mayora de ellos
con el fin de predecir el desempeo de motores
Disel alimentados con biocombustibles. Uno de
estos es el realizado por Ganapathy et al. [97]
en el 2009, que permite estudiar el desempeo
de un motor Disel mono-cilndrico de cuatro
tiempos, de inyeccin directa y enfriado por aire,
alimentado con biodisel de Jatropha. El mismo
ao Ganapathy et al. [98] utilizaron el enfoque
de optimizacin de Taguchi en combinacin
con un modelo de dos zonas para identificar
la configuracin ptima de parmetros de
entrada que maximizan el rendimiento trmico
de un motor Disel que usa biodisel Jatropha.
En 2010, Rajendra et al. [99] desarrollaron un
modelo de dos zonas con el objetivo de medir la
efectividad de la combinacin de dos mtodos
para mejorar el desempeo del motor y disminuir
las emisiones: el uso de biodisel de Jatropha
y el uso de materiales cermicos para aislar
la cmara de combustin (motor LHR). Gao
et al. [22] simularon la generacin de NOx y
material particulado en la combustin Disel
convencional, la combustin de alta dilucin
(HDC) y la combustin limpia de alta eficiencia
(HECC) en un motor Mercedes Benz 1.7 L,
usando un modelo de dos zonas. Por ltimo,
en 2011, Jagadish et al. [100], desarrollaron un

45

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

Modelos recientes para el estudio de la factibilidad del uso de biocombustibles

En la Tabla 1 se presenta una descripcin de
varios de los principales modelos que han
sido desarrollados en los ltimos aos con el
propsito de estudiar la factibilidad del uso de
biocombustibles en motores Disel. En dicha tabla
se presentan los objetivos y resultados de cada
trabajo, permitiendo as la comparacin de sus
respectivos alcances.
Como se puede apreciar, los modelos cerodimensionales
se
han
venido
utilizando
recientemente para comparar los parmetros de
desempeo de motores Disel alimentados con
biocombustibles, con los mismos parmetros
cuando dichos motores son alimentados con
Disel convencional, esto sin tener en cuenta
las emisiones contaminantes. Por otra parte los
modelos cuasi-dimensionales se han utilizado
para el mismo fin, pero incluyendo en los anlisis
las emisiones de contaminantes, esto debido a su
mayor poder predictivo a causa del uso de mltiples
zonas para los anlisis de la cmara de combustin.
Finalmente, los modelos multidimensionales se
han utilizado en casos en los cuales se necesita
comparar caractersticas espaciales especficas de
la inyeccin y del flujo de combustible, lo cual no
podra realizarse con modelos de otro tipo.

modelo que permite estudiar el desempeo de

un motor mono-cilndrico al ser alimentado con
varios tipos de biocombustibles. Una descripcin
ms detallada de estos modelos es presentada
en la Tabla 1.
En Colombia, en los ltimos aos, se han
desarrollado modelos de acuerdo a la metodologa
cuasi-dimensional en la Universidad Nacional y
la Universidad de Antioquia. En 2010, Ortiz [101]
elabor un modelo cuasi-dimensional con el fin de
analizar los efectos que tiene la modificacin de
la viscosidad del combustible sobre la liberacin
de energa, y Rodrguez [102] desarroll un
modelo de formacin de contaminantes de
combustin para la mezcla B5 que permite
determinar la incidencia de la estructura y
composicin del combustible en la generacin de
material particulado (MP), este ltimo se present
como una herramienta que permite obtener
resultados preliminares sobre el desempeo de
la mezcla antes de la implementacin en motores
de combustin interna. Mientras que en 2011
Villegas [103] utiliz un modelo de diagnstico
cuasi-dimensional para realizar un estudio
terico-experimental acerca de la incidencia de
las mezclas Disel-biodisel de aceite de palma
en la formacin de los xidos de nitrgeno y
material particulado.

Tabla 1. Comparacin de los modelos recientes para el estudio de biocombustibles.

Autores /
Pas / Ao

Zheng et
al. [25] /
China /
2005

Ramadhas
et al. [1] /
India / 2006

Tipo de Modelo

Estudio Realizado

Resultados

Multidimensional

Formacin de la mezcla de
gases, distribucin del NO en la
cmara, flujo de combustible en
el cilindro e incidencia del ngulo
de
inyeccin,
temperatura
inicial y composicin del gas
natural en el ngulo de ignicin,
presin en cmara y emisiones
contaminantes.

Retardar en el ngulo de inyeccin reduce

la emisin de NO. Incrementar la cantidad
de etano en el gas natural puede adelantar
la ignicin y por lo tanto aumentar las
emisiones de NO. Temperaturas iniciales
muy altas o muy bajas no son beneficiosas
para la salida de potencia y la emisin de
contaminantes.

Cero-dimensional

Efecto de la variacin de la
relacin aire combustible y la
relacin de compresin en la
eficiencia trmica al freno, la
presin pico y la temperatura
pico en el cilindro de un motor
alimentado con biodisel de
semilla de caucho, Disel
convencional y mezcla de Disel
con 20%vol de biodisel.

Se presenta un incremento en la
temperatura pico, la presin pico y la
eficiencia trmica al freno al aumentar la
relacin de compresin, y un decremento
en los mismos parmetros al aumentar la
relacin aire-combustible, tendencia que
fue seguida por todos los combustibles
analizados. Se observa la mayor presin
pico y eficiencia trmica al freno al
utilizarse Disel convencional, y la mayor
temperatura pico al utilizarse biodisel
puro.

46

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

Autores /
Pas / Ao

Rakopoulos
et al.
[32], [94]
/ Grecia /
2007

Agudelo
et al. [29] /
Colombia /
2009

Ganapathy
et al. [97] /
India / 2009

Gogoi y
Baruah [51]
/ India /
2010

Tipo de Modelo

Estudio Realizado

Resultados

Distribucin del dosado relativo,

temperatura, concentracin de
oxido ntrico y holln en el chorro
de combustible de un motor
Disel alimentado con aceite de
semilla de algodn, su biodisel
derivado, mezclas etanol-Disel
y Disel convencional.

Las zonas con dosados altos son ms

pequeas en los chorros de mezclas
etanol-Disel, aceite de semilla de
algodn y su biodisel derivado que en los
de Disel convencional, lo cual mejora la
combustin y mejora las emisiones.

Efecto del tipo de combustible,

la presin de inyeccin y la
presin del gas ambiente en
la penetracin de la punta
del chorro de combustible,
el dimetro medio de Sauter
(SMD) y la masa de combustible
evaporada.

Las propiedades del fluido afectan

significativamente los procesos de
atomizacin y vaporizacin y en menor
medida la penetracin del chorro.
Independientemente de las presiones de
inyeccin y del gas ambiente, el SMD se
incrementa con la viscosidad y la tensin
superficial mientras la tasa de evaporacin
se incrementa con la volatilidad del
combustible. Bajas presiones del gas
ambiente favorecen el proceso de
vaporizacin y la penetracin del chorro.
A medida que la presin de inyeccin se
incrementa el SMD disminuye y la tasa de
evaporacin aumenta.

Cuasidimensional

Incidencia
del
ngulo
de
inyeccin,
la
relacin
de
compresin
y
la
relacin
relativa aire/combustible en la
eficiencia trmica al freno, la
eficiencia mecnica, la eficiencia
volumtrica, la temperatura pico y
la presin media efectiva al freno,
de un motor Disel alimentado
con Disel convencional y
biodisel de Jatropha.

Para los dos combustibles se observa

que: cualquier cambio en el ngulo de
inyeccin, por encima o por debajo del
nominal, representa una disminucin en
la eficiencia trmina al freno y la eficiencia
mecnica; retardos en el momento de
la inyeccin implican disminuciones en
la temperatura pico; incrementos en la
relacin de compresin representan
incrementos en todos los parmetros
estudiados; e incrementos en la relacin
relativa aire combustible hacen que se
disminuya la presin media efectiva al
freno y la temperatura pico en el cilindro.

Cero-dimensional

Efecto de la velocidad y la
relacin de compresin en la
potencia al freno y la eficiencia
al freno de un motor Disel
mono cilndrico de 4 tiempos
alimentado con mezclas Disel/
Biodisel de aceite de Karanja.

Se observa un desempeo similar al

utilizar mezclas de Disel con el 20%vol
y 40%vol de Biodisel y un desempeo
superior al utilizar la mezcla Disel con
60%vol de Biodisel.

Cuasidimensional

Multidimensional

47

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

Autores /
Pas / Ao

Rajendra
et al. [99] /
India / 2010

Shi y Reitz
[27] / USA /
2010

Bueno et
al. [52] /
Brasil /
2010

Ortiz [101] /
Colombia /
2010

Colao
et al. [53]
/ Brasil /
2010

Tipo de Modelo

Estudio Realizado

Resultados

Efectividad del la combinacin

uso de biodisel de Jatropha y el
uso de materiales cermicos para
aislar la cmara de combustin
de un motor Disel (motor
LHR). Para esto se compar
el desempeo de las cuatro
combinaciones posibles: motor
LHR con Disel, motor LHR con
biodiesel, motor convencional
con Disel y motor convencional
con biodisel,

Se presenta una temperatura en cmara

mucho mayor cuando se utiliza biodisel
en un motor LHR, que cuando se utilizan
las dems combinaciones. El calor
liberado en la combustin y el trabajo
desarrollado son mayores al utilizarse
combustible Disel en el motor LHR. Las
prdidas de calor son mucho menores
en en un motor LHR que en un motor
convencional; mientras que las emisiones
de NO son mayores en un motor LHR
que en un motor convencional, lo cual
se cumple para los dos combustibles
analizados.

Multidimensional

Parmetros de inyeccin ptimos

de un motor de encendido por
compresin alimentado con
Disel, gasolina y combustible
E10, a condiciones de mediana
y alta carga.

Combustibles tipo gasolina presentan

mayor potencial para la combustin limpia
que el combustible Disel convencional.
En condiciones de carga media, la
efectividad del motor con gasolina se
ve afectada por el ngulo de la segunda
inyeccin y en el motor alimentado con
Disel se ve afectada por la cantidad
inyectada en la primera inyeccin.

Cero-dimensional

Impacto ocasionado por el uso

de biodisel de ster etlico de
soya en un motor Disel turbocargado por medio de un anlisis
energtico y exergtico. Para lo
cual se compar el desempeo
del motor alimentado con Disel
y mezclas de Disel con 5, 10, 15
y 20 y 30%vol de biodisel.

Se observa un mayor incremento en

la eficiencia trmica al freno al utilizar
mezcla de Disel con 20%vol de Biodisel
y un menor consumo especfico de
combustible al utilizar mezcla con 10%vol
de Biodisel.

Cuasidimensional

Efecto de la modificacin de
la viscosidad del combustible
en la liberacin de energa e
influencia de las propiedades
fsico-qumicas en la atomizacin
y evaporacin del combustible.

El fenmeno de rompimiento, donde se

manifiesta el efecto de la viscosidad de
lquido, del combustible incrementa el
rea superficial y disminuye el dimetro
medio de Sauter (SMD) de las gotas en
el chorro. No fue evidente la relacin
entre mayor rea, menor SMD y mayor
velocidad de evaporacin de las gotas.

Cero-dimensional

Perfil de temperatura del pistn

de un motor Disel alimentado
con Disel, biodisel y sus
mezclas.

El uso de diferentes mezclas Diselbiodisel no cambia el perfil de

temperaturas del pistn; sin embargo,
debido al menor poder calorfico de
biodisel, este ltimo hace que la
operacin del motor sea ms inestable.

Cuasidimensional

48

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

Autores /
Pas / Ao

Jagadish et
al. [100] /
India / 2011

Patil y
Akarte [55]
/ India /
2012

ipo de Modelo

Estudio Realizado

Resultados

Cuasidimensional

Efectos del uso de mezclas

Disel-biodisel y Disel-etanol
en el desempeo y emisiones
de un motor monocilndrico a
varias condiciones de carga y
porcentajes de recirculacin de
gases de escape (EGR).

El desempeo del motor es mejorado

al
utilizar
Disel
mezclado
con
biocombustibles en bajas cantidades.
Las emisiones de NOx se reducen
significativamente,
para
todos
los
combustibles, al utilizar el sistema EGR.
Las emisiones de holln se ven reducidas al
utilizar los biocombustibles considerados.

Cero-dimensional

Desempeo de un motor Disel

alimentado con mezclas Diselmetil ester de aceite de palma
(B0, B20, B60 y B100), por
medio de un anlisis del efecto
de la variacin de la relacin de
compresin en la presin pico,
la tasa de calor liberado y la
eficiencia trmica al freno, al usar
cada uno de los combustibles

Al incrementar la relacin de compresin

se presenta un incremento en la presin
pico y la eficiencia trmica al freno para
todas las mezclas analizadas. Para todos
los valores de relacin de compresin,
un incremento en la porcin de Biodisel
mezclado se traduce en un decremento
de la presin pico y un incremento de la
eficiencia trmica al freno.

Conclusin

matemtica de la forma de la ley de liberacin de

calor sin relacin directa con la fsica que controla el
proceso, mientras que otros describen y discretizan
el problema teniendo en cuenta la fsico-qumica
de los fenmenos ocurridos durante la combustin
de forma simplificada. En los ltimos aos se
han desarrollado varios modelos bajo este ltimo
enfoque con el fin de evaluar el desempeo de
motores Disel alimentados por biocombustibles
de una forma relativamente sencilla.
Por su parte los modelos cuasi-dimensionales,
que se presentan como un punto intermedio
entre los dos tipos anteriores, permiten calcular
los parmetros de desempeo del motor de un
modo mucho ms sencillo que con un modelo
multidimensional y, a diferencia de los modelos
cero-dimensionales, tiene el potencial para
calcular las emisiones del motor. Sin embargo,
presentan como inconveniente la necesidad de
ajustar algunos de sus coeficientes con resultados
experimentales en cada situacin, debido al
gran nmero de simplificaciones que requieren.
ltimamente la mayora de modelos del tipo cuasidimensional han sido modelos multizona, aunque
tambin se han presentado modelos de dos
zonas con el objetivo de predecir el desempeo
y emisiones de motores Disel alimentados con
biocombustibles.
En Colombia la mayor parte de los modelos de la
combustin han sido modelos cero-dimensionales,
muchos de ellos desarrollados en la universidad

En conclusin, los modelos de combustin Disel

se pueden clasificar en tres grandes grupos:
modelos multidimensionales, modelos cuasidimensionales y modelos cero-dimensionales.
Los primeros son los ms adecuados cuando se
busca tener una descripcin detallada del proceso
de combustin Disel, mientras los segundos y los
terceros son idneos cuando se desea obtener
una descripcin global de ste.
Los modelos multidimensionales han venido
evolucionando desde los aos 70 hasta llegar a los
actuales cdigos CFD comerciales, ellos tienen el
potencial para describir en detalle los fenmenos
ocurridos durante la combustin y tienen en
cuenta la distribucin espacial del proceso, pero
presentan la desventaja de ser sensibles a los submodelos utilizados, algunos de los cuales no estn
completamente definidos, adems de consumir
mayor tiempo en la simulacin del proceso, lo cual
los hace los menos eficientes desde el punto de
vista computacional.
Los modelos cero-dimensionales son los ms
sencillos, y por ende los ms eficientes a la hora
de calcular los parmetros de desempeo del
motor, sin embargo no describen en detalle los
subprocesos que ocurren durante el proceso
global de combustin Disel y son incapaces de
calcular las emisiones de contaminantes en el
motor. Algunos de ellos se basan en la descripcin

49

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

de Antioquia, aunque tambin se han elaborados

trabajos desde el enfoque cuasi-dimensional y
multidimensional.

Fluid Flow at All Speeds. Los Alamos, New

Mexico: Los Alamos National Laboratory;
1979. Report No. LA-8929-MS.
[10] Gosman AD, Johns RJR. Computer Analysis
of Fuel-Air Mixing in Direct-Injection Engines.
Warrendale, PA: SAE International; 1980
Feb. Report No. 800091.
[11] Syed SA, Bracco FV. Further Comparisons
of Computed and Measured DividedChamber Engine Combustion. Warrendale,
PA: SAE International; 1979 Feb. Report
No. 790247.
[12] Bracco FV. Applications of Steady-State
Spray Equations to Combustion Modeling.
AIAA J. 1974;12(11):153440.
[13] Harlow FH, Amsden AA. Numerical
calculation of multiphase fluid flow. J Comput
Phys. 1975;17(1):1952.
[14] Crowe CT. A computational model for the
gas-droplet flow field in the vicinity of an
atomizer. En: 11th JANAF Symposium; 1974
Oct 21-22; Northridge, Estados Unidos de
Amrica. Northridge: Combustion Institute;
1974. p. 18.
[15] Crowe C, Stock D, Sharma M. The
particle-source-in
cell/PSI-CELL/model
for gas-droplet flows. ASME J Fluids Eng.
1977;99:32532.
[16] Dukowicz JK. A particle-fluid numerical
model for liquid sprays. J Comput Phys.
1980;35(2): 22953.
[17] Gosman AD, Kralj C, Marooney CJ,
Theodossopouls P. The development of
the SPEED code for Diesel combustion
simulation. Florencia Italia: ATA-DMTI; 1992.
[18] Gosman AD, Kralj C, Marooney CJ,
Theodossopouls P. Development strategies
for Diesel combustion simulation using the
SPEED code. Institution of Mechanical
Engineers; 1992. Report No. Paper
C448/035.
[19] Fluent. FLUENT AND RAMPANT: Users
Guide. Fluent Incorporated; 1995.
[20] Amsden AA. KIVA-II: a computer program
for chemically reactive flows with sprays.
Los Alamos, NM: Los Alamos National
Laboratory; 1985.
[21] Amsden AA. KIVA-III: A KIVA Program
with Block-structured Mesh for Complex
Geometries.
Los
Alamos
National
Laboratory; 1993.
[22] Gao Z, Wagner RM, Sluder CS, Daw CS,
Green Jr. JB. Using a phenomenological

Agradecimientos
Los autores expresan sus agradecimientos a
las Vicerrectora de Investigacin y Extensin
de la Universidad Industrial de Santander por la
financiacin del proyecto Produccin de energa a
partir de residuos provenientes de la transformacin
de la biomasa (Proyecto No. 5451), del cual hace
parte este trabajo.
Referencias
[1]

[2]

[3]

[4]

[5]

[6]
[7]

[8]

[9]

Ramadhas AS, Jayaraj S, Muraleedharan

C. Theoretical modeling and experimental
studies on biodiesel-fueled engine. Renew
Energ. 2006;31(11):181326.
Fenollosa C. Descripcin del proceso
de combustin Disel. En: Modelado
fenomenolgico del proceso de combustin
por difusin Disel. Espaa: Reverte; 2005.
p. 7-29.
Colao MJ, Teixeira CV, Dutra LM.
Thermodynamic simulation and optimization
of Diesel engines operating with Disel and
biodiesel blends using experimental data.
inverse Probl Sci Eng. 2010;18(6):787812.
Heywood J. Modeling Real Engine Flow
and Combustion Processes. En: Internal
Combustion Engine Fundamentals. EE.UU:
McGraw-Hill; 1988. p. 748-816.
Ganesan V. Introduction. En: Computer
Simulation Of Compression-Ignition Engine
Processes. India: Universities Press; 2000.
p. 1-2.
Watson N. Turbocharging the Internal
Combustion Engine. John Wiley & Sons Inc;
1982.
Bracco FV. Introducing a New Generation of
More Detailed and Informative Combustion
Models. Warrendale, PA: SAE International;
1975 Feb. Report No. 751187.
Escalante H, Orduz J, Zapata H, Cardona M,
Duarte M. Atlas del Potencial Energtico de
la Biomasa Residual en Colombia. Colombia:
Ministerio de Minas y Energa; 2011.
Butler TD, Cloutman LD, Dukowicz J k.,
Ramshaw JD. CONCHAS: An Arbitrary
Lagrangian-Eulerian Computer Code for
Multi-component
Chemically
Reactive

50

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

computer model to investigate advanced

combustion trajectories in a CIDI engine.
Fuel. 2011;90(5):190718.
[23] Peng Krrholm F, Tao F, Nordin N. ThreeDimensional Simulation of Diesel Spray
Ignition and Flame Lift-Off Using OpenFOAM
and KIVA-3V CFD Codes. Warrendale, PA:
SAE International; 2008 Apr. Report No.
2008-01-0961.
[24] Gong Y, Kaario O, Tilli A, Larmi M, Tanner
FX. A Computational Investigation of
Hydrotreated Vegetable Oil Sprays Using
RANS and a Modified Version of the RNG
k - Model in OpenFOAM. Warrendale, PA:
SAE International; 2010 Apr. Report No.
2010-01-0739.
[25] Zheng QP, Zhang HM, Zhang DF. A
computational study of combustion in
compression ignition natural gas engine
with separated chamber. Fuel. 2005;84(12
13):151523.
[26] Jayashankara B, Ganesan V. Effect of fuel
injection timing and intake pressure on
the performance of a DI Diesel engine A
parametric study using CFD. Energ Convers
Manage. 2010;51(10):183548.
[27] Shi Y, Reitz RD. Optimization of a heavyduty compressionignition engine fueled
with Diesel and gasoline-like fuels. Fuel.
2010;89(11):341630.
[28] Gutirrez E. Modelado numrico de la
atomizacin y evaporacin de un chorro
Diesel mediante anlisis CFD (Tesis de
Maestra) Medelln, Colombia: Universidad
de Antioquia; 2008.
[29] Agudelo J, Agudelo A, Benjumea P. Estudio
de chorros Disel usando mecnica de
fluidos computacional. Rev fac ing univ
Antioquia. 2009;(49):619.
[30] Delgado A. Simulacin numrica de los
procesos al interior del cilindro en un motor
Disel mediante anlisis CFD (Tesis de
Maestra) Medelln, Colombia: Universidad
de Antioquia; 2010.
[31] Arrgle J, Lpez JJ, Garca
JM, Fenollosa C.
Development of a zero-dimensional Diesel
combustion model. Part 1: Analysis of the
quasi-steady diffusion combustion phase.
Appl Therm Eng. 2003;23(11):130117.
[32] Rakopoulos
CD,
Antonopoulos
KA,
Rakopoulos
DC.
Development
and
application of multi-zone model for
combustion and pollutants formation in

direct injection Diesel engine running

with vegetable oil or its bio-Diesel. Energ
Convers Manage. 2007;48(7):1881901.
[33] Chmela FG, Pirker GH, Wimmer A. Zerodimensional ROHR simulation for DI Diesel
engines A generic approach. Energ
Convers Manage. 2007;48(11):294250.
[34] Galindo J, Hernandez J. Modelado de
Motores. En: Motores de combustin interna
alternativos. Payri F, Desantes JM. Espaa:
Reverte: 2011. p. 746-796.
[35] Vibe II. Brennverlauf und Kreisproze von
Verbrennungsmotoren. VEB Verlag Technik;
1970.
[36] Ghojel JI. Review of the development
and applications of the Wiebe function: A
tribute to the contribution of Ivan Wiebe
to engine research. Int J Engine Res.
2010;11(4):297312.
[37] Dyechenko NK, Magidovich LY, Pugachiuv
BP. Approximation of heat release
characteristics in the cylinders of Diesel
engines.
Energomashinostroyeniye.
1969;310:73.
[38] Dyechenko NK, Magidovich LY, Pugachiuv
BP. Determination of main parameters of heat
release characteristics during combustion in
Diesel engines. Energomashinostroyeniye.
1970;316:54.
[39] Ghojel JI. Investigation of cycle processes of
direct injection Disel engines with cylindrical
piston bowl (Tesis Doctoral). Mosc, Rusia:
Moscow Automobile and Road Institute
(MADI); 1974.
[40] Ghojel JI. A Study of Combustion Chamber
Arrangements and Heat Release in D.I.
Diesel Engines. Warrendale, PA: SAE
International; 1982 Feb. Report No. 821034.
[41] Watson N, Pilley AD, Marzouk M. A
Combustion
Correlation
for
Disel
Engine Simulation. Warrendale, PA: SAE
International; 1980 Feb. Report No. 800029.
[42] Miyamoto N, Chikahisa T, Murayama T,
Sawyer R. Description and Analysis of
Diesel Engine Rate of Combustion and
Performance Using Wiebes Functions.
Warrendale, PA: SAE International; 1985
Feb. Report No. 850107.
[43] Witt H, Hassenforder M, Gissinger GL.
Modelling and Identification of a Diesel
Combustion Process with the Downhill
Gradient Search Method. Warrendale, PA:
SAE International; 1995 Feb. Report No.

51

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

950854.
[44] Bilcan A, Tazerout M, Le Corre O, Ramesh
A. Ignition Delay in Dual Fuel Engines: An
Extended Correlation for Gaseous Fuels.
ASME Conf Proc. 2001;2001(16699a):4955.
[45] Canova M, Garcin R, Midlam-Mohler S,
Guezennec Y, Rizzoni G. A control-oriented
model of combustion process in a HCCI
Diesel engine. En: Proceedings of the
American Control Conference; 2005 jun
8-10. Portland, Estados Unidos; Portland:
IEEE; 2005. p. 4446-51.
[46] Galindo J, Lujn JM, Serrano JR, Hernndez
L. Combustion simulation of turbocharger
HSDI Diesel engines during transient
operation using neural networks. Appl
Therm Eng. 2005;25(56):87798.
[47] Way RJB. Methods for Determination
of Composition and Thermodynamic
Properties of Combustion Products for
Internal Combustion Engine Calculations. P
I Mech Eng. 1976;190(1):68797.
[48] Hardenberg HO, Hase FW. An Empirical
Formula for Computing the Pressure Rise
Delay of a Fuel from Its Cetane Number
and from the Relevant Parameters of DirectInjection Diesel Engines. Warrendale, PA:
SAE International; 1979 Feb. Report No.
790493.
[49] Gardner TP, Henein NA. Diesel Starting: A
Mathematical Model. Warrendale, PA: SAE
International; 1988 Feb. Report No. 880426.
[50] Harris HD, Pearce F. A universal mathematical
model of Diesel engine performance. J Agr
Eng Res. 1990;47:16576.
[51] Gogoi TK, Baruah DC. A cycle simulation
model for predicting the performance of a
Disel engine fuelled by Diesel and bioDisel
blends. Energy. 2010;35(3):131723.
[52] Bueno AV, Velsquez JA, Milanez LF. Heat
release and engine performance effects of
soybean oil ethyl ester blending into Diesel
fuel. Energy. 2011;36(6):390716.
[53] Colao MJ, Teixeira CV, Dutra LM. Thermal
analysis of a Diesel engine operating
with
Diselbiodiesel
blends.
Fuel.
2010;89(12):374252.
[54] Payri F, Olmeda P, Martn J, Garca A. A
complete 0D thermodynamic predictive
model for direct injection Diesel engines.
Appl Energ. 2011;88(12):463241.
[55] Patil S, Akarte. Investigation on Effect
of Variation in Compression Ratio

on
Performance
and
Combustion
Characteristics of C.I Engine Fuelled With
Palm Oil Methyl Ester (POME) and its Blends
by Simulation. GJRE. 2012;12(2):3441.
[56]
Rakopoulos
CD,
Hountalas
DT.
Development and Validation of a 3-D MultiZone Combustion Model for the Prediction
of DI Diesel Engines Performance and
Pollutants Emissions. Warrendale, PA: SAE
International; 1998 Feb. Report No. 981021.
[57] Agudelo J, Corredor L, Hernandez J.
Simulacin del Proceso de Combustin en
Motores de Encendido por Compresin.
Energtica. 2000;23:10311.
[58] Agudelo A, Agudelo J, Benjumea P. Exergy
diagnosis of the combustion process in a
Diesel engine. Rev fac ing univ Antioquia.
2008;(45):4153.
[59] Barbosa J, Salcedo J. Herramienta educativa
en la web basada en la simulacin de los
procesos termofluidodinmicos de motores
de combustin interna, SIMTA 1.0. (Tesis
de Pregrado) Bucaramanga, Colombia:
Universidad Industrial de Santander; 2007.
[60] Rodrguez E. Software para la simulacin
de los ciclos termodinmicos de los motores
de combustin interna de cuatro tiempos
SICICLOTER 1.0 (Tesis de Pregrado)
Bucaramanga,
Colombia:
Universidad
Industrial de Santander; 2003.
[61] Chacn JL. Identificacin de un motor
Disel. En: IX Congreso Latinoamericano
de Control Automtico; 2000 nov 1-3.
Cali,Colombia; Cali: Universidad Autnoma
de Occidente; 2000. p. 1-5
[62] Probert RP. The influence of spray particle
size and distribution in the combustion of oil
droplets. Philos Mag. 1946;37(265):94105.
[63] Tanasawa PH. On the combustion rate of a
group of particles. Tech Report Tohoko Univ.
1953;18:6174.
[64] Austen AEW, Lyn W-T. Relation between
Fuel Injection and Heat Release in a
Direct-Injection Engine and the Nature of
the Combustion Processes. P I Mech Eng.
1960;14(1):4762.
[65] Lyn WT. Study of burning rate and nature
of combustion in Diesel engines. Symp (Int)
Combust. 1963;9(1):106982.
[66] Cook HA. Appraisal of Effects of Operating
Conditions and Engine Design on
Combustion in Reciprocating Engines by
Engine Cycle Calculations. Warrendale, PA:

52

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

SAE International; 1963 Jan. Report No.

630083.
[67] Cook HA. Diesel Engine Cycle Analysis
of Relationships of Fuel Injection to Fuel
Compression Ignition Characteristics and
Best Fuel Utilization. Society of Automotive
Engineers; 1965.
[68] Nagao F, Ikegami M, Oshima K. An Analysis
of Combustion Knock in a Diesel Engine.
Bull JSME. 1967;10(39):532-42.
[69] Shipinski J, Myers PS, Uyehara OA. A
Spray-Droplet Model for Diesel Combustion.
P I Mech Eng. 1969;184(10):2835.
[70] Shipinski J, Uyehara OA, Myers PS.
Experimental correlation between rate of
injection and rate of heat release in a Diesel
engine. 1968. Report No. 68-DGP-11.
[71] Whitehouse ND, Way R. Rate of Heat
Release in Diesel Engines and Its Correlation
with Fuel Injection Data. P I Mech Eng.
1969;184(10):1727.
[72] Whitehouse ND, Way RJB. A Simple Method
for the Calculation of Heat Release Rates in
Diesel Engines Based on the Fuel Injection
Rate. Warrendale, PA: SAE International;
1971 Feb. Report No. 710134.
[73] Whitehouse ND, Sareen BK. Prediction
of Heat Release in a Quiescent Chamber
Diesel Engine Allowing for Fuel/air Mixing.
Society of Automotive Engineers; 1974.
[74] Grigg HC, Syed MH. The Problem of
Predicting Rate of Heat Release in Diesel
Engines. P I Mech Eng. 1969;184(10):192
202.
[75] Khan IM, Greeves G, Probert DM. Prediction
of soot and nitric oxide concentrations
in Diesel engine exhaust. Institution of
Mechanical Engineers; 1971 p. 205 217.
[76] Abramovich GN. The Theory of Turbulent
Jets. Schindel L, editor. The MIT Press; 2003.
[77] Adler D, Lyn W-T. The Evaporation and
Mixing of a Liquid Fuel Spray in a Diesel Air
Swirl. P I Mech Eng. 1969;184(10):17180.
[78] Rife J, Heywood JB. Photographic and
performance studies of Diesel combustion
with a rapid compression machine. 1974.
[79] Chiu WS, Shahed SM, Lyn WT. A Transient
Spray Mixing Model for Diesel Combustion.
Warrendale, PA: SAE International; 1976
Feb. Report No. 760128.
[80] Hiroyasu H, Kadota T, Arai M. Development
and Use of a Spray Combustion Modeling
to Predict Diesel Engine Efficiency and

Pollutant Emissions: Part 1 Combustion

Modeling. Bull JSME. 1983;26(214):569-75.
[81] Dent JC, Mehta PS. Phenomenological
combustion model for a quiescent chamber
Diesel engine. 1981.
[82] Dent JC, Mehta PS, Swan J. A predictive
model for automotive DI Diesel engine
performance and smoke emissions.
Institution of Mechanical Engineers; 1982 p.
237 245. Report No. Paper C126/82.
[83] Kyriakides SC, Dent JC, Mehta PS.
Phenomenological
Diesel
combustion
model including smoke and NO emission.
Estados Unidos de Amrica: Society of
Automotive Engineers; 1986.
[84] Tinaut FV. Contribucin al estudio del
proceso de combustin en motores de
encendido por compresin de inyeccin
directa (Tesis Doctoral) Valencia, Espaa:
E.T.S. Ingenieros Industriales. Universidad
Politcnica de Valencia; 1986.
[85] Lapuerta M. Estudio fenomenolgico de la
combustin en motores Disel rpidos de
inyeccin directa (Tesis Doctoral) Valencia,
Espaa: E.T.S. Ingenieros Industriales.
Universidad Politcnica de Valencia; 1988.
[86] Desantes JM, Lapuerta M, Tinaut FV.
Combustion dans les moteurs Diesel
rapides injection directe. Modle bas
sur linteraction entre le jet et le mouvement
de lair dans le cylindre. Entropie.
1989;148:5161.
[87] Bazari Z. A DI Diesel Combustion and
Emission Predictive Capability for Use in
Cycle Simulation. Warrendale, PA: SAE
International; 1992 Feb. Report No.: 920462.
[88] Cui Y, Deng K, Wu J. A direct injection
Diesel combustion model for use in transient
condition analysis. P I Mech Eng D-J Aut.
2001;215(9):9951004.
[89] Kouremenos DA, Rakopoulos CD, Hountalas
DT. Multi-Zone Combustion Modelling for
the Prediction of Pollutants Emissions
and Performance of DI Diesel Engines.
Warrendale, PA: SAE International; 1997
Feb. Report No. 970635.
[90] Ottikkutti P, Van Gerpen J, Cui KR. Multizone
Modeling of a Fumigated Diesel Engine.
Warrendale, PA: SAE International; 1991
Feb. Report No. 910076.
[91] Shenghua L, Hwang JW, Park JK, Kim
MH, Chae JO. Multizone Model for DI
Diesel Engine Combustion and Emissions.

53

rev.ion. 2013;26(1):41-54. Bucaramanga (Colombia).

Warrendale, PA: SAE International; 1999

Aug. Report No. 1999-01-2926.
[92] Torkzadeh DD, Langst W, Kiencke U.
Combustion and Exhaust Gas Modeling of a
Common Rail Diesel Engine - an Approach.
Warrendale, PA: SAE International; 2001
Mar. Report No. 2001-01-1243.
[93] Bi X, Yang M, Han S, Ma Z. A Multi-Zone
Model for Diesel Spray Combustion .
Warrendale, PA: SAE International; 1999
Mar. Report No. 1999-01-0916.
[94] Rakopoulos
CD,
Antonopoulos
KA,
Rakopoulos DC, Hountalas DT. Multi-zone
modeling of combustion and emissions
formation in DI Diesel engine operating on
ethanolDiesel fuel blends. Energ Convers
Manage. 2008;49(4):62543.
[95] Komninos NP, Rakopoulos CD. Numerical
Investigation into the Formation of CO
and Oxygenated and Nonoxygenated
Hydrocarbon Emissions from Isooctaneand Ethanol-Fueled HCCI Engines. Energ
Fuel. 2010;24(3):165567.
[96] Komninos NP, Kosmadakis GM. Heat
transfer in HCCI multi-zone modeling:
Validation of a new wall heat flux correlation
under motoring conditions. Appl Energ.
2011;88(5):163548.
[97] Ganapathy T, Gakkhar P, Murugesan K.
An analytical and experimental study of
performance on jatropha biodiesel engine.
Therm Sci. 2009;13(3):6982.
[98] Ganapathy T, Murugesan K, Gakkhar RP.

Performance optimization of Jatropha

biodiesel engine model using Taguchi
approach. Appl Energ. 2009;86(11):247686.
[99] Rajendra Prasath B, Tamilporai P, Shabir
MF. Analysis of combustion, performance
and emission characteristics of low heat
rejection engine using biodiesel. Int J
Therm Sci. 2010;49(12):248390.
[100] Jagadish D, Kurmar R, K. M. Zero
dimensional simulation of Combustion
Process of a DI Diesel engine fuelled with
Biofuels. WASET. 2011;80:81925.
[101] Ortiz Valenzuela F. Modelado del efecto
de la viscosidad en la combustin de la
mezcla B20 de biodisel de palma africana
tenera en motores de combustin interna
encendidos por compresin / Modeling the
effect of viscosity on the combustion of
biodiesel B20 blend of African palm in ICECI (Tesis de Maestra) Bogot, Colombia:
Universidad Nacional de Colombia; 2010.
[102] Rodrguez Rodrguez AM. Modelacin de
la generacin de material particulado en
funcin de la composicin del combustible/
Modeling of particular metter generation in
function to the fuel composition (Tesis de
Maestra) Bogot, Colombia: Universidad
Nacional de Colombia; 2010.
[103] Villegas A. Estudio teorico-experimental de
las emisiones de NOx y material particulado
en motores Disel operando con biodiesel
(Tesis de Maestra) Medelln, Colombia:
Universidad de Antioquia; 2011.

54