本文仅供学习使用
本文参考：
《变分法基础-第三版》老大中
《变分学讲义》张恭庆
《Calculus of Variations of Optimal Control Theory》-变分法和最优控制论-Daneil Liberzon

---

1.3.2 向量场的通量和散度

设有向量场$\vec{a}$,在场内取一有向曲面$S$,$dS$是$S$上的曲面元素，在$dS$任取一点$M$,$\vec{n}$是点$M$处外法向方向的单位向量，若将$\vec{a}\cdot \vec{n}\mathrm{d}S$对曲面$S$积分：
$$\Phi ={\iint }_S\vec{a}\cdot \mathrm{d}\vec{S}={\iint }_S\vec{a}\cdot \vec{n}\mathrm{d}S={\iint }_S{a}_n\mathrm{d}S$$
则$\phi$称为向量$\vec{a}$沿外法线向量$\vec{n}$的方向通过曲面$S$的通量。

曲面$S$在三个坐标面上的投影分别为：
$$\{\begin{array}{l}\mathrm{d}y\mathrm{d}z=\cos \left(\vec{N},x\right)\mathrm{d}\vec{S}=\cos \alpha \mathrm{d}\vec{S}=l\mathrm{d}\vec{S}\\ \mathrm{d}z\mathrm{d}x=\cos \left(\vec{N},y\right)\mathrm{d}\vec{S}=\cos \beta \mathrm{d}\vec{S}=m\mathrm{d}\vec{S}\\ \mathrm{d}x\mathrm{d}y=\cos \left(\vec{N},z\right)\mathrm{d}\vec{S}=\cos \gamma \mathrm{d}\vec{S}=n\mathrm{d}\vec{S}\end{array}$$

有向曲面$\mathrm{d}\vec{S}$可表示为：$\mathrm{d}\vec{S}=\vec{n}\mathrm{d}S=\vec{i}\mathrm{d}y\mathrm{d}z+\vec{j}\mathrm{d}z\mathrm{d}x+\vec{k}\mathrm{d}x\mathrm{d}y$

向量$\vec{a}$通过曲面$S$的通量还可以写成以下形式：
$$\begin{gathered} \Phi =\iint \left[a_x\cos \left(\vec{N},x\right)+a_{\mathrm{y}}\cos \left(\vec{N},y\right)+a_{\mathrm{z}}\cos \left(\vec{N},z\right)\right]\mathrm{d}\vec{S} \\ =\iint \left(a_x\mathrm{d}y\mathrm{d}z+a_{\mathrm{y}}\mathrm{d}z\mathrm{d}x+a_{\mathrm{z}}\mathrm{d}x\mathrm{d}y\right) \end{gathered}$$
当曲面$S$是封闭曲面时，具有特别重要的意义。为了明显起见，采用在积分号上加圆圈的方法，即用：
$$\Phi =\iint a_{\mathrm{n}}\mathrm{d}S=\iint \vec{a}\mathrm{d}\vec{S}$$
来表示通过封闭曲面$S$的通量。

设有向量场$\vec{a}=\vec{a}\left(M\right)$ ,在场中一点$M$的某个邻域内作包含此点的任一封闭曲面$\Delta S$,$\Delta S$所围成的体积为$\Delta V$,计算向量$\vec{a}$穿过$\Delta S$面上的通量，当$\Delta V$以任意方式向点$M$收缩时，如果比式：
$$\underset{\Delta V\to 0}{\lim }\frac{\Delta \Phi }{\Delta V}=\underset{\Delta V\to 0}{\lim }\frac{\iint \vec{a}\cdot \mathrm{d}\vec{S}}{\Delta V}$$
的极限存在，则此极限称为向量函数$\vec{a}$在点$M$的散度，记作$\mathrm{div}\vec{a}$,$\mathrm{div}$为英文divergence的缩写，意为散度。当采用哈密顿算子时，可写作$\nabla \cdot \vec{a}$,即:
$$\begin{gathered} \mathrm{div}\vec{a}=\nabla \cdot \vec{a}=\underset{\Delta V\to 0}{\lim }\frac{\Delta \Phi }{\Delta V}=\underset{\Delta V\to 0}{\lim }\frac{\iint \vec{a}\cdot \mathrm{d}\vec{S}}{\Delta V} \\ \nabla \cdot \vec{a} \end{gathered}$$
显然点M的散度与所取的封闭曲面$\Delta S$的形状无关，它是一个不依赖于坐标系选取的数量，因此是一个标量。散度是通量对体积的变化率。散度为零的向量场称为无源场或无散场，也称为管状场。

虽然向量场的散度是一个不依赖于坐标系选取的标量，但在具体计算向量场的散度时，常需要用向量的分量来表示，而向量的分量和坐标的选择有关。因此需要讨论向量的散度与其分量的关系。

设直角坐标系向量：$\vec{a}=a_x\vec{i}+a_{\mathrm{y}}\vec{j}+a_{\mathrm{z}}\vec{k}$,微元体积$\Delta V$是以$\Delta x$、$\Delta y$和$\Delta z$为边长的微元六面体，并设微元六面体的中心坐标为$\left(x,y,z\right)$。作为一阶近似，在$x$方向存在如下关系：
$${a_x|}_{x+\frac{\Delta x}{2}}={a_x|}_x+\frac{\partial a_x}{\partial x}\frac{\Delta x}{2},{a_x|}_{x-\frac{\Delta x}{2}}={a_x|}_x-\frac{\partial a_x}{\partial x}\frac{\Delta x}{2}$$
在$x$轴方向流出微元六面体的净通量为：
$$\left({a_x|}_{x+\frac{\Delta x}{2}}-{a_x|}_{x-\frac{\Delta x}{2}}\right)\Delta y\Delta z=\frac{\partial a_x}{\partial x}\Delta x\Delta y\Delta z$$
同样，在$y$轴和$z$轴方向也可以到类似的关系，于是流出微元体积$\Delta V$的净通量为：
$$\iint \vec{a}\cdot \mathrm{d}\vec{S}=\frac{\partial a_x}{\partial x}\Delta x\Delta y\Delta z+\frac{\partial a_{\mathrm{y}}}{\partial y}\Delta x\Delta y\Delta z+\frac{\partial a_{\mathrm{z}}}{\partial z}\Delta x\Delta y\Delta z$$
将上式两端除以微元体积$\Delta V=\Delta x\Delta y\Delta z$并取极限，得：
$$\mathrm{div}\vec{a}=\underset{\Delta V\to 0}{\lim }\frac{\iint \vec{a}\cdot \mathrm{d}\vec{S}}{\Delta V}=\frac{\partial a_x}{\partial x}+\frac{\partial a_{\mathrm{y}}}{\partial y}+\frac{\partial a_{\mathrm{z}}}{\partial z}$$
又：
$$\nabla \cdot \vec{a}=\left(\frac{\partial }{\partial x}\vec{i}+\frac{\partial }{\partial y}\vec{j}+\frac{\partial }{\partial z}\vec{k}\right)\cdot \left(a_x\vec{i}+a_{\mathrm{y}}\vec{j}+a_{\mathrm{z}}\vec{k}\right)=\frac{\partial a_x}{\partial x}+\frac{\partial a_{\mathrm{y}}}{\partial y}+\frac{\partial a_{\mathrm{z}}}{\partial z}$$
故：
$$\mathrm{div}\vec{a}=\nabla \cdot \vec{a}=\frac{\partial a_x}{\partial x}+\frac{\partial a_{\mathrm{y}}}{\partial y}+\frac{\partial a_{\mathrm{z}}}{\partial z}$$
上式就是散度在直角坐标系中的表达式。在计算向量场的散度时，应用该式通常要比直接引用定义式更方便。一个向量函数a的散度是一个标量函数。在向量场中任一点，向量$\vec{a}$的散度等于它在各坐标轴上的分量对各自坐标变量的偏导数之和。

散度具有下列性质：

• 两个向量函数之和的散度等于各向量函数的散度之和
  $$\begin{gathered} \mathrm{div}\left(\vec{a}+\vec{b}\right)=\mathrm{div}\vec{a}+\mathrm{div}\vec{b} \\ \nabla \cdot \left(\vec{a}+\vec{b}\right)=\nabla \cdot \vec{a}+\nabla \cdot \vec{b} \end{gathered}$$
• 数量函数与向量函数之积的散度等于数量函数与向量函数的散度之积加上数量函数的梯度与向量函数之内积
  $$\begin{gathered} \mathrm{div}\left(\phi \vec{a}\right)=\phi \mathrm{div}\vec{a}+\vec{a}\cdot \mathrm{grad}\phi \\ \nabla \cdot \left(\phi \vec{a}\right)=\phi \nabla \cdot \vec{a}+\vec{a}\cdot \nabla \phi \end{gathered}$$

证明：（待补充）

求梯度的散度：$\nabla \cdot \nabla \left(\phi \psi \right)$
$$\nabla \cdot \nabla \left(\phi \psi \right)=\nabla \cdot \left(\phi \nabla \psi +\psi \nabla \phi \right)=\phi \nabla \cdot \nabla \psi +\nabla \phi \cdot \nabla \psi +\psi \nabla \cdot \nabla \phi +\nabla \psi \cdot \nabla \phi =\phi {\nabla }^2\psi +2\nabla \phi \cdot \nabla \psi +\psi {\nabla }^2\phi =\phi \Delta \psi +2\nabla \phi \cdot \nabla \psi +\psi \Delta \phi$$
式中，算子$\Delta =\nabla \cdot \nabla ={\nabla }^2$,作用在标量函数上的算子$\Delta$ 称为标量拉普拉斯算子，简称拉普拉斯算子或调和算子。在直角坐标系中，它可表示为：
$$\Delta =\nabla \cdot \nabla ={\nabla }^2=\frac{{\partial }^2}{\partial x^2}+\frac{{\partial }^2}{\partial y^2}+\frac{{\partial }^2}{\partial z^2}$$

1.3.3 高斯定理与格林公式

定理： 设向量$\vec{a}$在一封闭曲面$S$的内部有连续一阶偏导数，则$\vec{a}$在$S$上的通量等于该向量的散度对$S$所围成体积$V$上的积分，即：
$$\iint \vec{a}\cdot \mathrm{d}\vec{S}=\iint \vec{a}\cdot \vec{n}\mathrm{d}S={\iiint }_V\mathrm{div}\vec{a}dV={\iiint }_V\nabla \cdot \vec{a}dV$$
式中，$\vec{n}$为$S$上单位外法线向量；$\nabla$ 为哈密顿算子。

上式称为高斯-奥斯特罗格拉茨基定理，简称高斯定理，也称为散度形式的高斯公式或散度定理，它建立了空间区域V上连续函数的三重积分与其边界面S上的曲面积分之间的关系。
证明：（待补充）

• 高斯公式在理论研究和实际工作中有着广泛的应用，是一种有力的数学工具，许多公式都是从高斯公式出发而得到的。

设：
$$\begin{gathered} \vec{a}=P\left(x,y,z\right)\vec{i}+Q\left(x,y,z\right)\vec{j}+R\left(x,y,z\right)\vec{k} \\ \vec{n}=\cos \left(n,x\right)\vec{i}+\cos \left(n,y\right)\vec{j}+\cos \left(n,z\right)\vec{k}=\cos \alpha \vec{i}+\cos \beta \vec{j}+\cos \gamma \vec{k} \end{gathered}$$

将以上两式带入高斯公式，可得：
$$\iint \vec{a}\cdot \mathrm{d}\vec{S}=\iint \vec{a}\cdot \vec{n}\mathrm{d}S=\iint \left(P\cos \alpha +Q\cos \beta +R\cos \gamma \right)\mathrm{d}S={\iiint }_V\left(\frac{\partial P}{\partial x}+\frac{\partial Q}{\partial y}+\frac{\partial R}{\partial z}\right)dV$$
这个高斯公式是常见的一种形式。

设有平面向量场$\vec{a}$,在场内取一有向曲线$\Gamma$ ,$\mathrm{d}\Gamma$ 是$\Gamma$上的曲线元素，在$\mathrm{d}\Gamma$任取一点$M$,$\vec{n}$是点$M$处外法线方向的单位向量，若将$\vec{a}\cdot \vec{n}\mathrm{d}\Gamma$对于曲线$\mathrm{d}\Gamma$积分：
$$\Phi ={\int }_{\Gamma }\vec{a}\cdot \vec{n}\mathrm{d}\Gamma ={\int }_{\Gamma }\vec{a}\mathrm{d}\vec{\Gamma }$$
则称为$\Phi$ 向量场$\vec{a}$沿法向向量$\vec{n}$的方向通过曲线$\Gamma$的通量。在直角坐标系中，曲线$\Gamma$的法向单位向量$\vec{n}$、切向单位向量$\vec{\tau }$与该曲线存在下列关系：
$$\{\begin{array}{l}\vec{n}={\vec{n}}_x+{\vec{n}}_{\mathrm{y}}=n_x\vec{i}+n_{\mathrm{y}}\vec{j}=\cos \left(\vec{n},x\right)\vec{i}+\cos \left(\vec{n},y\right)\vec{j}=\cos \left(\vec{\tau },y\right)\vec{i}+\cos \left(\vec{\tau },-x\right)\vec{j}=\frac{\mathrm{d}y}{\mathrm{d}\Gamma }\vec{i}-\frac{\mathrm{d}x}{\mathrm{d}\Gamma }\vec{j}\\ \vec{\tau }={\vec{\tau }}_x+{\vec{\tau }}_{\mathrm{y}}={\tau }_x\vec{i}+{\tau }_{\mathrm{y}}\vec{j}=\cos \left(\vec{\tau },x\right)\vec{i}+\cos \left(\vec{\tau },y\right)\vec{j}=\cos \left(\vec{n},-y\right)\vec{i}+\cos \left(\vec{n},x\right)\vec{j}=-\frac{\mathrm{d}x}{\mathrm{d}\Gamma }\vec{i}+\frac{\mathrm{d}y}{\mathrm{d}\Gamma }\vec{j}\end{array}$$

设$\vec{a}=P\left(x,y\right)\vec{i}+Q\left(x,y\right)\vec{j}$，则通过曲线$\Gamma$的通量$\Phi$可写成：
$$\Phi ={\int }_{\Gamma }\vec{a}\cdot \vec{n}\mathrm{d}\Gamma ={\int }_{\Gamma }\vec{a}\mathrm{d}\vec{\Gamma }={\int }_{\Gamma }\left(P\mathrm{d}y-Q\mathrm{d}x\right)$$

设有平面向量场$\vec{a}=\vec{a}\left(M\right)$ ,在场中一点$M$的某个邻域内作包含此点的任一封闭曲线$\Delta \Gamma$ ,$\Delta \Gamma$所围成的面积为$\Delta D$,计算向量$\vec{a}$经过$\Delta \Gamma$线上的通量，当$\Delta D$以任意方式向点$M$收缩时，如果比式：
$$\underset{\Delta D\to 0}{\lim }\frac{\Delta \Phi }{\Delta D}=\underset{\Delta D\to 0}{\lim }\frac{{\oint }_{\Delta \Gamma }\vec{a}\cdot \vec{n}\mathrm{d}\Gamma }{\Delta D}$$
的极限存在，则此极限称为向量函数$\vec{a}$在点$M$的散度，记作$\mathrm{div}\vec{a}$,当采用哈密顿算子时，可写作$\nabla \cdot \vec{a}$，即：
$$\mathrm{div}\vec{a}=\nabla \cdot \vec{a}=\underset{\Delta D\to 0}{\lim }\frac{\Delta \Phi }{\Delta D}=\underset{\Delta D\to 0}{\lim }\frac{{\oint }_{\Delta \Gamma }\vec{a}\cdot \vec{n}\mathrm{d}\Gamma }{\Delta D}$$

显然点$M$的散度与所取的封闭曲线$\Delta \Gamma$的形状无关，它是一个不依赖于坐标系选取的数量，因此是一个标量。
令$\vec{b}=\varphi \vec{a}$,利用散度的性质,可得到以下公式：
$$\begin{gathered} {\iiint }_V\mathrm{div}\left(\varphi \vec{a}\right)\mathrm{d}V={\iiint }_V\nabla \cdot \left(\varphi \vec{a}\right)\mathrm{d}V={\iiint }_V\left(\varphi \nabla \cdot \vec{a}+\vec{a}\nabla \varphi \right)\mathrm{d}V \\ =\iint \varphi \vec{a}\cdot \mathrm{d}\vec{S}=\iint \varphi \left(\vec{a}\cdot \vec{n}\right)\mathrm{d}S \end{gathered}$$

定理： 设平面向量$\vec{a}$在一封闭曲线$\Gamma$ 的内部有连续一阶偏导数，则$\vec{a}$在$\Gamma$上的通量等于该向量的散度对$\Gamma$所围成面积$D$上的积分，即：
$${\oint }_{\Gamma }\vec{a}\cdot \vec{n}\mathrm{d}\Gamma ={\oint }_{\Gamma }{a}_n\mathrm{d}\Gamma ={\iint }_D\mathrm{div}\vec{a}\mathrm{d}D={\iint }_D\nabla \cdot \vec{a}\mathrm{d}D$$
式中，$\vec{n}$为$\Gamma$上外法线方向的单位向量；$\nabla$为哈密顿算子。上式称为格林公式或格林定理，它建立了平面区域$D$上的二重积分和沿该区域边界$\Gamma$的曲线积分之间的联系。

• 证明的方法与高斯定理的证明方法相同，只需将封闭曲面$S$换成封闭曲线$\Gamma$,并将$S$围成的体积$V$换成$\Gamma$围成的面积$D$即可。

设$\vec{a}=Q\left(x,y\right)\vec{i}-P\left(x,y\right)\vec{j}$,并注意到$\vec{n}=\frac{\mathrm{d}y}{\mathrm{d}\Gamma }\vec{i}-\frac{\mathrm{d}x}{\mathrm{d}\Gamma }\vec{j}$，则格林公式可写成：
$${\oint }_{\Gamma }\left(P\left(x,y\right)\mathrm{d}x+Q\left(x,y\right)\mathrm{d}y\right)={\iint }_D\left(\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\mathrm{d}x\mathrm{d}y$$

这里曲线积分沿着曲线$\Gamma$的正向，二重积分展布在区域$D$上。所谓曲线$\Gamma$的正向是这样的方向，设想一个人绕着曲线$\Gamma$行走，使得$\Gamma$围成的区域$D$总是保持在其左侧。

令$P=-y,Q=x$，则有：$D={\iint }_D\mathrm{d}x\mathrm{d}y=\frac{1}{2}{\oint }_{\Gamma }\left(x\mathrm{d}y-y\mathrm{d}x\right)$
上式为利用曲线积分计算平面区域$D$的面积公式

若$x=x\left(t\right),y=y\left(t\right)$，则上式可写成：$D={\iint }_D\mathrm{d}x\mathrm{d}y=\frac{1}{2}{\oint }_{\Gamma }\left(x\dot{y}-y\dot{x}\right)\mathrm{d}t$

令$P=u\frac{\partial u}{\partial y},Q=-u\frac{\partial u}{\partial x}$，则有：
$${\iint }_D\left(\frac{{\partial }^{2}u}{\partial x^2}+\frac{{\partial }^{2}u}{\partial y^2}\right)u\mathrm{d}x\mathrm{d}y+{\iint }_D\left[{\left(\frac{\partial u}{\partial x}\right)}^2+{\left(\frac{\partial u}{\partial y}\right)}^2\right]\mathrm{d}x\mathrm{d}y=-{\oint }_{\Gamma }\left(u\frac{\partial u}{\partial y}\mathrm{d}x-u\frac{\partial u}{\partial x}\mathrm{d}y\right)$$

设函数$\phi \text{、}\psi$ 和$\lambda$ 在区域$V$内有连续二阶导数，令$\vec{a}=\phi \lambda \nabla \psi$ ,根据散度的性质，有下面公式：
$$\nabla \cdot \vec{a}=\nabla \cdot \left(\phi \lambda \nabla \psi \right)=\phi \nabla \cdot \left(\lambda \nabla \psi \right)+\lambda \nabla \phi \cdot \nabla \psi$$
进而可知：
$$\vec{a}\cdot \vec{n}=\phi \lambda \nabla \psi \cdot \vec{n}=\phi \lambda \left|\nabla \psi \right|\vec{n}\cdot \vec{n}=\phi \lambda \frac{\partial \psi }{\partial n}$$
$$\begin{gathered} {\iiint }_V\nabla \cdot \left(\phi \lambda \nabla \psi \right)\mathrm{d}V={\iiint }_V\left[\phi \nabla \cdot \left(\lambda \nabla \psi \right)+\lambda \nabla \phi \cdot \nabla \psi \right]\mathrm{d}V \\ =\iint \phi \lambda \frac{\partial \psi }{\partial n}\mathrm{d}S=\iint \phi \lambda \left|\nabla \psi \right|\mathrm{d}S=\iint \phi \lambda \nabla \psi \cdot \vec{n}\mathrm{d}S \end{gathered}$$
若令上式中的$\lambda =1$,则可得到格林第一公式或奥斯特罗格拉茨基公式，也称为
格林第一定理或奥氏公式，即：
$$\begin{gathered} {\iiint }_V\nabla \cdot \left(\phi \nabla \psi \right)\mathrm{d}V={\iiint }_V\left(\phi \nabla \cdot \nabla \psi +\nabla \phi \cdot \nabla \psi \right)\mathrm{d}V={\iiint }_V\left(\phi \Delta \psi +\nabla \phi \cdot \nabla \psi \right)\mathrm{d}V \\ =\iint \phi \frac{\partial \psi }{\partial n}\mathrm{d}S=\iint \phi \lambda \left|\nabla \psi \right|\mathrm{d}S=\iint \phi \vec{n}\cdot \nabla \psi \mathrm{d}S \end{gathered}$$
将上式中的$\phi$与$\psi$易位，得：
$$\begin{gathered} {\iiint }_V\nabla \cdot \left(\psi \nabla \phi \right)\mathrm{d}V={\iiint }_V\left(\psi \Delta \phi +\nabla \psi \cdot \nabla \phi \right)\mathrm{d}V=\iint \psi \frac{\partial \phi }{\partial n}\mathrm{d}S \\ =\iint \psi \left|\nabla \phi \right|\mathrm{d}S=\iint \psi \nabla \phi \cdot \vec{n}\mathrm{d}S \end{gathered}$$
也是格林第一公式的另一种形式。将格林第一公式两种不同形式相减，便得到格林第二公式，也称为格林第二定理，即：
$${\iiint }_V\left(\phi \Delta \psi -\psi \Delta \phi \right)\mathrm{d}V=\iint \left(\phi \frac{\partial \psi }{\partial n}-\psi \frac{\partial \phi }{\partial n}\right)\mathrm{d}S=\iint \left(\phi \nabla \psi -\psi \nabla \phi \right)\cdot \vec{n}\mathrm{d}S$$
当$\phi =\psi$时，即可得到格林第三公式，也称为格林第三定理，即：
$${\iiint }_V\left[\phi \Delta \phi +{\left(\nabla \phi \right)}^2\right]\mathrm{d}V=\iint \phi \frac{\partial \psi }{\partial n}\mathrm{d}S=\iint \phi \nabla \phi \cdot \vec{n}\mathrm{d}S$$
令$\psi =1$，则：
$${\iiint }_V\Delta \phi \mathrm{d}V=\iint \phi \frac{\partial \phi }{\partial n}\mathrm{d}S=\iint \phi \nabla \phi \cdot \vec{n}\mathrm{d}S$$
在直角坐标系中，上式可写成：
$${\iiint }_V\left(\frac{{\partial }^2\phi }{\partial x^2}+\frac{{\partial }^2\phi }{\partial y^2}+\frac{{\partial }^2\phi }{\partial z^2}\right)\mathrm{d}V=\iint \left(\frac{\partial \phi }{\partial x}\mathrm{d}y\mathrm{d}z+\frac{\partial \phi }{\partial y}\mathrm{d}z\mathrm{d}x+\frac{\partial \phi }{\partial z}\mathrm{d}x\mathrm{d}y\right)$$
若$\Delta \phi =0$，则：
$$\begin{gathered} {\iiint }_V{\left(\Delta \phi \right)}^2\mathrm{d}V=\iint \phi \frac{\partial \phi }{\partial n}\mathrm{d}S=\iint \phi \nabla \phi \cdot \vec{n}\mathrm{d}S \\ \iint \frac{\partial \phi }{\partial n}\mathrm{d}S=\iint \nabla \phi \cdot \vec{n}\mathrm{d}S \end{gathered}$$
同理，只需将对封闭曲面$S$的积分换成对封闭曲线$\Gamma$ 的积分，将对$S$围成的体积$V$的积分换成对$\Gamma$围成的面积$D$的积分，并将三维的哈密顿算子和拉普拉斯算子换成二维的即可。

从历史上看，格林公式是独立提出来的，因而是一个原始公式。但由于它和高斯公式密切相关，故可认为高斯公式是平面格林公式在空间的推广，而格林公式是高斯公式在平面上的特殊情况。

• 若将格林第一公式中的$\phi$换成$\Delta \phi$，则有：
  $$\begin{gathered} {\iiint }_V\nabla \cdot \left(\Delta \phi \nabla \psi \right)\mathrm{d}V={\iiint }_V\left(\Delta \phi \nabla \cdot \nabla \psi +\nabla \Delta \phi \cdot \nabla \psi \right)\mathrm{d}V={\iiint }_V\left(\Delta \phi \Delta \psi +\nabla \Delta \phi \cdot \nabla \psi \right)\mathrm{d}V \\ =\iint \Delta \phi \frac{\partial \psi }{\partial n}\mathrm{d}S=\iint \Delta \phi n\cdot \nabla \psi \mathrm{d}S \end{gathered}$$
• 若将格林第一公式中的$\psi$换成$\Delta \psi$，则有：
  $$\begin{gathered} {\iiint }_V\nabla \cdot \left(\phi \nabla \Delta \psi \right)\mathrm{d}V={\iiint }_V\left(\phi \nabla \cdot \nabla \Delta \psi +\nabla \phi \cdot \nabla \Delta \psi \right)\mathrm{d}V={\iiint }_V\left(\phi {\Delta }^2\psi +\nabla \phi \cdot \nabla \Delta \psi \right)\mathrm{d}V \\ =\iint \phi \frac{\partial \Delta \psi }{\partial n}\mathrm{d}S=\iint \phi n\cdot \nabla \Delta \psi \mathrm{d}S \end{gathered}$$
  将上式中的$\phi$与$\psi$易位，得：
  $$\begin{gathered} {\iiint }_V\nabla \cdot \left(\psi \nabla \Delta \phi \right)\mathrm{d}V={\iiint }_V\nabla \left(\psi \nabla \cdot \nabla \Delta \phi +\nabla \psi \cdot \nabla \Delta \phi \right)\mathrm{d}V={\iiint }_V\nabla \left(\psi {\Delta }^2\phi +\nabla \psi \cdot \nabla \Delta \phi \right)\mathrm{d}V \\ =\iint \psi \frac{\partial \Delta \phi }{\partial n}\mathrm{d}S=\iint \psi n\cdot \nabla \Delta \phi \end{gathered}$$

进一步整理则有：
$$\begin{gathered} {\iiint }_V\psi {\Delta }^2\phi \mathrm{d}V=-{\iiint }_V\nabla \psi \cdot \nabla \Delta \phi \mathrm{d}V+\iint \psi \frac{\partial \Delta \phi }{\partial n}\mathrm{d}S \\ ={\iiint }_V\Delta \phi \Delta \psi \mathrm{d}V+\iint \left(\psi \frac{\partial \Delta \phi }{\partial n}-\Delta \phi \frac{\partial \psi }{\partial n}\right)\mathrm{d}S \end{gathered}$$

若标量函数$\phi$在封闭曲面$S$内具有连续的一阶偏导数，且$S$所围成的体积为$V$，则有下列梯度公式：
$$\iint \phi \mathrm{d}\vec{S}=\iint \phi \vec{n}\mathrm{d}S=\iint \vec{n}\phi \mathrm{d}S={\iiint }_V\mathrm{grad}\phi \mathrm{d}V={\iiint }_V\nabla \phi \mathrm{d}V$$
上式称为梯度形式的高斯公式或梯度定理。

证明：（待补充）

例： 设$v\Delta \left(p\Delta u\right)=u\Delta \left(p\Delta v\right)$，求证：$\nabla \cdot \left[v\nabla \left(p\Delta u\right)-u\nabla \left(p\Delta v\right)+p\Delta v\nabla u-p\Delta u\nabla v\right]=0$
令：$\phi =p\Delta u,\psi =p\Delta v$，根据散度的性质有：
$$v\Delta \left(p\Delta u\right)=v\Delta \phi =\nabla \cdot \left(v\nabla \phi \right)-\nabla v\cdot \nabla \phi$$
其中：$\nabla v\cdot \nabla \phi =\nabla \cdot \left(\phi \nabla v\right)-\phi \nabla \cdot \nabla v=\nabla \cdot \left(\phi \nabla v\right)-\phi \Delta v$，带入上式可得：$v\Delta \left(p\Delta u\right)=v\Delta \phi =\nabla \cdot \left(v\nabla \phi -\phi \nabla v\right)+\phi \Delta v$
将$u,v$易位，并由$\psi =p\Delta v$可得：$u\Delta \left(p\Delta v\right)=u\Delta \psi =\nabla \cdot \left(u\nabla \psi -\psi \nabla u\right)+\psi \Delta u$
将上述两式相减，并注意到有：$\phi \nabla v=\psi \Delta u=p\Delta v\Delta u$，得：
$$\begin{gathered} \nabla \cdot \left(v\nabla \phi -u\nabla \psi +\psi \nabla u-\phi \nabla v\right)=0 \\ \nabla \cdot \left[v\nabla \left(p\Delta u\right)-u\nabla \left(p\Delta v\right)+p\Delta v\nabla u-p\Delta u\nabla v\right]=0 \end{gathered}$$