CAPITULO I

CONDUCCION

INTRODUCCION

El flujo de calor conducido a través de un cuerpo es en realidad un flujo a

través de sus tres dimensiones y solamente en ciertos casos puede ser considerado
como un flujo Bidimensional o unidimensional.

El estado térmico de los cuerpos debido al flujo de calor intercambiado con el

medio externo, también puede ser influenciado por fuentes internas de calor; estos
efectos combinados pueden someter al cuerpo a un estado “estable” o “inestable de
flujo de calor”, que nos indica que el cuerpo se encuentra en un estado de energía
interna constante o variable, en este último caso el cuerpo estará calentándose o
enfriándose.

Estos efectos pueden ser expresados matemáticamente, mediante una

ecuación, denominada la “Ecuación General de conducción”, la cual podemos
expresarla en los diferentes sistemas coordenados; el estudio del Flujo de calor y la
distribución de temperaturas pueden ser hallados a través de estas ecuaciones, que
a su vez pueden particularizarse desde el estudio tridimensional al unidimensional,
des el estado inestable al estable de flujo de calor, en sistemas con o sin generación
interna y en las posibles combinaciones que puedan darse.

1 LA ECUACION GENERAL DE CONDUCCION

La ecuación general de conducción deberá establecerse considerando las tres

dimensiones espaciales y un punto del espacio podrá definirse usando
cualquiera de los tres tipos de sistemas de coordenadas que conocemos.

1.0.1 Descripción de los sistemas coordenados

- Sistemas de coordenadas cartesianas. – En este sistema un

Z punto espacial definido por las tres coordenadas lineales X, Y,
Z

P(X,Y,Z)

Y
 - Sistemas de coordenadas cilíndricas. – En este sistema un
punto queda definido por una coordenada angular , una
coordenada lineal Z y una radial r

Y ( r ,  , Z )

X


- Sistemas de coordenadas esféricas.-

En este sistema un punto queda definido

P(г,Φ,Z) por una coordenada radial r y dos
coordenadas angulares (  ,Y)

г
 1.0.2 Balance de calor

El balance energético lo realizamos considerando un elemento

diferencial dentro de un cuerpo sólido

- sistemas de coordenadas cartesianas

Calor que ingresa por Calor que sale Calor generado por Cambio de
conducción en las en las tres fuentes internas en energía
tres direcciones + direcciones + el = tiempo (de) interna en el
X,Y,Z por unidad de X,Y,Z en el tiempo de
tiempo (de) tiempo de

Matemáticamente:

. . .
( qx  q y  qz ).d  q m .dv.d 

. . .
(q ( x  dx)  q(Y  dY )  q( Z  dz) ).d  .dv.C.dt
 Donde:

. . .
qx ' qY ' Z = flujos de calor por unidad de tiempo que ingresan en las
direcciones X,Y,Z al sólido diferencial en (WATT).

. . . . .
.q x ' q x ' q Z q(Y  dY )' q( Z  dz)
Y

= Flujos de calor por unidad de tiempo que salen en las

direcciones (X,Y,Z) del sólido diferencial en (WaTT)

d = Lapso de tiempo en el cual se produce el análisis en (seg)

q ''' = Calor generado por unidad de volumen en el cuerpo diferencial en

(W/m3)

dV = Volumen del elemento diferencial en (m3)

 = Densidad del sólido diferencial en (kg/m3)

C = Calor específico del sólido en (J/kg.k)

dt = Variación de la temperatura durante el tiempo d en el sólido en

(K o °C)

Si la función en (X+dX) puede expresarse en (X) según la serie de

Taylor o el teorema del valor medio; de la siguiente manera:

q X  2qx
q( X  dX ) = q ( X ) + .dX  .dX 2  ...
X x 2

Utilizando solo 2 términos extensivo el análisis a las otras direcciones Y

y Z, se tiene:

q X
q( X  dX ) q x  .dx
x

qY
q(Y  dY )  qY  .dY
Z

q z
q( Z  dZ )  q Z  .dZ
Z

Reemplazando en la ecuación general de conducción esta se convierte en:

 q X q q
( .dX  Y .dY  Z .dZ ).d  q ''' .dv.d
X Y Z

- Sistema de coordenadas cilíndricas

En este caso el sólido diferencial puede expresarse de la siguiente manera:

 BALANCE DE CALOR

Cambio de
Calor que Calor generado Calor que sale energía
ingresa por por unidad de por conducción + interna en el
conducción en + volumen en el = en las tiempo (d )
las direcciones tiempo (d ) direcciones r, Z,
r, Z,  en el  en el tiempo
tiempo (d  ) (d )

MATEMATICAMENTE
. . . . .
(q r  q Z  q ) . d  + q .dv.d  (q r  dr  qZ  dZ  q  dZ  q  d ).d
. . . '''

 .dv.c.d 7
 y como en el caso anterior con el teorema del valor medio se obtiene:

q r q q
( .dr  Z .dZ  .d ), de  q m .dv.d   .dv.c.dT
r Z 

- Sistema de coordenadas esféricas

El sólido diferencial se expresa de la siguiente manera:

Sistemas de coordenadas esféricas

Calor que ingresa Calor Calor que sale Cambio de

por conducción generado por por conducción energía
en las direcciones + unidad de = en las + interna en
r,  y  en el volumen en el direcciones r, el tiempo
tiempo d tiempo (d ) y en el (d )
tiempo (d )

.
. . . .
(q r  q  q ).d  q .dv.d  (q r  dr  q  d  q  d  ).d  .dv.c.dtT
'''
1.1.0 LA CONDUCCION EN EL ESTADO ESTABLE

El estado estable de flujo de calor en todos lo sistemas es momentáneo, sin

embargo el régimen de funcionamiento de muchos sistemas es semejante al
estado estable de flujo de calor en consecuencia su estudio es de vital
importancia el cálculo y diseño térmico de muchos aparatos, así como
también en el conocimiento físico de algunos fenómenos térmicos.

1.1.1.1 ECUACION GENERAL DE CONDUCCION EN EL ESTADO ESTABLE DE

FLUJO DE CALOR

Un sistema en estado estable de fluido de calor, cede la cantidad de calor

que recibe y/o general en consecuencia su energía interna no cambia,
entonces:

.dV .C.dT  0

Y la ecuación general para el estado estable de flujo de calor es:

q X q q
( .dX  Y .dY  Z .dZ ).d  q ''' .dV .d
X Y Z

 q ''' .dV .d  0

De la misma manera sucede usando las otras coordenadas y los resultados se

muestran en la tabla (1.1.1)

TIPO DE ECUACION GENERAL DE CONDUCCION EN FLUJO

COORDENADAS ESTABLE DE CALOR

CARTESIANAS q X q q
( .dX  Y .dY  Z .dZ )d  q ''' .dv.d  0
X Y Z

CILINDRICAS q r q q
( .dr  Z .dZ  .d ).d  q ''' .dv.d  0
r Z 

ESFERICAS q r q q

( .dr  .d  .d ).d  q ''' .dv.d  0
r  

TABLA 1.1.1
1.1.2 CONDUCCION UNIDIMENSIONAL SIN GENERACION INTERNA

Como los sistemas de conducción unidimensionales mas sencillos no incluyen

la generación interna iniciaremos nuestro estudio con estos sistemas. Si el
flujo de calor es sólo unidimensional (coordenadas X), se tiene:

q y q Z
.dy  .dZ  0
y Z

Y si no existen fuentes internas de calor:

q' ' ' .dv.d  0

Y la ecuación de conducción unidimensional sin generación interna será:

q X
.dX .d  0
X

De igual manera sucederá en cualquier otra dirección, (Y ó Z) en cualquier

sistema coordenada, los resultados los apreciamos en las tabla 1.1.2
 ECUACIONES DE CONDUCCION UNIDIMENSIONAL SIN
GENERACION INTERNA

ECUACIONES
SISTEMAS DE DIRECCION DE FLUJO UNIDIMENSIONALES SIN
COORDENADAS DE CALOR GENERACION INTERNA

X q X
 .dX .d  0
X

CARTESIANAS Y qY
 .dY .d  0
Y

Z q Z
 .dZ .d  0
Z

r q r
 .dr.d  0
r

CILINDRICAS z q Z
 .dZ .d  0
Z

q
  d . .d  0


ESFERICAS q r
r  dr.d  0
r
q
  d .d  0


q
  d .d  0


Tabla 1.1.2
Para la aplicación de estas ecuaciones con la finalidad de hallar la distribución
de temperaturas deberá considerarse que la ecuación del flujo de calor
depende de la conductividad térmica y el área de transferencia de calor. Estos
parámetros no pueden siempre tomarse como valores constantes, ya que la
conductividad generalmente varía con la temperatura y el área de transferencia
depende de la geometría del cuerpo. La Figura. 1.2.1.1 nos muestra la forma
de variación de la conductividad térmica con la temperatura de algunas
sustancias y al final del capítulo presentamos tablas más completas de
conductividades térmicas.
Según lo expuesto aplicaremos la conducción unidimensional a los 3 modelos
básicos (placas planas, cilindros huecos y esferas huecas) y haremos
extensivo nuestros estudio a casos especiales que se presentan, daremos a
conocer también aplicaciones importantes de estos modelos en la ingeniería

Temperatura absoluta (K)

Variación en la conductividad térmica con la temperatura para distintos
materiales
1.1.2.1 CONDUCCION UNIDIMENSIONAL EN PLACAS PLANAS

Para una placa plana es apropiado el uso de coordenadas cartesianas y su

aplicación es fácil cuando la conductividad térmica (K) y el área (A) son
constantes; sin embargo cuando K y A son variable usaremos el método
diferencial que nos facilita el análisis.

Caso 1. Conductividad y área constantes

Sea una placa plana de espesor L y donde

T1  T2 , con un sistema coordenado en X:

L

T1 T2

X=0 X=L

Ecuación diferencial

La ecuación de flujo de calor es:

q x   K . A.(dT / dX )

La ecuación unidimensional es:

 (q x / x).dx.d   0

Reemplazando el flujo de calor:

 K . A.(dT / dX )
 .dx.d  0
X

Simplificaciones:

K = cte
A = cte
dx=integrador
d = No interviene como variable (Flujo Estable de calor).
Usando las simplificaciones correspondientes:
 ( 2T / x 2 ).k. A.dx.d  0

 T / x   0
2 2

Ecuaciones de distribución de temperaturas

 Integrando la ecuación diferencial

 Primera integración:

(T / X )  C1

 Segunda integración

T = C1. X  C2

 Condiciones de contorno

T=T1 para X = 0

T= T2 Para X = L

 Reemplazando en la segunda integración

T1  C 2
T2  C1 .L  C2  C1 .L  T1

 De las ecuaciones anteriores las constantes de integración son:

C1  (T2  T1 ) / L
C2  T1

 Reemplazando las constantes en la segunda integración:

T  (T2  T1 ).( X / L)  T1

 Ecuación para el flujo de calor

 El flujo de calor es:

qk  k. A.(T1  T2 ) / L
 CONCLUSIONES
 La distribución de temperaturas es lineal
 q k . Es el flujo de calor y su signo (+) nos indica el sentido de
izquierda a derecha.


Caso 2.- Conductividad y/o área variable. El método diferencial como

solución simplificada.

En este caso se presentan varias posibilidades que darán origen a

diferentes ecuaciones, por ejemplo:

 K = CONSTANTE, A = VARIABLE

(A / x).(T / x)  A.( 2T / x 2 )  0

 A= CONSTANTE, K= VARIABLE

(k / k ).(T / x)  k.( 2T / dx 2 )  0

 K=VARIABLE, A=VARIABLE

(k / k ).A  (A / x).k . T / x  ( 2T / x 2 ).k.A  0

Como se aprecia la solución de estas ecuaciones diferenciales es

bastante complicada, por lo tanto es mejor usar método diferencial.

 EL METODO DIFERENCIAL

Ilustraremos el uso del método en una placa plana de espesor “L” con
temperaturas T1  T2 de área constante y conductividad variable lineal:

k  k o .(1  T )

Donde:

k o  conductividad térmica del sólido a 0°C

 = Coeficiente de la dilatación térmica
 L

T1 T2

X=0 X=L

ECUACION DEL FLUJO DE CALOR

 De Fourier

q   K . A.(dT / dX )

 Ordenando convenientemente:

L T2
q  (dX / A)    K .dT
O T1

T2
(q.L) / A    k o .(1  T )dT
T1


 k o . T  (  T 2 ) / 2 T2
T1

 k .(To 2  T1 )  ( / 2).(T22  T12 ) 

 k o .1  ( / 2).(T1  T2 ).(T1  T2 )

 k m .(T1  T2 )

 En consecuencia el flujo de calor es:

q  k m .( A / L).(T1  T2 )

Nota: k m  k o .1  ( / 2).(T1  T2 )

 ** ECUACION DE DISTRIBUCION DE TEMPERATURAS

 De Fourier

X T
q. (dx / A)   k o .(1   .T ).dT
0 T1

 Integrando:

(q.X ) / A  ko .(T  T1 )  ( / 2).(T 2  T12 )

 Reemplazando el flujo de calo ya hallado se obtiene la
distribución de temperaturas, que ordenada
convenientemente es:

T  ( / 2).T 2  T1  ( / 2).T12  (k m / k o ).( X / L).(T1  T2 )

NOTA: Como se observa la distribución de temperaturas es

parabólica.

1.1.2.2 CONDUCCION UNIDIMENSIONAL SIN GENERACION INTERNA EN

CILINDROS HUECOS.

Siguiendo el mismo procedimiento que para la placa plana, se tiene:

CASO 1. CONDUCTIVIDAD CONSTANTE ( El Area es variable)

Para un cilindro hueco con radios r1 yr2 ytemperatu rasT1  T2 usando

ecuaciones en coordenadas cilíndricas:
 ECUACION DIFERENCIAL

○ La ecuación del flujo de calor según Fourier:

q   K . A.(dT / dr )

  K .2. .r.L(dT / dr )

○ La ecuación unidimensional en la dirección radial es:

 (qr / r ).dr.d  0

○ Reemplazando la ecuación de Fourier:

 k.2. .r.L(dT / dr )

 .dr.d  0
r

○ Simplificaciones

k ,2,  , L, d  Constantes

dr = Integrador

○ Usando las simplificaciones:

r.(dT / dr )
0
r

 ECUACION DE DISTRIBUCION DE TEMPERATURAS

○ Integrando la ecuación diferencial

- Primera integración:

(T / r )  C1 / r

- Segunda integración:

T  C1 .Int  C2
○ Condiciones de contorno

T  T1 parar  r1

T  T2 parar  r2

○ Reemplazando en la segunda integración:

T1  C1 .Inr1  C 2

T2  C1 .Inr2  C2

○ Las constantes de integración son:

C1  (T1  T2 ) / In(r1 / r2 )

C2  T1  (T1  T2 ).Inr1 / In(r1 / r2 )

○ Usando las constantes en la segunda integración se tiene:

T  T1  (T1  T2 ).In( / 1 ) / In(r1 / r2 )

 ECUACIONES DE FLUJO DE CALOR

○ De Fourier:

q   K .2. .r.L(dT / dr )

 k.2. .r.L(T1  T2 ) In(r2 / r1 )

○ Reemplazando el valor de C1:

q  2 .K .L.(T1  T2 ) / In(r2 / r1 )

CONCLUSION

 El flujo de calor es positivo

  La distribución de temperaturas es logarítmica

CASO 2.- CONDUCTIVIDAD VARIABLE

En este caso si tenemos un cilindro hueco con radios r1 y r2 ;

temperaturas T1  T2 , y conductividad variable lineal:

K  K o (1   .T )

 ECUACION PARA EL FLUJO DE CALOR

○ Usando el método diferencial:

q   KA.(dT / dr )

  K .2 .r.L.(d .T / dr )

○ Ordenado convenientemente e integrando:

r2 T2
q. (dr / r )  2. .L.  K o .(1  T )dT
r1 T1


q.In(r2 / r1 )  2. .L.K o . (T1  T2 ))( / 2).(T12  T22 
○ En consecuencia:

q  2. .L.k m .(T1  T2 ) / In(r2 / r1 )

  LA ECUACION DE DISTRIBUCION DE TEMPERATURAS

○ Con el mismo procedimiento:

r T
q. (dr / r )  2. .L.  k o .(1  T ).dT
r1 T1


q.In(r / r1 )  2. .L.k o . (T1  T )  ( / 2).(T12  T 2 ) 
○ Reemplazando el flujo de calor hallado y ordenando convenientemente se
tiene:

T  ( / 2).T 2  T1  ( / 2).T12  (k m / k o ).(T1  T2 ).In(r / r1 ) / In(r2 / r1 )

1.1.2.3 CONDUCCION UNIDIMENSIONAL EN ESFERAS HUECAS

Similarmente a los modelos anteriores se realiza el cálculo del flujo de calor y

de distribución de temperaturas y los resultados se muestran en la tabla
1.1.2.3.
 TABLA 1.1.2.3.

CONDUCCION UNIDIMENSIONAL. FLUJO ESTABLE DE CALOR. SIN GENERACION INTERNA

ECUACIONES DE DISTRIBUCION DE TEMPERATURAS, FLUJO DE CALOR Y RESISTENCIA TERMICA
CONDUCTIVIDAD TERMICA CONSTANTE CONDUCTIVIDAD TERMICA VARIABLE, LINEAL
K=KO.(1+  T)
MODELO BASICO K=cte
PLACA PLANA X  
T  T1  .(T1  T2 ) T .T 2  T1  .T12  (k m / k o ).( X / L).(T1  T2 )
L 2 2
K.A K m .A
q .(T1  T2 ) q .(T1  T2 )
L L

R  L /( K . A) R  L /( K m . A)

CILINDRO-HUECO
T1  T2   Ln(r / r1 ) k m
T  T1  .Ln(r / r1 ) T .T 2  T1  .T12  . .(T1  T2 )
Lm(r2 / r1 ) 2 2 Ln(r2 / r1 ) k o

2. .K .L 2. .k m .L
q .(T1  T2 ) q .(T1  T2 )
Ln(r2 / r1 ) Ln(r2 / r1 )

Ln(r2 / r1 ) Ln(r2 / r1 )
R R
2. .K .L 2. .k m .L
ESFERA-HUECA (T1  T2 ) 1 1
T  T1  .(  )
1 1 r1 r   r.(r2  r1 ) k m
(  ) T .T 2  T1  .T12  . .(T1  T2 )
r1 r2 2 2 r2 .(r  r1 ) k o
4. .k .r1 .r2 4. .k m .r1 .r2
q .(T1  T2 ) q .(T1  T2 )
(r2  r1 ) (r2  r1 )
1 1 1  1 1 1 
R .   R .  
4. .k  r1 r2  4. .k m  r1 r2 
1.1.2.4. CONDICION UNIDIMENSIONAL EN GEOMETRIAS DIFERENTES

Cuando se trata de encontrar el flujo de calor por conducción y la

distribución de temperaturas de estos sólidos cuya conductividad
podría ser constante, variable lineal o no lineal, se recomienda el uso
del método diferencial como método rápido y de fácil aplicación,
ilustraremos su uso a través de ejercicios numéricos.

EJERCICIO

El sólido que se muestra tiene su superficie lateral perfectamente

aislada. Calcule el flujo de calor a través del sólido para las siguientes
condiciones en grados centígrados y de metros:

 Conductividad térmica variable: K=50x (1+0.002T+0.003T 2)

 Distribución de temperaturas variable según la ley: T 2 = 1800X
 La temperatura en uno de los extremos es de: - °60C
 Identificación de las condiciones de contorno
 CALCULO DE LA LONGITUD DEL SÓLIDO

X=0 T2 = 1800 (0) T = T1 = 0

T1= 0 T2 = -60°C

X=L T22 =1800 (L) (-60)2 = 1800.L

 CALCULO DEL AREA VARIABLE EN FUNCION DE LA

DIMENSION X

 Del sistema coordenado.-

Ax = h x (El ancho = 1 mt)

 De la semejanza en el triángulo:

(10/2) = (a/X)

 Si

hx = 5+2.a = 5 + 10 . X

AX = 5 + 10.X

 FLUJO DE CALOR A TRAVES DEL SÓLIDO

 De Fourier:

q = - K.A. (dT/dX)

 Ordenando e integrando:

L T2
q. (dX / AX )    k.dT
O T1

 Reemplazando valores:

2 dX 60
q.    50(1  0.002T  0.003T 2 ).dT
0 (5  10 XI ) 0

.Ln(5  10 x2)  Ln(5)  50.(0  60)  0.001(0  3600)  0.001x(0  216,000)

1
q.
10
q x 0.16094 = 50 x 272.4

q = 84,627 Watt

EJERCICIO

El sólido tronco – cónico que se muestra en la figura está aislado en

forma perfecta en toda superficie lateral:

a. Calcule el flujo de calor

b. Determine la distribución de temperaturas

r1  0.5m
r2  0.75m
K  116(1  0.0006T

SOLUCION.-

 RELACIONES GEOMETRICAS

 Utilizando la coordenada Y se tiene

 La semejanza nos indica.

(r2  r1 ) / h  m.(h  Y )
m  (r2  r1 ).(h  Y ) / h

 Si
r  r1  m  r1  (r2  r1 ).(h  Y ) / h

 Derivando:

dr  (r1  r2 ).dY / h

dY  (h.dr ) /( r1  r2 )

 FLUJO DE CALOR A TRAVES DEL SÓLIDO

 De Fourier:

q   K . A.(dT / dY )

r1 T1
q. (dY / A)    k.dT
r2 T2

Si: A   .r 2

dY  h /( r1  r2 ).dr

k  116(1  0.0006T )

Reemplazando, simplificando y ordenando, se tiene:

r1dr 116. T
q. 2
 .(r1  r2 ). (1  0.0006).dT
r2 r h T2

 1 1   116.
q.    
.(r1  r2 ). (T1  T2 )  0.0003(T12  T22 ) 
 r1 r2  h

 116.
q
L

.(r1 xr2 ). (T1  T2 )  0.0003(T12  T22 ) 
 Reemplazando valores:


q  116. .(0.5x0.75). (538  93)  0.0003(5382  932 ) 
q  -136.7 x 529.33

q  -72,357 WATT

 CALCULO DE LA DISTRIBUCION DE TEMPERATURAS

Integrando desde r2 hasta r:

r T
q. (dr / r 2 )  (116 / h). .(r1  r2 ). (1  0.0006).dT
r2 T2


q. (1 / r )  (1 / r2 )  (116 / h). .(r1  r2 ). (T  T2 )  0.0003(T 2  T22 ) 
Utilizando los datos:

q  72,357 WATT; r1  0.5m; r2  0.75m y

r  r1  (r2  r1 ) / h.(h  Y )  0.75  0.25Y

Reemplazando y ordenando convenientemente:

T  0.0003T 2  902.7  (3032.4).Y /9.4  (3.1) xY 

ECUACION DE DISTRIBUCION DE TEMPERATURAS EN

EL SÓLIDO TRONCO-CONICO

1.1.2.4 APLICACIONES DIVERSAS

Las aplicaciones generalmente se dan usando un modelo básico para un

modelo físico determinado y como sabemos son: placas planas, cilindros
huecos y esferas huecas.

- Aplicaciones en placas planas

Cuando las otras dos dimensiones diferentes del espesor

(largo y ancho) son lo suficientemente grandes para no
considerar la bidimensionalidad del flujo de calor en una
pared o placa cualquiera, se usa este modelo básico, por
ejemplo:
 - Paredes simples o compuestas
- Aislamientos térmicos
- Conductores térmicos
En todos estos casos puede ser útil calcular cualquiera de los
parámetros típicos, por ejemplo: espesor de material, flujo de
calor, material (conductividad térmica o temperaturas de
operación.

EJERCICIO

Un horno de fierro con paredes blancas, tiene un espesor L1  0.025m y opera

a un tiempo completo (  8,760hr / año) con una temperatura interna
Tg  300C y una temperatura externa Tamb  26C ; el coeficiente combinado
interior es de h crgas =284 W/m2K y la conductividad del fierro es K1 = 62
W/mK.

El horno será aislado externamente con un aislante de conductividad k 2 =

0.087 W/mk, cuyo costo de instalación es de C = 1, 100 S/m 3, la depreciación
anual (mantenimiento, impuestos), es estimada en C1 = 15%, el costo de
combustible para la operación del horno es C2 = 0.8 x 10-5 S/W.hr, si el
coeficiente pelicular externo es 9W/m2.k, calcular:

a) El espesor económico del aislante (L2)

b) El ahorro neto anual por unidad de área

SOLUCION:
Tamb  26C
h amb  9W / m 2 .k
Tg = 300°C C  1.100 S / m 3
hg  284W / m 2 .k
C1  15%
C 2  0.8 x10 5 S / W .hr
q a''  calor perdido con aislante por unidad de área

qb''  calor perdido sin aislamiento por unidad de área

An  ahorro neto anual por unidad de área

 CALCULO DEL ESPESOR ECONOMICO DEL AISLANTE

- El ahorro neto anual está dado por:

An  (qb''  qa'' ) xxC 2  L2 xCxC1 ... (1)

Cantidad de calor

Ahorrado al poner Depreciación

Aislante por unidad

por área

- El calor perdido con aislante es:

(Tg  T)
q a'' 
1 L L 1
 1  2 
hg K1 K 2 h 

Reemplazando en la ecuación (1)

(Tg  T )
An  (qb''  ) xxC 2  L2 xCxC1
1 L1 L2 1
  
hg K1 K 2 h 

Derivando respecto al espesor aislante e igualando a cero.

(Tg  T ) xxC 2
An  (qb''  CxC1  0
1 L1 L2 1 2
(    ) xK 2
hg K1 K 2 h 
 (Tg  T ) xxC 2 1 L L 1
(  1  2  )
CxC1 xk 2 hg K1 K 2 h 

1/ 2
 k x(Tg  T) xxC 2  1 L 1 
L2   2   K 2 x  1      (2)
 CxC1   hg K1 h  

Derivando el ahorro neto por segunda vez:

d 2 An  2 x(Tg  T ) xxC 2
2
 3
0
dL2 1 L1 1 
   2
 xk 2
 hg K 1 h  

En consecuencia la ecuación 2 nos da un máximo

1/ 2
 0.087 x(300  26) x8,760 x0.8 x10 5   1 0.025 1 
L2     0.087    
 1,100 x0.15   284 62 9

L2  0.0906mts.

 CALCULO DEL AHORRO NETO

300  26
qb''  2,381W / m 2
1 0.025 1
 
284 62 9

(300  26)
q a''  236.93w / m 2
1 0.025 0.0906 1
  
284 62 0.087 9

Luego el ahorro neto es:

An  (2,381  23.9) x8,760 x0.8x105  0.0906 x1,100 x0.15

An  135.372S / m 2 xaño
- Aplicaciones en cilindros huecos y esferas huecas

 EL RADIO CRITICO DE AISLAMIENTO

Es el radio de un material considerado aislante que cuando se

coloca sobre una superficie curva, que puede ser cilíndrica o
esférica aumenta la transferencia de calor hasta un valor máximo
comparado con el caso en que no se colocase este material.
Ocurre generalmente cuando se aíslan tuberías o cables de
diámetros menores de 2 cm. y puede ocurrir si se aíslan superficies
curvas con materiales cuya conductividad térmica sea mayor que
los considerados aislantes.

ANALISIS MATEMATICO PARA EL CILINDRO HUECO

El flujo de calor a través del sistema

(Ti  Te )
q
Rcri  RK 1  RK 2  Rcre

O sea:

Ti  Te
q
1 Lnr2 r1 Lnr / r2 1
  
h cri . .Di .L 2. .K1 .L 2. .K .L h cre .2. .r.L
 Para que q sea máxima, el denominador tiene que ser mínimo por lo
tanto:

 ( RTOTAL )
0
r

( Rcri  RK 1  RK 2  Rcre
0
r

1 1 1 1
.( )  ( 2 )  0
2. .K 2 .L r 2. .h cre .L r

K2
r RADIO CRITICO AISLANTE
h cre
CILINDRICO.

TABLA DE RADIOS CRITICOS

FORMA SÓLIDA RADIO CRITICO

CILINDRICA r  K 2 / h cre

ESFERICA r  2.K 2 / h cre

- ANALISIS GRAFICO PARA EL CILINDRO HUECO
NOTA.- El procedimiento matemático y gráfico para la esfera similar.

EJERCICIO:

Un cable eléctrico de 1.27 cm. de diámetro va a ser aislado con hule

(K=0.1557 W/m.k); el cable estará instalado en el aire ( h cram  8.5W / .k ) al
21.1°C. investigar el efecto que el espesor del aislamiento produce sobre la
disipación de calor, suponiendo una temperatura de 65°C en la superficie
del cable.

SOLUCION.-
a) CALCULO DE RADIO CRITICO

k 0.1557(W / m.k )
rc   0.018317m
h cramb 8.5(W / m 2 K )

ANALISIS: El radio del conductor es 1.27/2 = 0.637 cm.

rconductor  rcrítico

En consecuencia el incremento de espesor de aislante inicialmente

aumentará la disipación de calor y después descenderá.

b) CALCULO DEL CALOR DISIPADO Y SUS ESPESORES DE

AISLANTE CORRESPONDIETE

Para el sistema el flujo de calor está dado por:

(T1  T )
q
In(r / r1 ) 1

2. .K ais .L h cramb .2. .r.L

q 2. .k ais .(T1  T )


L K ais
In(r / r1 ) 
h cramb .r

q 2 xx0.1557 x(65  21.1)


L In(r / 0.00635)  0.1557 /(8.5 xr )

CALOR DISIPADO SIN AISLANTE (r = 0.00635 m.)

Reemplazando este radio en la ecuación de flujo de calor se tiene:

q
 14.88W / m
L
 CALOR DISIPADO CON RADIO CRITICO DE AISLANTE

r  0.01831m)

Reemplazando el radio crítico en la ecuación de flujo de calor se

tiene:

q
 20.85W / m
L

RADIO A PARTIR DEL CUAL COMIENZA EL EFECTO AISLANTE

(q / L  14.88W / m)

Resolviendo la ecuación de calor con q/L = 14.88 W/m y dando

valores para el radio se tiene:

r q/L TERMINO f(r)

0.02 14.888 20.810

0.04 14.888 18.685
0.06 “ 16.834
0.08 “ 15.540
0.10 “ 14.608
0.095 “ 14.810
0.0935 14.888 14.884
 CONCLUSION

La disipación de calor vs radio del aislante se muestra en el gráfico

siguiente:

1.1.2.5 Diversas ecuaciones unidimensionales de aplicación

práctica

- Cilindros huecos

K . Am
qK  .(T1  T2 )
e

Donde:

Am = (A1 + A2)/2 Si A2/A1  1.5

Am = (A2 - A1)/Ln(A2/A1) Si A2 / A1  1.5

- ESFERAS HUECAS

K . Am
qK  .(T 1T2 )
e
 Donde:

Am = A1 .A2

e = r2 - r1

- Cajas huecas de paredes gruesas y espesor uniforme

(e)

K . Am
qK  .(T1  T2 )
e

En ellas se pueden presentar los casos siguientes:

 Caso 1.- Si las longitudes de las tres dimensiones interiores

son mayores que (el 5) y menores que (2.e)

Am  A1  (0.54).e  (1.2).e 2

A1  Area de la superficie interior

  Sumas de las aristas interiores

 Caso 2.- Si la longitud de una de las tres dimensiones

interiores son menores que (e/5)

Am  A1  (0.465).e.  (0.35).e 2

  suma de las ocho aristas interiores mayores que

(e/5)

 Caso 3.- Si la longitud de dos de las dimensiones interiores

son menores que (e/5)
 (2.78).(Ymax .e) (6.384).(Ymax .e)
Am  
log 10 ( A2 / A1 ) In( A2 / A1 )

Ymax = Dimensión interior mayor que (e/5)

 Caso 4.- Cuando las longitudes se las 3 dimensiones

interiores son menores que (e/5)

Am  (0.75). A1 . A2

1.1.2.7 LA CONDUCCION DE CALOR A TRAVES DE

TRAYECTORIAS CURVILINEAS

El flujo de calor en algunos casos puede seguir una

trayectoria no lineal y la ecuación de flujo de calor enunciada por
Fourier se puede adaptar perfectamente en estos caos dándole
la siguiente forma:

q   K . A.(dT / ds

Donde:

(dT/ds) = Gradiente de temperaturas en la

trayectoria curvilínea.

Ilustraremos su uso a través del siguiente ejercicio:

EJERCICIO

Para la pieza metálica proveniente de un cilindro hueco, con k

= 200 W/mk y aislada en su superficie interior y exterior,
calcular el flujo de calor por unidad de longitud Z a través de
la pieza: (considerar la longitud Z muy grande)
 

SOLUCION

 Como la pieza está aislada, los gradientes radial y longitudinal

son iguales
a0

(dT / dr )  (dT / dZ )  0

 La ecuación de Fourier se escribe

q   K . A.(dT / ds)

 Para una porción diferencial se tiene:

dq   K . A.(dT / ds)

 Se puede observar que:

dA = Z.dr y

dT d dT 1
(dT / ds)  .  .
d ds d r

 Reemplazando en la ecuación diferencial del Fourier

 dr dT
dq   K .Z . .
r d

 Integrando el calor entre r1 yr2

r2
q   K .Z . (dr / r ).(dT / d )
r1

q  K .Z .In(r1 / r2 ).(dT / d )

Calculo de (dT/d  ) usando las coordenadas cilíndricas y

considerando la variación de temperaturas en  se tiene:

(1 / r 2 ).( 2T /  2 )  0

+ Integrando:

(T /  )  C1

T  C1 .  C2 - - ( )

+ Condiciones de contorno

 0

  (3. ) / 4 para T  100C

+ Reemplazando en ( )

0  C1 .0  C2 
 C2  0

100  C1  (3. ) / 4  0 
 C1  42.44

+ Entonces:

(dT / d )  42.44K

 Reemplazando valores en el flujo de calor se tiene:

q  200 x1xIn(0.5 / 0.8) x42.44

q  3.989.3W / m
 NOTA.- El signo negativo indica que el flujo de calor va de
100 a 0°C

1.1.3 CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA DE

CALOR

Para la obtención de las ecuaciones unidimensionales con generación

interna bastará con incluir el término de dicha generación interna en las
ecuaciones unidimensionales de la tabla 1.2, el resultado de esta inclusión
se aprecia en la tabla 1.3.

La generación interna se presenta como consecuencia de:

 Paso de la corriente eléctrica en los conductores y

resistencias eléctricas.
 Flujo de Neutrones y otras partículas a través de las paredes
sólidas (reactores nucleares)
 Procesos de combustión (aparatos térmicos)
 Procesos metabólicos de transformación orgánica en calor
 ECUACIONES DE CONDUCCION UNIDIMENSIONAL CON
GENERACION INTERNA

DIRECCION
SISTEMAS DE DEL FLUJO ECUACIONES UNIDIMENSIONALES CON
COORDENADAS DE CALOR GENERACION INTERNA

X q X
 .dx.d  q ''' .dv.d  0
x
CARTESIANAS
qY
Y  .dY .d  q ''' .dv.d  0
X

Z q Z
 .dZ .d  q ''' .dv.d  0
Z

r q r
 .dr.d  q ''' .dv.d  0
r
CILINDRICAS
q Z
Z  .dZ .d  q ''' .dv.d  0
Z

 q
 .d .d  q ''' dv.d  0


q
  .d .d  q ''' .dv.d  0

ESFERICAS
r q r
 .dr.d  q ''' .dv.d  0
r

 q
 .d .d  q ''' .dv.d  0


Tabla 1.3
 El estudio de la generación interna lo realizaremos a través de los 3
modelos básicos no extendiendo aún el análisis a sistemas con otras
geometrías.

1.1.3.1 CONUDUCCION UNIDIMENSONAL CON GENERACION INTERNA

EN “PLACAS PLANAS”

Para la placa plana se usarán las coordenadas cartesianas, estudiando

dicha geometría con sus condiciones de contorno y analizando después
su funcionamiento dentro del medios sólidos y/o fluidos, que son los
que realmente condicionan su funcionamiento, es decir sus contornos y
sus parámetros internos.

- Ecuación generalizada en función de sus contornos

Consideramos una placa plana de espesor L y con temperaturas

superficiales T1 y T2 y que se encuentran generando q ''' (W / m 3 )

 LA ECUACION DIFERENCIAL

 La ecuación unidimensional es:

 (q x / X ).d  q ''' .dv.d  0 ………………..(1)

La ecuación de Fourier para el flujo de calor:

q   K . A.(dT / dX ) …………………………. ..(2)

 con (2) en (1) :

  ( K . A.(dT / dX )) / X .dX .d  q ''' . A.dX .d  0

 Ordenado convencionalmente:

 2T / X 2  q ''' / K

 ECUACION DE DISTRIBUCION DE TEMPERATURAS

 Integrando la Ecuación Diferencial

+ Primera integración:

(T / X )  (q ''' / K ).X  C1

+ Segunda integración:

 
T   q ''' /( 2.K ) . X 2  C1 . X  C2 …………………….(3)

+ Condiciones de contorno

T  T1 en X = 0

T  T2 en X = L

+ Reemplazando en la distribución de temperaturas:

T1  C 2

 
T2   (q ''' .L2 ) /( 2.K )  C1 .L  C2

+ Las constantes de integración son:


C1  (T2  T1 ) / L  (q ''' .L) /( 2.K ) 
C2  T1

 Reemplazando las constantes en la ecuación (3)

   
T   q ''' /( 2.K ) . X 2  (T2  T1 ) / L  (q ''' .L) /( 2.K ) . X  T1
* POSICION DE LA TEMPERATURA MAXIMA

Como la segunda derivada (d 2T / dX 2 )  0 entonces la función tiene un

máximo en (dT / dX )  0

Luego:

dT / dX   (q ''' / K ).X  C1  0

   
  q ''' / K . X  (T2  T1 ) / L  (q ''' .L) /( 2.K )  0

De aquí:

 
X T max  (T2  T1 ) / L  (q ''' .L) /( 2.k ) . k / q ''' 

* TEMPERATURA MAXIMA

 Reemplazando X Tmaz se tiene:

   
Tmax   q ''' /( 2.k ) . X T2 max  (T2  T1 ) / L  (q ''' .L) /( 2.k ) . X T max  T1

* FLUJO DE CALOR EN LOS CONTORNOS DE LA PLACA PLANA

+ en el plano X = 0

q X 0
 K . A.(dT / dX ) X 0

Si:

(dT / dX ) X 0 
 C1  (T2  T1 ) / L  (q ''' .L) /( 2.k ) 
Entonces:

q X 0 
  K . A. (T2  T1 ) / L  (q ''' .L) /( 2.K ) 
+ En el plano X = L

q X L   K . A.(dT / dX ) X L

Si:

(dT / dX ) X L 
  (q ''' .L) /( 2.K ) 

Con estas ecuaciones generalizadas, se pueden resolver situaciones más

sencillas de funcionamiento y algunas de estas gormas de funcionamiento
las hemos agrupado en los denominados casos particulares.

- CASOS PARTICULARES DE LA ECUACION GENERALIZADA

a) PLACA PLANA DE CONTORNOS IGUALES (SIMETRICA)

Sea una placa plana con contornos T1  T2 y de ancho “L” como se

muestra.

* ECUACION DE DISTRIBUCION DE TEMPERATURAS

 Con T1  T2 en la ecuación general

    
T   q ''' /( 2.K ) . X 2  (q ''' .L) /( 2.K ) . X  T1

* POSICION DE LA TEMPERATURA MAXIMA

 Con T1  T2 en la ecuación X T max se tiene:

X T max  L / 2

* ECUACION DE LA TEMPERATURA MAXIMA

 Con T1  T2 y X T max = L/2 se tiene:


Tmax  T1  (q ''' .L2 ) /(8.k ) 

* ECUACIONES DE FLUJO DE CALOR EN LOS

CONTORNOS.

+ En el plano X = 0

q X 0 
  K . A. (q . A.L) /( 2.K )
'''


q X 0  (q ''' . A.L) / 2  (q ''' .V ) / 2

+ En el plano X=L

q X L 
  K . A.  (q ''' .L) / K  (q ''' .L) /( 2.K ) 
q X L  (q ''' .L) / 2  (q ''' .v) / 2

CONCLUSION

+ La placa plana simétrica tiene distribución de temperaturas

parabólica.

+ La temperatura máxima se encuentra en el plano central.

 + El flujo de calor abandona ambos lados de la placa y en la
misma proporción.

b). PLACA PLANA CON EXTREMO ADIABATICO

Sea una placa con un extremo adiabático con en X = 0, con

contornos a temperaturas T1 y T2 y de ancho “L” como se muestra:

* CARACTERISTICAS

+ Extremo izquierdo adiabático, es decir:

En X = 0 se cumple (dT/dX) X 0  0, y q X 0 0

+ El Calor generado por la placa fluye a través del plano X = L:

CONCLUSION: T1  Tmax

* CONDICION PARA EL EXTREMO ADIABATICO (q X 0  0)

 Si

 
q X 0   K . A. (T2  T1 ) / L  (q ''' .L) /( 2.K )  0
  Entonces:

(T2  T1 ) / L  (q ''' .L) /(2.k )

* ECUACION DE DISTRIBUCION DE TEMPERATURAS

 Reemplazando la condición de extremo adiabático en la

ecuación de distribución de temperaturas (ecuación  )

Se obtiene:

 
T  T1  q ''' /( 2.k ) . X 2

* POSICION DE LA TEMPERATURA MAXIMA

Si ( 2T / X 2 )  0 entonces (T / dX )  0 tiene un máximo

  
(T / X )   (q ''' / k ).x  (T2  T1 ) / L  (q ''' .) /( 2.k )  0 
Despejando X

 
X T max  (T2  T1 ) / L  (q ''' .L) /( 2.k ) . k / q ''' 

Utilizando la condición de q X 0 0

X T max  0

* ECUACIONDE LA TEMPERATURA MAXIMA

Reemplazando X T max  0 en la ecuación de distribución de

temperaturas:

   
Tmax   q ''' /( 2.k ) .0 2  (T2  T1 ) / L  (q ''' .L) /( 2.k ) .0  T1

Tmax  T1
* ECUACIONES DE FLUJO DE CALOR EN LOS CONTORNOS

+ En el plano X = 0

q X  L   K . A.dT / dX  X 0  0

+ En el plano X = L

q X  L   K . A.dT / dX  X L


q X  L   K . A. (T2  T1 ) / L  (q ''' .L) /( 2.K ) 
Utilizando la condición para q X 0  0, se tiene:

q X L  q ''' .v

- ECUACION GENERALIZADA EN FUNCION DE LOS FLUIDOS EN

SUS CONTORNOS

En este caso el comportamiento ordinario de la placa plana es

cediendo calor a ambos medios fluidos. Sin embargo cuando la
generación interna es pequeña o cuando la relación entre algunos
parámetros (k,h,L) no es la usual puede suceder que la placa
generadora de calor a un solo fluido; este calor cedido corresponde
a la suma del calor generado más el que proviene del otro fluido.

Se la placa plana con generación interna, de espesor L en contacto

con 2 fluidos diferentes a T1 , h cr1 ; T 2 , h cr 2 , respectivamente.

Las dos situaciones físicas que se producirán son:

 CASO 1 CASO 2

En ambos casos, podemos escribir para el flujo de calor a cada lado de la

placa la ecuación de transferencia correspondiente:

* ECUACIONES DE TRANSFERENCIA DE CALOR

q X 0  K . A.(dT / dX ) X 0  h cr1 . A.(T1  T1 )

q X  L   K . A.(dT / dX ) X L  h cr 2 . A.(T2  T 2 )

NOTA.- El flujo de calor en plano X = 0 para la generación en ambos casos debe

modificar su signo para igualarlo al flujo de calor por convección.

* CALCULO DE LAS TEMPERATURAS SUPERFICIALES EN FUNCION

DEL AS TEMPERATURAS DE LOS FLUIDOS

 Desarrollando las ecuaciones con la generación interna


h cr1 . A.(T1  T1 )  K . A. (T2  T1 ) / L  (q ''' .L) /( 2.K ) 


h cr 2 . A.(T2  T 2 )   K . A. (T2  T1 ) / L  (q ''' .L) /( 2.k ) 

 Desarrollando y resolviendo para T1 yT2 se tiene:

q ''' L2 .(2.L.h cr 2  4.k )  (2.L.h cr1 .T1 ).(2.Lh cr 2  2.K )  4.L.k.h cr 2 .T 2 
T1 

(2.L.h cr1  2.k ).(2.L.h cr 2  2.K )  4k 2 

q ''' .L2 .(2.L.h cr1  4.k )  (2.L.h cr 2 .T 2 ).(2.L.hcr1  2.k )  4.L.k.h cr1 .T1 
T2 

(2.L.h cr1  2.k.).(2.L.h cr 2  2.k )  4.k 2 
 2.L3 .q ''' .(h cr1  h cr 2 )  4.L2 .h cr1 .h cr 2 .(T 2  T1 )
(T2  T1 ) 

(2.L.h cr1  2.k ).(2.L.h cr 2  2.k )  4k 2 

*** CASOS PARTICULARES

a) Placa Plana en contacto con fluidos iguales

 Si los fluidos son iguales se cumple que:

+ T1  T 2  T

+ T1  T 2  T

 Reemplazando estas condiciones en las en ecuaciones de

 
T1  T2  (q ''' .L) /( 2.h cr )  T

b) Placa plana con extremo adiabático y el otro en contacto con fluido

 Del balance de calor:

q X L  q ''' .v  h cr . A.(T2  T )


T2  (q ''' .L)h cr  T 

- RESUMEN DE ECUACIONES PARA LAS PLACAS PLANAS CON

GENERACION INTERNA DE CALOR.

Las fórmulas de aplicación son presentadas en las figuras 1.1.3.1.a y

1.1.3.1.b La primera nos indica las fórmulas para la ecuación
generalizada y los casos particulares.
La segunda tabla nos presenta los valores de las temperaturas T1 ,T2 y
T2  T1 ( en función de las temperaturas de los fluidos) y que pueden ser
reemplazadas en las ecuaciones de la Tabla 1.1.3.1.a. en caso de
necesitarse esta vinculación.
 CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA - PLACAS PLANAS: ECUACIONES EN FUNCION A SUS CONTORNOS

PLACA PLANA GENERICA    

T   q ''' /( 2.k ) .k 2  (T2  T1 ) / L  (q ''' .L) /( 2.k ) . X  T1


X T max  (T2  T1 ) / L  (q ''' .L) / 2.k ) . k / q '''  
   
Tmax   q ''' /( 2.k ) .k 2 T max  (T2  T1 ) / L  (q ''' .L) /( 2.K ) . X T max  T1


q X 0   K . A. (T2  T1 ) / L  (q ''' .L) /( 2.K ) 

q X 0   K . A. (T2  T1 ) / L  (q ''' .L) /( 2.K ) 
PLACA PLANA SIMETRICA   
T   q ''' /( 2.K ) .K 2  (q ''' .L) / 2.k . X  T1 

X T max  L / 2


Tmax  (q ''' .L2 ) /(8.K )  T1 
  
q X 0   (q ''' . A.L) / 2   (q ''' .V ) / 2 
  
q X  L  (q ''' . A.L.) / 2  (9 '''.V ) / 2 
PLACA PLANA DE EXTREMO ADIABATICO T  T  q /( 2.K ). X
1
''' 2

X T max  0
Tmax  T1
q X 0  0
T1  T2  q ''' /( 2.k ) condición adiabática
q X L  q ''' .V

TABLA 1.1.3.1.a.
 ECUACIONES DE TEMPERATURAS SUPERFICIALES EN FUNCIUON DE LAS TEMPERATURAS DE LOS FLUIDOS QUE LO RODEAN

PLACA PLANA GENERICA

q ''' .L2 .( L.h cr 2  2.K )  2.( L.h cr1 .T1 ).( L.h cr 2  K )  L.K .h cr 2 .T 2 
T1 

2. L.h cr1  K ).( L.h cr 2  K )  K 2 

q ''' .L2 .( L.h cr1  2.K )  2.( L.h cr 2 .T 2 ).( L.h cr1  K )  L.K .h cr1 .T1 
T2 

2. ( L.h cr1  K ).( L.h cr 2  K )  K 2 
q ''' .L.3 (h cr1  h cr 2 )  2.L.2.h cr1 .h cr 2 .(T 2  T1 )
(T2  T1 ) 

2. ( L.h cr1  K ).( L.h cr 2  K ).  K 2 
PLACA PLANA SIMETRICA

 
T1  T2  (q ''' .L) /( 2.h cr )  T

PLACA PLANA DE EXTREMO

ADIABATICO

 
T2  (q ''' .L) / h cr  T

TABLA 1.1.3.1.b.
EJERCICIO

Para la placa con generación interna calcular:

a) La posición de la temperatura máxima

b) La temperatura máxima
c) El calor disipado en los planos X = 0 y X = L

Si opera en dos condiciones diferentes como se muestran en la figura.

Considere en ambos casos para la placa plana.

K = 25 W/mK; q '''  2 x10 4 ; L  0.15mt

Largo = 1m; ancho = 1mt.

CONDICION DE OPERACIÓN CONDICION DE OPERACIÓN

N° 1 N°2

h cr1  h cr 2  5W / m 2 K h cr1  5W / m 2 K
T1  T 2  20C h cr 2  10W / m 2 K
T1  T 2  20C
SOLUCION:

* CONDICIONES DE OPERACIÓN N° 1

 Posición de la temperatura máxima

Como la placa es simétrica en su funcionamiento

X max  ( L / 2)  (0.15 / 2)  0.075mt

 Cálculo de la temperatura máxima:

Si:

  
T1  (q ''' .L) /( 2.h cr )  T  (2 x10 4 x0.15) /( 2.5)  20 

T1  320C

  
Tmax  (q ''' .L2 ) /(8.K )  T1  (2 x10 4 x0.152 ) /(8x25)  320 

Tmax  322.25C

 El flujo de calor en los extremos

+ En el plano X = 0

q X L (q ''' .V ) / 2  (2 x10 4 x1x1x0.15) / 2

q X 0  1,500Watt

+ En el plano X = L

q X L  1500Watt

* CONDICION DE OPERACIÓN N° 2

 Posición de la temperatura máxima

Si:

q ''' .L3 .(h cr1  h cr 2 )  2.L2 .h cr1 .h cr 2 .(T 2  T1 )

(T2  T1 ) 
2. ( L.h cr1  K ).( L.hcr 2  K )  K 2 
 2 x10 4 x0.153 x(5  10)  2 x0.15 2 x5 x10 x(20  20)
T2  T1 

2 x (0.15 x5  25).(0.15 x10  25)  25 2 
T2  T1  2.94C

Entonces:

 
X T max  (T2  T1 ) / L  (q ''' .L) /( 2.K ) . K / q ''' 

 
 (2.94 / 015)  (2 x10 4 x0.15) / 2 x25 . 25 / 2 x10 4 
X T max  0.0505m

 Cálculo de la temperatura máxima.

Si:

q ''' .L2 .( L.hcr 2  2.K )  2.( L.h cr1 .T ).( L.h cr 2  K )  L.K .hcr 2 .T 2
T1 

2. ( L.hcr1  K ).( L.h cr 2  K )  K 2 
2 x10 4 x0.15 2 x(0.15 x10  2 x25)  2 x(0.15 x5 x20).(0.15 x10  25)  0.15 x25 x10 x20


2 x (0.15 x5  25) x(0.15 x10  25)  25 2 
T1  215.4C

   
Tmax   q ''' /( 2.K ) . X T2 max  (T2  T1 ) / L  (q ''' .L) /( 2.K ) . X T max  T1

  
  2 x10 4 /( 2 x25) x0.052  (2.94 / 0.15)  (2 x10 4 x0.15) /( 2 x25) x0.05  21.8 
Tmax  216.4C
1.1.3.2. CONDUCCION UNIDIMENSIONAL CON GENERACION
INTERNA EN CILINDROS.

Función generalizada en función de sus contornos sea un

cilindro hueco de radios r1 y r2 y con temperaturas
superficiales T1 y T2 como se muestra en la figura.

* ECUACION DIFERENCIAL

 Usando coordenadas esféricas

 qr / r .dr.d  q ''' .dv.d  0

 La ecuación coordenadas esféricas

 qr / r .dr.d  q ''' .dv.d  0

 La ecuación de Fourier en la dirección radial

q   K . A.(dT / dr )

  K .2. .r.L.(dT / dr )

 Reemplazando:

  K .2. .r.L.(dT / dr ).dr  q ''' .2. .r.dr.L

  Simplificando y ordenando conveniente

 T / r  1/ r . T / r   q

2 2 '''
/K 

 ECUACION DE DISTRIBUCION DE TEMPERATURAS

 Integrando la ecuación diferencial

* Primera integración

  
r.T / r    q ''' / K . r 2 / 2  C1

* Segunda integración

 
T   (q ''' .r 2 ) /( 4.K )  C1 Inr  C2

* Condiciones de contorno:

r  r1 para T  T1

r  r2 para T  T2

* Las constantes de integración son:

 
C1  (T1  T2 )  q ''' (r12  r22 ) /( 4.K ) /In(r1 / r2 )

   
C2  T2  (q ''' .r22 ) /( 4.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .In(r2 ) / In(r1 / r2 )

* Reemplazando en la segunda integración:

   
T  T2  q ''' .(r 2  r22 ) /( 4.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .In(r / r2 ) / In(r1 / r2 )

 POSICION EN LA TEMPERATURA MAXIMA

Como  T / r   0
2 2

Entonces:
 T / r   0 tiene el valor máximo

T / r   (q ''' .r ) /(2.K ) (T1  T2 )  q ''' .(r12  r22 ) /(4.K ).1/(r.In(r1 / r2 ))  0

 
rT max  (2.K ) /( q ''' .In(r1 / r2 )) . (T1  T2 )  q ''' .(r12  r22 ) /( 4.K )   1/ 2

* TEMPERATUTA MAXIMA

   
Tmax  T2  q ''' .(rt 2max  r22 ) /( 4.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .In(rt max / r2 ) / In(r1 / r2 )

* FLUJO DE CALOR EN LOS CONTORNOS DEL CILINDRO

 en la superficie de radio r  r1

qr r1   K . A.(dT / dr ) r r1

   
q r r1  (2. .K .r1 .L). (q ''' .r1 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 /( r1 .In(r1 / r2 ))

- CASOS PARTICULARES DE LA ECUACION GENERALIZADA

a). Caso del cilindro macizo.

Sea un cilindro macizo que genera calor con radio exterior r2 y

temperatura de superficie exterior T2 como se muestra.

q r  r2
 T2

r2

 CARACTERISTICAS

   
 (q ''' .r1 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 / r1 .In(r1 / r2 )
q r r1 0  0 

r1  rt max  0

T1  Tmax

 ECUACION DE DISTRIBUCION DE TEMPARATURAS

o Reemplazando la condición de q r r1 0  0 y con

r1  0 se tiene:


T  T2  q ''' .(r 2  r22 ) /( 4.K ) 

 POSICION DE LA TEMPERATURA MAXIMA

rt max  r1  0

 TEMPERATURA MAXIMA
 o Usando la condición q r  r1 0  0, r1  0, rt max  0


Tmax  T2  (q ''' .r22 ) /( 4.K ) 

 FLUJO DE CALOR EN EL CILINDRO MACIZO

q r  r2   K . A.(dT / dr ) r r2

  
  K .2. .r2 .L.  (q ''' .r2 ) /( 2.K )  (T1  T2 )  q ''' .r12 ) /( 2.k.r2 ) 

o Con r1 = 0 se tiene:

q r  r2  q ''' .r22 . .L  q ''' . volumen del cilindro sólido

b). Otros casos.-

Además de este caso particular se ha realizado el estudio de otros

casos cuyas expresiones se encuentran tabuladas en el gráfico 1.1.3.2
 CONDUCCION UNIDIMENSIONAL CON GENERACION EN CILINDROS –ECUACION GENERALIZADA Y CASOS PARTICULARES

ECUACION GENERALIZADA    
T  T2  q ''' .(r 2  r22 ) /( 4.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .In(r / r2 ) / In(r1 / r2 )

  
rT max  (2.K ) /( q ''' .In(r1 / r2 )) . (T1  T2 )  q ''' .(r12  r22 ) /( r.K ) 
1/ 2

   
Tmax  T2  q ''' .(rt 2max  r22 ) /( 4.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .In(rt max / r2 ) / In(r1 / r2 )

   
q r r1  2. .k.r1 .L. (q ''' .r1 ) /( 2.k )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 /( r1 .In(r1 / r2 ))

   
q r r2  2. .K .r2 .L. (q ''' .r2 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 / r2 .In(r1 / r2 )

ECUACIONES PARA EL
CILINDRO MACIZO 
T  T2  q ''' .(r 2  r22 ) /( 4.K ) 

rt max  r1  0


Tmax  T2  (q ''' .r22 ) /( 4.K ) 
qr r2  q ''' . volumen cilindro macizo

   
 (q ''' .r1 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 / r1 .In(r1 / r2 )
q r r10  
Condición de flujo de calor en el Centro del cilindro macizo
 CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA EN CILINDROS
ECUACUON GENERALIZADA Y CASOS PARTICULARES

CILINDRO EN SUPERFICIE
INTERIOR ADIABATICA    
T  T2  q ''' .(r 2  r22 ) /( 4.K )  (q ''' .r12 ) /( 2.K ) .In(r / r2 )

rt max  r1

   
Tmax  T1  T2  q ''' .(r12  r22 ) /( 4.K )  q ''' .(r12  r22 ) /( 4.K ) .1 / r1 .In(r1 / r2 )

qr r2  q ''' . .L.(r22  r12

   
qr r1  0  (q ''' .r1 ).(2.K ) .  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 / 1 .In(r1 / r2 )
CONDICION DE SUPERFICIE INTERIOR ADIABATICA

CILINDRO CON SUPERFICIE    

T  T2  q ''' .(r 2  r22 ) /( 4.K )  (q ''' .r22 ) /( 2.K ) .Ln(r / r2 )
EXTERIOR ADIABATICA

rt max  r2

  
Tmax  T2  T1  (q ''' .r22 .In(r1 / r2 )) /( 2.K )  q ''' .(r12  r22 ) /( 4.K ) 
q r  r1  q ''' . .L.(r22  r12 )

   
  (q ''' .r2 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .1 / r2 .In(r1 / r2 )
qr r2  0 
CONDICION DE SUPERFICIE EXTERIOR ADIABATICA
TABLA 1.1.3.2
 CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA DE CILINDROS
ECUACIONES DE TEMPERATURAS SUPERFICIALES EN FUNCION DE LOS FLUIDOS

CILINDRO CON FLUIDO INTERIOR

K .C  B  ( A.T 2 )  K  A.

. E  C  ( D.T1 )
T1 

( K  A).( K  D)  K 2 
K .E  C  ( D.T1 )  K  D
. C  B  ( A.T 2 )
T2 

( K  A).( K  D)  K 2 
Ver A,B,C,D, al final de la tabla
CILINDRO MACIZO


T2  T 2  (q ''' .r2 ) /( 2.h cr 2 ) 

CILINDRO CON SUPERFICIE

INTERIOR ADIABATICA


T2  T 2  q ''' .(r22  r12 ) /( 2.r2 .h cr 2 ) 

TABLA 1.1.3.2.a
 CONDUCCION UNIDIMENSIONAL CON GNEERACION INTERNA EN CILINDROS
ECUACIONES PARA LAS TEMPERATURAS SUPERFICIALES EN FUNCION DE LOS FLUIDOS.

CILINDRO CON SUPERFICIE

EXTERIOR ADIABATICA


T1  T1  q ''' .(r22  r12 ) /( 2.r1 .h cr1 ) 

VALORES DE LAS EXPRESIONES A  h cr 2 .r2 .In(r1 / r2 )

A, B, C y D

 
B  (q ''' .r22 ) / 2 .In(r1 / r2 )

 
C  q ''' / 4 . r12  r22 
D  h c1 .r1 .In(r1 / r2 )

 
E  (q ''' .r12 ) / 2 .In(r1 / r2 )

NOTA: K= Conductividad térmica del sólido

TABLA 1.1.3.2.b
1.1.3.3 CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA
EN ESFERAS

- ECUACION GENERALIZADA EN FUNCION DE SUS CONTORNOS

Para la esfera hueca que se muestra si se encuentra generando calor

a razón de q ''' , con radios r1 , r2 y temperaturas superficiales T1 y T2 ,
se encontraran las ecuaciones típicas de distribución de
temperaturas, temperatura máxima, posición para la temperatura
máxima y flujos de calor.

AESFERA  4. .r 2 VESFERA  (4 / 3). .r 3

V  4. .r 2 .dr

* ECUACION DIFERENCIAL

o La conducción en la Dirección Radial

 dqr / dr .dr.d  q ''' .dV .d  0

o La ecuación de conducción de Fourier

q   K . A.(dT / dr )   K .4. .r 2 .(dT / dr )

 o Reemplazando el dV, en la ecuación de Fourier y simplificando:



  r 2 .(dT / dr )


q ''' .r 2
r K

* ECUACION DE DISTRIBUCION DE TEMPERATURAS

+ Primera Integración

 
r 2 .dT / dr    (q ''' .r 3 ) /(3.K )  C1

+ Segunda integración

 
T   (q ''' .r 2 ) /(6.K )  (C1 / r   C2

o Condiciones de contorno

T  T1 para r  r1

T  T2 para r  r2

o Usando estas condiciones se pueden hallar C1 y C2


C1  (r1 .r2 ) /( r1  r2 ). (T1  T2 )  q' ' '.(r22  r12 ) /(6.K ) 

  
C2  T1  (q ''' .r12 ) /(6.K )  r2 (r1  r2 ). (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) 

o Con las constantes; la ecuación de distribución de temperaturas

es:

   
T  T1  q ''' .(r12  r 2 ) /(6.K )  r2 /( r1  r2 ). (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) .1  (r1 / r )
 * Posición de la temperatura máxima

dT / dr   0  (q ''' .r ) /(3.K ) C1 / r 2 

o Reemplazando C1 y despejando r se tiene rt max

  
rt max  3.r1 .r2 .K . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) / q ''' .(r1  r2 ) 
1/ 3

* ECUACIONES PARA EL FLUJO DE CALOR

+ Flujo de calor en r = r1

q r  r1   K . A.(dT / dr ) r r1

   
(dT / dr ) r r1   (q ''' .r1 ) /(3.K )  2  (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r1 .(r1  r2 )

     
q r r1   K .4. .r12 . (q ''' .r1 ) /(3.K )  r2 . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r1 .(r1  r2 )

+ Flujo de calor en r = r2

q r  r2   K . A.(dT / dr ) r r2

   
(dT / dr ) r r2   q ''' .r2 ) /(3.K )  r1 . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r1 .(r1  r2 )

     
q r r2   k.4. .r22 . (q ''' .r2 ) /(3.K )  r1 . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r2 .(r1  r2 )
- CASOS PARTICULARES DE LA ECUACION GENERALIZADA.-

Ha realizado un análisis de los siguientes casos:

a). Esfera maciza

b). Esfera con superficie interior aislada
c). Esfera con superficie exterior aislada.

Los resultados se muestran en la Tabla 1.1.3.3.

- ECUCIONES EN FUNCION DE LOS FLUIDOS EN CONTACTO CON

LAS SUPERFICIES DE LA ESFERA

Se ha aplicado la ecuación generalizada, a la esfera maciza y a las

esferas con superficies adiabáticas.

Los resultados también se muestran en la tabla 1.1.3.3.

CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA EN ESFERAS – ECUACION GENERALIZADA Y CASOS PARTICULARES

ESFERA HUECA    
T  T1  q ''' .(r12  r 2 ) /(6.K )  r2 /( r1  r2 ). T1  T2 )  q ''' (r22  r12 ) /(6.K ) .1  (r1 / r )

 
rt max  3.r1 .r2 .K . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) / q ''' .(r1  r2 )  1/ 3

   
Tmax  T1  q ''' .(r12  rt 2max ) /(6.K )  r2 /( r1  r2 ). (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) .1  (r1 / rt max )

     
q r r1   K .4. .r12 . (q ''' .r1 ) /(3.K )  r2 . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r1 .(r1 / r1  r2 )

     
q r r2   K .4. .r22 . (q ''' .r2 ) /(3.K )  r1 . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r2 .(r1  r2 )

ESFERA MACIZA 
T  T2  q ''' .(r22  r 2 ) /(6.K ) 
rt max  0


Tmax  T1  T2  (q ''' .r22 ) /(6.K ) 

q r r 2  q ''' .V V = Volumen de la esfera

q r  r1 0   
 (q ''' .r1 ) /(3.K )  r2 . (T1  T2 )  q ''' .(r22  r12 ) /( 6.K ) /r1 .(r1 .(r1  r2 )
 0 
CONDICION DE FLUJO DE CALOR EN EL CAENTRO DEL CILINDRO

TABLA 1.1.3.3.

CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA EN ESFERAS – ECUACION GENERALLIZADA Y CASOS PARTICULARES
ESFERA CON SUPERFICIE
 INTERIOR ADIABATICA
 
T  T1  q ''' /(6.K ) . r12 .3  ((2.1 ) / r )  r 2 
rmax  r1

 
Tmax  T1  T2  q ''' .r12 .(r1  r2 ) /(3.r2 .K )  q ''' .(r22  r12 ) /(6.K ) 
qr r2 q ''' .4. .(r23  r13 ) / 3  q ''' .(V2  V1 )  VOLUMEN DE LA ESFERA HUECA

   
qr r1  0  (q ''' .r1 ) /(3.K )  r2 . (T1  T2 )  q ''' .(r22  r12 /(6.K ) /r1 .(r1  r2 )
CONDICION DE SUPERFICIE INTERIOR ADIABATICA

ESFERA CON SUPERFICIE

EXTERIOR ADIABATICA    
T  T1  q ''' /(6.K ) . r12  r 2  (2.r23 ) / r1 .(r1 / r )  1

rt max  r2

  
Tmax  T2  T1  q ''' .r22 (r1  r2 ) /(3.r1 .K )  q ''' .(r22  r12 ) /(6.K ) 
 
q r r1   q ''' .4. .(r23  r13 ) / 3  q ''' .(V2  V1 )  q ''' . Volumen de la esfera hueca

   
q r r2  0  (q ''' .r2 ) /(3.K )  r1 . (T1  T2 )  q ''' .(r22  r12 ) /(6.K ) /r2 .(r1  r2 )

CONDICION DE SUPERFICIE EXTERIOR ADIABATICA

TABLA 1.1.3.3.A
CONDUCCION UNIDIMENSIONAL CON GENERACION EN ESFERAS – ECUACIONES DE TEMPERATURAS SUPERFICIALES EN FUNCION
A LOS FLUIDOS
ESFERA HUECA CON 2 FLUIDOS

T1 
A  h .((q
cr 2
'''
 
.r1 ) / 3)  C  h cr1 .T1  B.  ((q ''' .r2 ) / 3)  D  h cr 2 .T 2 
A  h cr 2 . B.h cr1   B.A.


B  h . ((q .r ) / 3)  D  h .T  A.((q .r ) / 3)  C  h
cr1
'''
2 cr 2 2
'''
1 cr1 .T1 
A  h . B  h   B.A
T2
cr 2 cr1
Ver A, B, C, D al final de la tabla
ESFERA MACIZA CON FLUIDO
EXTERIOR


T2  T 2  (q ''' .r2 ) /(3.h cr 2 ) 

ESFERA CON SUPERFICIE

INTERIOR ADIABATICA


T2  T 2  q ''' .(r23  r13 ) /(3.r22 .h cr 2 ) 

TABLA 1.1.3.3.B
 CONDUCCION UNIDIMENSIONAL CON GENERACION INTERNA EN ESFERAS

ESFERA CON SUPERFICIE

EXTERIOR ADIABATICA


T1  T1  q ''' .(r23  r13 ) /(3.r12 .h cr1 ) 

VALORES DE LAS EXPRESIONES A, A  K .r1 /r2 .(r1  r2 )

B, C, D

B  K .r2 /r1 .(r1 .  r2 )

 
C  q ''' .r2 .(r1  r2 ) /6.r1 

 
D  q ''' .r1.(r1  r2 ) /6.r2 
TABLA 1.1.3.3 C
EJERCICIOS

La resistencia eléctrica hueca entrega calor a un fluido interior y exteriormente como

muestra la figura, la resistencia es de acero K  20W / m.k ). si la resistencia genera calor a

razón de: q '''  10 4 W / m 3 , calcular

a) El calor entregado por unidad de longitud en el interior y exterior.

b) La temperatura máxima y su posición.

CONDICIONES DEL FLUIDO

T1  T 2  30C

h cr1  2W / m 2 K

h cr 2  10W / m 2 K

r1  0.05mt

r2  0.1mt

SOLUCION:

* CALCULO DEL FLUJO DE CALOR

o Las constantes son:

A  h cr 2 .r2 .In(r1 / r2 )  10 x0.1xIn(0.05 / 0.1)  0.693

    
B  (q ''' .r22 ) / 2 .In(r1 / r2 )  10 4 / 4 . 0.052  0.12 / 2  xIn(0.05 / 0.1)  34.657

C  q / 4
'''
. r  r   10 / 4
1
2
2
2
. 0.05  0.1   18.75
4 2 2
 D  h cr1 .r1 In(r1 / r2 )  2 x0.05xIn(0.05 / 0.1)  0.0693

   
E  (q ''' .r12 ) / 2 .In(r1 / r2 )  (10 4 x0.052 ) / 2 .In(0.05 / 0.1)  8.664

o Las temperaturas T1 y T2 son:

k .C  B( A.T 2 )  K  A

. E  C  ( D.T1 )
T1 

( K . A).( K  D)  K 2 

20 x 18.75  34.657  0.693x30  20  0.693

.  8.66  18.75  0.069 x30


(20  0.693).(20  0.0693)  20 2 

T1  64.448C

T2  64.063C

o Flujo de calor en la superficie exterior:

   
q r r2  2. .K .r2 .L. (q ''' .r2 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) .

1/( r2 .In(r1 / r2 ))

q
L
r  r2   
 2 xx20 x0.1.  (10 4 x0.1) /( 2 x20)  (64.448  64.063)  10 4 x(0.05 2  0.12 ) /( 4 x20) 
1/(0.1xIn(0.05 / 0.1))

q
r  r2  214.085W / m
L

o El flujo de calor en la superficie interior es:

   
q r r1  2. .K .r1 L. (q ''' .r1 ) /( 2.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) 

1/(r1.In(r1 / r2 )

q
r  r1  21.621W /
L

* CALCULO DE LA POSICION DE LA TEMPERATURA MAXIMA

  
rt max  (2.K ) /( q ''' .In(r1 / r2 )) . (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) 
1/ 2

 
 (2 x20) /(10 4 xIn(0.05 / 0.1)) . (64.448  64.063)  10 4.

.(05 2  0.12 ) /( 4 x20)  1/ 2

rt max  0.05646mt

* CALCULO DE LA TEMPERATURA MAXIMA

  
Tmax  T2  q ''' .(rt 2max  r22 ) /( 4.K )  (T1  T2 )  q ''' .(r12  r22 ) /( 4.K ) . 

In(rt max / r2 ) / In(r1 / r2 )


Tmax  64.063  10 4 x(0.0564 2  0.12 ) /( 4 x20)  (64.448  67.063)
  10 4 x(0.052  0.52  0.12 ) /( 4.20) 
. In(0.05646 / 0.1) / In(0.05 / 0.1) / 

In(0.05 / 0.1) 
Tmax  64.459C

EJERCICIO

Un recipiente esférico de acero (K=20 W/mK), de Smm de espesor contiene

un líquido, cuya temperatura debe ser mantenida en 100°C, para lo cual el
recipiente generará una cantidad de calor necesaria.

a) Calcular el calor generado por el tanque esférico en las condiciones que

se muestran en la figura.
b) Calcular el calor perdido al ambiente.

* CALCULO DE T2


T2  T 2  (q ''' .r2 ) /(3.h cr 2 ) 

 20  (q ''' x0.805) /(3x5) 

T2  20  0.05366q '''

* CALCULO DE CALOR GENERADO

o En la ecuación de la temperatura máxima

  
Tmax  T1  T2  q ''' .r12 .(r1  r2 ) /(3.r2 .K )  q ''' .(r22  r12 ) /(6.K ) 

Tmax  20  0.5366q '''  q ''' x0.8 2 x(0.8  0.805) /(3x0.805x20) 

q '''
.(0.805 2  0.8 2 ) /( 6 x20) 

Tmax  20  0.05366q '''  6.6252 x10 5 xq '''  6.6875 5 q '''

100 = 20 + 0.05366 q’’’

W
q '''  1,490.86
m3

o Cálculo del calor perdido en el medio ambiente

q r  r2 
 q ''' .4. .(r23  r13 ) / 3 

 1,490.86 x4 xx(0.8053  0.83 ) / 3

q r  r2  60.32Watt