Moment preserving schemes for Euler equations

A high order accurate finite difference scheme is proposed for one-dimensional Euler equations. In the scheme a set of first three moments of each signal are preserved during the updating. The scheme is one of 5th order in space and 4th order in time. This feature is different from that in typical existing methods in which the use of the first three polynomials results in only 3rd order accuracy in space. The scheme has different features from the existing high order schemes, and the most noticeable are the simultaneous discretization both in space and time, and the use of moments of Riemann invariants instead of primitive physical variables. Numerical examples are given to show the accuracy of the scheme and its robustness for the flows involving shocks.