Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
10357258 | Journal of Computational Physics | 2005 | 23 Pages |
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/Ît) 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.
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Yingjie Liu,