SEARCH · Search NASA
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.
Two-dimensional, spatially growing, turbulent mixing layers are simulated numerically by a vortex method and the results are compared with those determined experimentally. The effects of artificial forcing on flow development are also studied. Many of the flow features which have been observed experimentally are reproduced, and good quantitative agreements between experiments and computations are obtained.
Research on combustion is being conducted at Lewis Research Center to provide improved analytical models of the complex flow and chemical reaction processes which occur in the combustor of gas turbine engines and other aeropropulsion systems. The objective of the research is to obtain a better understanding of the various physical processes that occur in the gas turbine combustor in order to develop models and numerical codes which can accurately describe these processes. Activities include in-house research projects, university grants, and industry contracts and are classified under the subject areas of advanced numerics, fuel sprays, fluid mixing, and radiation-chemistry. Results are high-lighted from several projects.
Current analytical models for large-scale air-sea interactions in the middle latitudes are reviewed in terms of known sea-surface temperature (SST) anomalies. The scales and strength of different atmospheric forcing mechanisms are discussed, along with the damping and feedback processes controlling the evolution of the SST. Difficulties with effective SST modeling are described in terms of the techniques and results of case studies, numerical simulations of mixed-layer variability and statistical modeling. The relationship between SST and diabatic heating anomalies is considered and a linear model is developed for the response of the stationary atmosphere to the air-sea feedback. The results obtained with linear wave models are compared with the linear model results. Finally, sample data are presented from experiments with general circulation models into which specific SST anomaly data for the middle latitudes were introduced.
Two-dimensional, spatially growing, turbulent mixing layers are simulated numerically by a vortex method and the results are compared with those determined experimentally. The effects of artificial forcing on flow development are also studied. Many of the flow features which have been observed experimentally are reproduced, and good quantitative agreements between experiments and computations are obtained.
Spanwise structures of 'nominally' two-dimensional (2D), spatially growing, turbulent mixing layers are simulated numerically by a three-dimensional (3D) vortex method. Small-amplitude 3D disturbances are introduced into an otherwise 2D flow field. Results show that a large-scale spanwise variation of the flow field is produced because of amplification of initial disturbances. Pairs of counterrotating streamwise vortices are formed as a result of stretching of primary spanwise vortices. The streamwise vortices are formed at a fixed spanwise location which depends on initial disturbances. The calculated magnitude of streamwise vorticity is close to that of the spanwise vorticity. The results also suggest that the presence of background disturbances may be essential for the formation of streamwise vortices. The results are in good qualitative agreement with experiments.
A three-dimensional (3-D) analytical approach to calculate the mixing and combustion efficiency along engine combustors that use sidewall fuel injectors is described. The analysis consists of three parts: (1) an empirical correlation for fuel penetration, (2) application of integral conservation laws following jet turning and pressure equilibration, and (3) an analytical solution for the simplified three-dimensional partial differential equations to describe the downstream mixing and combustion process. Numerical results are compared to some empirical and experimental work, and favorable agreements are shown. The importance of including anisotropic turbulence and three-dimensional effects are illustrated by comparison with data and other less comprehensive analyses. The analytical solution is suitable for rapid estimates of engine performance and parametric studies of combustor designs containing transverse or angled injector ports arranged in rows and/or columns.
Pursuant to a NASA-Langley development program for a scramjet HST propulsion system entailing the optimization of the scramjet combustor's fuel-air mixing and reaction characteristics, a numerical study has been conducted of the candidate parallel fuel injectors. Attention is given to a method for flow mixing-process and combustion-efficiency enhancement in which a supersonic circular hydrogen jet coflows with a supersonic air stream. When enhanced by a planar oblique shock, the injector configuration exhibited a substantial degree of induced vorticity in the fuel stream which increased mixing and chemical reaction rates, relative to the unshocked configuration. The resulting heat release was effective in breaking down the stable hydrogen vortex pair that had inhibited more extensive fuel-air mixing.
The Kelvin Helmholtz roll up of three dimensional, temporally evolving, plane mixing layers were simulated numerically. All simulations were begun from a few low wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity profile. The spanwise disturbance wavelength was taken to be less than or equal to the streamwise wavelength associated with the Kelvin Helmholtz roll up. A standard set of clean structures develop in most of the simulations. The spanwise vorticity rolls up into a corrugated spanwise roller, with vortex stretching creating strong spanwise vorticity in a cup shaped region at the vends of the roller. Predominantly streamwise rib vortices develop in the braid region between the rollers. For sufficiently strong initial three dimensional disturbances, these ribs collapse into compact axisymmetric vortices. The rib vortex lines connect to neighboring ribs and are kinked in the opposite direction of the roller vortex lines. Because of this, these two sets of vortex lines remain distinct. For certain initial conditions, persistent ribs do not develop. In such cases the development of significant three dimensionality is delayed. When the initial three dimensional disturbance energy is about equal to, or less than, the two dimensional fundamental disturbance energy, the evolution of the three dimensional disturbance is nearly linear (with respect to the mean and the two dimensional disturbances), at least until the first Kelvin Helmholtz roll up is completed.
A finite difference, three dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. Here, equations are solved in steadily rotating reference frames by using the steady state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two.
A finite-difference, three-dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudocompressibility approach and uses an implicit-upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. In this work, the equations are solved in steadily rotating reference frames by using the steady-state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two. Included in the appendix is a paper on incompressible viscous flow through artificial heart devices with moving boundaries. Time-accurate calculations, such as impeller and diffusor interaction, will be reported in future work.
A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. Unlike the Lagrangian method previously imposed which is valid only for supersonic flows, the present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multi-dimensional discontinuities with a high level of accuracy, similar to that found in one-dimensional problems.
The mixing performance in a mixing section of a rich burn/quick mix/lean burn (RQL) combustor was calculated using a 3-D numerical model in a non-reacting environment. The numerically calculated results were compared with the measured data reported by Hatch, Sowa, Samuelsen, and Holdeman, 1992. The numerical 3-D temperature fields qualitatively agree with the experimental data. Also, the development of the mixing flow and temperature non-uniformity trends throughout the mixing section for the numerically calculated results quantitatively agree with the measured data. The numerical model predicts less mixing and enhances the temperature gradients as compared to the measured data for the cases reported by Hatch et al. (1992) which include circular and slot orifice shapes (with different slant angles and aspect ratios). The predicted and measured results generally agree in the selection of the slanted slot orifice configuration yielding the best overall mixing performance (based on temperature uniformity) of all the configurations analyzed.
A hybrid grid technique for the solution of 2D, unsteady flows is developed. This technique is capable of handling complex, multiple component geometries in relative motion, such as those encountered in turbomachinery. The numerical approach utilizes a mixed structured-unstructured zonal grid topology along with modeling equations and solution methods that are most appropriate in the individual domains, therefore combining the advantages of both structured and unstructured grid techniques.
A unique formulation of describing fluid motion is presented. The method, referred to as 'extended Lagrangian method', is interesting from both theoretical and numerical points of view. The formulation offers accuracy in numerical solution by avoiding numerical diffusion resulting from mixing of fluxes in the Eulerian description. Meanwhile, it also avoids the inaccuracy incurred due to geometry and variable interpolations used by the previous Lagrangian methods. The present method is general and capable of treating subsonic flows as well as supersonic flows. The method proposed in this paper is robust and stable. It automatically adapts to flow features without resorting to clustering, thereby maintaining rather uniform grid spacing throughout and large time step. Moreover, the method is shown to resolve multidimensional discontinuities with a high level of accuracy, similar to that found in 1D problems.
The mixing performance in a mixing section of a rich burn/quick mix/lean burn (RQL) combustor was calculated using a 3-D numerical model in a non-reacting environment. The numerically calculated results were compared with the measured data reported by Hatch, Sowa, Samuelsen, and Holdeman, 1992. The numerical 3-D temperature fields qualitatively agree with the experimental data. Also, the development of the mixing flow and temperature non-uniformity trends throughout the mixing section for the numerically calculated results quantitatively agree with the measured data. The numerical model predicts less mixing and enhances the temperature gradients as compared to the measured data for the cases reported by Hatch et al. (1992) which include circular and slot orifice shapes (with different slant angles and aspect ratios). The predicted and measured results generally agree in the selection of the slanted slot orifice configuration yielding the best overall mixing performance (based on temperature uniformity) of all the configurations analyzed.
I propose to continue providing observers with basic data for interpreting spectra from stars, novas, supernovas, clusters, and galaxies. These data will include allowed and forbidden line lists, both laboratory and computed, for the first five to ten ions of all atoms and for all relevant diatomic molecules. I will eventually expand to all ions of the first thirty elements to treat far UV and X-ray spectra, and for envelope opacities. I also include triatomic molecules provided by other researchers. I have also made CDs with Partridge and Schwenke's water data for work on UV stars. The line data also serve as input to my model atmosphere and synthesis programs that generate energy distributions, photometry, limb darkening, and spectra that can be used for planning observations and for fitting observed spectra. The spectrum synthesis programs produce detailed plots with the lines identified. Grids of stellar spectra can be used for radial velocity-, rotation-, or abundance templates and for population synthesis. I am fitting spectra of bright stars to test the data and to produce atlases to guide observers. For each star the whole spectrum is computed from the UV to the far IR. The line data, opacities, models, spectra, and programs are freely distributed on CDs and on my Web site and represent a unique resource for many NASA programs. I am now in full production of new line lists for atoms. I am computing all ions of all elements from H to Zn and the first 5 ions of all the heavier elements, about 800 ions. For each ion I treat as many as 61 even and 61 odd configurations, computing all energy levels and eigenvectors. The Hamiltonian is determined from a scaled-Hartree-Fock starting guess by least squares fitting the observed energy levels. The average energy of each configuration is used in computing scaled-Thomas-Fermi-Dirac wavefunctions for each configuration which in turn are used to compute allowed and forbidden transition integrals. These are multiplied into the LS allowed and forbidden transition arrays. The transition arrays are transformed to the observed coupling to yield the allowed and forbidden line lists. Results are put on the web as they are finished. Provided I get funding,there will be more than 500 million lines. I will then compare ion by ion, to all the laboratory and computed data in the literature and make up a working line list for spectrum synthesis and opacity calculations with the best available data. As the laboratory spectrum analyses are improved, I will redo the calculations with the new energy levels. My original plan when I started the new calculations was to run through all the atoms using my old Cray programs from the 1980's that were limited to 1100 x 1100 arrays in the Hamiltonian for each J. Then I would go back and rerun the more complicated cases with 3000 x 3000 arrays so that I could include many more configurations and more configuration interactions. At present I am limited to 61 even and 61 odd configurations and I try to include everything up through n = 9. The current program runs on Alpha workstations. I decided to test the big program on Fe I and Fe II to see whether there was any great difference in the low configurations compared to those from the Cray program. Besides increasing the number of E1 lines by a factor of 6 to 7.7 million, there was an unexpected result: the electric quadrupole transitions were 10 times stronger than before because the transition integrals are weighted by r(exp 2) ---they become very large for high n, and because there are numerous configuration interactions that mix the low and high configurations. As a check I was able to reproduce Carstang's (1962) lower results by running his three configurations with my program. Since my model atom is still only a subset of a real Fe II ion, the true quadrupole A values are probably larger than mine. The magnetic dipole lines are affected by the mixing but the overall scale does not change. Because of this scovery I decided that there was no point in computing the small array cases. I have been running with as many configurations as I can and with thousands of parameters in the Hamiltonian. The computer runs take much longer to set up and produce than I had expected. I have concentrated on redoing the low iron group spectra, especially to get data for supernova modelers. I have done only Ca I -- Zn I, Ca II -- Zn II, CU I -- Cu XXIX, Zn I - Zn XXX, for practice at high stages of ionization, C I, C II, S I, and CL I and Ag I for people who were working on the laboratory spectra. Check my web site kurucz.harvard.edu for current additions. My latest calculations have been for carbon I and sulphur I, and silicon I is under way using the same elaborate approach as for C I, which took many months to do. These line lists greatly increase the number of lines in the ultraviolet, in the visible, and especially in the infrared. They will increase the opacity in A, F, and G stars. They will account for many unidentified lines in the sun.
A model of a temporal three-dimensional mixing layer laden with a multitude of drops of a fuel containing a large number of species is derived.
This paper presents a set of universal two-phase convection heat transfer correlations for modeling boiling heat transfer during pipe chilldown fit over the widest available range of cryogenic fluids and thermodynamic conditions. The correlations improve upon prior correlations that were developed separately for liquid nitrogen (LN2) and liquid hydrogen (LH2) pipe quenching datasets. The new correlations include equations to calculate the single-phase vapor heat transfer, film boiling heat transfer, transition boiling heat transfer, nucleate boiling heat transfer, single-phase liquid heat transfer, bulk vapor temperature during high quality film boiling, Leidenfrost temperature, critical heat flux, critical heat flux temperature, and the onset of nucleate boiling temperature. The correlations were validated against LH2, LN2, liquid methane, liquid oxygen, and liquid argon pipe quenching datasets. The eight datasets cover the following parameter ranges: pipe lengths of 0.1 to 6.5m; outer pipe diameters of 12.7 to 25.4 mm; pipe wall thicknesses from 0.51 to 1.64 mm; flow directions of upward, downward, and horizontal; gravity levels of 1g and 0g±0.01g. A numerical model assumes homogeneous mixing between the vapor and liquid and implicitly integrates the coupled energy equations for the pipe and fluid, as well as continuity for the fluid. The model was used to estimate the equilibrium quality and fluid mass accumulation along the pipes during chilldown from estimates of the fluid-pipe heat transfer extracted from temperature measurements. The correlations can be implemented in lumped parameter codes such as SINDA/FLUINT and the Generalized Fluid System Simulation Program (GFSSP) to improve the accuracy of chilldown time and chilldown boiloff mass predictions. Such predictions are useful for designing ground or in-space cryogenic liquid transfer systems.