It is shown by using a Dispersion-Relation-Preserving [Formula: see text] finite difference scheme that it is feasible to perform direct numerical simulation of acoustic wave propagation problems. The finite difference equations of the [Formula: see text] scheme have essentially the same Fourier-Laplace transforms and hence dispersion relations as the original linearized Euler equations over a broad range of wavenumbers (here referred to as long waves). Thus it is guaranteed that the acoustic waves, the entropy and the vorticity waves computed by the [Formula: see text] scheme are good approximations of those of the exact solutions of Euler equations as long as the wavenumbers are in the long wave range. Computed waves with higher wavenumber, or the short waves, generally have totally different propagation characteristics. There are no counterparts of such waves in the exact solutions. The short waves of a computation scheme are, therefore, contaminants of the numerical solutions. The characteristics of these short waves are analyzed here by group velocity consideration and standard dispersive wave theory. Numerical results of direct simulations of these waves are reported. These waves can be generated by discontinuous initial conditions. To purge the short waves so as to improve the quality of the numerical solution, it is suggested that artificial selective damping terms be added to the finite difference scheme. It is shown how the coefficients of such damping terms may be chosen so that damping is confined primarily to the high wavenumber range. This is important for then only the short waves are damped leaving the long waves basically unaffected. The effectiveness of the artificial selective damping terms is demonstrated by direct numerical simulations involving acoustic wave pulses with discontinuous wave fronts.
On the properties of numerical solutions of dynamical systems obtained using the midpoint method
