A priori error analysis of a mass-lumped midpoint finite element method with a structure-preserving solver for the Landau--Lifshitz--Gilbert equation with Dzyaloshinskii--Moriya interaction

技术摘要：带有 DMI 的质量集中中点有限元法求解 LLG 方程的先验误差分析

1. 问题陈述
本文研究了用于近似描述在居里点温度以下铁磁材料中磁化动力学的 Landau–Lifshitz–Gilbert (LLG) 方程的数值解。具体而言，研究重点在于包含体相 Dzyaloshinskii–Moriya 相互作用 (DMI)，这在有效场中引入了一阶项，并诱导了自然的手性边界条件。主要的数学挑战在于逐点单位长度约束（$|m|=1$）、解空间的非凸几何特性，以及方程极强的非线性结构。虽然仅含交换能的模型已被广泛研究，但 DMI 的引入改变了变分结构和边界行为，使得构建稳定、准确且保持结构的有限元方案变得尤为精细。

2. 方法论
作者提出并分析了一种全离散的质量集中中点有限元方案。该方法包含以下组成部分：

• 离散化： 使用最低阶保形拉格朗日有限元 ($P_1$) 对空间域进行离散化。时间积分采用隐式中点规则。
• 质量集中 (Mass-Lumping)： 使用了一种简化积分（质量集中）内积。这项技术至关重要，因为它允许在每个节点上精确强制执行节点单位长度约束，而标准的 $L^2$ 内积在处理非线性方案时往往无法保持这一特性。
• 有效场算子： 离散有效场结合了交换刚度和体相 DMI。DMI 项通过一个尊重手性边界条件的离散双线性形式来处理。为了便于分析，引入了一个偏移有效场算子，以确保离散双线性形式的一致强制性。
• 求解器策略： 由于中点方案在每个时间步都会产生非线性代数系统，本文提出了一种保持结构的固定点迭代法。与通用求解器不同，该迭代旨在在每次迭代中都保持节点的单位长度约束。
• 误差分析框架： 分析将数值解与精确解的椭圆投影进行比较。证明过程依赖于中点法固有的几何抵消、对一阶 DMI 项的仔细估计，以及从离散能量定律中导出的稳定性估计。

3. 主要贡献
本文对计算微磁学领域做出了三个主要贡献：

1. 公式化与稳定性： 作者为带有交换能和体相 DMI 的 LLG 方程制定了一种全离散质量集中中点方案。他们证明了该方案能精确保持节点单位长度约束并满足离散能量恒等式，从而在不需要几何网格限制（如锐角条件）或严格 CFL 类条件的情况下，为稳定性奠定了基础。
2. 最优先验误差估计（理想方案）： 在满足适当的正则性假设下，本文证明了理想方案（假设精确求解非线性系统）具有最优阶收敛估计。误差在能量范数 ($H^1$) 下以 $O(h + k^2)$ 的阶数收敛，代表空间一阶精度和时间二阶精度。据作者所知，这是针对质量集中中点法（即使在没有 DMI 的情况下）的第一个此类先验误差估计。
3. 不精确求解器的分析： 意识到在实践中必须通过迭代求解非线性系统，本文分析了求解器误差的影响。作者证明，如果基于特定的残差停止准则终止固定点迭代，只要求解器容差相对于离散参数足够小，所得的不精确方案仍能保持与理想方案相同的收敛速率。

4. 结果

• 理论界限： 主要理论结果（定理 3.7 和定理 4.8）确立了能量范数下的误差受 $C(h + k^2)$ 以及一个取决于求解器容差的项所限制。具体而言，对于统一的求解器容差 $\varepsilon_{fp}$，条件 $\varepsilon_{fp} = O(h)$ 足以维持最优收敛速率。
• 数值实验：
  • 收敛速率： 在光滑解上的数值测试证实了理论收敛速率：对于空间细化，在 $H^1$ 中约为一阶，在 $L^2$ 中约为二阶；对于时间细化，为二阶。
  • 求解器容差： 改变固定点容差的实验表明，误差随 $\varepsilon_{fp}/h$ 线性缩放，验证了关于求解器诱导误差的理论预测。
  • 结构保持： 涉及 Bloch 型磁天体（skyrmion）和高梯度“泡状”（bubbling）实验的模拟证实，该方法能将节点单位长度约束保持至机器精度，并能正确捕捉耗散能量行为。
  • DMI 效应： 实验展示了 DMI 符号 ($\kappa_d$) 对动力学的非对称影响。负的 $\kappa_d$ 加速了梯度集中（泡状形成），而正的 $\kappa_d$ 则倾向于抑制集中并促进弛豫，这与相互作用的手性性质一致。

5. 重要性与主张
本文声称提供了第一个针对应用于带体相 DMI 的 LLG 方程的质量集中中点有限元法的严谨先验误差分析。其意义在于弥合了该方法已知的稳定性属性（精确约束保持和能量耗散）与定量收敛速率之间的鸿沟。

作者强调，其分析并不依赖于限制性的网格条件（如角度条件）或严苛的时间步长限制（这些限制常困扰其他保持结构的方案）。此外，通过证明其实际使用的不精确求解器不会降低收敛速率，这项工作为使用这种计算高效且保持结构的算法来模拟复杂的如磁天体等手性磁现象提供了严谨的依据。结果支持了该方法在处理交换能、DMI 与手性边界条件之间微妙相互作用时的鲁棒性。