Growth performance and physiological parameters of the second generation selected and control groups of Pinctada martensii

1. Introduction

Ocean wave energy refers to the harnessing of the herculean power of ocean waves. Ocean waves hold a gargantuan amount of untapped energy, some of which we can use to power at least a portion of the world’s everyday electricity. The ocean wave energy has the advantages of being highly predictable, renewable and eco-friendly. An engineering device utilized for exploiting the ocean wave power is called a wave energy converter (WEC). In order to successfully design a wave energy converter, accurately simulating the random ocean waves in the WEC dynamic analysis process is of uttermost importance. However, until present, the majority of the people in the worldwide wave energy research community have applied simple linear irregular waves in their WEC dynamic simulation processes (Manuel et al., 2018; Sirnivas et al., 2016; Tom et al., 2016, 2018, 2019). The linear irregular wave model has the disadvantages that it can only generate unrealistic waves with horizontal symmetries, i.e., the generated waves have statistically symmetric wave crests and troughs. This linear wave model is only suitable for approximately simulating random waves from a very mild sea state in a very deep sea. However, real world ocean waves will become statistically asymmetric (i.e., having sharper and higher crests but smoother and shallower troughs) in a harsh deep sea or at a shallow water coastal site.

Fernandes and Fonseca (2013) has pointed out that most of the proposed WECs will be installed and operated in shallow water coastal sites where the water depths are less than 90 m. Because the ocean waves in these shallow water sites will become statistically asymmetric, the linearly simulated statistically symmetric waves obviously should not be used as the inputs in the analysis and design of most of the proposed WECs.

In order to generate shallow water random waves with statistical asymmetries, Lindgren (2015) and Wang (2019) applied quasi-linear wave models for simulating the movements of individual water particles. However, the random wave simulations presented in Lindgren (2015) and Wang (2019) assumed that wave energy is traveling in only one direction (considered the same direction as the wind). That is to say, Lindgren (2015) and Wang (2019) had respectively performed their stochastic wave simulations based on a uni-directional wave spectrum. Similarly, when Wang (2018a, b) and Wang and Wang (2018) used nonlinear wave models for studying the power performances of wave energy converters, they also applied uni-directional wave spectra during their stochastic simulation of asymmetric waves.

In the real world, however, wind-generated ocean wave energy does not necessarily propagate in the same direction as the wind; instead, the ocean wave energy usually spreads over various directions. Therefore, for an accurate description of random seas, it is necessary to clarify the spreading status of energy. The wave spectrum representing energy in a specified direction is called the directional spectrum, denoted by $S\left( {\omega ,\theta } \right)$. Obviously, the energy of a sea state characterized by the wave spectrum $S\left( {\omega ,\theta } \right)$ continuously varies as the wave frequency $\omega $ and direction $\theta $ changes. However, to the best of this author’s knowledge, in the current literature, there exists no research work that has performed nonlinear wave simulations based on a multi-directional spectrum during the power performance analysis of wave energy converters.

Motivated by the aforementioned facts, in this paper the power performances of a WEC operating in a nonlinear random sea characterized by a multi-directional spectrum $S\left( {\omega ,\theta } \right)$ will be rigorously investigated. The absorbed power of the WEC Power-Take-Off system will be predicted by incorporating a second order random wave model into a nonlinear dynamic filter. This will be a new approach, and, as the second order random wave model can be utilized to accurately simulate the nonlinear waves in an irregular multi-directional sea, avoids the inaccuracies resulting from using a first order linear wave model in the simulation process. The predicted results in this paper will be systematically analyzed and compared, and the advantages of using this new approach will finally be substantiated.

The reminder of this paper will be organized as follows: In Section 2 the theories behind the second order random wave simulation based on a multi-directional spectrum will be elucidated, and the measured wave elevation data from a multi-directional coastal sea will be utilized to validate the accuracies of the second order nonlinear random wave simulation method. In Section 3 the theoretical background of the nonlinear dynamic filter of a wave energy converter will be provided. In Section 4 the calculation results of some specific calculation examples will be presented and discussed, with concluding remarks finally summarized in Section 5.

2. The theories and validation of a second order random wave simulation method

2.1 The theories of the second order random wave simulation method

The fluid region is described by using the 3D Cartesian coordinates (x, y, z), with x the longitudinal coordinate, y the transverse coordinate, and z the vertical coordinate (positive upwards). Time is denoted by t. The location of the free surface is at z = η(x, y, t) at a specific time of t.

For real fluid flows in the sea surface, it is reasonable to neglect the effects of viscosity. Meanwhile, sea water has a low compressibility, and therefore normally sea water flows are incompressible. Considering the aforementioned facts, it is reasonably accurate to assume that the fluids on the sea surface are ideal (i.e., incompressible and inviscid). Furthermore, rotation of a fluid particle can be caused only by a torque applied by shear forces on the sides of the particle. Since shear forces are absent in an ideal fluid, the flow of ideal fluids is essentially irrotational. Then, the velocity potential $\varPhi \left( {x,y,z,t} \right)$ exists. If the water depth d at the sea bottom is constant, then for constant water depth d, the velocity potential $\varPhi \left( {x,y,z,t} \right)$ and the free surface elevation η(x, y, t) can be determined by solving the following boundary value problem:

In this study the following expansion is utilized to solve the system (1)–(4):

where ε is a small parameter which is typically proportional to the wave steepness. For an irregular sea state characterized by a specific wave spectrum ${S_{\eta \eta }}\left( \omega \right)$ in which $\omega $ denotes the angular frequency, it can be shown that a first order linear solution of the system (1)–(4) can be expressed as follows:

where $N$ is usually chosen to be a sufficiently large positive integer, i stands for the imaginary unit, ${c_n}$ denotes the sinusoidal wave component amplitude, ${\varepsilon _n}$ is the phase angle uniformly distributed in the interval [0, 2${\text π} $], ${\omega _n}$ denotes the radian frequency (${\omega _n} = 2{\text π} n/T$, T is the time interval) and ${k_n}$ denotes the wave number. ${\omega _n}$ and ${k_n}$ are related through the following linear dispersion relation:

where $d$ and $g$ are respectively the water depth and the gravitational acceleration.

Equation (7) describes an ideal linear irregular sea model and this model has been derived under the assumption that the wave heights are small compared to the wave length. Then, the surface elevation resulting from irregular ocean waves can be approximated as a superposition of multiple harmonic waves with different amplitudes and phases. Furthermore, in this ideal linear irregular sea model the random phase angles, ${\varepsilon _n}$ are assumed to be uniformly distributed between 0 and 2π. However, for modelling shallow water nonlinear waves, the above linear irregular sea model should be corrected by including second order terms as follows:

The terms $P({\omega _n},{\omega _m})$, ${r_{mn}}$ and ${q_{mn}}$ in Eqs (9) and (10) are called quadratic transfer functions, and there expressions are given as follows:

where the Kroenecker delta (${\text{δ} _{ - n,m}}$=1, if n+m = 0, 0 otherwise) is introduced to avoid a singular $P({\omega _n},{\omega _m})$. The wave surface elevations for the second order nonlinear waves can finally be obtained by combining Eq. (7) and Eq. (10) as follows:

The proposed second order random wave simulation method starts with taking a multi-directional spectrum $S\left( {\omega ,\theta } \right)$ and integrating the energy over all directions to give the total energy at each frequency. The obtained equivalent frequency spectrum $S\left( \omega \right)$ is then utilized to generate a nonlinear wave elevation time series by applying Eqs (8) and (14) at a specific site of the sea. Finally it should be pointed out that the coefficients ${c_n}$(or ${c_m}$) in Eq. (14) are determined by using a given spectrum $S\left( \omega \right)$ as follows:

where $\Delta \omega = {\omega _C}/N$, and ${\omega _C}$ is the upper cut off frequency beyond which the wave spectrum $S\left( \omega \right)$ may be assumed to be 0 for either mathematical or physical reasons.

As we know, a wave spectrum is actually obtained from the observed time series of water elevation. The transformation from water elevation to wave spectrum is based on the assumption that irregular waves can be treated as a combination of linear waves with different amplitudes, different frequencies and different phases. It would be better to give an explanation for how to reasonably consider the nonlinearity of waves in the inverse transformation (i.e., the transformation from the wave spectrum to the nonlinear water elevation). In the following, we will raise a specific calculation example of this inverse transformation.

This study shows as an example to use Eq. (7) to simulate the linear Gaussian wave time histories of a sea state with a JONSWAP spectrum with a significant wave height ${H_S}$= 6 m, a spectral peak period ${T_p}$=8 s, a peakedness factor $\gamma$=1 and corresponding to a water depth of 10 m.

Figure 1 shows our simulation of the linear part wave elevation time series ${\eta ^{\left( 1 \right)}}(x,t)$ which contains 150 wave elevation points (with an equal time distant of 0.1 s between two successive points). A wave elevation point has two coordinates (the first coordinate is time, and the second coordinate is the value of the wave elevation above the mean water level). The computation is based on Eq. (7) and the aforementioned JONSWAP spectrum.

Figure  1.  The simulation of the linear part $ {\eta ^{\left( 1 \right)}}(x,t) $ which contains 150 wave elevation points.

DownLoad: Full-Size Img PowerPoint

Figure 2 shows our simulation of the second order correction part wave elevation time series ${\eta ^{\left( 2 \right)}}(x,t)$ which contains 150 wave elevation points (with an equal time distant of 0.1 s between two successive points). The computation is based on Eq. (10) and the aforementioned JONSWAP spectrum. This is precisely an example for how to reasonably consider the nonlinear correction part of waves in the inverse transformation (i.e., the transformation from the wave spectrum to the nonlinear correction part water elevation).

Figure  2.  The simulation of the second order correction part ${\eta ^{\left( 2 \right)}}(x,t)$ which contains 150 wave elevation points.

DownLoad: Full-Size Img PowerPoint

In Fig. 3 the blue curve shows our simulation of the entire $\eta (x,t)$ (the linear part ${\eta ^{\left( 1 \right)}}(x,t)$ plus the second order correction part ${\eta ^{\left( 2 \right)}}(x,t)$) wave elevation time series which contains 150 wave elevation points (with an equal time distant of 0.1 s between two successive points). The computation is based on Eq. (14) and the aforementioned JONSWAP spectrum. In Fig. 3, the red curve shows our simulation of the linear part wave elevation time series ${\eta ^{\left( 1 \right)}}(x,t)$ which contains 150 wave elevation points based on the aforementioned JONSWAP spectrum. In Fig. 3, the green curve shows our simulation of the second order correction part wave elevation time series ${\eta ^{\left( 2 \right)}}(x,t)$ which contains 150 wave elevation points based on the aforementioned JONSWAP spectrum. We can see that the wave crests of the entire nonlinear waves $\eta (x,t)$ have become steeper and higher. This is a precise example for how to reasonably consider the nonlinearity of the entire waves in the inverse transformation (i.e., the transformation from the wave spectrum to the entire nonlinear waves).

Figure  3.  The simulations of the blue curve entire nonlinear waves $\eta (x,t)$, the red curve linear waves ${\eta ^{\left( 1 \right)}}(x,t)$ and the green curve second order correction waves ${\eta ^{\left( 2 \right)}}(x,t)$, respectively, of 150 points of wave elevation time series.

DownLoad: Full-Size Img PowerPoint

2.2 The validation of the proposed second order random wave simulation method

In this sub-section, the accuracy of the proposed second order random wave simulation method will be validated by a calculation example. Specifically, the proposed second order random wave simulation method will be applied for calculating the wave crest amplitude exceedance probabilities of a sea state with a multi-directional wave spectrum based on the measured surface elevation data at the coast of Yura. The measured surface elevation data at the coast of Yura were obtained at a location 3 km off Yura fishing harbor facing the Sea of Japan. The observations were carried out during the period from 11:10 am to 2:08 pm on November 24, 1987 by the Ship Research Institute, Ministry of Transport of Japan. Temporal sea surface elevations were measured with ultrasonic-type wave gages installed at three points in 42 m water depth. The sampling time interval during the measurement was 1 s. Based on these measured Yura coast surface elevation data, a multi-directional wave spectrum was estimated by using the Maximum Likelihood Method and is shown in Fig. 4. Fig. 5 shows a part of the measured wave elevation time series by the mid-point wave gage at this site.

Figure  4.  A multi-directional wave spectrum from the measured data at the coast of Yura.

DownLoad: Full-Size Img PowerPoint

Figure  5.  A part of the measured wave elevation time series at the coast of Yura.

DownLoad: Full-Size Img PowerPoint

In the field of ocean engineering, based on a specific wave spectrum, an empirical (or theoretical) model or numerical simulation methods can be used to calculate the exceedance probabilities of the wave crest amplitudes (which is directly related to the occurrence probability of the extreme waves). In the following, taking the multi-directional wave spectrum in Fig. 4 as a calculation example, this study will test the proposed second order random wave simulation method to calculate the exceedance probabilities of the wave crest amplitudes, and compare its accuracy with those of a theoretical model or a linear simulation method.

Figure 6 shows the calculated wave crest amplitudes exceedance probabilities based on the multi-directional spectrum shown in Fig. 4 and the measured Yura coast wave data. The continuous green curve in Fig. 6 represents the calculation results of the wave crest amplitudes exceedance probabilities directly obtained from the measured Yura coast wave data that contains 10700 wave elevation points. This continuous green curves based on the measured Yura coast wave data are used as the benchmark against the accuracy of the results from the various numerical simulation methods and from an existing theoretical wave crest amplitudes model which is checked.

Figure  6.  Comparison between the wave crest amplitudes exceedance probabilities from the linear simulation, from the Rayleigh distribution, from the nonlinear simulation, and the wave crest amplitudes exceedance probabilities results from the measured Yura coast wave data.

DownLoad: Full-Size Img PowerPoint

The continuous red curve in Fig. 6 represents the results of the wave crest amplitudes exceedance probabilities obtained from using the theoretical Rayleigh distribution model expressed as follows:

where h represents the wave crest amplitude, ${H_{\rm{{S}}}}$ represents the significant wave height, A_c represents crest amplitude. The wave crest amplitudes exceedance probabilities obtained from using the theoretical Rayleigh distribution model deviate a lot from the corresponding benchmark results obtained directly from the measured Yura coast wave data. This is not a surprise because the theoretical Rayleigh distribution is an ideal linear model and the measured Yura coast wave elevation series are obviously nonlinear. The continuous black curve in Fig. 6 represents the results of the wave crest amplitude exceedance probabilities obtained from using the linear simulation method based on the multi-directional spectrum shown in Fig. 4. The linear simulation process started with taking the multi-directional spectrum shown in Fig. 4 and integrating the energy over all directions to give the total energy at each frequency. The obtained equivalent frequency spectrum was then used to generate a wave elevation time series of 800000 points by applying Eqs (7) and (8). Next, the wave crest amplitudes time series were extracted from these 800000 wave elevation points. Finally, the wave crests amplitudes exceedance probabilities were obtained by statistical and mathematical processing the extracted wave crest amplitudes time series. Obviously, the wave crest amplitudes exceedance probabilities obtained from using the linear simulation method also deviate substantially from the corresponding benchmark results obtained directly from the measured Yura coast wave data.

In order to more accurately calculate the wave crest amplitudes exceedance probabilities, the proposed second order random wave simulation method was tried to use. In Fig. 6, a pink * represents the result of the wave crest amplitudes exceedance probability obtained from using the proposed second order random wave simulation method based on the multi-directional spectrum shown in Fig. 4. The proposed second order random wave simulation process started with taking the multi-directional spectrum shown in Fig. 4 and integrating the energy over all directions to give the total energy at each frequency. The obtained equivalent frequency spectrum was then utilized to generate a nonlinear wave elevation time series of 800000 points by applying Eqs (8) and (14). Next, the wave crest amplitudes time series were extracted from these 800000 wave elevation points. Finally, the wave crests amplitudes exceedance probabilities were obtained by statistical and mathematical processing the extracted wave crest amplitudes time series. Note that in the figure legend the phrase “from nonlinear simulation” exactly means “from second order random wave simulation”. The wave crest amplitudes exceedance probabilities obtained from using the proposed second order random wave simulation method fit quite well with the corresponding benchmark results obtained directly from the measured Yura coast wave data. The accuracy of the proposed second order random wave simulation method is therefore convincingly validated.

3. The theoretical background of the nonlinear dynamic filter of a wave energy converter

In the practice ocean engineering, the power performances prediction of a WEC is usually carried out by inputting simulated ocean waves in a nonlinear dynamic filter and performing subsequent time domain simulations. The mathematical equations of the WEC nonlinear dynamic filter are written as follows:

where ${{{M}}_{{\rm{RB}}}}$ is the WEC rigid body inertia matrix and ${{x}}\left( t \right)$ is the WEC position. Meanwhile, ${{A}}\left( \infty \right)$ is the WEC infinite-frequency added mass matrix and ${{\bar {{K}}}}\left( t \right)$ is the kernel function. ${{{P}}_{{\rm{wave}}}}\left( t \right)$, ${{{P}}_{{\rm{ext}}}}\left( t \right)$, ${{{P}}_{{\rm{visc}}}}\left( t \right)$ and $ - {{{P}}_{{\rm{hs}}}}$ are the wave excitation, external, hydrodynamic viscous and hydrostatic loads (forces and moments) respectively.

In the field of ocean engineering, the WEC hydrodynamic memory effects are captured in Eq. (17) by the convolution integral term that is a function of ${\dot{ x}}\left( \tau \right)$ and the kernel functions ${{\bar {{K}}}}\left( t \right)$ which are related to the WEC radiation damping. Meanwhile in Eq. (17), the external load ${{{P}}_{{\rm{ext}}}}\left( t \right)$ represent the Power-Take-Off system forces and moments.

However, up until present, the wave excitation load ${{{P}}_{{\rm{wave}}}}\left( t \right)$ in the WEC dynamic filter are typically calculated by using linearly simulated irregular waves as inputs. As pointed out in Section 1, the linear irregular wave model has the disadvantages that it can only generate unrealistic waves with horizontal symmetries. Even though some researchers (Lindgren, 2015; Wang, 2018a, b, 2019; Wang and Wang, 2018) had tried to use quasi-linear wave models or nonlinear wave models for simulating irregular waves with statistical asymmetries, all these researchers unfortunately performed their stochastic wave simulations based on some unrealistic uni-directional wave spectra. Therefore, Section 1 illustrated the extension of these research work is needed very much in order to study the power performances of a wave energy converter in a multi-directional random sea. Motivated by these facts, in the present study, the wave excitation load ${{{P}}_{{\rm{wave}}}}\left( t \right)$ was calculated by inputting nonlinearly simulated irregular waves based on a multi-directional wave spectrum. In the next section, some calculation examples regarding the power performances of a point absorber wave energy converter operating in a multi-directional random sea will be presented.

4. Some specific calculation examples and calculation result discussion

4.1 The point absorber wave energy converter

Our specific calculation examples will be carried out regarding a heaving two-body point absorber wave energy converter. Figure 7 shows this specific wave energy converter modeled in WEC-Sim (http://wec-sim.github.io/WEC-Sim/), an open source computer software for simulating WEC performances. The dimensions of this specific wave energy converter are shown in Fig. 8. This two-body point absorber WEC consists of a float and a spar/plate. The float has a diameter of 20 m, a thickness of 5 m and a mass of 727.01 t. The spar has a height of 38 m and a diameter of 6 m. The spar/plate has a mass of 878.3 t.

Figure  7.  The WEC-Sim model of the chosen heaving two-body point absorber.

DownLoad: Full-Size Img PowerPoint

Figure  8.  The main dimensions of the chosen heaving two-body point absorber. SWL represents still water line.

DownLoad: Full-Size Img PowerPoint

It should be emphasized that this two-body point absorber WEC is installed in a coastal sea area with a water depth of only 49.5 m. Obviously, the random waves in this shallow water area will be nonlinear according to the nonlinear random wave theory in Ochi (1998).

4.2 Some specific calculation examples and calculation result discussion

In most wave energy exploitation engineering projects that met in the real world, the only information knew beforehand regarding a sea state usually will be a specific wave spectrum. In Section 2, the accuracy of the proposed second order random wave simulation method was already demonstrated in generating irregular waves based on a measured multi-directional wave spectrum. In the following calculation examples, the proposed second order random wave simulation method was utilized for generating nonlinear irregular waves for calculating the wave excitation load ${{{P}}_{{\rm{ext}}}}\left( t \right)$ in Eq. (17). The absorbed power of the aforementioned WEC Power-Take-Off system was predicted for some specific sea states characterized with multi-directional JONSWAP spectra. Figure 9 shows one of these multi-directional JONSWAP spectra with a significant wave height ${H_{\rm{{S}}}}$=1 m, a spectral peak period ${T_p}$=12 s, a peakedness factor $\gamma $=1, and a cosine squared spreading function with the spreading parameter equal to 15. The mathematical expression of this multi-directional JONSWAP spectrum is as follows:

Figure  9.  A multi-directional JONSWAP wave spectrum with ${H_{\rm{{S}}}}$=1 m, ${T_p}$=12 s, $\gamma $=1, and a cosine squared spreading function with the spreading parameter equal to 15.

DownLoad: Full-Size Img PowerPoint

where $S\left( \omega \right)$ is a frequency spectrum and $D\left( \theta \right)$ is a spreading function.

The mathematical expression for the frequency spectrum $S\left( \omega \right)$ is as follows (Wang, 2014):

where $\text{δ}= 0.07,\;{\rm{if}}\; \;\omega < 2{\pi} /{T_p};$ $\text{δ} = 0.09, \; {\rm{if}}\; \;\omega \geqslant 2{\pi} /{T_p}.$ $\omega $ is the wave angular frequency.

The mathematical expression for the spreading function $D\left( \theta \right)$ in Eq. (18) is as follows:

where the spreading parameter s=15.

Figure 10 is a corresponding polar plot of the multi-directional JONSWAP spectrum as shown in Fig. 9.

Figure  10.  A corresponding polar plot of the multi-directional JONSWAP spectrum as shown in Fig. 9.

DownLoad: Full-Size Img PowerPoint

Next, the calculation results of the absorbed power of the aforementioned wave energy converter were placed in an ideal linear sea versus in a multi-directional nonlinear random sea. Figure 11 shows the predicted WEC absorbed power time series under the sea state of linear irregular waves based on a multi-directional JONSWAP wave spectrum with ${H_{\rm{{S}}}}$=12 m, ${T_p}$=12 s, $\gamma $=1 and a cosine squared spreading function with the spreading parameter equal to 15. Figure 12 shows the predicted WEC absorbed power time series under the sea state of nonlinear irregular waves based on the same multi-directional JONSWAP wave spectrum with ${H_{\rm{{S}}}}$=12 m, ${T_p}$=12 s, $\gamma $=1 and a cosine squared spreading function with the spreading parameter equal to 15. The calculation results in Figs 11 and 12 were obtained by solving the WEC nonlinear dynamic filter Eq. (17) in WEC-Sim. However, WEC-Sim does not have built-in functions for generating nonlinear irregular waves for calculating the wave excitation load ${{{P}}_{{\rm{ext}}}}\left( t \right)$ in Eq. (17). Therefore, this study externally generated nonlinear irregular waves by using the proposed second order random wave simulation method which imported into WEC-Sim. The nonlinear irregular wave simulation was started by taking the multi-directional JONSWAP wave spectrum (with ${H_{\rm{{S}}}}$=12 m, ${T_p}$=12 s, $\gamma $=1) and integrating the energy over all directions to give the total energy at each frequency. The obtained equivalent frequency spectrum was then utilized to generate a nonlinear wave elevation time series of 1200 points by applying Eqs (8) and (14). The generated nonlinear irregular waves time series were saved as a.mat file and then imported into WEC-Sim for calculating the wave excitation load ${{{P}}_{{\rm{ext}}}}\left( t \right)$ in Eq. (17) and subsequently obtaining the WEC absorbed power time series in Fig. 12. For comparison purpose, in this study, externally generated linear irregular waves by using Eqs (7) and (8) were also imported into WEC-Sim. The linear irregular wave simulation was started by taking the multi-directional JONSWAP wave spectrum (with ${H_{\rm{S}}}$=12 m, ${T_p}$=12 s, $\gamma $=1) and integrating the energy over all directions to give the total energy at each frequency. The obtained equivalent frequency spectrum was then utilized to generate a linear wave elevation time series of 1200 points by applying Eqs (7) and (8). The generated linear irregular waves time series were saved as a .mat file and then imported into WEC-Sim for calculating the wave excitation load ${{{P}}_{{\rm{ext}}}}\left( t \right)$ in Eq. (17) and subsequently obtaining the WEC absorbed power time series in Fig. 11. Having obtained the prediction results as shown in Figs 11 and 12, the statistical characteristic values were subsequently calculated based on the time series shown in these two figures. Our calculation results are summarized in the last row in Table 1.

Figure  11.  WEC absorbed power time series under the sea state of linear waves with ${H_{\rm{{S}}}}$=12 m.

DownLoad: Full-Size Img PowerPoint

Figure  12.  WEC absorbed power time series under the sea state of nonlinear waves with ${H_{\rm{{S}}}}$=12 m.

DownLoad: Full-Size Img PowerPoint

Table  1.  Predicted 1 200 s WEC absorbed power values at the sea states of ${T_p}$=16 s, $\gamma $=1, s=15

Statistical measures significant wave height/m	Mean value under sea state of linear waves/W	Mean value under sea state of second order nonlinear waves/W	Standard deviation value under sea state of linear waves/W	Standard deviation value under sea state of second order nonlinear waves/W	Sum value under sea state of linear waves/W	Sum value under sea state of second order nonlinear waves/W
1	138437.04116	138511.45857	147967.47721	148037.48374	3.59798×108	3.59991×108
2	183201.45543	183729.92515	219208.08886	219533.04823	4.76141×108	4.77514×108
3	199926.33733	199983.50274	249739.80807	251762.62593	5.19609×108	5.19757×108
4	274600.57560	276111.71603	403700.40894	408900.02558	7.13687×108	7.17614×108
5	315077.75210	316133.23968	454553.94521	454063.30863	8.18887×108	8.21630×108
6	433808.09451	444136.33628	726265.21959	769908.28318	1.12747×109	1.15431×109
7	490131.33153	503258.63700	776567.66295	809812.45690	1.27385×109	1.30797×109
8	594029.94690	621211.75487	978890.28927	1.02946×106	1.54388×109	1.61453×109
9	796963.82665	874346.20999	1.34261×106	1.53251×106	2.07131×109	2.27243×109
10	767280.02013	825077.57589	1.17164×106	1.30227×106	1.99416×109	2.14438×109
11	1.1492×106	1.2581×106	2.21824×106	2.68345×106	2.98680×109	3.26980×109
12	1.0986×106	1.1918×106	2.53972×106	2.64326×106	2.85525×109	3.09748×109

 | Show Table

DownLoad: CSV

By carefully comparing and analyzing the calculation results, the mean value of the 1 200 s WEC absorbed power under the sea state of ideal linear irregular waves is smaller than the corresponding mean absorbed power value when inputting nonlinear irregular waves. In order to investigate the influences of choosing different significant wave height values on the power performances of wave energy converters, the aforementioned WEC absorbed power values were sequently calculated under 11 other sea states characterized with a multi-directional JONSWAP wave spectrum with the following parameters respectively: ${H_S}$=1 m, ${T_{\rm{p}}}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=2 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=3 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=4 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_S}$=5 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=6 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=7 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=8 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=9 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=10 m, ${T_p}$=12 s, $\gamma $=1, s=15; ${H_{\rm{S}}}$=11 m, ${T_p}$=12 s, $\gamma $=1, s=15. The calculation results were also summarized in Table 1. By carefully comparing and analyzing these calculation results, in all these 11 sea states, the mean value of the 1 200 s WEC absorbed power under the sea state of ideal linear irregular waves is always smaller than the corresponding mean absorbed power value when inputting nonlinear irregular waves.

In order to further investigate the influences of using different random seed numbers for generating irregular waves on the power performances of wave energy converters, the aforementioned WEC absorbed power values were sequently calculated when inputting 11 different irregular wave elevation time series simulated based on the same multi-directional JONSWAP wave spectrum (${H_{\rm{S}}}$=12 m, ${T_p}$=12 s, $\gamma $=1, s=15) using different random seed numbers. The calculation results were summarized in Table 2. By carefully comparing and analyzing these calculation results, in all these 12 sea states, the mean value of the 1 200 s WEC absorbed power under the sea state of ideal linear irregular waves is always smaller than the corresponding mean absorbed power value when inputting nonlinear irregular waves. In one specific case (No. 11) the predicted WEC absorbed power value with the nonlinear irregular waves as inputs is 15.8% larger than that predicted with the linear irregular waves as inputs. All the aforementioned calculation results highlight the vital importance of using the nonlinear irregular waves simulated based on a multi-directional wave spectrum when studying the power performances of wave energy converters.

Table  2.  Predicted 1 200 s WEC absorbed power values at the sea states of ${H_S}$=12 m, ${T_p}$=16 s, $\gamma $=1, s=15

Statistical measures significant wave height/m	Mean value under sea state of linear waves/W	Mean value under sea state of second order nonlinear waves/W	Standard deviation value under sea state of linear waves/W	Standard deviation value under sea state of second order nonlinear waves/W	Sum value under sea state of linear waves/W	Sum value under sea state of second order nonlinear waves/W
12	1.18485×106	1.38075×106	2.01713×106	2.76839×106	3.07941×109	3.58857×109
12	1.12270×106	1.22233×106	2.04328×106	2.41229×106	2.91791×109	3.17685×109
12	1.15923×106	1.27871×106	2.22910×106	2.63901×106	3.01285×109	3.32336×109
12	1.47223×106	1.58350×106	2.69001×106	3.28096×106	3.82632×109	4.11550×109
12	1.18207×106	1.31431×106	2.10716×106	2.70589×106	3.07219×109	3.41588×109
12	1.06410×106	1.12036×106	1.74052×106	1.80526×106	2.76560×109	2.91183×109
12	1.32898×106	1.44481×106	2.55403×106	2.81290×106	3.45403×109	3.75506×109
12	1.12296×106	1.31056×106	1.97489×106	2.46568×106	2.91857×109	3.40615×109
12	1.40660×106	1.57163×106	2.67732×106	3.01967×106	3.65575×109	4.08466×109
12	1.10365×106	1.18066×106	1.76090×106	2.01755×106	2.86838×109	3.06853×109
12	1.27946×106	1.48130×106	2.27087×106	2.72554×106	3.32531×109	3.84990×109
12	1.04726×106	1.13799×106	1.62224×106	1.92767×106	2.72184×109	2.95764×109

 | Show Table

DownLoad: CSV

5. Concluding remarks

In the present study the power performances of a two body point absorber WEC operating in a nonlinear multi-directional random sea have been rigorously investigated. The absorbed power of the WEC Power-Take-Off system has been predicted by incorporating a second order random wave model into a nonlinear dynamic filter. This is a new approach that is uniquely proposed to ocean wave energy research community. It has been demonstrated in this paper that the second order random wave model can be utilized to accurately simulate nonlinear irregular waves in a multi-directional sea. This will help us to avoid the inaccuracies resulting from using a first order linear wave model in the WEC simulation process. The predicted results in this paper have been systematically analyzed and compared, and the advantages of using our proposed new approach have been convincingly substantiated. The research findings in this paper highlight the vital importance of using the nonlinear irregular waves simulated based on a multi-directional wave spectrum when studying the power performances of wave energy converters. In the future if enough funding is secured, the proposed second order random wave simulation method will also be compared with the related WEC model physical test data in order to fully verify the accuracy of the model.