《Scientific Reports》：Quadrature solution of fractional coupled burgers and plankton-oxygen dynamics under climate change
编辑推荐：
本研究旨在为数学建模与环境科学中的两个重要系统提供精确且计算高效的解：(1+1)维分数阶非线性耦合Burgers系统和分数阶耦合浮游生物-氧气模型动力学。尽管这些系统的复杂性不断增加，但现有的数值方法通常难以应对非局部算子的高计算成本。为了解决这一问题，研究人
本研究旨在为数学建模与环境科学中的两个重要系统提供精确且计算高效的解：(1+1)维分数阶非线性耦合Burgers系统和分数阶耦合浮游生物-氧气模型动力学。尽管这些系统的复杂性不断增加，但现有的数值方法通常难以应对非局部算子的高计算成本。为了解决这一问题，研究人员提出了一种鲁棒的混合框架，将分数阶微分求积法（FDQM）与Newton-Raphson（NR）迭代过程相结合。这些系统对于理解各种物理过程至关重要，包括湍流流体动力学以及水生环境中氧气与浮游生物的生物相互作用。这些模型使用了分数阶导数，其在捕获记忆和遗传特征方面提供了更大的灵活性，使其比传统的整数阶模型更适合正确模拟真实世界的过程。提出的工作使用混合了微分求积技术（DQM）版本的广义Liouville-Caputo分数阶模型来解决这些问题，该技术允许有效处理复杂的边界条件和空间导数。方程中的非线性通过Newton-Raphson迭代方法处理，这确保了解的稳定性和收敛性。所提出的系统在MATLAB中实现，并进行了完整的参数分析，以研究包括分数阶导数、氧气生成速率和浮游植物最大人均增长率在内的各种参数如何影响模型输出。本研究不仅证明了所建议技术的准确性、收敛性和效率，而且揭示了它们对重要参数的敏感性，使其更适用于真实环境情况。预计这项工作的结果将有助于为复杂系统建立更好的建模方法，最终使从环境科学到流体动力学领域的科研人员受益。
在数学建模与环境科学领域，分数阶偏微分方程被广泛应用于描述复杂流体力学与生态系统动力学过程。耦合Burgers系统在等离子体物理、流体力学及化学物理等众多学科中具有重要的应用价值。与此同时，全球变暖对海洋水温产生了深刻影响，水体温度作为主导变量之一，深刻影响着水中的氧气浓度，进而改变水生环境中初级物质的生产。浮游植物通过光合作用产生氧气，而浮游植物与浮游动物的捕食-被捕食动力学传统上往往缺乏对氧气浓度动态变化的充分考虑。由于分数阶算子具有非局部性，传统的数值方法在求解此类分数阶非线性系统时通常面临极高的计算成本和内存需求，难以保证解的稳定收敛。为解决这一计算瓶颈并更准确地模拟具有记忆和遗传效应的真实物理过程，研究人员开展了一项关于分数阶非线性系统数值求解的研究。该研究提出了一种将分数阶微分求积法（FDQM）与Newton-Raphson（NR）迭代技术相结合的鲁棒混合计算框架。通过引入广义Liouville-Caputo类型的分数阶导数，并利用DQM的全局加权系数矩阵，该框架能在使用显著较少节点数(N)的情况下实现高阶精度并有效处理非线性项。这项研究不仅验证了所提计算框架在处理复杂边界与高非线性耦合时的准确性与收敛性，还揭示了模型输出对关键参数的敏感性，为理解湍流动力学与水生生态系统的物理及生物相互作用提供了更真实的建模手段。该研究成果已发表在《Scientific Reports》上，为从环境科学到流体动力学等多个领域的复杂系统建模方法论发展做出了重要贡献。
研究人员为开展研究用到的主要关键技术方法包括：将分数阶微分求积法（FDQM）与Newton-Raphson（NR）迭代过程相结合构建混合数值求解框架。在空间与时间离散化过程中，采用了基于Lagrange插值多项式的多项式型微分求积法（PDQM）以及基于正则化Shannon核的离散奇异卷积微分求积法（DSCDQM-RSK）来近似各阶导数。模型采用了广义Liouville-Caputo类型的分数阶导数以捕获系统的记忆效应，并推导了对应的分数阶加权系数矩阵。针对方程的非线性特征，引入了计算Jacobian矩阵的Newton-Raphson迭代法以确保数值解的稳定收敛。研究在MATLAB环境中编写计算程序，通过对比已知解析解计算L∞误差来验证计算效率与精度，并开展不同参数下的数值模拟分析。
研究结果部分主要包含以下两个模型的数值验证与分析：
第一个问题（耦合Burgers系统）的研究结果：研究人员将微分求积法应用于耦合Burgers系统方程，得到非线性代数方程组，随后采用Newton-Raphson方法进行迭代求解。通过计算L∞误差对比了均匀与非均匀PDQM以及DSCDQM-RSK的性能。结果表明，在网格大小为20×20、时间参数为1至3秒的条件下，均匀PDQM的L∞误差在2.74×10
-4至5.24×10
-4之间；而使用基于Chebyshev-Gauss-Lobatto分布的非均匀网格，在15×15的网格下，误差显著降低至1.84×10
-4至2.99×10
-5，证明了非均匀PDQM比均匀方法及DSCDQM-RSK具有更高的精度和计算效率。研究还测试了不同网格大小、时间步长以及分数阶导数阶数α和β对结果的影响，发现W(x,t)和V(x,t)的数值解与解析解高度吻合，且随着α和β值的减小，计算结果略有下降。最终在21×21的网格下，α=1, β=1时，计算精度达到了1.75×10
-7，计算耗时仅为0.2773秒，充分证实了方法的稳定与高效。
第二个问题（分数阶浮游生物-氧气动力学模型）的研究结果：研究人员将该数值框架应用于描述气候变化背景下氧气(C)、浮游植物(U)与浮游动物(V)相互作用的非线性偏微分方程系统。研究将该系统分为三种情形进行仿真：情形1（时间动力学，即忽略扩散系数且分数阶参数φ和ρ取整）、情形2（分数阶导数作用下的模型动力学）以及情形3（空间与时间双向动力学，即考虑扩散且空间分数阶导数ε在1到2之间）。通过非均匀PDQM和DSCDQM-RSK的对比验证，发现结果在不同网格大小及氧气生产率(P)下均能良好收敛。参数分析显示，在情形1中，C(t)、U(t)和V(t)的动态变化与时间呈反比关系，直至100秒后趋于收敛；浮游植物捕获率参数(Q)的变化不会显著影响氧气浓度C(x,t)，但U(x,t)和V(x,t)的密度与Q值呈正比关系。在情形3中，C、U和V的分布与空间距离呈反比，直至x=80后趋于收敛。
总结讨论部分，研究人员翻译并总结研究结论如下：本研究成功地将基于Lagrange插值多项式与正则化Shannon核的微分求积技术与广义Liouville-Caputo型分数阶导数相结合，有效求解了具非线性的时间分数阶耦合Burgers系统及(1+1)维浮游生物-氧气动力学模型。通过与精确解的对比，验证了所提算法不仅具有极高的精确性和收敛性，同时在计算效率上表现优异。参数研究揭示了模型随时间和空间的动态演变规律，以及系统对关键生态参数的敏感性响应。研究发现，氧气浓度与氧气生成率及浮游植物的最大人均增长率呈正比，浮游动物密度也随浮游植物最大人均增长率呈正比增加。本研究的核心目标是揭示在气候变暖导致氧气匮乏等生态灾难潜在发生机制中的数学表征。通过调整参数Q的数值，模拟了氧气-浮游生物系统动态结构在气候改变过程中的演变轨迹。该工作不仅反映了当前系统在环境变化下的数学响应，也为未来在多维系统中引入风速、温度波动等更多复杂环境变量的综合研究奠定了坚实的数值基础。
