Hydrodynamic interactions and enhanced energy harnessing amongst many WEC units in large-size wave parks

: Interactions between wave energy converters (WECs) can significantly affect the overall energy-harnessing performance of a wave park. Although large-size wave parks with many WEC units are commonly considered in practical applications, it is challenging to simulate such parks due to huge computational costs. This paper presents a numerical model that uses the boundary element method (BEM) to simulate wave parks. Each wave energy converter (WEC) was modelled as a comprehensive system, including WEC buoys, power take-off, and mooring systems, with hydrodynamic interactions included. Two classical layouts for arranging 16 units were simulated using this numerical model. The energy-harnessing performance of these array layouts was analyzed for both regular waves and a selection of irregular sea state conditions with different wave directions, wave heights, wave periods and water depths. For each layout, three WEC separation distances were studied. An increase of up to 16% in the power performance of the WEC under regular waves was observed, which highlights the importance of interaction effects.


Introduction
Wave energy is a promising renewable energy source that has attracted increasing attention in the global market [1,2].It shows larger power density and round-the-clock characteristics compared with solar and wind energy [3].However, wave energy converter (WEC) technology has not converged to a certain design like solar panels or wind turbines, which explains the existence of many different kinds of WECs.Falcão [4] classified WECs into oscillating water columns, oscillating bodies and overtopping devices.Different working principles of WECs result in varied requirements for physical tests and numerical modelling.This situation hinders the development of WEC technology.
To harness wave energy cheaply, a typical solution is to install multiple WECs in an area, which constitutes a wave park (also called wave farm or WEC array), to reduce cost by sharing anchors, moorings, power cables, and other subsea installations.The side effect is that, as shown by the pioneering works of Budal and Falnes [5,6], the hydrodynamic interactions caused by diffracted and radiated waves between WECs can be constructive or destructive to the power output depending on the array's layout and sea state conditions, such as the wave height, wave period, wave encounter direction and water depth.In addition to the hydrodynamic interactions, there are also mechanical interactions caused by shared moorings, floaters and mechanically connected WEC components, which also influence the overall power performance of WECs in a wave park [7,8].Moreover, the influence of the interaction effects goes far beyond changes in the power output.Basin tests in [9] indicate a potential increase in mooring loads when installing multiple WECs in an array.
Various studies have been conducted on wave parks to address the influence of the interaction effect and guide wave park designs, with the number of WECs ranging from two to over one thousand.In the most simple wave park, a 2-WEC array, the interaction effects already play an important role in power performance and mooring fatigue damage.In [10], for medium separation distances ranging from 10 to 50 times the diameter of the WEC, the interaction can be low or high depending on the configuration, especially when the wave direction is not aligned with the array.For larger separation distances, such as 75 times the WEC diameter as described in [7], the hydrodynamic interaction has only a weak influence on the calculated power interaction.However, it results in a significant impact on mooring fatigue damage, ranging from an 80% reduction to a fourfold increase.
As the number of WECs increases, the complexity of the effects also builds up.Stallard et al. [11] analysed a 3-WEC rectilinear array, and the results indicate that if the waves align with the array, the performance of the middle WEC can be improved.Krivtsov and Linfoot [9] tested a 5-WEC array and found that the extreme peak mooring loads in the leading mooring line were approximately doubled in similar environmental conditions compared to those in a single device.Göteman [12] simulated a 6-WEC array with 1000 layouts.The study showed that the difference in power absorption between the best and worst layout is 23%, which indicates that the influence of the hydrodynamic interactions is non-negligible, and the interactions can play a positive role with well-designed layouts.Both [11,13] studied 9-WEC arrays with three rows and three columns.The results in [11] showed that the interactions are constructive along the front row but decay with distance through the array; moreover, the interaction factor is highly dependent on the settings of the power take-off (PTO) system.Lamprou et al. [13] varied the wave direction and found that the central WEC within the array has the largest positive interaction under 0-degree waves.However, no clear rules exist for other wave directions about whether the interaction effects are constructive or destructive among various WECs.
To further investigate the potential of wave parks, many studies went beyond 10-WEC arrays.For a 12-WEC array in [14], the interaction factors appear to vary with both incident wave conditions and device location within the array.The authors also found that there was row attenuation across the array in the direction of wave propagation in short-period waves, consistent with [11].Integrating multiple WECs with an offshore wind turbine platform has been a hot research topic in recent years [15][16][17][18][19][20].Interactions between WECs and wind turbine platforms can be seen as a continuation of the conventional WEC-WEC interaction problem, which also involves hydrodynamic and mechanical interaction.A hybrid system was studied in [15] with 24 WECs and four wind turbines.The multiple array configuration adopted by this wind-wave hybrid platform is beneficial when the wave direction is 0 degrees.However, this finding cannot be generalized to other wave directions, which is the same conclusion as in [13].
There are macroscopic studies of wave parks with more than 50 devices.The employment of such large-scale arrays lowers power fluctuation.As found in [21], doubling the number of WECs from 50 to 100 decreases the power variance by approximately 21%; doubling it from 100 to 200 decreases the power variance by approximately 19%.However, for large farms and realistic irregular sea states, the interaction effects are more likely to be destructive [22].
The ways to study interaction effects can be roughly classified into two groups: experiments and numerical simulations.Stallard et al. [11] conducted experiments with arrays of three, nine and twelve identical WECs.Weller et al. [14] studied a 12-WEC array (4 × 3) experimentally.Stratigaki et al. [23] presented their experimental study of a 25-WEC array.Similar experimental studies were conducted by Rahm et al. [24], da Fonseca et al. [25], Wolgamot et al. [26], and Ruiz et al. [27].However, studying interaction effects through experiments has two major drawbacks.Firstly, the results of an experiment usually correspond to a certain set of testing parameters and a specific WEC prototype, so they cannot be generalized to a different set of parameters or a different WEC prototype easily and reliably.This becomes even worse during the wave park design process when parameters are supposed to evolve continually.Moreover, experiments for large-scale wave parks carried out in wave tanks are usually scaled down, which introduces not-well-understood scaling effects.Secondly, the experiments' economic, time, and human resource costs are high.
Numerical simulations are becoming popular in the preliminary design of WECs and wave parks because of their high-speed and low-cost features.The most commonly used approach is the boundary element method (BEM), which has been adopted to analyse the interaction effects in many studies, such as Sinha et al. [28], Tay and Venugopal [29], Yang et al. [7], and Shao et al. [30].HydroD, WAMIT, and ANSYS Aqwa are three software programs based on BEM with a wide audience.The BEM method is based on linear potential theory and cannot consider severe sea conditions.Due to the limitation of the total amount of meshes in BEM software, the number of WECs is therefore restricted.Also, the detailed modelling of PTO and mooring systems is still uncommon because of their complexity and software limitations.The most feasible method for large-scale WEC arrays with over 50 devices is analytical multiple scattering theory, which can consider hydrodynamic interactions based on linear potential flow theory.However, all studies using this method have made compromises to reduce the computational cost.Göteman et al. [22] introduced an interaction cut-off distance, which eliminates the interaction effects between WECs when their distance is beyond a threshold.Göteman et al. [21] included interactions by radiated waves but neglected the diffraction between WECs; their model has not been validated or verified for large arrays, so their results only provide quantitative references for future WEC array designs.
From the above literature review, the challenges faced by wave park designers and researchers can be summarized as follows: (1) lack of general rules to guide the design process without limiting the WEC type, dimensions, and wave conditions.(2) lack of reliable and fast numerical tools that consider the complete system, which includes the WEC buoy, the PTO system and the mooring system.(3) Lack of efficient modelling methods to predict large-size parks with many units considering most of the primary hydrodynamic and mechanical interaction effects.(4) Lack of systematic analysis of various factors influencing a wave park's power performance.Specifically, the challenges (2)-( 4) will be partly addressed in this study.
In this paper, a BEM-based modelling method is developed to push the computation limit up to 16 WEC units within a park.This number has not been achieved for this type of method until now.Compared to the method proposed by Göteman et al. [21,22], the BEM-based modelling method takes into account all interaction effects.To validate this modelling approach, a general sensitivity based on CorPower Ocean's WEC, which is a representative type of bottom-fixed heaving point absorber, is carried out to investigate the power performance of two kinds of 16-WEC array layouts concerning the distance between WEC units, wave conditions and the shape of the arrays.The WEC model takes complete consideration of the buoy's hydrodynamic characteristics, including the viscous drag force and the end-stop feature of the PTO system.This study contributes to the understanding of how the power performance of wave parks is affected by multiple factors and provides valuable reference data for future wave park design.
The remainder of this paper is structured in the following way.Section 2 gives a general description of the theory behind the used method.Section 3 explains the numerical model emphasising the drag elements, the PTO system and the power performance calculation.Section 4 introduces the two 16-WEC layouts and the studied sea state conditions.Section 5 presents the results of the simulations with discussion.Section 6 summarizes the findings and makes conclusions.

Linear Potential Theory and the BEM
Linear potential theory assumes the fluid is inviscid and incompressible while the fluid motion is irrotational.Under this assumption, the fluid can be fully described by the Laplace equation and the Bernoulli equation with prescribed boundary conditions.BEM is a method to solve the velocity potential in the Laplace equation by deploying sources or a mixture of sources, sinks and dipoles over the mean wetted surface.After obtaining the velocity potential, the Bernoulli equation can be applied to obtain the pressure on the wetted surface.The total hydrodynamic forces can be obtained by integrating pressure over the wetted surface.Further details on linear potential theory and BEM can be found in [31,32].The BEM solver HydroD [33] was employed for simulations in the present study.

WEC Hydrodynamics and Interaction Effects
The hydrodynamic forces of a WEC unit are computed in the frequency domain.In this paper, the interaction effects refer to the hydrodynamic interactions caused by the radiation and diffraction between different WECs in the flow field.The total velocity potential for a flow field, including multiple interacting bodies, can be expressed under the framework of linear potential theory in the frequency domain as follows: where φ I and φ d are the velocity potentials of the incident and diffracted waves, φ rjm is the unit radiation potential due to the unit j-th motion of the m-th structure, X is the position coordinates of any point inside the flow field and x jm is the motion amplitude of that structure.The time-varying feature of the potentials is represented by e −iωt , where t is the time, and ω is the wave frequency.The total potential is solved first, and then the corresponding wave excitation force, added mass and damping for different wave frequencies ω are obtained.

WEC Dynamics
The WEC system can be represented as a rigid body dynamic system by integrating hydrodynamic forces, hydrostatic forces, mooring and power take-off (PTO) forces into the equation of motions: where M is the mass matrix and Ẍ is the acceleration.Equation ( 2) can be rewritten in a convolution integral form as follows: where A ∞ is the added mass at infinite frequency, c is the linear damping, K is the total stiffness and R is the retardation function.The hydrodynamic force is divided into two parts.On the left-hand side, the radiation force is represented by A ∞ Ẍ, which resembles an inertia force, and the convolution term The other part is the wave excitation force F excitation , which includes the force due to incident and diffracted waves.F other represents any other forces, for example, the viscous drag force.
The solving process of the WEC dynamics problem is conducted by DNV SESAM SIMA V4.2-00, which enables coupled mooring dynamic analysis and instantaneous hydrostatic force calculation [34].The fully coupled simulation procedure enables the simultaneous solution of the equation of motions for the WECs and the mooring lines.

Numerical Model
3.1.The Bottom-Fixed Heaving Point Absorber WEC Figure 1 illustrates the WEC prototype and the numerical model of CorPower Ocean's WEC.There are two coordinate systems: lowercase letters represent the local coordinates, while capital letters represent the global coordinates.The model consists of a wave buoy with a droplet shape, a stiff mooring line, a seabed anchor node and drag elements.The stiff mooring line provides negative axial stiffness, which amplifies the heave motion of the buoy with a value of −510 kN/m.This value was chosen based on the experiments and tuning of the WEC carried out by CorPower Ocean AB (Hägersten, Sweden).It is called a stiff mooring line because it has high bending stiffness so that the buoy moves only along the mooring line's axial direction.To avoid an extremely large force and violation of the assumption of linear potential theory, which requires small motions, the negative stiffness is set to 0 when the heave motion amplitude in the axial direction of the mooring line exceeds 2.5 m.As the weight of the WEC is far less than the weight of the water displaced by the WEC, the WEC is pre-tensioned by the stiff mooring line to balance extra buoyance.The seabed node is an anchor joint, which allows rotations around the x and y axes.The basic properties of the CorPower WEC system are listed in Table 1.1.4 × 10 6 1 The origin of the reference Cartesian coordinate system (x, y, z) is placed on the plane of the water surface at the geometric centre of the WEC buoy when it is in its unloaded neutral position.

Drag Elements
The drag elements are designed to include the viscous drag force, which is not included in the linear potential flow theory.The viscous force is assumed to be of comparable magnitude with the inertial force since the WEC dimension is close to wave amplitude values.In Figure 2, the location of the drag elements is shown in local coordinates.There are six cylindrical drag elements in total to represent the varying dimensions of the buoy along the vertical direction.The radius and height of each drag element are listed in Table 2.
The drag coefficients of heave and surge/sway, 0.25 and 0.50, respectively, were provided by CorPower Ocean AB through their experience.The viscous forces in surge, sway and heave directions are calculated as follows: where u x , u y and u z are the surge, sway and heave velocities.A 3 is the cross-section area of the drag element D3.A wet,z represents the project area of the wetted element in the z-direction, as follows:

PTO System
The PTO system of the WEC model is present inside the buoy.This PTO system is supposed to travel in an axial direction along the stiff mooring line.The PTO is simplified as a linear damper whose properties are introduced in Section 3.4.Moreover, when the device moves, it is necessary to ensure that the heave motions are restricted within the pre-defined limits to avoid damaging the device.This is taken care of by end-stop bumpers, which are modelled by introducing a non-linear force F es .This force varies according to the following curve: where z es is the end-stop distance beyond which the end-stop bumpers become active, and

Power Performance Calculation
The wave energy harnessed from this WEC unit primarily comes from the heave motions of the buoy.The absorbed hydrodynamic power is calculated as follows: Here, C PTO is the damping coefficient of the PTO, and Żrel is the heave velocity of the buoy in global coordinates.The power output of the WEC unit varies depending on the characteristics of the wave to which it is subjected.Therefore, to determine its hydrodynamic performance characteristics, the power output is calculated for different wave characteristics, represented as sea states in Table 3.The damping coefficient is defined as follows: The term R loss is the linear damping loss, including machinery losses, and is set to a constant value of 25 kNs/m.This value was provided by CorPower Ocean AB based on their experimental data.The term R u represents the useful damping consisting of various damping force elements, which change depending on the sea state the WEC is encountering.The useful damping also accounts for the damping from the generator.The values of R u under different sea state conditions are presented in Table 3.The power output is evaluated in terms of the averaged useful power during the simulation time, which is calculated as follows: In Section 5, the power refers to P u .

Study Cases
This section presents the cases studied in numerical simulations.In Section 4.1, two WEC array layouts are presented, which have been simulated and analysed for various WEC distances and sea state conditions.For all cases with regular wave conditions, the total simulation time was 1000 s; for irregular wave conditions, the total was 3 h.

Single WEC and Park Simulations
The two basic array layouts are shown in Figure 4.Each layout comprises 16 WECs arranged with a separation distance d between the WECs in the horizontal and vertical directions.The 4 × 4 layout has four WECs in each row and four in each column.The 8 × 2 layout consists of eight WECs in each of the two rows.This layout is wider and has more columns than the 4 × 4 layout.The intention behind these two layouts is to keep the same area but vary the number of rows and WECs in each row to analyse their interaction effects on the power generation capacity of WEC arrays.The regular square and rectangular shapes also reduce the potential challenges during the installation of multiple WECs and the deployment of electricity cables.
Table 3 presents the simulated sea state condition sets classified as regular waves (R) and irregular sea states (IR) where the JONSWAP spectrum with a peak enhancement factor of 2.4 was used.The regular waves are divided into sets of low (L) and high (H) waves for periods of 5 to 15 s, denoted by RL 1 to 11 and RH 1 to 11, whereas the three irregular sea states are denoted as IRL 1 to 3.These sea conditions cover the entire range of wave periods and wave heights that are of concern during the operation phase of WECs.The irregular sea states were defined by CorPower Ocean AB and reflect realistic sea state conditions at future installation sites.

Validation of the Single WEC Unit
Figure 5 presents the results from a numerical free decay test simulation.The simulation's purpose was to numerically obtain the WEC's resonance characteristics in heave and surge motions.The heave and surge resonant periods are observed as 8.5 and 17 s, respectively, and are aligned with the resonant periods of the real physical prototype of the WEC measured by CorPower Ocean AB.Therefore, the simulation model is validated in terms of its resonance characteristics.

Effects of Wave Periods and Wave Heights on Power Performance
Figure 6a shows the total useful power output from simulations of the 4 × 4 array layout with a 150 m distance between the WECs.Results from several regular wave conditions are presented for the wave heights 1 m and 4 m.There is a drop in the power output when the wave periods are between 7 and 10 s.This is caused directly by the motion pattern of the WEC, shown in Figure 7.The upward motion range can be up to +4 m, whereas the downward motion can only reach −3 m.A hypothesis for this abnormal motion response observed for all WECs in the array is that it is caused by parametric resonance.The basis for this hypothesis is that parametric resonance is usually observed when the excitation frequency is twice the natural frequency of the system.This is exactly the case as the power drop happens for wave periods between 7 and 10 s, which are around half of the surge resonance period 17 s.The parametric resonance results in a large surge motion, which restricts the range of the translation motion of the WEC along the stiff mooring line.However, validating this hypothesis is beyond the scope of the current study and is left for future work.In addition, the maximum total power output of the array from the simulations is 10.9 MW with a wave height of 4 m and only 0.8 MW with a 1 m wave height.Note that the wave height is increased four times from 1 to 4 m, while the power output is increased more than ten times.This demonstrates the considerable influence of wave height on the power output.An interaction factor was defined as the percentage of the power difference between an array and the same number of single units under the same sea state condition and with other parameters the same to evaluate the influence of interaction effects.The interaction factor of the 4 × 4 array layout is shown in Figure 6b.At the short wave periods of 5 s and 6 s, the power output is enhanced by interaction effects by at least 9% and up to 16%.For longer wave periods, from 12 to 15 s, the influence of interaction on the power performance is less than 1%.This indicates that waves with short wave periods are preferable for enhancing power performance, while waves with long periods have a negligible influence on power performance.

Effects of Water Depth on Power Performance
At the real installation site, the water depth may vary due to the seabed's changing topography, which changes the wave load on the WECs.The rod length that connects the buoy to the seabed should also be changed accordingly.These changes may affect the power performance.Therefore, two cases with 46 m and 55 m were simulated to study the influence of water depth on the power performance.Figure 8 shows the influence of water depth on the power output for the 4 × 4 array layout with a distance of 150 m.The simulations were carried out for the IRL 1 to IRL 3 sea states presented in Table 3.The results show a 1% increase in power performance for the deeper water depth.The reason for this minor positive effect is that a longer mooring line is needed for the 55 m water depth, which enables a larger range of global surge motion and, consequently, a larger range of translation along the stiff mooring line.Overall, the changing water depth has a negligible influence on the interaction effects since the interaction factors are similar for different water depths.

Effects of the Array layout on Power Performance
The influence of the array layout on the power output is presented in Figure 9.The simulations were carried out for the IRL 1 to 3 sea states presented in Table 3.The interaction factor is calculated for each IRL sea state and indicated by the percentage values.The 8 × 2 array layout performs better than the 4 × 4 array layout for IRL 1.This result can be explained by the results presented in Figure 10.For the 8 × 2 array layout with a 0-degree incident wave direction, all eight units in the first (upstream) row experience lower negative interaction effects than the 4 × 4 array layout.However, the second (downstream) row produces less power, probably due to the shielding effects of the first row.There is a significant power output drop for the 8 × 2 array in columns C5 and C6.The possible reason for this phenomenon is that their neighbours, C4 and C7, show large power increases due to positive interaction effects.For the 4 × 4 array layout, the downstream rows after the first-row experience higher negative interaction effects that reduce the motion and power output of the WECs.
Under IRL 2 and 3, the interaction effects are more favourable for the 4 × 4 array layout, while the 8 × 2 array layout experiences destructive interactions at IRL 3 and produces 4% less than the case without interaction.This is caused by the fact that the performance of the WECs in the second and third rows of the 4 × 4 array layout is enhanced by interaction effects while the performance of the first row of the 8 × 2 array layout is greatly reduced, as shown in Figure 11.

Effects of Wave Directions on Power Performance
The influence of the wave direction on the power output for the 4 × 4 array layout is presented in Figure 12.The red line represents the power at the 0-degree wave direction for the IRL 1 to 3 sea states (d = 150 m); the dark blue line represents the same sea states with the wave directions 4.7 degrees, 9.4 degrees and 14.0 degrees, respectively.The trend for the 0-degree case is upward, as is observed in the previous cases for the increased wave period.The percentages represent the increased positive power output due to interaction effects.The interaction is more significant in IRL 2 for the 9.4-degree wave direction as it produces 5% more power than the case without interaction at the 0-degree wave direction for IRL 2. However, normal waves are preferable for IRL 2 and IRL 3 as they enable higher power boosts of 10% and 5%, respectively, compared to their counterparts with incoming wave angles.

Effects of WEC Distances on Power Performance
Figure 13a shows the power output of the 4 × 4 array layout with different distances between the WECs for the RL 1 to 11 simulation cases.The case without interaction does not consider the radiation and diffraction effects from the WEC units in the numerical model.Its power output equals 16 isolated WECs under the corresponding sea state.It is included for comparison purposes in array layout simulations.When the wave period increases from 5 to 7 s, i.e., from RL 1 to 3, the power output is positively related to the wave period.However, from RL 4 to 9, the power output of all distances decreases due to the previously mentioned parametric resonance.Figure 13b shows the interaction factors that share the exact definition as shown in Figure 6b.It can be seen that the 150 m WEC distance is favourable since it enables a relatively large power output increase under RL 1 and 2 and avoids a possible large power output reduction under RL 3 to 11.

Conclusions
A numerical model is put forward to extend the margin of the computational capacity of BEM for wave park simulations.By developing new strategies for the numerical setup, this BEM-based model can model 16 instances at once, while also considering a detailed PTO system.Like the previous development in the family of BEM methods, the present BEM-based model considers all hydrodynamic and structural-dynamic interaction effects.
The bottom-fixed heaving point absorber WEC designed by CorPower Ocean was used to simulate and analyse interaction effects in relation to the array layout wave parks.Several factors varied, such as regular wave and irregular sea state conditions, water depth, wave encounter direction, and WEC distance in the array.
The influence of wave height on the hydrodynamic performance was found to be large.One example of simulation results showed that increased wave height in regular wave conditions with a factor of four increased the power output ten times.Water depth was studied for two cases, 46 m and 55 m, and there was a marginal difference in power output.The latter showed an increase in power output by 1% compared to the former.
The simulation results showed a clear influence of the interaction effects between WECs.They verify what previous researchers have found in experiments and numerical simulations.Each WEC type and array design must be studied carefully for many sea state conditions to find an appropriate distance between the WECs and array layout, to generate power outputs where the interaction factor is beneficial compared to single WEC installations.In this study, for example, the 8 × 2 array layout performed better for IRL 1, while the 4 × 4 array layout performed better for IRL 2 and 3. Regarding the distance between WECs, a 150 m distance seems to be the most favourable if all simulation results are included in the overall assessment since it results in the highest power production in total.

Figure 1 .
(a) The WEC prototype and (b) the numerical model.

Figure 2 .
Figure 2. Arrangement of the drag elements in local coordinates.
is the vertical position in local coordinates.The constants c and n are altered to obtain the desired results.H(•) is the Heaviside step function.The nonlinear stiffness force curves defined with different values of n are shown in Figure 3.The parameter n directly affects the degree of the equation: the larger the value of n, the higher the slope of the stiffness force.The value of n was set to 2 and the c value was set to 1 × 10 8 N/m 2 .The z es used in this WEC is 3.4 m.

Figure 3 .
Figure 3. Non-linear stiffness force defined with different n values for the end-stop damper.

Figure 5 .
Motions of the WEC model in the free decay test simulations: (a) heave and (b) surge.

Figure 6 .
(a) Power output of the 4 × 4 array layout (d = 150 m) for the regular wave conditions RL 1 to 11 and RH 1 to 11; see Table 3 for details.(b) Interaction factors for the cases presented in (a); see the text for details.

Figure 7 .
Figure 7. Illustration of an abnormal heave motion response of a WEC near resonance.Under regular waves, the wave height is 4 m, and the wave period is 10 s.

Figure 8 .
Figure 8. Power output of the 4 × 4 array layout (d = 150 m) for the water depths 46 m and 55 m and the irregular sea states IRL 1 to 3. The percentages are the interaction factors indicating the additional power generation due to interaction effects.

Figure 9 .
Figure 9. Power output of the 4 × 4 and 8 × 2 array layouts (d = 150 m) for the irregular sea states IRL 1 to 3. The percentages are the interaction factors indicating the additional power generation or power reduction due to interaction effects.

Figure 12 .
Figure 12.Power output of the 4 × 4 array layout (d = 150 m) for different wave directions and the irregular sea states IRL 1 to 3. The percentages are the interaction factors indicating the additional power generation due to interaction effects.

Figure 13 .
(a) Power output of the 4 × 4 array layout for different WEC distances for regular wave conditions RL 1 to 11.(b) Interaction factors for different WEC distances under regular wave conditions RL 1 to 11.

Table 2 .
Parameters of the drag elements.

Table 3 .
Regular wave and irregular sea state conditions in the numerical simulations.