电磁场学中亥姆赫兹定理证明过程的理解

什么样的矢量偏微分方程可以求解呢？
亥姆霍兹定理就是用来解决这个问题的。该定理表示：如果一个矢量场的散度、旋度和边界条件确定的话，那么就可以通过散度、旋度和边界条件来求解出该矢量场，而且该解唯一。
求解矢量场的思路是：

1. 将矢量场表示为$\delta (r)$函数与矢量的空间卷积。
2. 由于$\delta (r)$函数为标量函数$-\frac{1}{4\pi }\frac{1}{|r|}$的拉普拉斯运算$\delta (r)=-\frac{1}{4\pi }{\nabla }^2\frac{1}{|r|}$，因此矢量场就与散度和旋度产生了关系。
3. 通过矢量三重积公式${\nabla }^{2}A=\nabla \times (\nabla \times A)-\nabla (\nabla \cdot A)$，从而可以将矢量场表示为一无散场$\nabla \times A$的旋度与一标量场$\nabla \cdot A$的梯度之差。

在网上看了百度文库中的《亥姆霍兹定理的证明》【1】、百度文库中《格林公式证明亥姆霍兹定理》【2】、知乎中的《亥姆霍兹定理》【3】。文献【1】中的证明最简洁清晰；【2】中竟然对矢量函数应用格林定理，又对标量函数应用矢量三重积公式，我没明白是什么技巧，而且很明显多处公式错误，真有点怀疑作者和编辑；【3】中也存在对矢量函数应用格林定理的毛病（对公式4中最后一项就是用格林公式展开的）。
下面把自己对【1】的理解重新整理一下，更改了其中部分错误，优化局部的证明思路！

一、 $\delta (r)$函数和空间卷积

我认为这是亥姆霍兹定理的基点，类似信号系统中的冲激响应函数。有了该函数，矢量函数就可以表示为自己与$\delta (r)$的空间卷积了，即$\bar{F}(r)={\iiint }_V\bar{F}(r^{\prime })\delta (r-r^{\prime })d{V}^{\prime }$。
由于$\delta (r)=-\frac{1}{4\pi }{\nabla }^2(\frac{1}{|r|})$，所以，
$$\bar{F}(r)=-\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime }){\nabla }^2(\frac{1}{|r-r\prime |})d{V}^{\prime }$$ (1)
基于该卷积中的拉普拉斯算子，从而引导出来了积分、散度和旋度等对象，从而得到了亥姆霍兹定理：矢量场可以表示为旋度和梯度。

二、矢量三重积

1.拉普拉斯算子移到积分外

$$-\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime }){\nabla }^2(\frac{1}{|r-r\prime |})d{V}^{\prime }=-\frac{1}{4\pi }{\nabla }^2{\iiint }_V\bar{F}(r^{\prime })(\frac{1}{|r-r\prime |})d{V}^{\prime }$$（2）
上式关键难点是：式左边是拉普拉斯运算是对标量，式右边是拉普拉斯运算是对矢量。为什么相等？
首先拉普拉斯算子是对 r求导，而不是对r’求导；其次，上式是一个矢量相等公式，相当于三个标量等式。写详细一些，对上式右边应用拉普拉斯算子的定义公式$({\nabla }^{2}A{)}_x=\frac{{\partial }^2{A}_x}{\partial x^2}+\frac{{\partial }^2{A}_x}{\partial y^2}+\frac{{\partial }^2{A}_x}{\partial z^2}$展开即可得证：
− 1 4 π ∇ 2 ∭ V F ˉ ( r ′ ) ( 1 ∣ r − r ′ ∣ ) d V ′ = − 1 4 π ∭ V { e x ( ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 + ∂ 2 ∂ z 2 ) [ F ˉ x ′ ( r ′ ) ( 1 ∣ r − r ′ ∣ ) ] + e y ( ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 + ∂ 2 ∂ z 2 ) [ F ˉ y ′ ( r ′ ) ( 1 ∣ r − r ′ ∣ ) ] + e z ( ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 + ∂ 2 ∂ z 2 ) [ F ˉ z ′ ( r ′ ) ( 1 ∣ r − r ′ ∣ ) ] } d V ′ = − 1 4 π ∭ V F ˉ ( r ′ ) ∇ 2 ( 1 ∣ r − r ′ ∣ ) d V ′ -\frac{1}{4\pi}\nabla^2 \iiint_V \bar{F}(r^\prime) (\frac{1}{|r-r\prime|})dV^\prime = \\ -\frac{1}{4\pi}\iiint_V \{ e_{x} (\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2})[\bar{F}_{x^\prime}(r^\prime) (\frac{1}{|r-r\prime|})] \\+ e_{y} (\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2})[\bar{F}_{y^\prime}(r^\prime) (\frac{1}{|r-r\prime|})] \\+ e_{z} (\frac{\partial^2}{\partial x^2}+\frac{\partial^2}{\partial y^2}+\frac{\partial^2}{\partial z^2})[\bar{F}_{z^\prime}(r^\prime) (\frac{1}{|r-r\prime|})] \}dV^\prime \\ = -\frac{1}{4\pi} \iiint_V \bar{F}(r^\prime) \nabla^2(\frac{1}{|r-r\prime|})dV^\prime −4π1​∇2∭V​Fˉ(r′)(∣r−r′∣1​)dV′=−4π1​∭V​{ex​(∂x2∂2​+∂y2∂2​+∂z2∂2​)[Fˉx′​(r′)(∣r−r′∣1​)]+ey​(∂x2∂2​+∂y2∂2​+∂z2∂2​)[Fˉy′​(r′)(∣r−r′∣1​)]+ez​(∂x2∂2​+∂y2∂2​+∂z2∂2​)[Fˉz′​(r′)(∣r−r′∣1​)]}dV′=−4π1​∭V​Fˉ(r′)∇2(∣r−r′∣1​)dV′（3）

2. 将拉普拉斯运算转换为梯度偏导和矢量三重积之差

$${\nabla }^{2}A=\nabla (\nabla \cdot A)-\nabla \times (\nabla \times A)$$
上式为矢量三重积公式，即可以得到向量场的求解公式：
$\bar{F}(r)=-\nabla \Phi (r)+\nabla \times \bar{A}(r)$（4）
其中：
$\Phi (r)=\nabla \cdot (\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime })(\frac{1}{|r-r\prime |})d{V}^{\prime })$（5）
$\bar{A}(r)=\nabla \times (\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime })(\frac{1}{|r-r\prime |})d{V}^{\prime })$（6）

3. $\Phi (r)$的进一步推导

应用等式：$\nabla \cdot (\bar{F}(r^{\prime })\delta (r))=\bar{F}(r^{\prime })\cdot \nabla \delta (r)$ 以及等式：$\nabla (\frac{1}{|r-r^{\prime }|})=-{\nabla }^{\prime }(\frac{1}{|r-r^{\prime }|})$ ，可以得到：
$\Phi (r)=\nabla \cdot (\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime })(\frac{1}{|r-r\prime |})d{V}^{\prime })=\frac{1}{4\pi }{\iiint }_V\frac{\nabla \cdot \bar{F}(r^{\prime })}{|r-r^{\prime }|}d{V}^{\prime }-\frac{1}{4\pi }{\iint }_{S^{\prime }}\frac{\bar{F}(r^{\prime })\cdot d{S}^{\prime }}{|r-r^{\prime }|}$
（7）

4. $\bar{A}(r)$ 的进一步推导

这个推导难度较大。
首先要证明：$\nabla \times (\frac{\bar{F}(r^{\prime })}{|r-r^{\prime }|})=-\bar{F}(r^{\prime })\times \nabla (\frac{1}{|r-r^{\prime }|})$ （8）
由于 $\nabla \times (\frac{\bar{F}(r^{\prime })}{|r-r^{\prime }|})=\frac{\nabla \times \bar{F}(r^{\prime })}{|r-r^{\prime }|}-\bar{F}(r^{\prime })\times \nabla (\frac{1}{|r-r^{\prime }|})$ ，因为 $\nabla \times \bar{F}(r^{\prime })=0$，故上式成立！
另外一种证明方式是将叉乘写成行列式形式，求偏导时 可以视为常量，因此下面推导成立[3]：
$$\left|\begin{array}{ccc}{e}_x& e_y& e_z\\ \frac{\partial }{\partial x}& \frac{\partial }{\partial y}& \frac{\partial }{\partial z}\\ \frac{{\bar{F}}_{x^{\prime }}(r^{\prime })}{|r-r^{\prime }|}& \frac{{\bar{F}}_{y^{\prime }}(r^{\prime })}{|r-r^{\prime }|}& \frac{{\bar{F}}_{z^{\prime }}(r^{\prime })}{|r-r^{\prime }|}\end{array}\right|=\left|\begin{array}{ccc}{e}_x& e_y& e_z\\ \frac{\partial }{\partial x}\frac{1}{|r-r^{\prime }|}& \frac{\partial }{\partial y}\frac{1}{|r-r^{\prime }|}& \frac{\partial }{\partial z}\frac{1}{|r-r^{\prime }|}\\ {\bar{F}}_{x^{\prime }}& {\bar{F}}_{y^{\prime }}(r^{\prime })& {\bar{F}}_{z^{\prime }}(r^{\prime })\end{array}\right|=-\left|\begin{array}{ccc}{e}_x& e_y& e_z\\ {\bar{F}}_{x^{\prime }}& {\bar{F}}_{y^{\prime }}(r^{\prime })& {\bar{F}}_{z^{\prime }}(r^{\prime })\\ \frac{\partial }{\partial x}\frac{1}{|r-r^{\prime }|}& \frac{\partial }{\partial y}\frac{1}{|r-r^{\prime }|}& \frac{\partial }{\partial z}\frac{1}{|r-r^{\prime }|}\end{array}\right|$$
其次要证明：
${\iiint }_V\nabla \times \bar{F}(r^{\prime })d{V}^{\prime }={\iint }_S\overline{d{S}^{\prime }}\times \bar{F}(r^{\prime })=-{\iint }_S\bar{F}(r^{\prime })\times \overline{d{S}^{\prime }}$（9）
将等式左边的叉乘进行行列式展开，并且由于${\iiint }_V\frac{\partial F}{\partial x}dxdydz={\iint }_{S}Fdydz$，因此可推导得出下式成立：
$$\begin{array}{rl}\nabla \times \bar{F}(r^{\prime })d{V}^{\prime }& =\left|\begin{array}{ccc}{e}_x& e_y& e_z\\ \frac{\partial }{\partial x^{\prime }}& \frac{\partial }{\partial y^{\prime }}& \frac{\partial }{\partial z^{\prime }}\\ {\bar{F}}_x(r^{\prime })d{V}^{\prime }& {\bar{F}}_y(r^{\prime })d{V}^{\prime }& {\bar{F}}_z(r^{\prime })d{V}^{\prime }\end{array}\right|\\ & =\left|\begin{array}{ccc}{e}_x& e_y& e_z\\ \frac{\partial }{\partial x^{\prime }}d{V}^{\prime }& \frac{\partial }{\partial y^{\prime }}d{V}^{\prime }& \frac{\partial }{\partial z^{\prime }}d{V}^{\prime }\\ {\bar{F}}_x(r^{\prime })& {\bar{F}}_y(r^{\prime })& {\bar{F}}_z(r^{\prime })\end{array}\right|\\ & =\nabla d{V}^{\prime }\times \bar{F}(r^{\prime })\\ & =\overline{d{S}^{\prime }}\times \bar{F}(r^{\prime })\\ & =-\bar{F}(r^{\prime })\times \overline{d{S}^{\prime }}\end{array}$$（10）
$\bar{A}(r)$进一步展开的推导如下：
$$\begin{array}{rl}\bar{A}(r)& =\nabla \times (\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime })(\frac{1}{|r-r\prime |})d{V}^{\prime })\\ & =\frac{1}{4\pi }{\iiint }_V\nabla \times (\frac{\bar{F}(r^{\prime })}{|r-r^{\prime }|})d{V}^{\prime }\\ & =-\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime })\times \nabla (\frac{1}{|r-r^{\prime }|})d{V}^{\prime }\\ & =\frac{1}{4\pi }{\iiint }_V\bar{F}(r^{\prime })\times {\nabla }^{\prime }(\frac{1}{|r-r^{\prime }|})d{V}^{\prime }\\ & =\frac{1}{4\pi }{\iiint }_V(\frac{{\nabla }^{\prime }\times \bar{F}(r^{\prime })}{|r-r^{\prime }|}-{\nabla }^{\prime }\times (\frac{\bar{F}(r^{\prime })}{|r-r^{\prime }|}))d{V}^{\prime }\\ & =\frac{1}{4\pi }{\iiint }_V\frac{{\nabla }^{\prime }\times \bar{F}(r^{\prime })}{|r-r^{\prime }|}d{V}^{\prime }-\frac{1}{4\pi }{\iint }_S\frac{\overline{d{S}^{\prime }}\times \bar{F}(r^{\prime })}{|r-r^{\prime }|}\end{array}$$（11）
三、矢量场唯一性定理证明
由上可知，矢量场可以通过其散度、旋度和边界S上约束来计算。进一步可以证明矢量的散度、旋度和边界条件确定的矢量场是唯一的。使用反证法，假设不唯一，则存在$\delta \bar{F}={\bar{F}}_1-{\bar{F}}_2$，则：
$$\{\begin{array}{rl}& \nabla \cdot \delta \bar{F}=0\\ & \nabla \times \delta \bar{F}=0\end{array}$$
并且在边界曲面S上法向分量或切向分量给定，则：
$$\begin{array}{rl}& e_n\cdot \delta \bar{F}=0，或\\ & e_n\times \delta \bar{F}=0\end{array}$$，
假设 $\delta \bar{F}=\nabla \Phi$，即${\nabla }^2\Phi =0$。应用标量函数第一格林公式可得：
${\iiint }_V(|\nabla \Phi {|}^2+\Phi {\nabla }^2\Phi )d{V}^{\prime }={\iint }_S\Phi \frac{\partial \Phi }{\partial n}d{S}^{\prime }$（12）
${\iiint }_V(|\nabla \Phi {|}^2)d{V}^{\prime }={\iint }_S\Phi \frac{\partial \Phi }{\partial n}d{S}^{\prime }$（13）
如果$e_n\cdot \delta \bar{F}=0$，即式（13）左边项为0，从而 $\delta \bar{F}=0$。
如果在边界面S上$e_n\times \delta \bar{F}=0$，则在边界S上$\Phi$ 为常数。结合散度定理：${\iint }_V{\nabla }^2\Phi d{V}^{\prime }={\iint }_S\nabla \Phi \cdot \overline{d{S}^{\prime }}={\iint }_S\frac{\partial \Phi }{\partial n}d{S}^{\prime }$ ，可做如下推导：
$$\begin{array}{rl}{\iint }_S\Phi \frac{\partial \Phi }{\partial n}d{S}^{\prime }& ={\Phi }_s{\iint }_S\frac{\partial \Phi }{\partial n}d{S}^{\prime }\\ & ={\Phi }_s{\iiint }_V{\nabla }^2\Phi d{V}^{\prime }\\ & =0\end{array}$$
补充证明：如果在边界面S上$e_n\times \nabla \Phi =0$，则在边界S上$\Phi$为常数。
设边界曲面S的函数为：$\Psi (x,y,z)=0$ 。
因曲面S上： $e_n\times \nabla \Phi =0$，即$\frac{\nabla \Psi }{|\nabla \Psi |}\times \nabla \Phi =0$
所以， $\nabla \Psi \times \nabla \Phi =0$
所以， $\nabla \Psi =k\nabla \Phi$
所以曲面S上，$\Phi (x,y,z)=k\Psi (x,y,z)+Constant=Constant$ ，得证！