计算物理 ›› 2005, Vol. 22 ›› Issue (1): 13-18.

• 研究论文 • 上一篇    下一篇

双曲型守恒律的一种高精度TVD差分格式

郑华盛1,2, 赵宁1   

  1. 1. 南京航空航天大学空气动力学系, 江苏 南京 210016;
    2. 南昌航空工业学院信息与计算科学系, 江西 南昌 330034
  • 收稿日期:2003-12-17 修回日期:2004-06-27 出版日期:2005-01-25 发布日期:2005-01-25
  • 作者简介:郑华盛(1966-),男,江西景德镇,教授,博士,从事计算流体力学的研究.
  • 基金资助:
    航空科学基金(01A52003;02A52004);国防预研资助项目

A High Order Accurate TVD Difference Scheme for Hyperbolic Conservation Laws

ZHENG Hua-sheng1,2, ZHAO Ning1   

  1. 1. Department of Aerodynamics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China;
    2. Department of Information and Computational Science, Nanchang Institute of Aeronautical Technology, Nanchang 330034, China
  • Received:2003-12-17 Revised:2004-06-27 Online:2005-01-25 Published:2005-01-25

摘要: 构造了一维双曲型守恒律方程的一个高精度高分辨率的守恒型TVD差分格式.其主要思想是:首先将计算区域划分为互不重叠的小单元,且每个小单元再根据希望的精度阶数分为细小单元;其次,根据流动方向将通量分裂为正、负通量,并通过小单元上的高阶插值逼近得到了细小单元边界上的正、负数值通量,为避免由高阶插值产生的数值振荡,进一步根据流向对其进行TVD校正;再利用高阶Runge KuttaTVD离散方法对时间进行离散,得到了高阶全离散方法.进一步推广到一维方程组情形.最后对一维欧拉方程组计算了几个算例.

关键词: 双曲型守恒律, 通量分裂, 差分格式, 高阶精度, 欧拉方程组

Abstract: A high-order,high resolution,conservative TVD difference scheme is presented for one dimensional hyperbolic conservation equations.The basic idea is as follows.Firstly,the computation domain is divided into many non-overlapping subdomains,and then each subdomain is further subdivided into small cells according to the required accuracy; Secondly,by the flow direction,flux splitting is introduced,and high-order approximation in the subdomain are used to compute the positive/negative numerical fluxes at cell boundaries.Furthermore,TVD corrections are considered to prevent oscillations near discontinuities from the high-order interpolation.Moreover,by means of high-order TVD Runge-Kutta time discretization,a high-order fully discretization method is obtained.The extension to one dimensional systems is also carried out.Finally,numerical experiments on one dimensional Euler equations are given,and numerical results are satisfactory.

Key words: hyperbolic conservation laws, flux splitting, TVD difference scheme, high order accuracy, Euler equations

中图分类号: