Conceptual Study on Car Acceleration Strategies to Minimize Travel Time, Fuel Consumption, and CO 2 -CO Emissions

: A conceptual study was performed on intelligent driving acceleration strategies for vehicles equipped with internal combustion engines. Two archetypal acceleration scenarios of highway driving and urban driving were prescribed. Three trajectories were considered for each scenario. They involved (a) nearly constant acceleration, (b) fast acceleration first and slow acceleration later, and (c) slow acceleration first and fast acceleration later. The selected vehicle was a generic European small–medium passenger car. Engine inlet pressure and ignition time were optimized along each trajectory to minimize fuel consumption, CO, and CO 2 emissions, and travel time. The optimization process involved a methodological approach based on the higher-order singular value decomposition of the tensor form of the engine model. The optimized trajectories were analyzed and compared among themselves. Conceptual acceleration design guidelines for intelligent driving were provided that could be of interest when integrating vehicle/engine performance into the surrounding traffic flow.


Introduction 1.Overview
Reducing fuel consumption and pollutant emissions are critical issues for internal combustion engines in the context of car intelligent driving (including the case of battery electric-hybrid vehicles).The reason is that transition to fully electric or hydrogen-based transport systems may take longer than expected.During this transition, highly efficient optimized vehicles based on internal combustion engines may represent a feasible compromise.
The operation procedure of an internal combustion engine begins with a certain power demand that depends on the vehicle trajectory, and follows with an adjustment of the engine control parameters to satisfy that demand.This means that fuel consumption and pollutant emissions depend on both external (vehicle trajectory) and internal (engine variables) factors that, in turn, depend on a large number of parameters.
In this context of complex operation, it is doubtful that a human driver would be able to develop real time strategies to deal, optimally, with a non-linear problem that contains a large set of variables.Nowadays, the prevailing trend in the field is to move towards intelligent driving strategies by means of either using driver assisting engine control unit (ECU) software, or by trying to move to an autonomous concept vehicle.Intelligent driving approaches are very promising but require a large amount of R&D effort.The reason is that they involve (a) modelization of vehicle powerplant and associated mechanical and electrical/electronic systems, (b) development of a driving framework for a variety of drivers and trajectories, and interaction with other drivers, and (c) integration of said

Literature Review
Regarding the relation between acceleration and traffic flow, Choudhury and Islam, [5] have developed a "Latent Leader" approach to model the influence that acceleration decisions have in traffic streams with weak lane discipline (meaning that the driver may have several vehicles in front, and not necessarily a leader vehicle).In their conclusions, the authors point to the importance of space headway in a situation of vehicle deceleration.Goñi-Ros et al. [6] addressed traffic flow optimization at sags by means of acceleration control.Their results show that the optimal solution involves a deceleration-accelerationdeceleration-acceleration strategy.Ivanchev et al. [7] have studied available models of autonomous vehicle longitudinal acceleration in a context of mixed traffic (autonomous vehicles mixed with human drivers) to maximize the speed, safety, comfort, and efficiency of the vehicle population.For this aim, they formulated a multi-objective problem, and found the set of traffic/driving parameters that optimized a certain objective function.Interestingly enough, they state that the most important limitation of their study is the reliability of the human model.Manolis et al. [8] have presented an adaptive-cruise-controlbased strategy for traffic flow efficiency.This is achieved in real time, acting on the driving behavior (including acceleration strength).In their conclusions, increased acceleration at active bottlenecks was presented as one of the optimal strategies.Wang et al. [9] have studied how road curvature ahead and slope variations affect driver's deceleration and acceleration decisions.
Concerning the effect of acceleration on engine performance parameters, Ard et al. [10] present a car-following algorithm that controls acceleration strategies and reduces the engine fuel and/or electric consumption.In their conclusions, the authors report a 30% energy improvement as compared to a realistically calibrated human driver.It is important to note that, in addition to simulations, the authors used specifically adapted experimental vehicles: a gasoline Mazda CX-7, and an electric Nissan Leaf.Bishop et al. [11] developed engine maps of fuel consumption and emissions associated with some specific transient driving cycles.To this aim, they used data provided by the on-board diagnostic system in combination with a dynamic equation of the vehicle that included acceleration and a portable emissions measurement system.Some of their experimental vehicles were a Volkswagen Golf and a Chevrolet Cruze.Guo et al. [12] have developed a hybrid reinforcement learning eco-driving algorithm that accounts for longitudinal acceleration and deceleration as well as for lateral lane-changing maneuvers.The objectives of the algorithm were twofold: (a) to reduce fuel consumption, and (b) to maintain an acceptable travel time.Hibberd et al. [13] tested two haptic accelerator pedal systems, and one multimodal visual-auditory system to develop guidelines related to the fuel efficiency of the accelerator pedal position in the context of an eco-driving in-vehicle assistance system.The study used the driving simulator at Leeds University and considered several scenarios: deceleration (96 km/h to 64 km/h), cruise (64 km/h), and acceleration (64 km/h to 96 km/h).Sun et al. [14] developed a method to estimate vehicle emissions and energy consumption using mobile sensing data and specific vehicle trajectories (microflow).The considered trajectories accounted not only for the traffic state such as queuing and free-flowing, but also for the vehicle driving mode: cruise, idle, acceleration, and deceleration.Tsanakas et al. [15] proposed a method to estimate emissions and fuel consumption using the so-called virtual vehicle trajectories, which are generated out of partially observed traffic data.The engine/car data used in the study followed the previous work of Hausberger et al. [16].These authors experimentally categorized seven different gasoline cars and generated engine maps that included emissions information.The steady experimental maps were corrected, whenever applicable, with so-called transient emission correction factors relevant for transient driving cycles.In their conclusions, Tsanakas et al. [15] provided guidelines to improve the virtual vehicle trajectories technique, and stated a possible range of applications not necessarily related to emission-modelling.Recent experimental studies, performed at both test benches and actual traffic conditions, that relate emission characteristics to acceleration have been published by Kim et al. [17], Chadrashekar et al. [18], and Kuppili et al. [19].In this context, a very comprehensive semi-empirical methodology to link instantaneous acceleration and CO 2 emissions has been published by Suarez et al. [20].Also, it is relevant to refer to the work of Zhang et al. [21], which studied the influence of acceleration on vehicle emissions at intersections using the MOVES emission model (US Environmental Protection Agency, 2009) [22].
Finally, it is relevant to refer to the study presented by Blokpoel et al. [23].The authors focused on the problem of vehicle emissions when stopping and accelerating at traffic lights.The authors began the article by stating that this specific driving condition is one of the main contributing factors to vehicle emissions in urban environments.Then, they presented an adaptive control strategy to minimize CO 2 emissions at traffic lights.The outcome could be summarized as a series of interesting guidelines.First, emissions optimum ratio of delay time versus stops should be 1:165 at 50 km/h.Second, at 70 km/h speed limit, the ratio increases to 1:296.Since it was found that optimal free flow speed was 70 km/h, the maximum speed for optimal emissions with long distances between intersections and not many stops should be higher than 50 km/h.

The Present Study
In line with the references discussed above, the novelty of the present study is that it provides some conceptual guidelines for a somewhat idealized problem.The problem is when accelerating along a longitudinal trajectory, what is the best strategy to minimize fuel consumption and emissions while maintaining a reasonable travel time?To accelerate fast first and slow later?The opposite?Or something intermediate?An additional methodological novelty is the use of higher-order singular value decomposition (HOSVD) for generalization and densification of engine output data that, otherwise, may require an enormous computational effort.It should be remarked, nonetheless, that this study is of a conceptual nature.The aim is to first carry out an exploratory analysis based on a new methodology.As a second step in the near future, a series of additional practical aspects such as wear of certain components, quality of fuel, engine temperature, vehicle maintenance, shift behavior, gearbox, etc., will be considered.One example of the importance of these effects could be found in Crolla [24], where the influence of rolling resistance was highlighted in urban driving.
Regarding organization of the manuscript, the problem is described first (Section 2).Next, the methodology to solve the problem based on the minimization of a multi-objective function (Section 3) is described.Results are presented and discussed in Section 4, and, finally, conclusions and guidelines are given in Section 5.

Problem Definition
The problem consisted of the evaluation of different vehicle trajectories in two acceleration scenarios representative of highway and urban driving.The comparison was performed in terms of fuel consumption, CO 2 and CO emissions, and time needed to cover the trajectories.
Each acceleration scenario had its own constant transmission ratio.The first acceleration scenario (AC1) was high transmission ratio (highway driving), while the second (AC2) was low transmission ratio (urban driving).Three different trajectories (T1, T2, T3) were defined for each scenario.They involved nearly constant acceleration (T1), fast acceleration first and slow acceleration later (T2), and slow acceleration first and fast acceleration later (T3).
For each scenario and trajectory, two engine input parameters, inlet pressure and ignition timing, were optimized as a function of time along the trajectory to minimize four output parameters.The reason for this selection is that inlet pressure and ignition timing are two typical parameters that are acted upon by the engine control unit in a context of intelligent driving (they have a large influence on the thermodynamics cycle).The methodology is, however, generic and could accommodate other sets or input parameters.The output parameters were total fuel consumption, total emissions of both CO and CO 2 , and driving time along the trajectory.A restriction of the optimization process was to keep an acceptable risk of engine knocking.Then, the optimized trajectories were compared among themselves to extract generic guidelines regarding vehicle acceleration strategies.

The Car
A representative small-medium segment European type gasoline car was selected.Its relevant architectural and operational parameters were:

Acceleration Scenarios
Two acceleration scenarios were considered.The first (AC1) addressed acceleration at velocities of the order of 100 km/h along a horizontal plane with a high transmission ratio.The second (AC2) considered acceleration at velocities of the order of 30 km/h along an inclined plane with 5% slope with a low transmission ratio.There is, certainly, some level of arbitrariness in this selection since there is an enormous variety in possible acceleration trajectories in the context of either urban or highway traffic scenarios.However, what is intended is to select some trajectories that may illustrate the development of the methodology while being reasonably representative of the traffic scenarios.
Three different trajectories were defined for each acceleration scenario.

Acceleration Scenario 1
• The car accelerated from 20 m/s (72 km/h) up to 35 m/s (126 km/h) along a horizontal plane 1000 m long, see Figure 1; • Acceleration at both the beginning and end points was zero; • Three velocity trajectories (time laws) were considered.They were modelled via Bezier polynomials, and represented nearly constant acceleration (T1), fast acceleration first and slow acceleration later (T2), and slow acceleration first and fast acceleration later (T3).Each Bezier polynomial had 1000 points; • The transmission ratio was 23 km/h per 1000 rpm.
The idealized dynamic equation of the vehicle in this case (neglecting friction and rolling resistance) was: .
The term on the left, .
W engine , is the engine power.The first term on the right is the power needed to overcome the aerodynamic drag (ρ is the air density, v the vehicle velocity).The third and fourth terms represent the power required to accelerate the vehicle and the engine rotating parts, respectively (ω is the rotational velocity of the rotating parts, and I R the generalized moment of inertia).For each trajectory, v is known as a function of time (and ω is known, too, via the transmission ratio).Then, the acceleration is known, and the required engine power, .W engine , is obtained from Equation (1). Figure 1 shows the relevant parameters: velocity, acceleration, required engine power, and covered distance as a function of time, for the three trajectories T1, T2, and T3 of scenario AC1.
acceleration first and slow acceleration later (T2), and slow acceleration first and fast acceleration later (T3).Each Bezier polynomial had 1000 points;

•
The transmission ratio was 23 km/h per 1000 rpm.
The idealized dynamic equation of the vehicle in this case (neglecting friction and rolling resistance) was: The term on the left,  , is the engine power.The first term on the right is the power needed to overcome the aerodynamic drag ( is the air density,  the vehicle velocity).The third and fourth terms represent the power required to accelerate the vehicle and the engine rotating parts, respectively ( is the rotational velocity of the rotating parts, and  the generalized moment of inertia).For each trajectory,  is known as a function of time (and  is known, too, via the transmission ratio).Then, the acceleration is known, and the required engine power,  , is obtained from Equation (1). Figure 1 shows the relevant parameters: velocity, acceleration, required engine power, and covered distance as a function of time, for the three trajectories T1, T2, and T3 of scenario AC1.The value of the parameters at both starting and end points of the trajectories are given in Table 1 (they are the same for all trajectories within the acceleration scenario AC1).The value of the parameters at both starting and end points of the trajectories are given in Table 1 (they are the same for all trajectories within the acceleration scenario AC1).The outlook of scenario 2 is similar to that of scenario 1, but for three differences: • The car accelerated from 3 m/s (11 km/h) to 14 m/s (50 km/h) along an inclined plane (5% slope) 250 m long, see Figure 2; • The transmission ratio was 10.8 km/h per 1000 rpm; • An additional term was included in the right-hand side of Equation ( 1) to account for the gain of potential energy per unit time when climbing the slope: .
where g, m, and z are gravity acceleration, vehicle mass, and height.Figure 2 presents the relevant trajectory parameters of acceleration scenario 2 (AC2).The outlook of scenario 2 is similar to that of scenario 1, but for three differences: • The car accelerated from 3 m/s (11 km/h) to 14 m/s (50 km/h) along an inclined plane (5% slope) 250 m long, see Figure 2; The transmission ratio was 10.8 km/h per 1000 rpm; • An additional term was included in the right-hand side of Equation ( 1) to account for the gain of potential energy per unit time when climbing the slope: where g, m, and z are gravity acceleration, vehicle mass, and height.Figure 2 presents the relevant trajectory parameters of acceleration scenario 2 (AC2).Parameters at the beginning and end of the trajectories are presented in Table 2.  Parameters at the beginning and end of the trajectories are presented in Table 2.

Engine Model
The engine model and its validation are presented in Appendix A and Appendix B, respectively.In this subsection, only its main characteristics are outlined.

Engine Model
The engine model and its validation are presented in Appendices A and B, respectively.In this subsection, only its main characteristics are outlined.Vehicles 2024, 6 991

Tensor Mirror Model
Five 3D data tensors were generated using the engine model.The axes of the tensors represented the input parameters.The tensors themselves stored the output parameters.
Output parameters were engine power, fuel consumption, CO 2 and CO emissions, and risk of engine knocking.The total number of computed cases was 15 × 13 × 14 = 2730.The CPU time of each case was 3.4 min in a desktop computer (most of the computational load was associated to the chemical kinetics subroutine).Total CPU time for the 2730 cases was 155 h.

Densified Surrogate Tensor Model
The discretization described in Section 3.3 was coarse and, therefore, not practical for direct use on the solution of the problem under consideration.The reason for the coarseness was the high computational cost.The approach that was adopted to circumvent the difficulty was to perform a higher-order singular value decomposition (HOSVD) of the five computed tensors, followed by a densification of the tensors themselves based on cubic spline interpolation in the HOSVD modes.
HOSVD is a tensor decomposition technique developed by De Lathauwer et al. [25] and has been used extensively in many technical fields.For example, a practical case where HOSVD was used to densify an experimental database in which experiments were very expensive to perform was presented by Garcia-Magarino et al. [26].In the present study, the use of HOSVD presents two main advantages.The first is that densification by using this technique is not of a local nature.Instead, is it performed by interpolation on the HOSVD modes.This means that the process of generating new data is carried out according the global patterns (the HOSVD modes) of the database.The second is that if the optimization algorithm is of an algebraic nature (as in the present work), HOSVD preserves a tensor structure that facilitates its use in this algebraic context.

Generic Description of the Optimization Process
The velocity/acceleration paths defined in Sections 2.3.1 and 2.3.2,(see top subplots of Figures 1 and 2) yield the evolution of the vehicle RPM since the transmission ratio was prescribed for each acceleration scenario.Accordingly, the optimization path follows the RPM axis of the data tensor (see Figure 4).At each PA-IT plane of constant RPM it searches for the point at the next constant RPM plane that minimizes a given objective function.In Figure 4, point Q1 represents the vehicle operational point with data tensor coordinates (PA1, IT1, RPM1).Since the acceleration path is prescribed, the RPM at the next time step, RPM2, is known.Next point Q2 along the optimization path is located at horizontal plane RPM2, and its coordinates PA2 and IT2 are obtained via minimization of an objective function to be defined below.Then, the process goes to the next point Q3 and so on, until the end of the trajectory is reached.

Objective Function
The dimensionless objective function, Φ, to be minimized was defined as the summation of six different terms: The first term Φ is related to the prescribed engine power along the acceleration path (see bottom-left subplots of Figures 1 and 2).The definition of Φ is as follows: where  is the Matlab operator that rounds a certain number to the nearest integer, and  is the demanded engine power according to the corresponding acceleration paths (see Figures 1 and 2).Formula (4) penalizes values of engine power that deviate more than ±1% from the demanded power.For example, if  = 0.980  , then Φ = 100  (0.490 + |(0.98 − 1) 1 ⁄ |) = 100.On the other hand, if  = 0.995  then Φ = 100  (0.490 + |(0.995 − 1) 1 ⁄ |) = 0 .Formally, the engine power should be a restriction of the problem.However, in the present approach, it has been preferred to treat it as a separate term in the objective function.The second term, Φ , penalizes those data points whose risk of engine knocking, , is larger than 1: Third, fourth, and fifth terms evaluate the importance of fuel consumption, C (kg/s), and emissions of CO2 and CO in kg/s: where  ,  , and  are weighting factors. , 2 , and  are reference values that correspond to the end of the acceleration paths (i.e.: where the velocity is maximum and acceleration is zero, see Figures 1 and 2).
The sixth term, Φ , is related to the fact the engine operation has to be smooth.Then, the Euclidean distance inside the data tensor between one operational point and the next one that is being evaluated, see Figure 3, has to be as small as feasible.Thereby, Φ penalizes large distances inside the data tensors as follows:

Objective Function
The dimensionless objective function, Φ, to be minimized was defined as the summation of six different terms: The first term Φ 1 is related to the prescribed engine power along the acceleration path (see bottom-left subplots of Figures 1 and 2).The definition of Φ 1 is as follows: where round is the Matlab operator that rounds a certain number to the nearest integer, and P d is the demanded engine power according to the corresponding acceleration paths (see Figures 1 and 2).Formula (4) penalizes values of engine power that deviate more than ±1% from the demanded power.For example, if P = 0.980 P d , then Φ 1 = 100 round (0.490 + |(0.98 − 1)/1|) = 100.On the other hand, if P = 0.995 P d then Φ 1 = 100 round (0.490 + |(0.995 − 1)/1|) = 0. Formally, the engine power should be a restriction of the problem.However, in the present approach, it has been preferred to treat it as a separate term in the objective function.The second term, Φ 2 , penalizes those data points whose risk of engine knocking, Dd, is larger than 1: Third, fourth, and fifth terms evaluate the importance of fuel consumption, C (kg/s), and emissions of CO 2 and CO in kg/s: where λ 3 , λ 4 , and λ 5 are weighting factors.C re f , CO2 re f , and CO re f are reference values that correspond to the end of the acceleration paths (i.e.: where the velocity is maximum and acceleration is zero, see Figures 1 and 2).
The sixth term, Φ 6 , is related to the fact the engine operation has to be smooth.Then, the Euclidean distance inside the data tensor between one operational point and the next one that is being evaluated, see Figure 3, has to be as small as feasible.Thereby, Φ 6 penalizes large distances inside the data tensors as follows: where D is the Euclidean distance between two operational points inside the data tensor and D re f is the maximum distance between two points located at two constant RPM planes.For example, consider two points, 1 and 2, whose coordinates inside the data tensor are (PA, IT, RPM) = (0.4 bar, 4.5 deg, 3181 RPM), and (PA, IT, RPM) = (0.665 bar, −19 deg, 3225 RPM), respectively; these coordinates, in terms of positions along the tensor axes, are (21,60,256) and (74,13,261), respectively (recall the discretization of PA, IT, and RPM described at the beginning of Section 3.1).Then, the distance, D, between points 1 and 2 is Then, in this case, D/D re f = 71.01/184.46= 0.385.In the baseline optimization, the weight factors λ 1 , λ 2 , and λ 3 had the values 1, 1, and 1, respectively.

Detailed Description of the Optimization Process
The steps of the optimization process were as follows: 1.
The starting engine operating point, P0, in the (PA, IT, RPM) data tensor corresponds to the initial point of the acceleration path.This point has constant vehicle velocity and engine RPM's, zero acceleration, minimum fuel consumption, C, and minimum CO 2 and CO emissions.It provides the power needed to propel the vehicle at the prescribed constant velocity, and its risk of engine knocking is smaller than 1.

2.
The next point, P1, is selected along the corresponding acceleration path.In practical terms, this point, that is the next Bezier point that defines the path, prescribes vehicle velocity, engine power, and engine RPM (see Figures 1 and 2).

3.
Once the variable RPM is prescribed (it is called RPM1 for convenience), the objective function Φ, as defined in Equation ( 2), is computed for all points located in the plane of constant RPM1 of the (PA, IT, RPM) data tensor.The number of points where evaluation of Φ takes place is 141 × 121 = 17,061.

4.
Out from the 17,061 evaluations of Φ, the one that delivers the minimum value, Φ min , is selected.This point where Φ = Φ min corresponds to an associated pair of PA and IT (denoted with the subscript "1" for convenience).Then, the next optimum point in the data tensor of operation variables is (PA1, IT1, RPM1).

5.
The process is repeated iteratively until the end point of the acceleration path is reached.6.
Finally, integrals of the time evolution of fuel consumption, C, and CO 2 and CO emissions along the optimum path inside the data tensor are computed to obtain the total values of fuel consumed and pollutant emissions: C total , CO2 total , and CO total , respectively.

Acceleration Scenario AC1
Cumulative fuel consumption, C, and CO 2 and CO emissions along the three trajectories in scenario AC1 are presented in Figure 5.The integral values (total values) at the end of each trajectory, needed for comparative evaluation of the trajectories, are presented in Table 3.The integral values (total values) at the end of each trajectory, needed for comparative evaluation of the trajectories, are presented in Table 3.The following aspects could be observed: • The integrals of the CO 2 emissions alongside the 1 km long trajectories, T1, T2, T3, were similar.They were all in a band of 0.232 kg ± 6%.The integral of CO emissions varied between 0.631 × 10 −3 kg and 1.646 × 10 −3 kg.The ratio of the integral of CO 2 emissions to the integral of CO emissions had values between 362 and 134.• Total CO emissions in trajectory T3 nearly tripled those of trajectories T1 and T2.The reason is that the initial point of the trajectories had a high CO emission (see Table 1).Then, while T1 and T2 trajectories rapidly increase their velocity and acceleration, thereby moving away from the high CO emissions operating region, trajectory T3 keeps close to the initial operation point for a significant part of the trajectory.Thus, the integral of the CO emissions (that is the variable represented in Figure 5) is larger than its counterparts of trajectories T1 and T2.Then, T3 is discarded as a practical option when compared to trajectories T1 and T2.• Trajectories T1 and T2 present a similar lever of total CO emissions.Fuel consumption and CO 2 emissions of trajectory T2 are 8% larger than in trajectory T1.However, trajectory T2 saves time (by a factor of 11%) as compared to trajectory T1. • All this suggests that trajectories T1 and T2 are preferable in acceleration scenario AC1.Specific selection of each of them could be based on either minimizing emissions (T1) or reducing travel time (T2).In any case, all acceleration profiles that range from nearly constant acceleration to fast acceleration first and slow acceleration later might be acceptable in a practical driving situation.This provides a reasonable margin to the designer of intelligent driving approaches when these criteria are to be combined with other concerns related to the interaction of the vehicle with the surrounding traffic flow.
A 3D view of the evolution of the control parameters, PA, IT, and RPM for the three trajectories T1, T2, and T3, is presented in Figure 6.
The following aspects could be observed:

•
The integrals of the CO2 emissions alongside the 1 km long trajectories, T1, were similar.They were all in a band of 0.232 kg ± 6%.The integral of CO em varied between 0.631 × 10 −3 kg and 1.646 × 10 −3 kg.The ratio of the integral emissions to the integral of CO emissions had values between 362 and 134.

•
Total CO emissions in trajectory T3 nearly tripled those of trajectories T1 and T reason is that the initial point of the trajectories had a high CO emission (see Ta Then, while T1 and T2 trajectories rapidly increase their velocity and accele thereby moving away from the high CO emissions operating region, traject keeps close to the initial operation point for a significant part of the trajectory the integral of the CO emissions (that is the variable represented in Figure 5) is than its counterparts of trajectories T1 and T2.Then, T3 is discarded as a pr option when compared to trajectories T1 and T2.

•
Trajectories T1 and T2 present a similar lever of total CO emissions.Fuel con tion and CO2 emissions of trajectory T2 are 8% larger than in trajectory T1.Ho trajectory T2 saves time (by a factor of 11%) as compared to trajectory T1.

•
All this suggests that trajectories T1 and T2 are preferable in acceleration sc AC1.Specific selection of each of them could be based on either minimizing sions (T1) or reducing travel time (T2).In any case, all acceleration profiles tha from nearly constant acceleration to fast acceleration first and slow acceleratio might be acceptable in a practical driving situation.This provides a reasonabl gin to the designer of intelligent driving approaches when these criteria are combined with other concerns related to the interaction of the vehicle with th rounding traffic flow.
A 3D view of the evolution of the control parameters, PA, IT, and RPM for th trajectories T1, T2, and T3, is presented in Figure 6.Since in this study the air to fuel mass ratio is considered to be stoichiometr emissions depend directly on the gas temperature.During the time span between t of combustion and the beginning of exhaust in the expansion process (tce), chem netics (that depends strongly of temperature) still influence the chemical composi the gas.This exhaust temperature, in turn, depends on both ignition timing (IT) a gine load (inlet pressure, PA).Positive values of IT (ignition occurs before top dead c lead to lower exhaust temperatures because the time span, tce, in which the gas ex Since in this study the air to fuel mass ratio is considered to be stoichiometric, CO emissions depend directly on the gas temperature.During the time span between the end of combustion and the beginning of exhaust in the expansion process (tce), chemical kinetics (that depends strongly of temperature) still influence the chemical composition Vehicles 2024, 6 996 of the gas.This exhaust temperature, in turn, depends on both ignition timing (IT) and engine load (inlet pressure, PA).Positive values of IT (ignition occurs before top dead center) lead to lower exhaust temperatures because the time span, tce, in which the gas expands after combustion is larger.On the contrary, negative values of IT (ignition occurs after top dead center) lead to higher exhaust temperatures because the time span "tce" in which the gas expands is smaller.That is, negative values of IT lead to larger CO emissions (higher temperatures during expansion).This is observed in Figure 7, where the instantaneous CO emission (not the integral) is plotted as function of distance (top right), together with the exhaust temperature (bottom right) for the three trajectories T1, T2, and T3.In these two subplots, the direct relation between CO emissions and exhaust temperature is apparent.Subplot bottom-left of Figure 7 presents IT as a function of distance along the three trajectories.T3 (solid blue line) has negative values of IT for the first 300 m of the trajectory (approximately).This leads to higher gas temperatures during the "tce" time span and, therefore, to higher CO emissions as compared to the T1 and T2 trajectories.At the end of the 1000 m, IT is nearly the same for all three trajectories, but T3 shows a peak on engine load (PA), top left subplot, caused by the late power demand of this trajectory (see Figure 1).This increase in engine load leads, also, to higher CO emissions.
Vehicles 2024, 6 after combustion is larger.On the contrary, negative values of IT (ignition occurs a dead center) lead to higher exhaust temperatures because the time span "tce" in the gas expands is smaller.That is, negative values of IT lead to larger CO em (higher temperatures during expansion).This is observed in Figure 7, where the in neous CO emission (not the integral) is plotted as function of distance (top right), to with the exhaust temperature (bottom right) for the three trajectories T1, T2, and these two subplots, the direct relation between CO emissions and exhaust temper apparent.Subplot bottom-left of Figure 7 presents IT as a function of distance alo three trajectories.T3 (solid blue line) has negative values of IT for the first 300 m trajectory (approximately).This leads to higher gas temperatures during the "tc span and, therefore, to higher CO emissions as compared to the T1 and T2 trajecto the end of the 1000 m, IT is nearly the same for all three trajectories, but T3 shows on engine load (PA), top left subplot, caused by the late power demand of this tra (see Figure 1).This increase in engine load leads, also, to higher CO emissions.

Acceleration Scenario AC2
Cumulative fuel consumption, C, and CO2 and CO emissions along the three tories in scenario AC2 are presented in Figure 8.The integral values at the end of each trajectory are presented in Table 4.The integral values at the end of each trajectory are presented in Table 4. Broadly speaking, the results obtained go in line with those obtained for acceleration scenario AC1, meaning that trajectory T3 should be discarded because of its larger CO emissions.However, regarding trajectories T1 and T2, the latter (T2) appears to be preferable.It is true that its fuel consumption and CO 2 emissions are 5% larger than for T1, but CO emissions are 25% lower, and the travel time is 25% smaller as well.The evolution of the control parameters, PA, IT, and RPM is presented in Figure 9.
Vehicles 2024, 6 Broadly speaking, the results obtained go in line with those obtained for accel scenario AC1, meaning that trajectory T3 should be discarded because of its lar emissions.However, regarding trajectories T1 and T2, the latter (T2) appears to be able.It is true that its fuel consumption and CO2 emissions are 5% larger than for CO emissions are 25% lower, and the travel time is 25% smaller as well.The evolu the control parameters, PA, IT, and RPM is presented in Figure 9.

Conclusions
A conceptual study has been carried out regarding acceleration strategies for gent driving, with the objective being qualitative identification of optimized traje in view of minimizing fuel consumption, CO and CO2 emissions, and travel tim acceleration scenarios were selected that represent typical driving situations.The volved acceleration along a horizontal plane using a high transmission ratio (h driving).In the second, the vehicle accelerated along a short, inclined plane (5% with a low transmission ratio (urban driving).For each acceleration scenario, three ent trajectories were considered: (a) nearly constant acceleration (T1), (b) fast accel first and slow acceleration later (T2), and (c) slow acceleration first and fast accel later (T3).Engine parameters were typical of a generic European small-medium ger car (90 kW and 1150 kg).The engine model had two inputs (inlet pressure and i timing) and three outputs (fuel consumption, and CO2 and CO emissions).Trav was computed from the definition of the trajectories.Within each trajectory, time tion of the inputs was optimized so as to minimize the outputs (one restriction w

Conclusions
A conceptual study has been carried out regarding acceleration strategies for intelligent driving, with the objective being qualitative identification of optimized trajectories in view of minimizing fuel consumption, CO and CO 2 emissions, and travel time.Two acceleration scenarios were selected that represent typical driving situations.The first involved acceleration along a horizontal plane using a high transmission ratio (highway driving).In the second, the vehicle accelerated along a short, inclined plane (5% slope) with a low transmission ratio (urban driving).For each acceleration scenario, three different trajectories were considered: (a) nearly constant acceleration (T1), (b) fast acceleration first and slow acceleration later (T2), and (c) slow acceleration first and fast acceleration later (T3).Engine parameters were typical of a generic European small-medium passenger car (90 kW and 1150 kg).The engine model had two inputs (inlet pressure and ignition timing) and three outputs (fuel consumption, and CO 2 and CO emissions).Travel time was computed from the definition of the trajectories.Within each trajectory, time evolution of the inputs was optimized so as to minimize the outputs (one restriction was prescribed for risk of detonation effects).Then, optimized trajectories were compared among themselves.The optimization process involved a new methodological approach based on the higher-order singular value decomposition of the tensor form of the engine model.
The first conclusion (valid for both acceleration scenarios) is that trajectory T3 should be discarded because of the larger amount of CO emissions as compared to trajectories T1 and T2.In particular, T3 nearly tripled the emissions of trajectories T1 and T2 while having comparable CO 2 emissions.Then, whenever possible, this type of trajectory (slow acceleration first, fast acceleration later) should be avoided when programming intelligent driving strategies.
Trajectories T1 and T2 presented similar behavior regarding CO emissions in acceleration scenario 1 (they differ by a factor of 1% only).However, fuel consumption and CO 2 emissions of trajectory T2 are 8% larger while its travel time is smaller by a factor of 11%.Then, if fuel consumption and CO 2 emissions reduction are at a premium, trajectory T1 should be selected.If travel time reduction is preferred, the selection should be T2.
Regarding acceleration scenario 2, trajectory T2 is preferable over T1.Its fuel consumption and CO 2 emissions are larger by a factor of 5%, but both its CO emissions and travel time are smaller by a factor of 25%.
The selection of the trajectories becomes more apparent when results are analyzed in the context of the driving environment.For example, when driving along a highway, short acceleration intervals are interspersed between long intervals of constant velocity.In this situation, total travel time, say, from city A to city B, is, basically, associated with the long constant velocity intervals (not with the acceleration intervals), so the preferred trajectory for the acceleration phases should be T1.In the case of urban driving, short acceleration intervals occur very frequently, so the preferred trajectory should be T2.
The results obtained show that, in general, the class of trajectories situated in between T1 and T2 are preferable.This is important for the designers of intelligent driving systems.The reason is that there is reasonable margin in between T1 and T2 to accommodate further integration of optimized trajectories with algorithms that account for the state of traffic flow.
The global conclusion is that real-time adaptation of acceleration strategies to the driving environment may have a significant impact on the reduction in fuel consumption and emissions in internal combustion engines, without sacrificing travel time significantly.computed from the model (original tensor).Then, 600 cases times each of the five variables yields a total of 3000 ε values that are presented as a histogram in Figure A4.It could be observed that, for the vast majority of cases, discrepancy ε was smaller than 2 × 10 −5 , so tensor densification via HOSVD was appropriate for the problem under consideration.
where  _ ,  _ ,  _ , and  _ are variable "i" computed from the model, variable "i" interpolated vía HOSVD, and maximum and minimum values of variable "i" computed from the model (original tensor).Then, 600 cases times each of the five variables yields a total of 3000  values that are presented as a histogram in Figure A4.It could be observed that, for the vast majority of cases, discrepancy  was smaller than 2 × 10 −5 , so tensor densification via HOSVD was appropriate for the problem under consideration.

3 . Methodology 3 . 1 . 3 . 1 .
Overview of the Methodology The methodology has three main blocks: (a) engine model, (b) optimization algorithm, and (c) trajectories.The engine model has three components: (i) model equations, (ii) tensor mirror model, and (iii) densified surrogate tensor model.The optimization algorithm reads the model and the trajectories.It is supplied with the engine input parameters and it delivers the optimized engine output parameters.They, together with the trajectory's times, allow for a comparison of the different trajectories.An overview of the methodology in the form of a block diagram is presented in Figure 3. Overview of the Methodology The methodology has three main blocks: (a) engine model, (b) optimization algorithm, and (c) trajectories.The engine model has three components: (i) model equations, (ii) tensor mirror model, and (iii) densified surrogate tensor model.The optimization algorithm reads the model and the trajectories.It is supplied with the engine input parameters and it delivers the optimized engine output parameters.They, together with the trajectory's times, allow for a comparison of the different trajectories.An overview of the methodology in the form of a block diagram is presented in Figure 3.

Figure 3 .
Figure 3. Block diagram of the methodology.

•
Gasoline Otto cycle; • Zero-dimensional (time is the only independent variable); • Heat transfer to the outside is accounted for in all thermodynamic cycle processes; • Pressure losses through the intake and exhaust valves are modelled; • Combustion takes place with 10 species: CO2, H2O, OH, CO, H2, H, O2, O, N2, and N, considered to be chemical equilibrium; • Atmospheric engine (no compressor); • The model is made up of 35 integral-differential equations that involve all processes of the thermodynamic cycle.

Figure 3 .
Figure 3. Block diagram of the methodology.

•
Gasoline Otto cycle; • Zero-dimensional (time is the only independent variable); • Heat transfer to the outside is accounted for in all thermodynamic cycle processes; • Pressure losses through the intake and exhaust valves are modelled; • Combustion takes place with 10 species: CO 2 , H 2 O, OH, CO, H 2 , H, O 2 , O, N 2 , and N, considered to be chemical equilibrium; • Atmospheric engine (no compressor); • The model is made up of 35 integral-differential equations that involve all processes of the thermodynamic cycle.

Figure 4 .
Figure 4. Example of RPM plane inside the 3D tensor of a given output Q.

Figure 4 .
Figure 4. Example of RPM plane inside the 3D tensor of a given output Q.

Figure 5 .
Figure 5. Cumulative fuel consumption and emissions in scenario AC1.

Figure 5 .
Figure 5. Cumulative fuel consumption and emissions in scenario AC1.

Figure 7 .
Figure 7. Top left: engine load (inlet pressure) as a function of distance for trajectories T1 (bl (red), and T3 (blue).Top right: instantaneous CO emissions.Bottom left: ignition timing IT. right: exhaust gas temperature.

Figure 8 .
Figure 8. Cumulative fuel consumption and emissions in scenario AC2.

Figure 9 .
Figure 9. Cumulative fuel consumption and emissions in scenario AC2.

Figure 9 .
Figure 9. Cumulative fuel consumption and emissions in scenario AC2.

Table 1 .
Vehicle velocity, RPM, engine power, distance covered, fuel consumption, and pollutants emissions at the starting and end points of trajectories in acceleration scenario AC1.

Table 1 .
Vehicle velocity, RPM, engine power, distance covered, fuel consumption, and pollutants emissions at the starting and end points of trajectories in acceleration scenario AC1.

Table 2 .
Vehicle velocity, RPM, engine power, distance covered, fuel consumption, and pollutants emissions at the starting and end point of trajectories in acceleration scenario AC2.

Table 2 .
Vehicle velocity, RPM, engine power, distance covered, fuel consumption, and pollutants emissions at the starting and end point of trajectories in acceleration scenario AC2.

Table 3 .
Integral values (total values) at the end of each trajectory of fuel consumption, C, emissions, CO2 and CO, and trajectory time, in the acceleration scenario AC1.

Table 3 .
Integral values (total values) at the end of each trajectory of fuel consumption, C, emissions, CO 2 and CO, and trajectory time, in the acceleration scenario AC1.

Table 4 .
Integral values (total values) at the end of each trajectory of fuel consumption, C, emissions, CO 2 and CO, and trajectory time, in the acceleration scenario AC2.