Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters

An approach which closely maintains the non-dissipative nature of classical fourth- or higher-order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient, and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) switch of Harten. Spatially non-dissipative fourth- or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM switch is used to signal the appropriate amount of second- or third-order total variation diminishing (TVD) or essentially non-oscillatory (ENO) types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex?shock interactions, and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock?turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD, positive, or ENO schemes.

[1]  H. C. Yee,et al.  Linearized form of implicit TVD schemes for the multidimensional Euler and Navier-Stokes equations , 1986 .

[2]  P. Olsson.  Summation by parts, projections, and stability. II , 1995 .

[3]  H. C. Yee,et al.  Semi-implicit and fully implicit shock-capturing methods for nonequilibrium flows , 1989 .

[4]  Compressible mixing layer computations with high-order ENO schemes , 1997 .

[5]  Ami Harten,et al.  Implicit TVD schemes for hyperbolic conservation laws in curvilinearcoordinates , 1985 .

[6]  A. Bayliss,et al.  A fourth-order scheme for the unsteady compressible Navier-Stokes equations , 1985 .

[7]  H. C. Yee,et al.  Implicit TVD schemes for hyperbolic conservation laws in curvilinear coordinates , 1987 .

[8]  HI,et al.  The Effects of Numerical Viscosities I . Slowly Moving Shocks , 1996 .

[9]  Rosa Donat,et al.  Studies on error propagation for certain nonlinear approximations to hyperbolic equations: discontinuities in derivatives , 1994 .

[10]  A. Harten,et al.  The artificial compression method for computation of shocks and contact discontinuities. I - Single conservation laws , 1977 .

[11]  H. C. Yee,et al.  A class of high resolution explicit and implicit shock-capturing methods , 1989 .

[12]  Bertil Gustafsson,et al.  Fourth-order difference methods for hyperbolic IBVPs , 1995 .

[13]  J. Lambert.  Computational Methods in Ordinary Differential Equations , 1973 .

[14]  Parviz Moin,et al.  Interaction of isotropic turbulence with shock waves: effect of shock strength , 1997, Journal of Fluid Mechanics.

[15]  Gordon Erlebacher,et al.  High-order ENO schemes applied to two- and three-dimensional compressible flow , 1992 .

[16]  J. Butcher.  The Numerical Analysis of Ordinary Di erential Equa-tions , 1986 .

[17]  J. M. Watt.  Numerical Initial Value Problems in Ordinary Differential Equations , 1972 .

[18]  H. C. Yee,et al.  High-Resolution Shock-Capturing Schemes for Inviscid and Viscous Hypersonic Flows , 1990 .

[19]  H. C. Yee,et al.  A numerical study of a class of TVD schemes for compressible mixing layers , 1989 .

[20]  H. C. Yee,et al.  On symmetric and upwind TVD schemes , 1985 .

[21]  A. Jameson,et al.  Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes , 1981 .

[22]  Robert D. Moser,et al.  A numerical study of turbulent supersonic isothermal-wall channel flow , 1995, Journal of Fluid Mechanics.

[23]  Yanwen Ma,et al.  A High Order Accurate Difference Scheme for Complex Flow Fields , 1997 .

[24]  Brian E. Thompson,et al.  Accuracy evaluation of unsteady CFD numerical schemes by vortex preservation , 1995 .

[25]  H. C. Yee,et al.  Numerical simulation of shock wave diffraction by TVD schemes , 1987 .

[26]  Margot Gerritsen,et al.  Designing an efficient solution strategy for fluid flows. 1. A stable high order finite difference scheme and sharp shock resolution for the Euler equations , 1996 .

[27]  Chi-Wang Shu,et al.  Total variation diminishing Runge-Kutta schemes , 1998, Math. Comput..

[28]  M. Ciment,et al.  Higher order compact implicit schemes for the wave equation , 1975 .

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

[30]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[31]  Ami Harten,et al.  Self adjusting grid methods for one-dimensional hyperbolic conservation laws☆ , 1983 .

[32]  R. Hixon,et al.  High-Accuracy Compact MacCormack-Type Schemes for Computational Aeroacoustics , 1998 .

[33]  Neil D. Sandham,et al.  Direct simulation of turbulence using massively parallel computers , 1997, Parallel CFD.

[34]  H. C. Yee,et al.  Explicit and Implicit Multidimensional Compact High-Resolution Shock-Capturing Methods , 1997 .

[35]  P. Moin,et al.  Turbulence statistics in fully developed channel flow at low Reynolds number , 1987, Journal of Fluid Mechanics.

[36]  C. W. Gear,et al.  Numerical initial value problem~ in ordinary differential eqttations , 1971 .

[37]  Eli Turkel,et al.  Dissipative two-four methods for time-dependent problems , 1976 .

[38]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[39]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

[40]  Björn Sjögreen,et al.  High Order shock capturing methods , 1995 .

[41]  R. F. Warming,et al.  An implicit finite-difference algorithm for hyperbolic systems in conservation-law form. [application to Eulerian gasdynamic equations , 1976 .

[42]  H. C. Yee,et al.  Implicit Total Variation Diminishing (TVD) schemes for steady-state calculations. [in gas dynamics , 1983 .

[43]  Eleuterio F. Toro,et al.  Viscous Flux Limiters , 1992 .

[44]  H. C. Yee,et al.  Semi-implicit and fully implicit shock-capturing methods for hyperbolic conservation laws with stiff source terms , 1987 .

[45]  Mark H. Carpenter,et al.  The accuracy of shock capturing in two spatial dimensions , 1997 .

[46]  R. Hirsh,et al.  Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique , 1975 .

[47]  N. Sandham,et al.  A numerical investigation of the compressible mixing layer , 1989 .

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

[49]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[50]  Shi Jin,et al.  The Effects of Numerical Viscosities , 1996 .

[51]  G. Strang.  On the Construction and Comparison of Difference Schemes , 1968 .

[52]  D. Gottlieb,et al.  The Stability of Numerical Boundary Treatments for Compact High-Order Finite-Difference Schemes , 1993 .