基于動態參數響應模型的動力鋰離子電池循環容量衰減研究?

蔣躍輝艾亮賈明?程昀杜雙龍李書國

1)(中南大學冶金與環境學院,長沙 410083)

2)(湖南艾華集團股份有限公司,益陽 413002)

基于動態參數響應模型的動力鋰離子電池循環容量衰減研究?

蔣躍輝1)艾亮2)賈明1)2)?程昀1)杜雙龍1)李書國1)

1)(中南大學冶金與環境學院,長沙 410083)

2)(湖南艾華集團股份有限公司,益陽 413002)

(2016年12月10日收到;2017年3月10日收到修改稿)

結合鋰離子電池容量衰減主要機制,建立了基于動態參數響應的電化學熱耦合模型,研究磷酸鐵鋰電池循環壽命及電池內部的電化學行為.模型采用恒流恒壓充放電制度對電池進行循環充放電仿真計算.結果顯示:800次循環后,電池容量衰減率約為6.35%,電池內部固體電解質界面膜阻抗增大了15.6 m?·m?2.分別探討了充放電倍率、負極活性物質顆粒粒徑、負極固相體積分數對電池循環壽命的影響.研究表明:400次循環后,相較于1C(C表示充放電倍率)倍率下的容量衰減率3.31%,2C,3C,4C容量衰減率分別達到3.93%,4.69%和5.04%;負極活性顆粒粒徑為2和10μm時對應容量衰減率分別為2.89%,3.87%,差值接近1%;固相體積分數在[0.5,0.6]區間內能保持最好的電池循環穩定性和壽命發揮.

容量衰減,充放電倍率,體積分數,顆粒半徑

1 引言

近年來,隨著電動車的發展,鋰電池的發展制高點已逐漸從消費類電子產品領域轉向動力電池領域(純電動汽車、混合動力汽車),動力電池對性能和循環壽命要求也上升到了一個新的高度[1,2].以電池機理及大量實驗數據為構建基礎的數值模型,可以高效地研究電池循環容量的衰減情況.

根據已公開的研究工作,容量衰減的原因有電解液的氧化還原、電極在循環過程中的相變、負極表面固體電解質界面(SEI)膜的持續生成等[3?9],其中SEI膜的持續生成是電池在較高倍率下循環的主要衰減機制[10].負極表面溶劑與鋰離子發生的氧化還原反應是SEI膜增長的主要原因,且此反應與鋰電池在負載下鋰離子的嵌入/脫出同時發生.Ramadass等[11]第一次提出基于多物理場的鋰離子電池壽命模型,并指出溶劑在負極活性炭上發生的還原反應會不斷發生,從而SEI膜會持續增長,活性鋰也會被不斷地消耗.Phoehn等[12]和Safari等[13]發展了溶劑擴散模型來描述SEI膜在碳負極上的形成和增長.此外,大量研究集中在不同工況條件對鋰電池循環壽命的影響及如何準確地預測鋰電池的服役壽命[14?16].鋰離子電池在工作過程中伴隨著明顯的熱效應,電池溫度的實時變化對鋰離子電池內部化學反應的相關參數及物料傳輸過程有著重要影響.Shi等[17]詳細綜述了從微觀尺度到介觀尺度的相場模型,再到宏觀尺度上的多尺度仿真技術在鋰電池領域應用的國內外研究進展,在研究電池內部的電化學行為和熱行為方面,仿真技術具有傳統研究手段無法比擬的優勢.然而,在眾多針對性研究容量衰減的電化學模型,通常只是單純考慮電池的生熱情況及環境溫度對電池循環壽命的影響,并未將溫度的實時變化引起的參數變化耦合到電化學模型中.Ye等[18]提出了一個考慮電池充放電過程熱效應的容量衰減模型,使用經驗公式描述活性鋰損失及其對循環性能的影響,然而卻沒有從機理上對容量衰減的現象進行分析.Baek等[19]介紹了一種基于單顆粒模型的鋰離子電池模型.該模型仿真了由正極活性材料的損失及負極上SEI膜的形成導致的容量衰減及阻抗行為.采用單顆粒模型計算效率高,但該模型并沒有耦合熱模型,因而忽略了溫度依賴的動態參數變化,模型精確性有待提高.

基于鋰離子電池的基本理論和負極充電過程生成SEI膜的機理,并耦合鋰離子電池產熱模型,本文建立了具有電化學熱耦合特征的鋰離子電池容量衰減模型.模型涵蓋了鋰離子電池實現能量轉化過程中涉及的電場、熱場、濃度場等物理場,基于多物理場的協同作用,對電池循環過程進行仿真研究.另外,通過引入與溫度相關的動態參數,修正了溫度變化對電池行為的影響.利用所建立的模型,對磷酸鐵鋰電池進行了800次循環的仿真計算,詳細分析循環過程中鋰離子電池的電化學行為.此外,分別探討了充放電倍率、負極固相體積分數、活性物質顆粒粒徑等對電池循環壽命的影響.

2 模型建立

考慮到電池的溫度與電化學過程是相互影響的關系,而溫度對電池得很多參數諸如液相擴散系數,電導率等有重要影響,進而影響電池性能.因此,本工作基于Newman和Tiedemann[20]提出的電化學模型,并根據傳熱學基本原理,結合負極充電過程生成SEI膜的機理,建立電化學-熱耦合容量衰減模型.模型中熱耦合的實現主要通過將計算得到的電芯平均溫度實時引入到電化學模型中,并作為電化學反應的溫度.同時,由電化學模型計算得到的平均生熱率會反饋到熱模型,作為電芯的平均生熱率.

2.1 電化學模型建立過程

本文中磷酸鐵鋰電池的電化學模型是以Doyle和Newman[21]準二維模型為基礎,考慮了多孔電極理論、歐姆定律、濃溶液理論、鋰離子的固液相擴散過程及電極動力學.磷酸鐵鋰電池的工作原理如圖1所示.

圖1 (網刊彩色)磷酸鐵鋰電池的工作原理圖Fig.1.(color online)The schematic diagram of the lithium iron phosphate battery model.

2.1.1 電化學反應動力學

電池充放電過程中,鋰離子在以電解液為介質的正負極之間進行傳輸,實現電化學儲能.在電極顆粒表面的脫嵌過程通常用Butler-Volmer方程來進行描述:

式中:jLi為反應電流密度,i0為交換電流密度,η為過電勢;αa,αc分別為陽極、陰極電極反應轉化系數;R,T,F分別為通用氣體常數、電池實際溫度、法拉第常數.其中i0,η表達式如下:

式中,kj為液相電導率;分別為陽極電遷移數、陰極電遷移數;分別為最大固相鋰離子濃度、固相鋰離子濃度和液相鋰離子濃度.

2.1.2 質量守恒

鋰離子在電池內部鋰的物料傳輸主要體現在兩方面,分為電解液中的物料守恒和電極活性顆粒的物料守恒.前者主要描述鋰離子在電解質中的擴散和遷移,忽略對流傳質的影響;后者則采用Fick定律描述鋰離子在電極活性顆粒內的遷移過程.電解質和活性顆粒對應的相態分別為液相和固相,因此質量守恒方程分為固相和液相兩方面.

1)固相區域物質守恒

鋰離子在電極活性材料顆粒中脫嵌過程的質量守恒,用Fick第二定律在球形坐標下描述成

式中,cs,Ds分別表示鋰離子濃度和固相擴散系數;t,r代表時間和距球形顆粒中心的距離.在活性顆粒中存在著兩個邊界:顆粒中心鋰通量為零,顆粒表面鋰通量依據法拉第定律描述,因而邊界條件表示如下:

2)液相區域的質量守恒

式中x表示電池上位置到負極集流體靠近負極活性材料端的距離;Ln,Lsep,Lp分別表示負極活性涂層、隔膜、正極活性涂層的厚度.

2.1.3 電荷守恒

與質量守恒類似,建立的電荷守恒方程分為固相、液相兩部分.

1)固相區域的電荷守恒

其中εs表示電極活性材料的體積分數,rs為活性顆粒粒徑.在固相區域電荷守恒中存在著四個邊界:負極集流體接地,正極集流體上的電流密度則為正極極片平均電流密度,隔膜由于具有電子絕緣特性,因而沒有電流通過.上述四個邊界分別可用(13)—(16)式表示:

式中Ln,cc,Ln,Lsep,Lp,Lp,cc分別表示負極集流體、負極活性涂層、隔膜、正極活性涂層、正極集流體的厚度;φcc,σcc,keffs,Iapp分別代表集流體電勢、集流體電導率、固相有效擴散系數和電極平均電流密度.

2)液相區域電荷守恒

在電解液中主要的電荷傳遞是通過鋰離子的擴散和遷移實現的,依據Doyle等[22,23]的研究,鋰電池液相中的電流和電勢符合歐姆定律,可以得到電荷的守恒方程為

式中f代表活化系數.邊界條件一是忽略隔膜兩側形成的電勢降,隔膜處鋰擴散通量為零,另一邊界條件則是集流體上鋰的擴散通量為零.因此液相區域的邊界條件可以表示為

2.2 熱模型建立

鋰離子電池在充放電過程中生熱源有電極極化引起的活化極化熱(Qact)、反應過程焓變引起的反應熱(Qrea)、歐姆電勢引起的歐姆熱(QOhm),由于集流體上的電阻率遠遠小于電極的電阻率,所以本模型中忽略集流體上產生的熱.

總熱量Q用(20)式表示為

根據能量守恒,

依據牛頓冷卻定律,鋰離子電池熱模型的邊界條件可以用(22)式來描述:

式中ρ,Cp,λ,Tamb,h分別代表密度、比熱容、熱導率、環境溫度、電池表面換熱系數,具體數值見表1.

2.3 容量衰減機理描述

本文所建立模型的衰減部分,是以電解液發生氧化還原反應,在負極顆粒表面生成SEI膜,消耗活性鋰為主要衰減機制.

容量衰減的建模過程描述如下.

首先采用塔菲爾方程描述在碳負極表面上電解液與活性鋰發生的副反應:

在循環過程中由于容量損失,下面將逐一對隨著循環過程而不斷更新的參數進行敘述.首先是電池的SOC修正:

SEI膜厚度的增長量Lgrow由法拉第定律可用下式進行計算:

因而在完成一個循環后SEI膜的總厚度為

隨著SEI膜的不斷增厚,膜阻抗Rsei也將不斷增大,膜阻抗表達式如下所示:

副反應過電勢ηside表達式如下:

(30)式中j為總電流密度,它包括兩部分:

考慮SEI膜電導率后,整個負極的電導率可以認為是初始負極電導與SEI膜電導的串聯,其表達式表示為

(32)式中κs,n,κs,n,0,κsei分別表示生成SEI膜后的負極電導率、原始負極電導率、SEI膜電導率.隨著循環周期進行,負極鋰離子濃度表示為

負極SEI膜阻抗Rsei表示為

3 模型參數與有效性驗證

3.1 模型參數

參數的獲取是仿真研究的前提和基礎,鋰電池的建模涉及的參數多而繁,分為基本設計參數、電極動力學參數、電極熱力學參數.此外,鋰離子電池在工作過程中伴隨著明顯的熱效應,電池溫度的實時變化,對鋰離子電池內部化學反應的相關參數及物料傳輸過程有著重要影響,因此相關參數的修正顯得尤為重要.

3.1.1 基本設計參數

對于仿真研究來說,仿真對象的重要物性參數尤為重要.表1為電池設計參數、初始鋰離子濃度參數及動力學傳遞參數.表2詳細列出了生成SEI膜副反應的相關參數.

表1 電化學-熱耦合模型采用的電芯物化參數Table 1.The cell parameter values used in electrochemical-thermal model.

表2 電池副反應相關參數Table 2.Values of parameters used for the solvent reduction reaction.

3.1.2 電極動力學參數

本文通過建立電化學熱耦合模型研究容量衰減行為,由于電池溫度在工況下波動較大,因而有必要考慮動態參數的變化.下面對擴散系數和反應速率等動力學參數,采用Arrhenius公式進行修正,詳見(35)—(39)式.

液相擴散系數De,電導率ke及動力學因數ν表達式分別見(40)—(42)式,其中c為電解液鋰鹽濃度.

3.1.3 電極熱力學參數

在電池放電過程中,電化學系統的反應驅動力是兩電極材料的電勢差,而電極材料的電勢差決定于不同電池材料的工作電壓.工作電壓Ui表達式如下:

Uref,i為電極參比電勢,參比電勢的大小則隨著嵌鋰量而變化.依據Verbrugge和Koch[25]以及Yamada等[26]的工作,可以得到參比電勢與嵌鋰量的函數關系式.對正負極上的SOC定義如下:

正極磷酸鐵鋰和負極石墨材料的開路電勢與SOC的關系式為:

3.2 模型驗證

本工作建立一個基于電解液分解消耗活性鋰導致容量衰減的數值模型,并從電化學角度探討影響循環容量的因素.其中,SEI膜的增長是利用與鋰離子電池的充放電循環時的主反應結合的副反應過程描述.目前為止并沒有有效的實驗方案和測量方法對副反應發生過程進行研究.因此,通過副反應的實驗數據來驗證模型的有效性不太可能.由于本文主要研究鋰離子電池在循環過程中容量衰減的情況,故本文通過容量衰減的數據側面驗證模型的有效性.

3.2.1 容量衰減驗證

從圖2中可看出,仿真模型所展示出的容量衰減與實驗所獲得的容量衰減情況在整體上具有很好的符合度.但在實驗后期,模擬結果與實驗結果出現偏差,表現為實驗容量的衰減程度略大于仿真數據.這可能是仿真模型中著重研究的是以SEI膜生長為主要衰減機制,而在實際情況中,由于集流體與負極的接觸電阻增大,負極材料中的導電劑、黏結劑等的效果變差等一系列原因,在后期會對循環過程造成比較大的影響.但從電化學,熱場的角度研究循環衰減情況,本模型依然具有很好的應用價值.

圖2 循環過程中的容量衰減驗證Fig.2.Validation of capacity fading during the cycling process.

4 結果與討論

本文主要研究對鋰離子電池進行了800次循環計算后電池的容量衰減情況及電化學行為表現,此外探討了充放電倍率、負極固相體積分數、活性物質顆粒粒徑等對電池循環壽命的影響.電池容量衰減主要體現在活性鋰的消耗及SEI膜阻抗的增大,文中以方程(26)局部副反應電流密度控制方程描述副反應發生的速率.在討論電極設計參數對循環容量影響的分析過程中,以電流密度在循環過程中的變化為量度進行定性分析,進行400周的循環計算做定量分析.

研究中電池的充放電循環制度如圖3所示:

1)1C(C為充放電倍率)恒流放電,截止電壓為2.0 V;

2)在2.0 V電壓下進行恒壓放電,直至電流下降到0.1 A·m?2以下;

3)1C恒流充電,截止電壓為3.6 V;

4)在3.6 V電壓下進行恒壓充電,直至電流下降到0.1 A·m?2以下;

5)電池開路直至單次循環時間達到10000 s;

6)重復步驟1),直至指定循環次數.

圖3 (網刊彩色)鋰離子電池充放電過程中電壓、電流與時間的關系Fig.3.(color online)Battery potential and current versus time in discharge progress.

4.1 循環過程中的電化學分析

根據模型仿真電池循環過程所得到的數據,主要從電化學角度分析鋰離子電池在循環過程中電池內部副反應電流密度的分布,以及電池負極初始SOC與放電容量隨循環次數增加的變化情況.

4.1.1 副反應電流密度變化分析

如圖4所示,在第一個循環中,前期電池在1C恒定電流下進行放電,在整個放電過程中副反應并沒有發生,副反應電流密度為0,大約在3600 s結束放電過程.接著對電池進行1C恒流充電,副反應電流密度開始急劇增大,在7000 s附近電流密度達到峰值.在之后的恒壓充電過程中,副反應電流密度急劇減小,在充電過程完成時,電流密度變為0,副反應消失.對比第1,第400,第800循環時副反應的電流密度變化曲線,可以明顯看出第1,第400,第800循環的副反應電流密度峰值分別為3.13,2.8,2.54μA·m?2,呈現遞減的規律.此外,圖4中每次循環的副反應密度曲線所形成的陰影面積代表的是在一個循環中單位面積上副反應所消耗的活性鋰量.可以看出,隨著循環進行,單個循環所消耗的鋰呈遞減的規律.

圖4 (網刊彩色)不同循環周期的副反應電流密度變化Fig.4.(color online)The side-reaction current density versus cycles.

4.1.2 負極不同位置副反應電流密度分析

如圖5(a)所示,縱坐標方向表示電池充電初期在不同時間點的負極上副反應電流密度的值,可以看出在放電初期(4200—4300 s),電池負極電流密度逐漸減小.橫坐標表示電池在負極方向上不同位置副反應電流密度的分布情況,可以看出,從集流體側到隔膜側有一個遞減趨勢.圖5(b)所示為充電末期(6900—7000 s)負極上副反應電流密度分布情況,與圖5(a)所示不同的是隨著充電過程的進行,負極上副反應電流密度逐漸變大;但在負極上從集流體側到隔膜側副反應電流密度的分布規律與充電初期相同.說明在整個的循環周期內活性鋰在負極上不同位置消耗都不同,造成負極材料上的活性材料利用率不高,最終降低了電池容量的輸出.

4.1.3 容量保持量與循環次數關系的分析

表3列出的是鋰離子電池在1C倍率下1,300,500,800次循環時,鋰離子電池的SEI膜阻抗、負極初始SOC和放電容量.可以看出,隨著循環次數的增加,SEI膜阻抗不斷增大,負極初始SOC及放電容量都在不斷減小.800次循環后,SEI膜阻抗增量為15.6 ?·m?2,負極SOC減少了0.052,約占電池初始SOC的6.29%;電池的放電容量減少了1.21 A·h,約占放電總容量的6.35%,鋰離子電池的負極初始SOC與容量分別減小到93.71%和93.65%.

圖5 不同時刻負極充電過程副反應電流密度分布(a)充電初期;(b)充電末期Fig.5.The side reaction current density in the anode:(a)At an early stage of charge progress;(b)at an end stage of charge progress.

表3 SEI膜阻抗、SOC及放電容量與循環次數的關系Table 3.Values of SEI fi lm resistance,SOC and the discharge capacity versus cycle number.

4.2 循環過程的容量衰減因素

該部分主要從充放電制度和電極設計方面探討了循環過程容量衰減情況,分析了充放電倍率、體積分數、負極活性顆粒半徑等對電池循環容量衰減的影響.

4.2.1 充放電倍率的影響

由圖6(a)中可以看出,所得副反應電流密度變化曲線與數值建模原理相符:只發生在充電過程.在1C倍率下,大約于放電結束(約3600 s)后的充電過程,副反應電流密度快速增大,直至充電過程結束.另外,副反應電流密度的絕對最大值1C<2C<3C<4C,副反應電流密度曲線與Y軸0值線所圍成的面積大小Area1C

圖6 (網刊彩色)不同充放電倍率下的副反應電流密度及容量衰減(a)副反應電流密度;(b)循環容量衰減Fig.6.(color online)The evolution of side-reaction current density and cycling capacity fade at di ff erent charge and discharge rates:(a)The side-reaction current density;(b)the cycling capacity fade.

4.2.2 體積分數的影響

在負極材料的組成上,體積分數是對電極的固液比例的一個表征.基于(12)式分析,可以看出當εs增大時,比表面積隨之增大.所以在電流源不變的情況下,電流密度減小,負極上的固相擴散極化也會減弱,副反應所消耗的容量減小.如圖7(a)示,1C充放電倍率下,第1次循環的負極固相體積分數在[0.3,0.6]區間時與副反應電流密度大小呈正相關關系,對應的副反應的電流密度的峰值分別為0.5,0.36,0.28,0.23 mA·cm?2,副反應電流密度曲線與Y軸0值線所圍成的面積大小Area0.3>Area0.4>Area0.5>Area0.6.較于[0.3,0.5]區間的體積分數增量對副反應的電流密度影響,很明顯固相體積分數從0.5增大到0.6時的影響力驟減.固相體積分數在[0.3,0.4]區間內時,400次循環后,可以清楚地看到電池容量都有不同程度的下降,且下降容量與固相體積分數成反比關系,體積分數為0.3時,損耗容量達到0.65 A·h.所以,本文所研究電池的最佳固相體積分數區間應在[0.5,0.6]之間.

圖7 (網刊彩色)不同負極固相體積分數的副反應電流密度及容量衰減(a)副反應電流密度;(b)循環容量衰減Fig.7.(color online)The evolution of side-reaction current density and cycling capacity fade along with the changes of the volume fraction of the cathode:(a)The side-reaction current density;(b)the cycling capacity fade.

4.2.3 負極活性顆粒粒徑的影響

圖8(網刊彩色)不同負極活性顆粒半徑副反應電流密度及循環容量衰減圖(a)副反應電流密度;(b)循環容量衰減Fig.8.(color online)The evolution of side-reaction current density and cycling capacity fade along with the changes of the particles radius of the cathode:(a)The side-reaction current density;(b)the cycling capacity fade.

圖8 (a)所示為1C充放電倍率下,改變負極材料顆粒半徑時第1次循環的副反應電流密度的變化對比.該圖中電流密度曲線與Y軸0值線所圍成的面積代表一個循環內鋰離子的損耗量.因此,曲線所圍成的面積越大,表示一個循環內的鋰離子的消耗量越大.圖8(a)中Area2×10?6

5 結論

本文基于負極充電過程持續生成SEI膜的機理,建立電化學-熱耦合容量衰減模型.基于該模型,研究了磷酸鐵鋰電池的容量衰減行為.研究表明:

1)在電池充電過程中,負極上各處的鋰濃度不同,呈現出從集流體側到隔膜側的遞減趨勢.負極上各處的副反應速率與鋰濃度在負極上的分布規律一致;

2)1C恒流恒壓充放電的制度下,800次循環后,磷酸鐵鋰電池SEI膜阻抗增量為15.6 m?·m?2,電池的放電容量減少了1.21 A·h,約占放電總容量的6.35%;

3)目標電池在400次循環后,較于1C充放電循環的容量衰減,2C,3C,4C容量衰減增量為0.62%,1.38%,1.73%;在1C充放電循環制度下,得出負極固相體積分數在[0.5,0.6]區間內能保持最好的電池循環穩定性和壽命發揮;在1C下,電池容量衰減率隨負極活性顆粒粒徑的減小而減小.在研究區間[2,10]μm內,粒徑為2和10μm時對應的容量衰減率分別為2.89%,3.87%,差值達到1%.因此,充放電倍率、負極固相體積分數、負極活性顆粒粒徑等因素都對循環壽命具有較大的影響,是提高鋰電池循環壽命的重要途徑.

[1]Yan J D 2014 Acta Aeronaut.Astronaut.Sin.35 2767(in Chinese)[閆金定2014航空學報35 2767]

[2]Etacheri V,Marom R,Ran E,Salitra G,Aurbach D 2011 Energy Environ.4 3243

[3]Arora P,White R E,Doyle M 1998 ChemInform 145 3647

[4]Gang N,Haran B,Popov B N 2003 J.Power Sources 117 160

[5]Vetter J,Novák P,Wagner M R,Veit C,Moller K C,Besenhard J O,Winter M,Mehrens W,Vogler C,Hammouche A 2005 J.Power Sources 147 269

[6]Laresgoiti I,K?bitz S,Ecker M,Du S 2015 J.Power Sources 300 112

[7]Arora P,Popov B N,White R E 1998 J.Electrochem.Soc.145 807

[8]Tang Y W,Jia M,Li J,Lai Y Q,Cheng Y,Liu Y X 2014 J.Electrochem.Soc.161 E3021

[9]Cheng Y,Li J,Jia M,Tang Y W,Du S L,Ai L H,Yin B H,Ai L 2015 Acta Phys.Sin.64 210202(in Chinese)[程昀,李劼,賈明,湯依偉,杜雙龍,艾立華,殷寶華,艾亮2015物理學報64 210202]

[10]Groot J,Swierczynski M,Stan A I,K?r S K 2015 J.Power Sources 286 475

[11]Ramadass P,Haran B,Gomadam P M,White R,Popov B N 2004 J.Electrochem.Soc.151 A196

[12]Ploehn H J,Ramadass P,White R E 2004 J.Electrochem.Soc.151 A456

[13]Safari M,Morcrette M,Teyssot A,Delacourt C 2009 J.Electrochem.Soc.156 A145

[14]Honkura K,Takahashi K,Horiba T 2011 J.Power Sources 196 10141

[15]Ecker M,Gerschler J B,Vogel J,K?bitz S,Hust F,Dechent P,Sauer D U 2012 J.Power Sources 215 248

[16]Watanabe S,Kinoshita M,Nakura K 2014 J.Power Sources 247 412

[17]Shi S Q,Gao J,Liu Y,Zhao Y,Wu Q,Ju W W,Ouyang C Y,Xiao R J 2016 Chin.Phys.B 25 018212

[18]Ye Y,Shi Y,Tay A A O 2012 J.Power Sources 217 509[19]Baek K W,Hong E S,Cha S W 2015 Int.J.Automot.Tech.16 309

[20]Newman J,Tiedemann W 1975 Aiche.J.21 25

[21]Doyle M,Fuller T F,Newman J S 1993 J.Electrochem.Soc.140 1526

[22]Doyle M,Newman J,Reimers J 1994 J.Power Sources 52 211

[23]Doyle M,Newman J 1995 Electrochimica Acta 40 2191

[24]Li J,Cheng Y,Jia M,Tang Y W,Lin Y,Zhang Z A,Liu Y X 2014 J.Power Sources 255 130

[25]Verbrugge M W,Koch B J 2003 J.Electrochem.Soc.150 A374

[26]Yamada A,Koizumi H,Nishimura S,Sonoyama N,Kanno R,Yonemura M,Nakamura T,Kobayashi Y 2006 Nature Mater.5 357

[27]Srinivasan V,Wang C Y 2003 J.Electrochem.Soc.150 A98

[28]Gerver R E,Meyers J P 2011 Quatern Int.158 A835

[29]Doyle M,Newman J,Gozdz A S,Schmutz C N,Tarascon J M 1996 J.Electrochem.Soc.143 1890

PACS:82.80.Fk,81.70.Pg,82.47.Aa,88.05.DeDOI:10.7498/aps.66.118202

Cyclic capacity fading of the power lithium ion battery based on a numerical modelling with dynamic responses?

Jiang Yue-Hui1)Ai Liang2)Jia Ming1)2)?Cheng Yun1)Du Shuang-Long1)Li Shu-Guo1)
1)(School of Metallurgy and Environment,Central South University,Changsha 410083,China)
2)(Hunan Aihua Group Co.,LTD,Yiyang 413002,China)

10 December 2016;revised manuscript

10 March 2017)

It is one of the important issues for electric vehicle to utilize power batteries which have long lifetime and excellent performance.For optimizing electrochemical performance and lifetime of the lithium ion battery,an electrochemicalthermal model based on dynamic response is developed by COMSOL MULTIPHYSICS.The modeling theory is the reaction mechanism of lithium iron phosphate battery which also includes a parasitic reaction occurring in the constant current and constant voltage charging process.The model consists of three parts:electro-chemical model,thermal model and capacity fade model.A series of temperature-dependent parameters and lithium ion concentration-dependent parameters relevant to the reaction rate and Li+transport are employed in this model.Comparing with the results of the experimental test,the model shows high accuracy and reliability.The capacity losses and electrochemical behaviors of the battery in cyclic processes with di ff erent rates are investigated.The results show that when the battery is cycled at a rate of 1C,the capacity fading rate is about 6.35%,meanwhile the solid electrolyte interface membrane impedance of the battery is increased by 15.6 m?·m?2after 800 time cycle.In the charge process,the side reaction rate within the anode shows a decreasing trend along the direction from the cooper current collector to separator,which is consistent with the lithium concentration in the anode.Besides,the e ff ects of charge/discharge rate,negative active material particle radius and negative solid volume fraction on the battery cycle life are also discussed respectively.Compared with the fading rate of 3.31%after 400 time cycle with 1C rate,the capacity fading rates for 2C,3C,4C reach to 3.93%,4.69%and 5.04%respectively.When the average particle radii of the anode are 2m and 10m,corresponding capacity fading rates are 2.89%and 3.87%,showing a di ff erence of nearly 1%.The study for solid volume fraction demonstrates that the battery with a solid volume fraction varying in a range of[0.5,0.6]will keep a longest battery life.These results show that the model has great potential to optimize the design of the battery.

capacity fade,charge and discharge rate,volume fraction,particle radius

10.7498/aps.66.118202

?國家自然科學基金(批準號:51204211,51222403)、工信部工業轉型升級強基工程專項(批準號:0714-EMTC02-5271/6)和湖南省

戰略性新興產業與新型工業化專項資金(批準號:2015GK1045)資助的課題.

?通信作者.E-mail:jiamingsunmoon@aliyun.com

?2017中國物理學會Chinese Physical Society

http://wulixb.iphy.ac.cn

*Project supported by the National Natural Science Foundation of China(Grant Nos.51204211,51222403),the Special Foundation for Industrial Upgrading Transformation and Strengthen the Foundation of Ministry of Industry and Information Technology,China(Grant No.0714-EMTC02-5271/6),and the Foundation of Strategic Emerging Industrial Scienti fi c Project Research,China(Grant No.2015GK1045).

?Corresponding author.E-mail:jiamingsunmoon@aliyun.com