電雙層電容孔隙漸變多孔電極電化學特性分析

李澤宇，陳 威

(上海海事大學商船學院，上海 201306)

電化學超級電容器具有充電快、充放電周期長、工作溫度范圍寬等優點。電化學超級電容的數學模型可以分為三種類型：多孔介質微尺度建模[1-3]、等效電路模型[4-6]和多孔電極理論框架[7-9]。由于多孔電極中的變孔隙結構對電解質中離子的遷移和擴散有很大的影響。因此，許多學者將這種電極結構應用到液流電池和鋰離子電池的電極上[10-12]。這種多孔結構化設計使得相應的電化學性能有所提升。

本文基于多孔電極理論框架，將Chen 等[10]提出的新型電極結構布置在多孔電極上，并參考Subramanian 等[13]的研究修正了電解質和固相多孔電極的電導率；考慮到電解質溶液的擴散和遷移來分析多孔電極上的變孔隙結構對電雙層電容性能的影響，進一步補充了分級多孔電極的理論研究。通過在電雙層電容的電極上采用漸變孔隙結構，模擬分析了其結構對恒流恒壓及恒功率下的電化學性能的影響，得到了不同孔隙結構下的恒流充電恒壓放電特性，恒功率充電下電雙層電流源和電荷密度及恒功率放電下的電流電壓特性。

1 物理模型

如圖1(a)所示，電雙層電容器是由兩個活性炭電極、電解質溶液、隔膜等組成，隔膜將兩個多孔電極隔開，非水電解質LiPF6溶液充滿整個電雙層電容器。在電雙層中，由于電解液中的陰陽離子可以被穩定地吸附到具有活性材料的多孔電極上，因此，電荷是靜電存儲的。雙電層電容的電荷存儲是發生在電極材料上，在界面上沒有發生任何的電荷轉移。由于在多孔活性炭電極表面上所吸附的帶有正負電荷的離子被平衡，從而產生了電容。

圖1(b)為電雙層電容器幾何結構示意圖，藍色區域為具有多孔性的隔膜，里面充滿著自由電解質。左右黑色區域是含活性炭的多孔電極，電極孔隙內的溶液是電中性的。沿x軸依次排列的是50 μm 的活性炭電極，25 μm 的隔膜，50 μm的右側活性炭電極。在充電過程中電解質中的離子向兩側多孔電極擴散和遷移，鋰離子在多孔碳電極表面積累，形成雙電層；在放電過程中，雙電層被中和，其中儲存的能量也隨之釋放出來，從而達到放電的效果。

圖1 電雙層電容器原理示意圖(a)與幾何結構示意圖(b)

2 數學模型

2.1 模型假設

在模擬分析電雙層電容器的電化學特性時，做如下假設：

(1)該模型模擬的是一維的超級電容截面，忽略其長度和高度的邊緣效應。

(2)假設正負離子在dq+/dq=dq-/dq=-1/2 上是相等的。

(3)假定電容aC是恒定的。

(4)擴散過程的阻力無作用，忽略擴散阻力。

(5)認為是一維等溫的，忽略溫度的變化。

(6)當電極厚度為50 μm 時，假定多孔電極理論是有效的，因為輸運現象基本上是一維的。

2.2 描述方程和邊界條件

有效介質理論考慮了孔隙度和彎曲度對輸運參數值的影響。對液相采用了兩個唯象方程。首先，離子通過遷移和擴散進行輸運，一維通量關系由能斯特-普朗克方程提供：

式中：Nj為物質通量；Dj為擴散系數，m2/s；zj為電荷數；F為法拉第常數，C/mol；R為通用氣體常數，J/(mol·K)；T為絕對溫度，K；cj為濃度，mol/L；Φl為電解質液相電勢，V。

采用的第二種關系由于雙電層充放電而使物質j遠離電極表面的摩爾通量表示為：

式中：j為稀物質種類；C為電容，F；qj為電極表面電荷；Φs為固相多孔電極電勢，V。

固相多孔電極采用歐姆定律：

式中：is為固相多孔電極側的電流，A；σ為固相電導率，S/m。

考慮孔隙率和曲度，定義了多孔電極的有效電導率和有效擴散系數：

式中：ε為孔隙率；τ為曲度。

液相的歐姆定律呈濃度有關的非線性形式：

式中：il為液相電解質側的電流，A；k為電解質電導率，S/m；t+為電荷轉移數。

其中轉移數tj與離子擴散系數有關：

電解質的有效電導率與阿奇指數[14]有關：

式中：keff為電解質有效電導率，S/m。

通過能斯特方程計算遷移離子的遷移率：

中性液相中電荷守恒為：

式中：a為比表面積，cm2/cm3。

二元電解質的物質平衡產生的守恒方程：

其中多孔電極處的比表面積與孔隙率可以表示為[15]：

邊界條件如下所示。

在左集流體處（x=0）有：

在電極-隔膜邊界處（x=LA和x=LA+LS）均有：

在右集流體處（x=LA+LS+LB）有：

計算中所用參數和負載參數分別如表1 和表2 所示。

表1 計算中所用的參數

表2 負載參數

2.3 模擬計算過程和模型驗證

對所建立的模型使用COMSOL Multiphysics5.5 進行求解，其中在電化學模塊選擇“三次電流分布，Nernst-Plank”，以及使用數學中的事件來模擬由恒電流充電-恒電壓放電組成的負載循環。使用全局方程節點計算電流密度，在右側活性炭多孔電極邊界上指定恒功率充電，用靜息電位來確定初始值，采用參數掃描法模擬了不同的外加功率和靜息電位。選擇向后差分公式(BDF)方法用于時間步進器。求解使用的迭代方法是廣義最小殘差法，相對容差設定為10-4。

網格質量對數值計算的穩定性和準確性尤為重要。本文對模型網格進行了用戶控制網格，通過在兩側的多孔電極創建更高分辨率的網格，使計算結果準確穩定。自定義使用三種不同的網格數：普通網格(424 單元)、細網格(536 單元)、較細網格(1 948 單元)。分析結果如表3 所示。

表3 恒電流充電電壓范圍隨網格數變化

結果表明，網格大小質量對模擬結果產生的影響極小，為了節省計算所需的時間，本文采取較細化網格。將25 ℃下恒電流充電和電壓實驗所得到的電流電壓曲線，與數值模擬計算所得到的一個周期內25 ℃下的恒電流電壓曲線進行對比。實驗數據來自文獻[18]。從圖2 可以看出，模擬曲線與實驗曲線吻合度較高，由于模型建立過程中存在很多假設理想化以及部分物性參數選取有所差異，這個誤差是可以允許的，因此驗證了模擬的正確性。

圖2 電雙層電容恒電流電壓模擬結果與實驗[18]對比(25 ℃，100 A恒流充電，1.4 V恒壓放電)

3 結果與討論

3.1 漸變孔隙率結構對恒流充電恒壓放電性能的影響

多孔電極中不同孔隙結構下對應的電雙層電容器的恒流充電恒壓放電特性如圖3 所示。如圖3(a)所示，在孔隙率ε=0.67 時，電容器以100 A 恒流充電至1.8 V，然后恒壓5 s，靜息180 s，分別充電至最大電壓為1.8、2.0、2.2 和2.4 V。每個周期的電壓依次升高0.2 V，最終經過4 個周期充電電壓達到2.4 V。實驗中在很短時間內捕獲大電流充放電是較為困難的，一般恒壓充電2.2 V 放電的最大電流實際上大于2.4 V 放電時的最大放電電流。因此，通過數值模擬可以在短時間內得到準確的充放電特性。

通過采用不同孔隙結構可以提高充電電壓及充放電速率。通過圖3(a)和(b)的比較，可以發現采用較大孔隙率的多孔電極，電雙層電容的最大充電電壓有所提升。因為在較大孔隙的電極中，電解液的有效擴散系數較大，擴散阻力減小，而且遷移離子的遷移率也隨之增大。從圖3(a)和(c)可以看出，在多孔電極上使用漸變增加孔隙率結構可以提高最大充電電壓和充放電速率。從圖3(b)和(c)可以看出，對于較大孔隙的電極結構，采用漸變增加孔隙的電極結構一定程度上提高了系統的充放電速率。因此，相對于單純的較大孔隙率而言，漸變增加孔隙的多孔結構更具有電化學優勢。這是因為在漸變增加孔隙率的多孔電極中，電解質從電極內側可以較好地擴散至電極外側，內部電解質分布較為均勻。而且在線性增加孔隙率中，電極內側的孔隙率較大，使得此處具有較好的擴散率及較小的擴散阻力。然而，從圖3(c)和(d)中可以看出，電雙層電容充放電速率有所下降，因為漸變減小孔隙率的電極內側孔隙率較小，所以比表面積較大，從電雙層中釋放電解質較多，在電極內側形成電解質的累積。此時有效擴散系數減小，擴散阻力隨之增大。充放電速率有所下降。

圖3 不同孔隙結構下的恒流充電-恒壓放電電壓電流分布

3.2 漸變孔隙率結構對不同周期恒流充電恒壓放電響應的影響

電雙層電容重疊的電流電壓分布特性如圖4 所示，從圖中可以看出在恒流充電中每個周期增加的最大電壓值。從圖4(a)中可以看出，當孔隙率ε=0.67 時，電雙層電容處于恒流100 A 充電，恒壓1.4 V 放電的工況。當電壓達到2.5 V 的充電上限目標電壓時，一共要經歷3 個周期。在2.0～2.5 V 之間的電流與電壓曲線重疊，此時為電雙層電容的恒功率運行狀態。圖4(a)與(b)相比，可以看出在較大孔隙率的電雙層電容進行充電時到達各個周期的電壓所需時間相對縮短，提升了充電速率。從圖4(b)和(c)可以發現，采用漸變增加孔隙的電極結構，充電速率提升最為明顯，縮短了充電時間。此外，從圖4(d)中可以發現，對應的漸變減小孔隙結構的充放電周期加長，不利于系統性能的提升。由于漸變增加孔隙結構合理地規劃了具有較大比表面積的小孔隙率電極和具有較高擴散率和遷移率的大孔隙率電極，因此，漸變增加孔隙結構更具有電化學優勢，進一步提升了電雙層電容器的功率密度。

圖4 不同孔隙結構下的恒流充電和恒壓放電周期響應

3.3 漸變孔隙率結構對恒功率充電下電雙層電流源的影響

在恒功率充電下的電雙層電流源和電荷密度可以作為超級電容的電極利用率的指標。電化學電容在無量綱長度上的恒功率充電時的電雙層電流源分布特性如圖5 所示。從圖5(a)可以看出，當孔隙率ε=0.67且系統功率為0.5 kW、靜息電位為-2.1 V 時，電雙層電流源差值為2.5×106A/m3；隨著功率增加到0.85 kW 且靜息電位為-2.5 V 時，電雙層電流源差值為3.5×106A/m3。從圖5(b)中可以發現，對于較大孔隙率ε=0.77 且功率為0.5 kW 時，其電雙層電流源差為2.4×106A/m3；與此同時，當功率為0.85 kW 時，其電雙層電流源差為3.5×106A/m3。從圖5(c)中可以看出，對于漸變增加孔隙率ε=0.67～0.77 且功率為0.5 kW 時，其電雙層電流源差為2.6×106A/m3；隨著功率增加到0.85 kW，電雙層電流源差為4.1×106A/m3。在圖5(d) 中可以發現，漸變減小孔隙率ε=0.77～0.67且功率為0.5 kW，可以獲得2×106A/m3的電雙層電流源差；此外，功率增加到0.85 kW，與之對應的電雙層電流源差為3×106A/m3。

圖5 不同孔隙結構下恒功率充電下的雙電層電流源

在漸變增加孔隙率中隨著電極外側孔隙率的減小，漸變電極的平均孔隙率減小，電極的比表面積增大，引起電荷在電極表面積累的空間變大，電解液中的離子聚集增多，強化了電雙層電荷儲存量。因此，通過改變孔隙結構可以對電雙層電流源有所改進，其中漸變增加孔隙排列方式在較小功率0.5 kW 和較大功率0.85 kW，電雙層電流源都有較大的提升。

3.4 漸變孔隙率結構對恒功率充電下電荷密度的影響

恒功率充電時電雙層電容器的電荷密度分布如圖6 所示。從圖6(a)可以看出，當孔隙率ε=0.67 且功率為0.5 kW、靜息電位為-2.1 V 時，電荷密度變化為4.5×106C/m3；隨著功率增加到0.85 kW 且靜息電位為-2.5 V 時，對應的電荷密度變化為5×106C/m3。對于較大孔隙率ε=0.77 且功率0.5 kW 時，電荷密度變化增加了5×106C/m3；隨著功率增加到0.85 kW，對應的電荷密度變化增大了7×106C/m3。當漸變增加孔隙率ε=0.67～0.77 并且功率為0.5 kW，相應的電荷密度變化為4×106C/m3；隨著功率增加到0.85 kW，與之對應的電荷密度變化為5×106C/m3。當漸變減小孔隙率ε=0.77～0.67 并且功率為0.5 kW，相應的電荷密度變化為4.5×106C/m3；此外，當功率增大到0.85 kW 時，其電荷密度變化為6×106C/m3。

圖6 不同孔隙結構下的恒功率充電下電荷密度分布

從圖6(a)和(c)中可以發現，漸變增加孔隙率ε=0.67～0.77的電極在同一低功率下的電荷密度變化比常孔隙率ε=0.67的小，隨著功率的增加，在較高功率下的電荷密度與ε=0.67的變化相同。圖6(b)和(c)相比可以發現，采用漸變增加孔隙率電極結構的電容器，無論是高或低功率下，電荷密度變化幅度都有所減小。一般而言，電荷分布的均勻性隨充電功率的增大而減小，充電功率較低時對應的電極更均勻，電極利用率較高。漸變增加孔隙率可以在功率增大的同時保證電荷在電極上的均勻性。這是因為電極的比表面積增大，電解質溶液中的電荷更能均勻地分布在電極上。因此，與?？紫堵孰姌O結構相比，采用漸變增加孔隙結構的多孔電極可以提高電極利用率。

3.5 漸變孔隙率結構對恒功率放電的電流電壓的影響

圖7 顯示了不同功率參數下的放電電流和電壓分布。當常孔隙率ε=0.67 且功率為0.5 kW 時，對應的電壓為2.35 V；隨著功率增加到0.85 kW，相應的電壓為2.8 V。當?？紫堵师?0.77 且功率為0.5 kW 時，獲得的電壓為2.4 V；隨著功率增加至0.85 kW，與之對應的電壓為2.85 V。當漸變增加孔隙率ε=0.67～0.77 且功率為0.5 kW 時，相應的電壓為2.6 V；隨著功率增大到0.85 kW，與其對應的電壓為3.1 V。當漸變減小孔隙率ε=0.77～0.67 且功率為0.5 kW 時，對應的電壓為2.3 V；隨著功率增加到0.85 kW，所對應的電壓為2.75 V。在常孔隙率的電雙層電容中表現出充電功率越高，電壓變化越大。對比多孔電極上不同孔隙結構發現，采用漸變增加孔隙率在低功率和高功率下都可以較好地提高電雙層電容的放電電壓。這是因為漸變增加孔隙率可以提高固相電極的有效電導率。而且與常孔隙率相比，使用這種電極結構由于電極內孔隙率的合理布置，使得電雙層電容器在放電過程中的電壓損失減小。

圖7 不同孔隙結構下恒功率放電的電流電壓分布

4 結論

本文主要研究了一種漸變孔隙分布電極結構的電雙層電容器，研究結果表明，漸變增加孔隙率對其性能有較好的改進。得到的主要結論如下：

(1)在電雙層電容的恒流充電恒壓放電中，使用漸變增加孔隙結構相比于ε=0.67 來說，在第四個周期達到最大充電電壓時，所耗費時間從680 s 縮短到620 s，減小了60 s，因此充放電速率有所提升。

(2)在不同周期恒流充電恒壓放電響應中的恒功率運行下，漸變增加孔隙結構相比于ε=0.67 而言，恒功率運行狀態從15～27 s 縮短到6～12.5 s，達到恒功率時所需時間縮短了15 s，一定程度上提高了系統的功率密度。

(3)針對恒功率充電下雙電層電流源而言，與ε=0.67 相比，采用漸變增加孔隙率的電極結構在功率均為0.5 kW 時，電雙層電流源增加了0.1×106A/m3；當功率密度均為0.85 kW時，電雙層電流源增加了0.6×106A/m3。這種漸變電極結構強化了電雙層電荷儲存量。

(4)在恒功率充電電荷密度變化下，相對于ε=0.67 的電極結構來說，漸變增加孔隙率在功率均為0.5 kW 時，電荷密度變化減小了0.5×106C/m3；當功率增大到0.85 kW 時，電荷密度保持不變。漸變增加孔隙結構在小功率時減小了電荷密度的變化，在功率增加的同時保證了電荷在電極上的均勻性，一定程度上提高了超級電容的電極利用率。

(5)在恒功率放電下的電流電壓中，相比ε=0.67 的電極結構，采用漸變增加孔隙率在功率均為0.5 kW 時，放電電壓提升了0.25 V；在功率為0.85 kW 時，放電電壓提升了0.3 V。這種電極結構減小了放電過程中的電壓損失，適當地提高了超級電容的放電電壓。