【GAMES201学习笔记】快速多极展开

1. 相关概念

1.1 快速多极算法简介

快速多极展开算法（FMM）最初是用来处理大量粒子之间相互影响的问题，也可以用来计算基于 位势理论 的拉普拉斯方程的数值解，同时也可用来处理泊松方程。

同理，我们也可以用相关场的求解过程来辅助理解 FMM 求解泊松方程的这一过程，该处以多体系统的静电场为例。

1.2 静电场

1.2.1 公式的离散形式

对于空间 N 个带电粒子构成的系统，第 i 个粒子所在位置的静电势 $\Phi (x_i)$ ：

$$\Phi (x_i)=\sum _{j=1,j\ne i}^N\frac{m_j}{r_{ij}}$$

所受到的静电力 $E(x_i)$ ：

$$E(x_i)=\sum _{j=1,j\ne i}^N{m}_j\frac{{\mathbf{x}}_i-{\mathbf{x}}_j}{r_{ij}^3}=\sum _{j=1,j\ne i}^N{m}_j\frac{{\mathbf{x}}_i-{\mathbf{x}}_j}{\Vert {\mathbf{x}}_i-{\mathbf{x}}_j{\Vert }^3}$$

Note ：

• ${\mathbf{x}}_i$ ：表示第 i 个粒子所在位置；
• $m_j$ ：表示第 j 个粒子所占权重，与带电量 $q_j$ 成正比；
• $r_{ij}$ ：表示第 i 个粒子和第 j 个粒子之间的距离；
• $\Vert \cdot \Vert$ ：表示两点间的距离；
• 求解 $E(x_i)$ 相当于对 i 粒子与 i 粒子以外的所有粒子的相互力做一个累和，复杂度为 $O(N)$ ，同理，如果对全部 N 个粒子都使用该方式进行求解，复杂度高达 $O(N^2)$ 。

1.2.2 公式的积分形式

对于空间连续介质问题，假设粒子是连续分布在空间曲面 $S$ 上，则离散的多粒子系统可以过渡到连续介质系统。在点 $\mathbf{x}$ 处受到整个曲面的势为：

$$\Phi (\mathbf{x})={\int }_S\frac{m(\mathbf{y})}{|\mathbf{x}-\mathbf{y}|}\mathrm{d}S(\mathbf{y}),\mathbf{y}\in S$$

点 $\mathbf{x}$ 处受到整个曲面的力为：

$$E(\mathbf{x})={\int }_{S}m(\mathbf{y})\frac{\mathbf{x}-\mathbf{y}}{|\mathbf{x}-\mathbf{y}{|}^3}\mathrm{d}S(\mathbf{y}),\mathbf{y}\in S$$

Note ：

• $|\cdot |$ ：表示两点间的距离。

1.2.3 公式的复势形式

课程中为了式子尽可能的直观，主要是以该表达方式为例。

对于单个粒子产生的场

位于复平面上的点 ${\mathbf{z}}_j$ 的平面点电荷 $q_j$ 在点 $\mathbf{z}$ 的复势：

$$f_j(\mathbf{z})=\frac{-q_j}{2\pi {ϵ}_0}\ln (\mathbf{z}-{\mathbf{z}}_j)$$

其中的实部为势函数：

$${\varphi }_j(\mathbf{z})=\mathrm{R}\mathrm{e}{f}_j(\mathbf{z})=\frac{q_j}{2\pi {ϵ}_0}\mathrm{R}\mathrm{e}(\ln (\mathbf{z}-{\mathbf{z}}_j))$$

Note ：

• 虚部为通量函数，由于该处我们仅讨论势函数，故忽略。

对于 j 个粒子组成的系统产生的场

点 $\mathbf{c}$ 附近的 $m$ 个电荷对点 $\mathbf{z}$ 的势（直接累和）：

$$\varphi (\mathbf{z})=\frac{1}{2\pi {ϵ}_0}\mathrm{R}\mathrm{e}\sum _{j=1}^m{q}_j(\ln (\mathbf{z}-{\mathbf{z}}_j))$$

该处简化为：

$$\varphi (\mathbf{z})=\sum _{j=1}^m{q}_j\ln (\mathbf{z}-{\mathbf{z}}_j)\text{\hspace{1em}(1)}$$

将 $\ln (\mathbf{z}-{\mathbf{z}}_j)$ 在点 $\mathbf{z}$ 处做泰勒级数展开，即：

$$\ln (\mathbf{z}-{\mathbf{z}}_j)=\ln \mathbf{z}-\sum _{k=1}^{\infty }\frac{{\mathbf{z}}_j^k}{k{\mathbf{z}}^k}\text{\hspace{1em}(2)}$$

将公式（2）代入公式（1），得到：

ϕ ( z ) = ∑ j = 1 m q j ln ⁡ z − ∑ j = 1 m ∑ k = 1 ∞ q j