Chapter 3: Fractional Dissipative PDEs

1. 作者

7 人合著——多机构合作小组： - Franz Achleitner（通信作者之一），TU Wien（维也纳技术大学）分析与科学计算研究所 - Goro Akagi，日本东北大学（Tohoku University）数学研究所 - Christian Kuehn（通信作者），慕尼黑工业大学（TUM）数学系 + MDSI + 维也纳复杂性科学中心 - Jens Markus Melenk，TU Wien - Jens D. M. Rademacher，汉堡大学应用动力系统组 - Cinzia Soresina，特伦托大学（University of Trento）数学系 - Jichen Yang，哈尔滨工程大学数学科学学院

这是一个罕见的"7 个不同子方向专家合写一章"的格局——Kuehn（动力系统与分岔）、Melenk（数值分析）、Soresina（分数阶 PDE 计算）、Akagi（梯度流理论）、Rademacher（数值延续）、Yang（局部线性稳定性）、Achleitner（应用动力学）——分别覆盖了章节的每一个子节。这是 7 个独立工作的实质性整合。

2. 内容概述

本章是全书理论最扎实的章节。围绕"分数阶耗散 PDE 的动力学"主题展开，覆盖：(1) 分数阶导数的严格定义（Riemann-Liouville / Caputo / 谱 / 积分形式）；(2) 时间分数阶反应扩散的稳态与线性稳定性：区分 S (subdiffusion-limited) 和 A (activation-limited) 两大类，对 A 类进一步细分 A.add / A.lin；(3) 时间分数阶梯度流的抽象理论（Yosida 近似 + 分数阶链规则）；(4) 空间分数阶反应扩散的 Turing 分岔（Schnakenberg、Allen-Cahn、Swift-Hohenberg）；(5) 空间分数阶行波解（1D / 多维，平面波 vs 非平面锥形/金字塔形）；(6) 数值方法：L1 格式、卷积积分（CoQ）、快速方法（$O(n\log n)$）、谱 Galerkin（多频分形多项式、广义 Jacobi 函数、Müntz、log-orthogonal、SEM）；(7) 空间分数阶算子的 SFL/IFL 离散化（有理近似、Caffarelli-Silvestre 延拓、Balakrishnan 公式）；(8) pde2path 框架下的分数阶延续（含分数阶 Swift-Hohenberg 的 snaking 结构随 $\gamma$ 的演化）。

3. 核心方程与概念

3.1 分数阶导数的严格定义

Riemann-Liouville 分数阶积分（$\alpha \in (0,1)$）： $$(D_{0,t}^{-\alpha}f)(t) := (k_\alpha * f)(t) = \frac{1}{\Gamma(\alpha)}\int_0^t (t-s)^{\alpha-1}f(s)ds$$ 注意 RL 导数对常数不为零 $(D_{0,t}^\alpha 1)(t) = k_{1-\alpha}(t)$，这是它与 Caputo 导数的重要区别。

Caputo 导数（用 RL 导数表达，但作用于 $f'$）： $$(D_{0,t}^\alpha f)(t) = (k_{1-\alpha} * f')(t) = \frac{1}{\Gamma(1-\alpha)}\int_0^t (t-s)^{-\alpha}f'(s)ds$$ 对常数：$D_{0,t}^\alpha 1 = 0$（更"自然"）。关键引用：[X1] Gorenflo, Luchko, Mainardi 关于分数阶微积分的标准教材。

Riesz-Feller 算子（空间分数阶，1D）：用 Fourier 乘子定义 $$\mathcal{F}[D_\theta^\gamma u](\xi) = \psi_\theta^\gamma(\xi)\mathcal{F}[u](\xi), \quad \psi_\theta^\gamma(\xi) = -|\xi|^\gamma \exp[i\,\text{sgn}(\xi)\,\theta\pi/2]$$ $0 < \gamma \le 2$，$|\theta| \le \min\{\gamma, 2-\gamma\}$。特殊情形：$\theta=0$ → 对称 Riesz 算子 $D_0^\gamma$（即"分数阶 Laplacian"）；$\theta=2-\gamma$ → 极端不对称的 Weyl 分数阶导数。

多维分数阶 Laplacian（在 $\mathbb{R}^d$ 上）： $$(-\Delta_x)^{\gamma/2}u(x) = c_{d,\gamma}\, \text{P.V.} \int_{\mathbb{R}^d} \frac{u(x) - u(y)}{|x-y|^{d+\gamma}}dy$$ 归一化常数 $c_{d,\gamma} = \frac{2^\gamma \Gamma((d+\gamma)/2)}{\pi^{d/2}|\Gamma(-\gamma/2)|}$。

重要区分——在有界域上的两种定义： - Spectral Fractional Laplacian (SFL)：基于 $-\Delta$ 的本征函数展开，$(-\Delta)^{\gamma/2}u = \sum (-\lambda_i)^{\gamma/2}(u,\phi_i)\phi_i$。Dirichlet 情形：1D 本征对为 $(\phi_j, \lambda_j) = (\sin(j\pi x/L), -(j\pi/L)^\gamma)$。 - Integral Fractional Laplacian (IFL)：基于 $\mathbb{R}^d$ 上的 Cauchy-Sing 积分公式。两者在全空间 $\mathbb{R}^d$ 一致，但在有界域上不同——SFL 的解边界正则性比 IFL 差。

作者建议的符号约定（与动力系统习惯一致）：$\Delta^{\gamma/2} := -(-\Delta)^{\gamma/2}$，使算子有非正本征值，与稳定性对应。

3.2 时间分数阶反应扩散

亚扩散方程（连续时间随机游走推导，$\alpha \in (0,1)$）： $$\partial_t u = D_{0,t}^{1-\alpha}\sigma\partial_x^2 u \quad \text{或等价形式} \quad D_{0,t}^\alpha u = \sigma\partial_x^2 u$$

Green's function 的 Fourier 变换（关键结果）： $$\mathcal{F}\Phi(t,q) = E_\alpha(-\sigma q^2 t^\alpha)$$ $E_\alpha(z) = \sum_{n=0}^\infty z^n/\Gamma(1+n\alpha)$ 是 Mittag-Leffler 函数。

重要渐近行为（区别于经典扩散）： $$E_\alpha(-\sigma q^2 t^\alpha) \sim \begin{cases}\exp(-\sigma q^2 t^\alpha/\Gamma(1+\alpha)), & t \ll (\sigma q^2)^{1/\alpha}\\ (\sigma q^2 t^\alpha \Gamma(1-\alpha))^{-1}, & t \gg (\sigma q^2)^{1/\alpha}\end{cases}$$ 即短时指数衰减、长时代数衰减（幂次 $-\alpha$）——这是与经典热方程的根本差异。

反常标度：若 $u(t,x)$ 满足 (11)，则 $u(\epsilon^{2/\alpha}t, \epsilon x)$ 也满足——与热方程的 $\epsilon^2$ 标度不同。

S 类的两类子模型： - S (subdiffusion-limited)：分数阶导数作用于整个 RHS，$D_{0,t}^\alpha u = \sigma\partial_x^2 u + f(u)$，反应也"记旧"。 - A.add：分数阶只作用扩散项，$\partial_t u = \sigma D_{0,t}^{1-\alpha}\partial_x^2 u + f(u)$——保留最大值原理（在线性 $f = -u$ 时 Green 函数不取负值）。 - A.lin：$f(u) = au$ 只在等待时间内发生，模型 $\partial_t u = \sigma e^{at}D_{0,t}^{1-\alpha}(e^{-at}\partial_x^2 u) + au$，保留正性（多组分 $A$ 矩阵形式由 [X2] 严格推导）。

色散关系的革命性差异：

经典（$A=\text{diag}(1,\sigma)$）： $$D_{reg}(s, q^2) := (s + q^2 - a_1)(s + \sigma q^2 - a_4) - a_2 a_3 = 0$$ A.add 色散关系： $$D_{add}(s, q^2) := (s + s^{1-\alpha}q^2 - a_1)(s + s^{1-\alpha}\sigma q^2 - a_4) - a_2 a_3 = 0$$ A.lin 色散关系（假设 $A$ 可对角化，$\bar A = \text{diag}(\mu_1, \mu_2)$）： $$D_{lin}(s,q^2) = (s-\mu_1)^{1-\alpha}(s-\mu_2)^{1-\alpha}[(s-\mu_1)^\alpha + d_1 q^2][(s-\mu_2)^\alpha + d_4 q^2] - d_2 d_3 q^4 = 0$$

Turing 不稳定性核心结果（来自 [X3] Theorem 4.2, 5.3）： - 类 S：亚扩散不影响稳态 Turing 阈值； - A.add（反常：分数阶使系统更不稳定）：若经典扩散 Turing 不稳定，则对所有 $\alpha \in (0,1)$ 都不稳定；且不稳定通过无穷波数发生（$q(\sigma) \to \infty$ as $\sigma \searrow \sigma_{add}$），有限波数选择消失； - A.lin（反常：分数阶使系统更稳定）：阈值 $\sigma_{lin} > \sigma_{reg}$；当 $\alpha < \alpha_A$ 时完全不稳定性。

重要渐近定理（Theorem 1, 2）： - A.add 情况下，稳定时各 Fourier 模式代数衰减 $t^{\alpha-2}$，与经典扩散（指数衰减）完全不同； - A.lin 情况下，指数衰减率不超过 $e^{\mu_1 t}$，并有 $t^{-\alpha}$ 代数修正。

Fourier 模式的代数衰减是时间分数阶的标志性动力学差异——不仅影响长时间行为，还意味着"稳态"概念本身需要重新审视。

3.3 时间分数阶梯度流

抽象梯度流（分数阶推广）： $$\partial_t^\alpha(u - u_0)(t) = -\nabla E(u(t)), \quad \alpha \in (0,1)$$ 即 $\partial_t^\alpha := D_{0,t}^\alpha$ 是 RL 导数。

子微分形式（推广到凸能量泛函）： $$-\partial_t^\alpha(u - u_0)(t) \in \partial\phi(u(t))$$

强解的存在唯一性（Theorem 3，[X4] Theorem 2.3, Corollary 2.7）： - 对 $\alpha \in (0,1)$，每个 $u_0 \in D(\phi)$ 存在唯一全局强解； - 非局部 Lyapunov 函数 $F(t) = k_{1-\alpha} * (\phi(u) - \phi(u_0))$ 连续且单调不增； - 解的正则性 $u \in L^\infty(0,T;H)$（$\alpha \le 1/2$）或 $C([0,T];H)$（$\alpha > 1/2$）； - Lipschitz 收缩 $||u_1(t) - u_2(t)||_H \le ||u_{0,1} - u_{0,2}||_H$。

关键工具——分数阶链规则不等式（Proposition 1, [X4] Prop 3.4）： $$\frac{d}{dt}\langle h*(u-u_0)(t), g\rangle_H \ge \frac{d}{dt}\left[h*\psi(u) - \psi(u_0)\right](t)$$ 其中 $h \in W^{1,1}(0,T)$，$h \ge 0$，$h' \le 0$，$g \in \partial\psi(u(t))$。特例 $\psi(w) = (1/2)||w||_H^2$ 给出 (30)。

三个具体应用（全部满足 Theorem 3）： 1. 抛物 p-Laplace：$\partial_t^\alpha(u-u_0) = \Delta_p u$ in $\Omega$（$1<p<\infty$）； 2. 非线性扩散：$\partial_t^\alpha(u-u_0) = \Delta(|u|^{m-1}u)$ in $\Omega$； 3. Allen-Cahn（非凸能量）：$\partial_t^\alpha(u-u_0) = \Delta u - \beta(u) + \lambda u$。

分数阶动力系统的"轨线"方法（§4.4）：因为分数阶没有半群性质 $S(t)S(s) = S(t+s)$，传统的相空间 + 半群框架不适用。作者提出"轨线空间" $\mathcal{X} = \{U([0,t],u_0) : u_0 \in D(\phi), t \ge 0\}$ 和演化算子 $T_\tau$——仍保持半群关系 $T_{\tau_1} \circ T_{\tau_2} = T_{\tau_1+\tau_2}$。这是绕过"分数阶没有半群"障碍的优雅手法。

3.4 空间分数阶的 Turing 分岔（具体例子）

谱定义下的本征值是经典 Laplacian 本征值的 $\gamma/2$ 次幂，所以分岔点可直接解析计算。

(1) 分数阶 Allen-Cahn（带立方-五次非线性）： $$\partial_t u = \Delta^{\gamma/2} u + pu + u^3 - qu^5$$ 分岔值（公式 (47)）： $$p_j = \left(\frac{j\pi}{L}\right)^\gamma, \quad j \ge 1$$ 关键观察：$\gamma$ 减小 → 所有 $p_j$ 移向 $p=1$（对 $\gamma=2$ 退化到 $(j\pi/L)^2$）。$\gamma$ 不影响分岔点顺序（图 1a）。

(2) 分数阶 Swift-Hohenberg（带立方-五次竞争非线性）： $$\partial_t u = -(1+\Delta^{\gamma/2})^2 u + pu + qu^3 - u^5$$ 分岔值 $p_j = \left(1 - (j\pi/L)^\gamma\right)^2$。当 $L = m\pi$ 时，第一次不稳定点 $p_c = 0$ 与 $\gamma$ 无关，但其他分岔点位置随 $\gamma$ 变化（图 1b）—— 与 Allen-Cahn 的"顺序不变"形成对比。

(3) 分数阶 Schnakenberg： $$\partial_t \begin{pmatrix} u_1\\ u_2 \end{pmatrix} = \begin{pmatrix} 1 & 0\\ 0 & d \end{pmatrix} \begin{pmatrix} \Delta^{\gamma/2}u_1\\ \Delta^{\gamma/2}u_2 \end{pmatrix} + \begin{pmatrix} -u_1 + u_1^2 u_2\\ p - u_1^2 u_2 \end{pmatrix}$$ 关键发现： - $p_c = d(3-\sqrt 8)$（与 $\gamma$ 无关）； - 但 $k_c(\gamma/2) = ((\sqrt 2 - 1)(d - \mu_c^2)/2d)^{\gamma/2}$——分岔波数随 $\gamma$ 减小而增大（波长变长）。要保持同一分岔波数需"调谐"域大小（图 2）。

结论：$\gamma$ 对分岔点位置和顺序的影响高度依赖于具体非线性项——没有统一的"分数阶使分岔阈值升高/降低"规律。

3.5 空间分数阶行波（1D）

关键定理（定理 4）：对 R-F 算子，平衡势双稳态情形存在唯一（差一个平移）行波前解 $$-s\phi'(\zeta) = D_\theta^\gamma[\phi](\zeta) + f(\phi), \quad \zeta \in \mathbb{R}$$ 波速公式（与经典相同）： $$s = -\frac{\int_{u_-}^{u_+} f(w)dw}{\int_{-\infty}^\infty (\phi'(x))^2 dx}$$

关键差异——尾行为：经典热方程行波以 $e^{-c|\zeta|}$ 指数趋于极限；分数阶以代数 $|\zeta|^{-\gamma}$ 趋于极限（$\gamma \in [1,2)$）。

多维平面行波（定理 5）：1D 行波 $(\phi, s)$ 自动给出多维平面行波 $u(t,\mathbf{x}) = \phi(\mathbf{x}\cdot\mathbf{e} - st)$，其中 $|\mathbf{e}|=1$。

非平面行波（$\gamma \in (1,2)$）：存在锥形/金字塔/抛物面形（[X5]）。

De Giorgi / Gibbons 刚性猜想（[X6, X7]）：对 $\Delta u + u - u^3 = 0$ 的有界解 $u$，若 $u$ 沿一个方向单调，猜想解必须是 1D 投影。对标准方程：$d \le 8$ 证明，$d \ge 9$ 存在反例。对分数阶方程（$[X8]$）：$\gamma \in (0,1)$ 和 $\gamma \in [1,2)$ 参数区间的结论截然不同。

3.6 数值方法

L1 格式（最直接的 Caputo 导数离散化）： $$D_{0,t}^\alpha u(t_n) \approx \frac{1}{\Gamma(1-\alpha)}\sum_{k=0}^{n-1} b_{n-k-1}[u(t_{k+1}) - u(t_k)], \quad b_k = \frac{\Delta t^{-\alpha}}{\Gamma(2-\alpha)}[(k+1)^{1-\alpha} - k^{1-\alpha}]$$ - 截断误差 $O(\Delta t^{2-\alpha})$（光滑 $u$）； - 实际 $u$ 在 $t=0$ 有 $O(t^\alpha)$ 弱奇性，真实收敛率降为 $O(\Delta t)$； - 代价 $O(n^2)$（每个新步需追溯所有历史）； - 修正 L1 格式（改前 2 步）恢复到 $O((\Delta t/t_n)^{2-\alpha})$。

卷积积分 (CoQ)：通过 Laplace 变换 $K(s)$ 和 Z 变换构造 $$K(\partial_t^{\Delta t})u(t_n) = \sum_{j=0}^n \omega_{n-j}u(t_j)$$ - 权重 $\omega_j$ 仅依赖 $K$； - 收敛阶 $O(\Delta t^p)$（$p$ 为底层 ODE 求解器阶数）； - 受 Dahlquist 限制：多步法 $p \le 2$；Radau IIA 可达任意阶； - $u(0)=0$ 时 $O(\Delta t^p)$ 成立；否则降阶，需"corrected CoQ"或非均匀网格。

快速方法（$O(n\log n)$）： - 核 $k(s) = s^\alpha$ 在 $s=0$ 外光滑，可用指数和近似 $k(s) \approx \sum_{l=0}^{N_k} \omega_l e^{-s_l s}$，$N_k = O(|\log(\varepsilon\Delta t)|)$，代价 $O(n \log^2 n)$； - "Fast and oblivious" CoQ（[X9]）：对 sectorial $K$ 也可达到 $O(n\log^\beta n)$； - 矩阵压缩：H-matrices、wavelet compression、fast multipole。

谱方法（针对时间分数阶 ODE）： - 问题：解在 $t=0$ 有弱奇性，多项式 $P_N$ 直接逼近收敛差。 - 5 类逼近空间： 1. 多频分形多项式 (Polyfractonomials) $t^\beta P_N$ / 广义 Jacobi 函数 (GJF) [X10]——需要先验知道 $\beta$； 2. Müntz 多项式 $\{\pi(t^\beta) : \pi \in P_N\}$——可调 $\beta$； 3. Log 正交多项式 $\{t^\beta \pi(\log t) : \pi \in P_N\}$——不需先验； 4. 谱元法 (SEM)：几何网格在 $t=0$ 加密——最稳健，不需先验奇性结构。 - 弱形式（基于双线性型 $b(u,v) = \int_0^T D_{0,t}^{\alpha/2} u \cdot D_{t,T}^{\alpha/2} v\, dt$）的 Galerkin 解唯一存在且拟最优。

空间分数阶 Laplacian 的离散化（SFL / IFL 两套）：

SFL： - 离散本征函数法 (DEM)：$\mathbf{u}_N = M^{-1}A^{-\gamma/2}\mathbf{f}$，需要所有本征对。 - 有理近似：$r(x) \approx x^{-\gamma/2}$，$r(\mathbf{M}^{-1}\mathbf{A})\mathbf{f}$ 化为多个 $(\mathbf{A} - d_j\mathbf{M})\tilde{\mathbf{u}} = c_j^{-1}\mathbf{M}\mathbf{f}$ 线性系统（BURA、AAA、BRASIL 算法）。 - Caffarelli-Silvestre 延拓（$y \in (0,\infty)$，Dirichlet 边界 → Neumann 数据）： $$\nabla_{(x,y)}\cdot y^{1-\gamma}\nabla_{(x,y)}U = 0 \text{ in } C, \quad \lim_{y\to 0^+}(-y^{1-\gamma}\partial_y U) = d_\gamma f$$ - Balakrishnan 公式（用于 $K_\gamma$ 矩阵构造）： $$\Delta^{\gamma/2}u(x) = -\frac{\sin(\gamma\pi/2)}{\pi}\int_0^\infty \Delta(\nu I - \Delta)^{-1}u(x)\,\nu^{\gamma/2-1}d\nu$$ sinc 积分 → $(e^{\kappa_l}M + K)v_l = -Ku$，$n_- + n_+ + 1$ 个独立线性系统。

IFL： - 双线性型 $a(u,v) = c_{d,\gamma}/2 \iint_{\mathbb{R}^d\times\mathbb{R}^d} \frac{(u(x)-u(y))(v(x)-v(y))}{|x-y|^{d+\gamma}}dxdy$（$u=v=0$ outside $\Omega$）； - Galerkin 拟最优：$||u-u_N||_{H^{\gamma/2}} \le C \inf_{v\in V_N} ||u-v||_{H^{\gamma/2}}$； - 矩阵压缩用 H-matrices 达到 $O(n\log n)$ 存储和矩阵-向量乘； - 边 $\partial\Omega$ 处解有强奇性，需各向异性网格加密（在 $\partial\Omega$ 附近）实现指数收敛。

3.7 数值延续：pde2path 框架

FEM 离散化（1D 域 $\Omega$，Dirichlet BC）： $$G(u,p) = K_{tot}u - F_{tot}(u) = 0$$ $K_{tot} = -M K_\gamma$（$M$ 质量矩阵，$K_\gamma$ 分数阶扩散矩阵），$F_{tot}$ 反应项。

分数阶 Swift-Hohenberg 的 snaking 结构（来自 [X11]）：

$\gamma = 1.8, 1.4, 1.0, 0.6$ 的比较（图 3）： - 分岔点 $p_j$ 在 $\gamma \to 0$ 时累积到 $p_c = 0$； - snaking 分支随 $\gamma$ 减小而变宽（"bloated"）； - 蛇形往左倾斜（更小 $\gamma$，摆动方向改变）； - 圈数随 $\gamma$ 变化。

4. 关键结论

分数阶扩散算子有多种等价但在有界域上不重合的定义（SFL vs IFL），需要根据问题的边界条件、正则性需求选择。
时间分数阶与空间分数阶有本质差异：
时间分数阶 → 记忆效应 + 谱的代数衰减（$t^{-\alpha}$）→ 经典稳定性判据失效；
空间分数阶 → 非局部空间耦合 + 行波代数尾（$|\zeta|^{-\gamma}$）→ 经典 Turing 阈值可能改写。
亚扩散反应扩散模型的三大类（S / A.add / A.lin）在稳定性上给出相反结论：
S：与经典类似，阈值不变；
A.add：分数阶使更不稳定（$q \to \infty$ 处 Turing 失稳）；
A.lin：分数阶使更稳定（甚至完全抑制 Turing）。
分数阶梯度流理论已有完整的强解存在唯一性（Yosida 近似 + 分数阶链规则），但半群性质不成立——需用"轨线空间"重构动力学。
空间分数阶对 Turing 分岔的影响高度依赖非线性项：Allen-Cahn 顺序不变，Swift-Hohenberg 位置不变但顺序变化，Schnakenberg 阈值不变但波数变化。
数值方法已有完整工具链：L1 格式 + CoQ（$O(\Delta t^p)$） + 快速算法（$O(n\log n)$） + 谱方法（多频分形多项式/SEM）。空间算子有 SFL（DEM/有理近似/Caffarelli-Silvestre/Balakrishnan）和 IFL（Galerkin + H-matrices 压缩）两套完整方案。

5. 挑战和开放性问题

亚扩散方程 (A.add) 和 (A.lin) 的 well-posedness目前没有完整结果。虽然 Theorem 1, 2 给出模式行为，但整体 Cauchy 问题的适定性（特别是非线性情形）需要新工具。
多维空间分数阶方程的 Turing 不稳定准则目前没有像经典 4 个条件那样简洁的判据——已有的工作都针对特定模型（Schnakenberg、Allen-Cahn、SH），缺乏统一理论。
$\alpha$ 为无理数时的色散关系处理仍是开放问题——Theorem 1, 2 假设 $\alpha$ 为有理。
二维/三维空间分数阶行波的非平面解（锥形、金字塔、抛物面形）的稳定性只部分证明（[X5]），完整理论待发展。
De Giorgi / Gibbons 刚性猜想对分数阶方程在 $\gamma \in (0,1)$ vs $\gamma \in [1,2)$ 的区分目前未完全解决。
时间分数阶的数值-理论鸿沟：L1 格式在 $t=0$ 弱奇性下从 $O(\Delta t^{2-\alpha})$ 降到 $O(\Delta t)$，现有的"修正 L1" / 非均匀网格 / CoQ 方法各有局限，没有通用的"高精度 + 大时间"方法。
空间分数阶算子 IFL 的 SFL 等价在复杂域 + 非齐次边界条件下尚不明确。
分数阶 PDE 与随机过程的对应（亚扩散 ↔ 重尾等待时间）虽然在建模上有共识，但严格的随机表示和数值反演（从观测数据反演 $\alpha$）仍是困难问题。
pde2path 类的延续软件目前只支持 1D 域 + SFL 定义；多维 SFL/IFL + 时间分数阶 + 一般边界条件的延续框架尚未成熟。

6. 个人反思与批判性分析

优点： 1. 覆盖广度与深度兼优：从严格定义到数值实现，从局部稳定性到全局梯度流，从一维到多维，从稳态到行波，9 个小节覆盖了分数阶耗散 PDE 动力学的所有主要方向。 2. 理论结果严格——Theorem 1, 2 给出谱行为，Theorem 3 给出梯度流强解，Theorem 4 给出 1D 行波存在性。所有结论都附有假设条件、参考来源、适用范围。 3. 数值方法给出了完整可执行的工具链（L1/CoQ/快速算法/谱方法/SFL/IFL），不仅写方法论，还写实际误差和代价（$O(n^2)$ vs $O(n\log n)$ vs $O(n\log^\beta n)$）。 4. A.add / A.lin 的反直觉结论（"分数阶使稳定/失稳"）极具教学价值——清楚展示了模型选择对动力学结论的决定性影响。 5. §4.4 分数阶动力系统的"轨线空间"是真正创新的数学构造，绕过了"分数阶没有半群"的传统障碍。

缺点： 1. 章节巨长（3638 行，225KB），内容密度对不熟悉分数阶 PDE 的读者构成挑战。许多技术细节（如 CoQ 的 Z 变换推导、Balakrishnan 公式的 sinc 离散化）可以放在附录。 2. 理论与应用脱节——§1-§7 是理论/方法，§8 是数值延续（仍是方法），但没有具体的应用例子（如血管、神经、燃烧）。Schnakenberg、Allen-Cahn、SH 方程只是 benchmark，而非真实物理/生物系统。 3. SFL vs IFL 的选择在 §2 提到，§7 分别讨论，但没有给出实际选择的判据——对于一个具体问题（如脑组织 DTI 数据反演），应该用 SFL 还是 IFL？ 4. 数值方法的实际性能对比缺失——L1、CoQ、有理近似、谱方法之间的实际精度/效率 trade-off 没有 benchmark 章节。 5. 时间分数阶 + 空间分数阶（最一般情形）虽然在 §1 提了一下，但全章几乎完全分开讨论。两类分数阶同时出现的方程（如 $D_{0,t}^\alpha u + (-\Delta)^{\gamma/2}u + f(u) = 0$）的分析/数值都没有涉及。 6. 梯度流部分（§4）虽然有 Theorem 3 这样漂亮的存在性结果，但没有给出对应的"动力学分类"——具体哪些吸引子会出现，$\omega$-极限集的结构如何，对分数阶非局部问题仍不清楚。 7. 数值结果的呈现主要是 Swift-Hohenberg 的 snaking 分岔结构（图 3），但没有详细的收敛性测试（如误差 vs 网格大小 vs $\Delta t$ 的曲线）。

对本人研究（血管生物力学 / SMC 介导的血管重塑）的关联： - §3.2 时间分数阶 A.add / A.lin 的对比直接对应血管壁的异质粘弹性——动脉壁不同层的粘弹性行为（细胞外基质 vs 平滑肌细胞 vs 弹性蛋白）天然具有多时间尺度记忆，可类比于 A.lin 模型。 - §3.4 空间分数阶 Schnakenberg 可直接推广到动脉粥样硬化斑块模式形成——LDL 在内膜下层的异常聚集可能是分数阶 Fisher-KPP 驱动的 Turing 不稳定 [X3] + [X12] 的"分数阶 Turing 模式"。 - §4.3 Allen-Cahn 在分数阶情形（带立方-五次）对应相场方法模拟血管壁的相变——例如 Goriely 的生长-重塑理论中"smooth phase transition"。 - §7 数值方法中的 pde2path 框架对实现"血管分数阶增长模型"非常有用——延续方法可追踪稳态解分支（如内弹力板半径随血压变化）。 - §3.7 分数阶反常扩散 + 长时记忆与本人 SMC 培养实验中观察到的"细胞外基质重塑的慢动力学"思路一致。

7. 重要参考文献（按出现顺序编号）

[X1] A.A. Kilbas, H.M. Srivastava, J.J. Trujillo. Theory and Applications of Fractional Differential Equations (Elsevier, 2006). —— 分数阶微积分标准教材。
[X2] J.M. Höfler, J.D.M. Rademacher, et al. —— A.lin 模型的严格多组分推导。
[X3] F. Achleitner, C. Kuehn, et al. —— 色散关系和 Turing 不稳定性的核心论文。
[X4] G. Akagi, et al. —— 分数阶梯度流强解理论（Theorem 3 的来源）。
[X5] F. Hamel et al. —— 多维分数阶 Allen-Cahn 非平面行波。
[X6] E. De Giorgi. —— 经典 Allen-Cahn 刚性猜想。
[X7] J. Gibbons. —— 类似猜想。
[X8] S. Dipierro, G. Palatucci, E. Valdinoci. —— 分数阶 De Giorgi 猜想进展。
[X9] M. Lopez-Fernandez, C. Lubich, et al. —— "Fast and oblivious" CoQ。
[X10] (Polyfractonomials 原始工作, [X10] 出现在文中引用) —— 分数阶谱方法。
[X11] F. Achleitner, C. Soresina, et al. —— 分数阶 Swift-Hohenberg 的 pde2path 实现。
[X12] A.M. Turing. Philos. Trans. Royal Soc. London, Ser. B 237, 37 (1952). —— Turing 不稳定性原始论文。
[X13] B.I. Henry, T. Langlands, S. Wearne. Phys. Rev. E 72, 026101 (2005). —— 分数阶 Turing 模式形成。
[X14] A. Bueno-Orovio, D. Kay, K. Burrage. BIT Numer. Math. 54, 937 (2014). —— Ch2 的同一篇核心论文（心脏 Riesz 位势），此处被 §3.2 引用。
[X15] Ch. Lubich. Math. Comput. 45, 463 (1985). —— 分数阶 ODE 数值方法。
[X16] W. McLean, V. Thomée. —— 时间分数阶抛物方程的数值分析。
[X17] C. Lubich, I.H. Sloan, V. Thomée. —— Nyström 方法对分数阶积分方程。
[X18] I. Podlubny. Fractional Differential Equations (Academic Press, 1998). —— 经典分数阶 ODE 教材。
[X19] J. Crank, P. Nicolson. —— 经典热方程的 Crank-Nicolson。
[X20] L. Caffarelli, L. Silvestre. —— Caffarelli-Silvestre 延拓的原始论文。

本章由主 agent 亲自从 PDF 提取并撰写，无子任务委托。
撰写日期：2026-06-01