Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

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 1,567 records · Page 87

Numerical investigation of complex, transitional, and chemically reacting flows

This document is divided into two parts. Part A, on turbulent transport of chemical species in compressible reacting flows, includes the following topics: mixing of chemical species in reacting flows; chemically reacting turbulent flows; turbulent mixing of chemical species in reacting flows, which includes large eddy mixing, turbulent mixing, and molecular mixing; multiple time scale turbulence equations; species conservation equations for reacting flows; and a finite volume numerical method. Part B, on unsteady transitional flows over forced oscillatory surfaces, discusses unsteady transitional flow equations with moving boundaries and the application of the finite volume numerical solution method.

Kim, S.-W.↗

Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D

This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.

Pao, S. Paul↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗

CFD-Based Design Optimization for Single Element Rocket Injector

To develop future Reusable Launch Vehicle concepts, we have conducted design optimization for a single element rocket injector, with overall goals of improving reliability and performance while reducing cost. Computational solutions based on the Navier-Stokes equations, finite rate chemistry, and the k-E turbulence closure are generated with design of experiment techniques, and the response surface method is employed as the optimization tool. The design considerations are guided by four design objectives motivated by the consideration in both performance and life, namely, the maximum temperature on the oxidizer post tip, the maximum temperature on the injector face, the adiabatic wall temperature, and the length of the combustion zone. Four design variables are selected, namely, H2 flow angle, H2 and O2 flow areas with fixed flow rates, and O2 post tip thickness. In addition to establishing optimum designs by varying emphasis on the individual objectives, better insight into the interplay between design variables and their impact on the design objectives is gained. The investigation indicates that improvement in performance or life comes at the cost of the other. Best compromise is obtained when improvements in both performance and life are given equal importance.

Vaidyanathan, Rajkumar↗

LES of Temporally Evolving Mixing Layers by an Eighth-Order Filter Scheme

An eighth-order filter method for a wide range of compressible flow speeds (H.C. Yee and B. Sjogreen, Proceedings of ICOSAHOM09, June 22-26, 2009, Trondheim, Norway) are employed for large eddy simulations (LES) of temporally evolving mixing layers (TML) for different convective Mach numbers (Mc) and Reynolds numbers. The high order filter method is designed for accurate and efficient simulations of shock-free compressible turbulence, turbulence with shocklets and turbulence with strong shocks with minimum tuning of scheme parameters. The value of Mc considered is for the TML range from the quasi-incompressible regime to the highly compressible supersonic regime. The three main characteristics of compressible TML (the self similarity property, compressibility effects and the presence of large-scale structure with shocklets for high Mc) are considered for the LES study. The LES results using the same scheme parameters for all studied cases agree well with experimental results of Barone et al. (2006), and published direct numerical simulations (DNS) work of Rogers & Moser (1994) and Pantano & Sarkar (2002).

Hadjadj, A↗

Streamline‐Traced, External‐Compression Supersonic Inlets for Mach 2

A computational study was performed to explore the aerodynamic design and performance of streamline-traced, external‐compression (STEX) inlets for Mach 2.0. The performance metrics included inlet flow rates, total pressure recovery, and total pressure distortion. The study explored the use of round and flattened capture cross‐sections and a bleed slot. The design of the inlet and generation of the inlet geometry was performed using the Supersonic Inlet Design and Analysis (SUPIN) Tool. Computational grids were generated, and methods of computational fluid dynamics (CFD) were applied to solve the three dimensional, turbulent flow through the inlets using the Wind‐US CFD flow solver. It was found that the use of a bleed slot with about 5% bleed resulted in an inlet total pressure recovery of 95% at Mach 2 with acceptable radial and circumferential total pressure distortion. It was also found that there was only a slight decrease in performance between a round and flattened capture cross‐section. This suggests that the use of streamline tracing offers flexibility in shaping the inlet capture cross‐section for a more favorable integration with an aircraft wing or fuselage. This flexibility and the good performance of the inlet provides a novel inlet design for future commercial supersonic aircraft.

Computational Fluid Dynamics↗

Catalytic converters having non-linear flow channels

Disclosed is a honeycomb catalyst substrate core having geometrically non-linear flow channels. In an embodiment, the honeycomb catalyst substrate core includes helical flow channels. In another embodiment, the honeycomb catalyst substrate core includes sinusoidal flow channels. In yet another embodiment, the honeycomb catalyst substrate core includes helical plus sinusoidal flow channels. The honeycomb catalyst substrate core comprises a plurality of parallel non-linear flow channels formed along a longitudinal axis of symmetry of the catalyst substrate core, each non-linear flow channel configured such that a turbulent vortical flow occurs during engine exhaust gas flow. Also disclosed is a method for manufacturing a ceramic honeycomb having non-linear flow channels, comprising the steps extrusion soft ceramic material through a die whilst the die moves through six degrees of freedom along its axis of symmetry. Disclosure includes a method for manufacturing a ceramic honeycomb having non-linear flow channels using three-dimensional printing.

33 ADVANCED PROPULSION SYSTEMS↗

Study of effects of injector geometry on fuel-air mixing and combustion

An implicit finite-difference method has been developed for computing the flow in the near field of a fuel injector as part of a broader study of the effects of fuel injector geometry on fuel-air mixing and combustion. Detailed numerical results have been obtained for cases of laminar and turbulent flow without base injection, corresponding to the supersonic base flow problem. These numerical results indicated that the method is stable and convergent, and that significant savings in computer time can be achieved, compared with explicit methods.

Bangert, L. H.↗

Coupled Monte Carlo Probability Density Function/ SPRAY/CFD Code Developed for Modeling Gas-Turbine Combustor Flows

The success of any solution methodology for studying gas-turbine combustor flows depends a great deal on how well it can model various complex, rate-controlling processes associated with turbulent transport, mixing, chemical kinetics, evaporation and spreading rates of the spray, convective and radiative heat transfer, and other phenomena. These phenomena often strongly interact with each other at disparate time and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. Turbulence manifests its influence in a diffusion flame in several forms depending on how turbulence interacts with various flame scales. These forms range from the so-called wrinkled, or stretched, flamelets regime, to the distributed combustion regime. Conventional turbulence closure models have difficulty in treating highly nonlinear reaction rates. A solution procedure based on the joint composition probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices such as extinction, blowoff limits, and emissions predictions because it can handle the nonlinear chemical reaction rates without any approximation. In this approach, mean and turbulence gas-phase velocity fields are determined from a standard turbulence model; the joint composition field of species and enthalpy are determined from the solution of a modeled PDF transport equation; and a Lagrangian-based dilute spray model is used for the liquid-phase representation with appropriate consideration of the exchanges of mass, momentum, and energy between the two phases. The PDF transport equation is solved by a Monte Carlo method, and existing state-of-the-art numerical representations are used to solve the mean gasphase velocity and turbulence fields together with the liquid-phase equations. The joint composition PDF approach was extended in our previous work to the study of compressible reacting flows. The application of this method to several supersonic diffusion flames associated with scramjet combustor flow fields provided favorable comparisons with the available experimental data. A further extension of this approach to spray flames, three-dimensional computations, and parallel computing was reported in a recent paper. The recently developed PDF/SPRAY/computational fluid dynamics (CFD) module combines the novelty of the joint composition PDF approach with the ability to run on parallel architectures. This algorithm was implemented on the NASA Lewis Research Center's Cray T3D, a massively parallel computer with an aggregate of 64 processor elements. The calculation procedure was applied to predict the flow properties of both open and confined swirl-stabilized spray flames.

Source record↗

Significance of shock structure on supersonic jet mixing noise of axisymmetric nozzles

One of the key technical elements in NASA's high speed research program is reducing the noise level to meet the federal noise regulation. The dominant noise source is associated with the supersonic jet discharged from the engine exhaust system. Whereas the turbulence mixing is largely responsible for the generation of the jet noise, a broadband shock-associated noise is also generated when the nozzle operates at conditions other than its design. For both mixing and shock noise components, because the source of the noise is embedded in the jet plume, one can expect that jet noise can be predicted from the jet flowfield computation. Mani et al. developed a unified aerodynamic/acoustic prediction scheme by applying an extension of Reichardt's aerodynamic model to compute turbulent shear stresses which are utilized in estimating the strength of the noise source. Although this method produces a fast and practical estimate of the jet noise, a modification by Khavaran et al. has led to an improvement in aerodynamic solution. The most notable feature in this work is that Reichardt's model is replaced with the computational fluid dynamics (CFD) solution of Reynolds-averaged Navier-Stokes equations. The major advantage of this work is that the essential, noise-related flow quantities such as turbulence intensity and shock strength can be better predicted. The predictions were limited to a shock-free design condition and the effect of shock structure on the jet mixing noise was not addressed. The present work is aimed at investigating this issue. Under imperfectly expanded conditions the existence of the shock cell structure and its interaction with the convecting turbulence structure may not only generate a broadband shock-associated noise but also change the turbulence structure, and thus the strength of the mixing noise source. Failure in capturing shock structures properly could lead to incorrect aeroacoustic predictions.

Kim, Chan M.↗

Laminar and Turbulent Gaseous Diffusion Flames

Recent measurements and predictions of the properties of homogeneous (gaseous) laminar and turbulent non-premixed (diffusion) flames are discussed, emphasizing results from both ground- and space-based studies at microgravity conditions. Initial considerations show that effects of buoyancy not only complicate the interpretation of observations of diffusion flames but at times mislead when such results are applied to the non-buoyant diffusion flame conditions of greatest practical interest. This behavior motivates consideration of experiments where effects of buoyancy are minimized; therefore, methods of controlling the intrusion of buoyancy during observations of non-premixed flames are described, considering approaches suitable for both normal laboratory conditions as well as classical microgravity techniques. Studies of laminar flames at low-gravity and microgravity conditions are emphasized in view of the computational tractability of such flames for developing methods of predicting flame structure as well as the relevance of such flames to more practical turbulent flames by exploiting laminar flamelet concepts.

Faeth, G. M.↗

Center for modeling of turbulence and transition: Research briefs, 1993

This research brief contains the progress reports of the research staff of the Center for Modeling of Turbulence and Transition (CMOTT) from June 1992 to July 1993. It is also an annual report to the Institute for Computational Mechanics in Propulsion located at Ohio Aerospace Institute and NASA Lewis Research Center. The main objectives of the research activities at CMOTT are to develop, validate, and implement turbulence and transition models for flows of interest in propulsion systems. Currently, our research covers eddy viscosity one- and two-equation models, Reynolds-stress algebraic equation models, Reynolds-stress transport equation models, nonequilibrium multiple-scale models, bypass transition models, joint scalar probability density function models, and Renormalization Group Theory and Direct Interaction Approximation methods. Some numerical simulations (LES and DNS) have also been carried out to support the development of turbulence modeling. Last year was CMOTT's third year in operation. During this period, in addition to the above mentioned research, CMOTT has also hosted the following programs: an eighteen-hour short course on 'Turbulence--Fundamentals and Computational Modeling (Part I)' given by CMOTT at the NASA Lewis Research Center; a productive summer visitor research program that has generated many encouraging results; collaborative programs with industry customers to help improve their turbulent flow calculations for propulsion system designs; a biweekly CMOTT seminar series with speakers from within and without the NASA Lewis Research Center including foreign speakers. In addition, CMOTT members have been actively involved in the national and international turbulence research activities. The current CMOTT roster and organization are listed in Appendix A. Listed in Appendix B are the abstracts of the biweekly CMOTT seminar. Appendix C lists the papers contributed by CMOTT members.

Liou, William W.↗

Velocity-split Navier-Stokes solution procedure for incompressible high Reynolds number external flows

A method for solving the Navier-Stokes equations based on splitting the velocity vector into its rotational and irrotational parts was sucessfully applied to internal flow computations. The applicability of the method to external flows is examined by studying several model problems. The model problems are those of laminar and turbulent incompressible flow past a semi-infinite flat plate and laminar incompressible flow past a finite flat plate. For these problems, the procedure accurately reproduces the known solutions and is computationally very efficient even at high Reynolds numbers. Computational aspects of the method are discussed along with the possibility of using the procedure to retrofit a viscous capability into existing potential flow codes.

Dwoyer, D. L.↗

LOX manifold tee analysis

A 4000 Hz vibration phenomena was observed during the test firings of several space shuttle main engines (SSME). Experimental studies of this phenomena suggest that the problem might be associated with vortex shedding from the vanes within the LOX tee manifold. The objective of this study was to determine the unsteady, 3-D flow associated with these vanes by computational methods to identify and better understand the 4000 Hz vibration phenomena. A flow solver, FDNS, for the turbulent conservation equations was validated for predicting high frequency vortex dynamics and used to predict 2-D and 3-D flows within the LOX tee. 4000 Hz excitation oscillations were predicted for some flows and the entire 3-D flow structure was predicted for LOX tee flow. The complexity of the flow was revealed by this analysis, and computational methods for predicting these high frequency oscillations in future engine systems were established.

Chen, Yen-Sen↗

Reduction and Analysis of Phosphor Thermography Data With the IHEAT Software Package

Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.

Merski, N. Ronald↗

Conditionally-averaged structures in wall-bounded turbulent flows

The quadrant-splitting and the wall-shear detection techniques were used to obtain ensemble-averaged wall layer structures. The two techniques give similar results for Q4 events, but the wall-shear method leads to smearing of Q2 events. Events were found to maintain their identity for very long times. The ensemble-averaged structures scale with outer variables. Turbulence producing events were associated with one dominant vortical structure rather than a pair of counter-rotating structures. An asymmetry-preserving averaging scheme was devised that allowed a picture of the average structure which more closely resembles the instantaneous one, to be obtained.

Guezennec, Yann G.↗

Nonlinear global gyrokinetic delta- f turbulence simulations in a quasi-axisymmetric stellarator

We use the global gyrokinetic stellarator code XGC-S to look at turbulence physics in geometry of the proposed quasi-axisymmetric stellarator, QUASAR. In this article, the XGC-S code is first verified for nonlinear simulations of ion temperature gradient-driven turbulence by a successful benchmark with the established nonlinear tokamak code XGC1 in a CYCLONE-like tokamak case. We then apply the XGC-S code to QUASAR for nonlinear turbulence physics. We find that the highly poloidally localized linear mode gives way to a more evenly distributed nonlinear turbulent structure. Furthermore, the radial extent of the mode increases and occupies a large fraction of the volume, despite the initially localized temperature profile, as the profiles relax.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗