鋰離子電池仿生樹狀通道液冷板數(shù)值優(yōu)化

劉歡歡，任曉龍，張澤斌*

鋰離子電池仿生樹狀通道液冷板數(shù)值優(yōu)化

劉歡歡1,2，任曉龍2，張澤斌1*

（1.鄭州大學(xué) 機(jī)械與動(dòng)力工程學(xué)院，鄭州 450001；2.廣東佛山聯(lián)創(chuàng)工程研究生院，廣東 佛山 528300）

將鋰離子電池在5放電倍率、環(huán)境溫度298.15 K的工作條件下的最大溫度和溫差控制在適宜范圍內(nèi)。在鋰電池兩側(cè)插入帶通道的仿生樹狀液冷板，建立電池模塊液冷散熱模型。首先驗(yàn)證液冷數(shù)值模型的準(zhǔn)確性，然后通過(guò)正交試驗(yàn)設(shè)計(jì)研究分支流道角度、入口流速、分支流道寬度和液冷板厚度對(duì)冷板散熱性能的影響，并通過(guò)極差分析對(duì)4種影響因素進(jìn)行重要程度排序，確定最佳的參數(shù)組合。經(jīng)優(yōu)化后，鋰電池的最高溫度為302.4 K、溫差為3.4 K，滿足設(shè)計(jì)要求，但冷卻劑在流動(dòng)循環(huán)過(guò)程中存在泄漏風(fēng)險(xiǎn)。提出的樹狀液冷板滿足鋰電池?zé)峁芾硇阅艿囊?，?yōu)化后電池模塊的表面溫度在安全范圍內(nèi)，同時(shí)電池的溫度均勻性顯著提高。應(yīng)選擇導(dǎo)熱系數(shù)高的冷板材料，并注重冷板結(jié)構(gòu)設(shè)計(jì)，防止冷卻劑泄漏。

鋰離子電池；液冷板；樹狀結(jié)構(gòu)；正交試驗(yàn)

鋰離子電池具有能量密度高、循環(huán)壽命長(zhǎng)等特點(diǎn)，在儲(chǔ)能系統(tǒng)中得到了廣泛應(yīng)用。儲(chǔ)能系統(tǒng)中的鋰電池在放電、充電過(guò)程中會(huì)產(chǎn)生熱，并且難以及時(shí)散發(fā)熱量，易引起鋰電池組之間的熱量聚集，出現(xiàn)運(yùn)行溫差過(guò)大等現(xiàn)象，過(guò)高的溫度會(huì)導(dǎo)致鋰電池起火甚至爆炸[1]。分析近年來(lái)儲(chǔ)能系統(tǒng)的安全事故發(fā)現(xiàn)，火災(zāi)隱患多由鋰電池?zé)崾Э匾餥2]。由此可見，開發(fā)高效的鋰電池?zé)峁芾硐到y(tǒng)對(duì)于鋰電池模塊的性能和使用壽命至關(guān)重要，鋰電池?zé)峁芾硐到y(tǒng)應(yīng)保證鋰電池模塊之間的溫差在5 K以內(nèi)，最高溫度保持在293.15～313.15 K的最佳范圍內(nèi)[3]。目前，鋰電池?zé)峁芾硐到y(tǒng)有4種常用的冷卻方法：空氣冷卻、液體冷卻、相變材料冷卻和熱管冷卻[4]。研究人員還提出了一種混合熱管理系統(tǒng)，通過(guò)結(jié)合2種常用的冷卻方法來(lái)提高鋰電池模塊的冷卻效率[5]。

與其他冷卻系統(tǒng)相比，液冷系統(tǒng)需要水泵、管道、冷板等配件，導(dǎo)致其重量和成本增加，同時(shí)冷卻劑存在泄漏風(fēng)險(xiǎn)。由于液體冷卻系統(tǒng)具有較高的傳熱系數(shù)和較好的導(dǎo)熱性能，目前已成為鋰電池?zé)峁芾硌芯康闹攸c(diǎn)。劉霏霏等[6]分析了并行流道冷板的入口流速、通道數(shù)量、通道寬度等參數(shù)對(duì)鋰電池最高溫度和溫差的影響，結(jié)果表明，通道數(shù)量對(duì)冷板散熱性能的影響最顯著。付平等[7]通過(guò)比較傳統(tǒng)蛇形結(jié)構(gòu)冷板、外置式雙進(jìn)出水口蛇形冷板和內(nèi)嵌式雙進(jìn)出水口蛇形冷板的散熱性能后發(fā)現(xiàn)，雙進(jìn)出水口的內(nèi)置式結(jié)構(gòu)是一種較理想的冷板結(jié)構(gòu)。Li等[8]采用仿生波紋并行通道冷板，研究了放電速率、流道形狀、冷卻液質(zhì)量流量等因素對(duì)鋰電池最高溫度和溫差的影響，發(fā)現(xiàn)具有仿生波紋流道的冷板在冷卻效率上優(yōu)于并行直流道冷板。Wang等[9]采用正交試驗(yàn)設(shè)計(jì)方法研究了蛛網(wǎng)式通道冷板對(duì)鋰電池散熱性能的影響。在蛛網(wǎng)形通道的結(jié)構(gòu)參數(shù)中，通道寬度對(duì)冷板冷卻性能的影響最大，通道角度對(duì)冷板冷卻性能的影響最小。杜柏林等[10]建立了相變材料/液冷復(fù)合的鋰電池?zé)峁芾硐到y(tǒng)，并結(jié)合多目標(biāo)優(yōu)化算法得出了鋰電池組最高溫度和溫差最小的參數(shù)組合。

傳統(tǒng)的蛇形流道冷板存在功耗大的缺點(diǎn)，并行通道冷板存在溫度梯度大的缺點(diǎn)，因此在設(shè)計(jì)冷板結(jié)構(gòu)時(shí)需要綜合考慮散熱效果和功耗情況。自然系統(tǒng)存在的流動(dòng)結(jié)構(gòu)（如樹冠、葉脈、血管等）可為設(shè)計(jì)新的冷板提供思路，這些流道系統(tǒng)在很大程度上遵循默里定律[11]。默里定律給出了樹狀分支結(jié)構(gòu)在2個(gè)連續(xù)分支層次上分支直徑的最優(yōu)關(guān)系。

文中提出一種樹狀流道結(jié)構(gòu)冷板，首先通過(guò)單因素分析分支流道角度、寬度、冷板厚度和冷卻液流速對(duì)冷板散熱性能的影響，并采用正交試驗(yàn)設(shè)計(jì)、極差分析法對(duì)各因素進(jìn)行優(yōu)化分析，以得到各因素影響冷板散熱性能的重要程度及最佳參數(shù)組合。

1 模型和方法

1.1 幾何模型

選取8個(gè)方形鋰電池和9個(gè)液冷板組成鋰電池組模塊。鋰電池的長(zhǎng)度、寬度、高度分別為118、63、13 mm，鋰電池組模塊結(jié)構(gòu)如圖1a所示。在2個(gè)相鄰的鋰電池之間存在鋁制樹狀液冷板，液冷板的長(zhǎng)度和寬度分別等于鋰電池的長(zhǎng)度和寬度。由于鋰電池組采用幾何對(duì)稱排列方式，因此采用簡(jiǎn)化的鋰電池模塊。選擇2個(gè)1/2鋰電池和1個(gè)液冷板作為計(jì)算域，如圖1b所示，以降低計(jì)算成本，將頂部和底部表面設(shè)置為對(duì)稱面。樹狀冷板的二維結(jié)構(gòu)如圖1c所示。模型中使用的材料熱物性參數(shù)如表1所示。

圖1 鋰電池?zé)峁芾矸桨甘疽鈭D

表1 材料熱物性參數(shù)

Tab.1 Thermo-physical properties of materials

1.2 控制方程

冷卻劑的質(zhì)量、動(dòng)量、能量守恒方程[12]見式（1）—（3）。

冷卻劑雷諾數(shù)見式（4）。

式中：為冷卻通道水力直徑，m。

液冷劑的壓降計(jì)算見式（5）。

式中：為摩擦因數(shù)；為流道的長(zhǎng)度，m；ξ為局部損失系數(shù)。

液冷劑與鋰電池的換熱量的計(jì)算見式（6）。

式中：為對(duì)流換熱系數(shù)，W/(m2·K)；為換熱面積，m2；1、w分別為冷卻劑、壁面的溫度，℃。

1.3 初始和邊界條件

采用Ansys Fluent求解器，通過(guò)計(jì)算可知冷卻劑流動(dòng)入口的雷諾數(shù)均小于2 300，因此將冷卻劑流動(dòng)設(shè)置為層流模型條件。將環(huán)境溫度設(shè)定為298.15 K，鋰電池表面為自然對(duì)流，對(duì)流傳熱系數(shù)為5 W/(m2·K)。入口邊界為速度條件，出口邊界為零壓力條件。

1.4 產(chǎn)熱模型驗(yàn)證

采用Rao等[13]建立的鋰電池產(chǎn)熱模型，以恒定電流放電720 s來(lái)計(jì)算體積發(fā)熱率，放電速率為5（表示鋰電池的放電倍率），鋰電池的體積發(fā)熱率為240 000 W/m3。從圖2可知，得到的鋰電池最高溫度與文獻(xiàn)[13]中鋰電池的最高溫度吻合較好。由此可見，這里建立的鋰電池發(fā)熱模型合理，后續(xù)研究將以該模型為基礎(chǔ)。

2 數(shù)值仿真

2.1 條件假設(shè)

1）冷卻劑流動(dòng)是不可壓縮的。

2）鋰電池是一個(gè)均勻的產(chǎn)熱部件。

3）將液冷板的上下表面視為無(wú)厚度表面[14]。

圖2 產(chǎn)熱模型的驗(yàn)證

2.2 網(wǎng)格數(shù)量無(wú)關(guān)性驗(yàn)證

文中在對(duì)模型進(jìn)行網(wǎng)格劃分時(shí)采用四面體網(wǎng)格，為了檢驗(yàn)網(wǎng)格數(shù)量的無(wú)關(guān)性，采用6種網(wǎng)格數(shù)量（分別為25萬(wàn)、33萬(wàn)、42萬(wàn)、56萬(wàn)、90萬(wàn)、127萬(wàn)）進(jìn)行仿真分析。作為測(cè)定網(wǎng)格數(shù)量無(wú)關(guān)性的初始方案，將試驗(yàn)條件設(shè)置為冷卻劑入口流速0.1 m/s、分支流道角度45°、分支流道寬度2 mm、冷板厚度2 mm，后面在進(jìn)行仿真計(jì)算時(shí)若無(wú)特別說(shuō)明，則試驗(yàn)條件不變。計(jì)算結(jié)果如圖3所示，當(dāng)網(wǎng)格數(shù)量從56萬(wàn)增至127萬(wàn)時(shí)，溫度最大變化量約為0.1 K，相對(duì)誤差在0.03%內(nèi)。為了在確保計(jì)算精度足夠的同時(shí)減少時(shí)間成本，選擇網(wǎng)格數(shù)量56萬(wàn)進(jìn)行后面的仿真計(jì)算。

圖3 鋰電池最大溫度的網(wǎng)格獨(dú)立性試驗(yàn)結(jié)果

2.3 液冷數(shù)值模型驗(yàn)證

為了研究數(shù)值模型的準(zhǔn)確性，對(duì)鋁制微通道散熱器進(jìn)行了仿真模擬，F(xiàn)azeli等[15]進(jìn)行了相應(yīng)實(shí)驗(yàn)，模型幾何結(jié)構(gòu)如圖4a所示。在散熱器底面施加恒定的熱流，進(jìn)水流速為0.071 m/s，進(jìn)口冷卻液溫度為27 ℃。將模擬結(jié)果沿底面中心線的溫度與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行了比較，結(jié)果如圖4b所示?？梢钥闯?，最大相對(duì)誤差小于4%，因此文中的模擬結(jié)果可靠。

圖4 液冷數(shù)值方法驗(yàn)證

3 仿真結(jié)果分析

3.1 分支流道角度的影響

鋰電池的溫度和溫差隨時(shí)間的變化曲線如圖5所示。從圖5a可以看出，在前40 s內(nèi)，各方案中鋰電池的最高溫度基本相同，升溫速率在初始時(shí)刻最大，之后逐漸減小，這是因?yàn)槔鋮s劑流過(guò)冷板通道的每個(gè)角落并與鋰電池進(jìn)行對(duì)流換熱需要一定的時(shí)間。在40 s后，各方案中鋰電池的最高溫度開始出現(xiàn)差異，即分支流道角度對(duì)冷板散熱性能產(chǎn)生的影響開始顯現(xiàn)。

在分支流道角度為30°時(shí)，鋰電池的最高溫度、溫差分別為306.2、6.7 K，此時(shí)冷板的散熱性能最優(yōu)。在不同分支流道角度下，鋰電池的最高溫度均在理想工作范圍內(nèi)，但溫差大于5 K。當(dāng)分支流道角度從30°增至60°時(shí)，鋰電池的最高溫度增加了1.5 K，溫差增加了1.7 K。在分支流道角度大于60°后，鋰電池的最高溫度、溫差分別僅增加了0.1、0.2 K，可見在分支流道角度大于60°后對(duì)樹狀冷板散熱性能的影響較小。不同分支流道角度下鋰電池溫度分布云圖如圖6所示?？梢钥闯?，鋰電池溫度沿著冷卻劑的流動(dòng)方向不斷升高，這是因?yàn)槔鋮s劑不斷被鋰電池加熱，導(dǎo)致?lián)Q熱速率不斷減小。從圖5、圖6、表2中也可看出，隨著分支流道角度的增加，最高溫度、溫差和壓降都呈增大的趨勢(shì)。這是因?yàn)榱鲃?dòng)換熱面積隨著分支角的增大而減小，阻礙了冷卻劑與鋰電池的對(duì)流換熱。

圖5 不同分支流道角度下鋰電池的最高溫度和溫差隨時(shí)間的變化

圖6 不同分支流道角度下鋰電池溫度分布云圖

表2 不同分支流道角度下冷卻劑的壓降

Tab.2 Pressure drop of coolant under different branch runner angles

3.2 冷卻劑入口流速的影響

在不同流速下，鋰電池的最高溫度、溫差隨時(shí)間的變化過(guò)程如圖7所示，鋰電池的溫度分布云圖如圖8所示。可以看出，當(dāng)進(jìn)口速度為0.04 m/s時(shí)，鋰電池的最高溫度達(dá)到313.9 K，超出了鋰電池的理想工作溫度范圍。在其他5個(gè)速度下鋰電池的最高溫度均低于308.1 K，只有當(dāng)流速大于0.20 m/s時(shí)，鋰電池的溫差才能保持在5 K以下。隨著流速的增加，鋰電池的最高溫度不斷降低，這是因?yàn)榱魉俚脑黾訉?dǎo)致對(duì)流傳熱系數(shù)增大，在相同時(shí)間內(nèi)冷卻劑可以帶走更多的熱量。

當(dāng)流速?gòu)?.04 m/s增至0.16 m/s時(shí)，最高溫度和溫差顯著降低，最高溫度降低了9.2 K，溫差降低了8.4 K。當(dāng)入口流速?gòu)?.16 m/s變化到0.24 m/s時(shí)，最高溫度和溫差僅分別降低了1.3、1.1 K ，鋰電池的最高溫度和溫差緩慢下降。這是因?yàn)槔鋮s劑在流道入口段的局部傳熱系數(shù)較高，沿著流動(dòng)方向逐漸降低，并趨于穩(wěn)定。入口段的長(zhǎng)度與雷諾數(shù)呈正相關(guān)，當(dāng)入口流速?gòu)?.04 m/s增至0.16 m/s時(shí)，與入口流速?gòu)?.16 m/s變化到0.24 m/s時(shí)相比速度的變化范圍大，這意味著雷諾數(shù)增加的值較大，因此平均傳熱系數(shù)增加得較多。由式（6）可知，冷卻劑帶出的熱量較多，因此最高溫度和溫差變化明顯。由此可見，增大冷卻劑流速可以顯著改善鋰電池模塊的冷卻效果，但是進(jìn)一步增大入口流速會(huì)減緩鋰電池模塊的冷卻效果，同時(shí)會(huì)增大冷卻劑的壓降。如表3所示，當(dāng)冷卻劑的入口流速?gòu)?.04 m/s增至0.16 m/s時(shí)，速度增加了0.12 m/s，冷板的壓降增加了68.8 Pa。冷卻劑的入口流速?gòu)?.16 m/s增至0.24 m/s，雖然僅增加了0.08 m/s，然而冷板的壓降增大了67.4 Pa。從式（5）可知，壓降與摩擦因數(shù)和流體速度有關(guān)，并隨著速度呈指數(shù)變化，因此二者的壓降增加量相近。由以上分析可知，當(dāng)入口流速超過(guò)一定值時(shí)，增加入口流速并不能顯著改善液冷板的冷卻性能，而壓降會(huì)迅速增大，這意味著泵的功耗會(huì)迅速增大。

圖7 不同入口流速下鋰電池的最高溫度和溫差隨時(shí)間的變化

圖8 不同入口流速下鋰電池溫度分布云圖

表3 不同入口流速下冷卻劑的壓降

Tab.3 Pressure drop of coolant under different inlet flow rates

3.3 分支流道寬度的影響

最高溫度隨分支流道寬度變化的曲線如圖9a所示，隨著流道寬度的增加，鋰電池的最高溫度呈下降趨勢(shì)，在最初的100 s內(nèi)，鋰電池在不同分支流道寬度下的最高溫度基本相同。在100～400 s時(shí)，鋰電池的升溫速率逐漸減小，并趨于穩(wěn)定。在400 s后，最高溫度出現(xiàn)較長(zhǎng)時(shí)間的平臺(tái)期，這意味著由鋰電池產(chǎn)生的熱量可以及時(shí)排出，此時(shí)鋰電池的產(chǎn)熱率和散熱率處于平衡狀態(tài)。

當(dāng)流道寬度從1 mm增至2 mm時(shí)，鋰電池的最高溫度僅降低了0.2 K，當(dāng)流道寬度從2 mm增至3.5 mm時(shí)，鋰電池的最高溫度基本保持不變。這是因?yàn)樵黾恿鞯缹挾龋瑩Q熱面積增加得不大，因此對(duì)換熱量的影響較小。流道寬度對(duì)溫差的影響與最高溫度的影響相似，如圖9b所示，當(dāng)分支流道寬度從1 mm增至2 mm時(shí)，溫差減小了0.5 K；當(dāng)分支流道寬度大于2 mm后，溫差僅減小了0.2 K。由此可知，當(dāng)分支流道寬度大于2 mm后對(duì)冷板的散熱性能基本無(wú)影響。不同分支流道寬度下鋰電池的溫度分布云圖如圖10所示，可以看出不同流道寬度下溫度的分布相似，并呈現(xiàn)對(duì)稱性。雖然不同分支流道寬度下鋰電池的最高溫度均滿足要求，但溫差均大于5 K。流道寬度的變化主要會(huì)影響冷板的壓降，從表4可以看出，隨著寬度的增加，壓降呈負(fù)相關(guān)變化。

圖9 不同流道寬度下鋰電池的最高溫度和溫差隨時(shí)間的變化

圖10 不同分支流道寬度下鋰電池的溫度分布云圖

表4 不同分支流道寬度下冷卻劑的壓降

Tab.4 Pressure drop of coolant under different branch runner width

3.4 冷板厚度的影響

鋰電池的最高溫度和溫差隨著冷板厚度變化的曲線如圖11所示。隨著冷板厚度的增加，最高溫度呈現(xiàn)減小的趨勢(shì)。這是因?yàn)榱鞯栏叨扰c冷板厚度相同，隨著冷板厚度的增加，冷卻劑的流量隨之增大，在相同時(shí)間內(nèi)冷卻劑可以帶走更多的熱量。溫差隨時(shí)間變化的趨勢(shì)與最高溫度隨時(shí)間變化的趨勢(shì)相似，隨著冷板厚度的增加，溫差不斷減小，并且僅在板厚為4 mm時(shí)，鋰電池的最高溫度在適宜區(qū)間的同時(shí)溫差也小于5 K。

不同冷板厚度下鋰電池的溫度分布云圖如圖12所示。板厚從1 mm增至2.5 mm時(shí)，鋰電池的最高溫度、溫差分別降低了5.8、5.9 K，當(dāng)板厚大于2.5 mm后，鋰電池的最高溫度和溫差僅降低了1.8 K，鋰電池的最高溫度和溫差變化率不斷減小，這意味著冷板厚度對(duì)其散熱性能的影響逐漸減弱。由式（6）可知，冷卻劑帶出的熱量主要與對(duì)流換熱系數(shù)和冷卻面積有關(guān)，局部對(duì)流換熱系數(shù)在入口段較大，且沿著流動(dòng)方向逐漸減小。厚度增大，雷諾數(shù)增大，入口段長(zhǎng)度增大，這意味著平均換熱系數(shù)增大，同時(shí)換熱面積也增大。由于雷諾數(shù)和換熱面積的增加速率逐漸減小，因此厚度對(duì)冷板散熱性能的影響逐漸減弱。壓降的變化見表5，隨著冷板厚度的增加，壓降一直減小。上述結(jié)果表明，增加冷板厚度可以提高樹狀冷板的散熱性能，但增加到一定值后繼續(xù)增大冷板厚度對(duì)散熱性能的影響不會(huì)明顯增大。

圖12 不同冷板厚度下鋰電池溫度分布云圖

表5 不同冷板厚度下冷卻劑的壓降

Tab.5 Pressure drop of coolant under different cold plate thickness

4 正交試驗(yàn)分析

4.1 正交試驗(yàn)設(shè)計(jì)

上述研究表明，鋰電池的最高溫度（max）和溫差（Δ）與分支流道角度（因子A）、入口流速（因子B）、分支流道寬度（因子C）、冷板厚度（因子D）相關(guān)。為了使鋰電池達(dá)到更好的工作狀態(tài)，需進(jìn)一步優(yōu)化該冷卻模型。通過(guò)正交試驗(yàn)設(shè)計(jì)來(lái)安排和分析多因素試驗(yàn)，并對(duì)部分具有代表性的試驗(yàn)結(jié)果進(jìn)行分析，可以了解全部因素試驗(yàn)的情況，進(jìn)而找出各個(gè)因素水平的最佳組合[16]。這里使用具有4個(gè)因素和4個(gè)水平的正交試驗(yàn)來(lái)研究上述因素對(duì)鋰電池模塊散熱效果的影響，設(shè)計(jì)的正交試驗(yàn)及結(jié)果見表6。

表6 L16正交試驗(yàn)結(jié)果

Tab.6 L16 orthogonal test results

4.2 鋰電池的最高溫度、溫差及優(yōu)化結(jié)果分析

在正交試驗(yàn)中，為了評(píng)估各因素對(duì)樹狀通道冷板鋰電池?zé)峁芾硐到y(tǒng)最高溫度或溫差的影響，得到使鋰電池最高溫度或溫差最小的參數(shù)組合，這里采用極差分析法對(duì)表6中的正交試驗(yàn)結(jié)果進(jìn)行分析。在不同因素水平下，鋰電池的最高溫度曲線如圖13a所示，試驗(yàn)因素對(duì)應(yīng)的極差范圍越大，表示該因素對(duì)鋰電池最高溫度的影響越大[17]。因子A、B、C、D的極差范圍分別為1.45、5.33、1.20、7.01，因此這4個(gè)因素對(duì)鋰電池最高溫度的重要性排序?yàn)镈＞B＞A＞C。針對(duì)提出的樹狀通道冷板，使鋰電池的最高溫度最低的參數(shù)組合是D4B4A2C2。不同因素水平下鋰電池的溫差曲線如圖13b所示，可以看出使鋰電池的溫差最小的參數(shù)組合也是D4B4A2C2。

初始方案和正交優(yōu)化方案的數(shù)值模擬結(jié)果如圖 14所示，參數(shù)組合為D4B4A2C2時(shí)，鋰電池的最高溫度為302.4 K，溫差為3.4 K。與初始方案相比，鋰電池的最高溫度降低了4.4 K，溫差降低了4.1 K，可見D4B4A2C2方案為各個(gè)因素在研究范圍內(nèi)的最佳參數(shù)組合。

圖13 不同因素水平下鋰電池的最大溫度和溫差

圖14 不同方案下鋰電池最高溫度和溫差

5 結(jié)語(yǔ)

設(shè)計(jì)了一種新型的樹狀流道液冷板，通過(guò)數(shù)值模擬方法分析了分支流道角度、入口流速、分支流道寬度和冷板厚度對(duì)液冷板散熱性能的影響，主要結(jié)論如下。

1）冷板厚度對(duì)鋰電池散熱性能的影響最大，增加液冷板厚度可以顯著降低鋰電池的溫度和溫差，并減小壓降。增大分支流道角度，則鋰電池的最高溫度和溫差會(huì)增大，分支流道角度在大于60°后對(duì)冷板的散熱性能基本無(wú)影響。分支流道寬度主要影響冷板的壓降，對(duì)鋰電池溫度的影響較小。

2）在一定流速范圍內(nèi)，增大流速可以有效降低鋰電池的最高溫度，并改善其溫度均勻性。當(dāng)入口流速增加到一定值時(shí)，對(duì)鋰電池最高溫度和溫差的影響減弱，并且增大流速會(huì)使壓降迅速增加。

3）通過(guò)正交試驗(yàn)優(yōu)化后，鋰電池的最高溫度、溫差分別為302.4 K、3.4 K。與初始方案相比，鋰電池的最高溫度、溫差分別降低了4.4、4.1 K，有效降低了鋰電池的表面溫度，提高了溫度均勻性。

文中僅考慮了4個(gè)因素對(duì)液冷板散熱性能的影響，后續(xù)將考慮更多的影響因素。同時(shí)，應(yīng)對(duì)冷板進(jìn)行結(jié)構(gòu)設(shè)計(jì)，防止冷卻劑泄漏，并結(jié)合代理模型、多目標(biāo)優(yōu)化算法進(jìn)行全局優(yōu)化設(shè)計(jì)，以找出使冷板綜合性能最優(yōu)的參數(shù)組合。

[1] CHEN Zhao-liang, YANG Shu, PAN Min-qiang, et al. Experimental Investigation on Thermal Management of Lithium-Ion Battery with Roll Bond Liquid Cooling Plate[J]. Applied Thermal Engineering, 2022, 206: 118106.

[2] TAO Feng-bo, ZHANG Wei-jiang, GUO Dong-liang, et al. Thermofluidic Modeling and Temperature Monitoring of Li-Ion Battery Energy Storage System[J]. Applied Thermal Engineering, 2020, 181: 116026.

[3] ZHOU Hao-bing, ZHOU Fei, XU Li-peng, et al. Thermal Performance of Cylindrical Lithium-Ion Battery Thermal Management System Based on Air Distribution Pipe[J]. International Journal of Heat and Mass Transfer, 2019, 131: 984-998.

[4] 楊梓堙, 竺玉強(qiáng), 王亞平, 等. 基于PCM-液冷復(fù)合的鋰離子鋰電池?zé)峁芾硌芯縖J]. 電源技術(shù), 2022, 46(12): 1388-1392.

YANG Zi-yin, ZHU Yu-qiang, WANG Ya-ping, et al. Research on Thermal Management of Lithium-Ion Battery Based on PCM-Liquid Cooling Compound[J]. Chinese Journal of Power Sources, 2022, 46(12): 1388-1392.

[5] PATEL J R, RATHOD M K. Recent Developments in the Passive and Hybrid Thermal Management Techniques of Lithium-Ion Batteries[J]. Journal of Power Sources, 2020, 480: 228820.

[6] 劉霏霏, 陳洋洋, 胡建泰, 等. 雙冷板液冷鋰離子鋰電池溫升特性研究[J]. 汽車技術(shù), 2022(12): 15-22.

LIU Fei-fei, CHEN Yang-yang, HU Jian-tai, et al. Study on Temperature Rise Characteristics of Liquid-Cooled Lithium-Ion Battery with Double Cold Plates[J]. Automobile Technology, 2022(12): 15-22.

[7] 付平, 畢杰, 楊明飛. 液冷板結(jié)構(gòu)對(duì)鋰電池模組液冷性能的影響研究[J]. 電源技術(shù), 2021, 45(3): 313-316.

FU Ping, BI Jie, YANG Ming-fei. Influence Study of Liquid Cooling Plate Structure on Liquid Cooling Performance of Lithium Battery Module[J]. Chinese Journal of Power Sources, 2021, 45(3): 313-316.

[8] LI Bo, MAO Zhao-yong, SONG Bao-wei, et al. Study on Battery Thermal Management of Autonomous Underwater Vehicle by Bionic Wave Channels with Liquid Cooling[J]. International Journal of Energy Research, 2021, 45(9): 13269-13283.

[9] WANG Jian-feng, LIU Xiao-dong, LIU Fen, et al. Numerical Optimization of the Cooling Effect of the Bionic Spider-Web Channel Cold Plate on a Pouch Lithium-Ion Battery[J]. Case Studies in Thermal Engineering, 2021, 26: 101124.

[10] 杜柏林, 張甫仁, 李世遠(yuǎn), 等. CPCM/液冷復(fù)合鋰電池?zé)峁芾矸绞絻?yōu)化設(shè)計(jì)[J]. 包裝工程, 2022, 43(3): 202-209.

DU Bo-lin, ZHANG Fu-ren, LI Shi-yuan, et al. Optimization Design of Thermal Management Mode of CPCM/Liquid-Cooled Composite Battery[J]. Packaging Engineering, 2022, 43(3): 202-209.

[11] JING Da-lei, SONG Shi-yu, HE Lei. Reexamination of Murray's Law for Tree-Like Rectangular Microchannel Network with Constant Channel Height[J]. International Journal of Heat and Mass Transfer, 2019, 128: 1344-1350.

[12] ZUO Wei, ZHANG Yun-tian, JIAQIANG E, et al. Effects of Multi-Factors on Performance of an Improved Multi-Channel Cold Plate for Thermal Management of a Prismatic LiFePO4Battery[J]. Energy, 2022, 261: 125384.

[13] RAO Zhong-hao, WANG Qing-chao, HUANG Cong- liang. Investigation of the Thermal Performance of Phase Change Material/Mini-Channel Coupled Battery Thermal Management System[J]. Applied Energy, 2016, 164: 659-669.

[14] ZHAO Ding, LEI Zhi-guo, AN Chao. Research on Battery Thermal Management System Based on Liquid Cooling Plate with Honeycomb-Like Flow Channel[J]. Applied Thermal Engineering, 2023, 218: 119324.

[15] FAZELI S A, HOSSEINI HASHEMI S M, ZIRAKZADEH H, et al. Experimental and Numerical Investigation of Heat Transfer in a Miniature Heat Sink Utilizing Silica Nanofluid[J]. Superlattices and Microstructures, 2012, 51(2): 247-264.

[16] 王學(xué)章, 李科群. 鋰電池叉流流道液冷結(jié)構(gòu)設(shè)計(jì)及散熱特性分析[J]. 物理學(xué)報(bào), 2022, 71(18): 184702.

WANG Xue-zhang, LI Ke-qun. Liquid-Cooled Structure Design and Heat Dissipation Characteristics Analysis of Cross-Flow Channels for Lithium Batteries[J]. Acta Physica Sinica, 2022, 71(18): 184702.

[17] QI Wen-jie, HUANG Wen-qi, NIU Jun-tian, et al. Thermal Management of Power Battery Based on Flexible Swiss Roll Type Liquid Cooling Micro-Channel[J]. Applied Thermal Engineering, 2023, 219: 119491.

Numerical Optimization of Bionic Tree-like Channel Liquid Cold Plate on Lithium-ion Battery

LIU Huan-huan1,2, REN Xiao-long2, ZHANG Ze-bin1*

(1.School of Mechanical and Power Engineering, Zhengzhou University, Zhengzhou 450001, China; 2.Guangdong Shunde Innovative Design Institute, Guangdong Foshan 528300, China)

The work aims to achieve the goal that the maximum temperature and temperature difference of lithium-ion battery are within the suitable temperature range when the discharge rate is 5and the ambient temperature is 298.15 K. Bionic tree-like liquid cold plates with channels were inserted on both sides of the lithium-ion battery to establish the liquid cold plate heat dissipation model of the battery module. Firstly, the accuracy of the numerical model of liquid cold plate was verified, and then the effect of the angle of branch runner, the inlet velocity, the width of branch runner and the thickness of liquid cold plate on the cooling performance of cold plate was studied by orthogonal experimental design. The importance of the four factors was ranked by range analysis, and the optimal parameter combination was determined. The maximum temperature of the optimized scheme battery was 302.4 K and the temperature difference was 3.4 K, which met the design requirements. However, there was a risk of coolant leakage during the flow cycle. The tree-like liquid cold plate can meet the requirements of battery thermal management performance. The optimized battery module has a safe range of surface temperature, and the battery temperature uniformity is significantly improved. Cold plate materials with high thermal conductivity should be selected, and attention should be paid to the structural design of cold plate to prevent coolant leakage.

lithium-ion battery; liquid cold plate; bionic tree-like structure; orthogonal experiment

TK172；TM912.9

A

1001-3563(2023)19-0273-10

10.19554/j.cnki.1001-3563.2023.19.035

2023-04-14

國(guó)家自然科學(xué)基金面上項(xiàng)目（12272354）

責(zé)任編輯：彭颋