**Optimal Calibration Strategy of a Hybrid Electric Vehicle Equipped with an Ultra-Lean Pre-Chamber SI Engine for the Minimization of CO2 and Pollutant Emissions**

#### **Fabio Bozza, Vincenzo De Bellis, Enrica Malfi \*, Luigi Teodosio and Daniela Tufano**

Dipartimento di Ingegneria Industriale, Università di Napoli Federico II, 80125 Napoli, Italy; fabio.bozza@unina.it (F.B.); vincenzo.debellis@unina.it (V.D.B.); luigi.teodosio@unina.it (L.T.); daniela.tufano@unina.it (D.T.)

**\*** Correspondence: enrica.malfi@unina.it; Tel.: +39-081-7683264

Received: 15 June 2020; Accepted: 28 July 2020; Published: 3 August 2020

**Abstract:** The complexity of modern hybrid powertrains poses new challenges for the optimal control concerning, on one hand, the thermal engine to maximize its efficiency, and, on the other hand, the vehicle to minimize the noxious emissions and CO2. In this context, the engine calibration has to be conducted by considering simultaneously the powertrain management, the vehicle characteristics, and the driving mission. In this work, a calibration methodology for a two-stage boosted ultra-lean pre-chamber spark ignition (SI) engine is proposed, aiming at minimizing its CO2 and pollutant emissions. The engine features a flexible variable valve timing (VVT) control of the valves and an E-compressor, coupled in series to a turbocharger, to guarantee an adequate boost level needed for ultra-lean operation. The engine is simulated in a refined 1D model. A simplified methodology, based on a network of proportional integral derivative (PID) controllers, is presented for the calibration over the whole operating domain. Two calibration variants are proposed and compared, characterized by different fuel and electric consumptions: the first one aims to exclusively maximize the brake thermal efficiency, and the second one additionally considers the electric energy absorbed by the E-compressor and drained from the battery. After a verification against the outcomes of an automatic optimizer, the calibration strategies are assessed based on pollutant and CO2 emissions along representative driving cycles by vehicle simulations. The results highlight slightly lower CO2 emissions with the calibration approach that minimizes the E-compressor consumption, thus revealing the importance of considering the engine calibration phase, the powertrain management, the vehicle characteristics, and its mission.

**Keywords:** virtual engine calibration; ultra-lean combustion; hybrid vehicle; 0D-1D engine modelling

#### **1. Introduction**

In the 2013, the Diesel-gate scandal erupted in the USA, flooding into Europe and worldwide, involving more than 11 million of vehicles produced by various car manufacturers [1]. The related rumors hit the entire automotive sector, leading to the idea that the traditional internal combustion engine (ICE), fed by fossil fuels, would disappear within a few years. A full-electric mobility, based on battery electric vehicles (BEVs), has been deeply promoted by politicians and public opinion, as the unique available solution to fulfill more and more stringent regulations [2].

At the same time, however, alternative propulsion architectures, such as hybrid electric vehicles (HEVs) and plug-in HEV (PHEVs), due to their lower costs and greater autonomy, continuously gained new market shares. Today, they represent much more than a temporary solution, filling the gap between standard ICE-powered vehicles and full electric mobility.

As known, the choice for the optimal electrification level of a vehicle requires to consider its entire life cycle, including the CO2 formation during the vehicle manufacturing and the one associated to electrical energy production and distribution. Under this point of view, the overall CO2 emissions from the two extreme antagonist vehicles (ICEs and BEVs) become comparable and are strongly dependent on the so called "energy-mix" utilized for electrical energy production [3,4]. In particular, BEV or PHEV have consistent benefits only when renewable energy sources are used, leading to the conclusion that an electric mobility cannot be considered as the only practical alternative due to the limited and country-dependent availability of renewable sources. Additionally, as already remarked, the higher ownership cost of the BEV, combined with the high battery charging time and the low autonomy, strongly limits their market acceptance and diffusion [5]. In the light of those concerns, rather than a pure-electric future, a scenario should be expected in the years to come characterized by variegated technologies that are best suited to the contest in which they are employed. This means that ICE-based vehicles, HEVs, PHEVs, BEVs and fuel cell-based vehicles will coexist in the market for a long time, pushing car manufactures to improve the propulsion system efficiency [6] with the aim to respect the stringent regulation in terms of CO2 and pollutant emissions.

Among the various pathways to improve the thermal efficiency of ICEs, various knock mitigation measures were attempted, such as the variable compression ratio [7], the cooled exhaust gas recirculation (EGR) [8], and the water injection [9]. Benefits of fuel consumption highly depend on the considered operating condition [10]. An increasing research work is nowadays devoted to ultra-lean combustion systems since they have the potential of simultaneously reducing the cylinder-out NOx emissions and the fuel consumption in the whole engine operating domain. Among the various concepts analyzed in the literature to sustain a lean combustion [11,12], such as fuel stratification [13,14], HCCI [15] and SACI [16], the employment of a pre-chamber (PC) ignition system is nowadays considered the most practical solution. In this architecture, a small PC volume is connected to the main chamber (MC) through some orifices. The combustion process starts in the PC and propagates in the MC one, in the form of hot burned jets ejected from the pre-chamber. The high turbulence intensity of these jets ignites the MC mixture and ensures a stable flame propagation even under extremely lean mixtures [17,18].

The aim of this activity is to numerically analyze the calibration and the energy management strategy of a HEV equipped with an ultra-lean PC engine. The overall propulsion system under investigation is developed within the framework of the EAGLE H2020 project (https://h2020-eagle.eu/). The hybrid architecture and the components' sizing has been completely defined by the manufacturer and includes two electric motor/generator units, a battery, and two gearboxes. The system is designed to provide the maximum flexibility and a versatile series/parallel operation. A hybrid series-parallel architecture was chosen to combine the advantages of series powertrain in urban driving (including smooth takeoff) and of the parallel mode over rural/highway routes. Close-to-optimal operations for the ICE are guaranteed by a four-speed gearbox placed between the ICE and differential. A gearbox is preferred to a CVT power-split device to get a higher mechanical efficiency. The switch between series and parallel modes is handled by three clutches, which also help to decrease the mechanical losses when one of the motors is not used.

The thermal unit is equipped with an active pre-chamber ignition system and a flexible intake and exhaust VVT device. The very high air-flow demand of the ultra-lean architecture requires the adoption of a two-stage boosting system composed of a variable geometry turbocharger, coupled in series to an E-compressor (E-Comp). The presence of an electrically driven compressor implies that the energy input to the engine is not just the chemical energy given by the fuel. The energy input in the E-Comp also plays a role on the effective engine efficiency. The overall boost level can be then shared between the turbocharger compressor and the E-Comp, depending on the electrical energy input given to the second one. The latter is provided by the same battery pack of the hybrid propulsion system. In other words, the possibility to modify the overall boost level sharing between the LP compressor and the E-comp opens the way to two calibration variants. From one side, an intensive use of the E-Comp leads to a reduced backpressure and an improved brake thermal efficiency (BTE). A more frequent

battery discharge/recharge is, however, to be expected along a driving mission. The opposite situation occurs when the boost level is mainly provided by the standard turbocharger. The main objective of this paper is hence to provide some insights on the calibration and energy management strategy, which provides the minimum CO2 and pollutant emission along representative vehicle missions.

To this aim, the pre-chamber engine under investigation is schematized in a 1D modeling framework. An "in-house" quasi-dimensional model is employed to reproduce the main in-cylinder phenomena, such as mixture formation, combustion, turbulence, and emissions. Its reliability has been validated in a prototype single-cylinder unit in previous activities [19,20]. Here, the same model, with few enhancements, is applied to foresee the behavior of the multi-cylinder engine as well.

The engine control parameters are set by coupling the 1D model to an external optimizer. The latter is instructed to maximize the engine BTE or, alternatively, a properly defined overall BTE, which also considers the electrical energy absorbed by the E-Comp (strategies #1 and #2). Optimization results are obtained at two engine speeds over a load sweep. They are then utilized to define a number of heuristic rules implemented within several PID controllers in the 1D model, mimicking the mathematical optimization procedure. This rule-based (RB) calibration strategy is preliminary validated through a direct comparison with the results of the external optimizer. Then it is utilized to compute two sets of complete engine performance maps, summarizing the calibration strategies following the maximum BTE and maximum overall BTE (OBTE), respectively. Finally, these two sets are embedded in a vehicle simulation able to quantify the CO2 and pollutant emissions over two different driving cycles, namely the worldwide harmonized light vehicles test cycle (WLTC) and a real driving emission (RDE) cycle. This framework is schematized and summarized in the flowchart in Figure 1.

**Figure 1.** Flowchart of the methodology followed in the work.

The paper is organized as follows: firstly, the ultra-lean engine layout is described and the employed combustion model suitable for a pre-chamber engine is briefly recalled. Then, the numerical RB calibration procedure is illustrated and verified against the results of a multi-purpose optimizer. Finally, the energy management strategy developed by the authors is summarized [21] and the vehicle simulation outcomes are shown. In order to quantify the potential benefit of this HEV combined with this high-efficiency engine, the vehicle emissions are also compared against the current regulation limits (Euro 6d) and the CO2 European targets.

#### **2. Engine Description**

The SI engine under study, outlined in Figure 2, is a prototype power unit constituted by four cylinders and equipped with active pre-chambers. Main features of the investigated engine are summarized in Table 1. The PC features four orifices of 1mm [22], and each cylinder has two intake and exhaust valves. The engine presents a cam phaser for both the intake and the exhaust camshaft. Moreover, two different cams profiles can be selected on the intake camshaft, allowing for an early

valve closure (Standard Miller Lift-SML), or an extremely advanced valve closure (Extreme Miller Lift-EML) [22–24]. The latter is particularly suitable to limit the knock occurrence at high load.

**Figure 2.** Schematic engine layout.

**Table 1.** Engine main features.


The engine is fed with liquid gasoline through four port fuel injectors located just upstream the intake valves (InjMC), and four direct injectors located into each PC (InjPC). A two-stage boosting system provides the needed boosting level to match the prescribed load target. It is composed by a variable geometry turbocharger (LPC and LPT) and a high-pressure (HP) compressor (E-Comp), driven by an electric motor (EM). An electric generator (EG) mounted on the engine shaft provides the electrical energy to recharge the battery (Ba) and to power the EM. The electric fluxes are shown as a dashed line in Figure 2.

As experimentally and numerically demonstrated [18–20], the introduction of a pre-chamber allows to strongly extend the lean burn limit of a conventional SI engine. This engine is designed to operate with very lean air/fuel mixtures to improve the engine efficiency and reduce the NOX emission in a zone as wide as possible for the operating domain.

The presence of a high number of degrees of freedom makes the calibration of such a complex engine very challenging. In particular, the engine has nine control parameters, namely the relative air/fuel ratio (λ) in the chambers (labelled as λ*PC* and λ*MC*), the LPT rack position, the power input to the E-Comp, the intake and exhaust cam phases, the spark timing (or the combustion phasing, which can be represented by the angular position when half of the in-cylinder fresh charge is burned—50% of mass fraction burned, MFB50), the throttle valve position, and the SML/EML cam profile.

#### **3. D Engine Model Description and Validation**

Consistently with the engine geometry, a detailed 0D/1D scheme is developed in a commercial modeling framework. The description of the flow inside the intake and exhaust pipes is based on a 1D approach, whereas "in-house" 0D sub-models are used to reproduce the in-cylinder phenomena typical of this novel architecture, such as air/fuel mixture preparation, combustion, turbulence and heat transfer. Regarding PFI and DI injections, it is assumed that the 30% of the injected fuel instantaneously evaporates, without considering spray evolution and liquid wall film formation. An empirical correlation is used for the estimation of the mechanical friction losses, dependent on engine speed, combustion phasing, and in-cylinder pressure peak. The steady-state flow coefficients for both the intake and the exhaust valves have been measured on a similar single-cylinder research engine and are here assigned to predict the flow permeability through the cylinder head [19,20,22]. A standard map-based approach is employed to reproduce the boosting system.

Each cylinder is schematized by two 0D volumes, a constant one (PC) and a variable one (MC), connected by means of an equivalent orifice. In particular, the overall cross-sectional area of the four holes is assigned in the above orifice. The orifice model also takes into account a discharge coefficient, fixed to a constant value of 0.65, which gives the best agreement with the PC-MC pressure difference along the compression stroke. According to a classical filling-emptying technique, mass and energy equations are solved in both volumes to estimate the mass exchange between them, based on pressure difference and overall effective cross-sectional area. The combustion process is described by a rearranged quasi-dimensional fractal model, developed by the authors in the last years [25]. It is based on a two-zone (burned and unburned gases) approach, and it is applied in a similar way in both MC and PC. In this context, the overall burn rate can be expressed as the sum of "classical fractal flame propagation" and a "turbulent jet combustion," expressed by the Equation (1).

$$
\left(\frac{dm\_b}{dt}\right)\_{\text{overall}} = \left(\frac{dm\_b}{dt}\right)\_{\text{fractal}} + \left(\frac{dm\_b}{dt}\right)\_{\text{jet}}\tag{1}
$$

The first term of Equation (1), applied in both PC and MC, describes the burning rate through the fractal approach under the hypothesis that the flame front locally propagates at a laminar speed, *SL*, and that the combustion process is promoted by the turbulence-induced flame wrinkling, Σ, leading to the following expression:

$$
\left(\frac{dm\_b}{dt}\right)\_{fractal} = \rho\_{\text{\textquotedblleft}L} S\_L \Sigma \tag{2}
$$

where ρ*<sup>u</sup>* is the unburned gas density and *AL* is the laminar flame front area. The second term of Equation (1) models the burning rate contribution in the MC due to the turbulent jets, under the hypothesis that the jets entrain fresh charge (air and fuel) and rapidly burn it, in turn releasing heat. This is described by the equations below:

$$\left(\frac{dm\_b}{dt}\right)\_{jet} = \frac{dm\_{b,entr}}{dt} = \frac{m\_{entr} - m\_{b,entr}}{\tau}; \quad \tau = \frac{\Lambda\_T}{S\_L} \tag{3}$$

$$\frac{dm\_{\rm entr}}{dt} = \varepsilon\_{\rm jet} \dot{m}\_{\rm jet} \sqrt{\frac{\rho \text{pc}}{\rho \text{pC}}} \tag{4}$$

The jet-induced burning rate is assumed proportional to the difference between the current entrained mass (*mentr*) and its burned portion (*mb*,*entr*). The Equation (3) resembles the well-known eddy burn-up approach [26,27], for which τ is a characteristic time scale. This last is calculated as the ratio between the Taylor length scale, Λ*T*, and the laminar flame speed, *SL*. Finally, the total entrained mass, *mentr* is computed by the integration of its time derivative, according to Equation (4), which in turn is estimated by the semiempirical correlation reported in [28]. This entrained mass is a function of the mass flow rate from the PC, . *mjet*, on a tuning constant *cjet* ad on the density ratio between PC and MC. Equally, the burned entrained mass, (*mb*,*entr*), is estimated by the integration of Equation (3).

The summarized combustion model requires the evaluation of a few parameters—namely the flame wrinkling Σ and the Taylor scale Λ*T*—which depend on the turbulent field established in the combustion chamber. An "in-house" developed phenomenological K-k-T turbulence sub-model is employed to estimate those quantities [29]. The turbulence model takes into account the energy flux transferred from the mean flow kinetic energy, K, to the turbulent one, k. The model includes proper balance equations for tumble and swirl angular momentum, T and S. The model is applied to both PC and MC, and also considers the turbulence production induced by the incoming/outcoming flow through the orifices [20].

The description of the knock phenomenon is based on the calculation of the auto-ignition (AI) delay time of the air/fuel mixture in the unburned zone. In particular, AI is carried out with a tabulated approach [30], employing preliminary auto-ignition chemical kinetic simulations performed in a homogeneous reactor at a constant pressure. The kinetic scheme adopted for the estimation of the table includes five elements, 201 species, and 1548 reactions [31] and an additional skeletal sub-mechanism for toluene oxidation to handle a toluene reference fuel. The AI time, τ*AI*, is stored in the table as a function of pressure, temperature, equivalence ratio, and residual content. As known, due to the presence of the cyclic dispersion, faster-than-average pressure cycles may induce the presence of a stochastic soft-knock, even when the average cycle is working in a safe, non-knocking combustion. Since, in the adopted approach, the AI integral, Equation (5), is computed on the average cycle, a small safety margin is applied, and the knock event is considered to occur when the AI integral exceeds a tunable lower-than-unity threshold level.

$$\int \frac{dt}{\pi\_{AI}}\tag{5}$$

The model allows to estimate the regulated cylinder-out emissions, namely CO, HC, and NO. For the evaluation of CO and NO, a multizone-zone approach is applied in both PC and MC to estimate the burned gas stratification. Each burned parcel is compressed/expanded adiabatically according to the in-cylinder pressure. The estimated local zone temperatures are employed within the well-known extended Zeldovich mechanism for the estimation of the NO [32], whereas the CO is computed with a two-step reaction scheme [33]. The experimental findings show that in ultra-lean engines the noxious emissions are mainly composed of NO2, therefore in the adopted approach the NO production derived by the Zeldovich mechanism is assumed to completely oxidize the NO2 when the burned gases evolve along the exhaust pipes. Concerning the HC, the crevice contribution and the post-oxidation is considered only for the MC, neglecting the PC contribution. In the employed model, the HC filling and emptying within the crevice regions is followed [34]. The temperature in this volume is considered to be the same as the piston wall, while the pressure equals the one in the cylinder. When unburnt fuel is released from the crevices, a proper boundary layer temperature is computed, at which the released HC partly oxidize according to the kinetic rate proposed in [35].

The quasi-dimensional combustion model was tuned with reference to a research engine [19,20], having geometrical and architectural characteristics similar to the four-cylinder engine under study for which experimental data are not yet available. The model was tuned by a trial-and-error procedure aiming to obtain the best agreement with pressure signals in both PC and MC. A unique set of tuning constants was selected for all the operating conditions, characterized by very different speeds, loads, and air/fuel ratios (λ up to 2.4 in the main chamber) [19,20].

#### **4. Optimization Approach for Engine Calibration**

With the aim of calibrating the engine, a value for each of the nine control parameters listed in the "Engine Description" section must be specified in each operating condition. The purpose of this task is the optimization of the engine performance over the entire operating domain.

The calibration procedure here targets the maximum possible efficiency, fulfilling several constraints—i.e., maximum pressure inside MC and PC, knock intensity, turbo speed, boost level, E-Comp power, internal EGR amount, etc.—prescribed to confine thermal and mechanical stresses and operating reliability for the engine and for its components.

To find the optimal set of the control parameters at variable speeds and loads, an automatic optimization is performed, realized by using an external multi-purpose optimizer. The optimization is performed considering two different control strategies with differentiated fuel and electric consumptions. Particularly, the first control strategy (strategy #1) aims to maximize the conventional BTE, namely:

$$\text{BTE} = \frac{P\_{\text{ICE}}}{\dot{m}\_f \text{LHV}} \tag{6}$$

where *PICE* is the brake power at the engine shaft, *m˙ <sup>f</sup>* is the total injected fuel flow rate, and LHV is the lower heating value of the fuel.

The second one (strategy #2) is oriented to the optimization of the overall BTE:

$$\text{OBTE} = \frac{P\_{\text{ICE}} - P\_{\text{HPC}} / \eta\_{\text{HPC}}}{\dot{m}\_f \text{LHV}} \tag{7}$$

where *P*HPC is the mechanical power at the E-Comp shaft and ηHPC is the electro-mechanical efficiency of the E-Comp. Based on its definition, the OBTE also considers the electric energy required by the E-Comp and drained from the battery.

For both examined control strategies, the optimization is realized along a brake mean effective pressure (BMEP) sweep at two different engine speeds, namely 2000 and 3000 rpm. Low/medium engine rotational speeds are selected because of their relevance regarding the vehicle homologation cycles. The load sweep is automatically realized by alternatively maximizing or minimizing the BMEP, together with the BTE/OBTE, hence performing four different multi-objective problems for each engine speed. Each optimization problem hence identifies a single branch of the BTE/BMEP or OBTE/BMEP Pareto frontiers.

The selected optimization tool is the genetic MOGA-II algorithm, which represents the best choice for the case of a multi-variable/multi-objective problem. The nine control parameters previously listed, apart from the relative air-fuel ratio in the main chamber which is always metered to fix a value of two, represent the independent variables of the optimization problem. Each variable is changed within a prefixed allowable range. The selection of the intake valve lift profile is indeed handled by introducing a fictitious discrete variable, associating values of 0/1 to the Standard/Extreme Miller valve lifts, respectively.

The logical scheme of the optimization process was already reported in a previous author's work [36], and here it is partially modified by the addition of SML/EML and EVC (exhaust valve closure) variables. The optimizer algorithm, at each step with an iterative process, selects the eight engine control variables, passing them to the 1D engine model to evaluate the variables of interest. At the end of each simulation, the estimated values of the objective function are returned to the optimizer to start the next iteration. This process pursues as long as the optimal levels, being part of the Pareto frontier, are achieved.

As said, the optimizer evaluations are post-processed and filtered to check that some monitored variables do not exceed prescribed threshold levels. As an example, the maximum in-cylinder pressure is limited to 180 bar, the maximum plenum pressure to 4.5 bar, the maximum E-Comp power to 10 kW, the maximum AI integral to 0.8, etc.

#### **5. Assessment between RB and Optimizer Calibrations**

The previously described optimization procedure is too time consuming to be directly employed to compute the complete engine performance map. To speed up the numerical calibration, a set of heuristic rules were defined to reproduce as best as possible the behavior of the mathematical optimization (rule-based calibration). The RB was applied by including a network of logical switches, PID controllers, and "math functions" in the 0D/1D environment.

Preliminarily, a prefixed running line was chosen on the LPC map as a tradeoff between an appropriate surge margin and a suitably high-pressure ratio to avoid an excessive E-Comp shaft power absorption (below 10 kW). Under these operating conditions, the E-Comp running line was directly identified to reach the HP boost level necessary to fulfill a prescribed full load (FL) target, depicted in red in Figure 3. By changing the LPC running line, a different sharing of the overall boost level between the LPC and the E-Comp would have been accordingly identified. However, whatever the choice for the LPC running line, at very high engine speeds (usually over 4500 rpm) for the selected boosting devices, the constraints on the maximum E-Comp power and/or maximum boost level could not be fulfilled. For this reason, the only way to achieve the FL target was to inject more fuel in the main chamber, resulting in a λMC < 2. The L2 line in Figure 3 identifies the BMEP region below which a λMC = 2 operation can be sustained (light yellow domain), while the area above the L2 line defines the operating domain with λMC < 2 (highlighted in orange).

**Figure 3.** BMEP–rpm map containing the BMEP target (FL), and the L2 and BB lines.

In the low-load region, an additional line is defined corresponding to the BMEP levels reached with the LPT fully open and no power given to the E-Comp, named Base Boost Line (BB in cyan in Figure 3). In these conditions, the E-Comp is bypassed, acting on dedicated control valves. Above the BB line, the engine load is regulated by a progressive closing of the LPT and a proportional increase of the E-Comp power. Below the BB line, the LP and HP compressor settings continue to be fixed, while the load is only regulated by the progressive closure of the throttle valve.

Concerning the VVT position of intake and exhaust camshafts, simple rules are defined to select an advanced IVC (intake valve closure) at low load to have a pumping losses reduction. The intake valve advancing is, however, limited at very low load to avoid an internal EGR rate above 20%. In this way, a stable combustion process can be guaranteed in the real engine at the expense of a certain pumping losses increase.

In the mid-load range, the intake VVT is moved towards a later IVC to allow for increased air flow rates, while, at a very high load, an early IVC is once again preferred to reduce the effective compression ratio for knock control. The exhaust valve is accordingly controlled to limit the valve overlap and/or to increase the expansion work at a low load.

The actual positions of the BB and L2 lines also depend on the selection of the intake lift profile, namely the SML or the EML. To reduce pumping losses, the EML is always specified when the computed airflow rate is below 100 kg/h. By analyzing the results of the previous mathematical optimization, it was indeed found that for strategy #1 (max BTE), the SML is only selected below a specific BMEP limit (namely 6.5 bar), while it is always selected above 100 kg/h airflow rate in the strategy #2 (max OBTE). This can be easily explained by observing that the SML requires a lower boost level for the same load. For this reason, a lower E-Comp power is required, which improves the

OBTE. Moreover, to better reproduce the optimization results, a different running line is selected for the strategy #2, moving the boost sharing mainly on the LPC and reducing the E-Comp contribution.

The proposed rule-based methodology is quite refined to provide engine operations very close to the optimal BTE and OBTE, as confirmed in the following. In addition, the advantage of this RB strategy is related not only to its easy implementation in the used 1D software, but also to its much faster execution compared to a set of multi-objective optimizations.

To quantify the reliability of the two implemented calibration strategies, their outcomes were compared to the optimization results. The optimization outcomes are shown in the next figures with open symbols and label Opt and refer to the outcomes belonging to the Pareto frontier.

Figure 4 reports the comparisons of BTE and OBTE at 2000 rpm for the considered calibration strategies (#1 and #2). In both cases, a considerable agreement was realized along the BMEP sweep between Opt and RB procedures. Figure 4 also shows the potentiality of the analyzed engine architecture to reach a very high BTE level of about 50% at the higher load. In the latter conditions, lower OBTE is reached due to the relevant power request by E-Comp (see OBTE definition in Equation (5)). Indeed, OBTE assumes a maximum value around 42% at medium loads, with a slight reduction at the highest BMEP levels.

**Figure 4.** Brake thermal efficiency (BTE) and overall brake thermal efficiency (OBTE) comparison in a BMEP sweep at 2000 rpm for strategies #1 (**a**) and #2 (**b**).

Quite satisfactorily, RB/Opt agreements were found for the most relevant engine calibration variables depicted in Figures 5 and 6. The first assessment regards the intake valve strategy and the intake/exhaust valve timings, which are plotted in Figure 5. Considering the strategy #1, the Extreme Miller was mainly preferred all along the BMEP sweep, while in the strategy #2 at medium/high loads, the Standard Miller strategy was chosen.

On a general viewpoint at both low and high loads, IVC is selected near to the most advanced timing setting. This option, as said, is preferred for minimizing the pumping losses at low loads and mitigating the knock occurrence at high loads [37]. In the medium BMEP range, especially for strategy #2, IVC is slightly delayed, allowing for higher air flow rates and resulting in a better agreement with the Opt outcomes. EVC indeed exhibits an approximately constant trend with load.

**Figure 5.** Selected intake valve lift profile, IVC and EVC, comparison in a BMEP sweep at 2000 rpm for strategies #1 (**a**) and #2 (**b**).

**Figure 6.** Low/High pressure ratios and MFB50 comparison in a BMEP sweep at 2000 rpm for strategies #1 (**a**) and #2 (**b**).

Figure 6 shows the sharing of the overall engine boosting between the two compressors, depicting the HP and the LP pressure ratios. As expected, strategy #1 favored a higher boosting by the E-Comp instead of the LPC. Conversely, for strategy #2 the boosting was principally realized by the LPC; instead, the E-Comp worked by having a pressure ratio as low as possible. This choice involved an increased engine exhaust backpressure, but it resulted in a lower power absorbed by the E-Comp.

The results in Figure 6 also show that HP and LP pressure ratios assume unit values below 5 bar BMEP, suggesting a throttle-based load control. Above 5 bar, both LP and HP boosting are progressively regulated by a partial closing of the LPT rack and by simultaneously increasing the E-Comp power.

It is worth underlining that the RB strategy #1 between 4–10 bar shows some discrepancies if compared to the optimal outcomes. More specifically, the RB approach determined a delayed IVC, which is compensated by a lower HP boost pressure. Nevertheless, this reflects a slight BTE lowering compared to the optimizer output (less than 1%).

As a further assessment between RB and Opt methodologies, in Figure 6 the BMEP trends of MFB50 are reported, highlighting a satisfactory RB/Opt agreement. For this combustion parameter, the optimal value is usually around 6–10 CAD AFTDC. Indeed, in the case of later phasing, the combustion develops during the expansion stroke, leading to a less effective work exchange between the piston and the gas, while, in the case of too advanced MFB50, a relevant portion of the combustion process takes place during the compression stroke, increasing the compression work exerted by the piston on the fluid. An almost flat trend was obtained for strategy #1 in a wide load range (2–15 bar BMEP) with most of the points in the band 5–10 CAD AFTDC. A slightly delayed MFB50 at both higher and lower BMEP values was detected. At high loads, this is due to knock mitigation, while at low loads, it is because of the combustion lengthening and of the limitation on the maximum allowable spark advance. For the strategy #2, MFB50 showed a different trend with a more pronounced delay of the MFB50 at increasing loads (above 10 bar BMEP). This, of course, is the only possibility for knock control once the Standard Miller profile for the intake valves has been selected. This quite different behavior was very well captured by the RB strategies.

Although not discussed here for sake of brevity, another outcome of the optimization is that a close-to-stoichiometric mixture is preferred in the pre-chamber for whatever is the BMEP. A certain mixture leaning is required at high loads to limit the pressure peak in the PC.

To further prove the consistency of the RB calibration, Figures 7 and 8 present additional comparisons between RB and Opt strategies, referring to the engine speed of 3000 rpm. Once again, a satisfactory RB/Opt agreement for all the considered variables was reached all along the BMEP sweep. Figure 7 shows that, for this higher speed, the optimal control involved more relevant modifications of intake valve lift profiles and IVC/EVC timings at changing loads.

For strategy #1, different from the previous speed, the Standard Miller was selected in a restricted low/medium load range (4–5.5 bar BMEP) to enhance the effective compression ratio and hence the efficiency in the light of a minor knock tendency at this speed. However, the Extreme Miller was once again preferred as soon as the load increases (Figure 7). For strategy #2, an earlier transition from the Extreme Miller to the Standard Miller by increasing the BMEP was detected if compared to the results at 2000 rpm. Once again, at medium loads the IVC was delayed, allowing for increased airflow rates. EVC trends appear almost load insensitive regardless of the calibration strategy, and a modest discrepancy between RB and Opt procedures was detected especially at very low BMEPs (0–2 bar in Figure 7). This is probably due to a more refined calibration given by the optimizer, which tends to reduce the valve overlap with the aim to limit the internal EGR even more than the logics of the RB approach. Referring to the results plotted in Figure 8, similar considerations to the ones already carried out for the case at 2000 rpm still hold.

Concerning the combustion phasing differences between 2000 and 3000 rpm, a greater delayed MFB50 was observed at low loads for the engine speed of 3000 rpm, mainly ascribed to a pronounced combustion lengthening of its initial stage due to an excessive rich mixture in the pre-chamber at the spark event. This is due to the lower bound of the PC injection duration of 300 μs (see Table 1). At high loads, a less pronounced MFB50 delay occurred for the strategy #2 due to the lower knock tendency at 3000 rpm. In summary, the above discussed results highlight that the RB calibration demonstrates to be extremely trustworthy for various load levels and for two different speeds, allowing to reach efficiencies and control variables very close to the ones obtained by the optimizer-based methodology. In the light of the latter consideration, the developed RB calibration was broadened to the full engine speed/load range in order to investigate the close-to-optimal engine performance in the whole operating domain, as reported hereinafter.

**Figure 7.** Intake valve lift profile, IVC, and EVC comparison in a BMEP sweep at 3000 rpm for strategies #1 (**a**) and #2 (**b**).

**Figure 8.** Low/High pressure ratios and MFB50 comparison in a BMEP sweep at 3000 rpm for strategies #1 (**a**) and #2 (**b**).

#### **6. Map of Engine Performance**

According to the RB control strategy, the entire engine operating map was computed (consisting in 143 operating points—11 rpm × 13 BMEP) for the two different calibration strategies above discussed. In Figure 9, the BTE maps for calibration strategies #1 (a) and #2 (b) are shown. The maximum levels for both strategies occur at 2000 rpm, close to the full load region. The OBTE maps for the different strategies are shown in Figure 10. For these, the maxima occur at medium speeds (2000–2500 rpm) and medium–high loads (8–14 bar BMEP). According to the selected optimization targets, a 3% higher peak BTE was obtained in the strategy #1 (Figure 9a), while a 1% higher peak OBTE was obtained in the strategy #2 (Figure 10b).

**Figure 9.** Map of BTE (%) for strategies #1 (**a**) and #2 (**b**).

**Figure 10.** Map of OBTE (%) for strategies #1 (**a**) and #2 (**b**).

The gap between BTE and OBTE for the strategy #1 was about 7% in the OBTE peak zone, and gradually decreased at higher speeds and low loads. A similar behavior characterizes the strategy #2, having in the peak zone a gap of about 3–4%. This is due to the power absorbed by the E-Comp, whose maps are reported in Figure 11. For strategy #1, the E-Comp power gets to the upper bound of 10 kW, which is about the 8% of the power rated by the engine, while for strategy #2, it only gets a maximum of about 7 kW.

**Figure 11.** Map of E-Comp power consumption (kW) for strategies #1 (**a**) and #2 (**b**).

For the considered engine, featuring an ultra-lean air/fuel mixture, both strategies involved a very high plenum pressure, which is mandatory to achieve the load target. As an example, Figure 12a shows the boost pressure map for the strategy #1, which highlights a peak of 4.3 bar at the highest speeds and loads. The intake plenum pressure was conserved above 1.0 bar in most of the map (reduced pumping losses), and above 0.5 bar even in the area at lowest BMEP. The lambda target of 2 in the main chamber was reached in most of the operating domain, as shown in Figure 12b, except in the maximum power area. This is due to the attainment of the maximum E-Comp power, which obliges to reduce the mixture leaning (λMC < 2) to get the load. Similar boost pressure and lambda maps are also found for the strategy #2.

**Figure 12.** Map of intake plenum pressure (bar) (**a**) and (**b**) map of λMC (-) for strategy #1.

The predicted brake-specific emissions of NOx and HC are shown in Figure 13a,b, respectively, still referring to strategy #1. With the aim to limit the temperature peaks in the burned zone of the combustion chamber preventing the formation of the above pollutants, the ultra-lean combustion concept was used. The NOx production was very low in most parts of the operating map, except for the high speed and load zone. This is due to the lower mixture leaning in the above-mentioned zone, compared to the λMC = 2 target (see λMC contours in Figure 12b). Due to the higher burned temperature, if compared to the one in the MC, for this engine the NOx production was prevailing in the pre-chamber. Despite the large excess air, the HC production is not negligible over the whole engine domain (Figure 13b). This is due to the lean burning concept, which, compared to a conventional stoichiometric engine, reduces the in-cylinder temperatures and consequently makes the post-oxidation phenomenon less effective. Looking to Figure 13b, the shape of the HC iso-contours is justified by a faster post-oxidation at increasing speed (less time for wall heat losses) and load (lower percent wall heat losses).

**Figure 13.** Map of (**a**) brake specific NOx emission (g/kWh) and (**b**) map of brake-specific HC emission (g/kWh) for strategy #1.

#### **7. Vehicle Simulations**

#### *7.1. Hybrid Vehicle Architecture*

The main characteristics of the vehicle here investigated are reported in Table 2 [38]. The vehicle is a HEV, belonging to the C-class. It is characterized by a combined parallel/series powertrain, schematized in Figure 14. The powertrain is composed of an ICE, two electric motor/generator units, (EM and EG), a battery (Ba), three clutches (Cl1–3) and two gear boxes (GB1–2).


**Table 2.** Main characteristics of the tested HEV.

**Table 2.** *Cont*.


**Figure 14.** Powertrain schematic of the tested HEV.

#### *7.2. Simulation Platform*

The vehicle simulations were carried out by a "in-house developed" software, implemented in Fortran language (UniNa Vehicle Simulation-UNVS) [21]. It is based on a "backward-facing" (quasi-static) method [39]. All of the powertrain schematic components in Figure 14 were defined by a lumped-parameter method. The approach to describe the thermal unit was quasi-steady map based. In particular, the BSFC maps, derived from the calibration strategies #1 and #2, were implemented in the vehicle simulation, together with the curve of maximum shaft torque. In addition to the BSFC maps, the vehicle model handled the corresponding E-Comp power consumption maps. The electrical absorption of the E-Comp is considered as an additional load on the battery. The ICE thermal transient was neglected in the simulations, resulting in a null fuel consumption penalization at the cold start. The exhaust gas after-treatment system was not modeled, so the predicted emissions have to be considered as "engine-out" values. The electric unit efficiencies were described once again by a map-based approach, dependent on speed and shaft torque. The battery was treated by a conventional State of Charge (SoC) model [40]. In order to estimate the Joule-effect losses, an internal resistance was set, neglecting the possible variations in changing SoC and temperature. The mechanical losses in the gearboxes were evaluated assuming constant efficiencies (equal to 0.97).

#### *7.3. E*ffi*cient Thermal Electric Skipping Strategy*

The vehicle control strategy here adopted, if compared to the well assessed strategies (namely Pontryagin minimum principle [41], dynamic programming [42] or ECMS [43]), is designed to be simplified. The conventional power-split logics, typical of the abovementioned strategies, is substituted with an alternative utilization of the thermal and electric units for the vehicle driving (Efficient Thermal Electric Skipping Strategy–ETESS [21]). The basic concept behind ETESS is an alternate use of the electric units and thermal engine so to fulfill the power demand at the vehicle wheels. The choice between the two modes depends, at each time, on the comparison between the effective fuel rate of the thermal engine needed to fully satisfy the power demand at the vehicle wheels, and an equivalent fuel rate related to a pure electric driving of the vehicle. The ETESS is featured by a much-reduced computational effort if compared to most common optimization strategies, leading to

similar performance in terms of kilometric consumption and emission [21]. The results shown in the next section were obtained by the off-line variant of the ETESS, where the control is tuned to realize the battery energy balance between the beginning and the end of the driving mission [44].

#### *7.4. Simualtion Raster and Results Discussion*

The vehicle simulations were carried out for two driving cycles, namely the WLTC and a real driving cycle (RDC) [45]. The results of CO2 and pollutant emissions (grams per kilometer) are represented under the form of bar charts in Figures 15 and 16. Each couple of bars compares the results of the calibration strategy #1 (red) and strategy #2 (blue), reporting on the top percent difference.

**Figure 15.** Strategies assessment of CO2 emission and their percent difference along WLTC and RDC.

**Figure 16.** Strategies assessment of NOx emission (**a**) and of HC emission (**b**) and their percent difference along WLTC and RDC.

Starting from the CO2 results shown in Figure 15, a first consideration is that the emissions are below the target imposed by the European Union for the 2021 of 95 g/km, but above the one of the 2025, equal to about 81 g/km or whatever is the ICE control strategy. These results indicate that the HEV vehicle here analyzed, even equipped with a very efficient thermal unit, was not enough to match the EU 2025 target and a plug-in HEV variant is mandatory to get this. The pollutant emission outcomes, as stated above, refer to engine-out values compared to the tailpipe reference standards imposed by the legislation. The potential abatement due to an aftertreatment system is not explicitly considered. However, the comparison can also help to evaluate and design the most adequate exhaust after-treatment system (EATS) for the considered application. Looking at the values for the nitrogen oxides (Figure 16a), it can be seen that the emissions are within the bounds. The HC values (Figure 16b), for both the considered cycles and calibration strategies, are much higher than the limit imposed from the Euro 6d, highlighting the need for an oxidizing EATS for regulation compliance. Anyway, the technology of an oxidizing EATS at the current state-of-art is very robust, low-cost and reliable, not representing a substantial challenge for a real on-vehicle application. The engine calibration, either strategy #1 or #2, does not seem to exert a substantial influence on the pollutant emissions. It is worth noting that the ETESS control strategy was not conceived to minimize the pollutant emissions, but the fuel consumption and hence the CO2 emissions.

The comparison presented underlines that the strategy #2 allows slightly lower CO2 emissions for both the considered driving cycles. To clarify this result, the instantaneous trends of some representative control and performance parameters are plotted in Figure 17 with reference to the WLTC. The trends underline that the powertrain control is practically the same for both the ICE calibrations, resulting in almost superimposed traces of ICE power (b), EM power (c) and gear number (g). Some differences emerge for the trends of E-Comp power (d) and SoC (e), which reflect slight misalignments of fuel consumption (f) and battery power (h) profiles. Note that a positive value of battery power corresponds to an energy flux drained from the battery.

To better highlight the differences among the two engine calibrations, the plots of SoC, fuel rate and battery power are reproposed in Figure 18 with a magnified time scale in the most load-demanding portion of the driving cycle, namely between 1450–1750 s. During this period, the fuel rate is lower for the calibration strategy #1 if compared to the strategy #2, while the opposite occurs for the battery power consumption due to a higher E-Comp power absorption (Figure 17d). This reflects in a faster decrease of the SoC, which leads to reduced electric driving in the first portion of the cycle where the load demand is lower. For this reason, the ICE works for a longer time in the operating conditions, namely at low load/speed, where its efficiency is well below its maximum level. On the contrary, the calibration strategy #2 allows longer pure electric driving in the first portion of the cycle, avoiding for the thermal engine operations with a very reduced efficiency.

A similar behavior emerged for the RDC, as shown in Figure 19. This reports the vehicle speed profile (a) and the SoC trend (b). Note that the initial SoC level was chosen to get simultaneously the optimal control and battery energy balance along the cycle. Such a choice is a consequence of the simplified approach for the battery treatment, where the losses are not dependent on the instantaneous SoC level and on the possible stress when the extreme SoC bounds are attained.

Once again, during the most load-demanding portion of the cycle (after about 5400 s), the battery discharge for the calibration strategy #1 was faster, and this was compensated by a reduced use of electric driving in the first portion of the cycle. This occurrence globally penalizes the strategy #1 in comparison with the strategy #2.

Under a more general point of view, for a given battery capacity and vehicle mission, the engine calibration, which reduces the energy drain from the battery, allows for a more efficient choice of the electric driving phase along the cycle. In the presented cases, however, a limited benefit has been estimated. An explanation of such reduced fuel consumption differences can be drawn by analyzing the time evolution of the engine operating point. Referring to the RDC, the instantaneous engine operating points, sampled at each second, are shown in Figure 20. They are plotted over the map of OBTE difference between strategies #2 and #1. It can be observed that most of the samples are located in a speed range between 1500 and 2500 rpm and with loads below 8/9 bar BMEP, where the OBTE difference is lower than 0.6% and the E-Comp power is small (Figure 11). Changes between the two calibrations only emerge outside the above operating area where the percentage difference in terms of OBTE increases, contributing to reduce the fuel consumption of strategy #2. Similar considerations can be applied to the WLTC results. Different driving missions and more load demanding are expected to bring out more relevant differences between the considered strategies.

**Figure 17.** Instantaneous trends of vehicle speed (**a**), ICE power (**b**), EM power (**c**), E-Comp power (**d**), SoC (**e**), fuel rate (**f**), gear number (**g**), and battery power (**h**) along WLTC.

Even if the advantage of the strategy #2 cannot be considered yet as a general result due to the need of further verifications with different powertrain/vehicle types, it seems that a reduced battery consumption, even if the thermal engine will not give its maximum efficiency performance, is the path to improve the hybrid vehicle fuel economy. As stated above, further analyses will be performed in the next development of this activity to confirm the generality of this result.

**Figure 18.** Instantaneous trends of SoC (**a**), fuel rate (**b**), and battery power (**c**) along the last portion of WLTC.

**Figure 19.** Instantaneous trends vehicle speed (**a**) and SoC (**b**) along the RDC.

**Figure 20.** Instantaneous engine operating point along the RDC over the OBTE difference between strategy #2 and #1.

#### **8. Conclusions**

In this paper, a refined study of an innovative engine prototype was provided, aiming to estimate the impact of the engine calibration on CO2 and pollutant emission of a hybrid vehicle equipped with the considered power unit, consisting in a two-stage boosted ultra-lean pre-chamber SI engine.

As a first step, a physical model for the prediction of the burn rate for this engine architecture was recalled, able to sense its main geometrical parameters and operating conditions. The combustion model, validated in previous authors' activities, was integrated with advanced sub-models of turbulence, heat transfer, knock, and engine-out pollutant emissions formation.

The model was employed to estimate complete engine performance maps once a rule-based numerical pre-calibration had been set up. In particular, two variants of the calibration procedure were proposed, aiming to maximize either the BTE or the OBTE; this last option also considered the E-Comp power consumption. Despite the great number of control parameters, the RB approach was successfully verified against the results of a more complex and time-consuming procedure where the 1D engine model was coupled to an external optimizer built on a genetic algorithm (MOGA-II).

The estimated performance maps revealed the potential of this advanced engine architecture to reach very high efficiency levels. Depending on the calibration strategy, a 49% BTE peak was estimated at 2000 rpm (full load), and a 43% OBTE peak was also gained at a lower load. The analyses also underline low cylinder-out specific NOx emissions, thanks to the ultra-lean combustion concept and relevant cylinder-out specific HC emissions, due to an ineffective post-oxidation process.

The computed engine performance maps were employed in vehicle simulations to assess these strategies along two representative driving cycles, namely WLTC and an RDC. The comparisons underline that the strategy #2 allows slightly lower CO2 emissions for both the considered driving cycles. This is due to a more extended pure electric driving during the low-speed portion of the driving cycles, allowed by a reduced electric consumption in the high-speed driving mission. In this way, the time in which the engine operates with a reduced efficiency was shortened. Either control strategies do not substantially show a clear impact on the pollutant emissions.

The discussed methodology represents a potential powerful tool to forecast both engine and vehicle performance. Once tested and validated against experimental data, this procedure will contribute to supporting and driving the development phase of a high-efficient HEV.

**Author Contributions:** Conceptualization, F.B. and V.D.B.; Methodology, F.B. and V.D.B.; Software, E.M.; Validation, E.M., L.T. and D.T.; Data Curation, E.M., L.T. and D.T.; Writing—Original Draft Preparation, F.B., V.D.B., E.M., L.T. and D.T.; Writing–Review & Editing, V.D.B. and E.M.; Visualization, E.M., L.T. and D.T.; Supervision, F.B. and V.D.B.; Project Administration, V.D.B.; Funding Acquisition, F.B. and V.D.B. All authors have read and agreed to the published version of the manuscript.

**Funding:** This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No 724084.

**Conflicts of Interest:** The authors declare no conflict of interest.

#### **Definitions**/**Abbreviations**



#### **Symbols**



#### **Greeks**


#### **References**


© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
