广告
一种自动高效的发动机舱热管理建模方法
2020-08-29 23:41:02· 来源：汽车CFD技术之家
摘 要：汽车发动机舱热管理的CFD 分析需要对整车包含发动机舱内的众多零部件进行复杂的网格处理，而现有的工作流程主要依赖手工方式进行，工作量巨大，因此迫切
摘 要：汽车发动机舱热管理的CFD 分析需要对整车包含发动机舱内的众多零部件进行复杂的网格处理，而现有的工作流程主要依赖手工方式进行，工作量巨大，因此迫切需要一种自动化流程来实现。本文以我司某车型为例，利用FluentMeshing的自动化网格流程来完成表面网格处理和体网格生成，将原本需要1~2个月的网格处理工作缩短到1周之内完成，且人工干预工作大大减少，大幅提升了仿真工作的效率。该自动化网格流程内嵌最佳网格实践，同时也可根据我司不同车型几何质量及管理规范进行二次封装以适应多种车型（或改款）的自动高效建模需求。通过这套建模流程可以获得高质量体网格，模型计算结果（热态工况）与实验结果吻合较好。
引言
计算机仿真技术可以有效的缩短研发周期降低研发成本，在激烈的汽车市场竞争中发挥着重要作用。随着汽车工业的高速发展, 人们对汽车安全性、舒适性、节能性和环保性的要求也越来越高。发动机舱热管理仿真是汽车行业数值仿真一个重要方向。发动机舱里布置了包括前端冷却模块、发动机总成、进排气系统，电子部件等众多零部件，其中发动机和排气系统是高温热源，热源附近的温度敏感部件（如线束、橡胶件、电子零件等）受到高温影响可能会导致零部件加速老化甚至性能损害，因此通过对发动机舱内部的流动和温度分布的仿真，对机舱内关键部件进行热害分析，并为空间布局优化提供技术支持。CFD仿真可以缩短试验周期，降低试验成本。
发动机舱热管理分析需要对整车包含发动机舱内的众多零件进行网格处理，现有的CFD工作流程主要依赖手工包面的方式处理，这种方式需要包面之前对初始集合模型进行预处理，手动检测并修补发动机等零部件，手动预处理的工作量巨大，并且受人工影响因素大，表面网格质量很难保证，迫切需要一种自动化流程来实现。FluentMeshing的脚本驱动处理网格正是这样的自动化流程。它可以自动完成表面网格处理和体网格生成，需要人工干预的工作减少，大幅提升了仿真工作的效率。优化仿真中批量处理大规模模型的需求，更能体现出脚本处理的优势。该自动化网格流程适用于不同车型，可以积累最佳网格实践数据库，适应多种车型（或改款）的自动高效建模需求。
利用这套方法,Volvo S80车型仅用3天时间就可获得整车网格，并且这套一体化网格可以进行包括外气动，发动机舱热管理，共轭传热等计算[1]。这种全新高效的建模流程极大提升了Volvo热管理分析的效率，也更加印证了自动化网格流程的发展趋势。
下面对我司某车型利用自动网格流程进行发动机舱热管理仿真，并对关注部件进行热害预测的过程进行介绍。
1、自动化网格流程
自动化网格流程就是对初始几何模型运行一套已经准备好的脚本文件，就可以自动完成表面漏洞封闭、包面、面网格重构、体网格划分及网格质量提升的全套工作。
发动机舱热管理仿真的零部件众多，要仿真空气从格栅进入在整个发动舱内的流动和温度分布，需要通过包面的方式提取出流体域的边界。要对空气域进行包面，就需要把空气域处理成一个单连通的区域，这样就需要对发动机，乘员舱等主要区域的表面进行漏洞检测，封闭表面漏洞。传统的网格处理方式中漏洞封闭的过程全部需要工程师手动处理，工作量巨大，对于一个整车模型，网格处理需要一到两个月的时间。网格处理的工作时长和网格质量受工程师使用软件的熟练程度影响，差别较大。
1.1、自动化网格流程的优势
FluentMeshing的自动化网格处理可以依靠脚本驱动自动的实现发动机舱热管理这样复杂模型的处理，自动完成漏洞封闭、包面、重构、体网格划分及网格质量提升的工作，使用者只需要对初始网格模型进行模型分区命名等准备工作，调整合理的尺寸函数后，运行脚本文件即可。这样网格处理的时间可以缩短到1周之内完成，并且其中手动准备时间更少，其余都是自动运行的CPU时间，大幅提升了工作效率。
自动化网格流程的自动封闭漏洞功能非常强大。以发动机为例，通过参数设定，我们可以控制检测漏洞的尺寸，运行网格流程后自动封闭所有小于指定尺寸的漏洞，得到封闭的发动机外表面，免去了繁杂重复的漏洞检测和堵洞的手工操作。
自动化网格流程中可以根据车外造型自动生成加密区域，包括贴体加密区域，前端加密区域和尾流加密区域，并且可以通过调整参数控制加密区域的形状，尺寸和密度[2]。
自动化网格流程中的一个特色功能是自动缝隙封闭的功能，可以针对指定的对象按照定义的标准寻找并封闭缝隙，对于外气动或者发动机舱热管理分析中如外造型面上不同部件间的缝隙封补非常有用[2]。
此外，自动化网格工具还包含同时生成固体网格的功能，用于流固耦合传热计算。自动匹配交界面等工具，省去了发动机舱热管理仿真中众多交界面的选择和匹配工作。
对于发动机舱的优化仿真，往往需要进行批量的模型处理，这些模型如果在命名规则一致的前提下，就可以不再重复模型的尺寸函数等设定工作，直接运行脚本，实现批量处理。
使用自动化网格流程可以不断积累针对不同车型的最佳网格实践，同时也可根据我司不同车型几何质量及管理规范进行二次封装以适应多种车型以及改款车型的自动高效建模需求。
1.2、自动化网格流程
自动化网格流程使用的FluentMeshing脚本语言易读易用，使用者不需要精通编程语言就可以使用。使用这套网格流程，从包面到体网格生成再到网格质量提升都自动完成，网格处理时间更短，人工干预更少，网格质量更有保证。
脚本化网格的流程如图3所示：
2、数学模型
计算流体动力学（Computational Fluid Dynamics，简称CFD）是流体力学领域的重要分支学科。采用CFD方法研究问题时，常用方法是把空间区域离散成若干三维单元，以形成一个个立体网格( 或者格点)，然后在网格上应用相应的算法来解连续方程、运动方程（对于无粘流体求解欧拉方程，对于粘性流体求解纳维-斯托克斯方程(Navier-Stokes equations)）、能量方程及组分输运方程等，对于瞬态过程现象，还要添加物理变量对时间的离散项，在离散的时间步长进行计算。目前在工程领域CFD方法已经得到广泛的应用。
2.1、湍流模型介绍
本文发动机舱热管理仿真用到的湍流模型是工程仿真分析中使用最广泛的是k-ε系列湍流模型。k-ε模型是一种两方程湍流模型，它通过求解两个附加输运方程（湍动能k和湍动能扩散率ε），再以此为基础得到湍流粘度，从而对雷诺时均RANS方程中的附加雷诺应力项进行模化，使整个方程组封闭。标准k-e模型自它被提出以来就在诸多实际工程流动计算中扮演重要的角色，它以其良好的计算经济性、宽广的实用性以及合理的精度被许多工程技术人员的广泛使用。
湍动能k和湍动能耗散率ε从以下运输方程得到：
其中，Gk是由于平均速度梯度引起的湍动能k的产生项，Gb是由于浮力引起的湍动能k的产生项，YM代表课呀湍流中的脉动扩张的贡献，C1ε、C2ε和C3ε为经验常数，σk和σε分别是与湍动能k和耗散率ε对应的Prandtl数，Sk和Sε是用户定义的源项[3]。
2.2、传热和辐射模型介绍
传热在各工业领域都存在广泛的应用。称从空间中占据一个区域的物质的热能到空间中占据不同区域的物质称为热传递。传热可以通过三种主要方法进行：导热，对流和辐射。本文发动机舱热管理仿真要计算传热和辐射。
当流体运动时，自身携带热量，这称为对流。传热是跟流体流动耦合在一起的。当能量和流体流动方程都被激活就意味着计算了对流。当能量方程激活时，流体内的导热也进行了计算。
能量运输方程
其中，keff是有效导热率（k+kt，其中kt是根据所使用的湍流模型定义的湍流热导率），Jj是组分j的耗散通量。右侧的前三项分别表示由于传导，组分扩散和粘性耗散引起的能量传递[4]。
辐射传热是一种通过电磁波传递能量的能量传输模式。热辐射覆0.1~100um的电磁光谱部分。对于半透明物体 ( 如玻璃、燃烧废气), 辐射可以从体内发射，因此是一种体现象。对于不透明物体，辐射是一种表面现象，因为体内发射的电磁波都在本体内部被吸收。
在一个模拟中，如果辐射换热热流跟对流和导热传热处于相同数量级时需考虑辐射。
通常在高温条件下符合这种情况。但某些较低温度的情况下也符合，具体取决于应用对象。可先根据下式估计系统中导热或对流传热的量级，再将qrad与qconv进行比较。
Fluent中的辐射模型包括P1模型，DTRM(Discrete Transfer Method) 模型，S2S(Surface to surfacemodel) 模型，DO(Discrete ordinates model) 模型等。在选择辐射模型前，先要确定光学厚度：
光学厚度小意味着流体介质对辐射接近透明，辐射仅在计算域的边界上互相传递。光学厚度大意味着流体介质会吸收辐射，并再次发射辐射。
DO模型为对计算资源要求最高的辐射模型，同时适用性最广，精度最高。当光学厚度=0时，S2S具有与DO模型相当的精度，但所需计算资源更少。P1模型计算资源要求中等，精度尚可。DTRM模型比DO模型计算代价小，但不支持并行，因此很少使用。对于发动机舱热管理仿真，可以选择S2S模型。
3、仿真建模
3.1、几何模型
发动机舱热管理需要处理全车数模（乘员舱内部件除外），按照车身，动力总成，冷却模块和机舱内部件等对几何边界进行预处理归类命名等工作，为运行自动化网格流程做准备。初始几何模型外造型面,前端冷却模块和动力总成结构如图4-6所示：
3.2、网格模型
利用自动化网格流程包面并重构之后的面网格数量约1500万，其中换热器模块和风扇域边界等部件的最小尺寸为0.5mm。面网格如图7-9所示。
主流体域空间体网格类型为四面体网格和边界层，单独划分换热器网格和风扇域网格以及隔热罩固体网格，组合脚本自动化生成的流体域网格，体网格总数约6000万。体网格截面如图10所示 。
3.3、分析设定
计算工况为低速爬坡，车速50km/h，环境温度38℃。空滤吸气量为457.78kg/h， 风扇转速为2650rpm。换热器芯体定义多孔介质模型，发动机和排气系统设定温度热源，并考虑辐射影响。散热器采用Fluent的换热器模型定义。
3.3.1、换热器模型
换热器模型用于预测发动机舱环境中换热器组件的冷却效果，如预测辅助流体入口温度或者整体发动机放热量。ANSYSFluent中有两种可用的换热器模型：Macro模型和Dual-Cell模型。
Macro模型只需要创建一套网格单元，真实的单元用于主流体（冷却流体）一侧，macros用于辅助流体（热流体）一侧，在发动机舱热管理仿真中是普遍使用的模型。Dual-Cell模型需要创建重叠单元，一个代表主流体一个代表辅助流体，可以预测辅助流体的流动和特性[5]。
3.3.2、大规模并行计算
在并行计算规模不断增大的背景下，ANSYS Fluent 软件在HPC并行能力上持续发展，求解器优化，计算更快，并行效率更高。在硬件不变的条件下，软件并行效率的提升帮助我们在处理大规模计算时摆脱硬件的限制。Fluent目前是超级计算多核并行的世界纪录保持者[6]。ANSYS具有GPU并行的功能，对于计算量较大的模型还可以选择性得使用分解的CPU区域和GPU区域。
4、仿真结果分析
4.1、流场结果分析
流场分析可以得到发动机舱内的三维详细流场特点并分析其潜在影响。我可以输出的结果包含前端冷却模块(冷凝器、散热器、中冷器和风扇)的计算风量、机舱内关键部件周边的流场状况。
统计格栅的进风量，计算通过冷却模块的风量占格栅进风量的比重，可以辅助分析出冷却模块的冷却效果。例如可以通过增加导风板，调整冷却模块局部布局等方式提高进风量，从而提升冷却能力。
以下是Y向截面的速度分布：
发动机后方的排气系统属于高温部件，排气附近的风速较小会导致局部温度过高，因此可以调整发动机装饰罩等部件的位置及形状，增大发动机后方的空气速度，有助于降低热害风险。通过分析发动机舱内重点关注部件附近的流场分布，可以辅助发动机舱内的优化布局，再结合热态计算的结果，可更直观更准确的分析温度结果，预测潜在热害。
4.2、热态结果分析
以下分别是动力总成系统的温度分布和Y向截面的速度分布：
下图是散热器迎风面和背风面的温度分布：
以下是散热器通道分布，以及模拟出的水侧发热量分布：
5、结论及展望
发动机舱热管理仿真可以分析机舱内的流场分布, 如果结合一维软件比如GT－SUITE还可以进行一维三维耦合，计算换热器的匹配。温度场分析可以预测关注部件的潜在热害，为机舱的优化提供仿真依据。
自动化网格流程大幅缩短了仿真前处理工作的周期，对于提升仿真工作效率，缩短研发周期起到了重要的推进作用。本文的模型中包含了仅隔热罩的固体模型，我们后期会尝试增加更多的固体模型，开展流固耦合传热计算。后续我们会继续拓展自动化工作的流程，尝试将模型设定也通过脚本驱动来自动完成，进一步完善自动化建模流程。
基于脚本化的网格流程，ANSYS已经开发了ACT模块，它是基于最佳实践的快速网格生成模板，导入几何文件再输入部分管理数据，可对复杂的几何模型进行包面和体网格生成，支持发动机舱、气动外流场、舱室冷却的网格生成。ACT模块可以最小化人工干预设置，伴随生成过程自动生成Excel表格以进行输入数据的重用。我们在积累自动网格流程经验的基础上，进一步尝试ACT模块的应用。
作者：訾昌陆1，陈桂杰2，唐连伟2，钟修林2
- 浙江吉利汽车研究院有限公司
- 艾迪捷信息科技上海有限公司
广告 
广告
编辑推荐
最新资讯
-
荷兰Zepp氢燃料电池卡车-Europa
2024-12-22 10:13
-
NCACFE -车队油耗经济性报告（2024版）
2024-12-22 10:11
-
R54法规对商用车轮胎的要求（上）
2024-12-22 10:10
-
蔚来ET9数字架构解析
2024-12-22 09:53
-
4G/5G网络新时代的高效紧急呼叫系统NG-eCal
2024-12-20 22:33
广告
广告
