差分格式应用于非均匀网格的精度保持算法

Precision Preservation Algorithm for Difference Schemes Applied to Non-Uniform Grids

  • 摘要: 为解决有限差分格式在非均匀网格下精度下降的问题,本文通过研究发现一维非均匀网格中几何诱导误差的产生机制,构建仅含一维非均匀性的网格模型,理论分析与数值实验揭示了传统格式精度损失的机理。理论研究发现:现有差分格式如一阶迎风、MUSCL、WENO等,在非均匀网格下难以达到理论设计精度,其根源在于坐标变换引入的高阶导数项对截断误差的不可控影响。基于有限差分法对导数逼近的本质目标,创新性提出精度保持算法。该方法通过加权融合不同方向差商项,修正传统格式在非均匀网格上的离散策略,以重构导数逼近模型来有效消除几何诱导误差。通过线性流与平方流为初场数值验证表明,数值验证表明:改进后的一阶格式在非均匀网格中可严格保持线性流场(误差量级10–17),二阶格式进一步保持平方流场(误差量级10–16),本研究结果为后续非均匀网格误差消减算法的完善提供了新的思路。

     

    Abstract: To address the issue of accuracy degradation in finite difference schemes on non-uniform grids, this study investigated the mechanism behind geometry-induced errors in one-dimensional non-uniform grid configurations. A simplified grid model with only one-dimensional non-uniformity was constructed, and both theoretical analysis and numerical experiments were conducted to reveal the root causes of accuracy loss in traditional schemes. Theoretical findings indicated that existing difference schemes, such as first-order upwind, MUSCL, and WENO, failed to achieve their designed order of accuracy on non-uniform grids due to the uncontrolled influence of higher-order derivative terms introduced by coordinate transformations on the truncation error. Based on a re-examination of the fundamental goal of derivative approximation in finite difference methods, an innovative accuracy-preserving algorithm was proposed. This method reconstructed the derivative approximation model by combining directional difference quotients in a weighted manner, thereby correcting the discretization strategy of traditional schemes and effectively eliminating geometry-induced errors. Numerical validations using linear and quadratic flow fields as initial conditions demonstrated that the improved first-order scheme can strictly preserve linear flow fields (with error magnitude on the order of 10–17), and the second-order scheme can further preserve quadratic flow fields (with error magnitude on the order of 10–16). The findings of this study offer new insights for the development of error-reduction algorithms on non-uniform grids.

     

/

返回文章
返回