紧凑型LABFM：一种具有谱类分辨能力的无网格方法框架

《Journal of Computational Physics》：Compact LABFM: a framework for meshless methods with spectral-like resolving power

　　

无网格方法常用于偏微分方程（PDEs）系统的数值模拟，特别是涉及复杂几何或自由表面的问题。研究人员提出了一种基于局部各向异性基函数方法（LABFM）的新型紧凑格式，这是一种无网格方法，能够提供具有任意多项式一致性的空间算子近似。该方法通过使用隐式模板来优化每个算子的分辨能力，同时保持所得全局稀疏线性系统的对角占优性，从而模仿紧凑有限差分法。通过一系列收敛性测试和分辨能力分析，证明了新方法能提供改进的近似，随后计算了典型PDEs的解。观察到精度的显著提升，特别是对于包含高波数分量的解。研究人员的紧凑无网格方法为在复杂几何中实现具有谱类分辨能力的高阶PDEs模拟提供了途径，并有可能导致此类问题数值解精度的阶跃性变化。

该论文提出了一种基于局部各向异性基函数方法（LABFM）的紧凑无网格格式，旨在通过隐式离散化优化空间算子的分辨能力，从而在复杂几何中实现具有谱类精度的偏微分方程（PDEs）模拟。研究背景方面，传统有限差分法虽能通过紧凑格式（如Lele型）提升分辨能力，但受限于简单几何结构；而无网格方法（如SPH、RBF-FDs）虽可处理复杂边界，但存在收敛阶低或谱分辨率损失的问题。现有紧凑无网格方法（如RBF-HFDs、CMLS）尚未针对分辨能力进行系统优化。为此，研究人员开发了紧凑LABFM，其核心思想是在每个节点处构建隐式模板，通过优化系数最大化算子对宽波段波数的分辨能力，同时保持全局稀疏线性系统的对角占优性。研究通过收敛性测试、二维傅里叶分辨能力分析，以及求解粘性Burgers方程、泊松方程和Navier-Stokes泰勒-格林涡等典型问题，验证了该方法在精度上的显著提升，误差减少可达一个数量级，尤其适用于高波数分量。该工作发表在《Journal of Computational Physics》，为复杂几何中高精度谱类分辨率模拟开辟了新途径。

关键技术方法方面，主要包括：1. 基于LABFM的显式格式，使用Hermite多项式与Wendland C²平滑核函数构造各向异性基函数，通过求解局部线性系统获得权重，保证多项式一致性；2. 隐式模板选择，对于梯度算子，根据另一方向坐标偏差选择节点，对于拉普拉斯算子，取梯度模板的并集，最大节点数分别为9和17（含中心节点）；3. 系数优化，采用高斯函数形式设定隐式系数，通过迭代减小参数以最大化分辨能力，同时确保对角占优性和条件数良好，并通过二维傅里叶分析评估有效波数；4. 全局线性系统求解，使用GMRES算法和AMG预条件子（HYPRE库）进行稀疏迭代求解。

研究结果部分：

§4 分辨能力：通过二维傅里叶分析，比较了紧凑与显式LABFM的梯度算子和拉普拉斯算子的有效波数。结果显示，紧凑格式在宽波数范围内更接近精确微分，特别是拉普拉斯算子具有谱类精度。误差度量ε₁（梯度）和ε₂（拉普拉斯）表明，紧凑格式的误差更小，且较大的隐式模板通常带来更好的分辨能力。但紧凑格式的虚部误差略大，但影响较小。

§5 收敛性和稳定性测试：使用测试函数（傅里叶级数表示顶部函数）比较了梯度算子和拉普拉斯算子的L₂范数误差。紧凑格式在几乎所有分辨率下误差减少，二阶格式误差减少达80%，四阶格式减少60%以上。稳定性分析表明，梯度算子紧凑格式比显式更不稳定，但拉普拉斯算子所有格式均稳定，紧凑格式的振荡模式衰减快。

§6 应用：a) 粘性Burgers方程：在Re=100时，紧凑格式在平流主导阶段误差显著降低，二阶格式误差减少一个数量级，四阶格式减少80%。b) 泊松方程：在复杂几何（带孔方域和环形域）中，紧凑格式误差减少最高达一个数量级，且计算成本增加极小（仅需额外右端项计算）。c) 泰勒-格林涡：在固定网格上改变波数，紧凑格式误差减少最高达94%，且随波数增加改善持续。

讨论部分指出，该方法仅适用于光滑解，且系数形式和隐式模板选择可能非最优，但可扩展至其他无网格方法（如RBF-HFDs）。研究结论为：研究人员提出了一种新的紧凑无网格方法，基于Lele型格式优化分辨能力，在宽波数范围内显著改善，拉普拉斯算子具有谱类精度。在典型PDEs中误差减少达一个数量级。该方法可扩展至其他无网格方法，为下一代高精度无网格模拟提供基础，但仅适用于光滑解。