基于数字风洞的整车热管理性能开发.docxVIP

基于数字风洞的整车热管理性能开发

摘要：针对某商用车更换大功率发动机以满足国六排放和四阶段油耗法规要求，需要优化冷却系统，提升散热能力.采用格子玻尔兹曼方法（latticeBoltzmannmethod，LBM）建立热环境数字风洞，高精度预测整车热平衡与热保护性能，仿真与试验对标，冷却液温度预测值与试验测试值之间的误差小于1℃.基于此以散热器水温和中冷器出风温度作为优化目标，分析冷凝器、中冷器、散热器和风扇侵入量之间的交互效应，并据此进行优化设计，提出一套适用于发动机冷却系统虚拟标定的仿真流程.研究结果表明，风扇侵入量和中冷器高度对散热器散热性能影响较大.经过优化迭代后，散热器和中冷器的进气量分别提升了5.01%和7.87%，冷却模块的表面温度分布更加均匀，显著改善了发动机舱散热效率.

关键词：热平衡；数值模拟；玻尔兹曼方程；优化

中图分类号：U461文献标志码：A

随着运输行业的蓬勃发展，加之国家出台的油耗法规与排放法规，商用车开发时格外关注机舱热管理性能，这一性能的提升对满足日益严格的汽车油耗性能需求和排放标准至关重要.冷却系统的设计与布置直接影响整车的运行效率和可靠性，如果冷却系统设计不匹配或布置不合理会导致热风回流或逆流，进而升高发动机进出水温度，这不仅影响发动机充气效率和燃烧稳定性［1-3］，而且可能损害热敏元件，甚至引发安全事故.因此，提高整车散热性能成为当今汽车研发的难题［4-6］.

本文使用PowerFLOW软件进行仿真分析，建立数字风洞精确预测发动机水温和热敏元件表面温度.采用一维与三维耦合仿真技术，分析影响整车冷却系统散热性能的关键因素.通过多目标优化，探索冷却系统设计的优化空间，提升商用车散热性能和安全性，为商用车冷却系统设计提供参考［7-10］.

1数值方法

1.1用格子玻尔兹曼方法模拟流场

为了模拟车辆在热环境风洞中的非定常流动行为，本文采用PowerFLOW软件提供的格子玻尔兹曼方法（LBM）.与传统的计算流体力学（computationalfluiddynamics，CFD）方法相比，LBM方法基于数字物理技术，通过统计的方法描述真实流体运动，使用离散化的Boltzmann方程来完成流体的建模.这种方法在基本的动力学层面上对流体进行模拟，从而提供了一种更为精确和高效的流体分析手段.采用Bhatnagar-Gross-Krook（BGK）［11］形式的格子玻尔兹曼方程，求解流体分子之间的碰撞松弛达到平衡状态的过程，如式（1）所示.

式中：f≡f（x，ξ，t）描述了在时间t、位置x和速度ξ处的单粒子数密度f；τ为碰撞引起的松弛时间.

图1为Boltzmann方程的有限差分格式.使用D3Q19模型在三维立方晶格上离散化Boltzmann方程和BGK算子，由此产生的有限差分方程能够求解分布函数.

麦克斯韦-玻尔兹曼分布函数描述如式（2）所示.

式中：R为气体常数；D为空间的维数；ρ为质量密度；θ为归一化的标准温度，θ≡KBT/m；KB、T、m、u分别为玻尔兹曼常数、温度、粒子质量和流体流速.LBM跟踪分子在空间和时间中的运动来模拟气体和液体的流动.它遵循质量、动量和能量的守恒关系，宏观流体力学量是颗粒密度分布矩阵的直接结果，由式（3）给出.

式中：ε为内能密度.LBM中使用了3个微观过程来模拟流体行为，包括颗粒间相互作用、颗粒表面相互作用和平流颗粒运动.在基于数学模型使用的格子玻尔兹曼模型中，粒子存在于空间中的离散位置，如图2所示，粒子位于由体单元和曲面组成的格子上.其中曲面是物体的表面，体单元的格子与物体表面进行切割，形成用于计算的流体网格.在模拟过程中，粒子在离散的时间间隔内以离散的速度从一个体单元移动到另一个体单元时发生相互作用，当数百万个粒子在各体单元的每个方向上移动时，遵循粒子密度分布函数.相同体单元中的粒子碰撞并改变速度和方向，同时保持体单元中总质量、动量和能量的守恒.表面碰撞过程加强了质量守恒，同时与所需的边界条件交换动量和能量.动量的变化量对应于压力和摩擦力.在上述体单元到表面的流过程之后，聚集的粒子被反射并返回到流体体单元中［12］.

LBM的主要特点：它是一个高效的非定常流求解器，数值耗散低.它能够处理非常复杂的几何结构.它使用超大涡模拟（verylargeeddysimulation，VLES）湍流模型，包括子网格建模和壁面函数建模，直接求解各向异性湍流涡或非常大的漩涡.

浮力驱动的热对流通过CFD进行求解.热传导系数为：

CHT=q/（Twall-Tref）（4）

式中：q为热通量；Twall为表面温度；Tref为参考温度（近壁温度或特征温度）.

湍流建模模拟中