UNIVERSIDAD NACIONAL AUTÓNOMA DE

MÉXICO

Facultad de Química

6-10-2016
 Problema

Una barra de aluminio de longitud de 90.5 cm y 1.9 cm de diámetro, se calienta en uno de

los extremos con vapor de caldera a 2 kgf/cm2.
Encontrar la posición longitudinal en la barra donde las diferencias entre los calores de los
mecanismos conductivos y convectivos son despreciables. Reportar la eficiencia de la barra
en esta posición.

Laboratorio de Ingeniería Química II

Diagrama de
flujo
Profesor: Sergio Adrián García González

Práctica No. 4:

“Transferencia de calor de una barra cilíndrica”

Datos Grupo: 05
experimentales

Tabla 1. Datos
experimentales. Alumnos:

Tabla 2. Datos
experimentales de
 Chávez Díaz Óscar Iván
la barra, el vapor  Jacinto Quintana Paulina
y el ambiente.

Cuestionario
Semestre : 2017-1
1.- Elegir una
posición en la
barra y describir
lo que sucede con
la temperatura a
medida que pasa
el tiempo. ¿En qué tiempo se alcanza el estado estable o estacionario?
Explique por qué ha sido posible alcanzar ese estado.
Se eligió la posición de 30.5 cm en la barra, a medida que pasa el tiempo, los puntos se van
haciendo más cercanos entre sí, de esta manera se puede decir que se hacen estables. El
estado estacionario se alcanza a los 60 min, ya que es cuando ya no hay muchas variaciones
en el último segmento de la barra.
Este estado es posible alcanzarlo porque las pérdidas de calor por conducción y convección
obedecen un modelo que se describirá más adelante y que indica que la temperatura ya no
cambia ni respecto a la posición ni respecto al tiempo.

2.- Para una longitud x, ¿cómo es la temperatura desde el centro de la barra a la

superficie de la barra en la posición radial y en el estado estacionario?
Explique por qué.
La temperatura en el centro de la barra es mayor a la de la superficie, pero para fines
prácticos, podemos decir ambas temperaturas son iguales en su posición radial.

3.- Graficar los siguientes datos:

Gráfica 1. Las temperaturas en el régimen transitorio contra la posición x, teniendo como

parámetro el tiempo.
 Temperatura de la barra en función de la distancia
tiempo 0 min tiempo 10 min tiempo 20 min tiempo 30 min
tiempo 40 min tiempo 50 min tiempo 60 min
120

100

80
Temperatura de la barra (ºC)

60

40

20

0
0 10 20 30 40 50 60 70 80 90
distancia (m)

La gráfica 1 muestra los valores de la temperatura respecto a la longitud, se puede observar

que los valores aún no se mantienen constantes, por lo que se encuentran en el régimen
transitorio.

Gráfica 2. Las temperaturas en el régimen permanente contra la posición x.

 Temperatura de la barra en función de la distancia
120

100
Temperatura de la barra (ºC)

80
tiempo 40 min
60 tiempo 50 min
tiempo 60 min
40

20

0
0 10 20 30 40 50 60 70 80 90
distancia (m)

La gráfica 2 muestra los valores de temperatura a partir del minuto 40, se puede observar
que los valores se mantienen más constantes comparados con los valores de la gráfica 1, de
esta manera los valores se encuentran dentro del régimen permanente.

4.- Plantear un balance de energía en Kcal/h, en cualquier posición en la barra

aplicando el principio del balance de coraza a un elemento diferencial de volumen,
introduciendo la ecuación que define el coeficiente convectivo individual de
transferencia de calor Kcal/m2h°C y obtener la ecuación diferencial sin integrar.

Se realizará un balance de coraza para la determinación de la transferencia de calor en el

sistema:

∂T ∂T
kA −kA + hS ( T −T amb )=0
∂ x x+∆ x ∂x x

Donde S es el perímetro multiplicado por la fracción de barra donde se transfiere (S=PΔx),

∂T ∂T
kA −kA refiere al calor transmitido por conducción y hS ( T −T amb ) el
∂ x x+∆ x ∂x x
transferido por convección. Se divide entre Δx la ecuación anterior.

∂T ∂T
kA −kA
∂ x x+∆ x ∂x x
−hP ( T −T amb )=0
∆x
Se aplica el operador límite cuando Δx tiende a cero para pasar a la diferencial:

∂T ∂T
lim
∆ x→ 0
( kA
∂ x x+ ∆ x
∆x
−kA
∂x x
−hP ( T −T amb )=0 )
∂ ∂T
∂x ( kA
∂x )
−hP ( T−T amb ) =0

∂2 T hP hP
2
= ( T −T amb )=0 Donde se define m=
∂ x kA kA

5. En la ecuación diferencial lineal de segundo orden obtenida en el punto anterior,

defina a θ=Ts-Tamb. Ts es la temperatura de la superficie de la barra en cualquier
posición y Tamb es la temperatura del aire ambiente, incorpore esta nueva variable
en la ecuación para que ahora las temperaturas queden en función de θ.

Se realiza un cambio de variable de θ=T −T amb

dθ dT T amb d 2 θ d 2 T
= − ; =
dx dx dx dx 2 dx 2

d2θ
−mθ=0
dx 2

6. Aplicar las condiciones a la frontera a la ecuación.

7. Resuelva la ecuación diferencial en el estado estacionario y reporte las siguientes

condiciones:

7.1 El perfil de temperaturas:

Se sabe que θ=± √m= {− √ m,+ √ m }, por lo tanto, resolviendo la ecuación diferencial se
tiene que:
θ=C1 e−√ m x +C 2 e √ m x

Se tienen estas condiciones de frontera:

1 x=0 ; θ=θ =T −T ¿2 ¿ x=L ; dθ =0 ¿
b v a
¿ dx

Aplicándolas a la ecuación diferencial:

θ=θb=C 1 e−√m 0+C 2 e√ m 0

θ=θb=C 1 +C2
Aplicando la segunda condición:

dθ
=C 1 (− √ m ) e−√ m+ C2 ( √ m ) e √ m=0
dx
− mL mL
C 1 ( √ m ) e √ =C2 ( √ m ) e √

Resolviendo por sustitución

C 1=θb −C2

Se desarrolla y se despeja C2:

√ mL
( θ b−C 2) √ m e−√m L =C2 ( √ m ) e
θb √ m e−√ m L −C2 √ me−√ m L =C 2 √ m e√ m L
θb √ m e−√ m L =C2 ( √ me √ m L + √ me−√ m L )

θ b √m e−√ m L
C 2=
√ m ( e √ m L + e−√m L )
θ b e−√m L
C 2=
( e √m L +e−√ m L )

Sustituyendo el valor en C 1=θb −C2:

θb e−√ m L
C 1=θb − ( e √ m L + e−√m L )
Se desarrolla y despeja para obtener un valor de C1:

θ b ( e√ m L +e−√ m L ) −θb e−√ m L

C 1=
e √m L + e− √m L
θ b ( e√ m L +e−√ m L−e−√ m L )
C 1=
e √ m L +e−√ m L
θb e √ m L
C 1=
e√ m L +e−√ m L
Se sustituyen las constantes en la ecuación inicial para obtener:

θ b e√ m L θb e−√ m L
(
θ= √ m L −√ mL e
e +e
− √m x
)
+ √ m L −√ m L e√ m x
(e + e ) ( )
Se simplifica la ecuación de la siguiente manera:

e √ m L e−√ m L e−√ m L e√ m L
θ=θb
[( e √ m L +e−√ m L
θb
+ )(
( e√ m L +e−√ m L ) )]
θ= √m L −√ m L
( e√ m ( L− x ) +e √ m (x− L ))
e +e

Se usarán la siguiente relación trigonométrica para reducir aún más la ecuación:

e x +e− x
cosh ( h )=
2
Se modifica la ecuación y se sustituye la relación trigonométrica:

θ e √m ( L−x ) +e−√ m ( L−x ) 2

θb
=
e√ m L +e−√ m L 2 ()
θ cosh [ √ m ( L−x ) ]
=
θb cosh [ √ m L ]
cosh [ √m ( L−x ) ]
θ=θb
(cosh [ √ m L ] )
Se vuelve a la ecuación original recordando que θ=T −T amb

cosh [ √ m ( L−x ) ]
T =( T v −T amb )
{ cosh [ √m L ]
+T amb
}
7.2 La cantidad de calor transferido por convección desde la base de la barra:
Qb=√h̅PkATθbtanh(mL). ¿Cuál es el valor numérico en Kcal/h?

Qb= √ h́ Pk A T ∗θb tanh ( mL )

Q b=
(√ 10.639 mKcal2
h° C ) ( 0.05969m ) (197.169
Kcal
m° C)
2
( 0.0002835 m ) ( 89.2° C ) ¿

Kcal
Q b=16.7321
h

7.3La cantidad de calor que entra por la base de la barra utilizando el perfil de
temperaturas y obtener Qbo =√h̅PkATθbtanh(mL). ¿Cuál es el valor numérico en
kcal/h? Explique por qué Qb = Qbo

Q bo= √ h́ Pk A T ∗θb tanh ( mL )

Kcal Kcal
Qbo=
√( 10.639 2
m h° C ) (
( 0.05969 m ) 197.169
m° C )
( 0.0002835 m2) ( 89.2 ° C ) ¿

Kcal
Q bo=16.7321
h

Los valores de Qb y Qbo son iguales porque la transferencia de calor por convección es
exactamente igual tanto para el modelo teórico como el experimental.

7.4 La eficiencia de la aleta Ω. Definida como el calor real transferido por convección
desde la base de la barra entre el calor transferido por convección si la aleta estuviera
en toda su longitud a la temperatura uniforme del vapor Tv. ¿Cuál es el valor
numérico?

tanh ⁡(mL) tanh ⁡(3.05032)

Ω= =
mL 3.05032

Ω=0.326367
8. Grafique las temperaturas en el estado estacionario contra la posición x de la tabla
de datos experimentales a partir de la posición x = 10.5 cm y encuentre el modelo
matemático exponencial o potencial que ajuste mejor los datos. Reporte el modelo
encontrado.

T °C vs x metros
a t 40min Power (a t 40min) a t 50min
a t 60min Power (a t 60min)

f(x)
f(x) == 25.85
27.64 x^-0.54
x^-0.5
R²
R² == 0.98
0.99

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9

El modelo que mejor se ajusto fue el potencial cuyo coeficiente de correlación fue de
0.98529, el modelo es el siguiente:
T=27.636x-0.499

9. Derive el modelo encontrado en el punto anterior dT/dx y sustitúyalo en la ecuación

de conducción de Fourier Qx = -ATk dT/dx, obtenga el perfil puntual al graficar Qx
contra x utilizando como punto inicial el valor de 10.5 cm de la barra de aluminio.
Describir el significado del perfil obtenido.

T=27.636x-0.499

T=28.048x-0.499

dTdx=-13.79x-1.499

Qx=-kAT(-13.79x-1.499)

Qx=kAT(13.79x-1.499)
 Qx en función de x
25

20

15
Qx

10

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
x (cm)

En el gráfico se puede observar que a medida que la longitud aumenta la cantidad de calor,
generado por el mecanismo de conducción, a lo largo de la barra disminuye.

10. Obtenga el perfil puntual de la cantidad de calor transferido por convección desde
la base dela barra Qb(x) contra la posición x a partir del valor de 10.5 cm de la barra
de aluminio.
−hP θb
Qbx = senh [ m ( x −L ) ]
cosh ( mL ) m

Qbx en función de x
14

12

10

8
Qbx

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
x (cm)

La gráfica muestra que mientras la posición (x) aumenta el calor generado por el
mecanismo de convección disminuye.
11.En la expresión de la eficiencia de la aleta del punto 7.4, sustituya el valor al valor
L (longitud de la aleta) por los de las posiciones x. A partir de los valores de 10.5 cm y
grafique los valores.

Ω en función de x
1
0.9
0.8
0.7
0.6
0.5
Ω

0.4
0.3
0.2
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
x (cm)

En la gráfica se puede observar que al aumentar la posición de la barra (x), la eficiencia de

aleta disminuye.
A mayor longitud menor eficiencia.

12.Junte los perfiles de las gráficas 4 y 5 en una misma gráfica y reporte la posición
longitudinal en la barra donde las diferencias entre los calores de los mecanismos
conductivos y convectivos son despreciables y la eficiencia de la barra en esa posición.

Qbx & Qx en función de x

25

20

15
Qbx & Qx

10

0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
x (cm)