or perhaps with a forcing term g(u, x, t) on the right-hand side. Here u ϭ (u 1 , ..., u m ), f ϭ (f 1 , ..., f d ), x ϭ (x 1 , ..., x d )In this paper, we further analyze, test, modify, and improve the high order WENO (weighted essentially non-oscillatory) finite differand t Ͼ 0. involves a procedure of reconstructing point values from Euler systems of gas dynamics, we suggest computing the weights cell averages and could become complicated and costly for from pressure and entropy instead of the characteristic values to multi-dimensional problems. Later, Shu and Osher [14,15] simplify the costly characteristic procedure. The resulting WENO developed the flux version of ENO schemes which do not schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights and work well require such a reconstruction procedure. We will formulate for problems which do not contain strong shocks or strong reflected the WENO schemes based on this flux version of ENO waves. We also prove that, for conservation laws with smooth soluschemes. The WENO schemes of Liu et al. [9] are based tions, all WENO schemes are convergent. Many numerical tests, on the cell-averaged version of ENO schemes.including the 1D steady state nozzle flow problem and 2D shock For applications involving shocks, second-order schemes entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the are usually adequate if only relatively simple structures WENO scheme using the new smoothness measurement in resolv-are present in the smooth part of the solution (e.g., the ing complicated shock and flow structures. We have also applied shock tube problem). However, if a problem contains rich Yang's artificial compression method to the WENO schemes to structures as well as shocks (e.g., the shock entropy wave sharpen contact discontinuities. ᮊ 1996 Academic Press, Inc.interaction problem in Example 4, Section 8.3), high order shock capturing schemes (order of at least three) are more

Abstract. In this paper we further explore a class of high order TVD (total variation diminishing) Runge-Kutta time discretization initialized in a paper by Shu and Osher, suitable for solving hyperbolic conservation laws with stable spatial discretizations. We illustrate with numerical examples that non-TVD but linearly stable Runge-Kutta time discretization can generate oscillations even for TVD (total variation diminishing) spatial discretization, verifying the claim that TVD Runge-Kutta methods are important for such applications. We then explore the issue of optimal TVD Runge-Kutta methods for second, third and fourth order, and for low storage Runge-Kutta methods.

In this paper we design a class of numerical schemes that are higher-order extensions of the weighted essentially non-oscillatory (WENO) schemes of G.-S. Jiang and C.-W. Shu (1996) and X.-D. Liu, S. Osher, and T. Chan (1994). Used by themselves, the schemes may not always be monotonicity preserving but coupled with the monotonicity preserving bounds of A. Suresh and H. T. Huynh (1997) they perform very well. The resulting monotonicity preserving weighted essentially non-oscillatory (MPWENO) schemes have high phase accuracy and high order of accuracy. The higher-order members of this family are almost spectrally accurate for smooth problems. Nevertheless, they, have robust shock capturing ability. The schemes are stable under normal CFL numbers. They are also efficient and do not have a computational complexity that is substantially greater than that of the lower-order members of this same family of schemes. The higher accuracy that these schemes offer coupled with their relatively low computational complexity makes them viable competitors to lower-order schemes, such as the older total variation diminishing schemes, for problems containing both discontinuities and rich smooth region structure. We describe the MPWENO schemes here as well as show their ability to reach their designed accuracies for smooth flow. We also examine the role of steepening algorithms such as the artificial compression method in the design of very high order schemes. Several test problems in one and two dimensions are presented. For multidimensional problems where the flow is not aligned with any of the grid directions it is shown that the present schemes have a substantial advantage over lower-order schemes. It is argued that the methods designed here have great utility for direct numerical simulations and large eddy simulations of compressible turbulence. The methodology developed here is applicable to other hyperbolic systems, which is demonstrated by showing 1 that the MPWENO schemes also work very well on magnetohydrodynamical test problems.

No abstract

No abstract

No abstract

No abstract

No abstract