计算物理 ›› 2006, Vol. 23 ›› Issue (6): 717-720.

• 论文 • 上一篇    下一篇

柱坐标下Euler方程数值边界条件的一种处理方法

田保林   

  1. 北京应用物理与计算数学研究所, 北京 100088
  • 收稿日期:2005-06-24 修回日期:2006-01-17 出版日期:2006-11-25 发布日期:2006-11-25
  • 作者简介:田保林(1975-),男,山东成武,博士,从事计算流体力学方面的研究.
  • 基金资助:
    国家973项目(2005CB321700);中物院科技基金(20060648)资助项目

Numerical Boundary Condition of Euler Equations in Cylindrical Coordinate

Tian Bao-lin   

  1. Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
  • Received:2005-06-24 Revised:2006-01-17 Online:2006-11-25 Published:2006-11-25

摘要: 给出了柱坐标下Euler方程数值边界条件的一种处理方法.径向第一个网格点设在距离中心半点位置上.根据相应物理量的特性,在中心附近进行边界延拓,使得内点的高精度差分格式可以同样应用在网格中心附近,从而无需单侧差分格式,保持了一致的高阶精度.对于周向边界,也建立了一种周期延拓方法,使得在周向所有节点处都能够采用同样的高精度格式离散,并进行了数值试验.

关键词: 柱坐标, 极点奇异性, 高精度有限差分

Abstract: A method for the treatment of polar singularities of Euler equations is given. The first radial mesh point locates at a half-space away from the centerline. Based on the characteristics of the physical variables,the boundary near centerline is extended so that a high order finite difference scheme can be utilized as at inner mesh points.Similarly,in azimuthal direction the boundary is extended according to the periodicity.An uniform high-order precision is preserved during discretization of equations.

Key words: cylindrical coordinates, polar singularities, high-order finite difference schemes

中图分类号: