[论文解读] Index reduction of differential algebraic equations by differential algebraic elimination

Xiaolin Qin, Yang Lu+3 more|arXiv (Cornell University)|Apr 20, 2015

Modeling and Simulation Systems参考文献 20被引用 4

一句话总结

该论文提出了一种新颖的高指标微分代数方程（DAE）索引降阶方法，采用微分代数消去法，利用Dixon结式计算微分代数结式。通过引入变量束技术及启发式方法去除多余因子，该方法在保持结构特性的同时高效降低DAE索引，优于现有方法在基准测试实例中的表现。

ABSTRACT

High index differential algebraic equations (DAEs) are ordinary differential equations (ODEs) with constraints and arise frequently from many mathematical models of physical phenomenons and engineering fields. In this paper, we generalize the idea of differential elimination with Dixon resultant to polynomially nonlinear DAEs. We propose a new algorithm for index reduction of DAEs and establish the notion of differential algebraic elimination, which can provide the differential algebraic resultant of the enlarged system of original equations. To make use of structure of DAEs, variable pencil technique is given to determine the termination of differentiation. Moreover, we also provide a heuristics method for removing the extraneous factors from differential algebraic resultant. The experimentation shows that the proposed algorithm outperforms existing ones for many examples taken from the literature.

研究动机与目标

• 为解决物理与工程模型中出现的高指标DAE问题，这些方程因索引相关的刚性与病态性而难以数值求解。
• 将微分消去法与Dixon结式推广至多项式非线性DAE，实现系统的索引降阶。
• 设计一种结构化算法，在求导过程中保持DAE的原始结构特性。
• 引入变量束技术以确定最优的求导终止点，避免不必要的复杂性。
• 提出一种启发式方法，用于从微分代数结式中去除多余因子，提升解的准确性。

提出的方法

• 该方法采用微分代数消去法，系统地从原始DAE生成新方程组，利用Dixon结式计算扩展系统微分代数结式。
• 引入变量束技术以引导求导过程，确保在最优点终止，降低计算开销。
• 通过将导数视为新变量，计算微分代数结式，实现对扩展系统的代数化处理。
• 应用启发式策略识别并去除结式中的多余因子，提升降阶后系统的可靠性。
• 该方法旨在保持原始DAE的结构特性，如稀疏性与多项式非线性形式。
• 在文献中的基准DAE实例上实现并评估该方法，以检验其性能与鲁棒性。

实验结果

研究问题

• RQ1 如何将微分消去法推广至多项式非线性DAE，以实现有效的索引降阶？
• RQ2 可采用何种标准来确定最优求导次数，以避免过度求导？
• RQ3 如何有效识别并去除结式计算过程中引入的多余因子？
• RQ4 该方法在标准DAE测试用例上是否能实现优于现有索引降阶技术的性能？
• RQ5 变量束技术在多大程度上提升了索引降阶过程的效率与稳定性？

主要发现

• 所提出的算法通过系统性地应用微分代数消去法，并利用Dixon结式计算微分代数结式，成功实现了高指标DAE的索引降阶。
• 变量束技术有效确定了求导过程的终止点，防止了不必要的复杂性，提升了计算效率。
• 用于去除多余因子的启发式方法显著增强了降阶后系统的准确性与可靠性。
• 实验结果表明，该方法在一系列文献中的基准DAE实例上优于现有方法。
• 该方法保持了原始DAE的结构特性，如稀疏性与多项式形式，这有助于其高效性与鲁棒性。