Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Kolmogorov Behavior of Near-Wall Turbulence and Its Application in Turbulence Modeling

The near-wall behavior of turbulence is re-examined in a way different from that proposed by Hanjalic and Launder and followers. It is shown that at a certain distance from the wall, all energetic large eddies will reduce to Kolmogorov eddies (the smallest eddies in turbulence). All the important wall parameters, such as friction velocity, viscous length scale, and mean strain rate at the wall, are characterized by Kolmogorov microscales. According to this Kolmogorov behavior of near-wall turbulence, the turbulence quantities, such as turbulent kinetic energy, dissipation rate, etc. at the location where the large eddies become Kolmogorov eddies, can be estimated by using both direct numerical simulation (DNS) data and asymptotic analysis of near-wall turbulence. This information will provide useful boundary conditions for the turbulent transport equations. As an example, the concept is incorporated in the standard k-epsilon model which is then applied to channel and boundary flows. Using appropriate boundary conditions (based on Kolmogorov behavior of near-wall turbulence), there is no need for any wall-modification to the k-epsilon equations (including model constants). Results compare very well with the DNS and experimental data.

Shih, Tsan-Hsing↗
NASA NTRS1992

On turbulent spots in plane Poiseuille flow

Turbulence characteristics inside a turbulent spot in plume Poiseuille flow are investigated by analyzing a data base obtained from a direct numerical simulation. The spot is found to consist of two distinct regions - a turbulent area and a wave area. The flow inside the turbulent area has a strong resemblance to that found in the fully developed turbulent channel. Suitably defined mean and rms fluctuations as well as the internal shear layer structures are found to be similar to the turbulent counterpart. In the wave area, the inflexional mean spanwise profiles cause a rapid growth of oblique waves, which break down to turbulence. The breakdown process of the oblique waves is reminiscent of the secondary instability observed during transition to turbulence in channel and boundary layer flows. Other detailed characteristics associated with the Poiseuille spot are presented and are compared with experimental results.

Henningson, Dan S.↗
NASA NTRS1992

Systematic study of Reynolds stress closure models in the computations of plane channel flows

The roles of pressure-strain and turbulent diffusion models in the numerical calculation of turbulent plane channel flows with second-moment closure models are investigated. Three turbulent diffusion and five pressure-strain models are utilized in the computations. The main characteristics of the mean flow and the turbulent fields are compared against experimental data. All the features of the mean flow are correctly predicted by all but one of the Reynolds stress closure models. The Reynolds stress anisotropies in the log layer are predicted to varying degrees of accuracy (good to fair) by the models. None of the models could predict correctly the extent of relaxation towards isotropy in the wake region near the center of the channel. Results from the directional numerical simulation are used to further clarify this behavior of the models.

Demuren, A. O.↗
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
NASA NTRS1992

Modeling of the reactant conversion rate in a turbulent shear flow

Results are presented of direct numerical simulations (DNS) of spatially developing shear flows under the influence of infinitely fast chemical reactions of the type A + B yields Products. The simulation results are used to construct the compositional structure of the scalar field in a statistical manner. The results of this statistical analysis indicate that the use of a Beta density for the probability density function (PDF) of an appropriate Shvab-Zeldovich mixture fraction provides a very good estimate of the limiting bounds of the reactant conversion rate within the shear layer. This provides a strong justification for the implementation of this density in practical modeling of non-homogeneous turbulent reacting flows. However, the validity of the model cannot be generalized for predictions of higher order statistical quantities. A closed form analytical expression is presented for predicting the maximum rate of reactant conversion in non-homogeneous reacting turbulence.

Frankel, S. H.↗
NASA NTRS1992

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗
NASA NTRS1992

Low Reynolds number multiple-time-scale turbulence model and calculations of steady and pulsating shear layers

A low Reynolds number multiple-time-scale turbulence model (LMS) and its application to fully developed turbulent channel flows and pulsating pipe flows are presented. The LMS can describe the inequilibrium turbulence phenomena down to the viscous sublayer. The calculated fluid flow and turbulence fields for the channel flows are in better agreement with the direct numerical simulation (DNS) results than those obtained using a Reynolds stress turbulence model, and the calculated near-wall dissipation rates are in qualitatively correct agreement with the DNS results. The LMS also successfully predicts the rapidly varying phase-lead of the wall shearing stress that occurs in a narrow range of the dimensionless frequency omega (+) = (omega x nu)/(upsilon x tau(exp2)) for the pulsating pipe flows while various other turbulence models fail to predict this phenomenon, and the LMS yields significantly improved numerical results for a wide range of the dimensionless frequency compared with those obtained using a rapid distortion theory (RDT).

Kim, Sang-Wook↗
NASA NTRS1990

Turbulent thermal convection in a differentially rotating channel

Differentially rotating disks of gases and solids occur in several astrophysical systems, in particular in the inner parts of protostellar nebulae, of which our own solar system is thought to be a relic. The objectives of this paper are to: (1) study localized turbulence in circumstances approximating those found in accretion disks using previously existing expertise in performing direct numerical simulation of turbulent, incompressible channel flows with low Reynolds numbers; (2) determine the limitations of such calculations; and (3) extend the type of numerical simulation (e.g., to include density and stratification and compressibility effects and to accommodate higher Reynolds numbers with sub-grid scale modeling) so that the relevant physical effects are realistically captured.

Cabot, W.↗
NASA NTRS1990

Modeling the turbulent kinetic energy equation for compressible, homogeneous turbulence

The turbulent kinetic energy transport equation, which is the basis of turbulence models, is investigated for homogeneous, compressible turbulence using direct numerical simulations performed at CTR. It is shown that the partition between dilatational and solenoidal modes is very sensitive to initial conditions for isotropic decaying turbulence but not for sheared flows. The importance of the dilatational dissipation and of the pressure-dilatation term is evidenced from simulations and a transport equation is proposed to evaluate the pressure-dilatation term evolution. This transport equation seems to work well for sheared flows but does not account for initial condition sensitivity in isotropic decay. An improved model is proposed.

Aupoix, B.↗
NASA NTRS1990

Laminar flamelet modeling of turbulent diffusion flames

In modeling turbulent combustion, decoupling the chemistry from the turbulence is of great practical significance. In cases in which the equilibrium chemistry model breaks down, laminar flamelet modeling (LFM) is a promising approach to decoupling. Here, the validity of this approach is investigated using direct numerical simulation of a simple chemical reaction in two-dimensional turbulence.

Mell, W. E.↗
NASA NTRS1992

A new time scale based k-epsilon model for near wall turbulence

A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

Yang, Z.↗
NASA NTRS1991

Analysis of turbulence in the orthonormal wavelet representation

The usefulness of the wavelet transform for the analysis of turbulent flow fields is explored by examining the wavelet transform properties of a decomposition of turbulent velocity fields into modes that exhibit the localization in a wavenumber and physical space. The calculations are performed on 3D fields from direct numerical simulations of isotropic flow and homogeneous shear flow, and from measurements in two laboratory wind-tunnel experimental velocity signals (boundary layer and wake behind a circular cylinder). The analysis confirmed that there is strong spatial intermittency in nonlinear quantities; their mean spectral behavior results from a delicate balance between large positive and negative excursions. The wavelet analysis is a way to quantify these observations in a standardized fashion by using 'flow-independent eddies' to decompose the velocity field.

Meneveau, Charles↗
NASA NTRS1992

Low Mach number two-dimensional hydrodynamic turbulence - Energy budgets and density fluctuations in a polytropic fluid

Theory suggests that three distinct types of turbulence can occur in the low Mach number limit of polytropic flow: nearly incompressible flows dominated by vorticity, nearly pure acoustic turbulence dominated by compression, and flows characterized by near statistical equipartition of vorticity and compressions. Distinctions between these kinds of turbulence are investigated here by direct numerical simulation of two-dimensional compressible hydrodynamic turbulence. Dynamical scalings of density fluctuations, examination of the ratio of transverse to longitudinal velocity fluctuations, and spectral decomposition of the fluctuations are employed to distinguish the nature of these low Mach number solutions. A strong dependence on the initial data is observed, as well as a tendency for enhanced effects of compressibility at later times and at higher wave numbers, as suggested by theories of nearly incompressible flows.

Ghosh, S.↗
NASA NTRS1991

An eddy viscosity calculation method for a turbulent duct flow

The mean velocity profile across a fully developed turbulent duct flow is obtained from an eddy viscosity relation combined with an empirical outer region wake function. Results are in good agreement with experiments and with direct numerical simulations in the same flow at two Reynolds numbers. In particular, the near-wall trend of the Reynolds shear stress and its variation with Reynolds number are similar to those of the simulations. The eddy viscosity method is more accurate than previous mixing length or implicit function methods.

Antonia, R. A.↗
NASA NTRS1992

Instability mode interactions in a spatially developing plane wake

Direct numerical simulations are used to study the development of various instability modes in a spatially developing 2D wake. Five types of forcing of the inlet are investigated: fundamental mode, fundamental and one or two subharmonics, fundamental mode and random noise, and random noise only. Statistical analyses are carried out, and some numerical results are compared with experimental measurements.

Maekawa, Hiroshi↗
NASA NTRS1992

Secondary instability mechanisms in compressible, axisymmetric boundary layers

Secondary instability mechanisms in compressible, axisymmetric boundary layers are analyzed using spectrally-accurate mean-flow and stability codes. Results show that subharmonic disturbances are the most dangerous secondary disturbances in an environment with a low to moderate intensity of the primary disturbance. The relation between spatial and temporal analyses of the secondary disturbance is explored at Mach 1.6 along a flat plate and Mach 6.8 along a cone. Spatial direct numerical simulations are utilized to confirm the quantitative predictions from spatial secondary instability theory.

Ng, Lian L.↗
NASA NTRS1992

Secondary instabilities in compressible boundary layers

The secondary instability mechanisms in a two-dimensional (2-D) compressible boundary layer over a flat plate have been examined for a range of Mach numbers up to 4.5. For subsonic and low supersonic flows, fundamental resonance dominates when the amplitude of the (2-D) primary disturbance is high, while subharmonic resonance prevails in an environment with a low primary disturbance. At a high supersonic Mach number of 4.5, the secondary instability of the second-mode primary is stronger than that of the first-mode primary. Further, the subharmonic and the combination resonance modes, which are slightly detuned from the subharmonic, are the dominant instabilities. The influence of the propagation direction of the primary disturbance on secondary instabilities is investigated at Mach 1.6. Whereas the fundamental and the subharmonic disturbances propagate synchronously with a 2-D primary wave, this is not true for an oblique primary. A subset of the secondary instability results is verified against direct numerical simulations. Some comparisons between spatial and temporal secondary instabilities have been made at Mach 1.6.

Ng, Lian L.↗
NASA NTRS1991

Compressibility effects on the growth and structure of homogeneous turbulent shear flow

Direct numerical simulations of compressible homogeneous turbulent shear flow are used to provide insight into compressibility effects on turbulence. The simulations show a reduction in the growth rate of the turbulence compared to the incompressible case. Examination of the turbulent kinetic energy budget shows that the reduced growth rate is due to an increase in the dissipation rate due to the divergence of the velocity and to the pressure-dilatation correlation which acts to transfer energy between internal energy and kinetic energy. The structure of the turbulence is also examined. Visualizations of the flowfields reveal the presence of eddy shocklets. These shock structures are important contributors to the increased dissipation of compressible turbulence. A mechanism for the generation of the shocks is suggested.

Blaisdell, G. A.↗