Inicio Fenómenos de transporte Conducción del calor

Conducción del calor en una esfera

homogénea (I)
En esta página, se estudia la distribución de temperaturas de una esfera cuya temperatura
inicial T0 es uniforme y se sumerge en un baño térmico a la temperatura Ts.

Consideremos una esfera de radio R en la cual la distribución inicial de temperaturas y las

condiciones de contorno tienen simetría esférica. Las superficies isotérmicas son superficies
esféricas concéntricas y la temperatura es una función únicamente de la distancia radial r y del
tiempo t.

La esfera se calienta hasta una temperatura uniforme T0 (distribución inicial de temperaturas

T(r, 0)=T0. En el instante t=0, se sumerge en un recipiente grande de agua a temperatura Ts
que es continuamente agitada. La condición de contorno es por tanto, T(R, t)=Ts.

En el estado estacionario, después de un tiempo t→∞, la temperatura final de la esfera será

T(r,∞)=Ts, la temperatura del baño térmico.

La ecuación de la conducción del calor apropiada para resolver este problema es

1 ∂T 1 ∂ 2 ∂T
= (r )
α ∂t r2 ∂r ∂r

2
1 ∂T 2 ∂T ∂ T
= +
α ∂t r ∂r ∂r2

Solución de la ecuación de la conducción del calor

Definimos la función u(r,t)=T(r,t)-T(r,∞), en términos de esta nueva función, la ecuación de la
conducción del calor, las condición inicial en el instante t=0, y las condiciones de contorno en
r=R se escriben
 2
1 ∂u ∂ u 2 ∂u
= + 0 < r < R t > 0
α ∂t ∂r2 r ∂r

⋅ contorno ⇒ u(R, t) = 0

⋅ inicial ⇒ u(x,0) = T0 − Ts

Si hacemos la sustitución v(r,t)=r·u(r,t) nos queda la educación diferencial en derivadas

parciales
2
∂v ∂ v
= α
2
∂t ∂ r

Buscamos soluciones de la forma v(r, t)=F(r)·G(t) (variables separadas)

2
∂v(r,t) ∂ v(r,t)
= α
2
∂t ∂ r

2
1 1 dG(t) 1 d F (r)
2
= = −ω
α G(t) dt F (r) d2r

Integramos la primara ecuación diferencial

dG(t)
2
+ αω G(t) = 0
dt

2
G(t) = G(0) ⋅ exp(−αω t)

Integramos la segunda ecuación diferencial

2
d F (r)
2
+ ω F (r) = 0
2
d r

Es una ecuación diferencial similar a la de un Movimiento Armónico Simple. La solución es

F(r)=A·sin(ωr)+B·cos(ωr)

En primer lugar, F(0)=0 ya que u(r,t) tiene que ser finito cuando r=0, recuérdese que
u(r,t)=v(r,t)/r y v(r, t)=F(r)·G(t). Esto implica que B=0

Condiciones de contorno
u(R,t)=0, por lo que F(R)=A·sin(ωR)=0, es decir, ωn=nπ/R

La solución v(r,t) es la superposición de los productos Fn(r)·Gn(t). La solución de la ecuación

de la conducción del calor r·u(r,t)=v(r,t)
∞

2
r ⋅ u(r, t) = ∑ An sin(ωn r) ⋅ exp(−αωn t)

n=1

Condición inicial
La condición inicial es u(r,0)=T0-Ts, que nos permite determinar los coeficientes An

nπ
r ⋅ u(r,0) = ∑ An sin ( r)
n
R

Multiplicamos ambos miembros por sin(mπr/R) e integramos entre 0 y R, para ello, hacemos
el cambio de variable z=πr/R, dz=πdr/R

R R
∞
r r r
∫ r ⋅ u(r,0) sin (mπ ) dr = ∑ An ∫ sin (nπ ) sin (mπ ) dr
R R R
n=1
0 0

El resultado de la integral del segundo miembro cuando m≠n y cuando m=n es

 R π

r r R
∫ sin (mπ ) sin (nπ ) dr = ∫ sin (mz) sin (nz) dz
π
R R
0 0

2 π
n 1 n
(1 − )∫ sin (mz) sin (nz) dz = − sin (nz) cos (mz) + cos (nz) sin (nz)∣
∣ = 0
m 2 m m 2
0

π π
π
2 1 1 1 ∣ π
∫ sin (nz) dz = ∫ (1 − cos (2nz)) dz = (z − sin (2nz)) =
2 2 2n ∣ 2
0

0 0

Despejamos An

2 r
An = ∫ r ⋅ u(r,0) sin (nπ ) dr
R R
0

Teniendo en cuenta que u(r,0)=T0-Ts y efectuando el cambio de variable z=πr/R

π
n+1
2R 2R (−1)
An = (T0 − Ts ) ∫ z ⋅ sin (nz)dz = (T0 − Ts )
2
π π n
0

Solución completa
La temperatura T(r,t) en cualquier punto a una distancia r del centro de la esfera, en un
instante t, se compone de la suma de la temperatura en el estado estacionario Ts y de una
serie rápidamente convergente que describe el estado transitorio.

∞ n+1
2 2
(−1) r n π
2R
r ⋅ u(r, t) = (T0 − Ts ) ∑ sin (nπ ) exp (−α t)
π 2
n R R
n=1

T (r, t) = u(r, t) + T (r, ∞) =

∞ n+1
2 2
(−1) r n π
2R
T (r, t) = Ts + (T0 − Ts ) ∑ sin (nπ ) exp (−α t)
πr 2
n R R
n=1

Ejemplo
Creamos la función temperatura_1, que admite como parámetros: las temperatura del baño
Ts, la temperatura inicial, T0, el radio de la esfera R, el parámetro 1/α y el tiempo t en el que
queremos calcular la distribución de temperaturas en la esfera. Devuelve un vector x de
puntos a lo largo de la dirección radial y sus correspondientes temperturas en un vector T.
function [x,T]=temperatura_1(Ts,T0,R,a2,t)
x=linspace(eps,R,100);
if(t==0)
T=T0*ones(1,length(x));
return;
end
T=Ts*ones(1,length(x));
n=0;
v=1;
cte=pi^2/(a2*R^2);
k=2*R*(T0-Ts)./(pi*x);
while(v>0.01)
n=n+1;
v=exp(-cte*n^2*t);
T=T+v*(-1)^(n+1)*k.*sin(pi*n*x/R)/n;
end
end

Una serie infinita tenemos que interrumpirla en algún término de acuerdo con algún criterio.
Un posible criterio sería que la exponencial sea menor que un valor prefijado, por ejemplo
0.01. Nos surge otro problema en el origen x=0, se produce una indeterminación 0/0.
Calculamos la temperatura en el centro de la esfera mediante la expresión simplificada, o bien,
evitar el punto x=0, tomando un punto muy próximo.

Creamos un script que establece

La temperatura del baño térmico, Ts.

La temperatura inicial de la esfera, T0.
El radio de la esfera, R
El valor del coeficiente, 1/α que depende del material del que está hecho las esfera

Define el vector de los instantes t, en segundos, en el que queremos representar la

distribución de temperaturas a lo largo de la dirección radial.

Llama a la función temperatura_1 y representa mediante el comando plot cada una de las
graficas T(x) en la misma ventana.

Ts=100; %Temperatura, baño

T0=0; %Temperatura inicial
R=0.5; %Radio esfera
alfa=11352; %Coeficiente, 1/alfa

hold on
axis([0 R -5 100]);
for t=[2 100 200 300 400]
[x,T]=temperatura_1(Ts,T0,R,alfa,t);
plot(x,T,'displayName',num2str(t));
end
title('Evolución de la temperatura de una esfera')
xlabel('r')
ylabel('T')
legend('-DynamicLegend','location','northwest')
grid on
hold off

Temperatura en el centro de la esfera

Temperatura en el centro de la esfera cuando r→0

π
sin(nr )
R π
lim =
r→0 nr R

en función del tiempo t

 ∞ 2 2
n+1
n π
T (0,t) = Ts + 2(T0 − Ts ) ∑ (−1) exp(−α t)
2
R
n=1

Representamos la evolución de la temperatura en el centro de una esfera de aluminio de radio

R=0.5 m, cuya temperatura inicial es T0=0°C y se sumerge en un baño de agua hirviendo
Ts=100 °C. Como observamos en la gráfica, se tarda un tiempo, alrededor de 100 s, para que
la temperatura del centro de la esfera se vaya incrementando apreciablemente, por encima de
T0=0°C

Ts=100; %Temperatura, baño

T0=0; %Temperatura inicial
R=0.5; %Radio esfera
alfa=11352; %Aluminio, coeficiente, 1/alfa

t=linspace(10,1200,200);
T=Ts*ones(1,length(t));
cte=pi^2/(alfa*R^2);
for i=1:length(t)
v=1;
n=0;
while(v>0.01)
n=n+1;
v=exp(-cte*n^2*t(i));
T(i)=T(i)+v*2*(T0-Ts)*(-1)^(n+1);
end
end

plot(t,T);
title('Temperatura en el centro de una esfera')
xlabel('t')
ylabel('T')
grid on

Huevo duro
Una aplicación práctica es la estimación del tiempo que tarda un huevo sumergido en agua
hirviendo (100 °C) en hacerse duro, tal como se muestra en la fotografía.

En la figura, se muestra la evolución de la temperatura en en el centro de un huevo de radio R

supuesto esférico. La temperatura inicial es T0=25°C y se sumerge en un baño de agua
hirviendo Ts=100 °C. El parámetro α=1.53·10-7 m2/s para la clara y yema del huevo

Ts=100; %Temperatura, baño

T0=25; %Temperatura inicial

hold on
for R=(1:4)/100 %radios del huevo en cm
t=linspace(10,1400,200);
T=Ts*ones(1,length(t));
cte=1.53e-7*pi^2/R^2;
for i=1:length(t)
v=1;
n=0;
while(v>0.01)
n=n+1;
v=exp(-cte*n^2*t(i));
T(i)=T(i)+v*2*(T0-Ts)*(-1)^(n+1);
end
end
plot(t,T,'displayName',num2str(R));
end
line([0,1400],[85,85],'lineStyle','--','color','k')
hold off
title('Temperatura en el centro del huevo')
legend('-DynamicLegend','location','southeast')
xlabel('t')
ylabel('T')
grid on
En la serie infinita, tomamos el primer término n=1 y despreciamos los restantes
2
π
T (0,t) ≈ Ts + 2(T0 − Ts ) exp (−α t)
R2

2 Ts −T (0)
R
tF = − log ( )
π2α 2(Ts −T0 )

Se considera que un huevo está bien cocido cuando su centro alcanza la temperatura
T(0,t)=85°C. Sabiendo que la temperatura inicial del huevo es T0=25°C y la del baño
Ts=100°C, despejamos el tiempo tF=610 s, para un huevo de radio R=2 cm, que es un valor
similar al obtenido en la representación gráfica, donde se han sumado muchos términos de la
serie

Referencias
Unsworth J., Duarte F. J. Heat diffusion in a solid sphere and Fourier theory: An elementary
practical example. Am. J. Phys. 47 (11) November 1979, pp. 981-983

D Buay, S K Foong, D Kiang, L Kuppan, V H Liew. How long does it take to boil an egg?.
Revisited. Eur. J. Phys. 27 (2006) pp. 119-131

Energías Renovables Curso Interactivo de Física en Internet

Angel Franco García, Copyright © 2016