Engineering topics
Publications and source records attributed to Kim, Yeonjoon.
Computational predictions of vaporization properties aid the de novo design of green chemicals, including clean alternative fuels, working fluids for efficient thermal energy recovery, and polymers that are easily degradable and recyclable.
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.
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.
Bond dissociation energetics underpin the thermodynamics of chemical transformations where bonds are broken or formed and can also be used to predict reaction rates and selectivities.
One of the major biofuels used today is biodiesel, composed of fatty acid methyl esters. The fats and oils feedstocks used to make biodiesel are in relatively short supply. To meet increasing global demand, engineered microorganisms have been developed to catalyze conversion of sugar to fatty acid esters that contain unique ..beta..-hydroxy-esters. This study investigated the effect of a hydroxyl group on the combustion characteristics of methyl esters by comparing the chemical behavior of methyl hexanoate (MHx) and methyl 3-hydroxyhexanoate (M3OHHx) - used as surrogates for diesel-boiling-range esters with longer fatty acid chains. The oxidation of these esters was studied experimentally in a flow reactor at 0.84 and 10 bar, 600 to 1,100 K, and stoichiometric conditions; and in a constant volume combustion chamber at 5 and 10 bar, 600 to 900 K, and equivalence ratios of 0.3 and 0.6. MHx was more reactive in the constant volume chamber at temperatures below 800 K at 10 bar and equivalence ratio of 0.6. MHx also exhibited a higher indicated cetane number (16.4) than M3OHHx (8.1). We investigated this reactivity trend using an updated MHx kinetic model and M3OHHx model developed as part of this work. The kinetic models predicted that radicals formed from MHx were consumed by low-temperature chemistry reactions, whereas the M3OHHx reactivity was significantly governed by the ..beta..-radical (on the same carbon as the OH group) chemistry which mainly terminated via the less reactive chain propagation pathway to methyl-3-oxohexanoate + HO2. This study extends our understanding of structural effects on ester reactivity, which will allow for accurate surrogate formulation and performance simulations.
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.
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.
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.
Biofuels present opportunities for improving the performance and reducing emissions from internal combustion engines by incorporating oxygenated functional groups to the fuels. Among various oxygenates, ethers have been recognized as promising candidates for an alternative to conventional diesel fuel owing to their higher reactivity and lower sooting tendency. The detailed guidelines for designing ethers, however, have not been fully discussed, even though their combustion characteristics are sensitive to the molecular structure. This study was devoted to exploring the structure-property relationships, particularly focusing on the cetane number and yield sooting index, using five linear, branched, and cyclic ethers: di-amyl-ether, 4-butoxy-heptane, 3,3-dimethyl-oxetane, 2-ethyl-4-methyl-1,3-dioxolane, and 2-isopropyl-4-methyl-1,3-dioxolane. First, we examined the chemical kinetics underlying the sooting tendency of the test fuels. The combustion product distribution was measured from flow reactor experiments at 750-1100 K, F=3, at atmospheric pressure. As a result, it was revealed that the sooting tendency is closely related to the size of hydrocarbon intermediates in the high-temperature regime (>1000 K); that is, larger hydrocarbons lead to more soot precursor formation. The underlying chemistry determining the size of the hydrocarbons from the tested fuels was analyzed using reaction pathway analysis and quantum mechanics calculations, which showed that the branched and cyclic ether structures form abundant C3-C4 compounds. Moreover, the auto-ignition characteristics of the test fuels were studied using the flow reactor at low-temperature (400-700 K) and F=1. We found a clear difference in the combustion-product distribution from high and low reactivity fuels, which was then correlated to the systematic analysis of the key reaction energy barriers with the varying molecular structure.
This study developed the design principle of bio-derived ether fuels with low soot emission and high reactivity based on an experiment-theory combined approach.
In sustainable chemical engineering, catalytic upgrading of lignocellulosic biomass has recently gained attention for producing renewable platform chemicals. To achieve maximal biomass utilization, upgrading the underutilized lignin components is essential. Among various catalysts for lignin upgrading, supported vanadia (V2O5) catalysts are promising because of their cost-effectiveness and tunability of either dopant metals or catalyst supports. Here, computational studies are conducted to derive rational design guidelines of supported V2O5 catalysts for accomplishing the high catalytic activity of lignin upgrading to phenol, a key compound for producing bioplastics and biofuel blendstocks. Guaiacol was used as the model compound since it comprises the highest portion of depolymerized lignin. Computational mechanistic studies for the catalytic guaiacol conversion to phenol were performed for the V2O5 catalysts on Titania (TiO2) and silica (SiO2) to explain higher experimental phenol yields on V2O5/SiO2 than V2O5/TiO2. The hydrogen migration from the methoxy group to the aryl ring was identified as a rate-determining step, and the overall activation energies on the two catalysts were compared. A structural analysis was carried out for the catalysts and rate-determining transition states to gain further insights from mechanistic studies. It was concluded that the tilt angle of the aryl group in the hydrogen migration transition state is a key descriptor determining the catalytic activity of phenol formation. These features correlate well with activation energies and experimental phenol yields, indicating that they provide design guidelines for supported metal catalysts for lignin upgrading before experiments.
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.