Author Contributions: Tobias Nüesch and Alberto Cerofolini developed the theoretical framework for the control of real driving NOx emissions of a diesel HEV. In addition, they performed the simulations, they conducted the experimental validation, and they wrote the article. Giorgio Mancini played a key role in the setup and calibration of the internal combustion engine and the test bench. He was also partially involved in the drafting of the article. Christopher Onder was the technical supervisor, while Nicolò Cavina and Lino Guzzella were the project supervisors. All supervisors were involved in exchanging ideas on the development of the theoretical framework, its practical application, and/or the experimental validation. Moreover, all supervisors participated in the internal review of the article draft.

This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Motivated by the fact that the real driving NOx emissions (RDE) of conventional diesel vehicles can exceed the legislation norms by far, a concept for the control of RDE with a diesel parallel hybrid electric vehicle (HEV) is proposed. By extending the well-known equivalent consumption minimization strategy (ECMS), the power split degree of freedom is used to control the NOx emissions and the battery state of charge (SOC) simultaneously. Through an appropriate formulation of the problem, the feedback control is shown to be separable into two dependent PI controllers. By hardware-in-the-loop (HIL) experiments, as well as by simulations, the proposed method is shown to minimize the fuel consumption while tracking a given reference trajectory for both the NOx emissions and the battery SOC.

Light-duty diesel vehicles are known for their low fuel consumption, as compared to gasoline vehicles. However, due to legislative restrictions, vehicle manufacturers continuously have to make considerable efforts to reduce the pollutant emissions of diesel vehicles. Although the legislative limits have been continuously reduced over the last decade, the real driving emissions, which are the emissions emitted during every-day driving, can far exceed the legislative limits, even for Euro 6 certified light-duty vehicles, as shown in several studies [1–5]. One reason is that the homologation of the vehicles is performed on well-defined, but unrealistic driving cycles. The manufacturers unavoidably focus the optimization effort on such types of vehicle operating conditions. To reduce the discrepancy between the certified and the real-world pollutant emissions, the European commission is currently discussing measures to limit real driving emissions [5].

One option to cope with such a radical change would be to continuously monitor and control the pollutant emissions by an appropriate exhaust aftertreatment system. Another option is provided by electric hybridization of the vehicles, which not only offers a reduction of pollutant emissions, but also a simultaneous reduction of the CO2 emissions. Since hybrid electric vehicles (HEVs) have an additional degree of freedom for the control of the energy flows in the powertrain, the trade-off between fuel consumption and pollutant emissions can be further influenced.

Some studies can be found in the literature about the control of pollutant emissions for HEVs. For example, the authors of [6] propose a real-time rule-based strategy to optimize both fuel economy and pollutant emissions, taking into account cold-start emissions, by minimizing an overall normalized impact function. Similar approaches have been presented in [7–11], where an instantaneous optimization algorithm, with a similar structure to the well-known equivalent consumption minimization strategy (ECMS) [12,13], is built. In that case, the target is to minimize a weighted sum of multiple factors, for example fuel consumption and NOx, CO and CO2 emissions, while guaranteeing charge-sustaining conditions for various driving cycles. The weighting factors between the various components of the target cost function are constant and considered to be tuning parameters. For a diesel HEV equipped with a selective catalytic reaction (SCR) system, a noncausal extended ECMS is proposed in [14], including the minimization of tailpipe emissions, while considering the cold start behavior. A control framework with three state variables arises from the energy management extended with emissions management; these are the energy stored in the battery, the SCR catalyst temperature and total NOx tailpipe mass, resulting in a controller with an unstable co-state, which can be used in a fixed time window only.

Dynamic programming (DP) has also been applied to address the problem of building a supervisory control system for the fuel and emission reduction. Examples are given in [15,16] for a parallel HEV, where the gearshift strategy and the engine start/stop decision are optimized along with the torque split factor and, in [17], considering the power split as the only control input. In both studies, constant weighting factors for the multiple emission sources have been applied.

A general approach based on optimal control theory is proposed in [18], with a description of several possible extensions of the basic framework of ECMS. The authors describe how to include different pollutant components for an HEV, possibly taking into consideration thermal effects and aftertreatment systems. The authors claim that the solution of such a general problem is not available yet. Emissions can also be included in map-based ECMS approaches, as suggested in [19]. An experimental validation of a method based on a constant weighting factor for NOx emissions has been provided in [20] by means of hardware-in-the-loop (HIL) experiments.

Other studies focus on the control of transient emissions, especially regarding hybrid powertrains relying on a diesel combustion engine. An optimal energy management strategy is provided by the authors of [21] with a DP approach for constant weighting factors related to NOx and particulate emissions. A key idea for the energy management strategy of a diesel HEV is to use the electric motor for torque phlegmatization during transients, for example by adopting heuristic methods, as presented in [22,23], or model-based frameworks, as reported in [24–26].

However, a control strategy that includes online adaptation of the weighting factors for the objective pollutant emission, to take into account real-world driving conditions and a possible modification of the emissions target level, has not been demonstrated yet. Therefore, in this paper, an energy management strategy that allows for the tracking of a specific NOx emission level to respect real driving emission constraints is presented. Under these conditions, the strategy is constructed to minimize the fuel consumption while sustaining the battery state of charge (SOC).

The paper is structured as follows: after explaining the vehicle model in detail in Section 2; the energy management strategy, which takes into account the real driving emissions, is derived in Section 3; then, Section 4 presents an experimental validation of the method proposed, as well as a simulation case study that quantifies the fuel savings potential compared to a standard method not accounting for real driving emissions.

2.Vehicle Model

The vehicle under investigation is a fictitious executive class sedan. The powertrain architecture is of the pre-transmission parallel type, as illustrated in Figure 1. The powertrain consists of a seven-step automatic gearbox, a 40 kW electric motor, the power electronics, a 2 kW h battery, a clutch and a 170 kW diesel engine. For the simulation of the vehicle behavior, only the longitudinal dynamics of the vehicle are of interest, since the consideration of these dynamics is sufficient for energetic considerations [27]. The longitudinal dynamics are simulated using a so-called forward approach, in which the physical causality is respected. By contrast, for the energy management of the powertrain, a model-based approach is chosen, in which the powertrain behavior is predicted by a so-called backward approach, where the physical causality is inverted [27].

More details on the modeling of the vehicle are given in the following. All the values needed to parameterize the model are listed in Table 1.

2.1.Longitudinal Dynamics

The equation for the longitudinal dynamics of the vehicle is described by [27]:
(1)dv(t)dt=1m(t)⋅(Tw(t)rw−12ρaircdAv2(t)−crmvagcos(γ(t))−mvagsin(γ(t)))with:
(2)m(t)=mv+Θwrw2+ΘEMig2(t)rw2where v denotes the vehicle speed; Tw denotes the wheel torque; rw denotes the wheel radius; ρair denotes the air density; cd denotes the aerodynamic drag coefficient; A denotes the frontal area of the vehicle; cr denotes the rolling friction coefficient; mv denotes the nominal vehicle weight; Θw and ΘEM denote the inertia of the four wheels and the electric motor, respectively; ag denotes the gravitational acceleration; and γ denotes the road slope.

2.2.Gearbox Model

The rotational speed of the gearbox input shaft, ωg, is given by:
(3)ωg(t)=v(t)rw⋅ig(t)where ig is the gear ratio of the gear engaged, g ∈ {1,…, 7}. The gear engaged, g, is decided by the gear decision variable, ug, defined by the energy management. For simplicity, the gearshifts are assumed to occur instantaneously without any delays.

The torque delivered to the wheels, Tw, is calculated by:
(4)Tw(t)=Tg(t)⋅ig(t)⋅(ηg,0−ηg,1ωg,1⋅ωg(t))sign(Tg(t))with Tg being the input torque of the gearbox; and ηg,0, ηg,1, ωg,1 being the parameters to model the speed-dependent gearbox losses, which account for the increased friction at higher gearbox input speeds [28]. The values for the efficiency parameters are chosen based on expert knowledge.

The input torque of the gearbox is:
(5)Tg(t)=Tm(t)+uc(t)⋅Te(t)where Tm is the motor torque; Te is the engine torque; and uc is the clutch state command determined by the energy management. The value uc = 1 means that the clutch is closed, and the value uc = 0 means that the clutch is open. The clutch command is assumed to be realized instantaneously.

2.3.Electric Motor Model

The desired torque for the electric motor, determined by the energy management, is assumed to be realized instantaneously without any delays. The power of the electric motor, including the power electronics, is given by a steady-state efficiency map, depicted in Figure 2, so that:
(6)Pm(t)=fm(ωm(t),Tm(t))

The minimum and maximum torque values for the electric motor, as well as the minimum and maximum speeds are defined by:
(7)0≤ωm(t)≤ωm,max(8)Tm,min(ωm(t))≤Tm(t)≤Tm,max(ωm(t))

2.4.Engine Model

The desired torque of the engine is determined by the energy management. However, the rate of change of the engine torque is limited by 100 N m/s in order to prevent the formation of excessive soot and NOx emissions during transients [23,31].

The rotational speed of the engine, ωe, is defined by (for brevity, the notation of time is omitted):
(9)ωe={0ifue=uc=0,ωe,idleif(ue=1anduc=0)or if(ue=uc=1andωg<ωe,idle)ωgif(ue=uc=1andωg≥ωe,idle)or if(ue=0anduc=1)where ue is the engine on/off command, with 0 standing for “off” and 1 standing for “on”; and uc is the clutch state command, with 0 standing for “open” and 1 standing for “closed”.

The mass flow rate of the fuel consumed,
m*fuel, and NOx emissions,
m*NOx, are calculated using measured steady-state maps in the form of:
(10)m*fuel(t)=fe,fuel(ωe(t),Te(t))(11)m*NOx(t)=fe,NOx(ωe(t),Te(t))

The map for the fuel efficiency and the map for the NOx emissions are shown in Figure 3.

At all times, the engine is only allowed to operate within its limits, defined by maximum torque and speed range conditions, i.e.,
(12)ωe,min(t)≤ωe(t)≤ωe,max(13)0≤Te(t)≤Te(ωe(t))

2.5.Battery Model

The battery is modeled as an equivalent circuit with a constant open-circuit voltage, Voc, in series with a constant internal resistance, Ri [27]. Since the battery is assumed to be of the LiFePO4 type [32,33], the assumption of a constant open-circuit voltage is a valid approximation for the typical operating range of the battery [34–36].

The equation of the dynamics of the battery SOC, ξ, is given by:
(14)dξ(t)dt=−Ib(t)Q0with:
(15)Ib(t)=Voc(ξ(t))−Voc2(ξ(t))−4Ri(Pm(t)+Paux)2⋅Ri(ξ(t))where Ib is the battery current, and Paux is a constant power consumed by electric auxiliary units. Although in this paper, Voc and R are assumed to be constant, the notation for the dependency on the SOC is preserved to keep the following derivation of the causal controller more general.

The battery current and the battery SOC are constrained to:
(16)Ib,min≤Ib(t)≤Ib,max(17)ξmin≤ξ(t)≤ξmaxrespectively.

2.6.Driver

The driver model consists of a proportional-integral (PI) controller. The output of the model is a preliminary throttle position, θ̃, calculated by:
(18)θ˜(t)=kp,D⋅(vref(t)−v(t))+1Ti,D⋅∫0t(vref(τ)−v(τ))dτwith kp,D and TI,D being the proportional and integral controller parameters with the manually tuned values of 0.9 and 11 s, respectively. The final throttle position, θ, is then obtained by saturating the preliminary throttle signal as follows:
(19)θ˜(t)={100ifθ˜(t)>100θ˜if−100≤θ˜(t)≤100−100ifθ˜(t)<−100

Due to the saturation and the integral control part of the driver model, an anti-windup scheme is implemented [37].

The mapping of the throttle position, θ, to a desired torque command, Treq, is designed, such that:

θ = 100% corresponds to the maximum possible traction torque provided by both the engine and the motor at the current vehicle speed;

θ = 0% corresponds to a zero torque request; and

θ = −100% corresponds to the maximum possible brake torque.

The intermediate torque commands are obtained by linear interpolation of the throttle position.

2.7.Energy Management

The energy management defines the set points for the gear number to be engaged, ug, the clutch open/closed state, uc, the engine on/off state, ue, and the torque split, uts, between the engine and the motor.

The gear command, ug, is determined by a lookup table of the form:
(20)ug(t)=fg(v(t),θ(t))as shown in Figure 4. The figure shows the upshift and downshift lines depending on the vehicle velocity, v, and the throttle position, θ. For comfort reasons, after every gear change, the new gear has to remain in its state for at least 3 s before another gear can be engaged. This time-based hysteresis can only be overruled in the case the throttle is fully depressed.

The clutch command, uc, is defined as follows: if the engine is off, the clutch is assumed to be open. If the engine is on and the desired engine torque is larger than zero, the clutch is closed. Otherwise, if the engine is on and the desired engine torque is equal to zero, the clutch is open, and the engine is assumed to be running idle.

The torque split, uts, is defined as the ratio of the motor torque and the requested torque from the driver:
(21)uts=TmTreq

Then, the torques of the motor and the engine are obtained by:
(22)Tm=Treq⋅uts(23)Te=Treq⋅(1−uts)

Recharging during standstill is not considered here; however, this feature could be implemented easily.

The commands for the engine on/off, ue, and the torque split, uts, are determined by an online optimization method commonly known as the ECMS [12,13,34]. To account for the diesel exhaust emissions, the standard approach has been extended, as shown in Section 3 below.

3.Controller

This section presents a detailed description of the controller developed to optimize the fuel consumption of a diesel hybrid vehicle under a constraint for the NOx emissions. The description of the control system is given in Section 3.1, while the solution of the control problem, derived by means of optimal control theory, is provided in Section 3.2. Then, Section 3.3 proposes a derivation of a causal online controller. To implement such a controller, the dependency between the co-states must be identified, which is illustrated in Section 3.6. The final structure of the controller is shown in Section 3.7, and the assignment of the corresponding reference signal trajectories is described in Section 3.8.

3.1.System Description

The hybrid vehicle system can be described by a system of first-order ordinary differential equations:
(24)ddtx(t)=f(x(t),u(t))where x is the state vector and u the control input vector. The state and the control input vector are:
(25)x(t)=[ξ(t)mNOx(t)],u(t)=[ue(t)uts(t)]where ξ is the battery SOC, mNOx the cumulated NOx emissions, ue the engine on/off command and uts the torque split command, as previously introduced in Section 2.7. The control input for the gear number, ug, is decided independently of x and u, and the control input for the clutch state, uc, is defined by a known function of ue. Therefore, ug and uc are not subject to optimization and, hence, not included in u.

Notice that the vehicle speed is not considered to be a state variable, because the speed is assumed to be perfectly tracked. Based on the actual vehicle speed and the driver's torque request, the required torques and the shaft speeds can be determined using a backwards calculation [27].

Using these definitions, the model for the system dynamics is derived based on the description presented in the previous section. The state dynamics equations become:
(26)ddtx(t)=f(x(t),u(t))=[ddtξ(x(t),u(t))ddtmNOx(x(t),u(t))]=[−Ib(ξ(t),u(t))Q0m*NOx(u(t))]

The choice for the integral of the emissions as a state variable has an advantage compared to the specific emission level, as defined in legislation as:
(27)m¯NOx=mNOxdwith d being the distance driven. The derivative of the latter expression would be:
(28)ddtm¯NOx=ddtmNOx⋅d−mNOx⋅vd2=1d⋅(ddtmNOx−m¯NOx⋅v)with v being the vehicle speed. The derivative of the absolute cumulated emissions, mNOx, is:
(29)ddtmNOx=m*NOxwhere
m*NOx denotes the NOx mass flow rate. The following properties hold:
(30)∂(ddtm¯NOx)∂m¯NOx=−vd≠0(31)∂(ddtmNOx)∂mNOx=0

These properties play a role in the following derivation.

3.2.Solution Using Pontragin's Minimum Principle (PMP)

According to the general methodology introduced in [38], the problem can be defined as an optimal control problem with partially constrained final states. The problem is formulated as:
(32)minu{J(u)=∫0Tm*fuel(u(t))dt}such that:
(33)x(0)=x0(34)ddtx(t)=f(x(t),u(t))for allt∈[0,T](35)x(T)∈S

The closed final-time set, S, describes the constraints for the final state vector: mNOx has an upper constraint, since it has to be below a certain emission level, while ξ has to satisfy the charge-sustaining condition for the battery.

The optimal solution is found using Pontryagin's minimum principle (PMP) [38,39]. By defining the Hamiltonian function as:
(36)H(x(t),u(t),λ(t))=m*fuel(u(t))−λT(t)⋅f(x(t),u(t))with λ denoting the co-state vector, the optimal solution (denoted by o) must satisfy the following conditions:
(37)x(0)=x0(38)ddtxo(t)=∇λH|o=f(x0(t),uo(t))(39)ddtλo(t)=−∇xH|o=−∇xm*fuel(uo(t))−[∂f∂x(xo(t),uo(t))]T⋅λo(t)

The co-state vector must stay within the normal cone, T*, of the target set, S, at the final time for xo(T):
(40)λo(T)∈T*(S,xo(T))

The Hamiltonian of the system must be minimized for all times with respect to all admissible inputs u, such that:
(41)H(xo(t),uo(t),λo(t))≤H(x(t),u(t),λ(t))

The dynamics of the co-states can be rewritten, separating the battery SOC and the cumulated NOx co-states, as follows (where the dependencies on the inputs and the states are omitted):
(42)[ddtλξo(t)ddtλNOxo(t)]=[−∂m*fuel∂ξ−[∂f∂ξ]T⋅λo(t)−∂m*fuel∂mNOx−[∂f∂mNOx]T⋅λo(t)]

Since the fuel mass flow and the NOx emissions do not explicitly depend on the battery SOC, the equations hold:
(43)∂m*fuel(u(t))∂ξ=0(44)∂f(x(t),u(t))∂ξ=[−∂∂ξIb(ξ(t),u(t))Q00]

Moreover, since neither the system dynamics nor the fuel mass flow depend explicitly on the cumulated NOx emissions, the following expressions hold:
(45)∂m*fuel(u(t))∂mNOx=0(46)∂f(x(t),u(t))∂mNOx=[00]

Combining Equations (43)–(46) with Equation (42) leads to the following rewritten expression for the co-states dynamics:
(47)[ddtλξo(t)ddtλNOxo(t)]=[∂∂ξIbξ(t),u(t)Q0⋅λξ(t)0]

The importance of the choice of cumulated emissions is here demonstrated, since the term,
∂f∂mNOx, would not have vanished with the choice of specific NOx emissions as a state variable Equation (30). Further, considering Equation (47) and that
λNOxo(T)≥0, we have that:
(48)λNOxo=const.=λNOx≥0

The Hamiltonian in Equation (36) can be rewritten, having defined λNOx:
(49)H(⋅)=m*fuel(u(t))−λT(t)⋅f(x(t),u(t))=m*fuel(u(t))+λNOx⋅m*NOx(u(t))−λξ(t)⋅Ib(ξ(t),u(t))Q0=:H˜(⋅)

The latter new expression for the Hamiltonian in Equation (49) leads to a new Hamiltonian, H̃, for the system for the minimization of a weighted sum of fuel and emissions. This leads to a new problem definition:
(50)minu{J˜(u)=∫0T(m*fuel(u(t))+λNOx⋅m*NOx(u(t)))dt}such that:
(51)ξ(0)=ξ0(52)ddtξ(t)=−Ib(ξ(t),u(t))Q0for allt∈[0,T](53)ξ(T)=ξ0

For the reformulated optimization problem, λNOx is not a co-state but a weighting factor and, therefore, a given parameter. It quantifies the fuel equivalent of a given amount of NOx emissions, in a new Lagrangian reformulation. For this reason, the authors of [40] have defined the Equations (50)–(53) equivalent emissions minimization strategy (EEMS). The equivalence of the minimization Equations (50)–(53) to the previous Equations (32)–(35) is guaranteed, since the eliminated state does not appear in any of the other equations, and it is introduced only to enforce the limit on cumulated emissions with a final-state constraint. As a consequence, if the constant equivalence factor is known, the same optimal results will be achieved for the redefined problem. As a matter of fact, such value is not known a priori. For this reason, a causal controller based on the online calculation of the equivalence factor is proposed in the following section.

3.3.Emissions and Charge-Sustaining Causal Control

The objective of this section is to present a feedback controller derived for online control of the cumulative NOx emissions. Before deriving the proposed causal control framework, the cost functional in Equation (50) is conveniently rearranged by introducing the normalized NOx mass flow rate,
m*˜NOx:
(54)m*˜NOx(t)=m*NOx(t)⋅m*NOx,MAXm*fuelMAXsuch that:
(55)J¯(u)=∫0T(m*fuel(u(t))+λNOx⋅m*˜NOx(u(t)))dtwhere the normalization takes place using the maximum flow rates of fuel and NOx for the given experimental maps of the combustion engine. By means of the rewritten cost functional in Equation (55), the argument of the object function and the corresponding Hamiltonian are expressed by homogenous measuring units.

To achieve the goal of generating a charge and emissions sustaining strategy, two terms can be added to Equation (55) in order to penalize deviations from reference values for SOC, ξref, and deviations from the reference values for the normalized cumulated emission, m̃NOx,ref, thus obtaining a new formulation for the cost functional:
(56)J^(u)=∫0TL^(x(t),u(t))dt(57)=∫0T(m*fuel(u(t))+μ⋅(m˜NOx,ref(t)−m˜NOx(t)Δm˜NOx,norm)2q+β⋅(ξref(t)−ξ(t)Δξnorm)2p)dt

The extended cost Equation (57) leads to the extended Hamiltonian:
(58)H^(u)=m*fuel(u)+μ⋅(m˜NOx,ref(t)−m˜NOx(t)Δm˜NOx,norm)2q+β⋅(ξref(t)−ξ(t)Δξnorm)2p+λT⋅f(x,u)which must be minimized by the optimal input sequence:
(59)uo=argmin(H^(u))

Since the additional terms of the extended Hamiltonian do not explicitly depend on the control inputs, they will be minimized by the optimal policy, uo, as well.

To calculate the co-state dynamics, λ, the Hamilton–Jacobi–Bellman equations provide the following expression [39] for the optimal co-states vector, λo:
(60)λo(x,t)=∂C^o(x,t)∂xwhere Ĉo(x,t) denotes the optimal cost-to-go function associated with the cost function, ĵ(u). The optimal cost-to-go function for the state, x, at the time, t, is given by:
(61)C^o(x,t)=minu{∫tTL^(x(τ),u(τ))dτ|x(t)=x}

Since the optimal cost-to-go is not known a priori in a causal setting, the optimal cost-to-go function is estimated by a sub-optimal time-invariant function, formed by the sum of different independent cost indices, as follows [41]:
(62)C^(x)=C^(ξ,m˜NOx)==C^f1,NOx(m˜NOx)+C^f1,ξ(ξ,m˜NOx)+C˜f2+C˜ξ(ξ)+C˜m˜NOx(m˜NOx)

The five terms are explained in the following:

Ĉf1,ξ(ξ,m̃NOx): the additional fuel consumption caused by compensating for the current SOC deviation;

Ĉf1,NOx(m̃NOx): the additional fuel consumption caused by bringing the cumulated NOx close to the reference level;

C̃f2: a fuel consumption that is supposed to be independent of both the current SOC and the current emission level, needed to drive the rest of the driving mission with correct reference values;

C̃ξ(ξ): denotes the penalty for SOC deviations from the reference value;

C̃m(̃NOx)(m̃NOx): denotes the penalty of m̃NOx deviations from the reference value.

3.3.1.Cost of Sustaining the Battery SOC Ĉf1,ξ(ξ, m̃NOx)

Following the approach described in [41,42], the fuel energy used to compensate for the SOC deviations from the target value can be approximated by first estimating the energy stored in the battery at a certain ξ(t) with respect to the reference, ξref(t), as:
(63)EBT,Δξ=Q0∫ξ(t)ξref(t)Voc(ξ˜)dξ˜

Since in the future, this energy must be compensated for using the thermal path, a certain amount of fuel will be saved/consumed to discharge/charge the battery. Such a quantity will clearly depend on the future efficiencies of the engine and the electric path, which, in turn, depend on the future engine operating points. Moreover, the used engine operating points will also depend on the cumulated NOx emissions, since a second controller acts in parallel, modifying the choice of control inputs to track the desired emissions level. As a consequence, the average future charging/discharging overall efficiency, ηc, is a function of m̃NOx. The resulting cost to sustain the battery charge is:
(64)C^f1,ξ(ξ,m˜NOx)≈EBT,Δξηc(m˜NOx)⋅Hl=Q0ηc(m˜NOx)⋅Hl⋅∫ξ(t)ξref(t)Voc(ξ˜)dξ˜where Hl represents the lower heating value of the fuel.

3.3.2.Cost of Saving NOx Emissions Ĉf1,NOx(m̃NOx)

An expression for the additional fuel cost of saving NOx emissions can be found under the hypothesis that the battery SOC, ξ, has a much smaller time constant than the time horizon considered to bring the cumulated emissions to the reference value (e.g., minutes). If this holds, and the parameter, KFCN, identifies the relationship between fuel consumption and NOx emissions for the given engine, the associated cost-to-go can be expressed as:
(65)C^f1,NOx(m˜NOx)=Δmfuel=KFCN⋅Δm˜NOx=KFCN⋅(m˜NOx−m˜NOx,ref)

3.3.3.Cost of the Emissions Deviation Penalty C̃m̃NOx (m̃NOx)

The second term in Equation (57) stands for the costs for a future deviation of the actual NOx mass from its reference value. Under the hypothesis that the NOx controller is able to diminish the error between the actual and the reference value linearly with time, the cost of the emission penalty in the future can be estimated. The future evolution of the error between the actual and the reference NOx emissions at time τ ∈ [0, Th] is therefore calculated as:
(66)m˜NOx,ref−fut(τ)=m˜NOx−fut(τ)=(m˜NOx,ref(t)−m˜NOx(t))⋅(1−τTh)

The cost for the penalty is obtained by integrating this trajectory as follows:
(67)C˜m˜NOx(m˜NOx)=∫0Thμ⋅(m˜NOx,ref−fut(τ)−m˜NOx−fut(τ)Δm˜NOx,norm)2qdτ=μ⋅Th2q+1⋅(m˜NOx,ref(t)−m˜NOx(t)Δm˜NOx,norm)2q

3.3.4.Cost of the SOC Deviation Penalty C̃ξ(ξ)

Similarly to the treatise in Section 3.3.3, the costs of the third term in Equation (57) can be expressed by an equation equivalent to Equation (67). Here, the error between the actual SOC and the reference SOC is assumed to linearly diminish within the time, Tk. This leads to the following cost for the SOC deviations:
(68)C˜ξ(ξ)=∫0Tkβ⋅(ξref−fut(τ)−ξfut(τ)Δξnorm)2pdτ=β⋅Tk2p+1⋅(ξref(t)−ξ(t)Δξnorm)2p

3.3.5.Total Cost and Equivalence Factors

The total sub-optimal cost-to-go Equation (62) can now be expressed using Equations (67) and (68) as:
(69)C^(x)=C^(ξ,m˜NOx)==Q0ηc(m˜NOx)⋅Hl⋅∫ξ(t)ξref(t)Voc(ξ˜)dξ˜+KFCN⋅(m˜NOx(t)−m˜NOx,ref(t))+C^f2+β⋅Tk2p+1⋅(ξref(t)−ξ(t)Δξnorm)2p+μ⋅Th2q+1⋅(m˜NOx,ref(t)−m˜NOx(t)Δm˜NOx,norm)2q

Accordingly, the sub-optimal cost-to-go function is time-invariant, and so will be the sub-optimal co-states, which are given by:
(70)λ(x)=[λξ(x)λNOx(x)]=∂C^(x)∂x=[∂∂ξC^(x)∂∂m˜NOxC^(x)]

The partial derivatives of Equation (69) generate the following expressions for the two co-state Equation (70):
(71)λξ(x)=∂C^(x)∂x=−Q0⋅Voc(ξ)ηc(m˜NOx)⋅Hl−β˜⋅(ξref(t)−ξ(t))2p−1(72)λNOx(x)=∂C^(x)∂m˜NOx=KFCN+Q0Hl⋅∫ξ(t)ξref(t)Voc(ξ˜)dξ˜⋅∂∂m˜NOx(1ηc(m˜NOx))+μ˜⋅(m˜NOx,ref(t)−m˜NOx(t))2q−1where the following substitutions are adopted:
(73)μ˜=−2qμTk(2q+1)⋅Δm˜NOx,norm2p(74)β˜=2pβTh(2p+1)⋅Δξnorm2p

By analyzing Equations (71) and (72), the mutual relationship between the co-states is evident. In more detail, the first two terms of Equation (72) are formed by a theoretically constant term, KFCN, that instead depends on the operating points occurring in the period considered and another term that we suppose to be negligible, under the hypothesis that the dynamics of the average charging/discharging efficiency does not depend directly on the cumulated emissions. The simplified approach followed in this section is to replace the constant equivalence terms by a factor, λNOx,0, and to add an integrator, which is used to online adapt it during operation, to respect the average emissions target:
(75)λNOx(m˜NOx)=λNOx,0+μ˜⋅(m˜NOx,ref(t)−m˜NOx(t))2q−1+∫0tm˜NOx,ref(τ)−m˜NOx(τ)Ti,NOxdτ

Equation (75) represents a PI controller for the cumulative emissions level, which in online applications, will be measured by means of a dedicated sensor. The term, λNOx, can be directly implemented in the cost functional in Equation (55) that solves the equivalent Equations (50)–(53). Its corresponding Hamiltonian will be expressed by:
(76)H¯(x,u,λ)=[m*fuel(u)+λNOx(m˜NOx)⋅m*˜NOx(u)−λξ(ξ,m˜NOx)⋅Ib(ξ,u)Q0]

Equation (76) can be be expressed in terms of power by multiplying the whole Hamiltonian with the lower heating value, Hl, of the fuel to yield:
(77)H_(x,u,λ)=[Pfuel(u)+λNOx(m˜NOx)⋅PNOx(u)−Hl⋅λξ(ξ,m˜NOx)⋅Ib(ξ,u)Q0]=[Pfuel(u)+λNOx(m˜NOx)⋅P˜NOx(u)−s*(ξ,m˜NOx)⋅Pb(ξ,u)]with Pb representing the inner electrochemical battery power and the electrical energy equivalence factor [34]:
(78)s*(ξ,m˜NOx)=−λξ(ξ,m˜NOx)⋅HlVoc(ξ)⋅Q0

The combination of Equations (71) and (78) leads to a reformulated expression for the electrical energy equivalence factor, s*:
(79)s*(ξ,m˜NOx)=1ηc(m˜NOx)+HlVoc(ξ)⋅Q0⋅β˜⋅(ξref(t)−ξ(t))2p−1=s0*(m˜NOx,ξ)+β˜*⋅(ξref(t)−ξ(t))2p−1

Since the average conversion efficiency, ηc, will vary depending on the operating points of the components involved (engine, electric motor, battery) and the operating points will vary as a function of the driving cycle and of the NOx feedback controller already introduced,
s0* has to be adjusted during operation. The first adjustment depends directly on the actual normalized cumulative emission, m̃NOx, due to the action of the controller that online adapts λNOx, which will be clarified in the next section. The second adaptation is achieved using an integrator with integration time Ti,ξ as follows:
(80)s*(ξ,m˜NOx)=s0*(λNOx(m˜NOx))+β˜*⋅(ξref(t)−ξ(t))2p−1+∫0tξref(τ)−ξ(τ)Ti,ξdτ

Such a PI controller for the electrical energy equivalence factor was also proposed by the authors of [42,43].

Equations (75) and (80) represent the structure of the desired online causal emission and the battery charge controller. It is formed by two feedback PI controllers, linked together by a relationship between the constant values of the co-states,
s0*(λNOx(m˜NOx)). Since both PI controller outputs can be saturated, the controllers are extended with an anti-windup scheme [37,44].

The final structure of the controller is presented in Section 3.7. Since the controller is able to control the real driving NOx emissions (RDE) and since it is based on the ECMS, the controller is referred to as the RDE-ECMS.

3.4.Normalization of the Emissions Co-State

By normalizing the Hamiltonian function in Equation (77) by dividing H(·) by (λNOx + 1), the co-state, λNOx, is turned into a weighting factor, αNOx, to yield:
(81)H∘(x,u,λ)=1(λNOx+1)⋅[Pfuel(u)+λNOx⋅P˜NOx(u)+s0*(λNOx)⋅Pb(ξ,u)]=1(λNOx+1)⋅Pfuel(u)+λNOxλNOx+1⋅P˜NOx(u)+s0*(λNOx)λNOx+1⋅Pb(ξ,u)=αNOx⋅Pfuel(u)+(1−αNOx)⋅P˜NOx(u)+s0(αNOx)⋅Pb(ξ,u)

The introduction of the reformulated weighting factor, αNOx, for NOx emissions leads to the following possible cases:
(82)αNOx(t)={0NOxoptimal0<αNOx(t)<1NOx−fuel trade off1fuel optimal

The proportional gains of the PI controllers then become:
(83)kp.ξ=β˜*λNOx+1(84)kp,NOx=μ˜λNOx+1

3.5.Preventing Frequent Engine Starts and Stops

To prevent excessively frequent engine starts and stops that can arise due to the application of optimal control-based methods [12,13], a penalty, δ, for changing the engine on/off state is introduced [45,46]:
(85)H(x,u,λ)=αNOx⋅Pfuel(u)+(1−αNOx)⋅P˜NOx(u)+s0(αNOx)⋅Pb(ξ,u)+δ⋅Ie(x,u)

The function, ℐe, denotes an indicator function detecting a change request for the engine on/off state. If a change is requested, the indicator function is one and zero otherwise. By manual tuning, a value for the penalty, δ, of 0.1 × 10−3 kg × 43 MJ/kg proved to yield a reasonable performance. However, a penalty alone cannot ensure a minimum engine on/off dwell time, which is desired for comfort and emissions considerations.

Therefore, an additional heuristic engine on/off comfort function is implemented similarly to the one presented in [47]. In this comfort function, the desired engine on/off change request signal from the extended ECMS is not realized instantaneously. Instead, it has to remain in the same state, either on or off, for at least 1 s until it is transferred to the next level of a series of checks. At the next level, an engine on/off change request is only realized if the previously realized engine on/off state has remained for at least 5 s in its state. As such, the engine is either on or off for at least 5 s. This hysteresis can only be overruled if the throttle is fully depressed.

Due to these measures, the average number of engine starts and stops, on the four here considered driving cycles, could on average be reduced to a reasonable amount of 2.1 starts per minute compared to 4.6 starts per minute without any measure to prevent frequent starts and stops. The loss in fuel economy due to this comfort function amounts on average to 3.8% compared to a theoretical value for the fuel consumption obtained without any comfort function. The minimum engine on/off dwell time amounts to 5 s in almost any case for the driving scenarios considered.

3.6.Equivalence Factors Dependency Identification

The goal of this section is to describe the methodology applied to identify the dependency of s0 on αNOx. The key idea is to apply an optimal control method to the optimization problem described by the Hamiltonian function in Equation (85), such as, for example, DP or PMP with constant values for the weighting factors, αNOx and s0. In this case, PMP is adopted, since it is more suitable for the present application of a forward-facing vehicle model, including many input and state variables. This procedure is applied for several different values of αNOx, in order to identify, for each value, the corresponding unique constant equivalence factor, s0, that ensures a charge sustaining condition.

The methodology is applied to various driving scenarios, in this case to the four well-known driving cycles, New European Driving Cycle (NEDC), Federal Test Procedure 75 (FTP-75), Worldwide Harmonized Light Vehicles Test Procedure (WLTP) and California Unified Cycle (LA92). The simulation results of the identification procedure are depicted in Figure 5a. The relationship between αNOx and s0 is almost linear. However, a linear fit can lead to s0-values, which yield considerable deviations of the final SOC when simulating the vehicle on certain driving cycles without feedback. A quadratic fit, as indicated by the black dashed line, turned out to be more adequate for generating a unique relationship between the equivalence factors for the driving cycles considered.

Figure 5b shows the trade-off between fuel consumption and NOx emissions, using charge-sustaining s0-values for each αNOx ∈ {0,0.25,0.5,0.75,1}. For each driving cycle, the corresponding curve represents the optimal trade-off for the approach presented using Equation (85). Any causal method based on Equation (85) cannot yield results that are to the left or below the corresponding trade-off for a specific driving cycle.

3.7.Controller Structure

Based on the mathematical derivation of the controller presented in the previous sections, the desired controller, to be tested in simulation and experimental tests in the following sections, is illustrated in Figure 6.

3.8.Calculation of SOC and mNOx Reference Trajectories

The two PI controllers of Figure 6 require a proper definition of the reference trajectories for the respective controlled variables. The reference value for the SOC, ξref(t), could simply be a constant value representing the desired final SOC, which also coincides with the initial value, ξ0, to enforce a charge-sustaining constraint. Alternatively, the reference value can take into account that the current kinetic and potential energy of the vehicle can be recuperated in the future and stored as electrical energy with certain efficiencies, ηc,K, ηc,P, resulting in the following expression [48,49]:
(86)ξref(t)=ξ0−(ηc,K⋅12m(t)⋅v2(t)Q0Voc(ξ0)+ηc,P⋅m(t)⋅g⋅(h(t)−h0)Q0Voc(ξ0))where h(t) − h0 is the height difference between the current altitude and some reference altitude. The numerical values for ηc,K and ηc,P are both 0.65, which were obtained by manual tuning.

The reference cumulative emissions can be computed in a simple way by defining a specific emission level, m̄NOx, Equation (27) expressed in units of mg/km:
(87)mNOx,ref(t)=∫0tm¯NOx⋅v(τ)

The value for m̄NOx can reflect a potential limit imposed by legislation. However, since the controller is purely causal, the controller cannot guarantee to keep the NOx emissions below the limit. Therefore, a value for mNOx should be chosen that is below the actual limit value issued, for example, by legislation.

Further, Equation (86) can easily be extended to account for the characteristics of plug-in HEVs. In this case, the constant, ξ0, can be replaced by a driving distance-dependent reference trajectory, as proposed by the authors of [50]. Another possibility for both charge-sustaining and plug-in HEVs is to use predictive information to calculate a reference trajectory for the SOC [42,51], as well as for the NOx emissions. Even further, predictive methods could be applied to directly optimize the Equations (32)–(35). However, the application of these ideas to the problem at hand is left for future research.

4.Results

First, an experimental validation of the models and of the methodology, as presented in Sections 2 and 3, is shown. Then, a case study is presented in which the benefit of using an RDE-ECMS compared to an ECMS with a constant emission-related equivalence factor, αNOx, is analyzed.

4.1.Experimental Validation

The goal of this subsection is first to show that the RDE-ECMS presented in Section 3 works also in practice, and second, that the quasi-static modeling for the fuel consumption and the NOx emissions is sufficient.

4.1.1.Experimental Setup

For the experimental validation of the RDE-ECMS, the method presented in Section 3 was applied both in simulation and in HIL experiments. In the HIL experiments, only the engine was used in real hardware. The longitudinal dynamics, as well as the vehicle components were simulated on a computer. This setup allowed for the measurement of the real fuel consumption and the real NOx emissions without requiring the physical presence of the entire vehicle. For more details on HIL experiments, interested readers are referred to the literature [52–54].

In the HIL experiments here, the desired torque command, which is calculated by the energy management controller, is sent to the electronic control unit (ECU) of the engine, while the desired engine speed command is sent to the dynamometer of the engine test bench. The NOx emissions are measured using a VDO/NGK UniNOx sensor manufactured by Continental AG, Hanover, Germany. This sensor is likely to be employed also in real vehicles, due its low price and due to its ability to additionally measure the air-to-fuel ratio. The fuel consumption is taken by the ECU-internal estimation.

4.1.2.Experimental Results on the WLTP

We used the following setup to compare the simulation results to the results obtained with the HIL experiments:

Figure 7 shows the vehicle speed, the battery SOC, the normalized specific NOx emissions, the normalized specific fuel consumption and the equivalence factors, αNOx and s, for both the NOx emission and the battery power.

4.1.3.Comparison of the Simulation Results to the HIL Experiments

As can be seen from the figure, the vehicle speed trajectories of the simulation results and the HIL results are identical. The SOC trajectories of the simulation and the HIL results are similar; both are charge-sustaining at around the same SOC reference level. Furthermore, the NOx trajectories of both are very similar; after an initial transient phase, the trajectories become more stable, and they approach the desired reference level. A similar behavior is observed for the specific fuel consumption that, in addition, exhibits a visible offset that is explained below. Further, the dynamics of the equivalence factor, αNOx, are also similar for both the simulation and the HIL experiment. They both start from a value of zero and converge towards a value of one until the specific NOx emissions approach the reference level. Due to the highway part, which requires the engine to be used a lot, the equivalence factor, αNOx, is reduced to save NOx emissions. In the subsequent city driving part, αNOx again approaches a value of one until the next highway part, where the value of αNOx is reduced once more. Furthermore, the trajectories for the equivalence factor, s, are close to each other. The preliminary conclusion from the comparison of the simulation results to the experimental results is that in practice, the RDE-ECMS works the same as in the simulation.

However, as can be seen from the figure, the RDE-ECMS cannot guarantee maintaining the NOx emissions below the reference level since: (i) the strategy is purely causal; and (ii) the reference level can be too low for the controller to have a sufficient influence on the control action. Therefore, the reference trajectory has to be designed carefully.

The offset of the trajectories of the simulation and the HIL results is a consequence of some of the neglected dynamics in the engine model of the simulation. For example, the thermal dynamics are not considered in the simulation, although in practice, they can have an influence on the formation of pollutant emissions. In fact, such effects influence the behavior of the SOC and the NOx control, such that the trajectories of the SOC, NOx emissions, etc., become different for the results obtained with the simulation and the experiment. However, the offset between the trajectories is not a direct measure to quantify the modeling errors, because the trajectories show the actual, “uncorrected” emissions and the actual, “uncorrected” fuel consumption. For example, an offset of 5% in the trajectories for the fuel consumption does not mean that the true fuel consumption is 5% different, since there is also a certain offset in the SOC trajectories. Typically, a higher SOC means a also higher fuel consumption and also higher NOx emissions. For a fair comparison, the equivalent NOx emissions and the equivalent fuel consumption have to be calculated.

To make a fair comparison between the simulation and the experimental results, the NOx emissions and the fuel consumption have to be corrected to take into account the different levels of the SOC. Here, the correction is made based on the following equations:

Δmfuel,eq and ΔmNOx,eq are the equivalent fuel mass and the equivalent NOx mass as a function of the final battery charge, respectively [34]:
(90)Δmfuel,eq={ΔEbHl⋅1η¯b(d)⋅η¯m(d)⋅η¯eifΔEb≥0ΔEbHl⋅η¯b(c)⋅η¯m(c)⋅η¯eifΔEb<0(91)ΔmNOx,eq=∫0Tm*NOx(t)dt∫0Tm*fuel(t)dt⋅Δmf,eqwith η̄b, η̄m, η̄e denoting the averaged efficiencies of the battery, the motor and the engine, while the superscript, (c) denotes “charging phase” and (d) denotes “discharging phase”.

ΔEb is the amount of net energy stored in the battery at the end of the driving cycle:
(92)ΔEb=Q0∫ξ(0)ξ(T)Voc(ξ˜)dξ˜

Figure 8 shows a comparison of the SOC-corrected equivalent NOx emissions and corrected specific fuel consumption on the two driving cycles, NEDC and WLTP, for each, considering the two cases of a constant αNOx and a variable αNOx. The case with the constant αNOx-value refers to a simulation without feedback control of the NOx emissions. Instead, a constant αNOx = 0 is used throughout the simulation. By contrast, the case with a variable αNOx represents the case of the RDE-ECMS.

According to the figure, the results obtained with the simulation underestimate the results obtained with the HIL experiments. The error in all the simulations compared to the experimental data is well below 5%. Therefore, quasi-static models for the fuel consumption and the NOx emissions can be considered to be accurate enough for the simulation on the NEDC and WLTP driving cycles.

4.2.Simulation Case Study

By now, it was shown that the RDE-ECMS can simultaneously control the SOC and the NOx emissions. To show in addition that the RDE-ECMS yields a lower fuel consumption than a non-adaptive ECMS, these two strategies are compared on four different driving cycles, namely the NEDC, the FTP-75, the WLTP and the LA92.

Assume that the RDE have to be lower than a specific value, say 1.02 or 102% for the specific NOx emissions in this case. Two causal strategies are considered, which respect this NOx limit: (1) a non-adaptive PI-controlled ECMS, which is the RDE-ECMS, but with a fixed value for αNOx; and (2) the RDE-ECMS presented in Section 3 with a specific NOx reference value of 1.0.

The non-adaptive ECMS was tuned to respect the NOx emission limit on the worst case driving cycle, which, here, is LA92, while giving the lowest possible fuel consumption on all driving cycles. The optimized parameters of the non-adaptive PI-controlled ECMS are kp,ξ = 2, Ti,ξ = 480, p = 1 and αNOx = 0. For comparison, the parameters of the RDE-ECMS are kp,ξ = 2, Ti,ξ = 480, kp,NOx = 0.1, Ti,NOx = 75,000, q = 1 and αNOx(0) = 0. Note that, for simplicity, the values for the parameters, kp,ξ and Ti,ξ, of the RDE-ECMS were taken from the non-adaptive ECMS. The two other parameters, kp,NOx and Ti,NOx, were optimized on the NEDC and the WLTP driving cycles to yield an acceptable reference tracking, namely neither a too fast nor a too slow tracking. To ensure a fair comparison, the initial value for αNOx was chosen to be equal to that of the non-adaptive ECMS.

These two strategies were applied on each of five repetitions of the four different driving cycles, NEDC, FTP-75, WLTP and LA92. Figure 9 shows the normalized specific equivalent fuel consumption, the relative fuel consumption difference of the RDE-ECMS compared to the non-adaptive ECMS and the normalized specific equivalent NOx emissions for both the non-adaptive ECMS (“αNOx = 0”) and the RDE-ECMS (“αNOx var.”, where “var.” stands for “variable” in the meaning of “not fixed”).

As seen from Figure 9b, the fuel savings of the RDE-ECMS amount to 0%–7% compared to the non-adaptive ECMS. On all driving cycles, both strategies respect the NOx emission limit, indicated as a red line in the plot on the right-hand side. In the case of the LA92 driving cycle, the RDE-ECMS behaves the same as the non-adaptive ECMS, which is the reason that there are practically no fuel savings. On all other driving cycles, the RDE-ECMS provides a lower fuel consumption than the non-adaptive ECMS, since the RDE-ECMS increases the emission-related equivalence factor, αNOx, to move along the optimal NOx-fuel trade-off.

Figure 10 shows the performance of the RDE-ECMS and the non-adaptive ECMS in relation to the optimal trade-off between the fuel consumption and NOx emissions already shown in Figure 5b. The non-adaptive ECMS (“diamond marker”) achieves practically identical performance as the optimal non-causal solution for αNOx = 0 in terms of fuel consumption and NOx emissions. The RDE-ECMS (“star marker”) achieves a performance close to the optimal trade-off curve for the driving cycles, FTP-75, WLTP and LA92. For the NEDC, there is still a potential to reduce the fuel consumption by about 1.6% for the same amount of NOx emissions as calculated with the non-causal ECMS.

Overall, the RDE-ECMS proved to minimize the fuel consumption, while tracking a reference NOx emission level and sustaining the battery charge.

5.Conclusions

This paper presents an energy management strategy to account for NOx RDE of a diesel HEV. The method is based on the ECMS, which is extended with a state accounting for the NOx emissions. As demonstrated in simulation, as well as in HIL experiments, the strategy is able to minimize the fuel consumption, while following given reference trajectories for the NOx emissions and the battery SOC. By simulation, the strategy is shown to optimally adjust the trade-off between the fuel consumption and NOx emissions during operation. Compared to a conservative non-adaptive strategy, the advantages in terms of fuel consumption amount to more than 7% in favor of the proposed method. The strategy can be employed also in plug-in HEVs, without the need of adjusting the controller structure, by only modifying the reference trajectory for the battery SOC.

So far, the presented RDE-ECMS has been applied to warm engine conditions. A future evolution of this strategy will be to investigate the control of tailpipe NOx emissions for a diesel HEV equipped with a SCR system. This will require a model of the reduction efficiency of the aftertreatment system as a function of an additional state, represented by the thermal state of the SCR system. Moreover, a further potential is seen in the optimal design of the reference trajectories for the SOC and the NOx emissions.

A hybrid electric vehicle (HEV) architecture considered in this paper: pre-transmission parallel HEV.

Figure 2.

Efficiency map (in %) of the electric motor, including the inverter losses. The underlying measurement data were obtained for a 25-kW electric motor. The torque axis was then scaled linearly with the nominal power [29,30].

Figure 3.

Experimental data for the fuel efficiency and the NOx emissions of the diesel engine. The torque is artificially limited at 300 N m, due to the test bench limitations, while the actual maximum torque is 540 N m. The data covers a speed range of 1000–3000 rpm; higher speeds are not considered. (a) Measured normalized fuel efficiency of the diesel engine; and (b) measured normalized NOx emissions of the diesel engine (
g/kW hg/kW h).

Figure 4.

Heuristic gearshift table.

Figure 5.

Relationship between the equivalence factors and the trade-off between the fuel consumption and NOx emissions simulated for the New European Driving Cycle (NEDC), Federal Test Procedure 75 (FTP-75), and California Unified Cycle (LA92) driving cycles: (a) relationship between the equivalence factors, αNOx and s0; and (b) trade-off between the fuel consumption and NOx emissions as a function of αNOx.

Figure 6.

Controller structure.

Figure 7.

Comparison of results obtained with the simulation and with the hardware-in-the-loop (HIL) experiment on three repetitions of the Worldwide Harmonized Light Vehicles Test Procedure (WLTP).

Figure 8.

Comparison of the simulation results to the hardware-in-the-loop (HIL) experiments. The equivalent NOx emissions (m̄NOx,eq) and the equivalent fuel consumption (m̄fuel,eq) take into account the correction based on the SOC deviation between the initial and the final SOC. (a) NEDC, αNOx = 0; (b) NEDC, αNOx = var.; (c) WLTP, αNOx = 0; and (d) WLTP, αNOx = var.