記事ソース/div
メニュー現在 8 名がオンラインです。 最新の25件2023-12-12
2023-11-11
2023-11-06
2023-07-15
2022-09-14
2022-07-01
2022-06-12
2022-04-13
2021-12-03
2021-10-07
2021-08-12
2021-07-26
2021-06-30
2021-06-06
2021-05-02
2021-04-17
2021-03-20
2021-03-19
2021-03-11
|
記事ソース/div†これはrst2hooktailの記事ソース保存・変換用です（詳細）． コンバート公開・更新メニュー ▼▲記事ソースの内容
=======================================
div
=======================================
スカラー場の勾配を考えたとき、ベクトル微分演算子 $\nabla = \left( \frac{\partial}{\partial x_{1}},\frac{\partial}{\partial x_{2}},\frac{\partial}{\partial x_{3}} \right)$ というものを導入しました。そして、この $\nabla$ をスカラー関数に作用させたものを勾配（ ${\rm grad}$ ）と呼びました。 $\nabla$ をベクトル関数と内積を取る形で作用させたものを *発散* と呼びます。英語で発散を $divergence$ と言うので、記号 ${\rm div}$ を使う場合もあります。
<tex>
{\rm div} \bm{A} &= \nabla \cdot \bm{A} \\
&= \frac{\partial A_{1}}{\partial x_{1}}+\frac{\partial A_{2}}{\partial x_{2}}+\frac{\partial A_{3}}{\partial x_{3}} \tag{1}
</tex>
ただし、式中 $\bm{A}(x_{1},x_{2},x_{3})=(A_{1}(x_{1},x_{2},x_{3}),A_{2}(x_{1},x_{2},x_{3}),A_{3}(x_{1},x_{2},x_{3}))$ としました。単に $\nabla$ とベクトルの内積になっているだけなので、計算が特に難しいということは無いと思います。
発散のイメージ
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
ここでは、デカルト座標系を考えます。ベクトル場が、点 $M$ において $\bm{A}(x,y,z)$ と表現されるとき、 $M$ を一つの頂点に、各辺の長さが $\Delta x,\Delta y,\Delta z$ の直方体領域を考えます。
.. image:: Joh-Divergence002.gif
ベクトル場を、水（非圧縮流体）の流れだと考えると状況がイメージしやすいでしょう。この直方体領域は流れの中に置かれていますから、絶えず水が流れ込んだり出て行ったりしています。しかし、水は非圧縮流体だと仮定していますので、普通なら、入ってくる水量と出て行く水量は同じはずです。ところが、もしこの領域内に"温泉の湧き出し口"のように、絶えず水を噴出している穴が空いているとすればどうでしょうか？全体で、出て行く水の方が多くなりそうですね。逆に、水を吸い込んでしまう排水口のような穴が空いているとしたら、流入する水量の方が多くなりそうです。こうした点を *湧き出し* 、 *吸い込み* と呼びます。
.. [*] 湧き出し口や排水口に水道パイプはつながっているのか？などと、現実的なことを考えてはいけません。数学的に、空間のある一点が（上下水道なしでいきなり）水を噴出したり吸い込んだりしていると仮定しての話です。宇宙のどこかには、ブラックホールやホワイトホールと呼ばれる、湧き出しや吸い込みがあるようですが、そんなイメージでお願いします。
例として、 $x$ 軸方向に流入する量と流出する量の差を考えてみます。図のように、頂点 $A$ 付近に単位面積当たり流入する流体の量は $V_{x}$ で表わされ、頂点 $A'$ 付近では $V_{x}+\frac{\partial V_{x}}{\partial x}\Delta x$ と表わされます。その差は $\left( V_{x}+\frac{\partial V_{x}}{\partial x}\Delta x \right) - V_{x}= \frac{\partial V_{x}}{\partial x}\Delta x$ です。
.. image:: Joh-Divergence2.gif
ここまでは偏微分を定義通りに使っただけです。面倒なのは、正確に表現しようとすると『 $x$ 方向への変化』されも、 $y$ 軸や $z$ 軸に沿って変化することです。つまり、頂点 $B$ 付近で $x$ 方向の流入流は $V_{x}+\frac{\partial V_{x}}{\partial y}\Delta y$ と表現され、さらに $B'$ 付近ではこの値の $x$ 方向への変化を考え、流出流は次式のように表現されることになります。（ここで図に全て書き込むとごちゃごちゃしますので描きませんでした。大事なところですので、自分で似たような図でも描いて、しっかり考えてみて下さい。）
<tex>
V_{x}+\frac{\partial V_{x}}{\partial y}\Delta y + \frac{\partial}{\partial x}\left(V_{x}+ \frac{\partial V_{x}}{\partial y} \Delta y \right) \Delta x = V_{x}+\frac{\partial V_{x}}{\partial y}\Delta y +\frac{\partial V_{x}}{\partial x}\Delta x + \frac{\partial^2 V_{x}}{\partial x \partial y} \Delta y \Delta x \tag{2}
</tex>
ここで、直方体領域は十分に小さいとすれば、 $x$ 方向の流れは "ほぼ一様"だと考えて、二次の変化量（右辺第三項）を無視できます。結局、頂点 $B$ 付近での流量の差は $\left( V_{x}+\frac{\partial V_{x}}{\partial y}\Delta y \right) - \left( V_{x}+\frac{\partial V_{x}}{\partial y}\Delta y +\frac{\partial V_{x}}{\partial x}\Delta x \right) = \frac{\partial V_{x}}{\partial x}\Delta x $ と近似できることになります。これは頂点 $A$ 付近の流れと同じになります。 $z$ 方向に関しても同様に $\left( V_{x}+\frac{\partial V_{x}}{\partial z}\Delta z \right) - \left( V_{x}+\frac{\partial V_{x}}{\partial z}\Delta z +\frac{\partial V_{x}}{\partial x}\Delta x \right) = \frac{\partial V_{x}}{\partial x}\Delta x $ が言えますから、結局、面 $ABCD$ に流入する流体と、面 $A'B'C'D'$ から流出する流体の差は、単位面積あたり $\frac{\partial V_{x}}{\partial x}\Delta x$ と書いてしまって良いことになります。
.. important::
微小直方体領域で、 $x$ 軸方向に『 $($ 流出した量 $)-($ 流入した量 $)$ 』を、 $y$ 軸方向や $z$ 軸方向への二次的な変化を無視すると、単位面積あたり $\frac{\partial V_{x}}{\partial x}\Delta x$ と表わせる
ここで $x$ 軸方向の流れだけを考えて『 $($ 流出した量 $)-($ 流入した量 $)=($ 内部で湧き出した量 $)$ 』と考えてしまっては早計です。 $y$ 軸方向や $z$ 軸方向からの流体の出入りも考慮しなければならないからです。しかし、もし全方向で『 $($ 流出した量 $)-($ 流入した量 $)$ 』を考えれば、これは『 $($ 内部で湧き出した量 $)$ 』に等しくなります。
この直方体の $x$ 軸方向の断面積（上図の四角形 $ABCD$ の面積）が $\Delta y \Delta z$ であることに注意すると、 $x$ 軸方向の流量の差 $\frac{\partial V_{x}}{\partial x}\Delta x \Delta y \Delta z$ と表わされます。同様に、 $y$ 軸方向の流量の差は $\frac{\partial V_{y}}{\partial y}\Delta x \Delta y \Delta z$ 、 $y$ 軸方向の流量の差は $\frac{\partial V_{z}}{\partial z}\Delta x \Delta y \Delta z$ となり、この直方体全体から流出する流量はこれらを足して次式になります。
<tex>
\frac{\partial V_{x}}{\partial x}\Delta x \Delta y \Delta z + \frac{\partial V_{y}}{\partial y}\Delta x \Delta y \Delta z + \frac{\partial V_{z}}{\partial z}\Delta x \Delta y \Delta z \tag{3}
</tex>
最後に、これを直方体の体積 $\Delta x \Delta y \Delta z$ で割れば、 *単位体積当たりの湧き出し量* の表式を得ます。これは、もうお気づきのように ${\rm div} \bm{V}$ になっています。
.. important::
単位体積当たりの湧き出しは ${\rm div}\bm{V}=\frac{\partial V_{x}}{\partial x}+\frac{\partial V_{y}}{\partial y}+\frac{\partial V_{z}}{\partial z}$ で表わされます。
ベクトル場 $\bm{V}$ が水の流れ以外のベクトル場であっても、流れの物理的意味は少し異なるかも知れませんが、基本的には『単位体積当たりの流れの増減』を表わすのが発散だと考えることができます。このようなイメージを持っていれば、電磁気学で電場や磁場といった目に見えない場の計算に ${\rm div}$ が出てきても、慌てることは無いと思います。
また、 ベクトル解析奮闘記２_ では、以上の議論が直観的に分かりやすい絵と、味のある文で学習者の視点から説明されています。ぜひ、そちらも参照してください！
.. [*] 発散の意味の物理的意味は分かったでしょうか。流体力学で、湧き出しや吸い込みの無い非圧縮流体の流れは、 ${\rm div}\bm{u}=\bm{0}$ と表現されることになります。この式を見て、『非圧縮流で、湧き出しの無い流れだな』とすぐに読めるくらいになると素敵です。
発散の演算規則
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
発散をとる操作も線形演算ですので、以下のような算法がなりたちます。発散も線形ですね。
1. $\nabla \cdot (\alpha \bm{A}+\beta \bm{B})=\alpha \nabla \cdot \bm{A}+\beta \nabla \cdot \bm{B}$
2. $\nabla (\phi \bm{A})=(\nabla \phi)\cdot \bm{A}+\phi\nabla \cdot \bm{A}$
歴史
------------------------------------------------------------------
ベクトルは 四元数_ から発生したものですが、四元数の基底 $i,j,k$ を使ったベクトルの積は次式のようになります。（詳しくは 七次元の外積_ を参照下さい。）
<tex>
(iA_{1}+jA_{2}+kA_{3})(iB_{1}+jB_{2}+kB_{3})
&= -(A_{1}B_{1}+A_{2}B_{2}+A_{3}B_{3}) \\
& \ \ \ + i(A_{2}B_{3}-A_{3}B_{2}) + j(A_{3}B_{1}-A_{1}B_{3})+k(A_{1}B_{2}-A_{2}B_{1})
</tex>
この右辺第一項は、私達が知っている内積にマイナスがついた形をしています。ベクトル解析が生まれた頃は、どうしてもこの四元数の乗法が念頭にあったためか、内積にマイナスをつける流儀の人が多く、例えば電磁気学のマックスウェル( $\text{James Clerk Maxwell (1831-1879)}$ )は、四元数の基底を使って $\nabla = i \frac{\partial}{\partial x}+ j\frac{\partial}{\partial y}+k\frac{\partial}{\partial z}$ という記号を使っていました。この記号をベクトル $iA_{1}+jA_{2}+kA_{3}$ に作用させたときのスカラー項を $convergence$ と呼び、当時は ${\rm conv}$ などと書いていました。
<tex>
{\rm conv}\bm{A} = -\frac{\partial A_{1}}{\partial x}- \frac{\partial A_{2}}{\partial y}- \frac{\partial A_{3}}{\partial z}
</tex>
これは、ちょうど私達がいま使っている発散 $(divergence)$ と符号が逆ですから、物理的には吸い込みを意味することになります。 $\nabla$ を現在の形にし、ずっと扱い易くしたのはヘヴィサイド( $\text{Oliver Heaviside (1850-1925)}$ )の功績です。現在では ${\rm conv}$ なんて使う人は誰もいませんし、ベクトル解析の表記はとてもすっきりと自然な形にまとまっています。数学の進歩のためには美しい表記法が重要だというのは、微積分の記号を発明したライプニッツ( $\text{Gottfried Wilhelm Leibniz (1646-1716)}$ )の信念だったそうですが、 $\nabla$ や後述する $\triangle$ などは確かによく出来た記号で、こうした記号のお陰で、ベクトル解析の公式の中には芸術的とも言える美しいものがあるのです。
.. _七次元の外積: http://www12.plala.or.jp/ksp/vectoranalysis/SevenDCrossProd/
.. _四元数: http://www12.plala.or.jp/ksp/mathInPhys/quaternion/
.. _ベクトル解析奮闘記: http://www12.plala.or.jp/ksp/vectoranalysis/vecFuntou1/
.. _rotgrad=0,divrot=0: http://www12.plala.or.jp/ksp/vectoranalysis/ddw=0/
.. _ベクトル解析奮闘記２: http://www12.plala.or.jp/ksp/vectoranalysis/vecFuntou2/
@@author:Joh@@
@@accept: 2006-10-11@@
@@category: ベクトル解析@@
@@id: VectorDivergence@@
|
