Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactive capture”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Simulations of Low Mach Number Reactive Flows Coupled with Electric Fields: Preprint

Applying electric fields to flames can provide an effective, non-invasive control of the combustion process. In this work, an algorithm recently introduced to handle the interactions between reactive flows and an electric field is extended to handle multiple dimensions and adaptive mesh refinement (AMR), and applied to a series of laminar flames. Experimental data allows to validate the ability of the algorithm to capture the mechanisms driving the flame response to electric fields.

electric field↗
DOE OSTI2026

Development of Machine-Learned Interatomic Potentials to Predict Structure, Transport, and Reactivity in Platinum-Based Fuel Cells

Machine-learned interatomic potentials (MLIPs) have rapidly progressed in accuracy, speed, and data efficiency in recent years. However, training robust MLIPs in multicomponent systems remains a challenge. In this work, we train an MLIP to describe hydrated Nafion ionomers and platinum catalysts, which are important components of fuel cells, by constructing a diverse training set to describe the bulk polymer and interfacial catalyst–polymer interactions well. We use our trained MLIP to study the properties of the platinum–Nafion system, including polymer structure, proton mobility in a bulk Nafion polymer and near a platinum-Nafion interface, and reactions near and far from the interface, finding excellent results for structure and reactions contained within our training set. Transport seems to be well described, with both vehicular transport and Grotthuss hopping captured, although converged calculations of diffusivities were not computed because they require calculations of tens of nanoseconds that are challenging with current state-of-the-art MLIPs. The combined insights that this model provides can be leveraged to optimize fuel cell performance, and the approach can be applied to other chemical processes and devices where structure, transport, and reactivity all contribute to the overall observed performance.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2024

Reactive flash sintering and characterization of bulk high entropy nitrides

Over the past decade, numerous high-entropy ceramics have been synthesized, often displaying attractive properties. However, the study on facile preparation of bulk high entropy nitrides (HEN) are limited, despite its broad potential applications. This research demonstrates for the first time rapid fabrication (within ∼6 min) of bulk high-entropy nitrides, especially (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N, from binary nitride powder mixtures using a highly efficient reactive flash sintering (RFS) technique. X-ray diffraction (XRD) shows the HENs from RFS are near single-phase solid solutions with a rock salt crystal structure, while in situ synchrotron study carried out during RFS captured in real time the formation of HEN, which was preserved upon cooling, suggesting thermodynamic stability of the HEN phase, even up to extreme pressure (∼35.6 GPa). Microscopic analyses using SEM, STEM, and EDS reveal decent uniformity for HEN with no obvious segregation of elements, even to submicron scale. Some properties of the obtained bulk HENs are consistent with expectations. For example, their hardness and bulk modulus are close to estimates based on rule-of-mixture (ROM) values from the constituent binary nitrides. Meanwhile, some other measured properties seem to show surprises. For example, the fracture toughness for the HENs (e.g., 7.81 ± 1.40 MPa•m 1/2 or higher) turns out to be more than double of the expected ROM estimates. The significantly improved fracture toughness is attributed to the observed nano-layered structure of the HENs, despite the HEN’s cubic crystal structure and high hardness. In addition, the oxidation resistance shows improvement up till ∼800°C, possibly due to Ta doping that suppress oxygen vacancy formation in the oxide shell, while the 5-metal HEN of (Al 0.17 Nb 0.17 Ta 0.17 Ti 0.32 Zr 0.17 )N displays superconductivity (T c of ∼5–7 K from magnetism and resistivity measurements, slightly lower than ROM estimate), despite insulating property of starting AlN. Furthermore, future study combining experimental investigation using larger samples to confirm the observed increase in fracture toughness and oxidation resistance, theoretical modeling at different length scale, and more detailed structural/chemical characterization, especially at the atomic scale, are all needed to fully understand the inter-relationships between composition, processing, structure, and novel properties for these HENs and the development of related new materials for different applications.

Flash sintering↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Detonation performance experiments and modeling for the DAAF-based high explosive PBX 9701

Detonation performance experiments and modeling are reported for the explosive PBX 9701, which is composed of 97% 3,3$\acute{-}$diamino-4,4$\acute{-}$azoxyfurazan (DAAF) and 3% vinylidene uoride-chlorotriuoroethylene copolymer (Kel-F 800) binder by weight. PBX 9701 is a newly-developed reducedsensitivity explosive with increased performance relative to the triaminotrinitrobenzene or TATB-based PBX 9502 while still retaining low sensitivity to mechanical insult. The first detonation performance measurements for this formulation are presented, including front-curvature ratesticks and cylinder expansion test data. Prior shock initiation data is also reviewed. These data are used to develop programmed burn and reactive flow calibrations for existing commonly-used performance models which allow engineering calculations with PBX 9701. The calibration process involves several enhancements relative to conventional approaches including the use of an analytical scaling correlation to speed the calibration process and development of a new methodology to improve the consistency between the programmed burn and reactive model calibration process and the resulting calculations. Overall, PBX 9701 is seen to exhibit improved performance relative to insensitive explosives, with a trinitrotoluene (TNT) equivalence of 1.24. The detonation performance properties are found to be well captured by existing models.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

A theoretical kinetic study of ĊH 3 + ṄH 2 : From electronic structure to NH 3 /CH 4 combustion modelling implications

Carbon–nitrogen interaction reactions play an important role in governing the reactivity of ammonia blended fuels. However, there remains uncertainties regarding their detailed reaction pathways and rate constants, hampering the development of high-fidelity chemical kinetic models. In this study, the kinetics of ĊH 3 + ṄH 2 , a key C–N interaction reaction in ammonia/methane blend combustion have been investigated. The potential energy surface has been explored using the high-level ANL0F method, yielding highly accurate stationary point energies that agree with ATcT values within 0.1 kcal mol –1 . Variable reaction coordinate transition state theory is used to treat the barrierless association and decomposition reaction channels, based on directly sampled radical-radical interaction energies at the CASPT2-F12(2e,2o)/cc-pVTZ-F12 level of theory. The minimum transitional mode numbers of states obtained are then coupled with the RRKM/master equation to calculate temperature- and pressure-dependent rate constants. Our a priori calculations capture available experimental measurements from the literature very well. The calculated rate constants have been incorporated into an NH 3 /CH 4 chemical kinetic model currently under development at the University of Galway. The effect of the updated kinetic data for ĊH 3 + ṄH 2 on model predicted NH 3 /CH 4 fuel reactivity is elucidated.

ab initio↗
DOE OSTI2021

Screening candidate supplementary cementitious materials under standard and accelerated curing through time-series surface resistivity measurements and change-point detection

Time-series surface resistivity (SR) measurements are evaluated as an indicator of pozzolanicity in concretes containing 13 materials of both known and unknown pozzolanicity, under standard and accelerated curing conditions, and the results were validated against portlandite consumption obtained through thermogravimetric analysis. Change-point detection (CPD) analysis is applied to the data to provide a rigorous indicator of hydration age of the onset and rate of the pozzolanic reaction, if occurring. The SR measurements were able to capture the onset of the pozzolanic reactions in concretes curing at elevated temperature regimes, which allows determination of pozzolanicity weeks earlier than at standard curing temperatures. Using CPD along with slope analysis allows a comprehensive analysis of a material's reactivity and determination of pozzolanicity in as quickly as two days. Time-series SR measurements at elevated curing conditions, along with limits on material composition, are recommended for determination of pozzolanic reactivity in novel or variable materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Multi-reward reinforcement learning based development of inter-atomic potential models for silica

Abstract Silica is an abundant and technologically attractive material. Due to the structural complexities of silica polymorphs coupled with subtle differences in Si–O bonding characteristics, the development of accurate models to predict the structure, energetics and properties of silica polymorphs remain challenging. Current models for silica range from computationally efficient Buckingham formalisms (BKS, CHIK, Soules) to reactive (ReaxFF) and more recent machine-learned potentials that are flexible but computationally costly. Here, we introduce an improved formalism and parameterization of BKS model via a multireward reinforcement learning (RL) using an experimental training dataset. Our model concurrently captures the structure, energetics, density, equation of state, and elastic constants of quartz (equilibrium) as well as 20 other metastable silica polymorphs. We also assess its ability in capturing amorphous properties and highlight the limitations of the BKS-type functional forms in simultaneously capturing crystal and amorphous properties. We demonstrate ways to improve model flexibility and introduce a flexible formalism, machine-learned ML-BKS, that outperforms existing empirical models and is on-par with the recently developed 50 to 100 times more expensive Gaussian approximation potential (GAP) in capturing the experimental structure and properties of silica polymorphs and amorphous silica.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Prediction of the Thrust Performance and the Flowfield of Liquid Rocket Engines

In an effort to improve the current solutions in the design and analysis of liquid propulsive engines, a computational fluid dynamics (CFD) model capable of calculating the reacting flows from the combustion chamber, through the nozzle to the external plume, was developed. The Space Shuttle Main Engine (SSME) fired at sea level, was investigated as a sample case. The CFD model, FDNS, is a pressure based, non-staggered grid, viscous/inviscid, ideal gas/real gas, reactive code. An adaptive upwinding differencing scheme is employed for the spatial discretization. The upwind scheme is based on fourth order central differencing with fourth order damping for smooth regions, and second order central differencing with second order damping for shock capturing. It is equipped with a CHMQGM equilibrium chemistry algorithm and a PARASOL finite rate chemistry algorithm using the point implicit method. The computed flow results and performance compared well with those of other standard codes and engine hot fire test data. In addition, the transient nozzle flowfield calculation was also performed to demonstrate the ability of FDNS in capturing the flow separation during the startup process.

Wang, T.-S.↗
DOE OSTI2024

Diverse organic carbon dynamics captured by radiocarbon analysis of distinct compound classes in a grassland soil

Soil organic carbon (SOC) is a large, dynamic reservoir composed of a complex mixture of plant- and microbe-derived compounds with a wide distribution of cycling timescales and mechanisms. The distinct residence times of individual carbon components within this reservoir depend on a combination of factors, including compound reactivity, mineral association, and climate conditions. To better constrain SOC dynamics, bulk radiocarbon measurements are commonly used to trace biosphere inputs into soils and to estimate timescales of SOC cycling. However, understanding the mechanisms driving the persistence of organic compounds in bulk soil requires analyses of SOC pools that can be linked to plant sources and microbial transformation processes. Here, we adapt approaches, previously developed for marine sediments, to isolate organic compound classes from soils for radiocarbon ( 14 C) analysis. We apply these methods to a soil profile from an annual grassland in Hopland, California (USA), to assess changes in SOC persistence with depth (down to 1 m). We measured the radiocarbon values of water-extractable organic carbon (WEOC), total lipid extracts (TLEs), total hydrolyzable amino acids (AAs), and an acid-insoluble (AI) fraction from bulk and physically separated size fractions (< 2 mm, 2 mm–63 µm, and < 63 µm). Our results show that Δ 14 C values of bulk soil, size fractions, and extracted compound classes became more depleted with depth, and individual SOC components have distinct age–depth distributions that suggest distinguishable cycling rates. We found that AAs and TLEs cycle faster than the bulk soils and the AI fraction. The AI was the most 14 C-depleted fraction, indicating that it is the most chemically inert in this soil. Our approach enables the isolation and measurement of SOC fractions that separate functionally distinct SOC pools that can cycle relatively quickly (e.g., plant and microbial residues) from more passive or inert SOC pools (associated with minerals or petrogenic) from bulk soils and soil physical fractions. With the effort to move beyond SOC bulk analysis, we find that compound class 14 C analysis can improve our understanding of SOC cycling and disentangle the physical and chemical factors driving OC cycling rates and persistence.

58 GEOSCIENCES↗
DOE OSTI2025

Enabling BWR fuel rod analysis in the BISON fuel performance code

Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

MgO Nanostructures on Cu(111): Understanding Size- and Morphology-Dependent CO 2 Binding and Hydrogenation

To design and optimize cost-effective technologies for the capture, utilization, and storage of carbon dioxide (CO 2 ), we need a fundamental knowledge and control of chemical interactions associated with the capture and conversion of the molecule into high-value chemicals, minerals, and all kinds of materials. Bulk magnesium oxide (MgO) is frequently used for the trapping and storage of CO 2 by generation of magnesium carbonates. In this study, the growth and reactivity of MgO nanostructures on a Cu 2 O/Cu(111) substrate were investigated using scanning tunneling microscopy (STM) and synchrotron-based ambient-pressure X-ray photoelectron spectroscopy (AP-XPS). For extremely small concentrations of Mg (~ 0.01 monolayer (ML)), a well-ordered film of copper oxide with small clusters (0.2-0.5 nm in width, 0.4-0.6 Å in height) of embedded MgO was seen. At a coverage of 0.1 ML, MgO nanoparticles with a width of 0.4 to 1 nm and a height of ~ 1.5 Å were randomly distributed on the copper oxide. Further, random distribution was also observed when the MgO coverage was raised to 0.25 ML, with the width of the MgO particles increasing to 2-2.5 nm and the height reaching 2 Å. These oxide nanostructures displayed a high reactivity towards CO 2 and H 2 that is not seen for bulk MgO. Dissociation of H 2 was observed at room temperature with reaction of the H adatoms with CuO x and C-containing groups. On the small MgO nanostructures (< 1 nm in width), instead of plain carbonate formation, there was dissociation of CO 2 into CO and C species, opening reaction channels for the conversion of this harmful molecule into oxygenates and light alkanes.

03 NATURAL GAS↗
DOE OSTI2023

Microwave Regeneration and Thermal and Oxidative Stability of Imidazolium Cyanopyrrolide Ionic Liquid for Direct Air Capture of Carbon Dioxide

Abstract Understanding the oxidative and thermal degradation of CO 2 sorbents is essential for assessing long‐term sorbent stability in direct air capture (DAC). The potential degradation pathway of imidazolium cyanopyrrolide, an ionic liquid (IL) functionalized for superior CO 2 capacity and selectivity, is evaluated under accelerated degradation conditions to elucidate the secondary reactions that can occur during repetitive absorption‐desorption thermal‐swing cycles. The combined analysis from various spectroscopic, chromatographic, and thermal gravimetric measurements indicated that radical and S N 2 mechanisms in degradation are encouraged by the nucleophilicity of the anion. Thickening of the liquid and gas evolution are accompanied by 50 % reduction in CO 2 capacity after a 7‐day exposure to O 2 under 80 °C. To prevent long exposure to conventional thermal heating, microwave (MW) regeneration of the CO 2 ‐reactive IL is used, where dielectric heating at 80 and 100 °C rapidly desorbs CO 2 and regenerates the IL without any measurable degradation.

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

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Final Design for Thermal/Epithermal eXperiments using High 240 Pu Content Plutonium/Aluminum Zero Power Research Reactor Plates with Polyethylene Moderator (IER 520 Final Design CED-2 Report)

The US Department of Energy Nuclear Criticality Safety Program (NCSP) convened a multinational Thermal Epithermal eXperiments (TEX) meeting in July of 2011 to discuss the data and experimental needs of criticality safety practitioners. The number one and two priority integral experiment data needs were for 239 Pu and 240 Pu, with special emphasis on cross section performance in the intermediate energy range (from 0.625 eV to 100 keV). LLNL measured five critical configurations with LANL for the plutonium test bed (IER-184) and published the experiments as International Criticality Safety Benchmark Evaluation Project evaluation PUMET-MIXED-002. Modeling of the benchmark configurations using ENDF/B-VIII.0 nuclear data showed significant overprediction of reactivity for configurations that had a large percentage of fissions in the intermediate energy regime. This report documents a variation on the TEX plutonium test bed to provide a test of 240 Pu cross sections, with sensitivity of the configuration to 240 Pu radiative capture and fission cross sections a priority for the design.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2002

Quantum Dynamics Study of the Isotopic Effect on Capture Reactions: HD, D2 + CH3

Time-dependent wave-packet-propagation calculations are reported for the isotopic reactions, HD + CH3 and D2 + CH3, in six degrees of freedom and for zero total angular momentum. Initial state selected reaction probabilities for different initial rotational-vibrational states are presented in this study. This study shows that excitations of the HD(D2) enhances the reactivities; whereas the excitations of the CH3 umbrella mode have the opposite effects. This is consistent with the reaction of H2 + CH3. The comparison of these three isotopic reactions also shows the isotopic effects in the initial-state-selected reaction probabilities. The cumulative reaction probabilities (CRP) are obtained by summing over initial-state-selected reaction probabilities. The energy-shift approximation to account for the contribution of degrees of freedom missing in the six dimensionality calculation is employed to obtain approximate full-dimensional CRPs. The rate constant comparison shows H2 + CH3 reaction has the biggest reactivity, then HD + CH3, and D2 + CH3 has the smallest.

Wang, Dunyou↗
DOE OSTI2022

Reactive burn model calibration using high-throughput initiation experiments at sub-millimeter length scales

We report a first-of-its-kind model calibration was performed using Sandia National Laboratories’ high-throughput initiation (HTI) experiment for two types of vapor-deposited explosive films consisting of hexanitrostilbene (HNS) or pentaerythritol tetranitrate (PETN). These films exhibit prompt initiation, and they reach steady detonation at sub-millimeter length scales. Following prior work on HNS, we test the hypothesis of approximating these explosive films as fine-grained homogeneous solids with simple Arrhenius kinetics burn models. The model calibration process is described herein using a single-step as well as a two-step Arrhenius rate law, and it consists of systematic parameter sampling leading to a reduction in the model degrees of freedom. Multiple local minima are observed; results are given for seven different optimized parameter sets. Each model set is further evaluated in a two-dimensional simulation of the critical failure thickness for a sustained detonation. Overall, the two-step Arrhenius kinetics model captures the observed behavior for HNS; however, neither model produces a good fit to the PETN data. We hypothesize that the HTI results for PETN correspond to a heterogeneous response, owing to the smaller reaction zone of PETN compared to HNS (i.e., it does not homogenize the fine-grained hot spots as well). Future work should consider using the ignition and growth model for PETN, as well as other reactive burn models such as xHVRB, AWSD, PiSURF, and CREST.

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

Model predictive control of mixing controlled compression ignition operation for low reactivity fuels

Using gasoline or other low reactivity fuels with a pilot injection or port fuel injection in a compression ignition engine has shown great potential in reducing NOx emissions while keeping high thermal efficiency compared to diesel. However, excessive combustion noise is caused by a high maximum pressure rise rate in the cylinder due to the higher fractions of premixed charge of the low-reactivity fuel. This noise can result in structural damage to engine components and as such, combustion noise limits the range of the operating parameters and makes the control of such engines challenging. In this study, a simulation environment was built up in MATLAB/Simulink leveraging a physics-based zero-dimension combustion model to capture the in-cylinder pressure time traces as well as metrics relevant to thermal efficiency and combustion noise. Here, in order to also facilitate the control of emissions, machine learning models were investigated to capture NOx emissions. A kernel-based extreme learning machine (K-ELM) performed best and had a coefficient of correlation (R-squared) of 0.998. The combustion and NOx emission models are valid for not only conventional gasoline fuel but also oxygenated alternative fuel blends at three different pilot injection strategies. In order to track key combustion metrics while keeping noise and emissions within constraints, a model predictive control (MPC) was applied for a compression ignition engine operating with a range of potential fuels and fuel injection strategies. The MPC is validated under different scenarios, including a load step change, fuel type change, and injection strategy change, with proportional–integral (PI) control as the baseline. The simulation results show that MPC reduces about 26% of ringing intensity in the transient process and 17% at the steady state for E30. Generally, MPC can optimize the overall performance through modifying the main injection timing, pilot fuel mass, and exhaust gas recirculation (EGR) fraction.

42 ENGINEERING↗