Wave energy Extracting power from ocean waves Some basic principles Finn Gunnar Nielsen, Geophysical Institute, University of Bergen Finn.Nielsen@uib.no 14..17 1
Issues to be discussed A reminder on the nature of (deep water) gravity waves The Global picture Examples on wave energy converters Extracting wave energy by an oscillating system To absorb wave energy is about generating waves Theoretical efficiencies Two simple numerical examples 14..17
Classification of ocean waves 14..17 3
Linear, deep water waves (Airy waves) Dispersion relation: c As g tanh kd k k g d, tanh( kd) 1, c, kg k Kinetic = potential energy (per unit surface area): Energy flux: 1 H Ep Ev g 4 1 1 H 1 H g P cg Ep Ev g g k 4 k 14..17 4
Sf(f) (m sec) 18 14 1 1 8 6 4 Energy flux in a wave spectrum Peak frequency 16= 1/ (7.47 sec) Jonswap spectrum Mean energy frequency =1/ (6.78 sec) Zero upcrossing frequency = 1/ (5.88 sec] Hs = sqrt(m) = 3.5m.1..3.4.5 Frequency [1/sec] Spectral moments: n mn S d Significant wave height Average period Zero up-crossing period Energy mean period Hs T T T 1 1 4 m m m 1 m m m m 1 Energy flux: 1 1 1 P g m g T m g T H s 4 64 1 1 1 14..17 5
RESOURCES Waves (From World Energy Council) Average energy densities: 1 kw/m wave front 5 kw/m, 3% efficiency: 13 MWh/ym 1TWh/y: 76km 14..17 6
The potential exceeds the demands. Source: IPCC SRREN, 11 14..17 7
Main principles From: Babarit, Introduction to Ocean Wave Energy Conversion, 9 14..17 8
Illustration of main principles Source: Bedard(6)»overview: EPRI Ocean Energy Program» Presented to Duke Univerisity Global Change Centre. 14..17 9
Key principles for extracting wave energy Remember: Waves are not only quasi-static change in surface elevation Energy absorption requires a force working together with a velocity Falnes and Budal (1978): In order for an oscillating system to be a good wave absorber it should be a good wave generator. 14..17 1
The Tapchan project at Toftestallen, Øygarden (finished 1985) Principle. Test site 14..17 11
How is energy extracted? Key principles Linear considerations Linear oscillator interaction with waves Heaving buoy as example 14..17 1
Linear oscillator F(t) M X Dynamic equilibrium: Mx Bx Kx F() t x x x F( t) Harmonic oscillation, stationary solution. i t F t F cos t F Re[ e ] A A K M B = M K B Dynamic equation, frequency domain: M i B K x F 14..17 13
Linear oscillator - solution cos x xa t xare e F x M K ib i t 6 5 - -4 abs(x/x ) 4 3 (Deg) -6-8 -1-1 1-14 -16 1 3 4 5 / -18 1 3 4 5 / At resonance: Response controlled by damping! 14..17 14
Estimating the energy absorption by a heaving buoy Assumptions: Heave motion only Linearized analysis No viscous effects 14..17 15
The heaving point absorber, a linearized approach (1:3) Equation of motion (1DOF): mx Fex FR FP In frequency domain: m A i B B B k K x F 33 r l wl ex Wave radiation force: F A x B x k x Force due to power off-take: Heave motion: x R 33 r wl F Bx Kx P F Aex cos( t) m A k K B B B 33 wl r l 14..17 16
The heaving point absorber, a linearized approach (:3) Instantaneous power: P P t F x Kx Bx x P Integrated over one wave period: P P F B 1 Aex 1 m A33 k wl K B Br Bx A Maximum power at resonance 14..17 17
The heaving point absorber, a linearized approach (3:3) Optimum damping in power off-take: dpp / db Obtained for B B r Maximum mean power: P p _ max FAex 8B r 14..17 18
Theoretical limits for power extraction 3D axisymmetric body (deep water) Invoke Haskind relation (relates wave excitation force to wave radiation damping.) 3 Fex3 3D symmetric, heaving buoy: Br 33 3 g H Falnes and Budal, 1978: In order for an oscillating system to be a good wave absorber it should be a good wave generator. I.e. max power: P 3 H g 3 F Aex g p _ max 3 3 8Br 4 18 H T 3 Capture width : L cap 3 g 3 HT max power extraction 3 18 g 1 g TH 3 power in wave per meter 14..17 19
Heaving point absorber Theoretical versus Budal limit. - Or the implication of limited physical size. 1 C = = 7.9 kws/m 4 4 g 3 g C = =.44 kw/(m s 3 ) 3 18 Capture width (P A ): L=λ/π (Heave only) L = 3λ/π (Heave & surge) Resonance Semisubmerged sphere From Falnes (7) 14..17
E4 /TU 31.5.15 14..17 1
The idea of latching (1:) A. Babarit, G. Duclos, A.H. Clément: Comparison of latching control strategies for a heaving wave energy device in random sea 14..17
The idea of latching (:) Tout = Twa Max. W/o control Tout = 3Twa Source: Clement and Babarit (1) 14..17 1DOF system. Power absorbed in harmonic and first sub-harmonic latching modes, compared with uncontrolled, and ideal modes. Solid line, without control; squares, latching, Tout=Twa; circles, latching, Tout=3 Tin; dashed dotted line, max power 3
Main principles From: Babarit, Introduction to Ocean Wave Energy Conversion, 9 14..17 4
Interaction between waves and body - D Energy flux in incident wave: 1 g J cge A 4 Net absorbed wave power per unit length: g P A A A 4 D R T A T : amplitude of transmitted wave A R : amplitude of reflected wave 14..17 5
Interaction between waves and body D Incident wave plus heaving and surging body. 3 Radiated waves from a heaving D source 3 Radiated waves from a heaving plus surging source 1 1 3 Incident waves moving in positive x-direction wave elevation -1 - -1 - Wave due to heave -3-4 Heave + surge 1 wave elevation -3 - -15-1 -5 5 1 15 x 3 Radiated waves from a surging D source -5 - -15-1 -5 5 1 15 3 Incident plus radiated waves from a heaving and surging source -1-1 1-3 - -15-1 -5 5 1 15 x Incident wave wave elevation -1 - -3-4 Wave due to surge wave elevation -1 - -3-4 Incident + heave + surge -5 - -15-1 -5 5 1 15 x -5 - -15-1 -5 5 1 15 x Note importance of phasing 14..17 6
Example heaving buoy, motions INPUT Radius 1 m Draft 5 m Wave amplitude 1 m Mass 188.599 kg Water line stiffness, 31.5787 kn/m Power offtake stiffness 15.7894 kn/m Natural period 3.539 sec 8 Radiation damping at T 1947.69 N/(m/s) 7 Theoretical maximum energy absorbtion, at resonance (kw) 43.1835 Gam = B_ powerofftake / B_ radiation (T ) B_ powerofftake independent of frequency (not optimum) 6 5 Amplitude of motion gam =.1 gam =.5 gam = 1 gam = gam = 5 gam = 1 abs(x 3 )/A 4 3 x 3 1 x 1 1 3 4 5 6 7 8 9 1 Wave period (sec) 14..17 7
Heaving buoy - Power production Average power (kw) 45 4 35 3 5 15 1 Power production, wave amplitude 1 m gam =.1 gam =.5 gam = 1 gam = gam = 5 gam = 1 rage power (kw) 16 14 1 1 8 6 Budal limit with A= 1m Theoretical maximum versus "Budal limit" Theoretical max "Budal limit" 5 4 1 3 4 5 6 7 8 9 1 Wave period (sec) 4 6 8 1 1 14 Wave period (sec) 14..17 8
Power production off resonance - passive system INPUT Radius.5 m Draft 5 m Wave amplitude 1 m Mass 853.3117 kg Water line stiffness, 197.367 kn/m Power offtake stiffness 98.6835 kn/m Natural period 3.8998 sec Radiation damping at T 4977.4514 N/(m/s) Theoretical maximum energy absorbtion, at resonance (kw) 57.789 Gam = B_powerofftake / B_radiation (T) High damping in power off-take important to extract energy at periods above resonance Average power (kw) 7 6 5 4 3 Power production, wave amplitude 1 m gam =.1 gam =.5 gam = 1 gam = gam = 5 gam = 1 1 4 6 8 1 1 14 Wave period (sec) 14..17 9
Example: Gabriell Figure. Illustration of the buoy with a horizontal plate attached underneath. 14..17 3
Toftestallen In the 198 ies..and in 14 14..17 31
Wave power converters - Examples on installations in full / reduced scale (:) Bostrøm et al. (Sweden, 6-) Heaving buoy. Linear generator Fred Olsen, Buldra. (Norway 4-) Array of heaving buoys. Semisubmersible 14..17 3
Oscillating water column (OWC) LIMPET LIMPET 14..17 33
Relative motion device attenuator Pelamis (75 kw device) 14..17 34
Overtopping Wave Dragon (prototype kw, 4-7MW demonstrator) 14..17 35
New ideas - floating hose The ANACONDA concept. Ill. EPSRC Water inside the tubes are propagating in longitudinal direction. 14..17 36
Summary - Vaste amount of wave energy available. Technical availability uncertain. - No convergence on technical solutions - For an oscillating system to extract energy, it has to generate waves - Advanced control needed to enhance power offtake. - Survivability has shown up to be a critical issue. 14..17 37