Multi-Objective Grasshopper Optimization Based MPPT and VSC Control of Grid-Tied PV-Battery System

: This article presents the control of a three-phase three-wire (3P-3W) dual-stage grid-tied PV-battery storage system using a multi-objective grass-hopper optimization (MOGHO) algorithm. The voltage source converter (VSC) control of the presented system is implemented with adaptive kernel width sixth-order maximum correntropy criteria (AKWSOMCC) and maximum power point tracking (MPPT) control is accomplished using the variable step-size incremental conductance (VSS-InC) technique. The proposed VSC control offers lower mean square error and better accuracy, convergence rate and speed as compared to peer adaptive algorithms, i.e., least mean square (LMS), least mean fourth (LMF), maximum correntropy criteria (MCC), etc. The adaptive Gaussian kernel width is a function of the error signal, which changes to accommodate and filter Gaussian and non-Gaussian noise signals in each iteration. The VSS-InC based MPPT is provided with a MOGHO based modulation factor for better and faster tracking of the maximum power point during changing solar irradiation. Similarly, an optimized gain conventional PI controller regulates the DC bus to improve the power quality, and DC link stability during dynamic conditions. The optimized DC-link generates an accurate loss component of current, which further improves the VSC capability of fundamental load current component extraction. The VSC is designed to perform multi-functional operations, i.e., harmonics elimination, reactive power compensation, load balancing and power balancing at point of common coupling during diverse dynamic conditions. The MOSHO based VSS-InC, and DC bus performance is compared to particle swarm optimization (PSO) and genetic algorithm (GA). The proposed system operates satisfactorily as per IEEE519 standards in the MATLAB simulation environment.


Introduction
The grid-tied Photovoltaic (PV) system has become a natural choice for green energy, considering PV power's sharp cost curtailment [1]. The grid-tied PV systems have shown considerable commitment as a stable and reliable participant in the modern power grid [2,3].
Maximum power point tracking (MPPT) is crucial for achieving the desired efficiency and payback time of PV systems [4]. In [5,6], a detailed review of widely used traditional, intelligent and optimization based MPPT methods have been provided. In [7], an exhaustive review of artificial neural network (ANN) based MPPT has been provided along with ANN-based variable step-size MPPT is also proposed. The perturb and observe (P&O) and incremental conductance (InC) based MPPTs are among the widely implemented techniques due to their ease of implementation, reduced computational burden and adequate Numerous VSC control algorithms of various domains such as time domain, frequency domain, adaptive controls etc. for the grid-tied PV system have been investigated and implemented by the researchers. The synchronous reference frame (SRF), power balance theory (PBT), and many more [26] are among the widely used time-domain techniques, which are easy to implement and have less computational burden but produce more steadystate error. In [27], active reactive power control of instantaneous and average nature (IARC and AARC), balanced positive sequence compensation (BPSC), positive-negative sequence compensation (PNSC), etc. based on VSC controls have been proposed. The frequencydomain controls, i.e., Kalman filter-based controls (KFC) [28], Stockwell transformation (ST), etc. [29] have higher computational complexity and slow convergence speed. On the other hand, adaptive controls have a moderate computational burden and deliver reduced steady-state error with a faster convergence speed. Adaptive algorithms such as least mean square (LMS) [30], least mean fourth (LMF) [31], maximum correntropy criteria (MCC) [32] deliver faster convergence and robustness. The LMS control due to the fixed step size, produce a relatively higher steady-state error as compared to the LMF, which acts as a higher-order filter for the error signal. The MCC utilizes the Gaussian kernel which is effective against both Gaussian and non-Gaussian noise signals. Many advanced versions of LMS, LMF and MCC, i.e., hyperbolic cosine LMS (HCLMS) [33], hyperbolic tangent function (HTF) [34], LMS-LMF [35], adaptive neuro-fuzzy LMS (ANF-LMS) [36], zero attracting normalized LMF (ZAN-LMF) [37], normalized MCC (NMCC) [38], fixed forward prediction (FFP) [39], adaptive kernel width sixth-order MCC (AKWSOMCC) [40], etc. The HC-LMS offers minimal dynamic oscillations, HTF acts as a higher-order filter, which offers reduced mean square error (MSE), combined LMS-LMF reduces the shortcomings of LMS control by switching between second and fourth-order filter, but the performance of both controls reduces during heavily tailed Gaussian noises due to fixed step-size. The ANF-LMS delivers neuro-fuzzy tuned step-size for enhanced performance but increases the computational complexity. The ZAN-LMF, NMCC and FPP vouch for faster convergence rate and speed but have a high steady-state error due to fixed step size. The AKWSOMCC can handle Gaussian and non-Gaussian noises efficiently, while acting as a higher-order filter and delivering a better fundamental load current component. The function approximation and noise cancellation offered by AKWSOMCC are better than its peer adaptive controls due to adaptive kernel width, which also improves the convergence speed despite computational complexity. The proposed AKWSOMCC algorithms parameters are carefully tuned depending on their application ensuring the desired performance during the diverse dynamic condition. The optimized MPPT and DC bus further enhance the capability of the proposed control.
In the presented work, the MOGHO algorithm fulfils two objectives by providing the optimum δ for VSS-InC MPPT and PI controller gains (k p , k i ) for the V dc control, in a three-phase three-wire (3P-3W) grid-tied dual-stage PV-battery storage system. The VSC is controlled by the AKWSOMCC algorithm. During steady-state and dynamic states, the VSC control performs harmonics suppression, reactive power compensation, load balancing and active power balancing at the point of common coupling (PCC). The battery storage system maintains the power balance during dynamic conditions. The optimized DC bus generates an accurate loss current component, which further enhances the system's performance. The main attributes of the proposed work are as follows: a.
MOGHO based scaling factor (δ) optimization for VSS-InC MPPT to achieve faster tracking and reduced power oscillations at MPP. b.
MOGHO based PI controller gains (k p , k i ) optimization for the V dc control and generation of an accurate loss current component (i Loss ). c.
AKWSOMCC based VSC control for fundamental load current component extraction.

System Description
The proposed topology is presented in Figure 1. A dual-stage PV system of 33 kW along with battery storage is integrated with the 3P-3W grid. The PV array and battery design parameters are provided in Appendix A. The MOGHO algorithm optimizes the δ of VSS-InC based MPPT and PI controller gains k p &k i for V dc regulation. The boost converter is controlled by the gating signal delivered by the MPPT to stabilize the V PV . The battery current control during the charging and discharging phase is aided by the bi-directional buck-boost converter. The coupling capacitor at the DC link connects the DC and AC systems. The 3P-3W converter is attached to the PCC through interfacing inductors and ripple filters. The non-linear load of 16.5 kW and distribution grid of 415 V at 50 Hz are also attached to PCC.

System Description
The proposed topology is presented in Figure 1. A dual-stage PV system of 33 kW along with battery storage is integrated with the 3P-3W grid. The PV array and battery design parameters are provided in Appendix A. The MOGHO algorithm optimizes the of VSS-InC based MPPT and PI controller gains k &k for V regulation. The boost converter is controlled by the gating signal delivered by the MPPT to stabilize the V . The battery current control during the charging and discharging phase is aided by the bidirectional buck-boost converter. The coupling capacitor at the DC link connects the DC and AC systems. The 3P-3W converter is attached to the PCC through interfacing inductors and ripple filters. The non-linear load of 16.5 kW and distribution grid of 415 V at 50 Hz are also attached to PCC.

Implemented Research Methodology
The implemented research methodology is depicted in Figure 2. The V , I ,V , grid voltage magnitude (V ), source voltage ( ) and current ( ) are sensed from the simulated system. The derivative of PV power (P ) from the reference value (P * ) is utilized for the generation of first integral square error (ISE-1) as the first objective function (Obj-1) to be minimized by the MOGHO algorithm as per (1). The optimized VSS-InC MPPT regulates the boost converter. At the DC bus, the comparison of V and (V * ) is formed into second integral square error (ISE-2) as the second objective function (Obj-2) for MOGHO algorithm as per (2). The ISE as error indices is chosen, considering the decaying error with time. The MOGHO algorithm delivers the optimized and k & k for the DC bus control. The accurate current and feed-forward term (W ) are provided to the AKWSOMCC for generating the switching signals of VSC control.

Implemented Research Methodology
The implemented research methodology is depicted in Figure 2. The V PV , I PV ,V dc , grid voltage magnitude (V t ), source voltage (v Sabc ) and current (i Sabc ) are sensed from the simulated system. The derivative of PV power (P PV ) from the reference value (P * PV ) is utilized for the generation of first integral square error (ISE-1) as the first objective function (Obj-1) to be minimized by the MOGHO algorithm as per (1). The optimized VSS-InC MPPT regulates the boost converter. At the DC bus, the comparison of V dc and (V * dc ) is formed into second integral square error (ISE-2) as the second objective function (Obj-2) for MOGHO algorithm as per (2). The ISE as error indices is chosen, considering the decaying error with time. The MOGHO algorithm delivers the optimized δ and k p & k i for the DC bus control. The accurate i Loss current and feed-forward term (W PV ) are provided to the AKWSOMCC for generating the switching signals of VSC control.

Multi-Objective Grass-Hopper Optimization Algorithm and Implementation
The MOGHO is a nature-inspired swarm intelligence (SI) technique that is based o the natural foraging and swarming behaviour of the grasshoppers. The grasshopper lif cycle in two phases, i.e., nymph and adult are considered, where nymph does not hav wings and move in the wind direction, whereas adults have wings and can travel a large distance. The flow chart of the MOGHO is shown in Figure 3.
The grasshopper swarm is mathematically represented as (3), where is the fitnes function, is the social interaction term, is the gravity term, is the wind advectio term. The , , and are the random numbers within [0, 1], which introduces random ness.
The complete equation of the MOGHO is presented as per (4).

Multi-Objective Grass-Hopper Optimization Algorithm and Implementation
The MOGHO is a nature-inspired swarm intelligence (SI) technique that is based on the natural foraging and swarming behaviour of the grasshoppers. The grasshopper life cycle in two phases, i.e., nymph and adult are considered, where nymph does not have wings and move in the wind direction, whereas adults have wings and can travel a larger distance. The flow chart of the MOGHO is shown in Figure 3.
The grasshopper swarm is mathematically represented as (3), where P i is the fitness function, S i is the social interaction term, G i is the gravity term, A i is the wind advection term. The r 1 , r 2 , and r 3 are the random numbers within [0, 1], which introduces randomness.

MOGHO Based VSS-InC MPPT Control
The flow chart for the MOGHO based VSS-InC MPPT is shown in Figure 4. Th InC varies the D with ∆D as perturbation step instead of fixed perturbation in c tional InC algorithm. The MPP is achieved when the variance of output condu equals the negative of output conductance or the slope of the power-voltage (P-V reaches zero as in the conventional InC algorithm as per (10).
The upper and lower bounds for optimization are selected based on experi data presented in [7], which is 1-2.5% of I and likewise the D is also selecte optimization is performed offline to reduce the computational burden of MPP optimized VSS-InC MPPT offers better efficiency, faster MPP tracking and re steady-state oscillations as compared to conventional InC MPPT algorithm. The v and fixed step size MPPT operation are described in Figure 5. With a larger fixed ste the MPP tracking will be very fast but chances of oscillations across MPP will be ve during steady-state, which gets reflected in higher P oscillations. The smaller fix The complete equation of the MOGHO is presented as per (4).
The social interaction of grasshoppers is highly dependent on their social force and the distance between them. The social force (S(r)) is calculated as per (5), where f is the attraction intensity and l is the attractive length scale. The S(r) governs the social attraction and repulsion (or exploitation and exploration) of the grasshoppers. If f and l vary in the interval [0, 2.079] repulsion will occur and if f and l vary in the interval [2.079, 4] attraction will happen in grasshoppers. The grasshopper will be in its comfort zone when f and l remain exactly at 2.079. The d ij is the distance between i th and j th grasshopper andd ij is a unit vector showing the direction of the grasshopper movement.
The G i depending on gravitational constant (g) and unit vectorê g towards the centre of the earth. The A i term is a function of a constant u and unit vectorê w in wind direction. The grasshopper swarm behaviour of reaching and staying in the comfort zone can lead the iterative process towards local minima stagnation. For better chances of convergence, the modified equation of MOGHO is utilized as per (7).
The comfort zone coefficient (c = c 1 = c 2 ) calculated as per (8) will have the same value but a different impact on the fitness function calculation as per (7). The c 1 is similar to the inertial weight in PSO, which reduces the search area as the iterations increases to balance the exploration and exploitation phase. The c 2 term reduces the repulsion, comfort, and repulsion zone (or the search area around global optima) between grasshoppers linearly to increase the chances of convergence and reduce the local minima stagnation. c max and c min are chosen as 1 and 0.000001 for the current problem. At the end of the iterative process, the MOGHO delivers the optimized δ and k p , k i gains. The presented system is aided with the four dedicated controls, i.e., optimized VSS-InC based MPPT control, bi-directional converter control, optimized V dc control and AKWSOMCC based VSC control.

MOGHO Based VSS-InC MPPT Control
The flow chart for the MOGHO based VSS-InC MPPT is shown in Figure 4. The VSS-InC varies the D with ∆D as perturbation step instead of fixed perturbation in conventional InC algorithm. The MPP is achieved when the variance of output conductance equals the negative of output conductance or the slope of the power-voltage (P-V) curve reaches zero as in the conventional InC algorithm as per (10).
The upper and lower bounds for δ optimization are selected based on experimental data presented in [7], which is 1-2.5% of I −1 PV and likewise the D min is also selected. The δ optimization is performed offline to reduce the computational burden of MPPT and optimized VSS-InC MPPT offers better efficiency, faster MPP tracking and reduced steady-state oscillations as compared to conventional InC MPPT algorithm. The variable and fixed step size MPPT operation are described in Figure 5. With a larger fixed step size, the MPP tracking will be very fast but chances of oscillations across MPP will be very high during steady-state, which gets reflected in higher P PV oscillations. The smaller fixed step size increases the efficiency and reduces the steady-state oscillations, but the tracking speed reduces. The variable step size uses a larger step size during irradiation change to reach quickly near the MPP and gradually reduces the step size near MPP to avoid the oscillations.

MOGHO Based Control
The battery current control and V control is shown in Figure 6. The DC bus is regulated by a conventional PI controller, which is tuned with the MOGHO algorithm. The MOGHO algorithm minimizes the Obj-2 and delivers the optimized gains for better stabilization of the V and generation of an accurate current during dynamic conditions [13].

MOGHO Based Control
The battery current control and V control is shown in Figure 6. The DC bus is regulated by a conventional PI controller, which is tuned with the MOGHO algorithm. The MOGHO algorithm minimizes the Obj-2 and delivers the optimized gains for better stabilization of the V and generation of an accurate current during dynamic conditions [13].

MOGHO Based V dc Control
The battery current control and V dc control is shown in Figure 6. The DC bus is regulated by a conventional PI controller, which is tuned with the MOGHO algorithm. The MOGHO algorithm minimizes the Obj-2 and delivers the optimized gains for better stabilization of the V dc and generation of an accurate i Loss current during dynamic conditions [13]. Electronics 2021, 10, x FOR PEER REVIEW 9 of 26 The comparison of and battery current (I ) is fed to another PI controller, by which pulse width modulator (PWM) with logic gate, generates the adequate switching sequence. For the battery overcurrent protection, the I is limited during charging and discharging up to I = 30 A. The I control during the charging and discharging stage is performed by S7 and S8 switches of the bi-directional buck-boost converter.

AKWSOMCC Based VSC Control
The AKWSOMCC based compensator control is implemented by varying the kernel width as a function of the error signal as shown in Figure 7. The proper selection of the kernel width results in the smoothness of function approximation, better convergence rate and accuracy. The larger the kernel width, function smoothness will be more, but reduces sensitivity and may only behave as LMS. With a smaller kernel width, the locality of the training data is preserved, but the convergence rate will slow down. The proposed control with adaptive width brings the Gaussian and non-Gaussian noises of the error signal in the "observation window". The adaptive Gaussian kernel width is a function of error signals, try to cover all Gaussian and non-Gaussian noises under the Gaussian bell curve for filtration. With adaptive kernel width, the whole network matrix or observation window need not be analysed, hence convergence speed increases. The sixth-order MCC acts as a higher-order adaptive filter with sixth-order power optimization. The MSE reduction with higher-order filters produces better steady-state behaviour. The adaptive kernel width ( ) is a function of the error signal and fixed kernel width ( ) as per (11). The pre-decided kernel width ( = 10) will remain fixed during operation [27]. The error signal ( ) is generated as per (12). The comparison of i Loss and battery current (I Bat ) is fed to another PI controller, by which pulse width modulator (PWM) with logic gate, generates the adequate switching sequence. For the battery overcurrent protection, the I Bat is limited during charging and discharging up to I max Bat = 30 A. The I Bat control during the charging and discharging stage is performed by S 7 and S 8 switches of the bi-directional buck-boost converter.

AKWSOMCC Based VSC Control
The AKWSOMCC based compensator control is implemented by varying the kernel width as a function of the error signal as shown in Figure 7. The proper selection of the kernel width results in the smoothness of function approximation, better convergence rate and accuracy. The larger the kernel width, function smoothness will be more, but reduces sensitivity and may only behave as LMS. With a smaller kernel width, the locality of the training data is preserved, but the convergence rate will slow down. The proposed control with adaptive width brings the Gaussian and non-Gaussian noises of the error signal in the "observation window". The adaptive Gaussian kernel width is a function of error signals, try to cover all Gaussian and non-Gaussian noises under the Gaussian bell curve for filtration. With adaptive kernel width, the whole network matrix or observation window need not be analysed, hence convergence speed increases. The sixth-order MCC acts as a higher-order adaptive filter with sixth-order power optimization. The MSE reduction with higher-order filters produces better steady-state behaviour. The adaptive kernel width (σ n ) is a function of the error signal and fixed kernel width (σ 0 ) as per (11). The pre-decided kernel width (σ 0 = 10) will remain fixed during operation [27]. The error signal (e nx ) is generated as per (12).
The is the step size, , , are the unit templates or the in-phase components of source voltage generated as = V ⁄ (where x=a,b,c), where V is the voltage magnitude and will be the phase to ground source voltage. The W is the average weight signal calculated as per (14). The feed-forward term (W ) is introduced in the VSC control for the power balancing is calculated as per (15), where V and I are the PV voltage and current signal.
The W will be the overall weight generated as per (16), which is further utilized for the current reference generation ( * ; where x=a,b,c) as per (17). During specified power mode, the W will replace the W as per (18) and generate the reference currents ( * ) as per (19). The weight signals of each phase W pa (n + 1), W pb (n + 1) and W pc (n + 1) are generated as per (13).
The µ is the step size, µ pa , µ pb , µ pc are the unit templates or the in-phase components of source voltage generated as where V t is the voltage magnitude and v Sx will be the phase to ground source voltage. The W avg is the average weight signal calculated as per (14). The feed-forward term (W PV ) is introduced in the VSC control for the power balancing is calculated as per (15), where V PV and I PV are the PV voltage and current signal.
The W sp will be the overall weight generated as per (16), which is further utilized for the current reference generation (i * sx ; where x = a,b,c) as per (17). During specified power mode, the W sp f ix will replace the W sp as per (18) and generate the reference currents i * sx f ix as per (19).

Results and Discussion
The system is scrutinized under steady-state, irradiation variation mode, specified power mode, unbalanced grid voltage, abnormal grid voltage, unbalanced load, islanding, and re-synchronization mode, with optimized Duty cycle and V dc .

MOGHO Based VSS-InC MPPT Analysis
The V PV performance with MOGHO, PSO, and GA optimized δ for VSS-InC MPPT is compared to the non-optimized or initially selected fixed δ VSS-InC MPPT presented in Figure 8a. The non-optimized fixed δ is selected as 0.0001. During the steady-state at irradiation level of 1000 W/m 2 the V PV variations are minimum with MOGHO tuned MPPT, whereas with fixed δ VSS-InC variations are significant and continuous. The PSO and GA optimized MPPT also reduce the V PV steady-state variations and performed comparatively better than non-optimized MPPT. During irradiation change from 1000 W/m 2 to 600 W/m 2 and vice versa, the MOGHO tuned MPPT manages to track MPP more quickly and accurately to settle down to the reference V PV .

Results and Discussion
The system is scrutinized under steady-state, irradiation variation mode, specified power mode, unbalanced grid voltage, abnormal grid voltage, unbalanced load, islanding, and re-synchronization mode, with optimized Duty cycle and V .

MOGHO Based VSS-InC MPPT Analysis
The V performance with MOGHO, PSO, and GA optimized for VSS-InC MPPT is compared to the non-optimized or initially selected fixed VSS-InC MPPT presented in Figure 8a. The non-optimized fixed is selected as 0.0001. During the steady-state at irradiation level of 1000 W/m 2 the V variations are minimum with MOGHO tuned MPPT, whereas with fixed VSS-InC variations are significant and continuous. The PSO and GA optimized MPPT also reduce the V steady-state variations and performed comparatively better than non-optimized MPPT. During irradiation change from 1000 W/m 2 to 600 W/m 2 and vice versa, the MOGHO tuned MPPT manages to track MPP more quickly and accurately to settle down to the reference V .   The I PV comparison with MOGHO, PSO, GA optimized VSS-InC MPPT is compared with non-optimized VSS-InC MPPT in Figure 8b. During steady-state, MOGHO tuned MPPT maintains the I PV close to its reference value, while performing better than PSO and GA tuned MPPT. Even during irradiation variations, the I PV variations with MOGHO tuned MPPT remain minimum, which dies out more quickly due to faster and accurate tracking of MPP.
The P PV also performs better with MOGHO tuned VSS-InC MPPT in comparison with PSO and GA techniques as shown in Figure 8c. The P PV with adequate step size closely tracks P PV and deliver more average power, hence improving the MPPT efficiency. During irradiation variations, P PV variations are significantly reduced due to the accurate MPP tracking with MOGHO tuned VSS-InC MPPT.
The duty cycle (D) generated by the MOGHO tuned VSS-InC MPPT kept around the desired value, i.e., D = 0.137 as shown in Figure 8d. The D with MOGHO algorithm accurately tracks the MPP and takes very small steps to avoid the chances of oscillations around MPP. During dynamic state, the steps size increases to operate on a new load line and quickly track the shifted MPP for maximum power extraction. The D generated by MOGHO tuned MPPT during changing irradiations reduces the power oscillations to a greater extent and enhances the boost converter efficiency. Table 1 shows the comparative analysis of V PV , I PV settling time and boost converter efficiency with MOGHO, GA, PSO tuned MPPT and non-optimized VSS-InC MPPT with fixed δ. The VSC efficiency with MOGHO, GA, PSO tuned DC bus and non-optimized DC bus along with total dual-stage PV losses are also presented in Table 1.

MOGHO Algorithm Based DC Bus Analysis
The V dc is regulated by the MOGHO tuned conventional PI controller. The V dc should be maintained close to its reference value (V * dc ) to ensure the stability of the system throughout the steady and dynamic states. The MOGHO regulated V dc remains close to V * dc with minimum variations as compared to PSO regulated, GA regulated and non-optimized DC bus (with initially selected gains) during diverse dynamic conditions, i.e., initial transients, irradiation variation, unbalanced load, specified power mode and abnormal grid voltage as shown in Figure 9a. The optimized DC bus generates an accurate loss component of current which is further utilized in battery current control (I Bat ) as shown in Figure 6. The i Loss current is shown in Figure 9b, where the MOGHO based i Loss current have minimum variation as compared to PSO and GA. The MOGHO based accurate i Loss current reduces the I Bat ripples and improves the overall performance of bi-directional converter control. Table 2 presents the optimized design variables, i.e., k p , k i and δ with various optimization techniques and their initially selected or non-optimized values. The initially chosen or non-optimized design variables also serves as the initial condition for the implemented optimization techniques. Table 3 shows the comparative performance of V dc during diverse dynamic conditions with MOGHO, PSO and GA optimization techniques.

Steady-State Analysis
The steady-state is simulated by maintaining the solar irradiation level at 1000 W/m 2 . The total harmonics distortion (THD) of phase 'a' of grid voltage (v Sa ), grid current (i Sa ), load voltage (v La ), and load current (i La ) are satisfactory as per IEEE519 standards are shown in Figure 10a-d.

Steady-State Analysis
The steady-state is simulated by maintaining the solar irradiation level at 1000 W/m 2 . The total harmonics distortion (THD) of phase 'a' of grid voltage ( ), grid current ( ), load voltage ( ), and load current ( ) are satisfactory as per IEEE519 standards are shown in Figure 10a-d. Table 4. presents the comparative THD analysis of and during steady-state with MOGHO, PSO, GA and without optimization technique. With the MOGHO algorithm the proposed system delivers comparatively minimum THD.

Irradiation Variation Analysis
The insolation level is reduced from 1000 W/m 2 to 600 W/m 2 at 0.15 s of simulation time. The grid voltages and currents ( and ) remain in phase opposition, but the magnitude reduces with insolation level such as . The follows its reference currents ( * ) during the induced dynamic condition. The delivered to the grid reduces, as less power is being generated by the PV system while following its reference * . The power PV power ( ) and battery collectively are sufficient enough to satisfy the load requirements. The phase 'a' of load current ( ) and compensator current ( ) remains in phase with each other. The must be satisfied by the VSC itself leaving the sinusoidal nature of the source intact. The follows the * = 0, which is the reference

Irradiation Variation Analysis
The insolation level is reduced from 1000 W/m 2 to 600 W/m 2 at 0.15 s of simulation time. The grid voltages and currents (v Sabc and i Sabc ) remain in phase opposition, but the i Sabc magnitude reduces with insolation level such as P g . The i Sabc follows its reference currents i * Sabc during the induced dynamic condition. The P g delivered to the grid reduces, as less power is being generated by the PV system while following its reference P * g . The power PV power (P PV ) and battery collectively are sufficient enough to satisfy the load requirements. The phase 'a' of load current (i La ) and compensator current (i Ca ) remains in phase with each other. The Q g must be satisfied by the VSC itself leaving the sinusoidal nature of the source intact. The Q g follows the Q * g = 0, which is the reference reactive power exchange from the grid as shown in Figure 11a. The PV voltage (V PV ), current (I PV ) and power (P PV ) are also found following the reference PV voltage (V * PV ), current (I * PV ) and power (P * PV ). The V PV settles down to 600 V after some variations with a boost converter, while I PV and P PV change accordingly with a reduced irradiation level. The battery voltage (V Bat ) and current (I Bat ) depict the battery's charging and discharging rate. The I BAT follows its reference signal (I * Bat ). The V dc is maintained at the desired level of 700 V and keep matching with its reference DC link voltage (V * dc ) is also shown in Figure 11b.
(I * ) and power (P * ). The settles down to 600 V after some variations with a boost converter, while and change accordingly with a reduced irradiation level. The battery voltage ( ) and current ( ) depict the battery's charging and discharging rate. The follows its reference signal (I * ). The is maintained at the desired level of 700 V and keep matching with its reference DC link voltage (V * ) is also shown in Figure  11b.

Load Unbalancing Analysis
The load unbalancing is simulated from 0.35 s to 0.45 s by removing the phase 'a' of the load. The maintains the UPF with while remaining 180° out of phase. The exactly follows * leading towards a stable system. During load unbalancing, only two phases of the load are connected to the system. Hence the load requirement reduces to 66% of the actual load requirement. The delivered to the grid increases during unbalanced load. The exchange with grid remains zero in load unbalancing mode, as shown in Figure 12a. The , and remain as it is with a fixed insolation level. *

Load Unbalancing Analysis
The load unbalancing is simulated from 0.35 s to 0.45 s by removing the phase 'a' of the load. The i Sabc maintains the UPF with v Sabc while remaining 180 • out of phase. The i Sabc exactly follows i * Sabc leading towards a stable system. During load unbalancing, only two phases of the load are connected to the system. Hence the load requirement reduces to 66% of the actual load requirement. The P g delivered to the grid increases during unbalanced load. The Q g exchange with grid remains zero in load unbalancing mode, as shown in Figure 12a. The V PV , I PV and P PV remain as it is with a fixed insolation level. The I Bat becomes more negative and follows I * Bat as more power is available for battery charging. The V Bat also increases as the rate of battery charging increases. The V dc is kept at the required level of 700 V, such as V * dc , as depicted in Figure 12b.

Specific Power Mode Analysis
The fixed or specified power mode is simulated from 0.55 s to 0.65 s of simulation time, where the system has to supply the pre-decided amount of power to the grid for

Specific Power Mode Analysis
The fixed or specified power mode is simulated from 0.55 s to 0.65 s of simulation time, where the system has to supply the pre-decided amount of power to the grid for handling the peak demand. The P g varies along with i Sabc as per the pre-decided grid demand. The i Sabc keep following its reference signals i * Sabc while maintaining the UPF in phase opposition to v Sabc . The i La and i Ca maintain the in-phase relation. The Q g is also held around zero while following its reference Q * g as shown in Figure 13a. The V PV , I PV and P PV shows no variation as the irradiation level is kept at 1000 W/m 2 . The I Bat and V Bat change according to load requirement or charging/discharging phase. The PV and battery are collectively managing the peak demand occurring on the grid. The V dc is also maintained at 700 V, as V * dc , as shown in Figure 13b. handling the peak demand. The varies along with as per the pre-decided grid demand. The keep following its reference signals * while maintaining the UPF in phase opposition to . The and maintain the in-phase relation. The is also held around zero while following its reference * as shown in Figure 13a. The , and shows no variation as the irradiation level is kept at 1000 W/m 2 . The and change according to load requirement or charging/discharging phase. The PV and battery are collectively managing the peak demand occurring on the grid. The is also maintained at 700 V, as * , as shown in Figure 13b.

Abnormal Grid Voltage Analysis
The grid voltage abnormality is simulated by changing the magnitude of to 0.8 p.u. and 1.

Abnormal Grid Voltage Analysis
The grid voltage abnormality is simulated by changing the magnitude of v Sabc to 0.8 p.u. and 1.2 p.u. from 0.75 s to 0.85 s of simulation time. As the v Sabc changes i Sabc has to change to maintain the balance of power. The i Ca and i La remain in phase with each other. The P g changes abruptly with the abnormal grid voltage but settle down quickly. The Q g exchange with the grid is kept zero as Q * g as shown in Figure 14a. The V PV , I PV and P PV remain unchanged other than minor transients. The I Bat changes from charging to discharging mode as the grid voltage goes through sag and swell. The V dc is maintained at the desired level of 700 V as shown in Figure 14b.

Internal Control Signals Analysis
The internal control signals involving , * , W , W , and W are present Figure 15. The remains sinusoidal with unity magnitude throughout the opera The * magnitude changes under the induced diverse dynamic conditions. The W duces during load unbalancing, as the weight signal of phase 'a' diminishes to zero W also varies during abnormal grid voltage as weight signals of each phase chang

Internal Control Signals Analysis
The internal control signals involving µ pa , i * Sa , W avg , W PV , and W sp are presented in Figure 15. The µ pa remains sinusoidal with unity magnitude throughout the operation. The i * Sa magnitude changes under the induced diverse dynamic conditions. The W avg reduces during load unbalancing, as the weight signal of phase 'a' diminishes to zero. The W avg also varies during abnormal grid voltage as weight signals of each phase changes as per voltage sag and swell. The W PV tries to balance the active power flow through the DC link. The W sp varies as a function of W avg and W PV . Figure 16 shows the W avg variations with proposed AKWSOMCC, MCC, LMF and LMS control during diverse dynamic conditions. The AKWSOMCC control delivers W avg with the minimum steady-state and dynamic state error. The convergence speed of the proposed control is also comparable with LMS and LMF controls. As the W avg is having reduced ripples with the proposed control, which offers better VSC performance and reduced overall losses in the system. ditions. The AKWSOMCC control delivers W with the minimum steady-state and dynamic state error. The convergence speed of the proposed control is also comparable with LMS and LMF controls. As the W is having reduced ripples with the proposed control, which offers better VSC performance and reduced overall losses in the system.

Future Research Direction
The authors wish to extend the presented work by utilizing the optimization techniques for reducing the switching losses of the VSC and boost converter and validating it through hardware implementation.  ditions. The AKWSOMCC control delivers W with the minimum steady-state and dynamic state error. The convergence speed of the proposed control is also comparable with LMS and LMF controls. As the W is having reduced ripples with the proposed control, which offers better VSC performance and reduced overall losses in the system.

Future Research Direction
The authors wish to extend the presented work by utilizing the optimization techniques for reducing the switching losses of the VSC and boost converter and validating it through hardware implementation.

Future Research Direction
The authors wish to extend the presented work by utilizing the optimization techniques for reducing the switching losses of the VSC and boost converter and validating it through hardware implementation.

Conclusions
In the 3P-3W grid-tied dual-stage PV-battery storage system, a nature-inspired metaheuristic technique named MOGHO algorithm has been implemented for optimization of modulation factor for VSS-InC MPPT algorithm and PI controller gains for DC bus regulation. The VSS-InC with the optimized modulation factor offers faster and accurate maximum power point (MPP) tracking with reduced oscillations across MPP. With optimization implemented, the settling time of PV voltage and current have been greatly reduced from 13.2 ms to 2.1 ms just after sudden irradiation change along with oscillations. The boost converter efficiency has also swollen from 98.82% to 99.82% with optimized VSS-InC MPPT. The MOGHO tuned DC bus optimization offers enhanced system's stability during dynamic conditions and produces an accurate loss component of current. The accurate loss component of current further improves the AKWSOMCC based VSC control capability by extracting a precise fundamental component of current. The V dc rise time and settling time have been lessened to 0.8 ms and 3 ms with MOGHO in comparison with PSO, GA and non-optimized DC bus. Moreover, with MOGHO the steady-state and dynamic state error of V dc are obtained as low as 0.061% and 0.2% respectively. Due to the combined effect of optimized MPPT and DC bus, the VSC efficiency is grown from 98.47% to 99.23% and collective conduction losses of the boost converter and VSC have also been reduced from 875 W to 308 W. The parameter optimization is carried out on an offline basis, considering diverse dynamic conditions, to obtain the desired system's performance during online operation without slowing down the convergence speed of the system. The AKWSOMCC based VSC control, with adaptive kernel width, acts as a higher-order filter against both Gaussian and non-Gaussian noises, without affecting the convergence speed. The AKWSOMCC control performs multi-functional operations, i.e., harmonics elimination, load balancing, reactive power compensation and active power balancing during dynamic conditions. The presented system operates satisfactorily as per IEEE 519 standards.