SEARCH · Engineering Papers
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.
Papers are presented under the following subject areas: small scales; turbulence physics; compressible flow and modeling; and reacting flows and combustion.
The high total temperatures or total enthalpies required to duplicate the high-speed flight conditions in ground experiments often place stringent requirements on the material selection and cooling needs for the test articles and intrusive flow diagnostic equipment. Furthermore, for internal flows, these conditions often complicate the use of nonintrusive diagnostics that need optical access to the test section and interior portions of the flowpath. Because of the technical challenges and increased costs associated with experimentation at high values of total enthalpy, an attempt is often made to reduce it. This is the case for the Enhanced Injection and Mixing Project (EIMP) currently underway in the Arc-Heated Scramjet Test Facility at the NASA Langley Research Center. The EIMP aims to investigate supersonic combustion ramjet (scramjet) fuel injection and mixing physics, improve the understanding of underlying physical processes, and develop enhancement strategies and functional relationships between mixing performance and losses relevant to flight Mach numbers greater than 8. The experiments will consider a "direct-connect" approach and utilize a Mach 6 nozzle to simulate the combustor entrance flow of a scramjet engine. However, while the value of the Mach number is matched to that expected at the combustor entrance in flight, the maximum value of the total enthalpy for these experiments is limited by the thermal-structural limits of the uncooled experimental hardware. Furthermore, the fuel simulant is helium, not hydrogen. The use of "cold" flows and non-reacting mixtures of fuel simulants for mixing experiments is not new and has been extensively utilized as a screening technique for scramjet fuel injectors. In this study, Reynolds-averaged simulations are utilized (RAS) to systematically verify the implicit assumptions used by the EIMP. This is accomplished by first performing RAS of mixing for two injector configurations at planned nominal experimental conditions. The mixing parameters of interest, such as mixing efficiency and total pressure recovery, are then computed and compared to the values obtained from RAS under the true enthalpy conditions and using helium and hydrogen. Finally, the impact of combustion on mixing, often deemed small enough to neglect at hypervelocity conditions, is assessed by comparing the results obtained from the hydrogen-fueled reacting and non-reacting RAS. For reacting flows, in addition to mixing efficiency and total pressure recovery, the combustion efficiency and thrust potential are also considered. In all of the simulations, the incoming air Mach number and the fuel-to-air ratio are the same, while the total pressure, total enthalpy, and the fuel simulant vary depending on the case considered. It is found that under some conditions the "cold" flow experiments are a good approximation of the flight.
While hydrocarbon combustion and gas-surface interactions have been simulated for a variety of applications in continuum solvers, such processes have been rarely investigated in Direct Simulation Monte Carlo (DSMC). Thus Variable Hard/Soft Sphere (VHS/VSS)collision parameters are not available for most molecules of interest in the decomposition of methane, a significant pyrolysis gas product and potential atmospheric species for Titan re-entries. As these properties are essential for accurate modeling of gas transport and thermochemical effects, a method must be devised to compute them for use in DSMC. In this study, the methodology outlined in Stephani et al is utilized to compute the VHS/VSS parameters by fitting the Ω(1,1)and Ω(2,2) collision integrals [1]. These collision integrals are typically computed from intermolecular potentials and are used in CFD simulations to calculate transport properties such as viscosity, thermal conductivity, and diffusion. This approach ensures consistency between micro-scale collisions and macro-scale transport properties. To produce the necessary collision integrals needed for the fitting routine, we utilize the method described by Laricchiuta et al, where collision integrals are computed based on the constituent particles’ polarizability and effective electron number [2].Preliminary computations of pure-species Ω(1,1)and Ω(2,2)collision integrals produced from the coupled approach of Stephani et al and Laricchiuta et al for methane show good agreement with those computed via available ab-initio Potential Energy Surface Computations, as seen in Figure 1. Above 300K, no more than 12% deviation is observed between the two results over the temperature range of the fitted PES data [3].Current results of this study include full collision-specificVSS transport properties for a complex-species mixture involving major byproducts of pure methane decomposition. Targeted final results will include a comprehensive list of VSS parameters relevant toDSMC simulations of Titan atmospheric entry, based on the reduced kinetic model of Savajano et al [4].
A simple model of compressible reacting flow is studied. First, a dispersion relation is derived for the linearized problem making a distinction between frozen and equilibrium sound speed. Second, the stability of the Von Neumann-Richtmyer scheme applied to this model is studied. A natural generalization of the C.F.L. condition is found.
This report covers work done on the research project Analysis and Computation of Internal Flow Field in a Scramjet Engine. The work is supported by the NASA Langley Research Center (Computational Methods Branch of the High-Speed Aerodynamics Division) through research grant NAG1-423. The governing equations of two-dimensional chemically-reacting flows are presented together with the global two-step chemistry model. The finite-difference algorithm used is illustrated and the method of circumventing the stiffness is discussed. The computer program developed is used to solve two model problems of a premixed chemically-reacting flow. The results obtained are physically reasonable.
Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the nonequilibrium air results.
The possible means of simulating nonequilibrium reacting flows in hypersonic environments, and the required diagnostic techniques, are surveyed in two categories: bulk flow behavior and determination of chemical rate parameters. Flow visualization of shock shapes for validation of computational-fluid dynamic calculations is proposed. The facilities and the operating conditions necessary to produce the required nonequilibrium conditions, the suitable optical techniques, and their sensitivity requirements, are surveyed. Shock-tubes, shock-tunnels, and ballistic ranges in a wide range of sizes and strengths are found to be useful for this purpose, but severe sensitivity requirements are indicated for the optical instruments, which can be met only by using highly-collimated laser sources. Likewise, for the determination of chemical parameters, this paper summarizes the quantities that need to be determined, required facilities and their operating conditions, and the suitable diagnostic techniques and their performance requirements. Shock tubes of various strengths are found to be useful for this purpose. Vacuum ultraviolet absorption and fluorescence spectroscopy and coherent anti-Stokes Raman spectroscopy are found to be the techniques best suited for the measurements of the chemical data.
Fundamental aspects of spectral methods are introduced. Recent developments in spectral methods are reviewed with an emphasis on collocation techniques. Their applications to both compressible and incompressible flows, to viscous as well as inviscid flows, and also to chemically reacting flows are surveyed. The key role that these methods play in the simulation of stability, transition, and turbulence is brought out. A perspective is provided on some of the obstacles that prohibit a wider use of these methods, and how these obstacles are being overcome.
The present conference discusses low density aerothermodynamics, the drag of bodies in rarefied hypersonic flow, transitional hypervelocity aerodynamic simulation and scaling, high temperature kinetics and transport properties, electron-nitrogen molecule collisions in high temperature nonequilibrium air, theoretical studies of dissociative recombination, the transport properties of some atom/ion interactions in air, and the interaction energies, dipole transition moments, and transport cross sections of N(+)-N and O(+)-O. Also discussed are the results of studies of potential fluid mechanisms for enhanced stagnation heating, vortex-induced leeward heating on a biconic at Mach 6 and 10, the effects of surface discontinuities on convective heat transfer in hypersonic flow, computational convergence in chemical reacting flows, and three-dimensional viscous shock layer applications for the Space Shuttle Orbiter.
Fundamental aspects of spectral methods are introduced. Recent developments in spectral methods are reviewed with an emphasis on collocation techniques. Their applications to both compressible and incompressible flows, to viscous as well as inviscid flows, and also to chemically reacting flows are surveyed. The key role that these methods play in the simulation of stability, transition, and turbulence is brought out. A perspective is provided on some of the obstacles that prohibit a wider use of these methods, and how these obstacles are being overcome.
The construction of interferograms, schlieren, and shadowgraphs from computed flowfield solutions permits one-to-one comparisons of computed and experimental results. A method for constructing these images from both ideal- and real-gas, two- and three-dimensional computed flowfields is described. The computational grids can be structured or unstructured, and multiple grids are an option. Constructed images are shown for several types of computed flows including nozzle, wake, and reacting flows; comparisons to experimental images are also shown. In addition, the sensitivity of these images to errors in the flowfield solution is demonstrated, and the constructed images can be used to identify problem areas in the computations.
The construction of interferograms, schlieren, and shadowgraphs from computed flowfield solutions permits one-to-one comparisons of computed and experimental results. A method of constructing these images from both ideal- and real-gas, two and three-dimensional computed flowfields is described. The computational grids can be structured or unstructured, and multiple grids are an option. Constructed images are shown for several types of computed flows including nozzle, wake, and reacting flows; comparisons to experimental images are also shown. In addition, th sensitivity of these images to errors in the flowfield solution is demonstrated, and the constructed images can be used to identify problem areas in the computations.
An approximation method for evaluation of the caloric equations used in combustion chemistry simulations is described. The method is applied to generate the equations of specific heat, static enthalpy, and Gibb's free energy for fuel mixtures of interest to gas turbine engine manufacturers. Liquid-phase fuel properties are also derived. The fuels investigated include JP-8, synthetic fuel, and two blends of JP-8 and synthetic fuel. The complete set of fuel property equations for both phases are implemented into a computational fluid dynamics (CFD) flow solver database, and multiphase, reacting flow simulations of a well-tested liquid-fueled combustor are performed. The simulations are a first step in understanding combustion system performance and operational issues when using alternate fuels, at practical engine operating conditions.
Computational Fluid Dynamics (CFD) validation experiments at NASA Lewis are described. The material presented summarized the research in 3 areas: Inlets, ducts and nozzles; Turbomachinery; and Chemically reacting flows. The specific validation activities are concerned with shock boundary layer interactions, vortex generator effects, large low speed centrifugal compressor measurements, transonic fan shock structure, rotor/stator kinetic energy distributions, stator wake shedding characteristics, boundary layer transition, multiphase flow and reacting shear layers. These experiments are intended to provide CFD validation data for the internal flow fields within aerospace propulsion system components.
Computational Fluid Dynamics (CFD) validation experiments at NASA Lewis Research Center are described. The material presented summarizes the research in three areas: Inlets, Ducts and Nozzles; Turbomachinery; and Chemically Reacting Flows. The specific validation activities are concerned with shock-boundary layer interactions, vortex generator effects, large low speed centrifugal compressor measurements, transonic fan shock structure, rotor/stator kinetic energy distributions, stator wake shedding characteristics, boundary layer transition, multiphase flow and reacting shear layers. These experiments are intended to provide CFD validation data for the internal flow fields within aerospace propulsion system components.
Analyses and numerical procedures are presented to investigate the radiative interactions of absorbing-emitting species in chemically reacting supersonic flow in various ducts. Specific attention is directed in investigating the radiative contributions of H2O, OH, and NO under realistic physical and flow conditions. The radiative interactions in reacting flows are investigated by considering the supersonic flow of premixed hydrogen and air in a channel with a compression corner at the lower boundary. The results indicate that radiation can have significant influence on the flowfield and species production depending on the chemistry model employed.
Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.
Interest in circulating fluidized bed (CFB) boilers as a power generation technology has sky-rocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. Herin this paper summarizes the modeling effort on a 50 kW th CFB combustor designed, built, and operated at CanmetENERGY in Ottawa, Canada. The numerical model employs the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite. The MFiX-PIC model parameters for the simulation are tuned against cold-flow experiments from CanmetENERGY using olivine sand as the inert bed material. It is shown that for the relatively coarse fluid meshes and large parcel sizes necessitated by the scale of the simulation, filter size dependent corrections to the drag law must be incorporated to ensure accuracy of the simulation results. The validated cold flow model is extended to simulate reacting flow with torrefied hardwood as the feedstock and to validate the combustion reaction scheme. The species concentrations at the riser outlet are compared against CanmetENERGY’s experiments and show satisfactory agreement. The simulations demonstrate the ability of MFiX-PIC to accurately capture both the physics and chemistry of a CFB combustor at bench scales, which can be further extended to pilot- and industrial-scale systems.