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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
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.↗
NASA NTRS1991

Numerical study of axial turbulent flow over long cylinders

The effects of transverse curvature are investigated by means of direct numerical simulations of turbulent axial flow over cylinders. Two cases of Reynolds number of about 3400 and layer-thickness-to-cylinder-radius ratios of 5 and 11 were simulated. All essential turbulence scales were resolved in both calculations, and a large number of turbulence statistics were computed. The results are compared with the plane channel results of Kim et al. (1987) and with experiments. With transverse curvature the skin friction coefficient increases and the turbulence statistics, when scaled with wall units, are lower than in the plane channel. The momentum equation provides a scaling that collapses the cylinder statistics, and allows the results to be interpreted in light of the plane channel flow. The azimuthal and radial length scales of the structures in the flow are of the order of the cylinder diameter. Boomerang-shaped structures with large spanwise length scales were observed in the flow.

Neves, J. C.↗
NASA NTRS1991

Examination of various turbulence models for application in liquid rocket thrust chambers

There is a large variety of turbulence models available. These models include direct numerical simulation, large eddy simulation, Reynolds stress/flux model, zero equation model, one equation model, two equation k-epsilon model, multiple-scale model, etc. Each turbulence model contains different physical assumptions and requirements. The natures of turbulence are randomness, irregularity, diffusivity and dissipation. The capabilities of the turbulence models, including physical strength, weakness, limitations, as well as numerical and computational considerations, are reviewed. Recommendations are made for the potential application of a turbulence model in thrust chamber and performance prediction programs. The full Reynolds stress model is recommended. In a workshop, specifically called for the assessment of turbulence models for applications in liquid rocket thrust chambers, most of the experts present were also in favor of the recommendation of the Reynolds stress model.

Hung, R. J.↗
NASA NTRS1992

Probability density functions in turbulent channel flow

The probability density functions (pdf's) of the fluctuating velocity components, as well as their first and second derivatives, are calculated using data from the direct numerical simulations (DNS) of fully developed turbulent channel flow. It is observed that, beyond the buffer region, the pdf of each of these quantities is independent of the distance from the channel wall. It is further observed that, beyond the buffer region, the pdf's for all the first derivatives collapse onto a single universal curve and those of the second derivatives also collapse onto another universal curve, irrespective of the distance from the wall. The kinetic-energy dissipation rate exhibits log normal behavior.

Dinavahi, Surya P. G.↗
NASA NTRS1992

Modeling of turbulence and transition

The first objective is to evaluate current two-equation and second order closure turbulence models using available direct numerical simulations and experiments, and to identify the models which represent the state of the art in turbulence modeling. The second objective is to study the near-wall behavior of turbulence, and to develop reliable models for an engineering calculation of turbulence and transition. The third objective is to develop a two-scale model for compressible turbulence.

Shih, Tsan-Hsing↗
NASA NTRS1992

Near-wall reconstruction of higher order moments and length scales using the POD

An analysis of the near-wall behavior of the proper orthogonal decomposition (POD) eigenfunctions derived from direct numerical simulation (DNS) of channel flow is performed. Consistent with previous studies, a low order multi-mode reconstruction of the kinetic energy and Reynolds shear stress suffices. A similar reconstruction of the isotropic dissipation rate is shown to be insufficient, however. An analysis is performed of the multi-mode composition of the dissipation rate in the near-wall region, and it is shown that a significant number of higher-order modes are required to achieve the correct asymptotic consistency in the near-wall region. In an attempt to avoid this problem, a length scale definition is proposed in terms of an integration of the correlation tensor which factors in the presence of the wall. The wall is accounted for by only integrating out to 2y(+) and not over the entire domain. Viscous and inviscid estimates for the dissipation were used in the near-wall and core regions respectively, in conjunction with this length scale representation to obtain an estimate of the dissipation throughout the domain. The resulting dissipation exhibits the proper behavior near the wall and in the inertial layer. A 1 POD mode estimate of the length scale is computed and found to agree quite well with the length scale obtained when the entire correlation tensor is used.

Glauser, Mark N.↗
NASA NTRS1992

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗
NASA NTRS1992

Propagation velocity and space-time correlation of perturbations in turbulent channel flow

A database obtained from direct numerical simulation of a turbulent channel flow is analyzed to extract the propagation velocity V of velocity, vorticity, and pressure fluctuations from their space-time correlations. A surprising result is that V is approximately the same as the local mean velocity for most of the channel, except for the near-wall region. For y(+) is less than or equal to 15, V is virtually constant, implying that perturbations of all flow variables propagate like waves near the wall. In this region V is 55 percent of the centerline velocity U(sub c) for velocity and vorticity perturbations and 75 percent of U(sub c) for pressure perturbations. Scale-dependence of V is also examined by analyzing the bandpass filtered flow fields. Comprehensive documentation of the propagation velocities and space-time correlation data, which should prove useful in the evaluation of Taylor's hypothesis is presented. An attempt was made to explain some of the data in terms of our current understanding of organized structures, although not all of the data can be explained this way.

Kim, John↗
NASA NTRS1992

The transition prediction toolkit: LST, SIT, PSE, DNS, and LES

The e(sup N) method for predicting transition onset is an amplitude ratio criterion that is on the verge of full maturation for three-dimensional, compressible, real gas flows. Many of the components for a more sophisticated, absolute amplitude criterion are now emerging: receptivity theory, secondary instability theory, parabolized stability equations approaches, direct numerical simulation and large-eddy simulation. This paper will provide a description of each of these new theoretical tools and provide indications of their current status.

Zang, Thomas A.↗
NASA NTRS1992

The remarkable ability of turbulence model equations to describe transition

This paper demonstrates how well the k-omega turbulence model describes the nonlinear growth of flow instabilities from laminar flow into the turbulent flow regime. Viscous modifications are proposed for the k-omega model that yield close agreement with measurements and with Direct Numerical Simulation results for channel and pipe flow. These modifications permit prediction of subtle sublayer details such as maximum dissipation at the surface, k approximately y(exp 2) as y approaches 0, and the sharp peak value of k near the surface. With two transition specific closure coefficients, the model equations accurately predict transition for an incompressible flat-plate boundary layer. The analysis also shows why the k-epsilon model is so difficult to use for predicting transition.

Wilcox, David C.↗
NASA NTRS1993

Propagation velocity of perturbations in turbulent channel flow

A database obtained from direct numerical simulation of a turbulent channel flow is analyzed to extract the streamwise component of the propagation velocity V of velocity, vorticity, and pressure fluctuations from their space-time correlations. A surprising result is that V is approximately the same as the local mean velocity for most of the channel, except for the near-wall region. For y(+) less than 15, V is virtually constant, implying that perturbations of all flow variables propagate like waves near the wall. In this region V is 55 percent of the centerline velocity Uc for velocity and vorticity perturbations and 75 percent of U sub c for pressure perturbations. This is equal to U at y(+) = 15 for velocity and vorticity perturbations, and equal to U at y(+) = 20 for pressure perturbations, indicating that the dynamics of the nearwall turbulence is controlled by turbulence structures present near y(+) about 15-20. Scale dependence of V is also examined by analyzing the bandpass-filtered flow fields. This paper contains comprehensive documentation on the propagation velocities, which should prove useful in the evaluation of Taylor's hypothesis.

Kim, John↗
NASA NTRS1993

Spanwise scale selection in plane mixing layers

Direct numerical simulations of temporally evolving plane mixing layers undergoing as many as three pairings have been examined for evidence of spanwise scale change. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. The amplitude of the initial three-dimensional disturbances varied from infinitesimal to large enough to trigger a rapid transition to turbulence. The time required for a change of characteristic spanwise scale with infinitesimal three-dimensional disturbances was found to be very long, requiring three or more pairings to complete a doubling of the spanwise scale. Stronger three-dimensionality can produce more rapid scale changes, but it is also likely to trigger transition to turbulence. No evidence was found for a change from an organized array of rib vortices at one spanwise scale to a similar array at a larger spanwise scale.

Rogers, Michael M.↗
NASA NTRS1992

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

The presently obtained closed-form analytical expressions, which predict the limiting rate of mean reactant conversion in homogeneous turbulent flows under the influence of a binary reaction, are derived via the single-point pdf method based on amplitude mapping closure. With this model, the maximum rate of the mean reactant's decay can be conveniently expressed in terms of definite integrals of the parabolic cylinder functions. The results obtained are shown to be in good agreement with data generated by direct numerical simulations.

Madnia, C. K.↗
NASA NTRS1993

Effective equations and the inverse cascade theory for Kolmogorov flows

We study the two-dimensional Kolmogorov flows in the limit as the forcing frequency goes to infinity. Direct numerical simulation indicates that the low frequency energy spectrum evolves to a universal kappa (exp -4) decay law. We derive effective equations governing the behavior of the large scale flow quantities. We then present numerical evidence that with smooth initial data, the solution to the effective equation develops a kappa (exp -4) type singularity at a finite time. This gives a convenient explanation for the kappa (exp -4) decay law exhibited by the original Kolmogorov flows.

Weinan, E.↗
NASA NTRS1993

Perspective - 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 NTRS1993

The rollup of a vortex layer near a wall

The behavior of an inviscid vortex layer of non-zero thickness near a wall is studied, both through direct numerical simulation of the two-dimensional vorticity equation at high Reynolds numbers, and using an approximate ordinary nonlinear integro-differential equation which is satisfied in the limit of a thin layer under long-wavelength perturbations. For appropriate initial conditions the layer rolls up and breaks into compact vortices which move along the wall at constant speed. Because of the effect of the wall, they correspond to equilibrium counter-rotating vortex dipoles. This breakup can be related to the disintegration of the initial conditions of the approximate nonlinear dispersive equation into solitary waves. The study is motivated by the formation of longitudinal vortices from vortex sheets in the wall region of a turbulent channel.

Jimenez, Javier↗
NASA NTRS1993

A Reynolds stress model for near-wall turbulence

The paper formulates a tensorially consistent near-wall second-order closure model. Redistributive terms in the Reynolds stress equations are modeled by an elliptic relaxation equation in order to represent strongly nonhomogeneous effects produced by the presence of walls; this replaces the quasi-homogeneous algebraic models that are usually employed, and avoids the need for ad hoc damping functions. The model is solved for channel flow and boundary layers with zero and adverse pressure gradients. Good predictions of Reynolds stress components, mean flow, skin friction, and displacement thickness are obtained in various comparisons to experimental and direct numerical simulation data. The model is also applied to a boundary layer flowing along a wall with a 90-deg, constant-radius, convex bend.

Durbin, P. A.↗
NASA NTRS1993

Degrees of locality of energy transfer in the inertial range

Measured raw transfer interactions from which local energy transfer is argued to result are summed in a way that directly indicates the scale disparity (s) of contributions to the net energy flux across the spectrum. It is found that the dependence upon s closely follows the s exp -4/3 form predicted by classical arguments. As a result, it is concluded that direct numerical simulation measurements lend support to the classical Kolmogorov phenomenology of local interactions and local transfer in an inertial range.

Zhou, YE↗