一个广义的二阶保正性数值方法用于非自治动力系统及其应用

《Applied Mathematics and Computation》：A generalized second-order positivity-preserving numerical method for non-autonomous dynamical systems with applications

　　

在这项工作中，研究人员提出了一种广义的二阶非标准有限差分（NSFD）方法，用于由常微分方程（ODEs）定义的非自治动力系统。该方法将NSFD框架与右端函数的一种新的非局部近似相结合。该方法实现了二阶收敛，并且无条件地保持所有步长下解的正性。特别地，它避免了许多现有保正性二阶NSFD方法所要求的一个严格且不可或缺的条件。更准确地说，文献中常见的条件——ODEs的右端函数在任何时间步都不为零——被放宽了。因此，该方法易于实现且计算高效。数值实验，包括对SIR传染病模型的一个改进的NSFD格式，证实了理论结果。此外，研究人员展示了该方法对具有正解的非线性偏微分方程和边值问题的适用性，展现了其在实际建模中的多功能性。

在许多实际科学和工程问题中，非自治动力系统（由常微分方程（ODEs）定义）的解的正性是最常见且最重要的特征之一。现有的一阶非标准有限差分（NSFD）方法（由Mickens在1980年代引入）虽然能无条件保正，但通常精度不足。近年来，学者们发展了多种二阶保正性NSFD方法（如Alalhareth等提出的SOPESN方法、Conte等提出的α-NSFD方法等），这些方法在理论上是有效的，但共同面临一个严格且不可或缺的条件——右端函数在任何离散时间步都不为零（条件(2.10)）。该条件在现实中难以满足，尤其当数值解接近动力系统的零倾线时。因此，研究人员Manh Tuan Hoang和Matthias Ehrhardt开展了本研究，旨在构造一种无需此限制性条件的广义二阶保正性数值方法，以提升计算适用性和效率，并展示其在SIR传染病模型、非线性偏微分方程（PDEs）和边值问题（BVPs）中的实际应用价值。论文发表在《Applied Mathematics and Computation》。

研究人员采用的关键技术方法包括：基于表示定理将系统（1.1）转化为（1.5）形式（分解为f_i和g_i函数）；提出包含校正项A_i和B_i的新非局部近似NSFD格式（3.1）；通过泰勒展开推导二阶收敛条件（定理2），其中分母函数?_i需满足（3.5），校正项（A_i, B_i）取自集合D(v_i/(2κ_i))（v_i由（3.4）定义）；选用指数形式（3.9）或多项式形式（3.10）等分母函数；针对SIR传染病模型（4.1）具体构造了6种NSFD格式（2ndPNSFD1至2ndPNSFD6），并给出对应的?_i、φ_i选择（表1和表2）。未涉及具体试剂或培养操作，样本队列来源于经典的SIR流行病模型。

研究结果如下：
- “2. NSFD schemes based on well-known approaches”：通过分析现有二阶NSFD方法（2.1）、（2.5）和（2.8）的构造与收敛条件（引理1-3），明确指出它们均要求条件（2.10）（即F_i(t_k, y_k) ≠ 0），并通过一个二维示例（y₁' = t²(y₁-a)²+t⁴(y₂-b)⁴等）说明该条件在实际应用中的限制性。
- “3. Construction of new second-order positivity-preserving NSFD method”：提出新NSFD方法（3.1），包含校正项φ_i(A_i-y_i^k+1B_i/y_i^k)。定理1证明只要A_i, B_i ≥ 0，方法无条件保正。定理2给出二阶收敛条件：分母函数满足（3.5）且（A_i, B_i）∈ D(v_i/(2κ_i))，无需条件（2.10）。同时给出具体分母函数候选（3.9）-（3.11）。
- “4. Numerical simulation of an SIR epidemic model”：以可变系数的SIR模型（4.3）为测试问题（总人口守恒，仅考虑前两方程）。研究人员选择b(t)=1/(1+t), c(t)=2/(2+t)（精确解已知），构造6种NSFD格式（2ndPNSFD1-6），并与一阶NSFD格式（4.2）对比。表3-6显示所有二阶格式的全局误差随Δt²衰减（收敛阶接近2），而一阶格式仅呈一阶收敛。图1-3展示在Δt=2.0的大步长下，二阶格式仍保持正性并准确模拟动力学行为，而显式欧拉和梯形法则产生负解。
- “5. Second-order positivity-preserving NSFD method applied to nonlinear PDEs”：考虑一类非线性PDE（5.1）及其正性条件（定理3）。通过方法线（MOL）将PDE离散为ODE系统（5.10）或（5.13），再应用新NSFD方法。以PDE模型（5.14）为例（精确解u=e^-(x+t)），数值实验（表7）表明二阶PNSFD格式收敛阶接近2，一阶PNSFD格式收敛阶接近1，且二阶格式误差更小（图4-6）。
- “6. Second-order positivity-preserving NSFD method applied to BVPs”：考虑Bratu方程等非线性BVP（6.1）。利用解关于中点对称性，将BVP转化为初值问题（6.2），再改写为（1.5）形式（6.3），应用新NSFD方法求解。该方法通过求解一系列ODE问题得到正解，并借助射影法确定初始斜率。

在讨论部分，研究人员指出新方法可进一步推广为带可控参数α_i, β_i, γ_i的格式（7.1），以优化误差或保持其他动力学性质。同时提醒了方法线（MOL）在处理时变边界条件时可能出现的“刚性阶降”现象（如RK4退化为低阶），但本文数值示例中未观察到，将在未来工作中严格分析。研究结论部分翻译为：在本工作中，研究人员提出了一种简单且高效的方法，用于构造非自治动力系统的广义二阶保正性数值方法。该方法基于新的右端函数非局部近似与分母函数归一化。值得注意的是，所构造的方法避免了若干著名二阶保正性NSFD方法所要求的严格且不可或缺的条件，因此计算实现直接。还提供了所构造NSFD方法的重要应用，并通过数值实验支持并阐释了理论结果。因此，改进了文献[28]中针对SIR流行模型的NSFD格式。此外，介绍了并分析了所构造的二阶保正性NSFD方法在解实际中出现的一些PDE和BVP中的应用。