Difference between revisions of "Wave energy converters in coastal structures"
(→Application for wave energy converters) 
Dronkers J (talk  contribs) 

(26 intermediate revisions by 3 users not shown)  
Line 1:  Line 1:  
−  +  {style="float:right; marginleft: 1em; textalign:center; fontsize:85%"  
−  
−  {  
[[Image:Construction_of_a_coastal_structure.jpg300px]]  [[Image:Construction_of_a_coastal_structure.jpg300px]]  
    
'''Fig 1: Construction of a coastal structure.'''  '''Fig 1: Construction of a coastal structure.'''  
}  }  
+  == Introduction ==  
+  <p>  
Coastal works along European coasts are composed of very diverse structures. Many coastal structures are ageing and facing problems of stability, sustainability and erosion. Moreover climate change and especially sea level rise represent a new danger for them. [[Sea dikesCoastal dykes]] in Europe will indeed be exposed to [[waves]] with [[Wave heightheights]] that are greater than the dykes were designed to withstand, in particular all the structures built in shallow water where the depth imposes the maximal amplitude because of wave breaking.  Coastal works along European coasts are composed of very diverse structures. Many coastal structures are ageing and facing problems of stability, sustainability and erosion. Moreover climate change and especially sea level rise represent a new danger for them. [[Sea dikesCoastal dykes]] in Europe will indeed be exposed to [[waves]] with [[Wave heightheights]] that are greater than the dykes were designed to withstand, in particular all the structures built in shallow water where the depth imposes the maximal amplitude because of wave breaking.  
<p>  <p>  
Line 14:  Line 14:  
<br>  <br>  
<p>  <p>  
−  
−  
== Wave energy and wave energy flux ==  == Wave energy and wave energy flux ==  
<p>  <p>  
−  For regular water waves, the timemean [[waveswave]] energy density E per unit horizontal area on the water surface  +  For regular water waves, the timemean [[waveswave]] energy density <math>E</math> per unit horizontal area on the water surface [J/m<sup>2</sup>] is the sum of kinetic and potential energy density per unit horizontal area. The potential energy density is equal to the kinetic energy <ref name="ref1">Mei C.C. (1989) The applied dynamics of ocean surface waves. Advanced series on ocean engineering. World Scientific Publishing Ltd </ref> both contributing half to the timemean wave energy density E that is proportional to the wave height squared according to linear wave theory <ref name="ref1"/> (see [[Shallowwater wave theory]]): 
<p>  <p>  
<br>  <br>  
Line 25:  Line 23:  
(1)  (1)  
</div>  </div>  
−  <math>E= \frac{1}{8} \rho g H^2</math>  +  <math>E= \large\frac{1}{8}\normalsize \rho g H^2</math> 
</div>  </div>  
−  
<br>  <br>  
<p>  <p>  
−  g is the gravity and <math>H</math> the wave height of regular water waves. As the waves propagate, their energy is transported. The energy transport velocity is the group velocity. As a result, the timemean wave energy flux per unit crest length  +  <math>g</math> is the gravity acceleration and <math>H</math> the wave height of regular water waves. As the waves propagate, their energy is transported. The energy transport velocity is the group velocity. As a result, the timemean wave energy flux per unit crest length [W/m] perpendicular to the wave propagation direction, is equal to <ref name="ref1"/>: 
<p>  <p>  
<br>  <br>  
Line 38:  Line 35:  
(2)  (2)  
</div>  </div>  
−  <math> P=  +  <math> P = E \times c_{g}</math> 
</div>  </div>  
−  
<br>  <br>  
<p>  <p>  
−  with <math>c_{g}</math> the group velocity  +  with <math>c_{g}</math> the group velocity [m/s]. Due to the dispersion relation for water waves under the action of gravity, the group velocity depends on the wavelength <math>\lambda</math> [m], or equivalently, on the wave period <math>T </math>[s]. Further, the dispersion relation is a function of the water depth <math>h</math> [m]. As a result, the group velocity behaves differently in the limits of deep and shallow water, and at intermediate depths: 
<div style="textalign: center;">  <div style="textalign: center;">  
−  <math>(\frac{\lambda}{20} < h < \frac{\lambda}{2})</math>  +  <math>(\large\frac{\lambda}{20}\normalsize < h < \large\frac{\lambda}{2})\normalsize </math> 
</div>  </div>  
−  
<br>  <br>  
<p>  <p>  
Line 55:  Line 50:  
For regular waves in deep water:<br><div style="textalign: center;">  For regular waves in deep water:<br><div style="textalign: center;">  
<div style="float: right">(3)</div>  <div style="float: right">(3)</div>  
−  <math>c_{g} = \frac{gT}{4\pi} </math> and <math>P_{w1} = \frac{\rho g^2}{32 \pi} H^2 T</math>  +  <math>c_{g} = \large\frac{gT}{4\pi}\normalsize </math> and <math>P_{w1} = \large\frac{\rho g^2}{32 \pi}\normalsize H^2 T</math> 
</div>  </div>  
<p>  <p>  
<br>  <br>  
−  
The timemean wave energy flux per unit crest length is used as one of the main criteria to choose a site for wave energy converters.  The timemean wave energy flux per unit crest length is used as one of the main criteria to choose a site for wave energy converters.  
−  
<p>  <p>  
<br>  <br>  
−  +  For real seas, whose waves are random in height, period (and direction), the spectral parameters have to be used (for the definition of spectral wave parameters, see the article [[Statistical description of wave parameters]]). The spectral estimate of significant wave height <math>H_{m0} </math> is based on zeroorder moment of the spectral function as <math>H_{m0} = 4 \sqrt{m_0} </math>. Moreover the wave period is derived as follows <ref name ="ref2"> Vicinanza D., Cappietti L., Ferrante V. and Contestabile P. (2011) : Estimation of the wave energy along the Italian offshore, journal of coastal research, special issue 64, pp 613  617. </ref>.  
−  
−  For real seas, whose waves are random in height, period (and direction), the spectral parameters have to be used. <math>H_{m0} </math>  
<p>  <p>  
<br>  <br>  
<p><div style="textalign: center;">  <p><div style="textalign: center;">  
<div style="float: right">(4)</div>  <div style="float: right">(4)</div>  
−  <math>  +  <math>T_E = \large\frac{m_{1}}{m_0}\normalsize 
</math>  </math>  
</div>  </div>  
<p>  <p>  
<br>  <br>  
−  +  where <math>m_n</math> represents the spectral moment of order <math>n</math>. An equation similar to that describing the power of regular waves is then obtained <ref name="ref2"/> :  
−  
−  where <math>m_n</math>  
−  
−  
<p>  <p>  
<br>  <br>  
−  
−  
<br><div style="textalign: center;">  <br><div style="textalign: center;">  
<div style="float: right">(5)</div>  <div style="float: right">(5)</div>  
−  <math>P_{w1} = \frac{\rho g^2}{64 \pi} H_{m0}^2  +  <math>P_{w1} = \large\frac{\rho g^2}{64 \pi}\normalsize H_{m0}^2 T_E</math> 
</div>  </div>  
<p>  <p>  
<br>  <br>  
−  +  If local data are available (<math>H_{m0}^2, T_E </math>) for a sea state through insitu wave buoys for example, satellite data or numerical modelling, the last equation giving wave energy flux <math>P_{w1}</math> gives a first estimation. Averaged over a season or a year, it represents the maximal energetic resource that can be theoretically extracted from wave energy.  
−  +  If the directional spectrum of sea state variance <math>F (f,\theta)</math> is known with <math>f</math> the wave frequency [Hz] and <math>\theta</math> the wave direction [rad], a more accurate formulation is used:  
−  
−  If local data are available (<math>H_{m0}^2,  
−  If the directional spectrum of sea state variance F (f,  
<p>  <p>  
<br>  <br>  
Line 104:  Line 86:  
<p>  <p>  
<br>  <br>  
−  +  {style="float:right; marginright: 1em; textalign:center; fontsize:85%"  
−  {  
[[Image:Timemean wave energy flux along West European coasts.jpg280px]]  [[Image:Timemean wave energy flux along West European coasts.jpg280px]]  
    
Line 113:  Line 94:  
<p>  <p>  
−  From TOMAWAC model, the near shore wave atlas ANEMOC along the coasts of Europe and France based on the numerical modelling of wave climate over 25 years has been produced <ref name="Ref 2">Benoit M. and Lafon F. (2004) : A nearshore wave atlas along the coasts of France based on the numerical modeling of wave climate over 25 years, 29th International Conference on Coastal Engineering (ICCE’2004), Lisbonne (Portugal), pp 714726.</ref>. Using equation (  +  From TOMAWAC model, the near shore wave atlas ANEMOC along the coasts of Europe and France based on the numerical modelling of wave climate over 25 years has been produced <ref name="Ref 2">Benoit M. and Lafon F. (2004) : A nearshore wave atlas along the coasts of France based on the numerical modeling of wave climate over 25 years, 29th International Conference on Coastal Engineering (ICCE’2004), Lisbonne (Portugal), pp 714726.</ref>. Using equation (6), the timemean wave energy flux along West European coasts is obtained (see Fig. 2). This equation (6) still presents some limits like the definition of the bounds of the integration. Moreover, the objective to get data on the wave energy near coastal structures in shallow or intermediate water requires the use of numerical models that are able to represent the physical processes of wave propagation like the refraction, shoaling, dissipation by bottom friction or by wave breaking, interactions with tides and diffraction by islands. 
<p>  <p>  
The wave energy flux is therefore calculated usually for water depth superior to 20 m. This maximal energetic resource calculated in deep water will be limited in the coastal zone:  The wave energy flux is therefore calculated usually for water depth superior to 20 m. This maximal energetic resource calculated in deep water will be limited in the coastal zone:  
Line 131:  Line 112:  
# the converter survivability : the capacity of the converter to survive in extreme conditions. The survivability gives an estimate of the cost considering that the weaker are the extreme efforts in comparison with the mean effort, the smaller is the cost.  # the converter survivability : the capacity of the converter to survive in extreme conditions. The survivability gives an estimate of the cost considering that the weaker are the extreme efforts in comparison with the mean effort, the smaller is the cost.  
<p>  <p>  
−  Unfortunately, few data are available in literature. In order to determine the characteristics of the different wave energy technologies, it is necessary to class them first in four main families <ref name=  +  Unfortunately, few data are available in literature. In order to determine the characteristics of the different wave energy technologies, it is necessary to class them first in four main families <ref name=ref4/>. 
+  
−  { style ="marginleft: 1em  +  { style="marginleft: 1em; marginleft:auto; marginright:auto; width:650px; textalign:center; fontsize:85%" 
Valign="top" [[Image:Overtopping_with_lowhead_hydraulic_turbine.jpg250px]]  Valign="top" [[Image:Overtopping_with_lowhead_hydraulic_turbine.jpg250px]]  
Valign="top" Colspan="2"[[Image:Submerged_oscillating_bodies_with_hydraulic_motor%2C_hydraulic_turbine%2C_linear_electrical_generator.jpg350px]]  Valign="top" Colspan="2"[[Image:Submerged_oscillating_bodies_with_hydraulic_motor%2C_hydraulic_turbine%2C_linear_electrical_generator.jpg350px]]  
    
−  Valign="top" Overtopping with lowhead hydraulic turbine  +  Valign="top" Overtopping devices (OD) with lowhead hydraulic turbine. 
−  Valign="top" Colspan="2"Submerged oscillating bodies with hydraulic motor, <br> hydraulic turbine, linear electrical generator.  +  Valign="top" Colspan="2"Submerged oscillating bodies (SOB) with hydraulic motor, <br> hydraulic turbine, linear electrical generator. 
    
Valign="top"[[Image:Floating oscillating bodies with hydraulic motor, hydraulic turbine.jpg250px]]  Valign="top"[[Image:Floating oscillating bodies with hydraulic motor, hydraulic turbine.jpg250px]]  
Line 144:  Line 126:  
[[Image:Oscillating water column with water turbine 2.jpg180px]]  [[Image:Oscillating water column with water turbine 2.jpg180px]]  
    
−  Valign="top"  +  Valign="top" Point absorbers (heaving buoys, PA) with hydraulic <br> motor, hydraulic turbine, linear electrical generator. 
−  Valign="top" colspan="2"Oscillating water column with water turbine.  +  Valign="top" colspan="2"Oscillating water column (OWC) with water turbine. 
    
−  colspan ="3" '''Fig 3: The various wave energy technologies'''  +  colspan ="3" '''Fig 3: The various wave energy technologies.'''<p> 
}  }  
−  
<br>  <br>  
<p>  <p>  
−  An interesting result is that the maximum average wave power that a point absorber can absorb <math>P_{abs} </math>  +  An interesting result is that the maximum average wave power that a point absorber can absorb <math>P_{abs} </math>[W] from the waves does not depend on its dimensions <ref name ="ref4">De O. Falcão A. F. (2010) Wave energy utilization: A review of the technologies. Renewable and Sustainable Energy Reviews, Volume 14, Issue 3, April 2010, pp. 899–918. </ref>. It is theoretically possible to absorb a lot of energy with only a small buoy. It can be shown that for a body with a vertical axis of symmetry (but otherwise arbitrary geometry) oscillating in heave the capture (or absorption) width <math>L_{max}</math>[m] is as follows <ref name =ref4/>: 
−  
<br>  <br>  
<p>  <p>  
<div style="textalign: center;"><div style="float: right">(7)</div>  <div style="textalign: center;"><div style="float: right">(7)</div>  
−  <math>L_{max} = \frac{P_{abs}}{P_{w}} = \frac{\lambda}{2\pi}</math> or <math>1 = \frac{P_{abs}}{P_{w}} \frac{2\pi}{\lambda}</math>  +  <math>L_{max} = \large\frac{P_{abs}}{P_{w}} \normalsize = \large\frac{\lambda}{2\pi}\normalsize </math> or <math>1 = \large\frac{P_{abs}}{P_{w}} \frac{2\pi}{\lambda}\normalsize </math> 
</div>  </div>  
−  
<br>  <br>  
<p>  <p>  
−  {  +  {style="float:right; marginleft: 1em; textalign:center; fontsize:85%" 
[[Image:Upper limit of mean wave power absorption for a heaving point absorber.jpg300px]]  [[Image:Upper limit of mean wave power absorption for a heaving point absorber.jpg300px]]  
    
'''Fig 4: Upper limit of mean wave power <br>absorption for a heaving point absorber.'''  '''Fig 4: Upper limit of mean wave power <br>absorption for a heaving point absorber.'''  
}  }  
−  where <math>{P_{w}}</math> is the wave energy flux per unit crest length  +  where <math>{P_{w}}</math> is the wave energy flux per unit crest length [W/m]. An optimally damped buoy responds however efficiently to a relatively narrow band of wave periods. 
<p>  <p>  
−  Babarit  +  Babarit and Hals propose <ref name="ref5">Babarit A. and Hals J. (2011) On the maximum and actual capture width ratio of wave energy converters – 11th European Wave and Tidal Energy Conference Series (EWTEC’2011) – Southampton (UK).</ref> to derive that upper limit for the mean annual power in irregular waves at some typical locations where one could be interested in putting some wave energy devices. The mean annual power absorption tends to increase linearly with the wave power resource. Overall, one can say that for a typical site whose resource is between 2030 kW/m, the upper limit of mean wave power absorption is about 1 MW for a heaving WEC with a capture width between 3050 m. 
−  
−  
<p>  <p>  
+  In order to complete these theoretical results and to describe the efficiency of the WEC in practical situations, the capture width ratio <math>\eta</math> is also usually introduced. It is defined as the ratio between the absorbed power and the available wave power resource per meter of wave front times a relevant dimension <math>B</math> [m].  
<br>  <br>  
<p>  <p>  
<div style="textalign: center;"><div style="float: right">(8)</div>  <div style="textalign: center;"><div style="float: right">(8)</div>  
−  <math>\eta = \frac{P_{abs}}{P_{w}B} </math>  +  <math>\eta = \large\frac{P_{abs}}{P_{w}B}\normalsize </math> 
</div>  </div>  
<p>  <p>  
<br>  <br>  
−  +  The choice of the dimension <math>B</math> will depend on the working principle of the WEC. Most of the time, it should be chosen as the width of the device, but in some cases another dimension is more relevant. Estimations of this ratio <math>\eta</math> for different systems (see Fig. 3) are given <ref name = "ref5"/>: 33 % for oscillating water columns (OWC), 13 % for overtopping devices (OD), 929 % for point absorbers (heaving buoys, PA), 2041 % for submerged oscillating bodies (pitching devices, SOB). For energy converted to electricity, one must take into account moreover the energy losses in other components of the system.  
−  The choice of the dimension B will depend on the working principle of the WEC. Most of the time, it should be chosen as the width of the device, but in some cases another dimension is more relevant. Estimations of this ratio <math>\eta</math> are given <ref name = "ref5"/>: 33 % for OWC, 13 % for overtopping devices, 929 % for heaving buoys, 2041 % for pitching devices. For energy converted to electricity, one must take into account moreover the energy losses in other components of the system.  
<p>  <p>  
<br>  <br>  
Line 194:  Line 171:  
<br>  <br>  
<p>  <p>  
−  { style ="marginleft: 1em  +  { style ="marginleft: 1em" style="margin:0 auto; width:700px; textalign:center; fontsize:85%" 
[[Image:Oscillating water column configuration 1.jpg325px]]  [[Image:Oscillating water column configuration 1.jpg325px]]  
[[Image:Oscillating water column configuration 2.jpg325px]]  [[Image:Oscillating water column configuration 2.jpg325px]]  
Line 213:  Line 190:  
<br>  <br>  
<p>  <p>  
−  {  +  {style="float:right; marginleft: 1em; textalign:center; fontsize:85%" 
[[Image:Finistere area and locations of the six sites.jpg250px]]  [[Image:Finistere area and locations of the six sites.jpg250px]]  
    
Line 226:  Line 203:  
<p>  <p>  
<br style="clear:both;"/>  <br style="clear:both;"/>  
−  { style ="marginleft: 1em  +  { style ="marginleft:1em; textalign:center; fontsize:85%; width:700px; margin:0 auto" 
[[Image:Satelite_picture_Roskof.jpg219px]]  [[Image:Satelite_picture_Roskof.jpg219px]]  
[[Image:Satelite_picture_Molene.jpg225px]]  [[Image:Satelite_picture_Molene.jpg225px]]  
Line 252:  Line 229:  
<br>  <br>  
<p>  <p>  
−  { style ="marginleft: 1em  +  { style ="marginleft: 1em; margin:0 auto; width: 700px; textalign:center; fontsize:85%" 
[[Image:Waves measures at the Pierres Noires Lighthouse 1.jpg225px]]  [[Image:Waves measures at the Pierres Noires Lighthouse 1.jpg225px]]  
[[Image:Waves measures at the Pierres Noires Lighthouse 2.jpg435px]]  [[Image:Waves measures at the Pierres Noires Lighthouse 2.jpg435px]]  
Line 279:  Line 256:  
<p>  <p>  
== See also ==  == See also ==  
−  
−  
−  
−  
−  
−  
−  
* [[Wave energy converters]]  * [[Wave energy converters]]  
−  * [[  +  * [[Shallowwater wave theory]] 
−  * [[  +  * [[Hard coastal protection structures]] 
<p>  <p>  
<br>  <br>  
Line 300:  Line 270:  
<div align="center">  <div align="center">  
{ style="border:1px solid #abd5f5; background:#f1f5fc; margin:0em 0em 0em 0em;"  { style="border:1px solid #abd5f5; background:#f1f5fc; margin:0em 0em 0em 0em;"  
−  <center>The main author of this article is [http://www.coastalwiki.org/index.php?option=com_imis&module=person&Itemid=17&persid=11176 Sergent, Philippe]<br><small>With contributions by: François Bouttes, Bertrand Michard, Emmanuel Cosquer, Alain Clément, Aurélien Babarit, Virginie Baudry, Michel Benoit and Giovanni Mattarolo</small></center>  +  <center>The main author of this article is [http://www.coastalwiki.org/index.php?option=com_imis&module=person&Itemid=17&persid=11176 Sergent, Philippe]<br><small>With contributions by: François Bouttes, Bertrand Michard, Emmanuel Cosquer, Alain Clément, <br> Aurélien Babarit, Virginie Baudry, Michel Benoit and Giovanni Mattarolo</small></center> 
    
Line 306:  Line 276:  
}</div>  }</div>  
−  [[Category:  +  [[Category:Hard structures]] 
−  [[Category:  +  [[Category:Hydrodynamics]] 
−  
−  
−  
− 
Latest revision as of 16:32, 5 June 2020
Fig 1: Construction of a coastal structure. 
Contents
Introduction
Coastal works along European coasts are composed of very diverse structures. Many coastal structures are ageing and facing problems of stability, sustainability and erosion. Moreover climate change and especially sea level rise represent a new danger for them. Coastal dykes in Europe will indeed be exposed to waves with heights that are greater than the dykes were designed to withstand, in particular all the structures built in shallow water where the depth imposes the maximal amplitude because of wave breaking.
These structures need therefore to be modernized and adapted to climate change on one hand and to increase of maritime traffic and size of container carriers on the other hand.
This necessary adaptation will be costly but will provide an opportunity to integrate converters of sustainable energy in the new maritime structures along the coasts and in particular in harbours. This initiative will contribute to the reduction of the greenhouse effect. Produced energy can be directly used for the energy consumption in harbour area and will reduce the carbon footprint of harbours by feeding the docked ships with green energy. Nowadays these ships use their motors to produce electricity power on board even if they are docked. Integration of wave energy converters (WEC) in coastal structures will favour the emergence of the new concept of future harbours with zero emissions.
Wave energy and wave energy flux
For regular water waves, the timemean wave energy density [math]E[/math] per unit horizontal area on the water surface [J/m^{2}] is the sum of kinetic and potential energy density per unit horizontal area. The potential energy density is equal to the kinetic energy ^{[1]} both contributing half to the timemean wave energy density E that is proportional to the wave height squared according to linear wave theory ^{[1]} (see Shallowwater wave theory):
(1)
[math]E= \large\frac{1}{8}\normalsize \rho g H^2[/math]
[math]g[/math] is the gravity acceleration and [math]H[/math] the wave height of regular water waves. As the waves propagate, their energy is transported. The energy transport velocity is the group velocity. As a result, the timemean wave energy flux per unit crest length [W/m] perpendicular to the wave propagation direction, is equal to ^{[1]}:
(2)
[math] P = E \times c_{g}[/math]
with [math]c_{g}[/math] the group velocity [m/s]. Due to the dispersion relation for water waves under the action of gravity, the group velocity depends on the wavelength [math]\lambda[/math] [m], or equivalently, on the wave period [math]T [/math][s]. Further, the dispersion relation is a function of the water depth [math]h[/math] [m]. As a result, the group velocity behaves differently in the limits of deep and shallow water, and at intermediate depths:
[math](\large\frac{\lambda}{20}\normalsize \lt h \lt \large\frac{\lambda}{2})\normalsize [/math]
Application for wave energy converters
For regular waves in deep water:[math]c_{g} = \large\frac{gT}{4\pi}\normalsize [/math] and [math]P_{w1} = \large\frac{\rho g^2}{32 \pi}\normalsize H^2 T[/math]
The timemean wave energy flux per unit crest length is used as one of the main criteria to choose a site for wave energy converters.
For real seas, whose waves are random in height, period (and direction), the spectral parameters have to be used (for the definition of spectral wave parameters, see the article Statistical description of wave parameters). The spectral estimate of significant wave height [math]H_{m0} [/math] is based on zeroorder moment of the spectral function as [math]H_{m0} = 4 \sqrt{m_0} [/math]. Moreover the wave period is derived as follows ^{[2]}.
[math]T_E = \large\frac{m_{1}}{m_0}\normalsize [/math]
where [math]m_n[/math] represents the spectral moment of order [math]n[/math]. An equation similar to that describing the power of regular waves is then obtained ^{[2]} :
[math]P_{w1} = \large\frac{\rho g^2}{64 \pi}\normalsize H_{m0}^2 T_E[/math]
If local data are available ([math]H_{m0}^2, T_E [/math]) for a sea state through insitu wave buoys for example, satellite data or numerical modelling, the last equation giving wave energy flux [math]P_{w1}[/math] gives a first estimation. Averaged over a season or a year, it represents the maximal energetic resource that can be theoretically extracted from wave energy.
If the directional spectrum of sea state variance [math]F (f,\theta)[/math] is known with [math]f[/math] the wave frequency [Hz] and [math]\theta[/math] the wave direction [rad], a more accurate formulation is used:
[math]P_{w2} = \rho g\int\int c_{g}(f,h)F(f,\theta) dfd \theta[/math]
Fig 2: Timemean wave energy flux along West European coasts ^{[3]}. 
It can be shown easily that equations (5 and 6) can be reduced to (4) with the hypothesis of regular waves in deep water. The directional spectrum is deduced from directional wave buoys, SAR images or advanced spectral windwave models, known as thirdgeneration models, such as WAM, WAVEWATCH III, TOMAWAC or SWAN. These models solve the spectral action balance equation without any a priori restrictions on the spectrum for the evolution of wave growth.
From TOMAWAC model, the near shore wave atlas ANEMOC along the coasts of Europe and France based on the numerical modelling of wave climate over 25 years has been produced ^{[4]}. Using equation (6), the timemean wave energy flux along West European coasts is obtained (see Fig. 2). This equation (6) still presents some limits like the definition of the bounds of the integration. Moreover, the objective to get data on the wave energy near coastal structures in shallow or intermediate water requires the use of numerical models that are able to represent the physical processes of wave propagation like the refraction, shoaling, dissipation by bottom friction or by wave breaking, interactions with tides and diffraction by islands.
The wave energy flux is therefore calculated usually for water depth superior to 20 m. This maximal energetic resource calculated in deep water will be limited in the coastal zone:
 at low tide by wave breaking;
 at high tide in storm event when the wave height exceeds the maximal operating conditions;
 by screen effect due to the presence of capes, spits, reefs, islands,...
Technologies
According to the International Energy Agency (IEA), more than hundred systems of wave energy conversion are in development in the world. Among them, many can be integrated in coastal structures. Evaluations based on objective criteria are necessary in order to sort theses systems and to determine the most promising solutions.
Criteria are in particular:
 the converter efficiency : the aim is to estimate the energy produced by the converter. The efficiency gives an estimate of the number of kWh that is produced by the machine but not the cost.
 the converter survivability : the capacity of the converter to survive in extreme conditions. The survivability gives an estimate of the cost considering that the weaker are the extreme efforts in comparison with the mean effort, the smaller is the cost.
Unfortunately, few data are available in literature. In order to determine the characteristics of the different wave energy technologies, it is necessary to class them first in four main families ^{[5]}.
An interesting result is that the maximum average wave power that a point absorber can absorb [math]P_{abs} [/math][W] from the waves does not depend on its dimensions ^{[5]}. It is theoretically possible to absorb a lot of energy with only a small buoy. It can be shown that for a body with a vertical axis of symmetry (but otherwise arbitrary geometry) oscillating in heave the capture (or absorption) width [math]L_{max}[/math][m] is as follows ^{[5]}:
[math]L_{max} = \large\frac{P_{abs}}{P_{w}} \normalsize = \large\frac{\lambda}{2\pi}\normalsize [/math] or [math]1 = \large\frac{P_{abs}}{P_{w}} \frac{2\pi}{\lambda}\normalsize [/math]
Fig 4: Upper limit of mean wave power absorption for a heaving point absorber. 
where [math]{P_{w}}[/math] is the wave energy flux per unit crest length [W/m]. An optimally damped buoy responds however efficiently to a relatively narrow band of wave periods.
Babarit and Hals propose ^{[6]} to derive that upper limit for the mean annual power in irregular waves at some typical locations where one could be interested in putting some wave energy devices. The mean annual power absorption tends to increase linearly with the wave power resource. Overall, one can say that for a typical site whose resource is between 2030 kW/m, the upper limit of mean wave power absorption is about 1 MW for a heaving WEC with a capture width between 3050 m.
In order to complete these theoretical results and to describe the efficiency of the WEC in practical situations, the capture width ratio [math]\eta[/math] is also usually introduced. It is defined as the ratio between the absorbed power and the available wave power resource per meter of wave front times a relevant dimension [math]B[/math] [m].
[math]\eta = \large\frac{P_{abs}}{P_{w}B}\normalsize [/math]
The choice of the dimension [math]B[/math] will depend on the working principle of the WEC. Most of the time, it should be chosen as the width of the device, but in some cases another dimension is more relevant. Estimations of this ratio [math]\eta[/math] for different systems (see Fig. 3) are given ^{[6]}: 33 % for oscillating water columns (OWC), 13 % for overtopping devices (OD), 929 % for point absorbers (heaving buoys, PA), 2041 % for submerged oscillating bodies (pitching devices, SOB). For energy converted to electricity, one must take into account moreover the energy losses in other components of the system.
Civil engineering
Never forget that the energy conversion is only a secondary function for the coastal structure. The primary function of the coastal structure is still protection. It is necessary to verify whether integration of WEC modifies performance criteria of overtopping and stability and to assess the consequences for the construction cost.
Integration of WEC in coastal structures will always be easier for a new structure than for an existing one. In the latter case, it requires some knowledge on the existing coastal structures. Solutions differ according to sea state but also to type of structures (rubble mound breakwater, caisson breakwaters with typically vertical sides). Some types of WEC are more appropriate with some types of coastal structures.
Fig 5: Several OWC (Oscillating water column) configurations (by Wavegen – Voith Hydro). 
Environmental impact
Wave absorption if it is significant will change hydrodynamics along the structure. If there is mobile bottom in front of the structure, a sand deposit can occur. Ecosystems can also be altered by change of hydrodynamics and but acoustic noise generated by the machines.
Fig 6: Finistere area and locations of the six sites (google map). 
Study case: Finistere area
Finistere area is an interesting study case because it is located in the far west of Brittany peninsula and receives in consequence the largest wave energy flux along the French coasts (see Fig.2). This area with a very ragged coast gathers moreover many commercial ports, fishing ports, yachting ports. The area produces a weak part of its consumption and is located far from electricity power plants. There are therefore needs for renewable energies that are produced locally. This issue is important in particular in islands. The production of electricity by wave energy will have seasonal variations. Wave energy flux is indeed larger in winter than in summer. The consumption has peaks in winter due to heating of buildings but the consumption in summer is also strong due to the arrival of tourists.
Six sites are selected (see figure 7) for a preliminary study of wave energy flux and capacity of integration of wave energy converters. The wave energy flux is expected to be in the range of 1 – 10 kW/m. The length of each breakwater exceeds 200 meters. The wave power along each structure is therefore estimated between 200 kW and 2 MW. Note that there exist much longer coastal structures like for example Cherbourg (France) with a length of 6 kilometres.
Wave power flux along the structure depends on local parameters: bottom depth that fronts the structure toe, the presence of caps, the direction of waves and the orientation of the coastal structure. See figure 8 for the statistics of wave directions measured by a wave buoy located at the Pierres Noires Lighthouse. These measurements show that structures welloriented to West waves should be chosen in priority. Peaks of consumption occur often with low temperatures in winter coming with winds from East NorthEast directions. Structures welloriented to East waves could therefore be also interesting even if the mean production is weak.
Fig 8: Wave measurements at the Pierres Noires Lighthouse. 
Conclusion
Wave energy converters (WEC) in coastal structures can be considered as a land renewable energy. The expected energy can be compared with the energy of land wind farms but not with offshore wind farms whose number and power are much larger. As a land system, the maintenance will be easy. Except the energy production, the advantages of such systems are :
 a “zero emission” port
 industrial tourism
 test of WEC for future offshore installations.
Acknowledgement
This work is in progress in the frame of the national project EMACOP funded by the French Ministry of Ecology, Sustainable Development and Energy.
See also
References
 ↑ ^{1.0} ^{1.1} ^{1.2} Mei C.C. (1989) The applied dynamics of ocean surface waves. Advanced series on ocean engineering. World Scientific Publishing Ltd
 ↑ ^{2.0} ^{2.1} Vicinanza D., Cappietti L., Ferrante V. and Contestabile P. (2011) : Estimation of the wave energy along the Italian offshore, journal of coastal research, special issue 64, pp 613  617.
 ↑ Mattarolo G., Benoit M., Lafon F. (2009), Wave energy resource off the French coasts: the ANEMOC database applied to the energy yield evaluation of Wave Energy, 10th European Wave and Tidal Energy Conference Series (EWTEC’2009), Uppsala (Sweden)
 ↑ Benoit M. and Lafon F. (2004) : A nearshore wave atlas along the coasts of France based on the numerical modeling of wave climate over 25 years, 29th International Conference on Coastal Engineering (ICCE’2004), Lisbonne (Portugal), pp 714726.
 ↑ ^{5.0} ^{5.1} ^{5.2} De O. Falcão A. F. (2010) Wave energy utilization: A review of the technologies. Renewable and Sustainable Energy Reviews, Volume 14, Issue 3, April 2010, pp. 899–918.
 ↑ ^{6.0} ^{6.1} Babarit A. and Hals J. (2011) On the maximum and actual capture width ratio of wave energy converters – 11th European Wave and Tidal Energy Conference Series (EWTEC’2011) – Southampton (UK).
With contributions by: François Bouttes, Bertrand Michard, Emmanuel Cosquer, Alain Clément, Aurélien Babarit, Virginie Baudry, Michel Benoit and Giovanni Mattarolo
