Chapter 5: Traveling Waves in Fractional Models

1. 作者

Dmitry E. Pelinovsky，加拿大 McMaster University 数学系教授，Email: dmpeli@math.mcmaster.ca。Pelinovsky 是分数阶色散方程理论研究的领军人物之一，特别是在分数阶 KdV 方程的适定性、孤子稳定性、变分理论方面有突出贡献。McMaster 是加拿大分数阶 PDE 研究的重镇之一（与同日发布的分数阶 Schrödinger 团队并列）。本章是 Ch6（数值方法）和 Ch8（可积方程）的理论前奏——专注于纯理论方法（小振幅展开 + 不动点 + 变分）证明行波的存在与稳定性。

2. 内容概述

本章是全书理论性最强的章节，专门讨论分数阶 KdV/mKdV 方程的周期行波解的存在性与稳定性。两大模型： - 分数阶 KdV（二次非线性）：$u_t + uu_x = (-\Delta)^{\alpha/2} u_x$； - 分数阶 mKdV（三次非线性）：$u_t + u^2u_x = (-\Delta)^{\alpha/2} u_x$。

$\alpha = 2$ 退化到经典 KdV/mKdV，$\alpha = 1$ 对应 Benjamin-Ono (BO) 方程（用 Hilbert 变换）。本章用 3 种互补方法分析：(1) 小振幅展开（Lyapunov-Schmidt 约化）—— 得到 $\alpha$ 阈值 $\alpha_0 \approx 0.585$（KdV）和 $\approx 0.678$（mKdV），控制分岔方向；(2) 不动点方法（Krasnoselskii + Sturm 振荡理论）—— 证明正解的存在性（$k=2$ 情形，$k=3$ 仍开放）；(3) 变分方法—— 三种受限极小化（最小能量 + 动量 + 质量、能量 + 动量、二次能量 + 动量约束在三次能量和零均值上），对应约束最优化问题。最后给出谱稳定性判据（基于线性化算子 $M_c$ 的负本征值数），并显式给出 $\alpha = 1$（BO 方程）下的精确解（用双曲函数表示）。

3. 核心方程与概念

3.1 分数阶 KdV/mKdV 方程

分数阶 KdV（二次非线性，$\alpha \in (0, 2)$）： $$u_t + uu_x = (-\Delta)^{\alpha/2} u_x$$ 色散关系 $\omega(k) = c_1 k + c_{1+\alpha} |k|^\alpha k + O(|k|^{2\alpha} k)$ as $k \to 0$，对应非局部分数阶色散。

分数阶 mKdV（三次非线性）： $$u_t + u^2u_x = (-\Delta)^{\alpha/2} u_x$$

分数阶 Laplacian（Fourier 形式）在周期域 $\mathbb{T} = \mathbb{R}/(2\pi\mathbb{Z})$ 上： $$f(x) = \sum_{n\in\mathbb{Z}} f_n e^{inx}, \quad (-\Delta)^{\alpha/2} f(x) = \sum_{n\in\mathbb{Z}} |n|^\alpha f_n e^{inx}$$

守恒量（分数阶 KdV 类）：能量 $E(u) = \frac{1}{2}\int_\mathbb{T} ((-\Delta)^{\alpha/4} u)^2 dx - \frac{1}{6}\int_\mathbb{T} u^3 dx$；动量 $F(u) = \frac{1}{2}\int u^2 dx$；质量 $M(u) = \int u dx$。

能量空间：$H^{\alpha/2}(\mathbb{T})$，需要 $\alpha > 1/3$（KdV）或 $\alpha > 1/2$（mKdV）以保证非线性项的 Sobolev 嵌入。

适定性（关键定理）： - KdV：$H^s(\mathbb{R})$，$\alpha \in [1, 2]$，$s > (9-3\alpha)/4$（局部，Kenig-Ponce-Vega 1991）；全局 $L^2$ 适定性 $\alpha \in (1, 2]$（Herr-Ionescu-Kenig-Koch 2010）； - mKdV：全局 $H^{\alpha/2}(\mathbb{R})$，$\alpha \in (1, 2]$（Guo 2012）； - $\alpha = 1$（BO 方程）：临界动量 $F_0$ 处有限时间 blow-up（Martel-Pilod 2017）。

3.2 行波方程与 Galilean 变换

行波 ansatz $u(x,t) = \varphi_c(x-ct)$，$\varphi_c$ 满足 $$(-\Delta)^{\alpha/2}\varphi_c + c\varphi_c + b = \frac{1}{k}\varphi_c^k$$ $k = 2$（KdV）或 $3$（mKdV），$b$ 是积分常数。

Galilean 变换（消去 $b$，$k=2$ 情形）： $$\varphi_{c,b} - \sqrt{c^2 + 2b} = \phi_\omega, \quad \omega := \sqrt{c^2 + 2b}$$ 则 $\phi_\omega$ 满足 $$(-\Delta)^{\alpha/2}\phi_\omega - \omega\phi_\omega = \frac{1}{2}\phi_\omega^2$$ 重要结论：当 $c^2 + 2b < 0$ 时无周期解（因为 $\int(-\Delta)^{\alpha/2}\phi_c dx = -\int(b + c\phi_c - \phi_c^2/2)dx > 0$ 不可能为 0）。

对 $k=3$（mKdV），没有 Galilean 变换——$b$ 不可消去，需用两参数变分法。

3.3 小振幅展开

设 $\phi_\omega = a\phi_1 + a^2\phi_2 + a^3\phi_3 + O(a^4)$，$\omega = \omega_0 + a^2\omega_2 + O(a^4)$，$\omega_0 = 1$（$n=\pm1$ 处零本征值）。

逐步求解（$\phi_1 = \cos x$ 归一化）： - $O(a)$：$((-\Delta)^{\alpha/2} - 1)\phi_1 = 0$ —— $\phi_1 = \cos x$ ✓ - $O(a^2)$：$\phi_2 = -\frac{1}{4} + \frac{1}{4(2^\alpha - 1)}\cos(2x)$ - $O(a^3)$：$\omega_2 = \frac{1}{4} - \frac{1}{8(2^\alpha - 1)}$，$\phi_3 = \frac{\cos(3x)}{8(2^\alpha - 1)(3^\alpha - 1)}$

小振幅周期解（$k=2$，$b=0$）： $$\varphi_c(x) = 2c + a\cos(x) - \frac{a^2}{4} + \frac{a^2}{4(2^\alpha - 1)}\cos(2x) + \frac{a^3\cos(3x)}{8(2^\alpha - 1)(3^\alpha - 1)} + O(a^4)$$ $$c = 1 + \frac{2^{\alpha+1} - 3}{8(2^\alpha - 1)}a^2 + O(a^4)$$

关键阈值（来自 [X1]）： $$\alpha_0 = \frac{\log 3}{\log 2} - 1 \approx 0.585$$ $\alpha > \alpha_0$ 时 $c > 1$（解从 $\varphi = 2c$ 分岔到 $c > 1$ 侧），$\alpha < \alpha_0$ 时 $c < 1$——$\alpha < \alpha_0$ 存在 fold 分岔。

mKdV（$k=3$）阈值： $$\alpha_0 = \frac{\log 8 - \log 5}{\log 2} - 1 \approx 0.678$$ $m = 0$ 情形有 pitchfork 分岔（三种解共存）；$m \ne 0$ 时 pitchfork 破裂为 fold。

3.4 不动点方法

Green 函数（$(-\Delta)^{\alpha/2} + c$ 的逆）： $$G_c(x) = \frac{1}{2\pi}\sum_{n\in\mathbb{Z}} \frac{e^{inx}}{c + |n|^\alpha}$$ 关键性质（$G_c$ 偶函数，严格正，$x=0$ 处最大）：$\alpha > 1$ 时 $G_c(0) < \infty$；$\alpha \le 1$ 时 $G_c(0) = \infty$。

Krasnoselskii 不动点定理（[X2] Theorem 1.1）—— 正解存在性：定义正集 $$P_c := \left\{\varphi \in L^2(\mathbb{T}) : \varphi(x) \ge \frac{m_c}{M_c}||\varphi||_{L^2}\right\}$$ 其中 $m_c = \min_x G_c(x)$，$M_c = ||G_c||_{L^2}$。算子 $\mathcal{A}_c(\varphi) = \frac{1}{2}G_c * \varphi^2$ 在 $P_c$ 上有紧、自映射性，存在 $r_-, r_+ > 0$ 使得 $\mathcal{A}_c$ 在 $P_c \cap (B_{r_+} \setminus B_{r_-})$ 上有不动点。

Theorem 2（正解存在性）：固定 $\alpha \in (\alpha_0, 2]$，对每个 $c > 1$ 存在不动点 $\varphi_c \in H^{\alpha/2}(\mathbb{T})$，$\min_x \varphi_c(x) > 0$（前提是 $\text{Ker}(M_c) = \text{span}(\partial_x\varphi_c)$）。

关键假设——非退化性 $\text{Ker}(M_c) = \text{span}(\partial_x\varphi_c)$：意味着 $M_c := c + (-\Delta)^{\alpha/2} - \varphi_c$ 的零本征值是简单的。当 $\alpha < \alpha_0$ 时该条件被破坏（出现 fold 分岔），Theorem 2 不适用。

Sturm 振荡理论（来自 [X3]）：$M_c$ 在非退化性假设下只有一个负本征值——确保解是"正单瓣"。

$k=3$ 情形的开放性：正集 $P_c$ 的构造不直接（$\varphi^3$ 可负），解析证明仍未解决。

3.5 变分方法（三套框架）

变分问题 1（最直观的，§4.1）：约束极小化 $$\inf\{E(u) : F(u) = F_0, M(u) = M_0\}$$ 但 $E, F, M$ 与 $c, b$ 关系不直接，且 $C^1$ 平滑性在 $\alpha < \alpha_0$ 失败（fold）。

变分问题 2（[X4]，§4）：去掉 $M$ 约束，固定 $b$： $$\inf\{E(u) + bM(u) : F(u) = F_0\}$$ 对 $\alpha > 1/2$ 局部极小化存在（受孤子稳定性约束）。非退化 $\text{Ker}(M_c) = \text{span}(\partial_x\varphi_c)$ 在 $\alpha > 1/2$ 满足，$C^1$ 平滑性好。但未涵盖所有周期解。

变分问题 3（[X5]，最成功，eq. 24）：最优框架——极小化二次部分 + 动量，约束三次能量 + 零均值： $$\inf\left\{\int_\mathbb{T} \left[((-\Delta)^{\alpha/4}u)^2 + cu^2\right]dx : \int u^3 dx = 1, \int u\,dx = 0\right\}$$ Theorem 3（关键）：对每个 $\alpha \in (1/3, 2]$，$c \in (-1, \infty)$，存在全局极小化子 $\chi_c \in H^{\alpha/2}(\mathbb{T})$，具有单瓣剖面。这是本章最重要的存在性结果——对所有 $\alpha$ 一致。

参数 $b = b(c)$ 的确定： $$b(c) = \frac{1}{4\pi}\int_\mathbb{T}\varphi_c^2 dx = \frac{1}{2\pi}F(\varphi_c) > 0$$

mKdV 情形的变分（eq. 33）：双参数 $(c, m)$，$m$ 是周期解的平均值 $$\inf\left\{\int \left[((-\Delta)^{\alpha/4}u)^2 + cu^2\right]dx : \int u^4 dx = 1, \int u\,dx = 2\pi m\right\}$$ $m = 0$ 时存在 pitchfork 分岔；$m \ne 0$ 时 pitchfork 破裂为 fold。

3.6 谱稳定性

线性化方程（在移动坐标系下）： $$v_t = \partial_x M_{c,b} v$$ 其中 $M_{c,b} = (-\Delta)^{\alpha/2} + c - \varphi_{c,b}$（$k=2$ 情形）。

谱稳定性定义：若 $\partial_x M_{c,b}$ 在 $L^2(\mathbb{T})$ 上所有本征值 $\text{Re}(\lambda) \le 0$，则行波谱稳定；否则谱不稳定。

线性稳定（[X6]）：若满足三正交性 $\langle 1, v\rangle = \langle\varphi_{c,b}, v\rangle = \langle\partial_x\varphi_{c,b}, v\rangle = 0$，则 Lyapunov 泛函 $\langle M_{c,b} v, v\rangle$ 守恒 + 强制正定 → 线性稳定。

Theorem 5（稳定判据）：设 $\varphi_{c,b} \in H^\infty(\mathbb{T})$ 关于 $(c, b)$ 是 $C^1$。则 - 非退化 $\text{Ker}(M_{c,b}) = \text{span}(\partial_x\varphi_{c,b})$ ✓ - 雅可比 $\frac{\partial(F, M)}{\partial(c, b)}$ 的正本征值数 = $M_{c,b}$ 负本征值数 - 则行波线性稳定。

Theorem 6（$b = b(c)$ 情形，最简单稳定判据）：对问题 (24) 的局部极小化子 $\varphi_c$，$b'(c) > 0$ 线性稳定，$b'(c) \le 0$ 线性不稳定。

小振幅极限：$b'(c) = 2(2^\alpha - 1) + O(a^2)$，对 $\alpha \in (1/3, 2]$ 都 > 0——小振幅周期行波在所有 $\alpha$ 下都线性稳定。

mKdV 稳定判据（$k=3$）：与 $\partial_m b(c, m)$ 的符号相关——$\partial_m b < 0$ 对应简单负本征值，$\partial_m b > 0$ 对应两个负本征值。

$b'(c) = 0$ 阈值的"线性增长"失稳：行波不增长但扰动以 $O(t)$ 速度漂移（非指数不稳定）。

3.7 显式解：Benjamin-Ono 方程（$\alpha = 1$）

BO 方程周期解（$k = 2$）： $$\omega = \coth\gamma, \quad \phi_\omega(x) = \frac{2\sinh\gamma}{\cosh\gamma - \cos x} - 2\omega$$ $\gamma \in (0, \infty)$，小振幅极限 $\gamma \to \infty$ 对应 Theorem 1。

通过 Galilean 变换和零均值约束： $$\varphi_c(x) = \frac{2\sinh\gamma}{\cosh\gamma - \cos x} - 2, \quad c = \omega - 2, \quad b(c) = 2(\omega - 1)$$ 消元：$b(c) = 2(c+1)$ 对 $c \in (-1, \infty)$，$b$ 关于 $c$ 单调递增 → 由 Theorem 6 全部稳定。

3.8 数值结果概览

$b(c)$ 关系（图 4）： - $\alpha \in [1/2, 2]$：$b(c)$ 关于 $c$ 单调递增，全部稳定； - $\alpha \in (1/3, 1/2)$：$b(c)$ 非单调，$c \in (c_*, \infty)$ 处失稳。

mKdV 的三支解（图 5, 6，$\alpha = 1$）： - $m = 0$ 时存在半周期对称的 $\varphi_c(\pi - x) = -\varphi_c(x)$ 单支解； - 越过 $c_0$ 后出现 pitchfork 分岔，三支共存； - $m \ne 0$ 时 pitchfork 破裂为 fold，全局极小化子单独存在一支； - 数值验证：两支极小化子稳定，saddle 点失稳。

4. 关键结论

分数阶 KdV/mKdV 周期行波的存在性完整证明已通过 3 种互补方法（小振幅、不动点、变分）建立。最关键的存在性定理是 Theorem 3（变分问题 24 的全局极小化子对 $\alpha \in (1/3, 2]$, $c \in (-1, \infty)$ 全部存在单瓣解）。
谱稳定性完全由 $b'(c)$ 符号决定——这是分数阶行波理论的核心成果。判据形式简洁（$b'(c) > 0$ 稳定），与经典 KdV 的 $dP/dc > 0$ 判据同构。
$\alpha_0$ 阈值（$0.585$ for KdV, $0.678$ for mKdV）是 fold / pitchfork 分岔的临界值：$\alpha > \alpha_0$ 时单支稳定解普适存在，$\alpha < \alpha_0$ 时需要展开更精细的理论（多分支、saddle 点）。
Galilean 变换对 $k=2$ 情形消去了参数 $b$，使分析大幅简化；对 $k=3$ 失效，需双参数变分。
小振幅展开得到显式 Fourier 系数（$O(a^4)$ 精度），是任何数值方法的基准。
Benjamin-Ono 极限（$\alpha = 1$）的精确解是唯一可用解析表达的特殊情形，且 $b(c) = 2(c+1)$ 严格单调 → 全稳定。
变分方法 3（二次能量 + 动量约束在三次能量和零均值上）是最普适的框架——避免了 Theorem 2 的非退化性要求。

5. 挑战和开放性问题

$k=3$（mKdV）情形下正单瓣解的存在性（Theorem 2 的对应）仍未解决——$\varphi^3$ 不保持正性使得 $P_c$ 的构造不直接。
最小化子的非退化性 (29), (41) 是否普适成立——数值强烈暗示，但证明缺失。
周期解的唯一性——孤子情形（$c \to \infty$ 极限）已有唯一性（Frank-Lenzmann 2013），但周期情形仍开放。
$\alpha < 1/3$ 的适定性——分数阶 KdV 在 $H^{\alpha/2}$ 上的全局良定义性在 $\alpha < 1/3$ 时未知。
多支解的稳定性——saddle 点分支在 Theorem 5, 6 中已分析（$\partial_m b < 0$ 等），但完整的 $b(c, m)$ 曲面上稳定性分布图未绘制。
高阶非线性（quartic, quintic）的存在与稳定性——目前只处理二次和三次。
多维分数阶 KdV 的横向稳定性（[X7] 已开始研究孤子情形，周期行波情形未涉及）。
调制不稳定性（Benjamin-Feir instability）——分数阶 KdV 行波的旁支稳定性未研究。
mKdV 情形下，$b'(c)$ vs $\partial_m b(c, m)$ 的完整稳定性图谱未建立。
有限时间 blow-up 在 $\alpha = 1$（BO）临界动量 $F_0$ 处的实际可观测性——需要适定性理论在 $H^{\alpha/2}$ 中的精细刻画。

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

优点： 1. 理论性极强——是全书中唯一专注"用分析方法证明存在性和稳定性"的一章。三种方法（小振幅、不动点、变分）相互补充而非竞争，每种都贡献了独特的工具。 2. 变分问题 3（Theorem 3）的存在性定理异常漂亮——对所有 $\alpha \in (1/3, 2]$ 和所有 $c \in (-1, \infty)$ 同时给出存在性 + 单瓣性质。这比 Ch1 的"分数阶解只在某些参数下存在"严谨得多。 3. 谱稳定性判据 $b'(c)$ 的简洁性是数学之美的体现——把无限维算子的谱问题简化为标量函数的单调性。 4. 小振幅展开给出了显式 Fourier 系数，是数值方法的硬基准。 5. $\alpha_0$ 阈值的发现（$\log 3 / \log 2 - 1$）是出人意料的"自然常数"——$2^\alpha$ 的算术结构直接进入物理结论。 6. 数值结果与理论预测的对应——图 3-6 都很清晰，特别是 $b(c)$ 的非单调性预言了 $\alpha < 1/2$ 的失稳。

缺点： 1. 范围过窄——只讨论 KdV 和 mKdV 两个方程，但没有提及其他分数阶方程（如分数阶非线性 Schrödinger、分数阶 Boussinesq）。这种集中是力量也是限制。 2. 多维分数阶 KdV 的存在性/稳定性完全没有——只用了 §7 末一句"open question"提及。本章本可以拓展。 3. 变分方法 1（§4.1）作为引入给出了 $\alpha < \alpha_0$ 时 $C^1$ 平滑性失败的论述，但未深入分析 fold 后的多支解——直接跳到了变分方法 3。 4. 对分数阶 Laplacian 在 $H^{\alpha/2}$ 上的谱理论叙述过于浓缩——许多细节（如 $\alpha = 1$ 时 $G_c(0) = \infty$ 如何处理）未充分展开。 5. 数值方法只简短提到 Newton 和 Petviashvili，没有给出实际计算细节。与 Ch3, Ch6 形成对比：Ch5 是纯理论，Ch6 是纯数值。 6. 失稳机制——Theorem 5, 6 给出了稳定判据，但没有提供失稳模式（如最不稳定本征值的具体形式、增长率的渐近展开等）。 7. 与物理应用完全脱节——所有讨论都是"数学存在性 + 数学稳定性"，没有任何物理或生物学例子说明"为什么这些行波重要"。

对本人研究（血管生物力学 / SMC 介导的血管重塑）的关联： - 分数阶 KdV 的"长程相互作用"可类比血管壁层间弹性耦合——内弹力板、中膜、外膜之间的弹性传递本质上是非局部的（应力不会瞬时传递，而是通过长程相互作用）。分数阶 KdV 的 $\alpha$ 调节可对应"层间耦合强度"。 - Theorem 6 的"标量函数 $b'(c)$ 决定稳定性"思想可推广到血管壁——例如血管壁的稳定性指标可用某个"广义能量梯度"代替 $b'(c)$，判定血管是否进入失稳态（动脉瘤前兆）。 - mKdV 情形的 pitchfork 分岔对应血管壁的"双稳态重塑"——两种稳态（如正常壁 vs 增厚壁）共存，由某种参数（如血压）切换。分数阶算子调节分岔点位置。 - $\alpha_0$ 阈值的发现说明分数阶系统中存在"算术-物理"的深层对应（$2^\alpha = 3$ 处的临界），这种 $\alpha$ 阈值思想对研究血管粘弹性的幂律松弛可能有启发。 - 变分方法 3 的"约束优化"框架与本人血管生长-重塑理论（Goriely 2017, Ch17）的"约束极小化 + 拉格朗日乘子"思路同源。

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

[X1] P.G. Grinevich, P.M. Santini. —— 小振幅展开与 $\alpha_0$ 阈值发现。
[X2] M.A. Krasnoselskii. —— Krasnoselskii 不动点定理。
[X3] F. Steiner, P. Górka, et al. —— 分数阶 Sturm 振荡理论。
[X4] J. Angulo. Nonlinearity 31, 920 (2018). —— 变分方法 2（最小能量 + 约束动量）。
[X5] D.E. Pelinovsky, U. Le. —— 变分方法 3（最小二次能量 + 约束三次能量 + 零均值），Theorem 3 的来源。
[X6] V.I. Arnold. —— 经典 KdV 行波的线性稳定理论。
[X7] O. Riaño, S. Roudenko. arXiv:2210.09159 (2022). —— 多维分数阶 KdV 孤子横向稳定性。
[X8] R.L. Frank, E. Lenzmann. Acta Math. 210, 261 (2013). —— 分数阶 KdV 孤子唯一性。
[X9] Y. Martel, D. Pilod. Math. Ann. 369, 153 (2017). —— mKdV $\alpha = 1$ 的 blow-up。
[X10] C.E. Kenig, G. Ponce, L. Vega. J. Am. Math. Soc. 4, 323 (1991). —— KdV 经典适定性。
[X11] S. Herr, A.D. Ionescu, C.E. Kenig, H. Koch. Commun. PDE 35, 1827 (2010). —— 分数阶 KdV $L^2$ 全局适定性。
[X12] Z. Guo. J. Differ. Equ. 252, 2053 (2012). —— mKdV $H^{\alpha/2}$ 全局适定性。
[X13] L. Abdelouhab, J. Bona, M. Felland, J.C. Saut. Physica D 40, 360 (1989). —— 非局部 KdV 的奠基工作。
[X14] C. Klein, C. Sparber, P. Markowich. Proc. R. Soc. A 470, 20140364 (2014). —— 分数阶 Schrödinger 数值方法（Ch6 也引）。
[X15] F. Linares, D. Pilod, J.C. Saut. SIAM J. Math. Anal. 46, 1505 (2014). —— 分数阶 KdV 适定性扩展。
[X16] F. Linares, D. Pilod, J.C. Saut. Adv. Differ. Equ. 20, 853 (2015). —— 分数阶 KdV 孤子存在性。
[X17] L. Molinet, D. Pilod, S. Vento. Ann. Inst. H. Poincaré 35, 1719 (2018). —— $\alpha < 1$ 的适定性。
[X18] N. Laskin. Phys. Rev. E 66, 056108 (2002). —— 分数阶量子力学（Ch1 也引）。
[X19] A.M. Turing. Philos. Trans. Royal Soc. London, Ser. B 237, 37 (1952). —— Turing 不稳定（Ch2, Ch3 也引）。
[X20] J.D. Murray. Mathematical Biology II (Springer, 3rd ed., 2003). —— 行波在生物数学中的应用（Ch2 也引）。

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