《Mathematics》：The Variation Iteration Method Combined with the Natural Generalized Laplace Transform for Solving Fractional Burgers Equations
编辑推荐：
本研究采用新的自然广义拉普拉斯变换（Natural Generalized Laplace Transform, NGLT）与双自然广义拉普拉斯变换（Double Natural Generalized Laplace Transform, DNGLT）并结合
本研究采用新的自然广义拉普拉斯变换（Natural Generalized Laplace Transform, NGLT）与双自然广义拉普拉斯变换（Double Natural Generalized Laplace Transform, DNGLT）并结合变分迭代法（Variational Iteration Method, VIM）求解分数阶 Burgers 方程。首先，给出了自然变换、广义拉普拉斯变换以及 Caputo 分数阶导数的基本定义，为本研究提供理论基础。本文主要关注自然广义拉普拉斯变分迭代法（Natural Generalized Laplace Variational Iteration Method, NGLTVIM），该方法是求解常规问题的一种新途径。研究人员对所提方法的稳定性与收敛性进行了详细讨论，以证明其可靠性与有效性。通过若干数值算例，展示了该方法在求解一维分数阶耦合 Burgers 方程和奇异分数阶耦合 Burgers 方程中的有效性。结果表明，NGLTVIM 可成功应用于求解多种数学物理问题。
该论文发表于《Mathematics》，围绕分数阶 Burgers 方程的解析—半解析求解问题，构建了自然广义拉普拉斯变换（Natural Generalized Laplace Transform, NGLT）/双自然广义拉普拉斯变换（Double Natural Generalized Laplace Transform, DNGLT）与变分迭代法（Variational Iteration Method, VIM）耦合的新方法框架，并系统讨论了其稳定性与收敛性。研究背景在于，分数阶 Burgers 方程是流体动力学与扩散系统中的重要数学模型，其在经典 Burgers 方程中引入 Caputo 分数阶导数，用于刻画反常扩散与非局域效应，这类效应广泛存在于湍流、黏弹性材料等物理场景。现有方法如自然变换分解法、Adomian 分解法、拉普拉斯变分迭代法及若干推广形式，虽已用于多类偏微分方程与分数阶偏微分方程，但在处理非线性、耦合以及奇异型分数阶 Burgers 系统时，仍面临计算复杂、构造过程繁琐或统一性不足的问题。因此，有必要发展一种既保留变换法处理初值条件的优势，又能通过迭代机制有效处理非线性项，并具有严格稳定性与收敛性保证的新方法。
基于上述问题，研究人员提出自然广义拉普拉斯变分迭代法（NGLTVIM）以及相应的双重形式（DNGLTVIM），用于求解时间分数阶 Burgers 方程及其耦合系统。论文首先回顾自然变换、广义拉普拉斯变换和 Liouville–Caputo 偏导数的定义，并以 Banach 不动点定理和 Picard Λ-稳定性理论作为后续分析基础。随后，研究人员系统给出 NGLT 与 DNGLT 的定义、逆变换及若干基本性质，特别推导了 Caputo 分数阶偏导数在 DNGLT 下的变换公式，以及若干一阶、二阶偏导项的双自然广义拉普拉斯变换表达式。这些理论结果为将含分数阶导数和空间导数的 Burgers 方程转化为便于迭代求解的代数—积分形式奠定了基础。
方法上，研究人员针对二维 Burgers 方程建立了 DNGLTVIM 的一般步骤：先施加双自然广义拉普拉斯变换，再处理初始条件，继而取逆变换，将问题转化为关于辅助变量的关系式，对其求导后引入 VIM 构造校正泛函，最后得到可递推的迭代公式。通过变分理论确定拉格朗日乘子后，建立了新的迭代结构，并将解表示为收敛级数。进一步地，论文定义了相应算子并据此展开稳定性与收敛性分析。对于一般 Burgers 方程，研究人员还给出单重 NGLT 与双重 DNGLT 联合使用的求解框架，并结合 Adomian 多项式（用于展开非线性项）处理方程中的非线性耦合项，从而形成适用于一维耦合系统及奇异系统的统一算法流程。
主要技术方法概括如下：研究首先构建自然广义拉普拉斯变换（NGLT）与双自然广义拉普拉斯变换（DNGLT）的理论体系，推导分数阶导数及偏导项的变换性质；其次将变分迭代法（VIM）嵌入变换域与逆变换后的递推框架，通过拉格朗日乘子的变分确定获得迭代校正公式；再次利用 Banach 不动点定理与 Picard Λ-稳定性分析证明所建算子的稳定性，并在 Banach 空间中证明级数解的收敛性；最后通过两个数值算例验证方法对一维分数阶耦合 Burgers 方程与奇异分数阶耦合 Burgers 方程的适用性。论文未涉及样本队列来源。
以下结合论文各部分小标题，对主体内容作出凝练解读。
一、Introduction
引言部分指出，分数阶 Burgers 方程在流体流动演化、扩散系统建模中具有重要意义，Caputo 分数阶导数能够表征反常扩散和非局域行为。作者综述了自然变换分解法、Laplace 变分迭代法、同伦摄动法、微分变换法以及其他分数阶方程数值与解析方法在相关问题中的应用，说明现有研究已覆盖多类偏微分方程，但对于分数阶 Burgers 方程，尤其是耦合与奇异情形，仍有改进空间。基于此，研究人员明确提出将 VIM 与 NGLT 相结合，构建自然广义拉普拉斯变分迭代法（NGLVIM/NGLTVIM）来求解分数阶 Burgers 方程。
二、Properties of Natural Generalized Laplace Transform
本节建立方法论基础。研究人员定义了 NGLT 与 DNGLT，并给出逆变换形式；随后通过参数特化说明该变换框架可退化为双 Sumudu 变换、双 Laplace 变换以及 Laplace–Yang 变换，表明该方法具有较好的统一性。文中还给出若干基本函数在 NGLT 下的变换结果，以及 Gamma 函数参与的积分表达。更重要的是，作者通过定理推导了 Caputo 分数阶偏导数的 DNGLT 公式，以及若干空间偏导项在 DNGLT 下的变换表达式。这些结果说明，原方程中的分数阶时间导数和空间导数经变换后能够被系统地处理，从而为后续建立迭代格式提供了严格依据。
三、New Double Natural Generalized Laplace Variational Iteration Methods for Solving the 2D Burgers Equation
本节给出 DNGLTVIM 求解二维 Burgers 方程的核心构造。研究人员按照“变换—代入初值—逆变换—求导—变分迭代”的步骤，先将时间分数阶 Burgers 方程与初始条件变换到 DNGLT 域，再经逆变换得到包含源项和初值项的表达式；随后对辅助变量求导，并应用 VIM 构造递推关系。通过驻值条件确定拉格朗日乘子后，得到新的迭代公式。文中进一步指出，该迭代公式能够生成收敛级数解，且各项便于计算。研究人员还定义了相应算子及其分量形式，为后续稳定性与收敛性分析建立抽象框架。该部分的主要结论是：DNGLTVIM 能够将二维分数阶 Burgers 方程转化为可递推、可收敛的级数求解问题。
四、Examination of the Stability and Convergence of the NGLTVIM Method
本节从理论上验证方法的可靠性。研究人员将前述迭代过程写成 Banach 空间中的自映射问题，并证明该映射满足 Picard 稳定性条件。Theorem 6 表明，在给定线性算子与非线性算子满足相应有界性和压缩性条件时，所定义的自映射具有不动点，从而方法是 Picard Λ-稳定的。随后，Theorem 7 证明由迭代关系生成的部分和序列构成 Cauchy 序列，因此在 Banach 空间中收敛，进而相应级数解收敛到原问题的解。该节的核心贡献在于，论文不仅提出算法，还给出稳定性与收敛性的严格论证，增强了方法的数学可信度。
五、New Natural Generalized Laplace Variational Iteration Methods for Solving the Burgers Equation
在这一部分，研究人员将方法推广到更一般的 Burgers 方程求解。对于给定初始条件的时间分数阶 Burgers 方程，作者分别对控制方程和初值条件施加 DNGLT 与自然变换（NT），得到变换域中的代数关系；经逆变换和求导后，再利用 VIM 建立递推公式。文中还引入 Adomian 多项式来展开非线性项，如非线性乘积项及其耦合形式，从而使迭代公式能够逐项构造近似解。该节说明，NGLTVIM 不仅适用于二维情形，也适用于更一般的一维与耦合型 Burgers 系统，并具有统一处理线性项、非线性项与初始条件的优势。
六、Numerical Examples
本节通过两个算例验证方法的有效性。
Example 1
该例研究一维时间分数阶耦合 Burgers 方程。研究人员先对方程施加 NGLT，接着进行逆变换并对辅助变量求导，再依据 VIM 构造关于未知函数的递推关系；非线性项由前文定义的 Adomian 多项式处理。文中计算了前若干次迭代近似，并给出级数形式的近似解。进一步地，作者指出当分数阶参数 β=1 时，可恢复精确解。通过该例可见，所构造的解序列具有收敛与稳定特征。数值展示部分说明，不同分数阶阶次会影响解的演化轮廓，而当 β 取 1 时数值解与精确解一致，从而验证了方法的准确性。
Example 2
该例讨论奇异分数阶耦合 Burgers 方程组。由于方程含有奇异项，研究人员首先通过乘以 x 消除奇异性，再施加 NGLT 与 NT，借助前述定理将方程写成适合逆变换和迭代处理的形式。随后，经逆变换、求导和 VIM 构造校正泛函，得到递推关系，并逐步写出前几阶近似解。文中进一步给出一般第 n 阶近似解的表示，并说明当 n→∞ 时，所得级数收敛到该初值问题的精确解；当 β→1 时，精确解退化为经典情形下的结果。该例表明，所提方法不仅能处理标准分数阶 Burgers 方程，也能有效处理带奇异性的耦合系统。
七、Conclusions
结论部分指出，论文提出的自然广义拉普拉斯变换（NGLT）与双自然广义拉普拉斯变换（DNGLT）可作为求解复杂分数阶问题，尤其是分数阶 Burgers 方程的有效新技术。研究人员给出了相关定义、定理与实例，展示了这些方法的优势。所提方法通过将自然广义拉普拉斯变换与变分迭代法相结合，建立了新的解析求解方案，能够克服求解奇异分数阶微分方程时遇到的困难。数值算例证实，所提方法（DNGLTVIM）具有良好的精度、可靠性与高效性；所得结果表明，该方法可快速产生收敛级数解，并能高效应用于线性与非线性耦合系统。作者同时指出，未来可将该方法推广至更高维奇异偏微分方程系统及其他类型的非线性分数阶模型。
总体来看，该研究的主要意义在于：一方面，从理论上完善了 NGLT/DNGLT 作用于分数阶导数与偏导算子的性质，并建立了相应的不动点稳定性与收敛性证明；另一方面，从应用上展示了该方法对一维耦合、二维以及奇异分数阶 Burgers 方程的统一求解能力。论文的贡献集中体现在方法的新颖组合性、理论证明的完整性以及对典型数学物理模型的可操作性验证。
