Energy Storage Sizing Strategy for Grid-Tied PV Plants under Power Clipping Limitations

: This paper presents an analyses of an Energy Storage System (ESS) for grid-tied photovoltaic (PV) systems, in order to harness the energy usually lost due to PV array oversizing. A real case of annual PV power generation analysis is presented to illustrate the existing problem and future solutions. Three PV modeling techniques have been applied to estimate non-measured non-harnessed PV power to provide an ESS energy and power sizing strategy. Moreover, a control strategy to store or release power from the DC-link, without modifying the Maximum Power Point Tracking (MPPT) strategy, is presented. The results show an estimation of the annual power loss caused by oversizing the PV array. The ESS sizing strategy gives insight into not only the energy requirements, but also the power requirements of the system. Simulation results show that the proposed ESS control strategy is capable of harnessing the extra power without modifying the existing power converter of the PV plant nor its control strategy. presents a grid-tied central inverter PV plant with PV array oversizing, where the available PV power is truncated at the of the inverter (clipped), limiting the exported PV power. Both power curves were normalized to the inverter rating.


Introduction
Perpetration of renewable energy in electric markets has reached an impressive 26.5%, being wind, bio and solar power at the forefront of modern renewables development and integration to electric retail [1]. A crucial parameter when designing renewable energy plants is its load factor, also known as capacity factor or plant (load) factor, which corresponds to the ratio between the generated and rated energy of the plant during a certain amount of time. In UK, PV plants present annual load factors close to 10% [2], which are calculated considering the rated power of the converter. A common practice is to increase annual plant factor by oversizing the power rating of the PV array, with respect to the converter [3]. The ratio between PV array rated power and the inter AC rated output power is known as Inverter Loading Ratio (ILR) [4]; in places with high irradiation variability such as UK, PV array power ILR oversizing can reach as much as 40%, whereas, in places with lower irradiation variability, such as central Chile, oversizing is closer to 15%. Moreover, the continuous drop on PV module prices have encouraged the increase of ILR in PV plants [3]; some authors have even proposed ILR oversizing up to 80% [4].
When an oversized PV array reaches the power rating of the converter, the converter loses the ability to increase its current and therefore is unable to reduce the DC-link voltage and loses the ability to track the Maximum Power Point (MPP). This behavior is called clipping, and it forces the system to waste available PV power. Clipped power is the name assigned to this wasted power. Figure 1 presents a grid-tied central inverter PV plant with PV array oversizing, where the available PV power is truncated at the rating of the inverter (clipped), limiting the exported PV power. Both power curves were normalized to the inverter rating. Generation-demand matching paradigm, where intermittency and high variability of renewable resources play a major role, can be achieved by relying on other systems connected to the electrical network and/or by the addition of an Energy Storage System (ESS). The first solution is not suitable for harnessing clipped power, since it requires the clipped power to be transferred to the electric network through the inverter, which is already operating at its rated power. The latter solution presents a more promising alternative to enhance existing PV plants, enabling them to harness clipped power. This solution has been widely researched as an alternative not only to deal with generation-demand mismatch, but also as means for renewables to provide complementary services, such as load shifting [5], global maximum power point tracking [6] and peak-shaving [7]. Note that the standard location for ESSs is, as shown in Figure 1, beside the transformer (before [8] or after [9]).
Sizing the ESS is a fundamental part of designing a tailored solution to handle clipped power. ESS sizing strategies for PV applications have been previously proposed in the literature: in [10,11], sizing strategies to comply maximum power ramp rate regulation were proposed; in [12], a sizing strategy to provide support for household PV applications; in [5], a sizing strategy to balance the peak and off-peak electricity consumption; and, in [13], a sizing strategy for smoothing power output and storing clipped power at PV plant level. This latter sizing strategy consists on averaging the PV power beyond a certain power limit (clipping level), hence hiding power dynamics to the sizing process. Additionally, the analysis is based on a single sunny day. It must be noted that the power limitation is imposed by contract with the grid operator. In addition, this sizing strategy aims at providing a concentrated solution for a full PV plant, where energy storage is connected at the point of common coupling and inverter ratings are not a limitation for the power exceeding PV plant.
A much wider variety of ESS sizing can be found in the literature related to wind power applications. In [14], a sizing strategy to maximize service-hours per BESS unit cost is presented. The strategy forecasts power generation based in long term historic data and statistical noise; this prediction is then low pass filtered, allowing to obtain an ESS power reference curve, which is later processed by a cost function obtaining the ESS energy rating. Nevertheless, low pass filtering generates phase delay depending on frequency, consequently reshaping the power curve and leading to over or under sizing of the ESS. A sizing strategy to minimize penalties caused by not complying day-ahead power bidding is presented in [15]. The strategy generates 25 initial ESS power references by subtracting bid power from 43-hour-power generation forecasts. The initial references are then presented in a histogram, together with a compliance level, which can be used to generate the ESS sizing. The power generation forecast and bidding strategy are not described in the paper. Moreover, this method considers a 43-hour horizon, which is not ideal for PV systems' daily cycles. A hybrid ESS sizing strategy to comply with maximum power ramp rate regulation is presented in [16]. For this purpose, wind forecast and uncertain load behavior are subtracted, generating a power reference which is later transformed into frequency-domain by Discrete Fourier Transform (DFT). The result is later separated into low, medium and high frequencies, corresponding to the desired power output, the power reference for BESS and the power reference for Supercapacitors, respectively. However, DFT strategy decomposes the full signal into periodic sinusoids losing information regarding the time location of frequencies, therefore leading to a wrong sizing of the ESS. Another strategy to size a hybrid ESS while complying with maximum power ramp rate is proposed in [9]. Here, several historical datasets are filtered by wavelet discrete transform, generating a maximum power ramp rate compliant power curve. ESS power reference curve is obtained by averaging the differences between all original curves and their filtered version. The result is later filtered selecting high and low frequencies as supercapacitors and BESS power references, respectively. The strategy relies on averaging the results, hence masking some dynamic behaviors.
This document presents an analysis of the annual power generated by a PV plant. An analytical model was applied to estimate annual clipping losses. An ESS sizing strategy, based in historic data, was proposed; this strategy considers efficiency of the technology (energy storage technology and power electronics) and provides ESS energy and power sizing, required to recover a certain percentage of the annual clipped power. Additionally, configuration and control strategies were proposed to retrofit an existing PV plant, in order to handle clipped power without modifying the existing MPPT strategy. To validate, at power converter level, the technical feasibility of performing the clipping energy storage service, real PV system and power converters models including control strategies were simulated. The simulations shows specifically that existing central inverter based PV plants can be retrofitted to perform this service (without modifications to the central inverter topology and control). Moreover, the study provides an insight into the daily and seasonal behavior of PV power generation, hence suggesting the advantage of additional usage of ESS, as ancillary services, during idle hours.
To the best knowledge of the authors, the estimation of clipped power, the ESS sizing strategy, the proposed ESS configuration enabling fully usage of a Battery ESS (BESS) and the proposal of a control strategy to harness such power, are novel.
The document is arranged as follows: Section 2 presents a brief description of problem. Section 3 describes the PV model applied to estimate available MPP and a comparison between predicted power and empiric power measurements. Section 4 section presents the ESS sizing strategy. The selection of an Energy Storing Technologies (ESTs), capable of handling clipped power, is presented in Section 5. The control strategy, configuration and simulation of the ESS connected to the PV plant is presented in Section 6.

Problem Description
To emphasize the consequences of oversizing the PV array and highlight the effects of clipping, data from a PV plant located in UK with 39% PV array oversizing is presented in Figure 2. This PV plant has an empiric annual plant factor of 15.43%; if oversizing was neglected a yearly plant factor of 11.11% would have been obtained instead. A year of PV power generation from a central inverter grid-tied PV plant is shown in Figure 2a, where the installed capacity of the PV array is 2 MW, while the DC rating of inverter is 1.54 MW. Power measurements were taken at the DC side of the inverter once per minute during a full year, from 1 October 2016 to 30 September 2017. These measurements have been normalized with respect to the inverter rating. The surface presented in Figure 2b (lateral view of Figure 2a) is equivalent to overlapping all daily DC power generation curves, showing the power limitation, hereafter power clipping, caused by the power-oversized PV array reaching the power rating of the inverter. Seasonal behavior is shown in Figure 2c,d, where yearlong dawn, dusk and daily maximum power generation are presented. These characteristics were exploited to propose alternatives to obtain further usage and revenue from the ESS. Autumnal (Autumn) Equinox, Summer Solstice, Vernal (Spring) Equinox and Winter Solstice are identified in Figure 2c,d with the abbreviations AE, SS, VE and WS, respectively.
Irradiance, temperature and DC power measurements are the only data available for the PV plant located in UK. Since the power limitation applied to the system is imposed by the converter rating, power losses caused by clipping are neither estimated nor accounted for. Estimating those losses would result in a mandatory effort to assess clipping effects.

Proposed Sizing Strategy
The following steps correspond to the ESS sizing strategy applied to enable harnessing the PV plant clipped power.
1. PV modeling selection: Select a model to estimate clipped power. For this purpose, an analytical model depending on irradiance, module temperature and PV module parameters was chosen. The model predicts the MPP, which is later limited (clipped) at the rating of the converter (clipping level) and compared to the DC power measurements from the PV plant. Four error metrics and additional features were considered to show the accuracy of the model. 2. Clipped power estimation: Subtract the measured DC power from the selected model predicted DC power (not clipped) to estimate the clipped energy resulting from clipping the power curve. These calculations generate annual clipped power and energy curves with one minute spanning. 3. ESS sizing: Analyze annual clipped power and energy curves. For this purpose, a statistical analysis of daily clipped power and power-limited daily clipped energy is presented. Power and energy sizing are performed by setting a recovery ratio of the annual clipped energy; from here, ESS energy and power ratings are obtained.

PV Plant Model
There are several modeling techniques that can be used to estimate maximum power generation from a PV plant, such as single diode circuital model [17], double diode circuital model [18], artificial neural networks [19] and analytic model [20]. An analytical modeling method was considered to estimate power loss due to clipping in a grid-tied PV plant.
To validate the model, one year of data ( Figure 2) were compared to the power predicted by the model. Clipped DC power in kW, irradiance in kW/m 2 , module temperature in • K measurements and PV module data sheet parameters were available. Uniform irradiance and temperature conditions among PV modules were considered to estimate the output power. A description of the analytical modeling method and a quantitative comparison, with the empiric data, is given below.

Analytical Model
This mathematical model predicts the MPP (P mpp in W) as a function of the irradiance (G in W/m 2 ) and the module temperature (T in • K).
where k p , ∆T, A, η, η mppt , N ms and N sp are, respectively, the temperature coefficient of P mpp in %/ • K [21] (or maximum power correction factor for temperature [20]), module temperature difference between the module temperature T and the STC module temperatureT stc in • K (∆T = T − T stc ), area covered by PV cells in m 2 , STC module efficiency, MPPT efficiency [22], number of modules in series in each string and number of strings in parallel. A similar alternative is presented in [20] where A · η is replaced by P mpp stc /G stc .

Model Validation
The analytical model used to predict the clipped power was chosen, since it presents a low modeling error below clipping, presents a low model error compared data sheet stated STC, presents a low computational cost, does not require an optimization stage (neither parameter identification nor training) and is conceptually simple. The model prediction was clipped at the DC power rating of the converter to match the maximum DC power level (clipping level), and then the error between the clipped prediction and the measured DC power was calculated. The technical details applied to model the PV plant are presented in Table 1. PV modules correspond to the Jinko model JKM260-PP.
The error metrics applied to validate the PV plant model were Normalized Root Mean Squared Error (NRMSE), Normalized Mean Absolute Error (N MAE), Pearson linear correlation factor (Pearson) and Normalized Root Mean Squared Error Fitness (NRMSEF). For the first two error metrics, NRMSE and N MAE, the optimal value is 0%, while, in the second pair of error metrics, Pearson and NRMSEF, the optimum is 100%. Equations (2)-(5) correspond to the mathematical description of the metrics applied to analyze the error between the DC power measurement (X i ) and the clipped DC output power predicted by each model (X i ). Variables µ X , µX, σ X and σX in Equation (4) correspond, respectively, to the mean of X i andX i , and the standard deviation of X i andX i . N corresponds to the number of samples. Normalized metrics were measured respect to the clipping power level (P clip ).  Table 2 summarizes the error metrics applied to the model. The analytical model presents a low calculation time, a simple approach and small errors metrics. Therefore, this method was chosen to predict PV clipped power. It is also useful for predicting the clipped power in real time.
Aging and soiling effects were not considered in the previously described analytical model, nonetheless the model enables a straight forward update by fitting the STC efficiency term (η) in Equation (1).

ESS Sizing
ESS are composed by a bidirectional power converter and an EST, such as batteries, super capacitors, flywheels, CAES, HPS, among others [23]. Energy restrictions are imposed by EST, while charging and discharging power ratings depend on both EST and the bidirectional converter. Some ESTs, such as batteries, usually present different charge vs discharge power ratings. In this application, ESS charge power rating is limited by the clipped power.
PV power generation is strongly dependant on solar irradiance and the module temperature. The PV module temperature behaves as a low pass filter of the incident irradiance, with an equivalent time constant of a few minutes depending on the wind speed [24]. In addition, according to Vernica et al. [25], the PV plant output power can be modeled as a low-pass filtered version of the solar irradiance, where the cut-off frequency of the equivalent filter is determined by the area of the PV plant. Specifically, for the PV plant being analyzed in this work, the cut-off frequency of the equivalent filter is 0.01 Hz. Additionally, most grid code requirements related to power fluctuation (maximum power Ramp Rate) regulate power fluctuations per minute [26]. Moreover, standard PV plant available data range in sampling times between 1 and 30 min. Therefore, a data sampling rate of 1 min was selected to calculate the clipped energy and the ESS sizing strategy presented in this work.
Clipped power (P clipped ) is calculated according to Equation (6), where P mpp and P pv are, respectively, the predicted PV power generation from Equation (1) and measured PV power, both in kW. The clipped energy in kWh is estimated according to Equation (7), where dt corresponds to the sampling time in min (1 min).
ESS for clipping in PV plants are designed for a daily use cycle (dawn to dusk), which means there is no need to store energy for more than one day, and therefore the stored energy is completely depleted before a new day cycle.
The following efficiencies were considered in the sizing of the ESS: a single stage DC-DC converter with 97% efficiency [27] (94.09% round-trip efficiency), new commercially available batteries present a round-trip efficiency of 95% [28,29] and DC/AC inverter has a nominal efficiency of 97%. Therefore, the energy passing through the ESS (PV to ESS and ESS to Grid) would experience an efficiency of 86.70%. Figure 4a shows a histogram of the maximum recoverable daily-energy-loss (Ê del ) due to clipping per amount of days of occurrence during a year, and two overlapped curves showing the accumulated annual energy loss (due to clipping) and the total recoverable annual energy (including efficiency of the PV-ESS-grid system of 86.7%). Both curves are function of the maximum recoverable daily-energy-loss and were normalized respect to the annual energy loss (33 MWh).
The power rating analysis is shown in Figure 4b, for the case where four ESS power rating design criteria are depicted as a function of the recoverable-daily-energy-loss. All criteria include a 97% efficiency of the DC/DC power converter. The criteria C 1 to C 4 correspond respectively to maximum recoverable-daily-energy-loss in MWh (Ê del ), average recoverable daily-energy-loss in MWh (E del ), mean plus standard deviation of recoverable daily-energy-loss in MWh (E del + σ Edel ) and mean plus two times the standard deviation of recoverable daily-energy-loss in MWh (E del + 2 · σ Edel ).  As an example, to recover 80% of the annual energy lost due to clipping, an ESS of 600 kWh is required to store the maximum daily-energy-loss of (Figure 4a), which considering criteria C 1 (from Figure 4b) leads to a power rating of 200 kW.

Energy Storage Technology Selection
ESS are formed by an EST and its power converter. This section focuses on selecting a suitable EST alternative to store clipped energy. According to the type of energy conversion and the nature of the stored energy, ESTs may be classified as electric, chemical, mechanical and thermal [30], as shown in Figure 5. The main characteristics of ESTs, relevant for clipping, are summarized in Table 3 [31][32][33]. For a detailws description of each ESTs from Figure 5, please refer to [33,34].  Based on round-trip efficiency and maturity level, the best EST alternatives for handling clipped power are SC, LiIon batteries and FES. Nonetheless, energy cost of FES doubles the energy cost of SC and LiIon batteries. Moreover, installation costs are not included in the table, although they depend on the weight of the equipment that needs to be transported. In addition, considering energy density, we have have concluded that LiIon batteries are an adequate EST technology to be applied to handle clipped power. The following section shows a simulation of a LiIon BESS applied to a PV system, for validation of the sizing methodology.

Simulation Results
This section presents the control strategy, configuration and simulation of the ESS connected to the PV plant. BESS was selected as EST to be emulated during the simulation, though the previously described sizing strategy is valid for any EST. BESS was selected due to its modularity, which allows retrofiting each PV inverter at their DC side, and because BESS combine both energy and power density required for this application.

Configuration
The original PV system consists on a central inverter grid-tied PV plant, connected to the grid through a standard two level voltage source inverter (2LVSI). The addition of the ESS, formed by a battery pack and a single-stage DC/DC converter as in [35], merged to the DC-link of the PV system, enables clipped energy to be stored and released according to the system requirements and limitations. Figure 6 shows a simplified version of the full configuration; colored arrows illustrate the possible power paths. The grey arrow corresponds to a continuous power flow (from PV to the grid), while the green (PV to ESS) and purple (ESS to the grid) arrows correspond to excluding and non necessarily continuous power flows.
This configuration has a single stage DC/DC converter. However, as shown in Figure 6, the system can be implemented with a seconds DC/DC stage (marked as Optional). The difference between both alternatives lies in the percentage utilization of the capacity of the battery pack. A single DC/DC converter performing a standard battery charging strategy applies constant current (CC) charging mode until reaching a certain State of Charge (SoC), usually around 85%. The addition of a second DC/DC stage allows the control of the output voltage of the converter, hence enabling to transition into constant voltage (CV) charging mode, and therefore allowing full charging of the battery. Storing clipped power requires withdrawing power from the DC-link, since the MPP tracking voltage is imposed by the inverter, thus clipped-storable power must be controlled through current.
The simulation results are focused on the control strategy that enables performing power extraction and injection to the DC-link, without modifying the existing control strategy of the inverter. This strategy depends on the SoC of the battery pack, which was estimated through standard Coulomb counting method [36]. This estimation is independent of the battery model applied since it relies in the output current of the battery pack, enabling the utilization of an ideal circuital model [37], to emulate the behavior of the EST.

Control Schemes
The central inverter was considered to be controlled through a standard voltage oriented control strategy with MPPT reference, as depicted in the control scheme in Figure 7 [38]. Since the objective is to retrofit an existing PV system, the addition of the ESS is performed at the DC-side of the inverter and does not consider modifying the PV system or its control strategy. It must be noted that the addition of an ESS at the point of common coupling would not allow harnessing the power loss caused by power clipping, since the inverter power limitation would remain an issue, which is not the situation in other similar cases that can benefit from the same ESS sizing strategy, such as power curtailment, where the power limit is imposed by the control strategy and not by the rating of the converter.
Retrofitting forces the ESS control strategy to extract or inject power to the DC link without modifying the existing MPPT strategy of the PV inverter. For this purpose, a tailored ESS control strategy, shown in the flow chart in Figure 8, was designed. P pv , P inv , P ess , P mpp , P clip and P pre correspond, respectively, to PV plant output power (P pv = i pv · v pv ), inverter input power (P inv = i inv ·v pv ), ESS power (P ess = i ess · v pv ), estimated MPP according to Equation (1), clipping level, and pre-clipping level. This latter parameter is used as margin to perform power injection or subtraction to/from the DC link (P pre − P mpp ).P ess (< 0) andP ess (> 0) correspond to the minimum and maximum ESS power. Cases I-V match those in Figure 9. This control strategy can be adapted to handle any EST. For this purpose, estimation of the ESS SoC value (SoC) in the control scheme of Figure 8 must be replaced by an estimation of the available energy in the applied EST.
The ESS output current reference (i * ess ) is provided to a standard PI controller, generating the modulation index and switching pattern according to the power flow direction (storing or releasing energy). The SoC estimation was performed through standard Coulomb counting method [39].

Simulation
A simulation of the full system, formed by a 2LVSI operating as central inverter grid-tied PV system, and ESS formed by a Battery ESS (BESS) and an Isolated Bidirectional Boost Converter [35], is shown in Figure 9. The control strategies previously described in Section 6.2 were applied. It must be noted that the ESS was undersized to allow its control system to be tested in all possible scenarios and display those results in a single figure; this was performed considering the dynamics of the systems and an ESS operating range of 20-80% of the SoC. Several irradiance conditions were tested in the simulations, all with PV cell temperature of 298.15 • K.
The PV plant parameters correspond to those presented in Table 1, while Table 4 presents the parameters of the converters and ESS. To calculate per-unit values in Figure 9, the following base values were considered: G base = 1000 W/m 2 , P base = 1540 kW, i base = 3.767 kA and v base = 800 V.  From 0 to 0.03 s (Case I), the system is operating at a power rating below P clip (1540 kW) and the ESS is at a SoC level of 20%, hence ESS power reference is zero. From 0.03 to 0.14 s (Case II), a step in solar irradiance causes the PV power (P pv ) to reach P clip , since SoC < 80%, P * ess power reference is set to P * ess = P pre − P mpp (and i * ess = P * ess /v pv ), which reduces power flowing through the inverter (P inv ) enabling MPP tracking. From 0.14 to 0.20 s (Case III), a step down in solar irradiance causes the PV power (P pv ) to be lower than P clip , but, since P pre < P pv < P clip , the ESS power reference (P * ess ) is set to zero. From 0.2 to 0.25 s (Case II), a second step up in irradiance generates for the PV power to surpass P clip , the system behaves exactly as from 0.03 to 0.14 s. Once the ESS reaches 80% of the SoC, the ESS stops drawing power and the inverter losses momentarily the capability to track the MPP (Case IV from 0.25 to 0.28 s). Note that the ESS was undersized, to show the behavior of the system when reaching its maximum and minimum permitted SoC. From 0.28 to 42 s (Case V), a step down in solar irradiance causes the PV power to be below pre clipping level (P pv < P pre < P clip ), and the ESS releases power towards the DC-link (P * ess = P pre − P mpp ) and through the inverter into the grid. From 0.42 to 0.5 s, the ESS reaches 20% of the SoC, hence power flow from ESS towards the DC-link is stopped and the system goes back to operating in Case I; in this part, grid currents show high harmonic content due the comparison of power flowing to the grid and power required to generate voltage steps to perform MPPT. In a real case, the MPPT period is longer, hence harmonic content would be lower.
Note that the instantaneous power balance is given by Equation (8), where the corresponding terms are given by Equations (9)- (11). This can also be verified in power curves shown in Figure 9c,d. Variables v gd , v gq , i d and i q in Equation (9) correspond to the grid voltage and current in rotational coordinates.

Conclusions
The work described in this paper focused on the effects clipping, caused by oversizing a PV array with respect to the power rating of its converter. An insight into the benefits and drawbacks in terms of annual plant factor and energy losses of this commercial practice is presented. An analytical model was used to predict the annual power loss due to clipping; the validation of the model was performed by comparing the predicted power, limited at the clipping level, with real data over a time horizon of one year sampled every minute. An ESS sizing strategy based on recovering annual clipping losses is proposed. The strategy uses a power and energy approach, which considers statistical data to select the best fitting ESS ratings. For the analyzed PV plant, the sizing method determined that a ESS of 600 kWh per central inverter enables the retention of 80% of energy that would be lost due to clipping. From a power perspective, each central inverter, rated at 1.4 MW, requires being retrofitted with a 200 kW DC-DC converter for the ESS to enable the aforementioned clipped energy storage. In relation to the EST, LiIon based BESS was selected to perform a validation of the ESS and control system due to comparison of several criteria (LiIon have 10 times higher energy density than FC, about half the cost of FES among other values). Simulation results at power electronics level of the ESS, the central inverter, the grid connection, and their control show that the selected ESS can be retrofitted to existing central inverters, and provide clipping energy storage while still performing properly (perform MPPT, store/deliver energy, retain a controlled DC link and inject energy to the grid with high power quality).
The proposed methodology and analysis can be applied to determine, depending on their plant measurements and parameters, the size of a ESS to retrofit their PV plant for clipping energy storage. This information is necessary to perform an economic assessment, and how it can impact the levelized cost of energy (LCOE), and assist in the decision-making process. In addition, the sizing methodology can be adapted to perform new analysis on other ESS applications such as: load shifting, power curtailment, frequency and voltage regulation, base load generation, capacity firming, etc., providing further value for the ESS to the retrofitting of a PV plant. Furthermore, provided the proper model and data, it can also be extended to other renewable energy sources, such as wind energy and ocean energy.
Author Contributions: N.M. conceived, designed and applied the sizing methodology, and wrote the original draft; S.K. provided insight, resources, and supervision and reviewed the manuscript; P.Z. provided insight, resources, supervision and analysis; P.W. provided relevant key theoretical and technical suggestions and reviewed the manuscript; G.B contributed with mathematical modeling tools and insight; and F.G. provided experimental data, problem conceptualization and practical insight.

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

Abbreviations
The following abbreviations are used in this manuscript: