SEARCH · Search NASA

Results for “Eddy viscosity”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2010

A Nonlinear Dynamic Subscale Model for Partially Resolved Numerical Simulation (PRNS)/Very Large Eddy Simulation (VLES) of Internal Non-Reacting Flows

A brief introduction of the temporal filter based partially resolved numerical simulation/very large eddy simulation approach (PRNS/VLES) and its distinct features are presented. A nonlinear dynamic subscale model and its advantages over the linear subscale eddy viscosity model are described. In addition, a guideline for conducting a PRNS/VLES simulation is provided. Results are presented for three turbulent internal flows. The first one is the turbulent pipe flow at low and high Reynolds numbers to illustrate the basic features of PRNS/VLES; the second one is the swirling turbulent flow in a LM6000 single injector to further demonstrate the differences in the calculated flow fields resulting from the nonlinear model versus the pure eddy viscosity model; the third one is a more complex turbulent flow generated in a single-element lean direct injection (LDI) combustor, the calculated result has demonstrated that the current PRNS/VLES approach is capable of capturing the dynamically important, unsteady turbulent structures while using a relatively coarse grid.

Shih, Tsan-Hsing↗

NASA NTRS1998

Incompressible Turbulent Wing-Body Junction Flow

The overall objective of this study is to contribute to the optimized design of fan bypass systems in advanced turbofan engines. Increasing the engine bypass ratios have provided a major boost in engine performance improvement over the last fifty years. An engine with high bypass ratio (11-16:1) such as the Advanced Ducted Propulsion (ADP) is being developed and is expected to provide an additional 25% improvement in overall efficiency over the early turbofans. Such significant improvements in overall efficiency would reduce the cost per seat mile, which is a major government and Industry challenge for the 21th century. The research is part of the Advanced Subsonic Technology (AST) program that involves a NASA, U.S. Industry and FAA partnership with the goal of a safe and highly productive global air transportation system. The immediate objective of the study is to perform numerical simulation of duct-strut interactions to elucidate the loss mechanisms associated with this configuration that is typical of advanced turbofan engines such as ADP. However, at present experimental data for a duct-strut configuration are not available. Thus, as a first step a wing-body junction flow would be studied and is the specific objective of the present study. At the outset it is to be recognized that while duct-strut interaction flow is similar to that of wing-body junction flows, there are some differences owing to the presence of a wall at both ends of the strut. Likewise, some differences are due to the sheared inflow (as opposed to a uniform inflow) velocity profile. It is however expected that some features of a wing-body junction flow would persist. Next, some of the salient aspects of the complex flow near a wing-body junction, as revealed by various studies reported in the literature will be reviewed. One of the principle characteristics of the juncture flow, is the presence of the mean flow components in a plane perpendicular to the direction of the oncoming free-stream flow. The lateral curvature of the wing/strat causes the oncoming turbulent layer to skew about am axis (x-axis) parallel to the plane (xz-plane) of the mean shear. This is the principle mechanism for the generation of secondary flow. Such skew-induced secondary flows are slow to be attenuated by Reynolds stresses. Additional contribution to the generation of secondary flow comes from anisotropies in Reynolds stresses. Upstream of the strut, the mean-vorticity is directed span wise (along the y-direction). The presence of secondary flow in the vicinity of the strut causes the vorticity to stretch around the obstacle in a horse-shoe shape, with each leg having a vorticity of the opposite sense. The blockage effect of the strut imposes a severe adverse pressure gradient on the oncoming turbulent shear layer, causing boundary layer separation ahead of the leading edge, resulting in a vortex that rolls up and flows downstream into the juncture region. The separation vortices trailing in the wake of the wing can alter the lift or drag characteristics of the surfaces downstream of the wing-body juncture. Likewise, on submarines, the wake flow behind the appendage can degrade the performance of the propeller located downstream. The complex nature of this flow is caused by the presence of all six components of Reynolds stresses. Devenport and Simpson report that in the vicinity of the horse-shoe vortex there is intense recirculation with turbulent stresses being much larger than those normally observed in turbulent flows. These features contribute to making this flow a challenge to predict numerically. Some of the past studies provide useful insights into this flow that would guide our numerical efforts. In measurements reported by Shabaka and Bradshaw, the eddy viscosity tensor is seen to be non-isotropic and has negative components in certain regions. In an effort to evaluate the closure assumptions of various turbulence models, Devenport and Simpson used their own extensive measurements in juncture flows around the nose of a wing-body junction. Measured values of mean-velocity and/or turbulence kinetic energy was used to predict the magnitude of the shear stress vector. Algebraic stress models performed the best followed by Cebeci-Smith eddy viscosity model. The flow is reported to be dominated by a pressure field produced by the wing and the velocity field generated by the horseshoe vortex that is wrapped around the junction between the wing and wall. Kubendran et al. conclude from an experimental study that the shape of leading edge of the wing as characterized by its slenderness ratio is a major factor in determining the flow fields in the juncture region. The more thinner the leading edge of the juncture, the weaker the horseshoe vortex is. Also, with a slender leading edge, the secondary flow in the juncture would be due mainly to the cross-stream gradients of Reynolds stresses rather than due to a lateral skewing of the shear layer.

Krishnamurthy, R.↗

NASA NTRS1974

Polar symmetric flow of a viscous compressible atmosphere; an application to Mars

The atmosphere is assumed to be driven by a polar symmetric temperature field and the equations of motion in pressure ratio coordinates are linearized by considering the zero order in terms of a thermal Rossby number R delta I/(2a omega) sq where delta T is a measure of the latitudinal temperature gradient. When the eddy viscosity is greater than 1 million sq cm/sec, the boundary layer extends far up into the atmosphere, making the geostrophic approximation invalid for the bulk of the atmosphere. A temperature model for Mars was used which was based on Mariner 9 infrared spectral data with a 30% increase in the depth averaged temperature from the winter pole to the subsolar point. The results obtained for the increase in surface pressure from the subsolar point to the winter pole, as a function of eddy viscosity and with no-slip conditions imposed at the surface, are given.

Pirraglia, J. A.↗

NASA NTRS2003

On the Computation of Space-Time Correlations by Large-Eddy Simulation

Numerical comparisons in decaying isotropic turbulence suggest that there exist discrepancies in time correlations evaluated by DNS and LES using eddy-viscosity-type SGS models. This is consistent with the previous observations in forced isotropic turbulence. Therefore, forcing is not the main cause of the discrepancies. Comparisons among different SGS models in the LES also indicate that the model choice affects the time correlations in the LES. The multi-scale LES method using the dynamic Smagorinsky model on the small scale equation is the most accurate of the all models, the classic Smagorinsky model is the least accurate and the dynamic Smagorinsky model and spectral eddy viscosity model give intermediate results with small differences. The generalized sweeping hypothesis implies that time correlations in decaying isotropic turbulence are mainly determined by the instantaneous energy spectra and sweeping velocities. The analysis based on the sweeping hypothesis explains the discrepancies in our numerical simulations: the LES overpredicts the decorrelation time scales because the sweeping velocities are smaller than the DNS values, and underpredicts the magnitudes of time correlations because the energy spectrum levels are lower than the DNS ones. Since the sweeping velocity is determined by the energy spectra, one concludes that an accurate prediction of the instantaneous energy spectra guarantees the accuracy of time correlations. An analytical expression of sound power spectra based on Lighthill's theory and the quasi-normal closure assumption suggests that the sound power spectra are sensitive to errors in time correlations. Small errors in time correlations can cause significant errors in the sound power spectra, which exhibit a sizable drop at moderate to high frequencies accompanied by a shift of the peaks to lower frequencies. Based on the above analysis, two possible ways to improve the acoustic power spectrum predictions can be considered. The first is to construct better SGS models to improve the LES accuracy for time correlations. The second is to remedy the temporal statistics of the Lighthill stress tensor in order to "recover" the contribution from the unresolved scales in LES to time correlations.

He, Guo-Wei↗

NASA NTRS1987

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.↗

NASA NTRS1987

A low dimensional dynamical system for the wall layer

Low dimensional dynamical systems which model a fully developed turbulent wall layer were derived.The model is based on the optimally fast convergent proper orthogonal decomposition, or Karhunen-Loeve expansion. This decomposition provides a set of eigenfunctions which are derived from the autocorrelation tensor at zero time lag. Via Galerkin projection, low dimensional sets of ordinary differential equations in time, for the coefficients of the expansion, were derived from the Navier-Stokes equations. The energy loss to the unresolved modes was modeled by an eddy viscosity representation, analogous to Heisenberg's spectral model. A set of eigenfunctions and eigenvalues were obtained from direct numerical simulation of a plane channel at a Reynolds number of 6600, based on the mean centerline velocity and the channel width flow and compared with previous work done by Herzog. Using the new eigenvalues and eigenfunctions, a new ten dimensional set of ordinary differential equations were derived using five non-zero cross-stream Fourier modes with a periodic length of 377 wall units. The dynamical system was integrated for a range of the eddy viscosity prameter alpha. This work is encouraging.

Aubry, N.↗

NASA NTRS1987

Near-wall k-epsilon turbulence modeling

Mansour, N. N.↗

NASA NTRS1972

A finite difference method for predicting supersonic turbulent boundary layer flows with tangential slot injection

An implicit finite difference method has been applied to tangential slot injection into supersonic turbulent boundary layer flows. In addition, the effects induced by the interaction between the boundary layer displacement thickness and the external pressure field are considered. In the present method, three different eddy viscosity models have been used to specify the turbulent momentum exchange. One model depends on the species concentration profile and the species conservation equation has been included in the system of governing partial differential equations. Results are compared with experimental data at stream Mach numbers of 2.4 and 6.0 and with results of another finite difference method. Good agreement was generally obtained for the reduction of wall skin friction with slot injection and with experimental Mach number and pitot pressure profiles. Calculations with the effects of pressure interaction included showed these effects to be smaller than effects of changing eddy viscosity models.

Miner, E. W.↗

NASA NTRS1975

Calculation of shock-separated turbulent boundary layers

Numerical solutions of the complete, time-averaged conservation equations using several eddy-viscosity models for the Reynolds shear stress to close the equations are compared with experimental measurements in a compressible, turbulent separated flow. An efficient time-splitting, explicit difference scheme was used to solve the two-dimensional conservation equations. The experiment used for comparison was a turbulent boundary layer that was separated by an incident shock wave in a Mach 2.93 flow with a unit Reynolds number of 5.7 x 10 to the seventh power m. Comparisons of predicted and experimental values of surface pressure, shear stress along the wall, and velocity profiles are shown. One of the tested eddy-viscosity models which allows the shear stress to be out of equilibrium with the mean flow produces substantially better agreement with the experimental measurements than the simpler models. A tool is thereby provided for inferring additional information about the flow, such as static pressures in the stream, which might not be directly obtainable from experiments.

Baldwin, B. S.↗

NASA NTRS1982

Computation of incompressible, three-dimensional turbulent boundary layers and comparison with experiment

Incompressible three dimensional, turbulent boundary layer (3DTBL) experiments were simulated numerically by integrating the boundary layer equations together with an algebraic eddy viscosity turbulence model. For the flow treated, the downstream portion, where the crossflow was large, was not predicted with the present computational method; the flow was significantly influenced by elliptic flow field effects. Departures from the boundary layer concept are indicated. Calculations agreed reasonably well with the mean flow development up to separation. In one experiment the normal pressure gradients were found to be neligible in regions with large skewing and allowed testing turbulence models using the boundary layer equations. The simulation of this flow compared favorably with the experimental data throughout the flow field and suggested the applicability of algebraic eddy viscosity models for 3DTBLs.

Mueller, U. R.↗

NASA NTRS1967

Turbulent mixing /analysis/

Approximate methods for solving differential equations of motion to predict free turbulent mixing processes - single-form eddy viscosity model

EDDY VISCOSITY↗

NASA NTRS1997

Prediction of High-Lift Flows using Turbulent Closure Models

The flow over two different multi-element airfoil configurations is computed using linear eddy viscosity turbulence models and a nonlinear explicit algebraic stress model. A subset of recently-measured transition locations using hot film on a McDonnell Douglas configuration is presented, and the effect of transition location on the computed solutions is explored. Deficiencies in wake profile computations are found to be attributable in large part to poor boundary layer prediction on the generating element, and not necessarily inadequate turbulence modeling in the wake. Using measured transition locations for the main element improves the prediction of its boundary layer thickness, skin friction, and wake profile shape. However, using measured transition locations on the slat still yields poor slat wake predictions. The computation of the slat flow field represents a key roadblock to successful predictions of multi-element flows. In general, the nonlinear explicit algebraic stress turbulence model gives very similar results to the linear eddy viscosity models.

Rumsey, Christopher L.↗

NASA NTRS1986

Computation of turbulent boundary layer flows with an algebraic stress turbulence model

An algebraic stress turbulence model is presented, characterized by the following: (1) the eddy viscosity expression is derived from the Reynolds stress turbulence model; (2) the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale; and (3) the diffusion coefficients for turbulence equations are adjusted so that the kinetic energy profile extends further into the free stream region found in most experimental data. The turbulent flow equations were solved using a finite element method. Examples include: fully developed channel flow, fully developed pipe flow, flat plate boundary layer flow, plane jet exhausting into a moving stream, circular jet exhausting into a moving stream, and wall jet flow. Computational results compare favorably with experimental data for most of the examples considered. Significantly improved results were obtained for the plane jet flow, the circular jet flow, and the wall jet flow; whereas the remainder are comparable to those obtained by finite difference methods using the standard kappa-epsilon turbulence model. The latter seems to be promising with further improvement of the expression for the eddy viscosity coefficient.

Kim, Sang-Wook↗

NASA NTRS1995

Bleed Roughness Mechanism

Suction applied through wall bleed is expected to reduce boundary-layer thickness and lead to a profile that can resist separation. However, experiments with rows of suction holes or slots have shown that the benefits of suction are not realized under certain circumstances. Lee, Sloan and Paynter attributed this to a 'bleed roughness', and modeled the phenomena using a modified inner-length scale with convective lag for the eddy viscosity. Test information was required to determine the roughness values at each Mach number and bleed rate. The present work used numerical calculation, subsonic and supersonic experimentation, and physical modeling were used to explore bleed roughness. It was found that the phenomena can exist in the absence of bleed, that it is a function of geometry and Mach number, and that it does not appear in low subsonic flow. The inner-length scale/convective lag approach was found to be inconsistent with recent physical observation. Alternate steady models for the roughness effect disagreed with mass flux trends and the choked flow limit, and suggested that non-steady mechanisms are present. While new supersonic experiments did exhibit low-frequency plenum oscillations, high-frequency unsteadiness associated with the orifice flow scale appears to be the cause of bleed roughness. In fact, boundary-layer calculations with an oscillatory orifice flow produced the same effect in the time-averaged downstream profile for incompressibile flow. A new empirical turbulence-model was proposed which increased the eddy viscosity throughout the entire boundary-layer by a constant which was calibrated using the zero-bleed profile. Preliminary calculations agreed with the experimental data. The model significantly reduced the test information required to determine roughness values for turbulence modeling, and provided prediction capability for supersonic bleed flows. An abstract of a thesis entitled 'Boundary-Layer Bleed Roughness' is included.

Russell, David A.↗

NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗

NASA NTRS2015

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Simulation↗

NASA NTRS1990

A turbulent two-phase flow model for nebula flows

A new and very efficient turbulent two-phase flow numericaly model is described to analyze the environment of a protoplanetary nebula at a stage prior to the formation of planets. Focus is on settling processes of dust particles in flattened gaseous nebulae. The model employs a perturbation technique to improve the accuracy of the numerical simulations of such flows where small variations of physical quantities occur over large distance ranges. The particles are allowed to be diffused by gas turbulence in addition to settling under gravity. Their diffusion coefficients is related to the gas turbulent viscosity by the non-dimensional Schmidt number. The gas turbulent viscosity is determined by the means of the eddy viscosity hypothesis that assumes the Reynolds stress tensor proportional to the mean strain rate tensor. Zero- and two-equation turbulence models are employed. Modeling assumptions are detailed and discussed. The numerical model is shown to reproduce an existing analytical solution for the settling process of particles in an inviscid nebula. Results of nebula flows are presented taking into account turbulence effects of nebula flows. Diffusion processes are found to control the settling of particles.

Champney, Joelle M.↗

NASA NTRS1990

Three-dimensional compressible jet-in-crossflow calculations using improved viscosity models and adapted grid

Previous calculations of jet-in-crossflow problems have been sensitive to the turbulence and artificial viscosity models and to the grid. Consequently, the eddy viscosity model in the PARC3D code was modified to consider the turbulent jet by switching from the Baldwin-Lomax (1978) model to an axisymmetric jet model. A modified artificial viscosity model has been utilized and evaluated in this study as well. The new model includes cell size scaling and a directional dependence in the coefficients. Computational results from PARC3D demonstrate the effects of the viscosity models on the pressure distribution fore and aft of the jet and the ability of the adaptive grid scheme to adjust to the three-dimensional gradients around the jet.

Lytle, John K.↗