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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

CFD applications in pump flows

The objective of the proposed paper is to develop a computational procedure that solves incompressible Navier-Stokes equations for pump flows. The solution method is based on the pseudo-compressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. The equations are solved in steadily rotating reference frames and the centrifugal force and the Coriolis force are added to the equation of motion. As a benchmark problem, the flow through the Rocketdyne inducer is numerically simulated. A coarse grid solution is obtained with a single zone by using an algebraic turbulence model. In multi-zone fine grid computation, a one-equation Baldwin-Barth turbulence model is utilized. Numerical results are compared with experimental measurements and a good agreement is found between the two. The resulting computer code is then applied to the flow analysis inside a two-stage fuel pump impeller operating at 80 percent, 100 percent, and 120 percent of design flow.

Kiris, Cetin↗
NASA NTRS1992

Navier-Stokes analysis of blunt trailing edge airfoils

The flow around blunt trailing edge airfoils was studied by solving the Reynolds-averaged Navier-Stokes equations. The solution procedure combines a grid around the airfoil with a second grid for the wake so that the time advancement over the domain is fully implicit. This is not only very efficient for the algorithm but also allows implicit solutions of a one equation turbulence model appropriate for both boundary layers and wakes. An algebraic and two one-equation turbulence models are tested for a blunt RAE 2822 airfoil section and detailed comparisons with experimental data are presented in the trailing edge region.

Stanaway, Sharon↗
NASA NTRS2004

A Hybrid Approach To Tandem Cylinder Noise

Aeolian tone generation from tandem cylinders is predicted using a hybrid approach. A standard computational fluid dynamics (CFD) code is used to compute the unsteady flow around the cylinders, and the acoustics are calculated using the acoustic analogy. The CFD code is nominally second order in space and time and includes several turbulence models, but the SST k - omega model is used for most of the calculations. Significant variation is observed between laminar and turbulent cases, and with changes in the turbulence model. A two-dimensional implementation of the Ffowcs Williams-Hawkings (FW-H) equation is used to predict the far-field noise.

Lockard, David P.↗
NASA NTRS1985

On hairpin vortices as model of wall turbulence structure

A model of the hairpin vortex has been constructed and used in two distinct but related approaches. The first approach is kinematic in nature in which a synthesis procedure using hairpin vortices to provide a quantitative link between mean flow quantities and the statistical quantities of near wall turbulence has become developed. The second approach is dynamic in nature, and the evolution of an incipient 'representative' hairpin vortex as well as the distortion of a background laminar boundary layer flow, in which the hairpin vortex is immersed, has been simulated by numerical solution of the unsteady, three-dimensional Navier-Stokes equations.

Liu, N.-S.↗
NASA NTRS2001

Measurement of the Turbulence Kinetic Energy Budget of a Turbulent Planar Wake Flow in Pressure Gradients

Turbulence kinetic energy (TKE) is a very important quantity for turbulence modeling and the budget of this quantity in its transport equation can provide insight into the flow physics. Turbulence kinetic energy budget measurements were conducted for a symmetric turbulent wake flow subjected to constant zero, favorable and adverse pressure gradients in year-three of research effort. The purpose of this study is to clarify the flow physics issues underlying the demonstrated influence of pressure gradient on wake development and provide experimental support for turbulence modeling. To ensure the reliability of these notoriously difficult measurements, the experimental procedure was carefully designed on the basis of an uncertainty analysis. Four different approaches, based on an isotropic turbulence assumption, a locally axisymmetric homogeneous turbulence assumption, a semi-isotropy assumption and a forced balance of the TKE equation, were applied for the estimate of the dissipation term. The pressure transport term is obtained from a forced balance of the turbulence kinetic energy equation. This report will present the results of the turbulence kinetic energy budget measurement and discuss their implication on the development of strained turbulent wakes.

Liu, Xiao-Feng↗
NASA NTRS1989

Comparison of LDV measurements and Navier-Stokes solutions in a two-dimensional 180-degree turn-around duct

Results from an experimental and numerical investigation of turbulent subsonic flow inside a two-dimensional, strongly curved, 180 deg turn-around duct are presented. Data measured with a two-component, two-color laser Doppler velocimeter include profiles of mean axial velocity and local flow angle. Static pressure distributions are also measured. Results are obtained at a Mach number of 0.1, and at Reynolds numbers of 100,000 and a million based on channel height. Numerical calculations are performed using the incompressible Navier-Stokes equations with a Prandtl mixing length zero-equation turbulence model modified for internal flows. Theory and experiment are compared to evaluate the ability of the turbulence model to predict this class of internal flows with strong curvature.

Monson, Daryl J.↗
NASA NTRS1990

Embedded function methods for supersonic turbulent boundary layers

The development of embedded functions to represent the mean velocity and total enthalpy distributions in the wall layer of a supersonic turbulent boundary layer is considered. The asymptotic scaling laws (in the limit of large Reynolds number) for high speed compressible flows are obtained to facilitate eventual implementation of the embedded functions in a general prediction method. A self-consistent asymptotic structure is derived, as well as a compressible law of the wall in which the velocity and total enthalpy are logarithmic within the overlap zone, but in the Howarth-Dorodnitsyn variable. Simple outer region turbulence models are proposed (some of which are modifications of existing incompressible models) to reflect the effects of compressibility. As a test of the methodology and the new turbulence models, a set of self-similar outer region profiles is obtained for constant pressure flow; these are then coupled with embedded functions in the wall layer. The composite profiles thus obtained are compared directly with experimental data and good agreement is obtained for flows with Mach numbers up to 10.

He, J.↗
NASA NTRS2011

Low Boom Configuration Analysis with FUN3D Adjoint Simulation Framework

Off-body pressure, forces, and moments for the Gulfstream Low Boom Model are computed with a Reynolds Averaged Navier Stokes solver coupled with the Spalart-Allmaras (SA) turbulence model. This is the first application of viscous output-based adaptation to reduce estimated discretization errors in off-body pressure for a wing body configuration. The output adaptation approach is compared to an a priori grid adaptation technique designed to resolve the signature on the centerline by stretching and aligning the grid to the freestream Mach angle. The output-based approach produced good predictions of centerline and off-centerline measurements. Eddy viscosity predicted by the SA turbulence model increased significantly with grid adaptation. Computed lift as a function of drag compares well with wind tunnel measurements for positive lift, but predicted lift, drag, and pitching moment as a function of angle of attack has significant differences from the measured data. The sensitivity of longitudinal forces and moment to grid refinement is much smaller than the differences between the computed and measured data.

Park, Michael A.↗
NASA NTRS1986

Numerical studies of porous airfoils in transonic flow

A numerical tool is constructed to examine the effects of a porous surface on transonic airfoil performance and to help understand the flow structure of passive shockwave/boundary layer interactions. The porous region is located near the shock with a cavity underneath it. This study is composed of two parts. Solved in the first part, with an inviscid-flow approach, is the transonic full-potential equation associated with transpiration boundary conditions which are obtained from porosity modeling. The numerical results indicate that a porous airfoil has a wave drag lower than that of a solid airfoil. The observed lambda-shock structure in the wind-tunnel testing can be predicted. Furthermore, the lift could be increased with an appropriate porosity distribution. In the second part of this work, the modified version of either an interactive boundary layer (IBL) algorithm or a thin-layer Navier-Stokes (TLNS) algorithm is used to study the outer flow, while a stream-function formulation is used to model the inner flow in the shallow cavity. The coupling procedure at the porous surface is based on Darcy's law and the assumption of a constant total pressure in the cavity. In addition, a modified Baldwin-Lomax turbulence model is used to describe the transpired turbulent boundary layer in the TLNS approach, while the Cebeci turbulence model is used in the IBL approach. According to the present analysis, a porous surface can reduce the wave drag appreciably, but can also increase the viscous losses. As has been observed experimentally, the numerical results indicate that the total drag is reduced at higher Mach numbers and increased at lower Mach numbers when the angles of attack are small. Furthermore, the streamline pattern of passive shock/boundary layer interaction are revealed.

Chow, C. Y.↗
NASA NTRS1987

Computation of viscous transonic flow over porous airfoils

The viscous effects on transonic flow past an airfoil which contains a shallow cavity beneath a porous surface are studied numerically. The porous region occupies a small portion of the total airfoil surface, and is located near the shock. Both an interactive boundary layer (IBL) algorithm and a thin-layer Navier-Stokes (TLNS) algorithm have been modified for use in studying the outer flow, whereas a stream-function formulation has been used to model the inner flow in the small cavity. The coupling procedure at the porous surface is based on Darcy's law and on the assumption of a constant total presusre in the cavity. In addition, a modified Baldwin-Lomax turbulence model is used to consider the transpired turbulent boundary layer in the TLNS approach, and the Cebeci-Smith turbulence model is used in the IBL approach. According to the present analysis, a porous surface can reduce the wave drag appreciably, but it can also increase viscous losses. As has been observed experimentally, the numerical results indicate that the total drag is reduced at higher Mach numbers and increased at lower Mach numbers when the angles of attack are small. Furthermore, the streamline patterns of passive-shock and boundary-layer interaction are revealed in this study.

Chow, Chuen-Yen↗
NASA NTRS1986

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The objective of this investigation was the experimental evaluation of bulk turbulence and boundary thickness influence on the secondary flow development in a square, 90 deg turning duct. A three dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity in the large curved duct facility. The three dimensional development of the viscous shear layers in the curved duct has a strong influence on the complete flow field. Since ducted three dimensional flows are found in many engineering applications, including gas turbine engines, and contain high turbulence levels and high wall heat transfer rates, they present a difficult challenge to computational fluid mechanics codes. Turbulence modeling remains one of constraints to CFD advance due to inadequate physical understanding and experimental definition of turbulent shear flows. The results of this investigation expand the curved duct data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Crawford, Roger A.↗
NASA NTRS1990

Influence of bulk turbulence and entrance boundary layer thickness on the curved duct flow field

The objective of this investigation was the experimental evaluation of bulk turbulence and boundary thickness influence on the secondary flow development in a square, 90 deg turning duct. A three dimensional laser velocimetry system was utilized to measure the mean and fluctuating components of velocity in the large curved duct facility. The three dimensional development of the viscous shear layers in the curved duct has a strong influence on the complete flow field. Since ducted three dimensional flows are found in many engineering applications, including gas turbine engines, and contain high turbulence levels and high wall heat transfer rates, they present a difficult challenge to computational fluid mechanics codes. Turbulence modeling remains one of constraints to CFD advance due to inadequate physical understanding and experimental definition of turbulent shear flows. The results of this investigation expand the curved duct data base to higher turbulence levels and thicker entrance boundary layers. The experimental results provide a challenging benchmark data base for computational fluid dynamics code development and validation. The variation of inlet bulk turbulence intensity provides additional information to aid in turbulence model evaluation.

Shohadaee, Ahmad A.↗
NASA NTRS2013

Toward Verification of USM3D Extensions for Mixed Element Grids

The unstructured tetrahedral grid cell-centered finite volume flow solver USM3D has been recently extended to handle mixed element grids composed of hexahedral, prismatic, pyramidal, and tetrahedral cells. Presently, two turbulence models, namely, baseline Spalart-Allmaras (SA) and Menter Shear Stress Transport (SST), support mixed element grids. This paper provides an overview of the various numerical discretization options available in the newly enhanced USM3D. Using the SA model, the flow solver extensions are verified on three two-dimensional test cases available on the Turbulence Modeling Resource website at the NASA Langley Research Center. The test cases are zero pressure gradient flat plate, planar shear, and bump-inchannel. The effect of cell topologies on the flow solution is also investigated using the planar shear case. Finally, the assessment of various cell and face gradient options is performed on the zero pressure gradient flat plate case.

Pandya, Mohagna J.↗
DOE OSTI2026

Model Sensitivities of Biomass-Burning Aerosol Chemical Aging, Sulfate Formation, and Cloud Droplet Activation in the Southeastern Atlantic Using CESM and E3SM

Biomass-burning smoke drives large uncertainty in climate projections of the Earth's radiative balance. This is due to the chemical and physical evolution of smoke and its impact on clouds and radiation. Here we focus on the southeastern Atlantic region and its inflow of African biomass-burning smoke during August 2017. We evaluate smoke properties and processes in two coupled earth-system models, the Energy Exascale Earth System Model (E3SM) and Community Earth System Model (CESM). These are compared against in situ aircraft observations from two field campaigns, ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) and CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017 (CLARIFY-2017). Observations reveal an increase and subsequent decrease in smoke mean diameter, and a steady decrease in the mass ratio of organic aerosol (OA) to black carbon aerosol (BC) (OA:BC) over 4–12 days of aging, neither captured by the base models. Implementation of a photolytic loss scheme for secondary organic aerosol (SOA)—as a proxy for other heterogeneous volatilization chemistry—and a ∼1-day conversion for primary OA to SOA significantly improves the representation of this loss. In the boundary layer, both models show dimethyl sulfide driving a large increase in the sulfate aerosol mass fraction from the free troposphere, which is consistent with observations. Finally, models tend to underpredict cloud droplet number concentration partially due to weak modeled turbulent updraft strength, and model performance improves when the parameterized turbulent updraft strength is increased substantially. These results are expected to provide insights into future model development to reduce climate model uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Optimal white-noise stochastic forcing for linear models of turbulent channel flow

In the present study an optimisation problem is formulated to determine the forcing of an eddy-viscosity-based linearised Navier–Stokes model in channel flow at $Re_\tau \approx 5200$ ( $Re_\tau$ is the friction Reynolds number), where the forcing is white-in-time and spatially decorrelated. The objective functional is prescribed such that the forcing drives a response to best match a set of velocity spectra from direct numerical simulation (DNS), as well as remaining sufficiently smooth. Strong quantitative agreement is obtained between the velocity spectra from the linear model with optimal forcing and from DNS, but only qualitative agreement between the Reynolds shear stress co-spectra from the model and DNS. The forcing spectra exhibit a level of self-similarity, associated with the primary peak in the velocity spectra, but they also reveal a non-negligible amount of energy spent in phenomenologically mimicking the non-self-similar part of the velocity spectra associated with energy cascade. By exploiting linearity, the effect of the individual forcing components is assessed and the contributions from the Orr mechanism and the lift-up effect are also identified. Finally, the effect of the strength of the eddy viscosity on the optimisation performance is investigated. The inclusion of the eddy viscosity diffusion operator is shown to be essential in modelling of the near-wall features, while still allowing the forcing of the self-similar primary peak. In particular, reducing the strength of the eddy viscosity results in a considerable increase in the near-wall forcing of wall-parallel components.

Mechanics↗
NASA NTRS1997

Predicting NonInertial Effects with Algebraic Stress Models which Account for Dissipation Rate Anisotropies

Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coriolis-modified eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the turbulent model validation. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.

Jongen, T.↗
NASA NTRS1989

Further Improvements to Nozzle Boundary Layer Calculations in BLIMPJ

Further improvements made to advance the current Boundary Layer Integral Matrix Procedure - Version J (BLIMPJ) containing previously modeled simplified calculation methods by accounting for condensed phase, thick boundary layer and free stream turbulence effects are discussed. The condensed phase effects were included through species composition effect considered via input to the code and through particle damping effect considered via a turbulence model. The thrust loss calculation procedure for thick boundary layer effects was improved and the optimization of net thrust with respect to nozzle length was performed. The effects of free stream turbulence were approximately modeled in the turbulence model.

Praharaj, S. C.↗
NASA NTRS1986

Turbulent Boundary Layers in Oscillating Flows. Part 1: an Experimental and Computational Study

An experimental-computational study of the behavior of turbulent boundary layers for oscillating air flows over a plane surface with a small favorable mean pressure gradient is described. Experimental studies were conducted for boundary layers generated on the test section wall of a facility that produces a flow with a mean free stream velocity and a superposed nearly-pure sinusoidal component over a wide range of frequency. Flow at a nominal mean free stream velocity of 50 m/s were studied at atmospheric pressure and temperature at selected axial positions over a 2 m test length for frequencies ranging from 4 to 29 Hz. Quantitative experimental results are presented for unsteady velocity profiles and longitudinal turbulence levels obtained from hot wire anemometer measurements at three axial positions. Mean velocity profiles for oscillating flows were found to exhibit only small deviations from corresponding steady flow profiles, while amplitudes and phase relationships exhibited a strong dependence on axial position and frequency. Since sinusoidal flows could be generated over a wide range of frequency, studies at fixed values of reduced frequency at different axial positions were studied. Results show that there is some utility in the use of reduced frequency to correlate unsteady velocity results. The turbulence level u' sub rms was observed to vary essentially sinusoidally around values close to those measured in steady flow. However, the amplitude of oscillation and phase relations for turbulence level were found to be strongly frequency dependent. Numerical predictions were obtained using an unsteady boundary layer computational code and the Cebeci-Smith and Glushko turbulence models. Predicted quantities related to unsteady velocity profiles exhibit fair agreement with experiment when the Cebeci-Smith turbulence model is used.

Cook, W. J.↗