tesisUPV3034 PDF
tesisUPV3034 PDF
tesisUPV3034 PDF
Presentada por
CARLOS ALBERTO ROMERO PIEDRAHITA
En el
DEPARTAMENTO DE MQUINAS Y MOTORES TRMICOS de la
UNIVERSIDAD POLITCNICA DE VALENCIA
Para la obtencin del grado de
DOCTOR
Valencia, mayo de 2009
TESIS DOCTORAL
Realizada por:
Dirigida por:
TRIBUNAL CALIFICADOR
Presidente:
Secretario:
Vocales:
Vocales Suplentes:
Dr. D. Andrs Melgar Bachiller
Dr. D. Manuel Valds Del Fresno
Resumen
El trabajo desarrollado en la presente tesis est encaminado a estudiar los
fenmenos de la transmisin de calor y balance trmico en los motores de
combustin interna alternativos, tanto en el cilindro como en el sistema de
refrigeracin. El principal objetivo de la tesis ha sido la utilizacin y
extensin de un programa computacional de prediccin de temperaturas de
paredes del motor y de dos programas comerciales de modelado de sistemas
termohidrulicos para estudiar los flujos de calor dentro y fuera de la
estructura del motor, para evaluar el balance trmico de los motores en estado
estacionario y en condiciones transitorias y para modelar, simular y validar
diseos nuevos de sistemas de refrigeracin de motores para locomotoras.
Para lograr este objetivo se han realizado varias tareas que se explican a lo
largo de este documento, y que se sintetizan en una revisin de los modelos
de transferencia de calor utilizados en los clculos de flujos de calor en las
diferentes condiciones de contorno presentes en los motores de combustin;
la expansin de las capacidades de un modelo nodal para obtener de l los
flujos de calor discriminados entre las partes metlicas del motor y los calores
intercambiados entre los fluidos del sistema y las partes metlicas; el estudio
experimental de la evolucin de las temperaturas del motor durante el
calentamiento bajo un ciclo de conduccin con el fin de cuantificar el balance
trmico en este modo de operacin y proponer alternativas de bajo costo para
mejorar ese tiempo de calentamiento; el modelado, simulacin, validacin y
evaluacin del impacto sobre el calentamiento, el consumo y las emisiones
del motor, de alternativas propuestas en el diseo de los sistemas de
refrigeracin; y finalmente, la elaboracin y aplicacin de una metodologa
de prediseo, modelado y simulacin de sistemas de refrigeracin para
aplicaciones especiales.
Resum
El treball desenvolupat en la present tesi est encaminat a estudiar els
fenmens de transmissi de calor i el balan trmic en els motors de
combusti interna alternatius, tant en el cilindre com en el sistema de
refrigeraci. El principal objectiu de la tesi ha sigut la utilitzaci i extensi
dun programa computacional de predicci de temperatures de les parets del
motor i de dos programes comercials de modelat de sistemes termohidrulics
per estudiar els fluxos de calor dins i fora de lestructura del motor, per
ii
Abstract
In the thesis at hand the study of heat tranfer phenomena and thermal balance
inside the engine structure, as well as in the cooling system is addressed. It
has been the main goal of this thesis the use and extension of a model for the
prediction of combustion chamber wall temperatures, along with the
application of two commercial computational thermohydraulic programs, to
the study of heat flows inside and outside the engine structure in order to
assess the engine energy balance under steady and transitory operation, and
also to design or redesign new cooling systems for locomotive engines. To
accomplish the aforementioned goal some steps explained along this
document, were followed. Those tasks consisted on: a thorough review of the
heat transfer models used to calculate the heat flows in all the internal
combustion engine boundaries; the enhacement of a predictive wall
temperature model to detail all the heat flows between metallic parts and
boundary media; the experimental study of the engine warm-up under a
driving cycle to study its energy balance, and to propose low cost
modifications to the original cooling system intended to accelerate the engine
warm-up; the modelling and assessment of the proposed modifications in
terms of their impact on the engine warm-up time, as well as on the fuel
iii
iv
Amispadres
AValentinayCarlosAlberto
vi
vii
Agradecimientos
Esta tesis ha sido realizada durante mi trabajo como docente de la Escuela de
Tecnologa Mecnica de la Universidad Tecnolgica de Pereira, a ella debo
todo el apoyo econmico, y de ello estoy profundamente agradecido.
Ante todo, quiero agradecer a mis directores Antonio Jos Torregrosa y a
Pablo Csar Olmeda por la asesora brindada y por vincularme a su equipo de
trabajo. El Dr. Torregrosa me present en la primera charla que tuvimos, de
una manera insuperablemente clara, sinttica, precisa y completa, el problema
de la transferencia de calor entre los gases y las paredes de la cmara de
combustin en los motores, y este ha sido el inicio del trabajo cuyos
resultados presento, pero que no culmina ah, porque aunque la charla fue
breve, proyect todas las tareas sobre las que quisiera seguir trabajando y
para las que confo contar con su profundo conocimiento y su sabia visin.
Agradezco una vez ms al Dr. Torregrosa su excelente disposicin para
conmigo, su inters en mi trabajo, sus acertadas recomendaciones y tambin
por compartir sus ideas y pensamiento. A Pablo, incansable trabajador, le
agradezco grandemente el acompaamiento continuo, las siempre oportunas
observaciones y el rigor puesto en la lectura y correccin de esta tesis.
Agradezco a los responsables del Departamento de Mquinas y Motores
Trmicos, quienes me han facilitado los medios necesarios para llevar a cabo
los trabajos de investigacin que se presentan en la tesis, as como al personal
de secretara por ayudarme en todos los trmites administrativos.
Muy especialmente me gustara dar las gracias a Juanma por la ayuda
prestada en todos los temas informticos relacionados con la programacin en
C++. Por la grata compaa durante estos aos en el departamento quiero dar
las gracias a la gente de la lnea de Transferencia de Calor y, particularmente,
a Jaime Martin y Bart Degraeuwe. No puedo menos que agradecer la ayuda
de Benjamn Pla, Jos Manuel Lujan, Salvador Torr en la facilitacin de los
trabajos de medicin de temperaturas. Los trabajos realizados sobre los
sistemas de refrigeracin de las locomotoras no se habran podido llevar a
cabo sin la colaboracin de la empresa Vossloh y en particular de Toni y
Jordi. Gracias a los estudiantes Edouard, Juliene, Etienne por sus
contribuciones en este trabajo.
viii
Agradezco a todos mis amigos que no menciono porque son muchos, que me
han apoyado en esta aventura de venirme a hacer lo que me gusta. Es muy
importante para m reconocer el apoyo, an en la distancia, de mis
compaeros de trabajo en la UTP, pero sobre todo de Yamid, Hctor Fabio y
Gabriel. Gracias a mis hermanos por querer que emprendiera esta etapa y
sobre todo gracias a John Jairo y a Rosa.
Finalmente, lo ms importante, de la manera ms profunda, agradezco a mi
compaera Luz Adriana y a mis hijos por haberme aguantado durante estos
aos y por aceptar mi falta de atencin en ciertos momentos.
ix
ndice General
Lista de smbolos............... xv
Captulo 1. Introduccin..........
1.1. Motivacin......................
1.2. Objetivos.....
1.3. Antecedentes.......
1.4. Metodologa y planteamiento..
1.5. Bibliografa.....
Captulo 2. Estado del arte de la transferencia de calor en los
M.C.I.A..
2.1. Introduccin.....
2.2. La transferencia de calor en el motor......
2.3. Transferencia de calor entre los gases del cilindro y las paredes de
la cmara de combustin.........................
2.3.1. Modelos globales..
2.3.2. Modelos zonales y locales....
2.3.3. Modelos multidimensionales....
2.3.4. Transferencia de calor radiativa...............
2.4. Transferencia de calor en las pipas de admisin y escape...
2.5. Transferencia de calor del lado del refrigerante..
2.5.1. Modelos empricos...
2.5.2. Modelado CFD.
1
3
8
9
13
15
19
21
22
23
26
32
46
51
62
67
69
72
2.6. Transferencia de calor entre el aceite y el motor..... 73
2.6.1. Transferencia de calor entre el pistn y el aceite.. 73
2.6.2. Transferencia de calor entre el aceite y la estructura (bloque y
culata)......... 77
2.7. Transferencia de calor entre el motor y el ambiente.... 78
2.8. Conclusiones.... 81
2.9. Bibliografa.......... 83
95
97
99
100
102
103
105
107
107
117
118
123
130
134
135
136
136
138
141
143
143
145
xi
motor..
3.5.1. Prediccin de los flujos de calor al refrigerante...
3.5.2. Prediccin del comportamiento trmico del motor bajo
condiciones de calentamiento
3.5.3. Observaciones..
3.6. Conclusiones
3.7. Bibliografa..
Captulo 4. Consideraciones y herramientas de diseo y modelado
de los sistemas de refrigeracin de los motores
Resumen......
4.1. Introduccin.
4.2. Prediseo de los sistemas de refrigeracin de los motores..
4.2.1. Clculo preliminar de los componentes del sistema de
refrigeracin...
4.2.2. Metodologa de diseo de los sistemas de refrigeracin..
4.2.3. Clculo de los radiadores de los sistemas de refrigeracin..
4.3. Modelado de los sistemas de refrigeracin..
4.3.1. Modelos simplificados de la temperatura del motor a nivel de
sistema
4.3.2. Modelos a nivel de circuito del sistema de refrigeracin del
motor..
4.4. Herramientas de modelado unidimensional.
4.5. Conclusiones
4.6. Bibliografa..
Captulo 5. Modelado, simulacin y validacin del sistema de
refrigeracin de un motor HSDI para reducir su
tiempo de calentamiento....................................................
5.1. Introduccin.
5.2. Modelado trmico del motor
5.2.1. Descripcin y caracterizacin del sistema de refrigeracin.
5.2.2. Modelo 1D del sistema de refrigeracin original del motor.
5.3. Modelado de las modificaciones al sistema de refrigeracin del
151
151
154
158
159
161
165
167
167
169
170
184
187
196
197
199
208
212
213
217
219
221
221
226
xii
231
236
237
238
241
243
245
249
251
251
254
254
258
266
269
269
275
280
287
289
290
xiii
294
297
299
303
304
321
321
327
330
xiv
xv
LISTA DE SMBOLOS
Latinos
A
rea
Calor especfico a presin constante de un gas
Cp
cm
cu
D, d
f
F
h
h
i
k
K
K
L
M
m
m&
n
N
p
pme
Q
Q&
Dimetro
Factor de friccin
rea
Coeficiente de transmisin de calor
Entalpa
Nmero de cilindros del motor
Conductividad trmica
Conductancia trmica, constante
Coeficiente de momento, coeficiente de prdidas de presin
Longitud caracterstica
Par
Masa
Flujo msico
Velocidad de motor
Potencia
Presin
Presin media efectiva del motor.
Calor
S
T
t
u
xvi
U
V
V&
W
Griegos
Coeficiente de difusividad
Eficiencia
Viscosidad dinmica
Viscosidad cinemtica
Densidad
Constante de Stefan-Boltzmann
Velocidad angular
Incremento
Subndices
AC
Relativo al enfriador de aire
adm Relativo a la admisin
aire
Relativo al aire
bl
Relativo al bloque del motor
comb Relativo al combustible
e
Relativo a las condiciones efectivas de motor
en
Relativo a las condiciones de entrada
esc
Relativo al escape del motor
equiv Relativo a valores equivalentes
fr
Frontal
g
Relativo a las condiciones del gas
i
Relativo a las condiciones en el interior del cilindro
rad
Relativo al radiador, relativo a la radiacin
xvii
refr
sal
w
Relativo al refrigerante
Relativo a las condiciones de salida
Relativo a las condiciones en la pared
Relativo a condiciones basadas en corriente libre de un fluido
Superndices
en
Relativo a las condiciones de entrada
sal
Relativo a las condiciones de salida
Grupos adimensionales
Bi
Nmero de Biot
Fo
Nmero de Fourier
j
Nmero de Colbourn
Nu
Nmero de Nusselt
Pr
Nmero de Prandtl
Ra
Nmero de Rayleigh
Re
Nmero de Reynolds
St
Nmero de Stanton
Siglas
CFD
EDC+EUDC
EGR
EVC
EVO
HSDI
ITD
IVC
IVO
M.C.I.A.
NEDC
NTU
xviii
UPV
WAM
ACT
Captulo 1
Introduccin
1.1. Motivacin ................................................................................................ 3
1.2. Objetivos ................................................................................................... 8
1.3. Antecedentes ............................................................................................. 9
1.4. Metodologa y planteamiento.................................................................. 13
1.5. Bibliografa.............................................................................................. 15
1. Introduccin
1.1. Motivacin
1.1. Motivacin
La mayor parte de los esfuerzos por mejorar los sistemas de potencia de los
medios de transporte se ha enfocado a mejorar sus caractersticas dinmicas,
reducir las emisiones y reducir tambin el consumo de combustible mediante
desarrollos aplicados sobre todo a los sistemas de combustible, los sistemas
de intercambio de gases, los mecanismos de distribucin de gases, los
sistemas de inyeccin y los sistemas de control. Los incrementos en
ganancias logradas con los diferentes sistemas se han ido haciendo ms
pequeos y slo en los ltimos aos se ha vuelto la mirada hacia los sistemas
de refrigeracin, el eslabn olvidado en los sistemas de control de los
sistemas de potencia, como foco potencial de perfeccionamiento del motor y
ganancias traducibles en ahorros energticos y reduccin de emisiones.
Con los avances tecnolgicos aplicados a los sistemas de transporte, ya no es
suficientemente preciso hablar de sistemas de refrigeracin del motor, sino de
sistemas de gestin trmica por la complejidad de los fluidos que participan
de los intercambios de calor y la variedad de los componentes, algunos de
ellos no pertenecientes al grupo funcional del motor. Es por esto que a lo
largo de esta tesis se entender por sistema de gestin trmica al conjunto de
sistemas dedicado al intercambio de calor en los sistemas de potencia de los
vehculos, acoplados directa o indirectamente con el sistema bsico de
refrigeracin del motor.
Las exigencias de reduccin de emisiones han conducido a la implementacin
de sistemas de recirculacin de gases refrigerados; el aumento de la potencia
especfica de los motores ha hecho obligatorio en casi todos los motores
modernos el montaje de refrigeradores del aceite del motor y, desde hace
poco tambin, de la transmisin. La bsqueda de una mayor potencia ha
conducido a la aplicacin de tecnologas de turboalimentacin y
sobrealimentacin con al menos una etapa de refrigeracin del aire fresco
comprimido antes de entrar a los cilindros. Se cuentan as al menos tres
intercambiadores de calor que antes no existan y que se integran en los
sistemas convencionales de refrigeracin, haciendo necesaria la
implementacin de sistemas de control de flujos y temperaturas (de
refrigerantes, gases, aire de admisin y aceite), que se complementen con los
sistemas de control del motor, de las emisiones y de la transmisin.
1. Introduccin
1.1. Motivacin
1. Introduccin
1.1. Motivacin
Figura 1.2. Diagrama de bloques del modelado dinmico del vehculo con
nfasis en el modelado de la gestin trmica.
El tema de la presente tesis se inscribe en el estudio de la gestin trmica de
los sistemas de transporte y, ms especficamente, en el estudio del balance
energtico experimental del motor durante su calentamiento, y en el estudio
terico-experimental de los flujos de calor en estado estacionario y transitorio
dentro del motor con ayuda de un modelo trmico. Del estudio del balance
energtico del motor durante su calentamiento se ha concluido la necesidad
de proponer estrategias de refrigeracin que, sin afectar mucho la topologa
del sistema de refrigeracin, produzcan un calentamiento ms rpido del
1. Introduccin
1.2. Objetivos
Esta tesis ha sido realizada con el fin de estudiar los fenmenos de
transmisin de calor y el balance trmico en motores de combustin interna
alternativos, tanto en el cilindro como en sus sistemas de refrigeracin, en
condiciones estacionarias y transitorias, mediante la utilizacin y adaptacin
de un programa computacional de prediccin de temperaturas de paredes del
motor y de dos programas comerciales de modelado de sistemas
termohidrulicos, tambin utilizados para modelar, simular y validar diseos
nuevos de sistemas de refrigeracin de motores para locomotoras. Para lograr
este objetivo general se han considerado los siguientes objetivos especficos:
1.3. Antecedentes
1.3. Antecedentes
El modelado trmico del motor, visto este como fuente principal de
disipacin de calor en los medios de transporte, comprende el tratamiento de
seis subsistemas: las pipas de admisin y escape, la cmara de combustin o
el cilindro, las camisas de enfriamiento en el bloque y la culata, el aceite de
enfriamiento (aprovechado para enfriar la parte inferior del pistn y que
desempea un papel importante en la friccin entre la camisa y los anillos,
friccin que genera un flujo de calor adicional entre el pistn y la camisa) y,
el sexto y ltimo subsistema, las partes slidas que conforman el motor y
sirven para separar los otros cinco subsistemas. Estas partes se someten a
transferencia de calor inestable tridimensional. La inestabilidad es de dos
formas: las fluctuaciones rpidas en las superficies cilindro-gas y los
transitorios algo menos rpidos causados por cambios en la carga y la
velocidad. El constante almacenamiento y entrega de energa a ambos niveles
de frecuencia causa una prdida potencialmente grande de la energa
disponible, porque la energa es absorbida a elevadas temperaturas y liberada
luego a bajas temperaturas. Otros subsistemas importantes que entran en
juego en la transferencia de calor son el turbocompresor, los interenfriadores
y el sistema de tratamiento de emisiones.
Muchos de los trabajos relevantes relativos a la transferencia de calor y la
gestin trmica en los motores de combustin interna han sido divulgados en
las investigaciones de Borman y Nishiwaki [1.17], Finol y Robinson [1.20], y
Pang y Brace [1.1]. Algunos de los trabajos tericos y experimentales
relacionados con el estudio y mejora de los sistemas de refrigeracin son los
realizados por Luptowski et al. [1.22], Cortona [1.23], Chalgren y Traczyk
[1.3], Allen [1.19], Mohan et al. [1.24], entre otros.
Los modelos para valorar las temperaturas del motor pueden ser de un solo
nodo (globales), de tres nodos (tres nodos representados por la camisa, el
pistn y la culata), de nodos mltiples, y modelos de elementos finitos. Los
tres primeros se cuentan entre los modelos de capacitancias concentradas. El
mtodo de capacitancias concentradas es uno de los mtodos ms utilizados
1. Introduccin
10
1.3. Antecedentes
11
1. Introduccin
12
13
1. Introduccin
14
1.5. Bibliografa
15
1.5. Bibliografa
[1.1]
[1.2]
[1.3]
1. Introduccin
[1.4]
[1.5]
[1.6]
[1.7]
[1.8]
[1.9]
[1.10]
[1.11]
[1.12]
[1.13]
[1.14]
[1.15]
16
1.5. Bibliografa
[1.16]
[1.17]
[1.18]
[1.19]
[1.20]
[1.21]
[1.22]
[1.23]
[1.24]
[1.25]
17
1. Introduccin
18
Captulo 2
Estado del Arte de la Transferencia
de Calor en los M.C.I.A.
2.1. Introduccin ............................................................................................ 21
2.2. La transferencia de calor en el motor ...................................................... 22
2.3. Transferencia de calor entre los gases del cilindro y las paredes de
la cmara de combustin ................................................................................ 23
2.3.1. Modelos globales.............................................................................. 26
2.3.2. Modelos zonales y locales................................................................ 32
2.3.3. Modelos multidimensionales............................................................ 46
2.3.4. Transferencia de calor radiativa ....................................................... 51
2.4. Transferencia de calor en las pipas de admisin y escape ...................... 62
2.5. Transferencia de calor del lado del refrigerante...................................... 67
2.5.1. Modelos empricos ........................................................................... 69
2.5.2. Modelado CFD................................................................................. 72
2.6. Transferencia de calor entre el aceite y el motor .................................... 73
2.6.1. Transferencia de calor entre el pistn y el aceite ............................. 73
2.6.2. Transferencia de calor entre el aceite y la estructura (bloque y
culata) ......................................................................................................... 77
2.7. Transferencia de calor entre el motor y el ambiente ............................... 78
2.8. Conclusiones ........................................................................................... 81
2.9. Bibliografa.............................................................................................. 83
20
2.1. Introduccin
21
2.1. Introduccin
La transferencia de calor es un problema de gran inters para la comunidad
involucrada en el estudio de los motores de combustin interna (M.C.I.) por
las implicaciones que tienen las prdidas de calor hacia las paredes sobre el
rendimiento termodinmico del proceso de conversin de energa, por la
relacin que tiene sobre la elevacin de la temperatura de la superficie y con
ello la influencia sobre el comportamiento tribolgico y la durabilidad de los
componentes de la cmara de combustin y en particular del grupo cilindropistn. Adems de esto, la cintica de las reacciones qumicas que controlan
la formacin de emisiones, la formacin del holln y la aparicin del golpeteo
dependen en gran medida de la temperatura de los productos en la cmara, la
cual generalmente vara ampliamente dentro del espacio y al pasar el tiempo
dentro de la cmara de combustin. No es reciente ya la controversia entorno
al desempeo de los motores adiabticos y son muy actuales los desarrollos
en la gestin trmica de los vehculos. Todo esto motiva a repasar el estado
en el que se encuentra la temtica de la transferencia de calor en los motores
modernos.
En 1987 Borman y Nishiwaki [2.1] publicaron una vasta revisin del estado
de la transferencia de calor en los motores de combustin interna. Finol y
Robinson [2.84] presentaron tambin una revisin, pero sta se limit a
enumerar las principales correlaciones encontradas en la literatura del
coeficiente global de transferencia de calor en el cilindro del motor. Puede
adelantarse que las conclusiones del trabajo de Borman y Nishiwaki no han
variado mucho y lo que ha venido marcando la diferencia es el incremento de
trabajos basados en los mtodos computacionales. Se ha trabajado ms sobre
los modelos detallados del coeficiente de pelcula local, el procesado de
mediciones de temperaturas instantneas de pared y el refinamiento de
algoritmos matemticos para el procesamiento de esta informacin.
El estudio del estado del arte de la transferencia de calor en los M.C.I. en esta
tesis se ha abordado con el objetivo de explorar los modelos y trabajos
experimentales actuales empleados en el diagnstico y la prediccin de las
temperaturas de las paredes del motor y los flujos de calor en su interior, as
como la interrelacin trmica con los sistemas trmicos perifricos
relacionados con l, incluyendo, claro est, al ms convencional de ellos, al
22
Transferencia de
calor dentro del
cilindro
Transferencia de
calor en otros
sistemas del motor
Estudios analticos
Gestin trmica de
los motores de
combustin
Estudios
experimentales y
numricos
23
Las temperaturas de la estructura y de los fluidos son los factores que limitan
el incremento de la potencia especfica del motor (potencia por unidad de
volumen de trabajo del motor) y es por esto que es importante, no slo
continuar garantizando los niveles seguros de temperaturas segn las
limitaciones de tensiones trmicas de los materiales estructurales, las
tolerancias en los pares cinemticos, las exigencias tribolgicas y los lmites
impuestos por las propiedades de los fluidos, sino tambin avanzar en la
identificacin trmica de los motores, esto es, en la obtencin de modelos que
permitan predecir y controlar las temperaturas de los fluidos de las fronteras y
de las superficies de inters. Es importante poder predecirlas y controlarlas
para poder optimizar la combustin y reducir con ello el consumo y las
emisiones de contaminantes en todos los regmenes de operacin del motor,
en primer lugar, y luego para potenciar la recuperacin de la energa
evacuada con los gases de escape y con el refrigerante.
Para identificar trmicamente el motor se debe disponer de modelos que
permitan validar los intercambios de calor entre los gases de admisin y
escape y las pipas correspondientes, entre los gases del cilindro y las paredes
de la cmara de combustin, y as mismo entre el refrigerante, el aceite
lubricante y el aire circundante y las correspondientes superficies del motor.
En las partes siguientes de este apartado se comentarn los tratamientos de
modelado ms usados para estas interfases.
24
25
Figura 2.3. Coeficiente de transferencia de calor entre los gases del cilindro
y las paredes de la cmara de combustin. Las pruebas fueron realizadas en
un motor a gasolina Rover K16 [2.12].
Para cualquier zona de la cmara de combustin, lo correcto al determinar el
flujo de calor a travs de cualquier superficie es resolver la ecuacin de la
energa que, en su forma adimensional simplificada y con los supuestos de
que la conductividad y el calor especfico de los gases son constantes y de
que la presin es espacialmente uniforme, puede presentarse constituida por
sus trminos de almacenamiento, conduccin, conveccin, trabajo de presin
y fuente de calor (por conveccin y radiacin), como se describe en la
ecuacin (2.1):
1 T * dp*
T *
1
=
2T * v*T * +
+ q& * ,
(2.1)
*
Re Pr
p d
donde p*= p/p0, T*= T/T0 y v* = v/v0 son la presin, la temperatura y las
velocidades adimensionales, escaladas por los valores de referencia p0, T0 y
26
27
Nusselt en
(1923)
Briling
(1935)
Van Tyen
(1962)
Eichelberg,
(1939)
Pflaum
(1961)
Mtodo
Expresin y comentarios
Anlisis experimentales.
Recipiente
de
alta
presin y experimentos
en un motor real.
Originalmente obtenida
para predecir el flujo de
calor promediado en el
tiempo o estacionario.
Mejora de la ecuacin de
Nusselt
Mejora de la ecuacin de
Nusselt
Consider la radiacin.
Primera medicin del
flujo instantneo de
calor.
Basndose en la ecuacin
de Eichelberg, incorpora
el
efecto
del
supercargador,
empleando
diferentes
coeficientes
para
la
culata, el pistn y el
cilindro.
h = hc + hr =
= 5, 41103 ( 1 + 1, 24cm
)3
( T /100 ) ( Tw /100
4, 21104
T Tw
( 1/ g + 1/ w 1 )
4
h = ( 3,19 + 0,885cm
h = 2,1 3 cm
)3
)3
Annand
(1963)
LeFeuvre
(1969)
Annand y
Ma (1970)
Incorpora en la expresin
para el flujo instantneo
de calor el efecto del
swirl. Utiliza un nmero
de Reynolds local (por lo
que
estrictamente
hablando no es un
modelo global).
Extienden la frmula
inicial
de
Annand,
introduciendo
una
correccin
para
la
transitoriedad. Emplea
los resultados de un
motor
de
inyeccin
directa con swirl en la
admisin. Adopt la
p 2T
p 2T .
pT
h = K D , K , B pT f (cm ) ,
aplicando
la culata y el pistn;
, kW/m2K
p 2T +
q=a
k
(Re )0,7 (T Tw ) + b T 4 Tw4
D
),
0,35
0,8
q=a
k (t )
Re(t ) 0,8 Pr(t ) 0,33 {T (t ) T (t )} ,
r
r es un radio en el
Re =
r 2
(swirl); a = 0,047.
q=a
k
(Re )0,7 a(T Tw ) + a dT + b T 4 Tw4 , es la
D
dt
Dent
y
Suliaman
(1977)
Woschni
(1967)
Sitkei
(1962)
Zapf
(1966)
Nishiwaki
(1979)
Hohenberg
(1979)
velocidad
media
propuesta por Knight.
Mejora a la expresin de
LeFeuvre. Aplicaron un
anlisis de transferencia
de
calor
convectiva
turbulenta para flujo a lo
largo de una placa. Se
utiliza un nmero de
Reynolds local.
Trabaj sobre el balance
de energa , para
determinar el calor total
transferido a las paredes
de la cmara en cada
proceso
del
ciclo
completo de trabajo del
motor
Adopt la relacin para
la transferencia de calor
convectiva turbulenta de
Annand.
Aport un efecto del
swirl durante el proceso
de intercambio de gases
Nmero de Nusselt en las
carreras de admisin y
escape
Modific la ecuacin de
Woschni, para ajustar
mejor las predicciones a
los flujos instantneos de
calor medidos.
28
k r 2
q = 0,023
r
0,8
(T Tw )
VT
W = C1C m + C 2 s 1 p p0
p1V1
h = 2,36 104 (1 + b )
p 0,7Cm0,7
T 0,2 d e0,3
kW/m2K,
p en bar.
Oguri
(1960)
Mejora a la expresin
de Elser
S
(Re Pr )0,5
Nu = 6,51 + 0,5
cp
S
(Re Pr )0,5 2 + cos 20 o
Nu = 1,751 + 0,5
cp
)]
* En las expresiones de la tabla h es el coeficiente de transferencia de calor (kW/m2K); g, w son las emisividades
del gas y de la pared, respectivamente; cm es la velocidad del pistn, en m/s y p es la presin del gas en MPa; T y Tw
son las temperaturas del gas y de la pared, respectivamente; siendo k la conductividad trmica del gas; D el dimetro
del cilindro; en la frmula de Annand y Ma a, a y b son constantes obtenidas para el motor en cuestin y dependen
de las revoluciones del motor y de la relacin aire combustible. La derivada en el tiempo de la temperatura media
global se agrega en la frmula para compensar la naturaleza no estacionaria del flujo. El nmero de Reynolds se
calcula con la velocidad obtenida a partir del valor medio de la energa cintica del gas.
En la ecuacin de Woschni [2.16] W es la velocidad de referencia, la cual representa la velocidad media del gas que
afecta la transferencia de calor y est dada para cada proceso. Cm es la velocidad media del pistn. Para el proceso de
intercambio de gases C1 = 6,18, C2 = 0. Para el proceso de compresin C1 = 2,28, C2 = 0. Para los procesos de
combustin y expansin C1 = 2,28, C2 = 3,2410-3. p0 es la presin en MPa obtenida en arrastre, Vs es el volumen
barrido por el pistn, en m3. T1, V1 y p1 son determinados para un tiempo determinado en el que se conocen la
presin y la temperatura, por ejemplo, el cierre de la vlvula de admisin o el inicio del proceso de combustin.
Aunque fue escrita en la forma de transferencia de calor convectivo turbulento estacionario, la frmula de Woschni
29
incluye de manera concentrada todos los efectos de la conveccin y la radiacin. Posteriormente Woschni agreg un
trmino de swirl a su ecuacin, cambiando la forma del coeficiente C1 [2.115]. Para el proceso de intercambio de
gases C1 = 6,18 + 0,417Cu/Cm; para el resto del ciclo C1 = 2,28 + 0,308/Cm; Cu =DnD, y nD es la velocidad de
rotacin del rodete usado en el banco de pruebas de swirl.
En la expresin de Sitkei V es el volumen del cilindro, A es el rea que absorbe el calor, la constante adimensional b
se agreg en consideracin a velocidades de turbulencia adicionales y est entre 0 y 0,35, dependiendo del tipo de
cmara de combustin.
30
Visto el desarrollo que tienen los actuales sistemas de control de los motores,
hoy es posible utilizar algoritmos de control menos dependientes de tablas de
caracterizacin del motor experimentales y aprovechar la capacidad
computacional de los microcontroladores empleados en ellos para
implementar algoritmos basados en funciones. Para estos nuevos sistemas de
control es ya, entonces necesaria la utilizacin de modelos que todava
pueden corresponder a modelos globales para todos los procesos que ocurren
en el interior del cilindro. Un ejemplo de este tipo de aplicacin es el
tratamiento propuesto por Allmendinger et al. [2.23], quienes desarrollaron
un modelo predictivo para ser implementado en un programa para el control
en lazo cerrado, independiente de mapas, de un motor Diesel de inyeccin de
colector comn. Basndose en la relacin lineal aparente entre el coeficiente
de transferencia de calor y la temperatura de los gases, los autores
concibieron una expresin con coeficientes gausianos de la forma:
2
( E )
1
2
h ( , T ) = a ( ) + b ( ) T ; a , b ( ) =
e 2
(2.3)
2
2
31
32
33
34
Mtodo
Expresin y comentarios
1
2
Morel y
Keribar
(1985)
Basaron su trabajo en la
analoga
de
Colburn,
acoplada a un modelo de
dos zonas.
C f U eff C p P
2/3
1/ 2
2
2
;
U eff = U z + U y + 2 k
Ueff
es la
Bargende
(1990)
1/ 2
U eff
Vp (t )
= 2 K + 2 k +
, donde K y k corresponden a la
h=C
0,22
D
{
Char .length
Materialproperties
0,78
w
{
Gasvelocity
{
Combustion
Las constantes en la
ecuacin
fueron
determinadas para motores
a gasolina de inyeccin
directa.
35
0,22
1, 11V
0,073
2
+ Cp
;w =
dk
dt
2 k dV
3 V dt
1,5
L
1,5
ks
+ s
( V ) ; k
Un parmetro alternativo al
coeficiente de transferencia
de calor es la constante de
conductividad equivalente
de la capa lmite de espesor
Basada en la analoga de
Colburn. Emplean una ley
de potencia para definir el
perfil de velocidad:
St =
Cf
2
p c pU
( )
u + = C y + , igual que
Fiveland y
Assanis
(2002 )
h=
T T
y
= D
T Tw
*
*
0, 0153 F C pU
* 0,25
Aplicaron su modelo a un
motor HCCI. El modelo
tiene en cuenta la zona
gruesa (central) del gas o
zona de ncleo gaseoso
adiabtico, las regiones de
capa lmite y las regiones
intersticiales del grupo de
anillos.
, k eff es la
Chen y
Veshagh
(1989)
Tg Tw
0,25
T
Tw
0,4
El coeficiente de transferencia
h=
Cp
eff
Preff
eff
Preff
Pr
T
Pr
T
36
Modelos de transferencia de calor locales. Modelos de capa lmite con empleo de leyes
de pared.
Pinchon
[2.101],
Kleeman
[2.51],
referencian
a Launder
y Spalding
(1974).
q
h=
T T
w
; T
q
C p y
h =
Prt
u C T T
p w
(1 / k ) ln ( Ey + ) + P y
Prl
P = 9, 0
Prt
Prl
1
Prt
0,25
es
Han y
Reitz
(1997)
Modificaron la ley de la
pared. La turbulencia en las
regiones de nmeros de
Reynolds
elevados
la
modelan empleando un
modelo RNG k-. Se
consideran las variaciones
de
las
propiedades
termofsicas del fluido
dentro de la capa lmite.
qw =
C pu *T ln (T Tw ) 2,1 y + + 33,4 G / u *
( )
2,1ln y + + 2.5
k 0 u *
1
+ 0,0821 exp
w 0 ( )
p
kT0
p0
( ) dus (0. )
3
37
38
39
mientras que su espesor mnimo desciende a cerca de 0,3 mm cerca del punto
muerto superior. La variacin de la densidad en el espesor de la capa lmite
hace que la masa atrapada en la capa lmite en este punto corresponda a
aproximadamente el 30-35% de la masa total del cilindro.
Figura 2.10.
Variacin del
espesor de la capa
lmite trmica a lo
largo del ciclo de un
motor HCCI sobre
las superficies de
algunos nodos
seleccionados en la
camisa, el pistn y
la culata [2.13].
40
Figura 2.11. El mallado grueso no puede resolver la capa interna, por lo que
se requiere una relacin auxiliar que sirva de condicin de contorno, la
funcin de pared, para acoplar los flujos interno y externo.
El mtodo ms sencillo para modelar el flujo y los fenmenos asociados que
ocurren cerca de la pared ha sido el empleo de la ley de la pared (tratamiento
de funcin de pared), mediante la cual se suele expresar la velocidad media
paralela a la pared y las propiedades de turbulencia fuera de la subcapa
viscosa, en trminos de la distancia a la pared y de las condiciones de la pared
como son el esfuerzo cortante, el gradiente de presin, y la transferencia de
calor a la pared. As, las funciones de pared pueden usarse para proveer las
condiciones de contorno en las inmediaciones de la pared para las ecuaciones
de transporte de momento y turbulencia, en lugar de las condiciones en la
pared misma, de manera que no hay necesidad de resolver la subcapa viscosa
y se supera as la exigencia de un mallado muy fino. Este mtodo fue
propuesto por Launder y Spalding [2.91], quienes obtuvieron los perfiles de
velocidad y temperatura para las vecindades de la pared a partir de un anlisis
de flujo de Couette, asumiendo un flujo unidimensional estacionario e
41
42
u+ = y+
+
T = Pr u
, para y + 11, 63
(2.4)
En la subcapa inercial:
1
ln ( Ey + )
k
Pr
T + = Prt u + + P
Prt
u+ =
(2.5)
T+ =
+
y =
u C p (T Tw )
qw
yu
(2.6)
C0,25 k 0,5 y
43
44
0
t
p0 z 2
p0 z Prt z
14243 144424443
la min ar
donde 0 = k0 / C p .
(2.7)
turbulento
45
46
47
con sus partes mviles (pistn y vlvulas); y los algoritmos de solucin, cuya
funcin es resolver las ecuaciones algebraicas.
La mayora de los programas CFD comerciales emplean o una funcin de
pared o una modificacin de sta para describir el flujo en las cercanas de las
paredes en los clculos de transferencia de calor. Las ecuaciones promediadas
de Navier-Stokes para flujo compresible tridimensional se resuelven en un
mallado mvil, mientras que los flujos turbulentos son modelados bajo un
concepto de viscosidad turbulenta, empleando el modelo -.
Los programas CFD comerciales pueden acoplarse de manera iterativa a
programas de clculo de temperaturas de pared basados en elementos finitos,
FE (Fluent se puede acoplar con ANSYS, por ejemplo) y permitir el
clculo tambin de esfuerzos y deformaciones trmicas. En la primera parte
de la metodologa, el ciclo del motor se calcula con temperaturas constantes
asignadas para las diferentes partes de la cmara de combustin. Luego, en la
ejecucin del segundo ciclo, el programa CFD utiliza las temperaturas
calculadas por el programa FE, que ya son ms cercanas a las que se
asignaron en el primer ciclo [2.47].
Figura 2.15. Ubicacin de los sensores de flujo de calor en el pistn del motor
HPP [2.51].
48
Figura 2.16. Flujos de calor predichos en el pistn del motor HPP con el
modelo CFD, empleando una funcin de pared modificada [2.51].
49
50
51
52
Morel [2.7] formul un modelo de longitud de haz promedio simple para ser
usado en un programa de anlisis de ciclo. El modelo desprecia la radiacin
del gas, la radiacin pared-pared y la dispersin. La radiacin de material
particulado es tomada como gris. El flujo hacia una superficie con emisividad
i y temperatura Ti es entonces:
qi = i i si (Tr4 Ti 4 ) ,
(2.8)
donde Tr es la temperatura de radiacin del holln, si = 1-exp[1575(fm/s)TrLi] (2.10) y fm es la fraccin de la masa de holln, s/ es la
relacin entre las densidades del holln y el gas, Li es el recorrido radiativo
medio para la isima superficie. El recorrido radiativo medio se toma de
alguna manera arbitraria, relacionada estrictamente con la geometra de la
cmara, por ejemplo, para una superficie del cilindro se puede tomar el
dimetro de ste. La temperatura de radiacin fue definida arbitrariamente en
trminos de la temperatura de la zona quemada. La correlacin para fm se
obtuvo de los resultados empricos de Sitkei [2.22]. Borman y Nishiwaki
[2.1] estiman que es un modelo de trabajo razonable si se cuenta con
informacin para ajustar las frmulas de Tr y fm.
El modelo matemtico que describe el campo de radiacin en el cilindro del
motor es la ecuacin de transferencia radiativa (RTE) [2.98],[2.93]:
( ) I ( r , ) = ( anet + s ) I ( r , ) + kI b ( r ) +
s
S ( r, )
4
(2.9)
(2.10)
53
54
nmero de stos que han sido absorbidos por cada uno de los elementos
discretizados de la superficie. El flujo de calor radiativo QAj para la
superficie elemental Aj de la cmara de combustin puede definirse como:
kai
r
exp 0 kai rdr Tgi4 dVi dAj +
3
Aj V r
cosi cos j
r
+
exp 0 kai rdr siTsi4 dAi dAj
3
r
Aj V
QAj =
(2.11)
Ms
i =1
i =1
(2.12)
55
La Figura 2.20 muestra los flujos de calor radiativos en la parte superior del
pistn y en centro de la culata en un motor Diesel de inyeccin directa, junto
con los valores promedio en todas las superficies de la cmara de combustin
predichos con el mtodo multidimensional de Chen et al. [2.59], para cuya
parte radiativa se utiliz el mtodo de Monte Carlo.
56
57
58
n I ( n ) d
n<0
(2.13)
59
60
Tabla 2.3. Prdidas de calor para diferentes relaciones de swirl, (J) [2.48].
Combinada
Pistn
Camisa
Culata
Relacin de swirl
0,0
1,4
2,2
% de variacin
0,0 a 1,4
1,4 a 2,2
Convectiva
22,722
24,796 26,227
-3,60
5,77
Radiativa
7,221
7,147
10,119
-1,02
41,58
Total
32,942
31,943 36,346
-3,03
13,78
Convectiva
21,970
10,660
8,443
-51,48
-20,80
Radiativa
13,570
9,714
7,532
-28,42
-22,46
Total
35,540
20,374 15,975
-42,67
-21,59
Convectiva
38,751
34,477 29,685
-11,35
-13,90
Radiativa
11,89
10,490
6,543
-11,77
-37,626
Total
50,650
44,967 36,229
-11,22
-19,43
-19,10
-16,33
-18,34
-7,97
-11,54
-8,98
Convectiva
Radiativa
Total
61
Figura 2.25.
Temperatura pico del
pistn en funcin de la
fase de dispersin ()
[2.48].
62
63
(2.14)
Nu = 0, 023Re0,8 Pr n
(2.15)
(2.16)
(2.17)
Bauer et al.:
Nu = 0, 062Re
Nu = 0,140Re
0,73
0,66
Shayler et al.:
Nu = 0,135Re
0,713
[2.64]
(2.18)
Meisner y Sorenson:
Nu = 0, 0774Re
0,769
(2.19)
Wendland:
Nu = 0, 027Re
0,8
Pr
0,33
0,14
64
En esta ltima, el termino Centr tiene en cuenta los efectos de separacin del
fluido de las paredes al entrar al puerto. y s corresponden a la viscosidad
evaluada a las temperatura media del flujo y de la superficie,
respectivamente.
El nmero de Nusselt en el clculo del coeficiente de pelcula en las
pipas de admisin y escape emplea velocidades instantneas en lugar de
velocidades medias, para tener en cuenta las caractersticas pulsantes del
fluido, particularmente en el estudio del escape. Huber [2.107],
basndose en medidas de velocidad y temperatura realizadas en
conductos de escape, observ desfases temporales entre la seal de
velocidad y temperatura, por lo que propuso correlacionar el nmero de
Nusselt con una velocidad promedio mediante la ecuacin Nu = 1, 6 Re0,4 ,
en la que se pondera la historia de velocidades a su paso por el ducto.
Reyes [2.105] modific la expresin de Huber para dar cuenta la
dependencia que tiene el coeficiente de pelcula en el colector de escape
no slo de la velocidad instantnea del fluido en el instante considerado,
sino tambin de la historia de esa velocidad en instantes anteriores. En el
modelo de Reyes, la velocidad caracterstica se calcula como la suma de
la velocidad presente y las velocidades previas multiplicadas por sus
respectivos coeficientes de disipacin:
V (t ) =
k
c V ( t k t )
k =0
(2.21)
k =0
(2.22)
65
66
3/ 4
[2.67]
(2.23)
Nu = 0, 023Re Pr
0,8
0,4
D dU
1 0, 75 2
para la primera etapa,
U dt
(2.24)
0,8
D dU
0,8
0,4
018
Re
Pr
1
0,
75
Nu = 0,
correccin
no _ estable
U dt critical
(2.26)
67
68
69
0,49
k 0,79
c 0,45
f
p, f f
0,5 0,29
( h fg g )
f
F = (Twall Tsat )
0,24
(P
sat , wall
0,24
F S
Pcool )
0,75
(2.28)
70
S = 3, 4 0, 6 log ( Re )
Nu = 0, 027 Re
0,8
D
Pr bulk
wall
1
0,14
(2.30)
siendo bulk y wall la viscosidad del fluido a las temperaturas media del flujo y
de pared, respectivamente. Otros autores usan correlaciones similares con
correcciones especficas segn el motor.
Finlay et al. [2.72] demostraron que tanto la ebullicin nucleada como la
ebullicin de pelcula se presentan con frecuencia en las camisas de
refrigeracin de los motores. Norris et al. [2.69] realizaron mediciones de
temperaturas y un anlisis de transferencia con elementos finitos para
encontrar las regiones de conveccin pura, ebullicin nucleada y ebullicin
71
0,4
* 0, 091
Dh
1 + 23,99 Re 0,23 ( x / Dh )
19
1 ( x0 x )9 10
0,125
Re
0,363( Dh )
0,10
bulk
wall
0,14
(2.31)
72
73
74
Figura 2.32. Influencia del tipo de refrigeracin del pistn sobre las
direcciones del flujo de calor hacia la camisa y el aceite [2.112].
En dicho trabajo, se estim necesario introducir una correccin por velocidad
del motor:
hact
n
= href
n
ref
(2.32)
75
k
hpc = Apc Re0,8 Pr 0,33 Di0,33 0
D
(2.33)
4600
Para la corona
huc 2
Para la falda
0,35
N
= 240
4600
(2.34)
0,35
(2.35)
Nu
= 1 + a eb ( r / d )
Nu0
1 9
, Nu0 = 2, 67 Re
0,567
z
Pr 0
d
0.4
0,336
0,237
(2.36)
76
D
hoil = 68,17 ( rw ) n
b
(2.37)
d gal oil
Buscando excluir la dependencia en la expresin de las propiedades
termofsicas del aceite, durante el proceso de validacin, se busc un
'
m Con estas consideraciones los
coeficiente de la forma Cgal
= C gal koil / oil
valores encontrados fueron Cgal = 228,5 y m= 0,73. Los autores compararon
77
Para las partes de la culata en contacto con el aceite, Jarrier y Gentile [2.54]
utilizan una expresin Nu-Re para superficies planas en la forma:
Nu = 0, 664 Re0,5 Pr 0,33
(2.39)
Para el aceite que se escurre por las paredes del cilindro hacia el bloque,
Jarrier y Gentile utilizan una expresin establecida por Hausen:
Nu = 3, 66
0, 00668
1 + 0, 04
Dhyd
2L
Dhyd
2L
Re Pr
(2.40)
Re Pr
N 90
= hREF
2000 T
0.5
(2.41)
78
79
(2.42)
(2.43)
80
Nu = 1,86 Re Pr
basada en Sieder y
L s
Tate
Nu = 0,116 Re 125 Pr b
s
23
13
0,14
D1 2 3
1 +
L
13
13
0,14
Re > 10000
basada en Sieder y
Tate.
2.8. Conclusiones
81
Figura 2.34.
Temperatura en el
compartimento del motor
en el plano de simetra
(configuracin
condensador-radiadorventilador, RFM)
[2.114].
2.8. Conclusiones
En el presente captulo se ha realizado un anlisis del estado del arte en el
modelado de la transferencia de calor en los motores de combustin interna
con especial atencin a las expresiones empleadas para calcular los flujos
trmicos a travs de las principales interfases entre los elementos.
82
2.9. Bibliografa
83
En una segunda parte del captulo, se repasaron los modelos usados para
calcular el calor transmitido entre los gases de admisin y escape y sus
respetivas pipas. La mayora de las correlaciones encontradas valoran la
transferencia de calor a partir de la relacin Nu-Re-Pr, con ajuste para la
viscosidad del gas. Se destaca el trabajo de Depcik y Assanis [2.67], quienes
presentaron una correlacin universal y luego tambin presentaron una
propuesta para evaluar la transferencia de calor transitoria.
En una tercera parte de este captulo se ha explicado la complejidad del
proceso de transferencia de calor hacia el refrigerante y se han descrito
algunos de los modelos encontrados en la literatura, la mayora de ellos
basados en la ecuacin de Dittus-Boelter. Los trabajos ms importantes son
los trabajos de Norris [2.69] y de Robinson [2.70]. Los avances en los
modelos computacionales se estn aplicando a los flujos en la camisa de
refrigeracin y se espera que as como los modelos globales y zonales dentro
del cilindro han sido mejorados con los resultados de los modelos CFD, as
mismo sean mejorados los modelos para la interfaz paredes motor
refrigerante.
En el resto del captulo, y de una manera menos detallada se ha comentado
algunos de los tratamientos dados al clculo de los flujos de calor en las otras
interfases del motor, como entre el aceite y la estructura, entre el pistn y el
cilindro y entre la estructura del motor y el aire ambiente.
Los modelos de transferencia de calor se han estudiado dentro de la presente
tesis porque constituyen la base del modelado de las temperaturas y los flujos
de calor del motor. Se entiende as, que la adecuada integracin de las
expresiones para el clculo de los flujos de calor entre los componentes del
motor, en su parte interna y desde su superficie es la base de un modelo de
balance energtico del motor (modelo interno) que puede acoplarse al modelo
de gestin trmica del motor (modelo externo).
2.9. Bibliografa
[2.1] Borman, G., Nishiwaki, K. Internal combustin engine heat transfer.
Progress in Energy Combustion Sciences, 1987. Vol 13, pginas 1-46.
[2.2] Rakopoulos, C. D., Mavropoulos, G. C. Experimental instantaneous
[2.3]
[2.4]
[2.5]
[2.6]
[2.7]
[2.8]
[2.9]
[2.10]
[2.11]
[2.12]
[2.13]
84
2.9 Bibliografa
[2.14]
[2.15]
[2.16]
[2.17]
[2.18]
[2.19]
[2.20]
[2.21]
[2.22]
[2.23]
[2.24]
[2.25]
85
2001-01-1028.
Buttsworth, D. R. Transient response of an erodable heat flux gauge
using finite element analysis. Proc. Instn. Mech. Engrs, part D,
Journal of Automobile engineering, vol. 216, pp. 701-706, 2002.
R. Tatschl, B. Basara, J. Schneider, K. Hanjalic, M. Popovac, A.
Brohmer, J. Mehring. Advanced turbulent heat transfer modeling for
IC-engine applications using AVL FIRE. International
Multidimensional Engine Modeling Users Group Meeting. April 2,
2006, Detroit, MI.
Woschni, G. A Universally Applicable Equation for the
Instantaneous Heat Transfer Coefficient in the I.C. Engine. SAE
paper 670931.
Hohenberg, G. F. Advanced approaches for heat transfer
calculations, Diesel engine thermal loading. SAE SP-449, pp. 61-79,
1979.
Hiroyasu, H. Diesel engine combustion and its modeling.
International Symposium, COMODIA 85, pp. 53 -76, 1985.
Wilson, T.S., Bryanston-Cross, P. J., Chana, K. S. Dunkley, P.,
Jones, T. V., Hannah, P. High bandwidth heat transfer and optical
measurements in an instrumented spark ignition internal combustion
engine. SAE paper 2002-01-0747.
Heywood, J.B. Internal combustion engine fundamentals. McGrawHill, 1988.
Stone, R. Introduction to internal combustion engines. MacMillan
publishers, 1999.
Sitkei, G. Heat transfer and thermal loading in internal combustion
engines. Akademiai Kiado, 1973.
Allmendinger, K. Guzzella, L., Seiler, A., Loffeld, O. A method to
reduce the calculation time for an internal combustion engine model.
SAE paper 2001-01-0574.
Kavtaradze, A. I., Gaivoronskii, V. A. Fedorov, D. D. Onishchenko,
A. V. Shibanov. Calculation of radiative-convective heat transfer in
the combustor of Diesel engine. High Temperature, vol. 45, n 5, pp.
673680, 2007.
Schubert, C., Wimmer, A., Chmela, F. Advanced heat transfer model
[2.26]
[2.27]
[2.28]
[2.29]
[2.30]
[2.31]
[2.32]
[2.33]
[2.34]
[2.35]
[2.36]
[2.37]
86
2.9 Bibliografa
[2.38]
[2.39]
[2.40]
[2.41]
[2.42]
[2.43]
[2.44]
[2.45]
[2.46]
[2.47]
[2.48]
87
[2.49]
[2.50]
[2.51]
[2.52]
[2.53]
[2.54]
[2.55]
[2.56]
[2.57]
[2.58]
[2.59]
88
2.9 Bibliografa
[2.60]
[2.61]
[2.62]
[2.63]
[2.64]
[2.65]
[2.66]
[2.67]
[2.68]
[2.69]
[2.70]
89
[2.71]
[2.72]
[2.73]
[2.74]
[2.75]
[2.76]
[2.77]
[2.78]
[2.79]
[2.80]
[2.81]
[2.82]
90
2.9 Bibliografa
[2.83]
[2.84]
[2.85]
[2.86]
[2.87]
[2.88]
[2.89]
[2.90]
[2.91]
[2.92]
[2.93]
[2.94]
[2.95]
91
[2.96]
[2.97]
[2.98]
[2.99]
[2.100]
[2.101]
[2.102]
[2.103]
[2.104]
[2.105]
[2.106]
92
2.9 Bibliografa
[2.107]
[2.108]
[2.109]
[2.110]
[2.111]
[2.112]
[2.113]
[2.114]
[2.115]
93
94
Captulo 3
Caracterizacin trmica de los
M.C.I.A.
3.1. Introduccin ............................................................................................ 97
3.2. Proceso de combustin y flujo de calor a travs de las paredes de la
cmara de combustin de un motor alternativo ............................................. 99
3.2.1. El anlisis de Fourier...................................................................... 100
3.2.2. Transformada de Laplace ............................................................... 102
3.3. Balance trmico del motor .................................................................... 103
3.3.1. Balance energtico del motor en estado estacionario a carga
constante................................................................................................... 105
3.3.2. Balance energtico durante el calentamiento del motor................. 107
3.3.2.1. Consideraciones tericas y experimentales......................... 107
3.3.3. Balance energtico experimental ................................................... 117
3.3.3.1. Balance trmico del motor durante su calentamiento
bajo una carga y rgimen constantes, representativos del
ciclo de conduccin ............................................................... 118
3.3.3.2. Balance trmico del motor durante su calentamiento en
su proceso de calentamiento en el ciclo de conduccin
NEDC .................................................................................... 123
3.3.3.3. Energa recuperada en el vehculo durante el ciclo de
conduccin............................................................................. 130
96
3.1. Introduccin
97
3.1. Introduccin
El objetivo de este captulo es estudiar a la caracterizacin trmica de los
motores, entendiendo por caracterizacin la determinacin de las
temperaturas del motor, sus flujos de calor internos y externos y la inercia
trmica de sus paredes internas. Esta caracterizacin es particularmente
necesaria para alimentar programas de modelado de la combustin, de
operacin de sistemas de refrigeracin y tambin para derivar conclusiones
relativas a la utilizacin de la energa suministrada al motor. El trabajo
presentado en este captulo es el resultado de la combinacin de las
evaluaciones experimentales realizadas en un motor Diesel y la aplicacin y
extensin de un modelo trmico a la prediccin de temperaturas y flujos de
calor de este mismo motor.
Para lograr el objetivo mencionado se ha dispuesto la realizacin de ensayos
y adquisicin de informacin en el motor en las siguientes fases:
98
100
T
1 T
=
k
t c x x
(3.1)
Tw = Tm + [ An cos ( nt ) + Bn sin ( nt )]
(3.2)
n =1
(3.3)
N
k
T
T
k
+
( m l ) n ( An + Bn ) cos nt ( An Bn ) sin nt (3.4)
l
n =1
102
= hexp (T Tw,i ) Ai ,
dt exp
i
donde hexp =
(3.5)
qw
T Tw
( , s ) = s ( s )
103
e p e p( 2)
1 e 2 p
(3.6)
Ts
1
*)
* (
x 0 d
*
k
n 2 / ( * )
*
+
1
2
e
(3.7)
104
modelo detallado para calcular los transitorios trmicos del motor, pero como
puede inferirse de su trabajo, el modelo fue validado empleando un proceso
transitorio trmico controlado (la temperatura del motor se vari
manualmente entre 300 K y 380 K, realizando las mediciones cuando la
temperatura se estabilizaba). El trabajo de Morel no estaba dirigido a obtener
un balance energtico. Jarrier et al. [3.18] encerraron un motor en un
calormetro para realizar el balance trmico instantneo para diferentes
condiciones de operacin, identificando las rutas principales seguidas por los
flujos de energa. Aunque en el artculo se sugiere que la parte experimental
se ha realizado para simular ciclos de conduccin urbanos (ECE), no se
presentan en el trabajo resultados de balance energtico bajo el ciclo NEDC
completo. Jarrier et al. asumieron que durante el calentamiento del motor las
presiones instantneas de los gases en el cilindro eran idnticas a las del
estado estacionario, no considerando el hecho de que durante el
calentamiento la presin media indicada es afectada por la presin media
efectiva de friccin, la cual depende de la viscosidad del aceite (esto es, de la
temperatura). Shayler et al. [3.19] desarrollaron un modelo computacional
para predecir las caractersticas de calentamiento de un motor utilizando un
modelo de masas globales junto con correlaciones empricas para el clculo
disipado al refrigerante por el motor (correlacin de Taylor y Toong) y las
prdidas de friccin (correlacin de Patton). Los objetivos de otros trabajos
relacionados (Taymaz [3.20], Perez-Blanco [3.21], Lehner [3.22]) son
diferentes a la realizacin de un balance energtico. Adems en los trabajos
encontrados se hace uso de informacin obtenida bajo condiciones
estabilizadas.
En este apartado se hace mencin del balance energtico del motor de manera
genrica, a partir de la informacin disponible en la literatura, y luego, para el
motor estudiado, se calcula el balance energtico bajo condiciones de
calentamiento a carga y revoluciones constantes y tambin bajo el ciclo de
operacin NEDC. Tambin se evala la energa que realmente se recupera en
los dispositivos suplementarios del sistema de refrigeracin del motor,
representados por el enfriador de gases de escape recirculados, EGR, y el
calentador del habitculo. La informacin relativa a los gradientes de
temperatura y flujos de calor en condiciones de operacin transitorias puede
usarse para prever las tecnologas de recuperacin de energa que mejor se
ajustan a esos niveles de gradientes de temperaturas y flujos de calor.
105
106
140
120
100
R 2 = 0,87
80
60
40
20
0
0
50
100
150
200
107
108
109
dm fug
dEsis dU dQcomb dQrefr dW dmadm
dm
=
=
+
+
+
hadm + esc hesc +
h fug (3.8)
dt
dt
dt
dt
dt
dt
dt
dt
La variacin de la energa interna del sistema en el motor de combustin
interna es la resultante de la variacin por unidad de tiempo de la energa
suministrada al sistema con el combustible (Qcomb) y la carga de admisin
(madmhadm), el calor realizado en forma de trabajo mecnico incluyendo las
prdidas mecnicas (W), el calor disipado a travs de las paredes al
refrigerante (Qrefr), el calor equivalente de los gases de escape referidos a las
condiciones de admisin (meschesc) y el calor equivalente de los gases de la
cmara fugados a travs de los intersticios en los segmentos (mfughfug).
Entendido el motor como mquina de realizacin de trabajo, la Figura 3.7 da
una cabal representacin del sistema para ser modelado desde el punto de
vista termodinmico. No obstante esto, el motor no puede funcionar con su
solo volumen de trabajo y debe interactuar con el medio externo para renovar
110
111
mientras que en el esquema de la Figura 3.8 las prdidas por friccin s estn
presentes e incluidas en la energa del refrigerante.
Figura 3.8. Esquema para la descripcin del balance energtico del motor.
Las sustancias de entrada y salida son transportadoras de energa.
Basndonos en el esquema de la Figura 3.8, el balance energtico del motor
se realiza en este trabajo, de manera prctica, midiendo los flujos msicos de
aire y combustible, el refrigerante, los gases de escape y tambin midiendo
las revoluciones y el momento torsor del motor, con los que se calcula la
potencia efectiva.
Para efectos prcticos de clculo del balance energtico del motor, la
ecuacin (3.8) puede reescribirse como:
m& comb H i _ combcomb + m& air hair = Q& equiv + m& air hair = N e + Q& refr + Q& misc + Q& esc
(3.9)
donde, m& comb y m& air corresponden a los flujos msicos de combustible y aire,
es la energa desprendida por la
respectivamente, h denota la entalpa, Q&
equiv
112
(3.11)
Siendo m& refr y c p _ refr el flujo msico y el calor especfico del refrigerante, y
ent
sal
Trefr
y Trefr
las temperaturas del refrigerante a la entrada y la salida del motor,
Q& esc = (m& comb + m& air ) hesc (Tesc ) hadm (Tadm )
(3.12)
113
estructura, mientras que otra parte se transporta con los fluidos, refrigerante y
aceite, progresivamente segn la expresin:
dQdis
dT
= mbl c pbl bl + m& refr c p _ refr dTrefr + m& ac c p _ ac dTac
dt
dt
(3.13)
donde Qdis representa el calor total disipado por el motor a travs de las
paredes; Trefr y Tac son las temperaturas del refrigerante y del aceite, mbl y cpbl
representan la masa trmicamente sensible del bloque y el calor especfico del
material, m& refr y m& ac son los flujos msicos de refrigerante y aceite a travs de
la camisa de refrigeracin y las galeras de lubricacin, respectivamente (el
flujo del aceite es ms complejo y la inclusin del trmino en la expresin es
una aproximacin grosera que no contempla el volumen de aceite en el
crter), y c p _ ac representa el calor especfico del aceite.
Descripcin de la instalacin experimental
Relacin de Compresin
18
Potencia mxima
Revoluciones mximas
4750 rpm
114
------Valor
Disco seco
Relacin de transmisin
3,417
1,841
1,167
0,8238
0,6730
3,690
1700 kg
2,795 m
0,0140
0,410 kgm2
1,033 kgm2
1,033 kgm2
0,89
0,321 m
0,300
185 N
2,223 m2
1,202 kg/m3
115
116
117
118
Balance energtico
70
60
qe
qesc
qrefr
50
40
30
20
10
0
4
15
23
35
a)
b)
119
120
t h A
Tmt ( t ) = T0 + (Test T0 ) 1 exp g t
mc
t t
(3.14)
121
122
misc
123
100
80
18.58
12.66
20
32.83
35.93
40
100.00
60
Qesc
Qrefr
0
Qequiv
Qmisc
Ne
124
125
126
127
60
30
-30
0
400
800
1200
Tiempo (s)
128
Porcentaje de la energa
suministrada (%)
70.0
60.0
50.0
40.0
30.0
20.0
10.0
0.0
ECE-1
ECE-2
ECE-3
ECE-4
MVEG
Ciclo
total
qe
19.5
21.0
21.0
20.3
30.5
22.4
qesc
10.2
12.7
13.1
13.6
16.1
13.2
qrefr
7.4
10.0
10.0
11.6
22.9
12.4
63.0
56.3
55.9
54.5
30.5
52.0
qmisc
Figura 3.20. Reparto de los flujos energticos en las diferentes etapas del
ciclo.
Los porcentajes de las prdidas de calor hacia los gases de escape y el
refrigerante aumentan a medida que el vehculo se calienta y la demanda de
almacenamiento de calor en la estructura del motor y los fluidos disminuye,
tendiendo a estabilizarse el promedio una vez que se ha alcanzado la
temperatura de operacin regulada por el termostato.
129
130
131
Energa utilizada en el
calentador del habitculo (kW)
Q calentador_mean = 0,808 kW
2.5
2
1.5
1
0.5
0
0
400
800
1200
Tiempo (s)
2
Q egr_mean = 0,702 kW
1.5
0.5
0
0
400
800
1200
Tiempo (s)
132
133
Porcentaje de energa
recuperada (%)
50
25
0
0
400
800
1200
Tiempo (s)
134
Igual sucede con los sistemas modernos de enfriamiento del aire de admisin,
que llegan a disipar potencias trmicas del mismo orden que los radiadores
principales [3.8].
3.3.4. Observaciones
135
136
El modelo construido por Degraeuwe [3.13] para predecir los flujos de calor a
travs de las paredes de la cmara de combustin y las temperaturas de los
componentes se basa en la analoga existente entre los sistemas de flujo
elctricos y de transferencia de calor, como lo hicieron otros autores como
Bohac [3.45], Veshagh [3.46], o Gentile [3.18]. El motor es tratado como
una malla trmica, donde el calor puede transportarse entre capacitancias
concentradas a travs de conductancias concentradas (Reyes [3.44] las llama
globales). Para hacer corresponder el sistema real a su anlogo elctrico, el
pistn, la culata y el cilindro del motor se descomponen en trozos elementales
(nodos) de acuerdo al nmero de sensores de temperatura instalados en el
motor para validar el modelo y cuidando de observar el nmero de Biot para
asegurar el criterio de uniformidad de la temperatura en cada nodo. Estos
nodos se interconectan por resistencias trmicas de conduccin o conveccin,
estas ltimas asociadas a las condiciones de contorno (conveccin y contacto
variable) de los nodos. Las condiciones de contorno, se reitera, estn
representadas por la temperatura y coeficientes de conveccin de los gases
del cilindro, los gases de admisin y escape, el refrigerante y el aceite.
Adems cabe la condicin de contorno especial en el grupo pistnsegmentos-cilindro. La fundamentacin experimental, terica y la validacin
del modelo fueron ampliamente explicados por Degraeuwe [3.13] en su
trabajo de tesis doctoral.
3.4.1.1. Fundamentos del modelo nodal
137
camisa puede variarse automticamente, mientras que los nodos del pistn y
la culata son difciles de automatizar por lo complejo de la geometra de
ambas partes, como se ve en los modelos dibujados en las Figura 3.25. En el
modelo geomtrico el usuario es libre de discretizar estos elementos como
ms le convenga.
[C ] T = Q + [C ] T + H
[ t] [ ]
[K ] +
[ t +t ] [ ]
t
t
138
(3.14)
[K ] [T ] = [Q ] + [H ]
(3.15)
139
1
1
1
+ Ci 2 j .(x )Exp _ i 2 j
K i 2 j = K const
(3.16)
T gas =
Entre el gas y el
pistn o la culata:
1 720
Tgas ( )h( ).d
720 0
1 720
h( )d
720 0
1 720
h( )d
720 0
140
1 720
h( )d
720 0
Entre el plato de
vlvula y su asiento
K gas 2lin =
( )
( )
1 720
h( )Acyl _ n ( )d
720 0
( )
( )
z , = 0 si d pist head z
;
z , = 1 si d pist head > z
h (z ) =
1 720
h( ). ( ).d
720 0
Acyl _ n ( ) =
))
Entre el pistn y el
cilindro
1
1
1
K pis 2cyl = K seg
+ C pis 2cyl . S p Exp _ pis 2oil
Entre el cilindro
(culata) y el
refrigerante,
1
1
1
+ Ccyl 2cool .(Qcool )Exp _ cyl 2cool ;
K cyl 2cool = K cyl
basados en una
relacin Nu-Re-Pr,
1
K head 2cool = K head + Ccyl 2cool .(Qcool )Exp _ cyl 2cool
( )
0,8 0,4
.
Nu = 0, 023 Re D Pr
Klinoil = Aijhlinoil
Nu gal =
basados en una
relacin Nu-Re-Pr
Nu gal = C gal Re mgal :
hgal .d gal
koil
; Re gal =
S p .d gal
oil
( )
Entre superficies de
los puertos y los
gases:
Entre el pistn y el
cilindro, con Kseg
conductancia de los
anillos por unidad de
longitud:
t
D
K Pis _ i Lin _ j = con K seg Lin _ j
2
Tcycle
141
142
Cte_pis2oil
721,423
h_agua2cil
0,525
Exp_pis2oil
0,687
h_agua2cul_FD
1,219
Kpis2cil
3,876
h_cil2oil
864,869
h_agua2cul_Exh 0
143
144
Figura 3.28. Flujos de calor (en W) entre los nodos del modelo.
El flujo de calor recibido por el pistn diverge por dos caminos: por uno de
ellos el calor fluye hacia la camisa a travs de los segmentos (16 %), mientras
que por otro el calor fluye hacia el aceite a travs de la galera de lubricacin
y las superficies internas inferiores (84 %). Este calor recibido por el aceite se
transfiere al refrigerante y al medio circundante a travs del crter inferior
(para efectos del balance trmico, todo el calor recibido por el aceite se
considera transferido al refrigerante, como se coment en la segunda seccin
de este captulo). En el caso de la camisa, el calor proviene de los gases de
combustin y del pistn (a travs de los segmentos). Todo este calor se
145
146
Tabla 3.5. Resumen de los procesos transitorios de prueba del motor DV6.
Parmetro
Carga
(pme, bar)
Rango de variacin
Condiciones
Transicin de carga
Revoluciones
3,29 - 10,66
10,66 - 3,29
1000
3,9 - 16,86
16,86 - 3,85
1500
6,54 - 18,77
3,28 - 12,03
18,77 - 6,48
12,03 - 3,28
2000
4,18 - 12,64
12,64 - 4,14
2500
4,42 - 12,79
5,42 - 1,69
12,79 - 4,21
3000
4,22 - 12,22
12,22 - 4,21
3500
2,55 - 8,55
8,55 - 2,55
4000
Rgimen
(rpm)
2,41 - 8,81
8,81 - 2,37
4500
Transicin de revoluciones
Carga
1078 - 3013
3013 - 1078
4,04
1092 - 2710
1042 - 4001
2710 - 1091
4001 - 1029
8,08
147
Los resultados de las pruebas del modelo con este motor son muy
satisfactorios. A continuacin se describen algunos resultados para ilustrar la
validez del modelo.
El comportamiento transitorio de las temperaturas de la camisa del modelo
del motor se puede observar en las grficas de la Figura 3.29, en la que se
muestra la evolucin de la temperatura en dos puntos extremos a lo largo de
la camisa, distanciados 8,8 y 89,1 mm de la superficie de la culata, bajo dos
transitorios de carga: de 2,41 a 8,81 bar de presin media efectiva a un
rgimen constante de 4500 rpm y de 3,29 a 10,66 bar a 1000 rpm.
Figura 3.29.
Temperaturas
transitorias para un par
de nodos longitudinales,
a y B, extremos de la
camisa bajo dos
procesos transitorios: a)
4500 rpm (2 - 8 bar), b)
1000 rpm (3 - 10 bar).
148
149
150
Para algunos nodos del modelo las predicciones de temperaturas durante los
cambios de carga presentan grandes desviaciones con relacin a los valores
experimentales. Estas desviaciones estn asociadas tanto con las
incertidumbres (por ejemplo, el modelo asume que todo el calor proviene de
la combustin y no hace una discriminacin de los flujos de calor generados
por la friccin), como con las simplificaciones del modelo (reducido nmero
de nodos, correlaciones simplificadas y globales para valorar las condiciones
de contorno, etc.). En general, el modelo da una valoracin satisfactoria de la
evolucin de la temperatura y los flujos de calor durante los cambios de
operacin transitorios.
151
En el pargrafo 3.2 del presente captulo se coment que para este motor se
haban registrado las variables medias e instantneas operando bajo diferentes
condiciones de operacin estacionarias. Estas condiciones de operacin se
detallan en la tabla 3.7.
152
Rgimen
pme, bar
Rgimen
0
2
3
4
6
8
8,8
10
12
14
16
18
20
2500
2000, 3000
2500, 3000
1750, 2000, 2500
1750, 2000
2000, 2250
153
0 bar
2 bar
3 bar
4 bar
6 bar
20
8.8 bar
10 bar
10
12 bar
14 bar
16 bar
18 bar
Qrefr (kW)
40
30
0
700
1200
1700
2200
2700
3200
3700
4200
20 bar
(3.17)
154
Liner_40_79.2_0
Liner_40_26.4_2.125
Liner_40_61.6_2.125
Tm_cil_pred
Liner_40_26.4_8.5
Trefr_medida
Tac_medida
Tbl_medida
Tm_pist_pred
Tm_cul_ pred
155
Temperatura (C)
200
100
150
80
100
60
50
40
20
0
50
100
150
200
250
300
350
250
120
0
400
Tiempo (s)
156
200
Tm_cul_pred
Tac_medida
T_liner_int_wall44
Velocidad del vehculo
140
120
Temperatura (C)
160
100
120
80
60
80
40
40
Tm_pist_pred
Tm_cil_pred
Trefr_medida
Tblk_medida
20
0
0
200
400
600
800
1000
0
1200
Tiempo (s)
157
qcul_esc
qcul_adm
qw_refr
qw_aceite
qm_refr
5500
4500
3500
2500
1500
500
-500
200
400
600
800
1000
1200
Tiempo (s)
Figura 3.42. Componentes del balance energtico interno predicho del motor
durante el calentamiento bajo un ciclo de conduccin NEDC (considerado
para un cilindro).
Una comparacin de los flujos de calor al refrigerante durante el
calentamiento del motor a partir del modelo actualmente no es posible porque
158
3.6. Conclusiones
159
3.6. Conclusiones
En este captulo se ha realizado la caracterizacin trmica de un motor, a
partir de las pruebas experimentales y de la utilizacin de un programa
predictivo de temperaturas. Como resultado se tienen valores estimados de la
inercia trmica de la masa sensible interna del motor y del motor integrado
al sistema de refrigeracin, se dispone del mapa de disipacin de calor
registrado para un buen nmero de puntos de operacin y se cuenta con
valores de las temperaturas y los flujos de calor hacia las paredes, que sirven
como condiciones de contorno en otros modelos de gestin trmica del motor
y en modelos de prediccin y diagnstico de la combustin. Tambin se
cuenta con la caracterizacin hidrulica. Concretamente el conocimiento de la
inercia trmica del motor y los calores disipados a las paredes puede
emplearse en el diseo, modelado y simulacin de los sistemas de
refrigeracin.
La informacin de las predicciones tambin se ha usado para obtener una
correlacin lineal del calor disipado al refrigerante en funcin de la presin
media efectiva, la velocidad del motor y la temperatura del refrigerante. Un
resultado importante de este captulo ha sido tambin la determinacin
experimental del balance trmico del motor durante su calentamiento bajo
carga constante y bajo un ciclo de conduccin estndar, tras lo cual se han
sacado conclusiones referentes a los factores que influyen sobre el
calentamiento del motor. Ha quedado clara la influencia de la masa como
resistencia dinmica importante durante el calentamiento y se han calculado
los calores transportados y las temperaturas del refrigerante y los gases de
escape, como informacin inicial utilizable en los proyectos de recuperacin
de energa en los vehculos.
Las variables medias e instantneas adquiridas durante el calentamiento del
motor bajo carga constante y bajo ciclo de operacin pudieron emplearse en
160
valide con mayor resolucin las temperaturas del pistn, puede ser
empleando valores experimentales de un motor similar de
dimensiones y masas conocidas;
3.7. Bibliografa
161
3.7. Bibliografa
[3.1]
162
3.7. Bibliografa
163
164
[3.38] Reichelt, L., Meingast, U., Renz, U. Calculating transient wall heat
flux from measurements of surface temperature. Int. J. Heat Mass
Transfer, vol. 45, pp. 579-584, 2002.
[3.39] Morel, T., Keribar, R. Heat radiation in D. I. Diesel engines. SAE
paper 860445.
[3.40] Descombes, G., Maroteaux, F, Feidt, M. Study of the interaction
between mechanical energy and heat exchanges applied to IC engines.
Appl. Therm. Eng., vol. 23, issue 16, 2003, pp. 2061-20078.
[3.41] Alkidas, A. C. Performance and emissions acheivements with an
uncooled heavy_duty, single_cylinder Diesel engine. SAE paper
890144.
[3.42] Wolf-Heinrich, H. Aerodynamics of road vehicles. SAE, Warrandale,
1998.
[3.43] Woschni, G. A Universally applicable equation for the instantaneous
heat transfer coefficient in the I.C. engine. SAE paper 670931.
[3.44] Reyes M. Desarrollo de un modelo para el anlisis trmico de MCIA.
Seminario en el CMT, Universidad Politcnica de Valencia. Espaa.
Septiembre, 2003.
[3.45] Bohac, S.V., Baker, D.M., Assanis D.N. A global model for steadystate and transient S.I. engine heat transfer studies. SAE Paper
960073.
[3.46] Veshagh, A., Chen, C. A computer model for thermofluid analysis of
engine warm-up process. SAE paper 931157.
[3.47] Stone, R. Introduction to Internal Combustion Engines. MacMillan
publishers. 1999.
[3.48] Morel, T., Blumberg, P.N., Fort, E.F., Keribar, R. Methods for heat
analysis and temperature field analysis of the insulated Diesel. NASA
CR-174783, 1984.
[3.49] Lapuerta, M., Armas, O., Hernndez, J.J. Diagnosis of DI Diesel
combustion from in-cylinder pressure signal by estimation of mean
thermodynamic properties of the gas. Appl. Therm. Eng., vol. 19 pp.
513-529, 1999.
Captulo 4
Consideraciones y herramientas de
diseo y modelado de los sistemas de
refrigeracin de los motores.
Resumen ....................................................................................................... 167
4.1. Introduccin .......................................................................................... 167
4.2. Prediseo de los sistemas de refrigeracin de los motores ................... 169
4.2.1. Clculo preliminar de los componentes del sistema de
refrigeracin ............................................................................................. 170
4.2.2. Metodologa de diseo de los sistemas de refrigeracin................ 184
4.2.3. Clculo de los radiadores de los sistemas de refrigeracin............ 187
4.3. Modelado de los sistemas de refrigeracin ........................................... 196
4.3.1. Modelos simplificados de la temperatura del motor a nivel de
sistema...................................................................................................... 197
4.3.2. Modelos a nivel de circuito del sistema de refrigeracin del
motor ........................................................................................................ 199
4.4. Herramientas de modelado unidimensional .......................................... 208
4.5. Conclusiones ......................................................................................... 212
4.6. Bibliografa............................................................................................ 213
166
4.1. Introduccin
167
Resumen
En el captulo anterior, aparte del estudio de los flujos de calor y las
temperaturas en el interior del motor como componente, fue analizado el
balance trmico de ste acoplado a su circuito de refrigeracin, as como
tambin se pudo ver el tiempo que tarda el motor en calentarse y se
conocieron los factores de los que depende ese calentamiento. Pensadas estas
actividades como diagnstico de la respuesta de los sistemas de refrigeracin
modernos, cabe preguntarse qu trabajos pueden emprenderse para mejorar la
respuesta trmica de sistemas ya diseados o en operacin y qu herramientas
de computacin pueden integrarse a las tareas de rediseo, modelado y
simulacin de nuevos sistemas de refrigeracin. En los captulos quinto y
sexto de esta tesis se da respuesta a esos interrogantes, particularizando dos
tareas: la modificacin al sistema original del motor probado para mejorar su
calentamiento, y la concepcin, modelado y simulacin original del sistema
de refrigeracin de un motor. Como preparacin a esos trabajos, en el
presente captulo se describe el procedimiento de diseo de los sistemas de
refrigeracin, se introducen los fundamentos termohidrulicos que estn en la
base de los programas comerciales de modelado de estos sistemas y se
describe brevemente la utilizacin de los mismos.
4.1. Introduccin
El diseo funcional de los sistemas de refrigeracin parte de la necesidad de
proveer la suficiente capacidad de evacuacin de calor demandada por el
motor, y otras posibles fuentes, a intercambiadores de calor externos al motor,
como radiadores, recuperadores de calor, calentadores de cabinas, etc.,
conservando las temperaturas de los fluidos del sistema y las partes metlicas
dentro de unos valores de control propicios desde el punto de vista del
consumo de combustible y las emisiones, ajustndose a un balance energtico
entre fuentes y sumideros como se sintetiza en el cuadro de la Figura 4.1.
Al inicio slo se conoce la aplicacin, el motor que propulsa el vehculo o
medio de transporte (en general sin caractersticas trmicas), las opciones
para el posible trazado y ubicacin de su sistema hidrulico, el espacio y la
disposicin para los radiadores y los ventiladores y las condiciones
ambientales de operacin. Lo que se tiene que hacer es seleccionar los
radiadores y ventiladores adecuados para evacuar el calor disipado por las
168
169
170
(4.1)
171
2, 50 k g Re
0,7
Pr
0,3
4Vd n 2
D
n
4
4
T f Tw ) + 5, 98
(Tg Tw ) + 0,4341
(
2
V n 2
2
(4.2)
donde es la constante de Stefan-Boltzman; D es el dimetro del
cilindro; kg, Pr son la conductividad y el nmero de Prandtl de los
gases en el cilindro, respectivamente; Re es el nmero de Reynolds
utilizado en el clculo de la transferencia de calor entre los gases y las
paredes del cilindro; T g y T f son las temperaturas medias de los
gases y de la llama, respectivamente; y Tw es la temperatura de la
pared de la cmara de combustin.
Parish ha llegado a esta expresin partiendo de una extensin de la
formulacin de Annand [4.33] para el calor disipado por las paredes de la
cmara de combustin - Q = h A (T T ) + h A (T 4 T 4 ) + N -. La
refr
172
Tw = 110 C ), pueden
173
174
= 0,6638N
Q
refr
e
1200
R 2 = 0,86
600
= 0,2897N
Q
refr
e
R 2 = 0,846
800
1800
= 0,2037N
Q
refr
e
R 2 = 0,8167
400
0
0
1200
2400
3600
1200
2400
3600
175
(4.5)
donde Caire y aire son el calor especfico y la densidad del aire calculados
en
para una temperatura media del aire de entrada ambiente de Taire
= 40 C
(segn Zhelezko [4.35], sta es la temperatura promedio del aire a la entrada
al radiador).
El caudal de aire (en m3/s) para motores Diesel puede tomarse en el rango
(0,015...0,02)Ne [4.35], siendo Ne la potencia efectiva del motor. Segn un
proveedor de maquinaria pesada, para la refrigeracin adecuada de un motor
Diesel, se requiere un caudal mnimo de aire de 0,04 m3/s por cada kilovatio
de potencia neta del motor [4.34]. Otra posibilidad es fijar el valor del Taire
( Taire = 45 55 C ).
4. El gasto volumtrico de refrigerante que pasa a travs del radiador es:
176
(4.6)
(4.7)
(4.8)
(4.9)
(4.10)
177
(4.11)
Donde vaire es la velocidad del aire en el frente del radiador (vaire = 6...18 m/s
[4.4]) sin tener en cuenta la velocidad de movimiento del vehculo.
9. La profundidad del panel del radiador, en mm:
lrad = Frad
(F )
fr
rad
(4.12)
(4.13)
178
Figura 4.5. Recorrido convencional del aire del ambiente hacia el radiador
de un vehculo.
En los vehculos modernos el volumen de los radiadores y el diseo de los
pasajes de aire est comprometido por la configuracin de los motores, la
integracin de nuevos componentes como los requeridos para el tratamiento
de emisiones, la refrigeracin de la carga fresca, la refrigeracin del aceite, la
ubicacin del condensador del sistema de aire acondicionado del vehculo y
las caractersticas estticas y aerodinmicas del vehculo.
El mostrado es un orden para calcular los parmetros iniciales de diseo o
seleccin del radiador. La metodologa bien podra basarse estrictamente en
el espacio disponible solamente y partirse de all. Con los parmetros de
diseo iniciales, que se reducen a los gastos msicos de refrigerante y de aire,
179
Figura 4.6. Esquema simplificado de las resistencias al paso del aire por el
sistema de ventilacin del motor y los radiadores.
En general, para los sistemas de ventilacin de los motores, la presin de los
ventiladores est en el rango de pai = 600....1000 Pa.
La velocidad circular del aspa del ventilador, en m/s, en el borde externo de
se determina por [4.36]:
u = a pai / ai
(4.14)
donde a es un coeficiente que depende de la forma de las aspas: a = 2,2 ...
2,9 para aspas curvas, a = 2,8 ... 3,5 para aspas planas. Partiendo de
premisas de carcter acstico se toma u 110 m/s [4.36].
180
Vaire
'
vaire
(4.15)
fd =
2pvent
2Vvent
;
=
fd
v 2ps
Dv2 v 2ps
(4.18)
181
Figura 4.8. Relacin entre la velocidad de paso del aire y la velocidad del
vehculo para diferentes estados de encendido del ventilador en funcin de la
velocidad del vehculo [4.3].
182
V
Vrefer _ c = referl
(4.19)
183
(4.21)
184
Presin (bar)
0.25
DP_Sistema
0.2
Bomba_500
Bomba_780
0.15
Bomba_1200
Bomba_1800
0.1
Bomba_2200
0.05
0
0
10
20
30
40
50
60
70
80
90
100
185
186
187
188
189
11.32-737-SR
Unidades SI
Escalonado
Serpentn
0,737
0,1
11,32
0,138
0,004
18,72 mm
2,54 mm
0,446 aletas/mm
3,51 mm
0,102 mm
0,78
0,78
270
0,886 mm2/mm3
0,845
0,225
0,845
5,715 mm
0,001
0,254 mm
0,789
0,55
20,066 mm
13,97 mm
190
(4.22)
p ; Cmn = (mC
p ) mn ; Cmx = (mC
p ) mx
Siendo C = mC
sal
en
y Taire
son las temperaturas de salida y entrada del fluido
donde Taire
sal
en
y Trefr
son las temperaturas de salida y
fro, respectivamente, y Trefr
h / (siendo Dh el
Se calculan los nmeros de Reynolds, Re = mD
dimetro hidrulico y m el gasto msico del fluido correspondiente),
requeridos para determinar las caractersticas adimensionales de
transferencia de calor y de friccin de las superficies del radiador. A
continuacin se calculan los nmeros de Nusselt ( Nu = hDh k ) y
Colbourn ( j = St Pr 2 3 , con el nmero de Stanton dado por
St = Nu Re Pr ).
191
(4.24)
(4.25)
NTU mx =
192
Aa U a
p min )
(mc
(4.26)
C NTU 0,22
(4.27)
(4.28)
en
en
Trefr
y Taire
son las temperaturas de entrada del refrigerante y del aire,
respectivamente.
en
sal
Crefr (Trefr
Trefr
)
en
en
Cmn (Trefr
Taire
)
(4.29)
sal
sal
y Taire
. Si las temperaturas de salida son muy
Se obtienen as Trefr
(4.30)
193
(4.31)
194
Modificacin de los
parmetros geomtricos
( )
prefr = f Vrefr
paire = f Vaire
195
196
197
dT
= Qrefr Q r
dt
(4.32)
d T
= Q refr Q r
dt
(4.33)
198
Q refr =
Q refr
T
T +
Q refr
Q refr
h
(4.34)
refr
refr
refr
Q
Q
Q
en
(4.35)
Q r = r T + r m + enr Trefr
T
m
Trefr
En los sistemas de refrigeracin convencionales la variacin del caudal m
del refrigerante depende de la frecuencia de rotacin del rotor de la bomba y
del desplazamiento x del rgano de control (la vlvula de control
usualmente es el termostato).
en
Reemplazando esta expresin en la frmula (4.34), considerando Trefr
=0 y
Q refr Q r
Q refr
Q
d T
+ Fe T = r x +
h +
dt
x
h
(4.36)
199
200
(4.37,a)
(4.37,b)
m = 0
i =1
(4.39)
202
203
Primera zona
La primera zona se caracteriza por las dimensiones y la masa de la pared del
cilindro, cuya temperatura Tw puede obtenerse tras establecer el balance de
calor presentado en la ecuacin (4.40):
dTw Q g , w Q w,refr
=
dt
mwc pw
(4.40)
Q g , w = hg Aw (Tg Tw )
(4.41)
del lado del refrigerante, Tw , refr , del coeficiente de conveccin, hrefr y del rea
de conveccin, Ac :
k A
m
Q w,refr = w w (Tw Tw,refr ) = hrefr Ac (Tw,refr Trefr
)
(4.42)
204
sal
en
Trefr
+ Trefr
(4.43)
hc = c Re0,8 Pr 0,4
D
(4.44)
4 m cyl
refr refr D
(4.45)
cinemtica del
refrigerante, respectivamente.
Segunda zona
La segunda zona corresponde al refrigerante en la camisa de refrigeracin del
motor. La temperatura en esta zona puede asumirse aproximadamente igual a
la temperatura de salida del motor, con su variacin descrita por la ecuacin
diferencial [4.14]:
sal
dTrefr
dt
(4.46)
205
(4.47)
206
(4.48)
(4.51)
207
caracterstica.
Durante el calentamiento del motor (cerrado el paso por el radiador), los
flujos por los ramales principales, el del by-pass y del calentador del
habitculo se unen antes de entrar al motor, con lo que se puede calcular una
temperatura promedio:
en
Trefr
=
sal
sal
Tbysal m by + Texp.
tn mexp.tn + TEGR mhc
f (t )
m refr
(4.52)
208
209
210
Figura 4.17.
Variables calculadas
en los nodos de las
redes de los circuitos
termohidrulicos.
Para los trabajos de esta tesis se han utilizado los programas AMESim, que
basa sus soluciones en las tcnicas bond graph y multipuerto, y FlowMaster,
tambin multipuerto, pero que hace uso del mtodo de las caractersticas para
resolver las ecuaciones. Para el modelado del sistema de ventilacin del
motor estudiado en el captulo sexto se ha utilizado el programa PowerFlow.
En el apndice A se presenta un breve resumen de estos programas.
211
212
4.5. Conclusiones
El presente ha sido un captulo de transicin en el que, como parte sustancial
de la conceptualizacin de los sistemas de refrigeracin de los motores, se ha
abordado la composicin y prediseo de los componentes de estos sistemas,
desde una perspectiva convencional. Se ha retomado la valoracin del calor
evacuado al sistema de refrigeracin como factor clave para el prediseo del
sistema resumiendo el modelo aportado por Parish [4.1] y se ha sintetizado la
metodologa de clculo trmico de diseo de los radiadores fundamentada en
4.6. Bibliografa
213
4.6. Bibliografa
[4.1]
[4.2]
[4.3]
[4.4]
[4.5]
[4.6]
[4.7]
[4.8]
[4.9]
214
[4.10] Tafti, D.K. and Cui, J. Computations of flow and heat transfer in a
three dimensional multilouvered fin geometry. Int. J. Heat and
Fluid Flow, vol. 45, n 25, pp. 5007-5023, 2002.
[4.11] Kakac S., Liu H.: Heat exchangers selection, rating and thermal
design, CRC Press LLC, 2002.
[4.12] Oliet, C., Oliva, A., Castro, J. y Prez-Segarra, C.D. Parametric
studies on automotive radiators. Appl. Therm. Eng., vol. 27, pp.
2033-2043, 2007.
[4.13] Kaplan, J. A., Heywood, J. B. Modelling the spark ignition engine
warm-up process to predict component temperatures and
hydrocarbon emissions. SAE paper 910302.
[4.14] Cortona, E., Onder, C. H. Engine thermal management with
electronic cooling pump. SAE paper 2000-01-0965.
[4.15] Chen, C., Veshagh, A. A one-dimensional model for in-cylinder
heat convection based on the boundary layer theory. SAE paper
921733.
[4.16] Mohan, K. V, Arici, O., Yang, S., Johnson, J. H. A computer
simulation of the turbocharged Diesel engine as an enhancement
of the vehicle engine cooling system simulation. SAE paper
971804.
[4.17] Lehner, C., Parker, G., Arici, O. Johnson, J. Design and
development of a model based feedback controlled cooling
system for heavy duty diesel truck applications using a vehicle
engine cooling system simulation. SAE paper 2001-01-0336.
[4.18] Jarrier, L. Champoussin, J. C., Yu, R., Gentile, D. Warm-up of a
Diesel engine: experiment and modeling. SAE paper 2000-010299.
[4.19] Torregrosa, A. J., Olmeda, P., Degraeuwe, B., Reyes, M. A
concise wall temperature model for DI diesel engines. Appl.
Therm. Eng., vol. 26, 1320-1327, 2006.
[4.20] Incropera, F. P., DeWitt, D. P. Introduction to Heat Transfer. John
Wiley, New York, 1996.
[4.21] Strupp, N., Lemke, N., Tegethoff, W., Khler, J. MAC simulation
using Modelica/Dymola. 7th IIR Gustav Lorentzen Conf. Natural
Working Fluids, Trondheim, Norway, 2006.
4.6. Bibliografa
215
216
Captulo 5
Modelado, simulacin y validacin
del sistema de refrigeracin de un
motor HSDI para reducir su tiempo
de calentamiento
5.1. Introduccin .......................................................................................... 219
5.2. Modelado trmico del motor. ................................................................ 221
5.2.1. Descripcin y caracterizacin del sistema de refrigeracin. .......... 221
5.2.2. Modelo 1D del sistema de refrigeracin original del motor. ......... 226
5.3. Modelado de las modificaciones al sistema de refrigeracin del
motor para reducir el tiempo de calentamiento............................................ 231
5.4. Validacin experimental de las modificaciones propuestas y
valoracin de su influencia sobre las emisiones y el consumo de
combustible .................................................................................................. 236
5.4.1. Flujos de refrigerante experimentales de todas las
configuraciones. ....................................................................................... 237
5.4.2. Respuesta de temperatura del refrigerante. .................................... 238
5.4.3. Emisiones. ...................................................................................... 241
5.5. Conclusiones ......................................................................................... 244
5.6. Bibliografa............................................................................................ 246
218
5.1. Introduccin
219
5.1. Introduccin
El diseo de los sistemas de refrigeracin parte de la necesidad de proveer la
suficiente capacidad de disipacin de calor a plena carga, no obstante que la
mayor parte del tiempo los vehculos son conducidos a condiciones parciales
de carga. As, estando sobredimensionados, a cargas parciales de operacin
los sistemas de refrigeracin tradicionales son responsables de prdidas
innecesarias de calor, particularmente durante el calentamiento del motor. Las
medidas tecnolgicas tomadas para adaptar el motor a las condiciones
actuales de operacin pueden dar ganancias en la eficiencia del motor (es
decir, menor consumo de combustible), lo mismo que una reduccin de
emisiones nocivas. En trabajos publicados recientemente, como los de
Campbell et al. [5.5], Allen y Lasecki [5.6], y Pang y Brace [5.7], ya se han
venido demostrando los beneficios de las configuraciones de sistemas de
refrigeracin con control flexible. Las mejoras propuestas van desde la
introduccin de bombas controladas electrnicamente (Cortona y Onder [5.8],
Cho et al. [5.9]), ventiladores, persianas y termostatos (Chanfreau et al.
[5.10], Chalgren y Barron [5.11]), en los cuales se explota el potencial
ofrecido por las tecnologas de control computarizado, hasta la
implementacin de diferentes tcnicas como la refrigeracin de precisin en
Robinson et al. [5.12], la ebullicin nucleada en Campbell et al. [5.2], la
refrigeracin evaporativa en Porot et al. [5.14], o el aislamiento trmico de las
paredes, revisado por Jaichandar y Tamilporai [5.15]. Adicionalmente, la
importancia de la temperatura tanto del aire de admisin como del
refrigerante sobre las emisiones de un motor Diesel, bajo condiciones estables
de operacin en el motor caliente, han sido investigados recientemente por
Torregrosa et al. [5.16].
La introduccin de estrategias tendientes a reducir el tiempo de calentamiento
ha adquirido especial atencin tambin porque cuanto menos tiempo opere el
motor en regmenes ineficientes y con temperaturas de sus componentes por
debajo de los umbrales ptimos de operacin, menor ser el desgaste
experimentado por el motor, como lo anotan Broatch et al. [5.1].
El desarrollo de herramientas computacionales con capacidad para modelar y
analizar sistemas termohidrulicos con control integrado ha permitido en los
ltimos 10 aos avanzar mucho en el estudio de los sistemas de refrigeracin,
de manera particular. Recurriendo a la simulacin numrica son ostensibles
220
221
222
223
ste, mientras que en la Figura 5.3 se pueden observar las proporciones de los
flujos de refrigerante, V&refr , a travs de los diferentes ramales del circuito en
relacin con el flujo promovido por la bomba.
Figura 5.2.
Caracterstica presin
vs. flujo volumtrico de
la bomba del motor
ensayado a diferentes
revoluciones.
224
Prdida de
presin, p, Pa
12200
700
1800
700
1900
24200
23100
500
Aequiv,
mm2
61,75
257,77
160,75
257,77
156,46
87,10
185,13
833,38
52,0
56,5
125,0
400
20300
30000
968,97
148
268,96
200
500
800
1200
1600
3,92
4,86
5,43
5,88
4,02
5,08
5,93
6,7
4,05
5,16
6,16
6,99
4,1
5,28
6,23
7,18
4,16
5,31
6,32
7,24
400
500
600
225
6,61
7,15
7,63
8,02
9,13
9,98
8,43
9,67
10,72
8,67
9,98
11,13
8,84
10,26
11,37
1,5
2,5
4,0
5,5
7,0
7,5
23
25
26
27
31,5
34
37
39,5
40
43
50
55
46
51,5
58,5
67
50
58
67
74
52,5
62,5
73
82,5
En la tabla 5.4 se presentan los flujos de calor transferidos por los gases de
escape al refrigerante en el enfriador de EGR, obtenidos a partir de las
mediciones de calor transferido al refrigerante, la presin media efectiva
calculada y las correspondientes revoluciones del motor, en kW.
Tabla 5.4. Flujo de calor transferido por los gases de escape al refrigerante
en el enfriador de EGR, en kW.
Rgimen del motor, rpm
pme, bar
2,2
3,8
5,0
6,4
1635
1728
1890
2040
0,98
1,1
1,12
1,14
0,7
0,763
0,909
1,109
0,667
0,757
0,839
0,936
0,485
0,7
0,76
0,8
226
227
228
229
Figura 5.7.
Comparacin de las
temperaturas
experimental y
modelada del
refrigerante
durante el
calentamiento del
motor.
Como se puede ver en la Figura 5.7, el error del modelo comparando con la
respuesta experimental durante el calentamiento bajo condiciones de
operacin constantes es menor del 4 %, excepto para los primeros 10
segundos, tiempo durante el cual el motor an est operando a ralent, en el
ensayo experimental. El tiempo requerido para alcanzar los 80 C es de 400,5
segundos, mientras que el correspondiente tiempo experimental es de 395
segundos.
Una vez ajustado el modelo bajo condiciones de operacin estabilizadas, y
que la respuesta hidrulica y trmica del sistema bajo condiciones de carga
constante se consider satisfactoria, se cambi el modo de operacin del
modelo a condiciones de conduccin bajo ciclo de operacin NEDC. El calor
suministrado durante el ciclo de conduccin tiene el comportamiento
ilustrado en la Figura 5.8.
230
231
Figura 5.9.
Temperatura del
refrigerante a la
salida del motor
durante la operacin
de ste bajo un ciclo
de conduccin
NEDC.
232
233
M2
Cerrada
Cerrada
M3
Cerrada
Abierta
M4
Abierta
Cerrada
234
235
M1
M2
M3
M4
Flujo de refrigerante
236
Figura 5.11.
Temperatura del
refrigerante a la salida
del motor para las
cuatro configuraciones
modeladas del sistema
de refrigeracin, bajo
el ciclo de conduccin
NEDC.
237
238
239
10
90
Temperatura (C)
70
60
M1
M3
M1
M3
50
40
M2
M4
M2
M4
30
20
80
0
10
0
0
200
400
600
800
1000
-2
1200
Tiempo (s)
240
Experimental, (s)
CI 0.95 (s)
Con respecto a M1
Modelada
M1
M2
M3
M4
713,4 9,18
554,5 9,03
716,0 5,23
641,8 5,41
-22,3 %
0,4 % ~ 0 %
-10,0 %
721
539
730
655
241
242
243
0.5%
0.0%
-0.5%
-1.0%
-1.5%
-2.0%
-2.5%
-3.0%
-3.5%
-4.0%
-4.5%
-5.0%
M2
M3
M4
ECE 1
ECE 2
ECE 3
ECE 4
EUDC
TOTAL
-0.16%
-0.46%
-4.49%
-2.15%
-3.77%
0.03%
-2.58%
-1.38%
-0.61%
-0.39%
-1.62%
-0.68%
-1.59%
-1.39%
-0.84%
-0.45%
-0.15%
-0.66%
NOx
CO
HC
PM
Consumo de
combustible
M2
-1,28
-6,56
-11,41
11,09
-2,91
M3
4,97
0,81
- 0,56
1,21
-0,83
M4
1,08
-3,37
0,44
-13,41
-1,31
6,68
5,63
1,07
M1
5,85 4,34
244
5.5. Conclusiones
En este captulo se ha presentado el estudio del calentamiento de un motor
HSDI, cuyo balance trmico y caractersticas termohidrulicas fueron
determinados en el captulo tercero. Inicialmente se ha modelado el sistema
de refrigeracin con su configuracin original y se ha simulado su respuesta
bajo una carga y velocidad constantes del motor, con lo que se ha podido
acercar el comportamiento hidrulico y trmico al del sistema real. Se ha
explicado despus cmo el sistema original puede modificarse de manera
sencilla para que el motor caliente ms rpido, demostrando con la
simulacin y luego con la experimentacin esas posibilidades de mejora.
El modelado y la simulacin se han realizado con la ayuda del programa
computacional Amesim, los fundamentos de cuyos submodelos fueron
brevemente descritos en el captulo cuarto. En los modelos se han probado
cuatro configuraciones considerando slo aperturas y cierres totales de las
vlvulas, pero est claro que con, los modelos, las vlvulas bien pueden
modelarse para que sus secciones de paso se puedan variar de manera
continua.
Con ayuda del modelo desarrollado se estudiaron las respuestas trmicas de
las configuraciones propuestas para reducir el tiempo de calentamiento. Sin
afectar el diseo original del sistema de refrigeracin, slo incorporando en el
circuito original dos vlvulas de bajo costo, sin necesidad de estrategias
especiales y costosas de control, se ha concluido a partir de los resultados del
modelo la posibilidad de reducir el tiempo de calentamiento del motor tanto
bajo condiciones constantes de operacin, como durante la conduccin bajo
un ciclo de conduccin NEDC. Las modificaciones propuestas se
implementaron fsicamente en el motor real y se aprovech para evaluar las
5.5. Conclusiones
245
246
5.6. Bibliografa
[5.1]
[5.2]
[5.3]
[5.4]
[5.5]
[5.6]
[5.7]
[5.8]
[5.9]
5.6. Bibliografa
[5.10]
[5.11]
[5.12]
[5.13]
[5.14]
[5.15]
[5.16]
[5.17]
[5.18]
[5.19]
[5.20]
[5.21]
[5.22]
[5.23]
247
[5.24]
[5.25]
[5.26]
[5.27]
[5.28]
[5.29]
[5.30]
[5.31]
[5.32]
[5.33]
[5.34]
248
Captulo 6
Diseo, modelado y simulacin del
sistema de refrigeracin del motor
de una locomotora
Resumen ....................................................................................................... 251
6.1. Introduccin .......................................................................................... 251
6.2. Definicin de las especificaciones del motor y diseo preliminar
del sistema de refrigeracin.......................................................................... 254
6.2.1. Especificaciones de diseo del sistema de refrigeracin (tarea
tcnica) ..................................................................................................... 254
6.2.2. Clculo preliminar de los componentes del sistema de
refrigeracin ............................................................................................. 258
6.2.3. Disposicin de los componentes del sistema hidrulico y del
sistema de ventilacin .............................................................................. 266
6.3. Modelado y simulacin del sistema de refrigeracin............................ 269
6.3.1. Modelado y caracterizacin del sistema de ventilacin ................. 269
6.3.2. Caracterizacin de los componentes de los circuitos del
refrigerante ............................................................................................... 275
6.3.3. Modelado en FlowMaster del sistema de refrigeracin del
motor de la locomotora ........................................................................... 280
6.3.4. Alcance de resultados de los modelos............................................ 287
250
6.1. Introduccin
251
Resumen
En los captulos anteriores se ha estudiado la transferencia de calor en los
motores de automviles y se ha hecho uso de la informacin experimental de
operacin del motor para predecir las temperaturas y los flujos de calor en la
estructura del motor y en su circuito de refrigeracin, empleando un modelo
de prediccin de temperaturas y luego un programa de modelado 1D. La parte
conclusiva de esta tesis tiene que ver con el diseo asistido por herramientas
de modelado y simulacin de los sistemas de refrigeracin de motores Diesel.
A partir de las especificaciones del motor que propulsa una locomotora, la
caracterstica de la bomba y el punto de operacin de sta, se estima el calor
que debe evacuarse, se seleccionan el ventilador y el radiador adecuados a esa
capacidad trmica demandada y ajustados a las restricciones de espacio y
ubicacin dadas por la aplicacin particular, se modela el sistema as
ensamblado y se asegura el cumplimiento de las temperaturas del refrigerante
y del aire para las condiciones de plena carga del motor. Despus de esto, se
ensambla el sistema de refrigeracin real, se realizan las mediciones y se
comparan los resultados experimentales con los reales. El sistema de
refrigeracin diseado cumple con las especificaciones requeridas y la
metodologa de diseo se prueba como contrapropuesta a la prctica de
ensamble de componentes basada totalmente en el ensayo, dejando en
evidencia las ventajas del diseo integrado de las partes del refrigerante y del
aire asistido por las herramientas computacionales.
6.1. Introduccin
Actualmente slo la industria automotriz empieza a introducir metodologas
automatizadas para el diseo de sus sistemas refrigeracin [6.1]-[6.5], siendo
el diseo en aplicaciones para camiones y otros sistemas de transporte pesado
una actividad an no integrada o cerrada bajo una metodologa nica. Por
otro, lado la tendencia en el diseo de los motores de combustin alternativos
va haciendo que la casa matriz fabricante disponga la misma base estructural
compuesta por bloque, crter, culata(s), mecanismos manivela-deslizador y
de distribucin de gases y sistemas de alimentacin de aire y combustible
para un diverso mundo de aplicaciones y que slo los sistemas de control,
incluyendo el control de la temperatura del motor, dependan de la aplicacin
252
6.1. Introduccin
253
Figura 6.1.
Esquemas de
ventilacin
empleados en los
sistemas de
refrigeracin de las
locomotoras [6.7].
254
255
Enfriador de aceite
Configuracin
Disipacin estimada de calor, kW
Flujo estimado de refrigerante, l/s
Carcasa y tubos
439,21
57,72
Tanque de expansin
Capacidad, l
Flujo estimado de refrigerante, l/s
150
57,72
Compresor
Disipacin estimada de calor, kW
105,41
256
Cabeza (m)
50
40
& 2 + 308,03V
& + 54,22
H = -7936,4V
R 2 = 0,9733
n b = 950 rpm
30
20
10
0
0.02
0.04
0.06
0.08
Cabeza (m)
40
30
& 2 + 552,03V
& + 41,92
H = -46404V
R 2 = 0,994
n b = 950 rpm
20
10
0
0.01
0.02
3
257
258
259
35
7
35
12
85
85
35
35
1,015
1,015
1050
1050
1,18
1,18
en _ rad
Taire
, C
260
3,55
3,55
3500
0,0006
120
3500
0,0006
140
120
6
7
170
6
7
27
0,7
1000
0,75
2,5
27
0,7
1000
0,75
2,5
C
Temperatura media del aire,
Coeficiente global de transferencia de calor, Urefr,
W/m2K
rea de superficie de refrigeracin del radiador, Frad,
m2
rea frontal del radiador, Ffr, m2
819,32
19,55
31,40
313,10
7,47
7,00
81,50
79,00
52,50
52,50
99,23
131,19
434,57
2,79
144,64
1,07
261
222,34
78,47
1,46
1024
45,6
0,62
193,65
78,47
0,97
1549
19,9
0,21
700,00
5,87
700,00
5,87
262
263
264
265
800
700
600
500
400
300
200
100
0
30
50
70
90
Trefr_en_rad (C)
Figura 6.10. Variacin del calor disipado por el radiador a medida que el
motor calienta (la condicin de arranque con plena carga es hipottica,
pues el control del motor no permite su operacin a plena carga estando
an fro).
266
Figura 6.12.
Caracterstica del
ventilador
seleccionado.
En esta parte del proceso se define una configuracin provisional del circuito
con dimetros y longitudes tentativas de las tuberas y conexiones requeridas
segn la topologa y las limitaciones impuestas por los componentes.
267
268
Figura 6.14. Vista de ensamble parcial de las conexiones entre el motor y los
radiadores, en la que tambin se aprecian las bombas y el enfriador de
aceite.
269
6.3. Modelado
refrigeracin
simulacin
del
sistema
de
270
Los radiadores se modelan como un medio poroso con rea frontal Afr = 2,71
m2, con una disipacin de calor de 747,33 kW para el radiador del motor y
325,30 kW para el radiador del aftercooler .
En la Figura 6.17 y la Figura 6.18 se ilustran los campos de velocidades de
paso del aire antes y despus de la rejilla.
271
272
273
Los valores medios de las prdidas de presin del aire a su paso por el
sistema de ventilacin segn el esquema del sistema mostrado en la Figura
6.25 se resumen en la Tabla 6.4.
p, Pa
T, C
288
183,5
560
40,3
15,46
341
30,9
12,61
218,5
542,5
5-6 (ventilador1)
1307
5-6 (ventilador2)
1325
274
Los valores de las prdidas de presin sirven para calcular los coeficientes de
prdidas de presin, K (considerando la expresin p = KV& 2 ), que debern
introducirse en las caractersticas del sistema de ventilacin en el modelado
1D. Pueden calcularse los coeficientes individuales de los diferentes tramos,
discriminando las prdidas de presin por el radiador, o calcular un solo
coeficiente integral de prdidas y atriburselo al radiador. Esto hace ms
liviano el clculo y el error ser admisible porque las prdidas provienen en
rigor de los clculos CFD.
Los coeficientes de prdidas permiten obtener las curvas de prdidas de
presin del sistema de ventilacin que han de introducirse al programa de
modelado 1D y que se ilustran en la Figura 6.26.
275
Las bombas de los circuitos del motor (bomba 16G) y del enfriador de aire
(bomba D) se caracterizan con sus curvas presin - caudal, segn la
informacin dada por el fabricante. Graficando en las caractersticas de las
bombas los caudales nominales de operacin para los circuitos (0,0631 m3/s y
0,02524 m3/s para las bombas 16G y D, respectivamente), se obtienen las
276
277
278
279
3, 7 Re
donde Re es el nmero de Reynolds, hallado como Re = vDt , v es la
velocidad del flujo, v = V&refr At , es la rugosidad relativa. El coeficiente de
Colebrook permite hallar las prdidas de presin: p = Lt v 2 2 Dt y,
finalmente puede encontrarse el coeficiente de prdidas de presin
K = p V&refr 2 para los elementos equivalentes 1 8 en la Figura 6.29. El
coeficiente para todos los tubos tiene un valor de 0,018 0,02, en la
ecuacin de Colebrook.
Para la unin en Y, el coeficiente de prdidas de presin K3 se halla como
una prdida singular. Luego de conocer los valores de todos los coeficientes
de los componentes, se busca el coeficiente equivalente del sistema tras
sucesivas sumas de resistencias en serie y en paralelo como se explica en el
esquema de la Figura 6.30.
280
Los coeficientes de carga obtenidos son: Keq1 = 1309 bars2/m6, Keq2 = 1309
bars2/m6, Keq3 = 327 bars2/m6 y Kres = 1334 bars2/m6. En realidad estos
coeficientes no se necesitan sino para hallar un esquema equivalente ms
sencillo del circuito y en el programa de modelado 1D no se introducen, a
menos que en lugar de poner tuberas en el circuito reducido se pusieran
prdidas discretas (en tal caso los coeficientes hallados deberan dividirse por
la densidad y el inverso del cuadrado del rea: KDL = K/(/A2)).
Para asegurar la correspondencia del modelo reducido y de los componentes
equivalentes que lo integran, antes de componer el modelo en el programa 1D
se han procesado las caractersticas de carga del sistema y de presin de las
bombas en un proceso manual. En la Figura 6.31 se ven dichas caractersticas
para los subcircuitos del motor y del enfriador de aire.
Figura 6.32. Presentacin del modelo del circuito de refrigeracin del motor de la locomotora.
282
Los ramales hacia los radiadores tanto del motor como del enfriador de aire
tienen prcticamente la misma configuracin, excepto por las conexiones
hacia el compresor de aire y hacia el enfriador de aceite. A la salida de los
radiadores los flujos tanto de refrigerante como de aire confluyen y se dirigen
hacia la succin de las bombas y de los ventiladores, respectivamente.
En el modelo, en la manera como se obtienen las prdidas de presin, se logra
que los caudales que fluyen sean los adecuados y no hay mucha preocupacin
en relacin a la reparticin concentrada de las prdidas de presin, pues el
circuito se modela con un circuito equivalente de menor nmero de
accesorios. Entre los elementos equivalentes se han introducido prdidas
discretas de presin, que llevan la misma informacin hidrulica que un
bloque de motor o un intercambiador de calor. Es decir, en una conexin de
resistencias en serie en la que importa de forma neta la reparticin de
caudales, el mtodo es vlido. No lo es, si se requiere rigor en el anlisis de
las presiones en todos los componentes.
En el modelo simplificado se han conservado las tuberas de mayor longitud.
Las ms cortas y algunos elementos de unin y codos se han reemplazado por
elementos de prdidas concentradas como prdidas discretas. Las prdidas
discretas de presin se configuran con un rea de seccin de paso y los
coeficientes de prdidas de presin bajo flujo directo y reverso. Es comn la
utilizacin de estos elementos en los circuitos hidrulicos a manera de
potencimetros, para ajustar los caudales a los requeridos. Tambin pueden
usarse para ajustar las presiones en los sistemas y determinar las
caractersticas de componentes ausentes que deban incorporarse al sistema
fsico.
Las uniones en Y se configuran a partir de la geometra de los elementos
reales. Se les configuran los dimetros de los terminales, el ngulo entre los
ramales menores y los coeficientes de prdidas de presin para las
direcciones posibles y en funcin del tipo de combinacin en la que la unin
opera (dividiendo o combinando flujos, etc.). Para definir los coeficientes de
prdidas se utilizan las funciones de prdidas de presin de Miller [6.24].
Las tuberas del sistema son tuberas metlicas de acero laminado, con
rugosidad del orden de 0,025 mm. Al configurar su submodelo en el
programa, se entran el dimetro y su longitud, la opcin de clculo de la
283
284
los flujos msicos de refrigerante y de aire por unidad de rea frontal (figuras
5.28). La caracterstica que se introduce corresponde a condiciones estndar,
ya el programa automticamente se acopla a la variacin de las condiciones
segn evolucionen durante la simulacin.
Otros elementos presentes en el modelo de la Figura 6.32 son los
controladores y los medidores, a travs de los cuales se le pasa informacin a
los elementos principales y que se pueden configurar por programacin
mediante cdigo. Mediante controladores en el modelo se fijan los valores de
los calores a intercambiar en las fuentes trmicas y el estado y las
revoluciones de las bombas y los ventiladores. Se han colocado medidores
para conocer las temperaturas (variable de inters) a la salida del enfriador de
aire, del motor y de la bomba.
Para propsitos de ilustracin en la Figura 6.33 se destaca la conexin de los
dos ventiladores en paralelo en el sistema, el componente 17 simula la
atmsfera arriba de la locomotora; a este componente slo se le define el tipo
de fluido con el que opera, la presin y la temperatura. En la Figura 6.34 se
muestra el cuadro de dilogo de los parmetros que se ingresan del
ventilador.
285
286
287
Las grficas correspondientes para uno de los radiadores del enfriador de aire
se presentan en la Figura 6.38. Por los resultados del modelo se puede ver que
las temperaturas de los fluidos se mantienen en valores aceptables, segn las
condiciones establecidas en la tarea tcnica.
288
Q Q = m& m& + C p C p + T T
Si se reemplazan las desviaciones de flujo, calor y temperaturas del modelo
por las tolerancias de las variables en la ltima ecuacin, se obtendra la
contribucin discreta de las desviaciones de las variables a la desviacin total
del flujo de calor predicho, no considerando las interacciones. Lo correcto es
pensar que el modelo trabaja con unas tolerancias. Si se asumen las
desviaciones, por ejemplo desviaciones para las variables: m& = 3% ;
C p = 0, 012 ; T = 1,5 C , entonces la desviacin del flujo de calor predicho
puede alcanzar el valor de Q& = 214, 45kW .
289
de
290
291
292
Unid.
Elemento
Tipo
Funcin
Medidores de
caudal
Medidor de
caudal
13
Transductores
de presin
piezo-resistivos
13
Termopares
apantallados
Siemens
MAG 1100
DN80
Siemens
MAG 1100
DN25
Danfoss y
Cerobar.
rango de 010 bar
Tipo K
21
Termopares
apantallados
Tipo K
Sistema de
adquisicin de
datos de 30
canales
Registrador de
informacin y
osciloscopio de
16 canales
Fuente de
alimentacin
Columnas de
agua
Anemmetro
IO Tech
Personal
Daq/55
1
2
Medicin presin
refrigerante.
Medicin de temperatura
del
refrigerante
en
diferentes puntos del
sistema.
Medicin de temperatura
de aire en los radiadores,
ducto de ventilacin y
ambiente.
Registro de variables
medidas.
del
de
variables
Alimentacin
de
instrumentacin
Medicin de cada de
presin del aire.
Medicin de velocidad
del aire.
293
294
295
61,9
54,8
28,46
14,43
56,78
20,93
52,8
41,4
9,62
16,9
50
17,3
56,4
56,6
1,856
55,6
57,9
56,92
296
Ten_refr_AC = 52,8 C.
297
711,74
370
461,13
3,27
Balance
Experimental
trmico
radiador
del
Simulado
CFD
16,55
20,93
15,46
11
17,3
12,61
Las diferencias entre los caudales de aire medidos y los obtenidos por balance
trmico son bastante grandes, con lo que la evaluacin de los resultados del
modelo se hacen con los clculos derivados del balance trmico. Aunque
tambin puede argumentarse que el calor calculado con las temperaturas y
298
kW
-1423,4
-740
-264
2427,4
58,6
30,5
10,9
100
+1161
+802
+461,13
3,27
47,8
33,0
19,0
0,1
299
Simul.
64,26
56,37
7,89
14,4
60,58
46,18
29,49
Exper. Diferencia
2,36
61,9
54,8
1,57
7,1
0,79
0
14,4
56,78
3,8
42,38
3,8
28,46
1,03
16,89 16,55
0,34
797
711,44
85,56
% Error
3,81
2,86
11,13
0
5,63
8,97
3,62
2,05
12,04
Hay un calor evacuado por tuberas y dems que el modelo puede no estar
estimando bien. Aparte de esto estn las desviaciones en las magnitudes
medidas y simuladas (caudales y temperaturas).
Tabla 6.12. Comparacin de los valores medidos y los simulados en los
radiadores del aftercooler
Simul.
TEM
Entrada T (C) 53,69
Salida T (C)
41,42
Refrigerante
Enfriamiento (C) 11,69
Entrada T (C)
16,9
Salida T (C)
46,5
Aire
Calentamiento (C) 29,6
10,5
Flujo de refrigerante (kg/s)
Exper. Diferencia
% Error
52,8
41,4
11,4
16,9
50
33,1
9,62
0,9
0,02
0,29
0
3,5
3,5
0,88
0,92
0,048
2,54
0
7,0
3,02
9,14
13,44
11
2,44
22,18
408
370
38,16
10,31
300
Discusin. Las temperaturas del refrigerante al paso por las radiadores del
motor, predichas por el modelo, son ligeramente superiores que las
experimentales en trminos absolutos (C), con un error porcentual de
prediccin de temperatura medio de entrada y salida, del 3,3 % y un error en
el gradiente del refrigerante con respecto al experimental del 11 % (el
enfriamiento predicho por el modelo de diseo es mayor que el real
obtenido). La temperatura predicha para el aire tambin es mayor, con un
error del 3,8 %. Los flujos de refrigerante y de aire se predicen con una
diferencia menor del 1%. El calor disipado predicho es mayor en un 12 %
comparado con el experimental. Para la parte del enfriador de aire el anlisis
es similar, con menor error en la prediccin de las temperaturas, pero mayor
error en la prediccin de los flujos de refrigerante y de aire. El modelo
tambin predice un calor disipado superior en cerca de un 10 %.
301
cercanos para los flujos del refrigerante y del aire. La secuencia del proceso
se puede repetir las veces que sea necesario para optimizar la seleccin de los
radiadores y el diseo del sistema de ventilacin, pero antes debe pensarse en
realizar una evaluacin ms detallada de las tolerancias dentro de las cuales
se pueden realizar predicciones y mediciones en el sistema de refrigeracin,
pues son varias las fuentes de incertidumbre a pesar de estar en la direccin
adecuada.
El principal objetivo de cualquier proceso de diseo de ingeniera debe ser el
la definicin de un producto con las caractersticas funcionales definidas bajo
consideraciones claras de las incertidumbres y tolerancias involucradas. En el
diseo del sistema de refrigeracin se han realizado varias suposiciones desde
el inicio (partiendo de la adopcin de un valor de calor a disipar, que de
entrada no se conoce con mucha aproximacin). La acumulacin de estas
asunciones introduce variaciones en el diseo. Es por esto necesario entender
dnde se concentran las inexactitudes y tratar de eliminar todas las fuentes de
error posibles o, al menos, conocer con mayor certeza esas fuentes de error,
escogiendo las configuraciones en las que estos errores sean menos
probables. Para el caso concreto de este proyecto las principales
incertidumbres, aparte de las debidas al diseo constructivo de los
componentes, estn asociadas a la distribucin del flujo msico del aire a su
paso por las secciones de los radiadores. En esto las tcnicas experimentales
empleadas deben refinarse ms, pero adems debe contarse con las
especificaciones ms detalladas de los componentes fsicos. Como anotan
Taylor et al. [6.29], la cuantificacin de las incertidumbres en los problemas
de ingeniera trmica se ve dificultada porque los proveedores de los
componentes como bombas, ventiladores y accesorios, no proveen
informacin relacionada con la incertidumbre de sus caractersticas, y porque
tambin la informacin de ingeniera que reportan est basada generalmente
en criterios de semejanza y leyes de escalado. En la grfica de la figura 6.47
se ilustra la influencia que slo las tolerancias en las caractersticas integrales
del sistema de ventilacin, asumindolas de 5 % tanto en el flujo
volumtrico como en la presin, podran tener sobre el volumen de aire
conducido a travs de los radiadores.
Como se ve de la grfica, el caudal de aire puede variar de 55 a 62 m3/s, por
la sola variacin de las caractersticas constructivas de los componentes.
302
6.5. Conclusiones
En este captulo se ha seguido un procedimiento de diseo o rediseo de
sistemas de refrigeracin de motores descrito en el captulo cuarto. Conocidas
las especificaciones bsicas del motor y las particularidades de la aplicacin,
se ha estimado el calor a disipar, se han calculado los radiadores, se han
seleccionado los ventiladores, se ha perfilado un trazado de los circuitos del
sistema y se ha realizado su modelado para simular y comprobar que tanto el
circuito hidrulico ensamblado como el de ventilacin respondan a las
necesidades presentadas en la tarea tcnica. Se ha dedicado especial atencin
a ilustrar los resultados del modelado del flujo de aire a su paso por el
recorrido del sistema de ventilacin (sala de motor, radiadores, remanso,
6.6. Bibliografa
303
6.6. Bibliografa
[6.1]
[6.2]
[6.3]
[6.4]
[6.5]
[6.6]
[6.7]
[6.8]
[6.9]
[6.10]
[6.11]
[6.12]
[6.13]
304
Hill, N., Kollamthodi, S., Cross, S. Rail Diesel study, WP2 final
report: Technical and operational measures to improve the
emissions performance of Diesel rail, diciembre de 2005.
Page, R. W., Hnatczuk, W. y Kozierowski, J. Thermal
management for the 21st century improved thermal control and
fuel economy in an army medium tactical vehicle. SAE paper
2005-01-2068.
Sakai, T., Ishiguro, S., Sudoh, Y., Raab, G. y Hager, J. The
optimum design of engine cooling system by computer
simulation. SAE International Truck and Bus Meeting, Seattle,
1994.
Puntigam, W., Hrmann, T., Bernasch, J., Schierl, K., Wiesler, B.
y Hager J. Thermisches Management im Fahrzeug durch
Kopplung
unterschiedlicher
Simulationsmodelle
Wrmemangement des Kraftfahrzeuges IV, Essen, expert verlag,
2004 (consultado en http://www.vif.tugraz.at en diciembre de
2008).
PowerFlow users manual. Version 4.1.
Khlsysteme
of
the
company
Voith.
http://www.voithturbo.com/kuehlsysteme-produkteuehlsysteme.htm (consultada en enero de 2008).
Kays, W. M.; London, A. L. Compact heat exchangers. New
York, 1984.
Achachia, A., Cowell, T. A. Heat transfer and pressure drop
characteristics of flat tube and louvered plate fin surfaces.
Experimental Thermal and Fluid Science, vol. 1, pp. 147-157,
1988.
Incropera, F.P. and DeWitt, D.P. Fundamentals of heat and mass
transfer. New York, 1996.
Namai, K., Muramoto, H., Mochizuki, S. Flow visualization in
the louvered fin heat exchanger. SAE paper 980055.
Tafti, D.K. and Cui, J. Computations of flow and heat transfer in a
three dimensional multilouvered fin geometry. International
Journal of Heat and Fluid Flow, 2002.
Webb, R. L.,Trauger, P. Flow Structure in the Louvered Fin Heat
Exchanger Geometry. Experimental Thermal and Fluid Science,
6.6. Bibliografa
[6.14]
[5.15]
[6.16]
[6.17]
[6.18]
[6.19]
[6.20]
[6.21]
[6.22]
[6.23]
[6.24]
[6.25]
305
306
Captulo 7
Conclusiones y trabajos futuros
7.1. Conclusiones ......................................................................................... 309
7.2. Trabajos futuros..................................................................................... 316
308
7.1. Conclusiones
309
7.1. Conclusiones
Como fruto del trabajo realizado, se pueden extraer algunas conclusiones
principales, derivadas directamente de las exposiciones hechas a lo largo de
este escrito, y que se exponen a continuacin, segn el orden marcado por los
diferentes captulos. Relacionadas con el estudio del estado del arte del
modelado y las mediciones experimentales, de manera amplia, puede
afirmarse que:
310
7.1. Conclusiones
311
312
Extender las capacidades del modelo trmico, para predecir los flujos
de calor en el motor hacia los medios correspondientes a las
condiciones de contorno, en particular hacia el refrigerante y el aceite,
bajo condiciones de operacin estacionarias y transitorias.
Predecir las temperaturas y los flujos de calor entre los nodos del
modelo del motor, y obtener el mapa trmico de ste, junto con una
correlacin para cuantificar el calor disipado a travs de las paredes de
7.1. Conclusiones
313
314
7.1. Conclusiones
315
316
317
318
Apndice A
A.1. Acerca de AMESim
AMESim (sigla de Advanced Modeling Environment for performing
Simulations) es un programa multipuerto, en el que se intercambia
informacin de flujo de energa entre componentes a travs de los nodos de
conexin. AMESim posee unas libreras primarias o estndar y unas libreras
especializadas, como se puede apreciar en el cuadro de presentacin de
Amesim de la Figura A.1. Se tiene la posibilidad de crear modelos propios y
supercomponentes.
Apndice A
320
321
Anlisis de los resultados. En esta fase el usuario puede hacer anlisis de los
resultados mediante el procesamiento mediante anlisis de respuesta
dinmica, transformada rpida de Fourier, anlisis del lugar de las races,
posibilidad de realizar diseo de experimentos y optimizacin del modelo,
entre otras posibilidades. Como caracterstica importante vale la pena citar la
facilidad de realizar operaciones con los resultados y de editar los resultados
de esas operaciones.
Del grupo de libreras Thermal Fluid y Cooling System de AMESim, los
submodelos sobre los que se fundamenta el ensamble de un sistema de
refrigeracin convencional moderno son los mostrados en la Figura A.2.
Apndice A
322
Parmetros a definir
Archivo con velocidad de
referencia del motor y
velocidad del vehculo en
funcin de la relacin de
transmisin
Archivo con la posicin del
pedal de aceleracin en
funcin del tiempo
Archivo con la informacin
del par del motor en funcin
de la posicin del pedal y de
las revoluciones del motor.
Nombre del archivo con la
velocidad del vehculo en
funcin del tiempo
Nombre del archivo con las
relaciones de la caja de
velocidades en funcin del
tiempo
Nombre del archivo con la
inclinacin de la calzada en
funcin del tiempo
rea equivalente
Volumen de refrigerante
involucrado
ndice del fluido
Archivo con el calor
intercambiado entre el motor
y el refrigerante, en funcin
de la potencia del motor y
las revoluciones
323
Volumen de refrigerante
Masa del material (masa
trmicamente sensible)
rea equivalente
rea de transferencia de
calor de la camisa del motor
Coeficiente de transferencia
de calor
ndice del fluido
ndice del material del motor
Relacin de transmisin
entre la bomba y el motor
Volumen de refrigerante
asociado con la bomba
Eficiencia total de la bomba
Apndice A
324
expansin.
Constante de tiempo
ndice del fluido
Nombre del archivo dp =
f(Qv, velocidad del impeler)
325
Apndice A
326
327
Apndice A
328
A.3. Bibliografa
[A.1]
[A.2]
[A.3]
A.3. Bibliografa
329
[A.4]
[A.5]
[A.6]
[A.7]