非線形拡散方程式における質量保存則と解の有限時刻消滅について

4段階に分けて証明を行う.
Step:1 $u_0$は非負かつ$C_0^{\mathrm{\infty }}(\mathrm{\Omega })$と仮定する.
方程式の退化性または特異性を回避するため,
$$u_{0,n}=u_0+\frac{1}{n}$$
とした初期値境界値問題
$$\begin{array}{}\text{(7)}& \{\begin{array}{ll}{\partial }_t{u}_n=\mathrm{\Delta }{u}_n^m& \mathrm{i}\mathrm{n}\mathrm{\Omega }\times (0,\mathrm{\infty })=:Q,\\ u_n={\displaystyle \frac{1}{n}}& \mathrm{o}\mathrm{n}\partial \mathrm{\Omega }\times [0,\mathrm{\infty })=:\mathrm{\Sigma },\\ u_n=u_{0,n}& \mathrm{o}\mathrm{n}\mathrm{\Omega }\times \{t=0\}\end{array}\end{array}$$
を考える. 問題$\text{(7)}$は放物型方程式における最大値原理より
$$\begin{array}{}\text{(8)}& {\displaystyle \frac{1}{n}}⩽u_n(x,t)⩽M+{\displaystyle \frac{1}{n}}\mathrm{i}\mathrm{n}\bar{Q}\end{array}$$
であると予想される. ただし, $M=sup(u_0)$である. ここで, 問題$\text{(7)}$の方程式は
$${\partial }_t{u}_n=\mathrm{\nabla }\cdot (a(u_n)\mathrm{\nabla }{u}_n)$$
と書き換えられ, $u_n$が式$\text{(8)}$をみたしているとき$a(u_n)=m{u}_n^{m-1}$とする. 故に, $a(u_n)$は正値かつ滑らかで, $a(u_n)⩾c>0$である. このとき, 問題$\text{(7)}$はstandard quasilinear theory[cf. 1, Chapter 6]より古典解$u_n\in C^{2,1}(\bar{Q})\cap C^{\mathrm{\infty }}(Q)$を一意にもつことがわかる. よって, 放物型方程式における最大値原理を適用すると, 式$\text{(8)}$が成立する. 実際,
$${\partial }_t{u}_n=a(u_n)\mathrm{\Delta }{u}_n+\mathrm{\nabla }a(u_n)\cdot \mathrm{\nabla }{u}_n$$
であるので, $\hat{u}=u_n-(M+{\displaystyle \frac{1}{n}})$とおけば,
$${\partial }_t\hat{u}=a(u_n)\mathrm{\Delta }\hat{u}+\mathrm{\nabla }a(u_n)\cdot \mathrm{\nabla }\hat{u}$$
と$\hat{u}$に関する線型方程式になる. したがって, 最大値原理より${\hat{u}}_0⩽0$より$\hat{u}⩽0$すなわち$u_n⩽M+{\displaystyle \frac{1}{n}}$がしたがう. 同様にして$u_n⩾{\displaystyle \frac{1}{n}}$であるから, 古典解$u_n$は式$\text{(8)}$をみたす. さらに, $n$は任意なので, 同様に
$$\frac{1}{n+1}⩽u_{n+1}⩽M+\frac{1}{n+1}\mathrm{i}\mathrm{n}\bar{Q}$$
が成立する. ここで,
$$w(x,t)=u_n(x,t)-u_{n+1},$$
$$A(x,t)=m(m-1){\int }_0^1\{u_{n+1}(x,t)+\tau (u_n(x,t)-u_{n+1}(x,t)){\}}^{m-2}d\tau$$
とおくと,
$$\begin{array}{rl}{\partial }_{t}w(x,t)& =\mathrm{\nabla }\cdot (a_n(u_n)\mathrm{\nabla }{u}_n(x,t))-\mathrm{\nabla }\cdot (a_{n+1}(u_{n+1}(x,t))\mathrm{\nabla }{u}_{n+1}(x,t))\\ & =\mathrm{\nabla }\cdot (m{u}_n^{m-1}(x,t)\mathrm{\nabla }(u_n(x,t)-u_{n+1}(x,t)))\\ & +\mathrm{\nabla }\cdot ((m{u}_n^{m-1}(x,t)-m{u}_{n+1}^{m-1}(x,t))\mathrm{\nabla }{u}_{n+1}(x,t))\\ & =\mathrm{\nabla }\cdot (m{u}_n^{m-1}(x,t)\mathrm{\nabla }w(x,t))+\mathrm{\nabla }\cdot (A(x,t)w(x,t)\mathrm{\nabla }{u}_{n+1}(x,t))\\ & =:p(x,t)\mathrm{\Delta }w+\mathrm{\nabla }q(x,t)\cdot \mathrm{\nabla }w(x,t)+r(x,t)w(x,t)\end{array}$$
となる. ここで,
$$\begin{array}{rl}p(x,t)& :=m{u}_n^{m-1},\\ q(x,t)& :=m{u}_n^{m-1}+A(x,t)u_{n+1}(x,t),\\ r(x,t)& :=\mathrm{\nabla }\cdot (A(x,t)\mathrm{\nabla }{u}_{n+1})\end{array}$$
である. したがって, $p(x,t)>0$であるので, 最高階の微分の係数が退化しない$w(x,t)$に関する線型方程式に帰着できる. よって, $w_0=u_{0,n}-u_{0,n+1}$なので, 最大値原理より$w⩾0$ in $\bar{Q}$すなわち
$$u_{n+1}⩽u_n\mathrm{i}\mathrm{n}\bar{Q}$$
が得られる. 故に, $\{u_n\}$は下に有界な単調減少列であるので, 各点の意味で
$$u(x,t)=\underset{n\to \mathrm{\infty }}{lim}{u}_n,(x,t)\in \bar{Q}$$
となる$u$が存在し, $0⩽u⩽M$である. 一方, $\mathrm{\Sigma }$上では,
$$\underset{\mathrm{\Sigma }}{sup}u⩽\underset{\mathrm{\Sigma }}{sup}{u}_n={\displaystyle \frac{1}{n}}\to 0\mathrm{a}\mathrm{s}n\to \mathrm{\infty }$$
となる.

Step2: $L^{\mathrm{\infty }}(0,\mathrm{\infty };L^1(\mathrm{\Omega }))$上のCauchy列の構成.
$u\in L^{\mathrm{\infty }}(0,\mathrm{\infty };L^1({\mathbb{R}}^N))$となる$u$の存在性を証明するため, 古典近似解に対して次の$L^1$-contraction principleと呼ばれる不等式
$$\begin{array}{}\text{(9)}& \parallel u_n(t)-v_n(t){\parallel }_{L^1(\mathrm{\Omega })}⩽\parallel u_{n,0}-v_{n,0}{\parallel }_{L^1(\mathrm{\Omega })}\end{array}$$
が成立することを確かめる.
$p_k\in C^1(\mathbb{R})$を
$${\mathrm{s}\mathrm{i}\mathrm{g}\mathrm{n}}_0^{+}(s)=\{\begin{array}{ll}1& \mathrm{i}\mathrm{f}s⩾0,\\ 0& \mathrm{i}\mathrm{f}s<0\end{array}$$
の近似列とする. すなわち, $k\in \mathbb{N}$として
$$p_k(s)=\{\begin{array}{ll}1& \mathrm{i}\mathrm{f}s⩾{\displaystyle \frac{1}{k}},\\ 0& \mathrm{i}\mathrm{f}s<0,\end{array}0<p_k⩽1\mathrm{i}\mathrm{f}0<s<{\displaystyle \frac{1}{k}}$$
をみたす非負かつ非減少関数とする. 問題$\text{(7)}$の古典近似解を$u_n,v_n$とし,
$${\partial }_t(u_n-v_n)=\mathrm{\Delta }(u_n^m-v_n^m)$$
の両辺に$p_k(u_n^m-v_n^m)$をかけて$\mathrm{\Omega }$上で積分すると,
$$\begin{array}{rl}{\int }_{\mathrm{\Omega }}{p}_k(u_n^m(x,t)-v_n^m(x,t)){\partial }_t(u_n-v_n)(x,t)dx& ={\int }_{\mathrm{\Omega }}{p}_k(u_n^m(x,t)-v_n^m(x,t))\mathrm{\Delta }(u_n^m-v_n^m)(x,t)dx\\ & =-{\int }_{\mathrm{\Omega }}|\mathrm{\nabla }(u_n^m-v_n^m)(x,t)|p_k^{\prime }(u_n^m(x,t)-v_n^m(x,t))dx\\ & ⩽0\end{array}$$
となる. したがって,
$$p_k(u_n^m-v_n^m){\partial }_t(u_n-v_n)\to {\mathrm{s}\mathrm{i}\mathrm{g}\mathrm{n}}_0^{+}(u_n-v_n){\partial }_t(u_n-v_n)={\partial }_t[u_n-v_n{]}_{+}\mathrm{a}\mathrm{s}k\to \mathrm{\infty }$$
が成立するので,
$$\frac{d}{dt}{\int }_{\mathrm{\Omega }}[u_n-v_n{]}_{+}(x,t)dx$$
を得る. 故に, 両辺を$0$から$t$まで積分すれば,
$${\int }_{\mathrm{\Omega }}[u_n-v_n{]}_{+}(x,t)dx⩽{\int }_{\mathrm{\Omega }}[u_{n,0}-v_{n,0}{]}_{+}(x,t)dx$$
となるので, $u_n$と$v_n$を入れ替え同様の議論をすれば, 式$\text{(9)}$が成立する. これにより, 初期値にmollifier$J_{\epsilon }$を作用させ,
$$u_{0,n}(x)=(J_{1/n}\ast u_0)(x)+\frac{1}{n}$$
とした初期値境界値問題$\text{(7)}$の古典近似解と任意の$n,k\in \mathbb{N}$に対して
$$\begin{array}{rl}\parallel u_n(t)-u_k(t){\parallel }_{L^1(\mathrm{\Omega })}& ⩽\parallel u_{n,0}-u_{k,0}{\parallel }_{L^1(\mathrm{\Omega })}\\ & ⩽\parallel (J_{1/n}\ast u_0)-u_0{\parallel }_{L^1(\mathrm{\Omega })}+\parallel (J_{1/k}\ast u_0)-u_0{\parallel }_{L^1(\mathrm{\Omega })}+(\frac{1}{n}-\frac{1}{k})|\mathrm{\Omega }|\\ & \to 0\mathrm{a}\mathrm{s}n,k\to \mathrm{\infty }\end{array}$$
となる. よって, $\{u_n\}$は$L^{\mathrm{\infty }}(0,\mathrm{\infty };L^1(\mathrm{\Omega }))$上のCauchy列である. $L^{\mathrm{\infty }}(0,\mathrm{\infty };L^1(\mathrm{\Omega }))$はBanach空間であるから,
$$u_n\to u\mathrm{s}\mathrm{t}\mathrm{r}\mathrm{o}\mathrm{n}\mathrm{g}\mathrm{l}\mathrm{y}\mathrm{i}\mathrm{n}{L}^{\mathrm{\infty }}(0,\mathrm{\infty };L^1(\mathrm{\Omega }))\mathrm{a}\mathrm{s}n\to \mathrm{\infty }$$
となる$u\in L^{\mathrm{\infty }}(0,\mathrm{\infty };L^1(\mathrm{\Omega }))$が一意に存在する.

Step3: $L^1$-連続性の証明.
今得られた$u$に対しても,
$$\begin{array}{}\text{(10)}& \parallel u(t)-v(t){\parallel }_{L^1(\mathrm{\Omega })}⩽\parallel u_0-v_0{\parallel }_{L^1(\mathrm{\Omega })}\end{array}$$
が成立するので,
$$\begin{array}{rl}\parallel u(t)-u_0{\parallel }_{L^1(\mathrm{\Omega })}& ⩽\parallel u(t)-u_n(t){\parallel }_{L^1(\mathrm{\Omega })}+\parallel u_n(t)-u_0{\parallel }_{L^1(\mathrm{\Omega })}+\parallel u_{n,0}-u_0{\parallel }_{L^1(\mathrm{\Omega })}\\ & ⩽2\parallel u_0-u_{n,0}{\parallel }_{L^1(\mathrm{\Omega })}+\parallel u_n(t)-u_0{\parallel }_{L^1(\mathrm{\Omega })}\\ & \to 2\parallel u_0-u_{n,0}{\parallel }_{L^1(\mathrm{\Omega })}\mathrm{a}\mathrm{s}t\to 0,\\ & \to 0\mathrm{a}\mathrm{s}n\to \mathrm{\infty }\end{array}$$
を得る. 故に$t=0$での連続性が得られる. $t>0$の場合においても同様に計算すればよい. 以上より$u\in C([0,\mathrm{\infty });L^1(\mathrm{\Omega }))\cap L^{\mathrm{\infty }}(Q)$が証明された. $◻$

$u<v$ a.e. in $Q$と仮定する. このとき, $L^1$-contraction principleより
$${\int }_{\mathrm{\Omega }}[v(x,t)-u(x,t){]}_{+}dx⩽{\int }_{\mathrm{\Omega }}[v_0(x)-u_0(x){]}_{+}dx=0$$
となる. これは$u<v$に矛盾する. 故に$u⩾v$ a.e. in $Q$である. $◻$

$p_k\in C^1(\mathbb{R})$を前章と同様のものとし,
$$j_k(r)={\int }_0^r{p}_k(s)ds$$
とする. すなわち$j_k$は関数$s\mapsto [s{]}_{+}$の近似列である. さらに, ${\xi }_0\in C_0^{\mathrm{\infty }}({\mathbb{R}}^N)$を$0⩽{\xi }_0⩽1$,
$${\xi }_0(x)=\{\begin{array}{ll}1& \mathrm{i}\mathrm{f}|x|⩽1,\\ 0& \mathrm{i}\mathrm{f}|x|⩾2\end{array}$$
をみたすカットオフ関数とし, ${\xi }_n(x)={\xi }_0(x/n)$と定める. つまり
$${\xi }_n\to 1\mathrm{a}\mathrm{s}n\to \mathrm{\infty }$$
である. 方程式
$${\partial }_t(u_1-u_2)=\mathrm{\Delta }(u_1^m-u_2^m)$$
の両辺に$p_k(u_1^m-u_2^m){\xi }_n$をかけて${\mathbb{R}}^N$上で積分すると,
$${\int }_{{\mathbb{R}}^N}{p}_k(u_1^m(x,t)-u_2^m(x,t)){\xi }_n(x){\partial }_t(u_1(x,t)-u_2(x,t))dx={\int }_{{\mathbb{R}}^N}{p}_k(u_1^m(x,t)-u_2^m(x,t)){\xi }_n(x)\mathrm{\Delta }(u_1^m(x,t)-u_2^m(x,t))dx$$
となる. 右辺は
$$\begin{array}{rl}& {\int }_{{\mathbb{R}}^N}{p}_k(u_1^m(x,t)-u_2^m(x,t)){\xi }_n(x)\mathrm{\Delta }(u_1^m(x,t)-u_2^m(x,t))dx\\ & =-{\int }_{{\mathbb{R}}^N}{p}_k^{\prime }(u_1^m(x,t)-u_2^m(x,t))|\mathrm{\nabla }(u_1^m(x,t)-u_2^m(x,t)){|}^2{\xi }_n(x)dx-{\int }_{{\mathbb{R}}^N}{p}_k(u_1^m(x,t)-u_2^m(x,t))\mathrm{\nabla }(u_1^m(x,t)-u_2^m(x,t))\cdot \mathrm{\nabla }{\xi }_n(x)dx\\ & ⩽-{\int }_{{\mathbb{R}}^N}{p}_k(u_1^m(x,t)-u_2^m(x,t))\mathrm{\nabla }(u_1^m(x,t)-u_2^m(x,t))\cdot \mathrm{\nabla }{\xi }_n(x)dx\\ & =-{\int }_{{\mathbb{R}}^N}\mathrm{\nabla }{j}_k(u_1^m(x,t)-u_2^m(x,t))\cdot \mathrm{\nabla }{\xi }_n(x)dx\\ & ={\int }_{{\mathbb{R}}^N}{j}_k(u_1^m(x,t)-u_2^m(x,t))\mathrm{\Delta }{\xi }_n(x)dx\\ & ⩽\frac{1}{n^2}{\int }_{{\mathbb{R}}^N}(u_1^m(x,t)-u_2^m(x,t))\mathrm{\Delta }\xi (x/n)dx\end{array}$$
と評価できる. ここで, $m>1$のときは$u(t)\in L^1({\mathbb{R}}^N)$であれば$u^m(t)\in L^1({\mathbb{R}}^N)$であるので
$$\begin{array}{rl}\frac{1}{n^2}{\int }_{{\mathbb{R}}^N}(u_1^m(x,t)-u_2^m(x,t))\mathrm{\Delta }\xi (x/n)dx& ⩽\frac{1}{n^2}\parallel \mathrm{\Delta }\xi {\parallel }_{L^{\mathrm{\infty }}({\mathbb{R}}^N)}{\int }_{{\mathbb{R}}^N}|u_1^m(x,t)-u_2^m(x,t)|dx\\ & \to 0\mathrm{a}\mathrm{s}n\to \mathrm{\infty }\end{array}$$
となる. 一方, $0<m<1$のときは$(a-b{)}^m⩽2^{1-m}(a-b{)}^m$を用いると
$$\begin{array}{rl}\frac{1}{n^2}{\int }_{{\mathbb{R}}^N}(u_1^m(x,t)-u_2^m(x,t))\mathrm{\Delta }\xi (x/n)dx& ⩽\frac{2^{1-m}}{n^2}{\int }_{{\mathbb{R}}^N}(u_1(x,t)-u_2(x,t){)}^m\mathrm{\Delta }\xi (x/n)dx\\ & ⩽\frac{2^{1-m}}{n^2}{({\int }_{{\mathbb{R}}^N}|u_1(t)-u_2(t)|dx)}^m{({\int }_{{\mathbb{R}}^N}|\mathrm{\Delta }{\xi }_n(x){|}^{\frac{1}{1-m}})}^{1-m}\\ & \to 0\mathrm{a}\mathrm{s}n\to \mathrm{\infty }\end{array}$$
となる. 以上の計算を整理すると, $n,k\to \mathrm{\infty }$とすれば${\mathrm{s}\mathrm{i}\mathrm{g}\mathrm{n}}_0^{+}(u_1-u_2){\displaystyle \frac{d}{dt}}(u_1-u_2)={\displaystyle \frac{d}{dt}}[u_1-u_2{]}_{+}$であることから,
$$\frac{d}{dt}{\int }_{{\mathbb{R}}^N}[u_1(x,t)-u_2(x,t){]}_{+}dx⩽0$$
を得る. よって, 両辺を$0$から$t$まで積分すると,
$$\parallel [u_1(t)-u_2(t){]}_{+}{\parallel }_{L^1({\mathbb{R}}^N)}⩽\parallel [u_{0,1}-u_{0,2}{]}_{+}{\parallel }_{L^1({\mathbb{R}}^N)}$$
となるので, $u_1,u_2$を入れ替えて同様の議論をすれば補題の証明が完了する.$◻$