散度、旋度与三维向量场 | 多变量微积分 I | 自在学散度、旋度与三维向量场
这一章只问一个局部问题：一个向量场在某个很小的区域附近，是更像“有东西冒出来或流进去”，还是更像“让小转轮自转”？
前一章把向量场放到曲面上，计算穿过曲面的通量。本章把视角收回到一个点附近，定义两个从偏导数来的量：散度和旋度。它们不是新的积分技巧，而是描述向量场局部行为的语言。
源场净流出为正，汇场净流出为负，旋转场表现出明显的局部转动。
散度回答“净流出多少”，旋度回答“局部怎样转”。这两个问题看起来都很直观，但真正计算时都要用局部极限。只看一张箭头图，常常会被箭头长短或整体绕行骗过去。
三维向量场和算子
三维向量场把空间中的每个点对应到一个三维向量。通常写成
F(x,y,z)=P(x,y,z)i+
Q(x,y,z)j+
R(x,y,z)k
F=(P,Q,R) 其中 P、Q、R 分别是 x、y、z 方向的分量函数。比如流体速度场中，F(2,1,0) 表示空间点 (2,1,0) 处流体的速度方向和大小。
向量场的三个分量与 ∇ 算子对应的方向变化。
∇=⟨∂x∂,∂y∂,∂z∂⟩ 这个符号本身不是普通向量。它像一个“偏导数工具箱”：和标量函数配合时给出梯度，和向量场用点乘形式配合时给出散度，和向量场用叉乘形式配合时给出旋度。
把 ∇ 当成普通向量做形式运算，可以帮助记公式；但每一项最后都要回到偏导数。散度和旋度都是由局部分量的变化率组成的，不是从箭头图上直接数出来的。
散度：局部净流出
F=(P,Q,R) ∇⋅F=∂x∂P+∂y∂Q+∂z∂R 如果把 F 看成流体速度场，那么散度描述一个很小体积附近的净流出率。正散度表示流出多于流入，负散度表示流入多于流出，零散度表示局部没有净源汇。
微小体积 ΔV 中流出多于流入时，向量场在该点表现为正散度。
这里的“净”很重要。一个小盒子的某一面可以有很强的流出，但另一面也可能有同样强的流入。散度比较的是所有面的总效果，而不是某一个方向上的箭头大小。
∇⋅F=ΔV→0limΔV1∬∂BF⋅ndS 这里 B 是包围该点的一个很小体积，∂B 是它的边界曲面，n 是向外法向量。这个式子现在只作为直观定义使用；下一章会把它发展成闭曲面通量和体积分之间的定理。
一个最直接的判断
F=(x,y,z) ∇⋅F=∂x∂x+∂y∂y+∂z∂z=3 所以每个点附近都有正的局部净流出。它像均匀膨胀的流场：不是只有原点在“冒出东西”，而是空间每一点附近都在膨胀。
F=(−y,x,0) ∇⋅F=∂x∂(−y)+∂y∂x+∂z∂0=0 这个场绕着 z 轴转，但没有局部净源汇。一个小盒子某边流入的量，会被另一边流出的量抵消。
旋度：局部旋转轴
∇×F=(∂y∂R−∂z∂Q,∂z∂P−∂x∂R,∂x∂Q−∂y∂P) ∇×F=i∂x∂Pj∂y∂Qk∂z∂R 旋度是一个向量。它的方向表示局部旋转轴的方向，大小表示围绕这个方向旋转的强度。用右手法则判断方向：手指按局部旋转方向弯曲，拇指所指就是旋度方向。
小转轮直观展示局部旋转；按右手法则，逆时针转动对应向上的旋度方向。
如果把一个很小的转轮放到流场里，它的转轴方向就是旋度方向的一种直观读法。更准确地说，旋度在单位法向量 n 上的分量，来自该小面片边界的单位面积环流：
(∇×F)⋅n=ΔA→0limΔA1∮∂SF⋅Tds 这里 T 是边界的正向单位切向量。它说明旋度不是“整体路径绕了一圈”的代名词，而是点附近很小环路上的环流密度。
典型场的旋度
F=(−ωy,ωx,0) ∇×F=(0,0,2ω) 当 ω>0 时，从 z 轴正方向看，场逆时针旋转，旋度指向 +z 方向。这里的系数 2 不要忽略：在流体力学的常用解释中，旋度等于微小刚体旋转角速度的两倍。
F=(x,y,z) ∇×F=(0,0,0) 它到处往外发散，但不会让小转轮自转。这个例子提醒我们：散度和旋度测的是两种不同局部行为。
先算再解释
很多题目会同时要求散度和旋度。稳妥做法是先把 P,Q,R 写清楚，再分别代入两个公式，最后用语言解释结果。
例题：一个三维多项式向量场
F=(x2y,yz,z2) 求 ∇⋅F 和 ∇×F，并解释点 (1,−2,1) 附近的局部行为。
先确定三个分量。这里 P=x2y，Q=yz，R=z2。散度只取每个分量对自己方向的偏导数。
∇⋅F=∂x∂(x再按旋度公式逐项计算。第一分量比较 R 对 y 的变化和 Q 对 z 的变化，第二分量比较 P 对 的变化和 对 的变化，第三分量比较 对 的变化和 对 的变化。
把点 (1,−2,1) 代入散度。
∇⋅F(1,−2,把点代入旋度。
∇×F(1,−2,1)=(2,0,−1)这说明局部旋转轴朝向 。如果把很小的转轮放在这个点附近，最明显的转轴不是竖直方向，而是偏向正 、负 的方向。
散度的结果是标量，它只告诉你净源汇强弱；旋度的结果是向量，它还带有旋转轴方向。把这两类结果混成同一种量，是本章最常见的计算后解释错误。
例题：用结构快速判断
F=(ay,bx,0) ∇⋅F=∂x∂(ay)+∂y∂(bx)+∂z∂0=0 ∇×F=(0,0,b−a) 这个例子很适合训练直觉：ay 是水平速度却随 y 改变，bx 是竖直速度却随 x 改变，这种交叉变化会制造局部旋转；但它们不一定制造源汇。
与 Green 定理的关系
二维 Green 定理其实已经把“局部变化”和“边界累积”的关系讲过一遍。现在只是把语言换成散度和旋度。
F=(P(x,y),Q(x,y),0) ∇×F=(0,0,∂x∂Q−∂y∂P) ∮CPdx+Qdy=∬D(∂x∂Q−∂y∂P)dA 就是在说：边界上的总环流等于内部局部旋转的面积累积。
∇⋅F=∂x∂P+∂y∂Q ∮CPdy−Qdx=∬D(∂x∂P+∂y∂Q)dA 这句话对应另一种边界观察：边界上的向外通量等于内部局部源汇的面积累积。
从格林定理的二维边界累积，过渡到三维向量场中散度与旋度的局部解释。
本章不需要提前使用 Stokes 定理或散度定理。你只要先看懂 Green 定理的两种版本：环流对应旋度，通量对应散度。下一章的两个三维定理，就是把这个“边界累积等于内部局部量累积”的结构推广到曲面和空间区域。
常见误区
绕着走不一定处处有旋度；有发散箭头也要看闭区域内源与汇的净效应。
误区一：箭头变长就一定散度为正
散度看的是各分量沿自身方向的局部变化。箭头在某个方向变长，可能只说明速度大小变大，不一定说明小体积有净流出。
F=(0,x) ∇⋅F=∂x∂0+∂y∂x=0 误区二：绕着某点走就一定旋度非零
旋度看局部小环路的环流密度，不是看大路径有没有绕圈。某些场会沿着大圆环绕行，但在不包含奇点的区域内，局部小转轮并不自转。
F=(x2+y2−y,x2+y2x) 在原点以外，它的二维旋度分量为 0。但沿着包围原点的闭曲线，环流并不为 0。这里的关键是原点不在定义域内，局部判断和全局路径判断不能混为一谈。
误区三：散度为零就表示没有运动
散度为零只表示局部没有净源汇。流体仍然可以平移、剪切或旋转。刚体旋转场 F=(−y,x,0) 的散度为零，但旋度不为零。
误区四：旋度为零就一定能随便找势函数
在简单、没有洞的区域中，旋度为零通常和保守场联系在一起。但如果定义域有洞，局部旋度为零不一定能保证所有闭合路径的线积分都为零。遇到这类题目，要同时检查向量场的定义域。
散度和旋度是局部量，线积分和通量积分常常是区域或边界上的整体量。两者可以通过 Green 定理、Stokes 定理和散度定理连接，但不能在没有条件的情况下直接互相替代。
练习
练习一
F=(3x,−2y,z) 的散度和旋度，并解释它是否有局部净源汇和局部旋转。
散度为
∇⋅F=3−2+1=2旋度为
∇×F=(0,0,0)所以它有正的局部净流出，但没有局部旋转。这个场沿三个坐标方向伸缩，其中 x 和 z 方向的扩张总量大于 y 方向的压缩。
练习二
F=(−y,x,z) 散度为
∇⋅F=0+0+1=1旋度为
∇×F=(0,0,2)它同时有正的局部净流出和绕 z 轴正向的局部旋转。前两个分量造成平面旋转，第三个分量 z 造成竖直方向的伸长。
练习三
F=(ay,bx,0) 找出使 ∇×F=(0,0,3) 的 a,b 条件。它的散度是多少？
散度恒为
∇⋅F=0旋度为
∇×F=(0,0,b−a)所以需要
b−a=3例如 a=−23、b=23，得到刚体旋转场 。满足条件的组合不止这一组。
练习四
判断下列说法是否正确：若一个三维向量场的散度为 0，则它的旋度也一定为 0。
不正确。反例是
F=(−y,x,0)它的散度为
∇⋅F=0但旋度为
∇×F=(0,0,2)散度和旋度测量的是不同的局部行为，一个为零不能推出另一个为零。
练习五
F=(x2,xy,yz) 求 ∇⋅F 和 ∇×F。
散度为
∇⋅F=2x+x+y=3x+y旋度为
∇×F=(∂y∂(yz)−∂z∂(xy),计算时最容易漏掉的是第三个分量：∂x∂(xy)=y，而 。
小结
散度和旋度把三维向量场的局部行为拆成两个问题。散度是标量，描述小体积的净流出；旋度是向量，描述小转轮的旋转轴和强度。
∇⋅(P,Q,R)=Px+Qy+Rz ∇×(P,Q,R)=(Ry−Qz,Pz−Rx,Qx−Py) 理解时把它们放回边界思想：环流对应旋度，通量对应散度。下一章的 Stokes 定理和散度定理，会把这种局部和整体的对应关系推到三维曲面与空间区域上。
2
y
)
+
∂y∂(yz)+
∂z∂(z2)=
2xy+
z+
2z=
2xy+
3z
z
∇×F=(∂y∂z2−∂z∂(yz),∂z∂(x2y)−∂x∂z2,∂x∂(yz)−∂y∂(x2y))=(−y,0,−x2) 1
)
=
2(1)(−2)+
3(1)=
−1
该点附近的净效果是流入多于流出，因此它表现为弱的局部汇。
(2,0,−1)
F=(−23y,23x,0) ∂z∂(x2)
−
∂x∂(yz)
,
∂x∂(xy)
−
∂y∂(x2)
)
=
(z,0,y)
∂
y
∂(x2)
=
0
