A free-surface-only closure model for linear waves on general deep-water flows

Abstract

We present a novel, spatially two-dimensional (2D) set of equations to study the propagation of linear deep-water surface waves over a general steady three-dimensional (3D) background free surface flow. No assumptions of a flat background free surface, nor of an irrotational background flow, are made. The resulting model is not only a significant theoretical simplification, but also results in orders of magnitude faster computations. The linearized Euler equations are evaluated on the free surface of the base flow, and a closure condition is proposed to account for the vertical derivatives at the free surface. This generalizes a recent result for purely rotating background flows (Zuccoli, Brambley&Barkley, 2025, arXiv:2405.12078). The final model consists of five coupled first order partial differential equations (PDEs) to be solved on the 2D free surface, involving five unknowns: the horizontal and vertical perturbation velocities; the free surface perturbation height; and unexpectedly the gradient of perturbation pressure with depth at the surface. Two test problems are used to validate the model: a uni-directional vortical depth-varying base flow with a flat free surface; and perturbations to a travelling Gerstner wave solution. Eigenvalues and eigenfunctions of the linearized Euler equations are computed and compared with those of the model. Results show remarkable agreement between the two. Our study finds for the first time, to the best of our knowledge, that sufficiently steep two-dimensional Gerstner waves are unstable.