[论文解读] On FFT-based convolutions and correlations, with application to solving Poisson's equation in an open rectangular pipe

Robert D. Ryne|arXiv (Cornell University)|Nov 21, 2011

Particle accelerators and beam dynamics参考文献 6被引用 11

一句话总结

本文提出了一种基于FFT的方法，通过与格林函数的混合卷积和相关运算，求解开放矩形管道中的泊松方程，实现了在束流动力学模拟中高效且精确的势能计算。该方法通过适当的零填充和方向性变换，利用FFT将广泛使用的Hockney方法扩展至有界几何结构，支持标准格林函数与积分格林函数（IGF）技术，从而提升计算效率。

ABSTRACT

A new method is presented for solving Poisson's equation inside an open-ended rectangular pipe. The method uses Fast Fourier Transforms (FFTs) to perform mixed convolutions and correlations of the charge density with the Green function. Descriptions are provided for algorithms based on the ordinary Green function and for an integrated Green function (IGF). Due to its similarity to the widely used Hockney algorithm for solving Poisson's equation in free space, this capability can be easily implemented in many existing particle-in-cell beam dynamics codes.

研究动机与目标

• 开发一种在加速器物理中常见的开放矩形管道几何结构下，求解泊松方程的数值高效方法。
• 将广泛使用的基于FFT的Hockney自由空间泊松求解方法扩展至具有导电壁的有界几何结构。
• 在FFT框架内同时支持普通格林函数与积分格林函数（IGF），以提升精度和计算效率。
• 提供一种计算可扩展且模块化的解决方案，可轻松集成到现有的粒子网格（PIC）束流动力学代码中。

提出的方法

• 该方法使用三维快速傅里叶变换（FFTs）计算电荷密度与格林函数在开放矩形管道中的混合卷积和相关运算。
• 通过对电荷密度进行零填充，并对格林函数进行周期性延拓，利用FFT实现循环卷积，确保在物理区域内的结果正确。
• 通过方向性FFT处理混合边界条件：根据格林函数参数中的符号，在某些维度上执行正向变换，在其他维度上执行反向变换。
• 该方法同时支持标准格林函数与积分格林函数（IGF），后者可提升精度并降低计算成本。
• 该方法利用FFT不仅能计算卷积，还能通过在特定维度上反转傅里叶变换方向来计算相关运算的特性。
• 该实现设计为与现有基于Hockney方法的束流动力学代码兼容，仅需极少修改即可集成。

实验结果

研究问题

• RQ1 基于FFT的卷积与相关运算技术能否被适配以求解具有导电壁的开放矩形管道中的泊松方程？
• RQ2 如何利用格林函数将Hockney方法在效率与模块化方面的优势扩展至有界几何结构？
• RQ3 在管道几何结构中，基于FFT框架使用积分格林函数（IGF）在计算与精度方面有何优势？
• RQ4 如何通过方向性FFT在多维空间中高效实现混合卷积与相关运算？

主要发现

• 基于FFT的方法通过混合卷积与相关运算成功计算了势能，其计算复杂度与Hockney方法相同：O((2N)^d (log₂ 2N)^d)。
• 该方法正确处理了开放矩形管道的格林函数，其依赖于(x±x′, y±y′, z±z′)，通过方向性FFT考虑了参数中符号的变化。
• 通过零填充与周期性延拓，确保了FFT计算在未填充区域中恢复正确的物理卷积或相关运算结果。
• 该算法同时支持标准格林函数与积分格林函数（IGF），IGF在某些束流构型下可提升精度与效率。
• 该方法与现有基于Hockney的代码向后兼容，可实现无缝集成，仅需极少的实现工作量。