1. 引言
当环境温度低于−40˚C,不通过其他手段或措施,汽车无法靠本身的起动装置自行起动 [1] 。为了解决这类冷启动问题,有的在发动机外部布置一个锅炉加热水系统。该系统在发动机启动前将发动机内部冷却液从柴油机某合适部位引出,经过锅炉系统的循环水泵进入锅炉,经过锅炉加热后的冷却液返回柴油机内部循环流动加热发动机。该系统设计有几个关键点需要注意:一是加热速度快,最简单的办法是加大锅炉功率及循环水流量,但是这带来耗能多、布置空间紧张等的问题,因此要求锅炉热水尽可能多的流经发动机本体,尽量不要走发动机内部小循环通道流动;二是各缸加热尽量均匀;三是发动机正常启动后锅炉内部循环不要出现反向流动导致锅炉损坏;四是单向阀目前在温度−50˚C~150˚C大范围工作可靠性及耐久性差,尽量不用。
由于在低温状态下节温器关闭了大循环,因此计算不考虑大循环管路及水箱等部件,只是针对柴油机本体水套进行,用加热水经过发动机各缸的流量评价加热均匀性,用经过小循环通道被短路的流量占比来评价加热的快速性。
2. 锅炉加热水系统优化计算
2.1. 计算模型及方案
本次计算水套网格模型如图1所示,使用Hypermesh对几何模型进行处理,生成面网格,使用AVL Fire自动网格生成器FAME HYBRID生成体网格,最大网格尺寸4 mm,最小尺寸1 mm,总网格数约425万个。
节温器座小循环通道内在某些方案内安装有单向阀(如图2所示),针对有单向阀的情况计算时直接认为通道完全不通。
![](//html.hanspub.org/file/7-2810074x9_hanspub.png)
Figure 1. Cooling water jacket mesh model
图1. 水套网格模型图
计算方式为标准的求解N-S方程组迭代求解方式:给定初场,然后通过联立求解每个网络单元上的N-S方程组,不断迭代,直到计算达到收敛。计算中对于压力项和速度项之间的耦合关系采用SIMPLE算法,使用的湍流模型为k-ε两方程湍流模型,近壁面采用了标准壁面函数法。
计算方案具体见表1。方案5是在方案2基础上将节温器接口1管内径由原来的10 mm加大到14 mm,机体横向水道接口内径由原来的12 mm加大到18 mm。
2.2. 边界条件
采用稳态计算,进水口流量根据锅炉水泵性能按20 L/min计算,进口水温按试验温度−40˚C设置;出水口设置为压力:100,000 Pa;其余为壁面边界,考虑重力,不考虑传热。工质为water/glycol,体积分数各50%。
2.3. 计算结果与分析
模拟结果主要从压力、流速、流量等方面进行分析说明。
2.3.1. 压力计算结果
各方案进出口压差对比柱状图如图3所示。除了方案5加大进出口直径阻力非常小外,其余各方案总压差在16~21 kPa之间,相差不大。
![](//html.hanspub.org/file/7-2810074x10_hanspub.png)
Figure 2. Non-return valve assembly
图2. 节温器座小循环管道内安装单向阀的装配图
表1. 计算方案列表
2.3.2. 流速计算结果
各方案缸垫上水孔截面速度矢量图如图4所示,从图中可以看出方案3和方案6缸垫孔处的流速相对较高和较均匀,方案1各缸均匀性还可以,但是由于小循环短路流量占比大,通过机体缸盖的水流量小。
图5为方案6节温器座附近截面速度矢量图,从图中可以看出,在方案6节温器小循环口周围存在较大速度的垂直流动,产生的水墙阻碍水从小循环通道流动,这就是小循环流量低的原因。方案7稍微改动节温器部位的进水位置小循环流量就大幅度增加,因此一定要注意选择节温器接口1,不要选择接口2。
2.3.3. 流量计算结果
各方案对应各缸缸盖口流量占比如图6所示,从图中可以看出:方案6各缸流量相对较大,经过机体缸盖的流量占比达到92.9%,只有7.1%的流量经过小循环被短路;且方案6各缸比较均匀,第2缸流量最小,占比12.8%,第5缸流量最大,占比18.3%。方案6迹线图及对应支路流量占比结果见图7。其余方案结果在此不展示。
2.3.4. 各方案关键计算结果对比
表2为各方案关键结果评价,综合评价认为方案6最优。
3. 锅炉加热水系统试验验证
图8为各方案对应各缸缸盖口18分钟温升测试数据对比(方案5由于需要改变零件尺寸未验证,由于方案6结果已经满足整车厂要求,且台架资源紧张,方案7未再验证),从图中可以看出:方案6各缸出水口最低温升达到57.1˚C,比其他方案中较优的方案3 (40˚C)高17.1˚C;平均温升达到65.6˚C,比其他方案中较优的方案4 (49.5˚C)高16.1˚C。
4. 结语
通过AVL FIRE对某柴油机水套加热系统进行了多方案对比分析,优选出了各缸加热均匀、加热速
![](//html.hanspub.org/file/7-2810074x11_hanspub.png)
Figure 3. Pressure loss results of all cases
图3. 各方案进出口压差对比柱状图
![](//html.hanspub.org/file/7-2810074x12_hanspub.png)
Figure 4. Velocity distribution at cylinder gasket hole of all cases
图4. 各方案缸垫上水孔截面速度矢量图
![](//html.hanspub.org/file/7-2810074x13_hanspub.png)
Figure 5. Velocity distribution near thermostat of Case 6
图5. 方案6节温器座附近截面速度矢量图
![](//html.hanspub.org/file/7-2810074x14_hanspub.png)
Figure 6. Mass flow percentage at each cylinder inlet of all cases
图6. 各方案对应各缸缸盖口流量占比
![](//html.hanspub.org/file/7-2810074x15_hanspub.png)
Figure 7. Mass flow percentage at each cylinder inlet of Case 6
图7. 方案6迹线图及对应支路流量占比
![](//html.hanspub.org/file/7-2810074x16_hanspub.png)
Figure 8. Tested temperature rise at each cylinder inlet of all cases
图8. 各方案对应各缸缸盖口18分钟温升测试数据对比
![](Images/Table_Tmp.jpg)
Table 2. Simulation results assessment of all cases
表2. 计算方案关键结果评价
度高的设计方案,优化方案得到了试验验证。