A wave propagation method for three-dimensional hyperbolic conservation laws

A class of wave propagation algorithms for three-dimensional conservation laws and other hyperbolic systems is developed. These unsplit finite-volume methods are based on solving one-dimensional Riemann problems at the cell interfaces and applying flux-limiter functions to suppress oscillations arising from second-derivative terms. Waves emanating from the Riemann problem are further split by solving Riemann problems in the transverse directions to model cross-derivative terms. With proper upwinding, a method that is stable for Courant numbers up to 1 can be developed. The stability theory for three-dimensional algorithms is found to be more subtle than in two dimensions and is studied in detail. In particular we find that some methods which are unconditionally unstable when no limiter is applied are (apparently) stabilized by the limiter function and produce good looking results. Several computations using the Euler equations are presented including blast wave and complex shock/vorticity problems. These algorithms are implemented in the CLAWPACK software, which is freely available.

[1]  P. L. Roe,et al.  Optimum positive linear schemes for advection in two and three dimensions , 1992 .

[2]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[3]  E. Toro.  Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[4]  Philip L. Roe,et al.  Discontinuous solutions to hyperbolic systems under operator splitting , 1987 .

[5]  B. Vanleer,et al.  Towards the ultimate conservative difference scheme. IV. A new approach to numerical convection , 1977 .

[6]  Wai Sun Don,et al.  Numerical simulation of shock-cylinder interactions I.: resolution , 1995 .

[7]  R. F. Warming,et al.  Upwind second-order difference schemes and applications in unsteady aerodynamic flows , 1975 .

[8]  High-resolution finite volume methods for acoustics in a rapidly-varying heterogeneous medium , 1997 .

[9]  P. Lax.  Hyperbolic systems of conservation laws II , 1957 .

[10]  T. J. Barth,et al.  On unstructured grids and solvers , 1990 .

[11]  R. LeVeque.  High-resolution conservative algorithms for advection in incompressible flow , 1996 .

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

[13]  S. Zalesak.  Fully multidimensional flux-corrected transport algorithms for fluids , 1979 .

[14]  A. Bressan.  Hyperbolic Systems of Conservation Laws , 1999 .

[15]  P. Sweby.  High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[16]  R. LeVeque.  Wave Propagation Algorithms for Multidimensional Hyperbolic Systems , 1997 .

[17]  G. R. Shubin,et al.  An unsplit, higher order Godunov method for scalar conservation laws in multiple dimensions , 1988 .

[18]  Clinton P. T. Groth,et al.  Assessment of Riemann solvers for unsteady one-dimensional inviscid flows for perfect gases , 1988 .

[19]  J. Boris,et al.  Flux-Corrected Transport , 1997 .

[20]  Randall J. LeVeque,et al.  High-resolution finite-volume methods for acoustic waves in periodic and random media , 1999 .

[21]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

[22]  P. Colella.  Multidimensional upwind methods for hyperbolic conservation laws , 1990 .

[23]  I. Bohachevsky,et al.  Finite difference method for numerical computation of discontinuous solutions of the equations of fluid dynamics , 1959 .

[24]  S. J. Billett.  A class of upwind methods for conservation laws , 1994 .

[25]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[26]  Arthur Veldman,et al.  NUMERICAL METHODS FOR FLUID DYNAMICS 4 , 1993 .

[27]  P. Raviart,et al.  Numerical Approximation of Hyperbolic Systems of Conservation Laws , 1996, Applied Mathematical Sciences.

[28]  James J. Quirk,et al.  On the dynamics of a shock–bubble interaction , 1994, Journal of Fluid Mechanics.

[29]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[30]  J. Fromm.  A method for reducing dispersion in convective difference schemes , 1968 .

[31]  R. LeVeque.  Numerical methods for conservation laws , 1990 .

[32]  R. LeVeque,et al.  Adaptive Mesh Refinement Using Wave-Propagation Algorithms for Hyperbolic Systems , 1998 .

[33]  Philip L. Roe,et al.  Sonic Flux Formulae , 1992, SIAM J. Sci. Comput..

[34]  R. LeVeque,et al.  Three-dimensional Euler Computationsusing CLAWPACKJan , 1995 .

[35]  W. Grabowski,et al.  The multidimensional positive definite advection transport algorithm: nonoscillatory option , 1990 .

[36]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[37]  M. Crandall,et al.  The method of fractional steps for conservation laws , 1980 .

[38]  George Turkiyyah,et al.  ZPLCLAW: A Parallel Portable Toolkit for Wave Propagation Problems , 2000 .

[39]  B. P. Leonard,et al.  Positivity-preserving numerical schemes for multidimensional advection , 1993 .

[40]  P. Lax,et al.  Difference schemes for hyperbolic equations with high order of accuracy , 1964 .

[41]  Jay P. Boris,et al.  Flux-corrected transport. I. SHASTA, a fluid transport algorithm that works , 1973 .

[42]  J. Saltzman,et al.  An unsplit 3D upwind method for hyperbolic conservation laws , 1994 .

[43]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme III. Upstream-centered finite-difference schemes for ideal compressible flow , 1977 .

[44]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[45]  Lawrence Snyder,et al.  A programmer's guide to ZPL , 1999 .

[46]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .