Search NASA⌕ Search

SEARCH · Search NASA

Results for “kinetic simulation”

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

3D Simulation of Lower-Hybrid Drift Waves in Strong Guide Field Asymmetric Reconnection in Laboratory Experiments

The magnetic reconnection experiment has recently seen short wavelength (kp (e)(1) lower-hybrid waves near the electron diffusion region in strong guide field reconnection. Based on plasma parameters from the experiment, we perform a three-dimensional fully kinetic simulation in order to investigate the generation of the lower-hybrid waves and their effects on the reconnection process. We find that the low-beta regions around the reconnection site are unstable to the lower-hybrid drift instability propagating in the outflow direction, driven by the difference between the electron and ion outflows. The waves modify the electron distributions, leading to periodic opening and closing of gaps in electron velocity space, and provide a small contribution to the anomalous resistivity. Finally, the simulation results are discussed in the context of space observations and laboratory experiments.

J. Ng↗

Cycling phosphorus on the Archean Earth: Part II. Phosphorus limitation on primary production in Archean ecosystems

Several lines of evidence point to low rates of net primary production (NPP) in Archean oceans. However, whether Archean NPP was limited by electron donors or nutrients, particularly phosphorus (P), and how these factors might have changed over a billion years of recorded Archean history, remains contentious. One major challenge is to understand quantitatively the biogeochemical cycling of P on the early Earth. In Part I of this series (Hao et al., 2020), we estimated the weathering flux of P to the oceans as a function of temporally increasing continental emergence and elevation through Archean time. In Part II, we conduct thermodynamic and kinetic simulations to understand key processes of P cycling within the Archean ocean, including seafloor weathering, recycling of organic P, the solubility and precipitation of secondary phosphate minerals, and the burial diagenesis of P precipitates. Our calculations suggest low solubilities of apatite minerals in Archean seawater, primarily due to nearly neutral pH and high levels of Ca. This low solubility, in turn, implies a negligible contribution of apatite dissolution to P bioavailability in Archean seawater. We also simulate the solubility limits of common secondary P-bearing minerals, showing that vivianite would have been the least soluble P mineral in ferruginous Archean seawater (0.1–0.3 µM), even at moderate supersaturation states (Ω = 100 or 1000). If vivianite precipitation was kinetically favorable by microbial activities and mineral adsorption, the sinking flux of P as vivianite in Archean seawater could have reached the modern sinking flux, implying that vivianite precipitation was a potentially major sink for P in Archean oceans. During burial diagenesis, however, vivianite in porewater would have become less stable than Ca-phosphates of lower solubility. At elevated temperatures (>100 °C) associated with burial diagenesis and low-grade metamorphism, vivianite is predicted to react irreversibly with calcite to form apatite. Optimistic assumptions about the recycling efficiency of P on the Archean Earth lead us to estimate that by the end of the eon the total flux of P (continental weathering + recycling) could have supported NPP at levels up to 7% of the modern. The total flux of P would have been much lower on the early and middle Archean Earth, whereas fluxes of electron donors could have been higher, suggesting very low productivity and P-limitation of marine ecosystems during much of the eon. Comparing our estimates of NPP as limited by P supply with the estimate by Ward et al. (2019), in which NPP was limited by electron donors and metabolic efficiency, there could have been a transition between P-limited productivity in the early to middle Archean to electron donor-limitation closer to the eon’s end (assuming no oxygenic photosynthesis). Once oxygenic photosynthesis reached ecological significance, probably near the end of the Archean, our estimated flux of P would allow rapid oxidation of atmosphere.

Vivianite↗

Kinetic Monte Carlo Simulations of Rod Eutectics and the Surface Roughening Transition in Binary Alloys

In this paper we report on two related topics: Kinetic Monte Carlo simulations of the steady state growth of rod eutectics from the melt, and a study of the surface roughness of binary alloys. We have implemented a three dimensional kinetic Monte Carlo (kMC) simulation with diffusion by pair exchange only in the liquid phase. Entropies of fusion are first chosen to fit the surface roughness of the pure materials, and the bond energies are derived from the equilibrium phase diagram, by treating the solid and liquid as regular and ideal solutions respectively. A simple cubic lattice oriented in the {100} direction is used. Growth of the rods is initiated from columns of pure B material embedded in an A matrix, arranged in a close packed array with semi-periodic boundary conditions. The simulation cells typically have dimensions of 50 by 87 by 200 unit cells. Steady state growth is compliant with the Jackson-Hunt model. In the kMC simulations, using the spin-one Ising model, growth of each phase is faceted or nonfaceted phases depending on the entropy of fusion. There have been many studies of the surface roughening transition in single component systems, but none for binary alloy systems. The location of the surface roughening transition for the phases of a eutectic alloy determines whether the eutectic morphology will be regular or irregular. We have conducted a study of surface roughness on the spin-one Ising Model with diffusion using kMC. The surface roughness was found to scale with the melting temperature of the alloy as given by the liquidus line on the equilibrium phase diagram. The density of missing lateral bonds at the surface was used as a measure of surface roughness.

Bentz, Daniel N.↗

A moment-conserving discontinuous Galerkin representation of the relativistic Maxwellian distribution

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell–Jüttner distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here, we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

astrophysical plasmas↗

Kinetic corrections to heat-flow and Nernst advection for laser heated plasmas

Reduced models for approximating the impact of kinetic electron behavior on the transport of thermal energy and magnetic field are investigated. The thermal flux limiter has improved agreement with Vlasov–Fokker–Planck data when a harmonic form is used that adjusts the electron mean free path to account for electron–electron collisions; these results apply to both unmagnetized and magnetized plasmas. Once a magnetic field is incorporated, the mean free path should also be modified using the electron gyroradius. A flux limiter on Nernst advection of magnetic fields is also required; a form that limits Nernst by the same fraction as the thermal heat-flow best reproduces kinetic simulations. A flux limiter form for the cross terms (Righi–Leduc and cross-gradient-Nernst) is also suggested. Hohlraum simulations relevant to fusion experiments on the National Ignition Facility are found to be sensitive to all of these details.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The promising chemical kinetics for the simulation of propane-air combustion with KIVA-II code

The development of chemical kinetics for the simulation of propane-air combustion with the use of computer code KIVA-II since 1989 is summarized here. In order to let readers understand the general feature well, a brief description of the KIVA-II code, specially related with the chemical reactions is also given. Then the results of recent work with 20 reaction mechanism is presented. It is also compared with the 5 reaction mechanism. It may be expected that the numerical stability of the 20 reaction mechanism is better as compared to that of 5 reaction mechanism, but the CPU time of the CRAY computer is much longer. Details are presented in the paper.

Ying, S. J.↗

Application of a First-Principles Anomalous Transport Model for Electrons to Multiple Hall Thrusters and Operating Conditions

We have developed a physics-based model based on a pseudo-particle description of the electron cyclotron drift instability. A key improvement of the model with respect to previous work is that linear theory is not applied in the event of wave saturation and deviations of electrons or ions from a Maxwellian distribution function. In the acceleration region, the anomalous collision frequency is computed as the minimum value necessary to prevent the electron drift velocity from exceeding the thermal velocity. A functional based on the electron equilibration time is defined to control the transition from high to low resistivity regions. The model was previously applied to a single Hall thruster at its nominal operating condition, showing promising results that captured accurately the location of the thruster’s acceleration region. In this paper, we extend the use of this firstprinciples models to two additional thrusters, also considering multiple operating conditions for each of them. Numerical results are compared to experimental measurements obtained with non-invasive laser induced fluorescence. In general, the agreement between experiments and simulations is good. The model is able to predict the location of the acceleration region for all cases. We observe however that fine details, such as changes in the plasma potential gradient within the acceleration regions, are not captured. The model is also insensitive to changes in the magnetic field strength while experiments show that small shifts in location (of less than 5% of the acceleration channel length) occur. We plan to address the weaknesses of our method with the help of physical insight gained from kinetic simulations of the acceleration region.

Chaplin, Vernon H.↗

Particle Acceleration by Dense Impulsive Structures Moving in Ambient Magnetospheric Plasma. 3-D Hybrid Kinetic Modeling and MMS Observations

High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.

A. S. Lipatov↗

Off-Equilibrium Reactivity of Boron-Enriched Metal Diboride Surfaces in Electroreduction Conditions

Boron-based materials, featuring B-dependent reactivity and diverse phases, are emerging as promising catalyst systems. However, the catalytic mechanism on many borides remains poorly understood due to complex surface reconstructions under reaction conditions. Here, we investigate the MoB 2 surface in conditions of hydrogen evolution reaction in acidic media, using grand canonical global optimization, grand canonical density functional theory, ab initio molecular dynamics, free energy surface sampling, and an analytical model for electrochemical barrier evaluation. We propose a boron-enrichment strategy to tune the surface reactivity of the hexagonal face of MoB 2 . We reveal the dynamic nature of the B-enriched surface under H coverage and kinetic trapping of the system in the metastable regime with an extensive examination of the deactivation pathways. The metastable center B site on B-enriched surfaces, featuring buckled-up configuration and a usual relaxation effect, is found to be highly active toward HER via the Volmer–Heyrovsky mechanism. In conclusion, this work demonstrates how off-equilibrium behaviors can arise from the interplay between adsorbate coverage and surface reconstruction on a seemingly simple surface, and we present a theoretical framework and computational workflows to address these behaviors, along with other realistic complexities, in kinetics simulations.

Adsorption↗

Concerning the Charging of an Exploration Craft on and near a Small Asteroid

Introduction: An object immersed in an airless plasma environment will experience a natural process of surface charging in order to acheieve current balance, or zero net electric current to the object. It has been shown in recent computer simulations that the small-body plasma environment is very complex [1], considering effects of photoemission, topography, and formation of a plasma wake. For this work we consider an exploration craft (or astronaut) immersed within a plasma environment near an asteroid, which exhibits widely varying solar wind and photoelectric particle fluxes and continuously evolving illumination conditions. Objective: We aim to determine how an explo-ration craft or astronaut suit accumulates charge while located in the "nightside" asteroid wake where the particle fluxes are reduced, and in the dayside near-surface photoelectron sheath, by combining an object charging model [2] with kinetic simulations of a near-asteroid plasma environment [1]. We consider an astronaut floating near the asteroid while not in contact with the surface, as well as an astronaut moving along the surface using their hands/gloves to crawl along. Results: The modeling results suggest that remediation of triboelectric charge via accumulation of plasma currents is an important factor to consider when designing future NEA mission infrastructure, especially if repeated and frequent contact with the surface is planned. In shadowed regions such as the location shown in Fig. 1a, the plasma currents are so low (and the effective charge-remediation timescale so long, e.g. minutes to hours) that repeated contact with the surface tribocharges the glove in an uncontrollable fashion, as shown for two representative electron temperatures in Fig. 2a. The resulting buildup of significant negative charge would eventually initiate some other "current of last resort" [4] such as transport of positively-charged dust, field-emission from the glove, or significant alteration of environmental ion currents within the wake. In contrast, the few-meters-thick dayside photoelectron sheath in which the astronaut of Fig. 1b is immersed in is so rich in electrons (and hence so electrically conductive) that accumulated tribocharge dissipates almost instantaneously (e.g. in less than a ms) as shown in Fig. 2b. As our model astronaut orbits the NEA they would experience plasma currents and associated charge re-mediation times spanning many orders of magnitude, and the fusion between our numerical models provides a detailed understanding of the charging hazards possibly associated with contact-based NEA exploration.

Asteroid↗

Ohm’s Law, the Reconnection Rate, and Energy Conversion in Collisionless Magnetic Reconnection

Magnetic reconnection is a ubiquitous plasma process that transforms magnetic energy into particle energy during eruptive events throughout the universe. Reconnection not only converts energy during solar flares and geomagnetic substorms that drive space weather near Earth, but it may also play critical roles in the high energy emissions from the magnetospheres of neutron stars and black holes. In this review article, we focus on collisionless plasmas that are most relevant to reconnection in many space and astrophysical plasmas. Guided by first-principles kinetic simulations and spaceborne in-situ observations, we highlight the most recent progress in understanding this fundamental plasma process. We start by discussing the non-ideal electric field in the generalized Ohm’s law that breaks the frozen-in flux condition in ideal magnetohydrodynamics and allows magnetic reconnection to occur. We point out that this same reconnection electric field also plays an important role in sustaining the current and pressure in the current sheet and then discuss the determination of its magnitude (i.e., the reconnection rate), based on force balance and energy conservation. This approach to determining the reconnection rate is applied to kinetic current sheets with a wide variety of magnetic geometries, parameters, and background conditions. We also briefly review the key diagnostics and modeling of energy conversion around the reconnection diffusion region, seeking insights from recently developed theories. Finally, future prospects and open questions are discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Artificial correlation heating in PIC simulations

The Particle-in-Cell (PIC) method, a cornerstone in plasma modeling, is widely employed for its ability to simulate kinetic phenomena in device-scale domains. Part of what makes this possible is that computational macroparticles represent many physical particles. It converges under certain constraints, including a grid spacing that resolves the Debye length and a time step small enough to respect the Courant–Friedrichs–Lewy condition and plasma frequency stability limit. Here, we introduce a new constraint necessary to avoid Artificial Correlation Heating (ACH). This requires that the macroparticle coupling strength be smaller than one, Γ ω < 1, where Γ ω ≡ Γω 2/3 ⁠, Γ = Z 2 e 2 /(4πε 0 ak B T) is the physical coupling strength, and w is the macroparticle weight. This is particularly relevant to 3D simulations of dense plasmas, which are becoming common with modern computing power. If this condition is violated, the finite macroparticle weight artificially enhances the coupling strength and causes the plasma to heat until the macroparticle coupling strength is near unity, depending on the grid resolution. A comprehensive model of ACH is developed that incorporates electron density, temperature, macroparticle weight, and grid resolution. It is then tested using PIC simulations, delineating the boundaries of the method's applicability and offering a predictive framework for ACH. Moreover, the research explores a runaway heating process induced by ACH in the presence of ionization, which can lead to numerical instability. A conclusion of this study is that the onset of ACH can impose a more stringent constraint on the macroparticle weight and average number of macroparticles per cell than what is typically expected, particularly in 3D simulations of dense plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy Partition in Collisionless Counterstreaming Plasmas

Fast, counterstreaming plasma outflows drive magnetic field amplification, plasma heating, and particle acceleration in numerous astrophysical environments, from supernova remnant shocks to active galactic nuclei jets. Understanding how, in the absence of Coulomb collisions, energy is redistributed between the different plasma species remains a fundamental open question. We use three-dimensional, fully kinetic simulations to investigate energy partition in weakly magnetized counterpropagating plasmas. Our results reveal a complex interplay between different processes, where at early times the Weibel instability drives a first stage of magnetic field amplification and at late times the kinking of current filaments drives a second amplification stage via a dynamo-type mechanism. Electrons are heated primarily during the latter phase through magnetic pumping. By the time the flows thermalize, we observe that the final temperature ratio T e /T i and energy partition depend on the ion-to-electron mass ratio. For electron–proton flows, the electron thermal energy only reaches up to a few percent of the initial ion kinetic energy.

Marret, Alexis [SLAC National Accelerator Laborato↗

Electron-scale Magnetic Holes Generation Driven by Whistler-to-Bernstein Mode Conversion in Fully Kinetic Plasma Turbulence

Magnetic holes (MHs) are coherent structures characterized by a strong and localized magnetic field amplitude dip, commonly observed in the heliosphere. These structures come in different sizes, from magnetohydrodynamic to kinetic scales. Subion-scale MHs are usually sustained by an electron current vortex and exhibit a strong electron temperature anisotropy, with higher temperatures perpendicular to the background magnetic field. Magnetospheric multiscale observations (MMSs) have revealed electron-scale MHs to be ubiquitous in the turbulent Earth’s magnetosheath and the solar wind, potentially playing an important role in the energy cascade and dissipation. Despite abundant observations, the origin of electron-scale MHs is still unclear and debated. In this work, we use fully kinetic simulations to investigate the role of plasma turbulence in generating electron-scale MHs. We find that the turbulence spontaneously produces electron-scale MHs via the following mechanism: first, large-scale turbulent velocity shears produce regions with high electron temperature anisotropy; these localized regions become unstable, generating oblique electron-scale whistler waves; as they propagate over the inhomogeneous turbulent background, whistler fluctuations develop an electrostatic component, turning into Bernstein-like modes; the strong electrostatic fluctuations produce current filaments that merge into an electron-scale current vortex; the resulting electron vortex locally reduces the magnetic field amplitude, finally evolving into an electron-scale MH. We show that MHs generated by this mechanism have properties consistent with MMSs and nontrivial kinetic features with a “mushroom”-shaped electron velocity distribution function. Our results have potential implications for understanding the formation and occurrence of electron-scale MHs in astrophysical turbulent and space environments, such as the Earth’s magnetosheath and the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

A multiscale model to understand the interface chemistry, contacts, and dynamics during lithium stripping

A reversible Li-metal electrode, paired with a solid electrolyte, is critical for attaining higher energy density and safer batteries beyond the current lithium-ion cells. A stable stripping process may be even harder to attain as the stripping process will remove Li-atoms from the surface, and naturally reduce surface contact area, if not self-corrected by other mechanisms, such as diffusion and plastic deformation under an applied external stack pressure. Here, we capture these mechanisms occurring at multiple length- and time- scales, i.e., interface interactions, vacancy hopping, and plastic deformation, by integrating density functional theory (DFT) simulations, kinetic Monte Carlo (KMC), and continuum finite element method (FEM). By assuming the self-affine nature of multiscale contacts, we predict the steady-state contact area as a function of stripping current density, interface wettability, and stack pressure. We further estimate the exponential increase of overpotential due to contact area loss to maintain the same stripping current density. We demonstrate that a lithiophilic interface requires less stack pressure to reach the same steady-state contact area fraction than a lithiophobic interface. A “tolerable steady-state” contact area loss for maintaining stable stripping is estimated at 20 %, corresponding to a 10 % increase in overpotential. To constrain contact loss within the tolerance, the required stack pressure is 0.1, 0.5, and 2 times the yield strength of lithium metal for three distinct interfaces, lithiophilic Li/lithium oxide(Li2O), Li/lithium lanthanum zirconium oxide(LLZO), and lithiophoblic Li/lithium fluoride(LiF), respectively. The modeling results agree with experiments on the impact of the stack pressure quantitatively, while the discrepancy in stripping rate sensitivity is attributed to the simplifying interface interaction in our simulations. Overall, this multiscale simulation framework demonstrates the importance of electrochemical-mechanical coupling in understanding the dynamics of the Li/SE interface during stripping.

Feng, Min↗

Understanding the interplay between pilot fuel mixing and auto-ignition chemistry in hydrogen-enriched environment

The diesel-piloted dual-fuel compression ignition combustion strategy is well-suited to accelerate the decarbonization of transportation by adopting hydrogen as a renewable energy carrier into the existing internal combustion engine with minimal engine modifications. Despite the simplicity of engine modification, many questions remain unanswered regarding the optimal pilot injection strategy for reliable ignition with minimum pilot fuel consumption. The present study uses a single-cylinder heavy-duty optical engine to explore the phenomenology and underlying mechanisms governing the pilot fuel ignition and the subsequent combustion of a premixed hydrogen-air charge. The engine is operated in a dual-fuel mode with hydrogen premixed into the engine intake charge with a direct pilot injection of n-heptane as a diesel pilot fuel surrogate. Optical diagnostics used to visualize in-cylinder combustion phenomena include high-speed IR imaging of the pilot fuel spray evolution as well as high-speed HCHO* and OH* chemiluminescence as indicators of low-temperature and high-temperature heat release, respectively. Three pilot injection strategies are compared to explore the effects of pilot fuel mass, injection pressure, and injection duration on the probability and repeatability of successful ignition. The thermodynamic and imaging data analysis supported by zero-dimensional chemical kinetics simulations revealed a complex interplay between the physical and chemical processes governing the pilot fuel ignition process in a hydrogen containing charge. Hydrogen strongly inhibits the ignition of pilot fuel mixtures and therefore requires longer injection duration to create zones with sufficiently high pilot fuel concentration for successful ignition. Results show that ignition typically tends to rely on stochastic pockets with high pilot fuel concentration, which results in poor repeatability of combustion and frequent misfiring. In conclusion, this work has improved the understanding on how the unique chemical properties of hydrogen pose a challenge for maximization of hydrogen’s energy share in hydrogen dual-fuel engines and highlights a potential mitigation pathway.

33 ADVANCED PROPULSION SYSTEMS↗

Impact of gasoline composition on the effects of nitric oxide on autoignition and knock in a DISI engine

Modern spark-ignition engines use exhaust gas recirculation (EGR) to dilute the charge and suppress knock, enabling the use of higher compression ratios and/or more optimum combustion phasing for higher efficiency. The effectiveness of EGR is affected by the composition of the fuel and its chemical-kinetic interactions with combustion products. Among those, nitric oxide (NO) has been shown to strongly affect autoignition reactivity. However, the impact of fuel composition of the effect of NO on reactivity is not well-understood. Here, in this study, engine experiments were conducted to assess the impact of NO seeded to the intake on knock-limited operation of two gasoline fuels (high cycloalkane content, or HCA, and high olefin content, or HO). Results showed that compositionally-different fuels responded differently to NO. HCA, which was less knock-limited than HO for NO < 200 ppm, became more knock-limited for NO > 200 ppm. Moreover, it was found that differences in knock between fuels were caused by differences in autoignition chemistry and not in the sequential autoignition process of the end gas that occurs due to thermal stratification. Chemical kinetic simulations were performed to better understand the experimental results. For HCA, intermediate-temperature heat release had a greater impact on autoignition reactivity than low-temperature heat release, while the opposite was observed for HO. For both fuels, NO enhances the magnitude of low-temperature heat release via NO + HO 2 → NO 2 + OH. The effect of NO on reactivity was stronger for HCA because OH produced from NO helped to overcome the OH quenching effect of cyclopentane, a main species in HCA. In contrast, HO had relatively strong inherent low-temperature chemistry arising from iso-octane, which reduced the impact of NO on reactivity. For the range of NO mole fractions tested in this study, in-cylinder NO increased fuel’s knock propensity, especially for fuels with mild low-temperature chemistry.

Autoignition↗

Investigating the potential of a higher reactivity fuel to achieve faster heat-up of aftertreatment systems

Shortening the aftertreatment system heat-up period will be crucial to meet upcoming emissions regulations, which mandate significant reductions in tail-pipe emissions for on-road heavy-duty vehicles. This work studies the impact of increasing fuel reactivity on the combustion of post injections during the expansion stroke. The work demonstrates an operating strategy using a blend of di-butyl ether (a potential low-carbon biofuel component) and #2 diesel that achieves higher exhaust enthalpy, while maintaining similar engine-out emissions (NOx and combustion efficiency) in comparison to the baseline #2 diesel fuel (cetane number 41). On average, a 12% increase in exhaust enthalpy was observed for the blend (cetane number 60) relative to the baseline fuel. The ability to operate the higher reactivity fuel at a more retarded post-injection timing, where combustion with #2 diesel was found to be more unstable, and less complete was key to achieving higher exhaust enthalpy. However, a small fuel penalty (3%–5%) was also associated with the higher reactivity fuel with this operating strategy. The formation of fuel lean regions, due to the long ignition delays at late injection conditions, was predicted to be the leading contributor to the increased emissions of products of partial combustion. Here, chemical kinetic simulations using surrogates for the fuels used in the experiments, demonstrated the ability of a higher cetane fuel to outperform the baseline fuel under fuel lean conditions with equivalence ratios less than 0.7.

09 BIOMASS FUELS↗