ラプラス場

ベクトル場 \bm{A} が, \nabla \cdot \bm{A}= 0 を満たすとき管状ベクトル場, \nabla \times \bm{A}= \bm{0} を満たすとき層状ベクトル場と呼ぶのでした.この両方を満たすようなベクトル場,すなわち管状でもあり層状でもあるベクトル場を ラプラス場 と呼びます.

\nabla \cdot \bm{A} = 0, \ \ \ \nabla \times \bm{A} = \bm{0}   \tag{1}

ラプラス場が単連結領域で定義されている場合,この流れは ラプラス方程式 を満たすスカラーポテンシャルによって記述されます.

\bm{A} = \nabla \phi ,  \ \ \    s.t. \ \ \triangle \phi  = 0  \tag{2}

(2) で表現される \bm{A} が,式 (1) の条件を両方満たすことを確認してみて下さい.ラプラス方程式の解となる関数を,一般に 調和関数 と呼ぶのでした.調和関数には,以下に紹介するような重要な諸性質があります.ラプラス場の解析では,これらの性質が効いて来ます.

[*]以下の諸性質は,ベクトル場の問題というよりは調和関数の性質に関する話題です.物理の様々な場面で,ラプラス方程式を適当な境界条件下で解く問題は重要ですが,これは,ベクトル解析の問題というよりは,微分方程式論もしくはポテンシャル論の守備範囲になると思います.ですから,これらの事項を今すぐ覚える必要はないと思いますが,いずれ使う日が来ると思いますので,ここでざっと見ておいても損はないと思います.この記事では,あまり詳しい所までは取り上げません.

性質1

閉曲面 S 上で,調和関数 \phi の法線方向の方向微分の, S 全面に渡る面積分は常に 0 になります.

\int \int \limits _{S} \frac{\partial \phi}{\partial n}dS = 0  \tag{3}

この関係は, ガウスの発散定理\nabla \phi を代入して次のように変形すれば,左辺 =0 よりすぐに導けます.

\int \int \int \limits _{V} \triangle \phi dV = \int \int \limits _{S} \frac{\partial \phi}{\partial n}dS      \tag{4}

性質2

閉曲面 S 上で,二つの調和関数 \phi , \psi に対し グリーンの定理 \int \int \int \limits _{V} (\phi \triangle \psi - \psi \triangle \phi ) dV = \int \int \limits _{S}\left( \phi \frac{\partial \psi}{\partial n}  - \psi \frac{\partial \phi}{\partial n} \right) dS を適用すれば,左辺 =0 より次式を得ます.

\int \int \limits _{S}\left( \phi \frac{\partial \psi}{\partial n}  - \psi \frac{\partial \phi}{\partial n} \right) dS=0       \tag{5}

これは閉領域で定義された二つの調和関数が満たす関係式として重要です.

性質3

次の公式を グリーンの定理 の最後で求めました.

\phi (M) = -\frac{1}{4\pi} \int \int \int \limits _{V} \frac{1}{r}\triangle \phi dV -\frac{1}{4\pi}\int \int \limits _{S} \left[ \phi \frac{\partial }{\partial n} \left( \frac{1}{r} \right ) - \frac{1}{r}  \frac{\partial }{\partial n} \right] dS  \tag{6}

この右辺は \triangle \phi =0 を代入すると第一項が消えて簡単化します.式 (6) の表式に従う一般にはスカラーポテンシャル \phi では,領域内の点 M における値を決めるのに,領域における体積分と面積分の両方を考えなければならなかったのに対し, \phi が調和関数なら,面積分だけで良くなるわけです.

\phi (M) = -\frac{1}{4\pi}\int \int \limits _{S} \left[ \phi \frac{\partial }{\partial n} \left( \frac{1}{r} \right ) - \frac{1}{r}  \frac{\partial }{\partial n} \right] dS   \tag{7}

また,式 (7) の曲面 S を点 M を中心とする半径 \rho の球面 K_{g} とすると, \int \int \limits _{K_{g}} \frac{1}{r} \frac{\partial \phi}{\partial n}dS = \frac{1}{\rho} \int \int \limits _{K_{g}} \frac{\partial \phi}{\partial n}dK_{g} =0 が成り立ちます.(球は全方位に対称ですから,法線ベクトルを全部足したら 0 になりますね).また, \frac{\partial }{\partial n} \left( \frac{1}{r} \right) _{K_{g}}= \frac{\partial }{\partial n} \left( \frac{1}{r} \right) _{r=\rho }=-\frac{1}{\rho ^{2}} にも注意すると,式 (7) は次のように簡単化されます.

\phi (M) = -\frac{1}{4\pi \rho^{2}}\int \int \limits _{K_{g}} \phi dK_{g}      \tag{8}

球対称なポテンシャルの分布を考える問題は,割合によく物理学に出て来る重要な形ですので,式 (7) の形よりも式 (8) の簡単形で使う方が実用的かも知れません.式 (8) は『点 M における調和関数 \phi の値は, M を中心とする球内(半径によらない!!)の \phi の値の平均値に等しい』という主張だと読めます.

[†]これは著しく美しい調和関数の性質です.例えば,代表的な調和関数である三角関数を思い浮かべれば,周期的に変動する値を,平均値で表わすという意味が直観的に分かりやすいと思います.

性質4

調和関数 \phi が領域 D 内で定数関数ではないとすると, \phiD 内で最大値も最小値も取りません.

proof

例えば,もし D 内の点 M_{1}\phi が最大値を取ったとすると, D をすっぽり包むような M_{1} を中心とする半径 \rho (十分大きい)の球状領域 S を, \phi (M) < \phi (M_{1}) となるように取れます.ただし, MS 内の任意の点とします.一方,式 (8) より, \phi (M_{1}) = -\frac{1}{4\pi \rho^{2}}\int \int \limits _{S} \phi dS < \phi (M_{1}) (中辺は S 内の \phi の平均値の意味)が要請されます.これはおかしいので,結局, \phi は最大値を取らないことになります.同様にして,最小値も取れません.■

性質5【ディリクレ問題の定理】

theorem

領域 V の境界は閉曲面 S だとします.調和関数 \phi がもし S 上で一定値 c を取るなら, \phiV 内の至るところで一定値 c を取るということが言えます.

(式 (6)\phi=c と置けば示せます.)

性質6

theorem

ラプラス方程式 \triangle \phi =0 は,境界 S 上の境界条件下に解くとき,一意的な解を持ちます.

proof

もし \triangle \phi =0 が二つの解 \phi_{1},\phi_{2} を持つとすると,ラプラス方程式の線形性より,その差 \phi_{3} = \phi_{1}-\phi_{2} もやはりラプラス方程式を満たすはずです.ところで,境界 S 上では \phi_{1}\phi_{2} は同じ値を取るはずで, \phi_{3} = \phi_{1}-\phi_{2}=0 となります.すると,性質5より, V 内では至るところ \phi_{3}=0 となりますので,すなわち \phi_{1}=\phi_{2} が示されます.■

[‡]ラプラス方程式 \triangle \phi =0 に,境界条件として境界上の値 \phi (S) を与えて解く問題をディレクレ問題と言います.境界条件として境界上での方向微分 \frac{\partial \phi}{\partial n}_{S} を与えて解く問題をノイマン問題と言います.いずれも物理学で重要な微分方程式ですが,ここではその解法までを示すことはしません.興味のある人は微分方程式論やポテンシャル論といった題名の教科書を当たってみて下さい.

性質7【ノイマン問題の定理】

theorem

ラプラス方程式 \triangle \phi =0 が,境界 S 上で \frac{\partial \phi}{\partial n}=0 を満たすとき, V 内の至るところで \phi = const. となります.

proof

これは, グリーンの第一定理 の系 \int \int \int \limits _{V}( \phi \triangle \phi +( \nabla phi )^{2} )dV = \int \int \limits _{S} \phi \frac{\partial \phi}{\partial n}dS より, \int \int \int \limits _{V} (\nabla \phi)^{2}dV=0 を得ます.これより \nabla \phi =0 となり,積分して \phi = const. を得ます.■

性質8

theorem

ラプラス方程式 \triangle \phi =0 の二つの解 \phi_{1},\phi_{2} の方向微分が,境界 S 上で同一であるとき, \phi_{1}-\phi_{2}=const. となります.

proof

二解の差 \phi_{3}= \phi_{1}-\phi_{2} の方向微分を取って \frac{\partial \phi_{3}}{\partial n}= \frac{\partial \phi_{1}}{\partial n}-\frac{\partial \phi_{2}}{\partial n} となりますが,右辺は条件より 0 となります.これを積分すると, \phi_{1} - \phi_{2}=const. を得ます.■