2.1. Single Battery Modelling
Batteries are at the core of the CS system and this component plays a strategic role in the EV field. Therefore, there is a plethora of studies on the topic [
7,
17,
18,
19]. Some of the technologies employed include: Molten Salt (Na-NiCl2), Nickel Metal Hydride (Ni-MH), Lithium Ion (Li-Ion) and Lithium Sulphur (Li-S). In the study undertaken by Iclodean et al. a comparison among different type of batteries applied on the same vehicle is considered. The most commonly used commercially is the Li-Ion battery [
7], which has a small “memory effect” that reduces the initial capacity over time. Even though the energy density is higher, the tests demonstrated that, for the same conditions, the Li-Ion has less autonomy (battery duration at a specific load level) than Ni-MH, which are the also widely used because of their high energy/power density. Li-Ion is currently the most convenient battery technology, even with the associated strong drawbacks related to the autonomy and temperature sensitivity. Further, due to the fact that 25–50% of the entire cost of an EV concerns the battery, it is essential to be cognisant of the associated implications. In the context of a potential fast transition to EV, it is compulsory to implement technologies that are economically competitive in order to incentivize such transition. Consequently, the battery type that has been taken into account in this work is the Li-Ion since it has a strong competitive price and it is the battery technology that applies most satisfactorily to the relevant (technology) transition considerations. The Li-Ion batteries can be controlled by adjusting: voltage, current, temperature and the load supplied by the battery. As suggested in [
20], this makes them suitable for electrical vehicles and grid applications. Generally, the most important functions that define the state of the battery are Equations (
1) and (
2):
Among the different typologies to model battery effects, the most suitable for EV’s applications is the Equivalent Circuit Method (ECM). This method uses lumped electrical components and controlled voltage sources in order to model the battery’s physics and chemical reactions, therefore it is possible use to electrical circuits solvers, in particular the one suggested by Tremblay in [
10]. This model has the peculiarity that requires just three points from the manufacturer’s discharge
-Voltage curve to obtain the parameters. In addition to this, the SIMULINK
TM software uses this model in the SimScape library as the battery. This model is characterized by two components: a fixed internal resistance
that models the complex physical and chemical reactions that causes voltages drop and heat (modelled as Joule Effect). The latter is an ideal controlled voltage generator that represents the
(Open Circuit Voltage) of the battery as presented in Equation (
3). An extraction of the battery modelled characteristic is contained in
Table 1, starting from the battery’s rated capacity and voltage. Zhang X et al. considers an ECM model that has some limitations [
21]. It is necessary to be aware of the hypothesis used by the model in order to not overestimate the results of the algorithm. For the work presented here it is important to consider how the voltage behaves during the charging pattern because it is focused only on the unidirectional smart charging process. In order to compute the actual energy that is stored in the battery, it is necessary to know the voltage at the terminals as it is possible to see in Equation (
1).
Li-Ion batteries are extremely sensitive to temperature [
22]. Frequent charge/discharge cycles also decrease the life cycle of the batteries [
23,
24]. In addition, the internal resistance changes over time when the internal structure of the electrodes changes. The effects are more intense when the batteries are operating at high C-rates (the unit used to measure the speed at which a battery is fully charged or discharged). To avoid the aforementioned adverse effects caused by the fast charging, it is common to restrict the battery’s operation at their nominal C-rate, however Amietszajew et al. proved that it is possible to use the cells in a wider C-rate range [
15]. Since the batteries of an EV could be solicited to employ higher currents (for example, an EB), it is necessary to know the
% range that could resist fast charging. The first safety limit is to impose a lower limit on the
, since the voltage will drop drastically when the
is close to the lowest tolerable value, for instance, an
of 20%. When the battery reaches the highest values of permitted
, the internal resistance
of the battery tends to grow, increasing the over voltage, Joule Effect would result in a degradation of the cathode. Therefore, a safety limit is imposed for a
of 80% to limit the temperature, lithium planting and extreme polarization (which causes a voltage spike when the
). This is why a high Constant-Current (CC) charging for an
of less than 80%, and a Constant-Voltage (CV) charging with
greater than 80% are adopted in fast charging techniques. Even if this is the most common charging protocol, it is also the least time-efficient. The possibility of charging the batteries at higher rates when
is suggested by Amietszajew. To achieve a safe charging pattern, it is not necessary to impose a gradual reduction of the C-Rate, as the current magnitude will follow the Load demand, which increases and decreases naturally and gradually. In addition, outside the periods where there is an under load demand (and for high currents to facilitate some compensation), the charging current is less than or equal to nominal.
By considering just the polarization effects of the Tremblay model, a test on the voltage increase of a Li-Ion battery pack with 650 V as rated voltage, 660 Ah capacity
Q and a charging current of 600 A was carried out. As the polarization effect is proportional to the current Equations (
4) and (
5), it was selected as the worst case in order to verify the discrepancy between the rated and distorted voltage.
The result appreciable in
Figure 1 shows that the divergence of this voltage during the charging pattern is not significant since it is less than 5 V (0.77%) in the range bounded between 20% and 80%, demonstrating that it is possible to use the rated voltage curve as a reference. In addition, there is the possibility to consider the voltage constant during the charging process as the error incurred can be neglected. For a more qualitative result, the algorithm employs a (MATLAB
TM) polyfitted Voltage characteristic that is a function of the
for more precise results when the battery is close to 80%.
2.2. Single Charger and Linearization
The charger choice has an important influence on how the batteries are employed. The first difference is the possibility to have either an internal charger (On-Board Charger, in practice a controlled rectifier) in the EBs or one that is accessed externally [
25]. This method is adopted primarily in private EVs that have to be charged in domestic environments by connecting the EV directly to a low voltage AC supply via a plug-in receptacle [
26]. Therefore, the capacity of the internal charges are limited to small powers/capacity and it is not possible to have external control because the vehicle self-modulates the power according to its charging protocol. For the purposes of this paper, there will be an emphasis on external chargers, which are primarily connected to a Charging Station (CS). This charger typology can facilitate more power, and use higher voltage levels to reduce the current magnitude, and as a result, the associated charging losses.
There are two possibilities associated with external chargers and their control typologies. The first employs a combination of an inverter and a DC-DC converter per vehicle. The alternative one employs a DC-DC per vehicle to connect all the converter primaries to a common DC-Busbar and then a single larger inverter connects the vehicles to the grid. The first method is mostly for isolated charging columns and has the peculiarity of a relatively small inverter size. The typology is useful when there is the possibility to combine a small size of renewable distributed generation, with an EV. By controlling the DC-DC and the EV’s inverter, it is possible to track the power from the RDG (for example a domestic PV generation) and compensate for the power produced to absorb more uniform power from the distribution network [
27]. The effectiveness of this method depends on the quality of the power tracking system and the coordination between all the vehicles connected to the specific node.
As the number of vehicles supplied by the same CS increases, the coordination of the DC-DC and inverters becomes more complex, and the second configuration, as depicted in
Figure 2, is preferable, where the RDG is connected to the AC side and the vehicles are all controlled with DC-DC converters and then connected to the AC grid with a single inverter. With the second configuration, it is possible to modulate the current absorbed by each EV by controlling the DC-DC converters and then by regulating the larger inverter capacity to the grid’s requirements. The coordination in this case is simpler, but it is necessary to rely on a more expensive inverter. Moreover, there is no redundancy if the inverter becomes out-of-service. Since the energy that an EV could exchange is limited, the service provided in this approach is more effective as “Power intensive” rather than “Energy Intensive”, to accommodate any power fluctuations. The instantaneous value of the power generated by the RDG has to be tracked and communicated rapidly in order to have an effective control over the devices. Even if studies like [
27], proved the effectiveness of the method, the technology required is not yet available and the DN is not ready to host proficiently this type of service. Therefore, since the goal of the present study is to focus on a technique that is easily accessible from a technological and management perspective, mono-directional power-flow is prioritized. In this regard, the choice of the DC-DC converter could be further discussed in order to obtain a more efficient and economical service. The bi-directionality of the half-bridge is possible with the introduction of more components rather than a DC-DC Buck or Boost [
28]. Thereby, the efficiency is sacrificed in order to obtain a bi-directional power flow, which is not a requirement for the algorithm proposed here. So a DC-DC Buck converter is preferable both in terms of its efficiency (cost, complexity) and operability in the context of the inevitable time available before a full transition involving the full range of communications control is available.
The other problems involving EVs and the CS, include the uncertainty of the EV’s parking pattern; plugging problems; different typologies and sizes of the batteries, etc. [
27]. These aspects add randomness to the EV’s charge planning and management. EBs have: higher battery capacity, same size and typology, scheduled mission. These characteristics lead to a relatively higher involved power, the elimination of the diversified approach based on vehicle’s model, and a reduction in the randomness of the charging pattern. In order to have simpler relationships that involve constant parameters in the relative time-frame, it is necessary to modify the behavior of the inputs in Equation (
6).
Since the batteries are usually used in CC mode in the fast charging range, Equation (
1) can be expressed without the voltage terms and if the current is considered constant with discrete time steps, it is possible to simplify the
expression as Equation (
7). Subsequently, since this
formulation is related to the [
Ah] unit of measurement, it is necessary to modify it in terms of energy in [
Wh] to facilitate a comparison with the data from the DAEM. To do so Equation (
2) and more specifically the spline poly-fitted function, is employed. Therefore, energy, as described in Equation (
6), becomes Equation (
8).
With a view to achieving a process linearity in terms of power transfer, as it facilitates the EBs to be operated separately while still having a unified end impact on the final CS’s power, it is more useful to exploit the analytical expression to facilitate a forecast of the actual energy on subsequent discrete time steps. With Equation (
9), it is possible to find the new
level and the linearity of the formulation is clear. The computation of the new charge level is more accessible in
form, as provided in Equation (
10).
The application of Equation (
9) has been verified by simulating the charging process and including the Buck Converter’s switching effect, with different constant currents to confirm the linearity with a model built in SIMULINK
TM. By controlling the current with different values, the power, as illustrated in
Figure 3, and the energy, as portrayed in
Figure 4, at the DC/DC Buck terminals, follows a relatively linear and ‘constant’ pattern that is proportional to the current’s magnitude including the additional implications caused by the Buck’s switching.
In conclusion, the approach advocates for linear relationships in the consideration of battery energy and . In achieving this approach, a more complex model that includes the polarization effects, can adopt a linear energy/power with a step-wise current pattern, if the battery operation is restricted within the 20–80% range. Consequently, it is possible to associate the current and power quantities through suitable vectors (developed in MATLABTM) and compute the current required. This approach facilitates the opportunity to test iteratively the energy to be provided to the EBs in the time step until the entire energy of the whole CS matches the one required to improve the DN demand’s shape.
2.3. The Functionality of the Smart-Charging Algorithm
Some methods for managing electric buses or charging stations can be found in the literature. Zhuang P. and Liang H. proposed a stochastic method [
29], Han B. et al. a method that uses timetables and routing as constraints [
30], and Hasan M.M. et al. a method that improves electric motor efficiency [
31]. Gkiotsalitisin K. emphasized the necessity of regular charging periods for EBs in lowering passenger travel time [
32]. Zhang C. proposed a strategy for overall optimization of EB scheduling [
33].
The algorithm presented in this study intends to propose a solution that may be employed in a smart grid transition by leveraging readily available data (DAEM, data-sheets, etc.) to construct a CS load profile and the EB scheduling, to achieve a given demand shape in a restricted DN. Furthermore, exact currents and EB schedules are prepared every 15 min for general CS management purposes.
On the other hand, the algorithm could be used to forecast the effects of a given CS and EB on different DN configurations (loads, PV, etc.). It allows for the photovoltaic (PV) power and the penetration of self-consumption (allotment of PV power saved in batteries during the day to be used during the night to reduce power absorbed by the grid ) to be altered. To simulate alternative scenarios, the quota of the PV installed in the DN and the self-consumption as a proportion of its rated installed power can be changed. As the generated PV power saved for local consumption is not part of the usual load, it reduces the PV curve magnitude, thereby mitigating the Duck Shape (a sinking impact of the demand shape during the noon hours caused by increased PV generation, resulting in lower grid energy demand ) [
34]. The load is further decreased during the hours of darkness, when demand is already low. Furthermore, the produced power can be adjusted from 0 to 100% of the rated PV output to simulate varied weather situations.
The first operation is to polyfit the original demand curve in order to obtain Equation (
11). The algorithm requires the target level for the relevant day, which could be partial peak- shaving or total load-levelling, and produces a vector containing the objective demand per time step,
. The
vector is in discrete form, to compare the actual demand with the objective, and Equation (
11) is discretized by computing the average power in a
time period, consistent with how DAEM is presented. Thereafter, the equation represents its vectorial form
.
Thereafter, the energy difference between the actual demand and the aim is estimated for each time step
i long
, resulting in the
, Equation (
12) with
elements. The power absorbed by the charging station will be already incorporated in
by taking into account each EB continuously absorbing the rated current
. During the
, however, the vehicles will absorb a “smart current”
based on their
% status. The algorithm applies
, the current to each EB
j at each time step by computing the new
with Equation (
9) and then translating it to [kWh] with Equation (
10) to acquire the energy
. To start the smart charging process, the initial current
is enforced to produce
. Then, in order to compute the difference in energy between the normal operation and the smart process, the full CS’s energy transition with all m EBs is expressed in Equation (
13). The EBs will be charged with the same
during the smart charging process, with the exception of vehicles that, if charged with
, would exceed the
. As a result, the smart current for these vehicles will be charged with a current
calculated using Equation (
14).
The algorithm will first calculate the energy for each EB in order to obtain Equation (
13) for the
i time step, and then the objective function Equation (
15) will be tested. More specifically, it is checked if the chosen current
achieves the correct reduction/increase in CS energy at time step
i. If the criterion is not met, the process is repeated by increasing
by 0.1A until the condition is fulfilled. The
generated by Equation (
9) is then assigned to
for the next time steps
, and the method is continued until the last time step
for each element of
.
Figure 5 illustrates that it is the rise/reduction of energy between each time step that affects the ultimate load form, and not the energy already stored in the vehicles. In
Figure 5 case,
>
with
<<
, therefore the vehicles start at a specific
represented in yellow, the orange increment represents the additional charge that should be achieved if the vehicles were charged with rated current
. Thus, the blue level represents, in this case, the one that is required in order to achieve the DN requirements by charging the EBs with the smart current
(the initial demand lower than the required one).
The smart current will have four cases as in Equation (
16). The case that could present difficulty is when the difference between the demand and the objective power is greater than the one that achieves the proper CS (
>>
). In such a case the algorithm imposes
, so it turns off the entire CS and reduces the load as far as the CS can.
The
<<
example could take a long time to compute because the current only grows by 0.1A per iteration until it reaches the correct value or the maximum eligible current. To improve this issue, an adaptive current increment could be employed. The last unmarked instance with
might be implemented with a DC/DC Half-Bridge and is thus ignored. As a result, the currents (
,
) and
% are saved in Equation (
17) for each iteration at step I for each EB j. Vehicles that achieve
at time step
i are then connected with vehicles that achieve
at
, which is considered the best situation.
The tracking of vehicles that reach
is recorded in a vector
Equation (
18) that contains the number of EBs that have reached the limit,
. The entire fleet is then computed using Equation (
19) to achieve perfect load levelling. The
vector is used to create the time slots represented by the bars in
Figure 6, where the vehicles must be connected in order to have the sufficient power in the CS for the DN’s load-levelling. There is a need to have more EBs available to connect during the time steps where the number of ideal exchanged vehicles is higher. In general, Equation (
19) will generate numbers that have no physical or economic meaning (e.g., 70 EBs in the CS with a total fleet of 1000 EBs), hence the value is normalized by inputting the actual size of the available fleet. At this point, the method assigns the real number of EBs uniformly in the
slots and generates the
Equation (
20) vector, which will include
EBs as elements, with the remainder set automatically in idle mode (
) by Equation (
14).
In reality, some EBs with
will remain in the CS as there will be a practical constraint in the limitation on the number of EB available to exchange; however, vehicles in idle mode (with
) will be able to begin their mission. In other words, it is not necessary to have
m vehicles in the CS at all times. This vector might be used to build the CS timetables, and it is possible to see graphically which time steps are most relevant to actuate the exchange, in order to synchronise the mission with the electrical demand needs (as illustrated in
Figure 6). The blue bars in
Figure 6 appear only in some time steps and show the number of vehicles
that must be connected at time step
i in order to achieve an appropriate load-leveling (within the CS’s power restrictions).
Furthermore, in order to avoid impractical connection patterns, the starting
value must be accurately adjusted. If all of the EBs are set to
with a
of 20%, the connecting pattern will look like
Figure 6, where the vehicles reach the
value at once as it is possible to notice from the few high blue bars.
However, a uniform distribution of the initial
in the range of [20–80%] will yield a more realistic connection pattern (as illustrated in
Figure 7) and produce an accessible connection pattern, as it possible to notice from the number and the height of the blue bars. Since the EBs are generally all charged with the same
, by imposing this initial condition, the number of the disconnected or connected ones are well distributed over 24 h.
The ultimate mutual effects were computed in power terms and were compared to the starting power demand. Equation (
21) computes the power
of the CS when the EBs absorb
, and Equation (
22) computes
, the total power absorbed by the CS with
. The final load was then calculated using Equation (
23) and plotted to ensure the accuracy of the load levelling/peak shaving.
Holistically,
Figure 8 illustrates a representation of the whole algorithm in operation.
Color maps are another type of output produced by the method by utilizing the data stored in Equation (
17). They are useful for providing a graphical global picture of the entire charging operation during a 24-h period. Looking at the current’s color gradient can provide information on the DN’s demand shape. The maps are generated from two matrices of size
, which may be transposed into EXCEL
TM tables to generate a charging pattern plan for each EB at each time step. A graphical depiction is required to have a global view of what is going on in the CS. Furthermore, the amount of the absorbed current can be shown to represent the influence on the network. The
matrix in the ideal scenario and the
, with the charging current Is and SoC level (with the color) of each vehicle
j at time step
i.
Figure 9 depicts the illustration of Equation (
14) when the EBs are nearly fully charged. For example, at the row corresponding to 4am, it is possible to see that the charging current is 140 A for most of the EBs. However, for the vehicles connected to the slots 20 to 50, the current color is different, therefore those vehicles are charged with a lower current (from 20 A to 120 A) because they are close to the full charged state. That is, for the precise vehicle
j that will exceed the 80%
, the charging current will be tailored as
to achieve the exact full charge. The following square has the same color as the others since it represents a newly connected (discharged) EB that charges properly with
. Aside from the nearly charged vehicles, all of the others have the same charging current, which is why the color seems uniform during time step
i. It is possible to note that the color does not vary monotonously during the day, which is due to the network requirements: when
, the indicator turns red (high
) and progressively turns blue (low/zero
).
Since the current is represented in
Figure 9 as a percentage of the
, there are some specific cases:
On the other hand, the
map
Figure 10 can be used to depict the connection spots and the charging process velocity. When the gradient takes longer to transition from blue to red, the CS is absorbing less current, and vice versa. The map is handy for visually verifying the amount of time that the EBs would require to be fully charged by looking at the
between the top blue square and the last red square.