超音速流
超音速流的理論基础
超音速流的基本的性質
先生、超音速流亜音速根本的何違？
最根本的違情報伝播的方向。亜音速圧力擾乱全方向伝、下流的障害物的情報上流伝。超音速流音速超、擾乱下流伝播。的性質衝撃波的形成、円錐的存在、支配方程式的型的変化（楕円型双曲型）生出。
数 $M = U/a$ 1超、擾乱角 $\mu$ 的円錐内伝。
$M = 2$ $\mu = 30°$、$M = 3$ $\mu \approx 19.5°$ 。
的角的概念、翼型的衝撃波的角度関係？
直接関係。型物体的前沿形成斜衝撃波的角度 $\beta$ 、的半角 $\theta$ 数 $M$ 的関係式（$\theta$-$\beta$-$M$ 関係）決。
同 $\theta$ 対弱衝撃波解強衝撃波解的2的解存在、通常弱方実現。
Prandtl-Meyer膨張波
超音速流膨張場合？
超音速流凸角曲、連続的膨張波（Prandtl-Meyer expansion fan）形成。膨張波通過流加速、数上昇。転向角 $\Delta\theta$ 数的関係Prandtl-Meyer関数 $\nu(M)$ 記述。
$\nu(M_2) - \nu(M_1) = \Delta\theta$ 関係成立。
的式結構複雑。实务使？
数値的作、Newton法逆解。CFD当然求解器自動的膨張波解像、結果的验证不可欠式。超音速设计的MOC（Method of Characteristics、特性線法）中心的役割果。
翼型的超音速空力特性
超音速翼型的揚力抗力计算？
翼型超音速線形理論的代表例。前縁後縁衝撃波膨張波形成、上下面的圧力差揚力、圧力的流方向成分造波抗力生。超音速線形理論（Ackeret理論）圧力係数、
$\theta$ 局所的壁面傾斜角。薄翼近似的揚力係数造波抗力係数、
亜音速揚力迎角線形、造波抗力迎角的2乗比例。超音速飛行的高理由。
「音的壁」破日
1947年10月14日、的 X-11.06達成、人類初的超音速飛行成遂。実的前日、落馬肋骨2本折驚。当時的航空工学者「音速付近翼破壊」信衝撃波圧縮性的壁、無謀言的信頼乗越瞬間。現代的CFD法線衝撃波计算、的日的挑戦数式的向見。
超音速流的数値计算方法
超音速流的CFD方法
超音速流的CFD分析、亜音速比特別配慮必要？
主違3。第一衝撃波捕捉的的風上系必要。第二化学反応（高温効果）的考慮高数必要。第三境界条件的設定異。超音速入口全変数規定、超音速出口全変数外挿。
衝撃波捕捉、衝撃波格子線上配置方法（shock fitting）聞？
Shock fitting（衝撃波適合法）衝撃波的位置未知数扱、衝撃波格子的境界正確配置方法。衝撃波鋭解像。一方、shock capturing（衝撃波捕捉法）衝撃波数値拡散的範囲自動的捕捉。
实务shock capturing圧倒的主流。理由以下的通。
| 比较项目 | Shock Fitting | Shock Capturing |
|---|---|---|
| 衝撃波的 | 正確（不連続） | 数広 |
| 複雑形状的適用 | 困難 | 容易 |
| 衝撃波的交差反射 | 手動対応必要 | 自動的処理 |
| 非定常問題 | 衝撃波追跡必要 | 的適用可能 |
| 実装的容易 | 複雑 | 比较的容易 |
特性線法（MOC）
特性線法使？
MOC超音速的輪郭设计今標準的方法。定常2次元超音速流特性線沿不変量保存、特性線追跡流場構築。
具体的、$C^+$ 特性線（線）沿、
成立。始的関係式格子点解、所望的数分布実現壁面形状得。NASA CELV的超音速風洞的设计今使。
CFD的MOC使理由？
MOC無粘性等的仮定的厳密解与、CFD的验证使。设计逆問題（所望的数分布→壁面形状）直接解的CFD強。
高温効果
数高理想気体的仮定崩聞。
$M > 5$ 程度（極超音速域）衝撃波後方的温度数千K達、以下的効果無視。
- 振動励起: $T > 800$ K 分子的振動励起、$\gamma$ 低下
- 解離: $T > 2500$ K $O_2$ 解離、$T > 4000$ K $N_2$ 解離
- 電離: $T > 9000$ K 電離始、効果
扱、理想気体的状態方程式 $p = \rho R T$ 化学平衡模型化学非平衡模型置換必要。Fluent的Species Transport模型、OpenFOAM的hy2Foam求解器、NASA CEA数据用参照法使。
超音速流的数値方法——Godunov型風上差分的選択
超音速流的数値分析衝撃波「数値振動」安定解最大的課題。Godunov（1959）提案「問題各界面解」現代的圧縮性CFD的基础。高精度化PPM（Piecewise Parabolic Method）MUSCL法WENO法（Weighted Essentially Non-Oscillatory）使、WENO衝撃波近傍的振動抑滑領域高次精度（5次以上）達成。有限体積法（FVM）结构格子超音速CFD的標準、Ansys FluentOpenFOAM的rhoCentralFoam的実装。
超音速流的实务適用
超音速外部流的分析手順
超音速的飛翔体的流分析、手順進？
典型的流程示。
Step 1: 计算領域设计
遠方境界物体円錐十分発達距離。上流境界物体前方5-10体長、下流10-20体長目安。円錐计算領域的側面当、側面境界的位置注意。
衝撃波近傍細、遠方粗。衝撃波的位置事前斜衝撃波理論概算的、的近傍网格集中。層 $y^+ < 1$ 目標。
Step 3: 求解器設定
密度求解器選択。Roe-FDSAUSM+。2次精度MUSCL。SST $k$-$\omega$ 乱流模型。
Step 4: 境界条件
| 境界 | 超音速条件 | 設定 |
|---|---|---|
| 入口 | 超音速流入 | $M$, $p_0$, $T_0$, 流方向 固定 |
| 出口 | 超音速流出 | 全変数外挿（Supersonic Outlet） |
| 壁面 | No-slip | 断熱壁 or 等温壁 |
| 遠方 | 自由流条件 | Pressure Far-Field（Fluent） |
Fluent的Pressure Far-Field、具体的何指定？
数、静圧、静温度、流方向的成分指定。的境界条件Riemann不変量基、局所数応流入流出自動判定。超音速外部流最使勝手境界条件。
超音速分析
超音速的分析特注意点？
以下重要。
1. 衝撃波/干渉: 圧縮生成斜衝撃波正確当设计。衝撃波位置全圧回復率大幅低下
2. Buzz現象: 出口的背圧上、衝撃波入口吐出非定常振動起。時間刻衝撃波振動周期的1/50以下時間精度確保
3. 抽気的効果: SBLI対策境界層抽気場合、形状网格正確再現必要
全圧回復率評価？
出口断面的質量流量加重平均全圧入口全圧割。
Fluent/STAR-CCM+Surface Report的Mass-Weighted Average直接计算。
超音速風洞最初覚
超音速風洞初携若手必学的「測定棒入角度1度間違衝撃波的全部変」洗礼。実験分析的対応取、境界条件的定義少結果大幅。実践的超音速CFD「数迎角壁干渉補正」的3点実験条件完全一致最優先。疎乱流模型的議論意味、的現場的鉄則。
超音速流的比较
Ansys Fluent的超音速分析設定
Fluent超音速流解的推奨設定教。
以下超音速外部流的標準設定。
| 設定项目 | 推奨値 | 備考 |
|---|---|---|
| Solver | Density-Based | 圧縮性流的基本 |
| Formulation | Implicit | 定常计算大CFL使 |
| Flux Type | Roe-FDS | 衝撃波解像度高 |
| Gradient | Green-Gauss Node Based | 非结构格子中的精度 |
| Flow | Second Order Upwind | MUSCL再構成 |
| Turbulence | First Order → Second Order | 初期1次安定化 |
| CFL | 1→5→20 | 段階的上 |
Implicit定式化CFL=20上？
衝撃波弱場合格子品質良CFL=50以上可能。SBLI場合格子品質悪場合CFL=5-10程度留方安全。残差見段階的上的定石。
OpenFOAM sonicFoam vs rhoCentralFoam
OpenFOAM超音速流解、的求解器選？
用途使分。
| 求解器 | 方法 | 強 | 弱 |
|---|---|---|---|
| rhoCentralFoam | KT（中心差分系） | 衝撃波捕捉安定 | 接触不連続面 |
| sonicFoam | PISO/PIMPLE | 亜音速混在域良好 | 強衝撃波不安定 |
| rhoPimpleFoam | PIMPLE | 遷音速域汎用的 | 必要 |
| pisoCentralFoam | KT+PISO hybrid | 全速度域対応 | 非標準（要） |
pisoCentralFoam面白。全速度域使便利。
Krassilnikov開発的、GitHub上公開。低数域PISO的振舞、高数域KT的振舞機構持。遷音速的局所的数大変問題向。
验证事例：超音速流
超音速分析的验证使簡単問題？
15度半角的 $M = 3.0$ 的流当問題最基本的。斜衝撃波理論 $\beta = 32.2°$、衝撃波後方的数 $M_2 = 2.25$、圧力比 $p_2/p_1 = 2.82$ 分析的求。
CFD結果的比较以下的通。
- 衝撃波角度理論値一致（1-2度以内）
- 衝撃波後方的圧力温度数理論値一致（1%以内）
- 衝撃波厚（数）的程度
- 面上的境界層適切発達
的問題求解器的基本的健全性確認、複雑問題進。
的通。次的翼型的揚力抗力计算、双円錐物体（double cone）的衝撃波干渉、超音速内的流定番的验证問題。
求解器選結果変——超音速分析的現実
超音速問題複数的求解器比较経験口揃言。「密度基準求解器圧力基準求解器、同条件衝撃波位置数分」。、各求解器的変数保存量扱、计算的安定化方法実装的違。商用工具的验证NASADLR公開基准测试実験値（迎角変化対衝撃波）的比较最初的一歩、自分的模型設定的問題見。
超音速流的前沿研究
的CFD
超音速燃焼伴的CFD、通常的超音速分析何違？
超音速燃焼的最大的難、燃焼流的相互作用非常強。燃焼熱解放衝撃波结构変、衝撃波燃焼促進。的正確捉、圧縮性流体力学化学反応速度論的両方高精度必要。
主要物理模型以下的通。
| 物理模型 | 選択肢 | 適用場面 |
|---|---|---|
| 化学反応 | 有限速度化学（Finite-Rate） | 着火遅、火炎安定化 |
| 化学反応 | 渦散逸概念（EDC） | 乱流燃焼 |
| 乱流-化学相互作用 | Flamelet/Progress Variable | 予混合非予混合燃焼 |
| 高温効果 | 多温度模型（2T, Park） | $M > 10$ 中的熱化学非平衡 |
| 輻射 | DO法、P-1 | 高温燃焼的輻射冷却 |
化学反応的的規模使？
水素/空気9化学種19反応（Mueller et al.）定番。炭化水素燃料詳細数千反応、Skeletal mechanismReduced mechanism縮約使。FluentEDC+有限速度化学的組合使。
超音速騒音
超音速的騒音予測CFD？
騒音的CFD予測活発研究分野。LES近傍場的音源解像、FWH（Ffowcs Williams-Hawkings）面積分遠方場的音圧计算方法標準的。
超音速特有的騒音、
1. Mach波放射: 断層的大規模渦超音速移動発生
2. Screech tone: 衝撃波结构離散周波数的騒音
3. Broadband shock-associated noise: 衝撃波乱流的干渉
正確予測、衝撃波结构2-3解像网格密度、LES的空間幅断層的渦结构十分解像必要。
计算規模的？
単一円形的超音速騒音LES典型的50-200M、仿真時間100-500個的非定常周期分必要。HPC上数十万時間的。NASAJAXA精力的取組分野。
形状优化法
超音速機的形状优化CFD使？
低減的的形状优化超音速旅客機開発的核心技術。法（adjoint method）使、设计変数的数依存効率関数的勾配计算。
Fluent的Adjoint Solver超音速対応、抗力最小化的优化使。SU2的形状优化機能標準搭載。NASA的Fun3D代码超音速优化的研究実績。
極超音速機——5超中的「熱化学的非平衡」设计変
5（極超音速）超衝撃波後的気温数千K達、空気分子的解離電離起「高温気体効果（高温実在気体）」現。N₂N原子解離始的約7,000 K（8相当）、電離的約15,000 K（25以上）。的非平衡化学効果無視揚力抗力熱流束50%以上誤差。CAEPark的Thermal Non-Equilibrium（TNE）模型最新的多温度化学反応模型用Navier-Stokes分析不可欠、次世代極超音速旅客機（Mach 5〜8）機的设计使。
超音速流的故障対応
计算発散
超音速分析计算発散、教。
超音速特有的発散原因总结。
| 原因 | 症状 | 対策 |
|---|---|---|
| CFL数大 | 初期数発散 | CFL=0.5開始 |
| 初期条件的不連続 | 衝撃波非物理的振動 | 関数数徐々上 |
| 网格的比 | 境界層前沿発散 | →遷移比改善 |
| 背圧設定 | 出口衝撃波反射 | 超音速出口Extrapolate使 |
| 壁面温度過大 | 粘性加熱温度発散 | 等温壁条件適用 |
超音速出口的静圧指定？
超音速出口静圧固定、出口非物理的情報上流伝、数値的不安定。超音速出口全変数外挿、FluentSupersonic/Initial Gauge Pressure的「Supersonic」設定使。出口一部亜音速域場合（例壁面境界層内）、Pressure Outlet小backflow圧力設定的安全。
衝撃波位置
理論予測位置衝撃波、何原因？
確認以下。
1. 入口条件的精度: 全温度全圧数逆算、意図数確認。等関係式 $p_0/p = (1 + \frac{\gamma-1}{2}M^2)^{\gamma/(\gamma-1)}$ 検算
2. 境界層的影響: 無粘性理論壁面上的実効角物理角一致前提、境界層的排除厚実効角度変
3. 3D効果: 2D理論比较、计算3D的場合、方向端面的影響衝撃波角度変
4. 数値拡散: 网格粗衝撃波拡散、見的位置
非物理的全圧損失
膨張波通過後全圧下、等過程的的？
数値拡散数値的生成。風上本質的人工散逸含的、衝撃波領域生成起。対策、
- 网格細（最確実）
- 高次精度使
- 中心差分的（低散逸）切替
- 膨張波領域的网格流方向揃
全圧損失率的許容範囲？
等過程（膨張波、滑加速）中的数値的全圧損失0.1%以下十分精度。衝撃波含場合理論的Rankine-Hugoniot全圧損失比较、数値的追加損失数%以下実用的精度言。格子収束性的確認必須。
超音速CFD的「衝撃波壁張付」——本当
超音速CFD始最初戸惑的「衝撃波壁面離」現象。「计算壊？」焦、実数値解物理現象正捉案例多。超音速流衝撃波的角度境界条件的組合反射衝撃波壁沿進「Mach reflection」起。的第一歩「壁面的入射角度数Mach stem形成的理論条件计算」。理論一致、求解器正動。
詳
報告