双曲型偏微分方程式

同義語:双曲型方程式hyperbolic partial differential equationhyperbolic equation

概要

双曲型偏微分方程式(hyperbolic partial differential equation)とは、波動方程式を代表とし、時間について 2 階、空間について一様楕円型の 2 階の作用素をもつ $u_{tt}-\sum\partial_i(a^{ij}\partial_ju)+\cdots=f$ の形の方程式と、その 1 階の系である。係数が変わってもエネルギーの増え方は指数関数で抑えられ、初期境界値問題の解は一意でデータに連続に依存し、初期値の影響は $(a^{ij})$ の固有値の上限 $\Lambda$ について速さ $\sqrt\Lambda$ を超えて伝わらない。定数係数の系 $u_t=\sum A_j\partial_ju$ では $\sum\xi_jA_j$ の固有値が実数であることが適切性に必要で、対称双曲型または狭義双曲型なら解は $L^2$ ノルムで抑えられる。

$$\newcommand{C}[0]{\mathbb{C}} \newcommand{div}[0]{\mathbin{÷}} \newcommand{N}[0]{\mathbb{N}} \newcommand{Q}[0]{\mathbb{Q}} \newcommand{R}[0]{\mathbb{R}} \newcommand{Z}[0]{\mathbb{Z}} $$

前提知識: 偏微分方程式, 波動方程式, Gaussの発散定理, Fourier変換

双曲型偏微分方程式(hyperbolic partial differential equation)は、波のように有限の速さで伝わる現象を記述する発展方程式の型である。代表は 波動方程式 $u_{tt}=\Delta u$ で、2 階の線形方程式では、時間について 2 階、空間の変数について楕円型の 2 階の作用素をもつ
$$ u_{tt}-\sum_{i,j}\partial_i\bigl(a^{ij}(x,t)\,\partial_ju\bigr)+(\text{低階の項})=f $$
の形の方程式をいう。放物型偏微分方程式 では初期値の影響が瞬時に遠くまで及び、解は正の時刻で滑らかになるのに対し、双曲型の方程式では影響は有限の速さでしか伝わらず、初期値の特異性はそのまま運ばれる。
この記事では、まず 2 階の線形双曲型方程式を定義し、偏微分方程式 の定義「楕円型・放物型・双曲型」との関係を述べる。次に、係数が場所と時刻によって変わる場合に、エネルギー評価(解の一意性とデータへの連続な依存)と有限伝播速度を証明する。後半では 1 階の連立方程式(系)について双曲性・対称双曲性・狭義双曲性を定義し、定数係数の場合に Fourier 変換で初期値問題を解いて、双曲性が適切性に必要であることと、固有値が実数というだけでは足りないことを示す。

定義

以下、空間の変数を $x\in\mathbb{R}^n$、時間を $t$ とし、$\partial_i=\partial/\partial x_i$、$\nabla u=(\partial_1u,\dots,\partial_nu)$ と書く。添字の和は $1$ から $n$ までとる。

2 階の線形双曲型偏微分方程式

$\Omega\subset\mathbb{R}^n$ を開集合、$T>0$ とし、$\overline\Omega\times[0,T]$ の上の関数 $a^{ij}=a^{ji}$($C^1$ 級)と $b^1,\dots,b^n,b^0,q$(連続で有界)を与える。作用素
$$ Pu:=u_{tt}-\sum_{i,j}\partial_i\bigl(a^{ij}\partial_ju\bigr)+\sum_ib^i\partial_iu+b^0u_t+qu $$
が、定数 $0< \theta\le\Lambda$ について
$$ \theta\lvert\xi\rvert^2\le\sum_{i,j}a^{ij}(x,t)\xi_i\xi_j\le\Lambda\lvert\xi\rvert^2\qquad(\xi\in\mathbb{R}^n,\ (x,t)\in\overline\Omega\times[0,T]) $$
を満たすとき、方程式 $Pu=f$ を (一様)双曲型 の 2 階線形方程式という。$a^{ij}$ が $t$ によらず $a^{ij}=\delta_{ij}$ のとき、$P$ の主要部は波動作用素 $\partial_t^2-\Delta$ である。

$a^{ij}$ が $C^1$ 級なら $\sum_{i,j}a^{ij}\partial_i\partial_ju$ と $\sum_{i,j}\partial_i(a^{ij}\partial_ju)$ の差は 1 階の項 $\sum_{i,j}(\partial_ia^{ij})\partial_ju$ なので、発散形でない方程式 $u_{tt}=\sum a^{ij}\partial_i\partial_ju+\cdots$ も同じ定義に入る。変数 $(x,t)\in\mathbb{R}^{n+1}$ についての主要部の係数行列は、ブロックに分けて $\begin{pmatrix}-(a^{ij})&0\\0&1\end{pmatrix}$ であり、負の固有値が $n$ 個、正の固有値が $1$ 個ある。よってこの方程式は 偏微分方程式 の定義「楕円型・放物型・双曲型」の意味で各点で双曲型である。逆に、定数係数の 2 階の作用素が双曲型なら、主軸定理による直交変換で正の固有値(あるいは全体の符号を変えて負の固有値)の固有ベクトルの方向を時間の軸にとり、座標を拡大縮小すれば、主要部は $\partial_t^2-\Delta$ になる。定義では時間の方向をあらかじめ決め、その方向で初期値を与える問題を考える。

双曲型の方程式の例
  1. 媒質が一様でない波動方程式:密度 $\rho(x)>0$、弾性係数 $k(x)>0$($C^1$ 級で上下から正の定数で抑えられる)について $\rho u_{tt}=\nabla\cdot(k\nabla u)$。$\rho$ で割ると $u_{tt}=\sum_i\partial_i\bigl(\tfrac k\rho\partial_iu\bigr)-k\nabla\bigl(\tfrac1\rho\bigr)\cdot\nabla u$ となり、$a^{ij}=\frac k\rho\delta_{ij}$ である。点 $x$ での波の速さは $\sqrt{k(x)/\rho(x)}$ である。
  2. 電信方程式:$u_{tt}+2\gamma u_t=u_{xx}$($\gamma>0$)。$b^0=2\gamma$ の減衰の項をもつ。
  3. Klein–Gordon 方程式:$u_{tt}-\Delta u+m^2u=0$。$q=m^2$ である。

特性面と伝播の速さ:関数 $\varphi(x,t)$ の等位面 $\{\varphi=0\}$ で、各点で
$$ \varphi_t^2=\sum_{i,j}a^{ij}\varphi_{x_i}\varphi_{x_j} $$
となるものを $P$ の 特性面 という。主要部の係数行列で $\nabla_{(x,t)}\varphi$ を測った 2 次形式が $0$ になる面で、初期値を与えても 2 階の微分が方程式から決まらない面である。$a^{ij}=c^2\delta_{ij}$($c>0$ は定数)なら、円錐 $\lvert x-x_0\rvert=c\lvert t-t_0\rvert$ が特性面であり、波は速さ $c$ で伝わる。一般の係数でも、係数を 1 点の値で固定した定数係数の作用素について、方向 $\omega$($\lvert\omega\rvert=1$)に進む平面の特性面 $x\cdot\omega=\lambda t$ の速さは $\lambda=\sqrt{\sum a^{ij}\omega_i\omega_j}\le\sqrt\Lambda$ なので、$\sqrt\Lambda$ が伝播の速さの上限になる。これを下の thm-hpde-finite-speed で正確に述べる。

エネルギー評価

波動方程式では、運動エネルギーと弾性エネルギーの和 $\frac12\int(u_t^2+\lvert\nabla u\rvert^2)$ が保存される(偏微分方程式 の定理「波動方程式のエネルギー保存と一意性」)。係数が変わる場合や低階の項がある場合には保存はしないが、増え方が指数関数で抑えられる。以下 $\|\cdot\|$ は $L^2(\Omega)$ のノルムとし、定数
$$ \beta:=\sup\Bigl(\Bigl(\sum_i(b^i)^2\Bigr)^{1/2}+\lvert b^0\rvert+\lvert q\rvert\Bigr),\qquad\kappa:=\sup\bigl\|(\partial_ta^{ij}(x,t))_{i,j}\bigr\| $$
(上限は $\overline\Omega\times[0,T]$ でとり、$\kappa$ の中は行列の作用素ノルム)を使う。

エネルギー評価

$\overline\Omega$ を滑らかな境界をもつコンパクトな領域(Gaussの発散定理 の定義「滑らかな境界をもつコンパクトな領域」の意味)とし、$P$ を def-hpde-main の双曲型の作用素とする。$\overline\Omega$ を含む開集合 $U$ があって、$a^{ij}$ は $U\times[0,T]$ の上の $C^1$ 級の関数に延びているとする。$u$ は $U\times[0,T]$ の上で $C^2$ 級、$\partial\Omega\times[0,T]$ の上で $u=0$ とし、$f:=Pu$ とおく。
$$ E(t):=\frac12\int_\Omega\Bigl(u_t^2+\sum_{i,j}a^{ij}\partial_iu\,\partial_ju+u^2\Bigr)dx $$
とおくと、$C:=2+4\beta+(\beta+\kappa)/\theta$ について、$0\le t\le T$ で
$$ E(t)\le e^{Ct}\Bigl(E(0)+\frac12\int_0^t\| f(\cdot,s)\|^2ds\Bigr) $$
が成り立つ。

$E$ を微分して部分積分する

要点:$E'(t)$ を計算すると、主要部の寄与は部分積分で $\int u_t\bigl(u_{tt}-\sum\partial_i(a^{ij}\partial_ju)\bigr)$ にまとまり、境界の項は $u_t=0$ で消える。残りは $E$ と $\| f\|^2$ で抑えられ、Grönwall 型の議論で結論を得る。

詳しい証明を開く

被積分関数は $\overline\Omega\times[0,T]$ の上で $C^1$ 級なので、積分記号下の微分 により $E$ は $C^1$ 級で、$a^{ij}=a^{ji}$ を使うと

$$E'(t)=\int_\Omega\Bigl(u_tu_{tt}+\sum_{i,j}a^{ij}\partial_ju\,\partial_iu_t+\frac12\sum_{i,j}(\partial_ta^{ij})\partial_iu\,\partial_ju+uu_t\Bigr)dx$$

である。ベクトル場 $w_i:=u_t\sum_ja^{ij}\partial_ju$ は各時刻で $U$ の上の $C^1$ 級のベクトル場で、$\sum_i\partial_iw_i=\sum_{i,j}a^{ij}\partial_ju\,\partial_iu_t+u_t\sum_{i,j}\partial_i(a^{ij}\partial_ju)$ である。$\partial\Omega$ の上では $u=0$ がすべての時刻で成り立つので $u_t=0$、よって $w=0$ であり、Gaussの発散定理 から $\int_\Omega\sum_i\partial_iw_i\,dx=0$ である。したがって

$$E'(t)=\int_\Omega\Bigl(u_t\bigl(f-\textstyle\sum_ib^i\partial_iu-b^0u_t-qu\bigr)+\frac12\sum_{i,j}(\partial_ta^{ij})\partial_iu\,\partial_ju+uu_t\Bigr)dx$$

となる。$\lvert\nabla u\rvert^2\le\theta^{-1}\sum a^{ij}\partial_iu\,\partial_ju$ と $2\lvert yz\rvert\le y^2+z^2$ を使って各項を評価する:$\lvert u_tf\rvert\le\frac12(u_t^2+f^2)$、$\bigl\lvert u_t\sum b^i\partial_iu\bigr\rvert\le\frac\beta2(u_t^2+\lvert\nabla u\rvert^2)$、$\lvert b^0\rvert u_t^2\le\beta u_t^2$、$\lvert(1-q)uu_t\rvert\le\frac{\beta+1}2(u^2+u_t^2)$、$\frac12\bigl\lvert\sum(\partial_ta^{ij})\partial_iu\,\partial_ju\bigr\rvert\le\frac\kappa2\lvert\nabla u\rvert^2$。$u_t^2$ の係数の和は $1+2\beta$、$\sum a^{ij}\partial_iu\,\partial_ju$ の係数は $\frac{\beta+\kappa}{2\theta}$ 以下、$u^2$ の係数は $\frac{\beta+1}2$ で、いずれも $\frac C2$ 以下なので

$$E'(t)\le CE(t)+\frac12\| f(\cdot,t)\|^2$$

である。$F(t):=e^{-Ct}E(t)$ とおくと $F'(t)=e^{-Ct}(E'-CE)\le\frac12\| f(\cdot,t)\|^2$ なので、$0$ から $t$ まで積分して $F(t)\le E(0)+\frac12\int_0^t\| f\|^2ds$、すなわち結論を得る。$\square$

初期境界値問題の一意性とデータへの連続な依存

thm-hpde-energy の設定で、$u_1,u_2$ が同じ境界条件 $u=0$($\partial\Omega\times[0,T]$)を満たし、$Pu_k=f_k$ とする。$w:=u_1-u_2$ のエネルギーを $E_w$ とすると
$$ \sup_{0\le t\le T}E_w(t)\le e^{CT}\Bigl(E_w(0)+\frac12\int_0^T\| f_1-f_2\|^2ds\Bigr) $$
である。とくに $f_1=f_2$、$u_1(\cdot,0)=u_2(\cdot,0)$、$\partial_tu_1(\cdot,0)=\partial_tu_2(\cdot,0)$ なら $u_1=u_2$ である。

$P$ は線形なので $Pw=f_1-f_2$ であり、$w$ も境界で $0$ である。thm-hpde-energy を $w$ に使えばよい。後半では右辺が $0$ になるので $E_w=0$、被積分関数の $w^2$ の項から $w=0$ である($w$ は連続)。$\square$

$E(0)$ は初期値 $u(\cdot,0)$ の $1$ 階微分までと $u_t(\cdot,0)$ だけで決まる。評価は、初期値と右辺が $L^2$ の意味で少し変われば解も各時刻で $L^2$ の意味で($u$ と $\nabla u$ と $u_t$ について)少ししか変わらないことを表す。初期値問題がこのように解の存在・一意性・データへの連続な依存をもつことを 適切(well-posed)という。双曲型の方程式では、時間の向きを逆にした $v(x,t):=u(x,T-t)$ も同じ型の方程式($b^0$ の符号が変わる)を満たすので、同じ評価が過去に向かっても成り立つ。これは、時間を逆にたどると不適切になる放物型の方程式(放物型偏微分方程式 の反例の表の 1 行目)と対照的である。

有限伝播速度

全空間の解について、初期値がある球の上で $0$ なら、その球の上に立つ円錐の中で解は $0$ である。円錐の側面の傾きは $\sqrt\Lambda$ でよい。

有限伝播速度

$\Omega=\mathbb{R}^n$ とし、$P$ を def-hpde-main の双曲型の作用素、$c:=\sqrt\Lambda$ とする。$u$ は $\mathbb{R}^n\times[0,T]$ の上で $C^2$ 級とし、$x_0\in\mathbb{R}^n$、$R>0$ について
$$ K:=\{(x,t)\mid0\le t\le T,\ t\le R/c,\ \lvert x-x_0\rvert\le R-ct\} $$
とおく。$u(\cdot,0)$ と $u_t(\cdot,0)$ が閉球 $\overline B(x_0,R)$ の上で $0$ で、$K$ の上で $Pu=0$ なら、$K$ の上で $u=0$ である。

縮む球の上のエネルギー

要点:時刻 $t$ の球 $B(t):=B(x_0,R-ct)$ の上のエネルギー $E(t)$ を微分すると、球が縮む分のエネルギーの流出($-\frac c2\int_{\partial B}e$)が、部分積分で現れる境界の項を打ち消す。打ち消すことの鍵は、$a$ を内積とみた Cauchy–Schwarz の不等式 $\lvert\sum a^{ij}\partial_ju\,\nu_i\rvert\le c\,(\sum a^{ij}\partial_iu\,\partial_ju)^{1/2}$ である。残りは $E'\le CE$ となり、$E(0)=0$ から $E=0$ を得る。

詳しい証明を開く

$t_1:=\min(T,R/c)$ とし、$0\le t< t_1$ について $\rho(t):=R-ct$、$e:=u_t^2+\sum a^{ij}\partial_iu\,\partial_ju+u^2$、$E(t):=\frac12\int_{B(t)}e(x,t)\,dx$ とおく。

段 1(縮む球の上の積分の微分).$e$ は $C^1$ 級である。調和関数 の注意「この記事で使う積分の基本事項」の 1・2 により $G(r,t):=\int_{B(x_0,r)}e(x,t)\,dx=nV(B)\int_0^rs^{n-1}\int_Se(x_0+s\zeta,t)\,d\sigma(\zeta)\,ds$ であり($S$ は単位球面、$\sigma$ は正規化した表面積測度、$V(B)$ は単位球の体積)、$\partial_rG=\int_{\partial B(x_0,r)}e\,dS$、$\partial_tG=\int_{B(x_0,r)}e_t\,dx$(積分記号下の微分)はどちらも連続である。よって $E(t)=\frac12G(\rho(t),t)$ は微分可能で

$$E'(t)=\frac12\int_{B(t)}e_t\,dx-\frac c2\int_{\partial B(t)}e\,dS$$

である。

段 2(部分積分).$\frac12e_t=u_tu_{tt}+\sum a^{ij}\partial_ju\,\partial_iu_t+\frac12\sum(\partial_ta^{ij})\partial_iu\,\partial_ju+uu_t$ である。prf-hpde-energy と同じベクトル場 $w$ に 調和関数 の注意の 3(球の上の発散定理)を使うと、外向き単位法線を $\nu$ として

$$\int_{B(t)}\sum a^{ij}\partial_ju\,\partial_iu_t\,dx=-\int_{B(t)}u_t\sum\partial_i(a^{ij}\partial_ju)\,dx+\int_{\partial B(t)}u_t\sum a^{ij}\partial_ju\,\nu_i\,dS$$

である。$K$ の上で $Pu=0$ なので、prf-hpde-energy と同じ評価($f=0$)により内部の積分の和は $CE(t)$ 以下である。ここで $C$ の中の $\beta$、$\kappa$ は $K$ の上での上限とする($K$ はコンパクトなので有限である)。

段 3(境界の項).各点で $\langle\xi,\eta\rangle_a:=\sum a^{ij}\xi_i\eta_j$ は内積なので、Cauchy–Schwarzの不等式 と $\langle\nu,\nu\rangle_a\le\Lambda\lvert\nu\rvert^2=c^2$ から

$$\Bigl\lvert u_t\sum a^{ij}\partial_ju\,\nu_i\Bigr\rvert\le\lvert u_t\rvert\,\langle\nabla u,\nabla u\rangle_a^{1/2}\,c\le\frac c2\bigl(u_t^2+\langle\nabla u,\nabla u\rangle_a\bigr)\le\frac c2e$$

である。よって境界の項の和 $\int_{\partial B(t)}u_t\sum a^{ij}\partial_ju\,\nu_i\,dS-\frac c2\int_{\partial B(t)}e\,dS$ は $0$ 以下である。

段 4.以上から $[0,t_1)$ で $E'(t)\le CE(t)$ であり、$(e^{-Ct}E)'\le0$ と $E(0)=0$(初期値が $\overline B(x_0,R)$ で $0$ なので、そこで $u=u_t=0$、$\nabla u=0$)から $E=0$ である。被積分関数の $u^2$ の項から、$0\le t< t_1$ の各時刻で開球 $B(t)$ の上で $u=0$ である。$u$ は連続なので、$K$ の残りの点(各時刻の閉球の境界と、時刻 $t_1$ の切り口。$t_1=R/c$ なら頂点、$t_1=T< R/c$ なら閉球 $\lvert x-x_0\rvert\le R-cT$)でも $u=0$ である。$\square$

解の値を決める初期値の範囲

thm-hpde-finite-speed の設定で、2 つの $C^2$ 級の解 $u_1,u_2$ が $K$ の上で $Pu_1=Pu_2$ を満たし、初期値が $\overline B(x_0,R)$ の上で一致するなら、$K$ の上で $u_1=u_2$ である。つまり $K$ の点での解の値は、$\overline B(x_0,R)$ の上の初期値と $K$ の上の右辺だけで決まる。

差 $u_1-u_2$ に thm-hpde-finite-speed を使えばよい。$\square$

$a^{ij}=\delta_{ij}$ なら $c=1$ で、依存領域 の定理「波動方程式の依存領域(エネルギー評価)」と同じ主張である。証明はその定理と同じく縮む球の上のエネルギーを使うが、係数が変わる場合には、境界の項を $a$ の内積の Cauchy–Schwarz の不等式で抑えるところが要になる。点 $(x,t)$ での解の値に影響しうるのは、初期時刻で $\lvert y-x\rvert\le ct$ の範囲だけであり、影響は速さ $c=\sqrt\Lambda$ を超えて伝わらない。この性質から、全空間の解についても、初期値の台と右辺の台がコンパクトなら各時刻の解の台はコンパクトである(台から遠い点のまわりの球に thm-hpde-finite-speed を使う)。

弱解の存在

エネルギー評価は一意性と安定性を与えるが、解の存在は別に示す必要がある。初期境界値問題については、Hilbert 空間の方法で弱解の存在が示されている。

J.-L. Lions の定理の応用(双曲型の初期境界値問題)

$\Omega\subset\mathbb{R}^n$ を開集合、$a^{ij}=a^{ji}$、$\partial_ta^{ij}$、$a_0$、$\partial_ta_0$ を $\Omega\times(0,T)$ の上の有界な可測関数とし、ある $\theta>0$ について $\sum a^{ij}\xi_i\xi_j\ge\theta\lvert\xi\rvert^2$ とする。$f\in L^2(\Omega\times(0,T))$、$u_0\in H^1_0(\Omega)$、$v_0\in L^2(\Omega)$ なら、
$$ u_{tt}-\sum_{i,j}\partial_i\bigl(a^{ij}\partial_ju\bigr)+a_0u=f,\qquad u=0\ (\partial\Omega\times(0,T)),\qquad u(\cdot,0)=u_0,\quad u_t(\cdot,0)=v_0 $$
の弱解 $u\in C([0,T];H^1_0(\Omega))$、$u_t\in C([0,T];L^2(\Omega))$、$u_{tt}\in L^2(0,T;H^{-1}(\Omega))$ がただ 1 つ存在する。

証明の所在

Bre11 第 10 章の注 6 の Theorem 10.14(J.-L. Lions の定理)とその応用(pp. 345–346)に述べられている。この記事では証明しない。

$a^{ij}=\delta_{ij}$、$a_0=0$、$f=0$ で、$\Omega$ が $C^\infty$ 級で境界が有界のときは、Bre11 Theorem 10.7(p. 336)が Hille–吉田の定理 を使って、$u_0\in H^2\cap H^1_0$、$v_0\in H^1_0$ に対して強い意味の解を与えている。$H^1_0(\Omega)$ は境界で $0$ となる Sobolev 空間である(Sobolev空間)。初期値や係数がより滑らかで、境界での両立条件を満たせば解も滑らかになる(Bre11 Theorem 10.8、p. 336 と第 10 章の注 6、p. 346。この記事では証明しない)。双曲型の方程式は熱方程式のような平滑化をもたないので(Bre11 §10.3 の Remark 8、p. 339)、解の滑らかさは初期値の滑らかさで決まる。

1 階の系の双曲性

方程式を 1 階の連立方程式に直すと、双曲性は行列の固有値の条件として述べられる。以下、$A_1,\dots,A_n$ を $N$ 次の実正方行列、$B$ を $N$ 次の複素正方行列とし、$\mathbb{C}^N$ に値をとる未知関数 $u(x,t)$ についての定数係数の系
$$ u_t=\sum_jA_j\partial_ju+Bu\qquad(*) $$
を考える。$\xi\in\mathbb{R}^n$ について $A(\xi):=\sum_j\xi_jA_j$ とおく。

双曲型・対称双曲型・狭義双曲型の系
  1. すべての $\xi\in\mathbb{R}^n$ について $A(\xi)$ の固有値がすべて実数のとき、系 $(*)$ は 双曲型 であるという。
  2. $A_1,\dots,A_n$ がすべて対称行列のとき 対称双曲型 であるという。
  3. すべての $\xi\ne0$ について $A(\xi)$ が相異なる $N$ 個の実数の固有値をもつとき 狭義双曲型 であるという。

実対称行列の固有値は実数なので(Hermite行列)、対称双曲型の系も狭義双曲型の系も双曲型である。双曲型の判定には主要部 $\sum A_j\partial_j$ だけを使い、低階の項 $Bu$ は関係しない。係数が $(x,t)$ によって変わる系 $A_0(x,t)u_t=\sum_jA_j(x,t)\partial_ju+B(x,t)u$ についても、$A_0$ が正定値の対称行列で $A_j$ が対称なものを対称双曲型という。

波動方程式を対称双曲型の系に直す

$u$ が $u_{tt}=\Delta u$ の $C^2$ 級の解なら、$U=(U_0,U_1,\dots,U_n):=(u_t,\partial_1u,\dots,\partial_nu)$ は $C^1$ 級で
$$ \partial_tU_0=\sum_j\partial_jU_j,\qquad\partial_tU_j=\partial_jU_0\quad(j=1,\dots,n) $$
を満たす($N=n+1$)。$A_j$ は $(0,j)$ 成分と $(j,0)$ 成分だけが $1$ の対称行列なので、この系は対称双曲型である。$A(\xi)=\begin{pmatrix}0&\xi^{\mathsf T}\\\xi&0\end{pmatrix}$ の固有値は $\pm\lvert\xi\rvert$(固有ベクトル $(\pm\lvert\xi\rvert,\xi)$)と $0$(固有ベクトル $(0,\eta)$、$\eta\perp\xi$。重複度 $n-1$)なので、$n=1$ では狭義双曲型で、$n\ge2$ では狭義双曲型でない。$\lvert U\rvert^2=u_t^2+\lvert\nabla u\rvert^2$ はエネルギーの密度である。
同じようにして Laplace 方程式 $u_{tt}+u_{xx}=0$($n=1$)を $U=(u_t,u_x)$ の系に直すと $\partial_tU_0=-\partial_xU_1$、$\partial_tU_1=\partial_xU_0$ となり、$A=\begin{pmatrix}0&-1\\1&0\end{pmatrix}$ の固有値は $\pm i$ なので、双曲型でない。

定数係数の系の初期値問題

Fourier 変換は微分 $\partial_j$ を $2\pi i\xi_j$ を掛ける操作に変える(規格化は Fourier変換 の定義に合わせ、$\widehat g(\xi)=\int g(x)e^{-2\pi ix\cdot\xi}dx$ とする)。よって $(*)$ は各周波数 $\xi$ ごとの常微分方程式 $\partial_t\widehat u=G(\xi)\widehat u$、$G(\xi):=2\pi iA(\xi)+B$ になり、解の候補は
$$ u(x,t)=\int_{\mathbb{R}^n}e^{2\pi ix\cdot\xi}\,e^{tG(\xi)}\,\widehat g(\xi)\,d\xi\qquad(**) $$
である($e^{tG}$ は 行列の指数関数)。$\mathbb{C}^N$ のノルムを $\lvert\cdot\rvert$、$L^2(\mathbb{R}^n;\mathbb{C}^N)$ のノルムを $\|\cdot\|_2$ と書く。

定数係数の系の $L^2$ 評価

$g\in\mathcal S(\mathbb{R}^n;\mathbb{C}^N)$(成分が Schwartz空間 に属する)とし、$u$ を $(**)$ で定める。

  1. 系が対称双曲型なら、すべての $\xi$ と $t$ について $\| e^{tG(\xi)}\|\le e^{\| B\|\lvert t\rvert}$ であり、$B=0$ なら $e^{tG(\xi)}$ はユニタリ行列である。
  2. 系が狭義双曲型で $B=0$ なら、$A_1,\dots,A_n$ だけで決まる定数 $K$ があって、すべての $\xi$ と $t$ について $\| e^{tG(\xi)}\|\le K$ である。
  3. 1 または 2 の場合、$u$ は $\mathbb{R}^n\times\mathbb{R}$ の上で $C^1$ 級で $(*)$ と $u(\cdot,0)=g$ を満たし、各時刻で $u(\cdot,t)\in L^2$ である。1 の場合は $\| u(\cdot,t)\|_2\le e^{\| B\|\lvert t\rvert}\| g\|_2$($B=0$ なら等号)、2 の場合は $\| u(\cdot,t)\|_2\le K\| g\|_2$ である。
周波数ごとの行列の大きさを評価する

要点:1 は、$A(\xi)$ が実対称なので $2\pi iA(\xi)$ が $\lvert w\rvert^2$ を変えないことによる。2 は、相異なる固有値への射影を Lagrange の補間で書き、固有値の間隔を判別式で下から抑える。3 は、積分記号下の微分と Plancherel の定理による。

詳しい証明を開く

1:$v\in\mathbb{C}^N$ について $w(t):=e^{tG(\xi)}v$ とおくと $w'=Gw$ なので、Hermite 内積 $\langle\cdot,\cdot\rangle$ について $\frac d{dt}\lvert w\rvert^2=2\operatorname{Re}\langle Gw,w\rangle$ である。$A(\xi)$ は実対称なので $\langle A(\xi)w,w\rangle$ は実数で(Hermite行列)、$\langle2\pi iA(\xi)w,w\rangle$ は純虚数である。よって $\frac d{dt}\lvert w\rvert^2=2\operatorname{Re}\langle Bw,w\rangle\le2\| B\|\lvert w\rvert^2$ であり、$(e^{-2\| B\| t}\lvert w\rvert^2)'\le0$ から $t\ge0$ で $\lvert w(t)\rvert\le e^{\| B\| t}\lvert v\rvert$ である。$t\le0$ では $(e^{2\| B\| t}\lvert w\rvert^2)'\ge0$ から同様である。$B=0$ なら $\lvert w\rvert$ は一定で、$e^{tG}$ はノルムを保つ線形写像、すなわちユニタリ行列である(ユニタリ行列)。

2:$N=1$ なら $\lvert e^{2\pi itA(\xi)}\rvert=1$ なので $N\ge2$ とする。$S$ を単位球面とし、$\xi\in S$ について $A(\xi)$ の固有値を $\lambda_1(\xi)< \dots< \lambda_N(\xi)$ とする。$\lvert\lambda_k\rvert\le\| A(\xi)\|\le M:=\sum_j\| A_j\|$ である。固有多項式の 判別式 $D(\xi)=\prod_{k< l}(\lambda_k-\lambda_l)^2$ は固有多項式の係数の多項式なので $\xi$ の連続関数で、$S$ の上で正だから、コンパクト性により最小値 $d>0$ をもつ。1 組 $(k,l)$ を除く $\binom N2-1$ 個の因子はそれぞれ $(2M)^2$ 以下なので、$(\lambda_k-\lambda_l)^2\ge d/(2M)^{2(\binom N2-1)}=:\delta^2$ である。固有値が相異なるので $A(\xi)$ は対角化でき、固有空間への射影は $P_k=\prod_{l\ne k}\frac{A(\xi)-\lambda_lI}{\lambda_k-\lambda_l}$(Lagrange の補間)と書けて、$e^{2\pi itA(\xi)}=\sum_ke^{2\pi it\lambda_k}P_k$ である。$\| P_k\|\le(2M/\delta)^{N-1}$ なので、$\| e^{2\pi itA(\xi)}\|\le K:=N(2M/\delta)^{N-1}$ である。$r>0$ について $A(r\xi)=rA(\xi)$ の射影は $A(\xi)$ の射影と同じで、固有値は $r$ 倍になるので、同じ評価がすべての $\xi\ne0$ で成り立つ。$\xi=0$ では $e^{0}=I$ である。

3:1・2 の評価を $\| e^{tG(\xi)}\|\le K_T$($\lvert t\rvert\le T$)とまとめる。$(**)$ の被積分関数は $K_T\lvert\widehat g\rvert$ で、その $t$ による微分 $e^{2\pi ix\cdot\xi}G(\xi)e^{tG}\widehat g$ と $x_j$ による微分 $2\pi i\xi_je^{2\pi ix\cdot\xi}e^{tG}\widehat g$ は $\bigl(2\pi(M+1)\lvert\xi\rvert+\| B\|\bigr)K_T\lvert\widehat g\rvert$ で抑えられる。$\widehat g$ は Schwartz 関数なので(Schwartz空間 の定理「Schwartz空間上のFourier変換」)これらは可積分であり、積分記号下の微分 により $u$ は $C^1$ 級で、微分は積分の中で計算できる。すると $u_t-\sum_jA_j\partial_ju-Bu=\int e^{2\pi ix\cdot\xi}\bigl(G(\xi)-2\pi iA(\xi)-B\bigr)e^{tG}\widehat g\,d\xi=0$ である。$t=0$ では $(**)$ は $\widehat g$ の逆変換なので、Fourier変換 の定理「Fourier反転公式」から $u(\cdot,0)=g$ である。$h_t:=e^{tG}\widehat g$ は $L^1\cap L^2$ に属し、Plancherelの定理 の定理「Plancherelの定理($\mathbb{R}^n$)」の 3 により $u(\cdot,t)$ は $h_t$ の $L^2$ での逆 Fourier 変換である。逆変換はユニタリなので $\| u(\cdot,t)\|_2=\| h_t\|_2$ で、これは $\| e^{tG}\|$ の上限と $\|\widehat g\|_2=\| g\|_2$ の積以下である(成分ごとに使う)。$B=0$ の対称双曲型では $\lvert h_t(\xi)\rvert=\lvert\widehat g(\xi)\rvert$ なので等号が成り立つ。$\square$

対称双曲型の系では、低階の項 $Bu$ は解の大きさを $e^{\| B\|\lvert t\rvert}$ 倍までしか変えない。狭義双曲型の系では、各周波数で $A(\xi)$ を対角化する行列の大きさが周波数によらず抑えられることが評価の鍵である。

双曲性は適切性に必要である

$B=0$ とし、ある $\xi_0\in\mathbb{R}^n$ について $A(\xi_0)$ が実数でない固有値をもつとする。このとき、どの整数 $m\ge0$ と $T>0$、$C>0$ についても、すべての導関数が有界な $C^\infty$ 級の初期値 $g$ と、$\mathbb{R}^n\times[0,T]$ の上で有界な $(*)$ の $C^\infty$ 級の解 $u$ で $u(\cdot,0)=g$ かつ
$$ \sup_x\lvert u(x,T)\rvert>C\max_{\lvert\alpha\rvert\le m}\sup_x\lvert\partial^\alpha g(x)\rvert $$
となるものがある。つまり、有限回の微分まで込めた初期値の大きさで解の大きさを抑える評価は、どれも成り立たない。

周波数の高い平面波が指数関数的に増える

$A(\xi_0)$ は実行列なので、実数でない固有値は共役な組で現れる。虚部が負の固有値 $\mu$ と固有ベクトル $v$ をとる($A(0)=0$ なので $\xi_0\ne0$)。$k>0$ について
$$ u_k(x,t):=e^{2\pi ik(x\cdot\xi_0+\mu t)}v $$
とおくと、$\partial_tu_k=2\pi ik\mu u_k$、$\sum_jA_j\partial_ju_k=2\pi ik\,e^{2\pi ik(x\cdot\xi_0+\mu t)}A(\xi_0)v=2\pi ik\mu u_k$ なので $u_k$ は $(*)$ の解である。$\lvert u_k(x,t)\rvert=e^{2\pi kt\lvert\operatorname{Im}\mu\rvert}\lvert v\rvert$ なので $[0,T]$ の上で有界である。一方 $g_k:=u_k(\cdot,0)$ は $\lvert\partial^\alpha g_k\rvert\le(2\pi k\lvert\xi_0\rvert)^{\lvert\alpha\rvert}\lvert v\rvert$ を満たす。よって
$$ \frac{\sup_x\lvert u_k(x,T)\rvert}{\max_{\lvert\alpha\rvert\le m}\sup_x\lvert\partial^\alpha g_k\rvert}\ge\frac{e^{2\pi kT\lvert\operatorname{Im}\mu\rvert}}{\max(1,2\pi k\lvert\xi_0\rvert)^m}\to\infty\qquad(k\to\infty) $$
であり、$k$ を大きくとればよい。$\square$

Laplace 方程式の系(ex-hpde-wave-system の後半)では $\mu=\pm i$ で、解は $e^{2\pi k t}$ の程度で増える。これは 偏微分方程式 の例「Laplace 方程式の Cauchy 問題」の Hadamard の例と同じ現象である。

反例:条件を外すと成り立たないこと

固有値が実数であること(双曲性)は必要だが、それだけでは足りない。下の 2 行目と 3 行目は、双曲型ではあるが対称双曲型でも狭義双曲型でもない系($A(\xi)$ が対角化できない)である。

外す条件反例成り立たなくなること
双曲性(thm-hpde-fourier の仮定)Laplace 方程式の系 $A=\begin{pmatrix}0&-1\\1&0\end{pmatrix}$($n=1$、ex-hpde-wave-system)有限回の微分を込めた初期値の大きさによる解の評価(thm-hpde-necessity)
対称性・狭義双曲性(固有値は実数)$n=1$、$A=\begin{pmatrix}0&1\\0&0\end{pmatrix}$、$B=0$$\| u(\cdot,t)\|_2\le K\| g\|_2$
対称性・狭義双曲性(低階の項がある場合)$n=1$、同じ $A$ と $B=\begin{pmatrix}0&0\\1&0\end{pmatrix}$有限回の微分を込めた初期値の大きさによる解の評価
円錐の傾きが $c=\sqrt\Lambda$ であること(thm-hpde-finite-speed)$n=1$、$u_{tt}=c^2u_{xx}$ の解 $u=\varphi(x+ct)$、傾き $c'< c$ の円錐初期値が $\overline B(x_0,R)$ で $0$ なら円錐の上で $u=0$
表の反例の確かめ

2 行目:系は $\partial_tu_1=\partial_xu_2$、$\partial_tu_2=0$ で、$A(\xi)=\xi A$ の固有値は $0$ だけ(実数)なので双曲型である。解は $u_2=g_2$、$u_1=g_1+tg_2'$ であり、初期値の 1 階微分が解に現れる。

評価が成り立たないことの確かめを開く

$0$ でない $\varphi\in\mathcal S(\mathbb{R})$ をとり、$g=(0,\varphi(kx))$($k>0$)とおく。置換積分で $\| g\|_2=k^{-1/2}\|\varphi\|_2$、$\| u_1(\cdot,t)\|_2=t\,k\,\|\varphi'(k\,\cdot)\|_2=t\,k^{1/2}\|\varphi'\|_2$ なので、$t>0$ で $\| u(\cdot,t)\|_2/\| g\|_2\ge tk\|\varphi'\|_2/\|\varphi\|_2\to\infty$($k\to\infty$)である($\varphi'\ne0$。$\varphi'=0$ なら $\varphi$ は定数で、可積分なので $0$ になる)。一方 $\| u(\cdot,t)\|_2\le(1+\lvert t\rvert)(\| g\|_2+\| g'\|_2)$ は成り立つので、この系は微分を 1 回失う評価なら満たす。


3 行目:系は $\partial_tu_1=\partial_xu_2$、$\partial_tu_2=u_1$ で、$u_2$ は $\partial_t^2u_2=\partial_xu_2$ を満たす。主要部は 2 行目と同じなので双曲型である。
解が速く増えることの確かめを開く

$k>0$ について $\lambda:=\sqrt{2\pi k}\,e^{i\pi/4}$ とおくと $\lambda^2=2\pi ik$ である。$v:=(\lambda,1)$、$u(x,t):=e^{2\pi ikx+\lambda t}v$ とおくと、$\partial_tu_1=\lambda^2e^{2\pi ikx+\lambda t}=2\pi ik\,e^{2\pi ikx+\lambda t}=\partial_xu_2$、$\partial_tu_2=\lambda e^{2\pi ikx+\lambda t}=u_1$ なので解である。$\operatorname{Re}\lambda=\sqrt{\pi k}$ なので $\lvert u(x,T)\rvert=e^{T\sqrt{\pi k}}\lvert v\rvert$ である。初期値の $j$ 階の微分の大きさは $(2\pi k)^j\lvert v\rvert$ なので、どの $m$ についても $e^{T\sqrt{\pi k}}/(2\pi k)^m\to\infty$($k\to\infty$)となり、thm-hpde-necessity と同じ形の評価はどれも成り立たない。


対称双曲型の系では低階の項は解の大きさを $e^{\| B\|\lvert t\rvert}$ 倍までしか変えないが(thm-hpde-fourier の 1)、対角化できない系では低階の項が適切性を壊しうる。
4 行目:$\Lambda=c^2$ である。$\varphi$ を $C^2$ 級で、台が区間 $(R,R+\delta)$ に含まれ、$\varphi(R+\delta/2)\ne0$ となる関数とし、$x_0=0$ とする。$u=\varphi(x+ct)$ は左へ進む波で、初期値 $u(\cdot,0)=\varphi$、$u_t(\cdot,0)=c\varphi'$ は $[-R,R]$ の上で $0$ である。$t:=\frac{\delta}{2(c-c')}$ が $t< R/c'$ を満たすように $\delta$ を小さくとると、点 $x=R+\delta/2-ct$ は $\lvert x\rvert\le R-c't$ を満たし($x\le R-c't$ は $t$ の選び方から等号で成り立ち、$x\ge-(R-c't)$ も $\delta$ が小さければ成り立つ)、そこで $u=\varphi(R+\delta/2)\ne0$ である。

補足

  • 文献:弱解の存在は Bre11 第 10 章の注 6(pp. 345–346)、定数係数の波動方程式の初期境界値問題は同書 §10.3(pp. 335–340)にある。有限伝播速度とエネルギー評価、Fourier 変換による定数係数の系の扱いは、この記事では独自に証明した。
  • 解析的なデータ:係数・右辺・データが解析的なら、非特性の超平面の上の Cauchy 問題は型によらず局所的に解析的な解をもつ(Cauchy–Kovalevskayaの定理)。ただしデータへの連続な依存は与えないので、解析的でないデータを扱えるこの記事のエネルギー評価とは役割が違う。

関連項目

参考文献

[1]
Haim Brezis, Functional Analysis, Sobolev Spaces and Partial Differential Equations, Springer, 2011, §10.3(pp. 335–340、波動方程式の初期境界値問題。Theorem 10.7(p. 336)、Theorem 10.8(p. 336、解の滑らかさ)、Remark 8(p. 339、平滑化がないこと))、第 10 章の注 6 Theorem 10.14(J.-L. Lions の定理)とその応用(pp. 345–346、係数が変わる双曲型方程式の弱解)

Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する