FTF 2021

FÓRMULAS, TABLAS Y FIGURAS DE
TRANSFERENCIA DE CALOR
Juan Carlos Ramos González

Doctor Ingeniero Industrial
Raúl Antón Remírez

Doctor Ingeniero Industrial
Julio de 2020
Transferencia de Calor / Curso 2020-21
Fórmulas, Tablas y Figuras
ii
ÍNDICE
Fórmulas
Tema 1. Introducción a la transferencia de calor y a la conducción ...................................... 1

Tema 2. Conducción unidimensional en régimen estacionario ............................................. 3
Tema 3. Conducción bidimensional en régimen estacionario ............................................... 6
Tema 4. Conducción en régimen transitorio .......................................................................... 6
Tema 5. Introducción a la convección ................................................................................. 10
Tema 6. Convección forzada en flujo externo ..................................................................... 12
Tema 7. Convección forzada en flujo interno...................................................................... 15
Tema 8. Convección libre o natural ..................................................................................... 18
Tablas y Figuras
Tema 2. Conducción unidimensional en régimen estacionario ........................................... 20

Tema 3. Conducción bidimensional en régimen estacionario ............................................. 23
Tema 4. Conducción en régimen transitorio ........................................................................ 30
Tema 6. Convección forzada en flujo externo ..................................................................... 33
Tema 7. Convección forzada en flujo interno...................................................................... 38
Tema 8. Convección libre o natural ..................................................................................... 42
Tablas de propiedades termofísicas y de funciones matemáticas ........................................ 45
Alfabeto griego .................................................................................................................... 52
i
ii
TEMA 1. INTRODUCCIÓN A LA TRANSFERENCIA DE CALOR Y A LA

CONDUCCIÓN
 Calor o transferencia de calor o velocidad de transferencia de calor: q [J/s = W].

 Flujo calorífico o de calor: q  [W/m2].
dT
 Ley de Fourier: q x   k . q x  q x  A . En condiciones de régimen estacionario y con una
dx
dT T T T  T2 T
distribución lineal de temperaturas: q x  k  k 2 1  k 1 k .
dx L L L
 Conductividad térmica: k [W/mꞏK].
 Ley de enfriamiento de Newton: q x  h (T s  T ) .
 Coeficiente de transferencia de calor por convección: h [W/m2ꞏK].

 Potencia emisiva superficial: E [W/m2].
 Ley de Stefan-Boltzmann para un cuerpo negro: E b  Ts4 .
 Constante de Stefan-Boltzmann:  = 5,67ꞏ10-8 W/m2ꞏK4.

 El flujo de calor emitido por una superficie real a la misma temperatura que un cuerpo negro
siempre será menor y viene dado por: E  Ts4 , donde  es la emisividad, que puede variar
entre 0 y 1.
 Se llama irradiación, G, a la velocidad con la que la radiación incide sobre un área unitaria.
La proporción de la irradiación total que es absorbida por la superficie viene dada por la
absortividad,  (01), según la siguiente expresión: G abs  G . Irradiación de los
alrededores: G  Talr4 .
 Intercambio de radiación para una superficie gris y difusa ( = ):
 
q rad
q
A

 E b (Ts )  G   Ts4  Talr4 .  También se puede expresar como:
q rad
  hrad (Ts  Talr ) , siendo hrad el coeficiente de transferencia de calor por radiación:
 
hrad   ( Ts  Talr ) Ts2  Talr2 .
 Principio de conservación de la energía en un volumen de control formulado en un instante de

dE
tiempo (t): E ent  E gen  E sal  alm  E alm .
dt
1
 Principio de conservación de la energía en un volumen de control formulado en un intervalo

de tiempo (t): E ent  E gen  E sal  E alm .
 Principio de conservación de la energía en una superficie de control: E ent  E sal  0 .
   T  T  T    
 Ley de Fourier vectorial: q   kT   k  i  j k   i q x  j q y  k q z .
 x y z 
 Capacidad térmica volumétrica:  cp [J/m3ꞏK]. Mide la capacidad de un material para

almacenar energía térmica.
k
 Difusividad térmica:   [m2/s]. Mide la capacidad de un material para conducir energía
c p
térmica en relación con su capacidad para almacenarla.
 Ecuación de difusión de calor en coordenadas cartesianas:
  T    T    T  T  W 
k    k    k   q  c p .
x  x  y  y  z  z  t  m 3 
T
 Ecuación de difusión de calor vectorial: ꞏ(kT )  q  c p .
t
 En el caso de transmisión unidimensional en régimen estacionario y sin generación de
d  dT 
energía: k   0 . Teniendo en cuenta la ley de Fourier ( q x   k dT dx ), esta ecuación
dx  dx 
implica que el flujo de calor en la dirección de transmisión es una constante
( dq x / dx  0  q x  cte. ).
 Ecuación de difusión de calor en coordenadas cilíndricas (r radial,  angular o longitud, z

axial, elemento diferencial de volumen: drꞏrdꞏdz):
1   T  1   T    T  T
 kr  2  k    k   q  c p .
r r  r  r     z  z  t
 Ecuación de difusión de calor en coordenadas esféricas (r radial,  polar, cenital o colatitud, 

azimutal o longitud, elemento diferencial de volumen: drꞏrsendꞏd):
1   2 T  1   T  1   T  T
 kr  2  k   2  ksen   q  c p .
r r 
2
r  r sen      r sen   
2 2
  t
 Condición de contorno de primera clase o de Dirichlet: superficie mantenida a temperatura

constante, T(x = 0, t) = Ts.
2
 Condición de contorno de segunda clase o de Neumann: flujo de calor fijo o constante en la

T
superficie, q s( x  0)  k . Un caso especial es la superficie perfectamente aislada o
x x 0
T
adiabática,  0.
x x 0
 Condición de contorno de tercera clase o de Robin: corresponde a la transferencia de calor por

 , superficie  q conv
convección en la superficie, q cond  . Si el fluido está en contacto con la superficie
T
de la pared donde está el origen de coordenadas:  k  hT  T ( x  0, t ) . Si el fluido
x x 0
está en contacto con la superficie de la pared opuesta al origen de coordenadas:

T
k  hT ( x  L, t )  T  .
x xL
TEMA 2. CONDUCCIÓN UNIDIMENSIONAL EN RÉGIMEN ESTACIONARIO
Ts1  Ts 2 L
 Resistencia térmica de conducción para pared plana: Rt ,cond   .
qx kA
Ts  T 1
 Resistencia térmica de convección: Rt ,conv   .
q hA
Ts  Talr 1
 Resistencia térmica de radiación. Rt ,rad   .
qrad hrad A
T 1
 Coeficiente global de transferencia de calor, U: q x  UA T . Rtot   Rt   .
q UA
 Ley de Fourier expresada en forma integral para un sistema general en condiciones de
régimen estacionario sin generación de calor y con conducción unidimensional (en este caso,
x dx T
la transferencia de calor, qx, es una constante independiente de x): q x     k (T )dT .
x0 A( x) T0
(Ts1  Ts 2 ) ln(r2 / r1 )
 Resistencia térmica de conducción para una pared cilíndrica: Rt ,cond   .
qr 2Lk
1 1
 Resistencia térmica de convección para una pared cilíndrica: Rt ,conv   .
Ah 2rLh
3
 Resistencia térmica de conducción para una pared esférica:

(Ts1  Ts 2 ) 1 1 1
Rt ,cond     .
qr 4k  r1 r2 
1 1
 Resistencia térmica de convección para una pared esférica: Rt ,conv   .
Ah 4r 2 h
 El coeficiente global de transferencia de calor en una pared cilíndrica o esférica depende del
área en función de la cual se exprese: U 1 A1  U 2 A2  U 3 A3  ...  U i Ai   Rt  .

1
E gen  W 
 Generación de energía térmica por unidad de volumen: q  e gen  .
Vol  m 3 
d 2T  1 dAc  dT  1 h dAs 
 Ecuación de calor para una aleta:    (T  T )  0 .
dx 2  Ac dx  dx  Ac k dx 
 Distribución de temperaturas y transferencia de calor para aletas de área de sección

transversal uniforme:
 Caso A, con transferencia de calor por convección desde el extremo de la aleta (
dT
hAc T ( L)  T   kAc ):
dx xL
 ( x) cosh m( L  x)  (h / mk )senh m( L  x) senh mL  (h / mk ) cosh mL

 qf  M
b cosh mL  ( h / mk )senh mL cosh mL  (h / mk )senh mL
siendo  ( x)  T ( x)  T ,  b  Tb  T , m 2 
hP
kAc
, M   
hPkAc  b , P el perímetro y
Ac el área transversal.
d
 Caso B, extremo adiabático (  0 ):
dx x L
 ( x) cosh m( L  x)
 q f  M tanh mL
b cosh mL
 Caso C, extremo con temperatura establecida ((x = L) = L):

 ( x) ( L /  b )senh mx  senh m( L  x) cosh mL   L /  b 
 qf  M
b senh mL senh mL
 Caso D, aleta muy larga (L   y L  0, aplicable si mꞏL > 2,65):
 ( x)   b e  mx qf  M   
hPkAc  b
4
 La efectividad de una aleta se define como la razón entre la transferencia de calor de la aleta y
qf
la transferencia de calor que existiría sin la aleta:  f  , siendo Ac,b el área de la
hAc ,b b
sección transversal de la base de la aleta. El uso de aletas sólo se justifica cuando f  2.

b
 Resistencia térmica de una aleta: Rt , f  .
qf
 Teniendo en cuenta la resistencia térmica de convección de la base de la aleta, Rt,b = 1/hAc,b,

Rt ,b
se puede expresar la efectividad como:  f  .
Rt , f
 La eficiencia o rendimiento de una aleta se define como la razón entre el calor real transferido
por la aleta y el calor que transferiría si estuviera toda ella a la temperatura de la base:
qf qf
f   , siendo Af la superficie total de la aleta.
q máx hA f  b
 Teniendo en cuenta la ecuación que define la resistencia térmica de una aleta, se puede
1
expresar ésta en función de su eficiencia: Rt , f  .
hA f  f
 Para el caso de una aleta recta de sección transversal uniforme y con su extremo adiabático se
M tanh mL tanh mL
tiene:  f   .
hPL b mL
 Se puede emplear la expresión de la aleta con extremo adiabático para una aleta con extremo
activo, empleando una longitud de aleta corregida de la forma Lc = L+(t/2) para aleta
rectangular y Lc = L+(D/4) para aleta de aguja. Esta aproximación es válida cuando (ht/k) o
(hD/2k) < 0,0625.
 Aletas de sección transversal no uniforme. En las expresiones de la distribución de
temperaturas, la transferencia de calor y el rendimiento o eficiencia de este tipo de aletas
aparecen las funciones de Bessel modificadas de primera y segunda clase de orden 0 y orden
1 (I0, I1, K0 y K1) cuyos valores están tabulados en la Tabla H. En la Tabla 2.1 se muestran las
expresiones del rendimiento para distintos tipos de aletas de sección transversal no uniforme.
qt qt
 Dispositivo de varias aletas. Eficiencia global de la superficie:  o   , siendo qt
qmáx hAt b
la transferencia total de calor de la superficie total, At, que es la asociada a la superficie de las
5
aletas, Af, más la parte expuesta de la base, Ab. Es decir, At  NA f  Ab , siendo N el número
total de aletas.
 Este rendimiento también se puede expresar en función del rendimiento de una sola aleta:
NA f
q t  Nq f  q b  N f hA f  b  hAb b  q t   o hAt b   o  1  (1   f ) .
At
b 1
 Resistencia térmica efectiva del dispositivo de aletas: Rt ,o   .
qt hAto
TEMA 3. CONDUCCIÓN BIDIMENSIONAL EN RÉGIMEN ESTACIONARIO
 Factor de forma de conducción para sistemas bidimensionales, S: q  Sk (T1  T2 ) . Se obtiene

de la Tabla 3.1.
1
 Resistencia de conducción bidimensional: Rt ,cond ( 2 D )  .
Sk
 MDF: Para obtener la ecuación de diferencias finitas de un nodo aplicando el principio de
conservación de la energía a un volumen de control alrededor del nodo se supone que todo el
flujo de calor es hacia el nodo. Como estamos en régimen permanente la ecuación a emplear
es: E ent  E gen  0 . El término de energía entrante puede incluir calores de conducción o de
Tm 1,n  Tm , n
convección que se evalúan con la ley de Fourier ( q ( m 1,n )( m , n )  k ( yꞏ1) ) o la ley
x
de enfriamiento de Newton ( q (  )( m ,n )  hxꞏ1(T  Tm ,n ) .
TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO
Rt,int U ꞏLc Tint

 Definición general del número de Biot: Bi    .
Rt ,ext k Text
hLc Rt ,cond
 Número de Biot para un sólido con convección: Bi   .
k Rt ,conv
 Longitud característica: Lc = Vol/As. Para una pared plana de espesor 2L sometida a

convección simétrica en su superficie  Lc = L, y para un cilindro largo o una esfera de radio
ro  Lc = ro.
 ꞏt
 Número de Fourier: Fo  .
L2c
6
hLc
 El método de la resistencia interna despreciable es aplicable cuando Bi   0,1 .
k
 Distribución de temperaturas temporal en un sólido en el que se puede aplicar el método de la
resistencia interna despreciable:
 (t ) T (t )  T  hAs   t 
  exp  t   exp   .
 ini Tini  T  Volc p   t 
1
 Constante de tiempo térmica:  t  ꞏVolc p  Rt ꞏC t , siendo Rt la resistencia a la
hAs
transferencia de calor por convección y Ct la capacidad térmica del sólido.
 La transferencia total de energía que tiene lugar desde un sólido en el que se puede aplicar el
método de la resistencia interna despreciable durante un tiempo t será:
t t
Q(t )   qdt  hAs   (t )dt  Q (t )  Volc p ini 1  exp t  t   U ini 1  exp t  t  .
0 0
 Q(t )  U (0  t )  U (t )  U (0) .

 Solución analítica aproximada con el primer término (aplicable cuando Fo > 0,2) de la
distribución de temperaturas en una pared plana de espesor 2L sometida a convección:
T ( x, t )  T 180 180
 * ( x*, t*)   C1 exp( 12 Fo) cos(  1 x*)   o * cos(  1 x*) , siendo
Tini  T  
T ( x  0, t )  T
 o * ( x*  0, t )   C1 exp( 12 Fo) la temperatura del plano medio (x* = x/L =
Tini  T
0). Los valores de C1 y 1 se obtienen de la Tabla 4.1.
 Transferencia total de energía en función del tiempo desde o hacia una pared plana de espesor
 180 
sen 1 
Q(t )     * (0, t ) , siendo U = U = c Vol(T -
2L sometida a convección:  1 o ini p ini
1
o
Uo
T) = Qmáx la energía interna inicial de la pared referida a la temperatura del fluido o la
máxima cantidad de energía que se podría transferir al fluido o desde el fluido, según éste esté
a menor o mayor temperatura que la pared.
distribución de temperaturas en un cilindro largo de radio ro sometido a convección:
7
T (r , t )  T
 * (r*, t*)   C1 exp( 12 Fo) J 0 ( 1r*)   o * J 0 ( 1r*) , siendo
Tini  T
To  T
 o * (r*  0, t )   C1 exp( 12 Fo) la temperatura del eje central (r* = r/ro = 0) y J0
Tini  T
la función de Bessel de primera clase de orden cero cuyos valores se encuentran en la Tabla
G.
 Transferencia total de energía en función del tiempo desde o hacia un cilindro largo de radio
Q(t ) 2 * (0, t )
ro sometido a convección:  1 o J 1 ( 1 ) , siendo J1 la función de Bessel de
Uo 1
primera clase de orden uno cuyos valores se encuentran en la Tabla G y Uo = Uini =
cpVol(Tini - T) = Qmáx la energía interna inicial del cilindro referida a la temperatura del
fluido o la máxima cantidad de energía que se podría transferir al fluido o desde el fluido,
según éste esté a menor o mayor temperatura que el cilindro.
distribución de temperaturas en una esfera de radio ro sometida a convección:
T (r , t )  T 1 180 1 180
 * (r*, t*)   C1 exp( 12 Fo) sen(  1r*)   o * sen(  1r*) ,
Tini  T  1r *   1r * 
To  T
siendo  o * (r*  0, t )   C1 exp( 12 Fo) la temperatura del eje central (r* = r/ro =
Tini  T
0).
 Transferencia total de energía en función del tiempo desde o hacia una esfera de radio ro
Q(t ) 3 *  180 180 
sometida a convección:  1  o3 sen(  1 )   1 cos(  1 ) , siendo Uo = Uini =
Uo 1    
cpVol(Tini - T) = Qmáx la energía interna inicial de la esfera referida a la temperatura del
fluido o la máxima cantidad de energía que se podría transferir al fluido o desde el fluido,
según éste esté a menor o mayor temperatura que la esfera.
8
 En un sólido semiinfinito la condición de frontera interior es T(x, t) = Tini y la condición

inicial es T(x, 0) = Tini. Las soluciones analíticas para tres condiciones de frontera exterior
son:
Condición de frontera Distribución de temperaturas
T ( x, t )  Ts  x 
 erf  
Temperatura superficial Tini  Ts  2 t 
constante: T (0, t) = Ts T k (Ts  Tini )
q s(t )  k 
x x 0 (t )1 / 2
Flujo de calor superficial 2qo (t /  )1 / 2  x 2  qo x  x 
T ( x, t )  Tini  exp    erfc  
constante: q s  q o k  4t  k  2 t 
T ( x, t )  Tini  x 
Convección superficial:  erfc   
T  Tini  2 t 
T
k  hT  T (0, t )   hx h 2t    x h t 
x x0  exp  2  erfc   
k k 2 t k 
   
  
donde la función gaussiana de error, erf (), y la función complementaria de error, erfc (w) =
1 – erf (w), son funciones matemáticas estándar cuyos valores se encuentran en la Tabla E.
 Conducción multidimensional. Para las geometrías multidimensionales de la Tabla 4.2, la
solución multidimensional se expresa como un producto de soluciones unidimensionales que
corresponden a un sólido semiinfinito, una pared plana de espesor 2L o un cilindro infinito de
T ( x, t )  T T ( x, t )  T T ( r , t )  T
radio ro: S ( x, t )  ; P ( x, t )  ; C (r , t )  .
Tini  T Sólido Tini  T Pared Tini  T Cilindro
semiinfini to plana infinito
 MDF en régimen transitorio. Expresión en diferencias finitas de la velocidad de variación de

Tmp,n1  Tmp,n
la energía almacenada: E alm  c pV
 .
t
 Criterio de estabilidad en el MDF explícito: el coeficiente asociado con el nodo de interés en
el tiempo anterior (coeficiente de Tmp,n ) debe ser mayor o igual que cero. Así se obtiene un
 t
valor límite para Fo  , del que se obtiene el máximo valor permisible de t.
( x ) 2
9
TEMA 5. INTRODUCCIÓN A LA CONVECCIÓN
 Ley de enfriamiento de Newton: q   h (Ts  T ) ; q  h A(Ts  T ) .
 Coeficiente de transferencia de calor por convección local, h o promedio, h [W/m2ꞏK].

 Relación entre los coeficientes de convección local y promedio:
1
q   q dAs  (Ts  T )  hdAs  h As (Ts  T )
As As
 h
As 
As
hdAs . Para flujo sobre una
1 L
L 0
placa plana: h  hdx .
 Espesor de la capa límite de velocidad, (x): la y para la que u(y) = 0,99ꞏu.

 Espesor de la capa límite térmica, t(x): la y para la que (Ts - T(y))/(Ts - T) = 0,99.
 k f T / y y 0
 Relación del coeficiente de convección en la capa límite: h  .
Ts  T
u  x u  x
 Número de Reynolds: Re x   .
 
 Número de Reynolds crítico para el inicio de la turbulencia en flujo externo: Rex,c = 5ꞏ105.
 Expresión diferencial de la ecuación de conservación de la masa o de continuidad:
  ( u )  ( v)  kg 
  0  m 3 ꞏs  .
t x y
 Expresiones diferenciales de las ecuaciones de balance de la cantidad de movimiento o del
momento lineal:
 u u u  p    u 2  u v      u v  N kg 
  u v       2            X  3  2 2 
 x y t  x x   x 3  x y   y   y x  m m ꞏs 
 v v v  p    v 2  u v      u v 
  u  v       2            Y .
 x y t  y y   y 3  x y   x   y x 
 Expresión diferencial de la ecuación de conservación de la energía:

  T    T   ( pu )  ( pv)
k    k   q  Xu  Yv     
x  x  y  y  x y
, siendo V2 = u2 + v2.
   V2
    V 2
    V2
  W 
 u  e  gy     v e  gy       e  gy    3 
x   2  y   2  t   2   m 
 Expresión diferencial de la ecuación de conservación de la energía térmica para fluido
incompresible en flujo estacionario:
10
 T T    T    T 
c p  u v   k   k     q .
 x y  x  x  y  y 
 u v  2  u  2  v  2  
 Disipación viscosa:        2        .
 y x   x   y   
 Aproximaciones de capa límite: fluido incompresible ( constante), con propiedades

constantes (k, , etc.), fuerzas de cuerpo insignificantes (X=Y=0) y sin generación de energía (
u u v v T T
q  0 ).Además: u >> v y  , , en la capa límite de velocidad y  en
y x y x y x
la capa límite térmica.
u v
 Ecuación de conservación de la masa o de continuidad en la capa límite:   0.
x y
 Ecuaciones de balance de la cantidad de movimiento o del momento lineal en la capa límite:
u u 1 p  2u p
u v   2 y  0.
x y  x y y
2
T T  2T   u 
 Ecuación de conservación de la energía en la capa límite: u v   2    .
x y y c p  y 

 Número de Prandtl: Pr  .

hL T *
 Número de Nusselt: Nu   .
kf y * y * 0
 Las formas adimensionales de las soluciones de la capa límite adoptan la siguiente forma:
hL
Nu  f  x*, ReL , Pr  y Nu   f Re L , Pr  .
kf
 Relación entre los espesores de las capas límites hidrodinámica y térmica en régimen laminar:

 Pr 1 / 3 .
t
11
TEMA 6. CONVECCIÓN FORZADA EN FLUJO EXTERNO
 Temperatura de película es la temperatura media entre la del fluido y la de la superficie:

Ts  T
Tf  .
2
5 5x
 Espesor de la capa límite laminar:  lam ( x)   .
u  / x Re x
 Correlación de convección local para el flujo laminar sobre una placa plana con temperatura
hx x
superficial constante: Nu x   0,332 Re1x / 2 Pr 1 / 3 . Con propiedades a Tf y Pr  0,6.
k
 Relación entre los espesores de las capas límites de velocidad y térmica en régimen laminar:
 lam
 Pr 1 / 3 .
 t ,lam
 Correlación de convección promedio para el flujo laminar sobre una placa plana con
hx x
temperatura superficial constante: Nu x   0,664Re1x / 2 Pr 1 / 3 . Con propiedades a Tf y Pr 
k
0,6.
 Espesor de la capa límite de velocidad turbulenta:  turb  0,37 xRe x1 / 5 .
 Espesor de la capa límite térmica turbulenta:  t ,turb   turb .
 Correlación de convección local para el flujo turbulento sobre una placa plana con Ts = cte:
Nu x  0,0296 Re x4 / 5 Pr 1 / 3 . Con propiedades a Tf y 0,6 < Pr < 60.
 Para condiciones de capa límite mezclada (laminar y turbulenta) se trabaja con el coeficiente
1  xc
 0 hlam dx  x hturb dx  .
L
de convección promedio: h L 
L c 
 Correlación de convección promedio para capa límite mezclada (laminar y turbulenta) sobre
una placa plana con Ts = cte: Nu L  (0,037 Re L4 / 5  871) Pr 1 / 3 . Con propiedades a Tf y
 0,6  Pr  60 
 8
5ꞏ10  Re L  10  .
5
 Re x  5ꞏ10 5 
 c 
12
 Correlación de convección promedio para capa límite mayoritariamente turbulenta, es decir,

la longitud de la capa límite laminar es despreciable (L >> xc y ReL >> Rex,c), sobre una placa
 0,6  Pr  60 
 
plana con Ts = cte: Nu L  0,037 Re 4/5
L Pr 1/ 3
. Con propiedades a Tf y 5ꞏ10 5  Re L  10 8  .
 Re x  5ꞏ10 5 
 c 
 Correlación de convección local para flujo laminar sobre una placa plana que desprende un
flujo de calor superficial constante: Nu x  0,453 Re 1x / 2 Pr 1 / 3 , con Pr  0,6 y propiedades a Tf.
 Correlación de convección local para flujo turbulento sobre una placa plana que desprende un
flujo de calor superficial constante: Nu x  0,0308 Re x4 / 5 Pr 1 / 3 , con 0,6  Pr  60 y propiedades
a Tf.
 Flujo sobre una placa plana que desprende un flujo de calor superficial constante. La
q s
variación de la temperatura superficial local se obtiene con: Ts ( x)  T  .
hx ( x )
 Número de Reynolds para flujo cruzado sobre un cilindro: ReD = VꞏD/.
 Correlación de Hilpert para flujo cruzado sobre un cilindro: Nu D  CRe Dm Pr 1 / 3 . Los valores
de las constantes C y m se dan en la Tabla 6.3 en función de ReD. La Tabla 6.4 da los valores
de las constantes para cilindros no circulares. Las propiedades se evalúan a Tf. Válida para
fluidos con Pr  0,7.
1/ 4
 Pr 
 Correlación de Zhukauskas para flujo cruzado sobre un cilindro: Nu D  CRe Pr  m
D
n
 .
 Prs 
n  0,37 si Pr  10  0,7  Pr  500

Con   y  6 
. Los valores de las constantes C y m se dan en
n  0,36 si Pr  10  1  ReD  10 
la Tabla 6.5 en función de ReD. Las propiedades se evalúan a T, excepto Prs a Ts.
 Correlación de Churchill y Bernstein para flujo cruzado sobre un cilindro:
4/5
0,62 Re1D/ 2 Pr 1 / 3   ReD  5 / 8 
Nu D  0,3  1     . Con propiedades a Tf y ReDꞏPr > 0,2.
1  (0,4 / Pr ) 
2 / 3 1/ 4
  282.000  
13
 Correlación de Zhukauskas para flujo a través de un banco de tubos:
1/ 4  N L  20 
 Pr  1.000  Re 6
Nu D  CRe Dm,máx Pr 0,36   . Con  D , máx  2ꞏ10  . Las constantes C y m se dan
 Prs   0,7  Pr  500 
en la Tabla 6.6. Las propiedades se evalúan a T  (Tent  Tsal ) / 2 , excepto Prs a Ts. Para NL <
20 se aplica un factor de corrección tal que Nu D  C2 Nu D , donde C2 está dado en

N L  20 N L  20
la Tabla 6.7.
 ReD,máx se define en función de la velocidad máxima del fluido dentro del banco de tubos.
 ST es el espaciado transversal y SL el espaciado longitudinal (distancias entre centros de
tubos).
 Para la configuración alineada la velocidad máxima se da en el plano transversal entre dos
ST
tubos verticales y su valor es Vmáx  V.
ST  D
 Para la configuración escalonada se utiliza la misma expresión si la velocidad máxima se da
ST
en el plano transversal. Pero si se da en el plano diagonal la expresión es Vmáx  V.
2(S D  D)
La velocidad máxima ocurre en el plano diagonal si se cumple la siguiente condición (ver
1/ 2
 2  ST  2  ST  D
Figura 6.2): 2( S D  D)  ( S T  D)  S D   S L      .
  2   2
(Ts  Tent )  (Ts  Tsal )

 Diferencia de temperaturas media logarítmica: Tml  .
 Ts  Tent 
ln 
 Ts  Tsal 
Ts  Tsal  DNh 
 Cálculo de la temperatura de salida del flujo:  exp   , donde N es el
 
Ts  Tent  VN T S T c p 
número total de tubos y NT el número de tubos en el plano transversal.
 Cálculo de la transferencia de calor por unidad de longitud de tubo: q   h NDTml .
14
TEMA 7. CONVECCIÓN FORZADA EN FLUJO INTERNO
 El número de Reynolds para flujo interno se define en función del diámetro del tubo y de la
u m D u m D
velocidad media del fluido sobre la sección transversal del tubo: Re D   . Como
 
4m
m  u m Ac , para un tubo circular el número de Reynolds se puede expresar: ReD  .
D
 Número de Reynolds crítico para el inicio de la turbulencia en flujo interno: ReD,c = 2.300.
 xcd ,h 
 Longitud hidrodinámica de entrada para flujo laminar:    0,05Re D .
 D  lam
 xcd ,h 
 Longitud hidrodinámica de entrada para flujo turbulento: 10     60 .
 D  turb
 Expresión de la velocidad media en función del flujo másico integrado en la sección
m  u ( r , x) dAc 2 ro

Ac
transversal: u m ( x )    u ( r , x ) rdr .
Ac Ac ro2 0
 xcd ,t 
 Longitud de entrada térmica para flujo laminar:    0,05Re D Pr .
 D  lam
 xcd ,t 
 Longitud de entrada térmica para flujo turbulento:    10 .
 D  turb
 Temperatura media definida en función de la energía térmica transportada por el fluido:
U  E t   u (r )c v T ( r , x ) dAc  m c v Tm ( x )  Tm ( x ) 
 uc v TdAc 2 ro

Ac
 u (r )T ( r , x ) rdr .
Ac m c v u m ro2 0
 Bajo condiciones térmicas completamente desarrolladas se cumple:

  Ts ( x )  T ( r , x )  T dTs dTm
   0 . Además:    f ( r ) para q s = cte y
x  Ts ( x)  Tm ( x)  cd ,t x cd ,t dx cd ,t dx cd ,t
T Ts  T dTm
  f ( r ) para Ts = cte.
x cd ,t Ts  Tm dx cd ,t
 Al aplicar un balance de energía al flujo interno en un tubo de un gas ideal o de un líquido

incompresible se obtiene que la transferencia de calor por convección al fluido es igual a la
rapidez a la que aumenta la energía térmica del fluido: q conv  m c p (Tm , sal  Tm ,ent ) .
15
 Variación axial de la temperatura media para el caso de flujo de calor superficial constante:
q sP
Tm ( x)  Tm,ent  x.
m c p
 Variación axial de la temperatura media para el caso de temperatura superficial constante:
Ts  Tm ( x)  Px 
 exp  h.
Ts  Tm ,ent  m c 
 p 
 La transferencia total de calor se expresa en función de la diferencia de temperaturas media
Tsal  Tent (Ts  Tm , sal )  (Ts  Tm ,ent )
logarítmica: q conv  h As Tml ; Tml   .
ln(Tsal / Tent )  (Ts  Tm , sal ) 
ln
 (Ts  Tm ,ent ) 
 Caso general de un tubo con varias resistencias, por ejemplo, un tubo rodeado de un fluido
externo (convección interna y externa simultáneas):
Tsal T  Tm , sal  U As   1  T
  exp    exp   ; qconv  U As Tml  ml .
Tent T  Tm ,ent  m c   m c R  Rtot
 p   p tot 
 Correlación de convección local para flujo laminar en la región completamente desarrollada
hD 48
con flujo de calor superficial constante: Nu    4,36 . Propiedades a Tm.
k 11
 Correlación de convección local para flujo laminar en la región completamente desarrollada
hD
con temperatura superficial constante: Nu   3,66 . Propiedades a Tm.
k
 Número de Graetz: Gz D  ( D / x) Re D Pr .
 Correlación de Hausen para flujo laminar con longitud de entrada térmica (perfil de
velocidades desarrollado) y con temperatura superficial constante:
hD 0,0668 ( D / L ) Re D Pr
Nu D   3,66  . Propiedades a Tm  (Tm ,ent  Tm , sal ) / 2 .
1  0,04( D / L ) Re D Pr 
2/3
k
 Correlación de Sieder y Tate para flujo laminar interno con longitud de entrada combinada y
1/ 3 0 ,14
 Re Pr    
con temperatura superficial constante: Nu D  1,86 D    . Con propiedades a
 L/D   s 
 0,48  Pr  16.700 
Tm  (Tm ,ent  Tm , sal ) / 2 , excepto s a Ts y 0,0044    9,75 .
 
 s 
16
 Correlación de Dittus-Boelter para flujo turbulento interno completamente desarrollado,

válida tanto para flujo de calor como para temperatura superficial constante:
Nu D  0,023 Re D4 / 5 Pr n . Con n = 0,4 para calentamiento (Ts > Tm), n = 0,3 para enfriamiento
0,7  Pr  160
 
(Ts < Tm), las propiedades evaluadas a Tm y  ReD  10.000  .
 ( x / D)  10 
 Correlación de Sieder y Tate para flujo turbulento interno completamente desarrollado y con
grandes variaciones de las propiedades del fluido, válida tanto para flujo de calor como
0 ,14
  
temperatura superficial constante: Nu D  0,027 Re 4/5
Pr 1/ 3
  . Con propiedades a Tm,
 s
D

0,7  Pr  16.700
 
excepto s a Ts y  ReD  10.000  .
 ( x / D)  10 
 El número de Nusselt promedio en flujo turbulento para todo el tubo es igual al valor asociado
con la región completamente desarrollada, Nu D  Nu D ,cd , para valores de (L / D) > 60 y las
propiedades del fluido a Tm  (Tm ,ent  Tm , sal ) / 2 .
4 Ac
 Para tubos no circulares se trabaja con el diámetro hidráulico: Dh  , donde Ac es el
Pmojado
área de la sección transversal y Pmojado el perímetro mojado. Las expresiones del número de
u m Dh u m Dh 4m
Reynolds para el diámetro hidráulico son: Re Dh    .
  P
 Número de Nusselt local para flujo laminar completamente desarrollado en tubos no
circulares: Tabla 7.2.
 Correlaciones de convección para flujo turbulento completamente desarrollado en tubos no
circulares: Las mismas que para tubos circulares trabajando con el diámetro hidráulico.
17
TEMA 8. CONVECCIÓN LIBRE O NATURAL
g (Ts  T ) x 3
 Número de Grashof: Grx  .
2
 Relación entre la convección forzada y la convención libre: si GrL / Re L2  1 , la convección
libre se desprecia frente a la forzada; si GrL / Re L2  1 , la forzada se desprecia frente a la

libre.
 Soluciones de similitud para la convección libre laminar sobre una superficie vertical.
1/ 4
hx  Grx 
Número de Nusselt local: Nu x    ꞏ f ( Pr ) , siendo
k  4 
0,75Pr 1 / 2
f ( Pr )  . Número de Nusselt promedio:
0,609  1,221Pr 1/ 2
 1,238Pr 1/ 4
1/ 4
h L 4  GrL  4
Nu L     ꞏ f ( Pr )  Nu L .
k 3 4  3
g (Ts  T ) x 3
 Número de Rayleigh: Ra x  Grx Pr  .

 Transición entre la capa límite laminar y la turbulenta en placas verticales: Grx,c  109 
Rax,c / Pr  109.
 Correlación de Churchill y Chu para la convección libre sobre una superficie vertical a
temperatura constante aplicable para todo RaL:
2
h L  0,387 Ra1L/ 6 
Nu L   0,825  8 / 27 
. Propiedades calculadas a Tf = (Ts + T)/2.
k  
1  (0,492 / Pr ) 9 / 16  
 Correlación para la convección libre sobre una superficie vertical a temperatura constante
hL 0,670 Ra1L/ 4
aplicable al flujo laminar: Nu L   0,68  con RaL  109.
k 
1  (0,492 / Pr ) 9 / 16 4/9
Propiedades calculadas a Tf = (Ts + T)/2.

 Si la condición de la superficie es un flujo de calor constante en vez de una temperatura
uniforme, la diferencia de temperaturas (Ts - T) aumentará con x. Las correlaciones
anteriores son aplicables en este caso si Nu L y RaL se definen en términos de la diferencia de

temperaturas en el punto medio de la placa, TL / 2  Ts ( L / 2)  T . Como h  q s / TL / 2 es
18
necesario realizar un proceso iterativo para determinar Ts (L/2). Es posible obtener una
expresión para la temperatura en cualquier punto en función de la temperatura en el punto
1/ 5
 x
medio: T x  Ts ( x )  T  1,15  T L / 2 .
 L
 Los resultados anteriores también se pueden aplicar a cilindros verticales de altura L si el
espesor de la capa límite, , es mucho menor que el diámetro del cilindro, condición que viene
D 35
dada por:  .
L GrL1 / 4
 Para placas inclinadas (superficie superior de placa fría o superficie inferior de placa caliente)
se pueden emplear las correlaciones para placas verticales sustituyendo g por gꞏcos () para 0º
   60º ( se mide desde la vertical).
 Para placas horizontales se utiliza una longitud característica definida como el cociente entre
el área y el perímetro de la placa: Lc = As / P.
 Correlaciones de convección libre para la superficie superior de una placa horizontal caliente
o para la superficie inferior de una placa horizontal fría a temperatura constante:
Nu Lc  0,54 Ra1Lc/ 4 si 104  RaLc  107 y Nu Lc  0,15Ra1Lc/ 3 si 107  RaLc  1011. Propiedades
calculadas a Tf.
 Correlaciones de convección libre para la superficie inferior de una placa horizontal caliente o
para la superficie superior de una placa horizontal fría a temperatura constante:
Nu Lc  0,27 Ra1Lc/ 4 con 105  RaLc  1010. Propiedades calculadas a Tf.
 Correlación de Churchill y Chu para la convección libre sobre un cilindro largo horizontal:
2
h D  0,387 Ra1D/ 6 
Nu D   0,60  8 / 27 
con RaD  1012. Propiedades calculadas a Tf.
k  
1  (0,559 / Pr ) 9 / 16  
GrL
 Convección libre y forzada combinadas. Se produce cuando  1. Se utilizan las
ReL2
correlaciones convenientes corregidas con la siguiente expresión:
3
Nu combinada  Nu 3forzada  Nu libre
3
. El signo + se emplea cuando los dos flujos tienen el mismo
sentido o son perpendiculares y el signo - se emplea cuando los dos flujos tienen sentidos
opuestos.
19
TEMA 2. CONDUCCIÓN UNIDIMENSIONAL EN RÉGIMEN ESTACIONARIO.
Tabla 2.1. Eficiencia de formas comunes de aletas.
Descripción Esquema Dimensiones Eficiencia
A f  2 wL c
Aleta recta t tanh mLc
Lc  L  (t / 2)
de perfil f 
rectangular m  ( 2 h / kt ) 1 / 2 mLc
w
siendo w >> t
L
Aleta recta 
A f  2w L2  (t / 2) 2 
1/ 2
1 I1 (2mL)
de perfil m  ( 2 h / kt ) 1/ 2 f 
triangular w
mL I 0 (2mL)
t siendo w >> t
L
20
Tabla 2.1. Eficiencia de formas comunes de aletas (continuación).
Aleta recta

A f  w C1 L2  ( L2 / t ) ln(t / L  C1 ) 
2
de perfil
t
w 
C1  1  (t / L) 2 
1/ 2 f 
4(mL) 2
1 1/ 2
1
parabólico m  ( 2 h / kt ) 1 / 2
L
K1 (mr1 ) I1 (mr2c )  I 1 (mr1 ) K1 (mr2c )

Aleta A f  2 ( r22c  r12 )  f  C2
anular de t I 0 (mr1 ) K1 (mr2c )  K 0 (mr1 ) I 1 (mr2c )
r2 c  r2  (t / 2)
perfil (2r / m)
rectangular m  ( 2 h / kt ) 1 / 2 C2  2 1 2
L (r2c  r1 )
r1
r2
21
Tabla 2.1. Eficiencia de formas comunes de aletas (continuación).

Aleta de A f  DLc
tanh mLc
aguja D Lc  L  (D / 4) f 
cilíndrica mLc
m  ( 4 h / kD ) 1 / 2
L
D
Aleta de
aguja D
Af 
2
L2
 ( D / 2) 2 1/ 2
f 
2 I 2 (2mL)
mL I1 (2mL)
cónica m  ( 4 h / kD ) 1 / 2
L
L3  L 
Af  C3C4  ln(2DC4 / L)  C3 
Aleta de y = (D/2)ꞏ(1 - x/L)2 8D  2D 
2
C 3  1  2( D / L ) 2
f 
aguja D
4 / 9(mL) 2

1
1/ 2
1
parabólica

C 4  1  ( D / L) 
2 1/ 2
L m  ( 4 h / kD ) 1 / 2
22
TEMA 3. CONDUCCIÓN BIDIMENSIONAL EN RÉGIMEN ESTACIONARIO.
Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados

[q = Sk(T1 - T2)].
Descripción
Esquema Restricciones Factor de forma [m]
del sistema
T2
1.1. Esfera
enterrada z
2D
en un T1
z  D/2 S
medio 1  ( D / 4 z)
semiinfinito D
1.2. Esfera
aislado
enterrada
en un
z
medio
2D
T1
semiinfinito T2 T2
S
z  D/2
con 1  (D / 4z)
superficie D T2
aislada
(Bejan pág.
115)
1.3. Esfera T2
enterrada
T1
en un T2 T2
medio Ninguna S  2D
infinito D
(Holman T2
pág. 55)
1.4.
Conducción
entre dos T2
T1 2d
esferas en S
D d w/ D  3 d  ( D / 2w) 4  d
un medio 1  
infinito w
D  1  (d / 2w) 2  w
(Bejan y
Holman)
23
[q = Sk(T1 - T2)] (continuación).
Descripción
del sistema
1.5.
Cavidad aislado D aislado
hemisférica T1
en medio
semiinfinito T2 T2
Ninguna S  D
con T2
superficie
aislada
(Bejan)
T2 2L
2.1. L  D S
Cilindro de z arc cosh(2 z / D)
longitud L T1
enterrado L  D 2L
en un medio L S
z  3D / 2 ln(4 z / D)
semiinfinito D
2.2. 2
zD S 
Cilindro de T2 arc cosh(2 z / D)
longitud z 2
infinita T1 S 
enterrado   z > 2D ln 2 z / D   2 z D 2  1 
 
en un medio
semiinfinito D
(Rohsenow 2
z >> D S 
pág. 3-120) ln4 z / D 
2.3. T2
Cilindro
vertical de L 2L
longitud L T1
L  D S
enterrado ln(4 L / D)
en un medio D
semiinfinito
2.4.
Conducción
entre dos T2
T1 2L
cilindros L  D1 , D2 S
D d  4w 2  D 2  d 2 
paralelos de L  w arc cosh 
longitud L w  2 Dd 
en un medio
infinito
24
Descripción
del sistema
T2
2.5. Cilindro
de longitud L  z 
en medio de
planos
T1
z  D / 2 2L
S
paralelos de L L  z ln(8 z / D)
igual longitud D z
y ancho  
infinito T2
T2
2.6. Cilindro
de longitud L
centrado en 2L
T1 wD S
un sólido de w L  w ln(1,08w / D)
sección
cuadrada de D
igual longitud
T2
T1
2.7. Cilindro
excéntrico de d 2L
longitud L en z Dd S
 D 2  d 2  4z 2 
el interior de L  D arc cosh 
un cilindro de D  2 Dd 
igual longitud
2.8. Fila
infinita de
cilindros de
longitud 2
z>D S 
infinita en un ln(2 L / D)ꞏsenh(2z / L)
medio
semiinfinito
(Rohsenow)
25
Descripción
del sistema
2.9. Fila
infinita de
cilindros de
longitud 2
infinita en el z>D S 
plano medio ln(2 L / D)ꞏsenh(z / L)
de una placa
infinita
(Rohsenow)
3.1. Cubo
enterrado en
L
un medio Ninguna S  8, 24 L
infinito
(Holman)
4.1. 0 , 59
Paralelepípedo   z 
1,685L log1  
inmerso en un
Ninguna S   a 
medio z
0 , 078
semiinfinito  
(Holman) b
4.2. Agujero
de sección a
rectangular  5,7
muy largo en a>b S   2b
 3,5 z 
un medio ln 0, 25 0,75 
semiinfinito a b 
(Rohsenow)
26

Factor de forma
Descripción del sistema Esquema Restricciones
[m]
W
5.1. Pared plana con H

A
H  L/5 WH
superficies isotermas T1 S
T2 W  L/5 L
(Bejan)
T2
5.2. Esquina de dos paredes L

S  0,54W
W  L/5
contiguas
T1
T2
T2
L
5.3. Esquina de tres
T1
paredes contiguas con T2
diferencia de temperaturas Ninguna S  0,15 L
entre las superficies W
interior y exterior (Bejan) T2
T2
T2
D T1
6.1. Disco delgado sobre

Ninguna S  2D
medio semiinfinito T2
T1
6.2. Disco delgado
horizontal enterrado en un z 4,45D
T2 Ninguna S
medio semiinfinito (Kreith 1  ( D / 5,67 z )
pág. 112) D
27
Descripción
del sistema
6.3. Disco T1
delgado
horizontal z 2D
enterrado en T2 zD S
un medio ( / 2)  arc tan(D / 4 z )
D
semiinfinito
(Bejan)
6.4. Disco
delgado
horizontal aislado
enterrado en
un medio
T1 z T1
2D
T2 zD S
semiinfinito ( / 2)  arc tan(D / 4 z )
con D
superficie
aislada
(Bejan)
6.5. Dos
discos T1 T2
paralelos 2D
coaxiales en D L/ D  2 S
un medio ( / 2)  arc tan(D / 2L)
infinito L
(Bejan)
7.1. Placa 2L
horizontal z  L S
T1 L W ln(4 L / W )
delgada de
anchura W
z0 L
(dimensión z S
 al dibujo) T2 L W ln(4 L / W )
enterrada en
un medio L 2L
semiinfinito L  W S
(Bejan y z  2W ln(2z / W )
Holman)
28
Descripción del Factor de forma

Esquema Restricciones
sistema [m]
7.2. Placa vertical

delgada y larga según
0 , 24
la dimensión  al 1 z
  12 L
S   2,38ꞏ 
dibujo enterrada en 2 L z
un medio semiinfinito
(Rohsenow)
T1
7.3. Placa horizontal
delgada y larga según z 0 , 32
la dimensión  al 1 z
  12 L
T2 S   2,94ꞏ 
dibujo enterrada en 2 L z
un medio semiinfinito
L
(Rohsenow)
29
TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO.
Tabla 4.1. Coeficientes de la aproximación con un término de las soluciones de conducción

transitoria unidimensional.
Pared plana Cilindro infinito Esfera

Bi 1 (rad) C1 1 (rad) C1 1 (rad) C1
0,01 0,0998 1,0017 0,1412 1,0025 0,1730 1,0030
0,02 0,1410 1,0033 0,1995 1,0050 0,2445 1,0060
0,03 0,1732 1,0049 0,2439 1,0075 0,2989 1,0090
0,04 0,1987 1,0066 0,2814 1,0099 0,3450 1,0120
0,05 0,2217 1,0082 0,3142 1,0124 0,3852 1,0149
0,06 0,2425 1,0098 0,3438 1,0148 0,4217 1,0179
0,07 0,2615 1,0114 0,3708 1,0173 0,4550 1,0209
0,08 0,2791 1,0130 0,3960 1,0197 0,4860 1,0239
0,09 0,2956 1,0145 0,4195 1,0222 0,5150 1,0268
0,10 0,3111 1,0160 0,4417 1,0246 0,5423 1,0298
0,15 0,3779 1,0237 0,5376 1,0365 0,6608 1,0445
0,20 0,4328 1,0311 0,6170 1,0483 0,7593 1,0592
0,25 0,4801 1,0382 0,6856 1,0598 0,8448 1,0737
0,30 0,5218 1,0450 0,7465 1,0712 0,9208 1,0880
0,40 0,5932 1,0580 0,8516 1,0932 1,0528 1,1064
0,50 0,6533 1,0701 0,9408 1,1143 1,1656 1,1441
0,60 0,7051 1,0814 1,0185 1,1346 1,2644 1,1713
0,70 0,7506 1,0919 1,0873 1,1539 1,3225 1,1978
0,80 0,7910 1,1016 1,1490 1,1725 1,4320 1,2236
0,90 0,8274 1,1107 1,2048 1,1902 1,5044 1,2488
1,0 0,8603 1,1191 1,2558 1,2071 1,5708 1,2732
2,0 1,0769 1,1795 1,5995 1,3384 2,0288 1,4793
3,0 1,1925 1,2102 1,7887 1,4191 2,2889 1,6227
4,0 1,2646 1,2287 1,9081 1,4698 2,4556 1,7201
5,0 1,3138 1,2402 1,9898 1,5029 2,5704 1,7870
6,0 1,3496 1,2479 2,0490 1,5253 2,6537 1,8338
7,0 1,3766 1,2532 2,0937 1,5411 2,7165 1,8674
8,0 1,3978 1,2570 2,1286 1,5526 2,7654 1,8921
9,0 1,4149 1,2598 2,1566 1,5611 2,8044 1,9106
10,0 1,4289 1,2620 2,1795 1,5677 2,8363 1,9249
20,0 1,4961 1,2699 2,2881 1,5919 2,9857 1,9781
30,0 1,5202 1,2717 2,3261 1,5973 3,0372 1,9898
40,0 1,5325 1,2723 2,3455 1,5993 3,0632 1,9942
50,0 1,5400 1,2727 2,3572 1,6002 3,0788 1,9962
100,0 1,5552 1,2731 2,3809 1,6015 3,1102 1,9990
 1,5707 1,2733 2,4050 1,6018 3,1415 2,0000
Bi = hL/k para la pared plana y Bi = hro/k para el cilindro infinito y la esfera.
30
Tabla 4.2. Soluciones de conducción transitoria de sistemas multidimensionales expresadas

como producto de soluciones de sistemas unidimensionales.
Sistema Esquema Solución
S (x, t)
Sólido T ( x, t )  T
S ( x, t ) 
semiinfinito x Tini  T sólido
semiinfini to
P (x, t)
T ( x, t )  T
Pared plana P ( x, t ) 
x Tini  T pared
plana
2L1


C (r, t)
Cilindro T ( r , t )  T
C (r , t ) 
infinito r Tini  T cilindro
infinito
ro


S (x1, t)ꞏP (x2, t)
Placa
S ( x1 , t )ꞏP( x 2 , t )
semiinfinita x2
x1
2L2
31
Tabla 4.2. Soluciones de conducción transitoria de sistemas multidimensionales expresadas

como producto de soluciones de sistemas unidimensionales (continuación).
Sistema Esquema Solución

P (x1, t)ꞏP (x2, t)

Barra
rectangular x1
x2 P( x1 , t )ꞏP( x 2 , t )
infinita
 2L1
2L2
C (r, t)ꞏS (x, t)

Cilindro C ( r , t )ꞏS ( x , t )
semiinfinito r
x
ro
Barra P (x1, t)ꞏP (x2, t)ꞏS (x3, t)

rectangular S ( x 3 , t )ꞏP ( x1 , t )ꞏP ( x 2 , t )
semiinfinita x3
x2
x1
2L1
2L2
P (x1, t)ꞏP (x2, t)ꞏP (x3, t) x

3 2L3
Paralelepípedo P ( x1 , t )ꞏP ( x 2 , t )ꞏP ( x 3 , t )
rectangular x1
x2
2L1
2L2
C (r, t)ꞏP (x, t)
x
Cilindro corto 2L1 C ( r , t )ꞏP ( x , t )
r
ro
32
TEMA 6. CONVECCIÓN FORZADA EN FLUJO EXTERNO
Tabla 6.1. Tabla resumen de correlaciones para flujo externo sobre placa plana.
Correlaciones Transferencia de calor Condiciones

h x q x  hx Ts  T 
1* Nu x  x  0,332 Re 1x / 2 Pr 1 / 3 Placa a temperatura Ts constante. Régimen laminar. Valor local en x. Pr > 0,6.
k dq  hx Ts  T dAs
hx x q   h Ts  T  Placa a Ts constante. Régimen laminar. Valor promedio entre 0 y x (ó entre 0 y x = L).
2** Nu x   0,664 Re1x / 2 Pr 1 / 3
k q  h As Ts  T  Pr > 0,6.
hx x q x  hx Ts  T 
3* Nu x   0,0296Rex4 / 5 Pr 1 / 3 Placa a Ts constante. Régimen turbulento. Valor local en x. 0,6 < Pr < 60.
k dq  hx Ts  T dAs
hL L q   h Ts  T  Placa a Ts constante. Régimen mixto (parte laminar y parte turbulento). Valor
4** Nu L   (0,037 Re L4 / 5  871) Pr 1 / 3
k q  h As Ts  T  promedio entre 0 y x = L. Pr > 0,6. 5ꞏ105 < ReL < 108.
q   h Ts  T  Placa a Ts constante. Régimen predominantemente turbulento (parte laminar
hL L
5** Nu L   0,037 Re L4 / 5 Pr 1 / 3 despreciable  L >> xc y ReL >> Rex,c). Valor promedio entre 0 y x = L. Pr > 0,6.
k q  h As Ts  T  5ꞏ105 < ReL < 108.
q s Placa que desprende un flujo de calor uniforme. Régimen laminar. Valor local en x. Pr
6* Nu x  0, 453 Re 1x / 2 Pr 1 / 3 q s  cte. / Ts ( x)  T 
h x ( x) > 0,6.
q s
Placa que desprende un flujo de calor uniforme. Rég. turbulento. Valor local en x. 0,6
7* Nu x  0,0308 Re x4 / 5 Pr 1 / 3 q s  cte. / Ts ( x)  T 
< Pr < 60. h x ( x)
u x u x u L u L    c p c p
*: Rex     **: ReL     Condición de rég. turbulento para placa plana: Rex,c > 5ꞏ105 Número de Prandtl: Pr   
      k k
En todas las correlaciones las propiedades del fluido se calculan a la temperatura de película: T f  T s  T   / 2 ; Ts: Temperatura de la superficie [K]; T: Temp. del
flujo libre [K]; As: Área de transferencia de calor [m2]; ν: viscosidad cinemática [m2/s]; μ: viscosidad dinámica [N/m2ꞏs]; α: difusividad térmica [m2/s]; k:
conductividad térmica del fluido [W/mꞏK].
33
Tabla 6.2. Tabla resumen de correlaciones para flujo cruzado sobre cilindros.
Correlaciones para flujo cruzado sobre un cilindro Transferencia de calor Condiciones
Correlación de Hilpert. Los valores de las constantes C y m se dan en la
hD q   h Ts  T  Tabla 6.3 en función de ReD. La Tabla 6.4 da los valores de las constantes
1 Nu D   CReDm Pr 1 / 3
k q  h As Ts  T  para cilindros no circulares. Las propiedades se evalúan a Tf. Válida para
fluidos con Pr  0,7.
 n  0,37 si Pr  10 
Correlación de Zhukauskas. Con  n  0,36 si Pr  10  y
 
q  h Ts  T 
1/ 4
hD  Pr 
2 Nu D   CRe Dm Pr n  
  0,7  Pr  500
k  Prs  q  h As Ts  T   1  Re  10 6  . Los valores de las constantes C y m se dan en la Tabla
 D 
6.5 en función de ReD. Las propiedades se evalúan a T, excepto Prs a Ts
q  h Ts  T 
4/5
hD   Re D  5 / 8 
0,62 Re 1D/ 2 Pr 1 / 3 Correlación de Churchill y Bernstein. Con propiedades a Tf y ReDꞏPr >
3 Nu D  0,3  1    
k  
1/ 4
  282000  
1  (0,4 / Pr ) 2 / 3 q  h As Ts  T  0,2.
Correlaciones para flujo cruzado sobre un banco de N

Transferencia de calor Condiciones
cilindros
 N L  20 
1.000  Re  2ꞏ10 6
Correlación de Zhukauskas. Con  D , máx . Las
1/ 4
q  h Tml  0,7  Pr  500 
hD  Pr 
4 Nu D   CRe Dm, máx Pr 0 , 36  
 q  h DTml constantes C y m se dan en la Tabla 6.6. Las propiedades se evalúan a
k  Prs  T  (Tent  Tsal ) / 2 , excepto Prs a Ts. Para NL < 20 se aplica un factor de
q  q NL  h NDLTml
corrección tal que Nu D  C 2 Nu D , donde C2 está dado en la
N L  20 N L  20
Tabla 6.7
u  D u  D Vmax D Vmax D Ts  Tsal  DNh 
ReD   ReD , max   Cálculo de la temperatura de salida del flujo:  exp  
    Ts  Tent  VN S c 
 T T p
(Ts  Tent )  (Ts  Tsal ) ST: espaciado transversal; SL: espaciado longitudinal; NT: número de tubos en direc. transversal;
Diferencia de temperaturas media logarítmica: Tml 
 T  Tent  NL: número de tubos en direc. longitudinal; N = NT x NL: núm. total de tubos.
ln  s 
 Ts  Tsal 
ST ST
Config. alineada: Vmáx  V ; Config. escalonada: Vmáx  V si
ST  D ST  D
ST
2( S D  D)  (S T  D) ó Vmáx  V si 2(S D  D)  (S T  D) ;
2( S D  D )
34
Tabla 6.3. Coeficientes de la correlación de Hilpert para flujo cruzado sobre un cilindro
(Pr  0,7).
ReD C m
0,4 - 4 0,989 0,330
4 - 40 0,911 0,385
40 - 4.000 0,683 0,466
4.000 - 40.000 0,193 0,618
40.000 - 400.000 0,027 0,805
Tabla 6.4. Coeficientes de la correlación de Hilpert para flujo de aire cruzado sobre un
paralelepípedo.
Geometría Dibujo ReD C m
Cuadrado en V
D 5ꞏ103 - 105 0,246 0,588
diagonal
V
Cuadrado recto D 5ꞏ103 - 105 0,102 0,675
V 5ꞏ103 - 1,95ꞏ104 0,160 0,638

Hexágono recto D
1,95ꞏ104 - 105 0,0385 0,782
Hexágono en V
D 5ꞏ103 - 105 0,153 0,638
diagonal
Placa vertical D 4ꞏ103 - 1,5ꞏ104 0,228 0,731
Tabla 6.5. Coeficientes de la correlación de Zhukauskas para flujo de aire cruzado sobre
un cilindro.
ReD C m
1 - 40 0,75 0,4
40 - 1.000 0,51 0,5
103 - 2ꞏ105 0,26 0,6
2ꞏ105 - 106 0,076 0,7
35
Figura 6.1. Nu local para flujo de aire cruzado sobre un cilindro. (Incropera)
Tabla 6.6. Coeficientes de la correlación de Zhukauskas para el flujo cruzado sobre un

banco de tubos.
Configuración ReD,máx C m
2
Alineado 10 - 10 0,80 0,40
Escalonado 10 - 102 0,90 0,40
2 3
Alineado 10 - 10 Se aproxima como un cilindro único
Escalonado 102 - 103 Se aproxima como un cilindro único
3 5
Alineado (ST / SL > 0,7) 10 - 2ꞏ10 0,27 0,63
3 5 1/5
Escalonado (ST / SL < 2) 10 - 2ꞏ10 0,35(ST / SL) 0,60
Escalonado (ST / SL > 2) 103 - 2ꞏ105 0,40 0,60
Alineado 2ꞏ105 - 2ꞏ106 0,021 0,84
Escalonado 2ꞏ105 - 2ꞏ106 0,022 0,84
Para ST / SL < 0,7, la transferencia de calor es ineficiente y los tubos alineados no se deben usar.
36
Tabla 6.7. Coeficiente de corrección C2 de la correlación de Zhukauskas para el flujo

cruzado sobre un banco de tubos para NL < 20 y ReD > 103.
NL 1 2 3 4 5 7 10 13 16
Alineado 0,70 0,80 0,86 0,90 0,92 0,95 0,97 0,98 0,99
Escalonado 0,64 0,76 0,84 0,89 0,92 0,95 0,97 0,98 0,99
Figura 6.2. Disposición de los tubos en configuración alineada (a) y escalonada (b) en un
banco de tubos. (Incropera)
NL
NT
Figura 6.3. Esquema de un banco de tubos en flujo cruzado. (Incropera)
37
TEMA 7. CONVECCIÓN FORZADA EN FLUJO INTERNO
Figura 7.1. Número de Nusselt local en la región de entrada para flujo laminar en el
interior de un tubo circular con temperatura superficial uniforme. (Bejan)
Figura 7.2. Número de Nusselt local en la región de entrada para flujo laminar en el
interior de un tubo circular con flujo de calor superficial uniforme. (Bejan)
38
Diagrama 7.1. Metodología para seleccionar las correlaciones de convención forzada en flujo interno.
- Región de entrada  (x < xcd,t ó Gz-1 < 0,05)  Figuras 7.1. y 7.2. Siendo: xcd,t = 0.05ꞏDꞏReDꞏPr.
- Correlación local:
- q s  cte  NuD = 4,36.

- Región c. d.  NuD = cte - Ts = cte  NuD = 3,66.
(x > xcd,t ó Gz-1 > 0,05): - Tubo no circular  Tabla 7.1.
ReD < 2.300  Régimen Laminar:

- Problema de longitud de entrada térmica (si xcd,t >> xcd,h; Pr >> 1):
Correlación de Hausen.
- Reg. de entrada + c. d.:
- Problema de longitud de entrada combinada (si O(xcd,t)O(xcd,h); O(Pr)1):
Correlación de Sieder y Tate.
- Correlación promedio:
- q s  cte  Nu D  4,36 .
- Región c. d.  Nu D  cte : - Ts = cte  Nu D  3,66 .
- Tubo no circular  Tabla 7.2.
- Correlación de Dittus-Boelter.
- Correlación local
(en región c. d.: x/D > 10): - Correlación de Sieder y Tate:
(se usa con aceite siempre; con agua y con aire si hay grandes T).
ReD > 2.300  Régimen Turbulento:

- Correl. de Dittus-Boelter.
Tm ,ent  Tm , sal
- Correlación promedio con prop. a Tm  .
2
(condiciones c. d.: L/D > 60): - Correl. de Sieder y Tate.
* En el temario de este curso no se estudia la región de entrada en régimen turbulento.

39
Tabla 7.1. Tabla resumen de las correlaciones de convención forzada en flujo interno.
Correlaciones para tubos circulares Transferencia de calor Condiciones
hD q x  hTs  Tm  Tubo sometido a un flujo de calor superficial uniforme, q x  cte . Régimen laminar, correlación
1 Nu   4,36
k dq  hTs  Tm Ddx local, región completamente desarrollada. Propiedades calculadas a Tm.
hD q x  hTs  Tm  Tubo sometido a una temperatura superficial uniforme, Ts  cte . Régimen laminar, correlación
2 Nu   3,66
k dq  hTs  Tm Ddx local, región completamente desarrollada. Propiedades calculadas a Tm.
Correlación de Hausen. Tubo sometido a Ts  cte . Régimen laminar, correlación promedio,
hD 0,0668 ( D / L ) Re D Pr
3 Nu D   3,66  q  h PLTml región de entrada + c. d., problema de longitud de entrada térmica (perfil de velocidades
1  0,04( D / L ) Re D Pr 
2/3
k
desarrollado, xcd,t >> xcd,h, Pr >> 1). Propiedades calculadas a Tm  (Tm,ent  Tm,sal ) / 2 .
Correlación de Sieder y Tate. Tubo sometido a Ts  cte . Régimen laminar, correlación
1/ 3 0 ,14
hD  Re Pr     promedio, región de entrada + c. d., problema de longitud de entrada combinada (O(xcd,t) 
4 Nu D   1,86 D    q  h PLTml
k   O(xcd,h)). Propiedades calculadas a Tm  (Tm,ent  Tm,sal ) / 2 , excepto s a Ts. Rango de validez:
 L/D   s 
0,48 < Pr < 16.700 y 0,0044 < ( / s) < 9,75.
Correlación de Dittus-Boelter. Tubo sometido a qx  cte. o Ts  cte. Régimen turbulento,
correlación local, región completamente desarrollada. Con n = 0,4 para calentamiento (Ts > Tm)
hD q x  hTs  Tm 
5 Nu D   0,023Re D4 / 5 Pr n 0,7  Pr  160
k dq  hTs  Tm Ddx  
y n = 0,3 para enfriamiento (Ts < Tm). Propiedades a Tm. Rango de validez:  ReD  10.000  .
 ( x / D)  10 
Correlación de Sieder y Tate. Tubo sometido a qx  cte. o Ts  cte. Régimen turbulento,
correlación local, región completamente desarrollada. Grandes variaciones de las propiedades
q x  hTs  Tm 
0 ,14
hD   
6 Nu D   0,027 Re D4 / 5 Pr 1 / 3   0,7  Pr  16.700
k  s

 dq  hTs  Tm Ddx  
del fluido. Popiedades calculadas a Tm, excepto s a Ts. Rango de validez:  ReD  10.000  .
 ( x / D)  10 
hD q  h PLTml Mismas condiciones que correlación 5, pero correlación promedio para flujo completamente
7 Nu D   0,023ReD4 / 5 Pr n desarrollado, (L / D) > 60. Propiedades calculadas a Tm  (Tm,ent  Tm,sal ) / 2 .
k q  q / L q   q / PL
0 ,14
hD    q  h PLTml Mismas condiciones que correlación 6, pero correlación promedio para flujo completamente
8 Nu D   0,027 Re D4 / 5 Pr 1 / 3  

k  s  q  q / L q   q / PL desarrollado, (L / D) > 60. Propiedades calculadas a Tm  (Tm,ent  Tm,sal ) / 2 , excepto s a Ts.
u m D u m D (Ts  Tm ,ent )  (Ts  Tm , sal )
ReD   Las correlaciones 1 y 2 son válidas como promedio si L >> xcd,t Diferencia de temperaturas media logarítmica: Tml 
   Ts  Tm ,ent 
ln 
T T 
 s m , sal 
Tubos de sección no circular con régimen laminar, correlación local y región c. d.: Tabla 7.1.
Tubos de sección no circular con régimen turbulento y región c. d.: Correlaciones 5, 6, 7 u 8, pero trabajando con el diámetro hidráulico, Dh = 4ꞏAc / P. Ac: área de la sección
transversal. P: perímetro.
40
Tabla 7.2. Números de Nusselt para flujo laminar completamente desarrollado en tubos de
diferente sección transversal.
hDh
NuD 
k
b
Sección transversal q s uniforme Ts uniforme
a
Circular - 4,36 3,66
Rectangular
1,0 3,61 2,98
(a = altura, b =base)
Rectangular
1,43 3,73 3,08
Rectangular
2,0 4,12 3,39
Rectangular
3,0 4,79 3,96
Rectangular
4,0 5,33 4,44
Rectangular
8,0 6,49 5,60
Rectangular (a =
 8,23 7,54
altura, b =base)
Triangular - 3,11 2,47
41
TEMA 8. CONVECCIÓN LIBRE O NATURAL
Figura 8.1. Perfiles de velocidad y de temperatura para la capa límite laminar de

convección libre sobre una superficie vertical isoterma. (Incropera)
Figura 8.2. Vista lateral de los patrones de flujo de la convección libre sobre placas planas
inclinadas: Ts > T a la izquierda y Ts < T a la derecha.
42
Tabla 8.1. Tabla resumen de correlaciones de convención libre.
Transferencia de Representación
Correlación Condiciones
calor gráfica
1/ 4 Placa vertical con temperatura superficial constante, Ts = cte.
Nu x 
hx  Grx 
  ꞏ
0,75Pr 1 / 2 q x  hTs  Tm  Régimen laminar, Grx < 109. Correlación local.
1 k  4  0,609  1,221Pr 1 / 2  1,238Pr 1 / 4
dq  hTs  Tm dAz Ts > T Ts < T
Correlación promedio: Nu L 
hL 4
 Nu L
k 3
q  h Ts  T 
2
hL  0,387 Ra1L/ 6  Correlación de Churchill y Chu. Placa vertical con Ts = cte.
2 Nu L   0,825  9 / 16 8 / 27 
k  1  (0,492 / Pr )   q  h As Ts  T  Correlación promedio. Válida para todo RaL.
hL 0,670 Ra 1L/ 4 q  h Ts  T 

3 Nu L   0,68  Placa vertical con Ts = cte. RaL  109. Correlación promedio.
k 1  (0,492 / Pr ) 9 / 16 4 / 9 q  h As Ts  T 
h Lc Ts > T Ts < T
Nu Lc   0,54 Ra1Lc/ 4 Placa horizontal con Ts = cte. Superficie superior
4 k q  h Ts  T  de placa caliente o inferior de placa fría.
con 104  RaLc  107
h Lc q  h As Ts  T  Correlación promedio. Longitud característica
Nu Lc   0,15 Ra 1Lc/ 3 definida como el cociente entre el área y el
5 k perímetro de la placa: Lc = As / P.
con 107  RaLc  1011
Ts > T Ts < T
Placa horizontal con Ts = cte. Superficie inferior de
h Lc q  h Ts  T  placa caliente o superior de placa fría. Correlación
Nu Lc   0,27 Ra 1Lc/ 4
6 k q  h As Ts  T  promedio. Longitud característica definida como el
con 105  RaLc  1010 cociente entre el área y el perímetro de la placa: Lc
= As / P.
q  h Ts  T 
2
hD  0,387 Ra1D/ 6  Correlación de Churchill y Chu (promedio) para la
Nu D   0,60 
7 k  1  (0,559 / Pr )  
9 / 16 8 / 27
q  h As Ts  T 
Ts > T convección libre sobre un cilindro largo horizontal:
con RaD  1012.
g (Ts  T ) x 3 Ra x g (Ts  T ) L3 g (Ts  T ) D 3

Grx   ; RaL  GrL Pr  ; Ra D  ; Propiedades calculadas a Tf = (Ts + T)/2; Correlaciones 1 a 3: válidas para q x  cte si
2 Pr  
propiedades, Nu L y RaL se definen en función de la temperatura en el punto medio de la placa: Tf = (Ts(L/2) + T)/2; TL / 2  Ts ( L / 2)  T  h  qs / TL / 2 
Tx  Ts ( x)  T  1,15x / L TL / 2 ; Correlaciones 1 a 3: válidas para cilindros verticales de altura L si el espesor de la capa límite, , es mucho menor que el diámetro del cilindro
1/ 5
 D / L   35 / GrL1 / 4  ; Para placas inclinadas (superficie superior de placa fría o superficie inferior de placa caliente) se pueden emplear las correlaciones 1 a 3 sustituyendo g por
gꞏcos () para 0º    60º ( se mide desde la vertical).
43
Tabla 8.2. Metodología de resolución de problemas de convección según el tipo de

condición de contorno.
Flujo externo, placa plana:

Ts (local) q 𝑞" (local)
constante  ℎ 𝐴 𝑇 𝑇  ℎ 𝑇 𝑇
𝑞′′  𝑞" 𝐴  constante

𝑇 ℎ
Flujo externo, cilindro:

constante  ℎ 𝐴 𝑇 𝑇  ℎ 𝑇 𝑇
q′′  𝑞" 𝐴  constante

T h
Flujo externo, banco de tubos:

ℎ 𝐴 ∆𝑇 
constante  ó h T T
𝑚 𝑐 𝑇 , 𝑇 , 
T h
𝑇 es local (cambia en la dirección del flujo)
Flujo interno:
ℎ 𝐴 ∆𝑇 
constante  ó h T T
𝑚 𝑐 𝑇 , 𝑇 , 
T h
𝑇 es local (cambia en la dirección del flujo)
44
TABLAS DE PROPIEDADES TERMOFÍSICAS Y DE FUNCIONES MATEMÁTICAS
Tabla A. Propiedades termofísicas del aire a presión atmosférica.
T  cp ꞏ107 ꞏ106 kꞏ103 ꞏ106 Pr

(K) (kg/m ) (J/kgꞏK) (Nꞏs/m ) (m /s) (W/mꞏK) (m /s)
3 2 2 2
100 3,5562 1032 71,1 2,00 9,34 2,54 0,786

150 2,3364 1012 103,4 4,426 13,8 5,84 0,758
200 1,7548 1007 132,5 7,590 18,1 10,3 0,737
250 1,3947 1006 159,6 11,44 22,3 15,9 0,720
300 1,1614 1007 184,6 15,89 26,3 22,5 0,707
350 0,9950 1009 208,2 20,92 30,0 29,9 0,700
400 0,8711 1014 230,1 26,41 33,8 38,3 0,690
450 0,7740 1021 250,7 32,39 37,3 47,2 0,686
500 0,6964 1030 270,1 38,79 40,7 56,7 0,684
550 0,6329 1040 288,4 45,57 43,9 66,7 0,683
600 0,5804 1051 305,8 52,69 46,9 76,9 0,685
650 0,5356 1063 322,5 60,21 49,7 87,3 0,690
700 0,4975 1075 338,8 68,10 52,4 98,0 0,695
750 0,4643 1087 354,6 76,37 54,9 109 0,702
800 0,4354 1099 369,8 84,93 57,3 120 0,709
850 0,4097 1110 384,3 93,80 59,6 131 0,716
900 0,3868 1121 398,1 102,9 62,0 143 0,720
950 0,3666 1131 411,3 112,2 64,3 155 0,723
1000 0,3482 1141 424,4 121,9 66,7 168 0,726
1100 0,3166 1159 449,0 141,8 71,5 195 0,728
1200 0,2902 1175 473,0 162,9 76,3 224 0,728
1300 0,2679 1189 496,0 185,1 82 238 0,719
1400 0,2488 1207 530 213 91 303 0,703
1500 0,2322 1230 557 240 100 350 0,685
45
Tabla B. Propiedades termofísicas del aceite de motor a presión atmosférica.
T  cp ꞏ102 ꞏ106 kꞏ103 ꞏ107 Pr ꞏ103

(K) (kg/m ) (J/kgꞏK) (Nꞏs/m ) (m /s) (W/mꞏK) (m /s)
3 2 2 2
(K-1)
273 899,1 1796 385 4280 147 0,910 47000 0,70
280 895,3 1827 217 2430 144 0,880 27500 0,70
290 890,0 1868 99,9 1120 145 0,872 12900 0,70
300 884,1 1909 48,6 550 145 0,859 6400 0,70
310 877,9 1951 25,3 288 145 0,847 3400 0,70
320 871,8 1993 14,1 161 143 0,823 1965 0,70
330 865,8 2035 8,36 96,6 141 0,800 1205 0,70
340 859,9 2076 5,31 61,7 139 0,779 793 0,70
350 853,9 2118 3,56 41,7 138 0,763 546 0,70
360 847,8 2161 2,52 29,7 138 0,753 395 0,70
370 841,8 2206 1,86 22,0 137 0,738 300 0,70
380 836,0 2250 1,41 16,9 136 0,723 233 0,70
390 830,6 2294 1,10 13,3 135 0,709 187 0,70
400 825,1 2337 0,874 10,6 134 0,695 152 0,70
410 818,9 2381 0,698 8,52 133 0,682 125 0,70
420 812,1 2427 0,564 6,94 133 0,675 103 0,70
430 806,5 2471 0,470 5,83 132 0,662 88 0,70
46
Tabla C. Propiedades termofísicas del agua saturada.
T P  ifg cp ꞏ106 kꞏ103 Pr ꞏ106

(K) (bar) (kg/m ) (kJ/kg) (J/kgꞏK) (Nꞏs/m ) (W/mꞏK)
3 2
(K-1)
273,15 0,00611 1000 2502 4217 1750 569 12,99 -68,05
275 0,00697 1000 2497 4211 1652 574 12,22 -32,74
280 0,00990 1000 2485 4198 1422 582 10,26 46,04
285 0,01387 1000 2473 4189 1225 590 8,81 114,1
290 0,01917 999,0 2461 4184 1080 598 7,56 174,0
295 0,02617 998,0 2449 4181 959 606 6,62 227,5
300 0,03531 997,0 2438 4179 855 613 5,83 276,1
305 0,04712 995,0 2426 4178 769 620 5,20 320,6
310 0,06221 993,0 2414 4178 695 628 4,62 361,9
315 0,08132 991,1 2402 4179 631 634 4,16 400,4
320 0,1053 989,1 2390 4180 577 640 3,77 436,7
325 0,1351 987,2 2378 4182 528 645 3,42 471,2
330 0,1719 984,3 2366 4184 489 650 3,15 504,0
335 0,2167 982,3 2354 4186 453 656 2,88 535,5
340 0,2713 979,4 2342 4188 420 660 2,66 566,0
345 0,3372 976, 2329 4191 389 668 2,45 595,4
350 0,4163 973,7 2317 4195 365 668 2,29 624,2
355 0,5100 970,9 2304 4199 343 671 2,14 652,3
360 0,6209 967,1 2291 4203 324 674 2,02 697,9
365 0,7514 963,4 2278 4209 306 677 1,91 707,1
370 0,9040 960,6 2265 4214 289 679 1,80 728,7
373,15 1,0133 957,9 2257 4217 279 680 1,76 750,1
375 1,0815 956,9 2252 4220 274 681 1,70 761
380 1,2869 953,3 2239 4226 260 683 1,61 788
385 1,5233 949,7 2225 4232 248 685 1,53 814
390 1,794 945,2 2212 4239 237 686 1,47 841
400 2,455 937,2 2183 4256 217 688 1,34 896
410 3,302 928,5 2153 4278 200 688 1,24 852
420 4,370 919,1 2123 4302 185 688 1,16 1010
430 5,699 909,9 2091 4331 173 685 1,09
440 7,333 900,9 2059 4360 162 682 1,04
450 9,319 890,5 2024 4400 152 678 0,99
460 11,71 879,5 1989 4440 143 673 0,95
470 14,55 868,1 1951 4480 136 667 0,92
480 17,90 856,9 1912 4530 129 660 0,89
490 21,83 844,6 1870 4590 124 651 0,87
500 26,40 831,3 1825 4660 118 642 0,86
ifg: entalpía específica del cambio de fase entre líquido y gas.
47
Tabla D. Funciones hiperbólicas.

x senh x cosh x tanh x x senh x cosh x tanh x
0,00 0,0000 1,0000 0,00000 2,00 3,6269 3,7622 0,96403
0,10 0,1002 1,0050 0,09967 2,10 4,0219 4,1443 0,97045
0,20 0,2013 1,0201 0,19738 2,20 4,4571 4,5679 0,97574
0,30 0,3045 1,0453 0,29131 2,30 4,9370 5,0372 0,98010
0,40 0,4108 1,0811 0,37995 2,40 5,4662 5,5569 0,98367
0,50 0,5211 1,1276 0,46212 2,50 6,0502 6,1323 0,98661
0,60 0,6367 1,1855 0,53705 2,60 6,6947 6,7690 0,98903
0,70 0,7586 1,2552 0,60437 2,70 7,4063 7,4735 0,99101
0,80 0,8881 1,3374 0,66404 2,80 8,1919 8,2527 0,99263
0,90 1,0265 1,4331 0,71630 2,90 9,0596 9,1146 0,99396
1,00 1,1752 1,5431 0,76159 3,00 10,018 10,068 0,99505
1,10 1,3356 1,6685 0,80050 3,50 16,543 16,573 0,99818
1,20 1,5095 1,8107 0,83365 4,00 27,290 27,308 0,99933
1,30 1,6984 1,9709 0,86172 4,50 45,003 45,014 0,99975
1,40 1,9043 2,1509 0,88535 5,00 74,203 74,210 0,99991
1,50 2,1293 2,3524 0,90515 6,00 201,71 201,72 0,99999
1,60 2,3756 2,5775 0,92167 7,00 548,32 548,32 1,0000
1,70 2,6456 2,8283 0,93541 8,00 1490,5 1490,5 1,0000
1,80 2,9422 3,1075 0,94681 9,00 4051,5 4051,5 1,0000
1,90 3,2682 3,4177 0,95624 10,00 11013 11013 1,0000
Tabla E. Función gaussiana de error.
x erf (x) x erf (x) x erf (x)
0,00 0,00000 0,36 0,38933 1,04 0,85865
0,02 0,02256 0,38 0,40901 1,08 0,87333
0,04 0,04511 0,40 0,42839 1,12 0,88679
0,06 0,06762 0,44 0,46623 1,16 0,89910
0,08 0,09008 0,48 0,50275 1,20 0,91031
0,10 0,11246 0,52 0,53790 1,30 0,93401
0,12 0,13476 0,56 0,57162 1,40 0,95229
0,14 0,15695 0,60 0,60386 1,50 0,96611
0,16 0,17901 0,64 0,63459 1,60 0,97635
0,18 0,20094 0,68 0,66378 1,70 0,98379
0,20 0,22270 0,72 0,69143 1,80 0,98909
0,22 0,24430 0,76 0,71754 1,90 0,99279
0,24 0,26570 0,80 0,74210 2,00 0,99532
0,26 0,28690 0,84 0,76514 2,20 0,99814
0,28 0,30788 0,88 0,78669 2,40 0,99931
0,30 0,32863 0,92 0,80677 2,60 0,99976
0,32 0,34913 0,96 0,82542 2,80 0,99992
0,34 0,36936 1,00 0,84270 3,00 0,99998
2 w  u2
erf w 

 e du
0
erfc w  1  erf w
48
Tabla F. Primeras cuatro raíces de la ecuación trascendental, nꞏtan(n) = Bi, para

conducción transitoria en una pared plana.
hL
Bi  1 2 3 4
k
0 0,0000 3,1416 6,2832 9,4248
0,001 0,0316 3,1419 6,2833 9,4249
0,002 0,0447 3,1422 6,2835 9,4250
0,004 0,0632 3,1429 6,2838 9,4252
0,006 0,0774 3,1435 6,2841 9,4254
0,008 0,0893 3,1441 6,2845 9,4256
0,01 0,0998 3,1448 6,2848 9,4258
0,02 0,1410 3,1479 6,2864 9,4269
0,04 0,1987 3,1543 6,2895 9,4290
0,06 0,2425 3,1606 6,2927 9,4311
0,08 0,2791 3,1668 6,2959 9,4333
0,1 0,3111 3,1731 6,2991 9,4354
0,2 0,4328 3,2039 6,3148 9,4459
0,3 0,5218 3,2341 6,3305 9,4565
0,4 0,5932 3,2636 6,3461 9,4670
0,5 0,6533 3,2923 6,3616 9,4775
0,6 0,7051 3,3204 6,3770 9,4879
0,7 0,7506 3,3477 6,3923 9,4983
0,8 0,7910 3,3744 6,4074 9,5087
0,9 0,8274 3,4003 6,4224 9,5190
1,0 0,8603 3,4256 6,4373 9,5293
1,5 0,9882 3,5422 6,5097 9,5801
2,0 1,0769 3,6436 6,5783 9,6296
3,0 1,1925 3,8088 6,7040 9,7240
4,0 1,2646 3,9352 6,8140 9,8119
5,0 1,3138 4,0336 6,9096 9,8928
6,0 1,3496 4,1116 6,9924 9,9667
7,0 1,3766 4,1746 7,0640 10,0339
8,0 1,3978 4,2264 7,1263 10,0949
9,0 1,4149 4,2694 7,1806 10,1502
10,0 1,4289 4,3058 7,2281 10,2003
15,0 1,4729 4,4255 7,3959 10,3898
20,0 1,4961 4,4915 7,4954 10,5117
30,0 1,5202 4,5615 7,6057 10,6543
40,0 1,5325 4,5979 7,6647 10,7334
50,0 1,5400 4,6202 7,7012 10,7832
60,0 1,5451 4,6353 7,7259 10,8172
80,0 1,5514 4,6543 7,7573 10,8606
100,0 1,5552 4,6658 7,7764 10,8871
 1,5708 4,7124 7,8540 10,9956
49
Tabla G. Funciones de Bessel de primera clase.
x J0(x) J1(x)
0,0 1,0000 0,0000
0,1 0,9975 0,0499
0,2 0,9900 0,0995
0,3 0,9776 0,1483
0,4 0,9604 0,1960
0,5 0,9385 0,2423
0,6 0,9120 0,2867
0,7 0,8812 0,3290
0,8 0,8463 0,3688
0,9 0,8075 0,4059
1,0 0,7652 0,4401
1,1 0,7196 0,4709
1,2 0,6711 0,4983
1,3 0,6201 0,5220
1,4 0,5669 0,5419
1,5 0,5118 0,5579
1,6 0,4554 0,5699
1,7 0,3980 0,5778
1,8 0,3400 0,5815
1,9 0,2818 0,5812
2,0 0,2239 0,5767
2,1 0,1666 0,5683
2,2 0,1104 0,5560
2,3 0,0555 0,5399
2,4 0,0025 0,5202
50
Tabla H. Funciones de Bessel modificadas de primera y segunda clase.

x e-xꞏI0(x) e-xꞏI1(x) exꞏK0(x) exꞏK1(x)
0,0 1,0000 0,0000  
0,2 0,8269 0,0823 2,1408 5,8334
0,4 0,6974 0,1368 1,6627 3,2587
0,6 0,5993 0,1722 1,4167 2,3739
0,8 0,5241 0,1945 1,2582 1,9179
1,0 0,4658 0,2079 1,1445 1,6362
1,2 0,4198 0,2153 1,0575 1,4429
1,4 0,3831 0,2185 0,9881 1,3011
1,6 0,3533 0,2190 0,9309 1,1919
1,8 0,3289 0,2177 0,8828 1,1048
2,0 0,3085 0,2153 0,8416 1,0335
2,2 0,2913 0,2121 0,8057 0,9738
2,4 0,2766 0,2085 0,7740 0,9229
2,6 0,2639 0,2047 0,7459 0,8790
2,8 0,2528 0,2007 0,7206 0,8405
3,0 0,2430 0,1968 0,6978 0,8066
3,2 0,2343 0,1930 0,6770 0,7763
3,4 0,2264 0,1892 0,6580 0,7491
3,6 0,2193 0,1856 0,6405 0,7245
3,8 0,2129 0,1821 0,6243 0,7021
4,0 0,2070 0,1788 0,6093 0,6816
4,2 0,2016 0,1755 0,5953 0,6627
4,4 0,1966 0,1725 0,5823 0,6454
4,6 0,1919 0,1695 0,5701 0,6292
4,8 0,1876 0,1667 0,5586 0,6143
5,0 0,1835 0,1640 0,5478 0,6003
5,2 0,1797 0,1614 0,5376 0,5872
5,4 0,1762 0,1589 0,5280 0,5749
5,6 0,1728 0,1565 0,5188 0,5634
5,8 0,1697 0,1542 0,5101 0,5525
6,0 0,1667 0,1521 0,5019 0,5422
6,4 0,1611 0,1479 0,4865 0,5232
6,8 0,1561 0,1441 0,4724 0,5060
7,2 0,1515 0,1405 0,4595 0,4905
7,6 0,1473 0,1372 0,4476 0,4762
8,0 0,1434 0,1341 0,4366 0,4631
8,4 0,1399 0,1312 0,4264 0,4511
8,8 0,1365 0,1285 0,4168 0,4399
9,2 0,1334 0,1260 0,4079 0,4295
9,4 0,1305 0,1235 0,3995 0,4198
9,6 0,1278 0,1213 0,3916 0,4108
10,0 1,0000 0,0000
I n 1 ( x )  I n 1 ( x )  ( 2 n / x ) I n ( x )
51
ALFABETO GRIEGO
Mayúsculas Minúsculas Nombre

  alfa
  beta
  gamma
  delta
  épsilon
  seta o zeta
  eta
  zeta o theta
  iota
  kappa o cappa
  lambda
  my o mu
  ny o nu
  xi
  ómicron
  pi
  ro o rho
 ,  sigma
  tau
  ípsilon
  fi o phi
  ji
  psi
  omega
52

FTF 2021

Cargado por

Copyright:

Formatos disponibles

FTF 2021

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

FTF 2021

Cargado por

Copyright:

Formatos disponibles

FÓRMULAS, TABLAS Y FIGURAS DE

Juan Carlos Ramos González

Raúl Antón Remírez

Fórmulas, Tablas y Figuras

Fórmulas, Tablas y Figuras

Tema 1. Introducción a la transferencia de calor y a la conducción ...................................... 1

Tema 2. Conducción unidimensional en régimen estacionario ........................................... 20

Tablas de propiedades termofísicas y de funciones matemáticas ........................................ 45

Alfabeto griego .................................................................................................................... 52

Fórmulas, Tablas y Figuras

Fórmulas, Tablas y Figuras

TEMA 1. INTRODUCCIÓN A LA TRANSFERENCIA DE CALOR Y A LA

 Calor o transferencia de calor o velocidad de transferencia de calor: q [J/s = W].

 Coeficiente de transferencia de calor por convección: h [W/m2ꞏK].

 Constante de Stefan-Boltzmann:  = 5,67ꞏ10-8 W/m2ꞏK4.

 Intercambio de radiación para una superficie gris y difusa ( = ):

 Principio de conservación de la energía en un volumen de control formulado en un instante de

Fórmulas, Tablas y Figuras

 Principio de conservación de la energía en un volumen de control formulado en un intervalo

 Principio de conservación de la energía en una superficie de control: E ent  E sal  0 .

 Capacidad térmica volumétrica:  cp [J/m3ꞏK]. Mide la capacidad de un material para

 Ecuación de difusión de calor en coordenadas cilíndricas (r radial,  angular o longitud, z

 Ecuación de difusión de calor en coordenadas esféricas (r radial,  polar, cenital o colatitud, 

 Condición de contorno de primera clase o de Dirichlet: superficie mantenida a temperatura

Fórmulas, Tablas y Figuras

 Condición de contorno de segunda clase o de Neumann: flujo de calor fijo o constante en la

 Condición de contorno de tercera clase o de Robin: corresponde a la transferencia de calor por

está en contacto con la superficie de la pared opuesta al origen de coordenadas:

TEMA 2. CONDUCCIÓN UNIDIMENSIONAL EN RÉGIMEN ESTACIONARIO

Fórmulas, Tablas y Figuras

 Resistencia térmica de conducción para una pared esférica:

área en función de la cual se exprese: U 1 A1  U 2 A2  U 3 A3  ...  U i Ai   Rt  .

 Distribución de temperaturas y transferencia de calor para aletas de área de sección

 ( x) cosh m( L  x)  (h / mk )senh m( L  x) senh mL  (h / mk ) cosh mL

 Caso C, extremo con temperatura establecida ((x = L) = L):

 Caso D, aleta muy larga (L   y L  0, aplicable si mꞏL > 2,65):

Fórmulas, Tablas y Figuras

sección transversal de la base de la aleta. El uso de aletas sólo se justifica cuando f  2.

 Teniendo en cuenta la resistencia térmica de convección de la base de la aleta, Rt,b = 1/hAc,b,

Fórmulas, Tablas y Figuras

TEMA 3. CONDUCCIÓN BIDIMENSIONAL EN RÉGIMEN ESTACIONARIO

 Factor de forma de conducción para sistemas bidimensionales, S: q  Sk (T1  T2 ) . Se obtiene

TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO

Rt,int U ꞏLc Tint

 Longitud característica: Lc = Vol/As. Para una pared plana de espesor 2L sometida a

Fórmulas, Tablas y Figuras

 Q(t )  U (0  t )  U (t )  U (0) .

Fórmulas, Tablas y Figuras

Fórmulas, Tablas y Figuras

 En un sólido semiinfinito la condición de frontera interior es T(x, t) = Tini y la condición

 MDF en régimen transitorio. Expresión en diferencias finitas de la velocidad de variación de

Fórmulas, Tablas y Figuras

TEMA 5. INTRODUCCIÓN A LA CONVECCIÓN

 Ley de enfriamiento de Newton: q   h (Ts  T ) ; q  h A(Ts  T ) .

 Coeficiente de transferencia de calor por convección local, h o promedio, h [W/m2ꞏK].

 Espesor de la capa límite de velocidad, (x): la y para la que u(y) = 0,99ꞏu.

 Expresión diferencial de la ecuación de conservación de la energía:

Fórmulas, Tablas y Figuras

 Aproximaciones de capa límite: fluido incompresible ( constante), con propiedades