我们首先介绍矢量微分算子 ∇ (Nabla 或 Del 算子)。在三维笛卡尔坐标系中,它的定义是:
∇=i∂x∂+j∂y∂+k∂z∂
其中 i,j,k 分别是 x,y,z 方向的单位矢量。
下面我们将分别推导梯度、散度和旋度的计算公式及其物理意义。
一、梯度 (Gradient)#
梯度作用于一个标量场 f(x,y,z),结果是一个矢量场。
1、定义与推导#
考虑一个标量场 f(x,y,z)。我们想知道这个标量场在空间中某一点附近的变化情况。根据多元函数微分学的知识,标量场 f 的全微分(total differential)为:
df=∂x∂fdx+∂y∂fdy+∂z∂fdz这表示当位置从 (x,y,z) 移动一个微小位移 (dx,dy,dz) 时,f 值的微小变化量。
我们可以将上式看作两个矢量的点积。定义一个矢量,称为 f 的梯度,记作 ∇f 或 grad(f):
∇f=(∂x∂fi+∂y∂fj+∂z∂fk)再定义微小位移矢量 dr:
dr=dxi+dyj+dzk于是,全微分可以写成:
df=(∇f)⋅dr根据点积的定义,df=∣∇f∣∣dr∣cosθ,其中 θ 是 ∇f 和 dr 之间的夹角。
单位位移的变化率(即方向导数)为:
∣dr∣df=∣∇f∣cosθ这个变化率在 cosθ=1 时取得最大值,即当 θ=0 时。这意味着 dr 的方向与 ∇f 的方向一致。此时,最大的变化率为 ∣∇f∣。
所以,梯度 ∇f 的方向是标量场 f 变化率最大的方向,其大小是这个最大变化率。
计算公式 (笛卡尔坐标系):
∇f=∂x∂fi+∂y∂fj+∂z∂fk2、物理意义#
- 方向:梯度 ∇f 的方向指向标量场 f 在该点增加最快的方向。
- 大小:梯度 ∣∇f∣ 的大小表示标量场 f 在该点沿最快增加方向的单位距离变化率(即最大方向导数)。
- 等值面法线:对于一个等值面(或等值线)f(x,y,z)=C(常数),在该面上任意一点的微小位移 dr 都满足 df=0。由于 df=∇f⋅dr,所以 ∇f⋅dr=0。这意味着梯度 ∇f 与等值面在该点的切平面(或切线)上的任何矢量 dr 都正交。因此,梯度 ∇f 的方向是等值面的法线方向。
例子:
- 温度场 T(x,y,z):∇T 指向温度升高最快的方向,其大小表示温度变化率。热流密度通常与 −∇T 成正比(傅里叶热传导定律)。
- 电势场 V(x,y,z):−∇V 是电场强度 E。电场线指向电势降低最快的方向。
- 高度场 h(x,y):∇h 指向坡度最陡峭的向上方向。
二、散度 (Divergence)#
散度作用于一个矢量场 F(x,y,z)=Fxi+Fyj+Fzk,结果是一个标量场。
1、定义与推导#
散度的物理解释是矢量场在某一点的源的强度或汇的强度。一个更严格的定义是:
∇⋅F=ΔV→0limΔV1∫∫◯SF⋅dS其中 ΔV 是一个包含某点的微小体积元, S 是包围 ΔV 的闭合曲面,dS 是指向外侧的面积元矢量。这个积分表示通过闭合曲面 S 的矢量场 F 的总通量。散度就是单位体积的净流出通量。
我们考虑一个在点 (x,y,z) 附近,边长分别为 Δx,Δy,Δz 的微小长方体体积元 ΔV=ΔxΔyΔz。长方体的中心近似为 (x,y,z)。
我们计算通过这个长方体六个面的通量。
-
考虑沿 x 方向的通量:
- 右侧面 (在 x+2Δx 处,面积 ΔyΔz,法向 i):
通量 Φx,right≈Fx(x+2Δx,y,z)ΔyΔz
- 左侧面 (在 x−2Δx 处,面积 ΔyΔz,法向 −i):
通量 Φx,left≈−Fx(x−2Δx,y,z)ΔyΔz
- 沿 x 方向的净通量 ΔΦx=Φx,right+Φx,left:
ΔΦx≈[Fx(x+2Δx,y,z)−Fx(x−2Δx,y,z)]ΔyΔz
对于微小的 Δx,方括号内的差值可以近似为(根据中值定理或泰勒展开):
Fx(x+2Δx,y,z)−Fx(x−2Δx,y,z)≈∂x∂FxΔx
所以,
ΔΦx≈∂x∂FxΔxΔyΔz=∂x∂FxΔV
-
类似地,沿 y 方向的净通量:
- 上侧面 (在 y+2Δy 处,法向 j): Fy(x,y+2Δy,z)ΔxΔz
- 下侧面 (在 y−2Δy 处,法向 −j): −Fy(x,y−2Δy,z)ΔxΔz
- 净通量 ΔΦy≈∂y∂FyΔxΔyΔz=∂y∂FyΔV
-
沿 z 方向的净通量:
- 前侧面 (在 z+2Δz 处,法向 k): Fz(x,y,z+2Δz)ΔxΔy
- 后侧面 (在 z−2Δz 处,法向 −k): −Fz(x,y,z−2Δz)ΔxΔy
- 净通量 ΔΦz≈∂z∂FzΔxΔyΔz=∂z∂FzΔV
总的净通量
∫∫◯SF⋅dS=ΔΦx+ΔΦy+ΔΦz∫∫◯SF⋅dS≈(∂x∂Fx+∂y∂Fy+∂z∂Fz)ΔV根据散度的定义:
∇⋅F=ΔV→0limΔV1∫∫◯SF⋅dS=∂x∂Fx+∂y∂Fy+∂z∂Fz这也可以看作是 ∇ 算子与矢量场 F 的点积:
∇⋅F=(i∂x∂+j∂y∂+k∂z∂)⋅(Fxi+Fyj+Fzk)=∂x∂Fx+∂y∂Fy+∂z∂Fz计算公式 (笛卡尔坐标系):
∇⋅F=∂x∂Fx+∂y∂Fy+∂z∂Fz2、物理意义#
-
源与汇:散度描述了矢量场在某一点的“流出”或“汇入”的强度。
- ∇⋅F>0:该点是矢量场的“源”(source),表示有净流出。例如,流体在此处膨胀或有物质产生。
- ∇⋅F<0:该点是矢量场的“汇”(sink),表示有净流入。例如,流体在此处压缩或有物质消失。
- ∇⋅F=0:该点无源无汇,或者流入量等于流出量。这样的场称为无源场或螺线管场(solenoidal field)。例如,不可压缩流体的速度场(没有源或汇的情况下),磁场(∇⋅B=0,磁单极子不存在)。
-
高斯散度定理:散度与高斯散度定理密切相关,该定理指出:
∭V(∇⋅F)dV=∫∫◯SF⋅dS
这表示体积 V 内所有源(或汇)的总强度等于通过包围该体积的闭合曲面 S 的净通量。
例子:
- 流体速度场 v:∇⋅v 表示流体在某点的膨胀率或压缩率。对于不可压缩流体,∇⋅v=0。
- 电场 E:∇⋅E=ρ/ϵ0 (高斯定律的微分形式),其中 ρ 是电荷密度,ϵ0 是真空介电常数。电荷是电场的源。
- 磁场 B:∇⋅B=0 (高斯磁定律),表明不存在磁单极子。
三、旋度 (Curl)#
旋度作用于一个矢量场 F(x,y,z)=Fxi+Fyj+Fzk,结果是另一个矢量场。
1、定义与推导#
旋度描述了矢量场在某一点的“旋转”或“环流”的程度和方向。旋度矢量的方向是旋转轴的方向(遵循右手法则),其大小表示旋转的强度。
旋度 ∇×F 的任意方向 n^ 上的分量定义为:
(∇×F)⋅n^=ΔA→0limΔA1∮CF⋅dr其中 C 是一个以 n^ 为法线方向的微小闭合回路,ΔA 是该回路所围成的面积。积分 ∮CF⋅dr 称为矢量场 F 沿闭合路径 C 的环量。
我们推导旋度的 z 分量 (∇×F)z。考虑在 xy 平面内,以点 (x,y,z) 为中心的一个微小矩形回路,其顶点为 (x−2Δx,y−2Δy), (x+2Δx,y−2Δy), (x+2Δx,y+2Δy), (x−2Δx,y+2Δy)。面积 ΔA=ΔxΔy,法向为 k。我们沿逆时针方向计算环量。
- 路径 1 (底边): 从 (x−2Δx,y−2Δy) 到 (x+2Δx,y−2Δy)。dr=dxi。
∫1F⋅dr≈Fx(x,y−2Δy,z)Δx
- 路径 2 (右边): 从 (x+2Δx,y−2Δy) 到 (x+2Δx,y+2Δy)。dr=dyj。
∫2F⋅dr≈Fy(x+2Δx,y,z)Δy
- 路径 3 (顶边): 从 (x+2Δx,y+2Δy) 到 (x−2Δx,y+2Δy)。dr=−dxi。
∫3F⋅dr≈−Fx(x,y+2Δy,z)Δx
- 路径 4 (左边): 从 (x−2Δx,y+2Δy) 到 (x−2Δx,y−2Δy)。dr=−dyj。
∫4F⋅dr≈−Fy(x−2Δx,y,z)Δy
总环量 ∮CF⋅dr=∫1+∫2+∫3+∫4:
∮CF⋅dr≈[Fx(x,y−2Δy,z)−Fx(x,y+2Δy,z)]Δx+[Fy(x+2Δx,y,z)−Fy(x−2Δx,y,z)]Δy对于微小的 Δx,Δy:
- Fx(x,y−2Δy,z)−Fx(x,y+2Δy,z)≈−∂y∂FxΔy
- Fy(x+2Δx,y,z)−Fy(x−2Δx,y,z)≈∂x∂FyΔx
所以,
∮CF⋅dr≈(−∂y∂FxΔy)Δx+(∂x∂FyΔx)Δy=(∂x∂Fy−∂y∂Fx)ΔxΔy根据旋度分量的定义,(∇×F)z=(∇×F)⋅k:
(∇×F)z=ΔA→0limΔA1∮CF⋅dr=ΔxΔy(∂x∂Fy−∂y∂Fx)ΔxΔy=∂x∂Fy−∂y∂Fx通过对坐标进行轮换对称 (x→y→z→x),可以得到旋度的其他分量:
(∇×F)x=∂y∂Fz−∂z∂Fy(∇×F)y=∂z∂Fx−∂x∂Fz将这三个分量组合起来,得到旋度矢量:
∇×F=(∂y∂Fz−∂z∂Fy)i+(∂z∂Fx−∂x∂Fz)j+(∂x∂Fy−∂y∂Fx)k这可以形式上写成 ∇ 算子与矢量场 F 的叉积,通常用行列式表示:
∇×F=i∂x∂Fxj∂y∂Fyk∂z∂Fz计算公式 (笛卡尔坐标系):
∇×F=(∂y∂Fz−∂z∂Fy)i+(∂z∂Fx−∂x∂Fz)j+(∂x∂Fy−∂y∂Fx)k2、物理意义#
-
旋转性:旋度描述了矢量场在某一点的涡旋强度和旋转轴方向。
- 如果 ∇×F=0,则该场是有旋场 (rotational field)。想象在场中放置一个微小的桨叶轮,如果它开始旋转,则该点旋度非零。旋度矢量的方向是桨叶轮的旋转轴方向(由右手法则确定),其大小与旋转角速度成正比。
- 如果 ∇×F=0,则该场是无旋场 (irrotational field) 或保守场。这意味着场线不会形成闭合的小环路,或者说场中没有涡旋。对于无旋场,一定存在一个标量势 ϕ 使得 F=∇ϕ(如果场是单连通的)。
-
斯托克斯定理:旋度与斯托克斯定理密切相关,该定理指出:
∬S(∇×F)⋅dS=∮CF⋅dr
这表示通过一个开放曲面 S 的旋度通量等于矢量场 F 沿该曲面边界闭合路径 C 的环量。
例子:
- 流体速度场 v:∇×v 是流体的涡量 (vorticity) 的两倍,描述了流体微团的旋转。如果 ∇×v=0,流动是无旋流动。
- 磁场 B:∇×B=μ0J+μ0ϵ0∂t∂E (安培-麦克斯韦定律的微分形式),其中 J 是电流密度。电流和变化的电场是磁场的涡旋源。
- 静电场 Estatic:∇×Estatic=0,表明静电场是无旋场,因此可以表示为电势的梯度 Estatic=−∇V。
总结:
- 梯度 ∇f:标量场 → 矢量场。表示标量场变化最快的方向和速率。
- 散度 ∇⋅F:矢量场 → 标量场。表示矢量场的源强度或流发程度。
- 旋度 ∇×F:矢量场 → 矢量场。表示矢量场的旋转或环流程度。
这些算子是电磁学、流体力学、热力学等许多物理学分支中描述场的基本工具。