CUADERNO DE FÓRMULAS, TABLAS, FIGURAS Y PROBLEMAS PROBLEMAS DE TRANSFERENCIA DE CALOR
Juan Carlos Ramos González Doctor Ingeniero Industrial Raúl Antón Remírez Doctor Ingeniero Industrial Diciembre de 2009
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 Diciembre de 2009
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
ii
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Í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...................................................................................9 Tema 6. Convección forzada en flujo externo..................................................................... 11 Tema 7. Convección forzada en flujo interno...................................................................... 14 Tema 8. Convección libre o natural .....................................................................................17 Tema 9. Introducción a la radiación ....................................................................................19 Tema 10. Intercambio radiativo entre superficies ...............................................................22
Tablas y Figuras Tema 2. Conducción unidimensional en régimen estacionario ...........................................24 Tema 3. Conducción bidimensional en régimen estacionario.............................................27 Tema 4. Conducción en régimen transitorio........................................................................34 Tema 6. Convección forzada en flujo externo..................................................................... 37 Tema 7. Convección forzada en flujo interno...................................................................... 42 Tema 8. Convección libre o natural .....................................................................................46 Tema 9. Introducción a la radiación ....................................................................................48 Tema 10. Intercambio radiativo entre superficies ...............................................................50 Tablas de propiedades termofísicas y de funciones matemáticas........................................55 Alfabeto griego .................................................................................................................... 62
i
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
ii
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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]. • Ley de Fourier: q x′′ = −k
dT dx
. q x = q x′′ ⋅ A . En condiciones de régimen estacionario y con una
distribución lineal de temperaturas: q x′′ = − k
dT dx
T 2 − T 1
= − k
T 1 − T 2
= k
L
L
= k
∆T 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 = σ T s4 . • Constante de Stefan-Boltzmann: σ = 5,67·10-8 W/m2·K 4. • 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 = εσ T s4 , 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
4 alrededores: G = σ T alr .
• Intercambio ′′ = q rad
q A
de
radiación
para
una
4 ). = ε E b (T s ) − α G = εσ (T s4 − T alr
′′ = hrad (T s − T alr ) , q rad
(
superficie También
gris se
y puede
difusa
( α =
expresar
ε ):
como:
siendo hrad el coeficiente de transferencia de calor por radiación:
2 hrad = εσ (T s + T alr ) T s2 + T alr
).
• Principio de conservación de la energía en un volumen de control formulado en un instante de & & & tiempo (t ): E ent + E gen − E sal =
dE alm dt
& = E alm .
1
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• 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 . r r r ∂T r ∂T r ∂T r ′ ′ ′ ′ ′ ′ • Ley de Fourier vectorial: q = −k ∇T = −k i + j + k = i q x + j q y + k q z ′′ . x y ∂ ∂ ∂ z r
• Capacidad térmica volumétrica: ρ c p [J/m3·K]. Mide la capacidad de un material para
almacenar energía térmica. • Difusividad térmica: α =
k ρ c p
[m2/s]. Mide la capacidad de un material para conducir energía
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 • Ecuación de difusión de calor vectorial: ∇·(k ∇T ) + q& = ρ c p
∂T . ∂t
• En el caso de transmisión unidimensional en régimen estacionario y sin generación de
energía:
d dT 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 ( d q 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 ∂ ∂T ∂T k + k + q& = ρ c p . + 2 ∂r r ∂φ ∂φ ∂ z ∂ z ∂t
1 ∂ ∂T kr r ∂r
• 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 1 ∂ ∂T ∂ ∂T ∂T k + 2 2 . + 2 2 ksenθ + q& = ρ c p ∂r r sen θ ∂φ ∂φ r sen θ ∂θ ∂θ ∂t
1 ∂ 2 ∂T kr
r 2 ∂r
• Condición de contorno de primera clase o de Dirichlet: superficie mantenida a temperatura
constante, T ( x = 0, t ) = T s. 2
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Condición de contorno de segunda clase o de Neumann: flujo de calor fijo o constante en la ∂T . Un caso especial es la superficie perfectamente aislada o ∂ x x=0
superficie, q s′′( x = 0) = − k adiabática,
∂T =0. ∂ x x=0
• Condición de contorno de tercera clase o de Fourier: corresponde a la transferencia de calor ′′ , superficie = qconv ′′ . Si el fluido está en contacto con la por convección en la superficie, q cond ∂T = h[T ∞ − T ( x = 0, t )]. Si ∂ x x=0
superficie de la pared donde está el origen de coordenadas: − k
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 • Resistencia térmica de conducción para pared plana: Rt ,cond = • Resistencia térmica de convección: Rt ,conv = • Resistencia térmica de radiación. Rt ,rad =
T s − T ∞ q
T s − T alr qrad
=
=
1 hA
1 hr A
T s1 − T s 2
L
=
q x
.
kA
.
.
• Coeficiente global de transferencia de calor, U : q x = UA∆T . Rtot = ∑ Rt =
∆T q
=
1 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
x0
A( x)
la transferencia de calor, q x, es una constante independiente de x): q x ∫ • Resistencia térmica de conducción para una pared cilíndrica: Rt ,cond = • Resistencia térmica de convección para una pared cilíndrica: Rt ,conv =
3
T
= −∫ k (T )dT . T 0
(T s1 − T s 2 ) qr
1 Ah
=
=
1 2π rLh
ln(r 2 / r 1 ) . 2π Lk .
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Resistencia Rt ,cond =
térmica
(T s1 − T s 2 ) qr
=
de
conducción
para
una
pared
esférica:
1 1 1 − . 4π k r 1 r 2
• Resistencia térmica de convección para una pared esférica: Rt ,conv =
1 Ah
=
1 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
• Generación de energía térmica por unidad de volumen: q& = e& gen • Ecuación de calor para una aleta:
& E gen W . = Vol m 3
d 2T 1 dAc dT 1 h dA s dx 2
(T − T ∞ ) = 0 . + − 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
( hAc [T ( L) − T ∞ ] = − kAc θ ( x ) θ b
=
dT dx
): x = L
cosh m( L − x) + (h / mk )senh m( L − x) cosh mL + (h / mk )senh mL
siendo θ ( x) = T ( x) − T ∞ , θ b = T b − T ∞ , m 2 =
q f = M hP kAc
senh mL + (h / mk ) cosh mL cosh mL + (h / mk )senh mL
, M = hPkAc θ b , P el perímetro y Ac
el área transversal. • Caso B, extremo adiabático ( θ ( x ) θ b
=
d θ dx
= 0 ): x = L
cosh m( L − x) cosh mL
q f = M tanh mL
• Caso C, extremo con temperatura establecida ( θ ( x = L) = θ L): θ ( x ) θ b
=
(θ L / θ b )senh mx + senh m( L − x) senh mL
q f = M
cosh mL − (θ L / θ b ) senh mL
• Caso D, aleta muy larga ( L → ∞ y θ L → 0, aplicable si m· L > 2,65): θ ( x) = θ b e − mx
q f = M =
4
hPkAc θ b
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• La efectividad de una aleta se define como la razón entre la transferencia de calor de la aleta y
la transferencia de calor que existiría sin la aleta: ε f =
q f hAc ,bθ b
, siendo Ac,b el área de la
sección transversal de la base de la aleta. El uso de aletas sólo se justifica cuando ε f ≥ 2. • Resistencia térmica de una aleta: Rt , f =
θ b q f
.
• 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: ε f =
Rt ,b 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: η f =
q f q máx
=
q f hA f θ b
, siendo A f la superficie total de la aleta.
• Teniendo en cuenta la ecuación que define la resistencia térmica de una aleta, se puede
expresar ésta en función de su eficiencia: Rt , f =
1 hA f η f
.
• Para el caso de una aleta recta de sección transversal uniforme y con su extremo adiabático se
tiene: η f =
M tanh mL hPLθ b
=
tanh mL 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 ( I 0, I 1, K 0 y K 1) 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. • Dispositivo de varias aletas. Eficiencia global de la superficie: η o =
qt q máx
=
qt hAt θ b
, siendo qt
la transferencia total de calor de la superficie total, At , que es la asociada a la superficie de las 5
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
aletas, A f , 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: qt = Nq f + qb = N η f hA f θ b + hAbθ b ⇒ qt = η o hAt θ b ⇒ η o = 1 −
• Resistencia térmica efectiva del dispositivo de aletas: Rt ,o =
θ b qt
NA f
=
At
(1 − η f ) .
1 hAt η o
.
TEMA 3. CONDUCCIÓN BIDIMENSIONAL EN RÉGIMEN ESTACIONARIO • Factor de forma de conducción para sistemas bidimensionales, S : q = Sk (T 1 − T 2 ) . Se obtiene
de la Tabla 3.1. • Resistencia de conducción bidimensional: Rt ,cond ( 2 D ) =
1 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
convección que se evalúan con la ley de Fourier ( q ( m −1,n )→( m,n ) = k (∆ y·1)
T m−1,n − T m, n
∆ x
) o la ley
de enfriamiento de Newton ( q( ∞ )→( m,n ) = h(∆ x·1)(T ∞ − T m,n ) .
TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO • Definición general del número de Biot: Bi =
Rt,int Rt ,ext
=
• Número de Biot para un sólido con convección: Bi =
∆T int . ∆T ext hLc k
=
Rt ,cond Rt ,conv
.
• Longitud característica: Lc = Vol / A s. Para una pared plana de espesor 2 L sometida a
convección simétrica en su superficie ⇒ Lc = L, y para un cilindro largo o una esfera de radio r o ⇒ Lc = r o.
6
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Número de Fourier: Fo =
α ·t 2
Lc
.
• El método de la resistencia interna despreciable es aplicable cuando Bi =
hLc k
< 0,1 .
• Distribución de temperaturas temporal en un sólido en el que se puede aplicar el método de la
resistencia interna despreciable: θ (t ) θ ini
=
hA s = exp − T ini − T ∞ ρ Volc p
T (t ) − T ∞
t . τ t
t = exp −
• Constante de tiempo térmica: τ t =
1 hA s
·ρ Volc p = Rt ·C t , siendo Rt la resistencia a la
transferencia de calor por convección y C t 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 Q(t ) =
de t
∫0
la
resistencia
interna
despreciable
durante
un
tiempo
t
será:
t
∫0
qdt = hA s θ (t ) dt ⇒ Q (t ) = ρ Volc pθ ini [1 − exp(− t τ t )] = U ini [1 − exp(− t τ t )] .
• 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 2 L sometida a convección: θ * ( x*, t *) =
θ o
T ( x, t ) − T ∞ T ini − T ∞
* ( x* = 0, t ) =
= C 1 exp( −ζ 12 Fo) cos(
T ( x = 0, t ) − T ∞ T ini − T ∞
180 π
ζ 1 x*) = θ o
* cos(
180 π
ζ 1 x*) ,
siendo
= C 1 exp( −ζ 12 Fo) la temperatura del plano medio ( x* = x/ L =
0). Los valores de C 1 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
2 L sometida a convección: T ∞)
Q (t ) U o
180 ζ 1 π θ * (0, t ) , siendo U = U = ρ c Vol (T o ini p ini o
sen = 1−
ζ 1
= U má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. • Solución analítica aproximada con el primer término (aplicable cuando Fo > 0,2) de la
distribución de temperaturas en un cilindro largo de radio r o sometido a convección: 7
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras θ * ( r *, t *) =
θ o
T (r , t ) − T ∞
= C 1 exp( −ζ 12 Fo) J 0 (ζ 1 r *) = θ o * J 0 (ζ 1 r *) ,
T ini − T ∞
* (r * = 0, t ) =
T o − T ∞ T ini − T ∞
siendo
= C 1 exp( −ζ 12 Fo) la temperatura del eje central ( r * = r /r o = 0) y J 0
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 r o sometido
a convección:
Q(t ) U o
= 1−
2θ o * (0, t ) ζ 1
J 1 (ζ 1 ) ,
siendo J 1 la función de Bessel de
primera clase de orden uno cuyos valores se encuentran en la Tabla G y U o = U ini = ρ c pVol (T ini - T ∞) = U má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. • Solución analítica aproximada con el primer término (aplicable cuando Fo > 0,2) de la
distribución de temperaturas en una esfera de radio r o sometida a convección: θ * ( r *, t *) =
T (r , t ) − T ∞
= C 1 exp( −ζ 12 Fo)
T ini − T ∞
siendo θ o * (r * = 0, t ) =
T o − T ∞ T ini − T ∞
1
180 1 180 sen ( ζ 1 r *) = θ o * sen ( ζ 1 r *) , ζ 1r * π ζ 1 r * π
= C 1 exp( −ζ 12 Fo) la temperatura del eje central ( r * = r /r o =
0). • Transferencia total de energía en función del tiempo desde o hacia una esfera de radio r o
sometida a convección:
Q(t ) U o
= 1−
3θ o * ζ 13
180 180 sen( π ζ 1 ) − ζ 1 cos( π ζ 1 ) , siendo U o = U ini =
ρ c pVol (T ini - T ∞) = U má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 2 L o un cilindro infinito de radio r o: S ( x, t ) =
T ( x, t ) − T ∞ T ini − T ∞
; P ( x, t ) = Sólido semiinfinito
8
T ( x, t ) − T ∞ T ini − T ∞
; C (r , t ) = Pared plana
T ( r , t ) − T ∞ T ini − T ∞
. Cilindro infinito
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• En un sólido semiinfinito la condición de frontera interior es T ( x→∞, t ) = T ini y la condición
inicial es T ( x, 0) = T ini. Las soluciones analíticas para tres condiciones de frontera exterior son: Condición de frontera
Distribución de temperaturas T ( x, t ) − T s x = erf T ini − T s α t 2 k (T s − T ini ) ∂T = q s′′ (t ) = − k ∂ x x=0 (πα t )1 / 2
Temperatura superficial constante: T (0, t ) = T s
Condición de frontera Flujo de calor superficial constante: q s′′ = q o′′
Distribución de temperaturas x 2 q o′′ x 2qo′′ (α t / π )1 / 2 x − T ( x, t ) − T ini = exp − erfc − α k t k 4 α t 2 T ( x, t ) − T ini
x = erfc − T ∞ − T ini 2 α t hx h 2α t x h α t − exp + 2 erfc + k k k 2 α t
Convección superficial: ∂T − k = h[T ∞ − T (0, t )] ∂ x x=0
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 +
& la energía almacenada: E alm = ρ c pV
T m p,n1 − T m p, n
∆t
.
• Criterio de estabilidad en el MDF explícito: el coeficiente asociado con el nodo de interés en
el tiempo anterior (coeficiente de T m p,n ) debe ser mayor o igual que cero. Así se obtiene un valor límite para Fo =
α ∆t
(∆ x) 2
, del que se obtiene el máximo valor permisible de ∆t .
TEMA 5. INTRODUCCIÓN A LA CONVECCIÓN • Ley de enfriamiento de Newton: q ′′ = h(T s − T ∞ ) ; q = h A(T s − T ∞ ) . • Coeficiente de transferencia de calor por convección local, h o promedio, h [W/m2·K].
9
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Relación q=
∫
A s
entre
los
q ′′dA s = (T s − T ∞ )
placa plana: h =
1 L
∫
A s
coeficientes
de
hdA s = h A s (T s − T ∞ )
convección ⇒
h =
1 A s
local
∫
A s
hdA s
y
promedio:
. Para flujo sobre una
L
∫0 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 ( T s - T ( y))/(T s - T ∞) = 0,99. • Relación del coeficiente de convección en la capa límite: h = • Número de Reynolds: Re x =
ρ u ∞ x µ
=
u ∞ x
ν
− k f ∂T / ∂ y y =0 T s − T ∞
.
.
• Número de Reynolds crítico para el inicio de la turbulencia en flujo externo: Re x,c = 5·105. • Expresión diferencial de la ecuación de conservación de la masa o de continuidad: ∂ ρ ∂ ( ρ u ) ∂( ρ v) kg + + = 0 3 . ∂t ∂ x ∂ y m ·s • 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 kg N + v + = − + µ 2 − + + µ + + X 3 = 2 2 ∂ y ∂t ∂ x ∂ x ∂ x 3 ∂ x ∂ y ∂ y ∂ y ∂ x m ·s m ∂ x ∂v ∂v ∂v ∂ p ∂ ∂v 2 ∂u ∂v ∂ ∂u ∂v ρ u + v + = − + µ 2 − + + µ + + Y . ∂ y ∂t ∂ y ∂ y ∂ y 3 ∂ x ∂ y ∂ x ∂ y ∂ x ∂ x
ρ u
• Expresión diferencial de la ecuación de conservación de la energía: ∂ ∂T ∂ ∂T ∂( pu ) ∂( pv) − = k + k + q& − Xu − Y v + µ Φ − ∂ x ∂ x ∂ y ∂ y ∂ x ∂ y V 2 ∂ V 2 ∂ V 2 W ∂ + ρ v e + gy + + ρ e + gy + ρ u e + gy + 2 ∂ y 2 ∂t 2 m 3 ∂ x
, siendo V 2 = u2 + v2.
• 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 ∂ ∂T + v = k + k + µ Φ + q& . ∂ y ∂ x ∂ x ∂ y ∂ y ∂ x
ρ c p u
10
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
∂u ∂v 2 ∂u 2 ∂v 2 • Disipación viscosa: µ Φ = µ + + 2 + . ∂ x ∂ y ∂ y ∂ x • 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 en la capa límite de velocidad y en >> , , >> ∂ y ∂ x ∂ y ∂ x ∂ y ∂ x
( q& = 0 ).Además: u >> v y la capa límite térmica.
• Ecuación de conservación de la masa o de continuidad en la capa límite:
∂u ∂v + =0 . ∂ x ∂ y
• Ecuaciones de balance de la cantidad de movimiento o del momento lineal en la capa límite:
1 ∂ p ∂u ∂u ∂ 2u u +v =− + ν 2 ρ ∂ x ∂ x ∂ y ∂ y
y
∂ p = 0. ∂ y
∂T ∂T ∂ 2T ν • Ecuación de conservación de la energía en la capa límite: u +v = α 2 + c p ∂ x ∂ y ∂ y • Número de Prandtl: Pr =
ν α
• Número de Nusselt: Nu =
2
∂u . ∂ y
.
hL k f
=
∂T * . ∂ y * y*=0
• Las formas adimensionales de las soluciones de la capa límite adoptan la siguiente forma: Nu = f ( x*, Re L , Pr ) y Nu =
h L k f
= f ( Re L , Pr ) .
• Relación entre los espesores de las capas límites hidrodinámica y térmica en régimen laminar: δ δ t
≈ Pr 1 / 3 .
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: T f =
T s + T ∞
2
.
11
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
5
• Espesor de la capa límite laminar: δ lam ( x ) =
u ∞ /ν x
=
5 x Re x
.
• Correlación de convección local para el flujo laminar sobre una placa plana con temperatura
superficial constante: Nu x =
h x x k
= 0,332 Re x1 / 2 Pr 1 / 3 . Con propiedades a T f y Pr ≥ 0,6.
• Relación entre los espesores de las capas límites de velocidad y térmica:
δ lam δ t ,lam
≈ Pr 1 / 3 .
• Correlación de convección promedio para el flujo laminar sobre una placa plana con
temperatura superficial constante: Nu x =
h x x k
= 0,664 Re x1 / 2 Pr 1 / 3 . Con propiedades a T f y Pr ≥
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 T s = cte: 4/5
1/ 3
Nu x = 0,0296 Re x Pr
. Con propiedades a T f y 0,6 < Pr < 60.
• Para condiciones de capa límite mezclada (laminar y turbulenta) se trabaja con el coeficiente
de convección promedio: h L =
1 x
L ∫0 hlam dx + ∫ x hturb dx . L c
c
• Correlación de convección promedio para capa límite mezclada (laminar y turbulenta) sobre
una placa plana con T s = cte: Nu L = (0,037 Re L4 / 5 − 871) Pr 1 / 3 . Con propiedades a T f y 0,6 < Pr < 60 5 8 5·10 < Re L < 10 . Re x = 5·10 5 c
• 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 Re L >> Re x,c), sobre una placa 0,6 < Pr < 60 plana con T s = cte: Nu L = 0,037 Re L4 / 5 Pr 1 / 3 . Con propiedades a T f 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 x1 / 2 Pr 1 / 3 , con Pr ≥ 0,6 y propiedades a T f . 12
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• 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 T f . • 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 s ( x) = T ∞ +
q s′′
.
h x ( x )
• Número de Reynolds para flujo cruzado sobre un cilindro: Re D = 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 Re D. La Tabla 6.4 da los valores de las constantes para cilindros no circulares. Las propiedades se evalúan a T f . Válida para fluidos con Pr ≥ 0,7. Pr • Correlación de Zhukauskas para flujo cruzado sobre un cilindro: Nu D = CRe Pr Pr s m D
1/ 4
n
.
0,7 < Pr < 500 1 < Re < 10 6 . Los valores de las constantes C y m se dan en D
n = 0,37 si Pr ≤ 10 Con y n Pr = > 0 , 36 si 10
la Tabla 6.5 en función de Re D. Las propiedades se evalúan a T ∞, excepto Pr s a T s. • Correlación de Churchill y Bernstein para flujo cruzado sobre un cilindro: 1/ 2
Nu D = 0,3 +
• Correlación
1/ 3
0,62 Re D Pr [1 + (0,4 / Pr ) 2 / 3 ]1 / 4 de
Zhukauskas 1/ 4
0 ,36 Pr Nu D = CRe Dm,máx Pr Pr s
Re D 5 / 8 1 + 282 . 000
para
flujo
a
4/5
. Con propiedades a T f y Re D· Pr > 0,2. través
de
un
banco
de
tubos:
N L ≥ 20 . Con 1.000 < Re D,máx < 2·10 6 . Las constantes C y m se dan 0,7 < Pr < 500
en la Tabla 6.6. Las propiedades se evalúan a T = (T ent + T sal ) / 2 , excepto Pr s a T s. Para N L < 20 se aplica un factor de corrección tal que Nu D
N L < 20
= C 2 Nu D
N L ≥ 20
, donde C 2 está dado en
la Tabla 6.7. • Re D,máx se define en función de la velocidad máxima del fluido dentro del banco de tubos. • S T es el espaciado transversal y S L el espaciado longitudinal (distancias entre centros de
tubos). 13
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Para la configuración alineada la velocidad máxima se da en el plano transversal entre dos
tubos verticales y su valor es V máx =
S T S T − D
V .
• Para la configuración escalonada se utiliza la misma expresión si la velocidad máxima se da
en el plano transversal. Pero si se da en el plano diagonal la expresión es V máx =
S T
2( S D − D)
V .
La velocidad máxima ocurre en el plano diagonal si se cumple la siguiente condición (ver 2 S T 2 Figura 6.2): 2( S D − D ) < ( S T − D) ⇒ S D = S L + 2
• Diferencia de temperaturas media logarítmica: ∆T ml =
1/ 2
<
S T + D
2
.
(T s − T ent ) − (T s − T sal ) T s − T ent T T − s sal
.
ln • Cálculo de la temperatura de salida del flujo:
π DN h = exp − ρ VN T S T c p , donde N es el T s − T ent T s − T sal
número total de tubos y N T 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∆T ml .
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
velocidad media del fluido sobre la sección transversal del tubo: Re D =
ρ u m D µ
=
u m D
ν
& = ρ u m Ac , para un tubo circular el número de Reynolds se puede expresar: Re D = m
. Como
4m& π Dµ
.
• Número de Reynolds crítico para el inicio de la turbulencia en flujo interno: Re D,c = 2.300. x • Longitud hidrodinámica de entrada para flujo laminar: cd ,h ≈ 0,05 Re D . D lam x • Longitud hidrodinámica de entrada para flujo turbulento: 10 ≤ cd ,h ≤ 60 . D turb
14
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Expresión de la velocidad media en función del flujo másico integrado en la sección
transversal: u m ( x) =
& m
ρ Ac
∫ =
Ac
ρ u ( r , x) dAc ρ Ac
=
2 r o2
r o
∫0 u(r , x)rdr .
x • Longitud de entrada térmica para flujo laminar: cd ,t ≈ 0,05 Re D Pr . D lam x • Longitud de entrada térmica para flujo turbulento: cd ,t = 10 . D turb
• Temperatura media definida en función de la energía térmica transportada por el fluido: t = ∫ ρ u ( r )c v T ( r , x ) dAc
& = E & U
• Bajo
Ac
condiciones
∫ = m& c T ( x ) ⇒ T ( x ) =
Ac
v
térmicas
m
m
& cv m
completamente
dT s ∂ T s ( x) − T (r , x) ∂T 0 . Además: = = ∂ x T s ( x ) − T m ( x ) cd ,t ∂ x cd ,t dx T − T dT m ∂T = s ∂ x cd ,t T s − T m dx
ρ uc vTdAc
2 u m r o2
desarrolladas =
cd ,t
=
dT m dx
r o
∫0 u (r )T (r , x)rdr . se
cumple:
≠ f (r ) para q s′′ = cte y cd ,t
= f (r ) para T s = cte. 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: qconv = m& c p (T m, sal − T m,ent ) . • Variación axial de la temperatura media para el caso de flujo de calor superficial constante: T m ( x) = T m ,ent +
′ q s′ P x . & c p m
• Variación axial de la temperatura media para el caso de temperatura superficial constante: T s − T m ( x ) T s − T m,ent
Px h = exp − m& c p .
• La transferencia total de calor se expresa en función de la diferencia de temperaturas media
logarítmica: qconv = h A s ∆T ml ; ∆T ml =
(T − T m, sal ) − (T s − T m,ent ) ∆T sal − ∆T ent . = s ln(∆T sal / ∆T ent ) (T s − T m, sal ) ln (T s − T m ,ent )
15
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Caso de un tubo rodeado de un fluido externo (convección interna y externa simultáneas): U A ∆T sal T ∞ − T m, sal ∆T ml 1 q conv = U A s ∆T ml = ; . = = exp − s = exp − & c p & c p Rtot m m R ∆T ent T ∞ − T m,ent tot • Correlación de convección local para flujo laminar en la región completamente desarrollada
con flujo de calor superficial constante: Nu =
hD k
=
48 = 4,36 . Propiedades a T m. 11
• Correlación de convección local para flujo laminar en la región completamente desarrollada
con temperatura superficial constante: Nu =
hD k
= 3,66 . Propiedades a T m.
• 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 Nu D =
h D k
desarrollado) = 3,66 +
y
con
temperatura
superficial
constante:
0,0668( D / L) Re D Pr . Propiedades a T m = (T m ,ent + T m, sal ) / 2 . 1 + 0,04[( D / L) Re D Pr ]2 / 3
• Correlación de Sieder y Tate para flujo laminar interno con longitud de entrada combinada y 1/ 3
µ con temperatura superficial constante: Nu D = 1,86 L / D µ s Re D Pr
0 ,14
. Con propiedades a
0,48 < Pr < 16.700 . µ T m = (T m,ent + T m , sal ) / 2 , excepto µ s a T s y < < 0 , 0044 9 , 75 µ s
• 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 ( T s > T m), n = 0,3 para enfriamiento
0,7 ≤ Pr ≤ 160 (T s < T m), las propiedades evaluadas a T m y Re D ≥ 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
16
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
µ temperatura superficial constante: Nu D = 0,027 Re D Pr µ s 4/5
0 ,14
1/ 3
. Con propiedades a T m,
0,7 ≤ Pr ≤ 16.700 excepto µ s a T s y Re D ≥ 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 T m = (T m,ent + T m , sal ) / 2 . • Para tubos no circulares se trabaja con el diámetro hidráulico: Dh =
4 Ac P mojado
, donde Ac es el
área de la sección transversal y P mojado el perímetro mojado. Las expresiones del número de Reynolds para el diámetro hidráulico son: Re D = h
ρ u m Dh µ
=
u m Dh
ν
=
4m& µ 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.
TEMA 8. CONVECCIÓN LIBRE O NATURAL • Número de Grashof: Gr x =
g β (T s − T ∞ ) x 3
ν 2
.
• Relación entre la convección forzada y la convención libre: si Gr L / Re L2 <<1 , la convección
libre se desprecia frente a la forzada; si Gr L / Re L2 >> 1 , la forzada se desprecia frente a la libre. • Soluciones de similitud para la convección libre laminar sobre una superficie vertical.
Número
de
Nusselt
Gr = x Nu x = k 4 hx
local:
17
1/ 4
· f ( Pr ) ,
siendo
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
0,75 Pr 1 / 2 f ( Pr ) = . 1/ 4 1/ 2 (0,609 + 1,221 Pr + 1,238 Pr )
Número
de
Nusselt
promedio:
1/ 4
4 Gr 4 = L · f ( Pr ) = Nu L . Nu L = k 3 4 3 h L
• Número de Rayleigh: Ra x = Gr x Pr =
g β (T s − T ∞ ) x 3
να
.
• Transición entre la capa límite laminar y la turbulenta en placas verticales: Gr x,c ≈ 109 ⇒ Ra x,c / Pr ≈ 109.
• Correlación de Churchill y Chu para la convección libre sobre una superficie vertical a
temperatura
constante
aplicable
para
todo
Ra L:
2
0,387 Ra L1 / 6 h L Nu L = = 0,825 + . Propiedades calculadas a T f = (T s + T ∞)/2. 9 / 16 8 / 27 k [1 + (0,492 / Pr ) ] • Correlación para la convección libre sobre una superficie vertical a temperatura constante
aplicable al flujo laminar: Nu L =
h L k
= 0,68 +
0,670 Ra L1 / 4 9 / 16 4 / 9
[1 + (0,492 / Pr ) ]
con Ra L ≤ 109.
Propiedades calculadas a T f = (T s + 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 ( T s - T ∞) aumentará con x. Las correlaciones anteriores son aplicables en este caso si Nu L y Ra L se definen en términos de la diferencia de temperaturas en el punto medio de la placa, ∆T L / 2 = T s ( L / 2) − T ∞ . Como h = q s′′ / ∆T L / 2 es necesario realizar un proceso iterativo para determinar T s ( L/2). Es posible obtener una expresión para la temperatura en cualquier punto en función de la temperatura en el punto x medio: ∆T x = T s ( x) − T ∞ ≈ 1,15 L
1/ 5
∆T L / 2 .
• 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 dada por:
D L
≥
35 Gr L1 / 4
.
18
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• 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 = A s / 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 Ra L1 c/ 4 si
104 ≤ Ra Lc ≤ 107 y Nu L = 0,15 Ra L1 / 3 si 107 ≤ Ra Lc ≤ 1011. Propiedades c
c
calculadas a T f . • 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 Ra L1 c/ 4 con 105 ≤ Ra Lc ≤ 1010. Propiedades calculadas a T f .
• Correlación de Churchill y Chu para la convección libre sobre un cilindro largo horizontal: 2
0,387 Ra D1 / 6 h D Nu D = Ra D ≤ 1012. Propiedades calculadas a T f . = 0,60 + con 8 / 27 k [1 + (0,559 / Pr ) 9 / 16 ] • Convección libre y forzada combinadas. Se produce cuando
correlaciones
convenientes
3 3 3 = Nu forzada ± Nu libre Nu combinada .
corregidas
con
la
Gr L Re L2
≈ 1 . Se utilizan las
siguiente
expresión:
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 / r 2 = sen θ ·dθ ·dφ . Unidad: estereorradián (sr). El ángulo
sólido subtendido por el hemisferio sobre un diferencial de área
∫
d ω =
h
2π
π / 2
∫0 ∫0
senθ d θ d φ = 2π sr .
• Intensidad espectral emitida: I λ ,e (T , λ , θ , φ ) =
W . dA1 cos θ d ω d λ m 2 ·sr· µm
19
dqe
dA1
vale:
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras E λ (T , λ ) = q ′λ ′,e (T , λ ) =
• Potencia emisiva espectral: =∫
2π
0
dq e
=∫
2π
π / 2
∫0
0
dA1 d λ
I λ ,e (T , λ ,θ , φ ) cosθ d ω =
W I λ ,e (T , λ ,θ , φ ) cosθ senθ d θ d φ 2 m · µm
π / 2
∫0
.
∞
′ ,emit = ∫ E λ (T , λ )d λ [W/m2]. • Potencia emisiva total: E (T ) = q ′rad 0
• Emisor difuso: I λ ,e (T , λ ,θ , φ ) = I λ ,e (T , λ ) ⇒ E λ (T , λ ) = π I λ ,e (T , λ ) ⇒ E = π I e . • Irradiación
Gλ (λ ) = q λ ′′,i (λ ) =
espectral:
dq i dA1 ·d λ
=∫
2π
0
π / 2
∫0
I λ ,i (λ ,θ , φ ) cosθ senθ d θ d φ
[W/m2·µ m]. ∞
′ ,inc = ∫ Gλ (λ )d λ [W/m2]. • Irradiación total: G = q ′rad 0
• Radiación incidente difusa: I λ ,i (λ ,θ , φ ) = I λ ,i (λ ) ⇒ Gλ (λ ) = π I λ ,i (λ ) ⇒ G = π I i . • Radiosidad espectral: ′ ,emit + reflj (λ ) = J λ (λ , T ) = q ′rad
dq e + r dA1 ·d λ
=∫
2π
0
π / 2
∫0
I λ ,e + r (λ , T ,θ , φ ) cos θ senθ d θ d φ [W/m2·µ m].
∞
′ ,emit + ref = ∫ J λ (λ , T )d λ [W/m2]. • Radiosidad total: J (T ) = q ′rad 0
• Emisor y reflector difuso: I λ ,e+ r (λ , T ,θ , φ ) = I λ ,e+ r (λ , T ) ⇒ J λ (λ , T ) = π I λ ,e+ r (λ , T ) ⇒ J (T ) = π I e+ r (T ) .
2hc o2 , donde h = 6,6256·10 -34 J·s es la • Distribución de Planck: I λ ,b (λ , T ) = 5 λ [exp( hc o / λ 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 ) =
C 1
λ 5 [exp(C 2 / λ T ) − 1]
, donde C 1 y C 2 son la primera y segunda
constantes de radiación. • Ley de desplazamiento de Wien:
dE λ ,b (λ , T ) d λ
• Ley de Stefan-Boltzmann: E b = ∫
= 0 ⇒ λ máxT = C 3 = 2.898µ m·K . C 1
∞
0
λ [exp(C 2 / λ T ) − 1] 5
W/m2·K 4 es la constante de Stefan-Boltzmann. 20
d λ = σ T 4 ,
donde σ = 5,67·10-8
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras λ
• Emisión de banda: F (0 → λ ) =
λ
∫0 ∫ E , d λ = E , d (λ T ) = f (λ T ) . = 0 ∫0 σ T 5 σ T 4 λ E d ∫0 , E λ ,b d λ
λ
λ b
λ b
∞
λ b
• Emisividad superficial, ε : Emisividad (hemisférica) espectral: ε λ (λ , T ) =
E λ (λ , T ) E λ ,b (λ , T )
.
∞
Emisividad total (hemisférica): ε (T ) =
E (T ) E b (T )
ε (λ , T ) E , (λ , T ) d λ ∫ 0 = . En un emisor difuso λ
λ b
E b (T )
la emisividad direccional es independiente de la dirección. • Absortividad superficial, α : Absortividad (hemisférica) espectral: α λ (λ ) =
Gλ ,abs (λ ) Gλ (λ )
.
∞
Absortividad (hemisférica) total: α =
Gabs
α (λ )G (λ )d λ ∫ 0 = . α = α (λ ,θ , φ , sup .) ≠ α (T ) . ∫0 G (λ )d λ λ
λ
∞
G
λ
∞
• Para radiación solar ( T b = 5.800 K): α Sol
α (λ ) E , (λ , 5.800 K ) d λ ∫ 0 = . ∫0 E , (λ , 5.800 K )d λ λ
λ b
∞
λ b
• Reflectividad superficial, ρ : Reflectividad (hemisférica) espectral: ρ λ (λ ) =
Gλ , ref (λ ) Gλ (λ )
.
∞
Reflectividad (hemisférica) total: ρ =
G ref G
=
∫0 ρ (λ )G (λ )d λ . ∫0 G (λ )d λ λ
λ
∞
λ
• ρ λ + α λ + τ λ = 1 ⇒ ρ + α + τ = 1 . • Ley
de
E 1 (T s )
α 1 ε 1 α 1
=
=
ε 2 α 2
Kirchhoff
E 2 (T s )
α 2
(para
superficies
en
el
interior
de
un
recinto):
= ... = E b (T s ) . Como α ≤ 1 ⇒ E (T s) ≤ E b(Ts). También se cumple:
= ... = 1 ⇒ ε = α .
• Superficie difusa (emisora y receptora difusa): ε λ, θ y α λ ,θ son independientes de la dirección
(θ , φ ) ⇒ ε λ = α λ . independientes de λ . • Superficie gris: ε λ y α λ son
21
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Superficie gris difusa: ε λ, θ y α λ ,θ son independientes de λ (gris) y de la dirección (difusa) ⇒ ε
= α .
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
vista): F ij =
qi → j Ai J i
=
1 Ai
∫∫ Ai
cosθ i cosθ j π R 2
Aj
dAi dA j . F ji =
q j →i A j J j
=
1 A j
∫∫ Ai
cosθ i cosθ j π R 2
Aj
dAi dA j
.
Estas dos ecuaciones son válidas para superficies emisoras y reflectoras difusas y con radiosidad uniforme. • Relación de reciprocidad: Ai F ij = A j F ji . • Regla de la suma en un recinto de N superficies:
N
∑1 F = 1 . j =
ij
• Superficie plana o convexa: F ii = 0. n
• Para una superficie que se puede descomponer en la suma de varias, A j = ∑ Ak , se tiene que: k =1 n
F i ( j ) =
n
∑1 F
ik
k =
∑1 A F k
y F ( j )i =
ki
k =
n
∑1 A
.
k
k =
• Intercambio
neto
de
radiación
entre
dos
superficies
negras:
qij = qi → j − q j →i = Ai F ij E bi − A j F ji E bj = Ai F ijσ (T i 4 − T j4 ) .
• Transferencia neta de radiación desde la superficie i en un recinto con N superficies negras: qi =
N
∑1 A F σ (T 4 − T 4 ) . i
ij
i
j
j =
• Transferencia neta de radiación desde una superficie gris difusa en un recinto: qi = ( J i − Gi ) Ai = ( E i − α i Gi ) Ai =
E bi − J i
(1 − ε i ) / ε i Ai
, siendo (1 - ε i) /ε A i i, la resistencia radiativa
superficial. • Transferencia neta de radiación desde una superficie gris difusa en un recinto: qi =
N
∑1 j =
qij =
N
J i − J j
∑1 ( A F ) 1 , siendo ( A F )-1 la resistencia radiativa geométrica. j =
−
i
i ij
ij
22
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
• Intercambio
neto
q1 = − q 2 = q12 =
de
1 − ε 1 ε 1 A1
radiación
E b1 − E b 2
+
1 A1 F 12
+
1 − ε 2
en
un
recinto
de
dos
superficies:
.
ε 2 A2
• Intercambio neto de radiación entre dos superficies separadas por una cubierta de radiación: A1σ (T 14 − T 24 ) q12 = . 1 1 − ε 3,1 1 − ε 3, 2 1 Para el caso de N cubiertas de radiación con ε iguales + + + ε 1
ε 3,1
ε 3, 2
ε 2
(incluyendo las superficies extremas): (q12 ) N =
1 N + 1
(q12 ) 0 .
• Superficie rerradiante: superficie idealizada en la que la transferencia de calor neta por
radiación es cero: qi = 0 ⇒ J i = Gi = E bi.
23
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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 de perfil rectangular
t
Lc = L + (t / 2)
η f =
1/ 2
m = (2h / kt ) siendo w >> t
w
tanh mLc mLc
L
Aleta recta de perfil triangular
1/ 2
A f = 2 w[ L2 + (t / 2) 2 ] 1/ 2
w
t
m = (2h / kt ) siendo w >> t
L
24
η f =
1 I 1 (2mL) mL I 0 ( 2mL)
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 2.1. Eficiencia de formas comunes de aletas (continuación). Descripción
Aleta recta de perfil parabólico
Esquema
Dimensiones
Eficiencia
A f = w C 1 L2 + ( L2 / t ) ln(t / L + C 1 ) w
t
C 1 =
η f =
[1 + (t / L) 2 ]1 / 2
m = (2h / kt )
1/ 2
2
[4(mL)
2
1/ 2
+ 1]
+1
L
Aleta anular de perfil rectangular
A f = 2π ( r 22c − r 12 ) t
r 2 c = r 2 + (t / 2) m = (2h / kt )1 / 2
L r 1 r 2
25
η f = C 2
K 1 ( mr 1 ) I 1 ( mr 2 c ) − I 1 (mr 1 ) K 1 ( mr 2 c ) I 0 ( mr 1 ) K 1 ( mr 2 c ) + K 0 ( mr 1 ) I 1 ( mr 2 c ) C 2 =
(2r 1 / m) (r 22c − r 12 )
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla 2.1. Eficiencia de formas comunes de aletas (continuación). Descripción
Esquema
Dimensiones
Eficiencia
A f = π DLc
Aleta de aguja cilíndrica
η f =
Lc = L + ( D / 4)
D
m = ( 4h / kD)1 / 2
tanh mLc mLc
L
Aleta de aguja cónica
A f =
D
π D
2
[ L2 + ( D / 2) 2 ]1 / 2
m = ( 4h / kD)
η f =
1/ 2
2 I 2 (2mL) mL I 1 ( 2mL)
L
A f = Aleta de aguja parabólica
y = ( D/2)·(1 - x / L)2 D
π L3 L ln[(2 DC 4 / L) + C 3 ] C 3 C 4 − 8 D 2 D C 3 = 1 + 2( D / L ) 2 C 4 =
L
2 1/ 2
[1 + ( D / L) ]
m = (4h / kD)1 / 2
26
η f =
2
[4 / 9(mL) 2 + 1]1 / 2 + 1
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
TEMA 3. CONDUCCIÓN BIDIMENSIONAL EN RÉGIMEN ESTACIONARIO. Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)]. Descripción del sistema
Esquema
Restricciones
Factor de forma
T 2
1.1. Esfera enterrada en un medio semiinfinito 1.2. Esfera enterrada en un medio semiinfinito con superficie aislada (Bejan pág. 115)
z T 1
S =
2π D 1 − ( D / 4 z )
z > D / 2
S =
2π D 1 + ( D / 4 z )
D
aislado
z T 2
T 1
T 2
D
1.3. Esfera enterrada en un medio infinito (Holman pág. 55) 1.4. Conducción entre dos esferas en un medio infinito (Bejan y Holman)
z > D / 2
T 2
T 2 T 1
T 2
T 2
S = 2π D
Ninguna
D T 2
T 2
T 1
D
w / D > 3
d w
27
S =
2π d d ( D / 2w) 4 d 1 − − D 1 − ( d / 2 w) 2 w
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)] (continuación). Descripción del sistema 1.5. Cavidad hemisférica en medio semiinfinito con superficie aislada (Bejan)
Esquema
Factor de forma
Ninguna
S = π D
aislado
T 1 T 2
T 2 T 2
T 2
2.1. Cilindro de longitud L enterrado en un medio semiinfinito
S =
L >> D
z
2π L arc cosh(2 z / D)
T 1
L >> D z > 3 D / 2
L D
2.2. Cilindro de longitud infinita enterrado en un medio semiinfinito (Rohsenow pág. 3-120) 2.3. Cilindro vertical de longitud L enterrado en un medio semiinfinito 2.4. Conducción entre dos cilindros paralelos de longitud L en un medio infinito
D
aislado
Restricciones
S =
S ′ =
z ≈ D
T 2 z
T 1
S ′ = z > 2 D
2π L ln(4 z / D)
2π arc cosh(2 z / D ) 2π
2 ln (2 z / D ) + (2 z D ) − 1
D
z >> D
S ′ =
2π ln(4 z / D )
L >> D
S =
2π L ln(4 L / D)
T 2
L
T 1
D
T 2
T 1
D
d
L >> D1 , D2 L >> w
w
28
S =
2π L 4w 2 − D 2 − d 2 arc cosh Dd 2
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)] (continuación). Descripción del sistema
Esquema
Restricciones
Factor de forma
T 2
2.5. Cilindro de longitud L en medio de planos paralelos de igual longitud y ancho infinito
z T 1
z >> D / 2 L >> z
L D
S =
2π L ln(8 z / π D)
z
T 2
T 2
2.6. Cilindro de longitud L centrado en un sólido de sección cuadrada de igual longitud
T 1
w > D L >> w
w
S =
2π L ln(1,08w / D )
D
T 2
2.7. Cilindro excéntrico de longitud L en el interior de un cilindro de igual longitud
T 1 d z D
2.8. Fila infinita de cilindros de longitud infinita en un medio semiinfinito (Rohsenow)
29
D > d L >> D
S =
z > D
S ′ =
2π L D 2 + d 2 − 4 z 2 arc cosh Dd 2
2π ln[(2 L / π D)·senh (2π z / L)]
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)] (continuación). Descripción del sistema
Esquema
Restricciones
2.9. Fila infinita de cilindros de longitud infinita en el plano medio de una placa infinita (Rohsenow) 3.1. Cubo enterrado en un medio infinito (Holman)
z > D
L
Factor de forma
S ′ =
2π ln[(2 L / π D )·senh (π z / L)]
S = 8,24 L
Ninguna
4.1. Paralelepípedo inmerso en un medio semiinfinito (Holman)
Ninguna
4.2. Agujero de sección rectangular muy largo en un medio semiinfinito (Rohsenow)
b 1,685 L log1 + a S = 0 , 078 b c
a a > b
S ′ =
2b
3,5 z 0 , 25 0 , 75 a b
ln
30
+ 5,7
−0, 59
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)] (continuación). Descripción del sistema
Esquema
Restricciones
Factor de forma
W
5.1. Pared plana con superficies isotermas (Bejan)
H
A T 1
H > L / 5 W > L / 5
T 2
S =
WH L
L T 2
5.2. Esquina de dos paredes contiguas
L T 1
W > L / 5
S = 0,54W
Ninguna
S = 0,15 L
Ninguna
S = 2 D
T 2 W T 2
L
5.3. Esquina de tres paredes contiguas con diferencia de temperaturas entre las superficies interior y exterior (Bejan)
T 1
T 2 W
T 2
T 2
T 2
T 1
D
6.1. Disco delgado sobre medio semiinfinito
T 2
T 1
6.2. Disco delgado horizontal enterrado en un medio semiinfinito (Kreith pág. 112)
z T 2 D
31
Ninguna
S =
4,45 D 1 − ( D / 5,67 z )
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)] (continuación). Descripción del sistema 6.3. Disco delgado horizontal enterrado en un medio semiinfinito (Bejan) 6.4. Disco delgado horizontal enterrado en un medio semiinfinito con superficie aislada (Bejan)
Esquema
Factor de forma
T 1 z T 2
z > D
S =
2π D (π / 2) − arc tan( D / 4 z )
z > D
S =
2π D (π / 2) + arc tan( D / 4 z )
L / D > 2
S =
2π D (π / 2) − arc tan( D / 2 L)
D
aislado T 1
T 1
z T 2 D
6.5. Dos discos paralelos coaxiales en un medio infinito (Bejan) 7.1. Placa horizontal delgada de anchura W (dimensión al dibujo) enterrada en un medio semiinfinito (Bejan y Holman)
Restricciones
T 1
T 2 D L
T 1 z T 2 L
32
z >> L L > W
S =
z = 0 L > W
S =
L >> W z > 2W
S =
2π L ln(4 L / W ) π L
ln(4 L / W ) 2π L ln(2π z / W )
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 3.1. Factores de forma de conducción para sistemas bidimensionales seleccionados [q = Sk (T 1 - T 2)] (continuación). Descripción del sistema
Esquema
Restricciones
7.2. Placa vertical delgada y larga según la dimensión al dibujo enterrada en un medio semiinfinito (Rohsenow)
7.3. Placa horizontal delgada y larga según la dimensión al dibujo enterrada en un medio semiinfinito (Rohsenow)
Factor de forma
1 z < < 12 2 L
L S ′ = 2,38· z
0 , 24
1 z < < 12 2 L
L S ′ = 2,94· z
0 ,32
T 1 z T 2 L
33
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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.
Bi
0,01 0,02 0,03 0,04 0,05 0,06 0,07 0,08 0,09 0,10 0,15 0,20 0,25 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,0 2,0 3,0 4,0 5,0 6,0 7,0 8,0 9,0 10,0 20,0 30,0 40,0 50,0 100,0 ∞
Pared plana C 1 1 (rad)
0,0998 0,1410 0,1732 0,1987 0,2217 0,2425 0,2615 0,2791 0,2956 0,3111 0,3779 0,4328 0,4801 0,5218 0,5932 0,6533 0,7051 0,7506 0,7910 0,8274 0,8603 1,0769 1,1925 1,2646 1,3138 1,3496 1,3766 1,3978 1,4149 1,4289 1,4961 1,5202 1,5325 1,5400 1,5552 1,5707
1,0017 1,0033 1,0049 1,0066 1,0082 1,0098 1,0114 1,0130 1,0145 1,0160 1,0237 1,0311 1,0382 1,0450 1,0580 1,0701 1,0814 1,0919 1,1016 1,1107 1,1191 1,1795 1,2102 1,2287 1,2402 1,2479 1,2532 1,2570 1,2598 1,2620 1,2699 1,2717 1,2723 1,2727 1,2731 1,2733
Cilindro infinito C 1 1 (rad)
0,1412 0,1995 0,2439 0,2814 0,3142 0,3438 0,3708 0,3960 0,4195 0,4417 0,5376 0,6170 0,6856 0,7465 0,8516 0,9408 1,0185 1,0873 1,1490 1,2048 1,2558 1,5995 1,7887 1,9081 1,9898 2,0490 2,0937 2,1286 2,1566 2,1795 2,2881 2,3261 2,3455 2,3572 2,3809 2,4050
1,0025 1,0050 1,0075 1,0099 1,0124 1,0148 1,0173 1,0197 1,0222 1,0246 1,0365 1,0483 1,0598 1,0712 1,0932 1,1143 1,1346 1,1539 1,1725 1,1902 1,2071 1,3384 1,4191 1,4698 1,5029 1,5253 1,5411 1,5526 1,5611 1,5677 1,5919 1,5973 1,5993 1,6002 1,6015 1,6018
Esfera C 1 1 (rad)
0,1730 0,2445 0,2989 0,3450 0,3852 0,4217 0,4550 0,4860 0,5150 0,5423 0,6608 0,7593 0,8448 0,9208 1,0528 1,1656 1,2644 1,3225 1,4320 1,5044 1,5708 2,0288 2,2889 2,4556 2,5704 2,6537 2,7165 2,7654 2,8044 2,8363 2,9857 3,0372 3,0632 3,0788 3,1102 3,1415
1,0030 1,0060 1,0090 1,0120 1,0149 1,0179 1,0209 1,0239 1,0268 1,0298 1,0445 1,0592 1,0737 1,0880 1,1064 1,1441 1,1713 1,1978 1,2236 1,2488 1,2732 1,4793 1,6227 1,7201 1,7870 1,8338 1,8674 1,8921 1,9106 1,9249 1,9781 1,9898 1,9942 1,9962 1,9990 2,0000
Bi = hL/k para la pared plana y Bi = hr o/k para el cilindro infinito y la esfera.
34
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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 semiinfinito
S ( x, t ) =
x
T ( x, t ) − T ∞
P ( x, t )
Pared plana
P ( x, t ) =
x
T ini − T ∞
sólido semiinfinito
T ( x, t ) − T ∞ T ini − T ∞
pared plana
2 L1
C (r, t )
Cilindro infinito
C ( r , t ) =
r
T (r , t ) − T ∞ T ini − T ∞
r o
S ( x 1 , t )· P ( x 2, t )
Placa semiinfinita
S ( x1 , t )· P ( x 2 , t )
x 2 x 1
2 L2
35
cilindro infinito
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla 4.2. Soluciones de conducción transitoria de sistemas multidimensionales expresadas como producto de soluciones de sistemas unidimensionales (continuación). Sistema
Barra rectangular infinita
Esquema
Solución
P ( x 1 , t ) ·P ( x 2, t )
P ( x1 , t )· P ( x 2 , t )
x 2
x 1
2 L1 2 L2
C (r, t )· S ( x , t )
Cilindro semiinfinito
C ( r , t )·S ( x, t )
r x r o
Barra rectangular semiinfinita
P ( x 1 , t )· P ( x 2, t )· S ( x 3, t )
S ( x 3 , t )· P ( x1 , t )· P ( x 2 , t )
x 3 x 2
x 1
2 L1
2 L2
P ( x 1 , t ) ·P ( x 2, t )· P ( x 3, t ) x 3
Paralelepípedo rectangular
2 L3
x 2
x 1
P ( x1 , t )· P ( x 2 , t )· P ( x3 , t )
2 L1 2 L2 C (r, t )· P ( x , t )
Cilindro corto
x
2 L1
C ( r , t )· P ( x, t )
r r o
36
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
TEMA 6. CONVECCIÓN FORZADA EN FLUJO EXTERNO Tabla 6.1. Tabla resumen de correlaciones para flujo externo sobre placa plana. Correlaciones
1*
Nu x =
h x x
2**
Nu x =
h x x
3*
Nu x =
4**
h L L
Nu L =
5**
k
Nu L =
dq = h x (T s − T ∞ )dA s q ′′ = h (T s − T ∞ )
= 0,664 Re x1 / 2 Pr 1 / 3
k k
q x′′ = h x (T s − T ∞ )
= 0,332 Re x1 / 2 Pr 1 / 3
k
h x x
Transferencia de calor
q = h A s (T s − T ∞ ) q x′′ = h x (T s − T ∞ )
= 0,0296 Re x4 / 5 Pr 1 / 3
dq = h x (T s − T ∞ )dA s q ′′ = h (T s − T ∞ )
= (0,037 Re L4 / 5 − 871) Pr 1 / 3 h L L
q = h A s (T s − T ∞ ) q ′′ = h (T s − T ∞ )
= 0,037 Re L4 / 5 Pr 1 / 3
k
q = h A s (T s − T ∞ )
q s′′ = cte. / Nu x = 0,453 Re x Pr 1/ 2
6*
1/ 3
T s ( x) = T ∞ +
q s′′ h x ( x)
q s′′ = cte. / Nu x = 0,0308 Re x4 / 5 Pr 1 / 3
7* *: Re = x
ρ u ∞ x µ
=
u ∞ x
ν
**: Re = L
T s ( x) = T ∞ +
ρ u ∞ L µ
=
u ∞ L
ν
q s′′ h x ( x)
Condiciones
Placa a temperatura T s constante. Régimen laminar. Valor local en x. Pr > 0,6. Placa a T s constante. Régimen laminar. Valor promedio entre 0 y x (ó entre 0 y x = L). Pr > 0,6. Placa a T s constante. Régimen turbulento. Valor local en x. 0,6 < Pr < 60. Placa a T s constante. Régimen mixto (parte laminar y parte turbulento). Valor promedio entre 0 y x = L. Pr > 0,6. 5·105 < Re x,c < 108. Placa a T s constante. Régimen predominantemente turbulento (parte laminar despreciable ⇒ L >> xc y Re L >> Re x,c). Valor promedio entre 0 y x = L. Pr > 0,6. 5·105 < Re x,c < 108. Placa que desprende un flujo de calor uniforme. Régimen laminar. Valor local en x. Pr > 0,6. Placa que desprende un flujo de calor uniforme. Rég. turbulento. Valor local en x. 0,6 < Pr < 60.
Condición de rég. turbulento para placa plana: Re x,c > 5·105
Número de Prandtl: Pr =
ν α
=
µ ρ c p ρ k
=
µ c p k
En todas las correlaciones las propiedades del fluido se calculan a la temperatura de película: T = (T + T ∞ ) / 2 ; T s: Temperatura de la superficie [K]; T ∞: Temp. del flujo libre [K]; A s: Á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]. f
37
s
Tabla 6.2. Tabla resumen de correlaciones para flujo cruzado sobre cilindros. Correlaciones para flujo cruzado sobre un cilindro Nu D =
1
Nu D =
2
h D
Transferencia de calor q ′′ = h (T s − T ∞ )
= CRe Dm Pr 1 / 3
k
Pr = CRe Dm Pr n k Pr s
h D
q = h A s (T s − T ∞ )
1/ 4
q ′′ = h (T s − T ∞ ) q = h A s (T s − T ∞ )
Condiciones
Correlación de Hilpert. Los valores de las constantes C y m se dan en la Tabla 6.3 en función de Re D. La Tabla 6.4 da los valores de las constantes para cilindros no circulares. Las propiedades se evalúan a T f . Válida para fluidos con Pr ≥ 0,7. n = 0,37 si Pr ≤ 10 Correlación de Zhukauskas. Con y n = 0,36 si Pr > 10 0,7 < Pr < 500 1 < Re < 10 6 . Los valores de las constantes C y m se dan en la Tabla D
6.5 en función de Re D. Las propiedades se evalúan a T ∞, excepto Pr s a T s 3
Nu D =
h D k
= 0,3 +
0,62 Re D1 / 2 Pr 1 / 3 2 / 3 1/ 4
[1 + (0,4 / Pr ) ]
Re D 1 + 282000
5/8
4/5
q ′′ = h (T s − T ∞ )
Correlaciones para flujo cruzado sobre un banco de N cilindros
q = h A s (T s − T ∞ )
Correlación de Churchill y Bernstein. Con propiedades a T f y Re D· Pr > 0,2.
Transferencia de calor
Condiciones
Correlación de Zhukauskas. Con Nu D =
4
Pr = CRe Dm, máx Pr 0, 36 k Pr s
h D
q ′′ = h ∆T ml
1/ 4
q′ = h π D∆T ml
′ = h N π DL∆T ml q = q NL
N L ≥ 20 1.000 < Re 6 2 · 10 < , D máx . 0,7 < Pr < 500
Las
constantes C y m se dan en la Tabla 6.6. Las propiedades se evalúan a T = (T ent + T sal ) / 2 , excepto Pr s a T s. Para N L < 20 se aplica un factor de corrección tal que Nu D
N L < 20
= C 2 Nu D
N L ≥ 20
, donde C 2 está dado en la
Tabla 6.7 Re D =
ρ u ∞ D µ
=
u ∞ D
ν
Re D , max =
ρ V max D µ
=
V max D
Diferencia de temperaturas media logarítmica: ∆T ml =
Config. alineada: V máx =
S T S T − D
2( S D − D) > ( S T − D) ó V máx =
V ;
Cálculo de la temperatura de salida del flujo:
ν
(T s − T ent ) − (T s − T sal ) T − T ln s ent T s − T sal
Config. escalonada: V máx =
S T
2( S D − D)
V
si 2( S
D
S T S T − D
T s − T sal T s − T ent
π DN h = exp − VN S c ρ T T p
S T : espaciado transversal; S L: espaciado longitudinal; N T: número de tubos en direc. transversal; N L: número de tubos en direc. longitudinal; N = N T x N L: núm. total de tubos.
V
si
− D) < ( S T − D) ;
38
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 6.3. Coeficientes de la correlación de Hilpert para flujo cruzado sobre un cilindro (Pr 0,7). Re D
0,4 - 4 4 - 40 40 - 4.000 4.000 - 40.000 40.000 - 400.000
C
m
0,989 0,911 0,683 0,193 0,027
0,330 0,385 0,466 0,618 0,805
Tabla 6.4. Coeficientes de la correlación de Hilpert para flujo de aire cruzado sobre un paralelepípedo. Geometría Cuadrado en diagonal
Cuadrado recto
Hexágono recto
Hexágono en diagonal
Dibujo V
V
V
V
Re D
C
m
D
5·103 - 105
0,246
0,588
D
5·103 - 105
0,102
0,675
5·103 - 1,95·104
0,160
0,638
1,95·104 - 105
0,0385
0,782
5·103 - 105
0,153
0,638
4·103 - 1,5·104
0,228
0,731
D
D
D
Placa vertical
Tabla 6.5. Coeficientes de la correlación de Zhukauskas para flujo de aire cruzado sobre un cilindro. Re D
C
m
1 - 40 40 - 1.000 103 - 2·105 2·105 - 106
0,75 0,51 0,26 0,076
0,4 0,5 0,6 0,7
39
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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 Alineado Escalonado Alineado Escalonado Alineado ( S T / S L > 0,7) Escalonado ( S T / S L < 2) Escalonado ( S T / S L > 2) Alineado Escalonado
Re D,máx
C
2
m
10 - 10 0,80 0,40 2 10 - 10 0,90 0,40 2 3 10 - 10 Se aproxima como un cilindro único 2 3 10 - 10 Se aproxima como un cilindro único 3 5 10 - 2·10 0,27 0,63 3 5 1/5 10 - 2·10 0,35(S T / S L) 0,60 3 5 10 - 2·10 0,40 0,60 5 6 2·10 - 2·10 0,021 0,84 5 6 2·10 - 2·10 0,022 0,84
Para S T / S L < 0,7, la transferencia de calor es ineficiente y los tubos alineados no se deben usar.
40
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla 6.7. Coeficiente de corrección C 2 de la correlación de Zhukauskas para el flujo 3
cruzado sobre un banco de tubos para N L < 20 y Re D > 10 . N L Alineado Escalonado
1 0,70 0,64
2 0,80 0,76
3 0,86 0,84
4 0,90 0,89
5 0,92 0,92
7 0,95 0,95
10 0,97 0,97
13 0,98 0,98
16 0,99 0,99
Figura 6.2. Disposición de los tubos en configuración alineada (a) y escalonada (b) en un banco de tubos. (Incropera)
N L N T
Figura 6.3. Esquema de un banco de tubos en flujo cruzado. (Incropera)
41
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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)
42
Transferencia de Calor / Curso 2009-10 Fórmulas, Tablas y Figuras Diagrama 7.1. Metodología para seleccionar las correlaciones de convención forzada en flujo interno.
- Región de entrada ⇒ ( x < xcd ,t ó G z -1 < 0,05) ⇒ Figuras 7.1. y 7.2. Siendo: xcd,t = 0.05· D· Re D· Pr . - Correlación local: - Región c. d. ⇒ Nu D = cte ( x > xcd ,t ó G z -1 > 0,05): Re D < 2.300 ⇒ Régimen Laminar:
- q s′′ = cte ⇒ Nu D = 4,36. - T s = cte ⇒ Nu D = 3,66. - Tubo no circular ⇒ Tabla 7.1. - 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 :
- T s = 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): Re D > 2.300 ⇒ Régimen Turbulento:
- Correlación de Sieder y Tate: (se usa con aceite siempre; con agua y con aire si hay grandes ∆T ). - Correl. de Dittus-Boelter. con prop. a T m =
- Correlación promedio (condiciones c. d.: L/ D > 60): * En el
- Correl. de Sieder y Tate.
temario de este curso no se estudia la región de entrada en régimen turbulento.
43
T m,ent + T m, sal
2
.
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 7.2. Números de Nusselt para flujo laminar completamente desarrollado en tubos de diferente sección transversal. Nu D Sección transversal Circular Rectangular (a = altura, b =base) Rectangular (a = altura, b =base) Rectangular (a = altura, b =base) Rectangular (a = altura, b =base) Rectangular (a = altura, b =base) Rectangular (a = altura, b =base) Rectangular (a = altura, b =base) Triangular
b
hD h k
q s uniforme
T s uniforme
-
4,36
3,66
1,0
3,61
2,98
1,43
3,73
3,08
2,0
4,12
3,39
3,0
4,79
3,96
4,0
5,33
4,44
8,0
6,49
5,60
∞
8,23
7,54
-
3,11
2,47
a
45
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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: T s > T a la izquierda y T s < T a la derecha.
46
Tabla 8.1. Tabla resumen de correlaciones de convención libre. Transferencia de calor
Correlación 1/ 4
1
0,75 Pr 1 / 2 Gr Nu x = = x · 1/ 4 k 4 (0,609 + 1,221 Pr 1 / 2 + 1,238 Pr )
2
0,387 Ra L1 / 6 Nu L = = 0,825 + k [1 + (0,492 / Pr ) 9 / 16 ]8 / 27
hx
2
h L
Nu L =
3
h L k
= 0,68 +
Nu Lc =
4
h Lc k
4
5
h Lc k
dq = h(T s − T m )dA z
T s > T ∞
q ′′ = h (T s − T ∞ )
T s < T ∞
Condiciones
Placa vertical con temperatura superficial constante, T s = cte. Régimen laminar, Gr x < 109. Correlación local. h L 4 Correlación promedio: Nu L = = Nu L k 3 Correlación de Churchill y Chu. Placa vertical con T s = cte. Correlación promedio. Válida para todo Ra L.
q ′′ = h (T s − T ∞ )
Placa vertical con T s = cte. Ra L ≤ 109. Correlación promedio.
q = h A s (T s − T ∞ ) T s > T ∞
T s < T ∞
T s > T ∞
T s < T ∞
q ′′ = h (T s − T ∞ )
7
con 10 ≤ Ra Lc ≤ 10 Nu Lc =
q x′′ = h(T s − T m )
q = h A s (T s − T ∞ )
0,670 Ra L1 / 4 [1 + (0,492 / Pr )9 / 16 ]4 / 9
= 0,54 Ra L1c/ 4
Representación gráfica
q = h A s (T s − T ∞ )
= 0,15 Ra L1 c/ 3
con 107 ≤ Ra Lc ≤ 1011
Nu Lc =
6
h Lc 5
con 10
7
Gr x =
Nu D =
= 0,27 Ra Lc k ≤ Ra Lc ≤ 1010
q = h A s (T s − T ∞ )
0,387 Ra D1 / 6 = 0,60 + 9 / 16 8 / 27 k [ 1 ( 0 , 559 / ) ] Pr +
h D
g β (T s − T ∞ ) x 3
ν
q ′′ = h (T s − T ∞ )
1/ 4
2
=
Ra x Pr
; Ra L = Gr L Pr =
2
q ′′ = h (T s − T ∞ )
g β (T s − T ∞ ) L3
να
T s > T ∞
q = h A s (T s − T ∞ )
; Ra D =
g β (T s − T ∞ ) D 3
να
Placa horizontal con T s = cte. Superficie superior de placa caliente o inferior de placa fría. Correlación promedio. Longitud característica definida como el cociente entre el área y el perímetro de la placa: Lc = A s / P . Placa horizontal con T s = cte. Superficie inferior de placa caliente o superior de placa fría. Correlación promedio. Longitud característica definida como el cociente entre el área y el perímetro de la placa: Lc = A s / P .
Correlación de Churchill y Chu (promedio) para la convección libre sobre un cilindro largo horizontal: con Ra D ≤ 1012.
; Propiedades calculadas a T f = (T s + T ∞)/2; Correlaciones 1 a 3: válidas para q x′′ = cte si
propiedades, Nu L y Ra L se definen en función de la temperatura en el punto medio de la placa: T f = (T s( L/2) + T ∞)/2; ∆T L / 2 = T s ( L / 2) − T ∞ ⇒ h = q s′′ / ∆T L / 2 ⇒ 1/5
∆T x = T s ( x) − T ∞ ≈ 1,15( x / L ) ∆T L / 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
⇒ ( D / L ) ≥ (35 / Gr L ) ; 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). 1/ 4
47
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
TEMA 9. INTRODUCCIÓN A LA RADIACIÓN Figura 9.1. Distribución de Planck. Potencia Emisiva espectral del cuerpo negro.
10 10 10 10 10 10
Potencia emisiva espectral, E b 2 (W/m · m)
10 10 10 10 10 10 10 10
9 8
T = 5.800 K
λ máx·T = 2.898 m·K
7 6 5
T = 2.000 K
4
T = 1.000 K
3 2
T = 800 K
1
T = 300 K
0
T = 100 K
-1 -2
T = 50 K
-3 -4
0.1
0.2
0.4 0.6
1
2
4
6
10
20
40 60
longitud de onda, ( m) Tabla 9.1. Funciones de radiación de cuerpo negro.
T ( m·K)
200 400 600 800 1.000 1.200 1.400 1.600 1.800 2.000 2.200 2.400 2.600 2.800 2.898 3.000 3.200 3.400
I ,b ( , T ) F ( 0
5
T -1 ( m·K·sr)
)
0,3750·10-27 0,4903·10-13 0,1040·10-8 0,9911·10-7 0,1185·10 -5 0,5239·10 -5 0,1344·10 -4 0,2491·10 -4 0,3756·10 -4 0,4934·10 -4 0,5896·10 -4 0,6589·10 -4 0,7013·10-4 0,7202·10-4 0,7223·10 -4 0,7203·10-4 0,7060·10-4 0,6815·10-4
0,000000 0,000000 0,000000 0,000016 0,000321 0,002134 0,007790 0,019718 0,03934 0,06673 0,10089 0,14026 0,1831 0,2279 0,250108 0,2732 0,3181 0,3617 48
I ,b ( , T ) I ,b (
máx
, T )
0,000000 0,000000 0,000014 0,001372 0,016406 0,072534 0,18608 0,3449 0,5199 0,6831 0,8163 0,9122 0,9709 0,9971 1,000000 0,9971 0,9774 0,9436
100
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla 9.1. Funciones de radiación de cuerpo negro (continuación). T ( m·K)
3.600 3.800 4.000 4.200 4.400 4.600 4.800 5.000 5.200 5.400 5.600 5.800 6.000 6.200 6.400 6.600 6.800 7.000 7.200 7.400 7.600 7.800 8.000 8.500 9.000 9.500 10.000 10.500 11.000 11.500 12.000 13.000 14.000 15.000 16.000 18.000 20.000 25.000 30.000 40.000 50.000 75.000 100.000
F ( 0
I ,b ( , T )
)
T 5 -1 ( m·K·sr)
0,6504·10-4 0,6152·10-4 0,5781·10-4 0,5404·10-4 0,5033·10-4 0,4673·10-4 0,4331·10-4 0,4008·10-4 0,3706·10-4 0,3424·10-4 0,3164·10-4 0,2923·10-4 0,2701·10-4 0,2497·10-4 0,2310·10-4 0,2138·10-4 0,1980·10-4 0,1835·10-4 0,1703·10-4 0,1581·10-4 0,1469·10-4 0,1366·10-4 0,1272·10-4 0,1068·10-4 0,9015·10-5 0,7653·10-5 0,6533·10-5 0,5605·10-5 0,4833·10-5 0,4187·10-5 0,3644·10-5 0,2795·10-5 0,2176·10-5 0,1719·10-5 0,1374·10-5 0,9082·10-6 0,6233·10-6 0,2765·10-6 0,1405·10-6 0,4739·10-7 0,2016·10-7 0,4186·10-8 0,1358·10-8
0,4036 0,4434 0,4809 0,5160 0,5488 0,5793 0,6076 0,6337 0,6590 0,6804 0,7010 0,7202 0,7378 0,7541 0,7692 0,7832 0,7961 0,8081 0,8192 0,8295 0,8391 0,8480 0,8563 0,8746 0,8900 0,9031 0,9142 0,9237 0,9319 0,9400 0,9451 0,9551 0,9629 0,9700 0,9738 0,9809 0,9856 0,9922 0,9953 0,9980 0,9990 0,9997 0,999905
49
I ,b ( , T ) I ,b (
máx
, T )
0,9004 0,8517 0,8003 0,7481 0,6967 0,6470 0,5996 0,5549 0,5130 0,4741 0,4380 0,4047 0,3740 0,3457 0,3198 0,2960 0,2741 0,2541 0,2357 0,2188 0,2034 0,1891 0,1761 0,1478 0,1248 0,1060 0,09044 0,077600 0,066913 0,057970 0,050448 0,038689 0,030131 0,023794 0,019026 0,012574 0,008629 0,003828 0,001945 0,000656 0,000279 0,000058 0,000019
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
TEMA 10. INTERCAMBIO RADIATIVO ENTRE SUPERFICIES Tabla 10.1. Factores de forma radiativos para geometrías bidimensionales.
Geometría
Esquema
Factor de forma [(W i + W j ) 2 + 4]1 / 2 F ij = − 2W i
wi i
Placas paralelas centradas
L
[(W − W )
−
j
+ 4]
1/ 2
2W i
j
W i = wi / L W j = w j / L
w j
Placas paralelas inclinadas de igual anchura y una arista en común
2
i
j
w
α 2
F ij = 1 − sen i w j
Placas perpendiculares con una arista en común
1/ 2
1 + ( w j / wi ) − [1 + ( w j / wi ) 2 ] F ij = 2
w j
i wi
wk
Recinto de tres lados
w j k
F ij =
j
wi + w j − wk
2 wi
i wi
F ij = j i r i
Cilindros paralelos de radios diferentes
1 2 2 1/ 2 π + [C − ( R + 1) ] − 2π
− [C 2 − ( R − 1) 2 ]
1/ 2
r j
+
R 1 cos −1 − − 180º C C π R 1 cos −1 + − ( R + 1) 180º C C R = r j / r i S = s / r i C = 1 + R + S + ( R − 1)
s
π
j r
Cilindro y placa paralelos
L s1
i
F ij =
r s1 − s 2
π
s1
180º
L
−1 tan
− tan −1
s 2
s2
s
j D
Plano infinito y fila de cilindros
D 2 F ij = 1 − 1 − s
1/ 2
+
2 2 D π −1 s − D tan + 2 D s 180º
i
50
1/ 2
L
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Figura 10.1. Factor de forma radiativo para dos rectángulos paralelos alineados.
j
L
i
Y
X 1 0.8
Y / L = inf.
Y / L = 10 Y/L=4
0.5
Y/L=2 Y/L=1
0.3
Y / L = 0,6
0.2
Y / L = 0,4
j i F
0.1 0.08
Y / L = 0,2
0.05 Y / L = 0,1
0.03 0.02
0.01 0.1
Y / L = 0,05
0.2 0.3
0.5
1
2 X / L
51
3
5
10
20
40
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Figura 10.2. Factor de forma radiativo para dos discos paralelos coaxiales.
j
r j L
i
r i
1 0.9
8 6
5
3
4
0.8
rj / L = 2
0.7 1,5
0.6
1,25
j i 0.5 F
1
0.4 0,8
0.3 0.2
0,6
0.1
0,4 0,3
0 0.1
0.2
0.4
0.6 0.8 1 L / ri
52
2
4
6
8 10
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Figura 10.3. Factor de forma radiativo para dos rectángulos perpendiculares con una arista en común.
j Z
i
X
Y 0.5 Y / X = 0,02
0.45 0,05
0,1
0.4
0,2
0.35 0,4
0.3 0,6
j i 0.25 F
1
0.2
1,5
2
0.15 4
0.1
10
0.05 0 0.1
20
0.2
0.4
0.6 0.8 1 Z / X
53
2
4
6
8 10
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla 10.2. Intercambio neto de radiación en recintos especiales de dos superficies grises y difusas.
Geometría
Esquema
Planos paralelos grandes (infinitos)
plano 1 A2 , T 2 ,
Condiciones
A1 , T 1 ,
2
A1 = A2 = A 1
Objeto convexo pequeño en una cavidad grande
1 ε 1
A1 A2
=
r 1
1 ε 2
−1
r 2
q12
A1σ (T 14 − T 24 ) = 1 1 − ε 2 r 1
F 12 = 1
A1
r 1
A2 r 2
A1 , T 1 ,
=
r 12 r 22
F 12 = 1
A1
1
A2 A2 , T 2 ,
+
r 2
Cilindros concéntricos largos (infinitos)
Esferas concéntricas
q12 =
F 12 = 1
plano 2 r 1
Intercambio de radiación Aσ (T 14 − T 24 )
≈0
F 12 = 1 2
54
ε 1
q12 =
+
ε 2
r 2
A1σ (T 14 − T 24 )
1 ε 1
+
1 − ε 2 r 1 ε 2
2
r 2
q12 = A1 ε 1σ (T 14 − T 24 )
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
TABLAS DE PROPIEDADES TERMOFÍSICAS Y DE FUNCIONES MATEMÁTICAS Tabla A. Propiedades termofísicas del aire a presión atmosférica. T (K)
100 150 200 250 300 350 400 450 500 550 600 650 700 750 800 850 900 950 1000 1100 1200 1300 1400 1500
c p 3 (kg/m ) (J/kg·K)
3,5562 2,3364 1,7548 1,3947 1,1614 0,9950 0,8711 0,7740 0,6964 0,6329 0,5804 0,5356 0,4975 0,4643 0,4354 0,4097 0,3868 0,3666 0,3482 0,3166 0,2902 0,2679 0,2488 0,2322
1032 1012 1007 1006 1007 1009 1014 1021 1030 1040 1051 1063 1075 1087 1099 1110 1121 1131 1141 1159 1175 1189 1207 1230
3
7 6 k ·10 ·10 ·10 2 2 (N·s/m ) (m /s) (W/m·K)
71,1 103,4 132,5 159,6 184,6 208,2 230,1 250,7 270,1 288,4 305,8 322,5 338,8 354,6 369,8 384,3 398,1 411,3 424,4 449,0 473,0 496,0 530 557
2,00 4,426 7,590 11,44 15,89 20,92 26,41 32,39 38,79 45,57 52,69 60,21 68,10 76,37 84,93 93,80 102,9 112,2 121,9 141,8 162,9 185,1 213 240
55
9,34 13,8 18,1 22,3 26,3 30,0 33,8 37,3 40,7 43,9 46,9 49,7 52,4 54,9 57,3 59,6 62,0 64,3 66,7 71,5 76,3 82 91 100
6
·10 2 (m /s)
Pr
2,54 5,84 10,3 15,9 22,5 29,9 38,3 47,2 56,7 66,7 76,9 87,3 98,0 109 120 131 143 155 168 195 224 238 303 350
0,786 0,758 0,737 0,720 0,707 0,700 0,690 0,686 0,684 0,683 0,685 0,690 0,695 0,702 0,709 0,716 0,720 0,723 0,726 0,728 0,728 0,719 0,703 0,685
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla B. Propiedades termofísicas del aceite de motor a presión atmosférica. T (K)
273 280 290 300 310 320 330 340 350 360 370 380 390 400 410 420 430
c p 3 (kg/m ) (J/kg·K)
899,1 895,3 890,0 884,1 877,9 871,8 865,8 859,9 853,9 847,8 841,8 836,0 830,6 825,1 818,9 812,1 806,5
1796 1827 1868 1909 1951 1993 2035 2076 2118 2161 2206 2250 2294 2337 2381 2427 2471
3
k ·10 ·102 ·106 2 2 (N·s/m ) (m /s) (W/m·K)
385 217 99,9 48,6 25,3 14,1 8,36 5,31 3,56 2,52 1,86 1,41 1,10 0,874 0,698 0,564 0,470
4280 2430 1120 550 288 161 96,6 61,7 41,7 29,7 22,0 16,9 13,3 10,6 8,52 6,94 5,83
56
147 144 145 145 145 143 141 139 138 138 137 136 135 134 133 133 132
·107 2 (m /s)
Pr
·103 -1 (K )
0,910 47000 0,70 0,880 27500 0,70 0,872 12900 0,70 0,859 6400 0,70 0,847 3400 0,70 0,823 1965 0,70 0,800 1205 0,70 0,779 793 0,70 0,763 546 0,70 0,753 395 0,70 0,738 300 0,70 0,723 233 0,70 0,709 187 0,70 0,695 152 0,70 0,682 125 0,70 0,675 103 0,70 0,662 88 0,70
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla C. Propiedades termofísicas del agua saturada. T (K)
P (bar)
273,15 275 280 285 290 295 300 305 310 315 320 325 330 335 340 345 350 355 360 365 370 373,15 375 380 385 390 400 410 420 430 440 450 460 470 480 490 500
0,00611 0,00697 0,00990 0,01387 0,01917 0,02617 0,03531 0,04712 0,06221 0,08132 0,1053 0,1351 0,1719 0,2167 0,2713 0,3372 0,4163 0,5100 0,6209 0,7514 0,9040 1,0133 1,0815 1,2869 1,5233 1,794 2,455 3,302 4,370 5,699 7,333 9,319 11,71 14,55 17,90 21,83 26,40
i fg c p 3 (kg/m ) (kJ/kg) (J/kg·K)
1000 1000 1000 1000 999,0 998,0 997,0 995,0 993,0 991,1 989,1 987,2 984,3 982,3 979,4 976, 973,7 970,9 967,1 963,4 960,6 957,9 956,9 953,3 949,7 945,2 937,2 928,5 919,1 909,9 900,9 890,5 879,5 868,1 856,9 844,6 831,3
2502 2497 2485 2473 2461 2449 2438 2426 2414 2402 2390 2378 2366 2354 2342 2329 2317 2304 2291 2278 2265 2257 2252 2239 2225 2212 2183 2153 2123 2091 2059 2024 1989 1951 1912 1870 1825
4217 4211 4198 4189 4184 4181 4179 4178 4178 4179 4180 4182 4184 4186 4188 4191 4195 4199 4203 4209 4214 4217 4220 4226 4232 4239 4256 4278 4302 4331 4360 4400 4440 4480 4530 4590 4660
3
k ·10 ·106 2 (N·s/m ) (W/m·K)
1750 1652 1422 1225 1080 959 855 769 695 631 577 528 489 453 420 389 365 343 324 306 289 279 274 260 248 237 217 200 185 173 162 152 143 136 129 124 118
569 574 582 590 598 606 613 620 628 634 640 645 650 656 660 668 668 671 674 677 679 680 681 683 685 686 688 688 688 685 682 678 673 667 660 651 642
i fg : entalpía específica del cambio de fase entre líquido y gas.
57
Pr
·106 -1 (K )
12,99 12,22 10,26 8,81 7,56 6,62 5,83 5,20 4,62 4,16 3,77 3,42 3,15 2,88 2,66 2,45 2,29 2,14 2,02 1,91 1,80 1,76 1,70 1,61 1,53 1,47 1,34 1,24 1,16 1,09 1,04 0,99 0,95 0,92 0,89 0,87 0,86
-68,05 -32,74 46,04 114,1 174,0 227,5 276,1 320,6 361,9 400,4 436,7 471,2 504,0 535,5 566,0 595,4 624,2 652,3 697,9 707,1 728,7 750,1 761 788 814 841 896 852 1010
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla D. Funciones hiperbólicas. x
senh x
cosh x
tanh x
x
senh x
cosh x
tanh x
0,00 0,10 0,20 0,30 0,40 0,50 0,60 0,70 0,80 0,90 1,00 1,10 1,20 1,30 1,40 1,50 1,60 1,70 1,80 1,90
0,0000 0,1002 0,2013 0,3045 0,4108 0,5211 0,6367 0,7586 0,8881 1,0265 1,1752 1,3356 1,5095 1,6984 1,9043 2,1293 2,3756 2,6456 2,9422 3,2682
1,0000 1,0050 1,0201 1,0453 1,0811 1,1276 1,1855 1,2552 1,3374 1,4331 1,5431 1,6685 1,8107 1,9709 2,1509 2,3524 2,5775 2,8283 3,1075 3,4177
0,00000 0,09967 0,19738 0,29131 0,37995 0,46212 0,53705 0,60437 0,66404 0,71630 0,76159 0,80050 0,83365 0,86172 0,88535 0,90515 0,92167 0,93541 0,94681 0,95624
2,00 2,10 2,20 2,30 2,40 2,50 2,60 2,70 2,80 2,90 3,00 3,50 4,00 4,50 5,00 6,00 7,00 8,00 9,00 10,00
3,6269 4,0219 4,4571 4,9370 5,4662 6,0502 6,6947 7,4063 8,1919 9,0596 10,018 16,543 27,290 45,003 74,203 201,71 548,32 1490,5 4051,5 11013
3,7622 4,1443 4,5679 5,0372 5,5569 6,1323 6,7690 7,4735 8,2527 9,1146 10,068 16,573 27,308 45,014 74,210 201,72 548,32 1490,5 4051,5 11013
0,96403 0,97045 0,97574 0,98010 0,98367 0,98661 0,98903 0,99101 0,99263 0,99396 0,99505 0,99818 0,99933 0,99975 0,99991 0,99999 1,0000 1,0000 1,0000 1,0000
Tabla E. Función gaussiana de error. x
erf ( x )
x
erf ( x )
x
erf ( x )
0,00 0,02 0,04 0,06 0,08 0,10 0,12 0,14 0,16 0,18 0,20 0,22 0,24 0,26 0,28 0,30 0,32 0,34
0,00000 0,02256 0,04511 0,06762 0,09008 0,11246 0,13476 0,15695 0,17901 0,20094 0,22270 0,24430 0,26570 0,28690 0,30788 0,32863 0,34913 0,36936
0,36 0,38 0,40 0,44 0,48 0,52 0,56 0,60 0,64 0,68 0,72 0,76 0,80 0,84 0,88 0,92 0,96 1,00
0,38933 0,40901 0,42839 0,46623 0,50275 0,53790 0,57162 0,60386 0,63459 0,66378 0,69143 0,71754 0,74210 0,76514 0,78669 0,80677 0,82542 0,84270
1,04 1,08 1,12 1,16 1,20 1,30 1,40 1,50 1,60 1,70 1,80 1,90 2,00 2,20 2,40 2,60 2,80 3,00
0,85865 0,87333 0,88679 0,89910 0,91031 0,93401 0,95229 0,96611 0,97635 0,98379 0,98909 0,99279 0,99532 0,99814 0,99931 0,99976 0,99992 0,99998
erf w
2
w 0
e
u2
du
erfc w
1 erf w
58
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla F. Primeras cuatro raíces de la ecuación trascendental,
n·tan( n)
conducción transitoria en una pared plana. Bi =
hL k
0 0,001 0,002 0,004 0,006 0,008 0,01 0,02 0,04 0,06 0,08 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,5 2,0 3,0 4,0 5,0 6,0 7,0 8,0 9,0 10,0 15,0 20,0 30,0 40,0 50,0 60,0 80,0 100,0 ∞
1
0,0000 0,0316 0,0447 0,0632 0,0774 0,0893 0,0998 0,1410 0,1987 0,2425 0,2791 0,3111 0,4328 0,5218 0,5932 0,6533 0,7051 0,7506 0,7910 0,8274 0,8603 0,9882 1,0769 1,1925 1,2646 1,3138 1,3496 1,3766 1,3978 1,4149 1,4289 1,4729 1,4961 1,5202 1,5325 1,5400 1,5451 1,5514 1,5552 1,5708
2
3,1416 3,1419 3,1422 3,1429 3,1435 3,1441 3,1448 3,1479 3,1543 3,1606 3,1668 3,1731 3,2039 3,2341 3,2636 3,2923 3,3204 3,3477 3,3744 3,4003 3,4256 3,5422 3,6436 3,8088 3,9352 4,0336 4,1116 4,1746 4,2264 4,2694 4,3058 4,4255 4,4915 4,5615 4,5979 4,6202 4,6353 4,6543 4,6658 4,7124 59
3
6,2832 6,2833 6,2835 6,2838 6,2841 6,2845 6,2848 6,2864 6,2895 6,2927 6,2959 6,2991 6,3148 6,3305 6,3461 6,3616 6,3770 6,3923 6,4074 6,4224 6,4373 6,5097 6,5783 6,7040 6,8140 6,9096 6,9924 7,0640 7,1263 7,1806 7,2281 7,3959 7,4954 7,6057 7,6647 7,7012 7,7259 7,7573 7,7764 7,8540
4
9,4248 9,4249 9,4250 9,4252 9,4254 9,4256 9,4258 9,4269 9,4290 9,4311 9,4333 9,4354 9,4459 9,4565 9,4670 9,4775 9,4879 9,4983 9,5087 9,5190 9,5293 9,5801 9,6296 9,7240 9,8119 9,8928 9,9667 10,0339 10,0949 10,1502 10,2003 10,3898 10,5117 10,6543 10,7334 10,7832 10,8172 10,8606 10,8871 10,9956
= Bi , para
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
Tabla G. Funciones de Bessel de primera clase. x
J 0( x )
J 1( x )
0,0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1,0 1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 2,0 2,1 2,2 2,3 2,4
1,0000 0,9975 0,9900 0,9776 0,9604 0,9385 0,9120 0,8812 0,8463 0,8075 0,7652 0,7196 0,6711 0,6201 0,5669 0,5118 0,4554 0,3980 0,3400 0,2818 0,2239 0,1666 0,1104 0,0555 0,0025
0,0000 0,0499 0,0995 0,1483 0,1960 0,2423 0,2867 0,3290 0,3688 0,4059 0,4401 0,4709 0,4983 0,5220 0,5419 0,5579 0,5699 0,5778 0,5815 0,5812 0,5767 0,5683 0,5560 0,5399 0,5202
60
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras Tabla H. Funciones de Bessel modificadas de primera y segunda clase. x
- x e ·I 0( x )
0,0 1,0000 0,2 0,8269 0,4 0,6974 0,6 0,5993 0,8 0,5241 1,0 0,4658 1,2 0,4198 1,4 0,3831 1,6 0,3533 1,8 0,3289 2,0 0,3085 2,2 0,2913 2,4 0,2766 2,6 0,2639 2,8 0,2528 3,0 0,2430 3,2 0,2343 3,4 0,2264 3,6 0,2193 3,8 0,2129 4,0 0,2070 4,2 0,2016 4,4 0,1966 4,6 0,1919 4,8 0,1876 5,0 0,1835 5,2 0,1797 5,4 0,1762 5,6 0,1728 5,8 0,1697 6,0 0,1667 6,4 0,1611 6,8 0,1561 7,2 0,1515 7,6 0,1473 8,0 0,1434 8,4 0,1399 8,8 0,1365 9,2 0,1334 9,4 0,1305 9,6 0,1278 10,0 1,0000 I n+1 ( x ) = I n −1 ( x ) − ( 2n / x ) I n ( x )
- x e ·I 1( x )
x e ·K 0( x )
x e · K 1( x )
0,0000 0,0823 0,1368 0,1722 0,1945 0,2079 0,2153 0,2185 0,2190 0,2177 0,2153 0,2121 0,2085 0,2047 0,2007 0,1968 0,1930 0,1892 0,1856 0,1821 0,1788 0,1755 0,1725 0,1695 0,1667 0,1640 0,1614 0,1589 0,1565 0,1542 0,1521 0,1479 0,1441 0,1405 0,1372 0,1341 0,1312 0,1285 0,1260 0,1235 0,1213 0,0000
∞
∞
2,1408 1,6627 1,4167 1,2582 1,1445 1,0575 0,9881 0,9309 0,8828 0,8416 0,8057 0,7740 0,7459 0,7206 0,6978 0,6770 0,6580 0,6405 0,6243 0,6093 0,5953 0,5823 0,5701 0,5586 0,5478 0,5376 0,5280 0,5188 0,5101 0,5019 0,4865 0,4724 0,4595 0,4476 0,4366 0,4264 0,4168 0,4079 0,3995 0,3916
5,8334 3,2587 2,3739 1,9179 1,6362 1,4429 1,3011 1,1919 1,1048 1,0335 0,9738 0,9229 0,8790 0,8405 0,8066 0,7763 0,7491 0,7245 0,7021 0,6816 0,6627 0,6454 0,6292 0,6143 0,6003 0,5872 0,5749 0,5634 0,5525 0,5422 0,5232 0,5060 0,4905 0,4762 0,4631 0,4511 0,4399 0,4295 0,4198 0,4108
61
Transferencia de Calor / Curso 2009-10
Fórmulas, Tablas y Figuras
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
62
PROBLEMAS DE TRANSFERENCIA DE CALOR
Juan Carlos Ramos González Doctor Ingeniero Industrial Raúl Antón Remírez Doctor Ingeniero Industrial Diciembre de 2009
Transferencia de Calor / Curso 2009-10
Problemas
2
Transferencia de Calor / Curso 2009-10
Problemas
ÍNDICE Problemas Tema 1. Introducción a la transferencia de calor y a la conducción ............................. 1 Problemas Tema 2. Conducción unidimensional en régimen estacionario..................................... 7 Problemas Tema 3. Conducción bidimensional en régimen estacionario..................................... 15 Problemas Tema 4. Conducción en régimen transitorio ............................................................... 21 Problemas Tema 5. Introducción a la convección......................................................................... 29 Problemas Tema 6. Convección forzada en flujo externo ............................................................ 33 Problemas Tema 7. Convección forzada en flujo interno ............................................................. 37 Problemas Tema 8. Convección libre o natural ............................................................................ 41 Problemas Tema 9. Introducción a la radiación ............................................................................ 45 Problemas Tema 10. Intercambio radiativo entre superficies ....................................................... 51
i
Transferencia de Calor / Curso 2009-10
Problemas
2
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 1. INTRODUCCIÓN A LA TRANSFERENCIA DE CALOR Y A LA CONDUCCIÓN 1.
(2.7 del Incropera; Ley de Fourier) En el sistema mostrado en la figura se produce una conducción de régimen estacionario unidimensional sin generación de calor. La conductividad térmica es 25 W/m·K y el espesor L es 0,5 m.
T 1
T 2
L
x
Determine las cantidades desconocidas para cada caso de la tabla siguiente y dibuje la distribución de temperatura indicando la dirección del flujo de calor. Caso
T 1
T 2
dT /dx (K/m)
q x′′ (W/m2)
1 400 K 300 K 2 100 ºC -250 3 80 ºC 200 4 -5 ºC 4.000 5 30 ºC -3.000 2 Solución: 1) 200 K/m, -5.000 W/m ; 2) 498 K, 6.250 W/m 2; 3) -20 ºC, -5.000 W/m 2; 4) -85 ºC, -160 K/m; 5) -30 ºC, 120 K/m. 2.
(1.13 del Incropera; Convección) Un chip cuadrado isotérmico de lado 5 mm está montado en un sustrato de manera que sus superficies laterales e inferior están bien aisladas, mientras que la superficie frontal se expone a la corriente de un fluido refrigerante a 15 ºC. La temperatura del chip no debe sobrepasar los 85 ºC. Si el fluido refrigerante es aire ( h = 200 W/m2·K), ¿cuál es la potencia máxima admisible del chip? Si el fluido refrigerante es un líquido dieléctrico (h = 3.000 W/m 2·K), ¿cuál es la potencia máxima admisible del chip? Solución: 0,35 W y 5,25 W.
3.
(Radiación y balance de energía) Un antiguo alumno de la Escuela que trabaja en la ESA (Agencia Espacial Europea) nos ha transmitido la siguiente cuestión: Una sonda de exploración espacial cuyas placas de energía fotovoltaica tienen una superficie A p y una temperatura de fusión T p = 2.000 K es enviada en dirección al Sol. Calcular el radio de la órbita solar mínima ( Ro) a la que se podrá acercar la sonda al Sol. Datos: constante de 1
Transferencia de Calor / Curso 2009-10
Problemas Stefan-Boltzmann σ = 5,67·10 -8 W/m 2·K 4; temperatura de la superficie solar T s = 6.000 K; radio del Sol R s = 7·108 m; suponer que tanto el Sol como las placas se comportan como cuerpos negros (ε = α = 1). Solución: Ro = R s (T s/T p)2. 4.
(Convección y radiación) Una persona desvestida tiene una superficie de 1,5 m2 expuesta a un ambiente y a unos alrededores de 27 ºC. La temperatura de su piel es de 33 ºC y se puede considerar un emisor de radiación perfecto. Si el coeficiente de transferencia de calor por convección es de 9 W/m 2·K, hállese: a) Las pérdidas de calor por convección y por radiación. b) El gasto energético en kcal/día. Solución: a) qconv = Q& conv = 81 W, qrad = Q& rad = 56,8 W; b) 2.846 kcal/día.
5.
(2.6 del Incropera; Ley de Fourier) Para determinar el efecto de la dependencia de la temperatura de la conductividad térmica sobre la distribución de temperatura en un sólido, considere un material para el que esta dependencia puede representarse como k = k o + aT donde k o es una constante positiva y a es un coeficiente que puede ser positivo o negativo. Dibuje la distribución de temperatura de régimen estacionario asociada con la transferencia de calor en una pared plana para tres casos que corresponden a a > 0, a = 0 y a < 0.
6.
(2.11 del Incropera; Ley de Fourier) En el cuerpo bidimensional que se muestra en la figura se encuentra que el gradiente en la superficie A es ∂T /∂ y = 30 K/m. ¿Cuánto valen ∂T /∂ y y ∂T /∂ x en la superficie B?
Solución: ∂T /∂ y = 0; ∂T /∂ x = 60 K/m. 7.
(1.27 del Incropera; Balance de energía) Una placa de aluminio de 4 mm de espesor se monta en posición horizontal con su superficie inferior bien aislada. Se aplica a su superficie superior un recubrimiento que absorbe el 80% de cualquier radiación solar incidente y tiene una emisividad de 0,25. La densidad y el calor específico del aluminio son 2.700 kg/m 3 y 900 J/kg·K, respectivamente. a) Considere las condiciones para las que la placa está a una temperatura de 25 ºC y la superficie superior se expone súbitamente al aire ambiente a T ∞ = 20 ºC y a radiación solar que proporciona un flujo incidente de 900 W/m 2. El coeficiente de transferencia 2
Transferencia de Calor / Curso 2009-10
Problemas de calor por convección entre la superficie y el aire es h = 20 W/m2·K. ¿Cuál es la velocidad inicial de cambio de la temperatura de la placa? Suponga que no hay alrededores. b) ¿Cuál será la temperatura de equilibrio de la placa cuando se alcancen las condiciones de régimen estacionario? c) Represente mediante Excel ® una gráfica de la temperatura de régimen estacionario como función de la emisividad para 0,05 ≤ ε ≤ 1, para tres valores de la absortividad de la placa de 0,5, 0,8 y 1 con el resto de condiciones constantes. Si la finalidad es maximizar la temperatura de la placa, ¿cuál es la combinación más deseable de emisividad y absortividad de la placa? Solución: a) 0,052 K/s; b) T s = 321 K. 8.
(1.31 del Incropera; Balance de energía) En una etapa de un proceso de recocido, 1 hoja de acero inoxidable AISI 304 se lleva de 300 K a 1.250 K conforme pasa a través de un horno calentado eléctricamente a una velocidad de v s = 10 mm/s. El espesor y ancho de la hoja son t s = 8 mm y w s = 2 m, respectivamente, mientras que la altura, ancho y largo del horno son H o = 2 m, W o = 2,4 m y Lo = 25 m, respectivamente. La parte superior y cuatro lados del horno se exponen al aire ambiental y a alrededores a 300 K, y la temperatura de la superficie del horno, su emisividad y el coeficiente de convección respectivos son T s = 350 K, ε s = 0,8 y h = 10 W/m2·K. La superficie inferior del horno también está a 350 K y reposa en una placa de cemento de 0,5 m de espesor cuya base está a 300 K. Estimar la potencia eléctrica que se requiere suministrar al horno. Datos: k cemento (a 300 K) = 1,4 W/m·K. Propiedades termofísicas del acero inoxidable AISI 304: ρ = 7.900 kg/m 3. T (K) c p (J/kg·K) 600 557 800 582
Solución: 841 kW. 9.
(2.12 del Incropera; Ley de Fourier) Algunas secciones del oleoducto de Alaska están tendidas sobre tierra, sostenidas por columnas verticales de acero (k = 25 W/m·K) de 1 m de longitud y sección transversal de 0,005 m 2. En condiciones normales de operación se sabe que la variación de temperatura de un extremo a otro de la longitud de una columna se rige por una expresión de la forma 3
Transferencia de Calor / Curso 2009-10
Problemas T = 100 – 150 x + 10 x2 donde T y x tienen unidades de ºC y metros, respectivamente. Las variaciones de temperatura son insignificantes sobre la sección transversal de la columna. Evalúe la temperatura y la rapidez de conducción de calor en la unión columna-ducto ( x = 0) y en la interfaz columna-tierra ( x = 1). Explique la diferencia en las transferencias de calor. Solución: 18,75 W y 16,25 W. 10.
(2.17 del Incropera; Ley de Fourier) Un aparato para medir la conductividad térmica emplea un calentador eléctrico intercalado entre dos muestras idénticas de 30 mm de diámetro y 60 mm de longitud, prensadas entre placas que se mantienen a una temperatura uniforme T o = 77 ºC mediante la circulación de un fluido. Se pone grasa conductora entre todas las superficies para asegurar un buen contacto térmico. Se empotran termopares diferenciales en las muestras con un espaciado de 15 mm. Las caras laterales de las muestras se aíslan para que la transferencia de calor sea unidimensional.
a) Con dos muestras de acero inoxidable AISI 316 en el aparato, el calentador toma 0,353 A a 100 V y los termopares diferenciales indican ∆T 1 = ∆T 2 = 25,0 ºC. ¿Cuál es la conductividad térmica del material de la muestra de acero inoxidable y cuál la temperatura promedio de las muestras? Compare los resultados con los valores de la Tabla A.1 del Incropera. b) Calcular la conductividad térmica y la temperatura promedio de una muestra de hierro Armco puesta en lugar de la muestra inferior del acero AISI 316. En este caso el calentador toma 0,601 A a 100 V y los termopares diferenciales indican ∆T 1 = ∆T 2 = 15,0 ºC. c) ¿Cuál es la ventaja de construir el aparato con el calentador intercalado entre dos muestras en lugar de construirlo con una sola combinación muestra-calentador? ¿Cuándo resulta significativo el escape de calor por la superficie lateral de las muestras? ¿Bajo que condiciones esperaría que ∆T 1 ≠ ∆T 2? Datos: Propiedades termofísicas del acero inoxidable AISI 316: T (K) k (W/m·K) 300 13,4 400 15,2 4
Transferencia de Calor / Curso 2009-10
Problemas Propiedades termofísicas del Armco: T (K) k (W/m·K) 300 72,7 400 65,7 Solución: a) k = 15 W/m·K y T = 400 K; b) k = 70 W/m·K y T = 380 K. 11.
(2.21 del Incropera; Ecuación de calor) En una varilla cilíndrica de 50 mm de diámetro de combustible de un reactor nuclear ocurre generación interna de calor a q&1 = 5·10 7 W/m3, y en condiciones de régimen estacionario la distribución de temperatura es T (r ) = a+br 2, donde T está en grados Celsius y r en metros, mientras a = 800 ºC y b = -4,167·10 5 ºC/m2. Las propiedades de la varilla de combustible son k = 30 W/m·K, ρ = 1.100 kg/m 3 y c p = 800 J/kg·K. a) ¿Cuál es la velocidad de transferencia de calor por unidad de longitud de la varilla en r = 0 (línea central) y en r = 25 mm (superficie)? b) Si el nivel de potencia del reactor aumenta súbitamente a q& 2 = 10 8 W/m3, ¿cuál es la velocidad de cambio de temperatura en el tiempo inicial en r = 0 y en r = 25? Solución: a) q r ′ ( r = 0) = 0 y q ′r (r = 25) = 9,8 × 10 4 W/m ; b) 56,8 K/s.
12.
(2.24 del Incropera; Ley de Fourier, ecuación de calor y balance de energía) Un estanque solar poco profundo con gradiente salino consiste en tres capas fluidas distintas y se utiliza para absorber energía solar. Las capas superior e inferior están bien mezcladas y sirven para mantener las superficies superior e inferior de la capa central a temperaturas uniformes T 1 y T 2, donde T 2 > T 1. Considere condiciones para las que la absorción de la radiación solar en la capa central proporciona una generación no uniforme de calor de la − forma q& ( x ) = Ae ax , y la distribución de temperatura en la capa central es: T ( x ) = −
A ka
2
e
− ax
+ Bx + C
Las cantidades A (W/m3), a (1/m), B (K/m) y C (K) son constante conocidas, y k es la conductividad térmica que también es constante.
a) Obtenga expresiones para la rapidez a la que se transfiere calor por unidad de área de la capa inferior a la capa central y de ésta a la capa superior. 5
Transferencia de Calor / Curso 2009-10
Problemas
b) Determine si las condiciones son estacionarias o transitorias. c) Obtenga una expresión para la rapidez a la que se genera energía térmica en la capa central, por unidad de área superficial. Solución: a) q x′′ ( x = L ) = −
& ′′ = E gen
A a
A a
e
− aL
− Bk ; q x′′ ( x = 0) = −
(1 − e − ). aL
6
A a
− Bk ; b) Régimen estacionario; c)
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 2. CONDUCCIÓN UNIDIMENSIONAL EN RÉGIMEN ESTACIONARIO Pared plana 1.
(3.2 del Incropera) La ventana posterior de un automóvil se desempaña mediante el paso de aire caliente sobre su superficie interna. a) Calcular las temperaturas de las superficies interna y externa de una ventana de vidrio de 4 mm de espesor, siendo la temperatura del aire caliente T ∞,int = 40 ºC y su coeficiente de convección hint = 30 W/m 2·K y la temperatura del aire exterior T ∞,ext = 10 ºC y su coeficiente de convección hext = 65 W/m2·K. b) Evalúe cualitativamente la influencia de T ∞,ext y hext sobre las temperaturas. Datos: k vidrio (a 300 K) = 1,4 W/m·K. Solución: a) T int = 7,7 ºC y T ext = 4,9 ºC; b) Ambas disminuyen al aumentar hext y aumentan al aumentar T ∞,ext .
2.
(3.3 del Incropera) En la ventana posterior del automóvil del problema anterior se instala como sistema para desempañar su superficie interior un elemento de calentamiento consistente en una película transparente delgada con resistencias eléctricas. Al calentarse eléctricamente este dispositivo se establece un flujo de calor uniforme en la superficie interna. a) Calcular la potencia eléctrica por unidad de área de ventana necesaria para mantener la temperatura de la superficie interna a 15 ºC cuando la temperatura del aire interior es T ∞,int = 25 ºC y su coeficiente de convección hint = 10 W/m2·K. El aire exterior está en las mismas condiciones que en el problema anterior. b) Calcular la temperatura de la superficie externa de la ventana. c) Evalúe cualitativamente la influencia de T ∞,ext y hext sobre la potencia eléctrica.
′′ = 1,27 kW/m 2; b) T ext = 11,1 ºC; c) P elec ′′ aumenta al aumentar hext y Solución: a) P elec disminuye al aumentar T ∞,ext . 3.
(3.15 del Incropera) Una casa tiene una pared compuesta de madera, aislante de fibra y tablero de yeso, como se indica en el esquema. En un día frío de invierno los coeficientes de transferencia de calor por convección son hext = 60 W/m2·K y hint = 30 W/m2·K. El área total de la superficie es de 350 m2. Datos: Tablero de yeso: k (a 300 K) = 0,17 W/m·K. Propiedades termofísicas de la fibra de vidrio: T (K) ρ (kg/m3) k (W/m·K) 300 16 0,046 300 28 0,038 300 40 0,035 7
Transferencia de Calor / Curso 2009-10
Problemas Tablero de madera contraplacada: k (a 300 K) = 0,12 W/m·K.
a) Determine una expresión simbólica para la resistencia térmica total de la pared incluyendo los efectos de convección. b) Determine la pérdida de calor total de la pared. c) Si el viento soplara de manera violenta elevando hext a 300 W/m 2·K, ¿cuál sería el porcentaje de aumento relativo de la pérdida de calor? d) ¿Qué resistencia térmica influye en mayor medida sobre la pérdida de calor a través de la pared? Solución: b) 4.214 W; c) 0,45 %; d) La de la fibra de vidrio, que es el aislante y tiene la k menor. Resistencia de contacto 4.
(3.25 del Incropera) Un circuito integrado (chip) disipa 30.000 W/m 2 de calor eléctrico. El chip, que es muy delgado, se expone a un líquido dieléctrico en su superficie superior con hext = 1.000 W/m2·K y T ∞,ext = 20 ºC. En la superficie inferior se une a una tarjeta de circuitos de espesor Lb = 5 mm y conductividad k b = 1 W/m·K. La resistencia térmica de contacto entre el chip y la tarjeta es Rt ′′,c = 10-4 m2·K/W. La superficie inferior de la tarjeta se expone al aire ambiente para el que hint = 40 W/m2·K y T ∞,int = 20 ºC. a) Dibuje el circuito térmico equivalente señalando las resistencias térmicas, las temperaturas y los flujos de calor. b) ¿Cuál es la temperatura del chip para las condiciones de disipación de qc′′ = 30.000 W/m2? c) ¿Qué influencia tendría en la temperatura del chip el aumentar en un orden de magnitud la conductividad de la tarjeta de circuitos y en disminuir en un orden de magnitud la resistencia térmica de contacto entre el chip y la tarjeta? Solución: b) T c = 49 ºC; c) Prácticamente ninguna.
8
Transferencia de Calor / Curso 2009-10
Problemas Pared cilíndrica 5.
(3.37 del Incropera) Un calentador eléctrico delgado se inserta entre una varilla circular larga y un tubo concéntrico con radios interior y exterior de 20 y 40 mm. La varilla A tiene una conductividad térmica k A = 0,15 W/m·K y el tubo B k B = 1,5 W/m·K. La superficie externa está en contacto con un fluido a temperatura T ∞ = -15 ºC y un coeficiente de convección de 50 W/m 2·K. a) Determine la potencia eléctrica por unidad de longitud de los cilindros que se requieren para mantener la superficie externa del tubo B a 5 ºC. b) ¿Cuál es la temperatura en el centro de la varilla A? Solución: a) 251 W/m; b) 23,4 ºC.
6.
(3.44 del Incropera) Una corriente eléctrica de 700 A fluye a través de un cable de acero inoxidable que tiene un diámetro de 5 mm y una resistencia eléctrica de 6·10 -4 Ω/m. El cable esté en un medio que tiene una temperatura de 30 ºC y el coeficiente total asociado con la convección y la radiación entre el cable y el medio es aproximadamente 25 W/m 2·K. a) Si el cable está expuesto, ¿cuál es la temperatura de la superficie? b) Si se aplica un recubrimiento muy delgado de aislante eléctrico al cable, con una resistencia de contacto de 0,02 m2·K/W, ¿cuáles son las temperaturas superficiales del aislante y del cable? c) Si se usa un aislante de conductividad térmica 0,5 W/m·K, ¿cuál será el espesor de este aislante que dará el valor más bajo de la temperatura del cable? ¿Cuál es el valor de esa temperatura? Solución: a) T s,cable = 778,7 ºC; b) T s,cable = 1.153 ºC y T s,aislante = 778,7 ºC; c) e = 17,5 mm y T s,cable = 318,2 ºC.
7.
(3.45 del Incropera) Un tubo de acero de pared delgada de 0,20 m de diámetro y emisividad 0,8 se utiliza para transportar vapor saturado a una presión de 20 bar ( T sat = 485 K) en un cuarto para el que la temperatura del aire y de las paredes es 25 ºC y el coeficiente de transferencia de calor por convección en la superficie externa del tubo es 20 W/m2·K. a) ¿Cuál es la pérdida de calor por unidad de longitud del tubo expuesto (sin aislante)? b) Calcule la pérdida de calor por unidad de longitud del tubo si se añade una capa aislante de 50 mm de óxido de magnesio que también tiene una emisividad de 0,8. Calcular también la temperatura superficial exterior del aislante. c) El coste asociado con la generación del vapor saturado es de 4 €/10 9 J y el del aislante y su instalación de 100 €/m. Si la línea de vapor opera 7.500 horas al año, ¿cuánto tiempo se necesita para amortizar la instalación del aislante? Datos: Propiedades termofísicas del óxido de magnesio: T (K) k (W/m·K) 310 0,051 365 0,055 420 0,061 Solución: a) 3.702 W/m; b) q ′ ≈ 162 W/m y T s,ail,ext ≈ 30 ºC; c) ≈ 3 meses. 9
Transferencia de Calor / Curso 2009-10
Problemas Pared esférica 8.
(3.56 del Incropera) Una sonda esférica crioquirúrgica se incrusta en tejido enfermo con el propósito de congelarlo y destruirlo. La sonda tiene un diámetro de 3 mm y su superficie se mantiene a -30 ºC cuando se incrusta en tejido que está a 37 ºC. Se forma una capa esférica de tejido congelado alrededor de la sonda con una temperatura de 0 ºC en su superficie de contacto con el tejido normal. Si la conductividad térmica del tejido congelado es 1,5 W/m·K y el coeficiente de transferencia de calor por convección entre el tejido congelado y el normal es 50 W/m 2·K, ¿cuál es el espesor de la capa del tejido congelado? Resolución: En primer lugar se expone una manera de obtener una expresión de la resistencia térmica para una esfera. En un elemento diferencial de esfera la aplicación de la conservación de la energía implica que qr = qr+dr , es decir que la transferencia de calor es independiente del radio, para condiciones unidimensionales de régimen estacionario y sin generación interna de calor. La ecuación de Fourier para una esfera hueca cuyas superficies están en contacto con fluidos a temperaturas distintas y en condiciones de régimen estacionario sin generación de calor adopta la forma: qr = − kA
dT dr
= −k (4π r 2 )
dT dr
donde A = 4πr 2 es el área normal a la dirección de la transferencia de calor. Al integrar la ecuación anterior: qr 4π
r 2
dr
r 1
r 2
∫
T s 2
= − ∫ k (T )dT T s 1
Suponiendo k constante y al resolver para las condiciones de contorno de temperaturas conocidas en las superficies se obtiene: qr =
4π k (T s1 − T s 2 ) (1 / r 1 ) − (1 / r 2 )
La resistencia térmica para conducción adopta, por la tanto, la forma: Rt ,cond =
Y la de convección: Rt ,conv =
(T s1 − T s 2 ) qr
=
1 1
1 − 4π k r 1 r 2
1 4π r 2 h
Una vez visto esto se puede representar el circuito térmico equivalente al enunciado del problema: q
T s1
T s2 1 1 1 − 4π k r r + e
T ∞ 1 4π (r + e) 2 h
10
q
Transferencia de Calor / Curso 2009-10
Problemas La velocidad de transferencia de calor se puede expresar: q=
T ∞ − T s 2 1 4π h( r + e)
= 2
T s 2 − T s1 1 1 1 − 4π k r r + e
Al resolver se obtiene la siguiente ecuación en e2: e 2 + re −
k T s 2 − T s1 h T ∞ − T s 2
r = 0
Al resolver se obtiene: e = 5,34 mm. Generación interna de calor 9.
(3.73 del Incropera) El aire dentro de una cámara a T ∞,int = 50 ºC se calienta convectivamente con hint = 20 W/m2·K mediante una pared de 200 mm de espesor que tiene una conductividad térmica de 4 W/m·K y una generación de calor uniforme de 1.000 W/m3. Para prevenir que algo del calor generado se pierda hacia el exterior de la cámara, a T ∞,ext = 25 ºC con hext = 5 W/m2·K, se coloca un calentador de listón muy delgado sobre la pared exterior para proporcionar un flujo de calor uniforme, q o′′ .
a) Dibuje la distribución de temperaturas en la pared (T - x) para la condición de que no se pierde nada del calor generado dentro de la pared hacia el exterior de la cámara (es decir, quitando el calentador y aislando la superficie externa de la pared). b) ¿Cuáles son las temperaturas en las superficies externa e interna de la pared para esa condición? c) Determine el valor de qo′′ que debe suministrar el calentador de listón de modo que todo el calor generado dentro de la pared se transfiera al interior de la cámara. d) Si la generación de calor en la pared se cortara mientras el flujo de calor del calentador de listón permanece constante, ¿cuál sería la temperatura de la pared exterior en régimen permanente? Solución: b) T (0) = 65 ºC y T ( L) = 60 ºC; c) 200 W/m 2; d) 55 ºC. 10.
(3.83 del Incropera) Un elemento de combustible de reactor nuclear consiste en un núcleo cilíndrico sólido de radio r 1 y conductividad térmica k f . El núcleo de combustible está en buen contacto con un material de encamisado de radio externo r 2 y conductividad térmica k c. Considere condiciones de régimen estacionario para las que ocurre una generación de calor uniforme dentro del combustible a una razón volumétrica q& = e& gen y la superficie 11
Transferencia de Calor / Curso 2009-10
Problemas externa del encamisado se expone a un fluido refrigerante que se caracteriza por una temperatura T ∞ y un coeficiente de convección h. a) Obtenga expresiones para las distribuciones de temperatura en el combustible y en el encamisado, T f (r ) y T c(r ). b) Considere un núcleo de combustible de óxido de uranio para el que K f = 2 W/m·K y r 1 = 6 mm y un encamisado para el que K c = 25 W/m·K y r 2 = 9 mm. Si q& = e& gen = 2·108 W/m3, h = 2.000 W/m2·K y T ∞ = 300 K, ¿cuál es la temperatura máxima en el elemento de combustible? c) Evalúe cualitativamente la influencia de h sobre las temperaturas. ¿Es posible mantener la temperatura de la línea central del combustible por debajo de 1.000 K ajustando el flujo de refrigerante y, por tanto, el valor de h? Solución: b) T f (r = 0) = 1.458 K; c) Si h aumenta T f y T c disminuyen. No es posible. Superficies extendidas y aletas 11.
(3.109 del Incropera) Varillas de cobre circulares de diámetro D = 1 mm y longitud L = 25 mm se usan para reforzar la transferencia de calor de una superficie que se mantiene a T s1 = 100 ºC. Un extremo de la varilla se une a esta superficie (en x = 0) y el otro ( x = 25) se une a una segunda superficie que se mantiene T s2 = 0 ºC. El aire que fluye entre las superficies también está a una temperatura T ∞ = 0 ºC y tiene un coeficiente de convección h = 100 W/m2·K. a) ¿Cuál es la transferencia de calor de una sola varilla de cobre? b) ¿Cuál es la transferencia total de calor de una sección de 1 m x 1 m de la superficie a 100 ºC, si se instala una disposición de varillas separadas entre centros 4 mm? Datos: k cobre (a 300 K) = 401 W/m·K. Solución: a) q f = 1,51 W; b) qt = 103,8 kW.
12.
(3.114 del Incropera) A menudo se forman pasajes de aletas entre placas paralelas para reforzar la transferencia de calor por convección en núcleos compactos de intercambiadores de calor. Considere una pila de aletas de 200 mm de ancho y 100 mm de profundidad con 50 aletas de 12 mm de longitud. La pila completa está fabricada de aluminio (k = 240 W/m·K) de 1 mm de espesor. Las temperaturas máximas permisibles asociadas a las placas opuestas son T o = 400 K y T L = 350 K. El aire que fluye entre las placas tiene una h = 150 W/m2·K y una T ∞ = 300 K. ¿Cuáles son las disipaciones de calor de una aleta y del sistema de aletas en cada una de las placas?
12
Transferencia de Calor / Curso 2009-10
Problemas
Solución: q fo = 114,95 W; q fL = -88,08 W; qto = 5.972,5 W; qtL = -4.291,5 W. 13.
(3.131 del Incropera modificado, examen septiembre 2005) Se quiere disipar el calor generado en el interior de un transformador situando en una de sus paredes un dispositivo de aletas rectas. La pared del transformador tiene una conductividad térmica de 5 W/m·K y un espesor de 6 mm. Sobre ella se coloca un dispositivo de aletas de sección rectangular de aluminio (k al = 240 W/m·K). El soporte del dispositivo de aletas tiene un espesor de 4 mm. Entre la pared del transformador y el soporte de las aletas hay una resistencia de contacto − de valor Rt ′′,c = 10 4 m 2 ·K/W . Las aletas tienen una longitud de 25 mm, un espesor de 2 mm y la distancia entre ellas es de 2 mm. El calor generado en el transformador se puede asimilar a un flujo de calor uniforme sobre la pared de valor qi′′ = 10 5 W/m 2 . El aire
exterior está a 320 K y proporciona un coeficiente de convección de 100 W/m 2·K.
Pared del transformador k tra = 5 W/m·K
R" t,c
Soporte de k = 240 W/m·K las aletas al T al
T b
T int
t = 2 mm
= 2 mm
q" i
h = 100 W/m2·K T = 320 K
6 mm 4 mm
25 mm 13
Transferencia de Calor / Curso 2009-10
Problemas Problemas a) Dibuje el circuito térmico equivalente entre el interior del transformador y el aire exterior para la parte de pared que le corresponde a una aleta teniendo en cuenta que la dimensión perpendicular al dibujo es muy larga. b) Calcule los valores de las resistencias térmicas que aparecen en el circuito térmico anterior. c) Calcule la temperatura de la superficie interna del transformador, T int int . d) Calcule la temperatura de la superficie interna del soporte de aluminio (en contacto con la resistencia de contacto), T al al . e) Calcule la temperatura de la base de las aletas, T b. Solución: Solución: c) T int int = 532,3 K; d) T al al = 402,3 K; e) T b = 400,6 K.
14
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 3. CONDUCCIÓN BIDIMENSIONAL EN RÉGIMEN ESTACIONARIO Factores de forma 1.
(4.16 del Incropera) Con las relaciones de resistencia térmica desarrolladas en el tema 3 determine expresiones del factor de forma para las siguientes geometrías: a) Pared plana, capa cilíndrica y coraza esférica. b) Esfera hueca de superficie isotérmica de diámetro D diámetro D en en el interior de un medio infinito. Solución: Solución : a) A a) A// L, L, 2π L/ln( L/ln(r r 2/r 1), 4πr 1r 2/(r /(r 2-r 1); b) 2π D. D.
2.
(4.20 del Incropera) Un cable largo de transmisión de energía se entierra a una profundidad (distancia de la tierra a la línea central del cable) de 2 m. El cable está enfundado en un tubo de pared delgada de 0,1 m de diámetro y para hacer al cable superconductor (esencialmente cero disipación de energía), el espacio entre el cable y el tubo está lleno de nitrógeno líquido a 77 K. Si el tubo se cubre con un superaislante ( k i = 0,005 W/m·K) de 0,05 m de espesor y la superficie de la tierra ( k g = 1,2 W/m·K) está a 300 K, ¿cuál es la carga de enfriamiento por unidad de longitud de tubo [W/m] que debe suministrar un refrigerador criogénico para mantener el nitrógeno a 77 K? Solución: Solución: 9,89 W/m.
3.
(4.25 del Incropera) Por un tubo de cobre de pared delgada de 30 mm de diámetro fluye agua caliente a 85 ºC. El tubo está forrado de una capa cilíndrica excéntrica que se mantiene a 35 ºC y mide 120 mm de diámetro. La excentricidad, definida como la distancia entre los centros del tubo y la capa, es 20 mm. El espacio entre el tubo y la capa está lleno de un material aislante que tiene una conductividad térmica de 0,05 W/m·K. Calcule la pérdida de calor por unidad de longitud de tubo y compare el resultado con la pérdida de calor para una disposición concéntrica. Solución: Solución: 12,5 W/m y 11,33 W/m. Factores de forma con circuitos térmicos
4.
(4.28 del Incropera) Un fluido caliente pasa por tubos circulares de una plancha de hierro colado de espesor L espesor LA = 60 mm que está en contacto con unas placas de cubierta de espesor LB = 5 mm. Los canales tienen un diámetro D diámetro D = = 15 mm con un espaciado de línea central de L de Lo = 60 mm. Las conductividades térmicas de los materiales son k A = 20 W/m·K y k B = 2 75 W/m·K, y la resistencia de contacto entre los dos materiales es Rt ′′,c = 2·10 −4 m ·K/W. El 2
fluido caliente está a T i = 150 ºC y el coeficiente de convección es 1.000 W/m ·K. Las placas de cubierta se exponen al aire ambiental que está a 25 ºC y tiene un coeficiente de convección de 200 W/m2·K.
15
Transferencia de Calor / Curso 2009-10
Problemas Problemas
a) Determine la transferencia de calor de un solo tubo por unidad de longitud de la plancha en dirección dirección normal normal a la página, qi′ . b) Determine la temperatura de la superficie externa de la placa de cubierta, T s. c) Comente los efectos sobre q i′ y T s de un cambio en el espaciado de los canales. ¿Cómo afectaría a q i′ y T s aislar la superficie inferior? Solución: Solución: a) qi′ = 1.578,6 W/m; b) T s = 90,8 ºC; c) Si Lo aumenta, q i′ aumenta y T s disminuye. Si la superficie inferior está aislada, qi′ disminuye y T s permanece constante. 5.
(4.31 del Incropera) En el Tema 3 se supuso que cuando se une una aleta a un material base, la temperatura temperatura de la base no cambia. Lo que en verdad ocurre es que, si la temperatura del material de la base excede la temperatura del fluido, al colocar una aleta disminuye la temperatura de la unión, T j, por debajo de la de la base y el flujo de calor del material de la base a la aleta es bidimensional. Considere condiciones en las que una aleta larga circular de aluminio de diámetro D D = 5 mm se une al material de la base cuya temperatura lejos de la unión se mantiene a T b = 100 ºC. Las condiciones de convección en la superficie de la aleta son T ∞ = 25 ºC y h = 50 W/m2·K.
16
Transferencia de Calor / Curso 2009-10
Problemas a) Calcule la temperatura de la unión y la transferencia de calor cuando el material de la base es (i) aluminio aluminio (k (k = = 240 W/m·K) y (ii) acero inoxidable (k ( k = = 15 W/m·K). b) Repita los cálculos anteriores para el caso del aluminio si entre la unión de la aleta y el material de la base hay una resistencia térmica Rtc" = 3·10 −5 m2·K/W c) ¿Cómo influye el coeficiente de convección en la transferencia de calor?
Solución: a) (i) T j = 98 ºC, q f = Q& f = 4,8 W; (ii) T j = 78,4 ºC, q f = Q& f = 3,24 W; b) T j = 92 ºC, q f = Q& f = 4,08 W; c) Si h aumenta q f aumenta. 6.
(4.32 del Incropera) Se construye un iglú en forma de hemisferio con un radio interno de 1,8 m y paredes de nieve compactada de 0,5 m de espesor. En el interior del iglú el coeficiente de transferencia de calor por convección es 6 W/m2·K; en el exterior, en condiciones normales de viento, es 15 W/m 2·K. La conductividad térmica de la nieve compactada es 0,15 W/m·K. La temperatura de la capa de hielo sobre la que se asienta el iglú es de -20 ºC y tiene la misma conductividad térmica que la nieve compactada.
a) Suponiendo que el calor corporal de los ocupantes proporciona una fuente continua de 320 W dentro del iglú, calcule la temperatura del aire interior cuando la del aire exterior es -40 ºC. Considere las pérdidas de calor a través del suelo. b) ¿Cómo afecta a la temperatura interior el que el coeficiente de convección exterior se triplique debido al viento? ¿Y cómo afecta el doblar el espesor de las paredes?
Solución: a) T i = 1,2 ºC; b) T i = 0,8 ºC; T i = 20,8 ºC. 7.
(4.34 del Incropera) Un dispositivo electrónico en forma de disco de 20 mm de diámetro disipa 100 W cuando se monta sobre un bloque grande de aleación de aluminio (2024-T6) cuya temperatura se mantiene a 27 ºC. En la interfaz entre el dispositivo y el bloque hay − una resistencia de contacto Rtc" = 5·10 5 m2·K/W.
17
Transferencia de Calor / Curso 2009-10
Problemas
a) Calcule la temperatura que alcanzará el dispositivo suponiendo que toda la potencia que genera debe transferirse por conducción al bloque. b) Para aumentar la potencia del dispositivo se instala un sistema de aletas en la parte superior del dispositivo. Las aletas rectas de sección circular (aletas de aguja) están hechas de cobre (k = 400 W/m·K) y están expuestas a un flujo de aire a 27 ºC para el que el coeficiente de convección es 1.000 W/m 2·K. Para la temperatura del dispositivo que se calculó en el apartado a), ¿cuál es la potencia de operación permisible? Datos: Propiedades termofísicas de la aleación de aluminio 2024-T6:
T (K) k (W/m·K) 200 163 300 177 400 186 Solución: a) T d = 57 ºC; b) P eléct. = 138,65 W. Método de las diferencias finitas 8.
(4.41 del Incropera) Las superficies superior e inferior de una barra de conducción se enfrían convectivamente por acción de aire a T ∞, pero con h sup ≠ hinf . Los lados se enfrían manteniendo contacto con sumideros de calor a T o, a través de una resistencia térmica de contacto Rt ′′,c . La barra tiene conductividad térmica k y el ancho es el doble del espesor L. Considere condiciones de estado estacionario para las que se genera calor de manera uniforme a una tasa volumétrica q& debido al paso de una corriente eléctrica. Obtenga las ecuaciones en diferencias finitas para los nodos 1 y 13.
18
Transferencia de Calor / Curso 2009-10
Problemas
9.
(4.48 del Incropera) Las temperaturas de estado estacionario (en K) en tres puntos nodales de una varilla rectangular son como se muestra en la figura. La varilla experimenta una rapidez de generación de calor volumétrica uniforme de 5·10 7 W/m3 y tiene una conductividad térmica de 20 W/m·K. Dos de sus lados se mantienen a una temperatura constante de 300 K, mientras que los otros dos están aislados. a) Determine las temperaturas en los nodos 1, 2 y 3 resolviendo el sistema de 3 ecuaciones con 3 incógnitas que forman las ecuaciones nodales. b) Calcule la transferencia de calor por unidad de longitud de la varilla (W/m) a partir de las temperaturas nodales. Compare este resultado con la transferencia de calor calculada a partir del conocimiento de la generación volumétrica de calor y las dimensiones de la varilla.
Solución: a) T 1 = 362,4 K; T 2 = 390,2 K; T 3 = 369 K; b) q’ = 7.502,5 W/m; q’ = 7.500 W/m.
19
Transferencia de Calor / Curso 2009-10
Problemas 10.
(4.52 del Incropera) Una barra larga de sección transversal rectangular tiene 60 mm de ancho, 90 mm de largo, y una conductividad térmica de 1 W/m·K. Uno de sus anchos está sometido a un proceso de convección con aire a 100 ºC y un coeficiente de convección de 100 W/m2·K. El resto de los lados se mantiene a 50 ºC.
a) Con un espaciado de malla de 30 mm y mediante el método iterativo de Gauss-Seidel, determine las temperaturas nodales y la transferencia de calor (por unidad de longitud normal a la página) desde el aire a la barra. b) Utilizando Matlab® para resolver el sistema de ecuaciones (método de inversión de matrices), repita los cálculos con un espaciado de malla de 15 mm.
Solución: a) Empezando desde el lado sometido a convección las temperaturas nodales son: 81,7 ºC, 58,5 ºC y 52,1 ºC; q’ = 205 W/m; b) Temperaturas de los nodos a lo largo del ancho sometido a convección: 50 ºC, 80,33 ºC, 85,16 ºC, 80,33 ºC y 50 ºC; q’ = 156,27 W/m. 11.
(Basado en Ejemplo 5.1 del Chapman, 5ª edición) Se dispone de una varilla de hierro ( k = 50 W/m·K) de 1 cm de diámetro y 20 cm de longitud. La varilla se une en un extremo a una superficie calentada a 120 ºC y en el extremo libre se encuentra aislada. Su superficie lateral está en contacto con un fluido a 20 ºC para el que el coeficiente de transferencia de calor por convección es 10 W/m 2·K. a) Determine la distribución de temperaturas en la varilla resolviendo las ecuaciones nodales mediante el método iterativo de Gauss-Seidel y con un ∆ x = 5 cm. El número de iteraciones viene dado por un criterio de convergencia en la temperatura del extremo de un 1 %. Es decir, la diferencia relativa de la temperatura en el extremo en dos iteraciones sucesivas ha de ser inferior al 1 %.
E r (%) =
i i −1 − T extremo T extremo i T extremo
100 ≤ 1 %
b) Determine una aproximación a la pérdida de calor de la varilla a partir de la distribución discreta de temperaturas calculada en el apartado anterior. c) Compare la solución anterior con el calor perdido por la aleta calculado de manera exacta.
Solución: a) Comenzando desde la base la distribución de temperaturas es: 120 ºC, 84,8 ºC, 63,6 ºC, 52,8 ºC, 49,9 ºC; b) q faprox. = 3,24 W; c) q f = 3,32 W. 20
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 4. CONDUCCIÓN EN RÉGIMEN TRANSITORIO Método de la resistencia interna despreciable 1.
(5.8 del Incropera) Una bala esférica de plomo de 6 mm de diámetro se mueve aproximadamente a Mach 3. La onda de choque resultante calienta el aire alrededor de la bala a 700 K, y el coeficiente de convección promedio para la transferencia de calor entre el aire y la bala es 500 W/m 2·K. Si la bala sale de la escopeta a 300 K y el tiempo de vuelo es 0,4 s, ¿cuál es la temperatura en la superficie en el momento del impacto? Datos: Propiedades termofísicas del plomo a 300 K: k = 35,3 W/m·K; ρ = 11.340 kg/m 3; c p = 129 J/kg·K.
Solución: T = 351 K. 2.
(5.10 del Incropera) Una unidad de almacenamiento de energía térmica consiste en un canal rectangular largo, que está bien aislado en la superficie externa y encierra capas alternadas del material de almacenamiento y rejillas para el flujo.
Cada capa del material de almacenamiento es una plancha de aluminio de ancho W = 0,05 m que está a una temperatura inicial de 25 ºC. Considere condiciones en las que la unidad de almacenamiento se carga con el paso de un gas caliente a través de las rejillas, suponiendo que la temperatura del gas y el coeficiente de convección tienen valores constantes de T ∞ = 600 ºC y h = 100 W/m2·K a lo largo del canal. ¿Cuánto tiempo se tardará en alcanzar el 75 % del almacenamiento máximo posible de energía? ¿Cuál es la temperatura del aluminio en ese momento? Datos: Propiedades termofísicas del aluminio: ρ = 2.702 kg/m 3.
T (K) k (W/m·K) c p (J/kg·K) 300 237 903 400 240 949 600 231 1.033 Solución: t = 933,5 s ≈ 15,55 min y T = 456 ºC. 3.
(5.14 del Incropera) La pared plana de un horno se fabrica de acero al carbono simple ( k = 60 W/m·K; ρ = 7.850 kg/m 3; c p = 430 J/kg·K) y tiene un espesor de L = 10 mm. Para protegerla de los efectos corrosivos de los gases de combustión del horno, una superficie 21
Transferencia de Calor / Curso 2009-10
Problemas de la pared se cubre con una película delgada de cerámica que, para un área superficial unitaria, tiene una resistencia térmica de Rtf ′′ = 0,01 m 2·K/W. La superficie opuesta está bien aislada de los alrededores.
Al poner en funcionamiento el horno, la pared está a una temperatura inicial de T i = 300 K y los gases de combustión entran en el horno a T ∞ = 1.300 K, con lo que proporcionan un coeficiente de convección de 25 W/m2·K en la película cerámica. Suponiendo que la película tiene una resistencia térmica interna insignificante, ¿cuánto tiempo tardará la superficie interior del acero en alcanzar una temperatura de T si = 1.200 K? ¿Cuál es la temperatura T so de la superficie expuesta de la película cerámica en ese momento?
Resolución: Se dibuja el circuito térmico equivalente del sistema:
T so
T ∞
q
T si
1
Rtf ′′
hA
A
Como entre la pared y el fluido existe una película que aporta una resistencia térmica de contacto, para poder calcular el número de Biot y estudiar si se puede aplicar el método de la resistencia interna despreciable hay que trabajar con el coeficiente global de transferencia de calor, U : UA =
1
∑ R
⇒ U =
t
1 1 h + Rtf ′′
=
1 1 / 25 + 0,01
= 20 W/m 2 ·K
El número de Biot correspondiente será: Bi =
ULc k
=
20·0,01 60
= 0,0033
Al ser menor que 0,1 se puede aplicar el método de la resistencia interna despreciable: T si − T ∞
t = exp − T i − T ∞ τ 22
Transferencia de Calor / Curso 2009-10
Problemas τ = Rt C c = t = −τ ln
ρ Vc p
UA s
T si − T ∞ T i − T ∞
=
7.850·0,01·430 20
= 1.687,75·ln
= 1.687,75 s
1.200 − 1.300 300 − 1.300
= 3.886 s ≈ 65 min
A partir del circuito térmico equivalente, la temperatura de la superficie de la película vendrá dada por: h(T ∞ − T so ) = (T so − T si ) / Rtf ′′ ⇒ T so =
hT ∞ + T si / Rtf ′′ h + 1 / Rtf ′′
= 1.220 K
Efectos espaciales: análisis de semejanza 4.
(5.28 del Incropera) Considere la pared unidimensional que se muestra en el dibujo que inicialmente se encuentra a temperatura uniforme T i y se somete de pronto a la condición de frontera de convección con un fluido a T ∞.
Para una pared en particular, caso 1, la temperatura en x = L1 después de t 1 = 100 s es T 1( L1,t 1) = 315 ºC. Otra pared, caso 2, tiene diferentes condiciones de espesor y térmicas como se muestra en la siguiente tabla. Caso
L (m)
2
(m /s)
k (W/m·K)
T i (ºC)
T
(ºC)
2
h (W/m
·K)
-6
1 0,10 15·10 50 300 400 200 -6 2 0,40 25·10 100 30 20 100 ¿Cuánto tiempo tardará la segunda pared en alcanzar 28,5 ºC en la posición x = L2? Solución: t = 960 s = 16 min.
Conducción unidimensional: pared plana 5.
(5.32 del Incropera) Considere la unidad de almacenamiento de energía del problema 2, 3 pero con un material de mampostería de k = 0,70 W/m·K; ρ = 1.900 kg/m ; c p = 800 J/kg·K empleado en lugar del aluminio. ¿Cuánto tiempo se tardará en alcanzar el 75 % del almacenamiento máximo posible de energía? ¿Cuáles son las temperaturas máxima y mínima de la mampostería en ese momento? Solución: t = 1.174 s ≈ 19,6 min; T mín ( x = 0) = 412 ºC; T máx ( x = W /2) = 538 ºC. 23
Transferencia de Calor / Curso 2009-10
Problemas
Conducción unidimensional: cilindro 6.
(5.45 del Incropera) Una varilla larga de 40 mm de diámetro fabricada de zafiro (óxido de aluminio) e inicialmente a una temperatura uniforme de 800 K, se enfría de súbito con un 2 fluido a 300 K que tiene un coeficiente de transferencia de calor de 1.600 W/m ·K. Después de 35 segundos la varilla se envuelve en un aislante y no experimenta pérdidas de calor. ¿Cuál será la temperatura de la varilla después de un largo tiempo? Datos: Propiedades termofísicas del óxido de aluminio (zafiro) a T = 600 K: k = 18,9 3
W/m·K; ρ = 3.970 kg/m ; c p = 1.110 J/kg·K. Solución: T inf = 510 K.
Conducción unidimensional: esfera 7.
(5.48 del Incropera) En el tratamiento térmico para endurecer bolas de acero de rodamientos (k = 50 W/m·K; ρ = 7.800 kg/m3; c p = 500 J/kg·K) se desea aumentar la temperatura de la superficie por un tiempo corto sin calentar de manera significativa el interior de la bola. Este tipo de calentamiento se lleva a cabo mediante la inmersión súbita 2 de la bola en un baño de sal derretida con T ∞ = 1.300 K y h = 5.000 W/m ·K. Suponga que cualquier posición dentro de la bola cuya temperatura exceda 1.000 K se endurecerá. Estime el tiempo que se necesita para endurecer el milímetro externo de una bola de 20 mm de diámetro si su temperatura inicial es de 300 K. Solución: t = 3,4 s.
Sólido semiinfinito 8.
(5.61 del Incropera) Una grúa para levantar losas adheridas al suelo emplea un bloque de hierro que se mantiene a temperatura constante de 150 ºC mediante un calentador eléctrico empotrado. El bloque de hierro se pone en contacto con la losa para suavizar el adhesivo, lo que permite levantarla posteriormente. El adhesivo entre la losa y el suelo se suavizará lo suficiente si se calienta por encima de 50 ºC durante al menos 2 minutos, pero su temperatura no debe superar 120 ºC para evitar su deterioro. Suponga que la losa y el suelo tienen una temperatura inicial de 25 ºC y propiedades termofísicas equivalentes de k = 0,15 6 3 W/m·K; ρ ·c p = 1,5·10 J/m ·K. a) ¿Cuánto tiempo se tardará en despegar una losa de espesor 4 mm empleando el bloque de hierro? ¿La temperatura del adhesivo excederá 120 ºC? b) Si el bloque de hierro tiene un área superficial cuadrada de 254 mm de lado, ¿cuánta energía se eliminará de él durante el tiempo que se tarda en despegar la losa? Solución: a) t = 168,7 s; la temperatura es menor que 120 ºC; b) Q = 56.063,6 J.
Conducción multidimensional 9.
(5.75 del Incropera) Una punta cilíndrica de cobre de 100 mm de longitud y 50 mm de diámetro está inicialmente a una temperatura uniforme de 20 ºC. Las caras de los extremos se someten de pronto a una intensa rapidez de calentamiento que las eleva a una temperatura de 500 ºC. Al mismo tiempo, la superficie cilíndrica se somete a calentamiento por un flujo de gas con una temperatura de 500 ºC y un coeficiente de transferencia de calor de 100 W/m2·K. Determine la temperatura en el punto central del cilindro 8 segundos después de la aplicación súbita del calor. 24
Transferencia de Calor / Curso 2009-10
Problemas
Resolución:
En el problema en cuestión la longitud y el diámetro del cilindro son comparables por lo que se tiene una transmisión de calor bidimensional. La expresión de la ecuación de calor para este caso será: 1 ∂ ∂T ∂ 2T 1 ∂T = r + r ∂r ∂r ∂ x 2 α ∂t
Se puede demostrar que la solución ha esta ecuación se puede obtener por el método de separación de variables, llegando a la siguiente expresión: T ( x, r , t ) − T ∞ T i − T ∞
=
T ( x, t ) − T ∞ T i − T ∞
· Pared plana
T ( r , t ) − T ∞ T i − T ∞
Cilindro infinito
Es decir, la solución bidimensional se expresa como producto de las soluciones unidimensionales correspondientes a una pared plana y a un cilindro infinito. Para las soluciones unidimensionales se emplearán el método de la resistencia interna despreciable (si Bi < 0,1), la solución exacta o la aproximada con el primer término (si Fo > 0,2) según corresponda. En la Tabla 4.2 del Cuaderno de Fórmulas, Tablas y Figuras se presenta un resumen de soluciones para distintos sistemas multidimensionales como productos de las soluciones unidimensionales. Visto esto se resuelve este problema en concreto. Como se pide la temperatura central en t = 8 s, a partir de la ecuación vista antes se tiene que: T (0,0,8) − T ∞ T i − T ∞
=
T (0,8) − T ∞ T i − T ∞
· Pared plana
T (0,8) − T ∞ T i − T ∞
Cilindro infinito
A continuación se resuelve cada sistema unidimensional por separado. Se calcula el número de Biot para la conducción a través de una pared plana: Bi pp =
hLc k
=∞ 25
Transferencia de Calor / Curso 2009-10
Problemas Como para este caso la temperatura superficial está fijada esto es equivalente a tener una h infinita. Como el número de Biot es mayor que 0,1 no se puede emplear el método de la resistencia interna despreciable. Empleamos la solución aproximada con el primer término. A partir de la Tabla 4.1 del Cuaderno de Fórmulas, Tablas y Figuras se obtienen los coeficientes de la aproximación de un término: ξ 1 = 1,5707 y C 1 = 1,2733. Se
buscan las propiedades del cobre para una temperatura media 20 + 500 T = = 260 º C = 533 K ⇒ ρ = 8.933 kg/m3; k = 384 W/m·K; c p = 410 J/kg·K. 2
Se calcula el número de Fourier: Fo =
α t L2c
=
kt ρ c p L2c
=
384·8 8.933·410·0,05 2
de:
= 0,304 . Como Fo
> 0,2 la aproximación con el primer término es correcta. T (0,8) − T ∞ T i − T ∞
Pared plana
= C 1 exp(−ξ 12 Fo) = 1,2733·exp(−1,5707 2 ·0,304) = 0,601
A continuación se resuelve la conducción a través de un cilindro infinito. Se calcula el número de Biot: Bicil . =
hr o k
=
100·0,025 384
= 6,51·10 −3
Por ser menor que 0,1 se puede aplicar el método de la resistencia interna despreciable: T (0,8) − T ∞ T i − T ∞ τ = Rt C c =
Cilindro infinito
ρ Vc p hA s
T (0,8) − T ∞ T i − T ∞
Cilindro infinito
T (0,0,8) − T ∞ T i − T ∞
=
t = exp − τ
=
8.933·0,025·410 100·2
= 457,8 s
8 = exp − = 0,9827 457 , 8
T (0,0,8) − 500 20 − 500
=
T (0,8) − T ∞ T i − T ∞
· Pared plana
T (0,8) − T ∞ T i − T ∞
= 0,601·0,9827 = 0,5906 Cilindro infinito
T (0,0,8) = 216,5 ºC Lo correcto sería volver a calcular las propiedades del cobre para esta temperatura y repetir el problema. Método de las diferencias finitas 10.
(5.82 del Incropera) Un cilindro de material plástico (α = 6·10 -7 m2/s) está inicialmente a una temperatura uniforme de 20 ºC y está bien aislado a lo largo de su superficie lateral y en un extremo. En el tiempo t = 0 se le aplica calor en el extremo izquierdo de manera que T 0 aumenta linealmente con el tiempo a una razón de 1 ºC/s.
26
Transferencia de Calor / Curso 2009-10
Problemas
a) Con el método explícito obtenga las ecuaciones en diferencias finitas para los nodos 1, 2, 3 y 4. b) Con Fo = ½ construya una tabla que tenga como encabezados p, t y las temperaturas nodales de T 0 a T 4. Determine la temperatura de la superficie, T 0, cuando T 4 = 35 ºC. c) Resuelva el problema mediante el método implícito y el método de la inversión de matrices Solución: b) ∆t = 30 s; T 0 (t = 210 s) = 230 ºC. Resolución del apartado c): De las ecuaciones nodales explícitas se pueden deducir las ecuaciones nodales implícitas fácilmente: Nodo
1:
T 0 p +1 + T 2 p +1 − 2T 1 p +1 = Nodo 2: FoT 1 p
+1
1 Fo
(T +
p 1 1
− T 1 p ) ⇒ −(1 + 2 Fo)T 1 p +1 + FoT 2 p +1 = − FoT 0 p+1 − T 1 p
− (1 + 2 Fo)T 2 p+1 + FoT 3 p +1 = −T 2 p
Nodo 3: FoT 2 p +1 − (1 + 2 Fo)T 3 p +1 + FoT 4 p +1 = −T 3 p Nodo 4: 2 FoT 3 p
+1
− (1 + 2 Fo)T 4 p+1 = −T 4 p
Expresadas en forma matricial: Fo 0 0 − (1 + 2 Fo) T 1 p+1 − FoT 0 p+1 − T 1 p Fo p+1 p − (1 + 2 Fo) Fo 0 − T T 2 2 = + p 1 p T 3 − (1 + 2 Fo) 0 Fo Fo − T 3 − (1 + 2 Fo) T 4 p+1 0 0 2 Fo − T 4 p
[ A][· T ] = [C ] ⇒ [T ] = [ A]−1 ·[C ] Se inicia un proceso iterativo en que para cada tiempo p hay que evaluar el vector [C] antes de resolver las temperaturas en el tiempo siguiente p+1. Los resultados se pueden expresar en forma de tabla. Con Fo = 0,5 ⇒ ∆t = 30 s.
27
Transferencia de Calor / Curso 2009-10
Problemas
11.
p
t
T 0
T 1
T 2
T 3
T 4
0
0
20
20
20
20
20
1
30
50
2
60
80
3
90
110
(5.100 del Incropera) Se sueldan dos barras muy largas en la dirección normal a la página, las cuales tienen las distribuciones de temperaturas iniciales que se muestran en la tabla inferior. En el tiempo t = 0, la cara m = 3 de la barra de cobre hace contacto con la cara m = 4 de la barra de acero AISI 1010. La soldadura actúa como una capa interfacial de espesor insignificante y resistencia efectiva de contacto Rt ′′,c = 2·10 -5 m2·K/W.
a) Obtenga la ecuación en diferencias finitas explícita en términos de Fo y Bic =
∆ x , k · Rt ′′,c
para T 4,2 y determine el criterio de estabilidad correspondiente. b) Si Fo = 0,01, determine T 4,2 un intervalo de tiempo después de que se hace contacto. ¿Cuál es el ∆t ? ¿Se satisface el criterio de estabilidad? Datos: Propiedades termofísicas del acero AISI 1010 a 1.000 K: k = 31,3 W/m·K; ρ = 7.832 kg/m3; c p = 1.168 J/kg·K. Solución: a) Criterio de estabilidad: Fo ≤ 1/(4 + 2 Bic); b) ∆t = 1,17 s; T 4,2 (t = 1,17 s) = 806,3 K; sí se satisface el criterio de estabilidad. 28
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 5. INTRODUCCIÓN A LA CONVECCIÓN Coeficientes de transferencia de calor 1.
(6.2 del Incropera) Para la convección laminar libre de una superficie vertical caliente, el coeficiente de convección local se expresa como h x = Cx-1/4, donde h x es el coeficiente a la distancia x desde el inicio de la superficie y la cantidad C , que depende de las propiedades del fluido, es independiente de x. Obtenga una expresión para la razón h x / h x, donde h x es el coeficiente promedio entre el inicio ( x = 0) y la posición x. Dibuje la variación de h x y h x con x. Solución: h x / h x = 4/3.
2.
(6.5 del Incropera) Aire a una temperatura de flujo libre T ∞ = 20 ºC está en un flujo paralelo sobre una placa plana de longitud L = 5 m y temperatura T s = 90 ºC. Sin embargo, los obstáculos colocados en el flujo intensifican la mezcla al aumentar la distancia x desde el inicio, y la variación espacial de las temperaturas medidas en la capa límite están correlacionadas por una expresión de la forma T ( x, y) [ºC] = 20 + 70 exp(-600 xy), donde x e y están en metros. Determine y elabore una gráfica de la forma en que varía el coeficiente de convección local h con x. Evalúe el coeficiente de convección promedio h para la placa. Solución: h x = 600·k ·x [W/m2·K] = 15,44· x[W/m2·K]; h = 38,6 W/m 2·K. Perfiles de la capa límite
3.
(6.10 del Incropera) Agua a una temperatura T ∞ = 25 ºC fluye sobre una de las superficies de una pared de acero (AISI 1010) cuya temperatura es T s1 = 40 ºC. La pared es de 0,35 m de espesor y la temperatura de la otra superficie es T s2 = 100 ºC. Para condiciones de estado estacionario, ¿cuál es el coeficiente de convección asociado con el flujo de agua? ¿Cuál es el gradiente de temperatura en la pared y en el agua que está en contacto con la pared? Dibuje la distribución de temperaturas en la pared y en el agua contigua. Datos: Propiedades termofísicas del acero AISI 1010: ρ = 7.832 kg/m 3. T c p k (K) (J/kg·K) (W/m·K) 300 434 63,9 400 487 58,7 ∂T agua ∂T pared Solución: h = 700 W/m2·K; = -171,4 K/m; ∂ y ∂ y
= -17.222,22 K/m. y = 0
Transición de la capa límite 4.
(6.12 del Incropera) Considere un flujo de aire sobre una placa plana de longitud L = 1 m en condiciones para las que ocurre la transición en xc = 0,5 m con base en el número de 29
Transferencia de Calor / Curso 2009-10
Problemas Reynolds crítico Re x,c = 5·105. En las regiones laminar y turbulenta los coeficientes de convección local son, respectivamente: hlam ( x ) = C lam x
−0 ,5
hturb ( x) = C turb x
y
−0 , 2
donde C lam = 8,845 W/m3/2·K, C turb = 49,75 W/m1,8·K y x tiene unidades de m. a) Mediante la evaluación de las propiedades termofísicas del aire a 350 K determine la velocidad del flujo de aire. b) Desarrolle una expresión para el coeficiente de convección promedio hlam ( x ) , como función de la distancia desde el inicio de la placa, x, para la región laminar, 0 ≤ x ≤ xc. c) Desarrolle una expresión para el coeficiente de convección promedio de toda la placa h ( x) , como función de la distancia desde el inicio de la placa, x, para la región turbulenta, xc≤ x ≤ L. d) Dibuje una gráfica de hlam ( x) , hlam ( x ) , hturb ( x) y h ( x) para 0 ≤ x ≤ L. Solución: h ( x ) =
1 x
a)
[2C
V
=
u∞
=
20,92
m/s;
b)
hlam ( x ) = 2C lam x
−0 ,5
;
c)
]
x c0,5 + 1,25C turb x 0,8 − 1,25C turb x c0,8 .
lam
Ecuaciones de conservación y soluciones: flujo de Couette 5.
(6.18 del Incropera 5ª edición en inglés) Considere un cojinete hidráulico ligeramente cargado que usa aceite con propiedades constantes ρ = 800 kg/m 3, ν = 10 -5 m2/s y k = 0,13 W/m·K. El diámetro del eje es de 75 mm, el espacio hasta su apoyo es de 0,25 mm y trabaja a 3.600 rpm. a) Determine la distribución de temperaturas en la película de aceite suponiendo que no hay transferencia de calor hacia el interior del eje y que la superficie del apoyo del cojinete se mantiene a 75 ºC. b) ¿Cuál es la transferencia de calor del cojinete? 2
Solución: a) T ( y ) = T 0 − 6.
ρν ∂u y 2
− Ly ; b) q’( y = 0) = -1.508 W/m. K ∂ y 2
(6.21 del Incropera) Considere el flujo de Couette con transferencia de calor para el cual la placa inferior (placa móvil) se mueve con una velocidad de U = 5 m/s y está perfectamente aislada. La placa superior (placa estacionaria) está construida con un material de conductividad térmica k pe = 1,5 W/m·K y espesor L pe = 3 mm. La superficie externa se mantiene a T pe = 40 ºC. Las placas están separadas por una distancia Lo = 5 mm, que se llena con un aceite de motor de viscosidad µ o = 0,799 N·s/m 2 y conductividad térmica k o = 0,145 W/m·K.
30
Transferencia de Calor / Curso 2009-10
Problemas
a) Determine la distribución de temperaturas en la película de aceite, T o( y). b) Determine la temperaturas en los extremos de la película de aceite en contacto con las dos placas, T o( y = 0) y T o( y = Lo). c) Dibuje la distribución de temperaturas en la película de aceite y en la placa estacionaria. Solución: a) T ( y) = -2,755·10 6 ·y2 + 390 [K]; b) T o( y = 0) = 390 K y T o( y = Lo) = 321 K. 7.
(6.22 del Incropera) Un eje con un diámetro de 100 mm gira a 9.000 rpm en un cojinete hidráulico de 700 mm de longitud. Las propiedades del aceite lubricante son µ = 0,03 N·s/m2 y k = 0,15 W/m·K, mientras que el material del cojinete tiene una conductividad térmica de k c = 45 W/m·K.
a) Determine la disipación viscosa µ Φ [W/m3] en el lubricante. b) Determine la transferencia de calor [W] del lubricante suponiendo que no se pierde calor a través del eje. c) Si la cubierta del cojinete se enfría con agua, de modo que la superficie externa del cojinete se mantiene a 30 ºC, determine las temperaturas del cojinete y del eje, T c y T e. Solución: a) µ Φ = 6,66·10 7·W/m3; b) q = 14.646 W; c) T c = 81,3 ºC y T e = 303,3 ºC. Similitud y parámetros adimensionales 8.
(6.27 del Incropera, examen septiembre 2002) Un objeto de forma irregular tiene una longitud característica L1 = 1 m y se mantiene a una temperatura superficial uniforme T s,1 = 400 K. Cuando se coloca en aire atmosférico a una temperatura T ∞ = 300 K y se mueve con una velocidad V 1 = 100 m/s, el flujo promedio de calor desde la superficie al aire es 20.000 W/m2. Si un segundo objeto de la misma forma, pero con una longitud característica, L2 = 5 m, se mantiene a la misma temperatura superficial y se coloca en aire atmosférico a T ∞ = 31
Transferencia de Calor / Curso 2009-10
Problemas 300 K, ¿cuál será el valor del coeficiente promedio de convección si la velocidad del aire es V 2 = 20 m/s? 9.
(6.31 y 6.34 del Incropera) Considere las condiciones para las que un fluido con una velocidad de flujo libre V = 1 m/s fluye sobre una superficie con una longitud característica L = 1 m, lo que proporciona un coeficiente promedio de transferencia de calor por convección h = 100 W/m2·K. a) Calcule los parámetros adimensionales Nu L , Re L y Pr para los siguientes fluidos: aire, aceite de motor, mercurio y agua. Suponga que los fluidos están a una temperatura de 300 K. b) Dibuje la variación relativa entre el espesor de la capa límite hidrodinámica y la capa límite térmica para el flujo laminar sobre una placa plana para los cuatro fluidos anteriores. Suponga un coeficiente del número de Prandtl n = 0,35. Datos: Propiedades termofísicas del mercurio: k ·103 T c p Pr µ ·102 ν ·106 α ·107 ρ 3 2 2 2 (K) (kg/m ) (kJ/kg·K) (N·s/m ) (m /s) (W/m·K) (m /s) 300 13.529 0,1393 0,1523 0,1125 8.540 45,30 0,0248 Solución: a) aire: Nu L = 3.802; Re L = 62.932; Pr = 0,706; aceite de motor: Nu L = 690; Re L = 1.818; Pr = 6.403; mercurio: Nu L = 11,7; Re L = 8.888.889; Pr = 0,0248; agua: Nu L = 163; Re L = 1.166.092; Pr = 5,83; b) δ δ t
10.
δ δ t
= 0,885 ; aire
δ δ t
= 21,5 ; aceite
δ δ t
= 0,27 ; mercurio
= 1,85 . agua
(6.35 del Incropera) Se utiliza aire forzado a T ∞ = 27 ºC y V = 10 m/s para enfriar elementos electrónicos sobre una tarjeta de circuitos. Uno de tales elementos es un chip, de 4 mm por 4 mm, que se localiza a 120 mm desde el inicio de la tarjeta. Los experimentos revelan que el flujo sobre ésta es perturbado por los elementos y que la transferencia de calor por convección está correlacionada mediante una expresión de la forma: Nu x = 0,04 Re x0,85 Pr 1 / 3 Estime la temperatura superficial del chip si éste disipa 30 mW. Nota: Las propiedades del fluido se deben calcular a la temperatura media entre el fluido y la superficie.
Solución: T = 317 K = 44 ºC. 32
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 6. CONVECCIÓN FORZADA EN FLUJO EXTERNO Placa plana en flujo paralelo 1.
(7.3 del Incropera) Sobre ambas superficies de una placa plana de 1 m de longitud que se mantiene a 20 ºC fluye aceite de motor a 100 ºC y a una velocidad de 0,1 m/s. Determine: a) Los espesores de las capas límite de velocidad y térmica al final de la placa. b) El flujo local de calor al final de la placa. c) La transferencia de calor por unidad de ancho de la placa.
′′ = 1.304 W/m2; c) qtotal ′ = 5.216 W/m. Solución: Solución: a) δ L = 146,7 mm; δ t,L t,L = 14,3 mm; b) q L 2.
(7.11 del Incropera) Una placa plana de 1 m de ancho se mantiene a una temperatura superficial uniforme T s = 150 ºC mediante el uso de módulos rectangulares generadores de calor, controlados de manera independiente, de espesor a = 10 mm y longitud b = 50 mm. Cada módulo está aislado de sus alrededores, así como de su parte posterior. Aire atmosférico a 25 ºC fluye sobre la placa a una velocidad de 30 m/s. Las propiedades termofísicas del módulo son k = = 5,2 W/m·K, c p = 320 J/kg·K y ρ = = 2.300 kg/m 3.
a) Encuentre la generación de potencia que se requiere, q& [W/m3], en el módulo situado a 700 mm del inicio. b) Encuentre la temperatura máxima en ese módulo generador de calor. Solución: Solución: a) q& = 873.448 W/m 3; b) T máx máx = 158,4 ºC. 3.
(7.32 del Incropera) Un conjunto de componentes electrónicos disipadores de calor se monta sobre el lado inferior de una placa horizontal de aluminio de 1,2 x 1,2 m, mientras que el lado superior se enfría con un flujo de aire para el que u∞ = 15 m/s y T ∞ = 300 K. La placa se une a un recinto bien aislado de manera que todo el calor se debe transferir al aire. a ire. Además la placa de aluminio es lo suficientemente delgada para asegurar una temperatura casi uniforme en ella.
33
Transferencia de Calor / Curso 2009-10
Problemas
a) Si la temperatura del aluminio no debe exceder de 350 K, ¿cuál es la disipación de calor máxima permisible? b) Determine la disipación de calor máxima permisible si sobre la placa de aluminio se coloca un sistema de aletas longitudinales cuya longitud, espesor y espaciado (entre planos centrales de las aletas) son 25 mm, 5 mm y 10 mm, respectivamente. Datos: Propiedades termofísicas del aluminio: T c p k ρ 3 (K) (kg/m ) (J/kg·K) (W/m·K) 300 2.702 903 237 Solución: Solución: a) q = 2,14 kW; b) q = 12,71 kW. Flujo alrededor de un cilindro 4.
(7.42 del Incropera) Un tubo circular de 25 mm de diámetro exterior se coloca en un flujo de aire a 25 ºC y presión de 1 atm. El aire se mueve en flujo cruzado sobre el tubo a 15 m/s, mientras que la superficie externa del tubo se mantiene a 100 ºC. ¿Cuál es la transferencia de calor del tubo por unidad de longitud? Solución: Solución: q ′ = 527 W/m.
5.
(7.49 del Incropera) Aire a 27 ºC y a velocidad de 5 m/s pasa sobre la pequeña región A s (20 x 20 mm) en una superficie grande que se mantiene a T s = 127 ºC. Para estas condiciones se eliminan 0,5 W de la superficie A s. Con el fin de aumentar la velocidad de eliminación de calor, se sujeta una aleta de alfiler de acero inoxidable AISI 304 de diámetro 5 mm a A a A s, que se supone permanece a T s = 127ºC. Datos: Propiedades termofísicas del acero inoxidable AISI 304: T c p k ρ 3 (K) (kg/m ) (J/kg·K) (W/m·K) 300 7.900 477 14,9 400 7.900 515 16,6
34
Transferencia de Calor / Curso 2009-10
Problemas
a) Determine la velocidad máxima posible de eliminación de calor a través de la aleta. b) ¿Qué longitud de aleta proporciona esa velocidad máxima de eliminación de calor? c) Determine la efectividad de la aleta, ε f . d) ¿Cuál es el porcentaje de aumento en la transferencia de calor de A s debido a la instalación de la aleta? Solución: Solución: a) q fmáx = 2,3 W; b) L b) L = = 37,4 mm; c) ε f = 93,7; d) 455 %. 6.
(7.62 del Incropera) En un proceso de fabricación, una varilla larga recubierta de plástico ( ρ = = 2.200 kg/m 3, c p = 800 J/kg·K, k = = 1 W/m·K) de diámetro D diámetro D = = 10 mm está inicialmente a una temperatura uniforme de 25 ºC y súbitamente se expone a un flujo cruzado de aire a u∞ = 50 m/s y T ∞ = 350 ºC. ¿Cuánto tiempo le costará a la superficie de la varilla alcanzar 175 ºC, temperatura por arriba de la cual se curará el recubrimiento especial? Solución: Solución: t = = 8,27 s (Calculado con T f = ((25+175)/2 +350) /2 y correlación de Hilpert ⇒ ·K); Fo = = 0,1879. h = 240 W/m2·K); Fo Flujo a través de bancos de tubos
7.
(7.76 del Incropera) Considere un banco de tubos escalonado para el que el diámetro exterior del tubo es 16,4 mm y los espaciados longitudinal y transversal son S L = S T T = 20,5 mm. Hay siete líneas de tubos en la dirección del flujo de aire y ocho tubos por línea ( N ( N L = 7 y N T T = 8). La temperatura superficial de los tubos es de 70 ºC y la temperatura y velocidad del flujo de aire son 15 ºC y 6 m/s, respectivamente. Determine el coeficiente de convección y la transferencia de calor para el banco de tubos. Comience evaluando las propiedades del aire a T ∞ = 15 ºC. Solución: Solución: h = 235 W/m2·K; q’ = 28,5 kW/m.
8.
(7.80 del Incropera) Los componentes eléctricos montados en cada una de dos placas isotérmicas se enfrían al hacer pasar aire atmosférico entre ellas, y se usa un dispositivo en línea de aletas de alfiler de aluminio para aumentar la transferencia de calor al aire. Las aletas son de diámetro D diámetro D = = 2 mm, longitud L = L = 100 mm y conductividad térmica k = = 240 W/m·K. Los espaciados longitudinales y transversales son S L = S T T = 4 mm con un dispositivo cuadrado de 625 aletas ( N L = N = N T 25) montadas en placas cuadradas de ancho T = W = = 100 mm. El aire entra a una velocidad de 10 m/s y a una temperatura de 300 K.
35
Transferencia de Calor / Curso 2009-10
Problemas
a) Evaluando las propiedades del aire a 300 K, estime el coeficiente promedio de convección para el dispositivo de aletas. b) Suponiendo el coeficiente de convección anterior uniforme sobre las aletas y las placas determine la temperatura de salida del aire y la transferencia total de calor cuando las placas se mantienen a 350 K. Sugerencia: La temperatura de salida del aire está gobernada por una relación exponencial de la forma [( T placa – T sal ) / (T placa – T ent )] = & c p)], donde m& = ρ VLN TS exp[-( h At η o) / ( m T es el flujo de masa del aire que pasa a través de las aletas, At es el área superficial total de transferencia de calor (placas más aletas) y η o es la eficiencia superficial global del conjunto placas más aletas. Solución: a) h = 435 W/m 2·K; b) T sal = 321 K; qtotal = 2.467,3 W.
36
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 7. CONVECCIÓN FORZADA EN FLUJO INTERNO Consideraciones térmicas y balance de energía 1.
(8.12 del Incropera) Considere una varilla cilíndrica de combustible nuclear de longitud L y diámetro D que está encerrada en un tubo concéntrico. Por la región anular entre la & y la superficie exterior del varilla y el tubo fluye agua presurizada a un flujo másico m tubo está bien aislada. Hay generación de calor dentro de la varilla de combustible y se sabe que la rapidez de generación volumétrica varía senoidalmente con la distancia a lo largo de la varilla: q& ( x) = q& o sen(π x / L ) , siendo q& o [W/m3] una constante. Se puede suponer que existe un coeficiente de convección uniforme entre la superficie de la varilla y el agua.
a) Obtenga expresiones para el flujo de calor local q ′′( x) y para la transferencia total de calor q de la varilla de combustible al agua. b) Obtenga una expresión para la variación de la temperatura media del agua a lo largo del tubo, T m( x). c) Obtenga una expresión para la variación de la temperatura superficial de la varilla a lo largo del tubo, T s( x). Desarrolle una expresión para la posición x para la que se maximiza. D 2 Lq& o D 2 Lq& o π x Solución: a) q ′′( x) = ; b) T m ( x) = T m,ent + sen ; q = & c p 4 2 4m L m& c p q ′′( x) L . c) T s ( x ) = T m ( x ) + ; x máx = arctan − h π DLh Dq& o
2.
π x 1 − cos L ;
(8.15 del Incropera) Aire atmosférico entra en la sección caliente de un tubo circular con un flujo de 0,005 kg/s y una temperatura de 20 ºC. El tubo es de diámetro D = 50 mm y existen condiciones completamente desarrolladas con h = 25 W/m2·K sobre toda la longitud L = 3 m. a) Para el caso de un flujo de calor superficial uniforme, q s′′ = 1.000 W/m 2, determine la transferencia total de calor, q, y la temperatura media del aire que sale del tubo, T m,sal . 37
Transferencia de Calor / Curso 2009-10
Problemas ¿Cuál son los valores de las temperaturas superficiales del tubo en la entrada y en la salida, T s,ent y T s,sal ? Dibuje la variación con x de T m y T s. b) Si el flujo de calor superficial varía linealmente con x, q s′′( x) = 500 x [W/m2], ¿cuáles son los valores de q, T m,sal , T s,ent y T s,sal ? Dibuje la variación con x de T m y T s. c) Para cada tipo de proceso de calentamiento anterior, ¿qué flujos de calor se requieren para alcanzar en la salida una temperatura del aire de T m,sal = 125 ºC? Solución: a) q = 471,24 W; T m,sal = 113,6 ºC; T s,sal = 153,6 ºC; T s,ent = 60 ºC; b) q = 353,4 W; T m,sal = 90,2 ºC; T s,sal = 150,2 ºC; T s,ent = 20 ºC; c) q s′′ = 1.122 W/m2; q s′′( x ) = 748 x [W/m2]. Correlaciones de transferencia en tubos circulares 3.
(8.19 del Incropera) Aceite de motor a razón de 0,02 kg/s fluye por un tubo de 3 mm de diámetro y de 30 m de longitud. El aceite tiene una temperatura de entrada de 60 ºC, mientras que la temperatura de la pared del tubo se mantiene a 100 ºC por condensación de vapor sobre su superficie externa. a) Estime el coeficiente promedio de transferencia de calor para el flujo interno del aceite. b) Determine la temperatura de salida del aceite. c) Calcule la transferencia total de calor al aceite. Solución: a) h = 201 W/m 2·K; b) T m,sal = 89,6 ºC; c) q = 1.249 W.
4.
(8.24 del Incropera) En las etapas finales de producción se esteriliza un fármaco calentándolo de 25 a 75 ºC a medida que se mueve a 0,2 m/s por un tubo recto de acero inoxidable de pared delgada de 12,7 mm de diámetro. Un flujo de calor uniforme se mantiene mediante un calentador de resistencia eléctrica enrollado alrededor de la superficie externa del tubo. Si el tubo es de 10 m de longitud, ¿cuál es el flujo de calor que se requiere? Si entra fluido al tubo con un perfil de velocidad completamente desarrollado y un perfil uniforme de temperatura, ¿cuál es la temperatura superficial en la salida del tubo y a una distancia de 0,5 m desde la entrada? Las propiedades del fluido se pueden aproximar a: ρ = 1.000 kg/m 3; c p = 4.000 J/kg·K; µ = 2·10 -3 kg/s·m; k = 0,48 W/m·K; Pr = 10. Solución: q s′′ =12.700 W/m2; T s,sal = 152 ºC; T s ( x = 0,5) = 65,7 ºC.
5.
(8.39 del Incropera) El núcleo de un reactor nuclear de alta temperatura enfriado por gas tiene tubos de fluido refrigerante de 20 mm de diámetro y 1.500 mm de longitud. Entra helio a 600 K y sale a 1.000 K cuando el flujo es de 8·10 -3 kg/s por tubo. a) Determine la temperatura superficial uniforme de la pared del tubo para estas condiciones. b) Si el gas refrigerante es aire en vez de helio, determine el flujo másico y la temperatura de salida del aire si la rapidez de eliminación de calor y la temperatura superficial del tubo permanecen iguales.
38
Transferencia de Calor / Curso 2009-10
Problemas Datos: Propiedades termofísicas del helio: k ·103 T c p Pr µ ·107 ν ·106 α ·106 ρ 3 2 2 2 (K) (kg/m ) (J/kg·K) (N·s/m ) (m /s) (W/m·K) (m /s) 800 0,05924 5193 382 304 0,654 & aire ≈ 0,0474 kg/s. Solución: a) T s = 1.142,8 K; b) T m,sal ≈ 919 K; m 6.
(8.50 del Incropera) Considere un tubo de pared delgada de 10 mm de diámetro y 2 m de & = 0,2 kg/s y T m,ent = 47 ºC. longitud. Entra agua al tubo desde un recipiente grande a m a) Si la superficie del tubo se mantiene a una temperatura uniforme de 27 ºC, ¿cuál es la temperatura de salida del agua, T m,sal ? b) ¿Cuál es la temperatura de salida del agua si se calienta mediante el paso de aire a T ∞ = 100 ºC y V = 10 m/s en flujo cruzado sobre el tubo? Las propiedades del aire se pueden evaluar a una temperatura de película supuesta de 350 K. Solución: a) T m,sal ≈ 35-36 ºC; b) T m,sal ≈ 47 ºC.
7.
(8.53 del Incropera) Una tubería de acero (k = 60 W/m·K) que conduce agua caliente se enfría externamente mediante aire en flujo cruzado a una velocidad de 20 m/s y una temperatura de 25 ºC. Los diámetros interno y externo de la tubería son Dint = 20 mm y Dext = 25 mm, respectivamente. En cierta posición a lo largo de la tubería, la temperatura media del agua es 80 ºC. Suponga que el flujo dentro del tubo está completamente desarrollado con un número de Reynolds de 20.000. Encuentre la transferencia de calor al flujo de aire por unidad de longitud de tubería. Solución: q ′ = 490 W/m. Correlaciones en conductos no circulares
8.
(8.72 del Incropera) Aire a 3·10 -4 kg/s y 27 ºC entra en un conducto de sección rectangular de 1 m de longitud y con una sección de 4 por 16 mm. Se impone un flujo de calor uniforme de 600 W/m 2 sobre la superficie del conducto. ¿Cuáles son las temperaturas del aire y la de la superficie del conducto en la salida? Solución: T m,sal = 106,2 ºC; T s,sal ≈ 128,5 ºC.
9.
(8.80 del Incropera) Una placa fría es un dispositivo de enfriamiento activo que se une a un sistema generador de calor con el fin de disipar el calor mientras se mantiene el sistema a una temperatura aceptable. Normalmente se fabrica de un material de alta conductividad térmica, k pf , dentro del que se hacen canales por los que pasa el fluido refrigerante. Considere la placa fría de las dimensiones de la figura en la que las paredes laterales se pueden suponer aisladas y las paredes superior e inferior se mantienen en contacto con el dispositivo generador de calor a una temperatura constante T s. La velocidad y la temperatura media del refrigerante son um y T m,ent , respectivamente.
39
Transferencia de Calor / Curso 2009-10
Problemas
a) Suponiendo flujo turbulento completamente desarrollado a través de cada canal plantee todas las ecuaciones necesarias para calcular la temperatura de salida del refrigerante, T m,sal , y la transferencia total de calor a la placa fría, q. b) Resuelva el sistema con los siguientes datos: W = 100 mm, H = 10 mm, 10 canales cuadrados de w = h = 6 mm, δ = 4 mm. T s = 360 K y k pf = 400 W/m·K. El refrigerante es agua a um = 2 m/s y T m,ent = 300 K. Propiedades promedio del agua: ρ = 984 kg/m3, c p = 4.184 J/kg·K, µ = 489·10-6 N·s/m2, k = 0,65 W/m·K, Pr = 3,15. Solución: T m,sal = 305,3 K; q = 15.818 W.
40
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 8. CONVECCIÓN LIBRE Placas verticales 1.
(9.12 del Incropera) La puerta de un horno de 0,5 m de altura y 0,7 m de ancho alcanza una temperatura superficial promedio de 32 ºC durante el funcionamiento del horno. Estime la pérdida de calor al cuarto con aire ambiente a 22 ºC. Si la puerta tiene una emisividad de ε = 1 y los alrededores también están a 22 ºC, comente la pérdida de calor por convección libre con relación a la de radiación. Solución: qconv = 9,9 W; qrad = 21,4 W; qconv / qrad = 46 %.
2.
(9.17 del Incropera) El vidrio delgado de una ventana cuadrada de 1 m de lado separa aire quieto de una habitación a T ∞int = 20 ºC del aire ambiental exterior a T ∞ext = -20 ºC. Las paredes de la habitación y el medio externo están a la misma temperatura (T alr, int = 20 ºC, T alr, ext = -20 ºC). Si el vidrio tiene una emisividad de ε = 1, ¿cuál es su temperatura T ? ¿Cuál es la pérdida de calor a través del vidrio?
Solución: T ≈ 1 ºC; q = 174 W. 3.
(9.20 del Incropera) Un contenedor de paredes delgadas con agua caliente a 50 ºC se coloca en un baño de agua fría en reposo a 10 ºC. La transferencia de calor en las superficies interna y externa del contenedor se pueden aproximar por la convección libre en una placa vertical.
41
Transferencia de Calor / Curso 2009-10
Problemas Determine el coeficiente global promedio de transferencia de calor, U , entre el agua caliente del contenedor y el agua fría del baño. Solución: U ≈ 350 W/m2·K. Placas inclinadas 4.
(9.42 del Incropera) Una placa cuadrada de 1 m de lado inclinada un ángulo de 45º se expone a un flujo neto de calor por radiación de 300 W/m 2 en su superficie inferior. Si la superficie superior de la placa está bien aislada, estime la temperatura que alcanza la placa cuando el aire ambiente está en reposo y a 0 ºC. Solución: T = 61 ºC. Placas horizontales
5.
(9.34 del Incropera) Una parrilla circular de 0,25 m de diámetro y emisividad de 0,9 se mantiene a una temperatura superficial constante de 130 ºC. ¿Qué potencia eléctrica se requiere cuando el aire ambiente y los alrededores están a 24 ºC? Suponga que la parte inferior de la parrilla está aislada. Solución: q = 91 W.
6.
(9.29 del Incropera) Se tiene un transformador de potencia eléctrica de forma cilíndrica de diámetro 300 mm y altura 500 mm. Se desea mantener su temperatura superficial a 47 ºC mediante refrigeración por convección libre y radiación. La superficie tiene una emisividad de ε = 0,80. a) Determine cuanta potencia se puede eliminar de sus superficies lateral y horizontal superior (la inferior se supone aislada) cuando la temperatura ambiente y los alrededores están a 27 ºC. b) Se añaden 30 aletas verticales de 5 mm de espesor y 75 mm de longitud a lo largo de toda la superficie lateral. ¿Cuál es ahora la rapidez de eliminación de calor si las aletas tienen la misma emisividad de 0,80 y una conductividad térmica de 240 W/m·K?
Solución: a) q = 104,5 W; b) q = 531 W. Cilindros horizontales 7.
(9.60 del Incropera) En un tubo de pared delgada de 20 mm de diámetro circula un fluido caliente a una temperatura media de 45 ºC. El tubo se monta horizontalmente en aire en reposo a 15 ºC. Se enrolla una cinta delgada de calentamiento eléctrico sobre la superficie externa del tubo para evitar pérdidas de calor del fluido caliente al aire ambiente y mantener su temperatura constante. 42
Transferencia de Calor / Curso 2009-10
Problemas
a) Ignorando las pérdidas de calor por radiación, calcule el flujo de calor qe′′ que se debe suministrar mediante la cinta eléctrica. b) Calcule el flujo de calor teniendo en cuenta la radiación. La emisividad de la cinta es 0,95 y los alrededores también están a 15 ºC. c) Calcule el flujo de calor si se añade una capa de aislante (k = 0,050 W/m·K) de espesor 20 mm y emisividad 0,60. ¿Cuál será la temperatura superficial del aislante? Solución: a) qe′′ = 208 W/m 2; b) q e′′ = 388 W/m 2; c) q e′′ = 113,65 W/m 2; T aisl = 20 ºC. 8.
(9.51 del Incropera) Vapor saturado a 4 bar de presión con una velocidad media de 3 m/s fluye a través de una tubería horizontal de acero inoxidable AISI 302 cuyos diámetros interior y exterior son 55 y 65 mm, respectivamente. Se sabe que el coeficiente de transferencia de calor por convección para el flujo de vapor es 11.000 W/m 2·K. a) Si la tubería se cubre con una capa de 25 mm de espesor de aislante ( k = 0,051 W/m·K) y se expone a aire atmosférico a 25 ºC, determine la transferencia de calor por convección libre al cuarto por unidad de longitud de tubería. ¿Cuál es la temperatura de la superficie exterior del aislante? b) La radiación neta a los alrededores también contribuye a la pérdida de calor de la tubería. Si el aislante tiene una emisividad ε = 0,80 y los alrededores están a 25 ºC, ¿cuál será ahora la pérdida de calor de la tubería por unidad de longitud? Calcule la temperatura de la superficie exterior del aislante en estas condiciones. Datos: Acero inoxidable AISI 302: k (a 300 K) = 15,1 W/m·K. Solución: a) q ′ ≈ 50 W/m; T aisl ≈ 52 ºC.
43
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 9. INTRODUCCIÓN A LA RADIACIÓN Intensidad, potencia emisiva e irradiación 1.
(12.11 del Incropera) Con el fin de iniciar la operación de un proceso se emplea un sensor de movimiento infrarrojo (detector de radiación) para determinar la aproximación de una pieza caliente sobre un sistema transportador. La señal de salida del sensor es proporcional a la rapidez a la que la radiación incide sobre el sensor.
Para Ld = 1 m, ¿en qué posición x1 la señal S 1 del sensor es el 75 % de la señal correspondiente a la posición directamente debajo del sensor, S o ( x = 0)? Solución: x1 = 0,39 m. Radiación de cuerpo negro 2.
(12.17 del Incropera) Suponiendo que la superficie de la Tierra es negra, estime su temperatura si el Sol tiene una temperatura equivalente de cuerpo negro de 5.800 K. Los diámetros del Sol y de la Tierra son 1,39·10 9 y 1,29·10 7 m, respectivamente, y la distancia entre ellos es 1,5·1011 m. Solución: T T = 279 K = 6 ºC.
3.
(12.24 del Incropera) Un elemento de calentamiento radiante eléctrico en forma de anillo se comporta como un cuerpo negro y se mantiene a una temperatura de T h = 3.000 K y se usa en un proceso de producción para calentar una pequeña pieza que tiene un área superficial A p = 0,007 m 2. La superficie del elemento de calentamiento se puede suponer negra. Para θ 1 = 30º, θ 2 = 60º, L = 3 m y W = 30 mm, ¿cuál es la rapidez a la que la energía radiante emitida por el calentador incide sobre la pieza?
45
Transferencia de Calor / Curso 2009-10
Problemas
Solución: q = 278,4 W. Emisividad, absortividad, reflectividad y transmisividad 4.
(12.25 del Incropera) La emisividad hemisférica espectral del tungsteno se puede aproximar mediante la distribución que se describe más adelante. Considere un filamento cilíndrico de tungsteno de diámetro D = 0,8 mm y longitud L = 20 mm. El filamento se encierra en un bulbo al vacío y se calienta mediante una corriente eléctrica a una temperatura de estado estable de 2.900 K.
a) ¿Cuál es la emisividad hemisférica total cuando la temperatura del filamento es 2.900 K? b) Suponga que los alrededores están a 300 K. ¿Cuál es la rapidez inicial de enfriamiento del filamento cuando se desconecta la corriente? Suponga que el tungsteno se comporta como una superficie difusa. Propiedades termofísicas del tungsteno: ρ = 19.300 kg/m 3; c p (T = 2.000 K) = 167 J/kg·K; c p (T = 2.500 K) = 176 J/kg·K. Solución: a) ε = 0,352; b) dT / dt = -1.996 K/s. 46
Transferencia de Calor / Curso 2009-10
Problemas 5.
(12.39 del Incropera) Un pequeño objeto difuso y opaco a T s = 400 K se suspende en un horno grande cuyas paredes interiores están a T f = 2.000 K. Las paredes son difusas y grises y tienen una emisividad de 0,20. La emisividad hemisférica espectral para la superficie del objeto se da a continuación.
a) Determine la emisividad y la absortividad total de la superficie del objeto. b) Evalúe el flujo radiante reflejado y el flujo neto radiativo hacia la superficie. c) ¿Cuál es la potencia emisiva espectral a λ = 2 µ m? d) ¿Cuál es la longitud de onda λ 1/2 para la que la mitad de la radiación total emitida por la superficie está en la región espectral λ > λ 1/2?
′ = 362.880 W/m 2; q rad ′′ , neto = 543.594 W/m 2; c) E λ (λ Solución: a) ε s = 0,5; α s = 0,6; b) q ′ref = 2 µ m) = 0,1265 W/m2·µ m; d) λ 1/2 = 10,26 µ m. 6.
(12.55 del Incropera) Considere una superficie difusa opaca cuya reflectividad espectral varía con la longitud de onda como se muestra. La superficie está a 750 K, y la irradiación sobre un lado varía con la longitud de onda como se muestra. El otro lado de la superficie está aislado. ¿Cuáles son la absortividad total y la emisividad de la superficie? ¿Cuál es el flujo neto de calor radiativo hacia la superficie?
′ , neto = -11.757 W/m 2. Solución: α = 0,72; ε = 0,756; q ′rad Aplicaciones 7.
(12.74 del Incropera) Un procedimiento para medir la conductividad térmica de sólidos a temperaturas elevadas implica la colocación de una muestra en la parte inferior de un horno. La muestra tiene espesor L y se coloca en un contenedor cuadrado de lado W cuyos lados están bien aislados. Las paredes de la cavidad se mantienen a T w, mientras que la superficie inferior de la muestra se mantiene a una temperatura mucho más baja T c al hacer circular un refrigerante a través del contenedor de la muestra. La superficie superior de la 47
Transferencia de Calor / Curso 2009-10
Problemas muestra es difusa y gris con una emisividad ε s. Su temperatura T s se mide de manera óptica.
a) Ignorando los efectos de la convección, obtenga una expresión de la cual se pueda evaluar la conductividad térmica de la muestra en términos de las cantidades medidas y conocidas (T w, T s, T c, ε s y L). Las mediciones se realizan bajo condiciones de estado estacionario. Si T w = 1.400 K, T s = 1.000 K, T c = 300 K, ε s = 0,85 y L = 0,015 m, ¿cuál es la conductividad térmica de la muestra?
& c = 0,1 kg/s, ¿es b) Si W = 0,1 m y el refrigerante es agua con un flujo másico de m razonable suponer una temperatura uniforme T c de la superficie inferior? Solución: a) k = 2,93 W/m·K; b) ∆T agua = 3,3 K ⇒ Sí es razonable suponer T c uniforme. 8.
(12.78 del Incropera) Un termógrafo es un dispositivo que responde a la potencia radiativa que incide sobre su detector en la región espectral 9-12 µ m. El termógrafo proporciona una imagen, por ejemplo del lado de un horno, de la que se puede determinar la temperatura superficial.
a) Para una superficie negra a 60 ºC determine la potencia emisiva para la región espectral 9-12 µ m. b) Calcule la potencia radiante (en W) recibida por el termógrafo en el mismo rango (9-12 µ m) cuando ve, en una dirección normal, una pequeña área de una pared negra, 200 48
Transferencia de Calor / Curso 2009-10
Problemas mm2, a T s = 60 ºC. El ángulo sólido ω subtendido por la abertura del termógrafo cuando se ve desde el área objetivo es 0,001 sr. c) Determine la potencia radiante (en W) que recibe el termógrafo para la misma área de pared y ángulo sólido anterior cuando la pared es un material gris, opaco y difuso a T s = 60 ºC, con emisividad 0,70 y los alrededores son negros a T alr = 23 ºC. Solución: a) E b (9-12) = 144,77 W/m 2; b) qrad = 9,216·10 -6 W; c) qrad = 8,08·10 -6 W.
49
Transferencia de Calor / Curso 2009-10
Problemas
PROBLEMAS TEMA 10. INTERCAMBIO DE RADIACIÓN ENTRE SUPERFICIES Factores de forma 1.
(13.1 del Incropera) Determine F 12 y F 21 para las siguientes configuraciones con el uso de las relaciones básicas de los factores de forma. No use tablas ni gráficas.
Solución: a) F 12 = 1 y F 21 = 4/3π; b) F 12 = 0,5 y F 21 = 0,25; c) F 12 = 1 y F 21 = 2/ π; d) F 12 = 0,5 y F 21 = 1/ 2 ; e) F 12 = 0,5 y F 21 = 0; f) F 12 = 1 y F 21 = 1/8; g) F 12 = 0,5 y F 21 = 2/π. 2.
(13.4 del Incropera) Una capa metálica hemisférica delgada de diámetro D = 0,8 m se suspende dentro de un recinto cúbico de 1,5 metros de lado. Determine los factores de forma F 11, F 22 y F 33. 51
Transferencia de Calor / Curso 2009-10
Problemas
Solución: F 11 = 0,5; F 22 = 0; F 33 = 0,888. 3.
(13.7 del Incropera) Determine el factor de forma F 12 de los rectángulos perpendiculares que se muestran en la figura.
Solución: F 12 = 0,075. Intercambio de radiación de cuerpo negro 4.
(13.15 del Incropera) Un calentador tubular con una superficie interior negra de temperatura uniforme T s = 1.000 K irradia un disco coaxial. Determine la potencia radiante del calentador que incide sobre el disco, q s→1, y la irradiación sobre el disco, G1.
Solución: q s→1 = 15,6 W; G1 =7.938 W/m2. 52
Transferencia de Calor / Curso 2009-10
Problemas 5.
(13.32 del Incropera) El agua que fluye a través de un número grande de tubos largos circulares de pared delgada se calienta por medio de placas paralelas por encima y por debajo del dispositivo de tubos. El espacio entre las placas está al vacío y las superficies de las placas y tubos se pueden aproximar como cuerpos negros. Ignorando las variaciones axiales, determine la temperatura superficial de los tubos, T s, si el agua fluye a través de & = 0,2 kg/s y a una temperatura media T m = 300 K. cada tubo con un caudal másico m
Solución: T s ≈ 308 K. Intercambio de radiación entre superficies grises difusas 6.
(13.39 del Incropera) Considere dos placas paralelas muy largas con superficies grises difusas. Determine la irradiación y la radiosidad para la placa superior. ¿Cuál es el intercambio neto de radiación entre las placas por unidad de área?
′′ = 42.525 W/m2. Solución: G1 = 14.175 W/m2; J 1 = 56.700 W/m2; q12 7.
(13.51 del Incropera) Considere el recinto de tres superficies que se muestra en la figura. La placa inferior ( A1) es un disco negro de 200 mm de diámetro y se le suministra calor a razón de 10.000 W. La placa superior ( A2), un disco coaxial a A1, es una superficie gris difusa con ε 2 = 0,8 y T 2 = 473 K. Los lados, grises y difusos, entre las placas están perfectamente aislados. Suponga que la transferencia de calor por convección es insignificante. Determine la temperatura de operación de la placa inferior, T 1, y la temperatura de los lados aislados, T 3.
53
Transferencia de Calor / Curso 2009-10
Problemas
Solución: T 1 = 1.823 K; T 3 = 1.581 K. 8.
(13.57 del Incropera) El extremo de un tanque propulsor cilíndrico de líquido criogénico en el espacio se va a proteger de la radiación externa solar mediante la colocación de un escudo metálico delgado en la parte frontal del tanque. Suponga que el factor de forma F ts entre el tanque y el escudo es la unidad. Todas las superficies son difusas y grises y los alrededores están a 0 K. Encuentre la temperatura del escudo, T s, y el flujo de calor al extremo del tanque.
′′ = 25,4 W/m 2. Solución: T s = 338 K; q st 9.
(13.67 del Incropera) Un elemento cilíndrico largo de calentamiento de diámetro D1 = 10 mm, temperatura T 1 = 1.500 K y emisividad ε 1 = 1 se usa en un horno. El área de la parte inferior ( A2) es una superficie gris difusa con ε 2 = 0,6 y se mantiene a T 2 = 500 K. Las paredes lateral y superior están construidas con ladrillo refractario aislante que es difuso y gris con ε 3 = 0,9. La longitud del horno en dirección normal a la página es muy larga en comparación con el ancho w y el alto h. Ignorando la convección y tratando las paredes del horno como isotérmicas, determine la potencia por unidad de longitud que se debe proporcionar al elemento calentador para mantener condiciones de estado estacionario. Calcule la temperatura de la pared del horno.
54
Transferencia de Calor / Curso 2009-10
Problemas
Solución: q1′ = 8.542 W/m; T 3 = 733 K. Transferencia de calor multimodal 10.
(13.87 del Incropera) Una placa cuadrada, gris, difusa y opaca de 200 mm de lado y con una emisividad de 0,8 se coloca en la abertura de un horno y se sabe que está a 400 K en cierto instante. La parte inferior del horno, que tiene las mismas dimensiones que la placa, es negra y opera a 1.000 K. Las paredes laterales del horno están bien aisladas. La parte superior de la placa se expone al aire ambiente con un coeficiente de convección de 25 W/m2·K. El aire y los alrededores están a 300 K. Suponga que la convección en la superficie inferior de la placa es insignificante.
a) Evalúe la transferencia neta de calor radiativa para la superficie inferior de la placa. b) Si la placa tiene una masa de 2 kg y una calor específico de 900 J/kg·K, ¿cuál será el cambio en la temperatura de la placa con el tiempo, dT / dt ? Solución: a) q2, rad = 1.153 W (entrante); b) dT / dt = 0,57 K/s. 11.
(13.93 del Incropera) La absortividad espectral de una superficie difusa grande es α λ = 0,9 para λ < 1 µ m y α λ = 0,3 para λ > 1 µ m. La parte inferior de la superficie está bien aislada, mientras que la superior se puede exponer a dos condiciones diferentes.
55