Spatial Form of a Hamiltonian Dysthe Equation for Deep-Water Gravity Waves

An overview of a Hamiltonian framework for the description of nonlinear modulation of surface water waves is presented. The main result is the derivation of a Hamiltonian version of Dysthe’s equation for two-dimensional gravity waves on deep water. The reduced problem is obtained via a Birkhoff normal form transformation which not only helps eliminate all non-resonant cubic terms but also yields a non-perturbative procedure for surface reconstruction. The free surface is reconstructed from the wave envelope by solving an inviscid Burgers’ equation with an initial condition given by the modulational Ansatz. Particular attention is paid to the spatial form of this model, which is simulated numerically and tested against laboratory experiments on periodic groups and short-wave packets. Satisfactory agreement is found in all these cases.


Introduction
Because the cubic nonlinear Schrödinger (NLS) equation for water waves had been shown to have a limited range of applicability, Dysthe [1] extended this model to the next order in the perturbation analysis. The resulting equation is now widely used in the water wave community because of its efficiency and ability to describe realistic waves, in particular with moderately large steepnesses. Among the new features are effects from the wave-induced mean flow which lead to improvement in the stability properties of finite-amplitude waves. Due to its higher-order nonlinear terms, Dysthe's equation is also able to capture the asymmetric development of propagating wave packets. Originally derived in the context of gravity waves on deep water, it has since been extended to various other settings and has produced a large literature: e.g., finite depth [2], gravitycapillary waves [3], exact linear dispersion [4], dissipation [5], higher order [6], just to cite a few references. Interest in Dysthe's equation lies in the fact that it is more amenable to mathematical analysis and numerical simulation than the full water wave equations, while being able to capture salient features of weakly nonlinear wave packets. In particular, the associated computational cost is significantly less compared to that from direct (fully nonlinear) numerical solvers [7][8][9][10].
In this asymptotic regime, solutions are sought as mild modulations of monochromatic waves. Besides the wave steepness, another small parameter is associated with the spectral width relative to the dominant carrier wavenumber (the so-called narrowband approximation), and accordingly evolution equations are written for the complex wave envelope. Such models have typically been derived via the method of multiple scales from the full water wave equations [1,2,5], or as a reduction of the Fourier integro-differential equation that was first proposed by Zakharov [11] and later refined by Krasitskii [12] for resonantly interacting water waves in the small-amplitude limit [3,4,13]. These studies provide a formal derivation and analysis of Dysthe's equation, and it is only recently that rigorous results on its well-posedness have been obtained [14,15]. The rigorous justification of envelope equations for water waves has so far been proved for the NLS equation only [16,17]. This mathematical question remains an open problem for Dysthe's equation.
Despite these advances, early versions of Dysthe's equation share a fundamental shortcoming: unlike the NLS Equation [18], they are not Hamiltonian partial differential equations whereas the full water wave equations are. From a modeling point of view, it is desirable that a reduced model inherits important properties of the original system. Regardless of whether the modulational analysis is carried out from the full water wave equations or from Zakharov's equation, an envelope model is not guaranteed to be Hamiltonian unless there is a way to ensure that a symplectic structure is retained together with some conserved energy.
Progress in this direction has recently been made by Gramstad and Trulsen [19] and Craig et al. [20] who proposed Hamiltonian versions of Dysthe's equation. Gramstad and Trulsen [19] used Krasitskii's refined version of Zakharov's equation as a basis for their modulational analysis of three-dimensional gravity waves on arbitrary depth. Craig et al. [20] considered the two-dimensional problem of gravity waves on deep water and combined results from [21][22][23][24][25] to derive such models. Their formal approach involves a reduction of the basic Hamiltonian to Birkhoff normal form at fourth order, followed by a series of canonical transformations together with a procedure of homogenization. In particular, the normal form transformation helps eliminate all non-resonant cubic terms, which leads to a simple and elegant expression for the reduced Hamiltonian. These calculations are facilitated by the fact that the Dirichlet-Neumann operator has an explicit Taylor expansion, as devised by Craig and Sulem [7].
Another important step in the solution process for such approximations is the reconstruction of the free surface from the wave envelope, which is performed at a postprocessing stage but yet affects the predictions' accuracy. In classical (i.e., non-Hamiltonian) derivations of Dysthe's equation as well as in the alternate derivation by Gramstad and Trulsen [19], this reconstruction is carried out perturbatively in terms of a Stokes' expansion based on the Ansatz for the solutions' form. Contributions are included up to an order that is consistent with the order of approximation for Dysthe's equation. By contrast, Craig et al. [20] proposed a non-perturbative approach for surface reconstruction which requires solving a nonlinear partial differential equation (PDE) associated with the Birkhoff normal form transformation. This PDE for the surface elevation takes the generic form of an inviscid Burgers' equation. Higher-order harmonics of the surface wave spectrum are automatically generated from the carrier component through nonlinear interactions according to this PDE. As a consequence, the entire solution process in this approach fits within a Hamiltonian framework, meaning that both the envelope equation and the nonlinear PDE for surface reconstruction possess a well-defined Hamiltonian structure. Preliminary tests in [20] show very good agreement between this approximation and numerical simulations of the full water wave equations for the temporal evolution of Benjamin-Feir instability of Stokes waves on deep water. Furthermore, Craig et al. [20] provided an alternate spatial form of their Hamiltonian model in which the role of the spatial and temporal variables is interchanged. Spatial forms of Dysthe's equation have also received much attention in the literature as they are more suitable for comparison with time series recorded by fixed probes along the wave tank in laboratory experiments, or recorded by an array of buoys in field experiments [26][27][28]. Such a manipulation of the independent variables is applied to the envelope equation as well as to the procedure for surface reconstruction, which here requires a nontrivial adjustment of Burgers' equation associated with the Birkhoff normal form transformation.
In the present paper, after giving an overview of this Hamiltonian approach and reviewing recent results by Craig et al. [20], we propose an improved strategy for the surface computation. Focusing on the spatial dynamics (with the spatial and temporal variables switched), we further introduce a dimensionless rescaled form of the equations to facilitate comparison with laboratory experiments. We perform numerical simulations and test them under various wave conditions. On a related note, we also point out recent work by Dyachenko et al. [29][30][31] on a special reduction of Zakharov's equation for two-dimensional gravity waves on deep water (which has been referred to as the compact or super compact equation). This reduction is made possible in this case by taking advantage of explicit analytical expressions of 4-wave resonances and by identifying exact cancellation of terms. The resulting simpler form is similar to that for Dysthe's equation (while being free of the narrowband approximation) and still possesses a Hamiltonian structure. It is thus better suited than Zakharov's equation for analysis and computation, and numerical results have indicated that it is able to describe large-amplitude rogue waves and even waves reaching near the breaking point. A spatial Hamiltonian version of this compact equation can also be derived, as shown in [29,32].
The remainder of this paper is organized as follows. Section 2 introduces the basic mathematical formulation for two-dimensional gravity waves on deep water, including the Hamiltonian form of the governing equations and the Birkhoff normal form transformation to eliminate all non-resonant cubic terms. Section 3 describes the main steps in our Hamiltonian perturbation approach, including the modulational Ansatz to specify the special form of solutions in this asymptotic regime, the resulting Hamiltonian equation for temporal dynamics of the wave envelope, and the corresponding procedure for surface reconstruction. Section 4 presents the spatial version of this approximation together with its dimensionless rescaled form. Numerical simulations of our spatial model are then shown in Section 5 where it is compared to the classical formulation. Tests against laboratory experiments on periodic groups and short-wave packets are also presented, for which satisfactory agreement is obtained. Finally, concluding remarks are given in Section 6.

Hamiltonian Formulation
Consider a two-dimensional fluid of infinite depth bounded above by a free surface {y = η(x, t)}. The variables (x, y) denote the horizontal and vertical coordinates respectively, and the variable t denotes time. For an incompressible, inviscid and irrotational flow, the velocity potential ϕ(x, y, t) satisfies the boundary value problem where g is the acceleration due to gravity and ∇ = (∂ x , ∂ y ) is the spatial gradient. The Laplace problem (1)-(4) can be reduced from one posed inside the entire fluid domain to one posed on the free surface alone. As proposed by Craig and Sulem [7], a key ingredient is to introduce the Dirichlet-Neumann operator (DNO) which takes Dirichlet data ξ(x, t) = ϕ(x, η(x, t), t) at the free surface, solves Laplace's Equation (1) for ϕ with boundary condition (4), and returns the corresponding Neumann data (i.e., the surface normal velocity). Equations (1)-(4) then reduce to which constitute a closed Hamiltonian system for the two conjugate variables η and ξ in Zakharov's formulation of the water wave problem [7,11]. These can be cast into canonical form where the 2 × 2 matrix J represents the symplectic structure of the system, and the Hamiltonian coincides with the total energy and is an invariant for the temporal dynamics. All integrals are evaluated over R but this notation is omitted for convenience. Due to analyticity properties for a sufficiently regular free surface [33], the DNO can be written via a convergent Taylor series expansion about the quiescent state η = 0. Each term G j is homogeneous of degree j in η and is determined recursively [7]. The first three terms in this expansion are where D = −i∂ x (so that its Fourier symbol is k). The reader is referred to [9,21,24,[34][35][36] for other applications of this formulation to long-wave perturbation calculations as well as direct numerical simulations of nonlinear wave phenomena. With this observation, the Hamiltonian can be decomposed as represent the quadratic, cubic and quartic contributions respectively, in terms of (η, ξ). The first contribution H 2 describes the linear behavior while the next contributions H 3 and H 4 describe nonlinear mechanisms of 3-wave and 4-wave interactions, respectively. Our Hamiltonian approach amounts to applying a sequence of canonical transformations to the Hamiltonian (9) of the system, and to its symplectic structure (8). As a consequence, a Hamiltonian model is obtained at each level of approximation.

Normal Form Transformations
Because 3-wave resonances are not allowed for deep-water gravity waves, the Hamiltonian H(η, ξ) can be simplified by eliminating all non-essential cubic terms through a suitable canonical transformation [12]. This idea is inherent to the general theory of normal form transformations for Hamiltonian systems and has been successfully applied to the water wave problem, starting with the seminal work of Zakharov [11]. Birkhoff normal form up to a given order m results from a sequence of canonical changes of variables, so that the Taylor expansion of the transformed Hamiltonian up to order m contains only essential (i.e., resonant) terms.
In the two-dimensional problem of gravity waves on infinite depth [25], the normal form transformation for m = 3 can be defined as the flow at s = −1 subject to the vector field for some auxiliary Hamiltonian K 3 , and with initial condition at s = 0 being tied to the original physical variables (η, ξ). The variable s plays the role of an evolutionary variable for this auxiliary Hamiltonian system. More specifically, the normal form is defined through a Taylor expansion of the Hamiltonian about s = 0, where K 3 is homogeneous of degree 3, H n is homogeneous of degree n and the Poisson bracket is of degree n + 1. All cubic terms are eliminated if which is satisfied for Derivatives in the above expansion are related to Poisson brackets through by virtue of (11). The operator −i sgn(D) is the Hilbert transform in a Fourier multiplier form. Accordingly, it is convenient to introduce and The corresponding symplectic structure is given by J 0 = P 0 JP 0 = −J. This auxiliary system produces two equations: one for η, is an inviscid Burgers' equation, while the other for ξ, is its linearization along Burgers' flow. The new Hamiltonian then reads where the leading nonlinear contributions are now quartic terms, and follows from (12) together with the identity as shown in [25]. By construction, this normal form transformation is canonical and thus preserves the symplectic structure J in (8). More details on Birkhoff normal forms for the water wave problem can be found in [11,25,37,38].

Modulational Ansatz
Here we describe the special form of solutions to be sought in this asymptotic regime. We first make the change to complex symplectic coordinates z(x, t) and z(x, t) as defined by where the symbol. stands for complex conjugation, and System (8) then becomes after adjusting the symplectic structure J 1 = P 1 JP 1 . The next calculation involves the modulational Ansatz which implies that we are seeking solutions in the form of near-monochromatic waves with carrier wavenumber k 0 > 0 and with slowly varying complex envelope u(X, t) depending on X = εx. The dimensionless parameter ε = k 0 A 0 1 is a measure of the wave steepness (A 0 denotes a characteristic wave amplitude for the free surface). It is also a measure of the wave spectrum's narrowness around k = k 0 . The resulting equations of motion are where J 2 = εP 2 J 1 P 2 . The additional factor ε in this definition of J 2 reflects the change in symplectic structure from the spatial rescaling x → X = εx and plays a part in defining the relevant slow time scale for the resulting systems. More details on these canonical transformations in the context of Hamiltonian models for interfacial flows can be found in [22,24,39].

Envelope Equation for Temporal Dynamics
Naturally, the Hamiltonian (16) is also modified through the changes of variables (17) and (19). The first transformation (17) performs a diagonalization of the quadratic part in terms of normal modes (z, z) and the linear dispersion relation for deep-water gravity waves. The second transformation (19) sets the stage for the expansion of H in powers of ε. These calculations are facilitated by the fact that the DNO admits a Taylor expansion with an explicit form for its Taylor terms.
Due to the multiscale nature of this problem (fast oscillations in x and slow modulation in X), it is important to understand the action of Fourier multiplier operators on multiscale functions. The following result from [40] describes this action. Theorem 1. Let α(D) be a Fourier multiplier. For sufficiently smooth functions f (X), we have As an example, Furthermore, the presence of multiple scales needs to be appropriately dealt with in order to extract the nontrivial dynamics of the slowly varying wave envelope. The scale separation lemma of Craig et al. [21,23] is used to address this homogenization problem for our Hamiltonian system. As a result, terms with fast oscillations essentially homogenize to zero and thus do not contribute to the effective Hamiltonian, which implies that only 4-wave resonances (i.e., terms with the same number of u and u) are relevant among all possible quartic interactions in H 4 because the corresponding fast oscillations exactly cancel out. The following expressions for the leading terms in (16) are obtained, up to fourth order in ε. We draw the reader's attention to the extra term in H 4 (−i sgn(D)η, ξ) as compared to H 4 (η, ξ), which will be further discussed below. If we Taylor expand the dispersion relation ω(k 0 + εD X ) in ε, and retain terms of order up to O(ε 4 ), then the Hamiltonian (16) further reduces to where ω 0 = ω(k 0 ) = gk 0 and Initials 'c.c.' in (21) stand for the complex conjugate of all the preceding terms on the right-hand side of this equation, and Im(·) denotes the imaginary part.
We deduce from (20) that the evolution equation for u is given by which is a Hamiltonian counterpart to Dysthe's equation for deep-water gravity waves [1], associated with the Hamiltonian (21). The prefactor ε −1 in (20) should not be alarming because it suitably cancels out with the prefactor ε in (21). The operator D X is substituted with ∂ X = iD X wherever appropriate in (21) and (22). The cubic NLS equation is recovered if all O(ε 3 ) terms in (22) are neglected. The nonlocal term, which involves the Fourier multiplier in connection with the Hilbert transform, is a characteristic feature of Dysthe's equation and represents effects from the wave-induced mean flow [1,22]. By "mean flow", we refer to the lowest harmonic that is a slowly varying function of space. This nonlocal term naturally arises here from the expansion of H 4 (−i sgn(D)η, ξ) which is itself a direct consequence of the normal form transformation. It takes a simple explicit form, which differs from previous formulations of Dysthe's equation where it is coupled to the solution of an auxiliary Laplace problem for the unperturbed domain [4].

Reconstruction of the Free Surface
Because Equation (22) only solves for the wave envelope, another step is required in order to determine the actual shape of the free surface from knowledge of the wave envelope. In modulation theory, this reconstruction is typically a perturbative calculation in terms of a Stokes' expansion based on the modulational Ansatz. Contributions from up to the third harmonics of the wave spectrum are typically included in this expansion.
Our reconstruction procedure differs from the traditional method in that it inverts the transformations associated with our modulational Ansatz and with the third-order normal form. At any instant t, the original η is recovered by solving the inviscid Burgers' equation for s ∈ (−1, 0] with initial condition where u obeys (22) and η is related to η through (13). The choice of this initial condition is motivated by the special form of wave packets (17) and (19) we are modeling in this asymptotic regime. The final solution of (23) at s = 0 is meant to represent the original physical variable η. From (24), higher-order contributions from lower (i.e., mean flow) and higher harmonics are generated through nonlinear interactions according to (23). The short interval s ∈ (−1, 0] ensures that the smooth and small initial condition (24) cannot yet develop into a shock. As a consequence, the entire solution process including the surface reconstruction is carried out within the Hamiltonian framework. This approach was validated in [20] against direct numerical simulations of the full nonlinear Equations (6) and (7), as applied to the Benjamin-Feir instability of Stokes waves. The agreement was found to be slightly better than that obtained from predictions by the classical Dysthe's equation.

Envelope Equation
In view of comparing with laboratory experiments where time series of the surface elevation are typically recorded by probes at fixed locations along the wave tank, much consideration has also been given to the spatial version of Dysthe's equation as indicated by the numerous studies in the literature [26][27][28]. Such a model is established by switching the role of the Xand t-variables in the previous envelope equation, with X acting as an evolutionary variable for wave propagation along the wave tank.
Because we are unaware of any Hamiltonian formulation for the spatial dynamics associated with the basic Equations (1)-(4), we proceed by applying the standard method to (22). Let v(X, t) = u(X, t)e iω 0 t , so that Then using (26) recursively, all X-derivatives on its right-hand side can be replaced by t-derivatives up to a certain order. This implies and after inserting these expressions back in (26), we obtain which is the spatial version of (22). The coefficient ε −1 in (27) can be removed by reintroducing the original spatial variable x = X/ε, which yields The temporal counterpart to D is defined by D t = i∂ t so that its Fourier symbol is ω. The variable u is recovered by inverting (25). As is common for the spatial Dysthe's equation [26][27][28], there is no ∂ 3 t v term in (28). We note however an additional nonlinear term of the form v 2 ∂ t v and, partly for this reason, the Hamiltonian nature of (22) is not retained by (28). This result is not viewed as undesirable because, as stated earlier, it does not contradict any basic Hamiltonian formulation for the spatial dynamics of water waves. Furthermore, because wave-tank experiments are usually susceptible to dissipative effects [26,34,41], it is not clear whether a Hamiltonian model would be relevant in this context. Despite the absence of a conserved Hamiltonian, the L 2 norm (also called wave action) M = |v| 2 dt , is an invariant for the spatial dynamics of (28). This property can be proved by a direct calculation.

Reconstruction of the Free Surface
When using (28), it is also required to adjust the procedure for surface reconstruction accordingly. For (23), conversion of ∂ x η to ∂ t η is not straightforward given the complicated nature of (6) and (7). We use the approximation with η obeying simple advection at the linear phase speed ω 0 /k 0 . We choose (29) for right-moving solutions to be consistent with the advection part in (22). The new form of Burgers' equation is where For the initial condition (24), we first Taylor expand it as and then convert it to by virtue of (25) and (27). For simplicity, we denote a 0 = a(k 0 ), a 0 = a (k 0 ) and θ = k 0 x − ω 0 t. Parity conditions imply that a(k 0 ) = a(−k 0 ) and a (−k 0 ) = −a (k 0 ), with In this initial condition, the parameter x would be assigned the location of a probe along the wave tank. The symbol Re(·) denotes the real part.
The above expression of η(x, t) s=−1 relies on the scale separation between x and X, thus allowing the Fourier multiplier a −1 (D) to be expanded up to first order in ε. Equation (27) is invoked to approximate the X-derivative of v by t-derivatives. An alternate, perhaps simpler, approach consists in switching the xand t-variables directly in a −1 (D) by virtue of (29), yielding where v solves (28). In general, the latter option should be preferred because it does not involve any expansion of a −1 (D).

Rescaled Variables
We now consider unscaled versions of (28)- (31) where the perturbation parameter ε is absorbed back into the variables. To facilitate the comparison with wave-channel measurements, we change to a coordinate system moving at the linear group velocity via the dimensionless rescaled variables where γ is a scale factor that fits the computational domain into an interval of length 2π and B 0 is a characteristic (real-valued) amplitude for the wave envelope such that according to (17) and (19). This change of variables reintroduces ε in the same way as proposed by Lo and Mei [26], and will also allow us to compare our results with their model, which may be viewed as a spatial form of the classical (non-Hamiltonian) Dysthe's equation. Recall that the wave steepness is defined by ε = k 0 A 0 , hence k 2 0 B 2 0 = ε 2 a 2 0 /2, and the group velocity is given by ω 0 = g/(2ω 0 ) for gravity waves on deep water. To a fixed observer, τ denotes the negative time elapsed and χ represents the distance over which the wave group has advanced (i.e., the fetch). In terms of these new variables, the envelope equation (28) becomes where β = ω 0 a 2 0 /g = 1. The wave envelope is now represented by a complex function V(τ, χ) of τ and χ. For the surface reconstruction, Burgers' equation (30) takes the dimensionless form with initial condition up to order O(ε), where the phase function now reads θ = τ/(εγ) − χ/ε 2 . The alternate choice (31) reduces to after substituting with its definition (18). Note that the Fourier multiplier |D τ | 1/4 acts on all the functions inside the parentheses in (35). For the tests considered in the next section, we found that (34) and (35) produce similar results and we opted to present those given by (35).

Numerical Results
Validation of the Hamiltonian Dysthe's equation (22) was provided in [20] by comparing with time-domain simulations of the full Equations (6) and (7) and of the classical Dysthe's equation. In the present study, we focus on the spatial version (28) and test it against laboratory experiments by Keller [42] on periodic groups and by Su [43] on short-wave packets. These two different sets of experiments were discussed by Lo and Mei [26] and were used to validate their spatial version of the classical Dysthe's equation. We follow their procedure as described in the previous section, and take this opportunity to compare with their model as well. The latter takes the dimensionless form in terms of the same rescaled variables τ and χ as defined earlier. In this case (as is common in modulation theory), the free surface is reconstructed perturbatively via the Stokes' expansion up to the third harmonics, where

Numerical Methods
Equation (32) is discretized in τ by a pseudo-spectral method based on the fast Fourier transform F . The computational domain spans the interval 0 ≤ τ ≤ 2π with periodic boundary conditions and is divided into a regular grid of N collocation points. We denote the discrete Fourier coefficients U(κ) of a 2π-periodic function U(τ) by where κ represents the sideband wavenumbers associated with the variable τ. Integration over χ is carried out in the Fourier space so that the linear terms are solved exactly by the integrating factor technique. The nonlinear terms are integrated in χ by using a fourth-order Runge-Kutta scheme with constant step ∆χ.
The same numerical methods are used to solve Burgers' Equation (33) with initial condition (35) for the surface reconstruction. The reader is referred to [9,20,23,36] for more details on such simulations. In particular, Equation (33) is integrated over s by using the same step size ∆s = ∆χ. We remark that, although the surface reconstruction requires solving an additional nonlinear PDE, this calculation is relatively straightforward and is not necessarily performed at each downstream location χ (only when data on the free surface are needed at the post-processing stage, e.g., for plotting purposes). Furthermore, because this PDE is solved over a relatively short interval s ∈ (−1, 0], the associated cost is insignificant and does not undermine the overall efficiency of this modulational approach. The whole computation (i.e. solving for the wave envelope and determining the surface elevation) can be produced by a single self-contained computer program.
Equations (36) and (37) from Lo and Mei [26] are computed by the same numerical strategy, with exactly the same values for the numerical parameters.

Comparison with Experiments on Periodic Groups
For small amplitudes, the superposition of two progressive waves with slightly different frequencies gives rise to periodic groups whose envelopes move forward at the group velocity without change of shape. Keller [42] conducted laboratory experiments for such groups with maximum incident steepness ε up to 0.23. He found that each group tended to lean forward.
To mimic Keller's input data, Lo and Mei [26] prescribed the incident condition at χ = 0 for (36), so that | A 1 (±1, 0)| = 0.5 for the two parent sidebands (κ = ±1) of the first harmonic A 1 (τ, 0), as in the laboratory experiment. The initial values of the higher sidebands were small and are neglected in (38). The physical parameters are g = 9.81 m s −2 , k 0 = 8.865 m −1 and ε = k 0 A 0 = 0.23. The carrier frequency is estimated by the linear dispersion relation ω 0 = gk 0 . Similarly, we specify the incident condition for (32), so that based on the expression (35) of ζ(τ, χ) s=−1 . Here we deem the first harmonic V 1 (τ, χ) (and its sidebands) to be well approximated in this way because, as suggested by the scaling factor ε 2 of the nonlinear term, Equation (33) is expected to only produce small contributions to wave harmonics of the free surface over the short interval s ∈ (−1, 0]. We also point out that, because the variables A and V do not quite represent the same physical quantity, the numerical coefficients in (38) and (39) are not expected to be identical.
As shown in Figure 1, the computed surface displacements ζ compare favorably with the measurements (which are extracted from Figure 1 in [26]). These are limited sets of data for crests and troughs that were recorded at eight stations within 7.14 m from the wavemaker. Overall, the phase and asymmetry of the wave groups, which tend to lean forward, are well reproduced. Discrepancies are noticeable and may be attributed to experimental errors, physical effects that are not taken into account in the theory (e.g., friction, local wave breaking, undertow current in the wave tank), or uncertainty in the input data. Predictions of ζ based on our simulations of (32) and (36) are also compared in Figure 1. Recall that the free surface is reconstructed via (33) with (35) for our model (32), while it is determined by (37) for the spatial form (36) of the classical Dysthe's equation. The numerical parameters are γ = 0.229, N = 256 and ∆χ = 0.01. As expected, the two numerical solutions are indistinguishable (at least to graphical scale) near the wavemaker. They remain close together, although small discrepancies tend to develop around the dominant crest as the wave group travels down the tank. Model (32) seems to agree slightly better with the experiment, as indicated e.g., at x = 6 m (crest at τ 3.8) and at x = 6.57 m (trough at τ 3.9).   Figure 2a shows the evolution of the parent sidebands (κ = ±1) and higher sidebands κ = ±3 for the first harmonic of the free surface, again comparing predictions from (32) with Keller's measurements (which are extracted from Figure 2 in [26]). The agreement is good overall and is similar to that found for (36) as illustrated in Figure 2b. By contrast, predictions from the cubic NLS equation (by neglecting the higher-order nonlinear terms in Equation (32)) are plotted in Figure 2c and are clearly unsatisfactory. In particular, the corresponding sidebands κ = +1 and +3 behave quite differently from the measurements.  Figure 3 depicts the evolution of the same sidebands for a duration much longer than that of the experiment. Again, the similarity in results between (32) and (36) persists. We can recognize a near-recurrence phenomenon as often reported in such situations [44]. All these results on surface profiles and spectral amplitudes for periodic wave groups are consistent with Lo and Mei's observations (see Figures 1-3 in [26]). Finally, the evolution of the relative error on the rescaled wave action for (32) (U = V) and for (36) (U = A) is illustrated in Figure 4. The integral is computed via the trapezoidal rule over the periodic interval [0, 2π]. The reference value M 0 denotes the incident value of (40) at χ = 0. Overall, M is well conserved in both cases, despite a gradual loss of accuracy that is partly due to accumulation of numerical errors. Model (32) seems to perform slightly better in this regard, displaying a slightly lower error overall. Note that, although filtering was not used for the tests in Figures 1 and 2, it was needed for the tests in Figures 3 and 4 (using an exponential low-pass filter as in [9]) so that the numerical solution can be computed over the longer interval χ ∈ [0, 15] given the relatively large incident steepness ε = 0.23. Use of filtering likely contributes to the loss of accuracy as seen in Figure 4. We emphasize that the same filter together with the same values of numerical parameters were specified in the simulations of (32) and (36), and we point out that Lo and Mei [26] also used filtering.

Comparison with Experiments on Short-Wave Packets
Feir [45] performed laboratory experiments for packets of bell-shaped envelope and of roughly equal duration but different maximum amplitudes. For low wave slope, the envelopes remained symmetrical and flattened with distance. For relatively steep waves however, the envelopes first steepened forward and then split into two groups with the smaller group trailing behind. More extensive experiments in a long tank for initially square envelopes with maximum incident steepness ε = 0.09-0.28 were conducted by Su [43]. Surface displacements were recorded at eight stations within 106.7 m from the wavemaker. Su [43] confirmed Feir's results on the asymmetry of evolution and further revealed that the separated groups have different frequencies.
We consider a case that was investigated by Su [43] and, following Lo and Mei [26], the incident condition for (36) is taken to be where C = 0.5 and j a = 124N/256, j b = 134N/256. As so defined, the edges are slightly smoothed to avoid numerical instability. We specify the incident condition V(τ, 0) for (32) to be of the same form as (41) with the exception that C = 0.495. In this way, the maximum amplitude of ζ(τ, 0) s=0 via (33) coincides with the maximum amplitude of ζ(τ, 0) from (37). Figure 5 shows the computed free surface ζ at four stations along the wave tank for ε = 0.09 and ω 0 /(2π) = 0.96 Hz. The carrier wavenumber is estimated by inverting the linear dispersion relation k 0 = ω 2 0 /g. Wave profiles from our simulations of (32) and (36) are again compared in this figure. Discrepancies are barely noticeable for as far as x = 106.7 m. The numerical parameters are γ = 0.0868, N = 1024, ∆χ = 0.001 and filtering was not needed here. These results closely resemble the computations presented in Figure 5 from [26]. Note the slightly asymmetric envelope between x = 42.7 m and x = 91.5 m, which is also visible in Su's records. The wave period is longer in the lower group in front of the main group. The backward leaning of the main group is associated with nonlinear modification of the dispersion relation, as pointed out in [26]. A more quantitative analysis is provided in Table 1 where predictions of the maximum local steepness k 0 A m based on our simulations of (32) and (36) are compared to Su's measurements at all eight stations (which are reproduced from Table 1 in [26]). The agreement is good overall. The computations tend to overestimate the maximum local steepness as the wave travels down the tank, which is expected considering that dissipation is not accounted for in the theory. Indeed, dissipative effects on wave propagation are likely to be relevant for such a long tank as used in [43]. Compared to (36), predictions from (32) are found to be slightly closer to the experimental data at the three stations farthest downstream.

Conclusions
We have given an overview of the Hamiltonian modulational approach that has been advocated in a series of papers by Craig et al. [21][22][23][24][25] and which leads to a Hamiltonian version of Dysthe's equation for two-dimensional gravity waves on deep water, as recently shown in [20]. We have described the main steps in the derivation of the evolution equation for the complex wave envelope, together with the procedure to reconstruct the free surface from this wave envelope. Both operations are performed within the Hamiltonian framework, which ensures some consistency throughout the entire solution process. A key ingredient in this approach is a Birkhoff normal form transformation that eliminates all non-resonant cubic terms, and yields a simple and elegant expression for the reduced Hamiltonian. As a consequence of this reduction, the surface reconstruction is a non-perturbative calculation but requires solving a nonlinear PDE which, in this problem, takes the generic form of an inviscid Burgers' equation.
To facilitate comparison with laboratory experiments, an alternate spatial form of this Hamiltonian model can be written. This involves switching the role of the spatial and temporal variables in the equations governing the envelope evolution and the surface reconstruction. We have proposed an improved strategy for the latter computation, where a refined choice of the initial condition is specified to solve Burgers' equation. Focusing on the spatial dynamics and following Lo and Mei [26], we have introduced a dimensionless rescaled form of these equations and performed numerical tests against laboratory experiments by Keller [42] for periodic groups and by Su [43] for short-wave packets. Satisfactory agreement is obtained in all these tests. We have also validated our numerical results by comparing with the spatial form of the classical Dysthe's equation. Predictions from our new model and its classical counterpart are found to be similar in all the cases considered.
It would be of interest to extend the present methods and models to the finite-depth case as well as to the three-dimensional setting. Preliminary results in these directions based on the DNO can be found in [46]. These problems pose additional challenges and their study is envisioned for future work.
Author Contributions: Investigation, P.G., A.K. and C.S.; software, P.G. and B.X.; supervision, P.G. and C.S.; writing-review and editing, P.G., A.K., C.S. and B.X. All authors have read and agreed to the published version of the manuscript.