Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence 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 307 records · Page 17

Large-eddy simulation of a turbulent mixing layer

The three dimensional, time dependent (incompressible) vorticity equations were used to simulate numerically the decay of isotropic box turbulence and time developing mixing layers. The vorticity equations were spatially filtered to define the large scale turbulence field, and the subgrid scale turbulence was modeled. A general method was developed to show numerical conservation of momentum, vorticity, and energy. The terms that arise from filtering the equations were treated (for both periodic boundary conditions and no stress boundary conditions) in a fast and accurate way by using fast Fourier transforms. Use of vorticity as the principal variable is shown to produce results equivalent to those obtained by use of the primitive variable equations.

Mansour, N. N.↗

Simulation of electrostatic turbulence in the plasma sheet boundary layer with electron currents and bean-shaped ion beams

Plasma data from ISEE-1 show the presence of electron currents as well as energetic ion beams in the plasma sheet boundary layer. Broadband electrostatic noise and low-frequency electromagnetic bursts are detected in the plasma sheet boundary layer, especially in the presence of strong ion flows, currents, and steep spacial gradients in the fluxes of few-keV electrons and ions. Particle simulations have been performed to investigate electrostatic turbulence driven by a cold electron beam and/or ion beams with a bean-shaped velocity distribution. The simulation results show that the counterstreaming ion beams as well as the counterstreaming of the cold electron beam and the ion beam excite ion acoustic waves with a given Doppler-shifted real frequency. However, the effect of the bean-shaped ion velocity distributions reduces the growth rates of ion acoustic instability. The simulation results also show that the slowing down of the ion bean is larger at the larger perpendicular velocity. The wave spectra of the electric fields at some points of the simulations show turbulence generated by growing waves.

Nishikawa, K.-I.↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1996 Summer Program - Part 6

Topics considered include: New approach to turbulence modeling; Second moment closure analysis of the backstep flow database; Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers; Turbulent flame propagation in partially premixed flames; Ensemble averaged dynamic modeling. Also included a study of the turbulence structures of wall-bounded shear flows; Simulation and modeling of the elliptic streamline flow.

Source record↗

Numerical simulation of MHD turbulence in three dimensions

The evolution of Alfvenic turbulence in 3D spherical geometry can now be studied. In simulations, a fast stream is sandwiched between two slower streams. The inflow is both supersonic and superAlfvenic. Alfven waves entering the box are convected into the medium and interact nonlinearly with the velocity shear and with any structures (i.e., flux tubes) that might be present. These initial simulations suggest that velocity shear, even in spherical geometry, is able to drive a turbulent cascade which results in approximately Kolmogoroff-like power spectra.

Goldstein, M. L.↗

Advanced in turbulence physics and modeling by direct numerical simulations

The advent of direct numerical simulations of turbulence has opened avenues for research on turbulence physics and turbulence modeling. Direct numerical simulation provides values for anything that the scientist or modeler would like to know about the flow. An overview of some recent advances in the physical understanding of turbulence and in turbulence modeling obtained through such simulations is presented.

Reynolds, W. C.↗

Simulation of electrostatic turbulence due to sheared flows parallel and transverse to the magnetic field

A spatially two-dimensional electrostatic particle simulation code is used to examine the stability of a plasma equilibrium characterized by a localized transverse dc electric field and a magnetic-field-aligned electron drift for L much less than Lx, where Lx is the simulation length in the x direction and L is the scale length associated with the dc electric field. It is found that the dc electric field and the field-aligned current can together play a synergistic role to enable the excitation of electrostatic waves even when the threshold values of the field-aligned drift and the E x B drift are individually subcritical. The simulation results indicate that a broadband turbulence is associated with such an equilibrium.

Nishikawa, K.-I.↗

Critical comparison of second-order closures with direct numerical simulations of homogeneous turbulence

Recently, several models have been proposed for closing the second-moment equations, in which the velocity-pressure gradient tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so-called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behaviors in different flows have not been fully examined because they are unmeasurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the possibility to do this kind of study. With direct numerical simulation, we may use the solution to calculate exactly the values of these correlation terms and then directly compare them with the values from their modeled formulations. In this paper, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of 45 homogeneous flows, which were done by Rogers et al. (1986) and Lee and Reynolds (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones that give the best performance. The paper also describes the modeling procedure, model constraints, and the various evaluated models. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

Direct simulations of premixed turbulent flames with nonunity Lewis numbers

A principal effect of turbulence on premixed flames in the flamelet regime is to wrinkle the flame fronts. For nonunity Lewis numbers, Le is not equal to 1, the local flame structure is altered in curved regions. This effect is examined using direct numerical simulations of 3D isotropic turbulence with constant density, single-step Arrhenius kinetics chemistry. Simulations of Lewis numbers 0.8, 1.0, and 1.2 are compared. At the local level, curvature effects dominated changes to the flame structure while strain effects were insignificant. A strong Lewis-number-dependent correlation was found between surface curvature and the local flame speed. The correlation was positive for Le less than 1 and negative for Le greater than 1. At the global level, strain-related effects were more significant than curvature effects. The turbulent flame speed changed significantly with Lewis number, increasing as Le decreased. This was found to be due to strain effects that have a nonzero mean over the flame surface, rather than to curvature effects that have a nearly zero mean. The mean product temperature was also found to vary with Lewis number, being higher for Le greater than 1 and lower for Le less than 1.

Rutland, C. J.↗

A critical comparison of second order closures with direct numerical simulation of homogeneous turbulence

Recently, several second order closure models have been proposed for closing the second moment equations, in which the velocity-pressure gradient (and scalar-pressure gradient) tensor and the dissipation rate tensor are two of the most important terms. In the literature, these correlation tensors are usually decomposed into a so called rapid term and a return-to-isotropy term. Models of these terms have been used in global flow calculations together with other modeled terms. However, their individual behavior in different flows have not been fully examined because they are un-measurable in the laboratory. Recently, the development of direct numerical simulation (DNS) of turbulence has given us the opportunity to do this kind of study. With the direct numerical simulation, we may use the solution to exactly calculate the values of these correlation terms and then directly compare them with the values from their modeled formulations (models). Here, we make direct comparisons of five representative rapid models and eight return-to-isotropy models using the DNS data of forty five homogeneous flows which were done by Rogers et al. (1986) and Lee et al. (1985). The purpose of these direct comparisons is to explore the performance of these models in different flows and identify the ones which give the best performance. The modeling procedure, model constraints, and the various evaluated models are described. The detailed results of the direct comparisons are discussed, and a few concluding remarks on turbulence models are given.

Shih, Tsan-Hsing↗

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible, (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a non-dimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence.

Erlebacher, Gordon↗

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a nondimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence.

Erlebacher, Gordon↗

Tetrahedral-Mesh Simulations of Shock-Turbulence Interaction

Despite decades of development of unstructured mesh methods, direct numerical simulations (DNS) of turbulent flows are still predominantly performed on structured or unstructured hexahedral meshes with high-order finite-difference methods, weighted essentially nonoscillatory (WENO) schemes, or hybrid schemes formed by their combinations. Tetrahedral meshes offer easy mesh generation and adaptation around complex geometries and the potential of an orientation-free grid that would benefit the isotropic nature of small-scale dissipation, as well as the solution accuracy of intermediate scales. To advance the state of the art of unstructured-mesh simulation capabilities for shock/turbulence interaction, DNS using pure tetrahedral meshes are carried out with the space-time conservation element, solution element (CESE) method in this research. By its design, the CESE method is constructed based on a non-dissipative scheme and is a genuinely multidimensional numerical framework that is free from the use of an approximate Riemann-solver. The numerical framework also provides the ability to add numerical dissipation (the nondissipative scheme acts as the reference state like that of the reversible state in thermodynamics) when needed (with justification from mathematics/physics). The above-mentioned features along with the CESE method's consistent shock-capturing approach and strong enforcement of flux conservation in spacetime offers a novel method to accurately simulate turbulent flows and their interaction with shocks using tetrahedral meshes. Two canonical problems, namely, isotropic turbulence interaction with a normal shock and a Mach 2.9 turbulent boundary layer flow over a 24deg compression corner are investigated in this study. Computational results show reasonably good agreement with experimental data and results from structured-mesh, high-order simulations available in the literature. Successful validation of these canonical problems demonstrated here paves the way for future high-fidelity supersonic flow simulations involving complex-geometries.

Venkatachari, Balaji Shankar↗

A Piloted Simulation Study of Wake Turbulence on Final Approach

A piloted simulation study has been conducted in a research simulator to provide a means to estimate the effects of different levels of wake turbulence on final approach. A worst-case methodology was used to ensure conservative estimates. Fourteen airline pilots voluntarily participated in the study and flew almost 1000 approaches. The pilots rated the subjective severity of the disturbances using a special rating scale developed for this study. Several objective measures of the airplane/pilot response to the simulated wake turbulence were also made. All the data showed a large amount of variation between pilots and to a lesser extent for a given pilot. Therefore, the data were presented at 50, 70, 90 percentile levels as a function of vortex strength. The data allow estimates of the vortex strength for a given subjective or objective response and vice versa. The results of this study appear to be more conservative than the results of previous studies.

Stewart, Eric C.↗

Direct Numerical Simulation of Hypersonic Turbulent Boundary Layer inside an Axisymmetric Nozzle

As a first step toward a study of acoustic disturbance field within a conventional, hypersonic wind tunnel, direct numerical simulations (DNS) of a Mach 6 turbulent boundary layer on the inner wall of a straight axisymmetric nozzle are conducted and the results are compared with those for a flat plate. The DNS results for a nozzle radius to boundary-layer thickness ratio of 5:5 show that the turbulence statistics of the nozzle-wall boundary layer are nearly unaffected by the transverse curvature of the nozzle wall. Before the acoustic waves emanating from different parts of the nozzle surface can interfere with each other and undergo reflections from adjacent portions of the nozzle surface, the rms pressure fluctuation beyond the boundary layer edge increases toward the nozzle axis, apparently due to a focusing effect inside the axisymmetric configuration. Spectral analysis of pressure fluctuations at both the wall and the freestream indicates a similar distribution of energy content for both the nozzle and the flat plate, with the peak of the premultiplied frequency spectrum at a frequency of [(omega)(delta)]/U(sub infinity) approximately 6.0 inside the free stream and at [(omega)(delta)]/U(sub infinity) approximately 2.0 along the wall. The present results provide the basis for follow-on simulations involving reverberation effects inside the nozzle.

Huang, Junji↗

Simulation of the turbulent Rayleigh-Benard problem using a spectral/finite difference technique

The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.

Eidson, T. M.↗

Simulation of flat-plate turbulent boundary layers in cryogenic tunnels

The magnitudes of real-gas effects on flat-plate turbulent boundary layer simulations in a cryogenic nitrogen wind tunnel are investigated in order to determine the validity of the method used by Inger (1979) to estimate real-gas effects. Boundary layer solutions for real gases, ideal gases with a specific heat ratio of 1.6 and ideal diatomic gases (specific heat ratio 1.4) were obtained for the worst case conditions of maximum stagnation pressure (9 atm), minimum stagnation temperature (120 K) and Mach number of 1.2. Calculated boundary layer parameters such as friction coefficient and displacement thickness are shown to agree closely for the real gas and the ideal diatomic gas (specific heat ratio 1.4), while the ideal gas solution used by Inger is shown to differ from the real-gas values considerably. Results indicate that real-gas effects on a flat-plate turbulent boundary layer simulation in a cryogenic nitrogen tunnel are insignificant, and suggest the unlikelihood of the large real-gas effects reported by Inger for turbulent boundary layer shock interactions.

Adcock, J.↗

Efforts to Simulate Solar Wind Turbulence

A three-dimensional integration of the MHD equations in spherical coordinates has been developed that attempts to simulate a variety of solar wind conditions. These include the interaction of Alfven wave packets and the development of a turbulent cascade, the role of the heliospheric current sheet, the role of quasi-two-dimensional fluctuations in determining how magnetic field lines meander throughout the heliosphere, and the role of interstellar pickup ions in perturbing the solar wind in the outer heliosphere.

Goldstein, Melvyn L.↗

Studying Turbulence Using Numerical Simulation Databases - IX: Proceedings of the 2002 Summer Program

The ninth Summer Program of the Center for Turbulence Research was held during the period July 29th - August 23rd, 2002. The increase in number of participants, noted in the Preface to the Proceedings of the 2000 Program, continues: this year there were 50 participants from ten countries, and 30 hosts from Stanford and NASA-Ames. This Proceedings volume contains 32 papers that span a wide range of topics and an enormous range of physical scales. The papers have been divided into seven groups: Acoustics, RANS modeling, Combustion, Large-eddy simulation (LES), LES Numerics, Stratified Flows, and Fundamentals, In several cases, a paper could have fitted in more than one group so the classification is somewhat arbitrary.

Bradshaw, Peter↗