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.
The pyrolysis of general biomass materials is modeled via a superposition of cellulose, hemicellulose an lingnin kinetics. All three of the primary biomass components are modeled with multi-step kinetics involving both competetive primary pyrolysis and secodary tar decomposition reactions. Only typical (untreated) feedstocks are considered at atmospheric pyrolysis pressures.
The thermal decomposition of ammonium perchlorate (AP) and ammonium-perchlorate-based composite propellants is studied using the simultaneous thermogravimetric modulated beam mass spectrometry (STMBMS) technique. The main objective of the present work is to evaluate whether the STMBMS can provide new data on these materials that will have sufficient detail on the reaction mechanisms and associated reaction kinetics to permit creation of a detailed model of the thermal decomposition process. Such a model is a necessary ingredient to engineering models of ignition and slow-cookoff for these AP-based composite propellants. Results show that the decomposition of pure AP is controlled by two processes. One occurs at lower temperatures (240 to 270 C), produces mainly H2O, O2, Cl2, N2O and HCl, and is shown to occur in the solid phase within the AP particles. 200(micro) diameter AP particles undergo 25% decomposition in the solid phase, whereas 20(micro) diameter AP particles undergo only 13% decomposition. The second process is dissociative sublimation of AP to NH3 + HClO4 followed by the decomposition of, and reaction between, these two products in the gas phase. The dissociative sublimation process occurs over the entire temperature range of AP decomposition, but only becomes dominant at temperatures above those for the solid-phase decomposition. AP-based composite propellants are used extensively in both small tactical rocket motors and large strategic rocket systems.
Terrestrial meteorite impacts have been directly linked with multiple mass extinction events throughout Earth’s history. Among these impacts, those that land in sedimentary target rocks are thought to generate large quantities of CO 2 via decarbonation. This injection of CO 2 into the atmosphere has the potential to alter the climate and threaten terrestrial habitability. The magnitude of this change depends upon the net amount of CO 2 released, which is controlled by how much CO 2 is produced by the impact, how much CO 2 is removed by back-reactions after the impact, and other post-impact CO 2 sinks. To interrogate the behavior of CO 2 release into the impact atmosphere, we present carbon (δ 13 C), oxygen (δ 18 O), and clumped (Δ 47 ) isotope results from carbonate clasts preserved within the impact breccia of the Steen River Impact Structure (SRIS) in Alberta, Canada. These clasts exhibit a ~65‰ range in δ 13 C and a ~5‰ variation in δ 18 O. However, while δ 13 C and δ 18 O are positively correlated, Δ 47 unexpectedly has a negative relationship with δ 13 C and δ 18 O. Based on prior assumptions, there would either be (1) no relationship with carbonate Δ 47 and the bulk ratios because Δ 47 would be reset at extreme impact temperatures, or (2) Δ 47 would have a positive correlation with the bulk ratios, reflecting gradual Δ 47 ‘resetting.’ To reconcile the SRIS result with these expectations, we conducted a series of in vacuo heating experiments at temperatures above calcite decomposition. As predicted by Rayleigh fractionation, these heating experiments generated depletions in δ 13 C and δ 18 O that increased with reaction time. Mimicking the SRIS results, these experiments also produced concomitant increases in Δ 47 . Using an adaptation of a mechanistic model for decomposition fractionation (Hayles & Killingsworth, 2022) we hypothesize that these exotic isotope trends are caused by Rayleigh fractionation and a Δ 47 kinetic isotope effect that results from the disproportionation of O to CaO and CO 2 during thermal decomposition. This high Δ 47 CO 2 subsequently exchanges with the residual CaCO 3 . This work highlights a potential pathway for identifying and quantifying CO 2 generation by impacts and builds on the relatively limited literature characterizing the behavior of carbonate clumped isotopes at very high geologic temperatures.
The high-pressure carbon monoxide (HiPco) process, developed at Rice University, has been reported to produce single-walled carbon nanotubes from gas-phase reactions of iron carbonyl in carbon monoxide at high pressures (10-100 atm). Computational modeling is used here to develop an understanding of the HiPco process. A detailed kinetic model of the HiPco process that includes of the precursor, decomposition metal cluster formation and growth, and carbon nanotube growth was developed in the previous article (Part I). Decomposition of precursor molecules is necessary to initiate metal cluster formation. The metal clusters serve as catalysts for carbon nanotube growth. The diameter of metal clusters and number of atoms in these clusters are some of the essential information for predicting carbon nanotube formation and growth, which is then modeled by the Boudouard reaction with metal catalysts. Based on the detailed model simulations, a reduced kinetic model was also developed in Part I for use in reactor-scale flowfield calculations. Here this reduced kinetic model is integrated with a two-dimensional axisymmetric reactor flow model to predict reactor performance. Carbon nanotube growth is examined with respect to several process variables (peripheral jet temperature, reactor pressure, and Fe(CO)5 concentration) with the use of the axisymmetric model, and the computed results are compared with existing experimental data. The model yields most of the qualitative trends observed in the experiments and helps to understanding the fundamental processes in HiPco carbon nanotube production.
New CCT equations have been developed and optimized to simulate the start temperatures of the austenite decomposition process in low-alloyed steels using experimental CCT data published in the literature. Exceptionally, this optimization does not apply the nominal compositions of the steels, but the corresponding soluble compositions of the grain boundaries calculated using IDS software, depending on the reported austenitization treatments of the steels. These compositions, rather than the nominal ones, are expected to control the start of the austenite decomposition, which usually initiates at the grain boundaries. The new optimization treatment takes into account the solute microsegregation and the possible precipitate formation. Using IDS software, the new equations were validated with new experimental CCT data. Agreement was good not only for the austenite decomposition start temperatures, but also for the final phase fractions, indicating fairly reasonable predictions of phase transformation kinetics by the IDS. In addition, IDS simulations were compared with the experimental CCT data of five high-carbon steels, applying both the new equations based on grain boundary soluble compositions as well as the equations based on the nominal compositions. With the same experimental CCT data used in optimization, better agreement was obtained with the new equations, indicating the importance of determining the soluble compositions at the grain boundaries where the austenite decomposition process is likely to begin.
Competing microstructural evolution mechanisms can exist simultaneously when duplex stainless steels are operating for several decades in a high temperature service environment. Such competition between different microstructural evolution pathways can be difficult to ascertain using simple model alloy systems necessitating detailed structural and compositional analysis of phase transformation mechanisms in complex alloys Thus, in this study duplex stainless steels with complex but well understood chemistries were used to investigate the relative importance of different heterogeneous nucleation sites – specifically spinodal decomposition and Cu clustering – on precipitation of Ni-Si-Mn particles during long term thermal aging. Precipitation of Ni-Si-Mn particles in duplex stainless steels and ferritic steels during thermal aging and irradiation is known to have profound effect on mechanical properties. Using duplex stainless steels with custom modified compositions along with atom probe tomography characterization and first-passage kinetic Monte Carlo (FPKMC) simulations it is revealed that while the interface between Cr and Fe formed during spinodal decomposition can be a pathway for diffusion of Ni, Si, Mn, and Cu, it is not a site that encourages precipitation of Ni-Si-Mn rich particles. Instead, the presence of a higher concentration of Cu in these steels leads to formation of small Cu clusters with high energy interfaces that act as nucleation sites for Ni-Si-Mn particles because of its strong interaction particularly with Ni increasing the chemical driving force for precipitation. These results will help inform predictive models for the use of duplex stainless steels and other precipitation-hardened alloys for extended operation at high temperatures.
Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.
N-Nitrosamines are one of the environmentally significant byproducts from aqueous amine-based post-combustion carbon capture systems (CCS) due to their potential risk to human health. Safely mitigating nitrosamines before they are emitted from these CO 2 capture systems is therefore a key concern before widescale deployment of CCS can be used to address worldwide decarbonization goals. Electrochemical decomposition is one viable route to neutralize these harmful compounds. The circulating emission control waterwash system, commonly installed at the end of the flue gas treatment trains to minimize amine solvent emissions, plays an important role to capture N-nitrosamines and control their emission into the environment. The waterwash solution is the last point where these compounds can be properly neutralized before becoming an environmental hazard. In this study, the decomposition mechanisms of N-nitrosamines in a simulated CCS waterwash with residual alkanolamines was investigated using several laboratory-scale electrolyzers utilizing carbon xerogel (CX) electrodes. Hcell experiments revealed that N-nitrosamines were decomposed through a reduction reaction and are converted into their corresponding secondary amines thereby neutralizing their environmental impact. Batch-cell experiments statistically examined the kinetic models of N-nitrosamine removal by a combined adsorption and decomposition processes. The cathodic reduction of the N-nitrosamines statistically obeyed the first-order reaction model. Finally, a prototype flow-through reactor using an authentic waterwash was used to successfully target and decompose N-nitrosamines to below the detectable level without degrading the amine solvent compounds allowing them to be return to the CCS and lower the system operating costs. Furthermore, the developed electrolyzer was able to efficiently remove greater than 98% of N-nitrosamines from the waterwash solution without producing any additional environmentally harmful compounds and offers an effective and safe route to mitigate these compounds from CO 2 capture systems.
Abstract Biofuels are a promising ecologically viable and renewable alternative to petroleum fuels, with the potential to reduce net greenhouse gas emissions. However, biomass sourced fuels are often produced as blends of hydrocarbons and their oxygenates. Such blending complicates the implementation of these fuels in combustion applications. Variations in a biofuel’s composition will dictate combustion properties such as auto ignition temperature, reaction delay time, and reaction pathways. A handful of novel drop-in replacement biofuels for conventional transportation fuels have recently been down selected from a list of over 10,000 potential candidates as part of the U.S. Department of Energy’s (DOE) Co-Optimization of Fuels and Engines (Co-Optima) initiative. Diisobutylene (DIB) is one such high-performing hydrocarbon which can readily be produced from the dehydration and dimerization of isobutanol, produced from the fermentation of biomass-derived sugars. The two most common isomers realized, from this process, are 2,4,4-trimethyl-1-pentene (α-DIB) and 2,4,4-trimethyl-2-pentene (β-DIB). Due to a difference in olefinic bond location, the α- and β- isomer exhibit dramatically different ignition temperatures at constant pressure and equivalence ratio. This may be attributed to different fragmentation pathways enabled by allylic versus vinylic carbons. For optimal implementation of these biofuel candidates, explicit identification of the intermediates formed during the combustion of each of the isomers is needed. To investigate the combustion pathways of these molecules, tunable vacuum ultraviolet (VUV) light (in the range 8.1–11.0 eV) available at the Lawrence Berkeley National Laboratory’s Advanced Light Source (ALS) has been used in conjunction with a jet stirred reactor (JSR) and time-of-flight mass spectrometry to probe intermediates formed. Relative intensity curves for intermediate mass fragments produced during this process were obtained. Several important unique intermediates were identified at the lowest observable combustion temperature with static pressure of 93,325 Pa and for 1.5 s residence time. As this relatively short residence time is just after ignition, this study is targeted at the fuels’ ignition events. Ignition characteristics for both isomers were found to be strongly dependent on the kinetics of C 4 and C 7 fragment production and decomposition, with the tert-butyl radical as a key intermediate species. However, the ignition of α-DIB exhibited larger concentrations of C 4 compounds over C 7 , while the reverse was true for β-DIB. These identified species will allow for enhanced engineering modeling of fuel blending and engine design.
Spinodal decomposition is not the only way to make a bijel. Indeed, while spinodal decomposition produces structures with a potentially useful morphology, it can be challenging to make bijels using this method and the resulting systems can be hard to process and manipulate. Furthermore, exploiting the functional properties of the assembled particle monolayer is extremely challenging. Here, we show how the assembly of nanoparticle surfactants at the liquid–liquid interface can be used to kinetically trap liquids into a wealth of complex structures without using spinodal decomposition. We apply liquid three-dimensional printing and moulding methods, along with patterned substrates with controllable wetting properties, to build all-liquid devices with applications in chemical synthesis, separation, and purification. The functional properties of the assembled nanomaterials can be exploited to produce interfacially structured liquids that are plasmonically and magnetically responsive. Finally, we conclude by arguing that, while the field shows great promise, efforts need to be made to translate liquid bicontinuous systems out of the laboratory and into meaningful, real-world applications, as well applications in more ‘exotic’ disciplines, such as synthetic biology.
Polysulfone was shock loaded to pressures of 14.6–26.2 GPa in a series of gas gun-driven plate-impact experiments measuring material response with embedded electro-magnetic particle velocity gauges and optical velocimetry. The embedded electro-magnetic particle velocity gauges did not show a distinct two-wave structure but did show rounding that suggested a reaction but not a distinct separation of the reactants and product waves. In contrast, the transmission experiments fielded with optical velocimetry, with product pressures ranging from 21.0 to 26.2 GPa, showed well-defined two-wave structures due to shock-driven chemical decomposition of the polymer to products at a higher density. Distinct two-wave structures have not previously been published in literature during polymer compression; here, we observed these two-wave structures at both the polymer/lithium fluoride and polymer/polymethylmethacrylate interfaces. Hydrodynamic simulations of the experiments were performed using a newly constructed SESAME equation of state (EOS) for the polymer, a thermochemical EOS for the decomposition products, and an Arrhenius reaction rate model for the kinetics of the reaction. Finally, simulation results also demonstrated two-wave structures but were unable to quantitatively reproduce either the embedded gauge or optical velocimetry data.
To support the National Aerospace Plane project, the RPLUS3D CFD code has been developed at NASA Lewis. The code has the ability to solve three-dimensional flowfields with finite rate combustion of hydrogen and air. The combustion processes of the hydrogen-air system are simulated by an 18-reaction path, 8-species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though it uses a fully implicit numerical algorithm. A k-epsilon (two equation) turbulence model is incorporated into the RPLUS3D code.
The RPLUS3D code has been developed at the NASA Lewis Research Center to support the National Aerospace Plane (NASP) project. The code has the ability to solve three dimensional flowfields with finite rate combustion of hydrogen and air. The combustion process of the hydrogen-air system are simulated by an 18 reaction path, 8 species chemical kinetic mechanism. The code uses a Lower-Upper (LU) decomposition numerical algorithm as its basis, making it a very efficient and robust code. Except for the Jacobian matrix for the implicit chemistry source terms, there is no inversion of a matrix even though a fully implicit numerical algorithm is used. A k-epsilon turbulence model has recently been incorporated into the code. Initial validations have been conducted for a flow over a flat plate. Results of the validation studies are shown. Some difficulties in implementing the k-epsilon equations to the code are also discussed.
A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation, and a detailed wet-CO mechanism. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustor efficiency of an uni-element, tri-propellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.
A one-formula surrogate fuel formulation and its quasi-global combustion kinetics model are developed to support the design of injectors and thrust chambers of kerosene-fueled rocket engines. This surrogate fuel model depicts a fuel blend that properly represents the general physical and chemical properties of kerosene. The accompanying gaseous-phase thermodynamics of the surrogate fuel is anchored with the heat of formation of kerosene and verified by comparing a series of one-dimensional rocket thrust chamber calculations. The quasi-global combustion kinetics model consists of several global steps for parent fuel decomposition, soot formation, and soot oxidation and a detailed wet-CO mechanism to complete the combustion process. The final thermophysics formulations are incorporated with a computational fluid dynamics model for prediction of the combustion efficiency of an unielement, tripropellant combustor and the radiation of a kerosene-fueled thruster plume. The model predictions agreed reasonably well with those of the tests.
The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.
The spread of flames over the surface of solid combustible materials is affected by the environmental conditions. Variables such as flow velocity, oxygen concentration, ambient pressure, partial gravity or microgravity, may change the material flammability and influence the fire dynamics. This is an important fire safety issue for space exploration vehicles and space habitats which will very likely have different environments than those currently used on the International Space Station, such as reduced pressure and enriched oxygen concentration (SEA). However, testing experimentally the materials to be used and qualified for space exploration under these conditions is a difficult and expensive task. Here, numerical modeling is used to understand the dominant physico-chemical processes on the concurrent flame spread over thin fabrics under reduced ambient pressure (and in turn, buoyancy) and variable gravitational conditions. The numerical model uses the Fire Dynamics Simulator (FDS6) code with a single-step Arrhenius reaction rate for the solid phase decomposition. Different models are tested for the gas phase combustion kinetics. The model results are validated with experimental results obtained at similar reduced ambient pressure and flow conditions at 1 g. It is shown that as ambient pressure is reduced the flame spread rate over a thin fabric is also reduced, both experimentally and numerically. Numerical results are compared to an analytical approach previously developed to explain the experimental trends. Further interpretation of the model results provides information regarding the physics of the process and how they are affected by the lower pressure environments and gravity conditions. The results of this work provide guidance for potential on-earth testing for fire safety design in spacecraft and space habitats.
The general assumption that systems are in equilibrium is convenient, but imperfect. Instead, the partial or full expression of kinetic reaction steps may need to be determined. Using an isotope fractionation Monte-Carlo model under the harmonic approximation (Figure 1 envelopes), we find the range of kinetic isotope effects (i.e., resulting from unidirectional chemical reactions) contains and exceeds the range of accessible equilibrium isotope effects. For the 16O-17O-18O system, accessible kinetic isotope effects can yield Δ′17O variations as large 0.8‰ from a single kinetic fractionation step. Reactions that yield the largest deviations from equilibrium are those with the largest molecular mass difference between decomposition fragments, typically the two products, such as a gas and solid generated from the thermal decomposition of an initial solid. For example, previously published experiments on the thermal decomposition of calcite (CaCO3) and dehydroxylation of brucite (Mg(OH)2) are matched by density functional theory models of each kinetic isotope effect, thus revealing their origin in mass-dependent reaction. Beyond oxygen, our approach applies to any system with more than two isotopes.