COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用


  • 紐曼模型框架

     紐曼模型(Newman model)是用于描述鋰離子電池內部電化學和傳輸過程的一種數學模型。該模型以電池的正負極為基礎,通過一組偏微分方程來描述電池內部的電流、電壓和鋰離子濃度分布等關鍵參數。這個模型的主要目標是理解電池的性能和響應,以優化電池設計和管理。

以下是紐曼模型中的主要元素和方程:

  • 電極反應:模型考慮了正負極的電化學反應。在正極,鋰離子從電解質中遷移到正極材料,發生氧化反應。在負極,鋰離子從正極材料脫嵌并進入負極材料,發生還原反應。

  • 擴散:模型考慮了鋰離子在電解質中的擴散過程,其中Fick's第一定律用于描述鋰離子濃度梯度對擴散速度的影響。這包括正極和負極內的擴散。

  • 電解質導電性模型考慮了電解質的離子導電性,其中Ohm's Law用于描述電流與電場強度之間的關系。這部分描述了電池中的電流分布。

  • 極化模型還包括了由于電池材料的不完美和非均勻性而導致的極化效應。這些效應包括極化電阻、濃差極化等,會影響電池的性能和響應。

紐曼模型的具體數學方程因電池類型和設計而異,通常需要進行一些假設和簡化來使問題變得可行。紐曼模型中將復雜的電化學行為分為兩個相,液態電解質相和固態電極相,共由五個微分方程組成。這五個微分方程的作用示意圖如下。

在紐曼框架基礎上,可以耦合各種其他物理過程方程來擴展模型的能力(應對紐曼模型描述不了的場景)


COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖1


  • 電熱耦合

     電化學-熱耦合模型是基于電化學反應產熱而建立的電池模型,在紐曼模型的框架上耦合固體傳熱接口,主要用于模擬電池的溫度變化分布情況。鋰離子電池電化學-熱耦合模型由兩部分組成:研究電池內部化學反應的電化學模型以及描述電池溫度分布的熱模型。這兩個部分分工明確并相互耦合。首先,電化學模型計算出發熱功率,然后將發熱功率傳遞給熱模型,熱模型根據發熱功率計算出溫升,然后將此時電池溫度傳遞給電化學模型中受溫度影響的各參數,以此互相耦合實現電池的電壓和溫度模擬。電化學-熱耦合模型涉及的理論方程也分為兩部分,一部分是電化學模型所用 到的電荷守恒、質量守恒以及電極動力學,另一部分是熱模型構建所用的結合生熱、傳熱與散熱的能量守恒關系。兩部分相互耦合,使得模型能夠準確地反映出電池的電化學性能與熱性能,示意圖如下。?


COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖2


  • 電力耦合

     電化學-力耦合模型基于電化學插層反應而建立的電池模型,在紐曼模型的框架上耦合固體力學接口,主要用于模擬電池的內部應力變化分布情況。鋰離子電池電化學-力耦合模型由兩部分組成:研究電池內部化學反應的電化學模型以及描述電池應力分布的固體力學模型。鋰離子電池在電極的嵌入脫出過程稱之為插層反應,可類比于吸濕膨脹。對于吸濕膨脹,當固體進入一個潮濕的環境時,其中一部分固體有可能會吸收水分子。吸收和儲存水分子會導致固體膨脹,并使其面臨更大的應力和應變,這種效應被稱為吸濕膨脹。吸濕膨脹會產生非彈性應變,該應變與濃度和無應變參考濃度之間的差異成正比。在小變形理論中,吸濕膨脹的貢獻是疊加的,也就是說,非彈性應變是其他非彈性應變和吸濕應變的總和。吸濕膨脹系數是一個二階張量,可以被定義為各向同性、對角線或對稱性。因此,在不同的方向上,膨脹可以是不同的。在鋰電插層反應中,將石墨負極中嵌入鋰離子的濃度實時輸入進固體力學場中的應變接口中,即可計算電池充放電過程中,由于鋰離子濃度改變造成的膨脹和收縮效應,示意圖如下。


COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖3


  • 熱力耦合

     熱-力耦合模型基于電化學產熱反應而建立的耦合模型,在固體傳熱框架上耦合固體力學接口,主要用于模擬電池的內部由于溫度變化引起的應力變化。當固體材料的溫度上升時,其結構體積會因此而增加,這種現象稱為熱膨脹。受熱使得材料的動能增加,從而引發這一過程。固體分子通常是緊密排列的,因此固體具有一定的結構形狀。隨著溫度的上升,分子開始以更快的速度振動,并相互推擠。這一過程使相鄰原子間的距離增大,引起固體發生膨脹,進而使固體結構的體積增大。隨著結構體積的增加,固體單元會承受更高水平的應力。熱應力會對固體結構的強度和穩定性產生很大的影響,并可能使某些組件出現裂紋或斷裂。這些故障會破壞結構的整體設計,從而導致潛在的強度減弱和變形。此處注意的是,熱膨脹和插層膨脹的作用機理相似,插層膨脹系數與熱膨脹系數是影響這兩種膨脹的關鍵參數。


COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖4


  • 力相耦合

     熱-力耦合模型基于固體力學框架而建立的耦合模型,在固體力學上耦合相場損傷接口,主要用于模擬電池的內部由于應力變化引起的材料斷裂和失效。相場方法的核心思想是利用彌散的相邊界描述實際上較為尖銳的邊界,通過引入序參量,便可用連續函數描述斷裂模型,并通過相場控制方程控制序參量的演化,使得在模擬時不用顯式地追蹤裂紋面,而是通過序參量的自動演化獲取裂紋路徑及位置。


COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖5


  • 結果分析

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖6

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖7

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖8

    電池整體相場演變

    正極局部裂紋演變

    正極端部斷裂

    電池的斷裂情況收到應力影響,應力受到插層膨脹和熱膨脹兩方面的影響。

    正極的端部受到的不均因應力較大,顏色越深則破碎的概率越高。

    調整損傷因子后,正極端部呈現破碎狀態。相場斷裂指導極片設計。



    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖9


    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖10

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖11

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖12

    電池熱源分布

    電池溫度變化

    電池截面溫度分布

    電池內部卡靠近電解液的地方反應速率更快,產熱速率更快。

    電池隨著充電過程溫度整體上升。具體溫度變化需要進一步設置。

    電池的二位截面圖顯示了一個給定時間的溫度等高線分布。



    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖13


    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖14

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖15

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖16

    Von-Mises應力分布

    彈性應變能密度分布

    平均應力分布

    Von-Mises應力指導電極何處發生屈服斷裂,電極端部數值最高。

    彈性應變能密度是裂紋產生的驅動力,電極端部數值最高。

    平均應力是三個主應力之和,該應力不影響電極形變。



    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖17


    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖18

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖19

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖20

    鋰離子濃度變化

    電解質鹽濃度變化

    電流密度變化

    電池充電過程中,鋰離子從負極脫嵌進入正極。放電過程則相反。

    電池充電過程中電解質中的鋰鹽會參與到插層反應中。

    電流密度和離子通量關系密切,電極內部電流密度不均勻。



    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖21


    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖22

    鋰離子濃度和擴散通量。箭頭顯示離子傳輸方向。

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖23

    電解液濃度和擴散通量。箭頭大小顯示離子傳輸方向。

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖24

    電極電流密度,顯示電池內部的反應速率情況。

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖25

    Trescay應力分布,該應力體現電極的剪切破壞狀態。



    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖26


    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖27

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖28

    電池整體溫度分布

    電池內部鋰離子濃度分布

    電極和電解質產生的熱源使正極的溫度略高于負極的溫度。由于電極顆粒內部的溫度梯度幾乎為零,因此不同顆粒點處的溫度應力幾乎沒有變化,這對變形的影響不大

    通過建立異構鋰離子模型,體現出負極部分顆粒中與顆粒間的鋰離子濃度分布,隨著充電時間的增加,整體濃度上升,濃度梯度有所下降??拷娊庖禾帩舛容^高。

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖29

    COMSOL鋰電池技術仿真與應用(九)鋰電池電-熱-力-相全耦合模型搭建與應用的圖30

    正極顆粒徑向應力分布

    正極切向應力分布

    插入層狀材料中的鋰原子失去電子,并通過電解質向正極擴散。在粒子中,鋰離子濃度形成梯度沿半徑。由于離子的脫插,負極材料中的分子間力減小,間插層之間的空間減小,在宏觀上表現為粒子收縮。

    負粒子的鋰離子脫插過程(反應機理設置與陽極相同)和沿半徑的切向應力值。注意膨脹環向過程中的張應力是地表斷裂的驅動力,而徑向壓應力有利于保持結構穩定性。圖中顯示脫嵌機理。

    文章來源:中礦鋰電

登錄后免費查看全文
立即登錄
App下載
技術鄰APP
工程師必備
  • 項目客服
  • 培訓客服
  • 平臺客服

TOP

5
1
8