Laplacian と物理応用
mathvector-calculuslaplacianlecture
1導入
この講義では、Laplacian が局所平均からの偏差を特徴づけ、拡散・波動・ポテンシャルを記述する中心的な演算子であることを説明する。
2用語と定義
\Omega\subset\mathbb R^n を開集合とする。Laplacian は、f\in C^2(\Omega) であるスカラー場に対して
\Delta f=\nabla\cdot\nabla f
で定義される二階線形微分演算子である。\mathbb R^n の直交座標では
\Delta f=\sum_{i=1}^n\frac{\partial^2f}{\partial x_i^2}
となる。以下では Euclidean 空間のスカラー Laplacian を扱う。前講義では Gauss の発散定理から、\Delta u の体積積分が勾配 \nabla u の外向き境界流束に等しいことを確認した。この講義では、その局所的な意味を平均値との関係から解説する。
3方針
Laplacian は、一点の値と、その点を囲む小球面での平均との偏差の主要項を与える。この関係を Taylor 展開で導出してから、熱方程式・波動方程式・Laplace 方程式・Poisson 方程式へ接続する。
4局所平均からの差
f\in C^2(\Omega)、x\in\Omega とし、x を中心とする半径 r の閉球が \Omega に含まれるとする。その境界球面を \partial B_r(x) と書き、球面積で正規化した球面平均を
M_r f(x)=\frac{1}{|\partial B_r(x)|}\int_{\partial B_r(x)}f\,dS
と定義する。r\to0 で Taylor 展開すると、
M_r f(x)=f(x)+\frac{r^2}{2n}\Delta f(x)+o(r^2)
を得る。ここで o(r^2) は、r\to0 で o(r^2)/r^2\to0 を満たす剰余項を表す。一次の項と二階偏導関数行列 Hessian の非対角項は球面上の対称性で平均が 0 になり、対角項の和だけが Laplacian として残る。
したがって \Delta f(x)>0 なら、主要項が剰余項を支配する十分小さい r に対して、周囲の平均は中心の値より高い。\Delta f(x)<0 なら逆である。ただし、一点で \Delta f(x)=0 というだけでは、有限半径の平均と f(x) が等しいとは限らない。
一次元では、Laplacian は二階微分そのものであり、グラフの局所的な凹凸を特徴づける。高次元では、互いに直交する方向の二階変化を合計する。
\Delta f=0 が \overline{B_R(x)} の近傍で成立するとする。球面平均を単位球面 S^{n-1} 上の積分として微分し、Gauss の発散定理を適用すると、0<r<R で
\frac{d}{dr}M_r f(x)
=\frac{1}{|S^{n-1}|r^{n-1}}\int_{\partial B_r(x)}\nabla f\cdot\boldsymbol n\,dS
=\frac{1}{|S^{n-1}|r^{n-1}}\int_{B_r(x)}\Delta f\,dV=0
となる。したがって M_r f(x) は r に依存せず、r\to0 で M_r f(x)\to f(x) であるため、すべての 0<r<R において M_r f(x)=f(x) が成立する。これが調和関数の平均値性質である。
5代表例
以下では、均質・等方・定係数の最も単純な古典モデルと十分滑らかな古典解を考える。熱方程式では空間について C^2、時間について C^1、波動方程式では空間・時間について C^2、Laplace・Poisson 方程式では空間について C^2 を仮定する。
熱拡散率 \kappa>0 を持つ熱方程式は
u_t=\kappa\Delta u
であり、温度が周囲との差を均すように変化することを表す。
Laplace 方程式
\Delta u=0
は、内部に源のない平衡状態を表す。この解が調和関数である。
| 方程式 | Laplacian の役割 |
| 熱方程式 u_t=\kappa\Delta u | 周囲との差を均して温度を拡散させる |
| 波動方程式 u_{tt}=c^2\Delta u | 空間的な張力や曲がりを加速度へ変換する |
| Poisson 方程式 -\Delta u=f | ここではこの符号規約を採用し、源 f からポテンシャルを決定する |
| Laplace 方程式 \Delta u=0 | 源のない平衡を表す |
6調和関数の例
u(x,y)=x^2-y^2 では u_{xx}=2、u_{yy}=-2 であるため、\Delta u=0 である。このような関数を調和関数という。この関数は、その定義域内に閉包が含まれるすべての球で平均値性質を満たす。
別の例として u(x,y)=\log\sqrt{x^2+y^2} は、原点を除く領域で調和である。ただし原点に特異点があるため、領域を指定せずに全平面の調和関数として扱ってはならない。この例は Green 関数や基本解への入口になる。
7座標系による式の変化
二次元の直交座標では \Delta f=f_{xx}+f_{yy} である。r>0 の極座標では
\Delta f=f_{rr}+\frac{1}{r}f_r+\frac{1}{r^2}f_{\theta\theta}
となる。この式は x=r\cos\theta、y=r\sin\theta に連鎖律を適用し、f_{xx}+f_{yy} を整理して導出できる。r=0 では極座標そのものが特異なので、この表示をそのまま代入しない。同じ Euclidean Laplacian でも座標表示は変化し、1/r などの項は座標変換の結果として現れる。
8単位の確認
熱方程式で u の単位を \mathrm{K}、長さの単位を \mathrm{m} とすると、\Delta u の単位は \mathrm{K/m^2} である。\kappa の単位が \mathrm{m^2/s} なら、\kappa\Delta u は \mathrm{K/s} となり、u_t と一致する。
9混同しやすい点
scalar Laplacian はスカラー場へ作用する。\mathbb R^3 の直交座標では、C^2 級ベクトル場 \boldsymbol F の vector Laplacian を各成分へスカラー Laplacian を作用させて定義でき、
\Delta\boldsymbol F=\nabla(\nabla\cdot\boldsymbol F)-\nabla\times(\nabla\times\boldsymbol F)
が成立する。右辺を各成分で展開すると混合偏微分の項が相殺し、\Delta F_i だけが残る。曲線座標で「表示成分を単純に二階微分する」だけでは、基底の変化を落とすため同じ演算子にならない。
10どこまで成立するか
ここで扱ったのは Euclidean 空間の Laplacian と均質・等方な定係数モデルである。不均質または異方的な媒質では、たとえば \nabla\cdot(A(x)\nabla u) のような可変係数演算子が現れる。一般の Riemann 多様体では計量から定まる Laplace--Beltrami 演算子を用いる。また、平均値性質は一点の \Delta f=0 からではなく、球を含む領域での調和性から得られる。
11関連リンク
data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md
data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md
data/lecture/math/vector-calculus/vector-calculus-portal.lecture.n.md
data/lecture/physics/waves/wave-equation-basics.lecture.n.md
The Laplacian and Its Physical Applications
1Introduction
This lecture explains how the Laplacian characterizes deviation from a local average and serves as a central operator in the mathematical descriptions of diffusion, waves, and potentials.
2Terminology and Definition
Let \Omega\subset\mathbb R^n be open. For a scalar field f\in C^2(\Omega), the Laplacian is the second-order linear differential operator
\Delta f=\nabla\cdot\nabla f.
In Cartesian coordinates on \mathbb R^n,
\Delta f=\sum_{i=1}^n\frac{\partial^2f}{\partial x_i^2}.
This lecture concerns the scalar Laplacian in Euclidean space. The preceding lecture used the divergence theorem to show that the volume integral of \Delta u equals the outward boundary flux of \nabla u. Here, the local significance of this relation is explained through local averages.
3Strategy
The Laplacian supplies the leading term in the difference between the value at a point and the average over a small sphere surrounding that point. We first derive this relation from a Taylor expansion and then connect it to the heat, wave, Laplace, and Poisson equations.
4Deviation from a Local Average
Let f\in C^2(\Omega) and x\in\Omega, and suppose that the closed ball of radius r centered at x lies in \Omega. Define the spherical average by
M_r f(x)=\frac{1}{|\partial B_r(x)|}\int_{\partial B_r(x)}f\,dS.
A Taylor expansion as r\to0 gives
M_r f(x)=f(x)+\frac{r^2}{2n}\Delta f(x)+o(r^2).
Here o(r^2) denotes a remainder satisfying o(r^2)/r^2\to0 as r\to0. The linear terms and the off-diagonal entries of the Hessian have zero spherical average by symmetry. The sum of its diagonal entries remains and equals the Laplacian. Consequently, if \Delta f(x)>0, the leading term dominates the remainder and the surrounding average exceeds the central value for every sufficiently small r; if \Delta f(x)<0, the inequality is reversed. The condition \Delta f(x)=0 at one point alone does not imply equality with an average over a ball of finite radius.
In one dimension, the Laplacian is the second derivative and characterizes the local concavity or convexity of a graph. In higher dimensions, it sums second-order changes in mutually orthogonal directions.
Suppose \Delta f=0 on a neighborhood of \overline{B_R(x)}. Differentiating the spherical average after expressing it as an integral over the unit sphere, and then applying the divergence theorem, gives, for 0<r<R,
\frac{d}{dr}M_r f(x)
=\frac{1}{|S^{n-1}|r^{n-1}}\int_{\partial B_r(x)}\nabla f\cdot\boldsymbol n\,dS
=\frac{1}{|S^{n-1}|r^{n-1}}\int_{B_r(x)}\Delta f\,dV=0.
Thus M_r f(x) is independent of r. Since M_r f(x)\to f(x) as r\to0,
f(x)=M_r f(x)
for every 0<r<R. This identity is the mean-value property.
5Principal Equations
We consider the simplest homogeneous, isotropic, constant-coefficient classical models and solutions with sufficient regularity: C^2 in space and C^1 in time for the heat equation, C^2 in both space and time for the wave equation, and C^2 in space for the Laplace and Poisson equations.
For thermal diffusivity \kappa>0, the heat equation is
u_t=\kappa\Delta u.
It describes the evolution of temperature toward its surrounding average. The Laplace equation
\Delta u=0
describes an equilibrium without an interior source, and its solutions are harmonic functions.
| Equation | Role of the Laplacian |
| Heat equation u_t=\kappa\Delta u | Diffuses temperature by reducing deviations from surrounding values |
| Wave equation u_{tt}=c^2\Delta u | Converts spatial tension or curvature into acceleration |
| Poisson equation -\Delta u=f | Determines a potential from the source f under the sign convention used here |
| Laplace equation \Delta u=0 | Describes a source-free equilibrium |
6Examples of Harmonic Functions
For u(x,y)=x^2-y^2, one has u_{xx}=2 and u_{yy}=-2, hence \Delta u=0. Thus u is harmonic and satisfies the mean-value property on every ball contained in its domain.
Another example is
u(x,y)=\log\sqrt{x^2+y^2},
which is harmonic away from the origin. Because it has a singularity at the origin, it must not be described as harmonic on the entire plane. This example leads to Green functions and fundamental solutions.
7Dependence of the Formula on Coordinates
In two-dimensional Cartesian coordinates, \Delta f=f_{xx}+f_{yy}. In polar coordinates with r>0,
\Delta f=f_{rr}+\frac{1}{r}f_r+\frac{1}{r^2}f_{\theta\theta}.
This formula follows by applying the chain rule to x=r\cos\theta and y=r\sin\theta and simplifying f_{xx}+f_{yy}. Polar coordinates are singular at r=0, so this expression cannot be substituted there directly. The Euclidean operator itself is unchanged, but its coordinate expression changes; terms such as 1/r arise from the coordinate transformation.
8Dimensional Check
In the heat equation, suppose that u is measured in \mathrm K and length in \mathrm m. Then \Delta u has units \mathrm{K/m^2}. If \kappa has units \mathrm{m^2/s}, then \kappa\Delta u has units \mathrm{K/s}, which agrees with u_t.
9Common Confusions
The scalar Laplacian acts on scalar fields. For a C^2 vector field \boldsymbol F in Cartesian coordinates on \mathbb R^3, the vector Laplacian may be defined componentwise, and
\Delta\boldsymbol F=\nabla(\nabla\cdot\boldsymbol F)-\nabla\times(\nabla\times\boldsymbol F).
Expanding the right side componentwise cancels the mixed partial derivatives and leaves \Delta F_i in each component. In curvilinear coordinates, merely differentiating the displayed components twice omits the variation of the basis and therefore does not produce the same operator.
10Scope of Validity
The discussion above concerns the Euclidean Laplacian and homogeneous, isotropic, constant-coefficient models. In an inhomogeneous or anisotropic medium, a variable-coefficient operator such as \nabla\cdot(A(x)\nabla u) may appear. On a general Riemannian manifold, the metric determines the Laplace--Beltrami operator. Moreover, the mean-value property follows from harmonicity on a region containing the ball, not merely from the equation \Delta f=0 at its center.
11Related Lectures
data/lecture/math/partial-differential-equations/heat-wave-and-laplace-equations.lecture.n.md
data/lecture/math/vector-calculus/green-gauss-and-stokes-theorems.lecture.n.md
data/lecture/math/vector-calculus/vector-calculus-portal.lecture.n.md
data/lecture/physics/waves/wave-equation-basics.lecture.n.md