Transferencia de Calor

Operaciones de transferencia de Calor 3007820
Departamento de procesos y energía

Facultad de minas
Tarea N° 2 de Transferencia de Calor

Rosero Urbina Erla Anderlein.
Facultad de Minas, Universidad Nacional de Colombia, Medellín, Colombia.
earoserou@unal.edu.co
Febrero 15, 2018
Enunciado
Calcule el tiempo requerido y la energía para llevar un horno que quema gas natural de
25°C a la temperatura de operación. Use método de diferencias finitas. Qué hacer para
reducir el tiempo sin que el refractario tenga problemas. (Asuma diámetro y flujo de gas
para sostener la temperatura).
Desarrollo del Ejercicio
Suposiciones Y Datos Estándar:
 La geometría del horno de operación se ha tomado de forma cilíndrica y se

considera un solo aislamiento (concrax 1500).
 Para el horno de operación se tiene una temperatura a condiciones externa de 25°
C, el cual opera a una temperatura de 1500 °C
 El flujo másico de entrada es de 20 Kg/h
 La k es de 0.46 W/m2 °C y la capacidad calorífica de los gases es de 879 J/Kg °C
 El radio del material refractario es de 0.08 m y la altura del horno es de 3 m
 Se desprecia la transferencia de calor por radiación
 Se supone dirección unidireccional y propiedades constantes.
 Se asumen 6 nodos para el calculo por medio de diferencias finitas los cuales se
distinguen por estar en los bordes, en las esquinas y en los centros.
 El dr y dz son las condiciones para las iteraciones y para ambas se tiene un valor
de 0.01.
A continuación, se presentan las siguientes ecuaciones por el método de diferencias finitas,

ya que en cada uno de los nodos la distribución de calor puede cambiar según sea el caso.
Facultad de minas
 2≤𝑖 ≤𝑚−1𝑦2≤𝑗 ≤𝑛−1

2𝑖 + 1 2𝑖 − 1
𝑇𝑡 (𝑖, 𝑗) = Ϝ0 [( ) 𝑇𝑡 (𝑖, 𝑗 + 1) + ( ) 𝑇𝑡 (𝑖, 𝑗 − 1) + 𝑇𝑡 (𝑖 − 1, 𝑗) + 𝑇𝑡 (𝑖 + 1, 𝑗)] + (1
2𝑖 2𝑖
− 4Ϝ0 )𝑇𝑡 (𝑖, 𝑗)
 𝑗=𝑛𝑦2<𝑚 <𝑚−1
2𝑖 − 1 4𝑖 − 1
𝑇𝑡 (𝑖, 𝑗) = Ϝ0 [( ) 𝑇𝑡 (𝑖, 𝑗 − 1) + 𝑇𝑡 (𝑖 − 1, 𝑗) + 𝑇𝑡 (𝑖 + 1, 𝑗)] + [1 − ( ) Ϝ0 ] 𝑇𝑡 (𝑖, 𝑗)
2𝑖 𝑖
 2<𝑗 <𝑛−1𝑦𝑖 =1
2𝑖 + 1 2𝑖 − 1
𝑇𝑡 (𝑖, 𝑗) = Ϝ0 [( ) 𝑇𝑡 (𝑖, 𝑗 − 1) + ( ) 𝑇𝑡 (𝑖, 𝑗 + 1) + 2𝑇𝑡 (𝑖 + 1, 𝑗)] + (1 − 4Ϝ0 )𝑇𝑡 (𝑖, 𝑗)
2𝑖 2𝑖
 2<𝑗 <𝑛−1𝑦𝑖 =𝑚
2𝑖 + 1 2𝑖 − 1
𝑇𝑡 (𝑖, 𝑗) = Ϝ0 [( ) 𝑇𝑡 (𝑖, 𝑗 − 1) + ( ) 𝑇𝑡 (𝑖, 𝑗 + 1) + 2𝑇𝑡 (𝑖 − 1, 𝑗)] + (1 − 4Ϝ0 )𝑇𝑡 (𝑖, 𝑗)
2𝑖 2𝑖
 𝑖 =1𝑦𝑗 =𝑛
2𝑖 + 1 2𝑖 + 1
𝑇𝑡 (𝑖, 𝑗) = Ϝ0 [( ) 𝑇𝑡 (𝑖, 𝑗 − 1) + 2𝑇𝑡 (𝑖 + 1, 𝑗)] + {1 − [( ) + 2] Ϝ0 } 𝑇𝑡 (𝑖, 𝑗)
2𝑖 𝑖
 𝑖 =𝑚𝑦𝑗=𝑛
2𝑖 + 1 2𝑖 + 1
𝑇𝑡 (𝑖, 𝑗) = Ϝ0 [( ) 𝑇𝑡 (𝑖, 𝑗 − 1) + 2𝑇𝑡 (𝑖 − 1, 𝑗)] + {1 − [( ) + 2] Ϝ0 } 𝑇𝑡 (𝑖, 𝑗)
2𝑖 𝑖
El anterior sistema de ecuaciones se programó en Matlab, con el fin de determinar el tiempo

requerido para la operación total del horno, es decir hasta llegar a la temperatura máxima
de 1500°C. Como se puede resaltar en las ecuaciones, se tiene en cuenta el número de
Fourier por lo que se tuvo en cuenta la siguiente ecuación para el cálculo:
𝑘 × 𝑑𝑡
Ϝ0 =
𝑅𝑜 × 𝐶𝑝 × 𝑑𝑟
A continuación, se presenta el código que se utilizó en Matlab.

Facultad de minas
clear all
clc
%%
% Cantidades iniciales
Tf=1500;
Ti=25;
k=0.46;
r=0.08;
dr=0.01;
z=3.0;
dz=0.01;
n=(r/dr)+1;
m=(z/dr)+1;
cp=879;
RO=20;
dt=40;
%%
%Se definen las ecuaciones de flujo de entrada al horno y se declara una
%matriz para iterar las temperaturas:
Fo=((k*dt)/(RO*cp*dr));
T(1)= Tf;
Ta=zeros(m,n);
%%
%Iterando para llenar matriz de temperaturas:
for j=2:n
for i=1:m
Ta(i,j)=Ti;
end
end
Ta(:,1)=Tf;
%%
%Definiendo y ejecutando bucles de diferencias finitas:
time=0;
while Ta(m,n)<=849.9
time=time+dt;
for j=2:n
for i=1:m
if ((2<=i && i<=(m-1))&& (2<=j && j<=(n-1)));
Ta(i,j)=Fo*((((2*i)+1)/(2*i))*Ta(i,j+1)+(((2*i)-
1)/(2*i))*Ta(i,j-1)+Ta(i-1,j)+Ta(i+1,j))+(1-4*Fo)*Ta(i,j);
elseif (j==n && 2<=i && i<=m-1);
Ta(i,j)=Fo*((((2*i)-1)/(i))*Ta(i,j-1)+Ta(i-
1,j)+Ta(i+1,j))+(1-(((4*i)-1)/i)*Fo)*Ta(i,j);
elseif (2<=j && j<=n-1 && i==1);
Ta(i,j)=Fo*((((2*i)+1)/(2*i))*Ta(i,j-1)+(((2*i)-
1)/(2*i))*Ta(i,j+1)+2*Ta(i+1,j))+(1-4*Fo)*Ta(i,j);
elseif (2<=j && j<=n-1 && i==m);
Ta(i,j)=Fo*((((2*i)+1)/(2*i))*Ta(i,j-1)+(((2*i)-
1)/(2*i))*Ta(i,j+1)+2*Ta(i-1,j))+(1-4*Fo)*Ta(i,j);
elseif (i==1 && j==n);
Ta(i,j)=Fo*((((2*i)+1)/(i))*Ta(i,j-1)+2*Ta(i+1,j))+(1-
((((2*i)+1)/i)+2)*Fo)*Ta(i,j);
else (i==m && j==n);
Facultad de minas
Ta(i,j)=Fo*((((2*i)+1)/(i))*Ta(i,j-1)+2*Ta(i-1,j))+(1-
((((2*i)+1)/i)+2)*Fo)*Ta(i,j);
end
end
end
Fo=((k*dt)/(RO*cp*dr));
end
disp ('El tiempo requerido para que el horno alcance la temperatura de

operación es: ' )
disp(time)
Finalmente, el resultado que arroja Matlab es:

El tiempo requerido para que el horno alcance la temperatura de operación es: 8120
segundos
Lo que equivale a 2.25 horas
¿Qué hacer para reducir el tiempo sin que el refractario tenga problemas?
Aunque el tiempo de operación para este tipo de condiciones que se propusieron es
bastante bueno, se podría hacer una mejora en la reducción del tiempo de operación si en
el flujo de entrada se agregaran partículas de sílice debido a que tiene propiedades de
conductividad térmica alta además de tener estabilidad mecánica proporciona que se
reduzcan los tiempos de operación del fluido. Es importante señalar que la sílice también
es un gran material refractario por lo que permite el incremento de la producción al acortarse
los tiempos de calentamiento.

Transferencia de Calor

Cargado por

Copyright:

Formatos disponibles

Transferencia de Calor

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Transferencia de Calor

Cargado por

Copyright:

Formatos disponibles

Operaciones de transferencia de Calor 3007820

Departamento de procesos y energía

Tarea N° 2 de Transferencia de Calor

 La geometría del horno de operación se ha tomado de forma cilíndrica y se

A continuación, se presentan las siguientes ecuaciones por el método de diferencias finitas,

 2≤𝑖 ≤𝑚−1𝑦2≤𝑗 ≤𝑛−1

El anterior sistema de ecuaciones se programó en Matlab, con el fin de determinar el tiempo

A continuación, se presenta el código que se utilizó en Matlab.

disp ('El tiempo requerido para que el horno alcance la temperatura de

Finalmente, el resultado que arroja Matlab es:

También podría gustarte