Comparative advantages of high-order schemes for subsonic, transonic, and supersonic flows

This computational study aims to verify the order of accuracy of the COMPACT and weighted essentially non-oscillatory (WENO) finite difference schemes implemented in the AEROFLO software, and to identify the comparative advantages of these schemes relative to the low-order (MUSCL-based) schemes for a range of flow problems. The method of manufactured solutions was used to determine the order of accuracy of the spatial differencing schemes. The theoretical sixth-order of accuracy is verified for the COMPACT scheme for subsonic flows, while the theoretical fith-order WENO scheme exhibited a 3.5 order of accuracy for supersonic flows. The MUSCL scheme shows the theoretical secondorder accuracy for all flow regimes. The accuracy results were observed for both Cartesian and curvilinear grids. Several subsonic, transonic, and supersonic calculations were then used to evaluate the results from the highand low-order schemes. For the subsonic and transonic flow configurations, the high-order schemes generally require smaller CPU times, due to their ability to use larger time step sizes or their ability to generate better results with coarser grids as compared to the low-order schemes. For the supersonic flow configurations, both the highand the low-order schemes capture the shock locations very accurately, although the low-order schemes tend to exhibit significantly larger numerical noise in the regions behind the shocks.

[1]  Ruey-Hor Yen,et al.  Implicit Weighted Essentially Nonoscillatory Schemes for the Compressible Navier-Stokes Equations , 2001 .

[2]  Miguel R. Visbal,et al.  High-Order Schemes for Navier-Stokes Equations: Algorithm and Implementation Into FDL3DI , 1998 .

[3]  F. S. Billig,et al.  Shock-wave shapes around spherical-and cylindrical-nosed bodies. , 1967 .

[4]  T. Pulliam,et al.  A diagonal form of an implicit approximate-factorization algorithm , 1981 .

[5]  A. Roshko.  On the development of turbulent wakes from vortex streets , 1953 .

[6]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[7]  Patrick J. Roache,et al.  Symbolic manipulation and computational fluid dynamics , 1983 .

[8]  Christopher J. Roy,et al.  Verification of Euler/Navier–Stokes codes using the method of manufactured solutions , 2004 .

[9]  Viktoria Schmitt,et al.  Pressure distributions on the ONERA M6 wing at transonic Mach numbers , 1979 .

[10]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[11]  D. Fyfe,et al.  Economical evaluation of Runge-Kutta formulae , 1966 .

[12]  Miguel R. Visbal,et al.  On the use of higher-order finite-difference schemes on curvilinear and deforming meshes , 2002 .

[13]  Patrick Knupp,et al.  Code Verification by the Method of Manufactured Solutions , 2000 .

[14]  C. D. Harris.  Two-dimensional aerodynamic characteristics of the NACA 0012 airfoil in the Langley 8 foot transonic pressure tunnel , 1981 .

[15]  R. F. Warming,et al.  An Implicit Factored Scheme for the Compressible Navier-Stokes Equations , 1977 .

[16]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[17]  P. Spalart.  A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

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

[19]  P. Roe.  CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

[20]  D. Wilcox.  Turbulence modeling for CFD , 1993 .

[21]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[22]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .