Numerical Simulations of the Late Stages of Transition to Turbulence
This viewgraph presentation reviews direct numerical simulation in the late stages of the transition process to turbulence.
Sandham, N. D.↗
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.
Assessment of Turbulence-Chemistry Interaction Models in the National Combustion Code (NCC) - Part I
This paper describes the implementations of the linear-eddy model (LEM) and an Eulerian FDF/PDF model in the National Combustion Code (NCC) for the simulation of turbulent combustion. The impacts of these two models, along with the so called laminar chemistry model, are then illustrated via the preliminary results from two combustion systems: a nine-element gas fueled combustor and a single-element liquid fueled combustor.
Wey, Thomas Changju↗
The Effect of Two-Dimensional Forcing on Self-Similar Turbulent Plane Wakes
Three direct numerical simulations of time-evolving turbulent plane wakes with velocity deficit Reynolds numbers of about 2,000 have been simulated using a spectral numerical method with up to 600 x 260 x 160 modes. The initial conditions for the simulations are generated from direct numerical simulations of a turbulent boundary layer (momentum thickness Reynolds number of 670), and varying amounts of additional two- dimensional, forcing. In order to preserve the self-similar flow evolution, the forcing is implemented by multiplying all the two-dimensional modes in the initial condition by a constant factor. In the "natural" case no additional forcing is used; in the "forced" and "heavily forced" cases this factor is 5 and 20, respectively. The wake spreading rate Is increased by factors of 1.7 and 7.1 for the two forced cases. The Reynolds stresses are also increased by a similar or even larger factor. These results indicate that the plane wake is much more sensitive to initial forcing than the plane mixing layer. As in the plane mixing layer, two-dimensional forcing promotes more organized large-scale vortical flow structures and these structures axe sometimes separated by "braid regions" containing streamwise "rib" vortices, unlike in the unforced wake.
Rogers, M. M.↗
Turbulence structures associated with the bursting event
Turbulence structures in a wall bounded shear layer during the bursting event detected by a conditional sampling technique are investigated using data obtained from large eddy simulation of turbulent channel flow. Streamlines are constructed from the ensemble-averaged velocity field to illustrate the flow patterns associated with the bursting event. The spatting motions during the sweep event and the existence of a pair of counterrotating streamwise vortices during the ejection process are exhibited.
Kim, J.↗
Turbulence structures associated with the bursting event
Turbulence structures in a wall bounded shear layer during the bursting event detected by a conditional sampling technique are investigated using data obtained from large eddy simulation of turbulent channel flow. Streamlines are constructed from the ensemble-averaged velocity field to illustrate the flow patterns associated with the bursting event. The splatting motions during the sweep event and the existence of a pair of counterrotating streamwise vortices during the ejection process are exhibited.
Kim, J.↗
A Kinematically Consistent Two-Point Correlation Function
A simple kinematically consistent expression for the longitudinal two-point correlation function related to both the integral length scale and the Taylor microscale is obtained. On the inner scale, in a region of width inversely proportional to the turbulent Reynolds number, the function has the appropriate curvature at the origin. The expression for two-point correlation is related to the nonlinear cascade rate, or dissipation epsilon, a quantity that is carried as part of a typical single-point turbulence closure simulation. Constructing an expression for the two-point correlation whose curvature at the origin is the Taylor microscale incorporates one of the fundamental quantities characterizing turbulence, epsilon, into a model for the two-point correlation function. The integral of the function also gives, as is required, an outer integral length scale of the turbulence independent of viscosity. The proposed expression is obtained by kinematic arguments; the intention is to produce a practically applicable expression in terms of simple elementary functions that allow an analytical evaluation, by asymptotic methods, of diverse functionals relevant to single-point turbulence closures. Using the expression devised an example of the asymptotic method by which functionals of the two-point correlation can be evaluated is given.
Ristorcelli, J. R.↗
Annual Research Briefs, 2002: Center for Turbulence Research
Turbulent combustion remains the largest component of the CTR's core program. This program and several related activities at CTR are supported by NASA's Ultra Efficient Engine Technology Program. It is also intimately connected with the Department of Energy's ASCI program at Stanford which develops the technology for numerical simulation of realistic aircraft engines using state of the art massively parallel computers. In combustion modeling the attention has been directed to the modeling of higher levels of complexity such as spray dynamics, radiation and soot formation. Major aircraft engine manufacturers have shown considerable interest in this program; in particular, a significant active collaboration exists between CTR and the Pratt & Whitney Corporation. CTR's combustion program is essentially based on the large-eddy simulation technique, LES, which is actively being pursued at CTR for this and many other applications. Important accomplishments in LES included simulations with three-dimensional filters which result in grid independent calculations (that is why we call it "true" LES), and the development of the methodology for integration of LES and Reynolds Averaged computations. Optimization techniques are being studied and used for the important problem of wall boundary conditions for LES as well as for optimal shape design for aeroacoustic and aerodynamic performance gains.
ULTRA EFFICIENT ENGINE TECHNOLOGY↗
Closed-Loop HIRF Experiments Performed on a Fault Tolerant Flight Control Computer
ABSTRACT Closed-loop HIRF experiments were performed on a fault tolerant flight control computer (FCC) at the NASA Langley Research Center. The FCC used in the experiments was a quad-redundant flight control computer executing B737 Autoland control laws. The FCC was placed in one of the mode-stirred reverberation chambers in the HIRF Laboratory and interfaced to a computer simulation of the B737 flight dynamics, engines, sensors, actuators, and atmosphere in the Closed-Loop Systems Laboratory. Disturbances to the aircraft associated with wind gusts and turbulence were simulated during tests. Electrical isolation between the FCC under test and the simulation computer was achieved via a fiber optic interface for the analog and discrete signals. Closed-loop operation of the FCC enabled flight dynamics and atmospheric disturbances affecting the aircraft to be represented during tests. Upset was induced in the FCC as a result of exposure to HIRF, and the effect of upset on the simulated flight of the aircraft was observed and recorded. This paper presents a description of these closed- loop HIRF experiments, upset data obtained from the FCC during these experiments, and closed-loop effects on the simulated flight of the aircraft.
Belcastro, Celeste M.↗
An experimental study of the effects of rapid rotation on turbulence
Experiments, large eddy simulations and direct numerical simulations all show that rapid rotation has a striking effect on homogeneous nearly isotropic turbulence. The cascade process is effectively inhibited by rotation, and thus dissipation is greatly reduced. Some attempts have been made to incorporate this effect in turbulence models. Numerical simulation further showed the somewhat surprising result that anisotropic homogeneous turbulence subjected to rotation tended towards an isotropic state; however, the residual anisotropy was not zero. Reynolds (1989) performed a detailed analysis using Rapid Distortion Theory (RDT) and showed that a reduction in the anisotropy is indeed expected and if the anisotropy is produced by irrotational strain then the anisotropy tensor b(sub ij) is asymptotically driven to half its initial value. Our objective is to extend the work of the experiments mentioned above to lower turbulent Rossby numbers and to confirm some of the results obtained by Reynolds (1989) for anisotropic turbulence.
Veeravalli, S. V.↗
Dual spectra and mixed energy cascade of turbulence in the wavelet representation
The wavelet-transformed Navier-Stokes equations are used to define quantities such as the transfer of kinetic energy and the flux of kinetic energy through scale r at position x. Direct numerical simulations of turbulent shear flow reveal that although their mean spatial values agree with their traditional counterparts in Fourier space, their spatial variability at every scale is very large, exhibiting non-Gaussian statistics. The local flux of energy involving scales smaller than some r also exhibits large spatial intermittency, and it is negative quite often, indicative of local inverse cascades.
Meneveau, Charles↗
Turbulent transport in premixed flames
Simulations of planar, premixed turbulent flames with heat release were used to study turbulent transport. Reynolds stress and Reynolds flux budgets were obtained and used to guide the investigation of important physical effects. Essentially all pressure terms in the transport equations were found to be significant. In the Reynolds flux equations, these terms are the major source of counter-gradient transport. Viscous and molecular terms were also found to be significant, with both dilatational and solenoidal terms contributing to the Reynolds stress dissipation. The BML theory of premixed turbulent combustion was critically examined in detail. The BML bimodal pdf was found to agree well with the DNS data. All BML decompositions, through the third moments, show very good agreement with the DNS results. Several BML models for conditional terms were checked using the DNS data and were found to require more extensive development.
Rutland, C. J.↗
Turbulent Concentration of MM-Size Particles in the Protoplanetary Nebula: Scaled-Dependent Multiplier Functions
The initial accretion of primitive bodies (asteroids and TNOs) from freely-floating nebula particles remains problematic. Here we focus on the asteroids where constituent particle (read "chondrule") sizes are observationally known; similar arguments will hold for TNOs, but the constituent particles in those regions will be smaller, or will be fluffy aggregates, and are unobserved. Traditional growth-bysticking models encounter a formidable "meter-size barrier" [1] (or even a mm-cm-size barrier [2]) in turbulent nebulae, while nonturbulent nebulae form large asteroids too quickly to explain long spreads in formation times, or the dearth of melted asteroids [3]. Even if growth by sticking could somehow breach the meter size barrier, other obstacles are encountered through the 1-10km size range [4]. Another clue regarding planetesimal formation is an apparent 100km diameter peak in the pre-depletion, pre-erosion mass distribution of asteroids [5]; scenarios leading directly from independent nebula particulates to this size, which avoid the problematic m-km size range, could be called "leapfrog" scenarios [6-8]. The leapfrog scenario we have studied in detail involves formation of dense clumps of aerodynamically selected, typically mm-size particles in turbulence, which can under certain conditions shrink inexorably on 100-1000 orbit timescales and form 10-100km diameter sandpile planetesimals. The typical sizes of planetesimals and the rate of their formation [7,8] are determined by a statistical model with properties inferred from large numerical simulations of turbulence [9]. Nebula turbulence can be described by its Reynolds number Re = L/eta sup(4/3), where L = ETA alpha sup (1/2) the largest eddy scale, H is the nebula gas vertical scale height, and α the nebula turbulent viscosity parameter, and η is the Kolmogorov or smallest scale in turbulence (typically about 1km), with eddy turnover time tη. In the nebula, Re is far larger than any numerical simulation can handle, so some physical model is needed to extend the results of numerical simulations to nebula conditions.
protoplanetary nebula↗
Studies in turbulence
Various papers on turbulence are presented. Individual topics addressed include: modeling the dissipation rate in rotating turbulent flows, mapping closures for turbulent mixing and reaction, understanding turbulence in vortex dynamics, models for the structure and dynamics of near-wall turbulence, complexity of turbulence near a wall, proper orthogonal decomposition, propagating structures in wall-bounded turbulence flows. Also discussed are: constitutive relation in compressible turbulence, compressible turbulence and shock waves, direct simulation of compressible turbulence in a shear flow, structural genesis in wall-bounded turbulence flows, vortex lattice structure of turbulent shear slows, etiology of shear layer vortices, trilinear coordinates in fluid mechanics.
Gatski, Thomas B.↗
Wall-Modeled Lattice Boltzmann and Navier-Stokes Approaches for Separated Flows
Lattice Boltzmann (LB) and hybrid Reynolds-averaged Navier-Stokes/large eddy simulation (RANS/LES) methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's Revolutionary Computational Aerosciences (RCA) standard test cases for separated flows. A detailed comparison between the performance and accuracy of the two emerging numerical methodologies for turbulence resolving simulations, i.e. the LB and hybrid RANS/LES methods will be presented. This contribution addresses the RCA technical challenge to identify and down-select critical turbulence, transition, and numerical method technologies for 40% reduction in predictive error for standard turbulence separated flow test cases. Results for the 2D NASA wall-mounted hump and the axisymmetric transonic bump including time-averaged pressure coefficient, skin friction, and velocity pro les, as well as resolved and modeled Reynolds stresses for both numerical approaches will be presented and differences between LB and hybrid RANS/LES will be discussed.
Lattice Boltzmann↗
Multi-Fidelity Heating Prediction of Adaptable, Deployable Entry Placement Technology Vehicles
The objective of this work was to investigate a multi-fidelity modeling approach to accurately and efficiently predict the laminar and turbulent convective heating on adaptable, deployable entry placement technology vehicles in Mars entry. A previously developed co-Kriging based multi-fidelity modeling approach was used to model the laminar and turbulent convective heat fluxes at several surface locations along the vehicle, including the rib sections. The laminar convective heat flux multi-fidelity model was found to have a mean convective heat rate error of approximately 3% when compared to high-fidelity CFD simulations. The turbulent convective heat flux multi-fidelity model was found to have a mean convective heat rate error of approximately 8% when compared to high-fidelity CFD simulations. Compared to a single fidelity model, the multi-fidelity model required approximately one-third the number of high-fidelity model evaluations to obtain the same accuracy level. The computational cost of evaluating the multi-fidelity model was approximately five orders of magnitude less than one high-fidelity model simulation.
Mario Santos↗
Mixing and chemical reaction in sheared and nonsheared homogeneous turbulence
Direct numerical simulations were made to examine the local structure of the reaction zone for a moderately fast reaction between unmixed species in decaying, homogeneous turbulence and in a homogeneous turbulent shear flow. Pseudospectral techniques were used in domains of 64 exp 3 and higher wavenumbers. A finite-rate, single step reaction between non-premixed reactants was considered, and in one case temperature-dependent Arrhenius kinetics was assumed. Locally intense reaction rates that tend to persist throughout the simulations occur in locations where the reactant concentration gradients are large and are amplified by the local rate of strain. The reaction zones are more organized in the case of a uniform mean shear than in isotropic turbulence, and regions of intense reaction rate appear to be associated with vortex structures such as horseshoe vortices and fingers seen in mixing layers. Concentration gradients tend to align with the direction of the most compressive principal strain rate, more so in the isotropic case.
Leonard, Andy D.↗
Active layer model for wall-bounded turbulence
The active-layer model for wall-bounded turbulence hypothesizes that the nonlinear terms are large only in a thin layer near the wall, and hence the turbulence in the region outside the active inner layer can be modeled as a linear fluctuating flow driven by the active layer. This hypothesis is tested using data obtained from a direct simulation of turbulent channel flow. It is found that the nonlinear effects are the strongest near the wall with a maximum at around y sup + = 20 and, outside the near-wall region, these involve primarily the cascading mechanism leading to dissipation.
Landahl, M. T.↗
TPSAS-NF1676L-27659-DND
The Arctic is warming fast, a trend expected to continue. Rapid Arctic warming has far-reaching implications affecting the global energy budget, sea level rise, the global carbon cycle, and potentially the atmospheric and oceanic circulation. Our understanding of the forcing mechanisms driving Arctic climate change has rapidly evolved in the last decade. For 40 years, the melting of sea ice and snow has been implicated as the primary driver of amplified Arctic warming, however our results indicate a shift, whereby changes in the atmospheric temperature, humidity, and clouds have surpassed the warming effects of snow and sea ice loss in the latest CMIP5 model simulations. We examine the impacts of climate feedbacks on the Arctic surface energy budget: namely surface albedo, clouds, shortwave clear-sky feedbacks, longwave clear-sky feedbacks, ocean heat content, and surface turbulent fluxes. Simulations from 16 CMIP5 models running RCP8.5, a high-emission scenario are used. We employ the methodology used in Lu and Cai (2009) on a single CMIP3 model allowing an analysis of changes between CMIP3 and CMIP5.
Robyn C Boeke↗