跳转至

10.1.Describe a Vector Feild


1.1.梯度(grad)

对于数量集函数\(f(x, y)\),其\(\mathbf{l}\)方向上的方向导数为:

\[ \frac{\partial f}{\partial \mathbf{l}}=\mathbf{e_l}\cdot (\frac{\partial f}{\partial x},\frac{\partial f}{\partial y})^T \]

\((\frac{\partial f}{\partial x},\frac{\partial f}{\partial y})^T\)记作数量集函数\(f\)梯度

\[ \nabla f=(\frac{\partial f}{\partial x},\frac{\partial f}{\partial y})^T \]

有以下推论:

  1. 当单位向量\(\mathbf{e_l}\)与梯度\(\nabla f\)重合时,此方向上单位导数取得最大值,即\(\nabla f\)方向为函数值增长最快方向;
  2. \(\mathbf{e_l}\cdot\nabla f=0\)时,此方向上方向导数恒为\(0\),即\(\nabla f\)方向恒垂直于其等值线。

  3. 梯度与保守场

\(\nabla\)算子将数量集函数\(f\)映射到向量集函数\(\mathbf{F}=(P(x, y), Q(x, y))\),将向量场\(\mathbf{F}\)称为\(f\)的梯度场。

因此对于向量场\(\mathbf{F}\)分量\(P(x, y), Q(x, y)\),总能找到一个标量场\(f\),使得:

\[ P=\frac{\partial f}{\partial x}, Q=\frac{\partial f}{\partial y} \]

假设平面曲线\(L\)可利用\(x=\lambda(t), y=\mu(t)\)参数化表示,则在曲线\(C\)上:

\[ \frac{\text{d}f}{\text{d}t}=\frac{\partial P}{\partial x}\frac{\text{d}x}{\text{d}t}+\frac{\partial Q}{\partial y}\frac{\text{d}y}{\text{d}t} \]
\[ f(\mathbf{b})-f(\mathbf{a})=\int_{\mathbf{a}}^{\mathbf{b}}\mathbf{F}\text{d}\mathbf{s}=\int_L\mathbf{F}\text{d}\mathbf{s} \]

可以看出,梯度场\(\mathbf{F}\)\(A\)\(B\)的积分值与积分路径无关,只与积分的起点与终点有关,因此将\(\mathbf{F}\)称为保守场,将标量场\(f\)称为保守场\(\mathbf{F}\)

由以上性质,若对保守场\(\mathbf{F}\)中的环路\(C\)积分,则有:

\[ \oint_C \mathbf{F}\text{d}\mathbf{s}=0 \]
  • 判断保守场

若向量场\(\mathbf{F}(x, y)=(P, Q)\)是保守的,总能找到一个标量场\(f\),使得:

\[ P=\frac{\partial f}{\partial x}, Q=\frac{\partial f}{\partial y} \]

由于标量场\(f\)光滑且连续,因此:

\[ \frac{\partial P}{\partial y}=\frac{\partial ^2f}{\partial x\partial y}, \frac{\partial Q}{\partial x}=\frac{\partial ^2f}{\partial x\partial y} \]
\[ 即:\frac{\partial P}{\partial y}=\frac{\partial Q}{\partial x} \]

我们证明了其必要条件,现在证明充分条件:

定义\(\nabla\)算子:\(\nabla = (\frac{\partial}{\partial x_1}, \frac{\partial}{\partial x_2}, \cdot\cdot\cdot)\)

对于二维标量场\(f\)\(\nabla\)算子与其作点积将\(\mathbb{R}\to\mathbb{R}^2\),即:

\[ \nabla\cdot f=(\frac{\partial f}{\partial x}, \frac{\partial f}{\partial y}) \]

对于二维向量场\(\mathbf{F}\)\(\nabla\)算子与其作叉积将\(\mathbb{R}^2\to\mathbb{R}\),即:

\[ \nabla \times \mathbf{F}=\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y} \]

考虑到以下等式恒成立:

\[ \nabla \times (\nabla \cdot f)=0 \]

若对于二维向量场\(\mathbf{F}=(P, Q)\),有\(\frac{\partial Q}{\partial x}=\frac{\partial P}{\partial y}\),则:

\[ \nabla \times \mathbf{F}=0 \]

即必定存在一个标量场\(f\),使得\(\mathbf{F}=\nabla \cdot f\),即向量场\(\mathbf{F}\)是保守的,标量场\(f\)是该保守场的势。

因此\(\frac{\partial Q}{\partial x}=\frac{\partial P}{\partial y}\)\(\mathbf{F}=(P, Q)\)为保守场的充要条件。

关于\(\nabla\)算子的更多性质,会在下面的章节中介绍。


1.2.通量与散度(div)

向量集函数\(\mathbf{F}(x,y)\)表示一个平面流体的速度场,其通过曲线\(L\)通量可表示为:

\[ B=\int _L \mathbf{F}\cdot\mathbf{n}\text{d}s \]

其中\(\mathbf{n}\)垂直于曲线\(L\)一小段弧长\(s\)的切方向\(\mathbf{s}\)

当该速度场通过闭合曲线\(C\)时,规定向外为正方向,其通量可表示为:

\[ B=\oint_C\mathbf{F}\cdot \mathbf{n}\text{d}s \]

假定在平面向量场\(\mathbf{F}=(P, Q)\)中一点\((x_0, y_0)\)邻域内由\((x_0+\Delta x, y_0),(x_0, y_0+\Delta y),(x_0+\Delta x, y_0+\Delta y)\)确定的矩形\(D\)围成,通过该领域的总通量为:

\[ B=(\frac{\partial P_0}{\partial x}+\frac{\partial Q_0}{\partial y})\Delta x\Delta y \]

我们关心向量场\(\mathbf{F}\)在各区域的通量密度,因此定义向量场\(\mathbf{F}\)的散度(div)为:

\[ \text{div} \mathbf{F}=\nabla \cdot \mathbf{F}=\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y} \]

考虑矩形\(D\)的相邻区域\(D_1, D_2, D_3\)。流入相邻区域的通量等大且反向,因此对于更大的封闭区域C我们只需要关心边界的流入与流出情况。基于上述思路,我们不加证明地给出Green公式的通量形式:

\[ \oint_C\mathbf{F}\cdot \mathbf{n}\text{d}s=\iint _C \nabla \cdot \mathbf{F}\text{d}\sigma \]

1.3.旋度(curl)

向量集函数\(\mathbf{F}(x,y)\)表示一个平面流体的速度场,规定逆时针为正方向,其通过闭合曲线\(C\)旋量可表示为:

\[ \oint_C \mathbf{F}\cdot \mathbf{r}\text{d}s \]

基于上述思路,我们同样可以定义向量场\(\mathbf{F}=(P, Q)\)的旋度:

\[ \text{curl} \mathbf{F}=\nabla \times \mathbf{F}=\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y} \]

同样有Green公式的旋量形式:

\[ \oint_C\mathbf{F}\cdot \mathbf{r}\text{d}s=\iint_C \nabla \times \mathbf{F}\text{d}\sigma \]

根据亥姆霍兹分解(Helmholtz Decomposition),任何光滑且快速衰减的二维向量场\(\mathbf{F}\)(无调和)可分解为无旋的梯度场与无散的旋度场,即:

\[ \mathbf{F}=\nabla\cdot\phi+\nabla\times A \]
\[ \nabla\cdot\mathbf{F}=\nabla\cdot(\nabla\cdot \phi) \]
\[ \nabla\times\mathbf{F}=\nabla\times(\nabla\times A) \]

由此可见,任意封闭区域\(D\)的所有散度都由梯度场\(\nabla\cdot\phi\)贡献,而旋度所有都由旋度场\(\nabla\times A\)贡献。若对向量场\(\mathbf{F}\)作正交变换:

\[ \mathbf{F}^{\perp}=\nabla\times\phi+\nabla\cdot A \]

此时对于\(\phi\)\(A\),他们此时分别表示\(\mathbf{F}^{\perp}\)的旋度场与梯度场的标量势,因此对\(\mathbf{F}^{\perp}\)旋量形式的Green公式:

\[ \oint_C \mathbf{F}^{\perp}\cdot \mathbf{r}\text{d}s=\iint_C\nabla \times\mathbf{F}^{\perp}\text{d}\sigma \]

即:

\[ \oint_C\mathbf{F}\cdot\mathbf{n}\text{d}s=\iint_C\nabla\cdot \mathbf{F}\text{d}\sigma \]

因此,散度和旋度在正交变换下是对偶的,描述向量场\(\mathbf{F}\)的同一性质。