SIMPLE法
理論と物理
SIMPLE法の概要
先生、SIMPLE法って名前はよく聞くんですけど、何の略で、何をする手法なんですか？
SIMPLE は Semi-Implicit Method for Pressure-Linked Equations の略で、1972年にPatankarとSpaldingが発表した圧力-速度連成アルゴリズムだよ。非圧縮性Navier-Stokes方程式を有限体積法で解くとき、圧力と速度をどうやって整合的に求めるかという問題を解決するんだ。
圧力と速度の連成って、なぜそんなに厄介なんですか？
非圧縮性流れでは連続の式（質量保存）が圧力の独立した方程式を含まないんだ。密度が一定だから状態方程式で圧力を求められない。結果として、運動量方程式と連続の式を同時に満たす圧力場と速度場を見つける必要がある。これが圧力-速度連成問題だよ。
支配方程式
まず基本の方程式を教えてください。
非圧縮性流れの支配方程式は次のとおりだ。
運動量方程式（Navier-Stokes方程式）:
ここで $\mathbf{u}$ は速度ベクトル、$p$ は圧力、$\rho$ は密度、$\nu$ は動粘性係数だ。
未知数が4つ（u, v, w, p）で方程式も4つだから、原理的には解けるはずですよね？
そう。ただし圧力のPoisson方程式を直接導出して同時に解くのは計算コストが高い。SIMPLE法は「予測-補正」というアプローチで、この連成を効率的に分離して解くんだ。
SIMPLE法のアルゴリズム
具体的な手順を教えてください。
SIMPLE法の1反復は次の4ステップで構成される。
Step 1: 運動量方程式を仮の圧力 $p^*$ で解く
これで仮速度 $\mathbf{u}^*$ を得る。この速度場は一般に連続の式を満たさない。
Step 2: 圧力補正方程式を解く
連続の式の残差を消すための圧力補正 $p'$ を求める:
Step 3: 圧力と速度を補正する
ここで $\alpha_p$ は圧力の緩和係数だ。
Step 4: その他のスカラー方程式を解き、収束判定を行う
$a_P$ っていうのは運動量方程式の対角係数ですよね。圧力補正方程式の右辺が速度の発散ということは、連続の式の残差をゼロに近づけてるんですね。
その通り。反復を繰り返すと $\nabla \cdot \mathbf{u}^* \to 0$ となり、質量保存が満たされるんだ。
緩和係数の役割
$\alpha_p$ の緩和係数ってどういう意味ですか？
SIMPLE法では圧力補正で隣接セル係数（$a_N$項）の寄与を省略している近似があるため、1回の補正では過補正になりがちなんだ。だから圧力には $\alpha_p = 0.3$ 程度、速度には $\alpha_u = 0.7$ 程度の亜緩和（under-relaxation）を適用するのが一般的だよ。経験則として $\alpha_u + \alpha_p \approx 1$ とする流儀もある。
なるほど、近似の代償として緩和が必要になるわけですね。これがSIMPLEの「Semi-Implicit」たる所以ですか。
そうだ。完全に陰的（Fully Implicit）に解くなら連成ソルバーになるが、SIMPLE法は分離型（Segregated）ソルバーとして非常に効率的なんだ。
Patankarが1970年代の計算機でSIMPLE法を動かした話
SIMPLE法（Semi-Implicit Method for Pressure-Linked Equations）を開発したS.V. Patankarは、1972年当時の計算機でNavier-Stokes方程式を解くという挑戦に挑んだ。当時のメモリは現代のスマホの1万分の1以下。そこで彼が取った戦略が「圧力と速度を分離して交互に解く」という反復法——これがSIMPLEの本質。完璧な連立方程式を一気に解かなくていい、という発想の転換は、限られた計算資源を最大限活用するための知恵だった。現代でも数百万セルのCFD解析にSIMPLEが使われているのは、この「分離型反復」のシンプルさと実績の賜物です。
各項の物理的意味
- 時間項 $\partial(\rho\phi)/\partial t$：蛇口をひねった瞬間を思い浮かべてください。最初は水がバタバタと不安定に出て、しばらくすると安定した流れになりますよね？ この「変化している最中」を記述するのが時間項です。心臓の拍動で血流が脈打つのも、エンジンのバルブが開閉するたびに流れが変動するのも、すべて非定常現象。では定常解析とは？ 「十分時間が経って流れが落ち着いた後」だけを見る——つまりこの項をゼロにする。計算コストが大幅に下がるため、まず定常で解いてみるのがCFDの基本戦略です。
- 対流項 $\nabla \cdot (\rho \mathbf{u} \phi)$：川に落ち葉を落としたらどうなりますか？ 流れに乗って下流に運ばれますよね。これが「対流」——流体の動きが物を運ぶ効果です。暖房の温風が部屋の端まで届くのも、空気という「運び屋」が熱を対流で輸送しているから。ここが面白いところ——この項は「速度×速度」を含むため非線形です。つまり、流れが速くなるとこの項が急激に強くなり、制御が難しくなる。これが乱流の根本原因です。よくある勘違い：「対流と伝導は同じようなもの」→ 全然違います！ 対流は流れが運ぶ、伝導は分子が伝える。桁違いの効率差があります。
- 拡散項 $\nabla \cdot (\Gamma \nabla \phi)$：コーヒーにミルクを入れて放置したことはありますか？ かき混ぜなくても、しばらく経つと自然に混ざりますよね。あれが分子拡散です。では次の質問——ハチミツとお水、どちらが流しやすいですか？ 当然お水ですよね。ハチミツは粘性（$\mu$）が高いから流れにくい。粘性が大きいと拡散項が強くなり、流体は「もったりした」動きになります。レイノルズ数が小さい流れ（ゆっくり、ドロドロ）では拡散が支配的。逆にRe数が大きい流れでは対流が圧倒し、拡散は脇役になります。
- 圧力項 $-\nabla p$：注射器のピストンを押すと、液体が針先から勢いよく出ますよね？ なぜでしょう？ ピストン側が高圧、針先が低圧——この圧力差が流体を押す力になるからです。ダムの放水も同じ原理。天気図で等圧線がギュッと密になっている場所では？ そう、強風が吹きます。「圧力差があるところに流れが生まれる」——これがナビエ-ストークス方程式の圧力項の物理的意味。ここでの勘違いポイント：CFDの「圧力」は絶対圧ではなくゲージ圧のことが多い。圧縮性解析に切り替えたとたんに結果がおかしくなる場合、絶対圧/ゲージ圧の混同が原因かもしれません。
- ソース項 $S_\phi$：暖められた空気が上に昇る——なぜでしょう？ 周囲より軽く（密度が低く）なったから、浮力で押し上げられるのです。この浮力はソース項として方程式に追加されます。他にも、ガスコンロの炎で化学反応熱が発生する、工場の電磁ポンプで金属溶湯にローレンツ力がかかる…これらはすべて「外部から流体にエネルギーや力を注入する」作用であり、ソース項で表現します。ソース項を忘れるとどうなるか？ 自然対流の解析で浮力を入れ忘れると、流体は一切動かない——冬の部屋で暖房をつけたのに暖かい空気が上に行かない、という物理的にありえない結果になります。
仮定条件と適用限界
- 連続体仮定：クヌッセン数 Kn < 0.01（分子平均自由行程 ≪ 代表長さ）で成立
- ニュートン流体仮定：せん断応力と歪み速度が線形関係（非ニュートン流体では粘度モデルが必要）
- 非圧縮性仮定（Ma < 0.3の場合）：密度を一定として扱う。マッハ数0.3以上では圧縮性効果を考慮
- ブシネスク近似（自然対流）：密度変化を浮力項のみで考慮し、他の項では一定密度を使用
- 適用外ケース：希薄気体（Kn > 0.1）、超音速・極超音速流れ（衝撃波捕捉が必要）、自由表面流れ（VOF/Level Set等が必要）
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 速度 $u$ | m/s | 入口条件で体積流量から換算する際、断面積の単位に注意 |
| 圧力 $p$ | Pa | ゲージ圧と絶対圧の区別。圧縮性解析では絶対圧を使用 |
| 密度 $\rho$ | kg/m³ | 空気: 約1.225 kg/m³@20°C、水: 約998 kg/m³@20°C |
| 粘性係数 $\mu$ | Pa·s | 動粘性係数 $\nu = \mu/\rho$ [m²/s] との混同に注意 |
| レイノルズ数 $Re$ | 無次元 | $Re = \rho u L / \mu$。層流/乱流遷移の判定指標 |
| CFL数 | 無次元 | $CFL = u \Delta t / \Delta x$。時間刻みの安定性に直結 |
数値解法と実装
SIMPLE族アルゴリズムの比較
SIMPLE以外にもSIMPLECとかSIMPLERとか聞くんですけど、何が違うんですか？
SIMPLE法を改良したバリエーションがいくつかあるんだ。主要なものを比較しよう。
SIMPLEC（SIMPLE-Consistent）
Van Doormaalと Raithby（1984）が提案。圧力補正方程式で隣接セル係数の影響を部分的に考慮する改良版だ。
$a_P$ の代わりに $a_P - \sum a_N$ を使うことで、圧力の緩和係数を1.0に近づけられる。収束が速くなる場合が多いんだ。
SIMPLER（SIMPLE-Revised）
Patankar（1980）が改良版として提案。圧力補正の前に圧力方程式を別途解いて、より良い圧力場を推定する手法だ。
SIMPLECだと緩和係数を大きくできるのは嬉しいですね。実務ではどれを使うのが一般的ですか？
Rhie-Chow補間
有限体積法でSIMPLEを実装するとき、圧力と速度の格子配置ってどうなるんですか？
これは非常に重要なポイントだ。コロケーション格子（圧力と速度を同じ位置に配置）を使うと、圧力のチェッカーボードパターンが発生する問題がある。
有限体積法でSIMPLEを実装するとき、圧力と速度の格子配置ってどうなるんですか？
これは非常に重要なポイントだ。コロケーション格子（圧力と速度を同じ位置に配置）を使うと、圧力のチェッカーボードパターンが発生する問題がある。
Rhie-Chow補間（1983）はこれを解決するための手法で、セル面上の速度を計算するとき圧力の3次導関数的な項を追加することで、チェッカーボード振動を抑制する。
現在のほとんどの商用CFDコードはコロケーション格子＋Rhie-Chow補間を採用しているんだ。
スタガード格子（速度と圧力を別の位置に配置）では不要なんですか？
スタガード格子ではチェッカーボード問題は自然に回避される。ただし非構造格子への拡張が難しいため、現代のCFDコードではコロケーション格子が主流なんだよ。
線形方程式系の解法
各ステップで線形方程式系を解くわけですが、どんな手法が使われますか？
運動量方程式にはガウス-ザイデル法やILU前処理付きBiCGSTAB法が一般的だ。圧力補正方程式は楕円型でスペクトルが広いため、AMG（代数的マルチグリッド）前処理が有効なんだ。
| 方程式 | 推奨手法 | 内部反復回数の目安 |
|---|---|---|
| 運動量（u, v, w） | ガウス-ザイデル / ILU-BiCGSTAB | 3〜10 |
| 圧力補正 | AMG前処理付きCG / GAMG | 50〜200 |
| スカラー（k, epsilon等） | ガウス-ザイデル | 3〜5 |
圧力補正の方が反復回数が多いんですね。
圧力補正方程式はラプラシアン型で情報の伝播が遅いんだ。だからマルチグリッドが効くし、ここがボトルネックになることが多い。全体の計算時間の60〜80%を圧力補正に費やすこともある。
SIMPLECとSIMPLER——「C」と「R」の違いって何だっけ？
SIMPLE法には派生型がいくつかある。SIMPLECは「Consistent」の頭文字で、速度補正係数の計算を改良して収束を速めた版。SIMPLERは「Revised」で、圧力の初期推定精度を上げた版。実際のところ、どっちが速いかはケースバイケース。単純な管内流れではSIMPLECがSIMPLEより30%速く収束した事例がある一方、複雑な再循環流れでは差が出ないこともある。現場では「まずSIMPLEで様子を見て、収束が遅ければSIMPLECに切り替える」という判断プロセスが多い。大切なのはアルゴリズムの名前より緩和係数のチューニングだという意見も根強い。
風上差分（Upwind）
1次風上: 数値拡散が大きいが安定。2次風上: 精度向上するが振動のリスク。高レイノルズ数流れでは必須。
中心差分（Central Differencing）
2次精度だが、Pe数 > 2で数値振動が発生。低レイノルズ数の拡散支配流れに適する。
TVDスキーム（MUSCL、QUICK等）
リミッタ関数により数値振動を抑制しつつ高精度を維持。衝撃波や急勾配の捕捉に有効。
有限体積法 vs 有限要素法
FVM: 保存則を自然に満足。CFDの主流。FEM: 複雑形状・マルチフィジックスに有利。SPH等のメッシュフリー法も発展中。
CFL条件（クーラン数）
陽解法: CFL ≤ 1が安定条件。陰解法: CFL > 1でも安定だが、精度と反復回数に影響。LES: CFL ≈ 1を推奨。物理的意味: 1タイムステップで情報が1セル以上進まないこと。
残差モニタリング
連続の式・運動量・エネルギーの各残差が3〜4桁低下で収束と判断。質量保存の残差は特に重要。
緩和係数
圧力: 0.2〜0.3、速度: 0.5〜0.7が一般的な初期値。発散する場合は緩和係数を下げる。収束後は上げて加速。
非定常計算の内部反復
各タイムステップ内で定常解に収束するまで反復。内部反復数: 5〜20回が目安。残差がタイムステップ間で変動する場合は時間刻みを見直す。
SIMPLE法のたとえ
SIMPLE法は「交互に調整する」手法。まず速度を仮に求め（予測ステップ）、その速度で質量保存が満たされるよう圧力を補正し（補正ステップ）、補正された圧力で速度を修正する——このキャッチボールを繰り返して正解に近づく。2人で棚を水平にする作業に似ている：片方が高さを合わせ、もう片方がバランスを取り、これを交互に繰り返す。
風上差分のたとえ
風上差分は「川の流れに立って上流の情報を重視する」手法。川の中にいる人が下流を見ても水の出所は分からない——上流の情報が下流を決めるという物理を反映した離散化手法。精度は1次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実務での設定手順
先生、実際にCFD解析でSIMPLE法を使うとき、どんな手順で進めればいいですか？
定常解析でSIMPLE/SIMPLECを使う典型的なフローを示そう。
1. メッシュ生成と品質確認: 非直交性 < 70度、スキューネス < 0.85 が目標
2. 離散化スキームの選択: まず1次風上で安定に収束させ、そこから2次精度に切り替える
3. 緩和係数の初期設定: 圧力 0.3、速度 0.7、乱流量 0.5〜0.7
4. 残差モニタリング: 全変数の残差が $10^{-4}$ 以下、かつ監視点の物理量が安定
5. メッシュ感度解析: 3水準以上のメッシュで結果の依存性を確認
ソフト別のSIMPLE設定
各ソフトでの具体的な設定方法を教えてください。
Ansys Fluent での設定
Fluentでは Pressure-Based Solver を選択した上で、Coupled ではなく SIMPLE や SIMPLEC を Pressure-Velocity Coupling で選ぶんだ。
- Methods > Pressure-Velocity Coupling > SIMPLE or SIMPLEC
- 圧力離散化: Standard, PRESTO!, Second Order
- 運動量: Second Order Upwind（安定したら）
- Pseudo Transient オプションで擬似時間進行による安定化も可能
OpenFOAM での設定
OpenFOAMの定常ソルバー simpleFoam では fvSolution の SIMPLE エントリで設定する。
```
SIMPLE
{
nNonOrthogonalCorrectors 1;
residualControl
{
p 1e-4;
U 1e-4;
}
}
relaxationFactors
{
fields { p 0.3; }
equations { U 0.7; k 0.5; epsilon 0.5; }
}
```
nNonOrthogonalCorrectors はメッシュの非直交性が高い場合に増やすとよいんだ。
STAR-CCM+ での設定
STAR-CCM+では Segregated Flow モデルを選択すると SIMPLE アルゴリズムが使われる。Under-Relaxation Factors タブで緩和係数を設定できるよ。
緩和係数のチューニング指針
緩和係数を最適に設定するコツってありますか？
いくつかの実践的なガイドラインをまとめよう。
| 状況 | 圧力緩和 | 速度緩和 | 備考 |
|---|---|---|---|
| 標準的な流れ | 0.3 | 0.7 | デフォルト |
| 高Re数の乱流 | 0.2 | 0.5 | 安定性重視 |
| 自然対流 | 0.3〜0.5 | 0.5〜0.7 | 浮力項との連成 |
| 多孔質媒体 | 0.1〜0.2 | 0.3〜0.5 | ソース項が大きい |
| 収束が安定したら | 0.5 | 0.8 | 加速して仕上げ |
残差が振動するなら緩和を小さく、収束が遅いなら少しずつ大きくする。自動緩和機能を持つソルバーもあるよ。
「まず安全な設定で回して、安定したら攻める」という戦略ですね。
そのとおり。CFDは「急がば回れ」が鉄則だ。
緩和係数0.7/0.3はどこから来た？——SIMPLE実践の経験則
SIMPLEを使うとき「速度の緩和係数は0.7、圧力は0.3」という設定をよく見かける。これ、実は確固たる理論的根拠があるわけじゃない。Patankarの教科書や大量のベンチマーク事例から積み上げられた経験則だ。設定を0.9にすると最初は速く進む感じがするが、複雑な流れでは途中でブレて発散する。逆に0.5以下にするとひたすら遅い。ある自動車エンジンルーム冷却解析で0.7/0.3がなかなか収束しないとき、速度だけ0.5に下げたらスムーズに収束した——という微調整の経験談はCFDエンジニアなら誰しも持っているもの。緩和係数は「魔法の定数」ではなく、問題に応じて感触で調整するものです。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
主要CFDソルバーにおけるSIMPLE実装
SIMPLE法の実装って、ソフトによって違いがあるんですか？
基本原理は同じだが、各ソルバーが独自の改良や拡張を加えているんだ。
| ソルバー | 分離型アルゴリズム | 連成型の有無 | 特記事項 |
|---|---|---|---|
| Ansys Fluent | SIMPLE, SIMPLEC, PISO, Coupled | あり | Pseudo Transient, High Order Term Relaxation |
| Ansys CFX | 独自連成型のみ | 標準 | 分離型は非搭載。結合型AMGで高効率 |
| STAR-CCM+ | SIMPLE系 | あり | AMG前処理、自動緩和機能 |
| OpenFOAM | SIMPLE, SIMPLEC | 別途 | simpleFoam, rhoSimpleFoam, buoyantSimpleFoam |
Ansys Fluent
Fluentが最も選択肢が多そうですね。
FluentはPressure-Based SolverとDensity-Based Solverの2系統を持っていて、Pressure-BasedでSIMPLE/SIMPLEC/Coupled Algorithmを選べる。Coupled Algorithmは分離型ではなく圧力と速度を同時に解く連成型で、SIMPLEより収束が速いが、メモリを1.5〜2倍使うんだ。
Ansys CFX
CFXは元々AEA Technology（英国）が開発したコードで、設計思想として最初から連成型（Coupled）ソルバーを採用している。SIMPLEのような分離型アルゴリズムは搭載していないのが特徴だ。全変数を一つの行列系で解くため、SIMPLE系のような緩和係数チューニングが不要という利点がある。
OpenFOAM
OpenFOAMだと自分でアルゴリズムの中身を確認できるんですよね？
そうだ。simpleFoam のソースコード（applications/solvers/incompressible/simpleFoam/）を読めば、SIMPLE法の実装が100行程度のC++で書かれているのが分かる。教育目的には最適だよ。fvSolution で線形ソルバーや前処理を細かく制御できる。圧力にはGAMG（Geometric-Algebraic Multi-Grid）がデフォルトで使われることが多い。
ライセンスとコスト
コスト面ではどうですか？
| ソルバー | ライセンス | 年間コスト目安 | サポート |
|---|---|---|---|
| Ansys Fluent | 商用 | 数百万円〜 | 公式サポートあり |
| Ansys CFX | 商用 | Fluent同梱パッケージあり | 公式サポートあり |
| STAR-CCM+ | 商用 | 数百万円〜 | Siemens公式 |
| OpenFOAM | GPL v3 | 無償 | コミュニティまたは有償 |
OpenFOAMは無償だが、GUIやメッシャーは別途用意する必要がある。HELYX、simFlowなどのGUIフロントエンドを使う選択肢もある。商用ソルバーはワークフローの完成度が高い分、ライセンスコストがかかるんだ。
FluentのSIMPLEとOpenFOAMのsimpleFoam——「同じ」と思ったら大間違い
SIMPLEは汎用アルゴリズムだが、各ツールの実装にはクセがある。FluentのSIMPLEは速度補正に「Rhie-Chow補間」を自動適用しているが、OpenFOAMのsimpleFoamはfvSolutionの設定で明示的にコントロールできる。同じパイプ流れをFluentとOpenFOAMで解いて圧力損失が10%違った——と困っているエンジニアの原因の多くは、この補間方法の違いと緩和係数のデフォルト値の差。さらにFluentはコントロールアップデートの頻度を内部で調整するが、OpenFOAMは毎ステップ忠実にアルゴリズムを回す。「どのSIMPLE？」という確認は、ベンチマーク比較で必ず行うべき作業です。
選定で最も重要な3つの問い
- 「何を解くか」：SIMPLE法に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
分離型ソルバーの限界と連成型への発展
SIMPLE法って50年以上前のアルゴリズムですよね。今も現役なんですか？
現役どころか、CFDの中核だよ。ただし、強い圧力-速度連成がある問題（高速回転流、強い浮力流れ）では分離型SIMPLEの収束が遅くなる。そこで連成型（Coupled）ソルバーが注目されている。
連成型では運動量方程式と圧力方程式をブロック行列として同時に解く:
AMG前処理と組み合わせることで、SIMPLE比で5〜10倍の収束加速が得られる場合がある。
GPU対応と大規模並列化
GPU計算との相性はどうですか？
SIMPLE法はステップが逐次的なので、GPU並列化の恩恵を受けにくい面がある。ただし各ステップ内の線形ソルバー（特にAMG）をGPUで加速することは可能だ。Ansys Fluent 2023R1以降ではGPUネイティブソルバーが搭載され、数億セル規模の計算がGPUクラスターで実行可能になっている。
OpenFOAMでもAmgXやPETScを介したGPU対応が進んでいる。
機械学習との融合
AIでSIMPLE法を高速化するような研究はありますか？
いくつかの方向性があるよ。
- 初期値推定: ニューラルネットで良い初期圧力場を予測し、SIMPLE反復回数を削減
- 緩和係数の自動調整: 強化学習で問題に応じた最適な緩和係数を逐次決定
- 前処理の学習: 線形ソルバーの前処理行列をMLで近似し、AMGを高速化
- PINN（Physics-Informed Neural Networks）: N-S方程式の物理拘束をネットワークに埋め込み、メッシュなしで解を近似
特にPINNは最近論文をよく見ますね。
PINNはまだ複雑な産業問題には精度が不足しているが、パラメトリックスタディの高速サロゲートモデルとして有望だ。従来のSIMPLE法と組み合わせて、粗い初期推定をPINNで得てからSIMPLEで仕上げるハイブリッド手法も研究されている。
マルチフィジックスへの拡張
SIMPLE法は単相非圧縮性流れだけでなく、以下の分野にも拡張されている:
- VOF法との結合: 自由表面流れ（interFoam等）
- 共役熱伝達: 流体-固体の温度連成（chtMultiRegionFoam）
- 反応流: 燃焼解析での化学種輸送との連成
- 圧縮性低マッハ数流れ: rhoSimpleFoam による可変密度SIMPLE
SIMPLE法がこんなに広い分野で使われているとは知りませんでした。基本アルゴリズムをしっかり理解しておくのが大事ですね。
SIMPLEは「遅い」から消えるのか——Coupled Solverとの競争
最近、大手商用CFDツールは「Coupled Solver（連成ソルバー）」を前面に押し出してきている。圧力と速度を分離せず一気に解くから理論上は収束が速い——じゃあSIMPLEは時代遅れ？ そうでもない。Coupled Solverはメモリ消費がSIMPLEの4〜5倍になることがある。1000万セルのメッシュでCoupled Solverを使おうとしたら64GBRAMが必要で、HPC申請が通らなかった——という笑えない話も聞く。SIMPLEはメモリ効率が良く、並列化との相性も実績がある。「低負荷で十分な精度」が求められる設計最適化の大量ケース解析では、今もSIMPLEが現役です。
トラブルシューティング
よくあるトラブルと対処法
先生、SIMPLE法の計算が収束しなくて困ってます。何からチェックすればいいですか？
SIMPLE法の収束問題は大きく分けて4パターンあるんだ。順番に見ていこう。
1. 残差が振動して収束しない
症状: 圧力や速度の残差がジグザグに振動する
原因と対策:
- 緩和係数が大き過ぎる → 圧力を0.2、速度を0.5に下げる
- メッシュの非直交性が高い →
nNonOrthogonalCorrectorsを2〜3に増やす（OpenFOAM） - 離散化スキームが高次すぎる → 一旦1次風上に落として安定させる
2. 残差が下がらずフラットになる
症状: 残差が $10^{-2}$ 程度で停滞する
原因と対策:
- メッシュ品質不良 → セルのスキューネス、アスペクト比を確認
- 境界条件の不整合 → 質量流量の入口/出口バランスを確認
- 乱流モデルの初期値が不適切 → $k$ と $\epsilon$ の初期値を流速から見積もり直す
3. 残差が発散する
症状: 数ステップで残差が急増し NaN になる
原因と対策:
- CFL数が大きすぎる（擬似時間進行法使用時）→ 擬似時間刻みを小さく
- 境界条件で逆流が発生 → 出口をPressure Outletにし、backflow条件を確認
- メッシュに負体積セルがある → メッシャーでリメッシュ
4. 連続の式の残差だけ下がらない
速度の残差は下がるのに圧力が下がらないってことですか？
そう。これはSIMPLE法特有の問題で、圧力補正方程式の解法に起因することが多い。
- AMG前処理のサイクル数を増やす（V-cycle → W-cycle）
- 圧力の線形ソルバー反復回数を増やす（50→200）
- メッシュの格子数比（coarsening ratio）を確認
ソルバー別の典型的エラー
各ソフトで出やすいエラーメッセージを教えてください。
Divergence detected in AMG solver: pressure correction→ 圧力緩和を下げる、メッシュ品質確認Reversed flow in N faces→ 出口境界を流れから十分離す
GAMG: Solving for p, Initial residual = 1, Final residual = 1→ 圧力方程式が解けていない。境界条件の型を確認（fixedValue vs zeroGradient）Floating point exception→ 負体積セル、ゼロ除算。checkMeshを実行
Diverging solution - Loss of fluid cells→ セル品質チェック、Time Step Factor を下げる
デバッグの手順
収束しないときの系統的なデバッグ手順をまとめよう。
1. メッシュ品質を確認: checkMesh（OpenFOAM）、Mesh Quality Report（Fluent）
2. 境界条件を簡略化: まず全壁面でテスト
3. 1次精度で安定計算: 収束してから2次精度に切り替え
4. 緩和係数を保守的に設定: 圧力0.2、速度0.5から開始
5. 残差の変数別推移を観察: どの変数が問題か特定
6. 物理量の場を可視化: 非物理的な領域がないか確認
9割のトラブルはメッシュか境界条件に起因する。ソルバー設定を疑うのは最後だよ。
SIMPLEが収束しない——「圧力参照点なし」という見落とし
SIMPLEで「残差がずっとプラトーのまま下がらない」というトラブルでよくある原因が「圧力基準点（pressure reference point）の未設定」。SIMPLE法は圧力の絶対値ではなく相対値（勾配）を解くので、全壁面で囲まれた閉じた系では圧力のレベルが不定になる。OpenFOAMで "singular matrix" と怒られたり、Fluentで残差が振動し続けたりするのはたいていコレ。解決策は圧力参照点をドメインのどこか一点に設定するだけ。でも出口が開放されていれば不要なため「そんな設定が必要だとは知らなかった」という初心者あるあるのトラブルです。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——SIMPLE法の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告