Fundamentos estadísticos de la transformación adiabática

Una molécula de masa m se mueve con velocidad u, cuya componente horizontal es ux, choca con el émbolo de masa M que se mueve hacia la derecha con velocidad v.

Supondremos que las moléculas chocan elásticamente con el émbolo y las paredes del recipiente. Esta suposición es razonable, si consideramos que en un proceso adiabático no hay intercambio de energía en forma de calor a través de las paredes del recipiente con el entorno

Aplicamos la conservación del momento lineal y la energía cinética

{ m u x +Mv=m u x ' +Mv' 1 2 m u x 2 + 1 2 M v 2 = 1 2 m ( u x ' ) 2 + 1 2 M ( v' ) 2

Despejaremos las velocidades después del choque en este sistema de dos ecuaciones

{ m( u x u x ' )=M( v'v ) m( u x 2 ( u x ' ) 2 )=M( ( v' ) 2 v 2 )

Simplificamos la segunda ecuación

m( u x u x ' )( u x + u x ' )=M( v'v )( v'+v )

Resolvemos el sistema de dos ecuaciones con dos incógnitas

{ u x + u x ' =v'+v m( u x u x ' )=M( v'v )

El resultado es

v'= Mm M+m v+2 m m+M u x = 1 m M 1+ m M v+ m M 1+ m M u x u x ' =2 M m+M v Mm M+m u x =2 1 1+ m M v 1 m M 1+ m M u x

La masa de una molécula es mucho menor que la masa del émbolo, m<<M y también, la velocidad media de las moléculas es mucho mayor que la velocidad del desplazamiento del émbolo ux>>v

{ v'v u x ' 2v u x

Se cumple además

( u x ' ) 2 ( u x ) 2 4v u x

En cada impacto hay un cambio de energía.

Δ( 1 2 m u x 2 )= 1 2 m ( u x ' ) 2 1 2 m ( u x ) 2 = 1 2 m·4v u x =2mv u x

La energía cinética media de las moléculas de un gas ideal es

1 2 m u 2 = 1 2 m u x 2 + 1 2 m u y 2 + 1 2 m u z 2 = 3 2 kT

La energía E de un mol de gas se relaciona con la temperatura T

E= N A 1 2 m u 2 = 3 2 N A kT= 3 2 RT

NA es el número de Avogadro

En el intervalo de tiempo comprendido entre t y t+Δt chocan con el émbolo, las moléculas contenidas en el volumen

1 2 n( A u x Δt )

Siendo A el área del émbolo y n es el número de moléculas por unidad de volumen. El factor 1/2 se refiere a que de las moléculas contenidas en este volumen, la mitad (en promedio) se dirigen hacia el émbolo ux>0 y la otra mitad, se alejan del émbolo ux<0

Para un mol de gas, el número de moléculas por unidad de volumen es

n= N A Ax

En el intervalo de tiempo comprendido entre t y t+Δt, chocan con el émbolo

1 2 ( N A Ax )( A u x Δt )= 1 2 N A u x x Δt

El cambio de energía del gas ideal en el intervalo de tiempo comprendido entre t y t+Δt es el producto del número de choques por el cambio de energía en cada choque

ΔE=( 1 2 N A u x x Δt )Δ( 1 2 m u x 2 )= 1 2 N A u x x Δt( 2mv u x )=2 N A 1 2 m u x 2 v x Δt=2 N A ( 1 2 kT ) Δx x = N A kT Δx x

En el límite cuando Δt→0

3 2 N A kdT= N A kT dx x dT T = 2 3 dx x dT T = 2 3 dV V

El volumen del recipiente cilíndrico que contiene un mol de gas ideal, es V=A·x. Integramos

lnT+ln V 2/3 =cte T V 2/3 =cte

Gas ideal formado por una sola partícula

La pared fija está situada en x=0 y la pared móvil está situada inicialmente en la posición x0 en el instante t=0. La pared móvil se mueve con velocidad constante u hacia la pared fija.

La posición de la pared móvil en el instante t es x=x0-u·t. El movimiento de la partícula, cuya velocidad constante es v0, encerrada entre ambas paredes se muestra en la figura y se describe a continuación

Primer choque:

La partícula parte en el instante t=0, de la posición x0, se mueve hacia la pared fija con velocidad constante v0. Choca con la pared fija cambiando el sentido de su velocidad y se dirige al encuentro de la pared móvil, que se produce en el instante t1.

t 1 = 2 x 0 u+ v 0

La posición de la pared móvil en dicho instante es x1=x0-u·t1

Segundo choque:

Después del choque, velocidad de la partícula es v1=v0+2u y se dirige hacia la pared fija. Choca con la pared fija, cambiando el sentido de su velocidad y se dirige al encuentro de la pared móvil, que se produce en el instante t2.

t 2 = t 1 + 2 x 1 u+ v 1

La posición de la pared móvil en dicho instante es x2=x0-u·t2

Tercer choque:

Después del choque, la velocidad de la partícula es v2=v1+2u=v0+2(2u) y se dirige hacia la pared fija. Choca con la pared fija, cambiando el sentido de su velocidad y se dirige al encuentro de la pared móvil, que se produce en el instante t3.

t 3 = t 2 + 2 x 2 u+ v 2

La posición de la pared móvil en dicho instante es x3=x0-u·t3

Choque n+1:

Después del choque n la velocidad de la partícula es vn=vn-1+2u=v0+n(2u) y se dirige hacia la pared fija. Choca con la pared fija, cambiando el sentido de su velocidad y se dirige al encuentro de la pared móvil, que se produce en el instante tn+1.

t n+1 = t n + 2 x n u+ v n

La posición de la pared móvil en dicho instante es xn+1=x0-u·tn+1

La velocidad de la partícula va creciendo, con el número n de choques, alcanzando el valor

vn =v0+2nu

El desplazamiento de la pared móvil entre dos choques sucesivos que se producen en los instantes tn y tn+1 es

x n+1 x n = x 0 u( t n + 2 x n u+ v n )( x 0 u t n )= 2 x n u u+ v 0 +2nu x n+1 x n = 2 x n (2n+1)+ v 0 /u

De la formulación discreta a la continua.

El parámetro macroscópico observable es la posición x de la pared móvil o el volumen del recipiente que contiene la partícula. Expresaremos la fuerza que ejerce la pared sobre la partícula, el trabajo y la energía cinética de la partícula en función de x.

Como x es una función de n podemos escribir de forma aproximada

dx(n) dn = 2x(n) (2n+1)+ v 0 /u

Separando variables e integrando

x 0 x dx 2x = 0 n dn (2n+1)+ v 0 /u 1 2 ln x x 0 = 1 2 ln 2n+1+ v 0 /u 1+ v 0 /u

Despejando x

x= x 0 (1+ v 0 /u) 2n+1+ v 0 /u

El desplazamiento Δx de la pared móvil durante el intervalo de tiempo entre dos choques consecutivos es

x n+1 x n Δx= 2x 2n+1+ v 0 /u

La posición x de la pared móvil es función del número de choques n. La velocidad de la partícula v es también función del número de choques n, v=v0+2nu. Expresamos la velocidad v de la partícula en función de la posición x de la pared móvil

v=(u+ v 0 ) x 0 x u

Expresamos, igualmente el desplazamiento Δx en función de x

Δx= 2xu v+u = 2 x 2 u (u+ v 0 ) x 0

A partir de la expresión de x en función del número n de choques, despejamos n y derivamos respecto del tiempo.

2 dn dt = x 0 (1+ v 0 /u) x 2 dx dt

Teniendo en cuenta que dx/dt=-u, la frecuencia f de los choques (número de choques en la unidad de tiempo) es

f= dn dt = x 0 (u+ v 0 ) 2 x 2

El incremento de la frecuencia entre los choques se debe a dos causas:

Trabajo de la fuerza F que ejerce la pared móvil

El cambio de momento lineal en cada choque es  m(v+2u)-m(-v)=2m(u+v)

La fuerza F que ejerce la pared móvil se emplea en cambiar el momento lineal p de la partícula F=dp/dt. La fuerza F es el producto de la frecuencia de los choques por el cambio de momento lineal en cada choque.

F=2m(u+v)f=m (u+ v 0 ) 2 x 0 2 x 3

El trabajo de la fuerza F cuando la pared móvil se desplaza Δx es

F·Δx=m (u+ v 0 ) 2 x 0 2 x 3 · 2 x 2 u (u+ v 0 ) x 0 =2mu(u+ v 0 ) x 0 x

La variación de energía cinética en un choque es

Δ E k = 1 2 m ( v+2u ) 2 1 2 m v 2 =2mu(u+ v 0 ) x 0 x

El trabajo de de la fuerza F se invierte en incrementar la energía cinética de la partícula.

-F·ΔxEk

Relación entre las magnitudes microscópicas y los parámetros macroscópicos.

La energía cinética está relacionada con la temperatura, la fuerza con la presión y la posición de la pared móvil con el volumen del recipiente. Sustituimos:

-p·dV=k·dT

La ecuación de estado de los gases ideales relaciona, presión p, el volumen V y temperatura absoluta T.

nRT V dV=kdT dT T +(γ-1) dV V =0

Integrando miembro a miembro, obtenemos la ecuación de una transformación adiabática.

T V γ1 =cteP V γ =cte

donde γ al igual que k son constantes.

Actividades

Se introduce

Se pulsa en el botón titulado Nuevo

Se observa el movimiento de la partícula. En la parte superior izquierda, se proporciona el dato de la velocidad v de la partícula en cada instante t.

En la parte superior, se representan en dos gráficas la velocidad v de la partícula en función de la posición x del émbolo móvil.

Referencias

E N Miranda. Adiabatic reversible compression: a molecular view. Eur. J. Phys. 23 (2002) pp. 389–393

Vermillon, R. Compresing a classical particle in a 1-D box. Am. J. Phys. 40, September 1978, pp. 1340-1342