Wave breaking in the ocean affects the height of extreme waves, energy dissipation, and interaction between the atmosphere and upper ocean. Numerical modelling is a critical step in understanding the physics of wave breaking and offers insight that is hard to gain from field data or experiments. High-fidelity numerical modelling of three-dimensional breaking waves is extremely challenging. Conventional grid-based numerical methods struggle to model the steep and double-valued free surfaces that occur during wave breaking. The Smoothed Particle Hydrodynamics (SPH) method does not fall prey to these issues. Herein, we examine the SPH method's ability to model highly directionally spread overturning breaking waves by numerically reproducing the experiments presented in McAllister et al. (2019). We find that the SPH method reproduces the experimental observations well; when comparing experimental and numerical measurements we achieve coefficient of determination values of 0.92−0.95, with some smaller-scale features less well reproduced owing to finite resolution. We also examine aspects of the simulated wave's geometry and kinematics and find that existing breaking criteria are difficult to apply in highly directionally spread conditions.

Freak waves, abnormally large waves, that occur in the open-ocean can cause significant damage to offshore structures and vessels. In this paper, we attempt to numerically reproduce the experiments of McAllister et al., (2019, J. Fluid Mech. [1]), to investigate the potential properties of the Draupner freak wave [2] in more detail. We use a Smoothed Particle Hydrodynamics (SPH) method to solve the full-3D Navier-Stokes equations. This Lagrangian method is able to recreate wave breaking, and has the potential to fully reproduce these experiments with the aim of providing further insight into properties of the waves created such as their kinematics and geometry. We compare time histories of water surface elevation produced numerically using four different particle sizes with experimentally-obtained data. We find good agreement in the time domain, with r2 (coefficient of determination) values between experimental and numerical data of over 0.94 the error in maximum wave height was less than 5 % for the finest particle size (over 100 million particles). We also numerically reproduce wave breaking observed in the experiments, where jet formation and breaking phenomena are qualitatively similar in appearance.