Energy Optimization and Fuel Economy Investigation of a Series Hybrid Electric Vehicle Integrated with Diesel/RCCI Engines

Among different types of low temperature combustion (LTC) regimes, eactively controlled compression ignition (RCCI) has received a lot of attention as a promising advanced combustion engine technology with high indicated thermal efficiency and low nitrogen oxides ( NO x ) and particulate matter (PM) emissions. In this study, an RCCI engine for the purpose of fuel economy investigation is incorporated in series hybrid electric vehicle (SHEV) architecture, which allows the engine to run completely in the narrow RCCI mode for common driving cycles. Three different types of energy management control (EMC) strategies are designed and implemented to achieve the best fuel economy. The EMC strategies encompass rule-based control (RBC), offline, and online optimal controllers, including dynamic programing (DP) and model predictive control (MPC), respectively. The simulation results show a 13.1% to 14.2% fuel economy saving by using an RCCI engine over a modern spark ignition (SI) engine in SHEV for different driving cycles. This fuel economy saving is reduced to 3% in comparison with a modern compression ignition (CI) engine, while NO x emissions are significantly lower. Simulation results show that the RCCI engine offers more fuel economy improvement in more aggressive driving cycles (e.g., US06), compared to less aggressive driving cycles (e.g., UDDS). In addition, the MPC results show that sub-optimal fuel economy is achieved by predicting the vehicle speed profile for a time horizon of 70 s.


Introduction
Two-thirds of the oil consumption in the world is currently used in the transportation sector and half of that goes to passenger cars and light trucks [1].Prevalent consumption of the petroleum-based fuels leads to high greenhouse gas (GHG) emissions.The transportation sector in the U.S. and Europe accounts for around 25% of the total GHG emissions.In the U.S., the goal is to decrease the transportation GHG by 35% by 2025 [2] and in Europe the goal is to cut the transportation GHG by 67% by 2050 [3].In this context, automakers must reduce the GHG emissions by introducing advanced fuel-efficient technologies and also by using alternative fuels.low temperature combustion (LTC) engines include a family of internal combustion engines (ICEs) technologies, including premixed charge compression ignition (PCCI), homogeneous charge compression ignition (HCCI), and reactivity controlled compression ignition (RCCI) engines, which offer low engine-out NO x and low soot emissions [4], and peak net indicated thermal efficiency as high as 53% [5].Moreover, improvement in ICE fuel efficiency has the largest potential in improving hybrid electric vehicle (HEV) fuel economy improvement and reducing GHG emissions, compared to enhancement in other HEV component efficiency [6].Thus, integrating fuel-efficient LTC engines in HEVs has the potential to improve the vehicle fuel economy and decrease the GHG emissions.
However, utilizing the LTC engines in vehicles faces two major challenges: (i) limited engine operating range; and (ii) control complexity in mode transitions (e.g., SI ↔ LTC, CI ↔ LTC).To tackle these two challenges, this paper investigates the integration of an LTC engine with a series hybrid electric vehicle (HEV) by taking advantage of decoupling the ICE from the drivetrain.This allows the LTC engine to operate in a narrow operating range, removes the engine mode transients, and simplifies the LTC engine control.In addition, SHEVs are already available on the market, such as the Chevrolet Volt, Fisker Karma, and BMWi3.
Another important factor in improving the fuel economy of HEV powertrains is developing optimal energy management control (EMC) strategies to maximize fuel saving.Figure 1 lists the major EMC approaches in HEVs.These include (1) a rule-based controller (RBC) such as fuzzy [7] and on-off [8] strategies.These strategies are implementable in real-time applications due to their robustness and low computational cost.However, the RBC strategies do not offer the best HEV fuel saving [9,10] due to their offline design; (2) Global optimization strategies such as pontryagin minimum principle (PMP) [11][12][13] and dynamic programming (DP) [14][15][16][17].These strategies require the complete information of the driving cycle to determine the optimum EMC.While these global optimization strategies cannot be applied in real-time, the solution from these strategies can be used as a platform to find the ultimate fuel saving for evaluating other EMC strategies [9]; (3) Real-time optimization controllers such as stochastic dynamic programming (SDP) [18,19], equivalent consumption minimization strategy (ECMS) [20][21][22], and model predictive controllers (MPC) [23][24][25].These sub-optimal EMC strategies can be implemented in real-time.Prior HEV studies are divided into three groups according to the ICE type (Figure 1).In the first group, CI engines have been used in different HEV architectures [8,14,18,20,23,[26][27][28].The CI engines have been mostly used in sport utility vehicles (SUVs), trucks, and buses [7,18,[27][28][29].In the second group, SI engines have been used in HEV [11,12,30,31].Recently, SI engines with Atkinson cycle have become more popular in the market (e.g., Ford C-Max, Honda Accord PHEV, Lexus RX 450h, and Toyota Prius).In [32], an Atkinson cycle SI engine was used by Toyota R&D group to increase the fuel economy benefit of HEVs.In another study [33], a Honda Accord PHEV was redesigned based on the Atkinson cycle and 10% lower fuel consumption was reported compared to the SI engine.
Few studies are found in the literature that investigated integrating LTC engines with HEVs.Among different types of LTC engines, HCCI was the first type explored in hybrid electric powertrains.In the first study at Argonne National Laboratory in U.S., the effects of using a dual mode SI-HCCI engine in different vehicle electrification levels were analyzed [34].Their simulation results predicted from 6% to 15% fuel consumption reduction, depending on powertrain configurations and driving cycles.In another study [35], fuel economy improvement of a HCCI engine versus a SI engine for both mild and medium parallel HEVs was investigated.The authors reported a range from 17% to 35% fuel economy improvement with using a dual-mode SI-HCCI in comparison to a conventional SI non-hybrid powertrain.In both studies [34,35], a dual mode (SI-HCCI) engine was used with a rule-based EMC.Another study [21] was conducted to investigate the effect of utilizing a dual mode SI-HCCI engine on a power-split HEV acceleration performance.The authors quantified the trade-off between the vehicle fuel economy and the vehicle performance.In [36], we carried out the first study on utilizing a pure HCCI mode engine in a series hybrid powertrain and we found a 12.6% improvement in fuel economy in comparison with a SHEV running with a SI engine.In our next study [37], we investigated the impact of engine startup fuel penalty on SHEV HCCI-based powertrains.Our results show that the HCCI fuel economy improvement compared to the SI engine is independent of the engine start/up fuel penalty amount.
RCCI is the second type of LTC engine that has been recently investigated for HEVs.The study in [38] is the only study available in the literature for an RCCI-based HEV.In [38], researchers at the University of Wisconsin-Madison and Oak Ridge National Laboratory in U.S. used an RCCI engine in a series hybrid electric powertrain.The authors tested the vehicle for the US environmental protection agency (EPA) highway fuel economy test (HWFET) procedure to measure the vehicle fuel economy and emissions.Based on simulation results, the authors predicted that a series-parallel RCCI-HEV configuration will lead to a 12% fuel economy improvement for the Chevrolet Volt, which currently uses a SI engine.The study in [38] included preliminary results and no model-based EMC strategy was used for optimizing the energy balance between the battery and ICE energy sources.In addition, the results were presented only for the HWFET driving cycle.The EMC type and driving cycles will affect the RCCI-HEV performance and fuel saving results, as will be shown in this paper.
Moreover, in [39] we carried out the first study on integrating a multi-mode LTC engine with a SHEV architecture.The multi-mode LTC engine was able to switch between the HCCI, RCCI and conventional SI modes by incorporating a fuel penalty.Our results show a 9% to 13.1% fuel economy improvement, compared to an identical SHEV platform running with a single-mode SI engine.
To the best of the authors' knowledge, this is the first study undertaken to investigate the fuel economy benefit of integrating an RCCI engine in SHEV with advanced EMCs.The contribution from this study is threefold.First, it investigates the ultimate fuel saving of a dedicated RCCI-mode engine in SHEV configuration.Second, it investigates the effect of EMCs on the potential of fuel economy improvement for SHEV using an RCCI engine.Third, it studies the impact of driving cycle and also the battery initial SOC on the fuel economy of the RCCI-based SHEV.
This paper is organized as follows.In Section 2, the developed HEV model is presented and the main components of the HEV model are described.Section 3 explains three different EMC strategies, including on-off RBC, global optimization DP and MPC.In Section 4, the optimum number of engine operating points are selected with regards to the fuel economy.The SHEV results are presented in Section 5. Finally, Section 6 summarizes all findings from this paper.

SHEV Model Description
A SHEV with specifications listed in Table 1 is modeled in this section.The SHEV model encompasses different submodels including the RCCI engine, longitudinal vehicle dynamics (LVD), E-motor, and battery.The model is a forward-in-power HEV model developed in Matlab /Simulink to evaluate the EMC strategies and fuel economy benefits.Using the LVD model, the vehicle speed profile is calculated based on the available supplied traction torque by the E-motor after subtracting for drag and rolling resistance forces.The supervisory controller specifies the battery required power and the engine power based on the battery state of charge (SOC), and the driver power demand.Figure 2 shows the high level schematic of the model.The description for sizing and selection of HEV electrical components are found in [37].The SHEV model is briefly explained as follows.

IC Engines
Three different ICEs will be studied as part of the SHEV.These ICEs include RCCI, CI (diesel), and SI (gasoline) engines.The RCCI engine used in this work is based on a modified GM Z19DTH diesel engine.The diesel engine was converted to an RCCI engine by a group of scholars from the University of Wisconsin-Madison and Oak Ridge National Laboratory [40].The RCCI engine has a new piston design and its compression ratio is reduced from 17:1 to 15:1 [40].Table 2 lists the RCCI engine specifications and Figure 3 shows the RCCI engine Brake Specific Fuel Consumption (BSFC) map using the data from [40].The engine speed range is from 1000 rpm to 3500 rpm.The BSFC map of the engine (Figure 3) was created using the engine brake thermal efficiency (BTE) data from [40] and lower heating value (LHV) of fuels.The RCCI engine model is designed such that it takes the required brake mean effective pressure (BMEP) and the required engine power and then calculates the engine fuel consumption.For the comparative study, the BSFC model of a CI engine and a SI engine are included in the HEV model.The same base engine is used for both RCCI and CI (diesel) modes.Thus, the original GM Z19DTH diesel engine data from [40] is used for the CI engine model.Figure 4 shows the BSFC map for the CI engine.For a fair comparison, the SI engine is selected such that it has optimum power rating for the HEV size in this study.To this end, a GM A14XFL SI (gasoline) engine from a mid-size HEV on the market is chosen.Figure 5 shows the SI engine's BSFC map based on the data from [41].This figure demonstrates that the selected engine operating points (i.e., points # 1, 2, 3) are located in the lowest BSFC region of the SI engine.

Longitudinal Vehicle Dynamics
The LVD model calculates the vehicle actual speed V based on the vehicle traction and resistance forces and the dynamics associated with the vehicle.The following equation shows the governing dynamic model used for the LVD model: where M is the vehicle total mass; F gravity , F drag , and F roll are gravity force, aerodynamic drag force, and rolling resistance force, respectively.The resistive forces are calculated from the following equations: where θ is the road slope, C d is the vehicle drag coefficient, ρ is the air density, A is the vehicle frontal area, and f r is the rolling resistance coefficient.The main parameters of the LVD model are listed in Table 1.

E-Motor Model
A single gear ratio 75 kW E-motor [42] is utilized in this study.The E-motor requested power is calculated based on the supervisory controller's demanded power as the input to a map-based E-motor model.The E-motor model includes an efficiency map versus power request and rotational speed up to the maximum of 7500 rpm.Based on the E-motor operating point, the E-motor actual generated power is obtained from the efficiency map.Then, the E-motor traction torque is calculated using the electrical available power and the motor speed.The available electrical power consists of the E-motor generated power and the battery power by including a constant auxiliary load used for the entire driving cycle [43].The E-motor efficiency map covers both traction (η traction ) and regenerative braking (η braking ) modes.Equation (5) shows the torque relation of the E-motor: where, P available refers to the available electrical power in the battery.Consequently, the vehicle traction force (F traction ) and the E-motor speed (ω motor ) are calculated by: ) where, R w is radius of the wheels and Gr is the final gear ratio between the E-motor and the driveline.The Gr is designed to match the vehicle maximum designed speed with the E-motor maximum speed (i.e., the maximum vehicle speed is 100 mph and the maximum E-motor speed is 7500 rpm).
The actual vehicle speed is obtained from the LVD model.In addition, the E-motor was sized and selected based on the required traction torque and power for the vehicle gradeability and acceleration performances [37].Table 1 lists the main specifications of the E-motor.

Battery Model
Based on the E-motor operating point, the E-motor actual generated power is obtained from the efficiency map.The supervisory controller specifies the battery required power, as shown in Figure 2. A battery energy capacity and power is sized to meet the SHEV design requirements.The battery supplies 40 miles AER, which is a common range for SHEVs such as 2015 Cadillac ELR.A battery model is employed to estimate the battery SOC and the battery voltage.The battery model is from [44] and it is widely used in previous HEV and EV investigations.This model extracts a nonlinear model (experimental data [45]) for the battery characteristic discharge curve.Finally, the SOC is calculated by: where i is the battery current and Q is the battery nominal capacity.More discussions on the battery sizing and models are available in [37].

Vehicle Acceleration and Braking Performance
The HEV model is implemented in Matlab /Simulink.The vehicle acceleration performance was designed to reach 60 mph in 12.2 s, which is similar to a common mid-size HEV on the market [46].The vehicle braking performance was designed so as to be able to stop from 60 mph in less than 53.4 m, similar to that in [46].

Model Validation
The BSFC maps of the RCCI and CI engines are developed experimentally in [40] and the experimental data are used in this study.The models for two other HEV major components, including the E-motor and battery were validated against experimental data [37].The selected E-motor is UQM PowerPhase 75 and the E-motor efficiency model was validated against the manufacturer's experimental data [42] with RMSE less than 1.5% [37].Moreover, the battery model was validated against the experimental data from SAFT VL7P Li-ion battery [45] with the RMSE of 0.17 V [37].

Energy Management Controller (EMC) Design
The EMC strategies in this work aim to (i) keep the hybrid powertrain to operate in the charge sustaining (CS) mode; (ii) fulfill the driver power demand; and (iii) maximize the powertrain fuel saving [47].Here, three types of EMC strategies are designed for the SHEV-RCCI powertrain.These strategies include on-off RBC, DP, and MPC.A desired SOC window of 0.55 SOC 0.9 is used in all of the three EMCs.

EMC Type I: RBC
In on-off RBC, the EMC rules are designed heuristically without driving cycle information.The battery SOC is the only input to the control unit, which forces the battery to keep the SOC in the desired window (i.e., 0.55 SOC 0.9), by controlling the engine on/off status.The on-off RBCs cannot adapt their rules with changing driving cycles.This results in non-optimal efficiency for a wide operating range.Simple implementation for real-time EMC applications is the major advantage of on-off RBC strategies [48].Here, a heuristic RBC is designed to keep the battery SOC within the desired window.

EMC Type II: Global Optimization-DP
In a SHEV, the engine power profile is determined by an optimization-based EMC while minimizing a cost function for fuel/energy consumption.If the driving cycle is known, meaning that the driver power request is known, then a global optimal solution can be found [48].While this approach cannot be applied in real-time EMC applications, this method can provide an ideal baseline to assess different EMCs.In this work, the DP method is employed for structuring an optimization problem to find optimum strategies for minimizing a performance index.By doing the calculation backwards over the time horizon based on Bellman's Principle of Optimality [49], DP searches for the best control action among all the possible actions offline in time by assuming the entire driving cycle information is available.
The engine requested power is the control variable (u k ) and the battery SOC is the state variable in this formulation.The control variable is discretized by a time step of 1 s.An energy based performance index consisting of the fuel energy and the battery energy is formed in Equation ( 9): where, the index j represents the feasible transitions to the next time (t + where, η f path and η bpath represent the combined efficiency in the ICE path and the battery path for the transition between the kth states in the successive times.P bpath and P f path are the battery power and the engine produced power, respectively.m f ,start in Equation ( 10) is the fuel penalty for each engine start-up.The cost associated with the engine start-up is incorporated in Equation ( 10) by introducing the N kj constant which is equal to 1 when the engine is at start-up and it is equal to 0 during the rest of the engine operation.Q LHV is the gasoline fuel lower heating value.By finding the engine optimum requested power (u k ) to fulfill the driver power request, the battery produced power is calculated according to the following constraint at each state: In the backward DP, the optimal cost-to-go from the current time (t) to the end of the driving cycle is defined as: where, S * (t, u k ) is the optimal cost-to-go from the kth state at the current time t to the end of the driving cycle.S(t + 1, u k ) represents the optimal cost in the next time (i.e., t + 1) to the end of the driving cycle.At each time a state (k) which has the minimum cost among the different states as determined by the DP strategy as the optimal control variable.A nonlinear backward HEV model, which does not include the driver model, is used to increase flexibility for the real-time implementation of the optimal EMC model.This backward model assumes that the vehicle tracks exactly the driving cycle; thus, the vehicle power demand is directly calculated from the driving cycle.In addition, a high fidelity forward Simulink vehicle model is designed for assessing the EMC strategy and analyzing the HEV performance.The purpose for including DP results in this paper is to present the ultimate energy saving using the RCCI hybrid electric powertrain.The DP values serve as a benchmark for the comparison with RBC and MPC results with different time horizons.

EMC Type III: MPC
The MPC concept is deployed to form the EMC optimization problem on a moving receding horizon.Given that the vehicle speed can be predicted over a short time horizon (TH), the DP strategy is used over a short horizon to find a sequence of the sub-optimal control strategy.The MPC strategy is developed by assuming that the future driving cycle information over the TH is provided from the GPS data.The cost function (J) at the nth TH is shown in Equation (14).J is minimized over TH by selecting optimal ICE/generator power request.
subject to: P e,min < P e < P e,max P m,min where, the subscripts max and min denote maximum and minimum, respectively.The first term in Equation ( 14) refers to the fuel energy consumed by the ICE and the second term refers to the battery electrical energy consumed or recharged during the driving cycle.The subscripts e and m in Equation ( 15) denote engine and motor, respectively.The DP formulation in Section 3.2 is used over the TH to calculate the optimal ICE/generator power at time step n.A closed-loop MPC is designed to reject disturbances such as sudden changes in the estimated driving cycle data.However, this EMC strategy does not provide a globally optimal solution, but it can be used for real-time implementation.
In the current formulation, the time domain is discretized into one-second intervals.The solutions consist of local optimum control signal at each time step.

Selection of Engine Operating Points
In SHEVs, the engine can operate independently of the vehicle speed and the wheel torque.Thus, the engine has the opportunity to work most of the time at high efficiency operating points and low engine-out emission region.In this section, the engine operating points are selected for three different constant power levels as low, mid, and high to fulfill the vehicle power requirements and also to guarantee low engine BSFC and low engine-out emissions.The low BSFC and low engine-out emission constraints are discussed in subsections A and B.
A. Engine BSFC Constraint: The selected operating points are at three power levels designated as 10 kW, 20 kW, and 30 kW.Within a range of 10 percent for each power level, a search is done to find the lowest BSFC value.Figure 3 shows the three final selected operating points for the RCCI engine.In the single operating mode, the mid power level (i.e., point no. 2 in Figure 3) is selected since the engine can provide the mean power requirement for the vehicle at this operating point.For the 2-mode operation, a low level power and a mid level power are selected for the engine operating points (i.e., points no. 1 and 2 in Figure 3), and for the 3-mode operation all the three selected points are considered.

B. Engine Emissions Constraint:
The engine HC, CO, and NO x emissions are considered in selecting the engine operating points.CO and HC conversions in diesel oxidation catalysts (DOC) are a function of the exhaust gas temperature.Thus, the engine operating points are selected to meet the minimum DOC light-off temperature.The RCCI engine exhaust gas temperature for different engine speeds and torques are shown in Figure 6 based on the data from [40] along with the selected operating points.It is shown in Figure 6 that the exhaust gas temperature for all of the three engine operating points are above 290 • C.Moreover, in reference [50], it is shown that the DOC for the same RCCI engine achieved 90% and 100% HC and CO conversions, respectively, when the exhaust gas temperature in the RCCI engine was higher than 290 • C. Thus, the selected RCCI engine operating points meet the DOC light-off temperature to achieve low HC and CO emissions, similar to those in the CI and SI engines.In addition to the HC and CO emissions, the NO x emission is also taken into account for selecting the engine operating points.Table 3 lists the NO x emission associated with the selected operating points in both RCCI and CI engines.The table shows that the selected engine operating points in RCCI produce much less NO x compared to the selected operating points in the CI engine.This is a well-recognized fact in literature [5,40,51] and considered as one of the main advantages of RCCI engines.Thus, a smaller NO x after-treatment system is required in RCCI engines, compared to CI engines.Given that the SI and CI engines are production engines, the engines and exhaust after-treatment systems are optimized to meet HC, CO, and NO x constraint for the emission standard.In addition, the discussions above demonstrate that the RCCI engine produces low HC, CO, and NO x emissions in the three selected operating points.Thus, all the three engine operating points in the SI, CI, and RCCI engines have both low BSFC and low emissions advantages over the other engine operating regions at the designated power levels.
In this section, the effect of the number of the engine operating points on the HEV fuel economy is also discussed.Table 4 shows the effect of number of engine operating points and the engine start-up fuel penalty on the RCCI-SHEV fuel economy.The simulation results in Table 4 show that the RCCI fuel economy increases by increasing the number of engine operating points.This can be explained by improved overall ICE efficiency when three modes are utilized in RCCI for different driver power request levels.Given better fuel economy results by using the 3-mode EMC, the rest of the results in this paper are presented for the 3-mode EMC.

Results and Discussion
In this section, the results for the SHEV with the RCCI, CI, and SI engines are discussed.In Section 5.1, sensitivity of the driving cycle's prediction time horizon on the vehicle's fuel economy is studied.Moreover, the initial battery SOC effect on the vehicle's fuel economy is investigated in Section 5.2.Lastly, the effect of the type of the driving cycle is presented in Section 5.3.
In this work, a combined driving cycle (Figure 7) consisting of three standard driving cycles including UDDS (urban dynamometer driving schedule), HWFET (highway fuel economy test), and US06 is used to test the EMC strategies.The combined driving cycle is the base driving cycle for all the analyses in this paper except for Section 5.3, where different driving cycles are compared.Figure 8 describes the effect of the engine start-up fuel penalty (m f ,start ) on the fuel saving.The RCCI fuel economy improvement is constant versus the amount of m f ,start .This is because the power levels for each engine are similar; thus, the number of engine on/off switching will remain constant.This makes the RCCI fuel economy improvement independent of the m f ,start value.Figure 8 also shows the same trend for the RCCI fuel saving over the number of engine operating points and the fuel economy improves with increasing number of the engine operating points.Figure 9 illustrates the E-motor speed variation for the combined driving cycle.A single gear is designed to synchronize the E-motor and vehicle maximum speeds in order to expand the E-motor operating points to the whole E-motor operating region.In Figure 9, for the portion of the US06 driving cycle that the vehicle speed reaches to 80 mph, the E-motor speed is 6000 rpm.The E-motor operating points during the combined driving cycle are depicted in Figure 10.The E-motor efficiency points range from 70% to 90% in the traction and regenerative braking modes.The E-motor efficiency is higher at higher E-motor power (i.e., P e > 40 kW) and the E-motor efficiency improves with increasing the E-motor power.Figure 11 shows the vehicle tracking performance.The vehicle is able to follow the reference driving cycle with a root mean square error (RMSE) of 0.8 mph.This result also confirms that the HEV components' sizing meets the performance requirements during the driving cycle.Thus, the HEV model can be used as a testbed for evaluation of EMCs.  Figure 12 compares the engine on/off status and the battery SOC variation for both RBC and DP EMCs.The RBC strategy regulates the engine status command based on the SOC value.The RBC switches the engine on when the SOC reaches to its lowest allowed value (i.e., SOC min = 0.55).In the DP controller, the SOC variation is much less than that of the RBC.This allows the EMC to choose the most optimal control strategy over the driving cycle at each time.However, in the DP strategy, the engine switches on and off more frequently (Figure 12b).Figure 13 indicates the power distribution of the battery and the engine during the driving cycle for the DP control strategy.There are different operating modes, such as fully electric mode and opportunity charging mode based on Figure 13.For instance, the battery only supplies the driver power demand at t = 1282 s and runs the vehicle on the fully electric mode.For low power demands, at t = 3005 s, the vehicle runs in the opportunity charging mode, in which the engine/generator supplies the power demand and charges the battery simultaneously.In RBC, the engine is turned off during the braking, but in the DP strategy, the engine can either be on or off to minimize the cost function.In addition, at the low vehicle speeds, the mechanical braking assists the regenerating braking to supply the braking torque, however, at the high vehicle speeds, all the braking torque is supplied by the regenerative braking.The SOC variations over the driving cycle for the DP and MPC controllers are illustrated in Figure 14.The EMC is in the charge-sustaining mode and the EMC strategy is tuned to enforce the final SOC value to match the initial SOC value at the end of the driving cycle.Moreover, Figure 14 compares the SOC profile for each EMC.The SOC variation for different prediction TH is also shown.The results in Figure 14 show that all the designed EMCs can sustain SOC at the initial SOC by the end of the driving cycle.The designed DP EMC algorithm in this work is computationally efficient.
The simulation processing time on a 2.20 GHz Intel processor is about 72.6 s for the whole combined driving cycle (i.e., 11.5 ms per one second of the driving cycle).

Prediction Time Horizon Sensitivity
Figure 15 demonstrates the HEV fuel economy sensitivity with the driving cycle prediction TH for different engines.It is observed that the fuel economy increases by knowing more driving cycle information (i.e., larger TH) and it reaches a plateau, in which at TH = 120 s the HEV fuel economy merges to the DP results.There is a trade-off between the computational cost and fuel saving by the MPC strategy.Higher TH improves the HEV fuel economy at the cost of needing more computation time.In these simulations, THTH = 70 s offers a good compromise between computation load and fuel economy improvement.These results show that by knowing a short time horizon information of the vehicle's speed profile (i.e., TH = 70 s), 95% of the global optimal fuel economy (i.e., DP) is achieved.This figure also compares different engines' fuel economy variations with the EMC types.The RCCI engine offers the highest fuel economy compared to SI and CI engines.In the DP EMC, the RCCI engine has 5.9 and 1.1 MPG greater fuel economy than the SI and CI engines, respectively.

Initial SOC Sensitivity
In this section, the effect of running the RCCI engine in a lower initial SOC on the fuel economy is studied.The low SOC can represent operation in an extended range electric vehicle (E-REV).It is assumed the battery is operating at a lower initial SOC value (i.e., SOC initial = 0.3) during the charge-sustaining mode.Fuel economy of the HEV integrated with RCCI, SI, and CI engines for different EMCs are shown in Table 5. Comparing Figure 15 and Table 5 shows that the fuel economy is greater when the SOC initial is higher.In the lower SOC initial scenario, the battery has greater losses than that in higher initial SOC since the battery loss is greater in the low SOC region.Figure 16 shows the battery energy losses in both the scenarios.It can be seen that the SHEV with lower initial SOC has 31 kJ more battery energy losses than the larger SOC case in the combined driving cycle.The battery energy losses should be compensated by running the ICE for a longer time, which leads to lower fuel economy in the low SOC initial scenario.Figure 17 shows the engine produced power profile in both low SOC initial SOC and high SOC initial scenarios.The number of ICE on/off switching is identical in the two cases, but the engine runs for a longer time in the low initial SOC case.Therefore, the RCCI engine in this condition has more opportunity to work and consequently more fuel saving is achieved compared to the SI and CI engines.This is shown in Figure 18, where the SHEV with low SOC initial has higher fuel economy improvement (Fuel economy improvement is calculated by × 100 where FE HEV−RCCI is the fuel economy for the RCCI-HEV powertrain, and FE HEV is the fuel economy for SI-HEV or CI-HEV powertrain) compared to the high SOC initial scenario.

Driving Cycle Effect
Here, the effect of using various driving cycles is investigated for the RCCI fuel economy improvement for both SHEVs with low initial SOC and high initial SOC scenarios.Four different driving cycles including US06, HWFET, UDDS, and the combined driving cycle are studied.Figure 19 shows the effect of the various driving cycles on the HEV fuel economy improvement by incorporating the RCCI engine versus conventional SI and CI engines.The fuel economy improvement of the RCCI engine versus the SI engine varies from 13.1% to 14.2% depending on driving cycles.Moreover, the fuel economy improvement versus the CI engine varies from 1.8% to 3.0%.Table 6 lists the required average power of the HEVs in each of the studied driving cycles.By comparing Table 6 and Figure 19, it becomes clear that the average power of a driving cycle has a direct effect on the RCCI fuel economy saving over the SI and CI engines.In high power demand driving cycles (e.g., US06), the engine needs to run longer to charge the battery and compensate for the higher power demand.This gives the RCCI engine more opportunities to save more fuel since it is generally more fuel-efficient than conventional SI and CI engines.In addition, the initial SOC effect in Figure 19 is related to the difference in the length of the engine runtime, as previously mentioned in Section 5.2.

Conclusions
In this work, an RCCI engine was integrated with a SHEV powertrain.The fuel economy benefit of the RCCI-HEV powertrain was compared with the conventional ICE-HEV powertrains.Three different ICEs were studied, including a GM Z19DTH engine in CI and RCCI modes and a GM A14XFL SI (gasoline) engine.For a meaningful comparison among the different powertrains, the ICE produced the same amount of power for each powertrain.In selecting engine operating points, the engine-out emission constraints were considered.A high-fidelity forward-in-power SHEV model was developed in Matlab /Simulink with experimentally validated submodels.Three different types of EMC strategies, including RBC, DP, and MPC were developed to investigate the effect of the control strategy on the potential fuel saving from an RCCI-based HEV.The following summarizes the main findings from this work:

•
The simulation results showed that the RCCI engine offers significant potential for fuel saving in SHEV architecture.In the combined driving cycle, integrating an RCCI engine with a SHEV powertrain provided up to 12.6% higher fuel economy over a modern SI engine, while the improvement over the CI-based HEV was 2.2% for the combined driving cycle.

•
An increasing number of the engine operating points can increase the RCCI-HEV fuel economy (FE) improvement by 2.1% by utilizing more engine operating points.This study did not consider fuel penalty during ICE transients between engine modes.It is anticipated that excessively frequent mode switching is not desirable for fuel saving once a large number of engine operating points is used.

•
The 3-mode DP provides the best FE improvement among the EMCs studied.Using the 3-mode DP leads to 17.0% more FE improvement compared to the RBC EMC.This demonstrates the importance of designing optimal EMC strategies to maximize the FE improvement of the RCCI-HEV.

•
In the MPC strategy, more available information from the driving cycle (larger prediction time horizon) naturally leads to better fuel economy improvement and the MPC results eventually merge to the DP optimal controller results.A prediction horizon length of 70 s showed a compromise between the computational cost and fuel economy improvement in this work.

•
Among the four driving cycles studied, the driving cycle that has higher average power ( P) has higher RCCI fuel saving.Thus, the RCCI-HEV operating in the US06 driving cycle ( P = 19.1 kW) offered 14.2% fuel economy improvement versus 13.1% in the UDDS driving cycle ( P = 5.2 kW) in the SHEV.This was because the engine needs to run for a longer time to compensate for the higher power demand.This leads to more opportunity for the RCCI engine to save more fuel compared to the SI and CI engines.• RCCI-SHEV operation with a battery at low SOC in charge-sustaining mode resulted in a slightly higher fuel economy improvement for the RCCI engine in comparison to a higher initial SOC scenario.This is because the battery in a low SOC region has a larger energy loss which means that the engine is required to run for a longer time to compensate for the higher battery losses.This provides more opportunity for the RCCI engine to save more fuel in comparison to the SI and CI engines.
This study did not include a cost analysis of the HEV powertrain using the RCCI engine versus SI and CI engines in SHEV architecture.Cost analysis versus fuel saving will be the subject of future studies.

Figure 1 .
Figure 1.Prior studies categorized based on different types of internal combustion engines (ICEs) and energy management control (EMCs) incorporated in hybrid electric vehicle (HEVs).

Figure 2 .
Figure 2. A schematic diagram to illustrate the power and control flows among the main components of an RCCI-based SHEV.

Figure 3 .
Figure 3. BSFC map (g/kWh) of the RCCI engine in this study.The data from [40] was used to generate this figure.

Figure 4 .
Figure 4. BSFC map (g/kWh) of the GM Z19DTH CI (diesel) engine.The data from [40] was used to generate this figure.

Figure 5 .
Figure 5. BSFC map (g/kWh) of GM A14XFL SI (gasoline) engine.The data from [41] was used to generate this figure.

Figure 6 .
Figure 6.RCCI exhaust gas temperature map for selecting the RCCI engine operating points.The data from[40] was used to generate this figure.

Figure 7 .
Figure 7.The combined driving cycle used for the evaluation of the designed energy management control (EMCs).

Figure 8 .
Figure 8. Fuel economy (FE) improvement by using RCCI versus (a) SI and (b) CI engines as a function of number of engine operating points and engine start-up fuel penalty.(EMC: DP).

Figure 9 .
Figure 9. Traction E-motor speed during the combined driving cycle.

Figure 10 .
Figure 10.E-motor operating points, shown with 'o' symbol, plotted over the E-motor efficiency map.

Figure 11 .
Figure 11.Desired vehicle speed vs. actual vehicle speed with a root mean square tracking error of 0.8 mph.

Figure 12 .
Figure 12.Engine on/off status and the battery SOC pattern in the 3-mode DP and RBC strategies.(a) Engine on/off status and (b) the battery SOC pattern in the 3-mode DP and RBC strategies.

Figure 13 .
Figure 13.Power distribution in the SHEV during the combined driving cycle (EMC: 3-mode DP).

Figure 16 .
Figure 16.Comparison between battery energy loss for the low SOC initial and high SOC initial operating conditions (EMC: 3-mode DP).

Figure 17 .Figure 18 .
Figure 17.The RCCI engine power pattern in the SHEV for two initial SOC operating conditions (EMC: 3-mode DP).

Figure 19 .
Figure 19.HEV fuel economy (FE) improvement by using the RCCI engine versus the SI and CI engines in four different driving cycles (EMC: 3-mode DP).
[1,2,3,4] index k is the control variable indicator.k is a finite number and its size is equal to the number of possible values for the digitized control variable.In this work k ∈[1,2,3,4], where k = 1 represents the engine off situation, and k = 2, 3, or 4 represents the engine's three selected power levels (modes).Furthermore, ∆E f ) are the energy consumptions of the ICE path and the battery path respectively for the jth transition between two states.α is an equivalent factor to equate the electrical usage of the battery to virtual fuel consumption.α is a constant number and is sensitive to driving cycles.In this paper, j (t, u k ) α is specified offline to enforce the battery to operate in a charge-sustaining mode.∆E f j (t, u k ) and ∆E b j (t, u k ) are calculated by the following equations:

Table 3 .
[40]ne-out NO x emissions in RCCI and CI engines operating points.The data is from[40].

Table 4 .
RCCI engine fuel economy values in the SHEV architecture as a function of engine start-up fuel penalty and number of operating modes utilized in the engine (EMC: DP).

Table 5 .
Fuel economy comparison for different EMCs including the effect of MPC time horizon (TH) in the SHEV with SOC initial = 0.3 (EMC: 3-mode).

Table 6 .
Average power in the studied driving cycles.