[论文解读] Direct Methods and Symbolic Software for Conservation Laws of Nonlinear Equations

Willy Hereman, Paul Adams+3 more|ArXiv.org|Mar 1, 2008

Nonlinear Waves and Solitons参考文献 82被引用 36

一句话总结

本文提出用于计算非线性偏微分方程（PDEs）和微分-差分方程（DDEs）中守恒律的符号算法与软件工具。利用伸缩对称性，该方法通过待定系数法构造候选密度，借助变分法与同伦算子求解线性系统，并在 Mathematica 和 Maple 中实现高效计算，从而支持可积性分析及参数依赖型守恒律的检测，适用于 KdV、Boussinesq、sine-Gordon 和 Toda 晶格等系统。

ABSTRACT

We present direct methods and symbolic software for the computation of conservation laws of nonlinear partial differential equations (PDEs) and differential-difference equations (DDEs).The methods are applied to nonlinear PDEs in (1+1) dimensions with polynomial nonlinearities which include the Korteweg-de Vries (KdV), Boussinesq, and Drinfel'd-Sokolov-Wilson equations. An adaptation of the methods is applied to PDEs with transcendental nonlinearities. Examples include the sine-Gordon, sinh-Gordon, and Liouville equations. With respect to nonlinear DDEs, our methods are applied to Kac-van Moerbeke, Toda, and Ablowitz-Ladik lattices. To overcome the shortcomings of the undetermined coefficients method, we designed a new direct method which uses leading order analysis. That method is applied to discretizations of the KdV and modified KdV equations, and a combination thereof. Additional examples include lattices due to Bogoyavlenskii, Belov-Chaltikian, and Blaszak-Marciniak. The undetermined coefficient methods for PDEs and DDEs have been implemented in Mathematica. The code "TransPDEDensityFlux.m" computes densities and fluxes of systems of PDEs with or without transcendental nonlinearities. The code "DDEDensityFlux.m" does the same for polynomial nonlinear DDEs. Starting from the leading order terms, the new Maple library "discrete" computes densities and fluxes of nonlinear DDEs.

研究动机与目标

• 开发用于计算非线性 PDEs 和 DDEs 中多项式与超越型守恒律的直接、算法化方法。
• 通过创建高效、可扩展的软件工具，克服现有符号计算工具的局限性，以实现守恒律的高效计算。
• 通过非线性系统中守恒律层级的存在，实现可积性检测。
• 通过一种新颖的首项分析方法，将符号计算扩展至 DDEs，避免依赖伸缩不变性。
• 为研究人员提供开源软件工具，以分析非线性模型的物理与数学性质。

提出的方法

• 应用伸缩对称性分析，识别非线性 PDEs 和 DDEs 的尺度不变性。
• 通过待定系数法，将候选守恒密度表示为尺度不变单项式的线性组合。
• 利用连续变分导数（Euler 算子）推导系数的线性系统，根据非线性类型求解代数系统或常微分方程组。
• 采用连续同伦算子，通过分部积分法逆向总微分算子，以计算通量。
• 针对 DDEs，通过恒等算子的分解与基于首项计算的新型“分而治之”策略，对方法进行适配。
• 在 Mathematica（TransPDEDensityFlux.m、DDEDensityFlux.m）和 Maple（离散库）中实现算法，以支持自动化计算。

实验结果

研究问题

• RQ1 如何通过符号计算系统地计算具有多项式与超越型非线性的非线性 PDEs 的守恒律？
• RQ2 伸缩对称性在构造非线性 PDEs 和 DDEs 的候选守恒密度中起到何种作用？
• RQ3 同伦算子如何在 PDEs 和 DDEs 守恒律的背景下高效计算通量？
• RQ4 传统离散 Euler 算子与同伦算子在 DDEs 中存在哪些局限性？如何克服？
• RQ5 能否开发一种新方法，避免使用待定系数与伸缩不变性，以高效计算非多项式守恒律？

主要发现

• 该方法成功计算了经典 PDEs 的守恒律，包括 KdV、Boussinesq、Drinfel’d-Sokolov-Wilson、sine-Gordon、sinh-Gordon 和 Liouville 方程。
• 对于具有超越型非线性的 PDEs，该方法通过求解含待定函数系数的代数方程与常微分方程的混合系统，计算守恒律。
• 连续同伦算子通过分部积分法算法性地逆向总微分算子，实现通量的显式计算。
• 对于 DDEs，待定系数法在 Kac-van Moerbeke、Toda 和 Ablowitz-Ladik 晶格等多项式系统中表现有效。
• 针对 DDEs 的新型首项分析方法避免了不必要的项生成，无需依赖伸缩不变性或多项式性，即可高效计算守恒密度。
• Mathematica 中的 TransPDEDensityFlux.m 和 DDEDensityFlux.m 软件工具，以及 Maple 中的离散库，均已公开发布，并支持参数依赖型守恒律的检测。