Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coefficient of performance”

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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Space-Charge Aberrations in Single-Shot Time-Resolved Transmission Electron Microscopy

In this paper, we discuss effects of space-charge fields on imaging performance in a single-shot time-resolved transmission electron microscope. Using a Green’s function perturbation method, we derive analytical estimates for the effects of space-charge nonlinearity on the image formation process and the associated aberration coefficients. The dependence of these coefficients on the initial beam phase space distribution is elucidated. The results are validated by particle tracking simulations and provide fundamental scaling laws for the trade-off between temporal and spatial resolutions in single-shot time-resolved transmission electron microscopy.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Thermal WIMPs and the scale of new physics: global fits of Dirac dark matter effective field theories

We assess the status of a wide class of WIMP dark matter (DM) models in light of the latest experimental results using the global fitting framework GAMBIT. We perform a global analysis of effective field theory (EFT) operators describing the interactions between a gauge-singlet Dirac fermion and the Standard Model quarks, the gluons and the photon. In this bottom-up approach, we simultaneously vary the coefficients of 14 such operators up to dimension 7, along with the DM mass, the scale of new physics and several nuisance parameters. Our likelihood functions include the latest data from Planck, direct and indirect detection experiments, and the LHC. For DM masses below 100 GeV, we find that it is impossible to satisfy all constraints simultaneously while maintaining EFT validity at LHC energies. For new physics scales around 1 TeV, our results are influenced by several small excesses in the LHC data and depend on the prescription that we adopt to ensure EFT validity. Furthermore, we find large regions of viable parameter space where the EFT is valid and the relic density can be reproduced, implying that WIMPs can still account for the DM of the universe while being consistent with the latest data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Let’s Stick Together: Interplay of Moisture, Particle Size, and Anatomical Fraction in the Flowability of Corn Stover Crumbles

Handling of agricultural waste biomass feedstocks (such as corn stover) is a persistent challenge in the production of biofuels and chemicals at integrated biorefineries. Inconsistent flowability and high feedstock variability creates equipment downtime and processing challenges. This work investigates the influence of anatomical fractionation and moisture content on the bulk solids handling behavior of corn stover crumbles. A combination of rotational shear testing, tribological measurements, and surface and sorption characterization techniques (inverse gas chromatography (IGC), dynamic vapor sorption (DVS), and electron microscopy) were used to evaluate the contributions of particle composition, size, and surface properties to flow behavior. Shear and frictional testing revealed that moisture content has the greatest effect on flowability, significantly increasing unconfined yield strength (ƒ c ) and reducing the flow function coefficient (F F C), regardless of anatomical fraction. While differences in surface energy were observed between fractions, particularly in polar contributions (γ AB ), these did not correspond to meaningful differences in flow behavior. Coarser 6 mm particles exhibited poorer flow performance than 4 mm crumbles, likely due to broader particle size distributions and increased particle interlocking. Taken together, these results suggest that anatomical fractionation provides minimal benefit from a bulk handling perspective. Moisture content and particle size heterogeneity dominate handling behavior, reinforcing the utility of unfractionated corn stover in processing environments. These insights can inform the design of more efficient and reliable feedstock processing systems for agricultural biomass in integrated biorefineries.

09 BIOMASS FUELS↗
DOE OSTI2023

High Efficiency 2-Bladed Utility Wind Turbine Enhanced by CoFlow Jet Active Flow Control

The aerodynamic performance and power output of a high efficiency Co-Flow Jet (CFJ) 2-bladed Wind Turbine is studied in this paper. CFJ is a Zero-Net-Mass-Flux (ZNMF) active flow control method that dramatically increases airfoil lift coefficient and suppresses flow separation at a low energy expenditure. The 3D Reynolds Averaged Navier-Stokes (RANS) equations with one-equation Spalart-Allmaras (SA) turbulence model are solved to simulate the 3D flows of the wind turbines. The CFJ-Wind Turbine in this paper is modified from the NREL 3-bladed 3.4 MW reference wind turbine with the same RPM, blade length, distribution of chord, twist, and bending. The CFJ 2-bladed wind turbine utilizes the CFJ-NACA 64 series airfoils at most of the blade span from 31.88\% to the tip. The CFJ injection and suction slots are implemented along the full blade span to enhance the effectiveness and power output. The study shows that the 2-bladed CFJ turbine significantly increases power output compared with the 3-bladed baseline turbine at all wind speeds by 15.8\% at the wind speed of 9.8 m/s and by 25.9\% at 4 m/s respectively. The efficiency enhancement is substantially more at the low wind speed than at high wind speed. This would be very beneficial to increase the annual power output or capacity factor as wind turbines work at the speeds lower than the rated speed for most of the time. The results show that the flow field around the blade surface presents a significant enhanced suction peak near the blade leading edge, which results in a very high tangential force and thus increased wind turbine net power coefficient. Consistent with the previous study for the small NREL Phase VI research wind turbine, this paper demonstrates that CFJ 2-bladed wind turbine is an effective technology to extract more power from the wind at all speeds. This study is just an initial effort to demonstrate the advantage of the CFJ blade with no design optimization. It is believed that there is more potential to further enhance the efficiency improvement.

active flow control↗
DOE OSTI2024

Performance Improvements of the Griffin Solvers in FY24

The Griffin code is a MOOSE-based reactor physics application jointly developed by Idaho National Laboratory and Argonne National Laboratory under the Department of Energy Office of Nuclear Energy Nuclear Energy Advanced Modeling and Simulation Program. This fiscal year, we have made significant efforts to improve the performance of transport solver options and cross-section generation for the efficient use of Griffin in advanced reactor applications. For the HFEM-PN solver, the residual evaluations of HFEM kernels were optimized by utilizing the pre- computed averaged cross sections for individual elements. Numerical integration involving the evaluation of basis functions at quadrature points was bypassed by facilitating precomputed element mass matrices for response matrices. Red-black iterations were improved by introducing a new generalized minimum residual based solver. The memory usage of response matrix storage was significantly reduced by applying basis function rotations on interfaces and calculating volumetric odd-parity moments on the fly. Additionally, the adjoint flux and transient calculation capabilities of the HFEM-PN solver were successfully implemented and verified using the TWIGL benchmark problem. For the DFEM-SN solver, memory footprint and computation time were significantly reduced by not treating angular flux vectors as the MOOSE nonlinear system vectors. Specifically for IQS, scalar adjoint weighting was introduced to further eliminate angular adjoint flux storage in the MOOSE auxiliary system. It was demonstrated through the three-dimensional Advanced Burner Test Reactor core problem that the memory usage for transient calculations with the IQS method was reduced by over 7.5× compared to before the optimizations. For the self-shielding application programming interface, a new double-heterogeneity treatment method, named the Bell Function-Based Analytic Two-Region Slowing Down Method, was developed to efficiently flux-volume homogenize TRISO particles with the matrix. Additionally, optimizations were made to hyper- fine group (HFG) slowing down calculations by pretabulating collision probability coefficients and grouping isotopes, significantly reducing the computational time for calculating scattering sources per HFG. Lastly, the pin power reconstruction module was extended to account for temporal behavior in a microreactor analysis problem, specifically for a control drum transient. Verification tests for each of these improvements demonstrated significant performance enhancements and memory reduction.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

First-principles investigation of cerium and neodymium diffusion in BCC chromium and vanadium via vacancy-mediated transport

Lanthanide transport plays a crucial role in the performance and longevity of metallic nuclear fuels. This study examines the diffusion behavior of Ce and Nd—two major fission products—in body-centered cubic (BCC) Cr and V, which are potential liner or coating materials for mitigating fuel-cladding chemical interactions (FCCI). Using density functional theory (DFT) calculations and self-consistent mean-field (SCMF) analysis, the vacancy-mediated diffusion coefficients are evaluated. Our findings reveal that Ce and Nd act as oversized solutes and are strongly bound to vacancies in BCC Cr and V, with diffusivities in Cr significantly lower than in V and in hexagonal closed-packed (HCP) Zr, as investigated in our previous work. The activation energies for Ce and Nd diffusion are 3.39 and 3.32 eV, respectively, in BCC Cr, and 2.56 and 2.33 eV, respectively, in BCC V. Analysis of vacancy drag and partial diffusion coefficient ratios indicates a strong tendency for lanthanide enrichment at vacancy sinks in BCC Cr, and to a lesser extent in BCC V, with this effect persisting up to the melting point in Cr and remaining substantial for Nd in V at high temperatures. Under irradiation, the increase in vacancy concentration is expected to enhance lanthanide transport, potentially accelerating interactions at liner-cladding interfaces. Although BCC Cr exhibits relatively low lanthanide diffusivities under equilibrium conditions, the expected segregation tendencies under irradiation suggest that Zr liners may be a more favorable option. Further investigations using rate theory, cluster dynamics, and phase-field modeling are required to quantitatively assess the performance of these materials in reactor environments.

36 - MATERIALS SCIENCE↗
DOE OSTI2024

High Power Density Thermal Energy Storage With Phase Change Material in Enhanced Compact Heat Exchangers

Abstract Performance of a novel ultracompact thermal energy storage (TES) heat exchanger, designed as a microchannel finned-tube exchanger is presented. With water as the heating–cooling fluid in the microchannels, a salt hydrate phase change material (PCM), lithium nitrate trihydrate (LiNO3 · 3H2O), was encased on the fin side. To establish the hypothesis that small-length-scale encasement (<3 mm) of PCM substantially enhances heat transfer to yield very high power-density energy storage, heat exchanger designs with 10 and 24 fins/inch were considered. They were subjected to thermal cycling, or repeated heating (melting) and cooling (freezing), with inlet fluid flow mimicking diurnal variation between 42 °C and 25 °C (representing typical arid-region conditions) over an accelerated time period. By employing salt self-seeding to obviate subcooling during cooling or recrystallization, the TES was found to exhibit stable long-term (100 heating–cooling cycles) operation with very high PCM-side heat transfer coefficients (∼100–500 W/m2 K) and storage power density (∼160–175 kW/m3). In fact, with optimization of heating–cooling fluid flowrate for given charging–discharging time period and exchanger size, power density >300 kW/m3 can be achieved. The results clearly establish that highly compact heat exchangers used as TES units can provide very high-performance alternatives to conventional ones.

Engineering↗
DOE OSTI2023

The Jupiter High-240 Experiment

The Jupiter High-240 experiment was performed in May of 2019 by researchers from Los Alamos National Laboratory (LANL) at the National Criticality Experiments Research Center (NCERC) in the Device Assembly Facility (DAF) located at the Nevada National Security Site (NNSS). This experiment has previously been mentioned briefly in prior publication related to a collaborative effort with the Japan Atomic Energy Agency (JAEA) to assess lead void coefficients of reactivity in uranium- and plutonium-fueled systems with lead. This series of experiments supports JAEA’s research into the development of an accelerator driven transmutation system for spent nuclear fuel. The Jupiter High-240 experiment built upon the previous Jupiter experiment by incorporating plutonium fuel plates with higher 240 Pu content. Efforts to formally benchmark the original Jupiter experiment have continued for inclusion in the benchmark handbook of the International Criticality Safety Benchmark Evaluation Project (ICSBEP). Whereas there is much similarity between the two Jupiter experiments, there is a desire to also evaluate and benchmark this second experiment to further contribute towards the availability of lead-sensitive benchmarks. The components utilized in these two experiments have also been used to perform other subcritical and Rossi-α measurements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Temperature and illumination dependence of silicon heterojunction solar cells with a wide range of wafer resistivities

Abstract Recently, the significant improvements in the surface and contact passivation of silicon (Si) solar cells as well as their bulk quality have shifted their operating point to higher injections. Hence, they are less dependent on wafer doping. This shift opens an opportunity of using high‐resistivity wafers for practical photovoltaic applications, introducing a promising approach to push the cell efficiency towards the intrinsic limit and to improve the module reliability by increasing the cell breakdown voltage. Therefore, insights into the performance of Si solar cells using high‐resistivity wafers at various operating temperatures are of significant interest. In this study, we investigate the temperature‐ and illumination‐dependent performance of Si heterojunction (SHJ) solar cells using a wide range of wafer resistivities (between 3 and 1000 Ω⋅cm). Although a reduction in the passivation quality of the passivating contacts is observed at elevated temperature, the impact on the temperature coefficient of the open‐circuit voltage (TC V oc )—the dominant contributor to the temperature coefficient (TC) of the cell efficiency—is very limited. Their TC V oc are still dominated by the temperature dependence of the effective intrinsic carrier concentration. Furthermore, we also find that the investigated cells are more sensitive to temperature variation at lower illumination intensities. It is noteworthy that the efficiency of the cells fabricated using high‐resistivity wafers is comparable to that of the reference cells at any given temperature, highlighting the potential of using high‐resistivity wafers for solar cells.

14 SOLAR ENERGY↗
DOE OSTI2020

Modern Approaches to Exact Diagonalization and Selected Configuration Interaction with the Adaptive Sampling CI Method

Recent advances in selected configuration interaction methods have made them competitive with the most accurate techniques available and, hence, creating an increasingly powerful tool for solving quantum Hamiltonians. In this work, we build on recent advances from the adaptive sampling configuration interaction (ASCI) algorithm. We show that a useful paradigm for generating efficient selected CI/exact diagonalization algorithms is driven by fast sorting algorithms, much in the same way iterative diagonalization is based on the paradigm of matrix vector multiplication. We present several new algorithms for all parts of performing a selected CI, which includes new ASCI search, dynamic bit masking, fast orbital rotations, fast diagonal matrix elements, and residue arrays. The ASCI search algorithm can be used in several different modes, which includes an integral driven search and a coefficient driven search. The algorithms presented here are fast and scalable, and we find that because they are built on fast sorting algorithms they are more efficient than all other approaches we considered. After introducing these techniques, we present ASCI results applied to a large range of systems and basis sets to demonstrate the types of simulations that can be practically treated at the full-CI level with modern methods and hardware, presenting double- and triple-ζ benchmark data for the G1 data set. The largest of these calculations is Si$_2$H$_6$ which is a simulation of 34 electrons in 152 orbitals. We also present some preliminary results for fast deterministic perturbation theory simulations that use hash functions to maintain high efficiency for treating large basis sets.

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

High-pressure high-temperature synthesis and thermal equation of state of high-entropy transition metal boride

High-entropy transition metal boride (Hf 0.2 Ti 0.2 Zr 0.2 Ta 0.2 Mo 0.2 )B 2 sample was synthesized under high-pressure and high-temperature starting from ball-milled oxide precursors (HfO 2 , TiO 2 , ZrO 2 , Ta 2 O 5 , MoO 3 ) mixed with graphite and boron-carbide. Experiments were conducted in a large-volume Paris Edinburgh press combined with an in-situ Energy Dispersive X-ray Diffraction. The hexagonal AlB 2 phase with an ambient pressure volume V 0 = 27.93 ± 0.03 Å 3 was synthesized at a pressure of 0.9 GPa and temperatures above 1373K. High-pressure high-temperature studies on the synthesized high-entropy transition metal boride sample were performed to 7.6 GPa and 1873 K. The thermal equation of state fitted to the experimental data resulted in an ambient pressure bulk-modulus K 0 = 344 ± 39 GPa, dK/dT = –0.108 ± 0.027 GPa/K, and a temperature dependent volumetric thermal expansion coefficient α = α 0 + α 1 T + α 2 T –2 . The thermal stability combined with a high bulk-modulus establishes this high-entropy transition metal boride as an ultrahard high-temperature ceramic material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

High-resolution cryo-EM reconstructions in the presence of substantial aberrations

Here, an analysis is performed of how uncorrected antisymmetric aberrations, such as coma and trefoil, affect cryo-EM single-particle reconstruction (SPR) results, and an analytical formula quantifying information loss owing to their presence is inferred that explains why Fourier-shell coefficient-based statistics may report significantly overestimated resolution if these aberrations are not fully corrected. The analysis is validated with reference-based aberration refinement for two cryo-EM SPR data sets acquired with a 200 kV microscope in the presence of coma exceeding 40 µm, and 2.3 and 2.7 Å reconstructions for 144 and 173 kDa particles, respectively, were obtained. The results provide a description of an efficient approach for assessing information loss in cryo-EM SPR data acquired in the presence of higher order aberrations, and address inconsistent guidelines regarding the level of aberrations that is acceptable in cryo-EM SPR experiments.

3D reconstruction and image processing↗
DOE OSTI2024

Fourier Analysis and Loss Modeling for Inductive Wireless Electric Vehicle Charging with Reduced Stray Field

With the growth of electric vehicle (EV) popularity, different charging options to increase user convenience and reduce charging times are being considered and researched. Among these, inductive wireless power transfer (WPT) systems for EVs are being designed to meet specifications such as stray field, power level, efficiency, misalignment tolerance, and ground clearance, which are all heavily influenced by the coil geometry. The proposed Fourier Analysis Method (FAM) is an analytical method to directly design coil geometries to meet stray field and power level requirements through an optimization of Fourier basis function coefficients. The outputs of the optimization are complex, planar coil geometries that meet the power level and stray field constraints with minimized loss factors. Contours of these potentials determine the coil conductor paths and loss models predict the system efficiency and performance over misalignment. The Fourier representation of the geometry is used to conveniently calculate the coupling over misalignment, external proximity effect loss, and ferrite loss. A 6.6 kW prototype WPT system with low stray field and high efficiency is built from the optimization results to validate the models and showcase the usefulness of the FAM design approach.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗
DOE OSTI2022

Evaluation of Binders and Plasticizers in Kollidon VA 64-PEG Binder Systems

The binders, plasticizers, and dispersants in a polyvinylpyrrolidone/polyethylene glycol/glycerin binder system for PZT were evaluated. Kollidon VA 64 was investigated as a possible alternative binder to Kollidon 25 in a PZT powder system. The target amount of PEG300 in a Kollidon VA 64 system was predicted to be 15 to 30 wt.% PEG300 based on T g analysis by DSC. The compaction properties (slide coefficient, cohesiveness, green strength, etc.) were analyzed for Kollidon VA 64 – x PEG300 – glycerin systems. The properties in the range of x = 0 to 20 for systems without glycerin and x = 5 to 20 for systems with glycerin all exceeded the performance of the baseline Kollidon 25 system, of which VA 64 – 10 wt.% PEG300 – 5 wt.% glycerin with adsorbed moisture was the most promising composition due to a compact cohesiveness of 0.84 at 40 kpsi compared to a baseline of 0.44. The effect of dispersants on the compaction properties of a Kollidon 25 – PEG300 binder system was also analyzed, and the compaction properties were also compared to that of a Aquazol 200 – PEG6000 binder system. The powders with dispersant exhibited comparabl e per formance to the baseline, suggesting good compatibility. The compacts produce with the Aquazol 200 – PEG6000 binder exhibited decreased performance when compared to the baseline .

36 MATERIALS SCIENCE↗
DOE OSTI2021

ARC Software Validation Work for the FFTF Reactor

Extensive efforts have been carried out at ANL for the verification and validation of the Argonne Reactor Codes (ARC) software package currently used for the design of Versatile Test Reactor (VTR). The ARC software package consists of steady state neutronics and thermal hydraulics modeling capabilities which are being used by the VTR program to develop most of the VTR reactor design details which will be part of the licensing application. It is anticipated that this software will continue to be used for the design work and for initial operations although additional software may be introduced at a later time. The validation work was focused primarily on obtaining validation data consistent with VTR and usable for the ARC software. Because no critical facilities or operating fast spectrum reactors are available to do experiments for the VTR, the next best option is to identify historical experimental data that can be used as validation data. Early on in VTR, the ZPPR-15 set of experiments was identified as good validation data because of 1) the availability and quality of the data, 2) existing staff that are already familiar with the experimental machine and measurements, 3) most of the ZPPR-15 loadings of interest have already been processed into ARC models, and 4) a full uncertainty quantification has already been done for several loadings of ZPPR-15. The FFTF startup and operations data was identified as the most consistent reactor type that has validation data usable for VTR. Finally, the EBR-II fuel depletion measurements were identified as the best available validation data for VTR. It is important to note that both the FFTF and EBR-II reactors typically come with higher uncertainties than the ZPPR. In the frame of the discussed verification and validation efforts, the present document discusses the analysis of selected FFTF measurements included in the benchmark specifications of the International Reactor Physics Experiment (IRPhE) handbook. The FFTF reactor core configurations from the benchmark specification are presented in Section 2. The analysis is performed with the use of the ARC code suite available at ANL for fast reactor studies and is discussed in Section 3. The reactor parameters from the benchmark include criticality, neutron spectra, effective delayed neutron spectra, control rod worth, isothermal temperature coefficient and low energy gamma-ray spectra. The calculated values and the comparison with the experimental data is discussed in Sections 4 to 9 for each considered reactor parameter. Finally, conclusions are presented in Section 10.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

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↗