Crank–Nicolson跳跃格式在相场晶体模型中的二阶精度、无条件能量稳定性与误差估计

《Computers & Mathematics with Applications》：Stability and error estimate of the second-order Crank–Nicolson leap-frog scheme for the phase field crystal model

本文提出并严谨分析了一种针对相场晶体(PFC)模型的二阶Crank–Nicolson跳跃(CNLF)全离散数值格式。该工作证明了其满足质量守恒、无条件能量稳定（无条件能量稳定）和线性可解性，获得了时间上的二阶收敛精度，并进行了详细的误差估计。数值实验验证了其在长时间相粗化动力学模拟中的高效性与准确性。

　　

Section snippets

Numerical schemes

为便于后续分析，我们首先明确本文中使用的符号约定。Lebesgue空间L^p(Ω)与标准的L^p-范数相关联，定义为∥u∥_L^p:= (∫_Ω|u(x)|^pdx)^1/p， 对于 1 ≤ p < ∞。对于p = ∞的情况，我们使用本质确界范数，并定义空间L^∞(Ω)及相应的范数：∥v∥_L^∞:= ess sup_x∈Ω|v(x)|。特殊情况L²(Ω)与内积(·, ·)和范数‖·‖相关联。对于Sobolev空间，我们用……

Error analysis

本节在精确解具有温和正则性假设的前提下，推导了全离散格式2.2的严格误差界。

定义H_per²(Ω)-正交投影算子P_N: H_per²(Ω) → V_N，使得对于所有v ∈ H_per²(Ω)，

((Δ+β)(v - P_Nv), (Δ+β)φ_N)_Ω= 0 ? φ_N∈ V_N，

且满足∫_Ω(v - P_Nv) dx = 0。

关于这个投影算子的误差估计可以从定义域投影算子的标准结果中获得（见[27]），并可利用p-方法的估计技术得到（参见[28]、[29]）。

Lemma 5

对于所有v ∈ H_per^s(Ω) 满足l ≤……

Numerical results

在本节中，我们使用稳定化后的格式（Scheme 2.2）对PFC模型进行数值模拟，展示了所提出的Crank–Nicolson跳跃(CNLF)方法的精度、效率与能量稳定性。空间离散化采用Fourier谱方法，所有计算均通过快速傅里叶变换(FFT)加速。模拟在一个受周期性边界条件约束的方形区域[0, L]²上进行。在所有数值模拟中……

Conclusion

在这项工作中，我们为相场晶体(PFC)模型开发了一种高效、全离散的数值格式，实现了时间上的二阶精度和空间上的谱精度。所提出的方法将二阶Crank–Nicolson跳跃时间离散与Fourier伪谱空间近似相结合，确保了高精度和计算效率。我们的数值格式保证了三个基本理论性质：质量守恒、无条件能量……