Computational Considerations for the Simulation of Shock-Induced Sound

The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved because the chosen algorithm must also resolve discontinuities in the solution. The extent to which a high-order accurate shock-capturing method can be relied upon for aeroacoustics applications that involve the interaction of shocks with other waves has not been previously quantified. Such a study is initiated in this work. A fourth-order accurate essentially nonoscillatory (ENO) method is used to investigate the solutions of inviscid, compressible flows with shocks. The design order of accuracy is achieved in the smooth regions of a steady-state, quasi-one-dimensional test case. However, in an unsteady test case, only first-order results are obtained downstream of a sound-shock interaction. The difficulty in obtaining a globally high-order accurate solution in such a case with a shock-capturing method is demonstrated through the study of a simplified, linear model problem. Some of the difficult issues and ramifications for aeroacoustic simulations of flows with shocks that are raised by these results are discussed.

[1]  S. Osher,et al.  Propagation of error into regions of smoothness for non-linear approximations to hyperbolic equations , 1990 .

[2]  Chi-Wang Shu,et al.  A comparison of two formulations for high-order accurate essentially non-oscillatory schemes , 1994 .

[3]  Chi-Wang Shu.  Numerical solutions of conservation laws , 1986 .

[4]  Eckart Meiburg,et al.  A numerical study of the convergence properties of ENO schemes , 1990 .

[5]  Jay Casper,et al.  Using high-order accurate essentially nonoscillatory schemes for aeroacoustic applications , 1996 .

[6]  Mark H. Carpenter,et al.  Characteristic and Finite-Wave Shock-Fitting Boundary Conditions for Chebyshev Methods , 1994 .

[7]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[8]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[9]  H. L. Atkins,et al.  High-order ENO methods for the unsteady compressible Navier-Stokes equations , 1991 .

[10]  S. Osher,et al.  Regular ArticleUniformly High Order Accurate Essentially Non-oscillatory Schemes, III , 1997 .

[11]  Chi-Wang Shu,et al.  Non-oscillatory spectral Fourier methods for shock wave calculations , 1988 .

[12]  Jay Casper,et al.  Using high-order accurate essentially non-oscillatory schemes for aeroacoustic applications , 1995 .

[13]  Wai Sun Don.  Numerical study of pseudospectral methods in shock wave applications , 1994 .

[14]  S. Lele.  Compact finite difference schemes with spectral-like resolution , 1992 .

[15]  Wei Cai,et al.  Uniform high-order spectral methods for one- and two-dimensional Euler equations , 1991 .

[16]  A. Harten,et al.  Multi-Dimensional ENO Schemes for General Geometries , 1991 .

[17]  Chi-Wang Shu.  Numerical experiments on the accuracy of ENO and modified ENO schemes , 1990 .

[18]  David A. Caughey,et al.  Computing unsteady shock waves for aeroacoustic applications , 1994 .

[19]  A. Harten.  ENO schemes with subcell resolution , 1989 .

[20]  Stanley Osher,et al.  Propagation of error into regions of smoothness for accurate difference approximations to hyperbolic equations , 1977 .

[21]  Peter D. Lax,et al.  The computation of discontinuous solutions of linear hyperbolic equations , 1978 .