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