基于电化学-热耦合模型，以4节18650锂离子电池为研究对象，分析考虑温度效应时的对流换热系数对锂电池热扩散的影响及其程度。首先基于传热理论中的流体横掠顺排管束平均表面换热系数计算方法，计算得到不同温度和流速下锂离子电池表面对流换热系数，通过曲线拟合得到空气流速分别为0.05、0.1、0.2和0.3m/s时对流换热系数与温度的函数关系，得出对流换热系数与温度不完全呈线性变化；其次基于以上函数关系，通过数值模拟分析了考虑温度效应时的对流换热系数对锂电池热扩散的影响。

(相关资料图)

锂离子电池因具有高能量密度且循环寿命较长等优点给生活带来极大便利， 但近几年新能源汽车因锂离子电池屡次发生火灾甚至爆炸事件，众多研究人员对锂离子电池进行实验和仿真研究并分析热量来源与锂离子电池热失控的热特性。 Wang等利用加速绝热量热仪和电池循环仪研究锂电池的热响应，结果表明热量主要来源于放电和热失控，电池容量越大存在的风险越大；锂电池在高温的环境下工作，热稳定性明显降低。为降低锂电池热失控的发生概率，科研人员研究了锂电池热管理系统（BTMS），使锂电池的工作温度始终稳定在最佳温度范围内。刘晓东分别建立了集总模型和电-热耦合模型研究锂离子电池热管理系统，模拟结果表明集总热模型的优点在于模拟温度精度较高，而电-热耦合模型在模拟外部电压方面有优势。刘晓东基于集总模型设计冷板通道结构，选择效果较好的蛇形冷板通道并对其进行优化。王晓慧等建立强制风冷工况和自然冷却工况下的仿真模型，发现强制风冷下1C放电结束时温升约12.5℃，比自然冷却工况降低3℃，模型结果表明使用环境温度、0.5m/s的冷却空气时就可满足散热的需求并且单体电池温度能达到较好的一致性。大多研究在假设锂电池对流换热系数定值的前提下改变锂电池充放电倍率、不同荷电状态（state of charge，SOC）等参数研究锂电池热特性参数，而很少考虑对流换热系数的具体计算方法以及在受温度影响时对锂电池放电过程中温度场的影响。

对流换热系数和导热系数均是研究锂电池热力学仿真不可或缺的热物性参数。 冯旭宁等将锂电池简化为各向异性的长方体，建立锂离子电池热模型，利用“三线共点”图解法得出垂直于锂电池内部层叠方向的导热系数小于平行方向的导热系数。王翔等利用变对流换热系数（variable convective heat transfer coefficient，VCHTC）建立锂电池的生热优化模型研究不同环境温度和放电倍率下单体锂电池的温升特性，结果表明该模型结果与恒定对流传热系数模型在不同工况下的平均绝对百分误差分别为1.1%和6.9%。黄文才在高温热失控仿真模拟中通过设置不同对流换热系数研究该参数对锂电池热行为的影响，结果显示传热系数越大，锂电池发生热失控所需时间越短。以上研究表明了不同换热系数对锂电池散热的影响，但未定量给出对流换热系数考虑温度时对锂电池热扩散影响程度。选取某大型电池包中的4节18650锂电池作为研究对象，基于传热理论中的流体横掠顺排管束平均表面换热系数计算方法，模拟在强制风冷状况下锂电池4C放电的过程，研究对流换热系数考虑温度时对锂电池热扩散影响及其程度。

锂离子电池工作及产热原理

01锂离子电池工作原理

锂离子电池充放电过程是锂离子在正负极的脱嵌过程。放电时，锂离子从负极脱出穿过隔膜嵌入正极，电子通过外电路从负极到达正极；充电时，则反之。

通过下列反应式表示锂离子电池的工作原理：

总反应

正极反应

负极反应

式中，M表示Ni、Mn、CO等。

02锂离子电池产热模型

锂离子电池在正常工作中的产热主要包括电化学反应热 、焦耳热 、极化热 和副反应热 。

电化学反应热是锂电池在充放电过程中内部发生化学反应所吸收或释放的热量，充电过程中锂离子电池吸热， <0；反之锂电池放电时， >0。焦耳热是电流通过锂电池电阻所产生的热量，焦耳热在锂电池产热中具有较高的占比。极化热是电流通过极化内阻时由于焦耳效应所产生的热量，其产热方式与焦耳热相同。锂离子电池在正常工作时可忽略副反应热。

建模

01物理模型

仿真软件Comsol Multiphysics具有多物理场耦合的优点，本研究选用Comsol软件耦合集总电池接口、传热接口和流体流动接口来建立瞬态模型。

张晓光等建立相变冷却耦合空气冷却散热模型，通过模拟在不同单体间距下电池组的温度场来研究电池单体排布。结果表明，电池组内的温差随着单体电池间距的增大出现先减小后增大的变化趋势，单体电池的间距并非越大越好，在5mm时电池热管理性能最佳。本研究选取车载电池中的4节单体电池作为研究对象，间距均设为5mm，模拟在4C放电过程中锂电池温度等参数的变化。

实际应用中往往只能监测到锂电池表面温度，锂电池内部温度无法及时测量，这给锂电池的热管理增加了一定难度，且锂电池内部温度比表面温度高，因此监测锂电池的内部温度更具有研究价值。本研究在锂电池内部选取了6个点，研究其热物性参数的变化。图1中A、B、C、D、E、F的位置坐标为（0.000，0.000，0.0325）、（0.023，0.000，0.0325）、（0.009，0.000，0.0325）、（0.000，0.009，0.0325）、（0.0115，0.0115，0.0325）、（0.032，0.000，0.0325），位置坐标单位为m。该6点的位置坐标按照温度高低排列。后续分析均用1、2、3、4分别表示该模型中具体的锂电池，A、B、C、D、E、F分别表示其坐标位置。

图1 锂电池几何模型立体图

02传热模型

1、电化学-热耦合模型

考虑到锂电池内部实际结构的复杂性，进行如下假设：（1）锂电池的密度、比热容不随时间、温度变化而变化，忽略辐射散热；（2）锂电池内部受热均匀；（3）因锂电池组内各电池呈对称性分布，所以可将电池组四周视为绝热环境；（4）锂电池周围的空气假设为理想气体。

在多物理场耦合中选择电化学热接口即添加热源。电化学热接口基于上述假设，其锂电池导热能量守恒方程为

式中，ρ为锂电池平均密度，ρ=2523kg/ ； 为锂电池的平均比热容， =1145J/（kg·K）；T为微元体温度，K；t为时间，s； 、 、 为锂电池不同方向的导热系数，W/（m·K）；q为内热源，W/ ，r为单体锂电池的半径，m。

2、边界条件

本研究主要探究考虑温度效应时的对流换热系数对锂电池温度场的影响，因此选择传热学3类边界条件中的第三类边界条件，即选择Comsol Multiphysics热通量中的对流热通量，如下式：

式中， 为热通量，W/ ， 为锂电池外部空气的温度，℃。h采取流体横掠顺排管束平均表面换热系数计算方法。影响管束平均传热系数的影响因素主要有雷诺数Re、流体的普朗特数Pr。由于空气沿着进风方向的平均流速会发生变化，选定一个特征流速以计算Re。本研究分别选择0.1、0.3、0.5、0.7和0.9m/s作为平均流速。

100≤Re<1000时

1000≤Re<10000时

式中， ，表示动量扩散能力与热量扩散能力的一种量度；Nu为努塞尔数， ，表示壁面上流体的无量纲温度梯度；Re为雷诺数， ，表示惯性力与黏性力之比的一种量度；v为运动黏度，表示重力作用下流体内部流动阻力的量度， /s；α为热扩散系数， ，表示物体中某一点的温度的扰动传递到另一点的速率的量度， /s；h为对流换热系数，W/（ ·K），d为单体锂电池直径，m；λ为流体的导热系数，W/（m·K）；u为流体的流速（本文为空气流速），m/s， 为物性修正因子，本研究中 =1。具体计算结果如表1所示。

表1 对流换热系数计算结果

模拟结果及分析

01流速和环境温度对对流换热系数的影响

根据流体横掠管束的相关理论分别计算出空气流速为0.05、0.1、0.2、0.3m/s时各温度对应的对流换热系数，将计算结果与温度进行曲线拟合。影响对流换热的因素有很多，表征对流换热强弱的表面传热系数取决于流体速度v、黏度η、密度ρ、比热c，导热系数λ等多种因素的复杂函数，为方便数值模拟求解，需要对对流换热系数与温度进行曲线拟合，下列公式（10）~（13）分别是在不同空气流速下，对流换热系数与温度的拟合函数。图2~图5分别是在不同空气流速下对流换热系数与温度的拟合曲线图。

空气流速为0.05m/s时对流换热系数拟合多项式

空气流速为0.1m/s时对流换热系数拟合多项式

空气流速为0.2m/s时对流换热系数拟合多项式

空气流速为0.3m/s时对流换热系数拟合多项式

图2 空气流速0.05m/s时h与T的拟合曲线

图3 空气流速0.1m/s时h与T的拟合曲线

图4 空气流速0.2m/s时h与T的拟合曲线

图5 空气流速0.3m/s时h与T的拟合曲线

锂电池第一次充放电会在负极表面生成固体电解质界面（膜）（solid electrolyte interphase，SEI膜），SEI膜在80~100℃时会发生分解，失去保护负极的作用，本研究主要分析100℃以内对流换热系数随温度的变化。

图2、图3显示，空气流速为0.05、0.1m/s时的对流换热系数随着温度的升高而增大。温度在100℃内时，空气流速为0.05m/s时的对流换热系数随温度变化幅度较小，100℃时的对流换热系数与20℃的对流换热系数相比增长了3.60%。空气流速为0.1m/s时的对流换热系数随温度变化幅度较大，100℃时的对流换热系数与20℃的对流换热系数相比增长了11.11%。对流换热系数随着Nu、导热系数λ的升高呈增长趋势，Nu随着温度的升高而减小，空气的导热系数随着温度的升高而增大，对流换热系数与温度呈非线性相关（图4、图5）。空气流速为0.2m/s时的对流换热系数随着温度的升高呈现增—减—增的趋势。空气流速为0.2m/s时，100℃时的对流换热系数与20℃时的对流换热系数相比变化了0.79%。随着流速的增大，Nu数对对流换热系数的影响大于空气的导热系数，空气流速为0.3m/s时对流换热系数呈现出减—增—减的趋势。在该空气流速下，100℃时的对流换热系数与20℃时的对流换热系数相比变化了0.80%。

02对流换热系数对锂电池温度场的影响

将上述计算出的各空气流速对应下的对流换热系数温度函数应用在模型仿真中，空气温度、环境温度均为50℃，分别模拟计算出空气流速为0.05、0.1、0.2、0.3m/s，考虑温度效应时的对流换热系数对锂电池温度场的影响。图6~图9分别是空气流速为0.05、0.1、0.2、0.3m/s时，锂电池放电过程中对流换热系数考虑温度效应时对温度场的影响模拟结果。

图6 空气流速0.05m/s、考虑温度效应时锂电池温差变化

图7 空气流速0.1m/s、考虑温度效应时锂电池温差变化

图8 空气流速0.2m/s、考虑温度效应时锂电池温差变化

图9 空气流速0.3m/s、考虑温度效应时锂电池温差变化

空气流速为0.1m/s、锂电池放电至729s时，考虑温度效应时的对流换热系数的温差比对流换热系数为定值时降低了21.71%。图6~图9显示空气流速为0.1m/s时考虑温度效应的对流换热系数对锂电池温度场的影响较大，锂电池的温差随着对流换热系数的增大而减小，增大了模拟研究的准确性。空气流速为0.05、0.2、0.3m/s时考虑温度效应的对流换热系数使锂电池温差变化值均小于1.00%，对锂电池温度场影响较小。空气流速为0.05m/s时，温度对对流换热系数影响程度较小，在该流速下考虑温度效应的对流换热系数对锂电池温度场的影响可以忽略不计。随着空气流速的增大，强制对流散热效果远高于自然对流散热，因此，空气流速为0.2、0.3m/s时，考虑温度效应的对流换热系数对锂电池温度场影响也可以忽略。

由以上分析可知，不同流速下对流换热系数受温度的影响程度不同，为了减少误差，建议在数值模拟锂电池散热时考虑温度对对流换热系数的影响，若不同温度对流换热系数影响较大，则需要考虑温度对对流换热系数的影响，反之则可以忽略。

03锂电池内部不同点的温度分布

通过研究得到，空气流速为0.1m/s时，对流换热系数受温度影响最大。在此研究基础上，研究流速为0.1m/s时，对流换热系数考虑温度效应时锂电池内部温度的变化情况。根据模拟结果，得到锂电池1内部6个点放电过程中的温度变化（图10），同时还得到锂电池1、2放电至550s时在zx截面的温度云图（图11）。

图10中显示，A点温度一直最高，A点是1号电池的轴心位置，B点是2号电池的轴心位置。锂电池内外温差最大为2.4℃。根据相关文献得知，锂电池热量积聚发生热失控时会呈现出迅速且不易被察觉的特点，内部最高温度要比表面温度高200℃，因而锂电池内部温度比锂电池表面温度更具有研究价值，掌握锂电池内部温度能更好地对锂电池建立预警体系，便于在发生热失控之前及时将锂电池温度控制在安全范围内或者对锂电池采取紧急降温措施。

图10 电池内部点的温度变化

图11 550s时电池1、2的zx截面温度云图

结 论

在COMSOL软件中建立了电化学-热耦合模型，并在该模型中添加了空冷散热模型，模拟对流换热系数考虑温度时对锂电池热扩散影响及其程度，得到以下主要结论。

1）空气流速和温度均对对流换热系数有影响，不同流速下对流换热系数与温度并不完全呈线性变化，分别模拟了空气流速为0.05、0.1、0.2、0.3m/s的锂电池放电过程，其中空气流速为0.1m/s时对流换热系数随温度变化幅度最大，100℃时的对流换热系数相比与20℃的对流换热系数增长了11.11%。

2）空气流速为0.1m/s、在729s时对流换热系数考虑温度影响比设为常数时锂电池最大温差降低了21.71%；但空气流速为0.05、0.2、0.3m/s时对流换热系数的温度函数使锂电池的温差变化量均小于1.00%，在此流速下，对流换热系数受温度影响可忽略。

3）空气流速为0.1m/s、环境温度和空冷温度均为50℃、放电倍率为4C时锂电池内部温度比锂电池表面最大温度高2.4℃。

以上数值模拟研究发现，对流换热系数与温度并不完全呈线性变化，需要先计算温度对对流换热系数的影响，若计算得出温度对对流换热系数的影响小于3%时，则可以忽略因温度对对流换热系数影响所带来的模拟误差，否则需要考虑温度对对流换热系数所带来的模拟误差。

本文作者：韩润利、唐明云、王冬、张少杰

作者简介：韩润利，安徽理工大学深部煤矿采动响应与灾害防控国家重点实验室，安徽理工大学安全科学与工程学院，硕士研究生，研究方向为锂电池热管理系统；唐明云（通信作者），安徽理工大学深部煤矿采动响应与灾害防控国家重点实验室，安徽理工大学安全科学与工程学院，教授，研究方向为火灾防治。

原文发表于《科技导报》2023年第5期，欢迎订阅查看。