熱核

提供: miniwiki
移動先:案内検索


数学の特に熱伝導拡散の研究に現れる熱核(ねつかく、: heat kernel)とは、ある適切な境界条件を課された特定の領域上での熱方程式Heat equation)に対する基本解である。ラプラス作用素スペクトルの研究においても重要な道具の一つであり、したがって数理物理学の分野を通して有用な概念である。熱核は、境界がある特定の温度(通常はゼロ)に固定された領域内のある点に単位熱源が時間 t = 0 に置かれた際の、その領域全体での温度変化を表現するものである。

次のような d-次元ユークリッド空間 Rd の熱核が、最も有名である。

[math]K(t,x,y) = \frac{1}{(4\pi t)^{d/2}} e^{-|x-y|^2/4t}. \,[/math]

これは、すべての t > 0 および x,y ∈ Rd に対して、次の熱方程式

[math]\frac{\partial K}{\partial t}(t,x,y) = \Delta_x K(t,x,y)\,[/math]

の解となる。ただし初期条件は

[math]\lim_{t \to 0} K(t,x,y) = \delta(x-y)=\delta_x(y)[/math]

で与えられる。ここで δ はディラックのデルタ関数で、この極限はシュワルツの超函数の意味での極限である。すなわち、コンパクトなを持つなめらかなすべての函数 φ に対して

[math]\lim_{t \to 0}\int_{\mathbf{R}^d} K(t,x,y)\phi(y)\,dy = \phi(x) [/math]

が成り立つ。

Rd 内のより一般の領域 Ω について、上記のような陽的な定式化は一般に可能なものではない。次に述べる円板あるいは正方形領域の簡単な場合はそれぞれ、ベッセル函数テータ函数を含むものである。しかし、熱核は(ディリクレ問題に対しては)任意の領域上で t > 0 に対して依然として存在し、滑らかである。実際、境界が十分正則であるような任意のリーマン多様体上で、そのような熱核は存在する。より正確に言うと、そのような一般の領域におけるディリクレ問題に対する熱核は、次の初期値境界値問題の解として与えられる。

[math]\frac{\partial K}{\partial t}(t,x,y) = \Delta K(t,x,y) \rm{\ \ for\ all\ } t\gt 0 \rm{\ and\ } x,y\in\Omega[/math]
[math]\lim_{t \to 0} K(t,x,y) = \delta_x(y)\rm{\ \ for\ all\ } x,y\in\Omega[/math]
[math]K(t,x,y) = 0, \quad x\in\partial\Omega \rm{\ or\ } y\in\partial\Omega.[/math]

このような場合に、任意の領域上での熱核の正式な表現を導出することは、困難ではない。実際、ある連結領域(あるいは境界を伴う多様体)U でのディリクレ問題を考える。そのラプラシアンに対するディリクレ問題の固有値を λn とする。すなわち、

[math]\left\{ \begin{array}{ll} \Delta \phi + \lambda \phi = 0 & \mathrm{in\ }\ U\\ \phi=0 & \mathrm{on\ }\ \partial U \end{array}\right. [/math]

が成立する。ここで φn はそれらの固有値に対応する固有函数であり、L2(U) において正規直交化されている。このとき、逆ディリクレラプラシアン Δ−1コンパクトかつ自己共役であり、したがってスペクトル定理からそのような固有値は

[math]0 \lt \lambda_1 \lt \lambda_2\le \lambda_3\le\cdots,\quad \lambda_n\to\infty [/math]

を満たすものとなる。このとき、熱核は次のように表現される: テンプレート:NumBlk この級数を和の符号の下で形式的に微分すれば、熱方程式を満たすものであることが確かめられる。しかし、この級数の収束や正則性といった問題は、また慎重に考える必要がある。

熱核はまた、しばしば対応する積分変換と関連付けて考えられる。そのような積分変換は、コンパクトな台を持つなめらかな φ に対して

[math]T\phi = \int_\Omega K(t,x,y)\phi(y)\,dy [/math]

のように定義される。スペクトル写像定理によって、次のような T の表現を得ることが出来る。

[math]T = e^{t\Delta}.[/math]

関連項目

参考文献