# Modeling, Simulation and Analysis of On-Board Hybrid Energy Storage Systems for Railway Applications

Pablo Arboleya, Islam El-Sayed, Bassam Mohamed, Clement Mayet

<figure><img src="image_1.png" /></figure>

## To cite this version:

Pablo Arboleya, Islam El-Sayed, Bassam Mohamed, Clement Mayet. Modeling, Simulation and Analysis of On-Board Hybrid Energy Storage Systems for Railway Applications. Energies, 2019, 12 (11), pp.2199. 10.3390/en12112199. hal-02155600

HAL Id: hal-02155600

https://hal.science/hal-02155600

Submitted on 4 Sep 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Modeling, Simulation and Analysis of On-Board Hybrid Energy Storage Systems for Railway Applications †

Pablo Arboleya <fnref n="1" />,* ID, Islam El-Sayed <fnref n="1" />, Bassam Mohamed <fnref n="1" /> ID and Clement Mayet <fnref n="2" /> ID

Received: 30 April 2019; Accepted: 6 June 2019; Published: 10 June 2019

<figure><img src="image_3.png" /></figure>

check for
updates

**Abstract:** In this paper, a decoupled model of a train including an on-board hybrid accumulation system is presented to be used in DC traction networks. The train and the accumulation system behavior are modeled separately, and the results are then combined in order to study the effect of the whole system on the traction electrical network. The model is designed specifically to be used with power flow solvers for planning purposes. The validation has been carried out comparing the results with other methods previously developed and also with experimental measurements. A detailed description of the power flow solver is beyond the scope of this work, but it must be remarked that the model must be used with a solver able to cope with the non-linear and non-smooth characteristics of the model. In this specific case, a modified current injection-based power flow solver has been used. The solver is able to incorporate also non-reversible substations, which are the most common devices used currently for feeding DC systems. The effect of the on-board accumulation systems on the network efficiency will be analyzed using different real scenarios.

**Keywords:** rail transportation power systems; DC power systems; load flow analysis; power systems modeling; load modeling

## 1. Introduction

The use of power flow algorithms for planning traction networks is a widely-accepted technique [1,2]. However, the use of accurate models of the network and the trains may result in very complicated simulations [3]. In [4], the authors proposed a methodology to perform a fast estimation of the aggregated railway power system and traffic performance. In [5], for instance, the authors tried to reduce the computational burden of the whole simulation, proposing a compression technique that reduces the number of necessary simulations for evaluating the performance of the traction network. In [6], the energetic macroscopic representation (EMR) was used in order to simplify the development of the train and network mathematical models. This approach has been use extensively. For instance, in [7], the authors compared different traction substation models. The models presented using the EMR were extremely accurate and precise [8], and they may describe every single part of the system and are appropriate for applications in which the transient behavior of the trains is critical, such as hardware in the loop applications [9]. A very detailed model of an electric train fed by supercapacitors connected to fuel cells was presented in [10]. Pseudo-static models have been traditionally less accurate, but much faster [11]. They are usually combined with power flow solvers that must be able to manage the

<footnote n="1">¹ LEMUR Research Group, Department of Electrical Engineering, University of Oviedo, Campus of Gijón, 33204 Gijón, Spain; islam.aim@gmail.com (I.E.-S.); engbassam@gmail.com (B.M.)</footnote>

<footnote n="2">² SATIE—UMR CNRS 8020, Conservatoire National des Arts et Métiers, HESAM University, F-75003 Paris, France; clement.mayet@lecnam.net</footnote>

<footnote>* Correspondence: arboleyapablo@uniovi.es</footnote>

<footnote>† This paper is an extended version of our paper published in 2017 IEEE Vehicle Power and Propulsion Conference (VPPC), Belfort, France, 11–14 December 2017, pp. 1–5.</footnote>
complex non-smooth and non-derivable characteristics presented by the new power electronic devices installed within the railways [12].

In this paper, a very simple model for representing a train equipped with a hybrid energy storage system is presented. The combination of regenerative braking with the energy storage system is presented as one of the most effective ways to increase the overall efficiency of traction systems [13–15]. Schedule optimization together with these sets of technologies can be the key to very important energy savings [16–18]. In [19], the savings derived from the installation of an energy storage device in a diesel-electric locomotive were also evaluated. In [20], a method for optimization of the size of the storage devices in order to maximize the energy recovered by the regenerative braking system was presented. An accurate sizing of all the elements in the system was critical for a proper energy management within the system [21–23]. In other cases, the storage systems were installed off-board at the substation level [24]. Many applications, like the one presented in this paper, use hybrid energy storage systems based on batteries and ultracapacitors; usually, the battery covers a base load, and the ultracapacitors are used for feeding the peak power. In order to connect two different storage technologies to the same train, specific power electronic topologies are being investigated, the most common being the DC/DC converter that connects the DC traction network with the storage devices, which uses IGBT technology with a topology that allows the interleaved multi-channel buck-boost operation principle. This converter can adapt the voltage in the energy storage to the catenary voltage. This configuration provides the flexibility required by the operating conditions. Sometimes, series connection of the ultracapacitors or battery packs is needed in order to achieve a specific voltage level. Usually, there are parallel connected IGBT branches that can be fired at the same time or we can shift the firing to reduce the current ripple. In the case of hybrid energy storage systems, there are IGBT branches connected to the battery packs and other IGBT branches connected to the ultracapacitors. Other hybrid topologies use also a common connection; in these cases, the battery packs share the common bus of the ultracapacitors. Usually, buck mode is used for charging the ultracapacitors and the battery. In this situation, the converter drivers control the upper IGBTs, keeping open the devices placed at the bottom of the branch. The current flows from the train DC bus to the energy storage modules. The boost mode (for discharging) uses the IGBTs placed in the bottom, keeping open the ones at the top of the branch. This will produce a current flow from the energy storage devices to the train DC bus. The results obtained using the pseudostatic simulations proposed in this paper are a valuable input for the designers of the energy storage devices. However, a detailed description of the power electronic topology or the storage technology is far beyond the scope of this paper. In order to determine the size of the battery and the ultracapacitors, it is necessary to assess the performance of the system using an ideal integrated energy storage system like the one presented in this paper.

Wayside energy storage can be an alternative to the on-board accumulation systems and in some cases can also be better. In [12,25], the authors mentioned that for choosing the correct infrastructure/train configuration, multiple cases should be analyzed, and all possible technologies should be compared, like for instance reversible substations, wayside energy storage, and on-board energy storage. In order to take advantage of the power regenerated by the trains, the most common solution used in the past was: (1) to optimize the schedule, so when one train is accelerating, another train is braking (this solution is quite difficult to achieve without automatic driving systems); (2) the use of reversible substations; in this case, railway manufacturers are quite reluctant to change the conventional non-controlled diode-based substations because they are inexpensive and robust, so nearly all reversible substations are formed with an IGBT bridge in parallel with the diode bridge; the infrastructure keeps using the diodes to transfer power from AC to DC, and they use the IGBT to put power back into the DC system. Currently, with the drastic reduction in the prices of energy storage systems, it is not clear which is going to be the solution that will prevail, reversibility or accumulation. Probably, accumulation will be a much more reasonable solution in the near future since it helps to maintain the voltages in the DC traction network within the limits, but also increases the overall efficiency of the system. In [25], different configurations were compared, and we can
observe that on-board and wayside energy storage are better solutions in terms of efficiency. On-board accumulation is slightly better in terms of efficiency because the losses in the DC system are slightly lower, but wayside energy storage can be more cost-effective since a single accumulation device can be shared by all the trains in the system. Similar conclusions were obtained in [26]; in this work, the use of wayside energy storage improved the overall efficiency by more than 20% and also improved the controllability and the voltage profile of the network. Similar conclusions were obtained in [24,27] supporting the previous conclusions. In [24], the use of energy storage at the substation level reduced the annual cost of energy by 13%. We would like to point out that the main contribution of this paper is the mathematical model; a detailed comparison between the infrastructure alternatives considering wayside and on-board energy storage and reversible substations is beyond the scope of the work, which is focused specifically on providing an accurate and efficient model to simulate on-board energy storage.

What this work presents is a decoupled model that considers the train and the storage system as separate devices located at the same electrical point and interacting with each other and with the network. The authors did not consider the input filter, and the losses in this device can be estimated and added to the electromechanical and electrochemical device. This approach reduces the complexity of the mathematical model without reducing the accuracy when compared with previous proposed models [28]. The model has been embedded in a commercial software package that uses a modified current injection method to solve the power flow problem [12]; the solver allows combining regenerative trains with diode-based non-reversible substations. It must be remarked that even when there is the possibility of combining regenerative trains with Voltage Source Converter (VSC)-based reversible substations [29,30], in the vast majority of the cases, the existing infrastructures and also the new ones use non-controlled and non-reversible substations due to their robustness, reliability, and low cost [31]. The combination of the two technologies must be done carefully. It should be noticed that without a very precise coordination, the power regenerated by one train could be used to feed another train, but if there is no power demand when the train is regenerating, the voltage in the network will increase, and the substations will be blocked. In such cases, usually the voltage limit is reached, and most of the regenerated power must be burned in the rheostatic braking system. The use of on-board accumulation systems could alleviate this situation [32]. However, it will be demonstrated, the size and effectiveness of the storage system depend highly on the trains' power profile and the schedule.

The results obtained using this kind of simulations are very useful also to study how the accumulation systems are cycled in real applications. There are many on-going research projects that try to predict the life of the accumulation systems depending on the cycling. For instance, the common thought is that load leveling the charge and discharge of the battery packs can contribute to extending the life of the batteries. Recent studies revealed unexpected results. In [33], the authors cycled LiNiCoAl and LiFePO₄ modules with 18,650 cells 750 times in a period of six months, and they concluded that the degradation was higher with constant current loads than using dynamic pulse profiles. There is still much research to do in this direction, and every new storage technology requires these kinds of studies. The output of the pseudo-static simulation could be a valuable input for the researchers that try to reproduce real dynamic cycles in accumulation systems.

In the next section, the mathematical model of the train plus the energy storage device in different working modes is explained. Section 3 presents the results obtained using the proposed model in a real traction network. This section has four subsections: in the first three subsections, we describe the feeding infrastructure, the rolling stock, and the proposed scenarios. In the last subsection a deep analysis of the results is presented. Conclusions are stated in Section 4.

## 2. Mathematical Model of the Train and the Storage System

The proposed train model including the storage system has been simplified in order to lighten the computational burden. In [28], the storage system and the traction equipment of the train were considered as a single element; in this case, the storage system and the traction system separation
will substantially simplify the solving procedure. This simplification in the model will not reduce the accuracy because both models are equivalent. The time savings in simulation were around 20%. The solving procedure in this case allows the decoupling between the storage system and the traction equipment, and this is the main reason for this simplification in the mathematical model. From now on, we will use the term train when referring to the traction and breaking device, and storage will be reserved for the accumulation device. Both devices will be modeled as separate, but coupled models.

The storage control will consider the train behavior, as well as the catenary voltage, so we will assume that the train and the storage will be traveling along the grid separated, but connected to the same catenary point.

The joint power and current can be easily calculated as the sum of the current of the train and the storage. For the sign criteria, we will consider that the train will have positive current and power in traction mode when it is absorbing power from the electric grid. In the other case, the storage will have positive power and current when it is in discharging mode. In Figure 1a,b, the schematic representation of the train and storage mathematical models is depicted in traction and braking modes, respectively. We only consider two situations that are the most common ones, the train in traction mode with the storage system discharging and the train in braking mode with the storage system in charging mode.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Equivalent mathematical model of the train plus storage system (a) Traction mode. (b) Braking mode.</figcaption></figure>

## 2.1. Train Plus Storage Working in Traction Mode

The behavior of the train and the storage in traction mode is summarized in Figure 1a, and the functions $f_{t1}$, $g_{d1}$, and $g_{d2}$ are depicted in the curves presented in Figure 2a-c. In the next subsections, a detailed description of the train and the storage model will be presented.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Subfigures (a–c) represent the functions that model the behavior of the train and the storage in traction mode. Subfigures (d–f) represent the functions that model the behavior of the train and the storage in braking mode. (a) Train overcurrent protection in traction mode. (b) Maximum power that can be extracted from the storage system in traction mode as a function of the energy level. (c) Available discharging power in the storage system as a function of the voltage (overcurrent protection of the storage in traction mode). (d) Train squeeze control in braking mode. (e) Maximum power that can be injected into the storage system in braking mode as a function of the energy level. (f) Available charging power in the storage system as a function of the voltage (overcurrent protection of the storage system in braking mode).</figcaption></figure>

### 2.1.1. Train Behavior (Positive Power in Traction Mode)

The train model will read first the power reference provided by the external software; we will refer to this power as the mechanical reference power ($P_{mech,ref}^{train}$). It is possible to use coupled models that generate the mechanical power reference during the electrical simulation [34]. However, in most of the cases, the train mechanical model and the electrical problem are solved separately. To obtain the train electrical reference power ($P_{ref}^{train}$), the mechanical power will be divided by the efficiency of the electromechanical conversion that will include the motors, converters, and the rest of the equipment. This efficiency will be labeled as train efficiency in traction mode ($v_t^{train}$). This operation can be executed before launching the solver with all positive powers provided by the external software package, so the input file for the solver will already consider this efficiency according to the next expression:

$$P_{ref}^{train} = \frac{P_{mech,ref}^{train}}{v_t^{train}} \qquad (1)$$

It must be pointed out that this electric power reference ($P_{ref}^{train}$) is not the final electric power that the train will demand from the grid ($P_{train}$). The overcurrent protection can limit the power requested from the grid when the catenary voltage ($V_{cat}$) is very low for protection purposes.

Each train will have two configuration parameters to set the overcurrent protection ($V_{reg,t}^{train}$ and $V_{min}^{train}$); in Figure 2a, the function $f_{t1}$ is depicted, and the mathematical expression of such function is expressed in (2).
As can be observed, the actual power consumed by the train depends on the catenary voltage and the train electric power reference, as well as the parameters that define the overcurrent protection ($P^{train} = f_{t1}(P^{train}_{ref}, V^{train}_{reg,t}, V^{train}_{min}, V_{cat})$).

$$p^{train} = \begin{cases} 0 & V_{cat} \le V^{train}_{min} \\ \frac{V_{cat} - V^{train}_{min}}{V^{train}_{reg,t} - V^{train}_{min}} p^{train}_{ref} & V^{train}_{min} < V_{cat} < V^{train}_{reg,t} \\ p^{train}_{ref} & V_{cat} \ge V^{train}_{reg,t} \end{cases} \qquad (2)$$

### 2.1.2. Storage Behavior (Negative Power in Traction Mode)

The maximum amount of power that can be extracted from the primary active power source of the storage system (ultracapacitors, batteries, flywheel, or whatever other technology used) is labeled as $P^{STO}_{max}$ and depends on the level of charge of the storage system ($E^{STO}$) with the function $g_{d1}$ (see Figure 2b) that is expressed in (3).

$$P^{STO}_{max} = \begin{cases} P^{STO}_{rated,d} & E^{STO} \le E^{STO}_{reg,d} \\ \frac{E^{STO}}{E^{STO}_{reg,d}} P^{STO}_{rated,d} & E^{STO} > E^{STO}_{reg,d} \end{cases} \qquad (3)$$

The maximum discharging power will be zero if the energy stored is zero, and it will increase in a linear way until a specific energy level is achieved ($E^{STO}_{reg,d}$). The maximum discharging power will be constant for the interval between $E^{STO}_{reg,d}$ and the maximum energy that the system can store ($E^{STO}_{max}$). This maximum power ($P^{STO}_{max}$) can be calculated at each instant once the storage energy level is updated, but before launching the power flow algorithm, because it is independent of the catenary voltage and the train behavior.

The maximum power available from the storage primary source ($P^{STO}_{max}$) has to be multiplied by the electrochemical conversion process efficiency in discharging mode ($v^{STO}_{d}$) to obtain the maximum available electric power to be injected into the system ($P^{STO}_{max,av}$). This operation can be also executed before launching the power flow solver.

$$P^{STO}_{max,av} = v^{STO}_{d} P^{STO}_{max} \qquad (4)$$

Due to the storage system overcurrent protection, the maximum available power could be constrained if the voltage is too low by means of the function $g_{d2}$ (see Figure 2c). The output of this function is the actual available power that the storage device can inject into the system considering already the catenary voltage constraint ($P^{STO}_{av}$). To define the function $g_{d2}$, we need to define also the regulation parameters ($V^{STO}_{reg}$ and $V^{STO}_{min}$). The analytical expression of the function can be found in (5). Each train can have different regulation parameters. In the case of the storage system, these two parameters will be common for charging and discharging mode.

$$p^{train} = \begin{cases} 0 & V_{cat} \le V^{STO}_{min} \\ \frac{V_{cat} - V^{STO}_{min}}{V^{STO}_{reg} - V^{STO}_{min}} P^{STO}_{max,av} & V^{STO}_{min} < V_{cat} < V^{STO}_{reg} \\ P^{STO}_{max,av} & V_{cat} \ge V^{STO}_{reg} \end{cases} \qquad (5)$$

### 2.1.3. Coupling between the Train, Storage, and Network

Regarding the interaction between the train, the storage device, and the grid, many operational philosophies can be considered. In this particular case, one of the most popular working modes in real applications has been adopted, giving priority to the energy extraction from the storage system. When the train demands a given amount of power, if this power is available in the accumulation system, it is
going to be extracted from it. If not, the train will extract the available power from the storage system, and the rest will be imported from the network. With this control philosophy, the actual power that the storage system is going to inject into the system is going to be calculated following the next expression:

$$P^{STO} = \min(P^{train}, P_{av}^{STO}) \qquad (6)$$

Due to the employed power flow solving procedure [11], we will use the previous iteration of the catenary voltage to calculate first the power of the train ($P^{train}$), then we will calculate the storage power ($P^{STO}$). The catenary net power, representing the whole set (train plus storage device), can be calculated as $P^{cat} = P^{train} - P^{STO}$. Once the catenary power is calculated, a new iteration of the power flow will be launched to obtain a new value of the catenary voltage. The simulation will stop once the difference between all catenary voltages in all nodes in two successive iterations is lower than a specific tolerance.

Once the convergence is achieved, the energy level at the storage system must be updated using the following expression:

$$E_{t+1}^{STO} = E_{t}^{STO} - (1/v_{d}^{ACR}) \cdot \Delta t \cdot P^{ACR} \qquad (7)$$

## 2.2. Train Plus Storage Working in Braking Mode

An analogous procedure can be used for describing the train and the storage device working in braking mode. Such behavior is summarized in Figure 1b and the functions $f_{b1}$, $g_{c1}$, and $g_{c2}$ depicted in the curves represented in Figure 2d-f. In the next subsections, a detailed description of the train and the storage model will be presented.

### 2.2.1. Train Behavior (Positive Power in Traction Mode)

Again, the train model will read first the power reference provided by the external software ($P_{mech,ref}^{train}$). To obtain the train electrical reference power ($P_{ref}^{train}$), the mechanical power will be multiplied by the efficiency of the electromechanical conversion (see Equation (8)). This efficiency will be labeled as the train efficiency in braking mode ($v_{b}^{train}$).

$$P_{ref}^{train} = v_{b}^{train} P_{mech,ref}^{train} \qquad (8)$$

The squeeze control can limit the power requested from the grid when the catenary voltage ($V_{cat}$) is very high for protection purposes. Each train will have two configuration parameters to set the squeeze control ($V_{reg,b}^{train}$ and $V_{max}^{train}$). In Figure 2d, the function $f_{b1}$ is depicted. The mathematical expression of such a function is expressed in (9). As can be observed, the actual power that the train is going to regenerate depends on the catenary voltage and the train electric power reference, as well as the parameters that define the overcurrent protection ($P_{ref}^{train} = f_{b1}(P_{ref}^{train}, V_{reg,b}^{train}, V_{max}^{train}, V_{cat})$).

$$P^{train} = \begin{cases} P_{ref}^{train} & V_{cat} \le V_{reg,b}^{train} \\ \frac{V_{cat} - V_{max}^{train}}{V_{reg,b}^{train} - V_{max}^{train}} P_{ref}^{train} & V_{reg,b}^{train} < V_{cat} < V_{max}^{train} \\ 0 & V_{cat} \ge V_{max}^{train} \end{cases} \qquad (9)$$

### 2.2.2. Storage Behavior (Negative Power in Braking Mode)

The maximum amount of power that can be injected into the primary active power source of the storage system ($P_{max}^{STO}$) depends on the level of charge of the storage system ($E^{STO}$) with the function $g_{c1}$ (see Figure 2e) that is expressed in (10).

$$P_{max}^{STO} = \begin{cases} \frac{E^{STO} - E_{max}^{STO}}{E_{reg,c}^{STO} - E_{max}^{STO}} P_{rated,c}^{STO} & E^{STO} \le E_{reg,d}^{STO} \\ P_{rated,d}^{STO} & E^{STO} > E_{reg,d}^{STO} \end{cases} \qquad (10)$$
Again, this maximum power ($P_{max}^{STO}$) can be calculated at each instant once the storage energy level is updated, but before launching the power flow algorithm, because it is independent of the catenary voltage and the train behavior.

The maximum power that could be injected into the storage referring to the catenary side of the storage system ($P_{max,av}^{STO}$) has to be calculated using the electrochemical conversion process efficiency in charging mode ($v_c^{STO}$) according to (11). This operation can be also executed before launching the power flow solver.

$$P_{max,av}^{STO} = \frac{P_{max}^{STO}}{v_d^{STO}} \qquad (11)$$

The overcurrent protection of the storage system is considered also in braking mode using the function $g_{c2}$ (see Figure 2f). The output of this function is the actual available power that the storage device can extract from the system considering already the catenary voltage constraint ($P_{av}^{STO}$). To define the function $g_{c2}$, we need to define also the regulation parameters ($V_{reg}^{STO}$ and $V_{min}^{STO}$). The analytical expression of the function is similar to the one for the traction mode (see Equation (5)).

### 2.2.3. Coupling between the Train, Storage, and Network

To be consistent with the control philosophy previously described, the system will assign priority to the energy injection into the storage system. When the train regenerates a given amount of power, the storage system will try to use it (if possible) to increase its energy level. If not, the train will inject the maximum allowed power into the storage system, and the rest will be injected into the network. The actual power that the storage system is going to use for charging is going to be calculated following the next expression:

$$P^{STO} = \max(P^{train}, P_{av}^{STO}) \qquad (12)$$

The catenary net power, representing the whole set (train plus storage device) can be calculated as $P_{cat} = P^{train} - P^{STO}$. At each instant, the energy level of the storage system can be calculated as follows:

$$E_{t+1}^{STO} = E_{t}^{STO} - v_{c}^{STO} \cdot \Delta t \cdot P^{STO} \qquad (13)$$

## 3. Results' Analysis

In this section, we will analyze in depth the effect of the on-board accumulation system on the trains and network behavior in different scenarios. This section will have four subsections. In the first one, we will describe the basic feeding infrastructure, lines, and substations. In the second subsection, we will describe the rolling stock used in the simulations. The third subsection will be focused on the description of the different scenarios that will be analyzed in Section 4.

### 3.1. Feeding Infrastructure

The case study in this article will focus on the study of a real network consisting of two lines of 30.84 km and 36.93 km. The voltage level of the system was 3000 V. The simplified diagram of the network can be observed in Figure 3. Line blue is the longest; it had nine stops and four electrical nodes labeled as S1, S2, S3, and S4. The red line shared the first two electrical nodes with the blue line, and it had 17 stops and six electrical nodes labeled as S1, S2, S3, S4, S5, and S6. There were a total of eight electrical nodes and seven lines. Among the eight nodes, only three of them represented feeding substations; the rest were just nodes without any connection with the AC system. The three substations were placed in the nodes S3 and S5 of the red line and in the node S3 of the blue line. The three substations had the same characteristics. All of them were composed by a power transformer with rated power of 3 MW and a short circuit voltage of 5%; the no load output voltage of the rectifier was 3000 V, and the voltage at rated load (1000 A) was 2880 V. The equivalent impedance of each of the three substations in forward mode (AC to DC) was 270 mΩ. The equivalent impedance of the
overhead conductor and the rails (return circuit) were respectively 28.605 mΩ/km and 7 mΩ/km. In Table 1, the lengths of the different segments of the red and blue lines can be observed.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Schematic representation of the case study.</figcaption></figure>

<table><caption>**Table 1.** Length of the different line segments in km.</caption><thead><tr><th></th><th>S1 to S2</th><th>S2 to S3</th><th>S3 to S4</th><th>S4 to S5</th><th>S5 to S6</th></tr></thead><tbody><tr><td>Red Line</td><td>4.316</td><td>0.500</td><td>13.800</td><td>7.848</td><td>4.378</td></tr><tr><td>Blue Line</td><td>4.316</td><td>25.284</td><td>7.335</td><td>-</td><td>-</td></tr></tbody></table>

## 3.2. Rolling Stock

The train used in both lines was an electrical multiple unit (EMU). The whole unit was 2.940 m wide 4.265 m high and had a total length of 98.05 m, with an unladen weight of 157.3 t. The units were composed of five cars, the two ends having a driver's cabin and a normal floor. The middle car had a normal floor, while the other two cars had a low floor. The five cars were supported on two types of bogies, the trailer bogie and the tractor. The tractor bogie was always shared between two cars. The train was designed to use a standard Iberian track gauge (1668 mm) at a maximum speed of 120 km/h with almost 1000 passengers, although it could reach 160 km/h with minor modifications. The maximum total power of the train was 2.2 MW, and it had regenerative braking. In the base case, the trains were not going to be equipped with on-board energy storage systems. However, we considered the possibility of adding an on-board accumulator system based on a hybrid battery/ultracapacitor technology. The power profiles obtained during the simulation in the accumulation system would help the manufacturer to determine the percentage of energy that must be stored in the battery or ultracap parts. The electromechanical efficiency of the trains in traction and braking modes ($v_t^{train}, v_b^{train}$), as well as the electrochemical efficiency of the storage system in charging and discharging modes ($v_c^{STO}, v_d^{STO}$) were set to 0.95. The total accumulation system capacity ($E_{max}^{STO}$) was 7 kWh, and the on-board energy storage device rated charging and discharging power ($P_{rated,c}^{STO}, P_{rated,d}^{STO}$) was 1 MW.

Regarding the protection curves of the trains and the storage elements, the minimum and the regulation voltage of the train in traction mode ($V_{min}^{train}, V_{reg,t}^{train}$) were set to 1980 V and 2280 V. The same values have been selected for the minimum and regulation voltage of the energy storage system ($V_{min}^{STO}, V_{reg}^{STO}$). In braking mode, the regulation voltage and the maximum voltage of the squeeze control ($V_{max}^{train}, V_{reg,b}^{train}$) were 3300 V and 3600 V, respectively. In the cases in which the on-board accumulation system was activated, the system would be initialized with no charge.

In Table 2, the data summarizing the behavior of the trains in the different trips are collected. In the first column, we can observe the required mechanical power to complete the trip. It can be seen that the slope of the blue line was steeper because the difference between the power required for outward and return journeys was greater than on the red line. The average trip considering the two lines and both directions needed 229 kWh. The mechanical regeneration capacity in Column 2 is the available mechanical power that can be regenerated. Columns 3 and 4 contain the required electrical power and the electrical regeneration capacity considering already the efficiency of the electromechanical
conversion. The electrical regeneration capacity was usually around 40% of the required electrical power, except in the S1 to S4 trip of the blue line. In this case, because the train ascended a steep slope, the regeneration capacity was much lower, around 22% of the required electrical power. In the fifth column, we can see the minimum electric consumption. This consumption was calculated as the required electrical power minus the electrical regeneration capacity. Off course, this is a theoretical consumption that considers that we took advantage of all electrical power available for regeneration. This is not true, mainly because of two reasons. First, part of the power that was available to be injected into the catenary was burned in the rheostatic braking system when the squeeze control was activated in order to maintain the catenary voltage below the maximum level. In addition, if the train was equipped with on-board accumulation, the efficiency of the electrochemical conversion during the charging and discharging process also reduced the percentage of available regenerated power that could be reused. For these reasons, we uses these minimum consumption figures as a theoretical ceiling to compare the different solutions, but we must be aware of the fact that we will not reach this theoretical ceiling.

<table><caption>**Table 2.** Summary of the train behavior in the different trips; all data are in kWh.</caption><thead><tr><th>Trip</th><th>Required Mechanical Energy</th><th>Mechanical Regeneration Capacity</th><th>Required Electrical Energy</th><th>Electrical Regen. Capacity</th><th>Min. Electrical Consump. Theoretical</th></tr></thead><tbody><tr><td>S1 to S6 Red</td><td>245</td><td>112</td><td>258</td><td>106</td><td>151</td></tr><tr><td>S6 to S1 Red</td><td>240</td><td>107</td><td>253</td><td>102</td><td>151</td></tr><tr><td>S1 to S4 Blue</td><td>243</td><td>61</td><td>256</td><td>58</td><td>198</td></tr><tr><td>S4 to S1 Blue</td><td>187</td><td>99</td><td>197</td><td>94</td><td>103</td></tr><tr><td>Average Trip</td><td>229</td><td>95</td><td>241</td><td>90</td><td>151</td></tr></tbody></table>

### 3.3. *Description of the Selected Scenarios*

The authors developed eight different scenarios to study the influence of the accumulation system on the network, as can be observed in Table 3. There were four different paths for the trains, from S1 to S6 and from S6 to S1 for the red line and from S1 to S4 and from S4 to S1 for the blue line. Four different traffic densities were considered. Light traffic scenarios used a train headset of 50 min with 10 departures for each of the above-described routes. The medium traffic scenario considered a train headset of 35 min with 14 departures per route. The dense traffic scenario launched 24 trains per route with a headset of 20 min. Finally, the heavy traffic scenario launched 47 trains per route with a headset of 10 min. The simulation interval was very similar for all scenarios, and it went from eight hours and 18 min for the light traffic scenario to eight hours and 28 min for the heavy traffic scenario. Each of the four traffic densities were simulated without and with on-board accumulation systems without modifying the feeding infrastructure. The obtained results are presented in the next section and will clarify the effect of the on-board accumulation systems on the railway traction networks depending on the density of the traffic.

<table><caption>**Table 3.** Summary of the different proposed scenarios.</caption><thead><tr><th>Scenario Code</th><th>Traffic Density</th><th>On-Board Acc.System</th><th>Trains Headset (min)</th><th>Trains per Route</th><th>Number of Trains</th><th>Sim.Time (h:mm)</th></tr></thead><tbody><tr><td>L0</td><td>Light</td><td>No</td><td>50</td><td>10</td><td>40</td><td>8:18</td></tr><tr><td>L1</td><td>Light</td><td>Yes</td><td>50</td><td>10</td><td>40</td><td>8:18</td></tr><tr><td>M0</td><td>Medium</td><td>No</td><td>35</td><td>14</td><td>56</td><td>8:23</td></tr><tr><td>M1</td><td>Medium</td><td>Yes</td><td>35</td><td>14</td><td>56</td><td>8:23</td></tr><tr><td>D0</td><td>Dense</td><td>No</td><td>20</td><td>24</td><td>96</td><td>8:28</td></tr><tr><td>D1</td><td>Dense</td><td>Yes</td><td>20</td><td>24</td><td>96</td><td>8:28</td></tr><tr><td>H0</td><td>Heavy</td><td>No</td><td>10</td><td>47</td><td>188</td><td>8:28</td></tr><tr><td>H1</td><td>Heavy</td><td>Yes</td><td>10</td><td>47</td><td>188</td><td>8:28</td></tr></tbody></table>
### 3.4. Analysis of the Results

The above-described scenarios were simulated, and the obtained results are analyzed in this subsection. With the software used for simulation, we can obtain time-varying series of each electrical variable in the train or in the feeding network. An example of this detailed analysis can be observed in Figure 4, in which Train Number 4 on the red line route starting from S1 is represented in Scenario L1. It must be noticed how on seven occasions, the train burned part of the regenerated energy using the rheostatic system due to the high catenary voltage, even when the train was equipped with an energy storage system and the storage capacity was not full. It can be observed also how the power extracted/injected into the catenary differed from the train mechanical reference due to the efficiency of the electromechanical conversion process, but also because part of the power was provided by the energy storage system. In Figure 5, the fourth train of the blue line starting from S1 is represented, in this case for the heavy traffic scenario (H1). It has to be remarked that the voltage level was much lower and the network receptivity higher. The burned power in this scenario was nearly zero since the overvoltage protection was only activated on four occasions with a very short duration.

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Detailed behavior of the fourth train on the red line from S1 to S6 in the light traffic scenario with on-board energy storage. (a) Mechanical reference power; power extracted from the catenary and burned in the rheostatic system. (b) Catenary voltage. (c) Storage system power (in black) and state of charge of the energy storage system (in red) in (p.u.).</figcaption></figure>

The analysis of the time-varying curves was very interesting since it showed how electrical variables were correlated and it helped to understand the mathematical model proposed in the previous sections. However, in order to determine the effectiveness of the combination of the trains plus the infrastructure, another kind of analysis must be carried out. In Figure 6 is represented the power extracted from the AC network at substation S5 of the red line during the first four hours of simulation for all scenarios. Each subfigure represents two scenarios with the same train headset, but with and without the accumulation system. The correlation between the power extracted from the AC network was higher in the low traffic scenarios. Figure 7 depicts the energy extracted from the AC network at substation S5 in the eight possible scenarios. It must be noticed that the red solid line represents always the scenarios with energy storage. Obviously, scenarios with the same traffic were quite correlated. It should be pointed out that in the light, medium, and dense traffic scenarios, the energy consumption from the AC network in the cases with energy storage systems was a little bit lower. Paradoxically, that is not the case for the heavy traffic scenario in which the energy consumed
from the AC network was higher when the trains were equipped with energy storage systems. As will be observed, in heavy traffic scenarios, it was more efficient to share the energy surplus with other trains than storing it in the on-board accumulation system.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Detailed behavior of the fourth train on the blue line from S1 to S4 in the heavy traffic scenario with on-board energy storage. (a) Mechanical reference power; power extracted from the catenary and burned in the rheostatic system. (b) Catenary voltage. (c) Storage system power (in black) and state of charge of the energy storage system (in red) in (p.u.).</figcaption></figure>

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Power in MW obtained from the AC grid by substation S5 on the red line for all scenarios. (a) Light Traffic Scenario; (b) Medium Traffic Scenario (c) Dense Traffic Scenario; (d) Heavy Traffic Scenario.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Energy in MWh obtained from the AC grid by substation S5 on the red line for all scenarios.</figcaption></figure>

In Figure 8, we represent the Marey diagrams of each scenario just for the red line. The horizontal axis represents time; in this case, we represented the first 80 min of simulation. In the vertical axis, we represent the position of the train. Solid black lines represent trains from S1 to S6, while dashed-dotted lines represent trains circulating from S6 to S1. Vertical red lines represent the instants at which all the substations were blocked at the same time. It must be noticed that in light, medium, and dense traffic scenarios, the percentage of instants in which all substations were blocked at the same time dropped drastically when we added on-board energy storage systems, improving this percentage by more than 10%. However, in the heavy traffic scenario, the percentage of blocking instants was already very low (5%) in the case without on-board accumulation. The installation of on-board accumulation in the heavy traffic scenario produced the blocking of all substations in only 0.8% of the cases.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Marey diagrams representing the first 80 min of the schedule of the red line in the different scenarios. Solid lines represent the trains moving from S1 to S6, while dashed-dotted lines represent the trains moving from S6 to S1. The vertical red lines mark the instants at which all the substations in the system were blocked due to the high regenerated power surplus. The number in the top-left corner of each subfigure represents the percentage of instants at which all the substations were blocked at the same time for the whole simulation interval. The scenario is indicated in the top-right corner of each subfigure. (a) Light Traffic without accumulation; (b) Medium traffic without accumulation; (c) Dense traffic without accumulation; (d) Heavy traffic without accumulation; (e) Light Traffic with accumulation; (f) Medium traffic with accumulation; (g) Dense traffic with accumulation; (h) Heavy traffic with accumulation</figcaption></figure>

In the next paragraphs, we will analyze the aggregated results obtained for the eight scenarios from three points of view. First, we will present a full summary of the system. Second, a summary from the point of view of the trains will be analyzed, and finally, we will show a summary from the point of view of the substations. The use of aggregated data help us to measure the impact of the on-board storage solution in different traffic scenarios for the same feeding infrastructure. In Table 4, we summarize the most representative energies in the system; all the numbers are in MWh. In the first row, we have the electrical energy required by the trains. This energy is going to be the reference
energy, since it only depends on the number of trains and does not depend on the equipment of the trains. This energy was the same with and without accumulation. For the light traffic scenario, the 40 trains required a total of 9.62 MWh. For the heavy scenario, the 188 trains required 45.2 MWh, the relation between the number of trains and the required electrical energy being linear. Something similar happened with the second and the third row, which represent the regeneration capacity and the minimum consumption, respectively. The regeneration capacity considers the electromechanical conversion efficiency, and the numbers obtained represent the electrical energy that could be used for injecting into the catenary or charging the storage system. The minimum consumption is a theoretical concept that cannot be reached since it is obtained by subtracting the regenerated capacity from the required electrical energy. Up to now, these numbers only depended on the traffic, and did not vary whether or not the trains had on-board energy storage equipment. The next row (fourth row) represents the energy demanded by the train at the catenary level. It differs from the electrical energy required by the trains for one reason in the case of trains without an accumulation system and two reasons in the case of trains with an on-board accumulation system. In the first case, the train demand was different from the electrical energy requirement because the over-current protection prevented the absorption of too much power in the case of low catenary voltage. In the second case, part of the electrical energy provided by the train can be provided by the on-board accumulation system.

<table><caption>**Table 4.** Summary of the trains and the network behavior in the different scenarios.</caption><thead><tr><th rowspan="2">Scenario</th><th colspan="8">Energy (MWh)</th></tr><tr><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>Req.Electrical</td><td>9.62</td><td>9.62</td><td>13.5</td><td>13.5</td><td>23.1</td><td>23.1</td><td>45.2</td><td>45.2</td></tr><tr><td>Reg.Capacity</td><td>3.60</td><td>3.60</td><td>5.04</td><td>5.04</td><td>8.64</td><td>8.64</td><td>16.9</td><td>16.2</td></tr><tr><td>Min.Comps.</td><td>6.03</td><td>6.03</td><td>8.44</td><td>8.44</td><td>14.4</td><td>14.4</td><td>28.3</td><td>28.3</td></tr><tr><td>Train Demand</td><td>8.92</td><td>7.23</td><td>12.5</td><td>10.1</td><td>21.8</td><td>17.3</td><td>43.9</td><td>34.8</td></tr><tr><td>Train Inject.</td><td>1.61</td><td>0.55</td><td>2.56</td><td>0.81</td><td>5.3</td><td>1.60</td><td>13.1</td><td>3.99</td></tr><tr><td>Trains Net</td><td>7.31</td><td>6.68</td><td>9.94</td><td>9.31</td><td>16.5</td><td>15.7</td><td>30.8</td><td>30.8</td></tr><tr><td>Prov.Subs</td><td>7.71</td><td>6.98</td><td>10.3</td><td>9.65</td><td>16.7</td><td>16.1</td><td>31.0</td><td>31.2</td></tr><tr><td>Rheostatic</td><td>2.00</td><td>0.49</td><td>2.45</td><td>0.61</td><td>3.36</td><td>0.79</td><td>3.85</td><td>0.84</td></tr><tr><td>Non Supp.</td><td>0.62</td><td>0.16</td><td>0.66</td><td>0.16</td><td>0.51</td><td>0.14</td><td>0.27</td><td>0.09</td></tr><tr><td>Grid Losses</td><td>0.40</td><td>0.31</td><td>0.37</td><td>0.34</td><td>0.19</td><td>0.37</td><td>0.17</td><td>0.40</td></tr></tbody></table>

As can be observed, for the light traffic scenario without accumulation, the trains absorbed 0.72 MWh less than the required energy, so we would have delays in the system. For the heavy traffic scenario, the difference was 1.3 MW. However, in order to compare the scenarios, it is better to compare percentages with respect to the required electrical energy. These percentages can be found in Table 5. The authors are fully aware that this second table is redundant, but we think that it is important to analyze at the same time the data in MWh and in percentage with respect to the required electrical energy. The train demand in the scenario L0 was 92.6% of the required electrical energy, which means that 6.4% of the required energy cannot be provided. In order to distinguish the correlations between the different variables with the different traffic scenarios with and without energy storage and extract conclusions about the trends of the energy savings and consumption, the energetic data in percentage with respect to the total energy demanded by the trains are also represented in Figure 9. As can be observed, we had the same trends as the traffic increased for the cases with and without energy storage. However, even when the trend was the same, the values were very different, and it was necessary to add two different y-scales in order to compare the plots. The real energy demanded by the trains increased with the traffic, as well as the energy injected by the trains. That means that with denser traffic, the number of instants in which the overcurrent and overvoltage protections were activated was lower. We also reduced the energy provided by the substations and the grid losses when we increased the traffic.
<table><caption>**Table 5.** Summary of the trains and the network in the different scenarios in (%).</caption><thead><tr><th colspan="9">Energy in % with Respect to the Electrical Energy Required by Trains</th></tr><tr><th></th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>Req. Electrical</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td><td>100</td></tr><tr><td>Reg. Capacity</td><td>37.4</td><td>37.4</td><td>37.4</td><td>37.4</td><td>37.4</td><td>37.4</td><td>37.4</td><td>37.4</td></tr><tr><td>Min. Comps.</td><td>62.6</td><td>62.6</td><td>62.6</td><td>62.6</td><td>62.6</td><td>62.6</td><td>62.6</td><td>62.6</td></tr><tr><td>Train Demand</td><td>92.6</td><td>75.1</td><td>92.7</td><td>75.1</td><td>94.4</td><td>75.1</td><td>97.1</td><td>77.0</td></tr><tr><td>Train Inject.</td><td>16.7</td><td>5.7</td><td>19.0</td><td>6.0</td><td>22.9</td><td>7.0</td><td>29.0</td><td>8.8</td></tr><tr><td>Trains Net</td><td>75.9</td><td>69.4</td><td>73.7</td><td>69.0</td><td>71.4</td><td>68.1</td><td>68.1</td><td>68.2</td></tr><tr><td>Prov. Subs</td><td>80.1</td><td>72.5</td><td>76.5</td><td>71.6</td><td>72.3</td><td>69.7</td><td>68.5</td><td>69.1</td></tr><tr><td>Rheostatic</td><td>20.8</td><td>5.1</td><td>18.2</td><td>4.6</td><td>14.5</td><td>3.4</td><td>8.5</td><td>1.9</td></tr><tr><td>Non Supp.</td><td>6.4</td><td>1.7</td><td>4.9</td><td>1.2</td><td>2.2</td><td>0.6</td><td>0.6</td><td>0.2</td></tr><tr><td>Grid Losses</td><td>4.2</td><td>3.2</td><td>2.7</td><td>2.5</td><td>0.8</td><td>1.6</td><td>0.4</td><td>0.9</td></tr></tbody></table>

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Energy comparison of the different scenarios. In order to compare scenarios with a very different number of trains, all the energy values have been expressed in percentage referring to the total energy demanded by the trains in each scenario. This value can be found for each scenario in the first row of the Table 4. Black bars represent the cases without energy storage, and they are associated with the y-scale on the left. Red bars represent the cases with energy storage, and they are associated with the y-scale on the right.</figcaption></figure>

The non-supplied energy in the case of the heavy traffic scenario without energy storage was only 0.6%. The non-supplied energy is represented in Row 9 of Tables 4 and 5. This non-supplied energy in percentage is a good index of the network congestion, as it can be observed that the more trains in the system, the less the network congestion. In general, we can state that the installation of on-board energy storage always reduced the amount of non-supplied energy. In the worst scenario (light traffic), the non-supplied energy was reduced from 6.4–1.7% when we added storage to the trains. Row 5 represents the energy regenerated by the trains that was actually injected into the feeding system. Obviously, this energy increased with the traffic, but it was always reduced when we added energy storage to the trains. In Row 6, the trains’ net energy was obtained subtracting Row 5 (energy injected) from Row 4 (energy demanded). We can see that in the light, medium, and dense traffic scenario, the net energy was always lower when the trains were equipped with energy storage. However, the net
energy was nearly the same in the scenarios H0 and H1. Row 7 represents the energy provided by the substations. It is very interesting to remark that even in the worst scenario (L0), the substations provided only 80% of the energy demanded by the trains, the rest being provided by other trains in braking mode. Again, for the light, medium, and dense traffic scenarios, the energy provided by the substations was reduced when we added the energy storage. In the heavy traffic scenario, this trend was inverted, and the energy provided by the substations was slightly higher when we added energy storage to the trains. This is because when we had many trains and they were very close to each other, it was more efficient to use the regenerated energy in other trains than the stored energy in the on-board accumulation system. This is coherent with the data obtained representing the losses in the feeding system (Row 10). In the heavy traffic scenario, the losses were higher when the trains were equipped with an energy storage system. Finally, Row 8 represents the energy burned in the rheostatic system of the train, in all cases, it suffered a significant reduction when we added the on-board accumulation system.

In Table 6, we can observe the trains' average behavior in the different trips and in all scenarios, as well as the average trip considering all possible routes in all scenarios. The table is split into four blocks. In the first one, we represent the total energy obtained from the catenary. It should be noticed that the slope of the blue line was steeper than the slope in the red line. That is the reason why the energy obtained from the catenary in the outward and return journeys was so different, as well as the energy injected into the catenary. It must be noted that for the average trip, the energy obtained from the catenary increased when we increased the traffic because the voltage was more stable (within the limits), and the non-supplied energy decreased. The minimum consumption (theoretical) for the average trip defined in Table 2 was 151 kWh, and this number was obtained subtracting the electrical regeneration capacity from the required electrical energy. The best scenario in terms of average net consumption was H0 (heavy traffic without accumulation), with a net consumption of 164 kWh, only 5.8% above the theoretical consumption. Adding accumulation to the heavy traffic scenario increased the average trip net consumption by 1 kWh. In the rest of the scenarios, adding accumulation always improved the average trip net consumption. For instance, in the light traffic scenario, the net consumption for the average trip was 183 kWh without accumulation and 167 kWh with accumulation. It is very important to remark that when we considered the net consumption in the average trip for comparing with the theoretical limit, we were not considering the non-supplied energy that was quite high in the light, medium, and dense traffic scenarios without accumulation. When we added accumulation to the system, the non-supplied energy in the average trip was always lower than 2% with respect to the required electrical energy for the average trip; and in the case of heavy traffic, lower than 1%.

Table 7 contains the energy analysis of the substations for all scenarios; it is split into three blocks. The first block represents the energy injected into the DC traction system per each substation and the total. In the second block, we can observe the energy injected by each substation per train trip. Finally, this number is compared with the minimum theoretical consumption in Block 3. Again, we can observe that in the scenario H0, the total energy injected into all substations was only 9.3% above the minimum consumption.

Regarding the performance, both models (decoupled and integral) were equivalent from the point of view of the equations, and they provided the same result. However, the time savings on average when simulating the proposed decoupled model were around 20%. On average, the number of iterations of the decoupled model compared to the integral one was higher, from 6.5 iterations (average) to 7.9 iterations (average) per instant. The iterations using the decoupled model were faster. The average time to complete an iteration with the decoupled model was 0.39 ms, while the integral model spent 0.55 ms.
<table><caption>**Table 6.** Summary of the train behavior in the average trip for all lines and in all scenarios.</caption><thead><tr><th colspan="10">Energy Obtained from the Catenary (kWh)</th></tr><tr><th>Trip</th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>S1 to S6 Red</td><td>238</td><td>183</td><td>238</td><td>183</td><td>245</td><td>183</td><td>249</td><td>186</td></tr><tr><td>S6 to S1 Red</td><td>243</td><td>173</td><td>243</td><td>173</td><td>245</td><td>171</td><td>251</td><td>174</td></tr><tr><td>S1 to S4 Blue</td><td>230</td><td>222</td><td>230</td><td>221</td><td>233</td><td>223</td><td>246</td><td>233</td></tr><tr><td>S4 to S1 Blue</td><td>181</td><td>145</td><td>182</td><td>146</td><td>186</td><td>146</td><td>189</td><td>149</td></tr><tr><td>Average trip</td><td>223</td><td>181</td><td>223</td><td>181</td><td>227</td><td>181</td><td>234</td><td>186</td></tr></tbody><thead><tr><th colspan="10">Energy Injected into the Catenary (kWh)</th></tr><tr><th>Trip</th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>S1 to S6 Red</td><td>57</td><td>16</td><td>60</td><td>17</td><td>67</td><td>18</td><td>85</td><td>23</td></tr><tr><td>S6 to S1 Red</td><td>38</td><td>10</td><td>49</td><td>10</td><td>60</td><td>11</td><td>70</td><td>11</td></tr><tr><td>S1 to S4 Blue</td><td>9</td><td>7</td><td>13</td><td>7</td><td>28</td><td>12</td><td>34</td><td>13</td></tr><tr><td>S4 to S1 Blue</td><td>57</td><td>22</td><td>61</td><td>24</td><td>66</td><td>26</td><td>90</td><td>38</td></tr><tr><td>Average trip</td><td>40</td><td>14</td><td>46</td><td>15</td><td>55</td><td>17</td><td>70</td><td>21</td></tr></tbody><thead><tr><th colspan="10">Energy Burned in the Rheostatic Braking System (kWh)</th></tr><tr><th>Trip</th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>S1 to S6 Red</td><td>49</td><td>12</td><td>45</td><td>11</td><td>40</td><td>8</td><td>22</td><td>4</td></tr><tr><td>S6 to S1 Red</td><td>64</td><td>6</td><td>52</td><td>5</td><td>42</td><td>4</td><td>31</td><td>4</td></tr><tr><td>S1 to S4 Blue</td><td>50</td><td>13</td><td>45</td><td>12</td><td>29</td><td>7</td><td>24</td><td>5</td></tr><tr><td>S4 to S1 Blue</td><td>37</td><td>18</td><td>33</td><td>16</td><td>29</td><td>14</td><td>5</td><td>5</td></tr><tr><td>Average trip</td><td>50</td><td>12</td><td>44</td><td>11</td><td>35</td><td>8</td><td>21</td><td>5</td></tr></tbody><thead><tr><th colspan="10">Non-Supplied Energy (kWh)</th></tr><tr><th>Trip</th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>S1 to S6 Red</td><td>19</td><td>6</td><td>19</td><td>5</td><td>12</td><td>5</td><td>8</td><td>4</td></tr><tr><td>S6 to S1 Red</td><td>9</td><td>2</td><td>9</td><td>2</td><td>7</td><td>2</td><td>2</td><td>1</td></tr><tr><td>S1 to S4 Blue</td><td>9</td><td>4</td><td>24</td><td>5</td><td>21</td><td>4</td><td>10</td><td>2</td></tr><tr><td>S4 to S1 Blue</td><td>25</td><td>4</td><td>14</td><td>4</td><td>11</td><td>3</td><td>7</td><td>2</td></tr><tr><td>Average trip</td><td>15</td><td>4</td><td>16</td><td>4</td><td>13</td><td>3</td><td>7</td><td>2</td></tr></tbody></table>

<table><caption>**Table 7.** Summary of the substations' behavior in all scenarios.</caption><thead><tr><th colspan="10">Energy Injected into the DC Traction System (MWh)</th></tr><tr><th>Substation</th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>S5 (Red line)</td><td>2.35</td><td>2.03</td><td>3.04</td><td>2.76</td><td>4.88</td><td>4.64</td><td>8.99</td><td>9.02</td></tr><tr><td>S3 (Red line)</td><td>3.14</td><td>2.83</td><td>4.22</td><td>3.93</td><td>6.77</td><td>6.54</td><td>12.2</td><td>12.7</td></tr><tr><td>S3 (Blue line)</td><td>2.20</td><td>2.12</td><td>3.04</td><td>2.95</td><td>5.05</td><td>4.93</td><td>9.76</td><td>9.53</td></tr><tr><td>Total</td><td>7.71</td><td>6.98</td><td>10.3</td><td>9.65</td><td>16.7</td><td>16.1</td><td>31.0</td><td>31.2</td></tr></tbody><thead><tr><th colspan="10">Energy Injected into the DC Traction System per Train Trip (kWh/trip)</th></tr><tr><th>Substation</th><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>S5 (Red line)</td><td>59</td><td>51</td><td>54</td><td>49</td><td>51</td><td>48</td><td>48</td><td>48</td></tr><tr><td>S3 (Red line)</td><td>79</td><td>71</td><td>75</td><td>70</td><td>71</td><td>68</td><td>65</td><td>68</td></tr><tr><td>S3 (Blue line)</td><td>55</td><td>53</td><td>54</td><td>53</td><td>53</td><td>51</td><td>52</td><td>51</td></tr><tr><td>Total</td><td>193</td><td>175</td><td>184</td><td>172</td><td>174</td><td>168</td><td>165</td><td>166</td></tr></tbody><thead><tr><th colspan="10">Energy Injected into the DC Traction System per Train Trip above the Minimum Theoretical Consumption in %</th></tr><tr><th>L0</th><th>L1</th><th>M0</th><th>M1</th><th>D0</th><th>D1</th><th>H0</th><th>H1</th></tr></thead><tbody><tr><td>27.8</td><td>15.8</td><td>22.1</td><td>14.3</td><td>15.3</td><td>11.3</td><td>9.3</td><td>10.3</td></tr></tbody></table>
## 4. Conclusions

A simplified model of a storage system for power flow purposes was presented and tested. The proposed approach considered the train and the accumulation system as different devices placed at the same point of the railway traction network. Both devices were coupled and interacted with each other and with the electrical network. The proposed mathematical model reduced the computational burden of the simulation when it was compared with models that considered the train and the storage in an integral way. The model was embedded in a power flow solver, and extensive results were provided.

Analyzing the results obtained, it can be stated that the behavior of the system when adding this kind of technology cannot be easily predicted. It has been observed how reducing the train headset and increasing the traffic in the system, the non-supplied energy was reduced, as well as the number of instants at which all substations were blocked. In general, we could state that adding on-board accumulation systems to the trains reduced the burned energy in the rheostatic braking system due to the squeeze control activation and also the non-supplied power. However, this effect was highly correlated with the train headset. In order to perform comparisons with different traffic levels, all the energetic data must be normalized using as a rated value the total energy demanded by all trains in each specific scenario. In the cases studied in this work, we observed how the congestion in the network was inversely correlated with the traffic. That means that increasing the traffic, we alleviated the network, and the voltage profile was lower, so the activation of the overvoltage protection was less frequent. We observed also that with the heavy traffic scenario, the trains were able to absorb nearly 98% of the demanded energy when the on-board accumulation was installed and 77% without accumulation. The energy regenerated reached nearly 30% in the heavy traffic scenario with accumulation, when in the same situation without accumulation, it was only 9%, half of the injected regenerated power in the light traffic scenario with accumulation. As a general conclusion, it could be stated that up to a reasonable level, the increasing the traffic can alleviate the system, making it more efficient in relative terms. However, each specific infrastructure can have different behaviors, and it is important to make these kinds of studies during the infrastructure planning stage in order to make the correct decisions.

**Author Contributions:** Conceptualization, P.A.; methodology, B.M.; software, I.E.-S. and B.M.; validation, C.M. and P.A.; formal analysis, P.A., B.M. and C.M.; investigation, P.A., B.M., I.E.-S., C.M.; resources, P.A.; data curation, B.M. and I.E.-S.; writing—original draft preparation, P.A.; writing—review and editing, P.A. and C.M.; visualization, I.E.-S.; supervision, P.A.; project administration, P.A.; funding acquisition, P.A.

**Funding:** This research was partially funded by CAF TE and the Spanish Ministry of Science, Innovation and Universities through the CDTI strategic program CIEN (National Business Research Consortiums) under the granted project ESTEFI - FUO-EM-155-15

**Acknowledgments:** The authors would like to thank CAFTurnkey & Engineering, especially Peru Bidaguren and Urtzi Armendariz for their support during the development of these models.

**Conflicts of Interest:** The authors declare no conflict of interest. CAF TE provided the infrastructure data to perform the study. The Spanish Ministry of Science, Innovation and Universities had no role in the design of the study; in the collection, analyzes, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## Glossary

### Inputs

$P_{mech,ref}^{train}$ Mechanical power provided by an external software package; it is positive in traction mode and negative in braking mode.

### Parameters

$v_t^{train}, v_b^{train}$ Electromechanical efficiency of the train in traction and braking modes, respectively.
$v_d^{STO}, v_c^{STO}$ Electrochemical efficiency of the storage system in discharging and charging modes, respectively.

$E_{reg,d}^{STO}, E_{reg,c}^{STO}$ Storage system regulation energy in discharging and charging modes, respectively, used in the functions $g_{d1}$ and $g_{c1}$.
<table><tr><td>E<sup>STO</sup><sub>max</sub></td><td>Maximum amount of energy that can be stored; see functions g<sub>d1</sub> and g<sub>c1</sub>.</td></tr><tr><td>V<sup>train</sup><sub>reg,t</sub>, V<sup>train</sup><sub>min</sub></td><td>Train overcurrent protection parameters in traction mode; see function f<sub>t1</sub>.</td></tr><tr><td>V<sup>train</sup><sub>reg,b</sub>, V<sup>train</sup><sub>max</sub></td><td>Parameters for setting up the squeeze control of the train in braking mode; see function f<sub>b1</sub>.</td></tr><tr><td>V<sup>STO</sup><sub>reg</sub>, V<sup>STO</sup><sub>min</sub></td><td>Storage system overcurrent protection parameters; see functions g<sub>c2</sub> and g<sub>d2</sub>.</td></tr></table>

## Functions

$f_{t1}, f_{b1}$ Functions for obtaining the actual power of the train in traction and braking mode, respectively; this is the real power absorbed from the grid.

$g_{d1}, g_{c1}$ Functions for obtaining the maximum amount of power that can be extracted from or injected into the storage system, respectively; this function relates the maximum power in the physical storage system depending on the energy level.

$g_{d2}, g_{c2}$ Functions for obtaining the available discharging and charging power in the storage system, respectively; this power considers the storage system's overall efficiency and the network conditions.

## Variables

<table><tr><td>E<sub>STO</sub></td><td>Energy in the storage system at a specific instant.</td></tr><tr><td>p<sup>train</sup></td><td>Actual power that the train exchanges with the electrical traction network.</td></tr><tr><td>p<sup>STO</sup></td><td>Actual power that the storage system exchanges with the electrical traction network.</td></tr><tr><td>p<sup>cat</sup></td><td>Actual power that the set train plus storage system exchange with the electrical traction network; it is positive if it is absorbed from the grid.</td></tr><tr><td>p<sup>train</sup><sub>ref</sub></td><td>Train electrical reference power (considering the electromechanical efficiency); this power is obtained directly using the efficiency coefficients from the electrical machines and the input mechanical power; this power includes also the auxiliary power, and it refers to the catenary level.</td></tr><tr><td>P<sup>STO</sup><sub>rated,c'</sub></td><td>Energy storage device rated power in charging and discharging modes, respectively.</td></tr><tr><td>P<sup>STO</sup><sub>max</sub></td><td>Maximum power at each instant that can be extracted from the storage system; this power is at the accumulation system side; it does not consider the electrochemical efficiency of the process.</td></tr><tr><td>P<sup>STO</sup><sub>max,av</sub></td><td>Maximum power at each instant that can be extracted from or injected into the storage system; it only depends on the energy stored and the electrochemical efficiency, not on the network conditions.</td></tr><tr><td>P<sup>STO</sup><sub>av</sub></td><td>Real available power that the storage system can inject into or extract from the system in discharging and charging mode.</td></tr><tr><td>V<sub>cat</sub></td><td>Catenary voltage.</td></tr></table>

## References

1. Jabr, R.A.; Dzafic, I. Solution of DC Railway Traction Power Flow Systems including Limited Network Receptivity. *IEEE Trans. Power Syst.* **2017**. [CrossRef]

2. Arboleya, P.; Mohamed, B.; Gonzalez-Moran, C.; El-Sayed, I. BFS Algorithm for Voltage-Constrained Meshed DC Traction Networks With Nonsmooth Voltage-Dependent Loads and Generators. *Power Syst. IEEE Trans.* **2016**, *31*, 1526–1536. [CrossRef]

3. Takagi, R.; Amano, T. Optimisation of reference state-of-charge curves for the feed-forward charge/discharge control of energy storage systems on-board DC electric railway vehicles. *IET Electrical Syst. Transp.* **2015**, *5*, 33–42. [CrossRef]

4. Abrahamsson, L.; Söder, L. Fast Estimation of Relations Between Aggregated Train Power System Data and Traffic Performance. *Veh. Technol. IEEE Trans.* **2011**, *60*, 16–29. [CrossRef]
5. López-López, A.; Pecharromán, R.; Fernández-Cardador, A.; Cucala, A. Smart traffic-scenario compressor for the efficient electrical simulation of mass transit systems. *Int. J. Electrical Power Energy Syst.* **2017**, *88*, 150–163.
[CrossRef]

6. Mayet, C.; Horrein, L.; Bouscayrol, A.; Delarue, P.; Verhille, J.N.; Chattot, E.; Lemaire-Semail, B. Comparison of Different Models and Simulation Approaches for the Energetic Study of a Subway. *Veh. Technol. IEEE Trans.* **2014**, *63*, 556–565.
[CrossRef]

7. Mayet, C.; Delarue, P.; Bouscayrol, A.; Chattot, E.; Verhille, J.N. Comparison of Different EMR-Based Models of Traction Power Substations for Energetic Studies of Subway Lines. *IEEE Trans. Veh. Technol.* **2016**, *65*, 1021–1029.
[CrossRef]

8. Mayet, C.; Bouscayrol, A.; Delarue, P.; Chattot, E.; Verhille, J.N. Electrokinematical Simulation for Flexible Energetic Studies of Railway Systems. *IEEE Trans. Ind. Electron.* **2018**, *65*, 3592–3600.
[CrossRef]

9. Mayet, C.; Delarue, P.; Bouscayrol, A.; Chattot, E. Hardware-In-the-Loop Simulation of Traction Power Supply for Power Flows Analysis of Multitrain Subway Lines. *IEEE Trans. Veh. Technol.* **2017**, *66*, 5564–5571.
[CrossRef]

10. Li, Q.; Wang, T.; Dai, C.; Chen, W.; Ma, L. Power Management Strategy Based on Adaptive Droop Control for a Fuel Cell-Battery-Supercapacitor Hybrid Tramway. *IEEE Trans. Veh. Technol.* **2018**, *67*, 5658–5670.
[CrossRef]

11. Mohamed, B.; Arboleya, P.; González-Morán, C. Modified Current Injection Method for Power Flow Analysis in Heavy-Meshed DC Railway Networks With Nonreversible Substations. *IEEE Trans. Veh. Technol.* **2017**, *66*, 7688–7696.
[CrossRef]

12. Arboleya, P.; Mohamed, B.; El-Sayed, I. DC Railway Simulation Including Controllable Power Electronic and Energy Storage Devices. *IEEE Trans. Power Syst.* **2018**, *33*, 5319–5329.
[CrossRef]

13. Gelman, V. Energy Storage That May Be Too Good to Be True: Comparison Between Wayside Storage and Reversible Thyristor Controlled Rectifiers for Heavy Rail. *Veh. Technol. Mag. IEEE* **2013**, *8*, 70–80.
[CrossRef]

14. Gelman, V. Braking energy recuperation. *Veh. Technol. Mag. IEEE* **2009**, *4*, 82–89.
[CrossRef]

15. Arboleya, P.; Bidaguren, P.; Armendariz, U. Energy Is On Board: Energy Storage and Other Alternatives in Modern Light Railways. *Electr. Mag. IEEE* **2016**, *4*, 30–41.
[CrossRef]

16. Khodaparastan, M.; Dutta, O.; Saleh, M.; Mohamed, A.A. Modeling and Simulation of DC Electric Rail Transit Systems With Wayside Energy Storage. *IEEE Trans. Veh. Technol.* **2019**, *68*, 2218–2228.
[CrossRef]

17. Liu, H.; Zhou, M.; Guo, X.; Zhang, Z.; Ning, B.; Tang, T. Timetable Optimization for Regenerative Energy Utilization in Subway Systems. *IEEE Trans. Intell. Transp. Syst.* **2019**.
[CrossRef]

18. Liu, P.; Yang, L.; Gao, Z.; Huang, Y.; Li, S.; Gao, Y. Energy-Efficient Train Timetable Optimization in the Subway System with Energy Storage Devices. *IEEE Trans. Intell. Transp. Syst.* **2018**, *19*, 3947–3963.
[CrossRef]

19. Mayet, C.; Pouget, J.; Bouscayrol, A.; Lhomme, W. Influence of an Energy Storage System on the Energy Consumption of a Diesel-Electric Locomotive. *Veh. Technol. IEEE Trans.* **2014**, *63*, 1032–1040.
[CrossRef]

20. De la Torre, S.; Sánchez-Racero, A.J.; Aguado, J.A.; Reyes, M.; Martínez, O. Optimal Sizing of Energy Storage for Regenerative Braking in Electric Railway Systems. *IEEE Trans. Power Syst.* **2015**, *30*, 1492–1500.
[CrossRef]

21. Razik, L.; Berr, N.; Khayyamim, S.; Ponci, F.; Monti, A. REM-S-Railway Energy Management in Real Rail Operation. *IEEE Trans. Veh. Technol.* **2018**.
[CrossRef]

22. Lu, S.; Hillmansen, S.; Roberts, C. A Power-Management Strategy for Multiple-Unit Railroad Vehicles. *Veh. Technol. IEEE Trans.* **2011**, *60*, 406–420.
[CrossRef]

23. Khayyam, S.; Berr, N.; Razik, L.; Fleck, M.; Ponci, F.; Monti, A. Railway System Energy Management Optimization Demonstrated at Offline and Online Case Studies. *IEEE Trans. Intell. Transp. Syst.* **2018**, *19*, 3570–3583.
[CrossRef]

24. Graber, G.; Calderaro, V.; Galdi, V.; Piccolo, A.; Lamedica, R.; Ruvio, A. Techno-economic Sizing of Auxiliary-Battery-Based Substations in DC Railway Systems. *IEEE Trans. Transp. Electr.* **2018**, *4*, 616–625.
[CrossRef]

25. Mohamed, B.; El-Sayed, I.; Arboleya, P. DC Railway Infrastructure Simulation Including Energy Storage and Controllable Substations. In Proceedings of the 2018 IEEE Vehicle Power and Propulsion Conference (VPPC), Chicago, IL, USA, 27–30 August 2018; pp. 1–6.

26. Gee, A.M.; Dunn, R.W. Analysis of Trackside Flywheel Energy Storage in Light Rail Systems. *IEEE Trans. Veh. Technol.* **2015**, *64*, 3858–3869.
[CrossRef]
27. Alfieri, L.; Battistelli, L.; Pagano, M. Impact on railway infrastructure of wayside energy storage systems for regenerative braking management: A case study on a real Italian railway infrastructure. *IET Electr. Syst. Transp.* 2019. 10.1049/iet-est.2019.0005. [CrossRef]

28. Arboleya, P.; Coto, M.; Gonzalez-Moran, C.; Arregui, R. On board accumulator model for power flow studies in DC traction networks. *Electric Power Syst. Res.* **2014**, *116*, 266-275. [CrossRef]

29. Gomez-Exposito, A.; Mauricio, J.; Maza-Ortega, J. VSC-based MVDC railway electrification system. *IEEE Trans. Power Deliv.* **2014**, *29*, 422-431. [CrossRef]

30. Bai, Y.; Cao, Y.; Yu, Z.; Ho, T.K.; Roberts, C.; Mao, B. Cooperative Control of Metro Trains to Minimize Net Energy Consumption. *IEEE Trans. Intell. Transp. Syst.* **2019**. [CrossRef]

31. Kleftakis, V.A.; Hatziargyriou, N.D. Optimal control of reversible substations and wayside storage devices for voltage stabilization and energy savings in metro railway networks. *IEEE Trans. Transp. Electr.* **2019**. [CrossRef]

32. Clerici, A.; Tironi, E.; Castelli-Dezza, F. Multiport Converters and ESS on 3-kV DC Railway Lines: Case Study for Braking Energy Savings. *IEEE Trans. Ind. Appl.* **2018**, *54*, 2740-2750. [CrossRef]

33. Zhao, J.; Gao, Y.; Guo, J.; Chu, L.; Burke, A.F. Cycle life testing of lithium batteries: The effect of load-leveling. *Int. J. Electrochem. Sci* **2018**, *13*, 1773-1786. [CrossRef]

34. Arboleya, P. Heterogeneous multiscale method for multirate railway traction systems analysis. *IEEE Trans. Intell. Transp. Syst.* **2017**, *18*, 2575-2580. [CrossRef]

<figure><img src="image_13.png" /></figure>

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