تمرين محلول في الديناميكا الحرارية

غشاء شبه نفوذ: ما الكميات التي تتساوى؟

تمرين 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 يفرض ذلك تساوي الضغطين الكليين. درجة الحرارة النهائية مجددا T0T_0، من حصيلة الطاقة نفسها في السؤال 3. بما أن 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}.

يعطي الحسابان الضغط الكلي نفسه، كما يقتضي التوازن الميكانيكي للغشاء المحرر.