New NAG library software for first-order partial differential equations

New NAG Fortran Library routines are described for the solution of systems of nonlinear, first-order, time-dependent partial differential equations in one space dimension, with scope for coupled ordinary differential or algebraic equations. The method-of-lines is used with spatial discretization by either the central-difference Keller box scheme or an upwind scheme for hyperbolic systems of conservation laws. The new routines have the same structure as existing library routines for the solution of second-order partial differential equations, and much of the existing library software is reused. Results are presented for several computational examples to show that the software provides physically realistic numerical solutions to a challenging class of problems.

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

[2]  Blair Swartz.  Courant-Like Conditions Limit Reasonable Mesh Refinement to Order $h^2 $ , 1987 .

[3]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[4]  J. Verwer,et al.  A numerical study of three moving-grid methods for one-dimensional partial differential equations which are based on the method of lines , 1990 .

[5]  Randall J. LeVeque,et al.  A study of numerical methods for hyperbolic conservation laws with stiff source terms , 1990 .

[6]  Marek Machura,et al.  A Survey of Software for Partial Differential Equations , 1980, TOMS.

[7]  H. C. Yee.  Construction of explicit and implicit symmetric TVD schemes and their applications. [Total Variation Diminishing for fluid dynamics computation] , 1987 .

[8]  M. Berzins,et al.  An adaptive theta method for the solution of stiff and nonstiff differential equations , 1992 .

[9]  H. B. Keller.  A New Difference Scheme for Parabolic Problems , 1971 .

[10]  D. Howell Peregrine,et al.  Low frequency waves in the surf zone , 1993 .

[11]  Martin Berzins,et al.  Developing software for time-dependent problems using the method of lines and differential-algebraic integrators , 1989 .

[12]  Bram van Leer,et al.  On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe , 1984 .

[13]  S. P. Spekreijse,et al.  Multiple grid and Osher''s scheme for the efficient solution of the steady Euler equations Applied N , 1986 .

[14]  Norman L. Schryer,et al.  Designing software for one-dimensional partial differential equations , 1990, TOMS.

[15]  A. Benazzouz,et al.  A C 1 interpolant for codes based on backward differentiation formulae , 1986 .

[16]  B. Vanleer,et al.  On the Relation Between the Upwind-Differencing Schemes of Godunov, Engquist–Osher and Roe , 1984 .

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

[18]  R. W. Brankin,et al.  The stiff integrators in the NAG library , 1988, SGNM.

[19]  M. Berzins,et al.  Developments in the NAG library software for parabolic equations , 1990 .

[20]  Martin Berzins,et al.  A Method for the Spatial Discretization of Parabolic Equations in One Space Variable , 1990, SIAM J. Sci. Comput..

[21]  J. C. Mason,et al.  Scientific software systems: based on the proceedings of the international symposium on scientific software and systems, held at the Royal Military College of Science, Shrivenham, July 1988 , 1990 .

[22]  Philip L. Roe,et al.  Upwind differencing schemes for hyperbolic conservation laws with source terms , 1987 .

[23]  Graeme Fairweather,et al.  A Box Method for a Nonlinear Equation of Population Dynamics , 1991 .

[24]  P. L. Roe,et al.  A survey of upwind differencing techniques , 1989 .

[25]  S. Osher,et al.  Upwind difference schemes for hyperbolic systems of conservation laws , 1982 .

[26]  P. K. Sweby,et al.  “TVD” Schemes for Inhomogeneous Conservation Laws , 1989 .

[27]  Paul Glaister,et al.  Flux difference splitting for the Euler equations in one spatial co‐ordinate with area variation , 1988 .

[28]  Eli Turkel,et al.  On Numerical Boundary Treatment of Hyperbolic Systems for Finite Difference and Finite Element Methods , 1982 .

[29]  H. C. Yee,et al.  Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics , 1981 .

[30]  S. Spekreijse.  Multigrid solution of monotone second-order discretizations of hyperbolic conservation laws , 1986 .

[31]  H. C. Yee,et al.  Numerical wave propagation in an advection equation with a nonlinear source term , 1992 .

[32]  H. C. Yee,et al.  Boundary Approximations for Implicit Schemes for One-Dimensional Inviscid Equations of Gasdynamics , 1982 .

[33]  Stanley Osher,et al.  Convergence of Generalized MUSCL Schemes , 1985 .

[34]  A. Przekwas,et al.  A comparative study of advanced shock-capturing schemes applied to Burgers' equation , 1990 .

[35]  A. J. Przekwas,et al.  A comparative study of advanced shock-capturing shcemes applied to Burgers' equation , 1992 .

[36]  S. Osher,et al.  One-sided difference approximations for nonlinear conservation laws , 1981 .

[37]  Philip L. Roe,et al.  A comparison of numerical flux formulas for the Euler and Navier-Stokes equations , 1987 .

[38]  B. Koren.  Multigrid and defect correction for the steady Navier-Stokes equations , 1990 .

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