Search NASA⌕ Search

SEARCH · Search NASA

Results for “turbulence model”

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 109 records · Page 6

Turbulence modeling for computational aerodynamics

The status of turbulence modeling for external aerodynamic flows is reviewed, and closure concepts for the compressible form of the Reynolds-averaged Navier-Stokes equations are briefly outlined to establish a framework for comparison. The importance of experimental requirements for developing and verifying turbulence models is emphasized. Attention is then given to three important flow categories: attached flows, separating and reattaching flows, and trailing-edge flows. Examples of comparisons between experiments and computations for twoand three-dimensional flows are presented to illustrate the status of modeling. It is shown that, for most two-dimensional attached flows, eddy viscosity concepts will probably be adequate. For attached three-dimensional flows, however, eddy-viscosity and Reynolds-stress models are both deficient, failing to predict the proper surface shear-stress in rapidly skewing boundary layers, regardless of whether they are pressure driven or strain driven.

Marvin, J. G.↗

Atmospheric Turbulence Modeling for Aero Vehicles: Fractional Order Fits

Atmospheric turbulence models are necessary for the design of both inlet/engine and flight controls, as well as for studying coupling between the propulsion and the vehicle structural dynamics for supersonic vehicles. Models based on the Kolmogorov spectrum have been previously utilized to model atmospheric turbulence. In this paper, a more accurate model is developed in its representative fractional order form, typical of atmospheric disturbances. This is accomplished by first scaling the Kolmogorov spectral to convert them into finite energy von Karman forms and then by deriving an explicit fractional circuit-filter type analog for this model. This circuit model is utilized to develop a generalized formulation in frequency domain to approximate the fractional order with the products of first order transfer functions, which enables accurate time domain simulations. The objective of this work is as follows. Given the parameters describing the conditions of atmospheric disturbances, and utilizing the derived formulations, directly compute the transfer function poles and zeros describing these disturbances for acoustic velocity, temperature, pressure, and density. Time domain simulations of representative atmospheric turbulence can then be developed by utilizing these computed transfer functions together with the disturbance frequencies of interest.

Kopasakis, George↗

Turbulence modeling in noninertial frames of reference

The effect of an arbitrary change of frame on the structure of turbulence models is examined from a fundamental theoretical standpoint. It is proven, as a rigorous consequence of the Navier-Stokes equations, that turbulence models must be form invariant under arbitrary translational accelerations of the reference frame and should only be affected by rotations through the intrinsic mean vorticity. A direct application of the invariance property along with the Taylor-Proudman Theorem, material frame-indifference in the limit of two-dimensional turbulence and Rapid Distortion Theory is shown to yield powerful constraints on the allowable form of turbulence models. Most of the commonly used turbulence models are demonstrated to be in serious violation of these constraints and consequently are inconsistent with the Navier-Stokes equations in non-inertial frames. Alternative models with improved non-inertial properties are developed and some simple applications to rotating turbulent flows are considered.

Speziale, Charles G.↗

Implementation/validation of a low Reynolds number two-equation turbulence model in the Proteus Navier-Stokes code: Two-dimensional/axisymmetric

The implementation and validation of the Chien low Reynolds number k-epsilon turbulence model in the two dimensional axisymmetric version Proteus, a compressible Navier-Stokes computer code, are presented. The set of k-epsilon equations are solved by marching in time using a coupled alternating direction implicit (ADI) solution procedure with generalized first or second order time differencing. To validate Proteus and the k-epsilon turbulence model, laminar and turbulent computations were done for several benchmark test cases: incompressible fully developed 2-D channel flow; fully developed axisymmetric pipe flow; boundary layer flow over a flat plate; and turbulent Sajben subsonic transonic diffuser flows. Proteus results from these test cases showed good agreement with analytical results and experimental data. Detailed comparisons of both mean flow and turbulent quantities showed that the Chien k-epsilon turbulence model given good results over a wider range of turbulent flow than the Baldwin-Lomax turbulence model in the Proteus code with no significant CPU time penalty for more complicated flow cases.

Bui, Trong T.↗

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗

Towards a structure-based turbulence model

The author is developing a new type of turbulence model in which a new one-point quantity, the eddy structure tensor, carries information about the two-point structure of the turbulence. The model was motivated by the observation that conventional one-point turbulence models based only on the turbulent stresses do not predict the rapid changes in state that are found when anisotropic homogeneous turbulence is subjected to mean rotation, and hence are fundamentally incorrect for rotation. The model appears to give topologically correct predictions for the changes in stress state and structure state under all types of rapid distortions of homogeneous turbulence.

Reynolds, W. C.↗

Effect of Turbulence Models on Two Massively-Separated Benchmark Flow Cases

Two massively-separated flow cases (the 2-D hill and the 3-D Ahmed body) were computed with several different turbulence models in the Reynolds-averaged Navier-Stokes code CFL3D as part of participation in a turbulence modeling workshop held in Poitiers, France in October, 2002. Overall, results were disappointing, but were consistent with results from other RANS codes and other turbulence models at the workshop. For the 2-D hill case, those turbulence models that predicted separation location accurately ended up yielding a too-long separation extent downstream. The one model that predicted a shorter separation extent in better agreement with LES data did so only by coincidence: its prediction of earlier reattachment was due to a too-late prediction of the separation location. For the Ahmed body, two slant angles were computed, and CFD performed fairly well for one of the cases (the larger slant angle). Both turbulence models tested in this case were very similar to each other. For the smaller slant angle, CFD predicted massive separation, whereas the experiment showed reattachment about half-way down the center of the face. These test cases serve as reminders that state- of-the-art CFD is currently not a reliable predictor of massively-separated flow physics, and that further validation studies in this area would be beneficial.

Rumsey, Christopher L.↗

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

Near-wall k-epsilon turbulence modeling

The flow fields from a turbulent channel simulation are used to compute the budgets for the turbulent kinetic energy (k) and its dissipation rate (epsilon). Data from boundary layer simulations are used to analyze the dependence of the eddy-viscosity damping-function on the Reynolds number and the distance from the wall. The computed budgets are used to test existing near-wall turbulence models of the k-epsilon type. It was found that the turbulent transport models should be modified in the vicinity of the wall. It was also found that existing models for the different terms in the epsilon-budget are adequate in the region from the wall, but need modification near the wall. The channel flow is computed using a k-epsilon model with an eddy-viscosity damping function from the data and no damping functions in the epsilon-equation. These computations show that the k-profile can be adequately predicted, but to correctly predict the epsilon-profile, damping functions in the epsilon-equation are needed.

Mansour, N. N.↗

Verification and Validation of the k-kL Turbulence Model in FUN3D and CFL3D Codes

The implementation of the k-kL turbulence model using multiple computational uid dy- namics (CFD) codes is reported herein. The k-kL model is a two-equation turbulence model based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's equation is well suited for term-by-term mod- eling and displays useful features compared to other two-equation models. An important di erence is that this formulation leads to the inclusion of higher-order velocity derivatives in the source terms of the scale equations. This can enhance the ability of the Reynolds- averaged Navier-Stokes (RANS) solvers to simulate unsteady ows. The present report documents the formulation of the model as implemented in the CFD codes Fun3D and CFL3D. Methodology, veri cation and validation examples are shown. Attached and sepa- rated ow cases are documented and compared with experimental data. The results show generally very good comparisons with canonical and experimental data, as well as matching results code-to-code. The results from this formulation are similar or better than results using the SST turbulence model.

Abdol-Hamid, Khaled S.↗

Turbulence modeling in non-inertial frames of reference

The effect of an arbitrary change of frame on the structure of turbulence models is examined from a fundamental theoretical standpoint. It is proven, as a rigorous consequence of the Navier-Stokes equations, that turbulence models must be form invariant under arbitrary translational accelerations of the reference frame and should only be affected by rotations through the intrinsic mean vorticity. A direct application of the invariance property along with the Taylor-Proudman Theorem, material frame-indifference in the limit of two-dimensional turbulence and Rapid Distortion Theory is shown to yield powerful constraints on the allowable form of turbulence models. Most of the commonly used turbulence models are demonstrated to be in serious violation of these constraints and consequently are inconsistent with the Navier-Stokes equations in non-inertial frames. Alternative models with improved non-inertial properties are developed and some simple applications to rotating turbulent flows are considered.

Speziale, Charles G.↗

An Assessment of Turbulence Models for MSR Sample Retrieval Lander Aerothermal Design

Several cases from past entries in the T9 and T5 hypersonic test facilities are examined using various turbulence models and shock-layer thermochemistries with the objectives of verifying/validating turbulence models and determining appropriate multiplicative factors on predicted turbulent heating for use in aerothermal design of the heatshield of the Sample Retrieval Lander (SRL) of the Mars Sample Return (MSR) Flagship mission. The zero-equation turbulence models of Baldwin-Lomax and Cebeci-Smith are shown to replicate measurements to within experimental uncertainty for the T9 cases. For T5 cases, the zero-equation model of Baldwin-Lomax can also replicate measurements provided the shock-layer thermochemistry is modified. The modified thermochemistry resolves a dichotomy from past CFD that required different wall catalycity assumptions in the laminar and turbulent flow regimes. The present work shows that such assumptions are not necessary. Finally, an initial estimate for a multiplicative factor on turbulent heating is estimated for the SST turbulence model – the chosen model for SRL aerothermal design simulations.

Mars entry↗

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗

Algebraic turbulence modeling for unstructured and adaptive meshes

An algebraic turbulence model based on the Baldwin-Lomax model, has been implemented for use on unstructured grids. The implementation is based on the use of local background structured turbulence meshes. At each time-step, flow variables are interpolated from the unstructured mesh onto the background structured meshes, the turbulence model is executed on these meshes, and the resulting eddy viscosity values are interpolated back to the unstructured mesh. Modifications to the algebraic model were required to enable the treatment of more complicated flows, such as confluent boundary layers and wakes. The model is used in conjuction with an efficient unstructured multigrid finite-element Navier-Stokes solver in order to compute compressible turbulent flows on fully unstructured meshes. Solutions about single and multiple element airfoils are obtained and compared with experimental data.

Mavriplis, Dimitri J.↗

Advanced turbulence models for turbomachinery

Development and assessment of the single-time-scale k-epsilon turbulence model with different near-wall treatments and the multi-scale turbulence model for rotating flows are presented. These turbulence models are coded as self contained module decks that can be interfaced with a number of CFD main flow solvers. For each model, a stand-alone module deck with its own formulation, discretization scheme, solver and boundary condition implementations is presented. These satellite decks will take as input (from a main flow solver) the velocity field, grid, boundary condition specifications and will deliver turbulent quantities as output. These modules were tested as separate entities and, although many logical and programming problems were overcome, only wider use and further testing can render the modules sufficiently 'fool proof'.

Hadid, Ali H.↗

Progress and challenges in modeling turbulent aerodynamic flows

Progress in modeling external aerodynamic flows achieved by using computations and experiments designed to guide turbulence modeling is presented. The computational procedures emphasize utilization of the Reynolds-averaged Navier-Stokes equations and various statistical modeling approaches. Developments for including the influence of compressibility are provided; they point up some of the complexities involved in modeling high-speed flows. Examples of complementary studies that provide the status, limitations, and future challenges of modeling for transonic, supersonic, and hypersonic flows are given.

Marvin, Joseph G.↗

Recent Turbulence Model Advances Applied to Multielement Airfoil Computations

A one-equation linear turbulence model and a two-equation nonlinear explicit algebraic stress model (EASM) are applied to the flow over a multielement airfoil. The effect of the K-epsilon and K-omega forms of the two-equation model are explored, and the K-epsilon form is shown to be deficient in the wall-bounded regions of adverse pressure gradient flows. A new K-omega form of EASM is introduced. Nonlinear terms present in EASM are shown to improve predictions of turbulent shear stress behind the trailing edge of the main element and near midflap. Curvature corrections are applied to both the one- and two-equation turbulence models and yield only relatively small local differences in the flap region, where the flow field undergoes the greatest curvature. Predictions of maximum lift are essentially unaffected by the turbulence model variations studied.

Rumsey, Christopher L.↗

Heated Supersonic Axisymmetric Jet Cases for the NASA Turbulence Modeling Resource

Three new supersonic jet cases were built for the NASA Turbulence Modeling Resource (TMR). These cases were taken from the 6th AIAA Propulsion Aerodynamics Workshop (PAW) nozzle test problem, which utilized data taken at NASA GRC consisting of non-intrusive measurements of jet plume velocities and temperatures, including both mean values and turbulence statistics. The three jet cases all used a Mach 1.63 axisymmetric nozzle: a heated jet for each at on-design conditions, a temperature-matched on-design condition at Mach 1.63, and an off-design heated condition for the Mach 1.63 nozzle. Computational fluid dynamics (CFD) solutions are obtained using three established Reynolds-averaged Navier-Stokes (RANS) codes: Wind-US, FUN3D, and VULCAN-CFD, all using the Menter Shear Stress Transport k- turbulence model with vorticity source term (SST-V). The NASA TMR requires at least two CFD codes yielding essentially identical results to certify that the approaches are verified, meaning that they solve the posed CFD equations as intended. The three codes generally provided very close agreement with each other for jet plume quantities, with the only exception being static temperature. This discrepancy was determined to be the result of whether or not turbulent kinetic energy was considered in the definition of total energy. Validation of the SST-V turbulence model for this case (its level of agreement/disagreement with experimental data) is also addressed.

turbulance↗