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.
Presentation to be given at the FY 2020 SDRD Annual Meeting, to be held Sept. 23-24, on WebEx. Using CARS to Determine the Dynamic Temperature of Ejecta Particles Reacting with Surrounding Gas in a Shock Compression ExperimentSTL-042-20, Year 2
Various chemically reacting flow problems highlighting chemical and physical fundamentals rather than flow geometry are presently investigated by means of a comprehensive mathematical model that incorporates multicomponent molecular diffusion, complex chemistry, and heterogeneous processes, in the interest of obtaining sensitivity-related information. The sensitivity equations were decoupled from those of the model, and then integrated one time-step behind the integration of the model equations, and analytical Jacobian matrices were applied to improve the accuracy of sensitivity coefficients that are calculated together with model solutions.
Explore the source record for details and available documents.
Presentation to be presented during the Site-Directed Research and Development (SDRD) program FY 2021 review meeting (Webex), September 22–23, 2021.
The overall contractual effort and the theory and numerical solution for the Reacting and Multi-Phase (RAMP2) computer code are described. The code can be used to model the dominant phenomena which affect the prediction of liquid and solid rocket nozzle and orbital plume flow fields. Fundamental equations for steady flow of reacting gas-particle mixtures, method of characteristics, mesh point construction, and numerical integration of the conservation equations are considered herein.
All of the elements used in the Reacting and Multi-Phase (RAMP2) computer code are described in detail. The code can be used to model the dominant phenomena which affect the prediction of liquid and solid rocket nozzle and orbital plume flow fields.
A computer program written in conjunction with the numerical solution of the flow of chemically reacting gas-particle mixtures was documented. The solution to the set of governing equations was obtained by utilizing the method of characteristics. The equations cast in characteristic form were shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The characteristic directions for the gas-particle system are found to be the conventional gas Mach lines, the gas streamlines and the particle streamlines. The basic mesh construction for the flow solution is along streamlines and normals to the streamlines for axisymmetric or two-dimensional flow. The analysis gives detailed information of the supersonic flow and provides for a continuous solution of the nozzle and exhaust plume flow fields. Boundary conditions for the flow solution are either the nozzle wall or the exhaust plume boundary.
A users manual for the RAMP2 computer code is provided. The RAMP2 code can be used to model the dominant phenomena which affect the prediction of liquid and solid rocket nozzle and orbital plume flow fields. The general structure and operation of RAMP2 are discussed. A user input/output guide for the modified TRAN72 computer code and the RAMP2F code is given. The application and use of the BLIMPJ module are considered. Sample problems involving the space shuttle main engine and motor are included.
Steelmaking contributes 8% to the total CO 2 emissions globally, primarily due to coal-based iron ore reduction. Clean hydrogen-based ironmaking has variable performance because the dominant gas–solid reduction mechanism is set by the defects and pores inside the mm- to nm-sized oxide particles that change significantly as the reaction progresses. While these governing dynamics are essential to establish continuous flow of iron and its ores through reactors, the direct link between agglomeration and chemistry is still contested due to missing measurements. In this work, we directly measure the connection between chemistry and agglomeration in the smallest iron oxides relevant to magnetite ores. Using synthesized spherical 10-nm magnetite particles reacting in H 2 , we resolve the formation and consumption of wüstite (Fe 1-x O)—the step most commonly attributed to whiskering. Using X-ray diffraction, we resolve crystallographic anisotropy in the rate of the initial reaction. Complementary imaging demonstrated how the particles self-assemble, subsequently react, and grow into elongated “whisker” structures. Our insights into how morphologically uniform iron oxide particles react and agglomerate in H 2 reduction enable future size-dependent models to effectively describe the multiscale aspects of iron ore reduction.
A flame propagating through an iron-dust mixture can propagate in two asymptotic regimes. When the characteristic time of heat transfer between particles is much smaller than the characteristic time of particle combustion, the flame propagates in the continuum regime where the heat released by reacting particles can be modelled as a space-averaged function. In contrast, when the characteristic time of heat transfer is much larger than the particle reaction time, the flame can no longer be treated as a continuum due to dominating effects associated with the discrete nature of the particle reaction. The discrete regime is characterized by weak dependence of the flame speed on the oxygen concentration compared to the continuum regime. The discrete regime is observed in flames propagating through an iron dust cloud within a gas mixture containing xenon, while the continuum regime is obtained when xenon is substituted with helium.
There are indications (Greenberg et al., 1988), that fluffy (i.e., porous) particles are responsible for the observed 3.4 and 10 micron emissions of comet Halley. The absorption characteristics of small particles both solid and fluffy are needed in order to explain the Halley emissions. How isolated small solid particles react to an external radiation field is well known - the Rayleigh approximation. How these same small particles emit when assembled as fluffy aggregates in another question. To what degree are the emission spectra of isolated and aggregated particles comparable. In order to quantify the assertion that fluffy particles produce the observed Halley infrared emission features, the authors are performing calculations to determine the effect of porosity on the absorption characteristics of aggregates of interstellar grain-type particles. The calculations are based on an integral representation of the scattered electromagnetic field. Results are given with application to comet Halley.
The feasibility of adding AlCrTiV high entropy alloy (HEA) particles to molten aluminum to make metal matrix composites was studied. Equiatomic AlCrTiV HEA particles, produced by a solid state reduction process, were dispersed into pure aluminum and 7075 alloy (Al-5.8Zn-2.3Mg-1.4Cu, all weight %) melt and by a ultrasonic dispenser. It was found that the AlCrTiV particles reacted with molten aluminum, leading to the formation of Al 45 (Cr,V) 7 and Al 3 Ti intermetallic phases in the aluminum matrix composites. Although the AlCrTiV HEA particles did not survive in molten aluminum during composite processing, the final composite materials showed significantly improved mechanical properties, due to grain refinement and dispersion strengthening. Finally, the research suggests that the interactions between matrix aluminum alloys and the particles are key in controlling the final microstructure and properties of the composites.
The use of inert and redox-active particles for high-temperature energy storage requires the development of components that can efficiently transfer energy to high-pressure working fluids like supercritical carbon dioxide (sCO 2 ). Dilute flow reactors can enable high working fluid outlet temperatures and minimal parasitic losses compared to moving packed bed and fluidized bed reactors. This research uses both computational and experimental methods to explore the design trade-offs and practical challenges of a novel component for transferring energy from dilute flows of hot, reduced metal oxide (MO x ) particles to sCO 2 in tubes. A discretized thermal resistance network model, which accounts for particle hydrodynamics, multi-mode heat transfer, and reaction equilibrium, guides the design of a prototype device. This device is experimentally tested with a surrogate heat transfer fluids and inert particle temperatures up to 400°C and a heat duty exceeding 1 kW. The data are used to validate the thermal hydraulic sub-models, allowing for the simulation of reacting particle scenarios. Under nominal design conditions, the flow rate of reactive particles is predicted to be 30% lower than that of inert particles for the same energy recovered, with over 70% of the stored particle energy transferred to the sCO 2 . Furthermore, these findings can inform the design of more efficient energy recovery reactors for particle-based systems and can be integrated into system-level concentrated solar power models with thermal storage to optimize operating conditions.
Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO 2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. Here, in this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N 2 and CO 2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.
A numerical solution for chemically reacting supersonic gas-particle flows in rocket nozzles and exhaust plumes was described. The gas-particle flow solution is fully coupled in that the effects of particle drag and heat transfer between the gas and particle phases are treated. Gas and particles exchange momentum via the drag exerted on the gas by the particles. Energy is exchanged between the phases via heat transfer (convection and/or radiation). Thermochemistry calculations (chemical equilibrium, frozen or chemical kinetics) were shown to be uncoupled from the flow solution and, as such, can be solved separately. The solution to the set of governing equations is obtained by utilizing the method of characteristics. The equations cast in characteristic form are shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The particle distribution is represented in the numerical solution by a finite distribution of particle sizes.
A series of computer codes has been developed to predict gas-particle flows and resulting impingement forces, moments and heating rates to surfaces immersed in the flow. The gas-particle flow solution is coupled via heat transfer and drag between the phases with chemical effects included in the gas phase. The flow solution and impingement calculations are discussed. Analytical results are compared with test data obtained to evaluate gas-particle effects on the Space Shuttle thermal protection system during the staging maneuver.
A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.
A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.