亥姆霍兹分解定理(Helmholtz Decomposition Theorem),也称为亥姆霍兹定理(Helmholtz’s Theorem)或基本矢量分解定理,是矢量分析中的一个重要结论。它指出,在适当条件下,任何一个光滑的、快速衰减的三维矢量场都可以分解为**无旋场(保守场)和无散场(螺线管场)**两部分。
定理表述
设 \(\mathbf{F}\) 是一个定义在 \(\mathbb{R}^3\) 上的连续可微矢量场,且当 \(| \mathbf{r} | \to \infty\) 时,\(\mathbf{F}\) 衰减得足够快(例如,比 \(1/r\) 更快),那么 \(\mathbf{F}\) 可以表示为: \[ \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A} \] 其中:
- \(\phi\) 是一个标量场(标量势函数),满足 \(\nabla \phi\) 是无旋场(即 \(\nabla \times (\nabla \phi) = 0\))。
- \(\mathbf{A}\) 是一个矢量场(矢量势函数),满足 \(\nabla \times \mathbf{A}\) 是无散场(即 \(\nabla \cdot (\nabla \times \mathbf{A}) = 0\))。
物理意义
亥姆霍兹分解定理表明,任何矢量场都可以分解为:
- 无旋部分(纵场):由标量势 \(\phi\) 的梯度 \(\nabla \phi\) 描述,代表保守场(如静电场)。
- 无散部分(横场):由矢量势 \(\mathbf{A}\) 的旋度 \(\nabla \times \mathbf{A}\) 描述,代表螺线管场(如静磁场)。
数学推导
在数学上,亥姆霍兹分解可以通过傅里叶变换或格林函数法推导。其核心思想是求解泊松方程: \[ \nabla^2 \phi = -\nabla \cdot \mathbf{F}, \quad \nabla^2 \mathbf{A} = -\nabla \times \mathbf{F} \] 然后利用标量势和矢量势构造分解。
应用
亥姆霍兹分解在物理学中有广泛应用,例如:
- 电磁学:静电场和静磁场的分离(麦克斯韦方程组)。
- 流体力学:不可压缩流体的速度场分解(涡旋和无旋流动)。
- 量子力学:电磁势的规范选择(如库仑规范 \(\nabla \cdot \mathbf{A} = 0\))。
备注
- 在有限区域或有边界条件时,分解可能需要调整(如考虑调和场部分)。
- 该定理是更一般的**霍奇分解(Hodge Decomposition)**在三维欧氏空间的特例。
总之,亥姆霍兹分解定理提供了研究矢量场结构的强大工具,广泛应用于数学物理和工程领域。
1. 傅里叶变换法证明
假设矢量场 \(\mathbf{F}(\mathbf{r})\) 在 \(\mathbb{R}^3\) 上定义,且足够光滑(至少 \(C^2\)),并在无穷远处衰减足够快(如 \(|\mathbf{F}| \to 0\) 比 \(1/r\) 快)。我们可以利用傅里叶变换进行分解。
(1) 傅里叶变换表示
对 \(\mathbf{F}(\mathbf{r})\) 进行傅里叶变换: \[ \hat{\mathbf{F}}(\mathbf{k}) = \int_{\mathbb{R}^3} \mathbf{F}(\mathbf{r}) e^{-i \mathbf{k} \cdot \mathbf{r}} , d^3 \mathbf{r} \] 其逆变换为: \[ \mathbf{F}(\mathbf{r}) = \frac{1}{(2\pi)^3} \int_{\mathbb{R}^3} \hat{\mathbf{F}}(\mathbf{k}) e^{i \mathbf{k} \cdot \mathbf{r}} , d^3 \mathbf{k} \]
(2) 分解 \(\hat{\mathbf{F}}(\mathbf{k})\)
任意矢量场在傅里叶空间可以分解为: \[ \hat{\mathbf{F}}(\mathbf{k}) = \hat{\mathbf{F}}_L(\mathbf{k}) + \hat{\mathbf{F}}_T(\mathbf{k}) \] 其中:
- 纵向部分(无旋):\(\hat{\mathbf{F}}_L(\mathbf{k}) = \frac{\mathbf{k} (\mathbf{k} \cdot \hat{\mathbf{F}})}{|\mathbf{k}|^2}\)(平行于 \(\mathbf{k}\))
- 横向部分(无散):\(\hat{\mathbf{F}}_T(\mathbf{k}) = \hat{\mathbf{F}} - \hat{\mathbf{F}}_L = \hat{\mathbf{F}} - \frac{\mathbf{k} (\mathbf{k} \cdot \hat{\mathbf{F}})}{|\mathbf{k}|^2}\)(垂直于 \(\mathbf{k}\))
(3) 逆变换回实空间
- 纵向部分对应标量势 \(\phi\): \[ \mathbf{F}_L(\mathbf{r}) = -\nabla \phi, \quad \text{其中} \quad \phi(\mathbf{r}) = \frac{1}{(2\pi)^3} \int \frac{i \mathbf{k} \cdot \hat{\mathbf{F}}}{|\mathbf{k}|^2} e^{i \mathbf{k} \cdot \mathbf{r}} , d^3 \mathbf{k} \]
- 横向部分对应矢量势 \(\mathbf{A}\): \[ \mathbf{F}_T(\mathbf{r}) = \nabla \times \mathbf{A}, \quad \text{其中} \quad \mathbf{A}(\mathbf{r}) = \frac{1}{(2\pi)^3} \int \frac{i \mathbf{k} \times \hat{\mathbf{F}}}{|\mathbf{k}|^2} e^{i \mathbf{k} \cdot \mathbf{r}} , d^3 \mathbf{k} \]
因此,得到分解: \[ \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A} \]
2. 格林函数法(泊松方程法)
另一种方法是利用格林函数求解泊松方程,直接构造 \(\phi\) 和 \(\mathbf{A}\)。
(1) 标量势 \(\phi\) 的方程
我们希望: \[ \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A} \] 取散度: \[ \nabla \cdot \mathbf{F} = -\nabla^2 \phi \] 这是一个泊松方程,其解为: \[ \phi(\mathbf{r}) = \frac{1}{4\pi} \int_{\mathbb{R}^3} \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} , d^3 \mathbf{r}’ \]
(2) 矢量势 \(\mathbf{A}\) 的方程
取旋度: \[ \nabla \times \mathbf{F} = \nabla \times (\nabla \times \mathbf{A}) = \nabla (\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A} \] 选择库仑规范(\(\nabla \cdot \mathbf{A} = 0\)),则: \[ \nabla \times \mathbf{F} = -\nabla^2 \mathbf{A} \] 这也是泊松方程,其解为: \[ \mathbf{A}(\mathbf{r}) = \frac{1}{4\pi} \int_{\mathbb{R}^3} \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} , d^3 \mathbf{r}’ \]
(3) 最终分解
因此,我们得到: \[ \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A} \] 其中: \[ \phi(\mathbf{r}) = \frac{1}{4\pi} \int \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} , d^3 \mathbf{r}’, \quad \mathbf{A}(\mathbf{r}) = \frac{1}{4\pi} \int \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} , d^3 \mathbf{r}’ \]
总结
- 傅里叶变换法:在频域分解 \(\hat{\mathbf{F}}\) 为纵向和横向部分,再变换回实空间。
- 格林函数法:直接求解泊松方程得到 \(\phi\) 和 \(\mathbf{A}\)。
两种方法都证明了: \[ \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A} \] 即任何光滑、快速衰减的矢量场都可以分解为无旋部分(标量势梯度)和无散部分(矢量势旋度)。
================== Qwen3 =========================
一、定理内容(数学表达)
设 $\mathbf{F}(\mathbf{r})$ 是定义在三维空间中有足够光滑且衰减较快的向量场(例如在无穷远处趋于零),则可以唯一地将该向量场写成如下形式:
$$ \mathbf{F}(\mathbf{r}) = -\nabla \phi(\mathbf{r}) + \nabla \times \mathbf{A}(\mathbf{r}) $$
其中:
- $\phi(\mathbf{r})$ 是一个标量函数(scalar potential),
- $\mathbf{A}(\mathbf{r})$ 是一个矢量函数(vector potential)。
这个分解意味着:
-
$-\nabla \phi$ 是一个无旋场(irrotational),因为它的旋度为零: $$ \nabla \times (-\nabla \phi) = 0 $$
-
$\nabla \times \mathbf{A}$ 是一个无散场(solenoidal),因为它的散度为零: $$ \nabla \cdot (\nabla \times \mathbf{A}) = 0 $$
二、物理意义
亥姆霍兹分解揭示了向量场的两个基本性质:
- 散度源(由标量势描述):决定了向量场中“发散”的部分,比如电场中的电荷源。
- 旋度源(由矢量势描述):决定了向量场中“旋转”的部分,比如磁场或涡旋流体。
因此,亥姆霍兹分解广泛应用于:
- 流体力学(速度场分解为压缩分量与旋转分量)
- 电磁学(电场和磁场的基本结构)
- 等离子体物理
- 地球物理学等领域
三、亥姆霍兹分解的存在性条件
为了使亥姆霍兹分解成立,通常需要满足以下条件:
- 向量场 $\mathbf{F}(\mathbf{r})$ 在整个空间中是连续可微的;
- 向量场在无穷远处足够快地趋于零,使得积分收敛;
- 分解在一定的边界条件下是唯一的。
四、积分表示式(格林函数法)
在无限空间中,$\phi$ 和 $\mathbf{A}$ 可以通过如下积分公式显式写出:
-
标量势: $$ \phi(\mathbf{r}) = \frac{1}{4\pi} \int_{V} \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
-
矢量势: $$ \mathbf{A}(\mathbf{r}) = \frac{1}{4\pi} \int_{V} \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
其中 $\nabla’$ 表示对源点 $\mathbf{r}’$ 的梯度。
五、推广:斯托克斯分解(Stokes-Helmholtz 分解)
在有界区域上,若考虑边界条件,亥姆霍兹分解会涉及边界项,此时称为斯托克斯分解。这种情况下分解不是唯一的,除非附加边界条件,如 Dirichlet 或 Neumann 条件。
六、例子说明
示例 1:静电场
静电场 $\mathbf{E}$ 满足: $$ \nabla \times \mathbf{E} = 0 \Rightarrow \mathbf{E} = -\nabla \phi $$ 这说明静电场是一个纯无旋场(只有标量势)。
示例 2:静磁场
静磁场 $\mathbf{B}$ 满足: $$ \nabla \cdot \mathbf{B} = 0 \Rightarrow \mathbf{B} = \nabla \times \mathbf{A} $$ 这说明磁场是一个纯无散场(只有矢量势)。
七、总结
类型 | 数学形式 | 物理特性 |
---|---|---|
无旋场(标量势) | $-\nabla \phi$ | $\nabla \times \mathbf{F} = 0$ |
无散场(矢量势) | $\nabla \times \mathbf{A}$ | $\nabla \cdot \mathbf{F} = 0$ |
如果你有具体的应用场景(如电磁学、流体力学等),我也可以为你提供更详细的解释或推导!
任意定义在三维空间中的光滑向量场 $\mathbf{F}(\mathbf{r})$,若其散度和旋度在全空间有定义且衰减较快(例如在无穷远处趋于零),则可以唯一地写成一个无旋场与一个无散场之和: $$ \mathbf{F}(\mathbf{r}) = -\nabla \phi(\mathbf{r}) + \nabla \times \mathbf{A}(\mathbf{r}) $$
一、基本思路
我们希望将给定向量场 $\mathbf{F}$ 分解为两个部分:
- 标量势部分:$\mathbf{F}_\text{irrot} = -\nabla \phi$,满足 $\nabla \times \mathbf{F}_\text{irrot} = 0$
- 矢量势部分:$\mathbf{F}_\text{solen} = \nabla \times \mathbf{A}$,满足 $\nabla \cdot \mathbf{F}_\text{solen} = 0$
从而使得: $$ \mathbf{F} = \mathbf{F}_\text{irrot} + \mathbf{F}_\text{solen} $$
为此,我们需要构造合适的 $\phi$ 和 $\mathbf{A}$,使得这个等式成立。
二、使用格林函数法构造势函数
我们考虑无限空间中的一般情况,并利用格林函数法来构造 $\phi$ 和 $\mathbf{A}$。
格林函数
三维空间中拉普拉斯算子 $\nabla^2$ 的格林函数是: $$ G(\mathbf{r}, \mathbf{r}’) = \frac{1}{4\pi |\mathbf{r} - \mathbf{r}’|} $$
我们可以用它来求解泊松方程: $$ \nabla^2 \psi(\mathbf{r}) = f(\mathbf{r}) \Rightarrow \psi(\mathbf{r}) = \int G(\mathbf{r}, \mathbf{r}’) f(\mathbf{r}’) dV’ $$
三、构造标量势 $\phi$ 和矢量势 $\mathbf{A}$
设 $\mathbf{F}(\mathbf{r})$ 是给定的向量场。
定义:
- 散度源:$\rho(\mathbf{r}) = \nabla \cdot \mathbf{F}(\mathbf{r})$
- 旋度源:$\mathbf{J}(\mathbf{r}) = \nabla \times \mathbf{F}(\mathbf{r})$
我们尝试构造:
$$ \phi(\mathbf{r}) = \frac{1}{4\pi} \int_V \frac{\rho(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ = \frac{1}{4\pi} \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
$$ \mathbf{A}(\mathbf{r}) = \frac{1}{4\pi} \int_V \frac{\mathbf{J}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ = \frac{1}{4\pi} \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
四、验证分解成立
现在我们令:
$$ \mathbf{F}_\text{decomp}(\mathbf{r}) = -\nabla \phi(\mathbf{r}) + \nabla \times \mathbf{A}(\mathbf{r}) $$
我们想验证: $$ \mathbf{F}_\text{decomp}(\mathbf{r}) = \mathbf{F}(\mathbf{r}) $$
即:
$$ \boxed{ \mathbf{F}(\mathbf{r}) = -\nabla \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) + \nabla \times \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) } $$
为了验证这一点,我们可以从以下恒等式出发:
向量恒等式
对于任意光滑向量场 $\mathbf{F}(\mathbf{r})$,在无限空间中有如下恒等式(可由矢量分析推出):
$$ \mathbf{F}(\mathbf{r}) = -\nabla \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) + \nabla \times \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) $$
这正是我们要证明的结论!
五、唯一性说明
亥姆霍兹分解的唯一性依赖于边界条件。
如果我们要求:
- $\phi(\mathbf{r})$ 在无穷远处趋于零;
- $\mathbf{A}(\mathbf{r})$ 在无穷远处也趋于零;
那么这种分解是唯一的。否则,可能存在其他形式的势函数也能表示相同的向量场。
六、总结步骤
- 给定一个向量场 $\mathbf{F}(\mathbf{r})$。
- 计算它的散度 $\rho = \nabla \cdot \mathbf{F}$ 和旋度 $\mathbf{J} = \nabla \times \mathbf{F}$。
- 利用格林函数构造标量势 $\phi$ 和矢量势 $\mathbf{A}$: $$ \phi(\mathbf{r}) = \frac{1}{4\pi} \int_V \frac{\rho(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’, \quad \mathbf{A}(\mathbf{r}) = \frac{1}{4\pi} \int_V \frac{\mathbf{J}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
- 验证: $$ \mathbf{F}(\mathbf{r}) = -\nabla \phi(\mathbf{r}) + \nabla \times \mathbf{A}(\mathbf{r}) $$
- 若边界条件合适(如无穷远处趋于零),则此分解唯一。
七、补充说明
- 如果区域是有界的,还需要考虑边界上的通量或切向分量,此时称为斯托克斯分解。
- 在电磁学中,电场 $\mathbf{E}$ 可以分解为静电势产生的梯度项和磁感应变化引起的旋度项(见麦克斯韦方程组)。
- 在流体力学中,速度场也可以分解为压缩部分(标量势)和旋转部分(矢量势)。
如果你需要具体例子(如对某个具体向量场进行分解)或应用于特定物理问题(如电磁学、流体力学)的推导,我也可以继续为你详细展开!
对于任意光滑向量场 $\mathbf{F}(\mathbf{r})$,在三维无限空间中有: $$ \boxed{ \mathbf{F}(\mathbf{r}) = -\nabla \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) + \nabla \times \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) } $$
这个恒等式是亥姆霍兹分解的数学核心。我们将从矢量分析和格林函数的角度来推导它。
一、基本思想
我们希望将一个给定向量场 $\mathbf{F}(\mathbf{r})$ 表示为两个部分之和:
- 无旋部分(梯度项):由散度 $\rho = \nabla \cdot \mathbf{F}$ 决定;
- 无散部分(旋度项):由旋度 $\mathbf{J} = \nabla \times \mathbf{F}$ 决定。
这类似于求解泊松方程的方法,利用格林函数法来构造标量势和矢量势。
二、使用格林函数表示任意函数
在三维空间中,拉普拉斯算子 $\nabla^2$ 的格林函数为: $$ G(\mathbf{r}, \mathbf{r}’) = \frac{1}{4\pi |\mathbf{r} - \mathbf{r}’|} $$
我们可以用它来求解如下泊松方程: $$ \nabla^2 \psi(\mathbf{r}) = f(\mathbf{r}) \Rightarrow \psi(\mathbf{r}) = \int G(\mathbf{r}, \mathbf{r}’) f(\mathbf{r}’) dV’ $$
三、构造标量势与矢量势
定义:
-
标量势 $\phi(\mathbf{r})$ 满足: $$ \nabla^2 \phi(\mathbf{r}) = -\nabla \cdot \mathbf{F}(\mathbf{r}) \Rightarrow \phi(\mathbf{r}) = \frac{1}{4\pi} \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
-
矢量势 $\mathbf{A}(\mathbf{r})$ 满足: $$ \nabla^2 \mathbf{A}(\mathbf{r}) = -\nabla \times \mathbf{F}(\mathbf{r}) \Rightarrow \mathbf{A}(\mathbf{r}) = \frac{1}{4\pi} \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ $$
注意这里的积分变量是 $\mathbf{r}’$,而 $\nabla’$ 是对 $\mathbf{r}’$ 的梯度。
四、写出总场 $\mathbf{F}$ 的表达式
根据亥姆霍兹分解,我们定义: $$ \mathbf{F}_\text{decomp}(\mathbf{r}) = -\nabla \phi(\mathbf{r}) + \nabla \times \mathbf{A}(\mathbf{r}) $$
代入上面构造的 $\phi$ 和 $\mathbf{A}$,得到:
$$ \mathbf{F}_\text{decomp}(\mathbf{r}) = -\nabla \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) + \nabla \times \left( \frac{1}{4\pi} \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{|\mathbf{r} - \mathbf{r}’|} dV’ \right) $$
这就是我们要推导的目标恒等式!
五、关键恒等式的矢量分析推导
为了更严谨地验证这个恒等式成立,我们可以从以下恒等式出发:
向量恒等式(Green’s second identity for vector fields)
对于任意光滑向量场 $\mathbf{F}(\mathbf{r})$,在无限空间中满足一定衰减条件时,有:
$$ \mathbf{F}(\mathbf{r}) = \int_V G(\mathbf{r}, \mathbf{r}’) \nabla’^2 \mathbf{F}(\mathbf{r}’) dV’ $$
但更直接的方式是使用下面这个矢量分析恒等式:
核心恒等式(关键一步)
$$ \boxed{ \mathbf{F}(\mathbf{r}) = -\nabla \left( \int_V G(\mathbf{r}, \mathbf{r}’) (\nabla’ \cdot \mathbf{F})(\mathbf{r}’) dV’ \right) + \nabla \times \left( \int_V G(\mathbf{r}, \mathbf{r}’) (\nabla’ \times \mathbf{F})(\mathbf{r}’) dV’ \right) } $$
其中 $G(\mathbf{r}, \mathbf{r}’) = \frac{1}{4\pi |\mathbf{r} - \mathbf{r}’|}$
六、推导恒等式的过程(详细步骤)
我们考虑:
$$ \mathbf{F}(\mathbf{r}) = \int_V \delta(\mathbf{r} - \mathbf{r}’) \mathbf{F}(\mathbf{r}’) dV’ $$
由于 $\delta(\mathbf{r} - \mathbf{r}’) = -\nabla^2 G(\mathbf{r}, \mathbf{r}’)$,所以:
$$ \mathbf{F}(\mathbf{r}) = \int_V (-\nabla^2 G(\mathbf{r}, \mathbf{r}’)) \mathbf{F}(\mathbf{r}’) dV’ = -\nabla^2 \int_V G(\mathbf{r}, \mathbf{r}’) \mathbf{F}(\mathbf{r}’) dV’ $$
现在我们应用矢量拉普拉斯恒等式: $$ \nabla^2 \mathbf{F} = \nabla (\nabla \cdot \mathbf{F}) - \nabla \times (\nabla \times \mathbf{F}) $$
于是: $$ \mathbf{F}(\mathbf{r}) = -\nabla^2 \int_V G(\mathbf{r}, \mathbf{r}’) \mathbf{F}(\mathbf{r}’) dV’ = -\int_V G(\mathbf{r}, \mathbf{r}’) \left[ \nabla’^2 \mathbf{F}(\mathbf{r}’) \right] dV’ $$
再展开 $\nabla’^2 \mathbf{F} = \nabla’(\nabla’ \cdot \mathbf{F}) - \nabla’ \times (\nabla’ \times \mathbf{F})$,得:
$$ \mathbf{F}(\mathbf{r}) = -\int_V G(\mathbf{r}, \mathbf{r}’) \left[ \nabla’(\nabla’ \cdot \mathbf{F}) - \nabla’ \times (\nabla’ \times \mathbf{F}) \right] dV’ $$
整理后得到:
$$ \mathbf{F}(\mathbf{r}) = -\int_V G(\mathbf{r}, \mathbf{r}’) \nabla’(\nabla’ \cdot \mathbf{F}) dV’ + \int_V G(\mathbf{r}, \mathbf{r}’) \nabla’ \times (\nabla’ \times \mathbf{F}) dV’ $$
接下来,我们对第一个积分进行分部积分(矢量形式的高斯定理):
分部积分公式(矢量)
$$ \int_V G \nabla’ f dV’ = -\int_V f \nabla’ G dV’ + \text{边界项} $$
假设边界项趋于零(无穷远),则:
$$ \int_V G \nabla’(\nabla’ \cdot \mathbf{F}) dV’ = -\int_V (\nabla’ \cdot \mathbf{F}) \nabla’ G dV’ $$
然后我们把 $\nabla$ 提出来:
$$ -\int_V (\nabla’ \cdot \mathbf{F}) \nabla’ G dV’ = -\nabla \int_V (\nabla’ \cdot \mathbf{F}) G dV’ $$
类似地,对第二个积分:
$$ \int_V G \nabla’ \times (\nabla’ \times \mathbf{F}) dV’ = \nabla \times \int_V (\nabla’ \times \mathbf{F}) G dV’ $$
最终我们就得到了:
$$ \boxed{ \mathbf{F}(\mathbf{r}) = -\nabla \left( \int_V \frac{\nabla’ \cdot \mathbf{F}(\mathbf{r}’)}{4\pi |\mathbf{r} - \mathbf{r}’|} dV’ \right) + \nabla \times \left( \int_V \frac{\nabla’ \times \mathbf{F}(\mathbf{r}’)}{4\pi |\mathbf{r} - \mathbf{r}’|} dV’ \right) } $$
七、总结
我们通过以下步骤完成了推导:
- 利用格林函数表示单位脉冲;
- 将 $\mathbf{F}(\mathbf{r})$ 表示为积分形式;
- 应用矢量恒等式 $\nabla^2 \mathbf{F} = \nabla(\nabla \cdot \mathbf{F}) - \nabla \times (\nabla \times \mathbf{F})$;
- 使用分部积分将 $\nabla’$ 转移到 $G$ 上;
- 得出最终恒等式。
八、附注
- 如果 $\mathbf{F}(\mathbf{r})$ 定义在有限区域,则需要考虑边界项。
- 在物理中,例如电磁学和流体力学中,该恒等式被广泛用于构造电场、磁场、速度场的势函数。
如果你想要我用具体例子(如某个向量场)来验证这个恒等式是否成立,或者想看到在特定坐标系下的推导(比如球坐标或柱坐标),我也可以继续为你补充!
一、球坐标系简介
在球坐标系 $(r, \theta, \varphi)$ 中:
- $r$:径向距离
- $\theta$:极角(与 z 轴夹角)
- $\varphi$:方位角(绕 z 轴旋转)
单位基向量为: $$ \hat{\mathbf{e}}_r,\quad \hat{\mathbf{e}}_\theta,\quad \hat{\mathbf{e}}_\varphi $$
任意向量场可表示为: $$ \mathbf{F}(r, \theta, \varphi) = F_r(r, \theta, \varphi), \hat{\mathbf{e}}_r + F_\theta(r, \theta, \varphi), \hat{\mathbf{e}}_\theta + F_\varphi(r, \theta, \varphi), \hat{\mathbf{e}}_\varphi $$
二、亥姆霍兹分解回顾
我们希望将向量场 $\mathbf{F}$ 分解为两个部分:
$$ \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A} $$
其中:
- $\phi$ 是标量势(scalar potential),使得 $-\nabla \phi$ 是无旋场;
- $\mathbf{A}$ 是矢量势(vector potential),使得 $\nabla \times \mathbf{A}$ 是无散场。
三、球坐标系中的梯度、旋度和格林函数
1. 梯度(Gradient):
$$ \nabla f = \frac{\partial f}{\partial r}, \hat{\mathbf{e}}_r + \frac{1}{r} \frac{\partial f}{\partial \theta}, \hat{\mathbf{e}}_\theta + \frac{1}{r \sin \theta} \frac{\partial f}{\partial \varphi}, \hat{\mathbf{e}}_\varphi $$
2. 旋度(Curl):
对于矢量场 $\mathbf{A} = A_r \hat{\mathbf{e}}_r + A_\theta \hat{\mathbf{e}}_\theta + A_\varphi \hat{\mathbf{e}}_\varphi$,有:
$$ \nabla \times \mathbf{A} = \frac{1}{r \sin \theta} \left( \frac{\partial (A_\varphi \sin \theta)}{\partial \theta} - \frac{\partial A_\theta}{\partial \varphi} \right) \hat{\mathbf{e}}_r + \frac{1}{r} \left( \frac{1}{\sin \theta} \frac{\partial A_r}{\partial \varphi} - \frac{\partial (r A_\varphi)}{\partial r} \right) \hat{\mathbf{e}}_\theta + \frac{1}{r} \left( \frac{\partial (r A_\theta)}{\partial r} - \frac{\partial A_r}{\partial \theta} \right) \hat{\mathbf{e}}_\varphi $$
四、球坐标系下的亥姆霍兹分解公式
我们要构造 $\phi$ 和 $\mathbf{A}$,使得:
$$ \mathbf{F}(r, \theta, \varphi) = -\nabla \phi(r, \theta, \varphi) + \nabla \times \mathbf{A}(r, \theta, \varphi) $$
使用格林函数法构造势函数:
标量势 $\phi$ 的表达式:
$$ \phi(r, \theta, \varphi) = \frac{1}{4\pi} \int_{V’} \frac{\nabla’ \cdot \mathbf{F}(r’, \theta’, \varphi’)}{|\mathbf{r} - \mathbf{r}’|} r’^2 \sin \theta’, dr’ d\theta’ d\varphi’ $$
其中:
- 积分变量是 $(r’, \theta’, \varphi’)$
- 体积元 $dV’ = r’^2 \sin \theta’, dr’ d\theta’ d\varphi'$
- $|\mathbf{r} - \mathbf{r}’|$ 是两点之间的距离,在球坐标中为: $$ |\mathbf{r} - \mathbf{r}’| = \sqrt{r^2 + r’^2 - 2rr’(\sin\theta \sin\theta’ \cos(\varphi - \varphi’) + \cos\theta \cos\theta’)} $$
矢量势 $\mathbf{A}$ 的表达式:
$$ \mathbf{A}(r, \theta, \varphi) = \frac{1}{4\pi} \int_{V’} \frac{\nabla’ \times \mathbf{F}(r’, \theta’, \varphi’)}{|\mathbf{r} - \mathbf{r}’|} r’^2 \sin \theta’, dr’ d\theta’ d\varphi’ $$
五、完整形式的亥姆霍兹分解(球坐标系)
最终得到球坐标系下的亥姆霍兹分解表达式:
$$ \boxed{ \begin{aligned} \mathbf{F}(r, \theta, \varphi) &= -\nabla \left[ \frac{1}{4\pi} \int_{V’} \frac{\nabla’ \cdot \mathbf{F}(r’, \theta’, \varphi’)}{|\mathbf{r} - \mathbf{r}’|} r’^2 \sin \theta’, dr’ d\theta’ d\varphi’ \right] \\ &\quad + \nabla \times \left[ \frac{1}{4\pi} \int_{V’} \frac{\nabla’ \times \mathbf{F}(r’, \theta’, \varphi’)}{|\mathbf{r} - \mathbf{r}’|} r’^2 \sin \theta’, dr’ d\theta’ d\varphi’ \right] \end{aligned} } $$
六、物理意义与应用
这个公式说明了:
- 在球坐标系中,任何光滑且衰减良好的向量场都可以唯一地写成一个无旋场(由散度源决定)和一个无散场(由旋度源决定)之和。
- 它广泛应用于:
- 电磁学中电场和磁场的势描述;
- 天体物理中天体引力场或流体速度场的分解;
- 量子力学中角动量算符的处理;
- 流体力学中球对称流动问题。
七、特殊情形:轴对称场($\varphi$无关)
若向量场不依赖于方位角 $\varphi$,即具有轴对称性,则积分可以简化为关于 $r’$ 和 $\theta’$ 的二维积分,计算会更高效。
八、总结
操作 | 表达式 |
---|---|
标量势 $\phi$ | $\displaystyle \phi = \frac{1}{4\pi} \iiint \frac{\nabla’ \cdot \mathbf{F}}{ |
矢量势 $\mathbf{A}$ | $\displaystyle \mathbf{A} = \frac{1}{4\pi} \iiint \frac{\nabla’ \times \mathbf{F}}{ |
亥姆霍兹分解 | $\displaystyle \mathbf{F} = -\nabla \phi + \nabla \times \mathbf{A}$ |
如果你需要我用具体例子(例如某个球对称的向量场)来演示如何在球坐标系下进行亥姆霍兹分解,或者想看到数值实现方法,我可以继续为你详细展开!