SlideShare a Scribd company logo
03 June 2016
Study for the Hull Shape of a
Wave Energy Converter-Point Absorber
Filippos Kalofotias
Design Optimization & Modeling Improvement
Contents
 Background
 Problem Definition
 Research Objectives and Questions
 Model Description
 Design Optimization
 Modeling Improvement
 Model Validation
 Conclusions & Recommendations
 Questions?
03 June 2016Filippos Kalofotias 2
Background
03 June 2016Filippos Kalofotias 3
‘’Deciphering the title’’
Wave Energy = Energy transfer from wind waves
Converters = Devices transforming wave
energy to electrical
Point Absorber
Background
03 June 2016Filippos Kalofotias 4
Point Absorber
Problem Definition
03 June 2016Filippos Kalofotias 5
P. Wellens (2004) / S. Kao (2014)
Studied :
 Maximization of power extraction under irregular
waves
 Optimum Spring and Damper Configuration
 During sea state (Wellens)
 During smaller intervals of the sea state (Kao)
Did not study :
 Optimum dimensions of the buoy for the
studied wave climate
 Other shapes for the buoy than Cylinder
 Viscous Effects
 Wave force estimation with large motions
Control
Strategy
Design
Optimization
+
Physical
Modeling
Research Objectives & Questions
03 June 2016Filippos Kalofotias 6
Research Objectives
 Derive an efficient design for the
buoy (hull) of the Point Absorber
Shape
+
Dimensions
 Improve physical modeling of
the Point Absorber’s response to
waves
Viscous
+
Wave Force
Research Objectives & Questions
03 June 2016Filippos Kalofotias 7
Research Questions
What is the most efficient design?
 Influence of shape to efficiency?
 Influence of dimensions to efficiency?
 Efficiency comparison of different designs?
How to improve Point Absorber’s modeling?
 Viscous effects estimation?
 Wave force including large buoy’s motion?
 Influence of additions to the previous model?
Model Description
03 June 2016Filippos Kalofotias 8
Linear Mass-Spring-Damper System
Newton’s 2nd Law
Force Analysis
Restricted to up and down motion
Equation
Solution?
Coefficients?
Excitation Force?
Nonlinear Viscous
Effects?
Model Description
03 June 2016Filippos Kalofotias 9
Solution
Frequency Domain
Time Domain
 Irregular waves – JONSWAP
 Superposition of component
regular waves - RAO
 Linearity is a must
 Fast and valuable for statistics
 Time step solution of buoy’s
position and velocity
 Phase shifts needed
 Instantaneous power and
average over time
 Nonlinear effects such viscous
forces can be included
Model Description
03 June 2016Filippos Kalofotias 10
Coefficients
3D – Diffraction Theory / NEMOH
Spring, ksp and damper, β
 Radiation Damping, b(ω)
 Added Mass, a(ω)
 Predetermined and adjustable
 Resonance at peak frequency
 Optimum damper coefficient
Restoring, c
 Determined by waterline area
 Radius dependent
Model Description
03 June 2016Filippos Kalofotias 11
Excitation Force
3D – Diffraction Theory / NEMOH
 Force amplitude per frequency
 Per meter wave amplitude
 Doubling the wave amplitude
doubles the wave force
 Equilibrium position estimation
 Superposition principle
 For time dependent, Fexc (t)
phase shift is needed
NEMOH neglects viscous effects
Nonlinear Viscous
Effects?
Model Description
03 June 2016Filippos Kalofotias 12
Nonlinear Viscous Effects
Computational Fluid Dynamics / ComFLOW3
 Fully viscous numerical solution of
Navier-Stokes equations
 Numerical tank
 Boundary conditions
 Turbulence modeling
 Grid size
 Free surface
 Computationally expensive
 ComFLOW3 provides viscous
information for the Time Domain
model
Design Optimization
03 June 2016Filippos Kalofotias 13
Wave Climate
Scatter Diagram
Evaluation Criterion
 Hs ≤ 4.5m (95% of annual)
 Hs > 4.5m non operational
Efficiency
 Proportion of extracted power to
available
 Available depends on radius
Design Optimization
03 June 2016Filippos Kalofotias 14
Geometry
3 Shapes
Dimensions
 Max R=10m
 Max TD=15m
Design Optimization
Dimensioning
Frequency Domain Results
03 June 2016Filippos Kalofotias 15
Design Optimization
03 June 2016Filippos Kalofotias 16
Shape Evaluation
 3 Final designs, Cyl8, Bul6, Con6
 Similar efficiencies
 Viscous assessment is needed
How do we include viscous effects in
the Time Domain model?
Drag Force (Morison type)
Oscillating body in calm water
Oscillating body in waves / stagnation pressure
How do we find Cd?
Design Optimization
03 June 2016Filippos Kalofotias 17
Forced Oscillation Tests in ComFLOW3
Design Optimization
03 June 2016Filippos Kalofotias 18
Efficiency Comparison
 Cylinder has the largest drag
coefficients
 Bullet has smaller drag
coefficient for fast oscillations
than the Cone
 Drag force becomes important
at high velocities
 Drag coefficients vary according
to the flow conditions
 So, not steady during simulation
Drag coefficients, Cd results
Design Optimization
03 June 2016Filippos Kalofotias 19
Efficiency Comparison
Time Domain Results
with steady drag coefficients
Frequency Domain Results
no viscous forces
Spring coefficients per sea
state for Bullet and Cone
Design Optimization
03 June 2016Filippos Kalofotias 20
Final Design Selection
 Resonates further in the
scatter diagram
 Highest efficiency at the
most energetic sea states
 Less drag forces for fast
oscillations
 PTO damper optimization
can increase efficiency
Modeling Improvement
03 June 2016Filippos Kalofotias 21
Final Model 1
Parameterization of Cd
 Detailed inclusion of viscous effects
 Adjustment of Cd at every time step
Reynolds number
Keulegan-Carpenter number
Forced Oscillation Tests
Modeling Improvement
03 June 2016Filippos Kalofotias 22
Final Model 1
Assumption : Re and KC refer to the maximum velocity
Vm. In Time Domain model the maximum velocity is not
known before hand. The actual is used.
Results
What about PTO damping
coefficient optimization?
PTO damping coefficient selective optimization
Modeling Improvement
03 June 2016Filippos Kalofotias 23
Final Model 2
 Inclusion of the buoy’s motion in
excitation force estimation
 Calculation of Froude-Krylov force
at every time step
 NEMOH runs at 11 different
positions for the Diffraction force
 Estimation of Diffraction force at
every time step with interpolation
Results
Excitation force comparison
Model Validation
03 June 2016Filippos Kalofotias 24
Frequency vs Time Domain
Δω = 0.01 rad/s Δω = 0.001 rad/s
Model Validation
03 June 2016Filippos Kalofotias 25
Final Model 1 vs ComFLOW3
Forced Oscillation Test
Forced Oscillation Test
With waves
Conclusions & Recommendations
03 June 2016Filippos Kalofotias 26
Conclusions
What is the most efficient design?
 Within the 21 evaluated designs, Bul6 was qualified as the most efficient
 No guarantee that Bul6 is the optimum design
Influence of shape and dimensions?
 Coupled optimization problem
 Different shape leads to different optimum dimensions
 Resonance properties
 Excitation force
Efficiency comparison of different designs?
 Optimum dimensions per shape
 Viscous effects assessment
 Inclusion of viscous damping in PTO damper optimum configuration
Conclusions & Recommendations
03 June 2016Filippos Kalofotias 27
Conclusions
How to improve Point Absorber’s modeling?
 Final Model 1 increased the physical accuracy of the model.
 Final Model 2 proved to be computationally expensive without adding
significant physical information to the model
Viscous effects estimation?
 Stagnation pressure method for drag force
 Forced Oscillation Tests
 Drag coefficient parameterization
 Time step adjustment
Wave force including large buoy’s motion?
 It appear that for large buoys the assumption of estimating the excitation
force at equilibrium position is valid
 Linear wave theory and linear diffraction theory cannot add anymore
physical accuracy in the model
Conclusions & Recommendations
03 June 2016Filippos Kalofotias 28
Conclusions
Influence of additions to the previous model?
 More than 10% decrease in power extraction because of drag force
 Optimum PTO damper configuration shift as a result of the extra
damping
 Important for control strategies
Recommendations
 Nonlinear excitation force assessment
 Nonlinear restoring and radiation forces
 PTO modeling
 Point Absorber farms
03 June 2016Filippos Kalofotias 29
Questions?

More Related Content

PDF
Thesis Summary
PDF
Thesis Report
PDF
Question and answers webinar hydrodynamic modeling on the northwest european ...
PPTX
URC Poster
PDF
Webinar presentation ocean modelling and early-warning system for the gulf ...
PDF
Webinar ocean-modelling-and-gulf-of-thailand-questions-and-answers
PDF
Current research on simulations of flaoting offshore wind turbines
PDF
Simulation requirements and relevant load conditions in the design of floatin...
Thesis Summary
Thesis Report
Question and answers webinar hydrodynamic modeling on the northwest european ...
URC Poster
Webinar presentation ocean modelling and early-warning system for the gulf ...
Webinar ocean-modelling-and-gulf-of-thailand-questions-and-answers
Current research on simulations of flaoting offshore wind turbines
Simulation requirements and relevant load conditions in the design of floatin...

What's hot (20)

PPT
Performance studies on a direct drive turbine for wave power generation in a ...
PDF
Chemeca2015
PDF
Evaluation of the extreme and fatigue load measurements at alpha ventus
PDF
Towards the identification of the primary particle nature by the radiodetecti...
PDF
Webinar deltares hydrodynamic modeling on the northwest european shelf and no...
PDF
Weak and strong oblique shock waves
PDF
"Field characterization of location-specific dynamic amplification factors to...
PDF
Determination of Kite forces for ship propulsion
PDF
"Reduction of uncertainties associated to the dynamic response of a ship unlo...
PDF
Towards the identification of the primary particle nature by the radiodetecti...
PDF
Comparative Calibration Method Between two Different Wavelengths With Aureole...
PPTX
Voller_s4321724_presentation
PDF
Hydrostatics 1 n 2
PDF
Combustion and Mixing Analysis of a Scramjet Combustor Using CFD
PDF
4.11 radar cross section (rcs)
PDF
Methodology for Estimating Wave Power Potential in places with scarce instrum...
PPSX
Balancing of an air-bearing-based Acs Test Bed
PDF
Fluid discharge
PDF
WaReS Validation Report
PPTX
CALPUFF- Air Quality modelling
Performance studies on a direct drive turbine for wave power generation in a ...
Chemeca2015
Evaluation of the extreme and fatigue load measurements at alpha ventus
Towards the identification of the primary particle nature by the radiodetecti...
Webinar deltares hydrodynamic modeling on the northwest european shelf and no...
Weak and strong oblique shock waves
"Field characterization of location-specific dynamic amplification factors to...
Determination of Kite forces for ship propulsion
"Reduction of uncertainties associated to the dynamic response of a ship unlo...
Towards the identification of the primary particle nature by the radiodetecti...
Comparative Calibration Method Between two Different Wavelengths With Aureole...
Voller_s4321724_presentation
Hydrostatics 1 n 2
Combustion and Mixing Analysis of a Scramjet Combustor Using CFD
4.11 radar cross section (rcs)
Methodology for Estimating Wave Power Potential in places with scarce instrum...
Balancing of an air-bearing-based Acs Test Bed
Fluid discharge
WaReS Validation Report
CALPUFF- Air Quality modelling
Ad

Similar to Final presentation (20)

PPTX
Modelon JSME 2016 - Model Based Design for Fuel Cell Systems
PDF
20 Years and 20Mt, Statoil Storage Experience, Andrew Cavanagh - Geophysical ...
PPTX
Thesis presentation
PPTX
McGill Ozone Contactor Design Project
PDF
C012512230
PDF
Cfd study of a Pelton turbine runner by M.Minozzo R.Bergamin M.Merelli M.Galb...
 
PPTX
CRP presentation final
PDF
TSEP_2019_103_Original_V0.pdf
PDF
HYWEC (Hydrodynamics of Wave Energy Converters) Workshop at BCAM April 3-7, 2017
PPTX
Presentation CIE619
PDF
johan_malmberg_thesis_pres_final
PDF
Morphological model of the river rhine branches from the concept to the opera...
PDF
Project_ Report
PDF
CV_Vidal_s222447_CFD_Poster
PPTX
Determination of shock losses and pressure losses in ug mine openings (1)
PPTX
Determination of shock losses and pressure losses in ug mine openings
PPT
TU1.L10 - Arctic Sea Ice dynamics for Global Climate Models: Results from the...
PDF
3. Enhance DCM
PPTX
PRIME2 Model Evaluation
PPT
Development of Economical Analysis and Technical Solutions for Efficient Dist...
Modelon JSME 2016 - Model Based Design for Fuel Cell Systems
20 Years and 20Mt, Statoil Storage Experience, Andrew Cavanagh - Geophysical ...
Thesis presentation
McGill Ozone Contactor Design Project
C012512230
Cfd study of a Pelton turbine runner by M.Minozzo R.Bergamin M.Merelli M.Galb...
 
CRP presentation final
TSEP_2019_103_Original_V0.pdf
HYWEC (Hydrodynamics of Wave Energy Converters) Workshop at BCAM April 3-7, 2017
Presentation CIE619
johan_malmberg_thesis_pres_final
Morphological model of the river rhine branches from the concept to the opera...
Project_ Report
CV_Vidal_s222447_CFD_Poster
Determination of shock losses and pressure losses in ug mine openings (1)
Determination of shock losses and pressure losses in ug mine openings
TU1.L10 - Arctic Sea Ice dynamics for Global Climate Models: Results from the...
3. Enhance DCM
PRIME2 Model Evaluation
Development of Economical Analysis and Technical Solutions for Efficient Dist...
Ad

Final presentation

  • 1. 03 June 2016 Study for the Hull Shape of a Wave Energy Converter-Point Absorber Filippos Kalofotias Design Optimization & Modeling Improvement
  • 2. Contents  Background  Problem Definition  Research Objectives and Questions  Model Description  Design Optimization  Modeling Improvement  Model Validation  Conclusions & Recommendations  Questions? 03 June 2016Filippos Kalofotias 2
  • 3. Background 03 June 2016Filippos Kalofotias 3 ‘’Deciphering the title’’ Wave Energy = Energy transfer from wind waves Converters = Devices transforming wave energy to electrical Point Absorber
  • 4. Background 03 June 2016Filippos Kalofotias 4 Point Absorber
  • 5. Problem Definition 03 June 2016Filippos Kalofotias 5 P. Wellens (2004) / S. Kao (2014) Studied :  Maximization of power extraction under irregular waves  Optimum Spring and Damper Configuration  During sea state (Wellens)  During smaller intervals of the sea state (Kao) Did not study :  Optimum dimensions of the buoy for the studied wave climate  Other shapes for the buoy than Cylinder  Viscous Effects  Wave force estimation with large motions Control Strategy Design Optimization + Physical Modeling
  • 6. Research Objectives & Questions 03 June 2016Filippos Kalofotias 6 Research Objectives  Derive an efficient design for the buoy (hull) of the Point Absorber Shape + Dimensions  Improve physical modeling of the Point Absorber’s response to waves Viscous + Wave Force
  • 7. Research Objectives & Questions 03 June 2016Filippos Kalofotias 7 Research Questions What is the most efficient design?  Influence of shape to efficiency?  Influence of dimensions to efficiency?  Efficiency comparison of different designs? How to improve Point Absorber’s modeling?  Viscous effects estimation?  Wave force including large buoy’s motion?  Influence of additions to the previous model?
  • 8. Model Description 03 June 2016Filippos Kalofotias 8 Linear Mass-Spring-Damper System Newton’s 2nd Law Force Analysis Restricted to up and down motion Equation Solution? Coefficients? Excitation Force? Nonlinear Viscous Effects?
  • 9. Model Description 03 June 2016Filippos Kalofotias 9 Solution Frequency Domain Time Domain  Irregular waves – JONSWAP  Superposition of component regular waves - RAO  Linearity is a must  Fast and valuable for statistics  Time step solution of buoy’s position and velocity  Phase shifts needed  Instantaneous power and average over time  Nonlinear effects such viscous forces can be included
  • 10. Model Description 03 June 2016Filippos Kalofotias 10 Coefficients 3D – Diffraction Theory / NEMOH Spring, ksp and damper, β  Radiation Damping, b(ω)  Added Mass, a(ω)  Predetermined and adjustable  Resonance at peak frequency  Optimum damper coefficient Restoring, c  Determined by waterline area  Radius dependent
  • 11. Model Description 03 June 2016Filippos Kalofotias 11 Excitation Force 3D – Diffraction Theory / NEMOH  Force amplitude per frequency  Per meter wave amplitude  Doubling the wave amplitude doubles the wave force  Equilibrium position estimation  Superposition principle  For time dependent, Fexc (t) phase shift is needed NEMOH neglects viscous effects Nonlinear Viscous Effects?
  • 12. Model Description 03 June 2016Filippos Kalofotias 12 Nonlinear Viscous Effects Computational Fluid Dynamics / ComFLOW3  Fully viscous numerical solution of Navier-Stokes equations  Numerical tank  Boundary conditions  Turbulence modeling  Grid size  Free surface  Computationally expensive  ComFLOW3 provides viscous information for the Time Domain model
  • 13. Design Optimization 03 June 2016Filippos Kalofotias 13 Wave Climate Scatter Diagram Evaluation Criterion  Hs ≤ 4.5m (95% of annual)  Hs > 4.5m non operational Efficiency  Proportion of extracted power to available  Available depends on radius
  • 14. Design Optimization 03 June 2016Filippos Kalofotias 14 Geometry 3 Shapes Dimensions  Max R=10m  Max TD=15m
  • 15. Design Optimization Dimensioning Frequency Domain Results 03 June 2016Filippos Kalofotias 15
  • 16. Design Optimization 03 June 2016Filippos Kalofotias 16 Shape Evaluation  3 Final designs, Cyl8, Bul6, Con6  Similar efficiencies  Viscous assessment is needed How do we include viscous effects in the Time Domain model? Drag Force (Morison type) Oscillating body in calm water Oscillating body in waves / stagnation pressure How do we find Cd?
  • 17. Design Optimization 03 June 2016Filippos Kalofotias 17 Forced Oscillation Tests in ComFLOW3
  • 18. Design Optimization 03 June 2016Filippos Kalofotias 18 Efficiency Comparison  Cylinder has the largest drag coefficients  Bullet has smaller drag coefficient for fast oscillations than the Cone  Drag force becomes important at high velocities  Drag coefficients vary according to the flow conditions  So, not steady during simulation Drag coefficients, Cd results
  • 19. Design Optimization 03 June 2016Filippos Kalofotias 19 Efficiency Comparison Time Domain Results with steady drag coefficients Frequency Domain Results no viscous forces Spring coefficients per sea state for Bullet and Cone
  • 20. Design Optimization 03 June 2016Filippos Kalofotias 20 Final Design Selection  Resonates further in the scatter diagram  Highest efficiency at the most energetic sea states  Less drag forces for fast oscillations  PTO damper optimization can increase efficiency
  • 21. Modeling Improvement 03 June 2016Filippos Kalofotias 21 Final Model 1 Parameterization of Cd  Detailed inclusion of viscous effects  Adjustment of Cd at every time step Reynolds number Keulegan-Carpenter number Forced Oscillation Tests
  • 22. Modeling Improvement 03 June 2016Filippos Kalofotias 22 Final Model 1 Assumption : Re and KC refer to the maximum velocity Vm. In Time Domain model the maximum velocity is not known before hand. The actual is used. Results What about PTO damping coefficient optimization? PTO damping coefficient selective optimization
  • 23. Modeling Improvement 03 June 2016Filippos Kalofotias 23 Final Model 2  Inclusion of the buoy’s motion in excitation force estimation  Calculation of Froude-Krylov force at every time step  NEMOH runs at 11 different positions for the Diffraction force  Estimation of Diffraction force at every time step with interpolation Results Excitation force comparison
  • 24. Model Validation 03 June 2016Filippos Kalofotias 24 Frequency vs Time Domain Δω = 0.01 rad/s Δω = 0.001 rad/s
  • 25. Model Validation 03 June 2016Filippos Kalofotias 25 Final Model 1 vs ComFLOW3 Forced Oscillation Test Forced Oscillation Test With waves
  • 26. Conclusions & Recommendations 03 June 2016Filippos Kalofotias 26 Conclusions What is the most efficient design?  Within the 21 evaluated designs, Bul6 was qualified as the most efficient  No guarantee that Bul6 is the optimum design Influence of shape and dimensions?  Coupled optimization problem  Different shape leads to different optimum dimensions  Resonance properties  Excitation force Efficiency comparison of different designs?  Optimum dimensions per shape  Viscous effects assessment  Inclusion of viscous damping in PTO damper optimum configuration
  • 27. Conclusions & Recommendations 03 June 2016Filippos Kalofotias 27 Conclusions How to improve Point Absorber’s modeling?  Final Model 1 increased the physical accuracy of the model.  Final Model 2 proved to be computationally expensive without adding significant physical information to the model Viscous effects estimation?  Stagnation pressure method for drag force  Forced Oscillation Tests  Drag coefficient parameterization  Time step adjustment Wave force including large buoy’s motion?  It appear that for large buoys the assumption of estimating the excitation force at equilibrium position is valid  Linear wave theory and linear diffraction theory cannot add anymore physical accuracy in the model
  • 28. Conclusions & Recommendations 03 June 2016Filippos Kalofotias 28 Conclusions Influence of additions to the previous model?  More than 10% decrease in power extraction because of drag force  Optimum PTO damper configuration shift as a result of the extra damping  Important for control strategies Recommendations  Nonlinear excitation force assessment  Nonlinear restoring and radiation forces  PTO modeling  Point Absorber farms
  • 29. 03 June 2016Filippos Kalofotias 29 Questions?