[論文レビュー] A Numerical Formulation for Simulating Free-Surface Hydrodynamics

Thomas T. O’Shea, Kyle A. Brucker+2 more|arXiv (Cornell University)|Oct 7, 2014

Fluid Dynamics Simulations and Interactions参考文献 7被引用数 23

ひとこと要約

本稿では、船の周囲における自由表面流動、特にトランサム・ステン流動をシミュレートするために、インマージドボディ法と体積関数法（VOF）を組み合わせたカルテジアングリッド有限体積法を提案する。この手法は波の崩壊、スプレー、空気の混入を正確に捉えることができ、数値予測は複数の牽引速度におけるドラッグおよび波の断面測定値と強く一致している。

ABSTRACT

Cartesian-grid methods in combination with immersed-body and volume-of-fluid methods are ideally suited for simulating breaking waves around ships. A surface panelization of the ship hull is used as input to impose body-boundary conditions on a three-dimensional cartesian grid. The volume-of-fluid portion of the numerical algorithm is used to capture the free-surface interface, including the breaking of waves. The numerical scheme is implemented on a parallel computer. Recent improvements to the numerical scheme are discussed, including implementation of a new multigrid procedure and conversion to MPI communication. Numerical predictions are compared to laboratory measurements of a towed transom-stern model.

研究の動機と目的

• 波の崩壊や空気の混入を含む複雑な自由表面流れをシミュレートできる、堅牢で即戦可能なCFDフレームワークの開発。
• ボディ・フィットドグリッドの必要性を排除するため、インマージドボディ境界条件を用いたカルテジアングリッドを採用。
• 新しいマルチグリッドソルバーとMPIベースの通信を導入することで、大規模シミュレーションにおける数値的効率性と拡張性を向上。
• 異なる牽引速度におけるトランサム・ステンモデルの詳細な実験測定値と照合して、数値予測の妥当性を検証。
• 数十億グリッドポイントを含む大規模な海軍流体力学的問題の高解像度シミュレーションを可能に。

提案手法

• 圧縮性のないナビエ＝ストークス方程式を、分数ステップ時間積分スキームを用いて、カルテジアンステガラードグリッドで解く。
• 体積関数法（VOF）を用いて、セル単位の体積分率関数により自由表面界面を追跡し、鋭い液体・空気界面を表現するための滑らか化されたステップ関数を用いる。
• 表面パネル化を入力として用い、ボディ力法により船体の不透過性およびスリップなし境界条件をインマージドボディ境界条件として強制する。
• 分数ステップ法における圧力ポアソン方程式の収束を加速するために、新しいマルチグリッドソルバーを実装する。
• 高性能コンピューティングシステムにおける並列スケーラビリティを向上させるために、MPIベースの通信ルーチンを統合する。
• 乱流および粘性応力のための準解像度モデルを用い、SGS応力テンソルは動的手続きによりモデル化する。

実験結果

リサーチクエスチョン

• RQ1 VOFおよびインマージドボディ法を用いたカルテジアングリッド法は、船モデル周囲の波の崩壊およびトランサム・ステン流動を正確にシミュレートできるか？
• RQ2 異なる牽引速度におけるドラッグおよび波の断面プロファイルの数値予測は、実験測定値とどの程度一致するか？
• RQ3 新しいマルチグリッドソルバーは、自由表面流れの大規模CFDシミュレーションにおける収束速度をどの程度向上させるか？
• RQ4 プランジング波の崩壊に対する自由表面の高さおよび水粒子速度の数値予測は、実験データとどの程度一致するか？
• RQ5 高周波数の表面スペクトル特徴を解像する際の限界は何か。これはアリゼーションやモデル簡略化とどのように関連しているか？

主な発見

• トランサム・ステンモデルのドラッグ力に対する数値予測は、DTMBの実験測定値と4つの牽引速度すべてで優れた一致を示した。
• NFAコードによる波の断面プロファイルの予測は、全4速度および全横断位置で実験データと非常に良く一致しており、特に船尾領域で一貫した精度を示した。
• プランジング波の崩壊を模擬したシミュレーションにおける自由表面の高さは、形状および大きさの両面で実験測定値とよく一致した。
• このモデルは、波の崩壊、スプレー形成、空気の混入といった複雑な流れの特徴を、特に中程度の速度（7–8ノット）でうまく捉えた。
• R/Vアテナの後方のロースタール・テイル領域では、平均的な自由表面高さがよく一致したが、高周波数成分ではずれが生じており、これはアリゼーションや数値モデルに付属部品がないことが要因である可能性がある。
• 新しいマルチグリッドソルバーおよびMPI通信ルーチンの導入により、ソルバーの収束性と並列スケーラビリティが顕著に向上し、10億セルに近いグリッド解像度のシミュレーションが可能になった。