SEARCH · Engineering Papers

Results for “Direct numerical simulation”

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 577 records · Page 32

All sources NASA NTRS DOE OSTI

NASA NTRS1993

A research program for improving heat transfer prediction for the laminar to turbulent transition region of turbine vanes/blades

A program sponsored by NASA for the investigation of the heat transfer in the transition region of turbine vanes and blades with the objective of improving the capability for predicting heat transfer is described. The accurate prediction of gas-side heat transfer is important to the determination of turbine longevity, engine performance, and developmental costs. The need for accurate predictions will become greater as the operating temperatures and stage loading levels of advanced turbine engines increase. The present methods for predicting transition shear stress and heat transfer on turbine blades are based on incomplete knowledge and are largely empirical. To meet the objective of the NASA program, a team approach consisting of researchers from government, universities, a research institute, and a small business is presented. The research is divided into the areas of experiments, direct numerical simulations (DNS), and turbulence modeling. A summary of the results to date is given for the above research areas in a high-disturbance environment (bypass transition) with a discussion of the model development necessary for use in numerical codes.

Simon, Frederick F.↗

NASA NTRS1993

Parameterization of subgrid-scale stress by the velocity gradient tensor

The objective of this work is to construct and evaluate subgrid-scale models that depend on both the strain rate and the vorticity. This will be accomplished by first assuming that the subgrid-scale stress is a function of the strain and rotation rate tensors. Extensions of the Caley-Hamilton theorem can then be used to write the assumed functional dependence explicitly in the form of a tensor polynomial involving products of the strain and rotation rates. Finally, use of this explicit expression as a subgrid-scale model will be evaluated using direct numerical simulation data for homogeneous, isotropic turbulence.

Lund, Thomas S.↗

NASA NTRS1993

Application of incremental unknowns to the Burgers equation

In this article, we make a few remarks on the role that attractors and inertial manifolds play in fluid mechanics problems. We then describe the role of incremental unknowns for approximating attractors and inertial manifolds when finite difference multigrid discretizations are used. The relation with direct numerical simulation and large eddy simulation is also mentioned.

Choi, Haecheon↗

NASA NTRS1993

Reynolds stress closure modeling in wall-bounded flows

This report describes two projects. Firstly, a Reynolds stress closure for near-wall turbulence is described. It was motivated by the simpler k-epsilon-(v-bar(exp 2)) model described in last year's annual research brief. Direct Numerical Simulation of three-dimensional channel flow shows a curious decrease of the turbulent kinetic energy. The second topic of this report is a model which reproduces this effect. That model is described and used to discuss the relevance of the three dimensional channel flow simulation to swept wing boundary layers.

Durbin, Paul A.↗

NASA NTRS1993

The three-dimensional evolution of a plane wake

In the past three decades, linear stability analysis has led to a comprehensive understanding of the linear stages of transition in plane wakes. Our understanding of the nonlinear and turbulent stages is less developed. Nonlinear theory developed by Papageorgiou and Smith was used to study the long-wavelength regime in wakes. The nonlinear and turbulent stages were investigated experimentally, and few numerical studies examined the early nonlinear stages of forced wakes. The evolution of three dimensional disturbances in an incompressible wake is investigated using direct numerical simulations. The instantaneous three-dimaensional structures and corresponding statistics are presented.

Maekawa, H.↗

NASA NTRS1993

Mixing in a stratified shear flow: Energetics and sampling

Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed for Richardson numbers from 0 to 1 and for Prandtl numbers between 0.1 and 2. The results indicate that mixing efficiency R(sub f) varies with turbulent Froude number in a manner consistent with laboratory experiments performed with Prandtl numbers of 0.7 and 700. However, unlike the laboratory results, for a particular Froude number, the simulations do not show a clear dependence on the magnitude of R(sub f) on Pr. The observed maximum value of R(sub f) is 0.25. When averaged over vertical length scales of an order of magnitude greater than either the overturning or Ozmidov scales of the flow, the simulations indicate that the dissipation rate epsilon is only weakly lognormally distributed with an intermittency of about 0.01 whereas estimated values in the ocean are 3 to 7.

Ivey, G. N.↗

NASA NTRS1993

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud↗

NASA NTRS1993

Study and modeling of finite rate chemistry effects in turbulent non-premixed flames

The development of numerical models that reflect some of the most important features of turbulent reacting flows requires information about the behavior of key quantities in well defined combustion regimes. In turbulent flames, the coupling between turbulent and chemical processes is so strong that it is extremely difficult to isolate the role played by one individual physical phenomenon. Direct numerical simulation (hereafter DNS) allows us to study in detail the turbulence-chemistry interaction in some restricted but completely defined situations. Globally, non-premixed flames are controlled by two limiting regimes: the fast chemistry case, where the turbulent flame can be pictured as a random distribution of local chemical equilibrium problems; and the slow chemistry case, where the chemistry integrates in time the turbulent fluctuations. The Damkoehler number, ratio of a mechanical time scale to chemical time scale, is used to distinguish between these regimes. Today most of the industrial computer codes are able to perform predictions in the hypothesis of local equilibrium chemistry using a presumed shape for the probability density function (pdt) of the conserved scalar. However, the finite rate chemistry situation is of great interest because industrial burners usually generate regimes in which, at some points, the flame is undergoing local extinction or at least non-equilibrium situations. Moreover, this variety of situations strongly influences the production of pollutants. To quantify finite rate chemistry effect, the interaction between a non-premixed flame and a free decaying turbulence is studied using DNS. The attention is focused on the dynamic of extinction, and an attempt is made to quantify the effect of the reaction on the small scale mixing process. The unequal diffusivity effect is also addressed. Finally, a simple turbulent combustion model based on the DNS observations and tractable in real flow configurations is proposed.

Vervisch, Luc↗

NASA NTRS1992

On the local nature of the energy cascade

The local nature of the energy cascade in space and time is studied using direct numerical simulation of decaying and forced isotropic turbulence. To examine the concept that large scales evolve into smaller ones, we compute the Lagrangian correlation coefficient between local kinetic energy at different scales. This correlation is found to peak at a Lagrangian time-delay that increases with scale separation. The results show that, on average, the flow of energy to smaller scales is predominantly local in physical space and that the view of eddies decaying into smaller ones while transferring their kinetic energy appears to be, on average, quite realistic. To examine the spectral characteristics of the cascade under unsteady conditions, a pulse of large-scale energy is added to the large-eddy simulation of forced isotropic turbulence. As time progresses, the evolution of this pulse through bands of increasing wavenumbers is studied.

Meneveau, C.↗

NASA NTRS1992

A study of the fine scale motions of incompressible time-developing mixing layers

This work is an extension of a project conducted at the previous CTR summer program and was reported by Chen et al. (1990). In that program, the geometry and topology of the dissipating motions in a variety of shear flows was examined. All data was produced by direct numerical simulations (DNS). The partial derivatives of the velocity field were determined at every grid point in the flow and various invariants and related quantities were computed from the velocity gradient tensor. Motions characterized by high rates of kinetic energy dissipation and high enstrophy were of particular interest. Scatter diagrams of the invariants were mapped out and interesting and unexpected patterns were seen. Each type of shear layer produced its own characteristic scatter plot. In the present project, attention is focused on the incompressible plane mixing layer, and the scatter diagrams are replaced with more useful joint probability density contours. Comparison of the topology of the dissipating motions of flows at different Reynolds numbers are made. Also, plane mixing layers at the same Reynolds number but with different initial conditions are compared.

Soria, J.↗

NASA NTRS1992

Isotropy of small scale turbulence

The degree to which local isotropy is satisfied has been examined using direct numerical simulations for a fully developed channel flow. Attention is mainly given to the high wavenumber part of vorticity and temperature derivative spectra. The ratio of these spectra and their isotropic values depends on the particular quantity considered, the departure from isotropy being more pronounced for the temperature derivative than for the vorticity. When the Kolmogorov-normalized wavenumber is sufficiently large, isotropy is satisfied provided the (Kolmogorov-normalized) mean strain rate is sufficiently small. This result appears to be independent of the quantity considered and of the Reynolds number.

Antonia, R. A.↗

NASA NTRS1992

LES versus DNS: A comparative study

We have performed Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) of forced isotropic turbulence at moderate Reynolds numbers. The subgrid scale model used in the LES is based on an eddy viscosity which adjusts instantaneously the energy spectrum of the LES to that of the DNS. The statistics of the large scales of the DNS (filtered DNS field or fDNS) are compared to that of the LES. We present results for the transfer spectra, the skewness and flatness factors of the velocity components, the PDF's of the angle between the vorticity and the eigenvectors of the rate of strain, and that between the vorticity and the vorticity stretching tensor. The above LES statistics are found to be in good agreement with those measured in the fDNS field. We further observe that in all the numerical measurements, the trend was for the LES field to be more gaussian than the fDNS field. Future research on this point is planned.

Shtilman, L.↗

NASA NTRS1992

Investigation of the dilatational dissipation in compressible homogeneous shear flow

The dilatational dissipation rate within compressible homogeneous turbulent shear flow is studied using data from direct numerical simulations. It is found that the dilatational dissipation rate is mainly associated with large scale acoustic waves. Eddy shocklets are observed; however, they have little contribution to the average dissipation rate. A mechanism for the generation of eddy shocklets is shown to be the focusing of acoustic waves. Turbulence models for the dilatational dissipation rate are compared with data from the simulations. It is found that the formulation of Zeman used by Viegas & Rubesin to calculate a compressible mixing layer agrees well with the numerical simulation results for turbulent Mach numbers less than 0.3. However, it is also found that, for the Mach number range occurring in mixing layers, the model does not accurately represent the theory upon which it is based.

Blaisdell, G. A.↗

NASA NTRS1992

Evolution of the shock front and turbulence structures in the shock/turbulence interaction

The interaction of a weak shock front with isotropic turbulence has been investigated using Direct Numerical Simulation (DNS). Two problems were considered: the ability of the field equation (the equation for a propagating surface) to model the shock; and a quantitative study of the evolution of turbulence structure using the database generated by Lee et al. Field equation model predictions for front shape have been compared with DNS results; good agreement is found for shock wave interaction with 2D turbulence and for a single steady vorticity wave. In the interaction of 3D isotropic turbulence with a normal shock, strong alignment of vorticity with the intermediate eigenvector of the rate of strain tensor (S(sup *)(sub ij) = S(sub ij) - (1/3)(delta(sub ij))(S(sub kk))) is seen to develop upstream of the shock and to be further amplified on passage through the shock. Vorticity tends to align at 90 deg to the largest eigenvector, but there is no preferred alignment with the smallest eigenvector. Upstream of the shock, the alignments continue to develop even after the velocity derivative skewness saturates. There is a significant tendency, which increases with time throughout the computational domain, for velocity to align with vorticity. The alignment between velocity and vorticity is strongest in eddy regions and weakest in convergence regions.

Kevlahan, N.↗

NASA NTRS1992

A transport equation for eddy viscosity

A transport equation for eddy viscosity is proposed for wall bounded turbulent flows. The proposed model reduces to a quasi-homogeneous form far from surfaces. Near to a surface, the nonhomogeneous effect of the wall is modeled by an elliptic relaxation model. All the model terms are expressed in local variables and are coordinate independent; the model is intended to be used in complex flows. Turbulent channel flow and turbulent boundary layer flows with/without pressure gradient are calculated using the present model. Comparisons between model calculations and direct numerical simulation or experimental data show good agreement.

Durbin, P. A.↗

NASA NTRS1992

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗

NASA NTRS1992

Effect of finite-rate chemistry and unequal Schmidt numbers on turbulent non-premixed flames modeled with single-step chemistry

The interaction between a quasi-laminar flame and a turbulent flowfield is investigated through direct numerical simulations (DNS) of reacting flow in two- and three-dimensional domains. Effects due to finite-rate chemistry are studied using a single step global reaction A (fuel) + B (oxidizer) yields P (product), and by varying a global Damkoehler number, as a result of which the turbulence-chemistry interaction in the flame is found to generate a wide variety of conditions, ranging from near-equilibrium to near-extinction. Differential diffusion effects are studied by changing the Schmidt number of one reactive species to one-half. It is observed that laminar flamelet response is followed within the turbulent flowfield, except in regions where transient effects seem to dominate.

Chen, J. H.↗

NASA NTRS1992

Structure of turbulent non-premixed flames modeled with two-step chemistry

Direct numerical simulations of turbulent diffusion flames modeled with finite-rate, two-step chemistry, A + B yields I, A + I yields P, were carried out. A detailed analysis of the turbulent flame structure reveals the complex nature of the penetration of various reactive species across two reaction zones in mixture fraction space. Due to this two zone structure, these flames were found to be robust, resisting extinction over the parameter ranges investigated. As in single-step computations, mixture fraction dissipation rate and the mixture fraction were found to be statistically correlated. Simulations involving unequal molecular diffusivities suggest that the small scale mixing process and, hence, the turbulent flame structure is sensitive to the Schmidt number.

Chen, J. H.↗