Formulario Transferencia de Calor

INSTITUTO TECNOLÓGICO SUPERIOR
DE LIBRES
FÓRMULAS, TABLAS Y FIGURAS DE

TRANSFERENCIA DE CALOR
Alumno: Enrique Rodríguez Hernández
Docente: Ing. Guillermo Córdova

Fórmulas, Tablas y Figuras Transferencia de Calor
2
Enrique Rodríguez Hernández 5to B
Transferencia de Calor
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  k . q x  qx  A . En condiciones de régimen estacionario y con una
x dx
dT T2  T1 TT T
distribución lineal de temperaturas: qx  k  k k 1 2 k .
dx L L L
 Conductividad térmica: k [W/m·K].
 Ley de enfriamiento de Newton: qx  h(Ts  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 T 4s.
 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  Ts 4 , 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: Gabs  G . Irradiación de los
alrededores: G  Talr4 .
 Intercambio de radiación para una superficie gris y difusa ( = ):

q   E (T )  G   T 4  T 4  .
q
También se puede expresar como:
rad b s s alr
A
q rad  hrad (Ts  Talr ) , siendo hrad el coeficiente de transferencia de calor por radiación:
h   (T  T )T 2  T 2 .
rad s alr s alr
 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
 Principio de conservación de la energía en un volumen de control formulado en un intervalo
de tiempo (t): Eent  Egen  Esal  Ealm .
1
Transferencia de Calor Enrique Rodríguez Hernández 5to B
 Principio de conservación de la energía en una superficie de control: E& ent  E& sal  0 .
r  r T r T r T  r r r
 Ley de Fourier vectorial: q  kT    j   i q  jq   .
k i kq z
 x k 

y
y z
x
 Capacidad térmica volumétrica:  cp [J/m ·K]. Mide la capacidad de un material para

3
almacenar energía térmica.

k
 Difusividad térmica:   [m2/s]. Mide la capacidad de un material para conducir energía
 cp
térmica en relación con su capacidad para almacenarla.
 Ecuación de difusión de calor en coordenadas cartesianas:
  k T      kT   &q  cp T  W .
    kT 
x x y y z z t m3 
     
T
 Ecuación de difusión de calor vectorial: ·(kT )  q&  c .
t
p
 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  k dT dx ), esta ecuación

 x
 
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   kr T   1   k T     kT   q  cp T.
r r r r 2   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   kr 2 T   2 1 2   T  1   T  T
r 2 r r r sen    k    r 2 sen2  ksen   q&  cp .
  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 Fourier: corresponde a la transferencia de calor

por convección en la superficie, q cond ,superficie  qconv . Si el fluido está en contacto con la
 T
superficie de la pared donde está el origen de coordenadas: k  hT  T (x  0, t). Si
x x0
el fluido 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 hr A
T 1
 Coeficiente global de transferencia de calor, U: qx  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
la transferencia de calor, qx, es una constante independiente de x): q x     k (T )dT .
T
x0 A(x) T0
 Resistencia térmica de conducción para una pared cilíndrica: R (Ts1  Ts 2 ) ln(r2 / r1 )

  .
t ,cond
qr 2Lk
1 1
 Resistencia térmica de convección para una pared cilíndrica: R   .
Ah 2rLh
t ,conv
 Resistencia térmica de conducción para una pared esférica:
Rt ,cond  (Ts1  Ts 2 )  1  1  1 
q 4  .
r k  r1 r2 
3
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
 U3A 3  ...  Ui Ai   Rt  .
1
área en función de la cual se exprese: U 1A 1 U A2 2
 Generación de energía térmica por unidad de volumen: q&  e& E& gen  W 
 .
gen

Vol m  3 
 1 cdA  dT 
 Ecuación de calor para una aleta: d T2     1 sh(T  T
2
dA
  ) 0.
dx A dx dx A k dx  
 c   c 
 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
( hA T (L)  T   kA dT ):
c  c
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
hP
siendo  (x)  T (x)  T , b  Tb  T , m2  , M  hPkAc b , P el perímetro y Ac
kAc

el área transversal.
d
 Caso B, extremo adiabático (   0 ):
dx x L
 (x) cosh m(L  x) q  M tanh mL


b
f
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 q f  M  hPkAc b
 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:   , siendo A el área de la
hA 
f
c,b
c,b b
sección transversal de la base de la aleta. El uso de aletas sólo se justifica cuando f  2.
4

 Resistencia térmica de una aleta: Rt , f  q .
b
 Teniendo en cuenta la resistencia térmica de convección de la base de la aleta, Rt,b = 1/hAc,b,
se puede expresar la efectividad como:   Rt ,b .

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.
qmáx hAf 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  .
hAf  f
 Para el caso de una aleta recta de sección transversal uniforme y con su extremo adiabático se
tiene:   M tanh mL  tanh mL .

hPL b
f
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
aletas, Af, más la parte expuesta de la base, Ab. Es decir, At  NAf  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:
NAf
qt  Nq f  qb  N f hAf b  hAbb  qt  o hAtb   o  1  (1   f ) .
At
5
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 (2D)  .
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  k (y·1) ) o la ley
x
(m1,n)( m,n)
de enfriamiento de Newton ( q ()(m,n)  hx·1(T  Tm,n ) .
TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO
hLc Rt ,cond
 Número de Biot: Bi   .
k Rt ,conv
 Longitud característica: Lc = Vol/As. Para una pared plana de espesor 2L  Lc = L, y para un

cilindro largo o una esfera de radio ro  Lc = ro.
·t
 Número de Fourier: Fo  .
L2c
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 
 ini  Tini  T  exp t  exp    exp(Bi·Fo) .
 Volc p    t 
 

1
 Constante de tiempo térmica:  t  · Volcp  Rt ·Ct , siendo Rt la resistencia a la
hAs
transferencia de calor por convección y Ct la capacidad térmica del sólido.
6
 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á:
Q(t)  qdt  hA t
 (t)dt  Q(t)  Volc 
1 exp t    E .

 
t
0 s 0 p ini t alm
 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
 * (x*,t*)   C exp( 2 Fo) cos( x*)   * cos( x*) , siendo
Tini  T

1 1 1 o 1
T (x  0, t)  T
 * (x*  0, t)   C exp( 2 Fo) la temperatura del plano medio (x* = x/L =
o
Tini  T 1 1
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
2L sometida a convección: Q(t)  1  sen 1  * (0, t) , siendo Uo = Uini = cpVol(Tini - T) =

 
1
o
Uo
Umá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:
T (r, t)  T
 * (r*,t*)   C exp( 2 Fo)J ( r*)   * J ( r*) , siendo
Tini  T 1 1 0 1 o 0 1

 * (r*  0, t)  To  T  C exp( 2 Fo) la temperatura del eje central (r* = r/ro = 0) y J0
Tini  T
o 1 1
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 o * (0, t)  ) , siendo J1 la función de Bessel de
ro sometido a convección:  1 J 1 ( 1
1
 
Uo
primera clase de orden uno cuyos valores se encuentran en la Tabla G y Uo = Uini =
cpVol(Tini - T) = Umá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.
7
distribución de temperaturas en una esfera de radio ro sometida a convección:
T (r, t)  T 1 1
 * (r*,t*)   C exp( 2 Fo) sen( r*)   * sen( r*) , siendo
Tini  T 1 1
 
1 1 1 o 1
r * r *
 * (r*  0, t)  To  T  C exp( 2 Fo) la temperatura del eje central (r* = r/ro = 0).
Tini  T
o 1 1
 Transferencia total de energía en función del tiempo desde o hacia una esfera de radio ro
3 *
 1  o sen( )   cos( ), siendo Uo = Uini = cpVol(Tini -
Q(t)
sometida a convección:
13
 
1 1 1
Uo
T) = Umá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.
 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
semiinfinito plana infinito
 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 
T  T  erf 2 t 
Temperatura superficial constante: T ini s  
(0, t) = Ts T k (Ts  Tini )
qs (t)  k 
x x0 (t)1/ 2
Condición de frontera Distribución de temperaturas

Flujo de calor superficial 2q (t /  )1 / 2  x 2  q x  x ini
T (x, t)  T   exp   o erfc  
constante: qs  q o
o
 4t  k  2 t 
k    
T (x, t)  Tini  x 
 erfc
Convección superficial:   
T T  Tini  2 t 
k  hT T (0, t)
x x0   hx h2t   x h t 

  
exp erfc  
 k k 2    2 t k 
8
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.
 MDF en régimen transitorio. Expresión en diferencias finitas de la velocidad de variación de
T p1  Tm,n
p
la energía almacenada: E& alm  cpV m,n .

t
 Criterio de estabilidad en el MDF explícito: el coeficiente asociado con el nodo de interés en
p
el tiempo anterior (coeficiente de Tm,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
TEMA 5. INTRODUCCIÓN A LA CONVECCIÓN
 Ley de enfriamiento de Newton: q  h(Ts  T ) ; q  hA(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dA  (T  T ) hdA  h A (T  T )  h hdA . Para flujo sobre una

As s s  As s s s 
As 
As s
1

L
placa plana: h  hdx .
L 0
 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 

 Número de Reynolds: Re  u x  u x .
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)
0
kg 
  .
t x y m3 ·s 
 Expresiones diferenciales de las ecuaciones de balance de la cantidad de movimiento o del
momento lineal:
9
 u  v u  u  p   2 u 2 u v
u 
   u  v
  X  N
kg 
x y t x x 3 x y x y y x m  m 2 ·s 2 
3
         
 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.
   V 2     V 2     V 2  W 

 u e  gy    v e  gy     e  gy  2   m3 
x   2  y   2  t   
 Expresión diferencial de la ecuación de conservación de la energía térmica para fluido
incompresible en flujo estacionario:
 T 
 T T     k T    k    q& .
c  u  v 
  
 
 x y  x  x  y  y 
p
 u v 2  u 2  v  2  2  u v  2 
 Disipación viscosa:        2       3  x  y  .
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
T T
( q&  0 ).Además: u >> v y u  u , v , v 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   2 u p
u v y 0.
x y  x y 2 y
 2T   u 
2
 Ecuación de conservación de la energía en la capa límite: u T T

x  v y   y 2  c  y  .
 p  
 
 Número de Prandtl: Pr  .

 T *
hL
 Número de Nusselt: Nu  .
kf y * y*0
10
 Las formas adimensionales de las soluciones de la capa límite adoptan la siguiente forma:
hL
Nu  f x*, ReL , Pr  y Nu   f ReL , Pr  .
kf

 Relación entre los espesores de las capas límites hidrodinámica y térmica:   Pr1 / 3 .
t
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
T f .
2
5 5x
 Espesor de la capa límite laminar: lam (x)   .
u /x Rex
 Correlación de convección local para el flujo laminar sobre una placa plana con temperatura
hx x
superficial constante: Nu   0,332Re1 / 2 Pr1 / 3 . Con propiedades a Tf y Pr  0,6.
x x
k
 lam
 Relación entre los espesores de las capas límites de velocidad y térmica:  Pr1 / 3 .
 t ,lam
 Correlación de convección promedio para el flujo laminar sobre una placa plana con
h x
temperatura superficial constante: Nu x  x  0,664Re1 / 2 Pr1 / 3 . Con propiedades a Tf y Pr 

x
k
0,6.
 Correlación de convección local para flujo laminar de un fluido con número de Prandtl
pequeño, como los metales líquidos, en los que t   y se puede suponer velocidad
uniforme a través de la capa límite térmica: Nu x 0,565(Re ·Pr)
x
1/2
, con propiedades a Tf, Pr
 0,05 y Rex·Pr  100.


 Espesor de la capa límite de velocidad turbulenta: turb  0,37xRe1/
x
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,0296Re 4 / 5 Pr 1/ 3 . Con propiedades a Tf y 0,6 < Pr < 60.
x
 Para condiciones de capa límite1 mezclada (laminar y turbulenta) se trabaja con el coeficiente
  xc h dx  h dx .
de convección promedio: h
 
L
L
L 0 lam x c turb 
11
 Correlación de convección promedio para capa límite mezclada (laminar y turbulenta) sobre
una placa plana con Ts = cte: Nu L  (0,037Re4L/ 5  871)Pr1/ 3 . Con propiedades a Tf y
 0,6  Pr  60 
 5 8 
5·10 
 Re  5·10
Re  10 .

L
 xc 5
 Correlación de convección local para flujo laminar sobre una placa plana con longitud inicial
no calentada (la capa límite de velocidad comienza en x = 0 y la térmica en x = ):
Nux  0
Nux  , donde Nux  0 está dado por la expresión anterior para flujo laminar.
1  ( / x) 
3/ 4 1/ 3
 Correlación de convección local para flujo turbulento sobre una placa plana con longitud
inicial no calentada (la capa límite de velocidad comienza en x = 0 y la térmica en x = ):
Nux  0
Nux  , donde Nux  0 está dado por la expresión anterior para flujo
1  ( / x) 
9 / 10 1 / 9
turbulento.
 Correlación de convección promedio entre  y L sobre una placa plana con longitud inicial no
calentada: Nu L  Nu L
 0
L
L  

1  ( / L) (2 p 1) /( 2 p 2) 
2 p /( 2 p 1)
, siendo p = 1 para flujo laminar
y p = 4 para turbulento.
 Correlación de convección local para flujo laminar sobre una placa plana que desprende un
flujo de calor superficial constante: Nux  0,453Re1/x 2 Pr1/ 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: Nux  0,0308Re4x / 5 Pr1 / 3 , con 0,6  Pr  60 y propiedades
a Tf.
 Flujo sobre una placa plana que desprende un flujo de calor superficial constante. La
variación de la temperatura superficial local se obtiene con: T (x)  T  qs

s 
.
h x (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  CRemDPr1 / 3 . Los valores
de las constantes C y m se dan en la Tabla 6.1 en función de ReD. La Tabla 6.2 da los valores
de las constantes para cilindros no circulares. Las propiedades se evalúan a Tf.
12
Transferencia de Calor Fórmulas, Tablas y Figuras
1/4
 Pr 
 Correlación de Zhukauskas para flujo cruzado sobre un cilindro: Nu D  CReD Pr  Pr  .
m n
 s 
  0,7  Pr  500
Con n  0,37 si Pr  10
n  0,36 si Pr  10  1  Re  106 
y . Los valores de las constantes C y m se dan en
   D 
la Tabla 6.3 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:

D 
Nu D  0,3 D0,62Re 1 Pr1  
1/ 2 /3 4/5

Re
 
5/8
 . Con propiedades a Tf y ReDPr > 0,2.
1  (0,4 / Pr) 
2 / 3 1/ 4
  282.000  
 Correlación de Zhukauskas
para flujo a través de un banco de tubos:

1/ 4
 Pr   N L  20 

  6
Nu D  CRe D,máx Pr  Pr  . Con 1.000  Re D,máx  2·10 . Las constantes C y m se dan

m 0,36
 
 s     
 0,7 Pr 500
en la Tabla 6.4. 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

NL 20 NL 20
la Tabla 6.5.
 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 V má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 V máx  V.
2(S D  D)
La velocidad máxima ocurre en el plano diagonal si se cumple la siguiente condición (ver
  S 2 1 / 2 S  D
 D)  (S  D)  S  S   T  
2
Figura 6.2): 2(S  T .
 2  
D T D L
 2
 Diferencia de temperaturas media logarítmica: T ml  (Ts  Tent )  (Ts  Tsal ) .

 T  T 
ent 
ln s 
 s
T  T sal 
13
 Ts  Tsal  exp DNh 
 , donde N es el
 Cálculo de la temperatura de salida del flujo:
 

Ts  Tent  VNT ST 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 .
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
um D um D
velocidad media del fluido sobre la sección transversal del tubo: Re   . Como
D
 
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: Re D,c = 2.300.
 xcd ,h 
 Longitud hidrodinámica de entrada para flujo laminar:    0,05ReD .
 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
transversal: u m (x)  
Ac
Ac Ac ro2 0

u(r, x)rdr .
 xcd ,t 
 Longitud de entrada térmica para flujo laminar:    0,05ReD 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&  u(r)c T (r, x)dA  mc T (x)  T

ro
  ucvTdAc  2 u(r)T (r, x)rdr .


(x)
Ac
A
c
&vm m
c
t v m& cv u mro 2 0
 Bajo condiciones térmicas completamente desarrolladas se cumple:

 T   dTm
 s (x)  T (r, x)   0 . Además: T  dTs  f (r) para qs = cte y

x T (x)  T (x) 
 s m  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
14
 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 ) .
 Variación axial de la temperatura media para el caso de flujo de calor superficial constante:
 qsP
Tm (x)  Tm,ent m& c x .
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
 h A T ; T Tsal  Tent (Ts  Tm,sal )  (Ts  Tm,ent )

logarítmica: q   .
conv s ml ml
ln(Tsal / Tent ) (T s  Tm,sal ) 
ln  
 (Ts
 T )
m,ent 
 Caso de un tubo rodeado de un fluido externo (convección interna y externa simultáneas):

Tsal T  Tm,sal  UAs   1  Tml

 
 exp   
 exp & 
 ; qconv  UA s Tml  Rtot .
Tent T  Tm,ent  m& c p 

 mc p Rtot 
 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: GzD  (D / x)ReD Pr .
 Correlación de Hausen para flujo laminar con longitud de entrada térmica (perfil de
velocidades desarrollado) y con temperatura superficial constante:
hD

Nu D   3,66  0,0668(D / L)ReD Pr . Propiedades a T  (T T )/2.
1  0,04(D / L)ReD Pr 2 / 3

m m,ent m,sal
k
 Correlación de Sieder y Tate para flujo laminar interno con longitud de entrada combinada y
Pr    0,14
1/3
Re
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 
15
 Correlación de Dittus-Boelter para flujo turbulento interno completamente desarrollado:
NuD  0,023Re 4 / 5 Pr n. Con n = 0,4 para calentamiento (Ts > Tm), n = 0,3 para enfriamiento
D
0,7  Pr  160
(Ts < Tm), las propiedades evaluadas a Tm y Re D  10.000  .

 
 (x / D)  10 
 Correlación de Sieder y Tate para flujo turbulento interno completamente desarrollado y con
0,14 
 
grandes variaciones de las propiedades del fluido: Nu  0,027Re 4 / 5 Pr  1/  . Con
  s 

D 3
D
0,7  Pr  16.700
propiedades a Tm, excepto s a Ts y  Re  10.000 .
 D 
 (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  NuD,cd , para valores de (L / D) > 60 y las
propiedades del fluido a Tm  (Tm,ent  Tm,sal ) / 2 . 
 Para tubos no circulares se trabaja con el diámetro hidráulico: D  4 Ac , donde Ac es el

h
Pmojado
área de la sección transversal y Pmojado el perímetro mojado.

 Número de Nusselt local para flujo laminar completamente desarrollado en tubos no
circulares: Tabla 7.1.
 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.
TEMA 8. CONVECCIÓN LIBRE O NATURAL
g (T  T )x3
 Número de Grashof: Grx s  .
2
 Relación entre la convección forzada y la convención libre: si GrL / ReL2  1, la convección
libre se desprecia frente a la forzada; si GrL / ReL2  1, la forzada se desprecia frente a la

libre.
 Soluciones de similitud para la convección libre laminar sobre una superficie vertical.
hx Gr 1 / 4
Número de Nusselt local: Nu    x  f (Pr) ,

siendo
k  4 
x
16
1/2
0,75Pr
f (Pr)  . Número de Nusselt promedio:
0,609  1,221Pr1 / 2  1,238Pr1 / 4
 
4  Gr 1/ 4
 L  f (Pr)  4 Nu .
hL
Nu L   
3  4 
L
k 3
g (T  T )x 3
 Número de Rayleigh: Ra  Gr Pr s  .

x x

 Número de Rayleigh crítico para la transición de flujo laminar a turbulento en placas
verticales: Rax,c  109.
 Correlación de Churchill y Chu para la convección libre sobre una superficie vertical a
temperatura constante aplicable para todo RaL:

Nu L  h L   L 0,387Ra1 /6 2

0,825    . Propiedades calculadas a Tf = (Ts + T)/2.



1  (0,492 / Pr) 

k  9 / 16 8 / 27

 Correlación para la convección libre sobre una superficie vertical a temperatura constante
1/4
hL 0,670Ra

aplicable al flujo laminar: Nu L   0,68  L
con RaL/Pr  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
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
medio: Tx  Ts (x)  T  1,15 x  TL / 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 (superficies inferior y superior de placas frías y calientes,

respectivamente) se pueden emplear las correlaciones para placas verticales sustituyendo g
por g·cos () para 0º    60º ( se mide desde la vertical).
17
 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 L  0,54Ra1/L 4 si 10  RaL  10 y Nu L  0,15Ra1/L 3 si 10  RaL  10 . Propiedades

4 7 7 11
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 L  0,27Ra1L/ 4 con 105  RaL  1010. Propiedades calculadas a Tf.
 Correlación de Churchill y Chu para la convección libre sobre un cilindro largo horizontal:

Nu D  h D   D 0,387Ra1 /6 2

0,60    con RaD  1012. Propiedades calculadas a Tf.

8 / 27 

1  (0,559 / Pr) 

k  9 / 16


GrL
 Convección libre y forzada combinadas. Se produce cuando  1. Se utilizan las
Re2L
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.
TEMA 9. INTRODUCCIÓN A LA RADIACIÓN
 Ángulo sólido diferencial: d = dAn / r2 = sen·d·d. Unidad: estereorradián (sr). El ángulo

sólido subtendido por el hemisferio sobre un diferencial de área dA1 vale:
2  / 2
 d   
h 0 0
sendd  2 sr .
 Intensidad espectral emitida: I (T , ,  , ) 

dqe  W 
 ,e dA cosdd m2 ·sr·µm .
 
 E (T , )  q (T , ) 
1
 dq
e 2  / 2
 (T , , , ) cosd 
I
dA1d 0 0  W 
  ,e   ,e
 Potencia emisiva espectral: .
(T , ,  ,  ) cossendd  2 
 0 0 I  ,e
2  / 2 
 m ·µm 

 


 Potencia emisiva total: E(T )  q r ad ,emit   E (T , )d [W/m2].


0
18
 Emisor difuso: I  ,e (T , ,  ,  )  I  ,e (T ,  )  E (T , )  I  ,e (T , )  E  Ie .
 Irradiación espectral: G ()  q ()  dqi 2  / 2 (, , ) cossendd
I
dA1·d  0 

  ,i  ,i
0
[W/m2·m].

 Irradiación total: G  q r ad ,inc   G ()d [W/m2].
0
 Radiación incidente difusa: I  ,i (,  ,  )  I  ,i ()  G ()  I ,i ()  G  Ii .

 Radiosidad espectral: J ()  q ()  dqer 2  / 2 (, , ) cossendd
I
dA1·d  0 0

 rad ,emit reflj  ,er
[W/m2·m].

 Radiosidad total: J  q r ad ,emit ref   J  ()d [W/m2].
0
 Emisor y reflector difuso: I  ,er (, ,  )  I  ,er ()  J  ()  I ,er ()  J  Ier .
 2hc 2 -34
Distribución de Planck: I  ,b (,T )  o
, donde h = 6,6256·10 J·s es la
 5
 exp(hco / kT )  1

constante de Planck, co = 2,998·108 m/s es la velocidad de la luz en el vacío y k = 1,3805·10-23
J/K es la constante de Boltzmann. Como el cuerpo negro es un emisor difuso:
E ,b ( ,T )  I  ,b (  ,T )  C1 , donde C1 y C2 son la primera y segunda

5

 exp(C2 / T )  1

constantes de radiación.

dE ,b (,T ) 
 Ley de desplazamiento de Wien:  0  máx T  C3  2.898m·K .
d
C1
 Ley de Stefan-Boltzmann: E  
d  T 4 , donde  = 5,67·10-8
b 0 5 exp(C2 / T )  1
W/m2·K4 es la constante de Stefan-Boltzmann.
0E  ,b d  E d   E ,b d (T ) 
 Emisión de banda: F (0   )    ,b  
 0
f (T ) .

T 4 0 T 5  
 E  ,b d 0
 Emisividad superficial, : Emisividad direccional espectral (para una  y dirección

I (,  , ,T )
determinadas):  (, , ,T )   ,e . Emisividad direccional total:
 ,
I  ,b (,T )
I ( , ,T )
 ( ,  ,T )  e . Emisividad espectral (hemisférica):  (,T )  E (,T ) .
 
I b (T ) E ,b (,T )
19
  (,T )E  ,b (,T )d


Emisividad total (hemisférica):  (T ) 
E(T )

0

. En un emisor difuso
E b (T ) Eb (T )
la emisividad direccional es independiente de la dirección.
 Absortividad superficial, : Absortividad direccional espectral (para una  y dirección
I  ,i,abs (,  , )
determinadas):  (,  ,  )  . Absortividad (hemisférica) espectral:
 ,
I  ,i (,  , )

G ( )  ()G ()d
 ()   ,abs . Absortividad (hemisférica) total:   Gabs  0  
.

G() 
G
0 G (  )d 
   (, , , sup.)   (T ) . 
 ()E (, 5.800 K)d
 Para radiación solar (Tb = 5.800 K):  S 

0
  ,b
.

0
E  ,b ( , 5.800 K)d
 Reflectividad superficial, : Reflectividad direccional espectral (para una  y dirección
I  ,i,ref ( , , )
determinadas):  (, , )  . Reflectividad (hemisférica) espectral:
 ,
I  ,i (,  , ) 
G ( )   ()G ()d

. Reflectividad (hemisférica) total:   ref  0 
 ,ref G
 ()  
.

G () 
G
0 G ( )d
        1        1 .
 Ley de Kirchhoff (para superficies en el interior de un recinto):
 ...  E b (Ts ) . Como   1  E(Ts)  Eb(Ts). También se cumple:

E1 (Ts ) E2 (Ts )

1 2
1  2
  ...  1   = .
1  2
 Superficie difusa (emisora y receptora difusa): , y , son independientes de la dirección
(, )   = .
 Superficie gris:  y  son independientes de .
 Superficie gris difusa: , y , son independientes de  (gris) y de la dirección (difusa)  
= .
20
TEMA 10. INTERCAMBIO RADIATIVO ENTRE SUPERFICIES
 Factor de forma de radiación (también llamado de configuración, de apariencia, de visión o de

q
vista): F  i j  1 cosi cos j q cosi cos j
dA dA . F  ji  1 dA dA .
ij
Ai J i Ai Ai
A j R 2 i j ji
Aj J j Aj Ai
A j  R 2 i j
Estas dos ecuaciones son válidas para superficies emisoras y reflectoras difusas y con
radiosidad uniforme.
 Relación de reciprocidad: Ai Fij  Aj Fji .
N
 Regla de la suma en un recinto de N superficies:  Fij  1 .

j 1
 Superficie plana o convexa: Fii = 0.

n
 Para una superficie que se puede descomponer en la suma de varias, Aj   Ak , se tiene que:
k 1
n AF k ki
Fi ( j )   Fik y F( j )i  k 1
.
A
n
k 1
k
k 1
 Intercambio neto de radiación entre dos superficies negras:

q  q q AFE  A F E  A F  (T  T 4 ) .
4
ij i j ji i ij bi j ji bj i ij i j
 Transferencia neta de radiación desde la superficie i en un recinto con N superficies negras:
q   A F  (T 4  T 4 ) .
N
i i ij i j
j 1
 Transferencia neta de radiación desde una superficie gris difusa en un recinto:
q  (J  G ) A  (E   G ) A  Ebi  J i , siendo (1 - i) /iAi, la resistencia radiativa

(1   i ) /  i Ai
i i i i i i i i
superficial.
 Transferencia neta de radiación desde una superficie gris difusa en un recinto:
J i  J j , siendo (A F )-1 la resistencia radiativa geométrica.
q  N q   N
i j 1
ij ( A F )1 i ij
j 1 i ij
 Intercambio neto de radiación en un recinto de dos superficies:

 E E .
q1  q2  q12 1   b1 1 b2 1  
1 
 2
1 A1 A1 F12  2 A2
21
 Intercambio neto de radiación entre dos superficies separadas por una cubierta de radiación:
q 1 A1  (T2 4  T. 4 )
Para el caso de N cubiertas de radiación con  iguales
12
1 1  3,1  1  3,2 1
  
1  3,1  3,2 2
1
(incluyendo las superficies extremas): (q12 ) N  (q12 )0 .
N 1
 Superficie rerradiante: superficie idealizada en la que la transferencia de calor neta por
radiación es cero: qi = 0  Ji = Gi = Ebi.
22
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  2wLc
Aleta recta t tanh mLc
de perfil Lc  L  (t / 2) f
mLc
rectangular w m  (2h / kt)1 / 2
L
Aleta recta
de perfil

A f  2w L2  (t / 2) 2 
1/ 2
f 
1 I1 (2mL)
mL I0 (2mL)
triangular w m  (2h / kt)1 / 2
t
L
23
Tabla 2.1. Eficiencia de formas comunes de aletas (continuación).

A f  w C1 L  (L / t) ln(t / L  C 1)
2 2

Aleta recta 
C1  1  (t / L) 1 / 2  f  
2
de perfil
parabólico T
w 2
4(mL)  1
2
1/ 2
1
m  (2h / kt)1 / 2
L
A  2 (r 2  r 2 ) K (mr )I (mr )  I1 (mr1 )K 1 (mr2c )

Aleta f  C 2 1 1 1 2c
anular de t
f 2c 1 I (mr )K (mr )  K (mr )I (mr )
r2c  r2  (t / 2) 0 1 1 2c 0 1 1 2c
perfil (2r / m)
m  (2h / kt)1 / 2 C 2 1 2
rectangular L (r  r 2 )
2c 1
r1
r2
24
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
mLc
cilíndrica m  (4h / kD)1/ 2
L
D
Aleta de A  L
2
 (D / 2) 2  1/ 2
2 I 2 (2mL)
aguja D f
2 f 
mL I1 (2mL)
cónica m  (2h / kt)1 / 2
L
 
L3  L
 C 
Af    
 C3C 4 ln (2DC4 / L) 3

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  (2h / kt)1 / 2
25
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
del sistema
T2
1.1. Esfera
2D
enterrada z
en un T1
zD/2 S
1  (D / 4z)
medio
semiinfinito D
1.2. Esfera
aislado
enterrada
en un
z
medio
semiinfinito T2 T1 T2 2D
zD/2 S
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 26
T2
pág. 55)
1.4.
Conducción T2 4 d  (D /
entre dos T1
S
esferas en D d w/D3 1 
2w) 4  d D 1  (d /
un medio 
w 2 
infinito 
2w) w 
(Bejan)

= 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(2z / D)
longitud L T1
enterrado L  D 2L
en un medio L S
D
z  3D / 2 ln(4z / D)
semiinfinito
2.2. 2
Cilindro de T2 zD S
arc cosh(2z / D)
longitud z
infinita T1 2
S
enterrado   z > 2D ln 2z / D   2z D 2  1 

en un medio  
semiinfinito D 2
(Rohsenow z >> D S
pág. 3-120) ln4z / D
2.3. T2
Cilindro
vertical de L 2L
longitud L
T1
L  D S
ln(4L / D)
enterrado
en un medio D
semiinfinito
2.4.
Conducción
entre dos T2
T1 2L
L  D1 , D2 S
cilindros D d  4w2  D 2  d 2 
paralelos de L  w arc cosh 
longitud L w  2Dd 
en un medio
infinito
27

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(8z / D)
igual longitud D z
y ancho  
infinito T2
T2
2.6. Cilindro
de longitud L
centrado en un wD 2L
sólido de
T1
S
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
Dd S
longitud L en z
 D 2  d 2  4z 2 
el interior de L  D arc cosh 
un cilindro de D
 2Dd 
igual longitud
2.8. Fila
infinita de
cilindros de
longitud 2
S
ln(2L / D)·senh(2z / L)
z>D
infinita en un
medio
semiinfinito
(Rohsenow)
28

Descripción
del sistema
2.9. Fila
infinita de
cilindros de
longitud 2
S
ln(2L / D)·senh(z / L)
infinita en el z>D
plano medio
de una placa
infinita
(Rohsenow)
3.1. Cubo
enterrado en
un medio
L
Ninguna S  8,24L
infinito
(Holman)
4.1. 0,59
  b 
Paralelepípedo 1,685Llog1  
 a
S  0,078 
inmerso en un
Ninguna
medio
b
semiinfinito  
(Holman)  c 
4.2. Agujero a
 5,7
rectangular en
S2b
un medio a>b
 3,5z 
semiinfinito ln 0,25 0,75 
(Rohsenow) a b 
29

Descripción del sistema Esquema Restricciones Factor de forma

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

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

Ninguna S  2D
medio semiinfinito T2
T1
6.2. Disco delgado
z 4,45D
semiinfinito enterrado en S
T2 Ninguna
un medio semiinfinito 1  (D / 5,67z)
(Kreith pág. 112) D
30

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 / 4z)
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 / 4z)
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
z  L 2L
horizontal S
T1 LW ln(4L /W )
delgada de
anchura W
(dimensión z
z0 L
 al dibujo) T2 S
L W ln(4L /W )
enterrada en L
un medio 2L
L  W S
semiinfinito
(Bejan y z  2W ln(2z /W )
Holman)
31

Descripción del Factor de

Esquema Restricciones
sistema forma
7.2. Placa vertical

delgada y larga según
0,24
S  2,38· 
la dimensión  al 1 z L
  12
dibujo enterrada en un 2 L  z 
medio semiinfinito
(Rohsenow)
T1
7.3. Placa horizontal
delgada y larga según
z  L  0,32
la dimensión  al T2
1 z
  12 S  2,94 
dibujo enterrada en un 2 L ·  
medio semiinfinito  z 
L
(Rohsenow)
32
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.
33
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)
T (x, t)  T
Sólido S (x, t) 
semiinfinito x Tini  T sólido
semiinfinito

P (x, t)
T (x, t)  T
Pared plana P(x, t) 
x Tini  T pared
plana
2L1


C (r, t)
T (r, t)  T
Cilindro C(r, t) 
infinito r Tini  T cilindro
infinito
ro


S (x1, t)·P (x2, t)
Placa
S (x1 , t)·P(x2 , t)
semiinfinita x2
x1
2L2
34
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 x2 P(x1 , t)·P(x2 , t)
x1
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 (x3 , t)·P(x1 , t)·P(x2 , 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(x2 , t)·P(x3 , t)
rectangular x
x1 2
2L1
2L2
C (r, t)·P (x, t)
x
Cilindro corto 2L1 C(r, t)·P(x, t)
r
ro
35
TEMA 6. CONVECCIÓN FORZADA EN FLUJO EXTERNO
Tabla 6.1. Coeficientes de la correlación de Hilpert para flujo cruzado sobre un cilindro.
ReD C m
0,4 - 4 1,103 0,330
4 - 40 1,016 0,385
40 - 4.000 0,762 0,466
4.000 - 40.000 0,215 0,618
40.000 - 400.000 0,030 0,805
Tabla 6.2. 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.3. 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
36
Figura 6.1. Nu local para flujo de aire cruzado sobre un cilindro.
Figura 6.2. Disposición de los tubos en configuración alineada y escalonada en un banco de

tubos.
NL
NT
37
Tabla 6.4. Coeficientes de la correlación de Zhukauskas para el flujo cruzado sobre un

banco de tubos.
Configuración ReD,máx C m
Alineado 10 - 102 0,80 0,40
Escalonado 10 - 102 0,90 0,40
Alineado 102 - 103 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
Escalonado (ST / SL < 2) 103 - 2·105 0,35(ST / SL)1/5 0,60
3 5
Escalonado (ST / SL > 2) 10 - 2·10 0,40 0,60
Alineado 2·105 - 2·106 0,021 0,84
5 6
Escalonado 2·10 - 2·10 0,022 0,84
Para ST / SL < 0,7, la transferencia de calor es ineficiente y los tubos alineados no se deben usar.
Tabla 6.5. 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
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.
38
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.
Tabla 7.1. Números de Nusselt para flujo laminar completamente desarrollado en tubos de
diferente sección transversal.
hDh
Nu D 
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
39
CORRELACIONES DE CONVECCIÓN FORZADA EN FLUJO INTERNO
- Problema de longitud de entrada térmica (si xcd,t >> xcd,h).

- Región de entrada 
(x < xcd,t ó Gz -1 < 0,05)
 (Figuras 7.1. y 7.2.): - Problema de longitud de entrada combinada (si O(xcd,t)  O(xcd,h)).
- 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):
Correlación de Hausen.
- Reg. de entrada + c. d.:
- Problema de longitud de entrada combinada (si O(xcd,t)  O(xcd,h)):
Correlación de Sieder y Tate.
- Correlación promedio:
- q s  cte  Nu D  4,36 .
- Región c. d. (L >> xcd,t)  Nu D  cte : - Ts = cte  Nu D  3,66 .
- Tubo no circular  Tabla 7.1.
- 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.
40
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.
TEMA 9. INTRODUCCIÓN A LA RADIACIÓN
Figura 9.1. Distribución de Planck. Potencia Emisiva espectral del cuerpo negro.
9
10
8 T = 5.8 00 K
10
10
7 máx·T = 2.898 m·K
6
10
5 T = 2.0 00 K
10
4 T = 1.0 00 K
10
3
Potencia emisiva 10
2 T = 800 K
espectral, Eb 10
(W/m2·m) 10
1
0 T = 300 K
10
-1 T = 100 K
10
-2
10
-3 T = 50 K
10
-4
10
0.1 0.2 0.4 0.6 1 2 4 6 10 20 40 60 100
longitud de onda,  (m)
41
Tabla 9.1. Funciones de radiación de cuerpo negro.
I  ,b ( ,T ) I  ,b ( ,T )
T
(m·K)
F(0 ) T 5 I  ,b (máx ,T )
(m·K·sr)-1
200 0,000000 0,3750·10-27 0,000000
400 0,000000 0,4903·10-13 0,000000
600 0,000000 0,1040·10-8 0,000014
800 0,000016 0,9911·10-7 0,001372
1.000 0,000321 0,1185·10-5 0,016406
1.200 0,002134 0,5239·10-5 0,072534
1.400 0,007790 0,1344·10-4 0,18608
1.600 0,019718 0,2491·10-4 0,3449
1.800 0,03934 0,3756·10-4 0,5199
2.000 0,06673 0,4934·10-4 0,6831
2.200 0,10089 0,5896·10-4 0,8163
2.400 0,14026 0,6589·10-4 0,9122
2.600 0,1831 0,7013·10-4 0,9709
2.800 0,2279 0,7202·10-4 0,9971
2.898 0,250108 0,7223·10-4 1,000000
3.000 0,2732 0,7203·10-4 0,9971
3.200 0,3181 0,7060·10-4 0,9774
3.400 0,3617 0,6815·10-4 0,9436
3.600 0,4036 0,6504·10-4 0,9004
3.800 0,4434 0,6152·10-4 0,8517
4.000 0,4809 0,5781·10-4 0,8003
4.200 0,5160 0,5404·10-4 0,7481
4.400 0,5488 0,5033·10-4 0,6967
4.600 0,5793 0,4673·10-4 0,6470
4.800 0,6076 0,4331·10-4 0,5996
5.000 0,6337 0,4008·10-4 0,5549
5.200 0,6590 0,3706·10-4 0,5130
5.400 0,6804 0,3424·10-4 0,4741
5.600 0,7010 0,3164·10-4 0,4380
5.800 0,7202 0,2923·10-4 0,4047
6.000 0,7378 0,2701·10-4 0,3740
6.200 0,7541 0,2497·10-4 0,3457
6.400 0,7692 0,2310·10-4 0,3198
6.600 0,7832 0,2138·10-4 0,2960
6.800 0,7961 0,1980·10-4 0,2741
7.000 0,8081 0,1835·10-4 0,2541
7.200 0,8192 0,1703·10-4 0,2357
7.400 0,8295 0,1581·10-4 0,2188
7.600 0,8391 0,1469·10-4 0,2034
7.800 0,8480 0,1366·10-4 0,1891
42
Tabla 9.1. Funciones de radiación de cuerpo negro (continuación).
I  ,b ( ,T ) I  ,b ( ,T )
T
(m·K)
F(0 ) T 5 I  ,b (máx ,T )
(m·K·sr)-1
8.000 0,8563 0,1272·10-4 0,1761
8.500 0,8746 0,1068·10-4 0,1478
9.000 0,8900 0,9015·10-5 0,1248
9.500 0,9031 0,7653·10-5 0,1060
10.000 0,9142 0,6533·10-5 0,09044
10.500 0,9237 0,5605·10-5 0,077600
11.000 0,9319 0,4833·10-5 0,066913
11.500 0,9400 0,4187·10-5 0,057970
12.000 0,9451 0,3644·10-5 0,050448
13.000 0,9551 0,2795·10-5 0,038689
14.000 0,9629 0,2176·10-5 0,030131
15.000 0,9700 0,1719·10-5 0,023794
16.000 0,9738 0,1374·10-5 0,019026
18.000 0,9809 0,9082·10-6 0,012574
20.000 0,9856 0,6233·10-6 0,008629
25.000 0,9922 0,2765·10-6 0,003828
30.000 0,9953 0,1405·10-6 0,001945
40.000 0,9980 0,4739·10-7 0,000656
50.000 0,9990 0,2016·10-7 0,000279
75.000 0,9997 0,4186·10-8 0,000058
100.000 0,999905 0,1358·10-8 0,000019
TEMA 10. INTERCAMBIO RADIATIVO ENTRE SUPERFICIES
Tabla 10.1. Factores de forma radiativos para geometrías bidimensionales.
Geometría Esquema Factor de forma
wi Fij 

(W i  W j ) 2  4
1/ 2


i 2Wi
Placas paralelas centradas L (W  Wi ) 2  4 
1/ 2
 j
j 2Wi
wj
Wi  wi / L W j  w j / L
j
Placas paralelas inclinadas de w  
F 1  
igual anchura y una arista en  ij sen 
común i  2 
w
43
Tabla 10.1. Factores de forma radiativos para geometrías bidimensionales (continuación).
Geometría Esquema Factor de forma

j
 
Placas 1/ 2
perpendiculares wj 1  (w / w )  1  (w / w ) 2
con una arista Fij  j i j i
2
en común
i
wi
wk wj
wi  w j w
Recinto de tres k j Fij  k
lados 2wi
i
wi
  C 2  (R  1) 2 
1  1/ 2
F  

2 
ij
 
j 1/ 2
Cilindros  C 2  (R  1) 2 
i rj  1  R 1 
paralelos de ri
s  (R  1) cos  
radios 180º C C 
 
diferentes  1  R 1 
 (R  1) cos    
180º C C 
  
R  rj / ri S  s / ri C  1  R  S
j
r
r   1 s1 s2 
Cilindro y placa L Fij   tan  tan 1 
paralelos s1
s1  s 2 180º  L L 
i
s2
  D  2 
1/2
s j
Fij  1  1     
  s  
D
Plano infinito y
1/ 2
fila de cilindros  D    
2 2
s D
i   tan1  
 s  180º 
2
 D
44
Figura 10.1. Factor de forma radiativo para dos rectángulos paralelos alineados.
j
L
i Y
X
1
Y / L = inf. Y / L = 10
0.8
Y/L=4
Y/L=2
0.5
Y/L=1
0.3
Y / L = 0,6
0.2
Y / L = 0,4
Fij
0.1
0.08 Y / L = 0,2
0.05
Y / L = 0,1
0.03
0.02
Y / L = 0,05
0.01
0.1 0.2 0.3 0.5 1 2 3 5 10 20 40
X/L
45
Figura 10.2. Factor de forma radiativo para dos discos paralelos coaxiales.
j rj
L
i
ri
1
8
0.9
6 5 4 3
0.8
rj / L = 2
0.7
1,5
0.6
1,25
Fij
0.5
1
0.4
0,8
0.3
0,6
0.2
0,4
0.1
0,3
0
0.1 0.2 0.4 0.6 0.8 1 2 4 6 8 10
L / ri
46
Figura 10.3. Factor de forma radiativo para dos rectángulos perpendiculares con una
arista en común.
j
Z i
X
0.5
Y / X = 0,02
0.45
0,05 0,1
0.4
0,2
0.35
0,4
0.3
0,6
Fij
0.25
1
0.2
1,5
2
0.15
4
0.1
10
0.05 20
0
0.1 0.2 0.4 0.6 0.8 1 2 4 6 8 10
Z/X
47
Tabla 10.2. Intercambio neto de radiación en recintos especiales de dos superficies grises y
difusas.
Intercambio de
Geometría Esquema Condiciones
radiación
plano 1 A (T 4  T 4 )
q12 1 2 
Planos paralelos A1  A2  A
grandes A2, T2, 2 A1, T1, 1
F12  1 1 1
  1
(infinitos) plano 2  
1 2
r1 r
2
A1 r1 A  (T 4  T 4 )
Cilindros  q12  11 1 2
concéntricos A2 r2 1  2 r1

largos (infinitos) F12  1 1 2 r2
A r2 A  (T 4  T 4 )
1 1
r1 2
A2  r22
Esferas 1 1 12q 2
1 1   2  r1 

concéntricas r2
F12  1    r 
1 2  2 
A1, T1, 1 A1
Objeto convexo 0
pequeño en una A2 q  A   (T 4  T 4 )
12 1 1 1 2
cavidad grande F12  1
A2, T2, 2
48
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/m2)
3
(m /s) (W/m·K) (m2/s)
2
100 3,5562 1.032 71,1 2,00 9,34 2,54 0,786

150 2,3364 1.012 103,4 4,426 13,8 5,84 0,758
200 1,7548 1.007 132,5 7,590 18,1 10,3 0,737
250 1,3947 1.006 159,6 11,44 22,3 15,9 0,720
300 1,1614 1.007 184,6 15,89 26,3 22,5 0,707
350 0,9950 1.009 208,2 20,92 30,0 29,9 0,700
400 0,8711 1.014 230,1 26,41 33,8 38,3 0,690
450 0,7740 1.021 250,7 32,39 37,3 47,2 0,686
500 0,6964 1.030 270,1 38,79 40,7 56,7 0,684
550 0,6329 1.040 288,4 45,57 43,9 66,7 0,683
600 0,5804 1.051 305,8 52,69 46,9 76,9 0,685
650 0,5356 1.063 322,5 60,21 49,7 87,3 0,690
700 0,4975 1.075 338,8 68,10 52,4 98,0 0,695
750 0,4643 1.087 354,6 76,37 54,9 109 0,702
800 0,4354 1.099 369,8 84,93 57,3 120 0,709
850 0,4097 1.110 384,3 93,80 59,6 131 0,716
900 0,3868 1.121 398,1 102,9 62,0 143 0,720
950 0,3666 1.131 411,3 112,2 64,3 155 0,723
1.000 0,3482 1.141 424,4 121,9 66,7 168 0,726
1.100 0,3166 1.159 449,0 141,8 71,5 195 0,728
1.200 0,2902 1.175 473,0 162,9 76,3 224 0,728
1.300 0,2679 1.189 496,0 185,1 82 238 0,719
1.400 0,2488 1.207 530 213 91 303 0,703
1.500 0,2322 1.230 557 240 100 350 0,685
49
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/m2)
3 2
(m /s) (W/m·K) 2
(m /s) (K-1)
273 899,1 1.796 385 4.280 147 0,910 47.000 0,70
280 895,3 1.827 217 2.430 144 0,880 27.500 0,70
290 890,0 1.868 99,9 1.120 145 0,872 12.900 0,70
300 884,1 1.909 48,6 550 145 0,859 6.400 0,70
310 877,9 1.951 25,3 288 145 0,847 3.400 0,70
320 871,8 1.993 14,1 161 143 0,823 1.965 0,70
330 865,8 2.035 8,36 96,6 141 0,800 1.205 0,70
340 859,9 2.076 5,31 61,7 139 0,779 793 0,70
350 853,9 2.118 3,56 41,7 138 0,763 546 0,70
360 847,8 2.161 2,52 29,7 138 0,753 395 0,70
370 841,8 2.206 1,86 22,0 137 0,738 300 0,70
380 836,0 2.250 1,41 16,9 136 0,723 233 0,70
390 830,6 2.294 1,10 13,3 135 0,709 187 0,70
400 825,1 2.337 0,874 10,6 134 0,695 152 0,70
410 818,9 2.381 0,698 8,52 133 0,682 125 0,70
420 812,1 2.427 0,564 6,94 133 0,675 103 0,70
430 806,5 2.471 0,470 5,83 132 0,662 88 0,70
50
Tabla C. Propiedades termofísicas del agua saturada.
T P 1/·103 hfg cp ·106 k·103 Pr ·106

2
(K) 3
(bar) (m /kg) (kJ/kg) (J/kg·K) (N·s/m ) (W/m·K) (K-1)
273,15 0,00611 1,000 2.502 4.217 1.750 569 12,99 -68,05
275 0,00697 1,000 2.497 4.211 1.652 574 12,22 -32,74
280 0,00990 1,000 2.485 4.198 1.422 582 10,26 46,04
285 0,01387 1,000 2.473 4.189 1.225 590 8,81 114,1
290 0,01917 1,001 2.461 4.184 1.080 598 7,56 174,0
295 0,02617 1,002 2.449 4.181 959 606 6,62 227,5
300 0,03531 1,003 2.438 4.179 855 613 5,83 276,1
305 0,04712 1,005 2.426 4.178 769 620 5,20 320,6
310 0,06221 1,007 2.414 4.178 695 628 4,62 361,9
315 0,08132 1,009 2.402 4.179 631 634 4,16 400,4
320 0,1053 1,011 2.390 4.180 577 640 3,77 436,7
325 0,1351 1,013 2.378 4.182 528 645 3,42 471,2
330 0,1719 1,016 2.366 4.184 489 650 3,15 504,0
335 0,2167 1,018 2.354 4.186 453 656 2,88 535,5
340 0,2713 1,021 2.342 4.188 420 660 2,66 566,0
345 0,3372 1,024 2.329 4.191 389 668 2,45 595,4
350 0,4163 1,027 2.317 4.195 365 668 2,29 624,2
355 0,5100 1,030 2.304 4.199 343 671 2,14 652,3
360 0,6209 1,034 2.291 4.203 324 674 2,02 697,9
365 0,7514 1,038 2.278 4.209 306 677 1,91 707,1
370 0,9040 1,041 2.265 4.214 289 679 1,80 728,7
373,15 1,0133 1,044 2.257 4.217 279 680 1,76 750,1
375 1,0815 1,045 2.252 4.220 274 681 1,70 761
380 1,2869 1,049 2.239 4.226 260 683 1,61 788
385 1,5233 1,053 2.225 4.232 248 685 1,53 814
390 1,794 1,058 2.212 4.239 237 686 1,47 841
400 2,455 1,067 2.183 4.256 217 688 1,34 896
410 3,302 1,077 2.153 4.278 200 688 1,24 852
420 4,370 1,088 2.123 4.302 185 688 1,16 1.010
430 5,699 1,099 2.091 4.331 173 685 1,09
440 7,333 1,110 2.059 4.360 162 682 1,04
450 9,319 1,123 2.024 4.400 152 678 0,99
460 11,71 1,137 1.989 4.440 143 673 0,95
470 14,55 1,152 1.951 4.480 136 667 0,92
480 17,90 1,167 1.912 4.530 129 660 0,89
490 21,83 1,184 1.870 4.590 124 651 0,87
500 26,40 1,203 1.825 4.660 118 642 0,86
51
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 1.490,5 1.490,5 1,0000
1,80 2,9422 3,1075 0,94681 9,00 4.051,5 4.051,5 1,0000
1,90 3,2682 3,4177 0,95624 10,00 11.013 11.013 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
2w
erf w  e u du
2
 0
erfc w  1  erf w
52
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
53
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
54
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
In1 ( x)  In1 ( x)  (2n / x)In ( x)
55
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
56

Formulario Transferencia de Calor

Cargado por

Copyright:

Formatos disponibles

Formulario Transferencia de Calor

Cargado por

Información del documento

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

Formulario Transferencia de Calor

Cargado por

Copyright:

Formatos disponibles

INSTITUTO TECNOLÓGICO SUPERIOR

FÓRMULAS, TABLAS Y FIGURAS DE

Alumno: Enrique Rodríguez Hernández

Docente: Ing. Guillermo Córdova

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

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

almacenar energía térmica.

 En el caso de transmisión unidimensional en régimen estacionario y sin generación de

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

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

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

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

 Resistencia térmica de conducción para una pared cilíndrica: R (Ts1  Ts 2 ) ln(r2 / r1 )

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

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

 (x) cosh m(L  x) q  M tanh mL

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

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,

se puede expresar la efectividad como:   Rt ,b .

tiene:   M tanh mL  tanh mL .

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 enfriamiento de Newton ( q ()(m,n)  hx·1(T  Tm,n ) .

TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO

 Longitud característica: Lc = Vol/As. Para una pared plana de espesor 2L  Lc = L, y para un

0). Los valores de C1 y 1 se obtienen de la Tabla 4.1.

2L sometida a convección: Q(t)  1  sen 1  * (0, t) , siendo Uo = Uini = cpVol(Tini - T) =

T (x, t)  T T (x, t)  T T (r, t)  T

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

Condición de frontera Distribución de temperaturas

la energía almacenada: E& alm  cpV m,n .

TEMA 5. INTRODUCCIÓN A LA CONVECCIÓN

 Ley de enfriamiento de Newton: q  h(Ts  T ) ; q  hA(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.

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

 Ecuación de conservación de la energía en la capa límite: u T T

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:

 0,05 y Rex·Pr  100.

 Espesor de la capa límite térmica turbulenta:  t ,turb   turb .

variación de la temperatura superficial local se obtiene con: T (x)  T  qs

 Número de Reynolds para flujo cruzado sobre un cilindro: ReD = V·D/.

20 se aplica un factor de corrección tal que Nu D  C2 Nu D , donde C2 está dado en

 Diferencia de temperaturas media logarítmica: T ml  (Ts  Tent )  (Ts  Tsal ) .

TEMA 7. CONVECCIÓN FORZADA EN FLUJO INTERNO

 Expresión de la velocidad media en función del flujo másico integrado en la sección

m&  u(r, x)dAc

Ac Ac ro2 0

U&  E&  u(r)c T (r, x)dA  mc T (x)  T

  ucvTdAc  2 u(r)T (r, x)rdr .

 Bajo condiciones térmicas completamente desarrolladas se cumple:

 Variación axial de la temperatura media para el caso de temperatura superficial constante:

 h A T ; T Tsal  Tent (Ts  Tm,sal )  (Ts  Tm,ent )