Transferencia de Calor por Radiacion

TRANSFERENCIA DE CALOR  FUNDAMENTACION TEORICA: “TRANSFERENCIA DE CALOR POR RADIACION” Todos los cuerpos, cualquiera se

Views 52 Downloads 0 File size 2MB

Report DMCA / Copyright

DOWNLOAD FILE

Recommend stories

Citation preview

TRANSFERENCIA DE CALOR  FUNDAMENTACION TEORICA: “TRANSFERENCIA DE CALOR POR RADIACION” Todos los cuerpos, cualquiera sea su temperatura, emiten energía en forma continua desde sus superficies. Esta energía se denomina energía radiante y es transportada por ondas electromagnéticas, por este motivo, la energía radiante puede transmitirse aún en el vacío. La emisión continua de energía radiante por un cuerpo se denomina radiación. Como consecuencia de este fenómeno, dos cuerpos colocados en el vacío que están a diferentes temperaturas alcanzan el equilibrio térmico debido a que el de menor temperatura recibe energía radiante del otro cuerpo de mayor temperatura. Cuando la energía radiante es absorbida por un cuerpo, se transforma en calor; no obstante la energía radiante también puede ser reflejada (difundida) o refractada (propagada) por los cuerpos. Trataremos únicamente la energía radiante emitida por los sólidos y los líquidos, pues la emitida por los gases obedece a leyes muy diferentes. Hemos dicho que la energía radiante se transmite por ondas electromagnéticas, por lo tanto su velocidad de propagación será la de la luz (300.000 km/seg en el vacío). Las ondas electromagnéticas comprenden: radio ondas, ondas infrarrojas, luz visible, ondas ultravioletas y rayos X y γ; todas diferentes solamente en sus longitudes de ondas. Los cuerpos sólidos y líquidos emiten energía radiante que contiene ondas de todas las frecuencias, cuyas amplitudes dependen principalmente de la temperatura del cuerpo emisor y no del tipo de moléculas que lo formen. Si la radiación emitida por un cuerpo se hace incidir sobre un prisma, se descompone en radiaciones monocromáticas cuyo conjunto se denomina “espectro”. Cada radiación monocromática corresponde a una determinada longitud de onda λ, que está relacionada con la velocidad de propagación c por la ecuación: λ = c . T. D0onde T es el período correspondiente al fenómeno periódico al cual responde la radiación. Por otra parte, T = 1 / υ; siendo υ la frecuencia. El espectro se puede dividir en tres zonas: 1. zona infrarroja: constituida por radiaciones de longitud de onda superiores a 0,8 μ. 2. zona luminosa o visible, cuyas radiaciones poseen longitudes de onda comprendidas entre 0,4 y 0,8 μ. e impresionan la retina humana. 3. zona ultravioleta, cuyas longitudes de onda son inferiores a 0,4 μ. La energía radiante es emitida por toda la materia del cuerpo, pero en general, en su interior la energía emitida por cada punto es nuevamente absorbida por eso solamente se libera la energía correspondiente a una delgada capa de la superficie del cuerpo. No solo depende de la temperatura de la superficie sino también de su naturaleza. Distribución espectral de la energía radiante INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 1

TRANSFERENCIA DE CALOR Lumer y Pringssheim, efectuaron una serie de experimentos en los cuales tomaban las radiaciones emitidas a una cierta temperatura y medían su energía a distintas longitudes de onda. Así encontraron que la energía en las distintas longitudes de onda no era uniforme. Si E λ es la energía emitida con longitud de onda λ, la energía total a temperatura T está dada por:

Por lo tanto, el área encerrada por cada curva representa la energía total emitida a esa temperatura, será proporcional a la cantidad de calor transmitida por unidad de superficie y unidad de tiempo. Se puede observar que a temperaturas bajas, la energía emitida corresponde a radiaciones de longitudes de onda ubicadas en la zona infrarroja. A medida que el cuerpo aumenta su temperatura, emite radiaciones de longitud de onda cada vez menores, alcanzando la zona roja de luz visible y posteriormente al cubrir todo el espectro visible, la luz blanca. Por este motivo, los cuerpos a temperaturas elevadas presentan color rojo y también blanco. Radiación incandescencia Hemos visto anteriormente que la energía emitida por un cuerpo depende de su temperatura. La energía radiante recibida por un cuerpo, en general puede dividirse en tres partes: a) la energía transmitida o programada por el cuerpo sin absorberla; b) la energía reflejada o difundida según las leyes de la óptica y c) la energía que el cuerpo absorbe La cantidad de energía transmitida, reflejada o absorbida por un cuerpo, depende de la naturaleza del material, de la superficie y de la longitud de onda de la radiación. En INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 2

TRANSFERENCIA DE CALOR realidad no existen cuerpos totalmente permeables o impermeables. Por ejemplo, el vidrio es permeable a las radiaciones visibles pero absorbe las infrarrojas. Podemos imaginar la existencia de un cuerpo que sea absolutamente absorbente o sea un cuerpo que absorbiera todas las radiaciones que recibe. Un cuerpo teórico que cumple esta condición, se denomina cuerpo negro. Un cuerpo negro, se puede lograr casi perfectamente practicando un orificio pequeño, de superficie ΔS, en un recinto cerrado, opaco o recubierto de negro de humo, y mantenido a temperatura constante. La radiación absorbida o emitida por dicho sistema, es equivalente a la que correspondería a un cuerpo negro de superficie ΔS, a la misma temperatura. A unos 500 °C, la radiación que emite u cuerpo negro, comienza a tener radiaciones luminosas (rojo cereza). Midiendo la energía de dicha radiación, se puede medir la temperatura del cuerpo, procedimiento en que se basan los métodos ópticos de medición de temperatura en los hornos industriales. Poder emisivo y poder absorbente El poder emisivo o de emisión E, de un cuerpo, se define como la cantidad de calor emitida por unidad de superficie y por unidad de tiempo, en una dirección dad. El valor de E depende fundamentalmente del valor de λ y de T. En general se expresa relacionándola con el poder emisivo del cuerpo negro ideal. Se denomina coeficiente de emisividad e a la relación entre el poder emisivo del cuerpo E, y el poder emisivo del cuerpo negro ES, en iguales condiciones. O sea:

En la expresión anterior vemos que e debe ser un número, independiente de las unidades en que se mida el poder emisivo y cuyo valor está comprendido entre 0 y 1. Por ejemplo, entre 20 y 200°C, los valores aproximados de e son: MATERIAL

e

metal pulido

0,04 – 0,05

metal oxidado

0,80 - 0,90

madera lisa

0,80 - 0,90

material de construcción

0,90

vidrio liso

0,94

negro de humo

0,98

El poder absorbente o de absorción A, de un cuerpo, se define como la cantidad de calor absorbida por unidad de superficie y por unidad de tiempo. Su valor depende de λ y T. Se INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 3

TRANSFERENCIA DE CALOR denomina coeficiente de absorción a, a la relación entre el poder absorbente del cuerpo A y el poder absorbente correspondiente al cuerpo negro en las mismas condiciones AS.

Valor comprendido entre 0 y 1 Se deduce entonces que tanto aS correspondientes al cuerpo negro ideal, deben valer 1. Ley de KIRCHOFF Esta ley establece que la relación entre el poder emisivo y el coeficiente de absorción, es una constante para todas las superficies a valores de λ y de T dados. Si llamamos E1 y E2 a los poder emisivos de dos cuerpos cuyos coeficientes de absorción son a1 y a2, se deberá cumplir que: E1 = E2 también igual a

para el cuerpo negro

Como para el cuerpo negro, aS = 1; entonces el valor de la constante es igual a ES o sea el poder emisivo del cuerpo negro en las condiciones de λ y de T dadas.

Vemos entonces que todo cuerpo puede emitir radiación a una λ y T dadas, según el valor de su coeficiente de absorción. La radiación será mayor cuanto mayor sea el valor de a, en consecuencia el cuerpo negro es el que mayor radiación emite en tales condiciones, pues su valor de a es máximo, igual a 1. Según la ley de Kirchoff, para un cuerpo cualquiera, cuyo poder emisivo es E y su coeficiente de absorción es a, se debe cumplir que: Donde ES = poder emisivo del cuerpo negro Pero, según vimos antes, E = e; donde e es el coeficiente de emisividad del cuerpo, en consecuencia: Donde e = a Se deduce que para toda superficie, el coeficiente de emisividad es igual al coeficiente de absorción. Por tanto, si un cuerpo puede emitir una radiación λ a temperatura T, el mismo cuerpo es también capaz de absorberla en las mismas condiciones. Este fenómeno se conoce como inversión del espectro.

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 4

TRANSFERENCIA DE CALOR Resumiendo lo dicho, se puede establecer que la cantidad de calor transmitida por radiación y por unidad de tiempo, depende no solamente de la temperatura y de la naturaleza de la superficie del cuerpo sino también de las temperaturas y naturaleza de las superficies de los cuerpos circundantes. Cuerpos grises Son aquellos en los cuales el valor del coeficiente de emisividad e, permanece constante para todas las longitudes de onda y temperaturas. Como vimos que e = a, el coeficiente de absorción también debe ser constante. En la práctica no existen cuerpos grises, pues el valor de e no se mantiene constante, sin embargo, en la mayoría de los casos pueden considerarse grises a los cuerpos sin mucho error. El poder emisivo de un cuerpo gris será: E = e. ES Esta ecuación se considera válida para todas las longitudes de onda y en un intervalo dado de temperatura. Ley de Stefan Boltzmann Establece que la cantidad total de calor emitida (en todas las longitudes de onda), por unidad de tiempo y por unidad de superficie del cuerpo negro, es proporcional a la cuarta potencia de la temperatura absoluta del cuerpo. (Según lo visto antes, la cantidad de calor total emitida es proporcional al área encerrada por la curva de radiación: E λ = f(λ). Esta ley se puede expresar matemáticamente: Donde es el coeficiente de radiación total del cuerpo negro, que se puede definir como la radiación integral, para todas las direcciones y longitudes de onda transmitida por unidad de superficie del cuerpo negro, en la unidad de tiempo y por °K de temperatura. Es una constante universal. Par los cuerpos grises podemos aceptar que: δQ = σs S dτ T4 para un cuerpo gris : σ = e σs

δQ = e σs S dτ T4 o

Donde σs = 4,96 x 10-8 Kcal / m2 h ºK4

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 5

TRANSFERENCIA DE CALOR Calor transmitido por radiación Supongamos que un cuerpo 1 a temperatura T1 y de superficie S1 transmite calor a otro cuerpo 2 de temperatura T2 y superficie S2, considerando además que el medio que lo rodea no es absorbente. La cantidad de calor transmitida será igual a la cantidad de calor emitida por el cuerpo 1 a temperatura T1 menos la cantidad de calor reflejada por el cuerpo 2 y menos la cantidad de calor emitida por dicho cuerpo a T2 y absorbida por 1. Si el cuerpo 1 fuera gris y el 2 negro y rodeara totalmente a 1, las cantidades de calor serían: Ejemplo:



Calor emitido por 1



El cuerpo 2 por ser negro no refleja radiación



El cuerpo 1 absorbe de 2

ϕ1 = e1 σs S1 T14 ϕ2 = e1 σs S1 T24

Calor transmitido: ϕ = ϕ1 - ϕ2 = e1 σs S1 ( T14 - T24 ) Si t = T – 273,15

ϕ = αr S1 ( t1 - t2 )

Donde

αr = se denomina coeficiente de radiación. Calor transmitido por radiación y convección Si un cuerpo de temperatura t1 y superficie S1se encuentra dentro de un fluido a temperatura t2, siendo t1 > t2, transmite calor por convección y radiación. (Si estuviera apoyado, también transmitiría calor por conducción a través de los apoyos). El flujo total, transmitido por convección y radiación, según las ecuaciones ya vistas será: ϕ = α . S1 (t1 - t2) + αr . S (t1 - t2) donde α = coeficiente de convección ϕ = (α + αr ) . S1 (t1 - t2)

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 6

TRANSFERENCIA DE CALOR  PROBLEMAS – DESARROLLO:

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 7

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 8

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 9

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 10

TRANSFERENCIA DE CALOR

 CONCLUSIONES 



En el problema #1, a partir de la gráfica superior izquierda, εe disminuye con el aumento de la apertura, Da, como se esperaba. En el límite, cuando Da → 0, ε3 → 1 desde la cavidad se convierte en un recinto completo. En el gráfico superior derecho, el THC, la temperatura de la superficie semiesférica de re-radiante disminuye a medida que aumenta Da. En el límite, cuando Da → 0, T2 se acercará a la temperatura de la placa, Tp = 593 K. El efecto de la disminución de la placa de emisividad es para disminuir εe y disminuir T2. Las hipótesis que asumimos al comienzo, esto es cuerpo gris y superficies difusas así como despreciar la convección a través del medio, sirvió para la determinación de los resultados obtenidos

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 11

TRANSFERENCIA DE CALOR  REFERENCIAS BIBLIOGRAFICAS    

http://www.uco.es/~p32sebej/Energ%C3%A9tica/Radiaci%C3%B3n.pdf http://www.monografias.com/trabajos25/radiacion-termica/radiaciontermica.shtml http://es.scribd.com/doc/27408780/1/CONTENID-OO FUNDAMENTOS DE TRANSFERENCIA DE CALOR”/ Franck IncroperaDavid Dewitt/4°edicion/cap. 12-13

 ANEXOS Programación en matlab. PROBLEMA 1 %DATOS DEL PROBLEMA; e1=0.9;tp=593;ro=0.1;aro=ro/2;sigma=5.67*10^-8;e3=1;talr=300;pi=3.1416; %RELACIONES DE RECIPROCIDAD; syms F12 syms F32 F13=0;F11=0;F33=0;F31=0; F12=solve(F11+F12+F13==1) F32=solve(F31+F32+F33==1) a1=pi*(ro^2-(ro^2/4)); a3=pi*ro^2; q1=sigma*(tp^4-talr^4)/((1-e1)/(e1*a1)+(a1*F13+(1/(a1*F12)+1/(a3*F32))^1)^-1+(1-e3)/(e3*a3))%a) Eef=q1/(a3*sigma*tp^4)%b) Eb1=sigma*tp^4 j3=sigma*talr^4 syms j2 syms j1 syms t2 [j2,j1,t2]=solve(j2==sigma*t2^4,q1==(Eb1-j1)/((1-e1)/(e1*a1)),(j2j1)/(1/(a1*F12))+(j2-j3)/(1/(a3*F32))==0) t2=abs(t2)%c)

%d) %DATOS DEL PROBLEMA; tp=593;ro=0.1;sigma=5.67*10^-8;e3=1;talr=300;pi=3.1416; %RELACIONES DE RECIPROCIDAD; syms F12 syms F32 F13=0;F11=0;F33=0;F31=0; F12=solve(F11+F12+F13==1) F32=solve(F31+F32+F33==1) for e1=[0.5 0.6 0.7]; for aro=ro/8:ro/16:ro/2; a1=pi*(ro^2-(aro./2).^2); a3=pi*ro^2;

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 12

TRANSFERENCIA DE CALOR q1=sigma*(tp^4-talr^4)./((1e1)./(e1.*a1)+(a1.*F13+(1./(a1.*F12)+1/(a3*F32)).^-1).^-1)+(1e3)/(e3*a3); Eef=q1./(a3*sigma*tp^4) Eb1=sigma*tp^4; j3=sigma*talr^4; syms j2 syms j1 syms t2 j1=solve(q1==(Eb1-j1)./((1-e1)./(e1.*a1))); j2=solve((j2-j1)./(1./(a1.*F12))+(j2-j3)/(1/(a3*F32))==0); t2=solve(j2==sigma*t2^4); t2=abs(t2); plot(aro,t2,'*') hold on end end

PROBLEMA2 e1=0.8;t1=40+273;e2=0.9;t2=50+273;e3=0.7;t3=15+173;SIGMA=5.67*10^-8; %para 1 y 2 ; l=3;x=8;y=5; %luego buscamos los valores con los que tendremos que entrar en las tablas X1=x/l; Y1=y/l; F12=2*(log((1+X1^2)*(1+Y1^2)/(1+X1^2+Y1^2))^0.5+X1*(1+Y1^2)^0.5*atan(X1/( 1+Y1^2)^0.5)+Y1*(1+X1^2)^0.5*atan(Y1/(1+X1^2)^0.5)-X1*atan(X1)Y1*atan(Y1))/(3.1416*X1*Y1); F21=F12; %para 4 y 3 ; l=5;x=8;y=3; %luego buscamos los valores con los que tendremos que entrar en las tablas abscisa=x/l; ordenadas=y/l; X1=x/l; Y1=y/l; F34=2*(log((1+X1^2)*(1+Y1^2)/(1+X1^2+Y1^2))^0.5+X1*(1+Y1^2)^0.5*atan(X1/( 1+Y1^2)^0.5)+Y1*(1+X1^2)^0.5*atan(Y1/(1+X1^2)^0.5)-X1*atan(X1)Y1*atan(Y1))/(3.1416*X1*Y1); F43=F34; %para 1 y 3,4 ; z=5;x=8;y=3; %luego buscamos los valores con los que tendremos que entrar en las tablas abscisa=z/x; ordenadas=y/x; H=z/x;W=y/x; F13=(W*atan(1/W)+H*atan(1/H)(H^2+W^2)^0.5*atan(1/(H^2+W^2)^0.5)+0.25*log((W^2*(1*W^2+H^2)/((1+W^2)*(W ^2+H^2)))^(W^2)*(H^2*(1+H^2+W^2)/((1+H^2)*(H^2+W^2)))^(H^2)*(1+W^2)*(1+H^ 2)/(1+W^2+H^2)))/(3.1416*W);

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 13

TRANSFERENCIA DE CALOR F14=F13; %para 2 y 3,4 ; z=5;x=8;y=3; %luego buscamos los valores con los que tendremos que entrar en las tablas abscisa=z/x; ordenadas=y/x; H=z/x;W=y/x; F24=(W*atan(1/W)+H*atan(1/H)(H^2+W^2)^0.5*atan(1/(H^2+W^2)^0.5)+0.25*log((W^2*(1*W^2+H^2)/((1+W^2)*(W ^2+H^2)))^(W^2)*(H^2*(1+H^2+W^2)/((1+H^2)*(H^2+W^2)))^(H^2)*(1+W^2)*(1+H^ 2)/(1+W^2+H^2)))/(3.1416*W); F23=F24; %AREAS DEFINIMOS; X=5;Y=8;Z=3; A1=X*Y; A2=A1; A3=Y*Z; A4=A3; A5A=X*Z; A5B=A5A; %CONSIDERANDO AREAS 5A Y 5B COMO UN SOLO AREA A5 A5=A5A+A5B; %LUEGO TENEMOS POR SIMETRIA F32=A2*F23/A3; F31=A1*F13/A3; F41=F31; F42=F32; %ADEMAS TENEMOS QUE F15=1-F12-F13-F14; F25=1-F21-F23-F24; F35=1-F31-F32-F34; %LAS RELACIONES DE SIMETRIA F51=A1*F15/A5; F53=A3*F35/A5; F45=F35; F52=F51; F54=F53; %LAS POTENCIAS EMISIVAS EB1=SIGMA*t1^4; EB2=SIGMA*t2^4; EB3=SIGMA*t3^4; %determinamos los valores de J syms J1 syms J2 syms J3 syms J4 syms J5 [J1,J2,J3,J4,J5]=solve((EB1-J1)/((1-e1)/e1*A1)==(J1-J2)/(A1/F12)^-1+(J1J3)/(A1/F13)^-1+(J1-J4)/(A1/F14)^-1+(J1-J5)/(A1/F15)^-1,(EB2-J2)/((1e2)/e2*A2)==(J2-J1)/(A2/F21)^-1+(J2-J3)/(A2/F23)^-1+(J2-J4)/(A2/F24)^1+(J2-J5)/(A2/F25)^-1,(EB3-J3)/((1-e3)/e3*A3)==(J3-J1)/(A3/F31)^-1+(J3J2)/(A3/F32)^-1+(J3-J4)/(A3/F34)^-1+(J3-J5)/(A3/F35)^-1,0==(J4J1)/(A4/F41)^-1+(J4-J2)/(A4/F42)^-1+(J4-J3)/(A4/F43)^-1+(J4J5)/(A4/F45)^-1,0==(J5-J1)/(A5/F51)^-1+(J5-J2)/(A5/F52)^-1+(J5J3)/(A5/F53)^-1+(J5-J4)/(A5/F54)^-1)

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 14

TRANSFERENCIA DE CALOR

%hallamos el flujo de radiacion; Q1=A1*F12*(J1-J2)+A1*F13*(J1-J3)+A1*F14*(J1-J4)+A1*F15*(J1-J5) Q2=A2*F21*(J2-J1)+A2*F23*(J2-J3)+A2*F24*(J2-J4)+A2*F25*(J2-J5) Q3=A3*F31*(J3-J1)+A3*F32*(J3-J2)+A3*F34*(J3-J4)+A3*F35*(J3-J5)

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 15

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 16

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 17

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 18

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 19

TRANSFERENCIA DE CALOR

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 20

TRANSFERENCIA DE CALOR

 PROGRAMACION: Problema 1 %DATOS DEL PROBLEMA; e1=0.9;tp=593;ro=0.1;aro=ro/2;sigma=5.67*10^-8;e3=1;talr=300;pi=3.1416; %RELACIONES DE RECIPROCIDAD; syms F12 syms F32 F13=0;F11=0;F33=0;F31=0; F12=solve(F11+F12+F13==1) F32=solve(F31+F32+F33==1) a1=pi*(ro^2-(ro^2/4)); a3=pi*ro^2; q1=sigma*(tp^4-talr^4)/((1-e1)/(e1*a1)+(a1*F13+(1/(a1*F12)+1/(a3*F32))^-1)^-1+(1e3)/(e3*a3))%a) Eef=q1/(a3*sigma*tp^4)%b) Eb1=sigma*tp^4 j3=sigma*talr^4 syms j2 syms j1 syms t2 [j2,j1,t2]=solve(j2==sigma*t2^4,q1==(Eb1-j1)/((1-e1)/(e1*a1)),(j2-j1)/(1/(a1*F12))+(j2j3)/(1/(a3*F32))==0) t2=abs(t2)%c)

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 21

TRANSFERENCIA DE CALOR %d) %DATOS DEL PROBLEMA; tp=593;ro=0.1;sigma=5.67*10^-8;e3=1;talr=300;pi=3.1416; %RELACIONES DE RECIPROCIDAD; syms F12 syms F32 F13=0;F11=0;F33=0;F31=0; F12=solve(F11+F12+F13==1) F32=solve(F31+F32+F33==1) for e1=[0.5 0.6 0.7]; for aro=ro/8:ro/16:ro/2; a1=pi*(ro^2-(aro./2).^2); a3=pi*ro^2; q1=sigma*(tp^4-talr^4)./((1-e1)./(e1.*a1)+(a1.*F13+(1./(a1.*F12)+1/(a3*F32)).^-1).^-1)+(1e3)/(e3*a3); Eef=q1./(a3*sigma*tp^4) Eb1=sigma*tp^4; j3=sigma*talr^4; syms j2 syms j1 syms t2 j1=solve(q1==(Eb1-j1)./((1-e1)./(e1.*a1))); j2=solve((j2-j1)./(1./(a1.*F12))+(j2-j3)/(1/(a3*F32))==0); t2=solve(j2==sigma*t2^4); t2=abs(t2); plot(aro,t2,'*') hold on end end GRAFICAS -RADIO DE ABERTURA VS T2

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 22

TRANSFERENCIA DE CALOR

-RADIO DE ABERTURA VS Ee(emisividad efectiva)

Problema 2 e1=0.8;t1=40+273;e2=0.9;t2=50+273;e3=0.7;t3=15+173;SIGMA=5.67*10^-8; %para 1 y 2 ; l=3;x=8;y=5; %luego buscamos los valores con los que tendremos que entrar en las tablas X1=x/l; INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 23

TRANSFERENCIA DE CALOR Y1=y/l; F12=2*(log((1+X1^2)*(1+Y1^2)/(1+X1^2+Y1^2))^0.5+X1*(1+Y1^2)^0.5*atan(X1/(1+Y1^2)^ 0.5)+Y1*(1+X1^2)^0.5*atan(Y1/(1+X1^2)^0.5)-X1*atan(X1)-Y1*atan(Y1))/(3.1416*X1*Y1); F21=F12; %para 4 y 3 ; l=5;x=8;y=3; %luego buscamos los valores con los que tendremos que entrar en las tablas abscisa=x/l; ordenadas=y/l; X1=x/l; Y1=y/l; F34=2*(log((1+X1^2)*(1+Y1^2)/(1+X1^2+Y1^2))^0.5+X1*(1+Y1^2)^0.5*atan(X1/(1+Y1^2)^ 0.5)+Y1*(1+X1^2)^0.5*atan(Y1/(1+X1^2)^0.5)-X1*atan(X1)-Y1*atan(Y1))/(3.1416*X1*Y1); F43=F34; %para 1 y 3,4 ; z=5;x=8;y=3; %luego buscamos los valores con los que tendremos que entrar en las tablas abscisa=z/x; ordenadas=y/x; H=z/x;W=y/x; F13=(W*atan(1/W)+H*atan(1/H)(H^2+W^2)^0.5*atan(1/(H^2+W^2)^0.5)+0.25*log((W^2*(1*W^2+H^2)/((1+W^2)*(W^2+H^2 )))^(W^2)*(H^2*(1+H^2+W^2)/((1+H^2)*(H^2+W^2)))^(H^2)*(1+W^2)*(1+H^2)/(1+W^2+H^ 2)))/(3.1416*W); F14=F13; %para 2 y 3,4 ; z=5;x=8;y=3; %luego buscamos los valores con los que tendremos que entrar en las tablas abscisa=z/x; ordenadas=y/x; H=z/x;W=y/x; F24=(W*atan(1/W)+H*atan(1/H)(H^2+W^2)^0.5*atan(1/(H^2+W^2)^0.5)+0.25*log((W^2*(1*W^2+H^2)/((1+W^2)*(W^2+H^2 )))^(W^2)*(H^2*(1+H^2+W^2)/((1+H^2)*(H^2+W^2)))^(H^2)*(1+W^2)*(1+H^2)/(1+W^2+H^ 2)))/(3.1416*W); F23=F24; %AREAS DEFINIMOS; X=5;Y=8;Z=3; A1=X*Y; A2=A1; A3=Y*Z; A4=A3; A5A=X*Z; A5B=A5A; %CONSIDERANDO AREAS 5A Y 5B COMO UN SOLO AREA A5 A5=A5A+A5B; %LUEGO TENEMOS POR SIMETRIA F32=A2*F23/A3; F31=A1*F13/A3; F41=F31; F42=F32; INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 24

TRANSFERENCIA DE CALOR %ADEMAS TENEMOS QUE F15=1-F12-F13-F14; F25=1-F21-F23-F24; F35=1-F31-F32-F34; %LAS RELACIONES DE SIMETRIA F51=A1*F15/A5; F53=A3*F35/A5; F45=F35; F52=F51; F54=F53; %LAS POTENCIAS EMISIVAS EB1=SIGMA*t1^4; EB2=SIGMA*t2^4; EB3=SIGMA*t3^4; %determinamos los valores de J syms J1 syms J2 syms J3 syms J4 syms J5 [J1,J2,J3,J4,J5]=solve((EB1-J1)/((1-e1)/e1*A1)==(J1-J2)/(A1/F12)^-1+(J1-J3)/(A1/F13)^1+(J1-J4)/(A1/F14)^-1+(J1-J5)/(A1/F15)^-1,(EB2-J2)/((1-e2)/e2*A2)==(J2-J1)/(A2/F21)^1+(J2-J3)/(A2/F23)^-1+(J2-J4)/(A2/F24)^-1+(J2-J5)/(A2/F25)^-1,(EB3-J3)/((1e3)/e3*A3)==(J3-J1)/(A3/F31)^-1+(J3-J2)/(A3/F32)^-1+(J3-J4)/(A3/F34)^-1+(J3J5)/(A3/F35)^-1,0==(J4-J1)/(A4/F41)^-1+(J4-J2)/(A4/F42)^-1+(J4-J3)/(A4/F43)^-1+(J4J5)/(A4/F45)^-1,0==(J5-J1)/(A5/F51)^-1+(J5-J2)/(A5/F52)^-1+(J5-J3)/(A5/F53)^-1+(J5J4)/(A5/F54)^-1) %hallamos el flujo de radiacion; Q1=A1*F12*(J1-J2)+A1*F13*(J1-J3)+A1*F14*(J1-J4)+A1*F15*(J1-J5) Q2=A2*F21*(J2-J1)+A2*F23*(J2-J3)+A2*F24*(J2-J4)+A2*F25*(J2-J5) Q3=A3*F31*(J3-J1)+A3*F32*(J3-J2)+A3*F34*(J3-J4)+A3*F35*(J3-J5) Por Matrices: >> A=[-1.25 .0105 .0450 .045 0.55; 0.0467 -1.111 0.02 0.02 0.0245; .1286 .1286 -1.4286 .072 .0943; 0.3 .3 .18 -1 0.22; 0.29 0.29 0.176 0.176 -0.932] A= -1.2500

0.0105

0.0450

0.0450

0.5500

0.0467 -1.1110

0.0200

0.0200

0.0245

0.1286

0.1286 -1.4286

0.0720

0.0943

0.3000

0.3000

0.1800 -1.0000

0.2200

0.2900

0.2900

0.1760

0.1760 -0.9320

>> C=[-544.2; -617.2; -390.1; 0;0] INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 25

TRANSFERENCIA DE CALOR C= -544.2000 -617.2000 -390.1000 0 0 >> B=inv(A) B= -1.0206 -0.2475 -0.1400 -0.1780 -0.6650 -0.0637 -0.9331 -0.0287 -0.0366 -0.0736 -0.1501 -0.1635 -0.7428 -0.0971 -0.1910 -0.4515 -0.4916 -0.2365 -1.1483 -0.5744 -0.4510 -0.4911 -0.2374 -0.3019 -1.4473

>> J=B*C J= 762.7880 621.7881 472.3518 641.4502 641.1547

 DRIAGRAMAS DE FLUJO: Problema 1

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 26

TRANSFERENCIA DE CALOR

Problema 2

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 27

TRANSFERENCIA DE CALOR

 REFERENCIAS BIBLIOGRAFICAS:

 ZEMANSKY, MARK W.- CALOR Y TERMODINÁMICA. EDIT. AGUILAR S.A 1979

 SEARS, FRANCIS W.- TERMODINÁMICA. EDITORIAL REVERTÉ, S.A. 1969

 WILSON, JERRY D.- PHYSICS. EDIT.HEAT. SEGUNDA EDICIÓN, 1983

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 28

TRANSFERENCIA DE CALOR

 RESNICK Y HALLIDAY.- FISICA, EDITORIAL CECSA, PARTE I, 1990

 CEIT, UTN FACULTAD REGIONAL BUENOS AIRES.- APUNTE FISICA II B, CALOR Y TERMODINÁMICA, 1995

 FERNÁNDEZ Y GALLONI.- FISICA ELEMENTAL. EDITORIAL NIGAR. BUENOS AIRES.1980

 http://www.textoscientificos.com/fisica/transmision-calor/radiacion

 http://es.scribd.com/doc/27408780/Transferencia-de-Calor-por-Radiacion

INTERCAMBIO DE RADIACION, SUPERFICIES GRISES

Página 29