§E15.27波動方程式

最終更新

両端を固定した弦の振動では、初期変位と初速度から、その後の変位を求めることが問題となる。例えば、区間(0,π)(0,\pi)における初期変位sin⁡x\sin x、初速度零の振動は、波動方程式の解u(t,x)=cos⁡tsin⁡xu(t,x)=\cos t\sin xで表される。波動方程式は双曲型方程式の基本的な例であり、速度と変位の空間微分から定まるエネルギーが解の大きさを制御するが、その総量だけでは初期データの影響が届く範囲までは分からない。エネルギー解はこの量が有限であるという条件に基づく解の概念であり、非負な楕円型作用素の余弦・正弦作用素によって構成される。大域的なエネルギーとその空間的な分布を併せて考えることで、解の安定性と有限伝播速度を共通の枠組みで捉えることができる。本記事では、波動方程式のエネルギー解とその基本的な性質について解説する。

1 余弦・正弦作用素

定義 1.1 (余弦作用素と正弦作用素).HHを複素 Hilbert 空間、AAをHH上の非負自己共役作用素とする。内積は第一変数について線形とする。t∈Rt\in\Rとλ≥0\lambda\ge0に対して

ct(λ)=cos⁡(tλ),st(λ)={sin⁡(tλ)/λ,λ>0,t,λ=0c_t(\lambda)=\cos(t\sqrt\lambda),\qquad s_t(\lambda)= \begin{cases} \sin(t\sqrt\lambda)/\sqrt\lambda,&\lambda>0,\\ t,&\lambda=0 \end{cases}

とおき、可測汎関数計算により

CA(t)=ct(A),SA(t)=st(A)C_A(t)=c_t(A),\qquad S_A(t)=s_t(A)

と定める。CA(t)C_A(t)を 余弦作用素 (cosine operator)、SA(t)S_A(t)を 正弦作用素 (sine operator) という。V=D(A1/2)V=D(A^{1/2})には

∥v∥V2=∥v∥H2+∥A1/2v∥H2\|v\|_V^2=\|v\|_H^2+\|A^{1/2}v\|_H^2

を入れる。

命題 1.2.AAを複素 Hilbert 空間HH上の非負自己共役作用素、V=D(A1/2)V=D(A^{1/2})とする。このとき、任意のt∈Rt\in\Rに対して

∥CA(t)∥H→H≤1,∥SA(t)∥H→H≤∣t∣,∥A1/2SA(t)∥H→H≤1\|C_A(t)\|_{H\to H}\le1,\qquad \|S_A(t)\|_{H\to H}\le|t|,\qquad \|A^{1/2}S_A(t)\|_{H\to H}\le1

が成り立つ。SA(t):H→VS_A(t):H\to VとCA(t):V→VC_A(t):V\to Vは有界であり、

∥SA(t)h∥V≤1+t2∥h∥H,∥CA(t)v∥V≤∥v∥V\|S_A(t)h\|_V\le\sqrt{1+t^2}\|h\|_H,\qquad \|C_A(t)v\|_V\le\|v\|_V

が成り立つ。両作用素はHH上で強連続であり、SA(t)hS_A(t)hはVV値でも連続、CA(t)vC_A(t)vはVV値で連続である。さらに、HHにおける強微分として

ddtSA(t)h=CA(t)h(h∈H),ddtCA(t)v=−A1/2sin⁡(tA)v(v∈V)\frac{d}{dt}S_A(t)h=C_A(t)h\quad(h\in H),\qquad \frac{d}{dt}C_A(t)v=-A^{1/2}\sin(t\sqrt A)v\quad(v\in V)

が成り立つ。

証明.B=A1/2B=A^{1/2}とおく。スカラー評価

∣ct(λ)∣≤1,∣st(λ)∣≤∣t∣,∣λst(λ)∣=∣sin⁡(tλ)∣≤1|c_t(\lambda)|\le1,\quad |s_t(\lambda)|\le|t|,\quad |\sqrt\lambda s_t(\lambda)|=|\sin(t\sqrt\lambda)|\le1

と§E12.15 定理 1.2の積の定義域から、SA(t)H⊂D(B)S_A(t)H\subset D(B)と主張のノルム評価を得る。v∈Vv\in Vに対してBCA(t)v=CA(t)BvBC_A(t)v=C_A(t)Bvであるから、CA(t)C_A(t)はVV上でも縮小作用素である。

EAE_AをAAのスペクトル測度とし、μh(M)=(EA(M)h,h)H\mu_h(M)=(E_A(M)h,h)_Hとおく。汎関数計算のノルム公式

∥m(A)h∥H2=∫[0,∞)∣m(λ)∣2 dμh(λ)\|m(A)h\|_H^2=\int_{[0,\infty)}|m(\lambda)|^2\,d\mu_h(\lambda)

を用いる。有限な時間区間ではct,st,λstc_t,s_t,\sqrt\lambda s_tはそれぞれ11、時間区間の長さを含む定数、11で支配される。点ごとの連続性と優収束により、CA(t)hC_A(t)hとSA(t)hS_A(t)hはHH値で連続であり、BSA(t)hBS_A(t)hも連続である。BCA(t)v=CA(t)BvBC_A(t)v=C_A(t)Bvに同じ議論を適用すると、VV値の連続性を得る。

sts_tの時間差分商は絶対値が11以下であり、ctc_tに点ごとに収束する。ctc_tの時間差分商は絶対値がλ\sqrt\lambda以下であり、−λsin⁡(tλ)-\sqrt\lambda\sin(t\sqrt\lambda)に点ごとに収束する。前者には有限測度μh\mu_hを、後者には∫λ dμv<∞\int\lambda\,d\mu_v<\inftyを用いて優収束を適用すると、二つの強微分公式を得る。▨

2 エネルギー解

定義 2.1 (エネルギー解).AAを複素 Hilbert 空間HH上の非負自己共役作用素、V=D(A1/2)V=D(A^{1/2})、T>0T>0とする。V†V^\daggerをVVの連続反線形汎関数の空間とし、HHをh↦(v↦(h,v)H)h\mapsto(v\mapsto(h,v)_H)によりV†V^\daggerに埋め込む。A:V→V†\mathcal A:V\to V^\daggerを

⟨Au,v⟩=(A1/2u,A1/2v)H\langle\mathcal Au,v\rangle=(A^{1/2}u,A^{1/2}v)_H

で定める。u0∈Vu_0\in V、u1∈Hu_1\in H、f∈L1(0,T;H)f\in L^1(0,T;H)に対して、

u∈C([0,T];V)∩C1([0,T];H),u(0)=u0,ut(0)=u1u\in C([0,T];V)\cap C^1([0,T];H),\qquad u(0)=u_0,\quad u_t(0)=u_1

を満たし、u′′+Au=fu''+\mathcal Au=fがV†V^\dagger値の超関数として成立するuuを エネルギー解 (energy solution) という。すなわち、任意のv∈Vv\in Vとη∈Cc∞(0,T)\eta\in C_c^\infty(0,T)に対して

∫0T[(u(t),v)Hη′′(t)+(A1/2u(t),A1/2v)Hη(t)]dt=∫0T(f(t),v)Hη(t) dt\int_0^T\bigl[(u(t),v)_H\eta''(t) +(A^{1/2}u(t),A^{1/2}v)_H\eta(t)\bigr]dt =\int_0^T(f(t),v)_H\eta(t)\,dt

が成り立つ。時間の超関数と試験関数の組合せには共役を入れない。

定理 2.2.AAを複素 Hilbert 空間HH上の非負自己共役作用素、V=D(A1/2)V=D(A^{1/2})、T>0T>0とする。このとき、任意のu0∈Vu_0\in V、u1∈Hu_1\in H、f∈L1(0,T;H)f\in L^1(0,T;H)に対して、唯一のエネルギー解が存在し、

u(t)=CA(t)u0+SA(t)u1+∫0tSA(t−s)f(s) dsu(t)=C_A(t)u_0+S_A(t)u_1+\int_0^tS_A(t-s)f(s)\,ds

で与えられる。積分はVV値の Bochner 積分である。初期値と外力から解への写像は

V×H×L1(0,T;H)⟶C([0,T];V)∩C1([0,T];H)V\times H\times L^1(0,T;H) \longrightarrow C([0,T];V)\cap C^1([0,T];H)

の間の連続線形写像である。さらにu0∈D(A)u_0\in D(A)、u1∈Vu_1\in V、f∈C([0,T];V)f\in C([0,T];V)ならば

u∈C([0,T];D(A))∩C1([0,T];V)∩C2([0,T];H)u\in C([0,T];D(A))\cap C^1([0,T];V)\cap C^2([0,T];H)

が成り立ち、u′′(t)+Au(t)=f(t)u''(t)+Au(t)=f(t)がすべてのt∈[0,T]t\in[0,T]で成立する。D(A)D(A)にはグラフノルムを入れる。

証明.B=A1/2B=A^{1/2}とおき、外力項をw(t)w(t)と書く。命題 1.2と§E10.18 命題 4.3により

∥w(t)∥V≤1+T2∫0t∥f(s)∥H ds\|w(t)\|_V\le\sqrt{1+T^2}\int_0^t\|f(s)\|_H\,ds

である。強連続な作用素族を単関数近似に適用すると、積分の被積分関数は強可測である。時間差に対する優収束§E10.18 定理 5.1によりwwはVV値で連続である。SA(0)=0S_A(0)=0と差分商の評価から

w′(t)=∫0tCA(t−s)f(s) dsw'(t)=\int_0^tC_A(t-s)f(s)\,ds

を得る。実際、積分区間が共通の部分ではSAS_Aの差分商を∥f(s)∥H\|f(s)\|_Hで支配して優収束を適用する。新たに加わる長さ∣h∣|h|の部分の差分商のノルムは、その部分の∥f(s)∥H\|f(s)\|_Hの積分以下であり、h→0h\to0で零に収束する。w′w'の連続性もCAC_Aの強連続性と同じ優収束から従う。したがって公式のuuは指定の連続性と初期値を満たし、

ut(t)=−Bsin⁡(tA)u0+CA(t)u1+∫0tCA(t−s)f(s) dsu_t(t)=-B\sin(t\sqrt A)u_0+C_A(t)u_1 +\int_0^tC_A(t-s)f(s)\,ds

となる。

PN=1[0,N](A)P_N=1_{[0,N]}(A)とおく。PNP_NはH,VH,V上の縮小射影であり、HN=PNHH_N=P_NH上のAN=A∣HNA_N=A|_{H_N}は有界である。有界作用素のべき級数を微分すると、uN=PNuu_N=P_Nuは

uN′′+ANuN=PNfu_N''+A_Nu_N=P_Nf

をHNH_N値の超関数として満たす。積分とPNP_Nの交換は§E10.18 命題 5.2による。スペクトルノルム公式からPNui→uiP_Nu_i\to u_iがそれぞれV,HV,Hで成り立ち、優収束からPNf→fP_Nf\to fがL1(0,T;H)L^1(0,T;H)で成り立つ。伝播作用素の評価によりuN→uu_N\to uはCtVC_tV、uN′→u′u_N'\to u'はCtHC_tHで成立する。A:V→V†\mathcal A:V\to V^\daggerは有界であるから、試験関数と組み合わせた方程式の極限を取ることによりu′′+Au=fu''+\mathcal Au=fを得る。

zzを零初期値・零外力のエネルギー解とする。v∈HNv\in H_NはVVに属するから、弱方程式をvvで試すとzN=PNzz_N=P_NzはzN′′+ANzN=0z_N''+A_Nz_N=0を満たす。右辺−ANzN-A_Nz_Nは連続であり、超関数の方程式を二回積分して

zN(t)=−∫0t(t−s)ANzN(s) dsz_N(t)=-\int_0^t(t-s)A_Nz_N(s)\,ds

を得る。長さδ\deltaがδ2∥AN∥/2<1\delta^2\|A_N\|/2<1を満たす時間区間では、この式の最大ノルム評価からzN=0z_N=0である。零初期値を各区間の左端へ移して有限回繰り返すとzN=0z_N=0が[0,T][0,T]上で成立する。PN→IP_N\to Iは強収束するからz=0z=0であり、一意性を得る。

伝播作用素の評価と上のutu_tの公式から

sup⁡0≤t≤T(∥u(t)∥V+∥ut(t)∥H)≤CT(∥u0∥V+∥u1∥H+∥f∥Lt1H)\sup_{0\le t\le T}(\|u(t)\|_V+\|u_t(t)\|_H) \le C_T(\|u_0\|_V+\|u_1\|_H+\|f\|_{L^1_tH})

を得る。線形性により同じ評価を二つのデータの差へ適用することができるので、連続依存性が従う。

追加の正則性の仮定の下では、汎関数計算によりSA(t):V→D(A)S_A(t):V\to D(A)は有界であり、

ASA(t)v=sin⁡(tA)Bv,ACA(t)u0=CA(t)Au0AS_A(t)v=\sin(t\sqrt A)Bv,\qquad AC_A(t)u_0=C_A(t)Au_0

である。これらの式と優収束からu∈CtD(A)u\in C_tD(A)を得る。CA(t)C_A(t)はVV上で強連続であり、Bsin⁡(tA):D(A)→VB\sin(t\sqrt A):D(A)\to Vも強連続であるから、utu_tの公式はCtVC_tVに属する。方程式の右辺f−Auf-AuはHH値で連続であり、u′′=f−Auu''=f-Auからu∈Ct2Hu\in C_t^2Hを得る。▨

定理 2.3.AAを複素 Hilbert 空間HH上の非負自己共役作用素、V=D(A1/2)V=D(A^{1/2})、T>0T>0とし、uuを初期値(u0,u1)∈V×H(u_0,u_1)\in V\times H、外力f∈L1(0,T;H)f\in L^1(0,T;H)のエネルギー解とする。このとき、関数

E(t)=12(∥ut(t)∥H2+∥A1/2u(t)∥H2)E(t)=\frac12\bigl(\|u_t(t)\|_H^2+\|A^{1/2}u(t)\|_H^2\bigr)

は絶対連続であり、任意の0≤t≤T0\le t\le Tについて

E(t)−E(0)=ℜ∫0t(f(s),ut(s))H ds,E(t)-E(0)=\Re\int_0^t(f(s),u_t(s))_H\,ds,2E(t)≤2E(0)+∫0t∥f(s)∥H ds,∥u(t)∥H≤∥u0∥H+∫0t2E(s) ds\sqrt{2E(t)}\le\sqrt{2E(0)}+\int_0^t\|f(s)\|_H\,ds,\qquad \|u(t)\|_H\le\|u_0\|_H+\int_0^t\sqrt{2E(s)}\,ds

が成り立つ。特にf=0f=0ならばE(t)=E(0)E(t)=E(0)である。

証明.PN=1[0,N](A)P_N=1_{[0,N]}(A)、uN=PNuu_N=P_Nu、fN=PNff_N=P_Nfとし、ENE_NをEEの式でuuをuNu_Nに置き換えた量とする。有界なANA_Nの方程式からuN′u_N'は絶対連続である。内積を微分し、自己共役性を用いると、ほとんどすべてのttについて

EN′(t)=ℜ(fN(t)−ANuN(t),uN′(t))H+ℜ(ANuN(t),uN′(t))H=ℜ(fN(t),uN′(t))HE_N'(t)=\Re(f_N(t)-A_Nu_N(t),u_N'(t))_H +\Re(A_Nu_N(t),u_N'(t))_H =\Re(f_N(t),u_N'(t))_H

となる。積分後にN→∞N\to\inftyとする。定理 2.2の証明の収束からEN→EE_N\to Eは一様収束し、fN→ff_N\to fはLt1HL^1_tH、uN′→u′u_N'\to u'はCtHC_tHで収束する。したがって積分恒等式が成立し、その右辺が絶対連続であるからEEも絶対連続である。

ε>0\varepsilon>0に対し、ほとんどすべてのttで

ddt2E(t)+ε=ℜ(f(t),ut(t))H2E(t)+ε≤∥f(t)∥H\frac{d}{dt}\sqrt{2E(t)+\varepsilon} =\frac{\Re(f(t),u_t(t))_H}{\sqrt{2E(t)+\varepsilon}} \le\|f(t)\|_H

である。積分してε↓0\varepsilon\downarrow0とすると第一の評価を得る。u(t)=u0+∫0tut(s) dsu(t)=u_0+\int_0^tu_t(s)\,dsに§E10.18 命題 4.3を適用すると第二の評価を得る。f=0f=0のとき保存則は積分恒等式から従う。▨

例 2.4.Ω\Omegaを空でない有界 Lipschitz 開集合とし、Dirichlet Laplace 作用素の第一固有値をλ1D\lambda_1^D、対応するL2L^2ノルム11の固有関数をe1e_1とし、c>λ1Dc>\lambda_1^D、A=−ΔD−cIA=-\Delta_D-cI、α=c−λ1D\alpha=\sqrt{c-\lambda_1^D}とする。このとき、関数

u(t)=cosh⁡(αt)e1u(t)=\cosh(\alpha t)e_1

はutt+Au=0u_{tt}+Au=0と初期値(e1,0)(e_1,0)を満たす。固有関数の存在は§E15.19 定理 1.3による。実際、Ae1=−α2e1Ae_1=-\alpha^2e_1であり、utt=α2uu_{tt}=\alpha^2uである。形式的な保存量

12(∥ut∥22+(Au,u)2)=−α22\frac12\bigl(\|u_t\|_2^2+(Au,u)_2\bigr)=-\frac{\alpha^2}{2}

は正ではなく、指数的に増大する変位と速度を制御しない。非負性は定理 2.3の正のエネルギー評価に必要な仮定である。

3 固有関数展開

定理 3.1.Ω⊂Rn\Omega\subset\R^nを空でない有界 Lipschitz 開集合とし、AD=−ΔDA_D=-\Delta_D、AN=−ΔNA_N=-\Delta_Nをそれぞれ斉次 Dirichlet、斉次 Neumann 実現とする。B∈{D,N}B\in\{D,N\}に対して

VD=H01(Ω),VN=H1(Ω)V_D=H_0^1(\Omega),\qquad V_N=H^1(\Omega)

とおく。T>0T>0、(u0,u1)∈VB×L2(Ω)(u_0,u_1)\in V_B\times L^2(\Omega)、f∈L1(0,T;L2(Ω))f\in L^1(0,T;L^2(\Omega))とする。§E15.19 定理 1.3の完全正規直交固有系を(ejB)(e_j^B)、固有値を(λjB)(\lambda_j^B)とし、uij=(ui,ejB)2u_{ij}=(u_i,e_j^B)_2、fj(s)=(f(s),ejB)2f_j(s)=(f(s),e_j^B)_2と書く。このときエネルギー解は

u(t)=∑j=1∞uj(t)ejB,u(t)=\sum_{j=1}^\infty u_j(t)e_j^B,uj(t)=cos⁡(tλjB)u0j+st(λjB)u1j+∫0tst−s(λjB)fj(s) dsu_j(t)=\cos(t\sqrt{\lambda_j^B})u_{0j} +s_t(\lambda_j^B)u_{1j} +\int_0^t s_{t-s}(\lambda_j^B)f_j(s)\,ds

で与えられる。部分和はC([0,T];VB)C([0,T];V_B)で、時間微分した部分和はC([0,T];L2)C([0,T];L^2)で収束する。u0∈D(AB)u_0\in D(A_B)、u1∈VBu_1\in V_B、f∈C([0,T];VB)f\in C([0,T];V_B)ならば、さらにCtD(AB)C_tD(A_B)、Ct1VBC_t^1V_B、Ct2L2C_t^2L^2の各位相で収束する。P0P_0をker⁡AN\ker A_Nへの直交射影とすると、Neumann 解について

P0u(t)=P0u0+tP0u1+∫0t(t−s)P0f(s) dsP_0u(t)=P_0u_0+tP_0u_1+\int_0^t(t-s)P_0f(s)\,ds

が成り立つ。

証明. 形式q(u,v)=∫Ω∇u⋅∇v‾ dxq(u,v)=\int_\Omega\nabla u\cdot\overline{\nabla v}\,dxに§E15.24 命題 1.1を適用すると、D(AB1/2)=VBD(A_B^{1/2})=V_Bであり、グラフノルムはH1H^1ノルムと同値である。定理 2.2の公式をejBe_j^Bと内積に取ると、汎関数計算と§E10.18 命題 5.2により各係数の公式を得る。

QmQ_mを最初のmm個の固有関数への射影とする。§E15.19 定理 1.3によりQmu0→u0Q_mu_0\to u_0はVBV_Bで、Qmu1→u1Q_mu_1\to u_1はL2L^2で成立する。優収束によりQmf→fQ_mf\to fはLt1L2L^1_tL^2で成立する。連続依存性から、対応する解QmuQ_muは主張のエネルギー位相で収束する。

追加の仮定の下ではQmu0→u0Q_mu_0\to u_0はD(AB)D(A_B)、Qmu1→u1Q_mu_1\to u_1はVBV_Bで成立する。f([0,T])f([0,T])はVBV_Bのコンパクト集合である。QmQ_mはVBV_B上で一様有界かつ強収束するから、有限被覆による近似によりQmf→fQ_mf\to fはCtVBC_tV_Bでも成立する。定理 2.2の正則性の証明に現れた作用素評価をデータの差へ適用すると、追加の三つの収束を得る。零固有空間ではct(0)=1c_t(0)=1、st(0)=ts_t(0)=tであるから、同じ公式にP0P_0を適用して最後の式を得る。▨

例 3.2.Ω=(0,π)\Omega=(0,\pi)とする。このとき、Dirichlet 条件で初期値が(sin⁡x,0)(\sin x,0)、外力が零ならば、−∂x2sin⁡x=sin⁡x-\partial_x^2\sin x=\sin xより

u(t,x)=cos⁡tsin⁡xu(t,x)=\cos t\sin x

である。Neumann 条件で初期値が(0,1)(0,1)、外力が零ならば、定数関数は零固有関数であるから

u(t,x)=tu(t,x)=t

である。後者のエネルギーはE(t)=π/2E(t)=\pi/2で一定であるが、∥u(t)∥2=π t\|u(t)\|_2=\sqrt\pi\,tは無限に増大する。st(0)s_t(0)を零と置くと、後者の公式の初速度が零になり、初期条件を満たさない。

4 局所エネルギー

命題 4.1.Ω⊂Rn\Omega\subset\R^nを空でない有界C∞C^\infty開集合とし、a∈C∞(Ω‾;Rn×n)a\in C^\infty(\overline\Omega;\R^{n\times n})を時間に依存しない実対称行列とする。あるθ,c>0\theta,c>0に対して

θ∣ξ∣2≤a(x)ξ⋅ξ‾≤c2∣ξ∣2(x∈Ω‾, ξ∈Cn)\theta|\xi|^2\le a(x)\xi\cdot\overline\xi\le c^2|\xi|^2 \quad(x\in\overline\Omega,\ \xi\in\C^n)

とする。V=H01(Ω)V=H_0^1(\Omega)またはV=H1(Ω)V=H^1(\Omega)とし、

q(v,w)=∫Ω(a∇v)⋅∇w‾ dxq(v,w)=\int_\Omega(a\nabla v)\cdot\overline{\nabla w}\,dx

に対応する斉次 Dirichlet または conormal Neumann 実現をAAとする。このときAAは非負自己共役であり、D(A1/2)=VD(A^{1/2})=V、q(v,w)=(A1/2v,A1/2w)2q(v,w)=(A^{1/2}v,A^{1/2}w)_2が成り立つ。したがって定理 2.2と定理 2.3がこの実現に適用される。全空間Rn\R^nのA=−ΔA=-\Deltaでは同じ結論がV=H1(Rn)V=H^1(\R^n)について成り立つ。

証明. 形式qqは対称かつ非負であり、楕円性と上界からq(v,v)+∥v∥22q(v,v)+\|v\|_2^2は∥v∥H12\|v\|_{H^1}^2と同値である。§E15.16 定理 3.1によりqqは非負自己共役実現をもち、§E15.24 命題 1.1により平方根の主張を得る。弱形式§E15.16 定義 1.3から、V=H01V=H_0^1は斉次 Dirichlet 条件に、V=H1V=H^1は零 conormal 流束に対応する。全空間の場合には§E15.24 命題 1.4を用いる。▨

定理 4.2.Ω,a,A,V,θ,c\Omega,a,A,V,\theta,cを命題 4.1の有界領域の設定とするか、Ω=Rn\Omega=\R^n、a=Ia=I、A=−ΔA=-\Delta、V=H1(Rn)V=H^1(\R^n)、c=1c=1とする。T>0T>0、(u0,u1)∈V×L2(Ω)(u_0,u_1)\in V\times L^2(\Omega)、f∈L1(0,T;L2(Ω))f\in L^1(0,T;L^2(\Omega))に対するエネルギー解uuについて

e(t,x)=12(∣ut∣2+(a∇u)⋅∇u‾),J(t,x)=ℜ(ut‾ a∇u)e(t,x)=\frac12\bigl(|u_t|^2+(a\nabla u)\cdot\overline{\nabla u}\bigr), \qquad J(t,x)=\Re(\overline{u_t}\,a\nabla u)

とおく。このとき、滑らかな解には

∂te=div⁡J+ℜ(fut‾)\partial_te=\operatorname{div}J+\Re(f\overline{u_t})

が点ごとに成り立つ。任意の非負な有界 Lipschitz 関数χ:Ω‾×[0,T]→R\chi:\overline\Omega\times[0,T]\to\Rに対し、エネルギー解は任意の0≤r≤t≤T0\le r\le t\le Tについて

∫Ωχ(t)e(t) dx−∫Ωχ(r)e(r) dx=∫rt∫Ω[χse−∇χ⋅J+χℜ(fus‾)]dx ds\begin{aligned} \int_\Omega\chi(t)e(t)\,dx-\int_\Omega\chi(r)e(r)\,dx ={}&\int_r^t\int_\Omega\bigl[\chi_s e-\nabla\chi\cdot J +\chi\Re(f\overline{u_s})\bigr]dx\,ds \end{aligned}

を満たす。

証明. 滑らかな解に対し、aaが実対称かつ時間に依存しないことから

∂te=ℜ(uttut‾)+ℜ((a∇u)⋅∇ut‾)=ℜ((div⁡(a∇u)+f)ut‾)+ℜ((a∇u)⋅∇ut‾)=div⁡J+ℜ(fut‾)\begin{aligned} \partial_te &=\Re(u_{tt}\overline{u_t}) +\Re((a\nabla u)\cdot\nabla\overline{u_t})\\ &=\Re((\operatorname{div}(a\nabla u)+f)\overline{u_t}) +\Re((a\nabla u)\cdot\nabla\overline{u_t})\\ &=\operatorname{div}J+\Re(f\overline{u_t}) \end{aligned}

である。滑らかなχ\chiを掛け、§E15.12 定理 1.3と時間積分を適用する。境界流束はχℜ(ut‾(a∇u)⋅ν)\chi\Re(\overline{u_t}(a\nabla u)\cdot\nu)である。Dirichlet 条件ではutu_tのトレースが零、conormal Neumann 条件では(a∇u)⋅ν=0(a\nabla u)\cdot\nu=0であるから、どちらも境界流束は零である。

χ\chiが Lipschitz ならば、閉領域の近傍へ Lipschitz 拡張して平滑化することにより、χk→χ\chi_k\to\chiを一様に、∂sχk→∂sχ\partial_s\chi_k\to\partial_s\chiと∇χk→∇χ\nabla\chi_k\to\nabla\chiをほとんど至る所で、かつ一階微分を一様有界に保つ近似を得る。弱微分の存在と有界性は§E15.8 命題 2.3による。積の微分は§E15.8 命題 3.1による。平滑化した恒等式へ優収束を適用すると、Lipschitz なχ\chiについても恒等式を得る。

PN=1[0,N](A)P_N=1_{[0,N]}(A)とする。初期値にPNP_Nを適用し、外力を時間方向に平滑化してからPNP_Nを適用することにより、初期値はV×L2V\times L^2、外力はLt1L2L^1_tL^2で収束し、近似データがすべてのD(Am)D(A^m)に属する列をとることができる。対応する解もすべてのD(Am)D(A^m)に時間について滑らかな値をとる。有界領域では§E15.23 定理 5.1を反復適用すると、これらの解と時間微分は任意の整数次数の Sobolev 空間に属する。十分大きい整数kkに対し、§E15.10 系 4.2のp=2p=2の場合を適用すると、上の古典的な積分計算に必要なC2C^2級の空間正則性を得る。境界条件はD(A)D(A)の弱形式から従う。全空間では Fourier 変換をv^(ξ)=∫e−2πix⋅ξv(x) dx\widehat v(\xi)=\int e^{-2\pi ix\cdot\xi}v(x)\,dxとし、PNP_Nは4π2∣ξ∣2≤N4\pi^2|\xi|^2\le Nの周波数切断である。§E15.24 命題 1.4により任意の整数次数の Sobolev 正則性を得る。有界球に制限して§E15.10 系 4.2を適用すると、全空間でも局所的なC2C^2級の空間正則性を得る。全空間の積分計算には半径LLの空間切断を加える。その勾配をC/LC/Lで抑えると、余分な項はCL−1∫rt∥us∥2∥∇u∥2 dsC L^{-1}\int_r^t\|u_s\|_2\|\nabla u\|_2\,ds以下であり、L→∞L\to\inftyで零になる。

定理 2.2の連続依存性により、近似解はCtV∩Ct1L2C_tV\cap C_t^1L^2でuuに収束する。したがって近似解のeeとJJはCtL1C_tL^1で収束する。外力項も、外力のLt1L2L^1_tL^2収束と速度のCtL2C_tL^2収束によりL1L^1で収束する。重みとその微分は有界であるから、近似解の恒等式のすべての項で極限を取ることができる。▨

5 有限伝播速度と依存領域

定理 5.1.Ω,a,A,V,c\Omega,a,A,V,cを定理 4.2のいずれかの設定とし、uuを初期値(u0,u1)∈V×L2(Ω)(u_0,u_1)\in V\times L^2(\Omega)、外力f∈L1(0,T;L2(Ω))f\in L^1(0,T;L^2(\Omega))のエネルギー解とする。eeを定理 4.2のエネルギー密度とする。x0∈Rnx_0\in\R^n、R>0R>0に対し

Ds=Ω∩B(x0,R−cs)(0≤s<min⁡{T,R/c})D_s=\Omega\cap B(x_0,R-cs)\qquad(0\le s<\min\{T,R/c\})

とおく。このとき任意の0≤t<min⁡{T,R/c}0\le t<\min\{T,R/c\}について

2∫Dte(t,x) dx≤2∫D0e(0,x) dx+∫0t∥f(s)∥L2(Ds) ds\sqrt{2\int_{D_t}e(t,x)\,dx} \le\sqrt{2\int_{D_0}e(0,x)\,dx} +\int_0^t\|f(s)\|_{L^2(D_s)}\,ds

が成り立つ。さらに、u0=u1=0u_0=u_1=0がD0D_0上でほとんど至る所成り立ち、f=0f=0が{(s,x):0<s<t, x∈Ds}\{(s,x):0<s<t,\ x\in D_s\}上でほとんど至る所成り立つならば、u(s)=0u(s)=0が各0≤s≤t0\le s\le tのDsD_s上でほとんど至る所成り立つ。したがって二つの解がこの円錐内で一致するかどうかは、底面の初期データと円錐内の外力だけで決まる。

証明.z=a1/2∇uz=a^{1/2}\nabla uとおくと、a≤c2Ia\le c^2Iより

∣J∣≤∣ut∣ ∣a∇u∣≤c∣ut∣∣z∣≤c2(∣ut∣2+∣z∣2)=ce|J|\le|u_t|\,|a\nabla u| \le c|u_t||z| \le\frac c2(|u_t|^2+|z|^2)=ce

である。h:R→[0,1]h:\R\to[0,1]をh(q)=0h(q)=0(q≤0q\le0)、h(q)=qh(q)=q(0<q<10<q<1)、h(q)=1h(q)=1(q≥1q\ge1)で定める。0<ε<R−ct0<\varepsilon<R-ctに対して

χε(s,x)=h(R−cs−∣x−x0∣ε),Fε(s)=∫Ωχε(s,x)e(s,x) dx\chi_\varepsilon(s,x)=h\left(\frac{R-cs-|x-x_0|}{\varepsilon}\right),\qquad F_\varepsilon(s)=\int_\Omega\chi_\varepsilon(s,x)e(s,x)\,dx

とおく。距離関数の Lipschitz 定数は11であるから、§E15.8 命題 2.3と連鎖律により、ほとんど至る所で

∂sχε=−cεh′,∣∇χε∣≤1εh′\partial_s\chi_\varepsilon=-\frac c\varepsilon h',\qquad |\nabla\chi_\varepsilon|\le\frac1\varepsilon h'

である。h′≥0h'\ge0と∣J∣≤ce|J|\le ceから

(∂sχε)e−∇χε⋅J≤0(\partial_s\chi_\varepsilon)e-\nabla\chi_\varepsilon\cdot J\le0

を得る。定理 4.2によりFεF_\varepsilonは絶対連続であり、ほとんどすべてのs∈[0,t]s\in[0,t]について

Fε′(s)≤(∫Ωχε∣f∣2 dx)1/22Fε(s)F_\varepsilon'(s) \le\left(\int_\Omega\chi_\varepsilon|f|^2\,dx\right)^{1/2} \sqrt{2F_\varepsilon(s)}

となる。ここで重み付き Cauchy–Schwarz 不等式を用いた。δ>0\delta>0に対して2Fε+δ\sqrt{2F_\varepsilon+\delta}を微分し、積分してからδ↓0\delta\downarrow0とすると

2Fε(t)≤2Fε(0)+∫0t(∫Ωχε(s,x)∣f(s,x)∣2 dx)1/2ds\sqrt{2F_\varepsilon(t)}\le\sqrt{2F_\varepsilon(0)} +\int_0^t\left(\int_\Omega\chi_\varepsilon(s,x)|f(s,x)|^2\,dx\right)^{1/2}ds

を得る。ε↓0\varepsilon\downarrow0のとき、χε(s,x)\chi_\varepsilon(s,x)はDsD_sの特性関数へ点ごとに収束する。空間積分には優収束を、時間積分には支配関数∥f(s)∥2\|f(s)\|_2を用いて優収束を適用すると、主張の評価が従う。切断は領域全体で用いたので、球面と∂Ω\partial\Omegaの接触についての仮定は必要ない。

零データの仮定の下では、u0∈H1u_0\in H^1がD0D_0上で零であることから∇u0=0\nabla u_0=0もその開集合上でほとんど至る所成り立つ。局所評価によりus(s)=0u_s(s)=0が各s≤ts\le tのDsD_s上で成立する。0≤τ≤t0\le\tau\le tを固定すると、Dτ⊂DsD_\tau\subset D_s(s≤τs\le\tau)であるから、L2(Dτ)L^2(D_\tau)において

u(τ)∣Dτ=u0∣Dτ+∫0τus(s)∣Dτ ds=0u(\tau)|_{D_\tau}=u_0|_{D_\tau}+\int_0^\tau u_s(s)|_{D_\tau}\,ds=0

となる。二つの解の差へこの結論を適用すると、依存領域の主張を得る。▨

例 5.2. 全空間Rn\R^nでa=diag⁡(c12,…,cn2)a=\operatorname{diag}(c_1^2,\ldots,c_n^2)、cj>0c_j>0とする。T>0T>0、u0∈H1(Rn)u_0\in H^1(\R^n)、u1∈L2(Rn)u_1\in L^2(\R^n)、f∈L1(0,T;L2(Rn))f\in L^1(0,T;L^2(\R^n))とし、u∈C([0,T];H1(Rn))∩C1([0,T];L2(Rn))u\in C([0,T];H^1(\R^n))\cap C^1([0,T];L^2(\R^n))を初期条件u(0)=u0u(0)=u_0、ut(0)=u1u_t(0)=u_1と方程式

utt−∑j=1ncj2∂xj2u=fu_{tt}-\sum_{j=1}^nc_j^2\partial_{x_j}^2u=f

を超関数の意味で満たすエネルギー解とする。xj=cjyjx_j=c_jy_j、v(t,y)=u(t,c1y1,…,cnyn)v(t,y)=u(t,c_1y_1,\ldots,c_ny_n)とおく。このとき、超関数の意味で

∂yj2v=cj2∂xj2u,vtt−Δyv=f(t,c1y1,…,cnyn)\partial_{y_j}^2v=c_j^2\partial_{x_j}^2u,\qquad v_{tt}-\Delta_yv=f(t,c_1y_1,\ldots,c_ny_n)

である。変数変換の Jacobian は正の定数∏jcj\prod_jc_jであり、変位のH1H^1ノルム、速度のL2L^2ノルム、外力のLt1L2L^1_tL^2ノルムは変換前後で同値である。したがって定理 5.1の全空間の場合をvvに適用することができる。点(t,x)(t,x)の後方円錐の時刻ssの断面は、xx座標では

∑j=1n(zj−xj)2cj2<(t−s)2\sum_{j=1}^n\frac{(z_j-x_j)^2}{c_j^2}<(t-s)^2

という楕円体になる。この楕円体は半径(max⁡jcj)(t−s)(\max_jc_j)(t-s)の Euclid 球に含まれる。行列上界から得た球による速度評価は全方向の共通上界であり、座標変換による評価は各方向の係数を保持する。一次元の斉次方程式ではu(t,x)=g(x−ct)u(t,x)=g(x-ct)が任意のg∈Cc∞(R)g\in C_c^\infty(\R)に対して解であるから、進行方向の速度ccは実際に生じる。

前提記事

9 本の記事・単元を表示