Abstract

As an effective clean fuel, ethanol has the characteristics of improving antiknock quality and reducing emissions. It is an ideal antiknock additive for Homogeneous Charge Compression Ignition (HCCI) engines. The oxidation of gasoline-ethanol surrogates in HCCI engines is a very complex process which is dominated by the reaction kinetics. This oxidation process directly determines the performance and emissions of HCCI engines. Coupling the computational fluid dynamic (CFD) model with the gasoline-ethanol surrogate mechanism can be used for fuel design, so the construction of a reduced mechanism with high accuracy is necessary. A mechanism (278 species, 1439 reactions) at medium and low temperatures and experiments in a HCCI engine for the oxidation of gasoline-ethanol surrogates were presented in this paper. Directed relation graph with error propagation (DRGEP) method and quasi-steady-state assumption (QSSA) method were used in order to get a reduced model. Then, the kinetics of the vital reactions related to the formation and consumption of H and OH were adjusted. To validate the model, the HCCI experiments for the oxidation of gasoline-ethanol surrogates were conducted under different operating conditions. The verification result indicated that the present model can predict the oxidation process of gasoline-ethanol effectively.

1. Introduction

Due to rapid urbanization and industrialization, pollution levels are increasing at an alarming rate recently. Processing the fuel efficiently plays an important role in reducing the pollution due to fuel emissions, which can play an important role in improving the health of the citizens, especially in urban areas. Gasoline-ethanol is a multicomponent substance of thousands of macromolecular hydrocarbons and it is time-consuming and costly to develop a mechanism for the real fuel. Moreover, the application of the detailed model for complex gasoline-ethanol surrogate fuels in HCCI engine simulations is not practical with current computing resources, due to the large scale and the stiffness of the detailed mechanism. Therefore, the representative components of gasoline should be selected reasonably and the model of multicomponent gasoline surrogates should be reduced while maintaining its good performance.

Primary reference fuel (PRF), the two-component (iso-octane/n-heptane) mixture, is generally considered to be the most common surrogates for gasoline. In recent years, more and more experiments were conducted in HCCI engines under high pressure, medium and low temperatures, and low equivalent ratio conditions, which provided a basis for the application of PRF mechanism in HCCI engine simulations. The PRF oxidation process of “the first oxygen addition ⟶ the first isomerization ⟶ the second oxygen addition ⟶ the second isomerization” is the key section during the autoignition process. In addition, as a commonly used additive, ethanol has become an important component for gasoline surrogate fuels due to its good antiknock performance and low emissions. The chemical kinetic mechanism of ethanol-PRF coupled with CFD software helps to understand the oxidation phenomena of mixture such as autoignition, flame propagation, flameout, combustion stability, and emissions. This is of great importance for further improving combustion efficiency and reducing emissions.

The PRF models given by Halstead et al. [1] and Cox and Cole [2] were empirical models, which were still widely applied in the simulation of autoignition process. Then Li et al. [3] proposed a reduced model for predicting PRF oxidation behaviors, including ignition delay (τ), heat release rate (HRR), and molarity of vital species. This model can well predict the oxidation behaviors of PRF in the low and medium temperature, but the predictions at high temperature phase were not satisfactory. The reduced mechanism developed by Tanaka et al. [4] can be applicable for predicting τ, HRR, and knock in HCCI engines in a wide range, but it was difficult to predict the emission characteristics. However, the PRF mechanism constructed [5] by using the hierarchical expansion method can be used to calculate emissions of PHAs and other pollutants, although it was not accurate in predicting τ under the intake temperature (Tin) range of 300 K ∼ 434 K and the pressure (P) of 4.0 MPa.

The model of Curran et al. [6] can well predict the ignition process on a wider scale of Tins, Ps, and φs. According to the model [6], Ra and Reitz [7] proposed a reduced model, involving 41 species and130 reactions, which may predict in-cylinder pressure (P), τ, and HRR accurately. Then, in order to solve the problem of cross-reactions, Kirchen et al. [8] added the cross-reactions to the models of Tanaka et al. [4] and Marinov [9]. In 2013 and 2015, two mechanisms developed by Liu et al. [10] and Wang et al. [11] were also presented for the combustion of gasoline surrogates.

A detailed three-component (iso-octane/n-heptane/ethanol) model [12] may predict laminar flame speeds (SLs) accurately under high temperature. Then Zheng and Zhong [13] developed a reduced three-component model (50 species, 193 reactions). Its calculated τs were highly consistent with the experimental values. Based on this model, a three-component model [14] was proposed by adding some elementary reactions related to H and updating relevant kinetic parameters, which can predict SLs and τs accurately. Moreover, Lemaire et al. [15] analyzed the effect of the additive (ethanol) to gasoline on the formation of soot. They pointed out that adding 10% ∼ 30% (by volume) ethanol can reduce the production of soot precursor significantly; the amount of reduction for soot was 25% ∼ 81%. In 2019, Li et al. [16] developed a highly reduced four-component gasoline-ethanol model, which may predict the experimental data for PRF, toluene primary reference fuel (TRF), and PRF-ethanol surrogates.

In recent years, we have conducted in-depth studies on the HCCI test and chemical reaction kinetics of related mixtures (Energy, 2019, 169:572–579. Tehnički Vjesnik Technical Gazette, 2020, 27(5):1571–1578. (SCI); Acta Microscopica, 2020, 29(2):720–731. (SCI)). In the field of chemical reaction kinetics, the mechanism of linear alkane and iso-alkanes and the chemical reaction kinetics of paraffin fuels have been constructed. In addition, preliminary achievements have been made in the research of primitive reactions and active groups sensitive to ignition delay.

In summary, the combustion process of gasoline-ethanol blend has attracted more and more attention recently. Many gasoline-ethanol mixture mechanisms have been constructed. However, due to the stiffness caused by long simulation time scale, the existing models are too large in scale and have poor accuracy under the current computing resources. Furthermore, HCCI validation experiments fueled with a hydrocarbon blend or a hydrocarbon-oxygen blend are rare and more experimental data is needed to compare with the calculated value, in order to further verify the reduced model. Therefore, the objective of this paper is to perform HCCI experiments on the combustion of gasoline-ethanol surrogates and to develop a smaller size mechanism for the lean gasoline-ethanol surrogates (φ < 1) by implementing a reduction and merge scheme using DRGEP and QSSA methods. In order to carry out extensive validation of this reduced model, not only the calculation used by the proposed model should be compared with the results of the HCCI experiments, but also the new model should be compared with the previous literature models.

2. Kinetic Modeling

In this section, the processes of the initial mechanism construction, the automatic chemistry mechanism reduction for gasoline-ethanol surrogates, and the determination of the final mechanism are presented, as displayed in Figure 1. The reduced gasoline-ethanol mechanism was constructed by first reducing and then merging.

Firstly, in order to get a mechanism of smaller size, DRGEP method was used to eliminate the insignificant species efficiently, and QSSA method was used to identify the species that were in quasi-steady-state. Following the above process, the two submechanisms (the models of ethanol [9] and PRF [6]) were reduced. Taking the reduction of PRF submechanism as an example, according to the PRF reaction path given by Ra and Reitz [7], initial reactants (iso-octane/n-heptane/O2), intermediate components (C2H3, CH3, CH2CHO), and final products (CO2, H2O) were selected as target substances. The larger the error threshold was set, the smaller the reduced mechanism scale would be and the prediction accuracy would decrease. In the process of reduction, a smaller threshold is set and reduction is carried out several times to ensure the prediction accuracy. Subsequently, a powerful and accurate merge for the reduced submechanisms from three disparate fuels was conducted.

Secondly, several relevant reactions involving H and OH in the products or reactants were revised, and the parameters of these reactions were adjusted. Finally, the final mechanism (278 species and 1439 reactions) was proposed.

Thirdly, the repeated reactions and components of PRF and ethanol mechanisms are mainly small molecule reactions of C1∼C3, H, and O2. However, rate constants of the same reactions in these two mechanisms are different, leading to great changes in the generation or consumption rates of many free radicals (H, OH, H2O2, HO2, etc.). The rate constants of the small molecular reactions of C1∼C3, H, and O2 in PRF submechanism were selected.

2.1. Base Mechanism and Case Settings

In order to keep the consistency between predictions and experimental results on reaction rate, transport, and thermodynamic data, the detailed chemical kinetic mechanism of PRF (1034 species, 4236 reactions) [6] and the semidetailed ethanol model (57 species, 383 reactions) [9] were taken as base mechanisms.

Mechanism [6] has been used in the simulations of HCCI engines in [17, 18]. The ethanol submechanism has been validated by comparing with the experimental SLs and τs. The predictions by ethanol submodel were also compared with mole fractions of species measured in stirred and flow reactors. In this paper, by combining complex calculation with recently obtained rate constants, the ethanol submechanism was updated by adjusting the Arrhenius coefficients of some elementary reactions, as shown in Table 1 and [1922].

The reduction work of the mechanism for gasoline-ethanol surrogates was carried out on the conditions targeted for HCCI engines. HCCI engines generally operate at low φs, so the flame temperature should be lower compared with conventional internal combustion engines. To reduce the chemical model, 18 conditions were chosen at different temperatures (500 K∼1000 K) and different φs (0.25, 0.3 and 0.5). The three-component fuel (iso-octane/n-heptane/ethanol = 62%: 18%: 20% by volume) should be considered as a substitute for ethanol-gasoline fuels [23], and it was used as gasoline-ethanol surrogates in this paper.

All the related calculations were performed by using CHEMKIN package. Assuming that mass transport is infinitely fast, the gas phase reaction is controlled solely on the nature of each species, not by transport constraints. The whole domain has uniform thermodynamic and transport properties.

For nonadiabatic cases, heat transfer between the cylinder and the wall was noticed. Related parameters of heat transfer were set according to the Woschni formula [24]. The relevant data in the model were set according to the specification of the test engine bench and the running conditions. Specifications for the test engine were presented by Zhang and Wu [25].

2.2. Mechanism Reduction

The reduction of the mechanism was implemented by DRGEP and QSSA approaches.

The DRGEP method is used for initial reduction. For a detailed mechanism, when specific species (reactants, reaction products, and important intermediate species) are set as the target component, then a series of species are strongly coupled to the target component. When an error threshold is set, DRGEP method can identify unimportant components and thus remove those components and the reactions associated with them. The DRGEP method was implemented efficaciously by removing the species and reactions whose target variable error in the worst case exceeds the threshold. The worst-case error was considered to be the maximum relative error of the original and reduced mechanisms in the 24 target cases. When the worst-case error was within the threshold, the tolerance would be increased and the first stage reduction would continue. The reduction of the first phase was terminated until the worst-case error of simulation results exceeded the preset threshold (5%). At this stage, some species (ethanol, i-C8H18, n-C7H16, s-C2H4OH) was chosen to be the starting species. Relative tolerance of mole-fractions for CO, CO2, and H was 0.009, respectively, and threshold for iso-octane, s-C2H4OH, CH3CHO, n-heptane, and ethanol was set to 5%, respectively. After reduction, the first-stage model (296 species, 1691 reactions) was developed.

Subsequently, QSSA method was also applied. The QSSA method [2628] can be used to identify some intermediate species whose production and consumption rates were nearly equal. For these species, the change in concentration was almost negligible. After setting a threshold, these intermediate quasi-steady-state (QSS) species would be processed by a nonlinear algebraic system. In this paper, the relative tolerances of HRR and ignition delay were set to 0.1 compared with the initial mechanism; the relative tolerance of mole fractions for CO, iso-octane, and n-heptane were set to 7%. iC4H8 and C5H3 were identified under quasi-steady-state. It not only saved computational time, but also greatly reduced the stiffness of the ordinary differential equation system. After second-stage reduction, the model (277 species, 1437 reactions) was constructed.

3. Experimental Setting

In this paper, the HCCI experiments were performed to validate the mechanism of gasoline-ethanol surrogates and to provide more basic data on combustion characteristics of the test fuel. The three-component fuel (iso-octane/n-heptane/ethanol = 62%: 18%: 20% by volume) (a surrogate for 95 RON gasoline) was used as the test fuel in HCCI experiments. The selected HCCI operating conditions are shown in Table 2.

3.1. Experimental Setup

The test engine was retrofitted based on a water cooled, direct injection, naturally aspirated, original engine, CT2100Q. The first cylinder maintained the conventional diesel engine mode, while the second cylinder operated in HCCI mode.

In order to meet the requirements of HCCI operating mode, the intake system, exhaust system, and fuel system were modified. The details and schematic of the test HCCI engine system can be found in [25]. In order to control the intake temperature (Tin), an independent port-fuel-injection system and an electric heating system were installed on the second intake pipe. When the injection pulse width was bigger than 2.5 milliseconds, the test fuel was injected into the HCCI cylinder. The injection timing was set to 30°CA BTDC and duration can vary from 5 milliseconds to 8 milliseconds under different operating conditions. The P of the second cylinder was recorded by Kistler 6052A. This piezoelectric pressure makes it possible to record digital signals (TTL: >4, 5 V high, <1 V low level). These signals are then transmitted, either by ECU or statically from the test cylinder, and then it was analyzed by Kibox 283A. The accuracy for the instruments employed can be seen in [25].

3.2. Experimental Procedure

To start the test engine in HCCI mode smoothly, at the beginning, the engine worked in diesel mode. When the water temperature reached 95°C and oil temperature reached 85°C, the diesel supply to the first cylinder was stopped and the fuel injection for the other cylinder was started simultaneously. As a result, the operating mode was successfully switched. When the HCCI engine run stably, the in-cylinder pressures were recorded, averaged, and analyzed based on 100 consecutive cycles.

4. Results and Discussion

τ, SL, HRR, P, and molarity for vital species are the key parameters of the prediction.

4.1. Laminar Flame Speed

The simulation analysis was performed in flame speed reactor. Figure 2 gives SLs of three single-component fuels at Tin = 298 K and Pin = 0.1 MPa. The comparison of SLs between the results of calculation obtained by using the second-stage model and experimental results under different φs can be seen in Figure 2(a). It can be found that the trends of the calculated SLs with φs are consistent with the experimental values in [14]. That is, as φ increases from 0.6 to 1.3, the SLs of the three fuels first increase and then decrease. However, there is a significant gap between the two curves. In other words, the calculated SLs for the three single-component fuels cannot match with the experimental values. Therefore, the second-stage mechanism should be revised.

By analyzing the sensitivity of the base model on SL at selected operating conditions, the reactions that affect SL obvious were picked out. These highly sensitive main-chain branching reactions were revised and subsequently several reaction rate constants were also adjusted to improve the prediction. In addition, the two following reactions keeping highly sensitive were also added, based on the detailed discussion shown by Westbrook et al. [29]:

Any elementary reaction that produces hydrogen (H) radical increases the rate of branch reaction R9. R21 not only increases the production of H, but also affects the destruction of hydroxyl (OH) radical, so these two elementary reactions take vital parts in the autoignition process. Moreover, this stage reduced model cannot predict the reactions related to H accurately. As a result, the parameters of these reactions related to H require revision.

Based on detailed model developed by Mehl et al. [30], some Arrhenius coefficients were adjusted. The elementary reactions associated with H radical may accurately reproduce the characteristics of iso-octane flame, so vital reactions that affect SL were extracted, as shown below:

In addition, OH radical has an obvious effect on autoignition of gasoline-ethanol surrogates [14]; the parameters of the reactions involving H or OH in the products or reactants should be subjected to sensitivity analysis, and then the kinetics of the vital reactions need to be adjusted.

As a result, the final model including 278 species and 1439 reactions was proposed. To determine the prediction accuracy of this model, SLs of the three initial fuel components were calculated. Figure 2(b) shows the calculated value and the experimental data under different φs. It is obvious that final mechanism can predict SLs of the three initial components more accurately than second-stage model. This also indicates that the revision of the model is effective.

4.2. Pressure and Heat Release Rate

To further study the oxidation process of gasoline-ethanol surrogates in the HCCI engine, in-cylinder pressure (P) and heat release rate (HRR) in HCCI mode were calculated by coupling CHEMKIN with CFD software. The given HCCI operating conditions were as follows: n = 1200 r/min, Pin = 0.1 MPa, Tin = 423 K, and different φs (φ = 0.25, 0.30, 0.35, 0.40).

The HCCI combustion chamber model based on HCCI engine parameters is given in Figure 3.

Coupling CHEMKIN with CFD, calculated values and experimental results for in-cylinder pressures and HRRs were compared, when Pin = 0.1 MPa, n = 1200 r/min, φ = 0.25, 0.3, 0.35, 0.4, as shown in Figure 4.

As can be seen from Figure 4, curves of experimental values and the calculated data show the same trend. Firstly, under the above conditions, when equivalent ratio φ rises, the two peak values of P and HRR obtained by simulation show an increasing trend, and the corresponding timing of the peak value is advanced. Obviously, variation trend of the two parameters (P and HRR) obtained by simulation is consistent with the experimental data. Secondly, the calculated data are in good agreement with the experimental values.

P rises slightly as the piston moves up. After autoignition, P rises sharply until it reaches the peak in-cylinder pressure (Pmax), and then P gradually decreases as the piston moves down during the power stroke. Inevitably, there are two factors that may cause the predicted Ps to be higher than the experimental values. Firstly, the calculation is performed by a zero-dimensional model that ignores crevices and the inhomogeneity of mixture concentration and temperature in the cylinder. Secondly, the assumptions are closed, constant volume and adiabatic.

In addition, when Tin rises from 413 K to 433 K, the formation and combustion e are speed up, the formation of OH and related active groups are accelerated, so the ignition delay is shortened. However, the Pmax at high Tin (433 K) is smaller than that at low temperature (413 K). This is because Tin is too high, resulting in a reduction in the density of the mixture, which in turn reduces the quantity of injection fuel per cycle.

4.3. Heat Release Rate

The calculated HRR at two different φs (0.3, 0.4) using this model and the model by Li et al. [16] are compared with HCCI experimental data in Figures 5(a) and 5(b).

The given conditions were Pin = 0.1 MPa, Tin = 433 K and n = 1200 r/min.

The simulation curve (blue line) obtained in this study follows the trend: when φ increases, the heat release in the cylinder is more concentrated, and the shape of the HRR curve will change to a narrow and high trend. That is, when φ increases, the peak heat release rate increases accordingly, and the time of occurrence is advanced. This is mainly because, as φ increases, the number of activated molecules increases and the thermal energy in the cylinder is more sufficient, which not only leads to faster overall reaction, shorter ignition delay period, and concentrated heat release, but also increases the heat release in the cylinder.

Due to the inhomogeneity, the HRRs calculated by this model and the model of Li et al. [16] are significantly greater than the HCCI experimental data. As mentioned earlier, for the same reasons, the calculated Ps are also higher than the experimental results.

Above all, the predicted trends using this model are in consistency with the HCCI experimental results.

4.4. CA10 and CA50

In engine experimental research, τ is the duration from the opening of the injector needle valve until the moment when the P curve starts to separate from the pure compression curve at the compression process. This moment when P rises sharply refers to the CA corresponding to the heat release percentage being 10% (CA10). Moreover, when 50% of the heat has been released is usually considered as the midpoint of the combustion process per engine cycle, which was marked as CA50.

To further study the oxidation process in HCCI cylinder, CA10s and CA50s at different Tins by using the final model, the HCCI results are compared in Figures 6 and 7. As can be seen from Figure 6, when Tin decreases from 433 K to 393 K, CA10 is in advance about 6°CA ∼12°CA. Figure 7 shows that the midpoint of combustion process (CA50) at Tin of 393 K is 6°CA ∼12°CA later than at Tin of 433 K.

Above all, the results indicate that Tin play a vital part in oxidation process of the test fuel. Good agreements can be achieved between the calculation and the experimental data.

Comparison of the experimental data (scatters) [31] and calculated results (lines) for PRF100/PRF90/PRF0 in a rapid compression machine (RCM) can be seen in Figure 8. It shows the relationship between τ and φ. As expected, the higher the proportion of n-heptane in PRFs is, the shorter the τ is.

Moreover, the higher Tin and the higher T during the compression stroke cause the mixture to be more homogeneous. As a result, more active radicals may be generated and the reaction rate may be accelerated, which lead to a shorter τ.

Overall, a higher intake air temperature will advance the phase of the peak heat release rate, which means that the heat release of the fuel in the HCCI cylinder will be more concentrated; that is, when Tin is appropriately increased, CA10 and CA50 advanced and the HRR curve tends to be narrow and high.

4.5. Mole Fractions of the Vital Species

Comparison of the predicted mole fractions for iso-octane, heptane, and ethanol by our developed model and Li et al. model [16] at different Tins can be seen in Figure 9.

The following findings can be drawn from Figure 9. Firstly, under the same operating condition, the moments at which the three initial components begin to decrease sharply predicted by the two models are relatively close, as shown by the shaded area in Figure 9. This indicates that the τs predicted by the two mechanisms are roughly the same. Secondly, when Tin increases from 393 K to 433 K, the moments at which the three initial components begin to decrease sharply are advanced from −8°CA to −14°CA. This is because the formation of free radicals may be accelerated at high Tin, thus leading to the advance of autoignition.

In summary, the final model is predictive at the autoignition phase by comparing the predicted , P, SL, HRR, T, CA10, CA50, and mole fractions of vital species under selected conditions.

5. Conclusions

This work developed a mechanism (278 species, 1439 reactions) for gasoline-ethanol surrogates by implementing the reduction and merge scheme. DRGEP and QSSA methods were used to efficiently reduce the mechanism. Moreover, the kinetic parameters of the relevant reactions related to the formation and consumption of H and OH were adjusted.

HCCI experiments were conducted on the combustion of gasoline-ethanol surrogates. More HCCI experimental data were provided to validate the models for the oxidation.

The proposed mechanism was validated as well as the predictions of the previous literature model. Since the calculation is based on ideal assumptions, there is a gap between the simulated curve and the experimental curve. Overall, the prediction of this developed model was found satisfactory in terms of certain characteristic parameters involving SL, P, T, CA10, CA50, , HRR, and mole fractions of species under the selected HCCI conditions.

Based on the reduced mechanism and HCCI engine model developed in this paper, we will further analyze the influence of other boundary conditions (such as intake pressure, engine speed, and EGR rate) on the flame structure and combustion flow field of fuel combustion in the future.

Nomenclature

CA10:Timing at which 10% of the heat has been released (@ºCA)
CA50:Crank angle at which 50% of the heat has been released (@ºCA)
Tin:Intake temperature (K)
φ:Equivalence ratio

Abbreviation

BTDC:Before top dead center
CFD:Computational fluid dynamics
TRF:Toluene primary reference fuel.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by the Special Fund Chang’an University (300102228509, 300102228403, 300102228505, 300102229502, and 300102229202).