Computational study of some nonlinear shallow water equations

Transcription

1 Shiraz University of Technology From the SelectedWorks of Habibolla Latifizadeh 013 Computational study of some nonlinear shallow water equations Habibolla Latifizadeh, Shiraz University of Technology Available at:

2 Cent. Eur. J. Phys DOI: /s Central European Journal of Physics Computational study of some nonlinear shallow water equations Research Article Ali H. Bhrawy 1,, Mohamed A. Abdelkawy 1 Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 1589, Saudi Arabia Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef 6511, Egypt Received 09 December 01; accepted 7 February 013 Abstract: The shallow water equations have wide applications in ocean, atmospheric modeling and hydraulic engineering, also they can be used to model flows in rivers and coastal areas. In this article we obtained exact solutions of three equations of shallow water by using -expansion method. Hyperbolic and triangular periodic solutions can be obtained from the -expansion method. PACS 008: 0.30.p, 0.30.Jr, c, 0.70.Wz Keywords: -expansion method nonlinear physical phenomena nonlinear shallow water equations time dependent nonlinear system of shallow water Versita sp. z o.o. 1. Introduction The shallow water equations SWEs are a set of hyperbolic partial differential equations that describe the free-surface flows arising in shores, rivers, estuaries, tidal flats, coastal regions and also SWEs describe many other physical situations [1] such as storm surges, tidal fluctuations, tsunami waves, forces acting on off-shore structures and modeling the transport of chemical species. The SWEs, which are based on the conservation laws of mass and momentum, may be derived by integrating in depth the Reynolds averaged Navier-Stokes equations []. The melkawy@yahoo.com speed of shallow water waves SWWs is independent of wavelength or wave period and is controlled by the depth of water, also all SWW travel at the same speed [3]. In contrast, deep water waves of different length travel at different speeds the long ones faster than the short ones. The shallow water equations have been solved in many articles numerically and analytically. Variational iteration method RVIM is applied in [3] to compute the approximate solution for coupled Whitham-Broer-Kaup shallow water. While in [4], Benkhaldoun et al. proposed finite volume method for numerically solve shallow water equations on nonflat topography. In [5], the authors presented spectral alerkin method for the numerical solution of the shallow water equations in two dimensions. The shallow water equations over irregular domains with wetting and drying are discussed by second-order Runge-Kutta 1

3 Computational study of some nonlinear shallow water equations discontinuous alerkin scheme in [6]. Moreover, the exact solutions for shallow water equations are obtained in several articles suing different methods such as generally projective Riccati equation method [7], Bäcklund transformation and Lax pairs [8], the tanh-coth, Exp-function and Hirota s methods [9, 10]. We, firstly, will study time dependent nonlinear system of shallow water SW. The authors of [11 13] studied the coupled Whitham-Broer-Kaup CWBK equations that describe the dispersion of SWW with different scattering relations: u t + uu x + v x + βu xx = 0, v t + uv x + αu xxx βv xx = 0, 1.1 where u is the velocity and v the total depth, which can be used as a model for water waves. When α = 1 and β = 0, the CWBK equations are reduced to the approximate long-wave ALW equations in SW. When α = 0 and β = 0.5, the CWBK equations are reduced to the modified Boussinesq MB equations [14]. The second equation that we will study in this paper is the 1+1-dimensional SWE in the following form [15]: u xt + u xx u xxxt 4u x u xt u xx u t = The third one, we will study in this paper is the +1- dimensional SWE in the following form: 4u xt + u xxxy + 8u xy u x + 4u xx u y = A large variety of physical, chemical, and biological phenomena is governed by nonlinear partial differential equations NPDEs. The exact solutions of these NPDEs plays an important role in the study of nonlinear phenomena. In the past decades, many methods were developed for finding exact solutions of NPDEs as the inverse scattering method [16]- [18], Hirota s bilinear method [19], new similarity transformation method [0], homogeneous balance method [1, ], the sine-cosine method [3, 5], tanh function methods [6] [3], exp-function method [33] [36], Jacobi and Weierstrass elliptic function method [37] [39]... etc. In this article we find exact solutions of three SWEs using the expansion method [40] [44]. Kudryashov [45] proposed the -expansion method to find rational solutions of Schwarzian integrable hierarchies. Moreover, Wang et al. [46] developed this approach for the exact solution of nonlinear wave equations. The main advantages of the -expansion method, that the obtained solutions are more general solutions with some free parameters, while we can obtained other solution by taking different choice of the parameters. In addition 1. The -expansion method is a powerful technique to integrate nonlinear differential equations without need of boundary conditions or initial trial function. In contrast, in all finite difference, finite element and spectral methods, it is necessary to have boundary or initial conditions,. The -expansion method gives a general solution without approximation, in contrast, other methods give solutions in a series form and it becomes essential to investigate the accuracy and convergence of analytical and numerical methods, 3. In case of the Painlevè test of integrability fails, -expansion method can be used as a powerful technique to integrate the NEEs, 4. The rational solutions are computed by using - expansion method and they could not be obtained by other methods.. expansion method This section is devoted to the study of implementing the expansion method for a given partial differential equation u, u x, u y, u t, u xy,... = 0,.1 where ux, y, t is an unknown function in the independent variables x, y and t. In order to obtain the solution of Eq..1, we combine the independent variables x, y and t into one particular variable through the new variable ζ = x + y + t, ux, y, t = Uζ, where is the wave speed, using this variable, enables us to reduce Eq..1 to the following ordinary differential equation ODE U, U, U, U,... = 0.. The ODE is integrated as long as all terms contain derivatives in ζ, upon setting the constant of integration to zero, one is looking for although the non-zero case can be treated as well. At this stage, we search for the exact solutions which satisfy this ODE. For this purpose, Uζ can be expressed as a finite series of, ux, y, t = Uζ = N i=0 i a i, + µ + λ = 0..3

4 Ali H. Bhrawy, Mohamed A. Abdelkawy The parameter N can be determined by balancing the linear terms of highest order with the nonlinear ones, where N is a positive integer, so that an analytic solution in closed form may be investigated. Expressing the ODE. in terms of Eq..3 and comparing the coefficients of each power of in both sides, to get an over-determined system of nonlinear algebraic equations in the parameters, µ, λ and a i, i = 1,, N. "We next solve the overdetermined system of nonlinear algebraic equations by using Mathematica program. 3. Nonlinear system of shallow water In this section, we firstly implement the proposed method for solving the nonlinear system of SW: u t + uu x + v x + βu xx = 0, v t + uv x + αu xxx βv xx = 0, 3.1 as described in the introductory section, we perform a travelling wave reduction, ux, t = Uζ, vx, t = V ζ, ζ = x + t, 3. that converted 3.1 into a system of ODEs, U + UU + V + βu = 0, V + UV + αu βv = 0, 3.3 if we integrate the system 3.3 once, upon setting the constant of integration to zero, we obtain U + U + V + βu = 0, V + UV + αu βv = If we use the first equation in system 3.4 into the second one, we get only one ODE U + 3U + U 3 β + αu = Balancing the term U with the term U 3, implies to N = 1. Accordingly, the solution takes on the form U ζ = 1 a i i=0 i. 3.6 Substituting Eq. 3.6 into Eq. 3.5 and comparing the coefficients of each power of in both sides, getting an overdetermined system of nonlinear algebraic equations with respect to, µ, λ and a i ; i = 0, 1. Solving the over-determined system of nonlinear algebraic equations using Mathematica, we obtain two sets of constants: a µ = 0, a 0 = a 1 = ± α + β, and λ = 4 α + β, 3.7 3

5 Computational study of some nonlinear shallow water equations b a 0 = αµ + β µ a 1 = ± α + β, and λ = αµ + β µ α + β We find the following solutions of Eq. 3.5 For case a i α+β > 0 ii α+β < 0 iii = 0 U 1 = 1 4 U = 1 4 α+β α+β α+β α+β ζ + c cosh ζ + c sinh ζ + c cos ζ + c sin α+β α+β α+β α+β ζ, 3.9 ζ ζ, 3.10 ζ U 3 = ± c 1 α + β c 1 ζ + c For case b i α+β > 0 ii α+β < 0 iii = 0 U 4 = αµ + β µ ± 4 U 5 = αµ + β µ ± 4 α+β ζ + c cosh α+β ζ + c sinh α+β ζ + c cos α+β ζ + c sin α+β ζ, 3.1 α+β ζ α+β ζ, 3.13 α+β ζ U 6 = αµ + β µ ± c 1 α + β c 1 ζ + c Then the solutions of the system 3.1 are: 4 u 1 = 1 4 α+β x + t 4 α+β x + t + c sinh u = 1 α+β x + t + c cosh + c cos α+β x + t + c sin α+β x + t α+β x + t α+β x + t α+β x + t, 3.15,

6 Ali H. Bhrawy, Mohamed A. Abdelkawy u 3 = ± c 1 α + β, 3.17 c 1 x + c u 4 = α+β x + t + c cosh αµ + β µ ± α+β x + t + c sinh u 5 = αµ + β µ ± 4 where v in all above solutions given by the relation: α+β x + t + c cos α+β x + t + c sin α+β x + t, 3.18 α+β x + t α+β x + t, 3.19 α+β x + t u 6 = αµ + β µ ± c 1 α + β, 3.0 c 1 x + c V = U + U + βu. 3.1 t x a b Fig. 1 a Three-dimensional of the modulus of solitary wave solution u Eq where α = β = 0.5, µ = 1, = 1, b raphical representation of of the modulus of solitary wave solution u Eq where α = β = 0.5, µ = 1, = 1. The soliton solutions for system 3.1 are given in 3.15 and 3.18, while the triangular periodic solutions are given in 3.16 and Moreover, for the special values β = c = 0, α = 1, we obtain equivalent solutions for the system of the SWW equation that given in [47]. Three-dimensional of the modulus of solitary wave solutions u 1, Eq and u, Eq are displayed in Figures 1a and a respectively with values of parameters listed in their captions. Moreover, graphical representation for the same solutions are plotted in Figures 1b and b respectively dimensional shallow water equation We next consider the 1+1-dimensional SWE in the following form: u xt + u xx u xxxt 4u x u xt u xx u t = 0, 4.1 if we use the transformations ux, t = Uζ, ζ = x + t, 4. it carries Eq. 4.1 into the ODE + 1U U 6U U = 0, 4.3 5

7 Computational study of some nonlinear shallow water equations t x a b Fig. a Three-dimensional of the modulus of solitary wave solution u Eq where α = β = 0.5, µ = 1, = 1, b raphical representation of of the modulus of solitary wave solution u Eq where α = β = 0.5, µ = 1, = 1. where by integrating 4.3 once and taking the constants of integration to be zero we obtain + 1 U U 3 U = Considering the homogeneous balance between U and U, then N = 1, and U ζ = 1 a i i=0 Proceeding as in the previous case we obtain a b c i. 4.5 λ = µ 1, and a 1 =, µ = ±, λ = 0 and a 1 =, 4.7 λ = 1 + 4, µ = 0 and a 1 =. 4.8 For all above cases we find the same three solutions of Eq. 4.4: i > 0 U 1 = a 0 + ii < 0 U = a 0 + iii = 1 ζ + c cosh ζ + c sinh ζ + c cos ζ + c sin ζ, 4.9 ζ ζ, 4.10 ζ U 3 = a 0 + c 1 c 1 ζ + c Then the solutions of the Eq. 4.1 are: u 1 = a 0 + u = a 0 x + t + c cosh x + t, x + t + c sinh x + t 4.1 x + t + c cos x + t +, x + t + c sin x + t 4.13 c 1 u 3 = a c 1 x t + c The solution 4.1 is a soliton solution. Moreover the triangular periodic solution is obtained in For the special values c = 0, we give equivalent solutions for 1+1-dimensional shallow water equation to that given in [48]. 6

8 Ali H. Bhrawy, Mohamed A. Abdelkawy dimensional shallow water equation In this section, we study the +1-dimensional SWE in the following form: ii < 0 U = a 0 + ζ + c cos ζ ζ + c sin, 5.10 ζ 4u xt + u xxxy + 8u xy u x + 4u xx u y = 0, 5.1 if we use the following transformations iii = 0 U 3 = a 0 + c 1 c 1 ζ + c ux, t = U ζ, ζ = x + t, 5. inserting 5. into 5.1 yields an ODE. 4U + U + 8U U + 4U U = 0, 5.3 if we integrate the Eq. 5.3 once, upon setting the constant of integration to zero, we obtain Balancing the term U N = 1, and 4U + U + 6 U = U ζ = with the term 1 a i i=0 U, we obtain i. 5.5 Proceeding as in the previous case we obtain three sets of constants a b c λ = µ + 4, and a 1 = 1, µ = ±i, λ = 0 and a 1 = 1, 5.7 λ =,, µ = 0 and a 1 = For all above cases we find the same three solutions of Eq. 5.4: i > 0 U 1 = a 0 + ζ + c cosh ζ ζ + c sinh, 5.9 ζ Then the solutions of the equation 5.1 are: u 1 = a 0 + x + y + t + c cosh x + y + t x + y + t + c sinh, x + y + t 5.1 u = a 0 + x + y + t + c cos x + y + t x + y + t + c sin, x + y + t u 3 = a c 1 c 1 x + y + c We conclude that 5.1 is a soliton solution, while 5.13 is a triangular periodic solution. Moreover, we can obtain other soliton solutions that may be derived from the soliton solutions obtained in the last three sections, by setting c 1 = 0 or c = 0. The same procedure can be applied for triangular periodic solutions. 6. Conclusion We extend the -expansion method with symbolic computation to two nonlinear equations of SW and nonlinear system of SW. It is shown that hyperbolic solutions and triangular periodic solutions can be established by the - expansion method. The travelling wave solutions can be expressed as the hyperbolic functions, trigonometric functions and rational functions, and involved in two arbitrary parameters. As special values of parameters, the solitary waves are also derived from the travelling waves. The rational solutions are computed only by using -expansion method and they could not be obtained by other methods. Acknowledgments The authors are very grateful to the referees for carefully reading the paper and for their comments and suggestions which have improved the paper. 7

9 Computational study of some nonlinear shallow water equations References [1] R. Bernetti, V. A. Titarev, E. F. Toro, J. Comput. Phys. 7, [] Q. Liang, F. Marche, Adv. Water Resour. 3, [3] A. A. Imani, D. D. anji, H. B. Rokni, H. Latifizadeh, E. Hesameddini, M. H. Rafiee, Appl. Math. Modell. 36, [4] F. Benkhaldoun, I. Elmahi, M. Seaïd, Comput. Methods Appl. Mech. Engrg. 199, [5] T. Kröger, M. Lukáčovà-Medvid ovà, J. Comput. Phys. 06, [6]. Kesserwani, Q. Liang, Computers & Fluids 39, [7] Q. Wang, Y. Chen, B. Li, H. Zhang, Appl. Math. Comput. 160, [8] Y. Shang, Appl. Math. Comput. 187, [9] A. M. Wazwaz, Appl. Math. Comput. 01, [10] A. M. Wazwaz, Appl. Math. Comput. 0, [11]. B. Whitham, Proc. Roy. Soc. Lond. A 99, [1] L. J. F. Broer, Appl. Sci. Res. 31, [13] D. J. Kaup, Prog. Theor. Phys. 54, [14] Xie Fuding, Zhenya Yan, Phys. Lett. A 85, [15] A. M. Wazwaz, Appl. Math. Comput. 17, [16] S. hosh, S. Nandy, Nucl. Phys. B 561, [17] T. Koikawa, Phys. Lett. B 110, [18] R. A. Baldock, B. A. Robson, R. F. Barrett, Nucl. Phys. A 366, [19] A. M. Wazwaz, Appl. Math. Comput. 00, [0] A. N. Beavers Jr., E. D. Denman, Math. Biosci. 1, [1] M. L. Wang, Y. B. Zhou, Z. B. Li, Phys. Lett. A 16, [] J. F. Zhang, Phys. Lett. A 313, [3] A. M. Wazwaz, Appl. Math. Comput. 173, [4] M. Y. Moghaddam, A. Asgari, H. Yazdani, Appl. Math. Comput. 10, [5] S. Tang, C. Li, K. Zhang, Commun. Nonlinear Sci. Numer. Simul. 15, [6] E. Fan, Y. C. Hon, Z. Naturforsch. 57, [7] W. Malfliet, W. Hereman, Phys. Scr. 54, [8] W. Malfliet, W. Hereman, Phys. Scr. 54, [9] A. H. Khater, D. K. Callebaut, M. A. Abdelkawy, Phys. Plasmas 17, [30] A. H. Khater, D. K. Callebaut, A. H. Bhrawy, M. A. Abdelkawy, J. Comput. Appl. Math. 4, [31] E.. Fan, Phys. Lett. A 77, [3] A. M. Wazwaz, Commun. Nonlinear Sci. Numer. Simul. 11, [33] H. Naher, F. A. Abdullah, M. A. Akbar, Int. J. Phys. Sci. 6, [34] R. Sakthivel, C. Chun, Z. Naturforsch. A 65, [35] R. Sakthivel, C. Chun, J. Lee, Z. Naturforsch. A 65, [36] J. Lee, R. Sakthivel, Rep. Math. Phys. 68, [37] S. K. Liu, Z. T. Fu, S. D. Liu, Phys. Lett. A 89, [38] A. H. Bhrawy, M. A. Abdelkawy, A. Biswas, Commun. Nonlinear Sci. Numer. Simul. 18, [39] A. H. Bhrawy, A. Yildirim, M. M. Tharwat, M. A. Abdelkawy, Indian J. Phys. 86, [40] A. Malik, F. Chand, H. Kumar, S. C. Mishra, Indian J. Phys. 86, [41] A. Malik, F. Chand, S. C. Mishra, Appl. Math. & Comp. 16, [4] A. Malik, F. Chand, H. Kumar, S. C. Mishra, Pramana- J. Phys. 78, [43] A. Malik, F. Chand, H. Kumar, S. C. Mishra, Comp. & Math. Appl. 64, [44] F. Chand, A. K Malik, International Journal of Nonlinear Science, 14, [45] N. A. Kudryashov, A. Pickering, J. Phys. A: Math. en. 31, [46] M. L. Wang, X. Li, J. Zhang, Phys. Lett. A 37, [47] I. E. Inan, D. Kaya, Phys. Lett. A 355, [48] S. A. Elwakil, S. K. El-labany, M. A. Zahran, R. Sabry, Chaos. Soliton. Fract. 17,