COMSOL電化學(xué)耦合案例天花板!

電化學(xué)-熱耦合

COMSOL電化學(xué)耦合案例天花板!的圖1
COMSOL電化學(xué)耦合案例天花板!的圖2
COMSOL電化學(xué)耦合案例天花板!的圖3

COMSOL電化學(xué)耦合案例天花板!的圖4

本文主要介紹如何利用COMSOL Multiphysics 電化學(xué)-熱耦合模型對電動汽車動力電池的內(nèi)部溫度進(jìn)行數(shù)值模擬。

COMSOL電化學(xué)耦合案例天花板!的圖5
COMSOL電化學(xué)耦合案例天花板!的圖6
COMSOL電化學(xué)耦合案例天花板!的圖7

·

LiFePO4電池結(jié)構(gòu)可簡化為陽極極耳、含電解液的陽極集流體、陽極膜片、隔離膜、含電解液的陰極集流體和陰極極耳。電極隔膜和分離器由陽極膜片、陰極膜片和隔離膜組成。隔離膜中的聚合物骨架不導(dǎo)電子,電化學(xué)反應(yīng)只發(fā)生在固體活性物顆粒與電解液交界面處,其反應(yīng)方程式為:

陽極LixC6Li0C6xLi+xe-

陰極LixFe POC4xLi+xe-Fi0PO4

1、電化學(xué)模型

基于 Newman 的多孔性電極理論的電化學(xué)模型,其中描述正負(fù)電極顆粒表面電化學(xué)反應(yīng)過程的Buter-Volmer方程為:

COMSOL電化學(xué)耦合案例天花板!的圖8

j0為交換電流密度,單位為 A·cm-2;η 是局部過電位,單位為 V;αc和 αa是正負(fù)電極電化學(xué)反應(yīng)轉(zhuǎn)移系數(shù),取 0.5;F 為法拉第常數(shù),數(shù)值為96485 C·mol-1;R 為理想氣體常數(shù),數(shù)值為 8.314 J·mol-1· K-1

交換電流密度表達(dá)式為:

COMSOL電化學(xué)耦合案例天花板!的圖9

k0為反應(yīng)速率常數(shù);cs,max為材料最大固相鋰離子濃度;cs,surf為電極和電解液界面處鋰離子濃度。局部過電位表達(dá)式為:

COMSOL電化學(xué)耦合案例天花板!的圖10

?s為固相電勢;?1為液相電勢;Eeq為材料的平衡電位。

固相歐姆定律和液相歐姆定律分別為:

COMSOL電化學(xué)耦合案例天花板!的圖11

is和 i1分別代表固相電子電流和液相離子電流;σseff和 σ1eff分別代表固相和液相電勢;t+遷移數(shù)。

鋰離子在固相中的質(zhì)量守恒方程用菲克第二定律描述為:

COMSOL電化學(xué)耦合案例天花板!的圖12

cs代表固相鋰離子濃度;Ds為鋰離子的固相擴(kuò)散系數(shù);r 為球形顆粒的半徑。

鋰離子在液相中的質(zhì)量守恒方程用濃溶液理論描述為:

COMSOL電化學(xué)耦合案例天花板!的圖13

ε1為電極和隔膜中液相體積分?jǐn)?shù);D1eff為電解液中的等效擴(kuò)散系數(shù),用布拉格曼系數(shù)修正。

2、熱模型

鋰離子電池的熱模型的能量守恒方程為:

COMSOL電化學(xué)耦合案例天花板!的圖14

從左至右依次為熱積累、熱傳導(dǎo)和熱產(chǎn)生項(xiàng)。其中,電池電芯區(qū)域的產(chǎn)熱量主要包括可逆熱、極化熱與歐姆熱,依次為:

COMSOL電化學(xué)耦合案例天花板!的圖15

此外,極耳區(qū)域還存在由極耳的電阻所引起的歐姆熱,計(jì)算式為:

COMSOL電化學(xué)耦合案例天花板!的圖16

Itab為通過極耳處的電流;Atab為極耳的橫截面積,這里取極耳寬50 mm、厚0.3 mm;σtab為極耳自身材料的電導(dǎo)率,假定正極耳為鋁材料,負(fù)極耳為銅材料,電導(dǎo)率值參見表1電池電化學(xué)模型參數(shù)。經(jīng)過計(jì)算,求得正極耳的產(chǎn)熱率為 2.05×105W/m3負(fù)極耳的產(chǎn)熱率為 1.29×106W/m3

COMSOL電化學(xué)耦合案例天花板!的圖17

數(shù)

·

COMSOL電化學(xué)耦合案例天花板!的圖18

1、模型參數(shù)辨識

對模型中的主要參數(shù)進(jìn)行辨識。其中,對于導(dǎo)熱系數(shù)采取通用的熱阻串并聯(lián)的方法估算為:

法向?qū)嵯禂?shù)

COMSOL電化學(xué)耦合案例天花板!的圖19

展向?qū)嵯禂?shù)

COMSOL電化學(xué)耦合案例天花板!的圖20

Li為電芯單元中各部件的長度;Ki為各材料的導(dǎo)熱系數(shù)。

比熱容計(jì)算為:

COMSOL電化學(xué)耦合案例天花板!的圖21

根據(jù)表 1 中的具體數(shù)值,經(jīng)過計(jì)算得到法向導(dǎo)熱系數(shù)為 18.5 W·m-1·K-1,展向?qū)嵯禂?shù)為1.2 W·m-1·K-1,比熱容為1358 J·kg-1·K-1

表1 電池電化學(xué)模型參數(shù)

COMSOL電化學(xué)耦合案例天花板!的圖22

2、網(wǎng)格無關(guān)性驗(yàn)證

對模型進(jìn)行了網(wǎng)格無關(guān)性驗(yàn)證,消除網(wǎng)格數(shù)量對計(jì)算結(jié)果的影響。如圖 1 所示為放電電壓曲線,如圖 2 所示為電池表面溫度曲線。根據(jù)圖 1、圖 2,得出網(wǎng)格單元數(shù)為 44689、63479、75895 和108522,計(jì)算結(jié)果基本一致。綜合考慮計(jì)算量和準(zhǔn)確性,選擇網(wǎng)格單元數(shù)為 75895 的網(wǎng)格進(jìn)行后續(xù)計(jì)算。

COMSOL電化學(xué)耦合案例天花板!的圖23

圖1 網(wǎng)格無關(guān)性

(放電電壓曲線)

COMSOL電化學(xué)耦合案例天花板!的圖24

 圖2 網(wǎng)格無關(guān)性

(電池表面溫度曲線)

3、模型結(jié)果驗(yàn)證

COMSOL電化學(xué)耦合案例天花板!的圖25

圖3 電化學(xué)熱耦合模型

如 圖 3 所 示 為 磷 酸 鐵 鋰 電 池 一 維 模 型 。圖 3 從左至右依次為負(fù)極集流體、負(fù)極膜片、隔離膜、正極膜片和正集流體;其厚度分別為9 μm、59 μm、20 μm、9 μm 和 16 μm。

COMSOL電化學(xué)耦合案例天花板!的圖26

圖4 模型模擬的電池溫度與文獻(xiàn)結(jié)果的對比

圖 4(a)1D 模型在放電中后期溫度數(shù)據(jù)有較大偏差,這是由于 1D 模型假設(shè)內(nèi)部溫度一致且不考慮極耳對電芯的熱效應(yīng)而導(dǎo)致的。其中,3C 倍率偏差最大,最大平均偏差、絕對值平均偏差和均方根偏差分別為 7.4%、10.9%和14.3%。三者計(jì)算公式分別為

COMSOL電化學(xué)耦合案例天花板!的圖27

圖4(b)所示3D模型仿真數(shù)據(jù)與實(shí)驗(yàn)接近,其3C 倍率下最大平均偏差、絕對值平均偏差和均方根偏差分別為 3.8%、4.5%和 7.8%。因此,3D 模型可反應(yīng)電池內(nèi)部溫度分布,提供溫升曲線。

本文來自:COMSOL仿真交流

COMSOL電化學(xué)耦合案例天花板!的圖28

v

登錄后免費(fèi)查看全文
立即登錄
App下載
技術(shù)鄰APP
工程師必備
  • 項(xiàng)目客服
  • 培訓(xùn)客服
  • 平臺客服

TOP

5
7