封閉空間中不同壓力下火焰過孔板加速機理研究

趙健福, 周 磊, 鐘力嘉, 衛海橋

(天津大學 內燃機燃燒學國家重點實驗室, 天津 300072)

0 引 言

發動機小型強化技術被認為是最有前途的點燃式(Spark Ignition，SI)發動機節能減排技術措施之一[1]。然而，小型強化SI發動機熱負荷的增加，導致在其燃燒過程中更容易發生爆震(Knock)[2]、超級爆震(Super- Knock)[3- 4]等不正常燃燒現象，限制了小型強化SI發動機熱效率進一步提升。目前爆震產生的機理尚不明確，末端氣體自燃理論由于很好地解釋了光學實驗的結果而得到廣泛的認可。B?uerle[5]采用激光誘導熒光(Laser Induced Fluorescence，LIF)的方法檢測到內燃機末端氣體的熱點自燃并提出熱點的自燃分為誘發、起燃和傳播3個階段。Kawahara[6]在氫氣燃料壓縮膨脹機上通過高速攝像方法觀測到爆震過程中的自燃和壓力波，并提出爆震強度與末端氣體數量相關。與此不同，超級爆震是由“熱點”(Hot spot)先于火花點火自燃誘發的異常燃燒現象。Zahdeh[7]與Zaccardi[8]等通過內窺鏡觀測到了燃燒室內先于火花點火形成的“熱點”，并根據形成位置及分布排除了排氣門及火花塞電極等高溫表面點火的可能。潤滑油被認為是形成低速早燃“熱點”的可能物質之一，被稀釋后的潤滑油通過活塞環岸積累并竄入到燃燒室內發生蒸發和氧化，最終在火花點火之前發生自燃[9]，自燃后產生的火焰和壓力突變在傳播過程中誘導末端未燃混合氣再次發生自燃，并產生具有強壓力間斷的爆轟波[10- 11]，引發燃燒峰值壓力高達30MPa、壓力振蕩幅值超過20MPa的超級爆震[12]。綜上所述，無論是常規爆震還是超級爆震，本質特征都是未燃混合氣發生自燃，自燃后的湍流火焰- 壓力波相互作用最終導致了劇烈的壓力振蕩。對湍流火焰- 壓力波相互作用的相關研究是揭示爆震、超級爆震等非正常燃燒現象產生機理的關鍵。

從壓力波的產生機理來看，湍流火焰在傳播過程中不斷加速，壓縮前方未燃混合氣形成壓縮波，最終壓縮波互相疊加形成強壓力間斷面。因此，研究火焰加速產生壓力波的過程是火焰- 壓力波相互作用研究的基礎性前提。近年來,相關學者從實驗、理論和數值模擬等角度出發，開展了大量有關火焰加速過程及其機理的研究，其中光滑管和阻塞管均有所涉及[13- 23]。Shchelkin等[24]首次對光滑管內的火焰加速做出解釋，指出已燃氣體的熱膨脹在未燃混合氣中產生了流動并形成湍流，流動的不均勻又進一步造成火焰表面褶皺進而導致火焰加速。在此基礎上，Zeldovich等[25]進一步展開對光滑管內火焰加速過程的探究，提出火焰鋒面與上游不均勻的速度場相互作用導致火焰鋒面延展，這是導致火焰加速的主要原因，而湍流在火焰加速過程中只起到次要作用。根據Ciccarelli和Dorofeev[13]對火焰加速過程的綜述研究，層流火焰的光滑表面在發展過程中會受到水力學不穩定性(Darrieus- Landau instability)和熱擴散不穩定性(Thermal- Diffusion instability)的作用而產生褶皺，進而導致胞狀火焰的形成，增加火焰表面積。相比于光滑管而言，阻塞管中所設立的障礙物能在更大程度上促進湍流的形成和實現火焰加速[26- 31]。Bychkov等[32]提出由于障礙物的存在，在火焰傳播過程中會產生射流，射流會驅使火焰不斷加速。Petchenko等[33]基于障礙物封閉管道對火焰傳播過程進行了數值模擬，結果表明在強烈的Rayleigh- Taylor不穩定性作用下，聲波/壓力波的波動會使火焰鋒面產生很大程度的褶皺。除了Kelvin- Helmholtz不穩定性對火焰- 壓力波相互作用有影響之外，Richtmyer- Meshkov不穩定性也會使火焰- 壓力波相互作用對火焰鋒面產生擾動并提高能量釋放率[34]。肖華華等[35- 36]發現火焰在傳播過程中會與側壁面作用產生壓力波，火焰與壓力波或聲波的相互作用會導致火焰鋒面周期性振蕩以及形成不同的火焰結構。

上述研究從探究火焰加速機理的角度出發，揭示了火焰加速的本質。在此基礎上，相關學者對火焰加速過程的影響因素進行了探究，指出在不同發展階段火焰加速的影響因素也有所不同。Johansen等[19]指出在阻塞管內火焰發展的初始階段，火焰傳播速度低于反應物聲速時，火焰加速強度與障礙物的阻塞比表現出強烈的相關性。在這一階段，阻塞比的增大會導致湍流更加強烈，火焰表面積增大也會使體積燃燒速率升高，從而促進火焰加速。Ciccarelli等[13]則提出隨著火焰繼續向前發展，火焰加速后期主要受到火焰- 壓力波相互作用的影響，這一相互作用過程引發了Richtmyer- Meshkov不穩定性并伴隨有氣流收縮，造成火焰焰舌速度的波動。Valiev等[37]在研究中指出氣體壓縮作用也會降低火焰加速速率以及火焰焰舌速度峰值，在燃燒的早期階段抑制了指狀火焰的加速過程。他們還指出粘度、湍流以及障礙物幾何結構也是影響火焰加速的重要因素。然而,目前火焰過障礙物加速的機理仍不完善，氣流在火焰加速過程中的作用仍值得深入探討。

本文基于自主設計的定容燃燒彈，采用高速紋影技術捕捉火焰過孔板加速的全過程，探索火焰過孔板加速的機理；通過三維數值模擬方法，研究氣流在火焰加速過程中的作用。此外，通過對實驗和模擬結果對比，討論了初始壓力對火焰過孔板加速的影響。

1 實驗裝置及仿真方法

1.1 光學定容燃燒彈實驗平臺

整個實驗系統包括定容燃燒彈本體、高速紋影攝像系統、燃燒壓力采集系統、溫度控制系統、進排氣系統、高壓點火系統和同步控制系統，如圖1所示。定容燃燒彈的燃燒室為圓柱形的腔體，長度230mm，直徑100mm，容積2.32 L，彈體可承受極限壓力為10MPa。在燃燒室的前后壁面分別安裝直徑80mm的圓形石英玻璃視窗。孔板為3mm厚的鋼板，表面分布有不同孔徑的圓孔。Kistler 6113B型缸壓傳感器安裝于燃燒室上壁面，缸壓采集頻率為70 000Hz。Z字型布置的高速紋影系統用于捕捉火焰和沖擊波，高速相機的采樣頻率最高為210 000幀/s。但是，相機的采樣頻率越高，采樣區域越小，為捕捉火焰過孔板加速的全過程，本次研究采用的采樣頻率為90 000幀/s。定容燃燒彈的上下壁面分別安裝有1000W的加熱板，使整個彈體加熱到373K，防止燃燒產物中的水蒸氣冷凝成液滴，影響相機拍攝。

圖1 實現火焰- 沖擊波相互作用的可視化定容燃燒彈

Fig.1Schematicdiagramoftheexperimentalsetupfortheobservationofflameacceleration

本研究選用化學當量比的氫氣- 空氣混合氣作為測試燃料，并通過道爾頓分壓定律確定氫氣和空氣的配比。實驗開始前先將燃燒室抽真空，之后依次充入氫氣和空氣，待靜止2min后，同步觸發點火、圖像和缸壓采集。實驗條件如表1所示。

1.2 CFD仿真方法與驗證

使用美國Converge Science 公司開發的計算流體力學軟件Converge進行定容燃燒彈三維模型計算，如圖2所示。本模型使用的網格初始邊長為4mm。模擬過程中，點火區域在0～1ms的過程中經過了4層16倍的固定加密；點火階段，火焰前鋒的網格邊長為0.125mm。在孔板及其左右3mm, 也進行了4層16倍的固定加密。 此外，在整個火焰傳播過程中，根據當地溫度和H2、OH濃度對當地網格進行了最大4層16倍的動態自適應加密(Adaptive Mesh Refinement，AMR)。此舉可以保證火焰在層流傳播、射流傳播以及湍流傳播3個階段的鋒前網格尺寸為0.250mm。最后的結果顯示，在燃燒過程中，整個模型網格的數量在600～1000萬之間,時間步長在1×10- 8～1×10- 6s之間。

表1 實驗條件Table 1 Test conditions

圖2 定容燃燒彈幾何模型

考慮到實際所用定容彈的主體結構為圓柱體，為減小模擬復雜度，使用三維建模軟件Pro/E所搭建的圓柱模型近似模擬實驗所用的定容彈彈體(如圖2所示)，其總長230mm，直徑100mm。孔板的模擬也依據實驗中所用孔板完成，在圖2中用斷面處的小圓柱體表示。對于不同孔徑和不同位置的孔板，表現為模型中斷處小圓柱體的個數及直徑有所不同，到壁面的距離也有所改變。通過在UDF里添加監測點來實現對壓力的監測。模擬中，監測點的位置與實驗臺架壓力傳感器的位置保持一致并緊貼壁面。模擬火焰在定容燃燒彈中加速時，為了更好地模擬壁面散熱，采用了O′Rourke和Amsden模型，該模型能準確地模擬出熱量經壁面傳播散發的物理過程。由于燃燒進行的速度很快，壁面散熱對燃燒造成的影響較小，所以壁面采用恒溫模型。此外，本研究中采用大渦模擬(Large Eddy Simulation，LES)方法，即對大尺度渦團直接進行求解，對小尺度渦團則先建立亞網格模型再求解。在Converge三維數值模擬方法中，有2種亞網格模型，分別是零方程模型和一方程模型。零方程模型對于射流、混合層、邊界層等的流動比較簡單有效，但是其忽略了對流和擴散的影響，而一方程模型則彌補了這一缺點。本研究中選用了由Yoshizawa[38]和Menon[39]提出的一方程粘性模型。

燃燒模型方面，湍流火焰燃燒通過G方程來捕捉，其中層流火焰速度通過Gülder方程預測。本文采用將G方程和詳細化學反應機理模型SAGE結合的方法來模擬火焰傳播。使用G方程判定火焰面位置，并用化學平衡法計算火焰面處的反應，而在火焰內部和外部采用SAGE求解詳細化學反應機理。其中，SAGE采用CVODES方法(求解常微分方程的初始值問題)求解給定的化學反應機理中每一個單元反應的反應速率，并與輸運方程相耦合，可模擬寬泛的燃燒狀態，包括點火、預混燃燒和非層流燃燒。需要注意的是，采用CVODES方法時，模型的求解精度與化學反應機理直接相關。在化學反應動力學機理方面，采用了包含5元素10組分21反應的詳細反應機理。其中5個元素為C、 H、 O、 N、 Ar，10個組分為： H、 H2、 O、 O2、 OH、 H2O、 N2、 HO2、 H2O2、 Ar。

內燃機缸內流體運動數值的模擬由質量、動量和能量守恒方程以及組分輸運方程控制，控制方程如下所示：

(1)

(2)

(3)

(4)

其中，u是速度，ρ是密度，S是源項，p是壓力，σij為粘性應力張量，Ym是組分m的質量分數，D是質量擴散系數，e是比內能，K是導熱率，hm是組分m的焓值，T是溫度。

2 火焰過孔板加速

2.1 火焰過孔板加速過程

圖3(a)為初始壓力0.2MPa、孔板孔徑2mm下火焰過孔板過程的一系列紋影圖片；圖3(b)為孔板孔徑2 mm、初始壓力0.1、0.2和0.3MPa下火焰傳播速度的變化曲線。如圖3(a)所示，本文中火焰傳播速度是基于火焰前鋒位移對時間的微分計算所得。

如圖3所示，依據火焰傳播形態與速度，將封閉空間中火焰過孔板的加速過程分為3個階段：層流火焰階段、射流火焰階段和湍流火焰階段。

(a)

(b)

Fig.3(a)Chronologicalschlierenimagesofflameaccelerationpassingthroughtheperforatedplate; (b)Evolutionoftheflametipvelocity

層流火焰階段：在火花塞跳火之后，火焰呈球狀向前傳播，表面光滑。此后，由于受到水力學不穩定性和熱擴散不穩定性的影響，光滑的球形層流火焰表面出現胞狀結構[7]。在層流火焰階段，火焰傳播速度較低，約為20m/s左右，表現為先上升后下降的趨勢。在火焰發展初期，由于出現胞狀結構，導致火焰鋒面面積增加，火焰傳播速度略微增加。值得注意的是，本次實驗采用化學當量比的氫氣- 空氣混合氣，其Effective Lewis Number約為1，因此在層流火焰發展過程中，水力學不穩定性的作用會大于熱擴散不穩定性。后期火焰傳播速度下降是由于燃燒膨脹流場受到孔板阻礙作用的結果。如圖3(b)所示，對比不同初始壓力下的層流火焰傳播速度，發現在層流燃燒階段，火焰傳播速度與初始壓力呈負相關。

射流燃燒階段：在2.611ms，層流火焰穿過孔板并被孔板分割成數股射流火焰。過孔板后，火焰傳播速度迅速從16m/s上升至100m/s，增加約一個數量級。隨著射流火焰的發展，火焰傳播速度進一步升高，直至在距離火花塞約40mm處達到某一速度峰值，此后火焰傳播速度開始急劇下降。火焰過障礙物加速一直是一個熱點問題，并得到了廣泛研究。根據Bychkov的理論，在火焰過孔板加速的過程中，膨脹流場過孔板產生的射流驅動了火焰的加速。在過孔板前，膨脹流場沿火焰鋒面法向向四周傳播，當氣流穿過孔板之后，其傳播方向匯聚為一(與燃燒室軸線方向平行)，進而形成一股速度極高的強射流。火焰受到該強射流的驅動，傳播速度急劇增加。穿過孔板后，層流火焰被孔板分割成多束細小的火焰，使得火焰表面積增加，化學反應速率加快；此外，當火焰在障礙物空間傳播時，Rayleigh- Taylor不穩定性和Kelvin- Helmholtz不穩定性誘導了湍流的產生，使得火焰表面褶皺更加劇烈，表面分裂為更小的胞狀結構，進一步促進燃燒速率的加快[20, 40]。

湍流火焰階段：在2.767ms，多股射流火焰匯聚成一束湍流火焰，并在自加速機理的影響下，火焰傳播速度逐漸升高。如圖3(b)所示，過孔板之后，火焰傳播速度隨著初始壓力的增加而升高，與Chaudhuri等[41]對預混火焰爆炸的研究結果一致。

2.2 火焰過孔板加速過程中的流場分析

本研究中高速紋影攝像技術只能識別并拍攝密度梯度較大的射流現象，而無法拍攝到整個燃燒室內的其他流動。因此，為了獲得更多關于氣流速度的信息，通過三維數值模擬對火焰過孔板過程中火焰和氣流的傳播展開探究。

圖4(a)和(b)分別為溫度和速度的分布云圖；圖4(c)所示則為燃燒室中心線上溫度及速度的變化曲線。以火花塞為起點延長120mm作為研究范圍，依據火焰過孔板加速過程選擇5個時刻展開分析，分別為1.551、1.851、2.300、2.475和2.625ms。圖4(a)中，依據溫度變化可以看出紅色部分即為火焰。如圖4(b)所示，由于孔板的節流作用，可以明顯觀察到孔板后產生了強烈的射流現象。如圖4(c)所示，在1.551和1.851ms時刻，從溫度變化曲線可以看出火焰尚未穿過孔板，而在孔板附近可以觀察到明顯的氣流速度峰值，表明此時射流已經產生，并驅使火焰過孔板加速；在2.300、2.475和2.625ms時刻，火焰穿過孔板且火焰鋒面超過射流，此時可以觀察到火焰鋒前的氣流速度隨著與火焰鋒面距離的增加而降低，表明該階段下火焰不再受射流的驅動，相反，未燃區的氣體流動是在火焰鋒面的推動下而運動的。

(c)

圖4 初始壓力0.3MPa下火焰過孔板加速過程中的(a) 溫度分布；(b) 流場分布; (c) 速度分布

Fig.4Distributionof(a)temperature, (b)flowand(c)velocitywhenflamepassesthroughtheperforatedplateunderinitialpressureof0.3MPa

火焰在燃燒室中的加速傳播是不同因素綜合影響下的結果。強烈的射流驅使火焰過孔板加速，而當射流火焰超過射流時，氣體動力學對火焰傳播的影響則不再占據主導地位，此時湍流燃燒理論對火焰傳播起主要作用。湍流火焰傳播會受到火焰- 聲波/壓力波相互作用、Richtmyer- Meshkov不穩定性和湍流自加速的影響。根據湍流燃燒理論，當火焰與聲波/壓力波發生相互作用時，在Richtmyer- Meshkov不穩定性的作用下火焰表面積增加，促進了湍流火焰傳播。

按照圖4的處理方法，繪制初始壓力0.1和0.5MPa下燃燒室軸線上的溫度和速度變化情況，如圖5所示。分析0.1和0.5MPa下的速度分布情況，可以得到與0.3MPa下相同的結論。在火焰未穿過孔板之前，可以觀察到孔板附近有較高速度的強射流存在；當火焰接近孔板時，在該強射流的驅動下，火焰傳播速度急劇增加。而當火焰穿過孔板之后，可以觀察到未燃氣體區域的流場是在火焰的推動下而運動的，射流對火焰發展不再起作用。縱向對比0.1、0.3和0.5MPa的情況，可以發現燃燒溫度隨著初始壓力的升高而略有升高，同樣的，過孔板之后火焰鋒面的傳播速度也隨著初始壓力的升高而增加。有意思的是，在0.1MPa下發現，當火焰過孔板之后，已燃氣體區域的速度遠高于火焰鋒面處的速度，而在0.3和0.5MPa下，已燃區的速度在靠近鋒面時已經開始衰減，火焰鋒面處的速度是整個燃燒室中的最高速度。

(a) 0.1MPa

(b) 0.5MPa

Fig.5Distributionoftemperatureandvelocityunderinitialpressureof0.1and0.5MPa

3 不同初始壓力下火焰過孔板加速實驗與數值模擬結果對比

圖6所示為初始壓力分別為0.1、0.3和0.5MPa下封閉空間中火焰加速過程的實驗結果和數值模擬結果，其中彩色圖片代表實驗結果，黑白圖片表示數值模擬結果。如圖6所示，第一排是過孔板前層流階段的火焰傳播對比圖，第二排是過孔板后射流和湍流階段的火焰傳播對比圖。需要指出的是，圖6中所標注的時間為分別將層流和湍流階段的第一張圖片作為起始時刻下的相對時間。將實驗與數值模擬的火焰傳播圖片截取為相同的尺寸并保證孔板位置對齊，確保火焰鋒面位置進行對比時的準確性。

(a) 0.1MPa

(b) 0.3MPa

(c) 0.5MPa

圖6 初始壓力0.1、0.3和0.5MPa下火焰過孔板加速的實驗和數值模擬結果 (彩色圖片為實驗結果，黑白圖片為數值模擬結果)

Fig.6Comparisonofcombustionevolutionbetweenexperimentandsimulationunderinitialpressureof0.1,0.3and0.5MPa

對比實驗與數值模擬的結果，可以發現隨著初始壓力由低到高，火焰在過孔板前后的發展趨勢、結構形態和傳播速度上均呈現出高度的一致性。就火焰的發展趨勢而言，數值模擬結果與實驗結果一致，表現為層流火焰、射流火焰和湍流火焰3個發展階段。就結構形態而言，過孔板前的層流火焰階段，火焰表面呈現出球型結構，并在靠近孔板時受到射流的影響在火焰鋒面最前方形成了一種近似毛細現象的突起。觀察實驗結果發現，在層流燃燒階段火焰表面出現胞狀結構，這是火焰受水力學不穩定性和熱擴散不穩定性的作用而產生的[7]。而在模擬圖片中則未觀察到該現象，這是因為目前使用的LES方法無法對如此小尺度的結構做到精確的模擬，而繼續細化網格在LES中對于模擬胞狀結構而言收益較小。如果需要模擬出這種微觀結構，直接數值模擬DNS是一種更好的選擇。射流和湍流火焰階段，球形層流火焰鋒面最前方首先穿過孔板形成1束射流火焰，隨著火焰向前發展，考慮到不同初始壓力下孔板孔徑均為5mm，孔板上孔數及孔的分布是固定的，因此可以看到不同初始壓力下均形成3束射流火焰并最終發展合并為湍流火焰。注意到在湍流火焰階段，數值模擬圖片中可以明顯觀察到火焰表面出現褶皺，盡管與實驗結果相比褶皺程度較低。就火焰傳播速度而言，可以發現無論是層流還是湍流階段，在經過相同或相近的時間間隔之后實驗和數值模擬圖片中的火焰鋒面位置仍能保持接近，誤差不到5%，考慮到數值模擬時許多參數設置的局限性，可認為此誤差在正常范圍內。這說明本次模擬的火焰傳播速度與實驗結果呈現出較高的一致性。

圖7所示為初始壓力分別為0.1、0.3和0.5MPa下封閉空間中火焰加速過程的實驗與數值模擬的速度曲線，其中黑色曲線表示實驗得到的火焰傳播速度，紅色曲線表示數值模擬的火焰傳播速度。可以看出，在不同的初始壓力下，實驗與數值模擬的火焰傳播速度曲線在發展趨勢和數值上都呈現出較高的一致性。

在層流火焰階段，無論是哪一初始壓力，實驗和數值模擬的火焰傳播速度在數值和發展趨勢上都極為接近，前期均表現為上升趨勢，靠近孔板時受到孔板的阻礙作用開始下降。初始壓力為0.3MPa下，火焰過孔板后表現出“V”字型的變化趨勢，初始壓力為0.1和0.5MPa下的火焰過孔板后則表現出“S”型變化趨勢。在火焰速度達到峰值之后，由于射流作用減弱，火焰傳播速度下降；其后射流火焰發展為湍流火焰后由于湍流火焰的自加速作用，火焰傳播速度再次表現出上升趨勢。最后由于火焰和壓力波/壓縮波的相互作用傳播速度開始下降。

(a) 0.1MPa

(b) 0.3MPa

(c) 0.5MPa

圖7 初始壓力0.1、0.3和0.5MPa 下火焰傳播速度的實驗和數值模擬結果

Fig.7Comparisonofflametipvelocitybetweenexperimentandsimulationunderinitialpressureof0.1,0.3,0.5MPa

對不同的初始壓力橫向比較可以發現，射流階段的火焰傳播速度差異較小，火焰穿過孔板后其傳播速度從15m/s左右上升至60m/s左右；湍流階段的火焰傳播速度則表現出明顯的差異：0.1MPa下湍流火焰傳播速度峰值約為70m/s，0.3MPa下湍流火焰傳播速度峰值約為90m/s，0.5MPa下湍流火焰傳播速度峰值約為110m/s。

圖8所示為初始壓力分別為0.1、0.3和0.5MPa下封閉空間中火焰加速過程的缸壓曲線，其中黑色曲線表示實驗的壓力數據，紅色曲線表示數值模擬的壓力數據。需要指出的是，由于工作站計算能力有限，計算至壓力峰值出現的時刻所需時間過長，因此本次模擬僅給出部分壓力曲線。從圖中可以看出，壓力開始上升的時刻和壓力上升的趨勢在實驗和數值模擬之間表現出較高的相似性。初始壓力為0.1MPa下實驗和數值模擬壓力開始上升的時刻約為2.72ms；初始壓力為0.3MPa下實驗和數值模擬壓力開始上升的時刻約為2.56ms；初始壓力為0.5MPa下實驗和數值模擬壓力開始上升的時刻約為2.30ms。可以看出，隨著初始壓力的提高，壓力上升的時刻隨之提前。比較不同初始壓力下的情況，初始壓力為0.1和0.5MPa下實驗與數值模擬的壓力曲線吻合度較高，初始壓力為0.3MPa下實驗與數值模擬的壓力曲線則表現出一定的差異，但總體趨勢仍保持一致。

(a) 0.1MPa

(b) 0.3MPa

(c) 0.5MPa

圖8 初始壓力0.1、0.3和0.5MPa 下的缸內壓力實驗和模擬結果對比

Fig.8Comparisonofpressurebetweenexperimentandsimulationunderinitialpressureof0.1,0.3and0.5MPa

受限于工作站計算能力，模擬計算的缸壓未能表現整個燃燒過程，因此對不同初始壓力下的實驗結果進行對比。圖9為不同初始壓力下，缸內壓力及其對應的高通濾波壓力(帶寬為4kHz)隨時間的變化曲線，其中黑色曲線表示缸內壓力曲線，紅色曲線表示濾波壓力曲線，MAPO為最大壓力振蕩幅值。如圖所示，初始壓力為0.5MPa時對應的缸內壓力峰值為2.45MPa，高于初始壓力為0.3MPa下的1.37MPa和初始壓力為0.1MPa下的0.40MPa的壓力峰值。初始壓力為0.5MPa時對應的缸內最大振蕩幅值為0.07 MPa，同樣高于其他2個初始壓力。不難發現，無論是缸內壓力峰值還是缸內壓力的振蕩幅度，都隨著初始壓力的升高而升高。這是因為湍流火焰傳播速度隨著初始壓力的升高而升高，進而對末端未燃氣體產生強壓縮作用，形成較強的壓力波。

(a) 0.1MPa

(b) 0.3MPa

(c) 0.4MPa

圖9 初始壓力0.1、0.3和0.5MPa下的缸內壓力和壓力振蕩

Fig.9Experimentalresultsofpressureandpressureoscillationunderinitialpressureof0.1,0.3and0.5MPa

4 結 論

基于自主設計的定容燃燒彈和Converge三維數值模擬方法，對封閉空間中火焰過孔板加速機理及影響因素開展了研究，并討論了不同初始壓力對火焰過孔板加速過程的影響，得出以下結論：

(1) 依據火焰傳播形態與速度，將封閉空間中火焰過孔板的加速過程分為3個階段：層流火焰階段、射流火焰階段和湍流火焰階段。在過孔板前的層流火焰階段，火焰傳播速度呈現先上升后下降的發展趨勢：在火焰發展初期，由于受到水力學不穩定性和熱擴散不穩定性的影響，火焰表面誘導出現胞狀結構，導致火焰鋒面面積增加，火焰傳播速度略微增加。后期火焰傳播速度下降是由于燃燒膨脹流場受到孔板阻礙作用的結果。此后，層流火焰穿過孔板并被孔板分割成數股射流火焰，火焰傳播速度急劇增加約一個數量級。一方面，膨脹流場過孔板產生的強射流驅動了火焰的加速；另一方面，火焰在障礙物空間傳播時，Rayleigh- Taylor不穩定性和Kelvin- Helmholtz不穩定性使得火焰表面褶皺更加劇烈，表面分裂為更小的胞狀結構，促進燃燒速率加快。最后，多股射流火焰匯聚成一束湍流火焰，并在自加速機理的影響下火焰傳播速度逐漸升高。

(2) 通過數值模擬獲取了火焰過孔板過程中的流場變化情況，發現在火焰未達到孔板前，孔板附近存在速度較高的流場(強射流);當火焰接近孔板時，在該射流的驅動下，火焰穿過孔板后速度急速上升。然而當穿過孔板之后，火焰鋒面前的未燃氣體區域的流場速度沿著遠離火焰的方向逐漸下降，表明在該階段射流對火焰發展不再起主導作用，相反地，火焰開始驅動前方的未燃氣體運動。

(3) 對不同初始壓力下的實驗和數值模擬結果進行對比發現，火焰傳播速度的模擬結果不論從絕對值還是發展趨勢方面都能很好地吻合實驗結果，驗證了數值模擬結果的準確性。縱向對比不同壓力下的火焰過孔板過程，發現隨著初始壓力的提高，孔板前的層流火焰傳播速度略微降低，而孔板后的湍流火焰傳播速度隨著初始壓力的提高而升高，進而誘導較高強度的壓力波，導致較高的缸內壓力和壓力振蕩。