Universidad Nacional Mayor de San Marcos
Facultad de Ciencias Físicas
E.A.P. Física
“Trabajo Parcial de Análisis Numérico”
Alumno:
Quispe Padilla Branco Alexande
Código de estudiante:
16130020
Profesor:
Raul Felix Carita Montero
Curso:
Análisis Numérico
Lima-Perú
2021
Conducción del calor en una esfera homogénea
En este trabajo desarrollaremos 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
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, la condición inicial en el instante t=0,
y las condiciones de contorno en r=R se escriben
Si hacemos la sustitución v(r,t)=r·u(r,t) nos queda la educación diferencial
en derivadas parciales
Buscamos soluciones de la forma v(r, t)=F(r)·G(t) (variables separadas)
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)
Condición inicial
La condición inicial es u(r,0)=T0-Ts, que nos permite determinar los
coeficientes An
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
El resultado de la integral del segundo miembro cuando m≠n y
cuando m=n es
Teniendo en cuenta que u(r,0)=T0-Ts y efectuando el cambio de
variable z=πr/R
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.
Programación en Scilab:
Temperatura en el centro de la esfera
Temperatura en el centro de la esfera cuando r→0
en función del tiempo t
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')
xgrid
Temperatura en el centro de un huevo cocido
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 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
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',string(R));
end
xgrid
plot([0,1400],[85,85],'lineStyle','--','color','k')
title('Temperatura en el centro del huevo')
legend('-DynamicLegend','location','southeast')
xlabel('t')
ylabel('T')
En la serie infinita, tomamos el primer término n=1 y despreciamos los
restantes
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
