Experimental Modeling of an Isolated WECfarm Real-Time Controllable Heaving Point Absorber Wave Energy Converter

: To offer point absorber wave energy converters (WECs) as a bankable product on the marine renewable energy market, multiple WECs will be installed together in an array conﬁguration. The wave energy community (research and industrial) has identiﬁed the urgent need for available realistic and reliable data on WEC array tests in order to perform a better WEC array optimization approach and in order to validate recently developed (non-linear) numerical models. The ‘WECfarm’ project is initiated to cover this scientiﬁc gap on necessary experimental data. The ‘WECfarm’ experimental setup consists of an array of ﬁve generic heaving point-absorber WECs. The WECs are equipped with a permanent magnet synchronous motor (PMSM), addressing the need for WEC array tests with an accurate and actively controllable power take-off (PTO). The WEC array control and data acquisition are realized with a Speedgoat Performance real-time target machine, offering the possibility to implement advanced WEC array control strategies in the MATLAB-Simulink environment. The presented article describes the experimental setup, the performed tests and the results of the test campaign using a single, isolated ‘WECfarm’ WEC in April 2021 at the wave basin of Aalborg University (AAU), Denmark. A Coulomb and viscous friction model is determined to partly compensate for the drivetrain (motor, gearbox, rack and pinion) friction. A system identiﬁcation (SID) approach is adopted considering the WEC system to be composed of two single input single output (SISO) models, the radiation and the excitation model. Radiation tests yield the intrinsic impedance. Excitation tests yield the excitation frequency response function. Adopting an impedance matching approach, the control parameters for the resistive and reactive controller are determined from the complex conjugate of the intrinsic impedance. Both controllers are tested for a selection of regular wave conditions. The performed experimental test campaign using an isolated ‘WECfarm’ WEC allows a full evaluation of the WEC design prior to extending the setup to ﬁve WECs. Within the ‘WECfarm’ project, an experimental campaign with a ﬁve-WEC array in the Coastal and Ocean Basin (COB) in Ostend, Belgium, is under preparation.


Introduction
A point absorber wave energy converter (WEC) consists of a floating or submerged body to capture energy from different wave directions. The point absorber diameter should preferably be in the range of 5-10% of the prevailing wavelength [1]. Due to its ability to absorb energy from different directions, this WEC type is particularly suitable to put in arrays. In a WEC array, hydrodynamic interactions between the WECs occur through radiation and diffraction of waves. Both constructive and destructive interactions will occur between individual WECs within a WEC array called near-field interactions.
Göteman et al. define directions for future research necessary for a better WEC array optimization approach [2]. There is a need for available real life data for the validation of WEC array modeling and optimization. This is a research gap for the full wave energy sector, not only in the optimization of WEC arrays. In computational fluid dynamics-based numerical wave tanks (CNWT), the Power Take-Off (PTO) system is mostly modeled as a linear spring-damper system, not representing realistic PTO dynamics and inefficiencies and undermining the overall model fidelity [3]. To validate CNWT considering WECs, it is desirable to incorporate a realistic, nonlinear PTO model. As the computational power capabilities increase yearly, so do the numerical model's capabilities, stressing the need for experimental data to validate the model. However, publicly available databases from WEC array experiments are scarce. Vervaet et al. identified 17 experimental campaigns on point absorber WEC arrays, carried out during the last decades [4]. This limited number of experimental campaigns is due to the high cost of constructing and testing in wave basin facilities, as well as due to the complexity of the experiments and related instrumentation [5]. Therefore, the 'WECfarm' project aims to deliver a dataset to cover the research gap on the need for publicly available real life and reliable data to validate these new advanced numerical models. Vervaet et al. discuss the state of the art in physical modeling of point absorber WEC arrays and the identification of research gaps, resulting in design specifications of the 'WECfarm' experimental setup [4].
The 'WECfarm' experimental setup consists of an array of five heaving point-absorber WECs, designed as a unique test bench for future innovative WEC array research, able to address the current requirements and research gaps on physical WEC array testing. Given the limited number of five WECs, the WECfarm WEC array is not classified as a large WEC array, as the Manchester Bobber 25-WEC array [6], the PerAWaT project 24-WEC array [7] and the WECwakes project 25-WEC array [8]. Vervaet et al. discuss the features of the experimental setup, for which the most important ones are summarized below [4]. The WEC buoy is designed to be generic, being a truncated cylinder with a draft of 0.16 m and a radius of 0.30 m. The high diameter-to-draft ratio of 3.75 yields a flat Response Amplitude Operator (RAO) response and high resonance bandwidth, enhancing WEC-WEC interactions. The WECs are equipped with a Permanent Magnet Synchronous Motor (PMSM), addressing the need for WEC array tests with an accurate and actively controllable PTO. The air bushings linear guiding system excludes guiding friction in the power absorption measurements. The WEC array control and data acquisition are realized with a Speedgoat (Speedgoat, Köniz, Switzerland) Performance real-time target machine, offering the possibility to implement advanced WEC array control strategies in the MATLAB-Simulink environment. Wave basin testing with 'WECfarm' WEC arrays targets to include long-and short-crested waves and extreme wave conditions, representing real sea conditions. Within the 'WECfarm' project, two experimental campaigns have been performed at the Aalborg University (AAU) wave basin: (a) testing of the first WEC in April 2021, addressed in the presented article; (b) testing of a two WEC array in February 2022. An experimental campaign with a five WEC array, in the new wave basin; the Coastal and Ocean Basin (COB) in Ostend (Belgium) [9], is scheduled in 2023.
Friction characterization tests are performed to quantify the drivetrain (motor, gearbox, rack and pinion) introduced friction, whereafter a Coulomb and viscous based friction model for partial compensation of the drivetrain friction is implemented in the MATLAB-Simulink control model. Beatty et al. used a Proportional Integral (PI) force control to minimize the error between the target and measured forces for the physical model of a Wavestar WEC [10]. However, Bacelli et al. stress that closing the loop around a force sensor may induce negative consequences for the design of higher level control loops [11]. The presented friction compensation methodology provides an alternative for closing the PTO force feedback loop around the force sensor.
Coe et al. present a WEC control design based on the principle of impedance matching [12]. The control parameters yield from the complex conjugate of the intrinsic impedance, determined by radiation system identification (SID) tests, where the WEC is excited by a torque noise signal in calm water. Bacelli et al. performed SID tests for an isolated heaving point absorber WEC [13]. This article discusses the application of SID tests on an isolated 'WECfarm' WEC. The SID methodology will be extended to arrays with two to five WECs. Application of the impedance matching methodology for WEC arrays will yield valuable data and insights on WEC-WEC interactions and WEC array control optimization. Within the presented test campaign, the resistive and reactive controller are tested for a selection of regular wave conditions. This testing campaign allows a full evaluation of the WEC prior to extending the setup to five WECs.
Section 2 provides a detailed overview of the experimental setup, with a focus on the 'WECfarm' WEC, instrumentation and wave basin setup. The experimental results of the drivetrain friction model characterization tests, the SID tests and the power absorption tests are discussed in Section 3. A summary of the findings and conclusions are presented in Section 4.

WECfarm WEC and Instrumentation
For a detailed discussion on the design of the 'WECfarm' five-WEC array, the reader is referred to [4]. In this article, we consider a single, isolated 'WECfarm' WEC. Figure 1 shows a 3D rendering of the final design of the device. The used right-handed coordinate system has its origin at the intersection of the still water level (SWL) with the vertical axis through the center of the WEC buoy. This allows us to express displacements of the WEC buoy relative to the SWL. The x-axis corresponds with the positive wave propagation direction. The y-axis follows from the motor sign convention: a positive torque results in a downward motion of the WEC buoy. Therefore, the z-axis is pointed downwards to define positive forces, displacements, velocities and accelerations.   To exclude friction in the linear guiding, 40 mm OAV (OAV Air Bearings, Princeton, NJ, USA) air bushings are used. The air bushings are characterized by a load versus pressure curve, where one air bushing can cope with a maximum radial load of 720 N, for a nominal pressure of 5.5 bar [14]. A configuration of three OAV 40 mm air bushings guarantees a permanent layer of air between the guide shafts and the bushings for the most extreme wave conditions, resulting in zero-friction linear guiding on the condition of proper alignment. The PTO system of the WEC is designed as a PMSM connected to a gearbox powering a rack and pinion system. The pinion pitch circle radius R pinion is equal to 0.0212205 m. A Wittenstein (Wittenstein, Igersheim, Germany) single-stage gearbox 'NPR 025S-MF1-4 -2E1-1S' with ratio i = 4 is connected to a Beckhoff (Beckhoff Automation, Verl, Germany) PMSM 'AM8542-2E11-0000' with an inertia of 6.17 kg cm 2 , a rated torque of 3.97 Nm and a rated speed of 1200 RPM, for 230 V AC power supply [15,16]. The velocity on the pinion will be four times less than on the motor shaft, while the torque on the pinion will be four times more than on the motor shaft. The Beckhoff PMSM is powered and controlled by a Beckhoff motor drive type 'AX5103-0000-0212'. The hydrodynamic part of the WEC consists of an Acrylonitril-Butadieen-Styreen (ABS) thermofolded truncated cylindrical buoy, covered with a Polymethylmethacrylate (PMMA) plate. Figure 2a shows a 2D rendering of the WEC buoy with its dimensions. The WEC buoy is 0.32 m high and designed with a draft of 0.16 m. This draft corresponds with a submerged volume of 0.03683 m 3 . Therefore, the mass of the WEC buoy and hydrodynamically activated parts on top of it is 36.83 kg.  Figure 3 shows a scheme of the data acquisition and control flow for the 'WECfarm' setup with the isolated WEC, with a legend indicating the signal type. The MATLAB-Simulink real-time control model is built on the host PC and loaded on the Speedgoat Performance real-time target machine by Ethernet communication. This target machine runs the Simulink model and processes the input/output (I/O) at a sample frequency of 1000 Hz. In this context, the high sample frequency corresponds with the defined 'real-time' terminology. For each test, the various time series of each logged Simulink signal are saved within a single MATLAB structure.
A scheme on the bottom of Figure 3 shows the sensor input for the Speedgoat IO133 terminal board. The accelerometer ADXL335 (Analog Devices, Norwood, MA, USA) is used to measure the acceleration of the WEC buoy in the heave direction and is attached on top of the rack, the furthest position on the WEC from the water. The accelerometer has a linearity of ±0.3 % of the Full Scale Output (FSO) [17]. Three Tedea Huntleigh (Vishay Precision Group, Malvern, PA, USA) 50 kg load cells are placed between the hydrodynamic part (=the buoy) and the electromechanical part (=the motor) to measure the actual applied forces. The load cells with accuracy class C3 have a total error (per OIML R60) of 0.020% of the rated output [18]. A configuration of at least three load cells is required to avoid torsion and bending influencing the measurements. The mass of the three load cells together is equal to 0.682 kg, the mass above the load cells m top is equal to 27.610 kg and the mass below the load cells is equal to 8.534 kg, which results in a total hydrodynamically activated mass m of 36.83 kg. A TLE analog weight transmitter (Laumas Elettronica, Montechiarugolo, Italy) is used to amplify these three analog signals and to sum them to one analog signal. In case the WEC is locked, the wave heave excitation force F e can be measured. In case the motor is active, the load cells measure the PTO force F PTO . The upper micro switch and the lower micro switch are used as safety limit switches. It is necessary to limit the amplitude of the WEC buoy displacement to prevent the guiding system damaging the structure. The laser sensor Micro-Epsilon optoNCDT 1420-500 (Micro-Epsilon, Ortenburg, Germany) is installed as a backup for the motor encoder to measure the displacement of the WEC buoy relative to the SWL. The laser sensor has a linearity of ±500 µm, equivalent to ±0.1% of the FSO [19]. Moreover, the laser sensor can be used for displacement measurements for tests without the motor. Three pneumatic indicators, one for each air bushing, are used as a visual safety indicator in the pneumatic circuit. As long as the air bushing is provided with a certain air pressure, the red balloon in the indicator stays inflated, confirming the air bushings are pressurized. Figure 2b shows a picture of the WEC as installed at the AAU wave basin with indication of the location of the sensors.  The motor drive, the Speedgoat IO133 terminal board, DC power supply and loss current switches are centralized in the control cabinet. Figure 3 shows a picture of the inside of the control cabinet. During the experimental campaign, the VTI (VTI Instruments Corporation, Irvine, CA, USA) wave gauge sensor system was put inside the control cabinet to establish the connection of the VTI wave gauge sensor analog output with the Speedgoat IO133 terminal board analog input. The VTI wave paddle trigger is used for synchronization by providing a constant voltage signal from the moment the wave paddles are activited. The seismic accelerometer (PCB 393B04) is placed on top of the steel frame to quantify possible vibrations of this frame. Vibrations of the frame and resonance in particular are to be avoided, since these affect the measurement quality of the other sensors.
The torque request in the Simulink model is sent by the EtherCAT (Ethernet for Control Automation Technology) communication protocol to the Beckhoff motor drive as a Master Data Telegramm (MDT) process parameter. On the other hand, the Speedgoat target machine can receive by EtherCAT communication Amplifier Telegramm (AT) process parameters from the Beckhoff motor drive. The Beckhoff PMSM input and output signals are sent from and to the Beckhoff motor drive by the One Cable Technology (OCT), which allows to power the motor and process feedback. The motor drive receives the absolute position within one revolution at an 18 bit resolution from the single-turn absolute encoder. This encoder allows real-time determination of the state (position and velocity) of the WEC buoy. The drive provides the motor with a certain current, corresponding to a torque by multiplication with the torque constant of 1.91 Nm/A [16].
The uncertainty in the measured values with the given instrumentation is minimized by the selection of sensors with high resolution, high accuracy and low linearity error. It is important that the sensors are correctly calibrated, with a zero offset for the equilibrium position with draft 0.16 m, represented in Figure 2a. Besides the uncertainty related to sensor measurements, it is important to quantify the uncertainty of the physical testing results. Lamont-Kane et al. identified five distinct sources of uncertainty for physical testing of WEC-arrays [20]: Spatial variation of the wave-field within the wave basin; temporal variation of the wave-field from one repeat to another; the repeatability of model response for any single individual WEC; the reproducibility of model response between various nominally identical WECs (not applicable for the test campaign with a single, isolated WEC in the presented article) and the variation in the time-series of an incident irregular wave train. The quantification of these sources of uncertainty for the performed experimental campaign is not addressed in the presented article.

Wave Basin Setup
The experimental campaign took place at the AAU wave basin of the Ocean and Coastal Engineering Laboratory of the Department of the Built Environment [21]. The wave basin measures 14.60 m × 19.30 m × 1.50 m (length × width × depth) with an active test area of 8.00 m × 13.00 m (length × width). The wave generation system is 13 × 1.5 m (width × height) with 30 individually controlled wave paddles (snake type configuration). The system allows accurate generation of 3D waves due to narrow vertically hinged paddles (0.43 m segment width) with maximum wave height up to 0.45 m (at 3.0 s period) and typical maximum significant wave height H s in the range of 0.25-0.30 m. The wave basin is equipped with passive wave absorber elements. The AwaSys wave generation software is able to generate regular, irregular, solitary waves, execute 2-D and 3-D active wave absorption (reflection compensation) and generate 2nd order irregular unidirectional and multidirectional waves [22].
Seven resistive wave gauges (WGs) are installed in the wave basin to measure incident, diffracted, radiated and reflected waves during the tests. Figure 4a shows the planview layout of the wave basin with the central location of the WEC buoy, WGs, wave generation system and passive absorption. The interdistance between the different WGs and the WEC buoy are indicated in m and waves are generated from the bottom of Figure 4a. These seven WGs are also displayed in Figure 5a. The numbering is based on the used analog input ports of the Speedgoat IO133 terminal board. Figure 4b shows a 'bird's-eye' perspective of the experimental setup at the AAU wave basin. The WEC is attached to the bridge over the wave basin, complying with the spirit level requirements. Waves are generated from the top right corner in Figure 4b. The host PC for the Simulink control and the PC with the AwaSys software are located in line with the bridge, on the bottom right corner in Figure 4b. The water level is set equal to 1.010 m, which results in an equal positive and negative heave stroke with a magnitude of 0.25 m. This water level should be kept constant during the experimental campaign, since changing water level yields an offset for the position measurements from the laser sensor and motor encoder. A water density ρ of 1000 kg/m 3 is taken into account. Figure 5a shows a picture of the overview of the setup as schematized in Figure 3. Figure 5b shows the coordinate system and sign convention on a picture of the experimental setup, as adopted in Figures 1 and 2a. During this test campaign, waves are generated in the x-direction, corresponding to the direction of the greatest stiffness of the WEC.   Table 1 gives an overview of the different types of tests that are performed. Based on their underlying purpose, they are subdivided in three categories. The next sections discuss the different test categories.

Friction Model
The air bushings coefficient of friction is a function of air shear from motion, not from surface contact. Therefore, a friction coefficient of 0.00008 in log-scale due to the contribution from air molecules and gravitation can be taken into account [23]. To obtain the air gap between the shaft and the air bushing, a compressor with a two stage air filter supplies these bushings with clean and dry air under a nominal pressure of 5.5 bar. The gearbox break-away torque T 01 is expected to be around 0.20-0.40 Nm, with convergence towards lower values for a longer operational lifetime [24]. The motor static friction M R is reported to be equal to 0.02 Nm [15]. As a result, the total drivetrain friction experienced by the WEC will be mainly determined by the Coulomb friction attributed to the gearbox. Moreover, additional Coulomb and viscous friction attributed to the rack and pinion will occur. Based on the empirical characterization of the actual friction, a static friction model is constructed, discussed in Section 3.2.1 [25].

Friction Characterization Tests
The friction of the WEC with the motor and gearbox installed is characterized by a zero torque command on the motor. The WEC buoy is manually moved up and down from below the load cells with following elements in the test sequence: slowly at a targeted constant velocity close to 0.0 m/s and an accelerated motion with amplitude up to 0.40 m/s. Figure 6a shows a picture of a friction characterization test in the empty AAU wave basin.
The tests are repeated in the filled AAU wave basin to benefit from buoyancy forces. During the execution of the lifting procedure, the force acting on the loadcells F loadcells , the position z, the velocityż and the accelerationz are measured. The first element of the test sequence allows us to determine the Coulomb damping coefficient C Cou , while the second element allows us to determine the viscous damping coefficient C Vis .
The friction force F Friction to compensate for equals F loadcells reduced by the acceleration force F acc caused by m top and by the rotational inertia of the motor and the gearbox. To take this rotational inertia of the motor and the gearbox into account in F acc , it is expressed as a mass in the heave direction, i.e., the mass attributed to the motor and gearbox inertia m MGI : with J motor the motor inertia equal to 6.17 kg cm 2 [16] and J gearbox the gearbox inertia equal to 0.71 kg cm 2 [26]. It is anticipated that the actual inertia will be higher, as the inertia of the rack and pinion is not taken into account in Equation (1). The resulting F Friction is: In case the motor and the gearbox are not installed and the air bushings work properly, F acc should equal F loadcells , resulting in F Friction = 0 N. Figure 7a shows the z andż time series for the friction characterization Test_054. Figure 7b shows the time series of F loadcells , F acc , F Friction and the fitting of C Vis and C Cou . F loadcells exceeds F acc due to the addition of the drivetrain (motor, gearbox and rack and pinion) friction. Fitting a model to F Friction results in the friction compensation model F Comp , given by Equation (3).
Note that an additional minus sign is necessary to have F Comp in the same direction asż. The model is composed of a Coulomb part proportional to the sign of the velocity and a viscous part proportional to the velocity. A velocity boundaryż Bou is imposed to avoid rapid sign switches for the Coulomb friction term at low velocities. Consequently, no friction compensation is taken into account between −ż Bou andż Bou . These boundaries correspond approximately with the noise range on the velocity feedback from the motor encoder. Moreover, it is preferred to partly compensate for the friction to preserve a realistic PTO, expressed by the compensation factor C C . Figure 8a shows the F loadcells toż mapping for the friction characterization Test_054, which allows us to fit the parameters of F Comp . Table 2 provides the obtained values for C Vis , C Cou , C C , andż Bou .   Figure 7b shows the time series of F Comp . The proper functioning of the velocity boundaries can be observed as less sign switches occuring for F Comp (yellow curve) compared to the fitting of C Vis and C Cou (magenta curve). The difference between F Friction (green curve) and F Comp (yellow curve) is the remaining PTO friction.
The friction characterization tests are repeated with different lateral loading conditions, as an approximation for the surge or sway wave excitation force. In the wave propagation direction (x-direction) −20 kg is applied with a tension spring and rope around the WEC buoy, shown in Figure 6b (Test_057). This loading of 196 N corresponds approximately with a surge wave excitation force for a wave with period T = 1.0 s and wave height H = 0.40 m, according to linear potential flow simulations with the open-source software package openWEC [27], with the integration of the Boundary Element Method (BEM) code Nemoh. The F loadcells toż mapping in Figure 8b shows no increased F loadcells for Test_057 compared to Test_054, demonstrating the proper functioning of the air bushings. In the Simulink model, the uncompensated input motor torque τ motor,uncomp augmented with F Comp (R pinion /i) yields the compensated input motor torque τ motor,comp : Note that the friction compensation in Equation (4) is implemented as a feedback control structure, since the velocity output of the WEC is used. Commonly, a friction compensation feedforward control structure is adopted when tracking a reference is the objective, which is not the case here [28]. Table 3 provides an overview of all the executed friction characterization tests, with indication of the filling of the AAU wave basin, lateral loading force, lateral loading direction and applied air pressure. The Test_ID is the identification number of each performed test. To perform a SID of the WEC, a linear decomposition model formulation for the wave-WEC interaction is adopted. The problem is separated into radiation and excitation components. Distinct tests are performed for both the radiation and the excitation components to determine the WEC system response, whereafter radiation and excitation frequency response functions (FRFs) can be constructed from the discrete frequency components. At a high level of abstraction, the WEC can be considered as a system with two inputs: the surface elevation η and F PTO , as displayed in Figure 9. When the WEC system, consisting of the PTO and the WEC buoy, is assumed to be linear for small motions and small waves, superposition can be applied. The WEC system is considered to be composed of two single input single ouput (SISO) models, the radiation and the excitation model, as displayed in Figure 9 [12,13,29]: The quantities are expressed as a function of the angular frequency ω. The upper case indicates that these variables are all in the frequency-domain, while the hat symbol denotes that these variables are complex quantities. The control parameters for the impedance matching controller are determined from the intrinsic impedance Z i (ω), defined as the inverse of G(ω): Mechanical impedance is a measure of the opposition to motion from a source when a potential is applied, defined as the ratio of force (potential) to velocity (flow) in Equation (6). Note that Ẋ is used instead of Ż for the FRF ofż, to avoid confusion with the defined Z i . The excitation model is obtained by computing the ratio of the FRF of the output excitation force F e , measured as F loadcells , to the FRF of the input η at the location of the WEC, resulting in the excitation force coefficients H e (ω): In the first instance, the above-discussed SID tests are carried out in "open loop", which means no output feedback is considered, represented in Figure 9 by no feedback arrow. A "closed loop" is obtained at any time the WEC is controlled and F PTO is calculated based on an output measurement. In Figure 9, F PTO depends onż, as F PTO = Cż, with C representing the control system dynamics. In case of resistive damping, C is a negative constant in this formulation. As a result, the frequency-domain WEC equation of motion, as displayed in Figure 9 as the dual SISO model, is given by: By applying the superposition principle, the decomposition in the radiation and excitation model is given by:

Impedance Formulation and Radiation Tests
When the Fourier transform F differentiation property is applied to writez in terms ofż: the point absorber WEC equation of motion can be written in terms of Ẋ (ω) [12,30]: where m = 36.83 kg, A(ω) is the added mass coefficient and B(ω) is the hydrodynamic damping coefficient. The hydrostatic stiffness coefficient K is given by: where S is the cross-sectional area of the WEC buoy at the SWL equal to 0.283 m 2 , ρ = 1000 kg/m 3 and g the gravitational acceleration equal to 9.81 N/kg. Rearrangement of Equation (11) based on Equation (9), results in Z i (ω): Z i (ω) is experimentally determined by executing a forced oscillation test in the AAU wave basin, without waves generated by the wave paddles. In Equation (13), F e equals zero and F PTO is a chirp signal with a frequency spectrum covering the bandwidth of interest [13]. Since it is experimentally more convenient to use F PTO as an input and measure the outputż, Z i is obtained as the inverse of G, defined by Equations (5) and (6). Table 4 shows an overview of the performed radiation tests. The chirp-up PTO input torque is defined with an initial frequency of 0.0 Hz, a target time of 220 s and a frequency at target time of 4.0 Hz. The chirp-down PTO input torque is defined with an initial frequency of 4.0 Hz, a target time of 220 s and a frequency at target time of 0.0 Hz. The noise signal is defined with an amplitude of 1.0 Nm, multiplied by the defined PTO gain to control the absolute maximum heave amplitude of the WEC buoy z max . Most of the test were executed with C C = 0.0, without the friction compensation model discussed in Section 3.2 included. Test_112 and Test_113 consider C C = 0.6, which is the WEC system to build the resistive and reactive controller on.  The calculated Z i can be displayed as a bode plot with gain and phase, given in Figure 11 for Test_021, Test_059, Test_066, Test_112 and Test_113 [12]. Comparing Test_021 with Test_059 confirms that the chirp-up and chirp-down signal yields the same Z i (ω). Test_066 with a gain of 0.5 and z max = 0.041 m results in a higher identified Z i compared to Test_059 with a gain of 1.0 and z max = 0.116 m, stressing the importance of covering motion amplitudes representative for the WEC during operation. Test_112 and Test_113 demonstrate how the implementation of F Comp alters the WEC system dynamics. As expected,ż is more amplified when the friction is compensated. The WEC resonance frequency f n is equal at 0.84 Hz, corresponding to a natural period T n = 1.19 s.

Excitation Tests
H e (ω) is experimentally determined by locking the WEC in equilibrium position, as defined in Figure 2a, imposing a frequency rich η signal and measuring F loadcells . The WEC is fixed by deactivating the motor drive, resulting in an active holding brake. Application of Equation (7) results in H e (ω). Table 5 shows an overview of the performed excitation tests. The wave input generated by the wave generation system is a JONSWAP wave spectrum, defined by a peak enhancement factor γ of 3.3 and the mentioned significant wave height H s and peak period T p . In the case of Test_133 with regular waves, H and T are mentioned. Using a start signal from the wave paddles and recording the wave paddle motion, allow for a deterministic comparison of the same sea state over different control inputs and types.  Figure 12a displays an excitation test for the WECfarm WEC, where the incoming and diffracted waves can be observed. Figure 12b shows the F loadcells and η WG11 time series for Test_074. When a wave crest (negative η) passes the WEC buoy, the load cells are compressed, resulting in a positive F loadcells . When a wave trough (positive η) passes the WEC buoy, the load cells are under tension, resulting in a negative F loadcells . Figure 13a shows the resulting F loadcells , equivalent to F e , and η for Test_074 and Test_077. The calculated H e can be displayed with a gain and phase, given in Figure 13b for Test_074 and Test_077. Frequency smoothing has been performed on F loadcells , η and the gain and phase of H e . The Gaussian-weighted moving average over a window of 30 frequency intervals of 0.0042 Hz has been taken. Figure 13b confirms the higher H e for lower wave frequencies, as noticed in Figure 13a. To calculate H e (ω) based on the most straightforward conceptual definition, a WG is placed on the location of the WEC buoy when the WEC buoy is not present and then record the η time series. The WEC buoy is then put in place and the same wave time series are run again, this time measuring the force on the WEC buoy caused by the waves. For control purposes, this procedure is of little use, since it is clearly not possible to measure η at the point where the WEC is located, once the WEC is in place. Therefore, η obtained with WG 11 is used as an approximation for η on the location of the WEC buoy. It is shown by Bacelli et al. that if the distance between the WEC buoy and the WG is increased, the term describing the diffracted waves at the WG becomes small enough and the original model given by Equation (9) can be adopted [13].
To verify the assumption that the data of WG 11 of these excitation tests could be used for this purpose, at the end of the test campaign some tests with wave spectra are executed with the WEC buoy removed and WG11 moved to the position of the WEC buoy, equilinear with the other WGs. For these tests, the interdistance WG 6 to WG 11 is 1.31 m and the interdistance WG 11 to WG 8 is 0.98 m. Figure 14 shows a picture of this layout of the WGs. Three JONSWAP spectra and one regular wave are tested, defined in Table 6. The top Figure 15a shows the time series comparison between η measured by WG11 in excitation Test_074 and η measured by WG11 in Test_172. The bottom Figure 15a shows this comparison for Test_050 and Test_173. For both wave conditions, a minor difference between the two time series can be observed. Figure 15b shows η for Test_074 compared to Test_172 and η for Test_050 compared to Test_173. Figure 15a,b confirm that η obtained with WG 11, according to the layout of the WGs in Figure 4a, can be used as an approximation for η on the location of the WEC buoy.

Free Decay Tests
Free decay tests are executed to determine the decay response. The heave displacement z allows us to calculate the logarithmic decrement Λ, the corresponding damping ratio ξ d and the natural period T n . This T n should correspond with the one obtained from the radiation tests discussed in Section 3.3.2. Conceiving the WEC as a mass-spring-damper system, ξ d follows from Λ [31]: in which x 1 and x n are the values of the first and the n-th peak of the free decay motion z. For un underdamped system (0 < ξ d < 1), the oscillations of the WEC buoy fade exponentially over time and tend toward zero, yielding the envelope: in which z A is the amplitude of the free WEC response [31]. The angular natural frequency ω n of the WEC buoy can be calculated from the damped natural frequency ω d : Table 7 gives an overview of the executed free decay tests, with indication of the start position of the WEC buoy z start , C C , n used Equation (14), ξ d and T n . A positive z start corresponds to a submerged start position and a negative value corresponds to an elevated start position. The WEC buoy is submerged or elevated with a torque command. Once the torque is set equal to zero, the decay motion is initiated. Test_089 and Test_129 yield a T n equal to 1.19 s, corresponding to the value obtained by the bode plots of the radiation tests displayed in Figure 11. When the motor and gearbox with pinion are removed in Test_170, a T n equal to 1.00 s is obtained. This lower T n can be attributed to the removed inertia within the PTO drivetrain since removing inertia results in a lower T n . Test_089 and Test_129 yield a ξ d of 0.13 and 0.12, respectively, corresponding to an underdamped system. Given this low ξ d , ω n will closely approximate ω d . The limited draft in combination with the flat bottom will introduce important bottom slamming effects for tests where the WEC buoy re-enters the water after being lifted out [32]. These bottom slamming forces are assessed in Test_161 and Test_162, for which the start position of the WEC buoy is completely lifted out of the water. In addition, nonlinear viscous drag effects occur. The semi-empirical Morison equation describes the viscous force F vis in function of a viscous drag coefficient C D [33]: This C D can be identified from experiments or from fully viscous modeling methods based on the Navier-Stokes equations [34][35][36]. The experimental quantification of bottom slamming effects and nonlinear viscous drag effects is not addressed in the presented article.
Test_170 considers multiple free decay tests for the WEC buoy with the motor and gearbox with pinion removed. While z is obtained by both the laser and encoder for the other tests, only laser measurements are available when the motor is removed. The WEC buoy is manually pushed down or lifted from the measurement bridge. Figure 16a displays the WEC buoy displacing a water volume, equivalent to the volume of the WEC buoy, after being dropped from −0.202 m in Test_161. The concentric circular radiated waves can be observed. Bottom slamming forces are obtained from F loadcells . Figure 16b shows the free decay z time series and corresponding exponential decay envelopes according to Equation (15) for Test_089, Test_129 and Test_170. Test_089 with an elevated start position of -0.089 m results in a quasi-equivalent decay response as Test_078 with a submerged start position of 0.078 m. In order to allow a comparison of the symmetry of the WEC system, Test_089 is plotted with inversed sign. Figure 17 displays the four quadrants where the WEC PTO system acts as a motor or generator, according to the adopted sign convention [37]. When F PTO and the resultantż have an opposite sign, energy will be extracted from the waves and the PTO acts as a generator, occurring in quadrant II and IV. When F PTO and the resultantż have an identical sign, energy will be consumed and the PTO acts as a motor, occurring in quadrant I and III. The defined convention in the Simulink model is that a net positive absorbed power value P abs corresponds to mechanical power absorption and a net negative P abs corresponds to adding power to the system. Therefore, a minus sign is added to the multiplication of F PTO withż in order to obtain P abs : P abs = −F PTO ·ż (18) Note that this calculation results in instantaneous values. For regular waves, averaging over a number n of wave periods T results in the averaged absorbed power value P abs :

Power Absorption
3.4.1. Impedance Matching Impedance matching can be applied to enable maximum power transfer between two oscillatory systems, the ocean waves and the WEC PTO system. Since impedance is defined as a complex value, the PTO generally has a resistance component (real part) and a reactance component (imaginary part). The maximum power transfer theorem says that the maximum possible power is delivered to the PTO when the PTO impedance Z PTO (load impedance or input impedance) is equal to the complex conjugate (represented by *) of the impedance of the source Z i (intrinsic impedance or output impedance). For two impedances to be complex conjugates their resistances must be equal and their reactances must be equal in magnitude and opposite in sign. With the WEC impedance model given by Equation (13), an optimal PTO is obtained by Z PTO (ω) [12,38]: The frequency response of the Proportional Integral Derivative (PID) controller is given by [12]: K P , K I and K D are the proportional, integral and derivative gains of the controller, respectively. The rearrangement of Equation (21) results in the PID controller impedance: Given Equations (20) and (22), the proportional gain for the Proportional (P) controller, equivalent to the damping coefficient C PTO,P for the resistive control strategy, is given by: Given Equations (20) and (22), the proportional and integral gain for the Proportional Integral (PI) controller, equivalent to C PTO,PI and the spring coefficient K PTO,PI for the reactive control strategy, are given by: The presented study considers only P and PI control. Gu et al. present the implementation of a PID controller in a frequency domain model for a heaving point absorber WEC [39]. In this case the derivative controller acts on the acceleration term, corresponding to mass control.

Causal Impedance Matching P and PI Controller
For any causal system, the complex conjugate of the impedance will produce an acausal system. Since a WEC is a causal system, an acausal system is obtained for Z * i (ω). Therefore, to implement Equation (23) or Equations (24) and (25) perfectly across all wave frequencies, the future velocity of the WEC has to be known, resulting in an acausal controller implementation [12]. However, the controller does not need to operate perfectly at all wave frequencies simultaneously and can be designed to work well in a restricted range of frequencies. The energy in a typical sea state ranges over at most a single decade of frequencies, with the sea state changing appreciably only over the course of hours. This band-limited and slowly varying nature of ocean waves allows to utilize a causal realization of the impedance matching approach. Approximating Z * i (ω) in the peak frequency ω p of the design sea state is straightforward to implement in the Simulink model. The impedance of the causal impedance matching P controller is equal to C PTO defined in Equation (24) [12,38]: The impedance of the causal impedance matching PI controller is an interpolation of the impedance given by Equation (22) in ω p [12,38]: Based on the experimentally determined Z i in Test_112, the impedance of a causal impedance matching P controller and PI controller can be determined, according to Equations (26) and (27), respectively. Figure 18a shows the gain and phase of these controllers, for an interpolation point of T = 1.50 s. Since in the power absorption tests various wave conditions are considered, Figure 18b plots C PTO,P (ω) (Equation (23)), C PTO,PI (ω) (Equation (24)) and K PTO,PI (ω) (Equation (25)) in function of the wave frequency. The resulting control parameters for T = 1.50 s and T = 2.00 s are indicated. As expected, K PTO,PI (ω) = 0 kg/s 2 for T n = 1.19 s.

Resistive Control
The baseline control strategy is resistive control, equivalent to causal impedance matching P control as defined in Equation (26). F PTO,P is equal toż multiplied with a positive C PTO : A minus sign is added, since F PTO,P should opposeż. Tables 8 and 9 provide an overview of the performed resistive control tests, for regular and irregular waves, respectively. The applied C C is given. For the regular waves, characterized by T and H, a range of C PTO is tested during a single test. For the irregular waves, characterized by a JONSWAP wave spectrum with γ = 3.3, H s and T p , a fixed C PTO is applied.   Figure 19a shows a snapshot of the z andż time series for Test_100. Figure 19b shows the corresponding C PTO input, with a stepwise increase from 400 kg/s to 450 kg/s, and the P abs output.

Reactive Control
Reactive control, equivalent to causal impedance matching PI control as defined in Equation (27), aims to bring the WEC into resonance by the addition of a spring. F PTO,PI is composed ofż multiplied with a positive C PTO and z multiplied with negative K PTO .
Again, a minus sign is added, since F PTO,PI should opposeż. Table 10 gives an overview of the performed reactive control tests for regular waves, characterized by H and T. C C is fixed at 0.6 and a range of C PTO and K PTO is tested during each test.  Figure 20a shows a snapshot of the z andż time series for Test_138. The Simulink control model did not impose constraints on z,ż,z, allowing these high values for z andż. The stroke is mechanically limited to 0.25 m by the micro switches, displayed in Figure 3. Figure 20b shows the corresponding C PTO input, with a stepwise increase from 60 kg/s to 72 kg/s, and the P abs output. The net P abs is calculated as the period averaged difference between the integrated generated power (positive P abs ) and the integrated added power (negative P abs ), according to Equation (19) and the sign convention as defined in Figure 17. K PTO is equal to −1520 kg/s 2 for the displayed time window.

Power Absorption Comparison between the Resistive and Reactive Controller
Two regular wave conditions, characterized by H = 0.09 m, T = 2.00 s and by H = 0.09 m, T = 1.50 s, respectively, are considered to compare the presented resistive and reactive control strategy. P abs is calculated with n = 5 in Equation (19). For resistive control, C PTO,P is obtained according to Equation (23) and for reactive control, C PTO,PI and K PTO,PI are obtained according to Equations (24) and (25), respectively. Figure 18b displays the obtained control parameters for T = 2.00 s and T = 1.50 s. For C PTO,PI , lower values than experimentally indentified are adopted. For the regular wave with H = 0.09 m and T = 2.00 s, resistive control (Test_100) with C PTO = 450 kg/s yields P abs = 1.89 W and reactive control (Test_138) with C PTO = 72 kg/s and K PTO = −1520 kg/s yields P abs = 10.26 W. For the regular wave with H = 0.09 m and T = 1.50 s, resistive control (Test_103) with C PTO = 200 kg/s yields P abs = 3.09 W and reactive control (Test_147) with C PTO = 70 kg/s and K PTO = −1120 kg/s yields P abs = 11.75 W. Figure 21 visualizes these results. For the reactive controller, a significant increase in power absorption compared to the resistive controller is observed, which is according to the literature [40]. Figure 20a shows that resonance occurs as z max = 0.165 m is significantly higher than H/2 = 0.045 m. To properly assess the power performance of the reactive controller compared to the resistive controller, irregular waves should be considered, as discussed in Section 3.4.2.
The amplification on z requires the motor to add energy, which has a certain efficiency. Apart from larger peaks in F PTO and P abs , a major drawback of reactive control is the energy loss by dissipative processes inherent to the back-and-forth energy exchange between the PTO and the WEC buoy, especially when the magnitude of the exchanged energy is comparable to, or even significantly larger than, the net absorbed energy [41]. The presented research does not consider the PTO efficiency, nor the damping force, nor the reactive force. Strager et al. present a method to determine the optimal reactive control parameters for a given combination of non-ideal PTO efficiency and monochromatic wave frequency [42].

Discussion and Conclusions
Within the 'WECfarm' project, two test campaigns are performed at the AAU wave basin: (a) a testing of the first WEC in April 2021 and (b) a testing of a two-WEC array in February 2022, in preparation of five-WEC array tests. The main objective of the 'WECfarm' project is to cover the scientific gap on experimental data necessary for the validation of recently developed (non-linear) numerical models. The presented article discusses experimental testing of an isolated WEC, being the first test campaign within the 'WECfarm' project. The primary objective of evaluating the hydrodynamics, electromechanics, control platform, DAQ and structural performance of the WEC is to allow extending the setup to a five-WEC array.
The WEC buoy is a truncated cylinder with a high diameter to draft ratio to increase radiation and WEC-WEC interactions. The PTO system of the WEC is a PMSM connected to a gearbox powering a rack and pinion system. A configuration of three air bushings guarantees a permanent layer of air between the guide shafts and the bushings, resulting in zero-friction linear guiding. The WEC control and data acquisition are realized with a Speedgoat Performance real-time target machine, offering the possibility to implement advanced WEC array control strategies in the MATLAB-Simulink model. This unique PTO system in combination with the real-time target machine makes active and accurate PTO control possible.
The drivetrain (motor, gearbox, rack and pinion) friction is assessed by manually moving the WEC buoy up and down, from below the loadcells and under a zero torque command. Based on the empirical relationship between F loadcells andż, a simplified model based on Coulomb and viscous friction is determined. The resulting friction compensation model F Comp is implemented in Simulink as a torque augmentation on the torque command.
A SID approach is adopted, considering the WEC system to be composed of two SISO models, the radiation and the excitation model. The radiation tests in calm water, with F PTO chirp-up and chirp-down noise signals as input andż as output, yield the intrinsic impedance Z i . The excitation tests with the WEC buoy fixed, with as input various JONSWAP spectra and as output F loadcells , yield the wave excitation force coefficients H e (ω).
The assumption that WG measurements at sufficient distance from the WEC buoy may be used to characterize η on the location of the WEC is confirmed. Free decay tests characterize the WEC buoy decay response, confirming that T n = 1.19 s, obtained by the radiation tests.
Power absorption tests are executed with the resistive and reactive control strategy. The tests with regular waves are executed with 'real-time' tuning of the control parameters in the Simulink model. Adopting an impedance matching approach, the optimal C PTO,P for resistive control and the optimal C PTO,PI and K PTO,PI for reactive control are calculated from Z i . For the two selected regular wave conditions, characterized by H = 0.09 m, T = 2.00 s and by H = 0.09 m, T = 1.50 s, the reactive controller results in a significant higher averaged absorbed power P abs compared to the resistive controller. However, the presented research does not consider PTO efficiency, which is detrimental for reactive control.
The implementation of the friction compensation model F Comp proved to be a good methodology to partly compensate Coulomb and viscous friction attributed to the drivetrain. The experimentally determined radiation and excitation model yields a simple, though accurate, model of the WEC system. The intrinsic impedance, resulting from the radiation test, is used to design a causal impedance matching P and PI controller, equivalent to a resistive and reactive controller. Although this approach was only tested limitedly, mainly for regular waves, the successful extension to a multiple-WEC array, considering irregular waves, is confirmed. The testing of the isolated 'WECfarm' WEC proved to be successful and the extension of the setup to a five-WEC array will allow it to comply with the future research objectives.  Data Availability Statement: The datasets resulting from the WECfarm project will be made available in due time on https://www.awww.ugent.be, accessed on 23 August 2022.

Abbreviations
The following abbreviations are used in this manuscript: