UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL

Curso
FENÓMENOS DE TRANSPORTE

CAPITULO 10
ECUACIÓN GENERAL DE BALANCE
EN ESTADO ESTACIONARIO
FUNDAMENTOS DEL TRANSPORTE
TURBULENTO

DV 
Prof. [Link]átegui R
NR   2100
wreategui@[Link] 
2021 1
 2 Log10 ( N R f )  0.8
f
[Link]átegui 1
El objetivo de este tema es dar a conocer los principios del
transporte turbulento en masa, calor y cantidad de
movimiento, y que los estudiantes hagan uso de las diferentes
correlaciones adimensionales para resolver ejercicios
propuestos por el profesor.

[Link]átegui R 2
Introducción
• La mayoría de las aplicaciones prácticas de los fenómenos de transporte
son más complejos que los sistemas hasta hoy descritos.
• Por ejemplo, en un intercambiador de calor un fluido en movimiento debe
calentarse. En ésta aplicación se transporta calor y cantidad de movimiento.
Si se mantiene una velocidad moderada, el calor y la cantidad de
movimiento se transfieren por mecanismos paralelos del tipo molecular y
turbulento.

[Link]átegui R 3
[Link]átegui R 4
[Link]átegui R 5
[Link]átegui R 6
 • El transporte turbulento se debe al movimiento desordenado de
grandes moléculas o conjuntos de ellas. Estos grupos de
moléculas se denominan remolinos.
• La masa, el calor y la cantidad de movimiento pueden transferirse
por transporte turbulento en cualquier fluido.
• El transporte turbulento no puede describirse completamente por
una expresión matemática única, debido a que el mecanismo es
mucho más complejo que el transporte molecular.
• En estos casos es necesario combinar la teoría con la evidencia
experimental para obtener un cuadro adecuado de lo que es la
turbulencia.
• Una técnica convencional usada en la correlación de datos
experimentales es la del análisis dimensional, donde grupos
adimensionales apropiados como los números de Reynolds,
Nusselt y Prandtl se derivan de consideraciones dimensionales así
como funcionales.

[Link]átegui R 7
• La turbulencia consiste en un conjunto
de remolinos de varios tamaños que
coexisten en la corriente de flujo.
• Continuamente se forman remolinos
grandes, que se rompen en otros más
pequeños, que a su vez se transforman
en otros todavía menores. Los
remolinos más pequeños desaparecen.
• El diámetro del remolino más pequeño
es de 10 a 100 μm. Los remolinos
pequeños se destruyen.
• Debido a su gran tamaño de los
remolinos. el flujo dentro de un
remolino es laminar, los remolinos más
pequeños contienen cerca de 1012
moléculas agrupadas.
• Los remolinos son de tamaño
macroscópico, es decir el flujo
turbulento no es un fenómeno
molecular.

[Link]átegui R 8
• El flujo turbulento es muy importante en
muchas áreas de la ingeniería, la naturaleza
de la turbulencia ha sido objeto de una
investigación profunda.
• Las mediciones de fluctuaciones de velocidad
de los remolinos en un flujo turbulento han
ayudado a explicar la turbulencia.
• Para el flujo laminar existen soluciones
exactas, mientras que para el flujo turbulento
no existen soluciones exactas.
• La turbulencia puede generarse por contacto
entre dos capas de fluidos que se desplazan a
diferente velocidad o por una corriente de
flujo en contacto con un sólido como una
pared o una esfera.
• En flujo turbulento, la velocidad fluctúa en
todas direcciones. La siguiente gráfica
muestra la variación de la velocidad
instantánea Vx en dirección x en un punto
determinado del flujo turbulento

[Link]átegui R 9
[Link]átegui R 10
10.1Experimento de Reynolds
Se usa para demostrar las diferencias cualitativas entre el flujo laminar
y turbulento

colorante

Experimento 1. Patrón de flujo a velocidad media baja con

inyección de colorante: régimen laminar

[Link]átegui R 11
[Link]átegui R 12
 colorante

Experimento 2. Patrón de flujo a velocidad superior con inyección

de colorante: régimen intermedio

[Link]átegui R 13
[Link]átegui R 14
[Link]átegui R 15
 colorante

Experimento 3. Patrón de flujo a alta velocidad con inyección de

colorante: régimen turbulento

[Link]átegui R 16
Conclusiones:

1) Sobre cierta velocidad promedio para un sistema dado se forman

remolinos relativamente grandes que fluyen mediante una
conducta al azar.
2) Los remolinos son mayores y abundantes en el centro del tubo

3) Un aumento en la velocidad promedio del fluido amplía la zona

turbulenta hasta que el tubo se llena de remolinos.

[Link]átegui R 17
10.1Experimento de Reynolds
Se usa para demostrar las diferencias cualitativas entre el flujo laminar y
turbulento

Video: [Link]

[Link]átegui R 18
 N
Número de Taylor

R1
[Link]átegui R
R2 19
[Link]átegui R 20
 Glicerina Glicerina + Colorante Glicerina + colorante

μexp = 1.22 pa.s

N= 5 rps
Tiempo = 30 seg
ρ = 1210.8kg/m3 Ta= 19 (< 1700)

[Link]átegui R 21
 Glicerina Glicerina + Colorante Glicerina + colorante

μexp = 1.24 pa.s

N= 50 rps
Tiempo = 30 seg
ρ = 1210.8kg/m3

Ta= 1877 (> 1700)

[Link]átegui R Hay mezcla 22
 Flujo de Fluidos
[Link]
[Link]átegui R 23
10.2 Ecuaciones de Von Karman: Flujo turbulento

po pL
R Zona laminar

Zona intermedia

r
Zona turbulenta

FLUIDO x

[Link]átegui R 24
[Link]átegui R 25
 Para la zona laminar
U   y  , (1)
Donde : y   [0  5]
El parámetro (y ) es una posición generalizada para un fluido en
movimiento en una posición radial (r) particular.
El parámetro U  es una velocidad generalizada :
 V( r )
U  *
, ( 2)
U
Velocidad de fricción: U *
 Rx g c
U*  , (3)

 ( R  r ) U *
y  , ( 4)


[Link]átegui R 26
 Para la zona intermedia
U   3,05  5Ln ( y  ), (5)
Donde : y   [5  30]
Para la zona turbulenta
U   5,5  2,5 Ln ( y  ), (6)
y   [30  valor máximo (r  0)]
Definiciones adicionales :
R : radio del tubo
V(r) : velocidad puntual en r
r : posición radial
 : densidad del fluido
 .viscosidad
 Rx : tensión del fluido en la pared
[Link]átegui R 27
• Distribución de velocidad en una tubería, con flujo totalmente
desarrollado de un fluido newtoniano, para flujo laminar y turbulento
hasta NRe=10000

[Link]átegui R 28
 Distribución universal de velocidad, flujo turbulento de un fluido newtoniano
en un tubo liso.

[Link]átegui R 29
 Análisis matemático
Análisis matemático Acción combinada del flujo turbulentoy laminar puede
representarse de la siguiente manera :
Transporte molecular :
d V x ( r )
( rx g c ) m   , (7 )


dr
Transporte Turbulento :

( rx g c ) t   ET
dVx (r )
dr
, (8) 
( rx g c )T  ( rx g c ) m  ( rx g c ) t
dVx (r )
 rx g c  (  ET ) , (9)
dr
La difusivilidad de remolino o turbulenta varía desde cero correspondiente
a la pared ( flujo laminar) a valores elevados en la zona turbulenta

[Link]átegui R 30
 Balance de energía :
Caída de presión : P  p o - p L
P.S
S  R 2
Fuerza de resistencia debida a la tensión del fluido : Rx
 Rx . A
A  2RL
P.S   Rx . A, (10)
P
 Rx  ( ) R, (11)
2L
Para una posición cualquiera r :
P
 rx  ( )r , (12)
2L
r
 rx   Rx ( ), (13)
R

[Link]átegui R 31
Ejemplo: A través de un tubo de 2” DI, fluye agua con una caída de presión de 10 lbf/pie2. por pie de
tubo.
a) Determinar la difusividad de remolino en y+ 5, 30 y 600. en pie2/s
b) Determinar el espesor de la zona de transición. en pie
c) Determinar la difusividad máxima en pie2/s, así su posición (r).
d) Completar el siguiente cuadro y graficar (y+ -r), (y+ -ET)

y+ r/pie ET/pie2/s
0
5
10
20
30
50
100
200
400
800
[Link]átegui R 32
Cálculo de la tensión en la pared del tubo
(P ) R (10lbf / pie2 )(1 / 12 pie)
 Rx    0,417lbf / pie2
2L 2 x1 pie
Velocidad de fricción: U *
 Rx .g c 0,417x32,2
U*    0,464 pie / s
 62,4
Velocidad generalizada
Vr Vr
U  
U * 0,464
ANALISIS DE LA Ec. VON KARMAN
ZONA LAMINAR :
U   y
y   [0 - 5]
y   0  Vr  0
y   5  Vr  2,32pie/s
Espesor de la zona laminar :
y 5 x6,72x10 4
(R - r) y  5    1,16x10 4 pie
U * . 0,464x62,4

[Link]átegui R 33
ZONA de TRANSICION :
U   -3,05  5Ln(y )
y   [5 - 30]
y  5  Vr  2,32pie/s
y   30
U   3,05  5Ln (30)  13,96
Vr  13,96 * 0,464  6,48 pie / s
Espesor de la zona transición :
y 30x6,72x10 4
(R - r) y  30    6,96x10 4 pie
U * . 0,464x62,4
(R - r) y  30  (R - r) y  5  5,80x10 4 pie

[Link]átegui R 34
 ZONA de TRANSICION :
La difusivilidad de remolino
du  5 Transporte Turbulento :

  , (a)
dy y
dVx (r )
dV
du   *r , (b)
 rx g c  (  ET )
u dr
*
 dVx (r )
dy   (
u
)dr, (c)  rx g c    (  ET )
 dr
Dividiendo b y c :
ET , y   [ 5....]
du   dVr
 [ * 2 ] , (d )
dy 
(u )  dr
(a  d )
dVr 5(u * ) 2 5( rx g c )
  , (e)
dr y 
y 

[Link]átegui R 35
 En la Ec9 :
r 5 g
( Rx g c )  [  ET ]{ Rx  c }, ( f )
R y
ry 
  ET  ....
5 R
ry 
ET    .....(14)
5 R
Cálculo de la difusividad en y   5
ZONA LAMINAR
6,72 x10  4
  1,077 x10 5 pie 2 / s
62,4
Cálculo de la difusividad en y   5 Ec.14 :
E T  1,077 x10 5 pie 2 / s
Cálculo de la difusividad en y   30 : Ec..14
ry  0,082637 * 6,72 x10  4 * 30
  6,41x10 5 pie 2 / s
5 R 5 * (1 / 12) * 62,4
ET  5,33 x10 5 pie 2 / s
[Link]átegui R 36
Empleando la Ec.6 se obtiene :
ry 
  ET 
2,5 R
ry 
ET    ...(15)
2,5 R
y   30
r  0,082637pie
ET  11,74 x10 5 pie 2 / s
Cálculo de (E T )max en r  R/2
E T , max  385.59 x10 5 pie 2 / s
Cálculo de la difusividad en el centro del tubo : Ec.15
r0
y   3590
u   25.96
Vr 0  12,05 pie / s
E T  -1,08x10
[Link]átegui R
-5
pie 2
/s 37
Ejemplo: un fluido seudoplástico que obedece la ley exponencial y,
cuya densidad es de 961 kg/m3, fluye a través de un tubo circular liso
cuyo diámetro interior mide 0,0508m , a velocidad promedio de
6,10m/s . Las propiedades de flujo del fluido son n’=0,30 y
k’=2,[Link]’ /m2. Calcule la caída de presión por fricción en un tubo
de 30,3 m de largo.

[Link]átegui R 38
Solución:
Datos :
K '  2,744
n'  0,30
D  0,0508m
V  6,10m / s
  961kg / m 3
L  30,5m
Parámetro: Viscosidad generalizada
  K '*8n '1 , N .s n ' / m 2 ( SI )
  g c K '*8n '1 , lbm / pie.s 2 n ' ( sistema inglés)
D n 'V 2 n '  D n 'V 2 n '  (0.05080,30 )(6,101,17 )(961)
N Re, gen   
 K '*8n '1 2,744* (80, 70 )
N Re, gen  1,328x104 ( flujo turbulento)
Con : Nre, gen y n' , se obtiene : f  0,0032
L V2 30,5 (6,102 )
p  4f  4(0.0032)(961)( )  137,4 N / m 2
D 2 0,0508 2
p  2870lbf / pie2

[Link]átegui R 39
 TRANSFERENCIA DE CALOR
EN
REGIMEN TURBULENTO

[Link]átegui R 40
 Análisis matemático
dVx ( r )
 rx g c  (  ET ) , (9)
dr
La ecuación general para transferencia de calor en régimen
turbulento puede escribirsecomo :
d
  -(  E T ) , (16)
d
La densidad de flujo de calor total puede escribirsecomo :
Q Q Q
( )  ( ) c  ( ) T , (17)
A A A
q  qc  qT (17' )
Q dT
( ) c  k , (18)
A dy
Q d (CpT )
( ) c  (k / Cp ) , (18)
A dy
Q dT
( ) T  CpEH , (19)
A dy
Q d (CpT )
( ) T   EH , (19)
A dy
[Link]átegui R 41
 Balance de energía :
Donde :
Caída de presión : P  p o - p L
E H  Difusividad de remolino para P.S
transmisión de calor. S  R 2
Fuerza de resistencia debida a la tensión del fluido : Rx
Q dT
( )  - Cp[  E H ] , (20)  Rx . A
A dy A  2RL
Para la transferencia de calor es posible escribir P.S   Rx . A, (10)
P
una ecuación similar a la Ec.13 : [Link]  Rx  ( ) R, (11)
2L
r Para una posición cualquiera r :
q r  q R ( ), (21) P
R  rx  ( )r , (12)
2L
r dT r
q R ( )  - Cp[  E H ] , (22)  rx   Rx ( ), (13)
R
R dr

[Link]átegui R 42
 Balance de energía :
Transferencia de masa en flujo turbulento : Caída de presión : P  p o - p L
P.S
N dC
( A )  -[DAB  E M ] A , (23) S  R 2
A dy Fuerza de resistencia debida a la tensión del fluido : Rx
 Rx . A
Para un sistema cilindrico: A  2RL
P.S   Rx . A, (10)
r
 r   R ( ), (24) P
 Rx  ( ) R, (11)
R 2L
Para una posición cualquiera r :
r dCA P
( A ) R ( )  -[DAB  E M ] , (25)  rx  (
2L
)r , (12)

R dr  rx
r
  Rx ( ), (13)
R

[Link]átegui R 43
 Transporte Turbulento : Transferencia de masa en flujo turbulento :
NA dC
dVx (r ) ( )  -[DAB  E M ] A , (23)
 rx g c  (  ET ) A dy
dr Para un sistema cilindrico:
dV (r ) r
 rx g c    (  ET ) x  r   R ( ), (24)
R
dr
r dCA
( A ) R ( )  -[DAB  E M ]
ET , y   [ 5....] R dr
, (25)
Donde :
E H  Difusividad de remolino para
transmisión de calor.
Q dT
( )  - Cp[  E H ] , (20)
A dy
Para la transferencia de calor es posible escribir
una ecuación similar a la Ec.13 : [Link]
r
q r  q R ( ), (21)
R
r dT
q R ( )  - Cp[  E H ] , (22)
R dr

[Link]átegui R 44
 Ejercicio:. En una tubería de 2 pulg de diámetro, la velocidad r=R/2 es 40 pies /s,
para el aire a 80°F y 1 atm de presión. A esta velocidad existe un flujo turbulento.
Datos: ν=16.55x10 -5 pie 2/s , ρ = 0.0735 lb/pie 3.
Calcular:
a) La velocidad de fricción en pie/s.
b) El esfuerzo de corte en la pared (τRx ) y en r = R/2 en lbf/pie 2.
c) La caída de presión en 10 pies de tubería
d) El espesor de la zona turbulento y el porcentaje que representa con respecto al
diámetro..
e) La difusividad de remolino Et en y+=30 y en r=R/2.

ΔP
R

[Link]átegui R L 45
Para la zona laminar
r
U   y  , (1)
 rx   Rx ( ), (13)  rx  ( P )r , (12)
Donde : y   [0  5] R 2L
El parámetro (y  ) es una
posición generalizada
para un fluido en
movimiento en una Para la zona intermedia
posición radial (r) particular.
U   3,05  5 Ln ( y  ), (5)
El parámetro U  es una
velocidad generalizada :
Donde : y   [5  30]
 V( r ) Para la zona turbulenta
U  , ( 2)
U   5,5  2,5 Ln ( y  ), (6)
*
U
Velocidad de fricción: U *
y   [30  valor máximo (r  0)]
 Rx g c
U*  , (3) Definiciones adicionales :

 ( R  r ) U *
y  , ( 4)

[Link]átegui R 46
r  R/2
RU *
RU * X Fx
y   1.9 0.126806232
2 2(  /  )
1.91 0.00345964
1.9105 -0.002675596
Para la zona turbulenta
(Vr / U *)  5,5  2,5 Ln ( RU * /(2v))
1.9103 -0.000221866
F ( x)  0 1.9104 -0.001448792
U*  x
Fx  (Vr / x)  5,5  2,5 Ln ( Rx /(2v))  0
x...Fx.... Fx´....Fx / Fx´......( x  Fx / Fx´) P
 rx  ( )r , (12)
Vr  40 pie / s 2L
R  1 / 12 pie
v  16.55x10^ 5 pie^ 2 / s

[Link]átegui R 47
10.3 Fundamentos de los mecanismos de transferencia
A) Transferencia turbulenta: Cantidad de Movimiento

En el régimen laminar para un cilindro se tiene :

_
r
Vr  2 V[1  ( ) 2 ], (26)
R

dV 4V
 ( 2 )r , (27)
dr R
Ec.27 en la ecuación de la ley de la
viscosidad de Newton :

4 V
( ry g c ) r  (
2
)r , (28)
R
Tensión en la pared de la tubería :
 
4 V 8 V
( ry g c ) r  R  ( )( ), (29)
R D

[Link]átegui R 48
La Ec.29 establece que la transferencia
total de [Link]. en la pared, es igual a
la rapidez de [Link] a través del

8V
fluido por transferencia molecular ( )
D
• La Ec.29 establece una correlación
para dos mecanismos , la transferencia total
a la pared y la transferencia por difusión
[Link] mecanismos son iguales en
el flujo laminar, debido a que ningún otro
mecanismo actúa en el sistema.

[Link]átegui R 49
 Balance de fuerzas :
 Ry A  P.S , (30)
A  DL
S  D 2 / 4
P  po  pL
PDg c
( ry g c ) r  R  , (31)
4L
( Ec.31  Ec.29) :

PDg c 32V 
(4 ry g c ) r  R  ( ), (32)
L D
[Link]átegui R 50
• A medida que la velocidad media aumenta en el sistema, a cierta
velocidad, la transferencia de cantidad de movimiento se iniciará
mediante el mecanismo turbulento (Nr>2100).
• En un sistema en el cual está presente el flujo turbulento, el
mecanismo de éste debe combinarse con los dos representados
Ec.29.
2
mV y
[Link]ética( Ec ) 
2
V y
2
[Link]ética Ec
 
volumen V 2
 

Ec V 2
 , (33)
V 2
Flujo laminar :   0,50
Flujo tapon :   1,0

[Link]átegui R 51
 Debido a que la actividad de los remolinos es errática,
puede suponerse que la energía cinética de todos los
remolinos ( energía cinética media) dirigida en el sentido (y)

V 2
es proporcional a ( ).
2

V 2
El término ( ) es una descripción matemática global
2
del mecanismo de turbulencia y tiene las mismas unidades que
la Ec.32.

V 2
Al dividir la Ec.32 por ( ) se obtiene :
2
[Link]átegui R 52
 (8 ry g c ) r  R 2PDg c 
[ -
][ 
]  64[ 
], (34)
V 2
L V 2
D V
Cada término de esta ecuación, es una relación adimensional de dos mecanismos :
El término de la izquierda representa :
Transferencia total de [Link]/ Transferencia de [Link]. turbulento

Además define el factor de fricción: f

(8 ry g c ) r  R 2PDg c
f 
 
, (35)
V 2 V L 2

El término de la derecha representa :

Transferencia total de [Link] molecular/ Transferencia de [Link]. turbulento
y define el número de Reynolds (N R )
Régimen laminar :
64
f , (36)
NR
 
D V DQ DV
NR    , (37)
  v
f  F (NR )

[Link]átegui R 53
 Factores de fricción para el flujo en tubos

[Link]átegui R 54
 Expresión matemática Ecuación de Barr:
Ecuación de Colebrook-White:
1  K / D 5.1286
 2 Log10 ( ) 0.89 
1 K / D 2.51  f  3.7 N 
 2 Log10  
R

f  3.7 NR f 
Ecuación de Haaland:
Para el caso particular de tuberías lisas
la rugosidad relativa, es decir la
relación entre la rugosidad en las
1  K / D 1.1 6.9 
paredes de la tubería y el diámetro de  1.8Log10 ( )  
la misma, es muy pequeño con lo que
el término k/D es muy pequeño y
f  3.7 N R

puede despreciarse el primer sumando

situado dentro del paréntesis de la
ecuación anterior. Quedando en este k/D = rugosidad relativa total = D/ɛ
caso particular la ecuación del
Re = Número de Reynolds
siguiente modo:
f = factor de fricción
k/D0 D = diámetro interno de la cañería
1
 2 Log10 ( N R f )  0.8
f
[Link]átegui R 55
El diagrama de Moody es la representación gráfica en escala doblemente logarítmica del factor de
fricción en función del número de Reynolds y la rugosidad relativa de una tubería, diagrama hecho por
Lewis Ferry Moody.

[Link]átegui R 56
[Link]átegui R 57
[Link]átegui R 58
[Link]átegui R 59
V = 600 L/min = 0.0100 m3/s
Tubería acero comercial:
e = 0.000046 m DN 4", cédula 40 L = 28 m
D = 102.2604 mm A = 0.008213 m2 v = V / A = 1.2176 m/s

NRe = 879 * 1.2176 * 0.10226 / 0.000652 = 167861.73 Turbulento

D / e = 0.10226 / 0.000046 = 2223.0522

Ecuación de Darcy:
0,25
f 
2 L
2
  1 5,74 
log  
0,9 
Ftub  f
  3,7( D / e ) N Re  2g D

f = 0.018993 Ftub = 0.39337 m

[Link]átegui R 60
ECUACIONES DE FLUJO : • Ejemplo: Calcular la caída de
1 K / D 2.51  presión a través de 50pies de tubo
Flujo laminar :  2 Log10   
f  3.7 NR f  liso, para un aceite que fluye a
64
f , (38) una velocidad media de 8pie/s.
NR • Datos: D=4”, μ=5cp ,ρ=60
Flujo turbulento : lb/pie3

0,184
f 0,20
, (39) 1
 2 Log[ N R f ]  0,80
(N R ) f
5 x103  N R  2 x105 D V  
NR 
f  5,6 x103 
0,50
, (39) 
0,32
(N R ) 2PDg c
f 
3 x103  N R  3 x106  V  2 .L
1
 2 Log[ N R f ]  0,80, (40) F(x) = 0
f x =ΔP

[Link]átegui R 61
 D V   (4 / 12 pie) * (8 pie / s ) * (60lb / pie3 )
NR    47622.05
 3
3.360x10 lb / pie.s
2PDg c P(2 * (4 / 12) * 32.2) 5
f   5 .590 x10 P
 V  .L
2 2
8 * 60 * 50
NewtonRaphson : F ( x)  0 F ( x) 
1
 2 Log[ N R f ]  0,80
f
133.7468
 2 Log (356.052 P )  0.8
P
133.7468
F ( x)   2 Log (356.052 x )  0.8
x
X ......... F .......... F´......... F / F ´....... X  F / F´

[Link]átegui R 62
[Link]átegui R 63
• Factor de fricción de Fanning en función del número de Reynolds generalizado para fluidos
newtonianos y no newtonianos independientes del tiempo, fluyendo en tuberías lisas. Tomado de
D.W Dodge y A. B Metzner, A.I Ch. E. J, 5 189(1958). Reproducido con permiso.

[Link]átegui R 64
 Coeficiente de Transferencia para calor y masa en régimen turbulento :
Ecuación general :
d
  -[  Et] , (41)
dx
 : flujo de una propiedad
 : difusivilidad molecular(/ )
Et : difusibilidad de remolino
 : concentración de la propiedad
Sistema radial :
d
  -[  Et] , (42)
dr
Ademas :
 r   R (r / R ), (43)
d
 R (r / R )  -[  Et] , (44)
dr
Et  f (r )

[Link]átegui R 65
 r   R (r / R ), (43)
d
 R (r / R )  -[  Et] , (44)
dr
Et  f (r )
rdr
 R 0   R  R 0 d , (45)
R

[  Et]
 R  R 0 d   R ( R   0 ), (47)
 0 : Concentración en el centro del tubo
-
 : valor medio de la propiedad
- 1 VR
  0 dV , (48)
VR
V : flujo v olumétricoen un punto
VR : flujo v olumétricototal
dV  [Link]
dV  2[Link]
[Link]átegui R 66
R - 0  diferencia máxima de la -
4[  Et ]
propiedad transferente
e , (53)
R ( R   )

D
 R   R  0 d   , (49)
 Ec.53 en Ec.52
De la Ec. 45 

 R 0R
rdr

R
 0 rdr
R ( R A) R  e ( R   ) AR , (54)
- -
[  Et ] [  Et ] 
(R   )

R 2 R
, (50) ( R A) R   , (55)
2[  Et ]
- 1
Ec.50  Ec.49 eAR
R R 
-

 ( R   ), (51) ( R A) R : Rapidez de transferencia
2[  Et ] 
- (R   ) : Potencial de transferencia
- 4[  Et ] 
R  ( R   ), (52)
D 1
e : coeficiente de transferencia
: Resistencia de transferencia
eAR
[Link]átegui R 67
 Coeficiente de transferencia de calor en flujo turbulento en tubos
Ec.55 para transferencia de calor :

h 4(  Eq )
 eq  , (56)
c p  qD
  Q
  c pT Ts
A  Q
 
TR
Er  Eq
_ r _
[( c pT ) R  ( c pT )] T
QR   , (57)
c p
[ ]
hAR x
_
[(T ) R  (T )]
QR   , (58)
1
[ ]
hAR
LEY DE ENFRIAMIEN TO DE NEWTON
_
QR  [(T ) R  (T )]hAR , (59)
[Link]átegui R 68
 LEY DE ENFRIAMIEN TO DE NEWTON
Q  hA(Ts  T )

Aire :  ,  , T , h

Ts


[Link]átegui 69
Los mecanismos de transferencia de calor dependen de las condiciones de flujo.
Las condiciones de flujo dependen de la transferencia de cantidad de movimiento,
por lo que se establecen dos mecanismos: [Link] Molecular y turbulento.

Transferencia de calor :
Número de Nusselt : N Nu
hD
N Nu  , (60)
k
Número de Reynolds : N R
-
D V
NR  , (61)

Número de Prandtl : Npr
 cp
Npr   , (62)
 k
[Link]átegui R 70
FLUIDO QUE SE CALIENTA
N Nu  0,023N R
0 ,80
Npr 0,40 , (63)
FLUIDO QUE SE ENFRIA
N Nu  0,023N R
0 ,80
Npr 0,30 , (64)
LAS PROPIEDADE S DEL
FLUIDO SE EVALÚAN A
TEMPERATUR A PROMEDIO
Ecuación mejorada, e independiente
de si el fluido se calienta o enfría :

 0,14
N Nu  0,023N R
0 ,80
Npr 1/3 ( ) , (65)
R
N R  104
0,50  Npr  100

[Link]átegui R 71
[Link]átegui R 72
Ejemplo: A través de una tubería de 2” de diámetro interno y 10 pies de
longitud, fluye un fluido a razón de 400gal/min. La temperatura promedio es
de 120ºF. Calcular el coeficiente de transferencia de calor y el flujo de calor
si la temperatura de la pared del tubo es (a) 50ºF y (b) 170ºF.
A la temperatura promedio de 120ºF:
ρ=56,2 lb/pie3, cp = 0,465 btu/lb.ºF, k=0,070 Btu/[Link].ºF

T /ºF µ(cp)
50 100
120 20
170 6

LEY DE ENFRIAMIEN TO DE NEWTON

Q  hA(Ts  T )

[Link]átegui R 73
Solución :
  56,2 lb/pie3, cp  0,465 btu/lb.º F, k  0,070 Btu/[Link].º F
 Flujo 400 /(7.42 * 60)
V    41.18Pie / s
Area  (2 / 12) 2 / 4
Q  hA(Ts  T)
Número Reynolds : NR
 Q  194.68 * ( (2 / 12) *10) * (50  120)  71353.95Btu / h
D V
(2 / 12) * (41.18) * (56,2)
NR    28699 Enfriamiento!!!
 20 * 6.72x10-4
Transferencia de calor :
Número de Nusselt : N Nu
hD (2 / 12)h
N Nu    2.381h
k 0.070
Número de Prandtl : Npr

 cp 
0.465* (20 * 2.42)
Npr     321.514
 k 0.070

 20
  0.20
 R 100


N Nu  0,023N R
0 ,80
Npr 1/3 ( ) 0,14
R
2.381h  0.023* (28699) 0.80 * (321.514)1 / 3 * (0.20) 0.14
h  194.68Btu / h. pie 2 . F

[Link]átegui R 74
Ejemplo: A través de una tubería de 4” de diámetro interno y 20 pies de
longitud, fluye un fluido a razón de 450gal/min. La temperatura promedio es
de 120ºF. Calcular el coeficiente de transferencia de calor y el flujo de calor si
la temperatura de la pared del tubo es 170ºF.
A la temperatura promedio de 120ºF:
ρ=56,2 lb/pie3, cp = 0,465 btu/lb.ºF, k=0,070 Btu/[Link].ºF

T /ºF µ(cp)
50 100
120 20
170 6

LEY DE ENFRIAMIEN TO DE NEWTON

Q  hA(Ts  T )

[Link]átegui R 75
Ejercicio: Se desea calentar agua a razón de 50 lb/min, de 100 a 130°F en
un tubo de 2 pulg de diámetro, cuya superficie está a 220°. ¿qué longitud
deberá tener la tubería?
A la temperatura promedio Tp=115°F, las propiedades son:
K=0.366 Btu/[Link].°F, Cp=0.997 Btu/lb.°F, µ=1.57 Lb/pie.h, ρ=61.9
lb/pie3.

Transferencia de calor :
Número de Nusselt : N Nu
hD
N Nu  , (60)
k
Número de Reynolds : N R
-
D V
NR  , (61)

Número de Prandtl : Npr
 cp
Npr   , (62)
 k
COEFICIENTE DE TRANSFERENCIA DE MASA EN FLUJO
TURBULENTO DENTRO DE TUBOS.

De la Ec 53 y 55 :
  Dm
_ _
Er  E N
A  N a Coeficiente de transferencia de masa :
  ca _
4( D  e N )
K'c  e N  , (66)
ND

[Link]átegui R 77
Número de Schmidt : N Sc
 
N Sc   , (67)
Dm Dm
' Dm
kc  0,023N 0,83
Re N 0 , 33
Sc ( ), (68)
D
Para gases :
'

( ) R  K c ( paR  pa ), (69)
_
Na
A RT
_
Na
( ) R  k G ( paR  pa ), (70)
'

A
[Link]átegui R 78
• Ejemplo: Una columna de pared húmeda de 2pulg de diámetro interno contiene aire y
CO2 que fluyen a 3pies/s. La concentración de CO2 en el aire en un punto de la columna
es 0,1 fracción molar. En el mismo punto de la columna, la fracción de CO2 en el agua
en la interfase agua –aire es de 0,005 fracción molar. Determinar el coeficiente de
transferencia de masa y el flujo másico, sabiendo que la columna opera a 10 atm.y 25ºC.
• Dato: Dm=0,164cm2/s (1atm).

[Link]átegui R 79