USBTransTub: PROGRAMA COMPUTACIONAL QUE CALCULA LA
TRANSFERENCIA DE CALOR EN TUBERAS
Alejandra Sous Guevara
Miniproyecto de Ingeniera Qumica
dirigido por
Dosinda Gonzlez-Mendizbal
Departamento de Termodinmica y Fenmenos de Transferencia
RESUMEN
El programa USBTransTub, permite obtener la temperatura de salida de un fluido que circula por una tubera y
que intercambia calor con el exterior. El usuario puede especificar el tipo de fluido externo e interno, el material,
dimensiones y propiedades de la tubera, as como el tipo de conveccin y orientacin del fluido, teniendo la
posibilidad de incluir radiacin y aislantes. El programa cuenta con una base de datos que incluye las
propiedades para los fluidos y materiales de la tubera y aislantes. La aplicacin permite obtener la temperatura
de salida por dos vas: el caso integral y el caso diferencial, en este ltimo el usuario elige el nmero de
diferencias en las que quiere subdividir la tubera. USBTransTub, es un programa interactivo, rpido y de fcil
comprensin, que arroja resultados precisos. Fue desarrollado para Windows
98 (u otro ms actual), utilizando
como lenguaje de programacin Visual Basic
versin 6.0. Gracias al estilo de interfaz grfica Windows
, es
muy fcil de manipular por los usuarios permitindole adems, una mayor atencin al anlisis de resultados al no
emplear mucho tiempo en tareas repetitivas asociadas a este tipo de clculo.
INTRODUCCIN
A travs de la tecnologa de participacin utilizada en
Internet y el surgimiento de universidades y cursos virtuales, es
vital la bsqueda de propuestas con mayor interaccin entre la
teora de la transferencia de calor y su comprensin acadmica.
A travs de clculos como lo son, en este caso, la T
sal
en
tuberas expuestas a la prdidas o ganancia de calor por
procesos de enfriamiento o calentamiento de fluidos que
circulan por tuberas y la determinacin del aislante ptimo
para minimizar esta transferencia de calor. En la prctica, estos
clculos son muy tediosos y repetitivos ya que requieren datos
especficos de propiedades que dependen de la temperatura.
Por esta razn se propuso como objetivos para este trabajo:
Establecer un modo ms didctico para la comprensin de
la transferencia de calor en tuberas.
Futura inclusin en pginas web (cursos virtuales)
desarrollndose una mayor interaccin entre el curso y el
estudiante.
Agilizar el clculo de la temperatura de salida y espesor de
aislante ptimo.
Realizar tablas comparativas, grficas al exportar los
resultados a Excel
.
Para cumplir estos objetivos se propuso el programa
computacional que tuviese las siguientes caractersticas:
Programa computacional ejecutable, realizado en el
lenguaje de Microsoft Visual Basic
(versin 6.0)
Base de datos en Access
para la obtencin de las
propiedades como la conductividad y la viscosidad entre
otras. En la Fig. 1 se encuentra una de las grficas con su
respectiva lnea de tendencia utilizada (Polinmicas
cbicas)
AGUA
y = 1,996E-04x
3
- 2,191E-01x
2
+ 7,973E+01x - 5,380E+03
R
2
= 9,833E-01
3000
3500
4000
4500
5000
5500
6000
6500
7000
250 300 350 400 450 500 550 600 650
Temperatura [K]
C
p
[
J
/
k
g
K
]
Figura 1. Tendencia polinmica para las propiedades (Ej. Cp
del Agua).
FUNDAMENTOS TERICOS
Haciendo referencia a la tubera mostrada en la Fig. 3 un
balance de energa en estado estacionario:
Cambio de entalpa Calor intercambiado
del fluido interno = con los alrededores (1)
Suponiendo que el fluido interno no cambia de fase, se
puede escribir:
( )
D
= = -
Rt
T
Q T T Cp m
r sal e
&
(2)
donde Q
r
= Calor intercambiado con los alrededores.
Partiendo de la ec. (1), en la siguiente figura se observan los
flujos de calor que entran y salen de la tubera.
Q
conv int
Q
cond A
Q
cond C
Q
cond B
Tb
Q
conv ext
Q
entra
=m Cp (T
e
-T
sal
)
Q
rad
Figura 2. Balance de energa en un corte longitudinal de la
tubera.
Existe una analoga entre la difusin de calor y la carga
elctrica. De la misma manera que se asocia una resistencia
elctrica con la conduccin de electricidad, se asocia una
resistencia trmica con la conduccin de calor (Incropera y
De Witt, 1999).
Al considerar entonces la Fig. 3, se mostrar la expresin
que se utiliza fundamentalmente en este programa matemtico,
donde se desprecian las resistencias trmicas de contacto.
Siguiendo su nomenclatura, se tiene la siguiente expresin para
la transferencia de calor:
he k k k hi
sal
e
sal e
r
R R R R R
)
T T
T T
ln(
T T
Q
C B A
+ + + +
-
-
-
=
(3)
Cada una de las resistencias se irn describiendo a
continuacin. Es importante destacar que se est considerando
el diferencial medio logartmico de temperaturas como ms
preciso que el diferencial de temperaturas algortmico.
T
e
T
sal
T
b
T
sup1
T
sup4
T
sup3 T
sup2
T
Figura. 3 Distribucin de temperaturas para una pared
compuesta.
De la ec. (3) tenemos esta primera resistencia:
L r 2 h
1
R
1 1
hi
p
=
(4)
donde h
1
se calcula con las correlaciones de:
a) Flujo laminar (Re<2300) (Perry et al, 1992):
Hausen (Gz<100):
14 . 0
s
b
3 / 2
Gz 047 . 0 1
Gz 085 . 0
66 . 3 Nu
+
+ =
m
m
(4.a)
Sieder-Tate (Gz>100):
14 . 0
s
b 3 / 1
Gz 86 . 1 Nu
m
m
=
(4.b)
b) Flujo de transicin (2300<Re<10.000):
Hausen (Perry et al, 1992)
( ) A Pr 125 Re 116 . 0 Nu
3 / 1 3 / 2
D
- =
14 . 0
s
b
3 / 2
L
D
1 A donde
+ =
m
m
(4.c)
Sieder-Tate (Escartn y Ramrez, 1996)
( )
( )
[ ]( )
( )
[ ]
b
L
a
T
Pot
H
Pot
H H
j j j
2 1
= (4.d)
c) Flujo turbulento (Re>10.000):
Sieder-Tate (Perry et al, 1992):
14 . 0
s
3 / 1
D
D
D L
Pr Re
86 . 1 Nu
m
m
=
(4.e)
Colburn (Incropera y De Witt, 1999):
(4.f)
3 / 1 5 / 4
D D
Pr Re 023 . 0 Nu =
Dittus-Boelter (Incropera y De Witt, 1999):
(4.g)
n 5 / 4
D D
Pr Re 023 . 0 Nu =
donde n = 0.4 para calentamiento (T
s
> T
b
) y n = 0.3 para
enfriamiento (T
s
< T
b
) cumpliendo las siguientes condiciones:
10
D
L
000 , 10 Re
160 Pr 7 . 0
D
Todas estas correlaciones evalan las propiedades a la
temperatura de bulto (T
b
), que es un promedio aritmtico entre
la temperatura de entrada (T
e
) y temperatura de salida (T
sal
) del
fluido interno, excepto m
s
que est a la temperatura superficial
interna.
Volviendo a la ec. (2), luego de la resistencia convectiva
interna, descrita anteriormente, se tienen las resistencias
conductivas de:
R
, stas siguen como
expresin general, la siguiente:
C B A
k k k
R y R ,
Lk 2
) r r ln(
R
1 2
cond , t
p
= (5)
que como se mencion anteriormente, se deducen de la Ley de
Fourier en coordenadas cilndricas, es decir:
0
dr
dT
kr
dr
d
r
1
=
(6)
para las siguientes suposiciones:
Transferencia de calor unidimensional (direccin radial).
No hay generacin de energa.
Flujo de calor constante (Edo. Estacionario).
Por lo que, al usar como condiciones de borde para esta
ecuacin diferencial las respectivas T
superficiales
para cada caso,
manteniendo la nomenclatura de la ec. 2, se tiene la siguiente
expresin para la transferencia de calor:
) r r ln(
) T T ( Lk 2
Q
1 2
2 , s 1 , s
r
- p
=
(7)
y finalmente, tenemos la ltima resistencia incluida en la
ec. 2 que es:
L r 2 h
1
R
4 4
he
p
=
(8)
donde h
4
se calcula con una de las correlaciones que en
adelante se proponen, el usuario elegir segn el tipo de
conveccin y la orientacin del fluido y la correlacin a utilizar:
1. Para conveccin forzada y flujo perpendicular (Incropera y
De Witt, 1999):
Hilpert, propiedades a la temperatura de pelcula (T
f
):
3
(8.a)
/ 1 m
D D
Pr Re C
k
D h
Nu =
los valores de las constantes C y m se listan en la Tabla 1.
Tabla 1. Constantes para la ec. 8.a
Re
D
C M
0.4 4 0.989 0.805
4 40 0.911 0.618
40 4,000 0.683 0.466
4,000 40,000 0.193 0.385
40,000 400,000 0.027 0.805
Zhukauskas (1 < Re < 10
6
y 0.7 < Pr < 500), las
propiedades son, en este caso, a la temperatura de infinito
(T
), excepto Pr
s
que se evala a la temperatura superficial
externa (T
supe
):
4 / 1
s
n m
D D
Pr
Pr
Pr Re C Nu
=
(8.b)
Si Pr < 10, n = 0.37 y si Pr > 10, n = 0.36. Los valores de
las constantes C, m se listan en la Tabla 2.
Tabla 2. Constantes para la ec. 8.b
Re
D
C m
1 40 0.75 0.4
40 1,000 0.51 0.5
10
3
2 * 10
5
0.26 0.6
2 * 10
5
10
6
0.076 0.7
Churchill-Bernstein ( Re Pr > 0,2 ). Las propiedades
se evalan las propiedades de la temperatura de pelcula
(T
f
):
[ ]
5 / 4
8 / 5
D
4 / 1
3 / 2
3 / 1 2 / 1
D
D
000 . 282
Re
1 A donde
A
Pr) / 4 . 0 ( 1
Pr Re 62 . 0
3 . 0 Nu
+ =
+
+ =
(8.c)
2. Para conveccin forzada y flujo paralelo (Incropera y De
Witt, 1999):
Para Re<500.000, flujo laminar se tiene que:
3 / 1
L L
Pr Re 664 . 0 Nu = (9)
Para Re>500.000, flujo turbulento se tiene que:
(10)
3 / 1 8 . 0
L L
Pr ) 871 Re 037 . 0 ( Nu - =
2.1. En el caso de resolver la transferencia de calor en una
configuracin diferencial, se tiene una variacin en las
correlaciones que toma en cuenta el cambio en el
desarrollo del flujo (Re) con respecto a la longitud:
Para Re<500.000, flujo laminar la correlacin para el
coeficiente de transferencia de calor externo queda como:
) a b ( Pr
v
a b
664 . 0
h
2 / 1 2 / 1 3 / 1
2 / 1
-
m
r
-
=
-
(11)
Para 10
8
>Re>500.000, flujo turbulento quedar
como:
) a b ( Pr
v
a b
0164 . 0
h
5 / 4 5 / 4 3 / 1
5 / 4
-
m
r
-
=
(12)
donde a y b son los valores limites para el diferencial (Ver
Nomenclatura).
Para estas cuatro correlaciones se evalan las
propiedades a la temperatura de pelcula (T
f
).
3. Para conveccin libre y tubera en sentido vertical
(Incropera y De Witt, 1999):
Churchill y Chu (Se puede aplicar sobre todo
intervalo de Ra
L
) y se evalan las propiedades a la
temperatura de pelcula (T
f
):
[ ]
2
27 / 8
16 / 9
6 / 1
L
L
Pr) / 492 . 0 ( 1
Ra 387 . 0
825 . 0
+
+ = Nu
13.a)
Churchill y Chu corregida (Ra
L
<10
9
), se evalan las
propiedades a la temperatura de pelcula (T
f
) :
[ ]
9 / 4
16 / 9
4 / 1
L
L
Pr) / 492 . 0 ( 1
Ra 670 . 0
68 . 0 Nu
+
+ =
(13.b)
4. Para conveccin libre y tubera en sentido horizontal
(Incropera y De Witt, 1999):
Morgan (10
-10
< Ra
L
< 10
12
). Las propiedades se
evalan a la temperatura de pelcula (T
f
):
n
D D
CRa
k
D h
u N = =
(14.a)
donde los valores de las constantes C y n se listan en la Tabla 3.
Tabla 3. Constantes para la ec. 14.a
Ra
D
C m
10
-10
10
-2
0.675 0.058
10
-2
10
2
1.02 0.148
10
2
10
4
0.85 0.188
10
4
10
7
0.48 0.250
10
4
10
7
0.125 0.333
Churchill y Chu (Ra
L
<10
12
), se evalan las
propiedades a la temperatura de pelcula (T
f
):
[ ]
2
27 / 8
16 / 9
6 / 1
D
D
Pr) / 559 . 0 ( 1
Ra 387 . 0
60 . 0 u N
+
+ =
(14.b)
5. Para conveccin libre y tubera en sentido inclinado
(Incropera y De Witt, 1999):
Para este caso se trabajar slo con ngulos (a) menores de
40 desde la vertical, se usan las mismas ecs (13.a) y (13.b), con
la diferencia de utilizar la siguiente expresin para el clculo
del nmero de Rayleigh (Ra
L
):
f
3
sup
L
T
1
donde
L T T Cos g
Ra
= b
a n
- b a
=
(15)
En esta resistencia convectiva externa, se puede adicionar el
valor del coeficiente de transferencia de calor por radiacin,
llamado h
r
, en la Fig. 4 se presenta el balance en el caso que el
usuario quiera considerar la radiacin emitida por la tubera a
unos alrededores:
Figura 4. Intercambio de radiacin entre una superficie (tubera)
y los alrededores.
La velocidad total de transferencia de calor desde la
superficie es:
Q (16) ) T T ( ) T T ( hA Q Q
4
alr
4
s s rad conv
- es + - = + =
El coeficiente de transferencia de calor por radiacin hr que
se modela de forma similar a la conveccin:
(17) ) T T )( T T ( h
2
alr
2
s alr s r
+ - es
DESCRIPCIN DEL PROGRAMA
El programa USBTransTub est diseado con el siguiente
esquema de presentacin:
La organizacin del programa es estilo fichero, donde se
tienen dos carpetas principales llamadas: Configuracin del
sistema y Resultados, en la carpeta de Configuracin del
sistema se encuentran dos carpetas secundarias llamadas:
Lado externo e interno (Fig. 6.a), que incluye toda la
configuracin interna y externa en cuanto a fluidos se refiere y
a las restricciones anteriormente mencionadas. Cuando se elige
el tipo de conveccin y orientacin del fluido, las imgenes van
variando de acuerdo al tipo que se elige. En esta carpeta,
tambin est la opcin de elegir las correlaciones utilizadas para
el clculo de los coeficientes de la transferencia de calor, tanto
interno como externo, con todas las ecuaciones que se
presentaron en el captulo anterior.
Tubera (Fig. 6.b), donde se elige el tipo de material, la
longitud y dems componentes descritos a continuacin (Ver
Fig. 5):
Recubrimiento de aislante externo (se usa en la industria
como proteccin)
Recubrimiento de aislante interno. Ambos aislantes poseen
el mismo combo en lista de materiales, ver Tabla 4.
Recubrimiento de una coraza de acero inoxidable, aluminio
o cobre, con la opcin de poder personalizar otro material
incluyendo los coeficientes lineales de la conductividad.
Opcin de considerar radiacin en la coraza propuesta o
incluir una emisividad para una coraza personalizada.
Diversos materiales para la tubera.
Coraza
Aislante interno
Dimetro interno
Dimetro externo
Aislante externo
Radiacin
Emisividad
Figura 5. Corte longitudinal de la tubera
(a) Carpeta Lado externo e interno
(b) Carpeta Tubera
Figura 6. (a) y (b) Carpetas de la interfaz
Tabla 4. Materiales aislantes incluidos en la base de datos
(Incropera y De Witt, 1999)
Aislantes propuestos
Asbesto liso
Asbesto corrugado
Silicato de calcio
Vidrio celular
Lminas de vidrio
Magnesia (85%)
Granos de lana mineral
Lminas lana mineral
Poliestireno expandido
Poliuretano expandido
Corcho vegetal
Hoja de aluminio
Fibra asbesto
Fibra asbesto c/aglomerante
Fibra suelta de vidrio
Para la eleccin del material de la tubera, se defini una
sub categora en el tipo del material elegido, ste determina el
valor de la conductividad de la misma (Ver Tabla 5). Adems
se incluyeron especificaciones de dimetro nominal y catlogo
(Ver Tabla 6), que especifican el valor del espesor de la misma.
En este caso tambin se puede tener la oportunidad de
personalizar el material de la tubera incluyendo igualmente, los
coeficientes lineales de la conductividad.
Tabla 5. Sub tipo de materiales incluidos en la base de datos.
(Incropera y De Witt, 1999)
Material Sub Material
Acero Inoxidable
ANSI 302, ANSI 304, ANSI 316,
ANSI 347
Acero Comercial Carbn ordinario, Carbn-silicio,
Carbn-Mn-silicio, AISI 1010
Cobre Cobre puro, Bronce comercial,
Cartucho de latn, Constantan
Tabla 6. Tuberas incluidas de la base de datos (Normas ANSI, Ludwig, 1977) (Rodrguez-Alarcn y Guzmn, 1998)
Material Dimetro nominal Catlogo
Acero Inoxidable
1/8, , 3/8, , , 1, 1 , 1 , 2, 2,
3, 3 , 4, 5, 6, 8, 10, 12
5s, 10s, 40s, 80s
Acero Comercial , , 1, 1 , 1 , 2, 2 , 3, 3 , 4,
5, 6, 8, 10, 12, 14, 16, 18, 20, 24, 30
Tubera de pared estndar, Tubera de
pared extra fuerte, tubera de pared
doble extra fuerte, 20, 30, 40, 60, 80,
100, 120, 140, 160
Cobre , 3/8, , 5/8 Tipo K
a) Datos introducidos (fluidos)
Se propone una lista de fluidos interno y externos que estn
incluidos en la base de datos (Ver Tabla 7), adems de la
posibilidad de personalizar las propiedades de un fluido a
voluntad del usuario. (Ver Fig 8).
Valor de la temperatura de entrada de la tubera, como el
valor del caudal msico o velocidad (opcional).
Valor de la temperatura de infinito y velocidad externa, que
es automticamente cero cuando se estudian los casos de
conveccin libre.
Tabla 7. Fluidos incluidos en la base de datos (Incropera y De
Witt, 1999).
Fluidos propuestos
Agua [273,15<T<645]
Aire [100<T<3000]
Glicerina [273<T<320]
Etilenglicol [273<T<373]
Fren-12 [230<T<320]
Vapor de agua [380<T<850]
Figura 7. Interfaz en la personalizacin de propiedades.
b) Datos introducidos (restricciones)
Con la finalidad de que el programa sea ms didctico, se
consideraron algunos casos de clculo especiales denominado
Con restricciones (Ver Fig. 7.a). Por ejemplo, se pueden
incluir restricciones que reducen el tiempo de los clculos:
Coeficiente de transferencia de calor externo (h
e
).
Temperatura de superficie externa (T
supe
).
Coeficiente de transferencia de calor interno (h
i
).
Finalmente el clculo, como se mencion anteriormente, se
puede realizar integralmente o diferencialmente, en el primer
caso se efecta un solo clculo para toda la longitud de la
tubera, hallndose el valor de la temperatura de salida general.
A diferencia que para el caso diferencial, se van determinando
las temperaturas de salida para cada uno de los segmentos
fraccionados de la longitud total de la tubera, valor que el
mismo usuario puede elegir tenindose un mximo de 100
divisiones.
c) Estilo en la programacin
La programacin sigui un esquema modular, que es el
ptimo en la organizacin y distribucin de la data y desarrollo
general del algoritmo computacional.
Figura 8. Mdulos de programacin.
d) Diagrama de flujo
Una breve explicacin del algoritmo de clculo utilizado en
la programacin (Ver Fig. 9) comienza por igualar la
temperatura de salida con la temperatura de entrada. Este valor
slo sirve para poder comenzar la bsqueda en la base de datos
y posteriormente calcular las propiedades, de este modo se
determina el valor de los coeficientes de transferencia de calor,
de las resistencias y finalmente se encuentra un nuevo valor de
la temperatura de salida a travs de la ec. (18), que se compara
con la temperatura de salida inicial supuesta y con la cual se
repite de nuevo todo el proceso anterior.
+
-
=
T
)
Rt Cp m
1
exp(
T T
T
e
sal
&
(18)
Luego de cumplir con una tolerancia de 0.01, el programa
presenta resultados de todas las iteraciones realizadas, con todas
las propiedades y todos los valores, en general, empleados para
el clculo de la temperatura de salida.
e) Presentacin de los resultados
Estos se ubican en la carpeta llamada Resultados. Los
resultados se presentan de dos formas: los ms inmediatos se
presentan en el visor (Fig. 11) y luego se puede observar la
tabla con todos los valores de todas las iteraciones realizadas y
en sta se tiene la opcin de llevar a Excel para un
procesamiento posterior (Fig. 12).
Figura 9. Esquema del procedimiento.
Figura 10. Interfaz con visor y resultados.
Figura 11. Interfaz de la tabla de resultados.
RESULTADOS OBTENIDOS
Para observar la versatilidad del programa se trabaj con el
siguiente problema:
Una tubera de acero comercial de 100 m. de longitud,
dimetro nominal 2 in Catlogo 40, est cubierta por un aislante
de conductividad trmica k=0.1 W/m*K y espesor 2 cm. La
tubera se encuentra en un ambiente donde el aire sopla
transversalmente sobre la misma a una velocidad de 10 m/s, con
una temperatura promedio de 25 C. Si la tubera transporta
agua, la cual entra a 93 C, con una velocidad promedio de 1
m/s, calcule la temperatura de salida de la misma .
Resolucin:
1) Se introducen los siguientes datos y especificaciones
para los lados interno y externo (Ver Fig 6.a):
Fluido Externo: Aire
Tinf = 298 K
Velocidad = 10 m/s
Conveccin forzada y flujo perpendicular, usndose la
correlacin de Hilpert (predeterminada) para el he
Fluido interno: Agua
Te = 366 K
Velocidad = 1 m/s
2) Se introducen las especificaciones y valores
necesarios en la configuracin de la tubera (Ver Fig.
6.b):
Material: Acero comercial
Sub Tipo: Carbn ordinario
Dimetro nominal: 2
Catlogo: 40
Longitud: 100
Se personaliza el aislante: m (pendiente) = 0 y b
(coeficiente independiente) = 0.1, de modo que la
conductividad es independiente de la temperatura y
permanece constante.
3) Al oprimir el botn de Calcular (Ver Fig. 10)
aparecen los siguientes resultados en el visor:
Tsal = 365,1662 K
Fluido externo usado: Aire
Fluido interno usado: Agua
he = 47.0681 [W/m2 K]
hi = 6552.0652 [W/m2 K]
Q = 7322.1616 [W]
Sumatoria resistencias = 0.00923 [W/K]
A modo de comprobar las ventajas didcticas del programa,
se realiz el procedimiento descrito para las diferentes
configuraciones (ver Tabla 8). Es interesante observar el
incremento del calor en los casos de conveccin libre con
radiacin.
SINRADIACIN
Tipo de conveccin
yorientacindel fluidoenla tubera Tsal [K] He[W/mK] Q[W] # iteraciones
Conveccinforzada: flujo perpendicular 365,1662 47,0681 7322,1616 1
Conveccinforzada: flujo paralelo 365,1981 35,4059 7042,1618 2
Conveccinlibre: tubera vertical 365,6187 4,3449 3348,0823 3
Conveccin libre: tuberahorizontal 365,9992 0,0054 6,8942 1
Conveccin libre: tuberainclinada
1
365,6455 3,8568 3112,5473 3
CONRADIACIN
2
Tipo de conveccin
yorientacindel fluidoenla tubera Tsal [K] He[W/mK] Q[W] # iteraciones
Conveccinforzada: flujo perpendicular 365,0563 46,8723 8287,6854 2
Conveccinforzada: flujo paralelo 365,0535 35,4059 8312,3449 4
Conveccinlibre: tubera vertical 365,0543 1,0525 8304,7625 4
Conveccin libre: tuberahorizontal 365,0545 0,0018 8303,4254 4
Conveccin libre: tuberainclinada
1
365,0544 0,9235 8304,6037 4
1
=Angulo=40
2
=Coraza de Acero inox., espesor 1/2mmyTalr=150K
Tabla 8. Comparaciones para el caso ejemplo.
Igualmente este problema se plante tanto para el clculo
integral de la tubera como para un clculo diferencial de la
misma, obtenindose con este ltimo caso, grficas muy
interesantes sobre la disminucin de la temperatura con
respecto a la longitud. Obsrvese en la Fig. 12 la variacin de
la temperatura en el interior de la tubera donde se infiere que
en conveccin forzada existe mayor prdida de calor porque el
valor del coeficiente convectivo externo es mucho mayor que el
caso de conveccin libre, donde las prdidas son
aproximadamente despreciables:
365
365,2
365,4
365,6
365,8
366
366,2
0 50 100
Longitud [m]
T
[
K
]
Forzada: Flujo
Perpendicular
Forzada: Flujo
paralelo
Libre: Tubera
horizontal
Figura 12. Tendencia lineal de la temperatura de salida con
respecto a la distancia.
CONCLUSIONES
Los clculos realizados con USBTransTub han resultado
ser confiables ya que se corresponden completamente con
los clculos manuales.
La convergencia empleada en el mtodo descrito para la
obtencin de la temperatura de salida, ha sido efectiva y
reduce sustancialmente el nmero de iteraciones.
Se logran resultados inmediatos para clculos muy tediosos:
El uso de este programa reduce de manera significativa el
tiempo requerido para los clculos manuales y evita la
bsqueda en tablas y grficos al utilizar accesorios y
tuberas incluidos en la base de datos.
El programa permita la posibilidad de graficar y establecer
comparaciones entre las diversas configuraciones de la
tubera, lo que brinda al estudiante mayor versatilidad en la
materia.
La interfaz es muy fcil de manejar, agradable a asequible
al usuario.
FUTURAS MEJORAS
La programacin de dos casos: tubera enterrada y el caso
donde la tubera reciba un calor adicional.
Enriquecimiento de la base de datos.
La posibilidad de guardar archivos (.txt .dat).
Disear una ayuda mediante el sistema web (links).
Inclusin del caso de conveccin en cambio de fase.
NOMENCLATURA
a Longitud inicial del diferencial de tubera [m]
b Longitud final del diferencial de tubera [m]
D Dimetro de la tubera [m]
g Aceleracin de gravedad [m/s
2
]
h Coeficiente promedio de la transferencia de calor
[W/m
2
K]
hi Coeficiente convectivo del fluido interno [W/m
2
K]
heq Coeficiente convectivo del fluido externo [W/m
2
K]
k Conductividad [W/m K]
L Longitud [m]
m& Caudal msico [kg/s]
Nu Nmero de Nusselt
Pr Nmero de Prandtl [-]
Q Calor [W]
Ra Nmero de Rayleigh [-]
Re Nmero de Reynolds [-] definido como rVD/m
v Velocidad del fluido [m/s]
r Densidad [kg/m
3
]
m Viscosidad dinmica [Pa.s].
REFERENCIAS
Incropera, F., De Witt, D., 1999, Fundamentos de
Transferencia de Calor, Editorial Prentice-Hall, Captulos 1, 2,
3, 7, 8 y 9.
Perry, R., Green, D., Maloney, J., 1992, Manual del Ingeniero
Qumico, Editorial Mc Graw-Hill, Tomo 10. pp 10-14 y 10-16
Rodriguez-Alarcn, A., Guzmn, L., 1998, Desarrollo de un
programa para obtener perfiles de presin de fluidos
incompresibles en tuberas, Miniproyecto, Universidad Simn
Bolvar.
Escartn, T., Ramrez, G., 1999, Desarrollo de subrutinas en
Visual Basic para diseo y evaluacin de intercambiadores de
calor de doble tubo, tubo y coraza y placas empacadas,
Miniproyecto, Universidad Simn Bolvar.
