キャビテーション
理論と物理
概要
先生、キャビテーションって水中で泡ができる現象ですよね？
キャビテーションは、液体の局所圧力が飽和蒸気圧以下に低下したとき、液中に蒸気の空洞（キャビティ）が形成される現象だ。ポンプ羽根車、船舶プロペラ、弁の絞り部など高速流れ場で発生する。
沸騰とは違うんですか？
沸騰は温度上昇が駆動力、キャビテーションは圧力低下が駆動力だ。キャビテーションでは気泡が高圧域に移動すると急速に崩壊し、数千気圧の衝撃波と局所的な高温スポットが生じる。これがエロージョン（壊食）の原因になる。
支配方程式
CFDではどんなモデルを使うんですか？
まずキャビテーション数が基本パラメータだ。
CFDでは均質混合モデルが主流で、蒸気体積分率 $\alpha_v$ の輸送方程式にソース項を追加する。
ソース項のモデルにはどんな種類がありますか？
代表的なキャビテーションモデルを比較しよう。
| モデル | 基本概念 | 特徴 |
|---|---|---|
| Schnerr-Sauer | Rayleigh-Plesset方程式ベース | 気泡数密度 $n_0$ がパラメータ |
| Zwart-Gerber-Belamri | 簡略化RP方程式 | Fluent標準、調整係数で制御 |
| Singhal (Full Cavitation) | 質量輸送 | 非凝縮ガス（溶存空気）を考慮 |
| Kunz | 人工圧縮性ベース | 定常計算向き |
基礎となるRayleigh-Plesset方程式は球形気泡の成長を記述する。
第4項は粘性減衰、第5項は表面張力だ。CFDモデルでは2次項と粘性項を無視した簡略形から蒸発率を導出する。
気泡数密度 $n_0$ はどう決めるんですか？
一般水では $n_0 = 10^{13}$ /m³がデフォルト値だ。Schnerr-Sauerモデルでは $R_B = \left(\frac{3\alpha_v}{4\pi n_0}\right)^{1/3}$ で気泡半径を求め、蒸発率を計算する。
プロペラを壊す泡——キャビテーションが歴史を変えた瞬間
1893年、英国海軍の駆逐艦「Daring」は設計速度に全く届かず、プロペラに謎の損傷が続発しました。調査に当たったOsborne Reynoldsとその後継者たちが発見したのが「キャビテーション」現象です。局所圧力が水の蒸気圧（20℃で2.3 kPa）を下回ると蒸気泡が生成し、崩壊時に数百MPaの衝撃圧が発生してプロペラ材料を侵食します。この発見がきっかけとなりキャビテーション数σ = (p-pv)/(0.5ρu²)という無次元数が定義され、今日の水力機械設計の根幹を成しています。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
数値解法の詳細
キャビテーション解析の数値的なポイントは何ですか？
蒸気領域では混合体の音速が極端に低下し、圧縮性効果が顕著になる。水/蒸気の混合体音速は純水（約1500 m/s）より遥かに低く、数m/sまで下がることがある。
このため圧力-密度の結合が強く、圧力ベースソルバーではCoupled algorithmが推奨される。密度ベースソルバーの方が安定な場合もある。
乱流モデルの選択
乱流モデルは何を使えばいいですか？
標準k-εは乱流粘性を過大評価し、キャビティの非定常性を抑制する。Reboud修正が有効だ。
$n \approx 10$ で混合領域の乱流粘性を低減し、キャビティのシェディングを再現できる。SST k-ωもキャビテーション解析で良好な結果を示す。
非定常クラウドキャビテーションの詳細を捉えるにはDES、DDES、またはLESが必要になる。
ソルバー設定
| パラメータ | 推奨値 | 理由 |
|---|---|---|
| 圧力-速度連成 | Coupled | 圧力-密度の強い結合 |
| 空間離散化 | 2次精度以上 | キャビティ形状の解像 |
| 界面Courant数 | < 0.5 | 気泡成長/崩壊の捕捉 |
| 参照圧力 | 絶対圧基準 | 蒸気圧との比較のため |
OpenFOAMでの実装
OpenFOAMではどのソルバーを使いますか？
interPhaseChangeFoam がキャビテーション対応VOFソルバーだ。constant/transportProperties でモデルを指定する。SchnerrSauer、Kunz、Merkleから選択可能だ。
Fluentでの設定
Fluent側のポイントも教えてください。
Multiphase ModelでVOFまたはMixtureを選び、Cavitation Modelを有効化する。Zwart-Gerber-Belamriがデフォルトで、蒸発係数 $F_{evap} = 50$、凝縮係数 $F_{cond} = 0.01$ が標準値だ。非対称な係数は崩壊が蒸発より急速であることを反映している。
Schnerr-Sauer vs Zwart——キャビテーションモデル選択の実際
CFDでのキャビテーション解析で必ず議論になるのが質量輸送モデルの選択です。Schnerr-Sauerモデルは単一気泡の体積変化をRayleigh方程式から厳密に導出し、核密度の指定が不要という利点があります。一方、ZwartモデルはFluent標準で広く実績があり、蒸発・凝縮の非対称係数（Ce=0.02, Cc=0.01）によってヒステリシス挙動を再現できます。ポンプインデューサーの検証では、同一メッシュでも両モデルの予測するキャビテーション初生σが20%以上ずれることがあり、実験値との照合なしにモデルを選ぶのは危険です。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
実践ガイド
キャビテーション解析の実務的な手順を教えてください。
ポンプインペラの解析を例に説明しよう。
1. 単相定常解析: キャビテーションなしで流れ場を完全に収束
2. 蒸気圧設定: 動作温度での正確な飽和蒸気圧（25℃水で3170 Pa）
3. キャビテーション有効化: 単相解からリスタートしてモデルを有効化
4. NPSHの段階的低下: 入口圧力を下げてキャビテーション性能曲線を作成
5. 揚程3%低下点: NPSHrequiredを特定
メッシュ設計
キャビテーション解析に適したメッシュは？
| 領域 | 推奨サイズ | 理由 |
|---|---|---|
| 翼前縁 | 翼弦長の1/200 | キャビティ発生起点の解像 |
| 翼面吸い込み側 | 翼弦長の1/100 | キャビティ長さの再現 |
| 壁面プリズム層 | $y^+ \approx 1$ | 境界層剥離の正確な予測 |
| キャビティ後端 | 翼弦長の1/100 | re-entrant jetの捕捉 |
NPSHカーブって何ですか？
ポンプのキャビテーション性能を評価する標準的方法だ。
入口全圧を段階的に下げて各条件で定常計算を行い、揚程の変化をプロットする。揚程が3%低下した点がNPSH3%（required NPSH）になる。10〜15ケースの計算が必要で、前のケースからリスタートすると効率的だ。
エロージョン評価
壊食の予測もできますか？
定量的な壊食量予測は困難だが、壊食リスクの相対評価は可能だ。キャビティ崩壊時の壁面衝撃圧を指標にする。Fluentのcavitation erosion indicator、STAR-CCM+のerosion risk indexが利用可能だ。設計変更案の比較に有用だよ。
水力タービンの侵食予測——年間メンテコストを左右するCFD
ブラジルのItaipu水力発電所では、キャビテーション侵食による年間修繕コストが数十億円規模に達します。CFDによるキャビテーション予測は設計マージンを定量化するために使われますが、実務ではLES品質のシミュレーションが必要で、定常RANSでは気泡崩壊の非定常ダイナミクスが捉えられません。あるプロジェクトでは定常計算で「問題なし」と判断したランナーが、運転6か月で壊食が始まり、非定常計算をやり直したところキャビテーション雲の周期的崩壊が確認されたという事例があります。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
商用ツール比較
キャビテーション解析に対応しているツールを比較してください。
| ツール | モデル | 多相流手法 | 特徴 |
|---|---|---|---|
| Ansys Fluent | Zwart, Schnerr-Sauer, Singhal | VOF / Mixture | エロージョン指標あり |
| STAR-CCM+ | Schnerr-Sauer, Merkle | VOF | 音響解析連携（FW-H） |
| Ansys CFX | Rayleigh-Plesset | Homogeneous | 回転機械に強い |
| OpenFOAM | Schnerr-Sauer, Kunz, Merkle | VOF | 完全カスタマイズ可能 |
| NUMECA FINE/Turbo | 独自モデル | Homogeneous | ターボ機械専用 |
用途別推奨
用途によって使い分けるんですか？
STAR-CCM+のFW-H連携は何に使うんですか？
キャビテーションノイズの予測だ。キャビティ崩壊時の圧力変動を音源としてFfowcs Williams-Hawkings方程式で遠方場騒音を計算できる。潜水艦のステルス性能やソナー設計で重要だよ。
キャビテーション解析市場——水力機械メーカーが使うツールの実態
キャビテーション解析でのソフトウェア選択は業界によって明確に分かれます。水力タービン・ポンプメーカー（Voith、KSB等）ではANSYS CFX/Fluentが圧倒的シェアを持ち、CFXの結合ソルバーがキャビテーション収束安定性の面で評価されています。船舶・潜水艦分野ではFine/Marine（NUMECA）が流体-騒音連成解析ツールとして選ばれることが多く、キャビテーション騒音の周波数特性予測に強みがあります。自動車燃料噴射弁のキャビテーション解析ではマルチフィジクス対応のStar-CCM+が設計ループの速さから選ばれています。
選定で最も重要な3つの問い
- 「何を解くか」：キャビテーションに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
先端技術と研究動向
キャビテーション研究の最先端にはどんなトピックがありますか？
いくつかの方向性がある。
非定常クラウドキャビテーション
シート状キャビティの後端が周期的に崩壊してクラウド（蒸気塊）を放出する。Re-entrant jetや衝撃波の伝播による崩壊メカニズムが競合し、$St = fL_c / U_{\infty} \approx 0.2$〜0.3 の周期性を持つ。
LES/DESでこの非定常挙動を捕捉する研究が活発で、Schnerr（TUM）やBensow（Chalmers）のグループが先端的な成果を出している。
熱力学的効果
低温流体のキャビテーションは何が違うんですか？
液体水素や液体窒素では蒸発潜熱の吸収により局所温度が低下し、蒸気圧も下がる。この熱力学的効果がキャビテーションを抑制する。B-factorで評価する。
ロケットエンジンのターボポンプ設計（NASA、JAXA）で重要なテーマだ。Fluent、CFXに熱力学効果を追加したモデルがある。
気泡崩壊のDNS
壁面近傍での気泡崩壊では非球形変形が起き、壁面に向かうマイクロジェット（数百 m/s）が発生する。これを直接解くDNSが進んでおり、Tinguely（ETH）やJohnsen（Michigan）のグループが代表的だ。
機械学習の応用
キャビテーション発生の予測やエロージョンリスク評価にCNN/GANを適用する研究が増えている。CFDの非定常計算は計算コストが高いため、サロゲートモデルによるリアルタイム予測が期待されている。
熱力学的キャビテーション——冷水と熱水では別の物理
水温が80℃を超えると、キャビテーション気泡の成長に熱拡散が支配的になり、古典的な等温Rayleigh理論が破綻します。この「熱力学的キャビテーション」では蒸発に必要な熱が周囲液体から奪われるため、蒸発が自己抑制的になりキャビテーション強度が低下します。液体窒素や液体水素を扱う宇宙ロケットのターボポンプ設計では必ずこの効果を考慮する必要があり、NASAの研究では水温が100℃から120℃に上がるだけでキャビテーション初生σが0.3以上低下することが示されています。
トラブルシューティング
トラブルシューティング
キャビテーション解析でよくあるトラブルを教えてください。
順番に見ていこう。
1. キャビテーションが発生しない
症状: 実験では発生する条件なのにCFDで蒸気が生成されない。
対策:
- 飽和蒸気圧を正確に設定（25℃水で $p_v \approx 3170$ Pa）
- ゲージ圧と絶対圧の混同がないか確認
- 低圧部のメッシュを細分化
- 乱流粘性の過大評価をSST k-ωやReboud修正で抑制
2. キャビティが異常に大きい
対策:
- 凝縮係数を増やす（Zwartモデル: $F_{cond}$）
- 溶存空気を考慮するSinghalモデルに切り替え
- 出口境界が十分下流にあるか確認
3. 計算が発散する
発散の原因は何ですか？
蒸気領域での急激な密度変化が圧力方程式を不安定にする。対策としては単相解を完全収束→キャビテーション有効化の2段階アプローチ、Coupled solver使用、十分に小さいタイムステップが有効だ。
4. 非定常周期が実験と不一致
対策:
- RANSからDES/LESに切り替え
- キャビティ後端のメッシュ解像度向上
- タイムステップをキャビテーション周期の1/100以下に
- 出口境界の圧力反射を回避
5. ツール固有の注意点
| ツール | 注意点 |
|---|---|
| Fluent | Operating Pressureに注意。絶対圧基準が推奨 |
| CFX | Homogeneous multiphaseの凝縮ソース符号をバージョン確認 |
| STAR-CCM+ | VOFシャープネスとキャビテーションモデルの整合性確認 |
| OpenFOAM | SchnerrSauerの気泡数密度パラメータの感度が高い |
シミュレーションは発生しないのに実機では壊れる——キャビテーション診断の落とし穴
キャビテーションCFDで最も多い相談は「計算ではキャビテーションが出ないのに実機ではボコボコ音がする」です。原因の多くはメッシュの粗さにあります。キャビテーション初生は局所的な圧力最小点で起きるため、翼前縁や弁座コーナーのメッシュ解像度が不足するとその圧力谷が平均化されて消えてしまいます。目安はy+ < 1かつ翼厚方向に30層以上の境界層メッシュです。入口乱流強度の設定が5%から0.1%に変わるだけでキャビテーション初生位置が大きくずれる事例も報告されており、実機の入口条件を正確に再現することが不可欠です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——キャビテーションの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告