Assessment of the Worthwhileness of E ﬃ cient Driving in Railway Systems with High-Receptivity Power Supplies

: Eco-driving is one of the most important strategies for signiﬁcantly reducing the energy consumption of railways with low investments. It consists of designing a way of driving a train to fulﬁl a target running time, consuming the minimum amount of energy. Most eco-driving energy savings come from the substitution of some braking periods with coasting periods. Nowadays, modern trains can use regenerative braking to recover the kinetic energy during deceleration phases. Therefore, if the receptivity of the railway system to regenerate energy is high, a question arises: is it worth designing eco-driving speed proﬁles? This paper assesses the energy beneﬁts that eco-driving can provide in di ﬀ erent scenarios to answer this question. Eco-driving is obtained by means of a multi-objective particle swarm optimization algorithm, combined with a detailed train simulator, to obtain realistic results. Eco-driving speed proﬁles are compared with a standard driving that performs the same running time. Real data from Spanish high-speed lines have been used to analyze the results in two case studies. Stretches fed by 1 × 25 kV and 2 × 25 kV AC power supply systems have been considered, as they present high receptivity to regenerate energy. Furthermore, the variations of the two most important factors that a ﬀ ect the regenerative energy usage have been studied: train motors e ﬃ ciency ratio and catenary resistance. Results indicate that the greater the catenary resistance, the more advantageous eco-driving is. Similarly, the lower the motor e ﬃ ciency, the greater the energy savings provided by e ﬃ cient driving. Despite the di ﬀ erences observed in energy savings, the main conclusion is that eco-driving always provides signiﬁcant energy savings, even in the case of the most receptive power supply network. Therefore, this paper has demonstrated that e ﬀ orts in improving regenerated energy usage must not neglect the role of eco-driving in railway e ﬃ ciency. driving strategies for the commercial running time in the Puertollano–Ciudad Real section. Results are shown for di ﬀ erent values of the engine e ﬃ ciency and the linear resistance of the catenary. and


Introduction
Railways are considered an inexpensive, fast and safe transport mode. Furthermore, this transport mode is energy-efficient and its contribution to global warming is less severe compared with others [1]. Despite that, its increasing activity throughout the world [2] makes it necessary to continue working on the energy efficiency of rail transport to achieve a sustainable transport industry without losing sight of service quality improvements.
Eco-driving, also named speed profile optimization or train trajectory optimization, is one of the most important measures with which to achieve significant energy reductions in the energy consumption

Train Simulation Model
In this work, an accurate simulation model is used to calculate the speed profiles of a train and its consumption measured at the pantograph and at substations. It is divided into three modules: train, line and manual driving [47]. The train simulator is a software developed by Comillas Pontifical University researchers. This simulator is a desktop application for Windows coded in an object-oriented environment using the C++/CLI programming language. A validation against real on-board registered measurements was performed for high-speed trains on commercial services and night tests on the Madrid-Barcelona line, demonstrating its accuracy. Simulation results taken from this software showed an average error of 1.2% in running time and 0.4% in energy consumption compared with on-board measurements [47]. This accurate simulator has been applied for the design of efficient driving and tested in collaboration with Renfe and Adif.
The train module includes mass, length, running resistance, adhesion coefficient and the engine characteristics. These are introduced through the maximum traction and braking effort curves, which depend on the train speed. In this work, the traction and braking efficiencies are taken as equal and non-dependent on the train speed. This module includes also the energy consumption of the on-board auxiliary systems and a minimum constant braking effort that is applied along neutral zones in order to feed the auxiliary systems when no braking curve is affecting the train.
The physical characteristics of the line-track geometry (grades, grades transitions and track curves and their transitions), speed limits, tunnels, stopping points, stations, neutral zones, electrical substations and catenary features-are introduced in the line module.
The manual driving module calculates the traction and braking force needed at every time step, taking into account track data, speed limits, the position and speed of the train and stopping points. The driving module takes as input driving commands that can be easily interpreted by a human driver and generates speed profiles. In this work, the human driving is emulated by using a PI (proportional-integral) controller, in which the engine effort is modelled as a variable that represents the proportion of traction or brake that is demanded with respect to the maximum values for the current train speed.

Train Dynamics Model
At every time step the train position s(t), speed v(t) and acceleration a(t) are determined by using Newton's second law Equations (1)- (5), where m is the train mass, ρ is the dimensionless rotating mass factor, F r (v) is the running resistance and F g (s) are the forces that affect the train due to the track geometry. F m (t) is the engine traction force and F b (t) is the braking force. t 0 , t f , s 0 , s f , v 0 and v f are the initial and final times, positions and speeds of the train, respectively. The engine traction force F m (t) is bounded by a maximum traction curve that depends on the speed of the train. The braking force F b (t) has two components: the electrical or regenerative and the Energies 2020, 13, 1836 5 of 24 pneumatic brakes. The regenerative brake is limited by a maximum electrical braking curve; if the train needs to brake it will use pneumatic brake to complement the electrical when the electrical one is saturated. The manual driving module computes both forces at every simulation time step to control the train dynamics considering the driving commands while obeying speed limits and stopping the train at the stopping points. In neutral zones, a minimum braking effort is applied to feed the auxiliary systems using regenerative energy.
The running resistance F r (v) depends quadratically on the train speed and it is modelled by the Davis formula (Equation (6)): F r (v) = A + Bv + T f Cv 2 (6) where A, B and C are the Davis formula coefficients. T f is a factor that models the impact of tunnels on the running resistance and it modifies the quadratic term. It is equal to one when the train is outside a tunnel and greater otherwise. The considered forces on the train due to the track geometry, F g (s), are those due to the track grades and curves: where m is the mass of the train; g is the gravity acceleration; p(s) and R(s) are the averages of the grades and curve radius respectively affecting the train in the position s. K is a constant that depends on the gauge of the track which units are the same as R(s) (meters). The gauge in the considered lines is standard UIC (Internal Union of Railways) gauge, 1.435 m, and K = 600. Notice that the force due to curves is modelled as an equivalent grade.

Energy Consumption Model
The simulation model calculates the power consumed at the pantograph and at the electrical substations at every time step as a function of the train speed and the engine traction or braking effort by using Equations (8) and (9) where η t is the engine traction efficiency, η b is the engine regenerative brake efficiency, and both are considered constant. P aux is the power consumed by the auxiliary systems; F max (v) is the maximum regenerative brake force at speed v. r(s) =r|s − s ss | is the resistance of the catenary that depends on the distance between the train and the electrical substation, which is located at position s ss . In this article the influence of the lineal resistancer on the energy consumption is studied. V is the nominal line voltage and cosϕ is the power factor, and they are assumed to be constant. If a train is in a neutral zone, the power consumed at the electrical substations is zero. The energy consumption at the pantograph and substations can be obtained by integrating Equations (8) and (9) in time, respectively. This section may be divided by subheadings. It should provide a concise and precise description of the experimental results and their interpretation, and the experimental conclusions that can be drawn.

Driving Strategies and Commands
In this work the performances of two driving strategies are compared in terms of energy savings. In particular, a holding speed without braking and final coasting phase eco-driving strategy is compared against a standard driving strategy consisting in holding a constant speed with braking. All the regenerated energy at the pantograph is assumed to be returned to the AC power grid. The dependence Energies 2020, 13, 1836 6 of 24 of the energy consumption measured at the pantograph and electrical substations on the engine traction and brake efficiency and on the losses at the catenary is analyzed.
For the eco-driving strategy, the journey is divided in two phases. In the first phase, the driver must hold a certain speed by applying traction. If braking is necessary to maintain the constant speed, coasting is applied instead, and thus, the speed of the train will increase. In the second phase, the driver must apply coasting until the train reaches its braking curve to stop at the station.
For the standard strategy, a cruise speed is determined for the whole journey. The driver has to maintain that speed until reaching the braking curve. To drive at a constant speed, braking can be applied if necessary. Any driving strategy has to be described in easily-interpretable commands so that a driver can apply them. The driving commands are defined by using a command vector. It has two components. The first one contains the position s c that corresponds to the final point of the holding speed phase. The holding speed v c is introduced in the second component. The final coasting phase is performed starting at s c up to the braking curve in the holding speed without braking and the final coasting driving strategy. Notice that speed limits of the line are observed by the driving model braking the train when necessary.

Driving Model
A closed-loop PI controller is used to emulate the manual driving by controlling the train dynamics. It takes into account the maximum acceleration and service deceleration of the train, the jerk limitation, the maximum tractive effort, the speed limits and the energy that has to be regenerated by braking in neutral zones. The control loop minimizes the difference between an objective speed v obj and the train speed v. At every time step the objective speed is computed in terms of the driving commands and the track speed limits. To do so, the braking curves to the next ceiling speed reduction, v b (s), and to the next stop at a station, v f b (s), are calculated by using Equations (10) and (11).
where v next is the next ceiling speed of the track, d service is the service deceleration of the train, s is the position of the train, s next is the position of the next ceiling speed corrected with the train length if necessary and s end is the next stopping point. The objective speed v obj is therefore computed in terms of the maximum admissible speed v max , which is the minimum between the next braking curve and the ceiling speed at the position of the train v ceiling (s) The control variable U reg is defined as the traction or brake demand of the PI controller. To calculate it, at every time step the proportional and integral control components are calculated as described in Equations (14)- (18). For the standard driving the control variable is computed using Equation (17), while for the eco-driving strategy Equation (18) is used instead.
U P (t) = K P ·e(t) where e(t) is the difference between the objective and the train speeds, whose absolute value has to be minimized by the control. U P (t) and U I (t) are the proportional and integral components of the control variable U reg (t). K P and K I are the proportional and the integral gains, respectively. ∆t is the simulation time-step. If the control variable is positive, the control demands tractive effort, while if it is negative, it demands braking. If it is zero, the control demands coasting. In the eco-driving strategy (Equation (18)), the control demand is set to coasting when the train speed is below the maximum admissible speed and the control would demand braking under the current condition.
Once that the control variable has been obtained, the engine maximum traction and brake limits; comfort constraints-maximum acceleration and jerk; and the presence of neutral zones, have to be taken into account. Strong accelerations are considered to be unpleasant for passengers, so the maximum acceleration and deceleration have to be controlled. F maxAcc (s, v) and F maxDec (s, v) model the forces corresponding to the maximum admissible acceleration and deceleration respectively with corrections for the running and track geometry resistances. Equations (19) and (20) are used for calculating these forces.
where a max and d max are the maximum admissible acceleration and deceleration rates in absolute values, respectively. The rate of change of the train acceleration, jerk, can also affect passengers' comfort [49,[81][82][83]. The maximum and minimum admissible effort due to the jerk limitation, F maxJerk (t) and F minJerk (t), are calculated by using Equations (21) and (22).
where j max is the maximum admissible jerk. The total traction or braking effort demanded to the train after applying the comfort restrictions and engine maximum traction and brake effort, F TB (t), is computed by evaluating Equation (23) or Equation (24).
If the train is on a neutral zone, it has to use the regenerative brake to maintain the comfort and auxiliary systems; therefore, the total braking effort to be applied is are the initial and final points of the k-th neutral zone; and z is the total number of neutral zones in the considered line. Finally, the values of the tractive and brake forces to be applied in the train are obtained from the total traction or braking effort simply by considering that if F TB (t) is positive the train will apply traction and it will brake if F TB (t) is negative, as described in Equations (25) and (26).

Eco-Driving Optimization Algorithm
A multi-agent based optimization method called MOPSO (multi-objective particle swarm optimization algorithm) [78] is used to generate energy-efficient speed profiles [44]. Energy consumption/running time Pareto fronts are obtained for different values of the engine traction and brake efficiencies and the linear resistance of the catenary.
This nature-inspired algorithm imitates the behavior of a swarm of insects looking for food. These insects are considered particles that move iteratively within a search space. In this case, each particle corresponds to a speed profile and the running time and energy consumption have to be optimized. A fitness value is assigned to each particle. In biology terms, fitness stands for a degree of adaptability of an insect or how close it is to an abundant source of food. To optimize its fitness, each particle position is updated at every iteration by computing its velocity as a function of its past position with best fitness, pbest, and the global best position found by the whole swarm, gbest. The latter variable accounts for the information shared between the particle ensemble.
To initialize the algorithm, a set of particles with random positions and velocities is generated. At every iteration i, the non-dominated solutions, which are those particles for which there are not any calculated solutions with lower running time and energy consumption, are stored in an archive A. The non-dominated solutions in A are sorted in decreasing order by using a crowing distance (CD) operator [44]. This way, the low-density zones of the Pareto front are given more priority. pbest of each particle is updated and gbest is drawn randomly from the archive, giving priority to the solutions with higher crowding distance, which correspond to the solutions at the top of A.
The positions and velocities of the particles at the next iteration are generated by using Equations (27) and (28),v wherexp j ,vp j are the position and velocity of the j-th particle; w is an inertia constant that weights the previous velocities; c 1 and c 2 are two social factor constants that weight the distance to pbest j ; and gbest, r 1 and r 2 are random numbers drawn uniformly between 0 and 1.
After the position of all the particles in the current iteration has been updated, the new non-dominated solutions are included in the archive and all the newly dominated solutions are deleted. pbest is updated for all the solutions in the archive. The solutions in A are sorted in decreasing order of CD, and finally, the gbest is drawn randomly from A, giving a higher probability to those solutions with higher CD. This procedure is repeated until a certain number of iterations I is reached. Figure 1 shows the MOPSO algorithm flowchart.

Case Study and Results
In this section, the holding speed without braking with final coasting eco-driving strategy is compared against the standard holding speed driving, and the difference in energy consumption is quantified. In both cases, there is a single cruise speed and all the energy regenerated by the electrical brake is assumed to be returned to the power grid. The dependence of the energy consumption (measured at pantograph and at substations) on the engine efficiency and catenary linear resistance is analyzed. The nature-inspired optimization algorithm introduced in Section 5 is used for obtaining the Pareto front of speed profiles for the Talgo-Bombardier class 102 train running on two sections of Spanish high-speed lines. Track grades, curves, speed limits, tunnels, neutral zones and the positions of electrical substations are considered in the line models.
The first section analyzed is an 85.4 km long section between Calatayud and Zaragoza in the Madrid-Barcelona line with a 2 × 25 kV power supply system. The second section is a 38.9 km long line section between Puertollano and Ciudad Real in the Madrid-Sevilla line, with a 1 × 25 kV power supply system. The train has two 8 MW engines of 200 kN of maximum traction effort. The maximum traction curve is shown in Figure 2. The power consumed by the auxiliary systems is 325 kW and in any neutral zone the train has to regenerate enough energy to maintain these systems on, as presented in Equation (24). The running resistance curve is presented in Figure 3. The empty mass of the train is 324 t, and its length is 200 m. Furthermore, apart from the train's empty mass, 50 t-worth of passengers have been taken into account.

Case Study and Results
In this section, the holding speed without braking with final coasting eco-driving strategy is compared against the standard holding speed driving, and the difference in energy consumption is quantified. In both cases, there is a single cruise speed and all the energy regenerated by the electrical brake is assumed to be returned to the power grid. The dependence of the energy consumption (measured at pantograph and at substations) on the engine efficiency and catenary linear resistance is analyzed. The nature-inspired optimization algorithm introduced in Section 5 is used for obtaining the Pareto front of speed profiles for the Talgo-Bombardier class 102 train running on two sections of Spanish high-speed lines. Track grades, curves, speed limits, tunnels, neutral zones and the positions of electrical substations are considered in the line models.
The first section analyzed is an 85.4 km long section between Calatayud and Zaragoza in the Madrid-Barcelona line with a 2 × 25 kV power supply system. The second section is a 38.9 km long line section between Puertollano and Ciudad Real in the Madrid-Sevilla line, with a 1 × 25 kV power supply system. The train has two 8 MW engines of 200 kN of maximum traction effort. The maximum traction curve is shown in Figure 2. The power consumed by the auxiliary systems is 325 kW and in any neutral zone the train has to regenerate enough energy to maintain these systems on, as presented in Equation (24). The running resistance curve is presented in Figure 3. The empty mass of the train is 324 t, and its length is 200 m. Furthermore, apart from the train's empty mass, 50 t-worth of passengers have been taken into account.

Case Study and Results
In this section, the holding speed without braking with final coasting eco-driving strategy is compared against the standard holding speed driving, and the difference in energy consumption is quantified. In both cases, there is a single cruise speed and all the energy regenerated by the electrical brake is assumed to be returned to the power grid. The dependence of the energy consumption (measured at pantograph and at substations) on the engine efficiency and catenary linear resistance is analyzed. The nature-inspired optimization algorithm introduced in Section 5 is used for obtaining the Pareto front of speed profiles for the Talgo-Bombardier class 102 train running on two sections of Spanish high-speed lines. Track grades, curves, speed limits, tunnels, neutral zones and the positions of electrical substations are considered in the line models.
The first section analyzed is an 85.4 km long section between Calatayud and Zaragoza in the Madrid-Barcelona line with a 2 × 25 kV power supply system. The second section is a 38.9 km long line section between Puertollano and Ciudad Real in the Madrid-Sevilla line, with a 1 × 25 kV power supply system. The train has two 8 MW engines of 200 kN of maximum traction effort. The maximum traction curve is shown in Figure 2. The power consumed by the auxiliary systems is 325 kW and in any neutral zone the train has to regenerate enough energy to maintain these systems on, as presented in Equation (24). The running resistance curve is presented in Figure 3. The empty mass of the train is 324 t, and its length is 200 m. Furthermore, apart from the train's empty mass, 50 t-worth of passengers have been taken into account.  Regarding the comfort constraints, the maximum and the service decelerations are 0.4 / , the maximum acceleration is 0.67 / and the maximum admissible jerk is, in absolute value, 0.7 / . The power constant is considered to be = 1. The social factors and in the MOPSO algorithm are equal to 2 and the inertia constant is 0.2. Figure 4 shows the height profile of the line section between Calatayud and Zaragoza. As can be seen, the considered section is mostly downhill. There are three neutral zones, which are represented in magenta in the lower part of the figure. In this subsection, the impact of the efficiency ratio of the motors on the energy consumption is analyzed for both driving strategies. The traction and brake engine efficiencies are considered equal and constant within a range between 85% and 97%. The influence of the linear resistance of the catenary on the energy consumption is also studied. In the case of the Calatayud-Zaragoza section, the catenary configuration is 2 × 25 kV, and different values of its linear resistance, between 21 · 10 Ω/ and 37 · 10 Ω/ , are tested. The MOPSO algorithm was used to generate speed profiles for both standard and eco-driving strategies. For each driving strategy the non-dominated speed profiles in terms of energy consumption and running time define the Pareto front. Figure 4 shows Pareto fronts for the net consumption measured at the electrical substations for the extreme values of the considered interval for the efficiency ratio of the motors and for the linear resistance of the catenary. The Pareto fronts corresponding to the eco-driving strategy are represented in blue, while those corresponding to the standard driving strategy are represented in red.

Calatayud-Zaragoza Case Study
As shown in Figure 5, the consumed energy decreases as the engine efficiency increases. The energy consumption increases for higher values of the linear resistance, since more regenerated energy is dissipated at the catenary. For the fastest speed profilesm the difference in energy consumption is significantly lower than for the rest. These speed profiles are close to the flat-out speed profile; thus, the final coasting phase is very short, and therefore has little relevance in terms of energy-saving, so both driving strategies produce similar results in terms of energy consumption. It can also be observed that the eco-driving strategy produces faster and less energy-consuming speed profiles than the standard driving strategy, as all the solutions for the latter one are dominated by Regarding the comfort constraints, the maximum and the service decelerations are 0.4 m/s 2 , the maximum acceleration is 0.67 m/s 2 and the maximum admissible jerk is, in absolute value, 0.7 m/s 3 . The power constant is considered to be cosϕ = 1. The social factors c 1 and c 2 in the MOPSO algorithm are equal to 2 and the inertia constant w is 0.2. Figure 4 shows the height profile of the line section between Calatayud and Zaragoza. As can be seen, the considered section is mostly downhill. There are three neutral zones, which are represented in magenta in the lower part of the figure.

Calatayud-Zaragoza Case Study
Energies 2020, 13, x FOR PEER REVIEW 10 of 25 Regarding the comfort constraints, the maximum and the service decelerations are 0.4 / , the maximum acceleration is 0.67 / and the maximum admissible jerk is, in absolute value, 0.7 / . The power constant is considered to be = 1. The social factors and in the MOPSO algorithm are equal to 2 and the inertia constant is 0.2. Figure 4 shows the height profile of the line section between Calatayud and Zaragoza. As can be seen, the considered section is mostly downhill. There are three neutral zones, which are represented in magenta in the lower part of the figure. In this subsection, the impact of the efficiency ratio of the motors on the energy consumption is analyzed for both driving strategies. The traction and brake engine efficiencies are considered equal and constant within a range between 85% and 97%. The influence of the linear resistance of the catenary on the energy consumption is also studied. In the case of the Calatayud-Zaragoza section, the catenary configuration is 2 × 25 kV, and different values of its linear resistance, between 21 · 10 Ω/ and 37 · 10 Ω/ , are tested. The MOPSO algorithm was used to generate speed profiles for both standard and eco-driving strategies. For each driving strategy the non-dominated speed profiles in terms of energy consumption and running time define the Pareto front. Figure 4 shows Pareto fronts for the net consumption measured at the electrical substations for the extreme values of the considered interval for the efficiency ratio of the motors and for the linear resistance of the catenary. The Pareto fronts corresponding to the eco-driving strategy are represented in blue, while those corresponding to the standard driving strategy are represented in red.

Calatayud-Zaragoza Case Study
As shown in Figure 5, the consumed energy decreases as the engine efficiency increases. The energy consumption increases for higher values of the linear resistance, since more regenerated energy is dissipated at the catenary. For the fastest speed profilesm the difference in energy consumption is significantly lower than for the rest. These speed profiles are close to the flat-out speed profile; thus, the final coasting phase is very short, and therefore has little relevance in terms of energy-saving, so both driving strategies produce similar results in terms of energy consumption. It can also be observed that the eco-driving strategy produces faster and less energy-consuming speed profiles than the standard driving strategy, as all the solutions for the latter one are dominated by In this subsection, the impact of the efficiency ratio of the motors on the energy consumption is analyzed for both driving strategies. The traction and brake engine efficiencies are considered equal and constant within a range between 85% and 97%. The influence of the linear resistance of the catenary on the energy consumption is also studied. In the case of the Calatayud-Zaragoza section, the catenary configuration is 2 × 25 kV, and different values of its linear resistance, between 21·10 −6 Ω/m and 37·10 −6 Ω/m, are tested. The MOPSO algorithm was used to generate speed profiles for both standard and eco-driving strategies. For each driving strategy the non-dominated speed profiles in terms of energy consumption and running time define the Pareto front. Figure 4 shows Pareto fronts for the net consumption measured at the electrical substations for the extreme values of the considered interval for the efficiency ratio of the motors and for the linear resistance of the catenary. The Pareto fronts corresponding to the eco-driving strategy are represented in blue, while those corresponding to the standard driving strategy are represented in red.
As shown in Figure 5, the consumed energy decreases as the engine efficiency increases. The energy consumption increases for higher values of the linear resistance, since more regenerated energy is dissipated at the catenary. For the fastest speed profilesm the difference in energy consumption is significantly lower than for the rest. These speed profiles are close to the flat-out speed profile; thus, the final coasting phase is very short, and therefore has little relevance in terms of energy-saving, so both driving strategies produce similar results in terms of energy consumption. It can also be observed that the eco-driving strategy produces faster and less energy-consuming speed profiles than the standard driving strategy, as all the solutions for the latter one are dominated by those produced by the former.   In Table 1, two trends can be observed. The difference in energy consumption decreases as the linear resistance of the catenary decreases. It decreases even more significantly when the efficiencies of the engine brake and traction are higher. The energy saving at the electrical substations between the eco-driving strategy and the standard driving is 2.47%-4.69%.
Obviously, the measured energy consumption at the pantograph does not depend on the linear resistance of the catenary. Figure 6 shows the dependence of the energy savings (measured at the pantograph) on the energy efficiency.  In Table 1, two trends can be observed. The difference in energy consumption decreases as the linear resistance of the catenary decreases. It decreases even more significantly when the efficiencies of the engine brake and traction are higher. The energy saving at the electrical substations between the eco-driving strategy and the standard driving is 2.47-4.69%.
Obviously, the measured energy consumption at the pantograph does not depend on the linear resistance of the catenary. Figure 6 shows the dependence of the energy savings (measured at the pantograph) on the energy efficiency. It may be observed that for the considered engine efficiency interval, there is an approximately linear relation between the difference in energy consumption and the engine efficiency. The ecodriving strategy produces speed profiles with lower energy consumption than the standard driving strategy. Results shown in Figure 6 are especially relevant for train operators. The energy saving measured at pantograph of eco-driving with respect to standard driving is between 2.1% and 4.3%.
In Figure 7 the difference in energy consumption at the substations is analyzed for each linear resistance of catenary. Results are shown for different values of the engine efficiency and running time. A grey vertical dashed line is used to highlight the commercial running time, which is 24 min. It may be observed that for the considered engine efficiency interval, there is an approximately linear relation between the difference in energy consumption and the engine efficiency. The eco-driving strategy produces speed profiles with lower energy consumption than the standard driving strategy. Results shown in Figure 6 are especially relevant for train operators. The energy saving measured at pantograph of eco-driving with respect to standard driving is between 2.1% and 4.3%.
In Figure 7 the difference in energy consumption at the substations is analyzed for each linear resistance of catenary. Results are shown for different values of the engine efficiency and running time. A grey vertical dashed line is used to highlight the commercial running time, which is 24 min. Three different scenarios can be distinguished. Close to the flat-out running time, there is a large increase in the difference in energy consumption between the eco-driving and the standard driving strategies; this is due to the action of the final coasting phase in the former driving strategy. Next, there is a scenario in which the difference is approximately constant. For the slowest speed profiles, the difference in energy consumption is significantly lower, as can be seen in Figure 5, and the difference in energy consumption between the considered driving strategies becomes smaller as the running time increases, as the coasting phase will be shorter and therefore its impact will be less Three different scenarios can be distinguished. Close to the flat-out running time, there is a large increase in the difference in energy consumption between the eco-driving and the standard driving strategies; this is due to the action of the final coasting phase in the former driving strategy. Next, there is a scenario in which the difference is approximately constant. For the slowest speed profiles, the difference in energy consumption is significantly lower, as can be seen in Figure 5, and the difference in energy consumption between the considered driving strategies becomes smaller as the running time increases, as the coasting phase will be shorter and therefore its impact will be less significant. If the dependence of the results on the energy braking and traction efficiency is considered, it can be observed that for a fixed running time the difference in energy consumption seems to depend linearly on the difference in energy efficiency. For the same value of the engine efficiency, the difference in energy consumption increases slightly with the linear resistivity of the catenary.
Last, the speed profiles obtained for the commercial time with the two considered driving strategies are presented. Table 2 shows the command vectors for these two speed profiles. Notice that the considered line section ends at kilometric point (k.p.) 306.7 at Zaragoza station (stopping point). The coasting phase in eco-driving starts 18 km before the stopping point. In Figures 8 and 9, the speed profiles are shown. The train speed is represented in blue, the speed limits are represented by red lines and the track altitude is shown in green. The holding speed is represented by an orange dashed line.
Final position of the cruise phase (km) 288.7 Braking curve In Figures 8 and 9, the speed profiles are shown. The train speed is represented in blue, the speed limits are represented by red lines and the track altitude is shown in green. The holding speed is represented by an orange dashed line.   Figure 8 shows that in the holding speed regimes there are few points in which the train goes faster than the cruise speed in the eco-driving strategy, as there are few steep downhill grades. The most significant difference between the two speed profiles in Figures 8 and 9 is after position 288.7 km, in which the train is coasting in Figure 8, while in Figure 9 it is holding its speed until it reaches the braking curve. Near position 300 km the train has to brake due to the presence of a speed limit. In Figures 8 and 9, the speed profiles are shown. The train speed is represented in blue, the speed limits are represented by red lines and the track altitude is shown in green. The holding speed is represented by an orange dashed line.   Figure 8 shows that in the holding speed regimes there are few points in which the train goes faster than the cruise speed in the eco-driving strategy, as there are few steep downhill grades. The most significant difference between the two speed profiles in Figures 8 and 9 is after position 288.7 km, in which the train is coasting in Figure 8, while in Figure 9 it is holding its speed until it reaches the braking curve. Near position 300 km the train has to brake due to the presence of a speed limit.  Figure 8 shows that in the holding speed regimes there are few points in which the train goes faster than the cruise speed in the eco-driving strategy, as there are few steep downhill grades. The most significant difference between the two speed profiles in Figures 8 and 9 is after position 288.7 km, in which the train is coasting in Figure 8, while in Figure 9 it is holding its speed until it reaches the braking curve. Near position 300 km the train has to brake due to the presence of a speed limit.

Puertollano-Ciudad Real Case Study
Finally, results are shown for the second case study, in which eco-driving speed profiles were designed for the considered Talgo-Bombardier class 102 train in a 38.9 km long line section between Puertollano and Ciudad Real. This line section has a catenary configuration of 1 × 25 kV, and three different values of its linear resistance between 54·10 −6 Ω/m and 110·10 −6 Ω/m are considered in this study. In Figure 10 the height profile of the considered line section is represented in blue, while the two neutral zones are represented in magenta in the lower part of the figure. In general terms, this line section is flatter than the Calatayud-Zaragoza line section.

Puertollano-Ciudad Real Case Study
Finally, results are shown for the second case study, in which eco-driving speed profiles were designed for the considered Talgo-Bombardier class 102 train in a 38.9 km long line section between Puertollano and Ciudad Real. This line section has a catenary configuration of 1 × 25 kV, and three different values of its linear resistance between 54 · 10 Ω/ and 110 · 10 Ω/ are considered in this study. In Figure 10 the height profile of the considered line section is represented in blue, while the two neutral zones are represented in magenta in the lower part of the figure. In general terms, this line section is flatter than the Calatayud-Zaragoza line section. Again, speed profiles were generated for the two considered driving strategies by using the MOPSO algorithm. The obtained Pareto fronts are represented in Figure 11, in which they are sorted in terms of the net consumption at the substations. Results are shown for the extreme values of the considered engine efficiency interval, 85% and 97% and for catenary linear resistances 54 · 10 Ω/ , 70 · 10 Ω/ and 110 · 10 Ω/ . The standard and eco-driving speed profiles are represented in red and blue, respectively. The commercial time, which is 13 min, is highlighted by a vertical grey dashed line. Again, speed profiles were generated for the two considered driving strategies by using the MOPSO algorithm. The obtained Pareto fronts are represented in Figure 11, in which they are sorted in terms of the net consumption at the substations. Results are shown for the extreme values of the Energies 2020, 13, 1836 15 of 24 considered engine efficiency interval, 85% and 97% and for catenary linear resistances 54·10 −6 Ω/m, 70·10 −6 Ω/m and 110·10 −6 Ω/m. The standard and eco-driving speed profiles are represented in red and blue, respectively. The commercial time, which is 13 min, is highlighted by a vertical grey dashed line.
Energies 2020, 13, x FOR PEER REVIEW 16 of 25 Figure 11. Pareto fronts for the net consumption at the electrical substations for the extreme values of the engine traction and braking efficiency (85% and 97%) and the three considered values of the linear resistance of the catenary (54 · 10 Ω/ , 70 · 10 Ω/ and 110 · 10 Ω/ ). Figure 11 are similar to those in Figure 5, although this time the difference in energy consumption between two solutions with the same running time is usually larger than in the Calatayud-Zaragoza case. For a fixed running time, the energy consumption is larger when the engine is more inefficient and when the catenary has a larger linear resistance. The speed profiles in the holding speed without a braking Pareto curve dominate those in the holding speed with a braking Pareto curve. Table 3 shows the difference in energy consumption at the substations between the two considered driving strategies for the commercial running time in the Puertollano-Ciudad Real section. Results are shown for different values of the engine efficiency and the linear resistance of the catenary. Table 3. Difference (in %) in energy consumption at the electrical substations between the holding speed with braking and the holding speed without braking and final coasting driving strategies for different values of the energy traction and braking efficiency and the linear resistance in the Puertollano-Ciudad Real section for the commercial running time (13 min).   Figure 11 are similar to those in Figure 5, although this time the difference in energy consumption between two solutions with the same running time is usually larger than in the Calatayud-Zaragoza case. For a fixed running time, the energy consumption is larger when the engine is more inefficient and when the catenary has a larger linear resistance. The speed profiles in the holding speed without a braking Pareto curve dominate those in the holding speed with a braking Pareto curve. Table 3 shows the difference in energy consumption at the substations between the two considered driving strategies for the commercial running time in the Puertollano-Ciudad Real section. Results are shown for different values of the engine efficiency and the linear resistance of the catenary. Table 3. Difference (in %) in energy consumption at the electrical substations between the holding speed with braking and the holding speed without braking and final coasting driving strategies for different values of the energy traction and braking efficiency and the linear resistance in the Puertollano-Ciudad Real section for the commercial running time (13 min Once again, the lower the linear resistance is, the lower the difference in energy consumption is. Moreover, the difference in consumption between two speed profiles produced by using the two considered strategies is higher when the train engine is less efficient. The difference in energy consumption between the considered driving strategies for the commercial running time is between 9.70% and 12.55%. Figure 12 shows the energy consumption measured at the pantograph for the commercial time speed profiles. As can be observed, for an engine efficiency of 85%, the eco-driving strategy can be used for achieving a 10.8% energy saving with respect to the standard driving strategy, while for an engine efficiency of 97% the energy saving between the two strategies is 8.3%.

Results presented in
Energies 2020, 13, x FOR PEER REVIEW 17 of 25 Once again, the lower the linear resistance is, the lower the difference in energy consumption is. Moreover, the difference in consumption between two speed profiles produced by using the two considered strategies is higher when the train engine is less efficient. The difference in energy consumption between the considered driving strategies for the commercial running time is between 9.70% and 12.55%. Figure 12 shows the energy consumption measured at the pantograph for the commercial time speed profiles. As can be observed, for an engine efficiency of 85%, the eco-driving strategy can be used for achieving a 10.8% energy saving with respect to the standard driving strategy, while for an engine efficiency of 97% the energy saving between the two strategies is 8.3%. The difference in energy consumption between the considered driving strategies is analyzed in Figure 13 for the considered linear resistances of the catenary and engine efficiencies and for different running times. The difference in energy consumption between the considered driving strategies is analyzed in Figure 13 for the considered linear resistances of the catenary and engine efficiencies and for different running times.
The difference in energy consumption grows up to a maximum close to the commercial running time. The commercial running time (13 min, that is 780 s) is highlighted by means of a grey dashed line. The more-or-less constant scenario in Figure 7 is at this point concave; this indicates that the shape of the difference in energy consumption depends on the height profile of the track. This time, the difference in energy consumption is larger when considering different values of the linear resistivity, as there is a larger difference between those values. Again, the difference in consumption seems to decrease linearly with the engine efficiency. Energies 2020, 13, x FOR PEER REVIEW 18 of 25 The difference in energy consumption grows up to a maximum close to the commercial running time. The commercial running time (13 min, that is 780s) is highlighted by means of a grey dashed line. The more-or-less constant scenario in Figure 7 is at this point concave; this indicates that the In Table 4 the command vectors of the two driving speed profiles that satisfy the commercial running time are presented. The eco-driving final coasting phase is approximately 10 km long and its cruise speed is 17.4 km/h higher. Figures 14 and 15 show the commercial speed profiles obtained for the eco-driving and standard driving strategies, respectively. The train speed is represented in blue, the speed limits are represented by red lines and the track altitude is shown in green. The holding speed is represented by an orange dashed line. resistivity, as there is a larger difference between those values. Again, the difference in consumption seems to decrease linearly with the engine efficiency. In Table 4 the command vectors of the two driving speed profiles that satisfy the commercial running time are presented. The eco-driving final coasting phase is approximately 10 km long and its cruise speed is 17.4 km/h higher. Figures 14 and 15 show the commercial speed profiles obtained for the eco-driving and standard driving strategies, respectively. The train speed is represented in blue, the speed limits are represented by red lines and the track altitude is shown in green. The holding speed is represented by an orange dashed line.   In this case, the train coasts approximately in 25% of the total length of the line section. In comparison to Figures 8 and 9, in this section the role of the final coasting phase is more important. This explains the fact that the difference in energy consumption between the considered driving strategies in percentage is much higher for the Puertollano-Ciudad Real section. In Table 4 the command vectors of the two driving speed profiles that satisfy the commercial running time are presented. The eco-driving final coasting phase is approximately 10 km long and its cruise speed is 17.4 km/h higher. Figures 14 and 15 show the commercial speed profiles obtained for the eco-driving and standard driving strategies, respectively. The train speed is represented in blue, the speed limits are represented by red lines and the track altitude is shown in green. The holding speed is represented by an orange dashed line.   In this case, the train coasts approximately in 25% of the total length of the line section. In comparison to Figures 8 and 9, in this section the role of the final coasting phase is more important. This explains the fact that the difference in energy consumption between the considered driving strategies in percentage is much higher for the Puertollano-Ciudad Real section. In this case, the train coasts approximately in 25% of the total length of the line section. In comparison to Figures 8 and 9, in this section the role of the final coasting phase is more important. This explains the fact that the difference in energy consumption between the considered driving strategies in percentage is much higher for the Puertollano-Ciudad Real section.

Discussion
The studies developed in this paper have been carried out with the objective of assessing the effectiveness of eco-driving under different electrical scenarios and to answer the question about whether or not eco-driving is useful in high-receptivity networks. The results obtained suggest that the answer to the question is that it is useful to apply eco-driving. Besides, it has to be taken into account that the value of energy saving because of efficient driving will depend on the efficiency of the electrical chain (catenary loses and motor efficiency) and on the running time.
The cases studied in Section 6 are two cases of an alternate current power supply system with high receptivity. These systems are, by nature, much more effective using regenerate energy than direct current systems. In alternate current systems, the energy generated by train braking can be used by other trains or sent back to the distribution network for other uses. However, there is part of the regenerated energy that will be lost anyway because of the catenary losses and the motor's efficiency. Saving part of the lost energy is the main benefit of eco-driving, and the results support this. It has been shown that, the lower the energy lost in the motor and catenary, the lower the difference between efficient non-efficient driving.
The results obtained in Puertollano-Ciudad Real show energy savings produced by eco-driving according to the typical values between 10% and 15% in the literature [68]. On the other hand, the results obtained in Calatayud-Zaragoza present values that are further reduced. One reason to explain this is the track profile. Although the average grade in both cases is downhill, the difference in altitude between the end and start of the journey is higher in Calatayud-Zaragoza. The presence of steeper downhills leads to lower energy consumption results and lower differences between different speed profiles.
Catenary losses can also explain why the energy reduction figures obtained by the eco-driving strategy in the Calatayud-Zaragoza stretch are lower than those obtained in the Puertollano-Ciudad Real stretch (the 2 × 25 kV system is more efficient than the 1 × 25 kV system because of lower rate of power losses).
The influence of the commercial running time cannot be neglected. In both cases, it can be observed that there are three scenarios of energy saving due to eco-driving depending on the running time. Close to the flat-out running time, the energy saving that can be obtained by efficient driving is zero because there is no possibility of modifying the driving to obtain the minimum running time. From there, there is a small interval of running time wherein the energy saving provided by eco-driving increases as the running time increases. After that, there is a zone of running time values in which the energy saving rate is approximately constant. This zone typically contains the commercial running time for the journey. Finally, for the largest running times, the difference in energy consumption between efficient and non-efficient driving is gradually reduced.
Therefore, a second order reason for the reduced energy saving figures in Calatayud-Zaragoza is that the commercial running time of this journey is higher than usual and it is located in the running time interval in which the energy savings provided by eco-driving are reduced.
Despite this, the minimum energy saving provided by eco-driving is relevant (2.5%) in a scenario of 97% motor efficiency. This motor efficiency is very high compared with its typical values, which are between 85% and 90%. Therefore, taking into account the typical motor efficiency, the minimum energy saving provided by eco-driving can reach more that 3.5%.
With these results, it can be said that high network receptivity to braking energy can mitigate the effect of eco-driving, but it is still significant.
Obviously, results can vary in a field-test because of simulation error. However, as explained in Section 2, the value of simulation error in energy consumption is 0.4%. This error can slightly quantitatively modify the results obtained but cannot qualitatively modify them and the trends obtained.

Conclusions
The basis of eco-driving strategy is the substitution of the use of braking with coasting periods. Trains equipped with regenerative braking allow for recovering the kinetic energy during deceleration phases, which can be used by other trains or sent back to the utility grid in AC power supply systems or in DC systems equipped with reversible substations. In this context of high use of regenerated energy, it is necessary to analyze whether eco-driving design based on coast commands has an impact on energy savings. In this work, an eco-driving strategy has been compared to a standard driving strategy. The former one consists on holding a cruise speed without braking (if braking is needed to maintain that speed, it is not applied and the speed increases) and coasting before braking. The latter one consists of holding the cruise speed and braking if it is necessary to do so. A comparison has been performed in two real case studies by simulating a train on two high-speed Spanish lines. One of them is a section of the Madrid-Barcelona high-speed line, which is fed by a 2 × 25 kV power supply system. The other one is a section of the Madrid-Sevilla high-speed line, which is fed by a 1 × 25 kV Energies 2020, 13, 1836 20 of 24 power supply system. The influence of the regenerated energy losses at the catenary-which depend on its linear resistance, and the engine traction and braking efficiency, which are considered constant and equal-on the energy consumption of the speed profiles generated for both driving strategies has been studied. To generate the driving commands, a nature-inspired algorithm called MOPSO has been used for both driving strategies. Speed profiles have been obtained by simulating driving commands. For each driving strategy and for some values of the linear resistance and engine efficiency, the Pareto curves consisting of the non-dominated speed profiles in terms of energy consumption and running time have been computed. The difference in energy consumption between the considered driving strategies has been calculated for different values of the catenary resistance, engine efficiency and running time, with a special focus on the commercial time. For the commercial running time, the eco-driving strategy produces energy savings with respect to the standard driving strategy between 2.48% and 4.69% at the pantograph for the first case study and between 9.7% and 12.55% for the second.
For a specific running time, the holding speed without braking and final coasting eco-driving strategy produces solutions that consume less energy than the holding speed with braking standard driving strategy. The difference in energy consumption between the considered driving strategies is larger when the losses at the catenary are more significant and when the engine is less efficient.
It was found that the holding speed without braking and final coasting driving strategy produces speed profiles more energy-efficient than the standard driving strategy, mostly due to the final coast phase. For slower speed profiles, the energy consumption is much lower, as the cruise speed is lower; the final coast phase takes place more close to the brake curve at the end of the section; and then the difference in energy consumption between the considered driving strategies is smaller.
Finally, it can be concluded that, although the power supply network presents a high receptivity to braking energy, the eco-driving design can produce significant energy savings.