Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature reactivity coefficient”

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 109 records · Page 6

Ion correlations explain kinetic selectivity in diffusion-limited solid-state synthesis reactions

Establishing viable solid-state synthesis pathways for novel inorganic materials remains a major challenge in materials science. Previous pathway design methods using pairwise reaction approaches have navigated the thermodynamic landscape with first-principles data but lack kinetic information, limiting their effectiveness. This gap leads to suboptimal precursor selection and predictions, especially for reactions forming competing phases with similar formation energies, where ion diffusion is a critical influence. Here we demonstrate an inorganic synthesis framework by incorporating machine learning-derived transport properties through ‘liquid-like’ product layers into a thermodynamic cellular reaction model. In the Ba–Ti–O system, known for its competitive polymorphism, we obtain accurate predictions of phase formation with varying BaO:TiO2 ratios as a function of time and temperature. We find that diffusion–thermodynamics interplay governs phase compositions, with cross-ion transport coefficients critical for predicting diffusion-limited selectivity. This work bridges length scales and timescales by integrating solid-state reaction kinetics with first-principles thermodynamics and spatial reactivity.

Atomistic models↗

NCSP Integral Experiments at Sandia in FY21 [Slides]

This presentation provides details regarding integral experiments at Sandia National Laboratory for fiscal year 2021. The experiments discussed are as follows: IER 230: Characterize the Thermal Capabilities of the 7uPCX; IER 304: Temperature Dependent Critical Benchmarks; IER 305: Critical Experiments with UO 2 Rods and Molybdenum Foils; IER 306: Critical Experiments with UO 2 Rods and Rhodium Foils ; IER 441: Epithermal HEX Lattices with SNL 7uPCX Fuel for Testing Nuclear Data; IER 452: Inversion Point of the Isothermal Reactivity Coefficient; and IER 523: Critical Experiments with ACRR UO 2 -BeO Fuel.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗

Chemical kinetic interactions and sensitivity analyses for 2-ethylhexyl nitrate-doped PRF91 using a reduced mechanism

A numerical and experimental investigation about the chemical kinetic interactions between 2-ethylhexylnitrate (EHN) and PRF91 was performed in this study. Rapid compression machine experiments were conducted to investigate the effect of EHN on the autoignition reactivity of the fuel, and a reduced chemical kinetic mechanism was developed including an EHN sub-model. Experiments showed that the ignition delay decreases as the fuel is doped with EHN, but the effect of the doping level of EHN on the ignition is highly non-linear. Moreover, experiments showed that the EHN effectiveness is lowest during the transition between the low-temperature regime and the negative temperature coefficient (NTC) regime, and it rapidly increases as the temperature increases. Both detailed and (developed) reduced mechanisms were validated against the experimental results, allowing a more in-depth EHN-fuel chemistry study. Additionally, ignition delay sensitivity and EHN effectiveness sensitivity analyses were performed with the reduced mechanism to identify the reactions that control the interaction between EHN and the fuel. As the result, EHN thermal decomposition is only relevant for very low temperatures. The chemistry of EHN-doped fuel is more sensitive to intermediate temperature reactions than that of straight fuel, especially at lower temperatures, due to the effect of EHN on the NTC behavior of the fuel. Finally, the chemistry of EHN-doped fuel is very sensitive to NO 2 -to-NO reactions, especially at high temperatures, because these reactions transform the NO 2 generated by EHN into NO, which is a very effective fuel reactivity enhancer.

2-ethylhexyl nitrate↗

Fiber-Reinforced Reactive Nano-Epoxy Composites

An ultra-high-molecular-weight polyethylene/ matrix interface based on the fabrication of a reactive nano-epoxy matrix with lower surface energy has been improved. Enhanced mechanical properties versus pure epoxy on a three-point bend test include: strength (25 percent), modulus (20 percent), and toughness (30 percent). Increased thermal properties include higher Tg (glass transition temperature) and stable CTE (coefficient of thermal expansion). Improved processability for manufacturing composites includes faster wetting rates on macro-fiber surfaces, lower viscosity, better resin infusion rates, and improved rheological properties. Improved interfacial adhesion properties with Spectra fibers by pullout tests include initial debonding force of 35 percent, a maximum pullout force of 25 percent, and energy to debond at 65 percent. Improved mechanical properties of Spectra fiber composites (tensile) aging resistance properties include hygrothermal effects. With this innovation, high-performance composites have been created, including carbon fibers/nano-epoxy, glass fibers/nano-epoxy, aramid fibers/ nano-epoxy, and ultra-high-molecularweight polyethylene fiber (UHMWPE).

Zhong, Wei-Hong↗

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Electrodeposition of Magnonic V(tetracyanoethylene) 2 Thin Films

Molecule-based magnetic materials have been identified as promising candidates for application in magnonic technologies, owing not only to their solution processability but also because they can exhibit narrow ferromagnetic resonance (FMR) linewidths and low Gilbert damping coefficients—crucial prerequisites for the transmission of coherent magnons over macroscopic distances. In particular, V(TCNE) 2 , a compound with a threedimensional network structure composed of vanadium(II) centers linked by tetracyanoethylene (TCNE •− ) radical anions, displays magnonic properties comparable to yttrium iron garnet, the quintessential magnonic material in the field. However, existing solution and chemical vapor deposition methods for synthesizing V(TCNE) 2 require the use of highly reactive zero-valent molecular vanadium precursors, stymying research on this important material. Herein, we report a facile electrochemical method for the deposition of thin films of V(TCNE) 2 using readily obtainable and stable divalent vanadium precursors and TCNE •− anions generated by electrochemical reduction. Magnetization measurements reveal that the films exhibit ferrimagnetic ordering above room temperature, consistent with V(TCNE) 2 films synthesized via other methods. Moreover, the electrodeposited films exhibit narrow FMR linewidths as low as 17.5 G and a low Gilbert damping coefficient of 1.1 × 10 −3 , values that are on par with some currently integrated metallic magnonic materials. More generally, these results demonstrate that electrodeposition can provide a straightforward means of generating highperformance magnonic materials using readily available molecular precursors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tribological and microstructural comparison of HIPped PM212 and PM212/Au self-lubricating composites

The feasibility of replacing the silver with the volumetric equivalent of gold in the chromium carbide-based self-lubricating composite PM212 (70 wt. percent NiCo-Cr3C2, 15 percent BaF2/CaF2 eutectic) was studied. The new composite, PM212/Au has the following composition: 62 wt. percent NiCo-Cr3C2, 25 percent Au, 13 percent BaF2/CaF2 eutectic. The silver was replaced with gold to minimize the potential reactivity of the composite with possible environmental contaminants such as sulfur. The composites were fabricated by hot isostatic pressing (HIPping) and machined into pin specimens. The pins were slid against nickel-based superalloy disks. Sliding velocities ranged from 0.27 to 10.0 m/s and temperatures from 25 to 900 C. Friction coefficients ranged from 0.25 to 0.40 and wear factors for the pin and disk were typically low 10(exp -5) cu mm/N-m. HIPped PM212 measured fully dense, whereas PM212/Au had 15 percent residual porosity. Examination of the microstructures with optical and scanning electron microscopy revealed the presence of pores in PM212/Au that were not present in PM212. Though the exact reason for the residual porosity in PM212/Au was not determined, it may be due to particle morphology differences between the gold and silver and their effect on powder metallurgy processing.

Bogdanski, Michael S.↗

Tribological and microstructural comparison of HIPped PM212 and PM212/Au self-lubricating composites

The feasibility of replacing the silver with the volumetric equivalent of gold in the chromium carbide-based self-lubricating composite PM212 (70 wt percent NiCo-Cr3C2, 15 percent BaF2/CaF2 eutectic) was studied. The new composite, PM212/Au has the following composition: 62 wt percent NiCo-Cr3C2, 25 percent Au, 13 percent BaF2/CaF2 eutectic. The silver was replaced with gold to minimize the potential reactivity of the composite with possible environmental contaminants such as sulfur. The composites were fabricated by hot isostatic pressing (HIPping) and machined into pin specimens. The pins were slid against nickel-based superalloy disks. Sliding velocities ranged from 0.27 to 10.0 m/s and temperatures from 25 to 900 C. Frictions coefficients ranged from 0.25 to 0.40 and wear factors for the pin and disk were typically low 10(exp -5) cu mm/N-m. HIPped PM212 measured fully dense, whereas PM212/Au had 15 percent residual porosity. Examination of the microstructures with optical and scanning electron microscopy revealed the presence of pores in PM212/Au that were not present in PM212. Though the exact reason for the residual porosity PM212/Au was not determined, it may be due to practice morphology differences between the gold and silver and their effect on powder metallurgy processing.

Bogdanski, Michael S.↗

Pulse radiolysis and transient absorption spectra of aqueous solutions of sodium sulfamate

Chemical kinetics for the reactions of sulfamate ions (NH 2 SO 3 − ) with the primary products of water and nitric acid radiolysis were measured in aqueous solutions at ambient temperature. Using time-resolved electron pulse radiolysis techniques with a custom multichannel detection system, we examined the reactivity of NH 2 SO 3 − with the hydroxyl radical ( • OH), hydrogen atom (H • ), and nitrate radical (NO 3 • ). The sulfamate ion was found to react with • OH and H • with second-order rate coefficients of k • OH = (5.60 ± 0.04) × 10 6 M −1 s −1 and k H • = (7.96 ± 0.10) × 10 6 M −1 s −1 , respectively, and with NO 3 • with a rate coefficient of k NO 3 • = (1.67 ± 0.06) × 10 7 M −1 s −1 . The reactions of NH 2 SO 3 − with • OH and H • resulted in the formation of two transient radical species, one with maximum absorbance at 300 nm and a second with maxima at both 300 nm and 600 nm. These spectra are tentatively assigned to • NH 2 SO 3 and • NHSO 3 − , respectively. By measuring the absorbance of these radicals as a function of pH, the radical pK a was determined to be 9.5 ± 0.1. Overall, this work has implications for the longevity and performance of ferrous sulfamate, Fe(NH 2 SO 3 ) 2 , as a plutonium reductant in the reprocessing of used nuclear fuel.

Conrad, Jacy K. [Idaho National Laboratory (INL), ↗

Odd nitrogen in the mesosphere and thermosphere - Its chemistry and transport

The atomic nitrogen in the thermosphere, which is produced either in its ground term or in a highly reactive metastable term, reacts with molecular oxygen to form nitric oxide and atomic oxygen. However, nitric oxide is destroyed in reactions with atomic nitrogen in which molecular nitrogen and atomic oxygen is formed. Adopted reactions and rate coefficients for odd nitrogen chemistry are listed in a table. An analysis is conducted of the temperature distribution below 150 km and its effect on the odd nitrogen chemistry. The latitudinal distribution of nitric oxide is discussed. Satellite measurements have shown that nitric oxide concentrations at high latitudes are highly variable in both time and space. Their average concentrations are 3-4 times higher than at mid-latitudes. The measured variation of nitric oxide density as a function of latitude is shown in a graph.

Rusch, D. W.↗

Solvent selection for biphasic extraction of 5-hydroxymethylfurfural via multiscale modeling and experiments

We introduce a comprehensive conceptual framework for selecting solvents for reactive extraction in biphasic organic-water systems and demonstrate it for the separation of HMF (5-hydroxylmethylfurfural), a platform chemical produced in the acid-catalyzed dehydration of hexoses. We first perform in silico screening of ~2500 solvents, from the ADFCRS-2018 database using the ADF COSMO-RS implementation, and classification, based on the solvent partition coefficient. We then determine experimentally the partition coefficients for HMF, fructose, and products of HMF rehydration (levulinic acid (LA), and formic acid (FA)), the mutual water-organic solvent solubilities, and the separation factors in >50 select solvents spanning multiple homologous series at room temperature and a typical reaction temperature with in situ sampling. We find that COSMO-RS is excellent for screening purposes (typical error in most cases within a factor of ~2). Increased temperatures lead to significant reduction in partitioning, and room temperature measurements are clearly inadequate for solvent selection. Upon down selecting classes of solvents based on separation performance, we perform experimental thermal stability and reaction compatibility studies of a small set of solvents at relevant reactive-extraction temperatures. We discover that many substituted phenols exhibit an order-of-magnitude increase in partitioning compared to conventional solvents due chiefly to hydrogen bond interactions and show the necessary stability but retain a significant fraction of water and LA, factors that need to be considered in technoeconomic analysis. In contrast, anilines, aldehydes, and acids are good to excellent regarding separation but incompatible with this specific reaction medium. Furthermore, this multifaceted framework can be extended to other biomass-derived products and processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fitting Vibrational Energy Parameters to Falloff Experiments: CH3 + CH3 + He (0.6-2 Torr 200-298 K)

There are many reported experiments and theoretical calculations for the combination of methyl free radicals diluted in a non-reactive 'deactivator. At high pressure the ethane (C2H6) combination product is the sole product; with decreasing pressure the chemically activated C2H6 will decompose. The falloff in the observed rate coefficient is the result of the competition between collisional stabilization of the chemically activated CA and unimolecular decomposition. The dependence of the rate coefficient on pressure, temperature and collision properties is complex and can not be calculated from first principles. The understanding of this system is not only of fundamental importance but is relevant to the recent detection of methyl free radicals in the atmospheres of Saturn - and Neptune. The temperatures of these outer planet atmospheres are in the 140-200 K region with total pressures (predominately H2 and He) less than 0.2 Torr. Experimentally determined rate coefficients have been reported for this reaction at T = 296-906 K and T = 200-408 K mostly with argon as the deactivator. At T = 200 K only the high pressure rate coefficient has been determined. Complete falloff curves over a wide temperature range (200-1600 K) with a variety of weak collider models used to simulate argon as the deactivator have also been reported by Klippenstein and Harding (KH). More recently we have reported the experimental rate coefficients in the falloff region with helium as the deactivator at 200 and 298 K. In this paper we have used the calculated falloff curves reported by KH for argon to determine the average energy transferred per collision for helium in our recently reported experiments. Collision rates were converted using Lennard Jones parameters; the temperature dependence of this conversion factor is noted. The helium experiments were consistent with a down of approximately 100 cm (exp-1); the temperature dependence was slight. The magnitude of down and its temperature dependence will be compared to other systems with similar substrate complexity and reaction energetics.

Tardy, Dwight C.↗

TOUGHREACT V4.12

TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.

Spycher, Nicolas↗

Impact of carbamate formation on the surface tension of epoxy‐amine curing systems

Abstract Epoxy‐amine reactive systems, including their surface tension, have been studied due to their use and importance in composites and adhesives in laboratory and industrial applications. Understanding how the system's surface tension changes during curing could help target applications for tailored composites. Using the Wilhelmy plate method, surface tensions for epoxy‐curing agent systems were investigated under different environments. While the epoxy exhibited an expected negative temperature coefficient of −0.096 mN/m/°C, the curing agent exhibited a positive temperature coefficient of 0.073 mN/m/°C, which is due to the curing agent reacting with atmospheric CO 2 to form carbamate through a zwitterion pathway and functional groups accumulating at the surface of the reacting systems. Carbamate formation results in a surface tension below that of the individual constituents in the epoxy‐amine curing agent system. Fourier transform infrared (FTIR) spectroscopy was used to verify in situ formation of carbamate, and by reducing the system's temperature, it was demonstrated that the system surface tension profiles displayed micellar surfactant‐like characteristics. FTIR was used to explore the difference in composition between the surface and the bulk, in particular carbamate formation at the surface, resulting in a lowering of the surface tension of reacting systems.

Shenk, Timothy M.↗

Kinetics of Radiation-induced Cr( II ) and Cr( III ) Redox Chemistry in Molten LiCl-KCl Eutectic

Chromium (Cr) is a frequent constituent of the metal alloys proposed for molten salt nuclear reactor (MSR) applications, and is typically the least noble metal ion present. Consequently, chromium is preferentially corroded into molten salt solutions. Here, the redox poise and redox cycling of chromium ions in the salt can greatly influence its corrosivity towards structural alloys, ultimately impacting the longevity of MSR systems. Radiation-induced chemistry is expected to play a significant role in determining the chromium oxidation state distribution during MSR operations. In the present research, electron pulse radiolysis techniques were employed to characterize the reactivity of Cr(II) and Cr(III) ions with primary radiolysis products in molten lithium chloride–potassium chloride (LiCl–KCl) eutectic over a temperature range of 400–600 °C. Both chromium oxidation states were found to rapidly react with the primary products of molten chloride salt radiolysis, i.e., the solvated electron (e S - ) and the dichlorine radical anion (Cl 2 ˙ - ). For reactions with the e S - , second-order rate coefficients (k) of k = (4.1 ± 0.2) and (6.1 ± 0.3) × 10 10 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively, were determined. Temperature-dependent measurements allowed for the derivation of activation parameters for electron capture by Cr(II) and Cr(III). Both chromium ions also react with Cl 2 ˙ - , k = (7.2 ± 0.3) and (1.4 ± 0.1) × 10 9 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

Ignition delay times of NH 3 /DME blends at high pressure and low DME fraction: RCM experiments and simulations

Here, autoignition delay times of ammonia/dimethyl ether (NH 3 /DME) mixtures were measured in a rapid compression machine with DME fractions of 0, 2 and 5 and 100% in the fuel. The measurements were performed at equivalence ratios $\varphi$= 0.5, 1.0 and 2.0 and pressures in the range 10–70 bar; depending on the fuel composition, the temperatures after compression varied from 610 K to 1180 K. Admixture of DME is seen to have a dramatic effect on the ignition delay time, effectively shifting the curves of ignition delay vs. temperature to lower temperatures, up to ~250 K compared to pure ammonia. Two-stage ignition is observed at $\varphi$= 1.0 and 2.0 with 2% and 5% DME in the fuel, despite the pressure being higher than that at which pure DME shows two-stage ignition. At $\varphi$= 0.5, a reproducible pre-ignition pressure rise is observed for both DME fractions, which is not observed in the pure fuel components. A novel NH 3 /DME mechanism was developed, including modifications in the NH 3 subset and addition of the NH 2 + CH 3 OCH 3 reaction, with rate coefficients calculated from ab initio theory. Simulations faithfully reproduce the observed pre- ignition pressure rise. While the mechanism also exhibits two-stage ignition for NH 3 /DME mixtures, both qualitative and quantitative improvement is recommended. The overall ignition delay times for ammonia/DME mixtures are predicted well, generally being within 50% of the experimental values, although reduced performance is observed for pure ammonia at $\varphi$= 2.0. Simulating the ignition process, we observe that the DME is oxidized much more rapidly than ammonia. Analysis of the mechanism indicates that this ‘early DME oxidation’ generates reactive species that initiate the oxidation of ammonia, which in turn begins heat release that raises the temperature and accelerates the oxidation process towards ignition. The reaction path analysis shows that the low-temperature chain-branching reactions of DME are important in the early oxidation of the fuel, while the sensitivity analysis indicates that several reactions in the oxidation of DME, including cross reactions between DME and NH 3 species presented here, are critical to ignition, even at fractions of 2% DME in the fuel.

10 SYNTHETIC FUELS↗