[論文レビュー] A novel nonlocal potential solver based on nonuniform FFT for efficient simulation of the Davey-Stewartson equations

Norbert J. Mauser, Hans Peter Stimming+1 more|arXiv (Cornell University)|Sep 6, 2014

Numerical methods for differential equations参考文献 24被引用数 3

ひとこと要約

本稿では、デイビー＝スチュアートソン方程式における非局所ポテンシャルのための、非一様フーリエ変換（NUFFT）に基づく新規スムージング法を提案する。標準的なFFTでは、フーリエ乗数の不連続性に起因する精度制限が生じるが、これを克服する。極座標におけるポテンシャル積分の再定式化と、高周波成分と低周波成分への分割（それぞれ標準FFTと高次数ガウス求積により、NUFFTで高速化）により、$O(N\log N)$の計算量で、指数的減衰解では機械精度に達する大幅な精度向上が達成された。代数的減衰解に対しても顕著な改善が得られた。

ABSTRACT

We propose an efficient and accurate solver for the nonlocal potential in the Davey-Stewartson equation using nonuniform FFT (NUFFT). A discontinuity in the Fourier transform of the nonlocal potential causes accuracy locking if the potential is solved by standard FFT with periodic boundary conditions on a truncated domain. Using the fact that the discontinuity disappears in polar coordinates, we reformulate the potential integral and split it into high and low frequency parts. The high frequency part can be approximated by the standard FFT method, while the low frequency part is evaluated with a high order Gauss quadrature accelerated by nonuniform FFT. The NUFFT solver has $O(N\log N)$ complexity, where $N$ is the total number of discretization points, and achieves higher accuracy than standard FFT solver, which makes it a good alternative in simulation. Extensive numerical results show the efficiency and accuracy of the proposed method.

研究の動機と目的

• デイビー＝スチュアートソン方程式における非局所ポテンシャルの標準FFTスムージング法が、$\mathbf{k}=\mathbf{0}$で不連続となるフーリエ乗数 $k_x^2/(k_x^2 + k_y^2)$ のため精度が低下することを解消する。
• 微細な空間グリッドでも、零周波数近傍の不連続性の分解能が不足するため生じる「精度ロック」現象を克服する。
• 非局所ポテンシャルの高速かつ高精度で効率的な数値解法を構築し、$O(N\log N)$の計算量を維持しながら、低周波成分の分解能を向上させる。
• 非局所結合項を高精度に解釈することで、特に指数的減衰を示す解の高精度なシミュレーションを可能にする。
• 非線形偏微分方程式における他の非局所ポテンシャル問題へのNUFFTベースの手法の拡張に基盤を提供する。

提案手法

• フーリエ空間における非局所ポテンシャル積分を極座標に再定式化し、座標変換のヤコビアンのおかげで、乗数 $k_x^2/(k_x^2 + k_y^2)$ の不連続性が消失することを活用する。
• ポテンシャル積分を高周波成分と低周波成分に分割：高周波成分は一様グリッド上の標準FFTで計算する。
• 低周波成分は極座標における高次数ガウス求積により評価し、非一様フーリエ変換（NUFFT）を用いて散在する求積点での高速評価を実現する。
• NUFFTを用いて非一様求積点を一様グリッドに補間し、スペクトル評価を維持する。計算量は $O(N\log N)$ を保つ。
• 高周波成分と低周波成分の寄与を組み合わせ、より高い精度の非局所ポテンシャルを再構築する。
• 時間分割スペクトル法（TSFP）におけるデイビー＝スチュアートソン系全体の計算に、このスムージング法を応用し、自己結合項の計算にNUFFTポテンシャルスムージング法を用いる。

実験結果

リサーチクエスチョン

• RQ1 非局所ポテンシャルのフーリエ乗数における不連続性が、デイビー＝スチュアートソン方程式のスペクトルスムージング法における精度低下を引き起こす問題を、効果的に緩和できるか？
• RQ2 高周波成分には標準FFT、低周波成分にはNUFFTで高速化された高次数求積を組み合わせたハイブリッド手法が、$O(N\log N)$の計算量を維持しながら、優れた精度を達成できるか？
• RQ3 指数的減衰を示す解に対して、新規のNUFFTベースのスムージング法は、標準FFTに比べてどの程度精度が向上するか？
• RQ4 代数的減衰を示す解に対しては、切断誤差が支配的となるが、本手法の性能はいかがなものか？
• RQ5 NUFFTベースのポテンシャルスムージング法は、デイビー＝スチュアートソン系全体の時間分割スペクトル法に効果的に統合可能か？

主な発見

• 指数的減衰を示す波動関数に対して、NUFFTベースのスムージング法は機械精度に達する精度を実現し、標準FFTスムージング法を著しく上回る。
• 代数的減衰を示すラムプ解に対しても、NUFFTスムージング法は標準FFTに比べて誤差を数桁以上小さく抑えているが、依然として切断誤差が支配的である。
• 計算量は標準FFTと同等の $O(N\log N)$ を維持しながら、顕著な精度向上を達成しており、大規模シミュレーションにおいて極めて効率的である。
• 数値結果から、細かいグリッド（$h=1/8$）における波動関数およびポテンシャル誤差はNUFFTで $O(10^{-13})$ に低下するが、標準FFTでは $O(10^{-3})$ から $O(10^{-2})$ に plateau に達する。
• 初期条件や問題設定の違いに関わらず、本手法は一貫した精度向上を示し、有限時間 blowup を示すケースに対してもロバストである。
• 特に高精度が求められる非線形偏微分方程式における非局所ポテンシャルスムージング法として、標準FFTの強力な代替手段となり得る。