沈彥成，吳偉亮

● （上海交通大學(xué) 動(dòng)力機(jī)械和工程教育部重點(diǎn)實(shí)驗(yàn)室，上海 200240）

壓力式噴嘴在壓力脈動(dòng)下性能的研究

沈彥成，吳偉亮

● （上海交通大學(xué) 動(dòng)力機(jī)械和工程教育部重點(diǎn)實(shí)驗(yàn)室，上海 200240）

以機(jī)械壓力噴嘴為研究對(duì)象，應(yīng)用 VOF模型對(duì)脈動(dòng)壓力下的噴嘴內(nèi)部流場(chǎng)進(jìn)行了數(shù)值模擬，分析了噴嘴在脈動(dòng)壓力下的流量、霧化角等特性。結(jié)果表明：脈動(dòng)壓力下瞬時(shí)流量和霧化角呈現(xiàn)相同頻率的脈動(dòng)，且平均流量基本在設(shè)計(jì)流量附近。流量的振幅隨脈動(dòng)頻率的變化，呈現(xiàn)出先增大后減小的趨勢(shì)，并在400Hz時(shí)達(dá)到最大值；入口壓力與出口流量的相位差隨頻率的增大而增加，且近乎線性關(guān)系，脈動(dòng)壓力的幅值對(duì)相位差影響不大。

機(jī)械壓力噴嘴；壓力脈動(dòng)；霧化角；流量

0 引言

狹義上的機(jī)械壓力式噴嘴是指能將液體霧化的裝置，這類裝置在工農(nóng)業(yè)生產(chǎn)中獲得大量應(yīng)用。機(jī)械壓力式噴嘴包含多種種類，其中結(jié)構(gòu)較簡(jiǎn)單的是單路壓力噴嘴，因其壓力損失小、霧化效果好，在多種設(shè)備，特別是燃燒設(shè)備中得到廣泛應(yīng)用。本文討論對(duì)象就是這種結(jié)構(gòu)類型的噴嘴。從應(yīng)用情況看，迄今僅作為使液體霧化的器件噴嘴，在設(shè)計(jì)、計(jì)算中僅考慮其在給定壓力下流量、霧化角、霧化粒徑及分布等基本靜態(tài)特性。實(shí)際燃燒裝置中，由于各種因素可能導(dǎo)致噴嘴上下游壓力會(huì)發(fā)生變化，如燃油泵葉片旋轉(zhuǎn)作用、不穩(wěn)定燃燒、航空器機(jī)動(dòng)等，從而導(dǎo)致噴嘴工作性能發(fā)生改變，影響動(dòng)力設(shè)備的運(yùn)行，因此研究噴嘴在脈動(dòng)壓力下的性能對(duì)現(xiàn)實(shí)具有一定參考價(jià)值。特別是人們熟悉和掌握了噴嘴在脈動(dòng) 壓力下的性能后，不但可以以此改進(jìn)燃燒設(shè)備工作性能，還有可能將噴嘴作為控制器件，對(duì)燃燒過程實(shí)施有目的的控制，并發(fā)展出新的燃燒方式。工作于動(dòng)態(tài)壓力下的噴嘴在工程中是很常見的現(xiàn)象，如在內(nèi)燃機(jī)工作過程中燃油噴嘴就是以動(dòng)態(tài)形式工作，對(duì)其霧化過程人們開展了大量研究工作。但這類噴嘴工作最大特點(diǎn)是脈沖式工作，其不在本文研究范疇之內(nèi)。本文所指噴嘴脈動(dòng)壓力下的性能是指噴嘴在鍋爐、燃?xì)廨啓C(jī)等動(dòng)力裝置中上下游壓力持續(xù)脈動(dòng)下的性能，即噴嘴工作壓力由足夠的直流量和一定幅值交流量疊加而成，如圖1所示

圖1 噴嘴工作壓力變化

俄羅斯學(xué)者Bazarov首先對(duì)噴嘴在脈壓力下的性能進(jìn)行了的研究，其在忽略流體粘性的前提下，通過理論推導(dǎo)得到機(jī)械壓力噴嘴頻率特性方程，并分析了噴嘴各輸入、輸出參數(shù)振幅和相位角的的關(guān)系[3]。文獻(xiàn)[4]在繼承 Bazarov思想的基礎(chǔ)上考慮了粘性，分析了幾何參數(shù)和工況參數(shù)對(duì)機(jī)械壓力噴嘴相頻特性和和幅頻特性的影響。同時(shí)，也有一些學(xué)者通過實(shí)驗(yàn)的方法對(duì)噴嘴在脈動(dòng)壓力下的性能進(jìn)行了研究。文獻(xiàn)[5-6]測(cè)量了機(jī)械壓力噴嘴在壓力脈動(dòng)下的瞬時(shí)流量。通過測(cè)量噴嘴出口的軸向速度和和油膜厚度的方法得到瞬時(shí)流量，軸向速度在噴嘴出口壓力測(cè)得后通過理論推導(dǎo)的方法得到，油膜厚度通過電導(dǎo)法得到。但是，噴嘴在脈動(dòng)壓力下的性能試驗(yàn)研究由于流量脈動(dòng)頻率高，傳感器無法達(dá)到如此高的響應(yīng)靈敏度要求，因此直接進(jìn)行試驗(yàn)有很大的困難。近年來，機(jī)械壓力噴嘴的的數(shù)值計(jì)算逐漸發(fā)展成熟。如吳偉亮等人[7[8]運(yùn)用CFD研究了各個(gè)幾何參數(shù)對(duì)機(jī)械壓力噴嘴性能的影響進(jìn)行了分析，并與Abramovich理論及其優(yōu)化理論進(jìn)行了對(duì)比，突出了CCFD在機(jī)械壓力噴嘴設(shè)計(jì)上的作用。也有運(yùn)用CFD軟件對(duì)機(jī)械壓力噴嘴的在脈動(dòng)壓力下的性能進(jìn)行研究。如文獻(xiàn)[9-10]對(duì)二維軸對(duì)稱旋轉(zhuǎn)模型對(duì)單路壓力嘴在脈動(dòng)壓力下的性能進(jìn)行了數(shù)值模擬。

不難看出，由于機(jī)械壓力噴嘴的旋流室內(nèi)區(qū)域狹小，邊界層相對(duì)厚度大，液體粘性對(duì)噴嘴旋流室內(nèi)部流場(chǎng)存在很大影響，因此研究噴嘴動(dòng)力學(xué)特性時(shí)忽略液體粘性可能對(duì)計(jì)算結(jié)果帶來較大誤差。同樣，考慮到噴嘴內(nèi)流道結(jié)構(gòu)和尺度，采用軸對(duì)稱假設(shè)也會(huì)引入誤差。在本研究工作中，完全考慮液體粘性和流道細(xì)微結(jié)構(gòu)對(duì)流場(chǎng)的影響，對(duì)機(jī)械壓力噴嘴三維內(nèi)流場(chǎng)進(jìn)行了CFD計(jì)算，分析了噴嘴瞬時(shí)流量、霧化角隨同脈動(dòng)壓力頻率和振幅的變化規(guī)律，并對(duì)瞬時(shí)流量進(jìn)行了幅頻、相頻特性分析。

1 物理模型及網(wǎng)格

本文中的研究對(duì)象為如圖2所示的機(jī)械壓力噴嘴。流體由旋流槽1流入噴嘴，經(jīng)旋轉(zhuǎn)運(yùn)動(dòng)在旋流室2中形成緊貼壁面的旋轉(zhuǎn)液膜,最后通過噴口3噴出。由于旋轉(zhuǎn)液膜具有很高的切向速度，噴口軸線附近會(huì)形成局部低壓區(qū)，從而對(duì)噴嘴出口處的環(huán)境氣體形成抽吸作用，在噴嘴軸線附近形成氣錐。因此，機(jī)械壓力噴嘴內(nèi)的流動(dòng)為狹小空間內(nèi)的氣液兩相流動(dòng)。

圖2 機(jī)械壓力噴嘴的結(jié)構(gòu)

應(yīng)用網(wǎng)格劃分軟件對(duì)計(jì)算域進(jìn)行非結(jié)構(gòu)化網(wǎng)格的劃分。在噴嘴的噴口軸線附近區(qū)域進(jìn)行網(wǎng)格加密，使數(shù)值計(jì)算結(jié)果能較好的捕捉氣液分界面。計(jì)算域從噴嘴進(jìn)口向上游延伸10倍直徑長(zhǎng)度，這樣在上游給定總壓邊界條件下，在噴嘴入口處得到接近實(shí)際的速度分布。

圖3為噴嘴的計(jì)算域。

圖3 機(jī)械壓力噴嘴的計(jì)算域

2 數(shù)學(xué)模型

2.1 VOF方法

如前所述，壓力噴嘴內(nèi)的流動(dòng)為氣液兩相分層流動(dòng)，即氣液兩相之間存在明顯的邊界，對(duì)于此類問題，相界面的捕捉是數(shù)值求解的關(guān)鍵，本文在求解中使用VOF方法。

VOF方法適用于兩相間不相互滲透，兩相間具有明顯邊界且滑移可忽略的不可壓流動(dòng)。本文中，空氣的壓力變化不大，因此可以作為不可壓流體處理。

在VOF方法中，引入液相容積比f。

其中：νl和νg分別是計(jì)算網(wǎng)格中液相和氣相的體積。因此網(wǎng)格控制體中均為液體時(shí)，f=1，均為氣體時(shí)f=0，且0≤f≤1。因?yàn)閒是可以輸運(yùn)的守恒量，滿足無源輸運(yùn)方程，即：

此方程同基本的質(zhì)量、動(dòng)量守恒方程聯(lián)立求解。計(jì)算不考慮物理相變。

在VOF方法中，根據(jù)相函數(shù)f的分布構(gòu)造相界面，從而得到氣相和液相的分布。

2.2 湍流模型

本文采用RNGk-ε模型，該模型可更好處理流線彎曲程度較大的流動(dòng)，因此能較好地模擬本文問題中出現(xiàn)的強(qiáng)渦旋流動(dòng)。

2.3 表面張力模型

在CFD計(jì)算中，表面張力被當(dāng)做是控制單元表面上的動(dòng)量源項(xiàng)，控制單元表面上微弱的變形都會(huì)引起表面張力的變化。因此氣液分界面捕捉不準(zhǔn)確、曲率計(jì)算精度低等都會(huì)引起數(shù)值誤差的放大。如果要減小這種誤差，計(jì)算中應(yīng)使用雙精度，這樣會(huì)減小曲率計(jì)算中的舍入誤差及更好的捕捉氣液分界面，但是雙精度會(huì)降低計(jì)算的經(jīng)濟(jì)性。同時(shí)考慮到對(duì)于噴嘴內(nèi)部分層流動(dòng)，表面張力相對(duì)于慣性力等可忽略不計(jì)，因此在計(jì)算中忽略表面張力。

2.4 邊界條件和初始條件

計(jì)算模擬常溫下噴嘴在脈動(dòng)壓力下的工作性能。液體工作介質(zhì)為水，且計(jì)算中水和空氣的溫度為 25℃。噴嘴出口為靜壓邊界條件，壓力等于 1atm；噴嘴入口為總壓脈動(dòng)邊界，壓力大小為P，其表達(dá)式為：

其中：Pave為總壓脈動(dòng)的平均值，計(jì)算中取值為1.0MPa；δ為壓力脈動(dòng)量的幅值；f為壓力脈動(dòng)量的頻率。

為了給上述非定常過程提供初始流場(chǎng)，同時(shí)也為了給噴嘴的脈動(dòng)壓力下的性能提供一個(gè)對(duì)比分析的依據(jù)，因此先進(jìn)行定常計(jì)算，其入口為總壓邊界條件，壓力大小為定值1.0MPa。

3 結(jié)果與討論

圖4和圖5分別是定常計(jì)算結(jié)果。在噴口附近的沿軸線橫截面的絕對(duì)壓力分布圖和氣相-液相分布圖。圖片表明噴嘴內(nèi)部存在氣液分界面清晰的空氣錐。

圖4 軸向橫截面靜壓分布圖

圖5 軸向橫截面液相-氣相分布圖

計(jì)算表明，當(dāng)噴嘴工作壓力中存在正弦變化的脈動(dòng)分量時(shí)，噴嘴瞬時(shí)流量也出現(xiàn)同樣頻率的脈動(dòng)量，如圖6所示。

圖6 脈動(dòng)壓力幅值0.1MPa，頻率200Hz工況下瞬時(shí)流量

機(jī)械壓力噴嘴在不同的壓力脈動(dòng)頻率和幅值下瞬時(shí)流量將發(fā)生改變，為更好反映這種影響，引入平均流量這個(gè)概念。平均流量是指噴嘴在在壓力脈動(dòng)的條件下工作時(shí)，在壓力脈動(dòng)一個(gè)周期內(nèi)的平均流量，其表達(dá)式為：

其中：Gave為平均流量；T為脈動(dòng)周期；t為時(shí)間。

噴嘴在脈動(dòng)壓力工作條件下，其流量特性有可能改變，會(huì)與設(shè)計(jì)流量有一定的偏差，而這會(huì)顯著的影響燃燒的性能。因此引入相對(duì)平均流量，以與設(shè)計(jì)流量進(jìn)行對(duì)比分析：

其中G0為設(shè)計(jì)流量，也就是噴嘴在穩(wěn)定的額定工作壓力下的流量，在這里即為穩(wěn)態(tài)計(jì)算得到的流量。

由圖7知，脈動(dòng)壓力下噴嘴的平均流量與設(shè)計(jì)流量相差在5%之內(nèi)，壓力脈動(dòng)頻率和幅值對(duì)其平均流量存在一定的影響，且這種影響是非線性的。在200Hz附近，噴嘴的平均流量較多的偏離設(shè)計(jì)流量，且隨著脈動(dòng)壓力幅值的增大，噴嘴的平均流量更多得偏離設(shè)計(jì)流量。

圖7 噴嘴的相對(duì)平均流量

霧化角是除流量外噴嘴的又一重要特性參數(shù)。根據(jù)定義，機(jī)械壓力式噴嘴霧化角的計(jì)算公式為：

其中：Vθ、Vz分別為位于噴嘴出口貼壁處邊界層外的監(jiān)測(cè)點(diǎn)在柱坐標(biāo)下切向和軸向分速度；u、v、w分別為監(jiān)測(cè)點(diǎn)在直角坐標(biāo)系下的分速度，其中w為噴嘴的噴射方向。如圖8所示，當(dāng)噴嘴工作壓力中存在正弦變化的脈動(dòng)分量時(shí)，噴嘴瞬時(shí)霧化角也表現(xiàn)出同樣頻率的脈動(dòng)。

為更好了解噴嘴脈動(dòng)壓力下的響應(yīng)特性，了解噴嘴的流量振幅隨脈動(dòng)壓力頻率的變化，定義相對(duì)于設(shè)計(jì)流量的相對(duì)振幅為：

其中：GA為流量振幅。

圖8 脈動(dòng)壓力幅值0.1MPa，頻率200Hz下瞬時(shí)霧化角

圖9顯示，在不同脈動(dòng)壓力幅值下，機(jī)械壓力噴嘴呈現(xiàn)出相似的幅頻特性。在100Hz～400Hz頻率內(nèi)，機(jī)械壓力噴嘴的脈動(dòng)幅值迅速增大，并在400Hz附近達(dá)到最大值；超過 400Hz，隨著脈動(dòng)頻率的增大，脈動(dòng)幅值緩慢下降。流量脈動(dòng)幅值隨著壓力脈動(dòng)幅值的增加而單調(diào)增加。

圖9 機(jī)械壓力噴嘴的流量幅頻特性曲線

圖10 機(jī)械壓力噴嘴的流量相頻特性曲線

如圖10所示，噴嘴的入口壓力與出口流量之間存在相位差，具體是指同一時(shí)刻入口脈動(dòng)壓力與出口瞬時(shí)流量間相位角的差值，該差值反映了出口瞬時(shí)流量變化滯后于脈動(dòng)壓力變化的程度，其數(shù)值與管道長(zhǎng)度有關(guān)。對(duì)于多個(gè)噴嘴的燃燒裝置，可通過計(jì)算得到各個(gè)噴嘴的相位差，進(jìn)而調(diào)節(jié)和控制整個(gè)噴霧系統(tǒng)。由圖10還可知，相位差隨頻率的增大而增加，且近乎線性關(guān)系，脈動(dòng)壓力的幅值對(duì)相位差影響不大。

4 結(jié)論

1）在正弦脈動(dòng)壓力下，機(jī)械壓力噴嘴的瞬時(shí)流量和霧化角也呈現(xiàn)出相同頻率的正弦脈動(dòng)。

2）噴嘴的平均流量與設(shè)計(jì)流量相差在5%以內(nèi)。壓力脈動(dòng)頻率和幅值對(duì)其時(shí)均流量都存在一定的影響，且這種影響是非線性的。

3）在不同脈動(dòng)壓力幅值下，機(jī)械壓力噴嘴呈現(xiàn)相似的幅頻特性。在100Hz～400Hz頻率內(nèi)，脈動(dòng)幅值隨著頻率的增大而增大；當(dāng)頻率超過400Hz時(shí)，隨頻率的增加，脈動(dòng)幅值不斷減少。流量脈動(dòng)幅值隨壓力脈動(dòng)幅值的增加而單調(diào)增加。

4）出口流量與入口壓力的相位差隨頻率的增大而增加，且近乎線性關(guān)系，此時(shí)脈動(dòng)壓力的幅值對(duì)相位差影響不大。

[1]Huang Y,Yang V.Dynamics and stability of lean-premixed swirl-stabilized combustion [J].Progress in Energy and Combustion Science,2009,35(4):293-364.

[2]Lefebvre A H,Ballal D R,Bahr D W.Gas turbine combustion:alternative fuels and emissions[M].Boca Raton,FL:CRC Press,2010.

[3]Bazarov V G.Liquid Injector Dynamics [M].Moscow,1979.

[4]楊立軍,張向陽,高芳,等.液體噴嘴在脈動(dòng)壓力下的性能數(shù)值模擬[J].航空動(dòng)力學(xué)報(bào),2004 (19):866-872.

[5]Khil T,Kim S,Cho S,et al.Quantifying the variation of the mass flow rate generated in a simplex swirl injector by pressure fluctuation[J].AIAA Paper,2008,4849:2008.

[6]Khil T,Kim S,Cho S,et al.Quantification of the transient mass flow rate in a simplex swirl injector[J].Measurement Science and Technology,2009,20(7):075405.

[7]徐剛,吳偉亮.旋流式噴嘴內(nèi)流場(chǎng)的數(shù)值模擬 [J].能源技術(shù),2008,29(3):129-132.

[8]Hinckel J N,Villa Nova H F,Bazarov V G.CFD analysis of swirl atomizers [C]// 44th AIAA/ASME/SAE/ASEE Joint Propulsion Conference and Exhibit,2008.

[9]富慶飛,楊立軍,屈元元,等.敞口型液體機(jī)械壓力噴嘴在脈動(dòng)壓力下的性能數(shù)值模擬[J].Journal of Aerospace Power,2010,25(9):497-500.

[10]Rizk N K,Lefebvre A H.Internal flow characteristics of simplex swirl atomizers[J].Journal of Propulsion and Power,1985,1(3):193-199.

Research on Performance of Simplex Swirl Atomizer Under Pulsating Pressure Condition

SHEN Yan-cheng,WU Wei-liang
(Shanghai Jiao Tong University Key Laboratory of Power Machinery and Ministry Education,Shanghai 200240,China)

Numerical simulation of internal flow of the simplex swirl atomizer is carried out by using VOF model.The change of spray angle and flow mass is analyzed under pulsating pressure condition.The results show the spray angle and flow mass fluctuate at the same frequency with the pulsating pressure.The average flow mass approximates the designed flow mass.The amplitude of the mass flow first increases and then decreases with the increase of the frequency,and reach the maximum value at 400Hz.The phase difference between outlet flow mass and inlet pressure increases linearly with the pressure frequency increasing,and the fluctuating pressure amplitude has little effect on the phase difference.

simplex-swirl atomizer; pressure fluctuation ; spray angle; flow mass

V434

A

沈彥成（1988-），男，碩士研究生。研究方向：燃?xì)廨啓C(jī)燃燒器性能。