A Kinetic Modelling Study of Alcohols Operating Regimes in a HCCI Engine.
Modern society is facing multifaceted and complex energy-related challenges. Projection for the next 20-30 years foresee the combustion of fossil fuels to keep driving energy production, mostly due to higher energy demands for road, air and sea transport . While suitable alternatives to combustion exist for power generation, transportation requires high energy density sources, i.e. petroleum-derived liquid fuels (gasoline, diesel, kerosene, naphta). Moreover, environmental issues associated to combustion, pollution reduction targets  and other political and economic strategies endorse the use of biofuels (neat or in blends with commercial fuels) as the most promising near term alternative to fossil fuels for internal combustion engines.
First-generation fatty acid methyl-esters [3, 4] together with second-generation biofuels such as alcohols  and furans  have been one of the main topic of combustion chemistry research in the last decade. Particularly, large interest has been devoted to alcohols, due to the viable production pathways from biological matter [7, 8].
Starting from bioethanol, nowadays blended with gasoline up to 85% (E85 gasoline) and used in fuel flexible vehicles, an increasing interest in higher alcohols, with four or more carbon atoms has emerged because of favorable physical and thermodynamic properties . Energy density, heating value, viscosity, hygroscopicity and volatility of higher alcohols are in fact compatible with the requirements of existing distribution infrastructure (pipe lines, tanks etc.) and most importantly with modern engines. Therefore, the interest in [C.sub.4]-[C.sub.6] alcohols as fuel additives or as replacements for fossil fuels.
In the United States, the Octamix waiver  already allowed up to 16% butanol blends with gasoline as an equivalent to E10 gasoline. Many studies investigated the performances of biobutanol as a fuel or fuel additive. Rakopoulos et al.  as well as Siwale et al.  investigated blends of up to 24% v/v butanol in diesel, observing generally improved exhaust emission quality. Also a decrease in soot emissions was highlighted by Valentino et al. , due to the lower ignitability of butanol/diesel blends.
Pentanol is also an attractive second-generation bio-fuels and can be produced from renewable feedstock. Its higher energy density, higher heating value, higher viscosity, lower hygroscopicity and lower volatility  motivated the interest in pentanol isomers. Campos-Fernandez et al.  investigated power and fuel economy performance of diesel/pentanol blends (10-25% v/v of pentanol) in a direct injection compression ignition engine. Wei et al.  tested blends with up to 30% pentanol by volume. Li et al. experimentally investigated the emission performance, fuel economy and combustion characteristic of neat pentanol in a diesel engine . Overall these studies highlighted the potentials of the use of n-pentanol in diesel engines. Table 1 shows some important properties of several alcohols, gasoline and diesel fuels such as the lower heating value (LHV), research octane number (RON), motor octane number (MON) and cetane number (CN).
Correctly assessing the reactivity of a new fuel or a new gasoline or diesel formulation is largely a chemical kinetics problem. A better understanding of a specific chemical compound's effects on combustion performances (flame speed, auto-ignition etc.) and emission allow the design of a fuel or fuel blend for an existing technology, the tuning of an engine for an assigned fuel or the concerted development of fuels and engines . Typical surrogate mixtures such as primary reference fuels (PRFs) and toluene reference fuels (TRFs) have started to include different alcohols to investigate their possible impact on commercial fuels . From a chemical kinetic perspective, an accurate description of the combustion chemistry of every surrogate component (n-heptane, iso-octane, toluene, ethanol, butanol etc.) is imperative. The recent review of Sarathy et al.  discussed the state of the art understanding of alcohols combustion chemistry, synoptically pointing out many improvement margins and open questions. Following the inputs of Sarathy et al. , a joint research effort between the CRECK group at Politecnico di Milano (POLIMI) and the C3 group at National University of Ireland, Galway (NUIG) investigated the auto-ignition propensity of n-[C.sub.2]-[C.sub.5] alcohols in a twin-piston Rapid Compression Machine (RCM), with the aim of systematically revising alcohols low temperature chemistry.
Along with new fuel formulations to be used in conventional spark ignition (SI) and compression ignition (CI) engines, new engine technologies are also being investigated . Homogeneous Charge Compression Ignition (HCCI) engines have received great attention in the last few decades because of the possibility to obtain high power output with a cleaner combustion . In HCCI engines, fuel and air are fully premixed, as in conventional SI engines, and the charge is ignited through compression, relying on its auto-ignition characteristics, like in CI engines. Although the use of a premixed air-fuel charge ensures reduced NOx and soot emissions, it also causes High Heat Release (HRR) and Pressure Rise Rates (PRR) leading to ringing events , strongly limiting the power output. Charge dilution by means of lean equivalence ratio and Exhaust Gas Recirculation (EGR) allow preventing these phenomena, but excessive charge dilution leads to very low power output and unburned hydrocarbon emissions. Typical undesired phenomena such as knock, already limiting the efficiency of SI engine, are also an issue in HCCI engines, with the additional criticism derived from the lack of cycle control measures such as spark timing (SI engine) or fuel injection timing (CI engine).
Beside extensive experimental investigations of HCCI operating ranges, also including neat n-butanol or blends in the most recent literature [21,22, 23, 24, 25], some computational model providing a time- and cost-effective solution to analyze and optimize the engines, have been reported in the literature [26, 27, 28,29, 30, 31, 32]. In particular, stemming from the well-recognized fundamental role of chemical kinetics in HCCI combustion, and from the general lack of knowledge about the chemical phenomena ruling HCCI operability maps, Bissoli et al. (Energy & Fuels, 2017, Submitted) validated a multi-zone model to study the impact of different fuels and engine configurations (speed and boost) on the operability maps, from a chemical point of view, by means of advanced tools like Sensitivity Analysis and Rate of Production Analysis .
Starting from the description and validation of the updated POLIMI alcohols mechanism, this study investigates the chemical kinetics underlying the operability maps and efficiency of HCCI engines fueled with n-butanol, n-pentanol and a TRF/butanol gasoline surrogate .
The manuscript is organized as follows. Section 2 briefly discusses the experimental approach at NUIG also comparing alcohols ignition-delay time measurements with previous measurements from the literature. The kinetic mechanism of n-butanol and n-pentanol is then discussed, mostly focusing on the low temperature oxidation pathways. Comparison of experimental data with POLIMI model predictions for n-butanol, n-pentanol and PRF-TRFs/butanol mixtures is presented to prove reliability of the chemistry. Section 3 describes the main features of the multi-zone HCCI model and its validation method. Section 4 discusses operability and performance maps of HCCI engines charged with the same fuels for which ignition propensity is discussed in Section 2.
2. KINETIC MODEL OF N-BUTANOL AND N-PENTANOL LOW TEMPERATURE COMBUSTION
Modelling the autoignition in HCCI engines requires an accurate description of low and high temperature combustion kinetics of alcohols. Moving from the studies on propanol isomers by Frassoldati et al.  and of butanol isomers by Grana et al. , the high temperature kinetic mechanism has been recently revised and extended to pentanol isomers (n- and iso-pentanol) . A lumped low temperature mechanism to describe n-butanol oxidation was proposed by Pelucchi et al. . Starting from a systematic evaluation of the influence of different oxygenated functional groups on C-C and C-H bonds strength, a description of the main features and reaction pathways characterizing the reactivity of different oxygenated fuels (alcohols, aldehydes, ketones, methyl esters) was provided, with particular focus on intermediate and low temperature regimes . The same rules applied to n-butanol, were then extended to describe n-pentanol oxidation at low temperatures, according to the assumption that the effect of the functional group vanishes after the [beta] position [37, 38, 39], as schematically represented in Figure 1.
Stemming from the lack of a systematic investigation of alcohols auto-ignition chemistry under engine relevant conditions and from the scarcity of measurements at low temperatures and high pressures , ignition delay times were measured in the twin-piston rapid compression machine (RCM) at NUIG for n-butanol and n-pentanol/"air" stoichiometric mixtures at p=10-30 bar and T=700-925 K. The extended experimental targets for model validation guided a revision of the lumped low temperature mechanism, as discussed in Section 2.2.
New fuel formulations require the coupling of alcohols chemistry with standard PRFs and TRFs surrogates, representative of commercial fuels [33, 40]. The global POLIMI mechanism has been thoroughly validated for pure component such as n-heptane, iso-octane and toluene and their mixtures [41, 42, 43]. The mechanism uses a lumped description of the primary propagation reactions of larger species and primary intermediates [44, 45]. This approach, together with an extensive use of structural analogies and similarities within the different reaction classes, easily allows extension of the scheme to new species (e.g. n-hexanol), still maintaining a relatively low number of species (298 species, 11095 reactions). The kinetic mechanism, together with thermodynamic properties is provided in the Supporting Information.
2.1. Experimental Methodology
Ignition delay time measurements were carried out in the rapid compression machine described by Darcy et al.  at NUIG. n-butanol ([greater than or equal to]99%) and n-pentanol ([greater than or equal to]99%) were obtained from Sigma- Aldrich, while helium (99.9%), oxygen (99.5%), argon (99.9995%), nitrogen (99.95%), and carbon dioxide (99.5%) were supplied by BOC Ireland. Auto-ignition measurements for n-butanol and n-pentanol/"air" stoichiometric mixtures were conducted at temperatures of ~700-925 K, and pressures of 10-30 bar. Mixture compositions are reported in Table 2.
The experimental data thus obtained agree well with previous literature measurements in RCM and Shock Tubes (ST). Deviations between measurements at similar conditions have to be attributed to differences in the facility design and in diluent compositions ([N.sub.2]/Ar/C[O.sub.2]) in data from different authors. Deviations from ideal behavior caused by heat loss and complex fluid dynamics were taken into account for RCM simulations. Volume histories necessary for the correct simulation of the new data presented in this study are reported in the Supporting Information.
Figure 2a shows simulated pressure profiles for n-butanol and n-pentanol at Tc=830 K and 838 K, respectively. Experimental ignition delay times of the two alcohols at p=10 bar are compared in Figure 2b, highlighting the higher low temperature reactivity of n-pentanol, and a converging ignition propensity for T>~950 K.
2.2. Kinetic Mechanism
Prior to the systematic revision of n-butanol and n-pentanol kinetics, the thermodynamic properties of the fuels, fuel radicals (alkyl, peroxy, hydroperoxy alkyl etc.) and stable species (enols, ketohydroperoxides etc.) have been obtained based on the revised group contributions by Burke et al. . Properties of the lumped low temperature species have been obtained as a selectivity based average of the different possible isomers. It has to be noted that, according to the lumping procedure proposed by Ranzi et al. [44, 45] both the forward and the reverse rate constant in the low temperature pathways are explicitly assigned, nullifying the effect of updated thermochemistry of typical low temperature species. Figure 3 shows bond dissociation energies (BDE) for n-butanol C-H and C-C bonds, highlighting the weakening effect of the hydroxyl functional group on vicinal bonds. As summarized in Figure 1 such effect vanishes after the [beta]-positions.
Fuel consumption mostly occurs through H-abstraction to form fuel radicals. OH and H[O.sub.2] are the dominant abstracting radical over the whole temperature range of interest for HCCI combustion. The observed bond strength hierarchy directly impacts relative selectivity to the different H-abstraction channels. Rate constants for H-abstraction by OH, H[O.sub.2] and C[H.sub.3] on alcohols-specific positions have been adopted from Zhou et al. [48, 49, 50]. Alkane-like positions are treated according to the systematic approach of Ranzi et al. , as already reported in previous studies [34, 35, 36].
Figure 4 shows selectivity of H-abstractions by OH for both n-butanol and n-pentanol. As expected from the BDEs, the formation of [alpha] radical largely dominates, followed by the secondary alkane-like positions ([gamma] in n-butanol or [gamma] and [delta] in n-pentanol), and by the secondary [beta] position, whose C-H bond is ~1 kcal [mol.sup.-1] stronger than the alkane-like secondary positions. H-abstraction from the hydroxyl group only accounts for ~5 % of the overall selectivity.
Similarly to alkanes, the fuel radical can isomerize or decompose via [beta]-scission reactions . At lower temperatures alkyl radicals interact with oxygen forming peroxy radicals, activating the typical low temperature reaction pathways [5, 44, 53]. A schematic representation of the low temperature oxidation mechanism of alcohols is given in Figure 5, highlighting pathways that are particularly relevant or new pathways with respect to alkanes (thicker arrows).
The peculiarity of alcohols oxidation at low temperatures is that the [alpha]-hydroxyalkylradical (R-CH-OH) reacts with [O.sub.2] to rapidly form H[O.sub.2] and the parent aldehyde or ketone. This alcohol-specific propagation pathway strongly competes with the conventional alkane-like low temperature branching channel, largely justifying the relatively high RON and MON indices of alcohol fuels (ethanol, propanol and butanol isomers). As clear from the comparison between panel A and B of Figure 4, the importance of this reaction channel decreases with increasing chain length, making longer alcohols (pentanol, hexanol, etc.) gradually more similar to linear alkanes and therefore more suitable for diesel engines or new combustion technologies such as gasoline compression ignition (GCI) engines. As summarized by Kalghatgi , GCI engines could efficiently reduce N[O.sub.x] and soot by running CI engines using low octane gasoline (RON-70-85), spanning gasoline and diesel volatility ranges. n-pentanol and n-hexanol fall within such specifications both in terms of boiling temperatures (138 and 157 [degrees]C respectively) and ignition propensity.
The specific low temperature interactions of ethanol [alpha]-hydroxyalkyl radical with [O.sub.2] have been theoretically investigated by Zador et al.  and by da Silva et al. , clearly highlighting in the potential energy surface analysis the low lying pathway leading to the formation of H[O.sub.2] and acetaldehyde. Despite a more systematic theoretical evaluation of R-CH-OH + [O.sub.2] for a series of alcohols would be necessary in order to extrapolate meaningful rate rules, the high pressure limit rate constant of da Silva et al. , already adopted by Sarathy et al.  and Heufer et al. , has been corrected accounting for the extra ~4 kcal/mol reported by Zador  and co-workers at a higher level of theory. Figure 6 compares such rate constants.
The formation of a carbonyl compound and H[O.sub.2] is also considered in successive steps when internal isomerization reactions (R[O.sub.2][left and right arrow]QOOH) lead to the formation of a [alpha]-hydroxy-hydroperoxyalkyl radical, and its successive interactions with [O.sub.2] undergo similar pathways, as summarized in Figure 7.
The remaining fuel radicals produced by H-abstraction reactions follow the conventional low temperature pathways, whose rate constant are obtained from established alkanes rate rules  accounting for the different bond dissociation energies (i.e. C-H in [alpha] and [beta]) in isomerization and decomposition steps or from the mechanism of Sarathy  (e.g. tautomerization reactions, Waddington reactions etc.). The unconventional dehydration pathway of hydroperoxyalkyl (QOOH) radicals proposed by Welz et al.  is also included, and the rate constant is estimated taking into account the formation of a cyclic transition state and the calculated energy barrier (~13 kcal/mol). Further activity should be devoted to obtain a more accurate kinetic rate constant for this channel, whose steps are reported in Figure 8.
As clearly highlighted in the above discussion and schematically summarized in Figure 5, a large amount of aldehydes is produced in alcohols oxidation. In particular, the R-CH-OH+[O.sub.2] pathways directly produce the parent aldehydes R-(C=O)-H, whose high and low temperature oxidation have been included according to the work of Pelucchi et al. [57, 58].
The primary reaction in the lumped kinetic mechanism to describe the oxidation of n-butanol and n-pentanol is reported in Table 3.
Beyond the kinetic details of Table 3 it is possible to explain the higher reactivity of n-pentanol (Figure 2) through some rather simple observations. Beside the higher rate of consumption due to H-abstraction reactions by OH, the longer alkane-like carbon chain allows for enhanced isomerization steps in the low temperature branching pathways, as reported in Figure 9. In particular, Figure 9a shows forward (solid lines) and backward (dashed lines) R[O.sub.2][left and right arrow]QOOH isomerization lumped rate constants, whose ratio results in a ~20% higher equilibrium rate constant (Figure 9b) in n-pentanol oxidation compared to n-butanol.
This enhanced isomerization step translates into a more effective low temperature branching for n-pentanol, highlighted in terms of normalized ketohydroperoxides mole fraction in Figure 10. The ignition delay time also corresponds to the time at which complete consumption of ketohydroperoxides occurs.
2.3. Kinetic Mechanism Validation
The kinetic mechanism discussed in Section 2.2, has been validated over a wide range of experimental conditions: ignition delay times (T=670-1670 K, p=1-80bar, [PHI]=0.5-2.0), laminar flame speed (p=1 bar, T=343-473 K, [PHI]=0.7-1.7), speciation in jet stirred reactors (T=770-1100 K, p=1-10 bar, [PHI]=0.35-2.0). For sake of brevity only the ignition delay times are presented and discussed in the following section. Kinetic simulations have been performed using the OpenSMOKE++ framework of Cuoci et al. .
Figure 11 compares experimental ignition delay time [60, 61, 62, 63] for stoichiometric butanol/[O.sub.2]/Ar and butanol/air mixtures with kinetic model predictions, over a wide range of temperature and pressures. In a similar way, Figure 12 compares experimental [38, 64, 65] and calculated ignition delay time for n-pentanol. The kinetic mechanism is able to predict temperature and pressure dependence of ignition delay time for both fuels. Maximum deviations for n-butanol are within a factor of ~2 for the RCM data at 10 bar and the shock tube data at 8 bar (Figure 11, blue open squares and black open triangles). n-pentanol predictions systematically deviate from the 9 bar shock tube data  (Figure 12 green open diamond). Similar deviations were observed by Heufer et al.  and Sarathy et al.  (thin green line in Figure 12).
Sensitivity analyses have been carried out to investigate the chemistry responsible for the observed ignition trends. Figure 13 compares sensitivity coefficients at p=30 bar and T=650 K for n-butanol and n-pentanol. As expected from the fundamental similarities between the two alcohols and from the systematic development of the kinetic subsets to describe their combustion, the same hierarchy is observed within the important reaction pathways. The minor differences are due to more or less enhanced low temperature reactivity, resulting from different carbon chain length (i.e. longer alkane-like moiety). The sensitivity coefficients reported in Figure 13 have been scaled assuming a value of -1 for the most inhibiting reaction in both cases (OH+Fuel[left and right arrow]R[alpha]+[H.sub.2]O). At such temperature the most of the low temperature reactive flux produces ketohydroperoxides (KHYP) whose decomposition provides OH radicals to the H-abstraction reactions. The abstraction on the alkane-like moiety of the molecules largely increases reactivity, allowing the onset of the low temperature branching pathway. To a lesser extent this is also observed for the H-abstraction from the [beta] position. Clearly the H-abstraction from the [alpha]-site strongly decreases the reactivity as the only fate of such radical is to produce H[O.sub.2] and the corresponding aldehyde.
Additional interesting insights into the chemistry governing autoignition is given in Figure 14, showing the sensitivity coefficient of the most important fuel specific reactions as a function of temperature. For sake of clarity only the n-butanol case is reported and discussed, as one should expect analogous observations for n-pentanol.
Moving from the lowest temperature condition (T=650 K) already discussed, at T=850 K (blue bars of Figure 14) a rise in the relevance of H-abstractions by H[O.sub.2] is observed, despite this H-abstraction by OH still consumes the most of the fuel. H[O.sub.2] is mostly produced by H+[O.sub.2](+M[left and right arrow]H[O.sub.2](+M), and by HCO+[O.sub.2][left and right arrow]CO+H[O.sub.2] and by the same interaction of the [alpha] radical (R[alpha]) with [O.sub.2], discussed above (R-CH-OH+[O.sub.2][left and right arrow] H[O.sub.2]+ R-(C=O)-H). It has to be noted that such reaction does not show up in the sensitivity analysis as no alternative pathways exist for R[alpha]. In other words, this channel act as a "sink" of reactivity. H-abstractions by H[O.sub.2] directly compete with the termination reaction H[O.sub.2]+H[O.sub.2][left and right arrow][O.sub.2]+[H.sub.2][O.sub.2] strongly inhibiting auto-ignition. The competition between low temperature and intermediate temperature pathways clearly emerges at T=850 K. Observing the high pressure shock tube data of n-butanol in Figure 11 it is possible to observe a change of slope around this temperature. This phenomenon becomes even more pronounced in the case of n-pentanol, gradually approaching the typical negative temperature coefficient (NTC) behavior of n-alkanes. This is due to the increasing importance of alternative decomposition pathways of R[O.sub.2] and QOOH radicals and to the direct H-abstraction by [O.sub.2] on fuel radicals forming H[O.sub.2] and an aldehyde or an unsaturated alcohol ([O.sub.2]+R[left and right arrow]H[O.sub.2]+Aldehyde/Unsaturated Alcohol). Such alternative channels overcome the low temperature branching, simply propagating the radical chain reaction or, at the limit, inhibiting the reactivity producing less reactive radicals (e.g. H[O.sub.2]). These statements are confirmed by the increased sensitivity coefficient of the isomerization reaction R[O.sub.2][left and right arrow]QOOH and of the successive second addition to [O.sub.2] ([O.sub.2]+QOOH[left and right arrow]OOQOOH), playing a key role in defining the relative importance of the reactive flux giving branching or propagation.
At higher temperature (T=1200 K), the number of fuel specific reactions which are sensitive to ignition delay time determination largely decreases. In fact, the characteristic time of decomposition of fuel radicals (~[10.sup.9] [s.sup.-1]) does not allow for any interaction with molecular oxygen, resulting in dominant role of the [C.sup.0]-[C.sup.2] chemistry.
Recent experimental and kinetic modelling study focused on PRF and TRF mixtures blended with butanol [33, 40]. As an additional validation target of the proposed n-butanol and n-pentanol low temperature mechanism, comparison of the POLIMI mechanism with these experimental data are reported in Figure 15 and Figure 16.
3. MULTI-ZONE MODEL OF HCCI ENGINE
3.1. Engine Simulation Model
The multi-zone model of HCCI engine developed and validated by Bissoli et al. [32, 67] is briefly reiterated in this Section. Figure 18 shows the multi-zone configuration, conceived according to the " owiow-skiw" approach firstly introduced by Komninos et al. .
The simulation starts at the Intake Valve Closing (IVC) and ends at Exhaust Valve Opening (EVO), therefore it only describes the compression and expansion phases of a four-stroke engine. During each cycle the model evaluates system reactivity in each zone, together with heat and mass exchanged between adjacent zones. The thermal and composition stratification in the charge are accounted for by means of interactions between neighboring zones. A complete and adiabatic mixing between the exhaust gas trapped in the cylinder and the fresh EGR is assumed. EGR temperature and composition are updated at each cycle, thus affecting the initial conditions of the new cycle. In this work, the composition of the mixture at the beginning of the first cycle is assumed to be uniform and equal to the nominal one (i.e. without EGR).
The different assumptions underlying the model formulation are summarized herein:
1. In-cylinder mixture described as an ideal-gas,
2. Each zone is treated as an ideal reactor with uniform temperature, composition and time-variable volume,
3. Uniform pressure, except for the crevices (constant-volume, variable-mass).
Uniform pressure, except for the crevices (constant-volume, variable-mass).
Details on the heat transfer model, on laminar and turbulent contributions to heat and mass transfer have been discussed in the previous study of Bissoli et al. .
The multi-zone model was specifically conceived for simulations with detailed kinetic mechanisms and allows for conventional kinetic analyses such as Sensitivity Analysis and Rate of Production Analysis.
The validation of the model was previously reported for n-heptane, n-butanol, and methyl-esters  (+Ewergy & Fuels, 2017, Submitted). The multi-zone model was able to properly reproduce compression and expansion phases, pressure peak and the combustion phasing. Moreover, very good predictions of reactivity and emissions were obtained, based on CO, C[O.sub.2] and other intermediate species (aldehydes) profiles as a function of the compression ratio (CR).
The recent work of Bissoli et al. (Ewergy & Fuels, 2017, Submitted) focused on the experimental measurements obtained in a Ricardo E6 engine at Brunel University  for PRF80, PRF100 and ethanol. For the sake of comparison, the simulations performed in this study are carried out in the same engine configuration. Engine characteristics and specific model parameters are reported in Table 4. These parameters were obtained by performing simulations using 15 zones, with a crevice volume equal to 2% of the clearance volume. As already discussed , such number of zones constitutes a good compromise between computational efforts and the detailed description of the phenomena involved.
3.2. Operability and Performance Maps of HCCI Engines
Oakley and co-workers [69, 70] systematically explored the air-fuel ratio (AFR) versus EGR operating range of the Ricardo E6 single-cylinder engine at fixed speed and compression ratio (CR) for different fuels. Three zones limiting the stability of HCCI regime were identified: the first one at low loads, called partial burn, characterized by low engine efficiencies; the knocking region, observed at high loads, where high pressure rise rates occur, and the misfire zone, at high dilutions, where high cycle-to-cycle variations are registered. Figure 19 shows a typical operative map for a HCCI engine.
The features of these three different combustion regimes are highlighted in Figure 20. High loads (low [lambda]) and low EGR cause the ringing phenomena, characterized by a strong PRR. Partial burn is characterized by an incomplete combustion, with consequent reduction of combustion efficiency ([[eta].sub.comb]), defined as:
[[eta].sub.comb] = [[integral]HRR/[DELTA][H.sub.combustion]]
Misfire is observed at high EGR dilutions and represents an unstable condition, where the engine goes towards periodical oscillations. This regime is characterized by strong increases in cycle-to-cycle variability, quantified in terms of variation of the indicated mean pressure (CoV IMEP) defined as:
CoVIMEP = [[[sigma].sub.IMEP]/IMEP]
where the IMEP is calculated as:
IMEP = [[integral]pdV/[V.sub.disp]]
and [[sigma].sub.IMEP] is its standard deviation.
The operability maps presented in this study are generated by independently varying the EGR ratio (mass-basis) and [lambda], considered as the global in-cylinder air-to-fuel ratio (AFR) after mixing EGR with the fresh feed. For each map more than 150 EGR and [lambda] combinations are investigated, with 50 simulation cycles for each point to account for the EGR cycle-cycle variability. Different characteristics such as ignition timing, indicated mean effective pressure (IMEP) and exhaust emission are plotted inside the different maps averaging the last 10 simulation cycles. The limits of the stable combustion region are defined as the set of [lambda]-EGR satisfying the requirement of [[eta].sub.comb] [greater than or equal to] 80%, CoV IMEP [less than or equal to] 60% and PRR [less than or equal to] 6 bar/deg. Such threshold values are deduced from the experimental study of Oakley  and are comparable with those adopted in other studies [71, 72].
4. RESULTS AND DISCUSSION
Figure 21 compares the predicted HCCI operating regions for the investigated fuels. Ethanol and the TRF/butanol mixture  representative of a PR5801 gasoline (RON=95, MON=86.6), slightly extend the operability map toward the ringing region, due to the presence of highly anti-knocking components such as ethanol, iso-octane and toluene. n-pentanol instead shows a trend which is more similar to PRF100, with butanol lying in between. As previously observed (Ewergy & Fuels, 2017, Submitted), ringing phenomena are scarcely affected by the fuel type. n-butanol and n-pentanol show the highest flexibility to both EGR dilution and engine load, with the TRF/butanol surrogate showing a behavior in between ethanol and PRFs mixtures.
Figure 22 and Figure 23 exemplify the behavior of the different classes of fuels discussed in this paper, by comparing the performance of PRF80 with the TRF/butanol mixture. As expected, Figure 22 shows that engine load increases when air and EGR ratios decrease, since stoichiometric and undiluted conditions are approached. Both fuels show similar and quite narrow load ranges, with similar peak values of ~2.5 bar. TRF mixture shows a wider stability, allowing to reach stable combustion at lower loads (~1 bar).
A similar behavior is observed for the indicated thermal efficiency in Figure 23. PRF80 and TRF/butanol mixtures shows very similar trends, with efficiency increasing with fuel amount in the intake charge. On the contrary, air dilution and large EGR ratios reduce it. Both fuels show that thermal efficiency ranges from ~30% to ~36%, with TRF/butanol mixture showing a slightly large lower limit of efficiency (~28%), related with the larger operating region close to partial burn conditions.
The reason behind the large stability region of n-butanol and n-pentanol have to be referred to their ignition propensity. Beside the kinetic features of long chain alcohols, already discussed in Section 2, Figure 24 shows the effect of engine load and EGR on ethanol and n-pentanol constant volume batch reactor ignition delay times.
Similar responses for ethanol and n-pentanol are observed for varying EGR dilution (right panels) at low temperatures (T<900 K), justifying the very similar stability limits toward the misfire region. n-pentanol is much more sensitive to load variation (left panels), in particular at low temperatures, resulting in wider limits at low loads (high [lambda], lower fuel concentration). Sensitivity analysis from the multi-zone model allows to better understand such trends. Figure 25 compares the sensitivity coefficient of ethanol, n-butanol and n-pentanol at [lambda]~5 and EGR=20%. As already discussed by Bissoli et al. (Ewergy & Fuels, 2017, Submitted) ethanol chemistry is dominated by hydrogen peroxide ([H.sub.2][O.sub.2]) and hydroperoxy radical chemistry (H[O.sub.2]). The lack of a significant alkane-like moiety typical of longer alcohols, prevents any low temperature reactivity, resulting in extremely high yields of H[O.sub.2], mostly formed by the interaction of the [alpha] radical (R[alpha], C[H.sub.3]-CH-OH) with [O.sub.2], producing acetaldehyde. As reported in Figure 25, the typical low temperature reactions emerge as strongly sensitive for n-butanol and n-pentanol. The availability of such pathways, results in higher yields of OH, sustaining the reactivity (shorter ignition delay times, Figure 26) and enlarging heavier alcohols operability maps. These statements are confirmed by the sensitivity diagram. Moreover, comparable sensitivity coefficients are observed for H[O.sub.2]+H[O.sub.2]=[O.sub.2]+[H.sub.2][O.sub.2] and for [H.sub.2][O.sub.2] (+M) = 2OH(+M) for the different fuels, while butanol and pentanol are more sensitive to H-abstraction reactions by OH, whose concentration is sustained by the low temperature branching.
The addition of butanol to the TRF mixture investigated by Agbro et al.  is responsible for the extension of the stable region toward lower loads. Ignition delay times at [lambda]~5 and EGR=20% have been calculated for such mixture and are compared with those of ethanol, n-butanol and n-pentanol in Figure 26.
The TRF/butanol mixture is the slowest to ignite for T>900 K, while at low temperatures it behaves very similarly to n-butanol. These observations justify both the slightly narrower ringing region overlapping ethanol, and the lower loads at which the HCCI engine can operate with respect to PRF80 and PRF 100. To better investigate the acting chemistry, and the kinetic influence of n-butanol, sensitivity analyses have been carried out at both T=850 K and T=1150 K (Figure 27).
At the lowest temperatures the reactivity is dominated by H-abstraction by OH from n-butanol, producing R[alpha]. n-heptane and iso-octane low temperature reactions rule the ignition propensity, together with H-abstraction from the alkane-like moiety of n-butanol. At higher temperatures (T=1150 K) the reactivity is dominated by reactions belonging to the core [C.sup.0]-[C.sup.2] mechanism, together with the unimolecular intiation of iso-octane, the most abundant component in the fuel mixture, responsible for the overall radical chain initiation. Such intiation involves the formation of alkyl radicals such as iso-propyl (i[C.sub.3][H.sub.7]), iso-butyl, tert-butyl and neo-pentyl radical, whose successive decomposition reactions produce methyl radical and unsaturated species ([C.sub.3][H.sub.6], iso-[C.sub.4][H.sub.8]), sequentially leading to the formation of resonance-stabilized radicals (allyl, methylallyl). The abundance of the primary C-H site available for H-abstraction also justifies the lower ignitability of iso-octane at high temperatures.
Figure 28 shows predicted pressure traces at different dilutions and EGR ratios for n-butanol, n-pentanol and TRF/n-butanol mixture. These results are not filtered and highlight how the present model, thanks to the description of the interactions among neighboring zones is able to reproduce also the effects of charge stratification on the HCCI combustion progression, similarly to what observed in the experiments Bissoli et al. . Pressure traces at similar conditions (EGR=20, [lambda]~3) for the different fuels are also compared in the bottom right panel. Higher EGR (EGR=50) and lower loads ([lambda]~5) result in a smoother transition to the hot ignition in the case of alcohols, as highlighted by the different first derivatives at the inflection point (CAD ~-4.1 deg). The following analysis focuses on the two limiting cases: the high EGR dilution and low loads (EGR50, [lambda]=5) and the low EGR dilution and high loads (EGR20, [lambda]=3). Sensitivity analyses have been performed at the conditions of temperature, pressure and mixture composition marked by the black crosses in Figure 28, to investigate the relevant kinetics. Figure 29 and Figure 30 show the results.
Starting from the lower temperature, lower pressure cases (CAD=-30, T~900 K, p~9 bar) both mixtures are largely dominated by H[O.sub.2] chemistry, similarly to what observed for the cases of Figure 24. Namely, the reactivity is strongly inhibited by the termination step H[O.sub.2]+H[O.sub.2]=[O.sub.2]+[H.sub.2][O.sub.2] and enhanced by hydrogen peroxide decomposition ([H.sub.2][O.sub.2] (+M)=2OH(+M)). The reactivity is also strongly enhanced by H-abstractions by H[O.sub.2] from the weak C-H bond in [alpha]. In fact, this channel provides a large amount of [H.sub.2][O.sub.2], whose decomposition rules the reactivity. The importance of the correct definition of the relative weight of the channels involving QOOH (i.e. second addition to [O.sub.2], decomposition to H[O.sub.2] and unsaturated alcohols, backward isomerization to R[O.sub.2] etc.) are also highlighted in Figure 29. The intermediate temperature conditions move the focal point of the kinetics toward species preceding ketohydroperoxides formation (QOOH, OOQOOH). Once again, the relative selectivity of H-abstraction by OH from the different carbon site plays a major role.
Figure 30 highlights the important reactions at the inflection point, characterizing the transition to hot ignition. H-abstraction by H[O.sub.2] at the [alpha] position becomes the most important reaction in enhancing reactivity. Differently from what previously discussed for OH, where the H-abstraction from the same site generally reduces the reactivity, the direct formation of [H.sub.2][O.sub.2] associated with H-abstractions by H[O.sub.2] has a positive effect of reactivity. This is due to the extreme relevance of hydrogen peroxide decomposition: the formation of R[alpha], further propagates the reactivity producing additional H[O.sub.2] through the same interaction discussed above (R-CH-OH+[O.sub.2][left and right arrow] H[O.sub.2]+ R-(C=O)-H).
Concerning the higher first derivative at the inflection point observed for the EGR20 [lambda]=3 case, this is mostly due to a higher concentration of fuel radical in the inlet charge, speeding up fuel decomposition thus providing higher yields of radical species, including R[alpha].
The ignition of the TRF/butanol mixture is delayed of ~10 CAD. From a kinetic perspective, as previously discussed at the same EGR and [lambda] conditions in Figure 27, this has to be referred to the lower reactivity of iso-octane and toluene.
Figure 31 compares predicted formaldehyde (C[H.sub.2]O) and CO emissions for PRF80 and the TRF/butanol surrogates. C[H.sub.2]O and CO are typical pollutants largely released when incomplete combustion phenomena occur. PRF80 produces higher emissions of formaldehyde, found to be a factor of ~2-5 higher considering the same load and EGR conditions. Also CO emissions are higher of a factor of ~2-3 in the case of PRF80. In general formaldehyde and CO emissions increase for decreasing load and increasing EGR dilution. This phenomena is largely motivated by the decreasing thermal efficiency (Figure 23), typically observed when moving toward leaner and more diluted conditions.
At low temperatures, C[H.sub.2]O directly derives from ketohydroperoxides decomposition. Moving toward higher temperatures, C[H.sub.2]O production is ruled by interactions of H[O.sub.2] and C[H.sub.3] radicals producing methoxy radical (C[H.sub.3]O) through the reaction H[O.sub.2]+C[H.sub.3][left and right arrow]C[H.sub.3]O+OH. The further dehydrogenation of C[H.sub.3]O largely explains C[H.sub.2]O production. The larger H[O.sub.2] yields previously mentioned in the case of n-butanol addition, result in an enhanced H[O.sub.2]+H[O.sub.2]=[O.sub.2]+[H.sub.2][O.sub.2] channel, partly limiting formaldehyde formation. Further oxidative steps of formaldehyde involve the H-abstraction producing formyl radical (HCO), rapidly decomposed to produce CO.
A systematic experimental analysis of linear [C.sup.2]-[C.sub.5] alcohols ignition delay time at high pressure (10-30 bar) and low temperatures 700-925 K, allowed the revision of the POLIMI mechanism for alcohols oxidation. Firstly, this study presented and discussed the new experimental measurements for n-butanol and n-pentanol. The development of the kinetic mechanism is then discussed focusing on alcohol-specific reaction pathways in the low and intermediate temperature oxidation regime. Model comparison with a wide set of ignition delay time measurements from this and from previous studies proves reliability of the proposed kinetic model. Following the recent trends of fuel formulation, requiring the blending of biofuels such as alcohols to conventional hydrocarbon fuels, the kinetic model is also compared to PRF and TRF/n-butanol mixtures, further confirming the validity of the alcohol sub-mechanism. A broad kinetic discussion serves the goal of characterizing the ignition propensity of such fuels, highlighting key reaction steps.
The HCCI multi-zone model of Bissoli et al. , widely validated and recently used to highlight the promising performances of ethanol (Ewergy & Fuels, 2017, submitted), was used to investigate the operability maps of n-butanol, n-pentanol and the same TRF/n-butanol mixture, focusing on the kinetic features governing the performances of such fuels in HCCI engine. n-butanol and n-pentanol extend the operability maps of ethanol even further, mostly in terms of allowing lower loads, thus reducing the partial-burn region. Similarly, the addition of n-butanol to a TRF mixture representative of a PR5801 gasoline, significantly impacts the ignition chemistry, slightly extending the operating limits of PRF100 and PRF80 previously analyzed by Bissoli et al. (Ewergy & Fuels, 2017, submitted).
From a chemical kinetics perspective, an alcohol molecule can be divided into two moieties. The first one, alcohols specific, is strongly influenced by the presence of the hydroxyl moiety (R-OH), influencing C-H and C-C bond strengths and largely preventing alkane-like low temperature branching pathways as discussed in Section 2. The oxidation of the remaining moiety instead proceeds according to conventional pathways. Therefore, the longer the alkane-like portion the closer the [C.sub.n] alcohol ignition propensity to that of the [C.sub.n] alkane. As discussed by Pelucchi et al.  and by Heufer et al.  these chemical and kinetic features translate into a higher reactivity at higher temperatures and a lower reactivity at lower temperatures, when comparing alcohols with the homologous alkanes.
[1.] Exxonmobil, The Outlook for Energy: A View to 2040, 2015.
[2.] C2ES, Outcomes of the UN climate change conference in Paris. http://www.c2es.org/docUploads/cop-21-paris-summary-02-2016-final.pdf.
[3.] Coniglio L., Bennadji H., Glaude P.A., Herbinet O., Billaud F., Combustion chemical kinetics of biodiesel and related compounds (methyl and ethyl esters): experiments and modeling-advances and future refinements, Progress in Energy and Combustion Science 39 (2013) 340-382.
[4.] Shahabuddin M., Liaquat A., Masjuki H., Kalam M., Mofijur M., Ignition delay, combustion and emission characteristics of diesel engine fueled with biodiesel, Renewable and Sustainable Energy Reviews 21 (2013) 623-632.
[5.] Sarathy S.M., O[beta]wald P., Hansen N., Kohse-Hoinghaus K., Alcohol combustion chemistry, Progress in Energy and Combustion Science 44 (2014) 40-102.
[6.] Sudholt A., Cai L., Heyne J., Haas F.M., Pitsch H., Dryer F.L., Ignition characteristics of a bio-derived class of saturated and unsaturated furans for engine applications, Proceedings of the Combustion Institute 35 (2015) 2957-2965.
[7.] Atsumi S., Hanai T., Liao J.C., Non-fermentative pathways for synthesis of branched-chain higher alcohols as biofuels, Nature 451 (2008) 86-89.
[8.] Peralta-Yahya P.P., Keasling J.D., Advanced biofuel production in microbes, Biotechnology journal 5 (2010) 147-162.
[9.] EPA, Federal Register Environmental Protection Agency Federal Register, 2012.
[10.] Rakopoulos D., Rakopoulos C., Giakoumis E., Dimaratos A., Kyritsis D., Effects of butanol-diesel fuel blends on the performance and emissions of a high-speed DI diesel engine, Energy Conversion and Management 51 (2010) 1989-1997.
[11.] Siwale L., Kristof L., Adam T., Bereczky A., Mbarawa M., Penninger A., Kolesnikov A., Combustion and emission characteristics of n-butanol/diesel fuel blend in a turbo-charged compression ignition engine, Fuel 107 (2013) 409-418.
[12.] Valentino G., Corcione F.E., Iannuzzi S.E., Serra S., Experimental study on performance and emissions of a high speed diesel engine fuelled with n-butanol diesel blends under premixed low temperature combustion, Fuel 92 (2012) 295-307.
[13.] Campos-Fernandez J., Arnal J.M., Gomez J., Lacalle N., Dorado M.P., Performance tests of a diesel engine fueled with pentanol/diesel fuel blends, Fuel 107 (2013) 866-872.
[14.] Wei L., Cheung C., Huang Z., Effect of n-pentanol addition on the combustion, performance and emission characteristics of a direct-injection diesel engine, Energy 70 (2014) 172-180.
[15.] Li L., Wang J., Wang Z., Xiao J., Combustion and emission characteristics of diesel engine fueled with diesel/biodiesel/pentanol fuel blends, Fuel 156 (2015) 211-218.
[16.] Kalghatgi G.T., The outlook for fuels for internal combustion engines, International Journal of Engine Research, (2014) 1468087414526189.
[17.] Bergthorson J.M., Thomson M.J., A review of the combustion and emissions properties of advanced transportation biofuels and their impact on existing and future engines, Renewable and Sustainable Energy Reviews 42 (2015) 1393-1417.
[18.] Morgan N., Smallbone A., Bhave A., Kraft M., Cracknell R., Kalghatgi G., Mapping surrogate gasoline compositions into RON/MON space, Combustion and Flame 157 (2010) 1122-1131.
[19.] Yao M., Zheng Z., Liu H., Progress and recent trends in homogeneous charge compression ignition (HCCI) engines, Progress in Energy and Combustion Science 35 (2009) 398-437.
[20.] Epping, K., Aceves, S., Bechtold, R., and Dec, J., "The Potential of HCCI Combustion for High Efficiency and Low Emissions," SAE Technical Paper 2002-01-1923, 2002, doi:10.4271/2002-01-1923.
[21.] Komninos N., Rakopoulos C., Modeling HCCI combustion of biofuels: A review, Renewable and Sustainable Energy Reviews 16 (2012) 1588-1610.
[22.] Liu M.-B., He B.-Q., Zhao H., Effect of air dilution and effective compression ratio on the combustion characteristics of a HCCI (homogeneous charge compression ignition) engine fuelled with n-butanol, Energy 85 (2015) 296-303.
[23.] Han, X., Wang, M., and Zheng, M., "An Enabling Study of Neat n-Butanol HCCI Combustion on a High Compression-ratio Diesel Engine," SAE Technical Paper 2015-01-0001, 2015, doi:10.4271/2015-01-0001.
[24.] Xie, K., Yanai, T., Yang, Z., Reader, G. et al., "Emission Analysis of HCCI Combustion in a Diesel Engine Fueled by Butanol," SAE Technical Paper 2016-01-0749, 2016, doi:10.4271/2016-01-0749.
[25.] Mack J.H., Schuler D., Butt R.H., Dibble R.W., Experimental investigation of butanol isomer combustion in Homogeneous Charge Compression Ignition (HCCI) engines, Applied Energy 165 (2016) 612-626.
[26.] Martinez-Frias, J., Aceves, S., Flowers, D., Smith, J. et al., "HCCI Engine Control by Thermal Management," SAE Technical Paper 2000-01-2869, 2000, doi:10.4271/2000-01-2869.
[27.] Yelvington, P. and Green, W., "Prediction of the Knock Limit and Viable Operating Range for a Homogeneous-Charge Compression-Ignition (HCCI) Engine," SAE Technical Paper 2003-01-1092, 2003, doi:10.4271/2003-01-1092.
[28.] Yang, Y., Dec, J., Dronniou, N., Sjoberg, M. et al., "Partial Fuel Stratification to Control HCCI Heat Release Rates: Fuel Composition and Other Factors Affecting Pre-Ignition Reactions of Two-Stage Ignition Fuels," SAE Int. J. Ewgiwes 4(1):1903-1920, 2011, doi:10.4271/2011-01-1359.
[29.] Aichlmayr H., Kittelson D., Zachariah M., Miniature free-piston homogeneous charge compression ignition engine-compressor concept-Part II: modeling HCCI combustion in small scales with detailed homogeneous gas phase chemical kinetics, Chemical Engineering Science 57 (2002) 4173-4186.
[30.] Pinheiro, A., Vuilleumier, D., Kozarac, D., and Saxena, S., "Simulating a Complete Performance Map of an Ethanol-Fueled Boosted HCCI Engine," SAE Technical Paper 2015-01-0821, 2015, doi:10.4271/2015-01-0821.
[31.] Ortiz-Soto E., Assanis D., Babajimopoulos A., A comprehensive engine to drive-cycle modelling framework for the fuel economy assessment of advanced engine and combustion technologies, International Journal of Engine Research 13 (2012) 287-304.
[32.] Bissoli M., Frassoldati A., Cuoci A., Ranzi E., Mehl M., Faravelli T., A new predictive multi-zone model for HCCI engine combustion, Applied Energy 178 (2016) 826-843.
[33.] Agbro E., Tomlin A.S., Lawes M., Park S., Sarathy S.M., The influence of n-butanol blending on the ignition delay times of gasoline and its surrogate at high pressures, Fuel 187 (2017) 211-219.
[34.] Frassoldati A., Cuoci A., Faravelli T., Niemann U., Ranzi E., Seiser R., Seshadri K., An experimental and kinetic modeling study of n-propanol and iso-propanol combustion, Combustion and Flame 157 (2010) 2-16.
[35.] Grana R., Frassoldati A., Faravelli T., Niemann U., Ranzi E., Seiser R., Cattolica R., Seshadri K., An experimental and kinetic modeling study of combustion of isomers of butanol, Combustion and Flame 157 (2010) 2137-2154.
[36.] Nativel D., Pelucchi M., Frassoldati A., Comandini A., Cuoci A., Ranzi E., Chaumeix N., Faravelli T., Laminar flame speeds of pentanol isomers: An experimental and modeling study, Combustion and Flame 166 (2016) 1-18.
[37.] Pelucchi M., Cavallotti C., Ranzi E., Frassoldati A., Faravelli T., Relative Reactivity of [O.sub.x]ygenated Fuels: Alcohols, Aldehydes, Ketones, and Methyl Esters, Energy & Fuels 30 (2016) 8665-8679.
[38.] Heufer K.A., Sarathy S.M., Curran H.J., Davis A.C., Westbrook C.K., Pitz W.J., Detailed kinetic modeling study of n-pentanol oxidation, Energy & Fuels 26 (2012) 6678-6685.
[39.] Carstensen H.-H., Dean A.M., Development of detailed kinetic models for the thermal conversion of biomass via first principle methods and rate estimation rules, Computational Modeling in Lignocellulosic Biofuel Production, Oxford University Press, Oxford, UK 201 (2010).
[40.] Kumar K., Zhang Y., Sung C.-J., Pitz W.J., Autoignition response of n-butanol and its blends with primary reference fuel constituents of gasoline, Combustion and Flame 162 (2015) 2466-2479.
[41.] Ranzi E., Frassoldati A., Stagni A., Pelucchi M., Cuoci A., Faravelli T., Reduced kinetic schemes of complex reaction systems: fossil and biomass-derived transportation fuels, International Journal of Chemical Kinetics 46 (2014) 512-542.
[42.] Pelucchi M., Bissoli M., Cavallotti C., Cuoci A., Faravelli T., Frassoldati A., Ranzi E., Stagni A., Improved kinetic model of the low-temperature oxidation of n-heptane, Energy & Fuels 28 (2014) 7178-7193.
[43.] Pelucchi M., Development of Kinetic Mechanisms for the Combustion of Renewable Fuels, Department of Chemistry, Materials and Chemical Engineering "G. Natta", Politecnico di Milano, 2017, pp. 406.
[44.] Ranzi E., Dente M., Goldaniga A., Bozzano G., Faravelli T., Lumping procedures in detailed kinetic modeling of gasification, pyrolysis, partial oxidation and combustion of hydrocarbon mixtures, Progress in Energy and Combustion Science 27 (2001) 99-139.
[45.] Ranzi E., Frassoldati A., Granata S., Faravelli T., Wide-range kinetic modeling study of the pyrolysis, partial oxidation, and combustion of heavy n-alkanes, Industrial & engineering chemistry research 44 (2005) 5170-5183.
[46.] Darcy D., Nakamura H., Tobin C., Mehl M., Metcalfe W., Pitz W., Westbrook C., Curran H., A high-pressure rapid compression machine study of n-propylbenzene ignition, Combustion and Flame 161 (2014) 65-74.
[47.] Burke S., Simmie J., Curran H., Critical Evaluation of Thermochemical Properties of C1-C4 Species: Updated Group-Contributions to Estimate Thermochemical Properties, Journal of Physical and Chemical Reference Data 44 (2015) 013101.
[48.] Zhou C.-W., Simmie J.M., Curran H.J., Rate constants for hydrogen-abstraction by from n-butanol, Combustion and Flame 158 (2011) 726-731.
[49.] Zhou C.W., Simmie J.M., Curran H.J., Rate constants for hydrogen abstraction by HO2 from n-butanol, International Journal of Chemical Kinetics 44 (2012) 155-164.
[50.] Katsikadakos D., Zhou C.-W., Simmie J., Curran H., Hunt P., Hardalupas Y., Taylor A., Rate constants of hydrogen abstraction by methyl radical from n-butanol and a comparison of CanTherm, MultiWell and Variflex, Proceedings of the Combustion Institute 34 (2013) 483-491.
[51.] Ranzi E., Dente M., Faravelli T., Pennati G., Prediction of kinetic parameters for hydrogen abstraction reactions, Combustion science and technology 95 (1993) 1-50.
[52.] Curran H., Rate constant estimation for C1 to C4 alkyl and alkoxyl radical decomposition, International journal of chemical kinetics 38 (2006) 250-275.
[53.] Curran H., Gaffuri P., Pitz W.J., Westbrook C.K., A comprehensive modeling study of n-heptane oxidation, Combustion and flame 114 (1998) 149-177.
[54.] Zador J., Fernandes R.X., Georgievskii Y., Meloni G., Taatjes C.A., Miller J.A., The reaction of hydroxyethyl radicals with O 2: A theoretical analysis and experimental product study, Proceedings of the Combustion Institute 32 (2009) 271-277.
[55.] da Silva G., Bozzelli J.W., Liang L., Farrell J.T., Ethanol oxidation: Kinetics of the [alpha]-hydroxyethyl radical+ O2 reaction, The Journal of Physical Chemistry A 113 (2009) 8923-8933.
[56.] Welz O., Klippenstein S.J., Harding L.B., Taatjes C.A., Zador J., Unconventional peroxy chemistry in alcohol oxidation: the water elimination pathway, The journal of physical chemistry letters 4 (2013) 350-354.
[57.] Pelucchi M., Ranzi E., Frassoldati A., Faravelli T., Alkyl radicals rule the low temperature oxidation of long chain aldehydes, Proceedings of the Combustion Institute, (2016).
[58.] Pelucchi M., Somers K.P., Yasunaga K., Burke U., Frassoldati A., Ranzi E., Curran H.J., Faravelli T., An experimental and kinetic modeling study of the pyrolysis and oxidation of n-C 3 C 5 aldehydes in shock tubes, Combustion and Flame 162 (2015) 265-286.
[59.] Cuoci A., Frassoldati A., Faravelli T., Ranzi E., OpenSMOKE++: An object-oriented framework for the numerical modeling of reactive systems with detailed kinetic mechanisms, Computer Physics Communications 192 (2015) 237-264.
[60.] Vranckx S., Heufer K., Lee C., Olivier H., Schill L., Kopp W., Leonhard K., Taatjes C., Fernandes R., Role of peroxy chemistry in the high-pressure ignition of n-butanol-Experiments and detailed kinetic modelling, combustion and flame 158 (2011) 1444-1455.
[61.] Weber B.W., Kumar K., Zhang Y., Sung C.-J., Autoignition of n-butanol at elevated pressure and low-to-intermediate temperature, Combustion and Flame 158 (2011) 809-819.
[62.] Black G., Curran H.J., Pichon S., Simmie J.M., Zhukov V., Bio-butanol: Combustion properties and detailed chemical kinetic model, Combustion and Flame 157 (2010) 363-373.
[63.] Stranic I., Chase D.P., Harmon J.T., Yang S., Davidson D.F., Hanson R.K., Shock tube measurements of ignition delay times for the butanol isomers, Combustion and Flame 159 (2012) 516-527.
[64.] Heufer K., Bugler J., Curran H., A comparison of longer alkane and alcohol ignition including new experimental results for n-pentanol and n-hexanol, Proceedings of the Combustion Institute 34 (2013) 511-518.
[65.] Tang C., Wei L., Man X., Zhang J., Huang Z., Law C.K., High temperature ignition delay times of C5 primary alcohols, Combustion and Flame 160 (2013) 520-529.
[66.] Dagaut P., Togbe C., Experimental and modeling study of the kinetics of oxidation of butanol- n-heptane mixtures in a jet-stirred reactor, Energy & Fuels 23 (2009) 3527-3535.
[67.] Bissoli, M., Cuoci, A., Frassoldati, A., Faravelli, T. et al., "Detailed Kinetic Analysis of HCCI Combustion Using a New Multi-Zone Model and CFD Simulations," SAE Int. J. Ewgiwes 6(3):1594-1609, 2013, doi:10.4271/2013-24-0021.
[68.] Komninos, N., Hountalas, D., and Kouremenos, D., "Development of a New Multi-Zone Model for the Description of Physical Processes in HCCI Engines," SAE Technical Paper 2004-01-0562, 2004, doi:10.4271/2004-01-0562.
[69.] Oakley A.J., Experimental investigations on controlled autoignition combustion in a four-stroke gasoline engine, Brunel University School of Engineering and Design PhD Theses, 2001.
[70.] Oakley A., Zhao H., Ma T., Ladommatos N., Dilution effects on the controlled auto-ignition (CAI) combustion of hydrocarbon and alcohol fuels, (2001).
[71.] Oakley, A., Zhao, H., Ladommatos, N., and Ma, T., "Experimental Studies on Controlled Auto-ignition (CAI) Combustion of Gasoline in a 4-Stroke Engine," SAE Technical Paper 2001-01-1030, 2001, doi:10.4271/2001-01-1030.
[72.] Bhave, A., Kraft, M., Mauss, F., Oakley, A. et al., "Evaluating the EGR-AFR Operating Range of a HCCI Engine," SAE Technical Paper 2005-01-0161, 2005, doi:10.4271/2005-01-0161.
+39 02 02 2399 3243
Department of Chemistry, Materials and Chemical Engineering "G. Natta", Politecnico di Milano, Piazza Leonardo da Vinci 32, 20133, Milan, Italy
EGR - Exhaust Gas Recirculation
PRF - Primary Reference Fuels
TRF - Toluene Reference Fuels
RON - Research Octane Number
MON - Motor Octane Number
LHV - Lower Heating Value
CN - Cetane Number
SI - Spark Ignition
CI - Compression Ignition
HCCI - Homogeneous Charge Compression Ignition
GCI - Gasoline Compression Ignition
HRR - Heat Release Rate
PRR - Pressure Release Rate
RCM - Rapid Compression Machine
ST - Shock Tube
BDE - Bond Dissociation Energy
IVC - Intake Valve Closing
EVO - Exhaust Valve Opening
AFR - Air-to-Fuel Ratio
CR - Compression Ratio
CoV - Coefficient of Variation
IMEP - Indicated Mean Effective Pressure
[[eta].sub.comb] - Combustion Efficiency
TDC - Top Dead Center
Matteo Pelucchi, Mattia Bissoli, and Cristina Rizzo
Politecnico di Milano
Yingjia Zhang and Kieran Somers
National University of Ireland Galway
Politecnico di Milano
National University of Ireland Galway
Politecnico di Milano
Table 1. Properties of alcohols, gasolines and diesel fuels. Adapted from Sarathy et al.  and Kalghatgi . Fuel LHV (MJ/L) RON MON CN Gasoline ~30-33 88-98 80-88 n.d. Diesel ~35 n.d. n.d. 40-55 n-butanol 26.9 98 85 12 n-pentanol 28.5 80 74 20 n-hexanol 29.3 56 46 24 Table 2. [PHI]=1.0 alcohols / "air" mixture compositions (% mole fraction) tested in this study. n-butanol Mixture Fuel 02 N2 Ar C02 Mixl-C4,p=30bar 3.38 20.29 68.70 7.63 0.00 Mix2-C4,p=10bar 3.38 20.28 15.26 61.09 0.00 Mix3-C4,p=30bar 3.39 20.36 53.61 0.00 22.63 n-pentanol Mixture Fuel 02 N2 Ar C02 Mixl-C5,p=10bar 2.72 20.43 76.85 0.00 0.00 Mix2-C5,p=10bar 2.72 20.43 38.42 38.42 0.00 Table 3. Kinetic parameters of the lumped oxidation reactions of n-butanol and n-pentanol (units are mol, cm, s, cal). Lumped reactions n-butanol A n Ea R+[O.sub.2][right arrow][O.sub.2] 7.50E+12 0.0 0 R[O.sub.2][right arrow] R+[O.sub.2] 3.00E+13 0.0 30000 R[O.sub.2] [right arrow] Unsat Alcohols 0.80E+37 -7.5 39500 + H[O.sub.2] R[O.sub.2][right arrow]OH+CH2O+[C.sub.n-] 1.00E+10 0.0 22000 aldehyde (Waddington) R[O.sub.2] [right arrow] QOOH 8.36E+9 0.39 19621.3 QOOH [right arrow] R[O.sub.2] 3.55E+05 1.596 11124.06 [beta]-QOOH [right arrow] H[[??].sub.2] + 4.04E+07 1.823 23182.1 Unsat Alcohols [gamma]/[delta]-QOOH [right arrow] OH + 8.19E+15 1.013 23327.45 Olefins QOOH [right arrow] Cyclic Ether 2.78E+10 0.371 17120.7 + OH QOOH [right arrow] [H.sub.2]O + Cn-1 3.00E+10 0 13000 Aldehyde + HCO QOOH + [O.sub.2] [right arrow] 7.50E+12 0.0 0 OOQOOH OOQOOH [right arrow] QOOH + 3.00E+13 0.0 30000 [O.sub.2] aQOOH + [O.sub.2][right arrow]H[[??].sub.2] 3.00E+12 0.0 35000 + Aldehyde OOQOOH [right arrow] OQOOH 8.36E+9 0.39 19621.3 + OH OQOOH [right arrow] OH + 9.00E+15 0.0 41500 Products Lumped reactions n-pentanol A n Ea R+[O.sub.2][right arrow][O.sub.2] 1.00E+13 0.0 0 R[O.sub.2][right arrow] R+[O.sub.2] 4.00E+13 0.0 30000 R[O.sub.2] [right arrow] Unsat Alcohols 1.20E+37 -7.5 39500 + H[O.sub.2] R[O.sub.2][right arrow]OH+CH2O+[C.sub.n-] l.00E+10 0.0 22000 aldehyde (Waddington) R[O.sub.2] [right arrow] QOOH 2.16E+03 2.463 17204.05 QOOH [right arrow] R[O.sub.2] 5.97E+03 2.186 10839.18 [beta]-QOOH [right arrow] H[[??].sub.2] + 4.79E+12 0.48 27344.5 Unsat Alcohols [gamma]/[delta]-QOOH [right arrow] OH + 1.07E+17 1.335 23538.5 Olefins QOOH [right arrow] Cyclic Ether 2.78E+10 0.371 17120.7 + OH QOOH [right arrow] [H.sub.2]O + Cn-1 3.00E+10 0 13000 Aldehyde + HCO QOOH + [O.sub.2] [right arrow] 1.00E+13 0.0 0 OOQOOH OOQOOH [right arrow] QOOH + 4.00E+13 0.0 30000 [O.sub.2] aQOOH + [O.sub.2][right arrow]H[[??].sub.2] 3.00E+12 0.0 35000 + Aldehyde OOQOOH [right arrow] OQOOH 2.16E+03 2.463 17204.05 + OH OQOOH [right arrow] OH + 9.00E+15 0.0 41500 Products Table 4. Ricardo E6 . Engine characteristic and specific model parameters. Displacement 504 Speed [rpm] 1500 [cm.sup.3] Bore [mm] 76.2 Intake Valve Closing [[degrees]ATDC] 137 Stroke [mm] 111.1 Exhaust Valve Opening [[degrees]BTDC] 144 Rod Length [mm] 241.3 [Cu.sub.z] 0.12 Compression Ratio 11.5 [Cu.sub.w] 0.58
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||homogeneous charge compression ignition|
|Author:||Pelucchi, Matteo; Bissoli, Mattia; Rizzo, Cristina; Zhang, Yingjia; Somers, Kieran; Frassoldati, Ale|
|Publication:||SAE International Journal of Engines|
|Date:||Dec 1, 2017|
|Previous Article:||Functional Requirements to Exceed the 100 kW/l Milestone for High Power Density Automotive Diesel Engines.|
|Next Article:||Combustion Indexes for Innovative Combustion Control.|