[論文レビュー] An efficient solver based on low-rank approximation and Neumann matrix series for unsteady diffusion-type partial differential equations with random coefficients

Yujun Zhu, Min Li+2 more|arXiv (Cornell University)|Jan 16, 2026

Probabilistic and Robust Engineering Design被引用数 0

ひとこと要約

論文は、乱数係数を持つ拡散型の時変SPDEを効率的に解くために、確率剛度行列の一般化低ランク近似とNeumann級数を組み合わせたLRNS求解器を提案し、計算量とメモリを削減しつつ精度を維持する。

ABSTRACT

In this paper, we develop an efficient numerical solver for unsteady diffusion-type partial differential equations with random coefficients. A major computational challenge in such problems lies in repeatedly handling large-scale linear systems arising from spatial and temporal discretizations under uncertainty. To address this issue, we propose a novel generalized low-rank matrix approximation to represent the stochastic stiffness matrices, and approximate their inverses using the Neumann matrix series expansion. This approach transforms high-dimensional matrix inversion into a sequence of low-dimensional matrix multiplications. Therefore, the solver significantly reduces the computational cost and storage requirements while maintaining high numerical accuracy. The error analysis of the proposed solver is also provided. Finally, we apply the method to two classic uncertainty quantification problems: unsteady stochastic diffusion equations and the associated distributed optimal control problems. Numerical results demonstrate the feasibility and effectiveness of the proposed solver.

研究の動機と目的

• 乱数係数を持つ時系列拡散方程式の不確定性定量化を強健かつ効率的に行う動機づけ。
• MC-FEM離散化から生じる大規模な確率的線形系の計算負荷を抑える軽量ソルバーの開発。
• 提案LRNSアプローチの誤差解析を提供し、確率拡散問題や確率最適制御問題への適用可能性を示す。

提案手法

• グリニア低ランク近似(GLRAM)を用いRSVDからUとV_mを得て、圧縮比tauを所定に設定する新規な stochastic perturbation representation.
• 摂動行列の逆行列を切り捨てNeumann級数で近似する: (I + A^{-1}Ã_m^*)^{-1} ≈ sum_{r=0}^R (-A^{-1}Ã_m^*)^r.
• 大規模なM×L線形系を、低次元行列の掛け算と決定論解への低ランク補正の連続として転換。
• RMSREと摂動行列N = sum_m Ã_m Ã_m^Tの固有スペクトルを用いた誤差解析フレームワークを提供。
• RSVDベースのLRAとNeumann-seriesベースの逆算を組み合わせて、スケーラブルなLRNSソルバー（Algorithm 2）を得る。
• 時変確率拡散方程式と分散型確率最適制御問題への適用性を示す。

実験結果

リサーチクエスチョン

• RQ1 モンテカルロFEM設定で、低ランク表現は時系列SPDEの解をどのように加速するか。
• RQ2 RSVDベースのGENERALIZED低ランク近似とNeumann-series逆行列の組み合わせによる摂動系の精度と計算効率はどの程度か。
• RQ3 提案LRNSアプローチは、乱れた透磁率を持つ時変拡散方程式や関連する確率的制御問題を、精度を損なうことなく効率的に解けるか。

主な発見

• LRNSソルバーは高次元の行列反転を低次元の積と低ランク更新の列へ変換し、コストとストレージを削減。
• RSVDベースの一般化低ランク近似は、所定の圧縮比tauを用いて乱数付き剛度の摂動を効率的かつ正確に表現。
• Neumann級数の逆算は、切り捨てインデックスRと摂動サイズに応じて精度が改善する実用的な有限項近似を提供。
• 誤差解析はRMSREを摂動行列Nのスペクトル特性および選択したエネルギー捕捉e(tau)と結びつける。
• 本手法は、時変確率拡散方程式と分散型確率最適制御問題への適用性と効率性を示す。