葛子敬，臧孟炎，葉 鵬，謝金紅

(華南理工大學(xué)機(jī)械與汽車工程學(xué)院，廣東廣州 510640)

鋰離子電池組因其能量密度大、平均輸出電壓高、自放電率低等優(yōu)點(diǎn)廣泛應(yīng)用于電動(dòng)汽車領(lǐng)域。但是，電池組充放電過程中的大電流、電池組緊密的空間結(jié)構(gòu)以及惡劣的工作環(huán)境容易造成電池組溫升太高和溫度分布不均勻，嚴(yán)重時(shí)影響電池組的性能與壽命。因此，對電池組進(jìn)行散熱和溫度控制是保證電池性能的關(guān)鍵［1－2］。強(qiáng)制風(fēng)冷散熱系統(tǒng)［3］因其結(jié)構(gòu)簡單、質(zhì)量輕、成本低、能夠有效排除有害氣體、基本可以滿足電池組散熱要求等優(yōu)點(diǎn)廣泛應(yīng)用于電池組散熱系統(tǒng)。李騰等［4］詳細(xì)分析了鋰離子電池?zé)崮Ｐ偷难芯窟M(jìn)展，為設(shè)計(jì)鋰離子電池單體、模塊及熱管理系統(tǒng)提供依據(jù)。李哲等［5］著重研究了磷酸鐵鋰電池的容量特性、內(nèi)阻數(shù)值、開路電壓與環(huán)境溫度的關(guān)系，明確了磷酸鐵鋰電池的溫度特性。陳燕虹等［6］對電池組生熱和散熱溫度場進(jìn)行了仿真分析，得出電池組溫度場的分布情況，對電池組散熱系統(tǒng)的設(shè)計(jì)進(jìn)行了指導(dǎo)。劉皓等［7］通過計(jì)算流體動(dòng)力學(xué)(CFD)分析軟件對電池組的風(fēng)冷結(jié)構(gòu)進(jìn)行仿真分析，結(jié)果表明電池組在1C充放電條件下，電池組某一測點(diǎn)的溫度和溫差都處于正常范圍內(nèi)。

本文在上述研究的基礎(chǔ)上，采用風(fēng)冷散熱的方式對鋰離子電池組進(jìn)行冷卻，應(yīng)用仿真分析軟件FLUENT對電池組溫度場和流場進(jìn)行仿真分析，研究鋰離子電池組的風(fēng)冷散熱效果。

1 鋰離子電池生熱機(jī)理

鋰離子電池在進(jìn)行充放電過程中，會(huì)發(fā)生鋰離子和電子的嵌入與拖嵌，鋰離子和電子運(yùn)動(dòng)過程中產(chǎn)生的熱量稱為反應(yīng)熱Qr，鋰離子電池充電和放電過程為可逆反應(yīng)，反應(yīng)熱大小相等，符號相反;極化反應(yīng)產(chǎn)生極化熱Qp;鋰離子電池過充過放會(huì)引起副反應(yīng)，電解質(zhì)分解及自放電產(chǎn)生副反應(yīng)熱Qs;電池內(nèi)阻在充放電中產(chǎn)生焦耳熱Qj。由于電池管理系統(tǒng)的作用，Qs很小，可以忽略不計(jì)［8］。當(dāng)鋰離子電池溫度達(dá)到70～80℃時(shí)，反應(yīng)熱Qr占電池總熱量的較大比例;當(dāng)?shù)陀谏鲜鰷囟瘸浞烹姇r(shí)，焦耳熱Qj占較大比例。一般電池的正常工作溫度為－20～65℃［9］，因此鋰離子電池充放電過程中產(chǎn)生的熱量主要由極化熱和焦耳熱組成。

式中:Qtotal為鋰離子電池總的發(fā)熱量，J;Rp和Re分別為極化內(nèi)阻和電子流動(dòng)過程中內(nèi)阻，Ω;Rt為電池充放電過程中總內(nèi)阻，Ω;I為電池充放電過程中電流，A。

單位體積上電池自身的發(fā)熱功率即體熱源為:

式中:V為單體電池體積，m3。

2 計(jì)算模型的建立

2.1 鋰離子電池組模型

選用某款車載鋰離子電池為研究對象進(jìn)行建模，且在保證真實(shí)反映電池組內(nèi)部流場特性的前提下，對模型進(jìn)行一定簡化，并且只考慮由10個(gè)單體電池組成的電池組作為研究對象。單體電池尺寸為16mm×65mm×151mm，10個(gè)電池成單排排列，電池箱模型為如圖1所示的長方箱體結(jié)構(gòu)，空氣流道布置形式及空氣流道定義等如圖2所示(假設(shè)靠近入口端空氣流道為流道1，靠近入口端電池為電池1)，電池均勻排列且間距為6mm。使用ICEM軟件對鋰離子電池組幾何模型進(jìn)行網(wǎng)格劃分，對流體區(qū)域網(wǎng)格進(jìn)行加密處理，生成的結(jié)構(gòu)化網(wǎng)格如圖3所示，整個(gè)鋰離子電池組模型約有55萬個(gè)網(wǎng)格單元。

圖1 鋰離子電池組模型

圖2 電池箱縱斷面結(jié)構(gòu)

2.2 CFD仿真參數(shù)的設(shè)置

對鋰離子電池組風(fēng)冷散熱數(shù)值仿真時(shí)，提出如下假設(shè)條件［10］:1)在整個(gè)散熱過程中，流體看作不可壓縮的理想空氣，且不考慮空氣浮力;2)鋰離子電池設(shè)定為穩(wěn)定發(fā)熱熱源，入風(fēng)口為均勻進(jìn)風(fēng)條件;3)電池組輻射散熱影響小，不予考慮;4)流體與固體間無相對滑動(dòng)，接觸面設(shè)置成發(fā)熱量為0的耦合邊界條件。此外，假設(shè)電池單體材料質(zhì)地均勻，密度和比熱容不變;在同一方向上，電池的導(dǎo)熱系數(shù)相等且不受SOC和溫度的影響;電池充放電時(shí)，電池內(nèi)部電流密度均勻［11－12］。電池和空氣介質(zhì)的物性參數(shù)見表1。

圖3 電池組網(wǎng)格模型

表1 介質(zhì)物性參數(shù)

仿真分析參數(shù)設(shè)定:進(jìn)風(fēng)口空氣溫度為25℃，入口條件為velocity－inlet，初始風(fēng)速為5m/s，空氣入口端雷諾數(shù)為10 558，選用k－ε湍流模型求解;出口條件為outflow;計(jì)算域分為空氣流體域和電池固體域，采用標(biāo)準(zhǔn)壁面函數(shù)方法計(jì)算壁面與充分發(fā)展的湍流區(qū)域之間的黏性影響區(qū)域;空氣與電池表面進(jìn)行流固耦合換熱，電池箱外部進(jìn)行自然對流散熱，對流散熱系數(shù)為5W/(m2·K)。

2.3 CFD計(jì)算控制方程

鋰離子電池組數(shù)值仿真與模擬涉及到固體熱傳導(dǎo)、流體熱傳導(dǎo)以及流體與固體的耦合換熱，計(jì)算需要用到如下控制方程。

1)有內(nèi)熱源的電池組導(dǎo)熱控制方程。

式中:ρ1為電池密度;Cp為電池比熱容;T為電池溫度;t為溫度變化的時(shí)間;kx，ky，kz為電池在 x，y，z 3個(gè)方向上的熱導(dǎo)率;q為體熱源。

2)空氣冷卻控制方程。

包括以下3個(gè)方程:

①連續(xù)性方程。

②動(dòng)量方程。

③能量方程。

式中:ρ為空氣密度;u，v，w 為速度矢量 U 在 x，y，z 3個(gè)方向上速度分量;μ為空氣的動(dòng)力黏度;Su，Sv，Sw為動(dòng)量守恒方程的廣義源項(xiàng);k為流體的傳熱系數(shù);Cp為流體比熱容;ST為黏性耗散項(xiàng)。

3 仿真分析結(jié)果

3.1 風(fēng)冷散熱系統(tǒng)穩(wěn)態(tài)計(jì)算與瞬態(tài)計(jì)算比較

對電池組在1C充放電倍率(對應(yīng)電池組體熱源為20 000W·m－3)情況進(jìn)行強(qiáng)制風(fēng)冷散熱仿真分析。首先在FLUENT中進(jìn)行穩(wěn)態(tài)計(jì)算，穩(wěn)態(tài)計(jì)算迭代監(jiān)測曲線如圖4所示，得出仿真結(jié)果在計(jì)算210步時(shí)迭代收斂，此時(shí)電池組最高溫度為33.79℃，溫差為5.48℃。之后采用瞬態(tài)計(jì)算，迭代時(shí)間設(shè)置為3 600s，時(shí)間步長為1s，瞬態(tài)計(jì)算迭代監(jiān)測曲線如圖5所示，其中第10塊電池中心點(diǎn)p溫度監(jiān)測曲線如圖6所示，發(fā)現(xiàn)中心點(diǎn)溫度在3 600s之前已經(jīng)收斂，即瞬態(tài)計(jì)算在3 600s之內(nèi)已經(jīng)計(jì)算收斂，此時(shí)電池組最高溫度為33.70℃，溫差為5.59℃。瞬態(tài)計(jì)算結(jié)果和穩(wěn)態(tài)計(jì)算相比相差不大，但由于穩(wěn)態(tài)計(jì)算占用內(nèi)存少、計(jì)算時(shí)間短、計(jì)算效率高、收斂性好，因此采用穩(wěn)態(tài)計(jì)算方式仿真分析電池組風(fēng)冷散熱效果。

圖4 穩(wěn)態(tài)迭代收斂性曲線

圖5 瞬態(tài)迭代收斂性曲線

圖6 點(diǎn)p溫度檢測曲線

3.2 風(fēng)冷散熱系統(tǒng)仿真結(jié)果

電池組溫度云圖如圖7所示，電池組最高溫度為33.79℃，溫差為5.48℃，溫度云圖顯示電池左上半部分溫度較高，與預(yù)期結(jié)果吻合。各個(gè)單體電池最高溫度和溫差曲線如圖8所示，從分析結(jié)果可知，最高溫度與溫差曲線有很高的一致性，且靠近入口端電池溫度較低，遠(yuǎn)離入口端電池溫度較高。輸出各流道空氣流速的曲線如圖9所示，靠近入口端流道內(nèi)空氣流速較大，與電池耦合換熱帶走的熱量較多，從側(cè)面證實(shí)了靠近入口端電池溫度較低。但是總體來看各流道空氣流速相差不大，各流道較均勻的空氣流量分配使得各單體電池溫度均勻性較好，電池組最高溫度和溫差基本都處于正常工作范圍內(nèi)，即強(qiáng)制風(fēng)冷散熱系統(tǒng)可以滿足電池組散熱要求。

圖7 電池組溫度云圖

圖8 各單體電池溫度監(jiān)測曲線

圖9 各流道空氣流速監(jiān)測曲線

3.3 進(jìn)風(fēng)口空氣流速的影響

空氣掠過平板強(qiáng)制對流換熱系數(shù)h為:

式中:Nu為努塞爾數(shù);Re為雷諾數(shù);Pr為普朗特?cái)?shù);λ為流體熱導(dǎo)率;l為定型尺寸;ρ為流體密度;v為流體流速;L為特征長度;u為流體黏度;m，a，b為流動(dòng)狀態(tài)相關(guān)系數(shù)。

空氣掠過平板時(shí)，流體處于湍流狀態(tài)，此時(shí)系數(shù)a=0.8 ＜ 1，則h∝v0.8，即強(qiáng)制對流換熱系數(shù)h隨著流速v增大而增大，但是增大的幅度在不斷減小。

進(jìn)風(fēng)口空氣流速分別選取5m/s、10m/s、20m/s、30m/s和40m/s并建立相應(yīng)仿真模型進(jìn)行數(shù)值計(jì)算，電池組最高溫度和溫差與進(jìn)風(fēng)流速關(guān)系曲線如圖10、圖11所示。由曲線可知，隨著進(jìn)風(fēng)流速的增大，電池組最高溫度和溫差都在持續(xù)降低，電池組散熱效果越來越好。

圖10 不同進(jìn)風(fēng)流速下各單體電池最高溫度曲線

圖11 不同進(jìn)風(fēng)流速下各單體電池溫差曲線

電池組散熱效果提升率曲線如圖12所示。由圖12可知，當(dāng)流速低于20m/s時(shí)，電池組的散熱效果隨著流速的增大顯著增強(qiáng);當(dāng)流速高于20m/s后，流速的增大對電池組的散熱效果的提升率的影響越來越小。該仿真結(jié)果與式(9)得出的理論結(jié)果相吻合，符合真實(shí)性。

圖12 電池組散熱效果提升率曲線

雖然在空氣流速足夠大時(shí)，可以將電池充放電過程中產(chǎn)生的熱量有效地散出，但在實(shí)際應(yīng)用中，由于空氣流速的增大，電池組散熱效果提升率越來越小，電池管理系統(tǒng)效率越來越低，同時(shí)還提高了系統(tǒng)的能耗并增大了噪聲。因此在保證電池組良好散熱的條件下，將進(jìn)風(fēng)口空氣流速限定在一定范圍內(nèi)很有必要。

3.4 環(huán)境溫度的影響

電池組在充放電過程中，環(huán)境溫度的變化也會(huì)對電池組風(fēng)冷散熱效果產(chǎn)生影響。當(dāng)電池箱外環(huán)境溫度在－10～40℃之間變化時(shí)，電池組的最高溫度和溫差的變化曲線如圖13所示。由圖可知，環(huán)境溫度越低，電池組散熱效果越好，電池組最高溫度越低、溫差越小，與預(yù)測結(jié)果相符。

圖13 不同環(huán)境下最高溫度與溫差變化曲線

4 結(jié)束語

本文以計(jì)算流體動(dòng)力學(xué)方法對鋰離子電池組風(fēng)冷散熱結(jié)構(gòu)進(jìn)行建模和數(shù)值仿真分析。仿真結(jié)果表明進(jìn)風(fēng)口風(fēng)速與環(huán)境溫度對電池組散熱效果影響顯著，變化趨勢與理論結(jié)果相一致，對風(fēng)冷散熱結(jié)構(gòu)設(shè)計(jì)有一定指導(dǎo)作用;但同時(shí)也有一定局限性，需要相關(guān)實(shí)驗(yàn)驗(yàn)證仿真結(jié)果的準(zhǔn)確性。因而將在后續(xù)工作中盡可能地搭建風(fēng)冷散熱平臺，以進(jìn)一步完善風(fēng)冷散熱結(jié)構(gòu)的研究設(shè)計(jì)。

［1］ 劉振軍，林國發(fā)，秦大同，等.電動(dòng)汽車鋰電池組溫度場研究及其結(jié)構(gòu)優(yōu)化［J］.汽車工程，2012(1):80－84.

［2］ 平平.鋰離子電池?zé)崾Э嘏c火災(zāi)危險(xiǎn)性分析及高安全性電池體系研究［D］.合肥:中國科學(xué)技術(shù)大學(xué)，2014.

［3］ 付正陽，林成濤，陳全世.電動(dòng)汽車電池組熱管理系統(tǒng)的關(guān)鍵技［J］.公路交通科技，2005，22(3):119 －123.

［4］ 李騰，林成濤，陳全世.鋰離子電池?zé)崮Ｐ脱芯窟M(jìn)展［J］.電源技術(shù)，2009(10):927－932.

［5］ 李哲，韓雪冰，盧蘭光，等.動(dòng)力型磷酸鐵鋰電池的溫度特性［J］.機(jī)械工程學(xué)報(bào)，2011(18):115－120.

［6］ 陳燕虹，吳偉靜，劉宏偉，等.純電動(dòng)汽車電池箱的熱特性［J］.吉林大學(xué)學(xué)報(bào):工學(xué)版，2014(4):925－932.

［7］ 劉皓，楊凱，惠東，等.鋰離子電池組風(fēng)冷結(jié)構(gòu)設(shè)計(jì)［J］.電源技術(shù)，2011(10):1208－1210.

［8］ 姬芬竹，劉麗君，楊世春，等.電動(dòng)汽車動(dòng)力電池生熱模型和散熱特性［J］.北京航空航天大學(xué)學(xué)報(bào)，2014(1):18－24.

［9］ 張志杰，李茂德.鋰離子動(dòng)力電池溫升特性的研究［J］.汽車工程，2010(4):320－323.

［10］王福軍.計(jì)算流體動(dòng)力學(xué)分析:CFD軟件原理與應(yīng)用［M］.北京:清華大學(xué)出版社有限公司，2004.

［11］李騰，林成濤，陳全世.鋰離子電池?zé)崮Ｐ脱芯窟M(jìn)展［J］.電源技術(shù)，2009，33(10):927 －932.

［12］朱聰，李興虎，宋凌珺.電動(dòng)汽車用鋰離子電池生熱速率模型［J］.汽車工程，2014(2):174－180.