Mathematical Problems in Engineering

Volume 2018 (2018), Article ID 4276591, 12 pages

https://doi.org/10.1155/2018/4276591

## Three-Dimensional Coupled NLS Equations for Envelope Gravity Solitary Waves in Baroclinic Atmosphere and Modulational Instability

Correspondence should be addressed to Ruyun Wang

Received 3 August 2017; Accepted 30 October 2017; Published 8 January 2018

Academic Editor: Qin Yuming

Copyright © 2018 Baojun Zhao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Envelope gravity solitary waves are an important research hot spot in the field of solitary wave. And the weakly nonlinear model equations system is a part of the research of envelope gravity solitary waves. Because of the lack of technology and theory, previous studies tried hard to reduce the variable numbers and constructed the two-dimensional model in barotropic atmosphere and could only describe the propagation feature in a direction. But for the propagation of envelope gravity solitary waves in real ocean ridges and atmospheric mountains, the three-dimensional model is more appropriate. Meanwhile, the baroclinic problem of atmosphere is also an inevitable topic. In the paper, the three-dimensional coupled nonlinear Schrödinger (CNLS) equations are presented to describe the evolution of envelope gravity solitary waves in baroclinic atmosphere, which are derived from the basic dynamic equations by employing perturbation and multiscale methods. The model overcomes two disadvantages: (1) baroclinic problem and (2) propagation path problem. Then, based on trial function method, we deduce the solution of the CNLS equations. Finally, modulational instability of wave trains is also discussed.

#### 1. Introduction

According to media reported, “6.1 the Oriental Star cruise ship capsized event” was caused by the violence rainstorm attack which was brought by a sudden squall lines process. Unfortunately, the forecast for the severe weather such as squall lines is very difficult. As we know, the nonlinearity concentration of the gravity waves makes the energy assemble together and forms disastrous weather phenomena, such as squall lines and rainstorm. So constructing the theoretical model of gravity waves suits the real atmosphere condition and analyzing the formation mechanism of disastrous weather phenomena based on the theoretical model has significant scientific meaning and application value.

The first discovery of solitary waves was attributed to Russell [1]; then, the study of solitary wave continued to deepen. There was vast research literature devoted to the solitary waves. Solitary waves in the westerly shear flow were first found by Long [2]. Afterward, Benny [3] amplified the conclusions and got a conclusion that velocity and amplitude of solitary waves were related. Recently, a variety of equation models which described the solitary waves, such as ILW-Burgers equation and ZK-Burgers equation, were discussed by Yang et al. [4, 5]. Meanwhile the generation and evolution of solitary waves in different topography condition and different fluid depths were discussed. In recent years, in the solitary waves community, gravity solitary wave as a rising star had been paid more and more attention by many scientists. In the 60s, Long [2] proposed that the amplitude of atmospheric gravity wave satisfied KdV equation and got the solutions of amplitude. Then a lot of researchers obtained nonlinear KdV equations of gravity waves from the basic dynamic equation. From the original equation of two layers of -plane, under the shearing basic flow, the famous KdV equation was derived out by Li [6], and he pointed out that the nonlinear characteristic of gravity waves was the generation mechanism of squall lines. S. K. Liu and S. D. Liu [7] based on nonlinear atmospheric motion equation deduced out nonlinear gravity wave solutions and found that gravitational wave amplitude with propagation speed was proportional. Later, Boussinesq equation was derived in Luo [8] to describe the algebraic gravity solitary wave in atmosphere. Additionally, two-dimensional dissipative nonlinear Schrödinger equation [9] was also obtained to reflect the evolution of envelope solitary Rossby waves.

Therefore, for better understanding the characteristic of solitary waves deeply, it is desired to get the exact solutions [10, 11] of soliton equations. Up to now, people have obtained several sorts of nonlinear evolution equations and studied their properties, including the boundary value problem [12–14], Hamiltonian structure [15, 16], integrable systems [17, 18], and conservation law [19]. In the process of solving nonlinear evolution equations, many solution methods are found, like Darboux transformation method [20–22], Hirota method [23], homogeneous balance method [24], Jacobi elliptic function method [25], symmetry method [26], trial function method [27], the alternative variational-asymptotic method [28–31], and so on [32, 33]. Biswas et al. [34–36] have also found some other methods of solving the nonlinear partial equation, but we noticed that they have not been applied to solve the three-dimensional CNLS equations.

Through the analysis of previous research on the gravity wave we can get the following two points:

(1) In order to simplify the calculation, the former researches studied the propagation of gravity solitary waves in two-dimension space and barotropic atmosphere [37]. But for the envelope of gravity waves in the real atmosphere, it was not sufficient to consider only two dimensions. And the initial equation was a simple model in barotropic atmospherewhere , , , and . These variables were simplified. Therefore, the three-dimensional model in baroclinic atmosphere is more in line with the actual atmospheric conditions.

(2) Previous studies have focused on the solution of gravity waves and the dynamic characteristics of the wave equation mainly based on the numerical calculation from the original equation. As we all know, many of the variables in the original equation are simplified and dimensional; therefore, numerical simulation directly based on the original equation is usually imprecise and difficult. So it is necessary for us to find the appropriate model and solve it to study the evolution of the waves.

In this paper, overcoming the limitations of calculations and using the appropriate method we obtain a new model. The paper will be organized as follows: (1) using multiscale analysis and turbulence method, from the basic dynamic equations of multivariable in baroclinic environment, we derived out the gravity wave model in Section 2, which is a coupled nonlinear Schrödinger equation (CNLS). Not only is the model three-dimensional and more suited to describe the feature of two envelope gravity solitary waves in a plane, specially, but also it is a coupled model and can show the interaction process between two waves. (2) Based on the CNLS model, using the trial function method to solve the equations, the analytical solution was obtained in Section 3. By observing the structures of the solution, the evolution characteristics of gravity solitary waves was obtained in Section 4. (3) Finally, modulational instability of a uniform three-dimensional gravity waves trains was also discussed.

#### 2. Derivation of the Three-Dimensional CNLS Equations Group

Using the sum of disturbance pressure gradient force and buoyancy force, express the vertical pressure gradient force and gravity force, adopt Boussinesq approximation, and the basic dynamic equations of atmospheric motion are as follows:in the above equations, is the temperature of environmental flow field and is the density; . Each variable above will be dimensionless. Letwhere is the height of the homogeneous atmosphere, is the characteristics pressure of the ground, is the horizontal pressure changes, and is the pressure changes in the vertical direction. Substituting (3) into (2), get the dimensionless equations as follows:supposing in the above. Because the second term of the left side of the fourth formula is lesser, get the following approximate:introducing parameter (), where . By varying, (4) transforms to the following:

Introduction of multiscale variables (omit the sign at the top right corner of the variables)so long time and space scales are defined as in (6) were expanded according to the small parameter :where are the functions of , where is the speed of the basic flow, is the air pressure, and is the temperature field.

Substituting (8) and (9) into (6), the zero-order approximation of can be obtained:

It is shown that basic flow is geostrophic equilibrium and static balance. For a basic flow we can further assume thatin the first formula in (10), derivation of , and the second formula in (10), derivation of , so that

Further, take the first-order approximation of and introduce the new variables we have

After eliminating other variables in (14), we can get the equation about , as follows: where

Clearly, (15) is a variable separable equation; assume its solution isunder a certain definite solution condition, we can get ; further, all the solutions of (14) can be obtained:

Taking the second-order approximation of and introducing the new variables we havelet

After eliminating other variables in (20), the equation about can be obtained as follows: where

Similarly, (23) is also a variable separable equation; with comparison on both ends of (23), assume its basic solution is further, we can get all solutions of (23):

Similar to the above, introducing the new variables,the third-order approximation of transforms to

Let

Substitute (18) and (26) into (29) and collecting the secular-producing terms proportional to , we have

After eliminating other variables in (28), the equation about can be obtained as follows:

Introduce the following variables:There are conditions for solution of (31),where

#### 3. The Solutions of the Three-Dimensional CNLS Equations Group

In the past, Hui [38] gave a solitary wave solution of dimensional nonlinear Schrödinger equation with the help of a space coordinate transformation. In this chapter, this method will be applied. Based on trial function method, we will discuss the solution of the CNLS equations group. The equations group is

If we further introduce

(35) reduces towhere

By rescaling the spatial coordinate and both amplitudes, we can reduce the number of nontrivial coefficients of the CNLS system; the simplified system iswhere the new unknowns and coefficients are

The steadily oscillating solutions of this are solutions and their spatially periodical generalizations. Solitons and cnoidal waves are obtained by assumingwhere , and , are assumed to be real. Note that, in spite of the cross-modulation, the frequencies in the two modes are generally different. Substitute (41) into (39), a set of two coupled ODES equations’ group can obtained and can be called SCNLS equations group

By rescaling the spatial coordinate and both amplitudes, it is possible to reduce the number of nontrivial coefficients of the SCNLS system from six to three. The simplified system iswhere the new unknowns and coefficients are

To compute the solitary and cnoidal waves, it is always sufficient to solve this system of two ODES. However, to collide solitary waves it is necessary to shift the steadily oscillating solitary waves into a frame of reference such that they move at some translational velocity .

If , are solutions to the CNLS equations groupthe another solution iswhere the velocity is

which implies that the constants and must be proportional as indicated set and definewhere is a solution to

Substitute (48) into the ODE system (42) and collecting terms in and imposing the condition that the coefficients must be equal 2, the matrix system is obtained

The solution is

The limiting solitary waves for both branches are

So the solution of CNLS equations group iswhere are Landau coefficients and are interaction coefficients.

#### 4. Energy Variation Characteristics of Coupled Envelope Gravity Waves

In this section we will discuss two waves energy variation characteristics in the process of interaction; the first equation of (33) is multiplied by ; the conjugate equation is multiplied by , subtracting the two types. In the same way, the second equation of (33) is multiplied by ; the conjugate equation is multiplied by , subtracting the two types, we can obtain the following equations group: is integrated from to ; if and , we can getEquations (55) show that and are time invariants.

Obviously and are energy of two waves. This suggests that the energy of waves is conserved, when the coefficients satisfyor

Equation (33) has infinite conservation law. Except the two kinds of energy conversation and above, (33) have four movement invariants:

Equations (58) to (61) are the momentum of two waves, so we can confirm the conservation of momentum.

#### 5. Modulation Instability of a Uniform Gravity Waves Trains in 3-Dimension Space

We will discuss the modulation instability for three-dimensional gravity wave trains in this section. Based on (53) we can obtain the solution of a uniform wave train. Initial solution of the equation can take the following form:according to the studies of Yuen and Lake [39], the perturbed wave trains may be expressed as follows:where, in fact, denotes the amplitude of a uniform disturbance stream function, is the amplitude of a sideband disturbance and is small enough compared with to permit linearization, is eigenvalue to be determined, and its conjugate is . Informed by Craik [40], denotes the prescribed zonal and meridional wave number perturbations in the and directions. Here, we define , in which and are required. Linearization with respect to satisfied the following eigenvalue equation:

If and are defined, we can get if and satisfy the uniform gravity wave train is unstable.

When , the result resembles the instability diagram of Benjamin-Feir instability [41]. The case of has not been studied yet. To obtain more detailed results, for gravity solitary waves, we only study the case where the fundamental number is 1; that is,where is the latitude. We will discuss the influence of and latitude value on the instability regions. Taking and . In this case the instability growth rates 35 of a uniform gravity wave train at 30°N, 45°N, and 60°N for are shown in Figures 1, 2, and 3.