Ecodriving for Reduction of Bus Transit Emission with Vehicle ’ s Hybrid Dynamic Model

This paper formulates a global ecodriving optimal control to advise the green driving speed for bus transit to minimize the exhaust emission using Vehicle-to-Infrastructure (V2I) communication. Assuming communication between vehicles and infrastructure (V2I) and knowledge of traffic signal timings and waiting passengers at stations are known, an optimal driving speed is proposed to minimize the total vehicle emissions of the bus route.The dwell time of the bus transit at each station which includes two parts is proposed. A traffic lights timing model is employed as constraints to control the formation of the green wave band. Vehicle specific power (VSP)model is further applied to evaluate the exhaust emission level linkedwith the speed and acceleration of the bus transit. An approximate sixteen-kilometer traffic network including fourteen intersections and fifteen stations of Beijing bus transit line 1 in Chaoyang District, Beijing, is chosen to investigate the performance of the developed ecodriving approach.


Introduction
Traffic emissions are among the biggest sources of environmental pollution in cities.The emissions of vehicles contain several harmful substances, such as nitrogen oxides (NO  , i.e., nitrogen monoxide and nitrogen dioxide), hydrocarbons (HC), carbon monoxide (CO), carbon dioxide (CO 2 ), and fine particulate matter.Around the world, approximately 50% of the NO  emission and 90% of the CO emission come from traffic [1].Consequently, traffic emissions, as an environmentoriented traffic control objective, draw more attention [2][3][4][5][6][7].
Ecodriving in urban traffic networks using traffic signal information is a significant field of intelligent transportation system (ITS).Technology advancement in wireless communication, specifically Vehicle-to-Infrastructure (V2I) and Vehicle-to-Vehicle (V2V) communications, has the potential of data collection (road, vehicle, and passenger) at any space-time position of traffic networks.Knowledge of traffic signal timings has been proven to benefit the energy use of vehicles by reducing the sudden break and idling time at the intersections [3,8].
Hybrid systems have been an active research topic in the last few years and many results can be found in the control and computer science literature [9].As a special type of hybrid systems, switched systems have attracted increasing attention [10,11].Wei et al. study the problem of traction control, that is, how to stabilize a wheeled mobile robot (WMR) subject to wheel slippage to a desired configuration [12].Uthaichana et al. present a hybrid optimal control solution for the power management problem of a parallel hybrid electric vehicle (PHEV), based on a two-mode operation model of PHEV [13].
The green light optimal speed advisory (GLOSA) [14] can synchronize with green signals without stopping at intersections when vehicles are traveling at an optimal speed.In addition, GLOSA helps vehicles keep a driving speed as much as possible to reduce additional acceleration.In [15], authors provided a fast suboptimal solution that suggests an online speed advisory service to the driver for the energy minimization based on GLOSA system.A framework for vehicle's ecodriving was developed in [16] based on the V2I communication technology to enhance vehicle fuel consumption efficiency.In [17], Asadi and Vahidi propose the use of upcoming traffic signal information within the vehicle's adaptive driving control system to reduce stop time at stop lights and fuel consumption.To achieve this goal an optimization-based control algorithm is formulated that uses short range radar and traffic signal information to schedule an optimum velocity trajectory for the vehicle.More extensive discussion about how to find the energy-optimal velocity for a vehicle using V2I and V2V communication in urban traffic networks can be reviewed in the literature [18,19].
Here we are interested in optimal control of bus transit with minimizing the volume of pollutant of vehicle emissions, such as NO  , HC, and CO 2 .The vehicle emission model, vehicle specific power (VSP) model [20,21], is used as performance index to quantify the vehicle emissions, which is closely related to the vehicle dynamics (speed, acceleration) and the road grade.The referred trajectories induced from the traditional "bus priority" and "green wave" schemes of traffic signal timing are tracked by the driver in order to pass through the intersections without stopping.Using V2I communication technology, not only the information of traffic signal planning but also the estimated total delay at each station is transmitted to drivers for decision making.Physical topology graph is shown in Figure 1.Following from our discussion, we also obtain a suboptimal ecodriving trajectory for bus transit to pass through intersections in green light against the dwell time of bus transit at stations.
This paper is organized as follows: Section 2 modeling the vehicle's hybrid dynamic model with switching laws and guard conditions.The traffic signal timing scheme and the global optimal control model for ecodriving are introduced in Sections 3 and 4, respectively.Section 5 describes the iterative optimization method.In Section 6, we select an approximate sixteen-kilometer traffic route of Beijing bus transit line 1 for simulation to verify the validity of the vehicle's hybrid dynamic model through the numerical solution of the simulation.Conclusions and future works are development in Section 7.

Vehicle's Hybrid Dynamic Model
The movement of the vehicle is discontinuous, which is interrupted by the passengers regularly getting on and off at the stations by the bus transportation schedule.At the same time, the different vehicle dynamics, such as driving and dwell process, usually lead to the different emission levels of the pollutant.Therefore, it is convenient to build the movement of the vehicle with the hybrid dynamical model.Vehicle's hybrid dynamic model and switching laws among hybrid dynamics are formulated as follows.
Driving Mode.Using the vehicle kinematics, the dynamic of driving mode can be stated as the following equations: where () is the displacement of vehicle, V() is the speed of vehicle,  is the vehicle mass,   () is traction force of the engine, and   () means the road forces including aerodynamic drag, rolling resistance, and road grade forces [17].
The road forces   () can be described in detail as where  is the gravitational constant,  is the road grade angle,   is the rolling resistance coefficient,   is the mass density of air,  is the cross-sectional area, and   is the aerodynamic drag coefficient.
When the vehicle dynamics are in the driving mode the vehicle accelerates leaving the station, goes through the intersection with green light, and decelerates into the station.During the driving mode process the vehicle is driving at an expected speed if we neglect the transient acceleration and deceleration process.
Dwell Mode.The location of each station is known, the vehicle dynamics are in the dwell mode when () =   , where   is the location of station   .The dynamic of dwell mode can be stated as the following equations: where () is a timer used to record the vehicle's dwell time.When the vehicle leaves the station, the timer is cleared.
The dwell time of vehicle    at station   includes two parts; one is dwell time of passengers regularly getting on and off    (); the other one is control quantity   .  is controlled by the driver, meaning the dwell time after all the passengers have regularly finished getting on and off.
The dwell time of passengers regularly getting on and off    () can be estimated with the number of waiting passengers    () at station [23], in which where    () is the number of waiting passengers and (⋅) is a linear or nonlinear function about the number of passengers.
is a time independent control variable of vehicle's hybrid dynamic model.The upper bound value of   is   max , so the dwell time of vehicle    can be written as the following equation: where Dwell mode Driving mode Switching Laws.The vehicle's hybrid dynamic model contains two kinds of switching processes.The one switching process is switching from driving mode to dwell mode; it is autonomous switches; the other switching process is switching from dwell mode to driving mode; it is controlled switches.Figure 2 describes the switching laws and the guard conditions for the hybrid dynamic model of vehicle.The guard condition from the driving mode to the dwell mode is autonomous as the vehicle arriving stations, while the guard condition from the dwell mode to the driving mode is modeled as control variables of hybrid model such that the timer () is equal to the dwell time of vehicle    .

Traffic Signal Timing
Suppose there is an intersection between every two adjacent stations.As it is shown in Figure 3, vehicle leaves from any station   driving to the downstream adjacent station  +1 and goes through the intersection   with green light.Going through the intersection with green light can reduce vehicle acceleration and deceleration and vehicle brake numbers as well as vehicle delay and emissions.Therefore it is necessary to know the traffic signal timing of intersection in order to ensure the vehicle without stop at intersection.In our research, the traffic signal timing at intersections including cycle, split, and offset times is deterministic and given.Under the assumption that the information of signal timing of the traffic lights can be shared with the vehicle via the V2I communication, it is possible to formulate mathematically the state evolution of the traffic lights as follows: where   is the state of traffic light (red or green) at intersection   ,    is the time instant at which the vehicle arrives at the intersection   ,  is the cycle,    and   are the start instant and duration of the green phase, respectively, and  = 1, 2, . . . is the number of cycles.The offset of the neighboring traffic lights is set as where V is the "green wave" speed,   is the position of intersection   , and  +1 is the position of upstream intersection  +1 .The "green wave" timing scheme of adjacent traffic lights makes the vehicle pass through the next intersection without stop, which has been extensively studied in the field of traffic engineering to reduce delay.We restrict the vehicle to pass through the intersections in the driving mode of hybrid model without stop, for intersection   , satisfying where "∧" is on behalf of the logic operation "and."

Global Optimal Problem
Objective Function.Take the minimum of emissions (NO  , HC, CO, and CO 2 ) as the target.The vehicle specific power (VSP) model [20] is used as the performance index to measure the pollutants of vehicle emissions.VSP is a proxy variable for engine load and is defined as the source power per unit mass.The calculation formula of the VSP is as follows: where V is vehicle speed,  is vehicle acceleration,   is vehicle mass factor,  is acceleration of gravity,  is vertical rise/slope length,   is coefficient of rolling resistance,   is ambient air density,   is drag coefficient,  is frontal area of the vehicle,  is vehicle mass, and V  is wind speed.Actually the calculated value of VSP through ( 9) is dimensionless and has not any actual physical meaning.Table 1 is the VSP modes and the average modal emission rates of each.By querying Table 1, we can obtain the emissions.Emissions can be expressed as where  is the emissions and (⋅) is a piecewise function.

Constraint
residence times (4) , (5) signal timing (6) , ( 7) , ( 8) , where    =   is the time to arrive at station   .The constraint condition (8) for the "green wave" of the signal timing represents a set of disjointed time intervals at the same intersection.Therefore, the constraint conditions are nonconvex sets and the solution of the optimal ecodriving of vehicle has to be considered as a suboptimal way.

Iterative Optimization Method
Dwell time of passengers getting on and off is a disturbance to the global optimal control.It is difficult to predict the values of    () in advance for each station.In order to solve this problem we should update the dwell time    () through V2I communication at the moment when the vehicle arrives at each station.
Suppose there are  stations and  − 1 intersections.The space discretization is carried out according to the position of the station.All optimal solutions can be solved in the iterative optimization; for each step, the vehicle reaches to the next downstream station.Every step there are   stations and   inputs.  means how "far" we wish the future to be predicted for, where   ≤  − 1 is positive integer.  dictates the number of parameters used to capture the future control trajectory,   ≤   .There are  −   + 1 iterative optimization steps.   and    are input vectors, where   () is the total emission of th step iterative optimization and    and    are the emissions during the driving mode and dwell mode, respectively.
The performance indexes    and    of the vehicle's hybrid dynamics can be described as the following equations, respectively, with more details, for  = 1, 2, . . .,  − 1: where    and  +1  are the time instant at which the vehicle arrives at the stations   ,  +1 , respectively.
For th step the dimension of each input vector is   .We only implement the first element of these input vectors, that is,  1  () and  1  , while ignoring the rest of the elements (except the last step).When the next step arrives, the more recent measurement is taken from the state vector for calculation of the new sequence of input vector.This procedure is repeated in real time during the process of iterative optimization.The optimal emission for th step is  * (),  * () = (   +    ).Notice that when the vehicle arrives at station  −  +1 , that is, the last step, the number of stations at the downstream remains less than   .Then all elements of    The total optimal emissions  * total ,  * total = ∑ −  =1  * () +   ( −   + 1).
In the above iterative optimization method process, each step of the iterative optimization is an optimization problem.Main steps are including pruning algorithm, searching feasible path, and optimizing.More detailed steps are as follows.
Step 1. Determine the initial range of state variables.
Step 2. By pruning algorithm calculates feasible regions of each state variable.
Step 3. Determine feasible paths and calculate feasible regions of state variables for each path.
Step 4. Calculate feasible path's optimal solution through nonlinear programming.
Step 5. Repeat Step 4, until every feasible path has been calculated.
Step 6. Search for the optimal feasible path, that is, feasible path corresponding to the minimum value of the optimal index.
Remark 1.The classical model predictive control is a receding horizon optimization process.In our optimization method, the dynamic iterative process is receding along the spatial positions of the stations.

Experiments
A segment of bus transportation route of bus transit line 1 in Beijing is chosen as the simulation scene.This transportation route consists of fifteen stations, from the Yuquan Road station to the Ritan Road station, respectively, as seen in Figure 4, and fourteen traffic intersections with the twophase signal timing scheme.The spatial displacement of the stations and the traffic intersections are obtained from the Baidu Digital Map.
Bus transit line one is BJ6180C8CTD model of 18-meter LNG (liquefied natural gas) articulated vehicle made by FOTON motor incorporated company.The vehicle performance parameters are shown in Table 2.The unladen mass of vehicle is 17000 kilograms and the total mass is 28000 kilograms.Rated power  of the vehicle is 213 kilowatts.The physical dimension of the vehicle is 17.98-meter long, 2.55meter wide, and 3.1-meter high.See Table 2 for details.
Parameter values for calculating vehicle drag are shown in Table 3.The slope of the road is zero.The cross section of vehicle  is calculated by the height of the vehicle times the width of the vehicle.
The other model and simulation parameters, such as the station position, intersection position, speed limit of the vehicle, acceleration limit of vehicle, simulation times, and total number of stations are listed in Table 4.
Table 5 lists the parameters of traffic signal timing of the fourteen intersections.The cycle time is 120 seconds; the green time of the vehicle passing phase is 30 seconds.The offset of the neighboring intersections is 113-31-20-90-99-26-10-72-74-12-102-113-8 seconds, respectively.In this simulation, the start time  1  of the green phase of the first intersection is 0.
Simulation results of the total emissions  * total are 2420 grams.The total delay time at each station   , the maximum traction force    (i.e., the maximum value of each stage in Figure 7), the optimal speeds V * , the optimal acceleration  * , the optimal arriving time of each station    * , and the optimal arriving time of each intersection    * are shown in Table 6. Figure 5 gives the spatiotemporal diagram of optimal driving route for the public transits line one.Vehicle delays at   Figure 6 shows the variation curve of the optimal vehicle speed suggestion with time.The curve is periodic.In each period, the driving speeds are different.This is caused by the actual position of the stations and the total delay times.
Vehicle traction force of engine variation curves with time is shown in Figure 7.It is periodic too.In each period, the traction force of engine is zero at the station; it surges at the acceleration process and at last decreases to a constant value and maintains stability at driving mode.

Conclusions
According to the actual operation of the bus transit, vehicle's hybrid dynamic model is established.The dwell time of the vehicle at each station which includes two parts is proposed.Previous works have proved how the V2I has a positive impact on the energy consumption.In this work, further improvement of ecodriving to lower vehicle emission

Figure 4 :
Figure 4: Bus route of bus transit line one in Beijing city.

Table 1 :
[22]itions.Ensuring that the travel time vehicle must arrive at the terminal station   at the moment   ,    =   .The speed V  (), acceleration   (), and traction force of the engine    () values must be in the threshold range, V min ≤ V  () ≤ V max ,  min ≤   () ≤  max , and  min VSP modes and the average modal emission rates of each[22].

Table 3 :
Road forces calculation parameters.

Table 4 :
Model and simulation parameters.