2. Wave Energy Resource and the WEC Mechanism
2.1. Wave Energy Characteristics
In recent years, renewable energy sources have been hailed for their potential to lessen the environmental damage caused by the extraction of fossil fuels. In particular, Glendenning (1977)
[68][27] reported that wave energy is blossoming into a viable alternative power source and should be taken into account, due to its greater energy density and reliability than other sources. In order to make the most of this resource and to maximize energy harvesting, it is important to thoroughly study the temporal and spatial wave changes on the nearshore, as suggested by Neil and Hashemi (2013)
[69][28]. Over the last several years, many wave resource evaluations in various coastal regions across the globe have been conducted using the established approaches, such as the numerical modeling of the nearshore wave conditions
[13,70,71,72,73][13][29][30][31][32].
WECs are designed to operate in certain wave climates, and the peculiarities of the sea states also affect how much power they can produce
[74,75][33][34]. Therefore, the differences in the power range of WECs, based on the wave conditions should be evaluated along with the available energy resource to select the best locations for a wave energy plant. However, the majority of the analyses of wave energy resources are restricted to assessing the wave energy resources in a specific coastal region and identifying the sites with the maximum transformation
[76,77][35][36]. In light of these considerations, the innovative approaches for the wave energy evaluation must be created in order to pinpoint more precisely the locations with the most significant potential for collecting wave energy.
The different hydrodynamic modeling methods, such as Delft3D-Wave
[78][37], XBeach
[79][38], WaveWatch III
[80][39], and the standalone Simulating WAve Nearshore (SWAN)
[81[40][41][42],
82,83], are well-verified wind-wave models that may be used to analyze the wave energy resource and identify the nearshore hotspots. The best locations for harvesting wave energy are those where the wave energy is concentrated due to the combined effects of refraction and shoaling. The model inputs include offshore wave characteristics collected at a wave buoy and coastal bathymetry. The nearshore wave simulations using the SWAN model have been performed for some time and are anticipated to improve further.
The Delft University of Technology developed the spectral model SWAN
[83][42], based on third-generation wind-wave modeling, which forecasts the wave characteristics in shallow water. The wave action balancing equation is solved by SWAN using sources and sinks as follows:
where
N=N(σ,θ;x,y,t) denotes the action density;
σ is the relative radian frequency;
θ is the wave direction;
x and
y are the space variables;
t is the time;
cg,x, cg,y, cθ and
cσ are the group velocities in the geographical space and in the spectral space
(σ, θ), respectively.
Stot=Stot(σ,θ;x,y,t) is the sum of the physical processes of the wave energy generation, dissipation, and redistribution:
The parameters of Equation (2) are graphically illustrated in
Figure 1.
Figure 1.
Source terms of wind-wave models.
2.2. Wave Energy Converter Types
The first WECs were installed in 1973 by Salter
[84][43], who investigated the feasibility of wave power.
The amount of literature in this area has grown exponentially since the turn of the century. At present, there are at least six widely recognized categories of conversions, each of which is based on a different combination of factors, including the installation site, the PTO system type, the working regulations, whether the conversion is floating or submerged, and the number of degrees of freedom (DOFs) it possesses. These devices are categorized into different types, such as floating, fully submerged, partly submerged, and devices located on the sea floor. Six DOFs can be broken down into two sections. The first component includes the translational movement, such as heave, sway, and surge, while the second part carries the rotational degrees, such as roll, pitch, and yaw DOFs. WECs may be categorized based on the number of DOFs. The installation site for the conversions is divided into three sections. First, the WECs erected near the coastline or linked to man-made structures, such as breakwaters, are referred to as onshore systems. Second, the nearshore systems are positioned 500 to 2000 m from the coast at a depth of 10 to 25 m. Around a quarter of wavelengths are presented in the nearshore areas. Eventually, the WEC technologies vary from 40 to 100 m respective to the depth of the WEC installations, making adaptive conversions problematic due to the incoming waves
[85][44]. The various WEC type categories are briefly described below.
2.2.1. Point Absorbers
As a kind of WEC device, the point absorbers (PAs) have features and shapes thinner than the incident wavelength. The point absorbers may function in heave, pitch, or multiple DOFs, can be installed onshore or offshore, and can act as one- and two-bodies based on the number of bodies involved. While the latter (two-body) are absolute displacement references about the seabed or offshore structures, the majority of the former (one-body) are self-referenced. The standard shapes of one- and two-body PAs are shown in
Figure 2.
Figure 2.
Common point absorbers.
The PAs with one-body interact with the ocean waves using a single body that is either floating or submerged, and the body’s motion drives a stationary PTO device to generate the energy. The body’s natural frequency and response amplitude operator (RAO), an essential part of the dynamics, shows a low-pass behavior in the frequency domain and significantly impacts the body’s geometrical structure. The two-body point absorbers use two bodies’ relative motion and may operate offshore without requiring mooring devices since they do not need a fixed point of reference. The band-pass properties are unique to the two-body PAs, which allow the band-pass to be tuned to specific wave spectra
[86][45]. The two-body PAs have more complicated geometric forms and need more parameters to specify their geometries than the one-body PAs. Unlike the one-body PAs, which work as low-pass filters, the two-body PAs function as band-pass filters, with the appropriate bandwidths being very sensitive to the geometric forms.
2.2.2. Oscillating Water Columns
An open structure with an open inlet is used in an oscillating water device to trap air above the inner free surface and above the still water level. The trapped air is compressed and released by the wave action in alternating cycles by opening a valve that directs the airflow into a turbine connected to a power plant
[87][46]. The oscillating wave columns (OWCs) do not need an end-stop function since there is no oscillating linear motion present in them. The two types of OWCs are floating and fixed, whereby the fixed OWCs are connected to the breakwaters or the shore
[88][47].
Figure 3 illustrates examples of these OWCs. One significant benefit of fixed OWCs is that other functional needs, such as breakwaters
[89][48], may split the costs for the infrastructure, resulting in a reduced levelized cost of energy (LCoE). In addition, the testing conducted in open water on the permanently installed OWCs proves their reliability and resilience in extreme weather and waves
[24].
Figure 3.
Fixed (
left
) and floating (
right
) OWC.
2.2.3. Attenuators
Parallel to the direction of the waves, the attenuators are free-floating WEC devices made up of a series of bodies linked by hinges, such that the motion of one body relative to another may power a PTO. As seen in
Figure 4, the Pelamis
[90][49] is a prominent example of an attenuator-type WEC.
2.2.4. Terminators
The overtopping technologies, the oscillating wave surge converters (OWSCs), and the “duck-like” systems are all examples of terminators that are positioned perpendicular to the wave direction. Flaps, flaps with vanes, cylinder C-cell shapes, and flaps with arbitrary designs are all used in OWSCs to facilitate the interaction between the wave and the structure. The most important factors to consider while designing these forms are the flap size and depth to which they are submerged. There is some overlap between the categories of the PA and OWSC. However, in this context, the cylindrical devices are considered OWSCs since their working principle and optimization strategy are similar to those of other OWSCs. The overtopping devices, such as the WD and the sea slot-cone generator (SSG), as depicted in
Figure 5, rely on the angle, freeboard and draught, and the ramp shape to function at their best
[91,92,93][50][51][52].
Figure 5.
Terminator types of WECs.
2.3. Power Take-Off Mechanism
The mechanism that transforms the electrical energy generated by the primary conversion is known as the PTO. The direct conversion may be seen in a PA buoy or a water column resonator chamber. The efficiency with which the absorbed wave power is transformed into electricity is directly related to the PTO system’s influence on the wave energy converter’s mass, dimensions, and structural dynamics.
The PTO system has a direct effect on the LCoE
[94][53] since it affects the efficiency with which the wave energy is converted. The PTO system directly influences the effectiveness of the power conversion and, hence, the annual energy generation. The capital cost of a device is directly affected by the PTO system, which accounts for around 20–30% of the total cost
[95][54]. How well the PTO system functions determines on how much it costs to run and maintain, as well as how much energy it produces.
Figure 6 displays a schematic of the PTO’s impact on the LCoE
[96][55]. The PTO variables were studied by the Danish organization Partnership for Wave Energy to determine the impacts of improvements in the PTO efficiency and the PTO framework cost reduction on the LCoE.
Figure 6.
Parameters defining the LCoE for the WECs.
A decrease in the LCoE is associated with both an efficiency improvement and the PTO cost reduction. Although a boost in the PTO performance has a more noticeable impact on the LCoE, both variables significantly impact the LCoE, indicating the importance of the PTO system in a wave energy conversion.
3. Recent Advances in Optimizing the WEC Configuration
3.1. Optimization Approaches
Converging to the best solution directly in most actual engineering issues is difficult. Yet, putting up a loss function that gauge’s the solution quality and decreases its parameters to discover the best possible solution is not easy. Some of the research initiatives explore numerical approaches. This section thoroughly discusses several optimization methods applied in the WEC systematic design.
Noad and Porter
[106][56] used a multi-dimensional numerical optimization to assess the solution method’s correctness and to determine the appropriate device settings for the array optimization. The design of the WEC arrays was described by Ruiz et al.
[55][57], using four parameters, as well as looking at the GA, the glowworm swarm optimization (GSO), and the CMA-ES. In contrast to the GA and GSO, the CMA requires less processing power, as found by the researchers. Raju
[107][58] employed the Nelder–Mead (NM) search method to determine the best PTO parameters for each WEC and the derivative-free continuous optimization, the CMA-ES, to enhance the conversion location by minimizing the negative interactions.
Thomas et al.
[108][59] presented a shallow artificial neural network (ANN) to find the ideal latching times in irregular wave situations. According to the results, for particular sea conditions, the learnable WEC absorbs twice as much power as the WEC without latching and 30% more power when compared to the test wave’s best consistent latching time. Complex PTO control methods, such as the latching-declutching optimization, aim to improve the overall system performance. Due to the discontinuous nature of the loss function in this optimization issue, Feng and Kerrigan
[109][60] opted to adopt a novel derivative-free coordinate-search approach, developing a formulation based on historical wave data and forecasting future wave behavior. To demonstrate the algorithm’s effectiveness, it was compared to the simulated annealing (SA) algorithm, a global meta-heuristic technique without derivatives.
The bio-inspired algorithms use simple, naturally derived techniques to address complicated issues. Since many biological processes may be considered a local optimization, these techniques imitate nature. They fall under the scope of randomized algorithms and extensively use arbitrary choices. This method has gained much popularity because it can be used to tackle complicated issues in all major branches of computer science and offers solutions to a wide variety of issues. Bio-inspired algorithms have been divided into many categories, including ecology-based, multi-objective, many-objective, and EAs.
3.2. Layout-PTO-Geometry Optimization
Numerous studies investigated the prospect of enhancing the array’s capacity to absorb energy by considering various designs. Such arrangements may result in a 30% reduction in the array or even a growth of around 5%. Over the last decade, many researchers have examined the uncomplicated regular, and unpredictable patterns and found that the optimized patterns may be a solid option to impact extracting power positively.
A number of variables influence the placement of buoys in an array, which could be beneficial or destructive depending on the project’s ultimate objective. In particular, the closer WECs are to one another, the greater the adverse effects when less expensive cabling is used. The six most important elements in designing an array are discussed as follows. First, the number of WECs profoundly affects the overall design, whereby more than four transformations are needed for the polygonal and circular designs. Second, the array’s behavior may shift if each conversion is handled separately. Therefore, the power absorbed by an array would be the same as the power of the isolated powers, and the interaction effects of the array may even be negligible. Third, the power absorption characteristic is very sensitive to sea conditions. As a result, the wave interaction effects are applicable, and a lot of energy is consumed when the natural period of the system is similar to the wave period. In the case of regular waves, the wave’s direction and frequency also directly affect the excitation force. The interaction factor is susceptible to the strength of the excitation force used on the device. Fourth, the direction of the waves may have a significant role in determining the design. This may be tested in two ways: theoretically by rotating the array pattern or experimentally by employing different definitions at different points in the study. Despite the fact that most publications solely take into account unidirectional waves, the correlation of multidirectional waves also needs consideration. Fifth, the amount of energy harvested from each WEC is affected by its individual size. Finally, it is essential to consider how waves interact with the buoys when placing them to minimize the detrimental effects.
Table 1 presents the examples of the most recent literature that meet the requirements for selecting the optimal design of a WEC array, as part of their configurations, numbers of WECs, and types of conversion.
Table 1.
Summary of the recent publications on the WEC array optimization.
Author(s)–Year |
WEC Type |
WEC No. |
Objective Function |
Algorithm |
Ref. |
PA |
OWC |
Attenuator |
Terminator |
Deandres (2014) |
✓ |
|
|
|
2, 3, 4 |
q-factor |
GA + Parabolic Intersection |
[128][61] |
Baltisky (2014) |
|
|
✓ |
|
2, 3, 4, 5, 6 |
Mean AEP |
Global Control |
[129][62] |
Noad et al. (2015) |
|
|
|
✓ |
3, 5 |
Absorbed power |
Multi-Dimensional Optimization |
[106][56] |
Blanco et al. (2015) |
✓ |
|
|
|
2 |
Maximize power |
EA |
[130][63] |
Sharp & DuPont (2016) |
✓ |
|
|
|
5 |
Power and cost |
GA |
[131][64] |
Wu et al. (2016) |
✓ |
|
|
|
25, 50, 100 |
Computational cost |
EA and CMA-ES |
[54][65] |
Sarkar et al. (2016) |
|
|
|
✓ |
40 |
Maximize power |
GA and Monte Carlo |
[42][66] |
Ruiz et al. (2017) |
|
✓ |
|
|
>10 |
Maximize power |
CMA-ES, GA, GSO |
[55][57] |
Ferri (2017) |
✓ |
|
|
|
>50 |
Computational cost |
CMA-ES + MM |
[132][67] |
Giassi et al. (2017) |
✓ |
|
|
|
9, 12 |
Maximize power |
GA |
[112][68] |
Blanco et al. (2018) |
|
|
✓ |
|
2 |
Maximize power |
EA |
[119][69] |
Sharp and DuPont (2018) |
|
|
|
✓ |
5 |
Power and cost |
GA |
[113][70] |
Giassi et al. (2018) |
✓ |
|
|
|
4–14 |
Maximize power |
GA + Multiple Scattering |
[56][71] |
Fang et al. (2018) |
✓ |
|
|
|
3, 5, 8 |
Maximize power |
EA |
[133][72] |
Neshat et al. (2018) |
✓ |
|
|
|
4 |
Maximize power |
Meta-Heuristic Algorithm |
[134][73] |
Lyu et al. (2019) |
✓ |
|
|
|
3, 5, 7 |
Optimal control |
GA |
[110][74] |
Vatchavayi (2019) |
✓ |
|
|
|
4–9 |
Maximize power |
CMA-ES |
[107][58] |
Neshat et al. (2019) |
✓ |
|
|
|
16 |
Maximize power |
Neural Optimization + Analytical |
[135][75] |
Faraggiana et al. (2019) |
✓ |
|
|
|
1–3 |
Minimize LCoE |
PSO and GA |
[121][76] |
Neshat et al. (2020) |
✓ |
|
|
|
49, 100 |
Maximize power |
Multi-Strategy EAs |
[136][77] |
Neshat et al. (2020) |
✓ |
|
|
|
4, 16 |
Maximize power |
Cooperative EAs |
[137][78] |
Bosma et al. (2020) |
|
✓ |
|
|
5 |
Average power |
- |
[138][79] |