[論文レビュー] Compressive spectral method for the simulation of the water waves

Cihan Bayındır|arXiv (Cornell University)|Dec 19, 2015

Ocean Waves and Remote Sensing参考文献 11被引用数 3

ひとこと要約

本稿では、時間または周波数ドメインにおける水波のスパarsityを活用することで、スペクトルシミュレーションにおける計算コストを低減する圧縮スペクトル法を提案する。少ないスペクトル成分を用いた圧縮サンプリングを用いることで、古典的な高分解能スペクトル法とほぼ同一の精度を達成しながら、計算時間を顕著に短縮する。特に、KdV方程式に従うソリタリ波やJonswapスペクトルを有する線形海洋波の長時間シミュレーションにおいて顕著な効果を示す。

ABSTRACT

In this paper an approach for decreasing the computational effort required for the spectral simulations of the water waves is introduced. Signals with majority of the components zero, are known as the sparse signals. Like majority of the signals in the nature it can be realized that water waves are sparse either in time or in the frequency domain. Using the sparsity property of the water waves in the time or in the frequency domain, the compressive sampling algorithm can be used as a tool for improving the performance of the spectral simulation of the water waves. The methodology offered in this paper depends on the idea of using a smaller number of spectral components compared to the classical spectral method with a high number of components. After performing the time integration with a smaller number of spectral components and using the compressive sampling technique, it is shown that the water wave field can be reconstructed with a significantly better efficiency compared to the classical spectral method with a high number of spectral components, especially for long time evolutions. For the sparse water wave model in the time domain the well-known solitary wave solutions of the Korteweg-deVries (KdV) equation is considered. For the sparse water wave model in the frequency domain the well-known Airy (linear) ocean waves with Jonswap spectrum is considered. Utilizing a spectral method, it is shown that by using a smaller number of spectral components compared to the classical spectral method with a high number of components, it is possible to simulate the sparse water waves with negligible error in accuracy and a great efficiency especially for large time evolutions.

研究の動機と目的

• 時間または周波数ドメインにおける水波のスパarsityを活用することで、スペクトルシミュレーションにおける計算コストを低減すること。
• 少ないスペクトル成分を用いた圧縮サンプリングが、多数の成分を要する古典的スペクトル法と比較して高い精度を維持できることを示すこと。
• KdV方程式に従うソリタリ波およびJonswapスペクトルを有する線形海洋波の長時間シミュレーションにおける効率性の向上を図ること。
• 古典的スペクトル法との数値的比較を通じて、本手法の精度および計算コストの削減を検証すること。
• 津波モデリングのような計算負荷の高い応用分野における圧縮スペクトル法の可能性を検討すること。

提案手法

• 本手法は、時間（ソリタリ波）または周波数（線形海洋波）におけるスパarsityを活用して、少ないスペクトル成分から波場を再構成する圧縮サンプリングを用いる。
• ソリタリ波のシミュレーションでは、双曲正割二乗（sech²）解を有するKorteweg-de Vries（KdV）方程式を4次ルンゲ・クッタ時間積分法とFFTに基づくスペクトル微分を用いてシミュレートする。
• 線形海洋波のシミュレーションでは、周波数ドメインにおけるスパarsityをモデル化するためJonswapスペクトルを用い、深さ依存項を有するフーリエ級数で速度ポテンシャルを表現する。
• 圧縮スペクトル法は、アンダーサンプリングされたスペクトルデータから波場を再構成するためにl₁最小化を適用し、必要な成分数を削減する。
• 時間積分は4次ルンゲ・クッタ法を用い、スペクトル微分は高速フーリエ変換（FFT）で計算する。
• 本手法は、N=1024およびM=128/256の成分を用いた古典的スペクトル法と比較して、再構成された波場およびスペクトルを評価する。

実験結果

リサーチクエスチョン

• RQ1 圧縮サンプリングは、水波のスペクトルシミュレーションにおける計算コストを低減しつつ、精度を保持できるか？
• RQ2 長時間シミュレーションにおいて、圧縮スペクトル法の性能は古典的スペクトル法と比較して精度および計算時間の面でどのように異なるか？
• RQ3 時間（ソリタリ波）または周波数（線形海洋波）におけるスパarsityが、少ないスペクトル成分で効率的な再構成を可能にする程度はどの程度か？
• RQ4 スペクトル成分数（M）を変化させた場合、再構成精度および計算効率にどのような影響があるか？
• RQ5 圧縮スペクトル法は非線形波動および津波モデリングのシナリオに効果的に適用可能か？

主な発見

• N=1024、M=128の線形波シミュレーションにおいて、古典的手法の61.90秒から本手法では8.91秒にまで計算時間を短縮し、85.8％の削減を達成した。
• M=256、120,000ステップのシミュレーションでは、古典的手法の256.55秒から本手法では52.04秒にまで短縮され、79.7％の削減を達成した。誤差は無視できるほど小さく（RMS差 = 0.0070）、精度にほとんど影響がなかった。
• M=256の場合、古典的および圧縮スペクトル法の波形プロファイル間の平均二乗誤差はわずか0.0165にとどまり、ほぼ同一の精度を示した。
• 同じシミュレーションにおいて、120,000ステップで古典的手法と比較して79.7％の計算時間短縮が達成された。
• 長時間シミュレーションにおいて、本手法は計算コストの優位性がさらに顕著になり、時間ステップ数が増加するにつれて計算効率が向上した。
• M=128（N=1024の12.5％）のわずかなスペクトル成分数でも、RMS差 < 0.01という最小限の精度損失を維持した。