We present a Spectral Element Fully Nonlinear Potential Flow (FNPF-SEM) model developed for the simulation of wave-body interactions between nonlinear free surface waves and impermeable structures. The solver is accelerated using an iterative p-multigrid algorithm. Two cases are considered: (i) a surface piercing box forced into vertical motion creating radiated waves and (ii) a rectangular box released above its equilibrium resulting in freely decaying heave motion. The FNPF-SEM model is validated by comparing the computed hydrodynamic forces against those obtained by a Navier-Stokes solver. Although not perfect agreement is observed the results are promising, a significant speedup due to the iterative algorithm is however seen.
We present a depth-integrated Boussinesq model for the efficient simulation of nonlinear wave–body interaction. The model exploits a ‘unified’ Boussinesq framework, i.e. the fluid under the body is also treated with the depth-integrated approach. The unified Boussinesq approach was initially proposed by Jiang (2001) and recently analyzed by Lannes (2017). The choice of Boussinesq-type equations removes the vertical dimension of the problem, resulting in a wave–body model with adequate precision for weakly nonlinear and dispersive waves expressed in horizontal dimensions only. The framework involves the coupling of two different domains with different flow characteristics. Inside each domain, the continuous spectral/hp element method is used to solve the appropriate flow model since it allows to achieve high-order, possibly exponential, convergence for non-breaking waves. Flux-based conditions for the domain coupling are used, following the recipes provided by the discontinuous Galerkin framework. The main contribution of this work is the inclusion of floating surface-piercing bodies in the conventional depth-integrated Boussinesq framework and the use of a spectral/hp element method for high-order accurate numerical discretization in space. The model is verified using manufactured solutions and validated against published results for wave–body interaction. The model is shown to have excellent accuracy and is relevant for applications of waves interacting with wave energy devices.
We present a high-order nodal spectral element method for the two-dimensional simulation of nonlinear water waves. The model is based on the mixed Eulerian–Lagrangian (MEL) method. Wave interaction with fixed truncated structures is handled using unstructured meshes consisting of high-order iso-parametric quadrilateral/triangular elements to represent the body surfaces as well as the free surface elevation. A numerical eigenvalue analysis highlights that using a thin top layer of quadrilateral elements circumvents the general instability problem associated with the use of asymmetric mesh topology.We demonstrate how to obtain a robust MEL scheme for highly nonlinear waves using an efficient combination of (i) global L2 projection without quadrature errors, (ii) mild modal filtering and (iii) a combination of local and global re-meshing techniques. Numerical experiments for strongly nonlinear waves are presented. The experiments demonstrate that the spectral element model provides excellent accuracy in prediction of nonlinear and dispersive wave propagation. The model is also shown to accurately capture the interaction between solitary waves and fixed submerged and surface-piercing bodies. The wave motion and the wave-induced loads compare well to experimental and computational results from the literature.
For the assessment of experimental measurements of focused wave groups impacting a surface-piecing fixed structure, we present a new Fully Nonlinear Potential Flow (FNPF) model for simulation of unsteady water waves. The FNPF model is discretized in three spatial dimensions (3D) using high-order prismatic - possibly curvilinear - elements using a spectral element method (SEM) that has support for adaptive unstructured meshes. This SEM-FNPF model is based on an Eulerian formulation and deviates from past works in that a direct discretization of the Laplace problem is used making it straightforward to handle accurately floating structural bodies of arbitrary shape. Our objectives are; i) present detail of a new SEM modelling developments and ii) to consider its application to address a wave-body interaction problem for nonlinear design waves and their interaction with a model-scale fixed Floating Production, Storage and Offloading vessel (FPSO). We first reproduce experimental measurements for focused design waves that represent a probably extreme wave event for a sea state represented by a wave spectrum and seek to reproduce these measurements in a numerical wave tank. The validated input signal based on measurements is then generated in a NWT setup that includes the FPSO and differences in the signal caused by nonlinear diffraction is reported.
We present an arbitrary-order spectral element method for general-purpose simulation of non-overturning water waves, described by fully nonlinear potential theory. The method can be viewed as a high-order extension of the classical finite element method proposed by Cai et al. (1998)[5], although the numerical implementation differs greatly. Features of the proposed spectral element method include: nodal Lagrange basis functions, a general quadrature-free approach and gradient recovery using global L2projections. The quartic nonlinear terms present in the Zakharov form of the free surface conditions can cause severe aliasing problems and consequently numerical instability for marginally resolved or very steep waves. We show how the scheme can be stabilised through a combination of over-integration of the Galerkin projections and a mild spectral filtering on a per element basis. This effectively removes any aliasing driven instabilities while retaining the high-order accuracy of the numerical scheme. The additional computational cost of the over-integration is found insignificant compared to the cost of solving the Laplace problem. The model is applied to several benchmark cases in two dimensions. The results confirm the high order accuracy of the model (exponential convergence), and demonstrate the potential for accuracy and speedup. The results of numerical experiments are in excellent agreement with both analytical and experimental results for strongly nonlinear and irregular dispersive wave propagation. The benefit of using a high-order – possibly adapted – spatial discretisation for accurate water wave propagation over long times and distances is particularly attractive for marine hydrodynamics applications.