CAP.

1O: BALANDE DE ENERGIA EN LA ENVOLTURA Y DISTRIBUCION

DE TEMPERATURA EN SOLIDO Y EN FLUJO LAMINAR

10B.2 Calentamiento viscoso en flujo en una rendija. Encontrar el perfil de temperatura para
el problema de calentamiento viscoso que se muestra en la figura 10.4-2, cuando se
proporcionan las siguientes condiciones lmite: para

x 0, T T0

; para

x b, qx 0

Solucin:
Las ecuaciones 10.4-5 y 10.4-6 son todava vlidas para este problema.
2

10.4 5

dT

x C1
dx
b

en la que k es la conductividad calorfica del fluido.

2 2
C
x
T
1 x C2
10.4 6
k
k b 2
Las constantes de integracin
C.L.1:

C1

para x = 0

C2

se determinan a partir de las condiciones lmite:

T T0

dT
0
dx

qx k

C.L.2:
para x = b
Luego reemplazando en la ecuacin 10.4-5 la segunda condicin lmite se tiene:
2

0 b C1
1
b
2

b C1
b

Y reemplazando la primera condicin lmite en la ecuacin 10.4-6, obtenemos:

2
2
C
0
T0
1 0 C2
2
k
k b
2 T0 C2
Reemplazando las expresiones (1) y (2) por las constantes de integracin en 10.4-6



k b

x 2 C1
x C2
2
k

k b

x2
2

x T0
2
bk

T
T

2 x 2 2
x
2
bk
k 2b

T T0

2
x2 x

2b 2 b
k

T T0

T T0 x 1 x 2

2 / k b 2 b 2

T T0 x
1 x

2
/ k b
2 b

10B.18 Perfiles de temperatura en un reactor con densidad de flujo axial (figura 10B.18).

a) Demostrar que para una fuente de calor que depende linealmente de la temperatura, las
ecuaciones 10.5-6 a 10.5-14 tienen las soluciones (para
m m em em
e m m Z
I Z 1
2 m
2 m
m
e

m
e

10 B.18 1

m m

 10B.18 2

m e m e m Z m e m e m Z
m m
m2 e m m2 em

II Z

m2 m2 m m
e
2 m
2 m
m e m e

III Z 1

10 B.18 3
m

B
1 1 4N / B
2

Aqu
, donde. En la figura 10B.18 se muestran algunos
perfiles calculados a partir de estas ecuaciones.
b) Demostrar que, en el lmite cuando B tiende a infinito, la solucin anterior concuerda
con la de las ecuaciones 10.5-21 a 10.5-23.
c) Efectuar comparaciones numricas de los resultados en la ecuacin 10.5-22 y la figura
Z 0.0,0.5,0.9 1.0
N 2
10.B-18 para
en
y
.
d) Suponer que la ecuacin 9.6-9 es vlida. Demostrar que los resultados en la figura 10.B18 corresponden a un lecho catalizador de longitud L igual al dimetro de 4 partculas.

L / Dp
Debido a que en reactores industriales la relacin
rara vez es menor que 100, se
concluye que una suposicin razonable en clculos de diseo en estado estacionario es

ef , zz
despreciar
Solucin:
a) Las ecuaciones diferenciales en 10.5-6 a 10.5-8, y las condiciones lmite de 10.5-9 a
10.5-14 son todava vlidas, pero con forma lineal para la funcin F:
T T0
Sc1,
T1 T0

S c S c1

Sc1 T0
En el cual
y son constantes que describen la dependencia lineal de la velocidad de
reaccin de la temperatura.
I
2 I
dT d T

C
p
0
z
,ef,
z 0
10.5 6
dz
dz 2
Zona I:
II
2 II
dT d T S F

C
p
0
z
,ef,
c1
0 z L
10.5 7
dz
dz 2
Zona II:
III
2 III
dT d T

C
p
0
z
,ef,
z L
10.5 8
dz
dz 2
Zona III
Aqu hemos supuesto que podemos usar el mismo valor de la conductividad trmica
efectiva en las tres zonas. Estas tres ecuaciones diferenciales de segundo orden estn
sujetas a las seis condiciones lmite siguientes:
10.5 9
T I T1
z
C.L.1:
en
,

10.5 10

C.L.2:

en

z0

T I T II

 10.5 11
10.5 12
10.5 13
10.5 14

C.L.3:

z0

en

C.L.4:

C.L.5

zL

en

en

C.L.6:

zL
en

z ,ef,
,

T II T III

z ,ef,
,
z

dT I
dT II
z ,ef,
dz
dz

dT II
dT III
z ,ef,
dz
dz

T III finita

Las ecuaciones 10.5-10 a 10.5-13 expresan la continuidad de la temperatura y la

densidad de flujo de calor en los lmites entre las zonas. Las ecuaciones 10.5-9 y 10.514 especifican los requisitos en los dos extremos del sistema.
La solucin se facilita introduciendo las siguientes variables adimensionales:
10.5 15
Z z/L

10.5 16
10.5 17
10.5 18

T T0
T1 T0

C p0 L
z ,ef,

SC1 L

C p0 T1 T0

Ahora procedemos a adimensionalizar las ecuaciones diferenciales y las condiciones

lmite:
z
Z
dz 2 L2 dZ 2
dz LdZ

L
De:
T T0

T1 T0
T1 T0 d dT
T1 T0 d 2 d 2T

De:
i) Reemplazando equivalencias en la ecuacin diferencial 10.5-6:
dT I
d 2T I
C p0
z ,ef,
dz
dz 2

C p0 T1 T0 d I
T1 T0 d 2 I
z ,ef,
L
dZ
L2
dZ 2
L / C p0
Multiplicando la ecuacin anterior por
resulta:
z ,ef, d 2 I 1 d 2 I
d I

2
dZ
B dZ 2
C p0 L dZ
10.5 19
ii) Reemplazando equivalencias en la ecuacin diferencial 10.5-7:
dT II
d 2T II
C p0
z ,ef,
S c1 F
dz
dz 2

 C p0 T1 T0 d II
T T d 2 II S T T0
z ,ef, 1 2 0

c1
L
dZ
L
dZ 2
T1 T0
L / C p0 T1 T0

Multiplicando la ecuacin anterior por

resulta:
d II z ,ef, T1 T0
L
d 2 II

2
2
dZ
L
C p0 T1 T0 dZ
T T0
L

T

T
1 0 C p0 T1 T0

Sc1
d II 1 d 2 II

N II
dZ
B dZ 2

10.5 20

iii) Reemplazando equivalencias en la ecuacin diferencial 10.5-8:

dT III
d 2T III
C p0
z ,ef,
dz
dz 2

C p0 T1 T0 d III
T1 T0 d 2 III
z ,ef,
L
dZ
L2
dZ 2
L / C p0
Multiplicando la ecuacin anterior por
resulta:
III
2 III

d
d
1 d 2 III
z ,ef,

2
dZ
B dZ 2
C p0 L dZ
10.5 21

10.5 22
10.5 23
10.5 24
10.5 25
10.5 26
10.5 27

C.L.1:

en
Z

C.L.2:

en
Z

C.L.3:
C.L.4:

en
en

en

0
0
L
0
0
L

L
Z 1
L
Z

C.L.5

L
1
L

,
,

Z
L

I
,

T1 T0
1
T1 T0

I II
d I d II

dZ
dZ

II III
d II d III

dZ
dZ

III finita

C.L.6:
en
,
Las soluciones de las ecuaciones diferenciales homogneas lineales de segundo orden
con coeficientes constantes, son:
1 d 2I d I
d 2I
d I

B
0

B dZ 2
dZ
dZ 2
dZ
Regin I:
Ecuacin auxiliar:

r 2 Br 0

r 2 Br

Br

Por tanto la solucin general es:

I Z C1 C2 e BZ
10.5 28
d II 1 d 2 II

N II
dZ
B dZ 2

Regin II:

1 d 2 II d II

N II 0
B dZ 2
dZ
d 2 II
d II

B
BN II 0
2
dZ
dZ

m 2 am b 0;
Ecuacin auxiliar:

con a = B, b = BN
a a 4 1 b
2

m
m

B B 2 4 1 BN
2
B B 1 4N / B
2

B
1 1 4N / B
2

Luego la solucin general es:

II Z C3em3 Z C4 e m4 Z ,
10.5 29
Regin III:

10.5 30
10.5 31
10.5 32

para m m

Anloga a la regin I.
III Z C5 C6 e BZ

B
1 1 4N / B
2

B
1 1 4N / B
2

en las que:

Aplicando la condicin lmite de la ecuacin 10.5-14, se obtiene

C6 0

; por lo tanto, la
C1

temperatura en la zona III es constante. De la ecuacin 9.5-9 se deduce que

es la
unidad. Al aplicar las cuatro condiciones lmite restante se obtienen los siguientes
1 4N / B 0
resultados para el perfil de temperatura en las tres zonas (cuando
):
m
m
m m e e
e m m Z
I Z 1
2 m
2 m
m
e

m
e

10.5 33

Regin I:

 10.5 34

m e m e m Z m e m em Z
m m
m2 e m m2 e m

II Z
Regin II

10.5 35

m2 m2 m m
e
2 m
2 m
m e m e

III Z 1

Regin III:
b) En preparacin para tomar el lmite de B va al infinito, en primer lugar los siguientes
desarrollos en serie de Taylor
1
1
m (1 (1 (4 N / B) ...) B N O( B 1 )
2
2
1
1
m (1 (1 (4 N / B) ...) N O( B 1 )
2
2

N
B
Importante es sealar que el valor lmite de
para infinito no es cero, pero
.Entonces, las expresiones limitantes para los perfiles de temperatura en las tres
regiones son

m m e m ... ( m m ) Z
lim 1 ( Z ) lim 1 2 m
1 ( desde Z es negativo)
e
B
B
m e ...

m (m ...)em em Z m (m ...)e m e m Z
NZ

e
2 m
2 m
B
m
e

...
m
e

..

lim 1 ( Z ) lim

m2 ... m m
eN
e
B m 2 e m ...

lim 11 ( Z ) lim

estos resultados son consistentes con Eqs.10.5-21,22 y 25. Para obtener el segundo
forma de ecuacin Ec. 10,5-22, tenemos que substituye

F ( )

e integrar.

10C.l Calentamiento de un alambre elctrico cuyas conductividades elctrica y trmica

dependen de la temperatura. Encontrar la distribucin de temperatura en un alambre que se
caliente elctricamente, cuando las conductividades elctrica y trmica varan con la
temperatura corno sigue:
k
1 1 2 2 ...
10C.1 1 k0

ke
1 1 2 2 ...
10C.1 2 ke 0
Aqu

k0

ke 0

son los valores de las conductividades a la temperatura

T0

,y

(T T0 ) / T0

es una elevacin adimensional de la temperatura. Los coeficientes

y
son constantes.
Estos desarrollos en serie son tiles sobre intervalos de temperatura moderados.
a Debido al gradiente de temperatura en el alambre, la conductividad elctrica es una funcin
ke (r )
de la posicin.
, Por tanto, la densidad de la corriente tambin es una funcin de

r : I (r ) ke (r ) ( E / L)

Se ( r ) ke ( r ) ( E / L)

, y la fuente de calor el6ctrico tambin depende de la posicin:

. Entonces, la ecuacin para la distribucin de temperatura es

10C.1 3

1 d
dT

rk ( r ) ke ( r )
r dr
dr

Ahora, introducir las cantidades adimensionales

demostrar entonces que la ecuacin 10C.l-3 se vuelve

10C.1 4

r/R

B ke 0 R 2 E 2 / ke 0 L2T0

k
1 d k d
B e

d k0 d
ke 0

Cuando en esta ecuacin se insertan las expresiones en desarrollo en serie de potencias para
las conductividades, se obtiene
1 d
d

(1 1 2 2 ...) B(1 1 2 2 ...)

d
10C.1 5 d

sta es la ecuacin que debe resolverse para la distribucin adimensional de temperatura.

i
i
Empezar observando que s todas las
y las
fuesen cero (es decir, que ambas
conductividades fuesen constantes), entonces la ecuacin 10C.1-5 se simplificara a
1 d d

B
10C.1

d d

Cuando esta ecuacin se resuelve con las condiciones lmite de que

1
0
en
, se obtiene:
1
2
10C.1 7 4 B 1

finita

en

,y

sta es la ecuacin 10.2-13 en notacin adimensional.

Ntese que la ecuacin 10C.1-5 tendr la solucin de la ecuacin 10C.1-7 para valores

pequeos de
; es decir, para fuentes de calor dbiles. Para fuentes de calor ms
poderosas, postular que la distribucin de temperatura puede expresarse como una serie de
potencias en la intensidad adimensional

10C.1 8

Aqu las

de la fuente de calor:

1
B 1 2 1 B1 B2 ...
4

, son funciones de

pero no de

. Sustituir la ecuacin 10C.1-8 en la

B
ecuacin 10C.1-5 e igualar los coeficientes de potencias semejantes de
para obtener un
n
n 1, 2,3
conjunto de ecuaciones diferenciales ordinarias para las
, con
Estas

ecuaciones pueden resolverse con las condiciones lmite de que

n 0
1
en
De esta manera, obtener

10C.1 9

1
1
1

B 1 2 1 B 1 1 2 1 3 2 O( B 2 )
4
16
8

O( B 2 )
c

n finita 0

B2

donde
significa "trminos del orden de
y superiores"
Para los materiales descritos por la ley de Wiedemann-Franz-Lorenz (vase S9.5) la razn
k / k eT
es una constante (independiente de la temperatura). Por tanto,
k
k
0
10C.1 10 keT ke0T0

10C.1 11

Combinar esta ecuacin con las ecuaciones 10C.1-1 y 10C.1-2 para obtener
(1 1 2 2 ...) (1 1 2 2 ...)(1 )

Igualar los coeficientes de potencias semejantes de la temperatura adimensional para

i
i 1 1 1 2 1 2
obtener relaciones entre las
y las
:
,
, etc . Usar estas
relaciones para obtener

10C.1 12

1
1

B 1 2 1 B 1 2 1 2 2 ) O( B 2 )
4
16

Solucin:
Tenemos las cantidades adimensionales:
r R
dr Rd

De:
T T0
d 1 dT
dT
d

T0 .
T0
dr
T
dr
dr
Rd

0
De:
Ahora reemplazando estas cantidades en la ecuacin diferencial 10C.1-3, tenemos:

1 d
dT
E
rk (r ) ke (r )
r dr
dr
L

1 d

R Rd

k0T0
R2

k
ke 0 E
d
R k (r ) 0 T0
ke ( r )

k0
ke 0 L
R d

ke E 2
1 d k d

e0

d k0 d
k0 L

1 d k d
R 2 E ke
ke 0

d k 0 d
k 0T0 L k e 0

La anterior ecuacin se transforma fcilmente en 10C.1-4, teniendo en cuenta que

B ke 0 R 2 E 2 k0T0 L2 .

k
1 d k d
B e

d k0 d
ke 0

El uso de las ecuaciones. 10C.1.-1 y 2 lleva directamente a la ecuacin. 10C.1-5

1 d
d
2
2

(1 1 2 ...) B(1 1 2 ...)

d
d
b. Cuando la ecuacin 10C.1-8 se sustituye en la ecuacin. 10C.1-5 obtenemos
14 B(1 2 )(1 B1 B 2 2 ...)
En esta nueva ecuacin:
1 d
1

1 1 B 1 1 B1 ...
d
4

2 B 1 2 1 B1 ...
4

d
1

...
B 1 B1 ..

d 4

1
1

B 1 1 B 1 2 1 B1 ... 2
B 1 2 1 B1 ...

4
4

...

B1

Igualando trminos que contienen

, obtenemos
1 d
d 1

(1 2 ) B

d d 4
1 4 4 4 44 2 4 4 4 4 43
B1

1 d

d 1
1 d
d 1
1

(1 2 )1
(1 2 ) B 2 ...

1 4
4
4
d
d
1 4 4d 4 4 d4 4 4 4 4 4 4 44 2 4
4 4 4 4 4 4 4 4 4 4 4 43
1
B
{ 1 4 (1 ) B ...
1
4
4
4 2 4 4 43
1
B
2

B2

B2

1 d
d 1
B 1 2 B

d d 4

Lo cual, cuando se llevan a cabo las diferenciaciones, dar una identidad B = B, como debe ser.
B2
Luego equiparamos los trminos que contienen
, lo que da:

1 d d 1
1 d
d 1

B 1 2 B1
1 14 B 1 2
B 1 2
4
4

d d
d
d

1 B. 14 B 1 2

Divisin por

1
4

B2

y la realizacin de la diferenciacin en el trmino

luego da las
1
ecuaciones diferenciales para
multiplicando por (-1)

1 d d
1 d
d 1

1 2 1
1 1 2
1 2 1 1 2

d d
d
d

1 d d
1

1 2 1 1 1 2

d d

d
d
1 1 2

d
d

1
4

1
2

1 d d
1 d

1 2 1 1 1 2
1 1 2

d d
d
2

2 4

1 d d
1 d
2
2

1
d 1 2
1
d d

1
1 d d

1 2 1 1 1 2 1 2 4 3

d d
2

1 d d

1 2 1 1 1 2 1 1 2 2

d d

1 d d

1 2 1 1 1 2 2 1 1 2
d d

La realizacin de dos integraciones da:

1
1
1 2 1 1 2 2 4 1 4 2 4 C1 ln C2
8
16

C1
La constante de
debe ser igual a cero con el fin de satisfacer la condicin de contorno en el
eje del tubo, y la aplicacin de la condicin de contorno en la pared del tubo da:
1
C.L.1
1
3
0 1 1 C2
C2 18 1 163 1

8
16

1
2

Tomando la diferencia entre las dos ltimas ecuaciones da (despus dividiendo por
1 2 1 18 1 2 2 4 161 1 4 2 4 81 1 163 1

1
2

1
8

18 1 4 2 2 1 161 1 4 4 2 3

1 1 2 1 2
1 2

1
16

1 1 2 3 2
1 2

1
1
1 1 2 1 3 2
8
16
De acuerdo con la Ecuacin 10C.1-9
c. Entonces reordenando la ecuacion10.C-10 tenemos:
T T0 T

1
T0
T0

T
1
T0

k T
k
k
e
e 1
k0 ke 0 T0 ke 0
En esta ecuacin lo sustituimos por las ecuaciones 10C.1-1 y 10C.1-2, para obtener

1
1

... 1 1 2 2 ... 1

1 1 2 2 ... 12 2 3 ...
1 1 1 2 1 2 ...

1 1 1 2 2 1

Igualando los coeficientes de las potencias iguales de

nos da:
,
y as
sucesivamente. Entonces la ecuacin 10C.1-12 sigue directamente:
Partiendo de la ecuacin 10C.1-8:
14 B(1 2 )(1 B1 B 2 2 ...)
Y teniendo en cuenta que:
1 18 1 1 1 2 161 1 3 2
22 18 1 1 1 2 161 1 3 2

2 1 1 1 2 1 3 2
16
2 1 1 2 1 2 31 2
16

16

21 2 2 2 2 31 1 2
21 2 2 2 2 31 1 2
16

1 2 2 2 1 2
16
1 2 2 1 2
16

161 1 2 1 2 2
Luego multiplicando por B se tiene:

B1 161 B 1 2 1 2 2

1
1

B 1 2 1 B 1 2 1 2 2 ...
4
16

14 B(1 2 )[1 161 B( 1 2) ( 1 2) 2 ) O( B 2 )]

10C.2 Calentamiento viscoso con viscosidad y conductividad trmicas que dependen de la

temperatura (figuras 10.4-1 y 10.4-2). Considrese la situacin de flujo que se muestra en la
figura 10.42. Tanto la superficie estacionaria como la superficie mvil se mantienen a una
T0
p
k
temperatura constante
. Entonces, las dependencias de
y
respecto a la temperatura
estdn dadas por
k
1 1 2 2 ...
10C.2

1
k

0
ke

1 1 2 2 ...
10C.2

ke 0 0

Donde las

ai

"evaluado en

y las
T T0

1/
son constantes,

es la fluidez, y el subndice "0" significa

(T T0 ) / T0
. La temperatura adimensional se define como

a) Demostrar que las ecuaciones diferenciales que describen el flujo viscoso y la

conduccin de calor pueden escribirse en las formas

d d

0
10C.2 3 d 0 d
2

d k d
d

Br
0
0 d
10C.2 4 d k0 d

Donde

z / b x / b
,

Br 0b2 / k0T0

(el nmero de Brinkman)

b) La ecuacin para la distribucin adimensional de temperatura puede integrarse una vez

d / d C1.( / 0 )
C1
a fin de obtener
, donde
es una constante de integracin.
Luego, esta expresin se sustituye en la ecuacin de energa para obtener
d
d
2
2
2
1 1 2 ... BrC1 1 1 2 ... 0
10C.2

5
d

10C.2 6

Primero obtener los dos primeros trminos de una solucin en la forma:

10C.2 7

Br Br1 Br 2 2 ...

Br 0 Br1 Br 22 ...
Adems, se sugiere que la constante de integracin C1 tambin se desarrolle como una
serie de potencias en el nmero de Brinkman, para llegar a

10C.2 8

C1 Br C10 BrC11 Br 2C12 ...

y b
c) Repetir el problema, cambiando la condicin lmite en
especificar la temperatura).

qx 0

(en vez de

Respuestas:
1
Br 1 3 2 2 3 ...
12
b)
1
1
1
Br 2 Br 21 2 2 3 4
Br 2 1 2 2 3 4 ...
2
8
24

c)

1
Br 1 2 3 2 3 ...
6
1
1
1

Br 2 Br 21 4 2 4 3 4
Br 2 1 8 2 2 4 3 4 ...
2
8
24

Solucin:
a) El balance de energa y de impulso de momento est dado por:

WLex x WLex

x x

Dividiendo ambos miembros por

obtiene:

ex x ex

lim

x 0

WLx

y tomando el lmite cuando

x x

dex
0
dx

tiende a cero, se

d d z

0
dx
dx

d dT
d z
k
0
dx dx
dx

La primera ecuacin se multiplic por

b 2 / 0b

b 2 d d z
0
0b dx
dx
Donde utilizamos las siguientes variables adimensionales

z
b
b z
b d d z

x
b
b x
dx bd

d z b d

dx
bd

Reemplazo

b 2 d b d

0
0b b d b d

d
d

0
0 d

La segunda ecuacin se multiplica por

b 2 / 0b2

d dT
d z
k
0
dx dx
dx

Las variables adimensionales son:

z
b

b d d z

Br
;

0b 2
k0T0

dx b d
2

b 2 d dT
b 2 d z
k

0
0b 2 dx dx
0b 2 dx

Reemplazo

b2
d T0 d
b 2 b d
k

0
0b 2 bd b d
0 b 2 b d
2
2
b2
d T0 d
b 2 b d
k

0
0b 2 b d b d
0 b 2 b 2 d

b 2 T0 d
0b 2 b 2 d

Multiplicando por

d
d
k d d 0

0b 2

y dividiendo por

T0

k0

la ecuacin anterior queda:

k d
0b 2 d

0
0 K 0T0 d
k 0 d

d
d

d k d
0b 2 d

0
d k 0 d
0 k 0T0 d
2

k d
d

Br
0
0 d
k0 d

d
d

b. Cuando Eqs.10C.2-5,6,7 y 8 se utilizan en la ecuacin de la temperatura, y los coeficientes de

potencia iguales a Br. Se igualan para obtener los trminos que contienen la primera potencia de

d 2 1
C102 0
2
d
Resolviendo la ecuacin diferencial nos da:
d d
2

C 10
d d
d 1
2
C 10 d
d

d 1
C102 C3
d

C102 d C3 d C4

1
1 C102 2 C3 C4
2

1 0
Las constantes de integracin se determinan a partir de las condiciones de contorno que
a

0,1
1
0 C102 (0)2 C3 (0) C4
2
C4 0

1
1 C102 C3
2
1
0 C102 (1) C3
2
1
C3 C102
2

Reemplazando las constantes

1
1
1 C102 2 C102 (0)
2
2

De este modo se obtiene el siguiente resultado:

1
1 C102 2
2

Pasamos ahora a la variable

d d

0
d 0 d
d
0
0 d

d
C1
0 d

(junto con la Ecuacin 10C.2-2), que es:

C1
d
0

1 1 2 2 ...
0 0
d
C1 (1 1 2 2 ...)
d

Las expresiones de las ecuaciones 10C.2-6, 7 y 8

10C.2 7 10C.2 6

Br Br 1 Br 2 2 ...

10C.2 8 Br 0 Br1 Br 22 ...

C1 Br C10 BrC11 Br 2C12 ...
Reemplazamos:
d
0 Br 2 Br 2 2 ... C10 Br C11 Br 2 C12 ...

1 1 Br 1 Br 2 2 ...

2 Br 1 Br 2 2 ... ...
2

Igualamos Br=0
d0
C10 1 1 (0) 2 (0) ...
d
d0
C10
d
d0
C10
d

d C
0

10

0 C10 C2

as condiciones de contorno son:

Aplicando l

0 (0) 0 0 (1) 1
y

0 C10 (0) C2
C2 0

0 C10
1 C10
Y obtenemos:

0 (1) (0)

0
Por lo tanto tambin en:

1
1 C102 ( 2 )
2
1

1 2
1 2
2

1
2
2

A continuacin volvemos a la ecuacin de la energa y equiparamos los trminos

d 1 d
d
2 C11 1 1 0
1 1
d
d d

d 1 d 1
d 2 1 d 2 2
1 1

2C11 1 1 0
d d
d 2
d 2

d 2 2
d 2 d 2
d 1
1
1
2 C11 1 1
d 2
d

d 2

1
( 2 )
2

Como
d 1 1
1 2
d 2

Br 2

d 2 1
(2) 1
d 2 2
Reemplazamos:

d 2 2

1
1
1 1 2 6 3 4 2C11 1 2
2
d
4
2
2
d 22

1
1 1 6 6 2 2C11 1 2
2
d
4
2
d d 2
1

1 1 6 6
d d
4

2C

11

1
1 ( 2 )
2

d 2
1
1
2
2
1 1 6 6 d 2C11 d 1 d
d
4
2

d 2
1
1
1
1
1 3 2 2 3 d 2C11 d B1 2 3 d C3 d
d
4
2
3
2
2

1 1 2 3 1 4
1 1 3 1 4
2
C11 1 C3 C4
4 2
2
2 6
12

2 0 y 0;1.
Las condiciones de contorno son
10D.1 Prdida de calor desde una aleta circular (figura 10D.1).
a) Obtener el perfil de temperatura T(r) para una aleta circular de espesor 2B sobre un tubo
cuya temperatura en la pared exterior es T 0. Hacer las mismas suposiciones que se hicieron
en el estudio de la aleta rectangular en 10.7.
b) Deducir una expresin para la prdida de calor total desde la aleta.
Solucin:

r
a) El balance de energa se hace sobre un anillo de espesor
de la aleta circular. Se
encuentra en estado estacionario y la nica contribucin a la densidad de flujo de energa es
q. Por tanto, el balance de energa es:
2 r 2 B qr r 2 r 2 B qr r r 2h 2 r r T Ta 0
2 2B r

Dividiendo ambos miembros por

y tomando el lmite cuando
se obtiene:
rqr r qr r r
rh T Ta
lim

r 0
r
B

10 D.1 1

d rqr
dr

tiende a cero,

rh
T Ta
B
qr kdT / dr

Ahora insertamos la ley de Fourier (

del metal. Si k es constante, entonces se obtiene:
d dT
rh
r
T Ta

10 D.1 2
dr dr
Bk

), donde k es la conductividad trmica

Esta ecuacin debe resolverse con las condiciones lmite:

10 D.1 3
r R0 T T0
C.L.1:
en
,
dT
r R1 dr 0
10 D.1 4
C.L.2:
en
,
Ahora introducimos las siguientes cantidades adimensionales:
T Ta

T0 Ta
10 D.1 5
= temperatura adimensional
r

R0
10 D.1 6
= distancia adimensional
y procedemos a adimensionalizar la ecuacin diferencial 10D.1-2.

T Ta
T0 Ta

r
R0

De:

T0 Ta d dT ;

R0 d dr

De:
Reemplazando equivalencias en la ecuacin diferencial ordinaria 10D.1-2:
T T d R0 h T T
d
R0 0 a
0 a
R 0 d
Bk
R 0 d

R0 h
1 d T0 Ta d
T0 Ta


R0 d
d
Bk

T0 Ta

R0 h
1 d d
T0 Ta


R0 d d
Bk

R0

Multiplicando ambos miembros por

y dividiendo por , tenemos:
hR 2
1 d d

d d
Bk
10 D.1 7
Esta ecuacin tiene la siguiente solucin (con
1 d d
2


d d

2 hR0 2 / Bk

):

1 d d
d 2
2

d d
d 2

1 d d 2

2
d d 2
Esta ecuacin diferencial se puede acomodar para tener una ecuacin modificada de
Bessel,
Ecuacin
Solucin
2
2 2
2
y C1 I n (ax) C2 K n (ax)
x y '' xy ' (a x n ) y 0; n 0

d 2
d

2 2 0; n=0
2
d
d

( ) C1 I0 ( ) C2 K 0 ( )

Luego, las constantes son determinadas con las condiciones lmite:

1
1
C.L.1:
en
,
R1 / R0
d / d 0
C.L.2:
en
,
Aplicando la primera condicin lmite:
10 D.1 8 1 C1 I0 ( ) C2 K 0 ( )
Aplicando la segunda condicin lmite:

 10 D.1 9

0 C1 I1 ( R1 / R0 ) C2 K1 ( R1 / R0 )

Estas dos ecuaciones pueden ser resueltas simultneamente

Despejando C2 de la ecuacin 10D.1-9 nos da:
C I ( R1 / R0 )
C2 1 1
K1 ( R1 / R0 )
Reemplazando en la ecuacin 10D.1-8 tenemos:
C I ( R1 / R0 )
1 C1 I0 ( ) 1 1
K 0 ( )
K1 ( R1 / R0 )

K1 ( R1 / R0 ) C1 I 0 ( ) K1 ( R1 / R0 ) C1 I1 ( R1 / R0 ) K 0 ( )
K1 ( R1 / R0 ) C1 I0 ( ) K1 ( R1 / R0 ) I1 ( R1 / R0 ) K 0 ( )
C1

K1 ( R1 / R0 )
I0 ( ) K1 ( R1 / R0 ) I1 ( R1 / R0 ) K 0 ( )

Asimismo se procede para hallar la constante C2:

I1 ( R1 / R0 )
C2
I0 ( ) K1 ( R1 / R0 ) I1 ( R1 / R0 ) K 0 ( )
Por lo tanto el perfil de temperatura adimensional es:

( )

K1 ( R1 / R0 )
I ( )
I0 ( ) K1 ( R1 / R0 ) I1 ( R1 / R0 ) K 0 ( ) 0

( )

I1 ( R1 / R0 )
K ( )
I0 ( ) K1 ( R1 / R0 ) I1 ( R1 / R0 ) K 0 ( ) 0

K1 ( R1 / R0 ) I 0 ( ) I1 ( R1 / R0 ) K 0 ( )
I0 ( ) K1 ( R1 / R0 ) I1 ( R1 / R0 ) K 0 ( )

b) La prdida de calor total desde la aleta es:

dT
d
Q 2 R0 2 B k
4 Bk T0 Ta
dr r R0
d 1
I1 ( R1 / R0 ) K1 ( ) I1 K1 ( R1 / R0 )

4 Bk T0 Ta

I1 ( R1 / R0 ) K 0 ( ) I 0 ( ) K1 ( R1 / R0 )

UNIVERSIDAD NACIONAL JORGE BASADRE GROHMANN

FACULTAD DE INGENIERA
ESCUELA ACADMICO PROFESIONAL DE INGENIERA QUMICA

Tema: CAP. 1O: BALANDE DE ENERGIA EN LA ENVOLTURA Y

DISTRIBUCION DE TEMPERATURA EN SOLIDO Y EN FLUJO LAMINAR

Asignatura: Fenmenos de Transporte II

Docente: Mgr. Pedro Cornejo del Carpio
Grupo: 1
Integrantes:
-

Rossi Elizabeth Achata Merlin

Viviana Ximena Rodriguez Almanza
Carlos Roberto Costa Gil

3er ao

Tacna-Per
2015

2013-39329
2013-39343
2013-39321