スロッシング-構造連成
理論と物理
スロッシングの概要
スロッシングってどういう現象ですか？
容器内の液体が外部励振（地震、波浪、車両加速等）によって激しく揺動する現象だ。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タンカー設計では、この共振条件を見落としてタンク壁に繰り返し衝撃荷重がかかり、断熱材が剥がれる事故が多発した。以来「フルかほぼ空かの積載状態を保ち、共振レンジを通過させない」が業界のオペレーショナルルールとなっている。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\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にはどんな手法が使われますか？
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万粒子規模のシミュレーションを数時間で完了する環境も整ってきた。どちらが「正しい」かではなく、解析目的に応じた使い分けが実務の答えだ。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
実践ガイド
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タンカー建造に費やすスロッシング解析コストは数千万円規模にのぼる。
解析フローのたとえ
風船を膨らませたことがありますか？ あの瞬間、実は高度な流体-構造連成が起きています。内部の空気圧（流体）がゴム壁（構造）を押し広げ→広がった壁が内部の圧力分布を変え→変わった圧力がさらに壁を変形させる…このキャッチボールを計算ステップごとに繰り返すのがFSI解析です。
初心者が陥りやすい落とし穴
「片方向連成で十分でしょ？」——この判断ミスが連成解析で最も危険です。構造の変形が微小なら確かに片方向で足りますが、心臓弁の開閉のように変形が流路を大きく変える場合、片方向では全く話になりません。目安は「変形量が代表長さの1%を超えるか」。超えるなら双方向連成は必須です。片方向で済ませてしまった場合、結果が「もっともらしいけど実は大間違い」になる——これが最も怖いパターンです。
境界条件の考え方
連成界面のデータ交換は「国境の出入国管理」と同じです。各国（物理場）には独自の法律（支配方程式）がありますが、国境（界面）で人や物（力・温度・変位）のやり取りを正確に管理しないと、両国の経済（エネルギーバランス）が崩壊します。メッシュが一致していない場合の補間は「通訳」のようなもの——誤訳（補間誤差）が小さいほど良い結果が得られます。
ソフトウェア比較
ツール比較
スロッシング解析に使えるツールは？
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の特許管理が厳しいため独自設計タンクの普及は限られている。
選定で最も重要な3つの問い
- 「何を解くか」：スロッシング-構造連成に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
極低温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%以上変動することがある。原因の一つはカオス的な液面挙動——初期条件のわずかな差が液面形状を全く異なる状態に引き込む「バタフライ効果」が働く。もう一つは計算格子のサイズで、格子が粗いほどピーク圧力が過小評価される傾向がある。実務対応として「多数の解析ケースの統計的分布で評価する」「規則波ではなく不規則波を入力する」などのアプローチが標準化されつつある。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——スロッシング-構造連成の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
関連トピック
なった
詳しく
報告