Multi-Timescale Lookup Table Based Maximum Power Point Tracking of an Inverse-Pendulum Wave Energy Converter: Power Assessments and Sensitivity Study

: A novel, inverse-pendulum wave energy converter (NIPWEC) is a device that can achieve natural period control via a mass-position-adjusting mechanism and a moveable internal mass. Although the energy capture capacity of a NIPWEC has already been proven, it is still meaningful to research how to effectively control the NIPWEC in real time for maximum wave energy absorption in irregular waves. This paper proposes a multi-timescale lookup table based maximum power point tracking (MLTB MPPT) strategy for the NIPWEC. The MLTB MPPT strategy was implemented to achieve a theoretical “optimal phase” and “optimal amplitude” by adjusting both the position of the internal mass and linear power take-off (PTO) damping. It consists of two core parts, i.e., internal mass position adjustment based on a 1D resonance position table and PTO damping tuning based on a 2D optimal PTO damping table. Furthermore, power assessments and sensitivity study were conducted for eight irregular-wave sea states with diverse wave spectra. The results show that energy period resonance and the lookup table based PTO damping tuning have the highest possibility of obtaining the maximum mean time-averaged absorbed power. Additionally, both of them are robust to parameter variations. In the next step, the tracking performance of the MLTB MPPT strategy in terms of changing sea states will be studied in-depth.


Introduction
Wave energy is a carbon-free renewable resource that possesses the advantages of a high energy density and has outstanding availability and better predictability than wind or solar energy [1].Besides, the global wave power is about 29,500 TWh/year, according to the estimation from Ocean Energy Systems [2].If wave energy is fully utilized, the demand for global electricity consumption can be sufficiently satisfied.The data given by the International Energy Agency show that the global electricity consumption in 2019 reached 22,848 TWh [3].Therefore, wave energy has enormous development potential.
However, there still exist several technical and non-technical challenges for wave energy utilization (WEU) [1,4].The wave-by-wave, hour-by-hour, and site-by-site variations of the level of wave power represent one of the main obstacles restricting technological convergence and maturity.Therefore, many researchers focus on how to improve the performance of wave energy converters (WECs), especially wave conversion efficiency, under complex sea states [5].
WEC operation control has always been deemed a research hot spot since a successful control method can significantly enhance the wave conversion efficiency of an existing WEC with low additional investment [1].To date, plenty of control methods have been proposed based on different actuators and control algorithms, wherein the phase control methods (PCMs), natural period control methods (NPCMs), and maximum power point tracking (MPPT) are related to the research content of this paper.Hence, the aforementioned three control methods are introduced respectively.PCMs, e.g., latching control (LC) and declutching control (DeC), were first presented between the end of the 1970s and the beginning of this century.PCMs are usually based on PTO force tuning, which can compulsively synchronize WEC movement with the variations in wave excitation force to achieve an "optimal phase".The theoretical basis of PCMs was summarized by Falnes [6] as a universal theory on ocean waves and oscillating systems.Although the theory is linear, it qualitatively points out the core idea of most control methods, not only PCMs.Herein, the core idea means that the maximum absorbed power will appear when both "optimal phase" and "optimal amplitude" are achieved."Optimal phase" denotes a resonant state where the WEC velocity and wave excitation force are in phase.Meanwhile, "optimal amplitude" refers to damping matching; e.g., the linear damping of a power take-off (PTO) should be equal to radiation damping in regular waves.Nowadays, LC or DeC, with or without the knowledge of future excitation force, still attracts many researchers [7,8].
NPCMs are another solution for achieving "optimal phase".They alter the natural period/frequency of a WEC to be identical to the period/frequency of an incident wave.Three ways can be used to realize natural period control.The first is inertia adjustment [9,10], the second can be called stiffness adjustment (SA) [11][12][13][14][15][16], and the third is known as surface control [17].SA is the most common NPCM.It refers to adjusting the inherent stiffness term in a WEC.Marei et al. [11] changed compressed gas stiffness to achieve the resonance of an Archimedes wave swing.Additionally, Chen et al. [12] implemented the SA of a built-in horizontal axis pendulum to achieve resonance.Moreover, Refs.[13][14][15][16] pointed out that the hydrostatic restoring stiffness of an oscillating wave surge converter (OWSC), e.g., a solo duck WEC or an inverse-pendulum-type WEC, can be effectively regulated for resonance by using a moveable mass.
When compared to PCMs, NPCMs are mainly based on slow tuning processes.It means that energy capture performance can be effectively improved as long as a WEC adjusts its natural period to the wave characteristic period of an irregular wave within a large duration.In addition, NPCMs do not need to tune PTO force at a wide range, which reduces the control difficulty and improves control reliability.However, NPCMs are only able to realize "optimal phase".If "optimal amplitude" is required, PTO force tuning should be inevitably adopted.Therefore, how to comprehensively design and utilize two control mechanisms to achieve both "optimal phase" and "optimal amplitude" could become an issue worthy of in-depth research when applying NPCMs.
MPPT is a real-time search technology for the maximum power points of clean energy power, such as wind power, photovoltaics, thermoelectricity [18], etc.In 2009, Amon et al. [19] first introduced the concept of MPPT to WEU, conducting a perturbation and observation (P&O) algorithm on the duty cycle of a Buck circuit as well as the simulated resistance of a three-phase pulse width modulation (PWM) rectifier.Nowadays, the MPPT algorithm family includes a fixed-step P&O [19][20][21][22], a variable-step P&O [23,24], a segmental fixed-step P&O [25], a lookup-table-based (LTB) MPPT [26,27], diverse heuristic algorithms [28,29], etc.The fixed-step P&O is the most common algorithm and often serves as a reference for other algorithms.Besides, LTB MPPT is an ordinary but practical algorithm that has been applied to an oscillating water column-type WEC and a floatingpendulum WEC [26,27].Yue et al. [27] proved that a two-dimensional (2D) LTB MPPT algorithm, which contains an optimal duty cycle table derived from regular wave simulations, is a simple, reliable, and excellent solution for WEC output power enhancement in irregular waves.Refs.[30,31] constructed an emulator for a point absorber WEC (PAWEC) and reviewed the MPPT techniques or hybrid strategies of PAWECs.
Energies 2023, 16, 6195 3 of 25 This paper focuses on a novel inverse-pendulum wave energy converter (NIPWEC), which can conduct SA via a mass-position-adjusting mechanism (MPAM) and a moveable internal mass.The structural details of a NIPWEC can be seen in Section 2. Similar structures can be seen in Refs.[14][15][16]32].Although the energy capture capacity of a NIPWEC has already been proven via simulations in a numerical wave tank, the prototype tests in a physical wave flume and the annual average power assessments of a 1:16 fullscale prototype [14,15] for how to effectively control the NIPWEC in real-time is still a meaningful research topic.The inefficient implementation of a control method can lead to ordinary or even weakened energy capture performance.Hence, this paper researches an implementation strategy for SA and PTO force tuning, aiming to find a suitable control algorithm for NIPWECs for real irregular waves.
Dong et al. [32] proposed a frequency and amplitude control (FAC) strategy to realize the MPPT of a NIPWEC.Herein, frequency control means adjusting the internal mass position according to the peak frequency of a wave spectrum.Meanwhile, amplitude control refers to tuning the PTO force to find the inherent impendence of the NIPWEC, which is approximately calculated by a weighted average method.The premise of FAC is the fast Fourier transformation (FFT) and frequency-domain analysis of the wave excitation moment within the past hundreds of seconds at least.This may lead to the following two issues.
(1) The nonlinear factors, e.g., an endstop moment for the swing amplitude limitation and a sine function hydrostatic restoring moment, are ignored by frequency domain analysis.This may result in the deterioration of MPPT performance when the ideal NIPWEC inherent impendence is applied to PTO force tuning.
(2) The real-time tracking performance may be weakened since both frequency control and amplitude control are implemented every few hundred seconds.
In contrast to Ref. [32], this paper presents a simple and reliable control algorithm, i.e., the muti-timescale lookup table based maximum power point tracking (MLTB MPPT), for NIPWECs.This algorithm is constructed and conducted based on a similar idea that can be found in Ref. [27].MLTB MPPT consists of LTB internal-mass position adjustment and LTB PTO damping tuning, which are obtained via time domain analysis with consideration given to the above nonlinear factors.Additionally, MLTB MPPT is capable of tuning the PTO force wave by wave, which can reinforce its real-time tracking performance.Detailed descriptions of MLTB MPPT are illustrated in Section 2. Furthermore, the power assessments and sensitivity study of MLTB MPPT are comprehensively investigated.
The rest of this paper Is organized as follows.First, the structure and implementation procedure of MLTB MPPT is elaborated in Section 2. Next, a NIPWEC dynamic model is described in Section 3.Then, the parameter settings for the time domain simulations are displayed in Section 4. Afterwards, the maximum-power-point (MPP) analysis in regular waves, as well as the power assessments and sensitivity study for MLTB MPPT in irregular waves, is discussed in Section 5. Last, the conclusions, along with future work discussion, are revealed in the Section 6.

Three MPPT Implementation Methods for a NIPWEC
We focus on a NIPWEC with a mechanical power take-off.As shown in Figure 1, the NIPWEC is constrained on a fixed supporting shaft and swings along it driven by ocean waves.Next, a speed-increase mechanism accelerates the inverse-pendulum swing and converts it into the fast bidirectional rotation of a permanent magnet synchronous generator (PMSG).Then, the PMSG generates continuous but unstable three-phase alternating-current (AC) electricity.Afterwards, a three-phase PWM rectifier transforms the AC electricity into direct current (DC) and regulates the PTO force via the vector control of a PMSG.Last, the DC bus smooths DC voltage and balances power generation and consumption via an energy storage unit and diverse electric loads.The NIPWEC can automatically tune its natural period within a wide range by adjusting the vertical position of an internal mass, according to the changing characteristic period of irregular waves.Herein, the mass position-adjusting mechanism (MPAM) with a screw and a motor is adopted.The motor drives the screw to rotate forward/reverse in order to realize the up/down movement of the internal mass.Therefore, the NIPWEC is able to generate electricity at the resonant state for a long time, which leads to the effective improvement of the energy capture/generating capacity.
Falnes [6] pointed out that the single DoF motion in ocean waves can be assumed to be a forced linear oscillation system.Under this assumption, the maximum absorbed power can be obtained under the "optimal phase" condition, which means that the oscillation velocity and excitation force are in phase.However, Falnes' theory relies on a frequency domain model and ignores the influence of nonlinear factors, e.g., an endstop moment for the swing amplitude limitation and a sine function hydrostatic restoring moment [33].The two nonlinear factors are able to significantly affect the magnitude of the absorbed power under the "resonance" condition.This means that due to the existence of nonlinear factors, the optimal PTO damping max PTO R is not exactly equal to the radiation damping ω 55 ( ) r R at resonance wave frequency ω r .In addition, the actual ocean waves are irregular waves, for which the wave heights and wave periods change all the time.Therefore, it is meaningful to search for a simple and reliable MPPT algorithm based on a time domain model.
The basic goal of an MPPT algorithm is to achieve "optimal phase", i.e., resonance, and "optimal amplitude", i.e., optimal PTO damping, against the maximum absorbed power.For a NIPWEC, the above goal can be theoretically achieved via the following three implementation methods.Method 1 (M1): adjust the internal mass position X to achieve "optimal phase" and tune PTO damping PTO R to achieve "optimal amplitude".Herein, The NIPWEC can automatically tune its natural period within a wide range by adjusting the vertical position of an internal mass, according to the changing characteristic period of irregular waves.Herein, the mass position-adjusting mechanism (MPAM) with a screw and a motor is adopted.The motor drives the screw to rotate forward/reverse in order to realize the up/down movement of the internal mass.Therefore, the NIPWEC is able to generate electricity at the resonant state for a long time, which leads to the effective improvement of the energy capture/generating capacity.
Falnes [6] pointed out that the single DoF motion in ocean waves can be assumed to be a forced linear oscillation system.Under this assumption, the maximum absorbed power can be obtained under the "optimal phase" condition, which means that the oscillation velocity and excitation force are in phase.However, Falnes' theory relies on a frequency domain model and ignores the influence of nonlinear factors, e.g., an endstop moment for the swing amplitude limitation and a sine function hydrostatic restoring moment [33].The two nonlinear factors are able to significantly affect the magnitude of the absorbed power under the "resonance" condition.This means that due to the existence of nonlinear factors, the optimal PTO damping R PTOmax is not exactly equal to the radiation damping R 55 (ω r ) at resonance wave frequency ω r .In addition, the actual ocean waves are irregular waves, for which the wave heights and wave periods change all the time.Therefore, it is meaningful to search for a simple and reliable MPPT algorithm based on a time domain model.
The basic goal of an MPPT algorithm is to achieve "optimal phase", i.e., resonance, and "optimal amplitude", i.e., optimal PTO damping, against the maximum absorbed power.For a NIPWEC, the above goal can be theoretically achieved via the following three implementation methods.
Method 1 (M1): adjust the internal mass position X to achieve "optimal phase" and tune PTO damping R PTO to achieve "optimal amplitude".Herein, X = l M (t) − l M0 , where l M0 means the lower bound for l M (t), and l M (t) signifies the distance between the rotation center R and the adjustable internal mass center M. The PTO moment M PTO for M1 can be given by Energies 2023, 16, 6195 5 of 25 Method 2 (M2): adjust the PTO stiffness K PTO to achieve "optimal phase" and tune R PTO to achieve "optimal amplitude".M PTO for M2 is Method 3 (M3): M3 is similar to M2; the only difference from M2 is that the sine function hydrostatic restoring moment is considered when adjusting K PTO .M PTO for M3 is or "waves-by-waves", which means that the implementation time intervals of PTO damping tuning lie at the second or minute level.Furthermore, "look-up table" refers to the two key tables, i.e., the resonance position table and the optimal PTO damping table, derived from the MPP analysis for regular waves.Both of them are essential for guiding the MPPT in terms of regular waves.
When compared to other MPPT algorithms, MLTB MPPT has the following two features.
(1) MLTB MPPT has high reliability.It can effectively avoid search blindness, as the two look-up tables offer a meaningful reference for improving search accuracy.
(2) The principle of MLTB MPPT is simple.Complex wave prediction is unnecessary.Moreover, the two look-up tables can be easily obtained through a batch of regular wave simulations without regard given to irregular waves.

LTB Internal-Mass Position Adjustment
The internal mass position should be adjusted adaptively, along with the changeable real-time wave spectrum, to make sure that the NIPWEC can generate electricity at the resonant state most of the time.The core of the internal-mass position adjustment locates and keeps the natural period of the NIPWEC consistent with the characteristic period c T .Therefore, the adjusting process can be conducted as follows:  In terms of LTB internal-mass position adjustment, the precise adjustment of X can be realized via the MPAM with a screw and a controllable motor, e.g., a servo or stepper motor.The controllable motor drives the screw to rotate and to move the internal mass to a reference position X * , which can be given by querying the resonance position table according to the estimated characteristic period of a real-time wave spectrum.The real-time wave spectrum can be obtained through the FFT for the in situ wave elevation signal (WES) within the past hundreds of seconds.Herein, WES can be dependably achieved via the in situ wave-measuring buoy.Additionally, WEC motion can be utilized to identify wave excitation force/moment [14,34,35].Thereby, WES could be theoretically estimated via the identified wave excitation force/moment and the corresponding transfer function.However, few papers focus on the relevant issue.
For LTB PTO damping tuning, the accurate tuning of R PTO can be achieved by applying vector control to the PMSG using a three-phase PWM rectifier.The vector controller adopts a 0-d-axis-current control mode.Structurally, it only contains a current loop with the reference d-axis current i * d = 0.Meanwhile, the reference q-axis current i * q is given by a function of R PTOmax , which can be retrieved from the optimal PTO damping table according to the period and height of the last single wave or the statistical period and height of a previous multiple-wave group.
Herein, "multi-timescale" is reflected in two aspects.First, the required historical WES for the LTB internal-mass position adjustment and the LTB PTO damping tuning is different in terms of the signal length.The former needs a several-hundred-second WES to accurately calculate the characteristic period of a real-time wave spectrum, whereas the latter only needs the WES within the duration of a single wave or multiple waves.Second, there also exists a distinction between the implementation time intervals for the LTB internal-mass position adjustment and the LTB PTO damping tuning.Internal mass position can be adjusted at intervals of the kilo-second level since the obvious change of a real-time wave spectrum only occurs hourly [1,36].Meanwhile, PTO damping should be tuned wave-by-wave or "waves-by-waves", which means that the implementation time intervals of PTO damping tuning lie at the second or minute level.
Furthermore, "look-up table" refers to the two key tables, i.e., the resonance position table and the optimal PTO damping table, derived from the MPP analysis for regular waves.Both of them are essential for guiding the MPPT in terms of regular waves.
When compared to other MPPT algorithms, MLTB MPPT has the following two features.
(1) MLTB MPPT has high reliability.It can effectively avoid search blindness, as the two look-up tables offer a meaningful reference for improving search accuracy.
(2) The principle of MLTB MPPT is simple.Complex wave prediction is unnecessary.Moreover, the two look-up tables can be easily obtained through a batch of regular wave simulations without regard given to irregular waves.

LTB Internal-Mass Position Adjustment
The internal mass position should be adjusted adaptively, along with the changeable real-time wave spectrum, to make sure that the NIPWEC can generate electricity at the resonant state most of the time.The core of the internal-mass position adjustment locates and keeps the natural period of the NIPWEC consistent with the characteristic period T c .Therefore, the adjusting process can be conducted as follows: S1. Preset adjusting cycle, ∆T m .S2. Set sequence number, n t = 0. S3.Obtain the WES and save it.S4.Determine whether running time ≥ n t ∆T m .If the inequality is satisfied, go to S5.Otherwise, keep running S3.
S5. Implement FFT and filtering on the historical WES within the duration of dT m after the instant of n t ∆T m to obtain the real-time wave spectrum.
S6. Calculate the T c of the real-time wave spectrum.S7.Search the resonance position table and find out the resonance position X max corresponding to T c .If T c cannot be found in the resonance position table, an interpolation should be employed.
S8. Set X * = X max and input X * into the MPAM.S9.Move the internal mass to the given position X * via an MPAM and keep the position.S10.Determine whether the program needs to be stopped.If yes, terminate the program.Otherwise, n t = n t + 1 and return to S3.

LTB PTO Damping Tuning
If we assume that each single wave in an irregular wave sequence can be seen as a regular wave within the duration of a wave period, the corresponding R PTOmax given by an optimal PTO damping table can be applied wave-by-wave, i.e., every other wave period, to approach the real-time "optimal amplitude" condition.Although a single wave is not completely sinusoidal, the aforementioned assumption is essential and suitable for practical engineering.Furthermore, a 2D LTB MPPT algorithm with the same assumption has already been proposed and verified by Ref. [27].The only difference between the 2D LTB MPPT and the LTB PTO damping tuning here is that the control variable of the 2D LTB MPPT refers to the duty cycle of a buck-boost converter.
Similar to the 2D LTB MPPT, LTB PTO damping tuning can be implemented wave-bywave or waves-by-waves, i.e., at intervals of multiple wave periods.The first implementation strategy is named single-wave-period tuning (SWPT), whereas the second can be called multiple-wave-period tuning (MWPT).If frequent adjustment is not wanted, MWPT could be an alternative with comparable performance.Moreover, since there is not much difference in the periods/heights of the two adjacent single waves as well as the statistical periods/heights of the two adjacent multiple-wave groups, the corresponding R PTOmax of the last single wave or multiple-wave group is utilized for the present one here in order to avoid complex wave prediction.
The procedure of SWPT is described below.S1.Preset the single wave order number to n = −1.S7.Set i * q = f (R PTOmax ) and input it into the vector controller.Since the 0-d-axiscurrent control mode is adopted, f (R PTOmax ) can be expressed as where ω m is the PMSG angular velocity, n p is the number of pole pairs, φ f is the rotor flux, and k g is the acceleration ratio of a speed increase mechanism.S8.Tune R PTO to R PTOmax via the PMSG vector control and keep R PTO .S9. Determine whether the program needs to be stopped.If yes, terminate the program.Otherwise, return to S2.
In addition, the procedure of MWPT is as follows: S1.Preset the single wave order number to n = −1.
S2. Define the number of single waves k w (k w > 1).S3. Obtain the WES and save it.S4.Determine whether it reaches the n+1-th up-crossing zero point Z up,n+1 .If yes, n = n + 1 and go to S5.Otherwise, return to S3. S5.Determine whether n > 0. If yes, go to S6.Otherwise, return to S3. S6.Calculate the height H n and period T n of the n-th single wave.S7.Determine whether n is the multiple of k w .If yes, go to S8.Otherwise, return to S3. S8.Calculate the statistical height H k w and statistical period T k w of the multiple-wave group containing the last k w single waves.Herein, the mean wave height and period for the first one-third of the single waves of the multiple-wave group, which is rearranged in descending order of k w wave heights, is chosen as H k w and T k w .
S9. Search the optimal PTO damping table and find the R PTOmax corresponding to (H k w , T k w ).If (H k w , T k w ) cannot be exactly found in the optimal PTO damping table, an interpolation should be employed.S10.Set i * q = f (R PTOmax ) and input it into the vector controller.S11.Tune R PTO to R PTOmax via the PMSG vector control and keep R PTO .S12. Determine whether the program needs to be stopped.If yes, terminate the program.Otherwise, return to S3.

The Two Other PTO Damping Tuning Algorithms
The two other PTO damping tuning algorithms were also researched for a comparison with the LTB PTO damping tuning.The first one is a fixed-step P&O algorithm, and the second is a variable-step P&O algorithm, i.e., the admittance differentiation method (ADM).
The structure of a fixed-step P&O algorithm is shown in Figure 3. Herein, R PTO is perturbed with a fixed-step ∆R PTO in order to search the maximum absorbed power P PTO in real time.The update duration of R PTO is set as ∆T.This means R PTO is automatically updated at time intervals of ∆T.The perturbation direction of R PTO will be flipped once the P PTO value at present is lower than that at the previous moment.Similar algorithms were adopted in Refs.[19][20][21][22].

The Two Other PTO Damping Tuning Algorithms
The two other PTO damping tuning algorithms were also researched for a comparison with the LTB PTO damping tuning.The first one is a fixed-step P&O algorithm, and the second is a variable-step P&O algorithm, i.e., the admittance differentiation method (ADM).
The structure of a fixed-step P&O algorithm is shown in Figure 3. Herein, PTO R is perturbed with a fixed-step Δ PTO R in order to search the maximum absorbed power PTO P in real time.The update duration of PTO R is set as ΔT .This means PTO R is automatically updated at time intervals of ΔT .The perturbation direction of PTO R will be flipped once the PTO P value at present is lower than that at the previous moment.Similar algorithms were adopted in Refs.[19][20][21][22].ADM is based on the moving-least-square (MLS) power forecast and maximum capture width ratio (CWR) tracking.An analogous algorithm is described in Ref. [24].The structure of an ADM is displayed in Figure 4. Herein, PTO R is perturbed with a variable  ADM is based on the moving-least-square (MLS) power forecast and maximum capture width ratio (CWR) tracking.An analogous algorithm is described in Ref. [24].The structure of an ADM is displayed in Figure 4. Herein, R PTO is perturbed with a variable step size ∆R PTO_ADM .∆R PTO_ADM is given by where a means the step size coefficient, CWR the capture width ratio, R PTO_ADM,−1 or R PTO_ADM,−2 the PTO damping at the last moment or the moment before last, P PTO_ f orecast the absorbed power predicted by a MLS power forecast, and P wave as the input wave power P wave = J ir b.J ir is the wave power per unit width for irregular waves, and b is the NIPWEC dimension perpendicular to the wave propagation direction.
Energies 2023, 16, 6195 9 of 28 where a means the step size coefficient, CWR the capture width ratio, the PTO damping at the last moment or the moment before last,

Model: NIPWEC Dynamics
Based on Figure 5, the dynamic model of a NIPWEC can be expressed as the following Cummins' equation [37]. where

Model: NIPWEC Dynamics
Based on Figure 5, the dynamic model of a NIPWEC can be expressed as the following Cummins' equation [37].
where J P+M is the total moment of inertia, J 55∞ the additional mass at the wave frequency of +∞, θ is the swing angle, K(t) is the impulse response function of K(ω), F B is the buoyancy force, G P is the pendulum-hull gravity force, G M is the internal-mass gravity force, l B or l P is the distance between the rotation center R and the buoyancy center B or the pendulum-hull mass center P, M E is the wave excitation moment, M PTO is the PTO moment, M end is the endstop moment to prevent θ from exceeding its limits.

NIPWEC
We took the NIPWEC with a width of 5 m and a side area of 20 m 2 as a research case.The overall appearance and dimensions are illustrated in Figure 6.The detailed design parameters of this NIPWEC are listed in Table 1.In addition, the frequency-domain curves for the hydrodynamic parameters are shown in Figure 7. J P+M is the function of l M (t), and can be described by where J P is the pendulum-hull moment of inertia against R, J M0 is the internal mass moment of inertia against M, m M is the internal-mass mass.K(ω) can be represented as where R 55 (ω) and J 55 (ω) denote the radiation damping and additional mass against R, respectively.According to Ref. [38], M E can be written as where m is the number of wave frequency elements, ω j the j-th wave frequency element, f e5 and θ e5 are the amplitude and phase angle of the excitation-force coefficient, respectively, A is the wave amplitude, and rand() is a random number in (0,1).The expression of A has already been revealed in Ref. [27].This expression is based on the random interval division for wave frequency elements in order to avoid reproducing the WES or M E (t).
M PTO and M end can be given by M end = r end F end (11) where M g is the PMSG moment, and r end is the distance from R to the action point of the endstop force F end .Detailed formulas for M g and F end can be seen in Ref. [39].

Parameter Settings 4.1. NIPWEC
We took the NIPWEC with a width of 5 m and a side area of 20 m 2 as a research case.The overall appearance and dimensions are illustrated in Figure 6.The detailed design parameters of this NIPWEC are listed in Table 1.In addition, the frequency-domain curves for the hydrodynamic parameters are shown in Figure 7.

Irregular Wave Environments
Irregular wave simulations were implemented in MATLAB/Simulink.In order to study the applicability of MLTB MPPT for diverse irregular waves, the following three common wave spectrums are introduced to simulate irregular wave environments.
(1) JONSWAP spectrum [40] where s H and p T are the significant wave height and the peak period, respectively, γ (1 < γ < 10) means the peak elevation factor, σ signifies the peak shape parameter, and α is the energy scale parameter that ensures the validity of the following equation:

Irregular Wave Environments
Irregular wave simulations were implemented in MATLAB/Simulink.In order to study the applicability of MLTB MPPT for diverse irregular waves, the following three common wave spectrums are introduced to simulate irregular wave environments.
(1) JONSWAP spectrum [40] where H s and T p are the significant wave height and the peak period, respectively, γ (1 < γ < 10) means the peak elevation factor, σ signifies the peak shape parameter, and α is the energy scale parameter that ensures the validity of the following equation: wherein γ is positively correlated with the slenderness of a spectrum peak.When γ = 3.3, the wave spectrum can also be called a "standard JONSWAP spectrum".Moreover, σ is given segmentally in the frequency domain.The formula of σ is (2) Pierson-Moskowitz (P-M) spectrum [40]: where T z is the zero-crossing period.
(3) Ochi-Hubble spectrum [41]: where H s,i , T p,i , and λ i (i = 1 or 2) are the significant wave height, peak period, and shape parameter of the low/high frequency section, respectively.Herein, these six parameters can be separately given by a function of the significant wave height H s (H s = H 2 s,1 + H 2 s,2 ) for the Ochi-Hubble spectrum.
According to the above three wave spectra, eight sea states, i.e., SS1-SS8, can be defined in Table 2.Among them, SS1-SS5 are set as the standard JONSWAP spectra with different significant wave heights H s and peak periods T p in order to study the MPPT performance for the same wave-spectrum type but different H s or T p .Meanwhile, SS6-SS8 are defined as the JONSWAP spectrum, P-M spectrum, and Ochi-Hubble spectrum with the same H 2 s , respectively, for the sake of investigation into the MPPT performance of different wave-spectrum types.Moreover, the other investigated characteristic periods T c of each sea state are also shown in Table 2.The selection criteria of an investigated T c is based on whether it is common and easily available from the wave spectra of irregular wave environments, wherein the energy period T e , the mean period T m , and the zero-crossing period T z can all be conveniently obtained via spectral moments.The definition for a n-th spectral moment m n is Thereby T e , T m , and T z can be given by In addition, the wave spectrum density functions of the eight sea states are all displayed in Figure 8.

MPPT Implementation Methods
In order to examine the effectiveness of the three implementation methods, the NIP-WEC time-domain response in the regular wave with a wave period of 6 s and a wave height of 0.5 m was simulated.All the parameters are set for the optimal phase and amplitude.For M1,

MPPT 4.3.1. MPPT Implementation Methods
In order to examine the effectiveness of the three implementation methods, the NIP-WEC time-domain response in the regular wave with a wave period of 6 s and a wave height of 0.5 m was simulated.All the parameters are set for the optimal phase and amplitude.For M1, X = 1 m, R PTO =49,000 Nms.With regard to M2, K PTO = −1.055× 10 6 Nm, and R PTO = 79,000 Nms.In terms of M3, K PTO = −1.055× 10 6 Nm, and R PTO = 43,000 Nms.Moreover, the calculation formulas for the instantaneous input power P I N , absorbed power P PTO , PTO active power P R , and PTO reactive power According to M1, X max , R PTOmax , the maximum time-averaged absorbed power P PTOmax , and the maximum capture widths CW max corresponding to the diverse NIP-WEC MPPs were obtained by a series of time-domain simulations in the regular waves of different wave periods and wave heights.Herein, CW max is computed by where J r signifies the input wave power per unit width for a certain regular wave.

MLTB MPPT
The resonant states where the NIPWEC natural period is consistent with the four T c , i.e., T p , T e , T m , and T z , were comprehensively investigated for power assessments in terms of irregular waves.Herein, three prerequisites are set uniformly.First, the simulation duration of each case is set as 1000 s.Second, the optimal, fixed PTO damping R PTO_OFD against a resonant state is searched ergodically from 0 to 2 × 10 5 Nms (with an interval of 1 × 10 4 Nms) in order to find the maximum time-averaged absorbed power P PTO_OFD under the mode of fixed PTO damping.Third, for the sake of the sensitivity study, the resonant states of the NIPWEC natural period are equal to 0.95-1.05T c (with an interval of 0.025 T c ); these were researched.In addition, the sum δ var of the absolute values of the maximum positive/negative relative deviations was chosen as a unified sensitivity indicator for a certain variable var, e.g., P PTO_OFD or R PTO_OFD .A larger δ var means that var is more sensitive to a certain T c .
In order to study the time-averaged absorbed power and parameter sensitivity of different PTO damping tuning algorithms conveniently, the following three prerequisites are set to be consistent.First, the simulation duration of each case is set as 1000 s as well.Second, R PTO is restricted between 20,000 Nms and 200,000 Nms for whichever algorithm is adopted.Thirdly, all studied cases are based on the resonant state, where the natural period of a NIPWEC is identical to the energy period T e .
The selection criteria of an investigated key parameter for PTO damping tuning simulations are represented by whether the parameter is able to influence the MPPT tracking performance.
In terms of LTB PTO damping tuning, the investigated parameter is the number of single waves k w .A k w that is too large could weaken the real-time tracking performance and restrict the ability to obtain globally optimal solutions.However, if k w is too small, the negative power deviation for tuning processes could be accumulated quickly.This may also limit the capability of optimizing globally.Herein, the cases of k w = 1-10, i.e., the cases of the variable update duration mainly distributed within 4-90 s, are investigated for the comprehensive power assessments and sensitivity study.
The investigated parameters for a fixed-step P&O algorithm are the PTO-damping step size ∆R PTO and the update duration ∆T.A too large ∆R PTO and a too small ∆T will make the search process unstable.Meanwhile, a too small ∆R PTO and a too large ∆T will result in a reduction in the search speed.Therefore, the cases of ∆R PTO = 10,000 Nms and ∆T = 5-90 s (with an interval of 5 s), as well as the cases of ∆R PTO = 1000-20,000 Nms (with an interval of 1000 Nms) and ∆T = 10 s, were simulated in order.Moreover, the initial PTO damping was set to 20,000 Nms.
For an ADM, the step size coefficient a is chosen as an investigated parameter.Like the ∆R PTO for a fixed-step P&O algorithm, a too large a can make the search process unstable, whereas a too small a can result in a reduction in search speed.Therefore, the cases of a = 500-10,000 Nms 3 (with an interval of 500 Nms 3 ) were thoroughly researched.Moreover, the ADM sampling time was determined as 0.01 s, whereas its initial PTO damping was also set as 20,000 Nms.
In addition, the overall sensitivity value |Λ P | is proposed as a unified sensitivity indicator of the three PTO damping tuning algorithms.|Λ P | can be calculated as where m Λ is the total number of first-order difference quotients, P PTO is the time-averaged absorbed power against different PTO damping tuning algorithms, and p * the normalized value of a parameter.The normalization method of p * is given by A larger |Λ P | means that the relative fluctuation of the time-averaged absorbed power is more sensitive to a certain algorithm parameter.In general, M1 can achieve high active power without reactive power generation and the occurrence of an excessive PTO moment.Therefore, the subsequent MPPT research is based on M1.

MPP Analysis in Regular Waves
The resonance position max X , the optimal PTO damping (1) All three methods can achieve a basic in-phase situation in terms of the angular velocity .θ and the wave excitation moment M E , which results in the positive value of P I N most of the time (see Figure 9a-c).
(2) Both M1 and M3 can make the NIPWEC operate at the larger .θ, leading to the higher P R , when compared to M2.
(3) Both M2 and M3 generate enormous M PTO or P K , which substantially increases the difficulty of the PTO moment control.The M PTO amplitudes for M2 and M3 can reach approximately 10 times that of M E , while the P K amplitudes for them can attain more than five times that of P I N .However, M PTO for M1 has the same magnitude as M E .Moreover, the P K for M1 disappears all the time.
In general, M1 can achieve high active power without reactive power generation and the occurrence of an excessive PTO moment.Therefore, the subsequent MPPT research is based on M1.

MPP Analysis in Regular Waves
The resonance position X max , the optimal PTO damping R PTOmax , the maximum time-averaged absorbed power P PTOmax , and the maximum capture widths CW max against different MPPs are listed in Tables 3-6.It should be noted that the wave steepness of some regular waves, e.g., the regular wave of the wave period T = 3 s and the wave height H = 2.5 m, is greater than 0.1412.According to Schwartz's theory [42], the aforementioned regular waves do not exist since wave breaking happens.Therefore, the corresponding R PTOmax , P PTOmax , and CW max are all set as zero.Three key findings emerge.(1) X max is only related to T, but R PTOmax depends on both T and H.As T increases, R PTOmax shows a trend of first increasing and then decreasing.However, as H increases, R PTOmax only shows an increasing trend.
(2) P PTOmax is related to both T and H.Moreover, P PTOmax appears to have the same variation law as R PTOmax .
(3) CW max is also related to both T and H. Nevertheless, the variation law for CW max is different from that of R PTOmax or P PTOmax .As T increases, CW max increases first and then decreases.When H increases, CW max gradually decreases.
Moreover, both P PTOmax and CW max were compared with the corresponding solutions from the preliminary frequency-domain calculations in order to verify the correctness of the time-domain simulations (see Figure 10).The results indicate that P PTOmax or CW max is always less than or equal to the frequency-domain solution.When both T and H are small, the frequency-domain solution fits P PTOmax or CW max better than other cases.Moreover, the capture widths derived from the time/frequency-domain model are all below the theoretical upper bound.Hence, the time-domain simulation correctness is effectively verified.In order to further explain the variation law of MPPs, the regular waves of T = 6 s and H = 0.5-3 m, as well as the ones of T = 3-11 s and H = 1.5 m are taken as examples.The time-averaged absorbed power curves against different H or T are displayed in Figure 11.The results show that R PTOmax is usually larger than R 55 (ω r ).The difference between R PTOmax and R 55 (ω r ) becomes more obvious, along with the increase in H since R 55 (ω r ) is independent of H, whereas R PTOmax increases synchronously.In addition, the relationship between R PTOmax and P PTOmax presents an approximately quadratic correlation under the circumstance of the same T combined with a different H.Meanwhile, there exists a hysteresis phenomenon in the R PTOmax -P PTOmax relationship, under the circumstance of the same H combined with a different T.
the frequency-domain solution fits In summary, the MPPs in regular waves can be effectively gained via time-domain simulations.Moreover, the R PTOmax and R PTOmax -P PTOmax relationship is totally different from the frequency-domain solutions.

Power Assessments and Sensitivity Study for the LTB Internal-Mass Position Adjustment
The resonant states against different characteristic periods T c of the eight sea states were simulated.Moreover, the resonant states against 0.95-1.05T c were also considered for the sensitivity study.The statistical results of P PTO_OFD and R PTO_OFD are listed in Tables 7  and 8, respectively.The largest A P as well as the corresponding A R of each sea state is in bold and is underlined.Herein, A P signifies the arithmetic mean of P PTO_OFD , whereas A R means the arithmetic mean of R PTO_OFD .Moreover, Figure 12   (1) Although there is a certain possibility for all T c resonances to obtain the largest A P , it is better to keep resonant with T e because T e resonance can lead 50% of sea states to achieve the largest A P , as well as allowing other sea states to obtain a submaximal A P .
(2) When the largest A P is achieved, the unified sensitivity indicator δ P for P PTO_OFD is relatively low.This means that P PTO_OFD is not sensitive to the fine tuning of T c when P PTO_OFD appears around the largest A P .Moreover, T e resonance owns the minimum or subminimum δ P in six out of eight sea states.Hence, T e resonance is not sensitive to parameter changes in most sea states.
(3) Despite being in the same sea state, A R against different T c is quite different, and the change law of A R along with the variations in T c is elusive.Moreover, in contrast In summary, T e resonance is preferable in terms of power assessments when compared to other T c resonances.Moreover, the performance reliability of the T e resonance is high since its parameter sensitivity is low in most sea states.Hence, the subsequent simulations for three PTO damping tuning algorithms are all based on the T e resonance.

Power Assessments and Sensitivity Study for the LTB PTO Damping Tuning
Power assessments for the three PTO damping tuning algorithms, i.e., the fixed-step P&O algorithm and the ADM and LTB PTO damping tuning algorithms, are listed in Tables 9-11.The upper bound of the maximum or mean P T O P of three algorithms in each sea state is in bold and underlined.The results indicate the following three facts.
(  (1) In terms of the mean P PTO against the different parameters, LTB PTO damping tuning can obtain the maximum mean value in five out of eight sea states.Moreover, it can also achieve submaximum mean values in the other sea states.However, the fixed-step P&O algorithm can only get the maximum mean value in SS2 or SS3.Moreover, ADM has the best tracking performance for the SS8 of an Ochi-Hubble spectrum.
(2) For the maximum value of the time-averaged absorbed power against the different parameters, the obtaining probability of the three algorithms are not significantly different.
(3) With regard to |Λ P |, both the |Λ P | of the fixed-step P&O algorithm and the |Λ P | of the ADM are much larger than that of the LTB PTO damping tuning.Therefore, the former two algorithms are more sensitive to parameter changes than the LTB PTO damping tuning.Moreover, the fixed-step P&O algorithm is more sensitive to ∆T than to ∆R PTO .
Figure 13 visually displays the P PTO of the three algorithms in the form of a histogram, wherein the maximum P PTO and the minimum P PTO are represented by error bars.From Figure 13, it can be seen that the bar length, i.e., the mean P PTO , of the LTB PTO damping tuning is generally higher, and its error bar is the narrowest.Although the fixed-step P&O algorithm and ADM can obtain the maximum average absorbed power under certain sea states through parameter tuning, such as the fixed-step P&O algorithm in SS2, SS3, or SS7 (with the optimal {∆R PTO , ∆T} = {5000 Nms, 10 s}, {8000 Nms, 10 s} or {10,000 Nms, 80 s}, respectively), and the ADM in SS1 or SS8 (with the optimal a = 1500 Nms 3 or 4000 Nms 3 , respectively), the mean P PTO of them is generally low and owns a wider error bar.
under certain sea states through parameter tuning, such as the fixed-step P&O a in SS2, SS3, or SS7 (with the optimal { Δ P T O R , Δ T } = {5000 Nms, 10 s}, {8000 Nm {10,000 Nms, 80 s}, respectively), and the ADM in SS1 or SS8 (with the optimal Nms 3 or 4000 Nms 3 , respectively), the mean P T O P of them is generally low an wider error bar.Furthermore, we take SS2, SS5, SS6, and SS7 as examples.The variations in P PTO as a function of the algorithm parameters against the above sea states are displayed in Figure 14 for the sake of the visual indication of the specific impact of parameter changes on algorithm performance.It can be seen that P PTO fluctuates at a wide range with no obvious change pattern when the ∆T or ∆R PTO alters equidistantly.When a changes in an equal interval, the fluctuation range of the average absorption power of the ADM is smaller than that of the fixed-step P&O algorithm.However, there exist several sudden increases at certain values of a.When k w increases one by one, the variations in P PTO against the LTB PTO damping tuning is the smoothest without significant fluctuations.
In summary, the LTB PTO damping tuning is the least sensitive to parameter changes and has a high average absorbed power value.The overall reliability of this algorithm is high.Although the fixed-step P&O method and admittance differentiation method can search for larger power points against some specific parameters and sea states, the optimal parameters corresponding to the above power points have significant differences in different sea states, which makes it difficult to determine the universal optimal parameters.Moreover, the high sensitivity to parameter changes is another limitation.
equal interval, the fluctuation range of the average absorption power of the ADM is smaller than that of the fixed-step P&O algorithm.However, there exist several sudden increases at certain values of a .When In summary, the LTB PTO damping tuning is the least sensitive to parameter changes and has a high average absorbed power value.The overall reliability of this algorithm is high.Although the fixed-step P&O method and admittance differentiation method can search for larger power points against some specific parameters and sea states, the optimal parameters corresponding to the above power points have significant differences in different sea states, which makes it difficult to determine the universal optimal parameters.Moreover, the high sensitivity to parameter changes is another limitation.

Conclusions
This paper attempts to find a suitable MPPT technology for the optimized operation of a NIPWEC.First, a unique MPPT implementation method, i.e., adjusting the internalmass position to achieve "optimal phase" and tuning the linear PTO damping to achieve "optimal amplitude", was proposed.Then, MPPs in regular waves were discussed, and two lookup tables, i.e., a 1D resonance position table and a 2D optimal PTO damping table, were determined.Afterwards, an MLTB MPPT algorithm was put forward for the

Conclusions
This paper attempts to find a suitable MPPT technology for the optimized operation of a NIPWEC.First, a unique MPPT implementation method, i.e., adjusting the internalmass position to achieve "optimal phase" and tuning the linear PTO damping to achieve "optimal amplitude", was proposed.Then, MPPs in regular waves were discussed, and two lookup tables, i.e., a 1D resonance position table and a 2D optimal PTO damping table, were determined.Afterwards, an MLTB MPPT algorithm was put forward for the NIPWEC real-time operation in irregular waves, according to the aforementioned two lookup tables.MLTB MPPT consists of two core parts, i.e., the LTB internal-mass position adjustment and the LTB PTO damping tuning.Finally, the power assessments and sensitivity study for the LTB internal-mass position adjustment and the LTB PTO damping tuning were respectively conducted for eight simulated irregular-wave sea states of diverse wave spectra.The findings are as follows.
(1) It is preferable to choose M1, i.e., simultaneously adjusting the internal mass position and PTO damping, to realize MPPT.First, M1 has the same active power generation performance as the PTO impedance adjustment, with consideration given to the nonlinear effect of the hydrostatic restoring moment.Second, M1 does not generate excessive reactive power and PTO moment.
(2) For the LTB internal-mass position adjustment in irregular waves, energy period resonance possesses the highest possibility to obtain the maximum mean time-averaged absorbed power when compared with other characteristic period resonances.Moreover, the energy period resonance has high reliability since its parameter sensitivity is low in most sea states.
(3) When compared to the fixed-step P&O algorithm and a typical variable-step P&O algorithm, i.e., ADM, LTB PTO damping tuning is the most robust to parameter changes and has the highest mean time-averaged absorbed power in the majority of sea states.
In the future, the tracking performance of MLTB MPPT in changing sea states will be studied in-depth in order to further validate its applicability.
DC voltage and balances power generation and consumption via an energy storage unit and diverse electric loads.

Figure 1 .
Figure 1.Mechanical structure and electrical topology of the NIPWEC with a DC bus.

Figure 1 .
Figure 1.Mechanical structure and electrical topology of the NIPWEC with a DC bus.

Figure 2 .
Figure 2. Schematic diagram of the control structure for MLTB MPPT.

Figure 3 .
Figure 3. Block diagram of the fixed step P&O algorithm.

Figure 3 .
Figure 3. Block diagram of the fixed step P&O algorithm.

∞ 55 J
the additional mass at the wave frequency of ∞ + , θ is the swing angle, ( )Kt is the impulse response function of ω ( ) K , B F is the buoyancy force, P G is the pendulum-hull gravity force, M G is the internal-mass gravity force, B l or P l is the distance between the rotation center R and the buoyancy center B or the pendulum-hull mass center P, E M is the wave excitation moment, PTO M is the PTO moment, end M is the endstop moment to prevent θ from exceeding its limits.

Figure 5 .
Figure 5. Force analysis of the NIPWEC.Herein, F E , F R , G M , G P , and F B mean the wave excitation force, radiation force, internal-mass gravity force, pendulum-hull gravity force, and buoyancy force, respectively, M E , M R , M M , M P , and M B are the corresponding moments against R.

Figure 6 .
Figure 6.The appearance and dimensions of the researched NIPWEC.(a) A three-dimensional view.(b) The basic dimensions of the side profile.

Table 1 .
The design parameters of the researched NIPWEC.

Figure 6 .
Figure 6.The appearance and dimensions of the researched NIPWEC.(a) A three-dimensional view.(b) The basic dimensions of the side profile.

Figure 8 .
Figure 8. Wave spectrum density functions of eight sea states.

Figure 8 .
Figure 8. Wave spectrum density functions of eight sea states.

Figure 9 Figure 9 .
Figure 9 shows the NIPWEC time-domain response against three MPPT implementation methods.The findings are as follows.Energies 2023, 16, 6195 17 of 28

Figure 10 .Figure 10 .
Figure 10.Variations in the time-averaged absorbed power (a) and the capture width CW (b) as the function of T against different H . Note: solid lines: resonance and

Figure 10 .Figure 11 .
Figure 10.Variations in the time-averaged absorbed power (a) and the capture width CW (b the function of T against different H . Note: solid lines: resonance and visually displays the energy capture performance of different T c resonances in the form of a histogram.The maximum P PTO_OFD and the minimum P PTO_OFD are represented by error bars.The results demonstrate the following three things.

Energies 2023 ,
16, 6195 20 of 25 to P PTO_OFD , R PTO_OFD is still sensitive to the fine tuning of T c when P PTO_OFD appears around the largest A P .

A
against the largest P A of each sea state is in bold and underlined.

Figure 12 .
Figure 12.Histogram about the maximum time-averaged absorbed power _ P T O O F D P

)Figure 12 .
Figure 12.Histogram about the maximum time-averaged absorbed power P PTO_OFD at different resonant states against different characteristic periods, T c .Herein, the upper and lower bounds are displayed in the form of error bars.5.4.Power Assessments and Sensitivity Study for the LTB PTO Damping TuningPower assessments for the three PTO damping tuning algorithms, i.e., the fixed-step P&O algorithm and the ADM and LTB PTO damping tuning algorithms, are listed in Tables 9-11.The upper bound of the maximum or mean P PTO of three algorithms in each sea state is in bold and underlined.The results indicate the following three facts.

Figure 13 .Figure 13 .
Figure 13.Time-averaged absorbed power histogram against different PTO damping tu rithms.Herein, _ P T O O F D P

Figure 14 .
Figure 14.Variations in P T O P as a function of the parameters of different PTO damping tuning algorithms in four sea states.(a) Δ T of the fixed-step P&O algorithm.(b) Δ P T O R

Figure 14 .
Figure 14.Variations in P PTO as a function of the parameters of different PTO damping tuning algorithms in four sea states.(a) ∆T of the fixed-step P&O algorithm.(b) ∆R PTO of the fixed-step P&O algorithm.(c) a of the ADM.(d) k w of the LTB PTO damping tuning.

Figure 2. Schematic diagram of the control structure for MLTB MPPT.
S2. Obtain the WES and save it.S3.Determine whether it reaches the n+1-th up-crossing zero point Z up,n+1 .If yes, n = n + 1 and go to S4.Otherwise, return to S2. S4.Determine whether n > 0. If yes, go to S5.Otherwise, return to S2. S5.Calculate the height H n and period T n of the n-th single wave.S6.Search the optimal PTO damping table and find out the R PTOmax corresponding to (H n , T n ).If (H n , T n ) cannot be found exactly in the optimal PTO damping table, an interpolation should be employed. /

Table 1 .
The design parameters of the researched NIPWEC.

Table 2 .
Parameter settings for the simulated sea states.

Table 2 .
Parameter settings for the simulated sea states.

Table 4 .
Optimal PTO damping table.Herein background color is utilized to illustrate the value difference.The deeper red means the larger value, while the deeper blue refers to the smaller one.

Table 5 .
Maximum time-averaged absorbed power table.Herein background color is utilized to illustrate the value difference.The deeper red means the larger value, while the deeper blue refers to the smaller one.

Table 6 .
Maximum capture width table.Herein background color is utilized to illustrate the value difference.The deeper red means the larger value, while the deeper blue refers to the smaller one.

Table 7 .
Power assessments for resonant states against different characteristic periods T c .Herein, a sensitivity study is also conducted for 0.95-1.05T c .A P means the arithmetic mean of P PTO_OFD , and δ P is a unified sensitivity indicator for P PTO_OFD .The largest A P of each sea state is in bold and underlined.

Table 8 .
Analysis of the optimal, fixed PTO damping R PTO_OFD for the resonant states against different characteristic periods T c .Herein, a sensitivity study was also conducted for 0.95-1.05T c .A R means the arithmetic mean of R PTO_OFD , and δ R is a unified sensitivity indicator for R PTO_OFD .The A R against the largest A P of each sea state is in bold and underlined.

Table 9 .
Power assessments of the fixed-step P&O algorithm.Herein, a sensitivity study was also conducted for ∆R PTO = 10,000 Nms with ∆T = 5-90 s, as well as ∆T = 10 s, with ∆R PTO = 1000-20,000 Nms.The upper bound of the maximum or mean P PTO of three algorithms in each sea state is in bold and underlined.

Table 10 .
Power assessments of the ADM.Herein, a sensitivity study was also conducted for a = 500-10,000 Nms 3 .The upper bound of the maximum or mean P PTO of three algorithms in each sea state is in bold and underlined.

Table 11 .
Power assessments of the LTB PTO damping tuning.Herein, a sensitivity study was also conducted for k w = 1-10.The upper bound of the maximum or mean P PTO of three algorithms in each sea state is in bold and underlined.