[论文解读] A vertex-weighted-Least-Squares gradient reconstruction

Fan Zhang|arXiv (Cornell University)|Feb 15, 2017

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

一句话总结

本文提出了一种针对单元中心有限体积法的顶点加权最小二乘（VWLSQ）梯度重构方法，通过加权最小二乘法计算顶点梯度并将其平均以获得单元中心梯度，从而在非结构化三角形网格上提升精度与收敛性。VWLSQ(1)方案在各种网格条件下，尤其是畸变网格上，展现出更优的精度与鲁棒性，并可通过复用顶点梯度实现高效的黏性通量计算，无需额外的梯度重构。

ABSTRACT

Gradient reconstruction is a key process for the spatial accuracy and robustness of finite volume method, especially in industrial aerodynamic applications in which grid quality affects reconstruction methods significantly. A novel gradient reconstruction method for cell-centered finite volume scheme is introduced. This method is composed of two successive steps. First, a vertex-based weighted-least-squares procedure is implemented to calculate vertex gradients, and then the cell-centered gradients are calculated by an arithmetic averaging procedure. By using these two procedures, extended stencils are implemented in the calculations, and the accuracy of gradient reconstruction is improved by the weighting procedure. In the given test cases, the proposed method is showing improvement on both the accuracy and convergence. Furthermore, the method could be extended to the calculation of viscous fluxes.

研究动机与目标

• 提升单元中心有限体积方法在包含劣质单元的非结构化三角形网格上的梯度重构精度与收敛性。
• 解决传统节点平均格林-高斯（NA-GG）与基于单元的加权最小二乘（WLSQ）方案在处理网格偏斜与畸变时的局限性。
• 开发一种可复用顶点梯度以同时计算无黏与黏性通量的方法，提升计算效率。
• 在高长宽比网格上保持鲁棒性，同时消除对裁剪程序的需求。
• 为工业空气动力学模拟中现有的梯度重构技术提供更精确、更高效的替代方案。

提出的方法

• 通过基于顶点的加权最小二乘（VWLSQ）过程计算顶点梯度，求解由单元中心变量及其到顶点距离构成的超定方程组。
• 针对每个顶点，利用邻近单元中心值及其向量距离构建最小二乘系统，权重与距离成反比，以提升稳定性和精度。
• 通过计算单元所有邻近顶点的重构梯度的算术平均值，获得单元中心梯度。
• 加权方案通过强调距离更近、更相关的单元，提升最小二乘拟合的精度，降低对网格畸变的敏感性。
• 该方法避免了裁剪程序，并在高长宽比网格上保持鲁棒性。
• 通过复用相同顶点梯度，将该方法扩展至黏性通量计算，从而降低计算开销。

实验结果

研究问题

• RQ1 所提出的顶点加权最小二乘梯度重构方法在非结构化三角形网格上，与现有节点平均格林-高斯（NA-GG）和基于单元的WLSQ方案相比，在精度与收敛性方面表现如何？
• RQ2 VWLSQ(1)方案是否能在不使用裁剪程序的情况下，保持在高长宽比和畸变三角形网格上的高精度与鲁棒性？
• RQ3 复用顶点梯度同时用于无黏与黏性通量计算，是否相比传统方法提升了计算效率？
• RQ4 模板扩展与基于顶点的重构对非结构化有限体积方法中梯度精度的影响如何？
• RQ5 在数值耗散与熵产生方面，VWLSQ(1)方案与WLSQ(G)及其他先进方案相比表现如何？

主要发现

• 在所有测试方案中，VWLSQ(1)方案的数值阻力误差最低，NACA0012翼型案例的阻力系数为1.17×10⁻⁴。
• 在三角形网格上，VWLSQ(1)在精度与收敛性方面优于基于单元的WLSQ和节点平均GG方案，尤其在畸变与高长宽比区域表现更优。
• WLSQ(0)方案产生的阻力误差最小（2.26×10⁻⁵），但表现出显著的熵增加，表明存在高数值耗散，抵消了其精度优势。
• 熵分布显示，VWLSQ(1)的耗散低于WLSQ(0)，解释了其在精度与稳定性之间更优的平衡。
• VWLSQ(1)在精度与收敛性方面均优于WA-GG与PL-CLIP-GG，残差增长最小，且在混合网格上表现出稳定行为。
• 该方法在高长宽比网格上具有鲁棒性，无需裁剪，表明其在几何畸变下具备内在稳定性。