爆轰（爆轰波）
爆轰（爆轰波）的理論基础
概述
教授，爆轰與爆炸有什麼區別？
燃燒傳播有兩種形式。爆燃（缺陷燃燒）是亞音速火焰傳播，普通燃燒器火焰屬於此類。另一方面爆轰（爆轰波）是衝擊波與燃燒波結合形成的超音速傳播現象。傳播速度在甲烷/空氣中約1800 m/s，在氫氣/空氣中約2000 m/s。
衝擊波和燃燒一起傳播，是這樣嗎？
是的。先行衝擊波對未燃混合氣進行絕熱壓縮以升高溫度，高溫下化學反應迅速進行，反應能量維持衝擊波。這種自我維持機制是爆轰的本質。
Chapman-Jouguet理論
Chapman-Jouguet（CJ）理論是什麼？
CJ理論是熱力學地求出爆轰波傳播速度的理論。在衝擊波前後的Rankine-Hugoniot關係式中加入化學反應放熱量 $q$，應用CJ條件（爆轰波後流速恰好等於局部音速）。
用簡化形式，CJ爆轰速度大致可以寫成：
這裡 $\gamma$ 是比熱比，$q$ 是單位質量的放熱量 [J/kg]。
還有CJ馬赫數的概念吧？
CJ爆轰的馬赫數由以下公式給出：
氫氣/空氣（當量比1.0）的情況下 $M_{CJ} \approx 5.0$，甲烷/空氣為 $M_{CJ} \approx 5.2$。
ZND結構
爆轰波的內部結構是什麼樣的？
ZND（Zel'dovich-von Neumann-Doering）模型中，爆轰波具有三層結構。
1. 衝擊波面（von Neumann spike）：未反應氣體受到衝擊壓縮。壓力達到CJ值的約2倍
2. 誘導帶（Induction zone）：化學反應進行前的延遲區間。對應著火延遲時間
3. 反應帶：迅速進行化學反應，達到CJ狀態
誘導帶越短爆轰越穩定，是這樣嗎？
是的。誘導帶長度 $\Delta_i$ 直接關係到細胞大小 $\lambda$，有經驗公式 $\lambda \approx (10-30)\Delta_i$。當這個細胞大小相對於爆轰管直徑足夠小（管直徑 > 數個$\lambda$）時，穩定爆轰傳播才能維持。
爆轰理論是衝擊波力學和化學反應速率論的融合呢。
完全正確。用CFD處理需要既能準確捕捉衝擊波的數值方法，又能處理高溫高壓下化學反應速率。
爆轰——爆轰波以「音速5～10倍」傳播的原因
燃燒有「爆燃（缺陷燃燒：音速以下火焰傳播）」和「爆轰（爆轰波：衝擊波與燃燒一體化的超音速傳播）」兩種。爆轰波由Chapman-Jouguet（CJ）理論（1899～1905年）描述，燃燒氣體在音速下的流動對應「CJ面」，其速度是特性值。氫-空氣混合氣約2 km/s，天然氣-空氣約1.8 km/s的爆轰速度。CAE中將歐拉方程與化學反應模型結合，用高解像度數值方案分析爆轰波傳播、細胞結構、轉移（DDT：爆燃到爆轰的轉移）。
爆轰（爆轰波）的數值计算方法
數值方法的詳細
用CFD求解爆轰需要什麼數值方法？
爆轰的數值分析與普通RANS燃燒分析大不相同。需要準確捕捉衝擊波的高解像度方案，同時兼顧高溫高壓下的詳細化學反應。
空間離散化
捕捉衝擊波適合用什麼方案？
衝擊波捕捉常用高階精度TVD（Total Variation Diminishing）方案或WENO（Weighted Essentially Non-Oscillatory）方案。
| 方案 | 精度 | 特點 | 應用 |
|---|---|---|---|
| Roe + Minmod | 2階 | 穩定但數值擴散大 | 初期研究 |
| HLLC | 2階 | 分析接觸不連續面 | 通用 |
| WENO-5 | 5階 | 高精度但计算成本大 | DNS/高精度计算 |
| MUSCL-Hancock | 2階 | 成本-精度平衡良好 | 實用爆轰计算 |
WENO方案最理想，但實務上用MUSCL就足夠了嗎？
正是。MUSCL+HLLC黎曼求解器的組合是實務的折中方案。但網格大小要保證為爆轰細胞寬度 $\lambda$ 的1/20以下。氫/空氣（當量比1.0, 1 atm）時 $\lambda \approx 10$ mm，所以網格寬度需要0.5 mm以下。
時間積分
時間積分怎麼處理？
爆轰波傳播是微秒級現象，通常用顯式時間積分。但化學反應部分是Stiff的，所以用算子分裂（Strang splitting）分離流體輸送和化學反應。
我給出典型設定值。
| 參數 | 推薦值 | 備註 |
|---|---|---|
| CFL數 | 0.3-0.5 | 衝擊波捕捉需保守選擇 |
| 化學反應求解器 | CVODE (BDF) | Stiff系統必需 |
| 最小時間步 | $10^{-9}$ s | von Neumann尖峰分析 |
| 網格寬度 | $\lambda/20$ 以下 | 細胞結構分析最低條件 |
AMR（適應格子細分化）
均勻的細網格會導致计算成本巨大吧？
這時AMR就大顯身手了。只在爆轰波面附近網格細分，未燃和已燃區域保持粗網格。CONVERGE內置自動AMR，OpenFOAM也可用 dynamicRefineFvMesh 實現。AMR可將記憶體和计算時間減少1-2個數量級。
AMR的細分基準用什麼？
通常用溫度梯度 $|\nabla T|$ 或壓力梯度 $|\nabla p|$ 作為細分傳感器。OH質量分率梯度也有效用於追蹤反應帶。
爆轰的數值计算是衝擊波捕捉方案 + Stiff化學反應 + AMR的三位一體呢。
完全同意。少了其中任何一個都無法進行實用爆轰计算。
爆轰计算需要「超音速方案」的原因——CFL數0.3的壁障
爆轰數值计算最初面對的困難是離散化方案。爆轰波以超音速（馬赫5～10等效）傳播，普通燃燒CFD用的2階精度中心差分方案會立即發散。解決辦法是必須使用WENO或Roe法等衝擊波捕捉方案，CFL數也需限制在0.3以下。這意味著時間步比普通燃燒计算小一個數量級，计算時間增加10倍以上。「爆轰是CAE计算中特別難處理的一類」的原因在此。
爆轰（爆轰波）的實務應用
實踐指南
請教我實際進行爆轰分析的步驟。
爆轰分析的计算工作量比普通燃燒CFD大得多。讓我整理一下實用流程。
分析流程
1. 0D/1D預備计算 -- 用Cantera或SD Toolbox计算CJ條件（傳播速度、壓力比）。從ZND結構估出誘導帶長度和細胞大小
2. 網格设计 -- 細胞大小$\lambda$的1/20以下的分析度。使用AMR時設定最大細分級別
3. 爆轰開始條件 -- 設置Driver section（高壓高溫區）或用熱點直接著火
4. 時間積分設置 -- CFL = 0.3-0.5，總计算時間直到爆轰波到達管端
5. 後處理 -- x-t圖、細胞結構可視化、壓力溫度分佈
RDE（回轉爆轰引擎）的設定示例
最近受關注的RDE如何分析？
RDE中環形燃燒室的爆轰波沿周向傳播。從2D展開模型（將環形展開為平面）開始计算，再擴展到3D比较有效率。
給出RDE分析的典型設定。
| 项目 | 設定值 | 備註 |
|---|---|---|
| 燃料 | H2/空氣，當量比1.0 | CH4著火延遲太長，RDE不穩定 |
| 環形直徑 | 100-200 mm | 研究規模 |
| 網格 | 0.2-0.5 mm | 推薦使用AMR |
| 反應機構 | Li et al. 9種/21反應 | H2詳細機構已足夠緊湊 |
| 入口條件 | 總壓指定 + 混合氣組成 | 推薦choked inlet防止逆流 |
| 计算時間 | 數毫秒（5-10轉） | 確認定常旋轉狀態 |
H2較多是因為著火延遲短吧。
是的。甲烷著火延遲是氫的100倍以上，難以維持回轉爆轰波。實用化方面正在檢討甲烷/氫混合、乙烯（C2H4）等。
常見失敗與對策
| 症狀 | 原因 | 對策 |
|---|---|---|
| 爆轰不傳播 | 網格太粗 | 細分到$\lambda/20$以下 |
| 傳播速度與CJ值不符 | 數值擴散大 | 改用高階方案(MUSCL/WENO) |
| 看不到細胞結構 | 2D计算分析度不足 | 用0.1 mm以下網格重新计算 |
| RDE波消失 | 入口混氣供給跟不上 | 增加入口面積或總壓 |
爆轰分析是计算成本的戰爭呢。
是的。用2D AMR掌握概況，再對必要部分進行3D詳細计算是現實做法。
爆轰管爆炸事故——「爆燃到爆轰的轉移（DDT）」的恐怖
化學工廠管道中氫氧混合氣洩漏時，最初是溫和燃燒（爆燃），但火焰在管內傳播時會加速並轉化為爆轰（爆轰波），這稱為DDT（Deflagration-to-Detonation Transition），是配管设计人員最擔心的現象之一。爆轰轉化的瞬間，壓力會跳升到爆燃的10～20倍。爆轰CFD實踐中，為了進行安全評估，必須準備足夠細膩的網格來正確模擬DDT，並準確追蹤從點火開始的時間演變。
爆轰（爆轰波）的軟體比较
商用工具比较
能進行爆轰分析的CFD工具有哪些？
爆轰需要衝擊波捕捉和詳細化學反應兩個條件，所以對應工具比较有限。
| 工具 | 衝擊波方案 | AMR | 詳細化學反應 | 適性 |
|---|---|---|---|---|
| CONVERGE | PISO + ALE | 自動AMR | SAGE | 最優 |
| Ansys Fluent | Density-Based | 無（UDF外部） | Stiff Chemistry | 可能但需要技巧 |
| OpenFOAM (rhoCentralFoam) | KNP/KT | dynamicRefine | ode | 靈活，需編碼 |
| STAR-CCM+ | Coupled Flow | 無 | DARS | 有限 |
| 專用代碼 (AMROC, PeleC) | Godunov類 | Block-AMR | 詳細/縮約 | 最適研究 |
CONVERGE對爆轰最優出乎意料，印象中它是內燃機的工具呢。
CONVERGE的自動AMR和網格切割技術對爆轰波面追蹤非常相容。省去了大量網格生成工作。實際上RDE和脈衝爆轰引擎（PDE）研究中CONVERGE應用增加。
Ansys Fluent中的注意點
用Fluent進行爆轰计算時的注意點？
OpenFOAM中的實現
OpenFOAM用哪個求解器？
rhoCentralFoam 是可壓流的標準求解器，但需自行實現與化學反應耦合。社區版有 rhoReactingCentralFoam。用 dynamicRefineFvMesh 啟用AMR，以溫度、壓力梯度為細分基準是定石。
選定指南
最後選哪個比较好？
專用代碼的存在不知道。爆轰特殊要求的工具選擇很重要呢。
是的。直接套用普通燃燒CFD工具會因數值擴散導致爆轰細胞結構消失。要根據工具特性選擇。
爆轰CAE的「客戶」是防衛、太空、安全評估——商用工具支持滯後的原因
爆轰模擬市場很小。客戶主要限於防衛機構、太空機構、工廠安全評估，所以汎用商用CFD工具（Fluent、STAR-CCM+）的爆轰專用功能有限。ANSYS的Fluent中可選WENO方案，但爆轰設置示例在官方手冊中很少。實際上爆轰计算常用OpenFOAM的XiFoam或blastFoam（開源），或SANDIA、Lawrence Livermore自主開發代碼。商用工具「對爆轰投資多少」取決於市場規模——CAE世界也由經濟原理驅動。
爆轰（爆轰波）的前沿研究
前沿話題與研究動向
爆轰最前沿研究有哪些？
主要有三個方向。(1)回轉爆轰引擎（RDE）實用化，(2)爆轰驅動的壓力增益燃燒，(3)傾斜爆轰波引擎（ODWE）。
RDE（回轉爆轰引擎）
為什麼RDE受關注？
RDE的爆轰波在環形燃燒室連續迴轉，可提供定常推力。相對傳統爆燃燃燒器，熵增加較小，理論上熱效率提升5-15%。
| 機構 | 研究成果 | 燃料 |
|---|---|---|
| JAXA | 世界首次太空RDE實證（S-520火箭） | 乙烯/O2 |
| 美空軍研究所 (AFRL) | 與GE共同將RDE用於噴氣引擎燃燒器 | JP-8/空氣 |
| 名古屋大學 | 連續迴轉爆轰穩定化方法 | H2/空氣 |
| KAUST | LES+詳細化學反應的RDE數值分析 | H2/空氣 |
JAXA在太空進行實證實驗了呢。
2021年S-520-31號火箭實現了世界首次太空RDE作動。CFD對设计做出了巨大貢獻。
傾斜爆轰波引擎（ODWE）
ODWE是什麼？
這是為極超音速飛行體（馬赫5-10）設想的推進裝置，在楔形燃燒室內駐定傾斜爆轰波。是SCRAMJET之後的概念。駐波不需RDE的迴轉，但需針對馬赫數優化楔角度。
DNS（直接數值模擬）
進行爆轰DNS嗎？
近年Peta-FLOPS級HPC使得2D爆轰細胞結構可用DNS分析。使用PeleC或S3D（Sandia）。3D DNS還需Exascale級计算資源，現由美DOE的ECP（Exascale Computing Project）推進中。
機器學習的應用
爆轰研究也用機器學習嗎？
爆轰研究直接關係航空推進最前沿呢。
是的。特別是RDE作為實用化技術各國競相研究。CFD设计支持的重要性今後還會增加。
回轉爆轰引擎——爆炸「持續迴轉」的逆轉思維
連續運用爆轰的回轉爆轰引擎（RDE）從2000年代起研究加速。普通燃燒是「定壓燃燒」，爆轰接近「定容燃燒」，故熱效率理論上高10～15%。JAXA與各國防衛機構投入資金的原因。但實用化困難重重，在環形室內維持超音速爆轰波穩定運行難度極高。現在前沿研究試圖用機器學習對爆轰波動態行為進行即時控制。「控制爆炸」這一矛盾命題正成為該領域正面攻關的主題。
爆轰（爆轰波）的故障排除
故障排除
請教爆轰计算常見的故障。
爆轰计算比普通燃燒CFD敏感得多。整理主要故障。
1. 爆轰波不開始
可能原因：
- Driver section壓力溫度不足
- 網格太粗，衝擊波與反應帶無法分析
- 反應機構未能正確再現該條件著火延遲
對策：
- 設定Driver條件為 $p > 30\,p_1$, $T > 2500$ K
- 先從1D驗證開始條件再進行2D/3D擴展
- 用ZND计算確認所需分析度（誘導帶 $\Delta_i$ 至少20點）
2. 傳播速度與CJ值偏差
傳播速度與CJ理論值不符的情況怎麼辦？
速度低的情況：數值擴散過大，反應帶能量散逸。改用高階方案（MUSCL、WENO），細分網格。
速度高的情況：過驅動狀態，初始Driver影響仍存。讓爆轰波傳播足夠距離後再測（管直徑10倍以上）。
| 偏差 | 允許範圍 | 要因 |
|---|---|---|
| $D/D_{CJ} < 0.95$ | 需改善 | 數值擴散、網格不足 |
| $0.95 < D/D_{CJ} < 1.02$ | 良好 | -- |
| $D/D_{CJ} > 1.05$ | 需改善 | 過驅動、過渡狀態 |
3. RDE特有的問題
RDE计算特有的問題有嗎？
4. 计算成本巨大
對策：
- 活用AMR（波面附近高分析，其他粗網格）
- 先2D掌握概況再進行3D詳細计算
- H2情況用Li et al.的9種/21反應就足夠（無需縮約）
- 利用GPU並列化求解器（PeleC支援GPU）
爆轰计算是分析度和成本的權衡很嚴峻呢。
是的。所以用1D ZND分析估算必要分析度，再用AMR降低成本的步驟很重要。直接上3D均勻網格的話，计算可能數月不完成。
「爆轰波消失了」——爆轰CFD重啟困難的原因
爆轰计算故障排除的特有問題是「计算中途爆轰波消失，降級為爆燃（低速燃燒）」的現象。這種「quenching（消炎失速）」在實驗中也發生，但计算中多因數值黏性導致爆轰波人為衰弱。典型是網格太粗無法分析爆轰細胞結構（detonation cell），經驗法則「需要細胞寬度1/5以下的網格」。但實際爆轰細胞寬度只有數毫米至數厘米，實機規模全域細化會成本不現實，形成矛盾。
相關話題
價值
更詳細
錯誤