[論文レビュー] CN-Stream: Open-source library for nonlinear regular waves using stream function theory

Guillaume Ducrozet, Benjamin Bouscasse+3 more|arXiv (Cornell University)|Jan 29, 2019

Ocean Waves and Remote Sensing被引用数 10

ひとこと要約

CN-Stream は、改良されたストリーム関数理論を用いて、深海または浅海の非線形正規波を高精度に計算するオープンソースの Fortran ライブラリである。波動運動および自由表面プロファイルの高精度なシミュレーションを可能にし、古典的な Rienecker と Fenton の手法を拡張して、表面の傾きと速度ポテンシャルのための別々のモード数を可能にするとともに、目的の精度を満たすためにモード数を自動で最適化する。これにより、破壊限界に近い波への適用範囲が拡大される。

ABSTRACT

CN-Stream is a library for the computation of nonlinear regular ocean waves. The library is developed in order to be easily integrated with wave generation models in CFD solvers. It is based on the stream function theory and provides significant improvements regarding the applicability of the method for waves close to breaking (in deep or shallow water) compared to the classical implementation of Rienecker and Fenton [26]. The complete description of the wave field is available, including the free-surface evolution and the wave kinematics in the fluid domain. It is released as open-source, developed and distributed under the terms of GPL v3.

研究の動機と目的

• 任意の定常水深における非線形正規波を高精度にシミュレートできるオープンソースツールの開発。
• 深海または浅海における急波・破壊限界に近い波をモデル化する際の、古典的ストリーム関数法の限界を克服すること。
• 波構造物相互作用のシミュレーションに向けた CFD ソルバへのシームレスな統合を可能にすること。
• 最適な精度と計算効率を実現するための自動モード選択を提供すること。
• 物理的解析およびモデル検証のための詳細な波動運動および圧力場出力を提供すること。

提案手法

• 有限フーリエ級数展開を用いて、完全非線形ポテンシャル流れ問題をストリーム関数理論により解く。
• 自由表面境界条件を、ニュートン反復法を用いて解く非線形代数方程式系に還元する。
• 数値的精度を向上させるために、表面傾きのためのモード数（N1）と速度ポテンシャルのためのモード数（N2）を別々に設定する。
• ユーザーが定義した収束基準を満たすために、最適なコロケーション点数（モード数）を自動で決定する。
• 相対誤差と絶対誤差のしきい値（eps_inc, err_max）を組み合わせたハイブリッド収束戦略を採用し、収束性を確保する。
• 可視化および結合のための、モード振幅、速度場および圧力場、自由表面プロファイルを含む包括的な出力ファイルを生成する。

実験結果

リサーチクエスチョン

• RQ1 ストリーム関数理論をどのように拡張すれば、深海および浅海の両方において破壊限界に近い非線形正規波を高精度にモデル化できるか？
• RQ2 古典的な Rienecker と Fenton の定式化をどのように改善すれば、より高い波の傾きおよび広い水深範囲への適用範囲を拡大できるか？
• RQ3 自動モード選択は、波動シミュレーションにおける精度を著しく向上させるとともに、計算効率を維持できるか？
• RQ4 非線形波の運動および圧力場は、さまざまな波の傾きと水深において、線形予測とどのように異なるか？
• RQ5 CN-Stream は、CFD の検証および波構造物相互作用モデリングの信頼できる基準解として、どの程度の範囲で利用可能か？

主な発見

• CN-Stream は、波の傾きおよび相対的水深の広い範囲において、破壊に近い状態を含め、高精度で非線形正規波を計算できることを確認した。
• 表面傾きと速度ポテンシャルのための別々のモード数の使用により、標準的な実装と比較して数値的安定性および精度が向上した。
• 自動モード選択により、最小限のユーザー入力で最適な収束が保証され、発散のリスクが低減され、全体の堅牢性が向上した。
• 出力は、可視化および結合に適した形式（Tecplot など）で提供され、正確な波動運動、圧力場、自由表面プロファイルを生成した。
• CN-Stream は、ICARE、WCCH、OpenFOAM などの CFD ソルバに成功裏に統合され、波構造物相互作用のシミュレーションに使用された。
• 1回のシミュレーションあたり数秒で収束を達成しており、長期的な波動伝播研究およびモデル検証に適した性能を示した。