汪鳳山，毛曉芳，虞育松，張 榛，王 平

(1.北京控制工程研究所，北京100190;2.北京交通大學(xué)，北京100044)

雙組元離心式噴注器霧化性能的大渦模擬數(shù)值研究

汪鳳山1，毛曉芳1，虞育松2，張 榛1，王 平1

(1.北京控制工程研究所，北京100190;2.北京交通大學(xué)，北京100044)

基于氣液兩相體積混合分?jǐn)?shù)建立某型雙組元離心式噴注器內(nèi)部流場(chǎng)及霧化場(chǎng)的數(shù)學(xué)模型，并采用大渦模擬(LES)方法對(duì)其霧化性能進(jìn)行研究.計(jì)算結(jié)果表明，噴注器外路噴射速度低于內(nèi)路噴射速度，內(nèi)外兩路同時(shí)噴射時(shí)，兩路相互影響，霧化錐角增加;隨著噴注器出口直段長(zhǎng)度的增加，內(nèi)外路噴射速度和流量系數(shù)均呈單調(diào)下降趨勢(shì);隨著噴注壓降的增加，噴注器流量系數(shù)和出口噴射速度隨之增加，而霧化粒度不斷降低，霧化錐角不斷增大.

同軸離心式噴注器;混合分?jǐn)?shù)模型;大渦模擬;數(shù)值模擬;霧化粒度;霧化錐角

隨著航天科技的發(fā)展，新一代大型、長(zhǎng)壽命衛(wèi)星對(duì)雙組元姿控發(fā)動(dòng)機(jī)的穩(wěn)態(tài)及脈沖工作性能要求提出了更高的要求.液體火箭發(fā)動(dòng)機(jī)工作過(guò)程涉及到流動(dòng)、霧化、燃燒及傳熱等多個(gè)過(guò)程，其中推進(jìn)劑霧化是直接影響推力器比沖、可靠性與壽命的關(guān)鍵過(guò)程.為提高發(fā)動(dòng)機(jī)的工作性能，通常需要對(duì)推進(jìn)劑霧化單元(即噴注器)進(jìn)行良好設(shè)計(jì)，尤其要對(duì)推進(jìn)劑的霧化性能進(jìn)行詳細(xì)的研究.

液體霧化現(xiàn)象包括液膜形成、液膜破碎和霧滴形成，是一種復(fù)雜的物理過(guò)程，工質(zhì)物性、噴注器結(jié)構(gòu)設(shè)計(jì)、出口背壓條件及入口壓力等均對(duì)其霧化性能產(chǎn)生影響.一直以來(lái)，國(guó)內(nèi)外學(xué)者試圖通過(guò)采用流體力學(xué)原理建立對(duì)霧化過(guò)程的數(shù)學(xué)描述，并建立多種數(shù)學(xué)模型[1].目前，針對(duì)液體在入口壓力作用下的霧化現(xiàn)象，通常采用歐拉近似方法與歐拉—拉格朗日迭代方法(即ALE方法)兩種方數(shù)值方法進(jìn)行仿真分析.第一種方法完全忽略噴霧粒子的初速和在流場(chǎng)中的受力，液滴分布完全由氣流決定，通過(guò)平均計(jì)算，以擬流體的方法計(jì)算噴霧過(guò)程中的氣流和噴霧粒子的運(yùn)動(dòng)規(guī)律.該方法的優(yōu)點(diǎn)是計(jì)算工作量較小，速度快，能夠在較短的時(shí)間內(nèi)獲得噴霧粒子的運(yùn)動(dòng)規(guī)律，特別適用于計(jì)算粒徑較小、初速較低的噴霧輸運(yùn)現(xiàn)象，如飽和蒸氣運(yùn)動(dòng)、煤粉輸運(yùn)過(guò)程等.第二種方法中噴霧粒子的運(yùn)動(dòng)不僅與初速有關(guān)，而且還受氣動(dòng)力的影響，運(yùn)動(dòng)規(guī)律通過(guò)拉格朗日方法進(jìn)行描述，更符合噴霧過(guò)程的實(shí)際情況，能夠捕捉到霧化過(guò)程中的一些細(xì)節(jié)現(xiàn)象，精度也較高，因此在實(shí)際中也應(yīng)用較廣.

針對(duì)雙組元液體火箭發(fā)動(dòng)機(jī)，噴注器出口推進(jìn)劑霧滴速度通常較高，且在氣流中的受力情況復(fù)雜，因此，不宜采用歐拉近似方法對(duì)其霧化過(guò)程進(jìn)行仿真計(jì)算，通常采用歐拉-拉格朗日迭代方法.但由于噴霧粒子尺寸較小，受氣流場(chǎng)中湍流運(yùn)動(dòng)關(guān)系密切，而通常采用的雷諾應(yīng)力湍流模型空間分辨率較低，很難得到噴霧場(chǎng)的具體細(xì)節(jié).如20世紀(jì)90年代，Jeng等[2-4]運(yùn)用歐拉-拉格朗日迭代方法追蹤壓力旋流噴嘴的氣液分界面，同時(shí)運(yùn)用數(shù)值模擬和實(shí)驗(yàn)相結(jié)合的方法研究了結(jié)構(gòu)參數(shù)對(duì)液膜厚度和噴霧錐角的影響.周立新等[5]采用基于單流體模型及氣液兩相混合分?jǐn)?shù)概念對(duì)離心式噴嘴內(nèi)部流場(chǎng)進(jìn)行數(shù)值計(jì)算，并與Jeng的ALE方法進(jìn)行對(duì)比.但由于采用了雷諾應(yīng)力湍流模型，未能給出噴霧場(chǎng)的具體分布規(guī)律.

大渦模擬(LES)方法是通過(guò)將耗散尺度的湍流脈動(dòng)過(guò)濾掉，只求解大尺度湍流脈動(dòng)的一種新型納維-斯托克斯(N-S)方程數(shù)值求解方法.與傳統(tǒng)的雷諾時(shí)均模型(RANS)方法相比，具有分辨率高、普適性高、更能真實(shí)反映流動(dòng)微觀結(jié)構(gòu)等優(yōu)點(diǎn)，但計(jì)算量要比雷諾時(shí)均模型方法的要大.該方法能夠得到較為精細(xì)和準(zhǔn)確的流場(chǎng)細(xì)節(jié)，為通過(guò)采用數(shù)值計(jì)算的方法獲得噴霧場(chǎng)的細(xì)節(jié)提供了可能.隨著近年來(lái)計(jì)算機(jī)運(yùn)行速度和計(jì)算能力的迅速進(jìn)步，采用大渦模擬方法已經(jīng)可以處理部分工程問(wèn)題的湍流流動(dòng)，在國(guó)內(nèi)外噴霧學(xué)界已經(jīng)得到了越來(lái)越廣泛的應(yīng)用.

本文通過(guò)采用大渦模擬數(shù)值方法對(duì)某型雙組元離心式噴注器的噴霧流場(chǎng)進(jìn)行了數(shù)值仿真分析，得到了不同條件下噴注器的霧化性能參數(shù)和液滴分布特性，探索噴注器霧化特性數(shù)值求解的可行性與方法，所得結(jié)果可為未來(lái)雙組元噴注器的設(shè)計(jì)提供參考.

1 數(shù)值仿真模型與算法

1.1 氣液兩相流數(shù)學(xué)模型

針對(duì)噴注器出口噴霧區(qū)域內(nèi)同時(shí)存在氣液兩種不相容工質(zhì)，氣液兩相之間存在界面的情況.為對(duì)氣液兩相共存情況進(jìn)行描述，通過(guò)引入體積分?jǐn)?shù)的概念，可實(shí)現(xiàn)采用一組納維-斯托克斯(N-S)方程描述氣液兩相及相界面，這種方法被稱之為體積分?jǐn)?shù)法(VOF)，單元體積分?jǐn)?shù)εk的計(jì)算公式為

通過(guò)引入體積分?jǐn)?shù)，對(duì)于液體噴霧特性可用如下基于壓力、速度和體積分?jǐn)?shù)的N-S方程描述:△

對(duì)于牛頓流體，在斯托克斯假設(shè)下:

由公式可得:

式中，V為流體速度矢量，ρ為流體密度，p為流體各向同性壓力，F(xiàn)為體積力(表面張力)，g為重力，μ為流體動(dòng)力黏度，τ為剪應(yīng)力，I=δij為克羅內(nèi)克符號(hào)，S為速度梯度張量.對(duì)于牛頓流體，粘性系數(shù)為常數(shù)，則:

考慮到:

將方程(6)代入方程(5)，則有:△

相界面的跟蹤是不相溶多相流數(shù)值模擬的關(guān)鍵問(wèn)題之一，通過(guò)基于異相流體混合分?jǐn)?shù)的假設(shè)，VOF方法可以建立每一相流體體積分?jǐn)?shù)的輸運(yùn)方程，對(duì)相界面進(jìn)行跟蹤:

式中，Sεk是描述相間質(zhì)量轉(zhuǎn)移的物理量，由于霧化過(guò)程模擬不考慮相變，因此相間不存在質(zhì)量轉(zhuǎn)移，即Sεk=0.對(duì)于一定體積分?jǐn)?shù)下的密度和黏度滿足:

式中，下標(biāo)l和g分別代表液相和氣相.

氣、液相界面結(jié)構(gòu)在流場(chǎng)作用下變得非常復(fù)雜，精確求解三維空間的界面曲率和表面張力相當(dāng)困難.Brackbill等[6]利用較為簡(jiǎn)潔的連續(xù)表面張力模型(CSF)解決了表面曲率和表面張力求解的問(wèn)題，具體的表達(dá)式為:

1.2 方程的數(shù)值解法

N-S方程的基本特點(diǎn)是非線性和耦合性，在一般情況下不能用解析方法得到封閉形式的解，因此需要采用一定的數(shù)值解法對(duì)N-S方法進(jìn)行求解.在本文中，通過(guò)采用LES方法來(lái)完成對(duì)控制方程的數(shù)值求解.

LES方法求解N-S方程的思想是對(duì)控制方程進(jìn)行濾波處理，將流場(chǎng)中的大小湍流渦團(tuán)進(jìn)行分離.大尺度湍流渦由控制方程直接求解;小尺度湍流渦通過(guò)引入亞網(wǎng)格尺度(SGS)湍流模型進(jìn)行模擬求解.由于小尺度湍流渦結(jié)構(gòu)比大尺度湍流渦具有更好的各向同性，因此LES對(duì)SGS模型的要求不像雷諾應(yīng)力模型(RANS)那樣高.當(dāng)大渦模擬的網(wǎng)格間距逐漸減小時(shí)，SGS黏度趨于零，預(yù)測(cè)結(jié)果逐漸逼近直接數(shù)值模擬(DNS)結(jié)果.

在采用LES數(shù)值方法對(duì)N-S方程求解過(guò)程中，為了使方程封閉，需要引入模型對(duì)該亞網(wǎng)格尺度應(yīng)力項(xiàng)進(jìn)行計(jì)算.目前，使用最廣泛且最簡(jiǎn)單的亞網(wǎng)格尺度應(yīng)力模型是Smagorinsky渦黏性模型[7]，該模型通過(guò)建立未知亞網(wǎng)格尺度應(yīng)力與已知的變形速率張量之間的簡(jiǎn)單函數(shù)關(guān)系來(lái)實(shí)現(xiàn)對(duì)亞網(wǎng)格尺度應(yīng)力項(xiàng)的模擬.假設(shè)亞網(wǎng)格應(yīng)力張量正比于平均應(yīng)力張量，則:

式中，δij為克羅內(nèi)克符號(hào)，是已知的應(yīng)變張量，υt為湍流渦黏度(m2/s)，為濾波尺度(定義為(△χ△y△z)1/3).LES的湍流渦黏度υt∝uλ，其中，u

式中，k是表面曲率，σ是表面張力系數(shù)，n是垂直表面的法向單位向量，εl為液相體積分?jǐn)?shù)，其中:是流動(dòng)速度，λ為最大的耗散湍流渦泰勒尺度(Taylor尺度).

通過(guò)對(duì)上述微分方程組離散化，將連續(xù)變化形式描述的微分方程轉(zhuǎn)化為離散形式的代數(shù)方程，采用有限體積法在計(jì)算機(jī)上進(jìn)行數(shù)值求解.

2 仿真分析

2.1 噴注器的幾何結(jié)構(gòu)

本文研究的噴注器采用雙旋渦離心噴嘴結(jié)構(gòu)，該結(jié)構(gòu)具有霧化效果好，噴霧場(chǎng)分布均勻，以及便于在燃燒室壁面組織液膜冷卻等優(yōu)點(diǎn)，如圖1所示.為了盡量減小數(shù)值計(jì)算的網(wǎng)格數(shù)量，取網(wǎng)格平均間距為100μm，總網(wǎng)格單元數(shù)約5×106，噴嘴外區(qū)域的張角根據(jù)試算預(yù)先確定.

圖1 雙旋渦離心噴注器工作原理圖Fig.1 Principle of coaxial centrifugal injector

2.2 邊界條件

在噴注器入口處根據(jù)給定壓力入口邊界條件，湍動(dòng)能和耗散率根據(jù)經(jīng)驗(yàn)公式給出，出口設(shè)置壓力出口邊界條件，壁面設(shè)置無(wú)滑移固體邊界條件，近壁函數(shù)采用標(biāo)準(zhǔn)壁面函數(shù)法求解.

3 仿真結(jié)果與分析

圖2為噴注器在不同時(shí)刻的噴霧形態(tài)數(shù)值仿真結(jié)果，從圖中可以看出，液體離開噴嘴后，液膜迅速分裂為數(shù)根細(xì)長(zhǎng)的液絲，隨后連續(xù)液絲會(huì)分裂為較長(zhǎng)的離散液滴，從而完成整個(gè)霧化過(guò)程.對(duì)于外路噴嘴，由于出口液滴噴射速度相對(duì)較低(約4m/s)，其霧化過(guò)程相對(duì)較慢;而對(duì)于內(nèi)路噴嘴，由于出口液滴噴射速度較高(約9m/s)，其霧化過(guò)程進(jìn)度相對(duì)較快，但噴霧形態(tài)具有明顯的周向不對(duì)稱性，噴霧邊緣區(qū)域存在較長(zhǎng)的連續(xù)液絲;對(duì)于兩路同時(shí)噴射的情況，內(nèi)外噴嘴出口液滴相互之間發(fā)生強(qiáng)烈的相互作用，霧化速度介于內(nèi)外路獨(dú)立噴射之間，且基本消除了內(nèi)路單獨(dú)噴射時(shí)的周向不對(duì)稱現(xiàn)象.

圖2 噴注器不同時(shí)刻的噴霧形態(tài)Fig.2 Shape of the propellant spray varied with time

圖3 為噴注器在不同時(shí)刻噴霧錐角變化示意圖.從圖中可以看出，對(duì)于噴注器，大概經(jīng)過(guò)4ms左右，其噴霧錐已基本成形.由于外路速度較低，其霧化錐角與實(shí)際數(shù)值(110°)相差較大，而內(nèi)路霧錐角度卻略小于實(shí)際數(shù)值(90°);內(nèi)外兩路同時(shí)噴射時(shí)，由于同向旋流疊加的加強(qiáng)作用，其霧化角度比內(nèi)外路單獨(dú)噴射的時(shí)候均略有增加.

圖3 噴注器霧化錐角隨時(shí)間的變化Fig.3 Spray angle varied with time

圖4 為在0.5MPa壓差下，不同出口直段長(zhǎng)度下噴注器噴射速度與流量系數(shù)的變化示意圖.從圖中可以看出，隨著出口直段長(zhǎng)度的增加，內(nèi)外路噴射速度和流量系數(shù)均呈單調(diào)下降趨勢(shì).可見隨著出口直段長(zhǎng)度的增加，在摩擦力的作用下，噴注器總流阻不斷增加，從而導(dǎo)致噴嘴旋流速度不斷降低，進(jìn)而影響到噴注器噴霧特性.

圖4 內(nèi)外路噴射速度和流量系數(shù)隨出口直段長(zhǎng)度的變化Fig.4 Spray velocity ejected from the flow passage varied with the injector exit length

圖5 為不同噴注壓差下內(nèi)外路噴射速度和流量系數(shù)的變化.從圖中可以看出，隨著噴注壓差的增大，噴注器噴射速度和流量系數(shù)也隨之增大.

圖6為外路噴嘴霧化粒度和霧化錐角隨噴注壓降的變化.從圖中可以看出，隨著噴注壓降的增加，外路霧化粒度不斷降低，霧化錐角不斷增大.可見，隨著噴注壓降的提高，液體出口速度不斷升高，液膜與氣流之間的剪切作用更為強(qiáng)烈，液膜的不穩(wěn)定性增強(qiáng)，霧化粒度不斷降低;另外，由于噴注壓降的升高，增強(qiáng)了出口液膜的旋流速度，從而導(dǎo)致形成的霧錐角增大.

圖5 內(nèi)外路噴射速度和流量系數(shù)隨噴注壓差變化計(jì)算結(jié)果Fig.5 Spray velocity from the flow passage and flow coefficient varied with the injector pressure drop

圖6 外路霧化錐角和霧化粒度隨噴注壓差的變化Fig.6 Spray velocity and SMD ejected from the outer flow passage varied with the injector pressure drop

4 結(jié) 論

本文通過(guò)引入氣液兩相流體積分?jǐn)?shù)近似方法描述氣液兩相及相界面特性，并采用LES方法對(duì)某型雙組元離心噴注器的霧化性能進(jìn)行了數(shù)值研究，結(jié)果如下:

1)噴注器外路噴射速度低于內(nèi)路噴射速度，從而導(dǎo)致外路達(dá)到霧化穩(wěn)定的時(shí)間要比內(nèi)路更長(zhǎng);

2)內(nèi)外兩路同時(shí)噴射時(shí)，兩路相互影響，由于旋轉(zhuǎn)方向相同的緣故，噴注器霧化錐角較兩路均有所增加;

3)隨著噴注器出口直段長(zhǎng)度的增加，內(nèi)外路噴射速度和流量系數(shù)均呈單調(diào)下降趨勢(shì);

4)隨著噴注壓降的增加，噴注器流量系數(shù)和出口噴射速度隨之增加，而霧化粒度不斷降低，霧化錐角不斷增大.

對(duì)計(jì)算結(jié)果進(jìn)行分析，可見通過(guò)引入體積分?jǐn)?shù)(VOF)對(duì)氣、液兩相霧化流動(dòng)納維-斯托克斯(N-S)方程進(jìn)行簡(jiǎn)化，并采用大渦模擬(LES)的方法對(duì)其進(jìn)行數(shù)值求解，能夠較為準(zhǔn)確地捕捉到霧化過(guò)程中的細(xì)節(jié)，計(jì)算精度也較高，是一種值得推廣應(yīng)用噴注器霧化過(guò)程數(shù)值仿真方法.如能通過(guò)與實(shí)驗(yàn)測(cè)試數(shù)據(jù)進(jìn)行充分比較和分析，近一步優(yōu)化計(jì)算模型，將能為液體火箭發(fā)動(dòng)機(jī)噴注器的設(shè)計(jì)、優(yōu)化提供更準(zhǔn)確的參考依據(jù).

[1] Lefebvre A H.Atomization and sprays[D].Hemisphere Publishing.New York.1989

[2] Jeng SM，Joe M A，et al.Computational and experimental study of liquid sheet emanating from simple fuel nozzle[R].AIAA-97-0796

[3] Holtzclaw D，Sakman T，Jeng SM，et al.Investigation of flow in a simplex fuel nozzle[R].AIAA-97-2970

[4] Rizk N K，Chin J S.Comprehensive fuel nozzle model [R].AIAA-94-3278

[5] 周立新，雷凡培.離心式噴嘴內(nèi)流場(chǎng)特性的數(shù)值模擬[J].推進(jìn)技術(shù)，2002，23(6):480-484 Zhou L X，Lei F P.Numerical simulation of internal flow field of swirl nozzle[J].Journal of Propulsion Technology，2002，23(6):480-484

[6] Brackbill J U，Kothe D B，Zemach C.A continuum method for modeling surface tension[J].Journal of Computational Physics，1992，100:335-354

[7] Smagorinsky J.General circulation experiments with the primitive equations.I.the basic experiment[J]. Monthly Weather Review，1963，91(3):99-164

Num erical LES Study of Sp ray Perform ance in a Bi-Propellant Coaxial Centrifugal In jector

WANG Fengshan1，MAO Xiaofang1，YU Yusong2，ZHANG Zhen1，WANG Ping1
(1.Beijing Institute of Control Engineering，Beijing 100190，China; (2.Beijing Jiaotong Uniυersity，Beijing 100044，China)

A numerical study of spray performance in a bi-propellant coaxial centrifugal injector is proposed by using the large eddy simulation(LES)method，based on themultiphase volume fraction model of the gas and liquid.The results show that the spray velocity ejected from the inner flow passage is larger than the outer flow passage，and the spray angle increases slightly when both of the two flow passages simultaneously inject.The spray velocity and flow coefficient increase with the decrease of the injector exit length and increase of the injector pressure drop.The spray droplet size(SMD)decreases and the spray angle increases with the increases of the injector pressure drop.

coaxial centrifugal injector;multiphase fraction model;large eddy simulation(LES);numerical study;spray droplet size;spray angle

V43

A

1674-1579(2012)06-0013-05

汪鳳山(1981—)，男，高級(jí)工程師，研究方向?yàn)殡p組元推力器技術(shù);毛曉芳(1971—)，女，研究員，研究方向?yàn)楹教炱魍七M(jìn)技術(shù);虞育松(1979—)，男，副教授，研究方向?yàn)槿紵茖W(xué)與技術(shù);張 榛(1983—)，男，工程師，研究方向?yàn)殡p組元推力器技術(shù);王 平(1981—)，男，高級(jí)工程師，研究方向?yàn)殡p組元推力器技術(shù).

2011-08-10

DO I:10.3969/j.issn.1674-1579.2012.06.003