Publications

The formulation, implementation and usage of a numerical solution verification code is described. This code uses the Richardson extrapolation procedure to estimate the order of accuracy and error of a computational program solution. It evaluates multiple solutions performed in numerical grid convergence studies to verify a numerical algorithm implementation. Analyses are performed on both structured and unstructured grid codes. Finite volume and finite element discretization programs are examined. Two and three-dimensional solutions are evaluated. Steady state and transient solution analysis capabilities are present in the verification code. Multiple input data bases are accepted. Benchmark options are included to allow for minimal solution validation capability as well as verification.

The Detached Eddy Simulation (DES) and steadystate Reynolds-Averaged Navier-Stokes (RANS) turbulence modeling approaches are examined for the incompressible flow over a square cross-section cylinder at a Reynolds number of 21,400. A compressible flow code is used which employes a second-order Roe upwind spatial discretization. Efforts are made to assess the numerical accuracy of the DES predictions with regards to statistical convergence, iterative convergence, and temporal and spatial discretization error. Three-dimensional DES simulations compared well with two-dimensional DES simulations, suggesting that the dominant vortex shedding mechanism is effectively two-dimensional. The two-dimensional simulations are validated via comparison to experimental data for mean and RMS velocities as well as Reynolds stress in the cylinder wake. The steady-state RANS models significantly overpredict the size of the recirculation zone, thus underpredicting the drag coefficient relative to the experimental value. The DES model is found to give good agreement with the experimental velocity data in the wake, drag coefficient, and recirculation zone length.

Hypersonic transitional flows over a flat plate and a sharp cone are studied using four turbulence models: the one-equation eddy viscosity transport model of Spalart-Allmaras, a low Reynolds number κ-ε model, the Menter κ-ω model, and the Wilcox κ-ω model. A framework is presented for the assessment of turbulence models that includes documentation procedures, solution accuracy, model sensitivity, and model validation. The accuracy of the simulations is addressed, and the sensitivities of the models to grid refinement, freestream turbulence levels, and wall y+ spacing are presented. The flat plate skin friction results are compared to the well-established laminar and turbulent correlations of Van Driest. Correlations for the sharp cone are discussed in detail. These correlations, along with recent experimental data, are used to judge the validity of the simulation results for skin friction and surface heating on the sharp cone. The Spalart-Allmaras performs the best with regards to model sensitivity and model accuracy, while the Menter κ-ω model also performs well for these zero pressure gradient boundary layer flows. © 2001 American Institute of Aeronautics & Astronautics.

This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity - transport equation for non-turbulent fluctuation growth based on that proposed by Warren and Hassan (Journal of Aircraft, Vol. 35, No. 5) is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittence function based on the work of Dhawan and Narasimha (Journal of Fluid Mechanics, Vol. 3, No. 4). The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the spatial accuracy of selected predictions is analyzed.

This paper reports on the development of a unified one-equation model for the prediction of transitional and turbulent flows. An eddy viscosity--transport equation for nonturbulent fluctuation growth based on that proposed by Warren and Hassan is combined with the Spalart-Allmaras one-equation model for turbulent fluctuation growth. Blending of the two equations is accomplished through a multidimensional intermittency function based on the work of Dhawan and Narasimha. The model predicts both the onset and extent of transition. Low-speed test cases include transitional flow over a flat plate, a single element airfoil, and a multi-element airfoil in landing configuration. High-speed test cases include transitional Mach 3.5 flow over a 5{degree} cone and Mach 6 flow over a flared-cone configuration. Results are compared with experimental data, and the grid-dependence of selected predictions is analyzed.

Many Navier-Stokes codes require that the governing equations be written in conservation form with a source term. The Spalart-Allmaras one-equation model was originally developed in substantial derivative form and when rewritten in conservation form, a density gradient term appears in the source term. This density gradient term causes numerical problems and has a small influence on the numerical predictions. Further work has been performed to understand and to justify the neglect of this term. The transition trip term has been included in the one-equation eddy viscosity model of Spalart-Allmaras. Several problems with this model have been discovered when applied to high-speed flows. For the Mach 8 flat plate boundary layer flow with the standard transition method, the Baldwin-Barth and both k-{omega} models gave transition at the specified location. The Spalart-Allmaras and low Reynolds number k-{var_epsilon} models required an increase in the freestream turbulence levels in order to give transition at the desired location. All models predicted the correct skin friction levels in both the laminar and turbulent flow regions. For Mach 8 flat plate case, the transition location could not be controlled with the trip terms as given in the Spalart-Allmaras model. Several other approaches have been investigated to allow the specification of the transition location. The approach that appears most appropriate is to vary the coefficient that multiplies the turbulent production term in the governing partial differential equation for the eddy viscosity (Method 2). When this coefficient is zero, the flow remains laminar. The coefficient is increased to its normal value over a specified distance to crudely model the transition region and obtain fully turbulent flow. While this approach provides a reasonable interim solution, a separate effort should be initiated to address the proper transition procedure associated with the turbulent production term. Also, the transition process might be better modeled with the Spalart-Allmaras turbulence model with modification of the damping function f{sub v1}. The damping function could be set to zero in the laminar flow region and then turned on through the transition flow region.

A number of one- and two-equation turbulence models are examined for hypersonic perfect- and real-gas flows with laminar, transitional, and turbulent flow regions. These models were generally developed for incompressible flows, and the extension to the hypersonic flow regime is discussed. In particular, inconsistencies in the formulation of diffusion terms for one-equation models are examined. For the Spalart-Allmaras model, the standard method for forcing transition at a specified location is found to be inadequate for hypersonic flows. An alternative transition method is proposed and evaluated for a Mach 8 flat plate test case. This test case is also used to evaluate three different two-equation turbulence models:.a low Reynolds number k - ε model, the Menter k-ω formulation, and the Wilcox (1998) k -ω model. These one- and two-equation models are then applied to the Mach 20 Reentry F flight vehicle. The Spalart-Allmaras model and both k-ω formulations are found to provide good agreement with the flight data for heat flux, while the Baldwin-Barth and low Reynolds number k - ε models overpredict the turbulent heating rates. Careful attention is given to solution verification in the areas of both iterative and grid convergence. © protection in the United States.

The unsteady thin-layer Navier-Stokes equations for a perfect gas are solved with a linearized block Alternating Direction Implicit finite-difference solution procedure. Solution errors due to numerical dissipation added to the governing equations are evaluated. Errors in the numerical predictions on three different grids are determined where Richardson extrapolation is used to estimate the exact solution. Accurate computational results are tabulated for the hypersonic laminar flow over a spherical body which can be used as a benchmark test case. Predictions obtained from the code are in good agreement with inviscid numerical results arid experimental data.