[论文解读] An RBF-FD polynomial method based on polyharmonic splines for the Navier-Stokes equations: Comparisons on different node layouts

Gregory A. Barnett, Natasha Flyer+1 more|arXiv (Cornell University)|Sep 9, 2015

Computational Fluid Dynamics and Aerodynamics参考文献 22被引用 7

一句话总结

本文提出了一种新颖的RBF-FD方法，结合多焦样条（PHS）与高阶多项式，用于在各种节点布局上求解Navier-Stokes方程，实现无需调节形状参数的高阶精度。通过消除停滞误差并避免边界振荡，该方法在六边形和散乱节点布局上表现出卓越的稳定性和准确性。

ABSTRACT

Polyharmonic spline (PHS) radial basis functions (RBFs) are used together with polynomials to create local RBF-finite-difference (RBF-FD) weights on different node layouts for spatial discretization of the compressible Navier-Stokes equations at low Mach number, relevant to atmospheric flows. Test cases are taken from the numerical weather prediction community and solved on bounded domains. Thus, attention is given on how to handle boundaries with the RBF-FD method, as well as a novel implementation for the presented approach. Comparisons are done on Cartesian, hexagonal, and quasi-uniformly scattered node layouts. Since RBFs are independent of a coordinate system (and only depend on the distance between nodes), changing the node layout amounts to changing one line of code. In addition, consideration and guidelines are given on PHS order, polynomial degree and stencil size. The main advantages of the present method are: 1) capturing the basic physics of the problem surprisingly well, even at very coarse resolutions, 2) high-order accuracy without the need of tuning a shape parameter, and 3) the inclusion of polynomials eliminates stagnation (saturation) errors.

研究动机与目标

• 开发一种基于RBF-FD的高阶、无网格空间离散化方法，用于求解Navier-Stokes方程，采用多焦样条RBF与多项式。
• 通过将PHS RBF与高次多项式结合，消除传统RBF-FD方法中常见的停滞（饱和）误差。
• 评估不同节点布局（笛卡尔、六边形与散乱）对大气流模拟中解的精度、稳定性及涡旋结构的影响。
• 为RBF-FD在流体动力学问题中选择PHS阶次、多项式次数与带宽大小提供指导。
• 证明该方法在不同雷诺数与边界条件下均具有鲁棒性，且无需调整超粘性参数。

提出的方法

• 在最小二乘框架下，利用多焦样条RBF（ϕ(r) = r^{2m+1}）与高阶多项式（次数≥4）构建局部RBF-FD微分权重。
• 采用对称带宽策略，通过选择合适的节点集合确保数值稳定性和高阶收敛性，2D问题中使用13-和37-节点带宽。
• 提出一种新型边界处理策略，有效避免虚假振荡，尤其在非笛卡尔布局中表现优异。
• 采用超粘性（γ = 2^{-6}h^{-2k}）稳定模拟过程，且在雷诺数变化时无需调整参数。
• 采用与节点布局无关的代码结构：仅需修改一行代码即可从笛卡尔节点切换至六边形或散乱节点。
• 通过数值天气预报中的基准测试案例（如上升热泡、剪切流）验证不同布局下的精度与稳定性。

实验结果

研究问题

• RQ1 在RBF-FD对Navier-Stokes方程的离散化中，将高阶多项式与多焦样条RBF结合，如何影响停滞误差？
• RQ2 不同节点布局（笛卡尔、六边形与散乱）如何影响Navier-Stokes解的精度、稳定性及定性特征？
• RQ3 所提方法是否可在不调节形状参数的情况下实现高阶收敛？PHS阶次如何影响收敛速率？
• RQ4 在湍流区域中，该方法在不同节点分布下对细尺度涡旋与不稳定性解析的表现如何？
• RQ5 该方法是否可在不调整超粘性或时间步长的前提下，保持宽范围雷诺数下的稳定性和准确性？

主要发现

• 在RBF-FD基中引入多项式可有效消除停滞（饱和）误差，实现在粗网格下的稳定收敛。
• 在六边形节点布局下，该方法产生最精确且物理上更合理的解，尤其在边界附近振荡最小，优于笛卡尔节点。
• 散乱节点布局的性能与笛卡尔布局相当或更优，未出现精度或收敛性下降，且在解析细尺度涡旋方面优于笛卡尔布局。
• 对于2D问题，37-节点带宽可支持最高四次多项式（15个基函数），且为获得最佳稳定性，带宽大小应约为多项式基函数数量的两倍。
• 收敛速率主要由多项式次数决定，而非PHS阶次；更高PHS阶次仅略微改善收敛常数。
• 该方法在雷诺数增加六个数量级的范围内保持稳定与准确，且无需调整超粘性或时间步长，展现出在湍流区域的强鲁棒性。