Решённая задача по термодинамике

Полупроницаемая мембрана: какие величины выравниваются?

Задача 25 · Урок 6 — Второе начало термодинамики

  • равновесие
  • полупроницаемая мембрана
  • химический потенциал
  • парциальное давление
  • неподвижная перегородка
  • подвижная перегородка

Условие

Жёсткий изолированный сосуд объёмом 2V02V_0 содержит две смеси одноатомных идеальных газов A и B. Их разделяет диатермическая мембрана. Она пропускает A, но не B. Первоначально закреплённая и сделанная непроницаемой, она делит сосуд на равные объёмы V0V_0 при одинаковой температуре T0T_0. Слева находятся 2N02N_0 частиц A и N0N_0 частиц B; справа — N0N_0 частиц A и 2N02N_0 частиц B. Затем включают её избирательную проницаемость. Энергией мембраны и взаимодействиями между частицами пренебрегают.

В этой задаче принимается расширение термодинамического тождества на переменные числа частиц:

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 — химический потенциал на одну частицу компонента ss. Для каждой идеальной смеси также принимаем

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}.

Химический потенциал компонента s∈{A,B}s\in\{A,B\} равен

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

где p∘p^\circ — давление отсчёта.

  1. Нарисуйте сосуд и его содержимое. Сравните начальные полные давления, затем парциальные давления A. По химическим потенциалам предскажите начальное направление переноса A.
  2. Мембрана остаётся неподвижной. Запишите dStot\mathrm dS_{\mathrm{tot}} через разрешённые перераспределения dU1\mathrm dU_1 и dNA,1\mathrm dN_{A,1}, затем выведите условия стационарности. Следует требовать равенства полных давлений или μB\mu_B?
  3. Найдите конечную температуру, числа частиц A и полные давления. При температуре T0T_0 проверьте, что найденное распределение A даёт максимум энтропии относительно переноса A.
  4. Начиная из этого состояния, освободите мембрану. Она движется без трения и дополнительных сил, оставаясь диатермической и избирательно проницаемой. Какое условие добавляется? Найдите новое равновесие: объёмы, числа частиц A, температуру и общее полное давление.
  5. Численный расчёт: N0=NAn0N_0=\mathcal N_A n_0, где n0=1,0 moln_0=1{,}0\ \mathrm{mol}, V0=10 LV_0=10\ \mathrm L и T0=300 KT_0=300\ \mathrm K. Напоминаем соотношение R=NAkBR=\mathcal N_A k_{\mathrm{B}} и даём 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}}.

    Для неподвижной, затем для подвижной мембраны вычислите конечные давления и количества A с каждой стороны (в молях), а также конечные объёмы при подвижной мембране.

Подсказка

Подсказка
При постоянном объёме dS=dU/T−μA dNA/T\mathrm dS=\mathrm dU/T-\mu_A\,\mathrm dN_A/T, если B не проходит через мембрану. То, что получает одна сторона, теряет другая. Если мембрана подвижна, добавьте член (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1.

Подробное решение

Решение

Figure
Рисунок 1

Вопрос 1. Схема различает два компонента и показывает, что проходить может только A. Положим P∗=N0kBT0/V0P_*=N_0k_BT_0/V_0. Каждый отсек первоначально содержит 3N03N_0 частиц:

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_*.

Для парциального давления A учитываем только частицы A: 2N02N_0 слева и N0N_0 справа. Первоначально P1,i=P2,i=3P∗P_{1,i}=P_{2,i}=3P_*, но pA,1,i=2P∗p_{A,1,i}=2P_* и pA,2,i=P∗p_{A,2,i}=P_*. Поэтому μA,1−μA,2=kBT0ln⁡2>0\mu_{A,1}-\mu_{A,2}=k_BT_0\ln2>0: A вначале проходит слева направо. При dNA,1<0\mathrm dN_{A,1}<0 и равных температурах полное изменение энтропии (μA,2−μA,1) dNA,1/T0(\mu_{A,2}-\mu_{A,1})\,\mathrm dN_{A,1}/T_0 положительно. Равенства полных давлений недостаточно, чтобы остановить диффузию. В разности химических потенциалов два члена μA∘(T0)\mu_A^\circ(T_0) сокращаются; логарифмы объединяются в kBT0ln⁡(pA,1,i/pA,2,i)=kBT0ln⁡2k_BT_0\ln(p_{A,1,i}/p_{A,2,i})=k_BT_0\ln2.

Вопрос 2. Полная энергия и полное число A сохраняются; числа B с каждой стороны постоянны. Используя dU2=−dU1\mathrm dU_2=-\mathrm dU_1, dNA,2=−dNA,1\mathrm dN_{A,2}=-\mathrm dN_{A,1} и 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).

Объёмный член равен нулю, поскольку мембрана неподвижна. Складываем:

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}.

Группируя коэффициенты, находим

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}.

Теплопроводность и перенос A разрешены: их вариации можно проверять независимо. Поэтому стационарность требует 1/T1−1/T2=01/T_1-1/T_2=0, то есть T1=T2T_1=T_2, затем (μA,2−μA,1)/T=0(\mu_{A,2}-\mu_{A,1})/T=0, то есть μA,1=μA,2\mu_{A,1}=\mu_{A,2}. Она не требует ни P1=P2P_1=P_2 (мембрана закреплена), ни μB,1=μB,2\mu_{B,1}=\mu_{B,2} (B не проходит).

Вопрос 3. Начальная полная энергия равна 9N0kBT09N_0k_BT_0, полное число частиц остаётся 6N06N_0. Действительно, Ui=(3/2)(3N0+3N0)kBT0=9N0kBT0U_i=(3/2)(3N_0+3N_0)k_BT_0=9N_0k_BT_0. При общей конечной температуре TT имеем Uf=(3/2)(6N0)kBT=9N0kBTU_f=(3/2)(6N_0)k_BT=9N_0k_BT. Поэтому равенство Uf=UiU_f=U_i даёт T=T0T=T_0. Равенство химических потенциалов A принимает вид

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}.

Деление на kBT0k_BT_0 и потенцирование дают pA,1=pA,2p_{A,1}=p_{A,2}. Объёмы равны, поэтому множители kBT0/V0k_BT_0/V_0 сокращаются: NA,1=NA,2N_{A,1}=N_{A,2}. Поскольку NA,1+NA,2=3N0N_{A,1}+N_{A,2}=3N_0, каждое число равно 3N0/23N_0/2. Добавляя частицы 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}.

Таким образом,

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

Крепление мембраны уравновешивает эту разность давлений. Полагая a=NA,1a=N_{A,1} и NA,2=3N0−aN_{A,2}=3N_0-a при постоянных температуре и объёмах, из вопроса 2 получаем

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}.

Энергетический член сокращается, поскольку температуры равны и dU1+dU2=0\mathrm dU_1+\mathrm dU_2=0. Для второго дифференцирования запишем ln⁡(3N0−a)−ln⁡a\ln(3N_0-a)-\ln a: производная равна −1/(3N0−a)−1/a-1/(3N_0-a)-1/a. Поэтому снова получаем

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.

Равное распределение является единственным максимумом на 0<a<3N00<a<3N_0.

Вопрос 4. Стационарность по объёму добавляет P1=P2P_1=P_2: член, добавленный к dStot\mathrm dS_{\mathrm{tot}}, равен (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1. Его коэффициент должен исчезнуть; при T1=T2T_1=T_2 это требует равенства полных давлений. Конечная температура по тому же энергетическому балансу, что в вопросе 3, остаётся T0T_0. Поскольку pA,1=pA,2p_{A,1}=p_{A,2}, равенство полных давлений влечёт pB,1=pB,2p_{B,1}=p_{B,2}. B не проходит:

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

поэтому V2=2V1V_2=2V_1, затем 3V1=2V03V_1=2V_0: V1=2V0/3V_1=2V_0/3 и V2=4V0/3V_2=4V_0/3. Для A равенство парциальных давлений даёт NA,1/V1=NA,2/V2N_{A,1}/V_1=N_{A,2}/V_2, поэтому NA,2=2NA,1N_{A,2}=2N_{A,1}. Поскольку сумма равна 3N03N_0, получаем NA,1=N0N_{A,1}=N_0 и NA,2=2N0N_{A,2}=2N_0. Наконец, Pf=(N0+N0)kBT0/(2V0/3)=3P∗P_f=(N_0+N_0)k_BT_0/(2V_0/3)=3P_*. Итак, результаты:

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_*.

Равенство парциальных давлений B здесь следует из других условий и идеальной модели, а не задаётся обменом B. Количества A уже не равны: выравниваются концентрации, а значит и парциальные давления.

Вопрос 5. Для расчёта давлений переводим объёмы в единицы СИ: V0=10 L=10−2 m3V_0=10\ \mathrm L=10^{-2}\ \mathrm{m^3}. Постоянная NA\mathcal N_A переводит количество вещества в число частиц: N0=NAn0N_0=\mathcal N_A n_0 при n0=1,0 moln_0=1{,}0\ \mathrm{mol}. Сначала вычисляем

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}}.

При N0kB=n0RN_0k_B=n_0R масштаб давления, введённый в вопросе 1, равен

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},

поскольку 1 bar=105 Pa1\ \mathrm{bar}=10^5\ \mathrm{Pa}. Это давление служит для расчёта; само по себе оно не является полным давлением отсека.

Неподвижная мембрана. Каждый объём остаётся равным 10 L10\ \mathrm L. Вопрос 3 даёт NA,1=NA,2=3N0/2N_{A,1}=N_{A,2}=3N_0/2. Для перехода от частиц к молям делим на 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 не проходит через мембрану: слева остаётся 1 mol1\ \mathrm{mol} B, справа 2 mol2\ \mathrm{mol}. Поэтому полные количества равны 2,5 mol2{,}5\ \mathrm{mol} и 3,5 mol3{,}5\ \mathrm{mol}. При T0=300 KT_0=300\ \mathrm K уравнение P=nRT/VP=nRT/V даёт по отдельности

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}.

Полные давления различны; закрепление мембраны поддерживает это состояние. Парциальные давления A равны, поскольку A занимает равные объёмы при равных количествах и температурах.

Подвижная мембрана. Из вопроса 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.

Сумма остаётся 20 L20\ \mathrm L, объёмом сосуда. Числа частиц A теперь равны NA,1=N0N_{A,1}=N_0 и 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}.

Снова получаем полное количество 3 mol3\ \mathrm{mol} A. Добавляя B, находим соответственно 22 и 4 mol4\ \mathrm{mol} газа в отсеках. Используем точные объёмы до округления:

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}.

Оба расчёта дают одинаковое полное давление, как и требует механическое равновесие освобождённой мембраны.