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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

A basic code for the prediction of transient three-dimensional turbulent flowfields

A primitive pressure-velocity finite difference code has been developed to predict transient three-dimensional turbulent flow. The code is a simplified yet effective prediction procedure for use by persons with little experience in computational fluid dynamics, and into which user-oriented complexities can easily be added. The method is based on the transient two-dimensional Los Alamos SOLA prediction technique for laminar flows. Turbulence is simulated by means of a two-equation turbulence model; species diffusion and buoyancy are also included. Two applications of the code are presented to local destratification near the release structure of a reservoir and to the deflection of a jet entering normally into a uniform cross-flow. Predicted results exhibit good agreement with experimental data, showing that a useful characterization of fully three-dimensional flows is now available.

Busnaina, A. A.↗
NASA NTRS1984

Simulation Applications at NASA Ames Research Center

Aeronautical applications of simulation technology at Ames Research Center are described. The largest wind tunnel in the world is used to determine the flow field and aerodynamic characteristics of various aircraft, helicopter, and missile configurations. Large computers are used to obtain similar results through numerical solutions of the governing equations. Capabilities are illustrated by computer simulations of turbulence, aileron buzz, and an exhaust jet. Flight simulators are used to assess the handling qualities of advanced aircraft, particularly during takeoff and landing.

Inouye, M.↗
NASA NTRS1989

Computational simulation of transition to turbulence through inverse modeling

The present investigation has focused on a computational methodology for the fundamental case of transition in channel flow, in which recently published experimental data are utilized both as a stimulus and as a measure of merit of the method. The research has proceeded along three avenues in parallel. The first task has consisted of the development and verification of a computer code which calculates the mean evolution of flow in a channel similar to the one employed experimentally by Blair and Anderson. An analytical test case was created for the dual purposes of code verification and of highlighting the interactions between the Reynolds stress and the mean velocity profile. This test case generated a Reynolds stress by the residue in the momentum equation which is produced by a typical analytical velocity profile. By a substitution of this Reynolds stress into the appropriate code module, the correctness of the code may be verified, along with the accuracy of the computational method. The second task pursued has involved the development of a triple layer model for the Reynolds stress profile, which was suggested and derived from experimental velocity profiles. It is demonstrated that the innermost length scale is based on the local friction velocity, the intermediate layer corresponds to the usual logarithmic law of the wall region in which the normalized Reynolds stress is approximately unity, and the outermost layer is represented by a closed mathematical form depending explicitly on the velocity profile in the wake region. The third task was comprised of scrutiny of the excellent databases developed by Blair and others, and the planning of its incorporation into the transition analysis. These extensive measurements indicate that turbulent statistics in the transition regime may be considered to alternate between laminar and fully turbulent types, the proportions of which are quantified by a measured intermittency function.

Sepri, Paavo↗
NASA NTRS1992

A comparative study of turbulence models for overset grids

The implementation of two different types of turbulence models for a flow solver using the Chimera overset grid method is examined. Various turbulence model characteristics, such as length scale determination and transition modeling, are found to have a significant impact on the computed pressure distribution for a multielement airfoil case. No inherent problem is found with using either algebraic or one-equation turbulence models with an overset grid scheme, but simulation of turbulence for multiple-body or complex geometry flows is very difficult regardless of the gridding method. For complex geometry flowfields, modification of the Baldwin-Lomax turbulence model is necessary to select the appropriate length scale in wall-bounded regions. The overset grid approach presents no obstacle to use of a one- or two-equation turbulence model. Both Baldwin-Lomax and Baldwin-Barth models have problems providing accurate eddy viscosity levels for complex multiple-body flowfields such as those involving the Space Shuttle.

Renze, Kevin J.↗
NASA NTRS1994

Large Eddy Simulation of Homogeneous Rotating Turbulence

Study of turbulent flows in rotating reference frames has proven to be one of the more challenging areas of turbulence research. The large number of theoretical, experimental, and computational studies performed over the years have demonstrated that the effect of solid-body rotation on turbulent flows is subtle and remains exceedingly difficult to predict. Because of the complexities associated with non-homogeneous turbulence, it is worthwhile to examine the effect of steady system rotation on the evolution of an initially isotropic turbulent flow. The assumption of statistical homogeneity considerably simplifies analysis and computation; calculation of homogeneous turbulence is further motivated since it possesses the essential physics found in more complex rotating flows. The principal objectives of the present study have therefore been to increase our fundamental understanding of turbulent flows in rotating reference frames through an examination of the asymptotic state of homogeneous rotating turbulence; particularly as to the existence of an asymptotic state which is self similar. Knowledge of an asymptotic similarity state permits prediction of the ultimate statistical evolution of the flow without requiring detailed knowledge of the complex, and not well understood, non-linear transfer processes. Aside from examination of possible similarity states in rotating turbulence, of further interest in this study has been an examination of the degree to which solid-body rotation induces a two-dimensional state in an initially isotropic flow.

Squires, Kyle D.↗
NASA NTRS2004

Large-Eddy Simulation of Laminar-Turbulent Transition in a Swept-Wing Boundary Layer

A continuous adjoint approach for obtaining sensitivity derivatives on unstructured grids is developed and analyzed. The derivation of the costate equations is presented, and a second-order accurate discretization method is described. The relationship between the continuous formulation and a discrete formulation is explored for inviscid, as well as for viscous flow. Several limitations in a strict adherence to the continuous approach are uncovered, and an approach that circumvents these difficulties is presented. The issue of grid sensitivities, which do not arise naturally in the continuous formulation, is investigated and is observed to be of importance when dealing with geometric singularities. A method is described for modifying inviscid and viscous meshes during the design cycle to accommodate changes in the surface shape. The accuracy of the sensitivity derivatives is established by comparing with finite-difference gradients and several design examples are presented.

Huai, X.↗
NASA NTRS1990

Subgrid-scale backscatter in transitional and turbulent flows

Most subgrid-scale (SGS) models for large-eddy simulations are absolutely dissipative (that is, they remove energy from the large scales at each point in the physical space). The actual SGS stresses, however, may transfer energy to the large scales (backscatter) at a given location. Direct numerical simulations of turbulent channel flow and compressible isotropic turbulence are used to study the backscatter phenomena. In all flows considered roughly 50 percent of the grid points were experiencing backscatter when a Fourier cutoff filter was used. The backscatter fraction was less with a Gaussian filter, and intermediate with a box filter in physical space. Moreover, the backscatter and forward scatter contributions to the SGS dissipation were comparable, and each was often much larger than the total SGS dissipation. The SGS dissipation (normalized by total dissipation) increased with filter width almost independently of filter type and Reynolds number. The amount of backscatter showed an increasing trend with Reynolds numbers. In the near-wall region of the channel, events characterized by strong Reynolds shear stress correlated fairly well with areas of high SGS dissipation (both forward and backward). In compressible isotropic turbulence similar results were obtained, independent of fluctuation Mach number.

Piomelli, Ugo↗
NASA NTRS1990

On the Subgrid-Scale Modeling of Compressible Turbulence

A new sub-grid scale model is presented for the large-eddy simulation of compressible turbulence. In the proposed model, compressibility contributions have been incorporated in the sub-grid scale eddy viscosity which, in the incompressible limit, reduce to a form originally proposed by Smagorinsky (1963). The model has been tested against a simple extension of the traditional Smagorinsky eddy viscosity model using simulations of decaying, compressible homogeneous turbulence. Simulation results show that the proposed model provides greater dissipation of the compressive modes of the resolved-scale velocity field than does the Smagorinsky eddy viscosity model. For an initial r.m.s. turbulence Mach number of 1.0, simulations performed using the Smagorinsky model become physically unrealizable (i.e., negative energies) because of the inability of the model to sufficiently dissipate fluctuations due to resolved scale velocity dilations. The proposed model is able to provide the necessary dissipation of this energy and maintain the realizability of the flow. Following Zeman (1990), turbulent shocklets are considered to dissipate energy independent of the Kolmogorov energy cascade. A possible parameterization of dissipation by turbulent shocklets for Large-Eddy Simulation is also presented.

Squires, Kyle↗
NASA NTRSDate not supplied

Turbulence Effects on Shaped Booms: Propagation Simulations Using KZKFourier

Upcoming X-59 aircraft flight tests as part of NASA’s Quesst Mission are expected to occur in a range of atmospheric conditions. Sensitivity of ground waveform acoustic metrics to turbulent perturbations during propagation complicates the determination of noise levels. A series of simulations through turbulence was executed using the KZKFourier model of Stout et al. to develop a database of results for expanding the functionality of NASA tools for estimating turbulence effects on shaped-boom ground waveforms. For 45 cases covering a seven-factor design space, multiple realizations of turbulence were generated to characterize statistical turbulence effects on levels of six acoustic metrics. An individual simulation in the study produced over one thousand waveforms across a virtual microphone array. Different measures were used to evaluate refinement of results with increasing number of simulations. Data suggest that mean effects and variability were most strongly influenced by propagation distance and velocity fluctuation intensity. Although local increases in acoustic metrics were common, the overall average result in all cases was a reduction in metric levels. Sensitivity of acoustic metrics varied, with mean reductions of up to 2 dB in Perceived Level occurring across the 45 cases.

turbulence↗
NASA NTRS1991

Simulation and modeling of homogeneous compressible turbulence under isotropic mean compression

Compressible homogeneous turbulence subjected to isotropic mean compression is simulated numerically for high and low turbulent Mach numbers, at various compression rates. It is found that at low Mach numbers, the effects of viscosity variations on the development of the turbulent kinetic energy can be significant. A model consistent with the invariance of the Navier-Stokes equations to spherical compression, that takes into account variable viscosity effects is proposed. At high Mach numbers, the distribution of energy between the acoustic and solenoid fields is found to depend upon initial conditions. The simulation results are also used to evaluate two-equation compressible turbulence models.

Coleman, G. N.↗
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
NASA NTRS1974

Development and application of a non-Gaussian atmospheric turbulence model for use in flight simulators

A method is described for generating time histories which model the frequency content and certain non-Gaussian probability characteristics of atmospheric turbulence including the large gusts and patchy nature of turbulence. Methods for time histories using either analog or digital computation are described. A STOL airplane was programmed into a 6-degree-of-freedom flight simulator, and turbulence time histories from several atmospheric turbulence models were introduced. The pilots' reactions are described.

Reeves, P. M.↗
NASA NTRS1994

Numerical investigation of supersonic turbulent boundary layers with high wall temperature

A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Guo, Y.↗
NASA NTRS1992

Computed Turbulent Free Shear Flow Of Air

Standard k-epsilon model of turbulence yields fairly accurate results. Symposium paper discusses numerical simulation of turbulent free shear flow of nonreacting compressible fluid. Ability to compute such flows essential to advances in design.

Viegas, J. R.↗
NASA NTRS2000

Comment on "Heterodyne Lidar Returns in the Turbulent Atmosphere: Performance Evaluation of Simulated Systems"

The explanation for the difference between simulation and the zero-order theory for heterodyne lidar returns in a turbulent atmosphere proposed by Belmonte and Rye is incorrect. The theoretical expansion is not developed under a square- law-structure function approximation (random wedge atmosphere). Agreement between the simulations and the zero-order term of the theoretical expansion is produced for the limit of statistically independent paths (bi-static operation with large transmitter-receiver separation) when the simulations correctly include the large-scale gradients of the turbulent atmosphere.

Frehlich, Rod↗
NASA NTRS1990

Physical considerations governing computation of turbulent flows over bodies at large incidence

The physical aspects governing accurate numerical simulation of turbulent flows having large regions of crossflow separation are re-examined. Time-accurate, three-dimensional fine-grid Navier-Stokes solutions were obtained for turbulent subsonic and supersonic flows over a slender ogive-cylinder body of revolution at large angles of attack. These flowfields are complex and contain regions of crossflow separation and an organized leeward-side vortex structure. An algebraic eddy-viscosity turbulence model has been modified to correctly account for the effects of the vortices on the underlying viscous layers. The numerical results show that the vortical flow structure of both subsonic and supersonic high incidence flows are similar, and can be accurately treated by the turbulence model. As the angle of attack is increased, the effectiveness of the model increases since the influence of the vortices on the underlying boundary layer decreases.

Degani, David↗
NASA NTRSDate not supplied

An Improved Approach to the Predictability & Reliability of the Onset of Turbulence With Shocks

The construction of numerical schemes for (a) stable and accurate simulation of turbulence with strong shocks, and for (b) obtaining correct propagation speed of discontinuities in the presence of stiff source terms share one important ingredient – minimization of numerical dissipation while maintaining numerical stability. The dual requirements to achieve both numerical stability and minimal numerical dissipation are often conflicting since existing shock capturing schemes were designed mainly to be robust for rapidly developed turbulence-free flows and for shock waves without stiff source term. For the past two decades, Yee and collaborators have focused on an improved understanding of the nonlinear behavior of different high order shock-capturing methods. It was found that even very high order methods without proper nonlinear stability and numerical dissipation control can either numerically smear the onset of turbulence due to excess numerical dissipation, or induce (onset) numerical turbulence that is not physical turbulence due to lack of proper numerical dissipation to improve nonlinear stability for long time integration. Our approach is to combine (I) and (II) below for obtaining the physically correct onset of turbulence with shocks, including problems with stiff source terms: (I) Nonlinear dynamics is utilized to complement the traditional linearized stability theory (Yee & Sweby, Yee et al., Griffiths et al., Lafon & Yee, Yee, Wang et al., Kotov et al. 1990- 2015) in order to (i) Minimize numerically induced false transition to turbulence, (ii) Minimize numerical instability due to long time integration of turbulent flows, (iii) Minimize numerically induced standing wave solutions, and (iv) Minimize wrong propagation of speed of discontinuities due to the presence of stiff source terms. (II) Our recently developed physical preserving (structural preserving) high order methods with improved nonlinear stability & accuracy that are essential in minimizing spurious numerics are used.

HECC↗