A Hermite upwind WENO scheme for solving hyperbolic conservation laws

This paper proposes a new WENO procedure to compute problems containing both discontinuities and a large disparity of characteristic scales. In a one-dimensional context, the WENO procedure is defined on a three-points stencil and designed to be sixth-order in regions of smoothness. We define a finite-volume discretization in which we consider the cell averages of the variable and its first derivative as discrete unknowns. The reconstruction of their point-values is then ensured by a unique sixth-order Hermite polynomial. This polynomial is considered as a symmetric and convex combination, by ideal weights, of three fourth-order polynomials: a central polynomial, defined on the three-points stencil, is combined with two polynomials based on the left and the right two-points stencils. The symmetric nature of such an interpolation has an important consequence: the choice of ideal weights has no influence on the properties of the discretization. This advantage enables to formulate the Hermite interpolation for non-uniform meshes. Following the methodology of the classic WENO procedure, nonlinear weights are then defined. To deal with the peculiarities of the Hermite interpolation near discontinuities, we define a new procedure in order for the nonlinear weights to smoothly evolve between the ideal weights, in regions of smoothness, and one-sided weights, otherwise. The resulting scheme is a sixth-order WENO method based on central Hermite interpolation and TVD Runge-Kutta time-integration. We call this scheme the HCWENO6 scheme. Numerical experiments in the scalar and the 1D Euler cases make it possible to check and to validate the options selected. In these experiments, we emphasize the resolution power of the method by computing test cases that model realistic aero-acoustic problems.

[1]  R. Hindman On shock capturing methods and why they work , 1988 .

[2]  J. F. Mckenzie,et al.  Interaction of Linear Waves with Oblique Shock Waves , 1968 .

[3]  Gabriella Puppo,et al.  Compact Central WENO Schemes for Multidimensional Conservation Laws , 1999, SIAM J. Sci. Comput..

[4]  D. M. Bushnell,et al.  Numerical computations of turbulence amplification in shock wave interactions , 1984 .

[5]  Chi-Wang Shu Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[6]  T. Yabe,et al.  Exactly conservative semi-Lagrangian scheme for multi-dimensional hyperbolic equations with directional splitting technique , 2001 .

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

[8]  Sergio Pirozzoli,et al.  Conservative Hybrid Compact-WENO Schemes for Shock-Turbulence Interaction , 2002 .

[9]  Milo D. Dahl,et al.  Third Computational Aeroacoustics (CAA) Workshop on Benchmark Problems , 2000 .

[10]  S. Lele Compact finite difference schemes with spectral-like resolution , 1992 .

[11]  François Bouchut,et al.  A MUSCL method satisfying all the numerical entropy inequalities , 1996, Math. Comput..

[12]  Charles Hirsch,et al.  Numerical computation of internal & external flows: fundamentals of numerical discretization , 1988 .

[13]  Jianxian Qiu,et al.  On the construction, comparison, and local characteristic decomposition for high-Order central WENO schemes , 2002 .

[14]  Yuxin Ren,et al.  A characteristic-wise hybrid compact-WENO scheme for solving hyperbolic conservation laws , 2003 .

[15]  Mark H. Carpenter,et al.  Computational Considerations for the Simulation of Shock-Induced Sound , 1998, SIAM J. Sci. Comput..

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

[17]  A. Edelman,et al.  Nonnegativity-, monotonicity-, or convexity-preserving cubic and quintic Hermite interpolation , 1989 .

[18]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

[19]  P. Lax,et al.  On Upstream Differencing and Godunov-Type Schemes for Hyperbolic Conservation Laws , 1983 .

[20]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

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

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

[23]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method II: Two dimensional case , 2005 .

[24]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method: one-dimensional case , 2004 .

[25]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[26]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .