ベクトル解析 における発散 （はっさん、英 : divergence ベクトル場 の各点ごとの流入出 （英語版 ） ベクトル作用素 である。より技術的に言えば、対象点を含む近傍領域を定義しそこに出入りする流束 の総和と領域体積との比をとり、領域を無限小に近づけたときの極限である。
身近な例えでは、温度変化のある空気の各点の移動速度 ベクトル場をみる。一部の領域の空気を熱するとその膨張した空気は領域から全方向へ広がるから、領域の外側を向く速度場が生じる。このときの速度場の発散をとると、加熱された領域の内部で正値の分布であり、この領域は速度場全体にとっての流入（あるいは正の湧き出し、湧出、source ）域である。逆に空気が冷やされ収縮するならば、冷却される領域の発散は負値となり、その領域は流出（負の湧き出し、あるいは沈み込み、排出、sink ）域である。
物理的な言葉で言えば、三次元ベクトル場の発散は、各点においてそのベクトル場が流入や流出のような流動的振舞いをする度合を与える。これは、空間の無限小領域において入ってくるよりも出ていく方がどのくらい多いのかの度合いとしての「外向き度」を局所的に測るものである。発散がその点で零でないならば、その位置は湧出点か排出点でなければならない[1]
点 p におけるベクトル場 F Ω の滑らかな境界 bd(Ω) と交わる F Ω の体積 vol(Ω) で割ったものの、領域 Ω を一点 p に縮めるときの極限として定義される。これを式で書けば
(
div
F
)
p
:=
lim
Ω
→
{
p
}
1
vol
(
Ω
)
∮
bd
(
Ω
)
(
F
⋅
n
)
d
S
{\displaystyle (\operatorname {div} {\boldsymbol {F}})_{p}:=\lim _{\Omega \to \{p\}}{\frac {1}{\operatorname {vol} (\Omega )}}\oint _{\operatorname {bd} (\Omega )}({\boldsymbol {F}}\cdot {\boldsymbol {n}})\,dS}
となる。積分は境界面に直交する外向きの単位法ベクトル 場 n 面積分 である。各点において div F が得られて、これは点 p div F が F 流出密度 （発散密度 ）であることが分かる。
物理的な解釈から見れば、あらゆる点で発散がゼロとなるベクトル場は非圧縮性 (incompressible) あるいは管状 (solenoidal)
直感的に、全ての湧出量の和から全ての排出量の和を引けば領域から流れ出る正味の流れがわかるだろうと想像できる。これを精緻化したものが発散定理 である。
具体的な表示
編集
デカルト座標系での表示
編集
x , y , z ユークリッド空間 のデカルト座標系 とし、対応する単位ベクトル からなる基底 を i j k
連続的微分可能 なベクトル場 F U i V j W k
div
F
=
∇
⋅
F
=
∂
U
∂
x
+
∂
V
∂
y
+
∂
W
∂
z
{\displaystyle \operatorname {div} {\boldsymbol {F}}=\nabla \cdot {\boldsymbol {F}}={\frac {\partial U}{\partial x}}+{\frac {\partial V}{\partial y}}+{\frac {\partial W}{\partial z}}}
に等しい。これは座標で表されているけれども、物理的解釈が示唆する通り、この式の値は任意の直交変換 によって変わることはない。
しばしば用いられる発散の記法 “∇ · F ” は、中黒を点乗積 と見做して、∇ の成分（ナブラ の項を参照）と F 記号の濫用 の一種である。
連続的微分可能二階テンソル場 ε
div
(
ε
)
→
=
(
∂
ε
x
x
∂
x
+
∂
ε
x
y
∂
y
+
∂
ε
x
z
∂
z
∂
ε
y
x
∂
x
+
∂
ε
y
y
∂
y
+
∂
ε
y
z
∂
z
∂
ε
z
x
∂
x
+
∂
ε
z
y
∂
y
+
∂
ε
z
z
∂
z
)
{\displaystyle {\overrightarrow {\operatorname {div} ({\boldsymbol {\varepsilon }})}}={\begin{pmatrix}{\dfrac {\partial \varepsilon _{xx}}{\partial x}}+{\dfrac {\partial \varepsilon _{xy}}{\partial y}}+{\dfrac {\partial \varepsilon _{xz}}{\partial z}}\\[5pt]{\dfrac {\partial \varepsilon _{yx}}{\partial x}}+{\dfrac {\partial \varepsilon _{yy}}{\partial y}}+{\dfrac {\partial \varepsilon _{yz}}{\partial z}}\\[5pt]{\dfrac {\partial \varepsilon _{zx}}{\partial x}}+{\dfrac {\partial \varepsilon _{zy}}{\partial y}}+{\dfrac {\partial \varepsilon _{zz}}{\partial z}}\end{pmatrix}}}
になる。
円柱座標系
編集
a -方向の単位ベクトルを e a
F
=
e
r
F
r
+
e
z
F
z
+
e
θ
F
θ
{\displaystyle {\boldsymbol {F}}={\boldsymbol {e}}_{r}F_{r}+{\boldsymbol {e}}_{z}F_{z}+{\boldsymbol {e}}_{\theta }F_{\theta }}
に対し、その発散は、
div
F
=
∇
⋅
F
=
1
r
∂
∂
r
(
r
F
r
)
+
1
r
∂
F
θ
∂
θ
+
∂
F
z
∂
z
{\displaystyle \operatorname {div} {\boldsymbol {F}}=\nabla \cdot {\boldsymbol {F}}={\frac {1}{r}}{\frac {\partial }{\partial r}}(rF_{r})+{\frac {1}{r}}{\frac {\partial F_{\theta }}{\partial \theta }}+{\frac {\partial F_{z}}{\partial z}}}
と書ける[2]
球座標系
編集
球面座標系 において、天頂角を θ , z -軸周りの回転角を ϕ とすれば、発散は
div
F
=
∇
⋅
F
=
1
r
2
∂
∂
r
(
r
2
F
r
)
+
1
r
sin
θ
∂
∂
θ
(
sin
θ
F
θ
)
+
1
r
sin
θ
∂
F
ϕ
∂
ϕ
{\displaystyle \operatorname {div} {\boldsymbol {F}}=\nabla \cdot {\boldsymbol {F}}={\frac {1}{r^{2}}}{\frac {\partial }{\partial r}}{\bigl (}r^{2}F_{r}{\bigr )}+{\frac {1}{r\sin \theta }}{\frac {\partial }{\partial \theta }}(\sin \theta \,F_{\theta })+{\frac {1}{r\sin \theta }}{\frac {\partial F_{\phi }}{\partial \phi }}}
と書ける[3]
分解定理
編集
R 3 v r r → ∞v r 無回転成分 (irrotational part ) E r 無発散成分 (source-free part ) B r 回転 の項を参照）から明示的に決定される。即ち、無回転成分は
E
=
−
∇
Φ
(
r
)
,
Φ
(
r
)
=
∫
R
3
d
3
r
′
div
v
(
r
′
)
4
π
|
r
−
r
′
|
{\displaystyle {\boldsymbol {E}}=-\nabla \Phi ({\boldsymbol {r}}),\quad \Phi ({\boldsymbol {r}})=\int _{\mathbb {R} ^{3}}d^{3}{\boldsymbol {r}}'\,{\frac {\operatorname {div} {\boldsymbol {v}}({\boldsymbol {r}}')}{4\pi |{\boldsymbol {r}}-{\boldsymbol {r}}'|}}}
で与えられ、無発散成分 B スカラーポテンシャル Φ(r をベクトルポテンシャル A r −∇Φ を ∇ × A で、湧出密度 div v を循環密度 ∇ × v でそれぞれ置き換えた、
B
=
∇
×
A
(
r
)
,
A
(
r
)
=
∫
R
3
d
3
r
′
∇
×
v
(
r
′
)
4
π
|
r
−
r
′
|
{\displaystyle {\boldsymbol {B}}=\nabla \times {\boldsymbol {A}}({\boldsymbol {r}}),\quad {\boldsymbol {A}}({\boldsymbol {r}})=\int _{\mathbb {R} ^{3}}d^{3}{\boldsymbol {r}}'\,{\frac {\nabla \times {\boldsymbol {v}}({\boldsymbol {r}}')}{4\pi |{\boldsymbol {r}}-{\boldsymbol {r}}'|}}}
で与えられる。
この「分解定理」は電気力学 でも定常流に関する研究の副産物として得られた事実であり、三次元以外でも通用するもっと一般のヘルムホルツ分解 の特別の場合である。
以下の性質は、通常の微分積分学 における常微分の微分法則から導かれる。最も重要なことは、発散作用素が線型作用素 となること、つまり、
div
(
a
F
+
b
G
)
=
a
div
(
F
)
+
b
div
(
G
)
{\displaystyle \operatorname {div} (a{\boldsymbol {F}}+b{\boldsymbol {G}})=a\operatorname {div} ({\boldsymbol {F}})+b\operatorname {div} ({\boldsymbol {G}})}
が任意のベクトル場 F G 実数 a , b
積の微分法則 は以下の形で成立する。φ はスカラー場、F
div
(
φ
F
)
=
(
grad
φ
)
⋅
F
+
φ
(
div
F
)
{\displaystyle \operatorname {div} (\varphi {\boldsymbol {F}})=(\operatorname {grad} \varphi )\cdot {\boldsymbol {F}}+\varphi (\operatorname {div} {\boldsymbol {F}})}
∇ を用いた記法であれば、
∇
⋅
(
φ
F
)
=
(
∇
φ
)
⋅
F
+
φ
(
∇
⋅
F
)
{\displaystyle \nabla \cdot (\varphi {\boldsymbol {F}})=(\nabla \varphi )\cdot {\boldsymbol {F}}+\varphi (\nabla \cdot {\boldsymbol {F}})}
が成り立つ。二つの三次元ベクトル場 F G 交叉積 に対するもう一つの積の法則は、回転 curl を含む以下の形
div
(
F
×
G
)
=
(
curl
F
)
⋅
G
−
F
⋅
(
curl
G
)
{\displaystyle \operatorname {div} ({\boldsymbol {F}}\times {\boldsymbol {G}})=(\operatorname {curl} {\boldsymbol {F}})\cdot {\boldsymbol {G}}-{\boldsymbol {F}}\cdot (\operatorname {curl} {\boldsymbol {G}})}
∇ を用いた記法であれば、
∇
⋅
(
F
×
G
)
=
(
∇
×
F
)
⋅
G
−
F
⋅
(
∇
×
G
)
{\displaystyle \nabla \cdot ({\boldsymbol {F}}\times {\boldsymbol {G}})=(\nabla \times {\boldsymbol {F}})\cdot {\boldsymbol {G}}-{\boldsymbol {F}}\cdot (\nabla \times {\boldsymbol {G}})}
で書くことができる。
スカラー場 φ にラプラス作用素 を施したものは、φ の勾配の発散に等しい。即ち、
div
(
grad
φ
)
=
Δ
φ
{\displaystyle \operatorname {div} (\operatorname {grad} \varphi )=\Delta \varphi }
が成立する。任意の三次元ベクトル場の回転の発散は常に 0 に等しい。即ち、
∇
⋅
(
∇
×
F
)
=
0
{\displaystyle \nabla \cdot (\nabla \times {\boldsymbol {F}})=0}
が成り立つ。発散が 0 のベクトル場 F R 3 G F G R 3 G ポワンカレの補題 の項参照）。この主張が真でなくなる度合は、考える領域 U の複雑性を量化するのに適当な、鎖複体
{
scalar fields on
U
}
→
grad
{
vector fields on
U
}
→
curl
{
vector fields on
U
}
→
div
{
scalar fields on
U
}
{\displaystyle \{{\text{scalar fields on }}U\}\xrightarrow {\operatorname {grad} } \{{\text{vector fields on }}U\}\xrightarrow {\operatorname {curl} } \{{\text{vector fields on }}U\}\xrightarrow {\operatorname {div} } \{{\text{scalar fields on }}U\}}
のホモロジー によって測ることができる。こうしたことが、ドラームコホモロジー の起源および主な動機付けであった。
外微分との関係
編集
発散を外微分の特定の場合として表すことができて、これは R 3
j
=
F
1
d
y
∧
d
z
+
F
2
d
z
∧
d
x
+
F
3
d
x
∧
d
y
{\displaystyle j=F_{1}\,dy\wedge dz+F_{2}\,dz\wedge dx+F_{3}\,dx\wedge dy}
とする。これは局所速度 F (“stuff fluid” of density)
ρ
=
1
d
x
∧
d
y
∧
d
z
{\displaystyle \rho =1\,dx\wedge dy\wedge dz}
j の外微分 dj は
d
j
=
(
∂
F
1
∂
x
+
∂
F
2
∂
y
+
∂
F
3
∂
z
)
d
x
∧
d
y
∧
d
z
=
(
∇
⋅
F
)
ρ
{\displaystyle dj=\left({\frac {\partial F_{1}}{\partial x}}+{\frac {\partial F_{2}}{\partial y}}+{\frac {\partial F_{3}}{\partial z}}\right)dx\wedge dy\wedge dz=(\nabla \cdot {\boldsymbol {F}})\rho }
で与えられる。従ってベクトル場 F
∇
⋅
F
=
⋆
d
⋆
F
♭
{\displaystyle \nabla \cdot {\boldsymbol {F}}={\star }d{\star }{\boldsymbol {F}}^{\flat }}
と表すことができる。ここで上付きの ♭ は下げ同型 （英語版 ） ⋆ はホッジスター である。しかし、外微分は（曲線）座標系の変換と可換だが発散はそうではないので、流れ 2-形式自体を外微分とともに扱うほうが、ベクトル場と発散を扱うよりも容易であることに注意。
ベクトル場の発散の概念を、任意有限次元において一般に定義することができる。
F
=
(
F
1
,
F
2
,
…
,
F
n
)
{\displaystyle {\boldsymbol {F}}=(F_{1},F_{2},\dotsc ,F_{n})}
が R n
x
=
(
x
1
,
x
2
,
…
,
x
n
)
{\displaystyle {\boldsymbol {x}}=(x_{1},x_{2},\dotsc ,x_{n})}
d
x
=
(
d
x
1
,
d
x
2
,
…
,
d
x
n
)
{\displaystyle d{\boldsymbol {x}}=(dx_{1},dx_{2},\dotsc ,dx_{n})}
div
F
=
∇
⋅
F
=
∂
F
1
∂
x
1
+
∂
F
2
∂
x
2
+
⋯
+
∂
F
n
∂
x
n
{\displaystyle \operatorname {div} {\boldsymbol {F}}=\nabla \cdot {\boldsymbol {F}}={\frac {\partial F_{1}}{\partial x_{1}}}+{\frac {\partial F_{2}}{\partial x_{2}}}+\dotsb +{\frac {\partial F_{n}}{\partial x_{n}}}}
で与えられる。より複雑な曲線座標系 においても、適当な表示が得られる。
任意の自然数 n に対して、発散は線型作用素であり、積の法則
∇
⋅
(
φ
F
)
=
(
∇
φ
)
⋅
F
+
φ
(
∇
⋅
F
)
{\displaystyle \nabla \cdot (\varphi {\boldsymbol {F}})=(\nabla \varphi )\cdot {\boldsymbol {F}}+\varphi (\nabla \cdot {\boldsymbol {F}})}
を任意のスカラー場 φ に対して満足する。
あるいは発散の概念を、次元 n で体積要素 （あるいは密度 ）μ を持つ多様体、例えばリーマン多様体 やローレンツ多様体 に対して定義することもできる。R 3 X は (n − 1) -形式 j = iX μ X と μ との縮約によって定める。このとき、ベクトル場 X の発散は等式
d
j
=
div
(
X
)
μ
{\displaystyle dj=\operatorname {div} (X)\mu }
によって定められる。リー微分 に対する標準公式を用いれば、これを
L
X
μ
=
div
(
X
)
μ
{\displaystyle {\mathcal {L}}_{X}\mu =\operatorname {div} (X)\mu }
と書くこともできる。即ち、発散はベクトル場に沿って流すときの体積要素の膨張率を測るものである。
リーマン多様体やローレンツ多様体上で、計量体積要素に関する発散は、レヴィ・チヴィタ接続 ∇ を用いて計算することができて、
div
(
X
)
=
∇
⋅
X
=
X
a
a
{\displaystyle \operatorname {div} (X)=\nabla \cdot X={X^{a}}_{a}}
が成り立つ。真ん中は 1-形式値のベクトル場 ∇X とそれ自身との縮約であり、一番右は物理学者が使う従来の座標表示である（アインシュタインの和の規約 ）。
発散をテンソル に対しても一般化することができる。アインシュタインの和の規約に従って、反変ベクトル Fμ の発散は、
∇
⋅
F
=
∇
μ
F
μ
{\displaystyle \nabla \cdot {\boldsymbol {F}}=\nabla _{\mu }F^{\mu }}
で与えられる。ここで ∇μ は共変微分 である。同じことだが、文献によっては任意の混合テンソル の発散を、添字を上げる同型 “♯ ” を用いて
T が (p , q ) -テンソル（p -階反変かつ q -階共変）ならば、T の発散は
(
div
T
)
(
Y
1
,
…
,
Y
q
−
1
)
=
trace
(
X
↦
♯
(
∇
T
)
(
X
,
⋅
,
Y
1
,
…
,
Y
q
−
1
)
)
{\displaystyle (\operatorname {div} T)(Y_{1},\dotsc ,Y_{q-1})=\operatorname {trace} {\bigl (}X\mapsto \sharp (\nabla T)(X,\cdot ,Y_{1},\dotsc ,Y_{q-1}){\bigr )}}
(p , q − 1) -テンソル、即ち最初の二つの共変添字上の共変微分のトレースである。と定めるものもある。
参考文献
編集
Brewer, Jess H. (1999年4月7日). “DIVERGENCE of a Vector Field ”. Vector Calculus . 2007年9月28日 閲覧。 Theresa M. Korn; Korn, Granino Arthur. Mathematical Handbook for Scientists and Engineers: Definitions, Theorems, and Formulas for Reference and Review . New York: Dover Publications. pp. 157–160. ISBN 0-486-41147-8 関連項目
編集
外部リンク
編集
![{\displaystyle {\overrightarrow {\operatorname {div} ({\boldsymbol {\varepsilon }})}}={\begin{pmatrix}{\dfrac {\partial \varepsilon _{xx}}{\partial x}}+{\dfrac {\partial \varepsilon _{xy}}{\partial y}}+{\dfrac {\partial \varepsilon _{xz}}{\partial z}}\\[5pt]{\dfrac {\partial \varepsilon _{yx}}{\partial x}}+{\dfrac {\partial \varepsilon _{yy}}{\partial y}}+{\dfrac {\partial \varepsilon _{yz}}{\partial z}}\\[5pt]{\dfrac {\partial \varepsilon _{zx}}{\partial x}}+{\dfrac {\partial \varepsilon _{zy}}{\partial y}}+{\dfrac {\partial \varepsilon _{zz}}{\partial z}}\end{pmatrix}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf51f67128a160cd158588c54fc7b57640f70813)
