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.
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.