Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic parameters”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Chemical Preparation, Crystal Structure Reinvestigation and Vibrational Study of CoNa 3 P 3 O 10 ·12H 2 O and X-ray Characterization of the New Anhydrous Triphosphate CoNa 3 P 3 O 10

The triphosphate dodecahydrate of cobalt and sodium CoNa 3 P 3 O 10 ·12H 2 O was prepared by the ion exchange resin process. Its structure was studied by X-ray diffraction and determined in the monoclinic space group P12 1 /c1 with the unit-cell parameters a = 14.6650(5) Å, b=9.1916(3) Å, c = 15.0239(5) Å,β =90.2210(10)°, Z = 4 and V = 2025.13(12) Å 3 . The thermal dehydration of this compound was performed, leading to an anhydrous new form, CoNa 3 P 3 O 10 , which was characterized by X-ray diffraction. The obtained CoNa 3 P 3 O 10 crystallized in monoclinic space group P21/n with the unit-cell parameters a=15,3774 Å, b=7,6988 A°, c=14,2832A°, β=92,9115°. The characteristic IR wavenumbers of the P 3 O 10 5- ions observed in the vibrational spectra were calculated using isotopic substitutions, which confirms the existence of these groups in the studied compound. A comparison between the IR and Raman wavenumbers of CoNa 3 P 3 O 10 ·12H 2 O, CoNa 3 P 3 O 10 and Na 5 P 3 O 10 ·12H 2 was performed. A kinetic study was also made for CoNa 3 P 3 O 10 ·12H 2 O.

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

Structural Aspects of Neutron Survival Probabilities

The neutron survival probability (and related quantities including probabilities of extinction and initiation) is a central element of the broader stochastic theory of neutron populations and finds application in fields including reactor start-up, analysis of reactor power bursts and criticality accidents, and safeguards. In a full neutron transport formulation, the equation governing the single-neutron survival probability is a backward or adjoint-like integro-partial differential equation with the added complexity of being highly nonlinear. Analogous formulations of this equation exist in the context of many approximate theories of neutron transport, with the point kinetics formulation having received significant theoretical attention since the 1940s. This work continues this tradition by providing a novel analysis of the single-neutron survival probability equation using the tools of boundary layer theory. The analysis reveals that the “fully dynamic” solution of the single-neutron survival probability equation—and some key probability distributions derived from it—may be cast as a singular perturbation around the underlying quasi-static single-neutron probability of initiation. In this perturbation solution, the expansion parameter is the ratio of the neutron generation time to a macroscopic time scale characterizing the overall system evolution; this interpretation illuminates some of the fundamental structural aspects of neutron survival phenomena.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Raman and X-ray diffraction study of pressure-induced phase transition in synthetic Mg 2 TiO 4

Synthetic Mg 2 TiO 4 qandilite was investigated to 50 and 40.4 GPa at room temperature using Raman spectroscopy and X-ray diffraction, respectively. The Raman measurements showed that cubic Mg 2 TiO 4 spinel transforms to a high pressure tetragonal (I4 1 /amd, No.141) phase at 14.7 GPa. Owing to sluggish kinetics at room temperature, the spinel phase coexists with the tetragonal phase between 14.7 and 24.3 GPa. In the X-ray diffraction experiment, transformation of the cubic Mg 2 TiO 4 to the tetragonal structure was complete by 29.2 GPa, ~5 GPa higher than the transition pressure obtained by Raman measurements, owing to slow kinetics. The obtained isothermal bulk modulus of Mg 2 TiO 4 spinel is K T0 = 148(3) GPa when K T0 ’ = 6.6, or K T0 = 166(1) GPa when K T0 ’ is fixed at 4. The isothermal bulk modulus of the high-pressure tetragonal phase is calculated to be 209(2) GPa and V 0 = 270(2) Å 3 when K T0 ’ is fixed at 4, and the volume reduction on change from cubic to tetragonal phase is about 9%. The calculated thermal Grüneisen parameters (γ th ) of cubic and tetragonal Mg 2 TiO 4 phases are 1.01 and 0.63. Based on the radii ratio of spinel cations, a simple model is proposed to predict post-spinel structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Visualization and modeling of inhibition of IL-1β and TNF-α mRNA transcription at the single-cell level

Abstract IL-1β and TNF-α are canonical immune response mediators that play key regulatory roles in a wide range of inflammatory responses to both chronic and acute conditions. Here we employ an automated microscopy platform for the analysis of messenger RNA (mRNA) expression of IL-1β and TNF-α at the single-cell level. The amount of IL-1β and TNF-α mRNA expressed in a human monocytic leukemia cell line (THP-1) is visualized and counted using single-molecule fluorescent in-situ hybridization (smFISH) following exposure of the cells to lipopolysaccharide (LPS), an outer-membrane component of Gram-negative bacteria. We show that the small molecule inhibitors MG132 (a 26S proteasome inhibitor used to block NF-κB signaling) and U0126 (a MAPK Kinase inhibitor used to block CCAAT-enhancer-binding proteins C/EBP) successfully block IL-1β and TNF-α mRNA expression. Based upon this single-cell mRNA expression data, we screened 36 different mathematical models of gene expression, and found two similar models that capture the effects by which the drugs U0126 and MG132 affect the rates at which the genes transition into highly activated states. When their parameters were informed by the action of each drug independently, both models were able to predict the effects of the combined drug treatment. From our data and models, we postulate that IL-1β is activated by both NF-κB and C/EBP, while TNF-α is predominantly activated by NF-κB. Our combined single-cell experimental and modeling efforts show the interconnection between these two genes and demonstrates how the single-cell responses, including the distribution shapes, mean expression, and kinetics of gene expression, change with inhibition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

The high pressure dependence of x-ray induced decomposition of cadmium oxalate

The high proclivity of x rays to destabilize and distort molecular structures has been previously utilized in the synthesis of novel compounds. Here, we show that x-ray induced decomposition of cadmium oxalate induces chemical and structural transformations only at 0.5 and 1 GPa. Using x-ray diffraction and Raman spectroscopy, the synthesized product is identified as cadmium carbonate with cadmium oxalate remnants, which is stable under ambient conditions. At ambient and >1 GPa pressures, only degradation of the electronic density distribution is observed. The transformation kinetics are examined in terms of Avrami’s model, which demonstrates that despite the necessity of high pressure for efficient x-ray induced synthesis of cadmium carbonate, the rate and geometry of structural synthesis in the 0.5–1 GPa pressure range do not depend on the applied pressure. In addition, the possible role of intermolecular distance and molecular mobility in transformation yield is also discussed. Our experimental results indicate that x-ray induced photochemical synthetic pathways can be modulated and optimized by specific parameter selection such as high pressure.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Commensuration Effects in Layered Nanoparticle Solids.

We have developed a hierarchical nanoparticle transport simulator (HINTS) and adapted it to study commensuration effects in two classes of nanoparticle (NP) solids: (1) a bilayer NP solid (BNS) with an energy offset and (2) a BNS as part of a field effect transistor (FET). HINTS integrates the ab initio characterization of single NPs with the phonon-assisted tunneling transition model of the NP-NP transitions into a kinetic Monte Carlo based simulation of the charge transport in NP solids. First, we studied a BNS with an interlayer energy offset Delta, possibly caused by a fixed electric field. Our results include the following. First, in the independent energy-offset model, we observed the emergence of commensuration effects when scanning the electron filling factor (FF) across integer values. These commensuration effects were profound as they reduced the mobility by several orders of magnitude. We analyzed these commensuration effects in a five-dimensional parameter space, as a function of the on-site charging energy E-C, energy offset Delta, disorder D, the electron FF, and temperature k(B)T. We demonstrated the complexity of our model by showing that at integer FFs commensuration effects are present in some regions of the parameter space, while they vanish in other regions, thus defining distinct dynamical phases of the model. We determined the phase boundaries between these dynamical phases. Second, using these results as a foundation, we shifted our focus to the experimentally much-studied NP-FETs. NP-FETs are also characterized by an interlayer energy offset Delta, which, in contrast to our first model, is set by the gate voltage V-G and thereby related to the electron FF. We repeated many of our simulations and again demonstrated the emergence of commensuration effects and distinct dynamical phases in these NP-FETs. Notably, the commensuration effects in the NP-FETs showed many similarities to those in the independent energy-offset BNS.

Qu, Luman↗
DOE OSTI2022

Transport physics dependence of Bohm speed in presheath–sheath transition

We report the ion exit flow speed at the sheath entrance is constrained by the Bohm criterion, which is used as a boundary condition for simulations that do not resolve the sheath region. Traditional Bohm criterion analysis invokes the equation of state and, thus, ignores transport physics in the sheath transition problem. An expression for the Bohm speed away from the asymptotic limit is derived from a set of anisotropic plasma transport equations. The thermal force, collisional temperature isotropization, and heat flux enter into the evaluation of the Bohm speed. By comparison with kinetic simulation results, this expression is shown to be accurate in the presheath–sheath transition region rather than a single point at the sheath entrance over a broad range of collisionality.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Global Theory of Microtearing Modes in the Tokamak Pedestal

The pedestal of H-mode tokamaks displays strong magnetic fluctuations correlated with the evolution of the electron temperature. The microtearing mode (MTM), a temperature-gradient-driven instability that alters magnetic topology, is compatible with these observations. Here we extend the conventional theory of the MTM to include the global variation of the temperature and density profiles. The offset between the rational surface and the location of the pressure gradient maximum (μ) emerges as a crucial parameter for MTM stability. In conclusion, the extended theory matches observations on the Joint European Torus tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Effects of fuel composition and octane sensitivity on polycyclic aromatic hydrocarbon and soot emissions of gasoline–ethanol blend surrogates

The sooting propensity of a fuel is closely coupled with the fuel composition and chemistry. A detailed understanding of their effects is, therefore, needed to develop next-generation fuels which can minimize particulate emissions. With this overarching goal, the present work numerically investigates the effects of fuel composition and octane sensitivity (S) on polycyclic aromatic hydrocarbons (PAH) and soot emissions, for four-component gasoline-ethanol blend surrogates comprising isooctane, n-heptane, toluene, and ethanol. A partially-premixed counterflow flame is chosen as the canonical configuration for this study and simulations are performed using CHEMKIN-Pro-employing a kinetic mechanism developed by Park et al. (2017). In addition, a, detailed soot model based on the sectional method is used to capture the spatial characteristics of soot emissions. The kinetic mechanism and soot model are validated using available experimental data for various targets. A total of 86 TPRF-E mixtures, spanning a wide range of concentration of each component, and a wide range of S are analyzed. The effect of each non-paraffinic fuel component on the resultant PAH and soot emissions is investigated. PAH and soot emissions are found to vary significantly depending upon the blend composition. Additionally, based on the parametric sweeps, a regression analysis is carried out to identify global parameters that govern the formation of PAHs and soot. The analysis shows that both toluene content and S have a prominent effect on the formation of PAHs and soot, with toluene content having a stronger impact. Moreover, larger PAHs have higher dependency on toluene content and S. Furthermore, a detailed analysis is carried out to understand the physical and chemical phenomena associated with the observed trends of PAH and soot emissions.

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

Rate Coefficients for OH + NO (+N 2 ) in the Fall-off Regime and the Impact of Water Vapor

The termolecular, association reaction between OH and NO is a source of nitrous acid (HONO), an important atmospheric trace gas. Rate coefficients for the title reaction as recommended by evaluation panels differ substantially at the temperatures and pressures that prevail in the Earth’s boundary layer where the reaction is in the fall-off regime between low- and high-pressure limiting rate coefficients. Using pulsed laser methods for generation and detection of OH, we have reinvestigated the kinetics of the title reaction at pressures of 22–743 Torr (1 Torr = 1.333 hPa) and temperatures (273, 298, and 333 K) in pure N 2 and in N 2 –H 2 O bath gases. In situ optical absorption measurements were used to rule out any bias due to NO 2 or HONO impurities. Our rate coefficients (k 1 ) in N 2 bath gas are parametrized in terms of low-pressure (k o ) and high-pressure (k ∞ ) rate coefficients and a fall-off parameter (F C ) with k 1,0 N 2 = 7.24 × 10 –31 (T/300 K) –2.17 cm 6 molecule –2 s –1 , k 1,∞ = 3.3 × 10 –12 (T/300 K) –3 cm 3 molecule –1 s –1 , and F C = 0.53. Used with the “Troe” expression for termolecular reactions, these parameters accurately reproduce the current data in the fall-off regime and also capture literature rate coefficients at extrapolated temperatures. The presence of water vapor was found to enhance the rate coefficients of the title reaction significantly. The low-pressure limiting rate coefficient in H 2 O bath gas is a factor 5–6 larger than in N 2 , at room temperature (k 1,o H 2 O = 4.55 × 10 –30 (T/300 K) –4.85 cm 6 molecule –2 s –1 ) indicating that H 2 O is much more efficient in quenching the association complex HONO* through collisional energy transfer. Based on measurements in N 2 –H 2 O mixtures, a parametrization of k 1 including both N 2 and H 2 O as third-body quenchers was derived. Neglecting the effect of H 2 O results, e.g., in an underestimation of k 1 by >10% in the tropical boundary layer.

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

Schwinger-Keldysh effective action for hydrodynamics with approximate symmetries

We study hydrodynamic theories with approximate symmetries in the recently developed effective action approach on the Schwinger-Keldysh contour. We employ the method of spurious symmetry transformation for small explicit symmetry-breaking parameters to systematically constrain symmetry-breaking effects in the nonequilibrium effective action for hydrodynamics. We apply our method to the hydrodynamic theory of chiral symmetry in quantum chromodynamics at finite temperature and density and its explicit breaking by quark masses. We show that the spurious symmetry and the Kubo-Martin-Schwinger relation dictate that the Ward-Takahashi identity for the axial symmetry, i.e., the partial conservation of axial vector current (PCAC) relation, contains a relaxational term proportional to the axial chemical potential, whose kinetic coefficient is at least of the second order in the quark mass. In the phase where the chiral symmetry is spontaneously broken, and the pseudo-Nambu-Goldstone pions appear as hydrodynamic variables, this relaxation effect is subleading compared to the conventional pion mass term in the PCAC relation, which is of the first order in the quark mass. On the other hand, in the chiral symmetry restored phase, we show that our relaxation term, which is of the second order in the quark mass, becomes the leading contribution to the axial charge relaxation. Therefore, the leading axial charge relaxation mechanism is parametrically different in the quark mass across a chiral phase transition.

Goldstone bosons↗
DOE OSTI2022

Design Rules of the Mixing Phase and Impacts on Device Performance in High-Efficiency Organic Photovoltaics

In nonfullerene acceptor- (NFA-) based solar cells, the exciton splitting takes place at both domain interface and donor/acceptor mixture, which brings in the state of mixing phase into focus. The energetics and morphology are key parameters dictating the charge generation, diffusion, and recombination. It is revealed that tailoringthe electronic properties of the mixing region by doping with larger-bandgap components could reduce the density of state but elevate the filling state level, leading to improved open-circuit voltage ( V OC ) and reduced recombination. The monomolecular and bimolecular recombinations are shown to be intercorrelated, which show a Gaussian-like relationship with V OC and linear relationship with short-circuit current density ( J SC ) and fill factor (FF). The kinetics of hole transfer and exciton diffusion scale with J SC similarly, indicating the carrier generation in mixing region and crystalline domain are equally important. From the morphology perspective, the crystalline order could contribute to V OC improvement, and the fibrillar structure strongly affects the FF. These observations highlight the importance of the mixing region and its connection with crystalline domains and point out the design rules to optimize the mixing phase structure, which is an effective approach to further improve device performance.

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

Phase segregation mechanisms of small molecule-polymer blends unraveled by varying polymer chain architecture

As phase separation between the small-molecule semiconductor and the polymer binder is the key enabler of blend-based organic field-effect transistors (OFETs) fabricated by low-cost solution processing, it is crucial to understand the underlying phase separation mechanisms that determine the phase morphology, which significantly impacts device performance. Beyond the parameter space investigated in previous work, here we investigate the formation of blends by varying the branch architecture of the polymer binder and by shortening the solvent dry time using ultrasonic spray casting. The phase morphologies of the resulting blend films have been thoroughly characterized with a variety of techniques in three dimensions over multiple length scales, including AFM, energy-filtered transmission electron microscope, and neutron reflectivity, and have been correlated with electrical transport performance. From the results, we have inferred that the phase morphology is kinetically determined, limited by the inherent slow movement of polymer macromolecules. The kinetic picture, supported by molecular dynamics modeling, not only consistently explains our observations but also resolves inconsistencies in previous works. The achieved mechanistic understanding will guide further optimization of blend-based organic electronics, such as OFETs and organic photovoltaics.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Bottom-up synthesis of mesoscale nanomeshes of graphene nanoribbons on germanium

The synthesis of functional graphene nanostructures on Ge(001) provides an attractive route toward integrating graphene-based electronic devices onto complementary metal oxide semiconductor-compatible platforms. In this study, we leverage the phenomenon of the anisotropic growth of graphene nanoribbons from rationally placed graphene nanoseeds and their rotational self-alignment during chemical vapor deposition to synthesize mesoscale graphene nanomeshes over areas spanning several hundred square micrometers. Lithographically patterned nanoseeds are defined on a Ge(001) surface at pitches ranging from 50 to 100 nm, which serve as starting sites for subsequent nanoribbon growth. Rotational self-alignment of the nanoseeds followed by anisotropic growth kinetics causes the resulting nanoribbons to be oriented along each of the equivalent, orthogonal Ge < 110 > directions with equal probability. As the nanoribbons grow, they fuse, creating a continuous nanomesh. In contrast to nanomesh synthesis via top-down approaches, this technique yields nanomeshes with atomically faceted edges and covalently bonded junctions, which are important for maximizing charge transport properties. Additionally, we simulate the electrical characteristics of nanomeshes synthesized from different initial nanoseed-sizes, size-polydispersities, pitches, and device channel lengths to identify a parameter-space for acceptable on/off ratios and on-conductance in semiconductor electronics. The simulations show that decreasing seed diameter and pitch are critical to increasing nanomesh on/off ratio and on-conductance, respectively. With further refinements in lithography, nanomeshes obtained via seeded synthesis and anisotropic growth are likely to have superior electronic properties with tremendous potential in a multitude of applications, such as radio frequency communications, sensing, thin-film electronics, and plasmonics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Versatile stochastic model for predictive KMC simulation of fcc metal nanostructure evolution with realistic kinetics

Stochastic lattice-gas models provide the natural framework for analysis of the surface diffusion-mediated evolution of crystalline metal nanostructures on the appropriate time scale (often 10 1 –10 4 s) and length scale. Model behavior can be precisely assessed by kinetic Monte Carlo simulation, typically incorporating a rejection-free algorithm to efficiently handle the broad range of Arrhenius rates for hopping of surface atoms. The model should realistically prescribe these rates, or the associated barriers, for a diversity of local surface environments. However, commonly used generic choices for barriers fail, even qualitatively, to simultaneously describe diffusion for different low-index facets, for terrace vs step edge diffusion, etc. We introduce an alternative Unconventional Interaction–Conventional Interaction formalism to prescribe these barriers, which, even with few parameters, can realistically capture most aspects of behavior. Here, the model is illustrated for single-component fcc metal systems, mainly for the case of Ag. It is quite versatile and can be applied to describe both the post-deposition evolution of 2D nanostructures in homoepitaxial thin films (e.g., reshaping and coalescence of 2D islands) and the post-synthesis evolution of 3D nanocrystals (e.g., reshaping of nanocrystals synthesized with various faceted non-equilibrium shapes back to 3D equilibrium Wulff shapes).

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

Understanding cold electron impact on parallel-propagating whistler chorus waves via moment-based quasilinear theory

Earth's magnetosphere hosts a wide range of collisionless particle populations that interact through various wave-particle processes. Among these, cold electrons, with energies below 100 eV, often dominate the plasma density but remain poorly characterized due to measurement challenges such as spacecraft charging and photoelectron contamination. Understanding the contribution of these cold populations to wave–particle interaction is of significant interest. Recent kinetic simulations identified a secondary drift-driven instability, in which parallel-propagating whistler-mode chorus waves excite oblique electrostatic whistler waves near the resonance cone and Bernstein-mode turbulence. These secondary modes enable a new channel of energy transfer from the parallel-propagating whistler wave to the cold electrons. In this work, we develop a moment-based quasilinear theory of the secondary instabilities to quantify such energy exchange. Our results show that these secondary instabilities persist for a wide range of parameters and, in many cases, lead to nearly complete damping of the primary wave. Such secondary instability might limit the amplitude of parallel-propagating whistler waves in Earth's magnetosphere and might explain why high-amplitude oblique whistler or electron Bernstein waves are rarely observed simultaneously with high-amplitude field-aligned whistler waves in the inner magnetosphere.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

High-throughput multimodal exploration of a nanocrystalline Cu-Ag library

Sputter-deposited, nanocrystalline Cu-Ag thin films produced across a broad compositional and deposition-parameter space were evaluated to unravel the process-structure-property relationships important for creating hard, conductive electrical contacts and coatings. Combinatorial deposition involving pulsed direct current magnetron sputtering of elemental targets enabled swift examination of nearly the full range of alloy compositions and a relevant portion of deposition atomistics. Several high-throughput characterization modalities were employed to evaluate the chemistry, structure, and properties of the films. The resultant hardness, modulus, film density, crystal texture, and resistivity were analyzed in terms of key deposition characteristics (incident atom kinetic energy and incidence angle) predicted by binary-collision, kinematic Monte Carlo simulations. The study revealed improved hardness, parabolic resistivity dependence on composition, and compositional and process dependencies of film tarnishing. The results are discussed in the context of variations in microstructure and film density. Transmission electron microscopy and X-ray diffraction demonstrate several forms of compositional variation including solute segregation to grain boundaries as well as periodic, intragranular compositional modulations. Annealing of a Cu-rich alloy film exhibiting grain boundary segregation showed that this as-deposited, compositional variation is not stable above 100 °C. Finally, the Cu-Ag system is shown to have potential for hard, conductive, tarnish-resistant and room temperature-stable nanocrystalline thin films across the composition space.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗