ポンプキャビテーション
理論と物理
概要
ポンプのキャビテーションって、泡がボコボコ出て壊れる現象ですよね？
大まかにはそうだ。液体の静圧が飽和蒸気圧を下回ると蒸気泡が生成され、下流の高圧域で崩壊する。この崩壊時に局所的に数GPaの衝撃圧が発生し、インペラ表面を浸食するんだ。
数GPa！？ それは壊れますね…
性能面でも深刻だ。揚程低下、振動増大、騒音増加が起きる。だからキャビテーション回避はポンプ設計の最優先課題の一つだよ。
NPSHの定義
NPSHってよく聞きますが、正確な定義は何ですか？
NPSH（Net Positive Suction Head）は吸込み側で液体がどれだけ蒸気圧から余裕があるかを表す。
$p_{atm}$:大気圧、$p_v$:飽和蒸気圧、$z_s$:液面からポンプ中心までの高さ、$h_f$:吸込配管の損失水頭。これがシステム側のNPSH（Available）だ。
ポンプ側にもNPSHがあるんですよね？
NPSH_r（Required）はポンプ自体が必要とする最低限のNPSHで、揚程が3%低下する点として定義される。安全にはこうなる。
キャビテーション係数（Thoma数）でも表せる。
Rayleigh-Plesset方程式
気泡の成長と崩壊を記述する式ってありますか？
Rayleigh-Plesset方程式が基本だ。
$R$:気泡半径、$p_B$:気泡内圧、$p_\infty$:周囲圧力、$S$:表面張力。CFDのキャビテーションモデルはこれを簡略化した質量輸送方程式に落とし込んでいるよ。
キャビテーションが潜水艦を悩ませた歴史
ポンプキャビテーションの気泡力学（Rayleigh-Plesset方程式）が実用的な問題として重視されるようになったのは、第一次世界大戦中の潜水艦プロペラからです。スクリュー翼端付近で気泡が激しく発生・崩壊し、推進効率の低下と金属浸食という二重の問題が起きた。戦後の研究でNPSH（正味吸込ヘッド）という概念が体系化され、現代のポンプ設計に引き継がれています。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
均質混合モデル
CFDでキャビテーションをどうやって計算するんですか？
最も広く使われるのは均質混合（Homogeneous Mixture）モデルだ。液相と蒸気相を単一流体として扱い、蒸気体積分率 $\alpha_v$ の輸送方程式を解く。
$\dot{m}^+$ が蒸発（気泡生成）、$\dot{m}^-$ が凝縮（気泡崩壊）のソースタームだ。
ソースタームのモデルにはどんなものがありますか？
代表的な3つを比較しよう。
| モデル | 特徴 | 使用ソルバー |
|---|---|---|
| Zwart-Gerber-Belamri | 核生成点密度ベース、パラメータ調整しやすい | CFX（デフォルト）、STAR-CCM+ |
| Schnerr-Sauer | Rayleigh-Plessetベース、気泡数密度指定 | OpenFOAM（interPhaseChangeFoam）、Fluent |
| Singhal (Full Cavitation) | 非凝縮ガスも考慮、実用的 | Fluent |
CFXではZwartモデルが標準なんですね。
そう。蒸発係数 $F_{vap}=50$、凝縮係数 $F_{cond}=0.01$ がデフォルト値。核生成点体積分率 $\alpha_{nuc}=5 \times 10^{-4}$、気泡初期半径 $R_B=10^{-6}$ m が典型値だ。
数値設定のポイント
キャビテーション計算で収束させるコツはありますか？
難しい計算だから、いくつかポイントがある。
- タイムステップ: 非定常が必須。翼1枚通過時間の1/20～1/50が目安
- 収束基準: RMS残差 $10^{-5}$ 以上を目標（キャビテーション振動で$10^{-4}$にしか下がらないこともある）
- 初期条件: まず非キャビテーション状態で定常収束させ、そこからキャビテーションモデルをON
- 圧縮性: 蒸気-液体の密度比が大きいため、数値的に圧縮性を考慮する設定が必要
Zwart・Merkle・Singhal——3モデルの争い
CFDキャビテーション解析の現場でよく聞く「どのモデルを使えばいい？」という問い。Zwartモデルは核生成サイト密度を明示的に扱い、Merkleモデルは圧力差に直接反応する経験的な形式、Singhal（Full Cavitation）モデルは溶存ガスまで考慮する最も複雑な構成です。ベンチマーク研究ではケースによって順位が入れ替わることもあり、「これが一番」とは言い切れないのが実情。初期値・メッシュ・経験定数の設定が結果を大きく左右するため、自社実験との検証が欠かせません。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
NPSH特性曲線の取得手順
CFDでNPSH曲線を描くにはどうすればいいですか？
入口全圧を段階的に下げていく方法が一般的だ。
1. 基準計算: 十分高いNPSHa（キャビテーションなし）で定常解を取得
2. 入口圧力低下: 入口全圧を0.1～0.2 atm刻みで低下
3. 各点で非定常計算: キャビテーションモデルONで数回転分の非定常計算
4. 時間平均揚程を記録: 揚程が基準値から3%低下する点がNPSH_r
0.1 atm刻みは結構粗くないですか？
NPSH_r付近で揚程が急降下するから、まず粗く全体を把握して、3%低下付近を0.02～0.05 atm刻みで追い込むのが効率的だ。
可視化と評価
キャビテーションの結果はどうやって見ればいいですか？
蒸気体積分率 $\alpha_v = 0.5$ の等値面がキャビティ形状を表す。CFD-Postで以下を確認しよう。
| 可視化項目 | 確認ポイント |
|---|---|
| $\alpha_v$ 等値面 | キャビティの位置・大きさ・翼面への貼り付き |
| 翼面圧力分布 | 吸い込み面のどこで蒸気圧を下回るか |
| 翼面浸食リスク指標 | 気泡崩壊圧力の累積（CFXのErosion Model） |
| 径方向のキャビティ分布 | ハブ側 vs. シュラウド側の違い |
浸食リスクもCFDで評価できるんですか？
CFXにはCavitation Erosion Modelが実装されていて、気泡崩壊時のエネルギー密度から浸食リスクマップを生成できる。ただし定量的な寿命予測には実験との校正が必要だ。
設計改善の方向性
キャビテーション性能を改善するにはどんなアプローチがありますか？
主な設計パラメータの影響をまとめよう。
| パラメータ | NPSH_rへの影響 | トレードオフ |
|---|---|---|
| インデューサ追加 | 大幅低減 | 構造が複雑化、コスト増 |
| 翼入口角の最適化 | 入口衝撃角低減で改善 | off-design性能への影響 |
| 翼枚数増加 | 翼負荷低減で改善 | 摩擦損失増大 |
| 吸込口径拡大 | 流速低下で改善 | ポンプ寸法増大 |
NPSH特性試験はなぜ3%落ちで定義するのか
ポンプキャビテーションの実験で「NPSH₃」という言葉が出てきます。これは「揚程Hが3%低下した時点のNPSH値」という定義ですが、なぜ3%なのでしょうか？ 実は明確な物理的根拠があるわけではなく、1940～50年代に国際標準として合意された経験的な閾値です。CFDでNPSH特性を予測する際は「揚程が3%落ちる気相体積分率の条件」を再現することになりますが、これが数値的にも安定して検出できるかどうかが実務での重要なチェックポイントになります。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
主要ソルバーのキャビテーション機能比較
キャビテーション解析に適したソルバーはどれですか？
各ソルバーの対応状況を比較しよう。
| 機能 | Ansys CFX | Ansys Fluent | STAR-CCM+ | OpenFOAM |
|---|---|---|---|---|
| 均質混合キャビテーション | Zwart（標準） | Schnerr-Sauer, Singhal | Rayleigh-Plesset | interPhaseChangeFoam |
| 浸食モデル | Erosion Model | UDF拡張 | Cavitation Damage | なし（自作可） |
| ターボ専用メッシャー | TurboGrid | Turbo Meshing | 内蔵テンプレート | なし |
| 非定常回転計算 | Transient Rotor-Stator | Sliding Mesh | Rigid Body Motion | cyclicAMI |
CFXが一番機能が揃ってそうですね。
ポンプのキャビテーション解析に関してはCFXが最も実績がある。結合型ソルバーの安定性が密度比の大きいキャビテーション計算に向いているし、TurboGridとのワークフローも成熟している。
Fluent ではどうですか？
Fluent 2024R1以降ではTurbomachinery workflowが強化されて、キャビテーション付きポンプ解析も実用レベルになっている。特にPolyhedral meshとの組み合わせで自動メッシュ生成が楽なのは利点だ。
OpenFOAMでのキャビテーション
OpenFOAMでもキャビテーション計算はできますか？
interPhaseChangeFoamソルバーがSchnerr-SauerやKunzモデルを実装している。ただし回転機械との組み合わせはcyclicAMIの設定が手間だし、収束安定性も商用に劣る場面が多い。研究用途や予算制約時の選択肢だ。
計算コストの見積もり
キャビテーション計算はどのくらい時間がかかりますか？
非キャビテーション計算の5～10倍と考えておくといい。非定常が必須で、かつ蒸気-液体の界面変動で小さいタイムステップが要求されるからだ。
| モデル規模 | セル数 | コア数 | 回転数分の計算時間 |
|---|---|---|---|
| 単段ポンプ（1翼列） | 200万 | 32 | 4～8時間 |
| インデューサ付き | 500万 | 64 | 12～24時間 |
| 多段ポンプ | 1000万 | 128 | 2～4日 |
STAR-CCM+がポンプ業界に強い理由
キャビテーション解析ツールとしてSTAR-CCM+が水処理・化学プラント業界で多く選ばれる背景には、回転領域（MRF・Sliding Mesh）とキャビテーションモデルのワンストップ設定の使いやすさがあります。特に「VOF+キャビテーションモデル」の組み合わせが安定していて、NPSH特性曲線を自動的にバッチ計算するワークフローを組みやすい。ANSYS Fluentも同等の機能を持ちますが、セットアップのGUI体験という点でSTAR-CCM+を好む実務者が多いのが実情です。
選定で最も重要な3つの問い
- 「何を解くか」：ポンプキャビテーションに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
非定常キャビテーション現象
キャビテーションが非定常だとどんな問題が起きますか？
翼面上のキャビティが周期的に成長・崩壊を繰り返す「クラウドキャビテーション」は特に深刻だ。崩壊時の衝撃波が下流の翼面を叩き、浸食が加速する。
クラウドキャビテーションの周波数はどのくらいですか？
ストルーハル数 $St = fL/V \approx 0.2 \sim 0.4$ 程度が報告されている。$f$:振動周波数、$L$:キャビティ長さ、$V$:主流速度。翼弦長50mm、流速20m/sなら100～160Hz程度だ。
水中騒音の予測
キャビテーション騒音もCFDで予測できますか？
Ffowcs Williams-Hawkings（FW-H）方程式を使ったハイブリッド手法が使われる。CFDで近傍場の非定常圧力を計算し、音響アナロジーで遠方場の音圧を算出する。
Fluent と STAR-CCM+ にはFW-Hソルバーが内蔵されている。
キャビテーションがあるとないとで騒音はどう変わりますか？
キャビテーション発生時は広帯域騒音が10～20dB増加する。特に気泡崩壊に起因する高周波成分（1～100kHz）が顕著だ。ポンプの騒音規制が厳しい集合住宅用途などでは、NPSH余裕を十分に確保する設計が求められる。
流体-構造連成（FSI）
キャビテーションによるインペラの振動も解析できますか？
CFDの非定常圧力場をFEMの構造モデルに転写する一方向FSIが実用的だ。キャビティの崩壊による衝撃加振と翼の固有振動数の共振を確認する。Ansys System CouplingやSTAR-CCM+のCo-Simulationで実行可能だ。
完全な双方向FSIまで必要ですか？
ポンプインペラは剛性が高いから、ほとんどの場合一方向FSIで十分だ。ただし樹脂製インペラや薄肉翼では流体-構造の双方向連成が必要になることもある。
Coffee Break よもやま話
ロケットエンジンのポンプとキャビテーション
H-IIAロケットのターボポンプは液体水素を-253℃で取り扱い、毎秒数百Lの液体を超高圧に昇圧します。この極限環境でキャビテーションを制御するために、インペラ前段に「インデューサ」と呼ばれる低NPSHの予備翼を設けます。インデューサが一部キャビテーションを許容しながら圧力を上げ、メインインペラには安定した流れを供給する設計です。この設計思想はJAXAとIHIが長年の実験とCFDで磨き続けてきたもので、宇宙開発とポンプキャビテーション研究が深く交差した領域です。
キャビテーションによるインペラの振動も解析できますか？
CFDの非定常圧力場をFEMの構造モデルに転写する一方向FSIが実用的だ。キャビティの崩壊による衝撃加振と翼の固有振動数の共振を確認する。Ansys System CouplingやSTAR-CCM+のCo-Simulationで実行可能だ。
完全な双方向FSIまで必要ですか？
ポンプインペラは剛性が高いから、ほとんどの場合一方向FSIで十分だ。ただし樹脂製インペラや薄肉翼では流体-構造の双方向連成が必要になることもある。
ロケットエンジンのポンプとキャビテーション
H-IIAロケットのターボポンプは液体水素を-253℃で取り扱い、毎秒数百Lの液体を超高圧に昇圧します。この極限環境でキャビテーションを制御するために、インペラ前段に「インデューサ」と呼ばれる低NPSHの予備翼を設けます。インデューサが一部キャビテーションを許容しながら圧力を上げ、メインインペラには安定した流れを供給する設計です。この設計思想はJAXAとIHIが長年の実験とCFDで磨き続けてきたもので、宇宙開発とポンプキャビテーション研究が深く交差した領域です。
トラブルシューティング
よくある失敗
キャビテーション解析でありがちなミスって何ですか？
経験上、頻出する問題を挙げよう。
1. 飽和蒸気圧の設定ミス
水の場合、25℃で3170Pa、80℃で47400Paと温度で大きく変わる。温水ポンプなのに常温の蒸気圧で計算して「キャビテーションしない」と誤判断するケースが多い。
2. 定常計算でキャビテーションを評価
定常でも一応キャビティは出ますよね？
定常計算ではキャビティが「凍結」された形で出るが、実際のキャビテーションは本質的に非定常だ。特にNPSH_r付近では定常計算の揚程が非定常の時間平均と大きくずれることがある。NPSH_rの評価には必ず非定常計算を使うべきだ。
3. メッシュ不足
キャビティの界面はシャープな密度勾配を持つから、メッシュが粗いとキャビティ形状がぼやけて体積を過大評価する。翼面吸い込み側のメッシュを特に細かくする必要がある。
実験との比較検証
CFDの結果はどうやって検証すればいいですか？
段階的に検証するのが良い。
| 検証段階 | 比較対象 | 許容誤差 |
|---|---|---|
| 非キャビテーション | QH特性曲線 | 揚程±3%以内 |
| キャビテーション開始 | キャビティの目視観察 | 発生位置の一致 |
| NPSH_r | 3%揚程低下点 | ±0.5m以内 |
| 浸食位置 | 実機の浸食パターン | 定性的一致 |
NPSH_rで±0.5mって結構シビアですね。
そう。CFDでNPSH_rを±0.3m以内で予測できれば非常に良い精度だ。それ以上の精度を求めるなら、キャビテーションモデルの係数を実験データで校正する必要がある。
温度効果（サーモダイナミック効果）
高温の液体だとキャビテーションが変わるんですか？
液体窒素やLNGなど低温流体、あるいは高温水では、気泡生成時に潜熱で周囲液温が下がり蒸気圧が低下する「サーモダイナミック効果」が発生する。結果としてキャビテーションが抑制される。このため高温水ポンプのNPSH_rは冷水より低くなる。CFXではこの効果を考慮するオプションがあるよ。
「キャビテーションしていないのに翼が削れた」事件
水力機械の現場で稀に起きる現象として、「CFDでキャビテーションが出ないのに翼面にエロージョンが見つかる」ケースがあります。原因として疑われるのが「渦キャビテーション」——チップ漏れ渦の中心部で局所的に蒸気圧を下回る現象です。定常RANS計算ではこの渦コアを過大に拡散させてしまうため検出されにくい。LESや非定常計算を使って渦コアの低圧を解像しないと、エロージョン予測が抜け落ちる典型的な落とし穴です。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——ポンプキャビテーションの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告