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.
This presentation introduces our research on using droplet burning data to validate combustion kinetic mechanisms of transportation fuels: challenges and opportunities.
Not Available
Supercritical carbon dioxide (sCO 2 ) cycles are being investigated for the future of power generation and will contribute to a carbon-neutral future to combat the effects of climate change. These direct-fired closed cycles will produce power without adding significant pollutants to the atmosphere. For such cycles to be efficient, they will need to operate at significantly higher pressures (e.g., 300 atm for Allam Cycle) than existing systems (typically less than 40 atm). There is limited knowledge on combustion at these high pressures and with a high dilution of carbon dioxide. Also, various experimental and computational investigations and model developments have been performed by the University of Central Florida (UCF), Embry-Riddle Aeronautical University (ERAU), and Stanford University (S.U.) to improve the current knowledge base and support the design and development of sCO 2 combustors. This project's technical aspects include chemical kinetics development, fundamental sCO 2 combustion, combustion model and sub-model development, supercritical fluid injection characterization, and heat transfer characterization.
One of the major biofuels used today is biodiesel, composed of fatty acid methyl esters. The fats and oils feedstocks used to make biodiesel are in relatively short supply. To meet increasing global demand, engineered microorganisms have been developed to catalyze conversion of sugar to fatty acid esters that contain unique ..beta..-hydroxy-esters. This study investigated the effect of a hydroxyl group on the combustion characteristics of methyl esters by comparing the chemical behavior of methyl hexanoate (MHx) and methyl 3-hydroxyhexanoate (M3OHHx) - used as surrogates for diesel-boiling-range esters with longer fatty acid chains. The oxidation of these esters was studied experimentally in a flow reactor at 0.84 and 10 bar, 600 to 1,100 K, and stoichiometric conditions; and in a constant volume combustion chamber at 5 and 10 bar, 600 to 900 K, and equivalence ratios of 0.3 and 0.6. MHx was more reactive in the constant volume chamber at temperatures below 800 K at 10 bar and equivalence ratio of 0.6. MHx also exhibited a higher indicated cetane number (16.4) than M3OHHx (8.1). We investigated this reactivity trend using an updated MHx kinetic model and M3OHHx model developed as part of this work. The kinetic models predicted that radicals formed from MHx were consumed by low-temperature chemistry reactions, whereas the M3OHHx reactivity was significantly governed by the ..beta..-radical (on the same carbon as the OH group) chemistry which mainly terminated via the less reactive chain propagation pathway to methyl-3-oxohexanoate + HO2. This study extends our understanding of structural effects on ester reactivity, which will allow for accurate surrogate formulation and performance simulations.
A perturbation analysis modeling approach based on a genetic algorithm was used to identify possible reaction pathways to explain previous experimental observations of the strong acceleration of syngas auto-ignition by trimethylsilanol (TMSO). Organosilicon reactions were taken from an existing chemical kinetic mechanism for tetramethylsilane pyrolysis which also contained TMSO reactions. The experimental targets for optimization were ignition delay times and peak OH radical concentrations in the range 1010–1070 K, 5 atm with dilute (4%) ϕ = 0.1 syngas mixtures in N 2 /Ar containing 0, 200, and 1000 ppm TMSO. Hundreds of thousands of trials resulted in three models which accurately predicted the experimental auto-ignition times and peak OH concentrations. Rate of production analyses indicated that TMSO forms dimethylsilanediol, which then reacts creating two “catalytic loops.” The net effect of the loops is to accelerate the reaction H 2 + O = OH + H at a faster effective rate than H 2 /O 2 kinetics alone, thus enhancing the reactivity of syngas. The formation of the loops may be attributable to the high Si–C bond energies in siloxanols, reducing rates of thermal decomposition and allowing dimethylsilanediol to react with H 2 for a longer duration compared with its alkane analog. Here, the results of the current work provide valuable direction for fundamental studies of these new hypothesized reaction pathways.
Not Available
In contrast to the adage “Models are to be used, not believed”, combustion kinetics models have been intended to be predictive in nature. Theoretical chemical kinetics is now understood to provide a firm foundation for the reaction parameters, thereby facilitating predictive simulations of chemical reactivity, even in regimes that are poorly characterized by chemical kinetic and/or combustion experiments. In this study, we describe a theory-informed kinetics model (ThInK) for small molecule combustion chemistry (H 2 and C 1 – C 3 species) that is based on the prodigious use of theoretical predictions for reaction rate coefficients, thermochemistry, and transport parameters. The distinct features of this kinetics model, which was developed over the course of several decades, are illustrated through simulations of flame propagation and auto-ignition.
Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.
This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS), as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Gamma Technologies, LLC (GT or Participante), to incorporate the ability to access LLNL chemical kinetics technologies while using GT-SUITE, GT’s market leading engine simulation software. At the end of the project, LLNL has released Zero-RK version 3.5 with zero- and one-dimensional (0-D and 1-D) solver functionality that interfaces with GT’s GT-SUITE v2023 and later releases. GT has tested its product to assure their customers that the interface can provide reduction in chemistry solution time for detailed chemistry simulations. The process has also positioned GT to easily benefit from future improvements of the Zero-RK suite of tools developed under the DOE Vehicle Technologies Office.
This was a collaborative effort between Lawrence Livermore National Security, LLC (LLNS), as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Gamma Technologies, LLC (GT or Participante), to incorporate the ability to access LLNL chemical kinetics technologies while using GT-SUITE, GT’s market leading engine simulation software. At the end of the project, LLNL has released Zero-RK version 3.5 with zero- and one-dimensional (0-D and 1-D) solver functionality that interfaces with GT’s GT-SUITE v2023 and later releases. GT has tested its product to assure their customers that the interface can provide reduction in chemistry solution time for detailed chemistry simulations. The process has also positioned GT to easily benefit from future improvements of the Zero-RK suite of tools developed under the DOE Vehicle Technologies Office.
A method for designing a combustion system which emits less of at least one environmentally-harmful emission is presented. In a describing step, an injector which introduces a fuel into a combustion chamber is described via a CFD code. In a modeling step, combustion kinetics of the fuel are modeled via a pre-processing code as the fuel mixes and reacts with an oxidizer. In a first selecting step, at least one primary scalar is derived during the modeling of the combustion kinetics. In a performing step, a table look-up is performed to obtain at least one data from a look-up database based on the primary scalar. In a second selecting step, at least one secondary scalar is selected in addition to the primary scalar(s). In a specifying step, at least one chemical pathway of formation or destruction for the secondary scalar is specified via a chemistry manager wherein the secondary scalar is representative of the environmentally-harmful emission(s) of the chemical pathway(s). In a utilizing step, the data is utilized to evaluate the chemical pathway(s) to quantify the environmentally-harmful emission(s). In an identifying step, an improvement to the combustion system is identified which reduces the environmentally-harmful emission(s).
Catalytic Fast Pyrolysis (CFP) is a process that converts biomass into liquid intermediates suitable for transportation fuels by rapidly heating it in the presence of a catalyst, aiming to produce stable oils with reduced oxygen content. During CFP, the catalyst can become deactivated by the accumulation of coke, a carbon-rich deposit formed from the decomposition of biomass components. Unlike in petroleum refining, regenerating coked catalysts from biomass pyrolysis requires specific approaches due to the different chemical nature of the coke formed. An experimental technique, Temperature Programmed Oxidation (TPO), was used to study the de-coking process by gradually increasing temperature while monitoring the production of CO and CO2, which provides data for kinetic modeling. Utilizing data from TPO experiments, coke combustion kinetic model was developed to describe the rate of coke removal at different temperatures, allowing for simulation of regeneration processes. Then kinetic model is integrated into MFIX-Exa for the simulation of Biogenic Fluid Catalytic Cracker (BFCC) regenerator at different scales, enabling analysis of catalyst flow, temperature distribution, and regeneration efficiency under various operating conditions.
Abstract Syngas is a desirable fuel for combustion in the Allam-Fetvedt cycle, which involves combustion under supercritical-CO2 conditions. While some work has been conducted in collecting ignition delay times (IDT) at the extreme pressures required by these systems, significant model deficiencies remain. Additionally, considerable barriers in terms of nonideal gas dynamic effects have been shown for these experiments in shock tubes. Further investigation into the fundamental combustion kinetics of H2/CO/CO2 mixtures is required. Time-resolved speciation measurements for target species have been shown to better aid in improving the understanding of underlying chemical kinetics than global ignition delay time measurements. Therefore, laser absorption measurements of CO were measured behind reflected shock waves during combustion of syngas at 5 and 10 bar and temperatures between 1080 and 2100 K. The mixtures investigated utilized H2-to-CO ratios of 1:1 and 1:4, respectively, each at stoichiometric conditions, allowing for discussions of the effect of initial fuel composition. A ratio of fuel to CO2 of 1:2 was also utilized to represent commercially available syngas. The mixtures were diluted in helium and argon (20% He, 76.5% Ar) to minimize thermal effects and to expedite CO thermal relaxation during the experiment. The resulting CO time histories were then compared to modern chemical kinetics mechanisms, and disagreement is seen for this system, which is assumed to be fairly well known. This study elucidates particular chemistry that needs improvement in moving toward a better understanding of syngas combustion at elevated pressures.
Oxycombustion in supercritical carbon dioxide is an integral part of the Allam Cycle, a technology that enables carbon-neutral use of fossil-fuels and carbon-negative use of biofuels. We simulate oxycombustion in a realistic combustor geometry for two sets of fuel conditions: pure methane, and a 40 percent methane/60 percent carbon dioxide blend, both at 343.15 K. The fuel jet mixes with a preheated swirler of 20 percent oxygen and 80 percent carbon dioxide at 1005.35 K, with a 100 percent carbon dioxide coflow at 783.15 K. The entire system operates at a pressure of 300 bar, putting the entire system above the critical temperature and pressure of carbon dioxide. Simulations are performed using PeleC, a compressible block-structured adaptive-mesh refinement (AMR) reacting flow code. Direct comparison of results using thermodynamically self-consistent implementations of the ideal gas equation of state and the Soave-Redlich-Kwong (SRK) equation of state allow quantification of the real gas impacts on combustion kinetics, flow temperatures, and pollutant formation.
Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.
The reaction of unsaturated compounds with ozone (O 3 ) is recognized to lead to the formation of Criegee intermediates (CIs), which play a key role in controlling the atmospheric budget of hydroxyl radicals and secondary organic aerosols. The reaction network of two CIs with different functionality, i.e. acetaldehyde oxide (CH 3 CHOO) and glyoxal oxide (CHOCHOO) formed in the ozone-assisted oxidation reaction of crotanaldehyde (CA), is investigated over a temperature range between 390 K and 840 K in an atmospheric pressure jet-stirred reactor (JSR) at a residence time of 1.3 s, stoichiometry of 0.5 with a mixture of 1% crotonaldehyde, 10% O 2 , at an fixed ozone concentration of 1000 ppm and 89% Ar dilution. Molecular-beam mass spectrometry in conjunction with single photon tunable synchrotron vacuum-ultraviolet (VUV) radiation is used to identify elusive intermediates by means of experimental photoionization energy scans and ab initio threshold energy calculations for isomer identification. Addition of ozone (1000 ppm) is observed to trigger the oxidation of CA already at 390 K, which is below the temperature where the oxidation reaction of CA was observed in the absence of ozone. The observed CA + O 3 product, C 4 H 6 O 4 , is found to be linked to a ketohydroperoxide (2-hydroperoxy-3-oxobutanal) resulting from the isomerization of the primary ozonide. Products corresponding to the CIs uni- and bi-molecular reactions were observed and identified. A network of CI reactions is identified in the temperature region below 600 K, characterized by CIs bimolecular reactions with species like aldehydes, i.e., formaldehyde, acetaldehyde, and crotonaldehyde and alkenes, i.e., ethene and propene. The region below 600 K is also characterized by the formation of important amounts of typical low-temperature oxidation products, such as hydrogen peroxide (H 2 O 2 ), methyl hydroperoxide (CH 3 OOH), and ethyl hydroperoxide (C 2 H 5 OOH). Detection of additional oxygenated species such as alcohols, ketene, and aldehydes are indicative of multiple active oxidation routes. This study provides important information about the initial step involved in the CIs assisted oligomerization reactions in complex reactive environments where CIs with different functionalities are reacting simultaneously. It provides new mechanistic insights into ozone-assisted oxidation reactions of unsaturated aldehydes, which is critical for the development of improved atmospheric and combustion kinetics models.
Plasma stability in reactive mixtures is critical for various applications from plasma-assisted combustion to gas conversion. To generate stable and uniform plasmas and control the transition towards filamentation, the underlying physics and chemistry need a further look. Here, this work investigates the plasma thermal-chemical instability triggered by dimethyl-ether (DME) low-temperature oxidation in a repetitive nanosecond pulsed dielectric barrier discharge. First, a plasma-combustion kinetic mechanism of DME/air is developed and validated using temperature and ignition delay time measurements in quasi-uniform plasmas. Then the multi-stage dynamics of thermal-chemical instability is experimentally explored: the DME/air discharge was initially uniform, then contracted to filaments, and finally became uniform again before ignition. By performing chemistry modeling and analyzing the local thermal balance, it is found that such nonlinear development of the thermal-chemical instability is controlled by the competition between plasma-enhanced low-temperature heat release and the increasing thermal diffusion at higher temperature. Further thermal-chemical mode analysis identifies the chemical origin of this instability as DME low-temperature chemistry. This work connects experiment measurements with theoretical analysis of plasma thermal-chemical instability and sheds light on future chemical control of the plasma uniformity.
Blending petroleum fuels with oxygenates is a common approach to stem the depletion of crude oil while also mitigating the impact of their combustion on the environment. It has recently been considered that n-butyl acetate (BA, C6H12O2 , boiling point of 399K) could be a viable oxygenate additive to diesel fuel. In this application it is important to determine the influence of the fractional amount of BA on burning of the mixture. This presentation considers this problem from an experimental and computational approach using n-heptane (C7H14 , boiling point of 372K) as an essentially single component surrogate for diesel fuel. The burning configuration considered was ostensibly spherical symmetry as promoted by burning droplets under conditions where forced and natural convection effects were minimized. In this configuration the droplet and flame are concentric spheres and soot aggregates are trapped in a shell structure between the droplet and flame. Measurements were made of the droplet, flame and soot shell diameters through video imaging of the burning process. The experimental results were used to validate a detailed numerical model of the spherically symmetric droplet burning process that incorporated a consistent combustion kinetic mechanism for BA/heptane mixtures, comprised of 402 species and 16,872 reactions inclusive of soot chemistry. A model for soot formation was also included in the simulation to predict formation of the soot shell. The simulations agreed well with measured droplet and flame diameters. It was found both experimentally and computationally that adding BA to heptane had a minimal effect on the mixture burning rate, while the flame was positioned closer to the droplet surface and the simulated soot volume fraction decreased. The implications of these results are discussed for using BA as a potential additive to diesel fuel.