A generic point-absorbing wave energy converter is modeled in CFD as a vertical cylinder, moored with a single catenary chain that is fully coupled through a dynamic mooring code. The method of choice is very complete and takes much of the non-linearities in the highly coupled system of the moored body into account. The paper presents numerical results compared with experimental data for surge, heave and pitch motion in both decay tests and regular waves. Further, the wave motion response of the cylinder is computed using both a viscous and a non-viscous formulation as a first attempt to quantify viscous effects. Results show a good match between numerical and experimental results in heave, while the surge and pitch motion are more difficult to reproduce. The mooring load cycle appearance compares well with the experiments in shape but gives higher peak values. Although made at low Keulegan-Carpenter numbers, the simulations show vortical structures due to the heave motion, and the resulting motions are clearly affected by the inclusion or exclusion of viscosity. More test-cases and detailed experimental results are needed for further quantification of the viscous impact on floating point absorbers.
Accurate prediction of wave transformation is key in the design of coastal and nearshore structures which typically depends on numerical models. Turbulent and rotational effects call for the use of Computational Fluid Dynamics (CFD) solvers of which a large range of formulations including free surface treatments exists. Physical wave flume tests of wave propagation over a submerged bar with various levels of nonlinearity, regularity, and wave-breaking, dedicated to numerical model benchmarking or validation, were carried out in the Ocean and Coastal Engineering Laboratory of Aalborg University. Three fundamentally different CFD models each widespread within their category are benchmarked against the experimental data. The CFD models are based on (i) the Volume of Fluid (VoF) based interFoam solver of OpenFOAM, (ii) the sigma-transformation solver of MIKE 3 Waves Model FM, and (iii) the weakly compressible delta-SPH solver of DualSPHysics. Accuracy of the numerical models is assessed from surface elevation time series, evaluation metrics (averaged errors on surface elevations, amplitudes, phases, and wave set-up), and spectral analyzes to calculate the amplitude and phase contents of primary and higher-order components along the wave flume. Applicability is assessed from computational costs and ease-of-use factors such as the effort to configure the numerical models and achieve convergence. In general, the numerical models have high correlation to the physical tests and are as such suitable to model complex wave transformation with an accuracy sufficient for most coastal engineering applications. The VoF model performs more accurately under the turbulent conditions of breaking waves, increasing its relative accuracy in the prediction of downwave surface elevation. The sigma transformation model has simulation times one to two orders of magnitude lower than those of the VoF and SPH models.
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.
The paper presents incompressible Navier-Stokes simulations of the dynamics of a floating wave energy converter (WEC) coupled to a high-order finite element solver for cable dynamics. The coupled model has very few limiting assumptions and is capable of capturing the effects of breaking waves, green water loads on the WEC as well as non-linear mooring forces and snap loads, all of which are crucial for correct estimates of the extreme loads acting on the system in violent seas. The cable dynamics model has been developed as a stand-alone library that can be coupled to any body motion solver. In this study the open-source CFD package OpenFOAM has been employed. Preliminary test cases using incident regular Stoke's 5th order waves are presented, both for wave heights corresponding to operational conditions of the WEC as for a more severe condition in survival mode. It is illustrated that the coupled model is able to capture the complicated force propagation in the mooring cables.
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.
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 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.
Simulating the free decay motion and wave radiation from a heaving semi-submerged sphere poses significant computational challenges due to its three-dimensional complexity. By leveraging axisymmetry, we reduce the problem to a two-dimensional simulation, significantly decreasing computational demands while maintaining accuracy. In this paper, we exploit axisymmetry to perform a large ensemble of Computational Fluid Dynamics (CFDs) simulations, aiming to evaluate and maximize both accuracy and efficiency, using the Reynolds Averaged Navier–Stokes (RANS) solver interFOAM, in the opensource finite volume CFD software OpenFOAM. Validated against highly accurate experimental data, extensive parametric studies are conducted, previously limited by computational constraints, which facilitate the refinement of simulation setups. More than 50 iterations of the same heaving sphere simulation are performed, informing efficient trade-offs between computational cost and accuracy across various simulation parameters and mesh configurations. Ultimately, by employing axisymmetry, this research contributes to the development of more accurate and efficient numerical modeling in ocean engineering.
As more offshore wind energy projects are implemented, the risk of interactions between farms becomes more pronounced. While reduced surface roughness over water enhances airflow stability, it can also extend wake effects on downstream turbines. The study aims to enhance the understanding of wake interactions and efficiency variations based on the distance between neighboring farms. To assess the impact of neighboring farms across different scenarios and features, a methodology is developed to achieve computational optimality using an open-source Python-based library, PyWake, then verified by a well-established CFD software, Meteodyn. Then, the methodology is applied to a Brazilian offshore wind project currently under licensing as a reference point. The results indicate a 1–3% reduction in Annual Energy Production following the current Brazilian regulation for onshore projects of 20 times the blade tip height, as the minimum distance. This reduction translates to an approximate 3% increase in the Levelized Cost of Energy and a nearly 24% decrease in Net Present Value. These findings are crucial for offshore wind energy planning and its sustainable growth, indicating the need to define a minimum distance for the regulatory bodies. This would not only avoid future disputes but also enhance investor confidence.
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.