Çözümlü termodinamik alıştırması

Yarı geçirgen zar: hangi nicelikler eşitlenir?

Alıştırma 25 · Ders 6 — Termodinamiğin İkinci Yasası

  • denge
  • yarı geçirgen zar
  • kimyasal potansiyel
  • kısmi basınç
  • sabit bölme
  • hareketli bölme

Soru

2V02V_0 hacimli rijit ve yalıtılmış kap, tek atomlu ideal A veB gazlarının iki karışımını içerir. Isı geçiren bir zar bunları ayırır. A'yı geçirir, B'yi geçirmez. Başlangıçta sabit tutulan ve geçirimsiz yapılan zar, kabı aynı T0T_0 sıcaklığında eşit V0V_0 hacimlerine böler. Solda 2N02N_0 A ve N0N_0 B parçacığı; sağda N0N_0 A ve 2N02N_0 B parçacığı vardır. Sonra seçici geçirgenlik etkinleştirilir. Zarın enerjisini ve parçacık etkileşimlerini ihmal ediniz.

Bu alıştırmada termodinamik özdeşliğin değişken parçacık sayılarına genişletilmesini kabul ediniz:

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 türünün parçacık başına kimyasal potansiyelidir. Her ideal karışım için ayrıca,

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\} türünün kimyasal potansiyeli,

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

burada p∘p^\circ referans basınçtır.

  1. Kabı ve bileşenlerini çiziniz. Önce başlangıç toplam basınçlarını, sonra A'nın kısmi basınçlarını karşılaştırınız. Kimyasal potansiyellerle A'nın başlangıç aktarım yönünü öngörünüz.
  2. Zar sabit kalır. dStot\mathrm dS_{\mathrm{tot}} ifadesini izin verilen dU1\mathrm dU_1 ve dNA,1\mathrm dN_{A,1} yeniden dağılımları cinsinden yazıp durağanlık koşullarını bulunuz. Toplam basınçların veya μB\mu_B değerlerinin eşit olması gerekir mi?
  3. Son sıcaklığı, A parçacık sayılarını ve toplam basınçları bulunuz. T0T_0 sıcaklığında bulunan A dağılımının A aktarımına göre entropiyi maksimum yaptığını doğrulayınız.
  4. Bu durumdan başlayıp zarı serbest bırakalım. Sürtünmesiz ve ek kuvvetsiz hareket eder, ısı geçiren ve seçici kalır. Hangi ek koşul geçerlidir? Yeni dengeyi bulunuz: hacimler, A parçacık sayıları, sıcaklık ve ortak toplam basınç.
  5. Sayısal uygulama: N0=NAn0N_0=\mathcal N_A n_0, n0=1,0 moln_0=1{,}0\ \mathrm{mol}, V0=10 LV_0=10\ \mathrm L ve T0=300 KT_0=300\ \mathrm K. R=NAkBR=\mathcal N_A k_{\mathrm{B}} bağıntısını hatırlayıp şunu kullanınız: 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}}.

    Sabit ve ardından hareketli zar için son basınçları, her iki taraftaki A miktarını (mol olarak) ve zar hareketliyken son hacimleri hesaplayınız.

İpucu

İpucu
Sabit hacimde B zarı geçmiyorsa dS=dU/T−μA dNA/T\mathrm dS=\mathrm dU/T-\mu_A\,\mathrm dN_A/T. Bir tarafın aldığı aktarımlar diğerinden eksilir. Zar hareketliyken (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1 terimini ekleyiniz.

Ayrıntılı çözüm

Çözüm

Figure
Şekil 1

Soru 1. Şema iki türü ayırır ve yalnızca A'nın geçebildiğini gösterir. P∗=N0kBT0/V0P_*=N_0k_BT_0/V_0 tanımlayalım. Başlangıçta her bölüm 3N03N_0 parçacık içerir:

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'nın kısmi basıncı için yalnızca A parçacıklarını sayarız: solda 2N02N_0, sağda N0N_0. Başlangıçta P1,i=P2,i=3P∗P_{1,i}=P_{2,i}=3P_*, fakat pA,1,i=2P∗p_{A,1,i}=2P_* ve pA,2,i=P∗p_{A,2,i}=P_*. Böylece μA,1−μA,2=kBT0ln⁡2>0\mu_{A,1}-\mu_{A,2}=k_BT_0\ln2>0: A başlangıçta soldan sağa geçer. Eşit sıcaklıklarda dNA,1<0\mathrm dN_{A,1}<0 için toplam (μA,2−μA,1) dNA,1/T0(\mu_{A,2}-\mu_{A,1})\,\mathrm dN_{A,1}/T_0 entropi değişimi pozitiftir. Eşit toplam basınçlar difüzyonu durdurmaya yetmez. Kimyasal potansiyel farkında iki μA∘(T0)\mu_A^\circ(T_0) terimi yok olur; logaritmaları birleştirince kBT0ln⁡(pA,1,i/pA,2,i)=kBT0ln⁡2k_BT_0\ln(p_{A,1,i}/p_{A,2,i})=k_BT_0\ln2 bulunur.

Soru 2. Toplam enerji ve toplam A parçacık sayısı korunur; her taraftaki B parçacık sayıları sabittir. dU2=−dU1\mathrm dU_2=-\mathrm dU_1, dNA,2=−dNA,1\mathrm dN_{A,2}=-\mathrm dN_{A,1} ve dNB,1=dNB,2=0\mathrm dN_{B,1}=\mathrm dN_{B,2}=0 kullanılırsa,

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

Zar sabit olduğundan hacim terimi sıfırlanır. Toplayınca,

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

Katsayıları bir araya getirince,

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

Isı iletimi ve A aktarımı serbesttir: değişimleri bağımsız sınanabilir. Bu nedenle durağanlık 1/T1−1/T2=01/T_1-1/T_2=0, yani T1=T2T_1=T_2, ardından (μA,2−μA,1)/T=0(\mu_{A,2}-\mu_{A,1})/T=0, yani μA,1=μA,2\mu_{A,1}=\mu_{A,2} gerektirir. Ne P1=P2P_1=P_2 (zar sabittir) ne de μB,1=μB,2\mu_{B,1}=\mu_{B,2} (B geçemez) gerekir.

Soru 3. Başlangıç toplam enerjisi 9N0kBT09N_0k_BT_0 olup toplam parçacık sayısı 6N06N_0 kalır. Gerçekten Ui=(3/2)(3N0+3N0)kBT0=9N0kBT0U_i=(3/2)(3N_0+3N_0)k_BT_0=9N_0k_BT_0. Ortak son TT sıcaklığında Uf=(3/2)(6N0)kBT=9N0kBTU_f=(3/2)(6N_0)k_BT=9N_0k_BT. Bu yüzden Uf=UiU_f=U_i eşitliği T=T0T=T_0 verir. A'nın kimyasal potansiyellerinin eşitliği,

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 değerine bölüp üstelini alınca pA,1=pA,2p_{A,1}=p_{A,2}. Hacimler eşit olduğundan kBT0/V0k_BT_0/V_0 çarpanları sadeleşir: NA,1=NA,2N_{A,1}=N_{A,2}. NA,1+NA,2=3N0N_{A,1}+N_{A,2}=3N_0 olduğundan her biri 3N0/23N_0/2 olur. B parçacıklarını ekleyince,

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

Böylece,

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

Zarın desteği bu basınç farkını dengeler. Sabit sıcaklık ve hacimlerde a=NA,1a=N_{A,1} ve NA,2=3N0−aN_{A,2}=3N_0-a tanımlanırsa soru 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}.

Sıcaklıklar eşit ve dU1+dU2=0\mathrm dU_1+\mathrm dU_2=0 olduğundan enerji terimi yok olur. Yeniden türev almak için ln⁡(3N0−a)−ln⁡a\ln(3N_0-a)-\ln a yazalım: türevi −1/(3N0−a)−1/a-1/(3N_0-a)-1/a olur. Böylece,

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.

Eşit dağılım 0<a<3N00<a<3N_0 üzerinde tek maksimumdur.

Soru 4. Hacme göre durağanlık P1=P2P_1=P_2 ekler: dStot\mathrm dS_{\mathrm{tot}} ifadesindeki ek terim (P1/T1−P2/T2) dV1(P_1/T_1-P_2/T_2)\,\mathrm dV_1 olur. Katsayısı sıfır olmalıdır; T1=T2T_1=T_2 ile bu, toplam basınçların eşitliğini gerektirir. Soru 3'teki aynı enerji bilançosuyla son sıcaklık yine T0T_0. pA,1=pA,2p_{A,1}=p_{A,2} olduğundan eşit toplam basınçlar pB,1=pB,2p_{B,1}=p_{B,2} gerektirir. B geçemez:

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

dolayısıyla V2=2V1V_2=2V_1, ardından 3V1=2V03V_1=2V_0: V1=2V0/3V_1=2V_0/3 ve V2=4V0/3V_2=4V_0/3. A için eşit kısmi basınçlar NA,1/V1=NA,2/V2N_{A,1}/V_1=N_{A,2}/V_2, dolayısıyla NA,2=2NA,1N_{A,2}=2N_{A,1} verir. Toplam 3N03N_0 olduğundan NA,1=N0N_{A,1}=N_0 ve NA,2=2N0N_{A,2}=2N_0 bulunur. Son olarak Pf=(N0+N0)kBT0/(2V0/3)=3P∗P_f=(N_0+N_0)k_BT_0/(2V_0/3)=3P_*. Sonuçlar,

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 kısmi basınçlarının eşitliği burada izin verilen B alışverişinden değil, diğer koşullardan ve ideal karışım modelinden gelir. A miktarları artık eşit değildir: derişimleri ve dolayısıyla kısmi basınçları eşitlenir.

Soru 5. Basınçları hesaplamak için hacimleri SI birimlerine çevirelim: V0=10 L=10−2 m3V_0=10\ \mathrm L=10^{-2}\ \mathrm{m^3}. NA\mathcal N_A sabiti madde miktarını parçacık sayısına dönüştürür: n0=1,0 moln_0=1{,}0\ \mathrm{mol} ile N0=NAn0N_0=\mathcal N_A n_0. Önce,

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 ile soru 1'deki basınç ölçeği,

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

çünkü 1 bar=105 Pa1\ \mathrm{bar}=10^5\ \mathrm{Pa}. Bu basınç bir hesap ölçeğidir; bir bölümün toplam basıncı değildir.

Sabit zar. Her hacim 10 L10\ \mathrm L kalır. Soru 3 NA,1=NA,2=3N0/2N_{A,1}=N_{A,2}=3N_0/2 verir. Parçacık sayısından mole geçmek için NA\mathcal N_A değerine böleriz:

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 zarı geçemez: solda 1 mol1\ \mathrm{mol}, sağda 2 mol2\ \mathrm{mol} B kalır. Toplam miktarlar böylece 2,5 mol2{,}5\ \mathrm{mol} ve 3,5 mol3{,}5\ \mathrm{mol} olur. T0=300 KT_0=300\ \mathrm K sıcaklığında P=nRT/VP=nRT/V denklemi ayrı ayrı,

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

Toplam basınçlar farklıdır; zarın sabit tutulması bu durumu korur. A eşit sıcaklık ve miktarlarla eşit hacimleri doldurduğundan kısmi basınçları eşittir.

Hareketli zar. Soru 4'ten,

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.

Toplam, kabın hacmi 20 L20\ \mathrm L olarak kalır. A parçacık sayıları artık NA,1=N0N_{A,1}=N_0 ve 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}.

Toplam gerçekten 3 mol3\ \mathrm{mol} A'dır. B eklenince bölümler sırasıyla 22 ve 4 mol4\ \mathrm{mol} gaz içerir. Yuvarlamadan önce tam hacimleri kullanalım:

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

Serbest zarın mekanik dengesi gereğince iki hesap aynı toplam basıncı verir.