MPEMathematical Problems in Engineering1563-51471024-123XHindawi Publishing Corporation50939010.1155/2009/509390509390Research ArticleSoliton and Periodic Wave Solutions to the OsmosisK(2,2) EquationZhouJiangboTian LixinFanXinghuaLiaoShijunNonlinear Scientific Research CenterFaculty of ScienceJiangsu UniversityZhenjiang, Jiangsu 212013Chinaujs.edu.cn200902092009200907052009090720092009Copyright © 2009This 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.

Two types of traveling wave solutions to the osmosis K(2, 2) equation ut+(u2)x(u2)xxx=0 are investigated. They are characterized by two parameters. The expresssions for the soliton and periodic wave solutions are obtained.

1. Introduction

In 1993, Rosenau and Hyman  introduced a genuinely nonlinear dispersive equation, a special type of KdV equation, of the form

ut+a(un)x+(un)xxx=0,n>1, where a is a constant and both the convection term (un)x and the dispersion effect term (un)xxx are nonlinear. These equations arise in the process of understanding the role of nonlinear dispersion in the formation of structures like liquid drops. Rosenau and Hyman derived solutions called compactons to (1.1) and showed that while compactons are the essence of the focusing branch where a>0, spikes, peaks, and cusps are the hallmark of the defocusing branch where a<0 which also supports the motion of kinks. Further, the negative branch, where a<0, was found to give rise to solitary patterns having cusps or infinite slopes. The focusing branch and the defocusing branch represent two different models, each leading to a different physical structure. Many powerful methods were applied to construct the exact solutions to (1.1), such as Adomain method , homotopy perturbation method , Exp-function method , variational iteration method , and variational method [6, 7]. In , Wazwaz studied a generalized forms of (1.1), that is mK(n,n) equations and defined by

un-1ut+a(un)x+b(un)xxx=0,n>1, where a,b are constants. He showed how to construct compact and noncompact solutions to (1.2) and discussed it in higher-dimensional spaces in . Chen et al.  showed how to construct the general solutions and some special exact solutions to (1.2) in higher-dimensional spatial domains. He et al.  considered the bifurcation behavior of traveling wave solutions to (1.2). Under different parametric conditions, smooth and nonsmooth periodic wave solutions, solitary wave solutions, and kink and antikink wave solutions were obtained. Yan  further extended (1.2) to be a more general form

um-1ut+a(un)x+b(uk)xxx=0,nk1. And using some direct ansatze, some abundant new compacton solutions, solitary wave solutions and periodic wave solutions to (1.3) were obtained. By using some transformations, Yan  obtained some Jacobi elliptic function solutions to (1.3). Biswas  obtained 1-soliton solution of equation with the generalized evolution term

(ul)t+a(um)ux+b(un)xxx=0, where a,b are constants, while l,m, and n are positive integers. Zhu et al.  applied the decomposition method and symbolic computation system to develop some new exact solitary wave solutions to the K(2,2,1) equation

ut+(u2)x-(u2)xxx+uxxxxx=0, and the K(3,3,1) equation

ut+(u3)x-(u3)xxx+uxxxxx=0. Recently, Xu and Tian  introduced the osmosis K(2,2) equation

ut+(u2)x-(u2)xxx=0, where the positive convection term (u2)x means the convection moves along the motion direction, and the negative dispersive term (u2)xxx denotes the contracting dispersion. They obtained the peaked solitary wave solution and the periodic cusp wave solution to (1.7). In , the authors obtained the smooth soliton solutions to (1.7). In this paper, following Vakhnenko and Parkes's strategy [18, 19] we continue to investigate the traveling wave solutions to (1.7) and obtain soliton and periodic wave solutions. Our work in this paper covers and extends the results in [16, 17] and may help people to know deeply the described physical process and possible applications of the osmosis K(2,2) equation.

The remainder of this paper is organized as follows. In Section 2, for completeness and readability, we repeat [19, Appendix  A], which discusses the solutions to a first-order ordinary differential equaion. In Section 3, we show that, for traveling wave solutions, (1.7) may be reduced to a first-order ordinary differential equation involving two arbitrary integration constants a and b. We show that there are four distinct periodic solutions corresponding to four different ranges of values of a and restricted ranges of values of b. A short conclusion is given in Section 4.

2. Solutions to a First-Order Ordinary Differential Equaion

This section is due to Vakhnenko and Parkes (see [19, Appendix  A]). For completeness and readability, we state it in the following.

Consider solutions to the following ordinary differential equation:

(φφξ)2=ε2f(φ), where

f(φ)=(φ-φ1)(φ-φ2)(φ3-φ)(φ4-φ), and φ1, φ2, φ3, φ4 are chosen to be real constants with φ1φ2φφ3φ4.

Following  we introduce ζ defined by

dξdζ=φε, so that (2.1) becomes

(φζ)2=f(φ).

Equation (2.4) has two possible forms of solution. The first form is found using result 254.00 in . Its parametric form is

φ=φ2-φ1nsn2(wm)1-nsn2(wm),ξ=1εp(wφ1+(φ2-φ1)(n;wm)), with w as the parameter, where

m=(φ3-φ2)(φ4-φ1)(φ4-φ2)(φ3-φ1),p=12(φ4-φ2)(φ3-φ1),w=pζ,n=φ3-φ2φ3-φ1. In (2.5) sn(wm) is a Jacobian elliptic function, where the notation is as used in [22, Chapter  16]. Π(n;wm) is the elliptic integral of the third kind and the notation is as used in [22, Section  17.2.15].

The solution to (2.1) is given in parametric form by (2.5) with w as the parameter. With respect to w, φ in (2.5) is periodic with period 2K(m), where K(m) is the complete elliptic integral of the first kind. It follows from (2.5) that the wavelength λ of the solution to (2.1) is

λ=2εp|φ1K(m)+(φ2-φ1)(nm)|, where Π(nm) is the complete elliptic integral of the third kind.

When φ3=φ4, m=1, (2.5) becomes

φ=φ2-φ1ntanh2w1-ntanh2w,ξ=1ε(wφ3p-2tanh-1(ntanhw)).

The second form of the solution to (2.4) is found using result 255.00 in . Its parametric form is

φ=φ3-φ4nsn2(wm)1-nsn2(wm),ξ=1εp(wφ4-(φ4-φ3)(n;wm)), where m,p,w are as in (2.6), and

n=φ3-φ2φ4-φ2.

The solution to (2.1) is given in parametric form by (2.10) with w as the parameter. The wavelength λ of the solution to (2.1) is

λ=2εp|φ4K(m)-(φ4-φ3)(nm)|.

When φ1=φ2,   m=1, (2.10) becomes

φ=φ3-φ4ntanh2w1-ntanh2w,ξ=1ε(wφ2p+2tanh-1(ntanhw)).

3. Solitary and Periodic Wave Solutions to (<xref ref-type="disp-formula" rid="EEq1.7">1.7</xref>)

Equation (1.7) can also be written in the form

ut+2uux-6uxuxx-2uuxxx=0. Let u=φ(ξ)+c with ξ=x-ct be a traveling wave solution to (3.1), then it follows that

-cφξ+2φφξ-6φξφξξ-2φφξξξ=0, where φξ is the derivative of function φ with respect to ξ.

Integrating (3.2) twice with respect to ξ yields

(φφξ)2=14(φ4-4c3φ3+aφ2+b), where a and b are two arbitrary integration constants.

Equation (3.3) is in the form of (2.1) with ε=1/2 and f(φ)=(φ4-(4c/3)φ3+aφ2+b). For convenience we define g(φ) and h(φ) by

f(φ)=φ2g(φ)+b,where  g(φ)=φ2-4c3φ+a,f(φ)=2φh(φ),where  h(φ)=2φ2-2cφ+a, and define φL, φR, bL, and bR by

φL=12(c-c2-2a),φR=12(c+c2-2a),bL=-φL2g(φL)=a24-12c2a+c46-16(c3-2ac)c2-2a,bR=-φL2g(φL)=a24-12c2a+c46+16(c3-2ac)c2-2a. Obviously, φL, φR are the roots of h(φ)=0.

Without loss of generality, we suppose the wave speed c>0. In the following, suppose that a<c2/2 and a0 for each value c>0, such that f(φ) has three distinct stationary points: φL, φR, 0 and comprise two minimums separated by a maximum. Under this assumption, (1.7) has periodic and solitary wave solutions that have different analytical forms depending on the values of a and b as follows.

a<0

In this case φL<0<φR and f(φL)>f(φR). For each value a<0 and 0<b<bL (a corresponding curve of f(φ) is shown in Figure 1(a)), there are periodic loop-like solutions to (3.3) given by (2.10) so that 0<m<1, and with wavelength given by (2.12). See Figure 2(a) for an example.

The curve of f(φ) for the wave speed c=2. (a) a=-40, b=200; (b) a=-40, b=226.0424; (c) a=1.5, b=-0.05; (d) a=1.5, b=0; (e) a=16/9, b=-0.1; (f) a=16/9, b=0; (g) a=17/9, b=-0.24; (h) a=17/9, b=-0.1842.

Periodic solutions to (3.3) with 0<m<1 and the wave speed c=2. (a) a=-40, b=200 so m=0.8978; (b) a=1.5, b=-0.05 so m=0.6893; (c) a=16/9, b=-0.1 so m=0.8254; (d) a=17/9, b=-0.24 so m=0.8412.

The case a<0 and b=bL (a corresponding curve of f(φ) is shown in Figure 1(b)) corresponds to the limit φ1=φ2=φL so that m=1, and then the solution is a loop-like solitary wave given by (2.13) with φ2φ<φR and

φ3=12c2-2a+c6-13c2+3c4-2a,φ4=12c2-2a+c6+13c2+3c4-2a. See Figure 3(a) for an example.

Solutions to (3.3) with m=1 and the wave speed c=2. (a) a=-40, b=226.0424; (b) a=1.5, b=0; (c) a=16/9, b=0; (d) a=17/9, b=-0.1842.

0<a<4c2/9

In this case 0<φL<φR and f(φR)<f(0). For each value 0<a<4c2/9 and bL<b<0 (a corresponding curve of f(φ) is shown in Figure 1(c)), there are periodic valley-like solutions to (3.3) given by (2.10) so that 0<m<1, and with wavelength given by (2.12). See Figure 2(b) for an example.

The case 0<a<4c2/9 and b=0 (a corresponding curve of f(φ) is shown in Figure 1(d)) corresponds to the limit φ1=φ2=0 so that m=1, and then the solution can be given by (2.13) with φ3 and φ4 given by the roots of g(φ)=0, namely,

φ3=2c3-4c29-a,φ4=2c3+4c29-a. In this case we obtain a weak solution, namely, the periodic downward-cusp wave

φ=φ(ξ-2jξm),(2j-1)ξm<ξ<(2j+1)ξm,j=0,±1,±2,, where

φ(ξ)=(φ3-φ4tanh2(ξ4))cosh2(ξ4),ξm=4tanh-1φ3φ4. See Figure 3(b) for an example.

a=4c2/9

In this case 0<φL<φR and f(φR)=f(0). For a=4c2/9 and each value bL<b<0 (a corresponding curve of f(φ) is shown in Figure 1(e)), there are periodic valley-like solutions to (3.3) given by (2.5) so that 0<m<1, and with wavelength given by (2.8). See Figure 2(c) for an example.

The case a=4c2/9 and b=0 (a corresponding curve of f(φ) is shown in Figure 1(f)) corresponds to the limit φ3=φ4=φR=2c/3 and φ1=φ2=0 so that m=1. In this case neither (2.9) nor (2.13) is appropriate. Instead we consider (3.3) with f(φ)=(1/4)φ2(φ-(2c/3))2 and note that the bound solution has 0<φ<2c/3. On integrating (3.3) and setting φ=0 at ξ=0 we obtain a weak solution

φ=-2c3exp(-12|ξ|)+2c3, that is, a single valley-like peaked solution with amplitude 2c/3. See Figure 3(c) for an example.

4c2/9<a<c2/2

In this case 0<φL<φR and f(φR)>f(0). For each value 4c2/9<a<c2/2 and bL<b<bR (a corresponding curve of f(φ) is shown in Figure 1(g)), there are periodic valley-like solutions to (3.3) given by (2.5) so that 0<m<1, and with wavelength given by (2.8). See Figure 2(d) for an example.

The case 4c2/9<a<c2/2 and b=bR (a corresponding curve of f(φ) is shown in Figure 1(h)) corresponds to the limit φ3=φ4=φR so that m=1, and then the solution is a velley-like solitary wave given by (2.10) with φL<φφ3 and

φ1=c6-12c2-2a-13c2-3cc2-2a,φ2=c6-12c2-2a+13c2-3cc2-2a. See Figure 3(d) for an example.

4. Conclusion

In this paper, we have found expressions for two types of traveling wave solutions to the osmosis K(2,2) equation, that is, the soliton and periodic wave solutions. These solutions depend, in effect, on two parameters a and m. For m=1, there are loop-like (a<0), peakon (a=4c2/9), and smooth (4c2/9<a<c2/2) soliton solutions. For m=1,0<a<4c2/9 or 0<m<1,  a<c2/2, and a0, there are periodic wave solutions.

Acknowledgments

The authors are deeply grateful to the anonymous referees for careful reading of this paper and constructive comments. J. Zhou acknowledges funding from the Startup Fund for Advanced Talents of Jiangsu University (no. 09JDG013), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (no. 09KJB110003) and Jiangsu Planned Projects for Postdoctoral Research Funds. L. Tian was partially supported by the National Natural Science Foundation of China (no. 10771088). X. Fan was supported by the Postdoctoral Foundation of China (no. 20080441071), the Postdoctoral Foundation of Jiangsu Province (no. 0802073c) and the High-level Talented Person Special Subsidizes of Jiangsu University (no. 08JDG013).

RosenauP.HymanJ. M.Compactons: solitons with finite wavelengthPhysical Review Letters1993705564567EID2-s2.0-1204425319910.1103/PhysRevLett.70.564WazwazA. M.Compactons and solitary patterns structures for variants of the KdV and the KP equationsApplied Mathematics and Computation200313913754MR194937510.1016/S0096-3003(02)00120-0ZBL1029.35200HeJ.-H.Homotopy perturbation method for bifurcation of nonlinear problemsInternational Journal of Nonlinear Sciences and Numerical Simulation200562207208EID2-s2.0-4243080805HeJ.-H.WuX.-H.Exp-function method for nonlinear wave equationsChaos, Solitons & Fractals2006303700708MR223869510.1016/j.chaos.2006.03.020ZBL1141.35448HeJ.-H.WuX.-H.Construction of solitary solution and compacton-like solution by variational iteration methodChaos, Solitons & Fractals2006291108113MR221334010.1016/j.chaos.2005.10.100ZBL1147.35338HeJ.-H.Some asymptotic methods for strongly nonlinear equationsInternational Journal of Modern Physics B2006201011411199MR2220565ZBL1102.34039XuL.Variational approach to solitons of nonlinear dispersive K(m,n) equationsChaos, Solitons & Fractals2008371137143MR239706310.1016/j.chaos.2006.08.016ZBL1143.35361WazwazA. M.General compactons solutions and solitary patterns solutions for modified nonlinear dispersive equations mK(n,n) in higher dimensional spacesMathematics and Computers in Simulation2002596519531MR191782210.1016/S0378-4754(01)00439-6ZBL0996.35065WazwazA. M.Compact and noncompact structures for a variant of KdV equation in higher dimensionsApplied Mathematics and Computation200213212945MR191870410.1016/S0096-3003(01)00173-4ZBL1031.35128ChenY.LiB.ZhangH. Q.New exact solutions for modified nonlinear dispersive equations mK(m,n) in higher dimensions spacesMathematics and Computers in Simulation2004645549559MR203950010.1016/j.matcom.2003.10.005ZBL1073.35052HeB.MengQ.RuiW.LongY.Bifurcations of travelling wave solutions for the mK(n,n) equationCommunications in Nonlinear Science and Numerical Simulation2008131021142123MR241758110.1016/j.cnsns.2007.06.006YanZ. Y.Modified nonlinearly dispersive mK(m,n,k) equations. I. New compacton solutions and solitary pattern solutionsComputer Physics Communications200315212533MR1976396YanZ. Y.Modified nonlinearly dispersive mK(m,n,k) equations. II. Jacobi elliptic function solutionsComputer Physics Communications20031531116MR1979506BiswasA.1-soliton solution of the K(m,n) equation with generalized evolutionPhysics Letters A20083722546014602MR242428610.1016/j.physleta.2008.05.002ZhuY. G.TongK.LuT. C.New exact solitary-wave solutions for the K(2,2,1) and K(3,3,1) equationsChaos, Solitons & Fractals200733414111416MR231892610.1016/j.chaos.2006.01.090ZBL1137.35427XuC.TianL.tianlx@ujs.edu.cnThe bifurcation and peakon for K(2,2) equation with osmosis dispersionChaos, Solitons & Fractals2009402893901EID2-s2.0-2134444727310.1016/j.chaos.2007.08.042ZhouJ.TianL.Soliton solution of the osmosis K(2,2) equationPhysics Letters A20083724162326234MR245644710.1016/j.physleta.2008.08.053VakhnenkoV. O.ParkesE. J.Explicit solutions of the Camassa-Holm equationChaos, Solitons & Fractals200526513091316MR214931610.1016/j.chaos.2005.03.011ZBL1072.35156VakhnenkoV. O.ParkesE. J.Periodic and solitary-wave solutions of the Degasperis-Procesi equationChaos, Solitons & Fractals200420510591073MR202937910.1016/j.chaos.2003.09.043ZBL1049.35162ParkesE. J.The stability of solutions of Vakhnenko's equationJournal of Physics A1993262264696475MR125304910.1088/0305-4470/26/22/040ZBL0809.35086ByrdP. F.FriedmanM. D.Handbook of Elliptic Integrals for Engineers and Scientists19712ndBerlin, GermanySpringer358Die Grundlehren der Mathematischen Wissenschaften, vol. 67MR0277773AbramowitzM.StegunI. A.Handbook of Mathematical Functions1972New York, NY, USADover