In this research, a stand-alone bifacial PV- and battery energy storage system (BESS)-fed EV charging system is developed, as shown in Figure 1, using an ANN-controlled MPPT method. A single-diode equivalent circuit is used to model the bifacial PV array, accounting for temperature, series resistance, shunt resistance, solar irradiance, and rear-side irradiance. The PV output is supplied to the EV charging station through a DC–DC boost converter, while the BESS supports DC-bus voltage regulation. The Secretary Bird Optimization Algorithm (SBOA) is used to optimize the PI controller gains and filter parameters. The modeling of the system components and the proposed methodology are described in the following subsections.
Modeling of Components
1. Bifacial PV system
The PV modules help to achieve the required current and voltage output. In the present study, two monofacial PV modules were used to emulate the energy generation characteristics of bifacial arrays under different surface-albedo conditions9. The bifacial gain of irradiance BGg is determined by Equation (1)
(1)
Here, front irradiation is Gf, and rear irradiance is Gr. Maximum bifacial PV current is Imp, and voltage is Vmp. Finally, the output produced by the bifacial system with the chosen number of series Ns and parallel Np modules are given by Equation (2). The control of PV is given in Figure 2
(2)
2. Battery storage system (BSS)
Batteries are a crucial means of storing electrical energy18 . This work uses Li-ion batteries as the energy storage system. Among the available rechargeable battery options, this work focuses on Li-ion batteries, which are considered optimal for EVs due to their efficiency and power density. Equations (3) and (4) describe the Li-ion battery with constraints.
(3)
The specifications of the PV, storage battery, and converter are summarized in Table 2, while Table 3 presents the power-flow allocation among the PV array, energy-storage unit, and connected loads.
(4)
3. DC –AC converter
The Neutral-Point-Clamped (NPC) inverter is widely adopted in modern power conversion systems for renewable energy. Compared with the conventional two-level structure, the NPC configuration offers reduced total harmonic distortion (THD), lower electromagnetic interference, and lower voltage stress on semiconductor devices. The DC link of this topology is divided into three potential levels +Vdc/2, 0, and -Vdc/2, by means of paired capacitors and clamping diodes forming the neutral junction. By appropriately controlling the switching, the inverter generates three discrete voltage states at the output terminals, thereby improving waveform quality and power conversion efficiency. The neutral point of the DC bus in an NPC inverter plays a key role in balancing capacitor voltages and ensuring proper inverter operation. The inverter's control scheme is shown in Figure 3.
4. DC-DC Boost Converter
DC-DC converters step up the input DC voltage to a higher output voltage by controlling the duty cycle of a high-frequency switching device while maintaining efficiency. It has two operating modes. During the ON state, the diode is reverse-biased, and the input source supplies energy to the inductor, which stores it as a magnetic field, while the output capacitor supplies energy to the load. During the OFF state, the diode becomes forward-biased. Here, energy stored in the inductor is released through the diode to the capacitor connected at the output and load, combining with the input source to produce an output greater than the input voltage. By varying the PWM duty cycle, the converter's output voltage varies according to the relation
, where D is the duty cycle. It is used in solar, fuel cell, EV, and storage systems.
ANN control scheme for MPPT
In this work, an ANN-based MPPT is selected to collect the maximum power generated from the solar system17. Irradiation and temperature datasets used for ANN training were generated in MATLAB/Simulink under varying environmental operating conditions, corresponding to solar irradiance levels of 800–1000 W/m2 and temperature ranges of 20–25 °C, with a duty cycle (D). The ANN is trained to minimize the MSE between the obtained and required outputs (Op,
) to extract maximum power from the bifacial PV system. Here, the solar irradiation and temperature are used as inputs to the ANN, and the duty cycle is the output to control the boost converter. The structure of the developed ANN model for MPPT, as given in Figure 4, is considered in the work. Supplementary Table 1 presents the advantages and disadvantages of the proposed method compared with other standard methods.
An ANN consists of three main layers: an output layer (OPL), an input layer (IPL), and a hidden layer (HIL). This helps to transfer the data between the IPL and HIL. Subsequently, it is produced by the weights in the links connecting the IPL and HIL. In this context, computations are executed with a certain bias applied to the HL variable, and the resulting outcomes are accumulated in the OPL variable. Here, the LMBP-type ANN [17] is selected. The connection weights are tuned throughout training by measuring the error to achieve the target output. Here, LMBP is used to train ANNs with the MSE as the performance function. The LMBP algorithm utilizes the derived derivatives to update the weights, which exhibit the properties of effective learning and accelerated convergence
Every neuron in a multilayer perceptron network has a summation and an activation function. Nevertheless, there are certain numerical weights (wpk) that connect these neurons across the levels. When inputs are multiplied by weights as specified in Equation (5). The nonlinear sigmoid function is considered with MSE minimization given in Equation (6).
(5)
(6)
SBOA optimized the filter and control parameters of the proposed system
These days, metaheuristic algorithms play a key role in solving engineering problems. Figure 5 gives a classification of algorithms. SBOA algorithm23 mimics the behavior of the SB to survive in its natural habitat. The exploration and exploitation stages of the Secretary Bird's hunting behavior are modeled in SBOA. The algorithm's exploration phase mimics the SB behavior of catching snakes, while its exploitation phase mimics their behavior of dodging predators like eagles. The solution process begins by generating random values within their respective limits for each SB in the population, and by calculating the objective function value by running the Simulink model, treating each SB's value as a design parameter. Based on the objective value, each SB changes its position to represent a better solution. The SBOA mathematical modeling is discussed below:
Initial preparation phase
First, it is necessary to identify the initial solutions that were utilized to begin the search for a typical minimization of the objective function F(Y). In this case, the initial random population Y = [Y1, Y2, ......YN] of SB‘s is formed via N initial solutions. Equation (7) initializes the population and Yi in the represents a solution.
(7)
Where the choice variables' lower and upper bounds are denoted by lb and ub. In [0, 1], r is a random number. N Is the problem's dimension. Furthermore, the solution's fitness value Fi = F(Yi) is used to quantify the quality of Yi.
Hunting strategy of SB’s
Finding, eating, and assaulting the prey are the three main phases of SB hunting. The entire hunting process has been divided into three equal time intervals, t < 1/3T, 1/3T < t < 2/3T, and 2/3T < t < T. Here, t is the current iteration, and T is the maximum number of iterations. These intervals correspond to the three phases of the SB's predation: searching, consuming, and attacking prey. These divisions are based on the biological statistics of the bird's hunting phases and the time durations of each phase. Consequently, the following is how each SBOA phase is modeled:
Searching for prey
The SB must hunt prey from a safe distance during this phase. To gather sufficient data across the entire search region, the first step in optimization algorithms requires stronger exploration. By using the locations of the other two SB’s as a guide, the SB can investigate additional possible regions. To further increase algorithm diversity, differential mutation processes are thus introduced. When t is less than 1/3rd of T, Equation (8) is used to update the position of each individual. It will search for, consume, and assault prey. Consequently, the following is how each SBOA phase is modeled:
(8)
Where, Yi(t) is the current position of the ith SB, Ynewi(t) gives the updated solution, two individuals Yi(t), Yr2(t) were randomly chosen from the existing population. A random vector, for the parameter-selection strategy, R1 has 1 × N randomly chosen from [0, 1] for effective balance among exploitation and exploration phases. The updated solution is considered only if it results in an improved objective value compared to the previously obtained solution.
Consuming prey
SB uses nimble agility and moves to hover around the snake after spotting possible prey. The prey's patience will be used to lower its defenses by observing and enticing opponents as they circle. Other SB’s adjust their posture near the prey. Hunting success will be significantly increased with this approach. Equation (9) illustrates how the Brownian motion (B) is used to simulate the SB’s random movement when 1/3 T < t < 2/3 T.
(9)
Where the B is implemented in MATLAB using the standard normal distribution by the randomly generated vector r(1,M). SB’s then use Equation (10) to update their positions.
(10)
Where, Ybest(t) is the best solution and the exponential scaling factor exp((t/T)4) is chosen to increase the exploitation capacity for every iteration.
Attacking prey
After consuming continuously, the victim will run out of energy. The secretary birds should launch the assault now. Here, different attack methods, such as incessant steps /sporadic long jumps over a short time period, are simulated using the Lévy flight technique. The characteristics are described by Equations (11) and (12). Candidate options are near the optimal answer at the moment because SB’s swiftly approach the prey. This plan will be implemented when t > 2/3T.
(11)
(12)
Where, LF represents the Levy flight strategy and
, s= 0.01, and φ = 1.5 fixed values. μ and υ are arbitrary values that lie in [0, 1]. Where, τ denotes gamma and η = 0.5, the weighting factor 0.5 was adopted, as proposed in the original SBOA formulation, to regulate the step size and maintain a balance between large exploratory jumps and local refinement. Table 4 presents the values used for the ANN and the SB algorithm in the developed system.
Escape strategy for SB
When hunting other prey in the wild, SBs risk being hunted. Eagles, hawks, foxes, and jackals are the primary adversaries they must contend with. They must use a variety of evasion techniques to safeguard themselves or their prey when they perceive danger. This algorithm simulates the escape methods by modeling running modes (D2) and camouflage (D1). When facing adversaries, secretary birds initially blend into their surroundings to stay safe. Here, secretary birds strive to avoid local optima in the algorithms by updating their positions around the prey (the best individual). They will use flight or rapid-running tactics to stay safe if they are unable to evade the enemy. For reference, a random individual Yrand is chosen as the leader to avoid being restricted to a local optimum. Equations (13) and (14) are used by secretary birds to update their positions. The SBOA flowchart is shown in Figure 6.
(13)
(14)
Representation of design variables
The problem variables in this study include the design parameters, such as the PI controller gains of DC-DC converters, as well as the resistances and inductances of filters and converters. Equation (15) gives the representation of control variables. The control parameter bonds are listed in Supplementary Table 2.
(15)
Fitness function (FF)
The minimization of THD is selected as the objective function (Obj) in this study. The maximum value of FF is the minimum value of Obj, given in Equation (16) and (17).
(16)
(17)