[論文レビュー] Solution of boundary value and eigenvalue problems for second order elliptic operators in the plane using pseudoanalytic formal powers

Raúl Castillo-Pérez, Vladislav V. Kravchenko+1 more|arXiv (Cornell University)|Feb 5, 2010

Algebraic and Geometric Analysis参考文献 13被引用数 15

ひとこと要約

この論文は、平面における2階楕円型作用素の境界値問題および固有値問題を解くための擬解析的形式的べき乗法（MPFP）を導入する。擬解析関数理論を活用し、生成列の再帰的積分によってゼロ解の完全な系を構成することで、コロケーションを用いたメッシュレスな解法が可能となり、高い精度と高速収束性を示し、有限要素法と比較しても競争力を持つ。

ABSTRACT

We propose a method for solving boundary value and eigenvalue problems for the elliptic operator D=divpgrad+q in the plane using pseudoanalytic function theory and in particular pseudoanalytic formal powers. Under certain conditions on the coefficients p and q with the aid of pseudoanalytic function theory a complete system of null solutions of the operator can be constructed following a simple algorithm consisting in recursive integration. This system of solutions is used for solving boundary value and spectral problems for the operator D in bounded simply connected domains. We study theoretical and numerical aspects of the method.

研究の動機と目的

• 平面における2階楕円型作用素の境界値問題および固有値問題を解くためのメッシュレス数値法の開発。
• 作用素 $ D = \operatorname{div} p \operatorname{grad} + q $ のゼロ解の完全な系を、擬解析関数理論を用いて構成すること。
• 特に変数分離が不可能な領域において、従来の離散化に基づく手法（例：有限要素法）の代替を提供すること。
• 数値実験を通じて、本手法の効率性、収束速度および精度を示すこと。
• ゼロ解の完全な解系を用いて固有値を行列式の根の探索によって近似することで、固有値問題への応用を拡張すること。

提案手法

• 本手法は、Vekua方程式から導かれる擬解析的形式的べき乗を用いて、楕円型作用素 $ D = \operatorname{div} p \operatorname{grad} + q $ の完全な解系を構築する。
• 作用素 $ D u = 0 $ の特別な解 $ u_0 $ が必要であり、$ f = p^{1/2} u_0 $ が滑らかで消えないこと、かつ直交座標系で $ f = S(s)T(t) $ と変数分離可能であること。
• 最近の文献から得られるアルゴリズムに従い、生成列を再帰的に構築することで、解系の明示的構成が可能となる。
• 境界値問題の解は、これらの形式的べき乗の線形結合として近似され、境界上のコロケーションによって係数が決定される。
• 固有値問題では、行列 $ U = (u_k(z_j)) $（$ u_k $ は形式的べき乗、$ z_j \in \partial\Omega $）の行列式がゼロとなる $ \lambda $ を探索する。
• 本手法はラプラス方程式に対する調和多項式の手法を一般化し、$ p \equiv 1, q \equiv 0, u_0 \equiv 1 $ の場合に $ \operatorname{Re}(z-z_0)^n, \operatorname{Im}(z-z_0)^n $ に還元される。

実験結果

リサーチクエスチョン

• RQ1 有界な平面領域において、係数が変数分離可能であるという条件下で、作用素 $ D = \operatorname{div} p \operatorname{grad} + q $ のゼロ解の完全な系を明示的に構成できるか？
• RQ2 MPFPの収束速度および精度は、有限要素法などの従来の離散化手法と比較してどのように異なるか？
• RQ3 MPFPは、特にヘルムホルツ方程式に対して、固有値問題を効果的に解くために適応可能か？
• RQ4 関数 $ f = p^{1/2} u_0 $ の分離可能性が、形式的べき乗の構成と解の完全性を保証するために果たす役割は何か？
• RQ5 実際の応用において、境界コロケーション行列 $ U $ の行列式をどれほど信頼性を持って固有値を計算に用いることができるか？

主な発見

• MPFPは境界値問題において高い精度を達成しており、最大絶対誤差が $ N = 4 $ 時点で 0.079 から $ N = 32 $ 時点で $ 2.8 \times 10^{-10} $ に低下し、急速な収束が示された。
• 単位円板におけるヘルムホルツ方程式では、$ N = 21 $ で最初の5つの固有値を4桁の精度で計算可能であり、$ N = 23 $ で6つの固有値を同じ精度で得られた。
• 特に形式的べき乗の明示的構成が可能な問題において、有限要素法と比較して競争力のある性能を示した。
• 係数 $ p $ および $ q $ に関する所定の条件下で、擬解析的形式的べき乗の系は、解析関数論における $ (z - z_0)^n $ の完全性に類似した意味で完全である。
• 本手法は古典的な調和多項式の手法を一般化しており、$ p \equiv 1 $, $ q \equiv 0 $, $ u_0 \equiv 1 $ の場合に標準的な調和多項式に回復する。
• 複素数値の $ u_0 $ の場合に理論的困難が生じるものの、数値実験により、本手法は固有値問題に対しても有効であることが確認された。行列式 $ \det U $ は急速に減少するが、依然として高精度な固有値近似が得られた。