グリーンの恒等式のソースを表示


[[数学]]において'''グリーンの恒等式'''（グリーンのこうとうしき、{{Lang-en-short|Green's identities}}）とは、[[ベクトル解析]]に現れる三つの恒等式のことを言う。[[グリーンの定理]]を発見した数学者の[[ジョージ・グリーン]]の名にちなむ。

== グリーンの第一恒等式 ==

この恒等式は、ベクトル場 {{math|'''F''' {{=}} ''ψ''∇''φ''}} に対して[[発散定理]]と恒等式:

: <math>\nabla \cdot \psi \nabla \varphi = \psi \Delta \varphi + \nabla \psi \cdot \nabla \varphi</math>

を適用することで次のように得られる：{{mvar|φ}} と {{mvar|ψ}} をある領域 {{math|''U'' ⊂ '''R'''<sup>3</sup>}} 上で定義されるスカラー函数とし、{{mvar|φ}} は二階連続的微分可能、{{mvar|ψ}} は一階連続的微分可能とする。このとき次が成り立つ<ref name="strauss">{{cite book|last=Strauss|first=Walter|title=Partial Differential Equations: An Introduction|publisher=Wiley}}</ref>。

: <math>\int_U \left( \psi \Delta \varphi + \nabla \varphi \cdot \nabla \psi \right)\, dV  = \oint_{\partial U} \psi \left( \nabla \varphi \cdot \mathbf{n} \right)\, dS=\oint_{\partial U}\psi\nabla\varphi\cdot d\mathbf{S} </math>

ここで <math>\Delta</math> は[[ラプラス作用素]]、{{math|∂''U''}} は領域 {{mvar|U}} の境界、{{math|'''n'''}} は面素 {{math|''dS''}} に対する外向き法線ベクトル、{{mvar|d}}{{math|'''S'''}} は向き付けられた面素である。この定理は[[発散定理]]の特別な場合であり、{{mvar|ψ}} と {{mvar|φ}} の勾配をそれぞれ {{mvar|u}} と {{mvar|v}} で置き換えた[[部分積分]]の高次元版と本質的に同値である。

上述のグリーンの第一恒等式は、[[発散定理]]において {{math|'''F''' {{=}} ''ψ'''''Γ'''}} とすることで得られる次のより一般の恒等式の特別な場合である：

: <math>\int_U \left( \psi \nabla \cdot \mathbf{\Gamma} + \mathbf{\Gamma} \cdot \nabla \psi\right)\, dV  = \oint_{\partial U} \psi \left( \mathbf{\Gamma} \cdot \mathbf{n} \right)\, dS=\oint_{\partial U}\psi\mathbf{\Gamma}\cdot d\mathbf{S}. </math>

== グリーンの第二恒等式 ==

{{mvar|φ}} と {{mvar|ψ}} のいずれもが {{math|''U'' ⊂ '''R'''<sup>3</sup>}} 上で二階連続的微分可能であり、{{mvar|ε}} が一階連続的微分可能であるなら、{{math|'''F''' {{=}} ''ψ(ε''∇''φ'')}}と{{math|'''F''' {{=}} ''ψ(ε''∇''φ'')}} のそれぞれについて上の恒等式を立て, それらを辺々引くことで次が得られる：

: <math> \int_U \left[ \psi \nabla \cdot \left( \epsilon \nabla \varphi \right) - \varphi \nabla \cdot \left( \epsilon \nabla \psi \right) \right]\, dV = \oint_{\partial U} \epsilon \left( \psi {\partial \varphi \over \partial \mathbf{n}} - \varphi {\partial \psi \over \partial \mathbf{n}}\right)\, dS. </math>

{{math|''U'' ⊂ '''R'''<sup>3</sup>}} 上すべてで {{math|''ε'' {{=}} 1}} であるような特別な場合は、次が得られる：

: <math> \int_U \left( \psi \Delta \varphi - \varphi \Delta \psi\right)\, dV = \oint_{\partial U} \left( \psi {\partial \varphi \over \partial \mathbf{n}} - \varphi {\partial \psi \over \partial \mathbf{n}}\right)\, dS. </math>

上式において {{math|∂''φ''/∂'''n'''}} は {{mvar|φ}} の、面素 {{math|''dS''}} に対する外向き法線ベクトル {{math|'''n'''}} の方向での[[方向微分]]である。すなわち

: <math> {\partial \varphi \over \partial \mathbf{n}} = \nabla \varphi \cdot \mathbf{n}=\nabla_\mathbf{n}\varphi</math>

である。このことは特に、境界上で消失する函数に対する L2 内積において、ラプラシアンは[[自己共役作用素|自己共役]]であることを意味する。

== グリーンの第三恒等式 ==

グリーンの第三恒等式は、第二恒等式において {{math|''φ'' {{=}} ''G''}} とすることで得られる。ここで {{mvar|G}} は[[ラプラス作用素]]の[[基本解]]として定められる[[グリーン函数]]である。すなわち、次が成り立つ：

:<math> \Delta G(\mathbf{x},\mathbf{\eta}) = \delta(\mathbf{x} - \mathbf{\eta}).</math>

例えば {{math|'''R'''<sup>3</sup>}} において、ある解は次の形を取る：

:<math>G(\mathbf{x},\mathbf{\eta})= \frac{-1}{4 \pi \|\mathbf{x} - \mathbf{\eta} \|}.</math>

グリーンの第三恒等式は、{{mvar|ψ}} が {{mvar|U}} 上で二階連続的微分可能な函数である場合に、次式で与えられる：

: <math> \int_U \left[ G(\mathbf{y},\mathbf{\eta}) \Delta \psi(\mathbf{y})\right]\, dV_\mathbf{y} - \psi(\mathbf{\eta})=  \oint_{\partial U} \left[ G(\mathbf{y},\mathbf{\eta}) {\partial \psi \over \partial \mathbf{n}} (\mathbf{y}) - \psi(\mathbf{y}) {\partial G(\mathbf{y},\mathbf{\eta}) \over \partial \mathbf{n}} \right]\, dS_\mathbf{y}.</math>

{{mvar|ψ}} それ自身が[[調和函数]]、すなわち[[ラプラス方程式]]の解であるなら、{{math|∇<sup>2</sup>''ψ'' {{=}} 0}} よりこの式は次のように簡易化される：

:<math>\psi(\mathbf{\eta})= \oint_{\partial U} \left[\psi(\mathbf{y}) \frac{\partial G(\mathbf{y},\mathbf{\eta})}{\partial \mathbf{n}} - G(\mathbf{y},\mathbf{\eta}) \frac{\partial \psi}{\partial \mathbf{n}} (\mathbf{y}) \right]\, dS_\mathbf{y}.</math>

この積分の中の第二項は、問題が適切であるような領域 {{mvar|U}} の境界に対する[[グリーン函数]]として {{mvar|G}} が選ばれる（[[ディリクレ境界条件]]）とき、消える。すなわち

:<math>\psi(\mathbf{\eta}) = \oint_{\partial U} \psi(\mathbf{y}) \frac{\partial G(\mathbf{y},\mathbf{\eta})}{\partial \mathbf{n}} \, dS_\mathbf{y}</math>

となる。この形式はディリクレ境界条件の問題に対する解を構成する際に用いられる。[[ノイマン境界条件]]の問題に対する解を見つけるためには、境界上で法線微分が消失するようなグリーン函数を代わりに用いればよい。

さらに上述の恒等式は、 {{mvar|ψ}} が[[ヘルムホルツ方程式]]あるいは[[波動方程式]]の解で、{{mvar|G}} が適切なグリーン函数の解である場合にも適用される。そのような文脈において、この恒等式は[[ホイヘンス＝フレネルの原理]]の数学的表現となる。

== 多様体上での恒等式 ==

グリーンの恒等式はリーマン多様体上で成立する。この場合、初めの二つの恒等式は次のようになる。

:<math>\begin{align}
\int_M u\Delta v\, dV + \int_M \langle\operatorname{grad}\ u, \operatorname{grad}\ v\rangle\, dV &= \int_{\partial M} u N v d\widetilde{V} \\
\int_M \left (u \Delta v - v \Delta u \right )\, dV &= \int_{\partial M}(u N v - v N u) d \widetilde{V}.
\end{align}</math>

ここで {{mvar|u}} と {{mvar|v}} は {{mvar|M}} 上の滑らかな実数値函数、{{mvar|dV}} は計量と両立する[[体積形式]]、<math>d\widetilde{V}</math> は {{mvar|M}} の境界上で誘起された体積形式、{{mvar|N}} は境界上での向き付けられた単位法ベクトル場、{{math|∇<sup>2</sup>''u'' {{=}} div(grad ''u'')}} はラプラシアンである。

== グリーンのベクトル恒等式 ==

グリーンの第二恒等式は、二つのスカラー函数の二階微分と一階微分（の発散）との関係を表すものである。微分形式では、

:<math>p_{m}\Delta q_{m}-q_{m}\Delta p_{m}=\nabla\cdot\left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right)</math>

となる。ここで {{math|''p<sub>m</sub>''}} と {{math|''q<sub>m</sub>''}} は二つの任意の二階連続的微分可能なスカラー場である。この恒等式は、質量やエネルギーのようなスカラー場に対して連続の方程式が成立する担保となるため、物理学において非常に重要なものとなっている<ref>M. Fernández-Guasti. Complementary fields conservation equation derived from the scalar wave equation. ''J. Phys. A: Math. Gen.'', 37:4107–4121, 2004.</ref>。グリーンの第二恒等式はベクトル解析においてよく示されているが、多くの教科書ではスカラーの場合のみが扱われている。専門書においても、ベクトルの場合を見つけるのは容易ではない。ベクトル回折理論において、そのようなグリーンの第二恒等式の二つの場合が導入されている。一方は、[[クロス積]]の発散によるもの<ref>A. E. H. Love. The Integration of the Equations of Propagation of Electric Waves. ''Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character'', 197:pp. 1–45, 1901.</ref><ref>J. A. Stratton and L. J. Chu. Diffraction Theory of Electromagnetic Waves. ''Phys. Rev.'', 56(1):99–107, Jul 1939.</ref><ref>N. C. Bruce. Double scatter vector-wave Kirchhoff scattering from perfectly conducting surfaces with infinite slopes. ''Journal of Optics'', 12(8):085701, 2010.</ref>で、場における回転と回転に関する関係を表している：

:<math>\mathbf{P}\cdot\left(\nabla\times\nabla\times\mathbf{Q}\right)-\mathbf{Q}\cdot\left(\nabla\times\nabla\times\mathbf{P}\right)=\nabla\cdot\left(\mathbf{Q}\times\nabla\times\mathbf{P}-\mathbf{P}\times\nabla\times\mathbf{Q}\right).</math>

この方程式はラプラシアンを使って次のように書くことが出来る：

:<math>\mathbf{P}\cdot\Delta \mathbf{Q}-\mathbf{Q}\cdot\Delta \mathbf{P}+\mathbf{Q}\cdot\left[\nabla\left(\nabla\cdot\mathbf{P}\right)\right]-\mathbf{P}\cdot\left[\nabla\left(\nabla\cdot\mathbf{Q}\right)\right]=\nabla\cdot\left(\mathbf{P}\times\nabla\times\mathbf{Q}-\mathbf{Q}\times\nabla\times\mathbf{P}\right).</math>

しかし、項

:<math>\mathbf{Q}\cdot\left[\nabla\left(\nabla\cdot\mathbf{P}\right)\right]-\mathbf{P}\cdot\left[\nabla\left(\nabla\cdot\mathbf{Q}\right)\right],</math>

は発散によって書くことは出来ない。もう一方の場合は、二重ベクトルを導入する手法で、二価グリーン函数が必要となる<ref>W. Franz, On the Theory of Diffraction. ''Proceedings of the Physical Society. Section A'', 63(9):925, 1950.</ref><ref>Chen-To Tai. Kirchhoff theory: Scalar, vector, or dyadic? ''Antennas and Propagation, IEEE Transactions on'', 20(1):114–115, jan 1972.</ref>。ここで示すものはそれらの問題を避けるものである<ref>M. Fernández-Guasti. Green's second identity for vector fields. ISRN Mathematical Physics, 2012:7, 2012. Article ID: 973968. [http://www.isrn.com/journals/mp/2012/973968]</ref>。

グリーンの第二恒等式におけるスカラー場は、ベクトル場のデカルト成分であると考える。すなわち

:<math>\mathbf{P}=\sum_m p_{m}\hat{\mathbf{e}}_{m}, \qquad \mathbf{Q}=\sum_m q_{m}\hat{\mathbf{e}}_{m}</math>

とする。各成分に対する方程式を足し上げることで、次が得られる。

:<math>\sum_m \left[p_m\Delta q_m-q_m\Delta p_m\right]=\sum_m \nabla\cdot\left(p_m\nabla q_m-q_m\nabla p_m \right).</math>

この左辺は、ドット積の定義により、次のようなベクトル形式で書くことが出来る。

:<math>\sum_m \left[p_{m}\Delta q_{m}-q_{m}\Delta p_{m}\right]=\mathbf{P}\cdot\Delta\mathbf{Q}-\mathbf{Q}\cdot\Delta\mathbf{P}.</math>

右辺は、ベクトル作用素によって表現するにはやや複雑なものとなっている。発散作用素の加法についての分配性より、発散の和は和の発散に等しい。すなわち、

:<math>\sum_m \nabla\cdot\left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right)= \nabla\cdot\left(\sum_m p_{m}\nabla q_{m}-\sum_m q_{m}\nabla p_{m}\right) </math>

である。ドット積の勾配に対する次のベクトル恒等式を考える。

:<math>\nabla\left(\mathbf{P}\cdot\mathbf{Q}\right)=\left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}+\left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}+\mathbf{P}\times\nabla\times\mathbf{Q}+\mathbf{Q}\times\nabla\times\mathbf{P},</math>

これは、ベクトル成分について書くと次のようになる。

:<math>\nabla\left(\mathbf{P}\cdot\mathbf{Q}\right)=\nabla\sum_m p_{m}q_{m}=\sum_m p_{m}\nabla q_{m}+\sum_m q_{m}\nabla p_{m}.</math> 

この結果は、負の符号を除いて、ベクトル形式で表すことを所望していたものと似たものである。各項の微分作用素は、あるベクトル（<math>p_{m}</math>）あるいはもう一方（<math>q_{m}</math>）のいずれかに対して作用するため、各項は次のように表される：

:<math>\sum_m p_m \nabla q_m = \left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}+\mathbf{P}\times\nabla\times\mathbf{Q},</math>
:<math>\sum_m q_m \nabla p_m = \left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}+\mathbf{Q}\times\nabla\times\mathbf{P}.</math>

これらの結果は、[http://luz.izt.uam.mx/mediawiki/index.php/Green%27s_vector_identity ベクトル成分の評価]を介して厳密に証明することも出来る。以上より、右辺は次のようなベクトル形式で表される。

:<math>\sum_m p_{m}\nabla q_{m}-\sum_m q_{m}\nabla p_{m}=\left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}+\mathbf{P}\times\nabla\times\mathbf{Q}-\left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}-\mathbf{Q}\times\nabla\times\mathbf{P}.</math>

これら二つの結果を合わせることで、スカラー場に対するグリーンの定理と同様の結果が得られる：

:'''ベクトル場に対する定理''' 
::<math>\color{OliveGreen}\mathbf{P} \cdot \Delta \mathbf{Q} - \mathbf{Q} \cdot \Delta \mathbf{P} = \nabla \cdot \left[ \left( \mathbf{P}\cdot \nabla \right) \mathbf{Q} + \mathbf{P} \times \nabla \times \mathbf{Q}-\left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}-\mathbf{Q}\times\nabla\times\mathbf{P}\right].</math>

クロス積の回転は次のように書くことが出来る。

:<math>\nabla\times\left(\mathbf{P}\times\mathbf{Q}\right)=\left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}-\left(\mathbf{P}\cdot\nabla\right)\mathbf{Q}+\mathbf{P}\left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla\cdot\mathbf{P}\right);</math> 

グリーンのベクトル恒等式は、次のように書くことが出来る：

:<math>\mathbf{P}\cdot\Delta \mathbf{Q}-\mathbf{Q}\cdot\Delta \mathbf{P}= \nabla\cdot \left[\mathbf{P} \left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla\cdot\mathbf{P}\right)-\nabla\times\left(\mathbf{P}\times\mathbf{Q}\right)+\mathbf{P}\times\nabla\times\mathbf{Q}-\mathbf{Q}\times\nabla\times\mathbf{P}\right].</math>

回転の発散はゼロであるため、この第三項は消去され、次が得られる：

:'''グリーンのベクトル恒等式''' 
::<math>\color{OliveGreen}\mathbf{P}\cdot\Delta\mathbf{Q}-\mathbf{Q}\cdot\Delta\mathbf{P}=\nabla\cdot\left[\mathbf{P}\left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla\cdot\mathbf{P}\right)+\mathbf{P}\times\nabla\times\mathbf{Q}-\mathbf{Q}\times\nabla\times\mathbf{P}\right].</math>

同様の手順で、ドット積のラプラシアンはラプラシアンに関して次のように表すことが出来る：

:<math>\Delta\left(\mathbf{P}\cdot\mathbf{Q}\right)=\mathbf{P}\cdot\Delta\mathbf{Q}-\mathbf{Q}\cdot\Delta \mathbf{P} +2\nabla\cdot\left[\left(\mathbf{Q}\cdot\nabla\right)\mathbf{P}+\mathbf{Q}\times\nabla\times\mathbf{P}\right].</math>

この系として、整理されていない項はベクトルグリーン函数との比較によって発散に関して書くことが出来る：

:<math>\mathbf{P}\cdot \left[ \nabla \left(\nabla \cdot \mathbf{Q} \right) \right]-\mathbf{Q} \cdot \left[ \nabla \left( \nabla \cdot \mathbf{P} \right) \right] = \nabla\cdot\left[\mathbf{P}\left(\nabla\cdot\mathbf{Q}\right)-\mathbf{Q}\left(\nabla \cdot \mathbf{P} \right) \right].</math>

この結果は、右辺のベクトルのスカラー倍の発散を拡張することで確かめることが出来る。

== 関連項目 ==
* [[グリーン函数]]
* {{仮リンク|キルヒホッフの積分定理|en|Kirchhoff integral theorem}}

== 参考文献 ==
{{reflist}}

== 外部リンク ==
* {{SpringerEOM|title=Green formulas|urlname=Green_formulas}}
*[http://mathworld.wolfram.com/GreensIdentities.html] Green's Identities at Wolfram MathWorld

{{DEFAULTSORT:くりいんのこうとうしき}}
[[Category:ベクトル解析]]
[[Category:恒等式]]
[[Category:数学に関する記事]]
[[Category:数学のエポニム]]