key: cord-0781319-fhyiv07j
authors: Bukar, Abba Lawan; Tan, Chee Wei; Yiew, Lau Kwan; Ayop, Razman; Tan, Wen-Shan
title: A rule-based energy management scheme for long-term optimal capacity planning of grid-independent microgrid optimized by multi-objective grasshopper optimization algorithm
date: 2020-10-01
journal: Energy Convers Manag
DOI: 10.1016/j.enconman.2020.113161
sha: 6b1f1f4824b6286f06ee6172dd814184fdf37b0c
doc_id: 781319
cord_uid: fhyiv07j

Off-grid electrification of remote communities using sustainable energy systems (SESs) is a requisite for realizing sustainable development goals. Nonetheless, the capacity planning of the SESs is challenging as it needs to fulfil the fluctuating demand from a long-term perspective, in addition to the intermittency and unpredictable nature of renewable energy sources (RESs). Owing to the nonlinear and non-convex nature of the capacity planning problem, an efficient technique must be employed to achieve a cost-effective system. Existing techniques are, subject to some constraints on the derivability and continuity of the objective function, prone to premature convergence, computationally demanding, follows rigorous procedures to fine-tune the algorithm parameters in different applications, and often do not offer a fair balance during the exploitation and exploration phase of the optimization process. Furthermore, the literature review indicates that researchers often do not implement and examine the energy management scheme (EMS) of a microgrid while computing for the capacity planning problem of microgrids. This paper proposes a rule-based EMS (REMS) optimized by a nature-inspired grasshopper optimization algorithm (GOA) for long-term capacity planning of a grid-independent microgrid incorporating a wind turbine, a photovoltaic, a battery (BT) bank and a diesel generator ([Formula: see text]). In which, a rule-based algorithm is used to implement an EMS to prioritize the usage of RES and coordinate the power flow of the proposed microgrid components. Subsequently, an attempt is made to explore and confirm the efficiency of the proposed REMS incorporated with GOA. The ultimate goal of the objective function is to minimize the cost of energy (COE) and the deficiency of power supply probability (DPSP). The performance of the REMS is examined via a long-term simulation study to ascertain the REMS resiliency and to ensure the operating limit of the BT storage is not violated. The result of the GOA is compared with particle swarm optimization (PSO) and a cuckoo search algorithm (CSA). The simulation results indicate that the proposed technique’s superiority is confirmed in terms of convergence to the optimal solution. The simulation results confirm that the proposed REMS has contributed to better adoption of a cleaner energy production system, as the scheme significantly reduces fuel consumption, [Formula: see text] emission and COE by 92.4%, 92.3% and 79.8%, respectively as compared to the conventional [Formula: see text]. The comparative evaluation of the algorithms shows that REMS-GOA yields a better result as it offers the least COE (objective function), at $0.3656/kW h, as compared to the REMS-CSA at $0.3662/kW h and REMS-PSO at $0.3674/kW h, for the desired DPSP of 0%. Finally, sensitivity analysis is performed to highlight the effect of uncertainties on the system inputs that may arise in the future.

There has been a serious concern with the conventional energy systems owing to the continuous reduction and depletion of fossil fuels, as well as climate change and global warming issues. In this regard, various set-ups (renewable energy sources (RESs)) have been developed as an innovation to the conventional energy systems (traditional fossilbased systems), which are more viable and environmentally friendly, to substitute the inefficient conventional systems [1, 2] . RESs improve the security of energy supply, minimize greenhouse gas emission (GHG), expands energy supply markets, and reduce reliance on imported energy to remote regions [3] . In particular, solar and wind energy can be transformed via photovoltaic (PV) and wind turbines (WT) technology to provide an affordable and reliable electric power supply especially, for remote areas [4] . More importantly, the deployment of renewable energy (RE) systems can contribute to economic development, job creation, and various aspects of human well-being [5] . Accordingly, the RE-based microgrid Positioned at the proximity of the end-user has been endorsed as a substantial approach to resolving the problems involved in supplying electricity to isolated communities [4] . Such as the high cost of transmission due to difficult territory and long-distance, high power losses, and high GHG emissions.

Despite advantages of the RE-technologies, there are obstacles to their roll-out: the sporadic behavior of RESs (e.g., solar and wind energy) coupled with their reliance on climatological conditions, and often the RE technologies are expensive [7] . The cost of energy (COE) generation for a technology also relies on several factors. For the PV systems, COE is influenced by the environmental condition, the efficiency of the PV cell, availability of solar irradiance, and the system sizes [8] . Areas with high solar irradiance have a higher chance for electricity generation using the PV modules since high irradiance is received by the cells. PV cell efficiency is another factor that affects COE. Improved efficiency implies increased in electricity generation in the same condition leads to low COE. The efficiency of PV cells is influenced by operating temperature and the material used for designing the cells [9, 10] . The operating temperature is an additional factor because an increase in temperature reduces the efficiency of the PV cells and vice-versa. Moreover, the output power of the WT generator is influenced by the cut-in, cut-out, and rated speed specified by the manufacturer. Thus, it is recommended to hybridize various RESs in conjunction with an energy storage system (ESS) to form a RE-based microgrid [11] . The integration of RES and ESS requires a strategic and harmonized energy management scheme (EMS) to avoid unnecessary costs on the microgrid. Furthermore, as indicated earlier, the optimal capacity of the RES and the ESS has noticeable role in the overall cost of the microgrid. Thus, computing for the optimal capacity to achieve a minimum cost and reliably supply is demand preferred.

Several studies have focused on the investigation of the optimal design and capacity planning of the hybrid microgrid optimization problems. The optimization problem of microgrids is non-convex, nonlinear, multi-objective in nature, with discrete/integer variables and nonlinear/linear constraints [6] . The modality of the optimization problem is multimodal with many local optimums and a global optimum solution [10, 11] . This has precluded the application of exact mathematical (deterministic) optimizations to compute for the problem. Consequently, it is important to apply a superior optimization technique that guarantees computational tractability for solving such problems.

At the moment, research in artificial intelligence (AI) has resulted in a wide variety of procedures that can either be applied directly or adapted to solve complex problems in the energy sector (for example, in the energy system planning, design, operation, investment, and dynamic control and stability), which were thought to be unsolvable without several simplified assumptions [12] . Metaheuristic or heuristic optimization is a sub-division of AI that computes for global optimum solutions of non-deterministic polynomial time-hard (NP-hard) problems, which are not cooperative with exact mathematical optimization techniques, by using suitable computational variations of the natural systems [11] .

Before the ground-breaking application of particle swarm optimization (PSO) by Hakimi et al. [13] , deterministic or exact optimization approaches, such as mixed-integer linear programming [14] , dynamic programming [15] , linear programming [16] , mixed-integer nonlinear programming [17] , and decomposition-based optimization [18] , were applied to optimally design various forms of renewable and sustainable energy systems. These techniques are subject to some constraints on the derivability and continuity of the objective function; they often cannot track the dynamic changes and follow rigorous procedures. In this context several efforts have been made to harness the potential of the metaheuristic optimization algorithms to make a positive impact on the improvement of the quality of the solutions either by: (1) putting newly developed metaheuristics to test optimization problems using trial-anderror mechanism to determine a technique that will improve the solution quality as compared to techniques previously found to the efficient in solving a particular problem; or (2) developing optimization techniques that are well improved to a particular class of the problem at hand.

As an example of the first class of the works discussed, the application of newly-emerged metaheuristic algorithms has attracted much attention in the literature as compared to the improvement of the standard metaheuristic for RE-system design and capacity planning purposes. An evolutionary PSO (E-PSO) has been devised by Lorestani et al. [19] , and considering WT/microturbine (MT)/battery (BT)/photovoltaic thermal (PVT)/boiler off-grid microgrid, the performance of the devised E-PSO has been proven to be nearing global optimal solutions over harmony search (HS), differential evolution (DE), genetic algorithm (GA), and a standard PSO algorithm. Zhang et al. [20] have developed a novel algorithm by incorporating HS, simulated annealing (SA) and standard chaotic search (CS) algorithm, and name it hybrid HS-SA-CS. The authors have proven the outperformance of the developed approach in minimizing the life cycle cost (LCC) of a grid-independent microgrid as compared to HS and SA algorithms. Dong et al. [21] have applied an improved ant colony optimization (IACO) algorithm for the sizing of a remote microgrid system. The microgrid composed of PV, an electrolyzer (EL), a fuel cell (FC), WT, BT, and a hydrogen (H 2 ) tank. Mohseni et al. [6] have reported the use of a mothflame optimization algorithm (MFOA) to minimize the net-present cost (NPC) of a PV/BT/WT grid-independent microgrid. The best result obtained by the MFOA has been compared to those of the PSO, GA, and GA-PSO methods. Consequently, the authors concluded that MOFOA yields a better result than its counterparts. Likewise, Kaabeche et al. [22] have emphasized the merit of applying a firefly algorithm (FA) for optimal sizing of an autonomous microgrid. The authors consider electricity cost and the load dissatisfaction rate as indicators for computing the system costs and power supply reliability, respectively. Shang et al. [23] have employed improved PSO (IPSO) for the optimum sizing of a BT bank in a grid-independent microgrid. Zolfaghari et al. [24] have computed for the optimal capacity of a BT storage system that is to be combined into a hybrid autonomous microgrid comprising of PV,WT, BT, FC,MT, and a diesel generator (D ) gen using convex optimization. The authors have highlighted the pre-eminence of the convex optimization over PSO and GA, which are the two most wellknown algorithms used for optimal sizing of microgrids in the literature. In [25] , optimal sitting and sizing of a small-scale microgrid system which composed of WT, PV, a biodiesel generator (BD gen ), EL, FC, a H 2 tank reservoir, and a lead-acid BT have been designed using GA. The performance of the proposed approach has been compared with the exhaustive search algorithm (ESA). In [26] , a cost-effective microgrid comprising PV/biomass/wind was designed using artificial bee colony (ABC) algorithm and its superiority was validated over PSO. Maleki and Pourfayaz [27] demonstrated the excellent performance of the HS algorithm over simulated annealing (SA) algorithm, when applied to the microgrid sizing optimization problem. Sharma et al. [28] proposed the used grey wolf optimizer (GWO) for the optimal capacity planning problem of the microgrid BT storage system and have established its dominance over GA, PSO, and several other metaheuristics. Fetanat and Khorasaninenejad [29] applied ant colony optimization (ACO) to solve the size optimization problem of a PV/wind autonomous microgrid. The authors affirmed the supremacy of the proposed technique over GA and ABC. To minimize the total system cost of a hybrid RE system, Ghorbani provides a summary of the most recent work concerning the application of metaheuristic algorithms while computing for the optimal capacity of the grid-independent microgrids. The aforementioned studies have made valuable contributions to the optimal capacity planning of the RE-based microgrids energy system by improving the optimization protocol in the course of solving the microgrid problems. Finding a better metaheuristic algorithm for microgrid planning problems is crucial because of the high cost of its components. Therefore, a minor improvement in the metaheuristic algorithms will have a significant positive impact on the microgrid. These reasons have motivated researchers to examine maiden algorithms and compare them with the existing algorithm for microgrid capacity planning problems. The no-free-lunch theorem [31] has also emphasized on the need for comparisons, and it states that there cannot exist a metaheuristic algorithm that is explicitly best for computing all types of NP-hard optimization problems.

In view of the foregoing, this paper focuses on expanding the current knowledge of applying metaheuristic algorithms for long-term capacity planning of sustainable energy system by making a maiden attempt to study the potential advantages of applying a multi-objective grasshopper optimization algorithm (MOGOA) to optimize the operation of an EMS developed using a rule-based algorithm based on queuing in the planning phase. The GOA is a recently developed algorithm, and the major factor that distinguishes GOA from other algorithms such as PSO and GA is that, it balances the exploitation and exploration phase during the search process, and has proven the ability to avoid getting trapped in the local optimum [32, 33] . The rule-based algorithm is chosen to implement the EMS because it does not require any future journey profile for it to make a decision, thus making it is suitable for real-time application. To the best of our knowledge, no research work in the area has considered the combined rule-based EMS (REMS) and MOGOA to solve the optimal capacity planning problem of microgrids [33] . Furthermore, the review of studies performed indicates that researchers often do not implement and examine the (EMS) of the microgrid while computing for the capacity planning problem of microgrids. The proposed method is applied to optimally design a gridindependent microgrid that constitutes PV, WT, a BT bank, and D gen meant to meet the energy demand of an isolated community. The core contributions of this study can thus be summarized as follows:

• development of an informative mathematical model for a hybrid PV-WT-D gen -BT grid-independent microgrid to provide electricity for remote residential housing units at minimum cost.

• application of a rule-based algorithm based on queuing theory to develop an EMS to prioritize the utilization of RE sources and coordinate the power flow of the elements that form the microgrid.

• application of a state-of-art nature-inspired optimization algorithm to an optimal capacity planning problem of a grid-independent microgrid, and proof of its pre-eminence over PSO and the cuckoo search algorithm (CSA) in terms of nearing the optimal solutions.

• direct comparison of fuel consumption (FC), CO 2 emission, and COE from the proposed grid-independent microgrid under the control of the REMS to an equivalent conventional D gen used in remote communities.

• comprehensive performance analysis of the microgrid against future uncertainties that may arise from the system inputs.

The rest of the study is structured as follows: Section 2 defines the problem by presenting the topology and power flow of a conceptual grid-independent microgrid test case system. This section presents the mathematical description of the components that constitute the test case system, the proposed EMS deployed for the operation of the microgrid, and the description of the climatological conditions and energy demands of the study location. Section 3 describes the proposed optimal capacity planning method that integrates the EMS program. Section 4 details a case study through numerical simulations to validate the efficiency of applying the devised optimal capacity planning method on the proposed microgrid. Subsequently, the conclusion drawn from the study and area for future research is presented in Section 5. To sum up, a chart illustrating the simulation and connection between the sections of the article toward achieving the study objectives is depicted in Fig. 1 .

This section describes the problem by presenting the formation, components, and the power flow of the proposed grid-independent microgrid that is used as the test-case system to confirm the efficacy of the devised optimal capacity planning procedure. It also outlines the microgrid operational scheme that is used in the process of computing for the optimal capacities of the microgrid components. It then presents the projected climatological condition of the case study and power demand data streams, which were used as input data for the analysis. The structure of the proposed grid-independent microgrid is given in Fig. 2 . It incorporates WT, a BT bank, D gen , PV modules, a residential load, and converters.

In this section, the elements of the proposed grid-independent microgrid that are considered in the optimum capacity planning approach are modeled mathematically. The specifications, technical and economic parameters, and product models of components are also presented.

The output power of the WT generator (POW wt ) is a function of the wind velocity (V ) at time step t and the characteristics of the WT generator such as the cut-in speed

) cut out , and the rated speed V ( ) R . The POW wt at time step tcan be computed as follows [34] :

where N wt denote the number of WT. A EOcycle EO10 WT, with 3 blades rated at 1 kW, 48 V DC type is chosen for the analysis. Other details relating to the technical and economic aspects of the chosen WT can be found in Table 2 .

The output power of the PV generator (POW pv ) is a function of ambient temperature T ( ) amb and solar irradiance. The POW pv at time step tcan be computed as follows [34] :

where P r pv is the rated power of the PV panel W ( ), G represent solar irradiance (W m ) 2 , G ref represent the solar irradiance at reference conditions having a value of 1000(W m ) 2 , T cof denote the PV panel temperature coefficient, usually −3.7 × 10 -3 (1/°C) for poly and monocrystalline silicon [35] . T c denote the cell temperature, T ref denote temperature at standard test condition, normally set as 25°C, and N pv denote the optimum number of PV module determined at each iteration by the proposed capacity planning procedure. According to the PV model proposed by Markvar [36] , the cell temperature T c can be computed as follows:

where T noct denote the nominal operating cell temperature (°C), the value depends on the specification of the PV module specified by the manufacturer. A 36-cell monocrystalline PV module (STP275S-20/ Wem) rated 275W P [37] is considered in this study. The technical and economic parameters of the chosen PV module are given in Table 2 .

Due to the inconsistency of the wind and solar sources, the incorporation of BT to the grid-independent microgrid is inevitable. In this context, the state of charge (SOC) of BT bank is acquired based on Ref. [39] . The BT is in charging mode in case the power generated by the RESs (POW wt and POW pv ) is higher than the demand. The quantity of charge at time t is determined as:

where σ denote the self-discharge rate of the BT and is considered to be 0.007%/hour [40] . η B denote the efficiency of the BT during charging and discharging, it is assumed to be 85% in both cases [41] . The BT is said to be in discharge mode when the power generated by the RES is not adequate to fulfil the demand. The quantity of charge at time t is determined as:

The BT bank can supply the demand as long as the SOC t ( ) is greater than SOC minimum (SOC) _ . Similarly, the BT bank can be charged by surplus power generated until SOC t ( ) reaches SOC maximum ( − SOC). The BT capacity in Ampere-hour (Ah) is evaluated in accordance with the desired autonomy day (AD) and the demand using the following formulae [42] .

where E L denote the daily average demand and V s is the system voltage taken to be 48 V. The number of BTs (nBT ) ss in a series string is determined using Eq. (7) . The specifications of the BT used for the analysis are given in Table 3 .

where Max P _ s is the maximum number of parallel string. Accordingly, the number of BT's in each series string nBT ( ) ss , is computed using Eq.

where V B refers to BT bank voltage. The total number of BT required can be determined using Eq. (9).

TheD gen , act as back-up supply in the microgrid, its start ON when the minimum allowable DOD of the BT bank is reached. It is modeled according to its fuel consumption FC ( ) in (Litre/kW) and can be expressed by the following formulae [42] .

P DG and P r DG denote the generated power (kW) and the rated power 0.084 G , respectively [42] . The cost of fuel FUELC ( )throughout the useful lifetime of the system is given by Eq. (11) [42] .

where S f is the current price of diesel fuel per liter. CPV refers to the cumulative present value which translates all the costs throughout the useful life of the system to the initial moment of the investment. It can be computed using Eq. (12) [6, 43] .

where, r i denotes the real interest rate and n' denotes the project lifespan, typically equal to the life span of the PV panels. The CO 2 emitted by the D gen is determined based on the approach adopted by the Intergovernmental Panel on Climate Change (IPCC) which is outlined as follows [44] .

where NCV denotes the net calorific value of fuel (TJ/tons), EF denotes the emission factor (kg CO 2 /TJ fuel). The constants used for the analysis are obtained from Ref. [44] and are given in Table 4 .

The DC-AC converter tying the electrical loads to the grid-independent microgrid is modeled according to its efficiency (η inv ). It is Table 3 The specifications of the 6FM250D BT used for sizing. 

The sporadic nature of RESs is the reason which makes the EMS process complicated. Hence, it is not possible to rely solely on the RESs for meeting energy demand. In this context, the BT storage system and D gen can be integrated with the RESs to maintain a balance between demand and generation. The system operator, however, has to ensure the full utilization of the RES's and minimize the usage of the BT bank and D gen . This can only be achieved via a proper and well design EMS. The purpose of the EMS is to coordinate and efficiently manage the power flow coming from the microgrid subsystems [45] . The proposed EMS in this study is designed based on a cycle charging strategy and is implemented using a rule-based algorithm. The algorithm is mainly in the form of 'if' and 'then' descriptions [46] . Various scenarios are linked to the "if" statements and "then" statements then execute the operating modes. The various modes then execute the energy flow in the microgrid according to certain pre-defined rules. The main scenarios are defined based on the POW pv , POW wt , SOC, (

, and P l . By using the outputs of the previously developed sub-system models of the microgrid in Section 2.1, the algorithm first computes the energy generated by the RESs based on the measured climatological conditions G VandT ( , ) t a m b , and the initial value of the SOC. Afterward, the rulebased algorithm utilizes the data to decide on the system to switch ON or OFF. Regarding this work, four main operating modes are considered for the operation of the grid-independent microgrid and are described hereafter.

• Operating mode 1: If POW pv and POW wt is sufficient to fulfil energy demand, the remaining power will be channeled to charge the BT bank as long as <

. The algorithmic flowchart for this operation mode is shown in Fig. 3a .

• Operating mode 2: If

, the excess power will be discarded to a dump load.

• Operating mode 3: If

, then the energy stored in the BT will be used to fulfil demand requirement.

• Operating mode 4: If 

A case study was conducted to evaluate the optimal capacity of the conceptualized grid-independent microgrid for a remote community in Yobe State, Nigeria. The region is endowed with abundant RE resources (solar and wind energy) and consists of a cluster of an unelectrified residential housing unit. The houses are built 15 feet away from one another, they consist of a living room/dining, two bedrooms, two restrooms, a kitchen, a foyer, and a balcony. The topographical location of the case study region is identified with a latitude of 12.1871°N and a longitude of 11.7068°E. The climatological data were obtained from the data bank of the Nigerian Meteorological Agency (NIMET), for one year (01/01/2017 to 31/12/2017) [47] . The data indicate that the total daily solar radiation of the study location lies between 2, 618 and 7097kW h m day 2 , and the wind speed at an anemometer height of m 10 ranges between 1.3 and 14.4 m s. The heat map of the solar irradiance, wind speed, and ambient temperature extracted from the annual climatological data of the region is shown in Fig. 5 (a-c) . The accurate projection of the energy demand on the proposed grid-independent microgrid was beyond the scope of this work. The demand data used in the study are obtained based on the energy survey study conducted in remote communities by the power holding company of Nigeria [48] and Fig. 5 (d) depicts the seasonal energy demand. Both the climatological and energy demand data were used as input for the optimization.

The proposed optimal capacity planning of grid-independent microgrid consists of the formulation of the optimization problem which constitutes the objective function, constraints and design variables, and the state-of-art metaheuristic optimization algorithm, which are described in detail in the proceeding sub-sections.

Three main objective functions are of great importance in the design of microgrid, namely, the economic aspect objective (cost), the technical aspect objective (reliability), and the environmental aspect objective (CO 2 reduction). These objectives have conflicts with one another. For example, by increasing the reliability of the microgrid, the cost will exponentially increase. Equally, by decreasing the cost, the reliability will further decrease. As such, it is preferred to compute for sets of non-dominated objective functions called Pareto front. In this work, the objective functions that are to be minimized are COE which represents the economic aspect and the DPSP represents the technical aspect.

The COE ($ kW h) is the measure of the average NPC of the electricity generated, and it is expressed as follows [25] :

where TNPC denotes the total NPC. The TNPC encompassed the investment cost INVESTC ( ), the replacement cost REPLACEC ( ) , fuel cost FUELC ( ), and the O&M cost OPEMAINC ( ). Mathematically, the TNPC can be expressed as follows:

The INVESTC is the initial capital cost of system components; it is given by Eq. (17) .

where C WT , C PV , C DG , C CON , C BT refers to the cost of WT, PV, D gen , converter and BT, respectively. N PV , N BT , N DG , N WT , N CON denote the optimum number of PV modules, BT, D gen and WT, determined at each iteration. The OPEMAINC is determined as follows:

where C WT M , , C BA M , , C DG M , , and C PV M , refers to the cost of maintenance for the WT, BT bank, D gen and the PV panel, respectively. The maintenance costs of the converters are neglected. The REPLACEC of the system components can be determined as follow:

Since the WT and the PV panels have a lifespan of 25 years, theirREPLACEC is not considered for them. Hence only the converters, D gen and BT's are expected to be replaced. The expression for the FUELC has been derived previously in Eq. (11).

The reliability model for the grid-independent microgrid is based on DPSP. DPSP is a statistical index that shows the probability of power supply failure. It could be as a result of low energy generation from RESs or due to technical failure. DPSP value ranges from 0 to 1. DPSP equal to 0 implies that demand will be fulfilled in all conditions, while DPSP > 0 implies that demand will not be fully satisfied. It is computed by the following relationship [49] :

where N h denote the time interval and P generation denote the power generated.

Owing to the operational and physical limits of the target microgrid, various restrictions and limits are considered during the optimization process to evade undesired results. Thus, the following constraints are taken into account:

where N AD min 

For consideration of the DPSP, the constraint given in Eq. (23) is used during the optimization process.

where DPSP Desired is a user-defined and implies the allowable value of the DPSP usually 0-1.

In 

GOA is a swarm-based metaheuristic nature-inspired computational method developed by Saremi et al. in 2018 [33] . The algorithm aped the swarming activities of grasshoppers in the course of searching for food throughout their life-cycle. The mathematical equation developed to simulate the swarming behavior of the grasshoppers is given as follows [33] .

where X i denotes the position for the i th grasshopper in the search space. A i , S i , and G i denote the effect of the wind advection on the i th grasshopper, the social interaction of the i th grasshopper, and the gravity force for the i th grasshopper, respectively. The A i and G i components are expressed as follows: 

where d ij refers to the distance between the i th grasshopper and the j th grasshopper and can be computed as

is a unit vector of the distance from the i th grasshopper to the j th grasshopper. The s function denotes the strength of social forces, it is computed as follows:

where f denotes the intensity of attraction and l denotes the attractive length scale. Fig. 6 illustrates the influence of the s-function on the social interaction of the grasshopper (repulsion and attraction). An interval of [0, 15] is considered for the distance d ( ). It can be noticed that the repulsion force between the grasshoppers takes place at a distance between 0 and 2.079 units. After that, the grasshopper has entered its comfort zone, where there is no repulsion or attraction. The attraction of the grasshopper starts to increase after 2.079 until 4 and then further decreases gradually. A more clear demonstration is given in Fig. 7 , where different social behaviors can be observed in the artificial grasshoppers by altering values of l and f . In this work, 0.5 and 1.5 are chosen as the values of l and f during the optimization process. The conceptual model of the repulsion force, comfort force, and the attraction force between the grasshoppers is given in Fig. 8 .

Eq. (24) can be updated by substituting A i , S , i and G i components to solve the optimization problem. The mathematical model proposed by Mirjalili et al. [33] is thus given as follows:

where ub d and lb d denote the upper and lower bounds respectively, in the d th dimension. Cis a decreasing coefficient and is used to shrink the repulsion zone, comfort zone, and the attraction zone. T d  denotes the target value in the D th dimension, and is also referred to as the best solution found. The G i component (gravity force) is ignored (i.e. considered to be zero) in Eq. (29) and A i (wind force) is always in the direction of the best solution T d  . The adaptive parameter C was used twice to simulate the deceleration of the grasshoppers coming close to a source of food and eventually eating it. The first C has been used to minimize the search region toward the target grasshopper as the iteration increases. The second C is used to minimize the effect of repulsion, and attraction forces between grasshoppers proportional to the number iterations to shrink the attraction, repulsion, and comfort zones. The value of the parameter C is updated using Eq. (30) . The mechanism balances between exploitation and exploration of the GOA and proportionally reduces the comfort zone to the number of iteration. C min and C max denote the minimum and maximum values of the adaptive parameter C, t max refers to the maximum number of iterations, and t refers to the current iteration count.

GOA is was applied and Action 1 to 5 were followed to solve the optimization problem of the grid-independent microgrid. All the steps of the proposed approach are clearly explained in Section 3.4.1 and Fig. 9 . At the commencement of the simulation, like any other population-based algorithm, GOA places some random particles in the search landscape. The particles move in the search landscape in accordance with the governing equations of the algorithm to solve the optimization problem. Recently, GOA has been applied to address various RE-based energy system problems, some of the recent works can be found in ref. [50] [51] [52] [53] . 3.4.1. Application of the GOA for optimal capacity planning problem The computational steps of the devised GOA algorithm are as follows:

Step 1: Load input data. a) Load climatological database (ambient temperature, wind speed, solar irradiance). b) Load demand database. c) Load database of the economic indices and technical specifications of the microgrid components (according to Table 2 ).

Step 2: Initialization d) Set the GOA constants: g) Randomly generate an initial swarm of n grasshopper and apply to the objective function to evaluate the fitness of each candidate and determine the best one (T). h) Update adaptive parameter C using Eq. (30).

Step 3: Normalize distance between grasshoppers, update the position of the current search agent using Eq. (29) . Step 4: Update (T) if there is a better solution.

Step 5: Stopping criterion. If the number of iteration exceeds the maximum iterations then stop. The values of the control parameters used for the different algorithms and the overall flow chart of the proposed computational procedures are given in Table 5 and Fig. 9 , respectively.

This section centers on verifying the validity of the contributions introduced within the context of the developed optimal capacity planning method from different viewpoints and presents a cost-benefit analysis of the proposed microgrid based on the results obtained for the considered case study. The proposed microgrid presented in Fig. 2 (which can be considered to be a typical grid-independent microgrid for supplying electricity for remote residential housing) is employed to test and verify the effectiveness of the devised optimal planning capacity procedure. This is conducted by examining the impact of the contributions incorporated in the method compared to other approaches used in the planning studies of microgrid using numerical simulations. All simulations were performed using MATLAB software on an Intel Core i7-4770 CPU, 3.40 GHz PC.

This section presents the results of the resiliency of the REMS before optimizing it using the proposed capacity planning algorithm (GOA). The simulation process is aimed at examining the resiliency of the REMS over a long period and to ascertain that the BT operating limit is not violated despite the occurrence of seasonal variation. In this context, the maximum allowable DOD of the BT is set at 70% and the initial SOC is set at 80%. The PV array was rated 5.8 kW; the BT bank was rated 45.2kW h each 12 V/250 Ah; the WT was rated 5 kW; and D gen was rated 2 kW is used as an input for the REMS to examine the resiliency of the scheme and the rule-based algorithm used for implementing the scheme. Fig. 10(a-c) depicts the few samples of the electricity mix of the various system components that composed the microgrid and its corresponding BT capacity at time step t for the three seasonal variations (cold, hot and rainy season), whereas Fig. 10d depicts the results of the energy mix for a year. By referring to Fig. 10a , the rise, and fall of the P PV is reasonably smooth; this indicates that the irradiance is not disturbed by an abrupt change in the climatological conditions, such as a large cloud passing or rain. During the morning hours (0 h to 8 h), POW wt and POW pv are virtually zero; as such the BT discharges the energy stored in it to fulfil the energy demand (refer to the positive blue dotted line of Fig. 10a) . As the sun starts to rise (8 h-9 h), PV and partially WT starts to generate electricity to fulfil the demand, and BT starts to charge (refer to the negative blue dotted lines of Fig. 10a) . At the start of the sunset (15 h), BT, with the combined effort of WT and PV, starts to fulfil the demand. From there, BT and partially WT continue to supply the demand. Similarly, at 30 h, the BT has reached its minimum allowable limit to discharge (refer to Fig. 10a) . Thus, the D gen instantaneously switches ON to supply demand and starts to charge the BT. The D gen is switched OFF immediately when the sun starts to rise. Fig. 10d illustrates the annual electricity mix under the operation of the proposed REMS. The aforementioned facts that have been demonstrated for 50 h during the cold, hot and rainy seasons can be noticed from the annual energy mix. Despite the continuous deviation in all the quantities, the REMS can able to manage and coordinate the power flow of the microgrid subsystems.

For validation of the proposed GOA, four popular standard benchmark functions comprising unimodal and multimodal functions (Schwefel, Goldstein, Michalewicz, and Sphere functions) were made used. The mathematical formulation of the functions, and their dimension and range are described in detail in Table 6 . The optimum solution found by GOA for test functions were compared with the ones found by PSO and CSA. PSO and CSA were applied to assess the performance of the GOA because of their vast usage in the field of equipment capacity planning of microgrids. The optimum values of the test function were run 30 times by the GOA, PSO, and CSA. Table 6 provides the results of the validation.

Regarding the solutions given in Table 7 , the quality of solutions found by the GOA yields better results in terms of nearing the optimal solution. For the f 1 function, the best solutions found by the GOA, PSO, and CSA are similar, but with slight differences in the mean. For f f , 2 3 , and f 4 , the GOA yields better results compare to its counterparts. Figs. 11 and 12 present the shape of the test functions and their corresponding search history diagrams (convergence process) during the optimization process. It can be observed, that the multimodal test functions (refer to Fig. 11a, b, and 12a ) has many local optimums, which make them worthy for benchmarking of a metaheuristic algorithm. The unimodal function (refer to Fig. 12b ) has one global optimum with no local solutions. As per the search history diagram presented in Figs. 11 and 12 , the grasshoppers are inclined to the promising regions of the search landscape and eventually gather around the global optima (refer to the red dot in Figs. 11 and 12) . The pattern can be noticed in both multimodal and unimodal test functions. This result has demonstrated the capability of the GOA in balancing the exploitation and exploration to drive the grasshoppers near the global optima.

The proposed optimization algorithm is aimed at determining the optimal capacity of the grid-independent microgrid meant to supply the demand of remote residential housing unit, at DPSP Desired . In the literature, the DPSP value is interpreted to be valued from 0% to 1%. DPSP = 0% implies that the microgrid is autonomous full and has the provision of meeting an all-day demand without interruption.DPSP Desired > 0 implies that the microgrid is autonomous basic, i.e. the microgrid cannot fulfil the demand for 24 h a day. The analysis in this section was performed for a DPSP value of 0%. To ensure the reliability of the proposed REMS-GOA in computing the optimal capacity of the grid-independent microgrid required to supply the demand, REMS-CSA and REMS-PSO are programmed and applied to the optimization problem. Accordingly, the values of the different controlling parameters used for the algorithms are given in Table 5 .

The ideal convergence characteristic of the applied algorithms over 30 runs is depicted in Fig. 13 . The ordinate of the convergence curve represents the COE estimate and offers a platform for the comparison of the applied metaheuristic algorithms, whereas the abscissa represents the iteration count.

At the start of the iteration process, the COE values are $0.3951/ kW h, $0.3958/kW h, and $0.3857/kW h for REMS-GOA, REMS-CSA, and REMS-PSO, respectively. Afterward, the values decrease for all the applied algorithms during the iteration process. This means that the algorithms minimize the COE by moving towards the promising regions of the search landscape and eventually stop when the maximum iteration is reached. Consequently, a decrease in the objective function is vital because it leads to having extra information about the optimal values of the design variables. A thorough examination of the convergence curve indicates that REMS-GOA yields a better result as it offers the minimum objective function of $0.3656/kW h, followed by REMS-CSA with $0.3662/kW h and REMS-PSO with $0.3674/kW h. Nevertheless, the REMS-GOA convergence period in the minimum condition is 17 s and better than the REMS-CSA, which is 25 s. The REMS-PSO converges to the minimum condition at 7 s because PSO is prone to premature convergences [11] . Full details of the results obtained for the considered microgrid in the case ofDPSP Desired = 0% are presented in Table 8 . The minimum COE, which is regarded as the best choice, is $0.3656/kW h, as determined by employing the REMS-GOA method. Ranked second is convergence found by REMS-CSA, and Rank third is the one found by REMS-PSO. Thus, it can be concluded that the proposed REMS-GOA offers better performance in all aspects and can successfully be used to find the optimal solution of a complex microgrid design problem compared to its counterparts. The optimal system configuration found consists of a 7150 W capacity of PV array, 4000 W WT, 4.52 kW h BT bank capacity and D gen rated at 5 kW. For the optimal configuration achieved, we show the plot of the annual PV and WT output power in Figs. 14 and 15, respectively. Fig. 16 depicts the annual electricity mix, zoomed version of the electricity mix (at the top left), and the annual SOC of the BT (at the top right). The zoomed version of the SOC for the three different seasons is presented in Fig. 17 . One can notice in Fig. 17b that BT bank is exploited more intensively and is more vulnerable during the hot season, compared to the cloudy and rainy season, this is because of the high demand for energy during the season.

In the previous section (Section 4.3) the two objective functions, COE and DPSP were aggregated into one single objective to achieved DPSP = 0%. The optimization problem then becomes a mono-objective optimization problem. In other words, the COE becomes the main objective; while DPSP changed to a constraint. Hence, the approach has the disadvantage of determining only one single optimal solution (the best system configuration at a DPSP Desired level) that does not represent a trade-off between multiple objective functions. In this light, the current study further exploits the capability of a maiden multi-objective grasshopper optimization algorithm (MOGOA) embedded with Pareto optimal front to optimize the microgrid. MOGOA can simultaneously and independently optimize the objective functions (COE and DPSP), and perhaps the impact of the variation of the DPSP can only be realized by simultaneously considering both objectives. The multi-objective optimization problem (MOOP) is performed by linearly combining two objectives. To illustrate this fact comparatively, the problem is expressed as follows.

The equation has been applied to the two conflicting objectives, = F COE

. w in Eq. (32) denote a weighting factor randomly generated in the range of (0, 1). Its value is set at zero and progressively increased in the step of 0.05 up to 1. φ is a scaling factor chosen as 1000. The two objectives have different units in the MOOP, a penalty factor (PF) is considered to properly balance the objectives.

Two scenarios, case #1 and case #2 each corresponds to 5, and 10 residential housing units are considered to demonstrate the approach. MOGOA is run for 100 iterations and the solutions obtained for the microgrid configurations offer a remarkable and consistent distribution in the non-dominated front (refer to Figs. 18 and 19) . The results obtained not only offer the optimum solution, but a set of optimum tradeoff solutions (non-dominated solutions) in a Pareto front thereby giving the decision-makers several options to choose from. The Pareto-optimal front of the microgrid configuration for case #1 and case #2 is shown in Table 6 Standard benchmark test function. Table 9 .

Using the above technique, multiple sets of Pareto optimal solutions for the grid-independent microgrid system are obtained. For example, if solution #1 of case #1 is chosen by the decision-maker, the microgrid will constitute 6 WTs, an approximate of 20 units of PV modules, 3 days of system autonomy, etc., and the solution corresponds to REF of 90.50%, DPSP of 0.271%, and COE of $0.255/kW h.

Based on the capacity of the microgrid computed using the proposed REMS-GOA for = DPSP 0% Desired (refer to Section 4.3), an overall comparison between the proposed grid-independent microgrid operating under control of the REMS and a conventional D gen system is made; because the residential housing units presently depend solely on conventional D gen system. The comparison is made in terms of CO 2 emission, FC,and COE. Fig. 20a illustrates the ON/OFF operation of D gen under the control of the RfEMS (scenario 1) versus that of the conventional D gen considered as the reference for comparison of power supply for the remote residences (scenario 2). By using the IPCC method discussed in Section 2.1.4 [44] , the two scenarios can be examined in terms of CO 2 avoidance, COE andFC. The results of the comparison of the two scenarios is illustrated in Fig. 21 . The analysis reveals that the proposed microgrid has significantly reduced FC and CO 2 emission produced by the conventional D gen by 92.4% and 92.3% respectively. Additionally, the COE was reduced by 79.8%. Hence, it can be concluded that the proposed microgrid operating under the control of the REMS has a less environmental impact as compared to the conventional D gen power plant taken as a reference. For further examination, it can be observed that the ON/OFF of the D gen is well balanced with the SOC of the BT bank (Fig. 20b) . As seen in Fig. 20a the D gen operates for a longer period during the hot season, this effect can be observed in Fig. 20b , as the SOC of the BT bank, virtually becomes 100% as compared to that in the rainy and cold seasons. For this reason, the D gen is turned on for many hours during the hot season to compensate for the demand when the BT bank is depleted or when it reaches its minimum allowable SOC. 

The amount of energy flow based on running the program simulating the REMS of the microgrid using the optimum capacities of its component obtained by the proposed optimal capacity planning approach is shown in Fig. 22. Fig. 22(a and b) illustrate the total and percentage contribution of power produced by the energy sources in the microgrid during the cold, hot, and rainy seasons, respectively. Such analysis can be applied to other operational horizons (for example, for the whole project's life span) by simply changing and adapting the input parameters to the model. One other point that should not be ignored is the size of the microgrid element computed using the proposed optimal capacity approach for the same operational project lifespan when performing such an analysis. This mechanism will avert the occurrence of discrepancies between the demand and supply within the microgrid. As demonstrated in Fig. 22b , a total of 44% of the energy generated by the microgrid comes from the PV array, and 26% is from the WT. D gen and BT contributed 14% and 16%, respectively. A study was also conducted regarding the excess power produced by the microgrid. In this light, Fig. 23a portrays the plot of the annual excess power produced in each season against the energy demand, and Fig. 23b illustrates the percentage of excess power generated in each season. The results indicate that a total of 981.4 kW excess power would be produced annually, of which approximately 811.3 kW, 12.5 kW, and 171.5 kW would be produced during cold, hot, and rainy seasons respectively. One can note that virtually no amount of excess electricity is produced during the hot season (refer to Fig. 23a and b) , this is because of the high demand for energy during the hot season. Similarly, BT is extremely utilized during the hot season as shown, as illustrated previously in Fig. 20b ; therefore, it is an indication that a smaller amount of excess electricity will be produced during the hot season because BT is not allowed to fully charged during the hot season.

A breakdown of the microgrid subsystem procurement cost for the best combination of the capacities of its sub-system's element computed using the devised method, which implements the REMS and optimizes by GOA for the considered case study, is given in Fig. 24 . The contribution of the procurement cost for the WT accounts for approximately 25.3% of the total procurement cost of the microgrid, the PV array accounts for nearly 32.6%, and the converters account approximately 10%, followed by the BT, which contributes up to 24.2%. Remarkably, unlike other grid-independent microgrids, in which the cost of the BT bank is leading [54] , the BT comes third among the microgrid elements in terms of the initial cost in this project. This is due to the extremely complementary time behavior of the wind and solar resources in the study region. Moreover, the COE computed for the considered microgrid is found to be $0.36563/kW h assuming the proposed microgrid is to fully supply the residential energy demand. One can notice that the electricity price is much lower compared to the usage of standalone D gen , as the COE is found to be $1.81/kW h. Therefore, the calculated COE in this project is competitive to conventional D gen use by most off-grid communities. This makes the project economically viable for installation in Damaturu, Yobe State, Nigeria, since the required energy demand is currently being supplied by conventional D gen . It is expected that the deployment of the project will improve the reliability of power supply in the case study region, reduce GHG, and effectively contribute to the region's sustainable development goals.

This section presents the result of a sensitivity analysis implemented on the grid-independent microgrid. The sensitivity analysis was performed to examine the performance of the microgrid due to uncertainties on the system inputs that may arise in the future. Regarding the analysis, the influence of climatological conditions (wind speed, solar radiation, and temperature), diesel price, BT bank capacity, and SOC set-point is considered in this study. The parameters were carefully chosen based on the following motives:

• The effect of climatological conditions is taken into account to model a worse day scenario (a windy or cloudy sky condition).

• The price of crude oil has in recent times witnessed a substantial and unprecedented decline due to the COVID-19 pandemic, and it is expected to further decrease soon owing to the advent of the electric vehicle and 5G technologies. The result obtained for the optimal system configuration for theDPSP Desired = 0% is chosen to demonstrate the impact of uncertain changes on PV power output, diesel price, WT power output, and SOC set-point, details of the result is given in Table 8 . The sensitivity adjustment was made on a 10% increment/decrement. In this light, the sensitivity curves obtained due to the aforementioned uncertainties is presented in Fig. 25 . The "Base" point on the abscissa refers to the nominal values of the sensitivity adjustment (without increment or decrement). It can be observed that the diesel price has a tremendous effect on COE. As diesel prices increases, COE also increases. On the contrary, the COE value becomes lower for a higher value of PV and WT power and vice versa. This is because the D gen will operate less when the PV and WT generate more power. Equally, an increase in the SOC setpoint will cause a decrease in COE, but to the disadvantage of the BT bank. Finally, one can observe that at the base point, COE achieved is equal to the one computed by the GOA, i.e. $0.3656/kW h (refer to Table 9 ) is chosen as a reference to demonstrate the effect of changes on BT bank capacity. The sensitivity adjustment for the BT bank capacity was made from −50% to 150%. The sensitivity curves obtained are depicted in Fig. 26 . The "Base" point on the abscissa refers to the initial nominal BT capacity (45.2 kW h) chosen for the design. The right and left ordinates depict the DPSP and COE values for different BT capacities. Based on the curves, it can be observed that an increase in the level of BT capacity will reduce the DPSP and lead to an increment in the COE. Any decrease in the BT capacity will lower the COE, although to the detriment of DPSP. This analysis is important because it helps the designer to select the capacity of a BT bank based on the level of individual requirements.

This study considered a conceptual grid-independent microgrid to confirm the validity of the proposed EMS and optimal capacity planning method. The microgrid is meant to meet the energy demand of the residential housing unit in a remote community situated in Yobe State, Nigeria. The RE sources in the proposed microgrid can be updated depending on the RE potential of other regions. Accordingly, the PV and WT were considered as the main energy sources in the microgrid. The selection of the RE generation system, whose optimum sizes are under question, is context-dependent and requires prior viability studies to be conducted on the potentials for the utilization in the target area of the study. To put it differently, the devised method in this work is not principally targeted on and reserved for residential load only. The method can be practically applied to other energy projects with different load profiles (such as agricultural loads, industrial loads, and commercial loads), by simply modifying the corresponding variables and following the same approach.

Since the results of this work support the idea that implementing a proper EMS program while optimally sizing the elements of microgrid leads to the reduction of the overall system cost according to the reduced sizes of some of the elements. Therefore, REMS has been introduced to prioritize the utilization RE first and coordinate the power flow of the microgrid components whilst attempting to highlight the effectiveness of the GOA in solving for the considered problem. The numerical simulation results obtained have demonstrated the efficacy of the abovementioned advanced, cutting-edge approaches in contributing toward minimizing the cost of the grid-independent microgrid and the computational complexity of the problem-solving. Furthermore, since the optimal capacity planning problems of microgrids are nonlinear and non-convex combinational optimization problems, to make it amenable to other mathematical optimization methods, the performance of the proposed REM-GOA method is evaluated through a comparative study with REM-PSO and REMS-CSA in terms optimality of the solution and its convergence when applied to the optimal equipment capacity planning of a conceptualized microgrid for the considered case region. The present study shown has demonstrated that REM-GOA yields a better result as it offers the minimum objective function of $0.3656/kW h, compared to REMS-CSA and REMS-PSO which yields a COE of $0.3662/kW h and $0.3674/kW h respectively at = DPSP 0%

. The simulation results have also indicated that the proposed can REMS has significantly reduced the CO 2 emission, FC, and COE by 92.3%, 92.4%, and 79.8%, respectively.

To make the proposed microgrid smarter, it is recommended that future work should focus on the following key areas. The first is to develop a more sophisticated EMS considering microgrid pool (the interconnection of two or more microgrids) and incorporate it with the proposed optimal capacity method. The second key area is to focus on providing a platform for the consideration of harmonics originating from power conversion devices and nonlinear loads while designing the microgrid. The third area is to address the uncertainties inherent to various inputs variables (e.g., system losses, the inflation rate, and energy demand growth) of the microgrid network using uncertainty analysis methods, such as robust programming, chance-constrained, and Monte Carlo simulations. Moreover, the present study considered a single storage system. Future work can consider the integration of hydrogen storage, in which an EL can be used to produce H 2 . Modern equipment and other household appliances use dc voltage for their operation. Researchers have also discovered the benefit of a dc microgrid for localized load and the notion is to rewire homes to run on dc. Therefore, it is worthwhile to further explore the economic and technical viability of a dc microgrid. Algorithms like chaotic ant swarm [55] , memetic algorithm [56] , cultural Algorithm [57] , salp swarm optimization [58] , along with other nature-inspired approaches used in various optimization studies, are also likely to enrich the future research in RE-based energy system optimization studies. 

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper. 

A heuristic-based approach for optimizing a small independent solar and wind hybrid power scheme incorporating load forecasting

Harmony search optimization for optimum sizing of hybrid solar schemes based on battery storage unit

(a) The annual excess power produced versus energy demand (b) the seasonal percentage of excess power produced

Economic assessment of a pv/diesel/battery hybrid energy system for a non-electrified remote village in Nigeria

Techno-economic analysis of off-grid hybrid pv-diesel-battery system in Katsina State

Economic analysis of residential gridconnected photovoltaic system with lithium-ion battery storage

A demand response-centred approach to the long-term equipment capacity planning of grid-independent micro-grids optimized by the moth-flame optimization algorithm

Transportation in a 100% renewable energy system

Analysis of solar farm site selection based on TOPSIS approach

Analysis of stakeholder roles and the challenges of solar energy utilization in Iran

Energy and exergy efficiencies enhancement analysis of integrated photovoltaic-based energy systems

Economic viability assessment of sustainable hydrogen production, storage, and utilization technologies integrated into on-and off-grid microgrids: a performance comparison of different meta-heuristics

Artificial intelligence in power system optimization

Unit sizing of a stand-alone hybrid power system using particle swarm optimization. (PSO)

Optimal technology selection and operation of commercial-building microgrids

Optimal hydrogen storage sizing for wind power plants in day ahead electricity market

Distributed generation and renewable planning with a linear programming model

Fuel consumption minimization of a microgrid

Economic dispatch and optimal sizing of battery energy storage systems in utility load-leveling operations

Optimal sizing and techno-economic analysis of energy-and cost-efficient standalone multi-carrier microgrid

Sizing a stand-alone solar-wind-hydrogen energy system using weather forecasting and a hybrid search optimization algorithm

Optimal sizing of a stand-alone hybrid power system based on battery/hydrogen with an improved ant colony optimization

Firefly-inspired algorithm for optimal sizing of the renewable hybrid system considering reliability criteria

An improved particle swarm optimization algorithm applied to battery sizing for stand-alone hybrid power systems

Optimal sizing of battery energy storage systems in off-grid microgrids using convex optimization

Techno-economic optimization and environmental Life Cycle Assessment (LCA) of microgrids located in the US using genetic algorithm

Feasibility study of an islanded microgrid in rural area consisting of PV, wind, biomass and battery energy storage system

Sizing of stand-alone photovoltaic/wind/diesel system with battery and fuel cell storage devices by harmony search algorithm

Grey wolf optimisation for optimal sizing of battery energy storage device to minimise operation cost of microgrid

Size optimization for hybrid photovoltaic-wind energy system using ant colony optimization for continuous domains based integer programming

Optimizing a hybrid wind-PV-battery system using GA-PSO and MOPSO for reducing cost and increasing reliability

A new metaheuristic bat-inspired algorithm. Nature inspired cooperative strategies for optimization

Grasshopper optimisation algorithm: theory and application

Grasshopper optimization algorithm for multi-objective optimization problems

An optimal programming among renewable energy resources and storage devices for responsive load integration in residential applications using hybrid of grey wolf and shark smell algorithms

Optimal sizing of autonomous hybrid photovoltaic/wind/ battery power system with LPSP technology by using evolutionary algorithms

Optimal sizing of hybrid PV/diesel/ battery in ship power system

Central Bank of Nigeria (CBN). current interest rate

Optimal allocation and sizing of PV/ Wind/Split-diesel/Battery hybrid energy system for minimizing life cycle cost, carbon emission and dump energy of remote residential building

Design and optimal sizing of hybrid PV/wind/diesel system with battery storage by using DIRECT search algorithm

Increasing household solar energy penetration through load partitioning based on quality of life: the case study of Nigeria

Grey wolf optimizer for optimal design of hybrid renewable energy system PV-diesel Generator-battery: application to the case of Djanet city of Algeria

Multi-objective design of PV-wind-diesel-hydrogen-battery systems

Power management of a hybrid renewable system for artificial islands: a case study

A review on stand-alone photovoltaic-wind energy system with fuel cell: system optimization and energy management strategy

A rule-based energy management scheme for uninterrupted electric vehicles charging at constant price using photovoltaic-grid system

Switched model predictive control for energy dispatching of a photovoltaic-diesel-battery hybrid power system

A novel framework for optimal photovoltaic size and location in remote areas using a hybrid method: a case study of eastern Iran

A regional hybrid GOA-SVM model based on similar day approach for short-term load forecasting in Assam

Optimal voltage and frequency control of an islanded microgrid using grasshopper optimization algorithm

Placement and sizing of multiple distributed generation and battery swapping stations using grasshopper optimizer algorithm

Grasshopper optimization algorithm based two stage fuzzy multiobjective approach for optimum sizing and placement of distributed generations, shunt capacitors and electric vehicle charging stations

A feasibility study of a stand-alone hybrid solar-wind-battery system for a remote island

A novel hybrid chaotic ant swarm algorithm for heat exchanger networks synthesis

Clever algorithms: nature-inspired programming recipes

Review of recent trends in optimization techniques for solar photovoltaic-wind based hybrid energy systems

Salp Swarm Algorithm: a bio-inspired optimizer for engineering design problems

This work is supported by the Federal Government of Nigeria under the Petroleum Technology Development Fund (PTDF) scholarship number PTDF/ED/PHD/BAL/1200/17. They also acknowledge funding provided by UTMShine under vote Q.J130000.2451.09G32, and the facility provided by Universiti Teknologi Malaysia (UTM). The Nigeria Meteorological Agency, Kano is appreciated for providing raw wind, solar irradiance, and temperature data for the study. Lastly thanks to our colleague, Engr. Sara Ayup, who provided expertise and insight that greatly assisted the research.