李 瀟，王佔(zhàn)元，劉宗寬，周 磊，衛(wèi)海橋

（天津大學(xué)內(nèi)燃機(jī)燃燒學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津 300072）

0 概述

隨著環(huán)境污染問(wèn)題的日益嚴(yán)重和化石能源短缺問(wèn)題的不斷加劇，內(nèi)燃機(jī)領(lǐng)域的排放法規(guī)也逐漸嚴(yán)格。優(yōu)化燃燒過(guò)程和提高內(nèi)燃機(jī)的熱效率是降低污染物排放的有效手段，而燃油的霧化及混合氣的形成直接影響燃燒的進(jìn)行，對(duì)燃油在缸內(nèi)的擴(kuò)散混合過(guò)程進(jìn)行研究有重要意義。

在現(xiàn)代發(fā)動(dòng)機(jī)上，先進(jìn)燃燒模式和增壓技術(shù)的使用提高了熱效率，也使得缸內(nèi)溫度和壓力不斷增高。在某些高增壓發(fā)動(dòng)機(jī)上，燃油在噴射初始時(shí)刻缸內(nèi)壓力和溫度可達(dá)2.5 MPa 和700 K，已經(jīng)超過(guò)大部分烴類燃油的臨界壓力和溫度，例如正十二烷的臨界壓力和溫度分別為1.82 MPa 和658 K。由此可見(jiàn)，缸內(nèi)壓力和溫度超過(guò)燃料臨界點(diǎn)的狀況會(huì)大概率發(fā)生，對(duì)超臨界環(huán)境下的噴霧特性進(jìn)行研究顯得尤為重要。

超臨界流體具有獨(dú)特性質(zhì)，文獻(xiàn)［1］中通過(guò)試驗(yàn)證明超臨界流體密度類似液體，黏度類似氣體，而擴(kuò)散性質(zhì)介于液體和氣體之間。由于超臨界流體的特殊物性，燃油在超臨界環(huán)境下的噴霧也表現(xiàn)出不同于亞臨界環(huán)境的特有屬性。亞臨界環(huán)境下由于表面張力的存在，射流與周圍氣體之間存在明確的氣液界面，而在超臨界環(huán)境中的噴霧過(guò)程中，表面張力趨近于零導(dǎo)致氣液界面消失，擴(kuò)散過(guò)程靠湍流擴(kuò)散作用控制，射流表面各參數(shù)存在較大的梯度變化。關(guān)于對(duì)超臨界流體熱物性的預(yù)測(cè)，文獻(xiàn)［2］中指出超臨界流體分子間的相互作用不能被忽略，需使用真實(shí)流體狀態(tài)方程。超臨界流體的輸運(yùn)系數(shù)（如黏度系數(shù)和導(dǎo)熱系數(shù)）的計(jì)算也要考慮流體超臨界狀態(tài)的影響。

超臨界噴霧的試驗(yàn)研究已經(jīng)取得了較大進(jìn)展，文獻(xiàn)［3］中通過(guò)試驗(yàn)發(fā)現(xiàn)在超臨界環(huán)境下射流表面不會(huì)產(chǎn)生亞臨界環(huán)境下出現(xiàn)的液滴和液絲，而由液態(tài)條狀組成的梳子狀結(jié)構(gòu)代替。文獻(xiàn)［4］中利用平面激光誘導(dǎo)熒光（planar laser induced fluoresence，PLIF）技術(shù)，實(shí)現(xiàn)了在熱力學(xué)大范圍內(nèi)對(duì)流體密度和密度梯度的測(cè)量。許多學(xué)者也開(kāi)展了關(guān)于超臨界噴霧的數(shù)值模擬研究。文獻(xiàn)［5］中基于大渦模擬，利用描述真實(shí)流體行為的計(jì)算模型，對(duì)正庚烷噴入超臨界環(huán)境下的行為進(jìn)行了分析，發(fā)現(xiàn)密度梯度在噴霧過(guò)程中起著重要作用。文獻(xiàn)［6］中研究了不同超臨界壓力對(duì)低溫液氮噴霧的影響，指出超臨界下產(chǎn)生的密度梯度層會(huì)抑制射流徑向速度的振蕩，而促進(jìn)軸向速度的振蕩，且環(huán)境壓力的增加將導(dǎo)致射流提前過(guò)渡到自相似狀態(tài)。文獻(xiàn)［7］中利用直接數(shù)值仿真方法（direct numerical simulation，DNS）對(duì)超臨界射流中混合層的發(fā)展進(jìn)行了分析，證明密度梯度對(duì)混合層的穩(wěn)定性和湍流特性起關(guān)鍵作用。

從對(duì)超臨界噴霧的研究來(lái)看，目前的試驗(yàn)研究證明了在超臨界環(huán)境下噴霧形態(tài)的特殊性。由于超臨界流體的特殊性，亞臨界環(huán)境下的數(shù)值模擬策略不再適用，超臨界環(huán)境下的噴霧模型構(gòu)建具有挑戰(zhàn)性。關(guān)于密度梯度在超臨界噴霧中的作用原理及超臨界環(huán)境下射流擴(kuò)散演化過(guò)程的影響因素研究更是少見(jiàn)。

本研究基于上述超臨界噴霧的研究現(xiàn)狀，以開(kāi)源軟件OpenFoam 為計(jì)算平臺(tái)，構(gòu)建了適合超臨界環(huán)境的噴霧模型，分析了超臨界條件下渦的演化過(guò)程。模型采用真實(shí)流體的狀態(tài)方程及熱物性和輸運(yùn)系數(shù)的計(jì)算方法，結(jié)合大渦模擬通過(guò)對(duì)比超臨界和亞臨界環(huán)境下噴霧過(guò)程，重點(diǎn)研究了超臨界和亞臨界環(huán)境下燃油擴(kuò)散混合過(guò)程的區(qū)別與聯(lián)系，從根本上揭示了超臨界環(huán)境下射流擴(kuò)散的特殊性質(zhì)。

1 理論模型

1.1 大渦模擬基本方程

本文中采用大渦模擬的方法對(duì)超臨界噴霧過(guò)程進(jìn)行了模擬，大渦模擬是將流場(chǎng)在空間上進(jìn)行過(guò)濾，通過(guò)濾波函數(shù)使大渦和小渦分離，大尺度渦可以直接通過(guò)數(shù)值模擬求解，而小尺度渦（亞網(wǎng)格尺度）需要亞網(wǎng)格模型?；蠼?，因此大渦模擬能捕捉到雷諾平均無(wú)法捕捉到的流場(chǎng)細(xì)微結(jié)構(gòu)，得到的計(jì)算結(jié)果更接近真實(shí)噴霧形態(tài)。

控制方程為可壓縮形式的Navier-Stokes（NS）方程，包括連續(xù)方程、動(dòng)量方程、能量方程和組分方程，見(jiàn)式（1）～式（4）。黏性應(yīng)力張量表示見(jiàn)式（5）。

式中，ρ為密度；v為速度；p為壓力；hs為比焓；Yi為組分質(zhì)量分?jǐn)?shù)；τ為黏性應(yīng)力張量；Sρ、Sv、Sh和SY分別為密度、速度、焓和組分源項(xiàng)；μe和αe分別為有效動(dòng)力黏度和有效熱擴(kuò)散系數(shù)；(?v)T為速度梯度的轉(zhuǎn)置；I為單位矩陣。

液相與氣相間的傳質(zhì)及其他作用通過(guò)在控制方程中添加作用項(xiàng)實(shí)現(xiàn)，因噴霧過(guò)程中不考慮燃燒反應(yīng)，所以控制方程中未求解由于化學(xué)反應(yīng)產(chǎn)生的源項(xiàng)。

在大渦模擬中，上述控制方程經(jīng)過(guò)空間濾波將會(huì)產(chǎn)生亞網(wǎng)格尺度（subgrid scale，SGS）項(xiàng)，例如亞網(wǎng)格應(yīng)力張量τSGS、亞網(wǎng)格熱流量hSGS等無(wú)法直接求解的量，因此需要亞網(wǎng)格模型來(lái)封閉控制方程。本文選取的亞網(wǎng)格模型為動(dòng)態(tài)單方程模型［8］，其結(jié)構(gòu)簡(jiǎn)單且有足夠的精度，得到廣泛的應(yīng)用。

1.2 狀態(tài)方程及熱動(dòng)力學(xué)屬性

超臨界流體的部分參數(shù)會(huì)偏離亞臨界環(huán)境下的變化規(guī)律，為了準(zhǔn)確描述流體的物性參數(shù)變化，需要采用真實(shí)流體的狀態(tài)方程及熱物性和輸運(yùn)系數(shù)的計(jì)算方法。由于立方型狀態(tài)方程具有形式簡(jiǎn)單、精度高、計(jì)算成本低等優(yōu)點(diǎn)被大量使用，常用立方型狀態(tài)方程有Soave-Redlich-Kwong（SRK）和Peng-Robinson（PR）。文獻(xiàn)［9］中對(duì)比了SRK 和PR 兩種狀態(tài)方程對(duì)超臨界環(huán)境下噴霧預(yù)測(cè)的差異性，根據(jù)其研究結(jié)果，本文中選用的是PR 狀態(tài)方程，具體形式見(jiàn)式（6）～式（10）。

式中，T為流體的溫度；Mw為流體的物質(zhì)的量；pc和Tc分別為臨界壓力和臨界溫度；ω為流體的偏心因子；Ru為氣體常數(shù)；Spr、a和b為物質(zhì)性質(zhì)決定的常數(shù)。

超臨界流體黏度系數(shù)和導(dǎo)熱系數(shù)的計(jì)算選用文獻(xiàn)［10］中提出的計(jì)算方法，基本思想為利用特定的混合準(zhǔn)則計(jì)算多組分混合物在低壓下的黏度系數(shù)和導(dǎo)熱系數(shù)，之后利用修正函數(shù)得到高溫高壓狀態(tài)下的黏度系數(shù)和導(dǎo)熱系數(shù)。以黏度系數(shù)為例給出了計(jì)算過(guò)程。

式中，c、d為模型常數(shù)；η0為低壓黏度系數(shù)；η為高壓黏度系數(shù)；ηp為模型中間變量，無(wú)實(shí)際物理意義；Mwm、Vcm和Tcm分別為混合物的物質(zhì)的量、臨界體積和臨界溫度；k和?分別為勢(shì)能和玻爾茲曼常數(shù)；A6～A10為偏心因子和偶極矩的線性函數(shù)；Ω*、Fcm、T*和G2的具體計(jì)算方法見(jiàn)文獻(xiàn)［10］。

質(zhì)量擴(kuò)散系數(shù)采用文獻(xiàn)［11］中提出的高溫高壓修正方法，其利用Tr構(gòu)造了超臨界狀態(tài)下的修正函數(shù)，具體形式見(jiàn)式（16）。

式中，DH和DL分別為高壓和低壓下的質(zhì)量擴(kuò)散系數(shù)；Tr為溫度與臨界溫度的比值；A、B、C和E是隨壓力與臨界壓力的比值pr變化的模型經(jīng)驗(yàn)系數(shù)。

上述超臨界環(huán)境下物性參數(shù)計(jì)算模型的準(zhǔn)確性已從模型自身［10-11］和其在超臨界環(huán)境下的噴霧模擬中得到試驗(yàn)數(shù)據(jù)的驗(yàn)證，大量關(guān)于超臨界狀態(tài)下的模擬研究［1，12］均使用了此計(jì)算模型。

2 數(shù)值模擬

針對(duì)超臨界噴霧現(xiàn)象，本文中選取的模擬對(duì)象為美國(guó)Sandia 國(guó)家實(shí)驗(yàn)室關(guān)于正十二烷的噴射試驗(yàn)數(shù)據(jù)，該試驗(yàn)被大量用來(lái)進(jìn)行超臨界噴霧的研究。文獻(xiàn)［13］中的研究證明了此試驗(yàn)中存在超臨界狀態(tài)的油氣混合過(guò)程。根據(jù)已有試驗(yàn)數(shù)據(jù)建立噴霧模型，該試驗(yàn)將正十二烷通過(guò)直徑為0.09 mm 的噴孔噴入定容彈內(nèi)，由于不考慮燃燒過(guò)程，定容彈內(nèi)的環(huán)境氣體中不含O2，由N2、CO2和H2O 組成，密度約為22.8 kg/m3。初始溫度為363 K 的正十二烷在150 MPa 的壓力下噴入溫度900 K、壓力6.0 MPa 的環(huán)境中，而正十二烷的臨界溫度和臨界壓力分別為658 K 和1.82 MPa，在噴霧過(guò)程中會(huì)跨越燃油臨界點(diǎn)，使得燃油在超臨界環(huán)境中進(jìn)行霧化和混合。表1 給出了模擬中具體的邊界條件，其中N2、CO2和H2O 的物質(zhì)的量分?jǐn)?shù)分別為89.71 %、6.52 %和3.77 %。

表1 邊界條件

為探究超臨界和亞臨界環(huán)境下燃油擴(kuò)散混合的異同點(diǎn)，在上述模擬條件的基礎(chǔ)上，更改環(huán)境氣體的溫度和壓力，使整個(gè)噴霧過(guò)程都處于亞臨界環(huán)境中。本文中根據(jù)文獻(xiàn)［14］中的研究，環(huán)境溫度和壓力設(shè)為440 K 和2.9 MPa，并保證了環(huán)境氣體密度不變，均為22.8 kg/m3，表2 給出了兩種環(huán)境下的算例設(shè)置。

表2 不同環(huán)境下算例設(shè)置

本研究所采用的計(jì)算流體力學(xué)（computational fluid dynamics，CFD）軟件平臺(tái)為OpenFoam，該軟件平臺(tái)具有開(kāi)源、免費(fèi)等優(yōu)點(diǎn)，其自帶求解器類型齊全，泛用性強(qiáng)并可基于用戶需求自行擴(kuò)展，得到廣泛的應(yīng)用。本文中使用的SprayFoam 求解器基于歐拉-拉格朗日噴霧模型，求解器耦合PR 狀態(tài)方程和真實(shí)流體輸運(yùn)參數(shù)的計(jì)算方法，并結(jié)合大渦模擬策略，構(gòu)建了適用于超臨界環(huán)境下的噴霧模型。

計(jì)算域?yàn)橹睆?6 mm、高50 mm 的圓柱型區(qū)域，采用O-Block 對(duì)圓柱型計(jì)算域進(jìn)行劃分，并對(duì)局部網(wǎng)格進(jìn)行加密，總網(wǎng)格數(shù)量約為330 萬(wàn)，最小網(wǎng)格尺寸為50 μm。文獻(xiàn)［15］中研究了250.00 μm、125.00 μm、62.50 μm 和41.67 μm 網(wǎng)格對(duì)正十二烷噴射試驗(yàn)?zāi)M結(jié)果的影響，其指出大渦模擬采用62.50 μm 網(wǎng)格便能準(zhǔn)確捕捉噴霧特征，由此可知本文中使用的網(wǎng)格可以保證計(jì)算結(jié)果的準(zhǔn)確性。時(shí)間步長(zhǎng)由最大庫(kù)朗數(shù)控制，最大庫(kù)朗數(shù)為0.33，時(shí)間步長(zhǎng)約為3×10-8s。

3 結(jié)果分析

圖1 為超臨界和亞臨界環(huán)境下試驗(yàn)和計(jì)算的貫穿距隨噴射時(shí)間變化的對(duì)比圖。其中氣相貫穿距定義為燃油質(zhì)量分?jǐn)?shù)為0.1% 時(shí)對(duì)應(yīng)的位置與燃油噴嘴之間的最遠(yuǎn)距離，液相貫穿距定義為軸線方向上燃油質(zhì)量為已噴燃油90% 區(qū)域所占的軸向距離。從圖1 中可以看出計(jì)算模型能準(zhǔn)確捕捉不同環(huán)境下噴霧的液相和氣相貫穿距。在超臨界環(huán)境下的噴霧初期計(jì)算的液相貫穿距稍高于計(jì)算值，這也導(dǎo)致了氣相貫穿距的增加，亞臨界環(huán)境下也有類似的情況，但整體上計(jì)算結(jié)果和試驗(yàn)值基本吻合。本文中構(gòu)建的噴霧求解器可以在寬范圍的溫度和壓力下準(zhǔn)確預(yù)測(cè)噴霧特征。

圖1 不同環(huán)境下試驗(yàn)和計(jì)算噴霧貫穿距對(duì)比

為更好地證明超臨界環(huán)境下計(jì)算結(jié)果的準(zhǔn)確性，圖2 給出了超臨界環(huán)境下距離噴油位置處軸向距離X=17.85 mm 和25.00 mm 處沿徑向燃油質(zhì)量分?jǐn)?shù)分布的試驗(yàn)和計(jì)算結(jié)果對(duì)比圖。從試驗(yàn)和計(jì)算結(jié)果的對(duì)照中可以看出計(jì)算的燃油質(zhì)量分?jǐn)?shù)分布基本上與試驗(yàn)結(jié)果吻合，在噴霧邊緣的燃油質(zhì)量分?jǐn)?shù)稍微高出了試驗(yàn)值，這是由噴霧形態(tài)分布的不確定性所致。

圖2 不同軸向距離處沿徑向燃油質(zhì)量分?jǐn)?shù)分布

3.1 射流擴(kuò)散率

超臨界環(huán)境下燃油液滴的蒸發(fā)、熱物性參數(shù)和油氣混合層的發(fā)展演化與亞臨界環(huán)境下相比都表現(xiàn)出明顯的特殊性，這也使超臨界環(huán)境下燃油的射流擴(kuò)散過(guò)程不同于亞臨界。下面從燃油的分布入手研究?jī)煞N環(huán)境下的射流擴(kuò)散的差異。

圖3 和圖4 分別給出了0.6 ms 時(shí)刻兩種環(huán)境下燃油蒸氣的密度分布圖，后面的分析都基于此時(shí)刻進(jìn)行。從圖中可以明顯地看出亞臨界環(huán)境下噴霧有更長(zhǎng)的貫穿距，而超臨界環(huán)境下噴霧沿徑向充分發(fā)展。從燃油蒸氣密度的大小可以看出，超臨界環(huán)境下的燃油充分蒸發(fā)，霧化效果更好。

圖3 超臨界環(huán)境下密度分布

為了更細(xì)致地描述圖3 和圖4 中的密度分布現(xiàn)象，本文中利用射流擴(kuò)散率來(lái)量化燃油和環(huán)境氣體混合的程度。射流擴(kuò)散率是一個(gè)能衡量射流擴(kuò)散混合效率的指標(biāo)。射流擴(kuò)散率可以利用射流的速度場(chǎng)和密度場(chǎng)計(jì)算，需要根據(jù)應(yīng)用場(chǎng)景進(jìn)行選擇，本文中基于密度場(chǎng)的半寬半峰（half-width，half-maximum）規(guī)則［16］來(lái)計(jì)算射流擴(kuò)散率和擴(kuò)散角度。

圖4 亞臨界環(huán)境下密度分布

圖5 展示了基于射流密度場(chǎng)的半寬半峰規(guī)則計(jì)算射流擴(kuò)散率的具體過(guò)程，其中Δρ為密度增加量，Lρ為半寬半峰長(zhǎng)度。擴(kuò)散率定義為L(zhǎng)ρ沿射流軸線的線性回歸斜率。為了保證計(jì)算的準(zhǔn)確性，線性回歸必須在Lρ線性增加的區(qū)域進(jìn)行，在射流中通常是自相擬區(qū)域。根據(jù)半寬半峰規(guī)則得到射流擴(kuò)散角度，如式（17）所示。

圖5 基于密度場(chǎng)的半寬半峰規(guī)則

基于上述射流擴(kuò)散率的概念，表3 給出了兩種環(huán)境下的射流擴(kuò)散率和擴(kuò)散角度。超臨界環(huán)境下的擴(kuò)散率和擴(kuò)散角度均大于亞臨界，其燃油擴(kuò)散混合效果也優(yōu)于亞臨界，利用此方法得的超臨界環(huán)境下擴(kuò)散角度與文獻(xiàn)［17］中的試驗(yàn)結(jié)果一致性較好。

表3 擴(kuò)散率和擴(kuò)散角度

3.2 渦的分布

為了更加深入地研究?jī)煞N環(huán)境下擴(kuò)散率差異的原因，需要從更微觀的層面進(jìn)行分析，從根本上探究不同環(huán)境下射流擴(kuò)散混合過(guò)程的影響因素。在高速射流噴霧過(guò)程中，射流與周圍環(huán)境流體間的速度梯度非常大，初始噴油區(qū)域附近會(huì)產(chǎn)生剪切層，剪切層受Kelvin-Helmholtz 不穩(wěn)定性的影響，會(huì)發(fā)生渦的卷起、配對(duì)和分裂等級(jí)聯(lián)過(guò)程。隨著渦的不斷產(chǎn)生并向下游運(yùn)動(dòng)，周圍氣體會(huì)不斷地被卷吸到射流中，從而使燃油與環(huán)境氣體混合。由此可見(jiàn)渦的分布對(duì)噴霧過(guò)程中油氣混合有著關(guān)鍵的影響，根據(jù)渦的分布可以分析兩種環(huán)境下油氣混合狀態(tài)的差異。

圖6 和圖7 給出了兩種環(huán)境下噴油位置附近區(qū)域內(nèi)渦量的分布圖。超臨界和亞臨界環(huán)境下渦的分裂位置分別約為距離噴油位置軸向距離X=6.50 mm 和10.00 mm 處。可以看出超臨界環(huán)境下渦開(kāi)始產(chǎn)生的位置更靠近射流上游，超臨界環(huán)境下渦更早產(chǎn)生，加速了射流的破碎，使得更多的周圍氣體被卷吸，這種作用造成噴霧初始階段超臨界環(huán)境下渦量分布表面多褶皺，而亞臨界環(huán)境下表面較光滑，超臨界環(huán)境下射流的破碎趨勢(shì)更加明顯。超臨界下渦更早地產(chǎn)生促進(jìn)了油氣混合，同時(shí)也增強(qiáng)噴霧沿徑向的發(fā)展，擴(kuò)展率增大，相應(yīng)造成貫穿距降低，此結(jié)論與圖3 和圖4 的密度分布相一致。

圖6 超臨界環(huán)境下渦量分布

圖7 亞臨界環(huán)境下渦量分布

超臨界環(huán)境下渦更早產(chǎn)生的原因有待探究，需要從更深層次的渦的演化角度分析。

3.3 渦的演化

渦的分布直接影響著射流的破碎及其后續(xù)混合過(guò)程的進(jìn)行，為了進(jìn)一步分析渦的產(chǎn)生與發(fā)展，需要對(duì)渦的演化及其影響因素進(jìn)行探究，并量化各種影響因素所占的比例，從根本上分析超臨界和亞臨界環(huán)境下渦的分布不同的原因。根據(jù)文獻(xiàn)［18］中的研究，在不考慮其他外力的情況下，渦的演化方程為：

式中，d/dt代表物質(zhì)導(dǎo)數(shù)；右側(cè)第1 項(xiàng)表示渦的伸展，代表了速度梯度對(duì)渦的影響，它是產(chǎn)生復(fù)雜渦結(jié)構(gòu)的關(guān)鍵原因，并影響著湍流中動(dòng)能級(jí)聯(lián)過(guò)程；右側(cè)第2 項(xiàng)為渦的膨脹，這是由于流體體積的收縮和膨脹造成的，代表著流體的可壓縮性對(duì)渦的影響；右側(cè)第3 項(xiàng)為壓力梯度與密度梯度錯(cuò)位所產(chǎn)生的斜壓轉(zhuǎn)矩，是Rayleigh-Taylor 不穩(wěn)定性形成的原因；最后一項(xiàng)來(lái)自于由黏性應(yīng)力引起的渦的擴(kuò)散。

圖8 和圖9 分別為噴油位置附近區(qū)域內(nèi)兩種環(huán)境下不同項(xiàng)對(duì)渦的演化影響分布圖，自上而下依次為渦的伸展項(xiàng)、膨脹項(xiàng)、斜壓轉(zhuǎn)矩項(xiàng)和黏性應(yīng)力項(xiàng)。為便于分析，圖中僅給出了不同項(xiàng)垂直于平面方向的大小，在其他方向上各項(xiàng)具有同樣的趨勢(shì)。根據(jù)圖中各項(xiàng)的大小可以看出在兩種環(huán)境下渦的伸展的作用近乎相同，渦的膨脹作用也有同樣的趨勢(shì)，并且渦的伸展對(duì)渦的演化影響最大（達(dá)到1012水平）。在亞臨界環(huán)境下的斜壓轉(zhuǎn)矩和黏性應(yīng)力的影響較小，而斜壓轉(zhuǎn)矩的影響稍高于黏性應(yīng)力。在超臨界環(huán)境下斜壓轉(zhuǎn)矩的影響增強(qiáng)，與渦的膨脹作用相當(dāng)（均為1011），此時(shí)黏性應(yīng)力造成的影響遠(yuǎn)低于斜壓轉(zhuǎn)矩。

圖8 超臨界環(huán)境中不同項(xiàng)的影響

圖9 亞臨界環(huán)境中不同項(xiàng)的影響

整體對(duì)照兩種環(huán)境下不同項(xiàng)在渦的演化過(guò)程中的作用，可以明顯地看出斜壓轉(zhuǎn)矩的作用在超臨界環(huán)境中會(huì)大幅增強(qiáng)，其他3 項(xiàng)的影響基本保持不變，且3 項(xiàng)中速度梯度的作用最強(qiáng)。在高速射流中（例如本研究中的射流速度最高可達(dá)600 m/s），射流表面和周圍氣體之間會(huì)形成巨大的速度梯度，成為渦的演化過(guò)程中的決定因素，而渦的膨脹和黏性應(yīng)力的影響在溫度和壓力劇烈變化后仍保持不變，說(shuō)明周圍環(huán)境氣體的影響在高速射流中會(huì)相對(duì)降低，該現(xiàn)象是由于射流高速運(yùn)動(dòng)使得其與環(huán)境的相互作用時(shí)間縮短所致。

超臨界環(huán)境中斜壓轉(zhuǎn)矩的作用增強(qiáng)與其形成較大的密度梯度有關(guān)，圖10 和圖11 給出了兩種環(huán)境下的密度梯度分布。超臨界環(huán)境下的密度梯度明顯高于亞臨界，超臨界狀態(tài)下的大密度梯度使密度梯度和壓力梯度的錯(cuò)位所產(chǎn)生的斜壓轉(zhuǎn)矩增大，強(qiáng)化了渦的發(fā)展，進(jìn)而對(duì)混合層的演化過(guò)程造成重要的影響。由此可見(jiàn)，超臨界環(huán)境下產(chǎn)生的較大密度梯度是其提前出現(xiàn)渦的分離的關(guān)鍵影響因素，說(shuō)明密度梯度對(duì)射流混合過(guò)程起著重要作用，密度梯度的大小直接影響著超臨界環(huán)境噴霧的發(fā)展。

圖10 超臨界環(huán)境下密度梯度分布

圖11 亞臨界環(huán)境下密度梯度分布

由此可見(jiàn)超臨界條件下渦更早地產(chǎn)生的原因是超臨界下形成的較大密度梯度層，密度梯度的增大加大了壓力梯度與密度梯度錯(cuò)位所產(chǎn)生的斜壓轉(zhuǎn)矩，即渦的演化公式中?(1p) ×?p項(xiàng)增大，加速了渦的產(chǎn)生。超臨界環(huán)境下射流密度梯度的增大是由于射流進(jìn)入超臨界環(huán)境時(shí)其與周圍環(huán)境的質(zhì)量擴(kuò)散率接近于零［19］，射流與環(huán)境氣體的物質(zhì)交換受到抑制，使其形成一個(gè)高密度的混合層，混合層與環(huán)境氣體間產(chǎn)生較大的密度梯度。

綜上，欲強(qiáng)化渦的產(chǎn)生與演化以增強(qiáng)卷吸作用來(lái)強(qiáng)化燃料的噴霧混合：首先需要在噴油裝置允許的壓力條件下盡可能提高噴油壓力，提高射流與環(huán)境氣體間的速度梯度，加速渦的伸展；再者要提高環(huán)境溫度和壓力，使射流進(jìn)入超臨界狀態(tài)，以提高射流與環(huán)境間壓力梯度和密度梯度錯(cuò)位所產(chǎn)生的斜壓轉(zhuǎn)矩，增大Rayleigh-Taylor 不穩(wěn)定性。通過(guò)上述噴射策略，可大幅提高燃油霧化效率，使之與環(huán)境氣體充分混合，利于后續(xù)燃燒過(guò)程的進(jìn)行。

4 結(jié)論

（1）基于射流密度分布，超臨界環(huán)境下的射流擴(kuò)散率和擴(kuò)散角度均大于亞臨界。

（2）超臨界環(huán)境會(huì)促進(jìn)渦的形成，使周圍更多的氣體被卷吸入射流中，加速射流的破碎及后續(xù)的油氣混合過(guò)程。

（3）在高速射流中，速度梯度對(duì)渦的演化過(guò)程起主要作用，超臨界環(huán)境下由于密度梯度的增大會(huì)使斜壓轉(zhuǎn)矩的作用增強(qiáng)，加速渦的演化。