New Complex Hyperbolic Function Solutions for the ( 2 + 1 )-Dimensional Dispersive Long Water – Wave System

Abstract: In this paper, new algorithms called the “Modified exp( ́Ω)-expansion function method” and “Improved Bernoulli sub-equation function method” have been proposed. The first algorithm is based on the exp( ́Ω(ξ))-expansion method; the latter is based on the Bernoulli sub-Ordinary Differential Equation method. The methods proposed have been expressed comprehensively in this manuscript. The analytical solutions and application results are presented by drawing the twoand three-dimensional surfaces of solutions such as hyperbolic, complex, trigonometric and exponential solutions for the (2+1)-dimensional dispersive long water–wave system. Finally, a conclusion has been presented by mentioning the important discoveries in this manuscript.


Introduction
In recent years, studies on the nonlinear differential equations (NDEs) and systems (NDESs) have become a significant field among scientists.This is why many engineering problems can be represented by using NDEs and with the rapid development of computational algorithms.Therefore, many mathematical models have been proposed to understand of such problems in scientific and engineering such as the optical fiber communications, coastal and oceans engineering, fluid dynamics, plasma physics, chemical physics and many other scientific applications.Yong-Sik Cho, Dae-Hee Sohn and Seung Oh Lee have conducted a study on shallow-water equations for distant propagation of tsunamis [1].Furthermore, even some important diseases such as dengue epidemics, tuberculosis, HIV, AIDS, tsunamis, malaria, and cholera have been investigated by many scientists from all over the word [2][3][4][5][6][7].One of them is to obtain various solutions of coastal, oceans and fluid problems such as approximate, numerical, analytical and traveling wave solutions.Some important traveling wave solutions for nonlinear differential equations have been investigated by different authors [8][9][10][11][12][13][14][15].In the rest of this manuscript; we have explained the fundamental properties of the modified exp(´Ω(ξ))-expansion function method (MEFM) and Improved Bernoulli sub-equation function method (IBSEFM) in Section 2. We have studied to obtain some new analytical solutions such as hyperbolic, exponential, and complex hyperbolic function solutions by applying MEFM and IBSEFM to the (2+1)-dimensional dispersive long water-wave (DLW) system defined by [16]: u yt `uxxy ´2v xx ´`u 2 ˘xy " 0, v t ´vxx ´2 puvq x " 0, (1) where u = u(x, y, t) represents the horizontal velocity of water and v = v (x, y, t) gives the deviation height from the equilibrium position of the liquid [16].The solutions of (1) are very helpful for coastal scientists and engineers to apply the nonlinear water model to coastal and harbor design [17][18][19].Equation ( 1) was used to model nonlinear and dispersive long gravity waves traveling in two horizontal directions on shallow waters of uniform depth.Equation (1) appears in many scientific applications such as nonlinear fiber optics, plasma physics, fluid dynamics, and coastal engineering [16].
Step 3: According to different values of N and M, we can find various forms of Equation (6).By putting Equations ( 6) in (5), we can obtain a polynomial of exp(´Ω(ξ)).If we equal to zero all the coefficients of same power of exp(´Ω(ξ)), we can obtain a system of equations.When we solve this system with the aid of Wolfram Mathematica 9, this process gives different values of coefficients A 0 ,A 1 , . . .,A N , and B 0 ,B 1 , . . .,B M , ε.If we put these coefficients in Equation ( 6) by considering family conditions, we obtain many analytical solutions for Equation (2).

Improved Bernoulli Sub-Equation Function Method (IBSEFM)
IBSEFM formed by improving the Bernoulli sub-equation function method [25] will be given in this sub-section.Therefore, we consider the following steps.
Step 1.Let us consider the nonlinear partial differential equations system as the following: and take the wave transformation u px, y, tq " U pξq , ξ " kx `wy ´ct, v px, y, tq " V pξq , ξ " kx `wy ´ct, where w, k, c are both constants and non-zero and will be determined later.Substituting Equations ( 14) into (13) converts a nonlinear ordinary differential equation (NODE) as the following: where NODE is a polynomial of U, and, Step 2. Considering the trial equation of the solution in Equation ( 15), it can be written as the following: and, according to the Bernoulli theory, where F " F pξq is Bernoulli differential polynomial and F 1 " dF dξ " dF pξq dξ .Substituting the above relations in Equation ( 15) yields an equation of polynomial Ω pFq of F " F pξq According to the balance principle, we can get some values of n, m and M.
Step 3. Letting the coefficients of Ω(F) all be zero will yield an algebraic equations system Solving this system, we will determine the values of k, w, c, a 0 ,a 1 ,a 2 , . . .Step 4. When we solve nonlinear Bernoulli differential Equation (17), we obtain, according to b and d, two situations as the following: where ε is both an integration constant and non-zero.Using a complete discrimination system for polynomial of F(ξ), we solve Equation ( 16) with the help of using Wolfram Mathematica 9 and classify the exact solutions to Equation (16).For a better interpretations of results obtained in this way, we plot two-and three-dimensional surfaces of the solutions obtained by taking suitable parameters.

The Implementations of Techniques
In this subsection of the paper, we have obtained new analytical solutions to the DLW system (1) by using MEFM and IBSEFM.

Application of MEFM
Let us consider the travelling wave solutions of the DLW system (1) and we perform the transformation as the following; u px, y, tq " U pξq , ξ " kx `wy ´ct, v px, y, tq " V pξq , ξ " kx `wy ´ct. ( First of all, if we consider the following transformations for the DLW system (1), v px, y, tq " u y px, y, tq " Bu px, y, tq By , (23) we can obtain following partial differential equations: When we use the travelling wave transformations Equation (22), and using Equations ( 26) in (24), we obtain the following NODE; If we integrate Equation (27) twice by getting zero the integration constants, we obtain the following NODE ´cU ´k2 U 1 ´kU 2 " 0.
When we use Equations ( 26) in (25), after integrating twice by getting zero the integration constants we obtain the same NODE as Equation (28) below: When we rearrange to Equation (6), with the help of balance principle between U' and U 2 , we obtain the term for suitability: This relationship gives rise to various analytical solutions to the DLW system (1) as the following:

Application of IBSEFM
Let us consider the travelling wave solutions of the DLW system (1), and we perform the transformation as the following: u px, y, tq " U pξq , ξ " kx `wy ´ct, v px, y, tq " V pξq , ξ " kx `wy ´ct. ( First of all, if we consider the above transformations for the DLW system (1), we can obtain the same NODE as Equation (29): When we rearrange to Equation ( 16), with the help of balance principle between U 1 and U 2 , we obtain the term for suitability: This relationship allows us various analytical solution forms for the DLW system (1).

Case-1
If we choose M " n " 3 and m " 1, we can write the following equations for Equation ( 6): and where a 3 ‰ 0, b 1 ‰ 0 and F " F pξq .When we use Equations ( 52) and (53) in Equation ( 50), we get a system of equations for Equation (50) from the coefficients of polynomial of F. Solving this system with the help of Wolfram Mathematica 9 yields the following coefficients: Case 1.1: If Equations (54) along with ( 8) is put into (51) for u px, y, tq and in Equation ( 23) for v px, y, tq, new complex exponential function solution to the DLW system(1) can be obtained as the following: When we use coefficients Equations (55) along with (8) in (51) for u px, y, tq and in Equation ( 23) for v px, y, tq, we can obtain another new exponential function solution to the DLW system (1) as the following: b ?

Conclusions
The MEFM and IBSEFM have been applied to the (2+1)-dimensional dispersive long water-wave system Equation (1).Then, these methods which are newly submitted to literature in this paper have given new complex and exponential functional solutions such as Equations ( 39)-( 41) and ( 57)-( 59) for the DLW system Equation (1).We have shown that these analytical solutions are verified in Equation ( 1).Furthermore, we have already found the similar hyperbolic function solutions such as Equations ( 38), ( 42), ( 47), (48), and (58) with analytical solutions obtained by Wazwaz [16].Moreover, complex function solutions such as Equations ( 39)-( 41), (57), and (59) obtained in this paper by MEFM and IBSEFM are new analytical solutions when compared with analytical solutions obtained by Wazwaz [16].These results are very helpful for coastal and civil engineers to apply the nonlinear water model to coastal and harbor design [16].Figures 1-15 have been plotted by using Wolfram Mathematica 9 under the suitable values of parameters.It can be predicted that they are very useful for applications in engineering and science, especially, coastal and ocean engineering.To the best of our knowledge, applications of MEFM and IBSEFM to the DLW system (1) have not previously been submitted to literature.We think that these newly modified methods can also be applied to other models that result from engineering, science, coastal and ocean engineering.Case 1.5: Therefore, when we substitute coefficients Equations (33) along with Equation (8) in Equation (31) for   ,, u x y t and in Equation ( 23) for   ,, v x y t , we can obtain hyperbolic function solution to the DLW system (1) as the following, under the condition of Family-1; ,  23) for   ,, v x y t , we can obtain new complex hyperbolic function solution to the DLW system(1) as the following, under the condition of Family-1; where f         23) for  

x y t f x y t x y t i cw x y t
,, v x y t , we can obtain another new complex hyperbolic function solution to the DLW system (1) as the following, under the condition of Family-1; Taking Equation (37) along with Equation (8) in Equation (31) for  ,, u x y t and in Equation ( 23) for   ,, v x y t , we can obtain another new hyperbolic function solution to the DLW system(1) as the following, under the condition of Family-1; , where
where  

Application of IBSEFM
Let us consider the travelling wave solutions of the DLW system (1), and we perform the transformation as the following: First of all, if we consider the above transformations for the DLW system (1), we can obtain the same NODE as Equation (29): When we rearrange to Equation ( 16), with the help of balance principle betweenU  and 2 U , we obtain the term for suitability: This relationship allows us various analytical solution forms for the DLW system (1).

Case-1
If we choose 3 Mn and 1 m  , we can write the following equations for Equation ( 6): Case 1.2: Case 1.3: If Equation (54) along with Equation ( 8) is put into Equation (51) for   ,, u x y t and in Equation ( 23) for   ,, v x y t , new complex exponential function solution to the DLW system(1) can be obtained as the following: .

Conclusions
The MEFM and IBSEFM have been applied to the (2 + 1)-dimensional dispersive long water-wave system Equation (1).Then, these methods which are newly submitted to literature in this paper have given new complex and exponential functional solutions such as Equations ( 39)-( 41) and ( 57)-(59) for the DLW system Equation (1).We have shown that these analytical solutions are verified in Equation (1).Furthermore, we have already found the similar hyperbolic function solutions such as Equations ( 38), ( 42), (47), (48), and (58) with analytical solutions obtained by Wazwaz [16].Moreover, complex function solutions such as Equations ( 39)-( 41), (57), and (59) obtained in this paper by MEFM and IBSEFM are new analytical solutions when compared with analytical solutions obtained by Wazwaz [16].These results are very helpful for coastal and civil a n , b 0 ,b 1 ,b 2 , . . .,b m .

Figure 8 .
Figure 8.The 3D and 2D surfaces of Equation (42) by considering the values 0 where When we use Equations (43) and (44) in the Equation (29), we get a system of equations for Equation (29) from the coefficients of polynomial of with the help of Wolfram Mathematica 9 yields the following coefficients: Case 2.1:
46) When considering Equation (45) along with Equation (8) in Equation (31) for   ,, u x y t and in Equation (23) for   ,, v x y t , another new hyperbolic function solution to the DLW system(1) can be obtained as the following, under the condition of Family-1; 2 4 0
Substituting Equations (56) along with (8) in (51) for u px, y, tq and in Equation (23) for v px, y, tq, we can obtain another new complex exponential function solution to the DLW system(1) as the following: 40 R e ux,y,t