[論文レビュー] Mathematical and numerical analysis of time-dependent Ginzburg--Landau equations in nonconvex polygons based on Hodge decomposition

Buyang Li, Zhimin Zhang|arXiv (Cornell University)|Oct 14, 2014

Advanced Numerical Methods in Computational Mathematics参考文献 17被引用数 5

ひとこと要約

本稿は、磁気ポテンシャルがH¹正則性を欠く非凸多角形領域における時間依存Ginzburg–Landau方程式の適切な定式化を確立する。これにより、有限要素法による誤った解が生じる。Hodge分解を用いて、H¹適合未知数のみを含む等価な系に再定式化することで、最適収束率を達成する分離型・線形化有限要素法を構築した。数値的検証は特異性を持つL字型領域で行われ、標準的な有限要素法が失敗する状況でも有効であることが示された。

ABSTRACT

We prove well-posedness of time-dependent Ginzburg--Landau system in a nonconvex polygonal domain, and decompose the solution as a regular part plus a singular part. We see that the magnetic potential is not in $H^1$ in general, and the finite element method (FEM) may give incorrect solutions. To remedy this situation, we reformulate the equations into an equivalent system of elliptic and parabolic equations based on the Hodge decomposition, which avoids direct calculation of the magnetic potential. The essential unknowns of the reformulated system admit $H^1$ solutions and can be solved correctly by the FEMs. We then propose a decoupled and linearized FEM to solve the reformulated equations and present error estimates based on proved regularity of the solution. Numerical examples are provided to support our theoretical analysis and show the efficiency of the method.

研究の動機と目的

• 磁気ポテンシャルがH¹に属さない非凸多角形領域における時間依存Ginzburg–Landau系の適切な定式化を確立すること。これは、再入角を持つ領域で標準的な有限要素法が失敗する要因である。
• 磁気ポテンシャルがH¹に属さない場合、従来の有限要素法が破綻する問題に対処すること。これは再入角を持つ領域で一般的に生じる。
• Hodge分解を用いて、特異な磁気ポテンシャルの直接計算を回避し、主要未知数のH¹正則性を保証する再定式化された系を構築すること。
• 再定式化された系の正則性に基づいて、厳密な誤差推定が得られる分離型・線形化有限要素スケームを設計すること。
• L字型領域における数値的検証を通じて、標準的なFEMが失敗する状況でも本手法の収束性を確認すること。

提案手法

• Hodge分解を用いて、元の時間依存Ginzburg–Landau系を、磁気ポテンシャルを正則部と特異部に分離する形に再定式化する。
• H¹適合未知数のみを含む楕円型および放物型方程式の等価系を構築し、非H¹磁気ポテンシャルの直接近似を回避する。
• 再定式化された系に対して、分離型かつ線形化された有限要素法を適用し、安定かつ収束する数値解を得る。
• 時間方向には後退Euler積分法、空間方向には一次のピecewise線形有限要素法を用いる。
• 再定式化された系の未知数がH¹に属することを示し、その正則性に基づいて最適順序の誤差推定を証明する。
• 正確な解が与えられるL字型領域で数値実験を実施し、標準FEMと本手法の収束特性を比較する。

実験結果

リサーチクエスチョン

• RQ1 再入角を持つ非凸多角形領域において、磁気ポテンシャルがH¹に属さない場合でも、時間依存Ginzburg–Landau系は一意解をもつか？
• RQ2 磁気ポテンシャルがH¹に属さない場合、特にコーナー特異性を持つ領域では、標準的な有限要素法が収束する解を生成できるか？
• RQ3 特異な磁気ポテンシャルの直接計算を回避し、未知数のH¹正則性を保証するように、Ginzburg–Landau方程式を等価な系に再定式化できるか？
• RQ4 非凸領域における再定式化された系に有限要素法を適用した場合、厳密に示せる収束率は何か？
• RQ5 コーナー特異性の存在下で、本手法の収束性および精度は、標準FEMと比較してどのように異なるか？

主な発見

• 磁気ポテンシャルがH¹に属さない非凸多角形領域においても、解を正則部と特異部に分解することで、時間依存Ginzburg–Landau系は適切に定式化される。
• 元の系に直接適用した標準的な有限要素法は、メッシュを細かくしても収束せず、秩序パラメータでは約2.8×10⁻³、磁気ポテンシャルでは約7.8×10⁻²で誤差が停滞する。
• 本手法で用いられるHodge分解に基づく再定式化により、特異な磁気ポテンシャルの直接近似を回避でき、H¹適合有限要素法による正しく収束する解が得られる。
• 本手法で用いられる分離型・線形化有限要素法は、秩序パラメータでO(h¹·⁰⁹)、|ψ|でO(h¹·⁰⁷)、磁気ポテンシャルでO(h⁰·⁸⁷)の収束率を達成しており、定理3で示された理論的最悪ケースのO(h²·³)を上回る。
• 数値結果により、本手法は正確な解に収束するが、標準FEMはh = 1/256までメッシュを細かくしても収束しないことが確認された。
• 本手法は非凸な曲線境界を持つ多角形に対しても有効であり、解析は幾何学的形状の変更を伴わず、本質的に同様に適用可能である。