汽车工程学报

参数	液冷板（铝）	冷却液（水）
比热容c_p/［J/（kg·K）］	903	4 181.72
导热系数k/［W/（m·K）］	237	0.62
动力黏度μ/（Pa·s）		1.003E-3
密度ρ/（kg/m³）	2 702	997.561

参数	液冷板（铝）	冷却液（水）
比热容c_p/［J/（kg·K）］	903	4 181.72
导热系数k/［W/（m·K）］	237	0.62
动力黏度μ/（Pa·s）		1.003E-3
密度ρ/（kg/m³）	2 702	997.561

边界条件	参数值
入口温度/℃	25
出口压力/Pa	0
入口质量流量/（g/s）	1～5
热通量/（W/m²）	826

边界条件	参数值
入口温度/℃	25
出口压力/Pa	0
入口质量流量/（g/s）	1～5
热通量/（W/m²）	826

锂离子电池液冷散热通道拓扑优化设计研究

PDF下载

陈瑶 ¹ , 詹森 ¹ , 黄文姣 ² , 刘聪 ¹ , 许重阳 ¹

汽车工程学报 | 绿色低碳技术专栏 2025,15(4): 567-577

收起

汽车工程学报 | 绿色低碳技术专栏 2025, 15(4): 567-577

锂离子电池液冷散热通道拓扑优化设计研究

全屏

陈瑶¹, 詹森¹, 黄文姣², 刘聪¹, 许重阳¹

作者信息

¹ 重庆交通大学机电与车辆工程学院，重庆 400074

² 中国汽车工程研究院股份有限公司，重庆 401122

陈瑶（2000-），男，重庆市人，硕士研究生，主要研究方向为电池热管理系统设计与控制。 E-mail：cycy5369@163.com

通讯作者:

詹森（1986-），男，湖北十堰人，博士，讲师，主要研究方向为新能源汽车动力传动系统及控制。 E-mail：fishjason@163.com

Study on the Topology Optimization Design of Liquid Cooling Heat Dissipation Channels for Lithium-Ion Batteries

Yao CHEN¹, Sen ZHAN¹, Wenjiao HUANG², Cong LIU¹, Chongyang XU¹

Affiliations

¹ School of Mechatronics & Vehicle Engineering，Chongqing Jiaotong University，Chongqing 400074，China

² China Automotive Engineering Research Institute Co.，Ltd.，Chongqing 401122，China

出版时间: 2025-07-20 doi: 10.3969/j.issn.2095‒1469.2025.04.13

文章导航

摘要

收起

针对传统通道型电池液冷板压降高和均温性差的问题，采用多目标拓扑优化方法对液冷板进行优化设计。通过试验建立电池产热模型，基于变密度法构建液冷板拓扑优化模型，研究了不同进出口位置布置形式对拓扑优化后液冷板的性能影响，从中优选出性能最佳的液冷板，并将其与传统直通道液冷板对比分析。结果表明，不同进出口位置布置形式下得到的拓扑通道在温度和压降性能方面差异较大，当进出口布置在沿液冷板长边中心对称线时，拓扑优化后的液冷板具有最好的综合性能，其与直通道液冷板相比具有更强的流动换热性能，在入口流量为5 g/s时，最高温度、温度标准差和压降分别下降了1.38%、22.35%和28.36%。这种新型液冷板可为后续电池热管理系统散热设计提供新思路。

关键词

电池热管理系统 / 液冷板 / 拓扑优化 / 进出口位置

Abstract

收起

To address the issues of high pressure drop and poor temperature uniformity in traditional channel-type battery liquid cooling plates, a multi-objective topology optimization method was employed for the design optimization of the liquid cooling plate. An experimental model of battery heat generation was established, and a topology optimization model of the liquid cooling plate was constructed based on the variable density method. The impact of different inlet and outlet arrangements on the performance of the optimized liquid cooling plate was investigated, and the best-performing liquid cooling plate was selected and compared with the traditional straight-channel liquid cooling plate. The results indicate that the topology channels obtained under different inlet and outlet arrangements exhibit significant differences in temperature and pressure drop performance. When the inlet and outlet are arranged along the central symmetry line of the long edge of the liquid cooling plate, the topology-optimized liquid cooling plate demonstrates the best overall performance. Compared to the straight-channel liquid cooling plate, it exhibits stronger flow and heat transfer performance, with the maximum temperature, temperature standard deviation, and pressure drop reduced by 1.38%, 22.35%, and 28.36%, respectively, at an inlet flow rate of 5 g/s. This novel liquid cooling plate can provide new insights for the thermal design of future battery thermal management system.

Key words

battery thermal management system / liquid cooling plate / topology optimization / inlet and outlet position

引用本文

陈瑶, 詹森, 黄文姣, 刘聪, 许重阳. 锂离子电池液冷散热通道拓扑优化设计研究. 汽车工程学报, 2025 , 15 (4) : 567 -577 . DOI: 10.3969/j.issn.2095‒1469.2025.04.13

Yao CHEN, Sen ZHAN, Wenjiao HUANG, Cong LIU, Chongyang XU. Study on the Topology Optimization Design of Liquid Cooling Heat Dissipation Channels for Lithium-Ion Batteries[J]. Chinese Journal of Automotive Engineering, 2025 , 15 (4) : 567 -577 . DOI: 10.3969/j.issn.2095‒1469.2025.04.13

正文

收起

以电池为主要动力来源的电动汽车在解决能源与环境问题方面将发挥关键作用^[1]。锂离子电池以其高能量密度、低自放电率和较长的使用寿命等优点，成为电动汽车动力电池的首选^[2]。但是锂离子电池在充放电过程中焦耳热和电化学反应导致大量热量聚集，使锂离子电池容量和寿命衰减，极端情况下会发生热失控等安全性问题^[3-4]。因此，设计合理的电池热管理系统（BTMS）以确保锂离子电池工作在合理的温度范围内尤为重要，研究表明，锂离子电池最佳的工作温度范围为20～40 ℃，电池组各单元的最大温差不超过5 ℃^[5-6]。

根据冷却方式不同，BTMS大致可以分为风冷、液冷和相变材料冷却^[7]。其中，液冷BTMS因其优异的冷却性能、适中的制造和运行维护成本，目前已经被各大车企推广应用^[8]。液冷BTMS冷却性能的关键在于配置合适的液冷板。合理的液冷板通道结构设计可以改善电池组工作环境。BAO Yunhao等^[9]研究了一种超薄宽直流道液冷板，与蛇形、分叉和U形通道冷板相比，超薄宽直流道液冷板具有更低的最高温度、温差和运行功耗。陈雅等^[10]针对蛇形流道压降高、功耗大的问题设计了一种新型的二次流蛇形液冷板，结果表明，传统蛇形液冷板在增加了二次流结构后，液冷板进出口压降下降了90.69%。LIU Feifei等^[11]设计了一种仿生叶脉通道液冷板，并通过遗传算法对通道宽度、角度、数量和入口流量进行了优化，在3 C放电倍率下，优化后的液冷板可以将电池最高温度控制在33.34 ℃以内。此外，也有学者采用先进的拓扑优化方法对液冷板进行设计。陈帆^[12]讨论了换热量、出口流体焓值和固体域温度这3种不同目标函数对拓扑优化液冷板设计结果的影响，结果发现出口流体焓值为目标函数的液冷板的流动特性换热性能最佳。ZHANG Kezheng等^[13]针对多通道蛇形冷板中冷却剂流动路径长、压降高的问题，设计了一种多约束拓扑优化液冷板，相比于多通道蛇形冷板可使综合性能指标提高2～3倍。

传统液冷板通道结构设计很大程度上依赖于经验，手动优化结构参数提升液冷板性能具有较大局限性，而拓扑优化设计根据目标函数与优化算法自动求解符合特定性能要求的最佳通道结构，能设计出性能更优的液冷板通道结构。此外，不同进出口位置下拓扑优化后液冷板的通道结构不同，从而进一步影响液冷板温度和压降等性能。针对不同进出口位置对电池液冷板拓扑优化结果的影响进行的系统性研究较少，不能有效发挥拓扑优化的潜力。对此，本文以平均温度和流体耗散功为双目标函数，拓扑优化设计得到了8种典型位置下的二维通道拓扑形态，通过拉伸构建了三维模型，对所有模型进行了三维数值模拟，优选出了进出口位置性能最佳的液冷板，并和传统直通道液冷板进行了流动换热性能对比，优化设计结果能有效提高电池热管理系统性能。本文的研究框架如图1所示。

1 模型与方法

收起

1.1 电池产热模型

由Bernardi生热模型可知，电池放电过程中的产热与吸热分别可表示为^[14]：

Q 1 = I 2 R j + I T (∂ U o c v / ∂ T)

。（1）

Q 2 = m c p (d T / d t)

。（2）

式中：R_j为焦耳内阻；T为温度；U_ocv为开路电压；I为电流；m为质量；c_p为比热容。

在试验环境为绝热条件下时，电池的产热量等于吸热量，由式（1）和式（2）可得：

1 I d T d t = R j m c p I + 1 m c p T ∂ U o c v ∂ T

。（3）

针对同一个电池在较短时间的放电过程中，dT/dt和

T (∂ U o c v / ∂ T)

可以视为常数，式（3）可视为电池温度随时间变化的线性关系式，如图2所示，电池温升试验平台对容量为30 Ah，内阻为2 mΩ，质量为0.55 kg的软包锂离子动力电池在不同放电电流下的温升变化进行了测量。通过不同放电电流下温升试验数据建立了（1/I）

×

（dT/dt）与I的关系曲线，如图3所示。则式（3）可以表示为：

d T d t = 4.21 × 10 - 6 I 2 + 1.49 × 10 - 4 I

。（4）

式中：R_j/mc_p为拟合公式斜率。

因此，可计算电池比热容c_p为864 J/（kg·K），综合式（1）～（4）可得电池的产热公式为：

Q = 2 × 10 - 3 I 2 + 7.08 × 10 - 2 I

。（5）

最终可以得到电池的单位体积产热率为：

q = Q / V = 8.13 I 2 + 287.78 I

。（6）

由式（6）可知，电池在不同放电倍率下的产热量。

1.2 液冷板拓扑优化设计

1.2.1 设计域描述

液冷板在电池组中的布置方式如图4所示，为了简化模型和减少计算量，电池产生的热量被假设为均匀通过冷板，冷却液在冷板内循环流动，冷板表面与电池表面发生热交换，从而带走电池工作时所产生的热量。

本文基于变密度法对电池液冷板进行二维拓扑优化设计。变密度拓扑优化方法的基本原理是将设计域离散为多个设计微单元，赋予设计单元一个0～1的伪密度参数来描述该处是否有材料存在，当改变各个单元的密度值时可以表征出所需的拓扑结构，同时当密度为非0和1的状态时会出现不清晰的结构灰度单元，可通过密度过滤和投影来获得清晰的结构拓扑形态，最终优化问题是在目标函数下根据优化算法与约束条件寻求设计域内最佳的材料密度分布问题。拓扑优化设计域描述如图5所示，优化设计域的大小为200 mm×130 mm，为避免生成的流道与设计域边缘相交，在设计域的周围设置了2 mm的不可优化区域，设计域的总大小为134 mm×204 mm，其与电池表面大小一致，以此获得最大的热交换量。对于流体拓扑优化问题，设计域内被视为一个多孔介质区域，在进行拓扑优化时，多孔介质区域被离散为多个变量计算单元，当设计变量γ为1时，代表该区域为液体，当γ为0时，代表该区域为固体。

1.2.2 控制方程

在电池液冷板内，流体与冷板壁发生热交换，属于流固共轭传热问题，在稳态条件下满足流体力学三大控制方程。拓扑优化设计域中能量平衡方程可表示为：

ρ (γ) C p (γ) (u × ∇ T) = ∇ × (k (γ) ∇ T) + Q

。（7）

式中：ρ为密度；C_p为比热容；u为流速；T为温度；k为导热系数；Q为体积产热率。

方程中材料密度、比热容和导热系数均为设计变量γ的函数，通过函数形式将控制方程与设计变量进行耦合。

设计域连续性方程：

∇ × u = 0

。（8）

动量守恒方程：

ρ (γ) (u × ∇) u = - ∇ P + μ ∇ × (∇ u + ∇ u T) + F

。（9）

式中：P为压力；μ为动力黏度；F为体积力项，在动量方程中加入此项用于模拟多孔介质中流体流动所受到的摩擦阻力。

由Brinkman惩罚模型可知，多孔介质中流体流动的阻力与u相关，同时引入反渗透率

α

来控制多孔介质的局部渗透率，体积力可以表示为：

F = - α (γ) u

。（10）

当

α

为0时流体流动无障碍，而当其趋近于无穷大时，摩擦阻力无穷大，因此，改变

α

的值可以控制流体与固体的分布情况，同时

α

也是设计变量γ的函数。

1.2.3 材料插值

通过插值函数建立起反渗透率与设计变量之间的关系^[15]：

α (γ) = α m i n + (α m a x - α m i n) × q (1 - γ) q + γ

。（11）

α m a x = μ D a × L 2

。（12）

式中：q为惩罚因子，用于调节

α

的变化趋势，取值为0.05；

α

_max和

α

_min为仅存在固相和液相时的反渗透率；Da为达西数，用于描述多孔介质中黏性力和摩擦力的关系，其值越大则多孔介质渗透率越大，本文取值为1×10^-4；L为模型特征长度。

同理可建立材料导热率k、密度ρ和比热容c的插值函数，从而实现对固体与流体之间物理属性的转化：

k (γ) = k s + (k f - k s) × γ (1 + q k) q k + γ

。（13）

c (γ) = c s + (c f - c s) × γ (1 + q c) q c + γ

。（14）

ρ (γ) = ρ s + (ρ f - ρ s) × γ (1 + q ρ) q ρ + γ

。（15）

式中：s和f分别代表固体和流体。

1.2.4 过滤和投影

由于伪密度参数为人为假设，在拓扑优化后会存在密度为0～1的中间状态，这将导致结构不清晰、数值不稳和出现棋盘格分布现象等。为了消除网格依赖和棋盘格问题，采用了Helmholtz密度滤波器对设计变量进行过滤处理，其表达式为：

- r 2 ∇ 2 γ ˜ + γ ˜ = γ

。（16）

式中：

γ ˜

为过滤后的设计变量；r为过滤半径。

同时使用正切投影来消除滤波后的灰度单元以获得清晰的拓扑形态：

γ^= t a n h β (γ ˜ - γ β) + t a n h (β γ β) t a n h β (1 - γ β) + t a n h (β γ β)

。（17）

式中：

γ^

为投影后的设计变量；β和

γ β

为投影斜率和投影点，本文取值为16和0.5。

1.2.5 目标函数及约束

对于电池液冷板的设计要求一般包括散热性能指标和功耗性能指标，是一个多目标优化问题。本文考虑散热性能指标为使设计域的整体平均温度最小化：

Φ = ∫ Ω T d Ω / ∫ Ω 1 d Ω

。（18）

功耗目标函数由进出口总压降表示，其进一步可表示为设计域黏性耗散功：

Ψ = ∫ Ω 12 μ ∑ i, j ∂ u i ∂ x j + ∂ u j ∂ x i 2 + ∑ i α (γ) u i 2 d Ω

。（19）

总的目标函数可以表示为两个性能指标的加权和：

J = ω Φ Φ 0 + (1 - ω) Ψ Ψ 0

。（20）

式中：

Ψ 0

和

Φ 0

为两个目标函数的初始值，同时除以初始值以此实现目标函数归一化；ω为权重系数，本文取值为0.8。

为了使拓扑优化的结果具有实际工程意义，对设计域中流域体积分数进行了限制：

∫ Ω γ d Ω / ∫ Ω d Ω ≤ V f

。（21）

式中：V_f为流域占整个设计域的比例，本文取值为0.5。

1.3 拓扑优化设计结果

本文基于COMSOL平台实现拓扑优化算法的编译，采用SNOPT算法进行迭代求解，模型迭代终止条件为目标函数误差值低于1e^-6。拓扑优化模型的目标函数和流域体积分数的迭代过程如图6所示，目标函数值和体积分数在60步后逐渐趋于稳定，150步时达到收敛条件，停止迭代。

进出口在液冷板上的布置位置对液冷板的最终拓扑结果和综合性能具有较大的影响。对于进出口位置的选取原则，通过总结大量液冷板设计的参考文献发现，进出口位置设置一般位于冷板边长的中部或边缘位置，而进口与出口的相对位置一般为对向、对角或位于同一条边上，同时进口与出口不能相距太近以避免换热不够充分。对此，本文设计了一共8种不同进出口位置下的拓扑优化液冷板，包括沿着设计域的短边和长边布置的4种直进直出的对向位置，2种平行对角位置和2种进出口位于相同侧的位置，分别命名为Case1至Case8。8种不同进出口位置下的拓扑优化结果如图7所示。不同进出口位置下的流道拓扑形态均表现出由较宽的主流道和较细的次级流道构成，总体上呈现出自然界中的植物根系分布形态，这种结构将有利于强化散热和减少流体压降。本文将在后续三维仿真中具体分析这些不同进出口位置拓扑优化液冷板的性能。

2 三维流热耦合数值模型与验证

收起

2.1 边界条件

为了研究不同进出口位置下拓扑优化液冷板的三维热力学性能，对所有拓扑结果进行了纵向拉伸以获得三维数值模拟模型，拉伸的总高度为3 mm，流道高度为2 mm。在CATIA中将二维拓扑结果进行拉伸后，再导入流体仿真软件STAR-CCM+中进行数值模拟研究。液冷板和冷却液仿真的物理属性见表1。此外雷诺数可以衡量流体流动状态，其表示为：

R e = ρ f D h u i n μ

。（22）

式中：D_h为通道水力直径；u_in为入口流速；μ为流体动力黏度。

雷诺数小于2 300时可认为流体处于层流状态，本文所采用的入口最大质量流量为5 g/s，对应的雷诺数为419.18，其值小于2 300，因此，本文数值模拟中均采用层流模型。

为了简化仿真，作出如下假设：

1）冷却水在流动传热过程中是稳定的；

2）冷却水为不可压缩流动流体，液冷板的热物理性质不发生改变；

3）入口为质量流量入口，出口为压力出口；

4）忽略重力、黏性耗散和接触热阻的影响；

5）流体与固体交界面无滑移；

6）液冷板除发热面外的表面，其余表面为绝热。

由表2可知具体边界条件，其中热通量由电池热模型中式（6）计算出电池3 C放电倍率下的体积热再乘以电池厚度所得。

2.2 控制方程

三维数值模拟中质量守恒方程为：

∂ ρ f ∂ t + ∇ × ρ f u f = 0

。（23）

动量守恒方程为：

∂ ρ f u f ∂ t + ∇ × ρ f u f u f = - ∇ P + μ ∇ 2 u f

。（24）

能量方程为：

ρ f C p, f ∂ T ∂ t + ∇ × (u f T) = k f ∇ 2 T (F l u i d), ρ s C p, s ∂ T ∂ t = k s ∇ 2 T (S o l i d) 。

（25）

式中：u_f为冷却液流速；ρ_f和ρ_s分别为流体和固体密度；k_f和k_s分别为流体和固体导热率；C_p_，f和C_p_，s分别为流体和固体比热容。

电池热管理系统的设计要求是将电池最高温度和温差控制在合理的范围内，同时降低系统功耗。液冷板温度性能可以由表面最高温度T_max和表征温度均匀性的温度标准差T_σ定义，T_σ表示为：

T σ = ∫ A c T - T a v g 2 d A c ∫ A c d A c

。（26）

式中：T_avg为液冷板表面平均温度；A_c为液冷板受热面的面积。

平均换热系数定义为：

h a v g = Q A c (T a v g - T i n + T o u t 2)

。（27）

式中：T_in和T_out分别为冷却液的进口和出口温度。

努塞尔系数定义为：

N u = h a v g D h / k f

。（28）

流动阻力系数定义为：

f = Δ P ρ f u i n 2

。（29）

PEC指标用于评估液冷板换热与功耗的综合性能：

P E C = N u / N u r e f Δ P / Δ P r e f 1 / 3

。（30）

式中：Nu_ref和ΔP_ref分别为参照模型的努塞尔系数和进出口压降。

2.3 网格独立性分析

本文采用多面体网格对三维冷板模型进行了网格划分，同时为了保证计算的准确性，在流固边界层细化了网格数量，网格划分如图8所示。Case3的表面平均温度和压降随网格划分数量的变化如图9所示，网格范围在297 306～1 482 958个之间，当网格数达到778 628个时，平均温度和压降的变化范围均在0.2%以内，为了保证计算精度的同时节约计算时间，本文采用了网格数为778 628个进行数值模拟。

3 结果与讨论

收起

3.1 不同进出口位置的影响分析

3.1.1 温度分布

图10和图11分别为在入口质量流量为2 g/s时不同进出口位置拓扑优化液冷板的表面温度云图，以及最高温度、平均温度和温度标准差的具体数值对比。由表面温度云图可知，不同进出口位置布置下对拓扑优化后液冷板表面温度分布影响较大。其中Case2表面温度最高，达到了36.4 ℃，Case4、Case7和Case8均出现出一些极热区域，而Case3和Case6的表面温度梯度较小。由图11可知，Case3具有最低的平均温度为29.34 ℃，第2低的最高温度32.06 ℃和温度标准偏差1.5 ℃。而Case6具有最低的温度标准偏差1.47 ℃和最高温度31.70 ℃，第2低的平均温度29.42 ℃。以上结果表明，当进出口位置位于沿液冷板长边中心对称线和平行对角时，拓扑优化后液冷板表面具有较好的温度性能，采用这两种布置方式所形成的流道布局更均匀，冷却液在液冷板内充分均匀流动换热，从而降低了最高温度，提升了温度均匀性。

3.1.2 压力分布

图12和图13分别为不同进出口位置拓扑优化液冷板的中部截面压力云图和进出口压降数值对比。由压力云图可知，Case1和Case3的高压区域较少，而其他进出口类型液冷板的高压区域较多，其中Case8中的高压区域最多，在Case2和Case8这种进出口布置下，冷却液在液冷板内进行了较长的流动循环，导致较高的压力损失。由图13可知，Case1具有最小的压降为14.04 Pa，而Case3次之，压降为17.19 Pa。由于这两种进出口位置位于液冷板中心对称线上，其拓扑优化后形成的流道对称，且流动距离较短，使流体压降较低。而对于平行对角布置下的Case6，尽管其在温度性能方面表现优异，但由于这种布置下进出口的距离较远，导致拓扑优化后形成的流道路径较长，进而导致其压降较高，这将增加液冷板的能耗。因此，进出口位置的布置应当避免拓扑优化后形成较长的流动路径。

图14展示了以Case1为参照模型的所有类型液冷板的PEC值。由图可知，Case3具有最高的PEC值为1.32，这意味着Case3的综合性能最好，而尽管Case1压降低，但温度性能略差，Case5和Case6的温度性能虽然优于Case1，但压降性能略差。综合来看，进出口沿液冷板长边中心对称线布置下的Case3，在权衡温度和压降性能指标后是用于电池热管理系统的最佳选择。

3.2 流动与换热分析

3.2.1 流动特性

为了更深入地研究最优设计结果Case3的流动换热性能，选取1～5 g/s范围，间距为0.5 g/s的质量流量进行了仿真研究。同时为了对比拓扑优化液冷板与传统直通道液冷板的性能差异，建立了一块目前工业上常用的直通道电池液冷板。为了性能对比的公平性，其外形总尺寸大小与Case3相同，且流域体积分数同样保持为0.5，将其命名为RCP，如图15所示。

由图16a可知，Case3与RCP的压降随着质量流量的变化趋势。两种类型液冷板的压降都与质量流量成正相关。Case3的压降低于RCP，且随着质量流量的增加，差距变得明显，在5 g/s时，压降相比传统直通道液冷板降低了19.12 Pa（28.36%）。这是由于拓扑优化目标函数中有表征流动阻力的耗散功函数存在，在拓扑优化过程中会尽可能形成平滑的流道曲线以最小化流动阻力，其所形成的流道由截面积较宽和较细的分支流道组合而成，减少了速度梯度带来的能量损失，有效分散了流体压力。压降与液冷板泵送功率成正比，因此，拓扑优化液冷板可以有效减少泵功率，从而降低电池热管理系统的能耗。

由图16b可知，2种类型液冷板的流动阻力系数随质量流量的变化趋势。Case3的流动阻力系数小于RCP。传统直通道液冷板在拐角处为直角，这将导致较大的流动阻力，而拓扑优化通道在拐角处平滑过渡，减少了涡流的产生，且拓扑优化能降低流道壁面粗糙度，使拓扑优化形成的流道具有更小的流动阻力系数。因此，在流动性能方面，拓扑优化液冷板要优于传统直通道液冷板。

3.2.2 换热特性

由图17a可知，2种类型液冷板的最高温度随质量流量的变化趋势。Case3的最高温度低于RCP，拓扑优化目标函数中有温度性能指标，所形成的流道会尽可能地分布在液冷板的各个角落，使冷却液可以对所有受热位置进行冷却，从而减少了温度热点的出现，降低了最高温度。在流量较低时，两种液冷板的最高温度差距不明显，而随着质量流量的增加，拓扑优化液冷板的温度性能优势逐渐增加，在5 g/s时，最高温度相比于传统直通道液冷板降低了0.4 ℃（1.38%）。

由图17b可知，2种类型液冷板的温度标准差随质量流量的变化趋势。与最高温度类似，Case3的温度标准差也低于RCP，拓扑优化通过优化流道布局，寻找了给定热负荷下的最优流道分布，针对热源强度及分布，不断通过迭代寻优调整流道布局，最终形成了对称且均匀分布的流道，能有效提升温度均匀性，降低温度标准差。这种优势同样随着流量的增加而增加，在5 g/s时，温度标准差降低了0.19 ℃（22.35%）。这表明了拓扑优化液冷板能有效提升电池组的温度均匀性，延长电池寿命。

努塞尔系数可深入表征液冷板的对流传热能力。由图17c可知，2种类型液冷板的努塞尔系数随质量流量的变化趋势。Case3的努塞尔系数大于RCP，优化后的流道结构使流体流动更加均匀，改善了直通道存在流动死角或流速不均的问题，同时较复杂的流道形状增加了冷却液的接触面积，提高了对流换热的效率。此外，增加质量流量不能无限地提高液冷板的换热能力，在初始时努塞尔系数增幅较大，而随着流量进一步增加，趋势变缓，继续增加流量将会导致更高的压降。

入口质量流量为5 g/s时2种类型液冷板的温度、中部截面压力与速度云图对比如图18所示，Case3相比于RCP的高温区域要少，其表面温度分布更均匀。Case3的整体压力分布也比RCP更低，RCP在进口及转角处受到的阻力较大，而拓扑优化所形成的流道在拐角处近似于流线型，有效地降低了压降，从而减少泵功率。此外，均匀的流速分布能降低液冷板表面温度梯度，从而降低电池组温差，Case3的整体流速分布更均匀，而RCP中部流速明显高于两侧。综合前面的分析，Case3的流动与换热性能要远优于RCP，能有效提升电池热管理系统性能并降低系统功耗，证明了本文设计方法的有效性。

4 结论

收起

本文为了设计传热能力优异和流动阻力较低的电池液冷板，采用了双目标拓扑优化方法对液冷板进行了优化设计，研究了不同进出口位置布置下对拓扑优化后液冷板的性能影响，并将最优设计结果与传统直通道液冷板进行了对比。主要结论如下。

1）不同进出口位置布置下拓扑优化后的电池液冷板在温度和压降性能方面差异较大。进出口沿液冷板长边中心对称线布置的Case3的平均温度最低，而平行对角布置的Case6的最高温度和温度标准差最低。

2）当进出口布置在沿液冷板中心对称线时可以降低压降，而其余位置则导致较高的压降。进出口沿液冷板短边中心对称线布置的Case1具有最低的压降，Case3次之。

3）在权衡温度和压降性能指标下，Case3具有最好的综合性能，与传统直通道液冷板相比，其具有更强流动换热性能，在入口流量为5 g/s时，最高温度、温度标准差和压降分别下降了1.38%、22.35%和28.36%，有效提升了电池热管理系统性能。

基金

收起

重庆市教委青年项目(KJQN202000734)
机械传动国家重点实验室开放基金项目(SKLMT-MSKFKT-202112)
重庆市技术创新与应用发展重点项目(cstc2020jscx-dxwtBX0025)

参考文献

收起

文献

收起

参考文献引证文献

排序方式：

[1]

陈昊鹏，张天时，林晓东，等.应用于动力电池热管理的重力型热管仿真研究[J].汽车工程学报，2023，13（1）：54-63.

CHEN

Haopeng

， ZHANG

Tianshi

， LIN

Xiaodong

，et al.Simulation Study of Gravity Heat Pipe for Power Battery Thermal Management[J].Chinese Journal of Automotive Engineering，2023，13 （1）：54-63.（in Chinese）

[2]

李嘉鑫，李鹏钊，王苗，等.锂离子电池热管理技术研究进展[J].过程工程学报，2023，23（8）：1102-1117.

Jiaxin

， LI

Pengzhao

， WANG

Miao

，et al.Research Progress in Thermal Management Technology of Lithium-Ion Batteries[J].Journal of Process Engineering，2023，23 （8）：1102-1117.（in Chinese）

[3]

曾少鸿，吴伟雄，刘吉臻，等.锂离子电池浸没式冷却技术研究综述[J].储能科学与技术，2023，12（9）：2888-2903.

ZENG

Shaohong

， WU

Weixiong

， LIU

Jizhen

，et al.Review of Research on Immersion Cooling Technology for Lithium-Ion Batteries[J].Energy Storage Science and Technology，2023，12 （9）：2888-2903.（in Chinese）

[4]

李顶根，邹时波，徐鹏，等.不同热管理方案下锂离子电池模组温度特性分析[J].汽车工程学报，2020，10（2）：98-106.

Dinggen

， ZOU

Shibo

， XU

Peng

，et al.Analysis of Temperature Characteristics of Lithium-Ion Battery Modules Under Different Thermal Management Schemes [J].Chinese Journal of Automotive Engineering，2020，10 （2）：98-106.（in Chinese）

[5]

ZHAO

Gang

， WANG

Xiaolin

， NEGNEVITSKY

，et al.An Up-to-Date Review on the Design Improvement and Optimization of the Liquid-Cooling Battery Thermal Management System for Electric Vehicles[J].Applied Thermal Engineering，2023，219：119626.1-119626.21.

[6]

WANG

Jiao

， CHEN

Fan

， SHAO

Zhenyu

，et al.Study of the Influence of Objective Functions on the Topology Optimization Design of Battery Cold Plate[J].Applied Thermal Engineering，2023，226：120326.1-120326.12.

[7]

王毅军，周舟，李军.锂离子电池热管理技术[J].新能源进展，2023，11（1）：54-62.

WANG

Yijun

， ZHOU

Zhou

， LI

Jun

.Thermal Management Techniques for Lithium-Ion Batteries[J].Advances in New and Renewable Energy，2023，11（1）：54-62.（in Chinese）

[8]

SHENG

Lei

， ZHANG

Hengyun

， SU

Lin

，et al.Effect Analysis on Thermal Profile Management of a Cylindrical Lithium-Ion Battery Utilizing a Cellular Liquid Cooling Jacket [J].Energy，2021，220：119725.1-119725.12.

[9]

BAO

Yunhao

， SHAO

Shuangquan

. Numerical Study on Ultrathin Wide Straight Flow Channel Cold Plate for Li-Ion Battery Thermal Management [J].Journal of Energy Storage，2023，64：107263.1-107263.16.

[10]

陈雅，范立云，李晶雪，等.二次流蛇形通道锂离子电池散热性能[J].储能科学与技术，2023，12（6）：1880-1889.

CHEN

， FAN

Liyun

， LI

Jingxue

，et al.Heat Dissipation Performance of Secondary Flow Serpentine Channel Lithium-Ion Batteries [J].Energy Storage Science and Technology，2023，12 （6）：1880-1889.（in Chinese）

[11]

LIU

Feifei

， CHEN

Yangyang

， QIN

，et al.Optimal Design of Liquid Cooling Structure with Bionic Leaf Vein Branch Channel for Power Battery[J].Applied Thermal Engineering，2023，218：119283.1-119283.16.

[12]

陈帆.电池热管理系统中冷板流道的拓扑优化设计研究[D]. 镇江：江苏科技大学，2023.

CHEN

Fan

.Research on Topology Optimization Design of Cold Plate Flow Channel in Battery Thermal Manage-ment System[D]. Zhenjiang：Jiangsu University of Science and Technology，2023.（in Chinese）

[13]

ZHANG

Kezheng

， LI

Yang

， CHANG

S M

，et al. Hydraulic and Thermal Performance Enhancement for the Cold Plate Using Topology Optimization[J].Applied Thermal Engineering，2024，236：121829.1-121829.16.

[14]

DENG

Tao

， RAN

Yan

， YIN

Yanli

，et al.Multi-Objective Optimization Design of Thermal Management System for Lithium-Ion Battery Pack Based on Non-Dominated Sorting Genetic Algorithm II[J].Applied Thermal Engineering，2020，164：114394.1-114394.11.

[15]

XIA

Yang

， CHEN

， LUO

Jiwang

，et al.Numerical Investigation of Microchannel Heat Sinks with Different Inlets and Outlets Based on Topology Optimization [J].Applied Energy，2023，330：120335.1-120335.16.

2025年第15卷第4期

PDF下载

539

224

引用本文

BibTeX

文章信息

doi: 10.3969/j.issn.2095‒1469.2025.04.13

接收时间：2024-05-27
首发时间：2025-09-10
出版时间：2025-07-20

补充材料

相关文章

2种不同金属材料的力学参数

科 Family	属数 Number of genus	种数 Number of species	占总种数比例 Percentage of total species (%)	属 Genus	种数 Number of species	占总种数比例 Percentage of total species (%)
鹅膏菌科Amanitaceae	2	11	5.26	鹅膏菌属 Amanita	10	4.78
小菇科 Mycenaceae	2	12	5.74	丝盖伞属 Inocybe	5	2.39
多孔菌科 Polyporaceae	8	14	6.70	蜡蘑属 Laccaria	5	2.39
红菇科 Russulaceae	3	23	11.00	小皮伞属 Marasmius	6	2.87
				小菇属 Mycena	11	5.26
				光柄菇属 Pluteus	5	2.39
				红菇属 Russula	17	8.13
				栓菌属 Trametes	5	2.39

关闭全屏

BibTeX
EndNote
RefWorks
TxT