電火花點火式轉子發動機燃燒特性的數值模擬

0醫療模型及火炬的研究進展三角發動機是一種四向的旋轉發動機。與傳統的恢復發動機相比，它具有體積小、質量輕、速度快、功率密度高、運動拆卸少等優點。三角形發動機可在世界范圍內應用于無人機、軍用車輛、小型船只和坦克的輔助設備和空軍交通工具。中國的許多無人機和軍用車輛使用了直徑為1的三角形發動機作為動力和發電裝置。由于轉子發動機燃燒室內可燃物的單向流動,使其更適合燃燒天然氣、氫氣等清潔氣體燃料［2－3］,并且我國天然氣儲量豐富,因此研究天然氣轉子發動機具有重要的理論意義和實際應用價值［4］。國內外學者對天然氣轉子發動機進行了探索性的研究,研究表明為了克服天然氣燃燒速度慢、后燃嚴重等缺點,采用合理的點火位置和點火提前角能在很大程度上改善燃燒過程［5－6］。但是在燃燒過程中,點火提前角與缸內湍流的耦合作用,究竟如何影響火焰的傳播,國內外相關文獻未見報道。本文基于fluent軟件進行二次開發,實現了火花點燃式天然氣轉子發動機的動態模擬,在此基礎上,分析得到了點火提前角對缸內流場乃至缸內燃燒過程的影響。1發動機進排氣道區研究對象是文獻中所用發動機,其由218系列汽油轉子發動機改造而成,發動機進氣方式為周邊進氣,其工況為:節氣門全開,轉速3600r/min,缸內平均指示壓力為0.345MPa.甲烷和空氣混合氣預混后進入氣缸,其燃空比為1∶22.9(當量比為0.75),發動機基本參數如表1所示。如圖1所示,發動機內轉子的運動方向為逆時針運動,轉子轉動1圈,偏心軸轉動3圈,對應旋轉的角度為1080°,其中箭頭處為火花塞的安裝位置。結合轉子發動機的工作特點,進排氣道區是計算區域不變的區域,設置為靜態網格區域;轉子機汽缸和轉子之間的工作區域是計算區域隨時間變化的區域,設置為動態網格區域。由于動態網格的生成需要,采用了非結構的三角形網格。經過反復比較,確定的初始網格數量為39706,圖2為二維初始網格圖。2模型的構建和驗證2.1燃燒模型的建立轉子發動機缸內氣體為可壓縮黏性流動,一般都具有相當高的湍流度,而RNGκ-ε湍流模型嚴格考慮了流線型彎曲、漩渦、旋轉和張力的快速變化,它對于復雜流動有更高精度的預測潛力,適用于轉子發動機內部的復雜流動［8］。因此,湍流模型選擇RNGκ-ε湍流模型。該模型的表達式為式中:ρ為密度;t為時間;v為速度;x為坐標;μeff為有效粘度;Gk為層流速度梯度產生的湍流動能;Gb為由浮力而產生的湍流動能;YM為在可壓縮湍流中過渡擴散產生的波動;C1ε、C2ε和C3ε為經驗常數;ακ和αε分別為κ方程和ε方程的湍流Prandtl數;Sκ和Sε為源項。燃燒采用一步反應方程式,即為了綜合考慮缸內化學反應速率和湍流脈動對燃燒的影響,選擇層流有限速率/渦耗散結合的計算模型,其中Arrhenius和渦耗散反應速率都進行計算,凈反應速率取兩個速率的較小值［9］。進排氣口采用壓力邊界條件,設定進排氣口壓力為大氣壓。發動機缸壁溫度采用壁面邊界條件,溫度根據經驗值設定為400K.氮氧化物生成模型采用“熱力型”和“快速型”生成機理。2.2實驗結果分析首先,對流場計算結果進行了驗證,由于缺乏該機型流場的實驗數據,而DeFilippis等［10］利用激光多普勒測速儀對另外一種同樣采用周邊進氣的轉子發動機進行了測量,得到了缸內流場的變化過程,圖3為文獻實驗中激光多普勒測速儀的測點位置分布。所以利用已建立的動態模擬模型對文獻中發動機的流場進行數值模擬,模擬結果如表2所示。從表2中可以看出,計算得出的缸內流場分布規律與實驗結果十分相似,包括缸內流場的流速以及渦流的位置。比如曲軸轉角730°CA時發動機進氣口附近的流速,實驗結果和計算結果都約為20m/s,而渦流位置,實驗結果和計算結果都位于燃燒室前部區域。曲軸轉角1000°CA時由于燃燒室容積減小,渦流被壓縮而破碎成單向流,實驗結果和計算結果都顯示流場為單向流,并且流速都約為10m/s,由此說明選用的湍流模型對該類周邊進氣轉子發動機缸內流場的計算是準確的,對于其缸內流場具有高度的預測能力。同時為了驗證燃燒過程計算結果的正確性,采用所建立的計算模型,對文獻中由218系列改裝的預混天然氣轉子發動機進行了數值模擬,并與文獻中改裝后發動機的實驗數據進行對比。圖4和圖5分別是計算結果與文獻結果的示功圖對比和燃料已燃質量分數對比,發動機的點火提前角為42°CA,對比的范圍從點火時刻到上止點后60°CA.計算結果與文獻結果的氣缸壓力平均誤差約12%左右,燃料已燃質量分數平均誤差約9%左右。從圖4和圖5中可以看出:計算結果與文獻結果大體趨勢一致,吻合良好,只是在示功圖(見圖4)中壓力峰值處略有差別,說明所建立的計算模型是合理的。3計算與分析3.1點火后流場分布轉子發動機設計的點火提前角為42°CA,現分別取60°CA、56°CA、50°CA、47°CA、42°CA、36°CA作為點火提前角,并分析了相關計算結果。圖6給出了不同點火提前角時燃燒室內的流場,其中黑色圓點表示點火位置,標注的BTDC表示上止點前。由圖6可以看出,發動機缸內的流場可分為3個區域,即燃燒室后部的隨轉子運動方向的單向流區域、燃燒室中部的渦流區域、燃燒室前部的隨轉子運動方向的單向流區域。另外,隨著點火角度的推遲,渦流區域與點火中心的距離不斷縮小,其中60°CA、56°CA時二者距離較大,點火時點火位置完全處于燃燒室前部的隨轉子運動方向的單向流區域;當點火提前角分別為50°CA、47°CA、42°CA時,點火位置處于燃燒室中部渦流區與前部單向流區之間的過渡區域;而36°CA時,二者距離最小,點火中心已基本處于渦流區域中。根據點火時,點火位置所處的3種不同的區域,表3~表5給出了點火提前角分別為60°CA、47°CA、36°CA時發動機做功行程中流場和溫度場的演變過程。表中BTDC表示上止點前,ATDC表示上止點后。缸內的流場是隨著時間不斷變化的,當采用不同點火提前角時,點火時刻不同,火花塞附近氣體的流動速度和方向均不同,點火后,流場的不同還導致反應氣體和產物的傳熱傳質不同,直接影響了到火焰傳播的方向和速度,進而影響到發動機的經濟性、動力性及排放性能。表3~表5分別描述了點火后缸內流場和缸內火焰輪廓分布圖。燃燒初期,當點火提前角為60°CA時,點火位置完全處于燃燒室前部的隨轉子運動方向的單向流區域,火焰傳播主要受到燃燒室前部單向流的影響,如表3不同時刻溫度流場中曲軸轉角10°CA時,火焰中心形成后,火焰前鋒很快隨氣流的方向傳播到燃燒室的前部,引燃前部的未燃混合氣,而很難逆向地向燃燒室后部傳播;當點火提前角為47°CA時,點火位置已處于燃燒室中部渦流區與前部單向流區之間的過渡區域,如表4不同時刻溫度流場中曲軸轉角10°CA時,火焰前鋒除了隨燃燒室前部的單向流,快速地向前傳播外,還會由于燃燒室中部渦流的作用,把已燃的高溫混合氣帶回燃燒室的中后部,從而引燃中后部的未燃混合氣,這在很大程度上改善了后部未燃混合氣不易燃燒的狀況;當點火提前角為36°CA時,點火位置已基本位于中部渦流區域。火焰傳播主要受到中部渦流的作用,從表5可看出不同時刻溫度流場中曲軸轉角10°CA時,雖然渦流區域的流速很大,但只有少部分已燃混合氣會被渦流帶回燃燒室的中后部,這主要有兩個方面的原因:首先,火花塞點火時處于渦流區域中,此時其附近流場的速度很大,初期火焰中心已燃混合氣產生的熱量容易散失,不利于燃燒初期火焰中心的快速增大,所以初期已燃混合氣較少;其次,如表5不同時刻流場中曲軸轉角10°CA時,此時缸內流場中部的渦流已經開始破碎,高速渦流對已燃高溫混合氣的作用時間太短。燃燒中期,當點火提前角為60°CA時,如表3不同時刻流場中曲軸轉角20°CA時,轉子由于接近上止點,導致燃燒室空間狹小,從而燃燒室內部渦流消失,缸內整體流場為隨轉子運動方向的單向流,其進一步推動已燃的高溫混合氣繼續向燃燒室前端運動,如表3不同時刻溫度流場中曲軸轉角20°CA時,導致前部混合氣的快速燃燒,而后部未燃混合氣無法及時燃燒;當點火提前角為47°CA時,如表4不同時刻流場中曲軸轉角20°CA時,同樣由于燃燒室空間的狹小,缸內渦流消失成隨轉子運動方向的單向流,這同樣加速了燃燒室前部的未燃混合氣的燃燒,如表4不同時刻溫度流場中曲軸轉角20°CA時,而燃燒室尾部的未燃混合氣仍舊無法及時的燃燒;當點火提前角為36°CA時,如表5不同時刻流場中曲軸轉角20°CA時,缸內渦流消失形成的隨轉子運動方向的單向流,很快把中部已燃的高溫混合氣推向前部未燃混合氣區域,如表5不同時刻溫度流場中曲軸轉角20°CA時,所以整體上中部和前部的未燃混合氣可以及時燃燒,但后部的混合氣卻無法及時地被引燃。燃燒后期,當點火提前角為60°CA時,由于燃燒室容積增大,轉子凹坑內會出現一個強度不大的渦流,如表3不同時刻流場中曲軸轉角60°CA時,這個渦流可以在很大程度上加速燃燒室后部的燃燒速率。但是如表3不同時刻溫度流場中曲軸轉角60°CA時,由于此時中后部的混合氣沒有被引燃,火焰無法利用渦流對其燃燒速率的加速作用,所以整個燃燒過程持續時間長,燃燒速率較低;當點火提前角為47°CA時,如表4不同時刻流場中曲軸轉角60°CA時,同樣轉子凹坑內會出現一個強度不大的渦流,其加速了前部已燃高溫混合氣與后部未燃混合氣的混合速度,增加了反應速率,此時如表4不同時刻溫度流場中曲軸轉角60°CA時,燃燒室內的未燃混合氣基本燃燒干凈,整個燃燒過程縮短,燃燒效率較高;當點火提前角為36°CA時,如表5不同時刻溫度流場中曲軸轉角60°CA時,同樣由于此時后部的混合氣沒有被引燃,所以火焰無法利用此時流場中渦流對其燃燒速率加速作用,整體上燃燒過程較長,燃燒速率較低。綜上所述,當點火提前角在47°CA附近時,燃燒效率最高。這是因為點火提前角在47°CA時,火花塞處于渦流到單向流的過渡區域,并且其附近流場速度不是太大,這有利于燃燒初期火核的快速形成,燃燒速度的快速增加。在渦流消失時間一定的情況下點火提前角47°CA相比較原點火提前角42°CA,更加充分利用了渦流對火焰傳播所作用的時間。雖然點火提前角60°CA、56°CA、50°CA也充分利用了渦流的作用時間,但由于在這些時刻點火時,點火位置位于燃燒室前部隨轉子運動方向的單向流區域附近,火焰傳播主要受到燃燒室前部單向流的影響。火焰中心形成后,火焰前鋒很快隨氣流的方向傳播到燃燒室的前部,引燃前部的未燃混合氣,而很難逆向地向燃燒室中后部傳播,所以位于中部的渦流沒能起到加速火焰傳播的促進作用。總體上說:點火提前角以47°CA為一個分水嶺,即在渦流消失時刻一定的情況下,為了盡量充分利用渦流的作用時間,應該適當地增大點火提前角,以47°CA為界限,燃燒效率最高,點火提前角增加或降低時,燃燒速率均減小。3.2點火提前角對燃燒過程的影響從圖7不同點火提前角缸內的燃料已燃質量分數曲線可以看出:當點火提前角為47°CA時,燃燒過程中燃料已燃質量分數曲線的斜率最大,燃燒持續時間最短,即未燃混合氣燃燒速度最快,50°CA、42°CA次之,60°CA、56°CA、36°CA較差,所以燃燒速率以點火提前角47°CA時最高,并向兩邊遞減。由示功圖8可以看出:點火提前角為47°CA時,壓力峰值最大,50°CA、42°CA時次之,60°CA、56°CA、36°CA時較小,其中47°CA、50°CA時,比原點火提前角42°CA時大,即壓力峰值的整體趨勢與燃燒效率的變化趨勢相似,都是以點火提前角47°CA時最高,并向兩邊遞減。這主要是因為:當點火提前角為60°CA、56°CA時,點火位置在燃燒室前部隨轉子運動方向的單向流區域,導致燃燒室后部的混合氣不能及時燃燒,整個燃燒過程持續時間太長,所以其燃燒效率和壓力最大值小于其他點火提前角的情況。當點火提前角為50°CA、47°CA、42°CA時,火花塞點火時,其位于湍流到單向流的過渡區,點火后,火焰會同時向前部和后部傳播,缸內混合氣可以及時燃燒,整個燃燒過程持續時間短,所以它們的壓力曲線峰值要明顯高于其他點火提前角。其中當點火提前角為47°CA時點火位置處于湍流到單向流的過渡區,而且相比原點火提前角42°CA,其更充分利用了渦流對火焰傳播加速的作用時間,所以燃燒速率最大,壓力曲線峰值也最大。雖然點火提前角為50°CA時,相比原點火提前角42°CA也充分利用了渦流的作用時間,但其點火位置更接近于前部的單向流區,火花塞點火時,火焰前鋒很快隨氣流的方向傳播到燃燒室的前部,引燃前部的未燃混合氣,而很難逆向地向燃燒室中后部傳播,所以位于中部的渦流沒能起到加速火焰傳播的促進作用。當點火提前角為36°CA時,點火位置點火時,雖然其已完全位于渦流區域中,但點火位置處流場速度過大,不利于初期火焰的形成,并且在燃燒中期,渦流對已燃高溫混合氣的作用時間較短,所以整體燃燒速率仍然較低,同樣無法獲得較高的壓力峰值。由圖9缸內平均溫度曲線和圖10NOx生成曲線可以看出:NOx生成量在點火提前角為50°CA時最多,56°CA、47°CA、42°CA時次之,60°CA、36°CA時較少。點火提前角為47°CA時,燃燒速率最大,但其NOx生成量卻低于50°CA時,這主要是由于其燃料的快速燃燒導致其燃燒持續期相比點火提前角為50°CA時短,所以在發動機的膨脹行程行程初期其缸內平均溫度比點火提前角為50°CA時的缸內平均溫度大,但在