章嘉晶 張蘭春



摘要:建立漸變的冷卻液分配流道,利用計算流體力學的方法,分析進出水口布置方式、分配流道面積、分配流道上底與下底比例對液冷板散熱性能的影響規律。結果表明:異側布置方式比同側布置方式的液冷板具有更低的電池最高溫度和溫差;4種分配流道面積中,1 455 mm2分配流道的液冷板下電池的最高溫度最低,溫差略高,整體冷卻表現最好;隨著分配流道面積從1 455 mm2增大至2 425 mm2,電池的最高溫度上升,溫差稍有下降;隨著分配流道上底與下底比例的增加,電池最高溫度先下降后上升,溫差逐漸增大,上底與下底比例為3:7的液冷板下電池最高溫度為304.346 K,溫差為3.302 K,綜合冷卻性能最優。
關鍵詞:鋰離子電池;液冷散熱;液冷板;分配流道
中圖分類號:TM 912.9文獻標識碼:A文章編號:2095-7394(2021)06-0079-10
隨著世界經濟的高速發展,能源危機日趨嚴重。目前,各國都致力于推動與發展新能源產業,開發電動汽車被公認為是一種替代傳統內燃機技術的可持續解決方案[1-2]。作為電動汽車的關鍵部件之一,鋰離子電池其成本和性能在很大程度上決定了電動汽車的成本和可靠性。鋰離子電池具有重量輕、儲能大、功率大、無污染等特點,但其高倍率放電時,生熱迅速,需采用合適的冷卻方式來保證電池的工作性能和使用壽命[3-5];而液冷散熱系統的傳熱效率高,液流可控性強,是主流的散熱方案[6-7]。當前,方形電池的液冷散熱方式主要有液冷板和直接浸入冷卻液兩種[8],其中,平行流道液冷板結構簡單、實用性強,因此,國內外學者對其進行了多項研究。
Wang等人[9]研究對比了平行流道液冷板的四種串并聯排布方式,發現3塊液冷板并聯的排布方式散熱效果最優;Tang等人[10]提出了三種不同的冷卻結構,通過實驗和仿真發現,電池模塊的底部和兩側都有冷卻板的結構具有最佳的冷卻性能,而增加散熱面積能有效降低電池的最高溫度和溫差;Bladimir Ramos-Alvarado等人[11]提出了四種平行流道結構,并研究了加入橫向流道后四種結構的冷卻性能;安治國等人[12]研究了不同流道數量、不同流道截面形狀以及不同截面長寬比的平行流道液冷板模型;T. Amalesh等人[13]研究了平行流道的七種變形結構(矩形槽流道、方波流道、波紋管流道、正弦波流道、弧形流道、圓形槽流道、鋸齒形流道)的散熱性能;Liu等人[14]在平行流道間加入傾斜的分岔流道,使之變形為樹狀流道結構,采用仿真模擬和實驗驗證的方法測試該結構的散熱性能,通過與傳統平行流道比較,表明樹狀流道液冷板具有更好的散熱性能和溫度均勻性。
然而,上述研究方向往往集中于液冷板的排布方式、數量,散熱流道的幾何參數、結構設計等問題,而鮮有文獻討論冷卻液分配流道對液冷板冷卻效果的影響。本文建立漸變的分配流道,將進水口、出水口兩側的冷卻液分配流道設為直角梯形,用來調整散熱流道內冷卻液的分布及流速,從而提高液冷板的散熱效果;針對直角梯形分配流道,從進出水口布置方式、分配流道面積、分配流道上底與下底的比例三個方面,研究分配流道結構參數對散熱性能的影響規律。
1模型與數值計算
1.1模型建立
本文的研究對象為38 Ah三元功率型方殼鋰電池,鋰離子電池內部結構復雜,因此,在建模過程中將電池簡化為單一實體,其基本尺寸為148.2 mm×97 mm×27 mm。如圖1所不,單體電池由貼合的兩塊液冷板進行冷卻,液冷板流道分為豎直的矩形散熱流道和直角梯形的冷卻液分配流道兩個部分。矩形散熱流道等間距分布;直角梯形分配流道a為上底長度,b為下底長度;進出水口居中設于下底處。
1.2電池生熱機理及熱物性參數
鋰離子電池產熱的主要來源有四種方式[15]:正負電極間的電化學反應熱Qr;電解液發生分解反應等產生的副反應熱Qs;電池內部歐姆電阻產生的焦耳熱Qj;容差極化內阻生成的極化熱Qp。電池的總產熱量為Q,可表示為:
Q=Qr+Qs+Qj+Qp。(1)
對電池表面進行隔熱處理和充放電實驗測試,認為電池在工作中產生的熱量全部被自身吸收并轉化為電池表面溫度升高,則電池總產熱量由以下公式計算:
Q=cpmΔT。(2)
式中:Q為鋰離子電池總產熱量;cp為電池定壓比熱容;m為電池質量;ΔT為電池在實驗過程中的溫升。
由式(2)得到生熱功率的計算公式為:
式中:P為電池生熱功率;t為電池溫升所對應的工作時間。
參考傅家麒[16]的實驗數據,電池在3 C倍率下放電1 200 s的生熱功率為21.86 W,則單體電池體熱源可表示為:
其中:q為電池體熱源即單位體積的生熱率;V為單體電池體積。電池、液冷板和液體的熱物性參數如表1所示。
1.3網格無關性及邊界條件
使用ICEM軟件對液冷模型進行非結構性四面體網格劃分,并對流體區域的網格加密得到1 087 587、2 077 712、3 475 400、4 889 382、5 853 239、6 834 242、10 266 768 七種不同數量的網格模型。以溫差、壓降的計算結果作為衡量網格無關性的依據,結果如圖2所示。對比488萬數量與585萬數量網格模型的計算結果,溫差和壓降的計算精度均低于1.5%,通過綜合考慮仿真準確性和計算資源,最終選擇488萬數量的網格模型。具體網格劃分如圖3所示。
采用Fluent作為模擬仿真工具,選擇層流模型計算;設置環境溫度為300 K;冷卻液進口設為速度入口,速度為0.03 m/s;出口設為壓力出口,壓力為0。由于電池布置在電池箱中,密封性好而換熱能力較差,故對流換熱系數取1 W/(m2·K)。
2液冷板散熱特性仿真及分析
2.1進出水口布置方式對散熱性能的影響
如圖4(a)所示,進、出水口均布置在液冷板的右側,稱為同側布置方式;如圖4(b)所示,進、出水口布置在液冷板的左右兩側,稱為異側布置方式。為對比兩種布置方式散熱性能的優劣,選取分配流道面積1 445 mm2,上底與下底比例為3:7,對不同散熱流道數量下兩種布置方式的液冷模組進行研究。
如圖5所示,為不同流道數量下兩種布置方式電池的最高溫度和溫差。(1)隨著散熱流道數量的增加,在同側布置方式下,電池的最高溫度和溫差均有明顯升高。(2)隨著散熱流道數量的增加,在異側布置方式下,電池的最高溫度逐漸降低,溫差逐漸增大;當散熱流道數量從7增加到9時,電池最高溫度趨于穩定,溫差增大得十分有限。(3)不同散熱流道數量中,異側布置方式下電池的最高溫度和溫差均低于同側布置方式。
如圖6所示,為兩種布置方式下電池的溫度分布圖。由圖6可知:同側布置方式下電池的最高溫度明顯高于異側布置方式,且溫度分布不均,說明異側布置方式的液冷板具有更優的散熱性能。
如圖7所示,為兩種布置方式下液冷板流速圖。分析圖7可知:在同側布置方式下,冷卻液集中分布在靠近進出水口的3條散熱流道,而其余散熱流道內冷卻液分布少、流速低、接近斷流;在異側布置方式下,散熱流道內冷卻液分布均勻,流速較快。可見,異側布置方式的液冷板具有更好的冷卻性能及溫均性;故選擇異側布置方式的液冷板為研究對象,探究分配流道結構因素對散熱性能的影響。
2.2分配流道面積對散熱性能的影響
以異側布置方式、分配流道上底與下底之比為3:7的液冷板為前提條件,對比分析不同散熱流道數量下,4種分配流道面積(單個直角梯形分配流道面積分別為970 mm2、1 455 mm2、1 940 mm2、2 425 mm2)的液冷板散熱性能。
如圖8所示,為不同分配流道面積的液冷板下電池最高溫度和溫差。由圖8可知:970mm2的分配流道液冷板下電池的最高溫度和溫差最高;當分配流道面積從1 455 mm2增大至2 425 mm2時,電池的最高溫度上升,溫差略有下降。
如圖9所示,為970 mm2、1 455 mm2、2 425 mm2三種分配流道面積液冷板下電池的溫度分布圖。由圖9可知:970 mm2分配流道液冷板下電池的溫度最高,且高溫區域分布最多;1 455 mm2分配流道液冷板下電池的溫度分布與2 425 mm2液冷板下電池的溫度分布基本一致;2 425 mm2液冷板下電池的溫度云圖中出現一個明顯的高溫點。上述3 種分配流道面積液冷板下,電池的溫度分布與圖8 所示的最高溫度、溫差規律相符。
如圖10所示,為三種分配流道面積的液冷板流速圖。分析圖10可知:冷卻液從進水口流入,流經分配流道,并沿途流入散熱流道;970mm2的分配流道上底與下底長度短,分配流道內冷卻液流速較快,阻礙了其順利流入中間的5條散熱流道,造成冷卻液整體分布不均,中間3條散熱流道內冷卻液流速慢,散熱效果不佳;在2 425 mm2的分配流道中,下底長度長,冷卻液從靠近下底側的進水口流入后,流速稍慢,從而能更均勻地分流至各條散熱流道,則2 425 mm2的分配流道下電池的溫差小,7條散熱流道內冷卻液流速相對較慢,電池的最高溫度較高;1 455 mm2分配流道的上底與下底長度適中,中間3條散熱流道內冷卻液流速稍慢,其余流道內流速較快,冷卻液整體分布均勻,則該方案下電池的最高溫度最低,溫差略高。
進一步分析圖8、圖9可知:1 455 mm2分配流道液冷板下電池的溫度分布與2 425 mm2液冷板下電池的溫度分布基本一致,溫差小,且均在5℃以內,說明兩者溫均性相近;但2 425 mm2分配流道液冷板下電池的溫度云圖中有明顯的高溫點,即電池的最高溫度高于1 455 mm2的方案,可見1 455 mm2的分配流道液冷板散熱性能更優。綜合上述因素,1 455 mm2分配流道液冷板的綜合冷卻性能在4種分配流道面積中最優。隨著散熱流道數量的增加,4種不同分配流道面積的電池最高溫度均下降,當散熱流道數量從7增加到9時,最高溫度趨于穩定;當散熱流道數量為3時,溫差較高,而散熱流道數量增加到4時,溫差有較大幅度下降;當散熱流道數量繼續增加,溫差緩慢上升。因此,綜合考慮最高溫度和溫差,選取1 455 mm2、7條散熱流道的液冷板進行后續研究。
2.3分配流道上底與下底的比例對散熱性能的影響
在異側布置方式、分配流道面積1 455 mm2、7 條散熱流道的液冷板條件下,研究分配流道上底與下底的比例對散熱性能的影響。如圖11所示,為分配流道上底與下底不同比例的液冷板下電池最高溫度和溫差。由圖11可見:隨著上底與下底比例的增加,電池最高溫度先下降,后上升;當比例為3:7時,電池最高溫度最低,最高溫度為304.346 K,溫差為3.302 K;隨著上底與下底比例的增加,溫差逐漸增大。
如圖12所示,為3種不同上底與下底比例液冷板下電池的溫度分布圖。由圖12可知:比例為1:9的液冷板下電池溫度分布均勻,高溫區域集中;比例為3:7的液冷板下電池溫度分布均勻,高溫區域狹長;比例為7:8液冷板下電池的高溫區域大,高溫區域與其余部分對比明顯,溫度分布不均。
如圖13所示為不同上底與下底比例的液冷板流速圖。分析圖13可知:(1)當比例為1:9時,下底長度較長,導致從靠近下底一側流入散熱流道的冷卻液流速較低;上底較短,致使冷卻液不能充分流入靠近上底的散熱流道中,則靠近上底的散熱流道內冷卻液流速略低,從而有利于冷卻液在散熱流道內均勻分布,電池的溫差小;散熱流道內冷卻液流速稍慢,故該比例液冷板下電池最高溫度較高。(2)當比例為7:8時,下底的長度較短,冷卻液從進水口流入分配流道后流速快,致使較少的冷卻液能流入中間散熱流道;冷卻液流速慢,特別是中間三條散熱流道內冷卻液分布少,接近斷流,故該比例下電池的最高溫度和溫差最高。(3)當比例為3:7時,分配流道上底與下底的比例適中,靠近上底與下底的散熱流道內冷卻液流速稍快,中間流道內流速略慢,在保證流速的同時也能使冷卻液在散熱流道內分布均勻,故最高溫度最低,溫差略高。綜上,比例為1:9的液冷板與比例為3:7的液冷板下,電池溫度分布都均勻,兩者的溫差變化小;但比例為1:9的液冷板下電池的最高溫度明顯要高。因此,比例為3:7的液冷板整體冷卻效果最好。
3結論
本文采用計算流體力學的方法,分析進出水口布置方式、分配流道面積、分配流道上底與下底比例對液冷板散熱性能的影響規律,得出以下結論。
(1)異側進出水口布置方式的液冷板比同側布置方式具有更好的冷卻性能和溫均性。在異側布置方式下,散熱流道內冷卻液分布均勻,流速較快;而同側布置方式下,冷卻液分布不均,流速相差較大。
(2)4種分配流道面積中,1 455 mm2的分配流道下電池的最高溫度最低,溫差略高,綜合冷卻性能最優;970 mm2的分配流道液冷板散熱效果最差;隨著分配流道面積從1 455 mm2增大至2 425 mm2,電池的最高溫度上升,溫差稍有下降。
(3)隨著分配流道上底與下底比例的增加,電池的最高溫度先下降后上升,溫差逐漸增大。當比例為3:7時,分配流道上底與下底的比例適中,在保證流速的同時也能使冷卻液在散熱流道內均勻分布,該比例的液冷板下電池的最高溫度為304.346 K,溫差為3.302 K;其最高溫度在五種比例中最低,溫差略高。因此,比例為3:7的液冷板綜合冷卻性能最優。
參考文獻:
[1]中華人民共和國中央人民政府.中華人民共和國國民經濟和社會發展第十二個五年規劃綱要[EB/OL].(2011- 03-16)[2021-01-20].http://www.gov.cn/2011lh/content_ 1825838.htm.
[2]工業和信息化部.促進汽車動力電池產業發展行動方案[EB/OL].(2017-03-17)[2021-02-12].http://www.xinhua- net.com//auto/2017-03/17/c_1120643683.htm.
[3] ZHE L,ZHANG J,WU B,et al. Examining temporal and spatial variations of internal temperature in large-format laminated battery with embedded thermocouples[J]. Journal of Power Sources,2013,241(1):536-553.
[4] KIM U S,SHIN C B,KIM C S. Effect of electrode configuration on the thermal behavior of a lithium- polymer battery [J]. Journal ofPower Sources,2008,180(2):909-916.
[5] LIU J,KUNZ M,CHEN K,et al. Visualization of charge distribution in a lithium battery electrode[J]. Journal of Physical Chemistry Letters,2010,1(14):2120-2123.
[6] ZHAO J,RAO Z,LI Y. Thermal performance of mini-channel liquid cooled cylinder- based battery thermal management for cylindrical lithium- ion power battery[J]. Energy Conversion & Management,2015,103:157-165.
[7] GIULIANO M R,PRASAD A K,ADVANI S G. Experimental study of an air- cooled thermal management system for high capacity lithium- titanate batteries[J]. Journal of Power Sources,2012,216:345-352.
[8] VANGILS R W,DANILOV D,NOTTEN PHL,et al. Battery ther - mal management by boiling heat -transfer[J]. Energy Conversion and Management,2014,79:9-17.
[9] WANG J,LU S,WANG Y,et al. Effect analysis on thermal behavior enhancement of lithium-ion battery pack with different cooling structures[J].Journal of Energy Storage,2020,32:101800.
[10] TANG A,LI J,LOU L,et al. Optimization design and numerical study on water cooling structure for power lithiumbattery pack[J].Applied Thermal Engineering,2019,159:113760.
[11] RAMOS-ALVARADO B,LI P,HONG L,et al. CFD study of liquid- cooled heat sinks with microchannel flow field configurations for electronics,fuel cells,and concentrated solar cells[J]. Applied Thermal Engineering,2015,31 (14):2494-2507.
[12]安治國,丁玉章,劉奇,等.液冷管道對動力鋰電池組溫度場影響研究[J].電源技術,2018,42(6):795-798.
[13] THANGAPACHAIYAPPAN A,NARASIMHAN N L. Introducing new designs of minichannel cold plates for the cooling of Lithium-ion batteries[J].Journal of Power Sources,2020:479.
[14] LIU H L,SHI H B,SHEN H,et al. The performance management of a Li-ion battery by using tree-like mini-channel heat sinks:experimental and numerical optimization [J]. Energy,2019,189:116150.
[15]桂長清.動力電池[M].2版.北京:機械工業出版社,2012.
[16]傅家麒.電動汽車高功率動力電池液冷系統熱均衡性能研究[D].鎮江:江蘇大學,2019.
Study on the Influence of Structure Parameters of Liquid Cooling Plate Distribution Channel on Heat Dissipation Performance
ZHANG Jiajing,ZHANG Lanchun
(School of Automotive and Traffic Engineering,Jiangsu University of Technology,Changzhou 213001,China)
Abstract:A gradual coolant distribution channel is established,and the influence of inlet and outlet layout,distribution channel area,and the proportion of upper and bottom of distribution channel on the heat dissipation performance of liquid cooling plate is analyzed by computational fluid dynamics method. The results show that the maximum temperature and temperature difference of the battery with different side arrangement are lower than those with the same side arrangement. Among the four distribution channel areas,the highest temperature of the battery under the liquid cooling plate with 1455 mm2distribution channel is the lowest,the temperature difference is slightly higher,and the overall cooling performance is the best. With the increase of the distribution channel area from 1455 mm2 to 2425 mm2,the maximum temperature of the battery increases and the temperature difference slightly decreases. With the increase of the proportion of upper bottom and lower bottom of the distribution channel,the maximum temperature of the battery decreases first and then increases,and the temperature difference gradually increases. The maximum temperature of the battery under the liquid cooling plate with a ratio of 3:7 is 304.346 K,and the temperature difference is 3.302 K,which indicates that the comprehensive cooling performance is optimal.
Key words:lithium ion battery;liquid cooling;liquid cold plate;distribution channel