Balances Elementales de Membrana y Valores

en Variables de Procesos Estacionarios de Transferencia

Dr. Víctor Roberto Parra Sánchez

 Análisis Envolvente de Energía (1)
En palabras, el análisis envolvente de energía (láminas o envolturas normales a la
dirección de conducción de calor) para un sistema en estado estable:

Velocidad de Velocidad de Velocidad

entrada de salida de de
energía - energía + producción =0
calorífica calorífica de energía
calorífica

¿De qué manera puede entrar y salir energía calorífica del sistema?
¿De qué manera se puede generar energía calorífica en el sistema?
 Diferencias entre Transporte
Conductivo y Convectivo de Calor

Conducción: Generación de energía

Ley de Fourier calorífica:
• Degradación de energía
dT eléctrica.
qx   k • Disminución de velocidad de
dx
neutrones o fragmentos
Convección: nucleares durante la fisión.
• Degradación de energía
Transporte de energía
mecánica (disipación
calorífica debido al
viscosa).
movimiento global del
• Conversión de energía
fluido (calor sensible)
química a calorífica.
 Solución de Balances de Energía Calorífica
Tras el planteamiento del balance, el espesor del volumen de análisis puede hacerse que tienda
a cero, lo que conduce finalmente a una ecuación diferencial para la distribución de
temperatura.

Al integrar esta ecuación diferencial, aparecen constantes de integración que pueden evaluarse
a través de las condiciones en la frontera o condiciones límite:

a) La temperatura de la superficie conocida; T = T0

b) Densidad de flujo de calor en la superficie determinada, q = q0 (especificar el
gradiente de temperatura).
c) Para interfase sólido-fluido, la densidad de flujo de calor puede relacionarse con la
diferencia de temperaturas entre la superficie del sólido y la del fluido:

q  h T  T f 
d) En interfases sólido-sólido puede ocurrir la continuidad de temperatura y densidad de
flujo de calor.
Conducción del Calor Con un Manantial
Calorífico de Origen Eléctrico
r Características del sistema:
• Alambre de sección circular de radio R y conductividad
z eléctrica ke (ohm-1 cm-1).
• A través del alambre circula una corriente eléctrica I amp/cm2.
• Parte de la energía eléctrica se transforma en energía
calorífica (proceso irreversible), cuya rapidez de producción
por unidad de volumen está dada por:
I2
Se 
ke
(calor que se genera debido a la disipación eléctrica)
• No se tomará en cuenta variaciones en las propiedades del
material debido a aumento de la temperatura (poco
significativo).
• La temperatura de la superficie del alambre se mantiene a
T = T0.
• Diferencias de temperatura en dirección r.
 Balance de Energía

2rL qr r   2 r  r L qr r  r   2rrL Se  0

2rL qr r   2 r  r L qr r  r   2rrL Se

Dividiendo entre -2πΔrL

r  r q r r  r  rq   S r
r r

r
e
 Balance de Energía
Haciendo Δr  0:

lim 
  
 r  r qr r  r  rqr   S r

r
r
e
r 0
 
d rqr 
 Se r
dr
d rqr   S e rdr
2
Se r
rqr   C1
2
 Balance de Energía
despejando qr:
S e r C1
qr  
2 r
¿Cómo se evalúa C1?
en r = 0, qr no puede ser infinito, por lo tanto, C1 = 0, es decir:

r 0 qr  finito C1  0
Se r
qr 
2
 Balance de Energía
Para determinar el perfil de temperaturas en el alambre, se sustituye la ley de
Fourier:
dT S e r
k 
dr 2
Despejando:
Se
dT   rdr
2k
Integrando:
Se 2
T   r  C2 ¿Cómo se evalúa C2?
4k
 Balance de Energía
Se toma la segunda condición en la frontera, que dice que en la superficie
externa del cilindro, la temperatura es constante e igual a T0, es decir:

rR T  T0

Se 2
T0   R  C2
4k
Se 2
C2  T0  R
4k
 Balance de Energía
Sustituyendo en el perfil de temperaturas:

Se 2 Se 2 Se 2
T   r  C2 T   r  T0  R
4k 4k 4k

T  T0 
4k

Se 2 2
R r 
Se R 2
 r 2 2
Se R  r 
2

T  T0  1  2  T  T0  1    
4k  R  4k   R  
 Conducción del Calor Con un Manantial
Calorífico de Origen Eléctrico (2)
Se r
qr 
2
Se R 
2
r 
2

T  T0  1    
4k   R  

Incremento máximo de temperatura (r = 0):

2
Se R
Tmáx  T0 
4k
 Conducción del Calor Con un Manantial
Calorífico de Origen Eléctrico (3)
Incremento medio de temperatura (tarea):

2
T (r )  T rdrd
R

T  T0 
 
0 0
0
 2 R
  rdrd
0 0

2
Se R
T  T0 
8k
 Conducción del Calor Con un Manantial
Calorífico de Origen Eléctrico (4)
Flujo de calor en la superficie (para la longitud L del alambre):

Q r  R  2RL  qr r R

Se R
Q r R  2RL 
2

Q r  R  R LS e
2

¿Qué implica esta última ecuación?

 Comparación entre el Problema de Flujo Viscoso
y el Alambre Caliente

Flujo Viscoso en un Tubo Alambre Calentado

De la primera integración se obtiene la distribución de τrz qr
De la segunda integración se obtiene la distribución de vz T – T0
La condición límite para r = 0 es τrz = finito qr = finito
La condición límite para r = R es vz = 0 T – T0
Propiedad de transporte µ k
Término que corresponde al manantial (P0 – PL)/L Se
Suposiciones µ = constante k, ke = constante
 Problema
Un alambre de cobre tiene 2 mm de radio y 5 m de longitud. ¿Qué caída de
tensión dará lugar a un aumento de temperatura en el eje del alambre de 10 °C,
si la temperatura de la superficie externa del mismo es de 20 °C?
2
Se R
Tmáx  T0 
4k
2
I
Se 
ke
2 2
I R
Tmáx  T0 
4kke
 Problema
E
I  ke
L
 E 2 R 2  ke 
Tmáx  T0   2 
 
 4 L  k 

L k
E  2  T0 Tmáx  T0 
 R  keT0
 Problema
Para el cobre, el número de Lorentz k/keT0 es 2.23 x 10-8 vol2 K-2 . Así la
diferencia de voltaje necesaria para producir un aumento de temperatura de
10 °C es:

 5000mm  8  volt 
E  2  2.23 10   29310K
 2mm   K 


E  5000 1.49 10 4
54.1volt
E  40 volt
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias
Comúnmente se tratan aplicaciones en las cuales se lleva a cabo la transmisión
de calor a través de paredes constituidas por capas de varios materiales:

• Pared compuesta formada por tres

materiales de distinto espesor
x1 – x0 , x2 – x1 y x3 – x2
• Distintas conductividades caloríficas
k01, k12 y k23.
• En x = x0 la sustancia 01 está en
contacto con un fluido a la
temperatura Ta y la sustancia 23 está
en contacto con un fluido a la
temperatura Tb.
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (2)

• La transmisión del calor en las superficies externas (superficies) está dada por la
ley de enfriamiento de Newton, siendo h0 y h3 los coeficientes de transmisión de
calor.
• Aplicando un balance de energía para la lámina del material 01 con volumen
WHΔx, se obtiene:
q 01
x WH  q 01
x WH  0
x x  x

• Invirtiendo el orden de la entrada y salida y dividiendo entre WHΔx además de

hacer el espesor de la película lo más pequeño posible se tiene:
 q 01  q 01

lim  x 
x x
x  x
0
x 0 x 
 
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (3)

Por ser la definición de la derivada, se tiene:

dqx01
0
dx
Por lo tanto:

qx01  C1

Por la condición límite:

x  x0 qx01  q0
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (4)

En el estado estable, la densidad de flujo de calor debe ser la misma en las tres
regiones:

q  q  q  q0
01
x
12
x
23
x

Se sabe también por ley de Fourier que:

01
dT
q  k
01
x
01

dx
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (6)

Asimismo, por analogía:

01
dT
k 01
 q0
dx
12
dT
k 12
 q0
dx
23
dT
 k 23  q0
dx
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (7)

Integrando con los límites, x0, T0 y para x1, T1:

q0
dT 01
  01 dx
k
T0 q0 x0
   01  dx
01
dT
T1 k x1

T0  T1   q0
x0  x1 
k 01
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (8)

De la misma forma, para T12 y para T23:

T1  T2   q0
x1  x2 
k 12

T2  T3   q0
x2  x3 
k 23
Sobre la superficie, se tiene una transferencia de calor entre el sólido y un fluido que
puede ser explicada a través de la ley de enfriamiento de Newton:

q  h Ts  T f 
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (9)
Para las dos superficies del problema:
q0
Ta  T0 
h0
q0
T3  Tb 
h3
Sumando las ecuaciones del perfil de temperaturas 01, 12, 23 y las de los extremos:

 1 x1  x0 x2  x1 x3  x2 1 
Ta  Tb  q0   01
 12  23
 
 h0 k k k h3 
Conducción del Calor a través de Paredes Compuestas:
Suma de Resistencias (10)
Así:
Ta  Tb
q0 
1 3
xi  xi 1 1 
   i 1,i  
 h0 i 1 k h3 
Comúnmente se escribe la ecuación anterior de la forma de la ley de enfriamiento de
Newton:
q0  U Ta  Tb  Q0  U WH Ta  Tb 
Donde U se le conoce como coeficiente global de transferencia de calor:
1
1 3
xi  xi 1 1 
U     i 1,i  
 h0 i 1 k h3 
Conducción de Calor en una Aleta de Enfriamiento
Una aplicación sencilla, pero práctica, para la conducción de calor es el cálculo de la
eficacia de una aleta de enfriamiento (utilizadas para aumentar la superficie disponible
para la transmisión de calor en paredes metálicas y fluidos que son muy malos
conductores como los gases).
Conducción de Calor en una Aleta de Enfriamiento (2)
Características de los casos reales y aproximados por el modelo:

Caso Real Modelo

1.- T es función de z y x pero la 1.- T es sólo función de z.
influencia de z es más importante
2.- Por el extremo de la aleta (área 2.- No hay pérdida de calor por el
2BW) y los bordes (área 2BL + 2BL) se extremo y los bordes.
pierde una pequeña cantidad de calor.
3.- El coeficiente de transmisión de 3.- La densidad de flujo de calor en la
calor es función de la posición. superficie viene dada por q = h(T – Ta)
en la que h es una constante y T = T(z).
Conducción de Calor en una Aleta de Enfriamiento (3)
Aplicando un balance de energía a un segmento Δz de la aleta:

q z z  2 BW  q z z  z  2 BW  h 2Wz T  Ta   0

Dividiendo entre 2BWΔz, reordenando y tomando para Δz un valor muy pequeño:

 T  Ta 
dqz h

dz B
Introduciendo la ley de Fourier (qz = -kdT/dz), siendo k la conductividad calorífica del
metal que se considera como constante, se obtiene:
2
d T
 2 
h
T  Ta 
dz kB
Conducción de Calor en una Aleta de Enfriamiento (4)
Dado que es una ecuación diferencial de segundo orden se requieren las siguientes
condiciones límite para su resolución:

C.L.1 z0 T  Tw

dT
C.L.2 zL 0
dz
Se introducen las siguientes variables adimensionales:

T  Ta z N
hL2
 
Tw  Ta L kB
Coeficiente adimensional
Temperatura adimensional Distancia adimensional de transmisión de calor
Conducción de Calor en una Aleta de Enfriamiento (5)
De acuerdo a lo anterior, el problema puede replantearse de la siguiente manera:

d 
2
 N 2

d 2

Con:

C.L.1   0  1

d
C.L.2 0
d  1
Conducción de Calor en una Aleta de Enfriamiento (6)
Al integrar la ecuación y determinando las constantes de integración a través de las
condiciones límite se tiene la siguiente ecuación hiperbólica:

  cosh N  tanh N senhN

O de otra forma:

cosh N 1   

cosh N