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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

A numerical study of transition control by periodic suction-blowing

The applicability of active control of transition by periodic suction-blowing is investigated via direct numerical simulations of the Navier-Stokes equations. The time-evolution of finite-amplitude disturbances in plane channel flow is compared in detail with and without control. The analysis indicates that, for relatively small three dimensional amplitudes, a two dimensional control effectively reduces disturbance growth rates even for linearly unstable Reynolds numbers. After the flow goes through secondary instability, three dimensional control seems necessary to stabilize the flow. An investigation of the temperature field suggests that passive temperature contamination is operative to reflect the flow dynamics during transition.

Biringen, Sedat↗
NASA NTRS1987

Mechanisms by which heat release affects the flow field in a chemically reacting, turbulent mixing layer

The mechanisms by which heat release affects the fluid dynamics in a turbulent reacting mixing layer are studied by direct numerical simulation. In agreement with previous laboratory experiments, the heat release is observed to lower the rate at which the mixing layer grows and to reduce the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer are shown to produce changes in the flame vortex structure that act to produce more diffuse vortices than in the constant density case, resulting in lower rotation rates of fluid elements. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layers are shown to result from vorticity generation by baroclinic torques. The density reductions also lower the generation rates of turbulent kinetic energy and the turbulent shear stresses, resulting in less turbulent mixing of fluid elements. Calculations of the energy in the various wave numbers show that the heat release has a stabilizing effect on the growth rates of individual modes. A linear stability analysis of a simlified model problem confirms this, showing that low density fluid in the mixing region will result in a shift of the frequency of the unstable modes to lower wave numbers (longer wavelengths). The growth rates of the unstable modes decrease, contributing to the slower growth of the mixing layer.

Riley, J. J.↗
NASA NTRS1986

Numerical study of sink-flow boundary layers

Direct numerical simulations of sink-flow boundary layers, with acceleration parameters K between 1.5 x 10 to the -6th and 3.0 x 10 to the -6th, are presented. The three-dimensional, time-dependent Navier-Stokes equations are solved numerically, using a spectral method, with about one million degrees of freedom. The flow is assumed to be statistically steady, and self-similar. A multiple-scale approximation and periodic conditions are applied to the fluctuations. The turbulence is studied using instantaneous and statistical results. Good agreement with the experiments of Jones and Launder (1972) is observed. The two effects of the favorable pressure gradient are to extend the logarithmic layer, and to alter the energy balance of the turbulence near the edge of the boundary layer. At low Reynolds number the logarithmic layer is shortened and slightly displaced, but wall-layer streaks are present even at the lowest values of R(theta) for which turbulence can be sustained. Large quiescent patches appear in the flow. Relaminarization occurs at K = 3.0 x 10 to the -6th, corresponding to a Reynolds number R(theta) of about 330.

Spalart, Philippe R.↗
NASA NTRS1987

The structure of the vorticity field in homogeneous turbulent flows

The structures of the vorticity fields in several homogeneous irrotational straining flows and a homogeneous turbulent shear flow were examined using a database generated by direct numerical simulation of the unsteady Navier-Stokes equations. In all cases, strong evidence was found for the presence of coherent vortical structures. The initially isotropic vorticity fields were rapidly affected by imposed mean strain and the rotational component of mean shear and developed accordingly. In the homogeneous turbulent shear-flow cases, the roll-up of mean vorticity into characteristic hairpin vortices was clearly observed, supporting the view that hairpin vortices are an important vortical structure in all turbulent shear flows; the absence of mean shear in the homogeneous irrotational straining flows precludes the presence of hairpin-like vortices.

Rogers, Michael M.↗
NASA NTRS1987

Turbulence statistics in fully developed channel flow at low Reynolds number

A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of 3300, based on the mean centerline velocity and channel half-width, with about 4 million grid points. All essential turbulence scales are resolved on the computational grid and no subgrid model is used. A large number of turbulence statistics are computed and compared with the existing experimental data at comparable Reynolds numbers. Agreements as well as discrepancies are discussed in detail. Particular attention is given to the behavior of turbulence correlations near the wall. A number of statistical correlations which are complementary to the existing experimental data are reported for the first time.

Kim, John↗
NASA NTRS1987

Local intermodal energy transfer of the secondary instability in a plane channel

A mathematical technique for analyzing local energy-transfer rates among wave-vector triads is developed and applied to the data generated by direct numerical simulations of waves with various types of initial conditions. Starting the simulation with a primary two-dimensional wave and random noise produced structures very similar to those which evolved from three-dimensional center modes. The local transfer rates determined in this case help to explain the eventual deformation of the primary two-dimensional wave which was observed in wind-tunnel experiments. Including weak streamwise vortices in the initial flow field results in large amounts of energy being transferred to the K-type modes early in the simulation. The later development of the waves (and hence the energy-transfer rates) is similar to the previous case.

Singer, Bart A.↗
NASA NTRS1988

Annual Research Briefs, 1987

Lagrangian techniques have found widespread application to the prediction and understanding of turbulent transport phenomena and have yielded satisfactory results for different cases of shear flow problems. However, it must be kept in mind that in most experiments what is really available are Eulerian statistics, and it is far from obvious how to extract from them the information relevant to the Lagrangian behavior of the flow; in consequence, Lagrangian models still include some hypothesis for which no adequate supporting evidence was until now available. Direct numerical simulation of turbulence offers a new way to obtain Lagrangian statistics and so verify the validity of the current predictive models and the accuracy of their results. After the pioneering work of Riley (Riley and Patterson, 1974) in the 70's, some such results have just appeared in the literature (Lee et al, Yeung and Pope). The present contribution follows in part similar lines, but focuses on two particle statistics and comparison with existing models.

Moin, Parviz↗
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

Premixed turbulent flame calculation

The importance of turbulent premixed flames in a variety of applications has led to a substantial amount of effort towards improving the understanding of these flames. Although these efforts have increased the understanding, many questions still remain. The use of direct numerical simulation (DNS) in solving these questions is examined.

El-Tahry, S.↗
NASA NTRS1987

The decay of isotropic turbulence in a rapidly rotating frame

A direct numerical simulation of the decay of initially isotropic turbulence in a rapidly rotating frame was conducted. This 128 x 128 x 128 simulation was completed for a Reynolds number Re sub lambda = 15.3 and a Rossby number Ro sub lambda = 0.07 based on the initial turbulent kinetic energy and Taylor microscale. The numerical results indicate that the turbulence remains essentially isotropic during the major part of the decay (i.e., beyond the point where the turbulent kinetic energy has decayed to less than 10 percent of its initial value). The rapid rotation has the primary effect of shutting off the energy transfer so that the turbulence dissipation (and hence the rate of decay of the turbulent kinetic energy) is substantially reduced. Consequently, the anisotropy tensor remains essentially unchanged while the energy spectrum undergoes a nearly linear viscous decay (the same results that are predicted by Rapid Distortion Theory which is only formally valid for much shorter elapsed times. Surprisingly, no Taylor-Proudman reorganization of the flow to a two-dimensional state is observed. The implications that these results have on turbulence modeling are discussed briefly along with prospective future research.

Speziale, C. G.↗
NASA NTRS1987

Helicity fluctuations in incompressible turbulent flows

Results from direct numerical simulations of several homogeneous flows and fully developed turbulent channel flow indicate that the probability distribution function (pdf) of relative helicity density exhibits at most a 20 percent deviation from a flat distribution. Isotropic flows exhibit a slight helical nature but the presence of mean strain in homogeneous turbulence suppresses helical behavior. All the homogeneous turbulent flows studied show no correlation between relative helicity density and the dissipation of turbulent kinetic energy. The channel flow simulations indicate that, except for low-dissipation regions near the outer edge of the buffer layer, there is no tendency for the flow to be helical. The strong peaks in the relative helicity density pdf and the association of these peaks with regions of low dissipation found in previous simulations by Pelz et al.(1985) are not observed.

Rogers, Michael M.↗
NASA NTRS1988

Mixing control in a plane shear layer

An investigation of mixing processes in a plane shear layer by a direct numerical simulation with a high-order numerical technique, the isoparametric spectral element method, is presented. The computational domain includes a region with the splitter plate that allows to start the flow as a laminar (Blasius) boundary layers with different velocities and concentrations, merging in the shear layer and undergoing transition to highly-unsteady structures with intensive mixing. The augmentation of mixing as a function of a downstream location, the intensity and type of vortical structures in the layer is investigated by an introduction of disturbance at the high-velocity inflow. The character of the dynamic mixing is displayed along with time averaged profiles and statistical characteristics.

Korczak, K. Z.↗
NASA NTRS1988

Growth of correlation in compressible two-dimensional magnetofluid turbulence

Spectral transfer has been proposed as the primary mechanism for generating outward-propagating Alfven waves in the solar wind. This process has been investigated extensively for imcompressible magnetofluids, but the issue of whether it occurs in compressible magnetofluids such as the solar wind remains unresolved. The results of direct numerical simulations of nonisentropic-compressible two-dimensional MHD turbulence indicate that, for systems with finite initial cross helicity, the correlation between the fluctuating velocity field and the fluctuating magnetic field grows as a function of time. This growth of correlation can be interpreted as a turbulent process, as shown by examination of modal wavenumber spectra.

Dahlburg, R. B.↗
NASA NTRS1987

Transport of passive scalars in a turbulent channel flow

A direct numerical simulation of a turbulent channel flow with three passive scalars at different molecular Prandtl numbers is performed. Computed statistics including the turbulent Prandtl numbers are compared with existing experimental data. The computed fields are also examined to investigate the spatial structure of the scalar fields. The scalar fields are highly correlated with the streamwise velocity; the correlation coefficient between the temperature and the streamwise velocity is as high as 0.95 in the wall region. The joint probability distributions between the temperature and velocity fluctuations are also examined; they suggest that it might be possible to model the scalar fluxes in the wall region in a manner similar to the Reynolds stresses.

Kim, John↗
NASA NTRS1988

Ejection mechanisms in the sublayer of a turbulent channel

The structure of the vorticity field in the viscous wall layer of a turbulent channel is studied by examining the results of a fully resolved direct numerical simulation. It is shown that this region is dominated by intense three-dimensional shear layers in which the dominant vorticity component is spanwise. The advection and reproduction processes of these structures are examined and shown to be consistent with the classical generation mechanism for two-dimensional Tollmien-Schlichting waves. This process is fundamentally different from the usually accepted mechanism involving hairpin vortices.

Jimenez, Javier↗
NASA NTRS1989

Low-dimensional chaos in turbulence

Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Vastano, John A.↗
NASA NTRS1989

The effects of particle loading on turbulence structure and modelling

The objective of the present research was to extend the Direct Numerical Simulation (DNS) approach to particle-laden turbulent flows using a simple model of particle/flow interaction. The program addressed the simplest type of flow, homogeneous, isotropic turbulence, and examined interactions between the particles and gas phase turbulence. The specific range of problems examined include those in which the particle is much smaller than the smallest length scales of the turbulence yet heavy enough to slip relative to the flow. The particle mass loading is large enough to have a significant impact on the turbulence, while the volume loading was small enough such that particle-particle interactions could be neglected. Therefore, these simulations are relevant to practical problems involving small, dense particles conveyed by turbulent gas flows at moderate loadings. A sample of the results illustrating modifications of the particle concentration field caused by the turbulence structure is presented and attenuation of turbulence by the particle cloud is also illustrated.

Squires, Kyle D.↗
NASA NTRS1989

Transition to turbulence in laminar hypersonic flow

Progress in a recently started project aimed at the prediction of transition to turbulence in hypersonic flow is briefly discussed. The prediction of transition to turbulence is a very important issue in the design of space vessels. Two space vehicles currently under investigation, namely the aeroassisted transfer vehicle (AOTV) and the trans-atmospheric vehicle (TAV), suffer from strong aerodynamic heating. This heating is strongly influenced by the boundary layer structure. These aerospace vehicles fly in the upper atmospheric layer at a Mach number between 10 and 30 at very low atmospheric pressures. At very high altitudes the flow is laminar, but when the space vessel returns to a lower orbit, the flow becomes turbulent and the heating is dramatically increased. The prediction of this transition process is commonly done by means of experiments. The experimental facilities available nowadays cannot model the hypersonic flow field accurately enough by limitations in Mach and Reynolds number. These facilities also have a large free stream disturbance level which makes it very difficult to investigate transition accurately. An alternative approach is to study transition by theoretical means. Up to now numerical studies of hypersonic flow only discussed steady laminar or turbulent flow. This theoretical approach is extended to the study of transition in hypersonic flow by means of direct numerical simulations and additional theoretical investigations to explain the mechanisms leading to transition. A brief outline of how this research is to be performed is given.

Vandervegt, Jaap↗