UNIVERSIDAD NACIONAL DE

FACULTAD DE INGENIERA
ESCUELA ACADMICO PROFESIONAL DE
INGENIERA MECNICA
PRESENTADO POR:
RODRIGUEZ PEREZ, LUIS

PROBLEMA

CURSO:
TRANSFERENCIA DE CALOR

RADIACION

UNIDAD:
tr

III
DOCENTE:
ING. ELI GUAYAN HUACCHA
CICLO:
VII

TRUJILO PERU
2014
INDICE
I.

DE

PLANTEAMIENTO DEL PROBLEMA:

Considere una cavidad cilndrica cerrada en la parte inferior con una abertura en la parte
superior.
Para las siguientes condiciones, calcule la transferencia de radiacin a travs de la
abertura de la cavidad cuando la temperatura de los alrededores es 0K. Tambin
determine la emisividad efectiva de la cavidad, e.
(a) Todas las superficies son negras a 620K.
(b) La superficie inferior de la cavidad es difusa y gris con una emisividad de 0.6 a
600K, mientras que todas las superficies interiores son rerradiantes.
(c) Todas las superficies interiores son difusas y grises con una emisividad de 0.6 y
una temperatura uniforme de 620K
(d) Para las configuraciones de cavidad de 20, 40 y 80mm, grafique e, como funcin
de la emisividad de la superficie interior en un margen de 0.6 a 1.0. Todas las otras
condiciones permanecen iguales.

II.

ESQUEMAS DEL PROBLEMA:

III.

HIPOTESIS DE TRABAJO O SUPUESTOS TEORICO:

a) Para el caso (a) todas las superficies son negras.

b) Para el caso (b) la superficie inferior es difusa y gris, las dems reradiantes
(adiabtica).
c) Para el caso (c) todas las superficies son difusas y grises.
d) Las superficies tienen radiosidad uniforme.
e) Los alrededores (medio ambiente) se encuentra a 0K, de esta manera T1=0 y 1=1
IV.

ECUACIONES Y FORMULISMO ANALITICO:

El factor de visin para discos paralelos coaxiales:

Emisividad efectiva de la cavidad:

Un
Se refiere a la tasa de la potencia radiante dejando la cavidad para que a partir de
un cuerpo negro que tiene la misma rea de la abertura de la cavidad y a la
temperatura de la superficie interior de la cavidad es decir:

e=

1
Q
A1 T 4

Para una superficie negra, la transferencia neta de calor por radiacin desde
cualquier superficie i se determina:

Para una superficie es difusa y gris, la transferencia neta de calor por radiacin es:

V.

ANALISIS Y SOLUCIN DEL PROBLEMA:

Clculo de las reas:

2

A 1=

D1
=2.54469 x 104 m2
4

A 2=

D 22
4 2
=10.1788 x 10 m
4

A 3= D2 L=56.5487 x 104 m2

A 4 = A2 A 1=7.63411 x 104 m2

Clculo de los Factores de visin:

De las condiciones del problema sabemos que T1=0K, entonces el Eb1=J1, la
transferencia de calor por radiacin que pasa a travs de la superficie 1 (imaginaria) ser
debido al intercambio de calor por radiacin entre las otras superficies de la cavidad.
Los factores de visin sern (por simple inspeccin):

F1-1 = 0
F1-4 = 0
F2-2 = 0
F4-4 = 0

El clculo de F2-1 lo hacemos aplicando las frmulas para discos paralelos coaxiales:

R2=

36
=0.72
50

R1=

18
=0.36
50

S=1+

Tenemos:

1+0.362
=3.179
0.722

{ [

( ) ] }=0.0807

1
0.36
F21= 3.179 3.17924
2
0.72

2 1
2

Para el clculo de F2-(1,4) se utiliza la misma metodologa:

36
=0.72
50

R2=

R1,4 =

S=1+

36
=0.72
50
1+0.722
=3.929
0.722

Tenemos:

{ [

Por la regla de la superposicin:

F2-(1,4) = F2-1 + F2-4
F2-4 = F2-(1,4) - F2-1
Tenemos:

F2-4 = 0.1929

Por la regla de reciprocidad:

A2 F2-1 = A1 F1-2

F12=

A2
x F 21
A1

Tenemos:

F1-2 = 0.3228

Por la regla de la suma:

F1-1 + F1-2 + F1-3 + F1-4 = 1
F1-3 = 1 (F1-1 + F1-2 + F1-4)
Tenemos:

( ) ] }=0.2736

1
0.72
F2(1,4 )= 3.929 3.92924
2
0.72

F1-3 = 0.6772

2 1
2

Por la regla de reciprocidad:

A3 F3-1 = A1 F1-3

F31=

A1
x F13
A3

Tenemos:

F3-1 = 0.03047

Para el clculo de F2-3 aplicamos la regla de la suma:

F2-1 + F2-2 + F2-3 + F2-4 = 1
F2-3 = 1 (F2-1 + F2-2 + F2-4)
Tenemos:

F2-3 = 0.7264

Por la regla de reciprocidad:

A2 F2-3 = A3 F3-2

F12=

A2
x F23
A3

Tenemos:

F3-2 = 0.13075

Para el clculo de F4-2 aplicamos la regla de reciprocidad:

A2 F2-4 = A4 F4-2

F 42=

A2
x F24
A4

Tenemos:

F4-2 = 0.2572

Para el clculo de F4-3 aplicamos la regla de la suma:

F4-1 + F4-2 + F4-3 + F4-4 = 1
F4-3 = 1 (F4-1 + F4-2 + F4-4)
Tenemos:

F4-3 = 0.7428

Para el clculo de F3-4 aplicamos la regla de reciprocidad:

A3 F3-4 = A4 F4-3

F3 4=

A4
x F 43
A3

Tenemos:

F3-4 = 0.1003

Por ultimo hallamos el F3-3 con la regla de la suma:

F3-1 + F3-2 + F3-3 + F3-4 = 1
F3-3 = 1 (F3-1 + F3-2 + F3-4)
Tenemos:

F3-3 = 0.73848

Por tanto la matriz de los factores de visin Fij (1i4 y 1j4) nos queda:

0
0.3228 0.6772
0
0.0807
0
0.7264 0.1929
0.03047 0.13075 0.73848 0.1003
0
0.2572 0.7428
0

Solucin de la Parte (a):

Las superficies internas de la cavidad son negras a 620K, la tasa de calor que sale de A1:
Aplicando la ecuacin:

La transferencia de calor neta por radiacin que sale por la superficie 1:

1= A1 (T 41 T 4s )( F 11+ F 12 + F13 + F 1 4)
Q

 1=( 2.54469 x 104 )(5.67 x 108 )(04 620 4)(1)

Q
1=2.132W
Q
La emisividad efectiva de acuerdo a la definicin, queda:

1=1

Solucin de la Parte (b):

Para este caso tenemos que la superficie inferior de la cavidad es difusa y gris a 620K con
una emisividad

=0.6 , las dems son rerradiantes. El anlisis para para obtener

1
Q
en la cavidad lo podemos hacer con las ecuaciones planteadas en la parte IV para
superficies difusas y grises.

Para el rea A2:

E b 2J 2 J 2J 1
J J
J J
=
+ 2 3 + 2 4
1 2
1
1
1
A2 F21 A 2 F 23 A 2 F 24
2 A 2
Reemplazando:

8378.18J 2
J 2 0
J 2J 3
J 2J 4
=
+
+
10.6
1
1
1
4
4
4
0.6(10.1788 x 10 ) (10.1788 x 10 ) 0.0807 (10.1788 x 10 )0.7264 (10.1788 x 104 ) 0.1929

0.1929 J
( 20.1929 J 4 )
12567.271.5 J 2=0.0807 J 2+(0.7264 J 20.7264 J 3)+

12567.27=2.5 J 20.7264 J 30.1929 J 4

Para el rea A3:

0=

J 3J 1
J J
J J
+ 3 2 + 3 4
1
1
1
A 3 F 31 A3 F3 2 A 3 F 34

Reemplazando:

0=

J 30
J 3J 2
J 3J 4
+
+
1
1
1
4
4
(56.5487 x 10 )0.03047 (56.5487 x 10 )0.13075 (56.5487 x 104 )0.1003

0.13075 J
( 30.13075 J 2)+(0.1003 J 3 0.1003 J 4 )
0=0.03047 J 3+

0=0.26152 J 3 0.13075 J 20.1003 J 4

Para el rea A4:

0=

J 4J 1
J J
J J
+ 4 2 + 4 3
1
1
1
A 4 F 41 A 4 F 4 2 A 4 F 43

Reemplazando:

0=

J 4 0
J 4 J 2
J 4 J 3
+
+
1
1
1
4
4
(7.63411 x 10 )0 (7.63411 x 10 )0.2572 (7.63411 x 104) 0.7428

0=0+(0.2572 J 4 0.2572 J 2)+(0.7428 J 40.7428 J 3 )

0=J 40.2572 J 20.7428 J 3

Obtenemos por tanto un sistema de ecuaciones:

12567.27=2.5 J 20.7264 J 3 0.1929 J 4

0=0.26152 J 30.13075 J 2 0.1003 J 4
0=J 40.2572 J 20.7428 J 3

La solucin ser:

J2=7296.39 W/m2
J3=6107.62 W/m2
J4=6413.37 W/m2

Por tanto la transferencia de calor neta por radiacin que sale por la superficie 1:

= J 1J 2 + J 1J 3 + J 1J 4
Q
1
1
1
1
A1 F12 A 1 F 13 A 1 F 14
Reemplazando:

1=
Q

7296.39
6107.62
6413.37
+
+
1
1
1
4
4
2.54469 x 10 (0.3228) 2.54469 x 10 (0.6772) 2.54469 x 104 (0)

1=1.6518W
Q

La emisividad efectiva:

e=

1
Q
1.6518
=
=0.7748
4
4
A 1 T (2.54469 x 10 )(5.67 x 108 )(6204 )

Solucin de la Parte (c):

Para este caso tenemos todas las superficies son difusas y grises con

=0.6 y

temperatura 620K. Se procede a resolver de la misma manera que el caso (b):

Tenemos:
4

Eb 2=E b 3=Eb 4= T =5.67 x 10 ( 620 ) =8378.18W /m

Para el rea A2:

E b 2J 2 J 2J 1
J J
J J
=
+ 2 3 + 2 4
1 2
1
1
1
A2 F21 A 2 F 23 A 2 F 24
2 A 2
Reemplazando:

8378.18J 2
J 2 0
J 2J 3
J 2J 4
=
+
+
10.6
1
1
1
4
4
4
4
0.6(10.1788 x 10 ) (10.1788 x 10 ) 0.0807 (10.1788 x 10 )0.7264 (10.1788 x 10 ) 0.1929

0.1929 J
( 20.1929 J 4 )
12567.271.5 J 2=0.0807 J 2+(0.7264 J 20.7264 J 3)+

12567.27=2.5 J 20.7264 J 30.1929 J 4

Para el rea A3:

E b 3J 3 J 3J 1
J 3J 2
J 3 J 4
=
+
+
1 3
1
1
1
A3 F 31 A 3 F 32 A3 F3 4
3 A 3
Reemplazando:

8378.18J 3
J 30
J 3J 2
J 3J 4
=
+
+
10.6
1
1
1
4
4
4
0.6(56.5487 x 10 ) (56.5487 x 10 ) 0.03047 (56.5487 x 10 )0.13075 (56.5487 x 104) 0.1003

0.13075 J
( 30.13075 J 2)+(0.1003 J 3 0.1003 J 4 )
12567.271.5 J 3=0.03047 J 3 +

12567.27=1.76152 J 30.13075 J 20.1003 J 4

Para el rea A4:

E b 4 J 4 J 4 J 1
J 4J 2
J 4J 3
=
+
+
1 4
1
1
1
A 4 F 4 1 A 4 F 42 A 4 F 43
4 A 4

Reemplazando:

8378.18J 4
J 4 0
J 4 J 2
J 4 J 3
=
+
+
10.6
1
1
1
4
4
4
0.6(7.63411 x 10 ) (7.63411 x 10 ) 0 (7.63411 x 10 )0.2572 (7.63411 x 104) 0.7428

12567.271.5 J 4 =0+(0.2572 J 40.2572 J 2 )+(0.7428 J 4 0.7428 J 3)

12567.27=2.5 J 40.2572 J 20.7428 J 3

Obtenemos por tanto un sistema de ecuaciones:

12567.27=2.5 J 20.7264 J 30.1929 J 4

12567.27=1.76152 J 30.13075 J 20.1003 J 4
12567.27=2.5 J 4 0.2572 J 2 0.7428 J 3

La solucin ser:

J2=8050.55 W/m2
J3=8204.07 W/m2
J4=8292.74 W/m2

Por tanto la transferencia de calor neta por radiacin que sale por la superficie 1:

= J 1J 2 + J 1J 3 + J 1J 4
Q
1
1
1
1
A1 F12 A 1 F 13 A 1 F 14
Reemplazando:

1=
Q

8050.55
8204.07
8292.74
+
+
1
1
1
4
4
2.54469 x 10 (0.3228) 2.54469 x 10 (0.6772) 2.54469 x 104 (0)

1=2.0751W
Q
La emisividad efectiva:

e=

1
Q
2.0751
=
=0.9733
4
4
A 1 T (2.54469 x 10 )(5.67 x 108 )(6204 )

VI.

PROGRAMACIN EN MATLAB Y GRAFICAS:

a) PROGRAMA PARA RESOLVER EL PROBLEMA (Partes: a,b,c):

%% CAVIDAD CILINDRICA CERRADA CON ABERTURA AL MEDIO AMBIENTE
% Solucin de las partes (a)(b)(c)
format short g
% DATOS DEL PROBLEMA
Ta = 0
; % TEMP ALREDEDORES [K]
d1 = 18e-3 ; % [m]
d2 = 36e-3 ; %[m]
d4 = d2
; %[m]
L = 50e-3 ; %[m]
sigm = 5.67e-8; % [W/m2.K]

Cte Stefan-Boltzmann

r1 = d1/2;
r2 = d2/2;
r4 = d4/2;
% Clculamos el rea para cada superficie:
A1 = pi/4*d1^2
; % [m^2]
A2 = pi/4*d2^2
; % [m^2]
A3 = pi*d2*L
; %[m^2]
A4 = pi/4*d4^2 - A1; % [m^2]

%
%
%
%

CALCULO DEL NUMERO DE FACTORES DE FORMA:

N=nmero de superficies, N = 4 ;
N(N-1)/2 = 6 factores de forma
F11 + F12 + F13 + F14 = 1

% F21 + F22 + F23 + F24 = 1

% F31 + F32 + F33 + F34 = 1
% F41 + F42 + F43 + F44 = 1
% FACTOR DE FORMA PARA DISCOS COAXIALES: Fab
= Disco2Disco(ra,rb,L)
% Donde por simple inspeccin
F11
= 0
;
F14
= 0
;
F41
= 0
;
F22
= 0
;
F44
= 0
;
F21
= Disco2Disco(r2,r1,L);
F12
= F21 * A2/A1
;
F2_14 = Disco2Disco(r2,r4,L);
F24
= F2_14 - F21
; % F(2->1+4) = F(2->1) + F(2->4)
superposicin
F42
= F24 * A2/A4
;
F23
= 1 -(F21 + F22 + F24);
F32
= F23 * A2/A3
;
F43
= 1 -(F41 + F42 + F44);
F34
= F43 * A4/A3
;
F13
= 1 -(F11 + F12 + F14);
F31
= F13 * A1/A3
;
F33
= 1 -(F31 + F32 + F34);
% Matriz factores de forma
F = [F11 F12 F13 F14; F21 F22 F23 F24 ;...
F31 F32 F33 F34 ; F41 F42 F43 F44];
% emisividades
%
(a)
e = [ 1
1
1
1

(b)
1
0.6
NaN
NaN

(c)
1
; ...
0.6 ; ...
0.6 ; ...
0.6 ];

% A1
% A2
% A3
% A4

[n,m]=size(e);
% Temperaturas de las superficies
%
A1
A2
A3
A4
T = [0
620
620
620] ; % [K]
% Vector de areas
A = [A1 A2 A3 A4];
K=zeros(n);
C=zeros(1,n);
for p=1:m
for i=1:n
% Factor Ci = (1-ei)/ei ;
C(i) = (1 - e(i,p))/e(i,p) ;
% Factor Kii = 1 + Ci
K(i,i) = 1 + C(i) ;
I = eye(n);
if isnan(C(i))==1

for j=1:n
M(i,j)= I(i,j)-F(i,j) ;
end
B(i,1)= 0 ;

else

% Factor M(i,j)= -Ci*Fij + Kij

for j=1:n
M(i,j)= -C(i)*F(i,j) + K(i,j) ;
end
B(i,1)= sigm*T(i).^4 ;

end
end
% Columna de Constantes sigm*Ti^4: B
% Solucin: [J1 J2 J3 J4]
J = M\B;
% Colecin de soluciones
SJ(p,:) = J' ;
% Para todos los casos la emisividad efectiva: Ee = -q1/(A1*sigm*T^4)
% donde q1 = sum(A1 * F1j *(J1 - Jj))
q1 = 0 ;
for j=1:n
q1 = q1 + A1 * F(1,j) * (J(1)-J(j));
end
% Colecin de q1
Q1(p) = q1 ;
Ee(p) = -q1/(A1*sigm*T(2)^4);
end

fprintf('---------------------------------------------------------------\n');
fprintf('
J1
J2
J3
J4
Q1
Ee\n');
fprintf('----------------------------------------------------------------\n');
fprintf('(a)
%4.3f
%4.3f
%4.3f
%4.3f
%4.3f
%4.3f
\n',SJ(1,1),SJ(1,2),SJ(1,3),SJ(1,4),Q1(1),Ee(1));
fprintf('(b)
%4.3f
%4.3f
%4.3f
%4.3f
%4.3f
%4.3f
\n',SJ(2,1),SJ(2,2),SJ(2,3),SJ(2,4),Q1(2),Ee(2));
fprintf('(c)
%4.3f
%4.3f
%4.3f
%4.3f
%4.3f
%4.3f
\n',SJ(3,1),SJ(3,2),SJ(3,3),SJ(3,4),Q1(3),Ee(3));

Resultados:

b) PROGRAMA PARA LAS GRAFICAS (Parte d):

%% CAVIDAD CILINDRICA CON ABERTURA EN LA PARTE SUPERIOR
%Se variar la profundidad de la cavidad
%Parte (d) del ejercicio:

format short g
% DATOS DEL PROBLEMA
Ta = 0
; % Temperatura de alrededores [K]
d1 = 18e-3 ; % [m]
d2 = 36e-3 ; % [m]
d4 = d2
; % [m]
L = [20 40 80]*1e-3; % Valores de profundidad [m]
Ln=length(L);
E = 0.6:0.01:1; % Vector de emisividad;
Le = length(E);
sigm = 5.67e-8; % [W/m2.K] Cte Stefan-Boltzmann
r1 = d1/2;
r2 = d2/2;
r4 = d4/2;
% Clculo de el rea de cada una de las superficies
for l=1:Ln
for t=1:Le
A1 = pi/4*d1^2
; % [m^2]
A2 = pi/4*d2^2
; % [m^2]
A3 = pi*d2*L(l)
; % [m^2]
A4 = pi/4*d4^2 - A1; % [m^2]
% CALCULO DEL NUMERO DE FACTORES DE FORMA
%N=numero de superficies; N=4
% N(N-1)/2 = 6
% F11 + F12 + F13 + F14 = 1
% F21 + F22 + F23 + F24 = 1
% F31 + F32 + F33 + F34 = 1
% F41 + F42 + F43 + F44 = 1
% FACTOR DE FORMA ENTRE DISCOS COAXIALES: Fab
=
Disco2Disco(ra,rb,L)
% Por simple inspeccion:
F11
= 0
;
F14
= 0
;
F41
= 0
;
F22
= 0
;
F44
= 0
;
F21
= Disco2Disco(r2,r1,L(l));
F12
= F21 * A2/A1
;
F2_14 = Disco2Disco(r2,r4,L(l));
F24
= F2_14 - F21
; % F(2->1+4) = F(2->1) + F(2->4)
superposicin
F42
= F24 * A2/A4
;
F23
= 1 -(F21 + F22 + F24);
F32
= F23 * A2/A3
;
F43
= 1 -(F41 + F42 + F44);

F34
F13
F31
F33

=
=
=
=

F43 * A4/A3
;
1 -(F11 + F12 + F14);
F13 * A1/A3
;
1 -(F31 + F32 + F34);

% Matriz factores de forma

F = [F11 F12 F13 F14; F21 F22 F23 F24 ;...
F31 F32 F33 F34 ; F41 F42 F43 F44];
% emisividades de las superficies
%
(b)
e = [ 1
E(t)
NaN
NaN

(c)
1
;
E(t) ;
E(t)
;
E(t) ];

...
...
...
...

%
%
%
%

A1
A2
A3
A4

[n,m]=size(e);
% Temperaturas de las superficies
%
A1
A2
A3
A4
T = [0
620
620
620] ; %K
% Vector de areas
A = [A1 A2 A3 A4];
K=zeros(n);
C=zeros(1,n);
for p=1:m
for i=1:n
% Factor Ci = (1-ei)/ei ;
C(i) = (1 - e(i,p))/e(i,p) ;
% Factor Kii = 1 + Ci
K(i,i) = 1 + C(i) ;
I = eye(n);
if isnan(C(i))==1
for j=1:n
M(i,j)= I(i,j)-F(i,j) ;
end
B(i,1)= 0 ;
else
% Factor M(i,j)= -Ci*Fij + Kij
for j=1:n
M(i,j)= -C(i)*F(i,j) + K(i,j) ;
end
B(i,1)= sigm*T(i).^4 ;
end

end
% Columna de Constantes sigm*Ti^4: B
% Solucin: [J1 J2 J3 J4]
J = M\B;
% Colecin de soluciones
% Para todos los casos la emisividad efectiva: Ee = -q1/
(A1*sigm*T^4)
% donde q1 = sum(A1 * F1j *(J1 - Jj))
q1 = 0 ;

for j=1:n
q1 = q1 + A1 * F(1,j) * (J(1)-J(j));
end
% Colecin de q1
Ee(p,t,l) = -q1/(A1*sigm*T(2)^4); % emisividad efectiva
end

end

end
plot(E,Ee(1,:,1),'r',E,Ee(1,:,2),'k',E,Ee(1,:,3),'b')
xlabel('Emisividad de la superficie inferior');
ylabel('Emisividad efctiva')
grid on
figure(2);
plot(E,Ee(2,:,1),'r',E,Ee(2,:,2),'k',E,Ee(2,:,3),'b')
xlabel('Emisividad de las superficies internas');
ylabel('Emisividad efctiva')
grid on

0.95

Emisividad efctiva

0.9

20mm
40mm
80mm

0.85

0.8

0.75

0.7
0.6

0.65

0.7
0.75
0.8
0.85
0.9
Emisividad de la superficie inferior

0.95

1
0.995
0.99

Emisividad efctiva

0.985
0.98
0.975
20mm
40mm
80mm

0.97
0.965
0.96
0.955
0.95
0.6

VII.

0.65

0.7
0.75
0.8
0.85
0.9
Emisividad de las superficies internas

0.95

BIBLIOGRAFIA:

Transferencia de Calor y Masa Cengel 4Ed.

Fundamentos de transferencia de Calor y Masa Incropera 4Ed.