遠心ポンプ（Centrifugal Pump） — CAE用語解説
遠心ポンプとは
遠心ポンプって、要するにインペラを回して水を送る機械ですよね？ CFDで何を解析するんですか？
ざっくり言うと、インペラ（羽根車）が回転して流体に運動エネルギーを与えて、それがボリュート（渦巻き室）やディフューザで圧力エネルギーに変換される機械だ。CFDでは主に「どの流量でどれだけの揚程が出るか」——つまりH-Q曲線（揚程-流量特性）を予測するのが一番の基本だね。
ポンプのCFDって難しそうですけど、普通の外部流れの解析と何が違うんですか？
最大の違いは「回転する部品がある」ことだ。インペラは数百〜数千rpmで回ってるから、固定されたケーシングと回転する羽根の両方を同時に扱う必要がある。これを処理するためにMRF（Multiple Reference Frame）法やスライディングメッシュ法を使うんだ。あとは乱流が強いから、乱流モデルの選択もシビアだよ。
定義と基本原理
遠心ポンプの揚程って、物理的にはどう決まるんですか？ 回転数を上げればいくらでも上がるんですか？
理論的にはオイラーのターボ機械方程式で決まる。インペラ出入口での周方向速度成分の変化が揚程を生むんだ：
ここで $u$ はインペラの周速、$c_\theta$ は流体の絶対速度の周方向成分、添字1が入口、2が出口だ。回転数を上げると $u$ が増えるから揚程は上がるけど、損失も増えるし、キャビテーションのリスクも高まる。実際には相似則で「回転数の2乗に比例して揚程が上がる」と近似できるよ。
$u_2 c_{\theta 2}$ って、結局は羽根の出口角度で決まるってことですか？
その通り。出口羽根角 $\beta_2$ が大きい（後退翼）ほど $c_{\theta 2}$ が小さくなって揚程は下がるけど、H-Q曲線の傾きが急になって運転安定性が増す。逆に前進翼（$\beta_2$ が小さい）だと揚程は高くなるけどサージしやすい。実務では後退翼（$\beta_2 = 20^\circ \sim 40^\circ$）が一般的だね。この羽根角の最適化こそCFDの出番だ。
H-Q特性曲線とCFD
H-Q曲線をCFDで出すって、具体的にはどういう手順なんですか？ 流量を変えて何回も計算するんですか？
まさにそう。典型的には5〜8点の流量条件（締切運転0%から過大流量150%くらいまで）を設定して、それぞれ定常MRF解析を回す。入口に質量流量を与えて、出入口の全圧差から揚程を算出するんだ：
例えば口径200mmの汎用ポンプ（1450rpm）なら、設計点で流量 $Q = 0.03\,\mathrm{m^3/s}$、揚程 $H = 25\,\mathrm{m}$ みたいな値になる。CFDで出したH-Q曲線を実験データと比較して、設計点で揚程誤差 $\pm 3\%$ 以内なら「よく合ってる」と判断するのが業界の感覚だよ。
効率も同時に出せるんですか？ ポンプの効率って結構重要ですよね。
もちろん。水動力 $P_w = \rho g Q H$ をインペラのトルク $T$ と角速度 $\omega$ から出した軸動力 $P_s = T\omega$ で割れば水力効率が出る：
実務ではCFDの効率は実験より2〜5ポイントくらい高く出ることが多い。ディスクフリクション（インペラ背面の摩擦損失）や漏れ流量がCFDで完全には捉えられないからだ。だから「CFDの効率をそのまま信じない」のが鉄則。傾向比較やパラメトリックスタディに使うのが正しいアプローチだね。
インペラCFD解析の回転モデリング
MRF法って、メッシュは動かさないのに回転の効果を出せるんですか？ ちょっとイメージが湧かないです。
MRFでは、インペラを囲む領域を回転座標系として定義して、Navier-Stokes方程式にコリオリ力と遠心力の項を追加するんだ：
ここで $\mathbf{v}_r$ は相対速度、$-2\boldsymbol{\omega} \times \mathbf{v}_r$ がコリオリ力、最後の項が遠心力だ。メッシュ自体は固定のまま「回転座標系で見た流れ」を解いてるから、定常解析で済む。計算コストが小さいのが最大のメリットで、設計の初期段階で何十ケースも回すときに重宝するよ。
じゃあスライディングメッシュはどういうときに必要なんですか？ MRFで十分なら使わなくていい気もしますけど。
MRFは定常近似だから、インペラとボリュートの干渉（ブレード通過周波数の圧力脈動）や、舌部（カットウォーター）近傍の非定常流れは捉えられない。具体的にこんなケースではスライディングメッシュが必須になる：
- 圧力脈動・振動予測 — ブレード枚数 $z$ × 回転数 $n$ の周波数（BPF: Blade Passing Frequency）の圧力変動を評価したいとき
- 騒音解析 — 流体音源の時間波形が必要なとき
- 部分負荷運転 — 設計点から大きく外れた運転条件で旋回失速が起きるとき
- 過渡応答 — 起動・停止時のウォーターハンマーやサージの解析
計算コストはMRFの10〜100倍だけど、インペラ1回転あたり最低360ステップ（1度刻み）は必要で、数回転分は回さないと周期定常に達しない。現実的には3〜5回転分を計算して、最後の1回転のデータを使うことが多いよ。
MRFとスライディングメッシュの境界面って、どうやってデータを受け渡すんですか？
インターフェース面で補間するんだ。MRFの場合は「凍結ローター法」とも呼ばれて、回転側と固定側の境界面で速度・圧力を座標変換して受け渡す。スライディングメッシュでは、毎タイムステップごとにメッシュが相対的にずれるから、非整合メッシュ間の補間（GGI: General Grid Interface）が使われる。ANSYS FluentだとMixing Planeという周方向平均を取る方法もあって、段間の受け渡しに使えるよ。大事なのは、このインターフェース面のメッシュ解像度を回転側と固定側で揃えること。片方が粗いと補間誤差で揚程が数%ずれることがある。
キャビテーション予測
キャビテーションってポンプでよく問題になりますよね。インペラの吸込み側で泡がバチバチ弾けて壊食が起きるやつ。CFDで予測できるんですか？
できる。キャビテーション解析では多相流モデルを使って、液相と蒸気相の共存を計算する。一般的なのは均質混合モデル（Homogeneous mixture model）で、Rayleigh-Plesset方程式に基づいて気泡の成長・崩壊速度を計算するSchnerr-SauerモデルやZwart-Gerber-Belamiモデルが広く使われてる：
ここで $R_B$ は気泡半径、$p_v$ は飽和蒸気圧、$p$ は局所圧力だ。$p < p_v$ になると気泡が成長して、圧力が回復する領域で崩壊する。
NPSHっていう指標をよく聞くんですけど、あれもCFDで出せるんですか？
NPSHr（必要吸込ヘッド、required NPSH）をCFDで求めるのは実務的にかなり有用だ。手順としては、入口圧力を段階的に下げていって、揚程が3%低下した時点のNPSHをNPSHrと定義する：
例えば25℃の水だと $p_v \approx 3.17\,\mathrm{kPa}$ だから、入口全圧を130kPa、120kPa、110kPa…と下げていって、揚程の落ち始めを捉える。ポイントは、キャビテーション解析ではメッシュをかなり細かくしないと気泡の発生位置が正しく予測できないこと。特にインペラ前縁（リーディングエッジ）周りは、最低でも羽根高さ方向に20層以上のメッシュが必要だよ。
キャビテーションが起きると、CFDの収束性にも影響しますか？
めちゃくちゃ影響する。気液界面が激しく変動するから、定常解析だと収束しないことが多い。実務では非定常解析（スライディングメッシュ+キャビテーションモデル）にせざるを得ないケースが増える。あと厄介なのは、蒸気相の体積率が急変する領域で圧力が振動すること。対策としてはタイムステップを十分小さく（インペラ回転0.5度/ステップ以下）して、アンダーリラクセーション係数を下げるのが定石だね。OpenFOAMだと interPhaseChangeFoam ソルバーで扱えるよ。
実務上のポイント
遠心ポンプのCFDで、メッシュはどう作ればいいんですか？ インペラの形状って複雑で難しそうです。
ポンプCFDのメッシュ作成は確かに難関の一つだ。いくつかのベストプラクティスがある：
- 羽根間流路の構造格子 — TurboGridやANSYS BladeModelerで羽根間に構造格子（H-O-H型トポロジー）を切るのが理想。非構造格子よりy⁺の制御がしやすい
- 壁面メッシュ — SST k-ωを使うなら $y^+ \approx 1$ を狙いたい。壁関数で済ますなら $30 < y^+ < 300$ だけど、インペラ翼面の剥離を正しく捉えたいなら壁関数は避けたほうがいい
- ボリュートとインペラの接続 — インターフェース面でのメッシュ密度を揃えること。周方向に360分割が一つの目安
- 総メッシュ数 — 単段ポンプなら200〜500万セルが実務的な範囲。羽根1枚あたり最低30万セルは欲しい
乱流モデルは何がおすすめですか？ k-εとSST k-ωどっちがいいんですか？
ポンプCFDではSST k-ωが事実上のスタンダードだ。理由は明確で、インペラ翼面の逆圧力勾配領域で起きる境界層剥離の予測が、標準k-εよりずっと正確だからだ。特に部分負荷（設計点の50〜70%流量）での再循環流れの捕捉にはSSTが圧倒的に有利。ただし、設計点近傍のH-Q曲線を出すだけなら、Realizable k-εでもそこそこ合う。もし圧力脈動や騒音まで踏み込むなら、DESやSAS（Scale-Adaptive Simulation）も検討してみて。
最後に、ポンプCFDで「ここだけは注意しろ」ってポイントがあれば教えてください。
3つ挙げるなら：
- 入口・出口の延長管 — インペラ直前で境界条件を与えると非物理的な流れになる。吸込口から直径の5倍以上、吐出口から10倍以上の延長管をモデルに含めること
- 隙間（チップクリアランス）のモデリング — インペラとケーシングの隙間から漏れ流量が発生する。これを無視すると揚程を5〜10%過大に予測する。ウェアリング隙間は最低5層以上のメッシュで解像すること
- 収束判定 — 残差だけでなく、必ずモニタリングポイント（入口圧力・出口圧力・トルク）の時間変化を確認すること。残差が下がっても揚程が振動してたら収束してない
この3つを押さえるだけで、ポンプCFDの精度はぐっと上がるよ。
関連用語
- MRF（Multiple Reference Frame） — 回転機械の定常近似。インペラ領域を回転座標系で解く
- スライディングメッシュ — メッシュを実際に動かす非定常解法。動静干渉を正確に捉える
- キャビテーション — 局所圧力が飽和蒸気圧以下になり気泡が発生する現象
- H-Q曲線 — 揚程（Head）と流量（Flow rate）の関係を示す性能曲線
- NPSH — Net Positive Suction Head。ポンプ吸込側のキャビテーション余裕度
- BPF（Blade Passing Frequency） — ブレード枚数 $\times$ 回転数の圧力脈動周波数
- 比速度 $n_s$ — ポンプの形式分類に用いる無次元パラメータ。$n_s = n\sqrt{Q}/H^{3/4}$
CAE用語の正確な理解は、チーム内のコミュニケーションの基盤です。 — Project NovaSolverは実務者の学習支援も視野に入れています。
遠心ポンプCFDの実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
お問い合わせ（準備中）関連トピック
なった
詳しく
報告