Ecuación diferencial de conducción de calor en coordenadas

esféricas

I. Fundamentación:

Uno de los principales objetivos en un análisis de conducción de calor es determinar el campo

de temperatura en un medio que resulta de las condiciones impuestas sobre sus fronteras. Es
decir, deseamos conocer la distribución de temperaturas, que representa como varia la
temperatura con la posición en el medio. Una vez que se conoce esta distribución, el flujo de
calor por conducción en cualquier punto en el medio o en la superficie se calcula a partir de la
ley de Fourier.

La distribución de temperaturas para un sólido sirve para comprobar la integridad estructural

mediante la determinación de los esfuerzos térmicos, sus expansiones y deflexiones; así como
también, es útil para optimizar el espesor de un material aislante o para determinar la
compatibilidad de recubrimiento o adhesivos especiales que se usen con el material.

La ecuación diferencial de conducción de calor nos va a servir para determinar el campo de

temperaturas en un cuerpo o sistema. El fenómeno térmico se estudiara en un intervalo de
tiempo pequeño y en un volumen elemental del espacio, seguiremos la metodología de
aplicación de la conservación de la energía.

II. Deducción de la Ec. diferencial de calor para coordenadas esféricas:

2.1 Volumen de control:

Para la deducción definimos antes que nada un volumen de control infinitesimal pequeño
(diferencial), que para las coordenadas en que la deduciremos geométricamente es una cuña
esférica, como se muestra en la figura 1.
 Figura 1 Volumen diferencial de control para el análisis de conducción en coordenadas esféricas (r, ϕ, θ)

2.2 Hipótesis de trabajo:

- Medio homogéneo e isotrópico.

- Las variaciones de volumen debido al cambio de temperatura son despreciables.

- Propiedades físicas no se alteran o cambian, suponen constantes.

- Existe una fuente de calor uniformemente distribuida.

2.3 Principio de conservación de la energía:

Aplicaremos el principio de la conservación de la energía para un instante de tiempo “t”, el cual

enuncia lo siguiente: la velocidad de la energía calorífica neta impartida por la transferencia de
calor a través de todas sus superficies más la velocidad de la energía calorífica proporcionada
por las fuentes internas de calor es igual a la velocidad del incremento de la energía interna
considerada.

dU
E˙ent − E˙sal + Ė g=∆ E˙alm …….…balance entre las velocidades de energíadq 1+ dq2=
dt
…………………(I)

Dónde:

dq 1=velocidad neta de tranferencia de c alor por conducción [ W ] . ( E˙ent− E˙sal )

dq 2=velocidad de transferencia de calor liberada por las fuentes internas de calor [ W ] .( Ė g)

dU
=velocidad de incremento de la energía interna de la sustancia contenida en el volumen de control [ W ] . (∆
dt

2.4 Desarrollo:

2.4.1 Cálculo de áreas del elemento diferencial de volumen de control.- sea la siguiente
figura:
 Figura 2 Elemento diferencial de volumen en coordenadas esféricas

2.4.1.1 Según la dirección radial (r):

⇒ ⇒
❑ d Ar =r . dθ . r sin θ . dϕ ❑ d Ar +dr =( r +dr ) dθ .(r + dr). sinθ . dϕ
2.4.1.2 Según la dirección polar (ϕ):
⇒ ⇒
❑ d Ar =r 2 . sinθ . dθ . dϕ ❑ d Ar +dr =(r ¿¿ 2+2 r . dr). sin θ .dθ . dϕ ¿

⇒
dr
❑ d Aϕ =d A ϕ+dϕ= ¿)
2
⇒
❑ d Aϕ =d A ϕ+dϕ=r . dr . dθ

2.4.1.3 Según la dirección azimutal (θ):

⇒
dr ⇒
dr
❑ d Aθ = ¿) ❑ d Aθ +dθ = ¿)
2 2
⇒ ⇒
❑ d Aθ =r .sin θ dr . dϕ ❑ d Aθ +dθ=r . sin θ. dr . dϕ+r . cos θ . dr .dϕ . dθ
2.4.2 Cálculo de la velocidad neta de transferencia de calor introducida al elemento por
conducción:

Tenemos que:

dq 1=dq 1r + dq1 ϕ +dq1 θ …………………………(II)

Con:

dq 1 r=dq r−dqr +dr ………………………………(1)

dq 1 ϕ=dqϕ −dqϕ+dϕ ………………………….....(2)

dq 1θ =dqθ −dq θ+dθ………………………..…….(3)

Ahora, sabemos que el flujo de calor es:

dq ⇒
q '' = ❑ dq=q '' .dA
dA

Dónde:

q '' =flujo de calor [W /m2 ]

dq=velocidad de transferencia de calor por conducción [W ]

dA=área de conducción[m2 ]

- Para dq r:
dq r=q'r' . d A r=q'r' . r 2 . sin θ . dθ . dϕ
dq r=q'r' . r 2 . sin θ . dϕ dθ

- Para dq r+ dr: aplicando Serie de Taylor

∂ q'r'
(
dq r+ dr =q ''r +dr . d A r +dr = q ''r +
∂r )
. dr ¿

⇒
'' 2 '' ∂ q'r' 2
❑ dqr +dr=qr . r . sin θ . dθ . dϕ+¿ q r .2 r . sinθ . dr . dθ . dϕ+ .r . sin θ .dr . dθ . dϕ ¿
∂r

Reemplazando estos resultados en la ecuación (1):

dq 1 r=dq r−dqr +dr

'' 2 '' 2 '' ∂q ''r 2

dq 1 r=q . r . sin θ . dϕ dθ−q . r . sin θ . dθ . dϕ−¿ q .2r . sin θ . dr . dθ . dϕ−
r r r . r . sinθ . dr . dθ . dϕ ¿
∂r
⇒ ∂ q'r' 2
( ''
❑ dq1 r=− q .2 r . sin θ+
r
∂r )
. r . sin θ . dr . dθ . dϕ
- Para dq ϕ:
dq ϕ =q'ϕ' . d A ϕ =q ''ϕ . r . dr . dθ
⇒
❑ dqϕ =q ''ϕ . r . dr . dθ

- Para dq ϕ+dϕ: aplicando Serie de Taylor

∂ q'ϕ'
dq ϕ+dϕ =q'ϕ+dϕ
'
. d A ϕ+ dϕ= q 'ϕ' + ( ∂ϕ )
. dϕ ( r . dr .dθ )

⇒
'' ∂ q'ϕ'
❑ dqϕ+dϕ =q ϕ .r . dr . dθ+ . r . dr . dθ . dϕ
∂ϕ

Reemplazando estos resultados en la ecuación (2):

dq 1 ϕ=dqϕ −dqϕ+dϕ

'' '' ∂ q'ϕ'

dq 1 ϕ=q . r . dr . dθ−q .r . dr . dθ−
ϕ ϕ . r .dr . dθ . dϕ
∂ϕ
⇒ −∂ q'ϕ'
❑ dq1 ϕ = .r . dr . dθ . dϕ
∂ϕ

- Para dq θ:
dq θ=q 'θ' . d Aθ =q'θ' . r . sinθ dr . dϕ
⇒
❑ dqθ=q'θ' . r .sin θ dr . dϕ

- Para dq θ+dθ : aplicando Serie de Taylor

∂ q'θ'
(
dq θ+dθ =q'θ'+dθ . d A θ+dθ = q ''θ +
∂θ )
. dθ ( r . sin θ . dr . dϕ+r . cos θ . dr . dϕ .dθ )

⇒
'' '' ∂ q'θ'
❑ dqθ+ dθ =qθ . r . sinθ . dr . dϕ+q θ . r . cos θ . dr . dϕ . dθ+ .r . sin θ .dr . dϕ . dθ
∂θ

Reemplazando estos resultados en la ecuación (3):

dq 1θ =dqθ −dq θ+dθ

'' '' '' ∂ q'θ'

dq 1θ =q . r . sin θ dr . dϕ−q .r . sin θ . dr . dϕ−q . r . cos θ . dr . dϕ . dθ−
θ θ θ .r . sin θ .dr . dϕ . dθ
∂θ
⇒
❑ dq1 θ=−¿
 Finalmente, reemplazando estos últimos datos en la ecuación (II):

dq 1=dq 1r + dq1 ϕ +dq1 θ

∂ q'r' 2 ∂ q'ϕ' ∂q 'θ'

[ ''
r
''
dq 1= −( q .2 r .sin θ+ q . r . cos θ ) −
θ
∂r( . r . sin θ+
∂ϕ
.r+
∂θ )]
. r . sin θ dr . dϕ . dθ…...(III)

Pero, sabemos que por ley de Fourier:

∂T
q ''r =−k .

q '' =−k .
∂T
∂η
⇒

{
❑ q 'ϕ' =−k .
''
∂r
∂T
r . sinθ . ∂ ϕ
q θ =−k .
∂T
r . ∂θ

 Nota: en el denominador del gradiente de temperatura se toma la dimensión de la dirección

de coordenadas específicas de acuerdo al volumen de control, ver direcciones en figura 1.

Ahora, reemplazando estas ecuaciones en la ecuación (III):

[(
dq 1= − −k .
∂T
∂r
.2 r .sin θ−k .
∂T
r .∂θ )(
. r .cos θ − −r 2 . sin θ .
∂
∂r
k.(∂T
∂r
−r .
∂
)
∂ϕ
k.
∂T
(
r . sin θ .∂ ϕ )
−r . sin θ .
∂
∂θ
k.
∂T
r .∂ (
dq 1=
[(
∂T
∂r
.2 r .sin θ+
∂T
r . ∂θ )(
. r . cos θ + r 2 . sin θ .
∂ ∂T
( )
∂ r ∂r
+r.
∂
( ∂T
∂ ϕ r . sin θ . ∂ ϕ )
+r . sinθ .
∂ ∂T
(
∂θ r . ∂θ ))] k .dr . dϕ . dθ
⇒
∂ 2 ∂T 1 ∂ ∂T ∂ ∂T
[
❑ dq1= sinθ .
∂r (
r . + ) . + ( ) (
∂ r sin θ ∂ ϕ ∂ ϕ ∂θ
sinθ .
∂θ )]
k . dr . dϕ . dθ ………………..(a)

2.4.3 Cálculo de la velocidad de transferencia de calor “liberada” por las fuentes

internas de calor:

dq 2=q̇ f . dV

Dónde:

q̇ f =rapidez a la que se genera energia por unidad de volumendel medio [W /m 3 ]

dV =r 2 .sin θ . dr . dϕ. dθ

Por lo tanto,
 ⇒
❑ dq2=q̇ f . r 2 .sin θ . dr . dϕ. dθ ……………(b)

2.4.4 Cálculo de velocidad de incremento de la energía interna de la sustancia contenida

dentro del volumen de control:

Tenemos que:

dU dT
=dm . c p .
dt dt

dU dT
=ρ . dV . c p .
dt dt
⇒
dU dT
❑ = ρ. r 2 .sin θ . dr . dϕ .dθ . c p . ………………(c)
dt dt

Ahora, reemplazando ecuaciones (a), (b), (c) en ecuación (I):

dU
dq 1+ dq2=
dt

∂ 2 ∂T 1 ∂ ∂T ∂ ∂T
[ sin θ .
∂r(r . + ) . + ( ) (
∂r sin θ ∂ ϕ ∂ ϕ ∂ θ
sin θ .
∂θ )] k . dr . dϕ . dθ+q̇ . r . sin θ . dr . dϕ . dθ=ρ .r . sin θ .dr . dϕ . dθ . c
f
2 2
p

1 ∂ 2 ∂T 1 ∂ ∂T 1 ∂ ∂T dT
k.
[ .
r ∂r
2 (
r . + 2 ) . + 2( ) .
∂ r r . sin θ ∂ ϕ ∂ ϕ r . sin θ ∂θ
2
sin θ .
∂θ( )]
+ q̇ f =ρ . c p .
dt

k 1 ∂ 2 ∂T 1 ∂ ∂T 1 ∂ ∂T q̇ dT
. 2.
ρ .c p r ∂r[ ( )
r . + 2 . + 2 ( )
.
∂ r r . sin θ ∂ ϕ ∂ ϕ r . sin θ ∂θ
2
sinθ .
∂θ
+ f = (
ρ . c p dt )]
1 ∂ 2 ∂T 1 ∂ ∂T 1 ∂ ∂T q̇ dT
∴a.
[ ( )
r 2
.
∂ r
r .
∂r
+ 2 2
r .sin θ
.
∂ ϕ ∂ ϕ( )
+ 2
r .sin θ
.
∂ θ (
sin θ .
∂ θ
+ f =
ρ. )]
c p dt

Dónde:

k
a= =difusividad térmica[m2 / s]
ρ. c p

III. Conclusiones:

Como observamos el resultado obtenido es la forma general, en coordenadas esféricas de

la ecuación diferencial del calor, la misma que será de gran ayuda para obtener la
distribución de temperaturas T(r, ϕ, θ) como función del tiempo.
IV. Bibliografía:

1. INCROPERA, F.P.; DE WITT, D.P “Fundamentos de transferencia de Calor”, 4ta.

Ed. Edit. Mc GrawHill, Mexico 1999.
2. 2. YUNUS A. CENGEL “Transferencia de Calor”, 2da. Ed. Edit Mc GrawHill, Mexico
2004.
3. http://www.das.uchile.cl/~mhamuy/courses/AS42B/tema2.html
4. http://es.scribd.com/doc/22495144/La-ley-de-Fourier-y-la-ecuacion-de-calor
5. http://es.wikipedia.org/wiki/Coordenadas_esf%C3%A9ricas