Composite Schemes For Conservation Laws

Global composition of several time steps of the two-step Lax--Wendroff scheme followed by a Lax--Friedrichs step seems to enhance the best features of both, although it is only first order accurate. We show this by means of some examples of one-dimensional shallow water flow over an obstacle. In two dimensions we present a new version of Lax--Friedrichs and an associated second order predictor-corrector method. Composition of these schemes is shown to be effective and efficient for some two-dimensional Riemann problems and for Noh's infinite strength cylindrical shock problem. We also show comparable results for composition of the predictor-corrector scheme with a modified second order accurate weighted essentially nonoscillatory (WENO) scheme. That composition is second order but is more efficient and has better symmetry properties than WENO alone. For scalar advection in two dimensions the optimal stability of the new predictor-corrector scheme is shown using computer algebra. We also show that the generalization of this scheme to three dimensions is unstable, but by using sampling we are able to show that the composites are suboptimally stable.

[1]  Xu-Dong Liu,et al.  Solution of Two-Dimensional Riemann Problems of Gas Dynamics by Positive Schemes , 1998, SIAM J. Sci. Comput..

[2]  A. Seidenberg.  A NEW DECISION METHOD FOR ELEMENTARY ALGEBRA , 1954 .

[3]  Hoon Hong,et al.  Improvements in cad-based quantifier elimination , 1990 .

[4]  D. Gottlieb,et al.  Numerical stabilizers and computing time for second-order accurate schemes☆ , 1972 .

[5]  W. F. Noh.  Errors for calculations of strong shocks using an artificial viscosity and artificial heat flux , 1985 .

[6]  David L. Book,et al.  Flux-corrected transport II: Generalizations of the method , 1975 .

[7]  Bernd Einfeld.  On Godunov-type methods for gas dynamics , 1988 .

[8]  Well-Posed Problems and Stable Difference Operators , 1968 .

[9]  A. Tarski.  A Decision Method for Elementary Algebra and Geometry , 2023 .

[10]  A. Y. Le Roux,et al.  Convergence of an Antidiffusion Lagrange-Euler Scheme for Quasilinear Equations , 1984 .

[11]  A. Leroux,et al.  A new version of the two-dimensional Lax-Friedrichs scheme , 1994 .

[12]  Richard Liska,et al.  FIDE: a REDUCE package for automation of FInite difference method for solving pDE , 1990, ISSAC '90.

[13]  R. D. Richtmyer,et al.  Difference methods for initial-value problems , 1959 .

[14]  A. Harten,et al.  Self-adjusting hybrid schemes for shock computations , 1972 .

[15]  A. Kasahara,et al.  Nonlinear shallow fluid flow over an isolated ridge , 1968 .

[16]  James P. Collins,et al.  Numerical Solution of the Riemann Problem for Two-Dimensional Gas Dynamics , 1993, SIAM J. Sci. Comput..

[17]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[18]  Per Lötstedt,et al.  Nonlinear filters for efficient shock computation , 1989 .