グリーンの定理
定理
$u, v \in C^2( \bar{U})$としよう。それでは、以下の式が成り立つ。
(i) $\displaystyle \int_{U} \Delta u dx=\int_{\partial U} \dfrac{\partial u}{\partial \nu}dS$
(ii) $\displaystyle \int_{U} Dv \cdot Du dx = -\int_{U} u \Delta v dx+\int_{\partial U}\dfrac{\partial v}{\partial \nu}udS$
(iii) $\displaystyle \int_{U} (u\Delta v - v\Delta u )dx = \int_{\partial U} \left( \dfrac{\partial v}{\partial \nu}u - \dfrac{\partial u}{\partial \nu} v\right)dS$
これらをまとめて グリーンの公式Green’s formulaという。
- $\Delta$はラプラシアン
- $D$はグラディエント
- $\nu$は外向き単位法線ベクトル
証明
$u, v \in C^1(\bar{U})$とする。それでは、以下の式が成り立つ。
$$ \int_{U} u_{x_{i}}vdx = -\int_{U} uv_{x_{i}}dx + \int_{\partial U} uv\nu^{i} dS\quad (i=1,\dots , n) $$
(i)
部分積分公式において、$u$の代わりに$u_{x_{i}}$を、$v$の代わりに$1$を代入すると、以下の式を得る。
$$ \int_{U} u_{x_{i} x_{i}}dx = \int_{\partial U} u_{x_{i}}\nu^{i} dS \quad (i=1,\cdots , n) $$
全ての$i=1,\cdots, n$に対して加算すると、以下のようになる。
$$ \int_{U} (u_{x_{1} x_{1}}+\cdots +u_{x_{n} x_{n}} )dx = \int_{\partial U}( u_{x_{1}}\nu^{1} +\cdots u_{x_{n}}\nu^n)dS $$
ラプラシアンの定義と$\dfrac{\partial u}{\partial \nu}:=\boldsymbol{\nu}\cdot Du$により、次が成り立つ。
$$ \int_{U} \Delta u dx=\int_{\partial U} \dfrac{\partial u}{\partial \nu}dS $$
■
(ii)
部分積分公式において、$v$の代わりに$v_{x_{i}}$を代入すると、以下の式を得る。
$$ \int_{U} u_{x_{i}}v_{x_{i}}dx = -\int_{U} uv_{x_{i}x_{i}}dx + \int_{\partial U} uv_{x_{i}}\nu^{i} dS \quad (i=1,\cdots , n) $$
全ての$i=1,\cdots ,n$に対して加算すると、以下のようになる。
$$ \int_{U} (u_{x_{1}}v_{x_{1}}+\cdots +u_{x_{n}}v_{x_{n}} )dx = -\int_{U} u(v_{x_{1}x_{1}}+\cdots v_{x_{n} x_{n}})dx + \int_{\partial U} ( v_{x_{1}}\nu^1 +\cdots v_{x_{n}}\nu^n )udS $$
整理すると、以下のようになる。
$$ \int_{U} Du\cdot Dvdx = -\int_{U} u\Delta vdx + \int_{\partial U} \dfrac{\partial v}{\partial \nu}u dS $$
■
(iii)
(ii)で$u$と$v$の位置を入れ替えると、以下の式を得る。
$$ \int_{U} Du \cdot Dv dx = -\int_{U} v \Delta u dx+\int_{\partial U}\dfrac{\partial u}{\partial \nu}vdS $$
この式から(ii)を引くと、以下のようになる。
$$ 0= -\int_{U} ( v \Delta u -u\Delta v) dx+\int_{\partial U} \left( \dfrac{\partial u}{\partial \nu}v -\dfrac{\partial v}{\partial \nu}u \right)dS $$
整理すると、次を得る。
$$ -\int_{U} ( v \Delta u -u\Delta v) dx=\int_{\partial U} \left( \dfrac{\partial v}{\partial \nu}u -\dfrac{\partial u}{\partial \nu}v \right)dS $$
■