[论文解读] Interpolants of Lattice Functions for the Analysis of Atomistic/Continuum Multiscale Methods

Christoph Ortner, Alexander V. Shapeev|arXiv (Cornell University)|Apr 17, 2012

Advanced Mathematical Modeling in Engineering参考文献 10被引用 19

一句话总结

本文提出了一类适用于布拉维晶格上格点函数的（拟-）插值方法，实现了对离散原子数据的精确连续表示。该研究建立了关键的逼近性质——如节点插值的二阶精度和拟插值的三阶精度，这对分析原子/连续多尺度方法至关重要，并在索博列夫范数和勒贝格范数下给出了严格的误差界。

ABSTRACT

We introduce a general class of (quasi-)interpolants of functions defined on a Bravais lattice, and establish several technical results for these interpolants that are crucial ingredients in the analysis of atomistic models and atomistic/continuum multi-scale methods.

研究动机与目标

• 为多尺度建模中的离散格点函数开发连续插值的一般框架。
• 通过连续函数逼近离散原子数据，确保高阶精度，以支持数值分析。
• 在索博列夫空间和勒贝格空间中，为节点插值和拟插值建立严格的逼近误差估计。
• 通过稳定、可逆且具有局部紧支集的插值方法，支持原子/连续耦合方法的分析。
• 提供一个函数空间框架和技术工具，统一并拓展现有数值分析中关于多尺度方法的结果。

提出的方法

• 通过紧支集且属于 $W^{1,\infty}$ 的基函数 $\bar{\zeta}$ 定义节点插值 $\bar{v}(x) = \sum_{\xi \in \mathbb{Z}^d} v(\xi) \bar{\zeta}(x - \xi)$，并满足插值条件和多项式重构条件。
• 通过卷积引入拟插值 $\tilde{v}(x) = (\bar{\zeta} \ast \bar{v})(x)$，该方法保持更高正则性，并实现三阶精度。
• 对 $\bar{\zeta}$ 施加结构假设（Z1–Z4）：正则性、紧支集、仿射函数的多项式重构能力，以及在整数点的节点插值。
• 利用键密度引理和卷积恒等式，推导光滑函数被插值逼近的误差估计。
• 建立逆估计和局部范数等价性，将 $\|\nabla^j(\tilde{J}v - v)\|_{L^p}$ 用 $\|\nabla^k v\|_{L^p}$ 有界表示，其中 $\tilde{J}$ 是投影到拟插值空间的算子。
• 证明任意三次多项式均可精确表示为某个离散函数 $w \in \mathscr{U}$ 的拟插值 $\tilde{w}$，从而确保拟插值空间的一致性。

实验结果

研究问题

• RQ1 如何将离散格点函数一致且稳定地映射到连续函数，同时保持关键的物理和数学性质？
• RQ2 基函数 $\bar{\zeta}$ 需满足何种条件，才能保证节点插值的二阶精度和拟插值的三阶精度？
• RQ3 插值的逼近误差是否能在索博列夫范数下一致有界？其依赖关系如何随底层函数的正则性变化？
• RQ4 如何利用插值框架分析原子/连续多尺度方法中的收敛性和一致性？
• RQ5 在何种条件下插值操作是可逆的？这对多尺度耦合格式的稳定性有何影响？

主要发现

• 节点插值 $\bar{v}$ 属于 $W^{1,\infty}_{\rm loc}(\mathbb{R}^d; \mathbb{R}^m)$，确保局部利普希茨连续性和稳定性。
• 拟插值 $\tilde{v}$ 属于 $W^{3,\infty}_{\rm loc}(\mathbb{R}^d; \mathbb{R}^m)$，保证更高正则性，并对光滑函数实现三阶精度。
• 对任意 $v \in W^{k,p}_{\rm loc}(\mathbb{R}^d)$ 且 $k \leq 4$，误差 $\|\nabla^j(\tilde{J}v - v)\|_{L^p}$ 有界于 $C_{\tilde{J}} \|\nabla^k v\|_{L^p}$，证明了最优收敛速率。
• 任意三次多项式 $p$ 均可精确表示为某个 $w \in \mathscr{U}$ 的 $\tilde{w}$，确保拟插值空间的一致性。
• 插值算子 $\tilde{J}$ 是投影到拟插值空间的算子，其逼近误差由函数的 $k$ 阶索博列夫半范数控制。
• 当 $\bar{\zeta}$ 选择为 Q1 基函数（分段线性）时，$\bar{v}$ 对应张量积线性 B 样条，$\tilde{v}$ 对应三次 B 样条，验证了该框架的通用性及其与标准样条的联系。