渦励振（VIV）解析 — トラブルシューティングガイド
連成計算が発散する場合
VIVの連成計算を回したら発散してしまいました。何が原因でしょうか？
VIV-FSI解析の発散は主に以下の原因で起きる。
| 症状 | 原因 | 対策 |
|---|---|---|
| 1ステップ目で発散 | 界面マッピングの不整合 | 流体・構造の界面メッシュを確認。ノードの対応関係を検証 |
| 数ステップ後に発散 | added mass不安定性 | 強連成に切り替え。緩和係数を0.1〜0.3に下げる |
| 振幅成長後に発散 | メッシュ変形過大 | リメッシング頻度を上げる。overset meshを検討 |
| 特定流速でのみ発散 | ロックイン領域での共振 | 時間刻みを細かくし、サブイテレーション回数を増やす |
added mass不安定性って具体的にはどういうことですか？
流体の付加質量が構造質量と同程度以上（$m^* \lesssim 1$）になると、弱連成では人工的なエネルギーが界面に蓄積して発散する。Causin et al.（2005）の解析で理論的に示されている。対策は強連成＋緩和、またはRobin-Neumann条件の適用だ。
Strouhal数が文献値と合わない
計算で得られたStrouhal数が0.2から大きくずれるんですが。
チェックポイントを整理しよう。
- ブロッケージ比: 計算領域幅に対する円柱直径の比が5%以上だと壁の影響でStが変わる。$D/H < 0.03$ を確保する
- 時間分解能: FFTのサンプリングが不足すると周波数の同定精度が落ちる。Welch法で窓関数を適用し、十分な時間データを使う
- 境界条件: 出口境界が近すぎると後流が拘束される。$30D$ 以上離す
- 2D vs 3D: Re > 200 では3D効果が顕著。2D計算はStを過大評価する傾向がある
疲労評価への接続
VIV解析の結果を疲労寿命評価に使うにはどうすればいいですか？
時刻歴応力をレインフロー法でカウントし、S-N曲線とMinerの線形累積損傷則で疲労損傷度を算出する。
DNV-RP-C203のS-N曲線を使い、$D_{fatigue} < 1/DFF$（DFF: Design Fatigue Factor、通常3〜10）を満足させる。応力集中係数（SCF）の設定が結果を大きく左右するから注意が必要だ。
CFD-FSIの結果をそのまま疲労評価に流せるパイプラインがあると便利ですね。
AnsysではFluent→Mechanical→nCode DesignLifeの連携が可能だ。STAR-CCM+もFE-Safeとの連携インターフェースを持っている。ただし、CFDの時刻歴は非常に長いデータになるから、サイクルカウントの効率化が実務上の課題だよ。
心臓シミュレーション——究極のFSI問題
人間の心臓は1日に約10万回拍動し、血液を全身に送り出します。この過程は流体（血液）-構造（心筋・弁）-電気（刺激伝導系）の3場連成問題。心臓のデジタルツインの構築は連成解析の「聖杯」と呼ばれ、世界中の研究者が挑戦しています。実現すれば、手術のシミュレーションや薬の効果予測が患者ごとにカスタマイズできるようになります。
トラブル解決の考え方
デバッグのイメージ
連成解析のトラブルシューティングは「チームプレーの問題解決」に似ている。まず「どのチーム（物理場）に問題があるか」を切り分け、次に「チーム間の連携（データ転写）に問題がないか」を確認する。各物理場を単独で動かして問題がなければ、連成の設定が原因。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——渦励振（VIV）解析の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。
渦励振（VIV）解析の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →