Universidad Nacional de Trujillo

Escuela Profesional de Ingeniería Mecánica

Distribución de calor en una esfera en estado transitorio

Alumna:

Elías Díaz Lira

Curso:

Transferencia de Calor

Docente:

Ing. Jhon Haro

Semestre Académico:

2022-II
Distribución de calor en una esfera en estado transitorio

Hipótesis:

- Propiedades del material constantes

- Sin generación de calor
- Transferencia de calor unidimensional (radial)

Condiciones de contorno asumidas:

- 𝑈(𝑟, 0) = 𝑈0
- 𝑈(𝑅, 𝑡) = 𝑈𝑠
- 𝑈(𝑟, ∞) = 𝑈𝑠

Demostración:

Definimos 𝑇(𝑟, 𝑡) = 𝑈(𝑟, 𝑡) − 𝑈(𝑟, ∞)

Ecuación general de conducción de calor en coordenadas esféricas:

1 𝜕 2 𝜕𝑇 1 𝜕 2𝑇 1 𝜕 𝜕𝑇 𝑔̇ 𝑣 𝜌𝐶𝑝 𝜕𝑇
2
[𝑟 ] + 2 2 2
+ 2 [sin 𝜃 ] + =
𝑟 𝜕𝑟 𝜕𝑟 𝑟 sin 𝜃 𝜕𝜑 𝑟 sin 𝜃 𝜕𝜃 𝜕𝜃 𝑘 𝑘 𝜕𝑡

Por las hipótesis, obtenemos:

1 𝜕 2 𝜕𝑇 𝜌𝐶𝑝 𝜕𝑇
2
[𝑟 ]=
𝑟 𝜕𝑟 𝜕𝑟 𝑘 𝜕𝑡

1 𝜕𝑇 2
𝜕 2𝑇 𝜌𝐶𝑝 𝜕𝑇
2
(2𝑟 + 𝑟 2
)=
𝑟 𝜕𝑟 𝜕𝑟 𝑘 𝜕𝑡

2 𝜕𝑇 𝜕 2 𝑇 𝜕𝑇
+ 2 =𝛼 (1)
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑡

𝑣(𝑟, 𝑡) = 𝑟. 𝑢(𝑟, 𝑡)

𝜕𝑣 𝜕𝑇
=𝑟 (2)
𝜕𝑡 𝜕𝑡

𝜕𝑣 𝜕𝑇
= 𝑇(𝑟, 𝑡) + 𝑟
𝜕𝑟 𝜕𝑟
 𝜕 2𝑣 𝜕𝑇 𝜕 2𝑇
= 2 + 𝑟.
𝜕𝑟 2 𝜕𝑟 𝜕𝑟 2

𝜕 2𝑇 𝜕 2𝑢 𝜕 2𝑣
2. + 𝑟. = (3)
𝜕𝑟 2 𝜕𝑟 2 𝜕𝑟 2

Reemplazando (2) y (3) en (1)

𝜕𝑣 𝜕 2 𝑣
𝛼 = (4)
𝜕𝑡 𝜕𝑟 2

Supongamos que 𝑣 es de la forma: 𝑣 = 𝑀(𝑟)𝑁(𝑡)

𝜕𝑣
→ = 𝑀′ (𝑟)𝑁(𝑡)
𝜕𝑟

𝜕 2𝑣
→ = 𝑀′′ (𝑟)𝑁(𝑡)
𝜕𝑟 2

𝜕𝑣
→ = 𝑀(𝑟)𝑁′(𝑡)
𝜕𝑡

Reemplazando en (4), obtenemos:

𝑀′′ (𝑟)𝑁(𝑡) = 𝛼𝑀(𝑟)𝑁′(𝑡)

Dividiendo entre 𝑀(𝑟)𝑁(𝑡):

𝑀′′(𝑟) = 𝛼𝑁′(𝑡) = −𝜆2 (5)

Tomamos 𝜆2 > 0, debido a que los demás casos no resuelven el problema físico.

De la ecuación 5 obtenemos las siguientes 2 ecuaciones diferenciales:

𝑀′′ (𝑟) + 𝜆2 𝑀(𝑟) = 0 (6)

′ (𝑡)
𝜆2
𝑁 + 𝑁(𝑡) = 0 (7)
𝛼

De la ecuación (7), tenemos que:

𝑑𝑁 𝜆2
= − 𝑑𝑡
𝑁 𝛼
 𝜆2
ln(𝑁) = − 𝑡 + ln 𝐶1
𝛼

𝜆2
𝑁(𝑡) = 𝐶1 𝑒 − 𝛼 𝑡 (8)

La ecuación (6) tiene una solución análoga a la del MAS:

𝑀(𝑟) = 𝐶2 . sin 𝜆𝑟 + 𝐶3 . cos 𝜆𝑟

Usando condiciones de contorno, descartando soluciones triviales:

Tenemos que 𝑇(0, 𝑡) tiene que ser finito, como:

𝑣(𝑟, 𝑡)
𝑇(𝑟, 𝑡) =
𝑟

sin 𝜆𝑟
→ 0 𝑐𝑢𝑎𝑛𝑑𝑜 𝑟 → 0
𝑟

cos 𝜆𝑟
→ ∞ 𝑐𝑢𝑎𝑛𝑑𝑜 𝑟 → 0
𝑟

Por tanto, para que sea convergente, 𝐶3 = 0

𝑀(𝑟) = 𝐶2 . sin 𝜆𝑟

A su vez, tenemos que:

𝑇(𝑅, 𝑡) = 0, descartando la solución trivial, tenemos:

𝑀(𝑅) = 𝐶2 . sin 𝜆𝑅 = 0

𝑛𝜋
→ 𝜆𝑛 =
𝑅

Por el principio de superposición, obtenemos que 𝑣(𝑟, 𝑡):

∞
2
𝑣(𝑟, 𝑡) = 𝑟. 𝑢(𝑟, 𝑡) = ∑ 𝐴𝑛 sin(𝜆𝑛 𝑟) . 𝑒 −𝛼𝜆𝑛𝑡
𝑛=1

Tenemos que en 𝑡 = 0, 𝑇(𝑟, 0) = 𝑇0 − 𝑇𝑠 :

 ∞
𝜋𝑛
𝑟. 𝑢(𝑟, 0) = ∑ 𝐴𝑛 sin ( 𝑟)
𝑅
𝑛=1

𝑚𝜋
Multiplicando ambos miembros por sin 𝑟 e integramos entre 0 y R:
𝑅

𝑅 ∞ 𝑅
𝑚𝜋𝑟 𝑛𝜋𝑟 𝑚𝜋𝑟
∫ 𝑟. 𝑢(𝑟, 0). sin ( ) 𝑑𝑟 = ∑ 𝐴𝑛 ∫ sin ( ) . sin ( ) 𝑑𝑟
0 𝑅 0 𝑅 𝑅
𝑛=1

El resultado de la integral del segundo miembro cuando 𝑚 ≠ 𝑛 y cuando 𝑚 = 𝑛 es:

𝑛𝑟 𝑛
Hacemos el cambio de variable 𝑧 = , 𝑑𝑧 = 𝑅 𝑑𝑟
𝑅

𝑅
𝑚𝜋𝑟 𝑛𝜋𝑟 𝑅 𝜋
∫ sin ( ) . sin ( ) 𝑑𝑟 = ∫ sin(𝑚𝑧) . sin(𝑛𝑧) 𝑑𝑟
0 𝑅 𝑅 𝜋 0

𝜋 𝜋
𝑛2 1 𝑛
(1 − 2 ) ∫ sin(𝑚𝑧) sin(𝑛𝑧) 𝑑𝑧 = − sin(𝑛𝑧) . cos(𝑚𝑧) + 2 cos(𝑛𝑧) . sin(𝑛𝑧)| = 0
𝑚 0 𝑚 𝑚 0

𝜋 𝜋
1 𝜋 1 1 𝜋
∫ sin2 (𝑛𝑧) 𝑑𝑧 = ∫ (1 − cos(2𝑛𝑧))𝑑𝑧 = (𝑧 − sin(2𝑛𝑧))| =
0 2 0 2 2𝑛 0 2

Despejamos 𝐴𝑛 :

2 𝑅 𝑛𝜋𝑟
𝐴𝑛 = ∫ 𝑟. 𝑢(𝑟, 0) sin ( ) 𝑑𝑟
𝑅 0 𝑅

Entonces:

𝜋 (−1)𝑛+1
2𝑅 2𝑅
𝐴𝑛 = (𝑇 − 𝑇 ) ∫ 𝑧. sin(𝑛𝑧) 𝑑𝑧 = (𝑇 − 𝑇 ).
𝜋2 0 𝑠
0 𝜋 0 𝑠
𝑛

Por tanto, la temperatura 𝑈(𝑟, 𝑡) para cualquier punto a una distancia r del centro de la
esfera, en un instante 𝑡, se compone de la suma de la temperatura en el estado estacionario
𝑈𝑠 y de una serie convergente:
∞
2𝑅 (−1)𝑛+1 𝑛𝜋𝑟 𝑛2 𝜋 2
−𝛼. 2 𝑡
𝑈(𝑟, 𝑡) = 𝑇𝑠 + (𝑇 − 𝑇𝑠 ). ∑ . sin ( ).𝑒 𝑅
𝜋 0 𝑛 𝑅
𝑛=1