§E15.33熱核と波動の明示解

最終更新

熱方程式と波動方程式は同じ Laplace 作用素を用いるが、局所的な初期状態が後の時刻に影響する範囲は異なる。例えば、零でない非負の滑らかな初期値がコンパクトな台をもつ場合、熱解は、どの正の時刻にも全空間で正となる。このような違いを調べるには、解の存在や一意性に加えて、初期値が解に寄与する仕方を具体的に表す必要がある。解を核との積分として表す方法は、偏微分方程式の基本的な解析手法の一つであり、核の台・符号・拡大縮小から、熱の平滑化と波動の伝播を捉えることができる。本記事では、熱核と波動の明示解を通して、両者の性質と空間次元による伝播の違いを解説する。

1 Gauss 核と近似単位

命題 1.1.n≥1n\geq1を整数とし、t>0t>0に対してRn\R^n上の Gauss 核を

Gt(x)=(4πt)−n/2exp⁡(−∣x∣24t)G_t(x)=(4\pi t)^{-n/2}\exp\left(-\frac{|x|^2}{4t}\right)

と定める。このときGt>0G_t>0、∫RnGt(x) dx=1\int_{\R^n}G_t(x)\,dx=1であり、任意のs,t>0s,t>0について

Gt∗Gs=Gt+sG_t*G_s=G_{t+s}

が成り立つ。1≤p<∞1\leq p<\infty、f∈Lp(Rn;C)f\in L^p(\R^n;\C)ならば

lim⁡t↓0∥Gt∗f−f∥p=0.\lim_{t\downarrow0}\|G_t*f-f\|_p=0.

また、f∈L1(Rn;C)f\in L^1(\R^n;\C)ならば

∫Rn(Gt∗f)(x) dx=∫Rnf(x) dx.\int_{\R^n}(G_t*f)(x)\,dx=\int_{\R^n}f(x)\,dx.

証明. 正値性は式から従う。位相因子e−2πix⋅ξe^{-2\pi ix\cdot\xi}による Fourier 変換について、§E15.26 定理 4.1により

G^t(ξ)=e−4π2t∣ξ∣2,Gt∗h=etΔh(h∈L2(Rn))\widehat G_t(\xi)=e^{-4\pi^2t|\xi|^2},\qquad G_t*h=e^{t\Delta}h\quad(h\in L^2(\R^n))

が成り立つ。ξ=0\xi=0とすれば∫Gt=1\int G_t=1である。Gs∈L2G_s\in L^2に対するetΔe^{t\Delta}の Fourier 乗数はe−4π2t∣ξ∣2e^{-4\pi^2t|\xi|^2}なので、Gt∗GsG_t*G_sの Fourier 変換はe−4π2(t+s)∣ξ∣2=G^t+se^{-4\pi^2(t+s)|\xi|^2}=\widehat G_{t+s}である。したがってGt∗Gs=Gt+sG_t*G_s=G_{t+s}をL2L^2で得る。両辺は連続関数であるから、各点でも等しい。

1≤p<∞1\leq p<\inftyとする。§E14.7 定理 2.2の Bochner 積分表示とノルム評価により

∥Gt∗f−f∥p≤∫RnGt(y)∥f(⋅−y)−f∥p dy.\|G_t*f-f\|_p \leq\int_{\R^n}G_t(y)\|f(\cdot-y)-f\|_p\,dy.

§E14.7 定理 1.2により、任意のε>0\varepsilon>0に対して、∣y∣<δ|y|<\deltaならば∥f(⋅−y)−f∥p<ε\|f(\cdot-y)-f\|_p<\varepsilonとなるδ>0\delta>0が存在する。したがって

∥Gt∗f−f∥p≤ε+2∥f∥p∫∣y∣≥δGt(y) dy.\|G_t*f-f\|_p \leq\varepsilon+2\|f\|_p\int_{|y|\geq\delta}G_t(y)\,dy.

y=t zy=\sqrt t\,zと変数変換すると、最後の積分は∫∣z∣≥δ/tG1(z) dz\int_{|z|\geq\delta/\sqrt t}G_1(z)\,dzであり、t↓0t\downarrow0で零へ収束する。ε\varepsilonは任意なので、LpL^p収束を得る。

f∈L1f\in L^1のとき、∬Gt(y)∣f(x−y)∣ dx dy=∥f∥1\iint G_t(y)|f(x-y)|\,dx\,dy=\|f\|_1である。Fubini の定理により積分順序を交換すると、質量保存の等式を得る。▨

定理 1.2.n≥1n\geq1、1≤p≤q≤∞1\leq p\leq q\leq\infty、f∈Lp(Rn;C)f\in L^p(\R^n;\C)とし、α\alphaを多重指数とする。このとき、u(t,x)=(Gt∗f)(x)u(t,x)=(G_t*f)(x)は(0,∞)×Rn(0,\infty)\times\R^n上で滑らかであり、ut=Δuu_t=\Delta uを満たす。また、ttとffに依存しない定数Cn,p,q,αC_{n,p,q,\alpha}に対して

∥Dxαu(t)∥q≤Cn,p,q,αt−∣α∣/2−n2(1/p−1/q)∥f∥p(t>0)\|D_x^\alpha u(t)\|_q \leq C_{n,p,q,\alpha} t^{-|\alpha|/2-\frac n2(1/p-1/q)}\|f\|_p \qquad(t>0)

が成り立つ。ただし1/∞=01/\infty=0とする。

証明. 核の拡大縮小から

DαGt(x)=t−(n+∣α∣)/2(DαG1)(x/t)D^\alpha G_t(x)=t^{-(n+|\alpha|)/2}(D^\alpha G_1)(x/\sqrt t)

である。DαG1D^\alpha G_1は多項式と Gauss 関数との積なので、任意の1≤r≤∞1\leq r\leq\inftyに対して

∥DαGt∥r=t−∣α∣/2−n2(1−1/r)∥DαG1∥r.\|D^\alpha G_t\|_r =t^{-|\alpha|/2-\frac n2(1-1/r)}\|D^\alpha G_1\|_r.

p<∞p<\inftyのとき、§E14.7 定理 2.2により

∥(DαGt)∗f∥p≤∥DαGt∥1∥f∥p.\|(D^\alpha G_t)*f\|_p\leq\|D^\alpha G_t\|_1\|f\|_p.

1/p+1/p′=11/p+1/p'=1とおく。§E10.22 定理 4.2を各xxで適用すると、端点を含む任意のppについて

∣((DαGt)∗f)(x)∣≤∥DαGt∥p′∥f∥p.|((D^\alpha G_t)*f)(x)|\leq\|D^\alpha G_t\|_{p'}\|f\|_p.

p<q<∞p<q<\inftyのとき、h=(DαGt)∗fh=(D^\alpha G_t)*fに対して

∥h∥qq=∫∣h∣q−p∣h∣p≤∥h∥∞q−p∥h∥pp\|h\|_q^q=\int |h|^{q-p}|h|^p \leq\|h\|_\infty^{q-p}\|h\|_p^p

である。二つの端の評価を代入すると、主張のttの指数を得る。p=q=∞p=q=\inftyのときは

∥(DαGt)∗f∥∞≤∥DαGt∥1∥f∥∞\|(D^\alpha G_t)*f\|_\infty \leq\|D^\alpha G_t\|_1\|f\|_\infty

を積分から直接得る。

核の任意の時間・空間微分は、正時間では多項式と Gauss 関数との積である。ttが正のコンパクト区間を動き、xxがコンパクト集合を動くとき、y↦∂tkDxαGt(x−y)y\mapsto\partial_t^kD_x^\alpha G_t(x-y)とその差商は、Lp′L^{p'}に属する共通の減衰関数で支配される。p′=∞p'=\inftyの場合にも同じ減衰関数による一様評価が成り立つ。Hölder の不等式と優収束により、畳込みの下で任意階の微分を行うことができ、各微分は連続である。直接微分による∂tGt=ΔGt\partial_tG_t=\Delta G_tからut=Δuu_t=\Delta uを得る。▨

命題 1.3.0≤f∈L1(Rn;R)0\leq f\in L^1(\R^n;\R)が零関数でないならば、任意のt>0t>0、x∈Rnx\in\R^nについて(Gt∗f)(x)>0(G_t*f)(x)>0である。特に、零でない非負のコンパクトな台をもつ初期値から出発しても、正時間の熱解の台はRn\R^n全体である。

証明.f>0f>0となる集合は正の測度をもち、Gt(x−y)>0G_t(x-y)>0はすべてのyyで成り立つ。したがって非負関数Gt(x−y)f(y)G_t(x-y)f(y)の積分は正である。積分の有限性と連続性は定理 1.2による。▨

注意 1.4.1{t>0}Gt(x)\mathbf1_{\{t>0\}}G_t(x)は§E15.6 定理 3.1の熱作用素の基本解である。実数値の有界連続初期値については、Gt∗fG_t*fは唯一の有界古典解であり、初期値へ局所一様に収束する§E15.26 命題 4.2。初期値が一様連続ならば、その収束は全空間で一様である。

命題 1.5.n≥1n\geq1とし、KtK_tを§E15.32 定理 3.1のトーラスTn=(R/Z)n\mathbb T^n=(\R/\Z)^nの熱核、GtG_tを命題 1.1の全空間の Gauss 核とする。Q=[−12,12]nQ=[-\tfrac12,\tfrac12]^nの点をTn\mathbb T^nの点の代表と見なす。このとき、次の二つが成り立つ。

  1. nnだけに依存する定数CnC_nが存在し、0<t≤10<t\leq1とx∈Qx\in Qに対して 0<Kt(x)−Gt(x)≤Cn(4πt)−n/2e−1/(32t)0<K_t(x)-G_t(x)\leq C_n(4\pi t)^{-n/2}e^{-1/(32t)} が成り立つ。特に、sup⁡x∈Q∣Kt(x)−Gt(x)∣\sup_{x\in Q}|K_t(x)-G_t(x)|はt↓0t\downarrow0で零へ収束する。
  2. nnだけに依存する定数Cn′C_n'が存在し、t≥1t\geq1とx∈Qx\in Qに対して ∣Kt(x)−1∣≤Cn′e−4π2t,0<Gt(x)≤(4πt)−n/2|K_t(x)-1|\leq C_n'e^{-4\pi^2t},\qquad0<G_t(x)\leq(4\pi t)^{-n/2} が成り立つ。特に、t→∞t\to\inftyでKtK_tは定数11へ一様収束し、GtG_tは零へ一様収束する。また、u0∈L1(Tn)u_0\in L^1(\mathbb T^n)に対して(Kt∗u0)(x)=∫TnKt(x−y)u0(y) dy(K_t*u_0)(x)=\int_{\mathbb T^n}K_t(x-y)u_0(y)\,dyとおくと、∣(Kt∗u0)(x)−∫Tnu0 dy∣≤Cn′e−4π2t∥u0∥1|(K_t*u_0)(x)-\int_{\mathbb T^n}u_0\,dy|\leq C_n'e^{-4\pi^2t}\|u_0\|_1が成り立つ。

証明.主張 (1)を示す。§E15.32 定理 3.1によりx∈Qx\in Qに対してKt(x)=∑m∈ZnGt(x+m)K_t(x)=\sum_{m\in\Z^n}G_t(x+m)であるから、Kt(x)−Gt(x)=∑m≠0Gt(x+m)>0K_t(x)-G_t(x)=\sum_{m\ne0}G_t(x+m)>0である。m≠0m\ne0の最大ノルムをj=∣m∣∞≥1j=|m|_\infty\geq1とすると、∣x+m∣≥∣x+m∣∞≥j−12≥j/2|x+m|\geq|x+m|_\infty\geq j-\tfrac12\geq j/2であるからGt(x+m)≤(4πt)−n/2e−j2/(16t)G_t(x+m)\leq(4\pi t)^{-n/2}e^{-j^2/(16t)}である。0<t≤10<t\leq1ではj2/(16t)≥1/(32t)+j2/32j^2/(16t)\geq1/(32t)+j^2/32が成り立つ。∣m∣∞=j|m|_\infty=jを満たすmmの個数は(2j+1)n−(2j−1)n≤2n(2j+1)n−1(2j+1)^n-(2j-1)^n\leq2n(2j+1)^{n-1}であるから

Kt(x)−Gt(x)≤(4πt)−n/2e−1/(32t)∑j≥12n(2j+1)n−1e−j2/32K_t(x)-G_t(x)\leq(4\pi t)^{-n/2}e^{-1/(32t)}\sum_{j\geq1}2n(2j+1)^{n-1}e^{-j^2/32}

である。右辺の級数は収束するので、その値をCnC_nとすればよい。t−n/2e−1/(32t)→0t^{-n/2}e^{-1/(32t)}\to0により後半を得る。

主張 (2)を示す。§E15.32 定理 3.1によりKt(x)=1+∑k≠0e−4π2t∣k∣2e2πik⋅xK_t(x)=1+\sum_{k\ne0}e^{-4\pi^2t|k|^2}e^{2\pi ik\cdot x}である。k∈Zn∖{0}k\in\Z^n\setminus\{0\}とt≥1t\geq1に対してt∣k∣2≥(t−1)+∣k∣2t|k|^2\geq(t-1)+|k|^2であるから

∣Kt(x)−1∣≤e−4π2(t−1)∑k≠0e−4π2∣k∣2|K_t(x)-1|\leq e^{-4\pi^2(t-1)}\sum_{k\ne0}e^{-4\pi^2|k|^2}

となり、Cn′=e4π2∑k≠0e−4π2∣k∣2C_n'=e^{4\pi^2}\sum_{k\ne0}e^{-4\pi^2|k|^2}とおけばよい。Gt(x)≤Gt(0)=(4πt)−n/2G_t(x)\leq G_t(0)=(4\pi t)^{-n/2}は式から従う。最後の評価は、∫TnKt(x−y) dy=1\int_{\mathbb T^n}K_t(x-y)\,dy=1により

(Kt∗u0)(x)−∫Tnu0 dy=∫Tn(Kt(x−y)−1)u0(y) dy(K_t*u_0)(x)-\int_{\mathbb T^n}u_0\,dy=\int_{\mathbb T^n}\bigl(K_t(x-y)-1\bigr)u_0(y)\,dy

と書き、被積分関数の核の部分を上の評価で抑えれば従う。▨

例 1.6.s>0s>0を固定し、初期値をf=Gsf=G_sとする。このとき、命題 1.1により、熱解はu(t)=Gt+su(t)=G_{t+s}である。a=t+sa=t+sとおくと、Gauss 積分から

∥u(t)∥22=(8πa)−n/2.\|u(t)\|_2^2=(8\pi a)^{-n/2}.

∇Ga=−xGa/(2a)\nabla G_a=-xG_a/(2a)であり、各座標について部分積分すると

∫Rn∣x∣2e−∣x∣2/(2a) dx=na∫Rne−∣x∣2/(2a) dx.\int_{\R^n}|x|^2e^{-|x|^2/(2a)}\,dx =na\int_{\R^n}e^{-|x|^2/(2a)}\,dx.

したがって

∥∇u(t)∥22=n4a(8πa)−n/2,ddt∥u(t)∥22=−n2a(8πa)−n/2=−2∥∇u(t)∥22.\|\nabla u(t)\|_2^2=\frac n{4a}(8\pi a)^{-n/2},\qquad \frac{d}{dt}\|u(t)\|_2^2 =-\frac n{2a}(8\pi a)^{-n/2}=-2\|\nabla u(t)\|_2^2.

これはA=−ΔA=-\Deltaに対する熱のエネルギー等式§E15.26 定理 1.1と一致する。

例 1.7. 一次元でf=1(0,∞)f=\mathbf1_{(0,\infty)}とする。このとき、Φ(a)=∫−∞aG1(z) dz\Phi(a)=\int_{-\infty}^aG_1(z)\,dzとおくと

(Gt∗f)(x)=Φ(x/t).(G_t*f)(x)=\Phi(x/\sqrt t).

G1G_1は正で偶関数なので、Φ\Phiは狭義単調増加であり、Φ(0)=1/2\Phi(0)=1/2である。x<0x<0では誤差はΦ(x/t)<1/2\Phi(x/\sqrt t)<1/2、x>0x>0では誤差は1−Φ(x/t)<1/21-\Phi(x/\sqrt t)<1/2である。両側からx→0x\to0とすると誤差は1/21/2に収束するので

∥Gt∗f−f∥∞=12(t>0).\|G_t*f-f\|_\infty=\frac12\qquad(t>0).

したがって、近似単位のLpL^pノルム初期収束から有限指数という仮定を除くことはできない。

2 一次元の波動解

定理 2.1.f∈C2(R;C)f\in C^2(\R;\C)、g∈C1(R;C)g\in C^1(\R;\C)とする。このとき

utt−uxx=0,u(0,x)=f(x),ut(0,x)=g(x)u_{tt}-u_{xx}=0,\qquad u(0,x)=f(x),\qquad u_t(0,x)=g(x)

を満たす[0,∞)×R[0,\infty)\times\R上の古典解は唯一であり、

u(t,x)=f(x+t)+f(x−t)2+12∫x−tx+tg(y) dyu(t,x)=\frac{f(x+t)+f(x-t)}2 +\frac12\int_{x-t}^{x+t}g(y)\,dy

で与えられる。S1(t)h(x)=12∫x−tx+th(y) dyS_1(t)h(x)=\frac12\int_{x-t}^{x+t}h(y)\,dyとおけば、同じ公式は

u(t)=∂tS1(t)f+S1(t)gu(t)=\partial_tS_1(t)f+S_1(t)g

となる。

証明.r=x+tr=x+t、s=x−ts=x-tと変数変換すると、utt−uxx=−4ursu_{tt}-u_{xx}=-4u_{rs}である。urs=0u_{rs}=0を積分すると、各古典解はu(t,x)=a(x+t)+b(x−t)u(t,x)=a(x+t)+b(x-t)の形をもつ。初期条件からa+b=fa+b=f、a′−b′=ga'-b'=gである。よってa′=(f′+g)/2a'=(f'+g)/2、b′=(f′−g)/2b'=(f'-g)/2を積分して主張の式を得る。逆に、その式を微分すると方程式と両初期条件が成り立つ。

t>0t>0のときS1(t)S_1(t)の畳込み核は121[−t,t]\frac12\mathbf1_{[-t,t]}である。§E15.6 命題 4.2により、この核は因果的波動基本解の時刻ttにおける空間核と一致する。積分端を微分すると∂tS1(t)f=(f(⋅+t)+f(⋅−t))/2\partial_tS_1(t)f=(f(\cdot+t)+f(\cdot-t))/2を得る。▨

命題 2.2.f∈H1(R;C)f\in H^1(\R;\C)、g∈L2(R;C)g\in L^2(\R;\C)、T>0T>0とする。このとき、定理 2.1の式によって定めたuuは、

C([0,T];H1(R))∩C1([0,T];L2(R))C([0,T];H^1(\R))\cap C^1([0,T];L^2(\R))

に属する唯一のエネルギー解である。任意のT>0T>0とF∈C1([0,T]×R;C)F\in C^1([0,T]\times\R;\C)に対して、零初期値の古典解の外力項は

w(t,x)=12∫0t∫x−(t−s)x+(t−s)F(s,y) dy dsw(t,x)=\frac12\int_0^t\int_{x-(t-s)}^{x+(t-s)}F(s,y)\,dy\,ds

であり、wtt−wxx=Fw_{tt}-w_{xx}=Fを満たす。さらにF∈L1(0,T;L2(R))F\in L^1(0,T;L^2(\R))の場合も、この式の積分をL2L^2値の積分として定めると、得られる関数wwは零初期値のエネルギー解である。

証明.τah(x)=h(x+a)\tau_a h(x)=h(x+a)とおく。L2L^2における平行移動の等長性と連続性§E14.7 定理 1.2、および積分表示から

∥S1(t)h∥2≤t∥h∥2,∂xS1(t)h=τth−τ−th2,∂tS1(t)h=τth+τ−th2\|S_1(t)h\|_2\leq t\|h\|_2,\qquad \partial_xS_1(t)h=\frac{\tau_th-\tau_{-t}h}{2},\qquad \partial_tS_1(t)h=\frac{\tau_th+\tau_{-t}h}{2}

が成り立つ。後二式の空間微分は弱微分、時間微分はL2L^2ノルム微分である。空間微分の式は試験関数に対する積分端の計算から従い、時間微分の式は平行移動の連続性と積分の基本定理から従う。

H1H^1とL2L^2でそれぞれf,gf,gを滑らかなコンパクトな台をもつ関数で近似する。S1(t)S_1(t)とその時間・空間微分のL2L^2評価を近似解の差に適用すると、近似解はC([0,T];H1)∩C1([0,T];L2)C([0,T];H^1)\cap C^1([0,T];L^2)で主張の式による関数へ収束する。近似解の方程式を試験関数と組み合わせて極限を取ると、極限も弱方程式と初期条件を満たす。§E15.27 定理 2.2をA=−∂x2A=-\partial_x^2に適用すれば一意性を得る。

F∈C1F\in C^1の場合、積分端と時間を微分すると

wt(t,x)=12∫0t(F(s,x+t−s)+F(s,x−t+s)) ds,wtt(t,x)=F(t,x)+12∫0t(Fx(s,x+t−s)−Fx(s,x−t+s)) ds,wxx(t,x)=12∫0t(Fx(s,x+t−s)−Fx(s,x−t+s)) ds.\begin{aligned} w_t(t,x)&=\frac12\int_0^t\bigl(F(s,x+t-s)+F(s,x-t+s)\bigr)\,ds,\\ w_{tt}(t,x)&=F(t,x)+\frac12\int_0^t\bigl(F_x(s,x+t-s)-F_x(s,x-t+s)\bigr)\,ds,\\ w_{xx}(t,x)&=\frac12\int_0^t\bigl(F_x(s,x+t-s)-F_x(s,x-t+s)\bigr)\,ds. \end{aligned}

したがって方程式と零初期条件が成り立つ。F∈Lt1Lx2F\in L_t^1L_x^2のときはw(t)=∫0tS1(t−s)F(s) dsw(t)=\int_0^tS_1(t-s)F(s)\,dsがH1H^1値の Bochner 積分となる。t>0t>0のときS1(t)S_1(t)は核121[−t,t]\frac12\mathbf1_{[-t,t]}との畳込みであり、この核は§E15.6 命題 4.2により因果的波動基本解の時刻ttにおける空間核である。その Fourier 変換は

12∫−tte−2πixξ dx=sin⁡(2πtξ)2πξ\frac12\int_{-t}^te^{-2\pi ix\xi}\,dx=\frac{\sin(2\pi t\xi)}{2\pi\xi}

であり、これはA=−∂x2A=-\partial_x^2の Fourier 乗数λ=4π2ξ2\lambda=4\pi^2\xi^2に対するSA(t)S_A(t)の乗数sin⁡(tλ)/λ\sin(t\sqrt\lambda)/\sqrt\lambdaに等しい。有界な乗数が等しいからS1(t)=SA(t)S_1(t)=S_A(t)である。したがってw(t)=∫0tSA(t−s)F(s) dsw(t)=\int_0^tS_A(t-s)F(s)\,dsは§E15.27 定理 2.2の Duhamel 公式の外力項であり、この積分はエネルギー解である。▨

例 2.3.f∈Cc∞(R;R)f\in C_c^\infty(\R;\R)、g=−f′g=-f'とする。このとき、積分の基本定理と定理 2.1により

u(t,x)=f(x+t)+f(x−t)2−f(x+t)−f(x−t)2=f(x−t).u(t,x)=\frac{f(x+t)+f(x-t)}2 -\frac{f(x+t)-f(x-t)}2=f(x-t).

したがって

12∫R(∣ut(t,x)∣2+∣ux(t,x)∣2) dx=∫R∣f′(x−t)∣2 dx=∥f′∥22.\frac12\int_\R\bigl(|u_t(t,x)|^2+|u_x(t,x)|^2\bigr)\,dx =\int_\R|f'(x-t)|^2\,dx=\|f'\|_2^2.

3 球面平均と Kirchhoff の公式

補題 3.1.n≥1n\geq1を整数、h∈C2(Rn;C)h\in C^2(\R^n;\C)とする。σn−1=∣Sn−1∣\sigma_{n-1}=|S^{n-1}|とし、球面平均を

Mrh(x)=1σn−1∫Sn−1h(x+rω) dσ(ω)(r≥0)M_rh(x)=\frac1{\sigma_{n-1}}\int_{S^{n-1}}h(x+r\omega)\,d\sigma(\omega) \qquad(r\geq0)

と定める。n=1n=1ではS0={−1,1}S^0=\{-1,1\}上の数え上げ測度を用いる。このときr>0r>0において

(∂r2+n−1r∂r)Mrh=MrΔh=ΔxMrh.\left(\partial_r^2+\frac{n-1}{r}\partial_r\right)M_rh =M_r\Delta h=\Delta_xM_rh.

またM0h=hM_0h=h、∂rMrh∣r=0=0\partial_rM_rh|_{r=0}=0である。

証明.n≥2n\geq2のとき、§E15.6 補題 2.1をU=RnU=\R^n、f=hf=h、a=xa=xとして適用すると、0<r0<rに対して

σn−1rn−1∂rMrh(x)=∫B(x,r)Δh(y) dy\sigma_{n-1}r^{n-1}\partial_rM_rh(x)=\int_{B(x,r)}\Delta h(y)\,dy

が成り立つ。右辺を極座標で表すとσn−1∫0rsn−1MsΔh(x) ds\sigma_{n-1}\int_0^rs^{n-1}M_s\Delta h(x)\,dsであり、被積分関数はssについて連続であるから、rrで微分して

∂r(rn−1∂rMrh(x))=rn−1MrΔh(x)\partial_r\bigl(r^{n-1}\partial_rM_rh(x)\bigr) =r^{n-1}M_r\Delta h(x)

を得る。rn−1r^{n-1}で割れば最初の等式を得る。n=1n=1ではMrh(x)=(h(x+r)+h(x−r))/2M_rh(x)=(h(x+r)+h(x-r))/2の二回微分から同じ式を得る。xxに関する微分と球面積分の交換からΔxMrh=MrΔh\Delta_xM_rh=M_r\Delta hである。r=0r=0での微分は∇h(x)⋅σn−1−1∫Sn−1ω dσ=0\nabla h(x)\cdot\sigma_{n-1}^{-1}\int_{S^{n-1}}\omega\,d\sigma=0である。▨

補題 3.2.n≥1n\geq1とし、v∈C2([0,T]×Rn;C)v\in C^2([0,T]\times\R^n;\C)がvtt−Δv=0v_{tt}-\Delta v=0を満たすとする。このとき、あるx0∈Rnx_0\in\R^nとR>0R>0に対し、v(0)=vt(0)=0v(0)=v_t(0)=0がB(x0,R)B(x_0,R)上で成り立つならば

v(t,x)=0(0≤t≤T, ∣x−x0∣<R−t)v(t,x)=0\qquad(0\leq t\leq T,\ |x-x_0|<R-t)

である。特に、初期値が同じ二つの古典解は一致する。

証明.e=(∣vt∣2+∣∇v∣2)/2e=(|v_t|^2+|\nabla v|^2)/2、J=ℜ(vt‾∇v)J=\Re(\overline{v_t}\nabla v)とおくと、方程式から∂te=div⁡J\partial_te=\operatorname{div}J、∣J∣≤e|J|\leq eである。0≤t<min⁡{T,R}0\leq t<\min\{T,R\}に対し、球面での積分と Gauss の公式から

ddt∫B(x0,R−t)e(t,x) dx=∫∂B(x0,R−t)(J⋅ν−e) dσ≤0.\frac{d}{dt}\int_{B(x_0,R-t)}e(t,x)\,dx =\int_{\partial B(x_0,R-t)}(J\cdot\nu-e)\,d\sigma\leq0.

初期エネルギーは零なので、vt=0v_t=0は縮小する球の中で成り立つ。∣x−x0∣<R−t|x-x_0|<R-tならば、すべての0≤s≤t0\leq s\leq tでx∈B(x0,R−s)x\in B(x_0,R-s)である。よってv(t,x)=v(0,x)+∫0tvs(s,x) ds=0v(t,x)=v(0,x)+\int_0^tv_s(s,x)\,ds=0となる。端の時刻は連続性で扱う。同じ初期値をもつ解の差に適用し、各点についてR>tR>tを選べば一意性を得る。▨

定理 3.3.n=3n=3、f∈C3(R3;C)f\in C^3(\R^3;\C)、g∈C2(R3;C)g\in C^2(\R^3;\C)とする。このとき

u(t,x)=∂t[tMtf(x)]+tMtg(x)u(t,x)=\partial_t\bigl[tM_tf(x)\bigr]+tM_tg(x)

はutt−Δu=0u_{tt}-\Delta u=0、u(0)=fu(0)=f、ut(0)=gu_t(0)=gを満たす唯一の古典解である。t>0t>0に対して

S3(t)h(x)=tMth(x)=14πt∫∣y−x∣=th(y) dσ(y)S_3(t)h(x)=tM_th(x) =\frac1{4\pi t}\int_{|y-x|=t}h(y)\,d\sigma(y)

とおけば、u=∂tS3(t)f+S3(t)gu=\partial_tS_3(t)f+S_3(t)gである。

証明.h∈C2(R3)h\in C^2(\R^3)に対しv(t,x)=tMth(x)v(t,x)=tM_th(x)とおく。補題 3.1により、t>0t>0では

vtt=2∂tMth+t∂t2Mth=tΔxMth=Δxv.v_{tt}=2\partial_tM_th+t\partial_t^2M_th =t\Delta_xM_th=\Delta_xv.

同じ補題によりv(0)=0v(0)=0、vt(0)=hv_t(0)=h、vtt(0)=0v_{tt}(0)=0である。h=f∈C3h=f\in C^3の場合、vvの方程式を時間微分して∂tv\partial_tvも波動方程式を満たす。h=g∈C2h=g\in C^2の場合のvvと加えると、主張のuuはC2C^2級であり、方程式と初期条件を満たす。一意性は補題 3.2による。

t>0t>0のとき、tMttM_tの畳込み核は球面∣y∣=t|y|=t上の測度(4πt)−1dσ(y)(4\pi t)^{-1}d\sigma(y)である。§E15.6 命題 4.2により、この球面測度は因果的波動基本解の時刻ttにおける空間核と一致する。▨

注意 3.4. Kirchhoff の公式の各積分はコンパクトな球面上の積分なので、初期値が全空間でL2L^2に属する必要はない。T>0T>0と有界な空間領域Ω⊂R3\Omega\subset\R^3を固定する。[0,T]×Ω[0,T]\times\Omega上の解を求める際には、Ω\Omegaから距離TT以下の点を含む集合上で初期値を保つ滑らかな切断を取ることができる。切断後のデータはH1×L2H^1\times L^2に属し、公式はそのエネルギー解と一致する§E15.27 定理 2.2。切断の選択によらず局所解が一致することは、古典解については補題 3.2、エネルギー解については§E15.27 定理 5.1による。

4 Hadamard の降下法と伝播

定理 4.1.f∈C3(R2;C)f\in C^3(\R^2;\C)、g∈C2(R2;C)g\in C^2(\R^2;\C)とする。t>0t>0とh∈C(R2;C)h\in C(\R^2;\C)に対し

S2(t)h(x)=12π∫∣x−y∣<th(y)t2−∣x−y∣2 dyS_2(t)h(x)=\frac1{2\pi}\int_{|x-y|<t}\frac{h(y)}{\sqrt{t^2-|x-y|^2}}\,dy

と定める。このとき

u(t,x)=∂tS2(t)f(x)+S2(t)g(x)u(t,x)=\partial_tS_2(t)f(x)+S_2(t)g(x)

はutt−Δu=0u_{tt}-\Delta u=0、u(0)=fu(0)=f、ut(0)=gu_t(0)=gを満たす唯一の二次元の古典解である。t=0t=0での値は連続延長によって定める。

証明.hhを第三座標に依存しない関数H(x1,x2,x3)=h(x1,x2)H(x_1,x_2,x_3)=h(x_1,x_2)に延長する。半径ttの球面の上下半球をz=(z1,z2)z=(z_1,z_2)の円板へ射影すると、それぞれのグラフz3=±t2−∣z∣2z_3=\pm\sqrt{t^2-|z|^2}の面積要素は

dσ=tt2−∣z∣2 dzd\sigma=\frac{t}{\sqrt{t^2-|z|^2}}\,dz

である。したがって、二つの半球を加えることにより

S3(t)H(x,x3)=12π∫∣z∣<th(x+z)t2−∣z∣2 dz=S2(t)h(x).S_3(t)H(x,x_3) =\frac1{2\pi}\int_{|z|<t}\frac{h(x+z)}{\sqrt{t^2-|z|^2}}\,dz=S_2(t)h(x).

特異因子は可積分である。実際、z=twz=twと変数変換すれば

S2(t)h(x)=t2π∫∣w∣<1h(x+tw)1−∣w∣2 dw,∫∣w∣<1dw1−∣w∣2=2π.S_2(t)h(x)=\frac t{2\pi}\int_{|w|<1}\frac{h(x+tw)}{\sqrt{1-|w|^2}}\,dw, \qquad \int_{|w|<1}\frac{dw}{\sqrt{1-|w|^2}}=2\pi.

固定した円板上の表示では、データの許す次数まで積分の下で微分することができる。

延長した初期値に定理 3.3を適用する。得られる解は第三座標に依存しないので、三次元の Laplace 作用素は二次元の Laplace 作用素に等しい。よって主張の式が二次元の方程式と初期条件を満たす。一意性は補題 3.2による。第三座標に一定な延長を全空間のL2L^2データと見なす必要はない。有限な依存領域で延長を切断した場合にも、§E15.27 定理 5.1により同じ局所解を得る。▨

命題 4.2. 三次元の斉次波動方程式の初期値をf∈C3(R3)f\in C^3(\R^3)、g∈C2(R3)g\in C^2(\R^3)とする。このとき、t>0t>0、x∈R3x\in\R^3に対し、ffとggが球面{y:∣y−x∣=t}\{y:|y-x|=t\}のある開近傍で零ならば、解はu(t,x)=0u(t,x)=0を満たす。

証明. 球面はコンパクトであるから、あるε>0\varepsilon>0が存在して、∣r−t∣<ε|r-t|<\varepsilonならば球面∣y−x∣=r|y-x|=r上でf=g=0f=g=0となる。したがってMrf(x)=Mrg(x)=0M_rf(x)=M_rg(x)=0がそのrrの範囲で成り立つ。定理 3.3の二項は、ともに零である。▨

例 4.3. 各次元n∈{1,2,3}n\in\{1,2,3\}でf=0f=0とし、0≤g∈Cc∞(Rn;R)0\leq g\in C_c^\infty(\R^n;\R)、g≠0g\ne0を取る。supp⁡g⊂B(0,R)\supp g\subset B(0,R)とし、xxを固定する。このとき、t>∣x∣+Rt>|x|+Rならば、三次元では球面∣y−x∣=t|y-x|=tがsupp⁡g\supp gから離れるので

u(t,x)=14πt∫∣y−x∣=tg(y) dσ(y)=0.u(t,x)=\frac1{4\pi t}\int_{|y-x|=t}g(y)\,d\sigma(y)=0.

二次元では、同じ時刻にsupp⁡g⊂B(x,t)\supp g\subset B(x,t)であり、

u(t,x)=12π∫supp⁡gg(y)t2−∣x−y∣2 dy>0.u(t,x)=\frac1{2\pi}\int_{\supp g}\frac{g(y)}{\sqrt{t^2-|x-y|^2}}\,dy>0.

しかもt→∞t\to\inftyにおいて、台上で一様にt/t2−∣x−y∣2→1t/\sqrt{t^2-|x-y|^2}\to1となるので

t u(t,x)⟶12π∫R2g(y) dy>0.t\,u(t,x)\longrightarrow\frac1{2\pi}\int_{\R^2}g(y)\,dy>0.

一次元では、区間[x−t,x+t][x-t,x+t]が台全体を含むため

u(t,x)=12∫Rg(y) dy>0.u(t,x)=\frac12\int_\R g(y)\,dy>0.

いずれの次元でも、解に影響する初期値は閉球B(x,t)‾\overline{B(x,t)}の中にある。三次元ではさらに球面の近傍だけで十分であるという強い Huygens の原理が成り立つが、上の一次元と二次元の初期値は球面近傍で零であっても解を零にしない。

前提記事