Battery Energy Storage Systems in Microgrids: Modeling and Design Criteria

: Off-grid power systems based on photovoltaic and battery energy storage systems are becoming a solution of great interest for rural electriﬁcation. The storage system is one of the most crucial components since inappropriate design can affect reliability and ﬁnal costs. Therefore, it is necessary to adopt reliable models able to realistically reproduce the working condition of the application. In this paper, different models of lithium-ion battery are considered in the design process of a microgrid. Two modeling approaches (analytical and electrical) are developed based on experimental measurements. The derived models have been integrated in a methodology for the robust design of off-grid electric power systems which has been implemented in a MATLAB-based computational tool named Poli.NRG (POLItecnico di Milano—Network Robust desiGn). The procedure has been applied to a real-life case study to compare the different battery energy storage system models and to show how they impact on the microgrid design. several modeling approaches (i.e., empirical and electrical) have been compared and simulated. With respect to the related literature, the hallmark of this methodology is the capability to incorporate the uncertainties inherent to rural electriﬁcation processes: load consumption uncertainties, load evolution scenarios, unpredictable energy sources and proper BESS models—able to appropriately represent State of Charge (SoC) and State of Health (SoH) indicators—have been determined as pivotal features to correctly evaluate the robust design of off-grid systems. Once a case study is selected, the system operations are simulated over the entire plant lifetime. The evaluation of the battery replacement costs, the Loss of Load Probability (LLP), the Net Present Cost (NPC) and the levelized cost of electricity are proposed as performance indices to rank the possible microgrid conﬁgurations. Several lifetime load proﬁles are simulated to ﬁnd the optimal microgrid conﬁguration to fulﬁl the desired LLP with the minimum NPC (i.e., heuristic optimization method). The methodology has been proposed in the form of a computational tool in MATLAB R (cid:13) named Poli.NRG (POLItecnico di Milano—Network Robust desiGn) and applied to size a PV+BESS microgrid system to supply power to a rural village of Tanzania.


Introduction
Rural electrification represents one of the issues to be faced in order to provide energy services to 1.2 billion people living in Developing Countries (DCs) that do not have access to electricity [1]. Off-grid systems have an important role in this issue, and especially electrical power systems based on photovoltaic (PV) and battery energy storage systems (BESS). In the rest of this section a brief literature review of the methodologies and software that already address this issue is presented. Several gaps could be detected; this motivates the alternative design procedure discussed in this paper. Section 2 describes the features of the methodology for the robust design of off-grid electric power systems. Off-grid systems strongly rely on energy storage, consequently BESS models are investigated in Section 3. Finally, Section 4 introduces the study case in which the proposed methodology and BESS models are applied to a rural village of Tanzania. Results of the study case are presented in Section 5, where a detailed comparison of the different battery models evaluated aims at understanding their impact on the power system design.
The design process of off-grid systems is not straightforward since it means matching unpredictable energy sources with unknown or uncertain load demands and, at the end, providing the most favorable conditions in terms of reliability and costs. Within the scientific literature, there are three main methods to design off-grid power systems [2]. Each of them requires data about user's in a tool, namely Poli.NRG (POLItecnico di Milano-Network Robust desiGn). The tool is public available on [20].

•
Data inputs gathering and processing-All the necessary information regarding users' electric needs, fixed and variable equipment costs and weather data are collected and processed to obtain load and sources profiles. Specifically, lifetime load profiles are obtained by means of Poli.NRG (LoadProGen subroutine [21]) that can formulate different realistic daily load profiles starting from field data [22]. Renewable source profiles are formulated according to specific models obtained from weather stations or from databases [23,24].

•
System modeling and simulation-The operations of the specific off-grid power system are simulated over the plant lifetime according to specific components' and dispatch strategies models. In the case of a PV+BESS-based microgrid, the main components that must be modeled are PV generator, BESS and inverter [25,26]. The estimation of the PV energy output for each time-step k can be computed considering solar radiation profile, the size of the plant and the balance of system efficiency η BOS , which embraces all the losses not directly related to the sun energy conversion process. The inverter size is defined according to the power peak occurring within the load profile and considering the inverter efficiency (η Inv ). Storage system is considered crucial for the off-grid system design and it can be modeled with different degree of details; empirical and simplified electric models of BESS are hence included as options in the Poli.NRG tool. Section 3 is dedicated to detail the modeling approaches used for BESS.
For a given load curve LC and for given combinations of PV and BESS sizes, Poli.NRG provides techno-economic performance parameters (i.e., Loss of Load Probability, Net Present Cost and Levelized Cost of Energy) as output. Looking at the simulation process, for each time-step k, the difference between photovoltaic energy production EPV(k) and load consumption LC(k) represents the amount of energy that should be fulfilled by the BESS (charging or discharging): In real-life application, constraints are present which limit the BESS working conditions.
First, BESSs must respect technological constraints like the maximum allowable power (defined by its EPR max : maximum Energy Power Ratio). Consequently, the maximum energy that can be fluxed by the BESS in the given time-step k is computed as: For instance, if the EPR max is 0.5 and the BESS size is 1 kWh, the battery can provide or accept a maximum injection of 2 kW. According to this limit the energy that flow through the BESS is evaluated as: Secondly, BESSs must respect their physical limits. Depending on the BESS modeling approach adopted (detailed in Section 3), these limits can be the maximum/minimum SoC levels or the maximum/minimum voltage levels.
In discharge condition, BESS limitations can result in an amount of load that remains unsatisfied because the system is unable to supply it (assuming PV unavailability). This energy is represented by the Loss of Load (LL) indicator and it is computed as follows: For each simulation the system reliability is considered by computing the Loss of Load Probability LLP, which is the share of the electricity demand not fulfilled by the power system over its lifetime (LT) [27,28]: Then the Net Present Cost (NPC), which is defined as the present value of the sum of discounted costs that a system incurs over its lifetime, is calculated [29]: where the first term of the sum is the initial investment cost that is proportional to the installed capacity BESS size given the specific cost of BESS C BESS , the second term represents the sum of the net cash flows during each year y actualized with the expected discount factor r (O&M(y) are the operation and maintenance costs, while C R (y) accounts for replacement costs of BESS by taking into account the projected BESS cost at the specific year y), the third term RV(LT) represents the residual value of the assets (i.e., BESS) at the end of the expected plant lifetime LT.
The Levelized Cost of Energy (LCoE) is also computed since it is a convenient indicator for comparing the unit costs of different technologies over their life, and it is a reference value for the electricity cost that rural consumers would face [30,31]. Moreover, it has also been employed as objective function in several analyses that deal with renewable-based off-grid power systems (e.g., [32][33][34]): • Output formulation-A heuristic optimization method based on the imperialistic competitive algorithm [35] has been developed to manage the microgrid robust design. This method compares the techno-economic performances (LLP and NPC) of different combination of components' sizes in an iterative way. The optimal solution is the specific combination of components' sizes (PV opt ; BESS opt ) which has the minimum NPC value while fulfilling the desired level of LLP [25].
Once the optimal combination of components' sizes for a specific load profile is defined, the same heuristic procedure is repeated for other profiles within the same scenario. The new optimum points are likely to be different because of the different load profiles, hence the optimum points create a set of solutions instead of a single deterministic solution. The most robust solution (PV rbt ; BESS rbt ) within the set of solutions is computed as the most frequent solution among the obtained optimal solutions. New lifetime load profiles are tested till the following convergence criteria are satisfied: where (σ(PV rbt (n)); σ(BESS rbt (n))) indicates the standard deviations of the robust solution given the new n simulated profiles.

Novel BESS Models Proposed for a Proper Microgrid Design
For application of BESS within electrification approaches in DCs, selecting the best technology in terms of performance and durability is not the top priority: cost and availability on site play a relevant role in the final decision. The battery technology choice is usually defined by the economic capacity of the donor or investor. Thus, lead-acid batteries (in particular VRLA), which have benefited of years of development with inevitably cost reduction and global spread, still represent the most common choice for all the applications (Solar home systems, mid-size PV systems, microgrid, and grid-ties back-up systems) [36,37]. As regards lithium batteries, they are rarely seen in DCs applications despite leveraged kWh costs of Li-ion BESS are comparable to lead-acid batteries (when considering the higher depth of discharge, lifetime and number of cycles) [38]. The initial investment cost is a barrier for both investors and donor agencies. However, due to the current experienced reduction in the overall cost of lithium-ion technology [39], it is presumable to forecast that Li-ion BESS will have a relevant share in the sector of off-grid system in the near future.
In such a contest, in this work, Li-ion is chosen as technology of reference for the BESS modeling phase within the adopted methodology for the robust design of off-grid power systems. Two main tasks are identified for battery models: the estimation of the operating conditions (i.e., SoC estimation), the estimation of the ageing (i.e., State of health (SoH) estimation). Two different BESS modeling approaches are considered: empirical and electrical. Empirical models are based on an abstract vision of electrochemical devices. The battery is described by a few analytical equations which do not consider the electrochemical processes behind, but they are empirically fitted. These models usually focus on the evaluation of the SoC of the battery based on energy balances. Electrical models focus on the electrical properties of a battery. An equivalent electric circuit is used to represent the battery dynamics by reproducing voltage and current characteristics at the external terminals. These models could be very simple, with few circuital elements, or more complex, with each circuital element related to a precise physical phenomenon occurring in the cell. These models are computationally more expensive than empirical models, but they can reach a higher degree of accuracy.
Specifically, two different empirical models (model 1 and model 2) are compared with one electrical model (model 3). In all the cases, BESS models have to fulfil the energy balance ∆E(k) between production (PV) and demand (load) and update the State of Charge (SoC) and State of Health (SoH) indicators. The differences lie in the way in which the SoC and SoH are estimated (i.e., the modeling approach). Eventually, these differences will result in marked deviations in LLP and NPC estimations, impacting on the final design evaluation. The models have been developed within the framework of the collaboration between the Politecnico di Milano (DoE department) and CSEM-PV-Center (Swiss Center for Electronics and Microtechnology). All the measurement was performed at the Energy Storage Research Center (ESReC) located in Nidau (CH) on the Li-ion LNCO cell from Boston Power SWING5300 and then extended at the system level for the purpose of the analyses of the present work [40].

Empirical Models
In the empirical model, BESS is modeled as an ideal storage system. Main features are the SoC limits (SoC min and SoC max ), charge (η C ) and discharge (η D ) efficiencies, minimum SoH (SoH min ) and capacity fade index c f (i.e., the loss of SoH per cycle). Moreover, the energy stored in the battery needs to be updated based on the amount previously stored (E BESS (k − 1)) as follows: On this point, the two proposed empirical models are differentiated as regards of the efficiency values: • Model 1 assumes a fix value of round-trip efficiency η RT = 95% as claimed in the literature and/or manufacturers data for similar studies [41]. Charge and discharge efficiency can be derived considering symmetry in charge/discharge processes so that η C · η D = η RT . • Model 2 assumes a variable value of efficiency linked to the operating rate (In a battery system, the operating rate is usually represented by the C-rate that is the ratio between the actual current and the nominal capacity. For the purpose of this paper, C-rate is approximated by the E-rate=(∆E BESS (k)/∆t)/BESS size ) during the specific time-step k. The adopted function is derived from experimental tests and shown in Figure 1.
K Figure 1. Efficiency of the LNCO Li-ion cell BOSTON POWER SWING 5300 measured at ambient temperature with symmetric charge-discharge cycles at different C-rates (DoD = 100%) [42].
Then, the new state of charge (SoC) is computed based on the previous value: where SoH(k − 1) represents the assumed ageing state of the BESS during times step k (SoH computation is detailed in Section 3.3). In this modeling approach, the maximum energy that can be absorbed or provided by the battery ∆E BESS is limited by the SoC limits: SOC min and SOC max . If SOC(k) > SOC max a part of ∆E BESS in Equation (10) is not exploited to charge the battery. In the same way, if SOC(k) < SOC min a part of ∆E BESS generates a loss of load (Equation (4)).

Electrical Model
Electrical models can reproduce the full battery dynamics using equivalent electric circuits. There is a wide range of circuital models, with a different degree of complexity. The use of very complex electrical models like the one proposed in [43] or [44] within the proposed methodology is not suitable since it can bring to expensive computational effort. For this reason, simplified electrical models represent an appropriate compromise. Therefore, the proposed Model 3 is a simplified electrical model made by two circuital elements: the R-C series circuit of Figure 2. Specifically: • C BULK is a capacitor of big-variable capacitance called incremental, differential or intercalation capacitance [44,45]. The voltage drop at the terminals of this element represents the open circuit voltage (OCV) behavior and hence the SoC estimation [46]. Figure 3a shows the trend of C BULK in fresh cell condition which has been derived from OCV measurements (i.e., discharge curve at very low current: C/100) and adopted in the simulations. This curve is commonly called IC curve, referring to the Incremental/Intercalation Capacity. • R I NT is a resistor that models the cell over-potential in regime conditions. It represents the equivalent resistance of the cell at a characteristic stress frequency that depends on the typical BESS working conditions. The resistance has been measured in the lab at a frequency correspondent to the time-step of the procedure adopted. Poli.NRG is based on load and generation profile defined with a time-step equal to 1 min (∆t = 60s): the simulated BESS injections will have a fundamental frequency of 16 mHz. Consequently, the Resistance R I NT has been mapped at this frequency and different SoCs as shown in Figure 3 (fresh cell condition and quadratic fitting function have been assumed). The current flowing through the electrical cell model has been assumed equal to: where the ratio at the denominator represents the number of equivalent cells inside the BESS, obtained dividing the total capacity of the battery BESS size by the nominal capacity of a single cell E n,cell . The number of cells is used as the scaling factor to shift the calculations at cell level. The open circuit voltage V oc (k) and V(k) can be updated at instant k by accounting for the capacitance value C BULK and R I NT at the previous time-step: SoC(k) indicator is directly derived from Voc(k) as detailed in Figure 3c. Ageing is accounted by C BULK and R I NT parameters: SoH and State of Resistance (SoR) indicators are used to update both parameters during lifetime (SoH and SoR computation is detailed in Section 3.3).
In this modeling approach, the maximum energy that can be absorbed or provided by the battery ∆E BESS is limited by the voltage limits: V min and V max . If V(k) > V max , ∆E BESS in Equation (10) is not exploited to charge the battery. In the same way, if V(k) < V min , ∆E BESS generates a loss of load (Equation (4)). A more realistic behavior of the battery is considered compared to the limit of the empirical models.

Ageing Modeling
Beside the SoC estimation approaches for the considered models, BESS models should also be able to represent degradation during cycling (i.e., capacity fade and power fade effects). Understanding of the chemical and physical processes that lead to battery ageing is critical in optimizing battery usage, lifetime, chemistry and construction [47]. The issue is very challenging due to the large number of possible reactions and the interactions of these, along with the highly limited possibilities of collecting data during operation of the battery. In the literature, different degree of details is provided especially when dealing with capacity fade modeling. In [48,49] a cycling counting model is corrected by using a simplified power function which correct the cycle numbers regarding the depth of discharge (DoD) of each cycle. In [50] a more comprehensive ageing model is proposed which analytically correlate the capacity fade with the main ageing factors: cycling rate, DoD and temperature. In this work, all the previously presented BESS models (both empirical ones and electrical ones) accounts for ageing phenomena but with different degree of details: from a purely measurement-independent approach to a fully detailed model closely linked to the experimental results. Ageing studies are time intensive and usually only focus on one or two chemistries and a few ageing conditions. They are therefore not providing a complete understanding of the entire range of temperatures and state of charges. In [51] a wide review is provided on the Li-ion ageing research of commercially available batteries with different cathode and anode chemistries. In the proposed BESS models, ageing phenomena are taken into account in both empirical and electrical models. First, the cycles number registered during the time-step is computed as follows where cy(k) is the equivalent cycle during time-step k. The cumulate function of cy accounts for the cycles till time-step k. SoH can be then updated with the following equation: when SoH(K) = SOH min the battery is assumed to be exhausted and replaced and C R in Equation (6) is updated. As regards of the capacity factor (c f ) (i.e., the loss of capacity per cycle), the three models are differentiated: • Model 1 assumes constant c f that is based on the maximum number of cycles as claimed in the literature and/or manufacturers data for similar studies [52].
For instance, if cy max = 5000 and SoH min = 80%, c f will results in 0.004%/cycle. • Model 2 and Model 3 assume a variable value of c f that is linked to the operating condition, and in particular to the C-rate (Figure 4a).
In general, degradation is modeled differently among the proposed models: • Model 1 and Model 2: use SoH to update the available energy (i.e., capacity fade) as showed in Equation (11).
• Model 3: R and C parameters are updated each time-step. The IC curve of Figure 3a is scaled proportionally to SoH indicator accounting for the capacity fade: while the resistance function of Figure 3 is multiplied by the State of Resistance (SoR) indicator that account for the resistance increase at the time-step k (i.e., power fade): SoR(k) = SoR(k − 1) + cy(k) · r f (k) (20) where r f is the resistance increasing factor that depends on the operating condition during the specific time-step k. r f is hence computed each time-step by following the expression in Figure 4b.  Table 1 resumes the main assumptions of the three models adopted in the simulations. Table 1. Main characteristics of the BESS models adopted in the proposed methodology.

Model SoC Estimation SoH Estimation
Model 1 -Empirical (FIX) SoC limits

The Case Study
The methodology of Section 2 with the different BESS models detailed in Section 3 have been implemented in the tool Poli.NRG and they have been used to run different sets of simulations which are detailed in the following.
The procedure has been applied to design a microgrid, based on PV generators coupled with BESS lithium-ion technology, to supply power to the Ngarenanyuki Secondary School (Tanzania). This school was targeted by the project Energy4Growing [20] and it is a suitable test bed for the procedure since electrical data related to the school are available thanks to the monitoring system in place.
Thanks to a specific survey, electric needs of the Ngarenanyuki school have been collected. Data about the number and the type of appliances in use, their nominal power, and some qualitative information about users' behavior with respect to each appliance are organized and used as input for the Poli.NRG procedure (LoadProGen subroutine). Different possible realistic lifetime load curves are generated to be used in the simulations according to an evolution scenario that considers the load to increase with a constant rate (1%/year).
As regards of the solar resource, approach validated in [53] and the available dataset published in [23] have been used for the simulations. The method employed can take into accounts the variability of the solar resource throughout the year. For instance, lower values of irradiation are expected in February and November, in correspondence to the rainy seasons.
Technical and economical parameters in the simulations are reported in Table 2. Cost information about PV modules, batteries, and off-grid inverters are the result of a survey among Tanzanian local suppliers, while O&M, other investment costs and modeling parameters (efficiencies mostly) have been estimated based on experience. For each single lifetime profile and scenario, the optimal PV-BESS is defined as the one with the minimum NPC, respective a maximum LLP of 5%.

Simulations, Results and Discussion
In this section, the results obtained by different simulations of the study case with Poli.NRG tool are presented. Specifically, the discussions will be focused on two different layers of analysis: 1. Microgrid simulation with different BESS models. A detailed comparison between empirical (M1 and M2) and electrical (M3) BESS models is presented. As regards of the BESS empirical M1, two different values for the maximum number of cycles (cy max ) are tested: (i) M1-1 uses cy max = 10, 000 to reproduce the highest value claimed in the literature for Li-ion technology and (ii) M1-2 uses cy max = 3000 which represents a more reliable data that belongs to the Li-ion LNCO chemistry [40]. Simulations will demonstrate the high impact of BESS modeling approach on SoC and SoH estimations. 2. Microgrid robust design. The focus is on the microgrid robust design as proposed in the presented procedure, results will show how different BESS models lead to different optimal plant configuration.

Microgrid Simulation with Different BESS Models
The first level of analysis has been focused on the comparison of the different BESS models within the novel procedure for the robust design of off-grid systems. Different plant configurations have been selected and simulated with Poli.NRG. The PV size is fixed to 4 kW while the BESS is considered in three sizes: SMALL (1 kWh), MEDIUM (10 kWh) or LARGE (100 kWh). Figure 5 shows the impact of these hypotheses. A one-month winter scenario (February) is considered in the figure since it is characterized by high variability in the solar resources. The power to be delivered by the BESS is the same for the three simulated plants; however, the operating rate (E-rate) is higher for the SMALL BESS (blue line) and lower for the LARGE BESS (red line for MEDIUM and yellow line for the LARGE one). Higher E-rates will result in higher replacement costs making the hypothesis not convenient from a techno-economic point of view; vice-versa, the biggest the BESS capacity the higher the investment costs. This motivates the procedure proposed and implemented in Poli.NRG that looks at the best plant configuration that satisfy a predefined level of LoL at the minimum NPC.   Figure 6 shows the simulated SoC for the four different BESS models proposed within the three hypotheses of plant configuration. The two empirical M1 (black and blue lines), which represent BESS with fix efficiency and a fix number of maximum cycles cy max (M1-1:10,000 or M1-2:3000), are overlapped: differentiation will emerge due to ageing and hence simulating the plant for several years. The empirical M2 (green line), which is characterized by variable efficiency and lifetime according to the operating conditions, distinguishes from M1 in case of MEDIUM/LARGE plants. It results that M1 overestimates the discharge profile: the assumed constant efficiency is probably too severe, and it does not take into account the specific operating condition. This is a limitation of simplified models which use literature data and datasheet to set model's parameters. Special case in the one of SMALL plants in which M1-1, M1-2 and M2 saturate to the limits because of the high operating rate, then hiding the differences in SoC estimations. In general, this is a limitation of all empirical models which do not account for technological limitations (i.e., voltage limits). Conversely, electrical model M3 (red line) differs from empirical models. Charging phases stops at low level of SoC due to overcome voltage limits. The same happens in discharge phases thus affecting LoL estimation. This fact differentiates the simulation achievable with electrical model if compared to empirical models.  Figure 7 shows SoH trends along the plant lifetime. Sudden changes in the SoH indicator highlight replacements during years. In this case, M1-1 and M1-2 differentiate due to the different cy max . Higher cy max (M1-1) means less replacements. Empirical models that relate SoH estimation only on a maximum number of cycles as per Equation (17) without considering cycling conditions, may end up in wrong SoH estimations. This fact is confirmed in case of SMALL plants (undersized BESS) for M1-1: albeit operating rate is so stressful (high E-rate in Figure 5) it is not impacting on BESS replacements. Consequently, replacement is only based on reached maximum lifetime LT BESS,max . This is not the case of the model M1-2 where the assumption of cy max = 3000 significantly increases replacements. In other words, empirical models totally rely on assumptions about BESS aggregated parameters (i.e.,: efficiency and maximum cycle numbers), which should be carefully defined with respect to the final application. These problems are not present in the case of M2 and M3 because they account for the cycling rate in the ageing estimation. Apart from SMALL plant, these two models provide very similar SoH estimation (red and green line are almost overlapped): capacity fade is in fact derived from the same set of measurements ( Figure 4).
All the above findings, about SoC and SoH estimations, impact on the LoL estimations as showed in Table 3. Differences emerge in the simulation of MEDIUM and LARGE plants. M1-1 (cy max = 10,000) results in a lower LoL if compared to M1-2 (cy max = 3000) and this is due to higher capacity fade of the latter case: a lower capacity will indeed affect the LoL negatively. M3 results in the lowest values in all the three hypotheses. Given the higher accuracy provided by electrical model, it results that empirical models tend to overestimate LoL. This fact can be explained by the totally different approaches used to describe battery behaviors between empirical and electrical models: the first one makes use of energy/power balances to compute SoC, while the second one uses electrical quantities (i.e., voltages/currents) to derive SoC. This clearly affects the charge/discharge paths that in the end result into different LoL indicators.

Microgrid Robust Design
Given the results obtained by simulating different plant configurations, the second level of analysis has been focused on the robust design of the case study. Poli.NRG tool, improved with the new energy storage models, has been used to simulate the different load curves and to find the optimal plant configuration for each of them. Figure 8 and Table 4 show the obtained results in terms of sets of solutions. Each set of solution is related to a specific BESS modeling approach. In case of M1 models, around 90 LCs (each one detailed in a 20 years profile at 1 min time-step) are needed before reaching convergence, while M2 and M3 need less LCs: 42 and 62, respectively. The most frequent configuration within the simulated LCs indicates the most robust solution (PV rbt ; BESS rbt ). f rbt parameter indicates how many times the robust solution resulted to be the optimal solution for a specific model. PV optimal size is almost the same (4.10-4.15 kW) in all cases, being closely linked to the solar resources rather than the BESS models. BESS optimal size is instead influenced by the modeling approach. Results depict that BESS should range around 13.70-16.20 kWh with a maximum variation of 18%. Empirical models of type M1 tend to oversize the plant, especially if wrong assumptions about cycling conditions (see M1-1) are made. On the contrary, empirical model of type M2 and the electrical model M3 provide very similar results (i.e., the same PV size and the 2% of difference in between the BESS sizes). Economics are also influenced by the adopted BESS model. NPC and LCoE are computed as the average value among the f rbt simulated LCs that resulted in having the robust solution as the optimal one. Overall, NPC and LCoE for the robust plant configuration are around 34 ke and 0.47 e/kWh respectively. In general, LLP and NPV estimations of model M3 are surely preferable because capacity fade and power fade are the results of capacitance and internal resistance variations: BESS dynamic and ageing behaviors are linked to physical phenomena. But this comes at the price of a considerably higher computational effort: M3 is 14 times slower that M1 and 4 times slower than M2. Therefore, in planning tools like Poli.NRG, modified empirical models which account for capacity fade by using simplified analytical curves (Figure 1) can represent a suitable compromise to obtain reliable results in a reasonable time (Simulation time based on Intel R Core TM i7-4790 CPU 3.6GHz, RAM 16 GB.).

Conclusions
In this paper, a comprehensive methodology for the robust design of PV-BESS systems for rural electrification in Developing Countries (DCs) has been proposed. The attention has been dedicated to developing suitable BESS models to be included in the procedure: several modeling approaches (i.e., empirical and electrical) have been compared and simulated. With respect to the related literature, the hallmark of this methodology is the capability to incorporate the uncertainties inherent to rural electrification processes: load consumption uncertainties, load evolution scenarios, unpredictable energy sources and proper BESS models-able to appropriately represent State of Charge (SoC) and State of Health (SoH) indicators-have been determined as pivotal features to correctly evaluate the robust design of off-grid systems. Once a case study is selected, the system operations are simulated over the entire plant lifetime. The evaluation of the battery replacement costs, the Loss of Load Probability (LLP), the Net Present Cost (NPC) and the levelized cost of electricity are proposed as performance indices to rank the possible microgrid configurations. Several lifetime load profiles are simulated to find the optimal microgrid configuration to fulfil the desired LLP with the minimum NPC (i.e., heuristic optimization method). The methodology has been proposed in the form of a computational tool in MATLAB R named Poli.NRG (POLItecnico di Milano-Network Robust desiGn) and applied to size a PV+BESS microgrid system to supply power to a rural village of Tanzania.
The results showed that BESS models highly affect SoC and SoH estimations and thus LLP and NPC computations. Simplified empirical models, based on literature/manufacturers data might lead to oversized plant, especially if wrong assumptions about cycling conditions have been made. On the contrary, performance estimation made with electrical model has been shown to be preferable because degradation factors (capacity fade and power fade) are the direct results of model's parameters variations (capacitance and internal resistance variations). However, the higher accuracy comes at the expense of higher simulation time which increases by more than ten times. Therefore, modified empirical models which account for degradation factors (i.e., capacity fade) by including simplified analytical curves, instead of relying on manufacturers data, have been shown to represent a suitable solution to obtain reliable results with a viable computational effort. Overall, this paper has highlighted that attention is required when facing the design of off-grid power systems in DCs. Uncertainties on loads, generators and on the mathematical models adopted to simulate BESS behavior, could highly impact on the techno-economic evaluations.