# PV Based Standalone DC -Micro Grid System for EV Charging Station with New GWO-ANFIS MPPTs under Partial Shading Conditions

Authors: R. Ragul; N. Shanmugasundaram; Mariaraja Paramasivam; Suresh Seetharaman; Sheela L. Mary Immaculate
Journal: International Transactions on Electrical Energy Systems
 (2023)

---

## Abstract

The goal of this article is to use MPPTs (maximum power point trackers) to extort maximum power from best configuration or combine renewable resources and energy storage systems that all work together in off-grid for electric vehicle charging. The grey wolf algorithm (GWO) searches the MPP at partial shading condition (PSC) with following two consideration one is high oscillations around GMPPs, and other is that they are unable to track the new GMPPs after it has changed positions because the seeking agents will be busy around the previous GMPPs captured. Hence, in this paper, the proposed research objective is to find solutions to these two difficulties. The issue of oscillations around GMPPs was handled by combining GWO with ANFISs (adaptive Neuro-Fuzzy inference system) to gently tune output produced power at GMPPs. ANFISs are distinguished by their near-zero oscillations and precise GMPPs capturing. The second issue called they are unable to track the new GMPPs after it has changed positions is addressed in this work by using novel initialization by GWOs (Grey wolf Optimizations). In the MATLAB-Simulink and experiments demonstrate the effectiveness of the suggested GWO-ANFIS MPPTs based off-grid station for EVs (Electrical Vehicle) battery charging.


---

## Body

## 1. Introduction

With the growing environmental concerns, governments throughout the world have established a variety of carbon emission limits. EVs that are powered entirely by electricity have the potential to replace gasoline and diesel vehicles. EVs enhance energy efficiency through effective fuel economy, especially when the power is generated from renewable sources such as solar and wind. When it comes to EVs in smart grids, that are more than simply an electrical burden, but also a power resource [1, 2], several proposals have been made [3–7] presented an approach for modelling EVs as an additional load on the distribution network. Several scholars [8, 9] proposed a grid-connected PV system for EVs charging. In most cases, a backup battery is not required for a grid-connected system. There are two configurations for the grid-connected system. One is that the PV system works in conjunction with the utility grid to provide electricity to the load. As a result, the output of PV system is frequently less than load’s power. The other is that the PV system is designed to provide the load with the needed power, and the utility grid serves as a backup source during solar power variations in this arrangement. With the ever-increasing demand for electric vehicles, it became imperative to expand the number of charging stations. EVs require a high number of charging stations. As a result, EVs will be a substantial burden on the future distribution network system. The phrase “off-grid” of PV systems are not linked to main grids is known as a “stand-alone” PV systems. These stand-alone PV systems are ideal for electrifying small places in remote areas or small towns. They can also be used as charging points for EVs on long-distance routes as depicted in Figure 1. Many publications [10, 11] have discussed off-grid system design.Figure 1
Off-grid charging station for EVs.Due to seasonal weather conditions and the non-linear nature of solar irradiation, solar power output varies significantly, necessitating a hybrid application or backup systems [12]. Solar irradiance is not strongly connected across surrounding locations over short time periods, which is a critical factor in PV power output losses and swings [13]. MPPTs are considered significant components of PV systems for increased efficiency. Because of the non-linear features of PV arrays, it is unavoidable to design an effective maximum power point tracking system which is not only effective and improves the solar power system’s output power [14]. Furthermore, a variety of internal and external factors including series/parallel resistances, Diode factors, solar irradiations, PV array surfaces, internal temperatures, shadows and dirt impact on PV system’s output efficiencies. Because the dynamics of PV systems under PSCs change over time, MPPTs for PV powered systems must track GMPPs smoothly and steadily in a variety of situations including shades, degraded PVs, changes in PV arrays or PV characteristics.A number of algorithms for MPPTs have been suggested for enhancing PV system’s effectiveness, including hill climbing [12], P&O (perturb and observe) [15–17], incremental conductance [18], however these techniques produce oscillations at MPPTs, resulting in power losses. IC technique [18] was developed to decrease these oscillations and increase module efficiency, however it only lowered the oscillations. P&O and IC approaches failed in time spans specified by changing atmospheric conditions. A scanning technique was presented [19–21] to estimate panel’s maximum power-delivering potential at any given operational condition for PV systems with rapidly changing and PSCs and insolations.Soft computing methods employed for MPPTs in PV systems include PSOs (Particle Swarm Optimizations) [22], GWOs [23], CS, ACOs (Ant Colony Optimizations) [24], KHOs (Krill Herd Optimizations), FAs (Firefly algorithms) [25], ABCs (artificial bee colonies), MVOs (Multi-Verse Optimizations), ALOs (Ant Lion Optimizations), SCAs (Sine Cosine Algorithms), Dragonfly Algorithms, WOAs (Whale Optimization Algorithms), MFOs (Moth-Flame Optimizations). Several studies [26, 27] have examined and introduced all of these approaches. In shaded or non-shaded circumstances where most approaches could catch the GMPPs. Nonetheless, they are affected by two major issues in PV applications. The first difficulty is with dynamic or time-varying GMPPs location in P-V curves, where most of these techniques may catch the initial GMPPs and stay around it, but not if the GMPP positions changes over a period of time. The second difficulty is that random variables associated with all of these soft computing techniques create substantial power swings near steady states. To overcome these two challenges, novel and efficient ways have been applied in this study where the first issue is overcome by introducing re-initialization methods for soft computing methodologies while the second issue is resolved by using combination of soft computing techniques namely GWOs and ANFISs and where quick peak tracking and low oscillations around GMPPs differentiate ANFISs.


## 2. Related Work

Controller settings were optimised and adaptively modified using a predictive neural networks controller by Mohamed et al. [28]. The scheme anticipated control parameters by tracking grid current and dc-bus voltage mean square errors and eliminating them in a specified amount of time.SMCs were proposed by Pahari and Subudhi [29] for control discontinuity based on high frequency conversions for driving closed loop systems to reach and stay on planned sliding surfaces. This approach considerably enhanced PV system’s tracking speeds, however the step size of modulation depths of switching devices impacted system’s dynamic and steady-state properties. When ∆U increased, Although tracking speeds increased, variations in PV array’s output powers and voltages also increased. Pahari and Subudhi introduced Integral sliding mode controls to increase controller’s steady state performances.External voltage control, P&O, and an adaptive integral differential slip films were used to design their new sliding surface in Kihal et al. [30] and where derivative and integral terms eliminated overshoots during quick solar irradiation changes and reduced steady-state fluctuations.The study in [31] was a novel model of adaptive PID controllers based on ANFISs proposed to address concerns in MPPTs approaches. It aided in maximising the output of DC pump in terms of speed. Their proposed controller was also used to test performances of EVs. To get the maximum power from solar-powered pump many methods have been tried and ANFISs have also been used to optimise performances of these intelligent systems. However, the approaches have certain flaws, necessitating the development of a new paradigm. As a result, in this research, a technique is developed in which a PID controller is used, and the combination of ANFISs and PID improves the performance of MPPTs. Moreover, this work’s proposed combinations were tested on EVs.Padmanaban et al. [32] presented grid-based PV systems with MPPT control mechanisms. Their work’s ABC method converted ANFIS membership functions and their experimental study demonstrated that PV grid integrations were dependable and safe.The author of [33] recommended modifying INC (Incremental Conductance) to track GMPPs of PV systems with PSCs without the application of any nature-inspired intelligence approaches like PSOs or other comparable algorithms. The work identified p-v curve as a mixture of areas and monitored global peaks by moving operating points from one region to another. In addition, the system incorporated varying sample times for quicker global peak tracking under extreme PSCs caused by unpredictable PV module shading patterns. Their PV module was modelled in combination with SEPIC converters for a battery charging applications using MATLAB/SIMULINK platform. The proposed technique was put to the test under various shading patterns to see how quickly and accurately their tracker tracked the power points. Finally, the suggested algorithm’s idea was tested using PILs (Process-in-Loops) on TMS320F28027 LAUNCHPAD DSP board. Their result comparisons with recently published work in this sector showed that their suggested approach was more beneficial in terms of processing times.Hence Evolutionary algorithms have reduced steady state oscillations and have advantage of ease in implementation and simplicity of computation.


## 3. Design Charging Station Converters

The power converter is the most important technology associated with PV systems. Maximum power from PV module is dragged from converter and sent to the load. This should be accomplished in grid-connected systems with least harmonic content in the current and a PF value greater than 0.9. The output voltage in off-grid systems should be adjusted to the necessary value. The suggested charging station in Figure1 employs three DC-DC converter topologies. The SEPIC converter converts the fluctuating voltage from modules to a stable voltage, while the buck converter reduces the DC Bus voltage to voltage of the electric vehicle battery. The bidirectional converter is use d to control the battery bank’s charging and discharging. The following sections detail the design of three DC-DC converters.


### 3.1. DC-DC SEPIC Converter

Figure2 depicts Single Ended Primary Inductance Converter. [34] contains PV modules’ input source voltage VPV, SEPIC inductor LSEPIC−1andLSEPIC−2, Coupling Capacitor CC, DC-link capacitor CDC, a diode DMPPT, MOSFET (Metal–Oxide–Semiconductor Field- Effect Transistor) SMPPT as a switch and the output voltage VDC (DC Bus Voltage). The SEPIC converter’s output voltage can reduce or enhance the input PV voltage. Controlling or changing voltage may be done by adjusting the MOSFET’s duty cycle.Figure 2
The SEPIC converter diagram.As per the equation (34), the output voltage VDC is proportionate to input voltage VPV.(1)DMPPT=VDCVPV+VDC,here DMPPT is duty cycle.Equation (1) is use d to determine duty cycle of SEPIC converter. Solar energy ranges between 0 and 1000 W/m2 and temperature ranges from 25 to 75 degrees Celsius in the simulation.SEPIC converter duty cycle is determined for a cell temperature of 25°C and irradiance of 1000 W/m2 in non-shaded conditions. VPV = 116 V, as you may have observed. When VDC is equal to 110 V, DMPPTMIN at 110 V is 0.48.If the two PV panels are shaded, the SEPIC Converter Duty Cycle is calculated at 25°C and 1000 W/m2 at cell temperature of 25°C. As observed, VPVmin = 58 V and VDC = 110 V, resulting in DMPPTMax at 110 V of 0.65. Then inductance for continuous current, Coupling Capacitor of SEPIC converter is calculated below


#### 3.1.1. Design of SEPIC Inductor

In the SEPIC converter, the inductance for continuous current is computed as follows [34]:(2)LSEPIC−1&LSEPIC−2=VPVminDMPPTMax∆IoMax∗FSW=1.58mH≅1.6mH.The switching frequency,FSW, has been set at 25 kHz. The ripple current ∆IoMax is computed as follows: VIN is the converter input voltage = 58 V. (taken as 13 percent of the total current):(3)∆ILBoost=0.13∗IPV=0.13∗7.34=0.95.


#### 3.1.2. Design of SEPIC Coupling Capacitor

Coupling capacitorCc ripple voltage is:(4)CC=IPV1−DMPPTMax∆VCC∗FSW=7.351−0.6512∗25=8.5uF≅10uF.


### 3.2. DC-DC Buck Converter

Buck converter illustrated in Figure3 [35] is made up of DC-link voltage as an input source VDC, buck inductor LBUCK, buck capacitor CBUCK, a diode (D), MOSFET (S) as a switch and the output voltage VEV as an EVs battery voltage (350 V in Tesla S 100D cars).Figure 3
The buck converter diagram.The duty ratio is calculated as follows:(5)DBUCK=VEVVDC=60110=0.545.


#### 3.2.1. Design of Buck Inductor and Buck Capacitor

The buck converter’s inductance for continuous current is computed as follows [35]:(6)LBUCK=VDC−VEV∆ILBUCK∗FSW=110−601.51∗10,000=3.3mH,∆ILBUCK (The ripple current) is computed as (taken as 13% of the charge current):(7)∆ILBUCK=0.13∗ICGB=0.013∗11.66=1.51A.


### 3.3. The Bidirectional Converter

The storage battery is supplied as a backup service during energy outages to smooth power storage produced by off-grid systems. As a consequence of the advantages listed below, this article recommends a bidirectional buck-boost converter, as illustrated in Figure4 [23], regulate storage and delivery of electricity among PV system and battery bank:(i)
It provides the most cost-effective solution with the fewest external components.(ii)
It accomplishes voltage step-up and step-down using the fewest possible components.(iii)
It has a decreased duty cycle when in operation.(iv)
It has a high efficiency over a wide range of input and output voltages.(v)
It is less costly than the majority of converters.Figure 4
The bidirectional converter diagram.Based on power flow direction, the bidirectional converter may be configured in two operation modes: boost and buck [36] describes these modalities in detail.


#### 3.3.1. Design for the Bidirectional Converter

The inductance (L BI) can be represented in two different ways:In Buck mode:(8)LBI−Buck=VDC−VSSV∆IBI−BUCK∗FSW=110−602.08∗10,000=2.4mH,∆IBI−Buck (ripple current) is computed as (taken as 13% of charge current):(9)∆IBI−BUCK=0.13∗ICGS=0.13∗16=2.08A.In Boost mode:(10)LBI−BOOST=VDC−VSSV∆IBI−BOOST∗FSW∗VDC=110−601.65∗10,000∗110=0.02mH,here ∆IBI (ripple current) under Boost mode is computed as (taken as 13% of the charge current):(11)∆IBI−BOOST=0.13∗ICGSVDCVSSV=0.13∗6.95∗11060=1.65A,LBI−Buck is more than LBI−BOOST to make sure that inductor current can function during continuous conduction mode (CCM) [36], LBI−BUCK value is selected, that is, equal to 2.4 mH here DBI−BOOST is duty cycle under Boost mode and computed as:(12)DBI−BOOST=1−VSSVVDC=1−60110=0.455.


## 3.1. DC-DC SEPIC Converter

Figure2 depicts Single Ended Primary Inductance Converter. [34] contains PV modules’ input source voltage VPV, SEPIC inductor LSEPIC−1andLSEPIC−2, Coupling Capacitor CC, DC-link capacitor CDC, a diode DMPPT, MOSFET (Metal–Oxide–Semiconductor Field- Effect Transistor) SMPPT as a switch and the output voltage VDC (DC Bus Voltage). The SEPIC converter’s output voltage can reduce or enhance the input PV voltage. Controlling or changing voltage may be done by adjusting the MOSFET’s duty cycle.Figure 2
The SEPIC converter diagram.As per the equation (34), the output voltage VDC is proportionate to input voltage VPV.(1)DMPPT=VDCVPV+VDC,here DMPPT is duty cycle.Equation (1) is use d to determine duty cycle of SEPIC converter. Solar energy ranges between 0 and 1000 W/m2 and temperature ranges from 25 to 75 degrees Celsius in the simulation.SEPIC converter duty cycle is determined for a cell temperature of 25°C and irradiance of 1000 W/m2 in non-shaded conditions. VPV = 116 V, as you may have observed. When VDC is equal to 110 V, DMPPTMIN at 110 V is 0.48.If the two PV panels are shaded, the SEPIC Converter Duty Cycle is calculated at 25°C and 1000 W/m2 at cell temperature of 25°C. As observed, VPVmin = 58 V and VDC = 110 V, resulting in DMPPTMax at 110 V of 0.65. Then inductance for continuous current, Coupling Capacitor of SEPIC converter is calculated below


### 3.1.1. Design of SEPIC Inductor

In the SEPIC converter, the inductance for continuous current is computed as follows [34]:(2)LSEPIC−1&LSEPIC−2=VPVminDMPPTMax∆IoMax∗FSW=1.58mH≅1.6mH.The switching frequency,FSW, has been set at 25 kHz. The ripple current ∆IoMax is computed as follows: VIN is the converter input voltage = 58 V. (taken as 13 percent of the total current):(3)∆ILBoost=0.13∗IPV=0.13∗7.34=0.95.


### 3.1.2. Design of SEPIC Coupling Capacitor

Coupling capacitorCc ripple voltage is:(4)CC=IPV1−DMPPTMax∆VCC∗FSW=7.351−0.6512∗25=8.5uF≅10uF.


## 3.1.1. Design of SEPIC Inductor

In the SEPIC converter, the inductance for continuous current is computed as follows [34]:(2)LSEPIC−1&LSEPIC−2=VPVminDMPPTMax∆IoMax∗FSW=1.58mH≅1.6mH.The switching frequency,FSW, has been set at 25 kHz. The ripple current ∆IoMax is computed as follows: VIN is the converter input voltage = 58 V. (taken as 13 percent of the total current):(3)∆ILBoost=0.13∗IPV=0.13∗7.34=0.95.


## 3.1.2. Design of SEPIC Coupling Capacitor

Coupling capacitorCc ripple voltage is:(4)CC=IPV1−DMPPTMax∆VCC∗FSW=7.351−0.6512∗25=8.5uF≅10uF.


## 3.2. DC-DC Buck Converter

Buck converter illustrated in Figure3 [35] is made up of DC-link voltage as an input source VDC, buck inductor LBUCK, buck capacitor CBUCK, a diode (D), MOSFET (S) as a switch and the output voltage VEV as an EVs battery voltage (350 V in Tesla S 100D cars).Figure 3
The buck converter diagram.The duty ratio is calculated as follows:(5)DBUCK=VEVVDC=60110=0.545.


### 3.2.1. Design of Buck Inductor and Buck Capacitor

The buck converter’s inductance for continuous current is computed as follows [35]:(6)LBUCK=VDC−VEV∆ILBUCK∗FSW=110−601.51∗10,000=3.3mH,∆ILBUCK (The ripple current) is computed as (taken as 13% of the charge current):(7)∆ILBUCK=0.13∗ICGB=0.013∗11.66=1.51A.


## 3.2.1. Design of Buck Inductor and Buck Capacitor

The buck converter’s inductance for continuous current is computed as follows [35]:(6)LBUCK=VDC−VEV∆ILBUCK∗FSW=110−601.51∗10,000=3.3mH,∆ILBUCK (The ripple current) is computed as (taken as 13% of the charge current):(7)∆ILBUCK=0.13∗ICGB=0.013∗11.66=1.51A.


## 3.3. The Bidirectional Converter

The storage battery is supplied as a backup service during energy outages to smooth power storage produced by off-grid systems. As a consequence of the advantages listed below, this article recommends a bidirectional buck-boost converter, as illustrated in Figure4 [23], regulate storage and delivery of electricity among PV system and battery bank:(i)
It provides the most cost-effective solution with the fewest external components.(ii)
It accomplishes voltage step-up and step-down using the fewest possible components.(iii)
It has a decreased duty cycle when in operation.(iv)
It has a high efficiency over a wide range of input and output voltages.(v)
It is less costly than the majority of converters.Figure 4
The bidirectional converter diagram.Based on power flow direction, the bidirectional converter may be configured in two operation modes: boost and buck [36] describes these modalities in detail.


### 3.3.1. Design for the Bidirectional Converter

The inductance (L BI) can be represented in two different ways:In Buck mode:(8)LBI−Buck=VDC−VSSV∆IBI−BUCK∗FSW=110−602.08∗10,000=2.4mH,∆IBI−Buck (ripple current) is computed as (taken as 13% of charge current):(9)∆IBI−BUCK=0.13∗ICGS=0.13∗16=2.08A.In Boost mode:(10)LBI−BOOST=VDC−VSSV∆IBI−BOOST∗FSW∗VDC=110−601.65∗10,000∗110=0.02mH,here ∆IBI (ripple current) under Boost mode is computed as (taken as 13% of the charge current):(11)∆IBI−BOOST=0.13∗ICGSVDCVSSV=0.13∗6.95∗11060=1.65A,LBI−Buck is more than LBI−BOOST to make sure that inductor current can function during continuous conduction mode (CCM) [36], LBI−BUCK value is selected, that is, equal to 2.4 mH here DBI−BOOST is duty cycle under Boost mode and computed as:(12)DBI−BOOST=1−VSSVVDC=1−60110=0.455.


## 3.3.1. Design for the Bidirectional Converter

The inductance (L BI) can be represented in two different ways:In Buck mode:(8)LBI−Buck=VDC−VSSV∆IBI−BUCK∗FSW=110−602.08∗10,000=2.4mH,∆IBI−Buck (ripple current) is computed as (taken as 13% of charge current):(9)∆IBI−BUCK=0.13∗ICGS=0.13∗16=2.08A.In Boost mode:(10)LBI−BOOST=VDC−VSSV∆IBI−BOOST∗FSW∗VDC=110−601.65∗10,000∗110=0.02mH,here ∆IBI (ripple current) under Boost mode is computed as (taken as 13% of the charge current):(11)∆IBI−BOOST=0.13∗ICGSVDCVSSV=0.13∗6.95∗11060=1.65A,LBI−Buck is more than LBI−BOOST to make sure that inductor current can function during continuous conduction mode (CCM) [36], LBI−BUCK value is selected, that is, equal to 2.4 mH here DBI−BOOST is duty cycle under Boost mode and computed as:(12)DBI−BOOST=1−VSSVVDC=1−60110=0.455.


## 4. Design of Closed Loop System

The Grey wolf Optimizer combined with ANFISs are implemented to track the new GMPPs after it has changed positions are discussed follows.


### 4.1. Proposed Grey Wolf Optimizer Combined with ANFISs Global Maximum Power Point Tracking Techniques

Due the PSC the PV array generates the several peaks is generated with various LMPPs and a single GMPPs so tracking of maximum power location is obtain by GWO algorithm. The GWO enhance power quality and ensure the efficiency of PVCs (photovoltaic cells) under PSCs (partial shading conditions) by varying the duty cycle with low iteration. Mirjalili et al. originally proposed Grey Wolf Optimizer in 2014 [23], and it is considered one of the most recent heuristic optimization techniques. This approach is based on the way grey wolves hunt food in the wild, pursuing, attacking, and killing them. Grey wolves prefer to live in packs of 5 to 10 individuals. They follow a strict social dominance structure with four levels of leadership. Leaders, named as alpha (α), and subleaders, which they call beta (β), etc, are referred to as leadership pyramids, as illustrated in Figure 5, where wolf dominance rises from top to bottom [37]. The alpha wolf holds duty cycle of maximum power by corresponding random search duty cycle set. Accordingly the second and third maximum powers are named as beta and delta respectively. The rest of maximum power solution at the duty cycle is denoted as omega. The tracking global maximum power behaviour of grey wolves, we suppose that the alpha, beta and delta have the better knowledge about the global maximum power. Therefore, we save the first to three best powers obtained so far iterations and oblige the next iteration search duty cycle to update the duty cycle according to the best maximum power. GWOs are created to replicate grey wolves behaviour in the optimization domain. Grey wolf leadership hierarchy is established by assuming the leaders; wolves known as alpha (α), subleaders known as beta (β), lower rank wolves known as delta (δ), lowest rank wolves known as omega (ω).Figure 5
Leadership pyramid with four level (α, β, δ and ω).Grey wolves encircle food in the hunt, as previously stated. The following equations are presented to analytically model encircling behaviour [23]:(13)E⟶=C⟶.DP⟶t−D⟶t,DP⟶t+1=DP⟶t−A⟶.E⟶,here t indicates current iteration, A, C are coefficient vectors whose values has capability to save balance between exploration and exploitation in searching region, DP is prey’s position vector, D represents grey wolf’s position vector. Two coefficient vectors A and C are computed like:(14)A⟶=2a⟶.r1⟶−a⟶,C=2.r2⟶,here, coefficient a is decreasing linearly beginning 2 to 0, r1, r1 are arbitrary vectors with value [1, 0]. Grey wolves hunt by circling their prey, and pack should follow the commands of alpha wolf Dα as a high priority, and orders of the beta wolves Dβ and delta wolves Dδ as a lesser priority. The following equations can be used to model this leadership hierarchy numerically:(15)Eα⟶=C1⟶.Dα⟶−D⟶,Eβ⟶=C2⟶.Dβ⟶−D⟶,Eδ⟶=C3⟶.Dδ⟶−D⟶,D1⟶=Dα⟶−A⟶1.Eα⟶,D2⟶=Dβ⟶−A⟶2.Eβ⟶,D3⟶=Dδ⟶−A⟶3.Eδ⟶,DP⟶t+1=D1⟶+D2⟶+D3⟶3.The solution’s exploration and exploitation are determined bya and A values, wherein A≤1, the wolves incline to exploitation (converge to prey) and A≥1, wolves tend to exploration (diverge from prey as it might be one of LMPPs). As stated in introduction, hybrid GWO-ANFIS is suggested to combine benefits of both approaches and monitor variant GMPPs under variation PSCs with almost negligible oscillations around GMPPs power. To monitor the variation GMPPs, GWOs with two re-initialization approaches is presented. GWOs is re-initialized to seek for GMPPs at first. After GWOs has caught initial GMPPs, it is ended depending on first condition and the role is given to ANFISs [38]. ANFISs are employed to soften the initial GMPPs captured by GWOs and track precise value of GMPPs with nearly minimal fluctuation around global power. Inputs of ANFISs are VPV, IPV and TPV [28], which create the Duty (DMPPT output, whereas the ANFISs output is the DC-DC converter’s optimum duty ratio (SEPIC converter) [39, 40]. In contrast, if the PSCs changes, ANFISs are terminated, and GWOs must be re-initialized to search for and follow new GMPPs, as illustrated in Figure 6.Figure 6
PV energy system with hybrid GWO-ANFIS based MPPTs.This study proposes a combination of GWOs and ANFISs for taking full advantage of both techniques while overcoming their disadvantages. When it comes to capturing GMPPs of PV systems, GWOs are the most effective and have the fastest convergences, but it also has a lot of oscillations around GMPPs at steady states. ANFISs which are opposite of GWOs remain at local peaks with extremely low oscillations around GMPPs. As a result, GWOs were first used to quickly and effectively track GMPPs and subsequently ANFISs then start to behave like MPPTs with very low oscillations. As a result, GWOs are employed to follow the GMPPs and avoid local peaks at the start of MPPTs. Furthermore, at GMPPs, the drawback of major oscillations is avoided by pausing operations and enabling ANFISs to operate with extremely low oscillations. In order for GWOs to be terminated and ANFISs to complete the control of MPPTs on GMPPs that have already been captured, the following prerequisites must be met:(16)Pnew−Pold≤E1,here, Pnew is current power obtained from PV system, Pold is pervious iteration value of simulation, E1 is permissible limit before transferring control from GWOs to ANFISs, that is equivalent to 2% of generated power. Condition for terminating ANFISs and go back to GWOs to re-initialize the agents to search for new GMPPs under new PSCs via following conditions:(17)Pnew−Pold≥E2,here E2 is permissible limit before transferring control from ANFISs to GWOs, that equivalent to 5% of generated power.There are two suggested initialization approaches for re-initializing the GWOs to begin looking for a new GMPPs. As previously stated, the goal of combining GWOs and ANFISs is to take advantages of both approaches, with GWOs being rapid and reliable MPPTs approach for tracking GMPPs under PSCs while ANFISs have lowest GMPP oscillations. GWOs will be unable to catch new GMPPs when they alters their position since all wolves will be looking in the prior GMPPs search zone. As a result, if the controller detects a change in PSCs, the agents of GWOs should be re-initialized to scan the search region for new GMPPs. ANFISs will follow GWOs once they have identified GMPPs until the controller detects another change in the PSCs.If the condition stated in equation (27) is true, the controller will follow the PSCs modification. Figure7 shows a flowchart which summarizes processes of hybrid GWO-ANFIS with PSCs modification re-initialized.Figure 7
Flowchart of GWO-ANFIS based on PSCs change.


### 4.2. The Charge and Discharge Control

Power is transferred from and to the battery bank via the charging and discharging controller. The bidirectional converter in the stand-alone PV system with the battery bank illustrated in Figure8 has four distinct control modes depending on magnitude of power produced by PV array, that is dependent on temperature and irradiance.Figure 8
Four operating modes flowchart.Mode 1.
When obtained PV power is less than power needed by EVs battery, i.e. PPV < PDDh, and the charge status of SOC is greater than 40%, first control mode is used. The bidirectional converter enters buck mode in this mode, where the controller sends a control signal to make switch 1 (S1) detached while also sending a control signal to make switch 2 (S2) connect. This mode uses storage batteries to offer additional power.Mode 2.
When obtained PV power is less than power needed by EVs battery, i.e. PPV < PDDh, and the charge status of SOC is less than 40%, the second control mode is used. In this mode, the controller sends a control signal to both switches 1 (S1) and 2 (S2), causing them to be disconnected.Mode 3.
When obtained PV power exceeds power needed by EVs battery, i.e., PPV > PDDh and the charge status of SOC is less than 90% and more than 40%, the third control mode is used. PV array not only supplies power to EVs battery under maximum power point enabled management in this mode, but it also charges the battery bank with surplus power produced by PV array. In this mode, the controller sends a control signal to switch 1 (S1) to connect and a control signal to switch 2 (S2) to disconnect, causing the bidirectional converter to enter boost mode.Mode 4.
When the generated PV power exceeds the EVs battery’s power demand, i.e., PPV > PDDh, and the storage battery’s charge state (SOC) is greater than 90%, the fourth control mode is used. The PV array delivers power to EVs battery under maximum power point enabled management in this mode, and charge status of storage battery is maintained at more than 90% by constant voltage charging to prevent the battery from discharging. In this mode, the controller sends a control signal to both switches 1 (S1) and 2 (S2), causing them to be disconnected. The four operational modes are depicted in Figure8’s flowchart.


## 4.1. Proposed Grey Wolf Optimizer Combined with ANFISs Global Maximum Power Point Tracking Techniques

Due the PSC the PV array generates the several peaks is generated with various LMPPs and a single GMPPs so tracking of maximum power location is obtain by GWO algorithm. The GWO enhance power quality and ensure the efficiency of PVCs (photovoltaic cells) under PSCs (partial shading conditions) by varying the duty cycle with low iteration. Mirjalili et al. originally proposed Grey Wolf Optimizer in 2014 [23], and it is considered one of the most recent heuristic optimization techniques. This approach is based on the way grey wolves hunt food in the wild, pursuing, attacking, and killing them. Grey wolves prefer to live in packs of 5 to 10 individuals. They follow a strict social dominance structure with four levels of leadership. Leaders, named as alpha (α), and subleaders, which they call beta (β), etc, are referred to as leadership pyramids, as illustrated in Figure 5, where wolf dominance rises from top to bottom [37]. The alpha wolf holds duty cycle of maximum power by corresponding random search duty cycle set. Accordingly the second and third maximum powers are named as beta and delta respectively. The rest of maximum power solution at the duty cycle is denoted as omega. The tracking global maximum power behaviour of grey wolves, we suppose that the alpha, beta and delta have the better knowledge about the global maximum power. Therefore, we save the first to three best powers obtained so far iterations and oblige the next iteration search duty cycle to update the duty cycle according to the best maximum power. GWOs are created to replicate grey wolves behaviour in the optimization domain. Grey wolf leadership hierarchy is established by assuming the leaders; wolves known as alpha (α), subleaders known as beta (β), lower rank wolves known as delta (δ), lowest rank wolves known as omega (ω).Figure 5
Leadership pyramid with four level (α, β, δ and ω).Grey wolves encircle food in the hunt, as previously stated. The following equations are presented to analytically model encircling behaviour [23]:(13)E⟶=C⟶.DP⟶t−D⟶t,DP⟶t+1=DP⟶t−A⟶.E⟶,here t indicates current iteration, A, C are coefficient vectors whose values has capability to save balance between exploration and exploitation in searching region, DP is prey’s position vector, D represents grey wolf’s position vector. Two coefficient vectors A and C are computed like:(14)A⟶=2a⟶.r1⟶−a⟶,C=2.r2⟶,here, coefficient a is decreasing linearly beginning 2 to 0, r1, r1 are arbitrary vectors with value [1, 0]. Grey wolves hunt by circling their prey, and pack should follow the commands of alpha wolf Dα as a high priority, and orders of the beta wolves Dβ and delta wolves Dδ as a lesser priority. The following equations can be used to model this leadership hierarchy numerically:(15)Eα⟶=C1⟶.Dα⟶−D⟶,Eβ⟶=C2⟶.Dβ⟶−D⟶,Eδ⟶=C3⟶.Dδ⟶−D⟶,D1⟶=Dα⟶−A⟶1.Eα⟶,D2⟶=Dβ⟶−A⟶2.Eβ⟶,D3⟶=Dδ⟶−A⟶3.Eδ⟶,DP⟶t+1=D1⟶+D2⟶+D3⟶3.The solution’s exploration and exploitation are determined bya and A values, wherein A≤1, the wolves incline to exploitation (converge to prey) and A≥1, wolves tend to exploration (diverge from prey as it might be one of LMPPs). As stated in introduction, hybrid GWO-ANFIS is suggested to combine benefits of both approaches and monitor variant GMPPs under variation PSCs with almost negligible oscillations around GMPPs power. To monitor the variation GMPPs, GWOs with two re-initialization approaches is presented. GWOs is re-initialized to seek for GMPPs at first. After GWOs has caught initial GMPPs, it is ended depending on first condition and the role is given to ANFISs [38]. ANFISs are employed to soften the initial GMPPs captured by GWOs and track precise value of GMPPs with nearly minimal fluctuation around global power. Inputs of ANFISs are VPV, IPV and TPV [28], which create the Duty (DMPPT output, whereas the ANFISs output is the DC-DC converter’s optimum duty ratio (SEPIC converter) [39, 40]. In contrast, if the PSCs changes, ANFISs are terminated, and GWOs must be re-initialized to search for and follow new GMPPs, as illustrated in Figure 6.Figure 6
PV energy system with hybrid GWO-ANFIS based MPPTs.This study proposes a combination of GWOs and ANFISs for taking full advantage of both techniques while overcoming their disadvantages. When it comes to capturing GMPPs of PV systems, GWOs are the most effective and have the fastest convergences, but it also has a lot of oscillations around GMPPs at steady states. ANFISs which are opposite of GWOs remain at local peaks with extremely low oscillations around GMPPs. As a result, GWOs were first used to quickly and effectively track GMPPs and subsequently ANFISs then start to behave like MPPTs with very low oscillations. As a result, GWOs are employed to follow the GMPPs and avoid local peaks at the start of MPPTs. Furthermore, at GMPPs, the drawback of major oscillations is avoided by pausing operations and enabling ANFISs to operate with extremely low oscillations. In order for GWOs to be terminated and ANFISs to complete the control of MPPTs on GMPPs that have already been captured, the following prerequisites must be met:(16)Pnew−Pold≤E1,here, Pnew is current power obtained from PV system, Pold is pervious iteration value of simulation, E1 is permissible limit before transferring control from GWOs to ANFISs, that is equivalent to 2% of generated power. Condition for terminating ANFISs and go back to GWOs to re-initialize the agents to search for new GMPPs under new PSCs via following conditions:(17)Pnew−Pold≥E2,here E2 is permissible limit before transferring control from ANFISs to GWOs, that equivalent to 5% of generated power.There are two suggested initialization approaches for re-initializing the GWOs to begin looking for a new GMPPs. As previously stated, the goal of combining GWOs and ANFISs is to take advantages of both approaches, with GWOs being rapid and reliable MPPTs approach for tracking GMPPs under PSCs while ANFISs have lowest GMPP oscillations. GWOs will be unable to catch new GMPPs when they alters their position since all wolves will be looking in the prior GMPPs search zone. As a result, if the controller detects a change in PSCs, the agents of GWOs should be re-initialized to scan the search region for new GMPPs. ANFISs will follow GWOs once they have identified GMPPs until the controller detects another change in the PSCs.If the condition stated in equation (27) is true, the controller will follow the PSCs modification. Figure7 shows a flowchart which summarizes processes of hybrid GWO-ANFIS with PSCs modification re-initialized.Figure 7
Flowchart of GWO-ANFIS based on PSCs change.


## 4.2. The Charge and Discharge Control

Power is transferred from and to the battery bank via the charging and discharging controller. The bidirectional converter in the stand-alone PV system with the battery bank illustrated in Figure8 has four distinct control modes depending on magnitude of power produced by PV array, that is dependent on temperature and irradiance.Figure 8
Four operating modes flowchart.Mode 1.
When obtained PV power is less than power needed by EVs battery, i.e. PPV < PDDh, and the charge status of SOC is greater than 40%, first control mode is used. The bidirectional converter enters buck mode in this mode, where the controller sends a control signal to make switch 1 (S1) detached while also sending a control signal to make switch 2 (S2) connect. This mode uses storage batteries to offer additional power.Mode 2.
When obtained PV power is less than power needed by EVs battery, i.e. PPV < PDDh, and the charge status of SOC is less than 40%, the second control mode is used. In this mode, the controller sends a control signal to both switches 1 (S1) and 2 (S2), causing them to be disconnected.Mode 3.
When obtained PV power exceeds power needed by EVs battery, i.e., PPV > PDDh and the charge status of SOC is less than 90% and more than 40%, the third control mode is used. PV array not only supplies power to EVs battery under maximum power point enabled management in this mode, but it also charges the battery bank with surplus power produced by PV array. In this mode, the controller sends a control signal to switch 1 (S1) to connect and a control signal to switch 2 (S2) to disconnect, causing the bidirectional converter to enter boost mode.Mode 4.
When the generated PV power exceeds the EVs battery’s power demand, i.e., PPV > PDDh, and the storage battery’s charge state (SOC) is greater than 90%, the fourth control mode is used. The PV array delivers power to EVs battery under maximum power point enabled management in this mode, and charge status of storage battery is maintained at more than 90% by constant voltage charging to prevent the battery from discharging. In this mode, the controller sends a control signal to both switches 1 (S1) and 2 (S2), causing them to be disconnected. The four operational modes are depicted in Figure8’s flowchart.


## 5. Simulation Result

Simulating the suggested TPC Converter for performance assessment with the Sim Power-System Toolbox is done in MATLAB/Simulink setting. DC-microgrid settings provided in Table1 were used to obtain all simulation results in MATLAB 2021.Table 1
Parameters used in simulation.
Sl.noParameterValue1PV MPP powerPMPP213.15 W2PV MPP voltageVMPP29 V3PV MPP current (IMPP7.35 A4Series PV panel4No’s5Parallel PV panel1No’s6Storage system battery typeLead-acid7Storage system battery nominal voltageVSSV60 V8Storage system battery rated capacity16 Ah9SEPIC converter inductorsLSEPIC−1&LSEPIC−21.6 mH10SEPIC converter coupling capacitorCC10μF11Buck converter inductorLBUCK3.3 mH12Bidirectional converter inductorLBI2.4 mH13DC link capacitorCDC60μFThe solar irradiation fluctuates as per the following pattern in the simulation, while all other parameters remain constant. As a consequence of GWO-ANFIS MPPTs, the DC-microgrid is at steady state att = 0 s, and the four PV at 1000 wb/m2 create the current IPV of 6.8 Amps, the PV output voltage VPV of 125.8 Volt, and the PV power PPV of 856 watts as seen in Figure 9. There is no bypass diode conduct since all four PVs are operated in equal irradiance with just one global peak, as illustrated in Figure 10. As indicated in Figure 11, the EVs charging power PEV from the DC micro grid is 700 Watts, with a nominal dc link voltage VEV of 110 V and an EVs current IEV of 6.36 Amps. The remaining PV power PESS of 150 Watts is stored in the ESS battery at a current IESS of −2.49 A and a voltage of 60.2 Volt, increasing the battery SOC from 40% to 60%, as illustrated in Figure 12.Figure 9
Simulation results (a) PV voltage (b) PV current (c) PV power.
(a)(b)(c)Figure 10
Simulation results (a) bypass diode PV -1 (b) bypass diode PV -2 (c) bypass diode PV -3 (d) bypass diode PV -4.
(a)(b)(c)(d)Figure 11
Simulation results (a) EVs voltage (b) EVs current (c) EVs power.
(a)(b)(c)Figure 12
Simulation results (a) ESS voltage (b) ESS current (c) ESS SOC (d) ESS power.
(a)(b)(c)(d)Due to partial shade, the irradiation of Panels 3 and 4 drops to 300wb/m2 at t = 0.2 s. PV output voltage VPV falls to 59.5 Volts, PV current IPV reduces to 7.1 amps, and PV power PPV falls to 422.5 Watts. As illustrated in Figure 10, this global power is tracked using the GWO-ANFIS MPPTs algorithm by removing the darkened panel at this instant bypass diode conduction. Even when the PV lighting varies, the DC bus voltage remains constant at 110 V thanks to the battery converter. At t = 0 s, the battery is charging, but as the PV output power decreases, it drains. Because the DC bus voltage remains constant at 110 V, EVs power P EVs is kept constant at 700 W, as illustrated in Figure 11. The remaining power PESS of 277.5 Watts obtained by PV is stored in ESS battery at a current IESS of 4.64 A and a voltage of 59.7 Volt, lowering the battery SOC from 40% as seen in Figure 12The PV panel and diesel generator are started with a produced power of 422.5 watts and 790 watts, respectively, as the ESS battery discharges below 40% SOC att = 0.4 sec, as illustrated in Figures 9 and 13. To keep the EVs power PEV at 700 W, the DC bus voltage is kept constant, as illustrated in Figure 11. He remaining power PESS of 512.5 watts charged in the ESS battery, increasing the SOC as illustrated in the Figure 12.Figure 13
Simulation results (a) diesel generator voltage (b) diesel generator ESS current (c) diesel generator power.
(a)(b)(c)


## 6. Experimental Results

The proposed GWO-ANFIS tracking algorithm has been evaluated for its performance with P&O MPP tracking algorithms. The evaluation tracked maximum power from PV array system to DC using a simple boost converter. Table2 lists hardware parameter values of 4S configuration in Proposed GWO-ANFIS MPPT.Table 2
Hardware parameter values of 4S configuration.
ParameterValueBoost inductance4S2.36 mHDC filter capacitor4S200μFLoad resistor4S140 ohmsMOSFET 400 V, 10 A, N-channelIRF740Schottky diode 200 V, 10 AMUR 10200Maximum power PVC currentImpp2.90 AMaximum power PVC voltageVmpp17.5 VMaximum power PVC powerPmpp50 WOpen circuit PVC voltageVoc21.8 VSort circuit PVC currentIsc3.20To validate the effectiveness of the proposed MPPT, experiments were carried out on real PV array for both 4S. To create partial shading, transparent sheets of different shapes were placed on PV modules.Figures14(a) and 14(b) shows the PV voltage, PV current and PV power of GWO-ANFIS based MPPT control at 4S configuration with two different Pattern-1 and Pattern-2 respectively. The corresponding values are denoted in Table 3. Figure 15. Shows the DC Voltage, DC current and DC power of Proposed GWO-ANFIS MPPT algoritham for Pattern2.Figure 14
Experimental results for proposed GWO-ANFIS MPPT method for 4S configuration. (a) Pattern 1. (b) Pattern 2.
(a)(b)Table 3
Performance comparison of the proposed GWO-ANFIS MPPT method for 4S configuration.
Shading patternGP in wattsVPVIPVPPV (W)Track efficiency (%)Pattern 1108.2 W/75.5 V75.71.41106.7498.65Pattern 276.33 W/35.8 V35.12.1374.7697.94Figure 15
Experimental result of proposed GWO-ANFIS MPPT for pattern 2. (a) DC link voltage. (b) DC link current. (c) DC link power.
(a)(b)(c)


## 7. Conclusion

This article successfully built and simulated a charging station that use s a stand-alone PV system with an energy storage system to charge an EVs battery. The GWO-ANFIS method has been used with the MPPT controller. The P-V curve has several peaks due to partial shading; one GMPPs and many LMPPs. Heuristic approaches such as GWOs can readily capture the GMPPs before it starts hunting. The results showed that GWO-ANFIS with PSCs change re-initialization is the optimum approach for tracking the dynamic GMPPs. The combination of the ANFISs controller and the GWOs method significantly reduced output power oscillations. The charge and discharge phases of the battery energy storage were also effectively modelled and simulated employing the bidirectional converter. A bidirectional converter logic controller has also been created and simulated. The logic controller’s efficiency with the bidirectional controller was demonstrated by the findings. The suggested station’s design and power management are discussed and evaluated in MATLAB/Simulink using three alternative modes of operation. The experimental verification is done for PV 4S configuration using GWO-ANFIS MPPT with Boost converter which proves that the maximum power is tracked in PSC at an efficiency 98.65%.



---

*Source: 2073742-2023-03-15.xml*