Modern design for marine and coastal activities places increasing focus on numerical simulations. Several numerical wave models have been developed in the past few decades with various techniques and assumptions. Those numerical models have their own advantages and disadvantages. The proper choice of the most useful numerical tool depends on the understanding of the validity and limitations of each model. In the past years, REEF3D has been developed into an open-source hydrodynamic numerical toolbox that consists of several modules based on the Navier–Stokes equations, the shallow water equations and the fully nonlinear potential theory. All modules share a common numerical basis which consists of rectilinear grids with an immersed boundary method, high-order finite differences and high-performance computing capabilities. The numerical wave tank of REEF3D utilises a relaxation method to generate waves at the inlet and dissipate them at the numerical beach. In combination with the choice of the numerical grid and discretisation methods, high accuracy and stability can be achieved for the calculation of free surface wave propagation and transformation. The comparison among those models provide an objective overview of the different wave modelling techniques in terms of their numerical performance as well as validity. The performance of the different modules is validated and compared using several benchmark cases. They range from simple propagations of regular waves to three-dimensional wave breaking over a changing bathymetry. The diversity of the test cases help with an educated choice of wave models for different scenarios.

The experimental wave paddle signal is unknown to the numerical modellers in many cases. This makes it quite challenging to numerically reproduce the time history of free surface elevation for irregular waves. In the present work, a numerical investigation is performed using a computational fluid dynamics (CFD) based model to validate and investigate a non-iterative free surface reconstruction technique for irregular waves. In the current approach, the free surface is reconstructed by spectrally composing the irregular wave train as a summation of the harmonic components coupled with the Dirichlet inlet boundary condition. The verification is performed by comparing the numerically reconstructed free surface elevation with theoretical input waves. The applicability of the present approach to generate irregular waves by reconstructing the free surface is investigated for different coastal and marine engineering problems. A numerical analysis is performed to validate the free surface reconstruction approach to generate breaking irregular waves over a submerged bar. The wave amplitudes, wave frequencies and wave phases are modelled with good accuracy in the time-domain during the higher-order energy transfers and complex processes like wave shoaling, wave breaking and wave decomposition. The present approach to generate irregular waves is also employed to model steep irregular waves in deep water. The free surface reconstruction method is able to simulate the irregular free surface profiles in deep water with low root mean square errors and high correlation coefficients. Furthermore, the irregular wave forces on a monopile are investigated in the time-domain. The amplitudes and phases of the force signal under irregular waves generated by using the current technique are modelled accurately in the time-domain. The proposed approach to numerically reproduce the free surface elevation in the time-domain provides promising and accurate results for all the benchmark cases.

A comprehensive understanding of the marine environment in the offshore area requires phase-resolved wave information. For far-field wave propagation, computational efficiency is crucial, as large spatial and temporal scales are involved. For the near-field extreme wave events and wave impacts, high resolution is required to resolve the flow details and turbulence. The combined use of a computationally efficient large-scale model and a high-resolution local-scale solver provides a solution that combines accuracy and efficiency. This article introduces a coupling strategy between the efficient fully nonlinear potential flow (FNPF) solver REEF3D::FNPF and the high-fidelity computational fluid dynamics (CFD) model REEF3D::CFD within the open-source hydrodynamics framework REEF3D. REEF3D::FNPF solves the Laplace equation together with the boundary conditions on a sigma-coordinate. The free surface boundary conditions are discretized using high-order finite difference methods. The Laplace equation for the velocity potential is solved with a conjugated gradient solver preconditioned with a geometric multigrid provided by the open-source library Hypre. The model is fully parallelized following the domain decomposition strategy and the message passing interface protocol. The waves calculated with the FNPF solver are used as wave generation boundary conditions for the CFD-based numerical wave tank REEF3D::CFD. The CFD model employs an interface capturing two-phase flow approach that can resolve complex wave structure interaction, including breaking wave kinematics and turbulent effects. The presented hydrodynamic coupling strategy is tested for various wave conditions and the accuracy is fully assessed.