[論文レビュー] 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

ひとこと要約

本論文は、微分代数的消去法を用い、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の課題に対処すること。これらはインデックスに起因する剛性と悪条件性のため、数値的に解くのが困難である。
• 多項式的非線形DAEを扱えるように、Dixon結果式を用いた微分消去法を一般化し、体系的なインデックス低減を可能にすること。
• 微分の過程で元のDAEの構造的性質を保持する構造的なアルゴリズムを開発すること。
• 最適な微分終了点を特定するための変数ペンシル技術を導入し、不必要な複雑さを回避すること。
• 微分代数的結果式に含まれる余分な因子を効果的に同定・除去するためのヒューリスティクス手法を提供すること。

提案手法

• 本手法は、微分代数的消去法を用い、元のDAEから新たな方程式系を体系的に生成する。Dixon結果式を用いて、拡大された方程式系の微分代数的結果式を計算する。
• 微分プロセスを導くために変数ペンシル技術が導入され、最適な点で微分を終了させ、計算のオーバーヘッドを低減する。
• 導関数を新たな変数として扱うことで、拡張された方程式系の代数的取り扱いが可能になり、微分代数的結果式が計算される。
• 余分な因子を同定・除去するためのヒューリスティクス戦略が適用され、低減された方程式系の信頼性が向上する。
• 元のDAEの構造的性質（スパarsityや多項式的非線形性など）を保持するように設計されている。
• 本手法は、文献に掲載されたベンチマークDAE例を用いて実装・評価され、性能と頑健性が評価されている。

実験結果

リサーチクエスチョン

• RQ1 微分消去法を多項式的非線形DAEに一般化することで、効果的なインデックス低減が可能になるか。
• RQ2 過剰な微分を回避するための最適な微分回数の判定基準は何か。
• RQ3 結果式計算中に生じる余分な因子を効果的に同定・除去する方法は何か。
• RQ4 提案手法は、標準的なDAEテストケースにおいて、既存のインデックス低減技術を上回る性能を示せるか。
• RQ5 変数ペンシル技術は、インデックス低減プロセスの効率性と安定性をどの程度向上させるか。

主な発見

• 提案されたアルゴリズムは、微分代数的消去法を体系的に適用し、Dixon結果式を用いて微分代数的結果式を計算することで、高インデックスDAEのインデックス低減に成功している。
• 変数ペンシル技術は、微分の終了点を効果的に特定し、不必要な複雑さを防ぎ、計算効率を向上させている。
• 余分な因子を除去するためのヒューリスティクス手法は、低減された方程式系の正確性と信頼性を顕著に向上させている。
• 実験結果から、本手法は文献に掲載されたさまざまなベンチマークDAE例において、既存の手法を上回ることが示された。
• 本手法は、元のDAEの構造的性質（スパarsityや多項式形式など）を保持しており、これがその効率性と頑健性に寄与している。