智同生，王 順，郭 利

(1.南京工程學(xué)院，江蘇 南京 211167； 2.哈爾濱工業(yè)大學(xué)，黑龍江 哈爾濱 150001)

0 引言

脈動燃燒是指在一定聲學(xué)條件控制下周期性的燃燒過程，脈動燃燒技術(shù)的應(yīng)用是通過脈動燃燒器來實現(xiàn)的。相比于常規(guī)的穩(wěn)態(tài)燃燒技術(shù)，脈動燃燒具有眾多優(yōu)點：低過量空氣系數(shù)下保持高燃燒效率，較高的燃燒強度和傳熱系數(shù)，污染物排放低(尤其是NOx的排放)[1-3]。主流的脈動燃燒器種類分為三種：Schimidt型脈動燃燒器，Helmholtz型脈動燃燒器和Rijike型脈動燃燒器，其中Helmholtz型脈動燃燒器是目前國內(nèi)外應(yīng)用最廣泛的脈動燃燒器。目前已開展的關(guān)于Helmholtz脈動燃燒器的研究大多都是針對脈動燃燒的工程應(yīng)用問題展開的[4]，其中潘四普[5]研究了脈動燃燒器燃燒室外溫度場的溫度變化規(guī)律，為優(yōu)化加熱器內(nèi)部熱源結(jié)構(gòu)設(shè)計提供參考。陳輝[6]對脈動燃燒器混合室進氣特性進行了數(shù)值模擬研究，得出混合室夾角對脈動發(fā)動機的燃燒穩(wěn)定性、尾氣排放有著重要的影響。Xu yanying[7]對帶彎尾管的脈動燃燒器的換熱特定進行了研究，發(fā)現(xiàn)脈動氣流中產(chǎn)生的渦旋形狀和渦旋核心位置脈動燃燒器的對流傳熱特性有較大影響。對于Helmholtz脈動燃燒器整個燃燒過程的建立，尤其是實現(xiàn)穩(wěn)定點火的機理方面尚缺乏足夠的認識。而研究脈動燃燒的著火特性對控制脈動燃燒器的穩(wěn)定運行，擴展脈動燃燒器的運行范圍都有著重要指導(dǎo)意義，因此有必要對Helmholtz脈動燃燒器的的點火過程進行更加深入的研究。

本文基于課題組之前對Helmholtz型無閥自激燃燒器點火過程的實驗研究，利用數(shù)值模擬方法完善了對Helmholtz型無閥自激燃燒器點火過程的機理分析和探討，分析了從點火開始到實現(xiàn)穩(wěn)定脈動燃燒整個過程中燃燒器的溫度，壓力和氣體流動特性。研究結(jié)果對于說明脈動燃燒器的點火機理，從而更好的實現(xiàn)脈動燃燒器的穩(wěn)定點火和合理運行具有很好的指導(dǎo)意義。

1 實驗設(shè)備及實驗現(xiàn)象

1.1 實驗設(shè)備及系統(tǒng)

本課題組前期對Helmholtz型無閥自激燃燒器點火過程進行了實驗研究[8]，為說明本文模擬結(jié)果的可靠性，同時作為后續(xù)模擬結(jié)果的驗證依據(jù)，首先簡單介紹本課題組前期的實驗系統(tǒng)和實驗結(jié)果。實驗系統(tǒng)如圖1所示，實驗實物圖如圖2所示。實驗過程中，丙烷燃氣和空氣分別在氣瓶壓力和空氣泵的作用下進入脈動燃燒器，在預(yù)混室部位充分混合，然后由電火花點燃，在特制的石英玻璃脈動燃燒器中實現(xiàn)穩(wěn)定脈動燃燒。實驗利用高速相機記錄了整個點火過程的火焰狀態(tài)，并利用熱電偶和壓力傳感器等設(shè)備監(jiān)測了部分點位的壓力和溫度變化情況。

圖1 脈動燃燒器點火實驗研究系統(tǒng)圖1-電腦;2-虛擬示波器;3-壓力傳感器;4-點火電極;5-脈動燃燒器;6-熱電偶;7-高速相機;8-丙烷氣瓶;9-控制閥;10-質(zhì)量流量計;11-脈沖點火控制器;12-質(zhì)量流量計;13-控制閥;14-調(diào)節(jié)閥;15-空氣旁路;16-空氣過濾器;17-調(diào)節(jié)閥;18-鼓風(fēng)機

圖2 脈動燃燒器點火實驗研究實物圖

1.2 實驗現(xiàn)象

實驗過程中利用高速相機記錄了整個脈動燃燒器成功點火過程的火焰狀態(tài)，如圖3所示，從圖中可以看出，點火過程前期存在若干次點火爆燃，從燃燒器入口處點火，爆燃火焰向上傳播，然后熄滅，經(jīng)過若干次這樣的過程，燃燒器上端(出口部位)突然形成火焰，并向下傳播，最終不再熄滅，形成了穩(wěn)定的脈動燃燒火焰。

由于整個實驗所得數(shù)據(jù)有限，可監(jiān)測參數(shù)變化有限，因而，本文進一步利用數(shù)值模擬的方法對脈動燃燒器的整個點火過程進行了研究。

圖3 成功點火過程火焰狀態(tài)

2 數(shù)值模擬

2.1 網(wǎng)格劃分及模擬設(shè)置

物理模型的構(gòu)建和網(wǎng)格劃分是FLUENT進行數(shù)值模擬計算的第一步，本文分別利用CAD和Gambit軟件對脈動燃燒器進行了物理模型構(gòu)建和網(wǎng)格劃分。在網(wǎng)格劃分中，由于脈動燃燒器結(jié)構(gòu)方正，采用了正方形網(wǎng)格，網(wǎng)格質(zhì)量良好。經(jīng)過網(wǎng)格無關(guān)性驗證，最終選取網(wǎng)格數(shù)量為215 593。

同時，在模擬過程中選取了燃燒室入口，燃燒室中心，燃燒室出口，尾管中心四個典型位置作為壓力和速度的監(jiān)測分析位置，如圖4所示，分別命名為監(jiān)測點1，監(jiān)測點2，監(jiān)測點3和監(jiān)測點4。

圖4 脈動燃燒器網(wǎng)格劃分

脈動燃燒器中的燃燒過程是具有多組分化學(xué)反應(yīng)的湍流燃燒，模擬過程中采用標準k-epsilon湍流模型模擬湍流運輸，使用組分輸運方程(species transport)模擬燃燒過程，丙烷燃燒化學(xué)機理采用丙烷一步反應(yīng)模型，利用SIMPLE算法對各方程進行求解。為模擬實驗中的的電火花點火過程，在燃燒室入口的位置使用了Spark點火方式，從0 s開始，點火持續(xù)0.001 s，點火能量為1 J。

2.2 初始和邊界條件

基于脈動燃燒器運行的理論依據(jù)和前期實驗結(jié)果，對模擬中的初始和邊界條件進行了如下設(shè)置：

燃氣/空氣入口條件：入口采用速度進口邊界條件，給入的是丙烷和空氣混合氣，混合比例C3H8/O2為1∶5，進口流速為7 m/s。預(yù)混氣體的初始溫度為300 K壁面邊界條件：根據(jù)之前實驗所測得的燃燒器各處壁面溫度，模擬過程將燃燒器壁面分為三部分，分別是燃燒器混合室壁面，燃燒室壁面和尾管壁面，均設(shè)置為恒溫壁面，溫度分別為400 K，1 000 K和600 K。

尾管出口條件：壓力出口邊界條件。

3 數(shù)值模擬結(jié)果分析

脈動燃燒器點火成功的標志是實現(xiàn)穩(wěn)定持續(xù)的脈動燃燒，而燃燒器中壓力出現(xiàn)穩(wěn)定周期性變化是說明脈動燃燒成功建立的重要標志。因此，本文重點關(guān)注了點火過程中燃燒器中壓力的變化情況，輔以燃燒器內(nèi)溫度場的云圖分布情況來說明成功點火過程。同時為闡明脈動燃燒的點火機理，還進一步分析了煙氣的回流情況和爆燃現(xiàn)象對點火的影響。

3.1 點火過程壓力變化及溫度場的分布情況

整個點火過程中脈動燃燒器中壓力的變化情況如圖5所示，從圖中可以看出，脈動燃燒器點火過程基本可以分為三個階段：第一階段：多次爆燃過程，爆燃過程壓力振幅極大，爆燃頻率較小，稱為爆燃區(qū)；第二階段：由于存在點火延遲，爆燃之后有一個燃燒建立的過程，這個階段壓力幅度較小(大約持續(xù)0.05 s)，稱為過渡區(qū)；第三階段：實現(xiàn)穩(wěn)定脈動燃燒，壓力振幅平穩(wěn)，振蕩頻率穩(wěn)定，稱為穩(wěn)燃區(qū)。根據(jù)點火時燃燒器內(nèi)初始條件的不同，點火過渡過程將向著不同的方向發(fā)展。如果初始條件在燃燒器著火范圍內(nèi)，點火過渡過程將發(fā)展為成功著火；如果初始條件不在燃燒器著火范圍內(nèi)，燃燒器將不會著火，而可能重復(fù)爆燃過程或者轉(zhuǎn)變?yōu)闊o著火。這一結(jié)果和Dawson的研究結(jié)果基本吻合。Dawson[9]研究的是脈動燃燒器啟動過程中的壓力與熱量釋放，得出了脈動燃燒器的啟動時間在100 ms以內(nèi)，啟動過程可以分點火階段、脈動發(fā)展階段、循環(huán)振蕩階段。在脈動發(fā)展階段，由于強烈的非線性作用的影響，頻率和相位都沒有再現(xiàn)性。

圖6是對應(yīng)的點火成功后，穩(wěn)定脈動燃燒一個周期中燃燒器中溫度場的變化情況。從上到下依次為:(a)爆燃后煙氣回流點燃燃燒室內(nèi)新鮮燃料，燃燒室溫度升高;(b)燃燒產(chǎn)生的高溫?zé)煔馀蛎浥懦?(c)新鮮燃料進入，煙氣由于慣性繼續(xù)排出，溫度下降;(d)煙氣回流重新與新鮮燃料混合點燃，溫度升高。整個過程對應(yīng)圖5中穩(wěn)燃區(qū)內(nèi)一個壓力脈動周期。整個溫度場的數(shù)值模擬結(jié)果也和實驗結(jié)果(圖3)基本吻合。

圖5 點火過程的壓力變化情況

圖6 一個脈動周期內(nèi)燃燒器溫度分布

3.2 點火過程爆燃和煙氣回流情況

3.2.1 點火過程爆燃現(xiàn)象

無論是實驗結(jié)果，還是從圖5的壓力變化情況可以得出：脈動燃燒器前期點火過程存在著爆燃現(xiàn)象。爆燃是燃燒室內(nèi)積存的可燃混合物瞬間同時燃燒，從而使燃燒室壓力突然升高的一種現(xiàn)象。為了深入了解分析爆燃情況，進一步查看了圖5中壓力振幅最大區(qū)域的溫度云圖，如圖7(a)所示。最初在燃燒室進口處賦予能量進行點火，由于前期燃燒室燃料和空氣混合物在燃燒室內(nèi)的分布還不夠均勻，容易引起爆燃。爆燃瞬間，整個燃燒器內(nèi)溫度迅速上升，整個燃燒器內(nèi)的可燃氣體被瞬間點燃，產(chǎn)生大量高溫?zé)煔?，進而燃燒室內(nèi)壓力驟升，如圖5及圖7(b)所示。之后大量煙氣向后從尾管急速噴出，造成燃燒室內(nèi)壓力迅速下降，低于外界大氣壓，一些煙氣發(fā)生回流，且有部分外界冷空氣進入燃燒器，出現(xiàn)如圖7(c)所示的現(xiàn)象。

從圖5的圧力變化曲線中可以看出，本模擬中，點火過程脈動燃燒器發(fā)生了兩次爆燃現(xiàn)象，由于每次爆燃時燃燒的燃氣量不同，燃燒放出的熱量不同，進而燃燒室內(nèi)氣體的最大溫升和最高壓力也不同。其中，由于第一次爆燃是混合最不均勻的，其爆燃程度相較與后續(xù)的爆燃會更加劇烈。同時，需要指出的是，在前期實驗過程中發(fā)現(xiàn)，發(fā)生爆燃的次數(shù)并不確定，主要受混合程度的影響。

圖7 爆燃過程燃燒器溫度分布

3.2.2 燃燒器點火過程速度場的變化情況

燃燒器點火過程中，整個燃燒器內(nèi)的速度場分布對燃燒室內(nèi)的燃氣和空氣混合程度有較大的影響。同時，尾管的煙氣回流現(xiàn)象也對燃燒器的成功點火有重要作用。從圖8可以看出，在兩次爆燃過程中。煙氣的回流量和回流深度都很大。

整個燃燒器速度的發(fā)展規(guī)律和壓力基本類似，大致也可以分為3個區(qū)，爆燃區(qū)，過渡區(qū)和穩(wěn)燃區(qū)。爆燃過程燃燒室內(nèi)的氣體流速可以高達數(shù)百米每秒，實現(xiàn)穩(wěn)定脈動燃燒后，氣體流速為數(shù)十米每秒。尾管由于管徑較小，氣體流速要比燃燒室內(nèi)高很多?？偟膩碚f，燃燒器內(nèi)煙氣的速度變化與燃燒器內(nèi)的壓力變化有很大的關(guān)系，壓力絕對值大小影響煙氣的流速大小，壓力的正負影響煙氣的流動方向。壓力振蕩較大時，速度的振蕩也較大。燃燒室壓力由于爆燃急劇升高時，大部分煙氣向尾管方向流動；直至由于大量氣體排出，燃燒室壓力不斷降低低于大氣壓時，尾管煙氣開始發(fā)生回流，重新進入燃燒室。同時，對于整個脈動燃燒器，各個部位的速度變化規(guī)律基本相同。

圖8 點火過程燃燒器內(nèi)速度變化

脈動燃燒器內(nèi)速度的變化和壓力一樣對脈動燃燒器的成功點火和火焰穩(wěn)定性具有很大的影響?；亓鳠煔饽芗涌煨逻M燃氣和空氣的混合，并強化兩者間的換熱，縮短燃燒器時間。Valiev[10]等曾用一維和二維的模擬的方法研究了脈動火焰的穩(wěn)定性，得出結(jié)論：脈動燃燒器內(nèi)的火焰不穩(wěn)定性隨著入口流速和壓力振幅的增加而增加。本文的模擬結(jié)果也表明，煙氣回流量越大，越容易點火成功和后續(xù)新鮮燃氣進入實現(xiàn)穩(wěn)定脈動燃燒。

4 結(jié)論

本文利用FLUENT15.0軟件，在無特殊周期性邊界條件的影響下，采用標準的k-epsilon模型成功模擬了Helmholtz型無閥自激脈動燃燒器的燃燒過程。得到了以下結(jié)論：

(1)數(shù)值模擬結(jié)果與前期實驗結(jié)果基本吻合。說明了標準的k-epsilon模型對于模擬脈動燃燒的可靠性。

(2)燃燒室前期爆燃積攢的能量及高溫?zé)煔獾幕亓魇怯绊慔elmholtz型無閥自激脈動燃燒器的成功點火和穩(wěn)定燃燒的重要因素。