0% encontró este documento útil (0 votos)
61 vistas67 páginas

Unidad Iii - TC

El documento es un borrador del curso de Transferencia de Calor de la Universidad Nacional de Trujillo, centrado en la tercera unidad sobre la radiación y la condensación. Se explican conceptos fundamentales de la condensación, incluyendo su definición, tipos, y el modelo de Nusselt, así como el coeficiente local de transferencia de calor y la distribución de velocidad en la capa de condensado. Además, se abordan las ecuaciones relacionadas con el flujo másico de condensado y el calor liberado durante el proceso.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
61 vistas67 páginas

Unidad Iii - TC

El documento es un borrador del curso de Transferencia de Calor de la Universidad Nacional de Trujillo, centrado en la tercera unidad sobre la radiación y la condensación. Se explican conceptos fundamentales de la condensación, incluyendo su definición, tipos, y el modelo de Nusselt, así como el coeficiente local de transferencia de calor y la distribución de velocidad en la capa de condensado. Además, se abordan las ecuaciones relacionadas con el flujo másico de condensado y el calor liberado durante el proceso.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

ING ELI GUAYAN H.

BORRADOR RADIACION

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA

BORRADOR DEL CURSO DE:

TRANSFERENCIA DE CALOR
SECCION A

3RA UNIDAD: LA RADIACION


BORRADOR

Mg. ING. ELÍ GUAYÁN H.

TRUJILLO - PERÚ

2021

ING. ELÍ GUAYÁN HUACCHA Página 1


ING ELI GUAYAN H. BORRADOR RADIACION

LA CONDENSACIÓN (Borrador)

CONCEPTO
Es un fenómeno físico de cambio de estado, del estado gaseoso al estado líquido
La condensación ocurre cuando la temperatura de un vapor se reduce por debajo
de su temperatura de saturación del vapor, en el valor correspondiente de la
presión de trabajo. La energía latente del vapor se libera, y el calor se transfiere
a la superficie, y se forma la condensación.
Esquema:

Plano T-v

Tipos:
De película
Tipo gotas

Condensación de película sobre pared plana vertical: Modelo de Nusselt:


Esquema:

Hipótesis de trabajo:
La pared a temperatura uniforme Tp<Tv (menor que).
El vapor es estacionario y no ejerce ninguna fuerza de arrastre sobre el movimiento
del condensado.
El flujo de condensado es laminar.
No hay aceleración del fluido dentro de la capa de condensado
Propiedades termofísicas del fluido son constantes..
La transferencia de calor a través de la capa de condensado se realiza por conducción
pura.
La distribución de temperatura es lineal.

ING. ELÍ GUAYÁN HUACCHA Página 2


ING ELI GUAYAN H. BORRADOR RADIACION

EL COEFICIENTE LOCAL DE TRANSFERENCIA ℎ(𝑥) :

El coeficiente local de Transferencia de calor ℎ(𝑥) :

𝑑𝑄̇𝑐𝑜𝑛𝑣𝑒𝑐. = 𝑑𝑄̇𝑐𝑜𝑛𝑑𝑢𝑐.
𝑑𝑇 ∆𝑇
𝑑𝑄̇ = ℎ(𝑥) . 𝑑𝐴. ∆𝑇 = −𝑘𝑑𝐴 → 𝐸𝑠 𝑙𝑖𝑛𝑒𝑎𝑙 ∶
𝑑𝑦 ∆𝑦
(𝑇𝑣 − 𝑇𝑝 )
𝑑𝑄̇ = ℎ(𝑥) 𝑑𝐴(𝑇𝑣 − 𝑇𝑝 ) = 𝐾. 𝑑𝐴
𝛿𝑥
𝑘
ℎ(𝑥) = → (1)
𝛿𝑥
DISTRIBUCIÓN DE VELOCIDAD 𝑢(𝑦) : a través de la capa de condensado.

Considerando un elemento diferencial de volumen:

𝐹
𝜏= → 𝐹 = 𝜏. 𝐴
𝐴
En forma diferencial:

𝑑𝐹 = 𝜏. 𝑑𝐴
𝑝𝑒𝑠𝑜
𝛾= ⟶ 𝑝 = 𝛾. 𝑉𝑜𝑙
𝑣𝑜𝑙.
𝑝 = 𝜌. 𝑔. 𝑉𝑜𝑙
Diferencial de peso:

ING. ELÍ GUAYÁN HUACCHA Página 3


ING ELI GUAYAN H. BORRADOR RADIACION

𝑑𝑝 = 𝜌𝑙 . 𝑔. 𝑑𝑉 = 𝜌𝑙 . 𝑔 𝑑𝑥. 𝑑𝑦. 𝑑𝑧

1

EMPUJE (𝐸 ):

𝐸 = 𝛾𝑉 . 𝑉𝑜𝑙. 𝑑𝑒𝑠𝑝𝑙𝑎𝑧𝑎𝑑𝑜
𝐸 = 𝜌𝑉 . 𝑔. 𝑉𝑜𝑙.
En forma diferencial:

𝑑𝐸 = 𝜌𝑣 . 𝑔. 𝑑𝑥. 𝑑𝑦(1)
Sumando las fuerzas verticales: ∑ 𝐹𝑋 = 0 ↓ 𝑥(+)
𝜕𝜏
𝜌𝑙 . 𝑔. 𝑑𝑥. 𝑑𝑦 − 𝜏. 𝑑𝑥 + 𝜏. 𝑑𝑥 + 𝑑𝑥. 𝑑𝑦 − 𝜌𝑣 . 𝑔. 𝑑𝑥. 𝑑𝑦 = 0
𝜕𝑦
𝜕𝜏
Simplificando: 𝜌𝑙 𝑔 + − 𝜌𝑣 𝑔 = 0
𝜕𝑦

𝜕𝜏
= −(𝜌𝐿 − 𝜌𝑉 )𝑔
𝜕𝑦
𝜕𝑢
Además: 𝜏 = 𝜇𝑙 (Ley de la viscosidad de Newton)
𝜕𝑦

Reemplazando:
𝜕 𝜕𝑢
(𝜇𝑙 ) = −(𝜌𝑙 − 𝜌𝑣 )𝑔
𝜕𝑦 𝜕𝑦

𝜕 2𝑢 𝜕 2𝑢 𝜌𝑙 − 𝜌𝑣
𝜇𝑙 = −(𝜌𝑙 − 𝜌𝑣 )𝑔 → = −( )𝑔
𝜕𝑦 2 𝜕𝑦 2 𝜇𝑙
𝑑2 𝑢 𝜌𝑙 −𝜌𝑣
Ya es derivada total: = −( )𝑔
𝑑𝑦 2 𝜇𝑙

𝑑𝑢 𝜌𝑙 −𝜌𝑣
Integrando: = −( ) . 𝑔. 𝑦 + 𝐶1
𝑑𝑦 𝜇𝑙

(𝜌𝑙 −𝜌𝑣 )𝑔.𝑦 2


Integrando: 𝑢 = − + 𝐶1 𝑦 + 𝐶2
2𝜇𝑙

Condiciones de frontera:

Para 𝑦 = 0, 𝑢 = 0 (Sobre la superficie)


𝑑𝑢 𝜕𝑢
Para 𝑦 = 𝛿 → =0 (𝜏 = 𝜇𝑙 𝜕𝑦 = 0) en la interface
𝑑𝑦

Reemplazando: para 𝑦 = 0 → 𝑢 = 0 → 𝐶2 = 0
𝑑𝑢 (𝜌𝑙 −𝜌𝑣 )
Para 𝑦 = 𝛿 → = 0 → 𝐶1 = 𝑔𝛿
𝑑𝑦 𝜇𝑙

(𝜌𝑙 − 𝜌𝑣 )𝑔. 𝑦 2 (𝜌𝑙 − 𝜌𝑣 )


→ 𝑢=− + 𝑔. 𝛿. 𝑦
2𝜇𝑙 𝜇𝑙

ING. ELÍ GUAYÁN HUACCHA Página 4


ING ELI GUAYAN H. BORRADOR RADIACION

Factorizando:
𝜌𝑙 −𝜌𝑣 𝑦2
𝑢(𝑦) = ( ) 𝑔 (𝑦. 𝛿 − ) distribución de velocidad en la capa de condensado
𝜇𝑙 2

EL FLUJO MÁSICO DE CONDENSADO (𝑚̇)


Que pasa a través de la sección y que tiene profundidad unitaria
𝑚̇ = 𝜌𝐿 ∗ 𝑉𝑒𝑙𝑜𝑐 ∗ Á𝑟𝑒𝑎

ℎ𝑓𝑔 =calor latente de condensación

El diferencial de masa: 𝑑𝑚̇ = 𝜌𝐿 . 𝑢. 𝑑𝐴


𝑚̇ 𝛿
Integrando: ∫0 𝑑𝑚̇ = ∫0 𝜌𝐿 𝑢. 𝑑𝐴

𝛿 𝑦2
(𝜌𝑙 − 𝜌𝑣 )𝑔 (𝑦𝛿 − )
2
𝑚̇ = ∫ 𝜌𝐿 [ ] 𝑑𝑦 ∗ (1)
𝜇𝑙
0

𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑔 𝛿 3
→ 𝑚̇ = .
𝜇𝑙 3

Diferenciando con respecto a 𝛿 :

(𝜌𝐿 − 𝜌𝑣 )
→ 𝑑𝑚̇ = 𝜌𝑙 𝑔. 𝛿 2 . 𝑑𝛿
𝜇𝐿
Es la cantidad en que aumenta el condensado en la longitud 𝑋, cuando pasa de una
posición 𝑥 a (𝑥 + 𝑑𝑥), y el espesor de la capa de condensado aumenta de 𝛿 a
(𝛿 + 𝑑𝛿).

EL CALOR LIBERADO (𝒅𝑸̇):

ING. ELÍ GUAYÁN HUACCHA Página 5


ING ELI GUAYAN H. BORRADOR RADIACION

∆𝑇 = +(𝑇𝑣 − 𝑇𝑝 )

En 𝑥 = 0 → 𝛿 = 0

Durante la condensación:
∆𝑇 𝜕𝑇 ∆𝑇
𝑑𝑄̇ = ℎ𝑓𝑔 . 𝑑𝑚̇ = 𝑘𝐿 . 𝑑𝐴 ; ( ≅ )
𝛿 𝜕𝑌 𝛿
(𝜌𝑙 − 𝜌𝑣 ) ∆𝑇
𝑑𝑄̇ = ℎ𝑓𝑔 . 𝜌𝑙 . 𝑔. 𝛿 2 . 𝑑𝛿 = 𝑘𝑙 . (1)𝑑𝑥

𝜇𝐿 𝛿 𝑑𝐴
𝛿𝑥 𝑥
𝑘𝑙 . ∆𝑇. 𝜇𝑙
∴ ∫ 𝛿 3 . 𝑑𝛿 = ∫ 𝑑𝑥
0 0 ℎ𝑓𝑔 . 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑔
1⁄
4𝑘𝑙 . ∆𝑇. 𝜇𝑙 . 𝑥 4 𝐸𝑠𝑝𝑒𝑠𝑜𝑟 𝑑𝑒 𝑝𝑒𝑙í𝑐𝑢𝑙𝑎 𝑡𝑜𝑡𝑎𝑙
𝛿(𝑥) =[ ]
ℎ𝑓𝑔 . 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑔 𝑐𝑜𝑛 ∆𝑇 = (𝑇𝑣 − 𝑇𝑝 )
𝑘𝐿
Entonces de: ℎ(𝑥) =
𝛿(𝑥)

1⁄
ℎ𝑓𝑔 .𝜌𝑙 (𝜌𝑙 −𝜌𝑣 )𝑔.𝑘𝑙3 4
ℎ(𝑥) = [ ] coeficiente de película local
4.∆𝑇.𝜇𝑙 .𝑥

ℎ(𝑥) 𝑋
También: 𝑁𝑢(𝑥) =
𝑘𝐿

El coeficiente ℎ̅ = es una longitud 0 ≤ 𝑥 ≤ 𝐿


𝐿 1⁄
1 1 ℎ𝑓𝑔 . 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑔𝑘𝑙3 4 𝐿
ℎ̅ = ∫ ℎ(𝑥) 𝑑𝑥 = [ ] ∫ 𝑋 −1⁄4 . 𝑑𝑥
𝐿 𝐿 4. ∆𝑇. 𝜇𝑙 . 𝑥 0
0

1⁄
ℎ𝑓𝑔 . 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑔. 𝑘𝑙3 4
̅
→ ℎ = 0,943 [ ] 𝑇𝑒𝑜𝑟𝑖𝑐𝑜
𝜇𝑙 . ∆𝑇. 𝐿
4
̅
Óℎ = ℎ(𝑥) |𝑋=𝐿
3
𝐿 = Longitud de placa vertical

ING. ELÍ GUAYÁN HUACCHA Página 6


ING ELI GUAYAN H. BORRADOR RADIACION

NÚMERO DE JAKOB (𝐽𝑎):


Razón de energía sensible a latente absorbida durante el cambio de fase liquido-vapor

𝐶𝑝𝑙 (𝑇𝑠𝑎𝑡 − 𝑇𝑝 )
𝐽𝑎 =
ℎ𝑓𝑔

CALOR LATENTE MODIFICADO DE VAPORIZACIÓN:


Según Rohsenow:

ℎ𝑓𝑔 = ℎ𝑓𝑔 + 0,68. 𝐶𝑝𝐿 . (𝑇𝑠𝑎𝑡 − 𝑇𝑝 )

Donde: 𝐶𝑝𝐿 =Calor específico del líquido a la temperatura promedio de película


𝑇𝑝 =Temperatura de pared

O en función del número de Jakob: 𝐽𝑎



ℎ𝑓𝑔 = ℎ𝑓𝑔 + (1 + 0,68 𝐽𝑎)
CORRECCION Para placa Plana Vertical:
1⁄

ℎ𝑓𝑔 . 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )𝑔. 𝑘𝑙3 4
̅
ℎ = 0,943 [ ]
𝜇𝑙 . ∆𝑇. 𝐿

∆𝑇 = 𝑇𝑠𝑎𝑡 − 𝑇𝑝

CONDENSACIÓN DE PELÍCULA SOBRE UNA PLACA INCLINADA:

Con: 𝜃 ≤ 60°
1
ℎ̅ = ℎ̅𝑣𝑒𝑟𝑡𝑖𝑐𝑎𝑙 . (cos 𝜃) ⁄4
1⁄
𝑔.𝜌𝑙 (𝜌𝑙 −𝜌𝑣 )ℎ𝑓𝑔 .𝑘𝑙3 .cos 𝜑 4
Coeficiente local: ℎ(𝑥) = 1.2 [ ]
4𝜇𝑙 (𝑇𝑣 −𝑇𝑝 )𝑥

1⁄

𝑔.𝜌𝑙 (𝜌𝑙 −𝜌𝑣 )ℎ𝑓𝑔 .𝑘𝑙3 .cos 𝜑 4
Coeficiente promedio: ℎ(𝑥) = 0,943 [ ]
𝜇𝑙 (𝑇𝑣 −𝑇𝑝 )𝐿

Donde: 𝜌𝑙 = Densidad del líquido (𝐾𝑔⁄𝑚3 )

ING. ELÍ GUAYÁN HUACCHA Página 7


ING ELI GUAYAN H. BORRADOR RADIACION

𝜌𝑣 = Densidad del vapor (𝐾𝑔⁄𝑚3 )


ℎ𝑓𝑔 = Calor estante latente de condensación (𝐽⁄𝐾𝑔 ó 𝑊. 𝑠⁄𝐾𝑔)

𝑘𝑙 = Conductividad térmica del líquido (𝑊 ⁄𝑚. °𝐶 )


𝐿 = Longitud de la superficie (𝑚)
𝑇𝑣 = Temperatura del vapor (°𝐶)
𝑇𝑝 = Temperatura de la superficie de la pared (°𝐶)

𝜇𝑙 = Viscosidad del líquido (𝐾𝑔⁄𝑚. 𝑠)


CONDENSACIÓN SOBRE TUBOS HORIZONTALES (EXTERNO):
1⁄

𝑔. 𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )ℎ𝑓𝑔 . 𝑘𝑙3 4
ℎ̅ = 0,729 [ ]
𝜇𝑙 (𝑇𝑣 − 𝑇𝑝 )𝐷𝑒

Donde: 𝐷𝑒 = diámetro exterior de la tubería

CONDENSACIÓN SOBRE UN HAZ HORIZONTAL DE TUBOS:

1⁄
𝑔𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 ) ℎ∗𝑓𝑔 . 𝑘3𝑙 4
̅
ℎ = 0,729 [ ]
𝜇𝑙 (𝑇𝑣 − 𝑇𝑝 )𝑁. 𝐷𝑒

Donde: 𝑁 = N° de tubos horizontales dispuestos en


hilera vertical.

• Si tenemos: Tubería vertical de longitud L ℎ̅𝑡𝑣 = ℎ̅𝑡ℎ


Tubería longitud de diámetro ext. D.

• Si se tiene una tubería 𝐿 = 100𝐷 → ℎ̅𝑡ℎ = 2.44ℎ̅𝑡𝑣


PROPIEDADES FISICAS DEL CONDENSADO:
Las propiedades físicas del condensado se evalúan a la temperatura media
𝑇𝑚 (𝑘𝑙 , 𝜇𝑙 , 𝜌𝑙 )
𝑇𝑣 + 𝑇𝑝
𝑇𝑚 =
2

ING. ELÍ GUAYÁN HUACCHA Página 8


ING ELI GUAYAN H. BORRADOR RADIACION

Las propiedades del vapor se evalúan a 𝑇𝑠𝑎𝑡


Las propiedades del ℎ𝑓𝑔 se evalúan a 𝑇𝑠𝑎𝑡

REGIMENES DE FLUJO

Flujo laminar ondulado:


Kutateladze:
𝑅𝑒 𝑘𝐿 𝑔 1/3
ℎ= ( ) 30 < 𝑅𝑒 < 1800
1,08𝑅𝑒1,22 − 5,2 𝜈

Flujo turbulento sobre placas verticales:


Labuntsov:
𝑅𝑒 𝑘𝐿 𝑔 1/3
ℎ= ( ) 𝑅𝑒 > 1800
8750 + 58𝑃𝑟 −0,5 ( 𝑅𝑒 0,75 − 253) 𝜈

NÚMERO DE REYNOLDS DEL FLUJO DE CONDENSADO:


𝜌𝑙 . 𝑢𝑚 . 𝐷𝑒
𝑅𝑒 =
𝜇𝑙
Donde: 𝑢𝑚 = Velocidad media de la película de condensado
𝐷𝑒 = Diámetro hidráulico el flujo condensando
4𝐴 4. (á𝑟𝑒𝑎 𝑑𝑒 𝑙𝑎 𝑠𝑒𝑐𝑐𝑖ó𝑛 𝑡𝑟𝑎𝑛𝑠𝑣𝑒𝑟𝑠𝑎𝑙 𝑑𝑒𝑙 𝑓𝑙𝑢𝑗𝑜 𝑐𝑜𝑛𝑑𝑒𝑛𝑠𝑎𝑑𝑜)
𝐷𝑒 ≡ =
𝑝 𝑝𝑒𝑟í𝑚𝑒𝑡𝑟𝑜 ℎú𝑚𝑒𝑑𝑜
Para Placa Vertical:

𝑝=𝑏
𝐴 = 𝑏. 𝛿
𝐷𝑒 = 4𝛿
4. 𝛿. 𝑢𝑚 . 𝜌𝑙 4𝛤
∴ 𝑅𝑒 = ( )=
𝜇𝑙 𝜇𝑙
Donde: 𝛤 = 𝜌𝑙 . 𝑢𝑚 . 𝛿 = Tasa de flujo másico por unidad de anchura (𝐾𝑔⁄𝑚. 𝑠)
Para Tubo Vertical: (𝑅𝑒 en el extremo inferior)

ING. ELÍ GUAYÁN HUACCHA Página 9


ING ELI GUAYAN H. BORRADOR RADIACION

4 𝑚̇
𝑅𝑒 = ( )
𝜇𝑙 𝜋𝐷
𝑚̇: tasa de flujo másico de condensado en el extremo
inferior del flujo vertical (𝐾𝑔⁄𝑠)

Balance de energia:

𝑄̇ 𝑙𝑎𝑡𝑒𝑛𝑡𝑒 = 𝑄̇𝑐𝑜𝑛𝑣𝑒𝑐𝑐𝑖ó𝑛
𝑙𝑖𝑏𝑒𝑟𝑎𝑑𝑜

𝑚̇. ℎ𝑓𝑔 = ℎ̅. 𝐴𝑐𝑜𝑛𝑑𝑒𝑛𝑠. (𝑇𝑣 − 𝑇𝑝 )

𝑚̇. ℎ𝑓𝑔 = ℎ̅. 𝜋. 𝐷𝑒 . 𝐿(𝑇𝑣 − 𝑇𝑝 )

La tasa de flujo másico condensado (𝒎̇):

𝜋𝐷𝑒 . 𝐿. ℎ̅(𝑇𝑣 − 𝑇𝑝 )
𝑚̇ = ∗
ℎ𝑓𝑔

EJERCICIO DE CONDENSACION
Un tubo horizontal de 1m de longitud con una temperatura superficial de 70°C. se usa para
condensar vapor saturado a 1 atm.
Que diámetro se requiere para alcanzar un flujo de condensación de 125 Kg/h?
SOLUCION:
Según datos se tiene que el Tubo horizontal de 1 m de largo con temperatura superficial de
70 ° C que se utiliza para condensar vapor a 1 atm.
Se requiere hallar el diámetro para alcanzar un flujo de condensación de 125 Kg/h.
Esquema:

𝑇𝑠 = 70°𝐶 = 343°𝐾

𝑘𝑔
𝑚̇ = 125

D

Suponer:
• Condensación de película laminar.
• No condensables insignificantes en vapor.

ING. ELÍ GUAYÁN HUACCHA Página 10


ING ELI GUAYAN H. BORRADOR RADIACION

Según Tabla: vapor de agua (1atm):


𝑇𝑠𝑎𝑡 = 100°𝐶, 𝜌𝑣 = 0,5978 𝑘𝑔/𝑚3 y ℎ𝑓𝑔 = 2257 𝑘𝐽/𝑘𝑔
𝑇𝑠 +𝑇𝑠𝑎𝑡 70+100
Según Tabla: agua líquida ( 𝑇𝑓 = ( )= = 85°𝐶
2 2

𝜌𝑙 = 968,1 𝑘𝑔/𝑚3 , 𝑐𝑙,𝑝 = 4201 𝐽/𝑘𝑔. °𝐾, 𝑢𝑙 = 333𝑥10−6 𝑁. 𝑠/𝑚2 y


W
𝑘𝑙 = 0,673 m . °K

Calor latente modificado:



ℎ𝑓𝑔 = ℎ𝑓𝑔 + 0,68. 𝐶𝑝𝐿 . (𝑇𝑠𝑎𝑡 − 𝑇𝑝 )

ℎ𝑓𝑔 = 2257 + 0,68 ∗ 4,201(100 − 70) = 2342,7004𝐾𝐽/𝐾𝑔
Del balance:

𝑚̇. ℎ𝑓𝑔 = ℎ̅. 𝜋. 𝐷𝑒 . 𝐿(𝑇𝑣 − 𝑇𝑝 )
De la correlación del coeficiente de convección para condensación:
1⁄
𝑔𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )ℎ∗𝑓𝑔 . 𝑘3𝑙 4
̅
ℎ = 0,729 [ ]
𝜇𝑙 (𝑇𝑣 − 𝑇𝑝 ). 𝐷𝑒

Sustituyendo los valores en (I):


125
∗ 2342,7103 = ℎ̅. 𝜋. 𝐷𝑒 ∗ 1 ∗ (100 − 70)
3600
863,08
ℎ̅ =
𝐷𝑒
Luego reemplazando en:
1⁄
𝑔𝜌𝑙 (𝜌𝑙 − 𝜌𝑣 )ℎ∗𝑓𝑔 . 𝑘3𝑙 4
̅
ℎ = 0,729 [ ]
𝜇𝑙 (𝑇𝑣 − 𝑇𝑝 ). 𝐷𝑒
Donde:

𝐷 = 0,144m = 14,4cm

ING. ELÍ GUAYÁN HUACCHA Página 11


ING ELI GUAYAN H. BORRADOR RADIACION

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA

BORRADOR DEL CURSO DE:

TRANSFERENCIA DE CALOR
SECCION A

3RA UNIDAD: LA RADIACION


BORRADOR

Mg. ING. ELÍ GUAYÁN H.

TRUJILLO - PERÚ

2021

ING. ELÍ GUAYÁN HUACCHA Página 1


ING ELI GUAYAN H. BORRADOR RADIACION

LA RADIACIÓN (Borrador)

Concepto.- Es el tercer mecanismo de transferencia de calor y que se debe a la


propagación de ondas electromagnéticas (que viajan a la velocidad de la luz) emitidas por
un cuerpo o sustancia que está a determinada temperatura.

La radiación no necesita de un medio material para propagarse por que se da inclusive en


el vacío, la radiación es emitida por los cuerpos como consecuencia de las variaciones del
estado de cambio de los arreglos electrónicos y moleculares de la sustancia. En
aplicaciones de ingeniería se asume la radiación como un fenómeno superficial.

Radiacio n de
Radiación
los alrededde
ores
los alrededore s q rad , n et

Em isión
E misio de
n de
radiac ión
radiacion
su perficial
superf icial
T a lr
T S
VacíoVacio
Alrede dores
Alrededores

min = 0ºK absoluto


Tmin absoluto

Cresta

Valle

Espectro Electromagnético

Red Termica (0,1 - 100um)

-12 -10 -8 -6 -4 -2 0 2 4
10 10 10 10 10 10 10 10 10
(m)
Rayos X long.
Ondas de radio
onda
Infrarrojo TV, radar
Rayos

Rayos Ultra Luz visible (0,38 - 0,78um)


cosmicos violeta

ING. ELÍ GUAYÁN HUACCHA Página 2


ING ELI GUAYAN H. BORRADOR RADIACION

El espectro electromagnético nos da la característica de los distintos tipos de


radiación y nosotros nos referiremos a la radiación térmica.
1 micra = μm = 10-6 m
COLOR BANDA DE LONGITUDES DE ONDA
Violeta 0,40- 0,44 μm
Azul 0,44- 0,49 μm
Verde 0,49- 0,54 μm
Amarillo 0,54- 0,60 μm
Naranja 0,60- 0,67 μm
Rojo 0,63- 0,76 μm
Radiación de cuerpos negros.- Los cuerpos reales a pesar que pudieran tener la misma
temperatura no emiten la misma cantidad de radiación ya que esta depende de la naturaleza
del cuerpo, de su superficie, de la temperatura, por ello es que se ha creado un modelo de
comparación llamado el cuerpo negro, que es el que emite y absorbe la máxima cantidad de
radiación posible a cualquier temperatura y en cualquier longitud de onda, sin reflejarla.

El modelo de cuerpo negro proporciona una referencia muy útil para compararla con la
emisión radiante de cuerpos reales.

Qrad
Qrad real

Sup. Sup.
Negra Real
(T1) (T1)

LEY DE PLANCK:
Quien mediante los principios de la teoría cuántica elabora una relación matemática para el
poder emisivo o potencia emisiva de un cuerpo negro a una determinada longitud de onda.

Distribución de Planck:
𝐶1 𝑊
𝐸𝑏𝜆 = ; [ ]
𝜆 (𝑒 2⁄𝜆𝑇
5 𝐶 − 1) 𝑚2 . 𝜇𝑚
Donde: 𝐶1 = 1era constante de Planck (3,743𝑥108 𝑊. 𝜇𝑚4 ⁄𝑚2 )
𝐶2 = 2da constante de Planck (1,4387𝑥104 𝜇𝑚. °𝐾)
𝜆 = Longitud de onda (𝜇𝑚)
𝑇 = Temperatura absoluta del cuerpo negro (°𝐾)

ING. ELÍ GUAYÁN HUACCHA Página 3


ING ELI GUAYAN H. BORRADOR RADIACION

E b
2
(W/m um)
T 1

E b max
T 2

POTENCIA
EMISIVA T 3
ESPECTRAL
DEL CUERPO
NEGRO T 4

(u m)
ma x

Poder emisivo total del cuerpo negro (𝑬𝒃 ):


Se obtiene integrando la ecuación de Planck en todo el espectro electromagnético. (En todas
las longitudes de onda).

𝑊
𝐸𝑏 = ∫ 𝐸𝑏𝜆 . 𝑑𝜆 [ ]
0 𝑚2
∞ 𝐶 .𝑑𝜆
𝐸𝑏 = ∫0 𝜆5 (𝑒 𝐶12⁄𝜆𝑇 −1) , 𝐸𝑏 = Área bajo la curva

Ley de desplazamiento de Wein:


Derivando la distribución de Planck:
𝑑𝐸𝑏𝜆
=0
𝑑𝜆
Y con las hipótesis simplificadoras:
𝜆𝑚á𝑥 . 𝑇 ≅ 2897,8 𝜇𝑚. °𝐾
{
𝜆𝑚á𝑥 . 𝑇 ≅ 5215,6 𝜇𝑚. °𝑅

Ley de Steffan – Boltzman:

Establece que la máxima radiación que puede emitir un cuerpo negro (potencia emisiva
total) a una determinada temperatura por unidad de superficie y en todas las longitudes
de onda es:
𝑄̇𝑏 = 𝐴. 𝜎. 𝑇 4

ING. ELÍ GUAYÁN HUACCHA Página 4


ING ELI GUAYAN H. BORRADOR RADIACION

𝑄̇𝑏
= 𝜎. 𝑇 4
𝐴
𝐸𝑏 = 𝜎. 𝑇 4 [𝑊 ⁄𝑚2 ]
Donde: 𝜎 = 5,67𝑥10−8 𝑊 ⁄𝑚2 . °𝐾 4
Eb= potencia emisiva total
Emisión de banda:
Energía radiada desde un cuerpo negro en un intervalo de longitud de onda:

𝐸𝑏𝜆1−𝜆2 =?
𝜆1 𝐶1
Hasta 𝜆1 : 𝐸𝑏0−𝜆 = ∫0 𝑑𝜆
1 𝜆5 (𝑒 𝐶2 ⁄𝜆𝑇 −1)

𝜆2 𝐶1
Hasta 𝜆2 : 𝐸𝑏0−𝜆 = ∫0 𝑑𝜆
2 𝜆 (𝑒 ⁄𝜆𝑇 −1)
5 𝐶 2

En general:
𝜆
𝐶1
𝐸𝑏0−𝜆 = ∫ 𝑑𝜆
0 𝜆5 (𝑒 𝐶2⁄𝜆𝑇 − 1)
Dividiendo entre: 𝐸𝑏0−∞

𝜆 𝐶1
𝐸𝑏0−𝜆 ∫0 5 𝐶2⁄𝜆𝑇 𝑑𝜆
𝜆 (𝑒 − 1)
=
𝐸𝑏0−∞ ∞ 𝐶1
∫0 5 𝐶2⁄𝜆𝑇 𝑑𝜆
𝜆 (𝑒 − 1)
Pero:
𝜆 𝐶1
𝐸𝑏0−𝜆 ∫0 𝑑𝜆
𝜆 (𝑒 2⁄𝜆.𝑇
5 𝐶 − 1)
=
𝐸𝑏 𝜎. 𝑇 4
Acomodando:
𝐸𝑏0−𝜆 𝜆 𝐶1 .𝑇
= ∫0 𝑑𝜆 con: 𝑇. 𝑑𝜆 = 𝑑(𝜆𝑇)
𝐸𝑏 𝜎.𝜆5 .𝑇 4 .𝑇(𝑒 𝐶2 ⁄𝜆.𝑇 −1)

𝐸𝑏0−𝜆 𝜆𝑇
𝐶1 . 𝑇
=∫ 5 𝐶2 ⁄𝜆.𝑇 − 1)
𝑑(𝜆𝑇)
𝐸𝑏 0 𝜎. (𝜆. 𝑇) (𝑒

Poder emisivo para una Franja:

ING. ELÍ GUAYÁN HUACCHA Página 5


ING ELI GUAYAN H. BORRADOR RADIACION

𝐸𝑏0−𝜆 𝐸𝑏0−𝜆
2 1
𝐸𝑏𝜆1−𝜆2 = 𝐸𝑏0−∞ [ 𝐸 −𝐸 ] En tablas.
𝑏0−∞ 𝑏0−∞

Ejercicios de aplicación:
1. Determine el % de radiación emitida por un cuerpo negro a la temperatura de 300ºK que se
encuentra entre las longitudes de onda de 8 – 14 μm.
Solución:

Eb

T=300ºK

8 14

Datos:

T = 300 ºK

𝜆1 = 8 𝜇𝑚 𝜆1 𝑇 = (8𝜇𝑚)(300𝑜 𝐾) = 240 𝜇𝑚𝐾


𝜆2 = 14 𝜇𝑚 𝜆2 . 𝑇 = (14𝜇𝑚)(300𝑜 𝐾) = 4200 𝜇𝑚𝐾
𝐸𝑏𝜆1−𝜆2 =?

De tablas tenemos:

𝐸𝑏0−8 𝐸𝑏0−14
𝐸𝑏
= 0,140290 , 𝐸𝑏
= 0,516044

𝐸𝑏8−14 = 0,516044 − 0,140290


𝐸𝑏8−14 = 0,375754

Rpta: 37,57%

ING. ELÍ GUAYÁN HUACCHA Página 6


ING ELI GUAYAN H. BORRADOR RADIACION

Ejemplo:
Un elemento de tungsteno se calienta a la temperatura de 2000º K. ¿A qué longitud de
onda la cantidad de radiación emitida es máxima?
Solución:

E b

Eb max

T=200ºK

max

Datos:
𝑇 = 2000𝑜 𝐾
𝜆𝑚𝑎𝑥 =?
De la ley de Wein:
𝜆𝑚𝑎𝑥 𝑇 = 2897,8𝜇𝑚𝑜 𝐾
2897,8𝜇𝑚𝑜 𝐾
𝜆𝑚𝑎𝑥 =
2000𝑜 𝐾

𝝀𝒎𝒂𝒙 = 𝟏, 𝟒𝟒𝟗𝝁𝒎

ING. ELÍ GUAYÁN HUACCHA Página 7


ING ELI GUAYAN H. BORRADOR RADIACION

INTENSIDAD DE RADIACIÓN (𝑰): (Borrador)


La radiación que emite una superficie se propaga en todas las direcciones posibles y puede ser
no uniforme. Si se necesita conocer la radiación en una determinada dirección, entonces se
habla de la distribución direccional. La intensidad de radiación se refiere en concreto a la
radiación térmica emitida en una dirección en particular.

𝑑𝐴𝑠 = 𝑑𝐴1 . cos 𝜑


𝑑𝜔 = 𝑑𝑖𝑓. 𝑑𝑒 á𝑛𝑔𝑢𝑙𝑜 𝑠ó𝑙𝑖𝑑𝑜 (𝑠𝑟)
Donde:
𝜑: Ángulo cenital
𝜃: Ángulo Azimutal
Intensidad de Radiación Total: (𝑰𝑻 ):
Se define como la velocidad a la cual la energía radiante se emite en la dirección (φ, θ), por
unidad de área de la superficie emisora normal a dicha dirección y por unidad de ángulo sólido
en torno a esta misma dirección.
Intensidad de Radiación Espectral:
Se define como la razón a la que se emite energía radiante a la longitud de onda λ en la dirección
(φ, θ) por área de la superficie emisora normal a dicha dirección y por unidad de ángulo sólido
alrededor de ésta dirección y por intervalo de longitud de onda unitaria dλ; alrededor de λ.

Ángulo Plano Ángulo Sólido

𝑑𝐴𝑁
𝑑𝑆 = 𝑅. 𝑑𝜃 𝑑𝜔 = 𝑅2

ING. ELÍ GUAYÁN HUACCHA Página 8


ING ELI GUAYAN H. BORRADOR RADIACION

Intensidad espectral:

𝑑𝑄̇ 𝑊
𝐼𝜆, 𝑒 = ; [ 2 ]
𝑑𝐴1 . cos 𝜑 . 𝑑𝜔. 𝑑𝜆 𝑚 . 𝑠𝑟. 𝜇𝑚
Intensidad total:

𝑑𝑄̇ 𝑊
𝐼𝑒 = ; [ 2 ]
𝑑𝐴1 . cos 𝜑 . 𝑑𝜔 𝑚 . 𝑠𝑟

𝐼𝑇 = ∫ 𝐼𝜆 . 𝑑𝜆
0

𝑑𝑄̇ 𝑊
≡ 𝑑𝑄̇𝜆 ; [ ]
𝑑𝜆 𝜇𝑚

𝑑𝑄̇𝜆 = 𝐼𝜆, 𝑒 . 𝑑𝐴1 . cos 𝜑 . 𝑑𝜔 ………………………………………. (𝛼)

Hemisferio:

𝑑𝐴𝑁
𝑑𝜔 = …………….. (*)
𝑟2

⟹De la figura:

𝑑𝐴𝑁 = 𝑟 sen 𝜑𝑑𝜃. 𝑟𝑑𝜑


𝑑𝐴𝑁 = 𝑟 2 sen 𝜑𝑑𝜃. 𝑑𝜑; reemplazando en (*):

𝑟 2 .sen 𝜑.𝑑𝜃𝑑𝜑
𝑑𝜔 = ; 𝑑𝜔 = sen 𝜑𝑑𝜃. 𝑑𝜑
𝑟2

Por unidad de área de la superficie emisora

ING. ELÍ GUAYÁN HUACCHA Página 9


ING ELI GUAYAN H. BORRADOR RADIACION

𝑑𝑞̇ 𝜆 = 𝐼𝜆,𝑒 . cos 𝜑. sen 𝜑. 𝑑𝜑𝑑𝜃


2𝜋 𝜋 ⁄2
𝑞̇ 𝜆 = ∫0 ∫0 𝐼𝜆,𝑒 . cos 𝜑. sen 𝜑. 𝑑𝜑𝑑𝜃

Relación entre el poder emisivo (E) y la Intensidad de radiación (I):

La energía radiante emitida por dA en la dirección (phi, theta)

De (*): 𝑑𝑄̇𝜆 = 𝐼𝜆 . 𝑑𝐴. cos 𝜑. 𝑑𝜔

Acomodando:

𝑑𝑄̇𝜆 ̇
𝑑𝑄
⟹ = 𝐼𝜆, 𝑒 . sen 𝜑. cos 𝜑. 𝑑𝜃𝑑𝜑; 𝑑𝐴
𝜆
= 𝑑𝐸𝜆
𝑑𝐴

𝑑𝐸𝜆 = 𝐼 sen 𝜑. cos 𝜑. 𝑑𝜃𝑑𝜑

Integrando:

2π π⁄
Eλ = ∫0 ∫0 2I
λ,e sen φ . cos φ. dθdφ.

𝜃 = 0 → 2𝜋 ó 0 → 360𝑜
Limites:{ 𝜑 = 0 → 𝜋 ó 0 → 180𝑜
2

-En superfícies reales: 𝐼 = 𝐼(𝜑)

-En superfícies difusas (ideal):𝐼 = 𝑐𝑡𝑒

Para superfícies difusas:


Integrando:

𝐸𝜆 = 𝜋. 𝐼𝜆 ; 𝐸 = 𝜋. 𝐼
Para superfícies negras: (es emisor difuso: 𝐼𝑏 = 𝑐𝑡𝑒)

𝐸𝑏 = 𝜋. 𝐼𝑏

𝐸𝑏 = 𝜎. 𝑇 4


⇒ Relación: 𝐸 = ∫0 𝐸𝜆 . 𝑑𝜆

ING. ELÍ GUAYÁN HUACCHA Página 10


EGH BORRADOR

TABLAS DE RADIACION

EGH BORRADOR
EGH BORRADOR

EGH BORRADOR
EGH BORRADOR

EGH BORRADOR
ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA


SECCION A

BORRADOR DEL CURSO DE:

TRANSFERENCIA DE CALOR

3RA UNIDAD: LA RADIACION

Mg. ING. ELÍ GUAYÁN H.

TRUJILLO - PERÚ

2021

ING. ELÍ GUAYÁN H. Borrador Página 1


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

Ejemplo: (Borrador)
Considerando que el sol es un cuerpo negro con temperatura exterior 6000 K.
Determinar.
a) La longitud de onda λ para la cual la potencia emisiva espectral es máxima.
b) La fracción de energía radiante que se encuentra en el rango visible (0,38-0,78μm).
Datos:
𝑇 = 6000 °𝐾
𝜆 =?
𝐸𝑏𝜆 (𝑚𝑎𝑥) =?
𝜆𝑚𝑎𝑥

a) De la ley de Wien:
𝜆𝑚𝑎𝑥 ∗ 𝑇 = 2879,8 µm. K
2879,8 µm. K
𝜆𝑚𝑎𝑥 = = 0,483 µm
6000 𝐾

b) Rango:
𝐸𝑏 (0 − 0,38)
𝜆1 𝑇 = 0,38 ∗ 6000 = 2280 µm. K ⟶ = 0,116138
𝐸𝑏
𝐸𝑏 (0 − 0,78)
𝜆2 𝑇 = 0,78 ∗ 6000 = 4680 µm. K ⟶ = 0,590872
𝐸𝑏
La fracción de energía en el rango visible
= 0,590872 − 0,116138 = 0,474734 ó 47,47%
𝐸𝑏 (0,38 − 0,78) = 0,474734 ∗ 𝜎𝑇 4

Ejemplo: (Borrador)
Dos superficies pequeñas de área 0,9 cm2 se encuentran a 12 cm de distancia una de
la otra como en la figura. La superficie 1 es superficie negra y se mantiene a 1000°K.
calcule la velocidad con que la superficie 1 irradia a la superficie 2 (En W)

ING. ELÍ GUAYÁN H. Borrador Página 2


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

Supuesto: Consideramos las pequeñas áreas como elementos diferenciales de áreas

Se sabe que:
Para cuerpo negro el poder emisivo
está relacionado con la intensidad de
Radiación por:
𝐸𝑏 = 𝜋 ∗ 𝐼𝑏
𝐸𝑏 𝜎𝑇 4 (5,67 ∗ 10−8 )(1000)4
𝐼𝑏 = = =
𝜋 𝜋 𝜋
𝑊
𝐼𝑏 = 1,805𝑥104 ∗ 𝑐𝑡𝑒
𝑚2 . 𝑠𝑟
𝑑𝑄̇1−2 = 𝐼1 ∗ 𝑑𝐴1 ∗ cos ∅1 ∗ 𝑑𝜔1−2
El diferencial de ángulo sólido es:
𝑑𝐴𝑁 𝑑𝐴2 ∗ cos ∅2
𝑑𝜔1−2 = =
𝑟2 𝑟2
𝐼1 ∗ 𝑑𝐴1 ∗ cos ∅1 ∗ 𝑑𝐴2 ∗ cos ∅2
𝑑𝑄̇1−2 =
𝑟2
1,805 ∗ 104 ∗ (0,9𝑥10−4 ) ∗ 𝑐𝑜𝑠45° ∗ (0,9𝑥10−4 )𝑐𝑜𝑠60°
𝑑𝑄̇1−2 =
0,122
𝑑𝑄̇1−2 = 35,89 ∗ 10−4 𝑊

2) Ejemplo: La distribucion espectral de la radiacion emitida por una superficie


difusa se puede aproximar como en la figura. ¿Cuál es la potencia emisiva total?.

ING. ELÍ GUAYÁN H. Borrador Página 3


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

SOLUCION:

𝐸𝑏 = ∫0 𝐸𝑏𝜆 𝑑𝜆

Por tramos:

5 10 15 20 ∞
𝐸𝑏 = ∫0 0𝑑𝜆 + ∫5 100𝑑𝜆 + ∫10 200 𝑑𝜆 + ∫15 100𝑑𝜆 + ∫20 0𝑑𝜆

𝐸𝑏 = 100(5) + 200(5) + 100(5) = 500 + 1000 + 500

𝐸𝑏 = 2000 𝑊⁄𝑚2

ING. ELÍ GUAYÁN H. Borrador Página 4


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

LA IRRADIACIÓN (G)

La irradiación es la energía radiante por unidad de tiempo y por unidad de


superficie (receptor) que llega a un cuerpo.

𝑄̇𝑅𝑒𝑎𝑙

A
𝑄̇𝑅𝑒𝑎𝑙 𝑊
𝐺= [ 2]
𝐴 𝑚

La Irradiación (Radiación Incidente, i):

dw
dA

o
Radiación Incidente; Iλ,i .

Intensidad Espectral de Radiación Incidente:

Se define como la proporción en que la energía radiante de longitud de onda λ


incide en la dirección (φ, θ); por unidad de área de la superficie interceptora,
normal a esta dirección por unidad de ángulo sólido alrededor de ésta dirección,
y por intervalo de longitud de onda unitario dλ, alrededor de λ.

Relación entre la Intensidad de Radiación Incidente con la Irradiación G

a) Para la Irradiación Especial (𝐺𝜆 ) :

ING. ELÍ GUAYÁN H. Borrador Página 5


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

2𝜋 𝜋⁄ W
𝐺𝜆 (𝜆) = ∫0 ∫0 2 𝐼 (𝜆, 𝜑, 𝜃) cos 𝜑. sen 𝜑. 𝑑𝜑𝑑𝜃 ; [m2.μm]
𝜆,𝑖

Está definida como la rapidez a la que la radiación de longitud de onda λ, incide


sobre una superficie, por unidad de área de la superficie y por intervalo de
longitud de onda unitario dλ alrededor de λ.

b) Para la Irradiación Integral o Total (G) :

∞ W
G = ∫0 Gλ . dλ [ m2 ]

Convención :
Z
(r,o,o)
o

Y
o

RADIOSIDAD (J)

La radiosidad es la energía radiante por unidad de tiempo y de superficie, en total


evacuada de un cuerpo.

𝐽 = 𝐸 + 𝜌. 𝐺

Dónde:
J=radiosidad

ING. ELÍ GUAYÁN H. Borrador Página 6


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

G=Irradiación o radiación incidente por unidad de tiempo y unidad de área de la


superficie
E=poder de emisión del cuerpo negro
𝜌=reflectividad

Radiosidad Espectral (Iλ ):

2𝜋 𝜋⁄
2 𝑒 = 𝑒𝑚𝑖𝑠𝑖ó𝑛
𝐽𝜆 (𝜆) = ∫ ∫ 𝐼𝜆,(𝑒+𝑟) . (𝜆, 𝜑, 𝜃) cos 𝜑 . sen 𝜑 . 𝑑𝜑. 𝑑𝜃 {
0 0
𝑟 = 𝑟𝑒𝑓𝑙𝑒𝑥𝑖ó𝑛

Radiosidad Total: (𝑱𝑻𝑶𝑻 )



𝑊
𝐽 = ∫ 𝐽𝜆 (𝜆). 𝑑𝜆 ; [ ]
0 𝑚2

EMISIVIDAD (ε) (Borrador)

La emisión de energía radiante en superficies reales es menor que la emisión


de energía radiante de superficies negras (ideales)

Esquema:

𝑄̇ sup real
𝜀=
𝑄̇b (negro)
0 ≤ 𝜀 ≤1
La Emisividad de una Superficie Real ≠ cte; ε=f (T, λ,𝜑)

Esto quiere decir que hay varias emisividades para una superficie dependiendo
de los efectos a considerar.

A) Emisividad Direccional Espectral (ελ, φ):

ING. ELÍ GUAYÁN H. Borrador Página 7


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

A una temperatura dada:

𝐼(𝜆,𝑒) (𝜆, 𝜑, 𝜃,𝑇)


𝜀(𝜆,𝜑) (𝜆, 𝜑, 𝜃, 𝑇) =
𝐼(𝑏,𝜆) (𝜆,𝑇)

Es la razón entre la intensidad de radiación emitida por la superficie a una


longitud de onda específica, en una dirección específica, y la intensidad de
radiación emitida por un cuerpo negro , a la misma temperatura y la misma
longitud de onda.

B) Emisividad Direccional Total (εφ):

𝐼𝑒 (𝜑, 𝜃,𝑇)
𝜀𝜑 (𝜑, 𝜃, 𝑇) =
𝐼𝑏 (𝑇)

Se define como el cociente de intensidades de radiación emitida por un cuerpo


real y por un cuerpo negro, es decir usando intensidades totales de radiaciones
(integradas sobre totas las longitudes de onda).

Nota: En la práctica suele ser más conveniente trabajar con propiedades


relativas a la radiación promediadas sobre todas las direcciones y estas son las
llamadas propiedades hemisféricas.

C) La Emisividad Hemisférica Espectral (ελ):

𝐸𝜆 (𝜆,𝑇)
𝜀𝜆 (𝜆, 𝑇) =
𝐸(𝑏,𝜆) (𝜆,𝑇)

𝑬𝝀 : Potencia o poder emisivo espectral del cuerpo (real)

𝑬𝒃𝝀 : Potencia o poder emisivo espectral del cuerpo negro

d) Emisividad Hemisférica Total:

𝐸(𝑇)
𝜀(𝑇) =
𝐸𝑏 (𝑇)

Eb (T)= σT4: Ley de Stefan Bolzman o integral de la ley Plank.

Se define en términos de la energía de radiación emitida sobre todas las


longitudes de onda y en todas las direcciones posibles.

ING. ELÍ GUAYÁN H. Borrador Página 8


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

La emisividad hemisférica total, de una superficie a una temperatura dada


representa la razón entre la energía total de radiación emitida por la superficie y
la radiación emitida por un cuerpo negro de la misma área superficial a la misma
temperatura.

Tenemos:

𝐸𝜆 = 𝜀𝜆 𝐸𝑏𝜆
Total:

𝐸 = ∫ 𝐸𝜆 𝑑𝜆
0

∞ ∞
𝐸(𝑇) ∫0 𝜀𝜆 (𝜆,𝑇)𝐸𝑏𝜆 (𝜆, 𝑇)𝑑𝜆 ∫0 𝜀𝜆 𝐸𝑏𝜆 𝑑𝜆
𝜀(𝑇) = = ∞ =
𝐸𝑏 (𝑇) ∫0 𝜀𝑏𝜆 𝑑𝜆 𝜎𝑇 4

NOTA:
Superficie difusa: sus propiedas son independientes de su direccion.
Superficie Gris: Sus propiedades son independientes de la longitud de onda.

ING. ELÍ GUAYÁN H. Borrador Página 9


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

POR TRAMOS:
Para la emisividad espectral modelada por tramos:
𝜀1 = 𝑐𝑡𝑡𝑒 ,0 ≤ 𝜆 < 𝜆1
Si : 𝜀𝜆 = { 𝜀2 = 𝑐𝑡𝑡𝑒 , 𝜆1 ≤ 𝜆 < 𝜆2
𝜀3 = 𝑐𝑡𝑡𝑒, 𝜆2 ≤ 𝜆 < ∞

𝜆 𝜆 ∞
𝜀1 ∫0 1 𝐸𝑏𝜆 𝑑𝜆 𝜀2 ∫𝜆 2 𝐸𝑏𝜆 𝑑𝜆 𝜀3 ∫𝜆 𝐸𝑏𝜆 𝑑𝜆
1 2
𝜀(𝑇) = + +
𝐸𝑏 𝐸𝑏 𝐸𝑏

Relación util:

𝜋/2
𝜀𝜆 (𝜆 , 𝑇 ) = 2 ∫ 𝜀𝜆,𝜑 (𝜆 , 𝜑, 𝑇)𝑠𝑒𝑛𝜑 𝑐𝑜𝑠𝜑𝑑𝜑
0

EJERCICIOS DE APLICACIÓN:

1)¿Cuál es la emisividad hemisferica total cuando la temperatura de un filamento


de tugnsteno es 2900K?

ING. ELÍ GUAYÁN H. Borrador Página 10


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

SOLUCION:

𝐸(𝑇) ∫0 𝜀𝜆 𝐸𝑏𝜆 𝑑𝜆
𝜀(𝑇) = 𝐸 =
𝑏 (𝑇) 𝐸𝑏

En partes (por tramos):


2 ∞
𝜀1 ∫0 𝐸𝑏𝜆 𝑑𝜆 𝜀2 ∫2 𝐸𝑏𝜆 𝑑𝜆
𝜀𝑇𝑜𝑡𝑎𝑙 = +
𝐸𝑏 𝐸𝑏

𝜀1 𝐸𝑏 (0−2) 𝐸𝑏 (0−∞) 𝐸𝑏 (0−2)


𝜀𝑇𝑜𝑡𝑎𝑙 = +[ +− ]
𝐸𝑏 𝐸𝑏 𝐸𝑏

λT Tablas

0(2900) 0

2(2900) 0,720158

𝜀𝑇𝑜𝑡𝑎𝑙 = 0,45(0,720158) + 0.10[1 − 0.720158]

Rpta: 𝜺𝑻𝒐𝒕𝒂𝒍 = 𝟎, 𝟑𝟓𝟐𝟎𝟓𝟓𝟑

ABSORTIVIDAD, REFLECTIVIDAD Y TRANSMISIVIDAD (Borrador)


Esquema:

ING. ELÍ GUAYÁN H. Borrador Página 11


ING ELI GUAYAN H. BORRADOR SECCION A RADIACION

𝐺𝑎𝑏𝑠 𝐺𝑟𝑒𝑓 𝐺𝑡𝑟𝑎𝑛𝑠


𝛼= ; 𝜌= ; 𝜏=
𝐺𝑇 𝐺𝑇 𝐺𝑇
Por balance:

𝐺𝑎𝑏𝑠 + 𝐺𝑟𝑒𝑓 + 𝐺𝑡𝑟𝑎𝑛𝑠 = 𝐺𝑇𝑂𝑇

𝛼. 𝐺𝑇 + 𝜌. 𝐺𝑇 + 𝜏. 𝐺𝑇 = 𝐺𝑇
∴ 𝛼+𝜌+𝜏 =1
Para superficie “opaca”:
𝜏=0
→ 𝛼+𝜌≅1
ABSORTIVIDAD ESPECTRAL (𝜶𝝀 ):

𝐺𝜆 ∫0 𝛼𝜆 . 𝐺𝜆 . 𝑑𝜆
𝛼𝜆 = 𝑎𝑏𝑠 → 𝑇𝑜𝑡𝑎𝑙: 𝛼̅ = ∞
𝐺𝜆 ∫0 𝐺𝜆 . 𝑑𝜆

REFLECTIVIDAD ESPECTRAL (𝝆𝝀 ):



𝐺𝜆𝑟𝑒𝑓𝑙𝑒𝑗 ∫0 𝜌𝜆 . 𝐺𝜆 . 𝑑𝜆
𝜌𝜆 = → 𝑇𝑜𝑡𝑎𝑙: 𝜌̅ = ∞
𝐺𝜆 ∫0 𝐺𝜆 . 𝑑𝜆
TRANSMISIVIDAD ESPECTRAL (𝝉𝝀 ):

𝐺𝜆 ∫0 𝜏𝜆 . 𝐺𝜆 . 𝑑𝜆
𝜏𝜆 = 𝑡𝑟𝑎𝑛𝑠 → 𝑇𝑜𝑡𝑎𝑙: 𝜏̅ = ∞
𝐺𝜆 ∫0 𝐺𝜆 . 𝑑𝜆

LEY DE KIRCHOFF:
Total: 𝜀(𝑇) = 𝛼(𝑇)

Espectral: 𝜀𝜆 = 𝛼𝜆

ING. ELÍ GUAYÁN H. Borrador Página 12


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA

BORRADOR DEL CURSO DE:

TRANSFERENCIA DE CALOR
SECCION A

3RA UNIDAD: LA RADIACION

Mg. ING. ELÍ GUAYÁN H.

TRUJILLO - PERÚ

2021

ING. ELÍ GUAYÁN H. BORRADOR Página 1


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

LEY DE KIRCHHOFF: (Borrador)


A1<<A2

Recinto cerrado Isotérmico (Cavidad de cuerpo negro)

Considere un pequeño cuerpo, de emisividad e y coeficiente alpha, a la temperatura T,


encerrado en un recinto cerrado interno de temperatura T.
La irradiación G que experimenta el pequeño cuerpo en la cavidad difusa, es difuso e
igual a la emitida por un cuerpo negro a la temperatura T.
Y en condiciones de estado estable (Equilibrio térmico), la energía absorbida sería
igual a la energía emitida.
𝐺 = 𝐸𝑏 (𝑇)

𝐺𝑎𝑏𝑠 𝐴1 = 𝐸𝑒𝑚𝑖𝑡 𝐴1

𝛼 𝐺 𝐴1 = 𝜀 𝐴1 𝜎𝑇 4
𝛼𝐸𝑏 𝐴1 = 𝜀 𝐴1 𝜎𝑇 4
𝛼 (𝜎𝑇 4 ) 𝐴1 = 𝜀 𝐴1 𝜎𝑇 4

Total: 𝜀(𝑇) = 𝛼(𝑇)

Espectral: 𝜀𝜆 = 𝛼𝜆

ING. ELÍ GUAYÁN H. BORRADOR Página 2


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

EL FACTOR DE VISION Y EL INTERCAMBIO DE ENERGIA RADIANTE ENTRE


SUPERFICIES NEGRAS FINITAS (Borrador)
Esquema:
Con: 𝑇1 > 𝑇2

Donde: 𝑟̅ = vector
𝑛̅1 = normal (1)
𝑛̅2 = normal (2)
Supuestos:
-Radiación emitida difusa
-Superficies separadas por un medio no absorbente
-Superficies negras
El calor neto intercambiado por (1) y (2) por unidad de tiempo:
𝑄̇1,2 = 𝑄̇1−2 − 𝑄̇2−1 "𝑙𝑜 𝑞𝑢𝑒 𝑠𝑎𝑙𝑒 − 𝑙𝑜 𝑞𝑢𝑒 𝑙𝑙𝑒𝑔𝑎"

De (1) a (2): en forma diferencial (y usando la definición de intensidad)


𝑑𝐴2 . cos 𝜑2
𝑑𝑄̇1−2 = 𝐼1 . 𝑑𝐴
⏟ 1 . cos 𝜑1 . 𝑑𝜔2−1 = 𝐼1 . 𝑑𝐴1 . cos 𝜑 . ⏟ 𝑟2
á𝑟𝑒𝑎 𝑝𝑟𝑜𝑦𝑒𝑐𝑡.
𝑑𝑤2−1

𝐼1 . cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


𝑑𝑄̇1−2 =
𝑟2
Integrando:
cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2
𝑄̇1−2 = 𝐼1 . ∫ ∫
𝐴1 𝐴2 𝑟2

Dividiendo entre (𝐴1 . 𝐸𝑏1 ):

𝑄̇1−2 𝐼1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


= .∫ ∫
𝐴1 . 𝐸𝑏1 𝐴1 . 𝐸𝑏1 𝐴1 𝐴2 𝑟2

Y tenemos que : 𝐸𝑏1 = 𝜋. 𝐼1 (radiación de cuerpo negro)

ING. ELÍ GUAYÁN H. BORRADOR Página 3


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

𝑄̇1−2 𝐼1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


= .∫ ∫
𝐴1 . 𝐸𝑏1 𝐴1 (𝜋. 𝐼1 ) 𝐴1 𝐴2 𝑟2

𝑄̇1−2 1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


= .∫ ∫ → (𝑎𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙)
𝐴1 . 𝐸𝑏1 𝐴1 𝐴1 𝐴2 𝜋. 𝑟 2

1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


𝐹1−2 = .∫ ∫
𝐴1 𝐴1 𝐴2 𝜋. 𝑟 2

F1-2: Es el factor de forma o factor de visión.


Indica la fracción de energía radiante que sale de la superficie (1) y que es
interceptada por la superficie (2).

𝑄̇1−2 = 𝐴1 . 𝐸𝑏1 . 𝐹1−2 → [𝑊]

𝐸𝑏1 = 𝜎. 𝑇 4

Para la superficie (2): 𝑄̇2−1 = 𝐴2 . 𝐸𝑏2 . 𝐹2−1


1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2
𝐹2−1 = .∫ ∫
𝐴2 𝐴1 𝐴2 𝜋. 𝑟 2

F2-1: Factor de forma de 2-1.


Ahora, Igualando los integrales:
𝐴1 . 𝐹1−2 = 𝐴2 . 𝐹2−1
Se llama relación de reciprocidad.
El calor neto intercambiado por (1) y (2):

𝑄̇1,2 = 𝐴1 . 𝐸𝑏1 . 𝐹1−2 − 𝐴2 . 𝐸𝑏2 . 𝐹2−1

𝑄̇1,2 = 𝐴1 . 𝐸𝑏1 . 𝐹1−2 − 𝐸𝑏2 (𝐴1 . 𝐹1−2 )

𝑄̇1,2 = 𝐴1 . 𝐹1−2 (𝐸𝑏1 − 𝐸𝑏2 )

𝑄̇1,2 = 𝐴1 . 𝐹1−2 . 𝜎[𝑇14 − 𝑇24 ] → [𝑊]

FACTOR DE FORMA PARA DISCOS PARALELOS COAXIALES:

ING. ELÍ GUAYÁN H. BORRADOR Página 4


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

𝑟1
Relaciones: 𝑅1 = → 𝑎𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙
𝐿
𝑟2
𝑅2 = → 𝑎𝑑𝑖𝑚𝑒𝑛𝑠𝑖𝑜𝑛𝑎𝑙
𝐿

1+𝑅22
𝑆 =1+ 𝑅12

Factor:
1⁄
1 𝑅2 2 2
𝐹1−2 = {𝑆 − [𝑆 2 − 4 ( ) ] }
2 𝑅1

Ejemplo:
Dos discos concéntricos paralelos con diámetro D1=10 cm y diámetro D2=5 cm
distancia entre discos L= 10 cm uno del otro
a) Determine el factor de forma de 1 a 2
b) Si la superficie 1 está a 600°K y la superficie 2 está a 300°K ¿Cuál es el calor
emitido por 1 y que alcanza al cuerpo 2, suponiendo cuerpos negros?
Solución.
a)
𝑑1 = 10 𝑐𝑚 ⟶ 𝑟1 = 5 𝑐𝑚
𝑑2 = 5 𝑐𝑚 ⟶ 𝑟2 = 2,5 𝑐𝑚
𝐿 = 10 𝑐𝑚
𝑟1 𝑟2
𝑅1 = 𝐿
= 0,5 𝑅2 = 𝐿
= 0,25

1 + 0,252
𝑆 =1+ = 5,25
0,52
1 0,25 2 1/2
𝐹1−2 = [5,25 − (5.252 − 4( ) )]
2 0,5
𝐹1−2 = 0,0481 ⟶ 4,81%
b)

ING. ELÍ GUAYÁN H. BORRADOR Página 5


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

𝑑1 2
𝑄̇1−2 = 𝐴1 ∗ 𝐹1−2 ∗ 𝐸𝑏1 = 𝜋 ∗ 𝐹1−2 ∗ 𝜎𝑇1 4
4
(10 ∗ 10−2 )2
𝑄̇1−2 = 𝜋 ∗ 0,0481 ∗ 5,67 ∗ 10−4 ∗ 6004
4
𝑄̇1−2 = 2,776 𝑊

BALANCE DE RADIACION PARA SUPERFICIE CONCAVA CERRADA (INTERNA)


Sea un recinto de n partes

Por balance:
Para la superficie (1):

𝐴⏟1 . 𝐸𝑏1 = 𝐴
⏟1 . 𝐸𝑏1 . 𝐹1−1 + 𝐴
⏟1 . 𝐸𝑏1 . 𝐹1−2 + 𝐴1 . 𝐸𝑏1 . 𝐹1−3 + ⋯ + 𝐴1 . 𝐸𝑏1 . 𝐹1−𝑛
𝑡𝑜𝑑𝑜 𝑙𝑜 𝑞𝑢𝑒 𝑙𝑎 𝑝𝑎𝑟𝑡𝑒 𝑞𝑢𝑒 𝑙𝑎 𝑝𝑎𝑟𝑡𝑒 𝑞𝑢𝑒
𝑠𝑎𝑙𝑒 𝑑𝑒 (1) 𝑙𝑙𝑒𝑔𝑎 𝑎 (1) 𝑙𝑙𝑒𝑔𝑎 𝑎 (2)

1 = 𝐹1−1 + 𝐹1−2 + 𝐹1−3 + ⋯ + 𝐹1−𝑛

Para la superficie (2):


𝐴2 . 𝐸𝑏2 = 𝐴2 . 𝐸𝑏2 . 𝐹2−1 + 𝐴2 . 𝐸𝑏2 . 𝐹2−2 + 𝐴2 . 𝐸𝑏2 . 𝐹2−3 + ⋯ + 𝐴2 . 𝐸𝑏2 . 𝐹2−𝑛
1 = 𝐹2−1 + 𝐹2−2 + 𝐹2−3 + ⋯ + 𝐹2−𝑛

Por tanto, tenemos:


𝐹1−1 + 𝐹1−2 + 𝐹1−3 + ⋯ + 𝐹1−𝑛 = 1
𝐹2−1 + 𝐹2−2 + 𝐹2−3 + ⋯ + 𝐹2−𝑛 = 1
𝐹3−1 + 𝐹3−2 + 𝐹3−3 + ⋯ + 𝐹3−𝑛 = 1
.
.
.
{𝐹𝑛−1 + 𝐹𝑛−2 + 𝐹𝑛−3 + ⋯ + 𝐹𝑛−𝑛 = 1

ING. ELÍ GUAYÁN H. BORRADOR Página 6


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

Ejemplo:
Un tronco de cono tiene un diámetro 𝑑2 = 12 𝑐𝑚 𝑦 𝑑1 = 24 𝑐𝑚 y una altura de 12 cm.
a) Calcular el factor de forma entre la capa superior n°2 y la superficie lateral
b) Calcule el factor de forma de la superficie lateral consigo misma.

Factor de forma de 1 a 2
Relaciones:
𝑟1 12
𝑅1 = = =1
𝐿 12
𝑟2 6
𝑅2 = = = 0,5
𝐿 12
1 + 𝑅2 2 1 + 0,52
𝑆 =1+ =1+ = 2,25
𝑅1 12
1
1 2
0,5 2 2
𝐹1−2 = [2,25 − (2,25 − 4 ∗ ( ) ) ]
2 1

𝐹1−2 = 0,117217
Factor de forma de 2 a 1
De la relación de reciprocidad.
𝐴1 ∗ 𝐹1−2 = 𝐴2 ∗ 𝐹2−1

𝐴1 ∗ 𝐹1−2 𝜋(0,242 /4) ∗ 0,117217


𝐹2−1 = = = 0,468868
𝐴2 𝜋(0,122 /4)

Balance superficie 2

1 = 𝐹2−1 + 𝐹2−2 + 𝐹2−3

Pero 𝐹2−2 = 0 (Superficie plana perfecta)

1 = 𝐹2−1 + 𝐹2−3

1 = 0,468868 + 𝐹2−3

ING. ELÍ GUAYÁN H. BORRADOR Página 7


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

𝐹2−3 = 0,531132

Balance superficie 3

1 = 𝐹3−1 + 𝐹3−2 + 𝐹3−3

Balance superficie 1

1 = 𝐹1−1 + 𝐹1−2 + 𝐹1−3

Pero: 𝐹1−1 = 0

𝐹1−3 = 1 − 0,117217 = 0,882783

Factor de forma de 3 a 1:

Por reciprocidad
𝐴1 ∗ 𝐹1−3 = 𝐴3 ∗ 𝐹3−1
𝐴1
𝐹3−1 = ∗𝐹
𝐴3 1−3

de la figura:

g: generatriz

𝐴𝐿 = 𝜋 ∗ 𝑔(𝑅 + 𝑟)

𝑔 = √𝐿2 + (𝑅 − 𝑟)2

𝐴3 = 𝜋(𝑅1 + 𝑅2 )[(𝑅1 − 𝑅2 )2 + 𝐿2 ]0,5

𝐴3 = 𝜋(12 + 6)[(12 − 6)2 + 122 ]0,5

𝐴3 = 758,678 𝑐𝑚^2

𝜋 ∗ 242
𝐹3−1 = ∗ 0,882783 = 0,52639 0,531132 = 0,0791765
4 ∗ 758,678

Reciprocidad entre 2 y 3:
𝐴2 ∗ 𝐹2−3 = 𝐴3 ∗ 𝐹3−2

𝐴2 ∗ 𝐹2−3
𝐹3−2 =
𝐴3

De la ecuación de balance para la superficie 3

𝐹3−3 = 1 − 𝐹3−1 − 𝐹3−2

𝐹3−3 = 1 − 0,52639 − 0,0791765

ING. ELÍ GUAYÁN H. BORRADOR Página 8


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

𝐹3−3 = 0,394433

ING. ELÍ GUAYÁN H. BORRADOR Página 9


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA

BORRADOR DEL CURSO DE:

TRANSFERENCIA DE CALOR
SECCION A

3RA UNIDAD: LA RADIACION

Mg. ING. ELÍ GUAYÁN H.

TRUJILLO - PERÚ

2021

ING. ELÍ GUAYÁN H. BORRADOR Página 1


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

Ejercicio: (Borrador)
Determinar el factor de forma geométrica para un disco A1 muy pequeño y un
disco paralelo grande A2, colocado a una distancia L directamente arriba del más
pequeño.
De la fórmula del factor de forma:

1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


𝐹1−2 = .∫ ∫
𝐴1 𝐴1 𝐴2 𝜋. 𝑟 2

Consideración:

A1  menor, A1
dA1 = A1

𝐴1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴2


𝐹1−2 = .∫
𝐴1 𝐴2 𝜋. 𝑟 2

𝐿
Para el disco grande cos 𝜑1 = cos 𝜑2 =
𝑟

r =  2 + L2

𝑑𝐴2 = 𝜌 𝑑𝜃 𝑑𝜌
𝐴1 𝐿2 𝜌 𝑑𝜃 𝑑𝜌
𝐹1−2 = .∫ 2
𝐴1 𝐴2 𝑟 . 𝜋. 𝑟 2

ING. ELÍ GUAYÁN H. BORRADOR Página 2


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

𝑎
2𝜋
𝐴1 𝐿2 𝜌 𝑑𝜃 𝑑𝜌 𝐴1 𝐿2 𝑎 2𝜋𝜌𝑑𝜌
𝐹1−2 = ∫∫ 2 2 2
= ∫
𝐴1 0 𝜋(𝜌 + 𝐿 ) 𝐴1 𝜋 0 (𝜌2 + 𝐿2 )2
0

𝐴1 2𝜋. 𝐿2 𝑎 𝜌𝑑𝜌
𝐹1−2 = ∫ 2 2 2
𝐴1 𝜋 0 (𝜌 + 𝐿 )

𝐴1 1 𝑎
𝐹1−2 = . 2𝐿2 . ∫ (𝜌2 + 𝐿2 )−2 2𝜌𝑑𝜌
𝐴1 2 0

𝑎
𝐴1 2 (𝜌2 + 𝐿2 )−1
𝐹1−2 = .𝐿 . |
𝐴1 −1 0
𝑎
−𝐴1 2 1 −𝐴1 2 1 1
𝐹1−2 = .𝐿 .( 2 2
)| = .𝐿 [ 2 2
− 2]
𝐴1 𝜌 +𝐿 0 𝐴1 𝑎 +𝐿 𝐿
−𝐴1 2 𝐿2 − 𝑎2 − 𝐿2
𝐹1−2 = .𝐿 [ 2 ]
𝐴1 (𝑎 + 𝐿2 )(𝐿2 )

𝑎2
𝐹1−2 =
𝑎2 + 𝐿2

Rpta.

ING. ELÍ GUAYÁN H. BORRADOR Página 3


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

Configuraciones complejas reducibles a rectángulos perpendiculares con un


borde común: (Borrador)
1)
𝐸𝑛𝑒𝑟𝑔í𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑞𝑢𝑒 𝐸𝑛𝑒𝑟𝑔í𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑞𝑢𝑒 𝐸𝑛𝑒𝑟𝑔í𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑞𝑢𝑒
= +
𝑠𝑎𝑙𝑒 𝑑𝑒 (1) ℎ𝑎𝑐𝑖𝑎 (2 + 3) 𝑠𝑎𝑙𝑒 𝑑𝑒 (1) ℎ𝑎𝑐𝑖𝑎 (3) 𝑠𝑎𝑙𝑒 𝑑𝑒 (1) ℎ𝑎𝑐𝑖𝑎 (2)

𝐴1 . 𝐸𝑏1 . 𝐹(1)−(2+3) = 𝐴1 . 𝐸𝑏1 . 𝐹(1)−(3) + 𝐴1 . 𝐸𝑏1 . 𝐹(1)−(2)

𝐹(1)−(2+3) = 𝐹(1)−(3) + 𝐹(1)−(2)

Despejando:
𝐹(1)−(2) = 𝐹(1)−(2+3) − 𝐹(1)−(3)

Si queremos 𝐹(2)−(1) =?

Aplicar reciprocidad: 𝐴1 . 𝐹1−2 = 𝐴2 . 𝐹2−1

2)
𝐸𝑛𝑒𝑟𝑔í𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑞𝑢𝑒 𝐸𝑛𝑒𝑟𝑔í𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑞𝑢𝑒 𝐸𝑛𝑒𝑟𝑔í𝑎 𝑟𝑎𝑑𝑖𝑎𝑛𝑡𝑒 𝑞𝑢𝑒
= +
𝑠𝑎𝑙𝑒 𝑑𝑒 (1 + 3) ℎ𝑎𝑐𝑖𝑎 (2) 𝑠𝑎𝑙𝑒 𝑑𝑒 (1) ℎ𝑎𝑐𝑖𝑎 (2) 𝑠𝑎𝑙𝑒 𝑑𝑒 (3) ℎ𝑎𝑐𝑖𝑎 (2)
𝐴(1+3) . 𝐸𝑏 . 𝐹(1+3)−(2) = 𝐴1 . 𝐸𝑏 . 𝐹(1)−(2) + 𝐴3 . 𝐸𝑏 . 𝐹(3)−(2)

𝐴(1+3) . 𝐹(1+3)−(2) = 𝐴1 . 𝐹(1)−(2) + 𝐴3 . 𝐹(3)−(2)

Despejando:
𝐴1 . 𝐹(1)−(2) = 𝐴(1+3) . 𝐹(1+3)−(2) − 𝐴3 . 𝐹(3)−(2)
𝐴(1+3) . 𝐹(1+3)−(2) − 𝐴3 . 𝐹(3)−(2)
→ 𝐹(1)−(2) =
𝐴1

3)
Hallar: 𝐹1−2 =? ?

ING. ELÍ GUAYÁN H. BORRADOR Página 4


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

Por balance:
𝐴(1+3) . 𝐸𝑏 . 𝐹(1+3)−(2+4) = 𝐴1 . 𝐸𝑏 . 𝐹(1)−(2+4) + 𝐴3 . 𝐸𝑏 . 𝐹(3)−(2+4)

𝐴(1+3) . 𝐹(1+3)−(2+4) = 𝐴
⏟1 . 𝐹(1)−(2+4) + 𝐴3 . 𝐹(3)−(2+4) → (1)

Tenemos que:
Desarrollo de (*): 𝐴1 . 𝐹(1)−(2+4) = 𝐴1 . 𝐹(1)−(2) + ⏟
𝐴1 . 𝐹(1)−(4) → (2)

También: 𝐴(1+3) . 𝐹(1+3)−(4) = ⏟


𝐴1 . 𝐹(1)−(4) + 𝐴3 . 𝐹(3)−(4) → (3)
∗∗

Despejando (3): 𝐴1 . 𝐹(1)−(4) = 𝐴(1+3) . 𝐹(1+3)−(4) − 𝐴3 . 𝐹(3)−(4)

Y reemplazando en (2):
𝐴1 . 𝐹(1)−(2+4) = 𝐴1 . 𝐹(1)−(2) + 𝐴(1+3) . 𝐹(1+3)−(4) − 𝐴3 . 𝐹(3)−(4)

Reemplazando en (1):
𝐴(1+3) . 𝐹(1+3)−(2+4) = 𝐴1 . 𝐹(1)−(2) + 𝐴(1+3) . 𝐹(1+3)−(4) − 𝐴3 . 𝐹(3)−(4) + 𝐴3 . 𝐹(3)−(2+4)

1
∴ 𝐹(1)−(2) = [𝐴( ) . 𝐹(1+3)−(2+4) + 𝐴3 . 𝐹(3)−(4) − 𝐴(1+3) . 𝐹(1+3)−(4) − 𝐴3 . 𝐹(3)−(2+4) ]
𝐴1 1+3

METODO DEL ASPA:


se cumple:
𝐴1 . 𝐹1−3 = 𝐴2 . 𝐹2−4

ING. ELÍ GUAYÁN H. BORRADOR Página 5


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

Ejemplo:

Hallar el factor de forma 𝐹1−2 para la siguiente configuración.

Por balance.

𝐹1−2 = 𝐹1−(2+3) − 𝐹1−3

a) Factor de forma de 1 a (2+3)

4
𝑅1 = =1
4

De la gráfica el factor de forma 1 a (2+3)

𝐹1−(2+3) = 0,2

Hallamos 𝐹1−3
4
𝑅1 = = 1
4

2
𝑅2 = 4 = 0,5

De la gráfica el factor de forma de 1 a 3. 𝐹1−3 = 0,144

Por balance.

𝐹1−2 = 𝐹1−(2+3) − 𝐹1−3

𝐹1−2 = 0,2 − 0,144 = 0,056

b) Si 𝑇1= 500𝐾 y es cuerpo negro. Hallar el calor por radiación de 1 a 2.

𝑄̇1−2 = 𝐴1 ∗ 𝐹1−2 ∗ 𝐸𝑏1 = 16 ∗ 0,056 ∗ 5,67 ∗ 10−8 ∗ 5004

𝑄̇1−2 = 3175,2𝑊

ING. ELÍ GUAYÁN H. BORRADOR Página 6


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

Ejemplo:
Para la siguiente configuración hallar el factor de forma F1-2 ?

Supuesto: T1=T3
Eb1=Eb3
Hallar: F1-2
Solución: del balance
𝐴(1+3) . 𝐹(1+3)−(2) − 𝐴3 . 𝐹(3)−(2)
𝐹(1)−(2) =
𝐴1
-Calculando el factor de forma de (1+3) a (2)

De la gráfica:
𝐹(1+3)−(2) = 0,2
-Calculando el factor deforma de 3 a 2

De la gráfica:
𝐹3−2 = 0,29
Luego:

4 ∗ 4 ∗ 0,2 − 2 ∗ 4 ∗ 0,29
𝐹1−2 = = 0,11
2∗4

ING. ELÍ GUAYÁN H. BORRADOR Página 7


ING ELI GUAYAN H. BORRADOR RADIACION SECCION A

b) Si 𝑇1= 500𝐾 y es cuerpo negro. Hallar el calor por radiación de 1 a 2.

𝑄̇1−2 = 𝐴1 ∗ 𝐸𝑏1 ∗ 𝐹1−2 = 8 ∗ 5,67 ∗ 10−8 ∗ 5004 ∗ 0,11

𝑄̇1−2 = 3118,5 𝑊

Tarea de cuaderno:
Para la configuración hallar el factor de forma F1-2?

Solución:
Supuesto:
T1=T3
Eb1=Eb3
Por simetría
𝐹1−2 = 𝐹3−4 ; 𝐹3−2 = 𝐹1−4

A3=A1
Balance:

ING. ELÍ GUAYÁN H. BORRADOR Página 8


ING ELI GUAYAN H. BORRADOR RADIACION

UNIVERSIDAD NACIONAL DE TRUJILLO

FACULTAD DE INGENIERÍA

ESCUELA DE INGENIERÍA MECÁNICA

BORRADOR DEL CURSO DE:

TRANSFERENCIA DE CALOR
SECCION A

3RA UNIDAD: LA RADIACION

Mg. ING. ELÍ GUAYÁN H.

TRUJILLO - PERÚ

2021

ING. ELÍ GUAYÁN H. BORRADOR Página 1


ING ELI GUAYAN H. BORRADOR RADIACION

Factor de forma para dos rectangulos prependiculares (Borrador)

1 cos 𝜑1 . cos 𝜑2 . 𝑑𝐴1 . 𝑑𝐴2


𝐹1−2 = .∫ ∫
𝐴1 𝐴1 𝐴2 𝜋. 𝑟 2

𝑑𝐴1 = 𝑑𝑥1 𝑑𝑦1


𝑑𝐴2 = 𝑑𝑦2 𝑑𝑧2

𝑎 = √(𝑦1 − 𝑦2 )2 + 𝑥12 𝑟 2 = 𝑎2 + 𝑧22


𝑧2 𝑥1
cos 𝜑1 = cos 𝜑2 =
𝑟 𝑟

𝐿 𝐷 𝐷 𝑊
1 𝑥1 𝑧2 𝑑𝑥1 𝑑𝑦1 𝑑𝑦2 𝑑𝑧2
𝐹1−2 = ∫ ∫ ∫ ∫ 2
𝐷. 𝐿 0 0 0 0 2
𝜋 [(𝑦1 − 𝑦2 ) + 𝑥21 + 𝑧22 ]

ING. ELÍ GUAYÁN H. BORRADOR Página 2


ING ELI GUAYAN H. BORRADOR RADIACION

TRANSFERENCIA DE CALOR POR RADIACIÓN: SUPERFICIES GRISES Y


DIFUSAS (Borrador)
Radiosidad:J

Radiosidad (sale)

𝐽 = 𝐸 + 𝜌. 𝐺
2
J= Radiosidad (W/m )
E= Potencia emisiva (W/m2)
G=Irradiación (W/m2)
𝑄̇
𝐽 = 𝜀. 𝐸𝑏 + 𝜌. 𝐺 ; 𝑞̇ =
𝐴
𝜀 = Emisividad
Si: 𝜌+𝛼 =1; 𝜏≅0

𝐽 = 𝜀. 𝐸𝑏 + (1 − 𝛼). 𝐺
De la ley de Kirchhoff: 𝛼=𝜀

𝐽 = 𝜀. 𝐸𝑏 + (1 − 𝜀). 𝐺 …. (1)
FLUJO NETO DE CALOR POR RADIACION DESDE UNA SUPERFICIE
Haciendo balance de energía en la superficie:
(Radiación que sale – radiación que incide)

𝑄̇1 = 𝐴1 . 𝐽1 − 𝐴1 . 𝐺1

𝑄̇1 = 𝐴1 ( 𝐽1 − 𝐺1 ) … . . (2)
Despejando G de la ecuación (1) y remplazando en la ecuación (2):
𝐽1 − 𝜀1 𝐸𝑏1
𝑄̇1 = 𝐴1 ( 𝐽1 − )
1 − 𝜀1
𝐴1 𝜀1
𝑄̇1 = ( 𝐸 − 𝐽1 )
1 − 𝜀1 𝑏1
Acomodando:
( 𝐸𝑏1 − 𝐽1 )
𝑄̇1 =
1 − 𝜀1
𝐴1 𝜀1

ING. ELÍ GUAYÁN H. BORRADOR Página 3


ING ELI GUAYAN H. BORRADOR RADIACION

Por analogía eléctrica:


∆𝑉
𝐼=
𝑅𝑒
El numerador (𝐸𝑏 − 𝐽) se puede considerar como una diferencia de potencial.
1−𝜀
El denominador ( ) se puede considerar como una resistencia superficial a la
𝜀.𝐴
radiación.
“Resistencia de la superficie” a la radiación
1 − 𝜀1
𝑅1 =
𝐴1 𝜀1
Esquema:

1 − 𝜀𝑖 1−𝜀
𝑅𝑖 = 𝑅=
𝐴𝑖 . 𝜀𝑖 𝜀. 𝐴

FLUJO DE CALOR NETO ENTRE DOS SUPERFICIES:


Para dos superficies difusas grises y opacas:

Radiación Neta: rad. que sale – rad. que llega


𝑄̇1,2 = 𝐴1 . 𝐽1 . 𝐹1−2 − 𝐴2 . 𝐽2 . 𝐹2−1
Por la relación de reciprocidad:

𝐴1 . 𝐹1−2 = 𝐴2 . 𝐹2−1

𝑄̇1,2 = 𝐴1 . 𝐹1−2 (𝐽1 − 𝐽2 )


Acomodando:

ING. ELÍ GUAYÁN H. BORRADOR Página 4


ING ELI GUAYAN H. BORRADOR RADIACION

𝐽1 − 𝐽2
𝑄̇1,2 =
1
𝐴1 . 𝐹1−2
Por analogía eléctrica:

La “resistencia del espacio” a la radiación


1
𝑅1−2 = 𝐴 } Resistencia espacial
1 .𝐹1−2

PARA UN RECINTO DE 𝑁 SUPERFICIES:


La transferencia neta de calor desde la superficie i sea igual a la suma de transferencias
netas de calor desde la superficie i hacia cada una de las N superficies del recinto:

𝑁 𝑁 𝑁
𝐽𝑖 − 𝐽𝑗
𝑄̇𝑖 = ∑ 𝑄̇𝑖,𝑗 = ∑ 𝐴𝑖 . 𝐹𝑖−𝑗 (𝐽𝑖 − 𝐽𝑗 ) = ∑
𝑅𝑖−𝑗
𝑗=1 𝑗=1 𝑗=1

RECINTO CERRADO DE DOS SUPERFICIES (ENVOLTURAS GRISES DE DOS


COMPONENTES):

ING. ELÍ GUAYÁN H. BORRADOR Página 5


ING ELI GUAYAN H. BORRADOR RADIACION

Esquema:

1 − 𝜀1 1 1 − 𝜀2
𝑅1 = ; 𝑅1,2 = ; 𝑅2 =
𝜀1 . 𝐴1 𝐴1 . 𝐹1−2 𝜀2 . 𝐴2
Superficie Espacio Superficie

𝐸𝑏1 − 𝐸𝑏2 𝐸𝑏1 − 𝐸𝑏2


𝑄̇1,2 = ; 𝑄̇1,2 =
1 − 𝜀1 1 1 − 𝜀2 ∑𝑅
+ +
𝜀1 . 𝐴1 𝐴1 . 𝐹1−2 𝜀2 . 𝐴2
Donde: ∑ 𝑅 = sumatoria de resistencias
Envoltura(2)

cuerpo
gris(1)

Si la superficie interior es convexa: 𝐹1−1 = 0 → 𝐹1−2 = 1


Para grandes envolturas: 𝐴2 ≫ 𝐴1
𝐸𝑏1 − 𝐸𝑏2
𝑄̇ =
1 − 𝜀1 1 1 − 𝜀2
𝜀1 . 𝐴1 + 𝐴1 + 𝜀2 . 𝐴2

𝑄̇ = 𝜀1 . 𝐴1 (𝐸𝑏1 − 𝐸𝑏2 )

𝐸𝑏 = 𝜎. 𝑇 4
𝜎 = 5,67𝑥10−8 𝑊 ⁄𝑚2 . 𝐾 4

𝑄̇ = 𝜀1 . 𝐴1 . 𝜎(𝑇14 − 𝑇24 )

PLANOS PARALELOS INFINITOS (𝐴1 = 𝐴2 = 𝐴)

ING. ELÍ GUAYÁN H. BORRADOR Página 6


ING ELI GUAYAN H. BORRADOR RADIACION

1−𝜀1 1−𝜀2
1
𝜀1 𝜀2

Donde: 𝐹1−2 = 𝐹2−1 = 1


𝑄̇(1,2) 𝐸𝑏1 −𝐸𝑏2 𝑊
Calor neto por unidad de superficie: = 1 1 [ ]
𝐴 + −1 𝑚2
𝜀1 𝜀2

RECINTO CERRADO DE TRES SUPERFICIES (CUBIERTA GRIS DE TRES


COMPONENTES):

𝐸𝑏1 − 𝐽1 𝐽2 − 𝐽1 𝐽3 − 𝐽1
+ +
𝑅1 𝑅12 𝑅13
𝐽1 − 𝐽2 𝐸𝑏2 − 𝐽2 𝐽3 − 𝐽2
+ +
𝑅12 𝑅2 𝑅3

𝐽1 − 𝐽3 𝐽2 − 𝐽3 𝐸𝑏3 − 𝐽3
+ +
𝑅13 𝑅23 𝑅3

Generalmente se presentan dos situaciones prácticas:


1°) Una superficie componente (por ejemplo 𝐴3 ) es mucho mayor que los otros dos.
1−𝜀3
→ } → es una resistencia despreciable
𝜀3 .𝐴3

ING. ELÍ GUAYÁN H. BORRADOR Página 7


ING ELI GUAYAN H. BORRADOR RADIACION

2°) Una de las superficies se encuentra térmicamente aislada (𝑞̇ = 0 , emite radiación a la
misma rapidez con la que la absorbe)

𝑄̇1−3 − 𝑄̇3−2 = 0
Nodo 3: en equilibrio radiante
𝐴3 : es una superficie reirradiante

Superficie reirradiante:
Aquella que en condiciones estacionarias, pierde energía de radiación a igual ritmo
que la que gana y Qi = cero, se dice que la superficie vuelve a irradiar todo la energía
de radiación que recibe.

Ejercicio de aplicación

Para una esfera cuyo diámetro externo es de 30cm, situada dentro de una esfera de
60cm de diámetro interior, con 𝜀1 = 1 y 𝜀2 = 0,5 , determine el intercambio de calor 𝑄̇1,2 ,
si 𝑇1 = 527𝑜 𝐶 y 𝑇2 = 277𝑜 𝐶.

Solución

ING. ELÍ GUAYÁN H. BORRADOR Página 8


ING ELI GUAYAN H. BORRADOR RADIACION

3 1 =1 3 2 =0,5
2
R1
R2 R1=0,15m R =0,30m
1
T =800oK T =550 oK

𝑇1 = 527𝑜 𝑐 = 800𝑜 𝐾

𝑇2 = 277𝑜 𝑐 = 550𝑜 𝐾

𝐸𝑏1 − 𝐸𝑏2
𝑄̇1,2 =
1 − 𝜀1 1 1 − 𝜀2
𝜀1 𝐴1 + 𝐴1 𝐹1−2 + 𝜀2 𝐴2
Esquema:

E b1 J1 J2 Eb2

1- 13 1 1- 2
3
3 1 A1 A1 F1-2 3 2 A2

Resistencia R1 Resistencia R 1-2 Resistencia R2

a) Las potencias emisivas son:


𝐸𝑏1 = 𝜎. 𝑇14 = (5,67𝑥10−8 )(8𝑥102 )4 = 23224,32 𝑊/𝑚2
𝐸𝑏2 = 𝜎. 𝑇24 = (5,67𝑥10−8 )(5,5𝑥102 )4 = 5188,4 𝑊/𝑚2

b) Las resistencias del circuito son: (de superficie y espacial)

1−𝜀1 1−1
*) Resistencia (𝑅1 ): 𝑅1 = 𝜀1 𝐴1
= (1)[4𝜋(0,15)2 ] = 0

1−𝜀2 1−0,5
*) Resistencia (𝑅2 ): 𝑅2 = 𝜀2 𝐴2
= (0,5)[4𝜋(0,3)2 ] = 0,88419

1 1
*)Resistencia (𝑅1−2): 𝑅1−2 = 𝐴 𝐹
= [4𝜋(0,3)2 ].1 = 3,53677 ; 𝐹1−2 = 1
1 1−2

Entonces el calor por radiación:

𝐸𝑏1 − 𝐸𝑏2 23224 − 5188,4 18035


𝑄̇1,2 = = =
𝑅1 + 𝑅1−2 + 𝑅2 0 + 3,53677 + 0,88419 4,41867

𝑄̇1,2 = 4081,75 𝑊 𝑅𝑝𝑡𝑎.

ING. ELÍ GUAYÁN H. BORRADOR Página 9


BORRADOR

BORRADOR
BORRADOR

BORRADOR
BORRADOR

FACTOR DE FORMA
Hallar el factor de forma: F1-2

BORRADOR

También podría gustarte