Lemma数学ロードマップ

18 偏微分方程式論 · 第 5 章

発展方程式

目安 10〜12 時間定理など 9演習 6 問

この章の目標

  • バナッハ空間値関数の空間 Lp(0,T;X)L^p(0, T; X) と、時間についての弱微分を理解する
  • ガラーキン法により放物型方程式の弱解を構成し、エネルギー評価と一意性を示せる
  • 熱方程式を作用素半群として定式化し、ヒレ–吉田の定理の仮定を確かめられる
  • 放物型方程式の最大値原理、波動方程式のエネルギー保存と有限伝播速度を証明できる
  • シュレディンガー方程式のユニタリ群・分散評価・ストリッカーツ評価の意味を説明できる

前提:第2章〜第4章、10-functional-analysis 第7章(ストーンの定理・ヒレ–吉田の定理)、08-differential-equations 第3章・第7章

5.1 関数空間に値をとる関数

熱方程式 ut=Δuu_t = \Delta u を、時刻 tt ごとの「空間の関数」u(t)=u(⋅,t)u(t) = u(\cdot, t) の時間発展と見れば、関数空間 XX における常微分方程式 u′(t)=Δu(t)u'(t) = \Delta u(t) になる。Δ\Delta は非有界作用素なので、常微分方程式の理論(08-differential-equations 第3章)はそのままでは使えない。本章では二つの方法を扱う。

  1. エネルギー法とガラーキン法:有限次元の常微分方程式で近似し、エネルギー評価と弱コンパクト性で極限をとる。係数が時間に依存する問題や非線形問題にも使える。
  2. 作用素半群:u(t)=etΔu(0)u(t) = e^{t\Delta}u(0) を関数解析的に構成する(ヒレ–吉田の定理)。

どちらの方法でも、第4章の楕円型理論が土台になる。まず、関数空間に値をとる関数の空間を用意する。

定義 5.1(ボホナー空間)XX をバナッハ空間、T>0T > 0 とする。単関数の a.e. 極限として表される関数 u ⁣:[0,T]→Xu\colon[0, T] \to X(強可測関数)で、∥u∥Lp(0,T;X)=(∫0T∥u(t)∥Xpdt)1/p<∞\lVert u\rVert_{L^p(0,T;X)} = \bigl(\int_0^T\lVert u(t)\rVert_X^pdt\bigr)^{1/p} < \infty(p=∞p = \infty では ess sup⁡t∥u(t)∥X<∞\operatorname{ess\ sup}_t\lVert u(t)\rVert_X < \infty)となるものの全体を Lp(0,T;X)L^p(0, T; X) と書く。u,v∈L1(0,T;X)u, v \in L^1(0, T; X) について、すべての φ∈Cc∞(0,T)\varphi \in C_c^\infty(0, T) に対し ∫0Tφ′(t)u(t) dt=−∫0Tφ(t)v(t) dt\int_0^T\varphi'(t)u(t)\ dt = -\int_0^T\varphi(t)v(t)\ dt(XX におけるボホナー積分)が成り立つとき、vv を uu の弱導関数といい u′=vu' = v と書く。

XX がヒルベルト空間なら L2(0,T;X)L^2(0, T; X) もヒルベルト空間である。放物型方程式ではゲルファントの三つ組 H01(Ω)⊂L2(Ω)⊂H−1(Ω)H_0^1(\Omega) \subset L^2(\Omega) \subset H^{-1}(\Omega) を用いる:f∈L2(Ω)f \in L^2(\Omega) を汎関数 v↦(f,v)v \mapsto (f, v) とみて H−1H^{-1} の元と同一視し、H−1H^{-1} と H01H_0^1 の組を ⟨⋅,⋅⟩\langle\cdot, \cdot\rangle と書く。

定理 5.2 (1) u∈Lp(0,T;X)u \in L^p(0, T; X)、u′∈Lp(0,T;X)u' \in L^p(0, T; X)(1≤p≤∞1 \leq p \leq \infty)ならば、uu は C([0,T];X)C([0, T]; X) の元と a.e. で一致し、u(t)=u(s)+∫stu′(τ) dτu(t) = u(s) + \int_s^tu'(\tau)\ d\tau。 (2)(リオンス–マジェネス)Ω\Omega を有界開集合とし、u∈L2(0,T;H01(Ω))u \in L^2(0, T; H_0^1(\Omega))、u′∈L2(0,T;H−1(Ω))u' \in L^2(0, T; H^{-1}(\Omega)) とする。このとき uu は C([0,T];L2(Ω))C([0, T]; L^2(\Omega)) の元と a.e. で一致し、t↦∥u(t)∥L22t \mapsto \lVert u(t)\rVert_{L^2}^2 は絶対連続で、a.e. の tt で ddt∥u(t)∥L22=2⟨u′(t),u(t)⟩\frac{d}{dt}\lVert u(t)\rVert_{L^2}^2 = 2\langle u'(t), u(t)\rangle。

証明は Evans の第5章を参照。(2) は、u(t)u(t) 自身は H01H_0^1 にしか、u′(t)u'(t) は H−1H^{-1} にしかなくても、両者の組が「ちょうど釣り合って」L2L^2 ノルムの微分が計算できることを述べている。

5.2 放物型方程式の弱解:ガラーキン法

Ω\Omega を有界開集合、LL を第4章の仮定 (A) をみたす作用素(係数は tt によらない)、BB を対応する双線形形式とし、

ut+Lu=f  (Ω×(0,T)),u=0  (∂Ω×[0,T]),u=g  (t=0)(5.1)u_t + Lu = f \ \ (\Omega \times (0, T)), \qquad u = 0 \ \ (\partial\Omega \times [0, T]), \qquad u = g \ \ (t = 0) \tag{5.1}

を考える。データは f∈L2(0,T;L2(Ω))f \in L^2(0, T; L^2(\Omega))、g∈L2(Ω)g \in L^2(\Omega) とする。滑らかな解に v∈H01v \in H_0^1 を掛けて xx について部分積分すると (ut,v)+B[u,v]=(f,v)(u_t, v) + B[u, v] = (f, v) となる。各時刻で u(t)∈H01u(t) \in H_0^1 ならば B[u(t),v]B[u(t), v] が意味をもち、ut=f−Luu_t = f - Lu は H−1H^{-1} に属するのが自然である。

定義 5.3(弱解)u∈L2(0,T;H01(Ω))u \in L^2(0, T; H_0^1(\Omega))、u′∈L2(0,T;H−1(Ω))u' \in L^2(0, T; H^{-1}(\Omega)) が

  1. a.e. の t∈[0,T]t \in [0, T] とすべての v∈H01(Ω)v \in H_0^1(\Omega) について ⟨u′(t),v⟩+B[u(t),v]=(f(t),v)\langle u'(t), v\rangle + B[u(t), v] = (f(t), v)、
  2. u(0)=gu(0) = g(定理 5.2 (2) により u∈C([0,T];L2)u \in C([0, T]; L^2) なので意味をもつ)

をみたすとき、uu を (5.1) の弱解という。

定理 5.4(弱解の存在と一意性)(5.1) の弱解はただ一つ存在し、第4章の定数 α,γ\alpha, \gamma(定理 4.5)について

max⁡0≤t≤T∥u(t)∥2≤e(2γ+1)T(∥g∥2+∫0T∥f∥2dt),θ∫0T∥u∥H12dt≤e(2γ+1)T(∥g∥2+∫0T∥f∥2dt)\max_{0 \leq t \leq T}\lVert u(t)\rVert^2 \leq e^{(2\gamma+1)T}\Bigl(\lVert g\rVert^2 + \int_0^T\lVert f\rVert^2dt\Bigr), \qquad \theta\int_0^T\lVert u\rVert_{H^1}^2dt \leq e^{(2\gamma+1)T}\Bigl(\lVert g\rVert^2 + \int_0^T\lVert f\rVert^2dt\Bigr)

および ∥u′(t)∥H−1≤∥f(t)∥+α∥u(t)∥H1\lVert u'(t)\rVert_{H^{-1}} \leq \lVert f(t)\rVert + \alpha\lVert u(t)\rVert_{H^1}(a.e. tt)が成り立つ。ここで ∥⋅∥\lVert\cdot\rVert は L2(Ω)L^2(\Omega) のノルムである。

証明(詳しい概略). (i) ガラーキン近似. wkw_k を −Δ-\Delta のディリクレ固有関数(定理 4.10)とする。{wk}\lbrace w_k\rbrace は L2(Ω)L^2(\Omega) の正規直交基底で、(∇wk,∇wl)=λkδkl(\nabla w_k, \nabla w_l) = \lambda_k\delta_{kl}。um(t)=∑k=1mdmk(t)wku_m(t) = \sum_{k=1}^md_m^k(t)w_k を

(um′,wk)+B[um,wk]=(f(t),wk),dmk(0)=(g,wk)(k=1,…,m)(u_m', w_k) + B[u_m, w_k] = (f(t), w_k), \qquad d_m^k(0) = (g, w_k) \qquad (k = 1, \dots, m)

で定める。これは定数係数の線形系 d′+Md=F(t)d' + Md = F(t)(Mkl=B[wl,wk]M_{kl} = B[w_l, w_k]、Fk=(f(t),wk)∈L2(0,T)F_k = (f(t), w_k) \in L^2(0, T))であり、d(t)=e−tMd(0)+∫0te−(t−s)MF(s) dsd(t) = e^{-tM}d(0) + \int_0^te^{-(t-s)M}F(s)\ ds がただ一つの(絶対連続な)解である。

(ii) エネルギー評価. kk 番目の式に dmkd_m^k を掛けて和をとると 12ddt∥um∥2+B[um,um]=(f,um)\frac{1}{2}\frac{d}{dt}\lVert u_m\rVert^2 + B[u_m, u_m] = (f, u_m)。定理 4.5 の B[u,u]≥θ2∥u∥H12−γ∥u∥2B[u, u] \geq \frac{\theta}{2}\lVert u\rVert_{H^1}^2 - \gamma\lVert u\rVert^2 と 2(f,um)≤∥f∥2+∥um∥22(f, u_m) \leq \lVert f\rVert^2 + \lVert u_m\rVert^2 から

ddt∥um∥2+θ∥um∥H12≤κ∥um∥2+∥f∥2(κ=2γ+1)\frac{d}{dt}\lVert u_m\rVert^2 + \theta\lVert u_m\rVert_{H^1}^2 \leq \kappa\lVert u_m\rVert^2 + \lVert f\rVert^2 \qquad (\kappa = 2\gamma + 1)

e−κte^{-\kappa t} を掛けて 00 から tt まで積分し、∥um(0)∥≤∥g∥\lVert u_m(0)\rVert \leq \lVert g\rVert(ベッセルの不等式)を使うと、∥um(t)∥2+θ∫0t∥um∥H12≤eκT(∥g∥2+∫0T∥f∥2)\lVert u_m(t)\rVert^2 + \theta\int_0^t\lVert u_m\rVert_{H^1}^2 \leq e^{\kappa T}(\lVert g\rVert^2 + \int_0^T\lVert f\rVert^2)(グロンウォールの不等式)。時間微分については、v∈H01v \in H_0^1 と Pmv=∑k≤m(v,wk)wkP_mv = \sum_{k \leq m}(v, w_k)w_k について ∥Pmv∥≤∥v∥\lVert P_mv\rVert \leq \lVert v\rVert、∥∇Pmv∥2=∑k≤mλk(v,wk)2≤∥∇v∥2\lVert\nabla P_mv\rVert^2 = \sum_{k \leq m}\lambda_k(v, w_k)^2 \leq \lVert\nabla v\rVert^2(定理 4.10)だから ∥Pmv∥H1≤∥v∥H1\lVert P_mv\rVert_{H^1} \leq \lVert v\rVert_{H^1} で、

(um′,v)=(um′,Pmv)=(f,Pmv)−B[um,Pmv],∥um′∥H−1≤∥f∥+α∥um∥H1(u_m', v) = (u_m', P_mv) = (f, P_mv) - B[u_m, P_mv], \qquad \lVert u_m'\rVert_{H^{-1}} \leq \lVert f\rVert + \alpha\lVert u_m\rVert_{H^1}

固有関数を基底に選んだのはこの評価のためである。以上から (um)(u_m) は L2(0,T;H01)L^2(0, T; H_0^1) で、(um′)(u_m') は L2(0,T;H−1)L^2(0, T; H^{-1}) で有界である。

(iii) 極限. ヒルベルト空間の有界列なので、部分列について um⇀uu_m \rightharpoonup u(L2(0,T;H01)L^2(0, T; H_0^1))、um′⇀u♯u_m' \rightharpoonup u^\sharp(L2(0,T;H−1)L^2(0, T; H^{-1}))と弱収束し、弱導関数の定義で極限をとれば u♯=u′u^\sharp = u'。NN を固定し、v(t)=∑k=1Nαk(t)wkv(t) = \sum_{k=1}^N\alpha_k(t)w_k(αk∈C1[0,T]\alpha_k \in C^1[0, T])とすると、m≥Nm \geq N で ∫0T(⟨um′,v⟩+B[um,v]) dt=∫0T(f,v) dt\int_0^T(\langle u_m', v\rangle + B[u_m, v])\ dt = \int_0^T(f, v)\ dt。弱極限をとると uu についても同じ式が成り立ち、このような vv は L2(0,T;H01)L^2(0, T; H_0^1) で稠密なので、定義 5.3 の 1 が a.e. の tt で成り立つ。v(T)=0v(T) = 0 となる vv で時間について部分積分し、um(0)→gu_m(0) \to g(L2L^2)と比べると u(0)=gu(0) = g を得る。評価は弱極限でも保たれる(ノルムの弱下半連続性と、定理 5.2 (2) による uu 自身のエネルギー等式)。

(iv) 一意性. f=0f = 0、g=0g = 0 の弱解 uu について、定理 5.2 (2) と v=u(t)v = u(t) から ddt∥u∥2=−2B[u,u]≤2γ∥u∥2\frac{d}{dt}\lVert u\rVert^2 = -2B[u, u] \leq 2\gamma\lVert u\rVert^2。∥u(0)∥=0\lVert u(0)\rVert = 0 とグロンウォールの不等式より u≡0u \equiv 0。□\square

注意 5.5(正則性)∂Ω\partial\Omega が C2C^2 級、aij∈C1(Ω‾)a^{ij} \in C^1(\overline\Omega)、g∈H01(Ω)g \in H_0^1(\Omega) ならば、弱解は u∈L2(0,T;H2(Ω))∩L∞(0,T;H01(Ω))u \in L^2(0, T; H^2(\Omega)) \cap L^\infty(0, T; H_0^1(\Omega))、u′∈L2(0,T;L2(Ω))u' \in L^2(0, T; L^2(\Omega)) をみたす(第4章の境界正則性を各時刻に使う。Evans の第7章)。データが滑らかで適合条件をみたせば解は滑らかである。さらに放物型方程式には、gg が L2L^2 でしかなくても t>0t > 0 で解が滑らかになる平滑化効果がある(命題 5.6)。

5.3 半群による定式化

簡単のため L=−ΔL = -\Delta とし、Ω\Omega は有界で境界が C2C^2 級とする。X=L2(Ω)X = L^2(\Omega) 上の作用素 AA を

D(A)=H2(Ω)∩H01(Ω),Au=ΔuD(A) = H^2(\Omega) \cap H_0^1(\Omega), \qquad Au = \Delta u

で定める。境界条件は定義域 D(A)D(A) に組み込まれている。

命題 5.6(熱半群)AA は L2(Ω)L^2(\Omega) 上の縮小 C0C_0 半群 S(t)=etΔS(t) = e^{t\Delta} を生成する。−Δ-\Delta のディリクレ固有値・固有関数を λk\lambda_k, wkw_k(定理 4.10)とすると

S(t)g=∑k=1∞e−λkt(g,wk)wk,∥S(t)g∥≤e−λ1t∥g∥,∥ΔS(t)g∥≤1et∥g∥(t>0)S(t)g = \sum_{k=1}^\infty e^{-\lambda_kt}(g, w_k)w_k, \qquad \lVert S(t)g\rVert \leq e^{-\lambda_1t}\lVert g\rVert, \qquad \lVert\Delta S(t)g\rVert \leq \frac{1}{et}\lVert g\rVert \quad (t > 0)

証明. ヒレ–吉田の定理(10-functional-analysis 第7章 定理 7.18)の仮定を確かめる。D(A)⊃Cc∞(Ω)D(A) \supset C_c^\infty(\Omega) は稠密。λ>0\lambda > 0、f∈L2f \in L^2 に対し、λu−Δu=f\lambda u - \Delta u = f、u∈H01u \in H_0^1 はただ一つの弱解をもち(系 4.7 で c=λc = \lambda)、境界正則性(定理 4.16)により u∈H2u \in H^2、すなわち u∈D(A)u \in D(A) で (λ−A)u=f(\lambda - A)u = f。よって λ−A ⁣:D(A)→L2\lambda - A\colon D(A) \to L^2 は全単射で、v=uv = u を代入した λ∥u∥2+∥∇u∥2=(f,u)\lambda\lVert u\rVert^2 + \lVert\nabla u\rVert^2 = (f, u) から ∥(λ−A)−1f∥≤∥f∥/λ\lVert(\lambda - A)^{-1}f\rVert \leq \lVert f\rVert/\lambda。有界な逆をもつので AA は閉作用素である。したがって AA は縮小半群を生成する。wk∈D(A)w_k \in D(A)、Awk=−λkwkAw_k = -\lambda_kw_k なので e−λktwke^{-\lambda_kt}w_k は u′=Auu' = Au、u(0)=wku(0) = w_k の解であり、半群による解の一意性から S(t)wk=e−λktwkS(t)w_k = e^{-\lambda_kt}w_k。線形性と連続性から展開式を得る。AA は(グリーンの公式により)D(A)D(A) 上で対称なので、u∈D(A)u \in D(A) について (Au,wk)=−λk(u,wk)(Au, w_k) = -\lambda_k(u, w_k) であり

∥ΔS(t)g∥2=∑kλk2e−2λkt(g,wk)2≤(sup⁡λ>0λe−λt)2∥g∥2=1(et)2∥g∥2\lVert\Delta S(t)g\rVert^2 = \sum_k\lambda_k^2e^{-2\lambda_kt}(g, w_k)^2 \leq \Bigl(\sup_{\lambda > 0}\lambda e^{-\lambda t}\Bigr)^2\lVert g\rVert^2 = \frac{1}{(et)^2}\lVert g\rVert^2

第 2 の評価も同様である。□\square

∥ΔS(t)g∥≤(et)−1∥g∥\lVert\Delta S(t)g\rVert \leq (et)^{-1}\lVert g\rVert は平滑化効果の定量的な表現である。同様に ∥ΔkS(t)g∥≤(k/(et))k∥g∥\lVert\Delta^kS(t)g\rVert \leq (k/(et))^k\lVert g\rVert となり、∂Ω\partial\Omega が C∞C^\infty 級なら、境界正則性(定理 4.16)を繰り返し使って t>0t > 0 で S(t)g∈C∞(Ω‾)S(t)g \in C^\infty(\overline\Omega) がわかる。−A-A は非負の自己共役作用素なので、系 7.19(同第7章)から直接に半群の生成を結論してもよい。非斉次方程式 u′=Au+fu' = Au + f の解はデュアメルの公式 u(t)=S(t)g+∫0tS(t−s)f(s) dsu(t) = S(t)g + \int_0^tS(t - s)f(s)\ ds で与えられ、これは定理 5.4 の弱解と一致する。

5.4 放物型方程式の最大値原理

Ω\Omega を有界開集合、ΩT=Ω×(0,T]\Omega_T = \Omega \times (0, T] とし、放物境界 ΓT=ΩT‾∖ΩT=(Ω‾×{0})∪(∂Ω×[0,T])\Gamma_T = \overline{\Omega_T} \setminus \Omega_T = (\overline\Omega \times \lbrace 0\rbrace) \cup (\partial\Omega \times [0, T]) とおく。ΩT\Omega_T で u,∂iu,∂i∂ju,∂tuu, \partial_iu, \partial_i\partial_ju, \partial_tu が連続な関数の全体を C2,1(ΩT)C^{2,1}(\Omega_T) と書く。

定理 5.7(放物型の弱最大値原理)u∈C2,1(ΩT)∩C(ΩT‾)u \in C^{2,1}(\Omega_T) \cap C(\overline{\Omega_T}) が ΩT\Omega_T で ut−Δu≤0u_t - \Delta u \leq 0 をみたすならば、max⁡ΩT‾u=max⁡ΓTu\max_{\overline{\Omega_T}}u = \max_{\Gamma_T}u。

証明. まず ut−Δu<0u_t - \Delta u < 0 と仮定する。uu が最大値を (x0,t0)∈ΩT(x_0, t_0) \in \Omega_T でとったとすると、x0x_0 は u(⋅,t0)u(\cdot, t_0) の内部の最大点なので Δu(x0,t0)≤0\Delta u(x_0, t_0) \leq 0。また 0<t0≤T0 < t_0 \leq T で、t↦u(x0,t)t \mapsto u(x_0, t) は t0t_0 で(t0=Tt_0 = T なら左側から)最大なので ut(x0,t0)≥0u_t(x_0, t_0) \geq 0。よって ut−Δu≥0u_t - \Delta u \geq 0 となり矛盾。一般の場合は ε>0\varepsilon > 0 について v=u−εtv = u - \varepsilon t とすると vt−Δv≤−ε<0v_t - \Delta v \leq -\varepsilon < 0 なので、max⁡ΩT‾u≤max⁡ΩT‾v+εT=max⁡ΓTv+εT≤max⁡ΓTu+εT\max_{\overline{\Omega_T}}u \leq \max_{\overline{\Omega_T}}v + \varepsilon T = \max_{\Gamma_T}v + \varepsilon T \leq \max_{\Gamma_T}u + \varepsilon T。ε→0\varepsilon \to 0 とすればよい。□\square

系 5.8 u∈C2,1(ΩT)∩C(ΩT‾)u \in C^{2,1}(\Omega_T) \cap C(\overline{\Omega_T}) が ut−Δu=fu_t - \Delta u = f(ff は有界)をみたすならば max⁡ΩT‾∣u∣≤max⁡ΓT∣u∣+Tsup⁡∣f∣\max_{\overline{\Omega_T}}\lvert u\rvert \leq \max_{\Gamma_T}\lvert u\rvert + T\sup\lvert f\rvert。特に、初期値・境界値問題の古典解は一意であり、データに一様ノルムで連続に依存する。

証明. v=±u−tsup⁡∣f∣v = \pm u - t\sup\lvert f\rvert は vt−Δv≤0v_t - \Delta v \leq 0 をみたし、ΓT\Gamma_T 上で v≤max⁡ΓT∣u∣v \leq \max_{\Gamma_T}\lvert u\rvert。定理 5.7 を適用すればよい。□\square

1 次元の場合は 08-differential-equations 第7章 で扱った。さらに強最大値原理が成り立つ:Ω\Omega が連結で、uu が最大値を (x0,t0)∈ΩT(x_0, t_0) \in \Omega_T でとるならば、uu は Ω‾×[0,t0]\overline\Omega \times [0, t_0] で定数である(主張のみ;Evans の第2章・第7章)。したがって、g≥0g \geq 0、g≢0g \not\equiv 0、境界値 00 の熱方程式の解は、すべての t>0t > 0 で Ω\Omega 全体において正になる。熱は瞬時に全体に伝わる(無限の伝播速度)。

5.5 波動方程式

Ω\Omega を有界で境界が C1C^1 級の開集合とし、utt−Δu=fu_{tt} - \Delta u = f を考える。エネルギーを

E(t)=12∫Ω(ut2+∣∇u∣2)dxE(t) = \frac{1}{2}\int_\Omega\bigl(u_t^2 + \lvert\nabla u\rvert^2\bigr)dx

で定める(運動エネルギーと位置エネルギーの和)。

定理 5.9(エネルギー等式と一意性)u∈C2(Ω‾×[0,T])u \in C^2(\overline\Omega \times [0, T]) が utt−Δu=fu_{tt} - \Delta u = f(Ω×(0,T)\Omega \times (0, T) で)、u=0u = 0(∂Ω×[0,T]\partial\Omega \times [0, T] で)をみたすならば、E′(t)=∫Ωfut dxE'(t) = \int_\Omega fu_t\ dx であり、

E(t)1/2≤E(0)1/2+12∫0t∥f(s)∥ ds(0≤t≤T)E(t)^{1/2} \leq E(0)^{1/2} + \frac{1}{\sqrt{2}}\int_0^t\lVert f(s)\rVert\,ds \qquad (0 \leq t \leq T)

特に f=0f = 0 ならば E(t)=E(0)E(t) = E(0)(エネルギー保存)。初期値・境界値問題の古典解は一意である。

証明. 積分記号下で微分し、ガウス–グリーンの公式を使うと

E′(t)=∫Ω(ututt+∇u⋅∇ut) dx=∫Ωut(utt−Δu) dx+∫∂Ωut ∂νu dS=∫Ωfut dxE'(t) = \int_\Omega(u_tu_{tt} + \nabla u \cdot \nabla u_t)\,dx = \int_\Omega u_t(u_{tt} - \Delta u)\,dx + \int_{\partial\Omega}u_t\,\partial_\nu u\,dS = \int_\Omega fu_t\,dx

(境界で u=0u = 0 だから ut=0u_t = 0)。E′≤∥f∥∥ut∥≤∥f∥2EE' \leq \lVert f\rVert\lVert u_t\rVert \leq \lVert f\rVert\sqrt{2E} より、ε>0\varepsilon > 0 について ddtE+ε=E′2E+ε≤∥f∥2\frac{d}{dt}\sqrt{E + \varepsilon} = \frac{E'}{2\sqrt{E + \varepsilon}} \leq \frac{\lVert f\rVert}{\sqrt{2}}。積分して ε→0\varepsilon \to 0 とすれば評価を得る。一意性:2 つの解の差は f=0f = 0、初期値 00 の解なので E≡0E \equiv 0、よって ut=0u_t = 0、∇u=0\nabla u = 0 で、u≡0u \equiv 0。□\square

同じエネルギー評価を用いて、ガラーキン法により弱解が構成できる:g∈H01(Ω)g \in H_0^1(\Omega)、h∈L2(Ω)h \in L^2(\Omega)、f∈L2(0,T;L2(Ω))f \in L^2(0, T; L^2(\Omega)) に対し、初期条件 u(0)=gu(0) = g、ut(0)=hu_t(0) = h の弱解 u∈C([0,T];H01)∩C1([0,T];L2)u \in C([0, T]; H_0^1) \cap C^1([0, T]; L^2) がただ一つ存在し、エネルギー等式をみたす(証明は Evans の第7章)。半群の言葉では、U=(u,ut)U = (u, u_t) をエネルギー空間 H01×L2H_0^1 \times L^2 上の方程式 U′=AUU' = \mathcal{A}U、A(u,v)=(v,Δu)\mathcal{A}(u, v) = (v, \Delta u) とみると、A\mathcal{A} は歪自己共役で、etAe^{t\mathcal{A}} はエネルギーを保つ群になる(ストーンの定理)。熱方程式と違い、波動方程式には平滑化効果がなく、時間を逆向きにも解ける。

定理 5.10(有限伝播速度)u∈C2(Rn×[0,∞))u \in C^2(\mathbb{R}^n \times [0, \infty)) が utt−Δu=0u_{tt} - \Delta u = 0 をみたし、x0∈Rnx_0 \in \mathbb{R}^n、t0>0t_0 > 0 について B(x0,t0)B(x_0, t_0) 上で u(⋅,0)=ut(⋅,0)=0u(\cdot, 0) = u_t(\cdot, 0) = 0 とする。このとき円錐 K={(x,t)∣0≤t≤t0, ∣x−x0∣≤t0−t}K = \lbrace(x, t) \mid 0 \leq t \leq t_0,\ \lvert x - x_0\rvert \leq t_0 - t\rbrace 上で u=0u = 0。

証明. 0≤t<t00 \leq t < t_0 に対し、縮んでいく球の上の局所エネルギー e(t)=12∫B(x0,t0−t)(ut2+∣∇u∣2) dxe(t) = \frac{1}{2}\int_{B(x_0, t_0 - t)}(u_t^2 + \lvert\nabla u\rvert^2)\ dx を考える。半径が速さ 11 で減るので

e′(t)=∫B(x0,t0−t)(ututt+∇u⋅∇ut) dx−12∫∂B(x0,t0−t)(ut2+∣∇u∣2) dSe'(t) = \int_{B(x_0, t_0 - t)}(u_tu_{tt} + \nabla u \cdot \nabla u_t)\,dx - \frac{1}{2}\int_{\partial B(x_0, t_0 - t)}(u_t^2 + \lvert\nabla u\rvert^2)\,dS

第 1 項は部分積分により ∫∂But ∂νu dS\int_{\partial B}u_t\ \partial_\nu u\ dS に等しく(utt−Δu=0u_{tt} - \Delta u = 0)、∣ut ∂νu∣≤12(ut2+∣∇u∣2)\lvert u_t\ \partial_\nu u\rvert \leq \frac{1}{2}(u_t^2 + \lvert\nabla u\rvert^2) だから e′(t)≤0e'(t) \leq 0。e(0)=0e(0) = 0 なので [0,t0)[0, t_0) で e≡0e \equiv 0、すなわち KK(頂点を除く、連続性により頂点も含む)で ut=0u_t = 0、∇u=0\nabla u = 0。KK は連結で、t=0t = 0 で u=0u = 0 だから、KK 上で u=0u = 0。□\square

点 (x0,t0)(x_0, t_0) での解の値は、初期データの B(x0,t0)‾\overline{B(x_0, t_0)} 上の値だけで決まる(依存領域)。擾乱は速さ 11 を超えて伝わらない。これは第1章で見た波動作用素の基本解の台が円錐 ∣x∣≤t\lvert x\rvert \leq t に含まれることに対応する。フーリエ変換では解は u^=cos⁡(t∣ξ∣)g^+sin⁡(t∣ξ∣)∣ξ∣h^\hat u = \cos(t\lvert\xi\rvert)\hat g + \frac{\sin(t\lvert\xi\rvert)}{\lvert\xi\rvert}\hat h と表され、これらの乗数が指数型 tt の整関数に延長されることが、ペイリー–ウィーナー–シュワルツの定理(定理 2.18)を通じて有限伝播速度を与える。

5.6 シュレディンガー方程式

自由シュレディンガー方程式 i∂tu=−Δui\partial_tu = -\Delta u(x∈Rnx \in \mathbb{R}^n)を考える。H=−ΔH = -\Delta、D(H)=H2(Rn)D(H) = H^2(\mathbb{R}^n) はフーリエ変換により ∣ξ∣2\lvert\xi\rvert^2 を掛ける作用素にユニタリ同値なので、L2(Rn)L^2(\mathbb{R}^n) 上の非負の自己共役作用素である(実数値関数を掛ける作用素は、その最大定義域で自己共役)。ストーンの定理(10-functional-analysis 第7章 定理 7.16)により、U(t)=e−itHU(t) = e^{-itH} は強連続なユニタリ群で、g∈H2g \in H^2 なら u(t)=U(t)gu(t) = U(t)g は iu′=Huiu' = Hu、u(0)=gu(0) = g の解である。第2章の計算と合わせて次を得る。

命題 5.11(シュレディンガー群)U(t)g^(ξ)=e−it∣ξ∣2g^(ξ)\widehat{U(t)g}(\xi) = e^{-it\lvert\xi\rvert^2}\hat g(\xi) であり、

  1. すべての s∈Rs \in \mathbb{R} について ∥U(t)g∥Hs=∥g∥Hs\lVert U(t)g\rVert_{H^s} = \lVert g\rVert_{H^s}(特に L2L^2 ノルムの保存)。
  2. (分散評価)g∈L1∩L2g \in L^1 \cap L^2、t≠0t \neq 0 ならば sup⁡x∣U(t)g(x)∣≤(4π∣t∣)−n/2∥g∥L1\sup_x\lvert U(t)g(x)\rvert \leq (4\pi\lvert t\rvert)^{-n/2}\lVert g\rVert_{L^1}。
  3. 2≤p≤∞2 \leq p \leq \infty、1p+1p′=1\frac{1}{p} + \frac{1}{p'} = 1 ならば ∥U(t)g∥Lp≤(4π∣t∣)−n(12−1p)∥g∥Lp′\lVert U(t)g\rVert_{L^p} \leq (4\pi\lvert t\rvert)^{-n(\frac{1}{2} - \frac{1}{p})}\lVert g\rVert_{L^{p'}}。

証明. (1) ∣e−it∣ξ∣2∣=1\lvert e^{-it\lvert\xi\rvert^2}\rvert = 1 と HsH^s ノルムの定義による。(2) g∈Sg \in \mathcal{S} なら命題 2.20 であり、一般の場合は L1∩L2L^1 \cap L^2 での近似による((1) により L2L^2 で、したがって部分列が a.e. で収束する)。(3) は p=2p = 2((1))と p=∞p = \infty((2))の間のリース–トーリンの補間定理による(補間定理は本教材では証明しない)。□\square

熱方程式と違い、シュレディンガー方程式は HsH^s ノルムを保つので平滑化も減衰もしない(L2L^2 の意味で)。しかし (2)(3) の意味で解は空間に広がって各点で減衰する。この分散性を時空の積分ノルムで表したものがストリッカーツ評価である。

定理 5.12(ストリッカーツ評価, Strichartz estimates)2≤q,r≤∞2 \leq q, r \leq \infty、2q+nr=n2\frac{2}{q} + \frac{n}{r} = \frac{n}{2}、(q,r,n)≠(2,∞,2)(q, r, n) \neq (2, \infty, 2) とする(許容対)。定数 C=C(n,q,r)C = C(n, q, r) があって、すべての g∈L2(Rn)g \in L^2(\mathbb{R}^n) について

(∫R∥U(t)g∥Lr(Rn)q dt)1/q≤C∥g∥L2(Rn)\Bigl(\int_{\mathbb{R}}\lVert U(t)g\rVert_{L^r(\mathbb{R}^n)}^q\,dt\Bigr)^{1/q} \leq C\lVert g\rVert_{L^2(\mathbb{R}^n)}

(主張のみ。)指数の関係は、u(x,t)u(x, t) が解なら u(λx,λ2t)u(\lambda x, \lambda^2t) も解であるというスケーリングから必要になる(問題 5.6)。(q,r)=(∞,2)(q, r) = (\infty, 2) は L2L^2 ノルムの保存そのものである。q>2q > 2 の場合の証明は次の TT∗TT^{\ast} 論法による。Tg=U(⋅)gTg = U(\cdot)g とおくと ∥T∥2=∥TT∗∥\lVert T\rVert^2 = \lVert TT^{\ast}\rVert で、TT∗F(t)=∫U(t−s)F(s) dsTT^{\ast}F(t) = \int U(t - s)F(s)\ ds。命題 5.11 (3) により ∥U(t−s)F(s)∥Lr≤C∣t−s∣−2/q∥F(s)∥Lr′\lVert U(t - s)F(s)\rVert_{L^r} \leq C\lvert t - s\rvert^{-2/q}\lVert F(s)\rVert_{L^{r'}}(n(12−1r)=2qn(\frac{1}{2} - \frac{1}{r}) = \frac{2}{q})であり、時間変数についての ∣t∣−2/q\lvert t\rvert^{-2/q} との畳み込みはハーディ–リトルウッド–ソボレフの不等式により Lq′→LqL^{q'} \to L^q で有界である(0<2q<10 < \frac{2}{q} < 1 が必要)。端点 q=2q = 2(n≥3n \geq 3)はキール–タオ(1998)による。ストリッカーツ評価は、非線形シュレディンガー方程式 iut+Δu=±∣u∣p−1uiu_t + \Delta u = \pm\lvert u\rvert^{p-1}u の解を縮小写像の原理で構成するための基本的な道具である(第6章)。

まとめ

  • 発展方程式は関数空間に値をとる関数の常微分方程式とみなせる。弱解は u∈L2(0,T;H01)u \in L^2(0, T; H_0^1)、u′∈L2(0,T;H−1)u' \in L^2(0, T; H^{-1}) の枠で定式化され、u∈C([0,T];L2)u \in C([0, T]; L^2) となる(リオンス–マジェネスの定理)。
  • ガラーキン法:有限次元近似 → エネルギー評価(グロンウォール)→ 弱コンパクト性で極限、という手順で放物型方程式の弱解を構成した。一意性もエネルギー法による。
  • 熱方程式は、境界条件を定義域に組み込んだ A=ΔA = \Delta が縮小半群を生成すること(ヒレ–吉田の定理の仮定は第4章の楕円型理論で確かめられる)として捉えられ、∥ΔetΔg∥≤(et)−1∥g∥\lVert\Delta e^{t\Delta}g\rVert \leq (et)^{-1}\lVert g\rVert という平滑化効果をもつ。
  • 放物型方程式は最大値原理をみたし、熱は無限の速さで伝わる。
  • 波動方程式はエネルギーを保存し、擾乱は有限の速さでしか伝わらない。平滑化効果はない。
  • シュレディンガー方程式はユニタリ群 eitΔe^{it\Delta} で解かれ、HsH^s ノルムを保ちつつ ∣t∣−n/2\lvert t\rvert^{-n/2} で各点減衰する。これはストリッカーツ評価に一般化される。

演習問題

問題 5.1 ★ Ω\Omega を有界で境界が C1C^1 級とし、u∈C2(Ω‾×[0,T])u \in C^2(\overline\Omega \times [0, T]) が ut=Δuu_t = \Delta u、境界で u=0u = 0 をみたすとする。∥u(t)∥L2\lVert u(t)\rVert_{L^2} と ∥∇u(t)∥L2\lVert\nabla u(t)\rVert_{L^2} はともに tt について単調非増加であることを示せ。さらに Ω=(0,π)\Omega = (0, \pi) では ∥u(t)∥L2≤e−t∥u(0)∥L2\lVert u(t)\rVert_{L^2} \leq e^{-t}\lVert u(0)\rVert_{L^2} を示せ。

解答

12ddt∥u∥2=∫uΔu=−∥∇u∥2≤0\frac{1}{2}\frac{d}{dt}\lVert u\rVert^2 = \int u\Delta u = -\lVert\nabla u\rVert^2 \leq 0(境界項は u=0u = 0 で消える)。12ddt∥∇u∥2=∫∇u⋅∇ut=−∫Δu ut+∫∂Ωut ∂νu dS=−∫ut2≤0\frac{1}{2}\frac{d}{dt}\lVert\nabla u\rVert^2 = \int\nabla u \cdot \nabla u_t = -\int\Delta u\ u_t + \int_{\partial\Omega}u_t\ \partial_\nu u\ dS = -\int u_t^2 \leq 0(境界で ut=0u_t = 0)。Ω=(0,π)\Omega = (0, \pi) では問題 3.5 より ∥ux∥2≥∥u∥2\lVert u_x\rVert^2 \geq \lVert u\rVert^2 なので ddt∥u∥2≤−2∥u∥2\frac{d}{dt}\lVert u\rVert^2 \leq -2\lVert u\rVert^2、よって ∥u(t)∥2≤e−2t∥u(0)∥2\lVert u(t)\rVert^2 \leq e^{-2t}\lVert u(0)\rVert^2。減衰率 11 は第 1 固有値 λ1=1\lambda_1 = 1 である。

問題 5.2 ★ 系 5.8 の仮定のもとで、ut−Δu≤vt−Δvu_t - \Delta u \leq v_t - \Delta v(ΩT\Omega_T で)かつ ΓT\Gamma_T 上で u≤vu \leq v ならば ΩT‾\overline{\Omega_T} で u≤vu \leq v であることを示せ。これを用いて、0≤g≤10 \leq g \leq 1 を初期値、境界値 00 とする熱方程式の古典解は 0≤u≤10 \leq u \leq 1 をみたすことを示せ。

解答

w=u−vw = u - v は wt−Δw≤0w_t - \Delta w \leq 0 で ΓT\Gamma_T 上 w≤0w \leq 0 なので、定理 5.7 より w≤0w \leq 0。後半は v=1v = 1(vt−Δv=0v_t - \Delta v = 0、ΓT\Gamma_T 上 u≤1u \leq 1)と比べて u≤1u \leq 1、00 と比べて(−u-u に適用して)u≥0u \geq 0。

問題 5.3 ★★ 命題 5.6 の半群について、∥∇S(t)g∥≤(2et)−1/2∥g∥\lVert\nabla S(t)g\rVert \leq (2et)^{-1/2}\lVert g\rVert(t>0t > 0)を示せ。

解答

S(t)g∈H01S(t)g \in H_0^1 で、定理 4.10 の公式 ∥∇u∥2=∑kλk(u,wk)2\lVert\nabla u\rVert^2 = \sum_k\lambda_k(u, w_k)^2(u∈H01u \in H_0^1)を u=S(t)gu = S(t)g、(u,wk)=e−λkt(g,wk)(u, w_k) = e^{-\lambda_kt}(g, w_k) に使うと ∥∇S(t)g∥2=∑kλke−2λkt(g,wk)2≤sup⁡λ>0λe−2λt ∥g∥2\lVert\nabla S(t)g\rVert^2 = \sum_k\lambda_ke^{-2\lambda_kt}(g, w_k)^2 \leq \sup_{\lambda > 0}\lambda e^{-2\lambda t}\ \lVert g\rVert^2。λe−2λt\lambda e^{-2\lambda t} は λ=12t\lambda = \frac{1}{2t} で最大値 12et\frac{1}{2et} をとる。

問題 5.4 ★★ m≥0m \geq 0 とし、u∈C2(Rn×[0,∞))u \in C^2(\mathbb{R}^n \times [0, \infty)) がクライン–ゴルドン方程式 utt−Δu+m2u=0u_{tt} - \Delta u + m^2u = 0 をみたすとする。定理 5.10 の結論が同じ仮定のもとで成り立つことを示せ。また、u(⋅,t)u(\cdot, t) が各 tt でコンパクト台をもつとき、12∫(ut2+∣∇u∣2+m2u2) dx\frac{1}{2}\int(u_t^2 + \lvert\nabla u\rvert^2 + m^2u^2)\ dx が保存されることを示せ。

解答

e(t)=12∫B(x0,t0−t)(ut2+∣∇u∣2+m2u2) dxe(t) = \frac{1}{2}\int_{B(x_0, t_0 - t)}(u_t^2 + \lvert\nabla u\rvert^2 + m^2u^2)\ dx とおくと、定理 5.10 の証明と同様に

e′(t)=∫∂But ∂νu dS−12∫∂B(ut2+∣∇u∣2+m2u2) dS≤0e'(t) = \int_{\partial B}u_t\,\partial_\nu u\,dS - \frac{1}{2}\int_{\partial B}(u_t^2 + \lvert\nabla u\rvert^2 + m^2u^2)\,dS \leq 0

(内部の項は ut(utt−Δu+m2u)=0u_t(u_{tt} - \Delta u + m^2u) = 0 となって消える)。e(0)=0e(0) = 0 から KK 上 ut=∇u=0u_t = \nabla u = 0、u=0u = 0。保存則は、十分大きな球の上で同じ計算をすれば境界項が消えることによる。

問題 5.5 ★★(エネルギーの等分配)g,h∈Cc∞(R)g, h \in C_c^\infty(\mathbb{R}) とし、uu を utt=uxxu_{tt} = u_{xx}、u(x,0)=gu(x, 0) = g、ut(x,0)=hu_t(x, 0) = h の解とする。十分大きな tt について ∫ut2 dx=∫ux2 dx\int u_t^2\ dx = \int u_x^2\ dx であることを示せ。

解答

ダランベールの公式(08-differential-equations 第7章)より u=F(x+t)+G(x−t)u = F(x + t) + G(x - t)、F′=12(g′+h)F' = \frac{1}{2}(g' + h)、G′=12(g′−h)G' = \frac{1}{2}(g' - h) で、F′,G′F', G' はコンパクト台をもつ。ut=F′(x+t)−G′(x−t)u_t = F'(x + t) - G'(x - t)、ux=F′(x+t)+G′(x−t)u_x = F'(x + t) + G'(x - t) なので ∫ut2−∫ux2=−4∫F′(x+t)G′(x−t) dx\int u_t^2 - \int u_x^2 = -4\int F'(x + t)G'(x - t)\ dx。F′(⋅+t)F'(\cdot + t) の台は左に、G′(⋅−t)G'(\cdot - t) の台は右に速さ 11 で動くので、tt が大きければ交わらず、積分は 00 になる。

問題 5.6 ★★★ (1) 不等式 ∥U(t)g∥LtqLxr≤C∥g∥L2\lVert U(t)g\rVert_{L^q_tL^r_x} \leq C\lVert g\rVert_{L^2} が成り立つためには 2q+nr=n2\frac{2}{q} + \frac{n}{r} = \frac{n}{2} が必要であることを示せ。(2) 分散評価を用いて、すべての g∈L2(Rn)g \in L^2(\mathbb{R}^n) について t→∞t \to \infty で U(t)gU(t)g は L2L^2 で 00 に弱収束することを示せ(∥U(t)g∥L2\lVert U(t)g\rVert_{L^2} は一定であることに注意)。

解答

(1) gλ(x)=g(λx)g_\lambda(x) = g(\lambda x) とすると、g^λ(ξ)=λ−ng^(ξ/λ)\hat g_\lambda(\xi) = \lambda^{-n}\hat g(\xi/\lambda) から U(t)gλ(x)=(U(λ2t)g)(λx)U(t)g_\lambda(x) = (U(\lambda^2t)g)(\lambda x)。よって ∥U(t)gλ∥LtqLxr=λ−2q−nr∥U(t)g∥LtqLxr\lVert U(t)g_\lambda\rVert_{L^q_tL^r_x} = \lambda^{-\frac{2}{q} - \frac{n}{r}}\lVert U(t)g\rVert_{L^q_tL^r_x}、一方 ∥gλ∥2=λ−n/2∥g∥2\lVert g_\lambda\rVert_2 = \lambda^{-n/2}\lVert g\rVert_2。λ→0,∞\lambda \to 0, \infty を考えれば、指数が等しくなければ不等式は成り立たない。

(2) φ,g∈L1∩L2\varphi, g \in L^1 \cap L^2 ならば ∣(U(t)g,φ)∣≤sup⁡∣U(t)g∣ ∥φ∥1≤(4πt)−n/2∥g∥1∥φ∥1→0\lvert(U(t)g, \varphi)\rvert \leq \sup\lvert U(t)g\rvert\ \lVert\varphi\rVert_1 \leq (4\pi t)^{-n/2}\lVert g\rVert_1\lVert\varphi\rVert_1 \to 0。一般の g,φ∈L2g, \varphi \in L^2 については、ε>0\varepsilon > 0 に対し L1∩L2L^1 \cap L^2 の元 g′,φ′g', \varphi' を ∥g−g′∥2,∥φ−φ′∥2<ε\lVert g - g'\rVert_2, \lVert\varphi - \varphi'\rVert_2 < \varepsilon となるようにとる(Cc∞C_c^\infty は L2L^2 で稠密)。U(t)U(t) はユニタリなので ∣(U(t)g,φ)−(U(t)g′,φ′)∣≤ε(∥φ∥2+∥g′∥2)\lvert(U(t)g, \varphi) - (U(t)g', \varphi')\rvert \leq \varepsilon(\lVert\varphi\rVert_2 + \lVert g'\rVert_2) で、lim sup⁡t∣(U(t)g,φ)∣≤ε(∥φ∥2+∥g∥2+ε)\limsup_t\lvert(U(t)g, \varphi)\rvert \leq \varepsilon(\lVert\varphi\rVert_2 + \lVert g\rVert_2 + \varepsilon)。ε→0\varepsilon \to 0 とすればよい。質量は保存されるが、どの有界領域からも逃げていく。

この章を読み終えたら

「読了」にすると学習記録とロードマップに反映されます。演習の自己採点もお忘れなく。

この章の誤りを報告GitHub で見る