In this work, periodic lateral boundaries are developed in a time dependent mild-slope equation model, MILDwave, for the accurate generation of regular waves and irregular long and short crested waves in any direction. A single wave generation line inside the computational domain is combined with periodic lateral boundaries. This generation layout yields a homogeneous and thus accurate wave field in the whole domain in contrast to an L-shaped and an arc-shaped wave generation layout where wave diffraction patterns appear inside the computational domain as a result of the intersection of the two wave generation lines and the interaction with the lateral sponge layers. In addition, the performance of the periodic boundaries was evaluated for two different wave synthesis methods for short crested waves generation, a method proposed by Miles and a method proposed by Sand and Mynett. The results show that the MILDwave model with the addition of periodic boundaries and the Sand and Mynett method is capable of reproducing a homogeneous wave field as well as the target frequency spectrum and the target directional spectrum with a low computational cost. The overall performance of the developed model is validated with experimental results for the case of wave transformation over an elliptic shoal (Vincent and Briggs shoal experiment). The numerical results show very good agreement with the experimental data. The proposed generation layout using periodic lateral boundaries makes the mild-slope wave model, MILDwave, an essential tool to study coastal areas and wave energy converter (WEC) farms under realistic 3D wave conditions, due to its significantly small computational cost and its high numerical stability and robustness.

In this study, a series of modules is integrated into a wave-to-wire (W2W) model that links a Boundary Element Method (BEM) solver to a Wave Energy Converter (WEC) motion solver which are in turn coupled to a wave propagation model. The hydrodynamics of the WECs are resolved in the wave structure interaction solver NEMOH, the Power Take-off (PTO) is simulated in the WEC simulation tool WEC-Sim, and the resulting perturbed wave field is coupled to the mild-slope propagation model MILDwave. The W2W model is run for verified for a realistic wave energy project consisting of a WEC farm composed of 10 5-WEC arrays of Oscillating Surging Wave Energy Converters (OSWECs). The investigated WEC farm is modelled for a real wave climate and a sloping bathymetry based on a proposed OSWEC array project off the coast of Bretagne, France. Each WEC array is arranged in a power-maximizing 2-row configuration that also minimizes the inter-array separation distance d x and d y and the arrays are located in a staggered energy maximizing configuration that also decreases the along-shore WEC farm extent. The WEC farm power output and the near and far-field effects are simulated for irregular waves with various significant wave heights wave peak periods and mean wave incidence directions β based on the modelled site wave climatology. The PTO system of each WEC in each farm is modelled as a closed-circuit hydraulic PTO system optimized for each set of incident wave conditions, mimicking the proposed site technology, namely the WaveRoller® OSWEC developed by AW Energy Ltd. The investigation in this study provides a proof of concept of the proposed W2W model in investigating potential commercial WEC projects.

In this work, internal wave generation techniques are developed in an open source non-hydrostatic wave model (Simulating WAves till SHore, SWASH) for accurate generation of regular and irregular long-crested waves. Two different internal wave generation techniques are examined: a source term addition method where additional surface elevation is added to the calculated surface elevation in a specific location in the domain and a spatially distributed source function where a spatially distributed mass is added in the continuity equation. These internal wave generation techniques in combination with numerical wave absorbing sponge layers are proposed as an alternative to the weakly reflective wave generation boundary to avoid re-reflections in case of dispersive and directional waves. The implemented techniques are validated against analytical solutions and experimental data including water surface elevations, orbital velocities, frequency spectra and wave heights. The numerical results show a very good agreement with the analytical solution and the experimental data indicating that SWASH with the addition of the proposed internal wave generation technique can be used to study coastal areas and wave energy converter (WEC) farms even under highly dispersive and directional waves without any spurious reflection from the wave generator.

The weakly reflective wave generation is a wave generation and absorption method in phase-resolving models, based on the assumption that the waves propagating towards the wave generation boundary are small amplitude shallow water waves with direction perpendicular to the boundary. This assumption makes the method weakly reflective for dispersive and directional waves. The internal wave generation method was proposed by Vasarmidis et al. (2019b) as an alternative, for the non-hydrostatic wave model, SWASH, to avoid reflections. In this study, a comparison is made between the performance of the new internal wave generation method and the weakly reflective wave generation method. It is shown that using the internal wave generation leads to a significantly more accurate prediction of the resulting wave field in case of waves reflected back to the numerical boundary. Additionally, the internal wave generation method is extended to short-crested waves and SWASH is validated for the first time with experimental data for the cases of wave propagation over a shoal and wave diffraction around a wall. The proposed extended internal wave generation method increases the capability of SWASH towards the study of wave propagation of highly dispersive short-crested waves in coastal environments with minimal reflection from the boundaries.

In the recent years, the potential impact of wave energy converter (WEC) arrays on the surrounding wave field has been studied using both phase-averaging and phase-resolving wave propagation models. Obtaining understanding of this impact is important because it may affect other users in the sea or on the coastline. However, in these models a parametrization of the WEC power absorption is often adopted. This may lead to an overestimation or underestimation of the overall WEC array power absorption, and thus to an unrealistic estimation of the potential WEC array impact. WEC array power absorption is a result of energy extraction from the incoming waves, and thus wave height decrease is generally observed downwave at large distances (the so-called “wake” or “far-field” effects). Moreover, the power absorption depends on the mutual interactions between the WECs of an array (the so-called “near field” effects). To deal with the limitations posed by wave propagation models, coupled models of recent years, which are nesting wave-structure interaction solvers into wave propagation models, have been used. Wave-structure interaction solvers can generally provide detailed hydrodynamic information around the WECs and a more realistic representation of wave power absorption. Coupled models have shown a lower WEC array impact in terms of wake effects compared to wave propagation models. However, all studies to date in which coupled models are employed have been performed using idealized long-crested waves. Ocean waves propagate with a certain directional spreading that affects the redistribution of wave energy in the lee of WEC arrays, and thus gaining insight wake effect for irregular short-crested sea states is crucial. In our research, a new methodology is introduced for the assessment of WEC array wake effects for realistic sea states. A coupled model is developed between the wave-structure interaction solver NEMOH and the wave propagation model MILDwave. A parametric study is performed showing a comparison between WEC array wake effects for regular, long-crested irregular, and short-crested irregular waves. For this investigation, a nine heaving-point absorber array is used for which the wave height reduction is found to be up to 8% lower at 1.0 km downwave the WEC array when changing from long-crested to short-crested irregular waves. Also, an oscillating wave surge WEC array is simulated and the overestimation of the wake effects in this case is up to 5%. These differences in wake effects between different wave types indicates the need to consider short-crested irregular waves to avoid overestimating the WEC array potential impacts. The MILDwave-NEMOH coupled model has proven to be a reliable numerical tool, with an efficient computational effort for simulating the wake effects of two different WEC arrays under the action of a range of different sea states.