Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent flows”

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 973 records · Page 54

Simulation of turbulent shear flows at Stanford and NASA-Ames - What can we do and what have we learned?

The capabilities and limitations of large eddy simulation (LES) and full turbulence simulation (FTS) are outlined. It is pointed out that LES, although limited at the present time by the need for periodic boundary conditions, produces large-scale flow behavior in general agreement with experiments. What is more, FTS computations produce small-scale behavior that is consistent with available experiments. The importance of the development work being done on the National Aerodynamic Simulator is emphasized. Studies at present are limited to situations in which periodic boundary conditions can be applied on boundaries of the computational domain where the flow is turbulent.

Reynolds, W. C.↗

The structure of the vorticity field in turbulent channel flow. Part 2: Study of ensemble-averaged fields

Several conditional sampling techniques are applied to a data base generated by large-eddy simulation of turbulent channel flow. It is shown that the bursting process is associated with well-organized horseshoe vortices inclined at about 45 deg. to the wall. These vortical structures are identified by examining the vortex lines of three-dimensional, ensemble averaged vorticity fields. Two distinct horseshoe-shaped vortices corresponding to the sweep and ejection events are detected. These vortices are associated with high Reynolds shear stress and hence make a significant contribution to turbulent energy production. The dependency of the ensemble averaged vortical structures on the detection criteria, and the question of whether this ensemble-averaged structure is an artifact of the ensemble averaging process are examined. The ensemble-averaged pattern of these vortical structures that emerge from the analysis could provide the basis for a hypothetical model of the organized structures of wall-bounded shear flows.

Kim, J.↗

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.↗

Computation of turbulent supersonic flows around pointed bodies having crossflow separation

The numerical method developed by Schiff and Sturek (1980) on the basis of the thin-layer parabolized Navier-Stokes equations of Schiff and Steger (1980) is extended to the case of turbulent supersonic flows on pointed bodies at high angles of attack. The governing equations, the numerical scheme, and modifications to the algebraic eddy-viscosity turbulence model are described; and results for three cones and one ogive-cylinder body (obtained using grids of 50 nonuniformly spaced points in the radial direction between the body and the outer boundary) are presented graphically and compared with published experimental data. The grids employed are found to provide sufficient spatial resolution of the leeward-side vortices; when combined with the modified turbulence model, they are shown to permit accurate treatment of flows with large regions of crossflow separation.

Degani, D.↗

Developing flow in S-shaped ducts. 1: Square cross-section duct

Laser-Doppler velocimetry was used to measure the laminar and turbulent flow in an S-duct formed with two 22.5 deg sectors of a bend with ratio of mean radius of curvature to hydraulic diameter of 7.0. The boundary layers at the inlet to the bend were about 25% and 15% of the hydraulic diameter for the laminar and turbulent flows, respectively. Pressure-driven secondary flows develop in the first half of the S-duct and persist into the second half but are largely reversed by the exit plane as a consequence of the change in the sense of curvature. There is, however, a region near the outer wall of the second bend where the redistribution of the streamwise isotachs results in a reinforcement of the secondary flow which was established in the first half of the S-duct. The net redistribution of the streamwise isotachs is comparable to that occurring in unidirectional bends of stronger curvature. The wall pressure distribution was also measured for the turbulent flow and quantifies the expected large variations in the longitudinal pressure gradient distributions which occur at different radial locations.

Taylor, A. M. K. P.↗

Recent advances in PDF modeling of turbulent reacting flows

This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.

Leonard, Andrew D.↗

The effect of contraction on a homogeneous turbulent shear flow

In the present study, a homogeneous flow was generated and allowed to reach its asymptotic state. A two-dimensional contraction was inserted a short distance downstream. Both mean velocity and turbulence-stress profiles (normal and shear) were measured at several points along the flow direction. The pressure-strain terms were evaluated from the energy balance made along the wind-tunnel axis. On the basis of the obtained results, it is concluded that a uniformly sheared turbulent flow with approximate transverse homogeneity remains in its condition at all stages of its development through a two-dimensional contraction. Small deviations from homogeneity observed upstream of contraction become much less visible towards the end of contraction. The mean-flow development can be predicted adequately by the Bernoulli equation.

Sreenivasan, K. R.↗

Numerical calculations of turbulent reacting flow in a gas-turbine combustor

A numerical study for confined, axisymmetrical, turbulent diffusion flames is presented. Local mean gas properties are predicted by solving the appropriate conservation equations in the finite difference form with the corresponding boundary conditions. The k-epsilon two-equation turbulence model is employed to describe the turbulent nature of the flow. A two-step kinetic model is assumed to govern the reaction mechanism. The finite reaction rate is the smaller of an Arrhenius type of reaction rate and a modified version of eddy-breakup model. Reasonable agreement is observed between calculations and measurements, but to obtain better agreement, more work is needed on improvements of the above mathematical models. However, the present numerical study offers an improvement in the analysis and design of the gas turbine combustors.

Lin, Chin-Shun↗

A PDF closure model for compressible turbulent chemically reacting flows

The objective of the proposed research project was the analysis of single point closures based on probability density function (pdf) and characteristic functions and the development of a prediction method for the joint velocity-scalar pdf in turbulent reacting flows. Turbulent flows of boundary layer type and stagnation point flows with and without chemical reactions were be calculated as principal applications. Pdf methods for compressible reacting flows were developed and tested in comparison with available experimental data. The research work carried in this project was concentrated on the closure of pdf equations for incompressible and compressible turbulent flows with and without chemical reactions.

Kollmann, W.↗

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Asymmetry between deuterium and tritium turbulent particle flows

The asymmetry between deuterium (D) and tritium (T) turbulent particle fluxes in mixed D–T plasmas is studied with numerical simulations of nonlinear gyrokinetic turbulence in ion temperature gradient-dominated and trapped electron mode-dominated regimes. At 50–50 D–T concentration, the asymmetry, or flow separation, between D and T fluxes is such that the tritium is better confined than the deuterium in both regimes. To supplement the nonlinear simulations, an analytic quasilinear theory of the particle flux symmetry breaking is developed and is valid for general electron dynamics. Furthermore, this theory correctly predicts the numerically computed deviation in the ion density fraction (from 50–50), or ion density gradient (from the electron gradient), required to restore equal deuterium and tritium fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of Eddy Viscosity on Time Correlations in Large Eddy Simulation

Subgrid-scale (SGS) models for large. eddy simulation (LES) have generally been evaluated by their ability to predict single-time statistics of turbulent flows such as kinetic energy and Reynolds stresses. Recent application- of large eddy simulation to the evaluation of sound sources in turbulent flows, a problem in which time, correlations determine the frequency distribution of acoustic radiation, suggest that subgrid models should also be evaluated by their ability to predict time correlations in turbulent flows. This paper compares the two-point, two-time Eulerian velocity correlation evaluated from direct numerical simulation (DNS) with that evaluated from LES, using a spectral eddy viscosity, for isotropic homogeneous turbulence. It is found that the LES fields are too coherent, in the sense that their time correlations decay more slowly than the corresponding time. correlations in the DNS fields. This observation is confirmed by theoretical estimates of time correlations using the Taylor expansion technique. Tile reason for the slower decay is that the eddy viscosity does not include the random backscatter, which decorrelates fluid motion at large scales. An effective eddy viscosity associated with time correlations is formulated, to which the eddy viscosity associated with energy transfer is a leading order approximation.

He, Guowei↗

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.↗