Aproximación al equilibrio de dos gases contenidos en un recinto adiabático y separados por un émbolo (II).

En esta página, se estudia la evolución hacia el estado de equilibrio de un sistema aislado formado por dos gases ideales separados por una pared adiabática.

Un cilindro de sección S y longitud L está cerrado por ambos extremos, contiene un émbolo que divide el volumen en dos partes A y B, ambos contienen la misma cantidad de un gas ideal, un mol. El cilindro y el émbolo están adiabáticamente aislados. El émbolo está sujeto de modo que los volúmenes iniciales de cada parte son, respectivamente V₁₀=S·x₀, V₂₀=S(L-x₀). Las temperaturas iniciales del gas en cada una de las dos partes son T₁₀ y T₂₀.

En el instante t=0, se libera el pistón y suponemos que el émbolo se mueve sin rozamiento. Vamos a determinar el estado final de equilibrio.

El sistema alcanzará el equilibrio cuando las presiones de los gases separados por el émbolo sea la misma p_1f=p_2f=p_f

1.-Aplicamos la ecuación de los gases ideales a cada una de las partes

p_f·V_1f=nRT_1fp_f·V_2f=nRT_2f

2.-Como el sistema es aislado, la energía total permanece constante

U_1f+U_2f=U₁₀+U₂₀

La energía interna de un gas ideal solamente depende de la temperatura

T_1f+T_2f=T₁₀+T₂₀

3.-Los volúmenes de los gases cambian, pero el volumen total es constante e igual al inicial que ocupaban los gases.

V₁₀+V₂₀=V_1f+V_2f

Despejamos la temperaturas finales y la presión final

$\begin{array}{l} T_{1 f} = \frac{(T_{10} + T_{20}) V_{1 f}}{V_{1 f} + V_{2 f}} T_{2 f} = \frac{(T_{10} + T_{20}) V_{2 f}}{V_{1 f} + V_{2 f}} \\ p_{f} = n R \frac{(T_{10} + T_{20})}{V_{1 f} + V_{2 f}} = n R \frac{(T_{10} + T_{20})}{V_{10} + V_{20}} \end{array}$

Ahora bien, de estas ecuaciones no podemos despejar los volúmenes finales ni las temperaturas finales. Como veremos más adelante, la formulación de un modelo de gas ideal nos va a permitir describir la evolución desde el estado inicial al final de equilibrio.

Oscilaciones del émbolo

Supongamos que el émbolo se mueve sin rozamiento. El principio de conservación de la energía para el sistema aislado formado por los dos gases y el émbolo se escribe para un mol de gas.

$\frac{1}{2} M {(\frac{d x}{d t})}^{2} + c_{v} T_{1} + c_{v} T_{2} = c_{v} T_{10} + c_{v} T_{20}$

El primer término es la energía cinética del émbolo de masa M. El segundo, la energía interna de un mol de gas contenido en la parte izquierda y el tercero, la energía interna de un mol de gas contenido en la parte derecha.

Los dos gases experimentan una transformación adiabática

$\begin{array}{l} p V^{γ} = cte T V^{γ - 1} = cte \\ T_{1} x^{γ - 1} = T_{10} x_{0}^{γ - 1} T_{2} {(L - x)}^{γ - 1} = T_{20} {(L - x_{0})}^{γ - 1} \end{array}$

Recuérdese que c_p=c_v+R y γ=c_p/c_v

El principio de conservación de la energía se escribe.

$\frac{1}{2} M {(\frac{d x}{d t})}^{2} + \frac{R}{γ - 1} {T_{10} {(\frac{x_{0}}{x})}^{γ - 1} + T_{20} {(\frac{L - x_{0}}{L - x})}^{γ - 1}} = \frac{R}{γ - 1} (T_{10} + T_{20})$

El segundo término es la energía potencial E_p(x). que se representa en la figura. Para una energía total E, el émbolo oscila entre las dos posiciones x₀ y x₁ señaladas en la figura, con un determinado periodo

El émbolo parte de la posición x₀, en el instante t=0, y llega a la posición de máximo desplazamiento x₁ cuando su velocidad es dx/dt=0. Para calcular esta posición, es preciso resolver por procedimientos numéricos la ecuación trascendente

$T_{10} {(\frac{x_{0}}{x})}^{γ - 1} + T_{20} {(\frac{L - x_{0}}{L - x})}^{γ - 1} = T_{10} + T_{20}$

Ecuación del movimiento

Derivamos la ecuación de la conservación de la energía respecto con respecto del tiempo

$M \frac{d x}{d t} \frac{d^{2} x}{d t^{2}} + c_{v} \frac{d T_{1}}{d t} + c_{v} \frac{d T_{2}}{d t} = 0$

La ecuación del movimiento se escribe

$M \frac{d^{2} x}{d t^{2}} = R T_{10} \frac{x_{0}^{γ - 1}}{x^{γ}} - R T_{20} \frac{{(L - x_{0})}^{γ - 1}}{{(L - x)}^{γ}}$

Se resuelve esta ecuación diferencial de segundo orden por procedimiento numéricos con las siguientes condiciones iniciales: en el instante t=0, x=x₀, y dx/dt=0.

Actividades

Se pulsa el botón titulado Inicio

Se introduce

La temperatura inicial T₁₀ en ºC del gas situado en la parte izquierda del cilindro, actuando en la barra de desplazamiento titulada Temperatura 1
La temperatura inicial T₂₀ en ºC del gas situado en la parte derecha del cilindro, actuando en la barra de desplazamiento titulada Temperatura 2
La masa M del émbolo en kg, en el control de edición titulado Masa
Con el puntero del ratón se mueve el émbolo que separa los dos gases a la posición inicial x₀, determinando el volumen inicial de los gases proporcional a x₀ y (1.0-x₀), respectivamente.

Se pulsa el botón titulado Empieza

Se observa el movimiento oscilatorio del émbolo

Los termómetros marcan la temperatura de los dos gases en K: cuando el gas se comprime adiabáticamente se eleva su temperatura, cuando se expande su temperatura disminuye.

En la parte superior del applet, se representa la energía potencial E_p(x) del émbolo y la energía total mediante una recta horizontal, señalándose los dos puntos de intersección x₀ y x₁ que son los puntos de retorno en los que la velocidad del émbolo se anula. El mínimo de la curva E_p(x) señala la posición de equilibrio, la fuerza sobre el émbolo es nula, la velocidad del émbolo es máxima.

El contador de tiempo en la parte superior del applet, nos permite medir el periodo de las oscilaciones.

Mueva el émbolo con el puntero del ratón