ポンプCFD解析
理論と物理
概要
遠心ポンプのCFD解析って何を予測するんですか？
揚程（Head）、効率、軸動力の3つが基本だ。H-Q特性曲線を作ることがCFDの主目的になる。
揚程と効率の定義
揚程の式を教えてください。
ポンプの揚程は入口と出口の全水頭差だ。
CFDでは全圧差から直接計算するのが簡便だ。$H = (p_{t2} - p_{t1})/(\rho g)$。
効率は水力効率と全効率に分かれる。
$\tau$ はCFDから得られるインペラのトルク、$\omega$ は角速度だ。
水力効率と全効率の違いは？
水力効率はディスク摩擦や漏れを含まない流体力学的損失のみ。全効率はディスク摩擦、漏れ流れ、機械損失を全て含む。CFDで直接得られるのは水力効率で、ディスク摩擦と漏れはウェアリングとの隙間モデルを含めないと出てこない。
Euler揚程（理論揚程）
Euler方程式から理論揚程が出せるんですよね？
遠心ポンプでは入口に旋回がない（$C_{\theta 1}=0$）と仮定すると $H_{Euler} = U_2 C_{\theta 2}/g$ になる。スリップファクタ $\sigma_s$ を考慮した揚程は $H_{th} = \sigma_s \cdot H_{Euler}$ だ。CFDの揚程はこの理論揚程に水力損失を加えたものに対応する。
MRF法による定常解析
ポンプのCFDはMRF法が一般的ですか？
H-Q曲線の取得にはMRF法（定常）が標準だ。ボリュート付きの場合はFrozen RotorまたはSliding Meshを使う。ボリュートなしのガイドベーン付きポンプならMixing Planeも使える。
遠心ポンプ理論の礎——オイラーのポンプ方程式（1754年）と現代との繋がり
遠心ポンプの基本理論であるオイラーの「ターボ機械方程式」（U1*Vtheta1 - U2*Vtheta2 = g*H）は、Leonard Eulerが1754年に導いたものだ。回転体が流体に与えるエネルギーは入口と出口での周速と旋回速度成分の差に比例するという原理は、270年後の現代でもポンプ・タービン・コンプレッサの教科書の第1章に登場する不変の真理だ。Eulerは流体機械だけでなく剛体力学・数学・光学で巨大な業績を残した18世紀最大の数学者だ。驚くべきことに、現代CFDのNavier-Stokes方程式の元祖「オイラー方程式（粘性なし）」も彼の名を冠している——ターボ機械の基礎理論とCFDの数値解法の両方にEulerの名が刻まれているのは、流体力学がいかに彼の業績の上に成り立っているかを示す証だ。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
メッシュ生成
遠心ポンプのメッシュで注意すべき点は？
インペラはTurboGridで構造格子を生成するのが最も高品質だ。ボリュートは非構造のテトラ/ポリヘドラルメッシュを使う。
遠心ポンプのメッシュで注意すべき点は？
インペラはTurboGridで構造格子を生成するのが最も高品質だ。ボリュートは非構造のテトラ/ポリヘドラルメッシュを使う。
| 領域 | メッシュタイプ | セル数目安 | ツール |
|---|---|---|---|
| インペラ | 構造格子（H/J/L+O-grid） | 50～150万/ピッチ | TurboGrid |
| ボリュート | 非構造テトラ+プリズム | 100～300万 | Fluent Meshing, STAR-CCM+ |
| 吸込管 | 構造 or 非構造 | 20～50万 | 任意 |
| ウェアリング隙間 | 構造（六面体） | 10～30万 | 手動 |
ウェアリング隙間もモデルに入れるんですか？
漏れ流れの影響を評価するなら必須だ。隙間は0.2～0.5mmと非常に狭いから、径方向に最低10セル、軸方向に50セル以上が推奨される。
乱流モデルの選択
ポンプにはどの乱流モデルが適していますか？
SST k-omegaが定番だ。翼面の逆圧力勾配と剥離予測に優れている。ポンプの場合、翼枚数が少なく（5～7枚）翼負荷が高いため、k-epsilonでは剥離の予測が甘くなりがちだ。
壁関数と Low-Re のどちらを使うべきですか？
境界条件
ポンプの典型的な境界条件は？
ポンプの典型的な境界条件は？
出口を静圧固定にして入口流量を変えるのが最も安定する設定だ。
ポンプCFDの実践的設定——回転座標系とMRFの使い分けとメッシュ要件
遠心ポンプのCFD解析でまず決める設定は「MRF（移動参照系）vs スライディングメッシュ（SM）」だ。設計点付近の全揚程・効率の基本特性確認ならMRF（定常）が計算時間の点で圧倒的に有利で、スライディングメッシュの1/10〜1/50のコストで同等の精度が得られる。羽根通過周波数（BPF）の振動・騒音・疲労が問題になるポンプでは、非定常SMが必要だ。メッシュ要件はインペラ流路に翼面法線方向20〜30層のプリズム層（y+<1で低Re壁面処理使用時）と、流路内のセル数が最低50万/翼通路が実務上の目安。ボリュート（Volute）やスクロールも含めたフルモデル解析では全体セル数が300〜500万になるのが典型的だ。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
H-Q特性の計算手順
H-Q特性曲線はどうやって作りますか？
1. 設計流量 $Q_d$ で定常MRF(or Frozen Rotor)計算を収束
2. 流量を $0.2Q_d$ ～ $1.4Q_d$ の範囲で7～10点設定
3. 各点で入口質量流量を変更して再計算（前の点の結果をリスタート値に使用）
4. 各点で揚程H、トルクτ、効率ηを算出してプロット
低流量側で収束が悪くなるのはなぜですか？
低流量ではインペラ入口のインシデンス角が大きくなり翼面で大規模剥離が発生する。定常計算では非定常な剥離構造を一つの解に収束させようとして振動するんだ。0.3$Q_d$以下では非定常計算が必要になることが多い。
実験との比較
CFDの結果はポンプ試験とどのくらい合いますか？
BEPって何ですか？
Best Efficiency Point、最高効率点だ。ポンプの設計点はBEP付近に設定される。JIS B 8301やISO 9906に基づくポンプ試験のデータと比較する。
ディスク摩擦の評価
ディスク摩擦損失もCFDで計算できますか？
インペラの背面（シュラウド側・ハブ側）の隙間空間もモデルに含めれば計算される。含めない場合は経験式で補正する。
$C_M$ は摩擦モーメント係数で、隙間比 s/r とレイノルズ数の関数だ。
消防ポンプのCFD性能認証——日本消防検定協会の基準とCFD適用
消防用ポンプは人命に関わる装備のため、厳格な性能基準（吐出量・全揚程・効率）が法令（消防法施行規則）で定められており、国内では消防検定協会の型式認定が必要だ。CFDはこの認定取得プロセスの開発段階で、試作前の性能予測・設計修正コスト削減に活用されている。ある消防ポンプメーカーでは、CFD（RANS+スライディングメッシュ）で揚程と流量の関係（H-Q曲線）を事前評価し、4回の形状変更によって目標特性を達成した後に初めて試作した結果、従来比で開発期間を30%短縮した事例が紹介されている。CFD結果と実機試験の差は全揚程で±3%以内に収まっており、設計ツールとしての信頼性が確立されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
ソルバー比較
ポンプのCFDに最適なソルバーはどれですか？
CFXが一番実績あるんですね。
ポンプ業界ではCFXが事実上の標準だ。結合型ソルバーの安定性がポンプの低流量運転や二相流計算に向いている。Fluent は2024年以降のTurbo Workflowの強化で急速に追い上げている。
STAR-CCM+はどうですか？
ポリヘドラルメッシュの自動生成が魅力で、ボリュートやケーシングの複雑な形状に対してメッシュ工数が少ない。ただし翼間流路の構造格子ほどの精度は出ないから、TurboGridとの併用が理想的だ。
計算コストの見積もり
遠心ポンプ1台分の解析にどのくらい時間がかかりますか？
| 解析内容 | セル数 | コア数 | 所要時間 |
|---|---|---|---|
| 定常MRF（1運転点） | 200万 | 16 | 2～4時間 |
| H-Q特性（8運転点） | 200万 | 16 | 16～32時間 |
| 非定常Sliding Mesh（1運転点） | 200万 | 32 | 12～24時間 |
| キャビテーション評価 | 300万 | 64 | 24～48時間 |
H-Q特性なら1～2日で出せるんですね。
自動化スクリプトで流量点を順次計算すれば、夜間バッチで翌朝には結果が揃う。これがCFDの大きな利点だ。
ポンプCFDツール比較——ANSYS TurboSystemとSimericsの特徴と選定基準
ポンプCFD解析の商用ツールを比較すると、ANSYS TurboSystem（TurboGrid + CFX連携）は翼型設計と高品質構造格子生成のワンストップ環境が強みで、多段ポンプや混流ポンプの精度検証データが豊富だ。SimericsはPumpLinxから発展したポンプ特化ツールで、カルテシアンベースの自動メッシュ生成で形状変更時の再設定時間が最小化される。設計最適化ループを自動化する点でSimericsは商業製品の段階的設計に向く。OpenFOAMはpimpleFoam + AMIでポンプ解析が可能で、研究機関や中小企業のコスト制限下での活用が多い。2024年時点でANSYS OptiSLangとTurboSystemの連携最適化機能が強化されており、形状パラメータ感度解析の自動化でCFD最適化ループの工数が従来の1/5以下になったと報告されている。
選定で最も重要な3つの問い
- 「何を解くか」：ポンプCFD解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
多段遠心ポンプ
多段ポンプのCFDは難しいですか？
各段をMixing PlaneまたはFrozen Rotorで接続する。段間のリターンチャネル（逆ガイドベーン）の損失予測が精度を左右する。
リターンチャネルって何ですか？
1段目の出口旋回を除去して2段目の入口条件を整える流路だ。ここでの損失は全段効率の2～5%を占める。メッシュを粗くしがちな領域だが、丁寧にモデル化すべきだ。
スラリーポンプのCFD
固体粒子を含むスラリーのポンプ解析もCFDでできますか？
Eulerian-Eulerian法やEulerian-Lagrangian法で固液二相流を計算する。
| 手法 | 適用 | ソルバー |
|---|---|---|
| Eulerian-Eulerian | 高濃度スラリー（体積分率>10%） | CFX Multiphase, Fluent Mixture |
| Eulerian-Lagrangian (DPM) | 低濃度（<5%）、粒子軌跡評価 | Fluent DPM, STAR-CCM+ Lagrangian |
| DEM連成 | 大粒径、粒子間接触が重要 | STAR-CCM+ DEM, EDEM連成 |
摩耗予測もできますか？
Finnie やOkaの摩耗モデルで翼面の摩耗量分布を推定できる。粒子衝突速度と角度から摩耗率を計算する。FluentのDPM Erosion ModelやSTAR-CCM+のErosion Modelが利用可能だ。
ポンプ as タービン (PAT)
ポンプを逆転させてタービンとして使うPATもCFDで評価できますか？
できる。小水力発電で注目されるPATは、ポンプを逆流・逆回転で運転する。CFDでは回転方向と流れ方向を反転させた解析を行う。BEPがポンプモードとは異なる点に移動するから、CFDで性能マップを再構築する必要がある。
ポンプ空洞現象（キャビテーション）予測の最前線——気泡動力学とCFDの融合
ポンプのキャビテーション解析では、Rayleigh-Plesset方程式（気泡の膨張・収縮ダイナミクス）とCFDを結合した「フルキャビテーションモデル」が最先端だ。Schnerr-Sauerモデルはこの方程式を単純化した気泡数密度ベースの輸送方程式で、ANSYSのFluentやCFXに標準実装されている。気泡崩壊時の局所圧力ピーク（数GPa）は通常のRANS格子では解像不可能で、キャビテーション侵食のリスク評価には「Collapse Energy（崩壊時のエネルギー密度）」を空間積分したErosion Risk Indicator（ERI）の使用が有効だ。最先端研究ではボイド率の実験計測（X線CT）とCFDの高精度比較によるモデル改良が行われており、KAPLAN水車の吸い込み側キャビテーション予測精度の向上が世界中の水力機械研究者の注目を集めている。
トラブルシューティング
低流量での発散
低流量で計算が発散するのは何が原因ですか？
低流量ではインペラ入口で大きなプレスワール（逆旋回）が発生し、翼前縁で大規模な剥離が起きる。定常計算ではこの非定常現象を安定して捉えられない。
対策:
- 非定常計算に切り替える（Sliding Mesh）
- タイムステップを翼通過時間の1/30以下に設定
- 緩和係数を下げる（CFXのTimescale Factor: 0.5）
逆流境界条件
入口に逆流が出る場合はどうしますか？
低流量でインペラ入口のハブ付近に逆流が発生することがある。入口BCを通常のInletのままにすると逆流セルで非物理的な圧力が発生して発散する。
対策:
- 入口境界を上流に十分延長（パイプ直径の5倍以上）して、逆流が境界面に到達しないようにする
- CFXのOpening BCに変更
- Fluentの Pressure Inlet with backflow direction を設定
ボリュートとの干渉
ボリュートのカットオフ付近で圧力が振動します。
カットオフ（舌部）はポンプ最大の非定常力の発生源だ。Frozen Rotorだと翼位置固定で過大な干渉が出る。Sliding Meshに切り替えると時間平均値が安定する。
チェックリスト
ポンプCFDの最終チェック項目をまとめてください。
ポンプCFDで低流量域の全揚程が実測より大きい——リサーキュレーションの見落とし
遠心ポンプのCFDで低流量域（設計流量の50%以下）に全揚程を予測すると実測より5〜10%高くなる系統誤差の原因の多くは「入口リサーキュレーション（Inlet Recirculation）」の不適切なモデル化だ。低流量域ではポンプ吸い込み口で一部の流体が逆流して渦を形成し、実効的な流れ面積が減少する。この現象をCFD（定常MRF）で正確に捉えるには、吸い込み管を十分な長さ（最低5D）含めたモデルが必要で、吸い込み口を単純な均一入口条件にすると逆流が抑制されて全揚程が過大評価される。また低流量域では本質的に非定常性が高く（失速セル、旋回失速）、定常解析では物理的に解が存在しない条件が含まれる——非定常URANSに切り替えることが根本的な解決策になる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ポンプCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告