A high order compact scheme for hypersonic aerothermodynamics

A novel high order compact scheme for solving the compressible Navier-Stokes equations has been developed. The scheme is an extension of a method originally proposed for solving the Euler equations, and combines several techniques for the solution of compressible flowfields, such as upwinding, limiting and flux vector splitting, with the excellent properties of high order compact schemes. Extending the method to the Navier-Stokes equations is achieved via a Kinetic Flux Vector Splitting technique, which represents an unusual and attractive way to include viscous effects. This approach offers a more accurate and less computationally expensive technique than discretizations based on more conventional operator splitting. The Euler solver has been validated against several inviscid test cases, and results for several viscous test cases are also presented. The results confirm that the method is stable, accurate and has excellent shock-capturing capabilities for both viscous and inviscid flows

An implicit algorithm for the conservative, transonic full-potential equation with effective rotated differencing

A new differencing scheme for the conservative full potential equation which effectively simulates rotated differencing is presented. The scheme was implemented by an appropriate upwind bias of the density coefficient along coordinate directions. A fast, fully implicit, approximate factorization iteration scheme was then used to solve the resulting difference equations. Solutions for a number of traditionally difficult transonic airfoil test cases are presented

On central-difference and upwind schemes

A class of numerical dissipation models for central-difference schemes constructed with second- and fourth-difference terms is considered. The notion of matrix dissipation associated with upwind schemes is used to establish improved shock capturing capability for these models. In addition, conditions are given that guarantee that such dissipation models produce a Total Variation Diminishing (TVD) scheme. Appropriate switches for this type of model to ensure satisfaction of the TVD property are presented. Significant improvements in the accuracy of a central-difference scheme are demonstrated by computing both inviscid and viscous transonic airfoil flows

Three-dimensional multigrid algorithms for the flux-split Euler equations

The Full Approximation Scheme (FAS) multigrid method is applied to several implicit flux-split algorithms for solving the three-dimensional Euler equations in a body fitted coordinate system. Each of the splitting algorithms uses a variation of approximate factorization and is implemented in a finite volume formulation. The algorithms are all vectorizable with little or no scalar computation required. The flux vectors are split into upwind components using both the splittings of Steger-Warming and Van Leer. The stability and smoothing rate of each of the schemes are examined using a Fourier analysis of the complete system of equations. Results are presented for three-dimensional subsonic, transonic, and supersonic flows which demonstrate substantially improved convergence rates with the multigrid algorithm. The influence of using both a V-cycle and a W-cycle on the convergence is examined

Tailoring explicit time-marching schemes to improve convergence characteristics

Multi-stage time-stepping schemes, tailored to chosen spatial-differencing operators, are derived and tested. The schemes are constructed to give optimal damping of the high-frequency waves. They are ideal for use with multi-grid acceleration. The concept of characteristic time-stepping, necessary for the extension of the scalar analysis to systems of equations, is presented. The schemes show a marked improvement over Runge-Kutta schemes

The pros and cons of code validation

Computational and wind tunnel error sources are examined and quantified using specific calculations of experimental data, and a substantial comparison of theoretical and experimental results, or a code validation, is discussed. Wind tunnel error sources considered include wall interference, sting effects, Reynolds number effects, flow quality and transition, and instrumentation such as strain gage balances, electronically scanned pressure systems, hot film gages, hot wire anemometers, and laser velocimeters. Computational error sources include math model equation sets, the solution algorithm, artificial viscosity/dissipation, boundary conditions, the uniqueness of solutions, grid resolution, turbulence modeling, and Reynolds number effects. It is concluded that, although improvements in theory are being made more quickly than in experiments, wind tunnel research has the advantage of the more realistic transition process of a right turbulence model in a free-transition test