Tesis Transporte Calor Comsol

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 84

MODELACION DE FLUJO LAMINAR Y TRANSFERENCIA DE

CALOR EN HAZ DE TUBOS



MEMORIA PARA OPTAR AL TTULO DE INGENIERO CIVIL MECNICO




CARLOS FRANCISCO ESPINOSA BARRIOS



PROFESOR GUA:
RAMN FREDERICK GONZALEZ


MIEMBROS DE LA COMISIN:
WILLIAMS CALDERN MUOZ
ALEJANDRO ORTIZ BERNARDIN



SANTIAGO DE CHILE
JUNIO 2012



UNIVERSIDAD DE CHILE
FACULTAD DE CIENCIAS FSICAS Y MATEMTICAS
DEPARTAMENTO DE INGENIERA MECNICA
La conveccin es una de las 3 formas de transferencia de calor y se produce por
intermedio de fluidos que transportan el calor hacia otras zonas de distinta temperatura. Uno
de los casos de mayor importancia tecnolgica es la conveccin en el exterior de bancos de
tubos con flujo cruzado, debido a su importancia en procesos y aplicaciones industriales.
El coeficiente convectivo representado habitualmente por h, mide la influencia de las
propiedades de un fluido, de la superficie y su rgimen de flujo cuando se produce el fenmeno
de conveccin. Dado el alto costo de experimentos para predecir estos coeficientes, es una
alternativa viable la modelacin de este coeficiente mediante uso de Software de elementos
finitos y transferencia de calor.
En la presente Memoria se realiz la modelacin de flujo laminar y transferencia de calor en un
haz de tubos. El software de elementos finitos utilizado es COMSOL MULTIPHYSICS V3.5 que
est disponible en el departamento.
Dada la complejidad que presentan estos equipos en cuanto a su configuracin, se utilizaron
computadores potenciados para considerar grandes arreglos de tubos y de esta misma forma se
analiz el comportamiento e influencia entre tubos.
Se consider para este estudio un modelo de flujo laminar y transferencia de calor con bajos
nmeros de Reynolds, adems que el estudio se program en un rgimen transiente para ver la
evolucin de estos parmetros en el tiempo. Para una mayor rapidez de procesamiento se
restringi al caso de dos dimensiones.
El modelo bidimensional dependiente del tiempo se construy en dos versiones:
a) Para un tubo nico
b) Para un haz de ocho tubos, en cuatro configuraciones geomtricas.
En el rango de nmero de Reynolds entre 900 y 2400, los coeficientes convectivos para tubo
nico presentaron diferencias del orden del 10% con respecto a los resultados empricos. Con
respecto al modelo de haz de tubos, los errores fueron mayores, fluctuando entre un 5 y un
20%.
INDICE

1. Introduccin ........................................................................................................................................... 5
1.1. Antecedentes generales ................................................................................................................. 5
1.2. Motivacin ..................................................................................................................................... 7
1.3. Objetivos ........................................................................................................................................ 8
1.4. Alcances .......................................................................................................................................... 9
2. Antecedentes .......................................................................................................................................10
2.1 Consideraciones iniciales .............................................................................................................10
2.2 Flujos Contra y Co-corrientes .......................................................................................................14
2.3 Geometras Evaporadores ............................................................................................................14
2.3.1 Disposicin del arreglo .........................................................................................................14
2.3.2 Parmetros geomtricos ......................................................................................................15
2.4 Ecuaciones gobernantes ..............................................................................................................15
2.4.1 Ecuaciones principales .........................................................................................................15
2.4.2 Correlaciones conocidas de transferencia de calor .............................................................17
2.4.3 Cada de presin ...................................................................................................................23
2.5 Mtodo de elementos finitos .......................................................................................................25
2.6 Ecuaciones del Software ..............................................................................................................26
2.6.1 Weakly compresible Navier Stokes ..........................................................................................26
2.6.2 General Heat Transfer ..............................................................................................................27
3. Metodologa .........................................................................................................................................28
3.1 Geometra y condiciones de borde ..............................................................................................29
3.2 Parmetros de Ecuacin ...............................................................................................................31
3.3 Criterios de convergencia .............................................................................................................34
4. Anlisis de modelos ..............................................................................................................................36
4.1 Modelo bsico ..............................................................................................................................36
4.1.1 Sensibilidad de la simulacin al mallado (nmero de elementos de la malla) ....................36
4.1.2 Optimizacin del tiempo de simulacin ...............................................................................39
4.1.3 Validacin de los coeficientes convectivos ..........................................................................41
4.2 Pruebas de mejora Modelo Bsico ..............................................................................................45
4.2.1 Cambio en el mallado ...........................................................................................................45
4.2.2 Cambios geomtricos ...........................................................................................................48
4

4.2.3 Campos de temperatura y velocidad ...................................................................................51
4.2.4 Presin entorno al cilindro ...................................................................................................59
5. Modelo Haz de tubos ...........................................................................................................................60
5.1 Condiciones iniciales ....................................................................................................................60
5.2 Modelo St/D=Sl/D=1.5 .................................................................................................................61
5.3 Modelo St/D=2; Sl/D=1.5 .............................................................................................................69
5.4 Modelo St/D=1.5; Sl/D=2 .............................................................................................................72
5.5 Modelo St/D=Sl/D=2 ....................................................................................................................75
5.6 Anlisis Interaccin de tubos .......................................................................................................78
6. Conclusiones.........................................................................................................................................81
Bibliografa ...................................................................................................................................................84

























5

1. Introduccin
1.1. Antecedentes generales

La refrigeracin mecnica es un proceso en el cual se reduce la temperatura de un espacio
determinado. Este proceso se basa principalmente en las leyes de la termodinmica y depende
de los fenmenos de transmisin de calor. Su propsito es el enfriamiento de ambientes
habitables para la comodidad de sus ocupantes, as como obtener condiciones apropiadas para
la conservacin de sustancias.
Para la industria frutcola el almacenamiento refrigerado es un tema fundamental debido a
que evita el crecimiento de bacterias e impide reacciones qumicas no deseadas que se dan a
temperatura ambiente y por ende es un proceso obligatorio dentro de la conservacin de las
frutas.
El ciclo de refrigeracin ms usado en la industria es el ciclo de refrigeracin por compresin
de vapor, cuyo ciclo termodinmico se realiza mediante cuatro componentes principales: un
compresor, un condensador, una vlvula de expansin y un evaporador.
En este ciclo se logra la evaporacin del lquido refrigerante dentro de un intercambiador de
calor denominado evaporador. Para evaporarse, este lquido absorbe energa trmica sensible
del ambiente, y por ende lo enfra. Luego de este intercambio trmico, el vapor pasa por un
compresor que aumenta la presin y la temperatura del refrigerante para lograr su
condensacin dentro de otro intercambiador de calor (condensador). Finalmente el refrigerante
ya en estado lquido, pasa por una vlvula de expansin que le restituye la presin apropiada
para su ingreso al evaporador, permitiendo su evaporacin para un nuevo ciclo.








El evaporador de refrigeracin es un intercambiador de calor donde el refrigerante se
evapora al interior de tubos, y el fluido a enfriar pasa por exterior. Consideraremos en este
estudio que el fluido exterior es aire.
Figura 1.1: Ciclo de Refrigeracin
6

Principalmente, los evaporadores son diseados en haces o bancos de tubos, cuyos
arreglos son de dos tipos de geometra que corresponden a la rectangular y la triangular. Esto
corresponde a la forma en cmo se alinean los tubos. En la figura podemos observar una forma
tpica de los evaporadores. El flujo de aire pasa alrededor de los tubos y el refrigerante que pasa
por dentro de los tubos se evapora extrayendo energa del aire, que por ende se enfra. Estos
evaporadores tienen aletas por el lado del aire, las cuales abarcan todos los tubos (Figura 1.2).
El nmero de filas verticales de tubos en el arreglo es usualmente 4.











En esta memoria se construir un modelo de la geometra del evaporador, sin considerar las
aletas. El propsito del estudio es valores del coeficiente convectivo por el lado del aire.

Figura 1.2: Formas tpicas de evaporadores
7

1.2. Motivacin

En la presente memoria se pretende crear un modelo bidimensional de transferencia de
calor y flujo laminar en un banco de tubos. Considerando que ya se realiz una memoria que
antecede este estudio [1], una de las principales motivaciones es justamente profundizar el
anlisis, y abordar casos ms complejos. Dado que un modelo de un banco de tubos con aletas
requiere de un estudio tridimensional, que resulta imposible por la tecnologa con la cual se
cuenta actualmente en el departamento, este estudio se basa en el anlisis de un banco de
tubos sin aletas, el cual puede realizarse en geometra bidimensional.
La idea principal es representar bancos de tubos en forma realista, es decir con dimensiones
y espaciados que se dan en los actuales diseos de evaporadores. De esta forma, se puede
describir la influencia de parmetros geomtricos sobre el intercambio calrico. En la memoria
anterior, no fue posible realizar un anlisis con dimensiones reales, por que lo l estudio no
pudo arrojar los resultados que se buscaban.
El desarrollo de este proyecto involucra por supuesto dos de las principales disciplinas de la
Ingeniera Mecnica, como son la transferencia de calor y la mecnica de fluidos. Todo esto
sustentado por el uso de un Software de Elementos Finitos como lo es el Comsol Multiphysics
v3.5a que permite realizar el anlisis de los fenmenos fsicos.
El uso de este tipo de software se hace cada vez ms importante para el desarrollo
profesional. Son plataformas para realizar simulaciones donde se pueden analizar una enorme
variedad de escenarios y geometras a un bajo costo. Por esto mismo, es una alternativa viable a
todo tipo de experimentos.
El estudio de evaporadores es crucial para todas las aplicaciones de la industria frutcola
entre otras. Es conocido que las frutas se mantienen a temperaturas bajas para su conservacin,
por lo que estudios de esta ndole permiten desarrollar la eficiencia energtica, y ms an
permiten mejorar diseos dentro de los mismos equipos. Sin embargo, el estudio no se limitar
a las aplicaciones de la industria frutcola, sino que a diversas situaciones donde los fluidos
ambientes pueden ser tanto aire como agua, al igual que los lquidos refrigerantes.
El fenmeno que se aborda aqu es un caso clsico de estudio en mecnica de fluidos y
transferencia de calor. El flujo alrededor de un cilindro nico o a travs de un banco de tubos ha
sido descrito la mayor parte de las veces de manera experimental. No se dispone de gran
cantidad de informacin de origen numrico para este problema. Una de las posibles razones es
que el fenmeno a estudiar presenta las caractersticas de capa lmite con posibilidad de
separacin y transicin a turbulencia. Los elevados gradientes de velocidad y temperatura que
se producen en el flujo alrededor de un cilindro son los problemas mayores en el estudio, como
tambin la naturaleza no permanente de los flujos resultantes.


8

1.3. Objetivos

En el presente tema se pretende crear un modelo para el flujo y la transferencia de calor
alrededor de bancos de tubos. Se plantea bajo los siguientes objetivos:
General
Construir un modelo para la determinacin de coeficientes convectivos en el exterior
de intercambiadores de calor de haz de tubos con flujo cruzado.
Especficos:
Construccin de un modelo 2D y transiente para el anlisis de un tubo nico para la
validacin del coeficiente convectivo.
Creacin de un modelo de simulacin 2D y transiente del flujo de fluidos y
transferencia de calor en un intercambiador de flujo cruzado.
Obtener resultados de transferencia de calor con los modelos realizados. En
particular, analizar la interaccin entre tubos en funcin de los espaciamientos entre
ellos.

9

1.4. Alcances

En el estudio se realizan las siguientes simplificaciones:
Uso de nmeros de Reynolds bajos, dentro del rango laminar, puesto que la turbulencia es
difcil de simular con el software actual. La construccin del modelo geomtrico del banco de
tubos consistir en un arreglo de 6 a 9 tubos, limitndose a la capacidad del computador, y se
analizar la interaccin del flujo entre los mismos tubos, por lo que la distancia entre tubos ser
pequea.
Para el estudio se considerar el aire como un flujo incompresible de propiedades
constantes. A pesar que el software a usar realiza el acoplamiento de los fenmenos de flujo de
fluidos y transferencia de calor bajo la opcin de fluidos compresibles, no se usar esta opcin
dado que las velocidades de aire a usar son muy bajas.
A diferencia del estudio anterior, la capacidad del ordenador disponible es mayor, luego ser
posible realizar los estudios pertinentes de las interacciones de los tubos, y adems de un
nmero de arreglos de tubos mayor.
La temperatura de la pared exterior de los tubos se asume constante, por lo que no se
considera el fenmeno de ebullicin del refrigerante como parte del estudio. A la regin de
anlisis ingresar un flujo de aire de velocidad constante y temperatura uniforme.
El evaporador de refrigeracin es uno de los campos de inters para intercambiadores de
bancos de tubos, pero no ser el nico. Tambin se puede considerar casos en que el tubo est
a mayor temperatura que el aire, y otros fluidos distintos de aire, ya que con aire la
transferencia de calor es muy baja. Se considera en estos casos que los resultados de
transferencia de calor son vlidos tanto para calentamiento como para enfriamiento de los
fluidos que se desarrollan alrededor del haz de tubos.
10

2. Antecedentes
2.1 Consideraciones iniciales

El modelo de capa lmite explica la resistencia de un cuerpo sobre el fluido, pero claramente
cambia dependiendo la geometra del cuerpo en cuestin. Incluso para Reynolds grandes, el
acoplamiento con la corriente exterior es muy complejo. La forma del cuerpo afecta
considerablemente a la distribucin de presiones en el flujo, y por ende tiene una influencia
notable sobre la capa lmite.
Como se observa en la figura 2.1, en la parte frontal el dficit de cantidad de movimiento en la
capa lmite no es un problema debido a que la disminucin de presin (y por ende aumento de
velocidad) genera un gradiente favorable de presiones que empuja el fluido en la capa lmite.
Sin embargo, como se observa en la parte posterior la situacin es totalmente distinta. La
presin aumenta en la direccin del flujo apareciendo un gradiente adverso de presiones que
puede llegar a frenar el flujo en la capa lmite, eventualmente invertir su direccin, este es el
fenmeno de separacin de capa lmite [2].


Figura 2.1: Esquema del flujo externo alrededor de un cilindro as como de la capa limite desarrollada sobre su superficie

Luego, aparecen recirculaciones en la parte posterior del tubo, y en algunos casos, la capa lmite
se separa generndose en la parte posterior una estela pulsatoria que deflecta la corriente
principal.
11


Se puede realizar por lo tanto una descripcin cualitativa de la fenomenologa del problema.
Cuando el nmero de Reynolds es muy pequeo (Re<1), las fuerzas viscosas controlan el flujo, y
el patrn de este alrededor del cilindro es similar al de la figura 2.3. Para valores de Reynolds
mayores (Re=20), el flujo es laminar aguas arriba y en un punto se produce la separacin. Aguas
abajo aparece una regin con dos remolinos estacionarios (Figura 2.2).

Figura 2.3: Flujo alrededor de cilindro para Re<1
Para valores de Reynolds sobre 100 (Figura 2.4), los remolinos se inestabilizan y comienzan a
vibrar irregularmente para acabar desprendindose alternativamente a ambos lados de forma
peridica. Son los denominados vrtices de Von Karman.




Cuando los valores de Reynolds son tales que Re=10
4
, el flujo se separa del cilindro para formar
una estela simtrica aguas abajo (Figura 2.6). Cuando Re= 10
6
el coeficiente de arrastre se hace
constante y la transicin ocurre por delante del punto de separacin (Figura 2.5), que se retrasa
debido al carcter turbulento del flujo.

.
Figura 2.6: Flujo alrededor de cilindro para Re=10
4


Como se puede observar, la ocurrencia de una transicin de la capa lmite depende
directamente del nmero de Reynolds, que al mismo tiempo influye fuertemente en la posicin
del punto de separacin. A medida que el Reynolds crece, el punto de separacin se retrasa
(Figura 2.7), el ancho de estela se reduce, al igual que la fuerza de arrastre [3].
Figura 2.4: Flujo alrededor de cilindro para Re=100
Figura 2.2: Flujo alrededor de cilindro para Re=20
Figura 2.5: Flujo alrededor de cilindro para Re=10
6

12



Figura 2.7: Efecto de la turbulencia sobre la separacin

El coeficiente de arrastre es una funcin del nmero de Reynolds, como se observa en la figura
2.8. Al aumentar el Reynolds, el coeficiente de arrastre disminuye debido a la presin. Del orden
de Re=2x10
5
, la capa pasa de ser laminar a turbulenta, por lo que se retrasa la separacin, y por
ende el ancho de la regin de la estela.


Figura 2.8: Coeficiente de arrastre para un cilindro circular en flujo cruzado y para esfera.





13


En la figura 2.9, se observa resultados experimentales para la variacin del nmero de Nusselt
en funcin del ngulo. Para Re<10
5
, Nusselt disminuye a medida que el ngulo aumenta por el
desarrollo de la capa lmite. Alcanza un mnimo a los 80, para luego cambiar su tendencia, y
aumentar debido a la formacin de los vrtices. Por otro lado, para Re> 10
5
, se alcanzan dos
mnimos, el primero debido a la formacin de la capa lmite y el segundo debido al desarrollo de
la capa lmite turbulenta. El aumento entre ambos mnimos se debe a la transicin de capa
lmite a turbulenta. Finalmente, lo importante es que partir de esta figura, se puede ver la
tendencia del coeficiente convectivo en funcin de los ngulos. [3]



Figura 2.9: Nmero de Nusselt en cilindro circular para flujo cruzado
14

2.2 Flujos Contra y Co-corrientes

Existen dos tipos de arreglos para un evaporador en cuanto a su flujo, que corresponden a
un arreglo en co-corriente y contracorriente. De la misma forma, se puede trabajar con flujo
paralelo o cruzado, esto corresponde al ngulo de ataque del fluido al tubo, y como se
mencion anteriormente se trabajar con este ltimo. Esto se debe principalmente a que se
obtiene un mejor rendimiento en cuanto a la transferencia de calor para este tipo de arreglo,
mientras que la virtud de un sistema en co-corriente reside en la mantencin de una
temperatura relativamente constante entre los fluidos. [4]
A la salida del evaporador, el lquido refrigerante se encuentra tericamente en un estado
vaporizado sobrecalentado, pero en algunos casos se encuentra parcialmente vaporizado. En los
casos donde el lquido est parcialmente vaporizado es de suma importancia que est en
estado gaseoso antes de ingresar al compresor para de esta forma asegurar la integridad del
equipo.

2.3 Geometras Evaporadores

2.3.1 Disposicin del arreglo
Las geometras de los evaporadores se componen por arreglos de tubos triangulares y/o
rectangulares. Los parmetros que varan en esta geometra son el dimetro externo de los
tubos, la separacin vertical y horizontal de los tubos y el largo de los tubos. En este caso el
estudio consiste en la modelacin del banco de tubos en dos dimensiones por lo que no ser
relevante la variable del largo de los tubos, solo se basar en disposiciones alineadas y
alternadas como se puede observar en la figura 2.10. [5]

Figura 2.10: Flujo laminar en dos tipos de arreglos geomtricos, (a) Rectangulares y (b) triangulares

15

2.3.2 Parmetros geomtricos
Las principales variables que se manejan en este problema bidimensional corresponden al
dimetro de tubo y la distancia entre los centros de los tubos. El cambio en dichas variables
provoca justamente la caracterizacin del flujo y el movimiento del fluido. Como se quiere
trabajar con un fluido en rgimen laminar, el nmero de Reynolds ser bajo.
Las bobinas de los evaporadores usan generalmente tubos circulares, sin embargo, existen
tubos rectangulares y ovalados que se usan para aplicaciones especiales. Para un amplio rango
de aplicaciones, existen tpicos tamaos de tubos, donde el dimetro externo puede tomar
diversos valores tales como 5/16, 3/8, , 5/8, y 1 pulgada (7.9, 9.5, 12.7, 15.9, 19.1 y 25.4
mm).
De la misma forma, como vimos en el punto 2.3.1, los tubos pueden formarse en arreglos en
lnea como arreglos triangulares, siendo este primero, el mtodo ms comn. El espacio entre
tubos puede variar entre 0.6 a 2.5 pulgadas (16 a 64 mm). [4]
Una forma de caracterizar estos parmetros geomtricos, es crear una variable a partir de
estos dos componentes, es decir, se trabajara con L/D, donde L corresponde a la distancia entre
los tubos, y D al dimetro de estos. Todo esto sujeto a un nmero de Reynolds fijo.

2.4 Ecuaciones gobernantes

2.4.1 Ecuaciones principales
Las ecuaciones que gobiernan este sistema son principalmente las que vienen de las leyes de la
termodinmica, que son las de movimiento, continuidad y energa para flujo laminar
incompresible bidimensional transiente:

Balance de calor para flujo alrededor de un tubo colocado en un ducto [1]


Donde:
: Calor transferido por el tubo a la corriente de aire

: Calor especifico


: Flujo Masico

: Temperatura de entrada
16

: Temperatura de salida

Coeficiente convectivo medio alrededor del tubo



Donde:
: Coeficiente convectivo


: Diferencia entre temperatura de pared y entrada de aire
: rea de transferencia de calor

(Manto de tubo)

Nmero de Reynolds


Donde:
: Dimetro exterior del cilindro
: Viscosidad dinmica


: Velocidad de la corriente libre



Nmero de Prandtl


Donde:
: Conductividad trmica



Nmero de Nusselt


17


2.4.2 Correlaciones conocidas de transferencia de calor

a) Caso tubo nico con flujo cruzado
Se pueden obtener distintas correlaciones para nmeros de Nusselt local, sin embargo desde el
punto de vista de los clculos de ingeniera, es preferible utilizar correlaciones en condiciones
promedios globales. [6]

Hilpert



Donde las constantes C y m se listan en la tabla 2.1:
Tabla 2.1: Constantes para la ecuacin de Hilpert, para cilindro circular de flujo cruzado.
Re C m
0,4 - 4 0,989 0,33
4 - 40 0,911 0,385
40 - 4000 0,683 0,466
4000 - 40000 0,193 0,618
40000 - 400000 0,027 0,805

Zhukauskas (1972)




Donde C y m estn dados por la tabla 2.2
Tabla 2.2: Constantes para la ecuacin de Zhukauskas, para cilindro circular de flujo cruzado.
Re C m
4 - 40 0.75 0,4
40 - 1000 0.51 0.5
1000 - 200000 0.26 0.6
200000 - 2000000 0.076 0.7

18

En nuestro caso

se aproxima a 1 debido a que se utilizan propiedades constantes del


fluido.

Por otro lado, Churchill y Bernstein propusieron una sola ecuacin que cubre todo el rango de
Reynolds, pero que se recomienda principalmente para

Churchill Bernstein [7]



Whitaker present una correlacin ms general que toma en cuenta las contribuciones al
coeficiente promedio de transferencia de calor de la regin no separada de la capa lmite, de la
regin de la estela alrededor del cilindro y de los efectos de la temperatura.

Whitaker [7]



Donde

representa la viscosidad de la temperatura de la superficie de la pared.


Estas correlaciones se utilizan principalmente para calcular un Nusselt promedio para la
simulacin de un tubo con flujo cruzado.




19


En la figura 2.11 se puede observar los valores del coeficiente convectivo para las distintas
correlaciones en funcin de la velocidad del fluido que en este caso es aire, para un dimetro
D=0,5 cm.

Figura 2.11: Coeficiente convectivo en funcin de la velocidad para las distintas correlaciones

b) Caso arreglo de tubos con flujo cruzado
Para el caso de un haz de tubos, se desea conocer el coeficiente promedio de transferencia de
calor para todo el arreglo de tubos, es as como Grimison obtuvo una correlacin para un flujo
de aire, de la forma:
Grimison (1937) [3] [6]




Con



0,00
2,00
4,00
6,00
8,00
10,00
12,00
14,00
16,00
18,00
0 0,2 0,4 0,6 0,8 1 1,2
h

[
W
/
m
2
.
K
]

Velocidad [m/s]
Hilpert
Zhukauskas
Churchill Berstein
Whitaker
20



Donde para arreglo en lnea


Y para arreglo escalonado



Con


Donde
Velocidad de la corriente antes de interactuar con los tubos [m/s]

Distancia vertical entre centro de tubos.

Distancia horizontal entre centro de tubos.


Donde C1 y m se presentan en la tabla 2.3
Tabla 2.3: Constantes para el flujo de aire sobre un banco de tubos de 10 o ms lneas



St/D


1.25 1.5 2 3
Sl/D C1 m C1 m C1 m C1 m
Alineado

1.25 0.348 0.592 0.275 0.608 0.1 0.704 0.0633 0.752
1.5 0.367 0.586 0.25 0.62 0.101 0.702 0.0678 0.744
2 0.418 0.57 0.299 0.602 0.229 0.632 0.198 0.648
3 0.29 0.601 0.357 0.584 0.374 0.581 0.286 0.608

Escalonado

0.6 - - - - - - 0.213 0.636
0.9 - - - - 0.446 0.571 0.401 0.581
1 - - 0.497 0.558 - - - -
1.125 - - - - 0.478 0.565 0.518 0.56
1.25 0.518 0.556 0.505 0.554 0.519 0.556 0.522 0.562
1.5 0.451 0.568 0.46 0.562 0.452 0.568 0.488 0.568
2 0.404 0.572 0.416 0.568 0.482 0.556 0.449 0.57
3 0.31 0.592 0.356 0.58 0.44 0.562 0.428 0.574
21

En un arreglo en lnea, la velocidad mxima se produce en el pasillo de alta velocidad, es decir
en el plano trasversal A1 de la figura 2.10a, donde el fluido no choca con el tubo. Por otro lado,
en un arreglo escalonado, la mxima velocidad se puede dar en el plano transversal A1, o en su
defecto en el plano diagonal A2. Numricamente, ocurrir en A2 si las filas estn espaciadas de
modo que


Si Vmax no ocurre en A2, entonces se calcula nuevamente Vmax para arreglo en lnea.
Este resultado se puede extender a otros fluidos mediante la insercin del factor





De igual forma, si el arreglo es menor a 10 filas de tubos

se aplica un factor de
correccin tal que


Donde

esta dado por la tabla 2.4


Tabla 2.4: Factor de correlacin

para


N
l
1 2 3 4 5 6 7 8 9
Alineado 0.64 0.8 0.87 0.9 0.92 0.94 0.96 0.98 0.99
Escalonado 0.68 0.78 0.83 0.89 0.92 0.95 0.97 0.98 0.99

Zhukauskas (1988) [3] [6]

Zhukauskas propuso una correlacin de la forma:




Donde las constantes C y m se observan en la siguiente tabla:
Tabla 2.5: Constantes para ecuacin de banco de tubos con flujo cruzado
Configuracion Re max C m
Alineado 10-10
2
0.8 0.4
Escalonado 10-10
2
0.9 0.4
Alineado 10
2
- 10
3
Aproxima como cilindro unico
22

Escalonado 10
2
- 10
3
(Aislado)
Alineado 10
3
- 2x10
5
0.27 0.63
(St/Sl<0.7)

Escalonado 10
3
- 2x10
5
0.35(St/Sl)^(1/5) 0.6
(St/Sl<2)

Escalonado 10
3
- 2x10
5
0.4 0.6
(St/Sl>2)

Alineado
2x10
5
-
2x10
6
0.021 0.84
Escalonado
2x10
5
-
2x10
6
0.022 0.84

Si

se aplica un factor de correlacion tal que:


Donde

se encuentra en la table 2.6


Tabla 2.6: Factor de correlacin

para


N
l
1 2 3 4 5 7 10 13 16
Alineado 0.70 0.8 0.86 0.9 0.92 0.95 0.97 0.98 0.9
Escalonado 0.64 0.76 0.84 0.89 0.92 0.95 0.97 0.98 0.9

Para Re<1000,


Tabla 2.7: Factor de correccin para tubos escalonados con Re<1000 y Nl<20
N
L
1 2 3 4 5 7 10 13 16
C
2
0,83 0,88 0,91 0,94 0,95 0,97 0,98 0,99 1

Y la constante C
3
, se obtiene a partir de la tabla 2.8.
Tabla 2.8: Constantes para la ecuacin de Zhukauskas para Re<1000 para tubos escalonados
Re
d
C
3
m
1-40
0,75 0,4
40-1000
0,51 0,5
10
3
-2x10
5

0,26 0,6
2x10
5
-10
6

0,076 0,7


23

2.4.3 Cada de presin
El modelo de presin de Zhukauskas esta expresado por la siguiente ecuacin:


Donde f representa el factor de friccin y x el factor de correccin que se estn representados
grficamente para tubos escalonados en la figura:


Figura 2.12: Factor de friccin para modelo Zhukauskas

24


Figura 2.13: Factor de correccin para modelo Zhukauskas


25

2.5 Mtodo de elementos finitos

El mtodo de los elementos finitos (MEF) es un mtodo numrico general para la
aproximacin de soluciones de ecuaciones diferenciales parciales muy utilizado en diversos
problemas de ingeniera y fsica.
El MEF est pensado para ser usado en computadoras y se usa en el diseo y mejora de
productos y aplicaciones industriales, as como en la simulacin de sistemas fsicos y biolgicos
complejos. La variedad de problemas a los que puede aplicarse ha crecido enormemente, y
permite resolver ecuaciones diferenciales asociadas a un problema fsico sobre geometras
complicadas.
El MEF permite obtener una solucin numrica aproximada sobre un cuerpo, estructura o
dominio (medio continuo) dividindolo en un nmero elevado de subdominios denominados
elementos finitos. El conjunto de elementos finitos forma una particin del dominio tambin
denominada discretizacin. Dentro de cada elemento se distinguen una serie de puntos
representativos llamados nodos. Dos nodos son adyacentes si pertenecen al mismo elemento
finito; adems, un nodo sobre la frontera de un elemento finito puede pertenecer a varios
elementos. El conjunto de nodos considerando sus relaciones de adyacencia se denomina
malla.
Los clculos se realizan en esta malla, y a partir de estas relaciones de conectividad se calculan
el valor de un conjunto de variables incgnitas en cada nodo. Este conjunto se puede definir
como un sistema de ecuaciones lineales, siendo el nmero de ecuaciones proporcional al
nmero de nodos.
En nuestro caso, el programa de generador de mallas est incluido en el software que resuelve
estas ecuaciones, y se trata del programa Comsol Multiphysics v3.5.
Una importante propiedad del mtodo es la convergencia; si se consideran particiones de
elementos finitos sucesivamente ms finas, la solucin numrica calculada converge
rpidamente hacia la solucin exacta del sistema de ecuaciones. Por otro lado, si las relaciones
entre cada nodo no son continuas, o superan un cierto nmero considerado como la tolerancia
del modelo, el modelo tiende a diverger.




26

2.6 Ecuaciones del Software

El software escogido para la simulacin por elementos finitos es el Comsol Multiphysics v3.5.
Las ecuaciones que utiliza para el estudio corresponden a las ecuaciones del mdulo de
transferencia de calor, para una interaccin trmica del fluido en un anlisis transiente.

2.6.1 Weakly compresible Navier Stokes


Donde



Aplicacin a nuestro estudio:
No se considera dilatacin viscosa



No existen fuerzas externas =0

Por lo tanto la ecuacin queda de la siguiente forma:





27


2.6.2 General Heat Transfer


Donde




Aplicacin a nuestro estudio:
No se considera coeficiente de absorcin



No existe una fuente de calor externa

La ecuacin de transferencia de calor queda de la siguiente forma:






28

3. Metodologa

Para realizar la modelacin se utilizar el programa Comsol Multiphysics versin 3.5 que
permite realizar el anlisis mediante el uso de elementos finitos. Este modelo se construir de
tal forma que represente el problema formulado, tanto en la parte geomtrica como fsica (Fig.
3.2).
Se evaluar el sistema en rgimen transiente y el flujo ser laminar, las condiciones de
velocidad de entrada del aire sern uniforme al igual que la temperatura. La corriente de aire
ser bajo flujo cruzado y la temperatura de la pared de los cilindros se mantendr constante.
Se utilizarn de igual forma los modelos de Weakly compressible Navier Stokes y General
heat transfer ambos en rgimen transiente para el programa. Para construir el modelo se elige
la fsica en primer lugar, para luego construir la geometra del problema, en este caso, se debe
escoger para la transferencia de calor las opciones Fluid termal interaction y Non isothermal
Flow (Figura 3.1)
En lo que concierne la parte geomtrica, se deber generar mallas para discretizar esta,
que permitan la resolucin del sistema. Las mallas deben ser las del tipo boundary layer que
se afinan en la zona alrededor de los tubos, donde hay altos gradientes de velocidad y
temperatura.
La regin de anlisis es rectangular. A la entrada en x=0, se impone una velocidad dada, y la
salida corresponde a una condicin de flujo convectiva. Los bordes horizontales sin roce y con
condicin de simetra. En el cilindro, se usa condicin de no deslizamiento (velocidad nula).
La metodologa para la construccin del modelo es:
Fsica
Geometra
Propiedades y condiciones de borde
Mallado
Definicin del solver
Resolucin
La metodologa de avance ser la siguiente:
Simulacin en modelo bsico con flujo cruzado (1 tubo)
Ajuste de parmetros
Creacin de un modelo de haz de tubos (6 a 10)
Simulacin y resultados

La configuracin del software Comsol Multiphysics v3.5 para poder realizar este tipo de
simulaciones requiere de diversas variables que sern explicadas a continuacin.
29

En primer lugar los mdulos utilizados para la resolucin del problema son:
Weakly compressible Navier Stokes

General Heat Transfer
Estos mdulos aparecen automticamente del Model Navigator, una vez que se escoge la
fsica del problema que se desea resolver. En este caso, se escogi un mdulo de transferencia
de calor, para luego escoger una interaccin trmica del fluido en un rgimen transiente en 2
dimensiones.
En la figura 3.2 se puede apreciar la interfaz del Model Navigator.

Figura 3.1: Interfaz del navegador para encontrar el modulo indicado

3.1 Geometra y condiciones de borde

La geometra del problema consiste en un volumen de control rectangular que
representa un ducto, con un tubo circular dentro de este espacio permitiendo la interaccin
30

entre el fluido y el cilindro, todo esto en dos dimensiones. Este volumen de control es solo para
hacer posible la simulacin, puesto que lo que se desea estudiar es un tubo nico en una
corriente de aire sin otras paredes. Los valores de la regin de anlisis y del dimetro del tubo
variarn segn los resultados de las simulaciones, para este caso se utiliz un dimetro D=0.05
[m], un largo H=0.1 [m], y un ancho W=0.7 [m].
Las propiedades fsicas utilizadas en este problema son constantes. La figura 3.2 muestra las
etiquetas de cada uno de los bordes de la figura, luego, las condiciones de borde impuestas para
la simulacin de un tubo con flujo cruzado de aire son las siguientes:


Figura 3.2: Modelo de tubo con flujo cruzado con etiquetas en los bordes.

Flujo compresible
Flujo de aire en la entrada igual al de la corriente de aire en el eje x, u= cte., v=0 (borde
1)

Simetra en los bordes horizontales (borde 2 y 3)

Presin de salida igual a 0 (borde 4)

Condicin de no deslizamiento alrededor del cilindro. (borde 5,6,7,8)

31

Transferencia de calor
Aislacin de los bordes (bordes 2 y 3)

Temperatura de la entrada constante T=293.15K (borde 1)

Temperatura de la pared de los tubos constante T=253.15K (borde 5,6,7,8)

Flujo convectivo en la salida (borde 4)

3.2 Parmetros de Ecuacin

El segundo paso de programacin consiste en el seteo(configuracin) de los
parmetros para el solver. En primer lugar, se escoge un paso de tiempo pequeo (0.1) puesto
que el fenmeno que se est estudiando se desarrolla en un lapso corto, luego el fenmeno se
comporta de forma peridica. Este paso de tiempo indica el lapso de tiempo en el cual el
software grabar ese instante en su memoria ese, puesto que el verdadero paso de tiempo que
realiza el programa, se autodefine sujeto a las ecuaciones que est resolviendo. Mientras ms le
cueste al programa resolver las ecuaciones para su convergencia, menor ser el paso de tiempo.
32


Figura 3.3: Seteo del Solver Parameters.
Otra variable importante para la simulacin consiste en la tolerancia absoluta. Como se observa
en la figura 3.3, la tolerancia est programada para aceptar valores con un margen de hasta
1x10
-2
para las velocidades, y de 1x10
-1
para la temperatura. Este tipo de variables tiene gran
trascendencia en las simulaciones del modelo puesto que permite estabilizar estas ltimas.
33


Figura 3.4: Mtodo general para el Solver Parameters
El mtodo escogido para resolver las ecuaciones fue el mtodo BDF (Backward Differentiation
Formulas). Este mtodo (a diferencia del generalized Alpha) es conocido por su estabilidad, sin
embargo puede provocar severos efectos de Damping o amortiguamiento, que en este caso
son despreciables. El mtodo BDF es una familia de mtodos implcitos para la integracin
numrica de ecuaciones diferenciales ordinarias. Son mtodos lineales multipaso que para una
funcin y tiempo dado, aproximan la derivada de esa funcin utilizando la informacin de
tiempos ya computados, aumentando la precisin de la aproximacin.
Adems, este mtodo permite mayor velocidad en las simulaciones, lo que optimiza el tiempo
de simulacin y estabiliza la convergencia del modelo.



34

3.3 Criterios de convergencia

Se estableci como criterio de convergencia cuando la simulacin alcanza el tiempo final
programado. En la figura 3.5, se puede observar un grfico del paso inverso de tiempo en
funcin del tiempo para un modelo que converge. El software elige un paso de tiempo muy
pequeo en un principio, luego vara en funcin de la convergencia. Cuando la simulacin
converge, el paso de tiempo se hace estable despus de unas primeras iteraciones, y luego
alcanza un valor constante hasta alcanzar el tiempo final.

Figura 3.5: Paso de tiempo inverso en una simulacin

Generalmente en los casos donde las simulaciones divergen, el paso de tiempo se va
haciendo cada vez ms pequeo, hasta que su avance es prcticamente nulo. Como se trata del
paso de tiempo inverso, a medida que crece la curva, menor es el paso de tiempo, y viceversa
(Figura 3.6)


35


Figura 3.6: Paso de tiempo inverso en una simulacin que diverge
De esta misma forma, se puede apreciar el error en funcin del nmero de iteraciones. Para el
caso de una representacin que converge, el error flucta entre 10
-2
a 10
1
, y presenta
pendientes negativas en cada iteracin (Figura 3.7). En este caso las curvas estn tan cerca que
no se nota el efecto de varias pendientes negativas.

Figura 3.7: Error en funcin del nmero de iteraciones

36

4. Anlisis de modelos
4.1 Modelo bsico
4.1.1 Sensibilidad de la simulacin al mallado (nmero de elementos
de la malla)

El modelo bsico consiste en un tubo con flujo cruzado de aire. La idea principal de trabajar con
un tubo consiste en validar los resultados obtenidos por el programa a travs de una
comparacin con los resultados que se encuentran en la teora.
Una de las principales variables que maneja el programa es la cantidad de elementos o nodos
que se encuentran en la malla, que son determinantes en la precisin del clculo, y por ende en
el tiempo de simulacin de los modelos. Es por eso que es importante llegar a un nmero
ptimo de elementos, de manera de optimizar el tiempo de simulacin y llegar a una solucin
precisa.
Condiciones del tubo para la modelacin:


(Densidad del aire)



( Viscosidad dinmica del aire)



(Conductividad trmica del aire)



(Capacidad calrica del aire)



(Temperatura de entrada del aire)



(Temperatura de la superficie del tubo)



(Velocidad de entrada del aire)



(Dimetro del tubo)

(largo del rectngulo)

(Ancho del rectngulo)








37


Para analizar la cantidad de elementos a utilizar se realizaron simulaciones con 1302,
8320, 35858 y 67158 elementos. Se extrajo de estas simulaciones la variacin del
coeficiente convectivo en funcin del arco sobre el permetro del tubo (en la literatura
aparece tambin en funcin del ngulo). Como existen condiciones de simetra, la figura
4.1 muestra la variacin solo para la mitad del permetro expuesta al flujo.


Figura 4.1: Caso 1302 elementos, variacin del coeficiente convectivo en funcin del permetro del tubo.
Figura 4.2: Caso 8320 elementos, variacion del coeficiente convectivo en funcin del permetro del tubo

38



Figura 4.4: Caso 67158 elementos, variacin del coeficiente convectivo en funcin del permetro del tubo

Como se puede constatar en las figuras 4.1 a 4.4, a medida que aumenta el nmero de
elementos, el efecto sierra se ve disminuido. Esto se nota, ya que para una cantidad de
elementos relativamente baja como lo son las figuras 4.1 y 4.2, la curva presenta inestabilidad y
una variacin grande. A medida que se aumenta la cantidad de elementos (superior a 20000
nodos en las figuras 4.3 y 4.4), la curva no presenta estas variaciones, por lo que se concluye
que a partir de esa cantidad de elementos se pueden obtener soluciones de mayor precisin.
Figura 4.3: Caso 35858 elementos, variacin del coeficiente convectivo en funcin del permetro del tubo

39

Es importante notar que una mayor cantidad de nodos esta relacionado directamente con el
tiempo de estas simulaciones, por lo que es de mucha importancia optimizar el tiempo de
simulacin sin afectar la precisin del resultado.

4.1.2 Optimizacin del tiempo de simulacin

La optimizacin del tiempo es fundamental en este tipo de simulaciones, debido a que va ligado
directamente a la capacidad del ordenador que es una limitante para realizar las iteraciones, y
de la misma forma se pueden obtener resultados de una misma precisin con el manejo de
algunas variables en un menor tiempo.
Con la funcin Boundary integration en la pestaa de Post processing se obtuvo la curva del
flujo de calor (en [W/m]) en funcin del tiempo. La figura 4.5 muestra los primeros 30s de
simulacin con un paso de 0.1 para 1302 elementos.

Figura 4.5: Flujo total de calor en funcin del tiempo de 0 a 30s

Se puede observar que a partir de los 15s se alcanza una funcin peridica, por lo tanto para
obtener el valor promedio del flujo de calor en los primeros 30 segundos, se despreci los
valores obtenidos en el rango entre 0 a 15 s. De aqu se obtuvo un valor promedio de Q= 9.86
[W/m].
De igual forma, se simul para los siguientes 30s, con el mismo paso, e imponiendo en la
funcin solver parameters que se iniciara la simulacin con el resultado ya obtenido en los
40

primeros 30 segundos de simulacin (Current solution). Luego, se obtuvo la curva del flujo de
calor total en funcin del tiempo para los siguientes 30 segundos como se observa en la figura
4.6, que corresponde al estado final.


Figura 4.6: Flujo total de calor en funcin del tiempo de 30 a 60s

A diferencia de los primeros 30s, la curva es peridica en su totalidad, por lo tanto no es
necesario limpiar ciertos valores. Se obtuvo un valor promedio del flujo de calor Q=9.9 [W/m],
que se acerca bastante al anteriormente encontrado.
Se realiz este proceso para los distintos casos, 1302, 35858, y 67158 elementos que se
resumen en la tabla 4.1.
Tabla 4.1: Flujo de calor en funcin del tiempo para 3 casos, 1302, 35858,67158 Elementos.
Q [W/m]
Elementos 0 a 30s 30 a 60s
1302 9.86 [W/m] 9.9 [W/m]
35858 9.87 [W/m] 9.88 [W/m]
67158 9.86 [W/m] 9.85 [W/m]

Se concluye que en un rango de 15 segundos (a excepcin de la primera iteracin) con un paso
de 0.1 y con una malla superior a 20000 elementos se puede obtener resultados precisos
optimizando de esta forma el tiempo de simulacin.
41


4.1.3 Validacin de los coeficientes convectivos

La manera de obtener el coeficiente convectivo como output desde el software de Comsol,
consiste en la integracin de bordes sobre el permetro del tubo (Boundary Integration).
Como se vi en el punto anterior, el coeficiente convectivo es una funcin peridica, que
flucta en un rango de valores, por lo que la integracin de bordes sobre el permetro del tubo
en un tiempo especfico puede ser una mala aproximacin.
Por lo que, la forma de obtener el coeficiente convectivo promedio alrededor de los tubos
consiste en sacar el valor medio en todo el rango de tiempo estudiado, como se hizo en el punto
anterior.

Figura 4.7: Campo de velocidad para modelo de cilindro nico con flujo cruzado de aire

A partir de la ecuacin de Hilpert, se obtuvo que para una velocidad de entrada de 0.1 m/s, el
valor del coeficiente convectivo es de h=4.64 [W/m2.K]. Por otro lado a partir del mtodo
Boundary integration, y calculando el valor medio de la curva se obtuvo que h=3.14 [W/m2.K].
Los valores no son iguales entre la teora y la simulacin, sin embargo son del mismo orden. Este
error se puede producir por diversas razones, tanto en la geometra como en el mallado, por lo
que se requerir de un posterior anlisis para ver si se pueden mejorar los resultados. A
continuacin, se puede observar en la tabla 4.2, una tabla comparativa de los valores del
coeficiente convectivo calculados por el software Comsol y los calculados tericamente.
42

Tabla 4.2: Tabla Comparativa para valores entregados por Comsol y valores tericos
Velocidad [m/s] Re Q[w/m] h
1
[w/m^2*K] h terico Error %
2

0.1 316.22 9.85 3.14 4.64 32.48
0.2 632.43 14.04 4.47 6.41 30.33
0.3 948.65 17.85 5.68 7.75 26.67
0.4 1264.86 21.09 6.71 8.86 24.23
0.5 1581.08 24.04 7.65 9.83 22.16

Efectivamente, a medida que el valor de la velocidad de entrada aumenta, el error entre el
coeficiente convectivo terico y el que se obtuvo mediante Comsol disminuye. Este error
disminuye a un 22% con una velocidad de 0.5 m/s. Por lo que es de esperar, que para una
velocidad superior el error se pueda minimizar a cero. Este resultado es coherente, debido que
a mayor velocidad, mayor es el valor de h, por lo que el error cae a escalas menores (mientras
ms alto los valores de h, la diferencia se va haciendo despreciable).

Figura 4.8: Grfico de los coeficientes convectivos obtenidos por Hilpert y Comsol en funcin de la velocidad

Como se puede observar en la figura 4.8, el margen de diferencia es relativamente constante,
sin embargo a medida que los valores van creciendo, este error disminuye como se observa en
la figura 4.9.


1
Se calcula a partir de la formula




2
El error se calcula a partir de la formula


0,00
2,00
4,00
6,00
8,00
10,00
12,00
0,00 0,10 0,20 0,30 0,40 0,50 0,60
C
o
e
f
i
c
i
e
n
t
e

c
o
n
v
e
c
t
i
v
o

h
[
W
/
m
2
*
k
]

Velocidad [m/s]
h Comsol
h teorico
43


Figura 4.9: Porcentaje de Error en funcin de la velocidad

No obstante, una limitacin que se puede dar es la viabilidad de aumentar las velocidades. Esto
debido a que en un primer lugar, la simulacin se hace cada vez ms inestable, lo que hace que
el modelo no converja, y por otra parte, el tiempo de simulacin es demasiado largo (para una
velocidad de 0.5 m/s, el tiempo aproximado de simulacin fue de 36 horas).
Por otro lado, si observamos la figura 2.11, donde se grafican las correlaciones para calcular el
coeficiente convectivo en funcin de la velocidad, podemos notar que existe una correlacin
donde los valores se acercan an ms al obtenido mediante el mtodo elementos finitos, que es
la correlacin de Zhukauskas.
Realizando el mismo anlisis anterior, obtenemos la tabla 4.3:
Tabla 4.3: Tabla Comparativa para valores entregados por Comsol y valores tericos (Correlacin Zhukauskas)
Velocidad [m/s] Re Q[w/m] h [w/m^2*K] h Zhukaukas Error %
0.1 316.22 9.85 3.14 4.17 24.73
0.2 632.43 14.04 4.47 5.89 24.14
0.3 948.65 17.85 5.68 7.22 21.25
0.4 1264.86 21.09 6.71 8.33 19.42
0.5 1581.08 24.04 7.65 9.31 17.85

Efectivamente disminuye el error, los resultados son ms precisos, sin embargo esta precisin
radica en la utilizacin de una correlacin distinta. Notamos el mismo fenmeno anterior, a
medida que se aumenta la velocidad, el error disminuye hasta alcanzar un valor del orden del
17% para una velocidad de 0.5 m/s.



0,00
5,00
10,00
15,00
20,00
25,00
30,00
35,00
0,00 0,10 0,20 0,30 0,40 0,50 0,60
E
r
r
o
r

%

Velocidad [m/s]
error
44

Esto se puede apreciar mejor en las figuras 4.10 y 4.11 donde se observa la comparacin del
error entre el modelo creado, y las correlaciones Hilpert y Zhukauskas.

Figura 4.10: Comparacin coeficientes convectivos obtenidos por Hilpert, Zhukauskas y Comsol en funcin de la velocidad


Figura 4.11: Porcentaje de Error en funcin de la velocidad segn correlaciones

0,00
2,00
4,00
6,00
8,00
10,00
12,00
0,00 0,10 0,20 0,30 0,40 0,50 0,60
C
o
e
f
i
c
i
e
n
t
e

c
o
n
v
e
c
t
i
v
o

h
[
W
/
m
2
*
k
]

Velocidad [m/s]
h Comsol
Hilpert
Zhukauskas
0,00
5,00
10,00
15,00
20,00
25,00
30,00
35,00
0,00 0,10 0,20 0,30 0,40 0,50 0,60
E
r
r
o
r

%

Velocidad [m/s]
Hilpert
Zhukaukas
45

4.2 Pruebas de mejora Modelo Bsico

4.2.1 Cambio en el mallado

Debido a la falta de precisin en los resultados obtenidos, se busc los cambios necesarios para
ajustar los valores de h. En primer lugar, se not que no era necesario mallar la parte interior
del cilindro, lo que provoco cambios en el resultado final.

Figura 4.12: Campos de velocidad modelo sin mallado al interior del cilindro
A igual que el modelo anterior, este convergi hasta una velocidad de v=0.5 [m/s], para
velocidades mayores el modelo divergi debido a que alcanzo una singularidad (May have
reached a singularity). Se insisti cambiando la tolerancia de los parmetros correspondientes
(velocidad, temperatura, presin) para hacer el modelo ms estable, sin embargo no se
obtuvieron mejores avances.
Los resultados obtenidos para este primer modelo en lo que refiere al coeficiente convectivo
alrededor del tubo no se acercaron al valor terico calculado, sin embargo se apreci un cambio
notorio en su valor. En primer lugar, los valores de h pasaron a ser de menores en el primer
modelo estudiado (con malla interior) a muy superiores en este modelo sin mallado al interior
del tubo, como se observa en la tabla 4.4.



46


Tabla 4.4: Comparacin del valor h entre el modelo con y sin malla al interior del tubo.
Velocidad [m/s] Re Q[w/m] h* [w/m^2*K] h malla int.
0,10 316,22 39,44 6,28 3,14
0,20 632,43 56,18 8,94 4,47
0,30 948,65 69,87 11,12 5,68
0,40 1264,68 81,59 12,98 6,71
0,50 1581,08 92,21 14,68 7,65

Como se puede observar a primera vista, el coeficiente convectivo es casi el doble del primer
modelo. No encontramos explicacin a este cambio puesto que los dos modelos tienen como
condicin de borde temperatura constante alrededor del cilindro, sin embargo este resultado
nos permite variar parmetros para obtener mayor precisin en los resultados.
En la tabla 4.5 se observa el porcentaje de error del coeficiente convectivo obtenido del
programa comparado a las distintas correlaciones que se conocen en la literatura.
Tabla 4.5: Error entre coeficiente convectivo obtenido por Comsol para cilindro sin malla al interior v/s correlaciones
conocidas
Vel.
[m/s]
h*
[w/m^2*K] h Hilpert
Error
% Zhukauskas
Error
% Churchill B
Error
% Whitaker
Error
%
0,10 6,28 4,64 35,16 4,17 -50,67 4,71 -33,22 4,53 -38,69
0,20 8,94 6,41 39,40 5,89 -51,79 6,64 -34,71 6,63 -34,91
0,30 11,12 7,75 43,51 7,22 -54,12 8,13 -36,70 8,30 -34,02
0,40 12,98 8,86 46,56 8,33 -55,86 9,41 -38,01 9,74 -33,36
0,50 14,68 9,83 49,29 9,31 -57,56 10,54 -39,22 11,03 -33,09

Como se puede observar, el porcentaje de error va en aumento a medida que se aumenta la
velocidad, excepto para la correlacin de Whitaker, donde la diferencia sigue siendo grande, sin
embargo es constante del orden del 33%. En la figura 4.13 se observa este fenmeno.

47



Figura 4.13: h Comsol v/s h correlaciones
Se puede observar que las correlaciones empricas tienen un valor similar, y la curva que ms se
acerca a los resultados encontrados en el modelo Comsol es la curva de Whitaker. El valor de la
pendiente de las curvas de Comsol y Whitaker es prcticamente el mismo, por lo que se puede
concluir que existe un factor proporcional entre estos dos.

0,00
2,00
4,00
6,00
8,00
10,00
12,00
14,00
16,00
0,00 0,10 0,20 0,30 0,40 0,50 0,60
C
o
e
f
i
c
i
e
n
t
e

c
o
n
v
e
c
t
i
v
o

[
W
/
m
2
*
K
]

Velocidad flujo de aire [m/s]
h Comsol
h Hilpert
Zhukauskas
Churchill B
Whitaker
48

4.2.2 Cambios geomtricos

Para ajustar los resultados obtenidos, lo primero que se realiz fue cambiar la geometra del
problema. En este caso, el dimetro del tubo se mantuvo constante, sin embargo se vari el
volumen de control.
El largo del volumen de control se aument en 1.5 veces su tamao.

Figura 4.14: Modelo extendido
A pesar del cambio, los resultados fueron similares a los encontrados en la parte anterior.
Debido a esto, y adems que el modelo tiene mayor cantidad de elementos al ser ms grande,
se decidi volver al modelo anterior.
Para hacer los resultados ms precisos, se decidi cambiar la altura del volumen de control. En
este caso se hizo un modelo con un volumen de control de 0.2x0.7m como se observa en la
figura 4.15.
49


Figura 4.15: Modelo ancho 0.2x0.7 mts
En este modelo se not un cambio significativo en el coeficiente convectivo alrededor del tubo.
Efectivamente el hecho de tener un volumen de control angosto afectaba los valores de h. A
partir de este modelo, se realizaron simulaciones hasta una velocidad v=0,5 m/s, y los
resultados se muestran en la tabla 4.6.
Tabla 4.6: Error entre coeficiente convectivo obtenido por Comsol modelo ancho 0.2 mts v/s correlaciones conocidas
Vel.
[m/s]
h*
[w/m^2*K]
h
Hilpert
Error
% Zhukauskas
Error
% Churchill B
Error
% Whitaker
Error
%
0,10 5,17 4,64 -11,33 4,17 -24,17 4,71 -11,07 4,53 -15,54
0,20 7,49 6,41 -16,77 5,89 -27,10 6,64 -14,43 6,63 -14,59
0,30 9,39 7,75 -21,19 7,22 -30,08 8,13 -17,34 8,30 -15,08
0,40 9,98 8,86 -12,64 8,33 -19,79 9,41 -6,86 9,74 -2,92
0,50 11,05 9,83 -12,4 9,31 -18,63 10,54 -5,46 11,03 -0,46

Como se puede observar, los resultados se acercan mucho ms al valor terico de las
correlaciones, sin embargo, para las correlaciones de Hilpert, Zhukauskas y Churchill Berstein, el
error crece a medida que se aumenta la velocidad hasta 0,3 [m/s]. Por otro lado, y como se dio
para el modelo anterior, la correlacin de Whitaker mantiene un error constante del orden del
15% hasta una velocidad de 0,3 [m/s]. Para v=0,5 [m/s] obtenemos mejoras en el resultado
final, esto se debe principalmente a que se cambi la tolerancia absoluta de 0,005 a 0,01 para
las velocidades u y v. Este cambio se impuso como solucin puesto que el modelo no
converga a esas velocidades. Debido a que se observ una mejora en los valores de h, se
prob cambiando nuevamente el alto del volumen de control hasta llegar a un valor donde h
50

se haca invariante. Este valor se alcanz para una altura de 0,3 [m] con un dimetro de de tubo
de 0,05 [m], es decir, 6 veces el dimetro.
Tabla 4.7: Error entre coeficiente convectivo obtenido por Comsol para modelo ancho 0.3 mts v/s correlaciones conocidas
modelo ancho final
Vel.
[m/s]
h*
[w/m^2*K] h Hilpert
Error
% Zhukauskas
Error
% Churchill B
Error
% Whitaker
Error
%
0,10 5,03 4,64 -8,30 4,17 -20,73 4,71 -7,63 4,53 -12,10
0,20 7,31 6,41 -13,93 5,89 -24,05 6,64 -11,38 6,63 -11,54
0,30 9,13 7,75 -17,88 7,22 -26,60 8,13 -13,86 8,30 -11,60
0,40 9,57 8,86 -7,97 8,33 -14,82 9,41 -1,89 9,74 2,05
0,50 10,60 9,83 -7,82 9,31 -13,79 10,54 -0,62 11,03 4,59

Como se observa en este caso, con un ancho de 0,3 [m] se disminuye el error, aproximndose
al comportamiento de un medio infinito. Una vez ms la correlacin de Whitaker mantiene un
error constante hasta 0.3 [m/s], del orden del 11%. Luego, al igual que para el modelo de ancho
0,2 [m], se tuvo que cambiar la tolerancia lo que mejoro notoriamente los resultados, llegando
en algunos casos a errores del orden del -0,62%.
El cambio en la tolerancia para las simulaciones con velocidad superior a 0,3 m/s revel cambios
en los resultados finales. Surgi por lo tanto la inquietud, que esta variable pueda afectar los
valores de h para velocidades menores. Es por eso que se simul el mismo modelo, esta vez
con una tolerancia de 0,01 desde v=0,1 [m/s]. Los resultados se aprecian en la siguiente tabla:
Tabla 4.8: Comparacin h entre modelo Comsol y correlaciones para una tolerancia de 0,01
Vel.
[m/s]
h*
[w/m^2*K] h Hilpert
Error
% Zhukauskas
Error
%
Churchill
B
Error
% Whitaker
Error
%
0,10 5,03 4,64 -8,30 4,17 -20,73 4,71 -7,63 4,53 -12,10
0,20 7,30 6,41 -13,83 5,89 -23,95 6,64 11,27 6,63 -11,43
0,30 8,18 7,75 -5,51 7,22 -13,30 8,13 -0,57 8,30 1,69
0,40 9,25 8,86 -4,35 8,33 -10,97 9,41 1,96 9,74 5,90
0,50 10,53 9,83 -7,12 9,31 -13,06 10,54 0,11 11,03 5,33

Finalmente se puede concluir que los valores de h no cambian una vez que se vara la
tolerancia, pues mantienen un error del mismo grado de magnitud. La fluctuacin que se
produce es provocada por el sistema de ecuaciones del mismo software. Las correlaciones de
Hilpert, Churchill Berstein y Whitaker mantienen un error del orden del 10%, por lo que se
validan los resultados del coeficiente convectivo.




51

4.2.3 Campos de temperatura y velocidad

A partir de la representacin de los campos de temperatura y velocidad que se producen
alrededor del tubo se puede determinar qu tipo de rgimen es el que se forma alrededor del
cilindro. En la figura 4.16 se puede observar la temperatura en un instante t del modelo.

Figura 4.16: Distribucin de la temperatura en el modelo de velocidad v=0.3 m/s, t=30s
Este grfico es una fiel ilustracin de lo que sucede en la realidad, y puede ser muy til para ver
en qu tipo de rgimen se est desarrollando. Alrededor del cilindro encontramos bajas
temperaturas como se impuso en las condiciones de borde, y debido al flujo que atraviesa el
cilindro observamos la creacin de vrtices de Von Karman. Efectivamente, si comparamos esta
figura, con la figura 2.4, se puede ver inmediatamente que se est trabajando con Reynolds del
orden de 1x10
2
. Justamente en este modelo, se est trabajando con un Reynolds de 948 que
equivale a una velocidad de 0.3 m/s.
Los remolinos ya se muestran inestables, y se desprenden alternativamente a ambos lados de
forma peridica. Estos remolinos presentan un intercambio trmico con el cilindro, por lo que
tienen una temperatura menor a la ambiente, pero superior al tubo. En este caso se alcanzan
temperaturas alrededor de 280 K en los remolinos, siendo la temperatura ambiente de 293 K y
la temperatura alrededor del tubo de 253 K, por lo que se forma el fenmeno de intercambio de
calor para acondicionar el ambiente. De la misma forma, se puede representar la distribucin de
gradiente de temperatura alrededor del cilindro.
52


Figura 4.17: Gradiente de temperatura en el modelo de velocidad v=0.3 m/s, t=30s
En la figura 4.17, el comportamiento del gradiente de temperatura expresa el mismo fenmeno
que se puede observar en la figura 4.4, que es la curva que representa el coeficiente convectivo
en funcin del dimetro del tubo. En el punto de estancamiento frontal, existe un alto gradiente
de temperatura, por lo que el intercambio de calor est en su mximo valor (lneas juntas en la
parte frontal con colores rojos). A medida que se avanza en el dimetro del tubo, este gradiente
disminuye por lo que el intercambio de calor es menor, y el menor gradiente se produce en la
zona posterior del cilindro, donde se mezclan el aire enfriado con el aire caliente del ambiente.








53

Al igual que la figura anterior, se puede representar los campos de velocidad para este modelo
como se observa en la figura 4.18.

Figura 4.18: Campo de velocidad en el modelo de velocidad v=0.3 m/s, t=30s
El flujo presenta las caractersticas de capa lmite como en los anteriores grficos. Se puede
destacar que en el punto de estancamiento frontal, la velocidad es nula, luego la velocidad
aumenta por efecto de capa lmite (gradiente de presin favorable) en el permetro del tubo
hasta el punto de desprendimiento, donde la velocidad del fluido vuelve a ser baja. Los vrtices
e inestabilidad del fluido que se van desprendiendo del tubo presentan velocidades altas
alternativamente a ambos lados del cilindro en forma peridica.
Si observamos este fenmeno en lneas de flujo, el fluido incide sobre el cuerpo y tiende a
seguir su contorno. Luego de atravesar el tubo, en la parte posterior se puede ver como las
lneas de flujo ms juntas siguen el movimiento de los vrtices que se crean. En efecto, las lneas
de flujo ms alejadas del cilindro, no muestran ningn tipo de perturbacin pero a medida que
las lneas se acercan al cilindro estas presentan pequeas oscilaciones.
54


Figura 4.19: Lneas de flujo de velocidad de tubo simple con velocidad v=0.3 m/s, t=30s
Para ver lo que sucede en el permetro del cilindro, se procedi a obtener una imagen en
acercamiento del tubo como lo ilustra la figura 4.20.
55


Figura 4.20: Acercamiento de figura 4.19 alrededor del tubo.
Este tipo de grficos, permite determinar ngulos de desprendimiento de capa lmite sobre la
superficie del cilindro, adems de representar los fenmenos que suceden aguas abajo de este
mismo.







56

A partir de estos resultados, se obtuvo la figura 4.21, que representa el flujo de calor alrededor
del cilindro en funcin del tiempo.
Figura
4.21: Flujo de calor alrededor del cilindro en funcin del tiempo para v=0.3 m/s

A diferencia de la figura 4.6, la curva no muestra un comportamiento peridico, pero si
oscilatorio. Para este tipo de curvas, no es posible determinar frecuencia nica, sin embargo se
pueden llegar a notar algunos fenmenos particulares que analizaremos en la figura 4.22, en la
cual se observa un acercamiento de este grfico para los valores de q entre los tiempos [28s;
30s].
Entre los [20s, 30s] observamos 31 peaks superiores, por lo que el perodo T de la curva es
T=10/31= 0.3225 [s], con una frecuencia promedio f=1/T= 3.1 [Hz].
Este grfico es producto del aumento de velocidad del fluido de entrada, por lo tanto los flujos
pasan a ser ms turbulentos, y por eso las curvas pierden su carcter de peridicos.
57


Figura 4.22: Acercamiento figura 4.21, para valores de q en el rango [28s; 30s]

La figura 4.22 corresponde a un acercamiento de la curva del flujo de calor representada en la
figura 4.21. Se puede observar como mencionamos anteriormente que la curva no presenta
periodicidad, sin embargo, notamos eventos oscilatorios en perodos T=0.3s o t=0.4s. En este
rango de tiempo, observamos 6 peaks superiores, por lo que el perodo T=2/6= 0.33 [s], con
una frecuencia f= 1/T= 3 [Hz].








58

En la figura 4.23 se representa las lneas de flujo, para los perodos [28s; 28.5s].

(a) (b)

(c) (d)

(e) (f)
Figura 4.23: Secuencia correspondiente a las lneas de flujo para los instantes 28 a 28.5 segundos; (a) 28s; (b) 28.1s; (c) 28.2s;
(d) 28.3s; (e) 28.4s; (f)28.5s
El flujo posee una misma forma en cada instante, no obstante, se puede observar que la
separacin del flujo ocurre en distintos lugares de la superficie. En la figura 4.23 observamos un
ciclo y medio del fenmeno.
59

4.2.4 Presin entorno al cilindro

La presin es otro de los elementos que se puede analizar en este tipo de fenmenos. El flujo
incide en el tubo, provocando una presin alta en el punto de estancamiento frontal. Debido a
este gradiente de presin, el fluido se acelera a medida que recorre el permetro del tubo, y la
presin va disminuyendo hasta el punto de desprendimiento, donde la fuerza del fluido no es
capaz de superar la viscosidad del fluido, por lo que se llegan a formar gradientes de presin en
direccin opuesta. La figura 4.21 muestra este fenmeno en funcin del permetro del tubo.

Figura 4.24: Comportamiento de la presin alrededor del cilindro para Re=948


60

5. Modelo Haz de tubos

5.1 Condiciones iniciales

A partir del modelo de tubo simple con flujo cruzado se establecieron ciertos parmetros
para la simulacin, adems de ciertos lmites con respecto a la cantidad de nodos a
implementar. Al igual que para el modelo anterior, la simulacin requiere un paso progresivo
para aumentar la velocidad de entrada, por lo que es condicin necesaria modelar para una
velocidad e ir aumentando esa velocidad en 0.1 [m/s].
Las condiciones impuestas para el arreglo de tubos fueron las siguientes:


(Densidad del aire)



( Viscosidad dinmica del aire)



(Conductividad trmica del aire)



(Capacidad calrica del aire)



(Temperatura de entrada del aire)



(Temperatura de la superficie del tubo)



(Velocidad de entrada del aire)



(largo del rectngulo)

(Ancho del rectngulo)
















61

5.2 Modelo St/D=Sl/D=1.5

El primer modelo simulado consiste en 8 tubos formando un triangulo equiltero (Figura 5.1),
donde las medidas utilizadas estn definidas segn los valores tpicos de intercambiadores de
calor vistos en el captulo 2.3.2, adems de considerar el factor del cociente entre las distancias
tanto horizontal como vertical divididas por su dimetro como se observa en la tabla 2.3. El
dimetro escogido para este arreglo de tubos es de 1 pulgada. A continuacin se puede apreciar
su geometra:
(Dimetro del tubo)

(Distancia vertical entre centro de tubos)



(Distancia horizontal entre centro de tubos)



(Distancia plano diagonal entre tubos)



Figura 5.1: Esquema arreglo de tubos 53409 elementos

62

En la figura 5.2 se puede observar el campo de temperatura del fluido para el primer modelo de
arreglo de tubos estudiado. El flujo va de izquierda a derecha con velocidad uniforme de 0,5
[m/s] hasta que choca con la primera fila de tubos. Se puede notar que una vez que el fluido
choca con el tubo, se produce un intercambio de calor, donde los flujos interiores de la estela
poseen temperaturas inferiores a la del ambiente, y el contorno de color rojo representa el
lmite entre la temperatura del ambiente, y la que provoca el choque con los tubos. Posterior al
paso por los tubos, se pueden constatar la formacin vrtices que se producen principalmente
por la geometra impuesta en este modelo.


Figura 5.2: Campo de temperatura del fluido para modelo St/D=Sl/D=1.5






63

La figura 5.3 representa el campo de gradiente de temperatura del modelo. El aire entra con
una temperatura de 20[C] y el contorno del tubo se encuentra a -20[C]. Una vez que el aire
topa con los tubos, se aprecian cambios de temperatura que son producto de la turbulencia. El
aire caliente se queda estancado aguas abajo del tubo, provocando los mayores gradientes.
Luego se produce el frente de aire a la salida de la primera fila de tubos donde existe mayor
transferencia. Este gradiente de temperatura disminuye en la tercera fila, debido a que el aire
viene ms fro, por lo que el gradiente trmico es menor. En la zona central del arreglo, se
producen turbulencias, y el recorrido del flujo se hace impredecible.

Figura 5.3: Gradiente de temperatura St/D=Sl/D=1.5







64

En la figura 5.4 se aprecia el campo de velocidades del fluido. En el primer punto de contacto
del fluido y el tubo se puede constatar el punto de estancamiento en donde la velocidad tiende
a cero. Luego se produce el efecto capa lmite, a medida que el fluido recorre el permetro del
tubo se produce un aumento en la velocidad inducido por el aumento de presin que ejerce el
punto de estancamiento. En teora, la velocidad mxima se alcanza en los puntos donde la
presin es mnima, y esta se produce aproximadamente a unos 90 del cilindro, es decir, en un
cuarto de permetro donde se produce el desprendimiento de capa lmite. Luego, por el origen
de las vorticidades, las velocidades tienden a cambiar su direccin, por lo que se alcanzan
nuevos mnimos locales en la parte posterior de los tubos.


Figura 5.4: Campo de velocidad del fluido St/D=Sl/D=1.5
El flujo que incide sobre la segunda fila del arreglo, viene influenciado por la presencia de la
primera fila del cilindro, esto se observa directamente puesto que en la segunda fila, el punto de
estancamiento es ms difuso, y en la tercera fila, este punto ya casi no se observa. Esto se debe
a que el flujo ya no incide directamente sobre el tubo, solo en el cilindro central se puede
apreciar un pequeo punto de estancamiento. Se produce el mismo fenmeno ya explicado, a
medida que la presin es menor, la zona de aceleracin es menor por lo que el fluido sale a
bajas velocidades.
65

Finalmente otro concepto importante a notar, es que el conjunto de tubos, o banco de tubos
pareciera comportarse como un nico tubo gigante, donde se produce el efecto de capa lmite,
y a su alrededor se produce una aceleracin del fluido.


Figura 5.5: Campo de presin St/D=Sl/D=1.5

En la figura 5.5 se puede apreciar el fenmeno anteriormente explicado, en el punto donde
incide directamente el flujo se producen altas presiones (punto de estancamiento), y posterior a
los cilindros, se tienen zonas de ms baja presin.
Con respecto al anlisis de transferencia de calor, se pudieron obtener los siguientes resultados
que se grafican en la tabla 5.1, donde se compar directamente con la frmula de Zhukaukas
para analizar el comportamiento del coeficiente convectivo global. El flujo de calor q
representa la suma del flujo de calor alrededor de cada uno de los tubos del arreglo.




66

Tabla 5.1: Comparacin h global del arreglo St/d=Sl/D=1.5 v/s h Zhukauskas
V [m/s] Vmax[m/s] Sup. Remax Nusselt
h Zhukauskas
[w/m^2*K]
q comsol
[W/m]
h comsol
[w/m^2*K]
Error
%
0,2 0,6 A1 964 14,48 14,82 222,50 8,71 41,22
0,3 0,9 A1 1450 18,47 18,91 292,01 11,44 39,52
0,4 1,2 A1 1930 21,95 22,46 355,80 13,93 37,97
0,5 1,5 A1 2410 25,09 25,69 400,47 15,68 38,94


El coeficiente convectivo (h comsol) se obtuvo a partir del campo promedio entre los h de los 8
tubos del haz. Se calcul a partir de la siguiente ecuacin:





Estos resultados muestran un error del orden del 39% con respecto a la correlacin de
Zhukauskas para Remax>1000.
A partir de las frmulas vistas en el captulo 2, se obtuvo que la velocidad de referencia (Vmax)
pasa a travs de la seccin vertical A1 del arreglo, y luego a partir de esta velocidad se obtuvo
Remax, Nusselt y el coeficiente convectivo respectivamente. Vmax es por definicin una
velocidad de referencia para poder calcular la transferencia de calor en el banco de tubos, pero
no corresponde a la velocidad mxima dentro del arreglo.
De la tabla 5.1, se puede apreciar que el valor de Vmax es mucho mayor a la velocidad de
entrada. Para u=0.5 [m/s], el valor de Vmax es de 1.5 [m/s], siendo que en la figura 5.4 del
campo de velocidad, la velocidad mxima es de 0.929 [m/s], osea la de referencia debiera ser
an menor. Esa diferencia en los campos de velocidades parece ser la principal razn por la cual
existe ese margen de error entre los h.
Para entender porque se produce ese margen de error, se calcul la velocidad de referencia a
travs del Software. Como se observa en la figura 5.6, se obtuvieron los valores de la velocidad
en esa seccin del arreglo.

67


Figura 5.6: Campo de Velocidad en seccin A1
A travs de la funcin Cross-section plots parameters se obtuvo la figura 5.7. Estos valores se
exportaron para obtener la velocidad de referencia del programa.

Figura 5.7: Velocidad en funcin del arco de seccin A1



68

La velocidad de referencia se calcul mediante la ecuacin:


Donde
Sumatoria de las velocidades en casa punto de la seccin A1
: Paso elegido por el programa
C: Distancia vertical entre permetro de los tubos.
Finalmente, con estos resultados, se calcul el coeficiente convectivo mediante el mtodo de
Zhukauskas. En la tabla 5.2 se pueden observar los resultados.
Tabla 5.2: Comparacin h global del arreglo St/d=Sl/D=1.5 con velocidad de referencia de Comsol v/s h Zhukauskas
V [m/s]
Vmax
comsol[m/s] Sup. Re Nusselt
h Zhukauskas
[w/m^2*K]
q comsol
[W/m]
h comsol
[w/m^2*K]
Error
%
0,2 0,206 A1 331 8,21 8,40 222,50 8,71 -3,68
0,3 0,327 A1 525 9,87 10,10 292,01 11,44 -13,21
0,4 0,49 A1 787 12,08 12,37 355,80 13,93 -12,65
0,5 0,604 A1 971 13,36 13,68 400,47 15,68 -14,64

En primer lugar, hay que destacar que la velocidad de referencia obtenida mediante los datos
extrados del campo de velocidad es muy cercana a la velocidad de entrada. Debido a que el
Reynolds obtenido es inferior a 1000 para todas las velocidadas estudiadas, se utiliz una
segunda correlacin de Zhukauskas construida para casos con estos valores de Reynolds.
Efectivamente, los resultados se acercan mucho ms a los de Comsol, obteniendo mrgenes de
error del orden de 3 a 15%, lo que es comn dentro de este tipo de aplicaciones.









69

5.3 Modelo St/D=2; Sl/D=1.5

El segundo modelo consiste en un arreglo de 8 tubos, con la siguiente disposicin:
(Dimetro del tubo)

(Distancia vertical entre centro de tubos)



(Distancia horizontal entre centro de tubos)



(Distancia plano diagonal entre tubos)



Figura 5.8: Esquema arreglo de tubos 51509 elementos modelo St/D=2; Sl/D=1.5
70


Figura 5.9: Campo de temperatura modelo St/D=2; Sl/D=1.5
Este segundo modelo posee una mayor distancia vertical entre los tubos. Se produce el mismo
fenmeno que para el modelo anterior, pero se puede notar que al haber ms espacio entre los
tubos, el flujo que sale de cada tubo est ms definido, y la interaccin del flujo entre los tubos
parece no mezclarse en las filas 2 y 3 a diferencia de lo que ocurre en la figura 5.2. De la misma
forma, se puede observar que los vrtices que se crean abarcan un espacio mayor, pero las
temperaturas de estos vrtices parecen ser mayores al de la figura ya mencionada (del orden de
280 K), por lo tanto menor el enfriamiento.
71


Figura 5.10: Campo de velocidad modelo St/D=2; Sl/D=1.5
En la figura 5.10 se puede apreciar el campo de velocidad para el segundo modelo estudiado.
Aqu se pueden apreciar los puntos de estancamiento para prcticamente todos los tubos, a
excepcin del tubo central de la tercera fila, donde el fluido no incide directamente sobre el
tubo, sin tener un contacto frontal. Como existe un mayor distanciamiento vertical entre los
tubos, el flujo, al pasar alrededor del tubo de la segunda fila choca directamente con el tubo de
la tercera fila, a diferencia de la figura 5.4 donde el flujo no choca en la parte frontal, si no que
en la cara ms diagonal del tubo.
Tabla 5.3: Comparacin h global del arreglo St/d=2 Sl/d=1.5 v/s h Zhukauskas
V [m/s] Vmax[m/s] Sup. Re Nusselt
h Zhukauskas
[w/m^2*K]
q comsol
[W/m]
h comsol
[w/m^2*K]
Error
%
0,2 0,23 A1 369 7,73 7,91 215,60 8,44 -6,76
0,3 0,35 A1 562 9,53 9,75 285,40 11,18 -14,62
0,4 0,49 A1 782 11,51 11,78 349,50 13,69 -16,19
0,5 0,57 A1 916 12,66 12,95 378,99 14,84 -14,58




72

5.4 Modelo St/D=1.5; Sl/D=2

El tercer modelo consiste en un arreglo de 8 tubos, con la siguiente disposicin:
(Dimetro del tubo)

(Distancia vertical entre centro de tubos)



(Distancia horizontal entre centro de tubos)



(Distancia plano diagonal entre tubos)




Figura 5.11: Esquema arreglo de tubos 51509 elementos modelo St/D=1.5; Sl/D=2
73


Figura 5.12: Campo de temperatura modelo St/D=1.5; Sl/D=2

Este modelo posee una mayor distancia horizontal con respecto a sus tubos. Se observa un
fenmeno similar al de la figura 5.2, donde los vrtices que se forman entre la fila 2 y 3 se
mezclan, producindose grandes gradientes de temperatura.







74

Figura
5.13: Campo de velocidad modelo St/D=1.5; Sl/D=2
La figura 5.13 refleja el campo de velocidad que se produce para este modelo. Se observan 6
claros puntos de estancamiento en la parte frontal de los tubos. En primer lugar, los tres tubos
de la primera fila, luego los tubos de la segunda fila, y finalmente solo el tubo central de la
tercera fila. El flujo que sale de los tubos de la segunda fila parece no chocar con los tubos
externos de la tercera fila, efectivamente solo el flujo que va por el centro, tiene contacto
directo con el tubo central, donde se produce el efecto de capa lmite. El flujo no alcanza a
rodear completamente el tubo, y las lneas de flujo se vuelve a formar a una distancia larga
despus de chocar con los tubos, esto se debe principalmente a que se trabaja con velocidades
pequeas.
Tabla 5.4: Comparacin h global del arreglo St/d=1.5 Sl/d=2 v/s h Zhukauskas
V [m/s] Vmax[m/s] Sup. Re Nusselt
h Zhukauskas
[w/m^2*K]
q comsol
[W/m]
h comsol
[w/m^2*K] Error %
0,2 0,21 A1 337 7,91 8,10 220,10 8,62 -6,47
0,3 0,33 A1 530 9,72 9,95 289,90 11,35 -14,10
0,4 0,48 A1 771 11,62 11,89 350,50 13,73 -15,41
0,5 0,58 A1 938 12,76 13,07 388,99 15,23 -16,59

75

5.5 Modelo St/D=Sl/D=2

El cuarto modelo consiste en un arreglo de 8 tubos, con la siguiente disposicin:
(Dimetro del tubo)

(Distancia vertical entre tubos)



(Distancia horizontal entre tubos)



(Distancia plano diagonal entre tubos)





Figura 5.14: Esquema arreglo de tubos 50508 elementos St/D=Sl/D=2
76


Figura 5.15: Campo de temperatura modelo St/D=Sl/D=2

En la figura 5.15 se puede observar el campo de temperatura del modelo St/D=Sl/D=2. Este
modelo tiene un mayor distanciamiento tanto vertical como horizontal de los tubos. Al tener
una mayor rea, los vrtices que se crean abarcan un mayor espacio dentro del volumen de
control definido. La primera fila del arreglo tiene un comportamiento similar al de los otros
modelos, pero a partir de la segunda fila se produce el efecto de capa lmite para todos los
tubos como si fueran un slido en conjunto.
77


Figura 5.16: Campo de velocidad modelo St/D=Sl/D=2
En la figura 5.16 se puede apreciar el campo de velocidad para este modelo. Se observan 5
puntos de estancamiento en la parte frontal de los tubos. El flujo del centro flucta, pasando
tanto por abajo (como en la figura) como por arriba del tubo central de la tercera fila. Esta
fluctuacin se debe a que al espacio entre tubos es mayor, por lo tanto se observa el fenmeno
de vrtices entre tubos.

Tabla 5.5: Comparacin h global del arreglo St/d=Sl/d=2 v/s h Zhukauskas
V [m/s] Vmax[m/s] Sup. Re Nusselt
h Zhukauskas
[w/m^2*K]
q comsol
[W/m]
h comsol
[w/m^2*K]
Error
%
0,2 0,22 A1 353 7,62 7,80 212,30 8,31 -6,59
0,3 0,334 A1 537 9,43 9,65 283,50 11,10 -15,02
0,4 0,443 A1 710 10,94 11,20 332,78 13,03 -16,36
0,5 0,548 A1 880 12,16 12,45 367,55 14,39 -15,61

78

5.6 Anlisis Interaccin de tubos

El arreglo geomtrico se muestra de forma esquemtica en la figura 5.15. Las condiciones del
flujo estn dominadas, como vimos en el punto anterior, por los efectos de la separacin de
capa lmite y la interaccin entre las estelas, que por supuesto influyen directamente sobre la
transferencia de calor por conveccin.

Figura 5.17: Arreglo esquemtico para representacin de cada tubo.

A partir de los resultados obtenidos anteriormente, se analiz el coeficiente convectivo de cada
tubo del arreglo, para los cuatro casos distintos, a una velocidad v=0.5 [m/s]. Los resultados se
presentan a continuacin en la tabla 5.6.





79

Tabla 5.6: Valor del coeficiente convectivo para cada tubo del arreglo, (a) St/d=Sl/d=1.5; (b) St/d=2 Sl/d=1.5;(c) St/d=1.5
Sl/d=2; (d) St/d=Sl/d=2.

(a)

Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
1.1 57,02 17,72 2.1 59,76 18,58 3.1 40,5 12,59
1.2 56,8 17,66 2.2 58,9 18,31 3.2 38,5 11,97
1.3 57 17,72 3.3 39,9 12,40

(b)
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
1.1 52,84 16,55 2.1 56,1 17,44 3.1 41,2 12,81
1.2 52,6 16,48 2.2 55,2 17,16 3.2 31,5 9,79
1.3 52,12 16,33 3.3 40,3 12,53

(c)
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
1.1 54,68 17,13 2.1 55,5 17,39 3.1 36,5 11,44
1.2 55,42 17,36 2.2 54,93 17,21 3.2 39,6 12,41
1.3 54,93 17,21

3.3 36,1 11,31
(d)
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
Tubos
q
[W/m]
h
[W/m2.k]
1.1 50,29 15,76 2.1 52,54 16,46 3.1 39,87 12,49
1.2 49,77 15,59 2.2 52,02 16,3 3.2 33,42 10,47
1.3 50,47 15,81 3.3 39,25 12,3

En primer lugar se puede notar que el coeficiente de transferencia de calor est asociado
directamente por su posicin en el banco de tubos. El coeficiente para la primera lnea es similar
para los cuatro casos, pero lo ms notorio es que tiene el mismo valor que al de un solo tubo
con flujo cruzado, esto se debe a que se trata de la primera lnea del arreglo. La diferencia entre
los valores radica principalmente por la creacin de estelas, a medida que los tubos estn ms
cerca hay un mayor gradiente de temperatura por lo que la transferencia por conveccin es
mayor.
Se puede notar adems que los tubos que poseen un mayor coeficiente de transferencia de
calor estn asociados con los tubos entre lneas, y para este caso se trata de la segunda fila. Los
tubos de la primera lnea actan como una rejilla de turbulencia, aumentando el coeficiente
convectivo para las lneas siguientes. Sin embargo, en la tercera lnea apreciamos menores
valores de h, esto principalmente porque se trata de la ltima fila, por lo que no se forman
gradientes muy grandes posterior al choque del fluido con esos tubos.
En lo que concierne la comparacin entre los arreglos estudiados, se puede apreciar que para el
arreglo donde los tubos estn ms cercanos (a) se observa que los coeficientes de transferencia
80

de calor son mayores. Esto define un tipo de comportamiento, mientras ms cercano son los
tubos, mayor es la transferencia de calor.
Por otro lado, los arreglos (b) y (c), que corresponden a los arreglos con mayor espaciado
vertical y mayor espaciado horizontal respectivamente, se observa un comportamiento
relativamente similar, pero que las diferencias se explican bien a travs de las figuras vistas
anteriormente. Para la primera fila, en el caso (b) que corresponde al caso con mayor espaciado
vertical, tenemos un menor valor de h con respecto al caso (c), esto deriva directamente de la
diferencia de espaciado vertical entre los dos casos, a menor distancia, mayor transferencia.
En la segunda fila, los dos muestran un comportamiento muy similar, pero lo ms importante es
que los valores del coeficiente de transferencia de calor por conveccin son mayores en la
segunda fila para ambos casos, lo que corrobora lo mencionado anteriormente.
Para la tercera fila, la diferencia se aprecia en la las figuras 5.8 y 5.11. En la primera (caso (b),
espaciado vertical mayor), el flujo que sale de la segunda fila por efecto capa lmite incide
directamente sobre los tubos de la tercera fila, por lo tanto se aprecia un valor de h mayor. En
cambio en la figura 5.11 (caso (c) espaciado horizontal mayor), el flujo que sale de la segunda
fila solo golpea el tubo central, y por efecto de rebote el flujo incide de forma suave sobre los
tubos extremos.
Finalmente el caso (d), que corresponde al arreglo de tubos con mayor separacin tanto vertical
como horizontal, los valores de h son menores para las 3 filas. De aqu se puede concluir en
primer lugar que los valores del coeficiente convectivo esta sujetos a la posicin dentro del
arreglo, y en segundo lugar a la distancia que existe entre tubos.
Se debe mencionar de igual forma que estos valores estn sujetos a los mismos errores
encontrados en el estudio anterior, por lo que marcan una tendencia correcta, sin embargo sus
valores pueden diferir de valores encontrados mediante correlaciones empricas.









81

6. Conclusiones

En la presente memoria se logr crear un modelo para el flujo y la transferencia de calor
alrededor de bancos de tubos a travs del software de elementos finitos Comsol Multiphysics
v3.5.
En primer lugar, se analiz un modelo de un tubo simple con flujo de aire cruzado. Se impuso la
fsica del problema, seteandolo en el programa. Luego se cre la geometra del modelo, que
constituye el volumen de control y los tubos, para finalmente establecer las condiciones de
borde del problema. Entre las condiciones impuestas, se estableci trabajar con temperaturas
constantes en las paredes de los tubos y a la entrada del flujo, adems de velocidades
constantes. Esta simulacin se logr realizar para velocidades de 0.1 a 0.5 [m/s], siendo esta
ultima la ms alta permitida por el modelo, debido a que a mayores velocidades las
simulaciones divergan. Los resultados fueron mejorando a medida que se modificaban variables
dentro del estudio. Se pudo notar que existen variables determinantes dentro de este tipo de
problemas, tales como la cantidad de elementos que se usa, la tolerancia del solver, y la
geometra del modelo. Finalmente se logr la validacin del coeficiente convectivo para un tubo
nico con flujo cruzado, comparndolo con las principales correlaciones conocidas en la
literatura. Se llego hasta un 10% de error entre ambos valores, lo que es una muy buena
aproximacin.
El estudio del tubo simple con flujo cruzado, nos entrego una variacin oscilatoria del flujo de
calor con el tiempo, y un comportamiento muy similar al que se describe en la literatura con los
vrtices de Von Karman entre otros fenmenos.
Finalmente, este primer anlisis ayud para las futuras simulaciones. Se establecieron
parmetros fijos tales como las propiedades del fluido, y se permiti optimizar el tiempo en las
simulaciones sin perder precisin en los resultados gracias a la optimizacin de la cantidad de
nodos implementados.
En una segunda etapa, se cre un modelo de simulacin 2D transiente del flujo de aire y
transferencia de calor para un intercambiador de tubos de flujo de aire cruzado. Este modelo
consisti en cuatro casos, cada arreglo con 8 tubos, y un dimetro externo d=0.0254 [m] (1
pulgada). Los casos se diferenciaban en su distancia vertical u horizontal con respecto al centro
de los otros tubos lindantes.
Se pudo observar, en un primer instante, que las condiciones de cantidad de nodos impuesta
para la primera parte del estudio no eran suficientes, esto debido a que el problema tena
mayores gradientes de temperatura entre lneas, por lo que era importante un mallado fino en
esas zonas. De igual forma, se cre un volumen de control lo suficientemente grande para no
afectar los valores de h, y se dispuso en forma simtrica los tubos para que no haya variacin
entre los tubos en las extremidades. Una vez conseguido todos estos pasos, se procedi a la
resolucin del problema. Se consiguieron campos de velocidad, de temperatura y lneas de flujo
para los cuatro modelos simulados.
82

A diferencia del caso de tubo nico con flujo cruzado, los resultados obtenidos del software
para el haz de tubos y las correlaciones empricas de la literatura para arreglo de tubos con flujo
cruzado mostraron, en un principio, una diferencia mayor en el valor del coeficiente de
transferencia de calor por conveccin. Esto por el uso de una correlacin emprica (Zhukauskas)
basada en una ecuacin para el clculo de la velocidad de referencia. A travs de un clculo para
obtener la velocidad de referencia dentro del arreglo en el programa, se pudo mejorar
notablemente los resultados en la parte terica, y de esa forma convalidar los resultados
obtenidos por Comsol.
En cuanto al fenmeno y comportamiento del fluido en un haz de tubos, se puede decir que la
simulacin se asimila bastante a la teora de la literatura. En efecto, se produce el efecto de
capa lmite alrededor de los tubos, con un punto de estancamiento en la parte frontal del
cilindro, donde choca el fluido, debido a la alta presin que produce este mismo. Los gradientes
de temperatura crecen en las lneas internas del arreglo de tubos, y las estelas producidas por el
choque del fluido con el tubo inciden en los valores de h.
Finalmente se realiz el estudio de cada tubo dentro del arreglo, para cada uno de los casos con
velocidad de entrada del aire v=0.5 [m/s]. La primera lnea del arreglo tuvo el mismo
comportamiento que un tubo solo con flujo cruzado de aire, lo que se acerca a la realidad. De
igual modo, los valores de h de los tubos internos mostraron un aumento en cuanto al
coeficiente de transferencia de calor pero la tercera fila entrego valores bajos. En la teora, este
comportamiento es correcto puesto que la primera fila provoca la turbulencia, por lo que la
lneas siguientes presentan gradientes mayores de temperatura, que por ende conllevan a un
aumento de h, sin embargo no se produce eso para la tercera fila.
Los errores en la prediccin de los valores del coeficiente convectivo pueden deberse a diversas
razones. En primer lugar, se trabaja con propiedades constantes, siendo que en la realidad, las
propiedades fsicas del fluido cambian con la temperatura y presin. De igual forma, los errores
de simulacin provienen de las dificultades en la prediccin de altos gradientes de temperatura
alrededor del tubo. En efecto, se observ una mejor prediccin de h al alejar los bordes de la
regin. Esto implica que una mayor cercana del cilindro a otras superficies impide aproximar
bien el gradiente de temperatura sobre ste. Este hecho explica los mayores errores observados
en los coeficientes convectivos en el banco de tubos, en que la proximidad de los tubos dificulta
la correcta evaluacin del campo de temperatura entre stos. Otro punto observado, fue la
velocidad, y por ende el Reynolds con la cual se calculan los valores de Nusselt. Para las
correlaciones se utilizan las ecuaciones vistas en el captulo 2, donde a travs de una ecuacin
encontramos la velocidad mxima del fluido dentro del arreglo. Notamos que esta velocidad era
mayor a la observada en la simulacin, lo que finalmente arrastra una diferencia en los clculos
para obtener el valor de h.
Todas las comparaciones de los resultados de esta memoria con otros estudios se han hecho a
partir de correlaciones empricas que se encuentran en la literatura. Las correlaciones de
fenmenos de conveccin forzada tienen tpicamente errores del orden del 20%, derivado de
las condiciones del equipo experimental usado, de errores de medicin o simplemente de
arrastre o propagacin de errores en los clculos [3]. En efecto, los resultados de la
transferencia de calor para unas mismas condiciones se diferencian de un 20,3 hasta 48%, por lo
83

que esto puede repercutir en la sub o sobreestimacin del coeficiente de conveccin, y por lo
tanto, de todo el clculo trmico del equipo de intercambio de calor. [6]
A partir de este criterio, los resultados del estudio para tubo nico pueden considerarse vlidos.
Para el caso de banco de tubos hay mayores dificultades por lo que se requiere mayor precisin
para validar ciertos modelos, sin embargo se pudieron obtener diferencias del 5 al 20%.
Para poder juzgar la precisin de los resultados obtenidos con el modelo hubiese sido til
comparar con resultados de otras simulaciones numricas. Esto porque no hay suficientes
predicciones de coeficientes convectivos para bancos de tubos en rgimen laminar. Adems, la
ecuacin de Grimison no es vlida para Re<2000, el cual es prcticamente el lmite superior de
Reynolds en este estudio.
Parte de un estudio a futuro, consistira en simular estos modelos para velocidades mayores,
acercndose de esta forma a problemas que surgen en aplicaciones de la industria. Otro estudio
interesante seria realizar estas simulaciones en 3D, agregando variables tales como el largo del
tubo.

















84

Bibliografa

[1] Francisco Bozzo, Determinacin de coeficientes convectivos en el exterior de evaporadores de
refrigeracin mediante simulacin., 2010, vol. Memoria para optar al ttulo de Ingeniero Civil
Mecanico.
[2] R. Torres y J. Grau, Introduccin a la Mecnica de Fluidos y Transferencia de Calor con Comsol
Multiphysics. Barcelona, 2007.
[3] Frank P. Incropera, Fundamentos de Transferencia de Calor, Cuarta Edicin ed.
[4] sadik kaka, Boilers, evaporators and condensators., 1991.
[5] Prof. Jesus F. Ontiveros, "Operaciones unitarias II, Dimensionamiento de equipos de transferencia
de calor".
[6] Georgiy Polupan Yevgen Pysmennyys, "Estudio comparativo de los mtodos del calculo de
transferencia de calor en bancos de tubos," 2009.
[7] R.J. Goldstein S. Sanitjai, "Forced convection heat transfer from a circular cylinder in crossflow to air
and liquids," may 2004.
[8] R.P. Chhabra, V. Eswaran V.K. Mandhani, "Forced convection heat transfer in tube banks in cross
flow," september 2001.
[9] Harif Guerrero, Jesus Rojo Miguel Angel Diaz, "Simulacion y analisis de flujo transversal a traves de
bancos de tubos en lnea usando el programa CFX," 2008.
[10] Anthony F. Mills, Transferencia de calor.
[11] Tamotsu Igarashi Hajime Nakamura, "Variation of Nusselt number with flow regimes behind a
circular cylinder for Reynolds numbers from 70 to 30000".
[12] J.R. Culham, M.M Yovanovich W.A. Khan, "Convection heat transfer from tube banks in crossflow:
Analytical approach," march 2006.
[13] Discussion Forums Comsol.
http://www.comsol.com/community/forums/general/thread/13859/.
[14] Robert Parsons, ASHRAE Handbook., 2001.

También podría gustarte