スロッシング-構造連成
スロッシング-構造連成の理論基礎
スロッシングの概要
スロッシングってどういう現象ですか？
容器内の液体が外部励振（地震、波浪、車両加速等）によって激しく揺動する現象だ。LNGタンカーの貨物タンク、原子力発電所の使用済燃料プール、ロケットの推進剤タンクなどで重要になる。
支配方程式
スロッシングの数理モデルはどうなっていますか？
線形理論では速度ポテンシャル $\phi$ で記述する。自由表面のスロッシング固有振動数は矩形タンクの場合、
$L$ はタンク長さ、$h$ は液面高さ、$n$ はモード次数だ。
非線形スロッシング（大振幅、砕波、衝撃圧力を伴う）ではNavier-Stokes方程式をVOF法で解く。
タンク壁への衝撃圧力（スロッシングインパクト）は局所的に数MPaに達する。
構造との連成はどう扱うんですか？
タンク壁の弾性変形が液面挙動に影響する場合にFSIが必要になる。特に膜型LNGタンク（Mark III, NO96）では薄い波板構造がスロッシング衝撃で変形し、それが圧力分布を変える。
スロッシングの「固有振動数」——タンクの形が全てを決める
スロッシングで一番怖いのは「共振」だ。タンク内液体の固有振動数は、長方形タンクなら f₁ ≈ (1/2π)√(πg/L・tanh(πh/L)) で計算できる（Lはタンク長さ、hは液深）。液体が半分程度入った状態で、この固有振動数と船体のローリング周期が一致すると振幅が爆発的に増大する。1960年代の初期LNGタンカー設計では、この共振条件を見落としてタンク壁に繰り返し衝撃荷重がかかり、断熱材が剥がれる事故が多発した。以来「フルかほぼ空かの積載状態を保ち、共振レンジを通過させない」が業界のオペレーショナルルールとなっている。
スロッシング-構造連成の数値計算手法
数値手法
スロッシングFSIにはどんな手法が使われますか？
MPS法は日本発の手法なんですか？
越塚誠一教授（東京大学）が開発したMoving Particle Semi-implicit法だ。Prometech社のParticleworksとして商用化されている。自由表面の大変形や飛沫を安定して追跡でき、スロッシング解析で多くの実績がある。
衝撃圧力の評価
スロッシング衝撃圧力をどう評価するんですか？
衝撃圧力にはair pocket型（圧縮空気クッション）とflip-through型（直撃型）がある。
- Air pocket型: 液面と壁の間に空気が閉じ込められ、圧縮される。振動的な圧力波形。ピーク圧力は低いが持続時間が長い
- Flip-through型: 液面が壁面に沿って上昇し直撃する。非常に高いピーク圧力だが持続時間が短い
圧縮性を考慮する場合、空気のポリトロープ過程を含む二相流解析が必要だ。
SPHとVOF——スロッシング解析の「二大流派」
スロッシング数値解析では格子を使うVOF（Volume of Fluid）法と格子を使わないSPH（粒子法）の2つが長年競合している。VOFはOpenFOAMやSTAR-CCM+で実績が豊富で、気液界面の捉え方が直感的だ。一方SPHは液体の大変形（タンク壁への衝突やしぶきの飛散）を自然に扱えるのが強み。造船業界ではIHIや三星重工がSPHをスロッシングの衝撃圧力解析に積極活用している。最近はSPH計算をGPUで並列化し、100万粒子規模のシミュレーションを数時間で完了する環境も整ってきた。どちらが「正しい」かではなく、解析目的に応じた使い分けが実務の答えだ。
スロッシング-構造連成の実務適用
LNGタンクの解析手順
LNGタンカーのスロッシング解析はどう進めるんですか？
1. 耐航性解析で船体運動を計算（6自由度運動の時刻歴）
2. タンク形状と液位の設定（部分充填率が重要。20〜80%が危険領域）
3. CFD（VOF/SPH/MPS）でスロッシング解析
4. 壁面圧力の統計処理（短期・長期極値分布）
5. 構造応答解析（衝撃圧力をFEモデルに入力）
充填率が中間だと危険なのはなぜですか？
液面が低すぎると液体の質量が小さく、高すぎると液体の動きが制限される。40〜70%の中間充填率でスロッシングが最も激しくなる。IGC Code（国際ガス運搬船コード）では危険充填率範囲での航行を制限している。
衝撃圧力の統計処理
衝撃圧力のばらつきはどう扱うんですか？
スロッシング衝撃圧力は確率論的に非常に大きなばらつきを示す。同一条件でも衝撃ごとにピーク圧力が10倍以上変動する。
Gumbel分布またはWeibull分布で極値統計処理を行い、設計寿命に対する最大圧力を推定する。短期3時間の解析で数百回の衝撃イベントを取得し、極値分布のパラメータを推定する。
LNGタンカーのNo.96型タンク——スロッシング対策の結晶
世界のLNG輸送の主流を占めるメンブレン型タンクのうち、GTT社のNo.96型は複雑な断熱構造でスロッシング荷重を吸収する設計になっている。タンク容量14万〜17万立方メートルの大型LNGタンカーでは、スロッシング解析に数百ケースのシミュレーションが必要で、波浪条件・積載率・航行速度の組み合わせを総当たりで検証する。日本の大手造船所（三菱重工・川崎重工）では専用のスロッシング試験装置（6自由度振動台に1/50スケールモデルを乗せる）を持ち、CFD結果の検証に使っている。一隻のLNGタンカー建造に費やすスロッシング解析コストは数千万円規模にのぼる。
スロッシング-構造連成のソフトウェア比較
ツール比較
スロッシング解析に使えるツールは？
ParticleworksがLNG業界で使われている理由は？
MPS法は自由表面の大変形や砕波を安定して追跡できる。GUIでタンク動揺条件の設定が容易で、壁面圧力の自動統計処理機能を備えている。三菱重工やサムスン重工などの造船所での採用実績が豊富だ。
FSI連成まで含めるとどうですか？
LS-DYNAのSPH-FEM連成が最も直接的だ。ParticleworksやFLOW-3Dの結果を圧力データとしてエクスポートし、FEMソルバー（Nastran, Abaqus等）にmapped loadとして入力するone-way連成も実務では多い。
スロッシング解析ソフトのシェア争い——専用ツールvs汎用ツール
スロッシング解析のツール市場は面白い構造になっている。ANSYS Fluent・STAR-CCM+などの汎用CFDが大きなシェアを持つ一方、DNV Wasimのような船舶専用ツールも根強い支持を集める。海洋・石油メジャー（Shell、Total、BP）はSTAR-CCM+を採用するケースが多く、LNGタンカー専業の造船所はDNV-GT（GTTの技術評価ツール）との組み合わせを使うことが多い。国内では、国土交通省の「スロッシング評価手引き」がANSYSベースの解析事例を多く掲載しているため、行政提出案件はANSYS系が有利な状況だ。ただし、GTTの特許管理が厳しいため独自設計タンクの普及は限られている。
スロッシング-構造連成の先端研究
極低温FSI
LNGは-162°Cですよね。極低温の効果は考慮するんですか？
LNGの物性（密度、粘度、表面張力）は温度依存だ。さらに、LNGの蒸発（BOG: Boil-Off Gas）による気相の存在がスロッシングに影響する。気液二相のスロッシングでは、気相の圧縮性がair pocket型衝撃圧力を変化させる。
断熱材（ポリウレタンフォーム等）の極低温での脆性破壊もスロッシング衝撃との組み合わせで評価すべき問題だ。GTT社のMark IIIやNO96タンクでは、断熱材の衝撃応答が設計の鍵となっている。
比較スロッシング実験
実験との対応はどうなっていますか？
GTT社の大型スロッシング試験装置（6自由度動揺台+縮尺タンク）が業界標準の検証データを提供している。また、ISOPE（International Society of Offshore and Polar Engineers）のスロッシングベンチマークワークショップで複数コードの比較が定期的に行われている。
制振デバイス
スロッシングを抑制する方法はあるんですか？
バッフル（仕切り板）やTLD（Tuned Liquid Damper）が代表的だ。バッフルの形状最適化にCFDが活用されている。
$A_b$ はバッフル面積、$C_D$ は抗力係数だ。バッフルの開口率と位置の最適化がスロッシング抑制効果を左右する。
ロケット燃料タンクのスロッシング——アポロ計画の隠れた難問
スロッシングは海洋だけの問題ではない。アポロ計画でも液体酸素・液体水素タンク内の燃料スロッシングが飛行安定性を脅かした。ロケットの姿勢制御系が燃料の揺れに反応して過修正を繰り返し、制御不安定に陥るリスクがある。NASAはこれを「プロペラント・スロッシング」と呼び、タンク内にバッフル（仕切り板）を設けてスロッシングを抑制した。現在のSpaceXファルコン9でも大型推進剤タンクのスロッシング抑制は重要な設計課題で、CFD-FSI連成シミュレーションが姿勢制御アルゴリズムの開発と並行して行われている。
スロッシング-構造連成のトラブル対応
衝撃圧力の再現性問題
同じ条件で計算しても衝撃圧力が毎回違うんですが。
これはスロッシング衝撃の本質的な特性だ。微小な初期条件の違いが砕波の形状を変え、衝撃圧力に大きなばらつきをもたらす。対策は、
- 多数回の計算を実行して統計的に評価する
- 空間平均圧力（パッチ圧力）を使う
- Weibull/Gumbel分布で極値統計処理を行う
メッシュ密度と衝撃圧力の関係は？
スラミングと同様、ピーク圧力はメッシュ依存性が強い。力の積分値（衝撃力）の方がメッシュに対して安定する。構造応答を最終的な評価指標にすれば、構造の慣性フィルタ効果で高周波成分が減衰する。
液面追跡の精度
VOF法で液面が拡散してしまいます。
| 対策 | 詳細 |
|---|---|
| 界面圧縮スキーム | OpenFOAMのinterFoamのcAlpha=1で界面圧縮を有効化 |
| AMR | 液面近傍のメッシュを動的細分化 |
| HRIC/CICSAM | 高精度VOFスキームの選択 |
| Level Set + VOF | CLSVOF法で界面の鮮鋭性を改善 |
計算が非常に長時間かかる場合の対処法は？
スロッシング解析は長時間計算（数千秒以上の物理時間）が必要だ。GPUソルバー（ParticleworksはマルチGPU対応）やAMRによる動的メッシュ最適化で高速化を図る。また、等価不規則励振を短時間の規則励振に変換するequivalent regular wave法で計算時間を削減する手法もある。
スロッシング解析の「圧力ピーク再現性問題」
スロッシング衝撃圧力のCFD解析でよく起きる問題が「実験値に対してピーク圧力が大きくばらつく」という現象だ。同じ実験条件・同じソルバーで繰り返し計算しても、衝撃圧力のピーク値が±50%以上変動することがある。原因の一つはカオス的な液面挙動——初期条件のわずかな差が液面形状を全く異なる状態に引き込む「バタフライ効果」が働く。もう一つは計算格子のサイズで、格子が粗いほどピーク圧力が過小評価される傾向がある。実務対応として「多数の解析ケースの統計的分布で評価する」「規則波ではなく不規則波を入力する」などのアプローチが標準化されつつある。
関連トピック
なった
詳しく
報告