摘要：針對壓燃式發(fā)動(dòng)機(jī)中碳煙演化過程難以被實(shí)驗(yàn)準(zhǔn)確捕捉的問題，研究了高溫高壓環(huán)境中正十二烷噴霧燃燒焰中碳煙瞬態(tài)生成特性及動(dòng)力學(xué)行為，在實(shí)驗(yàn)方面，應(yīng)用消光輻射聯(lián)合技術(shù)以及燃燒成像測速技術(shù)，同步獲得高溫高壓環(huán)境中單孔噴油器正十二烷噴霧燃燒火焰中的碳煙體積分?jǐn)?shù)、碳煙溫度和碳煙速度場的瞬態(tài)分布；數(shù)值計(jì)算方面，基于OpenFOAM環(huán)境，發(fā)展歐拉噴霧模型和非穩(wěn)態(tài)火焰面進(jìn)度變量湍流燃燒模型的耦合模型，噴霧燃燒過程的數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)有高度的一致性。研究結(jié)果表明：相比于參考工況，低噴油壓力工況火焰浮起長度更短，碳煙在當(dāng)量比Φgt;2區(qū)間駐留時(shí)間更長，推動(dòng)了碳煙的生成，碳煙峰值質(zhì)量增加了80%；高溫環(huán)境工況碳煙在當(dāng)量比Φgt;2區(qū)間駐留時(shí)間沒有明顯差異，但更短的火焰浮起長度使得燃料在更高的當(dāng)量比開始燃燒，導(dǎo)致其碳煙峰值質(zhì)量增加了88%；富氧工況的燃燒溫度最高，碳煙粒子運(yùn)動(dòng)速度增大，使得碳煙在當(dāng)量比Φgt;2區(qū)間駐留時(shí)間明顯減少，同時(shí)火焰內(nèi)部當(dāng)量比較小，因而碳煙峰值質(zhì)量減少了42%。

關(guān)鍵詞：正十二烷；消光輻射聯(lián)合技術(shù)；歐拉噴霧模型；碳煙瞬態(tài)生成特性

中圖分類號：TK40 文獻(xiàn)標(biāo)志碼：A

DOI：10.7652/xjtuxb202403002 文章編號：0253-987X（2024）03-0015-10

A Study on Transient In-Flame Soot Evolution of n-Dodecane Spray

XUAN Tiemin1， SHI Zhizhao1， SHANG Weiwei2， HE Zhixia3， WANG Qian1

（1. School of Energy and Power Engineering， Jiangsu University， Zhenjiang， Jiangsu 212013， China; 2. College of

Vehicle and Traffic Engineering， Henan University of Science and Technology， Luoyang， Henan 471000， China;

3. Institute for Energy Research， Jiangsu University， Zhenjiang， Jiangsu 212013， China）

Abstract：Capturing the transient evolution of soot under diesel-like conditions through experimental methods is challenging. The transient formation characteristics and kinetic of soot in-flame of n-dodecane spray under high-temperature and high-pressure conditions is investigated in present study. Experimentally， the combined extinction and radiation （CER） methodology and combustion image velocimetry （CIV） are applied to capture the transient distributions of soot volume fraction， soot temperature， and soot velocity field of n-dodecane spray flames from a single-hole injector under high-temperature and high-pressure conditions. Numerically， a coupled model combining the Eulerian spray model and unsteady flame progress variable combustion model is developed in the OpenFOAM environment. The numerical results of the spray combustion process exhibit a high level of agreement with experimental data. The results show that， compared to the reference condition， in the lower injection pressure condition， the shorter lift-off length （LOL） and the longer soot residence time within the interval equivalence ratio Φgt;2 significantly propel the formation of soot， leading to an 80% increase in peak soot mass. In the higher-temperature condition， there is no significant difference in the soot residence time within the interval Φgt;2， while the shorter LOL leads the fuel to ignite at higher equivalence ratios， resulting in a significant increase in soot formation， with an 88% increase in peak soot mass. In the oxygen-rich condition， the flame temperature is the highest， which increase the particle velocity of soot. Consequently， the residence time of soot within the interval Φgt;2 is significantly reduced. Additionally， owing to the lower equivalence ratio in the flame， peak soot mass is reduced by 42%.

Keywords：n-dodecane; combined extinction and radiation imaging technique; Eulerian spray model; soot formation characteristics

壓燃式發(fā)動(dòng)機(jī)由于其優(yōu)良的經(jīng)濟(jì)性、動(dòng)力性將在相當(dāng)長的一段時(shí)間內(nèi)在商用車、工程機(jī)械、船舶等領(lǐng)域保持主要的動(dòng)力地位。特別是隨著生物柴油、醇類等可再生燃料的應(yīng)用，以及各種新型燃燒模式的大力推進(jìn)，為壓燃式發(fā)動(dòng)機(jī)的發(fā)展創(chuàng)造了新的活力。

壓燃式發(fā)動(dòng)機(jī)碳煙顆粒物的排放對人類健康和大氣環(huán)境都會產(chǎn)生有害影響，因此對燃燒火焰中碳煙形成過程進(jìn)行詳細(xì)深入理解是十分必要的。然而，由于碳煙的生成與氧化，涉及大量復(fù)雜的物理化學(xué)過程，目前還無法實(shí)現(xiàn)對形成過程進(jìn)行精準(zhǔn)的預(yù)測［1］。近些年，大量學(xué)者通過實(shí)驗(yàn)和數(shù)值計(jì)算方法對壓燃式發(fā)動(dòng)機(jī)噴霧火焰中的碳煙形成過程進(jìn)行了大量研究。通過研究發(fā)現(xiàn)，碳煙易形成于高溫低氧環(huán)境［2-4］，可以通過提供噴油壓力［2］、降低環(huán)境溫度［5-6］、減小噴嘴出口尺寸［7］、優(yōu)化噴油策略［8-9］和摻混含氧燃料［10-11］等形式降低火焰中碳煙的生成。目前針對壓燃發(fā)動(dòng)機(jī)噴霧火焰中碳煙的研究多數(shù)集中于準(zhǔn)穩(wěn)態(tài)火焰，對碳煙的瞬態(tài)生成過程，尤其是從碳煙動(dòng)力學(xué)的角度對碳煙生成特性的研究還十分有限。

本文將國際發(fā)動(dòng)機(jī)燃燒合作組織（engine combustion network， ECN）［12］中的A型單孔噴油器作為研究對象，以標(biāo)準(zhǔn)的正十二烷為燃油，通過結(jié)合實(shí)驗(yàn)和數(shù)值計(jì)算的方法，對一款定壓燃燒彈中的噴霧火焰中碳煙的生成特性展開研究。通過消光輻射結(jié)合技術(shù)（CER）和燃燒成像測速技術(shù)（CIV）獲得了碳煙的體積、碳煙溫度場和速度場的瞬態(tài)生成特性，并結(jié)合數(shù)值計(jì)算結(jié)果分析碳煙動(dòng)力學(xué)對碳煙的瞬態(tài)生成過程的影響，為推動(dòng)碳煙模型的發(fā)展提供了一定的理論依據(jù)。

1 實(shí)驗(yàn)裝置及方案

1.1 實(shí)驗(yàn)裝置

本文中所有的實(shí)驗(yàn)在一個(gè)定壓燃燒彈中進(jìn)行。此燃燒彈可以模擬發(fā)動(dòng)機(jī)上止點(diǎn)的熱力學(xué)工況，最高溫度可達(dá)1100K，最大壓力為15MPa。實(shí)驗(yàn)中應(yīng)用的噴油器為ECN中的A型單孔噴油器（210675#），噴孔直徑為90μm。有關(guān)此內(nèi)容詳細(xì)信息可參考文獻(xiàn)［6］。

本實(shí)驗(yàn)中應(yīng)用了高速紋影成像技術(shù)獲得氣相噴霧幾何輪廓的瞬態(tài)發(fā)展過程。另外，通過紋影圖像分析噴霧中灰度值增量變化進(jìn)而獲得了著火延遲期［13］。在相同工況下，應(yīng)用CER技術(shù)同步獲得了碳煙的體積分?jǐn)?shù)和溫度分布。高速紋影和消光輻射結(jié)合技術(shù)的示意圖如圖1所示。此外，在碳煙測試同時(shí)應(yīng)用OH*化學(xué)發(fā)光法獲得了OH*化學(xué)熒光，進(jìn)而獲得了相應(yīng)的火焰浮起長度。有關(guān)光路的詳細(xì)布置和圖片處理方法參見文獻(xiàn)［6］。

通過CER技術(shù)的消光圖像初步獲得的變量為碳煙的光學(xué)厚度KL值，其中K為一個(gè)正比于碳煙的體積分?jǐn)?shù)的參數(shù)，L為光學(xué)厚度路徑。參考文獻(xiàn)［14］中假設(shè)碳煙密度為1.8g/m3，則可以計(jì)算得到每個(gè)像素點(diǎn)光學(xué)路徑累計(jì)的碳煙質(zhì)量，再對整個(gè)碳煙云圖進(jìn)行積分則得到每個(gè)時(shí)刻下總體的碳煙質(zhì)量。為獲得碳煙質(zhì)量和碳煙體積分?jǐn)?shù)應(yīng)用的無量綱消光系數(shù)ke為7.59，ke通過Rayleigh-Debye-Gans （RDG）理論獲得［15］。

為了獲得碳煙生成量和溫度的空間分布信息，本文中還對碳煙的輻射光圖像以碳煙顆粒本身作為示蹤粒子，應(yīng)用類似PIV技術(shù)的圖片處理方式，獲得了碳煙顆粒的速度場，此方法稱為燃燒成像測速技術(shù)（CIV）［16-18］。CIV技術(shù)不需要額外光源即可直接根據(jù)碳煙輻射光強(qiáng)捕捉碳煙位置，由于碳煙有效密度在200～1800kg/m3范圍內(nèi)，碳煙聚合尺寸小于1μm，由斯托克斯阻力定律可得到碳煙形成時(shí)間尺度在毫秒數(shù)量級，因此用CIV計(jì)算速度場時(shí)可以忽略延遲誤差的影響。然而，由于碳煙時(shí)刻處于生成氧化過程，它的輻射光強(qiáng)都會隨著碳煙粒子大小、碳煙體積分?jǐn)?shù)和溫度等發(fā)生變化，此外CIV技術(shù)處理的光強(qiáng)為光學(xué)路徑累計(jì)的結(jié)果，因此通過CIV技術(shù)獲得的碳煙速度場很難獲得準(zhǔn)確的定量結(jié)果。但是，在激光測速技術(shù)難以應(yīng)用的情況下，此技術(shù)還是可以提供有效的半定量信息來分析流場特性。為了盡量滿足PIV查詢區(qū)的1/4原則，同時(shí)得到詳細(xì)流場信息，圖片處理時(shí)選取了雙通道查詢區(qū)，查詢區(qū)大小分別為64像素×64像素和32像素×32像素。

1.2 實(shí)驗(yàn)工況

本實(shí)驗(yàn)中應(yīng)用的燃油為ECN標(biāo)準(zhǔn)的正十二烷。對于每一個(gè)工況點(diǎn)設(shè)定激勵(lì)時(shí)間為3.5ms，其產(chǎn)生的實(shí)際噴霧持續(xù)時(shí)間約為5ms。實(shí)驗(yàn)中，保持相同環(huán)境密度ρg，以工況1為參考工況，其他3個(gè)工況通過變換噴油壓力Pinj、環(huán)境溫度Tg和O2的體積分?jǐn)?shù)構(gòu)成。不同實(shí)驗(yàn)工況參數(shù)信息見表1。在碳煙實(shí)驗(yàn)過程中，每個(gè)工況記錄40次噴油，高速紋影成像實(shí)驗(yàn)每個(gè)工況記錄15次噴油，最后求得樣本平均值以減少循環(huán)波動(dòng)影響。

2 數(shù)值模型及模型驗(yàn)證

2.1 數(shù)值模型

2.1.1 噴霧模型

基于歐拉方法的噴霧模型（Σ-Y）將噴射出具有大雷諾數(shù)和韋伯?dāng)?shù)的液/氣混合物視為具有單一速度場但密度可變的流體來描述射流噴霧［19-20］。對大尺度的流動(dòng)特性及小尺度下的霧化過程分別建模。用平均液體體積分?jǐn)?shù)來表征氣液混合相中的液相，氣液混合相中質(zhì)量加權(quán)平均液相質(zhì)量分?jǐn)?shù)定義為=ρlY/，其中ρl為液相密度；是氣液混合相密度；[AKY-]是液相平均體積分?jǐn)?shù)。假設(shè)氣液混合相中的氣液兩相不互溶，則平均液相質(zhì)量分?jǐn)?shù)與密度的關(guān)系可表示為

式中：μt為湍流黏度；Sc[KG-*3]t是湍流施密特?cái)?shù)。

小尺度的霧化通過求解給定時(shí)間和空間位置每單位體積的相間表面積密度Σ的輸運(yùn)方程來反映，如下式

式中：DΣ是擴(kuò)散系數(shù)；SΣevap是由于燃料蒸發(fā)導(dǎo)致相間表面變化項(xiàng)；CΣ是時(shí)間相關(guān)項(xiàng)；[AKΣ～]eq是驅(qū)動(dòng)局部表面密度的平衡或臨界表面密度；SΣinit是源項(xiàng)。

2.1.2 燃燒模型

在高Damkhler數(shù)（化學(xué)特征時(shí)間比物理特征時(shí)間小）的湍流噴霧燃燒中，局部火焰結(jié)構(gòu)對應(yīng)的射流火焰可用一組射流火焰面來描述湍流火焰［21-22］。應(yīng)用假定擴(kuò)散火焰面方法建立非穩(wěn)態(tài)火焰面/進(jìn)度變量（UFPV）模型。三維空間中擴(kuò)散火焰的任一參數(shù)可以用混合物分?jǐn)?shù)空間中4個(gè)獨(dú)立變量唯一表示，即平均混合分?jǐn)?shù)，平均混合分?jǐn)?shù)方差Z″2，平均標(biāo)量耗散率和進(jìn)度變量Yc。

進(jìn)度變量Yc隨時(shí)間單調(diào)遞增，由燃燒過程中關(guān)鍵組分質(zhì)量分?jǐn)?shù)線性組合得到，如下式

2.2 模型驗(yàn)證

基于實(shí)驗(yàn)中燃燒彈的幾何尺寸、噴油器噴孔尺寸等，以高度為108mm、半徑為25mm的圓柱體中的1/36的三維楔形區(qū)域作為計(jì)算區(qū)域。經(jīng)過網(wǎng)格無關(guān)性分析，將噴孔附近0.09mm內(nèi)劃分為 10個(gè)均勻網(wǎng)格，最小網(wǎng)格尺寸為0.01mm。湍流模型為標(biāo)準(zhǔn)k-ε模型，調(diào)整ε方程參數(shù)Cε1=1.60，以提高對燃油噴射貫穿距的計(jì)算精度，入口邊界上的湍流強(qiáng)度取噴射速度的5％，湍流長度尺度為孔口直徑的10%。

圖2給出了噴霧燃燒全局參數(shù)模型驗(yàn)證結(jié)果。本文遵循ECN對氣相貫穿距和液相長度的定義［12］：氣相貫穿距為從噴嘴出口到燃料混合分?jǐn)?shù)為0.001的最大軸向距離，其中混合分?jǐn)?shù)為燃料的質(zhì)量和混合物總質(zhì)量之比；液相長度為距離噴嘴具有95%液相質(zhì)量分?jǐn)?shù)的最大軸向距離。實(shí)驗(yàn)數(shù)據(jù)分別為15次紋影圖像獲得的氣相貫穿距與40次DBI獲得的液相貫穿距的平均值。從圖2中可以看出，在噴霧后期數(shù)值計(jì)算的氣相貫穿距比實(shí)驗(yàn)結(jié)果稍微偏大，而噴霧液相長度略小于實(shí)驗(yàn)結(jié)果，總體來說數(shù)值計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合很好。

數(shù)值計(jì)算的著火延遲期定義為從燃油噴射開始到出現(xiàn)最大溫度梯度（dTmax/dt）的時(shí)間間隔；數(shù)值計(jì)算的火焰浮起長度從開始噴嘴出口到OH的質(zhì)量分?jǐn)?shù)在該計(jì)算域中達(dá)到其最大值的14％的最小軸向距離，兩個(gè)定義都參照了ECN中關(guān)于著火延遲期和火焰浮起長度的標(biāo)準(zhǔn)。從圖2中可以看出，只有工況3時(shí)，CFD計(jì)算的著火延遲期與實(shí)驗(yàn)值有較為明顯的差別，其他工況下CFD獲得的著火延遲期和火焰浮起長度均與實(shí)驗(yàn)結(jié)果吻合較好。整體表明Σ-Y/UFPV模型能夠準(zhǔn)確地預(yù)測氣液相變和噴霧發(fā)展過程以及著火燃燒過程。

圖3顯示了由數(shù)值計(jì)算得到的噴霧軸線上的速度場和CIV技術(shù)得到碳煙的速度場的對比。從圖3中可以看出，實(shí)驗(yàn)值顯著小于數(shù)值計(jì)算值，這主要是由于碳煙在整個(gè)火焰內(nèi)并非均勻分布，具有較大的隨機(jī)性，通過循環(huán)平均后數(shù)值會大大降低。此外，CIV光學(xué)路徑累積結(jié)果的影響也會帶來較大誤差。在數(shù)值計(jì)算中各工況對噴霧軸線速度場的影響在趨勢上與實(shí)驗(yàn)結(jié)果吻合：工況4時(shí)速度值最高，工況1、3時(shí)速度值比較接近，工況2時(shí)速度值最小。因此，本文CIV計(jì)算得到的速度場可以用于后續(xù)的定性分析。

圖4所示為工況1時(shí)在噴霧發(fā)展的3個(gè)不同階段噴霧對稱面上碳煙體積分?jǐn)?shù)和溫度分布以及碳煙速度場分布，其中黑色曲線表征由紋影得到的噴霧輪廓。從圖4中可以看出，1ms時(shí)碳煙初始生成區(qū)域位于噴霧頭部，與噴霧一起受到空氣阻力作用，初始碳煙速度值較小，且在噴霧最前端由于空氣阻力耦合渦團(tuán)效應(yīng)使碳煙存在較大的徑向速度分量。隨著噴霧的發(fā)展，從2ms到4ms ，碳煙顆粒在噴霧上游區(qū)域高溫的噴霧內(nèi)部運(yùn)動(dòng)，速度值顯著增加，且在向下游軸線方向和噴霧邊緣徑向方向呈遞減趨勢，這與噴霧本身速度場分布趨勢一致。整體而言，高碳煙區(qū)域位于噴霧軸線富油燃燒區(qū)域，高溫區(qū)域靠近碳煙云邊緣，也是化學(xué)計(jì)量燃燒表面區(qū)域。

分別對圖4中碳煙云圖進(jìn)行積分轉(zhuǎn)換和求平均處理，得到火焰中碳煙的總體質(zhì)量、平均溫度T和平均速度值隨時(shí)間演變的過程，如圖5所示，從3個(gè)曲線都可以看出，碳煙分布經(jīng)歷了一些瞬態(tài)變化最后達(dá)到了準(zhǔn)穩(wěn)態(tài)。相對來說，工況3、4時(shí)碳煙達(dá)到穩(wěn)定狀態(tài)時(shí)間更早?？傮w而言，碳煙平均溫度和平均速度達(dá)到穩(wěn)定狀態(tài)略晚于碳煙質(zhì)量，如工況1時(shí)碳煙質(zhì)量約在2.5ms 達(dá)到穩(wěn)定狀態(tài)，而平均溫度和速度則約在 3ms 達(dá)到穩(wěn)定狀態(tài)。碳煙總體質(zhì)量都是先經(jīng)歷了一個(gè)峰值后回落到平穩(wěn)狀態(tài)，不同變量下峰值與穩(wěn)定值存在較大差別。穩(wěn)定階段的碳煙質(zhì)量隨著環(huán)境溫度升高和噴油壓力的降低而升高、隨氧含量的增加而降低。與絕熱火焰溫度趨勢一致，碳煙的平均溫度隨著環(huán)境溫度和氧含量的增加而增加，且氧含量變化對碳煙溫度產(chǎn)生顯著影響。理論上，噴油壓力不改變非燃燒工況下的油氣混合程度，絕熱火焰溫度在兩個(gè)噴油壓力下應(yīng)該保持一致。然而，較高的噴油壓力導(dǎo)致了更快的火焰發(fā)展，更早達(dá)到了穩(wěn)定狀態(tài)，使得在相同時(shí)刻下整體碳煙云的平均溫度略高于工況2的溫度。此外，由圖5（c）可以看出，碳煙平均速度與平均溫度趨勢一致，都隨噴油壓力、環(huán)境溫度和環(huán)境氧含量的增加而增加。定壓燃燒條件下，在相同的噴油壓力時(shí)，噴霧初始動(dòng)量相同，較高的局部溫度導(dǎo)致了更低的局部密度，因此高溫環(huán)境和氧含量導(dǎo)致了較快的碳煙運(yùn)輸速度。兩個(gè)不同噴油壓力所導(dǎo)致的平均速度不同的主要原因是初始動(dòng)量的不同。

3 結(jié)果與討論

本文接下來選取了碳煙演化過程中的兩個(gè)典型時(shí)刻，對瞬態(tài)火焰中碳煙的生成特性和動(dòng)力學(xué)特性進(jìn)行詳細(xì)研究。兩個(gè)時(shí)刻分別為碳煙初生時(shí)刻（圖5（a）中[LL]豎直虛線）、碳煙峰值時(shí)刻（圖5（a）中星號標(biāo)記位置）。

為了進(jìn)一步研究碳煙在空間特定位置動(dòng)力學(xué)特性，通過CFD計(jì)算結(jié)果應(yīng)用了一種混合跡線分析法［23］，以分析特定時(shí)刻下碳煙在運(yùn)動(dòng)軌跡流線上的熱力學(xué)特性和運(yùn)動(dòng)特性?；旌戏?jǐn)?shù)云圖混合跡線如圖6中灰色曲線所示。這些跡線通過對如下公式積分所得

式中：x和r分別為軸向和徑向位置；u和v分別為速度場軸向和徑向的對流項(xiàng)；udif和vdif分別為擴(kuò)散項(xiàng)。有關(guān)此混合跡線詳細(xì)介紹見文獻(xiàn)［23］。本文將以噴霧貫穿距70%、5%噴霧寬度所在跡線為例進(jìn)行分析。

3.1 碳煙初生時(shí)刻

圖7所示為碳煙初生時(shí)刻碳煙的體積分?jǐn)?shù)實(shí)驗(yàn)值和當(dāng)量比（Φ）數(shù)值計(jì)算值分布，定義初生時(shí)刻為碳煙體積分?jǐn)?shù)首次大于0.03×10－6對應(yīng)的時(shí)間。當(dāng)量比分布中黑色曲線表示的為上文所述的5%噴霧寬度所在跡線，豎直灰色直線表示瞬態(tài)時(shí)刻的火焰浮起長度。從圖7中可以看出，所有工況的碳煙初生位置全部位于噴霧的頭部，并且碳煙集中分布于兩側(cè)溫度較高和當(dāng)量比大于1的富油燃燒區(qū)域。整體來說，碳煙初生區(qū)域的當(dāng)量比梯度相對較大。工況1、2為相對較低的反應(yīng)工況，著火延遲期較長，碳煙初生時(shí)刻相對較晚，噴霧得到較長時(shí)間發(fā)展，噴霧貫穿距較長。工況2時(shí)由于相對較慢的噴霧發(fā)展，使得噴霧頭部當(dāng)量比略高于工況1，加上由于低速導(dǎo)致的相對較長的駐留時(shí)間加速了碳煙生成過程。工況3時(shí)噴霧發(fā)展的更慢，導(dǎo)致噴霧頭部當(dāng)量比更高，加上較短的火焰浮起長度處更濃的富油燃燒，也加速了碳煙的生成。工況4的高氧環(huán)境含量雖然導(dǎo)致了更早地生成碳煙，但是較高的環(huán)境氧含量使得其噴霧頭部的當(dāng)量比顯著降低，碳煙生成量并未出現(xiàn)明顯增高。為了進(jìn)一步研究4種工況在相同油氣混合程度下其他因素對碳煙生成的影響，得到了如圖8所示目標(biāo)跡線上相同當(dāng)量比下的平均溫度分布、乙炔的平均質(zhì)量分?jǐn)?shù)YC2H2以及駐留時(shí)間分布t*。其中，星號所標(biāo)記位置為瞬態(tài)火焰浮起長度處的當(dāng)量比。乙炔為碳煙生成的前驅(qū)物，本文的化學(xué)反應(yīng)機(jī)理中未耦合碳煙的生成機(jī)理，數(shù)值計(jì)算中用乙炔含量間接表征碳煙的生成量?；旌羡E線上特定當(dāng)量比的駐留時(shí)間通過下式積分獲得

式中：l為跡線上的運(yùn)動(dòng)距離。此駐留時(shí)間用于衡量在混合跡線上給定的當(dāng)量比下消耗的時(shí)間，此值越高越有利于碳煙的生成。

由圖7（b）可以看出，所有工況時(shí)噴霧頭部當(dāng)量比梯度較大，說明燃料在1lt;Φlt;3的物理空間緊密分布。從圖8中可以看出，相比較工況1，工況2時(shí)由于火焰浮起長度始于當(dāng)量比較大的上游區(qū)域，工況2的火焰浮起長度處Φ≈5.7，而工況1的Φ≈5，使得跡線上的溫度較高，加上較長的駐留時(shí)間，促進(jìn)了乙炔和碳煙的形成。

對于工況3的高溫環(huán)境，由于更短的火焰浮起長度，跡線溫度在高當(dāng)量比Φgt;3處更高，而其駐留時(shí)間略高于工況1，這主要由于碳煙初生時(shí)刻時(shí)工況1已經(jīng)得到更長時(shí)間的發(fā)展，5%徑向?qū)挾鹊嫩E線在高溫工況更靠近噴霧軸線，使其高當(dāng)量比范圍變大。工況4時(shí)高當(dāng)量比范圍顯著減小使得駐留時(shí)間也顯著減小，雖然在Φlt;3時(shí)的高溫導(dǎo)致了較多的乙炔生成，然后此范圍的絕對空間太小，使其對碳煙生成未產(chǎn)生明顯影響。

3.2 碳煙峰值時(shí)刻

圖9所示為碳煙峰值時(shí)刻碳煙的體積分?jǐn)?shù)實(shí)驗(yàn)值和當(dāng)量比數(shù)值計(jì)算值分布。此時(shí)，化學(xué)計(jì)量混合（Φ=1）表面尚未達(dá)到穩(wěn)定狀態(tài)，火焰頭部依然屬于富油燃燒區(qū)域，只是相對于碳煙初生時(shí)刻當(dāng)量比梯度顯著減小，因此高碳煙量區(qū)域依然位于火焰頭部，即在Φ=2等值線附近，文獻(xiàn)［24］表明碳煙一般在Φgt;2范圍內(nèi)產(chǎn)生。所有工況下的火焰浮起長度相對于碳煙初生時(shí)刻都不同程度的向上游移動(dòng)，基本達(dá)到了穩(wěn)定狀態(tài)。

圖10為碳煙峰值不同工況對混合跡線上特定當(dāng)量比時(shí)駐留時(shí)間、溫度和乙炔生成的影響。由圖10可以看出，在高當(dāng)量比Φgt;4范圍內(nèi)， 駐留時(shí)間和溫度分布與碳煙初生時(shí)刻相比，幾乎沒有發(fā)生變化。隨著噴霧發(fā)展，在2lt;Φlt;3范圍內(nèi)的駐留時(shí)間明顯增加，所有工況時(shí)隨著溫度的升高，乙炔在此當(dāng)量比范圍內(nèi)大量生成，并隨著向跡線下游進(jìn)行，即Φlt;2，乙炔逐漸被消耗，一部分演化成大分子環(huán)烴類和碳煙。

當(dāng)量比Φ=2和Φ=3的等值線呈紡錘體形態(tài)，說明噴霧在2lt;Φlt;3范圍得到充分發(fā)展，尤其對于工況1和2，在此范圍內(nèi)幾乎達(dá)到穩(wěn)定火焰狀態(tài)的分布，由混合控制理論可知，噴油壓力對油氣混合程度不會產(chǎn)生明顯影響［25］，上游由于不同火焰浮起長度導(dǎo)致的傳熱影響逐漸消失，此范圍內(nèi)的溫度分布趨于一致。相比工況1，低噴油壓力引起更長的駐留時(shí)間是碳煙生產(chǎn)量更顯著的主要原因。工況3時(shí)，當(dāng)量比2lt;Φlt;3的范圍也趨于一致，此范圍內(nèi)燃燒溫度略高于工況1，又由于相同的噴油壓力使得此范圍內(nèi)駐留時(shí)間也趨于一致，較高的乙炔和碳煙生成則主要由火焰浮起長度處更高的當(dāng)量比燃燒所致。工況4時(shí)在2lt;Φlt;3范圍的更高溫度和更大范圍產(chǎn)生了更多的乙炔，然而明顯降低的駐留時(shí)間，加速了其向下游低當(dāng)量比火焰面的傳輸過程，使得最后碳煙生成量相對工況1并未出現(xiàn)明顯變化。

4 結(jié) 論

本文以ECN組織A型單孔噴油器的正十二烷噴霧火焰為研究對象。通過CER和CIV技術(shù)獲得碳煙體積分?jǐn)?shù)、碳煙溫度和碳煙速度場的瞬態(tài)變化過程。此外，在開源計(jì)算流體軟件OpenFOAM的RANS環(huán)境下，應(yīng)用耦合歐拉噴霧模型（Σ-Y）和非穩(wěn)態(tài)火焰面進(jìn)度變量燃燒模型（UFPV）對此噴霧燃燒過程進(jìn)行數(shù)值計(jì)算。結(jié)合實(shí)驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果，分析了不同工況下正十二烷噴霧火焰碳煙的瞬態(tài)生成特性及動(dòng)力學(xué)特性，得到的主要結(jié)論總結(jié)如下：

（1）較低的噴油壓力導(dǎo)致較小的初始動(dòng)能，進(jìn)而使得碳煙運(yùn)動(dòng)速度更小；較高環(huán)境溫度下導(dǎo)致略高的燃燒溫度，引起的低密度環(huán)境使得碳煙速度略高于工況1；較大的氧含量工況下，燃燒溫度增高，引起的低密度環(huán)境也使得碳煙貫穿速度增加。

（2）低噴油壓力在碳煙初生時(shí)刻主要由于火焰浮起長度的不同，引起更高濃度富油燃燒，產(chǎn)生更高的火焰溫度，進(jìn)一步促進(jìn)了碳煙生成；在碳煙峰值時(shí)刻，低噴油壓力與工況1的當(dāng)量比分布和溫度分布多數(shù)趨于一致，更長的駐留時(shí)間對碳煙生成起到促進(jìn)作用。

（3）高溫環(huán)境下顯著減小了火焰浮起長度，提高了火焰初始位置的當(dāng)量比，使得火焰根部位置產(chǎn)生大量富油燃燒，同時(shí)產(chǎn)生了較高的火焰溫度，導(dǎo)致了大量碳煙生成。當(dāng)火焰充分發(fā)展后，火焰中下游處高溫環(huán)境工況并未對火焰溫度和油氣混合產(chǎn)生明顯影響，駐留時(shí)間和火焰溫度對碳煙生成影響較小。

（4）高氧含量工況下較短的著火延遲期使得碳煙初生時(shí)刻較早，然而高氧含量工況相對貧氧燃燒火焰溫度顯著提高，駐留時(shí)間較其他工況顯著減小，加上較窄的高當(dāng)量比范圍，最終碳煙生成質(zhì)量小于工況1。

參考文獻(xiàn)：

［1］SKEEN S A， MANIN J， PICKETT L M， et al. A progress review on soot experiments and modeling in the engine combustion network（ECN） ［J］. SAE International Journal of Engines， 2016， 9（2）： 883-898.

［2］XUAN Tiemin， PASTOR J V， GARCA-OLIVER J M， et al. In-flame soot quantification of diesel sprays under sooting/non-sooting critical conditions in an optical engine ［J］. Applied Thermal Engineering， 2019， 149： 1-10.

［3］PICKETT L M， SIEBERS D L. Soot in diesel fuel jets： effects of ambient temperature， ambient density， and injection pressure ［J］. Combustion and Flame， 2004， 138（1/2）： 114-135.

［4］賀鵬飛， 李鐵， 姜浩， 等. 柴油噴霧火焰中碳煙顆粒形貌及納觀結(jié)構(gòu)特性 ［J］. 內(nèi)燃機(jī)學(xué)報(bào)， 2019， 37（3）： 251-256.

HE Pengfei， LI Tie， JIANG Hao， et al. Effect of oxygen concentration on soot morphology and nanostructure in diesel spray flame ［J］. Transactions of CSICE， 2019， 37（3）： 251-256.

［5］ZHANG Ji， JING Wei， ROBERTS W L， et al. Soot temperature and KL factor for biodiesel and diesel spray combustion in a constant volume combustion chamber ［J］. Applied Energy， 2013， 107： 52-65.

［6］XUAN Tiemin， DESANTES J M， PASTOR J V， et al. Soot temperature characterization of spray a flames by combined extinction and radiation methodology ［J］. Combustion and Flame， 2019， 204： 290-303.

［7］張志強(qiáng)， 劉一澤， 鄔斌揚(yáng)， 等. 噴孔直徑對重型柴油機(jī)燃燒和排放的影響 ［J］. 內(nèi)燃機(jī)學(xué)報(bào)， 2022， 40（2）： 97-105.

ZHANG Zhiqiang， LIU Yize， WU Binyang， et al. Effect of nozzle diameter on combustion and emissions of a heavy duty diesel engine ［J］. Transactions of CSICE， 2022， 40（2）： 97-105.

［8］HADADPOUR A， JANGI M， PANG K M， et al. The role of a split injection strategy in the mixture formation and combustion of diesel spray： a large-eddy simulation ［J］. Proceedings of the Combustion Institute， 2019， 37（4）： 4709-4716.

［9］JEON J， PARK S. Effects of pilot injection strategies on the flame temperature and soot distributions in an optical CI engine fueled with biodiesel and conventional diesel ［J］. Applied Energy， 2015， 160： 581-591.

［10］PICKETT L M， SIEBERS D L. Non-sooting， low flame temperature mixing-controlled DI diesel combustion ［C］//SAE 2004 World Congress amp; Exhibition. Warrendale， PA， USA： SAE International， 2004： 2004-01-1399.

［11］瞿磊， 王忠， 李瑞娜， 等. 甲醇/生物柴油燃燒顆粒的結(jié)構(gòu)特征研究 ［J］. 西安交通大學(xué)學(xué)報(bào)， 2016， 50（5）： 81-86.

QU Lei， WANG Zhong， LI Ruina， et al. Structural characteristics of methanol/biodiesel combustion particles ［J］. Journal of Xi’an Jiaotong University， 2016， 50（5）： 81-86.

［12］National Technology and Engineering Solutions of Sandia. Engine combustion network ［DB/OL］. ［2023-04-08］. https：//ecn.sandia.gov/.

［13］DESANTES J M， GARCA-OLIVER J M， GARCA A， et al. Optical study on characteristics of non-reacting and reacting diesel spray with different strategies of split injection ［J］. International Journal of Engine Research， 2019， 20（6）： 606-623.

［14］CHOI M Y， MULHOLLAND G W， HAMINS A， et al. Comparisons of the soot volume fraction using gravimetric and light extinction techniques ［J］. Combustion and Flame， 1995， 102（1/2）： 161-169.

［15］KYL ， FAETH G M. Optical properties of overfire soot in buoyant turbulent diffusion flames at long residence times ［J］. Journal of Heat Transfer， 1994， 116（1）： 152-159.

［16］ZHA Kan， BUSCH S， WAREY A， et al. A study of piston geometry effects on late-stage combustion in a light-duty optical diesel engine using combustion image velocimetry ［J］. SAE International Journal of Engines， 2018， 11（6）： 783-804.

［17］DEMBINSKI H， ANGSTROM H E. Swirl and injection pressure effect on post-oxidation flow pattern evaluated with combustion image velocimetry， CIV， and CFD simulation ［C］//SAE/KSAE 2013 International Powertrains， Fuels amp; Lubricants Meeting.Warrendale， PA， USA： SAE International， 2013： 2013-01-2577.

［18］DEMBINSKI H W. The effects of injection pressure and swirl on in-cylinder flow pattern and combustion in a compression-ignition engine ［J］. International Journal of Engine Research， 2014， 15（4）： 444-459.

［19］VALLET A， BURLUKA AA， BORGHI R. Development of a Eulerian model for the “atomization” of a liquid jet ［J］. Atomization and Sprays， 2001， 11（6）： 619-642.

［20］VALLET A， BORGHI R. Modélisation eulerienne de l’atomisation d’un jet liquid ［J］. Comptes Rendus de l’Académie des Sciences-Series IIB-Mechanics-Physics-Astronomy， 1999， 327（10）： 1015-1020.

［21］WINKLINGER J F. Implementation of a combustion model based on the flamelet concept and its application to turbulent reactive sprays ［D］. Valencia： Universitat Politècnica de València， 2014.[LL]

［22］PREZ SNCHEZ E J. Application of a flamelet-based combustion model to diesel-like reacting sprays ［D］. Valencia： Universitat Politècnica de València， 2018.

［23］DESANTES J M， GARCIA-OLIVER J M， NOVELLA R， et al. A numerical study of the effect of nozzle diameter on diesel combustion ignition and flame stabilization ［J］. International Journal of Engine Research， 2020， 21（1）： 101-121.

［24］PICKETT L M， SIEBERS D L. Non-sooting， low flame temperature mixing-controlled DI diesel combustion ［C］//SAE 2004 World Congress amp; Exhibition. Warrendale， PA， USA： SAE International， 2004： 2004-01-1399.

［25］NABER J D， SIEBERS D L. Effects of gas density and vaporization on penetration and dispersion of diesel sprays ［C］//International Congress amp; Exposition.Warrendale， PA， USA： SAE International， 1996： 960034.

（編輯 武紅江 陶晴）

收稿日期：2023-08-14。

作者簡介：玄鐵民（1986—），男，副教授；何志霞（通信作者），女，教授，博士生導(dǎo)師。

基金項(xiàng)目：國家自然科學(xué)基金資助項(xiàng)目（52276116，51906087）。

網(wǎng)絡(luò)出版時(shí)間：2023-12-06網(wǎng)絡(luò)出版地址：https：∥link.cnki.net/urlid/61.1069.t.20231204.1851.008