“…Since the seminal work of Longuet-Higgins and Cokelet [75], increasingly efficient and generic numerical models were developed to simulate nonlinear ocean waves based on fully nonlinear potential flow (FNPF) equations, i.e., Euler equations for irrotational flows, in two-(2D), e.g., [2,13,19,21,22,35,37,38,40,43,44,46,47,58,65,66,81,84,88,111] or three-dimensions (3D), e.g., [5,11,25,26,34,55,82,83,94,114]. When simulating overturning waves in these models, as first proposed by [75], the free surface motion is typically described by an Eulerian-Lagrangian formulation.…”