High-fidelity simulations using computational fluid dynamics (CFD) for wave-body interaction are becoming increasingly common and important for wave energy converter (WEC) design. The open source finite volume toolbox OpenFOAM® is one of the most frequently used platforms for wave energy. There are currently two ways to account for moving bodies in OpenFOAM: (i) mesh morph-ing, where the mesh deforms around the body; and (ii) an overlooked mesh method where a separate body mesh moves on top of a background mesh. Mesh morphing is computationally efficient but may introduce highly deformed cells for combinations of large translational and rotational motions. The overlooked method allows for arbitrarily large body motions and retains the quality of the mesh. However, it comes with a substantial increase in computational cost and possible loss of energy conservation due to the interpolation. In this paper we present a straightforward extension of the spherical linear interpolation (SLERP) based mesh morphing algorithm that increases the stability range of the method. The mesh deformation is allowed to be interpolated independently for different modes of motion, which facilitates tailored mesh motion simulations. The paper details the implementation of the method and evaluates its performance with computational examples of a cylinder with a moonpool. The examples show that the modified mesh morphing approach handles large motions well and provides a cost effective alternative to overlooked mesh for survival conditions.
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.
Coupled mooring analysis using CFD with dynamic mooring models is becoming an established field. This is an important step for better predictions of responses of moored marine structures in extreme sea states and also for capturing the low-frequency response correctly. The coupling between the CFD and mooring solvers are most often carried out by exchanging the fairlead/anchor points and fairlead forces. In this paper we will discuss the effects of using (i) viscous fluid flow on a mooring component level (submerged buoys and clump weights) and (ii) the fluid-structure coupling between the viscous fluid solver and the mooring system.
Floating Power Plant (FPP) develops a hybrid floating wind and wave energy device. Pitching Wave Energy Converters (WECs) interact with the supporting structure, amplifying the motion of the WECs within the design wave frequency range. In this work we focus on the effect of the chamber geometry – without the WEC – in amplifying the waves inside the chamber. The simulations are carried out using two-phase Navier-Stokes simulations. We investigate the wave propagation and the interaction between waves and the fixed support structure. The simulations are compared to experimental tests performed in the wave basin at Aalborg University.
The International Energy Agency Technology Collaboration Program for Ocean Energy Systems (OES) initiated the OES Wave Energy Conversion Modeling Task, which focused on the verification and validation of numerical models for simulating wave energy converters (WECs). The long-term goal is to assess the accuracy of and establish confidence in the use of numerical models used in design as well as power performance assessment of WECs. To establish this confidence, the authors used different existing computational modeling tools to simulate given tasks to identify uncertainties related to simulation methodologies: (i) linear potential flow methods; (ii) weakly nonlinear Froude–Krylov methods; and (iii) fully nonlinear methods (fully nonlinear potential flow and Navier–Stokes models). This article summarizes the code-to-code task and code-to-experiment task that have been performed so far in this project, with a focus on investigating the impact of different levels of nonlinearities in the numerical models. Two different WECs were studied and simulated. The first was a heaving semi-submerged sphere, where free-decay tests and both regular and irregular wave cases were investigated in a code-to-code comparison. The second case was a heaving float corresponding to a physical model tested in a wave tank. We considered radiation, diffraction, and regular wave cases and compared quantities, such as the WEC motion, power output and hydrodynamic loading.
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.
This paper analyzes the nonlinear forces on a moored point-absorbing wave energy converter (WEC) in resonance at prototype scale (1:1) and at model scale (1:16). Three simulation types were used: Reynolds Averaged Navier-Stokes (RANS), Euler and the linear radiation-diffraction method (linear). Results show that when the wave steepness is doubled, the response reduction is: (i) 3% due to the nonlinear mooring response and the Froude-Krylov force; (ii) 1-4% due to viscous forces; and (iii) 18-19% due to induced drag and non-linear added mass and radiation forces. The effect of the induced drag is shown to be largely scale-independent. It is caused by local pressure variations due to vortex generation below the body, which reduces the total pressure force on the hole. Euler simulations are shown to be scale-independent and the scale effects of the WEC are limited by the purely viscous contribution (1-4%) for the two waves studied. We recommend that experimental model scale test campaigns of WECs should be accompanied by RANS simulations, and the analysis complemented by scale-independent Euler simulations to quantify the scale-dependent part of the nonlinear effects.
The wave loads and the resulting motions of floating wave energy converters are traditionally computed using linear radiation–diffraction methods. Yet for certain cases such as survival conditions, phase control and wave energy converters operating in the resonance region, more complete mathematical models such as computational fluid dynamics are preferred and over the last 5 years, computational fluid dynamics has become more frequently used in the wave energy field. However, rigorous estimation of numerical errors, convergence rates and uncertainties associated with computational fluid dynamics simulations have largely been overlooked in the wave energy sector. In this article, we apply formal verification and validation techniques to computational fluid dynamics simulations of a passively controlled point absorber.
The phase control causes the motion response to be highly nonlinear even for almost linear incident waves. First, we show that the computational fluid dynamics simulations have acceptable agreement to experimental data. We then present a verification and validation study focusing on the solution verification covering spatial and temporal discretization, iterative and domain modelling errors. It is shown that the dominating source of errors is, as expected, the spatial discretization, but temporal and iterative errors cannot be neglected. Using hexahedral cells with low aspect ratio and 30 cells per wave height, we obtain results with less than 5% uncertainty in motion response (except for surge) and restraining forces for the buoy without phase control. The amplified nonlinear response due to phase control caused a large increase in numerical uncertainty, illustrating the difficulty to obtain reliable solutions for highly nonlinear responses, and that much denser meshes are required for such cases.
Computational fluid dynamics (CFD) is becoming an increasingly popular tool in the wave energy sector, and over the last five years we have seen many studies using CFD. While the focus of the CFD studies have been on the validation phase, comparing numerically obtained results against experimental tests, the uncertainties associated with the numerical solution has so far been more or less overlooked. There is a need to increase the reliability of the numerical solutions in order to perform simulation based optimization at early stages of development. In this paper we introduce a well-established verification and validation (V&V) technique. We focus on the solution verification stage and how to estimate spatial discretization errors for simulations where no exact solutions are available. The technique is applied to the cases of a 2D heaving box and a 3D moored cylinder. The uncertainties are typically acceptable with a few percent for the 2D box, while the 3D cylinder case show double digit uncertainties. The uncertainties are discussed with regard to physical features of the flow and numerical techniques.
Floating wave energy converters (WECs) operating in the resonance region are strongly affected by non-linearities arising from the interaction between the waves, the WEC motion and the mooring restraints. To compute the restrained WEC motion thus requires a method which readily accounts for these effects. This paper presents a method for coupled mooring analysis using a two-phase Navier-Stokes (VOF-RANS) model and a high-order finite element model of mooring cables. The method is validated against experimental measurements of a cylindrical buoy in regular waves, slack-moored with three catenary mooring cables. There is overall a good agreement between experimental and computational results with respect to buoy motions and mooring forces. Most importantly, the coupled numerical model accurately recreates the strong wave height dependence of the response amplitude operators seen in the experiments.