Engineering Papers⌕ Search

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.

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

At least 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Numerical exploration of planetary arc dynamics

After analytically characterizing the response of a test particle to the combined effects of corotation and Lindblad resonances, due to the satellites of Neptune, as two coupled dynamical systems whose coupling depends on the collective effects in the arc, a formula is obtained for the energy imparted by the Lindblad resonance. This energy confines the arc material around the corotation resonance radius. The fact that the gradient of the torque density around the arc must be negative confirms the constraints defined by Lin et al. (1987) on the relative positions of the corotation and the Lindblad resonances if the arc is to be stable. These results are tested with a direct numerical simulation for two different arc configurations.

Sicardy, Bruno↗
NASA NTRS1991

Renormalization group analysis of anisotropic diffusion in turbulent shear flows

The renormalization group is applied to compute anisotropic corrections to the scalar eddy diffusivity representation of turbulent diffusion of a passive scalar. The corrections are linear in the mean velocity gradients. All model constants are computed theoretically. A form of the theory valid at arbitrary Reynolds number is derived. The theory applies only when convection of the velocity-scalar correlation can be neglected. A ratio of diffusivity components, found experimentally to have a nearly constant value in a variety of shear flows, is computed theoretically for flows in a certain state of equilibrium. The theoretical value is well within the fairly narrow range of experimentally observed values. Theoretical predictions of this diffusivity ratio are also compared with data from experiments and direct numerical simulations of homogeneous shear flows with constant velocity and scalar gradients.

Rubinstein, Robert↗
NASA NTRS1990

Nonlinear development of crossflow vortices

Nonlinear crossflow vortices in an incompressible three-dimensional boundary layer are computed by weakly nonlinear theory and direct numerical simulations. The parallel basic flow is defined by Falkner-Skan-Cooke similarity profiles. The temporal evolution of spanwise periodic, quasi-two-dimensional disturbances without variations along the vortex axis is considered. The nonlinear theory is based on the approach of Herbert (1980, 1983). The theory predicts the existence of the finite amplitude equilibrium states seen in earlier simulations. When the disturbance amplitudes are small, there is very good quantitative agreement in the fundamental disturbance velocity components between the theory and the simulations.

Singer, Bart A.↗
NASA NTRS1990

Instabilities in rotating channel flow

A combination of linear, weakly nonlinear, and secondary stability theories has been used to study the transition process in rotating channel flow. Results for parts of the theory have been confirmed by direct numerical simulations. Preliminary results indicate that the theory can help explain some experimental observations regarding the evolution of streamwise vortices.

Ng, Lian↗
NASA NTRS1991

Numerical investigation of instability and transition in rotating plane Poiseuille flow

Direct numerical simulations with initial small-amplitude random disturbances are used to investigate instability and transition in plane Pouseuille flow with spanwise system rotation, employing a large computational domain to facilitate the selection of the dominant spanwise wavelength of the observed streamwise vortices as well as the streamwise wavelength of the dominant secondary modes due to streamwise vortices. In all cases studied, the unstable primary three-dimensional waves persist in the nonlinear stage. Where there is primary three-dimensional wave decay, the finite-amplitude streamwise vortices can become unstable to two different three-dimensional secondary instability modes, depending on the Reynolds and the rotation numbers.

Yang, Kyung-Soo↗
NASA NTRS1990

Similarity of organized structures in turbulent shear flows

Recent analysis of databases generated by direct numerical simulations of homogeneous turbulent shear flows have revealed the presence of coherent structures similar to those in turbulent boundary layers. In this paper these findings and tentative conclusions on their significance are discussed.

Moin, Parviz↗
NASA NTRS1991

Dynamical eigenfunction decomposition of turbulent channel flow

The results of an analysis of low-Reynolds-number turbulent channel flow based on the Karhunen-Loeve (K-L) expansion are presented. The turbulent flow field is generated by a direct numerical simulation of the Navier-Stokes equations at a Reynolds number Re(tau) = 80 (based on the wall shear velocity and channel half-width). The K-L procedure is then applied to determine the eigenvalues and eigenfunctions for this flow. The random coefficients of the K-L expansion are subsequently found by projecting the numerical flow field onto these eigenfunctions. The resulting expansion captures 90 percent of the turbulent energy with significantly fewer modes than the original trigonometric expansion. The eigenfunctions, which appear either as rolls or shearing motions, possess viscous boundary layers at the walls and are much richer in harmonics than the original basis functions.

Ball, K. S.↗
NASA NTRS1991

Similarity between turbulent kinetic energy and temperature spectra in the near-wall region

The similarity between turbulent kinetic energy and temperature spectra, previously confirmed using experimental data in various turbulent shear flows, is validated in the near-wall region using direct numerical simulation data in a fully developed turbulent channel flow. The dependence of this similarity on the molecular Prandtl number is also examined.

Antonia, R. A.↗
NASA NTRS1991

Mixing transition and the cascade to small scales in a plane mixing layer

Direct numerical simulations of time-developing plane mixing layers have been performed using a variety of initial conditions. Up to two pairings of the dominant spanwise rollers have been simulated. When the flow is sufficiently three-dimensional, a pairing can cause the mixing layer to undergo a transition to small-scale turbulence. This small-scale transition is accompanied by increased mixing of a passive scalar, similar to observations of the mixing transition in experiments. As part of the transition process, thin vortex sheets are generated by vortex stretching and roll up in a two-dimensional mixing layer. This higher-order rollup is part of the cascade to small-scale turbulence. Estimates are obtained for the degree of three-dimensionality required for the pairing to initiate the transition.

Moser, R. D.↗
NASA NTRS1991

Preferential concentration of particles by turbulence

Direct numerical simulation of isotropic turbulence was used to investigate the effect of turbulence on the concentration fields of heavy particles. The hydrodynamic field was computed using 64-cubed points, and a statistically stationary flow was obtained by forcing the low-wave-number components of the velocity field. The particles used in the simulations were time advanced according to Stokes drag law and were also assumed to be much more dense than the fluid. Properties of the particle cloud were obtained by following the trajectories of 1,000,000 particles through the simulated flow fields. Three values of the ratio of the particle time constant to large-scale turbulence time scale were used in the simulations: 0.075, 0.15, and 0.52. The simulations show that the particles collect preferentially in regions of low vorticity and high strain rate. This preferential collection was most pronounced for particles with the instantaneous number density as much as 25 times the mean value. The fact that dense particles collect in regions of low vorticity and high strain in turn implies that turbulence may actually inhibit rather than enhance mixing of particles.

Squires, Kyle D.↗
NASA NTRS1991

Scalar dissipation and mixing in turbulent reacting flows

The mixing term that appears in probability density function-based formulations of the statistical behavior of turbulent reacting flows was studied with the use of direct numerical simulations of an irreversible, second-order, isothermal chemical reaction between initially segregated reactants in decaying, homogeneous turbulence. Three-dimensional, time-dependent pseudospectral calculations with 64-cubed Fourier modes were made as well as were Monte Carlo simulations of the joint concentration pdf equations for coalescence-dispersion (C/D) models. It was found that the development of the scalar microscales depends largely on the initial concentration distributions (both isotropic and nonisotropic distributions were used) and is relatively independent of the reaction rate constant. Joint concentration probability density functions were compared to predictions of D/D models that make use of parameters evaluated from the nonreacting or inert system.

Leonard, Andy D.↗
NASA NTRS1990

Numerical investigation of the spreading of self-excited stratified jets

The structure and evolution of self-excited subsonic periodic arrays of jets of constant and variable density are studied using spectral-element direct numerical simulations. The governing equation of motion is presented, and a method based on spectral element discretizations appropriate for simulating arbitrarily complex geometry jets and large density variations for subsonic flows is developed. Variable density fields are found to be more unstable than the corresponding uniform density fields with much higher rms values; as a result, their spreading is also considerably larger. There is a dramatic increase in spreading after a few pairings occur. Findings presented for low and high side-momentum flux reveal a shifting of the origin of instability from the near-field to the far-field, respectively, and suggest possible routes of stabilization.

Batcho, P. F.↗
NASA NTRS1991

Asymptotic behavior of curvature of surface elements in isotropic turbulence

The asymptotic behavior of the curvature of material elements in turbulence is investigated using Lagrangian velocity-gradient time series obtained from direct numerical simulations of isotropic turbulence. Several material-element ensembles of different initial curvatures and shapes are studied. It is found that, at long times, the (first five) moments of the logarithm of characteristic curvature and shape factor asymptote to values that are independent of the initial curvature or shape. This evidence strongly suggests that the asymptotic pdf's of the curvature and shape of material elements are stationary and independent of initial conditions. Irrespective of initial curvature or shape, the asymptotic shape of a material surface is cylindrical with a high probability.

Girimaji, S. S.↗
NASA NTRS1989

The topology and vorticity dynamics of a three-dimensional plane compressible wake

The three-dimensional aspects of transition in a low Mach number plane compressible wake are studied numerically. Comparisons are made between the topology of the velocity field and the vorticity dynamics of the flow based on results from direct numerical simulations of the full compressible Navier-Stokes equations. The velocity field is integrated to obtain instantaneous streamlines at different stages in the evolution. A generalized three-dimensional critical point theory is applied to classify the critical points of the velocity field.

Chen, Jacqueline H.↗
NASA NTRS1991

Motion of particles with inertia in a compressible free shear layer

The effects of the inertia of a particle on its flow-tracking accuracy and particle dispersion are studied using direct numerical simulations of 2D compressible free shear layers in convective Mach number (Mc) range of 0.2 to 0.6. The results show that particle response is well characterized by tau, the ratio of particle response time to the flow time scales (Stokes' number). The slip between particle and fluid imposes a fundamental limit on the accuracy of optical measurements such as LDV and PIV. The error is found to grow like tau up to tau = 1 and taper off at higher tau. For tau = 0.2 the error is about 2 percent. In the flow visualizations based on Mie scattering, particles with tau more than 0.05 are found to grossly misrepresent the flow features. These errors are quantified by calculating the dispersion of particles relative to the fluid. Overall, the effect of compressibility does not seem to be significant on the motion of particles in the range of Mc considered here.

Samimy, M.↗
NASA NTRS1991

Reynolds shear stress and heat flux calculations in a fully developed turbulent duct flow

The use of a modified form of the Van Driest mixing length for a fully developed turbulent channel flow leads to mean velocity and Reynolds stress distributions that are in close agreement with data obtained either from experiments or direct numerical simulations. The calculations are then extended to a nonisothermal flow by assuming a constant turbulent Prandtl number, the value of which depends on the molecular Prandtl number. Calculated distributions of mean temperature and lateral heat flux are in reasonable agreement with the simulations. The extension of the calculations to higher Reynolds numbers provides some idea of the Reynolds number required for scaling on wall variables to apply in the inner region of the flow.

Antonia, R. A.↗
NASA NTRS1991

On turbulent spots in plane Poiseuille flow

Turbulence characteristics inside a turbulent spot in plane Poiseuille flow are investigated by analyzing a database 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 r.m.s. 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 NTRS1990

Enstrophy budget in decaying compressible turbulence

The interaction of vorticity and shock waves in decaying compressible turbulence is investigated by direct numerical simulation of the Navier-Stokes equation using a spectral method with 64-cubed collocation points. The Taylor microscale Reynolds number decreases from 36 to 17 and the rms Mach number from 0.8 to 0.4 during the decay. It is observed that vorticity is created inside curved shocks through baroclinic interaction and intensified by the compression of fluid elements. Vortex stretching is, on the other hand, dominant outside shocks and is responsible for enhancing the strength of vorticity. A dipole of vorticity is typically created at the vertex of curved shocks, the structure of which is explained in terms of jump relations of velocity derivatives across a curved shock.

Kida, Shigeo↗