発散における基本定理
発散における基本定理は以下の式が成り立つことを主張している。
∫V(∇⋅v)dτ=∮Sv⋅da
この式は非常に重要なので、敬意を表してGauss の定理(Gauss’s theorem)、Green の定理(Green’s theorem)、発散定理(divergence theorem)という3つの名前が付けられている。他の基本定理のように、この定理もある領域(今は体積V)における導関数(今は発散)の積分は境界(今は体積Vの表面S)での関数の値に等しいということをあらわしている。境界項はそれ自身が積分(今は面積分)になっているということに注意せよ。このことは、ある線の境界が2つの端点で定まるのに対して、体積の境界が閉局面になるということからも丁度良いということが伺える。
幾何学的には、もしvが非圧縮性流体の流れをあらわすとき、vの流束((56)の右辺)は単位時間当たりに表面を通過する流体の総量をあらわしている。今、発散は流体が注がれる蛇口に値するようなある点からのベクトルの湧き出しを測っている。非圧縮性の流体で満たされた領域に沢山の蛇口があったとすると、その領域の境界から同量の液体が押し出されることになる。実際には、その量を計算する方法は、それぞれの蛇口からどれぐらい注がれるのかを数えてそれらを足し合わせるというやり方と、境界の周りの各点でどれだけ押し出されるのかを測って足し合わせるという2通りが存在する。これによって、発散定理のエッセンスは以下のような主張であると言える。
∫(考えている体積内の蛇口)=∮(それぞれの面の流れ)
問題
図1.29のように、各辺の長さが1の単位立方体を用いて、ベクトル関数v=y2ˆx+(2xy+z2)ˆy+(2yz)ˆzについて、発散における基本定理を確認せよ。
図1.29
解説
この場合、∇⋅v=2(x+y)であり、
∫V2(x+y)dτ=2∫10∫10∫10(x+y)dxdydz
∫10(x+y)dx=12+y, ∫10(12+y)=1, ∫101dz=1
となる。従って、発散定理の左辺は以下のように計算出来る。
∫V∇⋅vdτ=2
発散定理の右辺を計算するためには、各面で面積分を行って、その6面の結果を足し合わせればよい。
(1)
∫v⋅da=∫10∫10y2dydz=13
(2)
∫v⋅da=−∫10∫10y2dydz=−13
(3)
∫v⋅da=∫10∫10(2x+z2)dxdz=43
(4)
∫v⋅da=−∫10∫10z2dxdz=−13
(5)
∫v⋅da=∫10∫102ydxdy=1
(6)
∫v⋅da=−∫10∫100dxdy=0
故に、全流束は以下のように計算出来る。
∮Sv⋅da=13−13+43−13+1+0=2
これは期待した結果になっている。