Next Article in Journal
Defect Texts Mining of Secondary Device in Smart Substation with GloVe and Attention-Based Bidirectional LSTM
Previous Article in Journal
Redox Flow Batteries: A Literature Review Oriented to Automatic Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electro-Hydraulic Transient Regimes in Isolated Pumps Working as Turbines with Self-Excited Induction Generators

by
Filipe C. Madeira
1,
João F. P. Fernandes
2,*,
Modesto Pérez-Sánchez
3,
P. Amparo López-Jiménez
3,
Helena M. Ramos
4 and
P. J. Costa Branco
2
1
Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
2
IDMEC, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
3
Hydraulic and Environmental Engineering Department, Universitat Politècnica de València, 46022 Valencia, Spain
4
Civil Engineering, Architecture and Georesources Department, CERIS, Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Energies 2020, 13(17), 4521; https://doi.org/10.3390/en13174521
Submission received: 15 July 2020 / Revised: 3 August 2020 / Accepted: 27 August 2020 / Published: 1 September 2020
(This article belongs to the Section B: Energy and Environment)

Abstract

:
The use of pumps working as turbines (PATs) is a sustainable technical measure that contributes to the improvement of energy efficiency in water systems. However, its performance analysis in off-grid recovery systems is a complex task that must consider both hydraulic (PAT) and electrical machines (typically a self-excited induction generator-SEIG). Aside from several kinds of research that analyze the PAT-SEIG behavior under steady-state constant hydraulic and electrical conditions, this research focuses on the analysis of PAT-SEIG transient regimes, by analyzing their variation when a sudden change occurs in the hydraulic or electrical components. Analytical models were developed to represent the operation of SEIG, PAT, and the PAT-SEIG coupled system. Hydraulic and electromechanical experimental tests validated these models. An excellent fit was obtained when analytical and experimental values were compared. With these models, the impact on the operation of the PAT-SEIG system was examined when sudden change occurred in the excitation capacitances, resistive loads, or recovered head. With a sudden increase of resistive load, the hydraulic power and SEIG stator current remain almost constant. However, there is an increase of SEIG reactive power, decreasing the PAT-SEIG efficiency. Also, with a sudden increase of SEIG capacitors or PAT hydraulic head, the SEIG stator current increases once and not again, while PAT-SEIG efficiency decreases, but the induction generator can be overloaded. The development of this research is key to the advancement of future models which can analyze the coupling of micro-hydropower solutions.

1. Introduction

Micro-hydro power generation has a high potential to be used in many water sector applications where there is an excess of energy [1]. The sustainability improvement is focused on taking advantage of the hydraulic energy of water systems (i.e., pressurized or free surface) and transforming it into electric energy. This transformation needs the two main elements: a turbine and an electric generator. Both elements operate together, and therefore the analysis of the recovery system must be developed jointly for the water managers to maximize the efficiency improvement in their water systems [2]. One challenge regarding micro-hydro turbines is the cost of the turbine, which increases the initial investment, making it almost impracticable in small communities or off-grid locations [3]. One solution is the use of pumps working as turbines (PATs) since they can work in reversed mode, in a similar way as a turbine [4]. PAT is a hydraulic machine, and it requires an electrical machine to operate. Typically, commercial pumps use induction machines to convert electric power into mechanical power, and this machine is one of the most cost-effective available solutions in the market. Induction machines can also work in generator mode, connected to the grid or off-grid.
In this work, the use of PAT in rural areas where grid connection may not be possible due to physical or legal constraints will be analyzed. To work in off-grid mode, the induction generator requires an external source of reactive power. The reactive power source can be a battery coupled to the machine by an inverter or a set of capacitors connected directly to the machine. The latter solution is called a self-excited induction generator (SEIG). PAT and SEIG are coupled mechanically by a shaft operating as PAT+SEIG system (Figure 1). The analysis of this symbiosis is crucial since the global efficiency of the recovery system (PAT+SEIG) depends on each one of their efficiencies (i.e., fluid mechanical and electrical). When these efficiencies are analyzed, the maximum efficiency is located on different points. Therefore, the maximum efficiency in the PAT cannot be coincident with the maximum efficiency of the SEIG [5]. Another important aspect of these hydraulic-electromechanical systems is the investigation of transient behaviors due to sudden changes. The study of these types of influences in the system is of the utmost importance due to their effects on control stability. For example, in [6], a state feedback predictive control was applied to a hydro-turbine governing system to increase the control performance in the case of external disturbances. Also, the effect of a sudden drop in the prime mover speed of a SEIG (without the PAT) was analyzed in [7]. These effects, of sudden changes, may have an additional influence in PAT-SEIG systems during off-grid operation.
The modeling of water systems, as well as the PATs regulation, were considered in many research publications. Carravetta et al. [8] defined the different regulation types (i.e., hydraulic, electrical, and hydraulic-electrical modes). These regulation modes were compared in different water systems showing the variability of energy recovery according to the chosen strategy [9]. Different research works also focused on PAT operation and efficiency. The analysis of the tip clearance on pressure fluctuation intensity in PATs [10], cavitation performance and forces [11], or analysis of leakages [12], are also important factors to be considered. All these analyses are crucial to improve the management of the PATs. Accordingly, different optimization strategies were considered to locate the recovery systems in water systems. [13] summarized the different optimization procedures to maximize the recovered energy in the water sector. Other researches also considered the maximization of recovered energy in water systems [14], but few considered the singularity for PATs operating in off-grid mode. In this operation mode, the electrical efficiency acquires a greater significance, and the maximization of the recovered energy requires a detailed analysis due to the associated complexity.
Hence, this research focuses on analyzing the binomial hydraulic-electromechanical interaction system of Figure 1. It is composed of a pump working as a turbine (PAT) and an induction machine working as a self-excited induction generator (SEIG), which is connected to the PAT via a mechanical shaft. Inputs for the PAT+SEIG system are a differential pressure for the PAT imposed by upstream-downstream tanks, the electrical load, and the excitation capacitor bank of the SEIG. A hydraulic-electromechanical model was developed to simulate the behavior of the PAT and SEIG systems. The developed model was tuned with experimental data. The final model was validated, and it showed sufficiently accurate results for the overall system dynamics.
These analyses are representative of the system, and they allow operators to know the interaction between water and electrical machine through their different components of each part. This analysis is an important step forward to analyze the behavior of PAT+SEIG systems, which could be associated in serial or parallel in water pipe systems.

2. Self-Excited Induction Generator Model

The self-excited induction generator (SEIG) was modeled using a dq0 transformation [15] to capture the transient behavior of its electrical and magnetic quantities. Also, the change of the machine’s magnetizing inductance, L M , due to its non-linear B(H) characteristic was considered, as it highly influences the choice of the correct capacitance to self-excite the machine for a specific operating point [5] (different speeds and loads). The evolution of the magnetizing inductance as a function of the magnetization flux used here was determined experimentally in [5], for the induction machine used in the experimental results. The SEIG model was developed in Simulink, and it can be observed in Figure 2.
The dq induction machine model was developed using a stationary reference frame, aligned with phase A of the stator windings. Being all machine phases symmetric, the homopolar term will be null. The dq stator and rotor equations are shown in (1) to (4). In (1) and (2) u d s , u q s , i d s and i q s are the d and q axes stator voltages and currents, respectively, being its units V and A. Ψ d s and Ψ q s are the d and q axes stator linkage fluxes and R s is the stator phase resistance. In (3) and (4), i d r , i q r , Ψ d r and Ψ q r are the d and q axes rotor currents and linkage fluxes, respectively, R r is the rotor phase resistance (referred to the stator), and ω r is the rotor angular speed in rad/s [15].
u d s = R s i d s + d Ψ d s d t
u q s = R s i q s + d Ψ q s d t
0 = R r i d r + d Ψ d r d t ω r Ψ q r
0 = R r i q r + d Ψ q r d t + ω r Ψ d r
The relation between the stator and rotor linkage fluxes and the corresponding currents are presented in (5)–(8).
Ψ d s = L s i d s + L M i d r
Ψ q s = L s i q s + L M i q r
Ψ d r = L r i d r + L M i d s
Ψ q r = L r i q r + L M i q s
The electromagnetic torque, T e m , developed by the induction machine can be computed by (9), where n p p is the pair poles.
T e m = 3 2 n p p ( Ψ d s i q s Ψ q s i d s )
After the computation of the torque, the machine’s stator voltage and currents are converted into the three-phase quantities, a, b, and c, so that they can be compared with experimental values.

2.1. Magnetizing Inductance

As shown in [5], the SEIG model accuracy drastically increases when considering the variation of the magnetizing inductance, L M , with the magnitude of the magnetization flux, ϕ M . More specifically, for different resistive loads (RL = 120, 200 and 600 Ω), the deviation between experimental and numerical induction generator model decreased by about 37%, to an average deviation of 10.7%. The evolution of L M with the magnetization flux was obtained experimentally in [5] and was approximated by the polynomial Equation in (10). This equation shows the experimentally calculated magnetizing inductance as a function of ϕ m for the range of frequencies 20–60 Hz.
The effect of the magnetic flux in L M was included in the dq SEIG model through Equations (5)–(8). In each time instant, the model calculates the magnetic flux of the SEIG and chooses the correct value of L M . As the magnetic flux depends on the stator and rotor’s currents, the calculation of L M must be an iterative process in each time step. Fernandes et al. [5] defined the magnetizing inductance by the following equation:
L M ( ϕ M ) = 0.2192 ( ϕ M ) 3 0.8093 ( ϕ M ) 2 + 0.5531 ϕ M + 0.53
The magnetization flux can be computed using (11), where ϕ d m and ϕ q m are the d and q axes magnetization fluxes, which can be obtained using (12) and (13).
ϕ m = ϕ d m 2 + ϕ q m 2
ϕ d m = Ψ d s L s i d s
ϕ q m = Ψ q s L s i q s

2.2. Capacitor and Electrical Load Models

The SEIG required a set of capacitors (i.e., one per phase) to self-excite. The capacitors are responsible for providing the required reactive power for the induction machine to create the rotating magnetic field. However, to start the self-excitation process, the capacitors require a voltage to provide reactive power. This means that, with the induction machine rotating, no machine voltage is induced until provided with reactive power, and no reactive power is produced without applying a voltage to the capacitors. Therefore, to start the excitation process, it is required an initial voltage. This initial voltage can be stored in the capacitors, with an initial charge, or provided by the remnant voltage induced at the induction machine’s stator terminals due to asymmetries in its magnetic circuit. The residue voltage induced by the induction generator in our system was enough to start the self-excitation process. This remnant induced voltage was measured experimentally for different speeds (Figure 3). This shows the phase-to-neutral rms value of the remnant induced voltage as a function of its angular speed, N r . Each phase remnant voltage can be calculated from (14).
u r e m r m s = 0.00086 N r   ( R 2 = 0.999 )
Figure 4 shows the capacitor and load with the induction machine in parallel. The set of equations that expresses its behavior is given by (15). The capacitor bank receives the stator current of the induction machine and the current of the resistive load. It calculates the stator voltage, which is required for the induction machine and resistive blocks. In (15), u s a , b , c and i s a , b , c are the stator voltage and current of each machine’s phase, i R a , b , c are the load current of each load phase, and u r e m a , b , c are the induced remnant voltages in (14).
u s a , b , c = u C a , b , c = 1 C   i C a , b , c d t + u r e m a , b , c = 1 C ( i s a , b , c   +   i R a , b , c ) d t + u r e m a , b , c = 1 C   ( i s a , b , c + U s a , b , c R ) d t + u r e m a , b , c

3. Pump as a Turbine (PAT) Model

The pump as a turbine (PAT) model was also developed to couple with the SEIG model. The PAT model was based on the experimental curves obtained in [2,5], which were obtained for different speeds and hydraulic power. In these works, the relation between head and flow were acquired for different turbine speeds. The PAT characteristic curve can be defined by (16), according to [4,8,16], where H is the recovered head in meters water column (m w.c.), Q is the flow in l/s and A, B, and C are non-dimensional coefficients, which depend on the machine characteristics.
H   ( m   w . c . ) = A + B Q + C Q 2
In this case, the used machine was Etanorm 32-125 KSB. The calibrated characteristic curve is defined by Equation (17) [5]:
H   ( m   w . c . ) = 3.66 694.45 Q + 314560 Q 2
To obtain the new characteristic curves when the rotational speed changes, the affinity laws can be used, according to (18) [17], where α is the ratio between rotational speed N r and the nominal rotational speed ( N r e f ). In this research, the N r e f = 1050   rpm .
H = α 2 A + α B Q + C Q 2
Therefore, the used head curve to develop the hydraulic model was defined by (19), where α can be any value from 0.4 to 1.2. This range was defined considering the possibility of using the similarity laws [16]. Out of this range, the analytical model based on similarity laws, which defines the operational curve, may incur high deviations.
H = α 2 3.66 α 694.45 Q + 314560 Q 2
The PAT efficiency, η P A T , was also estimated using experimental results obtained in [2]. Figure 5 shows the efficiency surface using a linear interpolation of the experimental points for different operational points and rotational speeds.
The PAT was modeled through two main relations: its efficiency mapping is a function of PAT rotational speed and head curve (see Figure 5), and the Equation (19) to compute the flow as a function of the PAT head and its rotational speed. Considering the differential pressure (dPrh) imposed by the differential pressure tanks to the PAT terminals (PS1-PS2 in Figure 1), the PAT head can be computed using (20), where ρ is the water density (1000 kg/m3) and g is the standard acceleration due to gravity (9.81 m/s).
H = d P r h ρ g
Based on the head and SEIG-PAT speed, ω r , (resulting from the interaction between the PAT, SEIG, capacitors and electrical load), the hydraulic and mechanical torque, T h and T m e c can be computed using (21) and (22), respectively.
P h = ρ g H Q , T h = P h / ω r
P m e c = P h η P A T , T m e c = P m e c / ω r
Figure 6 shows the block diagram of the PAT model. At this point, it should be noted that the PAT model does not include the modeling of the water pipe system and the existing flow control tank, as shown in Figure 1. Therefore, the transient hydraulic dynamics related to the water circulation in the conveyance circuit and its effects on the PAT were neglected, since this analysis was developed in [17,18].
Finally, the connection between the SEIG and PAT systems are represented through the dynamic mechanical Equations (23) and (24), where T l o s s e s correspond only to the experimental values obtained for the friction and windage losses in SEIG. The overall diagram of the PAT-SEIG system is shown in Figure 7. The model inputs are the PAT head (H) and the values of the balanced resistive three-phase load and the capacitor bank.
J P A T + S E I G d ω r d t = T m e c T e m T l o s s e s
T l o s s e s = 1.05 × 10 5 N r

4. PAT-SEIG Model Validation

To validate the developed model, experimental setups were built. The error analysis between simulated and experimental data was done using different error indexes. In this research, to evaluate the fitness of the proposed model, four indicators were proposed:
  • I1. Nash–Sutcliffe index (NSI). This index is a fit indicator, which is used in temporal series. NSI value oscillates between − and 1. When values are below 0, the fit is considered poor. When the values are above 0, the model is considered good. Table 1 shows the used ranges to define the fit according to NSI values. NSI is defined by (25) where E i is the experimental value in each interval, E ¯ i is the average of the observed values and S i is the simulated value in each interval.
    N S I = 1 i = 1 N [ E i S i ] 2 i = 1 N [ E i E ¯ i ] 2
  • I2. Root Relative squared error (RRSE). It measures the error of the model by normalizing the variable. Perfect fits are defined when the RRSE value is zero. The efficiency of the simulation is better when the RRSE value is low. This index is defined by (26).
    R R S E = i = 1 N [ E i S i ] 2 i = 1 N [ E i E ¯ i ] 2
  • I3. Mean relative deviation (MRD). The index defines the significance of the error concerning variable value (27). The fit is good when MRD has values close to 0.
    M R D = 1 x | O i P i | / P i x
  • I4. Bias (BIAS). This index compares the tendency of the simulated values, determining if the simulated values are lower or higher than experimental data (28). The model overestimates if BIAS is negative. When BIAS is positive, the variable is underestimated by the model. The optimal value is zero when BIAS is analyzed.
    B I A S = i = 1 N [ O i P i ] i = 1 N [ O i ]
Table 1 shows the range values defining the goodness of the model when NSI, RRSE and BIAS are considered. These indexes were classified in very good, good, satisfactory, and unsatisfactory. The main indexes to measure the goodness of the model are the NSI, RRSE, and BIAS, according to [19,20]. The MRD value does not have a defined range for each classification. However, it provides additional information about the relative deviation between experimental and simulation results.

4.1. Self-Excited Induction Generator: d–q Model Validation

Experimental tests were done to validate the dq SEIG model under steady-state and transient conditions. For this, an experimental setup was developed in the Electrical Machines laboratory at Instituto Superior Técnico (IST), consisting of the SEIG driven by a DC motor, to impose the SEIG external torque, and a resistive electric load.
Figure 8 shows the “wye” squirrel-cage induction machine used as SEIG. Its rated values are listed in Table 2. The rated power and efficiency are P N = 0.55 kW and η N = 68%, being its nominal slip, s N , equal to 9%. Figure 9 shows the experimental setup. It is composed by two components. An electromechanical component associated with a DC motor coupled mechanically to the SEIG to act as an external torque source to the generator (Figure 9a) and an electric circuit component made by the variable capacitor and resistor banks, being both connected in parallel to the stator coils of the induction generator (Figure 9b). A general overview of the two experimental setup components is shown in Figure 9c. It shows the electric grid connected to an auto-transformer, which feeds the armature of the DC separated-excitation motor by a diode bridge rectifier. The DC machine is rigidly connected to the induction generator that feeds a capacitor and variable resistive bank.
The first series of tests were done to validate the transient dynamics of the induction generator during its initial self-excitation. Without the capacitor bank and resistive load, the generator is accelerated to a constant reference speed. During this operating condition and without the capacitor bank connected, the generator stator voltage is the remnant one, as indicated in Figure 3, for a speed of N r = 830 rpm, u r e m r m s = 0.71 V. Next, the capacitor bank is connected to begin the self-excitation process. Figure 10 shows the experimental (blue line) and simulated (red line) stator voltage results of the SEIG excitation process using two different capacitor bank values: (a) C = 50 μF and (b) C = 80 μF. Simulation results were obtained using the SEIG model developed in Section 2 (Figure 2). To replicate the laboratorial conditions, the SEIG model was simulated with an external mechanical load applied to its axis (provided from a DC motor, as seen in Figure 9). This was achieved using a similar model to the one shown in Figure 7, where the PAT was replaced by a DC motor torque-speed characteristic ( T d c = k ( U d c k ω r ) / R a , with k and R a equal to 1.05 and 1.6 Ω, respectively). Without capacitors to excite the induction generator, the DC motor torque was adjusted to meet the initial speed conditions, N r , of each experimental test. After reaching steady state, the capacitors values are added to starting the self-excitation process of the generator. Figure 10 shows the stator voltage in one phase of the SEIG during the excitation process for two different capacitor values: C = 50 μF in Figure 10a and C = 80 μF in Figure 10b.
After the establishment of the SEIG stator voltage, power losses occur inside the machine. These are composed mainly by the stator windings and rotor bars losses, P J , plus the magnetic core losses, P c o r e . The windings and rotor bars losses are due to the current injection from the capacitors to the SEIG, being the rotor bars losses caused by the induced rotor currents derived from the electromotive forces produced by the rotating magnetic field which is originated from the currents injected by the capacitor bank. Magnetic core losses are due to the rotating magnetic field created inside the SEIG. Therefore, from the mechanical point of view, the SEIG will demand a higher mechanical torque related to these losses, and thus, after the excitation, the SEIG speed decreases. The initial speed of the SEIG was 830 rpm (before excitation), and the final ones were 750 and 597 rpm (after excitation), respectively.
These error values were analyzed for these figures (Figure 10a,b) according to NSI, RRSE, MRD, and BIAS. These error indexes were described in Section 4 of this research. If stator voltage is analyzed when the capacitance was 50 μF, a very good fit was reached, according to Table 1. In this case, NSI, RRSE and BIAS were 0.955, 0.212 and 0.07, respectively. When the stator voltage was analyzed for a capacitance of 80 μF, the best fit was also very good, with NSI 0.93, RRSE 0.264, and BIAS 0.0357. Therefore, the developed model to simulate the self-excitation process can be considered correct. The deviations between experimental and simulation results are less than 0.044 (4.4% in absolute) for the final steady-state values (Table 3).
Once the validation of the self-excitation dynamic process was developed, the SEIG response to a change of electrical load needs to be analyzed. The same DC motor-SEIG model was used, now to simulate the behavior of the SEIG after a connecting of resistive load. The simulation was done in two steps: first, the SEIG was excited without electrical load and, after reaching a steady-state condition, the resistive load is connected. Figure 11 and Figure 12 show the experimental and simulation results of the SEIG: (a) stator voltage, (b) stator current, (c) active power, (d) reactive power, and (e) magnetizing flux; after connecting an RL of 600 Ω and 300 Ω, respectively. These resistive loads mimic an electric demand to the SEIG from a consumer. Since in the experimental setup, a DC motor provides the SEIG mechanical torque, its speed decreases with the increase of electric load. Notice that this is similar behavior to the PAT. When decreasing the resistive load, the electric current increases and, if the stator voltage remains (or has a very small change), the requested active power P a = 3 U s 2 R L would increase. However, Figure 11 shows that this is not true. When the requested active power increases, the speed of the DC motor decreases and thus also the SEIG electrical frequency and stator voltage. Decreasing the stator voltage will also decrease the stator current, Figure 11a,b, decreasing the active power, Figure 11c, and increasing the machine’s speed. Also, a reduction of electrical frequency will influence the reactive power, as shown in Figure 11d, provided by the capacitor to the SEIG, thus changing the SEIG operation point, per Figure 11e. In conclusion, this process will be iterative until a new steady state is obtained. The SEIG starts without electric load, with a capacitor bank of C = 35 μF, and the current at its rated value, I = 1.6 Arms. After the inclusion of the resistive load, R L = 600 Ω, the rms stator voltage and current drop, from 183 Vrms to 141 Vrms and from 1.6 Arms to 1.05 Arms, Figure 11a,b, respectively. This is due to the reduction of speed, N, and electrical frequency, f s , resulting into new active and reactive powers, shown in Figure 11c,d, and a new magnetizing flux inside the machine, E / f , shown in Figure 11e. As shown in Table 4, MRDs between the simulation and experimental results are less than 0.088 (8.8%). The highest deviation of 0.088 occurs for low magnetization points of the SEIG, E / f = 2.2. This is related to the magnetizing inductance, L M , obtained experimentally. For lower values of magnetic flux, the stator leakage cannot be neglected and, thus, the associated error to the determination of L M increases.
These errors were estimated to measure the goodness of fit. When the errors were calculated to the stator voltage (Figure 11a), the values of NSI, RRSE and BIAS were 0.629, 0.609 and −0.0434, respectively, reaching very good values for NSI and BIAS, and good for RRSE. Similar values were got when the stator current was analyzed between experimental and simulated values. Very goodness fit was reached for three indicators (NSI [0.687], RRSE [0.386] and BIAS [0.0105]). Excellent values were reached when active and reactive power were analyzed. The NSI, RRSE, and BIAS were 0.964, 0.189 and 0.055, respectively, for the active power. In contrast, the values were 0.933 (NSI), 0.258 (RRSE) and 0.0281 (BIAS) for the reactive power. Finally, when the magnetizing flux was studied, the values of NSI, RRSE, and BIAS were 0.784, 0.138, and −0.0317, respectively. Therefore, a very good fit was reached in this model when the SEIG was compared, connecting a resistive load equal to 600 Ω.

4.2. PAT Model Validation

A second set of experimental tests was performed to validate the PAT model represented by the PAT model in Figure 6. Figure 13a,b shows a scheme of the hydraulic circuit used with the PAT installed. The hydraulic circuit is composed of a pressure tank to constant guaranty pressure applied to the radial PAT and an Etanorm 32-125 KSB 4.8 PAT. This machine (PAT1) had a nominal rotational speed of 1020 rpm and a rated flow of 4.2 L/s, a flow control tank to eliminate possible turbulences in the hydraulic circuit, and a recirculating pump to close the hydraulic circuit. The specific speed was 51 rpm (m, kW), and the best efficiency point (BEP) was defined in the operation point 3.36 L/s and 4 m w.c. Figure 14 shows the experimental setup assembled at the laboratory (IST). Figure 14a highlights the PAT (1), the induction generator (2), the recirculating pump (3), the control flow tanks (4-5), and the pressure transducer (6). In Figure 14b, one highlights the measurement and the oscilloscope (7), the resistive load bank (8), the capacitor bank (9), the pressure tank (10), and the flowmeter (11).
A set of experimental tests were done to validate the PAT model in steady-state and transient conditions. Simulation results were now obtained using the PAT-SEIG model in Figure 7. Models’ inputs are the PAT head and the values of the resistive load and capacitors. All simulations start with a non-excited SEIG to obtain the PAT steady-state for a set of initial conditions. After reaching the steady-state hydraulic conditions, the capacitors are added to the SEIG to start with its excitation process, with and without a resistive load. Following the new steady-state conditions of the PAT-SEIG system (with and without resistive load), a sudden change of resistive load is applied to the SEIG, and the new transient behavior can be analyzed.
Starting without the capacitor bank and at no-load, a differential head of 5.82 m w.c. was applied to the PAT’s terminals, regulating the air vessel (Figure 14b element 10) to guarantee a uniform head value. After the establishment of the PAT steady-state conditions, the excitation capacitor bank was showing a value of 17.4 μF connected to the no-load SEIG to start its excitation process, as shown in Figure 15a,b. It shows the evolution of the no-load SEIG stator voltage and the electric frequency (experimental and simulation) until the magnetic saturation of the induction generator limited it, stabilizing around 250 V and 55 Hz. To verify the affection of the electric load in the excitation process, the previous starting process was repeated. In this case, an electrical load equal to 300 Ω was connected. The capacitor bank value was varied from C = 17.4 μF to C = 34.7 μF to guaranty the excitation of the induction generator. The results of this experimental test are shown in Figure 15c,d. When these figures are observed, the time to excite the SEIG increases when it starts with an electrical load and presents an oscillatory transient. Therefore, the best methodology is to start first the excitation without electrical load and to connect the electrical load only after the SEIG voltage is stabilized. If a comparison is developed, the analytical approximation is really good. It is confirmed when the errors are considered (Table 5). These show a very good fit for all analyzed errors applied to all variables. Therefore, the model correctly reproduced the behavior of the PAT-SEIG model.
Once the validation of the excitation process was developed, two additional tests were performed to verify the accuracy of the PAT-SEIG coupled models under transient conditions. After exciting the SEIG, a sudden change of electrical load or capacitor bank was applied to the PAT-SEIG system. Figure 16a,b shows the stator voltage and electric frequency transient behaviors before and after a reduction of the electric load. The PAT-SEIG system started with R L = 600 Ω, C = 17.4 μF and H = 7.5 m w.c. When the establishment of the SEIG steady-state was around 170 V of amplitude (120.2 Vrms), the electric load was decreased to R L = 300 Ω at 0.1s. With a lower resistance, the SEIG stator voltage drops, performing an oscillatory transient to the new steady-state regime around 130 V of amplitude (92 Vrms).
In Figure 16c,d, the load remained at R L = 300 Ω but an increase of the capacitor bank value was applied to the PAT-SEIG system at around 0.22 s. The PAT-SEIG started considering a C value of 17.4 μF, and it was increased to 34.7 μF. An overvoltage is seen in the SEIG stator voltage from 120 V to around 190 V, and a reduction of the electric frequency occurs from 58 Hz to 40 Hz.
In all tests, the simulation results follow the experimental ones. To measure the errors between them and to calibrate the model, the errors were defined as well as the goodness of fit. Table 6 shows a very good fit for all analyzed parameters. The error was very good in all analyses for different parameters, except for the stator voltage. For this parameter, the values were good and satisfactory when the NSI and RRSE were determined. The obtained values showed that the analytical model could be used in other hydraulic and electrical conditions, such as when these recovery systems operate in serial or parallel. MRD also showed satisfactory values since MRD values are under 0.1 (10%).

5. Impact of Electric and Hydraulic Perturbations in the PAT-SEIG Stability

The PAT-SEIG system is a complex hydraulic-electromechanical non-linear system, which has a high number of variables. These variables change over time rapidly, and this variation depends on the hydraulic and electric constraints (i.e., flow, head, electric load and capacitor value). The transient dynamics of the whole system results from the interaction behaviors of both PAT and SEIG subsystems. This means that a small sudden change in an electrical variable may influence the hydraulic one and vice-versa. In this section, a sudden change of electrical or hydraulic parameters is analyzed to show how this affects the system operation after reaching a steady-state condition. The tested scenarios consider sudden variations of (1) the value of the capacitor bank, (2) the resistive load, and (3) the hydraulic pressure.
Due to the limitation of the Hydraulic lab conditions in terms of the air-vessel and PAT hydraulic power, the SEIG was not able to reach its rated power within the tested range of speeds. Therefore, for simulation purposes, the PAT curve ( H   vs .   Q ) was reallocated to higher head values, as shown in Figure 17. This allows higher hydraulic power within the same range of speeds. The new proposal PAT was defined by the characteristic curve (29).
H   ( m   w . c . ) = 10.99 694.45 Q + 314560 Q 2
Figure 17 shows the hydraulic curves (H vs. Q and Ph vs. Q) for the original PAT (PAT1 in blue) and the new one (PAT2 in red). PAT2 is now able to produce enough hydraulic power to set the SEIG rated conditions (rated power of 550 W). The specific speed of PAT2 was 20 rpm (m,kW), and the best efficiency point (BEP) was defined in the operation point 3.48 L/s and 12.38 m w.c.

5.1. Variation of Excitation Capacitance

The first set of tests applied a sudden change in the excitation capacitor values, with the PAT-SEIG system already excited and under steady-state conditions. The initial steady-state conditions were reached for a differential PAT head, H , of 21.5 m w.c., an initial capacitor value, Ci, of 35 μF and an electric load of R L = 200 Ω. Next, the capacitors’ value was changed, and the new steady-state electromechanical and hydraulic conditions were recorded. Table 7 shows the final steady-state hydraulic and electric parameters after a sudden change of the capacitor’s value.
In this table, C is the applied variation of the capacitors in (%), P h is the resultant hydraulic power, N the PAT-SEIG speed, U s r m s and I s r m s are the SEIG stator voltage and current rms values, P a and Q r are the SEIG active and reactive power, and η s y s is the PAT-SEIG overall achieved efficiency. Note that the negative value in Q indicates that the reactive power is flowing from the capacitors to the SEIG.
Figure 18 is also possible to visualize the variation of the hydraulic, active and reactive powers and the system speed and efficiency. It is possible to observe that decreasing the capacitors’ value increases the system speed and efficiency. This follows the results obtained in [2,5], where an inverse relationship between the speed and capacitors’ values was obtained. The active power (red line in Figure 18a) provided from the SEIG shows a slight ± 15 % change of its value for the capacitor’s variation between −50% and +50%. In contrast, the reactive one (blue line) follows the variation of the capacitors. Due to mainly reactive power, the SEIG magnetic and windings losses decrease with the reduction of capacitors (red line in Figure 18b), thus requesting lower values of hydraulic power from the PAT for lower values of C. This shows that it is possible to control the PAT hydraulic power from changing the capacitor values in the SEIG.
Another important aspect is the behavior of the SEIG stator current with the variation of capacitors. Increasing the capacitors’ value increased the SEIG current. This can lead to an overload state of the SEIG that, if not considered, may significantly decrease the generator’s lifetime and/or damages it due to overheating of its electric insulation.

5.2. Variation of Resistive Load

The same study was done regarding a sudden change of the resistive load. Equal initial steady-state conditions were considered using a differential PAT head, H , of 21.5 m w.c., a capacitor value, C, of 35 μF and an initial electric load of R L i = 200 Ω. After reaching steady-state, a sudden change in the electrical resistance was applied.
Results in Table 8 and Figure 19 show that the change of SEIG’s resistive load highly influences the SEIG stator voltage and reactive power. Being the capacitor bank in parallel to the SEIG stator terminals, both share the same voltage. Therefore, a decrease of SEIG’s terminal voltage implies a decrease of reactive power stored at the capacitor bank. Despite the reactive evolution, the active power remains almost unchanged with the electric resistance, presenting a slight increase with the decrease of the resistance. For R L = −30%, the SEIG starts losing its excitation, getting near to its non-excited zone. The speed and PAT-SEIG system efficiency have similar behavior with the variation of the electrical resistance: they both decrease with the increase of the electric resistance. This is in agreement with the behaviors of the hydraulic, reactive and active power. Despite the active power remaining almost constant, with the increase of the reactive power, the magnetic power losses inside the SEIG increases, thus resulting in a decrease of the SEIG efficiency.
Another important aspect is the evolution of the current. When decreasing the electric resistance, the current presents a slow reduction. However, with the increase of the resistance, the current presents a higher increase. This must be taken into consideration, as an increase of the SEIG current to values higher than its nominal one, which may damage the SEIG in a medium or long term.

5.3. PAT Head Variation

The PAT head variation ( H ) is an important issue, as it is common to occur in industrial applications due to external influences. This case analyzed the impact of a sudden increase or decrease of the PAT head in the overall PAT-SEIG operation. The same initial steady-state conditions were considered with an initial differential PAT head, H i , of 21.5 m w.c., a capacitor value, C, of 35 μF and an electric load of R L = 200 Ω.
Results presented in Figure 20 and Table 9 shows there is a linear behavior of the hydraulic, active and reactive power and rotational speed with the variation of PAT head. Regarding the system efficiency, results indicated that only in an extreme case characterized by a significant decrease of head pressure (−50%) would the PAT-SEIG suffer a meaningful reduction of its efficiency. Another important effect is the increase of SEIG current with the PAT head, which may damage the SEIG if the current exceeds its rated value.

6. Conclusions

The study of the improvement of sustainable indicators in water systems has been an important issue in this last decade, especially in the context of sustainable development goals for 2030. The improvement of the efficiency in water systems can be leveraged by using PATs, which have been mainly studied in the last years under the hydraulic point of view. However, the PATs analysis cannot only be focused on hydraulic behavior, as PAT requires an electrical machine to operate. In the particular case of off-grid (absence of electrical grid) applications, using self-excited induction generators, the hydraulic-electric behavior becomes more complex.
This research focuses on delving into the study of the hydraulic-electrical transient stability interaction when possible perturbations occur in an off-grid PAT-SEIG system. Analytical models were developed to simulate the PAT and the SEIG behaviors, both individually validated through experimental tests. Then the coupled PAT-SEIG model also validated experimentally in IST hydraulic laboratory, was tested for different electrical loads. To analyze the accuracy of the models, error measures were considered, resulting in NSI values between 0.75 and 0.98 and BIAS values below 0.10.
The developed models can estimate the variation in SEIG voltage, current and frequency, as well as in the PAT flow and recovered head when a sudden change in the resistive load, capacitance and PAT head occurs. With a sudden increase of resistive load, the hydraulic power and SEIG stator current remain almost constant. However, there is an increase of the SEIG reactive power and a decrease in the PAT-SEIG efficiency. With a sudden increase of the SEIG capacitors or PAT head, the SEIG stator current increases and PAT-SEIG efficiency decreases. The increase of the SEIG current may be critical, as it may contribute to the machine overload.
The possibility to estimate these variations will allow operators to define new operating models, in which new recovery systems can be coupled between them. This analytical coupling will measure the injected perturbations in the water grid or electrical batteries when the electric and the hydraulic transients occur. These developed models are the basis for the real integration of the PATs in the water sector, where they can correctly operate under hydraulic and/or electric regulation conditions. The use of these models can be extrapolated to other off-grid electro-hydraulic systems, considering both electrical and hydraulic parameters of the used machines.

Author Contributions

Conceptualization, methodology and software: F.C.M. and J.F.P.F.; validation and formal analysis: J.F.P.F., M.P.-S., and P.J.C.B.; writing—original draft preparation, writing—review and editing, J.F.P.F., M.P.-S., P.J.C.B., and H.M.R.; supervision and lab concept, J.F.P.F., M.P.-S., P.J.C.B., and H.M.R.; final review H.M.R. and P.A.L.-J. All authors have read and agreed to the published version of the manuscript.

Funding

This research received some support from the project REDAWN (Reducing Energy Dependency in Atlantic Area Water Networks) EAPA_198/2016 from INTERREG ATLANTIC AREA.

Acknowledgments

This work was supported by FCT—Foundation for Science and Technology, I.P., through IDMEC, under LAETA, project UIDB/50022/2020. This research is supported by the “Program to support the academic career of the faculty of the Universitat Politècnica de València 2019/2020 in the project ‘A STEP AHEAD IN SUSTAINABILITY OF WATER SYSTEMS FOR THE ENERGY TRANSITION IN COMMUNITIES‘ of the Modesto Pérez-Sánchez. The authors wish to thank the project REDAWN (Reducing Energy Dependency in Atlantic Area Water Networks) EAPA_198/2016 from INTERREG ATLANTIC AREA Fluids PROGRAMME 2014–2020 and CERIS, the Hydraulic Laboratory, for the support in the conceptual developments and the experiments on PATs.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Postacchini, M.; Darvini, G.; Finizio, F.; Pelagalli, L.; Soldini, L.; di Giuseppe, E. Hydropower generation through pump as turbine: Experimental study and potential application to small-scale WDN. Water 2020, 12, 958. [Google Scholar] [CrossRef] [Green Version]
  2. Capelo, B.; Pérez-Sánchez, M.; Fernandes, J.F.P.; Ramos, H.M.; López-Jiménez, P.A.; Branco, P.J.C. Electrical behaviour of the pump working as turbine in off grid operation. Appl. Energy 2017, 208, 302–311. [Google Scholar] [CrossRef]
  3. Fecarotta, O.; Aricò, C.; Carravetta, A.; Martino, R.; Ramos, H.M. Hydropower potential in water distribution networks: Pressure control by PATs. Water Resour. Manag. 2014, 29, 699–714. [Google Scholar] [CrossRef]
  4. Ramos, H.; Borga, A. Pumps as turbines: An unconventional solution to energy production. Urban Water 1999, 1, 261–263. [Google Scholar] [CrossRef]
  5. Fernandes, J.F.P.; Pérez-Sánchez, M.; Silva, F.F.; López-Jiménez, P.A.; Ramos, H.M.; Branco, P.J.C. Optimal energy efficiency of isolated PAT systems by SEIG excitation tuning. Energy Convers. Manag. 2019, 183, 391–405. [Google Scholar] [CrossRef]
  6. Yi, Y.; Zhang, Z.; Chen, D.; Zhou, R.; Patelli, E.; Tolo, S. State feedback predictive control for nonlinear hydro-turbine governing system. J. Vib. Control 2018, 24, 4945–4959. [Google Scholar] [CrossRef]
  7. Khan, M.F.; Khan, M.R. Analysis of voltage build-up and speed disturbance ride-through capability of a self-excited induction generator for renewable energy application. Int. J. Power Energy Convers. 2016, 6, 157–177. [Google Scholar] [CrossRef]
  8. Carravetta, A.; del Giudice, G.; Fecarotta, O.; Ramos, H. PAT design strategy for energy recovery in water distribution networks by electrical regulation. Energies 2013, 6, 411–424. [Google Scholar] [CrossRef] [Green Version]
  9. Morani, M.C.; Carravetta, A.; del Giudice, G.; McNabola, A.; Fecarotta, O. A comparison of energy recovery by pATs against direct variable speed pumping in water distribution networks. Fluids 2018, 3, 41. [Google Scholar] [CrossRef] [Green Version]
  10. Liu, Y.; Tan, L. Tip clearance on pressure fluctuation intensity and vortex characteristic of a mixed flow pump as turbine at pump mode. Renew. Energy 2018, 129, 606–615. [Google Scholar] [CrossRef]
  11. Liu, Y.; Tan, L. Symmetrical and unsymmetrical tip clearances on cavitation performance and radial force of a mixed flow pump as turbine at pump mode. Renew. Energy 2018, 127, 368–376. [Google Scholar]
  12. Han, Y.; Tan, L. Dynamic mode decomposition and reconstruction of tip leakage vortex in a mixed flow pump as turbine at pump mode. Renew. Energy 2020, 155, 725–734. [Google Scholar] [CrossRef]
  13. Pérez-Sánchez, M.; Sánchez-Romero, F.J.; López-Jiménez, P.A.; Ramos, H.M. PATs selection towards sustainability in irrigation networks: Simulated annealing as a water management tool. Renew. Energy 2018, 116, 234–249. [Google Scholar] [CrossRef]
  14. García, I.F.; Nabola, A.M. Maximizing hydropower generation in gravity water distribution networks: Determining the optimal location and number of pumps as turbines. J. Water Resour. Plan. Manag. 2020, 146, 04019066. [Google Scholar] [CrossRef]
  15. Simões, M.G.; Farret, F.A. Modeling and Analysis with Induction Generators, 3rd ed.; CRC Press: Boca Raton, FL, USA, 2014; ISBN 9781482244670. [Google Scholar]
  16. Mataix, C. Turbomáquinas Hidráulicas; Universidad Pontificia Comillas: Madrid, Spain, 2009. [Google Scholar]
  17. Pérez-Sánchez, M.; López-Jiménez, P.A.; Ramos, H.M. PATs operating in water networks under unsteady flow conditions: Control valve manoeuvre and overspeed effect. Water 2018, 10, 529. [Google Scholar] [CrossRef] [Green Version]
  18. Carravetta, A.; Houreh, S.D.; Ramos, H.M. Pumps as Turbines, Fundamentals and Applications; Springer: Cham, Switzerland, 2018. [Google Scholar] [CrossRef]
  19. Pérez-Sánchez, M.; Sánchez-Romero, F.J.; Ramos, H.M.; López-Jiménez, P.A. Calibrating a flow model in an irrigation network: Case study in Alicante, Spain. Span. J. Agric. Res. 2017, 15, e1202. [Google Scholar] [CrossRef]
  20. Moriasi, B.; Arnold, J.; van Liew, M.; Binger, R.; Harmel, R.; Veith, T. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
Figure 1. PAT+SEIG scheme. Hydraulic and electrical circuit.
Figure 1. PAT+SEIG scheme. Hydraulic and electrical circuit.
Energies 13 04521 g001
Figure 2. SEIG model block diagram.
Figure 2. SEIG model block diagram.
Energies 13 04521 g002
Figure 3. Induced remnant voltage of the induction machine for different rotational speeds.
Figure 3. Induced remnant voltage of the induction machine for different rotational speeds.
Energies 13 04521 g003
Figure 4. Connection between induction machine, capacitors and resistive load.
Figure 4. Connection between induction machine, capacitors and resistive load.
Energies 13 04521 g004
Figure 5. PAT efficiency experimental curve.
Figure 5. PAT efficiency experimental curve.
Energies 13 04521 g005
Figure 6. PAT model.
Figure 6. PAT model.
Energies 13 04521 g006
Figure 7. SEIG-PAT model.
Figure 7. SEIG-PAT model.
Energies 13 04521 g007
Figure 8. Squirrel-cage induction machine used as SEIG.
Figure 8. Squirrel-cage induction machine used as SEIG.
Energies 13 04521 g008
Figure 9. Experimental setup to test the SEIG, isolated from the PAT: (a) DC motor coupled mechanically to the SEIG, (b) variable capacitor and resistor banks both connected in parallel to the stator coils of the induction generator, (c) scheme of two experimental setup components.
Figure 9. Experimental setup to test the SEIG, isolated from the PAT: (a) DC motor coupled mechanically to the SEIG, (b) variable capacitor and resistor banks both connected in parallel to the stator coils of the induction generator, (c) scheme of two experimental setup components.
Energies 13 04521 g009
Figure 10. SEIG stator voltage during the excitation process with: (a) C = 50 μF and (b) C = 80 μF. In blue, the experimental results. In red, the simulation ones.
Figure 10. SEIG stator voltage during the excitation process with: (a) C = 50 μF and (b) C = 80 μF. In blue, the experimental results. In red, the simulation ones.
Energies 13 04521 g010
Figure 11. SEIG stator voltage when connecting a resistive load of R L = 600 Ω with the SEIG excited with a capacitor bank having C = 35 μF: (a) stator voltage, (b) stator current, (c) active power, (d) reactive power and (e) magnetizing flux.
Figure 11. SEIG stator voltage when connecting a resistive load of R L = 600 Ω with the SEIG excited with a capacitor bank having C = 35 μF: (a) stator voltage, (b) stator current, (c) active power, (d) reactive power and (e) magnetizing flux.
Energies 13 04521 g011
Figure 12. SEIG stator voltage when connecting a resistive load of R L = 300 Ω with the SEIG excited with a capacitor bank having C = 35 μF.
Figure 12. SEIG stator voltage when connecting a resistive load of R L = 300 Ω with the SEIG excited with a capacitor bank having C = 35 μF.
Energies 13 04521 g012
Figure 13. Experimental hydraulic circuit: (a) circuit schematic and (b) PAT used (Etanorm 32-125 KSB 4.8).
Figure 13. Experimental hydraulic circuit: (a) circuit schematic and (b) PAT used (Etanorm 32-125 KSB 4.8).
Energies 13 04521 g013
Figure 14. Experimental hydraulic circuit developed for validation of the PAT model: in (a) the PAT-SEIG, control flow tanks and recirculating pump and, in (b) the pressure tank and the resistive load and capacitors.
Figure 14. Experimental hydraulic circuit developed for validation of the PAT model: in (a) the PAT-SEIG, control flow tanks and recirculating pump and, in (b) the pressure tank and the resistive load and capacitors.
Energies 13 04521 g014
Figure 15. PAT-SEIG stator voltage and frequency during excitation with: C = 17.4 μF, R L = Ω and H = 5.82 m w.c: (a) stator voltage and (b) electric frequency; and with C = 34.7 μF, R L = 300 Ω and H = 7.5 m w.c.: (c) stator voltage and (d) electric frequency.
Figure 15. PAT-SEIG stator voltage and frequency during excitation with: C = 17.4 μF, R L = Ω and H = 5.82 m w.c: (a) stator voltage and (b) electric frequency; and with C = 34.7 μF, R L = 300 Ω and H = 7.5 m w.c.: (c) stator voltage and (d) electric frequency.
Energies 13 04521 g015
Figure 16. PAT-SEIG stator voltage and electric frequency after a sudden change: (a,b) change of electric load from R L = 600 Ω to R L = 300 Ω, with C = 17.4 μF and H = 7.5 m w.c. and (c,d) change of capacitor from C = 17.4 μF to C = 34.7 μF with R L = 300 Ω and H = 7.5 m.w.c.
Figure 16. PAT-SEIG stator voltage and electric frequency after a sudden change: (a,b) change of electric load from R L = 600 Ω to R L = 300 Ω, with C = 17.4 μF and H = 7.5 m w.c. and (c,d) change of capacitor from C = 17.4 μF to C = 34.7 μF with R L = 300 Ω and H = 7.5 m.w.c.
Energies 13 04521 g016
Figure 17. Hydraulic characteristic curves of the original PAT (PAT1, blue color) and the new one (PAT2, red color), for N = Nref: (a) head curve and (b) hydraulic power curve.
Figure 17. Hydraulic characteristic curves of the original PAT (PAT1, blue color) and the new one (PAT2, red color), for N = Nref: (a) head curve and (b) hydraulic power curve.
Energies 13 04521 g017
Figure 18. Variation of hydraulic and electromechanical quantities after a change in the capacitor, C: (a) variation of hydraulic (black), active (red) and reactive (blue) powers, and (b) variation of speed (black) and system (red) efficiency.
Figure 18. Variation of hydraulic and electromechanical quantities after a change in the capacitor, C: (a) variation of hydraulic (black), active (red) and reactive (blue) powers, and (b) variation of speed (black) and system (red) efficiency.
Energies 13 04521 g018
Figure 19. Variation of hydraulic and electromechanical parameters after a change in the resistive load, R L : (a) variation of hydraulic (black), active (red) and reactive (blue) powers, and (b) variation of speed (black) and system (red) efficiency.
Figure 19. Variation of hydraulic and electromechanical parameters after a change in the resistive load, R L : (a) variation of hydraulic (black), active (red) and reactive (blue) powers, and (b) variation of speed (black) and system (red) efficiency.
Energies 13 04521 g019
Figure 20. Variation of hydraulic and electromechanical parameters after a change in the PAT head, H : (a) variation of hydraulic (black), active (red) and reactive (blue) powers, and (b) variation of speed (black) and system (red) efficiency.
Figure 20. Variation of hydraulic and electromechanical parameters after a change in the PAT head, H : (a) variation of hydraulic (black), active (red) and reactive (blue) powers, and (b) variation of speed (black) and system (red) efficiency.
Energies 13 04521 g020
Table 1. Classification of best fit, according to [19,20].
Table 1. Classification of best fit, according to [19,20].
Goodness FitNSIRRSEBIAS
Very GoodNSI > 0.60.00 ≤ RRSE ≤ 0.50 BIAS < ± 0.10
Good0.40 < NSI ≤ 0.600.50 < RRSE ≤ 0.60 ± 0.10 BIAS < ± 0.15
Satisfactory0.20 < NSI ≤ 0.400.60 < RRSE ≤ 0.70 ± 0.15 BIAS < ± 0.25
UnsatisfactoryNSI < 0.20RRSE > 0.70 BIAS > ± 0.25
Table 2. Nameplate data of the squirrel-cage induction machine used as SEIG.
Table 2. Nameplate data of the squirrel-cage induction machine used as SEIG.
Frequency50 Hz
Voltage400 V
Current1.6 A
Output Power0.55 kW
Power factor0.73
Speed910 rpm
Table 3. Experimental and model steady-state values after the excitation process.
Table 3. Experimental and model steady-state values after the excitation process.
ExperimentalModelMRD
C = 50 μFN (rpm)750758+0.010 (1.0%)
f s (Hz)35.235.0−0.006 (0.6%)
U s (Vrms)144145+0.007 (0.7%)
C = 80 μFN (rpm)597603+0.010 (1.0%)
f s (Hz)27.627.2−0.015 (1.5%)
U s (Vrms)113108−0.044 (4.4%)
Table 4. Comparison between steady-state experimental and model’s results after the application of resistive loads.
Table 4. Comparison between steady-state experimental and model’s results after the application of resistive loads.
ExperimentalModelMRD
InitialFinalInitialFinalInitialFinal
R L = 600 ΩN (rpm)839834842835+0.004 (+0.4%)+0.001 (+0.1%)
f s (Hz)41.040.040.939.3−0.000 (−0.0%) −0.018 (−1.8%)
U s (Vrms)183141191141.9−0.044 (−4.4%)−0.006 (−0.6%)
I s (Arms)1.61.051.531.12−0.023 (−2.3%)+0.023 (+2.3%)
E / f (V Hz 1 )4.463.534.673.61−0.005 (−0.5%)+0.003 (+0.3%)
R L = 300 ΩN (rpm)848843849851+0.001 (0.1%)+0.010 (+1.0%)
f s (Hz)41.240.341.540.0+0.007 (+0.7%)−0.007 (−0.7%)
U s (Vrms)18190184.684.8+0.002 (+0.2%)+0.058 (+5.8%)
I s (Arms)1.60.81.620.87+0.013 (+1.3%)+0.088 (+8.8%)
E / f (V Hz 1 )4.42.24.492.22+0.020 (+2.0%)+0.009 (+0.9%)
Table 5. Error analysis for the compared variables.
Table 5. Error analysis for the compared variables.
IndexesStator Voltage
(Figure 15a)
Frequency
(Figure 15b)
Stator Voltage
(Figure 15c)
Frequency
(Figure 15d)
NSI0.960 (VG)0.797 (VG)0.802 (VG)0.841 (VG)
RRSE0.200 (VG)0.451 (VG)0.441 (VG)0.399 (VG)
BIAS0.039 (VG)−0.071 (VG)0.095 (VG)−0.075 (VG)
MRD−0.0244 (−2.44%)0.062 (6.2%)0.010 (1.0%)0.0664 (6.64%)
(VG.—Very Good; G.—Good; S.—Satisfactory; P.—Poor).
Table 6. Error indexes for the compared variables.
Table 6. Error indexes for the compared variables.
IndexesStator Voltage
(Figure 16a)
Frequency
(Figure 16b)
Stator Voltage
(Figure 16c)
Frequency
(Figure 16d)
NSI0.569 (G)0.787 (VG)0.761 (VG)0.793 (VG)
RRSE0.648 (S)0.462 (VG)0.429 (VG)0.455 (VG)
BIAS0.063 (VG)−0.041 (VG)−0.034 (VG)−0.037 (VG)
MRD0.03600.0269−0.0540.281
(VG.—Very Good; G.—Good; S.—Satisfactory; P.—Poor).
Table 7. Steady-state hydraulic and electromechanical quantities after a sudden change in the capacitors’ value.
Table 7. Steady-state hydraulic and electromechanical quantities after a sudden change in the capacitors’ value.
C. P h   ( W ) N (rpm) U s r m s   ( V ) I s r m s   ( A ) P a   ( W ) Q r   ( Var ) η s y s
−50%10051365139.71.16331−39233.0%
−40%12231252152.61.34392−51532.0%
−30%13331173157.41.47416−59231.1%
−20%14011114158.81.58427−65830.5%
−10%14601057155.51.60403−66127.6%
0%15011010150.21.65384−66425.6%
+10%1512997157.01.83417−78327.6%
+20%1535969155.21.90412−81026.9%
+30%1553944154.41.95408−83326.3%
+40%1568922152.12.01403−85125.7%
+50%1581903149.42.07387−86624.4%
Table 8. Steady-state hydraulic and electromechanical parameters after a sudden change in the resistance value.
Table 8. Steady-state hydraulic and electromechanical parameters after a sudden change in the resistance value.
R L P h   ( W ) N (rpm) U s r m s   ( V ) I s r m s   ( A ) P a   ( W ) Q r   ( Var ) η s y s
−30%13111190119.11.54342−46626.1%
−20%14121105135.61.61387−56527.5%
−15%14391078139.61.60375−58826.0%
−10%14601056141.51.60366−60125.0%
0%15011010150.21.65384−66425.6%
+10%14991013153.51.66365−69824.3%
+15%15031007157.91.66360−72224.0%
+20%15071003159.21.69353−74923.4%
+30%15071003168.11.75359−82723.9%
Table 9. Steady-state hydraulic and electromechanical quantities after a sudden change in the PAT head
Table 9. Steady-state hydraulic and electromechanical quantities after a sudden change in the PAT head
H P h   ( W ) N (rpm) U s r m s   ( V ) I s r m s   ( A ) P a   ( W ) Q r   ( Var ) η s y s
−50%45788775.00.7496−14521.1%
−40%642912100.21.01163−26725.3%
−30%842942118.81.23221−38626.2%
−20%1057970132.91.41277−49926.2%
−10%1297995144.81.55334−60125.3%
0%15011010150.21.65384−66425.6%
+10%17821059170.21.92453−88525.5%
+20%20471086180.62.08521−101725.5%
+30%23201091181.32.11569−104025.3%

Share and Cite

MDPI and ACS Style

Madeira, F.C.; Fernandes, J.F.P.; Pérez-Sánchez, M.; López-Jiménez, P.A.; Ramos, H.M.; Costa Branco, P.J. Electro-Hydraulic Transient Regimes in Isolated Pumps Working as Turbines with Self-Excited Induction Generators. Energies 2020, 13, 4521. https://doi.org/10.3390/en13174521

AMA Style

Madeira FC, Fernandes JFP, Pérez-Sánchez M, López-Jiménez PA, Ramos HM, Costa Branco PJ. Electro-Hydraulic Transient Regimes in Isolated Pumps Working as Turbines with Self-Excited Induction Generators. Energies. 2020; 13(17):4521. https://doi.org/10.3390/en13174521

Chicago/Turabian Style

Madeira, Filipe C., João F. P. Fernandes, Modesto Pérez-Sánchez, P. Amparo López-Jiménez, Helena M. Ramos, and P. J. Costa Branco. 2020. "Electro-Hydraulic Transient Regimes in Isolated Pumps Working as Turbines with Self-Excited Induction Generators" Energies 13, no. 17: 4521. https://doi.org/10.3390/en13174521

APA Style

Madeira, F. C., Fernandes, J. F. P., Pérez-Sánchez, M., López-Jiménez, P. A., Ramos, H. M., & Costa Branco, P. J. (2020). Electro-Hydraulic Transient Regimes in Isolated Pumps Working as Turbines with Self-Excited Induction Generators. Energies, 13(17), 4521. https://doi.org/10.3390/en13174521

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop