Ejercicio resuelto de termodinámica

Membrana semipermeable: ¿qué magnitudes se igualan?

Ejercicio 25 · Lección 6 — El segundo principio de la termodinámica

  • equilibrio
  • membrana semipermeable
  • potencial químico
  • presión parcial
  • pared fija
  • pared móvil

Enunciado

Un recipiente rígido y aislado de volumen 2V02V_0 contiene dos mezclas de gases ideales monoatómicos A y B. Una membrana diatérmana las separa. Deja pasar A, pero no B. Inicialmente inmovilizada e impermeabilizada, divide el recipiente en volúmenes iguales V0V_0 a la misma temperatura T0T_0. A la izquierda hay 2N02N_0 partículas de A y N0N_0 de B; a la derecha, N0N_0 de A y 2N02N_0 de B. Después se activa su permeabilidad selectiva. Se desprecian la energía de la membrana y las interacciones entre partículas.

Para este ejercicio, se admite la extensión de la identidad termodinámica a números de partículas variables:

dU=T dS−P dV+μA dNA+μB dNB.\boxed{\mathrm dU=T\,\mathrm dS-P\,\mathrm dV +\mu_A\,\mathrm dN_A+\mu_B\,\mathrm dN_B.}

μs\mu_s es el potencial químico por partícula de la especie ss. También se admite, para cada mezcla ideal,

U=32(NA+NB)kBT,P=pA+pB,ps=NskBTV.U=\frac32(N_A+N_B)k_BT,\quad P=p_A+p_B,\quad p_s=\frac{N_sk_BT}{V}.

El potencial químico de la especie s∈{A,B}s\in\{A,B\} es

μs(T,ps)=μs∘(T)+kBTln⁡psp∘,\mu_s(T,p_s)=\mu_s^\circ(T)+k_BT\ln\frac{p_s}{p^\circ},

donde p∘p^\circ es una presión de referencia.

  1. Representar el recipiente y sus constituyentes. Comparar las presiones totales iniciales y después las presiones parciales de A. Utilizar los potenciales químicos para predecir el sentido inicial de transferencia de A.
  2. La membrana sigue fija. Expresar dStot\mathrm dS_{\mathrm{tot}} en función de las redistribuciones permitidas dU1\mathrm dU_1 y dNA,1\mathrm dN_{A,1} y deducir las condiciones de estacionariedad. ¿Deben ser iguales las presiones totales o los μB\mu_B?
  3. Determinar la temperatura final, los números de partículas de A y las presiones totales. A la temperatura T0T_0, comprobar que la distribución de A obtenida maximiza la entropía respecto a la transferencia de A.
  4. Partiendo de este estado, se libera la membrana. Se desplaza sin rozamiento ni fuerza adicional y sigue siendo diatérmana y selectiva. ¿Qué condición adicional aparece? Determinar el nuevo equilibrio: volúmenes, números de partículas de A, temperatura y presión total común.
  5. Aplicación numérica: N0=NAn0N_0=\mathcal N_A n_0, con n0=1,0 moln_0=1{,}0\ \mathrm{mol}, V0=10 LV_0=10\ \mathrm L y T0=300 KT_0=300\ \mathrm K. Recordar la relación R=NAkBR=\mathcal N_A k_{\mathrm{B}} y utilizar NA=6,022×1023 mol−1,kB=1,380649×10−23 J K−1.\mathcal N_A=6{,}022\times10^{23}\ \mathrm{mol^{-1}},\qquad k_{\mathrm{B}}=1{,}380649\times10^{-23}\ \mathrm{J\,K^{-1}}.

    Para la membrana fija y después móvil, calcular las presiones finales y las cantidades de A a cada lado (en moles), así como los volúmenes finales cuando la membrana es móvil.

Indicación

Pista
A volumen fijo, dS=dU/T−μA dNA/T\mathrm dS=\mathrm dU/T-\mu_A\,\mathrm dN_A/T si B no atraviesa la membrana. Las transferencias recibidas por un lado se pierden por el otro. Cuando la membrana es móvil, añadir el término (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1.

Solución detallada

Solución

Figure
Figura 1

Pregunta 1. El esquema distingue ambas especies y muestra que solo A puede atravesar la membrana. Sea P∗=N0kBT0/V0P_*=N_0k_BT_0/V_0. Cada compartimento contiene inicialmente 3N03N_0 partículas:

P1,i=(2N0+N0)kBT0V0=3P∗,P2,i=(N0+2N0)kBT0V0=3P∗.P_{1,i}=\frac{(2N_0+N_0)k_BT_0}{V_0}=3P_*,\qquad P_{2,i}=\frac{(N_0+2N_0)k_BT_0}{V_0}=3P_*.

Para la presión parcial de A solo se cuentan sus partículas: 2N02N_0 a la izquierda y N0N_0 a la derecha. Inicialmente P1,i=P2,i=3P∗P_{1,i}=P_{2,i}=3P_*, pero pA,1,i=2P∗p_{A,1,i}=2P_* y pA,2,i=P∗p_{A,2,i}=P_*. Así, μA,1−μA,2=kBT0ln⁡2>0\mu_{A,1}-\mu_{A,2}=k_BT_0\ln2>0: A se desplaza inicialmente de izquierda a derecha. Para dNA,1<0\mathrm dN_{A,1}<0 a temperaturas iguales, la variación de entropía total (μA,2−μA,1) dNA,1/T0(\mu_{A,2}-\mu_{A,1})\,\mathrm dN_{A,1}/T_0 es positiva. La igualdad de presiones totales no basta para detener la difusión. En la diferencia de potenciales químicos, los dos términos μA∘(T0)\mu_A^\circ(T_0) se cancelan; agrupando los logaritmos se obtiene kBT0ln⁡(pA,1,i/pA,2,i)=kBT0ln⁡2k_BT_0\ln(p_{A,1,i}/p_{A,2,i})=k_BT_0\ln2.

Pregunta 2. La energía total y el número total de partículas de A se conservan; los números de partículas de B a cada lado son constantes. Utilizando dU2=−dU1\mathrm dU_2=-\mathrm dU_1, dNA,2=−dNA,1\mathrm dN_{A,2}=-\mathrm dN_{A,1} y dNB,1=dNB,2=0\mathrm dN_{B,1}=\mathrm dN_{B,2}=0,

dSj=dUjTj−μA,jTj dNA,j(j=1,2).\mathrm dS_j=\frac{\mathrm dU_j}{T_j} -\frac{\mu_{A,j}}{T_j}\,\mathrm dN_{A,j} \quad (j=1,2).

El término de volumen se anula porque la membrana está fija. Sumando,

dStot=dU1T1−dU1T2−μA,1T1dNA,1+μA,2T2dNA,1.\mathrm dS_{\mathrm{tot}} =\frac{\mathrm dU_1}{T_1}-\frac{\mathrm dU_1}{T_2} -\frac{\mu_{A,1}}{T_1}\mathrm dN_{A,1} +\frac{\mu_{A,2}}{T_2}\mathrm dN_{A,1}.

Agrupando los coeficientes,

dStot=(1T1−1T2)dU1+(μA,2T2−μA,1T1)dNA,1.\mathrm dS_{\mathrm{tot}}= \left(\frac1{T_1}-\frac1{T_2}\right)\mathrm dU_1 +\left(\frac{\mu_{A,2}}{T_2}-\frac{\mu_{A,1}}{T_1}\right) \mathrm dN_{A,1}.

La conducción térmica y la transferencia de A están permitidas: sus variaciones pueden examinarse independientemente. La estacionariedad exige entonces 1/T1−1/T2=01/T_1-1/T_2=0, es decir, T1=T2T_1=T_2, y después (μA,2−μA,1)/T=0(\mu_{A,2}-\mu_{A,1})/T=0, es decir, μA,1=μA,2\mu_{A,1}=\mu_{A,2}. No exige P1=P2P_1=P_2 (la membrana está fija) ni μB,1=μB,2\mu_{B,1}=\mu_{B,2} (B no puede atravesarla).

Pregunta 3. Inicialmente, la energía total es 9N0kBT09N_0k_BT_0 y el número total de partículas sigue siendo 6N06N_0. En efecto, Ui=(3/2)(3N0+3N0)kBT0=9N0kBT0U_i=(3/2)(3N_0+3N_0)k_BT_0=9N_0k_BT_0. A una temperatura final común TT, Uf=(3/2)(6N0)kBT=9N0kBTU_f=(3/2)(6N_0)k_BT=9N_0k_BT. La igualdad Uf=UiU_f=U_i da entonces T=T0T=T_0. La igualdad de potenciales químicos de A se escribe

kBT0ln⁡pA,1p∘=kBT0ln⁡pA,2p∘.k_BT_0\ln\frac{p_{A,1}}{p^\circ} =k_BT_0\ln\frac{p_{A,2}}{p^\circ}.

Dividiendo por kBT0k_BT_0 y exponenciando se obtiene pA,1=pA,2p_{A,1}=p_{A,2}. Como los volúmenes son iguales, los factores kBT0/V0k_BT_0/V_0 se cancelan: NA,1=NA,2N_{A,1}=N_{A,2}. Como NA,1+NA,2=3N0N_{A,1}+N_{A,2}=3N_0, cada uno vale 3N0/23N_0/2. Añadiendo las partículas de B,

P1,f=(3N0/2+N0)kBT0V0,P2,f=(3N0/2+2N0)kBT0V0.P_{1,f}=\frac{(3N_0/2+N_0)k_BT_0}{V_0},\qquad P_{2,f}=\frac{(3N_0/2+2N_0)k_BT_0}{V_0}.

Así,

P1,f=52P∗,P2,f=72P∗.P_{1,f}=\frac52P_*,\qquad P_{2,f}=\frac72P_*.

El soporte de la membrana compensa esta diferencia de presión. Definiendo a=NA,1a=N_{A,1} y NA,2=3N0−aN_{A,2}=3N_0-a, a temperatura y volúmenes fijos, la pregunta 2 da

dStotda=μA,2−μA,1T0=kBln⁡pA,2pA,1=kBln⁡3N0−aa.\frac{\mathrm dS_{\mathrm{tot}}}{\mathrm da} =\frac{\mu_{A,2}-\mu_{A,1}}{T_0} =k_{\mathrm{B}}\ln\frac{p_{A,2}}{p_{A,1}} =k_{\mathrm{B}}\ln\frac{3N_0-a}{a}.

El término energético se anula porque las temperaturas son iguales y dU1+dU2=0\mathrm dU_1+\mathrm dU_2=0. Para derivar otra vez, se escribe ln⁡(3N0−a)−ln⁡a\ln(3N_0-a)-\ln a: su derivada es −1/(3N0−a)−1/a-1/(3N_0-a)-1/a. Así,

dStotda=kBln⁡3N0−aa,d2Stotda2=−kB(1a+13N0−a)<0.\frac{\mathrm dS_{\mathrm{tot}}}{\mathrm da} =k_{\mathrm{B}}\ln\frac{3N_0-a}{a},\qquad \frac{\mathrm d^2S_{\mathrm{tot}}}{\mathrm da^2} =-k_{\mathrm{B}}\left(\frac1a+\frac1{3N_0-a}\right)<0.

La distribución equitativa es el único máximo en 0<a<3N00<a<3N_0.

Pregunta 4. La estacionariedad respecto al volumen añade P1=P2P_1=P_2: el término adicional en dStot\mathrm dS_{\mathrm{tot}} es (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1. Su coeficiente debe anularse; con T1=T2T_1=T_2, esto exige la igualdad de presiones totales. La temperatura final vuelve a ser T0T_0, por el mismo balance energético de la pregunta 3. Como pA,1=pA,2p_{A,1}=p_{A,2}, la igualdad de presiones totales implica pB,1=pB,2p_{B,1}=p_{B,2}. B no puede atravesar la membrana:

N0V1=2N0V2,V1+V2=2V0,\frac{N_0}{V_1}=\frac{2N_0}{V_2},\qquad V_1+V_2=2V_0,

por lo que V2=2V1V_2=2V_1 y después 3V1=2V03V_1=2V_0: V1=2V0/3V_1=2V_0/3 y V2=4V0/3V_2=4V_0/3. Para A, la igualdad de presiones parciales da NA,1/V1=NA,2/V2N_{A,1}/V_1=N_{A,2}/V_2, de donde NA,2=2NA,1N_{A,2}=2N_{A,1}. Como la suma vale 3N03N_0, se obtiene NA,1=N0N_{A,1}=N_0 y NA,2=2N0N_{A,2}=2N_0. Finalmente, Pf=(N0+N0)kBT0/(2V0/3)=3P∗P_f=(N_0+N_0)k_BT_0/(2V_0/3)=3P_*. Los resultados son entonces

V1=23V0,V2=43V0,NA,1=N0,NA,2=2N0,Pf=3P∗.V_1=\frac23V_0,\quad V_2=\frac43V_0,\quad N_{A,1}=N_0,\quad N_{A,2}=2N_0,\quad P_f=3P_*.

La igualdad de presiones parciales de B se deduce aquí de las otras condiciones y del modelo de mezcla ideal, no de un intercambio permitido de B. Las cantidades de A ya no son iguales: se igualan sus concentraciones y, por tanto, sus presiones parciales.

Pregunta 5. Se convierten los volúmenes a unidades SI para calcular las presiones: V0=10 L=10−2 m3V_0=10\ \mathrm L=10^{-2}\ \mathrm{m^3}. La constante NA\mathcal N_A permite convertir una cantidad de sustancia en un número de partículas: N0=NAn0N_0=\mathcal N_A n_0 con n0=1,0 moln_0=1{,}0\ \mathrm{mol}. Primero se calcula

R=NAkB=(6,022×1023)(1,380649×10−23)≃8,314 J mol−1 K−1.R=\mathcal N_A k_{\mathrm{B}} =(6{,}022\times10^{23})(1{,}380649\times10^{-23}) \simeq8{,}314\ \mathrm{J\,mol^{-1}\,K^{-1}}.

Con N0kB=n0RN_0k_B=n_0R, la escala de presión introducida en la pregunta 1 es

P∗=N0kBT0V0=n0RT0V0≃1,0×8,314×30010−2=249420 Pa≃2,494 bar,P_* =\frac{N_0k_BT_0}{V_0} =\frac{n_0RT_0}{V_0} \simeq\frac{1{,}0\times8{,}314\times300}{10^{-2}} =249420\ \mathrm{Pa}\simeq2{,}494\ \mathrm{bar},

pues 1 bar=105 Pa1\ \mathrm{bar}=10^5\ \mathrm{Pa}. Esta presión es una escala de cálculo; no es por sí misma la presión total de un compartimento.

Membrana fija. Cada volumen sigue siendo 10 L10\ \mathrm L. La pregunta 3 da NA,1=NA,2=3N0/2N_{A,1}=N_{A,2}=3N_0/2. Para convertir partículas en moles, se divide por NA\mathcal N_A:

nA,1=nA,2=3N02NA=32×1 mol=1,5 mol.n_{A,1}=n_{A,2}=\frac{3N_0}{2\mathcal N_A} =\frac32\times1\ \mathrm{mol}=1{,}5\ \mathrm{mol}.

B no atraviesa la membrana: quedan 1 mol1\ \mathrm{mol} de B a la izquierda y 2 mol2\ \mathrm{mol} a la derecha. Las cantidades totales son entonces 2,5 mol2{,}5\ \mathrm{mol} y 3,5 mol3{,}5\ \mathrm{mol}. A T0=300 KT_0=300\ \mathrm K, la ecuación P=nRT/VP=nRT/V da por separado

P1,f=(1,5+1)×8,314×30010−2=623550 Pa≃6,24 bar,P2,f=(1,5+2)×8,314×30010−2=872970 Pa≃8,73 bar.P_{1,f}=\frac{(1{,}5+1)\times8{,}314\times300}{10^{-2}} =623550\ \mathrm{Pa}\simeq6{,}24\ \mathrm{bar}, P_{2,f}=\frac{(1{,}5+2)\times8{,}314\times300}{10^{-2}} =872970\ \mathrm{Pa}\simeq8{,}73\ \mathrm{bar}.

Las presiones totales difieren; la inmovilización de la membrana mantiene este estado. Las presiones parciales de A son iguales porque A ocupa volúmenes iguales con cantidades y temperaturas iguales.

Membrana móvil. Según la pregunta 4,

V1,f=23V0=23×10≃6,67 L,V2,f=43V0=43×10≃13,33 L.V_{1,f}=\frac23V_0=\frac23\times10\simeq6{,}67\ \mathrm L, \qquad V_{2,f}=\frac43V_0=\frac43\times10\simeq13{,}33\ \mathrm L.

La suma sigue siendo 20 L20\ \mathrm L, el volumen del recipiente. Las poblaciones de A son ahora NA,1=N0N_{A,1}=N_0 y NA,2=2N0N_{A,2}=2N_0:

nA,1=N0NA=1 mol,nA,2=2N0NA=2 mol.n_{A,1}=\frac{N_0}{\mathcal N_A}=1\ \mathrm{mol},\qquad n_{A,2}=\frac{2N_0}{\mathcal N_A}=2\ \mathrm{mol}.

El total es, en efecto, 3 mol3\ \mathrm{mol} de A. Añadiendo B, los compartimentos contienen respectivamente 22 y 4 mol4\ \mathrm{mol} de gas. Se utilizan los volúmenes exactos antes de redondear:

P1,f=(1+1)×8,314×300(2/3)×10−2=748260 Pa≃7,48 bar,P2,f=(2+2)×8,314×300(4/3)×10−2=748260 Pa≃7,48 bar.P_{1,f}=\frac{(1+1)\times8{,}314\times300}{(2/3)\times10^{-2}} =748260\ \mathrm{Pa}\simeq7{,}48\ \mathrm{bar}, P_{2,f}=\frac{(2+2)\times8{,}314\times300}{(4/3)\times10^{-2}} =748260\ \mathrm{Pa}\simeq7{,}48\ \mathrm{bar}.

Ambos cálculos dan la misma presión total, como exige el equilibrio mecánico de la membrana liberada.