[論文レビュー] Quasi-Monte Carlo finite element analysis for wave propagation in heterogeneous random media

M. Ganesh, Frances Y. Kuo+1 more|arXiv (Cornell University)|Apr 26, 2020

Mathematical Approximation and Integration参考文献 27被引用数 4

ひとこと要約

本稿では、有界領域における空間的に不均一な屈折率を有するヘルムホルツ方程式を解くための準モンテカルロ（QMC）有限要素法を提案する。最近開発された符号不変な変分形式を用いることで、標準的なヘルムホルツ形式の強制性の欠如を克服し、拡散問題と同等の誤差境界を有する収束解析を可能にする。

ABSTRACT

We propose and analyze a quasi-Monte Carlo (QMC) algorithm for efficient simulation of wave propagation modeled by the Helmholtz equation in a bounded region in which the refractive index is random and spatially heterogenous. Our focus is on the case in which the region can contain multiple wavelengths. We bypass the usual sign-indefiniteness of the Helmholtz problem by switching to an alternative sign-definite formulation recently developed by Ganesh and Morgenstern (Numerical Algorithms, 83, 1441-1487, 2020). The price to pay is that the regularity analysis required for QMC methods becomes much more technical. Nevertheless we obtain a complete analysis with error comprising stochastic dimension truncation error, finite element error and cubature error, with results comparable to those obtained for the diffusion problem.

研究の動機と目的

• ランダムで空間的に変化する屈折率を有する領域における波動伝播をシミュレートするための効率的な数値解法の開発。
• 複数の波長が存在する際、標準ガレルキン形式における強制性の欠如を解消すること。
• 従来、拡散問題に対して成功を収めた準モンテカルロ（QMC）法を、符号不変な変分形式を用いることで波動伝播に拡張すること。
• 確率的次元の切断、有限要素離散化、および立方則誤差を含む包括的な誤差解析を提供すること。
• 符号不変な形式に起因する複雑な確率的依存性を考慮に入れつつも、拡散問題と同等の収束率を達成すること。

提案手法

• ガネシュとモルゲンシュテルンが最近開発した、符号不変なヘルムホルツ方程式の変分形式を採用し、標準形式の符号不定性を回避する。
• 屈折率は、無限次元の一様確率ベクトル $\mathbf{y} \in [-1/2, 1/2]^\mathbb{N}$ でパrameter化された確率場としてモデル化され、平均成分と摂動成分に分解される。
• 符号不変な弱形式に確率的ガレルキン有限要素法を適用することで、安定かつ一貫性のある空間離散化が可能になる。
• 解の線形汎関数の期待値を計算するために、低不規則性列を活用した準モンテカルロ（QMC）立方則ルールが用いられる。
• 誤差は3つの成分に分解される：確率的次元の切断誤差、有限要素近似誤差、QMC立方則誤差。各成分は厳密に評価されている。
• 符号不変な形式がもたらす複雑な確率的依存性を扱うために、重み付き多重インデックスノルムを用いた技術的帰納法を用いた新しい正則性解析が開発された。

実験結果

リサーチクエスチョン

• RQ1 符号不変な変分形式を用いることで、ランダムかつ不均質な係数を有するヘルムホルツ方程式に準モンテカルロ法を効果的に適用できるか。
• RQ2 標準的ヘルムホルツ形式における強制性の欠如を克服し、QMC解析を可能にする方法は何か。
• RQ3 屈折率がランダムかつ空間的に不均一な場合、QMC有限要素法の収束挙動はどのようになるか。
• RQ4 この設定において、確率的、空間的、立方則的離散化を統合した誤差境界はどのように確立できるか。
• RQ5 無限次元確率入力に関する解の正則性を十分に分析でき、QMC収束率を支持できるか。

主な発見

• 提案手法は、符号不変な形式に起因する複雑さを考慮に入れても、拡散問題と同等の収束率を達成する。
• 最近開発された符号不変な変分形式を用いることで、標準形式の符号不定性を効果的に回避し、安定なガレルキン離散化を保証する。
• 確率的切断誤差、有限要素近似誤差、QMC立方則誤差の各成分に対して厳密な境界が確立され、総誤差は解の確率的微分の可積分性に依存した速度で減少する。
• 高階の確率的モーメントの成長を制御するための新しい技術的補題が確立され、QMC収束に必要な正則性解析が可能になった。
• 波長が領域サイズに比べて小さい場合でも、安定性と収束性が維持され、標準的手法が強制性の欠如により失敗する領域でも有効である。
• QMC法が標準モンテカルロ法に比べて収束速度が優れていることが分析で確認され、高次元確率的設定では顕著な計算コスト削減が期待できる。