什么是 旋度的旋度？
在向量微积分（Vector Calculus）中，旋度（Curl）是描述矢量场局部旋转性质的微分算子。当我们对一个矢量场先求旋度，再对结果求旋度时，就得到了“旋度的旋度”。这一运算在物理学，特别是电磁学、流体力学和弹性力学中扮演着核心角色。
对于任意二阶连续可微的矢量场 A，其旋度的旋度遵循以下著名的矢量恒等式：
其中：
- ∇ 是 Nabla 算子（矢量微分算子）。
- × 表示向量叉积（Cross Product）。
- · 表示向量点积（Dot Product）。
- ∇² 是矢量拉普拉斯算子（Vector Laplacian）。
这个公式不仅简洁优美，而且深刻揭示了矢量场的结构：一个矢量场的双旋度可以分解为该场的“散度梯度部分”与“拉普拉斯部分”之差。接下来，我们将通过两种主要方法进行严格推导。
旋度的旋度公式推导过程
推导 旋度的旋度 公式主要有两种方法：一种是基于矢量三重积的直观推导，另一种是基于爱因斯坦求和约定和 Levi-Civita 符号的严谨张量推导。我们将同时展示这两种方法，以确保理解的深度。
基于向量三重积恒等式的推导
我们首先回顾向量代数中的一个重要恒等式，即向量三重积公式：
A × (B × C) = B(A · C) - C(A · B)
现在，我们将向量 A 替换为微分算子 ∇，将 B 和 C 分别替换为 ∇ 和矢量场 A。注意，这里的 ∇ 是一个算子，它必须作用于其右侧的量。
步骤 1：形式上的替换
如果忽略 ∇ 的微分性质，直接套用公式，我们会得到：
∇ × (∇ × A) ≈ ∇(∇ · A) - A(∇ · ∇)
但这并不完全正确，因为 ∇ 是微分算子，必须使用莱布尼茨法则（乘积法则）。
步骤 2：应用莱布尼茨法则
在向量三重积 A × (B × C) 中，A 和 B 都作用于 C。在我们的例子中，第一个 ∇ 作用于第二个 ∇ 和 A。我们需要区分 ∇ 是只作用于 A 还是同时作用于 ∇ 和 A。
我们可以将 ∇ × (∇ × A) 写为两项之和：
- 项 1：第一个 ∇ 只作用于 A。
- 项 2：第一个 ∇ 只作用于第二个 ∇（这在数学上通常为零或需特殊处理，但在标准矢量恒等式推导中，我们主要关注 ∇ 作用于 A 的部分以及交叉项）。
更严谨的说法是，利用恒等式：
∇ × (A × B) = A(∇ · B) - B(∇ · A) + (B · ∇)A - (A · ∇)B
令 A = ∇，B = A（注意符号重载，这里 B 是矢量场 A）：
∇ × (∇ × A) = ∇(∇ · A) - A(∇ · ∇) + (A · ∇)∇ - (∇ · ∇)A
这里有一些项需要仔细辨析。由于 ∇ 是微分算子，(A · ∇)∇ 这一项通常被视为零，因为算子 ∇ 不能对自己求导（或者说，二阶导数的对称性使得某些交叉项抵消）。更标准的处理方式是直接使用分量法或 Levi-Civita 符号来避免歧义。
让我们采用更被广泛接受的简化逻辑：对于矢量场 A，我们有恒等式：
∇ × (∇ × A) = ∇(∇ · A) - ∇²A
这个结果可以通过张量分析严格证明，如下一个标签页所示。
基于 Einstein 求和约定与 Levi-Civita 符号
这是最严谨且通用的推导方法，适用于任意正交坐标系。
定义
- Nabla 算子的第 i 个分量：∂i
- 矢量场 A 的第 i 个分量：Ai
- Levi-Civita 符号：εijk
- Kronecker delta：δij
推导步骤
我们要计算矢量 B = ∇ × (∇ × A) 的第 i 个分量 Bi。
第一步，写出内层旋度的分量：
(∇ × A)_j = ε_jmn ∂_m A_n
第二步，写出外层旋度的第 i 个分量：
B_i = ε_ijk ∂_j (∇ × A)_k
代入内层结果：
B_i = ε_ijk ∂_j (ε_klm ∂_l A_m)
由于 εklm 是常数，可以提到微分符号外面：
B_i = ε_ijk ε_klm ∂_j ∂_l A_m
利用恒等式 εijk = εkij，重写为：
B_i = ε_kij ε_klm ∂_j ∂_l A_m
使用 Levi-Civita 的收缩恒等式：εkijεklm = δilδjm - δimδjl：
B_i = (δ_il δ_jm - δ_im δ_jl) ∂_j ∂_l A_m
展开这两项：
B_i = δ_il δ_jm ∂_j ∂_l A_m - δ_im δ_jl ∂_j ∂_l A_m
利用 Kronecker delta 的性质（δxy 仅在 x=y 时为 1，否则为 0）进行缩并：
- 第一项：δil 将 l 替换为 i，δjm 将 m 替换为 j。得到：∂j∂iAj
- 第二项：δim 将 m 替换为 i，δjl 将 l 替换为 j。得到：∂j∂jAi
因此：
B_i = ∂_i (∂_j A_j) - ∂_j ∂_j A_i
识别各项的物理意义：
- ∂jAj 是散度 ∇ · A。所以 ∂i(∂jAj) 是梯度 ∇(∇ · A) 的第 i 个分量。
- ∂j∂j 是拉普拉斯算子 ∇²。所以 ∂j∂jAi 是矢量拉普拉斯 ∇²A 的第 i 个分量。
最终得到：
B_i = [∇(∇ · A)]_i - [∇²A]_i
即：
∇ × (∇ × A) = ∇(∇ · A) - ∇²A
笛卡尔坐标系分量展开法
对于初学者，直接展开 x, y, z 分量可能更直观，虽然计算繁琐，但能避免抽象符号的混淆。
设 A = (Ax, Ay, Az)。
1. 计算 ∇ × A
(∇ × A)_x = ∂_y A_z - ∂_z A_y
(∇ × A)_y = ∂_z A_x - ∂_x A_z
(∇ × A)_z = ∂_x A_y - ∂_y A_x
2. 计算 ∇ × (∇ × A) 的 x 分量
根据旋度定义：
[∇ × (∇ × A)]_x = ∂_y (∇ × A)_z - ∂_z (∇ × A)_y
代入第一步的结果：
= ∂_y (∂_x A_y - ∂_y A_x) - ∂_z (∂_z A_x - ∂_x A_z)
展开微分：
= ∂_y ∂_x A_y - ∂_y ∂_y A_x - ∂_z ∂_z A_x + ∂_z ∂_x A_z
重组项，将含有 A_x 的项放在一起，其他项放在一起：
= - (∂_y ∂_y A_x + ∂_z ∂_z A_x) + (∂_x ∂_y A_y + ∂_x ∂_z A_z)
为了凑出完整的拉普拉斯算子 ∇²Ax = (∂_x² + ∂_y² + ∂_z²)Ax，我们在式子中加一项减一项 ∂_x ∂_x Ax：
= - (∂_x² A_x + ∂_y² A_x + ∂_z² A_x) + (∂_x² A_x + ∂_x ∂_y A_y + ∂_x ∂_z A_z)
提取公因式：
= - ∇² A_x + ∂_x (∂_x A_x + ∂_y A_y + ∂_z A_z)
识别括号内为散度 ∇ · A：
= - ∇² A_x + ∂_x (∇ · A)
这正是 [∇(∇ · A) - ∇²A]x 的第 x 分量。y 和 z 分量同理可证。
旋度的旋度的物理意义
理解公式的数学推导只是第一步，深入理解其物理内涵对于应用至关重要。
1. 场的分解
亥姆霍兹定理（Helmholtz Theorem）指出，任何充分光滑且衰减足够快的矢量场都可以分解为一个无旋场（梯度场）和一个无散场（旋度场）之和。旋度的旋度公式实际上展示了如何通过散度和拉普拉斯算子来重构或分析矢量场的结构。
2. 波动方程的核心
在电磁学中，麦克斯韦方程组是基础。通过对法拉第定律 ∇ × E = -∂B/∂t 取旋度，并利用 ∇ × (∇ × E) 恒等式，我们可以将电场 E 的方程转化为包含 ∇²E 的形式，从而推导出电磁波动方程，证明光是一种电磁波。
3. 粘性流体的扩散
在流体力学中，纳维-斯托克斯方程（Navier-Stokes Equations）描述了粘性流体的运动。其中的粘性项 ν∇²v 直接与 旋度的旋度 有关（对于不可压缩流体，∇·v=0，故 ∇×(∇×v) = -∇²v）。这表明粘性力本质上是速度场旋度的扩散。
网友们还关心的 旋度 周边知识
除了公式推导，用户通常还会关注 旋度 与其他矢量运算的关系，以及在实际工程中的具体计算技巧。以下是高频关注的热点话题。
矢量微积分算子对比表
| 运算类型 | 算子表示 | 输入 | 输出 | 物理意义 |
|---|---|---|---|---|
| 梯度 (Gradient) | ∇f | 标量场 | 矢量场 | 方向导数最大的方向，表示场的变化率 |
| 散度 (Divergence) | ∇ · F | 矢量场 | 标量场 | 源或汇的强度，表示场的发散程度 |
| 旋度 (Curl) | ∇ × F | 矢量场 | 矢量场 | 场的旋转程度，表示涡旋的强度 |
| 拉普拉斯 (Laplacian) | ∇²f | 标量场 | 标量场 | 函数值的平均值与中心值的偏差 |
| 矢量拉普拉斯 | ∇²F | 矢量场 | 矢量场 | 各分量的拉普拉斯算子组合 |
常见误区与注意事项
- 算子顺序： 在使用 ∇ 算子时，必须注意其微分属性。例如，∇(AB) ≠ A∇B + B∇A，而应遵循乘积法则 ∇(AB) = A(∇B) + B(∇A)（当 A,B 为标量时）。在 旋度的旋度 推导中，忽略这一点会导致错误。
- 坐标系选择： 上述恒等式在笛卡尔坐标系下最直观。在柱坐标或球坐标系下，虽然恒等式形式不变，但 ∇ 算子的具体表达式会变得非常复杂，计算时应谨慎。
- 无旋场与无散场： 如果 ∇ × A = 0，则 A 是保守场，可以表示为某个标量势的梯度（A = ∇φ）。此时 ∇ × (∇ × A) = 0。如果 ∇ · A = 0，则 A 是无散场，此时 ∇ × (∇ × A) = -∇²A。
历史沿革：矢量分析的发展
吉布斯与亥姆霍兹的贡献
约西亚·威拉德·吉布斯（Willard Gibbs）和奥利弗·亥姆霍兹（Helmholtz）等人独立发展了现代矢量分析方法，将麦克斯韦复杂的四元数方程简化为易于操作的矢量形式，确立了 ∇ 算子的地位。
麦克斯韦方程组的矢量化
麦克斯韦原始的方程组包含大量四元数运算。矢量分析的引入使得电磁理论变得简洁明了，旋度的概念成为理解电磁波传播的关键。
广义相对论与微分几何
在广义相对论中，旋度的概念被推广为外微分（Exterior Derivative）和微分形式。虽然形式更加抽象，但其核心思想——描述场的局部旋转和拓扑性质——依然保留。
关于 旋度的旋度 的常见问题
A: 这是因为梯度的旋度恒为零。从物理上看，梯度场是保守场，沿闭合路径的环量为零，因此没有“旋转”。从数学上看，∂²φ/∂x∂y = ∂²φ/∂y∂x（混合偏导数相等），导致旋度分量相互抵消。
A: 这是散度的旋度恒为零。旋度产生的矢量场总是无散的（无源无汇），因为旋度描述的是闭合的涡旋线，涡旋线不会开始或结束于某一点，因此净流出通量为零。
A: 通常不直接计算双旋度。而是先数值计算速度场的旋度（涡量 ω = ∇ × v），然后再对涡量场计算旋度。或者，直接利用恒等式，先计算散度和拉普拉斯项。有限体积法（FVM）或有限差分法（FDM）是常用的数值离散手段。
A: 在弯曲空间（黎曼流形）中，标准的 Nabla 算子需要替换为协变导数（Covariant Derivative）。此时，∇ × (∇ × A) 的表达式会引入黎曼曲率张量（Riemann Curvature Tensor）的项，不再简单地等于 ∇(∇·A) - ∇²A。曲率会干扰矢量的平行移动和旋转。
旋度的旋度公式 ∇×(∇×A) = ∇(∇·A) - ∇²A 是向量微积分中连接局部旋转、发散和扩散性质的桥梁。无论是通过矢量恒等式的直观推导，还是张量分析的严谨证明，该公式都展现了数学结构的统一美。在电磁学、流体力学等物理领域，它是推导波动方程和运动方程不可或缺的工具。掌握这一公式，对于深入理解矢量场理论具有重要意义。