[論文レビュー] Dynamical low-rank tensor approximations to high-dimensional parabolic problems: existence and convergence of spatial discretizations

Markus Bachmayr, Henrik Eisenmann+1 more|arXiv (Cornell University)|Jan 1, 2023

Tensor decomposition and applications被引用数 1

ひとこと要約

本稿は、高次元放物型PDEに対する動的低ランクテンソル近似の空間離散化の存在および収束を確立する。ヒルベルト空間において、空間的に離散化されたDLRA定式化の解が連続的DLRA解に収束することを証明し、従来の結果を、多変数テンソル系列および階層的形式にまで拡張し、多様体に基づく時間発展の厳密な解析を実施する。

ABSTRACT

We consider dynamical low-rank approximations to parabolic problems on higher-order tensor manifolds in Hilbert spaces. In addition to existence of solutions and their stability with respect to perturbations to the problem data, we show convergence of spatial discretizations. Our framework accommodates various standard low-rank tensor formats for multivariate functions, including tensor train and hierarchical tensors.

研究の動機と目的

• 高次元放物型問題におけるヒルベルト空間上での動的低ランクテンソル近似の存在および安定性を確立すること。
• 放物型PDEのDLRAにおける空間離散化の収束という、長年の未解決問題に取り組むこと。
• 変分的・多様体に基づく設定において、DLRAの理論的枠組みをテンソル系列および階層的テンソル形式へと拡張すること。
• 空間半離散化下での低ランクテンソル解の時間発展に対する厳密な解析を提供すること。
• 高次元PDEの文脈において、離散的数値スキームと連続的DLRA定式化の間のギャップを埋めること。

提案手法

• DLRA問題をヒルベルト空間設定下での低ランクテンソル多様体上での変分的時間発展として定式化する。
• テスト関数を多様体の接空間に制限したディラック＝フレンケルの原理（時間依存ガレルキン法）を適用する。
• 放物型PDEの弱定式化にGelfandの三重対を用い、データをL²(0,T;H⁻¹(Ω))に、初期データをL²(Ω)にとる。
• ヒルベルト空間上での変分的時間ステッピングスキームを用いて、最大時間区間における解の存在および一意性を証明する。
• 多線形パrametrizationおよび滑らかな曲線表現を用いて、テンソル系列多様体の接空間構造を解析する。
• 接空間への射影作用素のリプシッツ連続性および有界性を確立し、収束解析を可能にする。

実験結果

リサーチクエスチョン

• RQ1 高次元放物型PDEに対して、ヒルベルト空間設定下での動的低ランクテンソル近似の解が存在するか？
• RQ2 メッシュを細かくするに従って、空間離散化されたDLRAスキームは連続的DLRA解に収束するか？
• RQ3 テンソル系列および階層的テンソルなどの低ランクテンソル形式に対して、空間離散化の収束を厳密に証明できるか？
• RQ4 低ランクテンソル多様体の接空間構造は、DLRAスキームの安定性および収束にどのように影響するか？
• RQ5 滑らかでないデータを有する放物型問題の文脈において、DLRA解の存在および一意性を保証する条件は何か？

主な発見

• 本稿は、ヒルベルト空間フレームワーク下で、最大時間区間における放物型問題に対する動的低ランクテンソル近似の解の存在および一意性を証明する。
• 空間離散化の収束が連続的DLRA解に至ることを確立し、文献における重要な未解決問題を解決する。
• 最小限の正則性仮定の下で収束を示す：初期データはL²(Ω)、外力項はL²(0,T;H⁻¹(Ω))に属する。
• 一般的な多様体フレームワークを用いて、テンソル系列および階層的テンソルを含むさまざまな低ランクテンソル形式に適用可能である。
• 接空間への射影作用素が、コアテンソルの最小特異値σに対してO(1/σ)の定数でリプシッツ連続であることが示された。
• 固定ランクテンソルの多様体がヒルベルト空間のC¹部分多様体であることが証明され、解析に微分幾何学的道具を適用可能になる。