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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Direct numerical simulation of curved turbulent channel flow

Low Reynolds number, mildly curved, turbulent channel flow has been simulated numerically without subgrid scale models. A new spectral numerical method developed for this problem was used, and the computations were performed with 2 million degrees of freedom. A variety of statistical and structural information has been extracted from the computed flow fields. These include mean velocity, turbulence stresses, velocity skewness, and flatness factors, space time correlations and spectra, all the terms in the Reynolds stress balance equations, and contour and vector plots of instantaneous velocity fields. The effects of curvature on this flow were determined by comparing the concave and convex sides of the channel. The observed effects are consistent with experimental observations for mild curvature. The most significant difference in the turbulence statistics between the concave and convex sides was in the Reynolds shear stress. This was accompanied by significant differences in the terms of the Reynolds shear stress balance equations. In addition, it was found that stationary Taylor-Gortler vortices were present and that they had a significant effect on the flow by contributing to the mean Reynolds shear stress, and by affecting the underlying turbulence.

Moser, R. D.↗
NASA NTRS1986

Direct numerical simulations of a temporally evolving mixing layer subject to forcing

The vortical evolution of mixing layers subject to various types of forcing is numerically simulated using pseudospectral methods. The effect of harmonic forcing and random noise in the initial conditions is examined with some results compared to experimental data. Spanwise forcing is found to enhance streamwise vorticity in a nonlinear process leading to a slow, secondary growth of the shear layer. The effect of forcing on a chemical reaction is favorably compared with experimental data at low Reynolds numbers. Combining harmonic and subharmonic forcing is shown to both augment and later destroy streamwise vorticity.

Claus, Russell W.↗
NASA NTRS1987

Direct numerical simulation of buoyantly driven turbulence

Numerical simulations of homogeneous turbulence subject to buoyant forcing were performed. The presence of a mean temperature gradient combined with a gravitational field results in a forcing term in the momentum equations. The development of the turbulence was studied and compared to the decay of similar fields in the absence of gravity. In the buoyantly driven field, the vorticity is preferentially aligned with the intermediate eigenvector of the strain-rate tensor and the local temperature gradient is more likely to be aligned with the most compressive eigenvector. These relationships are qualitatively similar to those observed in previous shear flow results studied by Ashurst (1987). A tensor diffusivity model for passive scalar transport developed from shear flow results in Rogers, Moin, and Reynolds (1986) also predicts this buoyant scalar transport, indicating that the relationship between the scalar flux and the Reynolds stress is similar in both flows.

Ashurst, William T.↗
NASA NTRS1989

Direct numerical simulations of turbulent convection with a variable gravity and Keplerian rotation

Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Cabot, William H.↗
NASA NTRS1989

Direct numerical simulations of a temporally evolving mixing layer subject to forcing

The vortical evolution of mixing layers subject to various types of forcing is numerically simulated using pseudospectral methods. The effect of harmonic forcing and random noise in the initial conditions is examined with some results compared to experimental data. Spanwise forcing is found to enhance streamwise vorticity in a nonlinear process leading to a slow, secondary growth of the shear layer. The effect of forcing on a chemical reaction is favorably compared with experimental data at low Reynolds numbers. Combining harmonic and subharmonic forcing is shown to both augment and later destroy streamwise vorticity.

Claus, R. W.↗
NASA NTRS1989

The use of direct numerical simulation in the study of turbulent, chemically-reacting flows

Full turbulence simulations are used here to study the effects of chemical heat release on the large-scale structures in turbulent mixing layers. In agreement with laboratory results, it is found that the heat release lowers the rate at which the mixing layer grows and reduces the rate at which chemical products are formed. The baroclinic torque and thermal expansion in the mixing layer 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 the large-scale structures. Previously unexplained anomalies observed in the mean velocity profiles of reacting jets and mixing layer are shown to result from vorticity generation by baroclinic torques. Calculations of the energy in the various wavenumbers shows that the heat release has a stabilizing effect on the growth rate of individual modes. This methodology can be applied for Reynolds numbers less than several hundred and for Damkoehler numbers less than about ten.

Riley, J. J.↗
NASA NTRS1990

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗
NASA NTRS1992

Direct numerical simulation of a combusting droplet with convection

The evaporation and combustion of a single droplet under forced and natural convection was studied numerically from first principles using a numerical scheme that solves the time-dependent multiphase and multispecies Navier-Stokes equations and tracks the sharp gas-liquid interface cutting across an arbitrary Eulerian grid. The flow fields both inside and outside of the droplet are resolved in a unified fashion. Additional governing equations model the interphase mass, energy, and momentum exchange. Test cases involving iso-octane, n-hexane, and n-propanol droplets show reasonable comparison rate, and flame stand-off distance. The partially validated code is, thus, readied to be applied to more demanding droplet combustion situations where substantial drop deformation render classical models inadequate.

Liang, Pak-Yan↗
NASA NTRS1992

Direct numerical simulation of laminar breakdown in high-speed, axisymmetric boundary layers

The compressible Navier-Stokes equations are solved using spectral collocation and high-order compact-difference techniques to simulate the laminar breakdown in high-speed, axisymmetric boundary-layer flow. Mach 4.5 flow along a hollow cylinder and Mach 6.8 flow along a sharp cone are considered. Data obtained replicate two previously unexplained phenomena, namely, the appearance of so-called 'rope-like waves' and 'the precursor transition effect', in which transitional flow originates near the critical layer well upstream of the transition location at the wall. The numerical data also reveal that neither of these effects can be explained, even qualitatively, by linear stability theory alone. It is shown that rope-like appearance arises from secondary instability. Certain features of the precursor transition effect also emerge from secondary instability but its nature is revealed to be fundamentally nonlinear.

Pruett, C. D.↗
NASA NTRS1992

Direct numerical simulations of a reacting turbulent mixing layer by a pseudospectral-spectral element method

An account is given of the implementation of the spectral-element technique for simulating a chemically reacting, spatially developing turbulent mixing layer. Attention is given to experimental and numerical studies that have investigated the development, evolution, and mixing characteristics of shear flows. A mathematical formulation is presented of the physical configuration of the spatially developing reacting mixing layer, in conjunction with a detailed representation of the spectral-element method's application to the numerical simulation of mixing layers. Results from 2D and 3D calculations of chemically reacting mixing layers are given.

Mcmurtry, Patrick A.↗
NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗