柴油耦合噴嘴內(nèi)空化流動對近場噴霧影響的模擬

孫麗芝，何邦全，張更伽，張?巖

柴油耦合噴嘴內(nèi)空化流動對近場噴霧影響的模擬

孫麗芝1，何邦全1，張更伽1，張?巖2

(1. 天津大學內(nèi)燃機燃燒學國家重點實驗室，天津 300072；2. 中國北方發(fā)動機研究所，天津 300405)

柴油機共軌壓力的提高使空化現(xiàn)象快速發(fā)生而影響柴油近場噴霧的初次破碎，進而對后續(xù)燃燒和排放產(chǎn)生作用．因此搭建了單孔柴油噴油器內(nèi)外流耦合噴嘴三維數(shù)值模型，并進行了驗證．在此基礎上，模擬研究了針閥開啟初期不同噴油壓力下空化流動對近場噴霧特性的影響．模擬結(jié)果表明，噴孔內(nèi)空化發(fā)展經(jīng)歷了3個階段：單相流區(qū)(v＜0.3%)、空化發(fā)展區(qū)(0.3%＜v＜21%)及超空化區(qū)(v＞21%)，且空化發(fā)展中會伴隨云空化的產(chǎn)生．噴油壓力越大，空化初生時間越早，發(fā)展速度越快．進入超空化區(qū)后，空化氣泡潰滅最先引起近場噴霧上游的主液柱區(qū)的破碎，破碎區(qū)逐步向噴霧下游延伸．空化引起噴霧破碎粒徑減小，液滴數(shù)增加．空化使得噴霧徑向速度增大，徑向擴展能力增強，使得噴霧錐角增大．空化區(qū)渦值大，旋渦結(jié)構小且旋轉(zhuǎn)強度大，進入噴射腔后增強了對噴霧的擾動，使得噴霧液芯變細，促進了液芯外側(cè)的噴霧破碎．

噴嘴；空化流動；近場噴霧；破碎機理；數(shù)值模擬

柴油機的燃燒和排放特性與燃燒室內(nèi)形成的混合氣質(zhì)量息息相關，而影響混合氣質(zhì)量的一個重要因素就是噴霧特性[1-3]．噴霧特性受多種因素的影響，其中，噴嘴內(nèi)部的空化對近場噴霧的發(fā)展起著重要作用．目前國內(nèi)外學者就空化對近場噴霧的影響進行了諸多研究．He等[4]在噴油壓力為1.0MPa、背壓為0.1MPa的條件下，采用透明放大噴嘴進行可視化試驗，研究了不同噴孔形狀產(chǎn)生的空化對近場噴霧的影響，發(fā)現(xiàn)擴張型噴孔更易產(chǎn)生空化，空化又可促進近場噴霧錐角的增大．Desantes等[5]采用可視化方法試驗研究了噴油壓力為50MPa、背壓為2.5MPa時空化現(xiàn)象對噴霧初次破碎的影響，證明了噴霧錐角的增大與空化氣泡的潰滅有關．Kima等[6]采用透明噴嘴在背壓為0.2MPa下進行可視化試驗，研究了空化長度對噴霧錐角的影響，發(fā)現(xiàn)空化會造成噴霧錐角增大．Suh等[7]在背壓為大氣壓、噴油壓力為0.13～0.45MPa時進行可視化試驗，對不同長徑比下的孔內(nèi)空化及近場噴霧特性進行研究，結(jié)果表明，長徑比越小，空化現(xiàn)象越劇烈，噴霧霧化效果越好．Le?nik等[8]通過模擬計算孔內(nèi)空化與高速攝影試驗拍攝噴霧圖像相結(jié)合的方法，對比研究了柴油與生物柴油的空化程度和近場噴霧特性，結(jié)果表明，全負荷時柴油與生物柴油形成的空化現(xiàn)象差別不大．Lai等[9]采用模擬計算與試驗相結(jié)合的方法，研究了單孔和多孔噴油器的噴孔內(nèi)部流動及近場噴霧特性，發(fā)現(xiàn)多孔噴嘴內(nèi)部流動更劇烈，近場噴霧結(jié)構更不穩(wěn)定．Ghorbani?等[10]將噴孔近似簡化為微孔道，通過模擬計算微孔道內(nèi)部空化現(xiàn)象，再由高速攝影得到近場噴霧形態(tài)的方法，研究了大氣條件下噴油壓力對空化及近場噴霧的影響，研究發(fā)現(xiàn)，當噴油壓力超過5MPa后，液體射流轉(zhuǎn)變?yōu)樵旗F狀射流，且噴霧錐角隨噴油壓力的增大而增大．Hayashi等[11]亦采用模擬計算與可視化試驗相結(jié)合的方法研究了柴油噴嘴內(nèi)的空化流動對噴霧及燃燒的影響，結(jié)果表明，弦空化的發(fā)生可以增大噴霧錐角，增強噴霧破碎，從而提高混合氣質(zhì)量，優(yōu)化燃燒.因此，研究空化對近場噴霧的影響機理至關重要.

隨著車輛排放法規(guī)與使用需求的提升，對柴油機共軌壓力提出了更高的要求，在噴孔內(nèi)更易產(chǎn)生空化現(xiàn)象，尤其在噴油初期，進而影響近場噴霧[12]．目前國內(nèi)外仍主要集中在低壓條件下采用可視化試驗或孔內(nèi)模擬與孔外試驗相結(jié)合的方法對空化及近場噴霧特性進行研究，一方面對實際發(fā)動機運行工況的研究較少，另一方面孔內(nèi)外采用兩種研究方法增大了數(shù)據(jù)誤差，且孔外近場噴霧均采用試驗方式獲得，在空化對近場噴霧的影響機理研究方面存在局限．本文將建立柴油機單孔噴孔內(nèi)外耦合模型，在更接近發(fā)動機實際的高噴油壓力、高背壓的情況下，對針閥開啟初期、不同噴油壓力下的孔內(nèi)空化流動及近場噴霧進行同時模擬研究，以揭示噴油壓力對不同噴油時刻的瞬時空化發(fā)展狀態(tài)及對近場噴霧特性的影響機理．值得注意的是，綜合上述研究現(xiàn)狀，目前研究中常用噴嘴長徑比為3～5，但本文重點是內(nèi)外流耦合的模擬研究，所以長徑比不予考慮，按試驗室用真實噴嘴尺寸建模．

1?計算模型的構建

1.1?數(shù)學模型

柴油機噴孔尺寸極小，在高壓下孔內(nèi)燃油高速流動，因此采用RNG模型[13]描述噴油瞬態(tài)湍流運動過程．其中，湍動能的輸運方程為

關于湍流耗散率的輸運方程為

式中：C＝0.0845；Pr＝1.39；Pr＝1.39；C1＝1.55；C2＝1.68；C3＝－1.0；＝0.012；＝4.38．

燃油在噴嘴內(nèi)流動采用VOF模型，同時啟用基于Shields等[14]閃急沸騰理論的空化模型來計算孔內(nèi)空化現(xiàn)象．噴孔內(nèi)外的耦合計算采用ELSA模型[15]，模型使用兩個條件作為歐拉算法向拉格朗日算法的過渡標準，一是網(wǎng)格中的氣體體積分數(shù)，采用空隙率進行控制，二是網(wǎng)格中的液體表面積，采用比表面積進行控制．當網(wǎng)格空隙率小于0.8、液體燃油比表面積小于1.2時，網(wǎng)格中的燃油轉(zhuǎn)換成噴霧液滴，歐拉算法過渡為拉格朗日算法，其輸運方程為

本文采用KH-RT噴霧破碎模型[16-17]來描述湍流效應及空氣阻力引起的不穩(wěn)定波對近場噴霧的影響，采用NTC Collision模型[18]來計算噴霧過程中的油滴碰撞．

1.2?幾何模型

根據(jù)試驗室真實噴嘴尺寸圖建立了耦合噴嘴的三維幾何結(jié)構圖，如圖1所示，噴孔直徑()為0.16mm，孔長()為0.96mm．

圖1?噴嘴三維計算模型

為保證計算結(jié)果的準確性，對模型網(wǎng)格進行了合理設計．模型基礎網(wǎng)格大小為80μm．為捕捉針閥運動對流動的影響，針閥壁網(wǎng)格加密到20μm．壓力室和噴孔是燃油流動的重要邊界，網(wǎng)格分別加密到40μm和10μm．同時噴霧破碎也需精確計算，噴射腔網(wǎng)格加密到10μm．燃油物性及模型邊界參數(shù)設置如表1所示．

表1?燃油物性及邊界參數(shù)

Tab.1?Fuel physical properties and boundary parameters

2?試驗驗證

為保證模擬計算結(jié)果的準確性，對孔內(nèi)流動的空化模型及噴霧模型進行了驗證．

2.1?空化模型驗證

受試驗設備的限制，無法得到孔內(nèi)空化試驗圖像，因此采用Winklhofer等[19]在U型噴嘴中得到的空化圖像和質(zhì)量流量試驗數(shù)據(jù)進行空化模型的驗證．Winklhofer等人所做的空化試驗數(shù)據(jù)被公認為空化模型驗證的標準，其已廣泛應用于空化模型的驗證．圖2為噴孔質(zhì)量流量與空化形態(tài)的試驗與仿真結(jié)果對比圖．其中，藍色部分為氣相，表示空化區(qū)域；紅色部分為液相．由圖2(a)可見，不同壓差下噴孔質(zhì)量流量的模擬結(jié)果與試驗結(jié)果接近，且變化趨勢相同，質(zhì)量流量均隨壓差的提高而增大，當達到臨界空化壓差7MPa后，繼續(xù)增大壓差，質(zhì)量流量不再發(fā)生變化．由圖2(b)可見，不同壓差下孔內(nèi)空化的試驗與模擬現(xiàn)象接近，但存在細小差別，這是因為模擬所用燃料C12H26與試驗所用柴油理化性質(zhì)有差異，且模擬計算中假設壁面光滑，未考慮噴孔壁面的粗糙度.可見，本研究選用的空化模型可以正確描述噴孔內(nèi)的空化現(xiàn)象．

圖2?實驗與模擬結(jié)果對比

2.2?噴霧模型驗證

在實驗室的試驗臺架上得到了噴油壓力為90MPa、環(huán)境壓力為4.5MPa條件下的孔外近場噴霧試驗數(shù)據(jù)，以進一步地驗證耦合噴嘴模型及噴霧模型的準確性．圖3給出了模擬與試驗數(shù)據(jù)結(jié)果的對比圖．圖中現(xiàn)象顯示，在不同噴油后時刻(after start of injection，ASOI)，模擬與試驗得到的近場噴霧破碎長度與噴霧貫穿距接近，且噴霧破碎位置相同，這說明選用的噴霧模型可以較好地描述近場噴霧．模擬與試驗噴霧形狀存在細微差別是兩者所用燃油物性參數(shù)差異和定容彈內(nèi)空氣運動造成的．

圖3?近場噴霧的試驗與模擬圖像

3?結(jié)果與分析

3.1?不同噴油壓力下的空化流動

圖4給出了針閥開啟初期，各噴油壓力下孔內(nèi)空化程度隨ASOI的發(fā)展規(guī)律．其中，v為噴孔內(nèi)氣相體積分數(shù)，即氣相體積與噴孔總體積之比，藍色區(qū)域表示液相燃油，其余顏色為空化區(qū)域．圖中數(shù)據(jù)顯示，在ASOI為9μs之前，針閥剛剛開啟，孔內(nèi)燃油流速較低，各噴油壓力下噴孔內(nèi)均無空化現(xiàn)象發(fā)生．隨著針閥抬起，壓力室內(nèi)燃油壓力升高，孔內(nèi)燃油流速增加，在噴孔入口處發(fā)生流動分離，形成局部回流區(qū)．當回流區(qū)壓力低于燃油飽和蒸氣壓時，空化發(fā)生，孔內(nèi)氣相體積分數(shù)不斷增大．空化層隨燃油流動向噴孔出口移動，直至由噴孔出口流出，形成超空化現(xiàn)象．由于燃油的黏性阻力，空化層厚度由噴孔入口向噴孔出口不斷增加．依據(jù)圖5噴孔內(nèi)的燃油流動狀態(tài)，空化發(fā)展可分為3個時期：噴孔內(nèi)無空化現(xiàn)象的單相流區(qū)(v＜0.3%)、噴孔內(nèi)出現(xiàn)空化但空化未達到噴孔出口的空化發(fā)展區(qū)(v為0.3%～21%)、空化長度達到噴孔長度的超空化區(qū)(v＞21%)．此外，噴油壓力越大，氣相體積分數(shù)增加速率越快，空化在噴孔內(nèi)發(fā)展速度越快．在ASOI為9μs下，噴油壓力為120MPa時，孔內(nèi)無空化現(xiàn)象，處于單相流區(qū)；但噴油壓力為150MPa、180MPa與210MPa時，空化處于空化發(fā)展區(qū)．在ASOI為12μs時，各噴油壓力下的孔內(nèi)空化均處于空化發(fā)展區(qū)，噴油壓力越大，空化區(qū)域越大．ASOI為15μs時，在120MPa與150MPa噴油壓力條件下孔內(nèi)仍處于空化發(fā)展階段，而在噴油壓力為180MPa與210MPa時，燃油空化區(qū)已流出噴孔，進入超空化區(qū)．

圖5為噴油壓力為210MPa，噴油初期噴孔內(nèi)燃油的流動狀態(tài)變化．其中，黑線為燃油流動軌跡線．在ASOI為1μs時，針閥開啟瞬間，壓力室內(nèi)產(chǎn)生壓力降，出現(xiàn)了孔外空氣流入噴孔內(nèi)的現(xiàn)象，發(fā)生倒吸．在ASOI為9μs時，高壓燃油進入壓力室，噴孔上下為正壓梯度，噴孔內(nèi)產(chǎn)生空化現(xiàn)象，空化區(qū)顯示為明顯的回流區(qū)．在ASOI為13μs時，空化區(qū)超過噴孔長度的一半，此時由于燃油氣液界面存在速度差，出現(xiàn)了Kelvin-Helmholtz不穩(wěn)定現(xiàn)象，液相流對空化層存在剪切作用，使得空化界面產(chǎn)生波動，引起空化表面的氣泡破碎為細小結(jié)構，且空化尾部的回射流加速了空化脫落，在噴孔中形成氣液兩相流的云空化[20]，并隨燃油向噴孔下游移動．當ASOI為20μs時，空化現(xiàn)象占據(jù)整個噴孔長度，即進入超空化區(qū)，空化隨燃油流出噴孔．

圖4?噴油初期，不同噴油壓力下孔內(nèi)空化的發(fā)展

圖5 噴油壓力為210MPa時，噴油初期的孔內(nèi)燃油流動

3.2?空化對近場噴霧特性的影響

為了明確空化對近場噴霧形態(tài)的影響，給出了不同噴油壓力下距噴孔出口1mm內(nèi)平面的噴霧時空變化圖像，如圖6所示．為便于識別，根據(jù)孔內(nèi)氣相體積分數(shù)大小進行著色表示空化；空化潰滅使得噴霧破碎，破碎后的液滴脫離主射流，徑向運動趨勢增加，因此使用徑向速度V對噴霧圖像進行著色，隨著噴油后時間的變化即可以觀察破碎位置的動態(tài)變化，同時可以觀察粒子的運動方向．圖中現(xiàn)象顯示，ASOI為9μs時，在各噴油壓力下，孔內(nèi)均處于單相流區(qū)，無空化現(xiàn)象．燃油噴出噴孔后受周圍氣體的剪切阻力作用，形成傘狀頭部，頭部徑向速度V值增大，剪切邊界層外層噴霧在空氣卷吸的作用下破碎并形成細小的液滴群．隨噴油壓力增大，近場噴霧貫穿距加大．隨著噴油進行，ASOI為12μs時，不同的噴油壓力下孔內(nèi)出現(xiàn)不同程度的空化現(xiàn)象，均進入空化發(fā)展區(qū)，近場噴霧形態(tài)與上一時刻相同．ASOI為16μs時，噴油壓力inj≥150MPa的條件下，孔內(nèi)燃油流動進入超空化區(qū)．孔內(nèi)的空化氣泡隨燃油噴出噴孔，在高壓環(huán)境中潰滅，促進噴霧上游主液柱破碎．與上一時刻相比，近噴孔處噴霧徑向速度V增大，向兩側(cè)運動趨勢增大，產(chǎn)生破碎，使得主液柱變寬．值得注意的是，噴油壓力為120MPa時，孔內(nèi)未達到超空化流態(tài)，但孔內(nèi)云空化的發(fā)生使得近噴孔處噴霧主液柱出現(xiàn)了破碎現(xiàn)象．在ASOI為20μs時，各噴油壓力下孔內(nèi)空化均達到超空化流態(tài)，噴油壓力120MPa時，破碎集中在近場噴霧主液柱中上部，而噴油壓力210MPa時，整個噴霧貫穿長度均發(fā)生破碎，說明噴油壓力越大，噴霧破碎區(qū)域明顯增大，這主要是由于相同時刻，噴油壓力越大，空化氣泡潰滅提前，且高速噴霧攜帶空化氣泡運動更遠距離的原因．此外，當噴油壓力一定時，隨噴油時間的增加，空化氣泡隨近場噴霧向下運動，近場噴霧主射流破碎區(qū)域也逐步向下游發(fā)展．

圖6?不同噴油壓力下近場噴霧的徑向速度時空變化

為了進一步研究空化對破碎液滴大小的影響，提取了12μs、20μs時刻，各噴油壓力下距噴孔出口0.5mm截面處不同粒徑破碎液滴的數(shù)量，如圖7所示．結(jié)合圖6中的現(xiàn)象表明，在ASOI＝12μs時，孔內(nèi)流動處于空化發(fā)展區(qū)，空化未對近場噴霧產(chǎn)生影響，噴霧在空氣動力學作用下破碎，破碎粒子數(shù)隨粒子半徑呈現(xiàn)中間高、兩頭低的趨勢．破碎粒子的半徑集中在小于200μm的區(qū)域內(nèi)．噴油壓力越大，各粒子半徑下的液滴數(shù)越多．此外，噴油壓力分別為210MPa、180MPa、150MPa與120MPa時，破碎粒子數(shù)量的最大值分別對應半徑約為67.4μm、76.6μm、89.6μm和99.6μm，這是由于噴油壓力越大，噴霧與空氣的相對速度越大，環(huán)境氣體對燃油擾動增加，破碎粒子半徑減?。擜SOI為20μs時，孔內(nèi)燃油流動處于超空化區(qū)，空化對噴霧產(chǎn)生影響，相比于圖7(a)中ASOI為12μs，相同噴油壓力下破碎粒子數(shù)均約增加了兩倍，且各噴油壓力下破碎粒子的半徑減小，集中在粒徑小于70μm的區(qū)域．噴油壓力分別為210MPa、180MPa、150MPa與120MPa時，破碎粒子數(shù)量的最大值分別約在23.4μm、26.5μm、31.8μm、34.3μm處，較圖7(a)的粒子半徑的差異大大減小，這說明空化作用不僅增加了噴霧破碎數(shù)目，還增強了噴霧破碎程度．

圖7 距噴孔出口0.5mm截面處不同粒徑下的破碎粒子數(shù)

圖8給出了噴油初期的近場噴霧錐角變化圖，其中，錐角是對距噴孔出口6倍噴孔直徑長度的噴霧外輪廓進行最小二乘法擬合計算得到的．圖中數(shù)據(jù)顯示，噴油初期，噴油壓力越大，燃油噴出噴孔的時間越早；隨后各噴油壓力下的噴霧錐角均隨ASOI增大，達到峰值后，錐角趨于穩(wěn)定；噴油壓力越大，近場噴霧錐角增加的速度越快，同時刻錐角越大．以上現(xiàn)象主要是因為噴油壓力越大，噴孔出口速度越高，高壓下產(chǎn)生的高速射流與空氣的相互作用更強，促進了燃油破碎，使得噴霧錐角增大．此外，噴油壓力在210MPa、180MPa、150MPa與120MPa時，近場噴霧錐角達到峰值的時間分別為14μs、15μs、16μs和17μs，這說明噴油壓力越大，噴霧錐角峰值出現(xiàn)的時間越早．這是由于噴油壓力增大，一方面由于空氣動力學的作用增加了噴霧破碎，另一方面噴孔內(nèi)到達超空化所需時間越短，空化破碎引起近場噴霧主射流液相區(qū)發(fā)生破碎，增強了霧化．超空化現(xiàn)象穩(wěn)定后，噴霧錐角也逐漸趨于穩(wěn)定．

圖8?近場噴霧錐角隨ASOI的變化

圖9為噴油壓力在210MPa、180MPa、150MPa與120MPa時發(fā)生超空化時刻(14μs、15μs、16μs和17μs)的噴孔出口截面的徑向速度r與噴孔軸向速度V的比值沿噴孔直徑的分布．因為噴孔邊界速度矢量的方向?qū)婌F錐角有重要影響[21]，故r/V的比值大小可表征噴霧徑向擴展的能力，比值越大，噴霧擴展能力越強．圖中r/V的速度比值曲線說明，在ASOI為14μs時，噴油壓力為210MPa時，噴孔出口r/V呈現(xiàn)出中間低、邊緣高的特點，這是因為此時噴孔內(nèi)已達到超空化，噴孔中心為液相燃油，而噴孔邊界為空化區(qū)，空化區(qū)湍流強度大，渦動量大．噴油壓力為180MPa、150MPa與120MPa時，分別在15μs、16μs和17μs時刻達到超空化流態(tài)，r/V依次增大，噴霧徑向擴展能力增強，噴霧錐角增大．由此可見，空化現(xiàn)象對近場噴霧錐角的增大有促進作用．

圖9?噴孔出口Vr/Vz沿噴孔直徑的分布

3.3?空化對近場噴霧破碎的作用機理

氣液界面存在速度及黏度差，加之噴孔壁面的剪切作用，空化區(qū)的渦值增強[22]．因此對空化區(qū)的渦量值進行提取，噴油壓力為150MPa時，發(fā)現(xiàn)空化產(chǎn)生的渦值不低于1.1×107s-1，因此以這個值作為150MPa時空化產(chǎn)生的渦區(qū)邊界值，空化產(chǎn)生的渦區(qū)域定義為強渦區(qū)．圖10給出了噴油壓力為150MPa時，不同空化狀態(tài)下的近場噴霧破碎機理示意圖．其中，液芯是根據(jù)ELSA模型的過渡算法計算得到的液相噴霧．圖中現(xiàn)象顯示，在ASOI為12μs時，孔內(nèi)處于空化發(fā)展區(qū)，空化產(chǎn)生的強渦區(qū)只存在噴孔內(nèi)部，未對近場噴霧產(chǎn)生影響，液芯頭部受空氣阻力向上卷曲，液芯外側(cè)的破碎液滴群是由于空氣的強剪切和卷吸作用形成的．當ASOI為20μs時，孔內(nèi)到達超空化，強渦區(qū)隨著空化向近場噴霧中發(fā)展，造成噴霧主射流的破碎，液芯變細，且徑向速度V迅速增大，外側(cè)噴霧破碎增強．這是因為空化潰滅產(chǎn)生的強渦區(qū)對噴霧表面擾動增大，增強了噴霧主射流表面的不穩(wěn)定性，促進噴霧破碎．由此可見，超空化前，噴霧破碎主要受空氣動力學作用；孔內(nèi)達到超空化流態(tài)后，噴霧破碎受空化和空氣動力學的共同作用，其中，空化主導噴霧上游主射流的破碎．

圖10 噴油壓力為150MPa時，不同空化狀態(tài)時的近場噴霧破碎機理

圖11 以速度值著色的Q等值面(Q＝4e10)表示的空化及噴霧旋渦結(jié)構

4?結(jié)?論

建立了單孔柴油噴孔內(nèi)流與近場噴霧耦合計算模型，并用試驗數(shù)據(jù)進行了驗證．研究得出以下結(jié)論.

(1) 噴油壓力從120MPa提高到210MPa，噴孔內(nèi)空化初生時間提前．空化發(fā)展區(qū)，在噴油后相同時刻，噴油壓力越大，空化程度越強．進入超空化區(qū)后，空化程度不再變化．孔內(nèi)空化在不穩(wěn)定波及尾部回射流的作用下會脫落形成云空化．

(2) 孔內(nèi)達到超空化流態(tài)后，空化氣泡進入噴射腔引起噴霧主射流的破碎．噴孔外相同位置處，空化使得破碎液滴的半徑減小，且破碎液滴的數(shù)目大大增加．空化流使得近場噴霧破碎后徑向擴展能力增強，從而造成噴霧錐角增大．相同時刻，噴油壓力越大，噴霧破碎效果越好，噴霧錐角達到峰值的時刻越早且峰值越大．

(3) 空化區(qū)渦值大，旋渦結(jié)構小且旋轉(zhuǎn)強度大．超空化后，空化進入噴射腔，對噴霧的擾動增大，使得噴霧液芯變細，外側(cè)噴霧破碎增強．噴油后相同時刻，高噴油壓力下的噴油速度越高，空化區(qū)旋渦的旋轉(zhuǎn)強度越大，近場噴霧主液柱的破碎區(qū)域越大．

［1］ 徐?杰，張?力. 不同噴孔入口的柴油機噴嘴空化特性[J]. 燃燒科學與技術，2019，25(4)：304-308.

Xu Jie，Zhang Li. Cavitation characteristics of diesel engine nozzle with different nozzle entrances[J].，2019，25(4)：304-308(in Chinese).

［2］ 何志霞，鐘汶君，黃云龍，等. 針閥運動對柴油機噴嘴瞬態(tài)流動特性的影響[J]. 內(nèi)燃機學報，2012，30(4)：336-342.

He Zhixia，Zhong Wenjun，Huang Yunlong，et al. Investigation of transient behavior of cavitation flow in injector nozzles affected by the needle movement[J].，2012，30(4)：336-342(in Chinese).

［3］ 何邦全，張倓愷. 直噴汽油噴油器噴射過程近場噴霧特性[J]. 內(nèi)燃機學報，2013，31(5)：426-430.

He Bangquan，Zhang Tankai. Near-field characterization of direct injection gasoline sprays from a single-hole injector[J].，2013，31(5)：426-430(in Chinese).

［4］ He Zhixia，Guo Genmiao，Tao Xicheng，et al. Study of the effect of nozzle hole shape on internal flow and spray characteristics[J].，2016，71：1-8.

［5］ Desantes J，Payri R，Salvador F，et al. Influence of cavitation phenomenon on primary break-up and spray behavior at stationary conditions[J].，2010，89：3033-3041.

［6］ Kima B，Park S. Study on in-nozzle flow and spray behavior characteristics under various needle positions and length-to-width ratios of nozzle orifice using a transparent acrylic nozzle[J].，2019，143：118478.

［7］ Suh H，Lee C. Effect of cavitation in nozzle orifice on the diesel fuel atomization characteristics[J].，2008，29：1001-1009.

［8］ Le?nik L，Kegl B，Bombek G. The influence of in-nozzle cavitation on flow characteristics and spray break-up[J].，2018，222：550-560.

［9］ Lai Mingchia，Zheng Yi，Xie Xingbin，et al. Characterization of the near-field spray and internal flow of single-hole and multi-hole sac nozzles using phase contrast X-ray imaging and CFD[J].，2011，4(1)：703-719.

［10］ Ghorbani M，Sadaghiani A，Yidiz M，et al. Experimental and numerical investigations on spray structure under the effect of cavitation phenomenon in a micro-channel[J].，2017，31(1)：235-247.

［11］ Hayashi Tomohiro，Suzuki Masayuki，Ikemoto Masato. Effects of internal flow in a diesel nozzle on spray combustion[J].，2013，14(6)：646-654.

［12］ Ghiji M，Goldsworthy L，Garaniya V，et al. Effect of residual air bubbles on diesel spray structure at the start of injection[J].，2019，241：25-32.

［13］ Han Z，Reitz R D. Turbulence modeling of internal combustion engines using RNG κ-ε models[J].，1995，106(4/5/6)：267-295.

［14］ Shields B，Neroorkar K，SchmidtD P. Cavitation as rapid flash boiling[C] //. Venture，CA，2011.

［15］ Blanco P. Implementation and Development of an Eulerian Spray Model for CFD Simulations of Diesel Sprays[D]. València，Spain：Universitat Politècnica de València，2016.

［16］ Reitz R，Bracco F. Mechanisms of breakup of round liquid jets[J].，1986，3：233-249.

［17］ Xin J，Ricart L，Reitz R. Computer modeling of diesel spray atomization and combustion[J].，1998，137(1/2/3/4/5/6)：171-194.

［18］ Schmidt D，Rutland C. A new droplet collision algorithm[J].，2000，164(1)：62-80.

［19］ Winklhofer E，Kull E，Kelz E，et al. Comprehensive hydraulic and flow field documentation in model throttle experiments under cavitation conditions[C]//. Zurich，Swizerland，2001：574-579.

［20］ 潘森森，彭曉星. 空化機理[M]. 北京：國防工業(yè)出版社，2013：73-81.

Pan Sensen，Peng Xiaoxing.[M]. Beijing：National Defense Industry Press，2013：73-81(in Chinese).

［21］ Dong Pengbo，Inaba Takuya，Nishida Keiya，et al. Characteristics of the internal flow and the near-field spray of a single-hole injector and a multi-hole injector for diesel engines[J].，：，2016，230(5)：632-649.

［22］ Yu H，Goldsworthy L，Brandner P，et al. Modelling thermal effects in cavitating high-pressure diesel sprays using an improved compressible multiphase approach[J].，2018，222：125-145.

［23］ Hunt J C R，Wray A A，Moin P. Eddies，stream，and convergence zones in turbulent flows[C]//. Stanford，USA，1988：193-208.

［24］ Fu Wu-Shung，Lai Yu-Chih，Li Chung-Gang. Estimation of turbulent natural convection in horizontal parallel plates by the Q criterion[J].，2013，45：41-46.

Simulation of Effect of Cavitation Flow in a Diesel Injector Nozzle on Near-Field Spray

Sun Lizhi1，He Bangquan1，Zhang Gengjia1，Zhang Yan2

(1. State Key Laboratory of Engines，Tianjin University，Tianjin 300072，China；2. China North Engine Research Institute，Tianjin 300405，China)

With the increase of common rail pressure of a diesel engine，cavitation occurs rapidly，which affects the primary breakup of near-field spray and further affects combustion and emission. A coupled three-dimensional numerical model of single-hole diesel injector is established and verified. The effect of cavitation flow on near-field spray under different injection pressures at the beginning of needle valve opening is studied. The simulation results show that cavitation can be divided into three regions: single flow zone (v＜0.3%)，cavitation development zone (0.3%＜v＜21%) and super-cavitation zone (v＞21%). The flow of cavitation is accompanied by the occurrence of cloud cavitation. The higher the injection pressure，the earlier the cavitation initiation time and the faster the development speed. After the super-cavitation，the breakup occurs first in the main stream of the upstream of the near-field spray due to cavitation bubbles，and then the breakup area gradually extends to the downstream. Cavitation leads to decreased atomized particle size and increased droplet number. It also increases the radial velocity of spray and the tendency of radial spreading，thus increasing the spray cone. Cavitation area shows such characteristics as large vortex value，small vortex structure and large rotation intensity. When cavitation enters the injector cavity，the disturbance to the spray is enhanced，thus resulting in the thinner spray core and promoted break-up on the outside of the spray core.

nozzle；cavitation flow；near-field spray；breakup mechanism；numerical simulation

TK421

A

1006-8740(2021)01-0081-09

10.11715/rskxjs.R202003001

2020-03-02．

國家自然科學基金資助項目(51606175；51476151)．

孫麗芝（1994—??），女，碩士，lzsun@tju.edu.cn.

何邦全，男，博士，副教授，bqhe@tju.edu.cn.