Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms
Abstract
:1. Introduction
2. Numerical Methods
2.1. Spectral-Domain
2.2. Frequency-Domain
2.3. Comparison
3. WEC Types
3.1. Power Performance
3.2. Characteristic Dimensions
4. Wave Conditions
Irregular Waves
5. Comparison Metrics
5.1. Wave Field Disturbance
5.2. Wave Field Difference
5.3. Normalized, Mean-square, and Root-Mean-Square Difference
5.4. RMS Difference Region
5.5. Results Description
5.5.1. Wave Field Plots
- (a)
- SNL-SWAN wave field: This is the wave field produced by SNL-SWAN. Incident waves travel in the positive x direction, as indicated by the arrow. The plot is of the wave field disturbance coefficient, , which is equal to significant wave height at a field point normalized by incident significant wave height. The wave field coefficient represents the change in wave height due to the presence of the WEC: A value of 0 means the wave field is undisturbed, which is represented by green in the plot; a positive value means wave heights are higher than they would otherwise be (hot colors); and a negative value means there is a wave shadow (cool colors). It should be noted that because SNL-SWAN only removes energy from the wave field, the wave field coefficient values are always negative. The WEC or array are centered at . Profiles of the WECs are drawn on the plots.
- (b)
- WAMIT wave field: this is the wave field produced by WAMIT. Similar to the SNL-SWAN wave field, it is a plot of the wave disturbance field coefficient, and the WEC or array is centered at (0,0). However, since WAMIT models wave scattering (reflection) and radiation (wave produced by the motion of the device) as well as wave energy absorption, the plots show increases as well as decreases in the wave height.
- (c)
- Difference wave field: This is a plot of the difference between the WAMIT wave field and the SNL-SWAN wave field in terms of their wave field coefficient, . A value of 0 means there is no difference between the two wave fields (represented by the pink color). A positive value means that the WAMIT wave height is greater than the SNL-SWAN wave height at that point (represented by a light color). A negative value means that the WAMIT wave height is lower than the SNL-SWAN wave height at that point (represented by a dark color). In this plot, there are two blue circles centered around the WEC(s), referred to as the difference region. It is over these circles that quantitative values are given. The radii of these circles are:
- Single WEC: r = 5, 10 d;
- WEC Array: r = 30, 100 d;
where d is the characteristic diameter of the WEC (see Table 1). The radii of 5 and 10 d were used during the initial single WEC analysis, and larger regions are considered for the WEC array analysis.
5.5.2. RMS Difference/Disturbance Plots
5.5.3. Spectra Plots
- Single WEC: (−100,0), (50,0), (500,0), (500,100);
- WEC Array: (−20d,0), (20d,0), (100d,0), (86.6d,50d).
6. Single WEC Cases
6.1. Short Waves
6.2. Unidirectional Waves
6.3. Open Ocean Waves
6.3.1. Pitching Flap
6.3.2. Point Absorber
6.3.3. Hinged Raft Flap
7. WEC Arrays
7.1. Open Ocean Waves
7.2. Array Spacing
8. Conclusions
8.1. Single WEC
8.2. WEC Array
8.3. Guidance
- For arrays, SNL-SWAN shows differences of between 20% and 60% as compared to WAMIT, where the difference is normalized by the impact of the array on the wave field, which are between 1% and 5% at a distance of 30 characteristic diameters from the WECs;
- Directional wave spreading in the incident wave spectrum is important for SNL-SWAN accuracy; unidirectional waves are not modeled well;
- SNL-SWAN does not model wave reflection or scattering, which is a mechanism (in addition to power absorption) for the creation of a wave shadow. This causes larger errors in shorter waves. Preliminary results show the errors could be decreased by including reflections in SNL-SWAN;
- Truncating power absorption, because of RCW values greater than 1, worsens wave field modeling accuracy.
Author Contributions
Funding
Acknowledgments
Conflicts of Interest
References
- Folley, M. Numerical Modelling of Wave Energy Converters: State-of-the-Art Techniques for Single Devices and Arrays, 1st ed.; Elsevier: Amsterdam, The Netherlands, 2016. [Google Scholar]
- Lee, C.-H.; Newman, J.N. WAMIT 2017. WAMIT Software. Available online: http://www.wamit.com/ (accessed on 1 December 2019).
- Babarit, A. NEMOH 2017. NEMOH Software. Available online: https://lheea.ec-nantes.fr/logiciels-et-brevets/nemoh-presentation-192863.kjsp?RH=1489859622433 (accessed on 1 December 2019).
- Cummins, W.E. The Impulse Response Function and Ship Motions. 1962. Available online: https://books.google.com/books?id=GLLANwAACAAJ (accessed on 1 December 2019).
- TELEMAC. TOMAWAC 2017. TOMAWAC Software. Available online: http://www.opentelemac.org/index.php/modules-list/20-tomawac/ (accessed on 1 December 2019).
- TU Delft. SWAN 2019. SWAN (Simulating WAves Nearshore) Code. Available online: http://swanmodel.sourceforge.net/ (accessed on 1 December 2019).
- Millar, D.; Smith, H.; Reeve, D. Modelling analysis of the sensitivity of shoreline change to a wave farm. Ocean. Eng. 2007, 34, 884–901. [Google Scholar] [CrossRef]
- Smith, H.C.; Pearce, C.; Millar, D.L. Further analysis of change in nearshore wave climate due to an offshore wave farm: An enhanced case study for the Wave Hub site. Renew. Energy 2011, 40, 51–64. [Google Scholar] [CrossRef]
- Silverthorne, K.E.; Folley, M. A New Numerical Representation of Wave Energy Converters in a Spectral Wave Model. In Proceedings of the 9th European Wave and Tidal Energy Conference, Southampton, UK, 5–9 September 2011. [Google Scholar]
- Haller, M.C.; Porter, A.; Lenee-Bluhm, P.; Rhinefrank, K.; Hammagren, E.; Özkan Haller, T.; Newborn, D. Laboratory Observation of Waves in the Vicinity of WEC-Arrays. In Proceedings of the 9th European Wave and Tidal Energy Conference, Southampton, UK, 5–9 September 2011. [Google Scholar]
- Porter, A. Laboratory Observations and Numerical Modeling of the Effects of an Array of Wave Energy Converters. Master’s Thesis, Oregon State University, Corvallis, OR, USA, 2012. [Google Scholar]
- McNatt, J.C. Wave Field Patterns Generated by Wave Energy Converters. Master’s Thesis, Oregon State University, Corvallis, OR, USA, 2012. [Google Scholar]
- Rusu, E.; Onea, F. Study on the influence of the distance to shore for a wave energy farm operating in the central part of the Portuguese nearshore. Energy Convers. Manag. 2016, 114, 209–223. [Google Scholar] [CrossRef]
- Greenwood, C.; Christie, D.; Venugopal, V.; Morrison, J.; Vogler, A. Modelling performance of a small array of Wave Energy Converters: Comparison of Spectral and Boussinesq models. Energy 2016, 113, 258–266. [Google Scholar] [CrossRef] [Green Version]
- Venugopal, V.; Nemalidinne, R.; Vögler, A. Numerical modelling of wave energy resources and assessment of wave energy extraction by large scale wave farms. Ocean. Coast. Manag. 2017, 147, 37–48. [Google Scholar] [CrossRef]
- SNL-SWAN, 2019. SNL-SWAN (Sandia National Laboratories Simulating WAves Nearshore) Software. Available online: https://github.com/SNL-WaterPower/SNL-SWAN (accessed on 1 December 2019).
- Chang, G.; Ruehl, K.; Jones, C.A.; Roberts, J.; Chartrand, C. Numerical modeling of the effects of wave energy converter characteristics on nearshore wave conditions. Renew. Energy 2015, 89, 636–648. [Google Scholar] [CrossRef] [Green Version]
- Roberts, J.; Jones, C. Marine Renewable Energy Sediment Stability Evaluation Framework. In Proceedings of the 11th European Wave and Tidal Energy Conference, Nantes, France, 6–11 September 2015. [Google Scholar]
- Ruehl, K.; Porter, A.; Posner, A.; Roberts, J. Development of SNL-SWAN, a Validated Wave Energy Converter Array Modeling Tool. In Proceedings of the 10th European Wave and Tidal Energy Conference, Aalborg, Denmark, 2–5 December 2013. [Google Scholar]
- McNatt, J.C.; Porter, A.; Ruehl, K.; Roberts, J. Verification of the SNL-SWAN Spectral WEC-wave Model with Phase-resolved Linear Wave Fields. In Proceedings of the 12th European Wave and Tidal Energy Conference, Cork, Ireland, 27 August–1 September 2017. [Google Scholar]
- Dallman, A.R.; Neary, V.S. Characterization of U.S. Wave Energy Converter (WEC) Test Sites: A Catalogue of Met-Ocean Data; SAND2014-18206; Sandia National Laboratories: Albuquerque, NM, USA, 2014. [CrossRef]
- Chakrabarti, S.K. Hydrodynamics of Offshore Structures; Google-Books-ID: RNGW9CucxQsC; WIT Press: New Forest, UK, 1987. [Google Scholar]
SNL-SWAN | WAMIT | |
---|---|---|
Governing Principle | Conservation of spectral wave energy | Linear wave theory |
WEC Implementation | Extracts energy from spectrally defined wave field according to frequency-dependent definition of WEC power performance, resulting in reduction of wave height. | Computes WEC motions and wave fields via linear wave theory and superimposes the scattered and radiated wave fields on the incident. The net wave field shows interference patterns, particularly the wave shadow. |
Assumptions |
|
|
Intended Application | Modeling nearshore/far-field impacts of wave transformation over realistic bathymetry in a large (regional) domain | Modeling response of WECs and fluid domain subject to waves in a small (local) domain |
Pros |
|
|
Cons |
|
|
Device | Design | DOFs | PTO | Character | Dim 2 | Dim 3 | Character | Water |
---|---|---|---|---|---|---|---|---|
Name | Archetype | DOF | Width (w) | Diameter (d) | Depth (h) | |||
Pitching | Terminator | pitch | pitch | 24 m | 20 m | 1.5 m | 14 m | 25 m |
Flap (PF) | (beam) | (draft) | (thick) | |||||
Point | Point | surge, | heave | 10 m | 5 m | - | 10 m | 60 m |
Absorber | Absorber | sway, | (diameter) | (draft) | ||||
(PA) | heave | |||||||
Hinged | Attenuator | surge, | relative | 8 m | 90 m | 1.5 m | 16 m | 60 m |
Raft (HR) | sway, | pitch | (beam) | (length) | (draft) | |||
heave, | ||||||||
roll, | ||||||||
pitch, | ||||||||
yaw, | ||||||||
relative | ||||||||
pitch |
Single WEC | WEC Array | |
---|---|---|
Wave Periods | 4, 8, 12 s | 6, 8, 12 s |
Regular Waves | run | not run |
Incident Wave | 0, 40 | 0 |
Directions | ||
Spectra | Bretschneider, | Bretschneider |
JONSWAP | ||
Directional | Unidirectional | |
Spreading | s = 20 | s = 20 |
s = 10 |
The match between the SNL-SWAN and WAMIT wave field is perfect. There is no difference. | |
SNL-SWAN improves the fit of the wave field to WAMIT (compared to not modeling the WEC at all). Lower values are better. One could consider to represent the percentage difference between SNL-SWAN and WAMIT. That is, a value of means that there is a 20% difference between SNL-SWAN and WAMIT. | |
SNL-SWAN makes no net improvement in the wave field as compared to not modeling the WEC. A wave field without WEC effects scores . | |
SNL-SWAN predicts the wave field more poorly than if no WEC effects were present. This would occur when SNL-SWAN predicts a wave shadow where WAMIT does not predict a wave shadow. |
WEC | T | Spec | Dir | Spread | SWAN Diff | 5d | 10d | 5d | 10d |
---|---|---|---|---|---|---|---|---|---|
PF | 4 s | Bret | 0 | 10 | N | 0.78 | 0.79 | 0.07 | 0.03 |
PF | 8 s | Bret | 0 | UD | N | 0.99 | 0.91 | 0.05 | 0.03 |
PF | 8 s | Bret | 0 | 10 | N | 0.48 | 0.37 | 0.04 | 0.02 |
PA | 8 s | Bret | 0 | 10 | N | 0.45 | 0.38 | 0.02 | 0.01 |
HR | 8 s | Bret | 0 | 10 | N | 0.60 | 0.54 | 0.03 | 0.01 |
WEC | Rows | Space () | T | Spec | Spread | 30d | 100d | 30d | 100d |
---|---|---|---|---|---|---|---|---|---|
PF | 2 | 10 | 8 s | Bret | 20 | 0.34 | 0.34 | 0.04 | 0.01 |
PA | 2 | 5 | 8 s | Bret | 20 | 0.34 | 0.33 | 0.01 | 0 |
PA | 2 | 10 | 8 s | Bret | 20 | 0.36 | 0.35 | 0.01 | 0 |
© 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Share and Cite
McNatt, J.C.; Porter, A.; Ruehl, K. Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms. J. Mar. Sci. Eng. 2020, 8, 168. https://doi.org/10.3390/jmse8030168
McNatt JC, Porter A, Ruehl K. Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms. Journal of Marine Science and Engineering. 2020; 8(3):168. https://doi.org/10.3390/jmse8030168
Chicago/Turabian StyleMcNatt, J. Cameron, Aaron Porter, and Kelley Ruehl. 2020. "Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms" Journal of Marine Science and Engineering 8, no. 3: 168. https://doi.org/10.3390/jmse8030168
APA StyleMcNatt, J. C., Porter, A., & Ruehl, K. (2020). Comparison of Numerical Methods for Modeling the Wave Field Effects Generated by Individual Wave Energy Converters and Multiple Converter Wave Farms. Journal of Marine Science and Engineering, 8(3), 168. https://doi.org/10.3390/jmse8030168