ポンプCFD解析
ポンプ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も使える。
比速度という羅針盤——斜流ポンプ設計の出発点
ターボポンプ設計の最初の一歩は比速度（Specific Speed、Ns）の計算だ。Ns = n√Q / H^(3/4)（n：回転数、Q：流量、H：揚程）で定まるこの無次元数が、ポンプの形式（遠心・斜流・軸流）を決定する。斜流ポンプの設計領域（Ns = 400〜1200）では、理論的には1D設計（速度三角形）で基本寸法を決め、CFDで3D流れ詳細を検証するという2段階アプローチが定石だ。比速度の理論はアメリカの水力学者が20世紀初頭に体系化し、その後欧州・日本で独自の単位系による混乱が続いた。現代ではIEC規格によりm³/sとm単位での無次元比速度が国際標準化されている。
ポンプCFDの数値計算手法
メッシュ生成
遠心ポンプのメッシュで注意すべき点は？
インペラは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の数値的困難
斜流ポンプは比速度（Ns）が軸流ポンプと遠心ポンプの中間域（Ns = 400〜1200）に位置し、両方の流れ機構が混在するため数値解析が難しい。軸流成分では翼端渦と二次流れが支配的になり、遠心成分ではコリオリ力による流れの曲がりが効いてくる。この混在領域では乱流モデルの選択がポンプ効率予測に直接影響し、k-ω SSTモデルが通常のk-εより数%精度が良いことが多くの比較研究で示されている。メッシュ設計では軸方向・周方向・径方向それぞれの解像度バランスを取ることが重要で、一方向に細かくすれば良いわけではない。
ポンプCFDの実務適用
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
日本の地下放水路（首都圏外郭放水路など）や排水機場に使われる大型斜流ポンプは、突発的な大雨の際に毎秒200m³以上の流量を処理する能力が求められる。設計段階のCFD解析は全流量範囲（設計流量の50〜120%）にわたる効率・圧力特性の予測と、吸込み槽内の気泡巻き込みリスク評価が必須項目だ。特に吸込み渦の発生は実機で致命的なトラブルになるため、流域全体（吸込みピット形状を含む）をモデル化した自由表面CFD（VOF法）が納品前の必須検証となっている。モデルポンプ試験との相関は効率で±1%以内が国内メーカーの品質基準だ。
ポンプ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ツール選定の現実——汎用機vs専用機の攻防
斜流ポンプのCFD解析ツール選定では「ターボ機械専用機能の有無」が重要な評価軸だ。ANSYS TurboGridやNumeca AutoGrid5のような専用ブレードメッシュ生成ツールは、複雑な斜流翼形状でも数時間でHex-dominant高品質メッシュを生成できる。一方、汎用メッシャー（ANSYS Meshing、Pointwise）では同じ品質を得るのに熟練者でも1〜2日かかることがある。コスト面では専用ツールの年間ライセンスが数百万円になるため、年間解析案件数が少ない企業はOpenFOAMのturboBladeフォークを採用するケースが増えている。農業・治水系ポンプメーカーでは在来の1D設計ツール（Pump-Flo等）との連携を重視する傾向が強い。
ポンプCFDの先端研究
多段遠心ポンプ
多段ポンプの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で性能マップを再構築する必要がある。
斜流ポンプ最適化の最前線——位相解析で掴む不安定流れ
斜流ポンプの最先端研究は「回転失速（Rotating Stall）の早期検知」に集中している。回転失速は部分流量域で発生し、ポンプ特性曲線の折り返しを引き起こして運転不安定の原因になる。2010年代後半から固有直交分解（POD）や動的モード分解（DMD）をCFD結果に適用し、流れ場の支配モードを抽出する手法が論文発表されている。日本の農業用大型揚水ポンプ（揚程20〜50m、流量10m³/s超）のリニューアル設計では、このCFD + DMD解析が採用され、不安定流れの発生条件をそれまでの経験則より30%広い流量域まで回避できる設計を実現した。
ポンプCFDのトラブル対応
低流量での発散
低流量で計算が発散するのは何が原因ですか？
低流量ではインペラ入口で大きなプレスワール（逆旋回）が発生し、翼前縁で大規模な剥離が起きる。定常計算ではこの非定常現象を安定して捉えられない。
対策:
- 非定常計算に切り替える（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の最終チェック項目をまとめてください。
斜流ポンプのS字特性——設計者を悩ませる不安定曲線
斜流ポンプの特性曲線（H-Q曲線）に現れる「S字」は設計者の天敵だ。曲線が右肩上がりから折り返すS字形状を持つとき、複数の流量で同じ揚程が実現でき、並列運転時に不安定な流量配分が起きる。CFDでS字特性を捕捉するには、流量を設計点の20%程度まで絞った低流量域の解析が必要で、この領域ではRANSモデルの精度が著しく低下し大きな計算誤差が生じることがある。経験豊富なエンジニアは「低流量域はLES or DESで確認する」というルールを設けており、RANS結果だけを信じて設計を確定すると後工程で痛い目にあうことが業界で繰り返し起きている。
なった
詳しく
報告