Integrablepropertiesforavariable-coefficientBoussinesqequationfromweaklynonlineardynamicswithsymboliccomputation.doc

精品论文integrable properties for a variable-coefficient boussinesq equation from weakly nonlinear dynamics with symbolic computationya-xing zhang1 , tao xu1, juan li1, chun-yi zhang2, 3, hai-qiang zhang1 and bo tian11. school of science, p. o. box 122, beijing university of posts and telecommunications, beijing 100876, china2. ministry-of-education key laboratory of fluid mechanics and national laboratory for computational fluid dynamics, beijing university of aeronautics and astronautics, beijing 100083, china3. meteorology center of air force command post, changchun 130051, chinaabstractdescribing the weakly nonlinear dynamics of long waves embedded in marginally stable shear flows that vary in the streamwise direction, a variable-coefficient boussinesq equation is investigated in this paper. with symbolic computation, such a model is transformed into its constant-coefficient counterpart under certain constraints on the coefficient functions. by virtue of the obtained transformation, some integrable properties for this equation are derived, such as the auto-backlund transformation, nonlinear superposition formula, lax pair and dar- boux transformation. in addition, some soliton-like solutions are obtained from the integrable properties and the relevant physical applications are also pointed out.pacs numbers: 02.30.ik; 05.45.yv; 02.30.jr; 52.35.mw; 02.70.wz mailing address: 4精品论文1. introductionas an important nonlinear model in hydrodynamics, the boussinesq equation is widely employed to describe some physical and applied mathematic problems including the nonlinear kelvin-helmholtz instability 1, 2, nonlinear lattices 3, 4, nonlinear baroclinic instability 5,6, vibrations in a nonlinear string 7 and ion sound waves in plasmas 8. in addition, such an equation is also applicable to modelling the propagation of long waves in shallow water under gravity and in one-dimensional nonlinear lattices 10. applications of the constant- coefficient boussinesq equation can be seen from the following examples: (1) the originalboussinesq equation 11tt xx2xxxxxxu m u+ 3 (u2) + n u= 0 ,(1)governs the propagation of long waves in shallow water in a uniform environment, wherem and n are both constants, while x and t represent a “long” space and a “slow” time,respectively. (2) the “good” boussinesq equation 12 is written asutt = uxx uxxxx + (u2 )xx ,x r ,(2)which has been investigated in detail through the analytic and numerical methods 13. it has been pointed out in refs. 13 that eqn. (2) exhibits some interesting behavior: the small-amplitude solitons keep their shapes and velocities unchangeable in the interaction process, while the large-amplitude ones change into the so-called antisolitons after collision.in reality, the media may be inhomogeneous and the boundaries may be nonuniform, e.g., in plasmas 14, 15, superconductors 16, optical-fiber communications 17, 18, blood vessels and bose-einstein condensates 19. therefore, the variable-coefficient nonlinear evo- lution equations (nlees) are considered to be more realistic than their constant-coefficient counterparts in describing a large variety of real situations. for instance, in terrain-followingcoordinates, the second order boussinesq equation,%$ 2m () 2 m 0() 2f m () ftt + 3 m () f tt m ()f +ft + $tm 2()ft f = 0 ,(3)is a useful model in the presence of rapidly varying topographies 20, where the prime means differentiation with respect to , $ is the nonlinearity parameter, % is the dispersion parameter and m is the variable free surface coefficient. this model governs the velocitypotential along the bottom topography and when a topography is absent, it reduces to the standard boussinesq equation.on the other hand, the coefficient functions related to x and/or t bring about a greatamount of complicated calculations which are unmanageable manually, although the variable- coefficient nlees are of current importance in many fields of physics and engineering sci- ences 21, 22. fortunately, however, symbolic computation 14, 15, 23, as a new branch of artificial intelligence, makes it exercisable to deal with those nonlinear systems with variable coefficients.in this paper, we would like to investigate the following variable-coefficient boussinesq equation with the space-dependent coefficients 24utt + uxxx + + f (x) ux + u ux + g(x) ux = 0 ,(4)which can describe the long waves embedded in marginally stable shear flows varying in the streamwise direction and can govern new nonmodal or modal instabilities, where , and are constants defined in ref. 5; f (x) and g(x), related to the streamwise varying background flow, are both the analytic functions; f (x) uxx and g(x) ux contribute to theemergence of new (linear) modal or nonmodal instabilities. when the background flow is streamwise uniform, i.e., f (x) = 0 and g(x) = 0, eqn. (4) reduces to the constant-coefficient boussinesq equation, which possesses a lot of good properties such as the backlund transformation 25, 26 and lax pair 27, 28.as shown in refs. 18, 19 and references therein, it is a straightforward and effective method for investigating the variable-coefficient nlees to transform them into some known equations. thus, we organize our paper as follows. in section 2, with the aid of symbolic computation, the transformation from eqn. (4) to a constant-coefficient boussinesq equation is constructed under certain constraints. in section 3, taking advantage of the obtained transformation, a series of integrable properties are figured out when the relevant constraints are satisfied, such as the auto-backlund transformation, nonlinear superposition formula, lax pair and darboux transformation. in section 4, several soliton-like solutions, differing from the existing ones, are presented through the above properties and possible applications are also discussed. section 5 is the conclusions for this paper.2. transformation from eqn. (4) to a constant-coefficient boussinesq equationin this section, we will determine the conditions for eqn. (4) to be transformable into the following constant-coefficient boussinesq equationu u(u t tx x x x22)x x ux x = 0 .(5)similar to refs. 18, 19, the transformation from eqn. (4) to eqn. (5) is assumed to be of the formu = b(x, t) + a(t) u x (x), t (t) with x (x) = a x + b ,(6)where b(x, t), a(t), x (x) and t (t) are all analytic functions to be determined, while a andb are both constants.by virtue of symbolic computation, eqn. (4) with the substitution of format (6) becomesa(t) t 0(t)2 ut t + a4 a(t) ux x x x + a2 a(t) + a2 a(t) b + a2 a(t) f (x) ux xx+a2 a(t)2 u ux x + a2 a(t)2 u 2+ a a(t) g(x) + a a(t) f 0(x) + 2 a a(t) bx ux+ 2 a0(t) t 0(t) + a(t) t 00(t) ut + a(t) g0(x) + a00(t) + a(t) bxx u + b g0(x)x+btt + g(x) bx + f 0(x) bx + b2 + bxx + b bxx + f (x) bxx + bxxxx = 0 ,(7)where the prime notes differentiation with respect to independent variable, which is eqn. (5),provided thata4 + t 0(t)2 = 0 ,(8)(a2 + a2 b + a2 f (x) + t 0(x)2 = 0 ,(9)a2 a(t) + t 0(t)2 = 0 ,(10)a g(x) + a f 0(x) + 2 a bx = 0 ,(11)2 a0(t) t 0(t) + a(t) t 00(t) = 0 ,(12)a(t) g0(x) + a00(t) + a(t) bxx = 0 ,(13)xbtt + g(x) bx + f 0(x) bx + b2 + bxx + f (x) bxx+b g0(x) + bxx + bxxxx = 0 ,(14)from which we can geta(t) = a ,(15)t (t) = c t + ,(16)1b(x, t) = 2 where c = 0, a and are three constants.zf (x) +g(x) dx,(17)subsequently, with the help of symbolic computation, substituting eqn. (17) into eqn. (9)yields the constraints on the variable coefficients f (x) and g(x), as follows,zf (x) =g(x) dx 2 c2 a2 2 ,(18)f 0(x)2 + f 00(x) + f (x) f 00(x) + f (4) (x) = 0 .(19)integrating eqn. (19) with respect to x gives rise tof (x)2 f 00(x) + f (x) +2= 0 ,(20)which is an elliptic equation with the following exact analytic solutions 29a1 a2 secha 2 h a1 i2 (x x0)2 h a1 i,(a1 0, a2 f (x) 0, a2 f (x) 0)(21)where a1 = a2 sec3 and a2 = 1 .2 (x x0),(a1 0)therefore, the transformation from eqn. (4) to eqn. (5) turns out to beu = b(x) + a u x (x), t (t) ,(22)精品论文5withx (x) = a x + b ,(23)t (t) = c t + ,(24)1b(x) = 2 2zf (x) +g(x) dx,(25)where a2 = , a = and c2 = , while b and are both arbitrary constants, and 2 f (x) and g(x) obey constraints (18) and (19).it is worth noting that eqn. (5) is solvable by the inverse scattering method 7 and possesses many remarkable properties such as the backlund transformation 25, 26, painleve property 26, lax representation 27, 28, an infinite number of conservation laws 30 and darboux transformation 27, 31. accordingly, when constraints (18) and (19) are satisfied, eqn. (4) could be dealt with via soliton theory based on eqn. (5) (see section 3).3. integrable propertiesgenerally speaking, the variable-coefficient nlees are not integrable unless the coeffi- cient functions satisfy certain constraints. in view of the current importance of variable- coefficient nlees, many efforts have been devoted to studying their integrable properties by transforming them into the integrable constant-coefficient counterparts 32. in what fol- lows, via the obtained transformation in section 2, we will further investigate the integrable properties of eqn. (4) with constraints (18) and (19).3.1. auto-backlund transformation for eqn. (4)to eqn. (4), the auto-backlund transformation exists along with each family of those solutions, under the corresponding set of the aforementioned constraints on f (x) and g(x).by virtue of transformation (22), we will construct the auto-backlund transformation foreqn. (4).without loss of generality, we shall assume that = 1, = 6 and = 1 in eqn. (5) sincethe following equationis equivalent to eqn. (5).ut t ux x ux x x x 3 (u 2)x x = 0 ,(26)as we know, the auto-backlund transformation for eqn. (26) has been worked out inref. 29. taking u = wx , eqn. (26) turns out to bewt t wx x 6 wx wx x wx x x x = 0 ,(27)which has the auto-backlund transformation as belowwt + w0, t = + 2 (w w0) + (w w0)2 + (w w0 )3+ 3 + (w w0) (wx w0, x ) + (wx x w0, x x ) ,(28)1wt w0,t = + 3 (w w0) (wx w0,x ) + (wx x + w0,x x ) ,(29)精品论文3where w and w0 are two different solutions of eqn. (27), while 2 = 1 , and are botharbitrary constants.utilizing transformation (22) with eqns. (28) and (29), where u = wx and u = wx ,gives rise to the following auto-backlund transformation for eqn. (4),6a c 2a c (w w0)a2 c (w w0)3wt + w0,t =+ c a (w w0 ) +a +a2 c +a+ 3 c (w w0 ) a(wx w0,x ) + c (wxx w0,xx )2a2,(30)wt w0, t =c a a 3 c (w w0 )(wx w0,x) +c (2 b0 + wxx + w0,xx )a2 ,(31)where a, c, , a, and have been defined as above, while w(x, t) satisfieswtt + wxxxx + + f (x) wxx + wx wxx + g(x) wx = 0 .(32)3.2. nonlinear superposition formula for eqn. (4)furthermore, on the ground of the obtained auto-backlund transformation, we will get the nonlinear superposition formula for eqn. (4), which can be used to construct the multi- soliton-like solutions.via the seed solution w0(x, t) for eqn. (32), when 1 = 1 and 2 = 2, the solu-tions w1(x, t), w2(x, t) for eqn. (32) are respectively obtained. consequently, the nonlinearsuperposition formula for eqn. (32) is expressed in the formzw3 = 21b(x) dx w0 + a a 21 a w1+ a w2z( a2 + 2 a (w2,x w1,x)a )+ (a 1 a 2 + a w2 a w1)w1 + w2 2b(x) dx,(33)where is an arbitrary constant, from which, as u = wx, the multi-soliton-like solutions foreqn. (4) can be explicitly derived.3.3. lax pairas we have mentioned above, under constraints (18) and (19), eqn. (4) has been con- verted into the constant-coefficient one by virtue of transformation (22). actually, when constraints (18) and (19) hold, eqn. (4) can also be converted into the following boussinesq equationut t ux x + (u 2)x +13 ux x x x = 0 ,(34)2through transformation (22) with a = 2 , a2 = and c2 = , while b and are botharbitrary constants.3 3 精品论文it is known that the lax pair of eqn. (34) is9 3l = v +u 3 3+,(35)24 x2 x 3withm = u x 2,(36)33vx = 4 ux x 4 ut ,(37)where lis the operator for the spectral problem and mis the operator governing theassociated time evolution of the eigenfunctions 27. under constraints (18) and (19), we can directly obtain the lax pair for eqn. (4), i.e.,l = ,(38) t = m ,(39)where t = 0, meanwhile, and (x, t) denote the eigenvalue and eigenfunction respectively, while l and m are presented in the forms3 3 (u b) 13l = v 4 a x +2 a a+ x a3 x3 ,(40)c (u b)c2withm = a a2 x2 ,(41)33 avx = 4 a a (uxx b00) 4 a c ut .(42)then it is easy to prove that l and m satisfy the lax equationlt m l + l m = 0 ,(43) if and only if u(x, t) is a solution of eqn. (4), which suggests that eqn. (4) is integrable under constraints (18) and (19).3.4. darboux transformation for eqn. (4)in virtue of the lax representation for eqn. (4), we can further derive its darboux transformation in this subsection.let (x, t) be a solution of the lax pair, i.e., (38) and (39), and h(x, t) be a particularsolution with = 0. then the darboux transformation for this lax pair ishx = x 2 au = u +,(44)h hx ,(45)a2hx33 h32hxhxxhxxx x v = v 2 a a (b0 ux) + a3h3 h2 +h.(46)精品论文through direct calculation with the aid of symbolic computation, it can be proved that the identities of eqns. (44-46) satisfy the lax pairl = ,(47)t = m ,(48)where l and mare obtained from l and m by replacing u(x, t) and v(x, t) with u(x, t) andv(x, t) respectively. meanwhile, we can also verify that the relation between u(x, t) and v(x, t)coincides with the relation between u(x, t) and v(x, t) under the above transformations,namely33 avx = 4 a a (uxx b00) 4 a c ut .(49)therefore, transformations (44-46) constitute the darboux transformation for eqn. (4) andu(x, t) is a new solution of eqn. (4).as we have shown above, some integrable properties have been presented. additionally, based on transformation (22), some other integrable properties for eqn. (4) could also be obtained, including the hirota bilinear representation, similarity reductions to the ordinary differential equations of painleve type, backlund transformation in bilinear form, an infinite number of independent conservation laws and painleve property.4. soliton-like solutions and possible applicationsin this section, by virtue of the obtained properties, we will construct some exact analytic solutions for eqn. (4), which differ from those obtained in ref. 24.utilizing expressions (30) and (31) by taking w0 = r b(x) dx as the seed solution, we canget the following one-soliton-like solution for eqn. (4),u = b(x) + 2 a 2 sech2 (a x c r t + b r ) ,(50)witha2 = 6 , a =2, c2 = ,r = , =433 , = 3,(51)where b(x) has been defined in expression (17).on the other hand, we can also get the soliton-like solutions for eqn. (4) from the darboux transformation. it is easy to prove that b(x) is the solution for eqn. (4), so introducing u = b(x) into the lax pair of eqn. (4), we have13a3 xxx 4 a x = ,(52)ct = a2 xx .(53)the position of is responsible for different types of soliton-like solutions. consider thefollowing cases:4case 4.1 | = 1in this case, we haveau1 = b + 22log a1 e1 + a2 e2 + a3 e3 ,(54)a2where i = ki x cik2 t, aia2 x2are arbitrary constants, ki(i = 1, 2, 3) are the roots of thealgebraic equation3 a2 k 4 k3 + 4 a3 = 0 ,(55)while a, a, c and b have been defi