Equilibrio térmico

Equilibrio térmico

Consideremos dos bloques iguales de longitud L, fabricados con materiales de capacidades térmicas C1=m1·c1 y C2=m2·c2. Donde m1 y m2 son las masas de los bloques y c1 y c2 sus capacidades térmicas másicas en J/(kg·K).

Inicialmente cada bloque está a una temperatura T1 y T2, en recintos adiabáticos separados. En el instante t=0, se elimina la pared de separación y se ponen en contacto. Después de un tiempo, teóricamente infinito, los bloques adquieren una temperatura común de equilibrio Teq, ya que el calor cedido por el boque caliente Q es absorbido por el bloque más frío (se supone que no hay pérdidas de calor).

Q= C 1 ( T eq T 1 )Q= C 2 ( T eq T 2 ) T eq = C 1 T 1 + C 2 T 2 C 1 + C 2

Si los bloques son iguales, hechos del mismo material, la temperatura de equilibrio es Teq=(T1+T2)/2

Evolución hacia el estado de equilibrio

La ecuación diferencial que describe la conducción térmica es

T t =α 2 T x 2 α= K ρc

donde K, es la conductividad térmica, ρ la densidad, y c es la capacidad térmica másica o calor específico del material.

Las condiciones de contorno corresponden a extremos aislados

T(x,t) x | x=0 =0 T(x,t) x | x=2L =0

En la figura se representa el estado inicial T(x,0): las temperaturas iniciales T1 y T2 de los dos bloques de longitud L y el estado estacionario: la temperatura final de equilibrio T(x,∞)=(T1+T2)/2 .

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

Definimos la función u(x,t)=T(x,t)-T(x,∞), en términos de esta nueva función, la ecuación de la conducción del calor, las condiciones de contorno en x=0 y en x=2L y la condición inicial en el instante t=0, se escriben

1 α u t = 2 u x 2 0<x<2Lt>0 ·contorno u(x,t) x | x=0 =0 u(x,t) x | x=2L =0 ·inicialu(x,0)=T(x,0) T eq

Buscamos una solución de la forma u(x, t)=F(xG(t), variables separadas

1 α 1 G(t) dG(t) dt = 1 F(x) d 2 F(x) d x 2 = ω 2

Integramos la primera ecuación diferencial

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

Integramos la segunda ecuación diferencial

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

Es una ecuación diferencial similar a la de un Movimiento Armónico Simple, cuya solución es F(x)=Acos(ωx)+Bsin(ωx)

Condiciones de contorno

La solución completa u(x,t) es la superposición de los productos de las funciones Fn(xGn(x) para cada ωn=nπ/(2L)

u(x,t)= n=1 A n cos( ω n x )exp( α ω n 2 t )

Condición inicial

Solamente, queda por determinar los coeficientes An, identificando la solución para t=0 con la condición inicial u(x,0)

u(x,0) = n=1 A n cos( nπ x 2L ) u(x,0) =T(x,0) T eq

Multiplicamos ambos miembros por cos(mπx/(2L)) e integramos entre 0 y 2L

0 2L u(x,0) cos ( mπ 2L x ) dx= n=1 A n 0 2L cos( nπ 2L x )cos( mπ 2L x ) dx

Despejamos An

A n = 2L nπ ( T 1 T 2 )sin( nπ 2 ){ 0npar 2L nπ ( T 1 T 2 ) ( 1 ) n1 2 nimpar

Solución completa

u(x,t)= n=1 A n sin( nπ x 2L )exp( α n 2 π 2 4 L 2 t ) T(x,t)=u(x,t)+T(x,)

Así, la temperatura en cualquier punto de la varilla x, en un instante t, se compone de la suma de un término Teq, la temperatura de equilibrio y de una serie rápidamente convergente que describe el estado transitorio.

T(x,t)= T 1 + T 2 2 + 2L π ( T 1 T 2 ) n=1,3,5... (1) n1 2 n exp( α n 2 π 2 4 L 2 t ) cos( nπ 2L x )

Creamos la función temperatura, que admite como parámetros: las temperaturas iniciales T1 y T2 de los dos bloques iguales, la longitud L de cada bloque, el parámetro 1/α, y el tiempo t en el que queremos calcular la distribución de temperaturas a lo largo de los bloques de 0 a 2L. Devuelve un vector x de puntos a lo largo de los bloques y sus correspondientes temperturas en un vector T.

function [x,T]=temperatura_6(T1,T2,L,a2,t)
   x=linspace(0,2*L,100);
   T=ones(1,length(x))*(T1+T2)/2; 
   cte=pi^2/(a2*4*L^2);
   n=1;
   v=1;
   while(v>0.01)
        v=exp(-cte*n^2*t);
        T=T+2*L*(T1-T2)*v*cos(pi*n*x/(2*L))*(-1)^((n-1)/2)/(n*pi);
        n=n+2;
   end
end

Una serie infinita tenemos que interrumpirla en algún término de acuerdo con algún criterio. El que se ha empleado en este programa es el siguiente

exp( α n 2 π 2 4 L 2 t )<0.01

sea menor que un valor prefijado, por ejemplo 0.01.

Creamos un script que establece

Defina el vector de los instantes t, en segundos, en el que queremos representar la distribución de temperaturas a lo largo de los bloques.

T1=100; %temperatura del bloque izquierdo
T2=30; %temperatura del bloque derecho
alfa=11352; %Coeficiente, 1/alfa del aluminio
L=0.5; % longitud de la varilla
  
hold on
for t=[10 100 200 400, 800]
    [x,T]=temperatura_6(T1,T2,L,alfa,t);
    plot(x,T,'displayName',num2str(t));
end
title('Evolución de la temperatura de los bloques')
xlabel('x')
ylabel('T')
legend('-DynamicLegend','location','northeast')
grid on
hold off 

Establecimiento del equilibrio térmico

El equilibrio térmico se alcanza después de un tiempo teóricamente infinito. En la práctica, depende de la sensibilidad de los termómetros empleados para medir la temperatura. Diremos que se ha alcanzado el equilibrio térmico en el instante en el que la diferencia entre la temperatura del extremo izquierdo T(0,t) y la temperatura del extremo derecho T(2L,t) sea menor que un grado

Para calcular este instante, supondremos que después de un tiempo grande t la distribución de temperaturas viene dada esencialmente por el primer término del desarrollo en serie

T(x,t) T 1 + T 2 2 + 2L π ( T 1 T 2 )exp( α π 2 4 L 2 t )cos( π 2L x )

La diferencia de temperatura en los extremos es aproximadamente

T(0,t) T 1 + T 2 2 + 2L π ( T 1 T 2 )exp( α π 2 4 L 2 t ) T(2L,t) T 1 + T 2 2 2L π ( T 1 T 2 )exp( α π 2 4 L 2 t ) T(0,t)T(2L,t)= 4L π ( T 1 T 2 )exp( α π 2 4 L 2 t )

Añadimos al script la porción de código que calcula el tiempo t en segundos, a partir del cual consideramos que se ha establecido el equilibrio térmico

dT=1; %un grado, la diferencia de temperatura entre los extremos
tFin=-4*L^2*log(dT*pi/(4*L*(T1-T2)))*alfa/pi^2
tFin =   4.3672e+03

Referencias

F M S Lima. Approximate formulae for the thermalization of two blocks in an insulated system. Eur. J. Phys. 35 (2014) 055017