我最近在折騰氨燃料燃燒仿真這個課題前后跑了半年多的算例踩了不少坑也總結(jié)出一些還算靠譜的經(jīng)驗。項目代號就叫“主題083”實際內(nèi)容是一臺小型燃?xì)廨啓C(jī)燃燒室里燒氨氫混合燃料的全流程數(shù)值模擬。氨燃料近幾年在航運、發(fā)電、工業(yè)加熱領(lǐng)域的熱度很高大家都把它當(dāng)成零碳燃料的候選方案之一但真正能把這個燃燒過程仿明白的人并不多。原因很簡單氨的燃燒特性和甲烷、氫氣完全不同層流火焰速度低、自燃溫度高、可燃范圍窄再加上NOx生成路徑特別敏感稍有一點設(shè)置不當(dāng)計算結(jié)果就完全跑偏。這篇內(nèi)容就是圍繞氨燃料燃燒仿真從機(jī)理選擇、模型配置、網(wǎng)格策略到后處理判讀的完整實操記錄適合正在做燃燒室設(shè)計、排放預(yù)測或者實驗對照的工程師參考。1. 氨燃料燃燒仿真的核心挑戰(zhàn)與技術(shù)背景1.1 氨作為零碳燃料的兩面性氨分子里沒有碳燃燒產(chǎn)物理論上只有水和氮氣所以被看作理想的氫載體和直接燃料。但“理論上”三個字是重點實際燒起來遠(yuǎn)沒有這么干凈。氨的熱值低低位發(fā)熱量大約只有甲烷的一半左右也就是說同樣功率輸出需要更大的燃料流量。更麻煩的是它的化學(xué)反應(yīng)活性低常溫常壓下氨/空氣預(yù)混火焰的層流燃燒速度只有每秒十幾厘米的量級甲烷大約每秒四十厘米左右氫氣更是每秒兩三米?；鹧?zhèn)鞑ヂ馕吨鹧嫒菀妆淮迪ā⒒鼗痫L(fēng)險低也意味著在高速流場里穩(wěn)定火焰非常困難。這個特性直接決定了仿真思路。你用常規(guī)天然氣燃燒室的經(jīng)驗去套氨燃料幾乎不可能得到合理結(jié)果?；鹧驿h面位置、溫度分布、污染物生成量都會失真。所以做氨燃燒仿真之前先得把它的燃燒化學(xué)基礎(chǔ)盤清楚。氨的氧化路徑不是簡單的“燃料氧氣二氧化碳水”而是一整套包含氨基、亞氨基、硝酸基、一氧化氮、二氧化氮、一氧化二氮等多種中間產(chǎn)物的連鎖反應(yīng)網(wǎng)絡(luò)。仿真中如果機(jī)理沒選對后面所有內(nèi)容都是白搭。1.2 為什么必須用CFD而不是簡單零維計算氨燃燒的工程問題本質(zhì)上是三維湍流反應(yīng)流問題火焰穩(wěn)定、回火、吹熄、NOx生成這些現(xiàn)象都受到流場結(jié)構(gòu)、摻混程度、停留時間分布的多重影響。零維化學(xué)動力學(xué)計算只能給出理想條件下的著火延遲和產(chǎn)物趨勢無法回答“這個旋流器角度下火焰穩(wěn)不穩(wěn)”“壁面附近會不會形成高溫區(qū)導(dǎo)致熱力型NOx飆升”這類問題。CFD的價值就在于把流動、傳熱、化學(xué)反應(yīng)耦合在一起在給定幾何和工況下給出空間分布信息。我做的主題083項目就是典型例子。燃燒室入口有旋流葉片空氣和氨氫混合氣通過不同的通道進(jìn)入如果不做三維仿真你根本想象不到回流區(qū)中心會形成一個低溫富氨核心而高溫反應(yīng)區(qū)集中在剪切層附近。這種流場結(jié)構(gòu)對NOx生成的影響非常大零維模型完全捕捉不到。所以結(jié)論很明確只要涉及具體燃燒室結(jié)構(gòu)和運行工況評估CFD就是必需的關(guān)鍵是怎么把湍流模型、燃燒模型、化學(xué)機(jī)理以及數(shù)值格式這些環(huán)節(jié)組合好。2. 化學(xué)機(jī)理的選擇——氨燃燒仿真最重要的一道門檻2.1 三類機(jī)理怎么選氨燃燒化學(xué)機(jī)理的選擇是整條仿真鏈路里最影響結(jié)果、也最容易被新手忽略的環(huán)節(jié)。市面上能見到的氨氧化機(jī)理大致分三個檔位詳細(xì)機(jī)理、骨架機(jī)理和簡化機(jī)理。詳細(xì)機(jī)理包含幾十種組分和幾百步反應(yīng)描述精度最高覆蓋高溫點火和低溫氧化路徑但計算開銷極大。骨架機(jī)理是在詳細(xì)機(jī)理基礎(chǔ)上通過敏感性分析和路徑分析裁剪出來的組分?jǐn)?shù)量壓縮到二十種上下反應(yīng)步驟幾十步工程精度基本夠用。簡化機(jī)理進(jìn)一步壓縮有些甚至只有幾步集總反應(yīng)通常用于快速評估或者嵌入系統(tǒng)級仿真。我用過一個包含五十多種組分、四百多步反應(yīng)的詳細(xì)機(jī)理做零維反應(yīng)器驗證精度確實不錯著火延遲時間和實驗數(shù)據(jù)對得上但一旦放到三維網(wǎng)格里跑穩(wěn)態(tài)計算每步迭代光化學(xué)反應(yīng)源項就要消耗大量時間。后來換用某研究團(tuán)隊發(fā)布的骨架機(jī)理只有十九種組分、六十多步反應(yīng)計算速度提升了數(shù)倍關(guān)鍵參數(shù)如絕熱火焰溫度、層流火焰速度的偏差控制在百分之五以內(nèi)在工程判斷維度上完全夠用。機(jī)理類型典型組分?jǐn)?shù)量典型反應(yīng)步數(shù)適用場景計算開銷工程推薦度詳細(xì)機(jī)理50300基礎(chǔ)研究、機(jī)理驗證極高低骨架機(jī)理15-2550-80三維燃燒室仿真、排放預(yù)測中高高簡化機(jī)理5-128-30系統(tǒng)耦合、快速評估低中2.2 實際項目中的簡化策略與驗證方法選機(jī)理不能只看文獻(xiàn)推薦必須做驗證。我的做法是先用零維均質(zhì)反應(yīng)器把候選機(jī)理的著火延遲時間和層流火焰速度算一遍和實驗數(shù)據(jù)做對比。這里有個容易被忽略的坑氨的著火延遲時間在低溫段對壓力極其敏感如果你仿真的燃燒室工作壓力在十個大氣壓以上一定要確認(rèn)機(jī)理是否在高溫高壓條件下做過驗證很多公開發(fā)表的機(jī)理是基于常壓實驗數(shù)據(jù)標(biāo)定的外推到高壓時著火延遲能差一個數(shù)量級。在三維仿真里我強(qiáng)烈建議不要直接從詳細(xì)機(jī)理入手。將骨架機(jī)理導(dǎo)入求解器后先跑一個簡單的二維層流火焰驗證算例確認(rèn)火焰?zhèn)鞑ニ俣群拖ɑ鹄炻柿考壓侠砣缓笤俜诺饺S網(wǎng)格里跑。這樣可以快速排除機(jī)理文件格式錯誤、組分熱物性缺失這類低級問題。另一個經(jīng)驗是如果目標(biāo)只關(guān)注NOx排放趨勢而不追求絕對精度可以選擇帶簡化NOx子機(jī)理的燃燒模型在流場收斂后用后處理方式評估NO生成但對于燃料型NO占主導(dǎo)的氨燃燒這種后處理方式的誤差很大更可靠的辦法還是把氮氧化物相關(guān)組分包含在主反應(yīng)機(jī)理中進(jìn)行全耦合計算。3. 燃燒模型與湍流相互作用模型選型的關(guān)鍵判斷3.1 為什么渦耗散模型在氨燃燒上容易翻車很多做過甲烷燃燒仿真的工程師習(xí)慣性選用渦耗散模型或渦耗散概念模型來處理湍流燃燒這套思路在天然氣燃燒上表現(xiàn)不錯因為甲烷反應(yīng)活性高可以近似認(rèn)為化學(xué)反應(yīng)速率遠(yuǎn)大于混合速率燃燒由湍流混合控制。但氨燃燒完全不同它的化學(xué)反應(yīng)時間尺度和湍流混合時間尺度處于同一數(shù)量級甚至是反應(yīng)速率更慢。這時候如果還套用渦耗散模型火焰會被計算得過于“劇烈”預(yù)測的火焰長度偏短、溫度偏高NOx結(jié)果更是偏離實際。主題083項目里我專門做了對比實驗同一套幾何、同一套網(wǎng)格分別用渦耗散模型和考慮有限速率的模型計算前者的出口平均溫度比后者高出接近一百五十開爾文NO濃度差了兩倍以上。這個差別足以改變工程結(jié)論。所以做氨燃燒仿真燃燒模型必須能同時體現(xiàn)湍流混合和化學(xué)動力學(xué)控制兩方面的作用有限速率模型或者兼顧二者的渦耗散概念模型才是正確起點。3.2 EDC模型與有限速率模型的實測對比在項目主體計算中我最終采用了渦耗散概念模型配合骨架機(jī)理。這個模型的核心假設(shè)是化學(xué)反應(yīng)發(fā)生在精細(xì)結(jié)構(gòu)區(qū)域內(nèi)通過能量耗散率計算精細(xì)結(jié)構(gòu)質(zhì)量分?jǐn)?shù)和反應(yīng)時間尺度既能反映湍流對反應(yīng)的促進(jìn)又能保留化學(xué)機(jī)理的動力學(xué)特征。實際使用中它的計算量比純有限速率模型略高但穩(wěn)定性更好尤其是對氨氫混合燃料這種具有雙燃料特性的體系既能捕捉氫組分快速反應(yīng)放熱的特點也能體現(xiàn)氨組分緩慢氧化的特征。如果選擇直接有限速率模型對網(wǎng)格質(zhì)量的要求會更高反應(yīng)區(qū)附近需要足夠細(xì)的網(wǎng)格來分辨火焰鋒面結(jié)構(gòu)否則化學(xué)反應(yīng)源項完全不準(zhǔn)確。這在三維工程案例里通常是巨大的網(wǎng)格量和迭代開銷。我的建議是如果燃燒室的雷諾數(shù)高、網(wǎng)格量在三百萬以上優(yōu)先選渦耗散概念模型如果做的是小型實驗燃燒器或?qū)恿骰鹧婺M再考慮層流有限速率模型。另外要留意模型參數(shù)中的精細(xì)結(jié)構(gòu)體積分?jǐn)?shù)常數(shù)默認(rèn)值在大渦模擬和雷諾平均湍流模型下適用范圍不太一樣最好根據(jù)實驗點火距離做一次標(biāo)定。3.3 邊界與初始條件的設(shè)置參數(shù)氨燃燒仿真對邊界條件非常敏感特別是入口溫度和當(dāng)量比。氨的化學(xué)反應(yīng)活性低入口溫度差個五十開爾文可能導(dǎo)致著火位置移動很多甚至出現(xiàn)點火失敗。主題083項目中空氣入口溫度設(shè)定為六百五十開爾文燃料入口溫度三百開爾文這符合實際燃?xì)廨啓C(jī)燃燒室的溫度級配。壓力出口設(shè)定一點六兆帕和水蒸氣、氮氣、氧氣、氨、氫氣等組分的回流條件一起設(shè)置。必須提醒的一點是燃料入口的組分設(shè)置不要只填氨工程上為了穩(wěn)燃通常摻入體積分?jǐn)?shù)百分之二十到三十的氫氣這個比例直接影響反應(yīng)活性和火焰穩(wěn)定特征。湍流參數(shù)也要認(rèn)真對待。入口湍流強(qiáng)度我設(shè)的是百分之十左右水力直徑按實際通道尺寸計算。有人會隨手填個百分之五了事這在甲烷燃燒仿真里影響不大但在氨燃燒里可能改變回流區(qū)尺寸進(jìn)而影響火焰位置。初始條件方面如果不做點火模擬最好先做一個不帶化學(xué)反應(yīng)的冷態(tài)流場計算收斂后引入能量源項點火再切到反應(yīng)流計算。直接冷啟動反應(yīng)流很容易導(dǎo)致流場被化學(xué)源項帶偏產(chǎn)生非物理解。4. 實操過程虛擬燃?xì)廨啓C(jī)燃燒室氨氫燃燒仿真全流程4.1 幾何模型與網(wǎng)格策略主題083項目使用的是某型實驗燃燒室的三維模型幾何上包括旋流器、主燃區(qū)、摻混孔和過渡段。原始幾何是從實驗臺圖紙簡化的去掉了一些螺栓孔、臺階等不影響流場的細(xì)節(jié)。簡化幾何這一步非常關(guān)鍵保留太多細(xì)微結(jié)構(gòu)會讓網(wǎng)格量大增但保留太少又會丟失旋流特征。我的做法是保留旋流葉片角度和通道面積忽略圓角和倒角燃燒室主體只用結(jié)構(gòu)化六面體網(wǎng)格劃分旋流器區(qū)域用非結(jié)構(gòu)化網(wǎng)格過渡。網(wǎng)格策略上燃燒反應(yīng)區(qū)是加密重點。對于氨氫混合燃料火焰鋒面比純甲烷更薄對網(wǎng)格分辨率的要求更高。主流區(qū)網(wǎng)格控制在三到五毫米主燃區(qū)剪切層附近加密到一毫米以下近壁面使用壁面函數(shù)第一層網(wǎng)格高度按y在三十到八十之間控制?？偩W(wǎng)格量約四百二十萬這個規(guī)模在工程工作站上可以接受。網(wǎng)格無關(guān)性驗證我做了三套分別是兩百四十萬、四百二十萬和七百萬結(jié)果發(fā)現(xiàn)溫度和速度場差別不大但NO濃度在粗網(wǎng)格和中等網(wǎng)格之間有明顯變化說明氮氧化物對網(wǎng)格分辨率格外敏感最終取了中等偏細(xì)的版本作為基準(zhǔn)算例。4.2 物性參數(shù)與入口條件物性參數(shù)設(shè)定直接影響燃燒計算的準(zhǔn)確性。氨、氫、氮、水蒸氣、氧、一氧化氮等組分的比熱容采用溫度分段多項式擬合熱導(dǎo)率和黏度用動力學(xué)理論計算?；旌衔锏臒嵛镄园蠢硐霘怏w混合定律處理。這種細(xì)節(jié)很多新手不注意直接用常數(shù)比熱容去算結(jié)果火焰溫度偏差非常大。氨的比熱容隨溫度變化很明顯常溫到兩千開爾文區(qū)間內(nèi)數(shù)值變化超過百分之三十用常數(shù)近似會造成顯著誤差。邊界條件具體數(shù)值如下空氣入口質(zhì)量流量每秒零點三千克溫度六百五十開爾文組分按空氣標(biāo)準(zhǔn)比例燃料入口質(zhì)量流量按目標(biāo)熱功率換算氨和氫的摻混比例為體積比百分之七十比三十溫度三百開爾文。當(dāng)量比設(shè)定為零點九即燃料略偏富燃這是抑制NOx的常用手段之一。出口采用壓力出口表壓為零回流溫度為兩千米。壁面采用無滑移絕熱條件實際燃燒室有冷卻結(jié)構(gòu)但項目前期只做熱態(tài)性能評估絕熱假設(shè)可以接受后續(xù)如果要研究壁面熱負(fù)荷再單獨加共軛傳熱模型。4.3 求解設(shè)置與收斂判斷求解器我用的是某商業(yè)CFD軟件的密度基或壓力基求解器均可處理實際項目采用壓力基耦合算法。離散格式上動量方程用二階迎風(fēng)湍流方程用一階迎風(fēng)起步、穩(wěn)定后切換二階能量和組分方程始終用二階迎風(fēng)。化學(xué)反應(yīng)源項的雅可比矩陣更新頻率要控制好默認(rèn)更新間隔在反應(yīng)初期可能造成發(fā)散我調(diào)整每一到兩步迭代更新一次源項穩(wěn)定后放寬到每五步一次。收斂判斷方面分三個層次。第一層看常規(guī)殘差能量方程降到十的負(fù)六次方以下流動方程十的負(fù)四次方以下。第二層監(jiān)控燃燒室出口平均溫度和出口NO質(zhì)量分?jǐn)?shù)的變化曲線當(dāng)這兩條曲線在兩百步迭代內(nèi)波動小于百分之一時基本認(rèn)為統(tǒng)計穩(wěn)定。第三層檢查全場質(zhì)量守恒和能量守恒進(jìn)出口質(zhì)量流量偏差控制在百分之零點五以內(nèi)能量不平衡不超過總輸入功率的百分之二。這里要特別說明氨燃燒仿真中殘差曲線很容易出現(xiàn)振蕩尤其是NO組分方程因為它的生成速率對局部溫度極其敏感不必強(qiáng)求殘差單調(diào)下降只要監(jiān)控量穩(wěn)定且守恒達(dá)標(biāo)就可以接受。5. 常見問題排查與經(jīng)驗速查表5.1 點火失敗與火焰吹脫氨燃燒仿真里最常遇到的第一個坎就是點火失敗。表現(xiàn)為溫度場基本沒有反應(yīng)燃料和氧化劑只是單純混合出口組分和入口組分差不多。原因無非三類化學(xué)機(jī)理活性不足、點火源能量不夠、初始流場不適合反應(yīng)建立。排查順序是先用零維反應(yīng)器確認(rèn)機(jī)理在該溫度壓力下能著火再檢查三維模型點火區(qū)的溫度和組分是否達(dá)到著火條件最后看點火源區(qū)域網(wǎng)格分辨率夠不夠。第二個常見問題是火焰吹脫即反應(yīng)區(qū)被高速氣流吹向出口或火焰貼壁后難以維持。氨燃料本身火焰?zhèn)鞑ニ俣鹊瓦@個問題尤其突出。如果你算出來的火焰位置和實驗高速攝影結(jié)果對不上優(yōu)先懷疑湍流燃燒模型對反應(yīng)速率的估計偏高或偏低。我遇到過一次渦耗散概念模型標(biāo)定參數(shù)不當(dāng)導(dǎo)致火焰附著位置偏移的情況調(diào)整精細(xì)結(jié)構(gòu)體積分?jǐn)?shù)后結(jié)果立刻改善。另外在燃料中增加氫氣摻混比例可以顯著提升火焰穩(wěn)定性如果工程條件允許可以做一個摻混比例的參數(shù)掃描找到一個既穩(wěn)定又低排放的工況點。5.2 NOx預(yù)測偏差大的原因與校準(zhǔn)NOx預(yù)測偏差是氨燃燒仿真中最難處理的問題。氨燃燒的NOx主要來自燃料型路徑氨先氧化為NO一部分NO會通過再燃反應(yīng)還原為氮氣這個還原路徑對局部混合和溫度分布特別敏感。仿真中如果網(wǎng)格在反應(yīng)區(qū)不夠細(xì)或者湍流模型對混合強(qiáng)度的預(yù)測有偏差NO結(jié)果就會差得離譜。主題083項目中我對比了實驗測點數(shù)據(jù)發(fā)現(xiàn)主燃區(qū)下游的NO濃度實驗值和仿真值偏差約百分之三十五這是一個可以接受的工程精度但如果直接看出口NO濃度偏差很容易達(dá)到兩倍。校準(zhǔn)方法有幾個方向。一是查看溫度場的空間分布是否合理NO生成速率強(qiáng)烈依賴溫度溫度偏高往往導(dǎo)致NO偏高。二是檢查N2O和NO2的中間產(chǎn)物分布如果某個中間組分濃度異常高說明反應(yīng)路徑失衡。三是用不同當(dāng)量比做敏感性分析在當(dāng)量比零點八到一點零之間掃描看NO峰值是否落在富燃側(cè)趨勢正確說明機(jī)理和模型組合基本可靠。此外還有一個實操技巧計算分離渦或大渦模擬所需的亞格子化學(xué)源項投影時NO結(jié)果對反應(yīng)進(jìn)度變量的表查詢精度很敏感有條件的話可以提高查表分辨率。5.3 調(diào)試心得與參數(shù)敏感性參數(shù)敏感性分析這件事做氨燃燒仿真必須認(rèn)真對待這不是做給評審看的“面子工程”而是找到仿真可信邊界的重要方式。我重點看四個參數(shù)入口溫度、氫氣摻混比例、當(dāng)量比和湍流強(qiáng)度。入口溫度每增加五十開爾文火焰位置前移明顯NO生成量上升氫氣摻混比例從零增加到百分之三十火焰穩(wěn)定性和反應(yīng)區(qū)溫度同步上升但NO先降后升存在一個最優(yōu)區(qū)間當(dāng)量比從零點八變到一點零出口NO通常呈現(xiàn)先升后降的趨勢湍流強(qiáng)度加大后混合增強(qiáng)導(dǎo)致峰值溫度降低但局部高溫區(qū)面積不一定減小?，F(xiàn)象可能原因排查建議點火失敗機(jī)理活性不足或初始溫度過低檢查零維著火延遲提高點火區(qū)初始溫度火焰吹脫到出口燃燒模型反應(yīng)速率偏低調(diào)整精細(xì)結(jié)構(gòu)常數(shù)增大氫氣摻混比例出口NO過高網(wǎng)格分辨率不足或機(jī)理缺再燃路徑加密反應(yīng)區(qū)網(wǎng)格檢查N2O中間產(chǎn)物分布溫度持續(xù)振蕩源項更新頻率過高降低雅可比更新頻率增大時間步氨剩余濃度過高停留時間不足或當(dāng)量比偏富調(diào)整燃料分配檢查回流區(qū)結(jié)構(gòu)收斂但組分不合理入口組分定義錯誤復(fù)核各入口組分質(zhì)量分?jǐn)?shù)總和是否為1調(diào)試過程中我最大的體會是氨燃燒仿真必須頻繁回到“物理判斷”上來。殘差單調(diào)下降、收斂得很漂亮不代表結(jié)果可信。一個明顯不合物理的壁面高溫區(qū)或者奇怪的組分分布很可能對應(yīng)的是模型參數(shù)匹配問題。每次修改參數(shù)后先別急著跑完整算例用零維或一維預(yù)混火焰快速驗證一下反應(yīng)性趨勢再投入三維計算能節(jié)省大量機(jī)時。6. 后處理與工程判斷仿真的終點不在云圖6.1 關(guān)鍵指標(biāo)的提取方法仿真計算收斂后最忌諱的是做完一張溫度云圖就收工。氨燃燒仿真的價值在工程判斷必須把計算結(jié)果轉(zhuǎn)化成可對比的參數(shù)指標(biāo)。我通常提取四組數(shù)據(jù)出口截面面積平均溫度、出口NO與N2O排放濃度、燃燒效率、火焰長度。出口截面面積平均溫度用面積加權(quán)平均要避免直接用質(zhì)量加權(quán)因為燃燒室出口速度分布不均勻面積平均更能反映壁面耐溫需求。NO濃度換算成百萬分之一以百分之十五氧基準(zhǔn)折算方便和排放法規(guī)對照。燃燒效率的計算在氨氫混合燃燒時比純碳?xì)淙剂蠌?fù)雜不能只用二氧化碳生成量判斷。我采用燃料剩余法統(tǒng)計出口處氨和一氧化碳等未燃組分的化學(xué)能總和用總輸入化學(xué)能減去未燃化學(xué)能再除以總輸入化學(xué)能。主題083基準(zhǔn)工況下的燃燒效率約百分之九十九點二這符合預(yù)期。火焰長度通過提取反應(yīng)進(jìn)度變量或者OH自由基質(zhì)量分?jǐn)?shù)等值面來判斷這個指標(biāo)對燃燒室長度設(shè)計非常有用。6.2 數(shù)據(jù)可靠性判斷做完仿真后還要做數(shù)據(jù)可靠性審查我的習(xí)慣是把仿真結(jié)果和實驗數(shù)據(jù)或文獻(xiàn)經(jīng)驗放在一起對比。幾個硬性檢查點絕熱火焰溫度是否和化學(xué)平衡計算值接近如果不加入熱損失項主流區(qū)最高溫度不應(yīng)該明顯超過理論絕熱溫度出口氧氣濃度和入口當(dāng)量比是否匹配NO生成趨勢是否隨摻氫比例變化形成非單調(diào)曲線。每個檢查點都過了仿真結(jié)果才算真正可以用來指導(dǎo)工程決策。另外想多說一句氨燃燒仿真的數(shù)據(jù)量非常大尤其是組分多、反應(yīng)路徑復(fù)雜后處理時不要把二十個組分全都拉出來看。我一般重點關(guān)注溫度、速度、氨、氫、水蒸氣、NO、N2O、NO2這八個量其他中間產(chǎn)物只在問題排查時再看。看太多變量反而容易迷失重點工程問題需要的是關(guān)鍵矛盾的判斷不是科研級別的全路徑分析。7. 火候與分寸關(guān)于氨燃燒仿真的一些實在體會最后分享一點個人經(jīng)驗。氨燃料燃燒仿真的難點不在于某個單一環(huán)節(jié)而在于整個鏈條的匹配。機(jī)理、燃燒模型、網(wǎng)格、邊界條件每一項單獨看都有成熟方法但它們組合在一起時會產(chǎn)生獨特的耦合效應(yīng)。我在主題083項目里最大的收獲就是建立了一套“機(jī)理先驗證、模型再對比、網(wǎng)格做無關(guān)性、參數(shù)做敏感性”的完整流程每一步都不跳。這套流程聽起來繁瑣但它能讓你在結(jié)果不可信時迅速定位問題出在哪一環(huán)。還有一個小技巧氨燃燒仿真項目啟動時建議先找一個實驗工況點進(jìn)行全流程復(fù)算把計算值、實驗值放在一起校準(zhǔn)好模型參數(shù)再推廣到其他工況。直接上來就跑設(shè)計工況出現(xiàn)問題根本分不清是模型問題還是工況特點。把驗證工況跑順之后后面的參數(shù)掃描、優(yōu)化設(shè)計都會順很多。氨燃料的仿真還會持續(xù)發(fā)展新的機(jī)理、新的燃燒組織方式都在不斷出現(xiàn)核心能力不是記住某個模型而是能判斷什么條件下該信任仿真、什么條件下必須回看物理規(guī)律。