衝撃波管問題（Riemannソルバー） — よくある問題と対策

不連続面近傍のオーバーシュート/アンダーシュートは高次精度スキームの宿命だ。対策の選択肢を整理しよう。

Fluentでは「Solution Methods」でSpatial DiscretizationをFirst Order Upwindに一旦落として収束させ、その解を初期値にしてSecond Orderに切り替えるという手順が有効だ。また、Gradient LimiterのMethodをStandardからMultidimensional Limiterに変更すると振動が抑制されることもある。

Carbuncle現象は、強い垂直衝撃波の数値解が格子に依存した非物理的な不安定構造を示す問題だ。Roe系スキームで特に発生しやすい。衝撃波面が格子に平行に配置されたとき、衝撃波面上に不規則な凹凸が成長してしまう。

対策はいくつかある。

1. HLL系スキームへの切り替え: HLLやHLLCは接触不連続面の解像度が低い代わりにCarbuncleフリー

2. H-correction（Sanders et al., 1998）: Roeスキームに横方向の数値粘性を追加

3. Rotated Roe: フラックス計算の座標系を局所的に回転させる

4. AUSM+up: 低マッハ数補正付きAUSMはCarbuncleに対してロバスト

それも理由の一つだよ。AUSM系は衝撃波安定性に優れるから、汎用ソルバーのデフォルトには適しているんだ。

Godunov系スキームは低マッハ数で過度な数値拡散を生じる。これはフラックスの散逸項が $O(1/M)$ でスケーリングするためだ。対策として以下がある。

• プリコンディショニング: Turkel, Choi-Merkleなどの手法で固有値を修正し、低マッハ数での数値拡散を抑制。FluentではLow Speed Preconditionerとして実装されている
• AUSM+up: 圧力項に低マッハ数補正を追加したスキーム。全速度域で使用可能
• 圧力ベースソルバーへの切り替え: $M < 0.3$ が支配的なら非圧縮性ソルバーの方が効率的。FluentのCoupled Pressure-Based Solverは$M < 2$程度まで対応

All-speed schemeは活発な研究分野で、STAR-CCM+のCoupled Flowソルバーはそれに近い設計思想だよ。非圧縮極限と圧縮性を一つのフレームワークで扱える点が強みだ。