Nonlinear Boundary Dynamics in Reaction-Diffusion-Advection Systems with Nonlocal Delay: Bifurcation and Stability

技术摘要：具有非局部延迟的反应-扩散-平流系统中的非线性边界动力学

问题陈述
本文研究了一个定义在有界区域 $\Omega \subset \mathbb{R}^N$ 上的两物种 Lotka-Volterra 竞争模型的分岔与稳定性分析。该系统同时整合了三个具有数学挑战性的特征：

1. 平流（Advection）： 物种沿环境梯度（由势函数 $A_i(x)$ 和速度 $\eta$ 表示）的有向运动。
2. 非局部时间延迟（Nonlocal Time Delay）： 涉及空间卷积积分和成熟延迟 $\tau$ 的竞争项，模拟了幼体在成为成体进行竞争之前可能先进行扩散的过程。
3. 非线性边界条件： 生境边缘的密度依赖型迁出，通过 $\partial u_i / \partial n = \lambda u_i g_i(x, u_1, u_2)$ 进行建模，其中 $g_i$ 捕捉了净迁出率。

虽然这类系统的单个组成部分（例如带有线性边界的平流，或不含平流的非局部延迟）已有研究，但将这三者结合在一起的系统性分析框架仍付阙如。主要的挑战在于平流导致的自伴性丧失，以及非线性边界引入的代数复杂性，这使得无法直接应用标准的变分方法或中心流形约化。

方法论
作者开发了一个统一的分析框架，其核心在于两项主要的技术创新：

1. 加权变换与 Lyapunov-Schmidt 约化：

   • 为了处理平流，系统通过 $\hat{u}_i = e^{-\eta A_i} u_i$ 进行变换。这通过在加权 Sobolev 空间 $L^2(\Omega, e^{\eta A_i} dx)$ 中将平流-扩散算子转化为自伴形式，从而实现了主特征值的变分表征。
   • 对于稳态分岔，作者采用了 Lyapunov-Schmidt 约化而非中心流形理论。选择这一方法的原因是希望避免求解用于计算中心流形系数的偏微分方程。
   • 该方法将无穷维系统投影到线性化算子的有限维核（kernel）上。
     • 在非退化情况（$\lambda_1 \neq \lambda_2$）下，核是一维的，Crandall-Rabinowitz 定理确立了半平凡分支。
     • 在退化情况（$\lambda_1 = \lambda_2$）下，核是二维的。约化过程产生了一个有限维代数分岔方程 $G_j(z_1, z_2, \lambda) = 0$，其系数是涉及特征函数和核的显式积分。

2. 伴随延迟生成器的双线性型：

   • 为了分析稳定性和 Hopf 分岔，作者为线性化延迟生成器 $A_{\tau, \lambda}$ 构建了一个特定的双线性型 $\ll \cdot, \cdot \gg$。
   • 该形式推广了 Hale-Verduyn Lunel 双线性型（最初针对常微分方程），将其扩展到具有非局部延迟和平流的偏微分方程。
   • 至关重要的是，该构造在延迟区间内的双重积分中纳入了平流诱导的权重 $e^{\eta A_i}$，从而产生了一个指数因子 $\exp\{\eta A_k(\frac{\theta+\tau}{\tau}x - \frac{\theta}{\tau}y)\}$。这使得显式计算伴随算子并验证横截性条件成为可能。

主要结果

1. 稳态分岔：

   • 不同特征值 ($\lambda_1 \neq \lambda_2$)： 半平凡分支（其中一个物种持续存在而另一个灭绝）在主特征值处从平凡解处分岔。分岔方向由常数 $\kappa_{jj}$ 和 $w_{jj}$ 决定。
   • 重合特征值 ($\lambda_1 = \lambda_2$)： 一个共存分支出现。该分岔的性质（鞍节点分岔还是双鞍节点分岔）由从约化方程中导出的代数不变量决定，特别是 $\kappa_2 p_{11} - \kappa_1 p_{21}$ 的符号及相关的二次项。

2. Hopf 分岔判据：

   • 本文推导出了存在纯虚特征值（Hopf 分岔的前兆）的显式条件：
     $$\kappa_{11}^2(\eta) \geq w_{11}^2(\eta)$$
     此处，$\kappa_{11}$ 衡量内部非局部竞争，$w_{11}$ 衡量边界自我调节（由边界函数 $g$ 的导数导出）。
   • 判据的意义：
     • 在线性边界情形下（$w_{11}=0$），该条件始终成立，这意味着对于合适的延迟，Hopf 分岔是必然发生的。
     • 在非线性边界情形下，足够强的边界自我调节（$w_{11}^2 > \kappa_{11}^2$）可以完全抑制 Hopf 分岔，使系统在所有延迟值 $\tau \geq 0$ 下保持稳定。
   • 平流速度 $\eta$ 通过影响特征函数以及常数 $\kappa_{11}$ 和 $w_{11}$ 来改变临界延迟阈值 $\tau_0$。

3. 共存稳定性：

   • 对于退化情况（$\lambda_1 = \lambda_2$），共存分支的稳定性受两个代数不变量 $\Upsilon_1$ 和 $\Upsilon_2$ 控制。
   • 如果 $\Upsilon_1 > 0$ 且 $\Upsilon_2 > 0$，系统会在临界延迟 $\tau^*_z$ 处发生超临界 Hopf 分岔，从稳定的共存状态转变为振荡共存状态。

4. 对齐次核的应用：

   • 该理论应用于具有常系数和齐次核的一维河流模型。
   • 数值模拟证实了：
     • 平流将种群密度向下游移动。
     • 强边界自我调节（高迁出率）可以防止延迟诱导的振荡。
     • 当内部竞争占主导时，增加延迟会导致极限环（振荡共存）。

意义与主张
作者声称这项工作代表了首次针对同时包含非线性边界条件、非局部时间延迟和平流的反应-扩散-平流系统进行的系统性分岔与稳定性分析。

• 方法论贡献： 本文提供了一个统一的 Lyapunov-Schmidt 约化框架，该框架避免了使用复杂的泛函分析工具进行中心流形约化，从而提供了一种更直接的针对低维核的代数表述。
• 理论洞察： 一个主要贡献是发现非线性边界条件从根本上改变了稳定性格局。与线性边界问题中 Hopf 分岔是普遍现象不同，非线性边界可以作为一种稳定机制，如果边界自我调节足够强，它可以完全抑制振荡。
• 普适性： 该框架被呈现为适用于除种群生态学之外的其他具有相似数学结构的领域，如具有非局部连接性的流行病传播、神经场方程和趋化系统。

本文得出结论，内部竞争、边界自我调节以及平流之间的相互作用决定了种群是持续存在、灭绝还是表现出振荡动力学，并提供了可在受控系统中进行测试的显式判据。