Linear potential flow (LPF) models remain the tools-of-the-trade in marine and ocean engineering despite their well-known assumptions of small amplitude waves and motions. As of now, nonlinear simulation tools are still too computationally demanding to be used in the entire design loop, especially when it comes to the evaluation of numerous irregular sea states. In this paper we aim to enhance the performance of the LPF models by introducing a hybrid LPF-ML (machine learning) approach, based on identification of nonlinear force corrections. The corrections are defined as the difference in hydrodynamic force (viscous and pressure-based) between high-fidelity CFD and LPF models. Using prescribed chirp motions with different amplitudes, we train a long short-term memory (LSTM) network to predict the corrections. The LSTM network is then linked to the MoodyMarine LPF model to provide the nonlinear correction force at every time-step, based on the dynamic state of the body and the corresponding forces from the LPF model. The method is illustrated for the case of a heaving sphere in decay, regular and irregular waves – including passive control. The hybrid LPF model is shown to give significant improvements compared to the baseline LPF model, even though the training is quite generic.
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 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.
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.
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.
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.
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.
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.
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.
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.