Arxiv:1903.11909V1 [Physics.Flu-Dyn] 28 Mar 2019

Exact solution for progressive gravity waves on the surface of a deep fluid Nail S. Ussembayev∗ Computer, Electrical, and Mathematical Sciences and Engineering Divison King Abdullah University of Science and Technology, Thuwal 23955-6900, KSA Gerstner or trochoidal wave is the only known exact solution of the Euler equations for periodic surface gravity waves on deep water. In this Letter we utilize Zakharov’s variational formulation of weakly nonlinear surface waves and, without truncating the Hamiltonian in its slope expansion, derive the equations of motion for unidirectional gravity waves propagating in a two-dimensional flow. We obtain an exact solution of the evolution equations in terms of the Lambert W -function. The associated flow field is irrotational. The maximum wave height occurs for a wave steepness of 0.2034 which compares to 0.3183 for the trochoidal wave and 0.1412 for the Stokes wave. Like in the case of Gerstner’s solution, the limiting wave of a new type has a cusp of zero angle at its crest. Background.—The water waves problem even in the tial evaluated on the free surface. The expansion coef- completely idealized setting of a perfect fluid is notori- ficients simplify remarkably for two-dimensional waves ously difficult from a mathematical point of view. This propagating in the same direction and the infinite series is one of the reasons we have very few explicit solutions can be summed in a closed form if the wave steepness to the fully nonlinear free-surface hydrodynamics equa- is small. From the stationary periodic solution of the tions despite extensive research over hundreds of years. resulting nonlinear equations, an exact irrotational solu- The first exact solution was provided by Gerstner [1] in tion to the full system of hydrodynamic equations with a 1802 for pure gravity waves of finite height propagating free surface is then recovered. The solutions are provided on the surface of an inviscid and incompressible fluid of in Eulerian coordinates [10]. infinite depth. The practical importance of trochoidal Exact solutions play an important role not just because waves became apparent when in the 1860s Rankine and all the properties of the motion can be conveniently ex- Froude independently rediscovered them in connection pressed and analyzed without any approximation, but with the ship-rolling problem [2]. The fact that these also because they can provide a key insight into the na- waves are rotational and cannot be generated from the ture of more complex and physically more realistic flows, fluid at rest was seen as somewhat unsatisfactory. More- say, in which waves propagate under the combined influ- over, since the fluid particles move on a closed circular ence of gravity and surface tension. orbit, Gerstner’s solution does not predict the observed Let us finish this introduction with a note that two- mass transport in the direction of wave propagation [3]. dimensional flows are not of purely theoretical interest. In 1847 Stokes presented an approximate theory for ir- Their fundamental importance stems from the observa- rotational waves and showed that for deep water the wave tion that in many situations the motion of water waves profile coincides with a trochoid to the third order, but can be treated as quasi-two-dimensional. For example, not to the fourth [4]. He also conjectured the existence waves generated in the open ocean while coming ashore of the wave of greatest height distinguished by a sharp often appear as a wave-train of very long nearly paral- edge and argued that the enclosing angle must be 120◦. lel crests [11] and due to this symmetry can be more The mere existence of such waves was rigorously verified or less accurately studied and simulated by various two- only forty years ago [5], clearly illustrating the inher- dimensional models (see, e.g., [12]). ently complicated structure of the underlying nonlinear Theoretical formulation.—The potential flow φ(x; y; t) equations. The existence of extreme waves that are, as of an ideal incompressible fluid in the region Ω(t) = 2 hypothesized by Stokes, strictly convex between succes- (x; y) 2 R : − 1 < y ≤ ξ(x; t) bounded by the free sive crests was established very recently, but in spite of surface Γ(t) = fy = ξ(x; t)g satisfies the following set of some numerical evidence supporting the convexity of all equations: waves of extreme form their uniqueness remains an open problem [6]. Apart from Crapper’s exact solution for pure capillary ξt + ξxφx − φy = 0; on Γ(t) (1a) waves [7] found in 1957 and some explicit Gerstner-type arXiv:1903.11909v1 [physics.flu-dyn] 28 Mar 2019 1 2 2 solutions for equatorially trapped waves discovered lately φt + φ + φ + gξ = 0; on Γ(t) (1b) 2 x y [8], the field of exact solutions for water wave equations φ + φ = 0; in Ω(t) (1c) has not experienced significant advancement since the xx yy last century. In this Letter, we present an exact solution φy ! 0; as y ! −∞ (1d) for planar periodic waves of finite amplitude moving uni- formly and without change of form on a fluid of unlim- where g is the constant gravitational acceleration in the ited depth when only gravity acts as a restoring force. negative y direction and subscripts denote partial deriva- Following Zakharov [9], we use the power of Hamilto- tives. The effects of surface tension are neglected and nian formalism to write an asymptotic expansion of the the pressure at the free surface is taken to be zero. Once equations of motion for a pair of canonically conjugate ξ(x; t) and φ(x; y; t) are found, the velocity field and the variables – the surface elevation and the velocity poten- pressure within the fluid domain follow from 2 The major difficulty in applying the Hamiltonian formalism for water waves lies in expressing the kinetic energy in terms of free surface variables, ξ and , alone. v(x; y; t) = rφ, When the wave steepness is small, jξxj 1, one can write jvj2 the expansion of the Hamiltonian in a power series in the p(x; y; t) = −ρ gy + φt + ^ ^∗ ^ ^∗ 2 amplitudes of the waves. Let ξk = ξ−k and k = −k denote the Fourier components of ξ and , then upon solving the boundary-value problem for Laplace’s equa- where ρ is the fluid density. Since φ is harmonic in Ω(t), it tion the Hamiltonian (3) reads is uniquely determined by its trace on the free boundary, (x; t) = φj , and (1d). Zakharov ingeniously real- y=ξ(x;t) Z ized that equations (1a)-(1b) describing the wave motion 1 ^∗ ^ ^∗ ^ H = jkj 0 0 + gξ0 ξ0dk+ can be put in the form of Hamilton’s canonical equations 2 1 Z n+1 n+1 ! n+1 X (n+2) ^ ^ Y ^ X Y + E0;1;2;:::;n+1 0 1 ξiδ ki dkj: n=1 i=2 i=0 j=0 @ξ δH @ δH ρ = ; ρ = − (2) (4) @t δ @t δξ where we set ρ = 1 [13]. Computing the expansion ker- (n+2) with the Hamiltonian given as the sum of the potential nels E0;1;2;:::;n+1 relies on a series reversion associated energy measured with respect to the undisturbed fluid with expressing the velocity potential in terms of its value level y = 0 and the kinetic energy [9]: at the free surface [14]. These kernels can be obtained in principle to all orders, however, the iteration process quickly becomes cumbersome. A fully explicit recursion ! ρ Z 1 Z ξ relation for computing E(n+2) without a laborious H = dx gξ2 + jvj2dy : (3) 0;1;2;:::;n+1 2 −∞ −∞ series reversion has been proposed recently in [15] for an arbitrary depth and arbitrary space dimension. The Fourier transform is a canonical transformation, The Laplace equation (1c) and the decay condition at and hence the evolution equations (2) preserve their infinity (1d) represent the constraints on the dynamical Hamiltonian form in the canonically conjugate variables ^ ^ system whose total energy H is a constant of motion. ξk and k: ^ 1 Z n+1 n+1 ! @ξk δH ^ X (n+2) ^ Y ^ X = = jkj k + 2 E 1 ξiδ k − ki dk12:::n+1; (5a) @t ^∗ −0;1;2;:::;n+1 δ k n=1 i=2 i=1 ^ 1 Z n+1 n+1 ! @ k δH ^ X (n+2) ^ ^ Y ^ X = − = −gξk − nE 1 2 ξiδ k − ki dk12:::n+1 (5b) @t ^∗ 1;2;:::;n+1;−0 δξk n=1 i=3 i=1 where (4) was used to compute the variational derivative. The expansion kernels posses a remarkable property re- ξ± ± i ± = 0 (8a) vealed by the recursion relation: for any fixed n 2 N, it t x holds that ± 2 ± ± (i x ) t + gξ = ± (8b) 1 ∓ iξx (n+2) E Pn+1 = 0 (6) − i=1 ki;k1;k2;:::;kn+1 ± ± ± 1 where f = P f, P = 2 (I ± iH) are orthogonal self- adjoint projections onto positive and negative wavenum- if all wavenumbers have the same sign [15] and ber components, H is the Hilbert transform [16]. In sum- mary, the small-amplitude assumption, jξxj 1, in con- n+1 n+1 junction with the unidirectionality of wave propagation (n+2) σ Y E = − k (7) Pn+1 n=2 i for which relations (6)-(7) are valid in one dimension al- k1;k2;:::;kn+1;− i=1 ki n(2π) i=1 lows us to write apparently intangible equations (5) con- taining infinite sums in a concise form given in (8). with σ = −1 (σ = +1) if all wavenumbers are positive A few remarks are in order.

Load more