梯度、散度与旋度：场怎样在空间中变化

从一张地形图开始

地形高度 $f(x,y)$ 是标量场。梯度 $\nabla f$ 指向高度增长最快的方向；沿单位方向 $\mathbf n$ 的变化率为 $\nabla f\cdot\mathbf n$。因此梯度既给方向，也给变化强度。注意梯度的单位比原场多除以一个长度。

散度作用于向量场，问一个极小体积的净流出通量。旋度问一个极小回路的单位面积环流。它们都由同一个符号 $\nabla$ 组成，却回答不同的问题：

两个小场，把概念分开

先看 $\mathbf F=(x,y,0)$。箭头像从中心向外伸展，散度为 2，旋度为零。再看 $\mathbf G=(-y,x,0)$。箭头围绕中心转，散度为零，而旋度为 $(0,0,2)$。图像中的旋转和发散可以分别存在，也可以同时存在。

这两例只是无量纲教学场。真实磁场满足 $\nabla\cdot\mathbf B=0$，不能把第一例直接当成实际磁感应强度场。第二例可以帮助理解局部旋度，却也不能不检查来源与边界就当成一套完整物理解。

Laplacian 为什么像“和邻居比较”

一维二阶导数近似为 $[f(x+h)-2f(x)+f(x-h)]/h^2$。若某点高于两边，二阶导数为负；扩散方程便让它下降。若低于两边，它会上升。三维 Laplacian $\nabla^2 f$ 把三个方向的二阶导数加起来。

在 Cartesian 坐标和足够光滑的区域，这个恒等式能把两次旋度变成 Laplacian。磁扩散的推导依赖它；材料界面处应使用边界条件，不能直接把不连续系数当成常数移出微分。

一个带单位的计算

若温度 $T=300+2x^2$ K，且 $x$ 用米计，则在 $x=1$ m 处 $dT/dx=4$ K/m，$d^2T/dx^2=4$ K/m²。相同数学运算用于电磁场时，场的单位会改变，但“空间变化率”的含义不变。

把微分算子理解为局部实验

梯度考察向哪个方向变化最快，散度考察一个小体积的净流出，旋度考察小回路上的环流。三个概念分别处理标量到向量、向量到标量、向量到向量的变化。用量纲可以立即排除一些错误：若 B 的单位是 T，空间导数就是 T/m，二阶空间导数就是 T/m²。

对标量 $f=x^2+2y^2+3z^2$，直接求偏导得到

在点 (1,0,0)，向 x 方向移动一小步 h，f 增加约 2h。拉普拉斯算子则比较中心与周围的平均差异：一个局部凸起在扩散中趋于降低，一个凹陷趋于填平。

一个只有旋转、没有净流出的场

令 $\mathbf F=(-y,x,0)$，则

绕原点半径 r 的圆，场沿切向、大小为 r，因此环流为 $2\pi r^2$；旋度在圆面上的通量也为 $2\pi r^2$。这个例子帮助我们看到：散度为零绝不意味着场本身为零，更不意味着没有旋转。

扩散推导真正依赖的恒等式

磁通守恒把第一项变为零，才得到简洁的向量拉普拉斯形式。若数值网格没有保持离散磁通约束，额外的梯度项可能表现为非物理误差。这里各分量的拉普拉斯在笛卡尔坐标中可以分别求；换到球坐标时，基向量也随位置变化，不能把标量球坐标公式直接套给每一个向量分量。

研究拓展

从离散采样计算导数会放大噪声。网格更密不一定让实测差分更可靠，因为分母也变小。数值求解中应尽量保持离散散度、旋度之间的结构，而非任意拼接差分。

练习与解析