Answer

问题及解答

由多元函数的分部积分导出散度定理、Green第一恒等式、Green第二恒等式

Posted by haifeng on 2012-03-24 23:25:11 last update 2021-01-08 14:52:17 | Edit | Answers (1)

假设 $\Omega$ 是 $\mathbb{R}^n$ 中的一个有界开子集, 其边界 $\Gamma$ 是分段光滑的. 设 $u,v\in C^1(\bar{\Omega})$. 则有下面的分部积分公式.

\[\int_\Omega\frac{\partial u}{\partial x_i}vd\Omega=\int_\Gamma uv\nu_i d\Gamma-\int_\Omega u\frac{\partial v}{\partial x_i}d\Omega,\]

其中 $\hat{\nu}=(\nu_1,\nu_2,\ldots,\nu_n)$ 是 $\Gamma$ 的外法向量.

更一般的, 对于向量值函数 $\vec{V}=(v_1,v_2,\ldots,v_n)$, 有

\[\int_\Omega\nabla u\cdot\vec{V}d\Omega=\int_\Gamma(u\vec{V})\cdot\hat{\nu}d\Gamma-\int_\Omega u\nabla\cdot\vec{V}d\Omega.\]


注: (1) 令 $u\equiv 1$, 则得到 散度定理.

\[\int_\Omega\nabla\cdot\vec{V}d\Omega=\int_\Gamma\vec{V}\cdot\hat{\nu}d\Gamma,\]

\[\int_\Omega\text{div}(\vec{V})d\Omega=\int_\Gamma \vec{V}\cdot\hat{\nu}d\Gamma.\]

这里 $\text{div}(\vec{V})=\nabla\cdot\vec{V}$. 黎曼几何中流形 $M$ 上向量场 $X$ 的散度定义为

\[\text{div}X=\text{trace}(\xi\mapsto\nabla_{\xi}X)\]

(2) 若令 $\vec{V}$ 是一梯度场, 即存在函数 $v\in C^2(\bar{\Omega})$, 使得 $\vec{V}=\nabla v$, 则有

\[\int_\Omega \nabla u\cdot\nabla v d\Omega=\int_\Gamma(u\nabla v)\cdot\hat{\nu}d\Gamma-\int_\Omega u\Delta v d\Omega.\]

此即 Green 第一恒等式.


注:

(1) 有时也记 $\hat{\nu}=\vec{n}$. 或有时干脆写成 $n$. (这里 $n$ 不是指正整数, 而是 $\Gamma$ 的外法向量. 怀疑将 nu 简写为 n 了.)

(2) 由于 $(u\nabla v)\cdot\hat{\nu}=u(\nabla v\cdot\hat{\nu})=u\dfrac{\partial v}{\partial\hat{\nu}}=u\dfrac{\partial v}{\partial\vec{n}}$, 以及 $\Gamma=\partial\Omega$, 故 Green 第一恒等式也写为

\[\int_\Omega \nabla u\cdot\nabla v\mathrm{d}\Omega=\int_{\partial\Omega}u\frac{\partial v}{\partial\vec{n}}\mathrm{d}\partial\Omega-\int_\Omega u\Delta v\mathrm{d}\Omega.\qquad(*)\]

由于 $\Gamma=\partial\Omega$ 上的面积元写成 $\mathrm{d}\partial\Omega$ 不好看, 因此记之为 $\mathrm{d}S$. 既然面积元写为 $\mathrm{d}S$, 体积元为何不写为 $\mathrm{d}V$. 如此, 得到下面常见的形式:

\[\int_\Omega \nabla u\cdot\nabla v\mathrm{d}V=\int_{\partial\Omega}u\frac{\partial v}{\partial n}\mathrm{d}S-\int_\Omega u\Delta v\mathrm{d}V.\]

(3) 但是有些地方(书或其他资料)将 $\Delta v$ 写成 $\nabla^2 v$, (事实上应写为 $\Delta v=\mathrm{tr}\nabla^2 v$,  $\mathrm{tr}$ 表示对矩阵 $\nabla^2 v=\biggl(\dfrac{\partial^2 v}{\partial x_j\partial x_i}\biggr)_{n\times n}$ 取迹, 即对角线之和. 当然这里 $\nabla^2 v$ 也是指 $\sum\limits_{i=1}^{n}\dfrac{\partial^2 v}{\partial x_i^2}$.) 并将 $\Omega$ 上的积分写到一起, 从而又写为

\[\int_\Omega\Bigl(\nabla u\cdot\nabla v+u\nabla^2 v\Bigr)\mathrm{d}V=\int_{\partial\Omega}u\frac{\partial v}{\partial n}\mathrm{d}S.\qquad(1)\]

(4) 有些书或资料将 $\dfrac{\partial v}{\partial\vec{n}}$ 写为 $D_{\vec{n}}v$, 不要奇怪. 这是沿着 $\vec{n}$ 方向的方向导数. 方向导数在黎曼几何中被推广为联络, 也使用 $D$ 这个记号, 只不过下标是向量场.

(5) 在Green 第一恒等式中, 若令 $u=1$, (或等价地, 在散度定理中, 若令 $\vec{V}=\nabla v$) 则得

\[
\int_{\Omega}0\cdot\nabla v\mathrm{d}\Omega=\int_{\Gamma}(\nabla v)\cdot\hat{\nu}\mathrm{d}\Gamma-\int_{\partial\Omega}\Delta v\mathrm{d}\Omega,
\]

\[
\begin{split}
&\int_{\partial\Omega}\Delta v\mathrm{d}\Omega=\int_{\Gamma}(\nabla v)\cdot\hat{\nu}\mathrm{d}\Gamma\\
\Leftrightarrow\ &\int_{\partial\Omega}\Delta v\mathrm{d}V=\int_{\partial\Omega}\frac{\partial v}{\partial\vec{n}}\mathrm{d}S.
\end{split}
\]


 

Green 第二恒等式

根据 (1), 交换 $u$ 和 $v$ 的角色, 得到

\[\int_\Omega\Bigl(\nabla v\cdot\nabla u+v\nabla^2 u\Bigr)\mathrm{d}V=\int_{\partial\Omega}v\frac{\partial u}{\partial n}\mathrm{d}S.\qquad(2)\]

两式相减, 得

\[
\int_{\Omega}\Bigl(u\nabla^2 v-v\nabla^2 u\Bigr)\mathrm{d}V=\int_{\partial\Omega}\Bigl(u\frac{\partial v}{\partial n}-v\frac{\partial u}{\partial n}\Bigr)\mathrm{d}S.
\]

当然, 我们更喜欢写成

\[
\int_{\Omega}\Bigl(u\Delta v-v\Delta u\Bigr)\mathrm{d}V=\int_{\partial\Omega}\Bigl(uD_{\vec{n}}v-vD_{\vec{n}}u\Bigr)\mathrm{d}S.
\]

 

注: 由散度定理也可以直接推出 Green 第二恒等式. 只需在散度定理中令 $\vec{V}=u\nabla v-v\nabla u$, 计算即可推出.

 


References

http://en.wikipedia.org/wiki/Integration_by_parts

https://newikis.com/ja/グリーンの恒等式

https://www.math.ust.hk/~maklchan/ma4052/w12.pdf

1

Posted by haifeng on 2021-01-08 10:39:19

在散度定理

\[
\int_{\Omega}\mathrm{div}(\vec{V})\mathrm{d}\Omega=\int_{\Gamma}\vec{V}\cdot\hat{\nu}\mathrm{d}\Gamma
\]

中, 令 $\vec{V}=u\nabla v-v\nabla u$, 得

\[
\int_{\Omega}\mathrm{div}(\vec{V})\mathrm{d}\Omega=\int_{\Omega}\mathrm{div}(u\nabla v-v\nabla u)\mathrm{d}\Omega
=\int_{\Gamma}(u\nabla v-v\nabla u)\cdot\hat{\nu}d\Gamma
=\int_{\Gamma}(uD_{\hat{\nu}}v-vD_{\hat{\nu}}u)d\Gamma.
\]

若记 $D_i=\partial_{x_i}$, 则

\[
\begin{split}
\mathrm{div}\Bigl(u\nabla v-v\nabla u\Bigr)&=\mathrm{div}\bigl(u\cdot(D_1 v,\ldots, D_n v)\bigr)-\mathrm{div}\bigl(v\cdot(D_1 u,\ldots, D_n u)\bigr)\\
&=\Bigl[D_1(uD_1 v)+\cdots+D_n(uD_n v)\Bigr]-\Bigl[D_1(vD_1 u)+\cdots+D_n(vD_n u)\Bigr]\\
&=\Bigl[(D_1 u D_1v+\cdots D_n u D_n v)+(uD_1^2 v+\cdots uD_n^2 v)\Bigr]-\Bigl[(D_1 v D_1u+\cdots D_n v D_n u)+(vD_1^2 u+\cdots vD_n^2 u)\Bigr]\\
&=(uD_1^2 v+\cdots uD_n^2 v)-(vD_1^2 u+\cdots vD_n^2 u)\\
&=u\Delta v-v\Delta u,
\end{split}
\]

因此

\[
\int_{\Omega}(u\Delta v-v\Delta u)\mathrm{d}\Omega=\int_{\Gamma}(uD_{\hat{\nu}}v-vD_{\hat{\nu}}u)\mathrm{d}\Gamma.
\]

换个写法, 即

\[
\int_{\Omega}(u\Delta v-v\Delta u)\mathrm{d}V=\int_{\partial\Omega}(uD_{\vec{n}}v-vD_{\vec{n}}u)\mathrm{d}S.
\]