Optimization of Wave Energy Converter Arrays by an Improved Differential Evolution Algorithm

Abstract: Since different incident waves will cause the same array to perform differently with respect to the wave energy converter (WEC), the parameters of the incident wave, including the incident angle and the incident wave number, are taken into account for optimizing the wave energy converter array. Then, the differential evolution (DE) algorithm, which has the advantages of simple operation procedures and a strong global search ability, is used to optimize the wave energy converter array. However, the traditional differential evolution algorithm cannot satisfy the convergence precision and speed simultaneously. In order to make the optimization results more accurate, the concept of an adaptive mutation operator is presented to improve the performance of differential evolution algorithm. It gives the improved algorithm a faster convergence and a higher precision ability. The three-float, five-float, and eight-float arrays were optimized, respectively. It can be concluded that the larger the size of the array is, the greater the interaction between the floats is. Hence, a higher efficiency of wave energy extraction for wave energy converter arrays is achieved by the layout optimization of the array of wave energy converters. The results also show that the optimal layout of the array system is inhomogeneously distributed and that the improved DE algorithm used on array optimization is superior to the traditional DE algorithm.


Introduction
Consumption of fossil energy, such as oil, coal, and natural gas, has grown rapidly with the industrialization process.It is highly necessary to capture natural energy resources and convert them into electric energy or heat, going beyond non-renewable energy.Currently, the earth's surface area is 510 million km 2 , the ocean area of which is 361 million km 2 , accounting for 71% of the total area.It is estimated that the marine energy accounts for more than 70% of the world's total energy [1].China is a marine country whose sea area is vast and whose marine resources are abundant.Thus, the comprehensive utilization of marine energy has been proposed in China.
There is a large amount of renewable energy in the ocean, including wave energy, thermal energy, wind energy, salinity energy, tidal current energy, and tidal energy [2].Wave energy has gained considerable attention owing to its greater energy density compared to others.It is estimated that the wave energy reserve in the world is about 2.5 billion kW [3].As early as a century ago, people began to explore how wave energy can be used to generate electricity [4,5].There are various kinds of wave energy converters, such as point absorbers, attenuators, and terminators [6].
However, a single WEC has a low power efficiency, unstable power generation, and only small-scale generated electricity because it generally cannot absorb the ocean wave energy in different operator F used at an early stage prevents population diversity from local convergence and that a smaller F used later ensures accuracy with minimum error.The value of F is maintained in the best range for optimization results.The positions of each WEC in the array are determined by the improved DE algorithm, under more suitable incident wave parameters, so that greater wave energy is extracted from the ocean.The results show that the optimal layout is inhomogeneously distributed and that there are improved results after modification.

The Hydrodynamic Model of an Array System
In this paper, an improved oscillating buoy WEC system is adopted as shown in Figure 1.This system has a strong adaptability in nearshore and offshore regions.The float in the figure is connected with the reaction body through a pulley.The float can move up and down so that the wave energy can be converted to drive the permanent magnet synchronous generator (PMSG) continuously.The PMSG will then rotate in one direction under the automation control of the clutch and the bidirectional ratchet.The reaction body can maintain the cable tension and control the natural frequency of the WEC [27].
Energies 2018, 11, x FOR PEER REVIEW 3 of 18 smaller F used later ensures accuracy with minimum error.The value of F is maintained in the best range for optimization results.The positions of each WEC in the array are determined by the improved DE algorithm, under more suitable incident wave parameters, so that greater wave energy is extracted from the ocean.The results show that the optimal layout is inhomogeneously distributed and that there are improved results after modification.

The Hydrodynamic Model of an Array System
In this paper, an improved oscillating buoy WEC system is adopted as shown in Figure 1.This system has a strong adaptability in nearshore and offshore regions.The float in the figure is connected with the reaction body through a pulley.The float can move up and down so that the wave energy can be converted to drive the permanent magnet synchronous generator (PMSG) continuously.The PMSG will then rotate in one direction under the automation control of the clutch and the bidirectional ratchet.The reaction body can maintain the cable tension and control the natural frequency of the WEC [27].In this paper, it is supposed that the fluid is non-viscous and the fluid density ρ is invariable.The depth of the water is a finite value d, and the equation of wave free plane is z [28].The corresponding float model is cylindrical, and the state of it in seawater can be simplified as shown in Figure 2. In Figure 2, in the polar coordinate system, the submerging depth of the float is b, and the distance between the bottom of the float and the seafloor is h (h = d − b).The radius of the float is a.The size of the array is N (i.e., there are N floaters), and each float is numbered by j, j ∈ 1, 2, 3, …, N, In this paper, it is supposed that the fluid is non-viscous and the fluid density ρ is invariable.The depth of the water is a finite value d, and the equation of wave free plane is z [28].The corresponding float model is cylindrical, and the state of it in seawater can be simplified as shown in Figure 2.
Energies 2018, 11, x FOR PEER REVIEW 3 of 18 smaller F used later ensures accuracy with minimum error.The value of F is maintained in the best range for optimization results.The positions of each WEC in the array are determined by the improved DE algorithm, under more suitable incident wave parameters, so that greater wave energy is extracted from the ocean.The results show that the optimal layout is inhomogeneously distributed and that there are improved results after modification.

The Hydrodynamic Model of an Array System
In this paper, an improved oscillating buoy WEC system is adopted as shown in Figure 1.This system has a strong adaptability in nearshore and offshore regions.The float in the figure is connected with the reaction body through a pulley.The float can move up and down so that the wave energy can be converted to drive the permanent magnet synchronous generator (PMSG) continuously.The PMSG will then rotate in one direction under the automation control of the clutch and the bidirectional ratchet.The reaction body can maintain the cable tension and control the natural frequency of the WEC [27].In this paper, it is supposed that the fluid is non-viscous and the fluid density ρ is invariable.The depth of the water is a finite value d, and the equation of wave free plane is z [28].The corresponding float model is cylindrical, and the state of it in seawater can be simplified as shown in Figure 2. In Figure 2, in the polar coordinate system, the submerging depth of the float is b, and the distance between the bottom of the float and the seafloor is h (h = d − b).The radius of the float is a.The size of the array is N (i.e., there are N floaters), and each float is numbered by j, j ∈ 1, 2, 3, …, N, In Figure 2, in the polar coordinate system, the submerging depth of the float is b, and the distance between the bottom of the float and the seafloor is h (h = d − b).The radius of the float is a.The size of the array is N (i.e., there are N floaters), and each float is numbered by j, j ∈ 1, 2, 3, . . ., N, where the radius of float j is a j , and the quality is M j .When 0 r j a j , it represents the internal zone of float j.When r j a j , it represents the external zone of float j.The elastic effect of the cable can be equivalent to a spring with a coefficient of elasticity δ j and a damper with a damping coefficient of γ j [29].Further, it is assumed that the center coordinate O i of the float i is (0, 0), the center coordinate O j of the float j is (x j , y j ), the distance between O i and O j is L ij , the angle between the positive direction of the z-axis and the connecting line of float j and float i is θ ij , and β is the angle between the incident direction and the abscissa.The amplitude of the incident wave near field is H, the wave number is k 0 , the wave length is λ = 2π/k 0 , and the acceleration of gravity is g.The example of the two-float array in the Cartesian coordinate system is shown in Figure 3.
Energies 2018, 11, x FOR PEER REVIEW 4 of 18 where the radius of float j is aj, and the quality is Mj.When 0 ≪ rj ≪ aj, it represents the internal zone of float j.When rj ≫ aj, it represents the external zone of float j.The elastic effect of the cable can be equivalent to a spring with a coefficient of elasticity δj and a damper with a damping coefficient of γj [29].Further, it is assumed that the center coordinate Oi of the float i is (0, 0), the center coordinate Oj of the float j is (xj, yj), the distance between Oi and Oj is Lij, the angle between the positive direction of the z-axis and the connecting line of float j and float i is θij, and β is the angle between the incident direction and the abscissa.The amplitude of the incident wave near field is H, the wave number is k0, the wave length is λ = 2π/k0, and the acceleration of gravity is g.The example of the two-float array in the Cartesian coordinate system is shown in Figure 3.The analysis is under the theory of linearized potential flow, which has been validated in [30].The wave energy absorbed by the float can be calculated by the velocity potential φ, which can be decomposed into three components as [15,28,31] where φ I represents the velocity potential of the incident wave, φ S represents the scattered waves, φ R represents the radiated waves, φ D = φ I + φ S represents the velocity potential of the diffraction wave, and φ A represents the velocity potential of the incident wave near field.The φ A of the j-th float [32] is where the angular frequency ω satisfies that ω 2 = −km × g × tan(km × d), and km is the wave number of the evanescent wave.The superscript T in αj indicates the transpose, and the Ψj I and αj are represented as follows [14]: ) ) e , 0 ( ) The analysis is under the theory of linearized potential flow, which has been validated in [30].The wave energy absorbed by the float can be calculated by the velocity potential ϕ, which can be decomposed into three components as [15,28,31] where ϕ I represents the velocity potential of the incident wave, ϕ S represents the scattered waves, ϕ R represents the radiated waves, ϕ D = ϕ I + ϕ S represents the velocity potential of the diffraction wave, and ϕ A represents the velocity potential of the incident wave near field.The ϕ A of the j-th float [32] is where the angular frequency ω satisfies that , and k m is the wave number of the evanescent wave.The superscript T in α j indicates the transpose, and the Ψ j I and α j are represented as follows [14]: where J n is the n-order primal Bessel function, I n is the n-order modified primal Bessel function, and I j is the phase conversion factor of float j from the Cartesian coordinate system to the polar coordinate system.m (0 ≤ m ≤ +∞) is the discretization value of the z-axis in the Cartesian coordinate system.n (−∞ ≤ n ≤ +∞) is the discretization value of the θ axis in the polar coordinate system.
Energies 2018, 11, 3522 5 of 19 In the external zone of the j-th float, the scattered wave velocity potential ϕ S [31,33] is expressed as where A j T is a coefficient matrix of float j, and the matrix for coordinate transform T ji from Ψ j S to Ψ i I is represented as where H n is the n-order primal Hankel function, and K n is the n-order modified second Bessel function.
For the velocity potential of the incident wave near field, the value of the external zone of the float is different from that of the interior case.The internal scattering wave and the incident wave near field of the float are difficult to distinguish.Therefore, ϕ D is referred to as the diffraction wave field in the interior region of the float.In the interior region of the j-th float, ϕ D is expressed as where Ψ D j is calculated as I n (mπr j /h j ) where z = h j .Further, in the external region of the j-th float, the radiated wave field ϕ R can be represented as where R j (r j , z) is the radiation characteristic coefficient of the exterior region for float j, and Xj is the non-dimensional motion amplitude of float j, which is the movement amplitude of the float.
In the interior region of the j-th float, the radiated wave field ϕ R is represented as where R i (r i , z) is the radiation characteristic coefficient of the inner region for the j-th float.The The expression of the total velocity potential of incident wave is represented as A transformation matrix B j is adopted to show a relationship between the coefficient Ψ j I in Equation ( 11) and the coefficient Ψ j S in Equation ( 9) as Similarly, a transformation matrix B j could be introduced to represent the relationship between the coefficient Ψ j I in Equation ( 11) and the coefficient Ψ D j in Equation ( 7) within the float as Combining Equations ( 7), (10), and ( 13), the velocity potential inside the float j can be calculated, so it is easy to calculate the fluctuating force of the wave.
The Bernoulli equation is applied to integrate the bottom of the float S j .Then, in the positive direction of the z-axis, the hydrodynamic forces of float (F H j ) can be calculated as The elastic force F B j represents the buoyancy and gravity of a float in water, as follows where X j is the movement amplitude of float in the wave, which satisfies X j = H Xj .Furthermore, the float is pulled by the rope, and the elastic effect of the cable can be equivalent to a spring with an elastic coefficient of δ j and a damper with a damping factor of γ j .Therefore, the tension of a cable can be expressed as follows: where X j is the speed of the float in the wave, which satisfies X j = −iωH Xj .By combining Equations ( 15)-( 17), the following equation can be obtained: where X j is the acceleration of a float in the wave, which satisfies X j = −iω 2 H Xj , M j is the quality of the float, which satisfies M j = ρπa 2 j (d − h j ).A i and Xj are then given based on the above equations.Therefore, the wave energy extracted from each float can be calculated as [34] The energy gained by the WEC is the key to measuring the performance of a single WEC or a WEC array.However, whether the WEC array is better than a single WEC or what the criterion for measuring different kinds of WEC arrays should be are still open questions.Thus, an interaction coefficient q, in proportion to the total output power of the array system, is introduced to evaluate the quality of an array.The q is defined as follows: where P j is the energy gained by the float j in the array, P 0 is the energy gained by each isolated float.It has been found that the different incident waves will cause the same array to obtain different performances with respect to the wave energy converter.This is mainly due to the incident angle β and the incident wave number k 0 .Thus, the effect of these two parameters on the performance of the WEC will be studied in this paper.Its conclusion will be used to analyze whether the improved DE can be used to optimize the array layout.
In addition, in order to reflect the energy extraction of each float in an array, the interaction coefficient q j can also be used to reflect the operation condition of each float as

The Differential Evolution Algorithm
The DE algorithm, simple and powerful, can solve a problem by iterative optimization.It tries to improve the candidate solution and takes a given measure of quality into account, without making optimization problems differentiable.Therefore, the DE is widely used by people at home and abroad due to its accuracy, reliability, and robustness [35].Its efficiency and operation mainly depend on three procedures: mutation operation, crossover operation, and selection operation.The mutation scale factor F, population number Np, and crossover constant Cr are primary parameters that can control the calculation of the algorithm.Firstly, this algorithm generates randomly a population of Np candidate solutions, and each of them can be shown as x j i , where i ∈ 1, 2, 3, . . ., Np and j ∈ 1, 2, 3, . . ., D. Further, there are mutation operations and crossover operations to retrieve optimal values.In the end, a comparison is made between the previous generation and the individual with improved results [36,37].When the optimal solution is not obtained, the process is performed again.In this case, a maximum iteration number of Gm is defined.Note that each individual Np must carry on the mutation, crossover, and selection operations and does not directly replace the worst individual in the population.Therefore, the population size of each generation is constant, and each individual in the population will gradually approach the optimal value.For this reason, a threshold can be given.The optimal solution is not obtained until the fitness range of the population is over the threshold.The specific process of the differential evolution algorithm is presented as follows: (1) Initialization.Randomly generate the 0-th generation population X(0) = {X 1 (0), X 2 (0), . . ., X Np (0)}, where X i (0) = (x 1 i (0), x 2 i (0), . . ., x D i (0)).The initial population is chosen randomly under the given boundary constraints.It is generally assumed that all initialized populations satisfy the probability of uniform distribution.Set the bounds of the parameter variable as x j (L) < x j < x j (U) .Then where the value i is the integer between one and Np, j is the integer between one and D, and rand [0, 1] represents a series of random numbers between [0, 1].(2) Individual evaluation.Calculate every fitness value f (X i (G)) in the population.
The following mutation operations are performed for each X i (G) to generate a mutation vector (4) Crossover operation.Test vectors U i (G + 1) = (u 1 i (G + 1), u 2 i (G + 1), . . ., u D i (G + 1)) are obtained by the following crossover operation.
where randb(j) refers to the j-th estimation value of the stochastic number generator between [0, 1].rnbr(i), a stochastically selected sequence, is an integer between one and D. The function of rnbr(i) is to ensure that a value of V i (G + 1) can at least be obtained in U i (G + 1), like X i (G).( 5) Selection operation.Calculate the fitness value f (U i (G + 1)) of each test vector and compare them with f (X i (G)).Take the minimum, for example.There exists It is worth noting that each test vector competes only with the corresponding X i (G) rather than with each vector in the population.(6) Calculate the maximum and minimum for corresponding fitness value f (X(G + 1)) of the new population X(G + 1).Determine whether the difference between these two values is smaller than the threshold set in advance.If the calculation result of the difference is over the threshold and the number of iterations below its maximum (i.e., G < Gm), then repeat the above operation from Steps (2) to (6).

Improvement in the Differential Evolution Algorithm
A large amount of literature and experiments shows that, when scaling factor F ∈ [0, 2] is evaluated at [0.5, 1], the optimization results obtained will be superior [38].A larger F can produce a larger disturbance, which is conducive to maintaining the diversity of the population, but the efficiency of the search and the precision of the result are low.A smaller F with a strong local search capability is better able to obtain an optimal result.However, when F is small, it is easy for the population to lose its diversity quickly and precociously, which will lead to a local convergence.If a larger F can be used at an early stage to prevent the population diversity from local convergence, and a smaller F is used later to ensure the accuracy with minimum error, it can be predicted that the optimized results will be better than those where F is fixed.Therefore, the concept of an adaptive mutation operator is introduced as Equation (26) shows that the range of F is 2F 0 -F 0 .If F 0 = 0.5, the variation range is exactly 1-0.5, which not only satisfies the ideal range of F, but also can improve the result of Equation ( 23) as discussed above.
When crossover probability factor CR ∈ [0, 1] is smaller, a higher accuracy can be ensured.Theoretically, CR = 0.1 works best.However, because of its extremely slow convergence rate, it is not Energies 2018, 11, 3522 9 of 19 often used in practice.Increasing CR properly is beneficial to improve the convergence rate.It is often considered that CR = 0.9 is more effective in most cases [39].
Furthermore, the explanation of applying an improved DE algorithm on the optimization of an array is as follows: Firstly, the initialization is the same as the DE algorithm.Suppose the size of the array is N. Then there are N − 1 variables that need to be optimized.Each float's coordinate contains two components: the vertical and horizontal coordinates.Therefore, the dimension is D = 2 × (N − 1).Secondly, calculate the interaction coefficient of each float q x (x 2 , y 2 , . . ., x N , y N ), the optimum objective function of the array layout varying with the position of the float, at the initial position.Thirdly, set F 0 to 0.5.The optimal range of F is 1-0.5 in the mutation operation.Generate mutation vector V i (G + 1) = q best (G) + F × (q r1 (G) − q r2 (G)) in the form of DE/best/1/bin.Fourth, set CR to 0.9, and generate vectors U i (G + 1) by comparing the j-th estimation value with CR in the crossover operation.Fifth, calculate q u , the interaction coefficient of U i (G + 1), and compare it with q x .A better value is then selected.Finally, if G G m or the difference between the maximum and minimum of q u is larger than 0.001, end the procedure and determine the optimal value, or go on.

Parameters and Interaction Coefficient Analysis for WEC Arrays
It is well known that the size of the array and the distance between WECs will have a major impact on the output and efficiency of devices.However, the parameters of the incident wave, including the incident angle β and incident wave number k 0 , will cause the same array to perform differently.For an optimal analysis, the influence of β and k 0 on WEC output should be studied first.Setting incident wave parameters at optimal values will be beneficial for analysis.
The solution to the hydrodynamic model of the array system was calculated in MATLAB.All floats were calculated on a uniform scale to facilitate comparison with the isolated float.The values of the system parameters as shown in Table 1.Equation (20) shows that the interaction coefficient q is related to k 0 and β.In order to better optimize the array, the relationship among the interaction coefficient q, the incident wave number k 0 , and the incident angle β deserves to be investigated.

The Relationship between q and k 0
It was found that the effect of the wave number on q is related to the radius of the float.Thus, set the radius of the float to a constant and the wave number to a variable.The 2ak 0 is exactly the abscissa.Then, calculate the interaction coefficient q to obtain a graph when there are different sizes of float or different float spacings in MATLAB.Furthermore, set the wave number to a variable and the wave number to a constant to repeat the operation above.The optimal value of wave number can be chosen according to the calculation above.Figure 4 shows the relationship between q and k 0 .

The Relationship between q and β
Figure 5 shows the influences of the incident angle β on q with different amounts of floats and float spacing.As shown in Figure 5, the q values vary as the incident angle.Figure 5a-d, with the array size and float spacing set the same as in Figure 4, show that, if the array or the float spacing is different, a larger q can be obtained at β = 0. Therefore, β = 0 was selected in the following calculations as shown in Table 1.It can be seen that the value of q will change when the value of wave numbers is different.Optimizing the WEC array is the purpose here.However, it is impossible to optimize the array under each wave number.Figure 4a,b show q under the conditions of the same array size and two floats, but different float spacing, where the float spacing is equal to four times the radius in Figure 4a and eight times the radius in Figure 4b. Figure 4c,d show q under the conditions of the same float spacing and four times the radius, but different array sizes, with three floats in Figure 4c and four floats in Figure 4d.These four subfigures show that a larger q can be obtained in the vicinity of 2ak 0 = 0.8 with arrays of different sizes and float spacings.Therefore, the values of k 0 are fixed with the conditions of 2ak 0 = 0.8 in the next optimization, as shown in Table 1.

The Relationship between q and β
Figure 5 shows the influences of the incident angle β on q with different amounts of floats and float spacing.As shown in Figure 5, the q values vary as the incident angle.Figure 5a-d, with the array size and float spacing set the same as in Figure 4, show that, if the array or the float spacing is different, a larger q can be obtained at β = 0. Therefore, β = 0 was selected in the following calculations as shown in Table 1.
Figure 5 shows the influences of the incident angle β on q with different amounts of floats and float spacing.As shown in Figure 5, the q values vary as the incident angle.Figure 5a-d, with the array size and float spacing set the same as in Figure 4, show that, if the array or the float spacing is different, a larger q can be obtained at β = 0. Therefore, β = 0 was selected in the following calculations as shown in Table 1.Select 2ak 0 = 0.8 and β = 0, and take a two-float array as an to verify whether the k 0 and β discussed above are correct.As can be seen from Figures 4 and 5, when the float spacing is eight times the radius of the float (5 m), the q is larger than that when float spacing is four times the radius of the float.The two-float array is simple and accurate enough to show the difference.Thus, the two-float array was selected and the float spacing was set to 40 m. Figure 6 presents the corresponding graphics of various velocity potentials.
Figure 6a shows that the values and their distributions of the scattered wave velocity potential of Floats 1 and 2 are different.This indicates that these two floats will affect each other.Figure 6b shows that the velocity potentials of the radiation wave for Floats 1 and 2 are not the same.However, radiation waves are generated by the up and down movement of the float, which shows that the movements of each float are not the same with the interaction between floats.Therefore, it is necessary to calculate the movements of each float.At the same time, the differences between the floats in Figure 6a-c is further proof that the captured energy obtained in the array mode is different from that in isolated mode.This will be further verified with the following calculations.
radiation waves are generated by the up and down movement of the float, which shows that the movements of each float are not the same with the interaction between floats.Therefore, it is necessary to calculate the movements of each float.At the same time, the differences between the floats in Figure 6a-c is further proof that the captured energy obtained in the array mode is different from that in isolated mode.This will be further verified with the following calculations.

Simulation and Analysis of WEC Arrays
The simulation of WEC arrays is operated in MATLAB.Figure 7 shows a schematic of the array position.

Simulation and Analysis of WEC Arrays
The simulation of WEC arrays is operated in MATLAB.Figure 7 shows a schematic of the array position.As shown in Figure 7, the position of Float 1 is fixed at the position of (0, 0).The variable that needs to be optimized is the coordinates of the other float.Suppose the size of the array is N. Then there are N − 1 variables that need to be optimized, and the dimension is D = 2 × (N − 1).
The optimal results q depend on the required precision and the parameters of the improved DE algorithm, including the scaling factor, the crossover probability factor, and the size of the population.The required precision is set to 0.001, the range of the scaling factor is between 0.5 to 1, the crossover probability factor is set to 0.9, and the size of the population is set to 10.The given initial positions of the floats cannot influence the optimal layout and output of the WEC arrays.It will only affect the computation speed of the algorithm, whose difference is not so great.For example, if the initial position is closer to the optimal position, the computation speed of the algorithm will be faster and vice versa.The initial conditions are randomly set in the algorithm.Then, during the optimization of array, there will initially be a random population of float positions to determine a better position.It was found that the same optimal layout results can be achieved with different initial conditions.

Simulation of Three-Float, Five-Float, and Eight-Float Arrays
Three typical cases were investigated.Figure 8a shows the optimal result of a simulation of when there are three floats in the sea, and Figure 8b shows the individual fitness values of each generation, which reflects the maximum individual and minimum individual in each group.The blue line reflects the minimum individual in each generation, while the red line reflects the maximum.When the maximum and the minimum individual fitness value in the population is less than the set threshold, that is, when the red line and the blue line in the figure are close to coincidence, the result has converged to its optimal solution.As shown in Figure 7, the position of Float 1 is fixed at the position of (0, 0).The variable that needs to be optimized is the coordinates of the other float.Suppose the size of the array is N. Then there are N − 1 variables that need to be optimized, and the dimension is D = 2 × (N − 1).
The optimal results q depend on the required precision and the parameters of the improved DE algorithm, including the scaling factor, the crossover probability factor, and the size of the population.The required precision is set to 0.001, the range of the scaling factor is between 0.5 to 1, the crossover probability factor is set to 0.9, and the size of the population is set to 10.The given initial positions of the floats cannot influence the optimal layout and output of the WEC arrays.It will only affect the computation speed of the algorithm, whose difference is not so great.For example, if the initial position is closer to the optimal position, the computation speed of the algorithm will be faster and vice versa.The initial conditions are randomly set in the algorithm.Then, during the optimization of array, there will initially be a random population of float positions to determine a better position.It was found that the same optimal layout results can be achieved with different initial conditions.6.1.Simulation of Three-Float, Five-Float, and Eight-Float Arrays Three typical cases were investigated.Figure 8a shows the optimal result of a simulation of when there are three floats in the sea, and Figure 8b  As shown in Figure 7, the position of Float 1 is fixed at the position of (0, 0).The variable that needs to be optimized is the coordinates of the other float.Suppose the size of the array is N. Then there are N − 1 variables that need to be optimized, and the dimension is D = 2 × (N − 1).
The optimal results q depend on the required precision and the parameters of the improved DE algorithm, including the scaling factor, the crossover probability factor, and the size of the population.The required precision is set to 0.001, the range of the scaling factor is between 0.5 to 1, the crossover probability factor is set to 0.9, and the size of the population is set to 10.The given initial positions of the floats cannot influence the optimal layout and output of the WEC arrays.It will only affect the computation speed of the algorithm, whose difference is not so great.For example, if the initial position is closer to the optimal position, the computation speed of the algorithm will be faster and vice versa.The initial conditions are randomly set in the algorithm.Then, during the optimization of array, there will initially be a random population of float positions to determine a better position.It was found that the same optimal layout results can be achieved with different initial conditions.

Simulation of Three-Float, Five-Float, and Eight-Float Arrays
Three typical cases were investigated.Figure 8a shows the optimal result of a simulation of when there are three floats in the sea, and Figure 8b shows the individual fitness values of each generation, which reflects the maximum individual and minimum individual in each group.The blue line reflects the minimum individual in each generation, while the red line reflects the maximum.When the maximum and the minimum individual fitness value in the population is less than the set threshold, that is, when the red line and the blue line in the figure are close to coincidence, the result has converged to its optimal solution.Table 2 shows the energy capture coefficients of the WEC array.Similarly, a five-float array was optimized.Figure 9a shows the optimal result of the simulation diagram when there are five floats in the sea. Figure 9b shows the individual fitness values of each generation.Table 3 shows the corresponding energy capture coefficients of the WEC array.Similarly, a five-float array was optimized.Figure 9a shows the optimal result of the simulation diagram when there are five floats in the sea. Figure 9b shows the individual fitness values of each generation.Table 3 shows the corresponding energy capture coefficients of the WEC array.An eight-float array was then optimized.Figure 10 shows the optimal results for the eight-float array.An eight-float array was then optimized.Figure 10 shows the optimal results for the eight-float array.
Energies 2018, 11, x FOR PEER REVIEW 13 of 18 Table 2 shows the energy capture coefficients of the WEC array.Similarly, a five-float array was optimized.Figure 9a shows the optimal result of the simulation diagram when there are five floats in the sea. Figure 9b shows the individual fitness values of each generation.Table 3 shows the corresponding energy capture coefficients of the WEC array.An eight-float array was then optimized.Figure 10 shows the optimal results for the eight-float array.The optimal position result of each float when there are eight floats in the sea is shown in Figure 10a.The individual fitness values of each generation are shown in Figure 10b.Table 4 shows the corresponding energy capture coefficients of the WEC array.From the optimized WEC arrays with N = 3, N = 5, and N = 8, it can be seen that the larger the array is, the greater the interaction among the floats is.Hence, higher efficiency of the extracted wave energy is achieved.This means that the array composed of multiple floats can extract more wave energy compared with that of the single float operation mode.

Comparison with the Traditional DE Algorithm
In order to determine whether the improved DE algorithm is superior to the traditional DE algorithm, a comparison of results for these two algorithms was analyzed.Take a three-float array for example due to its simplicity and representativeness.The optimal position result of each float when there are three floats in the sea is shown in Figure 11a, and the individual fitness values of each generation is shown in Figure 11b.Table 5 shows the corresponding energy capture coefficients of the WEC array.
Energies 2018, 11, x FOR PEER REVIEW 14 of 18 The optimal position result of each float when there are eight floats in the sea is shown in Figure 10a.The individual fitness values of each generation are shown in Figure 10b.Table 4 shows the corresponding energy capture coefficients of the WEC array.From the optimized WEC arrays with N = 3, N = 5, and N = 8, it can be seen that the larger the array is, the greater the interaction among the floats is.Hence, higher efficiency of the extracted wave energy is achieved.This means that the array composed of multiple floats can extract more wave energy compared with that of the single float operation mode.

Comparison with the Traditional DE Algorithm
In order to determine whether the improved DE algorithm is superior to the traditional DE algorithm, a comparison of results for these two algorithms was analyzed.Take a three-float array for example due to its simplicity and representativeness.The optimal position result of each float when there are three floats in the sea is shown in Figure 11a, and the individual fitness values of each generation is shown in Figure 11b.Table 5 shows the corresponding energy capture coefficients of the WEC array.Based on the comparison between the improved DE algorithm and traditional DE algorithm, the total energy of the three-float array absorbed from the wave under the traditional DE algorithm is lower than that of the improved DE algorithm.The convergence rate is slower than that of the improved DE algorithm.Furthermore, the program running time of five floats under the improved DE algorithm is about 10 min and about 12 min under the traditional DE algorithm.The program running times of the three floats, under these two algorithms, are both no more than 2 min.If the array is larger, the difference will be more obvious.That is, the improved DE algorithm presented in this paper is effective for the layout optimization of an WEC array.The introduction of an adaptive mutation operator can obtain a better optimal result compared to the traditional DE algorithm.

Comparison with the Homogeneous Distribution of the Three-Float Array
It was found that the optimum of the three cases discussed above is not in a homogeneous distribution.Thus, it deserves further analysis compared with the array under the homogeneous distribution.Take a simple and representative three-float array for example.Three well-distributed arrays are analyzed as follows, whose distance and distribution are slightly similar to the layout results under the algorithm.Figure 12 shows the layout and velocity potential of a homogeneous distribution for a three-float array.Table 6 shows the energy capture coefficient of the WEC array (N = 3) under these two typical homogeneous distributions.Based on the comparison between the improved DE algorithm and traditional DE algorithm, the total energy of the three-float array absorbed from the wave under the traditional DE algorithm is lower than that of the improved DE algorithm.The convergence rate is slower than that of the improved DE algorithm.Furthermore, the program running time of five floats under the improved DE algorithm is about 10 min and about 12 min under the traditional DE algorithm.The program running times of the three floats, under these two algorithms, are both no more than 2 min.If the array is larger, the difference will be more obvious.That is, the improved DE algorithm presented in this paper is effective for the layout optimization of an WEC array.The introduction of an adaptive mutation operator can obtain a better optimal result compared to the traditional DE algorithm.

Comparison with the Homogeneous Distribution of the Three-Float Array
It was found that the optimum of the three cases discussed above is not in a homogeneous distribution.Thus, it deserves further analysis compared with the array under the homogeneous distribution.Take a simple and representative three-float array for example.Three well-distributed arrays are analyzed as follows, whose distance and distribution are slightly similar to the layout results under the algorithm.Figure 12 shows the layout and velocity potential of a homogeneous distribution for a three-float array.Table 6 shows the energy capture coefficient of the WEC array (N = 3) under these two typical homogeneous distributions.It can be seen that the homogeneous distribution cannot improve output compared with the layout calculated by the improved DE algorithm.In Figure 12, the distances between floats of the equicrural triangle layout and the linear layout are similar to the obtained optimal position calculated by the algorithm.However, their interaction coefficients q are far lower than the layout calculated by the algorithm.

Figure 2 .
Figure 2. Simplified diagram of the float in wave.

Figure 2 .
Figure 2. Simplified diagram of the float in wave.

Figure 2 .
Figure 2. Simplified diagram of the float in wave.

Table 1 .
Parameters of the float and the WEC array.

Figure 7 .
Figure 7.A schematic of the array position.

Figure 8 .
Figure 8. Optimal results for the three-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 7 .
Figure 7.A schematic of the array position.

Figure 7 .
Figure 7.A schematic of the array position.

Figure 8 .
Figure 8. Optimal results for the three-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 8 .
Figure 8. Optimal results for the three-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 9 .
Figure 9. Optimal results for the five-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 10 .
Figure 10.Optimal results for the eight-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 9 .
Figure 9. Optimal results for the five-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 9 .
Figure 9. Optimal results for the five-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 10 .
Figure 10.Optimal results for the eight-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 10 .
Figure 10.Optimal results for the eight-float array: (a) the optimal position; (b) individual fitness values for each generation.

Figure 11 .
Figure 11.Optimal results for the three-float array under the traditional DE algorithm: (a) the optimal position; (b) individual fitness values for each generation.

Figure 11 .
Figure 11.Optimal results for the three-float array under the traditional DE algorithm: (a) the optimal position; (b) individual fitness values for each generation.

Figure 12 .
Figure 12.Layout and velocity potential of homogeneous distribution for the three-float array: (a) equicrural triangle; (b) straight line.

Figure 12 .
Figure 12.Layout and velocity potential of homogeneous distribution for the three-float array: (a) equicrural triangle; (b) straight line.

Table 2
shows the energy capture coefficients of the WEC array.

Table 4 .
The situation of each float when N = 8.

Table 4 .
The situation of each float when N = 8.

Table 5 .
Energy capture coefficient of the WEC array under the traditional DE algorithm (N = 3).

Table 5 .
Energy capture coefficient of the WEC array under the traditional DE algorithm (N = 3).

Table 6 .
Energy capture coefficient of the WEC array (N = 3) under different homogeneous distributions.

Interaction Coefficient qj (pu) Interaction Coefficient q (pu)
It can be seen that the homogeneous distribution cannot improve output compared with the