Next Article in Journal
Grid-Following Mode Operation of Small-Scale Distributed Battery Energy Storages for Fast Frequency Regulation in a Mixed-Source Microgrid
Next Article in Special Issue
Influence of Oxymethylene Ethers (OMEn) in Mixtures with a Diesel Surrogate
Previous Article in Journal
How Does Digitization Succeed in the Municipal Water Sector? The WaterExe4.0 Meta-Study Identifies Barriers as well as Success Factors, and Reveals Expectations for the Future
Previous Article in Special Issue
Large-Eddy Simulation of Laser-Ignited Direct Injection Gasoline Spray for Emission Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Rapid Compression Machine Study of 2-Phenylethanol Autoignition at Low-To-Intermediate Temperatures

Department of Mechanical Engineering, University of Connecticut, Storrs, CT 06269, USA
*
Author to whom correspondence should be addressed.
Energies 2021, 14(22), 7708; https://doi.org/10.3390/en14227708
Submission received: 19 October 2021 / Revised: 9 November 2021 / Accepted: 16 November 2021 / Published: 17 November 2021
(This article belongs to the Special Issue Ignition and Combustion Characteristics of Automotive Fuels)

Abstract

:
To meet the increasing anti-knocking quality demand of boosted spark-ignition engines, fuel additives are considered an effective approach to tailor fuel properties for satisfying the performance requirements. Thus, screening/developing bio-derived fuel additives that are best-suited for advanced spark-ignition engines has become a significant task. 2-Phenylethanol (2-PE) is an attractive candidate that features high research octane number, high octane sensitivity, low vapor pressure, and high energy density. Recognizing that the low temperature autoignition chemistry of 2-PE is not well understood and the need for fundamental experimental data at engine-relevant conditions, rapid compression machine (RCM) experiments are therefore conducted herein to measure ignition delay times (IDTs) of 2-PE in air over a wide range of conditions to fill this fundamental void. These newly acquired IDT data at low-to-intermediated temperatures, equivalence ratios of 0.35–1.5, and compressed pressures of 10–40 bar are then used to validate the 2-PE model developed by Shankar et al. (2017). It is found that this literature model greatly overpredicts the current RCM data. The comparison of experimental and simulated results also provides insights into 2-PE autoignition behaviors at varying conditions. Further chemical kinetic analyses demonstrate that the absence of the O2-addition pathway of β- R . radical in the 2-PE model of Shankar et al. (2017) could account for the model discrepancies observed at low-to-intermediated temperatures.

Graphical Abstract

1. Introduction

The anti-knocking quality of a fuel is essential to spark-ignition (SI) engines for achieving higher compression ratios and thus higher thermodynamic efficiency. Fuels with high octane rating are desired for advanced SI engines because they not only can meet the increasing requirement on fuel’s anti-knocking quality but also exhibit reduced CO2 emissions [1]. As conventional gasoline will still be the fuel for running most of the SI engines, fuel additives become a viable way to tailor fuel properties for meeting those anti-knocking performance requirements. Moreover, with escalating global concerns about energy security and environmental sustainability, whether the candidate fuel additives can be derived from domestic biomass should be taken into consideration. Therefore, it is imperative to identify top bio-derived fuel additives that are suitable for further development and have the potential to enable more efficient engine design and operation.
2-Phenylethanol (2-PE) is an alcoholic compound from the aromatics family with high research octane number (RON = 110) and high octane sensitivity (S = 20) [2], which well suits the role of an anti-knocking additive especially at high engine load conditions [3]. 2-PE not only exhibits other desirable features such as low vapor pressure [4] and high energy density [2], but also has the potential to be produced via a low-cost feedstock known as lignocellulose [5,6]. As such, 2-PE has drawn researchers’ attention in recent years, both as alternate fuel and anti-knocking additive. Ratcliff et al. [7] blended several bio-derived oxygenates, including 2-PE, with gasoline and evaluated the knocking resistance and particle emissions of each candidate. It was found in [7] that 2-PE exhibits excellent anti-knocking performance. Pan et al. [8] investigated the effects of 2-PE addition on engine performance and emissions using a downsized gasoline SI engine. Compared with a gasoline reference fuel, 2-PE addition resulted in decreased in-cylinder pressure and heat release rates [8]. The combustion phasing was also affected as the crank angle at 50% burn (CA50) location was retarded with 2-PE addition, which reduced engine torque output and combustion efficiency under optimal SI timing [8]. Furthermore, it was reported in [8] that 2-PE addition affected engine emissions, with increased unburned hydrocarbon/CO and reduced nitric oxide in the exhaust gas. Those literature engine studies revealed both the great potential of 2-PE as an anti-knocking additive and the challenge of using 2-PE, namely the altered combustion phasing.
A key to resolving the combustion phasing issue is our fundamental understanding of the 2-PE autoignition characteristics under engine-relevant conditions and the underlying chemical kinetics, which however is very limited. Only a few fundamental works investigated the chemical kinetics of 2-PE. Sakai et al. [9] studied thermal decomposition of 2-PE using quantum mechanical calculations at the CBS-QB3 level. Shanker et al. [10] experimentally explored 2-PE autoignition via cooperative fuel research engine and shock tube (ST) measurements. In addition, a detailed chemical kinetic model describing the high temperature chemistry of 2-PE autoignition was developed in [10] and validated against their experimental data. Otten et al. [11] determined the oxidation reaction products of 2-PE via synchrotron photoionization and deduced the key oxidation pathways of 2-PE. Apart from the above-mentioned fundamental studies, both fundamental experimental data and chemical kinetic model at low-to-intermediate temperatures (<1000 K) are still lacking.
In view of this fundamental void, rapid compression machine (RCM) experiments are conducted over a wide range of engine-like conditions. A set of ignition delay time (IDT) data of neat 2-PE in air is acquired at low-to-intermediate temperatures, which provides the much-needed autoignition data for model validation. In addition, the high temperature 2-PE model developed by Shanker et al. [10] is used to simulate the newly acquired RCM data. Based on the current RCM IDT data, the model performance is evaluated, and the major discrepancy sources are analyzed and identified to provide guidance to future 2-PE model refinement and development.

2. Experimental and Computational Methodologies

2.1. Experimental Specifications

The RCM facility used for IDT measurements of 2-PE/air mixtures consists of a heated reaction chamber, a piston assembly, a hydraulic cylinder, a pneumatic cylinder, and a driving air tank, as shown in Figure 1. As the detailed descriptions of the RCM used in this work have been documented in [12,13], a brief overview is given here. This RCM simulates a single compression stroke of an internal combustion engine, and the fuel–oxidizer mixture in the reaction chamber is rapidly compressed by the piston assembly to reach the desired temperature and pressure at the end of compression (EOC). The piston assembly is driven by compressed air and stopped by hydraulic forces. At the EOC, the piston assembly is held in place which creates a constant volume reactor enclosed by the piston and the reaction chamber walls. To ensure homogeneity of the reactant mixture during and after the compression stroke, the piston is machined with a crevice to suppress the roll-up vortex as well as to contain the cold boundary layer gas. As such, the present RCM provides an idealized engine-relevant environment for investigating fuel autoignition behaviors. In-cylinder pressure histories are considered as the primary diagnostic of the RCM and are measured by a piezo-electrical thermal-shock resistant pressure transducer (Kistler 6125C-C20) along with a charge amplifier (Kistler 5010B). The EOC temperature and pressure (TC and PC, respectively) can be adjusted independently by varying the overall compression ratio, initial pressure (P0), and initial temperature (T0). The fuel–oxidizer mixtures are prepared in a pre-vacuumed mixing tank which is connected to the RCM reaction chamber by a heated manifold.
The raw experimental pressure traces are processed using a Python package called UConnRCMPy [14] to determine the resulting TC, PC, and IDT. Specifically, TC can be computed from the pressure profile obtained in the compression stroke by applying the “adiabatic core hypothesis” [15] via T 0 T C γ γ 1 d T T = ln P c P 0 , where γ is the temperature-dependent specific heat ratio. This process is carried out through the Python interface of Cantera 2.3.0 [16], from which TC is generated following the procedure detailed in Section 2.2.1. Figure 2a presents the representative experimental pressure trace obtained by the current RCM, as well as the definition of the total ignition delay time (τ). Namely, τ is defined as the time from the EOC to the maximum of the first-order time derivative of pressure (dP/dt). At least four consecutive runs are taken for each experimental condition to ensure repeatability, as shown in Figure 2b. Among repeated experiments, the one with τ closest to the average is selected as the representative pressure trace. The associated parameters (such as P0, T0) and outputs (such as PC, TC, and τ) of each representative experiment are reported instead of reporting all repeated experiments. The normal scatter of τ for repeated experiments is less than 15% of the reported value (while the standard deviations of repeated experiments are typically less than 6% of the reported values) and are plotted as error bars in the Arrhenius plots in the subsequent sections. To determine the machine-specified effect (i.e., the heat transfer effect [15]) on the reactive experiment as well as to confirm that there is no exothermicity during the compression stroke, the corresponding non-reactive experiment is conducted by replacing O2 in the reactive mixture with N2. Due to the similar specific heat values of O2 and N2, the change in the specific heat of the mixture is thus deemed negligible. Figure 2a demonstrates that the reactive pressure trace prior to the hot ignition is well-captured by its non-reactive counterpart.
High-purity (99.8%) 2-phenylethanol from Sigma-Aldrich, as well as high-purity (>99.99%) O2 and N2 from Airgas, are used to prepare the test mixtures in a pre-vacuumed stainless-steel mixing tank at room temperature. The liquid fuel (2-PE) is first injected into the mixing tank using a glass syringe. Subsequently, O2 and N2 are filled into the mixing tank consecutively based on barometric measurements. Adequate time interval is necessary between the fillings of two consecutive gases to minimize the temperature change (~5 K) in the mixing tank caused by filling the preceding gas. For “air”, the molar ratio of N2/O2 is kept at 3.76. A magnetic stirrer is placed at the bottom of the mixing tank to enhance mixing of the constituents. After filling the reactants, the heaters and magnetic stirrer are switched on and 4 h are allowed to completely vaporize the liquid fuel as well as stabilize the preheat temperature and ensure homogeneity of the reactant mixture. The preheat temperature used in this study is 420 K to prevent 2-PE condensation during experiments. The test conditions of this study, including equivalence ratio (ϕ), PC, TC, and mixture composition, are summarized in Table 1. The temperature regimes covered by the current RCM experiments are based on the consideration of avoiding excessively long τ and the constraints of machine setting capability.

2.2. Computational Specifications

2.2.1. Ignition Delay Time Simulations and TC Determination

Two types of IDT simulations are performed in this study. The first type is ‘RCM simulation’ which accounts for machine-specific effect by modeling it as change in reactor volume [17]. After confirming there is no exothermicity during the compression stroke by comparing the reactive pressure trace with its non-reactive counterpart, a combined pressure trace is then constructed and used to deduce the reactor volume history. The combined pressure trace uses the pressure profile during the compression stroke (before the EOC) from the reactive experiment, while the pressure profile after the EOC is taken from the corresponding non-reactive run. This is to ensure that the combined pressure trace faithfully captures the reactor condition during and after the compression stroke. Cantera [16] is used to convert the combined pressure trace to a volume trace by assuming that the mixture undergoes an isentropic compression during the compression stroke and an isentropic expansion after the EOC. With the generated volume trace as an input, the IdealGasReactor in Cantera [16] is used to yield the simulated reactive pressure trace, from which the IDT of RCM simulation is obtained following the same procedure of determining experimental τ described in Section 2.1. Moreover, by disabling the reactions (i.e., assigning the reaction multiplier to zero in Cantera) during RCM simulation, the temperature at the EOC is recorded as TC. More details about the RCM simulation and TC determination can be found in [14,18].
The second type of simulation is ‘CONV simulation’ that uses the IdealGasReactor in Cantera [16] assuming constant reactor volume. The CONV-simulated τ value is defined as the time required for the simulated temperature to increase by 400 K over the initial temperature, which is found to be close to that determined using the time derivative of pressure. The CONV simulations in this study are mainly used to simulate shock tube IDTs as well as perform model-based analyses.

2.2.2. Sensitivity Analyses

Besides IDT simulations, brute force sensitivity analyses of IDTs to the reaction rate coefficient of each reaction are carried out using CONV simulation in Cantera [16]. By perturbing the pre-exponential factor of the target reaction by a factor of 2 (0.5), the corresponding τ +   ( τ ) is obtained with the perturbed rate coefficient of k + = 2 k   ( k = 0.5 k ) , where k is the unperturbed reaction rate coefficient. The sensitivity coefficient for τ with respect to k is then defined as S τ = ln ( τ + / τ ) / ln ( k + / k ) . Thus, the reaction with a negative (positive) value of sensitivity coefficient promotes (retards) the overall reactivity.

3. Results and Discussion

3.1. Comparison with Literature Data

Figure 3 first compares the experimental τ of stoichiometric mixtures at 10 bar and 20 bar obtained from the current RCM experiments (2-PE/air) and the ST experiments (2-PE/Ar/O2) of Shankar et al. [10]. It is seen that the τ trends of the two studies in different temperature regimes have different slopes on the Arrhenius plot, with the RCM data at lower temperatures exhibiting steeper slopes than those of Shankar et al. [10] at higher temperatures. These slope differences can be attributed to two possible reasons: (1) different oxidizer compositions (i.e., air in RCM and highly diluted Ar/O2 in ST) and (2) different controlling autoignition chemistry in different temperature regimes. In addition, the pressure effects demonstrated by the two studies show different trends. While the highly-diluted ST data of Shankar et al. [10] display insignificant pressure effect on τ at high temperatures, a clear pressure effect is observed in the current RCM study as τ decreases noticeably with increasing pressure at low temperatures.
Figure 3 also shows that the CONV-simulated τ using the 2-PE model developed by Shankar et al. [10] agrees well with the ST data at high temperatures (>1000 K). However, comparing the RCM-simulated τ with the experimental results obtained in this work at lower temperatures (below 1000 K), the model of [10] significantly overpredicts τ up to a factor of 5. This suggests the change in controlling chemistry of 2-PE autoignition at different temperature regimes.

3.2. Reaction Pathway Analyses

As the 2-PE model of [10] displays significant overprediction against the RCM data at low-to-intermediated temperatures, the overall reaction pathways of this literature model are analyzed to identify the possible reasons causing such discrepancies. Figure 4 depicts the key reaction pathways/species included in the original model of [10] (marked in black). Note that the model of [10] did not include the O2-addition pathway on the β site. Consequently, the β- R . radicals mostly go through the styrene pathways, as shown in Figure 4. As the O2-addition pathways are generally insignificant at high temperatures, including or excluding those O2-addition pathways in the 2-PE model is not expected to affect the model performance at high temperatures. At low-to-intermediate temperatures, however, the O2-addition pathways become critical in the autoignition chemistry of hydrocarbon fuels due to their subsequent chain branching channels. Thus, the inclusion of the β site O2-addition pathway of 2-PE has the potential to improve the performance of the original model of [10] at low temperatures.
In view of the above, the original model in [10] is modified by adding four additional reactions: O2-addition on β- R . radicals, intramolecular H-abstraction (β- RO O . β- Q . OOH ), and two unimolecular decomposition reactions of β- Q . OOH . In Figure 4, the newly added reactions and species are marked with green arrows and shadows, respectively. These newly added pathways were proposed by Otten et al. [11] based on the experimentally identified reaction products of 2-PE oxidation. Due to the lack of theoretical calculations of those reaction rate parameters, they are determined herein via analogies to similar reactions of n-propylbenzene [19] and iso-octane [20]. Further, the thermodynamic properties of the two newly added species are acquired from the RMG database [21], which are estimated by using the group additivity method. This modified 2-PE model will be used to simulate the RCM/ST data and its performance will be compared with that of the original model of [10] in the following to illustrate the impact of the β site O2-addition pathway on the autoignition of 2-PE, especially at low-to-intermediate temperatures. The detailed information of the newly added species and reactions are provided in the Supplementary Material.

3.3. RCM Ignition Delay Times of 2-PE/Air Mixtures

Figure 5 compares the experimental τ obtained in this study at varying TC, PC, and ϕ. It is seen that τ monotonically decreases with increasing TC at all test conditions. Meanwhile, no negative temperature coefficient response is observed in the experiments and 2-PE only exhibits single-stage ignition behavior at all test conditions. The pressure effect on 2-PE autoignition is demonstrated in Figure 5 by varying PC at a given equivalence ratio, showing that increasing pressure monotonically reduces τ at all equivalence ratios investigated. For the stoichiometric condition shown in Figure 5c, although increasing PC from 10 bar to 20 bar leads to a large τ reduction, increasing PC from 30 bar to 40 bar only exhibits a very small τ reduction. At the fuel-lean conditions, while the ϕ = 0.5 case (Figure 5b) shows a somewhat similar pressure effect as the stoichiometric case, the extent of the pressure effect for the ϕ = 0.35 case is seen to remain unchanged with increasing PC. For the two pressures tested at the fuel-rich condition (ϕ = 1.5), Figure 5d illustrates that the τ variation with TC between PC = 20 bar and 30 bar changes minimally over the temperature range investigated here.
In addition to the pressure effect, the equivalence ratio effect is evaluated by varying ϕ at a given PC. As shown in Table 1, the mole fraction of O2 in the reactive mixture remains similar when varying ϕ from 0.35 to 1.5. Therefore, the present equivalence ratio effect represents the effect of fuel loading. In Figure 6, the experimental results show that increasing ϕ monotonically reduces τ at all PC investigated. It is also seen from Figure 6 that the extent of the equivalence ratio effect seems to be insensitive to the TC variation for a given PC.
Using RCM simulations, the performance of both the original model of [10] and the modified model is assessed by comparing the simulated τ with RCM data. It is seen from Figure 5 and Figure 6 that the original model (dashed lines) significantly overpredicts τ at all conditions typically in a factor of ~5, while the modified model (solid lines) agrees better with experimental results than the original model does at all tested conditions, with over-prediction of τ typically in a factor of ~2. This suggests that adding the β site O2-addition pathway alone in the 2-PE model of [10] greatly improves the model performance at low-to-intermediate temperatures as expected. Furthermore, both the original and modified models well capture the τ dependences on pressure and ϕ. Furthermore, note that both the original and modified models show gentler slopes on the Arrhenius plots than the experimental trends, especially at fuel-lean conditions as shown in Figure 5a,b. This is an important aspect to be improved in the future modeling effort of 2-PE autoignition by including all other pathways that are essential under low-temperature combustion conditions. In addition, although the modified model demonstrates better agreement with the RCM data compared to the original model, the discrepancies between the modified model and experiments are still large (a factor of 2), and hence further experimental and chemical kinetic modeling investigations are warranted.

3.4. Chemical Kinetic Analyses

To investigate the controlling chemistry of the original and modified models, brute force sensitivity analyses are further conducted for the stoichiometric 2-PE/air mixture to identify and compare the key reactions at two different initial temperatures of 850 K and 1500 K. As shown in Figure 7 for an initial pressure of 20 bar, the top ten promoting and retarding reactions in each model are listed. The original model of [10] is examined first. At 850 K, Figure 7a shows that the fuel H-abstraction reaction via O . H radical forming an α- R . radical exhibits a promoting effect while its β site counterpart displays an inhibiting effect. This is because the β- R . radicals in the original model mostly produce styrene, which is a “stable” species and will accumulate during the induction period. As the major destruction of styrene only occurs near hot ignition, the styrene produced during the induction period has no contribution to the growth of the radical pool, thereby leading to an inhibiting effect. By contrast, the α site pathway can route to the formation of C 6 H 5 C . H 2 , which can facilitate its subsequent chain-branching pathways in the toluene sub-model. Thus, the α site H-abstraction reaction via O . H radical in the original model exhibits a promoting effect.
In the modified model, however, the trend is reversed as seen in Figure 7a. The fuel H-abstraction reaction via O . H radical that forms an α- R . radical exhibits an inhibiting effect while its β site counterpart displays a promoting effect. This is because the newly added β site O2-addition pathway not only can bypass the styrene formation but also can link to the C 6 H 5 C . H 2 chemistry, as shown in Figure 4. The inclusion of the O2-addition pathway of β- R . radical shifts the overall effect of the β- R . radical formation reaction to promoting. As a result of the competition over fuel molecules, the overall effect of the α- R . radical formation reaction shifts towards inhibiting. In addition, the two newly added reactions mark in blue in Figure 7a demonstrate their high impact on 2-PE autoignition at 850 K. The unimolecular decomposition reaction of C6H5etohQOOH-B C6H5CHO + CH2O + OH exhibits a strong promotion effect due to its ability to produces O . H . The intramolecular H-abstraction reaction of C6H5etohO2R-B C6H5etohQOOH-B also displays a promoting effect. Due to the high activation energy of this intramolecular H-abstraction reaction, it is the controlling reaction in this channel.
At 1500 K, it is seen from Figure 7b that both the original and modified models show identical sensitivity coefficient trends, indicating the diminishing role of the newly added O2-addition pathway of β- R . radical at high temperatures. This is because the rate coefficient of the O2-addition reaction reduces as temperature increases while those of the competing high temperature pathways (high activation energy) increase rapidly with increasing temperature. As such, the branching percent of 2-PE that enters the O2-addition pathway on both α and β sites becomes negligible at 1500 K, thereby minimizing the impact of the O2-addition pathways on the overall reactivity. This is also reflected in the sensitivity analysis shown in Figure 7b, as the O2-addition reactions on the α and β sites are not in the top-10 list.
Figure 8a further compares the simulated results using the original and modified models with the ST data of [10]. It is found that the CONV-simulated τ values using the two models are identical and agree well with the experimental IDT data, showing that the inclusion of the O2-addition pathway of β- R . radical does not affect the model performance at high temperatures. Figure 8b presents the CONV-simulated τ using the original and modified models across different temperature regimes. The modified model is seen to be faster than the original model on the low temperature side. As temperature increases, the difference between the two models gradually reduces. Beyond 1100 K, the two models essentially yield the same results, which is consistent with the findings from the sensitivity analyses illustrated in Figure 7b. As the model-based analyses reasonably explain the discrepancies when validating the original and modified models against the RCM/ST IDT data, it is concluded that the inclusion of the β site O2-addition pathway is needed to improve the comprehensiveness of the model of [10] for predicting the autoignition of 2-PE under low-temperature combustion conditions. It must be pointed out that although the inclusion of the β site O2-addition pathway greatly improves the model performance at low-to-intermediate temperatures, a systematic construction of detailed low temperature reaction sets is required for rigorous model refinement and development. Obviously, low-temperature chemistry of 2-PE merits further experimental and computational investigations.

4. Conclusions

The autoignition experiments of 2-phenylethanol (2-PE) in air are performed in this study using a rapid compression machine (RCM) at varying compressed temperatures (TC), compressed pressures (PC), and equivalence ratios (ϕ). Experimental results show that the total ignition delay time (τ) of 2-PE monotonically decreases with increasing TC, PC, and ϕ. Moreover, only single-stage ignition behavior is exhibited for the conditions investigated and there is no negative temperature coefficient response observed in the experiments. While the literature ST study [10] showed an insignificant pressure effect on τ at high temperatures, the present RCM results demonstrate stronger influences of pressure on τ at lower temperatures. Comparison of the RCM data acquired in this work and the literature shock tube (ST) data in [10] with the simulated results using the 2-PE model in [10] shows that the model in [10] well predicts the high-temperature ST data of [10], but greatly overpredicts the current RCM results at low-to-intermediate temperatures. This suggests the change in controlling chemistry of 2-PE autoignition at different temperature regimes.
To examine the discrepancy sources of the 2-PE model in [10] at low-to-intermediate temperatures, a modified 2-PE model is compiled by adding the β site O2-addition pathway to the original model. While the modified model still overpredicts the current RCM data, it presents better agreement against experimental results than the original model does. This indicates that the inclusion of the β site O2-addition pathway can potentially account for the discrepancies of the original model. Additionally, the modified model still well predicts the ST data of [10] and demonstrates identical performance at high temperatures compared to the original model. Sensitivity analyses using the original and modified models illustrate that the inclusion of the β site O2-addition pathway shifts the overall effect of the β- R . radical formation reaction from inhibiting to promoting. As the importance of the O2-addition reactions greatly reduces with increasing temperature, the impact of the β site O2-addition pathway diminishes at high temperatures. Therefore, the discrepancies between the model in [10] and the current RCM data can be partially attributed to the absence of the O2-addition pathway of β- R . radical in the original 2-PE model.
This RCM study provides validation datasets for model development/refinement and insights into the 2-PE autoignition chemistry at low-to-intermediate temperatures. The potential issues in the 2-PE model in [10] are also investigated and identified to offer guidance for future development of a comprehensive 2-PE kinetic model that includes detailed low temperature reactions for 2-PE autoignition.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/en14227708/s1, (1) Information of the added species/reactions in the modified 2-PE model (Added species_reactions.txt). (2) RCM experimental pressure traces and deduced time-volume histories for RCM simulations.

Author Contributions

Conceptualization, C.-J.S. and R.F.; data acquisition, R.F.; simulation, R.F.; results analyses, R.F.; writing—original draft preparation, R.F.; writing—review and editing, C.-J.S.; supervision, C.-J.S.; project administration, C.-J.S.; funding acquisition, C.-J.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Lawrence Livermore National Laboratory via Standard Research Subcontract No. B630575.

Data Availability Statement

RCM experimental pressure traces and deduced time-volume histories for RCM simulations are available from the website of the Combustion Diagnostics Laboratory at the University of Connecticut https://combdiaglab.engr.uconn.edu/database/rcm-database (accessed on 15 November 2021).

Acknowledgments

The authors thank William J. Pitz (Materials Science Division, Lawrence Livermore National Laboratory, Livermore, CA 94551, USA) for supporting this work and providing valuable advice.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Nomenclature

2-PE2-phenylethanol
CA50crank angle at 50% burn
CONVconstant volume
EOCend of compression
IDTignition delay time
k reaction rate coefficient
PCpressure at the end of compression
P0initial pressure
RCMrapid compression machine
RONresearch octane number
Soctane sensitivity
SIspark ignition
STshock tube
S τ sensitivity coefficient for total ignition delay time
TCtemperature at the end of compression
T0initial temperature
ϕequivalence ratio
γspecific heat ratio
τtotal ignition delay time

References

  1. Wang, Z.; Liu, H.; Reitz, R.D. Knocking combustion in spark-ignition engines. Prog. Energy Combust. Sci. 2017, 61, 78–112. [Google Scholar] [CrossRef]
  2. Zhou, L.; Boot, M.D.; De Goey, L.P.H. The effect of the position of oxygen group to the aromatic ring to emission performance in a heavy-duty diesel engine. SAE Int. J. Fuels Lubr. 2012, 5, 1216–1239. [Google Scholar] [CrossRef]
  3. Szybist, J.P.; Splitter, D.A. Pressure and temperature effects on fuels with varying octane sensitivity at high load in SI engines. Combust. Flame 2017, 177, 49–66. [Google Scholar] [CrossRef]
  4. Emel’yanenko, V.N.; Dabrowska, A.; Verevkin, S.P.; Hertel, M.O.; Scheuren, H.; Sommer, K. Vapor pressures, enthalpies of vaporization, and limiting activity coefficients in water at 100 °C of 2-furanaldehyde, benzaldehyde, phenylethanal, and 2-phenylethanol. J. Chem. Eng. Data 2007, 52, 468–471. [Google Scholar] [CrossRef]
  5. McCormick, R.L.; Ratcliff, M.A.; Christensen, E.; Fouts, L.; Luecke, J.; Chupka, G.M.; Yanowitz, J.; Tian, M.; Boot, M. Properties of oxygenates found in upgraded biomass pyrolysis oil as components of spark and compression ignition engine fuels. Energy Fuels 2015, 29, 2453–2461. [Google Scholar] [CrossRef] [Green Version]
  6. Atsumi, S.; Hanai, T.; Liao, J.C. Non-fermentative pathways for synthesis of branched-chain higher alcohols as biofuels. Nature 2008, 451, 86–89. [Google Scholar] [CrossRef] [PubMed]
  7. Ratcliff, M.A.; Burton, J.; Sindler, P.; Christensen, E.; Fouts, L.; Chupka, G.M.; McCormick, R.L. Knock resistance and fine particle emissions for several biomass-derived oxygenates in a direct-injection spark-ignition engine. SAE Int. J. Fuels Lubr. 2016, 9, 59–70. [Google Scholar] [CrossRef] [Green Version]
  8. Pan, M.; Wei, H.; Feng, D.; Pan, J.; Huang, R.; Liao, J. Experimental study on combustion characteristics and emission performance of 2-phenylethanol addition in a downsized gasoline engine. Energy 2018, 163, 894–904. [Google Scholar] [CrossRef]
  9. Sakai, Y.; Ando, H.; Oguchi, T.; Murakami, Y. Thermal decomposition of 2-phenylethanol: A computational study on mechanism. Chem. Phys. Lett. 2013, 556, 29–34. [Google Scholar] [CrossRef]
  10. Shankar, V.S.B.; Al-Abbad, M.; El-Rachidi, M.; Mohamed, S.Y.; Singh, E.; Wang, Z.; Farooq, A.; Sarathy, S.M. Antiknock quality and ignition kinetics of 2-phenylethanol, a novel lignocellulosic octane booster. Proc. Combust. Inst. 2017, 36, 3515–3522. [Google Scholar] [CrossRef]
  11. Otten, A.; Wooten, M.; Medrano, A.; Fathi, Y.; Meloni, G. Investigation of oxidation reaction products of 2-phenylethanol using synchrotron photoionization. J. Phys. Chem. A 2018, 122, 6789–6798. [Google Scholar] [CrossRef] [PubMed]
  12. Mittal, G.; Sung, C.J. A rapid compression machine for chemical kinetics studies at elevated pressures and temperatures. Combust. Sci. Technol. 2007, 179, 497–530. [Google Scholar] [CrossRef]
  13. Das, A.K.; Sung, C.J.; Zhang, Y.; Mittal, G. Ignition delay study of moist hydrogen/oxidizer mixtures using a rapid compression machine. Int. J. Hydrog. Energy 2012, 37, 6901–6911. [Google Scholar] [CrossRef]
  14. Weber, B.W.; Sung, C.J. UConnRCMPy: Python-Based Data Analysis for Rapid Compression Machines. In Proceedings of the 15th Python in Science Conference, Austin, TX, USA, 10–16 July 2016; pp. 36–44. [Google Scholar]
  15. Lee, D.; Hochgreb, S. Rapid compression machines: Heat transfer and suppression of corner vortex. Combust. Flame 1998, 114, 531–545. [Google Scholar] [CrossRef]
  16. Goodwin, D.; Moffat, H.K.; Speth, R.L. Cantera: An Object-Oriented Software Toolkit for Chemical Kinetics, Thermodynamics, and Transport Processes; Version 2.3.0; Zenodo: Geneva, Switzerland, 2017; Available online: https://cantera.org/ (accessed on 18 January 2017).
  17. Sung, C.J.; Curran, H.J. Using rapid compression machines for chemical kinetics studies. Prog. Energy Combust. Sci. 2014, 44, 1–18. [Google Scholar] [CrossRef]
  18. Dames, E.E.; Rosen, A.S.; Weber, B.W.; Gao, C.W.; Sung, C.J.; Green, W.H. A detailed combined experimental and theoretical study on dimethyl ether/propane blended oxidation. Combust. Flame 2016, 168, 310–330. [Google Scholar] [CrossRef]
  19. Darcy, D.; Nakamura, H.; Tobin, C.J.; Mehl, M.; Metcalfe, W.K.; Pitz, W.J.; Westbrook, C.K.; Curran, H.J. A High-Pressure Rapid Compression Machine Study of n-Propylbenzene Ignition. Combust. Flame 2014, 161, 65–74. [Google Scholar] [CrossRef] [Green Version]
  20. Fang, R.; Kukkadapu, G.; Wang, M.; Wagnon, S.W.; Zhang, K.; Mehl, M.; Westbrook, C.K.; Pitz, W.J.; Sung, C.J. Fuel molecular structure effect on autoignition of highly branched iso-alkanes at low-to-intermediate temperatures: Iso-octane versus iso-dodecane. Combust. Flame 2020, 214, 152–166. [Google Scholar] [CrossRef]
  21. Gao, C.W.; Allen, J.W.; Green, W.H.; West, R.H. Reaction Mechanism Generator: Automatic construction of chemical kinetic mechanisms. Comput. Phys. Commun. 2016, 203, 212–225. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Rapid compression machine facility at the University of Connecticut.
Figure 1. Rapid compression machine facility at the University of Connecticut.
Energies 14 07708 g001
Figure 2. (a) Pressure trace of reactive experiment with stoichiometric 2-PE/air at PC = 30 bar and TC = 869 K, demonstrating the definition of the total IDT (τ). The non-reactive pressure trace by replacing O2 with N2 is also shown as a reference. (b) Representative and repeated experimental pressure traces with stoichiometric 2-PE/air at PC = 30 bar and TC = 869 K.
Figure 2. (a) Pressure trace of reactive experiment with stoichiometric 2-PE/air at PC = 30 bar and TC = 869 K, demonstrating the definition of the total IDT (τ). The non-reactive pressure trace by replacing O2 with N2 is also shown as a reference. (b) Representative and repeated experimental pressure traces with stoichiometric 2-PE/air at PC = 30 bar and TC = 869 K.
Energies 14 07708 g002
Figure 3. Experimental τ comparison of stoichiometric mixtures between the current RCM study (solid symbols) and the literature ST data of Shankar et al. 2017 [10] (hollow symbols). The simulated τ using the model of [10] are displayed as lines. Solid lines: CONV simulation. Dashed lines: RCM simulation.
Figure 3. Experimental τ comparison of stoichiometric mixtures between the current RCM study (solid symbols) and the literature ST data of Shankar et al. 2017 [10] (hollow symbols). The simulated τ using the model of [10] are displayed as lines. Solid lines: CONV simulation. Dashed lines: RCM simulation.
Energies 14 07708 g003
Figure 4. Reaction schematic of 2-PE autoignition proposed by Shankar et al. [10] (black arrow) and the newly added reactions (species) considered in the current investigation that are marked with green arrows (shadows).
Figure 4. Reaction schematic of 2-PE autoignition proposed by Shankar et al. [10] (black arrow) and the newly added reactions (species) considered in the current investigation that are marked with green arrows (shadows).
Energies 14 07708 g004
Figure 5. Experimental and simulated τ comparison for 2-PE/air mixtures at varying PC and (a) ϕ = 0.35, (b) ϕ = 0.5, (c) ϕ = 1.0, and (d) ϕ = 1.5. Experimental results are shown as symbols. Simulated results are displayed as lines. Solid lines: modified Shankar model. Dashed lines: original Shankar model [10].
Figure 5. Experimental and simulated τ comparison for 2-PE/air mixtures at varying PC and (a) ϕ = 0.35, (b) ϕ = 0.5, (c) ϕ = 1.0, and (d) ϕ = 1.5. Experimental results are shown as symbols. Simulated results are displayed as lines. Solid lines: modified Shankar model. Dashed lines: original Shankar model [10].
Energies 14 07708 g005
Figure 6. Experimental and simulated τ comparison for 2-PE/air mixtures at varying ϕ and (a) PC = 20 bar and (b) PC = 30 bar. Experimental results are shown as symbols. Simulated results are displayed as lines. Solid lines: modified Shankar model. Dashed lines: original Shankar model [10].
Figure 6. Experimental and simulated τ comparison for 2-PE/air mixtures at varying ϕ and (a) PC = 20 bar and (b) PC = 30 bar. Experimental results are shown as symbols. Simulated results are displayed as lines. Solid lines: modified Shankar model. Dashed lines: original Shankar model [10].
Energies 14 07708 g006
Figure 7. Brute force sensitivity coefficients of τ in the original and modified models for stoichiometric 2-PE in air based on CONV simulations at initial pressure of 20 bar and initial temperatures of (a) 850 K and (b) 1500 K. In panel (a), the newly added reactions in the modified model are marked in blue.
Figure 7. Brute force sensitivity coefficients of τ in the original and modified models for stoichiometric 2-PE in air based on CONV simulations at initial pressure of 20 bar and initial temperatures of (a) 850 K and (b) 1500 K. In panel (a), the newly added reactions in the modified model are marked in blue.
Energies 14 07708 g007
Figure 8. (a) Comparison of experimental and CONV-simulated τ for stoichiometric 2-PE/O2/Ar mixture at P5 = 10 and 20 bar per the test conditions of [10]. (b) Comparison of CONV-simulated τ between the original and modified models over a wide range of temperatures. Experimental results are shown as symbols. Simulated results are displayed as lines. Solid lines: modified Shankar model. Dashed lines: original Shankar model [10].
Figure 8. (a) Comparison of experimental and CONV-simulated τ for stoichiometric 2-PE/O2/Ar mixture at P5 = 10 and 20 bar per the test conditions of [10]. (b) Comparison of CONV-simulated τ between the original and modified models over a wide range of temperatures. Experimental results are shown as symbols. Simulated results are displayed as lines. Solid lines: modified Shankar model. Dashed lines: original Shankar model [10].
Energies 14 07708 g008
Table 1. Summary of test conditions for 2-PE/air mixtures in the current RCM experiments.
Table 1. Summary of test conditions for 2-PE/air mixtures in the current RCM experiments.
ϕPC (bar)TC (K)Molar Percent (%)
2-PEO2N2
0.3520, 30, 40886–9920.729920.855178.4150
0.520, 30, 40838–9701.039520.790078.1705
1.010, 20, 30, 40813–9372.057620.576177.3663
1.520, 30829–8933.055020.366676.5784
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fang, R.; Sung, C.-J. A Rapid Compression Machine Study of 2-Phenylethanol Autoignition at Low-To-Intermediate Temperatures. Energies 2021, 14, 7708. https://doi.org/10.3390/en14227708

AMA Style

Fang R, Sung C-J. A Rapid Compression Machine Study of 2-Phenylethanol Autoignition at Low-To-Intermediate Temperatures. Energies. 2021; 14(22):7708. https://doi.org/10.3390/en14227708

Chicago/Turabian Style

Fang, Ruozhou, and Chih-Jen Sung. 2021. "A Rapid Compression Machine Study of 2-Phenylethanol Autoignition at Low-To-Intermediate Temperatures" Energies 14, no. 22: 7708. https://doi.org/10.3390/en14227708

APA Style

Fang, R., & Sung, C. -J. (2021). A Rapid Compression Machine Study of 2-Phenylethanol Autoignition at Low-To-Intermediate Temperatures. Energies, 14(22), 7708. https://doi.org/10.3390/en14227708

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop