High-resolution difference methods with exact evolution for multidimensional waves

We consider the generalization of high-order upwind Strang methods for simulating waves. In $1+1$ dimensions the methods can be defined via the exact evolution over a single time step of an odd-order piecewise polynomial interpolant of the grid data. We construct a true multidimensional version for acoustic waves by applying the solution operator in integral form to the interpolant. We also examine the replacement of polynomials by bandlimited interpolation functions (BLIFs). Numerical experiments with turbulent wave fields are presented to verify the accuracy and stability of the multidimensional methods and to assess the relative effectiveness of the two interpolation techniques.