Optimized Compact Finite Difference Schemes with Maximum Resolution

Direct numerical simulations and computational aeroacoustics require an accurate finite difference scheme that has a high order of truncation and high-resolution characteristics in the evaluation of spatial derivatives. Compact finite difference schemes are optimized to obtain maximum resolution characteristics in space for various spatial truncation orders. An analytic method with a systematic procedure to achieve maximum resolution characteristics is devised for multidiagonal schemes, based on the idea of the minimization of dispersive (phase) errors in the wave number domain, and these are applied to the analytic optimization of multidiagonal compact schemes. Actual performances of the optimized compact schemes with a variety of truncation orders are compared by means of numerical simulations of simple wave convections, and in this way the most effective compact schemes are found for tridiagonal and pentadiagonal cases, respectively. From these comparisons, the usefulness of an optimized high-order tridiagonal compact scheme that is more efficient than a pentadiagonal scheme is discussed. For the optimized high-order spatial schemes, the feasibility of using classical high-order Runge-Kutta time advancing methods is investigated.

[1]  Z. Kopal,et al.  Numerical analysis , 1955 .

[2]  L. Collatz.  The numerical treatment of differential equations , 1961 .

[3]  B. Celikkol,et al.  Scattering of sound waves by a turbulent wake , 1974 .

[4]  H. Lomax.  Recent Progress in Numerical Techniques for Flow Simulation , 1975 .

[5]  A. Iserles.  Generalized Leapfrog Methods , 1986 .

[6]  Sanjiva K. Lele,et al.  Direct numerical simulation of compressible free shear flows , 1989 .

[7]  R. Clarksean.  Direct numerical simulation of a plane jet using the spectral-compact finite-difference technique , 1991 .

[8]  Scattering of sound waves by a compressible vortex , 1991 .

[9]  Parviz Moin,et al.  Direct computation of the sound from a compressible co-rotating vortex pair , 1992 .

[10]  Christopher K. W. Tam,et al.  Dispersion-relation-preserving schemes for computational aeroacoustics , 1992 .

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

[12]  Philip L. Roe,et al.  Development of non-dissipative numerical schemes for computational aeroacoustics , 1993 .

[13]  Yen Liu,et al.  Analysis of numerical approaches for acoustic equations , 1993 .

[14]  Duck-Joo Lee,et al.  Third Order UNO Schemes with an ACM for Linear Waves and Euler Equations , 1993 .

[15]  C. Tam,et al.  Dispersion-relation-preserving finite difference schemes for computational acoustics , 1993 .

[16]  C. Tam,et al.  A Study of the Short Wave Components in Computational Acoustics , 1993 .

[17]  J. Kim,et al.  Optimized Compact Finite Difference Schemes for Solving Wave Equations , 1996 .

[18]  J. Ekaterinaris.  Upwind Scheme for Acoustic Disturbances Generated by Low-Speed Flows , 1997 .

[19]  Tim Colonius,et al.  Numerically nonreflecting boundary and interface conditions for compressible flow and aeroacoustic computations , 1997 .