ベクトル積分定理(integral theorems of vector field)

残りの章ではベクトル積分定理とよばれる定理について学びます.これらの定理は簡単にいうと,閉領域の境界での関数の積分を,閉領域内部でのある微分演算子(勾配,発散,回転)を用いた積分で表わします.この関係の特殊な場合にすでに学んだ微分積分学の基本定理があります.つまり区間 $ [a,b]$ の端点での関数 $ f$ の値が,関数 $ f$ の導関数をその区間で積分したものと等しくなるというものです.つまり

$\displaystyle f(b) - f(a) = \int_{a}^{b} \frac{df}{dx} dx $

この場合境界値は $ a$$ b$ で微分演算子は導関数となります.

ここでは,2次元のベクトル場

$\displaystyle {\bf F}(x,y) = (P(x,y), Q(x,y)) $

を平面上の閉曲線 $ C$ がそれ自身で交わらない曲線上で積分します.このような閉曲線を 単一閉曲線(simple closed curve) といいます.

1828年に イギリスの数学者 George Green (1793-1841) によって示された 平面上のGreenの定理 とよばれている定理について考えます.

定理 8..1  

[Greenの定理] 区分的に滑らかな閉曲線 $ C$ を境界にもつ $ xy$ 平面上の有界閉領域を $ \Omega $ とする.またベクトル場 $ {\bf F} = (P(x,y),Q(x,y))$$ C^{1}$ 級とする.このとき

$\displaystyle \oint_{C}{\bf F}\cdot d{\bf r} = \iint_{\Omega}(\nabla \times {\bf F})\cdot \hat{\bf k}dx dy $

$\displaystyle \oint_{C}Pdx + Q dy = \iint_{\Omega}\left (\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y}\right ) dx dy $

が成り立つ.ただし, $ C$ の向きづけは進む方向に対して $ \Omega $ が左側にくるものとする.


図 8.4: Greenの定理
\begin{figure}\begin{center}
\includegraphics[width=5cm]{CALCFIG/Fig8-6-1.eps}
\end{center}\vskip -1cm
\end{figure}

証明

$\displaystyle {\bf F}\cdot d{\bf r} = (P(x,y),Q(x,y))\cdot (dx,dy) = Pdx + Qdy$

$\displaystyle (\nabla \times {\bf F})\cdot \hat{\bf k} = \frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y} $

より,

$\displaystyle \oint_{C}Pdx + Qdy = \iint_{\Omega}(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y})dx dy$

を示せばよいでしょう.

まず,領域 $ \Omega $ がV-simpleであり,またH-simpleでもあるとき標準領域とよぶことにします.つまり標準領域とは

$\displaystyle \Omega$ $\displaystyle =$ $\displaystyle \{(x,y) : a < x < b, y_{1}(x) < y < y_{2}(x) \}$  
  $\displaystyle =$ $\displaystyle \{(x,y) : c < y < d, x_{1}(y) < x < x_{2}(y) \}$  

のような領域です.そこで,まず $ \Omega $ が標準領域の場合について考えてみます.
$\displaystyle \iint_{\Omega}\frac{\partial P}{\partial y}dydx$ $\displaystyle =$ $\displaystyle \int_{a}^{b}\left[\int_{y_{1}(x)}^{y_{2}(x)}\frac{\partial P}{\partial y}dy\right] dx$  
  $\displaystyle =$ $\displaystyle \int_{a}^{b}[P(x,y_{2}(x)) - P(x,y_{1}(x))]dx$  
  $\displaystyle =$ $\displaystyle - \int_{a}^{b}[P(x,y_{1}(x)) - P(x,y_{2}(x))]dx$  

ここで曲線 $ y_{1}$ のグラフを左から右に進むようにパラメター$ t$ を用いて表わすと

$\displaystyle C_{1} : {\bf r}_{1}(t) = (t , y_{1}(t)),   t \in [a,b] $

次に 曲線 $ y_{2}$ のグラフを左から右に進むようにパラメター$ t$ を用いて表わすと

$\displaystyle C_{2} : {\bf r}_{2}(t) = (t , y_{2}(t)),   t \in [a,b] $

曲線 $ C$ は正の方向に進むので, $ C = C_{1} - C_{2}$ となる.よって
$\displaystyle \iint_{\Omega}\frac{\partial P}{\partial y}dydx$ $\displaystyle =$ $\displaystyle - \int_{a}^{b}[P(x,y_{1}(x)) - P(x,y_{2}(x))]dx$  
  $\displaystyle =$ $\displaystyle - (\int_{C_{1}}P(x,y)dx - \int_{C_{2}}P(x,y)dx ) = - \oint_{C}P(x,y)dx$  

となる.次に,
$\displaystyle \iint_{\Omega}\frac{\partial Q}{\partial x}dxdy$ $\displaystyle =$ $\displaystyle \int_{c}^{d}\left[\int_{x_{1}(y)}^{x_{2}(y)}\frac{\partial Q}{\partial x}dx\right] dy$  
  $\displaystyle =$ $\displaystyle \int_{c}^{d}[Q(x_{2}(y),y) - Q(x_{1}(y),y)]dy$  

となるが,ここで曲線 $ x_{1}$ のグラフを下から上に進むようにパラメター$ t$ を用いて表わすと

$\displaystyle C_{3} : {\bf r}_{3}(t) = (x_{1}(t),t),   t \in [c,d] $

次に 曲線 $ x_{2}$ のグラフを下から上に進むようにパラメター$ t$ を用いて表わすと

$\displaystyle C_{4} : {\bf r}_{4}(t) = (x_{2}(t),t),   t \in [c,d] $

曲線 $ C$ は正の方向に進むので, $ C = C_{4} - C_{3}$ となる.よって
$\displaystyle \iint_{\Omega}\frac{\partial Q}{\partial x}dxdy$ $\displaystyle =$ $\displaystyle \int_{c}^{d}[Q(x_{2}(y),y) - Q(x_{1}(y),y)]dy$  
  $\displaystyle =$ $\displaystyle \int_{C_{4}}Q(x,y)dy - \int_{C_{3}}Q(x,y)dy = \oint_{C}Q(x,y)dy$  

よって

$\displaystyle \oint_{C}P(x,y) dx + \oint_{C}Q(x,y) dy = \iint_{\Omega}\frac{\partial Q}{\partial x}dxdy - \iint_{\Omega}\frac{\partial P}{\partial y}dydx $

これより標準領域に対しては,Greenの定理が成り立つ.

次に領域 $ \Omega $ が一般の場合を考えます.この場合 $ \Omega $ を各部分が標準領域になるように分けると,

$\displaystyle \Omega = \Omega_{1} \cup \Omega_{2} \cup \cdots \cup \Omega_{n} $

このとき各部分領域 $ \Omega_{i}$ ではGreenの定理が成り立つので,その境界 $ \partial \Omega_{i}$

$\displaystyle \oint_{\partial \Omega_{i}}Pdx + Qdy = \iint_{\Omega_{i}}(\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y})dx dy $

それらを辺々くわえると右辺は

$\displaystyle \iint_{\Omega} (\frac{\partial Q}{\partial x} - \frac{\partial P}{\partial y})dx dy $

また左辺も $ \Omega $ の境界以外は打ち消しあって

$\displaystyle \oint_{C}Pdx + Qdy $

となるので, $ \Omega $ についてGreenの定理が成り立つ. $  \blacksquare$

例題 8..18  

ベクトル場 $ {\bf F} = (y, -x)$,曲線 $ C$ を原点を中心とする半径1の円とする.このとき曲線 $ C$ に沿っての線積分 $ \oint_{C}{\bf F}\cdot d{\bf r}$ を Greenの定理を用いて求めてみましょう.

曲線 $ C$ $ {\bf r} = {\bf r}(t) = (\cos{t},\sin{t})$ と表わせます.また $ \Omega $ は半径1の円の内部より線積分は

$\displaystyle \oint_{C} {\bf F}\cdot d{\bf r}$ $\displaystyle =$ $\displaystyle \oint_{C} ydx - xdy$  
  $\displaystyle =$ $\displaystyle \iint_{\Omega} \left(\frac{\partial (-x)}{\partial x} - \frac{\partial y}{\partial y} \right ) dx dy$  
  $\displaystyle =$ $\displaystyle \iint_{\Omega}(-2)dxdy = -2(半径1の円の面積) = -2\pi
\ensuremath{\ \blacksquare}$  

例題 8..19  

$ \displaystyle{{\bf F} = \frac{(-y,x)}{x^2 + y^2}}$ のとき, $ \displaystyle{\oint_{C}{\bf F}\cdot d{\bf r}}$ を求めてみましょう.ただし,曲線 $ C$ は例題8.6 と同じとします.

この問題のベクトル場 $ {\bf F}$$ C^{1}$ 級ではありません.よってGreenの定理は使えません.そこで線積分 $ \oint_{C}{\bf F}\cdot d{\bf r}$ の計算を直接行わなければなりません.

$\displaystyle {\bf r} = (x,y) = (\cos{t},\sin{t})$

より
$\displaystyle \oint_{C} {\bf F}\cdot d{\bf r}$ $\displaystyle =$ $\displaystyle \oint_{C}\frac{(-y,x)}{x^2 + y^2 }\cdot (dx,dy)$  
  $\displaystyle =$ $\displaystyle \oint_{C}\frac{-ydx + xdy}{x^2 + y^2}$  
  $\displaystyle =$ $\displaystyle \int_{0}^{2\pi}(\sin^{2}{t} + \cos^{2}{t}) dt = 2\pi
\ensuremath{ \blacksquare}$  

Stokesの定理

Greenの定理を アイルランドの数学者で物理学者の George Gabriel Stokes (1819-1903) が一般化したものをStokesの定理とよびます.まず, Stokesの定理を学ぶには,曲面の向きづけを行なう必要があります.

定義 8..4  

曲面 $ S$ 上の各点 $ (x,y,z)$ で法線ベクトル $ \hat{\bf n}(x,y,z)$ を適当に選び, $ \hat{\bf n}(x,y,z)$$ S$ 上で連続になるようにできるとき,曲面 $ S$向きづけられる曲面(orientable) という.また,このとき各点で選んだ $ \hat{\bf n}$ を,曲面の向きづけにより定まる法線単位ベクトルという.


向きづけられた曲面 $ S$ の境界の曲線 $ \partial S$ に沿っての線積分を, $ S$ 上での面積分に書き換える等式を与えるのがStokesの定理です.

定理 8..2  

[Stokesの定理] $ S:z = f(x,y), (x,y) \in \Omega$ をいくつかの区分的に滑らかな閉曲線を境界とする向きづけられた曲面とする.また,ベクトル場 $ {\bf F}$$ S$ 上で $ C^{1}$ 級とする.そのとき,

$\displaystyle \oint_{\partial S}{\bf F}\cdot d{\bf r} = \iint_{S}(\nabla \times {\bf F})\cdot \hat{\bf n} dS $

ただし, $ \partial S$$ S$ の境界を表わし,曲線 $ \partial S$ 上の線積分の向きは領域 $ S$ を左手にみるように $ \partial S$ を一周するものとする.つまり,法線単位ベクトル $ \hat{\bf n}$ に対して右手の法則に従う.


図 8.5: Stokesの定理
\begin{figure}\begin{center}
\includegraphics[width=5.1cm]{CALCFIG/Fig8-6-2.eps}
\end{center}\end{figure}

証明 まず, $ \iint_{S} [\nabla \times F_{1}\hat{\bf i}] \cdot \hat{\bf n}dS$ を考えよう.

$\displaystyle \nabla \times F_{1}\hat{\bf i} = \left\vert\begin{array}{ccc}
\ha...
...l F_{1}}{\partial z}\hat{\bf j} - \frac{\partial F_{1}}{\partial y}\hat{\bf k} $

より

$\displaystyle [\nabla \times F_{1}\hat{\bf i}] \cdot \hat{\bf n}dS = \left (\fr...
...f j} - \frac{\partial F_{1}}{\partial y}\hat{\bf n} \cdot \hat{\bf k}\right )dS$ (8.2)

$ {\bf r} = (x,y,z)$ を位置ベクトルとすると, $ {\bf r}_{y} = (0,1,z_{y}) $$ S$ の接線ベクトルとなるので,法線ベクトル $ \hat{\bf n}$ とは直交する.よって

$\displaystyle \hat{\bf n} \cdot {\bf r}_{y} = \hat{\bf n} \cdot \hat{\bf j} + z...
...box{または} \ \hat{\bf n} \cdot \hat{\bf j} = - z_{y}\hat{\bf n} \cdot \hat{\bf k}$

これを式8.2 に代入すると,

$\displaystyle [\nabla \times F_{1}\hat{\bf i}] \cdot \hat{\bf n}dS = - \left (\...
... + \frac{\partial F_{1}}{\partial y} \right ) \hat{\bf n} \cdot \hat{\bf k} dS $

ここで $ S$ 上では $ F_{1}(x,y,z) = F_{1}(x,y,f(x,y)) = H(x,y)$ とおけるので合成関数の微分法より

$\displaystyle \frac{\partial H}{\partial y} = \frac{\partial F_{1}}{\partial y} + \frac{\partial F_{1}}{\partial z}\frac{\partial z}{\partial y} $

となる.よって

$\displaystyle [\nabla \times F_{1}\hat{\bf i}] \cdot \hat{\bf n}dS = - \frac{\p...
...tial y}\hat{\bf n} \cdot \hat{\bf k} dS = - \frac{\partial H}{\partial y}dx dy $

を得る.これより

$\displaystyle \iint_{S}[\nabla \times F_{1}\hat{\bf i}] \cdot \hat{\bf n}dS = \iint_{\Omega}- \frac{\partial H}{\partial y}dx dy $

で表わせる.ここで, $ \Omega $$ S$$ xy$ 平面に正射影したものである.右側の積分は平面上の積分で,Greenの定理より $ \oint_{\partial \Omega}{\bf H}dx$ となる.$ \Omega $ の境界上の点 $ (x,y)$ での $ H(x,y)$ の値と, $ S$ の境界上の点 $ (x,y,z)$ での $ F_{1}(x,y,z)$ の値は等しく,また, $ dx$ はどちらの曲線でも同じなので,

$\displaystyle \oint_{\partial \Omega} {\bf H}dx = \oint_{\partial S}{F_{1}}dx $

または

$\displaystyle \iint_{S}[\nabla \times F_{1}\hat{\bf i}] \cdot \hat{\bf n}dS = \oint_{\partial S}{F_{1}}dx $

となる.同様にして,他の平面への正射影をとると,

$\displaystyle \iint_{S}[\nabla \times F_{2}\hat{\bf j}] \cdot \hat{\bf n}dS = \...
...bla \times F_{3}\hat{\bf k}] \cdot \hat{\bf n}dS = \oint_{\partial S}{F_{3}}dz $

となり,くわえると,

$\displaystyle \iint_{S}(\nabla \times {\bf F}) \cdot \hat{\bf n}dS = \oint_{\partial S}{\bf F}\cdot d{\bf r} $

となる. $  \blacksquare$

例題 8..20  

$ \displaystyle{{\bf F} = (-y,x,1),  S: z = (4 - x^2 - y^2)^{1/2}}$ のとき,Stokesの定理が成り立つことを示してみましょう.

$ S$ の境界 $ \partial S$ $ x^2 + y^2 = 4$ の円となります.よって位置ベクトル $ {\bf r} = (x,y,z) = (2\cos{t},2\sin{t},0)$ より線積分を求めると

$\displaystyle \oint_{\partial S}{\bf F}\cdot d{\bf r}$ $\displaystyle =$ $\displaystyle \oint_{\partial S}(-2\sin{t},2\cos{t},1)\cdot (-2\sin{t},2\cos{t},0)dt$  
  $\displaystyle =$ $\displaystyle \int_{0}^{2\pi}[4\sin^{2}{t} + 4\cos^{2}{t}]dt = 8\pi$  

次に面積分を求めてみます.

$\displaystyle \nabla \times {\bf F} = \left\vert\begin{array}{ccc}
\hat{\bf i} ...
...} & \frac{\partial}{\partial z}\\
-y & x & 1
\end{array}\right\vert = (0,0,2) $

次に法線単位ベクトルを求めると

$\displaystyle {\bf r}_{x} \times {\bf r}_{y} = \left\vert\begin{array}{ccc}
\ha...
...rac{y}{z}
\end{array}\right\vert = \left (\frac{x}{z}, \frac{y}{z}, 1 \right ) $

となるので

$\displaystyle \hat{\bf n} = \frac{( \frac{x}{z}, \frac{y}{z}, 1)}{\Vert( \frac{x}{z}, \frac{y}{z}, 1)\Vert } $

これより,
$\displaystyle \iint_{S}(\nabla \times {\bf F}) \cdot \hat{\bf n} dS$ $\displaystyle =$ $\displaystyle \iint_{S}(0,0,2) \cdot \frac{{\bf r}_{x} \times {\bf r}_{y}}{\Ver...
...}_{x} \times {\bf r}_{y} \Vert} \Vert{\bf r}_{x} \times {\bf r}_{y} \Vert dx dy$  
  $\displaystyle =$ $\displaystyle \iint_{\Omega}(0,0,2) \cdot {\bf r}_{x} \times {\bf r}_{y} dx dy$  
  $\displaystyle =$ $\displaystyle \iint_{\Omega}2 dx dy = 2(4\pi) = 8 \pi$  

よってStokesの定理が成り立つことが示せました. $  \blacksquare$

これまでに保存場では,ベクトル場はスカラー場の勾配と等しくなり,またベクトル場の回転は0になることをすでに学びました.では線積分との関係においては,どんなことが成り立つのか調べてみましょう.

定理 8..3  

ベクトル場 $ {\bf F}(x,y,z)$ では次の 3つの条件は同値である.

(1) $ {\bf F} = \nabla \phi$ となるスカラー関数 $ \phi(x,y,z)$ が存在する.(Fは保存場)

(2) いたるところ $ \nabla \times {\bf F} = {\bf0}$ が成り立つ.(渦なし)

(3) 任意の閉曲線 $ C$ について $ \oint_{C}{\bf F} \cdot d {\bf r} = 0$ が成り立つ (積分経路無関係).


証明 (1) $ \Rightarrow$ (2) $ \Rightarrow$ (3) $ \Rightarrow$ (1)を示す.

(1) $ \Rightarrow$ (2) $ \nabla \times {\bf F} = \nabla \times \nabla \phi = \left \vert \begin{array}{...
...z}\\
\phi_{x} & \phi_{y} & \phi_{z}
\end{array}\right \vert = (0,0,0) = {\bf0}$

(2) $ \Rightarrow$ (3) 閉曲線 $ C$ で囲まれた曲面 $ S$ を考えて,Stokesの定理を使うと

$\displaystyle \oint_{C} {\bf F} \cdot d {\bf r} = \iint_{S} (\nabla \times {\bf F}) \cdot \hat{\bf n} dS = \iint_{S} {\bf0} \cdot \hat{\bf n}dS = 0 $

(3) $ \Rightarrow$ (1) 定点P $ (x_{0},y_{0},z_{0})$ と動点Q$ (x,y,z)$ をむすぶ2つの曲線 $ C_{1},C_{2}$ をとり,Pから $ C_{1}$ を経てQに至り,Qから $ C_{2}$ を逆向きに通ってPに戻る道を $ C$ とすると,

$\displaystyle \int_{{\rm P}(C_{1}){\rm Q}}{\bf F} \cdot d{\bf r} - \int_{{\rm P...
...m Q}(C_{2}){\rm P}}{\bf F} \cdot d{\bf r} = \oint_{C}{\bf F} \cdot d{\bf r} = 0$

よって,

$\displaystyle \int_{{\rm P}(C_{1}){\rm Q}}{\bf F} \cdot d{\bf r} = \int_{{\rm P}(C_{2}){\rm Q}}{\bf F} \cdot d{\bf r} $

すなわち,PからQに至る線積分 $ \int_{{\rm PQ}}{\bf F} \cdot d{\bf r}$ は途中の経路に関係なく,終点Qの座標 $ (x,y,z)$ の関数で与えられる.よってこれを $ \phi(x,y,z)$ とすれば,


$\displaystyle \int_{\rm PQ}{\bf F} \cdot d{\bf r} = \phi(x,y,z) $


PからQに至る任意の曲線のベクトル方程式を $ {\bf r} = (x(s),y(s),z(s))$ とすると,

$\displaystyle \int_{{\rm PQ}}{\bf F} \cdot d{\bf r} = \int_{s_{0}}^{s} {\bf F} ...
... \left (F_{1}\frac{dx}{ds} + F_{2}\frac{dy}{ds} + F_{3}\frac{dz}{ds}\right )ds $

よって

$\displaystyle F_{1}\frac{dx}{ds} + F_{2}\frac{dy}{ds} + F_{3}\frac{dz}{ds} = \f...
...hi}{\partial y} \frac{dy}{ds} + \frac{\partial \phi}{\partial z} \frac{dz}{ds} $

曲線PQは任意,したがって $ x(s),y(s),z(s)$ も任意の関数でよいから,

$\displaystyle F_{1} = \frac{\partial \phi}{\partial x}, F_{2} = \frac{\partial \phi}{\partial y}, F_{3} = \frac{\partial \phi}{\partial z} $

これより,

$\displaystyle {\bf F} = (F_{1},F_{2},F_{3}) = (\frac{\partial \phi}{\partial x}...
...y}, \frac{\partial \phi}{\partial z}) = \nabla \phi
\ensuremath{ \blacksquare}$

この定理より,線積分をおこなうときに,ベクトル場がスカラーポテンシャルを持てば,積分をしなくとも答は0であることが分かります.

例題 8..21  

$ \displaystyle{\int_{C}(2x+yz,zx,xy) \cdot d{\bf r}}$ を求めてみましょう.ただし, $ C$ は点 $ (1,0,-1)$ から点 $ (2,-1,3)$ に至る曲線.

$\displaystyle \nabla \times (2x+yz,zx,xy) = \left\vert\begin{array}{ccc}
\hat{\...
...frac{\partial}{\partial z}\\
2x+yz & zx & xy
\end{array}\right\vert = (0,0,0) $

より,

$\displaystyle \int_{C}(2x+yz,zx,xy) \cdot d{\bf r} = 0
\ensuremath{ \blacksquare}$

次にドイツの数学者 Karl Friedrich Gauss (1777-1855) の名前をとってつけられた発散定理について学びます.

Gaussの発散定理

定理 8..4  

[Gaussの発散定理] ベクトル場 $ {\bf F}(x,y,z) = (F_{1},F_{2},F_{3})$ において,区分的に滑らかな閉曲面 $ S$ で囲まれた空間の領域を $ V$ とし, $ S$ の内部から外部に向かう法線ベクトルを $ \hat{\bf n}$ とすると,

$\displaystyle \iint_{S}{\bf F} \cdot \hat{\bf n}dS = \iiint_{V} {\rm div}{\bf F} dV $

$\displaystyle \iint_{S}(F_{1}dydz + F_{2}dzdx + F_{3}dxdy) = \iiint_{V} \left(\...
...artial F_{2}}{\partial y} + \frac{\partial F_{3}}{\partial z} \right) dx dy dz $

が成り立つ.


図 8.6: Gaussの定理
\begin{figure}\begin{center}
\includegraphics[width=6cm]{CALCFIG/Fig8-6-3.eps}
\end{center}\end{figure}

証明 まず, $ V$ が2つの曲面 $ S_{1},S_{2}$ で下と上からはさまれているとします.また, $ S_{1}$ $ z = f_{1}(x,y),  (x,y) \in \Omega$$ S_{2}$ $ z = f_{2}(x,y)  (x,y) \in \Omega$ で与えられているとします.このとき,

$\displaystyle \iiint_{V}\frac{\partial F_{3}}{\partial z} dV$ $\displaystyle =$ $\displaystyle \iiint_{V}\frac{\partial F_{3}}{\partial z} dzdydx = \iint_{\Omega}[\int_{z=f_{1}(x,y)}^{f_{2}(x,y)}\frac{\partial F_{3}}{\partial z} dz] dydx$  
  $\displaystyle =$ $\displaystyle \iint_{\Omega}\left[F_{3}\right ]_{z = f_{1}(x,y)}^{f_{2}(x,y)}dydx$  
  $\displaystyle =$ $\displaystyle \iint_{\Omega}[F_{3}(x,y,f_{2}(x,y)) - F_{3}(x,y,f_{1}(x,y)) ]dy dx$  

曲面 $ S_{2}$ において,曲線座標 $ (x,y)$ に対する法線単位ベクトル

$\displaystyle \frac{{\bf r}_{x} \times {\bf r}_{y}}{\Vert{\bf r}_{x} \times {\bf r}_{y}\Vert} $

$ \hat{\bf n}$ と一致しているが,曲面 $ S_{1}$ では法線単位ベクトルは $ -\hat{\bf n}$ に等しい.よって

$\displaystyle \iint_{\Omega}F_{3}(x,y,f_{2}(x,y)) dx dy = \iint_{S_{2}}F_{3} dx dy $

$\displaystyle -\iint_{\Omega}F_{3}(x,y,f_{1}(x,y)) dx dy = \iint_{S_{1}}F_{3} dx dy $

これより,

$\displaystyle \iiint_{V}\frac{\partial F_{3}}{\partial z} dV = \iint_{S_{2}}F_{3} dx dy + \iint_{S_{1}}F_{3} dx dy = \iint_{S}F_{3} dx dy $

が成り立つ.同様にして, $ S$ を他の平面に正射影することにより

$\displaystyle \iiint_{V}\frac{\partial F_{1}}{\partial z} dV = \iint_{S}F_{1} dy dz $

$\displaystyle \iiint_{V}\frac{\partial F_{2}}{\partial z} dV = \iint_{S}F_{2} dz dx $

を示すことができるので,これらをそれぞれくわえれば,

$\displaystyle \iiint_{V} \nabla \cdot {\bf F} dV = \iint_{S}{\bf F} \cdot \hat{\bf n}dS $

を得ることができる.

領域 $ V$ が一般な場合には,平面におけるGreenの定理の証明のように, $ V$ を部分領域に分割して証明すればよい. $  \blacksquare$

例題 8..22  

Gaussの発散定理により面積分

$\displaystyle \iint_{S} xz^2 dydz + (x^2 y - z^3)dzdx + (2xy + y^2 z)dxdy,  S : z = \sqrt{a^2 - x^2 - y^2}, z = 0$

を求めよ.

$ {\bf F} = (xz^2, x^2 y - z^3, 2xy + y^2 z)$ よりdivFを求めると

$\displaystyle {\rm div}{\bf F} = \nabla \cdot (xz^2, x^2 y - z^3, 2xy + y^2 z) = z^2 + x^2 + y^2 $

よって

$\displaystyle \iint_{S} xz^2 dydz + (x^2 y - z^3)dzdx + (2xy + y^2 z)dxdy = \iiint_{V}(x^2 + y^2 + z^2)dV $

ここで $ V$ は半径 $ a$ の上半球より球面座標変換 $ (\rho,\phi,\theta)$ を用いると
$\displaystyle \int_{\theta = 0}^{2\pi}\int_{\phi = 0}^{\frac{\pi}{2}}\int_{\rho = 0}^{a}\rho^2 \rho^2 \sin{\phi}d\rho d\phi d\theta$ $\displaystyle =$ $\displaystyle \int_{\theta = 0}^{2\pi} d\theta \int_{\phi = 0}^{\frac{\pi}{2}} \sin{\phi} \int_{\rho = 0}^{a}\rho^4 d\rho$  
  $\displaystyle =$ $\displaystyle 2\pi(1)(\frac{a^5}{5}) = \frac{2\pi a^5}{5}
\ensuremath{ \blacksquare}$  

演習問題


1.
次の線積分を求めよう.

(a) $ \displaystyle{\oint_{C}(x^2 - xy^3)dx + (y^2 - 2xy) dy}$, ただし, $ C$ は点 $ (0,0),(2,0),(2,2),(0,2)$ を頂点とする正方形.

(b) $ \displaystyle{\oint_{C}(2xy^3 - y^2 \cos{x})dx + (1 - 2y \sin{x} + 3x^2 y^2) dy}$, ただし, $ C$ は点 $ (0,0)$ から点 $ \displaystyle{(\frac{\pi}{2},1)}$ を放物線 $ 2x = \pi y^2$ に沿って進む.

(c) $ \displaystyle{\oint_{\partial S} -z^2 dx + xy^2 dy + z dz}$, ただし, $ \displaystyle{S : z = \sqrt{1 - (x^2 + y^2)}}$

2.
次の面積分を求めよう.

(a) $ \displaystyle{\iint_{S}(x^2 + y^2) dS}$, ただし, $ \displaystyle{S : z^2 = 3(x^2 + y^2), 0 \leq z \leq 3}$

(b) $ \displaystyle{\iint_{S}{\rm curl}{\bf F} \cdot \hat{\bf n} dS}$, ただし, $ \displaystyle{{\bf F} = (x^2 -x,-xy, 3z), S : z = 2 - \sqrt{x^2 + y^2}, z \geq 0}$

(c) $ \displaystyle{\iint_{S}{\bf F} \cdot \hat{\bf n} dS}$, ただし, $ \displaystyle{{\bf F} = (x-z,x^3 + yz, -3xy^3)}$, $ \displaystyle{S : z = 4 - y^2}$, $ x = 0, x=3, z = 0$

(d) $ \displaystyle{\iint_{S}x dydz + ydzdx + zdxdy}$, ただし, $ S$ は 円柱 $ \displaystyle{x^2 + y^2 = 9}$ と平面 $ z = 0, z= 3$ で囲まれた領域

3.
単連結な閉曲線 $ C$ で囲まれた領域の面積$ A$ は, $ \displaystyle{\frac{1}{2}\oint xdy - ydx}$ で与えられることを示そう.
4.
楕円 $ \displaystyle{x = a\cos{\theta}, y = b\sin{\theta}}$ の面積を求めよう.
5.
$ f,g$ をスカラー場とするとき次の式が成り立つことを示そう.

$ \displaystyle{\iint_{S}f \frac{\partial g}{\partial n} dS = \iiint_{V}(f \nabla^2 g + {\rm grad}f \cdot {\rm grad}g)dV }$ここで $ \displaystyle{\frac{\partial f}{\partial n}, \frac{\partial g}{\partial n}}$ はそれぞれ$ f,g$$ S$ における外向き法線方向の方向微分係数を表わす.