1055 words
5 minutes
From Fugacity to Raoult's Law
2024-11-29
TIP
  • 建議先看這個寶藏漫畫
  • Fugacity (逸壓):簡單想就是修正過的壓力 (非理想狀態,考慮分子間作用)

Why it matters?#

  1. 觀念適用於任何相:

    Fugacity=concentration×coefficient×pressure/idealfugacityFugacity = concentration \times coefficient \times pressure/ideal fugacity

    Where the coefficient describes the ideality of the species.

  2. 分析VLE時,若達平衡狀態:

    • 只要知道其中一相的濃度,就可以直接推出另一相的濃度
    • 計算各相含多少特定成分物質
  3. 常見應用 (=常考):分餾 (distilation)

    • 了解物質的混和、分離情形
    • 計算出料濃度 -> 反應器效率

Fugacity of Gases#

課本在CH7、9分別介紹了純物種和溶液中的fugacity,但總覺得不甚直覺 (可能是個人問題),所以想了一個我比較能理解的推導方式:

  1. 回想一下,我們是這樣定義化學位能 (chemical potential)

    • 純物種

      μ(T,P)=μ(T)+RTln(P)\mu(T,P) = \mu^\circ(T) + RT \ln(P)

    • 溶液 (物種 iiμi\mu_i)

      μi(T,P)=μi(T)+RTln(P)\mu_i(T,P) = \mu_i^\circ(T) + RT \ln(P)

    其中 μ(T)\mu^\circ(T)μi(T)\mu_i^\circ(T) 是標準狀態 (standard state, often P = 1 bar) 下的化學位能。

  2. 接著,在理想狀態下

    • 純物種

      GˉiIG(T,P)=μiIG(T,P)=μi(T)+RTln(P)\bar{G}_i^{IG}(T,P) = \mu_i^{IG}(T,P) = \mu_i^\circ(T) + RT \ln(P)

    • 溶液 (要考慮分壓:PiIGM=xiPP_i^{IGM} = x_i P)

      GˉiIGM(T,P,x)=μiIGM(T,P,x)=μi(T)+RTln(xiP)\bar{G}_i^{IGM}(T,P,x) = \mu_i^{IGM}(T,P,x) = \mu_i^\circ(T) + RT \ln(x_i P)

  3. 然而,現實世界中非理想的狀態下

    • 純物種

      Gˉi(T,P)GˉiIG(T,P)\bar{G}_i(T,P) \neq \bar{G}_i^{IG}(T,P)

    • 溶液

      Gˉi(T,P,x)GˉiIGM(T,P,x)\bar{G}_i(T,P,x) \neq \bar{G}_i^{IGM}(T,P,x)

  4. 於是定義 fugacity fˉi\bar{f}_i,作為修正過的壓力,取代原先假設理想狀態的壓力

    • 純物種

      GˉiIG(T,P)=μi(T)+RTln(P)\bar{G}_i^{IG}(T,P) = \mu_i^\circ(T) + RT \ln(P)

      Gˉi(T,P)=μi(T)+RTln(fˉi)\bar{G}_i(T,P) = \mu_i^\circ(T) + RT \ln(\bar{f}_i)

    • 溶液

      GˉiIGM(T,P,x)=μi(T)+RTln(xiP)\bar{G}_i^{IGM}(T,P,x) = \mu_i^\circ(T) + RT \ln(x_i P)

      Gˉi(T,P,x)=μi(T)+RTln(fˉi)\bar{G}_i(T,P,x) = \mu_i^\circ(T) + RT \ln(\bar{f}_i)

  5. 將修正前後的式子相減

    • 純物種

      GˉiGˉiIG=RT[ln(fˉi)ln(P)]\bar{G}_i - \bar{G}_i^{IG} = RT \left[ \ln(\bar{f}_i) - \ln(P) \right]

      GˉiGˉiIG=RTln(fˉiP)\bar{G}_i - \bar{G}_i^{IG} = RT \ln\left(\frac{\bar{f}_i}{P}\right)

      exp(GˉiGˉiIGRT)=fˉiP\exp\left(\frac{\bar{G}_i - \bar{G}_i^{IG}}{RT}\right) = \frac{\bar{f}_i}{P}

      fˉi=Pexp(GˉiGˉiIGRT)\bar{f}_i = P \exp\left(\frac{\bar{G}_i - \bar{G}_i^{IG}}{RT}\right)

    • 溶液

      GˉiGˉiIGM=RT[ln(fˉi)ln(xiP)]\bar{G}_i - \bar{G}_i^{IGM} = RT \left[ \ln(\bar{f}_i) - \ln(x_i P) \right]

      GˉiGˉiIGM=RTln(fˉixiP)\bar{G}_i - \bar{G}_i^{IGM} = RT \ln\left(\frac{\bar{f}_i}{x_i P}\right)

      exp(GˉiGˉiIGMRT)=fˉixiP\exp\left(\frac{\bar{G}_i - \bar{G}_i^{IGM}}{RT}\right) = \frac{\bar{f}_i}{x_i P}

      fˉi=xiPexp(GˉiGˉiIGMRT)\bar{f}_i = x_i P \exp\left(\frac{\bar{G}_i - \bar{G}_i^{IGM}}{RT}\right)

  6. 為了方便計算,試著只用 V,T,PV, T, P 表示 (假設 const. TT, VV: molar volume)

    dG=SdT+VdPdG = -S\,dT + V\,dP

    (GP)T=V\left( \frac{\partial G}{\partial P} \right)_T = V

    • 純物種

      G(T,P2)G(T,P1)=P1P2VdPG(T,P_2) - G(T,P_1) = \int_{P_1}^{P_2} V \, dP

      GIG(T,P2)GIG(T,P1)=P1P2RTPdP=RTlnP2P1G^{IG}(T,P_2) - G^{IG}(T,P_1) = \int_{P_1}^{P_2} \frac{RT}{P} \, dP = RT \ln\frac{P_2}{P_1}

      [G(T,P2)G(T,P1)][GIG(T,P2)GIG(T,P1)]=P1P2(VRTP)dP[G(T,P_2) - G(T,P_1)] - [G^{IG}(T,P_2) - G^{IG}(T,P_1)] = \int_{P_1}^{P_2} \left( V - \frac{RT}{P} \right) dP

    • 溶液

      Gi(T,P2,x)Gi(T,P1,x)=P1P2VidPG_i(T,P_2, x) - G_i(T,P_1, x) = \int_{P_1}^{P_2} V_i \, dP

      GiIG(T,P2,x)GiIG(T,P1,x)=P1P2ViIGdPG^{IG}_i(T,P_2, x) - G^{IG}_i(T,P_1, x) = \int_{P_1}^{P_2} V^{IG}_i \, dP

      [Gi(T,P2,x)Gi(T,P1,x)][GiIG(T,P2,x)GiIG(T,P1,x)]=P1P2(ViViIG)dP[G_i(T,P_2, x) - G_i(T,P_1, x)] - [G^{IG}_i(T,P_2, x) - G^{IG}_i(T,P_1, x)] = \int_{P_1}^{P_2} (V_i - V^{IG}_i) dP

  7. 假設 P = 0,此時任何流體都是理想狀態,因此 G(T,0)=GIG(T,0)G(T,0) = G^{IG}(T,0)

    • 純物種

      G(T,P)GIG(T,P)=0P(VRTP)dPG(T,P) - G^{IG}(T,P) = \int_0^{P} \left( V - \frac{RT}{P} \right) dP

    • 溶液

      Gi(T,P,x)GiIG(T,P,x)=P1P2(ViViIG)dPG_i(T,P, x) - G^{IG}_i(T,P, x) = \int_{P_1}^{P_2} (V_i - V^{IG}_i) dP

  8. 由以上推導可得結論 (我們通常只會被告知結論),對任何氣體

    • 純物種 (One Component)

      f=Pexp ⁣(GGIGRT)=Pexp ⁣(1RT0P(VRTP)dP)f = P \exp\!\left(\frac{G - G^{IG}}{RT}\right) = P \exp\!\left(\frac{1}{RT} \int_0^{P} \bigl(V - \frac{RT}{P}\bigr)\, dP \right)

      ϕ=fP=exp ⁣(GGIGRT)=exp ⁣(1RT0P(VRTP)dP)\phi = \frac{f}{P} = \exp\!\left(\frac{G - G^{IG}}{RT}\right) = \exp\!\left(\frac{1}{RT} \int_0^{P} \bigl(V - \frac{RT}{P}\bigr)\, dP \right)

    • 溶液 (Mixtures)

      fˉi=xiPexp ⁣(GˉiGˉiIGMRT)=xiPexp ⁣(1RT0P(VˉiVˉiIG)dP)\bar{f}_i = x_i P \exp\!\left(\frac{\bar{G}_i - \bar{G}_i^{IGM}}{RT}\right) = x_i P \exp\!\left(\frac{1}{RT}\int_0^P \bigl(\bar{V}_i - \bar{V}_i^{IG}\bigr)\, dP\right)

      ϕˉi=fˉixiP=xiPexp ⁣(GˉiGˉiIGMRT)=xiPexp ⁣(1RT0P(VˉiVˉiIG)dP)\bar{\phi}_i = \frac{\bar{f}_i}{x_i P} = x_i P \exp\!\left(\frac{\bar{G}_i - \bar{G}_i^{IGM}}{RT}\right) = x_i P \exp\!\left(\frac{1}{RT}\int_0^P \bigl(\bar{V}_i - \bar{V}_i^{IG}\bigr)\, dP\right)

    Where ϕ\phi, ϕˉi\bar{\phi}_i are the fugacity coefficients.

IMPORTANT

為了區別,通常會以 yiy_i 表示氣體的莫耳分率,液體則用 xix_i

Deriving Raoult’s Law#

  • 系統:容器中,液態溶液與其汽相達平衡 (皆假設理想狀態)
  • 主要應用:計算蒸氣壓

Fugacity of Liquids#

在前一段我們得到了 gas mixture 的 fugacity:

fˉiV=yiϕˉivP\bar{f}_i^V = y_i \bar{\phi}_i^v P

與氣體不同的是,液體、固體沒有所謂的分壓 1,也就是說壓力與fugacity不能直接寫成簡單的關係式。因此,以液體為例,我們必須多假設其在理想狀態下有一 fugacity filf_i^l,並套用同樣的觀念表示:

fˉil=xirifil\bar{f}_i^l = x_i r_i f_i^l

Where rir_i is the activity coefficient, describing the ideality of fluids.

NOTE

液體也可以用同樣的方式推導:

Gˉi=μi(T)+RTln(fˉi)\bar{G}_i = \mu_i^\circ(T) + RT \ln(\bar{f}_i)

若為理想溶液:

Gˉiid=μi(T)+RTln(xifi)\bar{G}_i^{id} = \mu_i^\circ(T) + RT \ln(x_i f_i)

GˉiGˉiid=RTln(fˉixifi)\bar{G}_i - \bar{G}_i^{id} = RT \ln(\frac{\bar{f}_i}{x_i f_i})

我們於是定義:

ri=fˉixifir_i = \frac{\bar{f}_i}{x_i f_i}

GˉiE=RTlnri\bar{G}_i^E = RT\ln{r_i} (Excess Gibbs Energy)

VLE (Vapor-Liquid Equilibrium)#

  1. 在溫度 TT 和蒸氣壓 PsatP^{sat} 時,GivGil=0G_i^v - G_i^l = 0 (飽和液體 -> 飽和蒸汽)

    fiv=fil=fisatf_i^v = f_i^l = f_i^{sat}

  2. 而在一飽和溶液中

    • 飽和蒸汽:fˉiV=yiϕˉivP\bar{f}_i^V = y_i \bar{\phi}_i^v P

    • 飽和液體:fˉil=xirilfil\bar{f}_i^l = x_i r_i^l f_i^l

  3. 於是我們可以列出

    yiϕˉivP=xirilfily_i \bar{\phi}_i^v P = x_i r_i^l f_i^l

    而在中低壓fil=fisat=ϕisatPisatf_i^l = f_i^{sat} = \phi_i^{sat} P_i^{sat} 2,因此:

    yiϕˉivP=xirilϕisatPisaty_i \bar{\phi}_i^v P = x_i r_i^l {\phi}_i^{sat} P_i^{sat}

  4. 假設汽相為理想氣體 ϕˉiv=1\bar{\phi}_i^v = 1ϕisat=1{\phi}_i^{sat} = 1;液體為理想溶液 ril=1r_i^l = 1,則:

    yiP=xiPisaty_i P = x_i P_i^{sat}

    此即著名的 拉午耳定律 (Raoult's Law)

Footnotes#

  1. Partial molar properties (待補)

  2. Why? (待補)

From Fugacity to Raoult's Law
https://princesswinnie1122.github.io/posts/2_fugacity/
Author
Winnie Sheng
Published at
2024-11-29