比速度と設計指針
理論と物理
概要
比速度ってターボ機械の型式選定に使うパラメータですよね？
そう。比速度は流量、揚程（圧力差）、回転数から計算される無次元パラメータで、最適な機械型式を決定するための最も基本的な指標だ。
比速度の定義
式を教えてください。
日本で広く使われるのはこの形だ。
$N$:回転数[rpm]、$Q$:流量[$m^3/min$]、$H$:揚程[m]。無次元化されていないが、実用上はこの有次元形が最も使われる。
国際的には無次元比速度（$\Omega_s$ or $\omega_s$）も使われる。
比速度の値でどう型式が決まるんですか？
| 比速度 $N_s$ (有次元) | 型式 | 特徴 |
|---|---|---|
| 100～300 | 遠心（ラジアル） | 高揚程・低流量 |
| 300～600 | 斜流（ミックスフロー） | 中揚程・中流量 |
| 600～1500 | 軸流 | 低揚程・大流量 |
Cordierダイアグラム
Cordierダイアグラムって何ですか？
無次元比速度 $\omega_s$ と無次元比直径 $\delta_s$ の最適線を示す図だ。
最適効率を実現する $\omega_s$-$\delta_s$ の関係がCordier線として描かれており、ポンプでもタービンでも成立する。設計の出発点として、まずCordier線上に乗るかどうかを確認する。
CFDを始める前にCordierダイアグラムで確認するんですね。
そう。比速度とCordier線から基本寸法（直径、翼幅）の初期値を決めてからCFDに進むのが効率的なワークフローだ。
比速度の概念の歴史——19世紀の水力工学者が生んだ設計類似則
比速度（Specific Speed、Ns）の概念は19世紀後半の水力機械工学者たちが類似の水車を比較するために経験的に発展させた。イタリアの水力工学者や独のVogelが独立に類似のパラメータを提案し、1910〜20年代に欧米で統一的な定義が整備された。当初は「水力比速度」として水車設計に使われ、後に遠心ポンプ・コンプレッサにも拡張された。無次元形式（Shape Number Omega = omega*sqrt(Q)/(g*H)^(3/4)）は1940年代以降に物理的意味を明確にする形で整備されている。日本の水力機械工学では今もNsの国内計算式（Ns = N*sqrt(Q)/H^(3/4)）が使われており、ISOの無次元式とは係数が異なるため国際論文との比較で混乱が生じやすい——使用している定義を明示することが実務の基本だ。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
1D設計とCFDの連携
比速度で型式を決めた後、CFDまでの設計フローはどうなりますか？
典型的な設計フローを示そう。
1. 仕様決定: 流量Q、揚程H（or圧力比）、回転数N
2. 比速度計算: $N_s$ を算出し型式を選定
3. Mean-Line設計: 速度三角形、翼角、子午面寸法を決定
4. Throughflow解析: 2Dの子午面流れ場を計算
5. 3D翼形状定義: BladeGen等で3D翼面を生成
6. 3D CFD: TurboGrid + CFXで詳細解析
7. 最適化: パラメトリックスタディまたは自動最適化
Mean-Line設計って何をするんですか？
スパン中央（Mid-Span）の速度三角形を設計する。入口/出口の絶対速度、相対速度、翼周速から翼角を決め、de Haller数やディフュージョンファクタで負荷の妥当性を確認する。
Throughflow解析
Throughflow解析って何ですか？
子午面（r-z平面）上で流れを軸対称に仮定して解く2D解析だ。翼列の影響をブレードフォースとして体積力で模擬する。スパン方向の速度・圧力分布が得られるので、翼根から翼端までの翼角分布を決める基礎になる。
| ツール | 手法 | 開発元 |
|---|---|---|
| Concepts NREC COMPAL/AXIAL | Streamline Curvature法 | Concepts NREC |
| AxSTREAM (SoftInWay) | Streamline Curvature法 | SoftInWay |
| NUMECA AutoBlade | 翼形状自動定義 | NUMECA |
| Vista CCD (Ansys) | 1D遠心設計 | Ansys |
Streamline Curvature法って何ですか？
子午面の流線曲率に基づいてスパン方向の圧力勾配を計算する手法で、ターボ機械のThroughflow解析の標準手法だ。数秒で結果が得られるから、パラメトリック設計に最適だ。
比速度とCFD結果の照合——η-Ns曲線の解析的検証と数値計算精度の関係
比速度Ns（Specific Speed）はターボ機械を流量・揚程・回転数で無次元化した設計類似パラメータで、Ns値から最適な翼形タイプが決まる。CFD解析後の効率は実験値と比較前に、まずNs-η（比速度-効率）の既知チャートと照合することが重要だ。Lomakin(1958)やKaplan(1935)等の実験集成から作られたNs-η最大効率曲線から外れていれば、設計問題かCFDモデルの問題がある。特に「CFD効率がNs-η曲線の上限を超えている（物理的に不可能な高効率）」場合は、損失のモデル化漏れ（ボリュート損失、軸受損失等）か境界条件の誤設定がほぼ確実だ。
風上差分（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留意点
比速度によってCFDのアプローチは変わりますか？
大きく変わる。型式ごとの特徴をまとめよう。
| 比速度範囲 | 型式 | CFDの主要課題 | 推奨メッシュ手法 |
|---|---|---|---|
| 低 (100-200) | 遠心（低流量） | ディフューザ損失、再循環 | TurboGrid + ボリュート非構造 |
| 中低 (200-400) | 遠心（標準） | ジェット/ウェイク構造 | TurboGrid |
| 中 (400-600) | 斜流 | 子午面曲率の影響 | TurboGrid (軸-径方向混合) |
| 高 (600-1000) | 軸流（ハブ比大） | チップ漏れ、二次流れ | TurboGrid |
| 超高 (>1000) | 軸流（ハブ比小） | 非圧縮流れ、失速 | TurboGrid or 非構造 |
斜流ポンプのCFDは特に難しいですか？
子午面の曲率が大きいため、翼間流路のメッシュが歪みやすい。TurboGridでJ型やL型トポロジを選択し、子午面曲率に追従させる設定が重要だ。
比速度と効率の関係
比速度によって到達可能な効率は変わりますか？
比速度と最高効率には明確な関係がある。
この経験式はAndersonの式と呼ばれ、$N_{s,opt}$付近で最高効率になる。CFDの結果がこの傾向から大きくずれる場合は、設計またはCFDの設定に問題がある可能性が高い。
CFDの結果を経験式と比較するのは良い検証方法ですね。
比速度チャートやCordierダイアグラムとの整合性確認は、CFD結果の sanity check として非常に有効だ。
下水処理場のポンプ比速度選定——揚程条件変化への対応とCFD活用
下水処理場の汚水ポンプは流入量が季節・降雨で大きく変動するため、広い流量範囲で安定した揚程特性が必要だ。比速度Nsを用いた設計では、Ns=100〜200程度の「混流ポンプ」が流量変動に対してH-Q曲線が安定したフラット特性を示すため選ばれることが多い。一方Ns>300の軸流ポンプは流量変化に対してH-Q曲線が急峻で、低流量域で不安定（サドル型）になる傾向がある。実際の設計では、最大流量・最小流量・設計流量の3点でCFD（RANS+スライディングメッシュ）を実施し、各条件での吸い込み振動・渦発生・キャビテーション余裕（NPSHavailable）を確認する手順が国内の下水道施設設計指針に記載されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
予備設計ツールの比較
1D/2D予備設計に使えるツールは何がありますか？
| ツール | 対象 | 手法 | CFD連携 |
|---|---|---|---|
| Concepts NREC COMPAL | 遠心圧縮機 | Mean-Line | BladeGen→TurboGrid |
| Concepts NREC PUMPAL | 遠心ポンプ | Mean-Line | BladeGen→TurboGrid |
| Concepts NREC AXIAL | 軸流機 | Streamline Curvature | BladeGen→TurboGrid |
| AxSTREAM (SoftInWay) | 全型式 | 1D+Throughflow | AxSTREAM→CFD各社 |
| Vista CCD (Ansys) | 遠心圧縮機 | 1D | BladeGen→TurboGrid直接 |
| Vista AFD (Ansys) | 軸流ファン | 1D | BladeGen→TurboGrid直接 |
| NUMECA AutoBlade | 全型式 | パラメトリック翼形状 | AutoGrid5→FINE/Turbo直接 |
AnsysのVista CCD/AFDはどういうツールですか？
Ansys Workbench内で1D設計からCFDまでシームレスに連携できる。Vista CCDで遠心圧縮機の子午面とインペラ寸法を決めると、そのままBladeGenで3D翼面を定義し、TurboGridでメッシュ生成、CFXで計算という一気通貫のワークフローだ。
設計最適化フロー
予備設計からCFD最適化まで一貫して行えるんですか？
Ansys Workbenchのパラメータ連携とoptiSLangの組み合わせで可能だ。1D設計パラメータ（翼角、子午面形状等）をDesign of Experiments（DoE）で振り、CFDの効率を応答面で近似して最適解を探索する。
何回くらいCFDを回す必要がありますか？
設計変数の数によるが、5～10変数なら50～200回のCFDが目安だ。Latin Hypercube Samplingで初期点を生成し、Krigingサロゲートモデルで応答面を構築する。
F1と空力の戦い
F1マシンは時速300kmで走ると、車重と同じくらいのダウンフォース（下向きの空力的な力）を発生します。つまり理論上、天井に貼り付けて走れる！ チームは数千CPU時間のCFDシミュレーションを毎週実行し、フロントウィングの角度を0.1°単位で最適化しています。F1はCAEの技術力がそのまま順位に直結する世界です。
選定で最も重要な3つの問い
- 「何を解くか」：比速度と設計指針に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
マルチポイント設計
設計点だけでなくオフデザインも考慮した設計は可能ですか？
多目的最適化でマルチポイント設計を行う。例えば以下の3点を同時に最適化する。
- 設計点: BEPでの最高効率
- 低流量点: 0.7$Q_d$での効率（部分負荷特性）
- 高流量点: 1.2$Q_d$での安定性（サージマージン）
パレートフロント上のトレードオフ解から最終設計を選択する。
パレートフロントって何ですか？
複数目的関数の間でトレードオフの境界を示す解の集合だ。設計点効率を上げると部分負荷特性が悪化する、といった関係がパレートフロントに現れる。設計者がトレードオフを見ながら最終判断する。
可変速運転
ポンプやファンの可変速運転はCFDでどう評価しますか？
インバータ駆動で回転数を変える場合、各回転数でのCFDを行うか、ファン法則でスケーリングする。
Re数が大きく変わる場合（回転数比2:1以上）はスケーリングの誤差が大きくなるから、CFDで直接計算するほうがよい。
AI/MLによる設計探索
機械学習を使った設計探索は実用化されていますか？
研究段階から実用段階に移りつつある。CFDの計算結果を学習データにしてサロゲートモデル（NN, Gaussian Process等）を構築し、設計空間を高速探索する。SoftInWayのAxSTREAMやNUMECA FINE/Designにはこうした機能が統合されつつある。
比速度を超える設計変数空間——非次元ターボ機械設計の最適化フロンティア
比速度Nsは単一設計点でのターボ機械特性を表すが、実際の運転範囲全体をカバーする最適設計には多点最適化（Multi-Point Optimization）が必要だ。特に可変速ポンプ・風力タービンのように広い流量範囲で高効率を維持する「フラット高効率曲線」設計は、単一Nsによる最適化では実現できない。現代のCFD最適化では「設計変数空間（翼弦長、積み重ね線、カンバー分布）の多次元最適化」をアジョイント法（感度計算の高効率化）と応答曲面（Surrogate Model）を組み合わせて実施する。Siemens Energyが公開した産業用コンプレッサーの最適化事例では、設計点効率を+1.5ポイント、サージマージンを+8%改善しながら部分流量域の効率低下を±0.5%以内に抑えた結果が発表されており、比速度以上の設計自由度の重要性を示している。
トラブルシューティング
比速度範囲外の設計
比速度の推奨範囲外で設計するとどうなりますか？
効率が大幅に低下する。例えば比速度100以下の極低比速度遠心ポンプでは、ディスク摩擦損失が水力仕事と同等レベルになり、効率30～50%に留まることがある。
そういう場合はどうすればいいですか？
設計変更の選択肢を検討する。
| 問題 | 対策 |
|---|---|
| 比速度が低すぎる | 多段化、回転数増加、直径縮小 |
| 比速度が高すぎる | 多吸込、回転数低下、直径拡大 |
| 型式不適合 | 遠心→斜流→軸流の型式変更 |
CFDで確認すべきオフデザイン挙動
CFDでオフデザイン性能を確認する際のポイントは？
遮断運転もCFDで計算するんですか？
安全評価として重要だ。バルブが閉まった状態でポンプが回り続けると、液温が急上昇する。CFDで内部再循環のパワーを求め、温度上昇率を推定する。
単位系の注意
比速度の単位系が文献によって違って混乱します。
| 単位系 | $N$ | $Q$ | $H$ | 日本の慣例値例 |
|---|---|---|---|---|
| 日本慣例 | rpm | $m^3/min$ | m | $N_s = 200$ |
| SI（ヨーロッパ） | rpm | $m^3/s$ | m | $n_q = 26$ |
| US慣例 | rpm | US gpm | ft | $N_s = 1032$ |
| 無次元 | rad/s | $m^3/s$ | J/kg | $\omega_s = 0.54$ |
同じポンプでも単位系で数値が全く異なるから、論文やカタログ値を比較する際は必ず単位系を確認すること。
比速度Nsの計算値が設計仕様と合わない——単位系と基準流量の定義ミス
比速度Nsの計算で「CFDの結果からNsを算出したが設計仕様書のNsと合わない」問題の多くは単位系の不一致だ。Nsには「工学単位式（Ns=N*sqrt(Q)/H^(3/4)）」と「SI無次元形式」があり、係数が全く異なる（水の場合の換算係数は約4.46）。日本の旧JIS規格はNs単位系をベースにしており、欧米のSI無次元式とは数値が10倍以上異なる。設計仕様書のNs値がどの定義を使っているかの確認が第一歩だ。また「流量Q」の定義——全流量か1翼列通過流量か——によっても差が出る（多段の場合）。CFD結果からNsを計算する際は必ず入口全圧・出口全圧・質量流量・回転数を明示し、計算式の分子・分母に入れる変数の定義を設計チームで統一することが実務でのトラブル防止の鉄則だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——比速度と設計指針の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告