藤田指数(Fujita exponent)とは、全空間 $\mathbb{R}^n$ の上の半線形熱方程式 $u_t=\Delta u+u^p$($p>1$)の非負の解について、有限時間の爆発と時間大域的な存在の境目になる指数 $p_F=1+2/n$ である。$1<p<p_F$ なら恒等的に $0$ でない非負の解はすべて有限時間で爆発し、$p>p_F$ なら熱核の小さな定数倍以下の初期値から出発した解は時間大域的に存在して熱方程式の解と同じ速さで減衰する。臨界の $p=p_F$ でも解は爆発する。1966 年に藤田宏が示した。
前提知識: 熱方程式, Banachの不動点定理, Jensenの不等式, Lebesgue積分
藤田指数(Fujita exponent)は、全空間 $\mathbb{R}^n$ の上の半線形熱方程式
$$
u_t=\Delta u+u^p\qquad(p>1)
$$
の非負の解について、「どんなに小さな初期値から出発しても有限時間で爆発する」場合と「熱核の小さな定数倍で押さえられる初期値なら時間大域的に解ける」場合の境目にあたる指数
$$
p_F:=1+\frac2n
$$
である。$1< p< p_F$ なら、恒等的に $0$ でない非負の解はすべて有限時間で爆発し、$p>p_F$ なら、熱核の小さな定数倍で押さえられる初期値に対しては解が時間大域的に存在して熱方程式の解と同じ速さで減衰する。境目の $p=p_F$ では、恒等的に $0$ でない連続な初期値からの解はすべて爆発する側に入る。
1966 年に藤田宏(Hiroshi Fujita)がこの二分法を証明した(Fuj66 Theorem 1・Theorem 2、p. 111)。藤田は序文で、Navier–Stokes 方程式の大域解の存在が空間の次元 2 と 3 で大きく異なることを挙げ、非線形の初期値問題で「空間の次元」と「非線形性の強さ」が組み合わさって大域的な振る舞いを決める様子を、この方程式で示すことを目的に掲げている(Fuj66 §1、pp. 109–110)。
$p_F$ が現れる理由は、次の 2 つの効果の競合として理解できる。空間の変化を無視した常微分方程式 $y'=y^p$ は、正の初期値から必ず有限時間で爆発する。一方、熱方程式は質量を保ったまま解を広げ、$\int u_0\,dx$ が有限な初期値に対して最大値は $t^{-n/2}$ の速さで減衰する。非線形項が解を押し上げる効果は、時間 $t$ の付近で $u^{p-1}\approx t^{-n(p-1)/2}$ に比例し、これが $t$ について可積分($n(p-1)/2>1$)なら拡散が勝ち、そうでなければ非線形項が勝つ。この記事では、mild 解の局所存在と一意性を示したあと、この競合を 2 つの不等式(prop-fe-lifespan と thm-fe-global)として正確に述べて証明する。
以下 $n\ge1$、$p>1$ とし、$\Phi(x,t)=(4\pi t)^{-n/2}e^{-|x|^2/(4t)}$ を熱核とする(熱方程式 の定義「熱核」)。$\mathbb{R}^n$ 上の有界な Borel 可測関数 $g$ について
$$
S(t)g(x):=\int_{\mathbb{R}^n}\Phi(x-y,t)\,g(y)\,dy\quad(t>0),\qquad S(0)g:=g
$$
とおき、$\|g\|_\infty:=\sup_x|g(x)|$ と書く。$g$ が連続なら $S(t)g$ は熱方程式の解である(熱方程式 の定理「熱方程式の Cauchy 問題の解」)。
$g,h$ を有界な Borel 可測関数、$s,t\ge0$ とする。
1 と 2 は $\Phi>0$ と $\int\Phi(y,t)\,dy=1$ から従い、4 は $\Phi(x-y,t)\le(4\pi t)^{-n/2}$ から従う。3 は Fubini の定理で積分の順序を入れ替え($\Phi$ の積分が $1$ で $g$ は有界なので絶対収束する)、熱方程式 の命題「熱核の性質」の 3($\int\Phi(x-z,t)\Phi(z-y,s)\,dz=\Phi(x-y,t+s)$)を使えばよい。$s=0$ または $t=0$ の場合は明らかである。$\square$
非線形項を熱半群で積分した形で方程式を書いたものが次の mild 解である。
$u_0$ を $\mathbb{R}^n$ 上の有界な Borel 可測関数で $u_0\ge0$ とし、$0< T\le\infty$ とする。Borel 可測関数 $u\colon\mathbb{R}^n\times[0,T)\to[0,\infty)$ が、各 $T'< T$ について $\mathbb{R}^n\times[0,T']$ で有界であり、すべての $(x,t)\in\mathbb{R}^n\times[0,T)$ で
$$
u(x,t)=S(t)u_0(x)+\int_0^t\!\!\int_{\mathbb{R}^n}\Phi(x-y,t-s)\,u(y,s)^p\,dy\,ds
$$
を満たすとき、$u$ を初期値 $u_0$ の mild 解($[0,T)$ 上の)という。右辺の第 2 項を $\int_0^tS(t-s)[u(s)^p](x)\,ds$ と書き、右辺全体を $\Gamma_{u_0}[u](x,t)$ と書く。
被積分関数は非負なので、2 重積分は Tonelli の定理により意味をもち、$u$ が Borel 可測なら $\Gamma_{u_0}[u]$ も Borel 可測である。初期値が 2 階までの有界連続な導関数をもつとき、有界な連続 mild 解は偏微分方程式の古典解であり、逆に有界な古典解は mild 解である(Fuj66 付録 Proposition A1(p. 118)・Proposition A4(pp. 121–124)。この記事では証明しない)。以下では解をすべて mild 解の意味で扱う。
次の恒等式は、時刻 $\tau$ の値を初期値として解き直しても同じ解が得られることを表す。
$u$ を初期値 $u_0$ の $[0,T)$ 上の mild 解、$0\le\tau< T$ とする。
1.$u(\tau)=S(\tau)u_0+\int_0^\tau S(\tau-s)[u(s)^p]\,ds$ に $S(r)$ を当て、Tonelli の定理で $S(r)$ を $s$ の積分の中に入れて lem-fe-semigroup の 3 を使う。
2 と 3.$t>\tau$ について、1 を $r=t-\tau$ で使うと
$$
S(t-\tau)[u(\tau)]+\int_0^{t-\tau}S(t-\tau-\sigma)[u(\tau+\sigma)^p]\,d\sigma=S(t)u_0+\int_0^tS(t-s)[u(s)^p]\,ds
$$
である(第 2 項で $s=\tau+\sigma$ と置換した)。2 では右辺が $u(t)$ なので左辺は $v(t-\tau)$ の mild 解の式になる。3 では同じ計算を、左辺の $u(\tau+\sigma)$ を $w(\tau+\sigma)=v(\sigma)$ に、右辺の $u$ を $w$ に置き換えて行う($s\le\tau$ では $w=u$ なので 1 が使える)。すると左辺は $v(t-\tau)=w(t)$、右辺は $\Gamma_{u_0}[w](t)$ になる。$t\le\tau$ では $w=u$ なので明らかである。$\square$
$u_0$ を有界な Borel 可測関数で $0\le u_0\le M$($M>0$)とし、$T_0:=\dfrac{M^{1-p}}{p\,2^p}$ とおく。
1.$X$ を、$\mathbb{R}^n\times[0,T_0]$ 上の Borel 可測関数で $0\le u\le2M$ を満たすもの全体とし、上限ノルムの距離を入れる。Borel 可測関数の一様極限は Borel 可測なので $X$ は完備である。$u\in X$ なら $\Gamma_{u_0}[u]$ は Borel 可測で、lem-fe-semigroup の 1・2 から $0\le\Gamma_{u_0}[u]\le M+T_0(2M)^p\le2M$ である($T_0(2M)^p=M/p< M$)。$a,b\in[0,2M]$ なら平均値の定理から $|a^p-b^p|\le p(2M)^{p-1}|a-b|$ なので、$u,v\in X$ について
$$
\bigl|\Gamma_{u_0}[u]-\Gamma_{u_0}[v]\bigr|(x,t)\le\int_0^tS(t-s)\bigl[p(2M)^{p-1}\|u-v\|\bigr]ds\le T_0\,p\,(2M)^{p-1}\|u-v\|=\tfrac12\|u-v\|
$$
である。Banachの不動点定理 の定理「Banach の不動点定理」により不動点 $u\in X$ がある。これが求める mild 解である。
2.$u,v$ をその 2 つの解とし、$T'< T$ を固定して、$\mathbb{R}^n\times[0,T']$ で $u,v\le K$ とする。$L:=pK^{p-1}$、$D(t):=\sup\{|u(x,s)-v(x,s)|: x\in\mathbb{R}^n,\ 0\le s\le t\}$ とおく。$t\le T'$ なら、上と同じ評価で $|u-v|(x,t)\le\int_0^tL\,D(t)\,ds=Lt\,D(t)$ なので $D(t)\le Lt\,D(t)$ であり、$t< 1/L$ なら $D(t)=0$ である。$\tau:=1/(2L)$ とすると $u=v$ が $[0,\tau]$ で成り立ち、lem-fe-shift の 2 により $u(\cdot+\tau)$ と $v(\cdot+\tau)$ は同じ初期値 $u(\tau)=v(\tau)$ の mild 解で、$[0,T'-\tau]$ で同じ $K$ で押さえられる。同じ議論を繰り返すと有限回で $[0,T']$ 全体で $u=v$ となる。$T'< T$ は任意である。$\square$
初期値 $u_0$ について、mild 解の存在する $T$ の上限を 最大存在時間 といい $T^*(u_0)\in(0,\infty]$ と書く。一意性(thm-fe-local の 2)と貼り合わせ(lem-fe-shift の 3)により、$[0,T^*)$ 上の mild 解 $u$ がただ 1 つある。これを単に初期値 $u_0$ の解という。$T^*=\infty$ のとき解は 時間大域的 であるといい、$T^*<\infty$ のとき解は 有限時間で爆発する という。
$T^*(u_0)<\infty$ なら $\displaystyle\lim_{t\to T^*}\ \sup_{x}u(x,t)=\infty$ である。
そうでなければ、$t_k\to T^*$ となる列と $M>0$ で $\sup_xu(x,t_k)\le M$ となるものがある。$T_0=T_0(M)$ を thm-fe-local のとおりとし、$t_k>T^*-T_0$ となる $k$ をとる。初期値 $u(\cdot,t_k)$ の $[0,T_0]$ 上の mild 解を lem-fe-shift の 3 で $u$ に貼り合わせると、$[0,t_k+T_0)$ 上の mild 解が得られ、$t_k+T_0>T^*$ は最大存在時間の定義に反する。$\square$
$u_0\equiv c>0$ とする。$y(t):=\bigl(c^{1-p}-(p-1)t\bigr)^{-1/(p-1)}$ は $y'=y^p$、$y(0)=c$ を満たし、$S(t)c=c$ なので、$u(x,t):=y(t)$ について
$$
\Gamma_{u_0}[u](x,t)=c+\int_0^ty(s)^p\,ds=y(t)
$$
である。よって $u$ は $[0,c^{1-p}/(p-1))$ 上の mild 解で、$t\to c^{1-p}/(p-1)$ で $y(t)\to\infty$ である。したがって $T^*(c)=c^{1-p}/(p-1)$ であり、$p$ と $n$ によらず爆発する。この初期値は可積分でなく、拡散によって減衰する余地がない。
次の不等式が爆発の側の議論の核である。証明では、$u(t)$ を熱核の重みで平均した量 $w(t)$ が常微分不等式 $w'\ge w^p$ を満たすことを使う。
$u$ を初期値 $u_0$ の $[0,T)$ 上の mild 解とし、$0< T<\infty$ とする。このとき、すべての $x_0\in\mathbb{R}^n$ について
$$
(p-1)\,T\,\bigl(S(T)u_0(x_0)\bigr)^{p-1}\le1
$$
である。とくに $T^*(u_0)>T$ なら、すべての $x_0$ でこの不等式が成り立つ。
$a:=S(T)u_0(x_0)$ とおく。$a=0$ なら示すことはないので $a>0$ とする。$0\le t< T$ について
$$
w(t):=\int_{\mathbb{R}^n}\Phi(x_0-y,T-t)\,u(y,t)\,dy=S(T-t)[u(t)](x_0)
$$
とおく。$u$ は $[0,t]$ で有界なので $w(t)$ は有限である。mild 解の式に $S(T-t)$ を当て、Tonelli の定理と lem-fe-semigroup の 3 を使うと
$$
w(t)=a+\int_0^tS(T-s)[u(s)^p](x_0)\,ds
$$
である。$y\mapsto\Phi(x_0-y,T-s)$ は確率密度で、$r\mapsto r^p$ は $[0,\infty)$ で凸なので、Jensenの不等式 の定理「Jensenの不等式」により $S(T-s)[u(s)^p](x_0)\ge w(s)^p$ である。よって
$$
w(t)\ge W(t):=a+\int_0^tw(s)^p\,ds\qquad(0\le t< T)
$$
である。$t'< T$ を固定すると、$w$ は $[0,t']$ で有界な可測関数なので $W$ は $[0,t']$ で Lipschitz 連続で、$W\ge a>0$、ほとんどすべての $t$ で $W'(t)=w(t)^p\ge W(t)^p$ である。$r\mapsto r^{1-p}$ は $[a,\infty)$ で Lipschitz 連続なので $G:=W^{1-p}$ も Lipschitz 連続で、ほとんどすべての $t$ で $G'=(1-p)W^{-p}W'\le-(p-1)$ である。Lipschitz 連続な関数は導関数の積分で表されるので、$0< G(t)\le a^{1-p}-(p-1)t$($t\le t'$)である。$t'< T$ は任意なので、すべての $t< T$ で $(p-1)t< a^{1-p}$、したがって $(p-1)T\le a^{1-p}$ である。$\square$
この命題から、初期値を大きくするとどの $p>1$ でも爆発することが分かる。
$g\ge0$ を有界な Borel 可測関数で、ほとんど至る所 $0$ というわけではないものとする。このとき $\lambda>\bigl((p-1)^{1/(p-1)}S(1)g(0)\bigr)^{-1}$ なら $T^*(\lambda g)\le1$ である。
$\Phi>0$ なので $S(1)g(0)>0$ である。$T^*(\lambda g)>1$ なら prop-fe-lifespan を $T=1$、$x_0=0$ で使って $(p-1)\lambda^{p-1}(S(1)g(0))^{p-1}\le1$ となり、$\lambda$ の仮定に反する。$\square$
$1< p< 1+\dfrac2n$ とし、$u_0\ge0$ を有界な Borel 可測関数で、ほとんど至る所 $0$ というわけではないものとする。このとき $T^*(u_0)<\infty$ である。すなわち、初期値がどれほど小さくても解は有限時間で爆発する。
$R>0$ を $m:=\int_{|y|\le R}u_0(y)\,dy>0$ となるようにとる。$T\ge1$ なら $|y|\le R$ で $e^{-|y|^2/(4T)}\ge e^{-R^2/4}$ なので
$$
S(T)u_0(0)\ge(4\pi T)^{-n/2}e^{-R^2/4}m=:c\,T^{-n/2}\qquad(c>0)
$$
である。$T^*=\infty$ と仮定すると、すべての $T\ge1$ で prop-fe-lifespan が使えて
$$
(p-1)\,c^{p-1}\,T^{\,1-n(p-1)/2}\le(p-1)\,T\,\bigl(S(T)u_0(0)\bigr)^{p-1}\le1
$$
である。$p< 1+2/n$ なので指数 $1-n(p-1)/2$ は正で、$T\to\infty$ とすると左辺は $\infty$ に発散し、矛盾する。$\square$
臨界の場合 $p=p_F$ には指数 $1-n(p-1)/2$ が $0$ になり、同じ議論からは初期値の質量が大きいときの爆発だけが従う。
$p=1+\dfrac2n$ とし、$u_0\ge0$ を有界な Borel 可測関数とする。$\int u_0\,dy=\infty$ であるか、または
$$
\frac2n\Bigl((4\pi)^{-n/2}\int_{\mathbb{R}^n}u_0\,dy\Bigr)^{2/n}>1
$$
なら、$T^*(u_0)<\infty$ である。
$p-1=2/n$ である。$T^{n/2}S(T)u_0(0)=(4\pi)^{-n/2}\int e^{-|y|^2/(4T)}u_0(y)\,dy$ を見ると、被積分関数は $T$ について単調に増加して $u_0$ に収束するので、Lebesgue積分 の定理「単調収束定理」により $T\to\infty$ で $(4\pi)^{-n/2}\int u_0\,dy$ に収束する。$T^*=\infty$ なら、すべての $T>0$ で prop-fe-lifespan から
$$
\frac2n\,\bigl(T^{n/2}S(T)u_0(0)\bigr)^{2/n}=\frac2n\,T\,\bigl(S(T)u_0(0)\bigr)^{2/n}\le1
$$
であり、$T\to\infty$ とすると仮定に反する。$\square$
実際には臨界の場合にも、初期値が恒等的に $0$ でない有界な非負の連続関数なら、解はすべて有限時間で爆発する(KST77 Theorem 2.1、pp. 409–410。非線形項 $f(u)=u^p$ は、$p=1+2/n$ のとき同定理の条件 (A.1)〜(A.3) を定数 $c=1$ で満たす。この記事では証明しない)。臨界の場合は藤田の論文の 2 つの定理のどちらにも含まれず、その後の研究で爆発の側に入ることが示された(KST77 §1、p. 408)。
$p>1+\dfrac2n$ とし、$q:=\dfrac{n(p-1)}2\ (>1)$、$\gamma>0$、
$$
\delta_0:=\Bigl(\frac{(q-1)(4\pi)^q\gamma^{q-1}(1-2^{1-p})}{p-1}\Bigr)^{1/(p-1)}
$$
とおく。有界な Borel 可測関数 $u_0$ が、ある $0<\delta\le\delta_0$ について $0\le u_0\le\delta\,\Phi(\cdot,\gamma)$ を満たすなら、$T^*(u_0)=\infty$ であり、解は
$$
0\le u(x,t)\le2\delta\,\Phi(x,t+\gamma)\qquad(x\in\mathbb{R}^n,\ t\ge0)
$$
を満たす。とくに $\sup_xu(x,t)\le2\delta\,(4\pi(t+\gamma))^{-n/2}$ で、解は熱方程式の解と同じ速さで減衰する。
要点:$U(t):=S(t)u_0$ とし、$V(t)=\lambda(t)U(t)$ の形の関数で $\Gamma_{u_0}[V]\le V$ となるものを作り、$0$ から始めた逐次近似を $V$ の下に閉じ込める。
上界関数:lem-fe-semigroup の 1・3 から $0\le U(t)\le\delta\Phi(\cdot,t+\gamma)\le m(t):=\delta(4\pi(t+\gamma))^{-n/2}$ である。
$$\int_0^\infty m(s)^{p-1}ds=\delta^{p-1}(4\pi)^{-q}\int_0^\infty(s+\gamma)^{-q}ds=\frac{\delta^{p-1}(4\pi)^{-q}\gamma^{1-q}}{q-1}$$
なので、$\delta\le\delta_0$ から $(p-1)\int_0^\infty m^{p-1}ds\le1-2^{1-p}$ である。よって $\lambda(t):=\bigl(1-(p-1)\int_0^tm(s)^{p-1}ds\bigr)^{-1/(p-1)}$ はすべての $t\ge0$ で定義され、$1\le\lambda\le2$、$\lambda'=\lambda^pm^{p-1}$ である。$V(x,t):=\lambda(t)U(x,t)$ とおく。
$\Gamma_{u_0}[V]\le V$:$0\le U(s)\le m(s)$ なので $V(s)^p=\lambda(s)^pU(s)^p\le\lambda(s)^pm(s)^{p-1}U(s)=\lambda'(s)U(s)$ である。lem-fe-semigroup の 1・3 から $S(t-s)[V(s)^p]\le\lambda'(s)S(t-s)U(s)=\lambda'(s)U(t)$ なので
$$\Gamma_{u_0}[V](t)\le U(t)+\int_0^t\lambda'(s)\,ds\;U(t)=\lambda(t)U(t)=V(t)$$
である。
逐次近似:$u^{(0)}:=0$、$u^{(k+1)}:=\Gamma_{u_0}[u^{(k)}]$ とおく。$u\le v$ なら $\Gamma_{u_0}[u]\le\Gamma_{u_0}[v]$ である($r\mapsto r^p$ の単調性と lem-fe-semigroup の 1)。$u^{(0)}=0\le u^{(1)}=U$ と $u^{(0)}\le V$ から始めて帰納法により、すべての $k$ で $u^{(k)}\le u^{(k+1)}\le V$ である($u^{(k+1)}=\Gamma[u^{(k)}]\le\Gamma[V]\le V$)。$u:=\lim_ku^{(k)}$ は Borel 可測で $0\le u\le V\le2\delta(4\pi\gamma)^{-n/2}$ であり、Lebesgue積分 の定理「単調収束定理」により $\Gamma_{u_0}[u^{(k)}]\to\Gamma_{u_0}[u]$ なので $u=\Gamma_{u_0}[u]$ である。よって $u$ は $[0,\infty)$ 上の mild 解で、一意性(thm-fe-local の 2)からこれが初期値 $u_0$ の解であり、$T^*=\infty$ である。最後に $V\le2\delta\Phi(\cdot,t+\gamma)$ である。$\square$
thm-fe-blowup と thm-fe-global は、それぞれ Fuj66 の Theorem 1 と Theorem 2(p. 111)を、有界な mild 解の枠組みで述べたものである。ここでの証明は原論文の議論と異なる。爆発の側は、原論文が随伴方程式の解 $\Phi(x,t-s+\varepsilon)$ で重みを付けた積分の微分不等式(Fuj66 Lemma 2.1、§2)を使うのに対し、mild 解の式と Jensen の不等式から直接導いた。大域存在の側は、原論文の重み付きノルムでの縮小写像(Fuj66 §3)の代わりに、$\lambda(t)S(t)u_0$ という比較関数を使った。いずれも、熱核の最大値が $t^{-n/2}$ で減ることと、$\int^\infty t^{-n(p-1)/2}dt$ が有限かどうかで話が分かれる。
| $p$ の範囲 | 小さい非負の初期値 | 大きい非負の初期値 | 根拠 |
|---|---|---|---|
| $1< p< 1+2/n$ | 0 でなければ必ず爆発 | 爆発 | thm-fe-blowup |
| $p=1+2/n$ | 0 でなければ必ず爆発 | 爆発 | prop-fe-critical-mass(質量が大きい場合)、KST77 |
| $p>1+2/n$ | $\delta\Phi(\cdot,\gamma)$ 以下なら大域解 | 爆発しうる | thm-fe-global、cor-fe-large |
次元が上がるほど $p_F=1+2/n$ は $1$ に近づき、大域解が存在する範囲が広がる。空間に広がる余地の大きい高い次元ほど、拡散が非線形項に勝ちやすい。$n=1,2,3$ では $p_F=3,\,2,\,5/3$ である。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| $p>1$(thm-fe-blowup) | $p=1$、$u(t)=e^tS(t)u_0$ | 有限時間の爆発 |
| $p< 1+2/n$(thm-fe-blowup) | $p>1+2/n$、$u_0=\delta\Phi(\cdot,1)$($\delta\le\delta_0$) | 有限時間の爆発 |
| $\delta\le\delta_0$(thm-fe-global) | $p>1+2/n$、$u_0=\lambda\Phi(\cdot,\gamma)$($\lambda$ が大きい) | 時間大域的な存在 |
各行は左の列の条件だけを破っている。
1 行目:$p=1$ の mild 解の式は $u(t)=S(t)u_0+\int_0^tS(t-s)u(s)\,ds$ で、$u(s)=e^sS(s)u_0$ を入れると lem-fe-semigroup の 3 から右辺は $S(t)u_0+\int_0^te^s\,ds\,S(t)u_0=e^tS(t)u_0$ である。これは各 $[0,T]$ で有界なので、解は時間大域的である。初期値の大きさや次元にはよらない。
2 行目:$p>1+2/n$ なので、thm-fe-global($\gamma=1$)により時間大域的な解がある。初期値は正で、ほとんど至る所 $0$ というわけではないので、破れている条件は $p< 1+2/n$ だけである。
3 行目:初期値は $\Phi(\cdot,\gamma)$ の定数倍という形を保ち、$p>1+2/n$ も満たすが、cor-fe-large により $\lambda$ が大きければ $T^*\le1$ である。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する