文献あり

非整数階時間微分を含む拡散方程式の初期値境界値問題

time-fractional-diffusion,PDE,weak solution

234

Initial-boundary value problems for time-fractional diffusion equations with time-independent coefficients

はじめに

近年, 異常拡散などの非局所拡散現象を記述する非整数階時間微分を含む拡散方程式の研究が盛んに行われている. 非整数階時間微分を含む拡散方程式とは
$(P)_{0}^{c} D_{t}^{α} u = L u$
である. ここで, $u = u (x, t)$ は未知関数, $L$ は一様楕円型作用素であり, $_{0}^{c} D_{t}^{α} u$ は $α$ 階Caputo微分と呼ばれ, $α \in (0, 1)$ と十分滑らかな $u$ に対して,
$_{0}^{c} D_{t}^{α} u := \int_{0}^{t} g_{α} (t - τ) u^{'} (τ) d τ, g_{α} (t) = \frac{t^{- α}}{Γ (1 - α)}$
と定義される. 方程式 $(P)$ は, 作用素 $L$ が空間変数にのみ依存する場合の解析は多く行われている. 実際, 有界領域における初期値境界値問題はLuchko4によって, Fourier級数を用いた解の表現が与えられている. また, 全空間における初期値問題は基本解の導出, 漸近挙動の解析がEidelman-Kochubei1によって与えられている. 一方, 作用素 $L$ が時間変数にも依存する場合は, $L^{2}$ 空間といったHilbert空間に属する解に対する評価しか行われておらず, 古典解の存在性などの基礎理論ですら構築されていない. その理由としては, 半群理論が適用できないことや, $t = 0$ での特異性のため, 常微分方程式の基礎理論が適用できないことが挙げられる.

上記の背景を基に, 今回は次の初期値境界値問題
$\begin{matrix} (1) & {\begin{cases} _{0}^{c} D_{t}^{α} u = L u + f & i n Ω \times (0, T) =: Q_{T}, \\ u = 0 & o n \partial Ω \times [0, T] =: Σ_{T}, \\ u = u_{0} & o n Ω \times {t = 0} \end{cases} \end{matrix}$
を考える. ただし, $L$ は
$L u = \sum_{i, j = 1}^{N} \partial_{i} (a_{i, j} (x) \partial_{j} u) + \sum_{j = 1}^{N} b_{j} (x) \partial_{j} u + c (x) u, \partial_{i} := \frac{\partial}{\partial x_{i}}$
であり, 任意の $x \in \overset{―}{Ω}$ に対して,
$\begin{matrix} (2) & λ | ξ |^{2} ⩽ \sum_{i, j = 1}^{N} a_{i, j} (x) ξ_{i} ξ_{j} ⩽ μ | ξ |^{2} \end{matrix}$
なる $λ, μ > 0$ が存在し, $a_{i, j} = a_{j, i}$ をみたす. このとき, 任意の $x \in \overset{―}{Ω}$ に対して $| a_{i, j} | ⩽ μ$ である. また, $b = (b_{1}, \dots, b_{N})$ , $A C [0, T]$ を $[0, T]$ 上で定義された絶対連続関数の空間とする. $a_{i, j}, b, c$ の正則性は後に与える. 以下の議論では, 3を基に行う(3では $L$ が時間変数にも依存している場合であるが, それについては次回述べる).

次に, 問題 $(1)$ の弱解の定義を次のように与える.

Weak solution

$u$ が次の $(i), (i i)$ をみたすとき, 問題 $(1)$ の弱解という.
(i) $u \in W (u_{0}, H_{0}^{1} (Ω), L^{2} (Ω) := {u \in L^{2} (0, T; H_{0}^{1} (Ω)); I^{1 - α} [u - u_{0}] \in_{0} H^{1} (0, T; H^{- 1} (Ω))}$ ,
(ii) 任意の $φ \in H_{0}^{1} (Ω)$ に対して,
$\begin{matrix} \frac{d}{d t} \int_{Ω} I^{1 - α} [u (x, \cdot) - u_{0} (x)] (t) φ (x) d x + \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} u (x, t) \partial_{i} φ (x) d x \\ (3) & = \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} u (x, t) φ (x) d x + \int_{Ω} c (x) u (x, t) φ (x) d x + \int_{Ω} f (x, t) φ (x) d x, a . e . t \in (0, T) \end{matrix}$
をみたす. ただし,
$_{0} H^{1} (0, T) = {f \in H^{1} (0, T); f (0) = 0}$
である.

問題 $(1)$ に対して, 次の定理が得られる.

Exitence of weak solution and energy estimate

$u_{0} \in L^{2} (Ω)$ , $f \in L^{2} (0, T; H^{- 1} (Ω))$ , 式 $(2)$ が成立し, $b, c \in L^{\infty} (Ω)$ とする. このとき, 問題 $(1)$ の弱解 $u = u (x, t)$ が一意に存在する. さらに,
$\begin{matrix} (4) & {‖ I^{1 - α} [u - u_{0}] ‖}_{H^{1} (0, T; H^{- 1} (Ω))} + ∥ u ∥_{L^{2} (0, T; H_{0}^{1} (Ω))} + ∥ u ∥_{H^{\frac{α}{2}} (0, T; L^{2} (Ω))} ⩽ C (∥ u_{0} ∥_{L^{2} (Ω)} + ∥ f ∥_{L^{2} (0, T; H^{- 1} (Ω))}) \end{matrix}$
が成立する. ここで, $C$ は $α, μ, λ, T, N, ∥ b ∥_{L^{\infty} (Ω)}, ∥ c ∥_{L^{\infty} (Ω)}$ に依存する定数である. また, $α \in (\frac{1}{2}, 1)$ のとき, $u \in C ([0, T]; H^{- 1} (Ω))$ である.

Approximate solutions

本章では, 通常の放物型方程式に対するGalerkin methodと同様に, 式 $(1)$ の近似解を構成することを目的とする. ${φ_{n} (x)}_{n \in N}$ を
${\begin{cases} - Δ φ_{n} = λ_{n} φ_{n} & i n Ω, \\ φ_{n} = 0 & o n \partial Ω \end{cases}$
をみたす $L^{2} (Ω)$ の正規直交基底, $H_{0}^{1} (Ω)$ の直交基底とし, 次の問題 $(1)$ の近似解
$\begin{matrix} (5) & u_{n} (x, t) = \sum_{k = 1}^{n} d_{n, k} (t) φ_{k} (x) \end{matrix}$
を求める. すなわち, 係数 $d_{n, k}$ を決定する. そのため, 次の初期値境界値問題
$\begin{matrix} (6) & {\begin{cases} _{0}^{c} D_{t}^{α} u_{n} = L u_{n} + f_{n} & i n Q_{T}, \\ u_{n} = 0 & o n Σ_{T}, \\ u_{n} = u_{n, 0} & o n Ω \times {t = 0} \end{cases} \end{matrix}$
を考える. ただし,
$u_{n, 0} (x) = φ_{k} (x) (\sum_{k = 1}^{n} \int_{Ω} u_{0} (y) φ_{k} (y) d y),$
$f_{n} (x, t) = (η_{1 / n} * f (x, \cdot)) (t), f_{n} \to f s t r o n g l y i n L^{2} (Q_{T}) a s n \to \infty$
である. 係数 $d_{n, k}$ は式 $(6)$ を ${φ_{1}, \dots, φ_{n}}$ による有限次元空間への射影を考えることによって決定する. すなわち, 問題 $(6)$ の $1$ 行目の方程式の両辺に $φ_{m}$ をかけて $Ω$ 上で積分をすると, 左辺は
$\begin{aligned} \int_{Ω} φ_{m} (x)_{0}^{c} D_{t}^{α} u_{n} (x, t) d x & = \int_{Ω} φ_{m} (x)_{0}^{c} D_{t}^{α} (\sum_{k = 1}^{n} d_{n, k} (t) φ_{k} (x)) d x \\ =_{0}^{c} D_{t}^{α} (\sum_{k = 1}^{n} d_{n, k} (t) \int_{Ω} φ_{m} (x) φ_{k} (x) d x) \\ =_{0}^{c} D_{t}^{α} d_{n, m} (t) \end{aligned}$
となり, 右辺第1項は部分積分より,
$\begin{aligned} \sum_{i, j = 1}^{N} \int_{Ω} φ_{m} (x) \partial_{i} (a_{i, j} (x) \partial_{j} u_{n} (x, t)) d x & = - \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} φ_{k} (x) \sum_{k = 1}^{n} d_{n, k} (t) \partial_{i} φ_{m} (x) d x \\ = - \sum_{i, j = 1}^{N} \sum_{k = 1}^{n} d_{n, k} (t) \int_{Ω} a_{i, j} (x) \partial_{j} φ_{k} (x) \partial_{i} φ_{m} (x) d x \end{aligned}$
となるので,
$\begin{matrix} _{0}^{c} D_{t}^{α} d_{n, m} (t) = - \sum_{k = 1}^{n} \sum_{i, j = 1}^{N} d_{n, k} (t) \int_{Ω} a_{i, j} (x) \partial_{j} φ_{k} (x) \partial_{i} φ_{m} (x) d x \\ + \sum_{k = 1}^{n} \sum_{j = 1}^{N} d_{n, k} (t) \int_{Ω} b_{j} (x) \partial_{j} φ_{k} (x) φ_{m} (x) d x + \sum_{k = 1}^{n} d_{n, k} (t) \int_{Ω} c (x) φ_{k} (x) φ_{m} (x) d x \\ (7) & + \int_{Ω} f_{n} (x, t) φ_{m} (x) d x \end{matrix}$
が得られる.
$d_{n} (t) = (d_{n, 1} (t), \dots, d_{n, n} (t)),$
$A_{m, k}^{n} = \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} ϕ_{k} (x) \partial_{i} ϕ_{m} (x) d x, A^{n} (t) = {A_{m, k}^{n}}_{k, m = 1}^{n},$
$B_{m, k}^{n} = \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} ϕ_{k} (x) ϕ_{m} (x) d x, B^{n} = {B_{m, k}^{n}}_{k, m = 1}^{n},$
$C_{m, k}^{n} = \int_{Ω} c_{n} (x) ϕ_{k} (x) ϕ_{m} (x) d x, C^{n} = {C_{m, k}^{n}}_{k, m = 1}^{n},$
$F^{n} (t) = (\int_{Ω} f_{n} (y, t) ϕ_{1} (y) d y, \dots, \int_{Ω} f_{n} (y, t) ϕ_{n} (y) d y),$
$d_{n, 0} = (\int_{Ω} u_{0} (y) ϕ_{1} (y) d y, \dots, \int_{Ω} u_{0} (y) ϕ_{n} (y) d y)$
と定義すると, 式 $(6)$ は,
$\begin{matrix} (8) & {\begin{cases} _{0}^{c} D_{t}^{α} d_{n} (t) = - \tilde{A} d_{n} (t) + F^{n} (t), \\ d_{n} (0) = d_{n, 0} \end{cases} \end{matrix}$
と表される. ここで, $- \tilde{A} = A^{n} - B^{n} - C^{n}$ である. 以上の議論から次の補題が成立する.

任意の $n \in N$ とに対して, 式 $(8)$ をみたす $d_{n} \in C^{1} ((0, T]) \cap C ([0, T])$ が一意に存在する. さらに, $t^{1 - α} d_{n}^{'} \in C [0, T]$ である.

方程式の両辺をLaplace変換すると, $α$ 階Caputo微分のLaplace変換は
$L [_{0}^{c} D_{t}^{α} d_{n}] (s) = s^{α} L [d_{n}] (s) - d_{n, 0} s^{α - 1}$
となるので, $L [d_{n}]$ について整理すると
$L [d_{n}] (s) = d_{n, 0} \frac{s^{α - 1}}{s^{α} + \tilde{A}} + L [F^{n}] (s) \frac{1}{s^{α} + \tilde{A}}$
となる. これは
$L [d_{n}] (s) = d_{n, 0} L [E_{α, 1} (- \tilde{A} t^{α})] (s) + L [F^{n}] (s) L [t^{α - 1} E_{α, α} (- \tilde{A} t^{α})] (s)$
を示唆している. したがって, 両辺をLaplace逆変換すると,
$\begin{matrix} (9) & d_{n} (t) = d_{n, 0} E_{α, 1} (- \tilde{A} t^{α}) + \int_{0}^{t} (t - s)^{α - 1} E_{α, α} (- \tilde{A} (t - s)^{α}) F^{n} (s) d s \end{matrix}$
となる. これが方程式の解になることは直接代入して確かめればよい. また, 一意性に関しては https://mathlog.info/articles/r2II7keeJFnfPUFwwp5C よりしたがう. 式 $(9)$ の両辺を微分すると,
$d_{n}^{'} (t) = - A t^{α - 1} E_{α, α} (- A t^{α}) + \int_{0}^{t} (t - s)^{α - 1} E_{α, α} (- \tilde{A} (t - s)^{α}) (F^{n})^{'} (s) d s + t^{α - 1} E_{α, α} (- \tilde{A} t^{α}) F^{n} (0)$
となる. よって, 両辺に $t^{1 - α}$ をかければ右辺は $[0, T]$ 上で連続関数となる. 以上で補題の証明が得られた.

この補題より, 次がしたがう.

$n \in N$ とする. このとき, 式 $(5)$ と式 $(8)$ で与えられる $u_{n}$ は, $m \in {1, \dots, n}$ に対して
$\begin{matrix} \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) φ_{m} (x) d x + \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} u_{n} (x, t) \partial_{i} φ_{m} (x) d x \\ (10) & = \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} u_{n} (x, t) φ_{m} (x) d x + \int_{Ω} c (x) u_{n} (x, t) φ_{m} (x) d x + \int_{Ω} f_{n} (x, t) φ_{m} (x) d x \end{matrix}$
をみたす. さらに, $x \in Ω$ と $β \in N^{N}$ に対して, $\partial Ω$ が十分滑らかなとき $\partial_{x}^{β} u_{n} (x, \cdot) \in A C [0, T]$ かつ $t^{1 - α} \partial_{x}^{β} u_{n}^{'} (x, \cdot) \in C ({\overset{―}{Q}}_{T})$ である.

Weak solutions

本章では, 問題 $(1)$ に対するEnergy評価を与える. そのため, 次の補題を証明する.

$w \in L^{2} (Q_{T})$ ,
$\begin{matrix} (11) & w (x, \cdot) \in A C [0, T] f o r x \in Ω, \end{matrix}$
かつ
$\begin{matrix} (12) & t^{1 - α} w_{t} \in L^{\infty} (Q_{T}) \end{matrix}$
と仮定する. このとき, 等式
$\begin{matrix} _{0}^{c} D_{t}^{α} ∥ w (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \frac{∥ w (t) - w (τ) ∥_{L^{2} (Ω)}^{2}}{(t - τ)^{α + 1}} d τ \\ (13) & + g_{α} (t) ∥ w (t) - w (0) ∥_{L^{2} (Ω)}^{2} = 2 \int_{Ω}_{0}^{c} D_{t}^{α} w (x, t) \cdot w (x, t) d x \end{matrix}$
が成り立つ.

形式的に計算する. Caputo微分の定義より,
$\begin{aligned} 2 \int_{Ω}_{0}^{c} D_{t}^{α} w (x, t) w (x, t) d x -_{0}^{c} D_{t}^{α} ∥ w (t) ∥_{L^{2} (Ω)}^{2} & = 2 \int_{Ω} w (x, t) \int_{0}^{t} g_{α} (t - s) w_{s} (x, s) d s d x - 2 \int_{0}^{t} g_{α} (t - s) \int_{Ω} w_{s} (x, s) w (x, s) d s \\ = 2 \int_{0}^{t} g_{α} (t - s) \int_{Ω} w_{s} (x, s) (w (x, t) - w (x, s)) d x d s \\ = - \int_{0}^{t} g_{α} (t - s) \frac{\partial}{\partial s} \int_{Ω} | w (x, t) - w (x, s) |^{2} d x d s \\ = - [g_{α} (t - s) \int_{Ω} | w (x, t) - w (x, s) |^{2} d x]_{s = 0}^{s = t} + \frac{α}{Γ (1 - α)} \int_{0}^{t} (t - s)^{- a - 1} \int_{Ω} | w (x, t) - w (x, s) |^{2} d x d s \\ = g_{α} (t) ∥ w (t) - w (0) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \frac{∥ w (x, t) - w (x, s) ∥_{L^{2} (Ω)}^{2}}{(t - s)^{a + 1}} d s \end{aligned}$
となる. 最後の等式第1項では $w$ に対して平均値の定理を用いることで得られる.

補題3より近似解に対してEnergy評価を与えることができる.

Energy Estimate for approximate solution

$u_{0} \in L^{2} (Ω)$ , $f \in L^{2} (0, T; H^{- 1} (Ω))$ , ${a_{i, j}}_{i, j = 1}^{N}$ が仮定 $(2)$ をみたし, $b, c \in L^{\infty} (Ω)$ とする. このとき, 各 $t \in [0, T]$ と $n \in N$ に対して, 近似解 $u_{n}$ は
$\begin{matrix} I^{1 - α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \int_{0}^{s} \frac{∥ u_{n} (s) - u_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d τ d s \\ + \int_{0}^{t} g_{α} (s) ∥ u_{n} (s) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} d s + λ \int_{0}^{t} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s \\ (14) & ⩽ C_{0} (∥ u_{0} ∥_{L^{2} (Ω)}^{2} + \int_{0}^{t} ∥ f_{n} (s) ∥_{H^{- 1}}^{2} d s) \end{matrix}$
をみたす. ここで, $C_{0}$ は $∥ b ∥_{L^{\infty} (Ω)}, ∥ c ∥_{L^{\infty} (Ω)}, λ, α, T$ に依存する定数である.

式 $(10)$ の両辺に $d_{n}, m$ をかけて $m = 1$ から $n$ まで和をとると,
$\begin{matrix} \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) u_{n} (x, t) d x + \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} u_{n} (x, t) \partial_{i} u_{n} (x, t) d x \\ = \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} u_{n} (x, t) u_{n} (x, t) d x + \int_{Ω} c_{n} (x) | u_{n} (x, t) |^{2} d x + \int_{Ω} f_{n} (x, t) u_{n} (x, t) d x \end{matrix}$
となる. 左辺は補題3と仮定 $(2)$ , 右辺にはHolderの不等式とCauchyの不等式より
$\begin{aligned} \frac{1}{2}_{0}^{c} D_{t}^{α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{2 Γ (1 - α)} \int_{0}^{t} \frac{∥ u_{n} (t) - u_{n} (s) ∥_{L^{2} (Ω)}^{2}}{(t - s)^{α + 1}} d s + \frac{1}{2} g_{α} (t) ∥ u_{n} (t) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} + λ ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} \\ ⩽ ∥ b ∥_{L^{\infty} (Ω)} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)} ∥ u_{n} (t) ∥_{L^{2} (Ω)} + ∥ c ∥_{L^{\infty} (Ω)} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ f_{n} (t) ∥_{H^{- 1} (Ω)} {(∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2})}^{\frac{1}{2}} \\ ⩽ \frac{1}{λ} ∥ b ∥_{L^{\infty} (Ω)}^{2} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{λ}{4} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ c ∥_{L^{\infty} (Ω)} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{1}{λ} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2} + \frac{λ}{4} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{λ}{4} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} \end{aligned}$
と評価できる. よって, まとめると
$\begin{matrix} _{0}^{c} D_{t}^{α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \frac{∥ u_{n} (t) - u_{n} (s) ∥_{L^{2} (Ω)}^{2}}{(t - s)^{α + 1}} d s + g_{α} (t) ∥ u_{n} (t) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} + λ ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} \\ (15) & ⩽ C ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{1}{λ} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2} \end{matrix}$
が得られる. ここで, $C$ は $∥ b ∥_{L^{\infty} (Ω)}, ∥ c ∥_{L^{\infty} (Ω)}, λ$ に依存する定数である. 次に $∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2}$ に対して評価を行う. 式 $(15)$ より
$\begin{matrix} (16) & _{0}^{c} D_{t}^{α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ C ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{1}{λ} ∥ f (t) ∥_{H^{- 1} (Ω)}^{2} \end{matrix}$
であるから,　式 $(16)$ の両辺を $α$ 階Riemann-Liouville積分すると, $∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} \in A C [0, T]$ であるので
$∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ ∥ u_{0} ∥_{L^{2} (Ω)}^{2} + C I^{α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{1}{λ} I^{α} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2}$
となる. したがって, Gronwallの不等式より
$∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ ∥ u_{0} ∥_{L^{2} (Ω)}^{2} \sum_{k = 0}^{\infty} \frac{t^{α k}}{Γ (α k + 1)} + \frac{1}{λ} \sum_{k = 0}^{\infty} C^{k} I^{α (k + 1)} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2}$
となる. この右辺の級数が絶対収束することを確かめる.
$\sum_{k = 0}^{\infty} \frac{t^{α k}}{Γ (α k + 1)} = E_{α, 1} (t^{α}) < \infty,$
$\begin{aligned} \sum_{k = 0}^{\infty} C^{k} I^{α (k + 1)} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2} & = \sum_{k = 0}^{\infty} C^{k} \int_{0}^{t} \frac{(t - s)^{α (k + 1) - 1}}{Γ (α (k + 1))} ∥ f_{n} (s) ∥_{H^{- 1} (Ω)}^{2} d s \\ ⩽ sup_{s \in [0, t]} ∥ f_{n} (s) ∥_{H^{- 1} (Ω)}^{2} \sum_{k = 0}^{\infty} C^{k} \frac{t^{α (k + 1)}}{Γ (α k + α + 1)} \\ = sup_{s \in [0, t]} ∥ f_{n} (s) ∥_{H^{- 1} (Ω)}^{2} t^{α} E_{α, α + 1} (C t^{α}) < \infty \end{aligned}$
と評価できる. よって, 式 $(16)$ の両辺を $0$ から $t$ まで積分し, Riemann-Liouville積分の加法性 $I = I^{1 - α} I^{α}$ を用いると
$\begin{matrix} I^{1 - α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \int_{0}^{s} \frac{∥ u_{n} (s) - u_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d τ d s \\ + \int_{0}^{t} g_{α} (t) ∥ u_{n} (s) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} d s + λ \int_{0}^{t} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s \\ ⩽ \frac{t^{1 - α}}{Γ (2 - α)} ∥ u_{0} ∥_{L^{2} (Ω)}^{2} + \frac{1}{λ} \int_{0}^{t} ∥ f_{n} (s) ∥_{L^{2} (Ω)}^{2} d s + C I ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} \end{matrix}$
となる. ここで, 先程の議論より
$\begin{aligned} C I ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} & ⩽ C ∥ u_{0} ∥_{L^{2} (Ω)}^{2} \sum_{k = 0}^{\infty} \frac{t^{α k + 1}}{Γ (α k + 2)} + \frac{1}{λ} \sum_{k = 0}^{\infty} C^{k + 1} I^{α (k + 1) + 1} ∥ f_{n} (t) ∥_{H^{- 1 (Ω)}}^{2} \\ = C ∥ u_{0} ∥_{L^{2} (Ω)}^{2} \sum_{k = 0}^{\infty} \frac{t^{α k + 1}}{Γ (α k + 2)} + \frac{1}{λ} \sum_{k = 0}^{\infty} \int_{0}^{t} \frac{C^{k + 1} (t - s)^{α (k + 1)}}{Γ (α (k + 1) + 1)} ∥ f_{n} (s) ∥_{H^{- 1} (Ω)}^{2} d s \\ ⩽ C ∥ u_{0} ∥_{L^{2} (Ω)}^{2} \sum_{k = 0}^{\infty} \frac{t^{α k + 1}}{Γ (α k + 2)} + \frac{1}{λ} \sum_{k = 0}^{\infty} \frac{C^{k + 1} t^{α (k + 1)}}{Γ (α (k + 1) + 1)} \int_{0}^{t} ∥ f_{n} (s) ∥_{H^{- 1} (Ω)}^{2} d s \\ ⩽ C ∥ u_{0} ∥_{L^{2} (Ω)}^{2} t E_{α, 2} (t^{α}) + \frac{C}{λ} t^{α} E_{α, α + 1} (C t^{α}) \int_{0}^{t} ∥ f_{n} (s) ∥_{H^{- 1} (Ω)}^{2} d s \end{aligned}$
と評価できる. さらに,
$sup_{n} sup_{t \in [0, T]} t E_{α, 2} (t^{α}) = sup_{n} T E_{α, 2} (T^{α}) =: D_{1},$
$sup_{n} sup_{t \in [0, T]} t^{α} E_{α, α + 1} (C t^{α}) = sup_{n} T^{α} E_{α, α + 1} (C T^{α}) =: D_{2}$
とする. ただし, $D_{1}, D_{2}$ は $∥ b ∥_{L^{\infty} (Ω)}, ∥ c ∥_{L^{\infty} (Ω)}, λ, T$ に依存する定数である. 以上で式 $(14)$ が得られ, 補題の証明が完了した.

定理1の証明

補題3より
$\begin{matrix} (17) & ∥ \nabla u_{n} ∥_{L^{2} (0, T; L^{2} (Ω))} + ∥ u_{n} ∥_{H^{\frac{α}{2}} (0, T; L^{2} (Ω))} ⩽ {\overset{―}{C}}_{0} \end{matrix}$
である. 次に $_{0}^{c} D_{t}^{α} u_{n}$ に対する評価を与える. $w \in H_{0}^{1} (Ω)$ のとき, $θ_{m}$ を定数として $w (x) = \sum_{m = 1}^{\infty} θ_{m} φ_{m} (x)$ である. $w_{n} (x) = \sum_{m = 1}^{n} θ_{m} φ_{m} (x)$ と表し, 式 $(10)$ の両辺に $θ_{m}$ をかけて $m = 1$ から $n$ まで和をとると
$\begin{matrix} \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) w_{n} (x) d x + \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} u_{n} (x, t) \partial_{i} w_{n} (x) d x \\ = \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} u_{n} (x, t) w_{n} (x) d x + \int_{Ω} c_{n} (x) u_{n} (x, t) w_{n} (x) d x + \int_{Ω} f_{n} (x, t) w_{n} (x) d x \end{matrix}$
となるので, Holderの不等式より
$\begin{matrix} | \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) w_{n} (x) d x | ⩽ μ N^{2} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)} ∥ \nabla w_{n} (x) ∥_{L^{2} (Ω)} + ∥ b ∥_{L^{\infty} (Ω)} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)} ∥ w_{n} ∥_{L^{2} (Ω)} \\ (18) & + ∥ c ∥_{L^{\infty} (Ω)} ∥ u_{n} (t) ∥_{L^{2} (Ω)} ∥ w_{n} ∥_{L^{2} (Ω)} + ∥ f_{n} (t) ∥_{H^{- 1} (Ω)} ∥ w_{n} ∥_{H_{0}^{1} (Ω)} \end{matrix}$
と評価できる. $u_{n} (x, \cdot) \in A C [0, T]$ より $_{0}^{c} D_{t}^{α} u_{n} = \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}]$ ,
${‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] (t) ‖}_{H^{- 1} (Ω)} = sup_{∥ w ∥_{H_{0}^{1} (Ω)} = 1} | \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) w_{n} (x) d x |$
である. よって, 式 $(18)$ の両辺は $L^{2} (0, T)$ の意味で $n$ について一様有界であるので
$\begin{matrix} (19) & sup_{n} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] ‖}_{L^{2} (0, T; H^{- 1} (Ω))} < \infty \end{matrix}$
が得られる. したがって, ${I^{1 - α} [u_{n} - u_{n, 0}]}$ は $_{0} H^{1} (0, T; H^{- 1} (Ω))$ の意味で $n$ について一様有界である. 式 $(17)$ , 式 $(19)$ より, 弱コンパクト性定理を用いると
$\begin{matrix} (20) & u_{n_{k}} \to u w e a k l y i n L^{2} (0, T; H_{0}^{1} (Ω)) a s k \to \infty, \end{matrix}$
$\begin{matrix} (21) & I^{1 - α} [u_{n_{k}} - u_{n_{k}, 0}] \to v w e a k l y i n L^{2} (0, T; H^{- 1} (Ω)) a s k \to \infty \end{matrix}$
となる部分列 ${u_{n_{k}}}$ と $u \in L^{2} (0, T; H_{0}^{1} (Ω)) \cap H^{\frac{α}{2}} (0, T; L^{2} (Ω))$ と $v \in_{0} H^{1} (0, T; H^{- 1} (Ω))$ が存在する. まず弱微分 $\frac{d}{d t} I^{1 - α} [u - u_{0}]$ が $L^{2} (0, T; H^{- 1} (Ω))$ の意味で存在することを示す. $ϕ \in C_{0}^{\infty} (0, T)$ , $φ \in H_{0}^{1} (Ω)$ とすると,
$\begin{aligned} \int_{0}^{T} ϕ (t) \int_{Ω} \frac{d}{d t} v (x, t) φ (x) d x d t & = lim_{k \to \infty} \int_{0}^{T} ϕ (t) \int_{Ω} \frac{d}{d t} I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) φ (x) d x d t \\ = lim_{k \to \infty} \int_{Ω} φ (x) \int_{0}^{T} ϕ (t) \frac{d}{d t} I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) d t d x \\ = - lim_{k \to \infty} \int_{Ω} φ (x) \int_{0}^{T} ϕ^{'} (t) I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) d t d x \\ = - lim_{k \to \infty} \int_{0}^{T} ϕ^{'} (t) \int_{Ω} I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) φ (x) d x d t \\ = - \int_{0}^{T} ϕ^{'} (t) \int_{Ω} I^{1 - α} [u (x, \cdot) - u (x)] (t) φ (x) d x d t \end{aligned}$
となる. よって,
$\frac{d}{d t} I^{1 - α} [u - u_{0}] = \frac{d}{d t} v \in L^{2} (0, T; H^{- 1} (Ω)) i n t h e w e a k s e n s e$
が示され, 得られた $u$ に対して式 $(4)$ が成立する. 次に, 得られた $u$ が等式 $(3)$ をみたしていることを確かめる. 稠密性より, $w (x) = \sum_{m = 1}^{K} θ_{m} φ_{m}$ に対して成立することを示せば十分である. 式 $(10)$ の両辺に $θ_{m}$ をかけて $m = K$ まで和をとる. このとき, $t_{0} \in (0, T)$ を固定し, mollifier $η_{ε} (t - t_{0})$ を両辺にかけて $0$ から $T$ まで積分すると,
$\begin{matrix} \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} \frac{d}{d t} I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) w (x) d x d t + \sum_{i, j = 1}^{N} \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} a_{i, j} (x) \partial_{j} u_{n_{k}} (x, t) \partial_{i} w (x) d x d t \\ = \sum_{j = 1}^{N} \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} b_{j} (x) \partial_{j} u_{n_{k}} (x, t) w (x) d x d t + \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} c (x) u_{n_{k}} (x, t) w (x) d x d t \\ + \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} f_{n} (x, t) w (x) d x d t \end{matrix}$
となる. それぞれの項において $k \to \infty, ε \to 0$ の極限を考える. 左辺第1項は
$\begin{aligned} \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} \frac{d}{d t} I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) w (x) d x d t & = - \int_{0}^{T} η_{ε}^{'} (t - t_{0}) \int_{Ω} I^{1 - α} [u_{n_{k}} (x, \cdot) - u_{n_{k}, 0} (x)] (t) w (x) d x d t \\ \to_{k \to \infty}^{} - \int_{0}^{T} η_{ε}^{'} (t - t_{0}) \int_{Ω} I^{1 - α} [u (x, \cdot) - u_{0} (x)] (t) w (x) d x d t \\ = \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} \frac{d}{d t} I^{1 - α} [u (x, \cdot) - u_{0} (x)] (t) w (x) d x d t \\ \to_{ε \to 0}^{} \int_{Ω} \frac{d}{d t} I^{1 - α} [u (x, \cdot) - u_{0} (x)] (t_{0}) w (x) d x f o r a . e . t_{0} \in (0, T) \end{aligned}$
となる. 以下, 同様に計算にする. $\partial_{i} w (x) η_{ε} (t - t_{0})$ は $Q_{T}$ 上で滑らかなので,
$\begin{aligned} \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} a_{i, j} (x) \partial_{j} u_{n_{k}} (x, t) \partial_{i} w (x) d x d t & \to_{k \to \infty}^{} \int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} a_{i, j} (x) \partial_{j} u (x, t) \partial_{i} w (x) d x d t \\ \to_{ε \to 0}^{} \int_{Ω} a_{i, j} (x) \partial_{j} u (x, t_{0}) \partial_{i} w (x) d x f o r a . e . t_{0} \in (0, T) \end{aligned}$
となる. 右辺は各項,
$\int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} b_{j} (x) \partial_{j} u_{n_{k}} (x, t) w (x) d x d t \to_{k \to \infty, ε \to 0}^{} \int_{Ω} b_{j} (x) \partial_{j} u (x, t_{0}) w (x) d x f o r a . e . t_{0} \in (0, T),$
$\int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} c (x) u_{n_{k}} (x, t) w (x) d x d t \to_{k \to \infty, ε \to 0}^{} \int_{Ω} c (x) u (x, t_{0}) w (x) d x f o r a . e . t_{0} \in (0, T),$
$\int_{0}^{T} η_{ε} (t - t_{0}) \int_{Ω} f_{n_{k}} (x, t) w (x) d x d t \to_{k \to \infty, ε \to 0}^{} \int_{Ω} f (x, t_{0}) w (x) d x f o r a . e . t_{0} \in (0, T)$
が得られる. したがって, 得られた $u$ が等式 $(3)$ をみたすことが確かめられた. 一意性は, $\tilde{u} = u_{1} - u_{2}$ とおき, 初期値境界値問題
${\begin{cases} _{0}^{c} D_{t}^{α} \tilde{u} = L \tilde{u} & i n Q_{T}, \\ \tilde{u} = 0 & o n Σ_{T}, \\ {\tilde{u}}_{0} = 0 & o n Ω \times {t = 0} \end{cases}$
を考え, 同様の議論をすれば $u_{1} = u_{2}$ a.e. が得られる. 最後に, $u$ の $H^{- 1} (Ω)$ 連続性を示す.
$u_{n} \to u s t r o n g l y i n L^{2} (0, T; H_{0}^{1} (Ω)) a s n \to \infty,$
$I^{1 - α} [u_{n} - u_{n, 0}] \to I^{1 - α} [u - u_{0}] s t r o n g l y i n L^{2} (0, T; H^{- 1} (Ω)) a s n \to \infty$
なる滑らかな近似列 ${u_{n}}$ をとると, $I^{1 - α} [u_{n} - u_{n, 0}] \in A C [0, T]$ であるので
$I^{α} \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] (t) = u_{n} (t) - u_{n, 0}$
となる. よって, Holderの不等式より
$\begin{aligned} ∥ u_{n} (t) - u_{n, 0} ∥_{H^{- 1} (Ω)} & = {‖ I^{α} \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] (t) ‖}_{H^{- 1} (Ω)} \\ ⩽ I^{α} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] (t) ‖}_{H^{- 1} (Ω)} \\ = \frac{1}{Γ (α)} \int_{0}^{t} (t - s)^{α - 1} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] (s) ‖}_{H^{- 1} (Ω)} d s \\ ⩽ \frac{1}{Γ (α)} {(\int_{0}^{t} (t - s)^{2 α - 2} d s)}^{\frac{1}{2}} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] ‖}_{L^{2} (0, T; H^{- 1} (Ω))} \\ = C_{α} t^{α - \frac{1}{2}} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] ‖}_{L^{2} (0, T; H^{- 1} (Ω))} \end{aligned}$
と評価できる. よって,
$∥ u_{n} - u_{n, 0} ∥_{C ([0, T]; H^{- 1} (Ω))} = max_{t \in [0, T]} ∥ u_{n} (t) - u_{n, 0} ∥_{H^{- 1} (Ω)} ⩽ C_{α, T} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] ‖}_{L^{2} (0, T; H^{- 1} (Ω))}$
が得られる. これより,
$\begin{aligned} ∥ (u_{n} - u_{n, 0}) - (u_{m} - u_{m, 0}) ∥_{C ([0, T]; H^{- 1} (Ω))} & ⩽ C_{α, T} {‖ \frac{d}{d t} I^{1 - α} [u_{n} - u_{n, 0}] - \frac{d}{d t} I^{1 - α} [u_{m} - u_{m, 0}] ‖}_{L^{2} (0, T; H^{- 1} (Ω))} \\ \to 0 a s n, m \to \infty \end{aligned}$
となるので, ${u_{n} - u_{n, 0}}$ は $C ([0, T]; H^{- 1} (Ω))$ 上のCauchy列となる. 故に, $u (0)$ はwell-definedになる. したがって, 先程の評価より
$max_{s \in [0, t]} ∥ u (t) - u_{0} ∥_{H^{- 1} (Ω)} ⩽ C_{α} t^{α - \frac{1}{2}} {‖ \frac{d}{d t} I^{1 - α} [u - u_{0}] ‖}_{L^{2} (0, T; H^{- 1} (Ω))} \to 0 a s t \to 0$
となるので, $u (0) = u_{0}$ である. 以上で定理1の証明が完了した.

Regular solutions

本章では, 次の定理を証明することを目的とする.

$u_{0} \in H_{0}^{1} (Ω)$ , $f \in L^{2} (0, T; L^{2} (Ω))$ , $b, c \in L^{\infty} (Ω)$ とし, $max_{i, j} ∥ \nabla a_{i, j} ∥_{L^{\infty} (Ω)} < \infty$ と仮定する. このとき, $I^{1 - α} [u - u_{0}] \in_{0} H^{1} (0, T; L^{2} (Ω))$ となる問題 $(1)$ の弱解 $u \in L^{2} (0, T; H^{2} (Ω)) \cap H^{\frac{α}{2}} (0, T; H_{0}^{1} (Ω))$ が一意に存在する. さらに,
$\begin{matrix} (22) & {‖ I^{1 - α} [u - u_{0}] ‖}_{H^{1} (0, T; L^{2} (Ω)} + ∥ u ∥_{L^{2} (0, T; H^{2} (Ω))} + ∥ u ∥_{H^{\frac{α}{2}} (0, T; H_{0}^{1} (Ω))} ⩽ C (∥ u_{0} ∥_{H_{0}^{1} (Ω)} + ∥ f ∥_{L^{2} (0, T; L^{2} (Ω))}) \end{matrix}$
が成立する. ここで, $C$ は $α, μ, λ, T, N, ∥ b ∥_{L^{\infty} (Ω)}, ∥ c ∥_{L^{\infty} (Ω)}, max_{i, j} ∥ \nabla a_{i, j} ∥_{L^{\infty} (Ω)}$ , ポアンカレ定数, $\partial Ω$ の $C^{2}$ ノルムに依存する定数である. さらに, $α \in (\frac{1}{2}, 1)$ のとき, $u \in C ([0, T]; L^{2} (Ω))$ である.

定理5を示すため, 次の補題を証明する.

$u_{0} \in H_{0}^{1} (Ω)$ , $f \in L^{2} (0, T; L^{2} (Ω))$ , $b, c \in L^{\infty} (Ω)$ とし, $max_{i, j} ∥ \nabla a_{i, j} ∥_{L^{\infty} (Ω)} < \infty$ と仮定する. このとき, 各 $t \in [0, T]$ と $n \in N$ に対して, 近似解 $u_{n}$ は
$\begin{matrix} I^{1 - α} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \int_{0}^{s} \frac{∥ \nabla u_{n} (s) - u_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d τ d s \\ + \int_{0}^{t} g_{α} (s) ∥ u_{n} (s) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} d s + λ \int_{0}^{t} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s \\ (23) & ⩽ \frac{t^{1 - α}}{Γ (2 - α)} ∥ \nabla u_{0} ∥_{L^{2} (Ω)}^{2} + \frac{4}{λ} \int_{0}^{t} ∥ f_{n} (s) ∥_{L^{2} (Ω)}^{2} d s + C_{1} (\int_{0}^{t} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s + \int_{0}^{t} ∥ u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s) \end{matrix}$
をみたす. ここで, $C_{1}$ は $α, λ, ∥ b ∥_{L^{\infty} (Ω)}, ∥ c ∥_{L^{\infty} (Ω)}, max_{i, j} ∥ \nabla a_{i, j} ∥_{L^{\infty} (Ω)}$ , $\partial Ω$ の $C^{2}$ ノルムに依存する定数である.

式 $(10)$ の両辺に $λ_{m} d_{n, m}$ をかけて $m = 1$ から $n$ まで和をとると,
$\begin{matrix} - \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) Δ u_{n} (x, t) d x - \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} u_{n} (x, t) \partial_{i} Δ u_{n} (x, t) d x \\ = - \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} u_{n} (x, t) Δ u_{n} (x, t) d x - \int_{Ω} c_{n} (x) u_{n} (x, t) Δ u_{n} (x, t) d x - \int_{Ω} f_{n} (x, t) Δ u_{n} (x, t) d x \end{matrix}$
となる. 左辺第1項は, $Δ u_{n} = 0$ , $_{0}^{c} D_{t}^{α} u_{n} = 0$ on $\partial Ω$ と $\nabla u_{n} (x, \cdot) \in A C [0, T]$ より
$\begin{aligned} - \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) Δ u_{n} (x, t) d x & = \int_{Ω}_{0}^{c} D_{t}^{α} \nabla u_{n} (x, t) \cdot \nabla u_{n} (x, t) d x \\ = \frac{1}{2}_{0}^{c} D_{t}^{α} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{2 Γ (1 - α)} \int_{0}^{t} \frac{∥ \nabla u_{n} (t) - \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2}}{(t - s)^{α + 1}} d s + \frac{1}{2} g_{α} (t) ∥ \nabla u_{n} (t) - \nabla u_{n, 0} ∥_{L^{2} (Ω)}^{2} \end{aligned}$
となる. 左辺第2項は, 3のProposition 9より
$\begin{aligned} - \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{j} u_{n} (x, t) \partial_{i} Δ u_{n} (x, t) d x & = \sum_{i, j = 1}^{N} \int_{Ω} \partial_{i} (a_{i, j} (x) \partial_{j} u_{n} (x, t)) Δ u_{n} (x, t) d x \\ ⩾ \frac{λ}{4} ∥ \nabla^{2} u_{n} (t) ∥_{L^{2} (Ω)}^{2} - C ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} \end{aligned}$
という評価が得られる. ここで, $C$ は $max_{i, j} ∥ \nabla a_{i, j} ∥_{L^{\infty} (Ω)}$ と $\partial Ω$ の $C^{2}$ ノルムに依存する定数である. 一方右辺は補題4と同様にして
$\begin{matrix} - \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} u_{n} (x, t) Δ u_{n} (x, t) d x - \int_{Ω} c_{n} (x) u_{n} (x, t) Δ u_{n} (x, t) d x - \int_{Ω} f_{n} (x, t) Δ u_{n} (x, t) d x \\ ⩽ \frac{2}{λ} ∥ b ∥_{L^{\infty} (Ω)}^{2} ∥ \nabla u_{n} (t) ∥_{L^{\infty} (Ω)}^{2} + \frac{2}{λ} ∥ c ∥_{L^{\infty} (Ω)}^{2} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} \\ + \frac{2}{λ} ∥ f_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{3 λ}{16} ∥ \nabla^{2} u_{n} (t) ∥_{L^{2} (Ω)}^{2} \end{matrix}$
と評価できるので, 以上をまとめると
$\begin{matrix} _{0}^{c} D_{t}^{α} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \frac{∥ \nabla u_{n} (t) - \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2}}{(t - s)^{α + 1}} d s \\ + g_{α} (t) ∥ \nabla u_{n} (t) - \nabla u_{n, 0} ∥_{L^{2} (Ω)}^{2} + \frac{λ}{8} ∥ \nabla^{2} u_{n} (t) ∥_{L^{2} (Ω)}^{2} \\ ⩽ \frac{4}{λ} ∥ f_{n} (t) ∥_{L^{2} (Ω)}^{2} + C (∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2}) \end{matrix}$
である. よって, 両辺を $0$ から $t$ まで積分すると, 式 $(23)$ が得られる.

定理5の証明

定理1と同様の議論をすることで
${‖ \frac{d}{d t} I^{1 - α} [u - u_{0}] ‖}_{L^{2} (0, T; L^{2} (Ω))} ⩽ C (∥ u_{0} ∥_{H_{0}^{1} (Ω)} + ∥ f ∥_{L^{2} (0, T; L^{2} (Ω))})$
と評価できる. したがって, 弱コンパクト性定理より
$u_{n_{k}} \to u w e a k l y i n L^{2} (0, T; H^{2} (Ω)) a s k \to \infty,$
$I^{1 - α} [u_{n_{k}} - u_{n_{k}, 0}] \to w w e a k l y i n H^{1} (0, T; L^{2} (Ω)) a s k \to \infty$
となる部分列 ${u_{n_{k}}}$ と $u \in L^{2} (0, T; H^{2} (Ω)) \cap H^{\frac{α}{2}} (0, T; H_{0}^{1} (Ω))$ と $w \in_{0} H^{1} (0, T; L^{2} (Ω))$ が存在する. 再度定理1と同様の議論をすることで得られた $u$ が等式 $(3)$ をみたしていること, 一意性と $L^{2}$ 連続性が示せる. 以上で定理の証明が完了した.

$L^{\infty}$ -Boundness of $L^{2} (Ω)$ -solutions

本章では, 次の定理を証明する.

$u_{0} \in L^{2} (Ω)$ , $f \in L^{\infty} (0, T; H^{- 1} (Ω))$ , $b, c \in L^{\infty} (Ω)$ とし, 式 $(2)$ が成立すると仮定する. このとき, 問題 $(1)$ の弱解
$u \in L^{\infty} (0, T; L^{2} (Ω)) \cap L^{2} (0, T; H_{0}^{1} (Ω)) \cap H^{\frac{α}{2}} (0, T; L^{2} (Ω))$
が一意に存在し,
$\begin{matrix} {‖ I^{1 - α} [u - u_{0}] ‖}_{H^{1} (0, T; H^{- 1} (Ω))} + ∥ u ∥_{L^{\infty} (0, T; L^{2} (Ω))} + ∥ u ∥_{L^{2} (0, T; H_{0}^{1} (Ω))} + ∥ u ∥_{H^{\frac{α}{2}} (0, T; L^{2} (Ω))} \\ (24) & ⩽ C (∥ u_{0} ∥_{L^{2} (Ω)} + ∥ f ∥_{L^{\infty} (0, T; H^{- 1} (Ω))}) \end{matrix}$
が成立する.
さらに, $u_{0} \in H_{0}^{1} (Ω)$ , $f \in L^{\infty} (0, T; L^{2} (Ω))$ のとき, $I^{1 - α} [u - u_{0}] \in W^{1, \infty} (0, T; H^{- 1} (Ω))$ となる弱解
$u \in L^{\infty} (0, T; H_{0}^{1} (Ω)) \cap L^{2} (0, T; H^{2} (Ω)) \cap H^{\frac{α}{2}} (0, T; H_{0}^{1} (Ω))$
が一意に存在し,
$\begin{matrix} {‖ I^{1 - α} [u - u_{0}] ‖}_{W^{1, \infty} (0, T; H^{- 1} (Ω))} + ∥ u ∥_{L^{\infty} (0, T; H_{0}^{1} (Ω))} + ∥ u ∥_{L^{2} (0, T; H^{2} (Ω))} + ∥ u ∥_{H^{\frac{α}{2}} (0, T; H_{0}^{1} (Ω))} \\ (25) & ⩽ C (∥ u_{0} ∥_{H_{0}^{1} (Ω)} + ∥ f ∥_{L^{\infty} (0, T; L^{2} (Ω))}) \end{matrix}$
が成立する.

まず式 $(24)$ を示す. 式 $(15)$ の両辺を $0$ から $t$ まで $α$ 階Riemann-Liouville積分すると,
$\begin{matrix} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (α) Γ (1 - α)} \int_{0}^{t} (t - s)^{α - 1} \int_{0}^{s} \frac{∥ u_{n} (s) - u_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d τ d s \\ + \frac{1}{Γ (α)} \int_{0}^{t} (t - s)^{α - 1} g_{α} (s) ∥ u_{n} (s) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} d s + \frac{1}{Γ (α)} \int_{0}^{t} (t - s)^{α - 1} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s \\ ⩽ ∥ u_{0} ∥_{L^{2} (Ω)}^{2} + C I^{α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{1}{λ} I^{α} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2} \end{matrix}$
となる. よって,
$\begin{aligned} C I^{α} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} & ⩽ ∥ u_{0} ∥_{L^{2} (Ω)}^{2} \sum_{k = 0}^{\infty} C^{k + 1} I^{α} (\frac{t^{α k}}{Γ (α k + 1)}) + \frac{1}{λ} \sum_{k = 0}^{\infty} C^{k + 1} I^{α (k + 2)} ∥ f_{n} (t) ∥_{H^{- 1} (Ω)}^{2} \\ ⩽ ∥ u_{0} ∥_{L^{2} (Ω)}^{2} \sum_{k = 0}^{\infty} \frac{C^{k + 1} t^{α (k + 1)}}{Γ (α (k + 1) + 2)} + \frac{1}{λ} ∥ f_{n} (t) ∥_{L^{\infty} (0, t; H^{- 1} (Ω))}^{2} \sum_{k = 0}^{\infty} \frac{C^{k + 1} t^{α (k + 2)}}{Γ (α (k + 2) + 1)} \\ = ∥ u_{0} ∥_{L^{2} (Ω)}^{2} C t^{α} E_{α, α + 2} (C t^{α}) + \frac{C}{λ} ∥ f_{n} (t) ∥_{L^{\infty} (0, t; H^{- 1} (Ω))}^{2} t^{2 α} E_{α, 2 α + 1} (C t^{α}) \end{aligned}$
と評価できるので,
$\begin{matrix} sup_{t \in [0, T]} ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α T^{α - 1}}{Γ (α) Γ (1 - α)} \int_{0}^{t} \int_{0}^{s} \frac{∥ u_{n} (s) - u_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d τ d s \\ + \frac{1}{Γ (α)} \int_{0}^{t} (t - s)^{α - 1} g_{α} (s) ∥ u_{n} (s) - u_{n, 0} ∥_{L^{2} (Ω)}^{2} d s + \frac{T^{α - 1}}{Γ (α)} \int_{0}^{t} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s \\ (26) & ⩽ C_{3} (∥ u_{0} ∥_{L^{2} (Ω)}^{2} + ∥ f ∥_{L^{\infty} (0, T; H^{- 1} (Ω))}^{2}) \end{matrix}$
が得られる. よって, $n_{k} \to \infty$ とすれば, $u \in L^{\infty} (0, T; L^{2} (Ω)) \cap L^{2} (0, T; H_{0}^{1} (Ω)) \cap H^{\frac{α}{2}} (0, T; L^{2} (Ω))$ が得られる.

次に, 式 $(25)$ を示す. 補題6の最後の不等式を $0$ から $t$ まで $α$ 階Riemann-Liouville積分すると,
$\begin{matrix} ∥ \nabla u_{n} (t) ∥_{H_{0}^{1} (Ω)}^{2} + \frac{α T^{α - 1}}{Γ (α) Γ (1 - α)} \int_{0}^{t} \int_{0}^{s} \frac{∥ \nabla u_{n} (s) - \nabla u_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d τ d s \\ + \frac{1}{Γ (α)} \int_{0}^{t} (t - s)^{α - 1} g_{α} (s) ∥ \nabla u_{n} (t) - \nabla u_{n, 0} ∥_{L^{2} (Ω)}^{2} d s + \frac{1}{α} \int_{0}^{t} (t - s)^{α - 1} ∥ \nabla^{2} u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s \\ ⩽ C (∥ \nabla u_{0} ∥_{L^{2} (Ω)}^{2} + \int_{0}^{t} (t - s)^{α - 1} ∥ \nabla u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s + \int_{0}^{t} (t - s)^{α - 1} ∥ u_{n} (s) ∥_{L^{2} (Ω)}^{2} d s + \int_{0}^{t} (t - s)^{α - 1} ∥ f_{n} (s) ∥_{L^{2} (Ω)}^{2} d s) \end{matrix}$
と評価でき, 右辺は式 $(26)$ より有界となる. したがって, $n_{k} \to \infty$ とすれば $u \in L^{\infty} (0, T; H_{0}^{1} (Ω)) \cap L^{2} (0, T; H^{2} (Ω)) \cap H^{\frac{α}{2}} (0, T; H_{0}^{1} (Ω))$ が得られる. 次に $_{0}^{c} D_{t}^{α} u_{n}$ の評価を行う. 定理1と同様にすると,
$\begin{matrix} | \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) w_{n} (x) d x | ⩽ μ N^{2} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)} ∥ \nabla w_{n} (x) ∥_{L^{2} (Ω)} + ∥ b ∥_{L^{\infty} (Ω)} ∥ \nabla u_{n} (t) ∥_{L^{2} (Ω)} ∥ w_{n} ∥_{L^{2} (Ω)} \\ (27) & + ∥ c ∥_{L^{\infty} (Ω)} ∥ u_{n} (t) ∥_{L^{2} (Ω)} ∥ w_{n} ∥_{L^{2} (Ω)} + ∥ f_{n} (t) ∥_{L^{2} (Ω)} ∥ w_{n} ∥_{H_{0}^{1} (Ω)} \end{matrix}$
となる. この右辺は, $L^{\infty} (0, T)$ の意味で一様有界である. よって,
${‖ \frac{d}{d t} I^{1 - α} [u - u_{0}] ‖}_{L^{\infty} (0, T; H^{- 1} (Ω))} ⩽ C (∥ u_{0} ∥_{H_{0}^{1} (Ω)}^{2} + ∥ f ∥_{L^{\infty} (0, T; L^{2} (Ω))})$
が得られる. 以上で定理の証明が完了した.

Improved regularity

$u_{0} \in H_{0}^{1} (Ω) \cap H^{2} (Ω)$ , $f \in H^{1} (0, T; L^{2} (Ω))$ , $b, c \in L^{\infty} (Ω)$ とし, 式 $(2)$ , $max_{i, j} ∥ \nabla a_{i, j} ∥_{L^{\infty} (Ω)} < \infty$ が成立すると仮定する. このとき, $I^{1 - α} [u - u_{0}] \in W^{1, \infty} (0, T; L^{2} (Ω)) \cap_{0} H^{1} (0, T; H_{0}^{1} (Ω))$ , $\partial_{t}^{α} (u - u_{0}) \in H^{\frac{α}{2}} (0, T; L^{2} (Ω))$ , $_{0}^{c} D_{t}^{α} (\partial_{t}^{α} (u - u_{0})) \in L^{2} (0, T; H^{- 1} (Ω))$ となる問題 $(1)$ の弱解 $u \in L^{\infty} (0, T; H^{2} (Ω))$ が一意に存在する. さらに
$\begin{matrix} {‖ I^{1 - α} [u - u_{0}] ‖}_{W^{1, \infty} (0, T; L^{2} (Ω))} + {‖_{0}^{c} D_{t}^{α} u ‖}_{L^{2} (0, T; H_{0}^{1} (Ω))} + {‖_{0}^{c} D_{t}^{α} u ‖}_{H^{\frac{α}{2}} (0, T; L^{2} (Ω))} \\ + {‖_{0}^{c} D_{t}^{α} (\partial_{t}^{α} (u - u_{0})) ‖}_{L^{2} (0, T; H^{- 1} (Ω))} + ∥ u ∥_{L^{\infty} (0, T; H^{2} (Ω))} \\ ⩽ C (∥ u_{0} ∥_{H^{2} (Ω)} + ∥ f ∥_{H^{1} (0, T; L^{2} (Ω))}) \end{matrix}$
が成立する. ここで, $C$ は $α, T, λ, μ, N, ∥ b ∥_{L^{\infty} (Ω)}, ∥ b ∥_{L^{\infty} (Ω)}$ に依存する定数である. また, $α \in (\frac{1}{2}, 1)$ のとき $I^{1 - α} [u - u_{0}] \in C^{1} ([0, T]; H^{- 1} (Ω))$ である.

$_{0}^{c} D_{t}^{α} u = \tilde{u}$ とする. 式 $(10)$ の両辺を $α$ 階Caputo微分し, $_{0}^{c} D_{t}^{α} d_{n, m}$ をかけて $m = 1$ から $n$ まで和をとると,
$\begin{matrix} \int_{Ω}_{0}^{c} D_{t}^{α} {\tilde{u}}_{n} (x, t) {\tilde{u}}_{n} (x, t) d x + \sum_{i, j = 1}^{N} \int_{Ω} a_{i, j} (x) \partial_{i} {\tilde{u}}_{n} (x, t) \partial_{i} u_{n} (x, t) d x \\ = \sum_{j = 1}^{N} \int_{Ω} b_{j} (x) \partial_{j} {\tilde{u}}_{n} (x, t) {\tilde{u}}_{n} (x, t) d x + \int_{Ω} c (x) | {\tilde{u}}_{n} (x, t) |^{2} d x + \int_{Ω} \tilde{f_{n}} (x, t) \tilde{u_{n}} (x, t) \end{matrix}$
となるので, 補題3より
$\begin{matrix} (28) & _{0}^{c} D_{t}^{α} ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (1 - α)} \int_{0}^{t} \frac{∥ {\tilde{u}}_{n} (t) - {\tilde{u}}_{n} (s) ∥_{L^{2} (Ω)}^{2}}{(t - s)^{α + 1}} d s + λ ∥ \nabla {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ C ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ {\tilde{f}}_{n} (t) ∥_{L^{2} (Ω)}^{2} \end{matrix}$
となる.
ここで, Holderの不等式より,
$\begin{aligned} ∥_{0}^{c} D_{t}^{α} f_{n} (t) ∥_{L^{2} (Ω)}^{2} & ⩽ {(I^{1 - α} ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)})}^{2} \\ = {(\int_{0}^{t} g_{α}^{\frac{1}{2}} (t - s) g_{α}^{\frac{1}{2}} (t - s) ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)} d s)}^{2} \\ ⩽ (\int_{0}^{t} g_{α} (t - s) d s) (\int_{0}^{t} g_{α} (t - s) ∥ f_{n}^{'} (s) ∥_{L^{2} (Ω)}^{2} d s) \\ ⩽ \frac{t^{1 - α}}{Γ (2 - α)} I^{1 - α} ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)}^{2} \end{aligned}$
と評価できる. よって,
$_{0}^{c} D_{t}^{α} ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ C ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{t^{1 - α}}{Γ (2 - α)} I^{1 - α} ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)}^{2}$
より両辺を $0$ から $t$ まで $α$ 階Riemann-Liouville積分すると, $I^{α} I^{1 - α} = I$ から
$∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ ∥ {\tilde{u}}_{n} (0) ∥_{L^{2} (Ω)} + C I^{α} ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + q (t) I ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)}^{2}$
であるので, Gronwallの不等式より
$∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} ⩽ ∥ {\tilde{u}}_{n} (0) ∥_{L^{2} (Ω)} \sum_{k = 0}^{\infty} \frac{C^{k} t^{α k}}{Γ (α k + 1)} + q (t) \sum_{k = 0}^{\infty} C^{k} I^{α k + 1} ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)}^{2}$
が得られる. したがって, 式 $(28)$ の両辺を $0$ から $t$ まで $α$ 階Riemann-Liouville積分すると,
$\begin{matrix} ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + \frac{α}{Γ (α) Γ (1 - α)} \int_{0}^{t} (t - s)^{α - 1} \int_{0}^{s} \frac{∥ {\tilde{u}}_{n} (s) - {\tilde{u}}_{n} (τ) ∥_{L^{2} (Ω)}^{2}}{(s - τ)^{α + 1}} d s + \frac{1}{Γ (α)} \int_{0}^{t} (t - s)^{α - 1} ∥ \nabla {\tilde{u}}_{n} (s) ∥_{L^{2} (Ω)}^{2} \\ ⩽ ∥ {\tilde{u}}_{0} ∥_{L^{2} (Ω)}^{2} + C I^{α} ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + q (t) I ∥ f_{n}^{'} (t) ∥_{L^{2} (Ω)}^{2} \end{matrix}$
となる.
$∥ {\tilde{u}}_{n} (0) ∥_{L^{2} (Ω)} ⩽ C ∥ u_{0} ∥_{H^{2} (Ω)}$
であり, 定理7と同様の議論をすれば, $n_{k} \to \infty$ より
$_{0}^{c} D_{t}^{α} u \in L^{\infty} (0, T; L^{2} (Ω)) \cap L^{2} (0, T; H_{0}^{1} (Ω)) \cap H^{\frac{α}{2}} (0, T; L^{2} (Ω))$
が得られる. 次に $u \in L^{\infty} (0, T; H^{2} (Ω))$ を示す.
$(10)$ の両辺に $λ_{m} d_{n, m}$ をかけて $m = 1$ から $n$ まで和をとり, $Δ u_{n} = 0$ on $\partial Ω$ より
$- \int_{Ω} L u_{n} (x, t) Δ u_{n} (x, t) d x = \int_{Ω}_{0}^{c} D_{t}^{α} u_{n} (x, t) Δ u_{n} (x, t) d x - \int_{Ω} f_{n} (x, t) Δ u_{n} (x, t) d x$
となる. よって, elliptic regularityより
$∥ u_{n} (t) ∥_{H^{2} (Ω)}^{2} ⩽ - \int_{Ω} L u_{n} (x, t) Δ u_{n} (x, t) d x + C ∥ u_{n} (t) ∥_{L^{2} (Ω)}^{2}$
となるので,
$∥ u_{n} (t) ∥_{H^{2} (Ω)}^{2} ⩽ C (∥ f_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ {\tilde{u}}_{n} (t) ∥_{L^{2} (Ω)}^{2} + ∥ u_{n} (t) ∥_{L^{2} (Ω)})$
が得られる. 右辺は $L^{\infty} (0, T)$ の意味で有界であるので, $n_{k} \to \infty$ とすれば $u \in L^{\infty} (0, T; H^{2} (Ω))$ を得る. さらに, 式 $(10)$ の両辺を $α$ 階caputo微分し, $θ_{m}$ をかけて $m = 1$ から $n$ まで和をとると,
$_{0}^{c} D_{t}^{α}_{0}^{c} D_{t}^{α} u \in L^{2} (0, T; H^{- 1} (Ω))$
を得る. 以上で定理の証明が完了した.

参考文献

[1]

Eidelman, Samuil D and Kochubei, Anatoly N., Cauchy problem for fractional diffusion equations, J. Differential Equations, 2004

[2]

Evans, Lawrence C., Partial differential equations, Graduate Studies in Mathematics, American Mathematical Society, Providence, RI, 2010, xxii+749

[3]

Kubica, Adam and Yamamoto, Masahiro, Initial-boundary value problems for fractional diffusion equations with time-dependent coefficients, Fract. Calc. Appl. Anal., 2018

[4]

Luchko, Yury, Some uniqueness and existence results for the initial-boundary-value problems for the generalized time-fractional diffusion equation, Comput. Math. Appl., 2010

投稿日：2023年8月29日

更新日：2023年11月2日

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

カメ

1387

大学院では非線形拡散方程式(主にFast Diffusion, Porous Medium), 非整数階時間微分を含む拡散方程式を専攻していました. 現在は非整数階時間微分を含む拡散方程式の可解性の研究をしています.

他の人のコメント

コメントはありません。

読み込み中

カメ

非整数階時間微分を含む拡散方程式の初期値境界値問題

Initial-boundary value problems for time-fractional diffusion equations with time-independent coefficients
はじめに
Approximate solutions
Weak solutions
Regular solutions
$L^\infty$-Boundness of $L^2(\om)$-solutions
Improved regularity
参考文献

非整数階時間微分を含む拡散方程式の初期値境界値問題

Initial-boundary value problems for time-fractional diffusion equations with time-independent coefficients

はじめに

Approximate solutions

Weak solutions

Regular solutions

L∞-Boundness of L2(Ω)-solutions

Improved regularity

参考文献

この記事を高評価した人

この記事に送られたバッジ

投稿者

コメント

他の人のコメント

$L^{\infty}$ -Boundness of $L^{2} (Ω)$ -solutions