Centralized AirflowControl to ReduceOutput Power Variation in a Complex OWC Ocean Energy Network

A centralized airflow control scheme for a complex ocean energy network (OEN) is proposed in this paper to reduce the output power variation (OPV). )e OEN is an integrated network of multiple oscillating water columns (OWCs) that are located at different geographical sites connected to a common electrical grid. )e complexity of the OWC-OEN increases manifolds due to the integration of several OWCs and design of controllers become very challenging task. So, the centralized airflow control scheme is designed in two stages. In control stage-1, a proportional-integral(PI-) type controller is designed to provide a common reference command to control stage-2. In control stage-2, the antiwindup PID controllers are implemented for the airflow control of all the OWCs simultaneously. In order to tune the large number of control parameters of this complex system, a fitness function based on integral squared error (ISE) is minimized using the widely adopted particle swarm optimization (PSO) technique. Next, the simulation results were obtained with random wave profiles created using the Joint North Sea Wave Project (JONSWAP) irregular wave model.)e OPV of the proposed OWC-OENwas reduced significantly as compared to the individual OWC. It was further observed that the OPV of the proposed scheme was lower than that achieved with uncontrolled and MPPT controlled OWC-OEN. )e effect of communication delay on the OPV of the proposed OWC-OEN scheme was also investigated with the proposed controller, which was found to be robust for a delay up to 100ms.


Introduction
Many countries around the world with a long coastal zone have been actively exploring ocean waves as an alternate energy source over the past five decades [1][2][3][4][5][6][7][8][9]. European countries have taken a leading role in the research and development of this sector [1][2][3]. Research in ocean wave energy also started later in countries such as Australia, Brazil, China, India, Indonesia, Japan, Korea, Maldives, Taiwan, and the United States [4][5][6][7][8][9]. e majority of the ocean wave power plants across the globe use oscillating water column technology (OWC) [10][11][12]. With a lot of consistent efforts by the research community, this technology has matured enough in recent years. However, due to high output power variation (OPV), hurdles remain in its grid integration. e OPV is defined as the difference between the maximum and minimum power generated from the OWC plant. Reducing OPV is one of the major challenges of OWC. It requires advanced control systems to minimize OPV as much as possible.

Literature Review of OWC Control.
Over the years, the OWC has been tried and tested by several control schemes . e initial stage of the controller design for OWC was primarily on phase control [13] and latching control [14]. Subsequently, in [15], a controller by swinging rotorblades of the wells turbine was proposed to optimize OWC performance. A reactive controller for a fixed two-dimensional OWC was developed in [16]. In [17], an airflow control scheme for OWC-Wells turbine using throttle and bypass valves was discussed. Justino and Falcao [18] presented the rotational speed controller for OWC-Wells turbine and electrical generator. An approach for obtaining the maximum output power from the OWC plant was described in [19]. In [20], a power electronic controller was designed for an OWC system with an induction generator-Wells turbine. In [21], the rotor side controller with constant frequency and variable speed operation of OWC was presented.
After the establishment of the OWC plant at the Mutriku harbor in Spain, a neural network network-based controller for OWC module was proposed in [22], which was supposed to change the external rotor resistance values to vary the induction generator slip. e PID-based control approaches for OWC plants were developed in [23] for airflow and rotor speed control. In [24], the control schemes proposed in [23] were validated on the experimental module. Neural network-based optimum speed reference for PI-based rotational speed control of the OWC system was proposed in [25]. In [26], an analysis was executed on fault-ride-through capability of OWC stations. It was equipped with doubly fed induction generator (DFIG) and an airflow-based controller. In [27], a sliding-mode control of OWC plants for maximizing output power of DFIG was designed.
In recent years, many advanced control techniques for OWC airflow as well as turbine rotor speed control techniques have been developed. In [28], a control technique was developed for efficiency enhancement of Wells turbine-OWC devices. A comprehensive survey of electrical controllers for ocean energy devices including OWC was discussed in [29]. In [30], the controller for merging local energy storage with OWC was proposed. A rotational speed optimization approach was applied in [31] to attain maximum power point tracking (MPPT) of OWC. In [32] as well, a latching control scheme was developed for MMPT of a floating OWC wave energy converter. An optimal speed controller based on the fuzzy-backstepping method was designed in [33].
ree Lyapunov-based nonlinear controllers were proposed in [34] for OWC plant emulation, rotational speed control, and DC link voltage control of DFIG's grid side converter (GSC), respectively. e flow controller was designed in [35] for Wells turbines for harnessing maximum wave power using OWC. An eventtriggered controller for the OWC energy plant was presented in [36] wherein the main objective was to minimize control updates when the controller and plant interact with the other via communication channel. In [37], a fuzzy gain scheduled and PI-type airflow controller for an OWC, which was an advanced version of the airflow controller, were presented earlier in [23,24]. Recently, some artificial intelligence-based airflow controller techniques have been proposed [38][39][40]. A self-adaptive global best harmony search algorithm-based airflow control of OWC was proposed by [38] wherein four variants of the harmony search algorithm were implemented and tested to optimize the control design of the PID controller in the airflow control scheme. In [39], another approach for airflow control was proposed which utilized symmetry-breaking concept to design the controller. Furthermore, an artificial neural network-based airflow controller was designed using surface elevation measurements [40]. is study considered real measured wave input data and generated power of the NEREIDA wave power plant. For more details on the control of OWC plants, please refer [41].

Main Contributions.
In the above-discussed control schemes, the output power of OWC can be maximized successfully using rotor speed regulation wherein actual rotor speed follows optimum reference speed. Also, output power can be controlled to some extent by using an air flow controller but the high OPV is still a huge challenge in OWC. To the best of the authors' knowledge this issue has not been addressed so far by forming a large network of OWC. e randomness in ocean wave input at the OWC plant is a serious concern. If the multiple OWCs operating at different geographical locations are integrated, then the same randomness of ocean waves might be helpful in reducing OPV as compared to that of an isolated OWC. e OPV can also be reduced further by suitably designing a centralized airflow controller for regulating the airflow of OWC chambers situated at different geographical sites. However, the design complexity of the controller for this system is increased.
Hence, in this paper, the issue of high OPV has been addressed by integrating multiple OWC sites located at different geographical locations. Hereafter, the integrated OWC sites would be referred to as OWC ocean energy network (OEN). Now, two control schemes would be designed for OWC-OEN: (i) localized MPPT control and (ii) centralized airflow control for reducing OPV. e localized MPPT control law would again be consisting of two sections: (i) an algorithm for optimum reference speed calculation and (ii) backstepping controller for rotor speed regulation. e centralized airflow controller would provide common commands to all individual OWC sites. e centralized airflow controller is a two stage control scheme. In control stage-1, a proportional-integral-(PI-) type controller is designed to provide a common reference command to control stage-2. In control stage-2, the antiwind PID controllers are implemented for the airflow control of all OWC sites. Suppose, there are N number of OWC sites included in OWC-OEN, two unknown parameters of the PI controller in control stage-1, and three unknown PID parameters of each OWC site in control stage-2. en, the 3N + 2 number of controller design parameters would be required to be properly selected so that OPV of OWC-OEN can be reduced considerably. In order to tune the 3N + 2 parameters, an integral squared error-(ISE-) type fitness function is defined. e fitness function is minimized with the help of particle swarm optimization (PSO) algorithm [42][43][44][45]. e PSO method is very instrumental in optimizing the parameters of complex systems.
2 Complexity e numerical simulations were performed on MATLAB ® platform using Joint North Sea Wave Project (JONSWAP) irregular wave model [46]. Several JONSWAP wave profiles are generated randomly using Wave Analysis for Fatigue and Oceanography (WAFO) MATLAB ® toolbox [47]. e performance of an arbitrarily selected OWC site has been analyzed for relevant OWC parameters such as rotor speed, control signal, Wells turbine flow coefficient, and DFIG output power of controlled/uncontrolled OWC. e performance of the proposed control system of OWC-OEN has been validated using a newly introduced performance index referred to as OPV index. e OPV index of all OWC sites, uncontrolled OWC-OEN, localized MPPT control of OWC-OEN, and centralized airflow control of OWC-OEN are compared. Also, the effect of transmission delay has been studied by introducing the delay in the OWC-OEN. e rest of the paper would be discussed with following sections. Wave-to-wire modelling and control of an individual OWC site in Section 2; design of OWC-OEN and centralized airflow controller in Section 3; discussion on simulation results in Section 4 ,and concluding remarks and future scopes in Section 5.

Wave-to-Wire Modelling and Control of an
Individual OWC Plant Figure 1 displays the schematic representation of OWC plant configuration and its control schemes. As shown in Figure 1, the OWC consists of a chamber, turbine, electrical generator, and controller [34]. e waves of the sea hit the front of the chamber, causing an increase and fall in the water. In the chamber, the air moves back and forth to a circular cavity above the chamber. us, the Wells turbine generator machine, which is mounted in the cavity, is operated by twoway air flow [16]. e controllers are certainly an essential component of OWC systems. ere are mainly three control configurations of OWC plants as uncontrolled, MPPT control, and airflow control. In uncontrolled OWC plant configuration, OWC operates without any control action. e output power generated in uncontrolled condition directly depends upon the sea wave input and the characteristics of turbine generator. In case of MPPT control, the rotor speed is regulated as per the reference rotor speed wherein the reference speed depends upon the sea wave input to the plant. is mechanism helps in extracting maximum output power from the OWC plant. Finally, in airflow control of the OWC plant, the control valve is used to restrict the airflow inside the OWC chamber in such a manner that desired electrical power is generated. Next, we will discuss OWC device modelling and its MPPT control.

Ocean Waves.
Most OWC experiments have chosen regular types of ocean waves for OWC output evaluation, but real sea wave actions and irregularities are highly unpredictable. In order to deal with this problem, researchers have therefore developed a number of practical sea wave models [46,47]. In this analysis, the irregular JONSWAP wave pattern is used to test the efficiency of the proposed controls. It is one of the commonly used models for marine technical research. e characterization parameters for JONSWAP are significant wave height (H m0 ), peak wave period (T p ), and peakedness parameter (c). e JONSWAP model gives best characterization in the range 3. 6 ���� H mo < T p < 5 ���� H mo and 1 < c < 7 [34]. e JONSWAP wave spectrum for the values of H m0 , T p , and c given in Table 1 is shown in Figure 2 wherein the highest frequency component is around 0.5 rad/s.

OWC Chamber Model.
e OWC consists of a closed room with four side walls with openings at top and bottom ( Figure 1). It is partly immersed in water and waves hit the lower section of the chamber. As the water level rises and falls, a two-way air flow is produced. e water level change is also referred to as the wave height. e air velocity mathematical definition is provided by [36] e air velocity definition established by (1) is the input for the Wells turbine. Here, it has been assumed that A valve � u v . It means that the control valve has unity transfer function for a simplified case or does not have any actuator delay.

Wells Turbine Model.
e Wells turbine [24,48] is a turbine of unique features and can be seen in Figure 3. is rotates in the same direction by input of two-way air flow,  Complexity 3 but has a restriction of stalling.
e power output falls considerably because of the stalling. A control unit is therefore necessary to prevent the stall.
Turbine torque description is given as follows [24]: where f t (ϕ) is stated as follows: where C t is the Wells turbine torque coefficient and depends on ϕ as displayed in Figure 3. e ϕ is given by As shown in Figure 4, ϕ ≤ ϕ th � 0.3 provides the highest torque coefficient to deliver the maximum torque and therefore the maximum output power. e turbine is attached to a DFIG. e coupling function of the turbine generator can thus be expressed as follows: where ω r has initial condition ω r0 .

DFIG Dynamics.
is analysis considered a directquadrature (dq) dynamic DFIG model. e dq dynamic model's benefit is that every three-phase voltage or current in the stationary frame can be interpreted as DC voltage or current in a synchronous revolving frame [15,16]. e state equations of dq-axis stator and rotor flux states of DFIG are given as follows [20]: where K � L s L r − L 2 m . ψ ds , ψ qs , ψ dr , and ψ qr are stator and rotor dq flux quantities. R s and R r are stator and rotor  4 Complexity resistances, whereas L s , L r , and L m are the stator, rotor, and mutual inductances. ω e is the stator supply frequency, v ds , v qs , v dr , and v qr are stator and rotor dq voltages. e flux states ψ ds , ψ qs , ψ dr , and ψ qr are subjected to initial conditions ψ ds0 , ψ qs0 , ψ dr0 , and ψ qr0 , respectively. e electromagnetic torque and output power expressions are where M � − (3/2)(p/2)(L m /K) and p is number of poles of DFIG.
e dq-axis voltages of the rotor are converted into abcaxis voltages of the rotor. Furthermore, the abc-axis voltage signals will be sent to pulse width modulation (PWM), the gate pulses of which will be sent to DFIG's DC/AC rotor side converter (RSC) [11,12,23].

Localized MPPT Control of an Individual OWC Plant.
e airflow control and rotor speed regulation blocks, as shown in Figure 1, are the two most important controllers for OWC control. As mentioned earlier, the rotor speed regulator is considered for localized MPPT control of the OWC plant.
e MPPT control can be achieved using linear or nonlinear controller in which the actual rotor speed is forced to track an optimum speed reference [24,25,27,31,33,34,36]. e optimum reference speed calculation algorithm is designed in such a manner that it would always lead to maximum output power. In this work, a backstepping-type rotor speed controller and optimum reference speed calculation algorithm have been designed using the approach described in [36]. e optimum reference speed calculation algorithm [36] is presented next. It is calculated using air velocity information in the following steps: (i) Step 1: calculate the value of V x using (1). (ii) Step 2: calculate peak values of V x as follows: if Step 3: pass V xp through zero-order-hold (ZOH) as follows: (iv) Step 4: for threshold value of the flow coefficient, φ th � 0.3, calculate the reference speed using (4) as follows: (v) Step 5: calculate ω 1d by limiting the minimum and maximum values of ω ref as follows: if where ω e is the minimum value of ω 1d and ω rp is the maximum value of ω 1d . (vi) Step 6: in order to avoid any sudden change in the reference value, pass ω 1d through a low-pass filter. e new reference value is calculated as follows: e frequency domain representation of h f is given as follows: As a result, the z 1d is the optimum reference speed to be used for proposed controllers. Next, the backstepping-type rotor speed regulator would be designed to force the actual rotor speed to track the optimum reference speed z 1d . e state space model from [36] is given below: where Furthermore, simplifying the above equation, we have where A step by step design approach is applied for the backstepping controller. For a second-order system expressed by (15) and (16), the controller is designed in two Complexity 5 stages. First, a virtual controller, z 2d , is to be designed. en, z 2d would be used for designing final control law u r . For designing the virtual control law, z 2d , let us choose the error component as follows: Next, the first order derivative of (17) is written as follows: Furthermore, by adding and subtracting k 2 z 2d term, (18) becomes e virtual control law z 2 d is chosen as follows: where σ 1 > 0. Define Substituting (20) and (21), (19) can be expressed as follows: _ Next, the derivative of (21) is given by For designing the control law which ensures the closed loop stability of systems (15) and (16), a Lyapunov function candidate, V lpf , is chosen as follows: e first order time derivative of (23) is expressed as follows: e control law u r gives where σ 2 > 0.
Here, (27) becomes negative semidefinite and _ V lpf does not contain any trajectories of error states z 1 and z 2 other than trivial trajectory z 1 � z 2 � 0. en, error states z 1 and z 2 converge asymptotically to zero. Hence, the system given by (15) and (16) with control law u r given in (26) is asymptotically stable.
In Section 3, the airflow control concept would be extended to the airflow control of OWC-OEN. As per the analysis of this work, the formation of OWC-OEN reduces the OPV as compared to the individual OWC plant. is could be further reduced by a centralized airflow controller to restrict the power generation as per the desired or reference power. e power output signal is sent by communication network to the controller. e controller output is again sent back to the OWC using the same channel. Needless to mention that the communication channel is a digital platform wherein the analog-to-digital (A/D) and digital-to-analog conversion (D/A) of plant and controller data also takes place simultaneously.

Formation of OWC-OEN and Centralized Airflow Controller Design
In this section, the concept depicted in Figure 1 is extended to a network of multiple OWC sites at different geographical locations.

Formation of OWC-OEN.
e OWC-OEN is proposed to be implemented in a smaller coastal region not greater than 50 km. OWC sites of OWC-OEN are independent from each other and are having at least 2 km distance between them. Typically, the distance covered by the wave front is a few hundred meters. As per the proposed configuration, two sites are separated by at least two kilometers and thus, the wave profiles at these two sites can be assumed to be uncorrelated. Next, each OWC site receives distinct sea wave input due to changed geographical location.
It is also assumed that each OWC site of OWC-OEN supplies the current to a common electrical grid, as shown in Figure 5, wherein the grid voltage and frequency are assumed to be constant. Hence, the total output power could be taken as a sum of output power of all individual OWC sites provided there is common grid voltage provided to each OWC. Suppose, there are N number of OWC sites connected to the grid network and are called as OWC-1, OWC-2, . . ., OWC-N. e output power from these OWC sites is denoted as P g1 , P g2 , . . . , P gN . e total output power, P gtotal , is given as follows: where T ek is the electromagnetic torque of kth OWC site; ω rk is the rotor speed of kth OWC site; V sk is the grid voltage of kth OWC site and V s1 � V s2 � · · · � V sN ; and I sk is the current supplied by kth OWC site to grid and I s1 ≠ I s2 ≠ · · · ≠ I sN . e centralized airflow control unit, as displayed by Figure 5, is located in the land area and is not very far (few tens of kilometers) from different OWC sites of OWC-OEN. Each OWC site is also equipped with MPPTcontrol at a local level so as to extract maximized output power. We consider the use of high speed fiber-optic cables for communication between the central airflow control unit and the OWC-OEN. 6 Complexity e output power data from OWC sites is transmitted to the controller and vice versa. e one-way delay time associated with fiber-optic cables for communication over large distances (several hundreds of kilometers) is typically around 100-150 ms [49]. As the distance between the central control unit and the OWC sites is only a few tens of kilometers, the communication delay would be much lower than 100 ms. In order to ensure a rigorous analysis of simulation results, three communication delays on higher side are considered as τ � 100 ms, 500 ms, and 1000 ms. ree types of configurations of OWC-OEN from the control point of view are considered in this study which are uncontrolled, localized MPPT control, and centralized airflow control. In uncontrolled configuration, the controllers are assumed to be absent at local (i.e., MPPT control) as well as central level (i.e., airflow control). erefore, the output power of each OWC site is delivered to the grid without performing any control action. In localized MPPT control of OWC-OEN, each OWC site is coupled with MPPTcontrol at the local level and the maximized power is extracted. en, the maximized output power from each OWC site is delivered to the grid. In centralized airflow control of OWC-OEN, each OWC site is integrated with control valves. e control valves are operated at the central level using the airflow controller that controls the amount of total power delivered to the grid.

Centralized Airflow Controller Structure.
e internal structure of central airflow control unit of OWC-OEN is shown in Figure 6. It consists of two control stages as control stage-1 and control stage-2. Control stage-1 consists of a PItype controller. e total output power, P gtotal , is compared with the reference total output power, P gtotal− ref , and error, e gtotal , is fed to the PI controller. e PI controller provides a common control command P g− ref to control stage-2. e internal structure of antiwindup PID-type controllers used for control stage-2 is shown in Figure 7.
e error e gk (t) in Figure 7 is given as follows: where k � 1, 2, . . . , N. P g− ref is the common reference command obtained from the central PI controller as depicted in Figure 6 and P gk is the output power obtained from kth OWC site via communication network. e suffix k is used to denote the PID k controller that belongs to kth OWC site. Similarly, K pk , T ik , and T dk are the proportional gain, integral time constant, and derivative time constant of the PID k , respectively. Next, the output of PID k controller, u vk , works as the input to the airflow control valve of kth OWC site. Antiwindup PID controller output characterises the effective valve area. erefore, the effective turbine duct area fluctuates according to the deviations in the valve area, and the required output power is achieved by adjusting the air flow rate. For simplicity, the transfer ratio between the valve input (control signal) and the valve output was taken as unity. e saturation effect is due to the valve that works as the final control actuator. erefore, the PID output is constrained to prevent the impact of valve saturation. e saturation range is taken between 0 and 1. e PID output has no effect on the valve reaction when the valve becomes saturated.
e valve goes to a steady state, but the error between the actual and desired values causes the amplitude of the PID integrator to increase continuously. It enhances the control effort without affecting the valve response. To avoid this, the antiwindup PID control scheme has been implemented wherein the integral action of PID is reduced via algebraic feedback, as shown in Figure 7. e antiwindup  PID controller scheme was also proposed earlier in [23,24] for the OWC plant.

OPV Index of Individual OWC Sites and OWC-OEN.
e OPV index is the difference between the maximum and minimum values of the actual output power of OWC divided by the total rated power of the OWC. e main objective of this paper is to reduce the deviations in the output power to make it suitable for electrical grid integration. In other words, the OPV index is the measurement of the suitability of the output power for the electrical grid, and it is desirable to keep this index as minimum as possible. Next, we define an OPV index as follows: where P g− max and P g− min are the maximum and minimum values of output power of an OWC, respectively. Now, for individual OWC sites, the OPV index is defined as follows: Next, the OPV for complete OWC-OEN is defined as follows: where P gk− rated is the rated power of the kth OWC site and P gtotal− rated is the rated power of OWC-OEN. e OPV index expressed in (10), (11), and (12) is very important for measuring the improvement in the power variation of OWC-OEN. e objective is to obtain the minimum possible value of OPV index by suitably designing the centralized airflow controller.

Centralized Airflow Controller Tuning Using PSO.
e complex controller tuning involves 3N + 2 unknown parameters to be determined. erefore, we choose PSObased optimization approach for tuning controller parameters. Now, we define an ISE type fitness function as follows: where e gtotal is defined as follows: erefore, the ISE can be written as follows: where T s is simulation run time of the OWC-OEN model. Next, the objective is to minimize ISE in such a way that all 3N + 2 parameters of OWC-OEN controllers are properly selected. So, the PSO algorithm has been applied for optimum selection of control parameters of OWC-OEN. e population-based swarm intelligence method is widely adopted to find the optimal solution when solving complex optimization problems with a large search area. e PSO approach [39,40] is a population-dependent search strategy that has a wide target field to address the problem of optimization. is approach represents the behavior of birds and fish schools. In a PSO system, an individual bird or fish is known as particles and each particle has its position and velocity. Now, particle travels in multidimensional search space depending on its own knowledge and the surrounding     [41,42]. However, there are many optimization algorithms available in the literature and anyone of them could be utilized for tuning purpose. e comparison of optimization algorithms is not taken into account in this study because the proposed research work is mainly focused on reduction in OPV using centralized airflow controller. e velocity and position update of each particle can be calculated using the current velocity and position as follows: Position update equation is given by where and k p− max � maximum number of iterations, k � iteration number, j � particle number, V k j � velocity of jth particle in kth iteration, W p � inertia weight factor, W p− max � Maximum inertia weight factor, W p− min � Minimum inertia weight factor, c p1 � cognitive acceleration factor, c p2 � social acceleration factor, rand () � random numbers uniformly distributed in the range (0, 1), X k j � position of jth particle in kth iteration, p best � local best position, and g best � global best position. e block diagram of PSO algorithmbased tuning of OWC-OEN controllers is given in Figure 8.

Discussion on Simulation Results
In this section, the centralized airflow controller has been validated through numerical simulations. e simulation parameters used for OWC-OEN airflow control and rotor speed regulation schemes are provided in Table 1. e number of OWC sites are chosen to be N � 20. e WAFO toolbox in MATLAB ® has been used for generating JONSWAP spectrum-based twenty distinct sea wave scenarios for twenty OWC sites located at different locations. e sea scenario SW1 (Figure 9) is chosen for demonstrating the MPPT outcome of OWC-1, whereas all sea scenarios, i.e., SW1, SW2, . . . , SW20, are selected for evaluating the performance of OWC-OEN with MPPT, uncontrolled OWC-OEN, and OWC-OEN with centralized airflow control. e wave height and corresponding airflow velocity for sea scenario, SW1, is shown in Figure 9.

Performance Analysis of MPPT Control.
e MPPT performance of OWC-1 at the local level has been evaluated in this section. e optimum reference speed, ω ref− OWC1 , is calculated first using the algorithm suggested in [36]. is reference is generated in a manner that the turbine flow coefficient always remains below the permissible threshold limit to escape the turbine stalling. If turbine stalling is avoided, then the output power can be maximized to its best possible level. Next, the backstepping-based nonlinear controller is implemented for rotor speed control that in turn provides the efficient tracking of reference speed by actual rotor speed. is is observed in Figure 10, wherein the actual rotor speed of OWC-1 very closely follows the optimum reference speed. e rotor speed for uncontrolled OWC-1 is clearly representing the turbine stalling case. e control effort needed to manipulate the rotor speed is also shown in Figure 10 which is having a range of ±150 volts and is required only when the reference value is changing.
Next, turbine flow coefficient and electrical power waveforms are shown by Figure 11. e turbine flow coefficient for MPPT-controlled OWC-1 always remains below the threshold level, i.e., 0.3, whereas for uncontrolled OWC-1, it exceeds the threshold and causes the stalling in Wells turbine. Similarly, the output power gets maximized in case of MPPT control as compared to uncontrolled OWC-1. However, the OPV is very high as seen in power waveform for OWC-1. e OPV index of all OWC units is provided in Start Randomly initialize the particles Run the OWC-OEN MATLAB model (.slx) calculate the fitness value (i.e., ISE) of each particle (i.e., OWC-OEN controller parameters) Find p best particles & fitness value of each p best particle  Table 2 which would be compared with the OPV index of OWC-OEN later. As shown in Figure 12, the total output power has been analyzed with and without MPPT control after forming OWC-OEN. For uncontrolled case, the output power obtained has a peak value of P gtotal− max � − 373.3 kW, whereas MPPT control provides the peak value of P gtotal− max � − 587.4 kW. Again, the OPV index of OWC-OEN is not at an acceptable level. erefore, the centralized airflow controller to be discussed next was required to deal with high OPV issue of OWC-OEN. e output power waveforms have negative amplitude which represents the generation mode of DFIG. In the generation mode, the negative peak of output power actually represents the P g− max . e value of output power that is close to time axis is treated as P g− min . In the motor mode of DFIG, the power waveforms would have positive amplitude.

Performance Analysis of Centralized Airflow Control.
e centralized airflow controller, as shown in Figure 6, is designed using PSO. For N � 20, there are 62 control parameters which have been tuned with an objective to minimize the ISE as given in (35). e ISE was minimized by applying PSO algorithm as per the flowchart shown in Figure 8. e PSO parameters for running optimization are given in Table 3. 50 Next, the optimized control parameters were obtained by running the optimization algorithm for 100 iterations. e tuned parameters of centralized airflow controller are provided in Table 4 e total output power is more concentrated in the lower range (i.e., below − 250 kW) and hence, the centralized airflow controller is found to be more effective in the lower range and the effectiveness gets reduced substantially at higher ranges. e total output power, P gtotal , and corresponding error waveforms for different references are shown in Figure 13. e total output power of OWC-OEN successfully tracks the reference power. Although for lower as well as higher  Table 5. It can be seen in the table that the centralized airflow controller leads to lowest OPV index among all configurations of OWC-OEN. e OPV index with MPPT controlled OWC-OEN is higher than that of uncontrolled OWC-OEN as provided in Table 5. e     12 Complexity  centralized airflow controller is also tested with different cases of communication delay in OWC-OEN in next section.

Performance Analysis of Centralized Airflow Control with
Communication Delay. e centralized airflow controller has been validated with different communication delays of τ � 100 ms, 500 ms, and 1000 ms. e total output power and corresponding error waveforms are shown in Figure 14.
e OPV index increases with increase in τ value as observed from the waveforms. So, it is required that the communication network should be able to transfer data with a very high speed in order to avoid any performance deterioration occurring due to delay. However, it is very much possible nowadays to transmit data with bandwidths in gigahertz range with advanced communication technology. It is also to be noted here that each OWC site may have different communication delays which needs further attention. It is observed that the total output power waveform in Figure 14(a) seems very similar to that of Figure 13 However, there is clear difference in performance of both ceases in terms of the OPV index, as given in Table 5.

Concluding Remarks and Future Challenges
is paper presented a novel concept of OWC-OEN and its centralized airflow control to mitigate the OPV issue of ocean energy systems. e centralized airflow controller consisting of two control stages was designed using PSO due to complexity of tuning parameters. e PSO helped in choosing large number of control parameters of a very complex OWC-OEN by minimizing an ISE-type fitness function. A comprehensive simulation analysis was presented to validate the proposed airflow control scheme. e JOWNSWAP-based irregular wave patterns were randomly selected as the input to individual OWC sites. Also, an OPV index was defined to compare the output power waveforms of individual OWC sites and OWC-OEN. e OWC-OEN with centralized airflow control delivers the lowest OPV index as compared to uncontrolled OWC-OEN and MPPT controlled OWC-OEN. e centralized airflow controller works reasonably well when OWC-OEN is under the effect of communication delay. As with the case of any other spatially distributed control network, large communication delays degrade the performance of the OWC-OEN. e communication delays considered in this study are assumed to be uniform for all OWC sites. e unequal communication delays for OWC sites located at different geographical locations would further increase the complexity of OWC-OEN. It is an interesting subject, which needs thorough investigation in further studies on OWC-OEN.

Data Availability
All the data are included within the article. Additional data can be provided by the corresponding author upon request (sunil.mishrafet@kiit.ac.in). 14 Complexity