An adaptive spectral/hp discontinuous Galerkin method for the two-dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non-conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p-1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h-type refinement, the parent element is subdivided into four similar sibling elements. The time-stepping is performed using a third-order Runge-Kutta scheme. The performance of the hp-adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p-adaptivity is more efficient than h-adaptivity with respect to degrees of freedom and computational time.
This paper presents the methods developed and key findings of the IWEC project performed by Ocean Harvesting Technologies AB (OHT). It aimed to reduce the levelized cost of energy (LCoE) of OHT’s wave energy converter InfinityWEC, by analysing how different key parameters impact cost and annual output using a model of a 100-MW array installation. Component-level cost functions were developed and mapped to key parameters and constraints of the system. A large number of system configurations were then evaluated with a numerically efficient 3 degree-of-freedom (DoF) nonlinear radiationdiffraction model in WEC-Sim along with OHT’s sea statetuned polynomial reactive control (PRC). The most promising configurations were identified and investigated in more detail. The configuration with the best LCoE were finally identified and analysed further, including estimation of the effect of changing the PRC to model predictive control, which resulted in 17-34% higher annual output and 12-23% lower LCoE. The final LCoE was found to be 93-162 EUR MWh at 100 MW installed capacity. An important finding from the study is that using simplified metrics such as CAPEX/ton was found to be irrelevant. Numerical wave tank testing, high-fidelity computational fluid dynamics (CFD), were used to tune the viscous drag of the 3 DoF WEC-Sim model. Applying verification and validation (V&V) techniques the CFD simulations showed a relatively large numerical uncertainty, but the average power and the motion responses were found to be sufficiently accurate.
The paper discusses the use of CFD simulations to analyse the parametric excitation of moored, full scale wave energy converters in six degrees of freedom. We present results of VOF-RANS and VOF-Euler simulations in OpenFOAM!R for two body shapes: (i) a truncated cylinder; and (ii) a cylinder with a smooth hemispherical bottom. Flow characteristics show large differences in smoothness of flow between the hull shapes, where the smoother shape results in a larger heave response. However the increased amplitude makes it unstable and parametric pitch excitation occurs with amplitudes up to 30". The responses in surge, heave and pitch (including the transition to parametric motion) are found to be insensitive to the viscous effects. This is notable as the converters are working in resonance. The effect of viscous damping was visible in the roll motion, where the RANS simulations showed a smaller roll. However, the roll motion was found to be triggered not by wave-body interaction with the incident wave, but by reflections from the side walls. This highlights the importance of controlling the reflections in numerical wave tanks for simulations with WEC motion in six degrees of freedom.