Global existence, blow-up behavior and numerical simulations for a class of chemotaxis-driven fish-mussel systems

技术摘要：一类趋化驱动的鱼类-贻贝系统的全局存在性、爆破行为及数值模拟

问题陈述
本研究探讨了一个描述营养物质（$u$）、鱼类种群（$v$）和贻贝种群（$w$）在有界区域 $\Omega \subset \mathbb{R}^n$ 内时空动力学的趋化驱动生态系统模型。该模型被构建为一个耦合偏微分方程（PDE）系统，并受齐次诺伊曼（Neumann）边界条件约束。该系统整合了营养物质扩散、鱼类向营养丰富区域的趋化运动，以及贻贝作为滤食者在营养循环中的作用。控制方程 (1.1) 包含一个参数 $\tau \in \{0, 1\}$，用以区分抛物型-抛物型-抛物型系统（$\tau=1$）与抛物型-椭圆型-抛物型系统（$\tau=0$）。本研究旨在建立经典解的全局存在性与有界性，分析导致有限时间爆破（特别是三维情况下）的条件，并开发一种数值方案以验证这些理论发现。

方法论
作者结合了分析技术与数值模拟：

1. 分析框架：

   • 局部存在性： 利用半群方法和巴拿赫不动点定理，基于椭圆型和抛物型正则性理论，建立了局部时间内唯一经典解的存在性。
   • 先验估计： 为了证明全局存在性，作者推导了 $L^\infty$ 和 $L^p$ 估计。对于 $\tau=1$ 的情况，使用了一个专门的测试函数 $\phi(x)$（通过引理 2.5 构建），以处理趋化项 $-\chi\nabla \cdot (v\nabla u)$。这使得在特定的趋化敏感度 $\chi$ 和营养物质界限约束下，能够推导出鱼类种群 $v$ 的 $L^p$ 估计。
   • 延展性判据： 全局存在性的证明依赖于证明解的 $L^\infty$ 范数保持有界，从而防止爆破。这通过结合 $L^p$ 估计与诺伊曼热半群的平滑性质来实现。
   • 爆破分析： 对于三维情况（$n=3$），作者推导了爆破时间 $T_{max}$ 的下界。这涉及构造一个辅助能量泛函 $\Psi_\tau$（定义于 1.8）并推导其演化的微分不等式。通过对该不等式进行积分，得到了一个显式的爆破时间下界。

2. 数值方法：

   • 离散化： 开发了一种有限元方法（FEM），使用分片线性（$P1$）四面体单元进行空间离散，并采用后向欧拉格式进行时间积分。
   • 实现： 使用 PETSc 中的 SNES Newton 线搜索求解器处理全耦合非线性系统，并利用经由 LU 分解预处理的 GMRES 求解线性化系统。模拟过程在 FEniCSx 中执行，并使用 MATLAB 进行可视化。
   • 验证： 通过网格收敛性研究以及将长期动力学与通过 MATLAB ode45 求解的空间均匀 ODE 模型进行对比，验证了数值方案的有效性。

主要贡献与结果

• 全局存在性与有界性（定理 1.1）： 本文确立了该系统在维度 $n \geq 3$ 下唯一经典解的全局存在性。具体而言，若初始数据满足一定的正则性条件，且参数满足以下不等式：
  $$\max \left\{ \tau\|u_0\|_{L^\infty(\Omega)}, \frac{M}{\alpha} \right\} < \min \left\{ \frac{\pi}{\chi}\sqrt{\frac{2}{n}}, \frac{\delta_2}{\sigma_2} \right\}$$
  则解 $(u, v, w)$ 在时间上是一致有界的。该结果对于 $\tau=0$ 和 $\tau=1$ 均成立。

• 爆破时间估计（定理 1.2）： 在解发生有限时间爆破的情况下（特别是在三维凸区域内），作者推导了一个可计算的爆破时间 $T_{max}$ 下界。该界限通过涉及能量泛函 $\Psi_\tau(0)$ 的初值及系统参数的积分来表示。此外，还提供了一个显式的简化下界：
  $$T_{max} \geq \frac{1}{2E_\tau \Psi_\tau(0)^2}$$

• 数值收敛性： 所提出的有限元方案在 $L^2(\Omega)$ 和 $L^\infty(\Omega)$ 范数下均表现出二阶收敛性，这与使用线性单元的结果一致。

• 数值模拟：

  • 稳定性： 在中等参数下的模拟表明，系统收敛至共存稳态，验证了由 ODE 分析预测的长期稳定性。
  • 爆破行为： 在具有高初始浓度和特定参数选择的三维立方体区域内的数值实验，展示了鱼类种群密度 $v$ 的快速、不受控增长，证实了有限时间爆破在理论上的可能性。模拟显示，尽管 $v$ 发生爆破，但营养物质 $u$ 和贻贝 $w$ 的密度保持有界，这与推导出的理论界限相符。

意义与主张
本文声称为理解涉及趋化鱼类和滤食贻贝的养殖生态系统的稳定性及潜在不稳定性提供了严谨的数学基础。通过将以往基于 ODE 的模型扩展到 PDE 框架，本研究强调了空间异质性和定向运动（趋化）在生态系统动力学中的关键作用。

作者将其工作定位为理论分析与计算建模之间的桥梁。他们证明了虽然在特定参数机制下可以保证全局存在性（从而防止由富营养化引起的崩溃），但如果违反这些条件，系统在更高维度下极易发生有限时间爆破。所开发的鲁棒有限元方案允许可视化这些复杂的动力学过程，为验证理论爆破判据以及模拟从稳定共存到种群爆发的转变提供了工具。本研究并非提出新的生态管理策略，而是提供了用于评估此类系统在不同环境条件下稳定性的分析与数值工具。