強(qiáng)迫擾動(dòng)下的射流撞擊霧化特性

李佳楠，雷凡培，楊岸龍，周立新

1. 西安航天動(dòng)力研究所 液體火箭發(fā)動(dòng)機(jī)技術(shù)重點(diǎn)實(shí)驗(yàn)室，西安 710100 2. 中國(guó)船舶集團(tuán)有限公司，北京 100097

燃燒不穩(wěn)定是液體火箭發(fā)動(dòng)機(jī)研制過程中經(jīng)常遇到的重大技術(shù)問題[1]，其重要特征表現(xiàn)為燃燒室中與聲學(xué)固有振型一致的大幅值背壓振蕩，振幅達(dá)到穩(wěn)態(tài)室壓的10%～100%，振蕩頻率從數(shù)百Hz到15 000 Hz以上[2]。燃燒不穩(wěn)定會(huì)導(dǎo)致發(fā)動(dòng)機(jī)極高加速度的機(jī)械振動(dòng)，對(duì)發(fā)動(dòng)機(jī)的機(jī)械系統(tǒng)造成破壞，還會(huì)破壞推力室內(nèi)的冷卻系統(tǒng)，導(dǎo)致超高熱流流向噴注器以及壁面，引起內(nèi)壁面的燒蝕，對(duì)發(fā)動(dòng)機(jī)造成毀滅性打擊。由于燃燒不穩(wěn)定自身復(fù)雜的耦合特性，以及推力室高溫、高壓、高熱流難以獲取數(shù)據(jù)的惡劣環(huán)境，人們對(duì)燃燒不穩(wěn)定的激發(fā)、維持機(jī)理仍沒有清晰的認(rèn)識(shí)，噴霧燃燒的各個(gè)物理化學(xué)子過程都有可能成為燃燒不穩(wěn)定激發(fā)與維持的關(guān)鍵子過程。

研究燃燒不穩(wěn)定的方法有自發(fā)激勵(lì)與強(qiáng)迫響應(yīng)2種方式[3]，自發(fā)激勵(lì)必須考慮燃燒室中所有的物理化學(xué)過程，強(qiáng)迫響應(yīng)則是通過在邊界施加擾動(dòng)，研究強(qiáng)迫擾動(dòng)下的噴霧燃燒過程。自發(fā)激勵(lì)的優(yōu)點(diǎn)是完全揭示了燃燒火焰與背壓振蕩之間的激發(fā)與反饋機(jī)理，強(qiáng)迫響應(yīng)的優(yōu)點(diǎn)則在于可以研究各種不同的子過程對(duì)于不同振幅、不同頻率的聲學(xué)壓力振蕩的響應(yīng)。相關(guān)研究表明霧化不僅影響穩(wěn)態(tài)的燃燒過程，而且與燃燒不穩(wěn)定的激發(fā)、維持有著千絲萬縷的聯(lián)系。撞擊式噴嘴的Hewitt準(zhǔn)則[4-6]表明，撞擊式噴嘴燃燒穩(wěn)定性由2個(gè)重要參數(shù)控制，一是射流直徑，二是噴射速度。通過改變噴射速度可以改變發(fā)動(dòng)機(jī)的穩(wěn)定性閾值，實(shí)現(xiàn)發(fā)動(dòng)機(jī)由穩(wěn)定向不穩(wěn)定的轉(zhuǎn)變，Hewitt準(zhǔn)則預(yù)測(cè)的穩(wěn)定性邊界與多種全尺寸發(fā)動(dòng)機(jī)的實(shí)際工作特性十分吻合[6]。Qin等[7]通過數(shù)值計(jì)算發(fā)現(xiàn)，液滴粒徑分布與空間分布對(duì)燃燒不穩(wěn)定產(chǎn)生重要影響，局部聚集的細(xì)小液滴會(huì)形成局部準(zhǔn)定容燃燒過程，并在推力室空間內(nèi)形成壓力尖峰，空間隨機(jī)產(chǎn)生的壓力尖峰最終會(huì)誘導(dǎo)出聲學(xué)振型，產(chǎn)生燃燒不穩(wěn)定，在這個(gè)過程中，噴霧是最為關(guān)鍵的環(huán)節(jié)。Bai等[8]發(fā)現(xiàn)，噴霧自激振蕩必然引起火焰振蕩以及燃料噴前壓力的振蕩，當(dāng)自激振蕩的頻率與燃燒室的固有頻率一致時(shí)，會(huì)誘發(fā)燃燒室的壓力振蕩，表明霧化的自激振蕩是發(fā)生不穩(wěn)定燃燒的非常關(guān)鍵的誘因。本文將采用強(qiáng)迫響應(yīng)的方法研究強(qiáng)迫擾動(dòng)條件下的射流撞擊霧化特性。

為了揭示燃燒不穩(wěn)定復(fù)雜的耦合機(jī)理，國(guó)內(nèi)外的研究機(jī)構(gòu)采用試驗(yàn)或者數(shù)值模擬的方法開展了聲學(xué)振蕩下的噴霧或者燃燒特性研究。Baillot等[9]試驗(yàn)研究了橫向聲場(chǎng)作用下的氣液同軸剪切噴嘴的霧化特性，當(dāng)射流位于速度波腹位置時(shí)，在非線性聲輻射壓力的作用下，射流由圓柱液體變?yōu)楸馄揭耗?，液核長(zhǎng)度減小，并且破碎更加劇烈。當(dāng)推力室內(nèi)多個(gè)噴嘴的霧化特性同時(shí)改變時(shí)，非穩(wěn)態(tài)燃燒過程會(huì)發(fā)生顯著變化，就有可能激發(fā)出燃燒不穩(wěn)定。Carpentier等[10]在Kundt管中研究了橫向聲學(xué)壓力振蕩對(duì)于單束射流的影響，當(dāng)射流位于壓力波腹位置時(shí)，射流的形態(tài)幾乎沒有發(fā)生變化。而當(dāng)射流位于速度波腹位置時(shí)，與Baillot等[9]觀察到的現(xiàn)象類似，射流變?yōu)楸馄揭耗ぁ．?dāng)射流位于壓力波腹與速度波腹之間時(shí)，射流在非對(duì)稱聲輻射壓力的作用下發(fā)生偏轉(zhuǎn)。Ficuciello等[11]開展了橫向聲學(xué)振蕩條件下氣液同軸射流的霧化試驗(yàn)，進(jìn)一步驗(yàn)證了聲學(xué)振蕩對(duì)于霧場(chǎng)的直接作用，發(fā)現(xiàn)了聲學(xué)耦合作用下的液滴聚集現(xiàn)象，并基于非線性聲學(xué)理論對(duì)試驗(yàn)觀察到的現(xiàn)象進(jìn)行了闡釋。Hardi等[12]開展了低溫液氧射流試驗(yàn)，當(dāng)存在橫向聲學(xué)振蕩時(shí)，液氧的液核長(zhǎng)度減小了近70%。為進(jìn)一步揭示聲波對(duì)于霧化過程的影響，基于氣液兩相流大渦模擬算法，Rutard等[13]開展了亞臨界條件下氣液同軸剪切噴嘴在橫向振蕩壓力場(chǎng)下的霧化特性數(shù)值模擬，同樣發(fā)現(xiàn)了射流扁平以及液核縮短的現(xiàn)象，并分析了橫向聲學(xué)振蕩對(duì)于噴霧動(dòng)力學(xué)以及液滴粒徑的影響。Dighe和Gadgil[14-16]研究了橫向聲波對(duì)于射流撞擊形成液膜破碎的影響，結(jié)果表明垂直于液膜的聲學(xué)壓力振蕩促進(jìn)了液膜的破碎，生成的液滴數(shù)目增多，液滴平均粒徑減小。Hakim等[17]采用大渦模擬研究了橫向聲學(xué)擾動(dòng)條件下低溫推進(jìn)劑的火焰特性，當(dāng)火焰位于速度波腹位置時(shí)，火焰長(zhǎng)度也會(huì)減小，并且變成扁平，同時(shí)在速度波腹位置出現(xiàn)周期性振蕩的現(xiàn)象，這與橫向壓力振蕩條件下觀察到的冷態(tài)噴霧場(chǎng)有些類似。

筆者團(tuán)隊(duì)[18]總結(jié)出了背壓振蕩影響撞擊霧化的兩種作用機(jī)制，背壓振蕩除了對(duì)噴霧燃燒產(chǎn)生直接作用外，還會(huì)通過改變噴注壓降影響噴射過程，進(jìn)而對(duì)后續(xù)的霧化、燃燒等子過程產(chǎn)生影響。當(dāng)推力室中出現(xiàn)縱向的不穩(wěn)定或者噴嘴位于橫向不穩(wěn)定的壓力波腹位置時(shí)，這種機(jī)制將會(huì)成為背壓振蕩影響霧化的主要作用機(jī)制，這也是本文重點(diǎn)關(guān)注的問題。當(dāng)噴注壓降周期性變化時(shí)，噴射速度也會(huì)發(fā)生周期性改變，此時(shí)，運(yùn)動(dòng)速度較快的流體就會(huì)追擊運(yùn)動(dòng)速度較慢的流體，從而使流體在空間呈現(xiàn)出堆積的現(xiàn)象，這就是Klystron效應(yīng)[19-20]。由系統(tǒng)動(dòng)力學(xué)分析[21-22]可以證明，燃燒室的壓力振蕩與管路中的壓力振蕩都會(huì)引起噴射速度的周期性變化，二者在幅值與相位上存在定量關(guān)系，在一定程度上可以等效處理，由于后端擾動(dòng)的研究開展較為困難，可以開展前端擾動(dòng)下的霧化試驗(yàn)。

Crane等[23-24]設(shè)計(jì)了機(jī)械擾動(dòng)裝置，研究了前端壓力擾動(dòng)對(duì)射流破碎的影響，建立了擾動(dòng)頻率、幅值與射流破碎長(zhǎng)度的關(guān)系，對(duì)低速射流破碎的Rayleigh分析作了改進(jìn)，發(fā)現(xiàn)在強(qiáng)迫擾動(dòng)的作用下，低速射流破碎機(jī)理由表面張力主導(dǎo)向Klystron效應(yīng)主導(dǎo)轉(zhuǎn)變。Chigier[25]應(yīng)用壓電式驅(qū)動(dòng)裝置對(duì)射流前端施加擾動(dòng)，同樣也發(fā)現(xiàn)了Klystron效應(yīng)，射流速度的周期性變化導(dǎo)致射流形態(tài)出現(xiàn)圓盤狀的結(jié)構(gòu)，在射流本征頻率附近的擾動(dòng)能增加液滴尺寸的均勻性?？抵覞萚26]采用水力擾動(dòng)裝置研究了有/無振蕩2種情況下氣液同軸離心式噴嘴自激振蕩的噴霧形態(tài)。在理論分析方面，Heister等[27]采用邊界元方法(BEM)研究了氣體縱向振蕩條件下的射流過程，縱向擾動(dòng)會(huì)使射流速度發(fā)生周期性變化，周期性擾動(dòng)的速度會(huì)導(dǎo)致射流形成Klystron效應(yīng)，通過不同頻率的擾動(dòng)速度分析發(fā)現(xiàn)，射流對(duì)其本征頻率處的擾動(dòng)響應(yīng)最強(qiáng)。在數(shù)值模擬方面，Srinivasan等[28-29]基于OpenFOAM，應(yīng)用VOF(Volume of Fluid)界面捕捉方法實(shí)現(xiàn)了射流破碎的數(shù)值模擬，通過在射流入口給定正弦的擾動(dòng)速度，研究了射流平均速度以及速度擾動(dòng)的頻率、幅值對(duì)射流破碎的影響。數(shù)值模擬研究表明，在相同的Strouhal數(shù)下速度擾動(dòng)振幅對(duì)射流破碎產(chǎn)生顯著影響，增加振蕩振幅和頻率都能夠增加射流的噴霧錐角。Yang和Turan[30]基于開源程序Gerris[31-32]，研究了低速與高速條件下，擾動(dòng)頻率、幅值對(duì)于射流霧化的影響。對(duì)于低頻與中頻擾動(dòng)而言，隨著擾動(dòng)頻率的增大，未受擾動(dòng)的射流長(zhǎng)度將會(huì)減小，液滴平均粒徑將會(huì)增大。對(duì)于高頻擾動(dòng)而言，霧化特性參數(shù)則變化較小，表明射流對(duì)高頻擾動(dòng)的響應(yīng)不敏感。擾動(dòng)幅值對(duì)于液滴數(shù)目以及液滴粒徑具有重要影響，擾動(dòng)幅值增大，液滴數(shù)目將會(huì)增大，而液滴平均粒徑將會(huì)減小。參照Bazarov等[33]設(shè)計(jì)的壓力擾動(dòng)裝置，Yang等[34-35]設(shè)計(jì)了前端壓力擾動(dòng)裝置，通過對(duì)噴前壓力施加擾動(dòng)，試驗(yàn)研究了噴射流量變化的Klystron效應(yīng)對(duì)撞擊霧化的影響，本文將在前人的基礎(chǔ)上繼續(xù)對(duì)強(qiáng)迫擾動(dòng)下的撞擊霧化過程開展系統(tǒng)深入的研究。

對(duì)霧化過程的研究，傳統(tǒng)上主要采用試驗(yàn)的方法開展，而與試驗(yàn)相比，數(shù)值模擬具有獨(dú)特的優(yōu)勢(shì)，可以對(duì)液膜、液絲的破碎過程進(jìn)行細(xì)致的觀察，捕獲有效的數(shù)據(jù)信息，還可以對(duì)試驗(yàn)現(xiàn)象進(jìn)行深入分析，揭示試驗(yàn)現(xiàn)象背后所蘊(yùn)含的物理本質(zhì)，本文將采用數(shù)值模擬與試驗(yàn)結(jié)合的方法開展研究。霧化過程屬于多相、多尺度的流動(dòng)過程，需要處理相界面追蹤、自由界面大變形等問題，對(duì)網(wǎng)格尺度、時(shí)間尺度要求很高，計(jì)算量很大。本文采用的數(shù)值模擬方法在處理霧化這種多相、多尺度的流動(dòng)過程中具有諸多優(yōu)勢(shì)：采用樹形自適應(yīng)加密算法對(duì)空間進(jìn)行離散，同時(shí)結(jié)合分段線性的VOF方法，在實(shí)現(xiàn)相界面精確捕捉的同時(shí)，與固定網(wǎng)格相比顯著降低了計(jì)算量，時(shí)間步長(zhǎng)也隨網(wǎng)格尺度相應(yīng)調(diào)整，實(shí)現(xiàn)了時(shí)間尺度自適應(yīng)；由高度函數(shù)結(jié)合連續(xù)表面力模型實(shí)現(xiàn)表面張力的精確求解，這在霧化過程的數(shù)值計(jì)算中尤為重要。

關(guān)于撞擊式噴嘴的霧化機(jī)理及霧化特性，國(guó)內(nèi)外已經(jīng)開展了大量的研究工作，基本掌握了其在穩(wěn)態(tài)條件下的工作特性，但對(duì)于其在非穩(wěn)態(tài)條件下的霧化特性目前還知之甚少，導(dǎo)致發(fā)動(dòng)機(jī)設(shè)計(jì)中無法全面準(zhǔn)確把握該種型式噴嘴的噴霧特性。當(dāng)發(fā)動(dòng)機(jī)在工況調(diào)節(jié)狀態(tài)下工作或出現(xiàn)燃燒不穩(wěn)定性時(shí)，因無法掌握其噴霧特性而給性能評(píng)估以及故障分析和定位造成嚴(yán)重阻礙。由于燃燒不穩(wěn)定是極其復(fù)雜的物理化學(xué)過程，本文將重點(diǎn)關(guān)注背壓振蕩影響撞擊霧化第一種機(jī)制，采用試驗(yàn)結(jié)合數(shù)值模擬的方法研究噴射速度的周期性變化對(duì)于撞擊霧化特性的影響，這對(duì)于進(jìn)一步認(rèn)識(shí)霧化在燃燒不穩(wěn)定中所起的作用以及揭示燃燒不穩(wěn)定機(jī)理都將起到積極的推動(dòng)作用。

1 試驗(yàn)系統(tǒng)

首先介紹一下試驗(yàn)系統(tǒng)，試驗(yàn)系統(tǒng)主要由高壓水路供應(yīng)系統(tǒng)、水力擾動(dòng)裝置、撞擊式噴嘴、脈動(dòng)壓力數(shù)采系統(tǒng)、高速攝影系統(tǒng)以及同步觸發(fā)裝置組成，其結(jié)構(gòu)示意圖如圖1所示。液體貯箱內(nèi)的工質(zhì)水經(jīng)過高壓氮?dú)鈹D壓之后分為兩路，一路流經(jīng)撞擊式噴嘴完成霧化過程，另一路流經(jīng)水力擾動(dòng)裝置排出。

水力擾動(dòng)裝置是產(chǎn)生噴前壓力擾動(dòng)的核心部件，其結(jié)構(gòu)示意圖如圖2所示，主要由定子與轉(zhuǎn)盤兩部分組成，定子與轉(zhuǎn)盤相互貼合。轉(zhuǎn)盤邊緣沿周向均勻布置一圈小孔，由電機(jī)帶動(dòng)轉(zhuǎn)動(dòng)，當(dāng)轉(zhuǎn)盤上的孔與定子的孔重疊時(shí)就有流體流經(jīng)擾動(dòng)裝置，當(dāng)定子上的孔被完全遮擋時(shí)，就沒有流體流經(jīng)

圖1 試驗(yàn)系統(tǒng)結(jié)構(gòu)簡(jiǎn)圖Fig.1 Sketch of test system

擾動(dòng)裝置，由此來產(chǎn)生噴前壓力的擾動(dòng)，通過調(diào)節(jié)電機(jī)轉(zhuǎn)速來改變壓力擾動(dòng)的頻率，試驗(yàn)達(dá)到的最高頻率約為3 563 Hz。

由Kistler脈動(dòng)壓力傳感器采集噴前的脈動(dòng)壓力，由XP5壓力傳感器記錄周期性變化的噴注壓降，采樣頻率均設(shè)定為20 480 Hz。由高速攝影對(duì)霧場(chǎng)進(jìn)行背光拍攝，設(shè)定高速相機(jī)的拍攝頻率為每秒拍攝20 480張圖片，圖像分辨率為512 pixel×512 pixel。通過同步觸發(fā)裝置實(shí)現(xiàn)脈動(dòng)壓力測(cè)量與霧場(chǎng)拍攝的同步，同步測(cè)量的原理為脈動(dòng)壓力傳感器對(duì)脈動(dòng)壓力不間斷采集，由信號(hào)發(fā)生器產(chǎn)生一個(gè)階躍信號(hào)觸發(fā)高速攝影工作。高速攝影觸發(fā)的時(shí)間精度為2 μs，DG535信號(hào)發(fā)生器的時(shí)間精度為ns量級(jí)，而數(shù)采系統(tǒng)采樣的時(shí)間精度為10 μs，根據(jù)木桶原理整個(gè)同步觸發(fā)裝置的時(shí)間精度應(yīng)為10 μs，目前試驗(yàn)所做的最高擾動(dòng)頻率約為3 563 Hz，10 μs的時(shí)間精度等級(jí)可以滿足試驗(yàn)需求。

試驗(yàn)采用的撞擊式噴嘴結(jié)構(gòu)示意圖如圖3所示，兩束圓柱射流軸線的交點(diǎn)定義為撞擊點(diǎn)，液體工質(zhì)加壓噴射之后在撞擊點(diǎn)處相互撞擊完成霧化

圖3 撞擊式噴嘴結(jié)構(gòu)示意圖Fig.3 Sketch of impinging jet injector

過程。圖3中標(biāo)注了撞擊式噴嘴的部分結(jié)構(gòu)參數(shù)，選用了射流直徑d為1 mm與0.8 mm兩種孔徑的噴嘴，撞擊角2θ=60°，噴嘴出口到撞擊點(diǎn)的噴射距離L約為6.9 mm。

2 數(shù)值模擬方案及算例驗(yàn)證

2.1 數(shù)值模擬方案

基于開源程序Gerris[31-32]實(shí)現(xiàn)兩股射流撞擊霧化過程的數(shù)值模擬，Gerris已在多種不同結(jié)構(gòu)形式的噴嘴霧化[36-38]中得到應(yīng)用，將采用的數(shù)值模擬方案描述如下。由八叉樹結(jié)構(gòu)形式的笛卡兒網(wǎng)格對(duì)空間進(jìn)行離散，結(jié)合自適應(yīng)加密算法實(shí)現(xiàn)特定區(qū)域的局部加密。由有限體積法直接數(shù)值求解不可壓Navier-Stokes方程組，由分段線性的VOF方法對(duì)氣液界面進(jìn)行重構(gòu)。表面張力采用Brackbill等[39]提出的CSF(Continuum Surface Force)方法，將相界面上一定寬度區(qū)域的表面張力等效為連續(xù)的體積力，在動(dòng)量方程中加入這一項(xiàng)。湍流的處理采用隱式大渦模擬方法[40-41]，由數(shù)值黏性充當(dāng)亞格子黏性模擬亞網(wǎng)格尺度的渦耗散。基于Linux系統(tǒng)，以MPI(Message Passing Interface)為并行平臺(tái)實(shí)現(xiàn)三維兩相流的并行計(jì)算。由于采用了動(dòng)態(tài)的自適應(yīng)網(wǎng)格，為了提高計(jì)算效率，采用了動(dòng)態(tài)負(fù)載平衡技術(shù)。

圖4 數(shù)值模擬計(jì)算域Fig.4 Computational domain of numerical simulations

求解的三維不可壓Navier-Stokes方程組為

(1)

(2)

(3)

采用經(jīng)典的VOF方法[42]對(duì)氣液界面進(jìn)行重構(gòu)，在流場(chǎng)中引入體積分?jǐn)?shù)α的概念，表征某一網(wǎng)格內(nèi)第一相流體體積與網(wǎng)格體積之比：α=1表示該網(wǎng)格內(nèi)充滿第一相流體；α=0表示網(wǎng)格內(nèi)充滿第二相流體；0<α<1表示網(wǎng)格內(nèi)存在兩相流體，故存在自由界面。由于體積分?jǐn)?shù)的引入，密度的對(duì)流方程式(1)可以等效為體積分?jǐn)?shù)α的對(duì)流方程

(4)

流體的密度及黏度系數(shù)由網(wǎng)格內(nèi)2種流體的體積分?jǐn)?shù)加權(quán)計(jì)算

ρ=ρlα+ρg(1-α)

(5)

μ=μlα+μg(1-α)

(6)

式中：ρl與ρg分別為液相與氣相的密度；μl與μg分別為液相與氣相的黏度系數(shù)。

VOF方法的一個(gè)優(yōu)點(diǎn)是不需要對(duì)界面的破裂、融合等現(xiàn)象作特殊的處理，因?yàn)檫@些拓?fù)浣Y(jié)構(gòu)的改變都隱式地包含在VOF方法中[43]。因此，VOF方法非常適合應(yīng)用于存在諸多破碎、聚合現(xiàn)象的霧化過程的計(jì)算。

采用經(jīng)典的時(shí)間分裂投影方法對(duì)控制方程進(jìn)行離散，離散后的方程為

(7)

(8)

(9)

(10)

同時(shí)需要求解泊松方程

(11)

2.2 算例驗(yàn)證

關(guān)于Gerris計(jì)算穩(wěn)態(tài)撞擊霧化的算例驗(yàn)證工作在相關(guān)文獻(xiàn)[36,44]中已經(jīng)開展過，本文將進(jìn)一步驗(yàn)證建立的數(shù)值模擬方案計(jì)算非穩(wěn)態(tài)撞擊霧化過程的有效性。對(duì)試驗(yàn)工況進(jìn)行數(shù)值模擬的關(guān)鍵是獲得噴嘴出口的周期性變化速度，以此作為數(shù)值計(jì)算的邊界條件?，F(xiàn)有的流量計(jì)都是穩(wěn)態(tài)流量計(jì)，無法獲得高頻脈動(dòng)的噴射速度，將通過流量壓降的關(guān)系計(jì)算高頻脈動(dòng)的噴射速度。驗(yàn)證算例采用直徑為0.8 mm的噴嘴，由XP5表壓傳感器記錄得到隨時(shí)間周期性變化的噴注壓降，由于測(cè)量點(diǎn)與噴嘴出口距離較近，忽略流體慣性產(chǎn)生的時(shí)滯，并假定噴嘴的流量系數(shù)不隨噴注壓降的變化而變化。由壓降與流量的計(jì)算公式(12)得到噴嘴的周期性變化流量。由式(13)換算成噴嘴出口的周期性變化速度u。要得到周期性變化的噴射速度還需要標(biāo)定噴嘴的流量系數(shù)Cd，由式(12)可以看出噴嘴流量的平方與噴注壓降呈線性關(guān)系，可以記錄一系列一一對(duì)應(yīng)的流量與壓降，通過曲線擬合得到流量系數(shù)。

(12)

(13)

XP5表壓傳感器記錄得到的周期性變化的噴注壓降如圖5所示，經(jīng)過曲線擬合可近似認(rèn)為噴注壓降為正弦形式的振蕩。噴嘴流量系數(shù)的標(biāo)定曲線如圖6所示，與上述分析一致，流量的平方與噴注壓降近似呈線性關(guān)系，經(jīng)過擬合可以得到噴嘴的流量系數(shù)Cd≈0.555。最終得到的噴嘴出口噴射速度隨時(shí)間的變化如圖7所示，噴嘴出口的周期性變化速度可以寫成式(14)的形式，其中平均噴射速度u0=16.1 m/s，速度脈動(dòng)幅值u′=0.186，擾動(dòng)頻率f=1 257 Hz。

u=u0(1+u′sin(2πft))

(14)

圖5 壓降信號(hào)與擬合曲線Fig.5 Pressure drop signals and fitted line

圖6 流量與壓降的對(duì)應(yīng)關(guān)系Fig.6 Relationship between mass flow rate and pressure drop

圖7 周期性變化的噴射速度Fig.7 Periodically varying injection velocity

由上述得到的周期性變化速度作為數(shù)值模擬的邊界條件，基于2.1節(jié)建立的數(shù)值模擬方案開展數(shù)值模擬。首先進(jìn)行網(wǎng)格無關(guān)性驗(yàn)證工作，分別設(shè)定網(wǎng)格最高加密等級(jí)為7級(jí)、8級(jí)與9級(jí)，對(duì)應(yīng)的最小網(wǎng)格尺度分別為234 μm、117 μm與58.6 μm，計(jì)算得到的撞擊式霧場(chǎng)分別如圖8所示。當(dāng)網(wǎng)格分辨率較低時(shí)，數(shù)值模擬捕捉到的液滴尺寸較大，無法準(zhǔn)確捕捉液膜的波動(dòng)破碎過程，同樣無法辨識(shí)霧化的Klystron效應(yīng)。隨著網(wǎng)格分辨率的提高，計(jì)算得到的霧場(chǎng)結(jié)構(gòu)越來越精細(xì)。當(dāng)網(wǎng)格最高加密等級(jí)為9級(jí)的時(shí)候，數(shù)值模擬能夠捕捉到霧場(chǎng)中液膜、液絲、液滴等精細(xì)結(jié)構(gòu)，同時(shí)也捕捉到了非穩(wěn)態(tài)霧化中周期性液滴群聚集現(xiàn)象。因此最高加密等級(jí)為9級(jí)可以滿足霧場(chǎng)空間發(fā)展對(duì)網(wǎng)格分辨率的要求，這與文獻(xiàn)[45]中網(wǎng)格無關(guān)性驗(yàn)證的結(jié)論一致，在后續(xù)的數(shù)值模擬中設(shè)定網(wǎng)格的最高加密等級(jí)為9級(jí)。

計(jì)算得到了兩束射流從噴射到撞擊形成液膜，液膜波動(dòng)破碎形成弓形液絲，弓形液絲進(jìn)一步收縮形成液滴的完整霧化過程，整個(gè)過程與Anderson等[46]建立的三步霧化模型完全一致。數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的對(duì)比如圖9所示，圖9(a)為霧場(chǎng)正面的對(duì)比，圖9(b)為霧場(chǎng)側(cè)面的對(duì)比。從宏觀上看，計(jì)算得到的霧場(chǎng)結(jié)構(gòu)與試驗(yàn)拍攝的霧場(chǎng)結(jié)構(gòu)是非常接近的。霧場(chǎng)正面形成以撞擊點(diǎn)為頂點(diǎn)的具有一定噴霧角的霧化區(qū)域，霧場(chǎng)側(cè)面液膜在撞擊波的作用下波動(dòng)破碎，并且在周期性變化的噴射速度作用下調(diào)制出Klystron效應(yīng)，霧場(chǎng)下游出現(xiàn)周期性脫落的大尺度液絲，在3.1節(jié)將重點(diǎn)分析強(qiáng)迫擾動(dòng)下的撞擊霧化特性。圖10展示了計(jì)算得到的非穩(wěn)態(tài)噴霧場(chǎng)結(jié)構(gòu)，根據(jù)試驗(yàn)以及數(shù)值模擬結(jié)果，其噴霧場(chǎng)可以大致劃分為2個(gè)區(qū)域，即一次霧化區(qū)域與二次霧化區(qū)域。一次霧化區(qū)域主要包含液膜破碎的區(qū)域，液膜初步破碎形成了大尺寸的液絲以及大顆粒的液滴，是相對(duì)粗糙的霧化區(qū)域。二次霧化區(qū)域則更靠近霧場(chǎng)的下游，液絲與大顆粒液滴在氣動(dòng)力以及表面張力的作用下進(jìn)一步收縮破碎形成更小尺度的球形液滴，霧化最終形成的液滴尺寸不僅取決于一次霧化形成的大顆粒液滴尺寸，還取決于二次霧化形成的小顆粒液滴尺寸。

圖9 數(shù)值模擬噴霧圖像與試驗(yàn)的對(duì)比Fig.9 Comparison of spray images between simulations and tests

圖10 數(shù)值模擬得到的非穩(wěn)態(tài)噴霧場(chǎng)Fig.10 Unsteady atomization field from numerical simulations

圖11 YZ平面的速度分布Fig.11 Distribution of velocity on plane YZ

霧場(chǎng)沿XY截面與YZ截面的體積分?jǐn)?shù)分布灰度圖分別如圖12(a)和圖12(b)所示，黑色區(qū)域的體積分?jǐn)?shù)為0，是氣相區(qū)域，白色區(qū)域表示體積分?jǐn)?shù)為1，為液相區(qū)域，氣液界面處的體積分?jǐn)?shù)0<α<1。 從兩幅圖中都可以明顯觀察到液膜的波動(dòng)特性，液膜向下游運(yùn)動(dòng)過程中，在Kelvin-Helmholtz不穩(wěn)定的作用下發(fā)生破碎。在圖12(a)中定義了液膜的破碎長(zhǎng)度與撞擊波的波長(zhǎng)λ，其中液膜破碎長(zhǎng)度是一個(gè)動(dòng)態(tài)的概念，液膜完全破碎之后，射流撞擊形成液膜的破碎長(zhǎng)度定義為撞擊點(diǎn)到液膜完全破碎位置的距離。從圖12(a)可以看出隨著與撞擊點(diǎn)距離的增加，液膜的波動(dòng)幅值在逐漸增大，達(dá)到破碎的臨界條件之后就開始破碎形成液絲。由于噴射速度是周期性變化的，速度較快的流體會(huì)追擊速度較慢的流體，數(shù)值模擬也捕捉到了射流中流體堆積形成的“鼓包”結(jié)構(gòu)。由圖12(b)可以看出，液膜下游黑色與白色區(qū)域交替出現(xiàn)，這正是由于液膜的波動(dòng)造成的，并且在越靠近下游的區(qū)域，弓形液絲的厚度也在減小，表明液膜在向下游運(yùn)動(dòng)過程中厚度是在不斷減小的，這也是造成液膜不穩(wěn)定性的另外一個(gè)因素[47]。在液膜破碎之后的下游區(qū)域，由于霧化的Klystron效應(yīng)出現(xiàn)了大尺度的弓形液體結(jié)構(gòu)。

圖12 不同截面的α灰度圖Fig.12 Grayscale images of α of different slices

圖13 液膜破碎區(qū)域速度分量的分布云圖Fig.13 Contour of velocity component of liquid sheet breakup region

圖14 液膜破碎區(qū)域渦量分量的分布云圖Fig.14 Contour of vortex component of liquid sheet breakup region

數(shù)值模擬統(tǒng)計(jì)液滴粒徑的方法是由液滴的體積進(jìn)行換算，由包含在氣液界面內(nèi)的網(wǎng)格體積與網(wǎng)格內(nèi)液相體積分?jǐn)?shù)相乘求和得到液滴的體積，假定液滴形狀為球形，由球形的體積計(jì)算公式計(jì)算得到液滴的直徑。統(tǒng)計(jì)不同網(wǎng)格分辨率條件下霧場(chǎng)液滴尺寸分布的概率密度函數(shù)(Probability Density Function, PDF)，其定義為某個(gè)粒徑區(qū)間內(nèi)的液滴概率，不同網(wǎng)格分辨率的液滴粒徑統(tǒng)計(jì)結(jié)果如圖15所示。當(dāng)網(wǎng)格加密等級(jí)較低時(shí)，數(shù)值模擬會(huì)捕捉到一些大尺寸液滴。隨著網(wǎng)格分辨率的提高，數(shù)值模擬捕捉到的液滴尺寸開始向小粒徑區(qū)域偏移。需要指出的是，當(dāng)網(wǎng)格分辨率不夠時(shí)，就會(huì)發(fā)生數(shù)值破碎而產(chǎn)生偽液滴，這種情況會(huì)隨著網(wǎng)格分辨率的提高而得到改善，但是并不能夠完全消除。如果非線性不穩(wěn)定也能產(chǎn)生飛濺的小液滴的話，由于網(wǎng)格分辨率的缺陷也會(huì)導(dǎo)致對(duì)這種小液滴捕捉不夠準(zhǔn)確，因此需要對(duì)這種小液滴進(jìn)行剔除。數(shù)值模擬的有效性應(yīng)體現(xiàn)在可以較為完整地展現(xiàn)整個(gè)霧化過程，并且數(shù)值破碎產(chǎn)生的小液滴所占的比例足夠小。假定小于4個(gè)網(wǎng)格捕捉到的液滴屬于偽液滴，使用Gerris自帶的RemoveDroplets函數(shù)將小于4個(gè)網(wǎng)格捕捉到的液滴全部剔除，因此統(tǒng)計(jì)曲線的起始點(diǎn)并不是零點(diǎn)。由圖15可以看出，計(jì)算得到的液滴粒徑分布范圍大致在0～350 μm，液滴粒徑概率密度分布的峰值在50～200 μm之間，概率密度數(shù)值基本在0.02以上。

圖15 不同網(wǎng)格分辨率的液滴尺寸分布Fig.15 Droplet diameter distributions of different mesh resolutions

3 結(jié)果分析

3.1 強(qiáng)迫擾動(dòng)下的射流撞擊霧化特性

在前面已經(jīng)提到，背壓振蕩與噴前壓力擾動(dòng)都會(huì)引起噴射速度的周期性變化，二者在幅值與相位上存在一定關(guān)系，可以通過開展噴前壓力擾動(dòng)試驗(yàn)，并結(jié)合系統(tǒng)動(dòng)力學(xué)分析將噴前壓力擾動(dòng)轉(zhuǎn)化為背壓振蕩對(duì)噴射的影響。采用水力擾動(dòng)裝置產(chǎn)生噴前壓力的擾動(dòng)，開展了前端擾動(dòng)條件下的射流撞擊霧化試驗(yàn)，試驗(yàn)采用孔徑為1 mm的噴嘴。首先對(duì)比強(qiáng)迫擾動(dòng)霧化與自然霧化之間的差別，然后著重探討強(qiáng)迫擾動(dòng)條件下的撞擊霧化特性。圖16為脈動(dòng)壓力傳感器記錄得到的典型噴前脈動(dòng)壓力p′，水力擾動(dòng)裝置迫使噴前壓力出現(xiàn)周期性振蕩，從而使噴注壓降發(fā)生周期性改變，最終導(dǎo)致噴射速度周期性變化，并調(diào)制出霧化的Klystron效應(yīng)。對(duì)該脈動(dòng)壓力信號(hào)作快速傅里葉變換(FFT)，可以得到該脈動(dòng)壓力的頻譜圖，如圖17所示，脈動(dòng)壓力的頻率為1 782 Hz。由高速攝影對(duì)存在噴前壓力擾動(dòng)與不存在噴前擾動(dòng)的噴霧場(chǎng)分別進(jìn)行拍攝，圖18(a)為不存在擾動(dòng)時(shí)的自然霧化圖像，圖18(b)為存在噴前壓力擾動(dòng)的霧化圖像。從宏觀上看，二液滴空間分布比較均勻，而對(duì)于存在擾動(dòng)時(shí)的噴霧場(chǎng)，在液膜下游則出現(xiàn)了大量液滴群聚集的弓形結(jié)構(gòu)，觀察不同時(shí)刻的高速攝影圖像可以發(fā)現(xiàn)弓形液滴群的出現(xiàn)表現(xiàn)出一定的周期性。Rayleigh準(zhǔn)則[48]指出了發(fā)生燃燒不穩(wěn)定的2個(gè)必要條件：① 存在周期性的釋熱波動(dòng)；② 壓力振蕩與釋熱波動(dòng)同相位。當(dāng)同時(shí)滿足這2個(gè)條件時(shí)，壓力振蕩就會(huì)在正反饋機(jī)制下持續(xù)放大。霧化的周期性必然導(dǎo)致周期性的燃燒釋熱，如果周期性的燃燒釋熱與周期性的壓力振蕩耦合在一起形成正反饋，則有可能調(diào)制出燃燒不穩(wěn)定，由此可以推斷，背壓振蕩調(diào)制出的周期性霧化有可能成為燃燒不穩(wěn)定的驅(qū)動(dòng)機(jī)制之一。

圖16 噴前脈動(dòng)壓力Fig.16 Pulsating pressures before injection

圖17 脈動(dòng)壓力的頻譜分布Fig.17 Spectral distribution of pressure fluctuations

圖18 撞擊式噴嘴的噴霧場(chǎng)Fig.18 Atomization field of impinging liquid jets

(15)

為獲得更多的數(shù)據(jù)信息，本文基于開源程序Gerris開展了數(shù)值模擬，重點(diǎn)關(guān)注的是噴射速度的周期性變化對(duì)于霧化特性的影響，因此，模擬強(qiáng)迫擾動(dòng)霧化的方法是按照式(14)在邊界給定隨時(shí)間周期性變化的噴射速度。設(shè)定射流的平均速度u0=31.5 m/s，速度的擾動(dòng)幅值u′=0.05，擾動(dòng)頻率f=1 782 Hz，其他結(jié)構(gòu)參數(shù)與試驗(yàn)所用的噴嘴參數(shù)一致。數(shù)值計(jì)算得到的自然噴霧場(chǎng)與強(qiáng)迫擾動(dòng)下的噴霧場(chǎng)對(duì)比如圖22所示，二者之間的差別與試驗(yàn)觀察到的霧場(chǎng)區(qū)別一致，自然噴霧場(chǎng)的液滴空間分布比較均勻，而強(qiáng)迫擾動(dòng)下的液膜下游出現(xiàn)了大尺度的弓形結(jié)構(gòu)，并且在霧場(chǎng)出現(xiàn)了局部液滴聚集的周期性現(xiàn)象。在一個(gè)周期T內(nèi)，霧場(chǎng)隨時(shí)間的變化如圖23所示，由一個(gè)周期內(nèi)的霧化過程可以對(duì)撞擊霧化的Klystron效應(yīng)形成機(jī)制進(jìn)行闡述。當(dāng)噴射速度出現(xiàn)周期性變化時(shí)，運(yùn)動(dòng)速度較快的流體追擊運(yùn)動(dòng)速度較慢的流體，在射流中出現(xiàn)“鼓包”的現(xiàn)象，如圖23(a)，之后“鼓包”向下游運(yùn)動(dòng)并在撞擊點(diǎn)位置相互撞擊，如圖23(d)。該現(xiàn)象類似于射流速度更高、直徑更大的兩束射流撞擊，撞擊之后的液膜厚度與寬度增加，如圖23(e)。 液膜變得更加不穩(wěn)定，大尺度的弓形液團(tuán)形成并在液膜下游脫落，如圖23(f)～圖23 (h)。液團(tuán)在向下游運(yùn)動(dòng)過程中快速破碎并形成大量在空間聚集的液滴群，整個(gè)物理過程在時(shí)間上表現(xiàn)出顯著的周期性特征，類似于射流的Klystron效應(yīng)，可以認(rèn)為是撞擊霧化表現(xiàn)出的Klystron效應(yīng)。

圖19 3個(gè)觀察窗的大小與位置Fig.19 Positions and sizes of three interrogation windows

圖20 自然霧化隨時(shí)間變化的灰度值與相應(yīng)的FFT結(jié)果Fig.20 Time varying grayscale values and corresponding FFT results for natural atomization

圖21 強(qiáng)迫擾動(dòng)霧化隨時(shí)間變化的灰度值與相應(yīng)的FFT結(jié)果Fig.21 Time varying grayscale values and corresponding FFT results for atomization with forced perturbations

圖22 數(shù)值模擬得到的不同噴霧場(chǎng)的對(duì)比Fig.22 Comparison of different atomization fields from numerical simulations

圖24 不同位置處速度分量隨時(shí)間的變化Fig.24 Time varying velocity component at different positions

與試驗(yàn)統(tǒng)計(jì)霧場(chǎng)信息的方法類似，在撞擊點(diǎn)下游某個(gè)位置處設(shè)定觀察窗，統(tǒng)計(jì)該區(qū)域的液滴信息。統(tǒng)計(jì)撞擊點(diǎn)下游10d與20d位置處，厚度為2d，寬度為30d區(qū)域內(nèi)的液滴數(shù)目以及Sauter平均直徑(SMD)隨時(shí)間的變化。自然霧化的液滴信息統(tǒng)計(jì)結(jié)果如圖25所示，可以看出液滴數(shù)目以及Sauter平均直徑的變化與試驗(yàn)觀測(cè)的灰度信息類似，呈現(xiàn)出隨時(shí)間無規(guī)律變化的特征。而當(dāng)噴射速度周期性變化時(shí)，觀察窗內(nèi)的液滴數(shù)目以及液滴的Sauter平均直徑則表現(xiàn)出周期性變化的特征，如圖26所示，并且周期性波動(dòng)的頻率與施加的強(qiáng)迫擾動(dòng)的頻率一致。對(duì)某個(gè)時(shí)刻兩個(gè)觀察窗內(nèi)的液滴粒徑分布進(jìn)行統(tǒng)計(jì)，統(tǒng)計(jì)結(jié)果分別如圖27(a)和圖27(b)所示。當(dāng)存在強(qiáng)迫擾動(dòng)時(shí)，小尺寸液滴所占的比例提高，大尺寸液滴所占的比例減小，表明在強(qiáng)迫擾動(dòng)的作用下，觀察窗內(nèi)的液滴平均尺寸有減小的趨勢(shì)。對(duì)一段時(shí)間內(nèi)觀察窗內(nèi)的液滴數(shù)目與液滴粒徑進(jìn)行時(shí)間平均，統(tǒng)計(jì)結(jié)果如表1所示。當(dāng)存在周期性擾動(dòng)時(shí)，觀察窗內(nèi)的液滴數(shù)目顯著增多，時(shí)間平均之后的Sauter平均直徑顯著減小，最大Sauter平均直徑也相應(yīng)減小。與自然霧化相比，強(qiáng)迫擾動(dòng)霧化的2個(gè)觀察窗內(nèi)液滴數(shù)目分別提高了6.0%與13.3%， 時(shí)間平均的Sauter平均直徑分別減小了8.8%與5.5%。由數(shù)值模擬結(jié)果可以看出，當(dāng)存在前端擾動(dòng)時(shí)，霧場(chǎng)空間局部區(qū)域生成的液滴數(shù)目顯著增多，并且液滴粒徑有減小的趨勢(shì)。由Qin等[7]的研究成果可以得知，燃燒不穩(wěn)定的激發(fā)與液滴粒徑有強(qiáng)的相關(guān)關(guān)系，液滴粒徑的減小與液滴數(shù)目的增多更有利于激發(fā)燃燒不穩(wěn)定，強(qiáng)迫擾動(dòng)下的霧化有可能成為燃燒不穩(wěn)定激發(fā)的重要中間環(huán)節(jié)。

圖25 自然霧化觀察窗內(nèi)的液滴數(shù)與SMD隨時(shí)間的變化Fig.25 Time varying droplet number and SMD of interrogation windows of natural atomization field

圖26 強(qiáng)迫擾動(dòng)霧化觀察窗內(nèi)的液滴數(shù)與SMD隨時(shí)間的變化Fig.26 Time varying droplet number and SMD of interrogation window of atomization field coupled with forced perturbations

圖27 觀察窗內(nèi)的液滴粒徑分布Fig.27 Droplet diameter distribution in interrogation windows

表1 自然霧化與強(qiáng)迫擾動(dòng)霧化觀察窗內(nèi)的液滴粒徑信息Table 1 Droplet information of two interrogation windows of natural atomization and forced atomization

3.2 擾動(dòng)頻率對(duì)霧化特性的影響

這一部分將分析擾動(dòng)頻率對(duì)于撞擊霧化的影響，通過改變電機(jī)轉(zhuǎn)速實(shí)現(xiàn)對(duì)頻率的調(diào)節(jié)，圖28給出了不同頻率的噴前脈動(dòng)壓力作快速傅里葉變換之后的頻譜圖，擾動(dòng)的頻率從1 347 Hz一直變化到3 563 Hz。圖29為對(duì)應(yīng)的不同擾動(dòng)頻率條件下的霧場(chǎng)圖像，從圖29可以看出在所研究的頻率范圍(1 347～3 563 Hz)內(nèi)，撞擊式噴嘴的霧化對(duì)擾動(dòng)都存在響應(yīng)，在強(qiáng)迫擾動(dòng)的作用下撞擊霧場(chǎng)都表現(xiàn)出周期性特征，表明撞擊式噴嘴的霧化在一個(gè)非常寬的頻帶范圍內(nèi)對(duì)擾動(dòng)都存在響應(yīng)。這也表明撞擊式噴嘴的工作狀態(tài)易受到外界的影響，并不是一種工作特別穩(wěn)定的噴注霧化單元，這種結(jié)構(gòu)形式的噴嘴非常容易產(chǎn)生燃燒不穩(wěn)定[20]。

圖28 不同脈動(dòng)壓力的頻譜分布Fig.28 Spectral distribution of different pressure fluctuations

從宏觀上看，隨著擾動(dòng)頻率的增大，空間局部聚集的弓形液滴群之間的間距在減小，表明液膜下游弓形液絲脫落的頻率也在提高，與強(qiáng)迫擾動(dòng)的頻率變化保持一致。

圖29 不同頻率下的強(qiáng)迫擾動(dòng)噴霧場(chǎng)Fig.29 Forced atomization field with different pressure fluctuation frequencies

圖30 同步測(cè)量的灰度與脈動(dòng)壓力信號(hào)Fig.30 Synchronous signals of grayscale values and fluctuating pressures

(16)

圖31 壓力振蕩與非穩(wěn)態(tài)釋熱之間的時(shí)滯Fig.31 Time delay between pressure fluctuations and unsteady heat release

假定相鄰2個(gè)脈動(dòng)壓力波形之間的相位差為360°，統(tǒng)計(jì)得到的不同頻率條件下霧場(chǎng)灰度與脈動(dòng)壓力之間的相位差如圖32所示。從圖中可以看出擾動(dòng)頻率對(duì)相位差影響很大，霧場(chǎng)相同位置處，擾動(dòng)頻率越高，灰度信息與擾動(dòng)壓力之間的相位差值也就越大。在相同的距離范圍內(nèi)，擾動(dòng)頻率越高，統(tǒng)計(jì)得到的相位差的變化區(qū)間也就越大，表明在相同的距離內(nèi)所包含的灰度變化周期更多，這也與從圖29中觀察到的弓形液滴群之間距離減小的現(xiàn)象一致。文中得到的是噴前脈動(dòng)壓力與霧場(chǎng)灰度之間的相位關(guān)系，通過系統(tǒng)動(dòng)力學(xué)分析可以得到產(chǎn)生相同流量脈動(dòng)的條件下，后端擾動(dòng)與前端擾動(dòng)之間的相位關(guān)系，這樣就可以間接得到后端壓力擾動(dòng)與脈動(dòng)灰度之間的相位關(guān)系，考慮其他子過程之后由圖31可以得到最終的相位關(guān)系。未來可進(jìn)一步考慮燃燒的影響，建立起基于周期性霧化的燃燒穩(wěn)定性分析理論體系。

圖32 不同頻率灰度值與振蕩壓力之間的相位角Fig.32 Phase angle between grayscale values and pressure fluctuations at different frequencies

3.3 擾動(dòng)幅值對(duì)霧化特性的影響

液體火箭發(fā)動(dòng)機(jī)推力室內(nèi)的室壓通常都存在小幅值的波動(dòng)，只有當(dāng)發(fā)生燃燒不穩(wěn)定時(shí)，擾動(dòng)壓力的幅值才會(huì)超出穩(wěn)態(tài)室壓的10%以上，擾動(dòng)壓力的幅值是甄別燃燒不穩(wěn)定是否發(fā)生的一個(gè)重要參數(shù)，這一部分將采用數(shù)值模擬的方法研究擾動(dòng)幅值對(duì)撞擊霧化的影響。當(dāng)推力室內(nèi)出現(xiàn)縱向振蕩時(shí)，噴嘴處于壓力波腹的位置，反壓振蕩影響霧化的主要作用機(jī)制是通過改變噴注壓降繼而改變噴射速度，從而對(duì)霧化產(chǎn)生影響，也就是本文所研究的作用機(jī)制。選取與某型推力室一階縱向不穩(wěn)定頻率接近的頻率作為參照頻率來研究擾動(dòng)幅值對(duì)于撞擊霧化的影響。根據(jù)式(14)，設(shè)定擾動(dòng)頻率為2 230 Hz，擾動(dòng)幅值的大小u′依次為0.01，0.05，0.1，0.2與0.5，計(jì)算得到不同擾動(dòng)幅值條件下的噴霧場(chǎng)分別如圖33所示。由圖33(a)可以看出，當(dāng)擾動(dòng)幅值很小時(shí)，霧場(chǎng)空間并沒有形成局部聚集的液滴群，強(qiáng)迫擾動(dòng)的噴霧場(chǎng)與自然噴霧場(chǎng)相比并沒有顯著差別，可以認(rèn)為噴嘴此時(shí)工作在準(zhǔn)穩(wěn)態(tài)條件下。只有當(dāng)擾動(dòng)幅值高于某一臨界值之后，霧場(chǎng)才會(huì)發(fā)生顯著改變。由圖33(b)～圖33(e)可以看出，隨著擾動(dòng)幅值的增大，振蕩壓力場(chǎng)對(duì)于霧化的影響加劇，空間局部聚集的液滴群更加明顯，并且液膜的變形更加劇烈。圖34為不同擾動(dòng)幅值條件下，霧場(chǎng)XY截面體積分?jǐn)?shù)分布灰度圖。從圖中可以看出，擾動(dòng)幅值的大小對(duì)于撞擊前射流的狀態(tài)影響較大。擾動(dòng)幅值增大，射流中液體堆積形成的“鼓包”更加明顯，表明霧場(chǎng)的Klystron效應(yīng)更加顯著。并且擾動(dòng)幅值對(duì)于射流撞擊之后的液膜波動(dòng)狀態(tài)、液膜的厚度分布以及破碎長(zhǎng)度都有影響，液膜的破碎機(jī)制也發(fā)生變化，破碎過程由撞擊波主導(dǎo)開始轉(zhuǎn)變?yōu)樽矒舨ㄅcKlystron效應(yīng)共同主導(dǎo)。

圖33 不同擾動(dòng)幅值下的強(qiáng)迫擾動(dòng)噴霧場(chǎng)Fig.33 Forced atomization field with different perturbation amplitudes

圖35 不同擾動(dòng)幅值下液膜的無量綱破碎長(zhǎng)度Fig.35 Dimensionless breakup lengths of liquid sheet with different perturbation amplitudes

圖36 不同擾動(dòng)幅值下撞擊點(diǎn)下游3d位置的速度分量Fig.36 Velocity component on point 3d downstream impingement point with different perturbation amplitudes

4 結(jié) 論

燃燒不穩(wěn)定的顯著特征是推力室內(nèi)周期性的背壓振蕩，背壓振蕩會(huì)引起噴射速度的周期性變化，從而調(diào)制出與自然霧化顯著不同的霧化特性。本文采用數(shù)值模擬結(jié)合試驗(yàn)的方法研究了強(qiáng)迫擾動(dòng)條件下的射流撞擊霧化特性，主要目的就是全面把握噴嘴的工作特性，進(jìn)一步認(rèn)識(shí)自然霧化與強(qiáng)迫擾動(dòng)霧化之間的差別以及霧化在燃燒不穩(wěn)定中所起的作用，得出的主要結(jié)論如下：

1)基于樹形自適應(yīng)網(wǎng)格與分段線性的VOF方法建立的霧化數(shù)值模擬方案可以細(xì)致刻畫出強(qiáng)迫擾動(dòng)條件下的撞擊霧化過程，模擬得到了強(qiáng)迫擾動(dòng)作用下的周期性噴霧場(chǎng)，可以為燃燒不穩(wěn)定正反饋機(jī)制的研究提供重要的參考，并為非穩(wěn)態(tài)燃燒的數(shù)值模擬提供準(zhǔn)確的初邊值條件。

2)強(qiáng)迫擾動(dòng)的噴霧場(chǎng)與自然噴霧場(chǎng)之間存在顯著差別，自然噴霧場(chǎng)的液滴在空間分布比較均勻，強(qiáng)迫擾動(dòng)條件下的射流撞擊噴霧場(chǎng)出現(xiàn)了弓形液滴群局部聚集的現(xiàn)象，并且在時(shí)間上表現(xiàn)出周期性特征，霧化頻率與強(qiáng)迫擾動(dòng)的頻率一致。在本文研究的頻率范圍(1 257～3 563 Hz)內(nèi)，撞擊式噴嘴的霧化對(duì)擾動(dòng)都有響應(yīng)，表明撞擊式噴嘴的霧化在一個(gè)非常寬的頻帶范圍內(nèi)對(duì)擾動(dòng)都存在響應(yīng)。

3)擾動(dòng)頻率主要影響霧場(chǎng)空間相鄰弓形液滴群之間的間距，以及霧場(chǎng)與振蕩壓力場(chǎng)之間的相位關(guān)系。擾動(dòng)頻率越高，霧場(chǎng)空間某一固定位置與擾動(dòng)壓力之間的相位差越大，并且在相同的距離范圍內(nèi)，相位差的變化區(qū)間范圍也越大。

4)擾動(dòng)幅值決定了霧化Klystron效應(yīng)的強(qiáng)度，只有當(dāng)擾動(dòng)幅值高于某一臨界值時(shí)，霧場(chǎng)才會(huì)表現(xiàn)出周期性特征。擾動(dòng)幅值增大，液膜的破碎長(zhǎng)度減小，撞擊點(diǎn)下游的流量波形由正弦波向陡峭前緣波轉(zhuǎn)變，表明流量由線性變化向非線性變化轉(zhuǎn)變。

5)強(qiáng)迫擾動(dòng)下噴霧場(chǎng)的許多特性都更加有利于激發(fā)燃燒不穩(wěn)定，比如霧場(chǎng)表現(xiàn)出的周期性特征，霧場(chǎng)出現(xiàn)局部聚集的液滴群，以及流量特性由線性向非線性的轉(zhuǎn)變等等，這些特性的出現(xiàn)有可能在正反饋機(jī)制中起到推波助瀾的作用，需要在非穩(wěn)態(tài)霧化的基礎(chǔ)上開展燃燒的研究來進(jìn)一步明確這些特征對(duì)燃燒特性的影響。

下一步工作：采用光學(xué)手段直接測(cè)量噴嘴出口周期性變化噴射速度，為數(shù)值模擬提供更加準(zhǔn)確的邊界條件；進(jìn)一步測(cè)量霧場(chǎng)的微觀信息，采用微距鏡頭結(jié)合圖像處理獲得瞬態(tài)液滴粒徑隨時(shí)間的變化規(guī)律；在霧化研究的基礎(chǔ)上進(jìn)一步建立燃燒響應(yīng)模型，從而構(gòu)建起基于周期性霧化的燃燒不穩(wěn)定分析理論體系。

致 謝

本項(xiàng)工作在國(guó)家超級(jí)計(jì)算天津中心的“天河一號(hào)”超級(jí)計(jì)算機(jī)上完成，感謝天津超算中心的大力支持。