An algebraic turbulence model for two- and three-dimensional separated flows is specified that avoids the necessity for finding the edge of the boundary layer. Properties of the model are determined and comparisons made with experiment for an incident shock on a flat plate, separated flow over a compression corner, and transonic flow over an airfoil. Separation and reattachment points from numerical Navier-Stokes solutions agree with experiment within one boundary-layer thickness. Use of law-of-the-wall boundary conditions does not alter the predictions significantly. Applications of the model to other cases are contained in companion papers.

Thin-layer approximation and algebraic model for separated turbulent flows

Flux vector splitting of the inviscid gasdynamic equations with application to finite-difference methods

Various artificial dissipation models that are used with central difference algorithms for the Euler equations are analyzed for their effect on accuracy, stability, and convergence rates. In particular, linear and nonlinear models are investigated using an implicit approximate factorization code (ARC2D) for transonic airfoils. Fully implicit application of the dissipation models is shown to improve robustness and convergence rates. The treatment of dissipation models at boundaries will be examined. It will be shown that accurate, error free solutions with sharp shocks can be obtained using a central difference algorithm coupled with an appropriate nonlinear artificial dissipation model. I. Introduction T HE solution of the Euler equations using numerical techniques requires the use of either a differencing method with inherent dissipation or the addition of dissipation terms to a nondissipative scheme. This is because the Euler equations do not provide any natural dissipation mechanism (such as viscosity in the Navier-Stokes equations) that would eliminate high frequencies which are caused by nonlinearitie s and especially shocks. A variety of numerical algorithms and computer codes for the Euler equations have been developed. Methods such as MacCormack's1 explicit

Artificial Dissipation Models for the Euler Equations

In 1977 and 1978, general purpose centrally space differenced implicit finite difference codes in two and three dimensions have been introduced. These codes, now called ARC2D and ARC3D, can run either in inviscid or viscous mode for steady or unsteady flow. Since the introduction of the ARC2D and ARC3D codes, overall computational efficiency could be improved by making use of a number of algorithmic changes. These changes are related to the use of a spatially varying time step, the use of a sequence of mesh refinements to establish approximate solutions, implementation of various ways to reduce inversion work, improved numerical dissipation terms, and more implicit treatment of terms. The present investigation has the objective to describe the considered improvements and to quantify advantages and disadvantages. It is found that using established and simple procedures, a computer code can be maintained which is competitive with specialized codes.

Recent improvements in efficiency, accuracy, and convergence for implicit approximate factorization algorithms

Two complementary procedures were developed to calculate the viscous supersonic flow over conical shapes at large angles of attack, with application to cones and delta wings. In the first approach the flow is assumed to be conical and the governing equations are solved at a given Reynolds number with a time-marching explicit finite-difference algorithm. In the second method the parabolized Navier-Stokes equations are solved with a space-marching implicit noniterative finite-difference algorithm. This latter approach is not restricted to conical shapes and provides a large improvement in computational efficiency over published methods. Results from the two procedures agree very well with each other and with available experimental data.

/pdf/calculation-of-supersonic-viscous-flow-over-delta-wings-with-4j1q0mbvps.pdf

Calculation of supersonic viscous flow over delta wings with sharp subsonic leading edges

A numerical method for solving the compressible form of the unsteady Navier-Stokes equations is described. This method was originally presented in 1970 and has since been modified during the development of computer programs at Ames for implementing models that account for the effects of turbulence in shock-induced separated flows. Although this paper does not describe the turbulence models themselves, a complete description of the basic numerical method is given with emphasis on the choice of a computational mesh for high Reynolds number flows, finite-difference approximations for mixed partial derivatives, extension of the Courant-Friedrichs-Lewy stability condition for viscous flows, mesh boundary conditions, and numerical smoothing for strong shock-wave calculations.

B. S. Baldwin

Papers

A numerical method for solving the Navier-Stokes equations with application to shock-boundary layer interactions