閻 凱，寧 智，趙 晉，呂 明，孫春華，3

(1.北京交通大學 機械與電子控制工程學院，北京 100044; 2.中國科學院電工研究所，北京 100190;3.上海交通大學 內(nèi)燃機研究所，上海 200240)

1 前 言

燃油霧化好壞關(guān)系內(nèi)燃機能否實現(xiàn)高效、清潔燃燒，而液體射流穩(wěn)定性對射流的破碎霧化具有決定性作用，其研究對于內(nèi)燃機燃燒技術(shù)的發(fā)展具有重要的理論意義和實際應用價值。

壓力旋流噴嘴具有可靠性高、霧化效果好、所需泵端壓力低等優(yōu)點[1]，已廣泛應用于航空和船用發(fā)動機、燃氣輪機及工業(yè)窯爐等噴射系統(tǒng)[2-3]。隨著技術(shù)的不斷發(fā)展，壓力旋流噴嘴也逐漸用在內(nèi)燃機，特別是汽油缸內(nèi)直噴(GDI)發(fā)動機的燃油噴射中[4]。

在內(nèi)燃機中，壓力旋轉(zhuǎn)噴嘴形成的圓環(huán)旋轉(zhuǎn)黏性液體射流與環(huán)境存在溫度差異。射流液體與周圍環(huán)境氣體間存在溫差時，溫度擾動的存在能導致液體射流界面穩(wěn)定性發(fā)生變化。雖然作者已研究了黏性、表面張力、旋流強度等對圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性的影響[5-7]，但并未考慮溫度擾動的影響。

近年來，國內(nèi)外一些學者采用射流線性穩(wěn)定性分析方法，對存在溫度擾動的液體射流穩(wěn)定性問題進行了研究，分析非均勻液體溫度、均勻環(huán)境溫度、周期性變化環(huán)境溫度等條件下的溫度擾動對液體射流穩(wěn)定性的影響，取得了重要研究成果[8-13]。研究表明，當射流表面存在溫度擾動時，射流表面將產(chǎn)生熱毛細力，進而產(chǎn)生Marangoni流動，從而對射流穩(wěn)定性產(chǎn)生影響；另外，射流液體的黏性、導熱率、比熱容等物性參數(shù)不同時，溫度擾動對射流穩(wěn)定性的影響會發(fā)生變化。目前，采用射流穩(wěn)定性方法對存在溫度擾動的液體射流穩(wěn)定性問題的研究，均是在比較簡單的條件下進行，如存在溫度擾動的二維平面液體射流或圓柱軸對稱液體射流穩(wěn)定性問題的研究等。存在溫度擾動的壓力旋流噴嘴或空氣助力旋流噴嘴等形成的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性問題的研究還未見報道。

本文采用射流穩(wěn)定性方法對存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性問題進行研究。建立同時存在速度擾動、壓力擾動以及溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性色散方程，并對色散方程進行空間模式下的數(shù)值求解。在此基礎上分析存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流的穩(wěn)定性，研究溫度擾動強度對圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性的作用，以及射流液體熱量傳遞特性對溫度擾動與圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性作用的影響。

2 色散方程的建立

本文討論的圓環(huán)旋轉(zhuǎn)黏性液體射流是具有兩個自由液面的中空環(huán)狀射流。圖1給出的是柱坐標系(r, θ, z)下圓環(huán)旋轉(zhuǎn)黏性液體射流示意圖。

圖1 圓環(huán)旋轉(zhuǎn)黏性液體射流示意圖Fig.1 Schematic diagram of annular swirl viscous liquid sheet

假定射流方向與z軸相反；氣液流體均為不可壓縮牛頓流體；圓環(huán)射流液體內(nèi)半徑為a，外半徑為b，厚度為h；液體兩側(cè)相位差為Φ；射流液體動力黏度為μ，射流周圍氣體為無黏流體；氣液界面處表面張力系數(shù)為σ；射流液體密度為ρl，圓環(huán)射流內(nèi)部氣體密度為ρi，外部氣體密度為ρo；射流液體初始溫度為Tl0，環(huán)形射流內(nèi)外部氣體具有相同溫度Tg，液體內(nèi)部具有雙線性溫度分布；液體射流旋轉(zhuǎn)速度分布采用固體渦核型旋轉(zhuǎn)速度分布Ar，A為液體旋轉(zhuǎn)角速度，r為半徑。

對射流控制方程組(N-S方程組)進行擾動，并對擾動后的射流控制方程組進行線性化和無量綱化處理，經(jīng)過推導得到無量綱形式的線性化射流擾動控制方程：

依據(jù)簡正模態(tài)法，設無量綱擾動速度、擾動壓力及擾動溫度有如下形式的解[14]：

將無量綱擾動速度、擾動壓力及擾動溫度代入射流擾動控制方程中，求解可得擾動量：

按照射流液體控制方程的求解過程，同樣可以求解得到氣體的擾動速度、擾動壓力及擾動溫度。

對射流擾動控制方程的運動學、動力學及熱力學邊界條件進行擾動、線性化及無量綱化處理，得到無量綱線性化擾動控制方程的邊界條件。

無量綱Marangoni數(shù)表征表面張力隨溫度的變化率，具有式(14)的形式。

圓環(huán)液膜內(nèi)自由液面動力學邊界條件(r=Ah)：

其中，WeT=We/(1-We?Ma)。

圓環(huán)液膜外自由液面動力學邊界條件(r=Bh)：

圓環(huán)旋轉(zhuǎn)黏性液體射流無量綱線性化擾動控制方程的熱力學邊界條件包括溫度連續(xù)性條件和熱流密度連續(xù)性條件。

無量綱化的溫度連續(xù)性條件：

無量綱化的熱流密度連續(xù)性條件：

通過整理得到12個無量綱射流邊界條件。將擾動速度、擾動壓力、擾動溫度代入12個無量綱射流邊界條件中，經(jīng)整理得到如下齊次方程組：

其中，X=[c11，c12，e11，e12，d11，d12，d21，d22，cT1，cT2，cg1，cg2]T；Λ(ω，k，m)是一個 12 階方陣。

式(23)為齊次線性代數(shù)方程組，有非零解的條件是方程組系數(shù)行列式為零；由此可得：

式(24)即為建立的描述存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流自由表面擾動發(fā)展的色散方程。

鑒于色散方程的嚴重非線性及復雜性，只能對其進行數(shù)值求解。以往研究者對該類問題的數(shù)學模型進行迭代求解時，大都采用拋物線法(亦稱 Müller法)。然而，與弦截法相比，雖然拋物線法的收斂速度略快，但其算法比較復雜，既要處理符號問題，還要考慮復數(shù)的運算。為此，本文采用弦截法對色散方程進行空間模式數(shù)值求解。色散方程是一個12階行列式，數(shù)值求解結(jié)果對初值十分敏感；初值選定不合適的話，有可能導致計算結(jié)果的不收斂。通過大量的數(shù)值計算發(fā)現(xiàn)，對于一定的擾動波數(shù)，求解擾動頻率比直接求解擾動增長率更容易收斂。因此，本文在數(shù)值求解擾動增長率時，首先確定一定擾動波數(shù)所對應的擾動頻率，然后再求解擾動增長率。

圖2給出的是Ma、Β、Λi和Λo、Pei和Peo為0時，射流擾動增長率的計算結(jié)果與文獻[14]給出的計算結(jié)果的比較。從圖2中可以看到，兩個計算結(jié)果相同。這從一定角度說明存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性色散方程以及色散方程數(shù)值求解方法的正確性。

圖2 射流擾動增長率的計算結(jié)果與文獻[15]中的結(jié)果對比Fig.2 Comparison of calculation results of disturbance growth rate with the results in the literature[15]

3 存在溫度擾動圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性

簡單條件下的溫度擾動對液體射流穩(wěn)定性影響的研究，如存在溫度擾動的二維平面液體射流或圓柱軸對稱液體射流穩(wěn)定性問題等，已經(jīng)證實溫度擾動的存在會導致液體射流界面的穩(wěn)定性發(fā)生變化；對存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性問題進行研究，研究過程中選擇的工質(zhì)為水與乳化劑的混合溶液，其動力黏度為 0.006 kg?m-1·s-1，初始入口速度為 6 m?s-1。

由于建立的描述存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性的色散方程中假設壓力擾動、速度擾動及溫度擾動皆具有相同的擾動頻率，因此下述分析中提到的擾動頻率是指壓力擾動、速度擾動以及溫度擾動的共同擾動頻率。

目前，射流線性穩(wěn)定性分析方法還無法直接對擾動強度與射流穩(wěn)定性之間的關(guān)系進行研究。但有研究認為，射流液體與周圍氣體間的溫差越大，溫度擾動強度越大[16]。因此，本文通過氣液溫差對射流穩(wěn)定性的影響來間接反映溫度擾動強度對射流穩(wěn)定性的作用，主要討論射流周圍氣體溫度高于射流液體溫度的情況。

最大擾動增長率反映的是射流的不穩(wěn)定程度。將最大擾動增長率對應的波數(shù)稱為最不穩(wěn)定波數(shù)，它與射流破碎液滴粒徑具有一定的對應關(guān)系[1]：

其中，kL為最不穩(wěn)定波數(shù)；

Ct為常數(shù)；f為與雷諾數(shù)和韋伯數(shù)有關(guān)的待定系數(shù)。

從式(25)中可以看到，最不穩(wěn)定波數(shù)越大，射流破碎液滴粒徑越小。

將射流最大擾動增長率對應的頻率稱為最不穩(wěn)定頻率，射流最不穩(wěn)定頻率與最不穩(wěn)定波數(shù)具有線性對應的關(guān)系。

不穩(wěn)定截斷頻率指的是不穩(wěn)定頻率的最大值，反映的是射流不穩(wěn)定頻率范圍和射流破碎液滴粒徑范圍。不穩(wěn)定截斷頻率越大，射流不穩(wěn)定頻率范圍越大、射流破碎產(chǎn)生的最小液滴粒徑越小、射流破碎液滴粒徑范圍越大。

通過最大擾動增長率、最不穩(wěn)定頻率以及不穩(wěn)定截斷頻率對存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性進行研究。為了表達上的簡潔，將軸對稱擾動下的射流穩(wěn)定性和非軸對稱擾動下的射流穩(wěn)定性簡稱為軸對稱穩(wěn)定性和非軸對稱穩(wěn)定性。

3.1 射流的軸對稱穩(wěn)定性與非軸對稱穩(wěn)定性

射流是否旋轉(zhuǎn)是射流軸對稱穩(wěn)定性和非軸對稱穩(wěn)定性的重要影響因素。

圖3和圖4分別是射流不旋轉(zhuǎn)和旋轉(zhuǎn)條件下，存在溫度擾動與不存在溫度擾動時圓環(huán)黏性液體射流的擾動增長率隨擾動頻率變化的比較。

圖3 射流不旋轉(zhuǎn)時溫度擾動對射流擾動增長率隨擾動頻率變化的影響Fig.3 Effects of temperature disturbance on disturbance growth rate without liquid swirling

圖4 射流旋轉(zhuǎn)時溫度擾動對射流擾動增長率隨擾動頻率變化的影響Fig.4 Effects of temperature disturbance on disturbance growth rate with liquid swirling

從圖3可以看到，射流不旋轉(zhuǎn)時，無論是軸對稱擾動還是非軸對稱擾動，存在溫度擾動時的擾動增長率隨擾動頻率的變化規(guī)律與不存在溫度擾動時的變化規(guī)律基本相同。然而，存在溫度擾動時，射流的最大擾動增長率、最不穩(wěn)定頻率以及不穩(wěn)定截斷頻率明顯大于不存在溫度擾動時的相應參數(shù)。說明溫度擾動的存在有利于圓環(huán)黏性液體射流的失穩(wěn)及破碎，也有利于射流破碎液滴粒徑的減小以及破碎液滴粒徑范圍的增大。

對于無旋轉(zhuǎn)圓環(huán)黏性液體射流來說，不論是否存在溫度擾動，軸對稱擾動射流的最大擾動增長率皆大于非軸對稱擾動射流。說明無旋轉(zhuǎn)圓環(huán)黏性液體射流的占優(yōu)模式為軸對稱擾動模式，溫度擾動不改變無旋轉(zhuǎn)圓環(huán)黏性液體射流的占優(yōu)模式。

從圖4中可以看到，與射流不旋轉(zhuǎn)時相比，在射流旋轉(zhuǎn)條件下，溫度擾動的存在對圓環(huán)黏性液體射流穩(wěn)定性產(chǎn)生更大的影響。

從圖4(a)中可以發(fā)現(xiàn)，在射流旋轉(zhuǎn)條件下，不存在溫度擾動時，不穩(wěn)定擾動頻率的變化范圍從零開始。而存在溫度擾動時，不穩(wěn)定擾動頻率的變化范圍則從一定擾動頻率開始。說明射流旋轉(zhuǎn)時，溫度擾動使得頻率低于一定值的擾動不會對射流的軸對稱不穩(wěn)定性產(chǎn)生影響，溫度擾動的存在有利于提高一定擾動頻率下(主要是擾動頻率較低情況下)圓環(huán)旋轉(zhuǎn)黏性液體射流的軸對稱穩(wěn)定性。這有可能是在較低的擾動頻率下溫度擾動與速度擾動及壓力擾動產(chǎn)生相互抵消的結(jié)果。另外，從圖4(a)中還可以看到，在射流旋轉(zhuǎn)條件下，存在溫度擾動時，不僅射流的最大擾動增長率、最不穩(wěn)定頻率以及不穩(wěn)定截斷頻率會顯著增大，而且擾動增長率隨擾動頻率的變化規(guī)律也會發(fā)生一定的變化。從圖4(b)中可以看到，在射流旋轉(zhuǎn)條件下，溫度擾動對非軸對稱擾動也有類似的影響。溫度擾動的存在同樣有利于射流旋轉(zhuǎn)時圓環(huán)黏性液體射流的失穩(wěn)及破碎、射流破碎液滴粒徑的減小以及破碎液滴粒徑范圍的增大。

對比圖4(a)和(b)可見溫度擾動時，非軸對稱擾動的不穩(wěn)定頻率范圍略大于軸對稱擾動的不穩(wěn)定頻率范圍，而不穩(wěn)定起始頻率和不穩(wěn)定截斷頻率皆小于軸對稱擾動的不穩(wěn)定起始頻率和不穩(wěn)定截斷頻率。

另外，對比圖4(a)和(b)還可以發(fā)現(xiàn)，不存在溫度擾動時，軸對稱擾動射流的最大擾動增長率小于非軸對稱擾動射流，非軸對稱擾動模式占優(yōu)。而存在溫度擾動時，軸對稱擾動射流的最大擾動增長率卻明顯大于非軸對稱擾動射流，軸對稱擾動模式占優(yōu)。說明在射流旋轉(zhuǎn)條件下，溫度擾動使圓環(huán)旋轉(zhuǎn)黏性液體射流的占優(yōu)模式發(fā)生改變，即由非軸對稱擾動模式轉(zhuǎn)變?yōu)檩S對稱擾動模式。這可能是由于溫度擾動與角向模數(shù)(反映擾動的類型)相互作用的結(jié)果。

3.2 射流溫差對射流穩(wěn)定性的作用

有研究認為，射流液體與周圍氣體間溫差在一定程度上反映了溫度擾動強度[13]。以無量綱參數(shù)B=βh/Tl0表征氣液溫差，Β越大，氣液溫差越大。

圖5和6分別給出的是射流不旋轉(zhuǎn)和旋轉(zhuǎn)時，圓環(huán)黏性液體射流的最大擾動增長率、最不穩(wěn)定頻率及不穩(wěn)定截斷頻率隨Β的變化。

圖5 射流不旋轉(zhuǎn)時擾動參數(shù)隨無量綱數(shù)Β的變化Fig.5 Variation of the disturbance frequency characteristics as a function of B without liquid swirling

從圖5可見，射流不旋轉(zhuǎn)時，Β對軸對稱擾動和非軸對稱擾動射流的最大擾動增長率具有基本相同的作用，即Β增加，最大擾動增長率增大；Β對軸對稱擾動的影響略大；不同Β時，軸對稱擾動射流的最大擾動增長率總是大于非軸對稱擾動射流。說明射流不旋轉(zhuǎn)時，氣液溫差，即溫度擾動強度增大有利于射流的失穩(wěn)及破碎。這與一般的理解基本一致，氣液溫差不改變射流的占優(yōu)模式。

射流不旋轉(zhuǎn)時，Β對軸對稱擾動射流和非軸對稱擾動射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率幾乎沒有影響。說明射流不旋轉(zhuǎn)時，氣液溫差，亦即溫度擾動強度對射流破碎液滴粒徑以及射流不穩(wěn)定頻率范圍和破碎液滴粒徑范圍沒有影響。

射流不旋轉(zhuǎn)時，軸對稱擾動射流的最不穩(wěn)定頻率和不穩(wěn)定截斷頻率總是大于非軸對稱擾動射流。說明不同氣液溫差，亦即不同溫度擾動強度時，軸對稱擾動射流的破碎液滴粒徑總是小于非軸對稱擾動射流，軸對稱擾動射流的破碎液滴粒徑范圍總是大于非軸對稱擾動射流。

通過圖6和圖5的比較可以看到，射流旋轉(zhuǎn)時，除了最大擾動增長率明顯增大以外，軸對稱擾動和非軸對稱擾動射流的最大擾動增長率隨Β的變化規(guī)律幾乎與射流不旋轉(zhuǎn)時完全相同。

圖6 射流旋轉(zhuǎn)時擾動參數(shù)隨無量綱數(shù)Β的變化Fig.6 Variation of the disturbance frequency characteristics as a function of B with liquid swirling

與射流不旋轉(zhuǎn)時Β對射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率幾乎沒有影響不同，射流旋轉(zhuǎn)時，隨著Β的增加，射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率增大。說明射流旋轉(zhuǎn)時，氣液溫差，亦即溫度擾動強度的增加有利于射流破碎液滴粒徑的減小及射流破碎液滴粒徑范圍的增大。為此，為獲得更好的霧化效果，在工程應用時可考慮增加圓環(huán)射流周圍氣體的溫度。

射流旋轉(zhuǎn)時，雖然軸對稱擾動和非軸對稱擾動射流的最不穩(wěn)定頻率和不穩(wěn)定截斷頻率皆隨Β的增加而增大，但軸對稱擾動射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率受Β的影響更大。

射流旋轉(zhuǎn)時，不同Β下軸對稱擾動射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率總是大于非軸對稱擾動射流，并且兩者間的差值隨著Β的增加而增大，說明隨著氣液溫差，亦即溫度擾動強度增加，軸對稱擾動射流破碎液滴粒徑小于非軸對稱擾動射流破碎液滴粒徑的程度以及破碎液滴粒徑范圍大于非軸對稱擾動射流破碎液滴粒徑范圍的程度增大。

3.3 貝克來數(shù)對射流穩(wěn)定性的作用

貝克來數(shù)Pe反映的是射流液體的熱量傳遞特性與溫度擾動對圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性作用的影響。

圖7給出的是射流不旋轉(zhuǎn)和旋轉(zhuǎn)時，圓環(huán)黏性液體射流最大擾動增長率隨貝克來數(shù)Pe的變化。

從圖7中可以看到，隨著Pe增加，不論是軸對稱擾動還是非軸對稱擾動，射流最大擾動增長率皆呈增大趨勢；說明射流液體的熱量傳遞特性越差，圓環(huán)黏性液體射流越容易失穩(wěn)。隨著Pe的增加，不旋轉(zhuǎn)射流最大擾動增長率的增大幅度略有減小，而旋轉(zhuǎn)射流最大擾動增長率的增大幅度略有增大。在不同Pe下，軸對稱擾動射流最大擾動增長率總是大于非軸對稱擾動射流，且隨著Pe的增加，軸對稱擾動與非軸對稱擾動射流最大擾動增長率之間的差值增大。說明射流液體的熱量傳遞特性對圓環(huán)黏性液體射流的占優(yōu)模式?jīng)]有影響，且射流液體熱量傳遞特性越差，軸對稱擾動模式占優(yōu)程度越高。

圖7 最大擾動增長率隨貝克來數(shù)Pe的變化Fig.7 Variation of the maximum disturbance growth rate as a function of Peclet number

圖8給出的是射流不旋轉(zhuǎn)和旋轉(zhuǎn)時，圓環(huán)黏性液體射流的最不穩(wěn)定頻率和不穩(wěn)定截斷頻率隨Pe的變化。從圖8可以看到，射流旋轉(zhuǎn)和不旋轉(zhuǎn)兩種情況下，Pe對圓環(huán)黏性液體射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率的影響表現(xiàn)出明顯不同特性。

圖8 最不穩(wěn)定頻率和不穩(wěn)定截斷頻率隨貝克來數(shù)Pe的變化Fig.8 Variation of the disturbance frequency characteristics as a function of Peclet number

需要說明的是，式(25)反映了最不穩(wěn)定波數(shù)與射流破碎液滴粒徑的關(guān)系，而射流最不穩(wěn)定頻率與最不穩(wěn)定波數(shù)具有線性對應的關(guān)系。所以本小節(jié)的敘述中將根據(jù)射流最不穩(wěn)定頻率的大小及范圍來反映射流破碎液滴粒徑及分布。從圖8(a)中可以看到，射流不旋轉(zhuǎn)時，隨著Pe的增加，軸對稱擾動射流最不穩(wěn)定頻率減小，且減小趨勢逐漸加快；非軸對稱擾動射流最不穩(wěn)定頻率減小，但減小趨勢逐漸變慢；軸對稱擾動射流最不穩(wěn)定頻率始終大于非軸對稱擾動射流。說明射流不旋轉(zhuǎn)時，射流液體的熱量傳遞特性越差，射流破碎液滴粒徑越大，而軸對稱擾動射流的破碎液滴粒徑總是小于非軸對稱擾動射流。

射流不旋轉(zhuǎn)時，隨著Pe的增加，軸對稱擾動和非軸對稱擾動射流的不穩(wěn)定截斷頻率皆基本保持不變，且軸對稱擾動射流不穩(wěn)定截斷頻率總是大于非軸對稱擾動射流。說明射流不旋轉(zhuǎn)時，射流液體的熱量傳遞特性對射流破碎液滴粒徑范圍沒有影響，軸對稱擾動射流的破碎液滴粒徑范圍總是大于非軸對稱擾動射流。

從圖8(b)中可以看到，射流旋轉(zhuǎn)時，最不穩(wěn)定頻率隨Pe的變化趨勢與射流不旋轉(zhuǎn)時完全不同。隨著Pe的增加，射流最不穩(wěn)定頻率增大，且軸對稱擾動射流最不穩(wěn)定頻率的增大有加快趨勢，而非軸對稱擾動射流最不穩(wěn)定頻率的增大有減慢趨勢；軸對稱擾動射流的最不穩(wěn)定頻率總是大于非軸對稱擾動射流。說明射流旋轉(zhuǎn)時，射流液體的熱量傳遞特性越差，射流破碎液滴粒徑越?。惠S對稱擾動射流的破碎液滴粒徑總是小于非軸對稱擾動射流，這一點與射流不旋轉(zhuǎn)時相同。

射流旋轉(zhuǎn)時，不穩(wěn)定截斷頻率隨Pe的變化與射流不旋轉(zhuǎn)時完全不同。隨著Pe的增加，不穩(wěn)定截斷頻率增大。軸對稱擾動射流不穩(wěn)定截斷頻率總是大于非軸對稱擾動射流。說明射流旋轉(zhuǎn)時，射流液體的熱量傳遞特性越差，破碎液滴粒徑范圍越大，而破碎液滴粒徑范圍總是大于非軸對稱擾動射流，這一點與射流不旋轉(zhuǎn)時相同。

4 結(jié) 論

(1) 推導出存在溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流擾動控制方程的解析解和邊界條件，建立了同時存在速度擾動、壓力擾動以及溫度擾動的圓環(huán)旋轉(zhuǎn)黏性液體射流穩(wěn)定性色散方程，對色散方程進行了空間模式求解，并進行了驗證分析。

(2) 溫度擾動的存在使得圓環(huán)旋轉(zhuǎn)黏性液體射流的最大擾動增長率、最不穩(wěn)定頻率和不穩(wěn)定截斷頻率增大，并使射流擾動增長率隨擾動頻率的變化規(guī)律改變。溫度擾動的存在有利于提高一定擾動頻率下射流的軸對稱擾動穩(wěn)定性，并使射流擾動的占優(yōu)模式發(fā)生變化。

(3) 對于旋轉(zhuǎn)射流，隨氣液溫差的增加，射流最大擾動增長率、最不穩(wěn)定頻率和不穩(wěn)定截斷頻率增大。軸對稱擾動射流最不穩(wěn)定頻率和不穩(wěn)定截斷頻率受氣液溫差影響更大。

(4) 在相同溫度擾動下，圓環(huán)旋轉(zhuǎn)射流的最大擾動增長率、最不穩(wěn)定頻率以及不穩(wěn)定截斷頻率皆隨Pe的增加而增大。軸對稱擾動射流的最大擾動增長率、最不穩(wěn)定頻率以及不穩(wěn)定截斷頻率總是大于非軸對稱擾動射流，且差值隨Pe的增加而增大。