制限関数 — CAE用語解説

主に、数値計算における物理量の振る舞いを制限する関数です。例えば、CFDで密度や圧力が負の値になったり、衝撃波の近くで解が激しく振動するのを防ぎます。具体的には、流束制限関数（Flux Limiter）が有名で、高次精度スキームの単調性を保証します。

良い質問だ。1次風上差分は単調性を保つが、数値拡散が大きく解が鈍る。一方、中心差分など高次精度スキームは振動（オーバーシュート）を起こす。制限関数は、この両者の流束を混合し、滑らかな領域では高次精度を、変動の激しい領域では1次風上に近づける「ブレンド係数」の役割を果たします。混合の度合いを決めるのが制限関数

いくつかあるが、よく使われるのはMINMOD、Superbee、Van Leerあたりだ。MINMODは最も控えめで安定性が高いが、解が少し鈍る。Superbeeは鋭い勾配を保持するが、過剰な平坦化を招くリスクがある。Van Leerはその中間で、多くの商用ソフト（Ansys Fluentの2nd-order upwindのデフォルトなど）で採用されている実用的な選択肢だ。

対流項の空間離散化の段階で組み込まれる。MUSCL（Monotone Upstream-centered Schemes for Conservation Laws）スキームが典型例だ。セル界面での変数値

原理は同じだが、実装は複雑になる。非構造格子では、セル界面を挟む「上流」と「下流」のセルが必ずしも整列していない。そのため、界面の法線方向の勾配を評価し、その方向に沿った上流・下流のセル中心値の差から

大きく影響する。制限関数は非線形なので、陰的ソルバーで線形化するとヤコビアンが複雑になる。そのため、完全な陰的解法では収束が遅くなることがある。実務では、制限関数の部分を「凍結」（前の反復の値を使う）して線形化したり、Ansys Fluentのように「高次項の再構成」を反復の途中で更新する「スキームの部分緩和」といった工夫がなされている。

まず、ソフトウェアのデフォルトから始めるべきだ。Ansys Fluentなら「Second Order Upwind」がデフォルトでVan Leer制限関数を使っている。これで振動や負の密度が生じたら、より控えめなMINMOD系（Fluentでは「First Order」に近い挙動になる「Modified HRIC」など）を試す。逆に、衝撃波の位置をシャープに捉えたい航空宇宙分野では、一時的にSuperbeeを試すこともあるが、発散リスクが高いので最終結果には使わないことが多い。

重要なチェックポイントは二つ。第一に、物理量の雲図やコンター図で、極端な勾配の部分（衝撃波、火炎面、自由表面）に「階段状」の不自然な段差が生じていないか。これは過剰な制限（過度な拡散）のサインだ。第二に、残差履歴が振動していないか。特に密度やエネルギーなどの残差が周期的に跳ねる場合、制限関数による非線形性が強く効いている可能性がある。その場合は、より滑らかな制限関数に切り替えるか、メッシュを改良する。

「Solution Methods」ペインの「Spatial Discretization」セクションにある。例えば「Momentum」のスキームで「Second Order Upwind」を選んだ状態で、その下の「Limiter」プルダウンメニューから「Default」「Minmod」「Superbee」「Van Leer」「QUICK」などを選択できる。ただし、デフォルトはVan Leerベースの独自実装だ。多相流のVOF法では「Interface Modeling」の下に「Sharpening/Diffusion」の設定があり、ここでHRICやCICSAMといった界面用の特殊な制限スキームを選ぶ。

Star-CCM+では、「Continua」→「Models」→「Implicit Unsteady」や「Segregated Flow」の設定内にある「High-Order Scheme」のオプションで制限を選択する。デフォルトは「Venkatakrishnan Limiter」という、勾配の大きさに応じて滑らかに1次精度に落とす高度な関数を使う。OpenFOAMでは、`fvSchemes`ファイルの`gradSchemes`や`divSchemes`の設定で直接指定する。例えば、`div(phi,U) Gauss limitedLinearV 1.0;` の最後の`1.0`が制限関数の係数（1.0は通常の制限、0は制限なし）にあたる。

直接「制限関数」という用語はあまり使わないが、類似の概念は存在する。例えば、陽解法のLS-DYNAでは、要素のひずみや体積変化が極端に大きくなるのを防ぐ「体積粘性（bulk viscosity）」や「砂時計モード抑制」が、制限関数と同じ「発散防止」の役割を果たす。Abaqus/Explicitの「体積粘性スケーリング係数」も同様だ。これらは、衝撃波や大変形領域での数値的不安定性を制限・抑制するという点で、CFDの制限関数と哲学的には同じだ。

第一の対策として有効だ。特に乱流モデル（k-εなど）の変数は正の値であるべきなので、制限が弱いと負値が生じやすい。Ansys Fluentでこのエラーが出たら、「Solution Methods」で「Turbulence」の離散化スキームを「First Order Upwind」に一時的に落とす（これは強力な制限と同等）。収束し始めたら「Second Order Upwind」に戻し、その際に「Limiter」を「Minmod」などより強力なものに設定する。根本的には、エラーが発生する領域のメッシュを細かくするのが最も確実だ。

可能性が高い。メッシュを細かくしても衝撃波の位置が数セル以上ずれる場合、使用している制限関数が過度に拡散的（MINMODなど）で、衝撃波面を鈍らせている。逆に、Superbeeのような鋭い制限関数は、粗いメッシュでは勾配を過大評価し、衝撃波を早く出現させすぎる。対策は、メッシュを十分に細かくした上で、標準的な制限関数（Van LeerやVenkatakrishnan）を使い、衝撃波の前後で10〜20セル以上を確保することだ。超音速流れのベンチマーク（NASAの楔問題など）で設定を検証するのが良い。

あり得る。特に剥離流れや後流の予測では、制限関数による過剰な数値拡散が乱流構造を減衰させ、実質的な粘性を増大させてしまう。その結果、剥離領域が小さくなり、抗力が過小評価される。これを検証するには、同じメッシュで異なる制限関数（例：デフォルト vs MINMOD）で計算し、後流の渦構造や剥離線の位置を比較する。差が大きい場合は、メッシュの解像度不足が根本原因であることがほとんどだ。DESやLESのような乱流モデルでは、制限関数の影響はより顕著になるので注意が必要だ。