旗・膜のFSI解析 — 数値解法と実装
数値手法の選択
旗のFSI解析ではどんな手法が使われますか？
IB法だと膜の厚みをゼロとして扱うんですか？
Peskin型IB法ではラグランジュ構造を1次元低い多様体として扱う。2Dの旗は1Dのファイバー、3Dの膜は2D面として埋め込む。力のスプレッドにDiracデルタ関数の正則化版を使う。
時間積分と安定性
膜が軽い場合、時間積分で問題が出ますか？
質量比が小さい系では陽的時間積分のCFL条件が非常に厳しくなる。半陰的スキーム（流体を陰的、構造を陽的に解く）や、完全陰的スキームが必要だ。
Causin型の安定性条件によると、弱連成の安定性は $\rho_s h / (\rho_f \Delta x)$ に依存する。この比が小さいと弱連成は不安定だ。
リバティ船の溶接割れ——連成問題の教訓
第二次世界大戦中、アメリカは「リバティ船」を溶接で大量生産し、戦争の物流を支えました。しかし約1,500隻のうち約400隻に船体の亀裂が発生。原因は溶接残留応力と低温脆性の連成——溶接時の急激な温度変化が残留応力を生み、北大西洋の冷たい海水で鋼材が脆くなり、亀裂が伝播したのです。現代の溶接シミュレーションは、この「温度→残留応力→破壊」の連鎖を予測できます。
離散化手法の詳細解説
空間離散化における手法選択が数値精度・安定性・計算コストに与える影響を詳述する。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
マトリクスソルバーの選定指針
問題規模と特性に応じた最適なソルバー選択のガイドライン。
| ソルバー種別 | 詳細・推奨条件 |
|---|---|
| ブロック対角前処理 | 各物理場の前処理を対角ブロックとして使用。実装が容易でスケーラビリティが良い。 |
| フィールド分割法 | 速度-圧力分割（流体）、変位-温度分割（熱-構造）等。物理的な意味に基づくブロック分割が安定性に寄与。 |
時間積分法と収束判定
ソルバー内部の制御パラメータと収束判定基準について記述する。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
数値解法の直感的理解
連成ソルバーのイメージ
モノリシック手法は「バンドの全員が同じ楽譜を見て同時に演奏する」——完璧な同期が取れるが、楽譜（連立方程式系）が巨大で複雑。パーティション手法は「各パートが別々に練習して合わせる」——個々の練習（ソルバー）は簡単だが、合わせ（データ交換）のタイミングと精度が課題。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
連成解析の安定性やデータ転写の精度は、マルチフィジックスの永続的な課題です。 — Project NovaSolverはこの課題に正面から取り組んでいます。
CAEの未来を、実務者と共に考える
Project NovaSolverは、旗・膜のFSI解析における実務課題の本質に向き合い、エンジニアリングの現場を支える道具づくりを目指す研究開発プロジェクトです。
プロジェクトの最新情報を見る →