§E15.26熱方程式

最終更新

熱の拡散に伴う温度分布の変化は、熱方程式によって記述される。例えば、区間(0,π)(0,\pi)の両端を温度零に保つと、初期温度分布sin⁡x\sin xは時刻ttにe−tsin⁡xe^{-t}\sin xとなるが、一般の初期値については、解の構成とともに、時間の経過が分布の滑らかさや極限に及ぼす影響を明らかにする必要がある。非負自己共役作用素の熱半群では、正のスペクトル値に対応する成分は指数的に減衰し、零固有空間の成分は変化しない。このスペクトル表示によって、解の滑らかさと長時間挙動を作用素の言葉で扱うことができる。熱方程式は放物型方程式の最も基本的な例の一つであり、高いスペクトル成分の減衰から得られる平滑化と、最大値原理による点ごとの制御に、その特徴が現れる。本記事では、熱半群を通して熱方程式の解とその基本的な性質を解説する。

1 熱半群と平滑化

定理 1.1.Ω⊆Rn\Omega\subseteq\R^nを開集合、H=L2(Ω;C)H=L^2(\Omega;\C)、AAをHH上の非負自己共役作用素とする。内積は第一変数について線形とする。u0∈Hu_0\in Hに対し

u(t)=e−tAu0(t≥0)u(t)=e^{-tA}u_0\qquad(t\geq0)

とおく。このときu∈C([0,∞);H)u\in C([0,\infty);H)、u(0)=u0u(0)=u_0であり、すべての整数m≥0m\geq0とt>0t>0について

u(t)∈D(Am),∥Ame−tA∥≤Cmt−m,C0=1,Cm=(m/e)m (m≥1)u(t)\in D(A^m),\qquad \|A^me^{-tA}\|\leq C_mt^{-m},\qquad C_0=1,\quad C_m=(m/e)^m\ (m\geq1)

が成り立つ。またu∈C∞((0,∞);H)u\in C^\infty((0,\infty);H)であり、整数k≥0k\geq0に対して

∂tku(t)=(−A)ku(t),12ddt∥u(t)∥22=−∥A1/2u(t)∥22\partial_t^ku(t)=(-A)^ku(t),\qquad \frac12\frac{d}{dt}\|u(t)\|_2^2=-\|A^{1/2}u(t)\|_2^2

が成り立つ。uuはC([0,∞);H)∩C1((0,∞);H)C([0,\infty);H)\cap C^1((0,\infty);H)に属し、正時間にD(A)D(A)に値をとる初期値問題v′+Av=0v'+Av=0、v(0)=u0v(0)=u_0の唯一の解である。

証明.EEをAAのスペクトル測度、μ(B)=∥E(B)u0∥22\mu(B)=\|E(B)u_0\|_2^2とおく。μ\muは[0,∞)[0,\infty)上の有限測度である。∣e−tλ∣≤1|e^{-t\lambda}|\leq1と優収束定理により、t↦e−tAu0t\mapsto e^{-tA}u_0は[0,∞)[0,\infty)で連続であり、t=0t=0でu0u_0に等しい。冪の評価と定義域への所属は§E12.25 命題 2.4による。

0<a<b<∞0<a<b<\inftyと整数k≥0k\geq0をとる。t∈[a,b]t\in[a,b]に対しλ2ke−2tλ\lambda^{2k}e^{-2t\lambda}はλ≥0\lambda\geq0上で一様に有界である。差商にも、平均値の積分表示からCa,kλk+1e−aλ/2C_{a,k}\lambda^{k+1}e^{-a\lambda/2}という上界を与えることができる。スペクトル積分で優収束定理を用いると、(−A)ke−tAu0(-A)^ke^{-tA}u_0は連続であり、その微分は(−A)k+1e−tAu0(-A)^{k+1}e^{-tA}u_0である。作用素積と乗数の一致には§E12.15 定理 1.2 (4)を用いた。特にu′=−Auu'=-Auである。

内積を微分し、スペクトル積分の等式を用いると

12ddt∥u(t)∥22=ℜ(−Au(t),u(t))2=−∫[0,∞)λ d∥E(λ)u(t)∥22=−∥A1/2u(t)∥22\frac12\frac{d}{dt}\|u(t)\|_2^2 =\Re(-Au(t),u(t))_2 =-\int_{[0,\infty)}\lambda\,d\|E(\lambda)u(t)\|_2^2 =-\|A^{1/2}u(t)\|_2^2

を得る。二つの解の差wwに同じ計算を適用すると、0<s<t0<s<tに対して∥w(t)∥2≤∥w(s)∥2\|w(t)\|_2\leq\|w(s)\|_2である。s↓0s\downarrow0とすればw(t)=0w(t)=0となる。▨

命題 1.2.Ω⊂Rn\Omega\subset\R^nを空でない有界C∞C^\infty開集合、ADA_Dをその Dirichlet Laplace 作用素とする。u0∈L2(Ω)u_0\in L^2(\Omega)に対しu(t)=e−tADu0u(t)=e^{-tA_D}u_0とおく。このとき、任意の整数k,m≥0k,m\geq0について、∂tku\partial_t^kuは正時間でH2m(Ω)H^{2m}(\Omega)に値をとる連続写像である。uuは(0,∞)×Ω‾(0,\infty)\times\overline\Omegaで時間と空間について滑らかであり、すべての整数j≥0j\geq0とt>0t>0についてTΔju(t)=0T\Delta^ju(t)=0が成り立つ。

証明.定理 1.1と§E15.23 系 6.1により、u(t)∈H2m(Ω)u(t)\in H^{2m}(\Omega)であり、0≤j<m0\leq j<mについてTΔju(t)=0T\Delta^ju(t)=0である。ADm∂tku=(−1)kADm+kuA_D^m\partial_t^ku=(-1)^kA_D^{m+k}uのL2L^2連続性と同じノルム同値により、∂tku\partial_t^kuはH2mH^{2m}で連続である。2m>r+n/22m>r+n/2となるmmを選び、高階 Sobolev–Hölder 埋め込み§E15.10 系 4.2を適用すると、すべての時間・空間微分が境界まで連続となる。▨

注意 1.3. 非負自己共役な楕円型実現AAに対しても、作用素冪への所属は定理 1.1により成り立つ。空間正則性への移行には楕円型評価を要する。例えば、整数k≥0k\geq0とu0∈L2(Ω)u_0\in L^2(\Omega)をとり、u(t)=e−tAu0u(t)=e^{-tA}u_0とおく。AAがL=−div⁡(a∇)+b⋅∇+cL=-\operatorname{div}(a\nabla)+b\cdot\nabla+cの斉次 Dirichlet 実現であり、Ω\Omegaが有界Ck+1,1C^{k+1,1}開集合、a∈Wk+1,∞a\in W^{k+1,\infty}、b,c∈Wk,∞b,c\in W^{k,\infty}で、一様楕円性を満たすとする。正時間のuuとAuAuに§E15.23 定理 5.1 (1)を00階から順に適用することによりu(t)∈Hk+2u(t)\in H^{k+2}を得る。各段階ではAu(t)Au(t)に対しても、さらに高い作用素冪への所属を用いて一段前の正則性を適用する。係数と境界の正則性が有限の場合には、この議論から無限階の空間正則性は従わない。

注意 1.4.H=L2(Ω;C)H=L^2(\Omega;\C)上の非負自己共役作用素AA、τ>0\tau>0、u0∈Hu_0\in H、可積分な外力f∈L1(0,τ;H)f\in L^1(0,\tau;H)をとる。このとき、非斉次方程式u′+Au=fu'+Au=f、u(0)=u0u(0)=u_0の軟解u(t)=e−tAu0+∫0te−(t−s)Af(s) dsu(t)=e^{-tA}u_0+\int_0^te^{-(t-s)A}f(s)\,dsの連続性、一意性、データの差に対する評価は、§E15.24 定理 2.1を熱の場合(β=1\beta=1)に適用すると得られる。

例 1.5.Ω=(0,π)\Omega=(0,\pi)の Dirichlet Laplace 作用素と初期値u0=1u_0=1をとる。このとき、定理 1.1によりu(t)=e−tAD1u(t)=e^{-tA_D}1はt↓0t\downarrow0にL2(0,π)L^2(0,\pi)で11に収束する。一方、正時間では命題 1.2によりu(t,0)=u(t,π)=0u(t,0)=u(t,\pi)=0である。したがって

sup⁡x∈[0,π]∣u(t,x)−1∣≥1\sup_{x\in[0,\pi]}|u(t,x)-1|\geq1

であり、閉区間上の一様収束は成立しない。

2 固有関数展開と長時間挙動

定理 2.1.Ω⊂Rn\Omega\subset\R^nを空でない有界 Lipschitz 開集合、B∈{D,N}B\in\{D,N\}、ABA_Bを Dirichlet または Neumann Laplace 作用素とする。L2(Ω;C)L^2(\Omega;\C)の内積を第一変数について線形とし、(λjB,ejB)(\lambda_j^B,e_j^B)を§E15.19 定理 1.3の完全正規直交固有系とする。このとき、u0∈L2(Ω)u_0\in L^2(\Omega)に対して

e−tABu0=∑j=1∞e−tλjB(u0,ejB)2ejBe^{-tA_B}u_0=\sum_{j=1}^\infty e^{-t\lambda_j^B}(u_0,e_j^B)_2 e_j^B

が成り立つ。級数はt≥0t\geq0でL2L^2収束し、t>0t>0ではすべての整数m≥0m\geq0に対してD(ABm)D(A_B^m)のグラフノルムで収束する。

証明.L2L^2での展開は§E12.25 定理 5.1 (1)による。cj=(u0,ejB)2c_j=(u_0,e_j^B)_2とおく。∑j∣cj∣2=∥u0∥22\sum_j|c_j|^2=\|u_0\|_2^2であり、固定したt>0t>0とm≥0m\geq0に対して

∑j(1+(λjB)2m)e−2tλjB∣cj∣2<∞\sum_j(1+(\lambda_j^B)^{2m})e^{-2t\lambda_j^B}|c_j|^2<\infty

である。実際、λ2me−2tλ\lambda^{2m}e^{-2t\lambda}は[0,∞)[0,\infty)で有界である。固有基底における作用素の表示§E15.19 定理 1.3を反復すると、D(ABm)D(A_B^m)は重み1+(λjB)2m1+(\lambda_j^B)^{2m}による二乗可算性で特徴づけられる。したがって級数の尾部はグラフノルムで零に収束する。▨

定理 2.2.Ω⊂Rn\Omega\subset\R^nを空でない有界 Lipschitz 開集合、B∈{D,N}B\in\{D,N\}、A=ABA=A_Bとし、P0P_0をker⁡A\ker Aへの直交射影、λ+\lambda_+をAAの最初の正固有値とする。このとき任意のu0∈L2(Ω)u_0\in L^2(\Omega)とt≥0t\geq0について

∥e−tAu0−P0u0∥2≤e−λ+t∥(I−P0)u0∥2\|e^{-tA}u_0-P_0u_0\|_2 \leq e^{-\lambda_+t}\|(I-P_0)u_0\|_2

が成り立つ。Dirichlet 条件ではP0=0P_0=0である。Neumann 条件では、各連結成分Ωℓ\Omega_\ell上で

(P0u0)∣Ωℓ=1∣Ωℓ∣∫Ωℓu0(x) dx(P_0u_0)|_{\Omega_\ell}=\frac1{|\Omega_\ell|}\int_{\Omega_\ell}u_0(x)\,dx

である。

証明. 固有値列は非負で無限大へ発散するから、零固有空間は有限次元であり、正固有値の最小値λ+\lambda_+が存在する。定理 2.1と Parseval の等式により

∥e−tAu0−P0u0∥22=∑λj>0e−2tλj∣(u0,ej)2∣2≤e−2tλ+∥(I−P0)u0∥22.\|e^{-tA}u_0-P_0u_0\|_2^2 =\sum_{\lambda_j>0}e^{-2t\lambda_j}|(u_0,e_j)_2|^2 \leq e^{-2t\lambda_+}\|(I-P_0)u_0\|_2^2.

Dirichlet 条件では、§E15.19 定理 2.1によりλ1D>0\lambda_1^D>0であるからker⁡A={0}\ker A=\{0\}であり、P0=0P_0=0である。Neumann 条件では、§E15.19 定理 2.1によりker⁡A\ker Aは各連結成分上で定数となるH1(Ω)H^1(\Omega)の関数全体であり、各成分の定数関数が核をなす。成分の指示関数を正規化した正規直交系への射影を計算すると、主張の平均値の式を得る。▨

例 2.3.Ω=(0,π)\Omega=(0,\pi)では−∂x2sin⁡x=sin⁡x-\partial_x^2\sin x=\sin xであり、sin⁡x\sin xは Dirichlet 条件を満たす。したがって初期値sin⁡x\sin xの解はe−tsin⁡xe^{-t}\sin xである。Neumann 条件では11とcos⁡x\cos xの固有値がそれぞれ00と11であるから、初期値1+cos⁡x1+\cos xの解は1+e−tcos⁡x1+e^{-t}\cos xであり、極限は平均値11である。

非負性を外すと減衰は成立しない。一般の空でない有界 Lipschitz 開集合で、c>λ1Dc>\lambda_1^DとしてA=AD−cIA=A_D-cIとおく。このとき、第一固有関数e1De_1^Dを初期値とする解は

e−tAe1D=e(c−λ1D)te1De^{-tA}e_1^D=e^{(c-\lambda_1^D)t}e_1^D

であり、そのL2L^2ノルムは指数的に増大する。

3 放物型最大値原理

定理 3.1.Ω⊂Rn\Omega\subset\R^nを空でない有界開集合、T>0T>0とし、実関数

u∈C(Ω‾×[0,T])∩C2,1(Ω×(0,T])u\in C(\overline\Omega\times[0,T])\cap C^{2,1}(\Omega\times(0,T])

がut−Δu≤0u_t-\Delta u\leq0を満たすとする。C2,1C^{2,1}の時間微分はTTでは左微分とする。放物型境界を

∂pQ=(Ω‾×{0})∪(∂Ω×[0,T])\partial_pQ=(\overline\Omega\times\{0\})\cup(\partial\Omega\times[0,T])

とおく。このとき、

sup⁡Ω‾×[0,T]u≤sup⁡∂pQu\sup_{\overline\Omega\times[0,T]}u\leq\sup_{\partial_pQ}u

が成り立つ。

証明.ε>0\varepsilon>0とし、v(x,t)=u(x,t)−εtv(x,t)=u(x,t)-\varepsilon tとおく。vt−Δv≤−εv_t-\Delta v\leq-\varepsilonである。コンパクト集合Ω‾×[0,T]\overline\Omega\times[0,T]上のvvの最大点が(x0,t0)∈Ω×(0,T](x_0,t_0)\in\Omega\times(0,T]にあると仮定する。空間変数についてΔv(x0,t0)≤0\Delta v(x_0,t_0)\leq0である。時間変数について、t0<Tt_0<Tならvt(x0,t0)=0v_t(x_0,t_0)=0、t0=Tt_0=Tなら左微分が非負である。したがってvt−Δv≥0v_t-\Delta v\geq0となり、vt−Δv≤−εv_t-\Delta v\leq-\varepsilonに反する。ゆえに最大点は放物型境界にある。v≤sup⁡∂pQuv\leq\sup_{\partial_pQ}uとu=v+εtu=v+\varepsilon tより、u≤sup⁡∂pQu+εTu\leq\sup_{\partial_pQ}u+\varepsilon Tである。ε↓0\varepsilon\downarrow0として結論を得る。▨

系 3.2.T>0T>0とし、u,vu,vをRn×[0,T]\R^n\times[0,T]上の有界な実連続関数で、正時間にはC2,1C^{2,1}に属するものとする。ut−Δu≤vt−Δvu_t-\Delta u\leq v_t-\Delta vかつu(x,0)≤v(x,0)u(x,0)\leq v(x,0)がすべてのxxについて成り立つとする。このとき、u≤vu\leq vである。特に、有界古典解の初期値問題は一意である。非負初期値の有界古典解は非負であり、実初期値u0u_0の有界古典解はsup⁡x,t∣u(x,t)∣≤∥u0∥∞\sup_{x,t}|u(x,t)|\leq\|u_0\|_\inftyを満たす。

証明.w=u−vw=u-vとし、M=sup⁡∣w∣M=\sup|w|とおく。ε>0\varepsilon>0に対し

wε(x,t)=w(x,t)−ε(∣x∣2+2nt+1)w_\varepsilon(x,t)=w(x,t)-\varepsilon(|x|^2+2nt+1)

とおくと、(∂t−Δ)wε≤0(\partial_t-\Delta)w_\varepsilon\leq0である。εR2>M\varepsilon R^2>MとなるRRをとれば、球BRB_Rの側面ではwε<0w_\varepsilon<0、初期面でもwε<0w_\varepsilon<0である。定理 3.1よりBR×[0,T]B_R\times[0,T]でwε≤0w_\varepsilon\leq0となる。RRは任意に大きくすることができるから、全空間でw≤ε(∣x∣2+2nt+1)w\leq\varepsilon(|x|^2+2nt+1)である。各(x,t)(x,t)を固定してε↓0\varepsilon\downarrow0とするとw≤0w\leq0を得る。二つの解の順序を交換すれば一意性が従う。零解との比較から非負性が、定数解±∥u0∥∞\pm\|u_0\|_\inftyとの比較からノルム評価が従う。▨

4 全空間の Gauss 核

定理 4.1.n≥1n\geq1とし、Fourier 変換をf^(ξ)=∫Rnf(x)e−2πix⋅ξ dx\widehat f(\xi)=\int_{\R^n}f(x)e^{-2\pi ix\cdot\xi}\,dxとする。L2(Rn;C)L^2(\R^n;\C)上の非負自己共役作用素A=−ΔA=-\Deltaを、Fourier 変換による乗数4π2∣ξ∣24\pi^2|\xi|^2、定義域{f∈L2:∣ξ∣2f^∈L2}\{f\in L^2:|\xi|^2\widehat f\in L^2\}によって定める。t>0t>0に対して

Gt(x)=(4πt)−n/2e−∣x∣2/(4t)G_t(x)=(4\pi t)^{-n/2}e^{-|x|^2/(4t)}

とおく。このとき

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

が成り立つ。最後の等式はL2L^2の等式であり、畳込みは正時間に滑らかな代表元を与える。

証明.§E15.3 命題 4.1にa=(4πt)−1a=(4\pi t)^{-1}を代入し、係数(4πt)−n/2(4\pi t)^{-n/2}を掛けると、主張の変換式を得る。特にGt≥0G_t\geq0かつ∫Gt=G^t(0)=1\int G_t=\widehat G_t(0)=1である。

f∈S(Rn)f\in\mathcal S(\R^n)のとき、Fubini の定理によりGt∗f^=G^tf^\widehat{G_t*f}=\widehat G_t\widehat fである。Fourier 反転§E15.3 定理 2.1と可測汎関数計算§E12.15 定理 1.2によりGt∗f=e−tAfG_t*f=e^{-tA}fを得る。一般のf∈L2f\in L^2に対してはS\mathcal Sによる近似をとる。Gt dxG_t\,dxは確率測度であるから、Cauchy–Schwarz の不等式と Fubini の定理により

∥Gt∗f∥22≤∫Rn∫RnGt(y)∣f(x−y)∣2 dy dx=∥f∥22.\|G_t*f\|_2^2\leq\int_{\R^n}\int_{\R^n}G_t(y)|f(x-y)|^2\,dy\,dx=\|f\|_2^2.

よって畳込みと半群の両方がL2L^2連続に延長され、等式はL2L^2全体で成り立つ。固定した正時間の任意の時間・空間微分は、GtG_tの微分とffの畳込みである。核の微分はL2L^2に属し、正時間のコンパクト区間ではL2L^2で連続である。Cauchy–Schwarz の不等式により微分積分交換と各微分の連続性が従う。▨

命題 4.2.u0∈Cb(Rn;R)u_0\in C_b(\R^n;\R)とする。u(x,t)=(Gt∗u0)(x)u(x,t)=(G_t*u_0)(x)(t>0t>0)、u(x,0)=u0(x)u(x,0)=u_0(x)とおく。このとき、uuは任意のT>0T>0に対しRn×[0,T]\R^n\times[0,T]上の唯一の有界古典解である。t↓0t\downarrow0のときu(⋅,t)→u0u(\cdot,t)\to u_0は局所一様に成り立つ。u0u_0が一様連続ならば全空間で一様に収束する。

証明.∫Gt=1\int G_t=1により∣u(x,t)∣≤∥u0∥∞|u(x,t)|\leq\|u_0\|_\inftyである。GtG_tの任意の時間・空間微分は、多項式と Gauss 関数の積であり、正時間のコンパクト区間上で共通の可積分な関数により抑えることができる。したがって畳込みを任意階に微分することができる。核を直接微分すると∂tGt=ΔGt\partial_tG_t=\Delta G_tであるからut=Δuu_t=\Delta uを得る。

コンパクト集合K⊂RnK\subset\R^nとδ>0\delta>0をとる。M=∥u0∥∞M=\|u_0\|_\inftyとおくと

sup⁡x∈K∣u(x,t)−u0(x)∣≤sup⁡x∈K∣y∣<δ∣u0(x−y)−u0(x)∣+2M∫∣y∣≥δGt(y) dy.\sup_{x\in K}|u(x,t)-u_0(x)| \leq\sup_{\substack{x\in K\\|y|<\delta}}|u_0(x-y)-u_0(x)| +2M\int_{|y|\geq\delta}G_t(y)\,dy.

第一項はu0u_0のコンパクト集合上の一様連続性によりδ↓0\delta\downarrow0で零となる。第二項の積分は変数変換y=t zy=\sqrt t\,zにより、固定したδ\deltaに対してt↓0t\downarrow0で零となる。δ\deltaを先に選び、次にttを選ぶと局所一様収束を得る。u0u_0が一様連続ならば、同じ評価でKKをRn\R^nに置き換えることができる。初期面での連続性と系 3.2により存在と一意性が従う。▨

例 4.3.

ϕ(t)={e−1/t2t>0,0t≤0,u(x,t)={∑k=0∞ϕ(k)(t)x2k(2k)!t>0,0t≤0\phi(t)=\begin{cases}e^{-1/t^2}&t>0,\\0&t\leq0\end{cases},\qquad u(x,t)=\begin{cases}\displaystyle\sum_{k=0}^\infty\phi^{(k)}(t)\frac{x^{2k}}{(2k)!}&t>0,\\0&t\leq0\end{cases}

とおく。このときu∈C∞(R2;R)u\in C^\infty(\R^2;\R)であり、ut=uxxu_t=u_{xx}、u(x,0)=0u(x,0)=0を満たすが、uuは零関数ではない。したがって、空間無限遠における有界性を課さない古典解の初期値問題には一意性がない。

証明.θ=1/8\theta=1/8とする。t>0t>0と∣z−t∣≤θt|z-t|\leq\theta tに対しz=t(1+w)z=t(1+w)と書くと∣w∣≤1/8|w|\leq1/8であり、

∣(1+w)−2−1∣≤2∣w∣+∣w∣2(1−∣w∣)2≤1749<12.|(1+w)^{-2}-1| \leq\frac{2|w|+|w|^2}{(1-|w|)^2} \leq\frac{17}{49}<\frac12.

したがってℜ(z−2)≥1/(2t2)\Re(z^{-2})\geq1/(2t^2)である。z↦e−1/z2z\mapsto e^{-1/z^2}はこの閉円板の近傍で正則であるから、§E5.6 系 4.2により任意の整数r≥0r\geq0について

∣ϕ(r)(t)∣≤r!(θt)−re−1/(2t2)|\phi^{(r)}(t)|\leq r!(\theta t)^{-r}e^{-1/(2t^2)}

を得る。

整数a,b≥0a,b\geq0とR≥1R\geq1を固定する。各項を時間にaa回、空間にbb回微分した絶対値は、2k≥b2k\geq bのとき

∣ϕ(k+a)(t)∣ ∣x∣2k−b(2k−b)!≤Ca,b(θt)−ae−1/(2t2)(k+1)a+bk!(R2θt)k(∣x∣≤R)\frac{|\phi^{(k+a)}(t)|\,|x|^{2k-b}}{(2k-b)!} \leq C_{a,b}(\theta t)^{-a}e^{-1/(2t^2)} \frac{(k+1)^{a+b}}{k!} \left(\frac{R^2}{\theta t}\right)^k\qquad(|x|\leq R)

である。実際、(k+a)!≤Ca(k+1)ak!(k+a)!\leq C_a(k+1)^ak!、(2k)!/(2k−b)!≤(2k)b(2k)!/(2k-b)!\leq(2k)^b、(2k)!≥(k!)2(2k)!\geq(k!)^2を用いればよい。2k<b2k<bの項の空間微分は零である。

非負整数ddとz≥0z\geq0に対し

∑k=0∞(k+1)dzkk!≤Cd(1+zd)ez\sum_{k=0}^\infty\frac{(k+1)^dz^k}{k!}\leq C_d(1+z^d)e^z

が成り立つ。この不等式は、∑zk/k!=ez\sum z^k/k!=e^zに作用素z∂z+1z\partial_z+1をdd回適用すると、左辺が非負係数多項式とeze^zの積になることから従う。よって微分級数の絶対値の和は

Ca,b(θt)−a(1+(R2θt)a+b)exp⁡(−12t2+R2θt)C_{a,b}(\theta t)^{-a} \left(1+\left(\frac{R^2}{\theta t}\right)^{a+b}\right) \exp\left(-\frac1{2t^2}+\frac{R^2}{\theta t}\right)

以下である。0<δ≤t≤T0<\delta\leq t\leq Tでは、各項の上界の和をt=δt=\deltaの逆冪とeR2/(θδ)e^{R^2/(\theta\delta)}によって抑えることができる。したがって、すべての微分級数は∣x∣≤R|x|\leq R、δ≤t≤T\delta\leq t\leq Tで一様収束し、項別微分が正当化される。

最後の指数評価は、t↓0t\downarrow0ですべてのttの冪より速く零に収束する。したがって正時間のすべての混合微分は、xxのコンパクト集合上一様に零へ収束する。各微分をt≤0t\leq0で零と延長すると連続である。時間微分のt=0t=0における差商も、さらにt−1t^{-1}を掛けた同じ評価によって零となる。空間微分はt=0t=0で恒等的に零である。微分の次数について帰納することによりu∈C∞(R2)u\in C^\infty(\R^2)を得る。

正時間では

ut=∑k=0∞ϕ(k+1)(t)x2k(2k)!,uxx=∑k=1∞ϕ(k)(t)x2k−2(2k−2)!u_t=\sum_{k=0}^\infty\phi^{(k+1)}(t)\frac{x^{2k}}{(2k)!},\qquad u_{xx}=\sum_{k=1}^\infty\phi^{(k)}(t)\frac{x^{2k-2}}{(2k-2)!}

であり、添字をずらすと両者は一致する。非正時間では両者は零である。u(0,t)=ϕ(t)>0u(0,t)=\phi(t)>0(t>0t>0)であるからuuは零解と異なる。系 3.2により、この非零解は任意のT>0T>0に対しR×[0,T]\R\times[0,T]上で有界ではない。▨

5 有界領域の Dirichlet 熱核

定理 5.1.Ω⊂Rn\Omega\subset\R^nを空でない有界C∞C^\infty開集合、ADA_Dを Dirichlet Laplace 作用素、(λj,ej)j∈N≥1(\lambda_j,e_j)_{j\in\NN}をL2(Ω;C)L^2(\Omega;\C)の完全正規直交固有系とする。内積は第一変数について線形とする。このとき

KD(t,x,y)=∑j=1∞e−tλjej(x)ej(y)‾K_D(t,x,y)=\sum_{j=1}^\infty e^{-t\lambda_j}e_j(x)\overline{e_j(y)}

はt>0t>0、x,y∈Ω‾x,y\in\overline\Omegaに対して定義される。任意のε>0\varepsilon>0に対し、この級数は任意階の時間・空間微分とともに[ε,∞)×Ω‾2[\varepsilon,\infty)\times\overline\Omega^2上で一様収束する。任意のf∈L2(Ω)f\in L^2(\Omega)、t>0t>0、x∈Ω‾x\in\overline\Omegaについて

(e−tADf)(x)=∫ΩKD(t,x,y)f(y) dy(e^{-tA_D}f)(x)=\int_\Omega K_D(t,x,y)f(y)\,dy

が成り立つ。左辺には正時間の滑らかな代表元をとる。また、KDK_Dはx∈∂Ωx\in\partial\Omegaまたはy∈∂Ωy\in\partial\Omegaならば零であり、s,t>0s,t>0に対して

KD(t,x,y)=KD(t,y,x)‾,KD(t+s,x,y)=∫ΩKD(t,x,z)KD(s,z,y) dzK_D(t,x,y)=\overline{K_D(t,y,x)},\qquad K_D(t+s,x,y)=\int_\Omega K_D(t,x,z)K_D(s,z,y)\,dz

が成り立つ。

証明.§E15.23 系 6.2により、各eje_jは境界まで滑らかである。任意の整数r≥0r\geq0に対して2Mr>r+n/22M_r>r+n/2を満たす整数MrM_rを選ぶと、同じ評価により、すべてのj∈N≥1j\in\NNについて

∥ej∥Cr(Ω‾)≤Cr(1+λj)Mr\|e_j\|_{C^r(\overline\Omega)}\leq C_r(1+\lambda_j)^{M_r}

を得る。ej∈H01(Ω)e_j\in H_0^1(\Omega)であるから、その境界値は零である。

時間微分の次数aaと空間多重指数α,β\alpha,\betaを固定する。微分後の第jj項の絶対値は、t≥εt\geq\varepsilonで

C(1+λj)Me−ελj≤CM,εe−ελj/2C(1+\lambda_j)^M e^{-\varepsilon\lambda_j} \leq C_{M,\varepsilon}e^{-\varepsilon\lambda_j/2}

以下である。§E15.19 系 3.4によりλj≥cΩj2/n\lambda_j\geq c_\Omega j^{2/n}であるから、右辺はjjについて可算和をもつ。実際、j2/n/log⁡j→∞j^{2/n}/\log j\to\inftyにより、十分大きいjjではe−εcΩj2/n/2≤j−2e^{-\varepsilon c_\Omega j^{2/n}/2}\leq j^{-2}である。Weierstrass の判定法により、すべての微分級数の一様収束を得る。

Ω\Omegaは有限測度でありf∈L1(Ω)f\in L^1(\Omega)であるから、核級数を積分の外へ出すことができる。したがって

∫ΩKD(t,x,y)f(y) dy=∑je−tλj(f,ej)2ej(x).\int_\Omega K_D(t,x,y)f(y)\,dy =\sum_j e^{-t\lambda_j}(f,e_j)_2 e_j(x).

右辺はxxについて一様収束し、定理 2.1によりL2L^2ではe−tADfe^{-tA_D}fに等しい。両辺の連続性により、等式は境界を含めて点ごとに成り立つ。

境界での消失と共役対称性は各項について成り立ち、一様収束によって和へ移る。半群則については二つの級数の積をzzについて積分する。両級数は絶対一様収束するから、積分と二重和を交換することができる。∫Ωej(z)‾ek(z) dz=δjk\int_\Omega\overline{e_j(z)}e_k(z)\,dz=\delta_{jk}により二重和は

∑je−(t+s)λjej(x)ej(y)‾=KD(t+s,x,y)\sum_j e^{-(t+s)\lambda_j}e_j(x)\overline{e_j(y)}=K_D(t+s,x,y)

となる。▨

前提記事