[論文レビュー] Numerical approach to $L_1$-problems with the second order elliptic operators

N. Limić, Mladen Rogina|ArXiv.org|Dec 21, 2007

Differential Equations and Boundary Problems参考文献 14被引用数 3

ひとこと要約

本稿は、境界条件がディリクレ型である第二順階楕円型作用素の $L^1$-問題を解くための数値的手法を提示する。有限差分スキームを用いて、コンpartimental（M行列に類似した）構造を持つグリッド解を構築し、それをハット関数に埋め込む。$p \in [1, d/(d-1))$ に対して、強い $L^p$ 収束および弱い $\.dot{W}^1_p$ 収束を証明する。主な貢献は、任意の次元において物理的構造を保存する、証明可能に収束する単調なスキームを提供し、厳密な誤差制御を実現することである。

ABSTRACT

For a second order differential operator $A(\msx) =- abla a(\msx) abla + b'(\msx) abla+ abla \big(\msb''(\msx) \cdot\big)$ on a bounded domain $D$ with the Dirichlet boundary conditions on $\partial D$ there exists the inverse $T(λ, A)= (λI+A)^{-1}$ in $L_1(D)$. If $μ$ is a Radon (probability) measure on Borel algebra of subsets of $D$, then $T(λ, A)μ\in L_p(D), p \in [1, d/(d-1))$. We construct the numerical approximations to $u =T(λ, A)μ$ in two steps. In the first one we construct grid-solutions ${\bf u}_n$ and in the second step we embed grid-solutions into the linear space of hat functions $u(n) \in \dot{W}_p^1(D)$. The strong convergence to the original solutions $u$ is established in $L_p(D)$ and the weak convergence in $\dot{W}_p^1(D)$.

研究の動機と目的

• 第二順階楕円型作用素の発散型表現を含む $L^1$-問題を解くための安定的で単調な数値スキームの開発。
• 離散問題の係数行列がコンpartimental構造（正の対角成分、非正の非対角成分、非負の列和）を持つことを保証し、M行列に類似した性質を再現すること。
• $p \in [1, d/(d-1))$ に対して、$L^p(D)$ での強い収束および $\dot{W}^1_p(D)$ での弱い収束を確立すること。
• リプシッツ境界を持つ $\mathbb{R}^d$ 内の一般領域および特異データ（レイン・メジャー）にまで既存の数値法を拡張すること。
• 次元 $d=2$ に対して構成的アルゴリズムを提供し、$d \geq 3$ に対して一般構成法を提示し、拡散テンソルに関する十分条件を提示すること。

提案手法

• 構造化グリッド上で第二順階楕円型作用素 $A(x)$ を有限差分法で離散化し、係数行列 $A_n$ が一貫性およびコンpartimental構造を持つように保証する。
• 離散線形系 $A_n u_n = \mu_n$ の解としてグリッド解 $u_n$ を構築する。ここで $\mu_n$ は離散レイン・メジャーである。
• グリッド解 $u_n$ を補間により連続な区分線形（ハット）関数 $u^{(n)} \in \dot{W}^1_p(D)$ に埋め込む。
• コンパクトネスおよび一貫性の議論を用いて、$u^{(n)}$ が真の解 $u = T(\lambda, A)\mu$ に対して $L^p(D)$ で強い収束および $\dot{W}^1_p(D)$ で弱い収束することを証明する。
• 連続および離散の変分的定式化を結びつけるために、双線形形式 $a(v,u)$ とグリーンの公式を用いる。
• $d \geq 3$ の場合、係数行列 $A_n$ のコンpartimental構造を保証するため、拡散テンソル $\{a_{ij}(x)\}$ に仮定 4.1 を課す。

実験結果

リサーチクエスチョン

• RQ1 第二順階楕円型作用素の発散型表現を含む一般の拡散テンソルを有する場合に、任意の次元において単調な有限差分スキームを構築できるか。
• RQ2 グリッド解がハット関数に埋め込まれた場合、$p \in [1, d/(d-1))$ に対して、$L^p(D)$ で強い収束および $\dot{W}^1_p(D)$ で弱い収束が成立するか。
• RQ3 拡散テンソルにどのような条件を課すと、係数行列 $A_n$ がコンpartimental構造を持つようになるか。
• RQ4 特異データ（例：デルタ関数）を有する問題に対して、提案スキームの数値的性能はいかがなものか。
• RQ5 $d=2$ において、直接的なアルゴリズム的構造を持つように実装可能か。

主な発見

• $d=2$ の場合、スキームはアルゴリズム的に明示的であり、直接適用可能であり、$W^1_2$ および $W^1_p$ 空間で収束が証明されている。
• $d \geq 3$ の場合、テンソル値関数 $\hat{a}$ が正定値であれば、補題 8.1 により $A_n$ のコンpartimental構造が保証される。
• 数値例 7.1 では、デルタ測度を有する問題に対して $\varepsilon_{1,\text{rel}} = 0.003371$ および $\varepsilon_{\infty,\text{rel}} = 1.6254$ が得られ、最大誤差は特異点付近に集中している。
• 例 7.2 では 210 回の反復後に $\varepsilon_{1,\text{rel}} = 8 \times 10^{-6}$ および $\varepsilon_{\infty,\text{rel}} = 4 \times 10^{-3}$ を達成し、高い精度が確認された。
• $d=2$ では回転を避けるために拡張ステンシルを用いることで、単調性が保たれ、解の表面に Saw-like のアーチファクトが生じない。
• 本手法は特異データ（レイン・メジャー）に対してもロバストであり、逆作用素 $T(\lambda, A)$ は $\mu \in R(D)$ を $L^p(D)$ および $\dot{W}^1_p(D)$ に写像する。ただし $p < d/(d-1)$ の場合に限る。