1. Introduction
Improving fuel efficiency and reducing the greenhouse emissions of vehicles has become crucial in the automotive industry. A wide range of technologies is used to this end. In [
1], simple models are developed to estimate the subsystems’ steady-state performances and efficiency of the entire vehicle. The models are combined to simulate the vehicle behavior for different driving cycles and assess the fuel economy for the effect of design variations and new technologies. In [
2], fuel saving was achieved by adjusting vehicle height, using an electromechanical system which results in reduced aerodynamic resistance of the vehicle.
An interest in electrification of the powertrain, i.e., using hybrid electric vehicles (HEV), as a means of improving fuel efficiency and minimizing greenhouse emissions is growing rapidly [
3]. By adding an electric motor (EM) in parallel to the internal combustion engine (ICE), it is possible to recover the braking energy, store it in the battery, and use it later during the traction phase. Furthermore, ref. [
4] introduced an ELPH (electrically peaking hybrid) control strategy that offers EM power usage for acceleration and deceleration phases, whereas ICE provides mean load power of the driving cycle. The main aim of this logic is to use ICE at high speeds with smooth torque demand.
A great deal of research and practical work has attempted to model the performance of the HEVs for optimizing their efficiency. In most of these works, the vehicle simulation models utilized the steady-state maps of the ICE to estimate fuel consumption [
3,
5]. It has been observed that measured fuel consumption data shows a higher value with respect to one calculated using the steady-state maps. This is due to the higher fuel consumption of the engine during transient phases. In [
6], the measured fuel consumption was compared to that calculated with steady-state fuel consumption maps as a function of measured torque and speed of the ICE. Testing two vehicles, such as a Chevrolet Silverado with a 4.3 L Ecotec engine and a Ford F-150 with a 2.7 L EcoBoost engine, they concluded that the difference, due to transient fuel consumptions, was in the range of 2–4%. Furthermore, they established that more aggressive driving conditions could increase this range due to increased transient contribution. Certainly, this difference depends not only on the driving cycle but also on the vehicle parameters and the engine type. Various engines with 3.8 Lto 8 L displacement were studied in [
7], demonstrating that on different drive cycles, the fuel enrichment due to transient phenomena varied by 2–8%.
A relatively wide range of literature describes different approaches to model the engine transient behavior that estimates precisely the fuel consumption under different operating conditions and for different applications [
6,
7,
8,
9,
10,
11,
12]. The authors of [
6] considered four transient events, such as powertrain torque management (during the gear upshift phase for shafts synchronization), tip-in demand, deceleration fuel cutoff, and cylinder deactivation. The influence of each event on the transient fuel rate has been analyzed using experimental data. The corresponding adjustments have been applied to the steady-state fuel consumption to obtain the transient one. In [
8], to produce an engine map that includes transient behavior, the engine torque and speed range was divided into equal bins. Then, the experimentally measured fuel consumption data was allocated to the bins to determine a representative value for each bin. The probability distribution function was analyzed and based on this data, the corresponding engine map was created. Fuel consumption prediction was improved relative to steady-state calculation in [
9] by applying the correction factor, which is function of torque, speed, and power change rates. In [
10], the vehicle fuel consumption was separated into two modes: cruising at a constant speed and acceleration operating modes. In each of these modes, the fuel consumption was calculated based on the instantaneous engine efficiency, approximated using an analytical function rather than the fuel consumption map. In [
11], a real-time prediction of the fuel consumption was evaluated as first order linear function that varies with VSP (vehicle specific power). The analysis in [
12] used vehicle speed (v) and acceleration (a) as dynamic variables to provide the transient correction to be added to the steady-state predicted value. The sign of VSP was used to choose a form of the correction factor.
The model proposed by [
12] was used as a reference in [
13] and further improved taking into account the computational time and accuracy of the model. In [
14], engine torque and rotation speed are used as dynamic variables. Four different cases were considered, depending on the sign of engine torque and rotation speed change rate. A correction factor was then applied to multiply the steady-state fuel consumption. The approach was convenient to compute the fuel consumption of a vehicle in transient conditions using a backward model.
The abovementioned transient corrections were performed on a vehicle with conventional powertrains to improve the accuracy of fuel consumption computation using steady-state maps. As to the authors’ knowledge, no attempt was made to study the improvement of the control strategy of the HEVs by integrating the transient correction of the engine fuel consumption.
In HEV applications, the improved accuracy of fuel consumption computation allows for further optimization of the torque split between EM and ICE. Equivalent fuel consumption minimization strategy (ECMS) is an online optimization-based strategy. Due to its online feature, the ECMS controller can define an optimal ratio of torque split at each instant of time, hence acting as a real-time controller. The ECMS was first applied by Paganelli et al. to HEV [
15], where a concept of conversion of the battery energy to equivalent fuel consumption was introduced. The equivalence factor had a substantial effect on the overall efficiency of the HEV. A genetic algorithm was used in [
16] for optimal selection of the equivalence factor. ECMS allows for penalizing battery use in case of its undesired state of charge (SOC), temperature, and degradation to maintain the desired operating conditions of the battery [
17,
18,
19]. However, the penalization of the ICE use in transient operations cannot be found in the literature.
To fill this gap, this paper proposes an approach to reduce the fuel consumption of HEVs in transient phases by imposing a transient correction factor (TCF) proposed in [
14] in the control strategy. The objective function to be minimized included TCF as a function of a torque changing. The electric machine can deliver torque quickly and with a minor efficiency variation over the operating range. Hence, the best mode to manage the torque transients is to rely more on the EM. An advantage of this approach is the possibility to reduce extra fuel consumption of the ICE in transient phases. Hence, the contributions of the paper comprise the model of the P2 HEV that includes the transient fuel correction and the developed control strategy to reduce fuel consumption by exploiting the e-motor more in transient phases.
The paper is organized as follows. First, the backward model of the HEV longitudinal dynamics and fuel consumption model is briefly described. Then, the transient fuel consumption model and its experimental validation based on experimental data available online [
20,
21] are reported. The TCF is introduced in the ECMS objective function. The results obtained with and without TCF are compared for the UDDS driving cycle.
2. Vehicle Backward Model
In this paper, analysis on transient fuel consumption was carried out based on the backward model of a Mazda CX9 powered by a conventional gasoline powertrain. The model was then validated with experimental data available from the Argonne National Laboratory (ANL) [
20]. The model of HEV was arranged by integrating ECMS, EM, and battery submodels using a P2 mild HEV configuration as shown in
Figure 1. The description of the Mazda CX9 conventional vehicle and HEV backward model validation was discussed in [
18]. The reader is kindly asked to refer to the work [
18] for a detailed description of the model and its submodels.
The backward model estimates the required torque from the vehicle speed imposed by the drive cycle [
3,
4]. Having the vehicle speed and its acceleration as inputs from the drive cycle, the required force on the vehicle tire patch
was calculated as:
where
M is the vehicle mass,
is the aerodynamic resistance force, and
is the tire rolling resistance. The road gradient forces were neglected as the road is supposed to be horizontal.
The effect of rotational masses, such as tire, crankshaft, and flywheel inertias, were taken into account in the model by adding the inertial torque values in corresponding sub-blocks.
Using vehicle configuration and transmission system parameters (
Table 1), the required torque
and angular speed
were calculated based on the force and linear speed of the vehicle.
3. Conventional Vehicle Fuel Consumption Model and Its Validation
Once the required torque and speed were evaluated, the engine fuel usage based on a steady-state map was estimated using a lookup table with the linear interpolation method. The steady-state map of the Mazda CX9 SkyActiv engine is shown in
Figure 2 based on the experimental data from [
21].
The results of the backward model and the experimental tests for fuel consumed during different drive cycles are reported in
Table 2. The results contain three cases of fuel consumption: (1) measured experimentally, (2) calculated using steady-state engine map and experimentally measured engine torque
and rotation speed
as an input, and (3) calculated using the steady-state engine map and model estimated values of torque and speed. The difference between fuel consumption obtained from the experiments (
) and the steady-state fuel consumption
based on measured engine torque (
) and angular speed (
) is in the range of 3.7–5.2%. This difference is due to transient characteristics of the engine that the steady-state maps do not cover. The steady-state fuel consumption
was determined by using the estimated torque and the speed in the backward model. The difference between
and
is in the range of 1–2.4%, which is mainly due to the error in estimating the engine torque and speed.
4. Transient Fuel Consumption Model and Its Validation
To improve the accuracy of the fuel consumption calculation from the backward model, TCF was used to calibrate the steady-state fuel consumption rate of the model (
) with the experimental one (
). As mentioned above in this paper, the transient model was developed based on the approach described in [
14]. In that work, three different representations for TCF were proposed. The most accurate fuel usage estimation was obtained when TCF was represented by the following expression [
14]:
Regression coefficients
can be estimated separately for the following use cases:
In [
14], it was established that coefficients
A,
B,
C and
D have an order of exponential magnitude −5, −5, −4 and −3 for all the cases, respectively. This means that coefficient
D is the dominant one. That indicates that the torque change rate mostly affects the transient fuel consumption increase. Therefore, TCF can be approximated with good accuracy using the following simplified form of Equation (2):
In the current work, coefficient
was set to one (
), so that that TCF was applied only if there was a torque change. Otherwise, TCF was 1, which means that the predicted fuel consumption was equal to the one computed using the steady-state map. Regression analysis was carried out using experimental data on UDDS, US06, and HW driving cycles available on the ANL database [
20]. The calculated values of
from the backward model was used to determine the torque change rate and fuel consumption
. Then, TCF was calculated as a ratio of measured
and steady-state fuel consumption
. Depending on the sign of the
changing rate, the data were split into 4 cases. The data sampling time during the experiments was 0.1 s and a moving average filter of order 10 was used for smoothing. Moreover, operating points with low distribution density were excluded from the analysis. In the case shown in
Figure 3a, the few operating points with a torque change rate higher than 200 Nm/s were omitted as outliers. The values of coefficient
D derived from this linear regression analysis are shown in
Figure 3 and summarized in
Table 3. These values of regression coefficient
D are aligned with those obtained in [
14].
The fuel consumption obtained from the steady-state map can then be corrected with TCF as:
Finally, the developed model with TCF was simulated for three drive cycles.
Figure 4 shows plots for
,
and
on the UDDS cycle. It can be seen that after introducing TCF, the fuel consumption result of
is more accurate with
compared with steady-state prediction.
Figure 5 demonstrates the cumulative sum of the fuel rate:
, and
over the cycle. The results show that overall fuel consumption with transient correction
and experimental total fuel usage (
were matched better. Thus, it can be concluded that the developed TCF model allows for effectively including the transient characteristics of the engine.
Table 4 summarizes the overall fuel consumption results for different driving cycles. It shows a good fit of
with
. In addition, the last column indicates the difference between
and
, the sum of all extra fuel usage during the transient phases.
5. ECMS with Steady-State Map
Equivalent consumption minimization strategy was applied to the system as a supervisory controller which defines the optimum ratio of energy split between EM and ICE. This controller is a real-time controller not requiring preliminary knowledge about the driving cycle. At each time step, the ECMS decides the power split depending on the current situation. Hence, it gives a local minimum because it minimizes the objective function at each instant. The objective function
is the sum of the fuel consumption rate of ICE
and the equivalent consumption rate of EM
as [
15,
17,
18]:
where,
can be derived by calculating engine fuel used to generate the same power provided by EM considering all the efficiencies of the power flow path.
Equivalent factor
:
where
LHV is the lower heating value of the fuel,
is the average efficiency of the ICE,
is the average efficiency of the EM,
is the average efficiency of battery and
is the average efficiency of the inverter. To reduce calculation time the equivalent factor
can be evaluated by using mean efficiencies. However, the fuel economy model can be improved even more by using specific efficiencies of each operation point.
The fuel consumption rate of the engine
is approximated as a polynomial with independent variables of engine torque and speed
from the engine steady-state map of
Figure 6.
The experimental data are approximated by a polynomial fit carried out by means of MATLAB curve fitting toolbox with a good quality
[
18].
To have a fair comparison of the fuel consumption of the HEV system the SOC, of the battery should maintain at the same level at the beginning and the end of the driving cycle. This is usually referred to as the charge sustaining operation mode of HEV. Therefore, the equivalent fuel consumption of the EM is penalized by an S-shaped correction function that varies with SOC. The expression suggested in [
14] was utilized in the objective function as shown in Equations (9) and (10). The ECMS controller with the penalization factor and the operational constraints are described in detail in Equations (8) and (11) [
18].
The objective function in the ECMS controller is described as:
where
is the penalization factor on
SOC [
17]:
And normalized battery
[
17]:
The torque split logic between the electric motor and the internal combustion engine must satisfy the following constraints:
The cost function in Equation (8) is two-variables (
and
) nonlinear optimization problem. However, the problem can be reduced to a single variable type by means of the first constraint described in Equation (11). Therefore, a built-in MATLAB function fminbnd solver can be used [
22]. The above control strategy is first implemented for comparison purposes with objective functions which consider the steady-state fuel consumption rate
without TCF. In
Figure 7 the ICE torque, EM torque, and requested torque at the gearbox input are plotted. The values of the torques at a time instant 375 s are indicated using data tips, which show that the sum of ICE and EM torques is equal to the requested torque.
Figure 7d shows that a charge sustaining mode is realized by the controller as the SOC at the beginning and end of the driving cycle are at the same level. For the UDDS driving cycle, the overall fuel consumption is about 620.7 g.
6. ECMS with Transient Correction Factor
As discussed above, in transient, the fuel consumption is corrected with the TCF of Equation (4). The introduction of the TCF in the objective function of ECMS penalizes the torque changing rate for ICE, which leads to the ICE operating more steadily than the ECMS based on steady-state fuel consumption only, as in the previous case.
Then the objective function with TCF can be introduced to Equation (8):
Figure 8 reports the values of the torques when TCF is considered and the charge sustaining condition is satisfied. The overall fuel consumption was reduced to 616.3 g, which corresponds to an additional 4.4 g of fuel savings with respect to ECMS with a steady-state map.
Figure 9 shows the torque values of the ICE and EM, and fuel consumption rate values for the two ECMS cases: with and without TCF. The zoom-in to a specific region is shown on the right plot. It can be seen that the engine torque
is smoother for the case with TCF (blue solid line). Due to the influence of TCF in the time interval of 55–57 s, the controller chooses a hybrid mode that minimizes the engine’s high transient operation. Instead, in the case when TCF is not considered (dashed orange line) the pure ICE traction mode is chosen by the controller in the same time range. On the last subplot (
Figure 9c), the fuel consumption rate shown highlights the reduction of the fuel consumption rate in the considered time interval. Hence, as the charge sustaining constraint is satisfied, the cumulative fuel consumption is also reduced to the values stated before.
The values of total fuel consumed on three driving cycles using ECMS controller with and without TCF are summarized in
Table 5. Because US06 is the most aggressive driving cycle in terms of transients, the largest benefit of 5.1 g fuel saving is achieved with it. In the HW driving cycle, which is characterized by steady driving scenarios instead, the fuel-saving is reduced to 3.2 g. In UDDS, the reduction is 4.4 g. This improvement can be increased further to 8.9 g by doubling battery capacity (14s24p). The increased utilization of electric traction leads to needs on a larger battery.