[論文レビュー] Power Grid Simulation using Matrix Exponential Method with Rational Krylov Subspaces

Hao Zhuang, Shih-Hung Weng+1 more|arXiv (Cornell University)|Sep 20, 2013

Matrix Theory and Algorithms参考文献 9被引用数 3

ひとこと要約

本稿では、電力網シミュレーションに適応的時間刻み法と有理的クリロフ部分空間を用いた行列指数法を提案する。低周波応答を活用することで、精度を損なわずに大幅な高速化を実現する。標準的なクリロフ基底に代えて、小さな固有値を強調する有理的形を用いることで、剛性制約を軽減し、1回のLU分解と前進/後退代入のみを維持したまま、大きな適応的時間刻みを可能にする。固定刻みの台形法と比較して最大18倍の高速化を達成する。

ABSTRACT

One well adopted power grid simulation methodology is to factorize matrix once and perform only backward forward substitution with a deliberately chosen step size along the simulation. Since the required simulation time is usually long for the power grid design, the costly factorization is amortized. However, such fixed step size cannot exploit larger step size for the low frequency response in the power grid to speedup the simulation. In this work, we utilize the matrix exponential method with the rational Krylov subspace approximation to enable adaptive step size in the power grid simulation. The kernel operation in our method only demands one factorization and backward forward substitutions. Moreover, the rational Krylov subspace approximation can relax the stiffness constraint of the previous works. The cheap computation of adaptivity in our method could exploit the long low frequency response in a power grid and significantly accelerate the simulation. The experimental results show that our method achieves up to 18X speedup over the trapezoidal method with fixed step size.

研究の動機と目的

• 長時間の電力網シミュレーションにおける固定刻み時間積分の非効率性、特に低周波応答に対する問題を解決する。
• 従来の行列指数法が時間刻みを制限する剛性の制限を克服する。
• 1回の因子分解と代入操作を維持することで、計算効率を保ちながら電力網シミュレーションにおける適応的時間刻みを実現する。
• 複雑で剛性の高いシステムを有する産業規模の電力分配網（PDN）に対して、シミュレーションの精度と性能を向上させる。
• 実世界のPDN設計に一般的に見られる特異な静電容量行列を、有理的クリロフ部分空間近似が効果的に扱えることを示す。

提案手法

• 電力網の修正ノード解析から導かれる微分代数系を解くために、行列指数法（MEXP）を用いる。
• 行列Aに基づく標準的クリロフ部分空間を、(I − γA)⁻¹を基底とする有理的クリロフ部分空間に置き換えることで、低周波動的をよりよく捉える。
• スケーリング不変な有理的クリロフ近似を適用し、高条件数を有する剛性系でも行列指数を正確に計算可能にする。
• 1回の(I − γA)のLU因子分解の後、前進/後退代入により有理的クリロフ部分空間を構築し、繰り返しの因子分解を最小限に抑える。
• 局所誤差推定に基づき時間刻みを動的に調整することで適応的時間刻みを実装し、低周波領域では大きなステップを許容する。
• 行列指数の定式化に、数値積分を用いずに、区分線形入力源を解析的に統合する。

実験結果

リサーチクエスチョン

• RQ1 有理的クリロフ部分空間は、行列指数に基づく電力網シミュレーションにおける剛性を効果的に低減し、より大きな時間刻みを可能にするか？
• RQ2 行列指数法における適応的時間刻みは、固定刻み台形積分と比較して、シミュレーション性能をどの程度向上させるか？
• RQ3 提案手法は、産業的電力網設計に一般的に見られる特異な静電容量行列をどのように処理するか？
• RQ4 実世界のPDNベンチマークにおいて、適応的MEXPと有理的クリロフ部分空間を用いることで、シミュレーション実行時間にどの程度の高速化が達成できるか？
• RQ5 低周波、中周波、高周波応答成分のすべてにおいて、業界標準のHSPICEシミュレーションと比較して、本手法は高い精度を維持できるか？

主な発見

• 産業用PDNベンチマークにおいて、固定刻み台形法と比較して平均17.9倍、最大18.1倍の高速化を達成。ノード数は最大740万ノードまで対応。
• 2500ノードの剛性RCメッシュでは、標準的MEXPの平均クリロフ部分空間次元115が、本手法では3.11にまで低下し、行列指数評価で224倍の高速化を達成。
• 従来のMEXP手法が正則化を必要とし、依然として剛性を示すのに対し、本手法は正則化なしに特異な静電容量行列を効果的に処理できる。
• シミュレーション結果では、HSPICEとの相対誤差最大7.33×10⁻⁴を示し、台形法の7.47×10⁻⁴と同等の高精度を確認。
• 長期間の低周波応答を活用した適応的時間刻みが可能で、固定刻み法とは異なり、安定領域では10 psよりはるかに大きな時間刻みが可能となる。
• 実行時間の内訳から、LU因子分解とDC解析の時間が総シミュレーション時間に比べて最小限であることが判明し、カーネル演算が効率的にアーモナイズされていることが確認された。