Analysis and Improvement of Upwind and Centered Schemes on Quadrilateral and Triangular Meshes

Second-order accurate upwind and centered schemes are presented in a framework that facilitates their analysis and comparison. The upwind scheme employed consists of a reconstruction step (MUSCL approach) followed by an upwind step (Roe's flux-difference splitting). The two centered schemes are of Lax-Friedrichs (L-F) type. They are the nonstaggered versions of the Nessyahu-Tadmor (N-T) scheme and the CE/SE method with epilson = 1/2. The upwind scheme is extended to the case of two spatial dimensions (2D) in a straightforward manner. The N-T and CE/SE schemes are extended in a manner similar to the 2D extensions of the CE/SE scheme by Wang and Chang for a triangular mesh and by Zhang, Yu, and Chang for a quadrilateral mesh. The slope estimates, however, are simplified. Fourier stability and accuracy analyses are carried out for these schemes for the standard 1D and the 2D quadrilateral mesh cases. In the nonstandard case of a triangular mesh, the triangles must be paired up when analyzing the upwind and N-T schemes. An observation resulting in an extended N-T scheme which is faster and uses only one-third of the storage for flow data compared with the CE/SE method is presented. Numerical results are shown. Other improvements to the schemes are discussed.

[1]  R. Courant,et al.  On the solution of nonlinear hyperbolic differential equations by finite differences , 1952 .

[2]  High accuracy capture of curved shock fronts using the method of space-time conservation element and solution elemen , 1998 .

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

[4]  High-resolution Non-oscillatory Central Schemes with Non-staggered Grids for Hyperbolic Conservation Laws Dedicated to Our Friend and Colleague , 1997 .

[5]  E. Tadmor,et al.  Non-oscillatory central differencing for hyperbolic conservation laws , 1990 .

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

[7]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[8]  Timothy J. Barth,et al.  The design and application of upwind schemes on unstructured meshes , 1989 .

[9]  Paul Arminjon,et al.  A Finite Volume Extension of the Lax-Friedrichs and Nessyahu-Tadmor Schemes for Conservation Laws on Unstructured Grids , 1998 .

[10]  Sin-Chung Chang,et al.  A 3-D non-splitting structured/unstructured Euler solver based on the space-time conservation element and solution element method , 1999 .

[11]  H. Huynh.  Accurate upwind schemes for the Euler equations , 1995 .

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

[13]  S. Osher.  Riemann Solvers, the Entropy Condition and High Resolution Difference Approximations, , 1984 .

[14]  S. Osher,et al.  High-Resolution Nonoscillatory Central Schemes with Nonstaggered Grids for Hyperbolic Conservation Laws , 1998 .

[15]  Meng-Sing Liou,et al.  An Accurate and Robust Flux Splitting Scheme for Shock and Contact Discontinuities , 1997, SIAM J. Sci. Comput..

[16]  Sin-Chung Chang,et al.  The Space-Time Conservation Element and Solution Element Method: A New High-Resolution and Genuinely Multidimensional Paradigm for Solving Conservation Laws. 1; The Two Dimensional Time Marching Schemes , 1998 .

[17]  Sin-Chung Chang.  The Method of Space-Time Conservation Element and Solution Element-A New Approach for Solving the Navier-Stokes and Euler Equations , 1995 .

[18]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

[19]  H. Huynh.  Accurate upwind methods for the Euler equations , 1995 .

[20]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[21]  Eitan Tadmor,et al.  Nonoscillatory Central Schemes for Multidimensional Hyperbolic Conservation Laws , 1998, SIAM J. Sci. Comput..

[22]  Sin-Chung Chang,et al.  A space-time conservation element and solution element method for solving the two- and three-dimensional unsteady euler equations using quadrilateral and hexahedral meshes , 2002 .

[23]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[24]  H. C. Yee,et al.  Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations. [in gas dynamics , 1983 .

[25]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .