[論文レビュー] Three-Dimensional Simulations of Deep-Water Breaking Waves

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

Ocean Waves and Remote Sensing参考文献 7被引用数 16

ひとこと要約

本稿では、移動計算フレーム内での完全非線形進行波を生成するために、新規の法線応力駆動技術を用いた体積関数法を用いた三次元数値シミュレーションフレームワークを提示する。主な発見は、弱いから非常に強いプラウンピングに至るまでの明確な破壊状態の同定、空気巻き込みとエネルギー散逸の定量的評価（強いイベントではピークエネルギーの50％以上が散逸）、およびレインズド平均化に比べてFavre平均化が自由表面追跡の精度を向上させることの実証である。

ABSTRACT

The formulation of a canonical deep-water breaking wave problem is introduced, and the results of a set of three-dimensional numerical simulations for deep-water breaking waves are presented. In this paper fully nonlinear progressive waves are generated by applying a normal stress to the free surface. Precise control of the forcing allows for a systematic study of four types of deep-water breaking waves, characterized herein as weak plunging, plunging, strong plunging, and very strong plunging.

研究の動機と目的

• 深海の破壊波を三次元的に研究するための標準的フレームワークの構築を目的とする。
• 弱いプラウンピングから非常に強いプラウンピングに至るまでのさまざまな破壊波タイプの体系的シミュレーションを可能とし、駆動力に対して正確な制御を可能とすることを目的とする。
• 破壊波における空気巻き込みおよびエネルギー散逸メカニズムの定量的評価を目的とする。
• 自由表面流れの統計的解析において、レインズド平均化とFavre平均化を比較することを目的とする。
• NFAコードを用いて50億～100億の未知数を有する大規模シミュレーションの実現可能性を評価することを目的とする。

提案手法

• 自由表面に法線応力駆動技術を適用し、移動計算フレーム内での完全非線形進行波の生成を実現する。
• 空気-水界面の追跡には体積関数法（VOF）を用い、各セル内の体積分率αが流体の占める割合を表す。
• 圧力および界面張力効果を組み込んだ保存則形式の運動量方程式を用いて、非圧縮性ナビエ＝ストークス方程式を解く。
• 有限体積法を用い、コロケートグリッドと圧力-速度結合アルゴリズムを採用する。
• グローバルな集合的処理を置き換えるために、1次元のサブ通信器を用いた変更済みMPI通信戦略を採用し、数千プロセッサへの効率的拡張を可能にする。
• 速度場および体積分率場にFavre平均化とレインズド平均化を適用し、自由表面ダイナミクスの捕捉における統計的性能を比較する。

実験結果

リサーチクエスチョン

• RQ1 法線応力駆動技術は、深海において完全非線形進行波をどのように制御的に生成可能にするか？
• RQ2 弱いプラウンピングから非常に強いプラウンピングに至るまでのさまざまな破壊波タイプに伴う、三次元的な流れ構造および渦度分布は何か？
• RQ3 破壊過程でエネルギーはどの程度散逸するのか？また、破壊の強度に応じてその散逸度はどのように変化するか？
• RQ4 さまざまな破壊状態において、自由表面下での空気脱ガスの速度はどのように比較されるか？
• RQ5 統計的平均化手法として、レインズド平均化とFavre平均化のどちらが、シミュレーションにおける真の自由表面位置とより良好に一致するか？

主な発見

• 強いおよび非常に強いプラウンピング破壊イベントでは、ピークエネルギーの50％以上が散逸しており、激しい破壊時に顕著なエネルギー損失が生じることを示している。
• 空気巻き込み直後に、主なプラウンピングイベントで流れ方向の強い渦度が形成されることが、三次元の渦度等値線図から明らかになった。
• 強いおよび非常に強いプラウンピングケースでは、初期の破壊イベントから長期間にわたり、顕著な量の空気が自由表面下に巻き込まれたまま残存している。
• 空気脱ガスの速度は、すべてのケースで線形的かつ同一であり、破壊強度に依存しない一貫した物理的メカニズムを示唆している。
• 速度場および体積分率場のFavre平均化は、レインズド平均化に比べて真の自由表面位置との一致度が高く、界面ダイナミクスの捉え込みにおいて優れている。
• 50億～100億の未知数を有するシミュレーションは現在、次世代ハードウェア上で最適条件の下でNFAコードが1時間に約100イテレーションを達成できるため、実現可能となった。