プロペラキャビテーションFSI
理論と物理
キャビテーションFSIの概要
プロペラのキャビテーションでFSIが必要になるのはなぜですか？
キャビテーション（空洞現象）によってプロペラ翼面に非定常圧力変動が作用し、翼の振動・エロージョン・騒音を引き起こす。翼の弾性変形がキャビテーションパターンを変化させるため、流体-構造連成が必要になる。
支配方程式
キャビテーションの数理モデルはどうなっているんですか？
均質混合流モデルが広く使われる。液相と気相の混合密度を用いてNavier-Stokes方程式を解く。キャビテーション数は、
$p_v$ は蒸気圧だ。相転移はRayleigh-Plesset方程式に基づく質量輸送モデルで記述する。Schnerr-SauerモデルやZwart-Gerber-Belamriモデルが代表的だ。
$\alpha$ は気相体積分率、$\dot{m}^+, \dot{m}^-$ は蒸発・凝縮の質量転移率だ。
構造側はどうモデル化するんですか？
プロペラ翼はシェル要素またはソリッド要素でモデル化する。複合材プロペラ（CFRP等）の場合は積層構成を考慮した異方性材料モデルが必要だ。FSI界面でのデータ転送はwet surface上で行う。
キャビテーション数σ——たった一つの無次元数が運命を決める
プロペラ設計者が最初に計算するのがキャビテーション数σ＝(p∞－pv)/(½ρU²)だ。σが1を下回るあたりから気泡が発生し始め、0.3を切るとプロペラ全面が蒸気に包まれる「スーパーキャビテーション」に突入する。面白いのは、スーパーキャビテーションをあえて利用した魚雷（ロシアの「シクヴァル」）が水中340km/hを叩き出したこと。通常の魚雷の4倍超の速度だ。民間船では絶対に避けるべき現象が、使い方次第で兵器の極意になる——この逆転の発想がキャビテーション理論の醍醐味でもある。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常の例】夏に線路のレールが伸びて隙間が狭くなる——温度上昇→熱膨張→応力発生の典型例。電子基板がはんだ付け後に反るのも、異なる材料の熱膨張率差による。エンジンのシリンダーブロックは高温部と低温部の温度差で熱応力が発生し、最悪の場合亀裂に至る。
- 流体-構造連成（FSI）項：流体圧力・せん断力が構造を変形させ、構造変形が流体領域を変化させる双方向の相互作用。【日常の例】強風で吊り橋のケーブルが振動する（渦励振）——風の力が構造を揺らし、揺れた構造が風の流れを変え、さらに振動が増幅する。心臓の血流と血管壁の弾性変形、航空機の翼のフラッタ（空力弾性不安定性）も典型的なFSI問題。片方向のみの連成で済む場合もあるが、変形が大きい場合は双方向連成が必須。
- 電磁-熱連成項：ジュール発熱 $Q = J^2/\sigma$ が温度上昇を引き起こし、温度変化が電気抵抗を変化させるフィードバックループ。【日常の例】電気ストーブのニクロム線は電流が流れると発熱（ジュール熱）して赤くなる——温度が上がると抵抗が変わり、電流分布も変化する。IHクッキングヒーターの渦電流発熱、送電線の温度上昇による弛み増加もこの連成の例。
- データ転写項：異なる物理場間のメッシュ不一致を補間で解決。【日常の例】天気予報で「気温のデータ」と「風のデータ」を合わせて体感温度を計算するとき、それぞれの観測地点が異なれば補間が必要——CAEの連成解析でも、構造メッシュとCFDメッシュは一般に一致しないため、界面でのデータ転写（補間）精度が結果の信頼性に直結する。
仮定条件と適用限界
- 弱連成仮定（片方向連成）：一方の物理場が他方に影響するが逆は無視可能な場合に有効
- 強連成が必要なケース：FSIでの大変形、電磁-熱連成での温度依存性が強い場合
- 時間スケールの分離：各物理場の特性時間が大きく異なる場合、サブサイクリングで効率化可能
- 界面条件の整合性：連成界面でのエネルギー・運動量保存が数値的に満たされることを確認
- 適用外ケース：3つ以上の物理場が同時に強く連成する場合、モノリシック手法が必要になることがある
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 熱膨張係数 $\alpha$ | 1/K | 鋼: 約12×10⁻⁶、アルミ: 約23×10⁻⁶ |
| 連成界面力 | N/m²（圧力）またはN（集中力） | 流体側と構造側で力の釣り合いを確認 |
| データ転写誤差 | 無次元（%） | 補間精度はメッシュ密度比に依存。5%以下が目安 |
数値解法と実装
回転体のFSI処理
回転するプロペラのFSIはどう処理するんですか？
滑り面（sliding mesh / rotating reference frame）を使って回転流体領域を設定する。FSI連成ではプロペラ翼面を界面とし、回転座標系で構造解析を実行する。
Ansys FluentではMRF（Moving Reference Frame）で定常計算後、sliding meshで非定常計算に切り替えるのが標準的だ。キャビテーションの非定常性を捕捉するには必ずsliding meshが必要だ。
キャビテーションモデルの設定
キャビテーションモデルのパラメータはどう設定するんですか？
Schnerr-Sauerモデルの場合、bubble number density $n_b$ がキーパラメータだ。デフォルトは $n_b = 10^{13}$ /m³だが、実験との比較でチューニングが必要なことがある。
| パラメータ | Schnerr-Sauer | Zwart-Gerber-Belamri |
|---|---|---|
| 核密度/核半径 | $n_b = 10^{13}$ /m³ | $R_b = 10^{-6}$ m |
| 蒸発係数 | 1.0 | $F_{vap} = 50$ |
| 凝縮係数 | 1.0 | $F_{cond} = 0.01$ |
| 飽和蒸気圧 | 温度依存 | 温度依存 |
時間刻みの設定は？
プロペラ1回転あたり360〜720ステップ（0.5°〜1.0°/step）が推奨される。キャビテーションの崩壊過程は極めて速いため、局所的にはCFL < 1が望ましい。
BEMとCFD——プロペラ解析の「二刀流」使い分け
プロペラの流体解析には境界要素法（BEM）とCFDの2つが長年併用されてきた。BEMはポテンシャル流理論に基づくため計算が爆速で、設計初期段階の何百ケースもの条件スイープに向いている。一方、CFDのLES（大渦シミュレーション）はキャビテーションの泡の生成・崩壊まで捉えられるが、1ケースあたり数百CPU時間かかる。実際の船舶設計では「BEMで粗探索→CFDで精査→水槽実験で検証」という3段構えが標準だ。最近はBEMの結果をCFDの初期条件にするハイブリッド手法も増えており、総解析時間を従来比30〜40%削減できる事例も出てきている。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
実践ガイド
解析手順
プロペラキャビテーションFSI解析の実務的な手順は？
1. プロペラ翼形状のCADデータ取得
2. 流体メッシュ生成（回転領域+静止領域）
3. 構造FEモデル作成（翼のみ。ハブは剛体）
4. wet-non-cavitating条件でFSIの動作確認
5. キャビテーションモデルを有効化して本計算
6. 推力・トルク係数、キャビテーションパターンの評価
メッシュはどのくらい細かくする必要がありますか？
翼面の $y^+$ は1以下が望ましい。特に翼端部とリーディングエッジ周辺はキャビテーション発生の起点になるため細密化する。
| 領域 | 要素サイズ目安 |
|---|---|
| 翼面境界層 | $y^+ < 1$、15層以上 |
| 翼端渦域 | 翼弦長/100 |
| 後流域 | 翼弦長/50 |
| 遠方 | 翼弦長/5 |
検証にはどんなベンチマーク問題がありますか？
PPTCプロペラ（Potsdam Propeller Test Case）がITTC（国際水槽会議）の標準ベンチマークだ。キャビテーションパターン、推力係数 $K_T$、トルク係数 $K_Q$ の実測値が公開されている。
大型タンカーのプロペラ交換——現場での「音」診断
ベテランの船舶エンジニアは「プロペラのキャビテーションは音でわかる」と言う。正常なら低くて穏やかな回転音だが、キャビテーションが起きると船尾から「バリバリ」「ガリガリ」という異音が聞こえてくる。実際、30万トン級のVLCC（超大型タンカー）では、プロペラ径が9m以上にもなり、回転数はたった毎分90回転ほど。それでもキャビテーションが発生すると、3か月で羽根表面に1cm以上の侵食孔ができることがある。シミュレーションで侵食量を予測し、ドック入り時期を最適化できれば、年間数千万円のコスト削減になる。
解析フローのたとえ
風船を膨らませたことがありますか？ あの瞬間、実は高度な流体-構造連成が起きています。内部の空気圧（流体）がゴム壁（構造）を押し広げ→広がった壁が内部の圧力分布を変え→変わった圧力がさらに壁を変形させる…このキャッチボールを計算ステップごとに繰り返すのがFSI解析です。
初心者が陥りやすい落とし穴
「片方向連成で十分でしょ？」——この判断ミスが連成解析で最も危険です。構造の変形が微小なら確かに片方向で足りますが、心臓弁の開閉のように変形が流路を大きく変える場合、片方向では全く話になりません。目安は「変形量が代表長さの1%を超えるか」。超えるなら双方向連成は必須です。片方向で済ませてしまった場合、結果が「もっともらしいけど実は大間違い」になる——これが最も怖いパターンです。
境界条件の考え方
連成界面のデータ交換は「国境の出入国管理」と同じです。各国（物理場）には独自の法律（支配方程式）がありますが、国境（界面）で人や物（力・温度・変位）のやり取りを正確に管理しないと、両国の経済（エネルギーバランス）が崩壊します。メッシュが一致していない場合の補間は「通訳」のようなもの——誤訳（補間誤差）が小さいほど良い結果が得られます。
ソフトウェア比較
ツール比較
プロペラキャビテーションFSIに使えるソフトウェアは？
| ツール | キャビテーションモデル | FSI | 特徴 |
|---|---|---|---|
| Ansys Fluent | Schnerr-Sauer, Zwart | System Coupling | 回転体FSIの実績豊富 |
| STAR-CCM+ | Schnerr-Sauer | 内蔵FSI | ポリヘドラルメッシュ。海事産業に強い |
| OpenFOAM (interPhaseChangeFoam) | Merkle, Kunz等 | preCICE連携 | OSS。カスタマイズ自由 |
| FINE/Marine (Cadence) | Merkle | FSI対応 | 旧NUMECA。船舶CFD専用 |
| HydroStar + NASTRAN | BEM(ポテンシャル) | モード重畳 | 簡易FSI。初期設計向き |
STAR-CCM+は海事分野で人気があるんですか？
Siemens Energy（旧CD-adapco）は海事CFDに長い歴史がある。STAR-CCM+のoverset meshと自動メッシングはプロペラ解析で使いやすい。また、Abaqus連携による強連成FSIも可能だ。
複合材プロペラの場合はどうですか？
異方性積層構造をモデル化できる構造ソルバーが必要だ。Abaqusの積層シェル要素（S8R等）やNastranのPCOMP/PCOMPG定義が適している。Ansys Mechanical ACP（Composite PrepPost）も積層定義に便利だ。
造船大国の解析ツール事情——韓国・中国・日本の使い分け
世界の造船シェアトップを争う韓国・現代重工業や中国・中国船舶集団は、プロペラキャビテーション解析にどのツールを使っているのか。韓国勢はStar-CCM+とANSYS Fluentを使い分け、水槽試験との相関取りに膨大な実績データを持つ。中国勢はオープンソースのOpenFOAMを改良した独自スタックも多い。日本の三菱・川崎・今治はFINE/Marineなど船舶特化ソルバーの信頼性を重視する傾向がある。面白いことに、同じキャビテーションモデル（Schnerr-Sauer法）を使っていても、メッシュ戦略や時間刻みの取り方で結果が10〜20%変わることがあり、各社はそのノウハウを門外不出にしている。
選定で最も重要な3つの問い
- 「何を解くか」：プロペラキャビテーションFSIに必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
複合材プロペラの適応変形
複合材プロペラが注目されているのはなぜですか？
CFRP（炭素繊維強化プラスチック）の異方性を利用して、荷重に応じて翼のピッチ角が変化するBend-Twist Coupling（BTC）を実現できる。これにより、広い運転範囲で高効率を維持するadaptive propellerが可能になる。
積層方向とBTC係数の関係は、
$D_{16}$ が非ゼロになる積層配列（例：[+45/-45]のアンバランス積層）でBTCが発生する。
キャビテーションエロージョン予測
キャビテーションによる翼面の侵食をシミュレーションで予測できますか？
キャビティ崩壊時の衝撃圧力を推定し、材料の耐エロージョン性と比較する。Patela（2016）のエロージョン強度指標を使うアプローチが広まりつつある。FSIを通じて翼の変形がキャビテーションパターンを変え、それがエロージョン分布に影響するという間接的な連成も重要だ。
騒音予測との統合
キャビテーション騒音の予測もFSIと関係しますか？
キャビテーションは船舶の水中放射騒音の主要源だ。FW-H音響アナロジーをCFD結果に適用して遠方場騒音を予測する。翼の変形がキャビテーション体積変動を変え、それが騒音スペクトルに影響する。Ansys FluentとSTAR-CCM+にはFW-Hソルバーが組み込まれている。
量子コンピュータとキャビテーション——次世代シミュレーションの最前線
キャビテーション-FSI連成の計算コストは並外れて高い。気泡1個の崩壊を忠実に追うDNS計算は、現代のスパコンでもプロペラ1回転に数週間かかることがある。そこで注目されているのが量子コンピュータを使った格子ボルツマン法の加速だ。IBMとDLRの共同研究では、量子回路で流体の位相空間を表現し、古典計算比で100倍の高速化を理論的に示した（2023年）。実用化はまだ先だが、「Navier-Stokes方程式を量子ビットで解く」という発想自体が、FSI解析の地平を根本から変えようとしている。
トラブルシューティング
キャビテーションが発生しない
計算しても実験で観測されるキャビテーションが再現できません。
チェックポイントを整理しよう。
| 確認項目 | 詳細 |
|---|---|
| 蒸気圧の設定 | 温度に対応した正しい$p_v$を設定しているか |
| 参照圧力 | operating pressureが正しく設定されているか |
| 乱流モデル | SST k-omegaを推奨。k-epsilonは粘度を過大評価しキャビテーションを抑制 |
| メッシュ密度 | 翼面の低圧部が十分に解像されているか |
| 時間刻み | 定常計算ではキャビテーションが安定しないことがある |
推力・トルクの不一致
$K_T$ や $K_Q$ が実験値から大きくずれます。
FSI連成が収束しない場合の対処法は？
回転体FSIでは遠心力による初期変形が大きいため、まず構造単体で遠心力解析を実行し、変形後の形状を初期形状としてFSI計算を開始する。また、連成の緩和係数を小さめ（0.2〜0.5）に設定して安定させてから徐々に上げるのが実用的だ。
「解析では問題なし」なのに船が振動する——よくある落とし穴
実務でよく聞く失敗談がある。「キャビテーション-FSI解析で船体振動は問題なし」と出たのに、就航後に船尾甲板が激しく振動し、船員が不快を訴えた事例だ。原因を調べると、解析では均一流入速度を仮定していたのに対し、実際は不均一な伴流（wake）が影響していた。プロペラの前方には船体があるため、流速は円周方向で最大30%もばらつく。この不均一性がキャビテーションの発生タイミングのずれを生み、周期的な圧力脈動を増幅させていた。「均一流入仮定」は計算を簡単にするが、振動問題では特に注意が必要だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——プロペラキャビテーションFSIの問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告