On the computation of sound in free and wall-bounded domains

A new Dispersion-Relation-Preserving (DRP) scheme has been developed using the Lax–Wendroff methodology. Two collocated grids are placed in a staggered formation and a staggered DRP scheme is used to calculate the spatial differentiation of the propagation and convection terms. A staggered filtering scheme of a six points stencil is developed to complete the transformation from one grid to another. Existing DRP Runge–Kutta schemes are used for the time marching. Stability limits and accuracy issues are investigated using a simple 1D advection equation. The new method is then tested for monopole and quadrupole radiation, diffraction effects of an aperture in a wall, and convection effects of shear flow. All demonstrate the good accuracy and numerical stability of the new method.