CFD解析の発散

数値的には、残差が単調に増加し、解が物理的に意味のない極端な値（例えば、局所的な圧力が10^30 Paなど）に飛んでいく現象です。物理的には、設定した条件が支配方程式（ナビエ-ストークス方程式）の解を安定に保つ範囲を超えていることを意味します。例えば、流速が速すぎて流れが剥離し、計算格子がその急激な変化を捉えきれなくなる状態です。

はい、重要な概念にCFL（Courant–Friedrichs–Lewy）条件があります。これは陽解法の時間刻み幅

理論的には陰解法は無条件安定と言われますが、実際の非線形問題ではそう単純ではありません。時間刻みが大きすぎると、非線形ソルバー（例えばニュートン法）の収束性が悪化し、結局発散に至ることがあります。Ansys Fluentのデフォルトの時間刻み設定は、現象の時間スケールを考慮して決める必要があります。例えば、シリンダー周りのカルマン渦を解く場合、ストローハル数から主流速度とシリンダー径で特徴時間を推定し、その1/20程度を初期の

初期計算や流れが複雑な領域では、高次精度スキームより「1次風上差分」から始めるのが鉄則です。これは数値的な拡散が大きく、解を安定させやすいためです。OpenFOAMでは`divSchemes`の設定で`Gauss linearUpwind`などが該当します。解が安定してきたら、2次精度の`Gauss linear`や`limitedLinear`に切り替えて精度を上げます。最初から`Gauss QUICK`のような高次スキームを使うと、発散のリスクが高まります。

大きな違いがあります。SIMPLE法は安定性が高い代わりに収束が遅く、定常解析の標準です。PISO法は非定常解析向きで、タイムステップあたりの修正回数を増やさないと発散しやすい。Ansys Fluentの「Coupled」ソルバーは、圧力と速度を同時に解くため、複雑な流れや浮力流れで強力ですが、メモリ消費が大きく、初期条件が悪いと一気に発散することもあります。実務では、複雑な内部流れにはまずSIMPLEで初期場を作り、その後Coupledに切り替える戦略も有効です。

発散しそうなときや発散した直後は、デフォルト値から下げます。例えばFluentでは、圧力の緩和係数はデフォルト0.3ですが、これを0.1や0.05まで下げます。運動量は0.7から0.3に下げます。これは更新量を小さくし、解の急激な変動を抑えるためです。ただし、収束は極端に遅くなるので、安定したら少しずつ元の値に戻していきます。逆に、収束が遅すぎるだけなら、緩和係数を上げる（ただし0.9が限度）ことで加速できます。

1) 計算を直ちに停止する。2) 最後に保存されたデータファイル（例えば、Fluentの`.dat`ファイル）から復元する。3) 「残差の履歴」と「領域内の最小/最大値（たとえば最小圧力）」を確認する。圧力が-10^5 Paなど物理的にあり得ない値になっているセルを特定します。4) そのセル周辺のメッシュ品質（歪み率、アスペクト比）をチェックします。これが最も一般的な初動です。

一般的なベストプラクティスとしては、歪み率（Skewness）は0.9未満（理想は0.5以下）、アスペクト比（Aspect Ratio）は100以下（境界層では1000を超えることもあるが注意）、直交性（Orthogonal Quality）は0.1以上（1に近いほど良い）を目安とします。Ansys MeshingやPointwiseではこれらの指標をレポートできます。特に、流速や圧力勾配が大きい領域のメッシュ品質は厳しくチェックします。

いくつかあります。第一に、先述の緩和係数を下げる。第二に、勾配の制限（Gradient Limiting）を有効にする。Fluentでは「Solution Limits」で最小絶対圧力（例：-1e5 Pa）や最大温度を設定し、解がその範囲を超えたらクリップできます。第三に、より頑健な離散化スキームに一時的に戻す。これらは根本解決ではありませんが、とりあえず計算を進めて問題領域を特定するのに役立ちます。

大きく異なります。FluentはGUI上で残差が赤く「Diverged」と表示され、最小/最大値の監視が容易です。ソルバーが自動的に「Turbulent Viscosity Limited」などの警告を出すことも。一方、OpenFOAMはコマンドライン出力で`time step continuity errors`が急激に増大（例えば1e+10など）したり、`solution singularity`というエラーで停止します。OpenFOAMは`controlDict`で`maxCo`（最大クーラン数）を設定し、これを超えると自動的に時間刻みを削減する機能があり、発散防止に効果的です。

強弱というより、アプローチが違います。COMSOLはデフォルトで「分離型ソルバー」を用い、各物理場を順番に解きます。非常に安定していますが、強い連成がある問題（例えば、化学反応を伴う高速流れ）では収束が難しいことがあります。その場合、「完全連成ソルバー」に切り替える必要があります。COMSOLの強みは、ソルバーの「定常化」設定で、自動的にダンピングをかけて発散を防ぐ試みをしてくれる点です。ただし、メッシュ依存性は他のソフトと同様にあります。

はい、「Implicit Unsteady」ソルバーにおける「Robustness Settings」が特徴的です。これには「Pressure Variation Damping」や「Momentum Damping」があり、発散しそうなときにソルバーが自動的にダンピング係数を増加させます。また、セル品質が悪い領域に対して自動的に数値的拡散を追加する「Accuracy/Stability Optimizer」機能もあります。初期のベンチマークでは、意図的に悪いメッシュを与えた場合、Star-CCM+は他のソフトより計算が落ちにくいという結果もありました。

これは非常に一般的な問題です。流入境界や壁面近傍で設定された乱流強度Iや乱流粘度比μ_t/μが大きすぎる、または初期場のkとεのバランスが取れていないのが原因です。具体的には、

界面の急激な変形は、体積分率の急勾配を生み、対流項の離散化が破綻しやすくなります。対策は、1) 界面付近のメッシュを十分に細かくする（少なくとも5〜10セルで界面を解像）。2) 界面捕捉スキームを、高次精度だが発散しやすい`Compressive`系から、まずは`Geometric Reconstruction`や`High Resolution Interface Capturing`（Fluent）など安定性重視のスキームで始める。3) 時間刻みを厳しく制御し、界面のクーラン数を0.25以下に保つ。OpenFOAMの`interFoam`では`alpha.water`の場に対する`nAlphaCorr`（反復回数）を増やすことも有効です。

原因は重なる部分もありますが、対処法が異なります。停滞は、ソルバーが局所的な解の「くぼみ」から抜け出せていない状態です。対策は、1) 初期推定値を変える（別の安定したケースから初期場を読み込む）。2) ソルバーを変える（SIMPLEからCoupledへ）。3) 離散化スキームの精度を一時的に下げて「滑らかに」する。一方、発散は解が制御不能に暴走している状態なので、緩和係数の低減や解の制限、メッシュの修正など、より強力な安定化策が必要です。停滞は「行き詰まり」、発散は「暴走」とイメージしてください。