Engineering Papers⌕ Search

Engineering topics

Kim, Seonah

Publications and source records attributed to Kim, Seonah.

High-Pressure Rate Rules for Ether Alkylperoxy Radical Isomerization

The first isomerization reaction of an alkylperoxy (RO 2 ) radical holds significant importance in low-temperature oxidation, as it governs the branching ratios of the hydroperoxyalkyl (QOOH) radicals, which influence the competition between the chain-propagation and chain-branching reactions. In this study, we systematically calculated high-pressure rate rules for the RO 2 isomerization reaction of monoethers, exploring 5-, 6-, 7-, and 8-membered ring transition states. Primary, secondary, and tertiary carbon sites, where both the abstracting peroxy group and the abstracted hydrogen are located, were considered, with particular emphasis on distinguishing between secondary carbons adjacent (alpha) and nonadjacent to the ether functional group. Using the G4//B3LYP/6-311++G(2df,2pd) level of theory and the transition state theory, we estimated the rate constants and the Arrhenius coefficient for over 120 possible isomerization reactions. We examined the effect of ring size and ring atoms, revealing that 6- and 7-membered ring isomerizations were generally the fastest. The impact of the ether functional group on transition states was investigated by comparing reactions with identical ring size, peroxy, and radical positions, but with the ether functional group positioned either outside (i.e., out) or inside (i.e., in) the transition state ring, leading to differences in the rate constants. When comparing to analogous alkane rate constants, differences of up to an order of magnitude were observed, underscoring the need for caution when assigning rate rules by analogy. We applied our rate constants in the di-iso-butyl ether kinetic model and evaluated their influence on low-temperature chemistry finding that they altered the branching ratios by up to a factor of 9, highlighting the significance of site-specific rate constants for more accurate low-temperature modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fuel property evaluation of unique fatty acid methyl esters containing β-hydroxy esters from engineered microorganisms

Unique fatty acid methyl esters (FAME) containing ..beta..-hydroxy esters were produced using an engineered microorganism by glucose fermentation. This study investigated the properties of the unique FAME mixture both neat and in blends with conventional diesel, as well as properties of ..beta..-hydroxy esters. The unique FAME blend contained relatively shorter-chain FAME (average fatty acid chain carbon number 14.6) with 58 % monounsaturated fatty acids and 9 % saturated and monounsaturated ..beta..-hydroxy acid chains. The unique FAME had significantly lower distillation T90 (321 °C versus 352 °C) and higher cetane number (56.7 versus 52) compared to soy biodiesel. Cloud points were within method repeatability. Unexpectedly (because of the lack of methylene-interrupted double bonds), the unique FAME had low oxidation stability (1.5 h) as determined by Rancimat induction period. Stability could be improved through addition of commonly used antioxidants. We speculate that monounsaturated ..beta..-hydroxy FAME may be the source of this instability. Blends with conventional diesel up to 50 vol% showed similar kinematic viscosity (within method repeatability) as blends of conventional FAME. The unique FAME had no effect on distillation T90 even at the 80% blend level. A 30 vol% blend into conventional diesel had a Rancimat induction period of only 2 h, very nearly the same as the neat unique FAME sample. The addition of antioxidants produced blends of acceptable stability. Based on an assessment of the properties of individual ..beta..-hydroxy FAME molecules, they have higher boiling point, higher cloud point, lower cetane number, and potentially lower storage stability than analogous FAME not having the ..beta..-hydroxy group. Removing them from the fuel product in the production process may result in a biodiesel product with superior properties to what is on the market today.

09 BIOMASS FUELS↗

A Machine Learning Model for Predicting Composition of Catalytic Coprocessing Products from Molecular Beam Mass Spectra

Demand for the development of an automated and integrated refining process for biofuels has increased in recent years due to the lack of generalized process inspection tools. In bio-oil upgrading processes, all process variables are maintained based on the offline specification of intermediates and products. A lack of real-time product specifications in batch-wise monitoring can cause process failure and wasted resources. Therefore, there is a need for a fast and accurate intermediates/product specification tool that can be used for real-time specification to reduce waste and mitigate the risk of process failure. Here, to address this gap, we developed a machine learning (ML) model for predicting speciated bio-oil composition, including paraffin, iso-paraffins, olefins, naphthene, and aromatics. The model is trained using the mass spectra from upgraded products collected in the vapor phase before condensation and predicts the composition of the condensed product. Training ML models using raw mass spectra is challenging due to numerous overlapped peaks originating from different parent compounds. With this in mind, we propose a protocol that (i) transforms raw mass spectra to chemistry-inspired predefined features and (ii) trains decision tree-based models using these features. Our results show that the random forest model was robust against overfitting and had the highest accuracy compared to other models. Moreover, a stochastic ablation method determined the eight most significant features while maximizing the accuracy. Our protocol facilitates real-time compositional analysis of upgraded bio-oils and thus real-time process monitoring. Additionally, this protocol enables the rational design of efficient catalysts and the determination of optimal process conditions.

09 BIOMASS FUELS↗

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗

Physics-informed graph neural networks for predicting cetane number with systematic data quality analysis

Designing alternative fuels for advanced compression ignition engines necessitates a predictive model for cetane number (CN). In this study, the physics-informed graph neural networks are introduced for a reliable CN prediction by considering molecular features pertinent to the physical properties of molecules that affect CN. The reliability of measured data is another key factor to consider for improving the predictive model. Various experimental instruments for measuring CN exist, including standard and non-standard methods. In this regard, a systematic data quality analysis was carried out for the total 630 CNs collected from literature and new measurements in this study using Advanced Fuel Ignition Delay Analyzer (AFIDA). The results from this data curation process were reflected in the model by imposing lower sample weights on the data coming from less reliable measurement techniques. This approach effectively maximized the prediction accuracy while incorporating data from all available sources. Using the sample weights decreased the mean absolute error (MAE) up to 0.8 CN units. The accuracy was also improved by introducing the CN-related physical properties (the number of hydrogen bond donors and acceptors); the test set MAE is 5.74 and 7.01 for the model with and without such properties, respectively. Investigating molecular structural effects on CN was also carried out to gain chemical insights into factors used to design new fuel candidates. The dimensionality reduction analysis of feature vectors showed a clear clustering in terms of functional groups and CN and the structural effect derived from the model was consistent with the physicochemical insights. Finally, this physics-informed model and data curation would be helpful for accurate CN prediction and inform rational fuel design.

97 MATHEMATICS AND COMPUTING↗

Enhancing $\phi$-sensitivity of ignition delay times through dilution of fuel-air mixture

The high $\phi$-sensitivity (η) of ignition delay time (τ IDT ) is one of the desirable fuel properties for the high-load extension of advanced compression ignition engine. Recent studies revealed the effectiveness of a high dilution rate (x D ) for enhancing η at the engine-relevant conditions. This study aims to quantify the effect of dilution on the η of isooctane using a combined experiment-simulation approach. The ignition delay of the isooctane/air mixture was measured with an Advanced Fuel Ignition Delay Analyzer (AFIDA) over the temperature range of 623 - 923 K at 10 bar pressure with global $\phi$ = 0.3 - 0.6, with and without 28.6% of additional N 2 dilution. For precise evaluation of experimental η, the facility effect of the AFIDA experiment was characterized with three-dimensional computational fluid dynamics (3-D CFD) simulation. The temperature in the combustion chamber from 3-D CFD indicated a substantial temporal dependency, varying up to ~52 K by charge-cooling of fuel injection and heat transfer from the wall. We introduced the dimensionless number θ(t) for characterizing the temporal profile of chamber temperature. Consideration of the facility effect using θ(t) resulted in better agreement between the experimental ..eta.. and zero-dimensional (0-D) kinetics simulation. The refined η were then further utilized to quantify the effectiveness of dilution to η The extent of η enhancement with dilution strategy was maximized at the low-temperature chemistry regime, increasing η by 77% with a 28.6% dilution rate. Further analysis on the dilution effect was carried out using 0-D kinetics simulation, revealing the critical dimensionless numbers relevant to the effectiveness of dilution to η enhancement. Here, this study is the first experiment-simulation combined research to quantify the effect of dilution on η, facilitating the kinetics model refinement for better reproduction of $\phi$-sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗

Sooting tendencies of terpenes and hydrogenated terpenes as sustainable transportation biofuels

Terpenes are a diverse group of molecules that are synthesized by plants and microorganisms through combining units of isoprene (2-methyl-1,3-butadiene). They typically contain rings and methyl branches, which gives them high energy densities and low freezing points and makes them appealing candidates for sustainable transportation biofuels. Between the original biosynthesis and upgrading options such as hydrogenation, they have a large degree of freedom of structures, e.g., different carbon skeletons, positions of double bonds, and functional groups. Therefore, structure-property data is needed to downselect potential fuel candidates. Here, we measured the sooting tendencies of 17 C10 monoterpenes and 7 of their hydrogenated analogues. The hydrogenated compounds were custom synthesized, so the quantities were too small for conventional smoke point measurements. Thus, the sooting tendencies were quantified with yield sooting index (YSI), which is based on the soot yield in a fuel-doped non-premixed methane flame. Derived smoke points (DSPs) were estimated from a correlation between YSI and smoke point for other hydrocarbons. The YSI of terpenes and their derivatives varies widely from 85.6 to 248.5. The YSI follows the trend: terpenes > dihydroterpenes > tetrahydroterpenes. The DSPs of all the tetrahydroterpenes and some dihydroterpenes are higher than that of a Jet-A fuel sample, suggesting that they offer soot reduction benefits. Further, the YSIs depend strongly on molecular structure; for example, α-pinene and β-pinene have identical carbon skeletons and differ only in the position of one carbon-carbon double bond, but the YSI of α-pinene is 34% higher than that of β-pinene. Detailed decomposition analysis via density functional theory (DFT) suggests that compared with β-pinene, α-pinene requires fewer steps to form the first aromatic ring and the process is more thermodynamically favorable. The YSI difference between the pinenes is mainly affected by the identity of the products from the dominant decomposition pathways.

09 BIOMASS FUELS↗

Connecting cation site location to alkane dehydrogenation activity in Ni/BEA catalysts

Ni-modified beta zeolite (Ni/BEA) catalysts activate carbon-hydrogen bonds in light alkanes, as demonstrated through isobutane reaction testing. Controlled synthesis of Ni/BEA allows for efficient introduction of ion-exchanged Ni sites at varying Ni loadings (0.43% - 1.8%). These catalysts exhibit site time yields (STY) for H2 production that increase with increasing Ni loading. A detailed analysis of secondary reactions and carbon deposition based on the relative molar flowrates of product C and H indicates that the observed increase in H2 STY with increasing Ni loading is likely attributed to both increasing alkane activation activity and increasing formation of hydrogen-deficient aromatic products retained within the catalyst pores. In situ diffuse-reflectance UV-visible-NIR absorbance and X-ray absorption spectroscopies indicate isolated, 4-coordinate Ni(2+) species for all loadings. Quantum mechanics/molecular mechanics modeling identifies two distinct Ni(2+) sites consistent with the structural characterization, but with differing relative stabilities due to their coordination environment. Computed reaction energetics for isobutane dehydrogenation demonstrate that the more stable Ni(2+) species at a six-membered 4Si-2Al ring, Ni-6MR, is less active for isobutane dehydrogenation than the less stable Ni(2+) at the five-membered 3Si-2Al ring, Ni-5MR. The differing local structure of the isolated cationic Ni sites in Ni/BEA offers a possible rationalization for the increased H2 STY observed at greater Ni loadings.

catalysts↗

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 1

We report di-isobutylene has received significant attention as a promising fuel blendstock, as it can be synthesized via biological routes and is a short-listed molecule from the Co-Optima initiative. Di-isobutylene is also popularly used as an alkene representative in multi-component surrogate models for engine studies of gasoline fuels. However, there is limited experimental data available in the literature for neat di-isobutylene under engine-like conditions. Hence, most existing di-isobutylene models have not been extensively validated, particularly at lower temperatures (< 1000 K). Most gasoline surrogate models include the di-isobutylene sub-mechanism published by Metcalfe et al. with little or no modification. The current study is undertaken to develop a detailed kinetic model for di-isobutylene and validate the model using a wide range of relevant experimental data. Part 1 of this study exclusively focuses on the low- to intermediate temperature kinetics of di-isobutylene. An upcoming part 2 discusses the high-temperature model development and validation of the relevant experimental targets. Ignition delay time measurements for the di-isobutylene isomers were performed at pressures ranging from 15 - 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in air and in the temperature range 650 - 900 K using two independent rapid compression machine facilities. In addition, measurements of species identified during the oxidation of these isomers were performed in a jet-stirred reactor and in a rapid compression machine. A detailed kinetic model for the di-isobutylene isomers is developed to capture the wide range of new experimental targets. For the first time, a comprehensive low-temperature chemistry submodel is included. The differences in the important reaction pathways for the accurate prediction of the oxidation of the two DIB isomers are compared using reaction path analysis. In conclusion, the most sensitive reactions controlling the ignition delay times of the DIB isomers under the pressure and temperature conditions necessary for autoignition in engines are identified.

09 BIOMASS FUELS↗

Predicting Catalytic Pyrolysis Aromatic Selectivity from Pyrolysis Vapor Composition Using Mass Spectra Coupled with Statistical Analysis

The behavior of fast pyrolysis (FP) and catalytic FP (CFP) of 20 renewable feedstocks was studied in a microscale reactor with molecular beam mass spectral analysis of products generated. A partial least-squares (PLS) model was constructed based on the FP vapor spectra that predicts the aromatic selectivity when upgrading over a ZSM-5 catalyst. Additionally, principal component analysis of both FP and CFP spectra was performed for comprehensive spectral analysis. This work highlighted the value of vapor-phase mass spectral screening to predict the subsequent feedstock performance and demonstrated that the quantity of coke deposited on the catalyst is not a reliable measure of catalyst deactivation when the feedstock type is varied.

09 BIOMASS FUELS↗