Soliton Turbulence in Approximate and Exact Models for Deep Water Waves

: We investigate and compare soliton turbulence appearing as a result of modulational instability of the homogeneous wave train in three nonlinear models for surface gravity waves: the nonlinear Schrödinger equation, the super compact Zakharov equation, and the fully nonlinear equations written in conformal variables. We show that even at a low level of energy and average wave steepness, the wave dynamics in the nonlinear Schrödinger equation fundamentally differ from the dynamics in more accurate models. We study energy losses of wind waves due to their breaking for large values of total energy in the super compact Zakharov equation and in the exact equations and show that in both models, the wave system loses 50% of energy very slowly, during few days


Introduction
Turbulence in nonlinear continuous media often accompanied by the appearance of localized nonlinear structures.In those cases, when the equations describing the medium have stable soliton solutions, the solitons are natural candidates for the role of such structures, and turbulence can be called soliton turbulence.
The study of soliton turbulence in various nonlinear models is of considerable interest that can be confirmed by numerous works.The problem of a statistical description of a large number of solitons in the framework of the Korteweg-de Vries equation was first considered in [1].The kinetic equation was obtained that describes the distribution function of soliton parameters for the case of a low density of solitons.The turbulence of rarified soliton gas plays the key role in the formation of wave field statistics (see the works [2,3]).With the use of direct numerical simulations methods, the collective behaviour of soliton ensembles was also studied in [4].The general case of dense soliton gas was considered in [5,6].
The Nonlinear Schrödinger Equation (NLS) is a special case of such models and a universal model for wave turbulence studying.However, it is an integrable equation, and turbulence in the framework of the NLS equation has essential peculiarities.The interaction of solitons in this model in some sense can be called simple.Collisions of the NLS solitons are completely "elastic", that is, there are no energy exchanges between them and their basic parameters, amplitudes and velocities, do not change.One of the main parameters that play an important role in the dynamics of solitons pairwise interactions is the relative phase of solitons at the moment of their collision.For example, the maximum amplitude amplification of the solitons in a collision is determined by the synchronization of their phases.Phase synchronization also plays an important role in the formation of extreme amplitude waves or rogue waves.In recent works [7,8] phase synchronization in multisoliton ensembles was studied analytically.Using the inverse scattering transform formalism, the authors of work [9] show that the dense soliton gas is the model of the asymptotic stage of the noise-induced modulation instability in the NLS equation.Rarefied soliton gas was also reproduced experimentally in an optical fiber ring resonator with the initial wavefield taken as arithmetic superposition of individual solitons [10] and in water wave tank [11].
It was really mentioned that the NLS equation is completely integrable.It makes it possible to construct the infinite amount of its periodic solution expressed in terms of Jacobi theta-functions.Some of these solutions described the nonlinear stage of modulation instability of the Stokes waves.The theory of these solutions is far elaborated in [12][13][14].Some of these solutions are beautiful.We plan to study all of them numerically in the closed future.
The soliton turbulence can appear as a result of the modulational instability of the homogeneous wave train.In the pioneer work [15] the evolution of a nonlinear wave train in deep water was studied experimentally and numerically.In particular, it was shown that at later stages of evolution, the wave train exhibits the Fermi-Pasta-Ulam recurrence phenomenon, which is consistent with the results of numerical simulation in the framework of the nonlinear Schrödinger equation.In a later work [16], it was shown that an arbitrary-shaped pulse ultimately decays into a set of solitons, the number of which is determined by the initial conditions.An experimental study of the spatial evolution of unidirectional waves on the water with the formation of nonlinear coherent structures coexisting with smaller stochastic waves was carried out in [17].The properties of extreme water waves developing as a result of modulation instability of initially uniform slightly modulated wave train in the framework of NLS equation and full Euler equations were examined in work [18].The authors identified a maximum possible amplitudes of non-breaking waves in the wide range of initial wave train parameters.
In this work, we investigate and compare soliton turbulence appearing as a result of modulational instability of the homogeneous wave train on the surface of deep water in the framework approximate and exact nonlinear hydrodynamic models.
The set of equations describing a two-dimensional potential flow of an ideal incompressible fluid of infinite depth with a one-dimensional free surface in a gravity field are well known and can be written as: Here x and y are the horizontal and vertical coordinates, t is time, g is the free-fall acceleration, η(x, t) is the shape of the surface, φ(x, y, t) is the hydrodynamic potential inside the fluid.The first equation in (1) is the Laplace equation for the hydrodynamic potential, while the second and third equations are kinematic and dynamic conditions at the fluid surface.As well known, a one-dimensional potential flow of an ideal fluid of infinite depth in the presence of gravity is a Hamiltonian system.As was shown by Zakharov [19], the surface elevation η(x, t) and the velocity potential at the surface ψ(x, t) = φ(x, y, t)| y=η of the fluid are canonically conjugated variables and satisfy the following equations: Here H is the Hamiltonian, i.e., total energy of the fluid: In the assumption of small wave steepness the Hamiltonian can be represented as the infinite series of powers of η and ψ (see [19]).We consider this series up to the fourth order term: Here the operator k means multiplication by |k| in k-space.
In 2011 (see papers [20,21]), Dyachenko and Zakharov assuming that all waves propagate in a single direction applied a canonical transformation to remove all cubic nonlinear terms and to drastically simplify fourth-order terms in the truncated Hamiltonian (4).This transformation is possible due to the unexpected cancellation of non-trivial four-wave interactions in the one-dimensional case [22].In 2017 analyzing the form of the coefficient of four-wave interactions, we found another the most optimal canonical transformation to new normal variables, which allows us to simplify the fourth-order term in the Hamiltonian further (see [23]).The water wave equation in terms of these variables has a remarkably simple form that we named "the super compact equation".The super compact equation includes only two nonlinear terms-wave term and an advection term, that can describe the initial stage of wave-breaking.Due to the specific type of cubic nonlinearity in that equation, it allows introducing an exact envelope for waves without the assumption of narrowness bandwidth (see [24]).This compact Dyachenko-Zakharov envelope equation reduces to the well-known Dysthe equation and further to the NLS equation in the limit of small steepness and assumption of the narrowband spectrum.
In this work, we observe soliton turbulence appearing as a result of modulational instability of the homogeneous wave train in the framework nonlinear hydrodynamic models for free surface waves: • the nonlinear Schrödinger equation (NLS), • the super compact Zakharov equation [23] (SCZ), • the fully nonlinear equations written in conformal variables [25] (RV).
The measure of the nonlinearity of wave at the surface of deep water is their average steepness µ =< ∇η >.We perform the series of numerical simulations in the framework of these models for a very long time (about hundred of thousands of characteristic wave period) varying initial measures of nonlinearity µ 0.04 − 0.08.We set the same initial conditions for all models, which allows us to compare the results and identify the applicability of approximate weakly nonlinear equations.We also take into account energy dissipation due to wave breaking in the super compact Zakharov equation and the fully nonlinear equations using the simple model of dissipation of such waves that we proposed earlier in [24].The results obtained in the framework of super compact Zakharov equation and the RV equations show that for moderate values of the average slope of the wave for a hundred thousand characteristic periods of the wave, the sea loses about 40-50% of its energy.

Nonlinear Models for Deep Water Gravity Waves
In this section, we very briefly derive three nonlinear hydrodynamic models for free surface waves: the super compact Zakharov equation based on the truncated Hamiltonian (4), the nonlinear Schrödinger equation and the fully nonlinear equations written in conformal variables.All the details of the derivation can be found in [23,26].

Super Compact Zakharov Equation and the Nonlinear Schrödinger Equation
We use canonical transformation from the pair of classical physical Hamiltonian variables η(x, t) and ψ(x, t) to normal complex variable c(x, t).This canonical transformation can be written as power series of c k up to the third order: ( and described in detail in [23].As a result of the transformation the Hamiltonian takes the following form in x-space: Here the operator Vk is multiplication by ω k k in Fourier space.And the equation of motion or the super compact Zakharov (SCZ) equation then is the following: Here the operator ωk is Fourier multiplier by the gk.The operator ∂ + x in the Fourier space is ikθ(k), where θ(k) is the Heaviside step function.And we also introduce the advection velocity U : The physical variables η and ψ can be recovered by the canonical transformation (5).An important property of this canonical transformation is that the linear and second order terms can be written in x -space, in a compact way: Here Ĥ-is the Hilbert transformation with eigenvalue isign(k) and operators kα act in Fourier space as multiplication by |k| α .Now if we suppose that Fourier spectrum of waves has a maximum at k = k 0 we can introduce envelope function where ω k 0 = gk 0 is the corresponding linear frequency.Using the Equation (10), which is actually canonical transformation, one can easily derive the exact equation for the envelope C that we named "Dyachenko-Zakharov envelope equation" [24].
The Dyachenko-Zakharov envelope equation written in the reference frame moving with the group velocity V 0 = ω k 0 2k 0 has the following form: The Dyachenko-Zakharov envelope Equation ( 11) is Hamiltonian, and the Hamiltonian is where the operator Vk has the following form in k-space: The Dyachenko-Zakharov envelope equation was derived without any assumptions about the spectral width of the wave packet.Moreover, the Equation ( 11) is the exact envelope form of the Equation ( 7) and has the same range of applicability.
One can extract the nonlinear Schrödinger (NLS) equation in the assumption of narrowband spectrum: ∂C ∂t The Hamiltonian has the following form:

Fully Nonlinear Conformal Equations
Following the work [26] we perform conformal mapping of the free surface liquid domain z = x + iy confined by a free boundary y = η(x, t) onto a half-plane of the new complex variable w = u + iv having a fixed boundary v 0 --see Figures 1 and 2. This transformation can be written in the following form: Then we use the following variables suggested by Dyachenko [25]: We define the functions U and B using the Hilbert operator and the projection operator P = 1 2 (1 + i H) as: Here and below the asterisks stands for complex conjugation.In the new variables ( 17) the fully nonlinear (RV) equations have the following form: with the boundary conditions:

Numerical Simulation of Soliton Turbulence. Initial Conditions
In this section we explain how we generate initial conditions for numerical experiments on soliton turbulence in the following nonlinear models for water waves: • the NLS Equation ( 14), • the SCZ Equation ( 7), • the RV Equation (18).
The pseudo-spectral Fourier method was used for the SCZ, and RV equations solving the fourth-order Runga-Kutta method was applied to calculate time evolution.The time integration step was chosen to ensure that the values of the integrals of motion-energy, momentum, and mass-are kept accurate to 10 decimal places.The FFTW3 library [27] was used for the fast Fourier transform procedure.The multiplication of grid functions was carried out in x-space and to calculate derivatives and nonlocal terms, the direct and inverse Fourier transform were used.The source codes for numerical simulation of soliton turbulence in the frame of SCZ and RV equations can be found by the following links: http://kachulin.itp.ac.ru/sceq/ and http://kachulin.itp.ac.ru/rveq/.
We perform simulations in periodic domain of length L = 10 km with gravity acceleration g = 9.81 m/s 2 .The initial conditions for numerical experiments in all three models are identical.They represent a homogeneous wave train with two harmonic perturbations having a random phase: Here, φ 0 and φ ±p are the random, but identical for all models, phases of monochromatic wave and the perturbations correspondingly.
We set the initial conditions in the SCZ equation in the form (20) and in the NLS equation directly using the transformation from c(x, t = 0) to the envelope function-see Equation (10).
Then we generate this initial conditions in conformal variables R(u, t = 0) and V(u, t = 0) for the fully nonlinear conformal Equation (18) with the following procedure.
First, we recover the surface profile η(x, t = 0) and the velocity potential on it ψ(x, t = 0) employing the canonical transformation (9).The conformal mapping described in the Section 2.2 satisfies the following equation: which can be solved by the iterative procedure: Here the index n stands for the step of the iteration.After the conformal mapping is found by the use of ( 22), we compute the function R only using its definition (17).Now in order to find the function V by its definition (17) we need to know the potential Φ(w) in conformal variables.The potential Ψ(u) in conformal variables can be computed using the Fourier transform of the known spectrum of the potential on the free surface ψ k as: where N is the number of grid points.
The carrier wave number k 0 in all numerical experiments is equal to 2π 100 , so that the characteristic wavelength is λ 0 = 100 m and the time period of the carrier wave is T 0 = 2π ω k 0 ≈ 8 s.The perturbations amplitude is 20 times less that the amplitude of the carrier wave δc 0 = c 0 20 , while c 0 is chosen so that in different experiments the average steepness take the values µ = 0.04, 0.06, 0.08.In addition for each average steepness µ we vary the dimensionless relative wave number of perturbation kp = L 2π k p so that: • kp = 1.In this case the harmonics adjacent to the carrier are perturbed.
• kp ∼the wave number of the maximum increment of modulation instability.For the average steepnesses µ = 0.04, 0.06 and 0.08 we set kp = 10, 15 and 20 correspondingly.

Numerical Simulations of Soliton Turbulence
In this section, we focus on the case of low wave energies performing numerical simulations in the three different nonlinear models for water waves.Figure 3 shows the dependence of total energy on time in these models.These dependencies were obtained in simulation with average wave steepness µ = 0.04 and the dimensionless perturbation wavenumber kp = 10 corresponding to the maximum increment of modulation instability.One can see from Figure 3 that in the NLS and the SCZ models, the total energy conserves throughout the entire simulation time.Meanwhile, in the fully nonlinear model, the formation of pre-breaking waves leads to the switching on of the damping at the time moment t damp ∼ 45 h-see Figure 3.We discuss the details of the damping model in the next section and now compare surface dynamics in different models until the time moment t damp .The modulation instability results in the formation of spatially periodic structures on the surface of the water-breathers.Besides, the breathers demonstrate an approximate periodic behavior in time-see the snapshots of the surface evolution during one such time period in Figure 4. We observe about 10 periods of the breather's evolution in each of three nonlinear models for water waves after what, and the free surface dynamics become more complex.However, for substantial evolution time, we report significant differences in the evolution of the free surface when the NLS equation compared to the more accurate models-see the discussion below.
Note that the simulations of the nonlinear Schrödinger equation and the super compact Zakharov equation were carried out in a reference frame moving with group velocity V 0 = 1 2 g k 0 , while simulations of the exact equations in the laboratory frame of reference.
Figure 5 shows how the NLS equation dynamics become more complex at large simulation times, t ≈ 20 h, exhibiting several different breather structures, though approximate periodic recurrences of the weakly modulated plane wave can be observed.At the large time we observe an almost complete reassemble of the initial condition comprised of 10 breathers structures on the fluid surface-see Figure 6.One of the key results of this work is that in both the super compact Zakharov and the fully nonlinear models the surface waves split into a set of interacting soliton structures at about the same time (t ≈ 20 h), after what no recurrence to the initial weakly modulated plane wave can be observed-see Figures 7 and 8  Dependence of maximum value of the surface amplitude on time turned out to be a crucial characteristic allowing to identify a fundamental difference between wavefield evolution in the three considered models-see Figure 9.We find that this function is almost periodic in time and does not exceed 3 m for the NLS equation model, meanwhile in more accurate models, it loses periodicity and starts to reach 4-6 m after t ≈ 20 h of simulations-see again Figure 9. Then we compare the statistical behavior of the studied models by computing the probability density functions of the surface amplitude PDF(η) averaged in time over the whole simulation-see Figure 10.We find that the NLS equation exhibits a low probability of high amplitude wave formation in comparison to the more accurate models.For example, the probability of formation of the 3-m waves in the NLS equation is from 10 to 100 times less than in the last ones.In contrast, there is a non-zero probability of formation of waves, exceed 3 m in the SCZ and the RV equations, while in the NLS equation, such events have not been detected.It is also important to note that the function PDF(η) corresponding to the NLS equation model has kinks at η ∼ 1 m due to the quasiperiodicity of the free surface dynamics discussed above.Now we study soliton turbulence in the different models at µ = 0.04 and kp = 1, i.e., in the case when the harmonics adjacent to the carrier are perturbed.Figure 11 shows the observed total wave energy dependencies on time.One can see that the total energy conserves throughout the whole simulation for the NLS equation and the SCZ equation.It was necessary to use the damping for the fully nonlinear model twice-at t ≈ 11.5 h and t ≈ 17 h.
Snapshots of surface dynamics in the three considered models are presented in Figure 12.At the initial stage in each model, the modulation instability leads to the appearance of a single wave of high amplitude on the free surface-see panels (a) in Figure 12.At this time, the spectra broaden significantly-see Figure 13, which shows the dependence of the spectra of the envelope function C(x) of the NLS equation and function c(x) of the SCZ equation.Then this wave splits into a set of coherent wave structures-see panels (b) in Figure 12.In the NLS model, these structures resemble breathers having similar velocities, approximately equal to the group velocity V 0 = 1 2 g k 0 corresponding to the characteristic wavelength 2π k 0 .In the more accurate nonlinear models, the arising objects mostly resemble localized soliton structures having different velocities and interacting with each other.As a result of their collision, large-amplitude waves can appear on the surface.In the RV equations, such waves reach the pre-breaking stage at t ≈ 11.5 h and t ≈ 17 h.Furthermore, therefore at these moments, the damping was switched on.The entire free surface dynamics, computed in different models, can be found by the following links: 1.
Similar to the case studied in the first part of this section, now we again observe that the surface dynamics demonstrate almost recurrence to the initial weakly modulated plane wave for the NLS equation.While the more accurate models exhibit fundamentally different complex dynamics of soliton turbulence-compare panel (d) in Figure 12I with panels (d) in Figure 12II,III.Found differences between the dynamics of the nonlinear Schrödinger equation and the more accurate models at µ = 0.04 and kp = 1 are also confirmed by the measured dependence of the maximum surface amplitude on time shown in Figure 14 and the PDF(η) shown in Figure 15.Not that the registered probability of large-amplitude wave formation in the fully nonlinear model is less than in the super compact Zakharov equation (see Figure 15) due to the fact that the PDFs are measured by averaging over the entire simulation time.Meanwhile, in the case of the fully nonlinear model, the total energy decreased due to the use of damping, and therefore large-amplitude waves appear less often at the end of the simulation.Also, Figure 14 shows that at large times the maximum amplitude of the surface in the super compact equation is almost always higher than in the full model.We conclude that even at a low level of energy and average wave steepness, the wave dynamics in the nonlinear Schrödinger equation fundamentally differs from the dynamics in the more accurate models.Meanwhile, the super compact Zakharov equation model demonstrates similar results when compared to the fully nonlinear equations.This is since in the soliton turbulence observed by us, even at low energy levels, high and steep waves appear on the water surface.The appearance of such waves is accompanied by a significant broadening of the spectrum-see again Figure 13.The NLS equation is a weakly nonlinear model, and, unlike the more accurate model, the SCZ equation, it was obtained using the additional assumption of the narrowband spectrum.Such events go beyond the applicability of the nonlinear Schrödinger equation, and therefore the surface dynamic is described incorrectly by this model.
In the next section, we compare numerical simulations of soliton turbulence for large values of total energy and averaged wave steepness only for the super compact Zakharov equation and the exact equations.We will pay special attention to the energy losses of wind waves due to their breaking.

Energy Losses Due to Wave Breaking at High Energy Levels
Recall that, modulation instability leads to the formation of high amplitude and/or large steepness waves that can break.Therefore we add damping providing dissipation of energy for pre-breaking waves to the super compact Zakharov equation and the exact equations in conformal variables.
Previously we studied in detail the initial stage of wave breaking and proposed model for damping, which removes only the pre-breaking waves, leaving the other part of the surface unchanged [24].The idea of our damping model is the following.Fourier spectrum exhibits exponential behavior at large wavenumbers shortly before the wave breaking.The slope of the exponent determines the distance from the real axis to the nearest singularity in the lower half-plane x + iy.This singularity moves closer and closer to the real axis at the pre-breaking stage.In order to prevent the wave breaking, we move the singularity closest to the x-axis away from it.More precisely, we multiply the Fourier spectrum by a special attenuation function that decreases exponentially at large wavenumbers k: We use breaking criterion based on advection velocity (8) in numerical simulation in the framework of the SCZ equation.Bjorkavag and Kalisch in the paper [28] showed that using a similar criterion for Boussinesq models gives good qualitative agreement with the experiments.If the value of the advection velocity exceeds half the group velocity, we stop the calculation and remove the pre-breaking wave.We "fix" the wave spectrum using the following damping function: and the coefficients were chosen as D = 400 and α = 0.75.Similarly, we took into account damping for the RV equations.Practically, in this case, the needed number of points N of the computational grid is about 2-4 times more than required for the SCZ equation for the same initial conditions due to specific properties of conformal mapping.Therefore we vary N in the numerical experiments with the exact model.Initially, the calculation was carried out on a grid with N = 16,384 points, so that the number of negative harmonics in the spectrum of R and V was 8192.When the spectrum of R (and V) became wider, and the average amplitude of 100 highest harmonics R av,N exceeded the threshold value R th equals to 10 times of the roundoff error, the number of computational nodes was doubled to Ñ = 32768.
We define the average amplitude R av,N as If the wave does not break and the width of the spectrum decreases, the number of computational points decreases to its initial value N =16,384.If the spectrum continues its expansion and the average amplitude of the 100 higher harmonics in the R spectra R av, Ñ with the number of the grid points Ñ =32,798 exceeds the roundoff error by a factor of 10 10 , then the simulation stops and we remove the pre-breaking wave.
Modulation instability at high "levels of nonlinearity" (µ = 0.06 and µ = 0.08 ) leads to the formation of an pre-breaking wave on the free surface.Figures 16 and 17 show the free surface profile η(x) and the spectrum of R immediately before and after the damping for numerical simulation in the framework of RV equations.One can see the spectrum dynamics in more details shortly before the damping right after that by following link: http://kachulin.itp.ac.ru/STmovies/spectrum-wb.avi, while corresponding dynamics of the free surface profile by http://kachulin.itp.ac.ru/STmovies/surface-wb.avi and http: //kachulin.itp.ac.ru/STmovies/surface-wb-zoom.avi.
The free surface profiles shown in Figure 16 demonstrate that damping acts only at the wave breaking point and does not affect the rest part of the water surface.In numerical experiments with initial parameters µ = 0.06, kp = 1 and µ = 0.08, kp = 1 the pre-breaking waves split into a set of interacting solitons in both models.The surface dynamics in these cases is similar to that which was discussed in the previous section.The difference is that at high energy levels, collisions of solitons more likely lead to the formation of large-amplitude waves.
Figures 18 and 19 show the total energy dependence on time in both models for the cases µ = 0.06, kp = 1 and µ = 0.08, kp = 1 correspondingly.
As one can see from Figures 18 and 19, the energy loss due to wave breaking is higher for the exact model in comparison with the approximate Zakharov equation.The differences in energy losses observed in two models considered can be explained by two following reasons.First, the super compact Zakharov equation is the approximate model corresponding to the Hamiltonian expanded in powers of η and ψ up to the fourth-order term.The expansion was done under the assumption of low wave steepness.Therefore, the dynamics of pre-breaking waves with large values of steepness will differ from the results obtained by using the exact RV equations.Second, the function c(x) in the SCZ equation and the function of R(u) in the exact equations differ in sense.c(x) is a complex normal Hamiltonian variable and related to the natural Hamiltonian variables η(x) and ψ(x) by canonical transformation (5), while the function of R(u) is related to conformal mapping by definition (17).Thus, although the fundamental idea of energy loss accounting is the same in these two models, it seems impossible to achieve complete agreement on damping procedure.However, in both models, the wave systems lose only about 40-50% of their energy for two days.

Conclusions
In this work, we investigate and compare soliton turbulence appearing as a result of modulational instability of the homogeneous wave train in three nonlinear models for surface gravity waves: the nonlinear Schrödinger equation, the super compact Zakharov equation, and the fully nonlinear equations written in conformal variables.

Figure 3 .
Figure 3. Dependence of the total energy E tot on simulation time for numerical experiments with µ = 0.04 and kp = 10 in different models.The black dots 1 corresponds to the NLS model, the blue solid line 2-to the SCZ equation and the red solid line 3 to the RV equations.

Figure 4 .Figure 5 .Figure 6 .
Figure 4. Snapshots of the free surface showing formation and evolution of periodic breathers in different nonlinear models.(I) corresponds to the NLS equation, figure (II) to the SCZ equation, and figure (III) to the RV equations.In each subfigure the panels (a) show the initial moment of time t = 0, while the panels (b), (c) and (d) correspond to the moments of time t ≈ 0.7 h, t ≈ 1.2 h and t ≈ 1.8 h. .

Figure 7 .Figure 8 .
Figure 7. Simulations in the framework of the SCZ model at large times.Different panels show snapshots of the free surface at: (a) t ≈ 20.5 h (b) t ≈ 22.2 h (c) t ≈ 38.3 h and (d) t ≈ 53.8 h.

Figure 9 .
Figure 9. Dependence of maximum surface amplitude on time in three different nonlinear models at µ = 0.04 and kp = 10.Black solid curve 1 corresponds to the NLS model, blue solid curve 2-to the SCZ equation and red solid curve 3-to the RV equations.

Figure 10 .
Figure 10.Probability density functions of surface amplitude PDF(η) in logarithmic scale for three three different nonlinear models at µ = 0.04 and kp = 10.Black solid curve 1 corresponds to the NLSe model, blue solid curve 2-to the compact Dyachenko-Zakharov envelope equation and red solid curve 3-to the fully nonlinear equations.

Figure 11 .
Figure 11.Dependence of the total energy E tot on simulation time for numerical experiments with µ = 0.04 and kp = 1 in different models.The black solid dots 1 corresponds to the NLS model, the blue solid line 2-to the SCZ equation and the red solid line 3 to the RV equations.

Figure 12 .
Figure 12.Snapshots of the free surface showing wave field evolution in different models at large times for µ = 0.04 and kp = 1. Figure (I) corresponds to the NLS equation, figure (II) to the SCZ equation, and figure (III) to the RV equations.In each subfigure the panels (a) show the moment of time t ≈ 2.8 h, while the panels (b), (c) and (d) correspond to the moments of time t ≈ 6.2 h, t ≈ 34.6 h and t ≈ 39 h.

Figure 13 .
Figure 13.Spectra of |C k | and |c k | in logarithmic scale at the moment of appearance of a single high amplitude wave on the free surface t ≈ 2.8 h.Black solid curve 1 corresponds to the NLS model, blue solid curve 2-to the SCZ equation.

Figure 14 .
Figure 14.Dependence of maximum surface amplitude on time in three different nonlinear models at µ = 0.04 and kp = 1.Black solid curve 1 corresponds to the NLS model, blue solid curve 2-to the SCZ equation and red solid curve 3-to the RV equations.

Figure 15 .
Figure 15.Probability density functions of surface amplitude PDF(η) in logarithmic scale for three different nonlinear models at µ = 0.04 and kp = 1.Black solid curve 1 corresponds to the NLSe model, blue solid curve 2-to the compact Dyachenko-Zakharov envelope equation and red solid curve 3-to the fully nonlinear equations.

Figure 17 .
Figure 17.Fourier spectrum of R in logarithmic scale before (solid purple curve) and after (solid green curve) the damping.

Figure 18 .Figure 19 .
Figure 18.Dependence of the total energy E tot on simulation time for numerical experiments with µ = 0.06 and kp = 1 in different models.The blue solid line 1 corresponds to the SCZ equation and the red solid line 3 to the RV equations.