Aalborg Universitet Intelligence-based battery management and economic analysis of an optimized dual-Vanadium Redox Battery ( VRB ) for a wind-PV hybrid system

This paper proposes an intelligent battery management system (BMS) implementing two large Vanadium Redox Battery (VRB) flow batteries in a master-slave mode to provide grid-level energy storage for a wind-solar hybrid power system. The proposed BMS is formulated to effectively meet a predetermined power dispatch formulated based on forecasted wind and solar data while incorporating features like peak shaving and ramp rate limiting. It is compared to a single battery module operated system to showcase the advantages of the proposed intelligent dual battery module in terms of appreciable reduction in battery size and costs while exhibiting improved lifecycle performance. The battery size is optimized based on heuristic optimization algorithms and modelled in Matlab/Simulink environment. An intelligent fuzzy-based BMS is used to control the dual VRB model to ensure optimized power sharing between batteries. The simulations were carried out and an in-depth economic analysis conducted to analyze the costs and other financial metrics of the hybrid project. Results proved the advantages of the dual battery with the proposed BMS and fortify that the introduction of time-based tariffs and other incentives will further make investments in VRB highly attractive for renewable applications.


Introduction
India is currently well placed on its path to address climate change with the Indian government announcing ambitious renewable projects.India aims to achieve 100 GW solar power and 60 GW of wind power by 2022.Wind energy generation capacity in India amounts to 67.9% of total renewable energy capacity by 2017.However, wind generation faces many issues and extensive investments are made for vast repowering activities for wind farms in India [1].Hybrid Renewable Energy Systems (HRES) integrating renewable power generating units with storage options will aim to increase the utilization of available land and help in achieving the solar targets.Multi-purpose repowering by integrating solar Photovoltaics (PV) with the existing wind turbines is being proposed with smart metering.Energy storage systems are increasingly being viewed as an effective solution in hybrid systems with higher renewable penetration to ensure system flexibility and reliability, thus replacing fossil driven dispatch able plants [2][3][4].Other solutions like Demand Side Management also allow for greater renewable penetration in modern micro grids but are often limited by available capacity and customer comfort in intelligent management of microgrids [5].
Battery Energy Storage Systems (BESS) are found to be the most widely implemented storage solutions for wind-PV systems [6][7][8].Many studies are now proposed on optimization [9] and economic studies [10] integrating batteries.Flow batteries prove to be good investment solutions for grid-level storage requirements with their large capacity and modularity in a detailed study conducted for Wind-PV-Battery systems [11].Advances in electrolyte and cell technology have brought down energy costs considerably.Unlimited use of recyclable electrolyte [12] provides umpteen number of operating cycles with little self-discharge and a long shelf life which makes it more attractive for grid level renewable applications.Vanadium Redox Battery (VRB) is a type of flow battery, which is predicted to capture a significant share of the future energy storage markets [13].VRBs have been implemented in many renewable systems in the world and are listed in the Department of Energy (DOE) database [14].
Energy Management Systems (EMS) are essential for batteries to ensure optimal operation and effective service lifetime.Control and management strategies influence the performance of storage systems as explored in the case of a battery integrated hybrid microgrid in [8,[15][16][17][18].Refs [19,20] presented a scheme to utilize two batteries in dual-topology to dispatch wind power where the two batteries were used interchangeably for charging and discharging to deliver wind power with better forecasting and increased reliability.Ensuring full charge/discharge cycle of batteries has proved to improve the operating lifetime of batteries [21].A control strategy for multiple BESS units with individual battery modules dispatched in a rank-based strategy ensured equal battery lifetimes and better performance [22].Intelligent controllers based on fuzzy logic can also ensure optimized and efficient control of hybrid systems integrating renewables and storage systems [23].The above ideas prove that control and management of storage in an HRES can be effectively planned to enhance system performance and storage life.However, modifying these strategies to be implemented with large scale flow battery systems is the core idea of the study.
In India, with the government aiming to encourage renewable penetration, it is high time to explore the concept of integrating and managing energy storage for renewables based on Indian climatic and market conditions.Modelling and operation of flow batteries is still being explored in current literature.The proposed study explores sizing and operation of flow batteries for Indian market conditions by modelling and integrating the VRB as single and dual battery modules in separate cases.The main objective of the paper is to study and analyze the performance of the VRB under two different topologies.The study methodology includes: (i) forecasting and scheduling of HRES output including peak shaving and ramp rate limiting for grid integration; (ii) battery sizing and energy management strategies and (iii) techno-economic analysis of the battery performance for both cases.Thus, this article aims to provide a complete solution pathway for VRB implementation with the HRES and also evaluates and justifies its economic investment case.
The paper is structured as follows.Section 2 includes modelling of system components, solar and wind forecasting and the scheduling of HRES power.In Section 3, an optimal size of the battery is evaluated using bat optimization algorithm and in Section 4, the battery is operated as per the two different cases.In the first case, the VRB is used as a single battery and in the second case; the battery is implemented as a dual-VRB system with a fuzzy logic controller, which manages power sharing and switching of the dual battery.Section 5 discusses the simulated results and battery performance under the two cases based on battery's State of Charge (SOC) and other economic parameters.All simulations are carried out in Matlab/Simulink.

System Modeling and Power Dispatch Scheduling
A test HRES is considered for the study, which consists of a 200 kW MICON M450-200 wind turbine and a 75 kW, PV panel supplying power to the local grid [11].Intermittencies in the wind and solar power make it highly unsuitable for power scheduling and grid integration, leading to heavy losses in the form of spilling and load shedding.Energy storage needs to be sized and operated to improve the HRES power dispatch, which is scheduled, based on the wind and solar forecast data.A VRB is selected due to its modularity in operation, high-energy capacity and high cycle life.
Figure 1a is the schematic diagram of the proposed hybrid system and Figure 1b shows the single and dual battery models.The wind turbine generates AC power and the solar panel and the VRB are connected to the AC grid through converters.Charge controllers used for the VRB in the simulation are bi-directional buck-boost converters.The energy management system performs the switching and controlling operations of the charge controllers to enable reliable operation of the HRES.
Energies 2018, 11, x FOR PEER REVIEW 3 of 18 battery models.The wind turbine generates AC power and the solar panel and the VRB are connected to the AC grid through converters.Charge controllers used for the VRB in the simulation are bidirectional buck-boost converters.The energy management system performs the switching and controlling operations of the charge controllers to enable reliable operation of the HRES.

Wind Turbine Model
The turbine starts generating power Pw(t) once wind speed v(t) exceeds cut-in speed vci up until wind speed reaches the rated speed vr as in the relation (1).As per Betz law, a power coefficient Cp governs the output power curve based on tip-speed ratio λ and pitch angle β.At rated speeds vr, the turbine generates rated power Prated-wt and continues until the cut-off speeds vco are reached, beyond which the turbine is stalled for protection [11,24].

Solar Model
Let Ppv be the power produced from a solar panel of rated power Ypv.Let GSTC and TSTC be the irradiance and temperature incident on the panel under standard test conditions, assumed to be 1000 W/m 2 at 25 °C respectively.Then the expression in (2) is used to evaluate the power generated from a solar panel at time t when Gc(t) is the irradiance in W/sq.m incident on the panel with an operating temperature Tc(t) °C.α is the temperature co-efficient of panel and fpv, the derating factor, normally accounts for conditions like soiling, shading, snow cover etc.Based on the environmental conditions of the site, α and fpv are assumed to be 0.04% and 94% respectively [11].
The cumulative power generated by the HRES Pgen(t) is found after considering efficiencies of wind and PV systems as shown in (3).

Wind Turbine Model
The turbine starts generating power P w (t) once wind speed v(t) exceeds cut-in speed v ci up until wind speed reaches the rated speed v r as in the relation (1).As per Betz law, a power coefficient C p governs the output power curve based on tip-speed ratio λ and pitch angle β.At rated speeds v r , the turbine generates rated power P rated-wt and continues until the cut-off speeds v co are reached, beyond which the turbine is stalled for protection [11,24].

Solar Model
Let P pv be the power produced from a solar panel of rated power Ypv.Let G STC and T STC be the irradiance and temperature incident on the panel under standard test conditions, assumed to be 1000 W/m 2 at 25 • C respectively.Then the expression in (2) is used to evaluate the power generated from a solar panel at time t when G c (t) is the irradiance in W/sq.m incident on the panel with an operating temperature T c (t) • C. α is the temperature co-efficient of panel and f pv , the derating factor, normally accounts for conditions like soiling, shading, snow cover etc.Based on the environmental conditions of the site, α and f pv are assumed to be 0.04% and 94% respectively [11].
The cumulative power generated by the HRES P gen (t) is found after considering efficiencies of wind and PV systems as shown in (3).

Battery Model
Since the 1980s, flow batteries storing electrolytes in separate tanks and circulating them into the battery have expanded the energy capacities of a battery and made the management of electrolytes easier.VRBs are capable of holding charge for long idle periods with almost zero self-discharge.Vanadium, a shiny soft metal used to improve the ductility of steel and soft iron became noted for its ability to gain and lose electrons and led to the exploration of a vanadium-based battery.The cost of vanadium was a major concern factor for the development of vanadium batteries.However, recycled vanadium from oil sludge and fly ash are being proposed with good results [25,26].Scalability and operating ease have made VRB capable of operating with high capacities and short response times.Hence, VRBs are being considered for varied grid level storage applications including black start, supply reserves, frequency regulation, renewable integration, load levelling, energy arbitrage etc. [14].A VRB is modelled as shown in Figure 2 below [27].

Battery Model
Since the 1980s, flow batteries storing electrolytes in separate tanks and circulating them into the battery have expanded the energy capacities of a battery and made the management of electrolytes easier.VRBs are capable of holding charge for long idle periods with almost zero self-discharge.Vanadium, a shiny soft metal used to improve the ductility of steel and soft iron became noted for its ability to gain and lose electrons and led to the exploration of a vanadium-based battery.The cost of vanadium was a major concern factor for the development of vanadium batteries.However, recycled vanadium from oil sludge and fly ash are being proposed with good results [25,26].Scalability and operating ease have made VRB capable of operating with high capacities and short response times.Hence, VRBs are being considered for varied grid level storage applications including black start, supply reserves, frequency regulation, renewable integration, load levelling, energy arbitrage etc. [14].A VRB is modelled as shown in Figure 2 below [27].The SOC of the battery is calculated as shown in ( 4) where, ( ) denotes the power to be charged/discharged to/from the battery in kW and is the VRB capacity in kWh.η is the efficiency of the VRB.The battery is modelled as a controlled voltage source by evaluating its cell stack voltage ( ) and n as the number of cell stacks in the battery architecture.Equivalent cell voltage is calculated as the product of the number of stack cells and the individual cell voltage ( ) which is set to be 1.5 V found from experimental analysis of VRB systems in [13,14] = + 2 1 − where = . ( Extrinsic losses are fixed losses occurring in the external circuitry and the pumps.Equivalent resistance is obtained using internal and parasitic resistance The current to the electrolyte circulating pumps (Ipump) is expressed in terms of stack current (Istack) and state of charge (SOC) of battery as per the below equation [13] = 1.011The SOC of the battery is calculated as shown in ( 4) In the above equation, P b (t) denotes the power to be charged/discharged to/from the battery in kW and E bess is the VRB capacity in kWh.η is the efficiency of the VRB.The battery is modelled as a controlled voltage source by evaluating its cell stack voltage (V stack ) and n as the number of cell stacks in the battery architecture.Equivalent cell voltage V eq is calculated as the product of the number of stack cells and the individual cell voltage (V cell ) which is set to be 1.5 V found from experimental analysis of VRB systems in [13,14] Extrinsic losses are fixed losses occurring in the external circuitry and the pumps.Equivalent resistance is obtained using internal and parasitic resistance The current to the electrolyte circulating pumps (I pump ) is expressed in terms of stack current (I stack ) and state of charge (SOC) of battery as per the below Equation [13] I pump = 1.011

Forecasting and Scheduling of Power Dispatch
Indian states generating wind power currently lack mature wind/solar forecasting systems [28], which makes managing power intermittencies difficult and may lead to voltage/frequency disruptions and grid failure.Wind forecasting methods can be broadly classified into two areas viz.statistical and intelligent techniques [29][30][31].A simple statistical method is employed here to forecast wind data using previously measured data with an error component e(t) [20].The forecast error is assumed to follow a random normal distribution with mean µ and standard deviation of σ.Forecasted wind speed v f (t) is calculated using actual wind speed v as in (8), where v r is the rated output wind speed of the wind turbine.
Solar data can also be predicted with statistical and intelligent learning approaches based on historical data [32].Here, solar irradiance is predicted on a day-ahead basis using multiple linear regressions on historical data, which aims at minimizing the mean square error and enables dynamic interactions between variables.The regression model is given by the below equation where G f is the response variable to be predicted based on the predictor variable and G c -the historical irradiation data.β f is the unknown parameter matrix and ε the error.Energy generated from renewables needs to be dispatched intelligently to the grid to ensure total evacuation.Many dispatch strategies have been proposed in the literature involving energy storage systems.Q. Li et al., [21] introduced a min-max strategy of power dispatch which was later improvised in [33] for achieving a better battery life.The Indian load dispatch centers follow a 15-min time block for dispatching the generated power.The scheduling is done every 30-min covering two time blocks.This is to ensure that any intermittency in the renewable power generated in this time period will be nullified by the VRB, to enable a steady HRES output.This leads to a better evacuation planning scenario for the grid operators.The determined schedule allows the VRB to save excess power generated at off-peak times and then discharge the saved power back to the grid at peak times.Let Dem(t) be the power to be dispatched from the HRES to be delivered at time t.Table 1 enlists the evaluation formulas for Dem(t) with peak shaving.PG av (i) is the average power generated from renewables (P gen ) for i th day based on forecasting.Next, dispatch power Dem is checked for ramp rate limit violations to exert a maximum power evacuation limit.A ramp rate limit of 15 kW every 10 min is imposed, and any excess energy is directed to be managed by the VRB [11].

Sizing Methodology
The VRB flow battery needs to balance the power difference P diff (t) occurring between the forecasted scheduled power Dem(t) and the actual power generated P gen (t): Optimization is used to find the smallest possible battery size (C vrb ) satisfying all constraints at minimum investment costs (IEEE Std.1184 has been followed for sizing).It is to be understood that the smaller the battery size, the more difficult it is for the battery to meet P diff ; thereby resulting in insufficient storage capacity and associated losses.Hence, the objective function is coined as in (8) to minimize battery cost while satisfying all other constraints listed in ( 16)- (19) [11,34].The addition of loss components in the objective function tends to balance the battery size for optimum capacity.The objective function includes cost coefficients for investments (χ) and maintenance costs (γ) with time factors (CRF) (See Nomenclature) where CRF = idr(1+idr) n (1+idr) n −1 is the capital recovery factor and In the above equations, idr, ifr, ir are the discount rate, inflation rate and interest rate on investment evaluated for n years of VRB battery life, LC denotes the summation of Shedding_Loss and Wind_spilling _loss calculated as below: Spilling losses = ∑ r i=1 P gen (t) − Dem(t) × ∆t × PR ws when Dem < P gen (14) Let the revenue earned by delivering P del (t) be evaluated based on the wind tariff of the region as Revenue Earned = ∑ r i=1 P del (t) × ∆t × PR ws (15) Constraints bounding the search algorithm are: SOC min ≤ SOC(t) ≤ SOC max where where The optimization problem is solved using bat optimization algorithm.Xin-She Yang formulated Bat optimization algorithm in 2010 [35].It is a nature-inspired optimization algorithm, which was derived from swarm intelligence wherein a group of microbats uses echolocation technique to find their prey.Many studies [36,37] have established that Bat algorithm can yield comparatively good search results.Initial conditions of the algorithm are set by defining the frequency (Q min and Q max ), amplitude (A i ) and rate (r i ) of sound pulses emitted by the participating microbat solutions.The best solution is selected by identifying the solution with minimum cost (x*).Velocities and positions of all bats are updated as in Equation ( 7) and search is continued until convergence.

Energy Management of Dual Battery
In this section, the battery thus sized is implemented with two different management strategies.VRB can be operated as modules connected in series and parallel to scale up the voltage and current respectively.In the first case, the VRB is operated as a single model and the second case the VRB is implemented in form of a dual battery model to where the two VRB modules are operated in complementary fashion as shown in Figure 1b.The simulations are run, and results are compared with the dual battery model.

Single Battery Model
The VRB is operated as a single model, thus absorbing energy when P gen > Dem and discharging otherwise as shown in Figure 3a.The battery SOC and the power limits are monitored to avoid any under or overcharging of the battery.At times of SOC violations, the battery is idle and power losses are calculated.

Single Battery Model
The VRB is operated as a single model, thus absorbing energy when Pgen > Dem and discharging otherwise as shown in Figure 3a.The battery SOC and the power limits are monitored to avoid any under or overcharging of the battery.At times of SOC violations, the battery is idle and power losses are calculated.

Dual Battery Model
VRB is implemented as a Dual model with two VRB modules viz.BESS 1 and BESS 2 and the EMS for the dual module is depicted in Figure 3b.Each battery module is capable of operating independently with BESS 1 in master-charging mode and BESS 2 set to be in slave-discharging mode and vice versa.They toggle their roles on reaching their predetermined SOC limits, i.e., when any one gets fully charged it is set to discharge and the other module charges respectively.The status of each battery can be indicated using select lines Sel 1 and Sel 2 , which when high indicates the corresponding battery in master/charging mode.At times when the power difference to be met is greater than the maximum power limits (P b,max ) of single module i.e., P diff (t) > P b,max , then the two batteries optimally share the power requirement based on their SOC.Hence, there arises a need to balance the sharing operation between the battery modules.This is implemented intelligently by a Fuzzy Logic Controller (FLC), which will decide the mutual sharing between modules based on their SOC and evaluate the power to be contributed by each battery module individually.Each fuzzy input is mapped to its member function (fuzzification) depending on its influence on the output.The output of each fuzzy set is inferred using fuzzy rules outlined in the controller.Defuzzification in the final step resolves the fuzzy outputs into a crisp number format [23,38].In this study, the FLC designed is a mamadani type fuzzy system built from four inputs namely, power difference component (dP), SOC of battery BESS 1 (SOC 1 ), SOC of battery BESS 2 (SOC 2 ) and BESS 1 status signal (Sel 1 ).The dP is the power that is to be charged/discharged by the battery bank.The FLC takes information of battery capacity at each interval through its inputs SOC 1 and SOC 2 to ensure the batteries do not overcharge/over discharge.The outputs from the fuzzy system are two parameters b 1 and b 2 which in turn give the selection information of batteries BESS 1 and BESS 2 respectively.
Figure 4a-d show the membership functions of inputs and outputs of FLC and trapezoidal functions are selected to represent them.SOCs are shown to be low (L) or high (H) or medium (M).The dP is shown to have four levels viz, Positive Big (PB), Positive (P), Negative Big (NB) and Negative (N).If dP is positive and small (< P b,max limits) it is positive and if more than P b,max , it is Positive Big (PB).Similarly, if dP is negative and small (< P b,max limits) it is N and if more than P b,max it is Negative Big (NB).The output parameters b 1 and b 2 of the FLC denote the selection status of batteries.The dP is checked for maximum battery power violations i.e., if abs (P diff (t) > P b,max ), then it indicates a condition when both the batteries need to operate together to meet P diff (t).The optimal sharing in this condition is implemented based on the output parameters of the FLC and Lagrange's multipliers.A total of 18 fuzzy rules are outlined based on the energy management strategy in Figure 3b and listed in Table 2 with all conditions and appropriate actions detailed.

Dual Battery Model
VRB is implemented as a Dual model with two VRB modules viz.BESS1 and BESS2 and the EMS for the dual module is depicted in Figure 3b.Each battery module is capable of operating independently with BESS1 in master-charging mode and BESS2 set to be in slave-discharging mode and vice versa.They toggle their roles on reaching their predetermined SOC limits, i.e., when any one gets fully charged it is set to discharge and the other module charges respectively.The status of each battery can be indicated using select lines Sel1 and Sel2, which when high indicates the corresponding battery in master/charging mode.At times when the power difference to be met is greater than the maximum power limits (Pb,max) of single module i.e., Pdiff(t) > Pb,max, then the two batteries optimally share the power requirement based on their SOC.Hence, there arises a need to balance the sharing operation between the battery modules.This is implemented intelligently by a Fuzzy Logic Controller (FLC), which will decide the mutual sharing between modules based on their SOC and evaluate the power to be contributed by each battery module individually.Each fuzzy input is mapped to its member function (fuzzification) depending on its influence on the output.The output of each fuzzy set is inferred using fuzzy rules outlined in the controller.Defuzzification in the final step resolves the fuzzy outputs into a crisp number format [23,38].In this study, the FLC designed is a mamadani type fuzzy system built from four inputs namely, power difference component (dP), SOC of battery BESS1 (SOC1), SOC of battery BESS2 (SOC2) and BESS1 status signal (Sel1).The dP is the power that is to be charged/discharged by the battery bank.The FLC takes information of battery capacity at each interval through its inputs SOC1 and SOC2 to ensure the batteries do not overcharge/over discharge.The outputs from the fuzzy system are two parameters b1 and b2 which in turn give the selection information of batteries BESS1 and BESS2 respectively.
Figure 4a-d show the membership functions of inputs and outputs of FLC and trapezoidal functions are selected to represent them as all inputs representing intervals.SOCs are shown to be low (L) or high (H) or medium (M).The dP is shown to have four levels viz, Positive Big (PB), Positive (P), Negative Big (NB) and Negative (N).If dP is positive and small (<Pb,max limits) it is positive and if more than Pb,max, it is Positive Big (PB).Similarly, if dP is negative and small (<Pb,max limits) it is N and if more than Pb,max it is Negative Big (NB).The output parameters b1 and b2 of the FLC denote the selection status of batteries.The dP is checked for maximum battery power violations i.e., if abs (Pdiff(t) > Pb,max), then it indicates a condition when both the batteries need to operate together to meet Pdiff(t).The optimal sharing in this condition is implemented based on the output parameters of the FLC and Lagrange's multipliers.A total of 18 fuzzy rules are outlined based on the energy management strategy in Figure 3b and listed in Table 2 with all conditions and appropriate actions detailed.

Simulation Results and Discussion
The solar and wind data for the sizing study were recorded from a wind farm in Southern Tamil Nadu over a period of 10 days at 10-min intervals and plotted in Figure 5a.Necessary data regarding wind turbine modelling is given in Table 3 [39].The simulations run for a period of 10 days.The Wind-PV hybrid system generates 23.205 MWh of energy during the simulation period.Based on the assimilated data, forecasting is done, and the power generated from the HRES is formulated as explained in Section 2. The mean absolute percent error for wind forecasting done based on statistical analysis was 9.3% and for solar using, multiple linear regression was 12%, which falls within acceptable standards [4].The power dispatch (Dem) is then scheduled as proposed in Section 2.4. Figure (i) shows the scheduled power dispatch (Dem) plotted against actual power generated by the HRES for the simulation period.Figure 5b (ii) is plotted for only 24 h to illustrate how the dispatch curve (Dem) follows peak shaving and ramp rate limiting while holding constant output for 30 min (two dispatch time blocks).Without any option for storage, huge power mismatches can be clearly seen here.18.54 MWh of energy is delivered to the grid with 3.64 and 4.64 MWh energy lost due to load shedding and wind power spilling respectively.The system suffers a total Loss of Power Supply Probability (LPSP) of 19% calculated as per the following equation

Single VRB Model Strategy
The VRB characteristics considered for the study are listed in Table 4.The size of battery system found using bat algorithm is 4600 Ah [11].Hence, three VRB modules, each of capacity 200 kWh, are operated in parallel as a single battery bank with a capacity of 5000 Ah.The simulation results prove that the VRB enables the HRES to effectively follow the scheduled dispatch curve with a LPSP of 3.06%.The battery SOC chart is plotted as in Figure 6a.The chart clearly shows underutilized battery capacity with many incomplete charging/discharging cycles, which may adversely affect the battery performance and lifetime.Figure 6b represents the power delivery curves for the single battery system against Pgen.

Single VRB Model Strategy
The VRB characteristics considered for the study are listed in Table 4.The size of battery system found using bat algorithm is 4600 Ah [11].Hence, three VRB modules, each of capacity 200 kWh, are operated in parallel as a single battery bank with a capacity of 5000 Ah.The simulation results prove that the VRB enables the HRES to effectively follow the scheduled dispatch curve with a LPSP of 3.06%.The battery SOC chart is plotted as in Figure 6a.The chart clearly shows underutilized battery capacity with many incomplete charging/discharging cycles, which may adversely affect the battery performance and lifetime.Figure 6b represents the power delivery curves for the single battery system against P gen .

Dual VRB BESS Module Strategy
The VRB bank is now implemented in the form of two separate VRB modules: BESS1 and BESS2, each 200 kWh 120 V (thereby two batteries of capacity 1666.66Ah); connected in parallel and operated for a cumulative capacity of 400 kWh (3333.33Ah).At the beginning of the simulation period, Sel1 of BESS1 is held high indicating that BESS1 is to be charging and BESS2 is for discharging.
The SOC charts of the dual battery modules are shown in Figure 7a,b respectively.Figure 8a shows the power delivered by the HRES plotted against the scheduled power.This figure shows that the dispatch is being delivered as scheduled with an LPSP ratio of 4.3%.The power mismatch curve is also plotted as Figure 7b where the negative quadrant shows losses due to spilling and positive for the load shed.Figure 9a shows the power that is charged/discharged from the two batteries.It was observed that the FLC effectively managed the two batteries so that optimal sharing of the load was

Dual VRB BESS Module Strategy
The VRB bank is now implemented in the form of two separate VRB modules: BESS 1 and BESS 2 , each 200 kWh 120 V (thereby two batteries of capacity 1666.66Ah); connected in parallel and operated for a cumulative capacity of 400 kWh (3333.33Ah).At the beginning of the simulation period, Sel 1 of BESS 1 is held high indicating that BESS 1 is to be charging and BESS 2 is for discharging.
The SOC charts of the dual battery modules are shown in Figure 7a,b respectively.Figure 8a shows the power delivered by the HRES plotted against the scheduled power.This figure shows that the dispatch is being delivered as scheduled with an LPSP ratio of 4.3%.The power mismatch curve is also plotted as Figure 7b where the negative quadrant shows losses due to spilling and positive for the load shed.Figure 9a shows the power that is charged/discharged from the two batteries.It was observed that the FLC effectively managed the two batteries so that optimal sharing of the load was At times, when the difference between dispatch and generation is beyond the maximum charge/discharge limit of a single VRB module, both batteries operate together.The voltages and the currents of the individual modules in the dual battery are plotted in Figure 9b,c respectively.achieved.At times, when the difference between dispatch and generation is beyond the maximum charge/discharge limit of a single VRB module, both batteries operate together.The voltages and the currents of the individual modules in the dual battery are plotted in Figure 9b,c respectively.

Inference of Results
By implementing the VRB in the dual-battery topology with effective energy, management of the size of the battery is now reduced to 400 kWh modules (two batteries of 200 kWh in parallel) as compared to the 450 kWh of the Single Battery model with higher power limits.In Figure 7c, the SOCs of all batteries are plotted for a period of 10 h to illustrate the behaviour of systems.It should be noted that single battery is not a sum of the two dual batteries but a much larger system.Comparing the SOC charts from Figure 6a and Figure 7a,b, it can be seen that in dual-model, each battery module undergoes fewer switching cycles and hence can deliver for a longer period.As a single module, VRB underwent nearly 17 charge-discharge cycles during the simulation period but it is clear that the battery capacity is underutilized as it has many incomplete charge/discharge cycles.In addition, while the VRB is used as dual-battery topology, the batteries did not experience more than 10-12 cycles for much-reduced size, implying a reduction of 30% of switching cycles.

5.4.Inference of Results
Through effective energy management of VRB in the dual-battery topology, the size of the battery is now reduced to 400 kWh modules (two batteries of 200 kWh in parallel) as compared to the 450 kWh of the Single Battery model with higher power limits.In Figure 7c, the SOCs of all batteries are plotted for a period of 10 h to illustrate the behaviour of systems.It should be noted that single battery is not a sum of the two dual batteries but a much larger system.Comparing the SOC charts from Figures 6a  and 7a,b, it can be seen that in dual-model, each battery module undergoes fewer switching cycles and hence can deliver for a longer period.As a single module, VRB underwent nearly 17 charge-discharge cycles during the simulation period but it is clear that the battery capacity is underutilized as it has many incomplete charge/discharge cycles.In addition, while the VRB is used as dual-battery, the batteries did not experience more than 10-12 cycles for much-reduced size, implying a reduction of 30% of switching cycles.
Thus, with reduced size, increased life and optimized cycle operation, the dual battery strategy is clearly a more efficient and better choice than the single module case.Figure 7d shows the power exchange between the battery systems and it is clear that the single battery has linear charging due to higher capacity but it needs to cater to all discharge/charge cycles and hence experiences frequent charge reversals.Whereas, in the dual case, both batteries share the charging/discharging operation and can reduce intermittent charge reversals and hence, operate better than the single battery (at instances of 37-40 h).The dual model results in a slightly higher LPSP ≈1% than the single module, which should be acceptable.A further understanding of the performance of the battery systems is explored by evaluating stress factors like the charge (CF), partial cycling values (PCV), time at low SOCs (TLS) and Ah Throughput of the battery; these are listed in Table 5 [40,41].The discharge rates are usually limited by limiting the P b,max for the batteries in the EMS.The results show that operating the battery modules in Single/Dual architecture has little impact on the CF or PCV values, which both prove optimal battery operation.However, the single model case suffered from very low TLS, showing underutilized capacity, which may affect the battery performance as the battery modules operate in limited charged condition characterized by low Ah Throughputs [40].Thus, with reduced size, increased life and optimized cycle operation, the dual battery strategy is clearly a more efficient and better choice than the single module case.Figure 7d shows the power exchange between the battery systems and it is clear that the single battery has linear charging due to higher capacity but it needs to cater to all discharge/charge cycles and hence experiences frequent

Economic Analysis
Flow batteries face major barriers regarding financial concerns for implementation in wind projects.[42,43] are some articles available in the literature, which attempt to explore the cost-effectiveness of battery technologies for applications like load shifting and utility scale integration.A similar approach has been attempted to explore the investment returns of the proposed system.Assuming net investment to be 600 $/kWh, the total investment for the two battery modules amounts to $240,000.An economic evaluation is carried out to understand the feasibility of investment and results tabulated in Table 5. Assume the power generated is sold to the utility at the rate of Rs. 2.75 per kWh, as per tariff fixed by the Tamil Nadu Electricity Regulatory Commission [44].The revenue lost due to wind power spilling and load shedding due to inadequate storage are evaluated as per (13)(14)(15)(16).The total revenue earned by operating the Wind-PV system (without storage) is $865 and the losses incurred is $1711, as maximum losses occur in the form of spilled renewable power which could not be accommodated in the grid.The same calculations are repeated after implementation of VRB in both strategies.Dual VRB HRES delivers 21.78 MWh of energy thus earning a total amount $1000.04.The losses now amount to $78.32.Thus, a total of about $1767 is saved in 10 days.Then the annual savings will amount to nearly $63,636.Thus, the battery investment can be earned back within five years of installation and it almost equals the payback for single case.If interest rates and inflation rates are to be considered (assume ir = 12%, ifr = 7% and γ = 3%) then the discounted payback period increases to about seven years.The total lifetime of the VRB, assuming it delivers 10000 cycles, is approximately 20 years (No. of battery cycles = 13 in 10 days).Thus, the system has a Net Present Value of $313,400 and a positive Benefit to Cost ratio of 1.949, both higher than the single battery case due to lesser investment costs.These values clearly prove that VRB investment for peak shaving in this Wind-PV system is highly profitable.

Conclusions
The proposed strategy for integration of VRB with the Wind-PV system has enabled the intermittent renewable generation unit to be considered as a dispatchable firm power-generating unit.The renewable system has been scheduled based on the forecasted wind speed data with strategic redistribution of generated power to meet peak loads using a VRB battery.The VRB size required to meet this scheduled dispatch has been evaluated using a meta-heuristic optimization algorithm.The VRB has been modeled for two cases as single and dual-battery models in Matlab/Simulink

Figure 1 .
Figure 1.(a) Wind-Photovoltaic-Vanadium Redox Battery (VRB) Hybrid Renewable Energy System (HRES) (b) Subsection showing architecture for Single and Dual VRB battery modules as proposed.

Figure 1 .
Figure 1.(a) Wind-Photovoltaic-Vanadium Redox Battery (VRB) Hybrid Renewable Energy System (HRES) (b) Subsection showing architecture for Single and Dual VRB battery modules as proposed.

Energies 2018 ,
11, x FOR PEER REVIEW 7 of 18 complementary fashion as shown in Figure1b.The simulations are run, and results are compared with the dual battery model.

Figure 3 .
Figure 3. Flow chart for Energy Management Strategy for (a) Single battery case (b) Dual battery model case.

Figure 3 .
Figure 3. Flow chart for Energy Management Strategy for (a) Single battery case (b) Dual battery model case.

Figure 4 .
Figure 4. Membership functions of inputs and outputs for fuzzy controller.Figure 4. Membership functions of inputs and outputs for fuzzy controller.

Figure 4 .
Figure 4. Membership functions of inputs and outputs for fuzzy controller.Figure 4. Membership functions of inputs and outputs for fuzzy controller.

Figure 5 .
Figure 5. (a) Wind speed and Solar irradiance (b) Forecasted Hybrid Renewable Energy System (HRES) power vs Actual HRES power generated (i) 10 days; (ii) for 24 h simulation period.

Figure 5 .
Figure 5. (a) Wind speed and Solar irradiance (b) Forecasted power vs. Actual power generated from the Hybrid Renewable Energy System (i) 10 days; (ii) for 24 h simulation period.

Figure 6 .
Figure 6.Single Battery Case Simulated Results (a) State of Charge (SOC) of Single VRB battery (b) Actual Power delivery by the HRES Pdel plotted against Power Dispatch Dem and Power Generated Pgen.

Figure 6 .
Figure 6.Single Battery Case Simulated Results (a) State of Charge of Single VRB battery (b) Actual Power delivery P del plotted against Power Dispatch Dem and Power Generated P gen .

Figure 7 .
Figure 7. (a) SOC of Dual module BESS1 (b) SOC of Dual module BESS2 (c) SOC of batteries in both cases for 10 h.(d) Power discharged from batteries in both cases for 10 h.

Figure 7 .
Figure 7. (a) State of Charge of Dual module BESS 1 (b) State of Charge of Dual module BESS 2 (c) State of Charge of batteries in both cases for 10 h.(d) Power discharged from batteries in both cases for 10 h.

Figure 8 .
Figure 8.(a) Actual Power delivered by the HRES P del plotted against Power Disptach Dem (b) Power mismatch Curve (Pdel-Dem).

Figure 9 .
Figure 9. (a) Power Charged or Discharged by Dual battery modules BESS 1 and BESS 2 .(b) BESS 1 and BESS 2 discharging voltage.(c) BESS 1 charging and discharging current (d) BESS 2 charging and discharging current for discharging.

Table 2 .
Rule base of fuzzy logic controller.

Table 4 .
Specification of VRB battery.

Table 5 .
Comparison of HRES based on simulated *Note: Simulation period for 10 days.