LiNi0.5Mn1.5O4電極的變固相擴散系數(shù)模擬

呂桃林，羅 英，晏莉琴，解晶瑩

( 1.哈爾濱工業(yè)大學化工學院，黑龍江 哈爾濱 150001； 2.上海動力儲能電池系統(tǒng)工程技術(shù)有限公司，上海 200245； 3.上?？臻g電源研究所，上海 200245 )

·理論研究·

LiNi0.5Mn1.5O4電極的變固相擴散系數(shù)模擬

呂桃林1，2，3，羅 英1，2，3，晏莉琴2，3，解晶瑩1，2，3

( 1.哈爾濱工業(yè)大學化工學院，黑龍江 哈爾濱 150001； 2.上海動力儲能電池系統(tǒng)工程技術(shù)有限公司，上海 200245； 3.上?？臻g電源研究所，上海 200245 )

通過實驗和參數(shù)辨識，確定了針對LiNi0.5Mn1.5O4體系電池所需要的模型參數(shù)，建立了適用于該體系的變固相擴散系數(shù)模型。該模型能很好地對LiNi0.5Mn1.5O4的放電過程進行模擬。通過對電池內(nèi)部的過電位進行分析，對電池性能的限制因素進行了說明。

LiNi0.5Mn1.5O4； 變固相擴散系數(shù)(VSSD)模型； 過電位； 鋰離子電池

LiNi0.5Mn1.5O4由于具有高比容量及較高的電壓平臺(4.7 V)，引起了人們的極大興趣[1]。通過合適的數(shù)值模型，可深入理解電池內(nèi)部物理化學過程及限制因素，為電池的優(yōu)化設計提供參考。目前應用較廣泛的模型是J.S.Dunning等[2]建立的P2D模型，以多孔電極理論和濃溶液理論為基礎，對鋰離子電池內(nèi)部的物理化學過程進行描述[3]。該模型將固相傳輸過程視為單相擴散過程，未對材料在工作過程中可能發(fā)生的相變過程進行分析。LiNi0.5Mn1.5O4在充放電過程中有兩個電位平臺，分別對應Ni的兩次價態(tài)變化。研究表明：鎳錳酸鋰材料在Li+嵌脫過程中經(jīng)歷了相轉(zhuǎn)變過程。需要對現(xiàn)有模型進行修正，以建立適用于該體系的數(shù)值模型。涉及相變的傳統(tǒng)材料有磷酸鐵鋰、鈷酸鋰等，針對這些材料的相變過程，研究者提出了收縮核模型、相場模型和固溶體模型[4]和變固相擴散系數(shù)(VSSD)模型等[5]。VSSD模型避免對相變過程的具體模擬，同時包含該過程的主要特征。

本文作者建立了適用于LiNi0.5Mn1.5O4體系的VSSD模型，同時確定了模型中所需要的參數(shù)。對材料工作中的內(nèi)部過程進行分析，同時，對LiNi0.5Mn1.5O4的過電位進行分解，分析電池放電過程中的限制因素，為材料的優(yōu)化設計提供參考。

1 模型的建立

研究對象是以金屬鋰作為對電極的二電極體系，在P2D模型中引入VSSD模型。VSSD模型用于模擬固相內(nèi)部的傳輸過程，P2D模型用于模擬電極過程[5-6]。

VSSD模型基于二元溶液熱力學對固相內(nèi)部的傳輸過程進行描述，將活性物質(zhì)視為空位與鋰位構(gòu)成的非理想二元溶液。用吉布斯自由能分析它們之間的相互作用力，同時忽略充放電過程中的體積變化，將顆粒內(nèi)部的物質(zhì)傳遞表示為[5-6]：

(1)

式(1)中Ds′是二元擴散系數(shù)，為常數(shù)；cs是固相Li+濃度；α為固相熱力學因子，由電極的平衡電位得到，表示為[5-6]：

(2)

式(2)中：γ是固相內(nèi)Li+活度系數(shù)，Us為平衡電位，y為嵌鋰量。如果用稀溶液理論表示，根據(jù)Fick第二擴散定律求解，物質(zhì)傳遞表達式為[5-6]：

(3)

式(3)中：Ds為菲克擴散系數(shù)(即αDs′)。

利用Butler-Volmer方程描述電化學反應過程，用反應物質(zhì)流表示顆粒表面Li+的脫嵌過程，則粒子表面的反應離子流密度與表面上的過電勢的關(guān)系為：

(4)

式(4)中：JJ為活性顆粒表面離子流通量；αa、αc為反應傳遞系數(shù)；j0為交換電流密度；η為表面過電位。交換電流密度表示為：

(5)

η=Φs-Φe-Us

(6)

式(6)中：Φs為固相表面電位，Φe為液相電位。

2 實驗

2.1 電池的制作

將LiNi0.5Mn1.5O4(上海產(chǎn)，電池級)與乙炔黑(上海產(chǎn)，電池級)、聚偏氟乙烯(廣州產(chǎn)，電池級)按質(zhì)量比8∶1∶1混勻，加入適量溶劑N-甲基吡咯烷酮(廣州產(chǎn)，AR)攪拌均勻，涂覆在16 μm厚鋁箔(上海產(chǎn)，電池級)上，在80～120 ℃下干燥24 h，再沖成直徑為12 mm的極片。以金屬鋰片(上海產(chǎn)，電池級)為負極，1 mol/L LiPF6/EC+EMC(體積比3∶1，深圳產(chǎn)，電池級)為電解液，在高純氬氣保護的手套箱中組裝CR2032型扣式電池。

2.2 性能測試

組裝好的電池在CT2001A 5 V/2 mA充放電測試儀(武漢產(chǎn))上，以0.10C在3.5～4.9 V循環(huán)20次，進行活化。以極小倍率(0.01C)對電池進行充放電，通過對充放電數(shù)據(jù)進行歸一化處理，取平均值作為LiNi0.5Mn1.5O4的平衡電位；在不同荷電狀態(tài)(SOC)下(以0.6 mA放電2 min，不同狀態(tài)之間的間隔為0.02 mAh)進行不同倍率的恒流間歇滴定，通過充放電前后的平衡電位差值，對固相擴散系數(shù)進行計算；在不同SOC進行小倍率脈沖，通過放電時刻的電位階躍，對電化學反應系數(shù)進行計算；材料的初始嵌鋰量通過測量平衡電位得到。此外，由于部分參數(shù)會在電池工作過程中發(fā)生變化或不方便測量，通過不同倍率放電測試，利用參數(shù)辨識的方法進行整體調(diào)整，以確定模型中所需要的參數(shù)。所選倍率分別為：0.05C(0.03 mA)、0.10C、0.20C、0.50C和1.00C。所有測試均在恒溫(25 ℃)下進行。

3 結(jié)果與討論

3.1 參數(shù)提取

為方便研究，在對鎳錳酸鋰的平衡電位進行測試后，對結(jié)果進行擬合，如圖1所示。

圖1 LiNi0.5Mn1.5O4平衡電位曲線Fig.1 Equilibrium potential curves of LiNi0.5Mn1.5O4

從圖1可知，實驗過程對平衡電位進行了較好的擬合。LiNi0.5Mn1.5O4在放電過程中存在兩個電壓平臺，分別為4.73 V與4.69 V，標記為Ⅰ、Ⅱ；存在兩次較明顯的電壓轉(zhuǎn)折，標記為A、B，可能對應著LiNi0.5Mn1.5O4的結(jié)構(gòu)轉(zhuǎn)變。

為了對VSSD模型進行驗證，通過恒流間歇滴定，測量固相擴散系數(shù)隨嵌鋰量的變化情況。恒流間歇滴定示意圖見圖2。

圖2 恒流間歇滴定示意圖Fig.2 Galvanostatic intermittent titration technique

從圖2可知，恒流間歇滴定過程中，有5個電位變化，其中，U1為起始平衡電位，U2為滴定開始時的電位，U3為滴定過程中的電位變化，U4為滴定結(jié)束時的電位，U5為滴定結(jié)束后的平衡電位。通過測試電極電位隨時間的變化，通過以下公式計算材料的固相擴散系數(shù)[7]：

(7)

式(7)中：t為放電時間，mB為活性物質(zhì)質(zhì)量，VM為活性物質(zhì)的摩爾體積，MB為活性物質(zhì)相對分子質(zhì)量，S為極片的面積，△Et為放電過程中的電壓階躍(即U2-U3)，△Es為穩(wěn)態(tài)電位變化(即U1-U5)。

通過在不同SOC下對電池進行恒流間歇滴定(GITT)，可以得出LiNi0.5Mn1.5O4在放電過程中的擴散系數(shù)的變化情況，測試過程如圖3所示。

圖3 恒流間歇滴定(GITT)測試過程Fig.3 Galvanostatic intermittent titration(GITT) test procedure

可根據(jù)式(2)，通過平衡電位曲線計算得到的材料的熱力學因子，并結(jié)合菲克擴散系數(shù)，得出LiNi0.5Mn1.5O4在放電過程中的擴散系數(shù)變化情況。通過GITT以及熱力學因子計算所得結(jié)果見圖4。

圖4 固相擴散系數(shù)隨嵌鋰量的變化情況Fig.4 The change of solid diffusion coefficient with lithium intercalation

通過對比恒流間歇滴定的測試結(jié)果與根據(jù)式(2)計算的結(jié)果可知，在整個放電區(qū)間內(nèi)，兩種方式所得的固相擴散系數(shù)基本吻合，因此認為VSSD模型中假設的變固相擴散系數(shù)反映了LiNi0.5Mn1.5O4材料的本身性質(zhì)，可用于進行模擬。在放電末期出現(xiàn)較大偏差，可能是放電末期電池很難達到平衡所致。

從圖4可知，LiNi0.5Mn1.5O4的固相擴散系數(shù)隨嵌鋰量呈現(xiàn)馬鞍狀變化。通過與平衡電位的變化情況對比可知，在平衡電位變化較快的階段，擴散系數(shù)較大；在平衡電位平臺區(qū)域，擴散系數(shù)明顯減?。辉陔娢黄脚_階躍區(qū)域，擴散系數(shù)有所抬升。這可能與材料在不同嵌鋰狀態(tài)下，Li+之間的相互作用力有關(guān)，也可能與材料在充放電過程中的晶格變化有關(guān)。

此外，可通過對P2D模型的簡化，方便地對材料充放電過程中的電化學反應系數(shù)進行研究。首先假設在小電流放電過程中電化學反應分布符合單粒徑模型，即滿足式(8)：

(8)

定義參考電位為Φe|x=0=0，則可得金屬鋰的電極電位為：

(9)

式(9)中Φs，Li為金屬鋰的電極電位，J為放電電流密度，j0，Li為金屬鋰的交換電流密度。

其次，在放電電流較小的情況下對Butler-Volmer方程(4)進行線性簡化，即：

(10)

然后，根據(jù)過電位公式(6)，可將階躍電位表示為：

(11)

由于沒有液相濃度極化且反應均勻分布，因此可對液相電位分布進行簡化：

(12)

(13)

再結(jié)合式(10)—式(12)，將LiNi0.5Mn1.5O4的交換電流密度表示為：

(14)

最后，結(jié)合交換電流密度公式(5)，可將LiNi0.5Mn1.5O4的電化學反應系數(shù)表示為：

(15)

可通過小電流脈沖測試，對LiNi0.5Mn1.5O4的電化學反應系數(shù)進行估算。對電池進行小倍率脈沖放電時，電位變化曲線與圖3相同，其中△U為階躍電位。脈沖測試中采用較小電流，結(jié)果如圖5所示。

圖5 脈沖測試過程Fig.5 Impulse test procedure

根據(jù)式(15)對LiNi0.5Mn1.5O4在反應過程中的電化學反應常數(shù)進行計算，所得結(jié)果如圖6所示：

圖6 不同嵌鋰量下LiNi0.5Mn1.5O4的電化學反應常數(shù)Fig.6 Rate constant of the LiNi0.5Mn1.5O4 electrode as a function of y

從圖6可知，LiNi0.5Mn1.5O4的電化學反應常數(shù)在放電過程中基本維持在同一數(shù)量級，因此可根據(jù)該實驗結(jié)果，對仿真結(jié)果進行調(diào)整優(yōu)化。針對電池模型計算所需要的參數(shù)見表1。

在電池仿真過程中，以實驗電流值作為輸入工況，同時根據(jù)電池工作前的擱置電位作為平衡電位，計算電極起始嵌鋰量，參數(shù)ys，0(0.05C)、ys，0(0.10C)、ys，0(0.20C)、ys，0

表1 模型參數(shù)Table 1 Model parameters

(0.50C)和ys，0(1.00C)的取值分別為0.1261、0.1265、0.1285、0.1314和0.1291。

3.2 實驗仿真

根據(jù)所建立的VSSD模型，對電池不同倍率的放電曲線進行仿真，并與實驗所得數(shù)據(jù)進行對比，如圖7所示。

圖7 實驗數(shù)據(jù)與仿真數(shù)據(jù)對比Fig.7 Comparison of experiment data and simulation data

從圖7可知，鎳錳酸鋰在放電過程中有兩個明顯的放電平臺，且在放電末期有一個不明顯的電壓拐點，與該材料的平衡電位對應。隨著放電倍率的增大，受極化作用的影響，平臺拐點逐漸消失。VSSD模型的結(jié)果與實驗數(shù)據(jù)吻合，說明所建立的VSSD模型可真實地反映鎳錳酸鋰的工作過程。

3.3 內(nèi)部分析

圖8 不同倍率放電下固相濃度分布仿真結(jié)果Fig.8 Concentration profiles of intercalated lithium within particle simulated by the VSSD model over the course of discharge

從圖8可知，在放電過程中，材料經(jīng)歷了兩次界面濃度階躍過程。在放電初始階段，當材料表面嵌鋰量達到0.15左右時，由于擴散系數(shù)急劇下降，使得Li+在顆粒表面累積，界面濃度迅速升高，形成與固相內(nèi)部的濃差界面。在之后的放電過程中，該界面逐漸向內(nèi)遷移。當表面嵌鋰量達到0.5左右時，由于擴散系數(shù)增加，使得濃差界面消除。隨著放電繼續(xù)，當表面嵌鋰量達到0.6時，由于固相擴散系數(shù)又一次減小，使得濃差界面再次形成，但第二平臺上的擴散系數(shù)較第一平臺上擴散系數(shù)大，使得濃差界面的遷移并不明顯。當放電接近拐點，即表面嵌鋰量達到0.9左右時，由于擴散系數(shù)急劇升高，使得濃差界面消失，固相內(nèi)部濃度梯度減小。隨著放電倍率的增大，材料表面嵌鋰量累積速度增加，達到擴散系數(shù)降低階段的時間逐漸縮短，因此形成的濃差界面的時間逐漸縮短，且濃差界面逐漸明顯。同時，隨著倍率的增加，兩次濃差界面之間的過渡區(qū)間逐漸減小，這主要是由于隨著倍率的增加，濃度變化速度變快引起的。

3.4 電位分解

利用數(shù)值模型，可進一步對電池工作過程中的電位發(fā)展進行分析。在放電過程中，電池的電位是由理想平衡電位以及各種極化電位組合而成的：

(16)

通過VSSD模型，對1.00C放電過程中的電位分解與各極化過電位分析如圖9和圖10所示。

圖9 電池放電電壓分解Fig.9 Battery discharge potential decomposition

圖10 電池放電過電位分析Fig.10 Battery discharge over potential analysis

從圖9和圖10可知，在各極化過電位中，反應極化過電位與固相濃差極化過電位為主導。反應極化過電位主要受材料反應極化系數(shù)的影響；固相濃差極化過電位在放電末期急劇增大，主要受與固相擴散有關(guān)的材料粒徑及擴散系數(shù)的影響。

基于電池電位的分解，對影響電池主要性能的反應極化系數(shù)、固相擴散系數(shù)及粒徑的敏感度進行分析，參數(shù)取值見表2。

表2 模型輸入?yún)?shù)范圍Table 2 The range of model parameters

1.00C放電下各參數(shù)的敏感度如圖11所示。

圖11 1.00 C放電下不同參數(shù)的敏感度分析Fig.11 Sensitivity of different parameters at 1.00 C discharge

從圖11可知，反應系數(shù)通過影響反應極化過電位主要對電池放電過程中的放電平臺產(chǎn)生影響；固相擴散系數(shù)與材料粒徑通過影響固相濃差極化過電位在放電末期急劇增大，主要對電池的放電容量產(chǎn)生影響。目前針對LiNi0.5Mn1.5O4的改性研究主要集中在控制形貌、摻雜離子以及包覆等。可通過仿真對材料的合成改性提供一定的幫助，如通過控制材料形貌與粒徑，可改善材料的容量性能；通過摻雜提高材料的反應能力，可改善材料的功率特性，提高材料的功率密度。

4 結(jié)論

通過對LiNi0.5Mn1.5O4材料參數(shù)的測定辨識，建立了適用于該材料體系的變固相擴散系數(shù)模型，對LiNi0.5Mn1.5O4的工作過程進行了仿真。對放電過程中顆粒濃度的分布情況進行仿真，由于擴散系數(shù)的改變，在放電過程中顆粒內(nèi)部發(fā)生了兩次濃度階躍。通過模擬不同參數(shù)對電池性能的影響，對材料的限制因素進行分析，為材料的改性方向提供了一定的指導。

[1] SANTHANAM R，RAMBABU B.Research progress in high voltage spinel LiNi0.5Mn1.5O4material [J].J Power Sources，2010，195(17)：5 442-5 451.

[2] DUNNING J S，BENNION D N，NEWMAN J.Analysis of porous electrodes with sparingly soluble reactants [J].J Electrochem Soc，1971，118(8)：1 251-1 256.

[3] 羅玲，宋文吉，林仕立，etal.鋰離子電池熱模型的研究現(xiàn)狀 [J].電池，2015，45(5)：280-283

[4] SRINIVASAN V，NEWANM J.Discharge model for the lithium iron-phosphate electrode [J].J Electrochem Soc，2004，151(10)：1 517-1 529.

[5] FARKHONDEH M，DELACOURT C.Mathematical modeling of commercial LiFePO4electrodes based on variable solid-state diffusivity [J].J Electrochem Soc，2012，159(2)：177-192.

[6] FARKHONDEH M，SAFARI M，PRITZKER M，etal.Full-range simulation of a commercial LiFePO4electrode accounting for bulk and surface effects：a comparative analysis [J].J Electrochem Soc，2014，161(3)：201-212.

[7] LIU P，WU H.Diffusion of lithium in carbon [J].Solid State Ionics，1996，92(1-2)：91-97.

[8] RAHIMIAN S K，RAYMAN S.WHITE R E.Extension of physics-based single particle model for higher charge-discharge rates [J].J Power Sources，2013，224(4)：180-194.

收稿日期：2017-01-04

Modeling of LiNi0.5Mn1.5O4electrodes based on variable solid-state diffusivity

LV Tao-lin1，2，3，LUO Ying1，2，3，YAN Li-qin2，3，XIE Jing-ying1，2，3

(1.SchoolofChemicalEngineering&Technology，HarbinInstituteofTechnology，Harbin，Heilongjiang150001，China；2.ShanghaiPowerEnergyStorageBatterySystemEngineeringTechnologyCo.，Ltd.，Shanghai200245，China；3.ShanghaiInstituteofSpacePowerSources，Shanghai200245，China)

By adopting experiment and parameter recognition method，it was confirmed that a larger number of parameters which needed for the model of new system could setup a suitable model that involved variable solid-state diffusivity (VSSD).The model was used to simulate performance of LiNi0.5Mn1.5O4electrodes.The model was validated that it could be used to better simulate the discharging processes of LiNi0.5Mn1.5O4，by analyzing the over potential of battery and explaining its limiting factors.

LiNi0.5Mn1.5O4； variable solid-state diffusivity(VSSD) model； over potential； Li-ion battery

2016-11-08

呂桃林(1986-)，男，河南人，哈爾濱工業(yè)大學化工學院博士生，研究方向：鋰離子電池仿真；

10.19535/j.1001-1579.2017.02.006

TM912.9

A

1001-1579(2017)02-0084-06

羅 英(1987-)，女，重慶人，哈爾濱工業(yè)大學化工學院博士生，研究方向：鋰離子電池電解液；

晏莉琴(1975-)，女，湖北人，上海動力儲能電池系統(tǒng)工程技術(shù)有限公司高級工程師，研究方向：鋰離子電池系統(tǒng)；

解晶瑩(1971-)，女，黑龍江人，上?？臻g電源研究所總研究師，研究員，研究方向：鋰離子電池，本文聯(lián)系人。