logo

グリーンの定理 📂偏微分方程式

グリーンの定理

定理

u,v∈C2(Uˉ)u, v \in C^2( \bar{U})としよう。それでは、以下の式が成り立つ。

  • (i) ∫UΔudx=∫∂U∂u∂νdS\displaystyle \int_{U} \Delta u dx=\int_{\partial U} \dfrac{\partial u}{\partial \nu}dS

  • (ii) ∫UDv⋅Dudx=−∫UuΔvdx+∫∂U∂v∂νudS\displaystyle \int_{U} Dv \cdot Du dx = -\int_{U} u \Delta v dx+\int_{\partial U}\dfrac{\partial v}{\partial \nu}udS

  • (iii) ∫U(uΔv−vΔu)dx=∫∂U(∂v∂νu−∂u∂νv)dS\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という。

証明

部分積分公式

u,v∈C1(Uˉ)u, v \in C^1(\bar{U})とする。それでは、以下の式が成り立つ。

∫Uuxivdx=−∫Uuvxidx+∫∂UuvνidS(i=1,…,n) \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)

部分積分公式において、uuの代わりにuxiu_{x_{i}}を、vvの代わりに11を代入すると、以下の式を得る。

∫Uuxixidx=∫∂UuxiνidS(i=1,⋯ ,n) \int_{U} u_{x_{i} x_{i}}dx = \int_{\partial U} u_{x_{i}}\nu^{i} dS \quad (i=1,\cdots , n)

全てのi=1,⋯ ,ni=1,\cdots, nに対して加算すると、以下のようになる。

∫U(ux1x1+⋯+uxnxn)dx=∫∂U(ux1ν1+⋯uxnνn)dS \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

ラプラシアンの定義と∂u∂ν:=ν⋅Du\dfrac{\partial u}{\partial \nu}:=\boldsymbol{\nu}\cdot Duにより、次が成り立つ。

∫UΔudx=∫∂U∂u∂νdS \int_{U} \Delta u dx=\int_{\partial U} \dfrac{\partial u}{\partial \nu}dS

■

(ii)

部分積分公式において、vvの代わりにvxiv_{x_{i}}を代入すると、以下の式を得る。

∫Uuxivxidx=−∫Uuvxixidx+∫∂UuvxiνidS(i=1,⋯ ,n) \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,⋯ ,ni=1,\cdots ,nに対して加算すると、以下のようになる。

∫U(ux1vx1+⋯+uxnvxn)dx=−∫Uu(vx1x1+⋯vxnxn)dx+∫∂U(vx1ν1+⋯vxnνn)udS \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

整理すると、以下のようになる。

∫UDu⋅Dvdx=−∫UuΔvdx+∫∂U∂v∂νudS \int_{U} Du\cdot Dvdx = -\int_{U} u\Delta vdx + \int_{\partial U} \dfrac{\partial v}{\partial \nu}u dS

■

(iii)

(ii)でuuとvvの位置を入れ替えると、以下の式を得る。

∫UDu⋅Dvdx=−∫UvΔudx+∫∂U∂u∂νvdS \int_{U} Du \cdot Dv dx = -\int_{U} v \Delta u dx+\int_{\partial U}\dfrac{\partial u}{\partial \nu}vdS

この式から(ii)を引くと、以下のようになる。

0=−∫U(vΔu−uΔv)dx+∫∂U(∂u∂νv−∂v∂νu)dS 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

整理すると、次を得る。

−∫U(vΔu−uΔv)dx=∫∂U(∂v∂νu−∂u∂νv)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

■

参照