1

UNIDAD I
MODOS BASICOS DE TRANSFERECNIA DE CALOR Y CONDCCION UNIDIMENSIONAL Y BI-
DIMENSIONAL

Modos básicos de transferencia de calor conducción, convección y radiación.

En Termodinámica se aprendió que la energía se puede transferir por medio de las

interacciones de un sistema con su alrededor, las mismas que se denominan trabajo y calor.

Pero la termodinámica no proporciona información con respecto a la naturaleza de esta

interacción o la rapidez con la que se produce.

¿Qué es la transferencia de calor?

Es la energía en tránsito debido a una diferencia de temperaturas.

 Siempre que exista una diferencia de temperaturas en un cuerpo o entre cuerpos, debe
ocurrir una transferencia de calor.

Procesos o Modos de transferencia de calor

Conducción:

Cuando existe una diferencia de temperaturas en un medio estacionario (puede ser sólido o un
fluido) se utilizara el termino conducción para referirnos a la transferencia de calor que se puede
producir a través del medio.

Convección:

Se refiere a la transferencia de calor que ocurrirá entre una superficie y un fluido en movimiento
cuando están a diferentes temperaturas.
2

Radiación Térmica:

Todas las superficies con temperatura finita emiten energías en forma de ondas
electromagnéticas.

Por tanto en ausencia de un medio, existe una transferencia neta de calor por radiación entre
dos superficies a diferentes temperaturas.

CONDUCCIÓN

Debemos inmediatamente referirnos a conceptos de actividad atómica y molecular, pues hay

procesos en estos niveles que sustentan este modo de transferencia de calor.

La conducción se considera como la transferencia de energía de las partículas más energéticas

a las menos energéticas de una sustancia debido a las interacciones entre las mismas.
3

Se habla de transferencia neta de energía debido al movimiento molecular aleatorio como una
difusión de energía.

La situación es muy similar en los líquidos aunque las moléculas están menos esparcidas y las
interacciones moleculares son más fuertes y frecuentes. De igual manera en un sólido la
conducción se atribuye a la actividad atómica en forma de vibraciones reticulares,
modernamente se atribuye la transferencia de energía a ondas reticulares inducidas por el
movimiento atómico (en un no conductor se da esto). En un conductor además, tenemos el
movimiento de traslación de electrones libres.

Es posible cuantificar los procesos de transferencia de calor en términos de ecuaciones o

modelos apropiados, estos sirven para calcular la cantidad de energía que se transfiere por
unidad de tiempo.

Ley de Fourier

Para pared plana unidimensional:

El signo – es consecuencia de que el calor que se transfiere en dirección

de la temperatura decreciente.

K: conductividad térmica (W/m*K); característica del material

𝒒′′ 𝒙 (𝑾/𝒎𝟐 ): Velocidad con que se transfiere el calor por unidad de área perpendicular a la
dirección de transferencia y es proporcional al gradiente de temperatura.
4

Para condiciones de estado estable el gradiente se transforma:

Esta ecuación proporciona un flujo de calor, es decir velocidad del calor transferido por unidad
de área.

El calor transferido por unidad de tiempo 𝒒𝒙 (𝑾)a través de una pared plana de área A es:

𝒒𝒙 = 𝒒′′ 𝒙 ∗ 𝑨

Ejemplo:

La pared de un horno industrial se construye con ladrillo de arcilla refractaria de 0.15 m d

espesos, que tiene una conductividad térmica de 1.7 (W/m*K). Mediciones realizadas durante
la operación en estado estable revelan temperaturas de 1400 ºK y 1150 ºK en la superficie
interna y externa.

¿Cuál es la velocidad de perdida de calor a través de una pared que tiene 0.5 m por 3 m de lado?
5

Problemas de conducción:

1. Si un flujo de calor de 3 KW se conducen a través de una sección de un material aislante

de área de sección transversal 10𝑚2 y espesor 25 cm. Si la temperatura de la superficie
interna (caliente) es de 415 ºC y la conductividad térmica del material es 0.2W/m*K
¿Cuál es la temperatura de la superficie externa?

Datos:

q= 3Kw

k=0.2 W/m*K

A=10𝑚2

e=0.25m

T=415 ºC=688 ºK
6

∆𝑇
𝑞′′𝑥 = −𝑘
𝐿
𝑞 = 𝑞′′𝑥 ∗ 𝐴
3000(𝑊)
𝑞′′𝑥 =
10𝑚2
𝑊
𝑞 ′′ 𝑥 = 300
𝑚2
𝑇2 − 𝑇1
𝑞′′𝑥 = −𝑘
𝐿
𝑞′′𝑥 𝐿 = −0.2𝑇2 + 0.2(688)

300 ∗ 0.25 = −0.2𝑇2 + 137.6

𝑇2 = 313ºK = 40ºC
2. Se determina que el flujo de calor a través de una tabla de madera de 50 mm de espesor,
cuyas temperaturas sobre las superficies interna y externa son 40 ºC y
20 ºC respectivamente es de 40 W/𝑚2
¿Cuál es la conductividad térmica de la madera?

Datos:

𝑞′′𝑥 = 40 W/𝑚2

e=50mm

Ti=40 ºC

Te=20 ºC

∆𝑇
𝑞′′𝑥 = −𝑘
𝐿
𝑞′′𝑥 ∗ 𝐿
− = −𝑘
∆𝑇
0.05𝑚
𝑘 = −40 ∗
(20 − 40)𝐶
𝑊
𝑘 = 0.1
𝑚∗𝐾
7

3. La temperatura de las superficies interna y externa de una ventana de vidrio de 5mm de

espesor son 15 y 5 ºC ¿Cuál es la perdida de calor a través de una ventana que mide
𝑊
1*3m de lado? La k del vidrio es 1.4𝑚∗𝐾

Datos:

e=5mm

Ti=15 ºC=288 ºK

Te=5 ºC=278 ºK
𝑊
k=1.4
𝑚∗𝐾

h=3m

a=1m

𝑊 (278 − 288)ºK
𝑞′′𝑥 = −1.4 ∗
𝑚𝐾 0.005𝑚
𝑊
𝑞′′𝑥 = 2800
𝑚2
𝑊
𝑞 = 2800 ∗ 1 ∗ 3𝑚2
𝑚2
𝑞 =8.4 KW

Ejercicios sugeridos: 1.2-1.4-1.5-1.6-1.7

8

CONVECCIÓN

Se compone de dos mecanismos: movimiento molecular aleatorio (difusión) y movimiento

global o macroscópico del fluido, este se da en presencia de un gradiente de temperatura que
contribuye a la transferencia de calor.

La transferencia se da por una superposición de transporte de energía por el movimiento

aleatorio de las moléculas y por el movimiento global del fluido.

Convección transporte acumulado del fluido.

La convección ocurre entre un fluido en movimiento y una superficie limitante cuando estos
tienen diferentes temperaturas.

La velocidad varía desde cero sobre la superficie a un valor finito µ∞ asociado con el flujo.

Esta región se la denomina capa límite hidrodinámica o de velocidad.

Además si las temperaturas de la superficie y del fluido difieren, habrá una región del fluido a
través del cual la temperatura varía de 𝑇𝑠 en y=0 a 𝑇∞ en el flujo exterior. Esta región
denominada capa limite térmica puede ser más pequeña, más grande o del mismo tamaño que
aquella en la que varía la velocidad, de todos modos si 𝑇𝑠 > 𝑇∞ ocurrirá la transferencia de calor
por convección entre la superficie y el flujo exterior.

Flujo de aire 𝑻∞ 𝒒′′

𝑇𝑠

Caloportador L

𝑇𝑠 > 𝑇∞
9

Se analiza en relación al fluido caloportador mas no al sólido, existen dos tipos de convección.

Convección Forzada: cuando el caloportador es obligado a pasar con una velocidad o caudal.

Convección Natural: cuando el caloportador está en reposo

La convección se da por el movimiento del fluido inducido por las burbujas de vapor generadas
en el fondo de una cacerola en la que está hirviendo agua.
10

Ley de Enfriamiento de Newton

𝑞′′𝑐 = 𝑞 ′′ = ℎ(𝑇𝑠 − 𝑇∞ )

𝑞′′𝑐 = flujo de calor por convección (𝑊/𝑚2 )

𝑇𝑠 =temperatura de la superficie

𝑇∞ =temperatura del fluido

h= coeficiente de transferencia por convección, es relativo al caloportador.

𝑞𝑐 = 𝑞′′𝑐 ∗ 𝐴

A= área de convección ( 𝑚2 )

RADIACIÓN

La radiación térmica es la energía emitida por la materia que se encuentra a una temperatura
finita. Aunque la atención será la radiación de superficies solidas, también existe radiación
térmica de líquidos y gases. La energía emitida entre superficies a diferente temperatura se da
por ondas electromagnéticas.

Si la radiación se produce en vacío la

Fluido transferencia es más efectiva

A B

𝑇𝐴 = 𝑇𝐵

La radiación que la superficie emite se origina a partir de la energía

térmica de la materia limitada por la superficie y la velocidad a la
que libera energía por unidad de área (𝑊/𝑚2 ).

Se denomina Potencia Emisiva Superficial

Ley de BOLTZMANN.

𝑬𝒃 = 𝝈 𝑻𝒔 𝟒 : La superficie se denomina cuerpo negro o radiador

ideal de la superficie

𝝈: Constante de BOLTZMANN
𝑾
𝝈 = 𝟓. 𝟔𝟕 ∗ 𝟏𝟎−𝟖
𝒎𝟐 𝑲𝟒
𝑻𝒔 : Temperatura absoluta (°K)

Si el flujo de calor emitido por una superficie real es menor que la de un cuerpo negro:
11

𝑬𝒃 = 𝜺 ∗ 𝝈 ∗ 𝑻𝒔 𝟒
𝜺: emisividad 0≤ 𝜺 ≤1

Es una propiedad que proporciona la eficiencia con la que una superficie emite energía en
relación con un cuerpo negro.

La emisividad depende del material, de la superficie y del acabado.

𝜺=α, área A y temperatura 𝑻𝒔

La radiación puede incidir sobre una superficie desde sus alrededores. La radiación se origina
desde una fuente especial, como el sol, o de otras superficies a la que se expone la superficie de
interés.

Sin tomar en cuenta la fuente designamos la velocidad a la que tarda esa radiación incide sobre
un área unitaria de la superficie como la irradiación G.

La velocidad a la que la energía radiante es absorbida por área superficial unitaria se la evalúa a
partir de una propiedad radiativa de la superficie denominada absortividad.

donde:

α: absortividad

G: irradiación G= 𝝈 ∗ 𝑻𝒔 𝟒

0≤α≤1

Si α<1 y la superficie es opaca, partes de la irradiación se reflejan.

Si la superficie es semitransparente partes de la irradiación se transmiten.

Si suponemos que la superficie es tal que α=ε(superficie gris), la velocidad neta de transferencia
de calor por radiación desde la superficie expresada por unidad de área de la superficie es:
12

Diferencias Energía térmica liberada debido a la emisión por radiación y la que se gana debido a
la absorción por radiación.

Existen muchas aplicaciones en las que conviene expresar el intercambio neto de calor por
radiación en la forma.

Las superficies pueden transferir simultáneamente calor por convección a un gas contiguo;
donde la velocidad total de transferencia de calor desde la superficie es:

Ejemplo 1.2

Una tubería de vapor sin aislamiento pasa a través de un cuarto en el que el aire y las paredes
están a 25°C. El diámetro exterior de la tubería es 70 mm y la temperatura superficial y
emisividad son 200°C y 0.8 respectivamente.

¿Cuánto vale la potencia emisiva de la superficie y la irradiación? Si el coeficiente asociado con

la transferencia de calor por convección libre de la superficie al aire es 15 𝑊/𝑚2 °𝐾

¿Cuál es la velocidad de perdida de calor de la superficie por unidad de longitud de la tubería?

Datos:

Aire

𝑇∞ = 25°𝐶

h=15 𝑊/𝑚2 °𝐾

𝑇𝑠 = 200°𝐶
ε=0.8 𝑇∞
=
∅=70mm

CEE: perdida de calor por convección

𝐺 = 𝜎 ∗ 𝑇𝑎𝑙𝑟 4 irradiación

Potencia emisiva
𝑊
𝐸 = 𝜀 ∗ 𝜎 ∗ 𝑇𝑠 4 = 0.8 (5.67 ∗ 10−8 ) (200 + 273)4 °𝐾 4
𝑚2 𝐾 4
𝑊
𝐸 = 22700
𝑚2
13

𝑊
𝐺 = 5.67 ∗ 10−8 ∗ (25 + 273)4 °𝐾 4
𝑚2 𝐾 4
𝑊
𝐺 = 447
𝑚2
Perdida de calor en la tubería

Por convección con el aire del cuarto y por intercambio de radiación con las paredes

𝐴 = 𝜋𝐷𝐿

𝑞 = ℎ𝐴(𝑇𝑠 − 𝑇∞ ) + 𝜀𝐴𝜎(𝑇𝑠 4 − 𝑇𝑎𝑙𝑟 4 )

𝑊 𝑊
𝑞´ = 15 2
(𝜋 ∗ 0.01𝑚)(473 − 298) + 0.8(𝜋 ∗ 0.01𝑚) ∗ 5.67 ∗ 10−8 2 4 (4734
𝑚 °𝐾 𝑚 𝐾
− 2984 )𝐾 4
𝑾
𝒒´ = 𝟗𝟗𝟖
𝒎
hr = εσ(𝑇𝑠 + 𝑇𝑎𝑙𝑟 )(𝑇𝑠 2 + 𝑇𝑎𝑙𝑟 2 )
𝑊
hr = 0.8 ∗ 5.67 ∗ 10−8 (473 + 298)°𝐾(4732 + 2982 )
𝑚2 𝐾 4
𝑾
𝐡𝐫 = 𝟏𝟏
𝒎𝟐 °𝑲
𝑞𝑟 = ℎ𝑟(𝑇𝑠 − 𝑇𝑎𝑙𝑟 )
𝑊
𝑞𝑟 = 11 (200 − 25)°𝐾
𝑚2 °𝐾
𝑊
𝑞𝑟 = 1925
𝑚2
EJERCICIOS: 1.1/1.4/1.5/1.6/1.7/1.9/1.11/1.17

Balance de Energía

Para realizar el análisis de balance de energía debemos establecer al igual que en termodinámica
el volumen de control.

La primera ley se aplica una vez determinado e identificado el volumen de control, una región
del espacio limitado por una superficie de control a través de la cual pueden pasar la energía y
la materia.
14

SISTEMA

COND

CONV

RAD
T
COND

CONV

RAD

Aplicamos conservación de Energía al Volumen de control

𝑑𝐸𝑎𝑙𝑚
𝐸𝑒𝑛𝑡 + 𝐸𝑔 − 𝐸𝑠𝑎𝑙𝑒 = = 𝐸𝑎𝑙𝑚
𝑑𝑡

En un instante estos términos incluyen la velocidad a la que la energía térmica y mecánica entra
y sale de la superficie de control.

Ealm: Este tipo de energía se almacena dentro del sistema cuando la temperatura en su interior
se incrementa o decrementa con el tiempo.

𝑑𝑇
𝐸𝑎𝑙𝑚 = 𝜌 ∗ 𝑉 ∗ 𝐶𝑝 ∗
𝑑𝜃

𝑬𝒂𝒍𝒎: Energía Almacenada [BTU/h]

𝝆: Densidad del material [lb/p3]

𝑽: Volumen p3

𝑪𝒑: Capacidad calorífica [BTU/lb°K]

𝒅𝑻
𝒅𝜽
: Diferencial de temperatura en función del tiempo

𝜽: Tiempo

Eg: Este tipo de energía se la obtiene cuando otras formas de energía en el interior del sistema
se convierten en energía térmica.
15

𝐸𝑔 = 𝑔 ∗ 𝑉

𝒈: Tasa de generación por unidad de volumen [w/m3]

𝑽: Volumen [m3]

Ejemplo: fricción, electricidad (resistencia) Reacción combustible.

En el balance solo hay energías térmicas, no hay energías mecánicas ni trabajo.

Balance de energía en una superficie

Con frecuencia aplicaremos el requerimiento de conservación de energía a la superficie de un

medio. En este caso especial la superficie de control no incluye masa o volumen.

La Eg y Ealm no son relevantes y es necesario

tratar el fenómeno superficial

𝐸𝑒𝑛𝑡 − 𝐸𝑠𝑎𝑙𝑒 = 0

Tenemos 3 formas de transferencia de calor

para la superficie de control, conducción desde
el medio hacia la superficie de control q”cond;
convección desde la superficie al fluido e
intercambio de radiación desde la superficie
hacia los alrededores q”rad.

 El balance toma la forma

q”cond – q”conv – q”rad = 0

16

Ejercicio 1.1 , Ejemplo 1.5

Los gases calientes de combustión de un horno se separan del aire ambiental y sus alrededores,
que están a 25°C mediante una pared de ladrillos de 0,15m de espesor. El ladrillo tiene una
conductividad térmica de 1,2 W/(m°K) y una emisividad superficial de 0,8. Se mide una
temperatura de la superficie externa de 100°C en CEE. La transferencia de calor por convección
libre del aire contiguo a la superficie se caracteriza por un coeficiente de convección de h=20
w/m2°K. ¿Cuál es la temperatura de la superficie interior del ladrillo?

1. CEE
Solución
𝐸𝑛𝑡 − 𝐸𝑠𝑎𝑙𝑒 = 0

𝑞"𝑐𝑜𝑛𝑑 – 𝑞"𝑐𝑜𝑛𝑣 − 𝑞"𝑟𝑎𝑑 = 0

𝑇1 −𝑇2 4
𝐾 𝐿
= ℎ( 𝑇2 − 𝑇∞ ) + 𝜀𝜎(𝑇24 − 𝑇𝑎𝑙𝑟 )

𝑊 (𝑇1 − 373)
1,2
𝑚𝐾 0,15
𝑊 𝑊
= 20 2
(373 − 298) + 0,8 (5,67 ∗ 10−8 2 ) (3734 − 2984 )
𝑚 𝐾 𝑚 𝐾

T1 = ? 𝑊 𝑊
= 1500 𝑚2 + 520 𝑚2 = 2020

0,15𝑚 𝑊
𝑇1 = 373𝐾 + (2020 2 ) = 625𝐾
𝑊 𝑚
1,2 𝑚𝐾

= 352 ℃

Ejercicios 1.7/1.9/1.23/2.22/2.25

Ejemplos 1.6/1.7
17

Ejemplo 1.6

El recubrimiento sobre una placa se cura exponiendo ésta a la acción de una lámpara infrarroja
que proporciona una irradiación de 2000 W/m2. El recubrimiento absorbe 80% de la irradiación
y tiene una emisividad de 0,50 también expuesta a un flujo de aire y a amplios alrededores para
los cuales las temperaturas son 20°C y 30°C respectivamente.

1. Si el coeficiente de convección entre la placa y el aire ambiente es 15 W/m2K. ¿Cuál es

la temperatura de curación de la placa?
2. Las características finales del recubrimiento, incluidos uso y durabilidad se sabe que
dependen de la temperatura a la que ocurre la curación. Un sistema de flujo de aire es
capaz de controlar la velocidad del aire sobre la superficie cerrada, pero el ingeniero de
procesos necesita saber en que forma depende la temperatura del coeficiente de
convección.
Proporcione la información deseada con el cálculo y aplicación de la temperatura de la
𝑊
superficie como función de h 2 ≤ ℎ ≤ 200 𝑚2 𝐾
. Que valor de h proporcionara una
temperatura de 50°C.

𝐸𝑒𝑛 = 𝐸𝑠𝑎𝑙𝑒

(∝ 𝐺)𝑙𝑎𝑚𝑝 − 𝑞𝑐𝑜𝑛𝑣 − 𝑞𝑟𝑎𝑑 = 0

4
∝ 𝐺𝑙𝑎𝑚𝑝 − ℎ(𝑇 − 𝑇∞ ) − 𝜀𝐺(𝑇 4 − 𝑇𝑎𝑙𝑟 )=0

𝑊 𝑊 𝑊
0,8 ∗ 2000 2
− 15 2 (𝑇 − 293)𝐾 − 0,5 ∗ 5,67 ∗ 10−8 2 (𝑇4 − 3044 )𝐾 4
𝑚 𝑚 𝐾 𝑚 𝐾

𝑇 = 377°K → 104℃
18

Ejemplo 1.23

Se conecta un resistor eléctrico a una batería como se muestra en el esquema. Después de una
breve fluctuación transitoria la resistencia toma una temperatura de EE casi uniforme de 95℃,
mientras que la batería y los alambres de conexión permanecen a la temperatura ambiente de
25℃.

No tome en cuenta la resistencia eléctrica de los alambre de conexión.

a) Considere el resistor como un sistema alrededor del cual se coloca una superficie de
control y se aplica 𝐸𝑒𝑛𝑡 + 𝐸𝑔 − 𝐸𝑠𝑎𝑙𝑒 = 𝐸𝑎𝑙𝑚. Determine los valores
correspondientes de estas energías.
Si se coloca una superficie de control alrededor del sistema entero. ¿Cuáles son los
valores de Eent, Eg, Esale, Ealm.
b) Si se disipa energía eléctrica de manera uniforme dentro del resistor que es un cilindro
de ∅ = 60 𝑚𝑚 y longitud L = 25 cm. ¿Cuál es la velocidad de generación de calor
volumétrica 𝑞̇ (W/m3) (tasa de generación).
c) Sin tener en cuenta la radiación del resistor. ¿Cuál es el coeficiente de convección?

a.

𝐸𝑒𝑛𝑡 + 𝐸𝑔 − 𝐸𝑠𝑎𝑙𝑒 = 𝐸𝑎𝑙𝑚

𝑃 = 𝑉𝐼 = 24𝑉 ∗ 6𝐴 = 144𝑊

Eg, Ealm

𝐸𝑖𝑛 = 0 Ein

V.C Esale

𝐸𝑜𝑢𝑡 = 144𝑊
19

𝐸𝑔𝑒𝑛 = 144𝑊

𝐸𝑎𝑙𝑚 = 0

Eg, Ealm

Ein

Esale

𝐸𝑎𝑙𝑚 = 0
𝐸𝑔 = 0

Eentra = 144W
Esale = 144W

b.

𝐸𝑔 = 𝑔 ∀

𝜋𝐷 2
𝑉= ∗𝐿
4

𝜋(0,06)2
𝐸𝑔 = 144 = 𝑔 ( ) ∗ 0,25
4

𝑊
𝑔 = 2,04 ∗ 105
𝑚3

c.
20

AT 
𝜋𝐷 2
= 𝐴1
4

𝐴 3 = 𝜋𝐷𝐿

𝜋𝐷 2
= 𝐴2
4

𝜋𝐷 2
𝐴 𝑇 = 𝜋𝐷𝐿 + 2
4

𝐸𝑜𝑢𝑡 = ℎ 𝐴𝑠 (𝑇𝑠 − 𝑇∞)

0,062 2
144 = ℎ [𝜋 ∗ 0,06𝑚 ∗ 0,25𝑚 + 2 (𝜋 𝑚 )] (95 − 25)
4

144 = ℎ [0,0471 + 0,0057]𝑚2 (95 − 25)

𝑊
ℎ = 39
𝑚2 °K
21

Conducción en Estado Estable

∆𝑇
𝑞 = 𝐾𝐴
∆𝑋

El signo (-) es necesario porque el calor se

transfiere en la dirección que la
temperatura decrece.
𝑞𝑥 𝑑𝑇
𝑞" = = −𝑘
𝐴 𝑑𝑋

El flujo de calor siempre será normal hacia una superficie de temperatura constante
denominado isotérmica.

Las superficies isotérmicas son planos normales a la dirección de X.

El flujo de calor es una cantidad vectorial.

𝜕𝑇 𝜕𝑇 𝜕𝑇
𝑞" = −𝑘∇𝑇 = −𝑘 ( 𝑖+ 𝑗+ 𝑘)
𝜕𝑥 𝜕𝑦 𝜕𝑧

∇: 𝑛𝑎𝑏𝑙𝑎

donde;

nabla (∇): es el operador tridimensional

T(x, y, z): es el campo escalar de temperatura

En coordenadas cartesianas
22

𝑞" = 𝑞"𝑥 𝑖 + 𝑞"𝑦 𝑗 + 𝑞"𝑧 𝑘

𝜕𝑇
𝑞"𝑥 = −𝑘
𝜕𝑥

𝜕𝑇
𝑞"𝑦 = −𝑘
𝜕𝑦

𝜕𝑇
𝑞"𝑧 = −𝑘
𝜕𝑧

Propiedades Térmicas de la materia

𝑞"𝑥 Propiedad
𝑘≡− física
𝜕𝑇
𝜕𝑥

El flujo de calor por conducción aumenta con el incremento de conductividad térmica.

La 𝑘𝑆 > 𝑘𝐿 > 𝑘𝐺

PROPIEDAD TERMODINAMICA

Se relacionan con el estado de equilibrio.

La densidad (𝜌) y el calor específico (Cp)

𝝆 ∗ 𝑪𝒑[𝑱/𝒎𝟑 °𝑲] Capacidad térmica volumétrica

Mide la capacidad de un material para almacenar energía térmica

23

la razón de la conductividad térmica a la capacidad térmica es una importante propiedad

denominada difusividad térmica.
𝒌
𝜶=
𝝆 ∗ 𝑪𝒑
Mide la capacidad de un material conducir energía térmica en relación con la capacidad de
almacenar.

CONDUCCION UNIDIMENSIONAL EN ESTADO ESTABLE.

Pared Plano

𝑻𝟏 > 𝑻𝟐
Balance de Energía

𝑬𝒆 + 𝑬𝒈 − 𝑬𝒔 = 𝑬𝒂𝒍𝒎

𝑬𝒆 = 𝑬𝒔
𝐴𝑑𝑇 𝐴𝑑𝑇
−𝑘 = −𝑘
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=𝐿
𝐴𝑑𝑇 𝐴𝑑𝑇
−𝑘 +𝑘 =0
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=𝐿
Si existe energía generada en el interior de la pared

𝑄𝑥 = 0
𝑄𝑥 = 𝐿
𝑬𝒆
24

𝐴𝑑𝑇 𝐴𝑑𝑇
−𝑘 + 𝐸𝑔 + 𝑘 =0
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=𝐿
𝐸𝑔 = 𝜌 ∗ 𝑉

𝜌: densidad

V: volumen

𝑉 = 𝐴 ∗ ∆𝑥
𝐸𝑔 = 𝜌 ∗ 𝐴 ∗ ∆𝑥
𝐴𝑑𝑇 𝐴𝑑𝑇
−𝑘 + 𝜌 ∗ 𝐴 ∗ ∆𝑥 + 𝑘 =0
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=𝐿
𝑑𝑇 𝑑𝑇
−𝑘 + 𝜌 ∗ ∆𝑥 + 𝑘 =0
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=𝐿
𝑑𝑇 𝑑𝑇
𝑘 −𝑘 + 𝜌 ∗ ∆𝑥 = 0
𝑑𝑥 𝑥=𝐿 𝑑𝑥 𝑥=0
𝑑𝑇 𝑑𝑇
𝑘 −𝑘
𝑑𝑥 𝑥=𝐿 𝑑𝑥 𝑥=0
+𝜌=0
∆𝑥
Si k=cte.

𝑑2 𝑇 𝜌
+ =0
𝑑𝑥 2 𝑘
Modelo matemático de la temperatura con generación

𝑑2 𝑇
=0
𝑑𝑥 2
Modelo matemático de la temperatura sin generación

Clasificación de las condiciones de Borde

Primera clase

𝑇1 𝑇2

𝑥=0 𝑥=𝐿
Se la obtiene cuando se define un valor numérico de temperatura o una función de temperatura
en los límites del sistema.
25

𝑇1 = 𝑓1 (𝑥)
𝑇(0,𝑡) = 𝑇𝑠

𝑇2 = 𝑓2 (𝑥) Temperatura superficial constante

𝑥=0 𝑥=𝐿

Segunda clase

Se la obtiene cuando se define un valor numérico del calor o función de calor en los límites del
sistema.

𝑑𝑇 𝑑𝑇
−𝑘 = 𝑞0 −𝑘 = 𝑞𝐿
𝑑𝑥 𝑑𝑥 𝑥=𝐿

𝑥=0 𝑥=𝐿
𝑞0 = 𝑞1 (𝑥)
𝑞𝐿 = 𝑞2 (𝑥)

Tercera clase

Se la obtiene cuando los límites del sistema están sometidos a transferencia de calor por
convección

Fluido 1
𝑑𝑇
ℎ𝐴(𝑇∞1 − 𝑇𝑥=0 ) = −𝑘 𝑑𝑥
𝑥=0

𝑇2
26

𝑑𝑇
−𝑘 𝑑𝑥 = 𝑥ℎ2 𝐴(𝑇𝑥=𝐿 − 𝑇∞2 )
𝑥=0

De:

𝑑2 𝑇
=0
𝑑𝑥 2
𝑑2 𝑇
∫ =0
𝑑𝑥 2

∫ 𝑑2 𝑇 = − ∫ 𝑑𝑥 2

𝑑𝑇
∫ = 𝐶1 𝑝𝑎𝑟𝑎 𝑐𝑜𝑛𝑑𝑢𝑐𝑐𝑖𝑜𝑛 𝑢𝑛𝑖𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙 𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟 𝑒𝑠 𝑢𝑛𝑎 𝑐𝑡𝑒.
𝑑𝑥
𝑇 = 𝑥𝐶1 + 𝐶2

Si x=0 T1=C2

Si x=L T2=LC1+T1
𝑇2 − 𝑇1
𝐶1 =
𝐿
𝑥
𝑇= (𝑇2 − 𝑇1) + 𝑇1
𝐿
𝑻𝟐 − 𝑻𝟏
𝑻= 𝒙 + 𝑻𝟏
𝑳
Ecuación de distribución de temperatura

Tarea: demostrar la ecuación de la difusión de calor a partir de esta la ecuación de distribución

de temperatura.
𝛿𝑞𝑥
𝑞𝑥 + 𝑑𝑥 = 𝑞𝑥 + 𝑑𝑥
𝛿𝑥
𝛿𝑞𝑦
𝑞𝑦 + 𝑑𝑦 = 𝑞𝑦 + 𝑑𝑦
𝛿𝑦
𝛿𝑞𝑧
𝑞𝑧 + 𝑑𝑧 = 𝑞𝑧 + 𝑑𝑧
𝛿𝑧
𝐸̇𝑒 + 𝐸̇𝑔 − 𝐸̇𝑠 = 𝐸̇𝑎𝑙𝑚

𝛿𝑇
𝐸̇𝑔 = 𝑞̇ 𝑑𝑥 𝑑𝑦 𝑑𝑧 𝐸̇𝑎𝑙𝑚 = 𝜌 𝐶𝑝 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝛿𝑡

𝛿𝑇
𝑞𝑥 + 𝑞𝑦 + 𝑞𝑧 + 𝑞̇ 𝑑𝑥 𝑑𝑦 𝑑𝑧 − 𝑞𝑥 + 𝑑𝑥 − 𝑞𝑦 + 𝑑𝑦 − 𝑞𝑧 + 𝑑𝑧 = 𝜌𝐶𝑝 𝑑𝑥 𝑑𝑦 𝑑𝑧
𝛿𝑥
𝛿𝑇 𝛿𝑇 𝛿𝑇
𝑞𝑥 = −𝑘𝑑𝑦 𝑑𝑧 𝑞𝑦 = −𝑘𝑑𝑥 𝑑𝑧 𝑞𝑧 = −𝑘𝑑𝑥 𝑑𝑦
𝛿𝑥 𝛿𝑦 𝛿𝑧
27

Ecuación de difusión de calor

𝜹 𝜹𝑻 𝜹 𝜹𝑻 𝜹 𝜹𝑻 𝜹𝑻
(𝒌 ) + (𝒌 ) + (𝒌 ) + 𝒒̇ = 𝝆 𝑪𝒑
𝜹𝒙 𝜹𝒙 𝜹𝒚 𝜹𝒚 𝜹𝒛 𝜹𝒛 𝜹𝒕

𝜹𝟐 𝑻 𝜹𝟐 𝑻 𝜹𝟐 𝑻 𝒒̇ 𝟏 𝜹𝑻
+ + + =
𝜹𝒙𝟐 𝜹𝒚𝟐 𝜹𝒛𝟐 𝒌 𝜶 𝜹𝒕
Si en condición de estado estable no existe cambio de energía almacenada
𝛿 𝛿𝑇 𝛿 𝛿𝑇 𝛿 𝛿𝑇
(𝑘 ) + (𝑘 ) + (𝑘 ) + 𝑞̇ = 0
𝛿𝑥 𝛿𝑥 𝛿𝑦 𝛿𝑦 𝛿𝑧 𝛿𝑧
Si la transferencia de calor es unidimensional
𝑑 𝑑𝑇
(𝑘 ) = 0
𝑑𝑥 𝑑𝑥

EJERCICIO

La distribución de temperatura a través de una pared de 1m de espesor en cierto instante está

dada como 𝑇(𝑥) = 𝑎 + 𝑏𝑥 + 𝑐𝑥 2 , donde T en °C y x en m. a=900°C, b=-300°C, c=-50°C/𝑚2 . Una
generación de calor uniforme g=1000 W/𝑚3 , está presente en la pared de área 10𝑚2 cuya 𝜌 =
1000 𝑘𝑔/𝑚3 y si conductividad térmica es de k=40 W/m°K.

Determinar:

1. La rapidez de calor que entra en la pared (x=0) y sale de la pared (x=Lm)

2. La rapidez de cambio del almacenamiento de energía en la pared Ealm=??
3. Velocidad respecto al tiempo del cambio de temperatura en x=0; 0.25; 0.5 m

Datos:

T(x)= 𝑎 + 𝑏𝑥 + 𝑐𝑥 2

a=900°C

b=-300°C

c=-50°C/𝑚2

g=1000 W/𝑚3

A=10𝑚2

𝜌 = 1000 𝑘𝑔/𝑚3

k=40 W/m°K.

1:

T(x)= a + bx + cx 2
dT
= b + 2cx
dx
28

dT
qentra = qx (0) = −kA
dx x=0
qentra = −𝑘𝐴(𝑏 + 𝑠𝑐𝑥)𝑥=0
°𝐶 𝑊
𝑞𝑒𝑛𝑡𝑟𝑎 = −𝑘𝐴𝑏 = 300 ∗ 40 ∗ 10𝑚2
𝑚 𝑚°𝐾
𝒒𝒆𝒏𝒕𝒓𝒂 = 𝟏𝟐𝟎 𝒌𝑾

𝑑𝑇
𝑞𝑠𝑎𝑙𝑒 = 𝑞𝑥=𝐿 = −𝑘𝐴
𝑑𝑥 𝑥=𝐿
𝑞𝑠𝑎𝑙𝑒 = −𝑘𝐴(𝑏 + 2𝑐𝑥)𝑥=𝐿

𝑞𝑠𝑎𝑙𝑒 = −(𝑏 + 2𝑐𝐿)𝑘𝐴

°𝐶 𝐶 𝐶 𝑊
𝑞𝑠𝑎𝑙𝑒 = − [300 + 2 (−50 2 ) ∗ 1𝑚 + 2 (−50 2 ) ∗ 1] ∗ 40 ∗ 10𝑚2
𝑚 𝑚 𝑚 𝑚°𝐾
𝒒𝒔𝒂𝒍𝒆 = 𝟏𝟔𝟎𝒌𝑾

2:
𝐸𝑒 + 𝐸𝑔 − 𝐸𝑠 = 𝐸𝑎𝑙𝑚

𝐸̇𝑔 = 𝑞𝑔̇ ∗ 𝐴 ∗ 𝐿
̇ = 𝑞𝑒𝑛𝑡𝑟𝑎 + 𝑞𝑔𝑒𝑛 − 𝑞𝑠𝑎𝑙𝑒
𝐸𝑎𝑙𝑚
𝑊
̇ = 120𝑘𝑤 + 1000
𝐸𝑎𝑙𝑚 ∗ 10𝑚2 ∗ 1𝑚 − 160𝑘𝑊
𝑚2
̇ = −𝟑𝟎𝒌𝑾
𝑬𝒂𝒍𝒎
Ejercicios: 1.9/1.11/1.17/1.22/1.23/2.4/2.21/2.23

Ecuación de Distribución de temperatura

(𝑻𝟏 − 𝑻𝟐)
𝑻=− 𝒙 + 𝑻𝟏
𝑳
Ecuación de Fourier
𝒅𝑻
∫ = 𝑪𝟏
𝒅𝒙
𝑇2 − 𝑇1
𝐶1 =
𝐿
𝒅𝑻
𝑸 = −𝒌𝑨
𝒅𝒙
Ecuación de Fourier
29

𝑻𝟐 − 𝑻𝟏
𝑸 = −𝒌𝑨[− ]
𝑳
𝒌𝑨
𝑸=− 𝑳
(𝑻𝟏 − 𝑻𝟐)Calor constante no depende de x

𝑻𝟏 − 𝑻𝟐
𝑸=
(𝑳⁄𝒌𝑨)

Analogía Eléctrica

R: resistencia térmica a la conducción

𝑇1 − 𝑇2
𝑞=
𝑅

R: resistencia térmica a la convección

1
R=ℎ𝐴
𝑐

𝑇1 − 𝑇2
𝑞=
𝑅

Pared compuesta
30

𝑇1 − 𝑇2 𝑇2 − 𝑇3 𝑇3 − 𝑇4
1𝑞2 = 2𝑞3 = 3𝑞4 =
𝑅1 𝑅2 𝑅3

1𝑄2 = 2𝑄3 = 3𝑄4

𝑇1 − 𝑇4 𝑇1 − 𝑇4
𝑞= →𝑞=
∑𝑅 𝑅1 + 𝑅2 + 𝑅3
31

Convección

𝑇∞, 2
𝑇∞,1 𝑇1 𝑇2 𝑇3 𝑇4

R= (L1/k1A) ∞
∞ ∞ ∞ ∞ ∞
𝑞𝑐 = ℎ𝐴𝑐 ∆𝑇
∞
∞ ∞ ∞ ∆𝑇 𝑇∞ − 𝑇 ∞ 𝑇∞ − 𝑇𝑐 ∞
𝑞𝑐 = = = Escriba aquí la ecu
1 1
Escriba aquí la ecuación. (
Escriba aquí laEscriba
ecuación. )Escriba𝑅aquí la ecuación.
)aquí (la ecuación. Escriba aquí la ecuación.
ℎ𝐴 𝑐 ℎ𝐴 𝑐

𝑇∞ 𝑇1

1
𝑅𝑐 =
ℎ𝐴𝑐
∞ 𝑇∞1 − 𝑇∞2
∞
𝑞=
𝑅 + 𝑅1 + 𝑅2 + 𝑅3 + 𝑅𝑐2
∞ 𝑐1 ∞
Escriba aquí la ecuación.
Escriba aquí la ecuación.

Ejercicio 3.11

Se desea mantener a 5°C el interior de un refrigerador cuyas dimensiones en la base son de

45cm * 45cm y la altura de 1,2m, las paredes del refrigerador están construidas por dos láminas
de acero de 0,318 cm de espesor con 5cm de aislante de fibra de vidrio entre ellas. Si la
temperatura ambiente en la cocina es de 30°C determinar:

El flujo de calor que debe extraer el refrigerador para mantener las condiciones especificadas,
determine además las temperaturas en sus paredes interiores.

𝑊
Fibra de vidrio ℎ1 = 10 𝑚2 °𝐶
𝑊
ℎ2 = 15 𝑚2 °𝐶
𝑊
𝑘𝑎𝑐𝑒𝑟𝑜 = 40 𝑚°𝐶
32

𝑊
acero acero 𝑘𝑣𝑖𝑑𝑟𝑖𝑜 = 0.04 𝑚°𝐶

𝑇∞2 − 𝑇∞1
𝑄=
𝑅𝑐1 + 𝑅1 + 𝑅2 + 𝑅3 + 𝑅𝑐2

Convección
1
𝑅𝑐1 =
ℎ1𝐴𝑐
1
𝑅𝑐1 =
10 ∗ 0.45 ∗ 1.2
°𝐶
𝑅𝑐1 = 0.185
𝑊
Conducción
𝐿1
𝑅1 =
𝑘1𝐴
0.00318 𝑚
𝑅1 =
𝑊
40 ∗ 0.45 ∗ 1.2 𝑚°𝐶 ∗ 𝑚2

°𝐶
𝑅1 = 0.00015
𝑊
𝐿2
𝑅2 =
𝑘2𝐴
0.05 𝑚
𝑅2 =
𝑊
0.04 ∗ 0.45 ∗ 1.2 𝑚°𝐶 ∗ 𝑚2

°𝐶
𝑅2 = 2.31
𝑊
𝑅1 = 𝑅3
33

30°𝐶 − 5°𝐶
𝑄=
°𝐶
(0.185 + 0.0015 + 2.31 + 0.015 + 0.1235) 𝑊

𝑄 = 9.55 𝑤𝑎𝑡𝑡𝑠 → 1 𝑝𝑎𝑟𝑒𝑑

4 paredes

𝑄𝑇 = 9.55 ∗ 4 = 38.19 𝑊 calor que debe remover la unidad de refrigeración

Temperaturas
𝑇∞2 − 𝑇4
𝑄=
𝑅𝑐2
𝑇4 = −𝑄𝑅𝑐2 + 𝑇∞2
𝑻𝟒 = 30°𝐶 − [9.55 ∗ 0.1235] = 𝟐𝟖. 𝟖𝟐°𝑪

𝑇4 − 𝑇3
𝑄=
𝑅3
𝑇3 = −𝑄𝑅3 + 𝑇4
𝑻𝟑 = 28.82°𝐶 − [9.55 ∗ 0.00015] = 𝟐𝟖. 𝟖𝟏°𝑪

𝑇3 − 𝑇2
𝑄=
𝑅2
𝑇2 = −𝑄𝑅2 + 𝑇3
𝑻𝟐 = 28.81°𝐶 − [9.55 ∗ 2.31] = 𝟔. 𝟕𝟓°𝑪

𝑇2 − 𝑇1
𝑄=
𝑅1
𝑇2 = −𝑄𝑅1 + 𝑇2
𝑻𝟐 = 6.75°𝐶 − [9.55 ∗ 0.0015] = 𝟔. 𝟕𝟒°𝑪

Ejercicio

La pared plana de un tanque con área de 5 m2 está cubierta con un aislamiento térmico de 2
capas. La pared del tanque es de acero con espesor L1=8mm y R1=46,5 W/m℃, la primera capa
de aislamiento en contacto con el acero está hecha de un refractario (NOVOASBOZURITA) con
34

espesor L2=50mm y R2= 0,144 W/m℃, la segunda capa de aislamiento mas externo es un
enlucido con espesor L3=10mm y R3= 0,698 W/m℃. Si la temperatura de la pared interna de
tanque es T1=250℃ y de la superficie más externa del aislamiento T4=50℃. Determinar la
cantidad de calor que se transmite a través de las paredes del tanque y las temperaturas internas
en las superficies de contacto de los aislamientos.

Datos

L1=8mm R1=46,5

L2=50mm R2=0,144

L3=10mm R3=0,698

T1=250℃

T4=50℃

Circuito Q

T1 T2 T3 T4

R1 R2 R3

𝐿1 0,008
𝑅1 = 𝑅1 =
𝑅1 𝐴 (46,5)(5)
𝐿2
𝑅2 =
𝑅2 𝐴

𝑅1 = 0,000034 ℃/𝑊
35

0,05 0,01
𝑅2 = 𝑅3 =
(0,144)(5) (0,698)(5)
𝐿3
𝑅3 =
𝑅3 𝐴

𝑅2 = 0,069 ℃/𝑊 𝑅3 = 0,00286 ℃/𝑊

𝑇1 − 𝑇4
𝑞𝑇 =
𝑅1 + 𝑅2 + 𝑅3

𝑞𝑇 = 2781,87

𝑇1 − 𝑇2 𝑇2 = 𝑇1 − 𝑄𝑅1
1𝑄2 =
𝑅1

𝑇2 = 250 − (2920,95)𝑊 ∗ 𝑅1

𝑇2 = 249,9 ℃

𝑇2 − 𝑇3 𝑇3 = 𝑇2 − 𝑄𝑅2
2𝑄3 =
𝑅2

𝑇3 = 249,9℃ − (2920,95)𝑊 ∗ 𝑅2

𝑇2 = 57,956 ℃
36

EL CILINDRO

Para condiciones de EE sin generación de calor la forma apropiada de la ecuación de calor es:
1𝑑 𝑑𝑇
(𝑘𝑟 ) = 0
𝑟 𝑑𝑟 𝑑𝑟

Consideración:

Donde k es una variable por el momento.

La rapidez a la que se conduce la energía a través de cualquier superficie cilíndrica se expresa:

𝑑𝑇 𝑑𝑇
𝑞𝑟 = −𝑘𝐴 = −𝑘(2𝜋𝑟𝐿)
𝑑𝑟 𝑑𝑟

𝑑𝑇
en 1 [𝑘𝑟 𝑑𝑟 ] es independiente de r

Ecuación de Fourier
37

𝑑2 𝑇 1 𝑑𝑇
+ =0
𝑑𝑟 2 𝑟 𝑑𝑟

𝑑𝑇
𝑝=
𝑑𝑟

𝑑𝑝 1
+ 𝑝=0
𝑑𝑟 𝑟

𝑑𝑝 1
∫ = − ∫ 𝑑𝑟
𝑑𝑟 𝑟

ln 𝑝 = − ln 𝑟 + ln(𝐶1 )

𝐶1
ln 𝑝 = ln ( )
𝑟

𝐶1
𝑝=
𝑟

𝑑𝑇 𝐶1
=
𝑑𝑟 𝑟

𝐶1
∫ 𝑑𝑇 = ∫ 𝑑𝑟 ⇒ 𝑇 = 𝐶1 ln 𝑟 + 𝐶2
𝑟

Distribución de temperatura

Fourier
Suponemos 𝑇(𝑟1) = 𝑇𝑆,1 y 𝑇(𝑟2) = 𝑇𝑆,2 (Condiciones
de frontera)

𝑇𝑆,1 = 𝐶1 ln 𝑟1 + 𝐶2 (3)

𝑇𝑆,2 = 𝐶1 ln 𝑟2 + 𝐶2 (4)

(3) – (4)
38

(𝑇𝑆,1 − 𝑇𝑆,2 ) = 𝐶1 ln 𝑟1 − 𝐶1 ln 𝑟2

𝑟1
(𝑇𝑆,1 − 𝑇𝑆,2 ) = 𝐶1 ln ( )
𝑟2

𝑟2
(𝑇𝑆,1 − 𝑇𝑆,2 ) = − 𝐶1 ln ( )
𝑟1

(𝑇𝑆,1 − 𝑇𝑆,2 )
𝐶1 = − 𝑟 (5)
ln 𝑟2
1

(5) en (3)

(𝑇𝑆,1 − 𝑇𝑆,2 )
𝑇𝑆,1 = − 𝑟 ln 𝑟1 + 𝐶2
ln 2
𝑟1

(𝑇𝑆,1 − 𝑇𝑆,2 )
𝐶2 = 𝑇𝑆,1 + 𝑟 (ln 𝑟1 ) (6)
ln 𝑟2
1

(5) y (6) en la distribución de temperatura 𝑇 = 𝐶1 ln 𝑟 + 𝐶2

(𝑇𝑆,1 − 𝑇𝑆,2 ) (𝑇𝑆,1 − 𝑇𝑆,2 )

𝑇= − 𝑟2 ln 𝑟 + 𝑇𝑆,1 + 𝑟 ln 𝑟1
ln ln 2
𝑟1 𝑟1

(𝑇𝑆,1 − 𝑇𝑆,2 ) 𝑟
𝑇(𝑟) = 𝑇𝑆,1 − 𝑟2 ∗ ln ( ) (7)
ln 𝑟 𝑟1
1
39

(𝑇𝑆,1 − 𝑇𝑆,2 ) 𝑟
𝑇(𝑟) = 𝑇𝑆,2 − 𝑟2 ∗ ln ( ) (8)
ln 𝑟 𝑟2
1

𝑑𝑇
𝑞"𝑟 = −𝑘 𝑠𝑖 𝐴𝑟 = 2𝜋𝑟𝐿
𝑑𝑟

𝑑𝑇
𝑞𝑟 = −𝑘 𝐴(𝑟)
𝑑𝑟

𝑑𝑇 (𝑇𝑆,1 − 𝑇𝑆,2 ) 1
=− 𝑟 ∗
𝑑𝑟 ln 𝑟2 𝑟
1

(𝑇𝑆,1 − 𝑇𝑆,2 ) 1
𝑞𝑟 = 𝑘 (2𝜋𝑟𝐿) [ 𝑟 ∗ ]
ln 𝑟2 𝑟
1

(𝑇𝑆,1 − 𝑇𝑆,2 )
𝑞𝑟 = 2𝜋𝑘𝐿 ∗ [ 𝑟 ] (9)
ln( 2)
𝑟1
40

𝑇∞,1 − 𝑇∞,4
𝑞𝑟 =
𝑅𝑐1 + 𝑅1 + 𝑅2 + 𝑅3 + 𝑅𝑐2

Coeficiente global de T.C : U

𝑞 = 𝑈 ∗ 𝐴 ∗ ∆𝑇

𝐴 = 2𝜋𝑟1 𝐿

1
𝑈=
∑ 𝑅𝑇
41

1
𝑈=
1 𝑟 𝑟 𝑟 𝑟 𝑟 𝑟 𝑟 1
+ 1 ln (𝑟2 ) + 𝐾1 ln (𝑟3 ) + 𝐾1 ln (𝑟4 ) + 𝑟1
ℎ1 𝐾𝐴 1 𝐵 2 𝐶 3 4 ℎ4

Ejemplo:

La posible existencia de un espesor de aislamiento óptimo para sistemas radiales lo sugiere la

presencia de efectos que compiten asociados con un aumento en este espesor.

En particular, aunque la resistencia de conducción aumenta al agregar un aislante, la resistencia

de convección disminuye debido al aumento del área de la superficie exterior.

Por ello puede existir un espesor de aislamiento que minimice la perdida de calor al maximizar
la resistencia total a la transferencia de calor. Resuelva este problema considerando el siguiente
sistema:

1. Un tubo de cobre con pared delgada de radio 𝑟𝑖 se usa para transportar un fluido
refrigerante de baja temperatura y está a una temperatura 𝑇𝑖 , que es menor que la del
aire del medio a 𝑇∞ alrededor del tubo. ¿Hay un espesor óptimo asociado con la
aplicación de aislante al tubo?

𝑊
𝑘 = 0,055
𝑚°𝐾

𝑟
ln( )
𝑟𝑖 1 R térmica total
𝑅𝑇𝑜𝑡 = 2𝜋𝐾
+ 2𝜋𝑟ℎ
42

Transferencia de
𝑇∞ − 𝑇𝑖 calor por unidad de
𝑞′ =
𝑅𝑇𝑜𝑡 L

𝑑𝑅𝑇𝑜𝑡
=0
𝑑𝑟

1 1 1
∗ 𝑟 −
𝑟𝑖 ( ) 2𝜋𝑘 2𝜋𝑟 2 ℎ
𝑟𝑖

1 1
− =0
2𝜋𝑟𝑘 2𝜋𝑟 2 ℎ

𝑘
𝑟=
ℎ

𝑊 𝑊
con ℎ = 5 𝑚2 𝐾
y 𝑘 = 0,055 𝑚𝐾

calcule el radio crítico

0,055
𝑟𝑐𝑟 = = 0,011𝑚
5

Ejercicio 3.49 pág. 143

El vapor que fluye a través de un tubo largo de pared delgada mantiene la pared del tubo a una
temperatura uniforme de 500 °K. El tubo está cubierto con una manta aislante compuesta con
dos materiales diferentes A y B.

Se supone que la interfaz entre los dos materiales tiene una resistencia de contacto infinita, y
que toda la superficie externa está expuesta al aire, para el cual 𝑇∞ = 300 °K y h = 25 W/𝑚2 °K.
43

a) Dibuje el circuito térmico del sistema. Usando los símbolos precedentes, marque todos
los nodos y resistencias pertinentes.

b) Para las condiciones que se establecen, ¿Cuál es la pérdida total de calor del tubo?
¿Cuáles son las temperaturas de la superficie externa Ts,2(A) y Ts,2(B)?

1
𝑅𝑐𝑜𝑛𝑣(𝐴) = 𝑅𝑐𝑜𝑛𝑣(𝐵) =
𝜋𝑟2 ℎ

𝑟
ln (𝑟2 )
1
𝑅𝑐𝑜𝑛𝑑(𝐴) =
𝜋𝑘𝐴

𝑟
ln (𝑟2 )
1
𝑅𝑐𝑜𝑛𝑑(𝐵) =
𝜋𝑘𝐵

𝑞𝑇𝑜𝑡𝑎𝑙 = 𝑞𝐴 + 𝑞𝐵
44

𝑚°𝐾 𝑚°𝐾
𝑅𝑐𝑜𝑛𝑣 = [𝜋 ∗ 0,1 ∗ 25]−1 [ ] = 0,1273
𝑊 𝑊

0,1
ln 𝑚°𝐾
0,05
𝑅𝑐𝑜𝑛𝑑𝐴 = = 0,1103
𝜋∗2 𝑊

0,1
ln 𝑚°𝐾
0,05
𝑅𝑐𝑜𝑛𝑑𝐵 = = 0,88254
𝜋 ∗ 0,25 𝑊

𝑇𝑠,1 − 𝑇∞ 𝑇𝑠,1 − 𝑇∞
𝑞𝑇𝑜𝑡𝑎𝑙 = +
𝑅𝑐𝑜𝑛𝑑(𝐴) + 𝑅𝑐𝑜𝑛𝑣(𝐴) 𝑅𝑐𝑜𝑛𝑑(𝐴) + 𝑅𝑐𝑜𝑛𝑣(𝐵)

(500 − 300) (500 − 300)

𝑞𝑇𝑜𝑡𝑎𝑙 = +
(0,1103 + 0,1273) (0,8825 + 0,1273)

𝑊
𝑞𝑇𝑜𝑡𝑎𝑙 = 841,75 + 198,05 = 1039,8
𝑚

𝑇𝑠,1 − 𝑇𝑠,2 𝑞𝐴 𝑅𝑐𝑜𝑛𝑑(𝐴) = 𝑇𝑠,1 − 𝑇𝑠,2

𝑞𝐴 =
𝑅𝑐𝑜𝑛𝑑𝐴

𝑇𝑠,2 = 𝑇𝑠,1 − 𝑞𝐴 𝑅𝑐𝑜𝑛𝑑(𝐴)

𝑇𝑠,2 = 500 − 841,75 ∗ 0,1103

𝑇𝑠,2 = 407,15 °𝐾
45

𝑇𝑠,1 − 𝑇𝑠,2 𝑞𝐵 𝑅𝑐𝑜𝑛𝑑(𝐵) = 𝑇𝑠,1 − 𝑇𝑠,2

𝑞𝐵 =
𝑅𝑐𝑜𝑛𝑑𝐵

𝑇𝑠,2 = 𝑇𝑠,1 − 𝑞𝐵 𝑅𝑐𝑜𝑛𝑑(𝐵)

𝑇𝑠,2 = 500 − 198,05 ∗ 0,8825

𝑇𝑠,2 = 325,22 °𝐾

𝑇𝑠,1 − 𝑇∞
𝑞𝑇𝑜𝑡𝑎𝑙 =
𝑅𝑒𝑞𝑢𝑖𝑣𝑎𝑙𝑒𝑛𝑡𝑒

Ejercicios 3.54 pág. 144

Una esfera hueca de aluminio, con un calentador eléctrico en el centro, se utiliza en pruebas
para determinar la conductividad térmica de materiales aislantes. Los radios interior y exterior
de la esfera son 0.15 y 0.18 m respectivamente, y la prueba se hace en condiciones de estado
estable, en las que la superficie interna del aluminio se mantiene a 250°C. En una prueba
particular, una capa esférica de aislante se funde sobre la superficie externa de la esfera y
alcanza un espesor de 0 .12 m.

El sistema está en un cuarto para el que la temperatura del aire es 20°C, y el coeficiente de
convección en la superficie externa del aislante es 30 W/m2 °K. Si se disipan 80 W por el
calentador bajo condiciones de estado estable. ¿Cuál es la conductividad térmica del aislante?

Continuación del ejercicio de la hoja 39-A

T1 𝑅𝐶𝑜𝑛𝑑 𝑎𝑖𝑠𝑙 𝑇∞

𝑅𝐶𝑜𝑛𝑑 𝐴𝑙 𝑅𝐶𝑜𝑛𝑣 q = 80W

1 1
−
𝑟1 𝑟2
𝑅𝑐𝑜𝑛𝑑 𝐴𝑙 = 1 1
4𝜋𝑘𝐴𝑙
𝑟2 − 𝑟3
𝑅𝑐𝑜𝑛𝑑 𝑎𝑖𝑠𝑙 =
4𝜋𝑘𝑎𝑖𝑠𝑙
46

1
𝑅𝑐𝑜𝑛𝑣 =
ℎ4𝜋𝑟32

𝑇1 − 𝑇∞
𝑞𝑇𝑜𝑡𝑎𝑙 = = 80
𝑅𝑇𝑜𝑡𝑎𝑙

250 − 20
80 =
1 1 1 1
( − ) − 1
0,15 0,18 0,18 0,30
+ 4𝜋𝑅 + 2
4𝜋(230) 𝑎𝑖𝑠𝑙 30(4𝜋)(0,3)

0,177 230
3,84 ∗ 10−4 + + 0,029 = = 2,875
𝑅𝑎𝑖𝑠𝑙 80

𝑊
𝑅𝑎𝑖𝑠𝑙 = 0,062
𝑚°𝐾

LA ESFERA

Consideremos una esfera hueca para analizar la conducción

Volumen diferencia de control

La conservación de energía es 𝑞𝑟 = 𝑞𝑟 + 𝑑𝑟

Para condiciones unidimensionales de Estado Estable sin generación interna de calor:

Fourier:
𝑑𝑇
𝑞𝑟 = −𝑘𝐴
𝑑𝑟
𝐴 = 4𝜋𝑟 2
𝑑𝑇
𝑞𝑟 = −𝑘4𝜋𝑟 2
𝑑𝑟
Como 𝑞𝑟 es una constante independiente de r
47

𝑇𝑠,2
𝑞𝑟 𝑟2 𝑑𝑟
∫ = ∫ 𝑘(𝑇)𝑑𝑇
4𝜋 𝑟1 𝑟 2
𝑇𝑠,1

Como k es constante

Ejercicios propuestos 3.40/3.46/3.59

Ejemplos 3.4/3.3

Ejemplo 3.5

Un contenedor metálico esférico de pared delgada se utiliza para almacenar nitrógeno líquido a
77°K. el contenedor tiene un diámetro de 0.5m y está cubierto de un aislante reflector al vacío
compuesto de polvo de dióxido de silicio. El aislante tiene un espesor de 25mm y la superficie
externa se expone al aire del ambiente a 300°K. Se sabe que el coeficiente de convección es 20
𝑊⁄ . La entalpia de vaporización y la densidad del nitrógeno líquido son 2∗ 105 𝐽/𝐾𝑔 y 804
𝑚2 °𝐾
kg/𝑚3 respectivamente.

1. Cuál es la transformación de calor al nitrógeno liquido

2. Cuál es la velocidad de evaporación del nitrógeno

Propiedades: polvo de dióxido de silicio al vacío (300°K) k=0.0017 W/m°K

Análisis
48

El balance de energía para una superficie de control alrededor del nitrógeno

Se asocia con la perdida de energía latente debido a la evaporación

Pérdida por día

49

Volumétrica:

Observaciones:

Conducción con generación de energía

La sección anterior se analizó la distribución de temperatura en un medio; mediante condiciones

de frontera.

Se analizaran distribuciones de temperatura de procesos que pueden ocurrir dentro del medio.
En particular se analizara situaciones para las que la energía térmica se genera debido a la
conversión de alguna otra fuente de energía.

Un proceso común es la conversión de energía eléctrica a térmica en un medio conductor de

corriente.

La razón a al que se genera E al pasar una corriente I a través de un medio de

re resistencia eléctrica Re.

Si la generación de potencia (W) ocurre de manera uniforme a lo largo del medio de volumen V
la razón de generación volumétrica es 𝑊/𝑚3

La Pared Plana

Consideraciones:

Generación de E uniforme por unidad de volumen (𝑞̇ 𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒) y las superficies se mantienen
a Ts,1 y Ts,2 para una k=constante.

a) Condiciones de frontera asimétricas.

50

b) Condiciones de frontera simétricas

𝑑𝑇
−𝑘 = ℎ(𝑇𝑠 − 𝑇∞ )
𝑑𝑥 𝑥=𝐿

𝑞̇ 𝐿 = ℎ(𝑇𝑠 − 𝑇∞ )

c) Superficie adiabática en el plano medio

𝑞̇ 𝐿
𝑇𝑠 = 𝑇∞ + ℎ

Por las condiciones de frontera se establece

51

Distribución de temperatura flujo de calor en cualquier punto de la pared

Con generación de 𝑞̇ el flujo de calor ya no es independiente de x

Si ambas superficies se mantienen a una misma temperatura común

La temperatura máxima se tiene en el plano medio

Temperatura adicional generada por 𝐸̇𝑔 o 𝑞̇

𝑞̇ 𝐿2
𝑇𝑔 =
2𝑘
Para este caso la distribución de temperatura es

Para el plano simétrico el gradiente de temperatura es cero

𝑑𝑇
(𝑑𝑥) =0 No existe transferencia de calor a través del
𝑥=0
plano y se representa con la superficie
adiabática c)

Pared con generación de energía considerando convección

52

𝑇∞ 𝑇∞

x=-L x=0 x=L

∞ ∞

∞ 𝑑𝑇
−𝑘𝐴 = ℎ𝐴(𝑇𝑥=𝐿 − 𝑇∞) ∞
𝑑𝑥
Escriba aquí la ecuación. Escriba aquí la ecuación.

ℎ𝐿
𝐵𝑖 =
𝑘
Número de Biot

Una pared plana se compone de dos materiales Ay B. la pared del material A tiene una
generación de calor uniforme 𝑞̇ = 1.5 ∗ 106 𝑊⁄ 3 , 𝑘𝐴 = 75 𝑊⁄𝑚°𝐾 y un espesor 𝐿𝐴 =
𝑚
50𝑚𝑚. El material b de la pared no tiene generación y su 𝑘𝐵 = 150 𝑊⁄𝑚°𝐾

y un espesor 𝐿𝐵 = 20𝑚𝑚. La superficie interior del material A esta bien aislada, mientras que
la superficie exterior del material B se enfría con un flujo de agua con 𝑇∞ = 30°𝐶 y ℎ =
1000 𝑊⁄ 2 .
𝑚 °𝐾
1. Dibuje la distribución de temperatura que existe en el compuesto bajo condiciones de EE.
2. Determinar la temperatura To de la superficie aislada y la temperatura T2 de la superficie
enfriada.

1)
53

b) pendiente cero
a) parabólica
𝑅 150
d) cambio de pendiente 𝑅𝐵 = 75
=2
𝐴
c) pendiente lineal
e) el gradiente para el agua es grande

Balance de energia

𝐸̇𝑒𝑛𝑡𝑟𝑎 + 𝐸̇𝑔 − 𝐸̇𝑠𝑎𝑙𝑒 = 𝐸̇𝑎𝑙𝑚

𝐸̇𝑒𝑛𝑡𝑟𝑎 = 𝐸̇𝑠𝑎𝑙𝑒

Para relacionar T2
No existe generación de calor en B
Tenemos Estado Estable y para la superficie de control área unitaria el flujo de calor hacia
el material x=𝐿𝐴 es igual al flujo de calor desde el material en x=𝐿𝐴 + 𝐿𝐵 debido a la
convección

Volumen de control sobre el material A

𝐸̇𝑒𝑛𝑡𝑟𝑎 + 𝐸̇𝑔 − 𝐸̇𝑠𝑎𝑙𝑒 = 𝐸̇𝑎𝑙𝑚
𝐸̇𝑔 = 𝐸̇𝑠𝑎𝑙𝑒
Para una area superficial unitaria
𝑞̇ 𝐿 = ℎ(𝑇𝑠 − 𝑇∞ )
𝑞̇ 𝐿𝐴 = ℎ(𝑇2 − 𝑇∞ ) → 𝑇2

Como To es la temperatura máxima en el plano medio

Aplicamos analogía eléctrica en B porque el material A esta aislado

54

𝑇1 − 𝑇∞
𝑞¨ =
𝑅𝑒𝑞𝑢𝑖𝑣
Para un área superficial

𝐿
𝑅𝑐𝑜𝑛𝑑 = 𝑅𝐵
𝐵
1
𝑅𝑐𝑜𝑛𝑣 =
ℎ
55

Sistemas Radiales

La razón a la que se que genera

calor dentro del cilindro debe ser
igual a la rapidez con que se igual
a la rapidez con que se transmite
calor por convección de la
superficie del cilindro a un fluido
en movimiento.

Esta condición permite que la

temperatura de la superficie se mantenga en un valor 𝑇𝑆
1𝑑 𝑑𝑇 𝑞̇
(𝑟 ) + = 0 (1)
𝑟 𝑑𝑟 𝑑𝑟 𝑘
Se integra separando variables y suponiendo generación uniforme
𝑑𝑇 𝑞̇
𝑟 = − 𝑟 2 + 𝐶1 (2)
𝑑𝑟 2𝑘
Repitiendo el procedimiento: Solución general para la distribución de temperatura
𝑞̇ 2
𝑇(𝑟) = − 𝑟 + 𝐶1 ln 𝑟 + 𝐶2 (3)
4𝑘
Para 𝐶1 𝑦 𝐶2 se aplica las condiciones de frontera
𝑑𝑇
| =0 𝑦 𝑇(𝑟0 ) = 𝑇𝑆
𝑑𝑟 𝑟=0
En condiciones de frontera simétrica estamos ubicados en el plano medio de una pared por lo
que

𝑟 = 0 ⇒ 𝐶1 = 0
𝑞̇ 2
𝐶2 = 𝑇𝑠 + 𝑟 ⇒ 𝑇𝐶
4𝑘 0
Se determina entonces que la distribución de temperatura es: Se puede evaluar la T-C
𝑞̇ 𝑟02 𝑟2 para cualquier radio r en el
𝑇(𝑟) = (1 − ) + 𝑇𝑠 cilindro.
4𝑘 𝑟02

Evaluando la temperatura en la línea central se obtiene:

To es la
temperatura de la
𝑞̇ 𝑟12 línea central
𝑇(0) = 𝑇0 = + 𝑇1
4𝑘
𝑇(𝑟) − 𝑇𝑠 𝑟 2
=1−( )
𝑇0 − 𝑇𝑠 𝑟0
Para relacionar la temperatura de la superficie 𝑇𝑠 , con la temperatura del fluido frio 𝑇∞ , se usa
un balance de energía en la superficie o un balance global de energía.
56

𝑞̇ (𝜋𝑟02 𝐿) = ℎ(2𝜋𝑟0 𝐿)(𝑇𝑠 − 𝑇∞ )

𝑞̇ 𝑟
𝑇𝑠 = 𝑇∞ +
2ℎ

Factor de Forma del Calor para varias geometrías (Conducción)

Pared plana
𝐾𝐴
𝑞= (𝑇 − 𝑇2 )
𝐿 1
𝐴
𝐹=
𝐿

𝑞 = 𝐾𝐹(𝑇𝑠,1 − 𝑇𝑠,2 )

Cilindro

𝑇𝑠,1 − 𝑇𝑠,2
𝑞 = 2𝜋𝐾𝐿 [ 𝑟 ]
ln (𝑟2 )
1

2𝜋𝐿
𝐹= 𝑟
ln (𝑟2 )
1

𝑞 = 𝐾𝐹(𝑇𝑠,1 − 𝑇𝑠,2 )

Esfera
4𝜋
𝐹=
1 1
𝑟1 − 𝑟2

𝑞 = 𝐾𝐹(𝑇𝑠,1 − 𝑇𝑠,2 )
57

Otras geometrías

Cilindros Excéntricos
2𝜋𝐿
𝐹= 1+𝜌2 −𝜀2
𝑐𝑜𝑠ℎ −1 ( )
2𝜌

𝑟 𝑒
𝜌 = 𝑟1 𝜀=𝑟
2 2

𝑞 = 𝐾𝐹(𝑇𝑠,1 − 𝑇𝑠,2 )

Cilindro circular en un ducto hexagonal

2𝜋𝐿
𝐹= 𝑟2
ln ( ) − 0,10669
𝑟1

𝑞 = 𝐾𝐹 (𝑇𝑠,1 − 𝑇𝑠,2 )

Cilindro circular en un ducto cuadrado

2𝜋𝐿
𝐹= 𝑟2
ln ( ) − 0,27079
𝑟1

𝑞 = 𝐾𝐹 (𝑇𝑠,1 − 𝑇𝑠,2 )
58

Cilindro infinito en un medio semi infinito

2𝜋𝐿
𝐹=
𝑙
𝑐𝑜𝑠ℎ−1 (𝑟)

𝑞 = 𝐾𝐹 (𝑇𝑠,1 − 𝑇𝑠,2 )

Resistencia Térmica a la Radiación

1
𝑅=
ℎ𝐴𝐶

𝐴𝐶 = 𝐴𝑅

ℎ ⇒ 𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛

1 1 1 1
𝑅𝑅 = = +
ℎ𝑅 𝐴𝑅 𝑅𝑒𝑞 𝑅𝑅 𝑅𝐶

1
𝑅𝑒𝑞 =
ℎ𝑅 𝐴 + ℎ𝐴
59

𝑇1 − 𝑇0
𝑞=
𝐿
(𝑅𝐴) + 𝑅𝑒𝑞

𝑞𝑟 = 𝜀𝐴𝜎(𝑇14 − 𝑇24 )

𝑞𝑟 = 𝜀𝐴𝜎(𝑇12 + 𝑇22 )(𝑇1 + 𝑇2 )(𝑇1 − 𝑇2 )

𝑇1 + 𝑇2 𝑇13 + 3𝑇12 𝑇2 + 3𝑇1 𝑇22 + 𝑇23

𝑇𝑚 = ≡ 4𝑇𝑚3 =
2 2

(𝑇12 + 𝑇22 )(𝑇1 + 𝑇2 ) = 4𝑇𝑚3

𝑞𝑟 = 𝜀𝐴𝜎4𝑇𝑚3 (𝑇1 − 𝑇2 )

ℎ𝑟 = 𝜀𝜎𝑇𝑚3

𝑇1 − 𝑇2
𝑞𝑟 =
1
( )
ℎ𝑅 ∗ 𝐴

𝑇2 + 𝑇0 3
ℎ𝑅 = 4𝜀𝜎 ( )
2

Por lo general se considera que ℎ𝐶 ≡ ℎ𝑅 .

Valores de emisividad
Superficie 𝜺
Aleación de Al 0,035
Al negro anodizado 0,80
Tierra 0,94
Pintura blanca acrílica 0,90
Pintura negra esmaltada 0,78
Asfalto 0,88
Concreto 0,90
Vidrio Pirex 0,80
Acero inoxidable 0,30
60

Conducción y convección combinadas en un elemento estructural.

Sistemas aleteados

Si 𝑇𝑆 es fija dos formas en las que es posible aumentar la transferencia de calor. El h podría
aumentarse incrementando la velocidad del fluido y/o podría reducirse la temperatura del fluido
𝑇∞ .
Existe una opción adecuada incrementando el área de la superficie a través de la cual ocurre la
convección. Esto se logra con el empleo de aletas que se extienden desde la pared al fluido
circundante.

La conductividad térmica del material de la aleta tiene fuerte efecto sobre la distribución de
temperatura a lo largo de la aleta y por tanto influye en el grado al que la transferencia de calor
aumenta.

Idealmente el material de la aleta debe tener una conductividad térmica grande para minimizar
variaciones de temperatura desde la base hasta la punta.
61

Formas de aletas (Configuraciones)

a) Aleta recta de sección transversal uniforme.

b) Aleta recta de sección transversal no uniforme.
c) Aleta anular.
d) Aleta de aguja.

Modelo matemático
Área superficial del
elemento diferencial

Área de la sección
transversal varia con x

𝑞𝑥 = 𝑞𝑥+𝑑𝑥 + 𝑑𝑞𝑐𝑜𝑛𝑣 (1) 𝐵𝑎𝑙𝑎𝑛𝑐𝑒 𝑑𝑒 𝑒𝑛𝑒𝑟𝑔í𝑎

62

𝑑𝑇
𝑞𝑥 = −𝐾𝐴𝐶 (2) 𝐹𝑜𝑢𝑟𝑖𝑒𝑟
𝑑𝑥
La conducción de calor en x + dx es:
𝑑𝑞𝑥
𝑞𝑥+𝑑𝑥 = 𝑞𝑥 + 𝑑𝑥 (3)
𝑑𝑥
𝑑𝑇 𝑑 𝑑𝑇
𝑞𝑥+𝑑𝑥 = −𝐾 𝐴𝐶 − 𝐾 (𝐴𝐶 )
𝑑𝑥 𝑑𝑥 𝑑𝑥

𝑑𝑞𝑐𝑜𝑛𝑣 = ℎ𝑑𝐴𝑆 (𝑇 − 𝑇∞ ) (4)

Fourier
dAs: Área superficial del
𝑑𝑇 elemento diferencial
1. 𝑞𝑥 = −𝑘𝐴𝐶
𝑑𝑥

2. 𝑑𝑞𝑐𝑜𝑛𝑣 = ℎ𝑑𝐴𝑆 (𝑇 − 𝑇∞ )

sustituyendo (3) y (4) en (1)

𝑑 𝑑𝑇 ℎ 𝑑𝐴𝑆
(𝐴𝐶 ) − (𝑇 − 𝑇∞ ) = 0 (5)
𝑑𝑥 𝑑𝑥 𝐾 𝑑𝑥

𝑑2 𝑇 1 𝑑𝐴𝐶 𝑑𝑇 1 ℎ 𝑑𝐴𝑆
3. +( ) −( ) (𝑇 − 𝑇∞ ) = 0 (6)
𝑑𝑥 2 𝐴𝐶 𝑑𝑥 𝑑𝑥 𝐴𝐶 𝐾 𝑑𝑥

Para resolver esta expresión acudimos primero a aletas simples  aletas rectangulares rectas.

Aletas de área de sección transversal uniforme

Cada aleta se une a una superficie de base de temperatura 𝑇(0) = 𝑇𝑏 y se extiende en un fluido
de temperatura 𝑇∞ .

𝑑𝐴𝑆
Para estas aletas 𝐴𝐶 = 𝑐𝑡𝑒 y 𝑑𝑥
= 𝑃𝑥 , donde 𝐴𝑆 es el área de la superficie de la base a x y P es
el perímetro.
𝑑𝐴𝑆 𝑑𝐴𝑆
⇒ =0 𝑦 =𝑃
𝑑𝑥 𝑑𝑥
Entonces (6) se reduce a:
63

𝑑2 𝑇 ℎ𝑝
− (𝑇 − 𝑇∞ ) = 0 (7)
𝑑𝑥 2 𝑘𝐴𝐶

Para simplificar la ecuación transformamos la variable dependiente definiendo:

Definimos un exceso de temperatura 𝜃 y transformamos la variable dependiente

𝜃(𝑥) ≡ 𝑇(𝑥) − 𝑇∞ (8)

como 𝑇∞ = 𝑐𝑡𝑒

𝑑𝜃 𝑑𝑇 𝑑2 𝜃 𝑑2 𝑇
= ; =
𝑑𝑥 𝑑𝑥 𝑑𝑥 2 𝑑𝑥 2
(8) en (7) y cambiando de variables dependientes

𝑑2 𝜃
− 𝑚2 𝜃 = 0 (9)
𝑑𝑥 2

de la (9) haciendo una analogía con (7)

ℎ𝑃
𝑚2 = (10)
𝑘𝐴𝐶

la ecuación (9) es una ecuación diferencial lineal de segundo orden homogénea con
coeficiente cte y su solución general es:

Solución de

𝑑2 𝜃
− 𝑚2 𝜃 = 0
𝑑𝑥 2
64

𝜃(𝑥) = 𝐶1 𝑒 𝑚𝑥 + 𝐶2 𝑒 −𝑚𝑥 (11)

Para conducir 𝐶1 𝑦 𝐶2 se aplican condiciones de frontera:

1. Condición en la base de la aleta (x=0)

𝜃0 = 𝑇𝑏 − 𝑇∞ ≡ 𝜃𝑏 (12)

2. Condición en el extremo de la aleta (x=L) implica cuatro situaciones física:

Caso A

Transferencia de calor por convección desde el extremo de la aleta

𝑑𝑇
ℎ𝐴𝐶 [𝑇(𝐿) − 𝑇∞ ] = −𝑘𝐴 | (13)
𝑑𝑥 𝑥=𝐿

𝑑𝜃
ℎ𝜃(𝐿) = −𝑘 | (14)
𝑑𝑥 𝑥=𝐿

La rapidez a la que la energía se transfiere haca el fluido por convección desde el extremo debe
ser igual a la rapidez a la que la energía alcanza el extremo por conducción a través de la aleta.

(11) en (12) y (14)

(7) en (8) y (10)

𝜃𝑏 = 𝐶1 + 𝐶2 (15)

ℎ(𝐶1 𝑒 𝑚𝐿 + 𝐶2 𝑒 −𝑚𝐿 ) + 𝑘𝑚(𝐶2 𝑒 .𝑚𝐿 − 𝐶1 𝑒 𝑚𝐿 ) (16)

ℎ
𝜃 cosh 𝑚 (𝐿 − 𝑥) + ( ) 𝑠𝑒𝑛 ℎ 𝑚 (𝐿 − 𝑥)
= 𝑚𝑘 (17)
𝜃𝑏 ℎ
cosh 𝑚𝐿 + ( ) 𝑠𝑒𝑛 𝑚𝐿
𝑚𝑘
65

Al aumentar x la gradiente de
temperatura disminuye.

Es consecuencia de la reducción de Transferencia de calor por conducción q(x) con el aumento

de x debido a las perdidas por convección continuas en la superficie de la aleta.

qp se puede evaluar utilizando distribución de temperaturas pero se puede usar Fourier en la

base de la aleta.

𝑑𝑇 𝑑𝜃
𝑞𝑓 = 𝑞𝑏 = −𝑘𝐴𝐶 | = −𝑘𝐴𝐶 | (18)
𝑑𝑥 𝑥=0 𝑑𝑥 𝑥=0

ℎ
𝑠𝑒𝑛ℎ 𝑚𝐿 + ( ) cosh 𝑚𝐿
𝑞𝑓 = √ℎ𝑃𝑘𝐴𝐶 𝜃𝑏 𝑚𝑘 (19)
ℎ
cosh 𝑚𝐿 + ( ) 𝑠𝑒𝑛ℎ 𝑚𝐿
𝑚𝑘

𝑇 = 𝑇∞ + (𝑇𝑏 − 𝑇∞ )𝑒 −𝑚𝑥 (20)

𝑄 = 𝐾𝐴𝐶 𝑚(𝑇𝑏 − 𝑇∞ )−𝑚𝑥 (21)

66

Ejemplo 3.8

Una varilla muy larga de 5 mm de diámetro tiene un extremo que se mantiene a I00°C.

La superficie de la varilla se expone al aire ambiente a 25 C con un coeficiente de transferencia

de calor por convección de 100 W/m2°K.

Determinar:

1. Las distribuciones de temperaturas a lo largo de varillas construidas de cobre puro,

aleación de aluminio 2024 y acero inoxidable tipo AISI 316. ¿Cuáles son las pérdidas de
calor correspondientes de las varillas?
2. Calcule el largo de las varillas para que la suposición de una longitud infinita de una
estimación exacta de la pérdida de calor.

1. Condiciones de estado estable.

2. Conducción unidimensional a lo largo de la varilla.
3. Propiedades constantes.
4. Intercambio de radiación insignificante con los alrededores.
5. Coeficiente convectivo constante y uniforme.
6. Varilla infinitamente larga.
67

Tabla A1. Incropera

Cobre
𝑇𝑏 + 𝑇∞
𝑇=( ) = 62,5℃ ≈ 335°𝐾
2
𝑊
𝑘 = 398
𝑚°𝐾

Aluminio (2024)(T=335°𝐾); k = 80 W/m°𝐾

Acero inoxidable AISI 316 (T=335°𝐾); k=14 W/m°𝐾

𝜃
= 𝑒 −𝑚𝑥
𝜃𝑏

𝜃 = 𝑇 − 𝑇∞

𝜃𝑏 = 𝑇𝑏 − 𝑇∞

𝑇 − 𝑇∞
= 𝑒 −𝑚𝑥
𝑇𝑏 − 𝑇∞

𝑇 = 𝑇∞ + ( 𝑇 − 𝑇∞ ) 𝑒 −𝑚𝑥

1
ℎ𝑃 2
𝑚=( )
𝑘𝐴𝐶

𝑃 = 𝜋𝐷

𝜋𝐷 2
𝐴𝐶 =
4

1
4ℎ 2
𝑚=( )
𝑘𝐷
68

𝑚 = 14,2 𝑚−1 𝐶𝑢

𝑚 = 21,2 𝑚−1 𝐴𝑙

𝑚 = 75,6 𝑚−1 𝐴𝑐𝑒𝑟𝑜 𝐴𝐼𝑆𝐼 316

𝑥 = (50 − 100 − 150 − 200 − 25 − 300

69

La pérdida de calor

𝑞𝑓 = √ℎ𝑃𝑘𝐴𝐶 𝜃𝑏

1
𝑊 𝑊 𝜋 2
𝑞𝑓 = [100 2 ∗ 𝜋 ∗ 0,005𝑚 ∗ 398 ∗ (0,005)2 𝑚2 ] (100 − 25)℃
𝑚 °𝐾 𝑚°𝐾 4

𝑞𝑓𝐶𝑢 = 8,3 𝑊

𝑞𝑓𝐴𝑙 = 5,6 𝑊

𝑞𝑓𝑎𝑐 = 1,6 𝑊

b) Como no hay pérdida de claro en el extremo de una varilla infinitamente larga, podemos
utilizar:

Caso B  la pérdida de calor en el extremo es insignificante

𝑞𝑓 = √ℎ𝑃𝑘𝐴𝐶 𝜃𝑏 tanh 𝑚𝐿
Caso B = Caso D

qf = Mtanh mL = M
𝑞𝑓 = √ℎ𝑃𝑘𝐴𝐶 𝜃𝑏
qf = qf

tanh mL = 1
𝑞𝑓 = 𝑞𝑓 = tanh 𝑚𝑙 = 1

arctanh(1) = 𝑚𝐿

mayor que -1 (-1 y 1) y menor que 1

Si:

tanh 𝑚𝑙 ≥ 0,99
70

𝑚𝐿 ≥ 2,65

1
2,65 𝑘𝐴𝐶 2
𝐿 ≥ 𝐿∞ ≡ = 2,65 [ ]
𝑚 ℎ𝑃

𝐿∞𝐶𝑢 = 0,19𝑚

𝐿∞𝐴𝑙 = 0,13𝑚

𝐿∞𝐴𝑐 = 0,04𝑚

2,3
Por teoría no tiene sentido extender la aleta mas allá de 𝐿 =
𝑚

2,3 2,65
≤𝐿≤
𝑚 𝑚

Desempeño de una aleta

Las aletas se utilizan para aumentar la transferencia de calor de una fuente porque aumentan el
área efectiva de la superficie. Sin embargo la aleta misma produce una resistencia a la
transferencia de calor.

 No ∋ seguridad que la transferencia aumente.

Evaluamos la efectividad de la aleta 𝜀𝑓

𝜀𝑓 se define como la razón de la T.C de la aleta a la TC que existirán sin la aleta.

𝑞𝑓
𝜀𝑓 =
ℎ𝐴𝑐,𝑏 𝜃𝑏
71

donde, 𝐴𝑐,𝑏 área de la sección transversal en la base de la aleta .

Para cualquier diseño 𝜀𝑓 debe ser tan grande como sea posible.

Para una aleta infinitamente larga

1
𝑘𝑃 2
𝜀𝑓 = ( )
ℎ𝐴𝐶

 La efectividad de una aleta aumenta por la elección de un material de alta conductividad

térmica.
Aleaciones de Al y Cu son las adecuadas pero el Al tiene bajo costo y peso.

 La efectividad aumenta al aumentar la razón del perímetro al área de la sección

transversal.
Por esto se prefiere el uso de aletas delgadas pero poco espaciadas, con la salvedad de
que el hueco de la aleta no se reduzca de tal manera que se reduzca el coeficiente de
convección.

1
𝑘𝑃 2
𝜀𝑓 = ( )
ℎ𝐴𝐶

Aquí el caso de aleta se justifica si h es pequeño, h es pequeño para gases  la necesidad

de aletas en este tipo de fluido es grande.

Si tenemos gas y liquido las aletas van del lado del gas.

El radiador de un carro las aletas van en la superficie externa del tubo.

Si 𝜀𝑓 > 2 𝑒𝑠 𝑐𝑟𝑖𝑡𝑒𝑟𝑖𝑜 𝑝𝑎𝑎 𝑢𝑠𝑎𝑟 𝑎𝑙𝑒𝑡𝑎𝑠

1
𝑘𝑃 2
Para 𝜀𝑓 (ℎ𝐴 ) >4
𝐶

Proceso convección libre h (W/m2°𝑲)

Gases 2 – 25
Líquidos 50 – 1000
C.F gases 25 – 250
C.F líquidos 50 – 20000
72

Para el caso B 98% se alcanza con mL=2.3b por lo que no tiene sentido extender las aletas más
allá de L=2.3/m
2.3 2.65
≤𝐿≤
𝑚 𝑚
También en función de resistencia térmica se puede cuantificar

Resistencia de aleta
𝜃𝑏
𝑅𝑡,𝑓 = Diferencia de temperaturas entre la base y el fluido
𝑞𝑓

Este resultado es útil cuando tenemos superficies con aletas mediante un circuito térmico

La resistencia en la base expuesta debido a la convección

1
𝑅𝑡,𝑏 =
ℎ𝐴𝑐,𝑏
𝑅𝑡,𝑏
𝐸𝑓 = La efectividad en una razón de resistencias
𝑅𝑡,𝑓

En función de las temperaturas

El impulso máximo para la convección es la diferencia de temperaturas en la base (x=0) y el

fluido Ɵb=𝑇𝑏 − 𝑇∞ . Si toda la aleta estará a la temperatura de la base tendríamos la máxima
rapidez para disipar energía.

Pero siempre existe un gradiente de temperatura a lo largo de la aleta

Af: es el área de la superficie de la aleta

Cuando la perdida en el extremo es insignificante B

𝑛𝑓 : Se aproxima a máx. 1 y min 0

Si L→0 o L→∞

Lc=L+ (t/2) para una aleta rectangular

Lc=L+ (D/4) para una aleta recta alfiler.

73

Los errores asociados a la aproximacion son insignificantes

En una aleta rectangular si 𝑤 ≫ 𝑡

𝑤:ancho

𝑡:espesor

Se aproxima el perímetro a: P=2w

Eficiencia global:

𝒒𝒕 es la transferencia de calor total del área de

la superficie 𝑨𝒕 asociada con las aletas y la parte
expuesta de la base (denominada superficie primaria).
Si existe N aletas cada una de las áreas superficiales 𝑨𝒇
y el área de la superficie primaria 𝑨𝒃

𝑨𝒕 = 𝑵𝑨𝒇 + 𝑨𝒃

La transferencia total por convección de las aletas y la superficie principal (sin aleta) se expresa
como:

h: equivalente para las superficies principal y con aletas y 𝒏𝒇

74

en la equivalencia de una sola aleta.

Por resistencia térmica

𝑞𝑡 = ℎ𝐴𝑡 𝑛𝑜 𝜃𝑏

Aletas como parte integral

Aletas adheridas a la base

75

𝑅𝑡,𝑐 ≪ 𝑅𝑡,𝑓

Para ingresar en el grafico 3.18 y 3.19

ℎ 1/2 𝑟2𝑐
𝐿𝑐 3/2 ( )
𝑘𝐴𝑝 𝑟1

Calculo de 𝒏𝒇 (método gráfico)

1)

Ap determinantes

Lc

2) calculamos

ℎ 1/2 𝑟
𝐿𝑐 3/2 (𝑘𝐴𝑝) 2𝑐
𝑟1

3) determinamos 𝒏𝒕 de las gráficas 3.18 y 3.19 (pág. 123)

4) calculamos calor ideal

a) aleta rectangular recta

b) aleta triangular recta

c) aleta circunferenciales

a) 𝑞 = ℎ(𝑃 ∗ 𝐿𝑐)(𝑇𝑏 − 𝑇∞ )

b) 𝑞 = ℎ(𝑃 ∗ 𝐿𝑐)(𝑇𝑏 − 𝑇∞ )

c) 𝑞 = ℎ 2𝜋(𝑟2𝑐 2 − 𝑟1 2 )(𝑇𝑏 − 𝑇∞ )

5) 𝑞𝑟 = 𝑛𝑓 ∗ 𝑞
76

Ejemplo:

Se colocan aletas circunferenciales de acero de 3/8 “de largo y t=1/4” en un tubo de 1” de

diámetro, con el fin de disipar calor desde la superficie del tubo que se encuentra a 400°F al
ambiente que se encuentra a 70°𝐹. Calcular el calor que disipa cada aleta considerando un
ℎ𝑒𝑥𝑡𝑒𝑟𝑛𝑜 = 20 𝐵𝑇𝑈/𝑝𝑖𝑒 2 °𝐹 y 𝑘𝑎𝑐𝑒𝑟𝑜 = 30 𝐵𝑇𝑈/𝑝𝑖𝑒°𝐹

L=3/8 “→1/32’

r1=1/2”→1/24’

t=1/4”→1/48’

Tb=400°𝐹 h=20 𝐵𝑇𝑈/𝑝𝑖𝑒 2 °𝐹

T∞=70°𝐹 k=30 𝐵𝑇𝑈/𝑝𝑖𝑒°𝐹

Lc=L+t/2

Lc=1/32 + (1/48)/2

Lc=1/24’

r2c=r1+Lc

r2c=1/24 + 1/24

r2c=1/12’

Ap=Lc*t

Ap=1/24’+1/48’

Ap=8.6805∗ 10−4
3⁄
3/2
ℎ 1/2 1 2 20 1/2
𝐿𝑐 ( ) =( ) ( ) = 0.2357
𝑘𝐴𝑝 24 30 ∗ 8.6805 ∗ 10−4
𝑟2𝑐 1 ∗ 24
= =2
𝑟1 12 ∗ 1
Grafico 3.19

𝑛𝑓 = 93%

𝑞𝑓 = ℎ 2𝜋(𝑟2𝑐 2 − 𝑟1 2 )(𝑇𝑏 − 𝑇∞ )

𝐵𝑇𝑈 12 12
𝑞𝑓 = 20 2𝜋([ ] − )𝑝𝑖𝑒 2 (400 − 70)°𝐹
𝑝𝑖𝑒 2 °𝐹 12 24
𝑞𝑓 = 215.98 𝐵𝑇𝑈

𝑞𝑔 = 𝑛𝑓 ∗ 𝑞𝑓 = (0.93) + 215.98 = 200.86 𝐵𝑇𝑈

77

Eficiencia global

Aletas son parte integral de la base

Aletas adheridas a la base

Ejemplo 3.9

El cilindro del motor de una motocicleta está fabricado de aleación de aluminio 2024T6 y tiene
una altura H=0.15m y un diámetro exterior D=50mm. Bajo condiciones de operación típicas la
superficie exterior del cilindro está a una temperatura de 500°𝐾 y se expone al aire ambiental a
300°𝐾, con un coeficiente de convección de 50𝑊/𝑚2 °𝐾.

Unas aletas anulares están fundadas integralmente con el cilindro para aumentar la
transferencia de calor a los alrededores. Considere cinco de estas aletas, de espesor t=6mm,
longitud L=20mm e igualmente espaciadas.

¿Cuál es el aumento en la transferencia de calor debido al uso de aletas?

78

Con las aletas colocadas la transferencia de calor es

r2c=r1+Lc

Lc=L+t/2

Lc=0.02+0.003

Lc=0.023

r2c=0.025+0.023

r2c=0.048

𝐴𝑡 = 𝑁𝐴𝑓 + 2𝜋𝑟1(𝐻 − 𝑁𝑡)

𝐴𝑡 = 5 ∗ 0.0105𝑚2 + 2𝜋(0.025)(0.15 − 5 ∗ 0.06)𝑚2

𝐴𝑡 = 0.0716𝑚2
𝑟2𝑐 0.048
= = 1.92
𝑟1 0.025
Ap=Lc*t

Ap=0.023*0.006

Ap=1.38∗ 10−4 𝑚2

ℎ 1/2 3 50 1/2
𝐿𝑐 3/2 ( ) = (0.023) ⁄2 ( ) = 0.1539
𝑘𝐴𝑝 186 ∗ 1.38 ∗ 10−4
Grafico 3.19

𝑛𝑓 = 95%
79

Sin aletas la transferencia de calor por convección seria:

q=236 W

Ejemplo 3.10

La transferencia de un transistor se puede aumentar insertándolo en una base de aluminio 𝑘 =

200 𝑊⁄𝑚𝐾 que tiene 12 aletas longitudinalmente fabricadas integralmente sobre su superficie
externa. El radio del transistor y la altura son r1=2mm y H=6mm, respectivamente, mientras que
las aletas n son de longitud L=r3-r2=1mm y espesor uniforme t=0.7mm. El espesor de la base de
la manga es r2-r1=1mm, y la resistencia de contacto de la interfaz base-transistor es 𝑅"𝑡,𝑐 =
2
10−3 𝑚 𝐾⁄𝑊. Aire a 𝑇∞ = 20𝐶 fluye sobre la superficie de la aleta, lo que proporciona un
coeficiente de convección aproximadamente uniforme de h=25𝑊⁄ 2 .
𝑚 𝐾
1. Suponiendo una transferencia de calor unidimensional en la dirección radia, dibuje el
circuito equivalente para la transferencia de calor de la capa del transistor (r=r1) al aire.
2. Evalué cada una de las resistencias en el circuito anterior. Si la temperatura del
transistor es T1=80C ¿Cuál es la rapidez de la transferencia de calor de la base?

1) Condición de estado estable El circuito debe explicar la

2) La transferencia de calor es insignificante en las resistencia de contacto,
superficies inferior y superior. conducción en la base,
3) Conducción radial unidimensional convección en la base expuesta y
conducción/convección en las
aletas
80

4) Radiación insignificante.

Para una sola aleta

𝜃𝑏 𝜃𝑏
𝑞𝑓 = 𝑅𝑡,𝑓 =
𝑅𝑡,𝑓 𝑞𝑓
81

Para una sola aleta:

ℎ
𝜃𝑏 cosh 𝑚𝑙 + ( ) 𝑠𝑒𝑛ℎ 𝑚𝑙
𝑅𝑡,𝑓 = = 𝑚𝑘
𝑞𝑓 [𝑠𝑒𝑛ℎ 𝑚𝑙 + ( ℎ ) cosh 𝑚𝑙]()1/2
𝑚𝑘
cosh(0.1997) + (6.25 ∗ 10−3 )𝑠𝑒𝑛ℎ (0.1997)
𝑅𝑡,𝑓 =
[𝑠𝑠𝑒𝑛ℎ (0.1997) + 6.25 ∗ 10−3 cosh(0.1997)](0.0168)
𝑅𝑡,𝑓 = 293.13𝐾/𝑊

Para 12 aletas:
293.13
R t,f(12) = K/W
12
𝑅𝑡,𝑓(12) = 24.427𝐾/𝑊

Para la base expuesta:

Ejercicio 3.113

Un arreglo experimental para medir la conductividad térmica de materiales solidos implica el

uso de las varillas largas que son equivalentes en todos los aspectos excepto que una está
fabricada de un material estándar de kA conocida mientras que el otro está fabricado con el
material cuya kB queremos conocer. Ambas varillas se unen en un extremo a una fuente de calor
de temperatura fija Tb, se expone a un fluido de temperatura T∞ y se instrumentan con
termopares para medir la temperatura a una distancia fija x1 de la fuente de calor. Si el material
estándar es aluminio en kA=200(W/mK) y las medidas revelan valores de TA=100C y Tb=60C a
x1 para Tb=100C y T∞=25C ¿Cuál es el calor de kB?
82

Varillas muy largas

𝜃
= 𝑒 −𝑚𝑥
𝜃𝑏
𝑇 − 𝑇∞
= 𝑒 −𝑚𝑥
𝑇𝑏 − 𝑇∞
𝑇−𝑇∞ ℎ𝑃 1/2
ln (𝑇 ) = 𝑚𝑥 = [𝑘𝐴] 𝑥
𝑏 −𝑇∞

𝑇𝐴 − 𝑇∞ ℎ𝑃 1/2
ln ( ) = 𝑚𝑥 = [ ] 𝑥 (1)
𝑇𝑏 − 𝑇∞ 𝑘𝐴 𝐴

𝑇𝐵 − 𝑇∞ ℎ𝑃 1/2
ln ( ) = 𝑚𝑥 = [ ] 𝑥 (2)
𝑇𝑏 − 𝑇∞ 𝑘𝐵 𝐴
(1)/ (2)
𝑇 −𝑇
ln (𝑇𝐴 − 𝑇∞ ) 𝑘𝐵 1/2
𝑏 ∞
=[ ]
𝑇 −𝑇 𝑘𝐴
ln ( 𝑇𝐵 − 𝑇∞ )
𝑏 ∞

𝑘𝐵 = 56.6 𝑊⁄𝑚𝐾
83

Ejercicio 3.124

Un calentado de aire consiste en un tubo de acero (k=20W/mK) con radios interno y externo de
r1=13 mm y r2=16 mm, respectivamente y ocho aletas longitudinales fabricadas integralmente,
cada una de espesor t=3mm. Las aletas se extienden a un tubo concéntrico que tiene radio
r3=40mm y aislado en la superficie externa. Agua a T∞= 90 C fluye a través del tubo interno
mientas que aire a 𝑇∞,0 = 25 𝐶 fluye a través de la región anular formada por el tubo
concéntrico más grande

a) Dibuje el circuito térmico equivalente del calentador y relacione c/resistencia térmica

con los parámetros apropiados del sistema
b) Si hi=5000𝑊⁄ 2 y ho=200𝑊⁄ 2 ¿Cuál es la transferencia de calor por unidad de
𝑚 𝐾 𝑚 𝐾
longitud?

a)

𝑛𝑓 = 0.49
84

𝑁 𝐴𝑓
𝑛𝑜 = 1 − (1 − 𝑛𝑓 )
𝐴𝑡
8 ∗ 0.048
𝑛𝑜 = 1 − (1 − 0.49)
0.4602
𝑛𝑜 = 0.5744
1 1
𝑅𝑡,𝑜 = =
𝑛𝑜 ℎ𝐴𝑡 0.5744 ∗ 200 ∗ 0.4602
𝑚𝐾
𝑅𝑡,𝑜 = 18.91 ∗ 10−3
𝑊

𝑇∞,𝑖 − 𝑇∞,𝑜 30 − 25 𝑊
𝑞= = = 2824
𝑅𝑡𝑜𝑡𝑎𝑙 ∑𝑅 𝑚
Ejercicios: 3.126/3.127/3.128/3.129

3.142/3.143/3.144/3.145

7.2 considere aire atmosférico a 25 C en un flujo paralelo a 5m/s sobre ambas superficies de
una placa de 1m de longitud que se mantiene a 75C.

a) determine el espeso de la capa límite de velocidad

b) el coeficiente local de convección hL

c) el flujo de calor al final de la placa

7.3 sobre ambas superficies de una placa plana de 1m de longitud se mantiene a 20C, fluye
aceite de motor a 100C y a una velocidad de 0.1 m/s

a) Espesor de la capa límite de velocidad y capa límite térmica al final de la placa

b) El flujo local de calor hL y q”x

85

Ejercicio 3.130

Se instala aletas anulares de aluminio de 2mm de espesor y 15mm de longitud sobre un tubo de
aluminio de 30mm de diámetro. Se sabe que la resistencia de contacto térmico entre una aleta
2
y el tubo es 𝑅"𝑡,𝑐 = 2 ∗ 10−4 𝑚 𝐾⁄𝑊. Si la pared del tubo está a 100C y el fluido contiguo está
a 25C con h=75𝑊⁄ 2 ¿Cuál es la transferencia de calor de 1 sola aleta?
𝑚 𝐾
¿Cuál sería la transferencia de calor si la resistencia de contacto pudiera eliminarse?

𝑇𝑡 − 𝑇∞
𝑞𝑓 =
𝑅𝑡,𝑐 + 𝑅𝑓
𝑅"𝑡,𝑐
𝑅𝑡,𝑐 =
𝐴𝑏
Ab=2π (0.015*0.002)

Ab=1.88495∗ 10−4 𝑚2
2
2 ∗ 10−4 𝑚 𝐾⁄𝑊
𝑅𝑡,𝑐 =
1.88495 ∗ 10−4 𝑚2

𝑅𝑡,𝑐 = 1.06 𝐾⁄𝑊

𝑞𝑓 𝑞𝑓
𝑛𝑓 = =
𝑞𝑚𝑎𝑥 ℎ 𝐴𝑓 𝜃𝑏
𝜃𝑏
𝑅𝑡,𝑓 = 𝑞𝑓 = 𝑛𝑓 ∗ 𝑞𝑚𝑎𝑥 = 𝑛𝑓 ∗ ℎ 𝐴𝑓 𝜃𝑏
𝑞𝑓
1
𝑅𝑡,𝑓 =
𝑛𝑓 ∗ ℎ 𝐴𝑓

𝐴𝑓 = 2𝜋(𝑟2,𝑐 2 − 𝑟1 2 )
86

𝑟2,𝑐 = 𝑟2 + 𝑡/2
2𝑚𝑚
𝑟2,𝑐 = 30𝑚𝑚 +
2
𝑟2,𝑐 = 0.031𝑚
𝑡
𝐿𝑐 = 𝐿 +
2
2𝑚𝑚
𝐿𝑐 = 15𝑚𝑚 +
2
𝐿𝑐 = 0.016𝑚
𝐴𝑝 = 𝐿𝑐 ∗ 𝑡
𝐴𝑝 = 0.016 ∗ 0.002

𝐴𝑝 = 3.2 ∗ 10−5 𝑚2

ℎ 1/2
𝐿𝑐 3/2 ( ) = 0.20
𝑘𝐴𝑝
𝑛𝑓 = 0.94
1
𝑅𝑓 =
0.94 ∗ 75 ∗ 2𝜋(0.0312 − 0.0152 )
𝐾
𝑅𝑓 = 3.07
𝑊
𝑇𝑡 − 𝑇∞
𝑞𝑓 =
∑𝑅
100 − 25
𝑞𝑓 = = 18.2 𝑊
(1.06 + 3.07)

Sin resistencia de contacto

𝜃𝑏 100 − 25
𝑞𝑓 = =
𝑅𝑓 3.07

𝑞𝑓 = 24.4 𝑊
87

CONDUCCION TRANSITORIA

kg 𝐽
ρ ∗ Cp = capacidad termica volumetrica ρ [𝑚 3 ] 𝐶𝑝[𝑘𝑔°𝐾]

CTV = ρ ∗ Cp [𝐽⁄ 3 ]
𝑚 °𝐾
α= difusividad térmica
𝑘 2
α = 𝐶𝑝 [𝑚 ⁄𝑠] Mide la capacidad de un material para conducir energía térmica en relación con
su capacidad de almacenar energía térmica.

Cuando un sólido experimenta un cambio súbito en su ambiente térmico, se considera una pieza
forjada de metal caliente que se encuentra a una
temperatura inicial uniforme Ti y que se templa
por inmersión en un líquido que se encuentra a
una temperatura más baja T∞<Ti, si el templado
comienza en t=0, la temperatura del solido
disminuye para un t > 0, hasta llegar a T∞. Esto
se da debido a que la transferencia de calor por
convección en la interfaz solido-liquido.

Ti: temperatura uniforme del cuerpo (inicial)

Suponemos que la temperatura del solido es espacialmente uniforme en cualquier instante

Los gradientes de temperatura del solido son insignificantes. Método de la resistencia interna
despreciable.

Si no tenemos gradientes de temperatura no podemos calcular la difusión de calor.

Para resolver el problema usamos la temperatura transitoria, para esto realizamos un balance
de energía global en el sólido.

Debemos relacionar la velocidad de perdida de calor en la superficie con la rapidez de cambio

de la energía interna.

1)
̇
−𝐸𝑠𝑎𝑙𝑒 = 𝐸̇𝑎𝑙𝑚
2)
𝑑𝑇
−ℎ𝐴𝑠(𝑇 − 𝑇∞ ) = 𝜌𝑉𝐶𝑝
𝑑𝑡
3)

𝜃 ≡ 𝑇 − 𝑇∞ Separando variables e
𝜌𝑉𝐶𝑝 𝑑𝜃 integrando desde la
= −𝜃 condición inicial para la
ℎ𝐴𝑠 𝑑𝑡
que t=0 y T(0)=Ti
88

𝜃
𝑡
𝜌𝑉𝐶𝑝 𝑑𝜃
∫ = − ∫ 𝑑𝑡
ℎ𝐴𝑠 𝑑𝑡 𝑜
𝜃𝑖

4)
𝜃𝑖 = 𝑇𝑖 − 𝑇∞
5)
𝜌𝑉𝐶𝑝 𝜃𝑖
ln =𝑡
ℎ𝐴𝑠 𝜃
6)
𝜃𝑖 𝑇 − 𝑇∞ ℎ𝐴𝑠
= = 𝑒𝑥𝑝 [− ( ) 𝑡]
𝜃 𝑇𝑖 − 𝑇∞ 𝜌𝑉𝐶𝑝

Constante térmica de tiempo

𝜌𝑉𝐶𝑝
τt =
ℎ𝐴𝑠
7)
1
τt = [ ] (𝜌𝑉𝐶𝑝) = 𝑅𝑡 𝐶𝑡
ℎ𝐴𝑠
𝑅𝑡 : Resistencia a la transferencia de calor por convección

𝐶𝑡 : Resistencia interna despreciable del solido

NOTA: cualquier aumento en 𝑅𝑡 y 𝐶𝑡 ocasionaría que un sólido responda más lentamente a

cambios en su ambiente térmico.

Transferencia total de calor Q

𝑡 𝑡

𝑄 = ∫ 𝑞 𝑑𝑡 = ℎ𝐴𝑠 ∫ 𝜃𝑑𝑡
0 0

Sustituyendo θ de la ecuación 6) e integrar se tiene

89

𝑡
𝑄 = (𝜌𝑉𝐶𝑝)𝜃𝑖 [1 − 𝑒𝑥𝑝 (− )]
𝜏𝑡
Además Q está relacionada con el cambio de energía interna de salida

−𝑄
≡ ∆𝐸𝑎𝑙𝑚 𝑝𝑎𝑟𝑎 𝑒𝑙 𝑡𝑒𝑚𝑝𝑙𝑎𝑑𝑜 𝑄 𝑒𝑠 𝑝𝑜𝑠𝑖𝑡𝑖𝑣𝑎, 𝑐𝑢𝑎𝑛𝑑𝑜 𝑒𝑙 𝑠𝑜𝑙𝑖𝑑𝑜 𝑠𝑒 𝑐𝑎𝑙𝑖𝑒𝑛𝑡𝑎 𝑄 𝑒𝑠 𝑛𝑒𝑔𝑎𝑡𝑖𝑣𝑎

Consideración de Conducción en estado estable a través de una pared plana de área A.

Esta superficie se mantiene a 𝑇𝑠,1 y la otra se

expone a un fluido de temperatura 𝑇∞ < 𝑇𝑠,1 por
lo que 𝑇𝑠,2 será un valor intermedio 𝑇∞ < 𝑇𝑠,2 <
𝑇𝑠,1

Curva a): Bi << 1

Curva b): Bi ≈ 1

Curva c): Bi >> 1

𝑇∞ < 𝑇𝑠,1

𝑇∞ < 𝑇𝑠,2 < 𝑇𝑠,1

𝑘𝐴
(𝑇 − 𝑇𝑠,2 ) = ℎ𝐴(𝑇𝑠,1 − 𝑇∞ )
𝐿 𝑠,1
𝐿
𝑇𝑠,1 − 𝑇𝑠,2 ( ⁄𝑘𝐴) 𝑅𝑐𝑜𝑛𝑑 ℎ𝐿
= = =
𝑇𝑠,1 − 𝑇∞ (1⁄ ) 𝑅𝑐𝑜𝑛𝑣 𝑘
ℎ𝐴

La resistencia a la conducción dentro del solido es mucho menor que la resistencia a la

convección a través de la capa límite del fluido.

9)
ℎ𝐿
𝐵𝑖 ≡
𝑘
Biot  proporciona la variación de temperatura en el sólido / variación de temperatura entre la
superficie y el fluido.
90

La T.C se la ∆𝑇 Muy
𝑇(𝑥,𝑡) ≈ 𝑇(𝑡) relaciona entre grande la T.C
Variación de el sólido y el ∆𝑇 en el interior
temp y el fluido, Significante del solido es
tiempo internamente muy grande
la ∆𝑇
significante

es
Se debe satisfacer la siguiente condición:
insignificante.
ℎ𝐿
𝐵𝑖 ≡ < 0.1
𝑘
Longitud característica Lc

Lc -> relación entre el volumen del sólido y el área de superficie

𝐿𝑐 ≡ 𝑉/𝐴𝑠
Para una pared plana de 2L  Lc=L

Cilindro largo: 𝐿𝑐 = 𝑟𝑜 /2

Esfera: 𝐿𝑐 = 𝑟𝑜 /3

Para un cilindro o esfera calentada en forma simétrica Lc sería igual al radio real 𝑟𝑜 en lugar de
𝑟𝑜 /2 o 𝑟𝑜 /3.

Con Lc=V/As
ℎ 𝐴𝑠 𝑡 ℎ𝑡 ℎ 𝐿𝑐 𝑘 𝑡
= = = ∗ 2
𝑉 𝜌 𝐶𝑝 𝐿𝑐 𝑘 𝜌𝑐 𝐿𝑐
ℎ𝐿𝑐 𝛼
= 𝑘
∗ 𝐿𝑐 2 Fo es un tiempo sin dimensión que junto con Biot son
ℎ 𝐴𝑠 𝑡 𝛼𝑡 características de conducción transitoria
= 𝐵𝑖 ∗ 𝐹𝑜 𝑑𝑜𝑛𝑑𝑒 𝐹𝑜 ≡ 2
𝜌 𝑉 𝐶𝑝 𝐿𝑐
91

𝜃𝑖 𝑇 − 𝑇∞
= = exp(−𝐵𝑖 ∗ 𝐹𝑜)
𝜃 𝑇𝑖 − 𝑇∞
Ejemplo 5.1

Una unión termopar cuya forma se aproxima a una esfera, se usara para la medición de la
temperatura en un flujo de gas. Se sabe que el coeficiente de convección entre la superficie de
unión y el gas es h= 40𝑊⁄ 2 , y que las propiedades termofisicas de la unión son
𝑚 °𝐾
𝑘 = 20𝑊⁄ 𝐽 3
𝑚°𝐾 c=400 ⁄𝑘𝑔°𝐾 y 𝜌 = 8500 𝑘𝑔/𝑚 . Determine el diámetro de la unión esta a
25°C y se coloca en un flujo de gas que está a 200°C ¿Cuánto tiempo tardara la unión en alcanzar
199°C?

Ejemplos 5.2/5.4
92

Ejercicios 5.6/5.7/5.9/5.31

Análisis General del Método de Resistencia interna Despreciable

Resistencia Interna

Aunque la conducción transitoria en un sólido normalmente se inicia mediante la transferencia

de calor por convección hacia o desde un fluido contiguo, otros procesos tal vez introduzcan
condiciones térmicas transitorias dentro del sólido, por ejemplo:

Un sólido se separa de sus alrededores mediante un gas o un vacío. Si las temperaturas del sólido
y los alrededores difieren el intercambio de radiación ocasiona que cambie la energía térmica
interna y por ello la temperatura del sólido.

Aplicando conservación de
energía para cualquier instante t:

Ecuacion diferencial ordinaria de

primer orden no lineal y no
homogenea no es posible
integrara para obtener una soluccion exacta.

Pero si no existe flujo de calor impuesto o generado y convección son insignificantes frente a la
radiación, se reduce a:

Al separar variables e integrar desde la condición inicial hasta cualquier tiempo t se tiene:

Al evaluar las integrales y reacomodar con el tiempo que se requiere para alcanzar T
93

Pared plana con convección

Solución EXACTA

𝑥 ∗ Es una coordenada adimensional espacial

𝑥
𝑥∗ =
𝐿
𝑡 ∗ Tiempo adimensional
𝛼𝑡
𝑡∗ = = 𝐹𝑜
𝐿2
L: es la mitad del espesor de la pared.

Si el espesor es pequeño en relación al ancho y la altura de la pared es razonable suponer que

la conducción ocurre en la dirección X

Cilindro Infinito o Esfera

Si la pared esta inicialmente a una temperatura uniforme

𝑇(𝑥,0) = 𝑇𝑖 y se la sumerge súbitamente en un fluido 𝑇∞ ≠ 𝑇𝑖

las temperaturas se encuentran resolviendo la ecuación:

𝑑2 𝜃 ∗ 𝑑𝜃 ∗ 𝜽 Relación adimensional de la variable

= 𝜽∗ ≡
𝑑𝑥 ∗2 𝑑𝐹𝑜 𝜽𝒊 𝑇−𝑇∞
dependiente𝜃 ∗ =
𝑇𝑖−𝑇∞

Para las condiciones:

𝜃 ∗ (𝑥 ∗ , 0) = 1
𝑑𝜃 ∗
=0
𝑑𝑥 ∗ 2 𝑥=0
𝑑𝜃 ∗
= −𝐵𝑖𝜃 ∗ (1, 𝑡 ∗ )
𝑑𝑥 ∗ 𝑥∗ =1
La convección para las superficies en 𝑥 ∗ = ±1 son las mismas, la distribución de temperaturas
en cualquier instante debe ser simétrica alrededor del plano medio (𝑥 ∗ = 0)
94

La solución exacta:

1)
∞

𝜃 = ∑ 𝐶𝑛 exp(−𝜁𝑛2 𝐹𝑜)(𝑐𝑜𝑠𝜁𝑛 𝑥 ∗ )
∗

𝑛=1

Cn ^ 𝜁 se encuentra en la tabla 5.1 pág. 227

2)
𝛼𝑡
𝐹𝑜 =
𝐿2
3)
4𝑠𝑒𝑛 𝜁𝑛
𝐶𝑛 =
2𝜁 + 𝑠𝑒𝑛(2𝜁𝑛)
Valores característicos de 𝜁𝑛 son las raíces positivas de la ecuación

4) 𝜁𝑛 tan 𝜁𝑛 = 𝐵𝑖

Solución aproximada:

Para valores de Fo > 0.2

5)
𝜃 ∗ = 𝐶1 exp(−𝜁12 𝐹𝑜)(𝑐𝑜𝑠𝜁1 𝑥 ∗ )
6) La dependencia de la temperatura con
∗ ∗
𝜃 = 𝜃0 cos( 𝜁1 𝑥 ) ∗ respecto al tiempo x en cualquier lugar
dentro de la pared es la misma que la de la
temperatura del plano medio
7)
𝑇0 − 𝑇∞ Representa la temperatura en el plano
∗
𝜃 = medio 𝑥 ∗ = 0
𝑇1 − 𝑇∞
8)

𝜃 ∗ = 𝐶1 exp(−𝜁12 𝐹𝑜)

TRANSFERENCIA TOTAL DE ENERGIA

En algunas situaciones es útil conocer la energía total que

disminuye en la pared en cualquier tiempo t eh el proceso
transitorio.

Se aplica para condiciones de tiempo limitado por la

condición inicial (t=0) y cualquier t > 0
9)
95

Energía interna inicial de la pared. Cantidad máxima de transferencia de energía que podría
ocurrir para un t=∞

10) 𝜃0 ∗ 𝑑𝑒 𝑙𝑎 𝑒𝑐𝑢𝑎𝑐𝑖𝑜𝑛 8)

C1 ^ 𝜁1 tabla 5.1

Cilindro infinito

𝜃 ∗: Energía total transferida de la pared

𝐽1 , 𝐽0 : Funciones de Bessel de primera clase Tabla B-4

Esfera

Donde:
𝛼𝑡
𝐹𝑜 =
𝑟𝑜 2

Donde 𝜁𝑛 son las raíces positivas de la ecuación trascendental 1 − 𝜁𝑛 cot 𝜁𝑛 = 𝐵𝑖

Ejemplo 5.2/5.4 Propuestos 5.5/5.6/5.7/5.31

Ejemplo 5.3 pág. 231

Considere una tubería de acero (AISI 1010) que tiene 1m de diámetro interno y una pared con
espesor de 40 mm. La tubería está fuertemente aislada en el interior y antes del inicio del flujo
las paredes de la tubería se encuentran a una temperatura uniforme de -20 C. con el inicio del
flujo se bombea aceite caliente a 60 C por la tubería, con lo que se crea una condición conectiva
de superficie que corresponde a ℎ = 500 𝑊⁄ 2 en la superficie interior de la tubería.
𝑚 𝐾
1) Cuáles son los números de Biot y Fourier apropiados, 8 minutos después de iniciado el
flujo
2) A t=8min ¿Cuál es la temperatura de la superficie exterior cubierta por aislante?
3) ¿Cuál es el flujo de calor q” 𝑊⁄ 2 a la tubería desde el aceite en t= 8 min?
𝑚
4) ¿Cuánta energía por metro de longitud de tubería se ha transferido del aceite en t= 8
min?
96

𝐵𝑖 = 0.313 > 0.1 No se puede usar resistencia interna

𝐹𝑜 = 5.64 > 0.2 Podemos usar convección en una pared plana

Temperatura

Del plano medio

TABLA 5.1

𝐵𝑖 = 0.313; 𝐶1 = 1.047; 𝜁 = 0.531 𝑟𝑎𝑑

INTERPOLAR

Bi 𝜁𝑖 C1
0.30 0.5218 1.0450
0.313 0.531 1.047
0.40 0.5932 1.0580

Para 8 min tenemos 𝑇(0,8min)

𝑇0 − 𝑇∞
𝜃∗ =
𝑇𝑖 − 𝑇∞
97

𝑇0 = 𝑇∞ + 𝜃 ∗ (𝑇𝑖 − 𝑇∞ )

𝑇0 = 60 + 0.213(−20 + 60)
𝑇0 = 42.92 𝐶

3) la transferencia de calor a la superficie interna en x=L se da por convección y para cualquier

tiempo

Ley de Newton

t=8min=480 s

El flujo de calor en t=8 min

Transferencia total de energía

98

INTRODUCCION A LA CONVECCION

Hemos analizado la transferencia de calor por conducción y se ha analizado la convección solo

hasta el punto en que proporciona una posible condición de frontera en problemas de
conducción.

En esta sección analizaremos la convección como la transferencia de energía entre una

superficie y un fluido que se mueve sobre esta  comprender los mecanismos físicos que
fundamentan la transferencia de calor por convección.

 Determinar los métodos adecuados para realizar los cálculos de transferencia de calor
por convección

El flujo de calor local se expresa como: 𝑞" = (𝑇𝑠 − 𝑇∞ )ℎ (1)

como las condiciones de flujo varían de punto a punto sobre la
superficie q” y h también varían a lo largo de la superficie.

 𝑞 = ∫𝐴𝑠 𝑞"𝑑 𝐴𝑠 (2)

 𝑞 = (𝑇𝑠 − 𝑇∞ ) ∫𝐴𝑠 ℎ 𝑑 𝐴𝑠 (3)

Definiendo un coeficiente de convección promedio ℎ̅ para toda la superficie el calor

transferido:

𝑞 = ℎ̅𝐴𝑠(𝑇𝑠 − 𝑇∞ ) (4)
Igualando (3) y (4)

1
ℎ̅ = ∫ ℎ 𝑑 𝐴𝑠 (5)
𝐴𝑠 𝐴𝑠

El flujo sobre una placa plana de manera especial h varía con la distancia x desde la primera
orilla por lo que la ecuacion (5) se reduce a:

1 𝐿
ℎ̅ = ∫ ℎ 𝑑𝑥 (6)
𝐿 0

Tarea ejemplo 6.1

Capa límite de velocidad o hidrodinámica

99

Considere el flujo sobre una placa plana.

Cuando las partículas entran en contacto
con la superficie adquieren una velocidad
cero. Estas partículas actúan para retardar
el movimiento de la capa contigua al fluido
por lo que también retarda el movimiento
de las partículas en la capa. Esto se repite
una y otra vez hasta llegar a una distancia 𝑦 = 𝛿 de la superficie.

Esta desaceleración del movimiento del fluido se asocia con los esfuerzos cortantes 𝝉 que actúan
en planos paralelos a la velocidad del fluido.

La cantidad 𝜹  espesor de capa límite y se define como el valor de 𝒚 para el que 𝜇 = 0.99𝜇∞
(velocidad del fluido).

Este perfil de velocidad se refiere a la forma en la que 𝝁 varía con 𝒚 a través de la capa límite.

 existe dos regiones distintas -> capa fluida delgada =, gradiente de velocidad y
esfuerzos cortantes grandes
->una región fuera en la que 𝜇∞ y 𝝉 son insignificantes.

Capa limite térmica

Si las temperaturas del flujo sobre la superficie difieren se produce una capa limite térmica. Al
inicio de la placa el perfil de temperatura es uniforme con 𝑇(𝑦) = 𝑇∞ . Pero las partículas del
fluido hacen contacto con la placa alcanzando el equilibrio térmico a la temperatura de la
superficie de la placa. A su vez intercambian energía estas partículas, con las de la capa
adyacente del fluido y se producen en el fluido gradientes de temperatura.
𝑇𝑠 −𝑇
𝑑𝑡 = 𝑦 En donde se cumple = 0.99
𝑇𝑠 −𝑇∞

El flujo de calor se puede obtener en y=0 a través de Fourier

(6)

En la superficie no existe movimiento del fluido y la transferencia de calor ocurre solo por
conducción.

Combinando (1) y (6)

100

(𝑇 − 𝑇∞ ) = 𝑐𝑡𝑒

𝑑𝑡 → 𝑖𝑛𝑐𝑟𝑒𝑚𝑒𝑛𝑡𝑜 𝑒𝑛 𝑥
𝑑𝑇
𝑑𝑖𝑠𝑚𝑢𝑛𝑢𝑦𝑒 𝑎𝑙 𝑎𝑢𝑚𝑒𝑛𝑡𝑎𝑟 𝑥
𝑑𝑦 𝑦=0

𝑞"𝑠 𝑦 ℎ 𝑑𝑖𝑠𝑚𝑢𝑛𝑢𝑦𝑒𝑛 𝑎𝑙 𝑎𝑢𝑚𝑒𝑛𝑡𝑎𝑟 𝑥

Convección

Capa limite

Es necesario determinar si el flujo es laminar o turbulento.

La transferencia de calor por convección depende de estas condiciones

Región turbulenta

Capa de amortiguamiento

Subcapa laminar

Laminar: Re ≤ 2100

Transición: Re < 3 ∗ 106

Turbulento: Re ≥ 3 ∗ 106

 en la capa limite laminar el movimiento del fluido es altamente ordenado y es posible

identificar líneas de flujo a lo largo de las cuales se mueven las partículas. El movimiento
del fluido a lo largo de una línea de flujo se caracteriza por componentes de la velocidad
en X y Y.
El componente de la velocidad en la dirección normal contribuye significativamente a la
trasferencia de movimiento, energía o especias a través de la capa limite.
El movimiento del fluido en la capa limite turbulenta es altamente irregular y se
caracteriza por fluctuaciones de velocidad; estas aumentan la transferencia de
movimiento, energía y especies por lo que aumenta la fricción de la superficie así como
la transferencia por convección.
 Se caracteriza por gradientes de velocidad y esfuerzos cortantes.
 101

En la región completamente turbulenta las condiciones

se caracterizan por un movimiento tridimensional
aleatorio de porciones grandes de fluido => aumentos en
el espesor de la capa límite, esfuerzo cortante de la
pared y coeficiente de convección.

El espesor 𝛿 y el
coeficiente h
varían

La transición comienza en alguna posición Xc la que se determina mediante el agrupamiento de

variables llamado Número de Reynolds.
𝜌 𝜇∞ 𝑥 x: es la distancia característica
Reynolds es una relación de fuerzas 𝑅𝑒𝑥 =
𝜇
de inercia y fuerzas viscosas. 𝑅𝑒𝑥 : # de Reynolds crítico y para
𝜇∞ 𝐿 𝑉 ∗ 𝐿
Si Re ↑ => fuerzas de inercia son↑ 𝑅𝑒𝐿 = = el flujo sobre una placa plana
𝑟 𝑟 varia de 105 ª 3 ∗ 106
Si Re peq => fuerzas viscosas son↑ dependiendo de la aspereza de la
superficie y del nivel de
turbulencia del flujo libre

𝜌 𝜇∞ 𝑥𝑐
𝑅𝑒𝑥,𝑐 = = 5 ∗ 105
𝜇
Valor representativo

Las propiedades del fluido varían con la temperatura a través de la capa límite
𝑇𝑠 + 𝑇∞
𝑇𝑓 = 𝑡𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑑𝑒 𝑝𝑒𝑙𝑖𝑐𝑢𝑙𝑎
2
𝑇𝑠 : temperatura de la superficie

𝑇∞ : temperatura de flujo libre

Variables Independientes adimensionales

𝑥 𝑦
𝑥∗ ≡ 𝑦∗ ≡
𝐿 𝐿
L: longitud característica para la superficie de interés (Ej.: longitud placa plana)
𝑢 𝑢
𝑢∗ ≡ 𝑣∗ ≡
𝑉 𝑉
𝑉: Velocidad a contracorriente
102

Flujo externo

𝑁𝑢𝑥 = 𝑓4 (𝑥 ∗ , 𝑅𝑒𝑥 , Pr) 𝑓 𝑒𝑠 𝑓𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑓𝑟𝑖𝑐𝑐𝑖ó𝑛

̅̅̅̅̅̅𝑥 = 𝑓5 (𝑅𝑒𝑥 , Pr)
𝑁𝑢
𝐶𝑝 𝑈 𝑣
𝑃𝑟 = =
𝑅 𝛼
Flujo laminar

Espesor para capa limite

5 5𝑥
𝛿= =
√𝜇∞ /𝑣𝑥 √𝑅𝑒𝑥
𝛿
𝑃𝑟1/3 = 𝑒𝑠𝑝𝑒𝑠𝑜𝑟 𝑐𝑎𝑝𝑎 𝑙𝑖𝑚𝑖𝑡𝑒 𝑡𝑒𝑟𝑚𝑖𝑐𝑎
𝛿𝑡
𝜏𝑠 = 0.332𝜇∞ √𝜌 𝜇 𝜇∞ /𝑥

 Coeficiente local de fricción es:

𝜏𝑠,𝑥 1
𝐶𝑓,𝑥 = 2 = 0.664𝑅𝑒𝑥 −2
𝜇
𝜌 ∞
2
Utilizando el coeficiente local de convección
ℎ𝑥 𝑥 1
𝑁𝑢𝑥 = = 0.332 ∗ 𝑅𝑒𝑥 1/2 ∗ 𝑃𝑟 3 𝑃𝑟 ≥ 0.6
𝑘

Nusselt global
̅̅̅
ℎ𝑥 𝑥 1
̅̅̅̅̅̅𝑥 =
𝑁𝑢 = 0.664 ∗ 𝑅𝑒𝑥 1/2 ∗ 𝑃𝑟 3 𝑃𝑟 ≥ 0.6
𝑘
̅̅̅̅̅̅𝑥 = 2𝑁𝑢𝑥
𝑁𝑢
̅̅̅
ℎ𝑥 = 2ℎ𝑥
Para el flujo laminar sobre una placa isotérmica
1
0.3387 ∗ 𝑅𝑒𝑥 1/2 ∗ 𝑃𝑟 3
𝑁𝑢𝑥 =
[1 + (0.0468/Pr)2/3 ]1/4
103

Flujo Turbulento

Se asocia con la existencia de fluctuaciones aleatorias en el fluido. Es beneficioso ya que

proporciona aumento de transferencia de calor. Sin embargo el movimiento es complicado y
difícil describirle teóricamente.

𝑅𝑒𝑥 ≤ 107

𝛿 = 0.37𝑥𝑅𝑒𝑥 −1/5

𝐶𝑓,𝑥 = 0.059𝑅𝑒𝑥 −1/5 𝑅𝑒 ≤ 107

𝑁𝑢𝐿
𝑁𝑢𝑥 = 𝑆𝑡 𝑅𝑒𝑥 𝑃𝑟 = 0.0296𝑅𝑒𝑥 4/5 𝑃𝑟 1/3 𝑆𝑡 = 𝑅𝑒
𝐿 𝑃𝑟

0.6 < 𝑃𝑟 < 60

Número de Stanton (St)
𝑁𝑢𝐿
𝑆𝑡 =
𝑅𝑒𝐿 𝑃𝑟
Condiciones de capa limite mezclada
4
̅̅̅̅̅
𝑁𝑢𝐿 = (0.037𝑅𝑒𝐿 5 − 𝐴)𝑃𝑟 1/3
4
𝐴 = 0.037𝑅𝑒𝑥,𝑐 5 − 0.664𝑅𝑒𝑥,𝑐 1/2

En la capa limite mezclada:

 si se supone un numero de Reynolds de transición representativo

𝑅𝑒𝑥,𝑐 = 5 ∗ 105 : Se utiliza:
4
 ̅̅̅̅̅
𝑁𝑢𝐿 = (0.037𝑅𝑒𝐿 5 − 871)𝑃𝑟 1/3

0.6 < Pr < 60

5
⌊ 5 ∗ 10 < 𝑅𝑒𝐿 < 108 ⌋ Rango de aplicación
𝑅𝑒𝑥,𝑐 < 5 ∗ 105

Coeficiente de fricción
0.074 1742
̅̅̅̅
𝐶𝑓,𝑙 = −
1/5 𝑅𝑒𝐿
𝑅𝑒𝐿
8
5 ∗ 105 < 𝑅𝑒𝐿 ≤ 10
[ 5 ]
𝑅𝑒𝑥,𝑐 = 5 ∗ 10
Para estas condiciones:
104

4/5
Si 𝐿 ≫ 𝑥𝑐 (𝑅𝑒𝐿 ≫ 𝑅𝑒𝑥,𝑐 ); 𝐴 ≪ 0.037 𝑅𝑒𝐿
4
 ̅̅̅̅̅
𝑁𝑢𝐿 = 0.037𝑅𝑒𝐿 5 ∗ 𝑃𝑟 1/3

𝐶𝑓,𝑙 = 0.074𝑅𝑒−1/5
̅̅̅̅̅

Ejemplo 7.1

Aire a presión de 6 𝑘𝑁⁄𝑚2y a una temperatura de 300 C fluye con una velocidad de 10 𝑚⁄𝑠
sobre una placa plana de 0.5m de longitud estime la velocidad de enfriamiento por unidad de
ancho de la placa, necesaria para mantenerla a una temperatura superficial de 27 C.

Un gas a diferentes presiones

Flujo laminar
105

7.1 Propuesto

Considere los siguientes fluidos a una temperatura de película de 300 °K en un flujo paralelo sobre
una placa plana con velocidad de 1 m/s: aire atmosférico, agua, aceite de motor y mercurio.

a) Para cada fluido, determine los espesores de la capa límite de velocidad y de la capa
térmica a una distancia de 40 mm desde el borde o inicio de la placa.

𝛿𝑡

𝑚 𝛿
𝜇∞ = 1 𝑠

𝑇𝑓 = 300°𝐾

0,04 m
x

Propiedades Tabla A.4

Aire

𝑚2
(300°𝐾 1𝑎𝑡𝑚) 𝑉 = 15,89 ∗ 10−6 𝑃𝑟 = 0,707
𝑠

Tabla A.6

Agua

−6 𝑁. 𝑆
𝜇 855 ∗ 10 𝑚2 𝑚2
𝑉= = = 0,858 ∗ 10−6 𝑃𝑟 = 5,83
𝜌 𝑘𝑔 𝑠
997 3
𝑚

Tabla A.5
106

Aceite de motor

𝑚2
𝑉 = 550 ∗ 10−6 𝑃𝑟 = 6400
𝑠

Tabla A.5

Mercurio

𝑚2
𝑉 = 0,113 ∗ 10−6 𝑃𝑟 = 0,0248
𝑠

5𝑥 𝛿
𝛿= 1 𝛿𝑡 = 1
𝑅𝑒𝑥2 𝑃𝑟 3

𝑚
𝜇∞ 𝑥 1 𝑠 ∗ 0,04𝑚
𝑅𝑒𝑥 = =
𝛾 𝛾

Fluido 𝑅𝑒𝑥 𝛿(𝑚𝑚) 𝛿𝑡 (𝑚𝑚)

Aire 2517 3,99 4,48
Agua 4,66*104 0,93 0,52
Aceite de motor 72,7 23,5 1,27
Mercurio 3,54*105 0,34 1,17

7.18 Propuesto agua

107

7.14

𝑞′
𝑇𝑠 = 40℃
𝑚
𝜇∞ = 0,6 𝑠

𝑇∞ = 4℃

L = 1,5 m

a) Usando 𝑇𝑓 evaluar 𝑞 ′ y las propiedades del agua.

Tabla A.6

4 + 40
𝑇∞ = 4℃ = 277°𝐾 𝑇𝑓 = 295°𝐾 = ( ) + 273°𝐾
2

→ 𝑇𝑓 = 295°𝐾 (𝐴𝑔𝑢𝑎 𝑇𝑓 = 295°𝐾)

𝑚2 𝑤
𝛾 = 0,961 ∗ 10−6 𝐾 = 0,606 Pr = 6,62
𝑠 𝑚°𝐾

→ (𝐴𝑔𝑢𝑎 𝑇𝑠 = 40℃ = 313°𝐾)

𝑁. 𝑆
𝜇 = 657 ∗ 10−6
𝑚2

→ (𝐴𝑔𝑢𝑎 𝑇∞ = 4℃ = 277°𝐾)
108

𝑘𝑔 𝑁. 𝑆 𝜇∞ 𝑚2
𝜌∞ = 1000 ; 𝜇∞ = 1560 ∗ 10−6 ; 𝛾∞ = = 1,560 ∗ 10−6
𝑚3 𝑚2 𝜌∞ 𝑠

𝑤
𝐾∞ = 0,577 𝑃𝑟 = 11,44
𝑚°𝐾

a) 𝑞 ′ = ℎ̅ 𝐿 (𝑇𝑠 − 𝑇∞ ) ℎ̅ = ?

𝑚
𝜇∞ ∗ 𝐿 0,6 𝑠 ∗ 1,5𝑚
𝑅𝑒𝐿 = = = 9,365 ∗ 105
𝛾 𝑚2
0,961 𝑠

0,6 < 𝑃𝑟 < 60

⇒
5 ∗ 105 < 𝑅𝑒𝐿 ≤ 108 Capa limite Mezclada

4
𝑅𝑒𝑥,𝑐 = 5 ∗ 105 ̅̅̅̅̅ ⁄ 1
𝑁𝑈𝐿 = [0,037 𝑅𝑒𝐿 5 − 871] 𝑃𝑟 ⁄3

4 1
𝑁𝑈𝐿 = [ 0,037 (9,365 ∗ 105 ) ⁄5 − 871] ∗ 6,62 ⁄3 = 2522
̅̅̅̅̅

̅̅̅̅̅
𝑁𝑈𝐿 ∗ 𝐾 2522 ∗ 0,606 𝑤
̅̅̅𝐿 =
ℎ = = 1019 2
𝐿 1,5𝑚 𝑚 𝐾

𝐾𝑤
𝑞 ′ = 1019 ∗ 1,5 ∗ (40 − 4) = 55
𝑚
109

a) Con la temperatura de flujo libre (𝑇∞)

𝜇∗𝐿
𝑅𝑒𝐿∞ =
𝜇∞

𝑚
0,6 𝑠 ∗ 1,5 𝑚
𝑅𝑒𝐿∞ =
𝑚2
1,560 ∗ 10−6 𝑠

𝑅𝑒𝐿∞ = 5,769 ∗ 105

⁄ 4 1
̅̅̅̅̅
𝑁𝜇𝐿 = [0,037 𝑅𝑒𝐿 5 − 871] 𝑃𝑟 ⁄3

1
̅̅̅̅̅
𝑁𝜇𝐿 = [0,037 ∗ 5,769 ∗ 105 − 871] 11,44 ⁄3

̅̅̅̅̅
𝑁𝜇𝐿 = 1424

𝑁𝜇𝐿 ∗ 𝐾 1424 ∗ 0,057 𝑤

̅̅̅
ℎ𝐿 = = = 575
𝐿 1.5 𝑚°𝐾

𝑘𝑤
𝑞 ′ = 575 ∗ 1,5 ∗ (40 − 4) = 31,1
𝑚

Cuando se genera una capa limite turbulenta desde el inicio de la placa (alambre fino o
generador de turb.) se tiene:

⁄4 1
≫ 𝑥𝑐 ̅̅̅̅̅
𝑁𝜇𝐿 = 0,037 𝑅𝑒𝐿 5 𝑃𝑟 ⁄3

4 1
𝑁𝜇𝐿 = 0,037 (9, 365 ∗ 105 ) ⁄5 6,62 ⁄3 = 4157
̅̅̅̅̅
110

𝑁𝜇𝐿 ∗ 𝐾 4157 ∗ 0,606 𝑤

̅̅̅𝐿 =
ℎ = = 1679 2
𝐿 1,5 𝑚 °𝐾

𝑤 𝑘𝑤
𝑞 ′ = 1679 ∗ 1,5(40 − 4) = 90,7
𝑚°𝐾 𝑚

5 ∗ 105 ∗ 𝜇 𝜇 𝑚2
𝑥𝑐 = [ ]
𝜌 ∗ 𝜇∞ 𝜌 𝑠

𝜇. 𝑠
5 ∗ 105 ∗ 15,60 ∗ 106
𝑥𝑐 = 𝑚2
𝑚 𝑘𝑔
0,6 1000
𝑠 𝑚3

𝑥𝑐 = 0,0013 𝑚

𝐿 ≫ 𝑥𝑐

⁄ 4 1
̅̅̅̅̅
𝑁𝜇𝐿 = 0,037 𝑅𝑒𝐿 5 𝑃𝑟 ⁄3

Ejemplo 8.4

Un concepto que se utiliza para colectar energía solar implica la colocación de un tubo en el punto
local de un reflector parabólico y hacer pasar un fluido por el tubo.
111

El efecto neto de este arreglo se puede aproximar al de crear una condición de calentamiento
uniforme en la superficie del tubo. Es decir, se puede suponer que el flujo de calor resultante al
fluido q" es una constante a lo largo de la circunferencia y del eje del tubo. Considere la operación
con un tubo de diámetro D = 60 mm en un día soleado para el que 𝑞𝑠 "= 2000 W/𝑚2 .

1. Si entra agua presurizada al tubo a 𝑚̇ = 0.01 kg/s y Tm = 20 °C, ¿qué longitud de tubo L
se requiere para obtener una temperatura de salida de 80°C?

2. ¿Cuál es la temperatura superficial en la salida del tubo, donde se puede suponer que
existen condiciones completamente desarrolladas?

𝑚̇𝐶𝑝 (𝑇𝑚.𝑜 − 𝑇𝑚.𝑖 )

𝐴𝑠 = 𝜋𝐷𝐿 =
𝑞𝑠"

𝑚 𝐶̇ 𝑝
𝐿=
𝜋𝐷𝑞𝑠"

𝑘𝑔⁄ 𝐽
0,01 𝑠 ∗ 4181 ⁄𝑘𝑔 ∗ 𝐾
𝐿= (80 − 20)℃ = 6.65 𝑚
𝜋 ∗ 0.060𝑚 ∗ 2000 𝑊⁄ 2
𝑚
112

𝑞𝑠"
𝑇𝑠.𝑜 = + 𝑇𝑚.𝑜
ℎ

𝑘𝑔
4 𝑚̇ 4 ∗ 0.01 ⁄𝑠
𝑅𝑒𝐷 = = = 603
𝜋𝐷𝜇 𝜋 ∗ 0.060 𝑚 ∗ 352 ∗ 10−6 𝑁 𝑠⁄ 2
𝑚

ℎ𝐷
𝑁𝜇𝐷 = = 4.36
𝑘

𝑘 0.670 𝑊⁄𝑚. 𝐾
ℎ = 4.36 = 4.36 = 48.7 𝑊⁄ 2
𝐷 0.06 𝑚 𝑚 .𝐾

2000 𝑊⁄ 2
𝑇𝑠.𝑜 = 𝑚 + 80℃ = 121℃
𝑊
48.7 ⁄ 2
𝑚 .𝐾

Casos especiales

Todas las expresiones anteriores del número de Nusselt se restringen a situaciones para las que
la temperatura superficial Ts es uniforme.

 Una excepción común implica la existencia de una longitud inicial no calentada (Ts = T∞),
corriente arriba de una sección calentada (Ts ≠T∞)
113

El crecimiento de la capa limite de velocidad empieza en x = 0 mientras que la capa limite térmica
empieza en x = 𝜉.

 ∄ transferencia de calor para 0 ≤ 𝑥 ≤ 𝜉

Si 𝑞𝑠" = 𝑢𝑛𝑖𝑓𝑜𝑟𝑚𝑒

1⁄ 1⁄
𝑁𝜇𝑥 = 0,453 𝑅𝑒𝑥 2 𝑃𝑟 3 Pr ≥ 0,6

Longitud no calentada

𝑞𝑠"
𝑇𝑠(𝑥) = 𝑇∞ + 𝑞 = (𝑇𝑠 − 𝑇∞ ) ∫ ℎ 𝑑∆𝑠
ℎ𝑥

𝑞 = 𝑞𝑠" ∗ ∆𝑠

̅̅̅̅̅̅̅̅̅̅̅̅̅ 1 𝐿 𝑞𝑠" 𝐿 𝑥
(𝑇𝑠 − 𝑇∞ ) = ∫ (𝑇𝑠 − 𝑇∞ )𝑑𝑥 = ∫ 𝑑𝑥
𝐿 0 𝐿 0 𝑘𝑁𝜇𝑥

̅̅̅̅̅̅̅̅̅̅̅̅̅ 𝑞𝑠" 𝐿
(𝑇𝑠 − 𝑇∞ ) =
𝐾𝑁𝜇𝑥

⁄ 1
1
̅̅̅̅̅
𝑁𝜇𝑥 = 0,680 𝑅𝑒𝐿 2 𝑃𝑟 ⁄3
114

Ejemplo 6.5

Las pruebas experimentales sobre una parte del álabe de turbina que se muestra indican un flujo
de calor hacia la hoja de q" = 95000 W/𝑚2 ; para mantener una temperatura superficial en estado
estable de 800℃, se elimina el calor que se transfiere al álabe haciendo circular un fluido
refrigerante dentro del mismo. Determinar:

1. El flujo de calor que llega al álabe si la temperatura se reduce a 700℃ al aumentar el flujo
de fluido refrigerante.

2. El flujo de calor en la misma posición adimensional para un álabe turbina similar que tiene
una longitud de cuerda L= 80 mm, cuando el álabe opera en un flujo de aire a 𝑇∞ =
1150°C y V = 80 m/s, con Ts = 800°C.

Condición 1 (inicial)

𝑞𝑠" = ℎ1 (𝑇∞ − 𝑇𝑆1 )

𝑞1"
ℎ1 =
(𝑇∞ − 𝑇𝑆1 )

Solo se reduce la temp. En la superficie

𝑞1" = ℎ40 (𝑇∞ − 𝑇𝑆1 )

115

𝑞1"
ℎ40 =
(𝑇∞ − 𝑇𝑆1 )

Para que h = cte

𝑞2" = ℎ2 (𝑇∞ − 𝑇𝑆2 )

𝑞2"
ℎ2 =
(𝑇∞ − 𝑇𝑆2 )

ℎ = 𝑐𝑡𝑒

ℎ1 = ℎ2

𝑞2" = ℎ80 (𝑇∞ − 𝑇𝑆2 )

𝑞2"
ℎ80 =
(𝑇∞ − 𝑇𝑆2 )

ℎ40 = ℎ80

No cambia

𝑞1" 𝑞2"
=
(𝑇∞ − 𝑇𝑆1 ) (𝑇∞ − 𝑇𝑆2 )

𝑇∞ − 𝑇𝑆2
𝑞2" = 𝑞1" ∗
𝑇∞ − 𝑇𝑆1

(1150 − 700)
𝑞2" = 95000 ∗
(1150 − 800)
116

𝑤
𝑞2" = 122000
𝑚2

Para que h = cte con las condiciones de cambio de temperatura en la superficie.

𝑉2 𝐿2
𝑅𝑒𝐿2 =
𝛾

160 ∗ 40
𝑅𝑒𝐿1 =
𝛾

6400
𝑅𝑒𝐿1 =
𝛾

𝑉1 𝐿1
𝑅𝐿1 =
𝛾

80 ∗ 80
𝑅𝐿2 =
𝛾

6400
𝑅𝐿2 =
𝛾
117

Para el álabe más largo y el flujo de aire reducido

⇒ 𝑅𝑒𝐿1 = 𝑅𝑒𝐿2

Pr = 𝑐𝑡𝑒

⇒ 𝑁𝜇1 = 𝑁𝜇2

ℎ1 𝐿1 ℎ2 𝐿2
=
𝑘1 𝑘2

𝐿1
ℎ2 = ℎ
𝐿2

𝑞1" = ℎ1 (𝑇∞ − 𝑇𝑆 )

𝑞1"
ℎ1 =
(𝑇∞ − 𝑇𝑆 )

𝑞1" 𝐿1
ℎ2 =
(𝑇∞ − 𝑇𝑆 ) 𝐿2

𝑞2" = ℎ2 (𝑇∞ − 𝑇𝑆 ) ⇒

𝑞1" 𝐿1
𝑞2" = ∗ (𝑇∞ − 𝑇𝑆 )
(𝑇∞ − 𝑇𝑆 ) 𝐿2

𝐿1
𝑞2" = 𝑞1" ∗ 𝐿2

𝑤 0,04
𝑞2" = 95000 ∗
𝑚2 0,08

𝑤
𝑞2" = 47500
𝑚2

Flujo Interno

Flujo laminar en un tubo circular de radio ro donde el fluido entra con una velocidad
uniforme
 118

𝜇𝑚 : Velocidad media del fluido sobre la sección transversal del tubo

𝐷: Diámetro del tubo

Reynolds critico en un flujo completamente desarrollado

Reynolds critico al inicio del tubo

Para flujo interno es
importante conocer la
extensión de la región
de entrada que
depende de si el flujo es
laminar o turbulento

Para condiciones completamente turbulentas 𝑅𝑒𝐷 ≈ 10000

Flujo laminar turbulentas 𝑅𝑒𝐷 ≤ 2300

La longitud hidrodinámica de entrada se da por la expresión

Suponemos que el fluido entra

al tubo desde una boquilla por
lo que el perfil de velocidad es
casi uniforme
119

𝜋 𝐷2
Para un flujo en tubo circular 𝐴𝑐 = 4

4 𝑚̇
𝑅𝑒𝐷 =
𝜋𝐷𝜇

Gradiente de presión y factor de fricción en un flujo completamente desarrollado

El ing. se encuentra siempre con caídas de presión y este parámetro determina los
requerimientos de potencia de bombeo o ventilación.

Para determinar la caída de presión se trabaja con el factor de fricción de Moody (o de Darcy).

Parámetro adimensional

𝑑𝑝
− ( ⁄𝑑𝑥 )𝐷
𝑓𝑓𝑎𝑐𝑡𝑜𝑟 𝑑𝑒 𝑓𝑟𝑖𝑐𝑐𝑖ó𝑛 = No se confunde con el
𝜇2 𝑚
𝜌 coeficiente de fricción.
2

Factor de fricción de Fanning

𝑇𝑠
𝐶𝑓 =
𝜇2 𝑚
𝜌 2

𝑇𝑠 = − 𝜇 (𝑑𝑢⁄𝑑𝑟)
𝑟=𝑟0

𝑓
𝐶𝑓 =
4

Para un flujo completamente desarrollado (Laminar)

64
𝑓=
𝑅𝑒𝐷
120

Para un flujo turbulento completamente desarrollado, el análisis es más complicado, y debemos

depender finalmente de resultados experimentales. Se utiliza el diagrama de Moody de la figura
8.3, se presentan factores de fricción para un intervalo del numero de Reynolds.

El factor de fricción es una función de la condición de la superficie del tubo. Es un mínimo para
superficies suaves y se incrementa al aumentar la rugosidad de la superficie, “e”.

−1⁄
4
𝑓 = 0,316 𝑅𝑒𝐷 𝑅𝑒𝐷 ≲ 2 ∗ 104

−1⁄
5
𝑓 = 0,184 𝑅𝑒𝐷 𝑅𝑒𝐷 ≳ 2 ∗ 104

Para números grandes de Reynolds (Expresión Petukhov)

𝑓 = (0,790 ln 𝑅𝑒𝐷 − 1,64)−2 3000 ≤ 𝑅𝑒𝐷 ≤ 5 ∗ 106

El factor de fricción y dp/dx = cte

𝑑𝑝
𝑓, 𝑦 𝑑𝑥
= 𝑐𝑡𝑒 para la región completamente desarrollada.

𝑑𝑝
( )𝐷
𝑑𝑥
⇒ ∆𝑃 = 𝑃1 − 𝑃2 caída de presión de la expresión 𝑓 ≡ 𝜇 𝑚2
𝜌
2

Para una posición axial desde 𝑥1 𝑎 𝑥2

𝑃2
𝜌𝜇 𝑚2 𝑥2 𝜌𝜇 𝑚2
∆𝑃 = − ∫ 𝑑𝑝 = 𝑓 ∫ 𝑑𝑥 = 𝑓 (𝑥2 − 𝑥1 )
𝑃1 2𝐷 𝑥1 2𝐷

64
𝑓=
𝑅𝑒𝐷

La potencia P en W requerida para vencer la resistencia al flujo asociada con la caída de

121

presión:

Ρ = (∆𝑝)∀̇

∀̇  Flujo volumétrico
𝑚̇
∀̇= 𝜌
[Fluido incompresible]

Entra fluido al tubo a una temperatura uniforme 𝑇(𝑟,0) que es menor que la temperatura de la
superficie Ts  𝑞𝑐𝑜𝑛𝑣 y empieza a producir una capa limite térmica.

Ahora si en a superficie del tubo nos imponemos una Ts=cte o flujo de calor uniforme 𝑞𝑠" = 𝑐𝑡𝑒
 se alcanza una capa térmica completamente desarrollada .

La forma del perfil de temperatura completamente desarrollada 𝑇(𝑟,𝑥) es variable según se

mantenga una Ts=cte o 𝑞𝑠" = 𝑐𝑡𝑒

Para el flujo laminar la longitud de entrada térmica se puede expresar como:

𝑋𝑐𝑑,𝑡
( ) ≈ 0,05 𝑅𝑒𝐷 𝑃𝑟
𝐷 𝑙𝑎𝑚

Si Pr > 1 La capa limite hidro se desarrolla mas rápido que la capa limite térmica
122

𝑋𝑐𝑑ℎ < 𝑋𝑐𝑑𝑡

Si Pr > 1 lo inverso es cierto

Pr en aceites es extremadamente grande Pr ≳ 1000 𝑋𝑐𝑑ℎ ≪ 𝑋𝑐𝑑𝑡

⇒ El perfil de velocidad se desarrolla completamente en la región de entrada térmica.

Consideraciones generales

Como el flujo en un tubo está encerrado se puede explicar un balance de energía para determinar
cómo varia la temperatura media 𝑇𝑚(𝑥) con la posición a lo largo del tubo, y como esta
relacionada la transferencia total de calor por convección 𝑞𝑐𝑜𝑛𝑣 con la diferencia de temperaturas
en la entrada y salida del tubo.

La temperatura media 𝑇𝑚 , es una temperatura de referencia conveniente para flujos internos

que desempeña una función muy similar a la de la temperatura de flujo libre 𝑇∞ para flujo
externo.

⇒ La ley de enfriamiento de Newton:

(1) 𝑞𝑠" = ℎ (𝑇𝑠 − 𝑇𝑚 )

𝑇𝑚 debe variar en la dirección del flujo

𝑑 𝑇𝑚
≠ 0 𝑝𝑎𝑟𝑎 𝑞𝑢𝑒 𝑜𝑐𝑢𝑟𝑟𝑎 𝑙𝑎 𝑡𝑟𝑎𝑛𝑠𝑓𝑒𝑟𝑒𝑛𝑐𝑖𝑎 𝑑𝑒 𝑐𝑎𝑙𝑜𝑟
𝑑𝑥

Si 𝑇𝑆 > 𝑇𝑚 ⇒ 𝑆𝑒 𝑡𝑟𝑎𝑛𝑠𝑓𝑖𝑒𝑟𝑒 𝑐𝑎𝑙𝑜𝑟 𝑎𝑙 𝑓𝑙𝑢𝑖𝑑𝑜 𝑦 𝑇𝑚 𝑎𝑢𝑚𝑒𝑛𝑡𝑎 𝑐𝑜𝑛 𝑥

Si 𝑇𝑆 < 𝑇𝑚 ⇒ 𝐿𝑜 𝑜𝑝𝑢𝑒𝑠𝑡𝑜 𝑒𝑠 𝑐𝑖𝑒𝑟𝑡𝑜

123

𝐸𝑛𝑡𝑟𝑎 + 𝐸𝑔 − 𝐸𝑠 ≡ 𝐸𝑎𝑙𝑚

𝑑(𝐶𝑣 𝑇𝑚 + 𝑃𝑣)
(2) 𝑑𝑞𝑐𝑜𝑛 + 𝑚̇ (𝐶𝑣 𝑇𝑚 + 𝑃𝑣) − [𝑚̇ (𝐶𝑣 𝑇𝑚 + 𝑃𝑣) + 𝑚̇ 𝑑𝑥 ] = 0
𝑑𝑥

(3) 𝑑𝑞𝑐𝑜𝑛𝑣 = 𝑚̇𝑑 (𝐶𝑣 𝑇𝑚 + 𝑃𝑣)

Si suponemos que el fluido es un gas ideal

𝑃𝑣 = 𝑅𝑇𝑚 ; 𝐶𝑝 = 𝐶𝑣 + 𝑅

la Ecuación (3) queda:

(4)𝑑𝑞𝑐𝑜𝑛𝑣 = 𝑚̇𝐶𝑝𝑑𝑇𝑚

Una forma especial de la ecuación (4) se relaciona con las condiciones para todo el tubo. En
particular al ingresar de la entrada i, a la salida O del tubo:
124

(5) 𝑞𝑐𝑜𝑛𝑣 = 𝑚 𝐶𝑝 (𝑇𝑚, ̇ 𝑜 − 𝑇𝑚, 𝑖)

∃ 3 𝑖𝑚𝑝𝑜𝑟𝑡𝑎𝑛𝑡𝑒𝑠 𝑣𝑎𝑟𝑖𝑎𝑙𝑏𝑙𝑒𝑠 𝑡𝑒𝑐𝑛𝑖𝑐𝑎𝑠:

𝑞𝑐𝑜𝑛𝑣 ; 𝑇𝑚, 𝑜; 𝑇𝑚, 𝑖

Es una expresión que se aplica independientemente de las condiciones térmicas de la superficie

o de las condiciones de flujo, de la ecuación (4) 𝑑𝑞𝑐𝑜𝑛𝑣 = 𝑚̇𝐶𝑝 𝑑 𝑇𝑚

La transferencia de calor para cada elemento diferencial

(6) 𝑑𝑞𝑐𝑜𝑛𝑣 = 𝑞𝑠" 𝑃 𝑑𝑥; 𝑃 = 𝜋𝐷 (𝑡𝑢𝑏𝑜)

𝑑𝑇𝑚 𝑞𝑆" 𝑃 𝑃
= = ℎ (𝑇𝑠 − 𝑇𝑚)
𝑑𝑥 𝑚̇𝐶𝑝 𝑚̇𝐶𝑝

→ 𝑞𝑐𝑜𝑛𝑣 = 𝑚̇𝐶𝑝 𝑑 𝑇𝑚 = 𝑞𝑠" 𝑃 𝑑𝑥 𝑞𝑠" = ℎ (𝑇𝑠 − 𝑇𝑚)

𝑑𝑇𝑚 𝑞𝑠" 𝑃 (𝑇𝑠 − 𝑇𝑚)𝑃

(7) = =ℎ
𝑑𝑥 𝑚̇𝐶𝑝 𝑚̇𝐶𝑝

(4) La expresión (7) es extremadamente útil.

La solución de (7) para Tm(x) depende de la condición térmica de la superficie.

 1: Flujo de calor superficial cte.

2: Temperatura superficial cte.

FLUJO DE CALOR SUPERFICIAL CTE 𝒒"𝒔 𝑷 = 𝒄𝒕𝒆

Para flujo de calor superficial cte observamos que es sencillo determinar 𝑞𝑐𝑜𝑛𝑣 . Como 𝑞𝑠" es
125

independiente de X

⇒ 𝑞𝑐𝑜𝑛𝑣 = 𝑞𝑠" (𝑃 ∗ 𝐿) (8)

Se puede utilizar la ecuación (5) para determinar el cambio en la temperatura del fluido Tm,o –
Tm,i

Para 𝑞𝑠" = 𝑐𝑡𝑒 se puede usar la ecuación (7)

𝑑𝑇𝑚 𝑞𝑠" 𝑃
= ≠ 𝑓(𝑥)
𝑑𝑥 𝑚̇𝐶𝑝

Al separar variables e interpretar desde x=0

Condiciones
𝑞𝑠" 𝑃 superficiales
(9) 𝑇𝑚(𝑥) = 𝑇𝑚, 𝑖 + 𝑥 𝑞𝑠" = 𝑐𝑡𝑒 constantes
𝑚̇𝐶𝑝

En consecuencia la temperatura media, varia en forma lineal con x a lo largo del tubo.

Ts – Tm inicialmente es pequeña (debido

al valor grande de h en la entrada). Al
aumentar x el valor de h disminuye
porque se desarrolla la capa limite.

Pero en la R.C.D h es indepediente de x.

Por lo tanto Ts – Tm es independiente de

x en la R.C.D.
126

TEMPERATURA SUPERFICIAL CONSTANTE Ts = cte

Los resultados para T.C total y distribución axial de la temperatura son completamente diferentes
para Ts=cte.

Definido ∆𝑇 = 𝑇𝑠 − 𝑇𝑚 la ecuación (7) se puede expresar como:

𝑑 𝑇𝑚 𝑑(∆𝑇) 𝑃
(10) = − = ℎ ∆𝑇
𝑑𝑥 𝑑𝑥 ̇
𝑚 𝐶𝑝

Separando variables e integrando desde la entrada hasta la salida del tubo:

∆𝑇𝑜 𝐿
𝑑(∆𝑇) 𝑃
∫ = − ∫ ℎ𝑑𝑥
∆𝑇𝑖 ∆𝑇 𝑚̇𝐶𝑝 0

𝐿
∆𝑇𝑜 𝑃 1
ln ( )= − (∫ ℎ𝑑𝑥 )
∆𝑇𝑖 𝑚̇𝐶𝑝 0 𝐿

El coeficiente promedio de convección se expresa como:

1 1 𝐿
ℎ̅ = = ∫ ℎ𝑑𝐴𝑠; ℎ̅ = ∫ ℎ𝑑𝑥
∆𝑠 ∆𝑠 𝐿 0

∆𝑇𝑜 𝑃𝐿
(11) ln ( )= − ̅̅̅
ℎ
∆𝑇𝑖 𝑚̇𝐶𝑝 𝐿

∆𝑇𝑜 𝑇𝑠 − 𝑇𝑚, 𝑜 𝑃𝑥 ̅
(− ℎ)
(12) = = 𝑒 𝑚̇𝐶𝑝 𝑇𝑠 = 𝑐𝑡𝑒
∆𝑇𝑖 𝑇𝑠 − 𝑇𝑚, 𝑖

(5) ℎ̅ es el valor promedio desde x=0 (entrada del tubo) hasta un valor x
127

⇒ 𝑇𝑠 − 𝑇𝑚 disminuye exponencialmente.

La expresión para la T.C se complica por la naturaleza exponencial de la disminución de la

Temperatura de la ecuación (5) 𝑞𝑐𝑜𝑛𝑣 = 𝑚̇ 𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖).

Reacomodada:

𝑞𝑐𝑜𝑛𝑣 = 𝑚̇ 𝐶𝑝 [(𝑇𝑠 − 𝑇𝑚, 𝑖) − (𝑇𝑠 − 𝑇𝑚, 𝑜)] = 𝑚̇ 𝐶𝑝 (∆𝑇𝑖 − ∆𝑇𝑜)

𝑞𝑐𝑜𝑛𝑣 = 𝑚̇𝐶𝑝 (∆𝑇𝑖 − ∆𝑇𝑜)

(13) 𝑞 𝑐𝑜𝑛𝑣 = ℎ̅ 𝐴𝑠 ∆𝑇𝑚𝑙 𝑇𝑠 = 𝑐𝑡𝑒

𝑑𝑜𝑛𝑑𝑒 𝐴𝑠 = 𝑃 ∗ 𝐿

∆𝑇𝑚𝑙 = 𝑇𝑒𝑚𝑝𝑒𝑟𝑎𝑡𝑢𝑟𝑎 𝑚𝑒𝑑𝑖𝑎 𝑙𝑜𝑔𝑎𝑟𝑖𝑡𝑚𝑖𝑐𝑎

∆𝑇𝑜 − ∆𝑇𝑖
(14) ∆𝑇𝑚 =
∆𝑇𝑜
ln ( ∆𝑇𝑖 )

Para algunas implicaciones en la temperatura del fluido externo la que se fija en lugar de la
temperatura de la superficie del tubo y el valor de ℎ̅ se reemplaza por 𝑈
̅.
128

𝑈̅ 𝐴𝑠
∆𝑇𝑜 𝑇∞ − 𝑇𝑚, 𝑜 (− ̇
𝑚 𝐶𝑝
)
(15) = = 𝑒
∆𝑇𝑖 𝑇∞ − 𝑇𝑚, 𝑖

̅𝐴𝑠 ∆𝑇𝑚𝑙
(16) 𝑞 = 𝑈

Ejemplo 8.2

Un sistema para calentamiento de agua desde una temperatura de entrada Tm,i = 20°C a una
temperatura de salida Tm,o = 60°C implica hacer pasar el agua por un tubo de pared delgada que
tiene diámetros interno y externo de 20 y 40 mm. La superficie externa del tubo está bien aislada
y el calentamiento eléctrico dentro de la pared proporciona una generación uniforme de 𝑞̇ =
106 𝑊⁄ 3 .
𝑚

a) Para un flujo de masa de agua 𝑚̇ = 0 l kg/seg, ¿que tan largo debe ser el tubo para alcanzar
la temperatura de salida que se deseada?
a) Si la temperatura de la superficie interna del tubo es Ts=70℃ en la salida; cual es el
coeficiente local de transferencia de calor por convección en la salida.
129

𝑇𝑚, 𝑖 + 𝑇𝑚, 𝑜 (20 + 60)℃

̅̅̅̅
𝑇𝑚 = ; ̅̅̅̅
𝑇𝑚 = ; ̅̅̅̅
𝑇𝑚 = 40℃; ̅̅̅̅
𝑇𝑚 = 313°𝐾
2 2

𝐾𝐽
𝐶𝑝𝑎𝑔𝑢𝑎 = 4,179
𝑘𝑔 °𝐾

∄ 𝑇. 𝐶 𝑎 𝑙𝑎 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒 𝑒𝑥𝑡𝑒𝑟𝑛𝑎 𝑝𝑜𝑟 𝑒𝑙 𝑐𝑎𝑙𝑒𝑛𝑡𝑎𝑚𝑖𝑒𝑛𝑡𝑜

𝐸𝑒𝑛𝑡𝑟𝑎 − 𝐸𝑠 + 𝐸𝑔 ≡ 𝐸𝑎𝑙𝑚

𝐸𝑔 = 𝐸𝑠

′′′
𝐸𝑔 = 𝑞𝑐𝑜𝑛𝑣

𝜋
(6) 𝐴 = (𝐷𝑜 2 − 𝐷𝑖 2 )
4

̇
𝐸𝑔 = 𝑞𝑐𝑜𝑛𝑣 ∗𝑉

𝜋
𝐸𝑔 = 𝑞̇ (𝐷𝑜 2 − 𝐷𝑖 2 ) ∗ 𝐿
4

𝑞𝑐𝑜𝑛𝑣 = 𝑚̇𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖) 𝑞𝑆" = 𝑐𝑡𝑒

𝐸𝑔 = 𝑞𝑐𝑜𝑛𝑣

𝜋
𝑞̇ (𝐷𝑜 2 − 𝐷𝑖 2 ) ∗ 𝐿 = 𝑚̇𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖)
4

4 𝑚̇ 𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖)
⇒ 𝐿=
𝜋 (𝐷𝑜 2 − 𝐷𝑖 2 )𝑞̇
130

𝐿 = 17,74 𝑚

b) 𝑞𝑆" = ℎ𝑜 ∆𝑇 Para la salida del tubo; ho es coeficiente local en la salida.

En la salida x=L.

∆𝑇 en la superficie interna de la salida del tubo

∆𝑇 = (𝑇𝑠, 𝑜 − 𝑇𝑚, 𝑜) 𝑞𝑆" = ℎ𝑜 (𝑇𝑠, 𝑜 − 𝑇𝑚, 𝑜)

𝜋
𝐸𝑔 = 𝑞̇ (𝐷𝑜 2 − 𝐷𝑖 2 )𝐿 [𝑊]
4

𝑊
𝑞𝑆" = [ ]⇒
𝑚2

𝜋
𝐸𝑔 𝑞̇ 4 (𝐷𝑜 2 − 𝐷𝑖 2 )𝐿
𝑞𝑆" = 𝑞𝑆" =
𝐴𝑟𝑒𝑎 𝜋 𝐷𝑖 𝐿

𝑞̇ (𝐷𝑜 2 − 𝐷𝑖 2 ) 𝑊
𝑞𝑆" = ; 𝑞𝑆" = 1,5 ∗ 104
4 𝐷𝑖 𝑚2

𝑞𝑆" = ℎ𝑜 ( 𝑇𝑠, 𝑜 − 𝑇𝑚, 𝑜)

𝑞𝑆" 𝑞̇ (𝐷𝑜2 − 𝐷𝑖 2 )
ℎ𝑜 = ; ℎ𝑜 =
𝑇𝑠, 𝑜 − 𝑇𝑚, 𝑜 𝐷𝑖 (𝑇𝑠, 𝑜 − 𝑇𝑚, 𝑜)

𝑊
ℎ𝑜 = 1500
𝑚2 °𝐾

Ejemplo 8.3

La condensación de vapor sobre la superficie externa de un tubo circular de pared delgada de 50

mm de diámetro y 6 m de longitud mantiene una temperatura superficial uniforme de 100°C. Por
el tubo fluye agua a razón de 𝑚̇= 0.25 kg/s, y sus temperaturas de entrada y de salida son Tm,i =
15°C y Tm,o = 57°C. ¿Cuál es el coeficiente promedio de convección asociado con el flujo de agua?
131

15 + 57
̅̅̅̅̅ =
𝑇𝑚 = 36℃ → 309°𝐾
2

𝐽
𝐶𝑝𝑎𝑔𝑢𝑎 = 4178
𝑘𝑔°𝐾

𝑞𝑐𝑜𝑛𝑣 = 𝑚̇ 𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖)

𝑞𝑐𝑜𝑛𝑣 = ℎ̅ 𝐴𝑠 ∆𝑇𝑚𝑙 → 𝑝𝑎𝑟𝑎 𝑇𝑠 = 𝑐𝑡𝑒

𝑚̇𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖) = ℎ̅ 𝐴𝑠 ∆𝑇𝑚𝑙

𝑚̇𝐶𝑝 (𝑇𝑚, 𝑜 − 𝑇𝑚, 𝑖)

ℎ̅ = 𝐴𝑠 = 𝜋𝐷𝐿
𝐴𝑠 ∆𝑇𝑚𝑙

∆𝑇𝑜 − ∆𝑇𝑖
∆𝑇𝑚𝑙 =
∆𝑇𝑜
ln ( ∆𝑇𝑖 )

(𝑇𝑠 − 𝑇𝑚, 𝑜) − (𝑇𝑠 − 𝑇𝑚, 𝑖) (100 − 57) − (100 − 15)

∆𝑇𝑚𝑙 = ∆𝑇𝑚𝑙 =
𝑇𝑠 − 𝑇𝑚, 𝑜 100 − 57
ln ( 𝑇𝑠 − 𝑇𝑚, 𝑖 ) ln ( )
100 − 15

∆𝑇𝑚𝑙 = 61,63℃

𝐽
4178 ⁄𝑘𝑔°𝐾 ∗ 0,25 𝑘𝑔 ∗ (57 − 15)
ℎ̅ =
5
132

𝑊
ℎ̅ = 755,22
𝑚2 °𝐾

(7) FLUJO LAMINAR EN TUBOS CIRCULARES

ANÁLISIS TÉRMICO Y CORRELACIONES DE CONVECCIÓN

Es necesario conocer los coeficientes de convección, por lo que se explicará la manera teórica
de obtenerlos para el flujo laminar.

REGIÓN COMPLETAMENTE DESARROLLADA

En cualquier punto del tubo se suponen aplicables las aproximaciones de capa limite.

Combinando la ley de enfriamiento de Newton.

𝑑𝑇𝑚 𝑞𝑆" 𝑃
(3) 𝑞𝑆" = ℎ (𝑇𝑠 − 𝑇𝑚) 𝑦 𝑙𝑎 𝑠𝑜𝑙𝑢𝑐𝑖ó𝑛 𝑑𝑒 (4) =
𝑑𝑥 𝑚̇𝐶𝑝

48 𝐾 11 𝑞𝑆" 𝐷
ℎ= 𝑇𝑚 − 𝑇𝑠 = (5)
11 𝐷 48 𝐾

Si 𝑞𝑆" = 𝑐𝑡𝑒

ℎ𝐷
(6) 𝑁𝑢𝐷 = = 4,36
𝐾

Por lo tanto en un tubo circular caracterizado por flujo de calor superficial uniforme y condiciones
completamente desarrolladas 𝑁𝑢𝐷 es una constante, independiente de 𝑅𝑒𝐷 , 𝑃𝑟 y la posición
axial.
133

 Para condiciones laminares completamente desarrolladas con una temperatura

superficial constante las aproximaciones de capa limite de velocidad o hidrodinámica se
𝜕2 𝑇
satisfacen otra vez exactamente y la aproximación de capa limite térmica (𝜕𝑥 2 ) ≪
𝜕2 𝑇
( ) a menudo es razonable.
𝜕𝑟 2

El perfil de velocidad se expresa como:

𝜇(𝑟) 𝛾 2
= 2 [1 − ( ) ]
𝜇𝑚 𝛾𝑜

Para Ts = cte

𝑑𝑇𝑠
=0
𝑑𝑥

𝜕𝑇 (𝑇𝑠 − 𝑇) 𝑑𝑇𝑚
| = |
𝜕𝑥 𝑐𝑑,𝑡 (𝑇𝑠 − 𝑇𝑚) 𝑑𝑥 𝐶𝑑,𝑡

Combinando estas expresiones:

Sii Ts = cte:

1 𝑑 𝑑𝑇 2 𝜇𝑚 𝑑𝑇𝑚 𝛾 2 𝑇𝑠 − 𝑇
(𝛾 ) = ( ) [1 − ( ) ]
𝛾 𝑑𝛾 𝑑𝛾 𝛼 𝑑𝑥 𝛾𝑜 𝑇𝑠 − 𝑇𝑚

Para resolver esta ecuación se puede utilizar un procedimiento iterativo que implica hacer
aproximaciones sucesivas al perfil de temperaturas.

El resultado no se lo puede describir mediante una simple expresión algebraica. Sin embargo se
puede mostrar que el número de Nusselt que se obtiene es:

Sii Ts = cte

ℎ𝐷
𝑁𝑢𝐷 = = 3,66
𝐾
134

Revisar  Tarea ejemplo 8.4

REGIÓN DE ENTRADA

Para resolver la Ecuación de Energía:

𝜇𝜕𝑇 𝑉𝜕𝑇 𝛼 𝜕 𝛾𝜕𝑇

+ = ( )
𝜕𝑥 𝜕𝛾 𝛾 𝜕𝛾 𝜕𝛾

En la región de entrada es complicado, ya que la velocidad y temperatura dependen de x así como

𝜕𝑇
también de 𝛾. El gradiente de temperatura 𝜕𝑥 tal vez ya no se pueda simplificar mediante

𝜕𝑇 𝑑𝑇𝑚
| = | 𝑞𝑆" = 𝑐𝑡𝑒
𝜕𝑥 𝑐,𝑑𝑡 𝑑𝑥 𝑐,𝑑𝑡

𝜕𝑇 𝑇𝑠 − 𝑇 𝑑𝑇𝑚
| = | 𝑇𝑠 = 𝑐𝑡𝑒
𝜕𝑥 𝑐,𝑑𝑡 𝑇𝑠 − 𝑇𝑚 𝑑𝑥 𝑐,𝑑𝑡

A pesar de esto se obtienen dos diferentes soluciones de longitud de entrada. La solución mas
sencilla es para el problema de longitud de entrada térmica, y se basa en la suposición de que las
condiciones térmicas se generan en presencia de un perfil de velocidad completamente
desarrollado.

Esta situación existirá si la posición a la que comienza la T.C. estuviera precedida por una longitud
inicial no calentada o para fluidos con Pr muy alto como los aceites
135

Los números de 𝑁𝑢𝐷 en principio don infinitos para x = 0 y disminuyen a sus valores asintóticos
(completamente desarrollados) al aumentar x.

Conforme Pr  ∞ las condiciones completamente desarrolladas se alcanzan para

𝑥
(𝐷 )
[ ] ≈ 0,05
𝑅𝑒𝐷 𝑃𝑟

Para la condición de Ts = cte se puede determinar ̅̅̅

ℎ𝑜 a partir de la correlación de Hausen.

Longitud de
𝐷 entrada térmica
0,0668 ( 𝐿 ) 𝑅𝑒𝐷 𝑃𝑟
̅̅̅̅̅̅
𝑁𝑢 𝐷 = 3,66 + 2⁄
𝐷 3
1 + 0,04 [( 𝐿 ) 𝑅𝑒𝐷 Pr]

Para la longitud de entrada combinada (Térmica y de velocidad). La correlación adecuada es la

de Sieder y Tate.

1⁄
3
𝑅𝑒𝐷 𝑃𝑟 𝜇 0,14
̅̅̅̅̅̅
𝑁𝑢 𝐷 = 1,86 [ ] ( )
𝐿 𝜇𝑆
( )
𝐷

𝑇𝑠 = 𝑐𝑡𝑒

0,48 < Pr < 16,700

𝜇
0,0044 < ( ) < 9,75
𝜇𝑆

𝐿
10 ≲ ≤ 30
𝐷

𝑇𝑚, 𝑖 + 𝑇𝑚, 𝑜
𝜇 𝑒𝑙𝑒𝑣𝑎𝑑𝑎 𝑎 𝑇𝑚 =
2

𝜇𝑆 𝑒𝑙𝑒𝑣𝑎𝑑𝑎 𝑎 𝑇𝑠
136

(9) CORRELACIONES DE CONVECCIÓN PARA FLUJO TURBULENTO EN TUBOS CIRCULARES

El análisis es mucho más complicado, por lo tanto se hace mas énfasis en correlaciones empíricas.

Una expresión clásica para calcular 𝑁𝑢𝐷 local para F.T.C.D (hidrodinámica y térmicamente) en un
tubo circular suave se debe a Chilton-Colburn

2
𝐶𝑓 𝑓 𝑁𝑢𝐷 2
= = 𝑆𝑡 𝑃𝑟3 = 𝑃𝑟 3
2 8 𝑅𝑒𝐷 𝑃𝑟

Si 𝑓 = (0,790 ln 𝑅𝑒𝐷 − 1,64)−2

4 1
𝑁𝑢𝐷 = 0,023 𝑅𝑒𝐷5 𝑃𝑟3

De manera general Dittus-Boelter proponen la siguiente correlación

4
𝑁𝑢𝐷 = 0,023 𝑅𝑒𝐷5 𝑃𝑟𝑛

Donde:

n= 0,4 para calentamiento (𝑇𝑠 > 𝑇𝑚)

n= 0,3 para enfriamiento (𝑇𝑠 < 𝑇𝑚)

Experimentalmente confirmado para las siguientes condiciones

0,7 ≤ Pr ≤ 160
137

𝑅𝑒𝐷 ≳ 10000

𝐿
≳ 10
𝐷

Esta expresión se debe usar solo para diferencias de temperaturas de pequeñas a moderada
(Ts – Tm) con todas las propiedades evaluadas a Tm.

Para variaciones grandes de las propiedades se puede utilizar la correlación de Sieder y Tate.

4 1
𝜇 0,14
𝑁𝑢𝐷 = 0,027 𝑅𝑒𝐷5 𝑃𝑟3 ( )
𝜇𝑆

0,7 ≤ Pr ≤ 16700
Todas las propiedades
evaluadas a Tm excepto
𝑅𝑒𝐷 ≥ 10000 𝜇𝑆

𝐿
≳ 10
𝐷

Con estas expresiones se podrían en determinados casos incidir en errores al 25%.

Mediante el uso de correlaciones mas recientes se puede reducir los errores en al menos un
10%.

Por ejemplo Petukhov propone:

𝑓
(8) 𝑅𝑒𝐷 𝑃𝑟
𝑁𝑢𝐷 = 1 2
𝑓 2
1,07 + 12,7 (8) (𝑃𝑟3 − 1)

Valida para:
138

0,5 < Pr < 2000

104 < 𝑅𝑒𝐷 < 5 ∗ 106

f se puede obtener del diagrama de Moody o para tubos suaves 𝑓 = (0,790 ln 𝑅𝑒𝐷 − 1,64)−2

Para numros de Reynolds pequeños Gnielinsky modifica la correlación de Petakhov y propone:

𝑓
(8) (𝑅𝑒𝐷 − 1000)𝑃𝑟
𝑁𝑢𝐷 = 1 2
𝑓 2
1 + 12,7 (8) (𝑃𝑟3 − 1)

Valida para:

0,5 < Pr < 2000

3000 < 𝑅𝑒𝐷 < 5 ∗ 106

Tarea
(10) Ejercicios 8.38 y 8.23 En clase
Ejemplos 8.4/8.5

Propuestos 8.38 b
Tarea
8.39
Ejemplos 8.4/8.5
8.44 a

Propuestos 8.34
Consultar
8.38 (b)
Correlaciones de convección en tubos
8.39
no circulares.
8.41

8.43
Ejemplo 8.6
8.54

8.78
8.45  8.38

8.46

Ver la correspondencia en Incropera.

139

Ejemplo 8.6

Consulta

Correlaciones de convección en tubos no circulares

𝟒𝑨𝒄
𝑫𝒏 ≡
𝒑

Ejemplo 8.6

Anillos en tubos concentricos