Central schemes on overlapping cells

  • Authors:
  • Yingjie Liu

  • Affiliations:
  • School of Mathematics, Georgia Institute of Technology, Atlanta, GA 30332, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2005

Quantified Score

Hi-index 31.47

Visualization

Abstract

Nessyahu and Tadmor's central scheme [J. Comput. Phys. 87 (1990)] has the benefit of not using Riemann solvers for solving hyperbolic conservation laws. But the staggered averaging causes large dissipation when the time step size is small compared to the mesh size. The recent work of Kurganov and Tadmor [J. Comput. Phys. 160 (2000)] overcomes this problem by using a variable control volume and results in semi-discrete and fully discrete non-staggered schemes. Motivated by this work, we introduce overlapping cell averages of the solution at the same discrete time level, and develop a simple alternative technique to control the O(1/@Dt) dependence of the dissipation. The semi-discrete form of the central scheme can also be obtained to which the TVD Runge-Kutta time discretization methods of Shu and Osher [J. Comput. Phys. 77 (1988)] or other stable and sufficiently accurate ODE solvers can be applied. This technique is essentially independent of the reconstruction and the shape of the mesh. The overlapping cell representation of the solution also opens new possibilities for reconstructions. Generally speaking, more compact reconstruction can be achieved. In the following, schemes of up to fifth order in 1D and third order in 2D have been developed. We demonstrate through numerical examples that by combining two classes of the overlapping cells in the reconstruction we can achieve higher resolution.