Comparison of coordinate-invariant and coordinate-aligned upwinding for the Euler equations

A floating-shock fitting method for the Euler equations has been developed that uses one-sided spatial differences along and across streamlines. This method has been applied to unsteady shocked flow in a duct with a ramp, to supercritical flow over a circular cylinder, and to transonic flow over airfoils. Compared to methods using coordinate-aligned upwind differencing, the coordinate-invariant upwinding generally required fewer grid points to produce crisp shocks and shears. For transonic airfoils, coordinate-invariant upwind differencing advances the agreement between computations and experiment, and it reduces the grid dependency of the computed results for strong shocks. The computational expenditure is comparable to that of coordinate-aligned upwind methods for the Euler equations.