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 CO
2 formation during the vehicle manufacturing and the one associated to electrical energy production and distribution. Under this point of view, the overall CO
2 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 CO
2 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 CO
2 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.
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 CO
2 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,
23,
24]. The latter is particularly suitable to limit the knock occurrence at high load.
The engine is fed with liquid gasoline through four port fuel injectors located just upstream the intake valves (Inj
MC), and four direct injectors located into each PC (Inj
PC). 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,
19,
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 NO
X 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 and ), 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).
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,
, and that the combustion process is promoted by the turbulence-induced flame wrinkling, Σ, leading to the following expression:
where
ρu 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:
The jet-induced burning rate is assumed proportional to the difference between the current entrained mass (
) and its burned portion (
. 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,
, and the laminar flame speed,
. Finally, the total entrained mass,
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,
, on a tuning constant
ad on the density ratio between PC and MC. Equally, the burned entrained mass, (
, 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
—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,
, 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.
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 NO
2, therefore in the adopted approach the NO production derived by the Zeldovich mechanism is assumed to completely oxidize the NO
2 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:
where
PICE is the brake power at the engine shaft,
ṁf 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:
where
PHPC 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).
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.
Quite satisfactorily, RB/Opt agreements were found for the most relevant engine calibration variables depicted in
Figure 5 and
Figure 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 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 MFB
50 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 MFB
50, 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 MFB
50 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, MFB
50 showed a different trend with a more pronounced delay of the MFB
50 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,
Figure 7 and
Figure 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 MFB
50 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 MFB
50 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.
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).
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.
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.
The predicted brake-specific emissions of NO
x 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 NO
x 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 NO
x 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).
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.