Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Combustion modeling”

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.

At least 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Core-Shell Oxidative Aromatization Catalysts for Single Step Liquefaction of Distributed Shale Gas (Final Technical Report)

The objective of this project was to design and demonstrate a core-shell structured multifunctional catalyst to convert the light (dry) components of shale gas into liquid aromatic compounds (primarily benzene and toluene) in a single step. Operated in a modular oxidative aromatization system (OAS) under a cyclic redox scheme, the novel catalyst and process can significantly improve the value and transportability of distributed shale gas. Since the project started, each quarter addressed a different set of tasks related to the completion of the milestone detailed in the project award. The yearly summaries of these tasks are summarized below: Q1-Q4: • Conducted project planning and literature search. • Investigated a number of SHC redox catalysts using thermogravimetric analysis and fixed-bed reactor experiments. • Initiated process modeling towards generating two process models for the methane DHA base case and OAS process. • Developed DHA catalysts capable of producing >500 g/kg-cat-hr aromatics at 80% or greater aromatics selectivity at 700°C. Q5-Q8: • Developed alternative approaches with sequential bed configurations to enhance the aromatic yields based on OCM+DHA • Improved the zeolite synthesis efficiency by using the microwave-assisted technique and investigated the synthesis conditions on the zeolite yield, crystalline structure and morphology • Constructed a set of Aspen Plus process models with significant energy savings for OAS as compared to the base case non-oxidative DHA. • Adapted conventional hydrothermal method to be applicable to the microwave synthesizer unit for more efficient catalyst synthesis. • Studied the structure of the OCM catalyst and the dispersion of the carbonate in the redox reactions and in methane flow with Raman Spectroscopy. Q9-Q12: • Scaled up the catalyst synthesis with the microwave synthesis method. Based on its performance, procedural characterizations and catalytic performance testing were further conducted for the new microwave synthesized catalysts with the newly-developed product analysis procedure. • Developed the reaction system setup for the C2-DHA or OCM+DHA reaction product and achieved a better product collection-analysis method for the aromatic products with an improved carbon balance. The product from the OCM reaction exhibited complicated effects on the DHA catalyst. • Conducted additional OCM catalyst characterization using Near Ambient Pressure X-ray Photoelectron Spectroscopy and in situ Raman characterization • Validated the significant energy savings for OAS as compared to the base case non-oxidative DHA. Successfully set up the simulation model for the OCM+DHA+SHC reaction system based on the updated experimental results from NCSU. Q13-End of project: • Synthesized new zeolite catalysts by the microwave method, conducted characterizations (XRD, SEM, and TEM) and catalytic behavior testing. • Explored the “wet” C 2 H 6 and C 2 H 4 DHA reactions with using steam co-feed. A subsequent reduction as the regeneration step can regenerate the DHA catalyst and recover 99% activity of the fresh performance. • Achieved a 15.3% single-pass aromatic yield from methane by rationally combining the OCM and DHA at different temperatures. • Conducted a 105-hour stability test with an improved regeneration procedure, with an average aromatic yield of 13.8%. • Developed new catalyst and achieved a record-high 23.2% yield.

03 NATURAL GAS↗
DOE OSTI2022

Matrimid substrates with bicontinuous surface and macrovoids in the bulk: A nearly ideal substrate for composite membranes in CO 2 capture

Scalable fabrication of thin-film composite (TFC) membrane for post-combustion carbon capture is often limited by the availability of high-performance nanoporous substrate. Ideally, the substrate should allow for fast gas transport at the selective layer/substrate interface as well as in the bulk of the substrate. In this study, highly permeable Matrimid substrates were prepared via vapor- and nonsolvent-induced phase separations. The phase separation mechanism of the Matrimid®/N-methyl-2-pyrrolidone casting solution was modulated by the addition of LiCl to control the solution thermodynamic stability, and nanoporous Matrimid substrates were formed with a bicontinuous surface and macrovoids in the bulk. This hierarchically optimized structure resulted in a CO 2 permeance of 2.60 × 10 5 GPU (1 GPU = 1 × 10 –6 cm 3 (STP) cm –2 s –1 cmHg –1 = 3.349 × 10 –10 mol m –2 s –1 Pa –1 ), which was ca. 11 times more permeable than a benchmark polyethersulfone substrate with cellular pores. The benefit of using this new substrate was demonstrated by coating a 170-nm amine-containing polymer to form a TFC facilitated transport membrane. The membrane exhibited a CO 2 permeance of 932 GPU at 57°C, which was 72 GPU higher than the counterpart coated on the benchmark substrate. Meanwhile, the CO 2 /N 2 selectivity was remained at 158. This permeance improvement could be well explained by the resistance-in-series model, where the improved permeance was attributed to the reductions in substrate and lateral diffusion mass transfer resistances. The upper bound (UB) analysis indicates that the Matrimid substrate is nearly ideal for the state-of-the-art polymers for CO 2 /N 2 separation. Here, the substrate improvement could also reduce the parasitic energy associated with the membrane process due to the reduced membrane cost and flue gas compression requirement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Heat transfer enhancement due to cold cap motion from bubbling in a waste glass melter

In this study, a computational fluid dynamics (CFD) model was developed to model the motion of a solid cold cap in a waste glass melter. Forced convection bubblers at the base of the melter release air into the molten glass, which forms large bubbles that travel upwards to the cold cap and augment heat transfer from the glass to the cold cap. The CFD model employs the Navier- Stokes equations to solve for the fluctuating flowfield using a rigid body motion dynamic fluid body interaction module. This allows for movement of the floating body in response to the bubbling forces calculated at each time step. Further, the heat flux delivered to the cold cap by the convective bubbling is studied as a function of the bubbling rate. Results for the moving cold cap are compared with the computed heat flux trends for a stationary cold cap. The heat flux delivered to the cold cap from the glass is 25% higher for the case with the moving cold cap. The heat flux was found to be proportional to v0.6 as opposed to v0.9 (where v is the normalized bubbling rate) for the stationary cold cap.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Industrial emission inventories of iron and phosphorus for past, present and future

Global emission inventories of aerosol iron by species, and of phosphorus, were developed for inclusion in Earth system modeling that includes atmospheric transport. The time period covered was 1850-2010 to evaluate human influence. A variety of approaches to compare measured and modeled values of aerosol iron was synthesized to estimate how well anthropogenic emissions are known. Anthropogenic emissions are constrained within a factor of two and if there is a bias, it is an overestimation. The exception is shipping emissions, which may be underestimated. Direct radiative forcing by iron aerosol is small in the global average (0.02 W m -2 ) but reaches 0.25 W m -2 around source regions. Regions of high direct radiative forcing are poorly constrained by observations. The central value of change in oceanic net primary productivity due to iron deposition is about 0.26 PgC yr -1 , also highly localized. This leads to a decrease in atmospheric CO2 accumulation that corresponds to a forcing of 0.008 W m -2 since 1850 (0.003-0.18 W m -2 ). Emissions of particulate matter to 2050 were estimated in response to different policy scenarios, including the need for renewable energy installations under decarbonization scenarios. Under the business-as-usual scenario, the total metal demand is 195 Mt/yr in 2020, peaking at 270 Mt/yr in years 2040-2045; under the Rapid Decarbonization scenario, total metal demand is nearly doubled at its peak: 480 Mt/yr in year 2040. Primary particulate matter emissions attributable to processing metals for renewable energy are also amplified under rapid decarbonization, but total atmospheric emissions are much lower because of the shift away from combustion.

58 GEOSCIENCES↗
DOE OSTI2022

Collaborative Research: Unravelling the Physics Associated with the Production of Extremely Dense Plasma States of Microscale (Final Report) Nanosecond-pulsed Discharges

The aim of this project is to study and establish the physical mechanisms that contribute to the formation of anomalously dense plasmas in high-pressure nanosecond-pulsed discharges. These discharges have a broad range of applications such as plasma-assisted combustion, plasma flow actuators, biomedical sterilization and exotic materials synthesis. The structure and formation of these discharges, producing high plasma densities of ~10 14 -10 15 cm -3 , are well-studied and understood. Fast-pulsed microscale high-pressure discharges can be driven to even higher densities of > 10 19 cm -3 , approaching warm dense matter conditions. The mechanisms that generate these plasmas have not been understood. Analysis of the warm dense matter state under laboratory conditions is an expensive and non-trivial endeavor. For instance, dense plasmas can be generated by electrical explosions of metal foils and wires. Plasmas generated after the explosion have a short lifetime and often present difficult conditions for diagnostics. Generation of dense plasmas was also achieved during high-voltage nanosecond pulsed discharges when the so-called explosive electron emission is obtained. Unfortunately, this process is very difficult to control for the studies of warm dense matter. In our recent study, we have shown that additional heating of plasma by lasers can further increase the density of plasma and even lead to the fully ionized state. This method, potentially, allows better control of the plasma parameters. In this work, we studied a second stage laser-heated micro-discharge using a self-consistent one-dimensional particle-in-cell Monte Carlo-collision (1D PIC-MCC) model coupled with Maxwell’s equations. We predicted the generation of a fully ionized plasma on the picosecond time scale. However, this model considered the plasma as an ideal gas despite the high pressure and the nearly fully ionized state. The ideal plasma model assumes that the dilute gas approximation is valid, where the inter-particle interactions are negligible. For charged particles this assumption holds as long as the shielded Coulomb potential assumption is valid. For very high plasma densities, this concept breaks down since the Debye sphere surrounding each charged particle no longer contains enough electrons to statistically provide the shielding of the single particle Coulomb interaction potential. At such densities, the plasma can no longer be described as ideal and non-ideal coupling effects need to be considered. In this report, we elucidate our recent work of developing a PIC-MCC model with improvements for non-ideal plasma conditions due to Coulomb coupling at high densities. In particular, we study the interaction of green light radiation and a dense microplasma, and explore the non-ideal plasma effects in this interaction. In this computational model, we implement the two most important non-ideal effects: ionization potential depression (IPD) and enhanced collision cross sections. Our primary goal is to study the physics associated with electromagnetic (EM) wave heating, also called the second-stage wave-heating, and establish the role of plasma non-ideality in this phenomenon. Our secondary goal is to improve the chemistry mechanism of the 1D PIC-MCC model by including a more detailed excited species collision treatment. At high pressures, stepwise ionization from excited species might play an important role in the ionization process. Previously, this ionization mechanism was neglected due to the excitation collision cross section of xenon being smaller than that of ionization. However, a preliminary study showed that the excited species density in the initial microplasma was an order of magnitude higher than the electron density. Therefore, my aim is to determine the significance of this additional ionization pathway to the plasma generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Zeolite Membrane Reactor for Pre-Combustion Carbon Dioxide Capture (Final Scientific/Technical Report)

Water-gas shift (WGS) reaction followed by carbon dioxide (CO 2 ) separation is a critical step in the integrated gasification combined cycle (IGCC) process for fossil-fuel-fired electrical power generation with CO 2 capture. To intensify the IGCC process hydrogen-permselective zeolite membrane reactor offers promise to replace the conventional energy-intensive fixed-bed reactors and solvent-based CO 2 capture units. The objectives of this project were to develop a bench-scale zeolite membrane reactor (total membrane area: 932 cm 2 for a 21-tube membrane bundle) for the water-gas-shift reaction of raw syngas from an oxygen-blown coal-gasifier for H 2 production with simultaneous CO 2 separation at the capability of about 2 kilograms H 2 per day (equivalent to 2 kW IGCC power plant) and to demonstrate significant progress toward achieving overall performance goal of 90% CO 2 capture rate with 95% CO 2 purity at the cost of electricity 30% less than the baseline carbon capture approaches. This report summarizes results obtained in this project on scaling up the zeolite membrane reactor by a factor of 200 in membrane area, tests of the bench-scale zeolite membrane reactor for the water-gas-shift reaction at high temperature and high-pressure, and techno-economic analysis of the integration of the zeolite membrane reactor in IGCC power plant for the electrical generation with CO 2 capture. With effective pore modification by catalytic cracking deposition of MDES (methyldiethoxysilane), fabrication of MFI-type zeolite membranes was successfully scaled-up from a lab-scale disk type to bench-scale multiple-tube bundles on low-cost alumina supports. A Co-Mo based sour shift catalyst was evaluated and used in the zeolite membrane reactors for water-gas-shift reaction. The reaction kinetic and gas-permeation equations were developed and employed in mathematic models to guide/predict experiments/performance of zeolite membrane reactors for water-gas shift reaction. Multiple-tube zeolite membrane bundles and reactors were designed, fabricated and tested for gas separation and water-gas-shift reactions with real raw syngas from a coal-fired gasifier operated at high temperature and pressure. The zeolite membrane reactors demonstrated good long-term thermal and chemical stability and constant H 2 permeance (>300 GPU) together with a Knudsen selectivity (~4.7) for H 2 over CO 2 in the field test (cumulative time >28 hours) with high-sulfur coal-derived syngas. The zeolite-membrane-reactor integrated IGCC process was designed using the performance experimentally measured by the University of Cincinnati team with a single-tube zeolite membrane reactor that offers a CO conversion >98% with more than 90% of CO 2 and H 2 captured. With the above integrative approaches, a techno-economic analysis for a cost-benefit comparison was performed to uncover features that determine the power output, capital expenditure, operating expenditure, cost of electricity and cost of CO 2 capture in a 550-MW zeolite-membrane-reactor integrated IGCC process. The integration of the zeolite membrane reactor in IGCC could provide a significant reduction of 80% and 27% in the power consumption for Selexol™ Acid Gas Removal and CO 2 compression, respectively, which lowers the total auxiliary power consumption by 12.5%. However, the low pressure required at permeate stream for maintaining the driving force of hydrogen permeation through the zeolite membrane costs a huge power in permeate compressor, compensating the power consumption reduction achieved with the membrane reactor. Thus, for coal-fired IGCC for electricity generation with 90% CO 2 captured, the integration of the membrane reactor could provide a CO conversion ~99% and a significant drop in cost-of-electricity using zeolite membrane with H 2 permeance >600 GPU and the H 2 /CO 2 selectivity over 70.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Numerical Investigation of a Central Fuel Property Hypothesis Under Boosted Spark-Ignition Conditions

In the present work, a central fuel property hypothesis (CFPH), which states that fuel properties are sufficient to provide an indication of a fuel’s performance irrespective of its chemical composition, was numerically investigated. In particular, the objective of the study was to determine whether Research Octane Number (RON) and Motor Octane Number (MON), as fuel properties, are sufficient to describe a fuel’s knock-limited performance under boosted spark-ignition (SI) conditions within the framework of CFPH. To this end, four TPRF-bioblendstock surrogates having different compositions but matched RON (=98) and MON (=90), were first generated using a non-linear regression model based on artificial neural network (ANN). Additionally, three unconventional bioblendstocks were included in the analysis: di-isobutylene (DIB), isobutanol, and Anisole. Skeletal reaction mechanisms were generated for the TPRF-DIB, TPRF-isobutanol, and TPRF-anisole blends from a detailed kinetic mechanism. Thereafter, numerical simulations were performed for the fuel surrogates using the skeletal mechanisms and a virtual cooperative fuel research (CFR) engine model, under a representative boosted operating condition. In the computational fluid dynamics (CFD) model, the G-equation approach was employed to track the turbulent flame front and the well-stirred reactor model combined with the multi-zone binning strategy was used to capture auto-ignition in the end-gas. In addition, laminar flame speed (LFS) was tabulated for each blend as a function of pressure, temperature, and equivalence ratio a priori, and the lookup tables were used to prescribe laminar flame speed as an input to the G-equation model. Parametric spark timing sweeps were performed for each fuel blend to determine the corresponding knock-limited spark advance (KLSA) and 50% burn point (CA50) at the respective KLSA timing. It was observed that despite same RON, MON, and engine operating conditions, the TPRF-anisole blend exhibited markedly different knock-limited performance from the other three blends. This deviation from the octane index (OI) expectation was shown to be caused by differences in laminar flame speed. However, it was found that relatively large fuel-specific differences in LFS (>20%) would have to be present to cause any appreciable deviation from the OI framework. Otherwise, RON and MON would still be robust enough to predict a fuel’s knock-limited performance.

42 ENGINEERING↗
DOE OSTI2021

Integrated modeling methodology for ash agglomeration in poly-disperse fluidized beds using particle population framework

This article discusses a unique agglomeration modeling methodology developed based on binary collisions to combine the effects of heterogeneity in ash chemistry and granular physics. A simple population balance is defined to find changes in the particle size distribution (PSD) of a fluidized bed. Thermodynamic equilibrium calculations and a computational fluid dynamics (CFD) code, are used to obtain hydrodynamic parameters. A method to calculate and use a distribution of collision frequencies in fluidized beds, corresponding to a poly-disperse particle size distribution, was developed in order to incorporate the particle level heterogeneities. A distribution of collision frequencies obtained for poly-dispersed particles showed a three orders of magnitude higher collision frequency amongst the smaller particles than the coarser ones, at the initiation of agglomeration. Ash agglomeration occurred when the slag amount (binder) was less than 10 wt% at temperatures <850 °C for Pittsburgh No. 8 coal.

42 ENGINEERING↗
DOE OSTI2025

Computational screening of fly ash zeolite sorbents for boric acid removal

In the United States, many impoundments at coal-fired power plants contain elevated contaminants like arsenic, boron, barium, and selenium. Zeolites synthesized from fly ash show promise as sorbents for these contaminants. However, optimizing sorption capacity is challenging due to numerous possible topologies, silicon to aluminum (Si/Al) ratios, and cation types. In this study, molecular simulations are used to design cationic zeolites for boric acid adsorption. Force field models based on quantum mechanical calculations (PBE + D2) for Na-, Ca-, Mn-, and Fe-exchanged chabazite and LTA are presented. The new D2FF force fields reproduce DFT energies with about half the error of UFF. Zeolite performance depends on Si/Al ratio and cation type, with low Si/Al ratio chabazite (CHA) and phillipsite (PHI) zeolite frameworks exchanged with Ca 2+ or Na + /Ca 2+ mixtures showing the highest adsorption. In conclusion, these findings suggest tailored fly ash-derived zeolites could provide effective boron removal from leachate ponds.

CCR impoundment↗
DOE OSTI2025

5G-TSN Architecture Capable of Providing Real-time Situational Awareness to Fossil-Energy (FE) Generation Systems (Final Technical Report)

This final report highlights the comprehensive achievements of the project focused on developing and validating a 5G-Time Sensitive Networking (TSN) architecture tailored for real-time operational awareness in fossil energy systems. The initiative successfully advanced through a series of technical milestones, including the integration of EMI-aware network models, deployment of advanced simulation frameworks, and real-world performance characterization at key sites such as UTEP and Fabens. Through the strategic use of NetSim® software, the team created and validated network configurations for wired and wireless environments, tested under varying congestion conditions, and verified network slicing implementations for URLLC-specific applications. Major accomplishments include the migration of simulation tools to the latest NetSim® version to support accurate modeling of TSN and network slicing, extensive EMI measurement campaigns, and the development of a robust simulation model for end-to-end SCADA system integration. Simulations compared both TDD and FDD duplexing modes, revealing insights into their performance under congested conditions. The wireless network was benchmarked for throughput, jitter, and delay metrics, aligning with 3GPP Release 15/16 and IEEE 802.1-TSN standards. A peer-reviewed conference paper was accepted and published, contributing to the broader academic and industrial discourse on 5G-TSN integration in energy systems, in addition to a journal article. Despite minor delays due to software limitations, the project achieved its objectives and delivered validated architecture ready for deployment in advanced energy network environments.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2021

Comparison of temperature adaptive calibration methods for laser induced fluorescence based fuel-in-oil instrument

Fuel dilution of engine oil (or fuel-in-oil, FiO) is an important issue as multiple and late-cycle fuel injection, integral to many combustion efficiency and emissions improvements (e.g., downsized boosted gasoline engines and catalyst thermal management) increases FiO rate. In addition to causing general wear and corrosion in engine due to decreased oil viscosity and pH buffering, FiO is also believed to cause destructive low speed pre-ignition (or super knock) in boosted SI engines. To understand the effects of engine operating conditions on the FiO rate, an optical diagnostic capable of measuring transient FiO on minute timescales was recently developed and demonstrated on a modified GM Ecotech engine system (Neupane et al., Applied Spectroscopy 2021). The measurement is based on adding a dye to the fuel and monitoring for its presence in oil via laser-induced fluorescence (LIF). The measured LIF signal is related to FiO concentration via pre-determined calibration factors using a multivariate classical least square (CLS) method.Since fluorescence quantum yield is a function of temperature, measured LIF intensity not only depends on FiO concentration but also oil temperature. To expand the applicability of the FiO diagnostic to transient oil-temperature conditions (e.g., cold start in practical engines), this study develops a method to account for oil-temperature variations. The effect of oil temperature (20°C - 95°C) on the LIF spectra of eight FiO samples ranging from ~0.8-15% was investigated. LIF intensities of the FiO samples decreased linearly with increasing temperature; the reductions being more significant at dye peaks. We develop a new calibration model (T-adaptive CLS) incorporating the temperature (T) effects on LIF intensity that enables simultaneous calculation of FiO and oil temperature. The improved FiO diagnostic with T-adaptive calibration is more robust, and applicable to varying oil-temperature conditions. For example, when strategies such as multiple/late fuel injections are applied to overcome cold-start instability due to use of low vaporization bio-based fuels such as ethanol, FiO rate is expected to be very high; the improved T-adaptive FiO diagnostic is hence relevant for engine- and fuel-system calibration and optimization. The diagnostic could also provide validation data for flow-field and spray interaction CFD models, further broadening the diagnostic’s utility for advancing engine technology and efficiency.

Neupane, Sneha↗
DOE OSTI2025

Preliminary Kinetic Analysis of Non-Equilibrium Plasma- Assisted Methanol Pyrolysis and Oxidation Experiments

Efforts to enhance power generation efficiency and reduce emissions have driven interest in novel combustion techniques, including non-equilibrium plasma (NEP) ignitors. NEP ignitors show promise in improving energy conversion efficiency, fuel reforming, emission control, and lean-flammability limits. However, their adoption is hindered by a limited understanding of the interplay between plasma-enhanced combustion and thermal chemistry, particularly for complex fuels under engine-relevant conditions. Developing experimentally validated kinetic mechanisms is therefore critical. Additionally, the increasing interest in renewable biofuels like ethanol and methanol, coupled with the desirable qualities of NEP ignitors, presents a compelling opportunity for study. Therefore, this work acts as an extension of a previous work (Bopaiah et al., 2023) pertaining to the experimental results of NEP-assisted methanol pyrolysis and oxidation. Experiments were performed with a custom-built plasma flow reactor at 0.5 atm and temperatures from 523-1203 K. All reactive mixtures are extremely diluted to minimize exothermicity due to reactivity, allowing isothermal assumptions and the isolation of plasma chemistry from thermal chemistry. A dielectric barrier discharge plasma, at 14 kV and 15 ns full-width half maximum, was applied to the reactive mixture at varying frequencies to maintain the number of pulses with increasing temperature. Steady-state product speciation was performed downstream of the reactor with ex-situ GC/MS diagnostics. The attained experimental results were examined through in-depth analysis performed by means of an in-development plasma-coupled kinetic mechanism. As discussed in the previous work, the plasma significantly accelerates methanol pyrolysis, increasing stable intermediate production, including oxygenated and nitrile species. Plasma-assisted oxidation shows even faster fuel consumption compared to pyrolysis and a 200 K ignition shift compared to thermal oxidation. For plasma-assisted pyrolysis, the model demonstrates that accelerated fuel consumption stems from dissociative quenching of excited N2 states with fuel and H2, generating H radicals that react to rapidly form CH3 and CH2OH radicals. At low temperatures, these radicals recombine to produce oxygenates, while CH3 drives nitrile and hydrocarbon formation at higher temperatures. While the model captures pyrolysis trends well, discrepancies in methane, ethylene, and ethanol predictions are present. Similarly, the model faces challenges in accurately representing plasma-assisted oxidation, predicting a much steeper fuel gradient and ignition 100 K earlier than the experiment. While a similar scheme to pyrolysis is nested in the reaction pathway, the enhancement of the O and H radical fluxes and their initiation of the OH and HO2 radical pools dominate fuel and intermediate oxidation. The overestimation of these processes is shown to be responsible for the divergence of model from experiment. While the modelling predictions of this preliminary mechanism are not perfect, they serve as a valuable starting point. Primarily, they elicited new reaction pathways that are not otherwise possible in thermal chemistry induced reaction kinetics. The results also provide a basis for the future work that should be performed. For example, theoretical and experimental studies on excited nitrogen species and fuel/fuel radical interactions, quantification of the NOx production, and the kinetics behind the slow ignition observed in oxidation should be emphasized.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗
DOE OSTI2020

Computational fluid dynamics analysis of char conversion in Sandia’s pressurized entrained flow reactor

Design and analysis of practical reactors utilizing solid feedstocks rely on reaction rate parameters that are typically generated in lab-scale reactors. Evaluation of the reaction rate information often relies on assumptions of uniform temperature, velocity, and species distributions in the reactor, in lieu of detailed measurements that provide local information. This assumption might be a source of substantial error, since reactor designs can impose significant inhomogeneities, leading to data misinterpretation. Spatially resolved reactor simulations help understand the key processes within the reactor and support the identification of severe variations of temperature, velocity, and species distributions. In this work, Sandia’s pressurized entrained flow reactor is modeled to identify inhomogeneities in the reaction zone. Tracer particles are tracked through the reactor to estimate the residence times and burnout ratio of introduced coal char particles in gasifying environments. The results reveal a complex mixing environment for the cool gas and particles entering the reactor along the centerline and the main high-speed hot gas reactor flow. Furthermore, the computational fluid dynamics (CFD) results show that flow asymmetries are introduced through the use of a horizontal gas pre-heating section that connects to the vertical reactor tube. Computed particle temperatures and residence times in the reactor differ substantially from the idealized plug flow conditions typically evoked in interpreting experimental measurements. Furthermore, experimental measurements and CFD analysis of heat flow through porous refractory insulation suggest that for the investigated conditions (1350 °C, <20 atm), the thermal conductivity of the insulation does not increase substantially with increasing pressure.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Diol isomer revealed as a source of methyl ketene from propionic acid unimolecular decomposition

Carboxylic acids are important combustion intermediates, especially regarding the combustion of oxygenates such as ethyl esters. The purpose of this study is to discern the unimolecular decomposition pathways of one such carboxylic acid, propionic acid, using microreactor flow experiments with line-tunable photoionization mass spectrometry and infrared spectroscopy, along with high-level electronic structure calculations (CCSD(T)/cc-pV∞Z//M06-2X/cc-pVTZ level of theory) and master equation theory. Here, microreactor experiments were performed at 300–1500 K, pressures decreasing from roughly 300 Torr to vacuum pressure, and around 100 μs residence times. Primary products are methyl ketene, ketene, ethylene, and methyl radical. Theory suggests the two lowest energy bond fissions are active at these conditions and responsible for the production of ketene, methyl radical, and some ethylene. Importantly, theory revealed the significance of the isomerization of propionic acid to propene-1,1-diol, and the subsequent dehydrogenation reaction of the diol as an alternative explanation for the unimolecular formation of methyl ketene. This is predicted to be the dominant thermal decomposition pathway at lower temperatures (up to ~850 K). Line-tunable VUV photons allowed for isomer resolution of the products and showed no evidence of the diol of propionic acid surviving until detection, suggesting propene-1,1-diol decomposes either to methyl ketene and water rapidly, fragments upon ionization to a distonic ion inseparable from methyl ketene through our detection methods, or never stabilizes as propionic acid well-skips directly to methyl ketene. Infrared spectroscopy showed no evidence of the decarboxylation of propionic acid at these conditions. Updated rate constants and branching ratios for propionic acid decomposition are calculated and provided for future modeling studies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Randomized Algorithms for Scientific Computing (RASC)

Randomized algorithms have propelled advances in artificial intelligence (AI) and represent a foundational research area in advancing AI for Science. Future advancements in DOE Office of Science priority areas such as climate science, astrophysics, fusion, advanced materials, combustion, and quantum computing all require randomized algorithms for surmounting challenges of complexity, robustness, and scalability. Advances in data collection and numerical simulation have changed the dynamics of scientific research and motivate the need for randomized algorithms. For instance, advances in imaging technologies such as X-ray ptychography, electron microscopy, electron energy loss spectroscopy, or adaptive optics lattice light-sheet microscopy collect hyperspectral imaging and scattering data in terabytes, at breakneck speed enabled by state-of-the-art detectors. The data collection is exceptionally fast compared with its analysis. Likewise, advances in high-performance architectures have made exascale computing a reality and changed the economies of scientific computing in the process. Floating-point operations that create data are essentially free in comparison with data movement. Thus far, most approaches have focused on creating faster hardware. Ironically, this faster hardware has exacerbated the problem by making data still easier to create. Under such an onslaught, scientists often resort to heuristic deterministic sampling schemes (e.g., low-precision arithmetic, sampling every nth element) and sacrifice potentially valuable accuracy. Dramatically better results can be achieved via randomized algorithms, reducing the data size as much as or more than naive deterministic subsampling can achieve, while retaining the high accuracy of computing on the full data set. By randomized algorithms we mean those algorithms that employ some form of randomness in internal algorithmic decisions to accelerate time to solution, increase scalability, or improve reliability. Examples include matrix sketching for solving large-scale least-squares problems (see Figure 1) and stochastic gradient descent for training machine learning models. We are not recommending heuristic methods but rather randomized algorithms that have certificates of correctness and probabilistic guarantees of optimality and near-optimality. Such approaches can be useful beyond acceleration, for example, in understanding how to avoid measure zero worst-case scenarios that plague methods such as QR matrix factorization.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Connected and Learning Based Optimal Freight Management for Efficiency

The management of the future heterogenous fleet is a complex decision-making problem. The heterogenous fleet is emerging as decarbonization technologies are deployed by fleets toward lowering the freight operation emissions in Medium and Heavy-duty vehicles. Traditionally, in fleets characterized by a homogeneous Diesel Internal Combustion Engine (ICE) powertrain, the process of fleet planning and operational optimization unfolds sequentially without the necessity to account for powertrain and vehicle-specific characteristics during dispatch decisions. Fleets with trucks less than 5 years old tend to maintain stable vehicle efficiency with minimal operational reliability risks for fleet managers. However, the landscape changes with the incorporation of emerging powertrain technologies, which lack extensive operational data and service experiences. This includes technologies like hybrid, Electric, Fuel Cell, or alternative fuel ICE. Operational decisions for fleets featuring heterogeneous powertrain technologies and facing limited access to alternative fueling and charging stations become intricate, requiring careful consideration and optimization at each dispatch. The difference in efficiency characteristics of emerging technologies, their range limitations, and the restricted availability of charging/alternative fueling infrastructure, coupled with sensitivity to driving conditions (e.g., EV range reduction in low temperatures) and their impact on component aging (such as batteries), become pivotal factors influencing the reliable and efficient freight transportation. To make the path toward low emission freight transportation efficient and reliable, an AI-assisted fleet management software is developed in this project to help fleet managers in optimizing both adoption of emerging powertrain decarbonization, connected and automated technologies and also operating the fleet after such technologies are deployed as schematically. Freight transportation requirements are different depending on the cargos to be shipped, customer requirements and regions of operations. This further highlights the need for software and digital solutions to tailor deployment and operation of emerging powertrain, connectivity, and automation technologies toward the specific fleet operation requirements. The fleet management optimizer was also integrated with a model of the fleet to simulate the operation of the fleet over 1 year of the baseline fleet operation (250,000+ shipments) indicating the significance of day-to-day variations on emissions and energy consumption of a freight transportation fleet. The results demonstrate ≥20% improvement in freight efficiency in terms of WTW CO2 per ton-mile of cargo shipments while all fleet operation constraints are enforced, and the cost (CapEx and OpEx) is minimized.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗