In this paper we present high-order I-stable centered difference schemes for the numerical simulation of viscous compressible flows. Here I-stability refers to time discretizations whose linear stability regions conta...In this paper we present high-order I-stable centered difference schemes for the numerical simulation of viscous compressible flows. Here I-stability refers to time discretizations whose linear stability regions contain part of the imaginary axis. This class of schemes has a numerical stability independent of the cell-Reynolds number Re, thus allows one to simulate high Reynolds number flows with relatively larger Re, or coarser grids for a fixed Re. On the other hand, Re cannot be arbitrarily large if one tries to obtain adequate numerical resolution of the viscous behavior. We investigate the behavior of high-order I-stable schemes for Burgers' equation and the compressible Navier-Stokes equations. We demonstrate that, for the second order scheme, Re ≤ 3 is an appropriate constraint for numerical resolution of the viscous profile, while for the fourth-order schemes the constraint can be relaxed to Re ≤ 6.0ur study indicates that the fourth order scheme is preferable: better accuracy, higher resolution, and larger cell-Reynolds numbers.展开更多
基金Research supported by the National University of Singapore grant No. R-151-000-016-112. Email address: bao@cz3.nus.edu.sg.
文摘In this paper we present high-order I-stable centered difference schemes for the numerical simulation of viscous compressible flows. Here I-stability refers to time discretizations whose linear stability regions contain part of the imaginary axis. This class of schemes has a numerical stability independent of the cell-Reynolds number Re, thus allows one to simulate high Reynolds number flows with relatively larger Re, or coarser grids for a fixed Re. On the other hand, Re cannot be arbitrarily large if one tries to obtain adequate numerical resolution of the viscous behavior. We investigate the behavior of high-order I-stable schemes for Burgers' equation and the compressible Navier-Stokes equations. We demonstrate that, for the second order scheme, Re ≤ 3 is an appropriate constraint for numerical resolution of the viscous profile, while for the fourth-order schemes the constraint can be relaxed to Re ≤ 6.0ur study indicates that the fourth order scheme is preferable: better accuracy, higher resolution, and larger cell-Reynolds numbers.
