Energy Storage Power and Energy Sizing and Specification Using HSSPFC

The intermittent nature of renewable sources requires the integration of Energy Storage Systems (ESSs) with appropriate power and energy densities. One of the applications of Hamiltonian Surface Shaping and Power Flow Control (HSSPFC) is to size ESSs for power and energy densities by employing them as sole actuators of Microgrid (MG) systems. This Article provides a comprehensive yet simplified example of utilization of HSSPFC to size ESSs of inverter-based three-phase MG systems under hierarchical control. Here, the distributed Hamiltonian controller is expanded for control of parallel ESSs and power sharing metrics are defined to distribute power between hybrid storage systems according to their power and energy density capabilities. Simulated hybrid ESSs comprising battery and flywheel systems are used as examples to demonstrate the behaviour of the expanded control, verify the power sharing criteria and illustrate ESS design and specification by utilizing HSSPFC.


Introduction
From a design viewpoint, Energy Storage (ES) technologies need to be carefully chosen so that a suitable bandwidth of operation versus power transients is ensured. Long-term transients such as generation fluctuations over time-span of hours or days can be alleviated by using storage systems with high energy density such as batteries [1]. On the other hand, short-term transients such as load step changes may specifically require storage systems with high power density such as flywheels or super-capacitors [2]. The inherent benefits of these technologies can be combined to create a more capable hybrid storage system [3].
The Hamiltonian Surface Shaping and Power Flow Control (HSSPFC) method employs a multi-disciplinary approach to address stability and performance problems by providing solutions for power flow and energy transfer in various physical systems under control [4]. This work aims to illustrate how HSSPFC can be used for specification of battery and flywheel hybrid ESSs. One of the applications of the HSSPFC method is to obtain ES capacity and bandwidth requirements for dc and ac Microgrid (MG) systems [5][6][7][8]. The Hamiltonian control, unified in its original form is a centralized control method, however, decentralized and distributed solutions have been offered in recent works such as References [6,7,9] for dc, and Reference [8] for ac MG systems. From ES control and design perspectives, while Reference [6,8] provide solutions to obtain Energy Storage System (ESS) requirements, the Device Level Controllers (DLCs) have not yet been specified. In other words, it has not yet been shown how the end results interpret into sizing specific storage systems for capacity (energy density) and bandwidth (power density). Moreover, in References [6,8], a series combination of source and the ES element will have its specific limitations. Furthermore, ESS design and control ESSs was reformed according to the HSSPFC method. For the rest of this work, the combination of source, ESS, inverter and its filter/transmission line is referred to as the DGU or the local MG system. The state-space representation of the system demonstrated in Figure 1 is where, v d v q = βλv dc cos(φ) sin(φ) .
The compact HSSPFC representation is where, and, This system can further be analyzed for stability and performance using the HSSPFC [4][5][6]. A single DGU connected to a load was presented for the sake of simplicity, however, there are typically multiple DGUs in a MG system. In Section 4, a d-q droop control [22] will be be leveraged for sharing the load power between parallel DGUs in a distributed form. In the next sections, the overall control systems will be derived based on the model that was presented in Figure 1.

Hamiltonian Control Derivations for Parallel Hybrid Storage Systems
The reference state-space of the system in Figure 1 is obtained as below, This expression consists of reference, nominal and measured values. The feed-forward control expressions for ES and inverter command controls are obtained from the steady-state solution of this reference state-space system. In this section, the ES feed-forward and feedback, and inverter feed-forward control laws for parallel sources and ESSs are obtained.
The series combination of the ES element u with the source voltage v might not always be equivalent to parallel form since at some operation points internal currents may flow between the two components. In this case, the parallel equivalent model would not be feasible. Moreover, in the series topology, ES element might limit the current supply of the source, and introduce more complexity such as requirement of high-side driving for converter control. In this section, the parallel form of the combined source and ESS is realized from the baseline series form. The general aim is to develop control laws for the ES element to satisfy specific constraints. Ideally, the effects of each constraint on ES element can be identified and eventually stacked using superposition principle to obtain the overall ES requirements. Here, the parallel ESS actuator is constrained to include zero steady-state response for the ES element and regulate the inverter dc-input voltage v dc . The inverter dc-input reference voltage denoted by v dc,re f will be used to enable the Hamiltonian controller to determine the ES requirements.
Considering (9), corresponding error state is For the system in (1), the error is defined as Therefore, the ES PI feedback is and the ES feed-forward is Hence, the overall ES feed-forward and feedback control law for the series form in (1) is The series topology of v and u can be obtained from the parallel form and vice versa according to Reference [6] where, v v and u u represent the energy source and storage device for the parallel form. i v , i u , R v and R u are the corresponding currents and resistances. The objective is to keep i u at zero level at steady-state then superimpose the requirements to maintain a constant v dc versus the change in the input source v v .
To have a zero i u , u u should satisfy from which the series form response for u can be obtained as Assuming v v,avg is the average value of the source voltage v v , the dc error e dc can be defined as Considering (15), in series form, this error can be scaled to To remove and compensate for e dc , the ES element u in (14) can be set to From above, u 1 and u 2 are the required responses for u to satisfy the constraints in (16) and (19). Using superposition, from (17) and (20), the overall lumped u should satisfy To enforce the above equation, from (13), v dc,re f can be set to For overall ES control, from (15), u u is controlled according to or i u can be controlled according to Hence, for a current-controlled storage system, the overall ES command is where, u is obtained from (14) and v dc is the measured dc bus voltage.

Primary Control System: DQ Droop Control
The d-q droop control [22] as the primary controller for the source MG in Figure 1 is given as where, V * d,i and V * q,i are the droop voltage settings. R d,i and R q,i are the droop slope settings. I d,i and I q,i are the corresponding currents injected into d-q buses. V db and V qb are the measured bus d-q voltages.
If I d and I q are determined, the real and reactive powers, P and Q can be obtained according to The droop settings V * d,i , V * q,i , R d,i and R q,i are obtained from a secondary controller to improve regulation of V db and V qb which represent the three-phase bus voltage.

Secondary Control System: Power Flow Calculations
Droop control in (26) can operate without communication and load information however, the bus voltage may suffer for lack of communication for long duration of time. A secondary controller can update the droop settings in (26) periodically to adjust the bus voltage level with respect to the load [8]. If I db and I qb represent the load currents, the injected currents from DGUs are shared according to Considering (26), for a nominal bus voltage represented by V * db and V * qb , V d,i and V q,i can be set as so that P and Q are shared according to For larger MG systems with more lines and loads, the droop settings can be obtained from the solution of conventional power flow calculations. In this case, the aim of the power flow problem is to set nominal values for local bus electrical levels and define the power sharing weightings.
It is important to note that the secondary control for Figure 1 is significantly simple since the distributed DGUs measure the same local bus. However, for more complex systems, power flow calculations are typically required. For the expanded system in Figure 2, the simplified secondary control calculations in dq coordinates are where, R line,i and L line,i represent the parameters of the transmission lines that connect corresponding DGUs to the Points of Common Coupling (PCC); hence one major assumption is that the values of the RL lines are known. From above, (31a) to (31d) represent a four equation and four variable system which can be solved to obtain I db and I qb , from which, I d,i and I q,i can be calculated to obtain V d,i , V q,i , R d,i and R q,i . It is important to note that two other equations need to be added to (31) so that the specific droop control method such as slope changing or curve shifting or both are determined [23].
Applying additional constraints will open the grounds for optimization and augmentation of additional loops [24,25] to the secondary controller which are out of the scope of this work. Here, it is assumed that the secondary controller is always present hence the values of two of the four droop settings can be chosen arbitrarily. Throughout the next sections, an assumption is that the secondary control provides local references appropriate for the droop slope changing method; hence, load current and line information are used to define droop settings for bus voltage regulation and current sharing accuracy.

Hybrid Battery and Flywheel System Control and Specification
A hybrid storage system, with its series and parallel battery and flywheel cells requires some form of power sharing scheme. Here, the hybrid system consists of parallel battery and flywheel systems where each have their respective series and parallel cells. The battery system is considered as the primary storage system and the flywheel system compensates when the battery cannot effectively track the reference control signal in (25). The reference signals for individual battery and flywheels cells are where, N p,batt and N p, f w are the number of parallel cells for battery and flywheel systems respectively. i u,re f ,total is the reference current for the overall hybrid system and i u,batt,meas is the measured current injected by the overall BESS. The equivalent hybrid ESS [11] can be represented as where, i batt and i f w represent the injected currents by the BESS and FESS respectively. w cut−o f f ,batt and w cut−o f f , f w represent the estimated minimum cut-off frequencies of the BESS and FESS. (33) presents first-order linear state-space estimations of BESSs and FESSs. Although there are draw-backs of using ESS models in such low fidelity forms (e.g., uniformity in charge and dis-charge dynamics); however, in next sections (e.g., Section 7) it will be shown that higher order systems under their corresponding controls can be optimized to mimic such first-order behavior.
To avoid undesired excursions for the BESS, an additional filter can be added to the system in (33) and the scaling expressions in (32). Hence, the overall system is represented as where, ω cut−o f f , f ilter and i f ilter,re f are the cut-off frequency and the output of the LPF respectively. i ESS is the overall hybrid ESS supplied current. For this system, the overall input is i u,re f ,total . Considering (34), for the overall reference current i u,re f ,total , an ESS designer can choose the cut-off frequencies ω cut−o f f ,batt and ω cut−o f f , f w according to the accepted system tolerance for the dc-link voltage v dc . With the cut-off frequencies identified, one can use it to estimate the battery and flywheel systems as well as the controller parameters. Since ESS devices and the sources are in parallel and v dc is a common entity, the i re f = i u,re f ,total command can be substituted with P re f /v dc . In the next section, P re f will be used in order to determine power and energy capacities of ESSs.

ES Power and Energy Sizing
Specific measures are needed for design of ESSs for power and energy capacities. Typically, there are two aspects to consider; power and energy density requirements. For a single ES element, the ESS should meet both of these requirements, that is having high power and energy densities simultaneously. Should the single ESS fail to meet these requirements, power quality issues may occur and the local MG system may fail. Generally, a hybrid ESS is more flexible in the sense that they can be controlled according to their merits. The disadvantage in this case would be the complexity of the control system compared to a single ESS operation, however, there is more potential for further optimization. A hybrid storage is capable to reduce battery degradation since faster power fluctuations can be compensated by the flywheel system.
Here, it is assumed that the average source power minus the losses meets the load demands. And, the net of power and energy of the ESS over the operation cycle is considered to be zero. In other words, the average power for individual BESS and FESSs over one cycle is zero. In this section, the overall approach is to share the ESS power according to BESS and FESS bandwidth support capabilities by defining appropriate power sharing metrics. The power sizing of the ESSs is determined from the power spectrum of the reference power signal using the Power Spectral Density (PSD) of the ESS response [26]. On the other hand, the energy sizing is performed by analysis of the time-domain requirements.

Power Requirements
If P re f is the reference ESS command signal with an average value of zero, the required RMS value of the continuous power signal is According to the conservation of energy principle, Parseval's theorem states where, S r is the PSD. Therefore, for discrete one-sided power spectrum representation, assuming the magnitude of each frequency content is the corresponding RMS value, the expression in (36) can be re-written as where, P S [n] is the power spectrum corresponding to frequency component n. Here, N is the length of the spectrum. Ideally, one can split the power spectrum by choosing N 0 such that where, P re f ,RMS,ESS is the required RMS of power signal of the ES device for the case the overall PSD is portioned. N 0 (≤ N) represents the maximum frequency that the device can support and it is directly related to the ESS cut-off frequency. γ is a probabilistic metric that should meet (0,1). N 0 and γ are design variables and can be chosen based on the capability, topology and size of the ES system. Considering a system with only one ES element, for a small value of N 0 (or γ), a smaller portion of the overall power can be supplied. On the other hand, for a larger value of N 0 (or γ), a more capable ESS is required. In the case that the device cannot sufficiently support a large N 0 (and the corresponding γ), electrical levels may degrade and the power quality may suffer. LPFs split the power spectrum around a cut-off frequency in a weighted manner rather than a clean split such as in (38). For a more precise approach, (38) should be re-written versus the weighting applied by the filter transfer function. Assuming a hierarchical control architecture, if the amalgamation of the ESS and the filter is represented as H(jn), (38) yields Here, the cut-off frequency of H(jn) and γ are the design variables. If H(jn) with a specific cut-off frequency is present, assuming it is an accurate estimation of the lumped system, γ can be found. On the other hand, if a power share of say 90% is desired, (39) can be used to calculate the required cut-off frequency.
It is important to note that the above expression deals with the power sharing based on the bandwidth of the ESS rather than sharing based on the overall energy capacity. In other words, the energy sharing is dictated by the power sharing. Hence, energy is allocated by the control system based on the priority that is given to each storage device according to its bandwidth. Now, assuming H(jn) is the BESS with the LPF in (34), for a battery/flywheel hybrid storage system, the power sharing can be defined as where, P re f ,RMS,Batt and P re f ,RMS, f w represent the power requirement from the battery and flywheel devices respectively. Here, respective to the BESS, the flywheel device may have higher bandwidth of operation but it is still band-limited and may fail to support the full requested power. Hence, the overall ESS must meet (39) with a sufficiently high γ.
A first-order LPF might not be an accurate representation of the BESS. As a result a filter matching leakage may occur. In Section 7, it will be shown how a more detailed BESS (such as in Reference [11]) can be optimized so that its frequency response matches a first-order LPF.

Energy Requirements
Energy of the overall ESS can be obtained from, The range of the processed energy for a cycle is hence, The required ES range is This expression represents the required capacity of the hybrid ESS. This value is eventually shared when the control system shares the overall power between the battery and flywheel sub-systems according to the their bandwidth of operation. Hence, the above expression can be re-written for each ESS sub-system. Considering a battery system, the required energy range is Assuming the initial energy of the battery is known, keeping the battery charge and discharge range between 20% and 80% of the overall State of Charge (SOC) yields where, (Ah) capacity is the overall energy capacity for a battery. Hence, the minimum overall capacity of the battery can be obtained from (Ah) capacity ≥ 1 0.6 E r,batt .
The above can then be used to determine energy density of BESS. Considering a flywheel system, the required energy range is Assuming the Depth of Discharge (DOD) for the flywheel system over a cycle is 75%, Hence, the flywheel energy capacity should meet

Battery Energy Storage System and First-order LPF Matching and Estimation
In previous sections, BESS and FESS systems were described as first-order LPFs such as in (34). However, the BESS in Reference [11] is clearly a higher order system. Mismatching the frequency responses of both systems may lead to substantial inaccuracies in power sharing based on frequency response. Moreover, choosing a large ESS cut-off frequency may lead to an over-sized BESS from the power density point of view.
The general approach is to choose the order of the LPF system in (34) so that the filter can be a more accurate representation of the BESS. Here however, for the sake of simplicity, the filter is kept as a first-order LPF and the BESS control is optimized to mimic a first-order LPF. In this case, the battery control system PI gains can be considered as optimization variables. The optimization problem is set up as

Minimize
F(k p,batt , k i,batt , k p,u , k i,u ) : where, H BESS (jw) and H f ilter (jw) represent the frequency responses of BESS and the first-order LPF respectively. Generally, H BESS (jw) is obtained from battery with its corresponding Device Layer Control (DLC) system such as in Reference [11]. ω cut−o f f represents the cut-off frequency for both BESS and LPF. ω 0 and δω represent initial frequency and the frequency deviation respectively. kδω represent the step in frequencies for which the frequency responses are matched. (50b) shows the inequality constraints for the PI gains of the nested loops of the BESS model. From (50), the objective function F can be minimized so that an accurate estimate of the LPF is obtained. The overall approach from above is equivalent to fitting the BESS and LPF Bode gain plots. It is important to note that the BESS cut-off frequency can be chosen relatively larger than the LPF so that the minimum power contribution of the BESS shown in (40) is ensured. However, (50) yields an accurate estimation of the LPF. The design of higher order LPFs are left for future iterations of this work.

Hybrid Battery and Flywheel Storage Discussion
Considering ideal first-order band-limited BESS and FESS in (34), hybrid ESS of each DGU in (2) can be sized and specified by following these steps, 1. The system with an ideal ESS controlled under (25) is run and overall P re f , voltage and current levels are obtained. Time domain analysis is performed to obtain the overall energy capacity. 2. From voltage and current levels, the number of series and parallel battery and flywheel cells and individual battery P re f 's are obtained. 3. (37) is used to declare RMS of the power signal in frequency domain as one-sided power spectrum. 4. If H(jn) represents the frequency response of the BESS with the first-order LPF as shown in (34), its cut-off frequency is obtained using (39) and (40). The FESS is similarly sized for appropriate frequency response (power density). 5. Considering (34), system with the hybrid ESS under the specified power sharing is run and time domain analysis (i.e., in Section 6.2) is used to size BESS and FESS for the corresponding energy capacities.
It is important to note, for the reduced-order system in Section 2, under the control in Section 3, and using the hybrid BESS and FESS (e.g., in Reference [11]), the overall energy trade with BESS and FESS will not be zero due to losses. Compensation and sizing ESSs versus system losses can be done by modifying (34), however, it is out of the scope of this work.

Illustrative Examples
The overall system in Figure 2 with parameters in Table 1 is simulated in Wolfram Mathematica, SystemModeler [27] and Modelica [28]. The system includes two DGUs that are connected to a common bus and an RLC load through RL transmission lines. Two examples are presented in this section. For both cases, the ESSs aim to compensate for source fluctuations by enforcing the references that are fed by the local Hamiltonian control system. The arbitrary nominal value for the bus voltage is 240V pk−pk and the secondary control scheme is defined so that the bus voltage always returns to this nominal value. This is a hard constraint for the operation of the system causes the overall MG system mimic a constant power load from source-side perspective. The aim is to maintain the bus voltage and the load power using the sources and the ESSs. It is assumed that the average power value of the source is sufficient to maintain the load and the ESS is utilized only to smooth the power to the average value. For the load voltage of 240V pk−pk , the corresponding nominal d-q voltages V * db and V * qb , are 147V and 0V respectively. Voltage sources v v,1 and v v,2 have randomly chosen dc values around 300V with superimposed uniform random white noise with amplitude of less than 20V pk−pk . A practical constraint to consider is that the source voltages cannot have values less than the corresponding dc bus voltages v dc,1 and v dc,2 . The random voltage component aims to introduce fluctuations so that the contributions of ESSs with different bandwidths are presented for a worst-case scenario for a field deployed MG.
In the first example, load power is solely supported by one DGU and droop control is set so that DGU 2 does not contribute to the load power. This baseline example aims to demonstrate the efficacy of the combination of the Hamiltonian control in (25) and demonstrate the trade-off between the ESS response versus the electrical level regulation of the system. The results are shown in Figures 3-13 for a single ESS element with first-order model in (33a). In the second example, the system in Figure 2, under primary droop control in (26), secondary control in (31) and Hamiltonian control in (25) is simulated. In this example, the steps in Section 8 are followed for load and source profiles shown in Figure 14.
Step 5 however, is re-iterated to utilize the matched BESS and FESS models in Reference [11]; hence, the PI gains of the cascaded BESS controller in Reference [11] are optimized according to (50) to match the first-order LPF.
It is important to note that the aim of ESS sizing in these examples is not to determine BESS or FESS device component parameters. Instead, the intent is to analyze the power split using the sizing criteria in (40) and observe the subsequent energy allocation. Hence, the power and energy sizing are done from the perspective of the Point of Common Coupling (PCC) between the DC-link and individual BESS and FESS subsystems. Sizing an energy storage device for its specific components and parameters under various DLCs and Energy Management Systems (EMSs) requires device-specific design measures, rigorous analysis and addressing multiple trade-offs at device and control levels which are out of the capacity of this work. Hence, the simulation cases are designed to serve the objective of the proposed control rather than to dictate detailed device-specific component sizing.

Single DGU with Constant Load Example
A 6 s simulation of the system is performed for a fixed load value of R b = 2 Ω. The source voltage with superimposed sampled variable voltage is shown in Figure 3. The results for a band-limited ES element with cut-off frequency of ω cut−o f f = 1000 rad/s are presented in Figures 4-8. The maintained three-phase bus voltages and the injected load currents are shown in Figure 4a,b respectively. Figure 5 shows the dc bus voltage and Figure 6 demonstrates the current of the band-limited ES element. Figure 7a,b demonstrate the bus voltage and injected load current amplitude fluctuations respectively. Direct and quadrature current components, i d and i q , represent the injected current to the bus/load and are shown in Figure 7c,d.
The overall dc power of the source and the ESS is presented in Figure 8a. Figure 8b shows dc power is delivered to the load with some reduction due to the transmission line losses. The significant loss is expected since the system is under load of 2 Ω at 240V pk−pk bus voltage. The contributing ESS dc power is demonstrated in in Figure 8c.  Figure 12. It can be seen that there is a trade-off between the ESS size and bandwidth of operation and the voltage fluctuations in dc bus. Figure 13 demonstrates the required normalized capacity against the bandwidth sweep. The normalized capacity is obtained by dividing the required capacity at specific frequency support over the maximum required capacity when an ideal ES element is present. It is important to understand that in this example, the aim is to put the most possible stress on the ESS. Typically, filters such as dc reservoir capacitors are chosen significantly larger than the ones chosen for system in Figure 2. Bus fluctuations that have higher frequency contents can generally be attenuated with such capacitors, however, here such capacitors are down-sized to highlight the behaviour of the control.

Parallel DGUs Example
In this example, the two DGU system in Figure 2 with parameters in Table 1 and under droop control with weighting of 1:2 is considered. The source and load profiles are shown in Figure 14. Considering the guidelines in Section 8, first, the system is run with an ideal ESS model. The corresponding reference powers, P ESS,1 and P ESS,2 are shown in Figure 15. From (43), the overall required energy capacities, for DGU 1 and DGU 2 are 5159.9J and 4529.9J respectively. For BESS power-share of more than 95%, the overall RMS of the power signals are obtained as 2042RMS and 2408RMS for DGU 1 and DGU 2 respectively. Moreover, according to (39), with γ = 95%, the cut-off frequencies are found to be 12.56 rad/s and 18.84 rad/s for filters corresponding to BESS 1 and BESS 2 . With the individual cut-off frequencies known, the control gains of BESSs from Reference [11] are optimized to match first-order LPFs as shown in Figures 16 and 17. The obtained gain values are shown in Table 2.     The simulation is re-run and the BESS with specified control gains are used. The system performance is shown in Figures 18-22. Figure 18 demonstrates individual hybrid current contributions for BESSs and FESSs. It can be seen that significant low-frequency portions of the currents are allocated to BESSs, and the FESSs are requested only for high power and fast fluctuations. As mentioned before, for energy density sizing, the overall injected power (or current) from the ESS device is considered rather than detailed specification of battery and flywheel system components. The obtained BESS RMS contributions are 1957.3RMS and 2306.6RMS corresponding to 95.84% and 95.78% of the total RMS of the power signals. Using (41) and (43), the total required energy from BESS 1 and BESS 2 are 4836.1 J and 4372.7 J respectively which match the expected results with an acceptable accuracy. Figure 19 presents the source powers and overall ESSs dc powers. Individual dc-link voltages, v dc,1 and v dc,2 are shown in Figure 20a,b respectively. It can be seen that ESSs regulate dc voltages versus the source fluctuations rather than the load which is consistent with the results of the previous example. However, in constant load periods, v dc variations are relatively smaller which demonstrate the efficacy of the hybrid ESSs for power bandwidth support. Figure 20c,d demonstrate the dynamics of individual inverter commands (λ, φ) for the two DGUs.

BESS Parameters FESS Parameters
Considering the load-side bus, Figure 21 shows real and reactive power sharing under the droop control. γ 1 and γ 2 from (28) are chosen to be 0.33 and 0.66 respectively. It can be seen that the power sharing is effectively maintained versus the load changes. The maintained three-phase bus voltage is demonstrated in Figure 22a. It is evident that as the sources and load fluctuate, the three-phase ac bus amplitude (denoted by phase a voltage v a ) is kept at 240V pk−pk . Figure 22b shows the corresponding load current while Figure 22c demonstrates the shared three-phase currents which both are consistent with the specified 1:2 power sharing weightings.

Conclusions and Future Work
In this paper, the HSSPFC method derivations, under hierarchical control and for parallel topology of source and ESS were presented. Specific constraints were defined and added to the control law to put the parallel ESS at zero output conditions and compensate for source-side fluctuations. A simple power splitting method for hierarchical ESS control was defined and utilized to size hybrid BESSs and FESSs for the corresponding power densities. The results attest if the LPF and the band-limited BESS frequency responses are matched, the described power sharing based on power spectrum yields accurate results. The performance of a MG system with hybrid ESSs and two DGUs under d-q droop control versus variable sources and load was demonstrated.
Most of the assumptions in this article can be challenged for future work. The control law defined in (25) will further be modified to account for constraints that address various control scenarios such as power smoothing [13], MPPT [15], peak shaving [14] and power curtailment [12]. Moreover, detailed device controllers such as in Reference [29] will be used to specify and integrate series battery module control. Furthermore, design of more detailed Power Management Systems (PMSs) based on Wavelet transform methods [30] can also be investigated. Funding: This research was funded in part by the US National Science Foundation grant number #1541148.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this Article: