邢鵬

摘要：內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬分析是優(yōu)化冷卻水系統(tǒng)、提高冷卻效率的重要手段。結(jié)合RPI模型為基礎(chǔ)，通過歐拉模型、沸騰氣泡行為子模型、壁面沸騰傳熱模型、湍流模型等創(chuàng)建，及時(shí)進(jìn)行試驗(yàn)研究，總結(jié)出內(nèi)燃機(jī)冷卻水腔數(shù)值波動情況，以此達(dá)到冷卻系統(tǒng)優(yōu)化的目的。

Abstract： Numerical simulation analysis of cooling water chamber of internal combustion engine is an important means to optimize cooling water system and improve cooling efficiency. Based on RPI model， Euler model， sub model of boiling bubble behavior， wall boiling heat transfer model and turbulence model were established. The numerical fluctuation of cooling water cavity of internal combustion engine was summarized to optimize the cooling system.

關(guān)鍵詞：內(nèi)燃機(jī);冷卻系統(tǒng);湍流模型;歐拉模型

Key words： internal combustion engine;cooling system;turbulence model;euler model

中圖分類號：TK124? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號：1674-957X（2020）23-0062-02

0? 引言

基于RPI模型的內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬分析，主要對冷卻系統(tǒng)在運(yùn)行過程中，缸蓋鼻梁區(qū)出現(xiàn)換熱方式不當(dāng)?shù)葐栴}加以解決，特別是內(nèi)燃機(jī)冷卻系統(tǒng)在新技術(shù)應(yīng)用基礎(chǔ)上，運(yùn)行始終處于高負(fù)荷狀態(tài)，這樣一來就需要提高換熱效率。但是在實(shí)際冷卻中，因?yàn)槭軣崃悴考L時(shí)間運(yùn)行，熱負(fù)荷變化明顯，尤其是氣缸蓋以及活塞等部位，造成內(nèi)燃機(jī)冷卻系統(tǒng)出現(xiàn)異常，運(yùn)行效率明顯下降，甚至影響到內(nèi)燃機(jī)運(yùn)行安全。在這種情況下必須加大對內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬分析，結(jié)合RPI模型，總結(jié)出更有效的處理與優(yōu)化方案。

1? 內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬分析現(xiàn)狀

內(nèi)燃機(jī)作為工業(yè)生產(chǎn)的重要設(shè)備，內(nèi)燃機(jī)運(yùn)行冷卻處理至關(guān)重要。內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬分析一直是國內(nèi)外學(xué)者研究的對象，特別是其中的內(nèi)換熱狀態(tài)相關(guān)實(shí)驗(yàn)。針對內(nèi)燃機(jī)鼻梁區(qū)以及活塞部件等方面頻繁受到高熱負(fù)荷的影響出現(xiàn)運(yùn)行異常，董非，苑天林等對此進(jìn)行核態(tài)沸騰研究，并且合理規(guī)劃核態(tài)沸騰區(qū)域，為后續(xù)冷卻處理明確方向[1]。內(nèi)燃機(jī)運(yùn)行中，缸蓋受到水腔結(jié)構(gòu)的影響，加上缸蓋本身內(nèi)部結(jié)構(gòu)復(fù)雜，封閉處理特殊，單純采用缸蓋冷卻水研究的方式對沸騰換熱問題解決無法實(shí)現(xiàn)，并且消耗較高成本?；诖宋渲緜ソY(jié)合冷卻水腔內(nèi)沸騰換熱影響展開模擬研究，利用CFD數(shù)值模擬的方式轉(zhuǎn)換出Chen模型，除此之外還得出BDL模型[2]。陳雪飛，秦國良等人就BDL模型進(jìn)行完善優(yōu)化，很大程度上將計(jì)算精度提高[3]。模型優(yōu)化中涉及大量半經(jīng)驗(yàn)公式以及數(shù)據(jù)等，不斷改善模型分析問題，打破分析局限，尤其是RPI兩相流模型的應(yīng)用，采取氣泡成核密度以及分離直徑等理論，準(zhǔn)確分析兩相流運(yùn)行狀態(tài)，通過流動狀態(tài)的總結(jié)，對內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬精確分析。根據(jù)兩相流沸騰研究，對其中氣泡形狀加以計(jì)算，對兩相流分布有效統(tǒng)計(jì)，科學(xué)改善內(nèi)燃機(jī)冷卻問題。

2? 創(chuàng)建內(nèi)燃機(jī)冷卻水腔內(nèi)沸騰傳熱數(shù)學(xué)模型

內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值分析中，必須及時(shí)創(chuàng)建沸騰傳熱數(shù)學(xué)模型，具體涉及到歐拉模型、壁面沸騰傳熱模型、沸騰氣泡行為子模型、湍流模型等。

2.1 歐拉模型分析? 基于內(nèi)燃機(jī)冷卻問題剖析，積極創(chuàng)建歐拉模型，并且根據(jù)能量方程，找到其中的連續(xù)性。其中連續(xù)性方程（1），動量方程（2），能量方程（3）。

此次計(jì)算中，涉及到速度矢量vq;相同條件下冷卻系統(tǒng)的質(zhì)量轉(zhuǎn)移pq、qp;壁面耗散力以及浮升力Fwl，q、FliFt，q;虛擬質(zhì)量力與外部體積力Fvm，q、Fq;質(zhì)量源相與應(yīng)力應(yīng)變張量Sq、q;湍流離散力與壓力、熱通量Ftd，q、p、qq;原項(xiàng)與相間焓、比定壓熱熔Sq、Qpq=-Qqp、hq;相間熱交換強(qiáng)度Qpq。

2.2 沸騰氣泡行為子模型分析? 通過對內(nèi)燃機(jī)冷卻系統(tǒng)的研究，配合歐拉模型，積極研究沸騰氣泡行為子模型，其中涉及到脫離直徑、密度、脫離頻率、等待時(shí)間等元素。其中氣泡脫離直徑模型的創(chuàng)建，主要選擇Tolubinsky研究所得的公式（4）。

計(jì)算公式中主要涉及到氣泡的參考直徑d0，取值0.0006m，同時(shí)還需要掌握飽和溫度，由此得出氣泡的脫離直徑。

氣泡成核密度模型涉及到的公式（5），主要涉及到兩個(gè)數(shù)值，n=1.805，m=210。因?yàn)槠鋵儆诮?jīng)驗(yàn)關(guān)聯(lián)式，所以在具體應(yīng)用中還要參考內(nèi)燃機(jī)冷卻具體情況。脫離頻率模型公式為（6），涉及到液相密度、氣相密度，具體為pl，pg，單位為kg/m3。

除此之外還涉及到氣泡等待時(shí)間的模型設(shè)計(jì)，主要選擇Kurl推薦公式（7）。

2.3 壁面沸騰傳熱模型分析? 對于壁面沸騰傳熱模型的創(chuàng)建，由RPI模型衍生而來，基本理論主要參照J(rèn)udd等的研究。結(jié)合內(nèi)燃機(jī)冷卻系統(tǒng)運(yùn)行中，流道受到溫度變化影響，壁面會不定期出現(xiàn)沸騰現(xiàn)象，這期間便需要計(jì)算總熱通量Qw，涉及到的計(jì)算條件主要包括蒸發(fā)項(xiàng)熱通量，由Qe表示，流項(xiàng)熱通量，由Qc表示，激冷項(xiàng)熱通量，由Qq表示。具體計(jì)算關(guān)系見式（8）、式（9）、式（10）、式（11）。

計(jì)算公式中，作為計(jì)算必備條件，必然會涉及到液相比熱容Cpl，無量綱溫度T，氣相成核密度n”，液相密度ρl，常數(shù)FA，無量綱距離Y+，液相溫度Tl取值為250，氣泡脫離頻率f，氣泡脫離直徑dw，氣泡等待時(shí)間tw，液相汽化替熱hlg，氣泡分離接觸分度Tquench。

2.4 湍流模型分析? 基于RPI模型的內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬分析中，湍流模型分析主要因?yàn)閮?nèi)燃機(jī)運(yùn)行期間，冷卻系統(tǒng)的水腔位置熱負(fù)荷較大，本身冷卻道設(shè)計(jì)便比較繁瑣，一旦冷卻系統(tǒng)運(yùn)行速度流速加快的情況，整個(gè)系統(tǒng)就會處于湍流狀態(tài)，這種情況就會影響到腔內(nèi)數(shù)值變化，因此必須創(chuàng)建湍流模型。以Realizable k-ε為基礎(chǔ)，得到內(nèi)燃機(jī)冷卻系統(tǒng)的湍流能k湍流耗散率ε，具體公式為式（12）、式（13）。

計(jì)算公式中，為了滿足計(jì)算要求，涉及到冷卻系統(tǒng)的浮力以及平均流速下冷卻運(yùn)行湍流能，分別為Gb、Gk，還涉及到源項(xiàng)Sk、Sg，常數(shù)Cg3，液相動力粘度為μ。在此基礎(chǔ)上還包含如下數(shù)值：σg=1.3、σk=1.0、μt=ρlCt、Cμ=0.09。

其中湍動能計(jì)算，公式為kd=Ck，湍流年度計(jì)算則需要根據(jù)方程（14）。

3? 內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬試驗(yàn)計(jì)算

內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬試驗(yàn)計(jì)算，首先試驗(yàn)對象為冷卻系統(tǒng)的缸蓋，及時(shí)創(chuàng)建溫度場模型，選擇最常用的增壓水冷四缸柴油機(jī)類型，在初始狀態(tài)下對缸蓋玩電腦讀測量記錄，隨后在沸騰傳熱情況下進(jìn)行仿真處理與計(jì)算，具體結(jié)構(gòu)、性能等參數(shù)詳見表1。

仿真試驗(yàn)中，結(jié)合柴油機(jī)具體情況設(shè)置硬度塞測溫點(diǎn)，數(shù)量為20個(gè)，按照1-20進(jìn)行編號，并且繪制局部實(shí)物圖。根據(jù)數(shù)值模擬要求，設(shè)計(jì)冷卻流道網(wǎng)格以及缸蓋模型，以STAR CCM+幾何模型為參考，要求網(wǎng)格尺寸必須控制到4mm，并且邊界厚度設(shè)定在1.2mm。因?yàn)閮?nèi)燃機(jī)冷卻腔內(nèi)道復(fù)雜，所以設(shè)計(jì)網(wǎng)絡(luò)數(shù)量基數(shù)較大，約83.5萬個(gè)。利用計(jì)算軟件設(shè)定模擬分析系數(shù)，并且繪制傳熱系數(shù)曲線，統(tǒng)計(jì)瞬時(shí)溫度變化，由此推斷出柴油機(jī)運(yùn)行情況，同時(shí)還能夠了解熱邊界變化，及時(shí)對模擬結(jié)果進(jìn)行分析，計(jì)算出平均環(huán)境溫度與平均傳熱系數(shù)。根據(jù)統(tǒng)計(jì)數(shù)值生成氣泡分布云圖、冷卻水腔內(nèi)流線分布云圖。及時(shí)對在內(nèi)燃機(jī)冷卻系統(tǒng)進(jìn)行優(yōu)化處理，提高冷卻效率。其中內(nèi)燃機(jī)冷卻系統(tǒng)中的缸蓋模擬數(shù)值計(jì)算，結(jié)合水腔壁面體積數(shù)值生成云圖，發(fā)現(xiàn)缸蓋的主要換熱方式為對流換熱，雖然會出現(xiàn)局部沸騰換熱情況，則主要集中在高溫區(qū)域。因?yàn)閾Q熱期間會出現(xiàn)大量氣泡，這些氣泡受到缸蓋結(jié)構(gòu)的影響會附著在壁面位置，試驗(yàn)研究中的缸蓋氣道，氣泡最密集區(qū)域?yàn)闅獾辣橇何恢?，密集度達(dá)到27.7%。結(jié)合這種情況制定流動環(huán)境改善的方案，改善氣泡聚集情況，防止因?yàn)閾Q熱惡化導(dǎo)致缸蓋氣泡出現(xiàn)膜態(tài)沸騰現(xiàn)象，進(jìn)而導(dǎo)致缸蓋出現(xiàn)裂痕。通過對冷卻水腔內(nèi)流線云圖結(jié)構(gòu)分析，發(fā)現(xiàn)氣道鼻梁位置因?yàn)榻Y(jié)構(gòu)原因，截面積相對來講比較大，所以流量比較充足，并且監(jiān)測平均流速發(fā)現(xiàn)能夠達(dá)到2m/s?？墒窃趯?shí)際流動中，因?yàn)榫o鄰排氣道，加上預(yù)熱塞對鼻梁區(qū)所造成的影響，限制了水腔壁面的流動速度，導(dǎo)致壁面的曲率出現(xiàn)異常，冷卻液流動遇到死區(qū)。針對這種現(xiàn)象必須及時(shí)對幾何結(jié)構(gòu)進(jìn)行調(diào)整，調(diào)整冷卻液平均流速，保證制冷系統(tǒng)的正常運(yùn)行。及時(shí)對比試驗(yàn)中的溫度與流動數(shù)值，特別是缸蓋的鼻梁位置，作為測量的重點(diǎn)區(qū)域，頻繁出現(xiàn)冷卻液流動異常情況，因此，需要定期對相關(guān)數(shù)值進(jìn)行計(jì)算比較，利用仿真實(shí)驗(yàn)的方式，總結(jié)出有效的改善措施，將計(jì)算誤差降到最小，提高內(nèi)燃機(jī)冷卻系統(tǒng)運(yùn)行穩(wěn)定性，并且保證冷卻系統(tǒng)數(shù)據(jù)統(tǒng)計(jì)精度。

4? 結(jié)束語

綜上所述，結(jié)合內(nèi)燃機(jī)冷卻系統(tǒng)腔內(nèi)數(shù)值模擬分析，選擇出最適合的冷沸騰換熱模型，并且進(jìn)行模擬數(shù)值分析，計(jì)算出冷卻系統(tǒng)優(yōu)化相關(guān)數(shù)值，嚴(yán)格控制模擬計(jì)算誤差。冷沸騰換熱模型對缸蓋溫度控制效果理想，并且為內(nèi)燃機(jī)冷卻系統(tǒng)優(yōu)化提供準(zhǔn)確參考。

參考文獻(xiàn)：

[1]董非，苑天林，武志偉，等.基于RPI模型的內(nèi)燃機(jī)冷卻水腔內(nèi)數(shù)值模擬研究[J].化工學(xué)報(bào)，2019，070（S2）：250-257.

[2]武志偉.冷卻介質(zhì)對內(nèi)燃機(jī)冷卻水腔內(nèi)沸騰換熱影響的模擬研究[D].2020.

[3]陳雪飛，秦國良，賈誠，等.基于氣腔模型和節(jié)流閥模型的離心壓縮機(jī)喘振數(shù)值研究（英文）[J].風(fēng)機(jī)技術(shù)，2019（5）.