白曙

（上海柴油機股份有限公司，上海200438）

沸騰對柴油機冷卻系統(tǒng)數值模擬的影響

白曙

（上海柴油機股份有限公司，上海200438）

計算對比了不考慮沸騰和考慮沸騰2種冷卻系統(tǒng)數值模擬計算，得出結論：考慮沸騰傳熱對內燃機冷卻水腔內流動與壓力的分布影響不明顯，而對冷卻水腔內傳熱過程的影響是很大的。若只考慮純對流傳熱，計算結果可能與實際情況存在很大的差異。因此，在對強化內燃機進行流動與傳熱問題的研究時必須考慮沸騰傳熱的因素，以獲得更為真實、準確的結果。

內燃機 沸騰傳熱 冷卻水腔 流-固耦合 有限元分析

1 前言

柴油機冷卻水腔內的傳熱是一個極為復雜的過程。在水冷式柴油機氣缸蓋中，冷卻水腔內的傳熱有對流傳熱和沸騰傳熱2種方式[1]，在對強化內燃機進行流動與傳熱問題的研究時，國內外對沸騰傳熱的研究都比較少。本文希望通過對比計算，檢驗在對強化內燃機進行傳熱問題研究時，是否有必要考慮沸騰問題。

2 三維組合模型的建立

計算模型建立的準確性、有效性直接關系到計算的精確性，是有限元模型的基礎。本文的計算對象主要是由氣缸蓋、氣缸墊、機體、氣缸套、機體冷卻器等零部件裝配而成。通過三維造型，再現柴油機的實際結構。由于氣缸蓋、機體上一些尺寸較小的螺釘孔、銷孔和測試孔等對有限元計算結果的影響很小，而且在對這些小孔劃分網格時容易產生不良單元，因此，這些結構在建模過程中省略。如圖1所示，圖中(a)為機體；(b)為氣缸套；(c)為氣缸墊；(d)為氣缸蓋；(e)為機體冷卻器。將這些零部件裝配起來，得到組合模型，如圖2所示。

將裝配好的模型以STL格式導入到SC/Tetra軟件的前處理模塊中。由于SC/Tetra軟件具有自動識別模型區(qū)域的功能，即只要將該柴油機裝配模型中冷卻水腔的進出口封閉起來，使其形成一個封閉的空間，SC/Tetra軟件就自動建立了該柴油機冷卻水腔的模型，如圖3所示。

3 有限元模型的建立

柴油機熱分析的有限元模型建立過程包括計算網格的劃分，邊界條件的確定，材料屬性的設定。

來稿日期：2012-12-12

圖1 某柴油機零部件三維模型

圖2 某柴油機三維組合模型

圖3 水腔幾何模型

3.1 計算網格的劃分

在劃分流固耦合交界面網格時，采整體劃分的方法同時生成固體域和流體域的網格，這樣就保證了2個區(qū)域的網格在流-固交界面上完全一致。

圖4是固體域網格模型。對于機體、缸蓋的流體域網格構成則比較復雜：在遠離壁面的區(qū)域采用的是非結構化四面體網格，在近壁處插入由五面體的三棱柱和金字塔網格構成的附面層網格。圖5是流體區(qū)域網格模型。由于簡化的缸墊結構外形較為規(guī)則且較?。▔壕o后的厚度僅為1.2mm），為了減小計算規(guī)模，缸墊本體及其水孔均采用了較為規(guī)則的六面體網格。

圖4 某柴油機零部件固體域網格模型

圖5 某柴油機零部件流體域網格模型

將這些網格模型裝配在一起形成了某型號柴油機機體-缸蓋-缸套-缸墊-冷卻水腔的流-固耦合計算模型。整個模型網格總數為6 213 765個，網格節(jié)點為1 710 816個。

3.2 邊界條件

3.2.1 流動邊界條件

冷卻液入口設定為入口質量流量為2.9 kg/s和入口溫度為78℃；冷卻液出口壓力設為相對靜壓邊界，取值為0 Pa。流體壁面設定為無滑移邊界，固體外壁面設定為自由滑移壁面。

3.2.2 缸內燃氣側的傳熱邊界條件

由Boost軟件仿真得到，某型號柴油機一個循環(huán)周期內缸內等效燃氣溫度和平均傳熱系數αgm分別為1 043 K和927 W/(m2·K)。

（1）缸蓋火力面燃氣側的傳熱邊界條件

根據文獻推薦，確定了柴油機缸蓋火力面的分區(qū)方式及各子區(qū)域的局部傳熱系數，如圖6和表1所示。

圖6 燃燒室表面分區(qū)示意圖

表1 各分區(qū)的傳熱邊界條件

對缸套內壁面進行了分區(qū)處理，并將各區(qū)域的平均傳熱系數和等效燃氣溫度施加到對應的缸套內表面上，具體數值如表2所示。

3.2.3 進、排氣道的傳熱邊界條件

根據文獻推薦，取進氣道內空氣的換熱系數250 W/(m2·K)，環(huán)境空氣溫度為325 K；排氣道表面?zhèn)鳠嵯禂?50 W/(m2·K)，廢氣溫度為795 K，將這些傳熱邊界條件分別施加于進、排氣道內表面單元上，如圖7所示。

3.2.4 其它傳熱邊界條件

機體、缸蓋的外壁面取平均傳熱系數為20 W/ (m2·K)，周圍的環(huán)境溫度為323 K；螺栓表面和螺栓孔內表面近似按絕熱邊界處理；氣缸墊與機體、缸蓋有接觸的表面均按絕熱邊界處理。

3.3 材料屬性

某型號柴油機所使用的冷卻介質為清潔的天然軟水，其熱物性參數隨溫度變化關系見表3。機體、氣缸蓋等材料參數見表4。

3.4 BDL沸騰模型

計算流動沸騰傳熱時，其中沸騰部分的傳熱系數計算公式可以借用靜止液體沸騰(容積沸騰)的計算公式。本文采用的BDL模型，是在Chen模型[3]的基礎上提出的一種單向流沸騰模型。Chen模型中S的取值，依賴于全局雷諾數的大小，但是全局雷諾數描述的是整體流，因此Chen模型忽略局部流動參數的影響，損失了部分計算精度。BDL模型的優(yōu)點就在于其克服了Chen模型的缺點，針對系數S進行了修改，考慮了局部流速以及長度尺寸等因素。

表2 氣缸套內壁面換熱邊界條件

圖7 進、排氣道表面網格模型

3.5 求解與收斂判斷

為了判斷計算是否收斂，一般在計算開始之前都要設置殘差監(jiān)控器和某個特定位置處的變量監(jiān)控器。本文中設置計算殘差為10-4，選定的監(jiān)測面為某個上水孔截面。計算后的殘差變化曲線如圖8所示，某個上水孔的質量流量變化如圖9所示。

在求解時只有當計算殘差滿足收斂要求以及所監(jiān)視的截面上的變量趨于穩(wěn)定時方可結束本次計算。在相同的收斂準則、計算步長（2000步）和計算CPU數目（4CPU）下，通過計算發(fā)現，不考慮沸騰傳熱的計算則大約需要43.8 h，而考慮沸騰傳熱則需要約50.2 h。

表3 冷卻水主要參數[2]

表4 發(fā)動機主要材料參數

圖8 計算殘差曲線

圖9 上水孔質量流量的變化曲線

4 計算結果分析

4.1 沸騰對冷卻水腔的流動與壓力分布的影響

為研究沸騰傳熱對冷卻水腔內冷卻水的流動以及壓力分布的影響，取如圖10所示截面位置上的計算結果。如圖11所示，考慮沸騰與否，所示截面的速度分布和平均流速分布趨勢基本一致。如圖12所示，考慮沸騰與否，整個水腔的壓力分布結果基本一致。所以沸騰傳熱對冷卻水腔內的流動以及壓力分布影響不明顯。

圖10 缸蓋水腔截面示意圖

綜上可知，沸騰傳熱對冷卻水腔內流動與壓力分布影響不明顯。

4.2 沸騰對冷卻水腔傳熱的影響

如圖13所示，考慮沸騰與否，流固交界面上的溫度場變化很大，可見沸騰傳熱對冷卻水腔內傳熱過程的影響是很大的。

由圖14可見，考慮沸騰傳熱后計算得到的火力面最高溫度為521.7 K，比由單相對流傳熱計算得到的結果下降了8.9 K，而且1缸與4缸火力面

最高溫度的差值也從單相對流傳熱計算時的6.3 K減小到2.1 K。

圖11 考慮沸騰與不考慮沸騰傳熱計算的各個截面速度分布對比圖

圖12 冷卻水腔壓力分布云圖

圖13 機體流-固交界面溫度分布云圖

綜上可知，沸騰傳熱對冷卻水腔內傳熱過程的影響是很大的，尤其是對于熱負荷較高的如噴油器與缸蓋底板交界處、缸蓋鼻梁三角區(qū)等區(qū)域，不能忽略沸騰傳熱的影響。

5 結論

本文通過數值模擬計算，詳細分析了考慮沸騰傳熱與否對計算結果產生的影響：考慮沸騰傳熱對冷卻水腔內流動與壓力的分布影響不明顯，而對冷卻水腔內傳熱過程的影響很大。若只考慮純對流傳熱，計算結果可能與實際情況存在很大的差異。因此，在對強化內燃機進行流動與傳熱問題的研究時必須考慮沸騰傳熱的因素，以獲得更為真實、準確的結果。

The Influence of Boiling on Diesel Engine Cooling System in the Numerical Simulation

Bai Shu
（Shanghai Diesel Engine Co.,Ltd.,Shanghai 200438,China）

This study focuses on the comparison between the simulation results with or without boiling. The result of study shows that boiling mainly affects the heat transfer in cooling jacket,but does not affect the flow and pressure of the coolant in the cool jacket.The computational result is very different from the fact if convective heat transfer is considered only.So it is important to consider boiling in research in the heat transfer of internal combustion engine of high power-to-volume ratio.The will give more exact results.

internal combustion engine,boiling heat transfer,cooling water jacket, fluid-solid couple,FEA

白曙（1985-），女，碩士，主要研究方向為CFD-計算流體力學。

10.3969/j.issn.1671-0614.2013.01.001