# Nonlinear evolution equations for broader bandwidth wave packets in crossing sea states.

1. IntroductionThe study of evolution of weakly nonlinear surface gravity waves in crossing sea states has attracted considerable interest in recent years. The reason is that the freak waves (freak waves are extreme oceanic waves, also known as rogue waves or giant waves or killer waves) are frequently observed in crossing sea states, characterized by two different wave-train systems. The extraordinarily large amplitude freak waves generated in crossing sea states may even lead to ship accidents (see [1, 2]). Several mechanisms (see [1, 3, 4])-dispersion enhancement of transient wave groups, effect of variable bathymetry, wave-current interaction, nonlinear modulational instability, statistical approach with the assumption of Gaussian wave field, and so forth--have been suggested for the possible formation of such killer waves. The problem of weakly nonlinear interaction between two wave systems propagating in two different directions has been studied by some authors ([5-7]) as a possible mechanism resulting in freak wave generation. Onorato et al. [5] derived two coupled nonlinear Schrodinger equations that describe the nonlinear interaction of two weakly nonlinear, obliquely propagating wave packets. By the use of these equations they obtained an expression for the growth rate of instability for unidirectional perturbation. Shukla et al. [6] extended their result by plotting instability growth rate diagrams for bidirectional perturbation. In both papers, stability analysis has been performed with two coupled evolution equations, correct up to third order in wave steepness. In Onorato et al. [5] and also in Shukla et al. [6], it is concluded that freak waves canbeformeddue to weaklynonlinearinteraction in crossing sea states. But it is well known that an evolution equation, which is correct up to fourth order in wave steepness, gives a more appropriate description of the dynamical behavior of waves with sufficiently small wave steepness. The reason is that, at fourth order, the wave-induced mean flow terms appear and these terms modify considerably the growth rate of instability ([8]).

Keeping this in view, Gramstad and Trulsen [9] have derived two coupled modified nonlinear Schroodinger (hereafter referred to as CMNLS) equations that describe the evolution of two two-dimensional narrow band wave systems with different directions of propagation. The corresponding coupled evolution equations for broader bandwidth wave packets in crossing sea state condition are derived in the present paper. The importance of the assumption of broader bandwidth wave packets lies in the fact that, in reality, the bandwidth of the ocean wave spectra often exceeds the range of the narrow bandwidth condition. The evolution equation for a single broader bandwidth wave packet was first derived by Trulsen and Dysthe [10]. The order of the bandwidth, for narrowband wave packet, is taken as [epsilon], [epsilon] being a measure of weakness of nonlinearity. Trulsen and Dysthe [10] chose the order of bandwidth to be of O([[epsilon].sup.1/2]) and derived the evolution equation, correct up to O([[epsilon].sup.7/2]) terms. The instability regions of a uniform Stokes wave in the perturbed wavenumber plane obtained from this equation are found to be in better agreement with the exact results of McLean et al. [11]. An equation, which is one order higher than the equation derived by Trulsen and Dysthe [10], was later derived by Debsarmaand Das [12], which is correct up to O([[epsilon].sup.4]) terms. In this paper, we have presented also the fourth order nonlinear evolution equations in crossing sea states for narrow band wave packets. These equations are in a form different from those derived by Gramstad and Trulsen [9], in the sense that the two equations do not involve the wave-induced mean flow velocity potential explicitly. The two evolution equations of Gramstad and Trulsen [9] can however be rewritten as two equations that do not involve mean flow potential terms explicitly by substituting mean flow potential terms [[phi].sub.a] and [[phi].sub.b] as given in (11) of Gramstad and Trulsen [9]. The evolution equations of Gramstad and Trulsen [9] also have certain advantages because of their Hamiltonian structure. It is well known that the Hamiltonian structure provides a systematic account of symmetries inherent to the problem and of corresponding conservation laws (see [13]).

Using the two coupled nonlinear envelope equations for both narrow and broader bandwidth wave packets, we have investigated the stability of two obliquely propagating uniform wave trains. It is observed that the growth rate of instability in crossing sea states is much higher than that for a single wave packet. Expectedly, it is found that the growth rate of instability of one wave train becomes higher with the increase of the amplitude of the second wave train. We have also plotted figures showing a comparison between the growth rate of instability for narrow and broader bandwidth wave packets. It is also observed that the growth rate of instability becomes higher when the angle of interaction between two wave packets becomes smaller. Shrinkage is observed in the unstable region in the perturbed wavenumber plane when crossing sea states with broader bandwidth wave packets is considered.

2. Derivation of Evolution Equations

Taking xy-plane as the undisturbed free water surface, we consider propagation of two wave packets centered at wavenumbers [[??].sub.a] = (k, l) and [[??].sub.b] = (k, -l), respectively. Both waves have the same linear frequency [[omega].sub.c] = [square root of g[k.sub.c]], where [k.sub.c] = [absolute value of [[??].sub.a]] = [absolute value of [[??].sub.b]] and g is the acceleration due to gravity. The two wave packets are meeting at an angle 2[theta]. To derive envelope equations we have followed the standard methodology that starts with Zakharov [14] integral equation in the spatial dimensions as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where B([??], t) is related to the free surface elevation [zeta]([??], t) by the following equation:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

Here c.c. denotes complex conjugate and [B.sup.*] = c.c. of B. The kernel function [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] was first given by Zakharov [14] and, later on, its expression that corresponds exactly to the Hamiltonian was derived by Krasitskii [15]. Because of the presence of the factor [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] in (1) a nonzero contribution to the integral in (1) with k = [[??].sub.a] is obtained for the three cases; namely, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. So, for the derivation of the envelope equation of first wave packet having carrier wavenumber vector [[??].sub.a] we take [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Clearly, anonzero contribution to the integral in (1) is obtained when [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Similarly, for the derivation of the second envelope equation we take [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. It is assumed that the spectral widths [absolute value of [??]]/[k.sub.c], [absolute value of [[??].sub.j]]/[k.sub.c] are of order O(e) for narrow band wave packet and of order O([[epsilon].sup.1/2]) for broad band wave packet, e being the order of smallness of wave steepness. The process of derivation of envelope equations adapted here is similar to the process of derivation of evolution equations for a surface gravity wave packet in the presence of another wave packet described in Debsarma and Das [16].

Finally, the envelope equations are written in terms of a1 and [a.sub.2] which are the complex amplitudes of the two wave envelopes given by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

Hereafter, the equations are written using the following dimensionless quantities:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

and finally (~) notation has been dropped out.

The two CMNLS equations for broader bandwidth wave packets, correct up to order O([[epsilon].sup.7/2]) in wave steepness, are the following:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (6)

Here [H.sub.x] and [H.sub.y] represent, respectively, x and y components of two-dimensional Hilbert transform operator defined by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (7)

The coefficients [[beta].sub.i]'s, [[lambda].sub.i]'s, and [[mu].sub.i]'s are given in Appendix A. When the terms within third bracket on left-hand side of (6) are omitted, the resulting equations become CMNLS equations for narrow band wave packets, correct up to fourth order in wave steepness. The coefficients [[beta].sub.1], [[beta].sub.2], [[beta].sub.3], [[beta].sub.4], [[beta].sub.5], [[lambda].sub.1], and [[mu].sub.1] are in agreement with the corresponding coefficients in Shukla et al. [6].

Setting [a.sub.2] = 0 in the envelope equation (4) we can recover the single envelope equation for a single surface gravity wave packet with carrier wavenumber [[??].sub.a].

3. Stability Analysis

Two obliquely propagating Stokes wave trains as obtained from envelope equations (6) are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (8)

where [a.sub.10] and [a.sub.20] are real constants.

In order to study the effect of fourth order nonlinearity on the stability properties of the above two interacting wave trains, we introduce the following infinitesimal perturbation:

[a.sub.1] = [a.sup.(0).sub.1] (1 + [a'.sub.1]), [a.sub.2] = [a.sup.(0).sub.2] (1 + [a'.sub.2]). (9)

We now substitute (9) in the envelope equations (6) and linearize with respect to perturbed quantities. Assuming [a'.sub.1] = [a'.sub.1r] + [ia'.sub.1i] and [a'.sub.2] = [a'.sub.2r] + [ia'.sub.2i], where [a'.sub.1r], [a'.sub.1i], [a'.sub.2r], and [a'.sub.2i] are real, and then separating into real and imaginary parts we get four coupled equations in [a'.sub.1r], [a'.sub.1i], [a'.sub.2r], and [a'.sub.2i]. Then taking Fourier transform of these four equations with respect to space variables x and y defined according to

[??] = 1/2[pi] [integral] [[integral].sup.[infinity].sub.-[infinity]] f(x, y) [e.sup.-i(Lx+My)] dx dy (10)

and assuming time dependence of the perturbed quantities [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] to be of the form exp (-i[OMEGA]t), we obtain finally the following system of four equations:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (11)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (12)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (13)

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (14)

The condition of nontrivial solution of the above system gives the following biquadratic dispersion relation:

[[([OMEGA] - [W.sub.+]).sup.2] - [S.sub.+]][[([OMEGA] - [W.sub.-]).sup.2] - [S.sub.-]] = P([[OMEGA].sub.-]) ([OMEGA] - [W.sub.-]) - [Q.sub.+] ([OMEGA] - [W.sub.+]) - [Q.sub.-]([OMEGA] - [W.sub.-]) + R. (15)

The coefficients [A.sub.[+ or -]], [B.sub.[+ or -]], [C.sub.[+ or -]], D, [E.sub.[+ or -]], [F.sub.[+ or -]], P, [Q.sub.[+ or -]], R, [S.sub.[+ or -]], and [W.sub.[+ or -]] are given in Appendix B. The subscripts ([+ or -]) are used just for convenience of writing the coefficients.

We have plotted growth rate of instability [G.sub.r] = Im([OMEGA]) both for unidirectional and bidirectional perturbations. In Figures 1-4, growth rate of instability [G.sub.r] is plotted against L, L being perturbation wavenumber along x-direction. In Figure 1, it is observed that the growth rate of instability for two obliquely interacting wave packets is much higher than that for a single wave packet. This is true at third order and also at fourth order. Figure 2 shows that there is an increase in the growth rate of instability of one wave packet when the amplitude of the other wave packet increases. Figure 3 shows a comparison between the growth rate of instability for narrowband wave packets and that of broader bandwidth wave packets. From Figure 3 it is revealed that the growth rate of instability in crossing sea states is much higher than that for a single wave packet having broader bandwidth. In Figure 4 we can observe the variation in [G.sub.r] with respect to [theta] for broader bandwidth wave packets when [a.sub.10] = 0.18 and [a.sub.20] = 0.15. The figure shows that the growth rate of instability is higher when the angle of interaction between the two wave packets is smaller.

Solving the dispersion relation (15) numerically, we have plotted instability growth rate curves in perturbed wavenumber plane. In Figures 5, 6, 7, and 8, we have shown contour plots for instability growth rate for bidirectional perturbations. In Figures 5-7, amplitudes of the two uniform wave trains are considered equal, while in Figure 8, the amplitudes are different. It is known from the investigations of Trulsen and Dysthe [10] and Debsarma and Das [12, 17] that the stability analysis made from higher order evolution equation for broader bandwidth wave packets gives more accurate results. Comparing the upper row figures with the corresponding lower row figures in Figures 5-8, we see that there is little shrinkage in the instability region for broader bandwidth wave packets.

4. Conclusion

Keeping in view the importance of higher order evolution equations valid for broad band wave packets, we have derived here two coupled nonlinear evolution equations for broader bandwidth surface gravity wave packets under the crossing sea states assumption. The evolution equations derived here in crossing sea states can be considered as an extension of the evolution equation for a single broader bandwidth wave packet derived by Trulsen and Dysthe [10]. Using the evolution equations derived here we have also carried out stability analysis of two obliquely interacting uniform wave trains. Results obtained here are in agreement with the corresponding results of Shukla et al. [6] at third order with narrow band wave assumption. Modifications in the instability growth rates with the assumption of broader bandwidth wave packets are shown graphically in this paper. It is found that the growth rate of instability in a situation of crossing sea states is much higher than that for a single wave packet even for broader bandwidth wave packets. So, generation of freak waves due to nonlinear interaction of broader bandwidth wave packets in crossing sea states is confirmed. It is also found that the growth rate of instability is higher when the angle of interaction between the two wave packets is smaller. Comparing the stable-unstable regions for broader bandwidth wave packets with those for narrow bandwidth wave packets, we see that there is little shrinkage in the instability region for broader bandwidth wave packets.

Appendices

A. Coefficients of the Evolution Equations

Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (A.1)

B. Coefficients of the Dispersion Relation (15)

Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (B.1)

http://dx.doi.org/10.1155/2014/597895

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

[1] C. Kharif and E. Pelinovsky, "Physical mechanisms of the rogue wave phenomenon," European Journal of Mechanics B: Fluids, vol. 22, no. 6, pp. 603-634, 2003.

[2] A. Toffoli, J. Monbaliu, J. M. Lefevre, and E. Bitner-Gregersen, "Dangerous sea-states for marine operations," in Proceedings of the 14th International Offshore and Polar Engineering Conference (ISOPE '04), pp. 85-92, Toulon, France, May 2004.

[3] K. B. Dysthe, "Modelling a "rogue wave"--speculations or a realistic possibility ?" in Proceedings of the Rogue Waves, M. Olagon and G. Athanassoulis, Eds., pp. 255-264, Ifremer, France, 2002.

[4] B. S. White and B. Fornberg, "On the chance of freak waves at sea," Journal of Fluid Mechanics, vol. 355, pp. 113-138, 1998.

[5] M. Onorato, A. R. Osborne, and M. Serio, "Modulational instability in crossing sea states: a possible mechanism for the formation of freak waves," Physical Review Letters, vol. 96, no. 1, Article ID 014503, 2006.

[6] P. K. Shukla, I. Kourakis, B. Eliasson, M. Marklund, and L. Stenflo, "Instability and evolution of nonlinearly interacting water waves," Physical Review Letters, vol. 97, no. 9, Article ID 094501, 2006.

[7] F. E. Laine-Pearson, "Instability growth rates of crossing sea states," Physical Review E, vol. 81, Article ID 036316, pp. 1-7, 2010.

[8] K. B. Dysthe, "Note on a modification to the nonlinear Schrodinger equation for application to deep water waves," Proceedings of the Royal Society of London A, vol. 369, pp. 105-114, 1979.

[9] O. Gramstad and K. Trulsen, "Fourth-order coupled nonlinear Schrodinger equations for gravity waves on deep water," Physics of Fluids, vol. 23, no. 6, Article ID 062102, 2011.

[10] K. Trulsen and K. B. Dysthe, "A modified nonlinear Schrodinger equation for broader bandwidth gravity waves on deep water," Wave Motion, vol. 24, no. 3, pp. 281-289, 1996.

[11] J. W. McLean, Y. C. Ma, D. U. Martin, P. G. Saffman, and H. C. Yuen, "Three-dimensional instability of finite-amplitude water waves," Physical Review Letters, vol. 46, no. 13, pp. 817-820, 1981.

[12] S. Debsarma and K. P. Das, "A higher-order nonlinear evolution equation for broader bandwidth gravity waves in deep water," Physics of Fluids, vol. 17, no. 10, Article ID 104101, 2005.

[13] T. B. Benjamin and P. J. Olver, "A higher order nonlinear evolution equation for broader bandwidth gravity waves in deep water," Journal of Fluid Mechanics, vol. 125, pp. 137-185, 1982.

[14] V. E. Zakharov, "Stability of periodic waves of finite amplitude on the surface of a deep fluid," Journal of Applied Mechanics and Technical Physics, vol. 9, no. 2, pp. 190-194, 1968.

[15] V P. Krasitskii, "On reduced equations in the Hamiltonian theory of weakly nonlinear surface waves," Journal of Fluid Mechanics, vol. 272, pp. 1-21, 1994.

[16] S. S. Debsarma and K. P. Das, "Fourth-order nonlinear evolution equations for a capillary-gravity wave packet in the presence of another wave packet in deep water," Physics of Fluids, vol. 19, no. 9, Article ID 097101, 2007.

[17] S. Debsarma and K. P. Das, "A higher-order nonlinear evolution equation for broader bandwidth gravity waves in deep water," Physics of Fluids, vol. 17, no. 10, Article ID 104101, 2005.

S. Debsarma, (1) S. Senapati, (2) and K. P. Das (1)

(1) Department of Applied Mathematics, University of Calcutta, 92, A.P.C. Road, Kolkata 700 009, India

(2) Department of Mathematics, Abhedananda Mahavidyalaya, Sainthia, Birbhum 731234, India

Correspondence should be addressed to S. Debsarma; suma_debsarma@rediffmail.com

Received 14 February 2014; Accepted 2 May 2014; Published 9 June 2014

Academic Editor: Leonard Pietrafesa

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Debsarma, S.; Senapati, S.; Das, K.P. |

Publication: | International Journal of Oceanography |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 3016 |

Previous Article: | Noisy reverberation suppression using AdaBoost based EMD in underwater scenario. |

Next Article: | Sediment chemistry and meiofauna from the northern Gulf of Mexico continental shelf. |

Topics: |