PISO法 — 先端応用とDNS/LESへの展開
PISO法のDNS/LESにおける役割
DNSやLESでPISO法はどのように使われていますか？
有限体積法ベースのDNS/LESでは、PISO法が圧力-速度連成の標準手法だ。特にOpenFOAMのdnsFoamやpisoFoamは学術研究で広く使われている。
DNSでは数値誤差が物理的なスケールに直接影響するため、PISO法の各ステップで非常に厳しい収束基準（残差 $10^{-8}$ 以下）が要求される。
分数ステップ法との関係
分数ステップ法（Fractional Step Method）とPISO法は関係があるんですか？
密接に関係している。分数ステップ法（Chorin, 1968; Temam, 1969）は:
1. 圧力項を除いた運動量方程式で仮速度を求める
2. 圧力ポアソン方程式を解く
3. 速度を補正して発散フリーにする
PISO法はこの分数ステップ法をFVM上で実装したものと見なせるんだ。FEM系のコード（例えばNek5000）では分数ステップ法と呼ばれ、FVM系（OpenFOAM, Fluent）ではPISO法と呼ばれることが多い。本質的には同じ考え方だよ。
高次精度FVMとの組み合わせ
空間精度をさらに上げる手法はありますか？
通常のFVMは2次精度が標準だが、以下の手法で高次精度化が可能だ:
- MUSCL再構築: 3次精度の勾配再構築。OpenFOAMのlinearUpwindが近い
- WENO（Weighted Essentially Non-Oscillatory）: 3〜5次精度。衝撃波を含む問題にも適用可能
- Compact Scheme: 高精度だが並列化が困難
LESでは空間フィルタ幅がメッシュサイズで決まるため、高次精度スキームを使うとフィルタのカットオフが急峻になり、解像可能な渦の範囲が広がるんだ。
スペクトル要素法との比較
FVMのPISO以外に、非定常CFDで使われる手法はありますか？
スペクトル要素法（SEM）が高精度DNS/LESで有力だ。代表的なコードにNek5000やNekRSがある。多項式次数Nを上げることで指数関数的に精度が向上する（p-refinement）。
ただしSEMは複雑形状のメッシュ生成が困難で、産業用途ではFVMのPISO/PIMPLE系が依然主流だ。最近はGPU対応のNekRS（Argonne国立研究所）が10億格子点規模のDNSを実現している。
今後の展望
PISO法の将来的な発展方向としては:
- GPU Native実装: PETScやAmgX経由ではなく、ソルバー全体をGPU最適化
- 適応メッシュ精密化（AMR）との統合: 渦構造に追従して動的にメッシュを細密化
- エクサスケール計算: 数兆セル規模のLES/DNSに向けた通信最適化
PISO法の基礎を理解しておけば、これらの先端技術も追いかけやすいですね。
レイノルズの実験（1883年）——乱流発見の瞬間
オズボーン・レイノルズは、管内の水にインクを流す実験で「層流から乱流への遷移」を発見しました。流速を上げていくと、インクの線がある瞬間にグチャグチャに乱れる。この劇的な瞬間を、レイノルズは数学的に $Re = \rho uD/\mu$ という無次元数で表現した。100年以上経った今も、CFDエンジニアが最初に確認するのはこのレイノルズ数です。
先端技術を直感的に理解する
この分野の進化のイメージ
CFDの最先端は「天気予報の進化」に似ている。かつての天気予報（RANS）は平均的な傾向しか分からなかったが、最新の数値天気予報（LES/DNS）は個々の雲の動きまでシミュレーションできる。AIとの融合により「数秒で近似予測」も可能になりつつある。
なぜ先端技術が必要なのか — PISO法の場合
従来手法でPISO法を解析すると、計算時間・精度・適用範囲に限界がある。例えば、設計パラメータを100通り試したい場合、従来手法では100回の解析が必要だが、サロゲートモデルを使えば数回の解析結果から100通りの予測が可能になる。「全部試す」から「賢く推測する」への転換が先端技術の本質。
CFDメッシュの品質管理や乱流モデルの選定に悩む時間を、もっと創造的な設計作業に使えたら。 — Project NovaSolverはそんな実務者の声から生まれました。
Project NovaSolver — CAE実務の課題に向き合う研究開発
「PISO法をもっと効率的に解析できないか？」——私たちは実務者の声に耳を傾け、既存ワークフローの改善を目指す次世代CAEプロジェクトに取り組んでいます。具体的な機能はまだ公開前ですが、開発の進捗をお届けします。
進捗通知を受け取る →