A high-order compact numerical algorithm for supersonic flows

A dissipative compact two-four scheme (second-order time, fourth-order space) similar to the original MacCormack scheme has been developed, that exhibits greater accuracy than conventional fourth-order schemes. The dissipative nature of the scheme allows it to resolve weak discontinuities without artificial damping. A derivation of the scheme is presented, as well as the theoretical stability characteristics. The temporal scheme is then generalized into a steady-state formulation which achieves fourth-order spatial accuracy at steady-state. Several test problems are used to show that the scheme is more accurate than the traditional MacCormack scheme, and is nearly as efficient.