プロペラCFD解析
理論と物理
プロペラ性能の基礎
先生、船舶プロペラの性能をCFDで予測するって、どういうパラメータを見るんですか？
プロペラ性能は推力係数 $K_T$、トルク係数 $K_Q$、効率 $\eta_0$ で評価する。
ここで $T$ は推力、$Q$ はトルク、$n$ は回転数、$D$ はプロペラ直径、$J = V_A/(nD)$ は前進係数だ。オープンウォーター特性（$K_T$-$J$曲線と$K_Q$-$J$曲線）がプロペラの基本性能を表す。
CFDではどんな結果と比較するんですか？
ITTC（国際水槽試験会議）やSMPF（船舶性能推進委員会）の標準プロペラ（例：DTMB 4119、KCS）の水槽試験データと比較する。$K_T$ と $K_Q$ でCFDと実験の差が2%以内なら良好、5%以内なら実用レベルだ。
数値モデリングのアプローチ
プロペラCFDにはどんなモデリング手法がありますか？
主に3つある。(1) MRF（Moving Reference Frame）: プロペラ周りを回転座標系で定常計算。最もコストが低い。(2) Sliding Mesh（Rigid Body Motion）: プロペラを実際に回転させる非定常計算。船体-プロペラ干渉が正確。(3) Overset Mesh（Chimera法）: 背景メッシュと重ね合わせで運動を扱う。STAR-CCM+で広く使われる。
どれを使うべきですか？
オープンウォーター特性の算出ならMRFで十分。船尾伴流中のプロペラ（自航状態）ではSliding MeshかOversetが必要。キャビテーション解析では非定常のSliding Meshが標準だ。
船舶プロペラ理論の歴史——ランキン-フルードの運動量理論（1865年）
プロペラ推力の最初の理論的記述はWilliam Rankine(1865)とW.J.M. Froude(1878)による「運動量理論（Actuator Disk Theory）」だ。プロペラを仮想の無限薄ディスクとして流れに運動量を与える理想化で、推力T = ρA(V+va)×2vaという単純な式（vaは誘導速度）を導いた。この理論は推力を発生させる物理の本質を捉えながら計算が単純で、今も設計初期段階の出力・効率見積もりに使われている。プロペラ効率の理論上限（Betz限界 η=1/(1+va/V)）もこの理論から導かれ、風力タービンのBetz係数（η_max≈59.3%）と同じ数学的構造を持つ。130年後のCFDはBladeを詳細にモデル化した完全な翼型解析を可能にしたが、Rankineの見通しの良さは今でも教育現場で語り継がれる。
各項の物理的意味
- 時間項 $\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$。時間刻みの安定性に直結 |
数値解法と実装
メッシュ設計
プロペラのメッシュ生成で気をつけることは何ですか？
翼面近傍の境界層メッシュが最重要だ。$y^+ \approx 1$ のprism layerを配置し、SST k-ωモデルで壁面を直接解像する。翼前縁・後縁には局所的な細分化が必要で、前縁曲率半径に対して少なくとも10セルを配置する。翼端渦を捕捉するにはtip近傍のメッシュも細かくする必要がある。
セル数の目安は？
1翼あたり50万〜200万セルが標準。周期境界条件を使って1翼分だけ計算すればコストを1/$Z$（$Z$は翼枚数）に削減できる。ただし船尾伴流中の自航計算では非一様流入なので全翼モデルが必須だ。STAR-CCM+のTrimmed MeshやFluentのPolyhedral Meshが品質と自動化のバランスに優れている。
キャビテーション解析
キャビテーションのCFDモデルはどうなっていますか？
Schnerr-Sauer modelやZwart-Gerber-Belamri modelが主流だ。VOF法で気液二相流を解き、局所圧力が蒸気圧を下回る領域で気泡の生成（蒸発）を、圧力が回復する領域で消滅（凝縮）をモデル化する。
FluentではMultiphase > VOF > Schnerr and Sauerで有効化。STAR-CCM+ではCavitation ModelをPhysicsで設定する。OpenFOAMではinterPhaseChangeFoam（Schnerr-Sauerモデル内蔵）を使う。
キャビテーション解析で注意すべき点は？
(1) 時間刻みが非常に小さい必要がある（$\Delta t \sim 10^{-5}$〜$10^{-4}$s）。(2) メッシュがキャビティ厚さ（数mm）を解像できるほど細かいこと。(3) 乱流モデルはSST k-ωにReboud correction（密度補正）を加えるとシート・クラウドキャビテーションの挙動が改善する。
キャビテーションモデルの「Schnerr-Sauer」はなぜ人気なのか
プロペラCFDでキャビテーションを解析する際、OpenFOAMユーザーに最もよく使われるのがSchnerr-Sauerモデルだ。Zwart-Gerber-Belamriモデルなど複数の選択肢があるが、Schnerr-Sauerは蒸発・凝縮の係数チューニングが不要で、泡核密度さえ設定すれば動く手軽さが支持される理由。ただし実際のキャビテーション崩壊（インプロージョン）の精細な再現にはLESや高密度メッシュが必要で、10億セル超の計算になることも珍しくない。現場では「全翼面均一6層プリズム層」が崩れないメッシュ品質管理が最初の難関になる。
風上差分（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次だが、流れの方向を正しく捕捉するため安定性が高い。
実践ガイド
自航解析（Self-Propulsion）
自航解析ってどういうシミュレーションですか？
船体とプロペラを同一の計算領域に配置し、プロペラが船を推進する状態をシミュレーションする。推力と抵抗が釣り合う自航点（Self-Propulsion Point）を求めることで、実船の必要馬力を予測する。
計算手順はどうなりますか？
(1) まず船体のみの抵抗計算（裸船体抵抗）を行う。(2) 次にプロペラのオープンウォーター特性を別途計算する。(3) 船体+プロペラのSliding Mesh計算で、回転数を調整しながら推力=抵抗となる点を探す。STAR-CCM+のPropeller Performance機能は自動的にself-propulsion pointを探索してくれる。
船尾伴流の影響はどうやって評価するんですか？
プロペラ面での伴流分布（wake fraction）$w = 1 - V_A / V_S$ を計算する。公称伴流（nominal wake、プロペラなし）と有効伴流（effective wake、プロペラあり）の両方を評価する。伴流の非一様性がプロペラの変動荷重、振動、キャビテーションの原因になるので、周方向分布を詳細に分析することが重要だ。
検証のためのベンチマーク
プロペラCFDの検証に使えるデータはありますか？
代表的なベンチマークデータを挙げよう。
| プロペラ | 翼数 | 特徴 | データ元 |
|---|---|---|---|
| DTMB 4119 | 3翼 | 非キャビテーション検証用 | ITTC |
| PPTC VP1304 | 5翼 | キャビテーション/圧力変動 | SVA Potsdam |
| KCS | 5翼 (KP505) | 自航解析用。KCS船体と組合せ | NMRI/SIMMAN |
これらのデータはSIMMANワークショップやITTC Benchmark Casesとして公開されている。
LNG船の省エネプロペラ——バラスト条件と満載条件の両立設計
LNG（液化天然ガス）輸送船はバラスト状態（空荷）と満載状態で喫水が大きく変わり、推進効率の要求が全く異なる。従来は満載条件に最適化したプロペラ設計が主流だったが、国際海事機関（IMO）のCII（Carbon Intensity Indicator）規制強化で両条件での燃費効率が重要になった。CFD（RANS + 自由表面VOF法 + プロペラ回転モデルMRF）による最適化では、ピッチ分布とキャンバーを両条件で妥協点を探るMulti-Point最適化が行われる。ある造船メーカーの設計事例では、CFDによる形状最適化で年間燃料消費量を3.8%削減し、CO₂排出量削減とCII等級改善を同時に達成した結果が造船学会で発表されている。
解析フローのたとえ
CFDの解析フローは「水族館の水槽を設計する」感覚で考えてみてください。まず水槽の形を決め（計算領域）、水の入り口と出口を設計し（境界条件）、ポンプの強さを設定する（流量条件）。魚がどう泳ぐか見たければ粒子追跡。水温が気になれば熱解析を追加。…どうですか？ 意外と直感的ではありませんか？
初心者が陥りやすい落とし穴
「y+って何ですか？」——この質問が出たら要注意。壁面近くのメッシュ解像度を表すy+は、CFDの結果精度を左右する最重要パラメータの1つ。壁関数を使うなら30〜300、壁を完全に解像するなら1以下。これを確認せずに「摩擦抵抗が合わない！」と悩む人がとても多い。体温計の先端をちゃんと脇に挟まないで「熱がないのに37.5度って出た！」と慌てているようなものです。
境界条件の考え方
入口の境界条件は「蛇口をどのくらい開けるか」と同じ。ちょろちょろ出すか（低速）、全開にするか（高速）。でもCFDではもう一つ——「どのくらい暴れた水を出すか」（乱流強度）も指定する必要があります。蛇口の開け方を間違えると、下流のシンク全体の流れが変わりますよね？ CFDでも入口条件のミスは下流全体に波及します。
ソフトウェア比較
STAR-CCM+でのプロペラ解析
STAR-CCM+はプロペラ解析に強いと聞きました。
STAR-CCM+はMarine業界で非常に広く使われている。Rigid Body Motion（Sliding Mesh相当）やOverset Meshが安定しており、DFBI（Dynamic Fluid Body Interaction）との統合が優れている。プロペラ解析の推奨設定は、Segregated Flow + SST k-ω + VOF（キャビテーション時）+ Rigid Body Motion。Polyhedral Mesh + Prism Layerの組合せが品質面で推奨される。
Ansys Fluentでのプロペラ解析
Fluentの場合は？
FluentではSliding Mesh（Mesh Motion）でプロペラ回転を扱う。MRF（Frame Motion）でオープンウォーター特性を求め、Sliding Meshで自航解析を行う二段階アプローチが効率的。圧力ベースCoupled solverにSST k-ωの組合せが標準。キャビテーションはMultiphase VOF + Schnerr-Sauer。Fluentのpolyhedral meshing（Fluent Meshing内）でプロペラ翼面の高品質メッシュが生成できる。
OpenFOAMでのプロペラ解析
OpenFOAMではどうですか？
MRFならsimpleFoamまたはpimpleFoamにMRFZone設定を追加。Sliding MeshはpimpleFoam + dynamicMeshDictのsolidBodyMotionFvMeshで回転を指定する。メッシュはsnappyHexMeshで生成できるが、プロペラのような複雑形状では商用メッシャー（Pointwise、ANSA等）からインポートするほうが品質が安定する。
キャビテーション解析はOpenFOAMでも可能ですか？
interPhaseChangeFoamがSchnerr-Sauer modelを内蔵している。ただし商用ソルバーと比べるとキャビテーション解析の安定性がやや劣ることがあり、時間刻みやメッシュ品質への感度が高い。OpenFOAM v2212以降では改善が進んでいるよ。
船舶プロペラCFDツール比較——MPUF-3A vs StarCCM+ vs OpenFOAMのKTとKQの精度
プロペラの推力係数KT・トルク係数KQの予測精度でCFDツールを比較すると、海洋分野の専用コード（MPUF-3A等の揚力面法）は計算が速いが翼端渦やキャビテーションの詳細な再現は不得意。StarCCM+はMRF（複数参照系）とSliding Meshを組み合わせたプロペラ解析が成熟しており、ITTCの比較テストでKT誤差±3〜5%、KQ誤差±5〜8%程度が報告されている。OpenFOAMのpimpleFoam + AMI（Arbitrary Mesh Interface）は研究機関でのオープンな検証が進んでおり、プロペラ開水面特性のbenchmark計算では商用ツールと遜色ない精度が公開論文で示されている。ただし設定の手間と自動化のしやすさでは商用ツールが依然優位だ。
選定で最も重要な3つの問い
- 「何を解くか」：プロペラCFD解析に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
水中放射騒音（URN）
プロペラの騒音予測もCFDでやるんですか？
IMO（国際海事機関）の水中騒音低減ガイドラインの影響で、プロペラの水中放射騒音（Underwater Radiated Noise）予測が重要になっている。騒音の主因はキャビテーションで、キャビティの崩壊時に高圧パルスが発生する。
CFDではFfowcs Williams-Hawkings（FW-H）方程式を使って遠方場の音響を予測する。FluentとSTAR-CCM+にはFW-Hソルバーが内蔵されており、非定常プロペラ計算の結果から騒音スペクトルを算出できる。ただしキャビテーション騒音の正確な予測にはメッシュ解像度と時間分解能の両方が極めて高い必要がある。
エネルギー節減デバイス（ESD）
EEDI（Energy Efficiency Design Index）対応でプロペラ周辺のデバイスが注目されていますね。
プロペラ前方のPre-Swirl Device（ステーター、ダクト）やプロペラ後方のPost-Swirl Device（PBCF：Propeller Boss Cap Fin、Hubvane）が代表的だ。これらはプロペラ周辺の流れを最適化して推進効率を2〜6%向上させる。
ESDのCFD評価はどうやりますか？
ESD装着時と未装着時の自航計算を行い、必要馬力の差を定量化する。ESDの形状最適化にはパラメトリックスタディやアジョイント法が使われる。STAR-CCM+のDesign ManagerやFluentのAdjoint Solverで、ESD翼角度やコード長の最適値を探索できる。
電動推進とリムドライブ
電動推進のプロペラ設計にもCFDは使われますか？
使われるよ。アジマススラスター（ABBのAzipod、Rolls-RoyceのMermaid等）やリムドライブ（rim-driven thruster）のCFDでは、モーター部の隙間流れや電磁気力との連成が追加的な課題になる。ポッド（nacelle）やストラットの形状最適化もCFDの領域だ。
プロペラキャビテーション侵食の予測——CFD-侵食モデル連成の最先端
船舶プロペラのキャビテーション侵食は、気泡崩壊時の局所衝撃圧（数GPa）が金属表面を繰り返し叩くことで生じる。CFD（Reynolds平均+キャビテーションモデル）で気泡生成・崩壊の空間分布を求め、崩壊エネルギーを侵食リスク指標（Erosion Risk Indicator: ERI）に変換する連成手法が最先端だ。特に「ジェット気泡崩壊（Microjet）」による局所圧力の方向性（壁面法線方向か否か）が侵食強度を決め、単純なキャビテーション体積分率だけでは侵食部位を誤予測する。DNV GLが公開したベンチマーク試験では、CFDキャビテーションモデルの侵食予測が実験の侵食跡と空間的に一致するまでに5世代の改良を要した歴史がある。
トラブルシューティング
$K_T$ が実験値と10%以上ずれる
オープンウォーター計算で $K_T$ が実験値と合いません。
まず以下を確認しよう。(1) 翼面メッシュが前縁曲率を十分に解像しているか（前縁に10セル以上）。(2) $y^+ < 1$ を達成しているか。特にSST k-ωでは $y^+$ が5〜30の中途半端な値だと精度が悪化する。(3) 翼端の後流が十分に解像されているか。メッシュ感度分析で収束を確認する。
MRFとSliding Meshで結果が違うんですが、どちらが正しいですか？
オープンウォーター条件（一様流入）ならMRFとSliding Meshの定常平均値はほぼ一致するはず。差が5%以上あるなら、MRFのインターフェース面でのメッシュ不整合が原因の可能性がある。FluentではMesh InterfaceのInterpolation方法をConservativeに設定し、orphan cellsがないか確認しよう。
Sliding Meshで圧力が振動する
Sliding Mesh計算で圧力が大きく振動して発散しそうです。
時間刻みが大きすぎる可能性がある。プロペラが1ステップで回転する角度を1度以下にするのが目安。$n = 10$ rpsなら $\Delta t < 1/(360 \times 10) \approx 2.8 \times 10^{-4}$s。また、sliding interface上のメッシュサイズが回転域と固定域で大きく異なると補間エラーが出る。interface面のセルサイズを揃えるのが重要だ。
キャビテーション解析が不安定
キャビテーション解析が数回転で発散します。
キャビテーションモデルの追加は数値的に不安定になりやすい。(1) まず非キャビテーション条件で十分に収束した流れ場を初期値にする。(2) 時間刻みをさらに小さくする（$\Delta t \sim 10^{-5}$s）。(3) VOFのunder-relaxationを0.3程度に下げる。(4) Courant数を0.5以下に抑える。STAR-CCM+のAutomatic CFL Ramping機能は初期の不安定期を自動的に乗り越えるのに有効だよ。
プロペラCFD推力が実験と15%ズレる——インペラモデルの設定と回転領域の境界処理
プロペラCFDで実験と15%以上の推力誤差が出る典型的原因は2つある。第一は「MRF（移動参照系）領域の設定ミス」——回転領域の境界をブレード先端に近づけすぎると、回転座標系から静止座標系への速度変換で誤差が増幅する。MRF境界は最低でもプロペラ半径Rの1.5倍以上離すことが推奨される。第二は「スリップ速度の扱い」——プロペラ後流の強い旋回（スワール）成分が出口境界まで伝播する場合、出口境界を後流域の3〜5D以上下流に設置しないと逆流が発生し推力計算が不安定になる。自由表面を含む場合、水中の換気（Air Ingestion）でプロペラ回転数が大きく変動する場合があり、VOF法との連成が必須になる。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——プロペラCFD解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告