华南理工大学学报(自然科学版) ›› 2009, Vol. 37 ›› Issue (11): 151-156.

• 机械工程 • 上一篇    

上风式LRPIM求解对流占优问题的影响因素分析

张祖军 阮锋 刘毅 刘成军   

  1. 华南理工大学 机械与汽车工程学院, 广东 广州 510640
  • 收稿日期:2008-11-25 修回日期:2009-03-26 出版日期:2009-11-25 发布日期:2009-11-25
  • 通信作者: 张祖军(1980-),男,博士生,主要从事模塑成型及其计算机技术研究. E-mail:xz800110@yahoo.cn
  • 作者简介:张祖军(1980-),男,博士生,主要从事模塑成型及其计算机技术研究.
  • 基金资助:

    广东省教育部产学研结合项目(2007A090302028)

Analysis of Factors Influencing the Solving of Convection-Dominated Problems by Means of Upwinding LRPIM

Zhang Zu-jun  Ruan Feng  Liu Yi  Liu Cheng-jun   

  1. School of Mechanical and Automotive Engineering, South China University of Technology, Guangzhou 510640, Guangdong, China
  • Received:2008-11-25 Revised:2009-03-26 Online:2009-11-25 Published:2009-11-25
  • Contact: 张祖军(1980-),男,博士生,主要从事模塑成型及其计算机技术研究. E-mail:xz800110@yahoo.cn
  • About author:张祖军(1980-),男,博士生,主要从事模塑成型及其计算机技术研究.
  • Supported by:

    广东省教育部产学研结合项目(2007A090302028)

摘要: 为提高对流扩散问题中对流项占优时的数值稳定性,采用局部上风式权函数,在局部径向基点插值无网格法(LRPIM)基础上构建了上风式局部径向基点插值无网格法(ULRPIM).将其与LRPIM方法的对流占优问题计算结果进行比较,发现ULRPIM能够得到更好的结果.通过对比不同参数时ULRPIM计算结果的相对误差,研究了各个参数对数值结果稳定性的影响规律.结果表明:权函数的偏置量需要随着对流强弱而变化,径向基函数(RBF)参数q、支持域尺寸对计算结果的影响比较大,RBF参数ac、Gauss积分子域数对计算结果的影响相对较小.

关键词: 对流占优, 上风式, 局部径向基点插值法, 稳定性

Abstract:

In order to improve the numerical stability during the solving of convection-dominated convection-diffusion problems, an upwinding local radial point interpolation method (ULRPIM) is proposed based on the local radial point interpolation method (LRPIM) and the local upwinding weight function. The proposed method is then compared with the LRPIM through some numerical calculations, finding that ULRPIM achieves better results. Mo- reover, the influences of different parameters on the stability of calculation results are investigated by comparing the relative errors of ULRPIM results at different parameters. It is found that the offset of weight function should change with the Peclet number, and that both the parameter q of radial basis functions (RBF) and the size of the support domain have a greater influence on the calculation results than the RBF parameter ac and the quantity of Gauss integral points.

Key words: convection domination, upwinding, local radial point interpolation method, stability