Tesis Transporte Calor Comsol
Tesis Transporte Calor Comsol
Tesis Transporte Calor Comsol
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
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
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]
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
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
Donde
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:
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:
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)