[论文解读] 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）的新方法，用于求解Davey-Stewartson方程中的非局部势能，克服了标准FFT因傅里叶乘子在$\mathbf{k}=\mathbf{0}$处不连续而导致的精度限制。通过在极坐标下重写势能积分，并将其分解为高频与低频分量——分别采用标准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求解器在Davey-Stewartson方程非局部势能计算中，由于傅里叶乘子$k_x^2/(k_x^2 + k_y^2)$在$\mathbf{k}=\mathbf{0}$处不连续而导致的精度下降问题。
• 克服因零频附近不连续性解析能力差而引起的“精度锁定”现象，即使在空间网格非常精细时亦然。
• 开发一种快速、精确且高效的非局部势能数值方法，保持$O(N\log N)$复杂度的同时，提升低频分量的分辨率。
• 通过精确解析非局部耦合项，实现对具有指数衰减特性的Davey-Stewartson系统高精度模拟。
• 为将基于NUFFT的方法推广至其他非局部势能问题（如非线性PDE中的类似问题）提供基础。

提出的方法

• 在傅里叶空间中使用极坐标重写非局部势能积分，由于坐标变换的雅可比行列式，傅里叶乘子$k_x^2/(k_x^2 + k_y^2)$的不连续性被消除。
• 将势能积分分解为高频与低频分量：高频部分在均匀网格上通过标准FFT计算。
• 低频部分采用极坐标下的高阶高斯求积法计算，并利用非均匀快速傅里叶变换（NUFFT）加速在散乱求积点上的高效求值。
• 使用NUFFT将非均匀求积点插值到均匀网格上，以实现谱精度计算，同时保持$O(N\log N)$复杂度。
• 将高频与低频分量的贡献合并，重构出完整非局部势能，实现更高精度。
• 将该求解器集成到时间分裂谱方法（TSFP）中，用于求解完整的Davey-Stewartson系统，利用NUFFT势能求解器计算自耦合项。

实验结果

研究问题

• RQ1 在Davey-Stewartson方程的谱求解器中，非局部势能傅里叶乘子的不连续性所导致的精度下降问题是否能被有效缓解？
• RQ2 结合标准FFT处理高频分量与NUFFT加速的高阶求积处理低频分量的混合方法，是否能在保持$O(N\log N)$复杂度的同时实现更优的精度？
• RQ3 与标准FFT相比，该新型NUFFT求解器在指数衰减解上的精度提升程度如何？
• RQ4 对于代数衰减解，该求解器表现如何，尤其在截断误差可能占主导地位的情况下？
• RQ5 基于NUFFT的势能求解器能否有效集成到完整Davey-Stewartson系统的分时间步谱方法中？

主要发现

• 对于具有指数衰减的波函数，NUFFT求解器可达到机器精度，显著优于标准FFT求解器。
• 对于具有代数衰减的孤立子解，即使截断误差仍占主导，NUFFT求解器仍将势能误差降低了数个数量级。
• 该方法保持了$O(N\log N)$复杂度，与标准FFT相当，同时实现了显著更高的精度，使其在大规模模拟中极具效率。
• 数值结果表明，使用NUFFT时，波函数和势能误差在细网格（$h=1/8$）下可降至$O(10^{-13})$，而标准FFT的误差则稳定在$O(10^{-3})$至$O(10^{-2})$之间。
• 该求解器在不同初始条件和问题设置下均表现出鲁棒性，包括有限时间爆破情形，且精度提升一致。
• 所提出的方法是非线性PDE中非局部势能求解器的强有力替代方案，尤其适用于对高精度有严格要求的应用场景。