1. Introduction
The concept of energy transition has gained significant momentum in recent years as reported by numerous authors. Potrč et al. [
1] and Erat et al. [
2] explore the gradual transition to a renewable energy system in the coming decades, while Griffiths [
3], Haas [
4], and Zhang [
5] offer a political and economic perspective on the energy transition challenge. This increased attention is largely due to the urgent need to reduce greenhouse gas emissions, especially those originating from fossil fuel-powered plants. Muryshev et al. [
6] analyse the correlation between changes in global temperature and atmospheric CO
2 content and Mohammed et al. [
7] study the effects of anthropogenic greenhouse gas emissions. The rising concern about global warming highlights the importance of addressing this issue promptly.
Exploring various renewable energy solutions, wind power emerges as a leading option due to its sustainability and vast deployment potential. The recent growth of wind energy production has been remarkable, leading to significant advancements in wind turbine technology. To name a few, Chen et al. [
8] investigate the reliability of offshore wind farms, Sharma et al. [
9] propose an open-source CFD model for wind turbines, Hassanpour et al. [
10] and Ma et al. [
11] suggest optimal design methods for wind farms, and finally Huang et al. [
12] investigate the effect of the wind turbine installation on the power characteristics. Moreover, regarding case studies, Li et al. [
13] provide an overview of offshore wind turbine installations in China, while Boudounit and Saifaoui [
14] propose a wind farm design approach based on the Dakhla Site in Morocco. Finally, Poudyal [
15] analyses potential strategies for addressing the energy crisis in Nepal.
Wind turbines capture the wind’s kinetic energy and convert it into electricity driving an AC generator. The generated AC current is then transformed into DC by a converter, suitable for battery storage. An inverter subsequently converts the DC electricity back to AC to supply power to the grid or directly power AC loads. Pan et al. [
16] propose a battery/supercapacitor hybrid energy storage, while Reddy et al. [
17] study a RBFN controller grid-tied wind energy conversion system.
This configuration facilitates the attainment of nominal frequency tolerances exclusively through solid-state power electronics while also ensuring a degree of power supply continuity by utilising battery storage. However, this conventional approach faces inherent challenges like the energy dissipation in solid-state power components and the costs and environmental concerns surrounding battery storage systems. To address these limitations, this study presents an alternative system configuration that eliminates the battery system or significantly reduces its size. A pumped hydro storage system employing a reversible hydraulic pump (referred to as a Pump as Turbine) is utilised for energy storage. Regarding the
PAT model, Marre et al. [
18] and Stefanizzi et al. [
19] provide experimental studies for performance prediction, while Liu et al. [
20] offer a theoretical model. On the other hand, Rossi et al. [
21] describe a case study of an aqueduc
PAT application. In this system configuration, when the wind energy exceeds the demand, the
PAT, in pump mode, stores this excess in the storage reservoir. Conversely, when wind speed is insufficient, the
PAT can operate in turbine mode to assist in meeting the wind turbine’s external load. This alternative plant layout offers numerous advantages like decreased installation and operational costs, prolonged
PAT lifespan, and diminished waste management expenses. Nevertheless, this alternative plant configuration necessitates a suitable sizing of the reversible hydraulic machine to maximize the wind turbine’s performance under fluctuating wind speeds and external loads.
A comprehensive sensitivity study of a multi-objective analysis applied to a wind farm integrated with hydro storage is presented in this paper. In [
22], Fornarelli et al. describe a multi-objective sensitivity analysis of a thermal energy storage system, while in [
23] they propose a multi-objective optimal design method for a thermal energy storage system. The optimisation problem seeks to maximise three key objectives:
the net change in the mass of water (potential energy) stored by the PAT in the pumped hydro system, ();
the captured wind energy, ();
the extra torque provided by the PAT in hydraulic turbine mode in order to handle heavier loads ().
However, since wind speed and external electric generator load are inherently random (stochastic), the optimisation results should consider a variety of temporal patterns for both variables. Therefore, in the present study, deterministic time series data for wind speed and external load generated based on Weibull distributions are first used since they provide the best representation of the frequency distribution of both wind speed and load. Akgül et al. [
24], Wais [
25] and Shin et al. [
26] provide alternative statistical methods for the Weibull distribution applied to wind energy using different approaches, while Jahan et al. [
27] analyse Weibull and truncated Weibull distribution of wind speed in Bangladesh. Stochastic components are then superimposed on these time series to evaluate their influence on the system’s behaviour. To achieve this goal, a multi-objective optimisation approach is employed not only to the initial specific time trends but also to the perturbed ones.
The wind power system dynamics was simulated using the MATLAB/Simulink software suite R2023b, which allowed the construction of a block diagram of the wind turbine connected to the
PAT via a constant ratio gearbox as illustrated in
Figure 1 and
Figure 2. It is noteworthy that the electric generator model was intentionally omitted from the simulated power scheme for simplicity. Nevertheless, the power plant model incorporates a time-varying torque load to represent the electric generator’s demand on the wind turbine shaft. Due to its reversible operational mode, the Pump-As-Turbine (
PAT) system was simulated both in hydraulic turbine mode and in centrifugal pump mode.
The forthcoming paragraphs delineate such an analysis scheme, and subsequently, in the Results section, the outcomes of the conducted analyses are presented.
2. Wind System Theoretical Model
A comprehensive model of a wind power system was developed using the MATLAB/Simulink software suite. As depicted in
Figure 2, wind speed and wind turbine angular speed serve as inputs to the wind turbine component, which generates a driving torque directed to the shaft block.
The shaft block ensures that the wind turbine system remains in dynamic equilibrium. This is achieved by accounting for the driving torque generated by the wind turbine blades, the resistant torque that represents the load from the electric generator, and the torque originating from the
PAT. The generator load torque was modelled with a time-varying trend based on a Weibull distribution to simulate varying user demands, as supported by Akgül et al. [
24], Wais [
25], and Shin et al. [
26].
When the available wind energy exceeds the demand, the
PAT, in pump mode, stores this excess (in the form of potential energy) by pumping a specific quantity of water into an elevated reservoir. Instead, if the external load exceeds the maximum torque achievable by the wind turbine, the
PAT changes its operation and functions as a hydraulic turbine. Conversely, if the load approaches the maximum power point (
), then the
PAT is mechanically decoupled from the turbine. Of course, in order to implement such a
PAT management strategy, an algorithm capable of tracking the maximum power point of the wind turbine is required. However, in this work, the algorithm is merely postulated and not implemented, primarily because the main objective of this research is to evaluate which
PAT systems exhibit sustained performance in response to perturbations introduced to temporal patterns of the wind velocity and the AC generator torque. However, S. Yan et al. [
28,
29] developed some interesting works based on the Memory-Event-Triggered Fuzzy control algorithm for maximum power point tracking in wind turbines. Finally, it is worth noting that in the described approach, a constant value of the wind turbine to
PAT speed ratio (design variable) is chosen.
2.1. Wind Turbine
A three-bladed horizontal-axis wind turbine (HAWT) is selected for analysis purpose in this study. The Blade Element Method (
BEM) to describe the aerodynamic turbine blade element interactions is adopted. This method discretises the flow field into stream-tubes to assess both axial and tangential momentum conservation as detailed in Hidayanti et al. [
30] and Laalej et al. [
31]. By referencing the airfoil aerodynamic properties, such as lift and drag coefficients, as can be found in Wolfe et al. [
32] and Coder et al. [
33], the Blade Element Method (BEM) determines thrust and torque forces acting on each infinitesimal blade element (
) at any radial distance (
r) [
34]:
in the context provided,
denotes air density, while
a,
and
indicate the induced axial and tangential velocity coefficients and inflow angle, respectively. Moreover,
A represents the swept area, while
denotes the total blade area; the solidity ratio
is given by
. Finally,
denotes the undisturbed wind speed and
represents the rotor angular velocity. Thus, from Equations (
1), the torque (
) and power (
P) generated by the wind turbine can be expressed as follows [
34]:
where
e
are the minimum and maximum blade radius values, respectively. Using Equation (
2), the power coefficient [
34],
, of the wind turbine can be determined in the following manner:
The following are the main data of the wind turbine used in this study: diameter
m, moment of inertia
, and maximum tip speed ratio (defined as
)
.
Figure 3 shows the power coefficient,
, versus tip speed ratio,
, based on BEM theory.
2.2. Pump as Turbine
The Pump as Turbine is a centrifugal pump capable of operating as a hydraulic turbine by reversing the direction of rotation of the impeller and consequently the fluid flow direction. Reversible centrifugal pumps, called Pump as Turbine, offer a unique advantage. By simply reversing the impeller’s rotation, they can be converted into hydraulic turbines, changing the direction of fluid flow. The adoption of
PATs offers lower costs and reduced maintenance requirements compared to typical hydraulic turbines. However, their implementation, especially in turbine mode, is made more difficult due to the lack of data of
PAT manufactures. As previously observed, Marre et al. [
18] and Stefanizzi et al. [
19] provide an experimental-based predictive performance analysis, while Liu et al. [
20] offer a theoretical one. Furthermore, Liu et al. [
35] and Stefanizzi et al. [
36] propose two different optimal design approaches for
PAT in more recent works for specific applications. One of the simplest methods to predict the performance of such hydraulic machinery when used as hydraulic turbines, also considered in this study, utilises flow rate
, efficiency
and head
as a key parameters at the Best Efficiency Point (
BEP) of the reversible centrifugal pump (pump mode) and then calculating the corresponding parameters when operating as a hydraulic turbine. Several studies report semi-empirical correlations capable of approximating these associations. Conversely, other studies have proposed alternative relationships based on the pump specific speed as can be found in Binama et al. [
37] and in Bogdanović et al. [
38].
This work adopts the approach described in Rossi et al. [
39] to evaluate the operational characteristics of the reversible hydraulic machine in turbine mode through empirical formulas, applicable not only at the maximum efficiency point (BEP), but also for various operating conditions. This procedure requires as the initial step to select the design working point (
BEP) of the
PAT operating as a pump. Moreover, to determine the hydraulic turbine(
PAT) curve, it is essential to have data regarding the corresponding centrifugal pump characteristic map. Fortunately, this requirement is often met since the pump manufacturers typically include this data. In this specific study, a generic centrifugal pump map is considered. This allows the derivation of the operating maps of any pump simply by scaling the values of the design flow rate, the design head, the efficiency and speed value curves according to their corresponding
BEP values. This approach facilitates the execution of a multi-objective sensitivity analysis for four different centrifugal pump families, each characterised by distinct specific speeds.
The operational maps of the reversible hydraulic machine (modeled in Matlab-Simulink) are presented in
Figure 4 for the pump and the turbine mode, respectively. The four examined
PAT families are listed in
Table 1.
3. Multi-Objective Optimisation
Numerous authors have discussed the theoretical framework of multi-objective optimisation analysis and the concept of Pareto optimality (non-dominated solutions) as reported by Periaux et al. [
40] and Caramia et al. [
41]. Single-objective optimisation involves finding the absolute minimum or maximum value of a function with one or more variables. This function represents the single objective that needs to be optimised. In this case, there is a clear best solution that optimises the single objective. On the other hand, multi-objective optimisation deals with problems where there are multiple, often conflicting, objectives to be optimized simultaneously. These objectives are represented by a mathematical function with multiple variables. Since achieving the best outcome for all objectives simultaneously is often impossible, the goal is to detect the Pareto front, which represents the set of non-dominated solutions. These solutions represent the best trade-offs between the different objectives.
Considering
as the vector of design variables and
as the vector of corresponding output performance variables, we let
and
be two distinct design variable vectors with corresponding objective function vectors
and
. For simplicity, let us consider the case where all independent objective functions must be maximised [
42]. It follows that if
the design vector
is identified as dominated by
. Conversely, if
then the solution
is labelled as non-dominated by
. Finally, if
the two solutions do not dominate each other. Equations (
4) and (
5) enable the identification of non-dominated solutions by those dominated, while Equation (
6) concerns the dominance relation between two non-dominated solutions. Employing Equations (
4)–(
6), dominance for all the solutions under consideration can be assessed, and therefore it is possible to determine a first approximation of the
Pareto front. It represents the set of optimal, non-dominated solutions, i.e., those solutions for which there is no solution that is simultaneously better for all the objectives considered. However, there are more efficient strategies rather than an exhaustive evaluation of dominance relations (Equations (
4)–(
6)) across all conceivable design variable values which allow to identify the
Pareto front. One of the most established and effective methods is Design of Experiments (
Doe). It facilitates the determination of output performances
by systematically varying the independent input variables (
) following a predetermined scheme, like a Cartesian grid approach as adopted by Dambrosio et al. [
43] and Dambrosio et al. [
44].
However, applying dominance checks (Equation (
4)) only to
DoE points within the design variable space provides an initial estimate of the Pareto frontier. To achieve a more accurate representation of the Pareto front, a multi-objective optimisation approach is required. This process is focused on finding, within the design variable domain, the directions that lead to improvements in performance variables based on the specific optimisation problem, as in the methodology proposed by Dambrosio et al. [
23]. Nevertheless, since wind speed and external electric load are inherently random (stochastic), the optimisation results should account for various temporal patterns in both variables. Therefore, this investigation first considers fixed temporal patterns for both the stochastic variables (illustrated in
Figure 5), as previously implemented by Dambrosio, but limited to four
PATs [
45]. Based on these reference profiles, random perturbations of
with respect to the baseline value, characterised by a pseudo-random number with normal distribution, are added to simulate realistic operating conditions. It is important to emphasise that a single perturbation of the wind velocity and external load is not sufficient to adequately represent the complexity of the operating conditions. For this reason, three distinct variations for both variables are examined, appropriately compared with the initial reference cases. The temporal variations of reference wind speed and external load, along with their perturbed counterparts, are graphically illustrated in
Figure 5,
Figure 6,
Figure 7 and
Figure 8. The reference profiles are represented by purple lines, while the newly generated profiles are shown in blue.
4. Results
This paragraph summarises the key findings from a multi-objective approach aimed at evaluating the robustness of the Pareto front for a wind power unit incorporating a pumped hydro storage system with a
PAT system (depicted in
Figure 1 and
Figure 2). This analysis, conducted using the MATLAB/Simulink software suite, aimed to maximise three wind system performances: the net change in the mass of water (potential energy) stored in the pumped hydro system, the captured wind energy, and the extra torque provided by the
PAT in hydraulic turbine mode. However, the main objective of this research was to evaluate which
PAT systems exhibit sustained performance in response to perturbations introduced to temporal patterns of the wind velocity and the AC generator torque that can be considered as stochastic variables. In the current study, the design variables considered during the analysis were the pump (
PAT)-specific speed (
PAT geometry), the wind turbine to
PAT speed ratio
, and the reversible hydraulic machine head. Although the electric generator is not included in the MATLAB/Simulink model, a representative torque load was implemented to simulate the wind system’s behaviour. Furthermore, concerning the time history profiles for both the stochastic variables (the wind velocity and the AC generator torque), the current study initially considered reference time trends, illustrated in
Figure 5. Then, based on these reference profiles, random perturbations were added to simulate realistic operating conditions. For this reason, three distinct variations (illustrated in
Figure 6,
Figure 7 and
Figure 8) for both variables were examined. The present discussion is inspired by a previous work of Dambrosio [
45] concerning the performance analysis of the wind power plant, but limited to four
PATs. The present analysis considered twenty
PATs organized into four families, each consisting of five elements; each
PAT family is characterised by a given specific speed value. The design conditions for these twenty systems are provided in
Table 2 below.
An individual analysis of each
PAT system was conducted to determine the design variable values that correspond to the non-dominated solutions on the Pareto front (Equations (
4)–(
6)). This approach enabled the identification of the best approximating Pareto front among the analysed
PATs. This was accomplished employing a
BOTTOM-UP approach, from a singular study of each
PAT to a global analysis among all
PAT families. Below, for brevity, only the results pertaining to the overall dominance study across the four
PAT families are presented.
From
Figure 9,
Figure 10,
Figure 11 and
Figure 12, it can be noted that
PAT Family 2 is completely dominated, while a significant number of non-dominated solutions belongs to Family 3 (yellow dots). In other words, considering the reference wind speed and external load time history (
Figure 5), the most performing
PAT family is the third one. These results suggest focusing further investigation only on Family 3 elements; specifically, the additional analysis aimed to evaluate the
PAT performances as the wind speed and external load profiles assume the trends shown in
Figure 6,
Figure 7 and
Figure 8.
Then, the dominance studies were performed for each
PAT and globally for the entire third family. For the sake of brevity, only the comparison between the best design solutions of the third
PAT family obtained in the four study cases (reference profiles
Figure 5 and perturbed profiles
Figure 6,
Figure 7 and
Figure 8) was reported. Finally, following further dominance analysis, the new Pareto front obtained (pink dots) is reported in
Figure 13,
Figure 14,
Figure 15 and
Figure 16.
From the dominance analysis, it can be seen that the new Pareto front is primarily composed of operating conditions obtained following the implementation of the first and second perturbations of the wind velocity and the AC generator torque signals into the model.
In the following, the results of the variations experienced by the third
PAT family non-dominated solutions related to the reference profiles of the wind speed and external load are presented. The dominance analysis also highlights situations where a solution initially appears non-dominated after a particular perturbation. However, if another perturbation is applied, this solution might become dominated. In such cases, it is important to determine the extent of deviation the solution underwent in the performance domain. This allows understanding of how a specific hydraulic machine reacts to the perturbations taken into consideration. It is quite important to identify the best solutions, i.e., those design conditions that remain non-dominated even when subjected to perturbations in wind speed and external load.
Table 3 summarises these best solutions. The left part identifies the design conditions, the central part shows the absolute deviations from the closest non-dominated condition (in the performance domain); the performance of two operating conditions is compared by calculating the distance (norm) between them in the space defined by the relevant performance parameters (
;
;
). The right-hand side of
Table 3 displays the percentage differences in results compared to those obtained using the reference wind velocity and AC generator torque profiles (
Figure 5).
The operational stability is highlighted by the smallest relative deviations that a particular
PAT undergoes in the performance space relative to the reference signals case. Specifically, ten
PATs emerge that continue to be particularly performing, as reported in
Figure 17.
Two operating conditions for PAT3.3 exhibit zero deviation from the best design solutions identified in the four considered cases. This is due to the fact that these operating conditions belong simultaneously to the Pareto fronts obtained for all the analysed perturbations.
5. Conclusions
This research introduces a novel wind system scheme that eliminates or significantly reduces the need for battery storage by introducing a new storage system. The proposed design incorporates a reversible hydraulic machine, known as a Pump as Turbine (PAT). The PAT can operate as either a pump or a turbine by simply reversing the impeller’s rotation, enabling bidirectional fluid flow. In this plant scheme, when the wind energy exceeds the demand, the PAT, in pump mode, stores this excess (in the form of potential energy) by pumping a specific quantity of water into an elevated reservoir. However, in cases of insufficient wind speed, the PAT, operating as a turbine, contributes to overcoming the wind turbine external load. This alternative plant configuration offers numerous advantages including lower installation and operational costs, a longer lifespan for the energy storage system, and diminished waste management expenses.
This research analyses four PAT families, each comprising five PAT configurations, for a total of twenty PATs. The primary objective is to assess the robustness of the Pareto front for a wind power system integrating a pumped hydro storage system equipped with a PAT using a multi-objective approach. This analysis, conducted using MATLAB/Simulink, focuses on maximising three key performance metrics of the wind system: net water mass change in the pumped hydro storage, captured wind energy, and torque output from the PAT in turbine mode. On the other hand, the design variables considered for this analysis include PAT-specific speed (determining PAT geometry), the wind turbine to PAT speed ratio, and the reversible hydraulic machine head. The main objective of this research is to evaluate the impact of stochastic variations in wind velocity and AC generator torque on the performance of non-dominated PAT solutions, identifying those that exhibit robust performance under these conditions. This latter aspect of the analysis focuses on the third PAT family, identified as the top performer, and therefore, this family is the primary subject of this study.
The analyses conducted yield promising results, demonstrating that certain PAT configurations excel at handling perturbations in wind speed and external load, thereby enhancing operational stability. The ten best-performing PAT configurations are identified. Finally, it is reasonable to assume that these ten PATs remain viable candidates, even when subjected to various perturbations in external load and wind speed profiles. This assertion is supported by their consistently strong performance across four different case studies.
Regarding future developments, the first development concerns of the electrical generator model into the Simulink model of the wind system. Additionally, the installation and operating costs could be included in the set of variables to be optimised in order to compare this type of plant with those equipped with battery storage systems. A further deepening of the sensitivity analysis requires considering the wind speed and the AC electric generator external load as stochastic processes based on a Weibull frequency distribution. Finally, as mentioned previously, to ensure greater realism, a control system can be introduced to achieve the wind turbine maximum power point tracking.