Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence modeling”

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.

At least 199 records · Page 11

An Initial Exploration of Improved Numerics within the Guidelines of the Negative Spalart-Allmaras Turbulence Model

A simple modification to the negative Spalart-Allmaras turbulence model is suggested so that when the turbulence working variable, ~v, is negative, diagonal dominance is increased, as is the tendency for the time-advancement scheme to push e toward positive values. Owing to the fact that the modification is only active when ~v is less than zero, the physical model is left unchanged. Using the proposed modification with a strong implicit solver based on Newton's method, convergence rates can be somewhat improved, with typical reductions in iterations and computer time on the order of 15-50%. The benefits are realized primarily when second- or higher-order accuracy is used for discretizing the convective terms in the turbulence model because of large overshoots that can occur with these schemes at the edges of boundary layers and wakes. For flowfields with few regions of negative ~v, or on very fine meshes where ~v is always greater than zero, little or no benefits should be expected.

Anderson, William K.↗

Turbulence modeling in hypersonic inlets

A study is conducted to analyze the performance of different turbulence models when applied to flow through a Mach 7.4 hypersonic inlet. The analysis, which is two-dimensional, is done by comparing computational results from a Parabolized Navier-Stokes code and a full Navier-Stokes code, with experimental data. The McDonald-Camarata (MC) and Baldwin-Lomax (BL) models were the two zero-equation models used in the study. The Turbulent Kinetic Energy (TKE) model was chosen as a representative higher order model. The MC model, when run with transition of flow, provides a solution which compares excellently with the data. Transition has a first order effect on the overall solution provided by the code. The BL model predicts separation of flow in the inlet, which contradicts experimental findings. The TKE model does not perform any better than the MC and BL models, despite the fact that it is a higher order turbulence model. The BL and TKE models predict transition in the inlet at a location which is much earlier than observed in the experiment. This may be attributed to the empirical constants used to determine the point of transition.

Ng, W. F.↗

NASA Trapezoidal Wing Computations Including Transition and Advanced Turbulence Modeling

Flow about the NASA Trapezoidal Wing is computed with several turbulence models by using grids from the first High Lift Prediction Workshop in an effort to advance understanding of computational fluid dynamics modeling for this type of flowfield. Transition is accounted for in many of the computations. In particular, a recently-developed 4-equation transition model is utilized and works well overall. Accounting for transition tends to increase lift and decrease moment, which improves the agreement with experiment. Upper surface flap separation is reduced, and agreement with experimental surface pressures and velocity profiles is improved. The predicted shape of wakes from upstream elements is strongly influenced by grid resolution in regions above the main and flap elements. Turbulence model enhancements to account for rotation and curvature have the general effect of increasing lift and improving the resolution of the wing tip vortex as it convects downstream. However, none of the models improve the prediction of surface pressures near the wing tip, where more grid resolution is needed.

Rumsey, C. L.↗

Comparison of turbulence models for powered-lift flow fields

The development of a practical turbulence model for the complex flows with strong streamline curvature associated with powered-lift aircraft was investigated. The approach was to modify existing models that were compatible with the widely used kappa-epsilon model already available in an existing 3D Navier-Stokes code. The state-of-the-art of current turbulence models was examined. Several models that provide curvature corrections to the kappa-epsilon model were numerically investigated. An algebraic Reynolds stress model, which provides additional physics, was also examined. Test cases including an impinging jet, an impinging jet in a crossflow, and two adjacent impinging jets with the resultant upwash fountain were calculated. The first two test cases could be modeled reasonsbly well with some of the models. The upwash fountain could not be modeled with sufficient accuracy, which is consistent with the findings of other investigators.

Roberts, Donald W.↗

Evaluation of Full Reynolds Stress Turbulence Models in FUN3D for Axisymmetric Jet Flow

Two full seven-equation turbulence models available in the FUN3D code are evaluated for their ability to improve the computation of challenging mixing flows encountered in aerospace propulsion. These models are the SSG/LRR and Wilcox full second-moment Reynolds stress models. They solve equations for the six components of the Reynolds stress and a seventh equation for the turbulent length scale. Two standard eddy viscosity models are also evaluated for comparison, the Spalart-Allmaras (SA) one-equation model and the Menter Shear Stress Transport (SST-V) two-equation turbulence model. Flow through an axisymmetric reference nozzle is examined at three flow conditions: subsonic unheated, subsonic heated, and near sonic unheated. Centerline profiles of velocity and turbulent kinetic energy and radial profiles of velocity, turbulent kinetic energy and turbulent stresses are examined. Results showed that the SA model did well at predicting the jet potential core length, but over-mixed the downstream flow, whereas the SST-V model over-predicted the potential core length. The Wilcox-model significantly over-predicted the potential core length and under-predicted the mixing and was not well-suited for the jet flows evaluated, however the SSG/LRR Reynolds stress model did well at predicting the mixing rate and mean velocity for all cases examined.

mixing layers↗

The influence of turbulence models on computer-simulated aircraft landing

Three models of simulated atmospheric turbulence were used to determine their respective influence on a computer-simulated, DC-8 type, aircraft landing. Statistical results based on a sample of 20 simulated landings illustrate the influence of turbulence simulated by the different models. For the same spectrum function, a non-Gaussian turbulence model produces the higher degree of touchdown variability. When different spectrum functions are used the von Karman model induces a smaller standard deviation in position error than the Dryden model. The turbulence energy contained by the higher frequency fluctuation has little effect on the landing position.

Frost, W.↗

A Modified Mixing Length Turbulence Model for Zero and Adverse Pressure Gradients

The modified mixing length (MML) turbulence model was installed in the Proteus Navier-Stokes code, then modified to make it applicable to a wider range of flows typical of aerospace propulsion applications. The modifications are based on experimental data for three flat-plate flows having zero, mild adverse, and strong adverse pressure gradients. Three transonic diffuser test cases were run with the new version of the model in order to evaluate its performance. All results are compared with experimental data and show improvements over calculations made using the Baldwin-Lomax turbulence model, the standard algebraic model in Proteus.

Conley, Julianne M.↗

Effects of turbulence models on the prediction of transonic wing flows

An investigation of the effects of turbulence models on the prediction of transonic wing flows is performed. The turbulence models used in this study are the equilibrium model of Baldwin and Lomax, and the original and modified models of Johnson and King. Comparisons with experimental data are presented which show clearly that the modified Johnson-King model works much better than the equilibrium model.

Abid, Ridha↗

A two-equation turbulence model for two-phase flows

A two-equation turbulence model has been developed for predicting two-phase flows. The two equations describe the conservation of turbulence kinetic energy and dissipation rate of that energy for the carrier fluid in a two-phase flow. They have been derived rigorously from the momentum equations of the carrier fluid. Closure of the time-mean equations is achieved by modeling the turbulent correlations up to third order. The new model eliminates the need to simulate in an ad hoc manner the effects of the dispersed phase on turbulence structure. Preliminary testing indicates that the model is successful in predicting the main features of a round gaseous jet laden with uniform-size solid particles.

Elghobashi, S. E.↗

Turbulence Modeling: Progress and Future Outlook

Progress in the development of the hierarchy of turbulence models for Reynolds-averaged Navier-Stokes codes used in aerodynamic applications is reviewed. Steady progress is demonstrated, but transfer of the modeling technology has not kept pace with the development and demands of the computational fluid dynamics (CFD) tools. An examination of the process of model development leads to recommendations for a mid-course correction involving close coordination between modelers, CFD developers, and application engineers. In instances where the old process is changed and cooperation enhanced, timely transfer is realized. A turbulence modeling information database is proposed to refine the process and open it to greater participation among modeling and CFD practitioners.

Marvin, Joseph G.↗

Turbulence Modeling Validation, Testing, and Development

The primary objective of this work is to provide accurate numerical solutions for selected flow fields and to compare and evaluate the performance of selected turbulence models with experimental results. Four popular turbulence models have been tested and validated against experimental data often turbulent flows. The models are: (1) the two-equation k-epsilon model of Wilcox, (2) the two-equation k-epsilon model of Launder and Sharma, (3) the two-equation k-omega/k-epsilon SST model of Menter, and (4) the one-equation model of Spalart and Allmaras. The flows investigated are five free shear flows consisting of a mixing layer, a round jet, a plane jet, a plane wake, and a compressible mixing layer; and five boundary layer flows consisting of an incompressible flat plate, a Mach 5 adiabatic flat plate, a separated boundary layer, an axisymmetric shock-wave/boundary layer interaction, and an RAE 2822 transonic airfoil. The experimental data for these flows are well established and have been extensively used in model developments. The results are shown in the following four sections: Part A describes the equations of motion and boundary conditions; Part B describes the model equations, constants, parameters, boundary conditions, and numerical implementation; and Parts C and D describe the experimental data and the performance of the models in the free-shear flows and the boundary layer flows, respectively.

Bardina, J. E.↗

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work.

Turbulent Boundary Layer↗

Nonequilibrium turbulence modeling effects on transonic vortical flows about delta wings

The Johnson-King turbulence model that is a viable method for calculating two dimensional transonic separated flows was extended into three dimensions. The implementation was done for Navier Stokes flow solvers written in general curvilinear coordinates. The present approach used in turbulence modeling is based on streamwise integration of an ordinary differential equation (o.d.e.) that governs the maximum Reynolds shear stress behavior. Streamwise integration of the o.d.e. approach was found to offer great mathematical simplicity and economy for three dimensional Navier Stokes methods. Thus, the new method is quick, simple, and very cheap. The new method was first checked against the data of a well known transonic axisymmetric bump experiment, and a good agreement was obtained. Later, the new method was used to compute the flow around a low aspect ratio wing in a transonic wind tunnel. Finally it was employed to study the nonequilibrium turbulence effects on the transonic vortical flows about a 65 deg sweep round leading edge delta wing.

Unver Kaynak↗

A Study of Two-Equation Turbulence Models on the Elliptic Streamline Flow

Several two-equation turbulence models are compared to data from direct numerical simulations (DNS) of the homogeneous elliptic streamline flow, which combines rotation and strain. The models considered include standard two-equation models and models with corrections for rotational effects. Most of the rotational corrections modify the dissipation rate equation to account for the reduced dissipation rate in rotating turbulent flows, however, the DNS data shows that the production term in the turbulent kinetic energy equation is not modeled correctly by these models. Nonlinear relations for the Reynolds stresses are considered as a means of modifying the production term. Implications for the modeling of turbulent vortices will be discussed.

Blaisdell, Gregory A.↗

A k-Omega Turbulence Model for Quasi-Three-Dimensional Turbomachinery Flows

A two-equation k-omega turbulence model has been developed and applied to a quasi-three-dimensional viscous analysis code for blade-to-blade flows in turbomachinery. the code includes the effects of rotation, radius change, and variable stream sheet thickness. The flow equations are given and the explicit runge-Kutta solution scheme is described. the k-omega model equations are also given and the upwind implicit approximate-factorization solution scheme is described. Three cases were calculated: transitional flow over a flat plate, a transonic compressor rotor, and transonic turbine vane with heat transfer. Results were compared to theory, experimental data, and to results using the Baldwin-Lomax turbulence model. The two models compared reasonably well with the data and surprisingly well with each other. Although the k-omega model behaves well numerically and simulates effects of transition, freestream turbulence, and wall roughness, it was not decisively better than the Baldwin-Lomax model for the cases considered here.

Chima, Rodrick V.↗

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗

Modeling turbulent, reacting flow

Several of the approximations or models involved in the development of a numerical combustor flow code are examined. In the first section, the importance of numerical accuracy is illustrated, and the impact that improved-accuracy schemes have on slowing convergence is demonstrated. Solution algorithms that can speed convergence are discussed and some performance features of these algorithms are illustrated. A sample calculation displaying the importance of boundary conditions on a three-dimensional numerical prediction is presented. The inaccuracy of a current turbulence model in highly turbulent (nonequilibrium) regions is described. Finally, the surprisingly good performance of a six-flux model in describing radiation heat transfer is displayed. In all the areas examined, continued research is still needed, but valuable engineering tools are available today.

Claus, Russell W.↗

On the validation of a code and a turbulence model appropriate to circular control airfoils

A computer code for calculating flow about a circulation control airfoil within a wind tunnel test section was developed. This code is being validated for eventual use as an aid to design such airfoils. The concept of code validation being used is explained. The initial stages of the process were accomplished. The present code was applied to a low subsonic, 2-D flow about a circulation control airfoil for which extensive data exist. Two basic turbulence models and variants thereof were successfully introduced into the algorithm, the Baldwin-Lomax algebraic and the Jones-Launder two equation models of turbulence. The variants include adding a history of the jet development for the algebraic model and adding streamwise curvature effects for both models. Numerical difficulties and difficulties in the validation process are discussed. Turbulence model and code improvements to proceed with the validation process are also discussed.

Viegas, J. R.↗