氨燃料燃燒仿真尤其是這個(gè)編號(hào)“主題083”的案例讓我最近一段時(shí)間既踩了不少坑也攢下了大量一手經(jīng)驗(yàn)。如果你正在做氨、氨氫混合燃料的燃燒模擬或者準(zhǔn)備從傳統(tǒng)碳?xì)淙剂限D(zhuǎn)向零碳燃料仿真這篇文章應(yīng)該能幫你少走很多彎路。氨燃料作為航運(yùn)、發(fā)電領(lǐng)域的熱門替代方案燃燒特性非常特殊常規(guī)仿真套路常常失靈需要我們從機(jī)理選擇、湍流模型到求解設(shè)置全鏈路重新思考。我下面就把這套“主題083-氨燃料燃燒”項(xiàng)目從零到一的完整過(guò)程拆開(kāi)講包括方案選型、實(shí)操參數(shù)、遭遇的坑和最終排查方法全部基于真實(shí)計(jì)算記錄可以直接參考。1. 氨燃料燃燒仿真的整體設(shè)計(jì)與思路拆解1.1 氨燃料作為零碳燃料的燃燒特性“反轉(zhuǎn)”把燃料從甲烷換成氨原本順手的CFD燃燒仿真流程會(huì)立刻變得別扭。氨是一種零碳燃料分子式NH3儲(chǔ)存和運(yùn)輸基礎(chǔ)設(shè)施相對(duì)成熟熱值雖然不高但作為低碳轉(zhuǎn)型的重要選項(xiàng)已經(jīng)被廣泛討論。問(wèn)題在于氨的燃燒特性和傳統(tǒng)碳?xì)淙剂喜町愄蟪Ｒ?guī)燃燒模型很難直接套用。從參數(shù)上看氨的最低點(diǎn)火能量高層流火焰速度只有甲烷的1/5左右可燃極限窄絕熱火焰溫度比甲烷低約100K。直觀感受就是氨火焰“點(diǎn)不著、燒不穩(wěn)、燒得慢”。在三維燃燒仿真里這些特性會(huì)造成幾個(gè)連鎖反應(yīng)點(diǎn)火階段需要更高的初始能量火焰前鋒厚度更厚湍流與化學(xué)相互作用的時(shí)間尺度變化很大NOx生成路徑也更加復(fù)雜。根本不是簡(jiǎn)單地替換一下燃料物性就能搞定的。我在項(xiàng)目設(shè)計(jì)階段反復(fù)權(quán)衡過(guò)一個(gè)問(wèn)題仿真到底要解決什么當(dāng)時(shí)定下三個(gè)目標(biāo)第一復(fù)現(xiàn)實(shí)驗(yàn)中的火焰形態(tài)和溫度場(chǎng)第二準(zhǔn)確預(yù)測(cè)出口氮氧化物NOx濃度第三為一個(gè)摻氫燃燒改造方案提供趨勢(shì)參考。三個(gè)目標(biāo)共同決定了后續(xù)所有選型的方向偏向任何一方都容易失衡。1.2 為什么拒絕簡(jiǎn)化全局反應(yīng)選擇組分輸運(yùn)和有限速率模型很多做甲烷燃燒仿真的朋友習(xí)慣用兩步或四步簡(jiǎn)化反應(yīng)比如“甲烷-空氣全局反應(yīng)”計(jì)算快工程上夠用。但氨燃料這里行不通。氨燃燒過(guò)程中NOx的生成非常依賴中間體像NH、NH2、HNO、N2O這些全局反應(yīng)完全無(wú)法描述。如果強(qiáng)行用簡(jiǎn)化機(jī)理出口NOx數(shù)值可能差一個(gè)數(shù)量級(jí)溫度場(chǎng)趨勢(shì)也會(huì)失真。于是“主題083”項(xiàng)目采用組分輸運(yùn)模型與有限速率化學(xué)反應(yīng)的耦合方式也就是說(shuō)在計(jì)算流動(dòng)和傳熱的同時(shí)每個(gè)網(wǎng)格內(nèi)的組分輸運(yùn)方程都帶有詳細(xì)化學(xué)反應(yīng)源項(xiàng)。這相當(dāng)于把一個(gè)化學(xué)反應(yīng)器嵌入到每一個(gè)控制體里用詳細(xì)機(jī)理計(jì)算局部反應(yīng)速率。這樣做的優(yōu)點(diǎn)是很明顯的上游的流動(dòng)變化比如回流區(qū)、渦脫落會(huì)即時(shí)影響組分濃度進(jìn)而影響反應(yīng)放熱再反饋給流場(chǎng)能捕捉局部熄火和再燃這類非線性現(xiàn)象。作為對(duì)比我也測(cè)試過(guò)火焰面模型flamelet model?；鹧婷婺Ｐ驮陬A(yù)混程度高、火焰厚度薄的情況下表現(xiàn)不錯(cuò)但氨/空氣火焰的層流火焰速度低火焰面容易被拉伸局部熄火風(fēng)險(xiǎn)高火焰面假設(shè)會(huì)失真。而且氨的Lewis數(shù)偏離1較大擴(kuò)散效應(yīng)特殊火焰面模型難以處理。最終選擇有限速率模型代價(jià)是計(jì)算量顯著上升但為了NOx預(yù)測(cè)的可信度值得。2. 核心細(xì)節(jié)解析與實(shí)操要點(diǎn)2.1 化學(xué)機(jī)理的取舍不能迷信“越大越好”氨燃燒的詳細(xì)反應(yīng)機(jī)理很多有些包含幾百個(gè)組分、數(shù)千個(gè)反應(yīng)有些簡(jiǎn)化到十幾步。我在項(xiàng)目初篩階段把可用機(jī)理大致分了三類小尺寸骨架機(jī)理幾個(gè)組分、中等尺寸還原機(jī)理幾十個(gè)組分、詳細(xì)機(jī)理上百個(gè)組分。三類算下來(lái)計(jì)算時(shí)間差異能到一個(gè)數(shù)量級(jí)而NOx預(yù)測(cè)精度差異同樣顯著。這里我的經(jīng)驗(yàn)是實(shí)際工程仿真千萬(wàn)別一上來(lái)就選最大機(jī)理。一個(gè)含50個(gè)組分、300個(gè)反應(yīng)的機(jī)理在三維網(wǎng)格上跑起來(lái)每步化學(xué)反應(yīng)求解耗時(shí)巨大尤其是網(wǎng)格量到幾百萬(wàn)以后整個(gè)算例可能根本推進(jìn)不動(dòng)。比較務(wù)實(shí)的路徑是先用理想反應(yīng)器模型進(jìn)行預(yù)篩在0D/1D層面把多個(gè)機(jī)理的預(yù)測(cè)結(jié)果和實(shí)驗(yàn)數(shù)據(jù)做對(duì)比看哪個(gè)在目標(biāo)溫度、壓力、當(dāng)量比區(qū)間內(nèi)表現(xiàn)穩(wěn)定再把它帶進(jìn)三維計(jì)算。在“主題083”里我最終選了一個(gè)中等尺寸的公開(kāi)氨氧化機(jī)理主要篩選標(biāo)準(zhǔn)有三個(gè)層流火焰速度在0.8-1.3當(dāng)量比范圍內(nèi)誤差小于15%NOx轉(zhuǎn)折趨勢(shì)與實(shí)驗(yàn)一致點(diǎn)火延遲時(shí)間數(shù)量級(jí)不偏差。實(shí)測(cè)下來(lái)這套機(jī)理在富燃側(cè)和貧燃側(cè)都能保持合理行為后續(xù)三維結(jié)果的出口NOx趨勢(shì)與實(shí)驗(yàn)吻合度也可以接受。注意機(jī)理驗(yàn)證階段不能只看火焰速度。一定要對(duì)比NOx隨當(dāng)量比的變化趨勢(shì)否則很可能選出一個(gè)火焰速度很準(zhǔn)、NOx卻完全跑偏的機(jī)理。2.2 NOx生成路徑怎么解析、怎么開(kāi)啟氨燃料的NOx生成和傳統(tǒng)燃料有個(gè)關(guān)鍵區(qū)別燃料本身含氮所以“燃料型NOx”是主要貢獻(xiàn)者而且其他路徑也不能忽略。實(shí)際計(jì)算時(shí)要關(guān)注三條路徑熱力型NOx高溫下N2與O2反應(yīng)生成NO主要出現(xiàn)在溫度超過(guò)1800K的高溫區(qū)域。氨火焰溫度即便偏低局部熱點(diǎn)仍可能觸發(fā)熱力型路徑。燃料型NOx來(lái)自燃料中氨氮的轉(zhuǎn)化主要是NH3通過(guò)脫氫變成NH2、NH、N基團(tuán)再與O、OH反應(yīng)生成HNO、NO或者通過(guò)N2O中間體再生成NO。這條路徑在氨燃燒中貢獻(xiàn)最大。瞬發(fā)型NOx碳?xì)淙剂现谐Ｒ?jiàn)氨燃料中占比相對(duì)低但在富燃區(qū)域仍有影響機(jī)理中會(huì)自動(dòng)包含。在CFD軟件中NOx處理有兩種方式后處理方式和在線耦合方式。所謂后處理就是先算出溫度場(chǎng)和組分場(chǎng)再單獨(dú)求解NOx輸運(yùn)方程好處是計(jì)算負(fù)擔(dān)小壞處是NOx的生成對(duì)溫度、OH和N2O濃度非常敏感如果溫度場(chǎng)本身就是用簡(jiǎn)化機(jī)理算出來(lái)的后處理結(jié)果自然不準(zhǔn)。我采用在線耦合方式讓NOx相關(guān)組分直接參與化學(xué)反應(yīng)源項(xiàng)計(jì)算這樣局部溫度波動(dòng)和中間體濃度能實(shí)時(shí)反饋到NOx生成上。實(shí)際仿真時(shí)還特別關(guān)注到N2O路徑。N2O不僅是溫室氣體也是低溫燃燒中生成NOx的一個(gè)高效跳板。如果只開(kāi)啟熱力型NO和燃料型NO常常會(huì)漏掉N2O轉(zhuǎn)化帶來(lái)的那部分NO導(dǎo)致排放值被低估。為了說(shuō)清楚這幾種路徑的貢獻(xiàn)我抽取了火焰中幾個(gè)代表性位置對(duì)比各基元反應(yīng)的反應(yīng)速率整理出一張簡(jiǎn)表NOx生成路徑主要發(fā)生條件關(guān)鍵中間體對(duì)氨燃燒的貢獻(xiàn)熱力型溫度大于1800K、高O2濃度N原子中等受高溫區(qū)影響燃料型氨脫氫后的含氮自由基氧化NH2、NH、HNO非常大N2O路徑低溫、富燃、淬熄區(qū)N2O、NH、NO明顯且常被忽略從表中可以看出如果只算熱力型路徑對(duì)氨燃料來(lái)說(shuō)是嚴(yán)重失真。這也是為什么“主題083”把NOx在線耦合作為必選項(xiàng)而不是可選項(xiàng)。3. 實(shí)操過(guò)程與關(guān)鍵環(huán)節(jié)實(shí)現(xiàn)3.1 幾何簡(jiǎn)化與網(wǎng)格劃分的針對(duì)性調(diào)整“主題083”模擬對(duì)象是一個(gè)實(shí)驗(yàn)室尺度的旋流燃燒器結(jié)構(gòu)不算復(fù)雜但細(xì)節(jié)如果原樣建模網(wǎng)格量會(huì)非常夸張而且許多小圓角、螺栓孔對(duì)燃燒結(jié)果幾乎沒(méi)有影響。所以第一步做幾何清理保留燃燒室主體、旋流器葉片區(qū)域和排氣管去掉內(nèi)部支撐件、法蘭倒角及測(cè)點(diǎn)探針。這個(gè)過(guò)程看起來(lái)簡(jiǎn)單實(shí)際很關(guān)鍵旋流器幾何直接決定流場(chǎng)旋轉(zhuǎn)強(qiáng)度不能隨意簡(jiǎn)化而下游測(cè)量段則可以適當(dāng)拉長(zhǎng)因?yàn)榘比紵磻?yīng)慢火焰和NOx的生成會(huì)在下游持續(xù)很長(zhǎng)時(shí)間。網(wǎng)格劃分上我先用六面體核心網(wǎng)格劃分整個(gè)計(jì)算域旋流葉片附近用局部加密捕捉剪切層和回流區(qū)燃燒區(qū)域加密尺度設(shè)定為1mm左右火焰上游至出口方向逐漸放大到5mm。近壁面處理采用增強(qiáng)壁面函數(shù)第一層網(wǎng)格高度約0.2mm保證最終y接近1這套網(wǎng)格的總量約為300萬(wàn)。這里要特別提醒氨火焰層流火焰速度低火焰厚度比甲烷大所以燃燒區(qū)域的網(wǎng)格密度要求并沒(méi)有想象中那么極限但沿軸向方向需要延伸得更長(zhǎng)。否則火焰后燃段落入粗網(wǎng)格區(qū)NOx守恒就會(huì)出現(xiàn)問(wèn)題。我在網(wǎng)格無(wú)關(guān)性驗(yàn)證中做了粗、中、細(xì)三套網(wǎng)格軸向溫度剖面和出口NOx在中等網(wǎng)格和細(xì)網(wǎng)格之間偏差已小于3%于是把中等網(wǎng)格作為正式計(jì)算網(wǎng)格。3.2 邊界條件、湍流模型與反應(yīng)流求解設(shè)置在邊界條件上燃料入口和空氣入口分開(kāi)設(shè)置。燃料為常溫氨氣模擬時(shí)還疊加了一部分氫氣用來(lái)研究摻氫影響入口質(zhì)量流量根據(jù)目標(biāo)當(dāng)量比0.9計(jì)算?？諝馊肟诓捎铆h(huán)形進(jìn)風(fēng)設(shè)置旋流數(shù)0.6速度方向按旋流器出口角給定。入口湍流強(qiáng)度設(shè)為5%水力直徑通過(guò)實(shí)際截面計(jì)算。求解器選擇上壓力基穩(wěn)態(tài)求解器是首選但氨燃燒反應(yīng)慢穩(wěn)態(tài)很難一步到位。我的做法是先用不帶反應(yīng)的等溫流場(chǎng)算收斂獲得穩(wěn)定的旋流流場(chǎng)后再激活反應(yīng)模型并使用偽瞬態(tài)耦合方法逐步推進(jìn)讓火焰有充分時(shí)間發(fā)展。湍流模型選用SST k-omega它在旋流和回流區(qū)的表現(xiàn)優(yōu)于標(biāo)準(zhǔn)k-epsilon而且對(duì)近壁處理更友好?；瘜W(xué)反應(yīng)第一步激活時(shí)不要直接上完整機(jī)理。我采取“分階段加載”策略先只激活主要組分和基礎(chǔ)反應(yīng)跑上千步判斷溫度場(chǎng)有沒(méi)有上升趨勢(shì)確認(rèn)點(diǎn)火成功后再啟用完整機(jī)理和所有NOx路徑。這樣能顯著避免啟動(dòng)階段振蕩發(fā)散。求解器開(kāi)啟雙精度壓力-速度耦合使用SIMPLE算法亞松弛因子控制為壓力0.3、動(dòng)量0.5、能量0.9、組分0.8時(shí)間推進(jìn)時(shí)先給1e-5秒量級(jí)的時(shí)間步穩(wěn)定后逐步放大到1e-4秒。這里面有個(gè)容易被忽略的細(xì)節(jié)氨的比熱容隨溫度變化明顯高溫下與常溫差異很大所以物性設(shè)置不能使用常比熱假設(shè)必須采用分段多項(xiàng)式或NASA多項(xiàng)式擬合。我踩過(guò)一次坑用常比熱算出來(lái)的火焰溫度比實(shí)驗(yàn)整整低了100K后面果斷換成變比熱才對(duì)得上。3.3 后處理溫度場(chǎng)、火焰面與NOx統(tǒng)計(jì)的驗(yàn)證方法計(jì)算收斂后后處理不只是截圖顏色好看更關(guān)鍵的是提取可量化的指標(biāo)。第一軸向溫度分布沿著燃燒室中心線提取和實(shí)驗(yàn)熱電偶測(cè)量結(jié)果對(duì)比判斷火焰駐定位置是否一致。氨火焰一般會(huì)比天然氣火焰更長(zhǎng)、更靠下游如果模擬溫度提前上升說(shuō)明反應(yīng)速率偏高或點(diǎn)火位置有誤。第二用化學(xué)當(dāng)量混合分?jǐn)?shù)等值面來(lái)顯示火焰面位置配合OH質(zhì)量分?jǐn)?shù)觀察反應(yīng)區(qū)厚度。OH在氨火焰中濃度較低需要選擇合適的閾值否則視野里什么都看不到。第三出口NOx按面積加權(quán)平均再換算成ppm和排放指數(shù)與排氣測(cè)量值對(duì)比。我還做了一步額外的數(shù)據(jù)挖掘把軸向各截面的NO、NO2、N2O質(zhì)量分?jǐn)?shù)分別導(dǎo)出觀察NOx在下游的演化。氨燃燒中NO在高溫區(qū)生成后下游低溫區(qū)往往會(huì)通過(guò)含氮中間體部分還原回N2所以出口截面如果離火焰太近測(cè)到的NO會(huì)虛高相反如果離火焰太遠(yuǎn)計(jì)算域太大數(shù)值耗散又會(huì)掩蓋部分NO還原。最終我把統(tǒng)計(jì)截面定在距燃燒器出口約10倍直徑的位置保證了和實(shí)驗(yàn)采樣點(diǎn)的一致性。4. 常見(jiàn)問(wèn)題與排查技巧實(shí)錄4.1 點(diǎn)火失敗或火焰被吹熄三個(gè)常用“藥方”氨燃料仿真中最讓人抓狂的問(wèn)題就是怎么點(diǎn)火都點(diǎn)不著?，F(xiàn)象是溫度場(chǎng)沒(méi)有上升組分殘差一直高迭代幾千步仍然是冷態(tài)。主要原因在于氨的最低點(diǎn)火能量高、反應(yīng)速率低單純依靠入口高溫來(lái)點(diǎn)燃并不現(xiàn)實(shí)。我嘗試過(guò)幾種處理辦法最終有效的是下面三種組合。第一種凍結(jié)流場(chǎng)點(diǎn)火法。先在不激活反應(yīng)的條件下收斂流場(chǎng)然后在一個(gè)小體積區(qū)域內(nèi)設(shè)置高溫區(qū)溫度給到1800K左右把高溫區(qū)域作為點(diǎn)火源同時(shí)激活反應(yīng)。這樣做的前幾百步高溫區(qū)附近會(huì)開(kāi)始放熱火焰逐漸傳播到整個(gè)剪切層。第二種加入能量源項(xiàng)法。在點(diǎn)火區(qū)域加入一個(gè)隨時(shí)間衰減的能量源項(xiàng)相當(dāng)于火花塞能量注入等火焰自行維持后再把源項(xiàng)降到零。第三種逐步放寬法。初始時(shí)把入口速度降低50%等反應(yīng)穩(wěn)定后再慢慢調(diào)高到目標(biāo)值。這個(gè)技巧對(duì)低火焰速度的氨非常有效模擬里“吹熄”就跟你端著蠟燭走太快結(jié)果火苗熄滅一樣。在“主題083”中真正穩(wěn)定下來(lái)的流程是先凍結(jié)流場(chǎng)再加高溫點(diǎn)火區(qū)和短時(shí)能量源項(xiàng)雙管齊下大約迭代300步后燃燒室溫度開(kāi)始穩(wěn)定上升。需要注意的是高溫點(diǎn)火區(qū)的位置必須放在回流區(qū)邊緣那里有高溫產(chǎn)物回流能把著火核心“拉住”如果放在主流區(qū)很可能直接被吹走。4.2 NOx數(shù)值異常偏高偏低怎么定位出口NOx偏差大是氨燃料仿真里最常見(jiàn)的問(wèn)題而且方向和原因都不一樣。先看偏高的情況。通常是因?yàn)闇囟葓?chǎng)偏高比如沒(méi)開(kāi)輻射模型火焰溫度被高估50-80KNOx就對(duì)溫度非常敏感稍有偏差排放可能差出一大截。氨火焰盡管溫度不如甲烷高但一旦超過(guò)1800K熱力型NOx生成速率指數(shù)上升。解決辦法是開(kāi)啟輻射模型離散坐標(biāo)法配合灰氣體加權(quán)平均吸收系數(shù)能有效降低峰值溫度。同時(shí)檢查壁面邊界是絕熱還是等溫實(shí)際燃燒室壁面有冷卻不能一概用絕熱。再看偏低的情況。氨燃燒的NOx主要來(lái)自燃料型路徑如果機(jī)理中沒(méi)有包含含氮中間體的完整路徑或者沒(méi)有在線激活N2O中間體NOx自然偏低。另外入口空氣中如果含有微量濕度或雜質(zhì)會(huì)影響基元反應(yīng)速率不能完全忽略。還有一個(gè)容易漏掉的因素是計(jì)算域下游長(zhǎng)度。NH3燃燒慢NOx在下游被N2O和NH基團(tuán)還原的過(guò)程會(huì)持續(xù)很長(zhǎng)距離出口截面太近統(tǒng)計(jì)值會(huì)偏低出口截面太遠(yuǎn)又可能讓數(shù)值擴(kuò)散過(guò)頭。我最終采用的方法是沿出口軸線下游布置多個(gè)監(jiān)測(cè)面觀察NOx濃度變化曲線是否平穩(wěn)取平臺(tái)區(qū)域作為統(tǒng)計(jì)值。4.3 發(fā)散與不收斂化學(xué)反應(yīng)剛性的處理方法三維詳細(xì)機(jī)理燃燒仿真中發(fā)散幾乎是每個(gè)項(xiàng)目都要經(jīng)歷的事。氨燃料算例尤其容易因?yàn)榉磻?yīng)機(jī)理剛性很強(qiáng)不同基元反應(yīng)的時(shí)間尺度差異極大從納秒到秒都有。反應(yīng)源項(xiàng)更新稍有不慎能量方程和組分方程就會(huì)振蕩最終導(dǎo)致殘差飆升。我總結(jié)了幾條針對(duì)性措施。使用雙精度求解器單精度在化學(xué)反應(yīng)源項(xiàng)計(jì)算中誤差太大殘差很難降到理想水平。降低組分方程的亞松弛因子尤其對(duì)NOx相關(guān)組分先給到0.5以下等溫度場(chǎng)穩(wěn)定后再逐步調(diào)高。如果穩(wěn)態(tài)計(jì)算反復(fù)振蕩果斷切換到非穩(wěn)態(tài)時(shí)間推進(jìn)。氨燃燒火焰發(fā)展時(shí)間長(zhǎng)穩(wěn)態(tài)算法對(duì)初始場(chǎng)依賴太強(qiáng)時(shí)間推進(jìn)反而能更真實(shí)地展現(xiàn)火焰建立過(guò)程。時(shí)間步長(zhǎng)起步設(shè)為1e-5秒每500步嘗試提升1.5倍如果殘差出現(xiàn)反彈立刻退回上一步長(zhǎng)。整個(gè)算例用了大約5500步迭代加上中間多次調(diào)整設(shè)置總計(jì)算時(shí)長(zhǎng)比同等規(guī)模甲烷燃燒仿真多了將近一倍。這也是氨燃料的真實(shí)代價(jià)要有心理準(zhǔn)備。5. 實(shí)用工具選型與流程規(guī)范5.1 工具鏈怎么配才順手關(guān)于具體軟件我不建議迷信任何單一品牌而是以“機(jī)理驗(yàn)證、三維計(jì)算、后處理”三段式工具鏈為主線。機(jī)理驗(yàn)證階段我會(huì)用理想反應(yīng)器工具計(jì)算點(diǎn)火延遲、層流火焰速度和反應(yīng)路徑貢獻(xiàn)三維階段使用主流壓力基CFD軟件但必須確保支持詳細(xì)化學(xué)反應(yīng)機(jī)理導(dǎo)入后處理階段則關(guān)注軸向剖面、等值面和組分散點(diǎn)圖等能力。工具不需要最貴關(guān)鍵是能方便地導(dǎo)入自定義機(jī)理并且支持雙精度求解。“主題083”項(xiàng)目里我專門做了一套腳本化的邊界條件管理把所有工況參數(shù)統(tǒng)一放在一個(gè)文本配置里這樣批量掃描時(shí)不會(huì)弄混。這個(gè)習(xí)慣很推薦后續(xù)回過(guò)頭查數(shù)據(jù)、復(fù)盤異常工況都方便得多。5.2 計(jì)算流程規(guī)范與數(shù)據(jù)記錄燃燒仿真包含大量中間環(huán)節(jié)每組網(wǎng)格對(duì)應(yīng)什么機(jī)理、什么邊界條件如果只靠腦子記一個(gè)項(xiàng)目做完基本就亂了。我的方法是每個(gè)算例獨(dú)立文件夾里面至少存放四類文件網(wǎng)格文件、機(jī)理文件、邊界條件配置、后處理腳本。命名上明確標(biāo)注版本號(hào)比如“mesh_v3”“chem_okafo_v2”。方案調(diào)整后保持原文件只讀新版本另行復(fù)制。這樣做的一個(gè)直接好處是后續(xù)想追蹤某個(gè)NOx異常數(shù)據(jù)是誰(shuí)引入的可以在半小時(shí)內(nèi)定位到具體文件。另外一個(gè)容易被忽略的點(diǎn)是計(jì)算監(jiān)視。不能只看全局殘差我在燃燒室內(nèi)部布置了若干個(gè)溫度監(jiān)視點(diǎn)、速度監(jiān)視點(diǎn)和組分監(jiān)視點(diǎn)。全局殘差降到1e-4不代表火焰已經(jīng)穩(wěn)定只有監(jiān)視點(diǎn)數(shù)值持續(xù)平穩(wěn)才是真正的收斂。特別是氨火焰發(fā)展慢監(jiān)視點(diǎn)溫度可能還在緩慢上升貿(mào)然認(rèn)為收斂是對(duì)后續(xù)結(jié)果不負(fù)責(zé)的。6. 后續(xù)擴(kuò)展與個(gè)人體會(huì)6.1 從單工況到工況掃描如何用仿真數(shù)據(jù)做燃燒優(yōu)化“主題083”基礎(chǔ)工況完成后我立刻做了變工況掃描。主要變化的參數(shù)是當(dāng)量比、入口溫度和摻氫比例。掃描結(jié)果是典型的氨燃燒規(guī)律當(dāng)量比在1.05附近出現(xiàn)NOx峰值往貧燃側(cè)和富燃側(cè)走NOx都有下降趨勢(shì)富燃側(cè)雖然NOx低但氨泄漏量上升燃燒效率下降。摻氫后火焰速度提升著火穩(wěn)定性變好但NOx排放并不是單調(diào)上升而是非線性的存在一個(gè)最佳摻混點(diǎn)。這類結(jié)論僅靠實(shí)驗(yàn)成本會(huì)很高仿真最大的價(jià)值就在于快速把趨勢(shì)找出來(lái)再集中實(shí)驗(yàn)驗(yàn)證幾個(gè)關(guān)鍵點(diǎn)。如果你要做類似擴(kuò)展建議每次變工況時(shí)保持其他所有條件不變并且從已收斂的鄰近工況作為初始場(chǎng)比如從當(dāng)量比0.9做完后再算1.0時(shí)直接以上一個(gè)溫度場(chǎng)作為初始場(chǎng)收斂速度能加快一半。這也是氨燃燒仿真提效的一個(gè)實(shí)用技巧。6.2 親手做過(guò)一遍后我的三個(gè)核心體會(huì)關(guān)于氨燃料燃燒仿真我的感受可以用三句話概括。第一句話氨和甲烷的仿真差距比大多數(shù)人想象的大得多不是換一個(gè)化學(xué)反應(yīng)文件那么簡(jiǎn)單?；鹧嫠俣?、點(diǎn)火延遲、NOx路徑、下游后燃每個(gè)環(huán)節(jié)都在逼著你重新思考設(shè)置。第二句話機(jī)理驗(yàn)證是大樓地基無(wú)論三維算得多炫機(jī)理選錯(cuò)一切都白搭。我強(qiáng)烈建議你至少花兩三天時(shí)間做0D/1D層面的機(jī)理對(duì)比這一筆投入是最劃算的。第三句話氨燃燒仿真很適合用“先冷態(tài)、后反應(yīng)、先簡(jiǎn)化、后詳細(xì)”的漸進(jìn)策略每一步可控遇到問(wèn)題能準(zhǔn)確定位而不是推倒重來(lái)。最后再分享一個(gè)小技巧調(diào)試階段把壁面溫度設(shè)為略高于實(shí)際值比如450K能顯著幫助點(diǎn)火等火焰穩(wěn)定后再調(diào)回真實(shí)冷卻條件算出來(lái)的排放值才靠譜。希望這篇“主題083”的經(jīng)驗(yàn)總結(jié)能讓你在氨燃料仿真這條路上少踩幾個(gè)我踩過(guò)的坑。