[論文レビュー] Semi-explicit solutions to the water-wave dispersion relation and their role in the nonlinear Hamiltonian coupled-mode theory

Theodosios K. Papathanasiou, Christos Papoutsellis+1 more|arXiv (Cornell University)|Feb 22, 2018

Coastal and Marine Dynamics参考文献 85被引用数 6

ひとこと要約

本稿では、特に減衰モードを含む水波分散関係の根を計算するための高精度で準明示的な解析的公式を提示している。これらはハイブリッドPicard-Householder反復スキームに基づくもので、3回以内の反復で機械精度の根の特定が可能である。この公式により、任意の水深における完全非線形水波シミュレーションに向けたハミルトニアン結合モード理論（HCMT）の効率性と安定性が著しく向上する。

ABSTRACT

The Hamiltonian coupled-mode theory (HCMT), recently derived by Athanassoulis and Papoutsellis [1], provides an efficient new approach for solving fully nonlinear water-wave problems over arbitrary bathymetry. In HCMT, heavy use is made of the roots of a local, water-wave dispersion relation with varying parameter, which have to be calculated at every horizontal position and every time instant. Thus, fast and accurate calculation of these roots, valid for all possible values of the varying parameter, are of fundamental importance. In this paper, new, semi-explicit and highly accurate root-finding formulae are derived, especially for the roots corresponding to evanescent modes. The derivation is based on the successive application of a Picard-type iteration and the Householder's root finding method. Explicit approximate formulae of very good accuracy are obtained, which are adequate to support HCMT for many types of applications. In most demanding cases, e.g. very steep, deep-water waves, machine-accurate determination of the required roots is achieved by no more than three iterations, using the explicit forms as initial values. Exploiting this root-finding procedure in the HCMT, results in an efficient, numerical solver able to treat fully nonlinear water waves over arbitrary bathymetry. Applications to demanding nonlinear problems demonstrate the efficiency and the robustness of the present approach.

研究の動機と目的

• 空間的・時間的各点で水波分散関係を解く必要が生じるハミルトニアン結合モード理論（HCMT）における計算ボトルネックを解消すること。
• 初期推定値が不適切な深海条件において、ニュートン・ラプソン法などの標準的根の特定法が不安定・非効率になる問題を克服すること。
• あらゆる水深領域で機械精度を保証する最小限の反復回数で収束する、安定的かつ高次の反復法を開発すること。
• 波数の正確な閉形式初期近似を提供することで、複雑な水深分布上での完全非線形水波の安定的で長時間にわたるシミュレーションを可能にすること。
• 再利用可能で信頼性の高い根の特定公式を提供することで、HCMTおよび類似のマルチモード手法の広範な応用を促進すること。

提案手法

• まずPicard型反復により良好な初期推定値を生成し、次にハウスホルダーの高次根の特定法を適用する二段階プロセスを用いて、局所的水波分散関係の根の準明示的近似を導出する。
• 分散関係方程式にハウスホルダー反復を適用し、波数 $k_n$ に関して非線形であるが、周波数パラメータ $\mu(x,t)$ は空間的・時間的に変化する。
• すべての水深条件、特に減衰モードを含む、$k_n(\mathbf{x},t)$ の初期近似の明示的・閉形式表現を構築する。
• Picardに続くハウスホルダーの複合スキームは、ハウスホルダーの収束次数を保持するが、誤差定数が顕著に低減されることを示す。
• 得られた根の特定手順を、係数 $A$, $B$, $C$ がこれらの高精度な $k_n$ 値を用いて計算されるHCMTフレームワークに統合する。
• 急峻な水深変化上での孤立波の進化、波の散乱、波の崩壊といった完全非線形波問題におけるHCMTシミュレーションを通じて、手法の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 すべての水深領域、特に減衰モードを含む状況において、水波分散関係の根を効率的かつ高精度に計算する方法は何か？
• RQ2 ニュートン・ラプソン法が深海で失敗するような場合でも、収束性と高精度を保証する反復戦略は何か？
• RQ3 低次と高次反復法を組み合わせたハイブリッド手法は、標準的手法に比べて優れた収束特性を示せるか？
• RQ4 準明示的初期近似を用いることで、機械精度の根の特定に必要な反復回数はどの程度削減できるか？
• RQ5 提案手法の根の特定法は、完全非線形水波の長時間シミュレーションにおけるHCMTの安定性と効率性をどの程度向上させるか？

主な発見

• 提案された2次および3次準明示的手法は、3回以内の反復で波数計算において機械精度に達する。
• 理論的収束次数がより高いにもかかわらず、2次手法が3次手法よりも計算効率に優れている。
• 本手法はあらゆる水深条件で安定しており、初期推定値を慎重に選ばない限りニュートン・ラプソン法が発散する深海領域でも発散を回避する。
• $k_n$ の閉形式初期近似は、反復ソルバーの出発点として十分に精度が高く、安定的で長期的なシミュレーションを可能にする。
• 新規の根の特定手順を用いたHCMTシミュレーションは、急激な水深の増加上での孤立波の崩壊や、正弦波リップル上でのブラッグ散乱といった複雑な非線形現象を的確に再現した。
• 結果として、本手法により、任意の水深上での完全非線形水波の効率的で次元削減された高精度なシミュレーションが可能となり、直接数値解法と同等の精度を達成するが、著しく低い計算コストで実現できることが示された。