SEARCH · Engineering Papers

Results for “Kinetic modeling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 775 records · Page 43

All sources NASA NTRS DOE OSTI

DOE OSTI2020

An experimental investigation of flame and autoignition behavior of propane

We report autoignition delay time data are one important means to develop, quantify, and validate fundamental understanding of combustion chemistry at low temperatures (T<1200 K). However, low-temperature chemistry often has higher uncertainties and scatter in the experimental data compared with high-temperature ignition data (T>1200 K). In this study, autoignition properties of propane and oxygen mixtures were investigated using the University of Michigan rapid compression facility in order to understand the effects of ignition regimes on low-temperature ignition data. For the first time for propane, autoignition delay times were determined from pressure histories, and autoignition characteristics were simultaneously recorded using highspeed imaging of the test section through a transparent end-wall. Propane mixtures with fuel-to-O 2 equivalence ratios of φ = 0.25 and φ = 0.5 and O 2 -to-inert gas molar ratios of 1:3.76 were studied over the pressure range of 8.9 to 11.3 atm and the temperature range of 930 – 1070 K. The results showed homogeneous or strong autoignition occurred for all φ = 0.25 experiments, and inhomogeneous or mixed autoignition occurred for all φ = 0.5 experiments. While a limited temperature range is covered in the study, importantly the data span predicted transitions in autoignition behavior, allowing validation of autoignition regime hypotheses. Specifically, the results agree well with strong-autoignition limits proposed based on the Sankaran Criterion. The autoignition delay time data at the strong-ignition conditions are in excellent agreement with predictions using a well-validated detailed reaction mechanism from the literature and a zerodimensional modeling assumption. However, the experimental data at the mixed autoignition conditions were systematically faster than the model predictions, particularly at lower temperatures (T< ~970 K). The results are an important addition to the growing body of data in the literature that show mixed autoignition phenomena are important sources of the higher scatter observed in the low-temperature autoignition data for propane and other fuels. The results are discussed in terms of different methods to capture the effects of pre-autoignition heat release associated with mixed autoignition conditions and thereby address some of the discrepancies between kinetic modeling and experimental measurements.

42 ENGINEERING↗

DOE OSTI2025

Transforming Aeration Energy in Water Resource Recovery Facilities (WRRFs) through Suboxic Nitrogen Removal (Final Report)

The objective of this project was to advance two key technological components—aeration control strategies and process design methodologies—to support the development and broader adoption of suboxic biological nitrogen removal (SBNR). The project focused on achieving the following three goals: • Enhance Model Predictive Control (MPC) Technology: Advance the DO/Nmaster MPC platform from its initial 2018 pilot deployment at the Chico Water Resource Recovery Facility in California to full-scale integration. This included partnering with a blower technology commercialization partner and incorporating machine learning (ML) capabilities to enable nationwide deployment. • Bridge Knowledge Gaps in SBNR Process Design: Address fundamental gaps in SBNR process understanding through controlled pilot-scale testing at a dedicated pilot facility. These efforts supported the development of robust kinetic models to inform reliable SBNR control, operational strategies, and design frameworks. • Demonstrate Full-Scale Implementation of Low DO/SBNR with ML: Transition low dissolved oxygen (DO)/SBNR coupled with ML from pilot-scale trials to full-scale demonstration in flow-through biological nutrient removal (BNR) systems, with the goal of enabling scalable, nationwide adoption in activated sludge treatment processes. The project included demonstration of SBNR at the pilot scale as performed by Hampton Roads Sanitation District (HRSD) and at the full-scale as performed by the Los Angeles County Sanitation Districts' (LACSD) Pomona Water Reclamation Plant (POWRP).

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1989

On the driving force of PAH production

The kinetic factors affecting the production of polycyclic aromatic hydrocarbons (PAH) in high-temperature pyrolysis and combustion environments are analyzed. A lumped kinetic model representing polymerization-type growth by one irreversible step and two reversible steps is considered. It is shown that at high temperatures, PAH growth is controlled by the superequilibrium of hydrogen atoms; at low temperatures and low H2 concentrations, the PAH growth rate is proportional to the rate of the H-abstraction of a hydrogen atom from aromatic molecules; while at low temperatures and high H2 concentrations, it is controlled by the thermodynamics of the H-abstraction and the kinetics of acetylene addition to aromatic radicals. The presence of oxygen mainly affects the small-molecule reactions during the induction period.

Frenklach, Michael↗

DOE OSTI2021

Non-Boltzmann Effects in Chain Branching and Pathway Branching for Diethyl Ether Oxidation

Low-temperature (LT) engine applications have several potential benefits, including reduced emissions and increased efficiency. Attaining these benefits requires accurate kinetic modeling of LT chain branching, which depends heavily on ketohydroperoxide (KHP) decomposition. For diethyl ether (DEE), a promising biofuel, current estimates of the KHP decomposition rate constant are largely based on empirical fits to data. In this study, we investigate the most important KHP isomer in DEE LT oxidation by applying variable reaction coordinate transition state theory to the main pathway for KHP decomposition: OO bond fission to produce •OH and a keto-alkoxy radical, •OQ'O. We also use ab initio kinetics methods to investigate the decomposition of •OQ'O, where we find dominant branching to acetic acid, with the remaining flux going to CH 3 C(O)OCHO. Additionally, new time-resolved measurements of DEE and acetic acid concentrations during LT (450–600 K) DEE oxidation are obtained in a laser photolysis flow reactor coupled with multiplexed photoionization mass spectrometry. These new experimental data, along with jet-stirred reactor data in the literature, are compared with the predictions of a recent DEE mechanism (Tran et al. Proc. Comb. Inst. 2019, 37, 511-519) that was modified with the newly calculated ab initio rate constants for KHP and •OQ'O decomposition. The predictions of the modified mechanism are quite poor when compared to the experimental data; this is primarily due to the new KHP ⇌ •OQ'O + •OH rate constant, which is 1–2 orders of magnitude slower than empirical values employed in recent mechanisms. To reconcile the new KHP rate constant and the experimental data, we explore and quantify the possible role of non-Boltzmann (nB) reaction sequences. The nB reactions have a substantial effect on both the overall mechanism reactivity and the •OQ'O branching to acetic acid. We also provide guidance on the proper implementation of nB reactions in kinetic mechanisms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Corrosion-Resistant Non-Carbon Electrocatalyst Supports for PEFCS

Insufficient durability of commercial catalysts is a persistent issue for PEFC applications.1 Herein, Pt on Nb-doped-TiO 2 is reported as a highly durable electrocatalyst on which the direct 4-electron reduction of oxygen to water is more facile compared to Pt/C. All performance metrics are reported comparing Pt/Nb-TiO2 with a 15% Pt loading, against a commercial Pt/C catalyst (46.5% Pt loading, Tanaka, K. K.). Nb-doped-TiO 2 with high surface area and high electronic conductivity was synthesized using the supercritical fluid method. Initially, the durability of the catalyst was characterized using accelerated stability tests (ASTs) involving 10,000 high potential cycles (DOE/FCCJ protocol) and the Pt/Nb-TiO 2 was found to retain 78% of its initial electrochemically active surface area (ECSA) compared to the 57.6 % retained by Pt/C. These observations were in excellent agreement with previous reports that the Pt particle size of Pt/C grew from 2nm to 8nm during the AST protocol along with severe corrosion and amorphization of the carbon surface.2, 3 TEM and XPS studies of the Pt-Nb-doped-TiO 2 catalyst showed that the Pt particle size grew from 3nm to 6nm and the Nb(IV) and Ti(III) in the support were oxidized to Nb(V) and Ti(IV) after the durability test. Thus, the improvement in Pt/Nb-TiO 2 ECSA retention was attributed to the lower extent of particle growth and lack of oxidative support loss upon oxidation as compared to Pt/C. The oxygen reduction reaction (ORR) performance was characterized by linear polarization using a rotating disk electrode (RDE). The electrochemical surface areas of Pt/Nb-TiO 2 and Pt/C were found to be 48m2·g -1 and 83m2·g -1 respectively, and the mass activity for the ORR at 0.9V vs. RHE were found to be 150 mA·mg-1Pt and 124 mA·mg-1Pt respectively. The improved mass activity on Pt/Nb-TiO 2 was attributed to strong metal support interaction (SMSI) between the Nb-TiO 2 support and the Pt catalyst based on the 625 meV decrease in the binding energy of the Pt4f x-ray photo-electron spectroscopy (XPS) peaks of Pt/Nb-TiO 2 compared to Pt/C. To quantify the impact of the SMSI, a kinetic model was applied to calculate the elementary reaction rate constants for the various steps of the ORR on both catalysts. The reaction rate constant (k1) for the direct 4-electron transfer pathway to produce H 2 O was significantly larger in Pt/Nb-TiO 2 as compared to Pt/C. Thus, the reduction in the electron binding energy as observed in the XPS was found to aid the facile filling of the higher energy 2p orbitals of O 2 thereby leading to improved 4-electron transfer kinetics and improved overall activity.

08 HYDROGEN↗

DOE OSTI2022

Simulation of the photodetachment spectra of the nitrate anion (NO 3 ₋ ) in the B ~ 2 E' energy range and non-adiabatic electronic population dynamics of NO 3

The photodetachment spectrum of the nitrate anion (NO 3 ₋ ) in the energy range of the NO 3 second excited state is simulated from first principles using quantum wave packet dynamics. The prediction at 10 K and 435 K relies on the use of an accurate full-dimensional fully coupled five state diabatic potential model utilizing an artificial neural network. The ability of this model to reproduce experimental spectra was demonstrated recently for the lower energy range. Analysis of the spectra indicates a weaker Jahn–Teller coupling compared to the first excited state. The detailed non-adiabatic dynamics is studied by computing the population dynamics. An ultra-fast non-statistical radiationless decay is found only among the Jahn–Teller components, which is followed by a slow statistical non-radiative decay among the different state manifolds. The latter is reproduced perfectly by a simple first order kinetics model. The dynamics in the second excited state is not affected by the presence of a conical intersection with the first excited state manifold.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Electron Energy Partition across Interplanetary Shocks. III. Analysis

An analysis of model fit results of 15,210 electron velocity distribution functions (VDFs), observed within ±2 hr of 52 interplanetary (IP) shocks by the Wind spacecraft near 1 au, is presented as the third and final part on electron VDFs near IP shocks. The core electrons and protons dominate in the magnitude and change in the partialto-total thermal pressure ratio, with the core electrons often gaining as much or more than the protons. Only a moderate positive correlation is observed between the electron temperature and the kinetic energy change across the shock, while weaker, if any, correlations were found with any other macroscopic shock parameter. No VDF parameter correlated with the shock normal angle. The electron VDF evolves from a narrowly peaked core with flaring suprathermal tails in the upstream to either a slightly hotter core with steeper tails or much hotter flattop core with even steeper tails downstream of the weaker and strongest shocks, respectively. Both quasi-static and fluctuating fields are examined as possible mechanisms modifying the VDF, but neither is sufficient alone. For instance, flattop VDFs can be generated by nonlinear ion acoustic wave stochastic acceleration (i.e., inelastic collisions), while other work suggested they result from the combination of quasi-static and fluctuating fields. This three-part study shows that not only are these systems not thermodynamic in nature; even kinetic models may require modification to include things like inelastic collision operators to properly model electron VDF evolution across shocks or in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2020

Electron Energy Partition Across Interplanetary Shocks. III. Analysis

Solar wind↗

DOE OSTI2023

Low-temperature oxidation pathways are critical to thermal incineration of PFAS-laden materials

With growing desire to destroy per- and poly-fluoroalkyl substances (PFAS) now known to be detrimental to human health, a sound understanding of fluorocarbon combustion chemistry is important to efficient thermal destruction within incinerators. While most fluorocarbon combustion models and the sets of reactions contained within them were originally developed for the high temperatures encountered in flame suppression applications, they have often been used to assess PFAS destruction in incinerators, which emphasize a lower range of temperatures. We present results that demonstrate that low-temperature fluorocarbon oxidation pathways—not yet known to play a role in fluorocarbon combustion—impact key incinerator performance metrics, including: PFAS surrogate mole fractions, products of incomplete destruction, and waste destruction efficiencies. The results further point to the utility of NO as a potential additive. The present results show the influence of these pathways for CF 3 O 2 , for which some data are available, but analogous pathways would also occur for other fluoroalkylperoxy radicals, for which little is known. Furthermore, the results demonstrate the need for future work to identify and characterize low-temperature pathways more generally, consider such pathways in kinetic model development, and experimentally probe intermediate temperature conditions to better understand, design, and control thermal destruction technologies for improved PFAS management.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2012

Role of Transport and Kinetics in Growth of Renal Stones

Renal stone disease is not only a concern on earth but could conceivably pose as a serious risk to the astronauts health and safety in Space. In this paper, a combined transport-kinetics model for growth of calcium oxalate crystals is presented. The model is used to parametrically investigate the growth of renal calculi in urine with a focus on the coupled effects of transport and surface reaction on the ionic concentrations at the surface of the crystal and their impact on the resulting growth rates. It is shown that under nominal conditions of low solution supersaturation and low Damkohler number that typically exist on Earth, the surface concentrations of calcium and oxalate approach their bulk solution values in the urine and the growth rate is most likely limited by the surface reaction kinetics. But for higher solution supersaturations and larger Damkohler numbers that may be prevalent in the microgravity environment of Space, the calcium and oxalate surface concentrations tend to shift more towards their equilibrium or saturation values and thus the growth process may be limited by the transport through the medium. Furthermore, parametric numerical studies suggest that changes to the renal biochemistry of astronauts due in space may promote development of renal calculi during long duration space expeditions.

Kassemi, Mohammad↗

NASA NTRS1974

Development of a global pollution model for CO, CH4, and CH2O

The current status of a global pollution model for carbon monoxide, methane, and formaldehyde is described. The physico-chemical action is considered of these three pollutants in the troposphere. This geographic restriction is convenient since the tropopause provides a natural boundary across which little transport occurs. The data on sources and sinks for these pollutants is based on available information and assumptions relative to the major man-made and natural contributions. The distributions and concentrations of methane, formaldehyde, and carbon monoxide in the atmosphere are interrelated by the chemical reactions in which they participate. A chemical kinetic model based on the pseudo-steady state approximation for the intermediate species was developed to account for these reactions. The numerical procedure used to mathematically describe the pollution transport is a mass conservative scheme employing an integral flux approach.

Peters, L. K.↗

DOE OSTI2021

Ultrafast conformational dynamics of Rydberg-excited N -methyl piperidine

We have observed the ultrafast conformational dynamics of electronically excited N-methyl piperidine (NMP) using time-resolved Rydberg fingerprint spectroscopy. Optical excitation at various wavelengths ranging from 212 nm to 229 nm leads to the 3s or 3p Rydberg states and induces coherent oscillatory motions with periods of about 700 fs. Furthermore, these coherent motions survive the internal conversion from 3p to 3s but then dephase on a time scale of a few oscillations. Intramolecular vibrational energy redistribution on a picosecond time scale leads to an equilibrium between two conformeric structures that are separated in binding energy by 0.09 eV. Model calculations using the DFT-SIC method are in excellent agreement with the experiments and identify the conformers as the chair and twist structures of NMP. The analysis of the equilibrium parameters at long time delays as a function of excitation wavelength allows for the extraction of thermodynamic parameters for the conformeric transformation. We derive an enthalpy of the chair to twist reaction in the 3s excited state of 62 meV with an entropy of 19.70 J mol –1 K –1 . An activation energy of 276 meV is also obtained with a kinetic model.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2021

Preliminary Thermoluminescent Dosimeter Glow Curve Analysis with Automated Glow Peak Identification for LiF Mg,Ti

When appropriately analyzed, thermoluminescent dosimeter glow curve analysis allows for improved quantification of thermoluminescent material behavior while flagging abnormalities. The mathematical separation of a glow curve into contributions from energetically unique trap states, or glow curve analysis, may be used to remove undesired effects of signal fading for complex materials. A generalized glow curve analysis software for the separation of glow curves is presented in this paper. Written in C ++ , the software uses the first-order kinetics model with automatic peak identification. The automatic identification of peaks is achieved through a unique peak-finding algorithm. Here, the program was performance tested using experimental glow curve data from LiF:Mg,Ti, and comparative results are presented.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Deterministic modelling of the SPERT IV reactor transients using the multi-physics capability of WIMS

The ANSWERS WIMS reactor physics code is being developed for whole core multi-physics modelling. The established neutronics capability for lattice calculations has recently been extended to be suitable for whole core modelling of Small Modular Reactors (SMRs). A whole core transport, SP3 or diffusion flux solution is combined with fuel assembly resonance shielding. An integrated thermal hydraulic solver permits temperature and density variations to feedback to the neutronics calculation. This capability can be applied to both steady state and time dependent transient problems. Nuclear reactor design and safety case development requires assessment of a range of reactor transients, to inform both normal operation limitations and accident scenario analysis. This paper presents new methodology developed in WIMS to couple the core neutronics to the integrated core thermal hydraulics solver for the simulation of reactor transients in whole core models. To support the validation of the multi-physics capability of WIMS, this capability has been applied to the reactivity insertion transient experiments performed in the plate fuelled SPERT-IV reactor. This study employs WIMS using the whole core solver MERLIN, which calculates the time-dependent flux distribution and magnitude, coupled to ARTHUR, which solves for the time-dependent thermal hydraulics solution, and provides thermal feedback via the cross sections generated by GEOM, which performs resonance shielding calculations and generates cross section data for the plate geometry. This spatial kinetics model, with dynamic cross section generation, allows for variations in the temperature and neutron flux profile with time.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRS2001

Production, Outflow Velocity, and Radial Distribution of H2O and OH in the Coma of Comet C/1995 O1 [Hale-Bopp] from Wide Field Imaging of OH

Observations of OH are a useful proxy of the water production rate (Q(sub H2O)) and outflow velocity (V(sub out)) in comets. We use wide field images taken on 03/28/1997 and 04/08/1997 that capture the entire scale length of the OH coma of comet C/1995O1 (Hale-Bopp) to obtain Q(sub H2O) from the model-independent method of aperture summation. We also extract the radial brightness profile of OH 3080 angstroms out to cometocentric distances of up to 10(exp 6) km using an adaptive ring summation algorithm. Radial profiles are obtained as azimuthal averages and in quadrants covering different position angles relative to the comet-Sun line. These profiles are fit using both fixed and variable velocity two-component spherical expansion models to determine VOH with increasing distance from the nucleus. The OH coma of Hale-Bopp was more spatially extended than in previous comets, and this extension is best matched by a variable acceleration of H2O and OH that acted across the entire coma, but was strongest within 1-2 x 10(exp 4) km from the nucleus. This acceleration led to VOH at 10(exp 6) km that was 2-3 times greater than that obtained from a 1P/Halleytype comet at 1 AU, a result that is consistent with gas-kinetic models, extrapolation from previous observations of OH in comets with Q(sub H2O) > 10(exp 29)/s, and radio measurements of the outer coma Hale-Bopp OH velocity profile. When the coma is broken down by quadrant, we find an azimuthal asymmetry in the radial distribution that is characterized by an increase in the spatial extent of OH in the region between the orbit-trailing and anti-sunward directions. Model fits to this area and comparison with radio OH measurements suggest greater acceleration in this region, with VOH UP to 1.5 times greater at 10(exp 6) km radial distance than elsewhere in the coma.

Harris, Walter M.↗

DOE OSTI2024

Intrinsic rotation modulation by diffusive neutral particles in tokamaks

Abstract The modulated transport model, a model kinetic ion transport equation for the pedestal and scrape-off layer (SOL), is generalized to self-consistently include effects of a single neutral particle species. The neutrals contribute additional transport terms, modifying the v ∥ -dependent orbit-averaged ion diffusivities of the original work and the resulting predicted intrinsic rotation of the ions. After making simplifying assumptions of the neutral transport, in particular taking the continuous transport limit via a short charge-exchange step expansion, we derive relatively simple analytic expressions that capture the diffusive neutral physics. Within the scope of the model’s validity, the neutral-driven intrinsic rotation can compete with the turbulence-driven intrinsic rotation. However, for physically motivated parameters, the neutral-driven intrinsic rotation appears negligible, either on a term-by-term basis or due to a strong cancellation between the neutral-driven momentum diffusion and pinch terms. It appears that a treatment containing finite charge-exchange steps is necessary to capture neutral transport of strong flow momentum into the confined region from the SOL.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Control of solvent adsorption in porous liquids through solvent-solvent interactions

Type 3 Porous Liquids (PLs) are a new class of materials that offer transformative potential for gas capture and utilization. These PLs are comprised of a bulky solvent with a suspended empty nanoporous material that enables selective and high-capacity capture from dilute and complex gas steams. The relationship between the nanoporous material, the structure of the solvent molecules, and the time dependence of solvent adsorption into the nanoporous material governs long term adsorption properties. Herein, molecular dynamics (MD) calculations evaluated the solvent dynamics of nine neat solvents with a broad range of chemical identities and experimental density measurements probed the time-dependent adsorption of these solvents into the ZIF-8 metal-organic frameworks (MOFs) to form PL dispersions over three weeks. Two of the nine dispersion, using glyceryl triacetate and 2’-hydroxyacetopheneone as solvents, presented sufficiently low solvent infiltration—characterized by solvent adsorption less than 40% of the ZIF-8 pore volume—to be viable as PLs. The experimental solvent adsorption data for ZIF-8 dispersed in the four aromatic solvents (acetophenone, methylbenzoate, 2-isopropylphenol, and 2’-hydroxyacetophenone) was fit to a kinetic model to quantify rate of solvent adsorption into ZIF-8. The rates of adsorption of these four similar-sized solvents demonstrated that solvent adsorption is slowed and limited by solvent clustering dynamics that are, in turn, driven by intermolecular hydrogen bonding. Ultimately, experimental solvent sorption measurements indicated that glyceryl triacetate and 2’-hydroxyacetophenone with ZIF-8 formed PLs with the most stable porosity and the greatest potential for gas capture. In conclusion, the combined experimental-computational materials exploration framework used here reveals a novel relationship between molecular scale solvent dynamics and macroscale solvent adsorption kinetics critical for the discovery and rapid evaluation of Type 3 PL compositions.

Hydrogen bonding↗

DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗