MEC-311: Equipos de Procesos Térmicos TAREA N°2 Semestre I, 2009
Las figuras muestran un diagrama esquemático de un convertidor Pierce-Smith utilizado en la fundición y
refinería de cobre ACME S.A. (se usa para la conversión de metal blanco a cobre blister). El reactor es
esencialmente cilíndrico, con un diámetro interior de 4.5 m y una longitud interior de 13 m. El casco
Tarea Nº3
metálico del convertidor tiene un espesor de 25 mm en todas las paredes, y es de acero con una
conductividad térmica de 50 W/mK. Sobre este caso de acero se encuentra montado el revestimiento
refractario.
Más especificaciones:
MEC278: Equipos de Procesos Térmicos – C. Rosales
manto: ladrillos de alta alúmina, espesor total 380 mm; emisividad a la temperatura típica de
trabajo: 0.45; conductividad térmica 1.8 W/mK.
culata norte: ladrillos de Cr-Mg, espesor total 440 mm; emisividad a la temperatura típica de
La Figuras muestran trabajo:
esquemas simplificados
0.35; conductividad térmica de un reactor para especies gaseosas. El espacio interior
2.8 W/mK.
es esencialmente un volumen cilı́ndrico. El manto cilı́ndrico está dividido en dos zonas. La zona A tiene
culata sur: concreto refractario, espesor total 420 mm; emisividad a la temperatura típica de
paredes calentadas mediante resistencias eléctricas homogéneamente distribuidas. La zona B tiene paredes
trabajo: 0.28; conductividad térmica 1.48 W/mK.
que cuentan con un sistema de enfriamiento, después de la capa de material refractario. En el lado Sur, hay
además una aberturaLa que
plantacircular,
de ACME está coaxial
situada con
a nivelladelcámara cilı́ndrica,
mar, y el interior que conecta
del convertidor hacia
se encuentra un ducto de salida de
a presión
atmosférica. La temperatura ambiente en el lugar es de 15 °C. La emisividad de las superficies exteriores
productos. del casco es 0.7 y existe además enfriamiento por convección natural. El coeficiente de convección exterior
El diámetro interior es en
se estima 3.3[m],
10 W/m2K. mientras que la longitud de la cámara es 9.0[m]. La longitud de la Zona
A es 3.5[m] y el diámetro de laetapas
Durante ciertas abertura de este
del proceso salida
equipoesvacía
0.8[m].
completamente su carga de cobre blister hacia un
horno de refino,
En el instante mostrado y debe permanecer
el reactor está vacı́ovacío (solo
durantehay
variasaire
horastransparente
esperando la siguiente
en carga de metal Para un análisis
el interior).
blanco para procesar. Durante el tiempo que permanece vacío, el convertidor debe mantenerse caliente
simplificado, el interior
(paradel reactor
evitar esfuerzosse puede
térmicos dividir
sobre en las Para
los refractarios). siguientes superficies,
esto se enciende para
un quemador lasnatural
de gas cuales la temperatura
se puede considerar (situado
aproximadamente
en la culata norte), eluniforme:
cual opera con un 5% de exceso de aire y produce gases de combustión cuya
distribución de temperatura se puede suponer homogénea en el interior del reactor.
lado Sur
manto
Sector 2
Sector 1
lado
Norte
Sector 1 Sector 2
3.5 m
9m
Lado Norte: Temperatura: 1080◦ C; emisividad: 0.45
Manto cilı́ndrico en Zona A: Temperatura: 1140◦ C; emisividad: 0.52
B:1 Temperatura: 850◦ C;Sector
Manto cilı́ndrico en ZonaSector emisividad:
2 0.40
Lado Sur: Temperatura: 780◦ C; emisividad: 0.38
Abertura: La temperatura equivalente del medio que se “ve” a través de la abertura es 60◦ C.
1. Determine el calor (kW) transferido desde las paredes calentadas, el calor transferido a las paredes
con enfriamiento y el calor
3.5 mque se pierde por la abertura hacia el ducto.
9m
2. Para mejorar laAppendix
eficiencia,
C se ha hecho algunas modificaciones al reactor, como se ilustra más abajo.. 883
La abertura fue cerrada y en su lugar se colocó un panel radiante (con la misma geometrı́a de la
abertura anterior), que tiene una emisividad 0.6 y que envı́a un flujo neto de radiación de 39[kW/m2 ].
Sector
8 Sector 1 Two 2finite rectangles of same length, having one
Por otra parte, se instaló aislación térmica en el resto de la superficie del lado Sur, de modo que esa
common edge and having an angle of 90° to each
parte se puede considerar adiabática.
ℓ other.
En estas nuevas condiciones las nuevas temperaturas de las h otras superficies
w son:
H= ; W=
A2 l l
Lado Norte: 1045◦ C
h 90° ◦C
Manto cilı́ndrico en Zona
3.5 m A: 1080
A1 1 1 1
F1− 2 = W tan W + H tan H
−1 −1
Manto cilı́ndrico en Zona
w B: 850◦ C 9m πW
Determine las nuevas temperaturas de equilibrio del lado 1
Sur y del nuevo 1 radiante.
panel (1 + W 2 )(1 + Calcule
H 2)
− H 2 + W 2 tan −1 + ln
además el calor transferido a la pared con enfriamiento que quedó. H 2 + W 2 4 1 + W + H 22
W 2 (1 + W 2 + H 2 )
W2
H 2 (1 + H 2 + W 2 )
H2
× 2
× 2
(11 + W )(W + H ) (1 + H )( H + W )
2 2 2 2
Sector 1 Sector 2
9 Plane element dA1 to circular disk in plane parallel to
element; normal to element passes through center of
disk.
dA1 r2
Fd1− 2 = 2 2
h +r
3.5 m
9m
h
A2
Indicaciones: Para calcular los factores de forma entre las 5 superficies utilice las siguientes fórmulas, y
álgebra de factores de forma: r
10
Factor de forma entre dos discos paralelos con centrsoParallel
sobre circular
la mismadisksnormal:
with centers along the same
normal.
A2 r r
R1 = 1 ; R2 = 2
h h
r2
1 + R22
X = 1+
R12
h 2
1 R
F1− 2 = X − X2 − 4 2
2 R1
r1
A1
Factor de forma desde el manto de un cilindro a si mismo:
1 p 2
F1−1 = 1 − 4z + 1 − 1
2z
donde K10138.indb 883 r 8/18/10 10:45:27 AM
z=
L
con r el radio del cilindro y L su longitud.
clc
clear all
r1 = 1.65;
r2 = 0.4;
L1 = 3.5;
L2 = 9;
A1 = pi* r1^2;
A2 = 2*pi*r1*L1;
A3 = 2*pi*r1*(L2-L1);
A5 = pi* r2^2;
A6 = A1;
A7 = A1;
A4 = A7-A5;
A8 = 2*pi*r1*L2;
R1 = r1/L2;
R2 = r1/L2;
X = 1 + (1+ R2^2)/(R1^2);
F11 = 0;
F17 = (X-sqrt(X^2 -4*(R2/R1)^2))/2;
F18 = 1-F17;
%%%%%%
R1 = r1/L1;
R2 = r1/L1;
X = 1 + (1+ R2^2)/(R1^2);
F16 = (X-sqrt(X^2 -4*(R2/R1)^2))/2;
F12 = 1-F16;
F13 = F18-F12;
%%%%%%
R1 = r1/L2;
R2 = r2/L2;
X = 1 + (1+ R2^2)/(R1^2);
F15 = (X-sqrt(X^2 -4*(R2/R1)^2))/2;
F14 = F17-F15;
%%%%%%%%
R1 = r1/(L2-L1);
R2 = r1/(L2-L1);
X = 1 + (1+ R2^2)/(R1^2);
1
F67 = (X-sqrt(X^2 -4*(R2/R1)^2))/2;
F63 = 1-F67;
R1 = r1/(L2-L1);
R2 = r2/(L2-L1);
X = 1 + (1+ R2^2)/(R1^2);
F65 = (X-sqrt(X^2 -4*(R2/R1)^2))/2;
F64 = F67-F65;
r = r1;
L = (L2-L1);
z = r/L;
F33 = 1- (sqrt(4*z^2 +1)-1)/(2*z);
F46 = F64*(A6/A4);
F43 = 1-F46;
F34 = F43*(A4/A3);
F56 = F65*(A6/A5);
F53 =1-F56;
F35 = F53*(A5/A3);
%%%%%%
F51 = F15*(A1/A5);
F41 = F14*(A1/A4);
F48 = 1-F41;
F58 = 1-F51;
F42 = F48-F43;
F52 = F58-F53;
F31 = F13*(A1/A3);
F32 = 1-F31-F33-F34-F35;
%%%%
r = r1;
L = L1;
z = r/L;
F22 = 1- (sqrt(4*z^2 +1)-1)/(2*z);
F21 = F12*(A1/A2);
F23 = F32*(A3/A2);
F24 = F42*(A4/A2);
F25 = F52*(A5/A2);
F44 =0;
F45 =0;
F54 = 0;
F55 = 0;
2
F = [F11 F12 F13 F14 F15;...
F21 F22 F23 F24 F25;...
F31 F32 F33 F34 F35;...
F41 F42 F43 F44 F45;...
F51 F52 F53 F54 F55];
A = [A1;A2;A3;A4;A5];
currentFolder = pwd;
filePath = fullfile(currentFolder, 'T3 Termica', 'FyA.mat');
save(filePath, 'F', 'A');
Published with MATLAB® R2024b
3
n=5;
currentFolder = pwd;
filePath = fullfile(currentFolder, 'T3 Termica', 'FyA.mat');
load(filePath);
A;
F;
C = zeros(n, n);
e = [0.45 0.52 0.4 0.38 1] ;
T = [1353.15 1413.15 1123.15 1053.15 333.15];
sigma = 5.669*10^(-8);
for j = 1:n
for k = 1:n
if j == k
delta = 1;
else
delta = 0;
end
C(j, k) = (1 / e(k)) * (delta - (1 - e(k)) * F(j, k));
end
end
E = zeros(5,1);
for j = 1:n
E(j)=sigma*e(j)*T(j)^4;
end
b=zeros(5,1);
for j = 1:n
sum=0;
for k = 1:n
sum = sum + F(j,k)*E(k);
end
b(j) = E(j)-sum;
end
q = C\b;
Q = zeros(5,1);
for j = 1:n
Q(j) = q(j)*A(j);
end
1
Q
Comprobacion = Q(1)+Q(2)+Q(3)+Q(4)+Q(5);
disp("El Q transferido desde el lado Norte y el Manto A: ")
disp(Q(1)+Q(2)+ "[W]")
disp("El Q transferido al lado Sur y al Manto B: ")
disp(Q(3)+Q(4)+ "[W]")
disp("El Q que se pierde por la abertura: ")
disp(Q(5)+ "[W]")
Q =
-36233.0655805435
519771.812597762
-385839.025170773
-72632.4340895746
-25067.2877568707
El Q transferido desde el lado Norte y el Manto A:
483538.747[W]
El Q transferido al lado Sur y al Manto B:
-458471.4593[W]
El Q que se pierde por la abertura:
-25067.2878[W]
Published with MATLAB® R2024b
2
n=5;
currentFolder = pwd;
filePath = fullfile(currentFolder, 'T3 Termica', 'FyA.mat');
load(filePath);
A;
F;
q = [1 1 1 0 39000];
C = zeros(n, n);
e = [0.45 0.52 0.4 0.38 0.6] ;
T = [1318.15 1353.15 1123.15];
sigma = 5.669*10^(-8);
for j = 1:n
for k = 1:n
if j == k
delta = 1;
else
delta = 0;
end
C(j, k) = (1 / e(k)) * (delta - (1 - e(k)) * F(j, k));
end
end
E = zeros(3,1);
for j = 1:3
E(j)=sigma*e(j)*T(j)^4;
end
%Ahora se empezará a plantear el sistema de ecuaciones%
d = zeros(5,1);
d(1) = E(1) - C(1,5)*q(5) - F(1,2)*E(2) - F(1,3)*E(3);
d(2) = E(2) - C(2,5)*q(5) - F(2,1)*E(1) - F(2,2)*E(2) - F(2,3)*E(3);
d(3) = E(3) - C(3,5)*q(5) - F(3,1)*E(1) - F(3,2)*E(2) - F(3,3)*E(3);
d(4) = -C(4,5)*q(5) - F(4,1)*E(1) - F(4,2)*E(2) - F(4,3)*E(3);
d(5) = -C(5,5)*q(5) - F(5,1)*E(1) - F(5,2)*E(2) - F(5,3)*E(3);
z = [C(1,1) C(1,2) C(1,3) F(1,4) F(1,5);...
C(2,1) C(2,2) C(2,3) F(2,4) F(2,5);...
C(3,1) C(3,2) C(3,3) F(3,4) F(3,5);...
C(4,1) C(4,2) C(4,3) -1 0;...
C(5,1) C(5,2) C(5,3) 0 -1];
x = z\d;
1
%El vector x contiene q1, q2, q3, Eb4 y Eb5 respectivamente
T4 = (x(4)/(sigma*e(4)))^(1/4);
T5 = (x(5)/(sigma*e(5)))^(1/4);
Q = zeros(5,1);
for j = 1:3
Q(j) = x(j)*A(j);
end
Q(4) = 0;
Q(5) = 39000*A(5);
disp("La nueva temperatura de equilibrio del lado sur es: ")
disp(T4+ "[K]")
disp("La temperatura de equilibrio del panel radiante es: ")
disp(T5+ "[K]")
disp("El calor transferido al manto B, que es la única pared con
enfriamiento: ")
disp(Q(3)+ "[W]")
Q =
-21492.2620452246
376139.103495198
-374250.379608373
0
19603.5381584003
La nueva temperatura de equilibrio del lado sur es:
1227.7689[K]
La temperatura de equilibrio del panel radiante es:
1353.6181[K]
El calor transferido al manto B, que es la única pared con enfriamiento:
-374250.3796[W]
