[論文レビュー] Space-Time Least-Squares Isogeometric Method for Parabolic Problems

Monica Montardini, Matteo Negri+2 more|arXiv (Cornell University)|Sep 26, 2018

Advanced Numerical Analysis Techniques被引用数 1

ひとこと要約

本稿では、高次スムーズなスプラインを用いた放物型問題に対する空間時間最小二乗等階的配置法を提案する。テンソル積構造を活用して、高速対角化を用いてシルベスタ型方程式を解く、堅牢なプリコンディショナを設計する。この手法は、多項式次数に依存せず、自由度にほぼ比例する計算コストを保証し、高次等階的配置離散化において非常に効率的かつ並列化可能である。

ABSTRACT

In this paper, we propose a space-time least-squares isogeometric method to solve parabolic evolution problems, well suited for high-degree smooth splines in the space-time domain. We focus on the linear solver and its computational efficiency: thanks to the proposed formulation and to the tensor-product construction of space-time splines, we can design a preconditioner whose application requires the solution of a Sylvester-like equation, which is performed efficiently by the fast diagonalization method. The preconditioner is robust w.r.t. spline degree and mesh size. The computational time required for its application, for a serial execution, is almost proportional to the number of degrees-of-freedom and independent of the polynomial degree. The proposed approach is also well-suited for parallelization.

研究の動機と目的

• 放物型問題に対する高次等階的配置法の計算非効率性を解消する。これは、悪条件な線形系に起因する。
• 空間と時間の両方で滑らかで高次スプラインを用いた、高精度な解を得るための空間時間定式化を構築する。
• 多項式次数やメッシュサイズの変化に対しても性能を維持する、堅牢で効率的なプリコンディショナを設計する。
• プリコンディショナの適用コストが自由度にほぼ線形に比例し、スプライン次数に依存しないように保証する。
• 大規模問題に対する解法プロセスを、テンソル積構造を活用して効率的に並列化可能にする。

提案手法

• 高次Bスプラインを用いた等階的配置法により、空間時間領域における放物型問題を定式化する。
• 空間時間弱形式に最小二乗変分法を適用し、システム行列の対称性と正定値性を保証する。
• 空間時間スプラインのテンソル積構造を活用して、システムをシルベスタ型行列方程式に分解する。
• テンソル構造に基づくプリコンディショナを設計し、線形システムの解法をシルベスタ方程式の解法に還元する。
• 空間的および時間的作用素の固有値分解を活用して、高速対角化法を用いてシルベスタ方程式を効率的に解く。
• 多項式次数やメッシュの細分化にかかわらず、プリコンディショナが有効かつ計算的に効率的であることを保証する。

実験結果

リサーチクエスチョン

• RQ1 高次スムーズスプラインを用いた空間時間最小二乗等階的配置法は、放物型問題に対して堅牢な性能を発揮できるか？
• RQ2 得られる線形システムの解法コストを最小限に抑えつつ、高次精度を維持できるか？
• RQ3 提案されたプリコンディショナは、多項式次数やメッシュサイズの変化に対しても有効かつ効率的か？
• RQ4 テンソル積構造を活用して、解法プロセスを効率的に並列化できるか？
• RQ5 自由度および多項式次数の観点から、プリコンディショナの適用のスケーラビリティと計算複雑度はいかほどか？

主な発見

• 提案手法は、空間および時間の両方で高次スムーズスプラインを用いることで、優れた近似特性を実現する高次精度を達成する。
• プリコンディショナの適用コストは、自由度に対してほぼ線形に増加し、多項式次数に依存しない。
• 高速対角化法により、テンソル積構造に起因するシルベスタ型方程式を効率的に解くことが可能になる。
• プリコンディショナは、メッシュの細分化およびスプライン次数の増加に対し、堅牢で、広範な離散化において効率性を維持する。
• 空間時間システムの分離構造のおかげで、並列実行に適した手法である。
• 計算実験により、ソルバの性能がほぼ最適であり、多項式次数の増加に伴い著しい劣化がないことが確認された。