[论文解读] 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

一句话总结

本文在非凸多边形区域中建立了时变Ginzburg–Landau方程的适定性，其中磁势不具有H¹正则性，导致标准有限元方法产生错误的解。通过应用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.

研究动机与目标

• 在具有凹角的非凸多边形区域中建立时变Ginzburg–Landau系统的适定性，标准有限元方法因磁势奇点而失效。
• 解决当磁势不属于H¹空间时传统有限元方法失效的问题，这在具有钝角的区域中是常见现象。
• 利用Hodge分解重构系统，避免直接计算奇异磁势，确保关键未知量具有H¹正则性。
• 设计一种基于重构系统正则性的解耦且线性化的有限元格式，并给出严格的误差估计。
• 在L形区域中通过数值实验验证该方法，证明其在标准FEM失效时仍能实现收敛。

提出的方法

• 通过Hodge分解重写原始的时变Ginzburg–Landau系统，将磁势分解为正则部分与奇异部分，实现解耦。
• 构建一个等价的椭圆与抛物型方程组，仅涉及H¹-符合的未知量，避免对非H¹磁势的直接逼近。
• 对重构系统应用解耦且线性化的有限元方法，实现稳定且收敛的数值解。
• 采用向后欧拉时间离散化方法，并结合分片线性有限元进行空间逼近。
• 基于重构系统未知量的正则性，证明最优阶误差估计，这些未知量已被证明属于H¹空间。
• 在具有精确解的L形区域上实施数值实验，比较标准FEM与所提方法的收敛行为。

实验结果

研究问题

• RQ1 在具有钝角的非凸多边形区域中，即使磁势不具有H¹正则性，时变Ginzburg–Landau系统是否仍存在唯一解？
• RQ2 当磁势不属于H¹空间时，标准有限元方法是否仍能产生收敛解，特别是在具有角奇点的区域中？
• RQ3 是否可以将Ginzburg–Landau方程重写为一个等价系统，避免直接计算奇异磁势，并确保未知量的H¹正则性？
• 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的精细网格下也无法收敛。
• 该方法在非凸曲线多边形中依然有效，分析可无需本质修改地推广至此类几何形状。