SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
A recently developed multiple-scale turbulence model which attempts to circumvent the deficiencies of earlier models by taking nonequilibrium spectral energy transfer into account is presented. The model's validity is tested by predicting the confined swirling coaxial jet flow in a sudden expansion. It is noted that, in order to account for anisotropic turbulence, a full Reynolds stress model is required.
This study focuses on the prediction of airfoil characteristics, including lift and drag over a range of Reynolds numbers. Two different turbulence models, which represent two different types of models, are tested. The first is a standard isotropic eddy-viscosity two-equation model, and the second is an explicit algebraic stress model (EASM). The turbulent flow field over a general-aviation airfoil (GA(W)-2) at three Reynolds numbers is studied. At each Reynolds number, predicted lift and drag values at different angles of attack are compared with experimental results, and predicted variations of stall locations with Reynolds number are compared with experimental data. Finally, the size of the separation zone predicted by each model is analyzed, and correlated with the behavior of the lift coefficient near stall. In summary, the EASM model is able to predict the lift and drag coefficients over a wider range of angles of attack than the two-equation model for the three Reynolds numbers studied. However, both models are unable to predict the correct lift and drag behavior near the stall angle, and for the lowest Reynolds number case, the two-equation model did not predict separation on the airfoil near stall.
Supersonic jet plumes were studied using a two-equation turbulence model employing corrections for compressible dissipation and pressure-dilatation. A space-marching procedure based on an upwind numerical scheme was used to solve the governing equations and turbulence transport equations. The computed results indicate that two-equation models employing corrections for compressible dissipation and pressure-dilatation yield improved agreement with the experimental data. In addition, the numerical study demonstrates that the computed results are sensitive to the effect of grid refinement and insensitive to the type of velocity profiles used at the inflow boundary for the cases considered in the present study.
The Johnson-King turbulence model as extended to three-dimensional flows was evaluated using finite-difference boundary-layer direct method. Calculations were compared against the experimental data of the well-known Berg-Elsenaar incompressible flow over an infinite swept-wing. The Johnson-King model, which includes the nonequilibrium effects in a developing turbulent boundary-layer, was found to significantly improve the predictive quality of a direct boundary-layer method. The improvement was especially visible in the computations with increased three-dimensionality of the mean flow, larger integral parameters, and decreasing eddy-viscosity and shear stress magnitudes in the streamwise direction; all in better agreement with the experiment than simple mixing-length methods.
The nonequilibrium algebraic turbulence model of Johnson and King has been applied successfully to both airfoils and wings. However, in wing applications a different expression for the outer 'equilibrium' eddy viscosity has been used because of apparent difficulties in determining the boundary-layer thickness. It is shown in this paper, that this expression does not give eddy viscosity levels for two-dimensional flows that are consistent with Clauser's correlation for equilibrium turbulent boundary layers. And as a consequence, results in inferior predictions when used in airfoil calculations. The objective of this paper is to develop a formulation applicable to wing flows which reduces to the proper form in the two-dimensional case. This necessitates, in part, finding a better way of estimating the boundary-layer thickness from the computed flowfield.
The flow in the base region of launch vehicles is characterized by flow separation, flow reversals, and reattachment. Computation of the convective heat flux in the base region and on the nozzle external surface of Space Shuttle Main Engine and Space Transportation Main Engine (STME) is an important part of defining base region thermal environments. Several turbulence models were incorporated in a CFD code and validated for flow and heat transfer computations in the separated and reattaching regions associated with subsonic and supersonic flows over backward facing steps. Heat flux computations in the base region of a single STME engine and a single S1C engine were performed using three different wall functions as well as a renormalization-group based k-epsilon model. With the very limited data available, the computed values are seen to be of the right order of magnitude. Based on the validation comparisons, it is concluded that all the turbulence models studied have predicted the reattachment location and the velocity profiles at various axial stations downstream of the step very well.
Navier-Stokes solutions about a supercritical airfoil with aileron deflection have been computed using the CFL3D code coupled with the Baldwin-Lomax, Johnson-King, Baldwin-Barth, and Spalart-Allmaras turbulence models. Computations were made at a Mach number of 0.716 and chord Reynolds numbers of 5, 15, and 25 million. The airfoil was analyzed with both 0 deg and 2 deg (TED) aileron deflections. Comparisons over a range of angles-of-attack showed that solutions obtained using the Baldwin-Barth turbulence model presented the best agreement with experimental pressures and sectional lift coefficients. However, Reynolds number trends in sectional lift coefficient and in aileron effectiveness were not predicted consistently.
A two-dimensional-dependent calculation of a forced turbulent mixing layer is reported. The approach employs a k-epsilon turbulence model to represent the effects of motions which are three dimensional or out of phase with the forcing cycle. Calculations that correspond to conditions of the measurements of Weisbrot (1984) and Weisbrot and Wygnanski (1988) indicate that the main features of the mixing layer development are captured by the computation. The effect of grid refinement on the computed result indicates the importance of adequate resolution of the near field in correctly reproducing the flow evolution. The inlet boundary condition, similarly, is shown to play a decisive role in determining subsequent development.
A k-epsilon turbulence model suitable for compressible flow, including the new wall function formulation, has been incorporated into an existing compressible Reynolds-averaged Navier-Stokes code, F3D. The low Reynolds number k-epsilon model of Chien (1982) was added for comparison with the present method. A number of features were added to the F3D code including improved far-field boundary conditions and viscous terms in the streamwise direction. A series of computations of increasing complexity was run to test the effectiveness of the new formulation. Flow over a flat plate was computed by using both orthogonal and nonorthogonal grids, and the friction coefficients and velocity profiles compared with a semi-empirical equation. Flow over a body of revolution at zero angle of attack was then computed to test the method's ability to handle flow over a curved surface. Friction coefficients and velocity profiles were compared to test data. All models gave good results on a relatively fine grid, but only the wall function formulation was effective with coarser grids. Finally, in order to demonstrate the method's ability to handle complex flow fields, separated flow over a prolate spheroid at angle of attack was computed, and results were compared to test data. The results were also compared to a k-epsilon model by Kim and Patel (1991), in which one equation model patched in at the wall was employed. Both models gave reasonable solutions, but improvement is required for accurate prediction of friction coefficients in the separated regions.
Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.
Computational uncertainties in turbulence modeling for three dimensional inlet flow fields include flows approaching separation, strength of secondary flow field, three dimensional flow predictions of vortex liftoff, and influence of vortex-boundary layer interactions; computational uncertainties in vortex generator modeling include representation of generator vorticity field and the relationship between generator and vorticity field. The objectives of the inlet flow field studies presented in this document are to advance the understanding, prediction, and control of intake distortion and to study the basic interactions that influence this design problem.
Four different turbulence models are used to compute the flow over a three-element airfoil configuration. These models are the one-equation Baldwin-Barth model, the one-equation Spalart-Allmaras model, a two-equation k-omega model, and a new one-equation Durbin-Mansour model. The flow is computed using the INS2D two-dimensional incompressible Navier-Stokes solver. An overset Chimera grid approach is utilized. Grid resolution tests are presented, and manual solution-adaptation of the grid was performed. The performance of each of the models is evaluated for test cases involving different angles-of-attack, Reynolds numbers, and flap riggings. The resulting surface pressure coefficients, skin friction, velocity profiles, and lift, drag, and moment coefficients are compared with experimental data. The models produce very similar results in most cases. Excellent agreement between computational and experimental surface pressures was observed, but only moderately good agreement was seen in the velocity profile data. In general, the difference between the predictions of the different models was less than the difference between the computational and experimental data.
Four turbulence models are described and evaluated for transonic flows using the upwind code CFL3D and the central-difference code TLNS3D. In particular, the effects of recent modifications to the half-equation model of Johnson-King are explored in detail, and different versions of the model are compared. This model can obtain good results for both two-dimensional (2D) and three-dimensional (3D) separated flows. The one-equation models of Baldwin-Barth and Spalart-Allmaras perform well for separated airfoil flows, but can predict the shock too far forward at the outboard stations of a separated wing. The equilibrium model of Baldwin-Lomax predicts the shock location too far aft for both 2D and 3D separated flows, as expected. In general, all models perform well for attached or mildly separated flows.
Three-dimensional Navier-Stokes simulations have been done for transonic and low supersonic flow past a nonaxisymmetric nozzle typical of those advocated for advanced fighter airplanes. The jet exhaust is included in the calculations. The investigation compares the performance of the unmodified Baldwin-Lomax turbulence model with its performance when enhanced by the Degani-Schiff and the Goldberg modifications. Solutions are presented for Mach numbers of 0.80, 0.94, and 1.20 at 0-deg angle of attack and a Reynolds number of 20 x 10 to the 6th. The numerical results, which are compared to the wind-tunnel data, show that the three turbulence models predict considerably different shock locations, separated-flow regions, and flowfields.
We formulate a new second-order closure turbulence model by employing a recent closure for the pressure-temperature correlation at the equation level. As a result, we obtain new heat flux equations that avoid the long-standing issue of a finite critical Richardson number. The new, structurally simpler model improves on the Mellor-Yamada 1982 and Galperin et al. 1988 models; key feature includes enhanced mixing under stable conditions facilitating agreement with observational, experimental and high-resolution numerical data sets. The model predicts a planetary boundary layer height deeper than predicted by models with low critical Richardson numbers, as demonstrated in single column model runs of the GISS ModelE general circulation model.
A computer code for the evaluation and/or optimization of the predicative potential of second order turbulent closure models in simple two dimensional flow configurations is discussed. A procedure for the numerical solution of the steady constant property Navier-Stokes equations are described together with algebraic, one dimensional and two dimensional equations of turbulence closure models. Four turbulence models are compared with several sets of experimental data. The effects of initial conditions and boundary conditions are also described. The effects of purely numerical parameters, such as mesh size, boundary locations, and convergence criteria are presented.
The importance of reducing jet noise in both commercial and military aircraft applications has made jet acoustics a significant area of research. A technique for jet noise prediction commonly employed in practice is the MGB approach, based on the Lighthill acoustic analogy. This technique requires as aerodynamic input mean flow quantities and turbulence quantities like the kinetic energy and the dissipation. The purpose of the present paper is to assess existing capabilities for predicting these aerodynamic inputs. Two modern Navier-Stokes flow solvers, coupled with several modern turbulence models, are evaluated by comparison with experiment for their ability to predict mean flow properties in a supersonic jet plume. Potential weaknesses are identified for further investigation. Another comparison with similar intent is discussed by Barber et al. The ultimate goal of this research is to develop a reliable flow solver applicable to the low-noise, propulsion-efficient, nozzle exhaust systems being developed in NASA focused programs. These programs address a broad range of complex nozzle geometries operating in high temperature, compressible, flows. Seiner et al. previously discussed the jet configuration examined here. This convergent-divergent nozzle with an exit diameter of 3.6 inches was designed for an exhaust Mach number of 2.0 and a total temperature of 1680 F. The acoustic and aerodynamic data reported by Seiner et al. covered a range of jet total temperatures from 104 F to 2200 F at the fully-expanded nozzle pressure ratio. The aerodynamic data included centerline mean velocity and total temperature profiles. Computations were performed independently with two computational fluid dynamics (CFD) codes, ISAAC and PAB3D. Turbulence models employed include the k-epsilon model, the Gatski-Speziale algebraic-stress model and the Girimaji model, with and without the Sarkar compressibility correction. Centerline values of mean velocity and mean temperature are compared with experimental data.
Two dimensional quasi-three dimensional Navier-Stokes solvers were used to predict the static and dynamic airload characteristics of airfoils. The following three turbulence models were used: the Baldwin-Lomax algebraic model, the Johnson-King ODE model for maximum turbulent shear stress, and a two equation k-e model with law-of-the-wall boundary conditions. It was found that in attached flow the three models have good agreement with experimental data. In unsteady separated flows, these models give only a fair correlation with experimental data.