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

En esta página, se estudia la evolución hacia el estado de equilibrio de un sistema aislado formado por dos gases ideales contenidos en un cilindro de sección S separados por una pared adiabática. Inicialmente el émbolo está fijo. Las características de los dos gases son las siguientes:

En el instante t=0, se libera el émbolo. El sistema alcanzará el equilibrio cuando las presiones de los gases separados por el émbolo en reposo sean las mismas p1f=p2f=pf

  1. Aplicamos la ecuación de los gases ideales a cada una de las partes
  2. pf·V1f=n1RT1f

    pf·V2f=n2RT2f

  3. Como el sistema es aislado, la energía total permanece constante
  4. U1f+U2f=U10+U20

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

    n1cv1T1f+n2cv2T2f= n1cv1T10+n2cv2T20

  5. Los volúmenes de los gases cambian, pero el volumen total es constante e igual al inicial que ocupaban los gases.
  6. V10+V20=V1f+V2f

Ahora bien, de estas ecuaciones no podemos despejar los volúmenes finales ni las temperaturas finales. Es preciso hacer ciertas hipótesis acerca del movimiento del émbolo y su interacción con las moléculas del gas.

Movimiento del émbolo

El émbolo de masa m se mueve cuando hay una la diferencia de presión entre los dos gases que separa.

(p1-p2)S

Supondremos además, que el émbolo experimenta una fuerza de rozamiento proporcional a su velocidad y de sentido contrario a ésta.

f r =e dx' dt

donde e es una constante de proporcionalidad y x’=x-x0 es la posición del émbolo respecto de la posición de equilibrio x0. Cuando se alcanza el estado de equilibrio el volumen final del gas 1 es V1f.

La ecuación del movimiento del émbolo se escribe

m d 2 x' d t 2 =( p 1 p 2 )Se dx' dt x'= V 1 V 1f S m d 2 V 1 d t 2 =( p 1 p 2 ) S 2 e d V 1 dt d 2 V 1 d t 2 = S 2 m ( k 1 ( V 1 ) γ 1 k 2 ( V T V 1 ) γ 2 ) e m d V 1 dt

Donde hemos empleado la ecuación de la transformación adiabática

p 1 V 1 γ 1 = k 1 p 2 V 2 γ 2 = k 2 V T = V 1 + V 2

Hemos obtenido la ecuación diferencial que nos da el volumen V1 del gas 1 en función del tiempo.

k1 y k2 serían constantes si no se tuviese en cuenta el rozamiento, pero no pueden serlo si lo hay como vamos a ver a continuación.

Trabajo de la fuerza de rozamiento

El trabajo realizado por la fuerza de rozamiento se transforma en calor que se reparte entre los dos recipientes. En el tiempo dt se genera el calor

f r dx'=e dx' dt dx'= e S 2 ( d V 1 dt )d V 1

Una fracción α de este calor se transfiere al gas 1 y el resto, (1-α) al gas 2.

Vamos ahora a estudiar el proceso por el cual el gas en el intervalo de tiempo Δt, cambia su volumen de V a VV. Este proceso lo podemos dividir por conveniencia en dos etapas:

En la transformación adiabática A→B

p V γ = p B ( V+ΔV ) γ =k

donde k es la constante que caracteriza el proceso adiabático

En la transformación isócora B→C,

ΔQ=n c v ( T C T B )ΔQ=α e S 2 ( dV dt )ΔV

El trabajo de la fuerza de rozamiento se transforma en calor, una parte del cual se transfiere al gas, que incrementa su energía interna, es decir su temperatura.

Al finalizar el proceso la nueva constante k’ que caracteriza la etapa adiabática del  siguiente proceso que comienza en el instante tt  valdrá

p C ( V+ΔV ) γ =k'

Vamos a relacionar k’ con k teniendo en cuenta que ΔV<<V

T C T B = ΔQ n c v ( p C p B )(V+ΔV)= R c v ΔQ ( k'k ) ( V+ΔV ) γ (V+ΔV)= R c v ΔQ Δk=α eR c v S 2 ( dV dt ) ( V+ΔV ) γ1 ΔV

En el límite, cuando el intervalo de tiempo Δt→0

dk dt =α eR c v S 2 ( dV dt ) 2 V γ1

Sistema de ecuaciones diferenciales

El comportamiento del sistema formado por los dos gases y el émbolo de acuerdo con el modelo adoptado, se describe por el siguiente sistema de ecuaciones diferenciales:

d 2 V 1 d t 2 = S 2 m ( k 1 ( V 1 ) γ 1 k 2 ( V T V 1 ) γ 2 ) e m d V 1 dt d k 1 dt =α eR c v1 S 2 ( d V 1 dt ) 2 ( V 1 ) γ 1 1 d k 2 dt =(1α) eR c v2 S 2 ( d V 1 dt ) 2 ( V T V 1 ) γ 2 1

Resolvemos el sistema de ecuaciones diferenciales por procedimientos numéricos, con las siguientes condiciones iniciales.

En el instante t=0.

V 1 = V 10 ( d V 1 dt )=0 k 1 = p 10 V 10 γ 1 k 2 = p 20 V 20 γ 2

Como el sistema formado por los dos gases y el émbolo es aislado la energía total permanecerá constante e igual a la inicial.

n 1 c v1 T 1 + n 2 c v2 T 2 + 1 2 m v 2 = n 1 c v1 T 10 + n 2 c v2 T 20

donde v=dx’/dt es la velocidad del émbolo

n 1 c v1 T 1 + n 2 c v2 T 2 + 1 2 m S 2 ( d V 1 dt ) 2 = n 1 c v1 T 10 + n 2 c v2 T 20

La constancia de la energía del sistema en cada instante es el criterio que nos permitirá probar la adecuación del procedimiento numérico para la resolución del sistema de ecuaciones diferenciales.

Actividades

Se introduce

Se pueden seleccionar cuatro modos de visualización:

Se pulsa el botón titulado Nuevo

Comprobar que si no hay rozamiento, el émbolo oscila permanentemente, no se alcanza el estado de equilibrio.

Si hay rozamiento, el émbolo describe una especie de oscilación amortiguada, al cabo de un cierto tiempo se alcanza el equilibrio igualándose las presiones de los dos gases a ambos lados del émbolo.


||
||
||
||

Referencias

Anacleto J., Anacleto J. A. C. Thermodynamical interactions: subtleties of heat and work concepts. Eur. J. Phys. 29 (2008) 555-566