TRAVELING WAVE SOLUTIONS TO THE HYPERELASTIC ROD EQUATION†
TRAVELING WAVE SOLUTIONS TO THE HYPERELASTIC ROD EQUATION†
Journal of Applied Mathematics & Informatics. 2015. May, 33(3_4): 261-273
• Received : January 08, 2015
• Accepted : March 17, 2015
• Published : May 30, 2015
PDF
e-PUB
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
BYUNGSOO MOON

Abstract
We consider the hyperelastic rod equation describing nonlinear dispersive waves in compressible hyperelastic rods. We investigate the existence of certain traveling wave solutions to this equation. We also determine whether two other equations(the b -family equation and the modified Camassa-Holm equation) have our solution type. AMS Mathematics Subject Classification : 35Q35, 37K45.
Keywords
1. Introduction
This paper is concerned with the hyperelastic rod equation
PPT Slide
Lager Image
which was obtained as a model for nonlinear waves in cylindrical hyperelastic rods with u ( t, x ) representing the radial stretch relative to a pre-stressed state, and the physical parameter γ ranging from −29.4760 to 3.1474 [6 , 7 , 8] . From the mathematical view point, we regard γ as a real number. Among (1), there are two other important equations.
When γ = 1, it recovers the standard Camassa-Holm(CH) equation [3]
PPT Slide
Lager Image
arises as a model for the unidirectional propagation of shallow water waves over a flat bottom [3] , as well as water waves moving over an underlying shear flow [12] . The CH equation is completely integrable systems with a corresponding Lax pair formulation, a bi-Hamiltonian structure, and an infinite sequence of conservation laws [3] . The CH equation also admits peaked solitary waves or ”peakons” [3] : u ( t, x ) = ce −|xct| , c ≠ 0, which are smooth except at the crests, where they are continuous, but have a jump discontinuity in the first derivative. These peakons are shown to be stable [5] . Moreover, it admits the multi-peakon solutions (cf. [3] ).
If γ = 0, then (1) becomes the regularized long wave equation, a well-known equation for surface wave in a channel [1] . All solutions are global and the solitary waves are smooth. Despite having a Hamiltonian structure, the equation is not integrable and its solitary waves are not solitons [9] .
For general γ
PPT Slide
Lager Image
, mathematical properties of (1) have been studied further in many works. Dai and Huo [8] observed that (1) was formally shown to admit smooth, peaked, and cusped traveling waves. Subsequently, Lenells [13] used a suitable framework for weak solutions to classify all weak traveling waves of the hyperelastic rod equation (1). Constantin and Strauss [6] also investigated the stability of a class of solitary waves for the rod equation (1) on the line.
Note that if p ( x ) :=
PPT Slide
Lager Image
, x
PPT Slide
Lager Image
, then u =
PPT Slide
Lager Image
= p m , where m := u uxx and ∗ denotes the convolution product on
PPT Slide
Lager Image
, given by
PPT Slide
Lager Image
This formulation allows us to define a weak form of (1) as follows :
PPT Slide
Lager Image
In this paper, we are interested in a special traveling wave (i.e. function of x ct for some c > 0) and in fact that solution will also be a solitary wave solution (i.e. function of x ct , decays to zero as | x | → ∞). Our solution will be bounded in x
PPT Slide
Lager Image
for each fixed t > 0. In fact, it will be bounded by
PPT Slide
Lager Image
. The solution will be valid for any fixed value of c > 0. This will be a weak solution, namely u ( t, x ) will not be differentiable at the ”peak” point. At any other points the solution will have all the derivatives.
We are motivated by the paper [4] , where they exhibited a similar solution to the following equation
PPT Slide
Lager Image
where 0 ≤ β < 1, but their result is also applicable to
PPT Slide
Lager Image
for any α 2 > 0 and β ∈ [0, 1). By using transformation x cx , we have
PPT Slide
Lager Image
Choose c such that c =
PPT Slide
Lager Image
. Then we obtain
PPT Slide
Lager Image
Our case is then obtained by choosing α > 0 such that
PPT Slide
Lager Image
3). Our case is
PPT Slide
Lager Image
with β ∈ [0, 1). In our case, we can use result in ]undefined [4] to get a solution for our equation
PPT Slide
Lager Image
when γ ∈ (1, 3]. Solution in [4] covers our equation when γ ∈ (1, 3], but not when γ > 3. In our case, we also cover γ > 3. In fact, our solutionwill be valid when γ > 1.
This paper is organized as follows. In Section 2, we give the needed results to pursue our goal. In Section 3, we shall show the existence of traveling wave solutions for (2) under the condition γ > 1, by using an analogous analysis in [4] . Finally, in Section 4 we determine whether the b -family equation [11] and the modified Camassa-holm equation [10] have our solution type.
2. Preliminaries
In this section, we present the following basic technical Lemmas which play a key role to obtain our main results.
Lemma 2.1. Suppose γ > 1 and
PPT Slide
Lager Image
be a sequence of positive real numbers defined
PPT Slide
Lager Image
Define
PPT Slide
Lager Image
Then, there exist R > 0 such that φ ( x ) is well-defined on | x | < R, where R is the radius of the convergence for φ ( x ).
Proof . Let
PPT Slide
Lager Image
be a sequence of positive real numbers defined
PPT Slide
Lager Image
Set
PPT Slide
Lager Image
It is obviously that αn ξn , for n = 1, 2, .... Note that
PPT Slide
Lager Image
Solving the equation (5), we have
PPT Slide
Lager Image
This implies ϱ ( x ) is well-defined on
PPT Slide
Lager Image
Hence, φ ( x ) is well-defined, at least, on
PPT Slide
Lager Image
which means there is
PPT Slide
Lager Image
such that φ ( x ) is well-defined on | x | < R . This completes the proof of Lemma2.1. ΋
Lemma 2.2. Consider the following initial value probem ( IVP )
PPT Slide
Lager Image
The solution of (6) can be implicitly expressed as (11).
Proof . Notice that (6) is independent of z . By assuming that
PPT Slide
Lager Image
:= ᴪ( φ ), we reach
PPT Slide
Lager Image
Solving the first-order ordinary differential equation (7), we get
PPT Slide
Lager Image
By (8) and the transformation z x with x := ez used in (15), we obtain
PPT Slide
Lager Image
Using the transformation φ η with φ =
PPT Slide
Lager Image
we have
PPT Slide
Lager Image
Solving (10) to obtain an implicit analytic formula for function φ :
PPT Slide
Lager Image
for x > 0. This completes the proof of Lemma 2.2. ΋
Under the same assumption of Lemma 2.1, we have the following important result to get (39).
Lemma 2.3. Assume γ > 1. Let
PPT Slide
Lager Image
and φ ( x ) be defined in (3) and (4). Then φ ( x ) has a continuous extension to x = R such that φ ( R ) =
PPT Slide
Lager Image
, or lim x R φ ( x ) =
PPT Slide
Lager Image
.
Proof . By Lemma 2.1, we see tht φ ( x ) is well-defined on | x | < R . In order to show that φ ( x ) has a continuous extension such that φ ( R ) =
PPT Slide
Lager Image
, or lim x R φ ( x ) =
PPT Slide
Lager Image
, let us define
PPT Slide
Lager Image
Since
PPT Slide
Lager Image
we obtain
PPT Slide
Lager Image
Notice that for | x | < R
PPT Slide
Lager Image
By (12) and (13), we have
PPT Slide
Lager Image
Since our goal is to show lim x R φ ( x ) =
PPT Slide
Lager Image
, we may assume that x > 0. Using the transformation z x with x :=
PPT Slide
Lager Image
, we have
PPT Slide
Lager Image
Substituting (15) into (14), we obtain the ordinary differential equation
PPT Slide
Lager Image
which is independent of z . By Lemma 2.2, we obtain implicitly expression of solution for φ
PPT Slide
Lager Image
for x > 0. Let
PPT Slide
Lager Image
By direct calculation we have
PPT Slide
Lager Image
By the Implicit Function Theorem, we conclude that φ ( x ) is well-defined in | x | < R , where
PPT Slide
Lager Image
and
PPT Slide
Lager Image
Since αn ≥ 0, n = 1, 2, ..., it is easy to check that
PPT Slide
Lager Image
This completes the proof of Lemma 2.3. ΋
3. Traveling wave solutions
In this section, we discuss the existence of traveling wave solutions to the hyperelastic rod equation (2) with γ > 1.
Theorem 3.1. Let γ > 1. For every c > 0, the functions of the form
PPT Slide
Lager Image
where α n and R are defined in (3) and (20), is a solution to (2) in the weak sense .
Proof . We need to show that u ( t, x ) satisfies (2). Notice that
PPT Slide
Lager Image
PPT Slide
Lager Image
Using (23), (24), and p ( x ) =
PPT Slide
Lager Image
for x
PPT Slide
Lager Image
, we calculate from (2) that
PPT Slide
Lager Image
where
PPT Slide
Lager Image
with Bn =
PPT Slide
Lager Image
When x > ct , by using
PPT Slide
Lager Image
we write K 2 defined in (25) as K 2 = I 1 + I 2 + I 3 , where
PPT Slide
Lager Image
We directly compute I 1 as follows:
PPT Slide
Lager Image
In a similar manner,
PPT Slide
Lager Image
and
PPT Slide
Lager Image
Plugging (27)-(29) into (26), we deduce that for x > ct
PPT Slide
Lager Image
Using (25) and (30), we obtain for x > ct
PPT Slide
Lager Image
While for the case x ct , we split second term( K 2 ) of (25) into the following three parts:
PPT Slide
Lager Image
For J 1 , a direct computation gives rise to
PPT Slide
Lager Image
Similarly, one obtains
PPT Slide
Lager Image
and
PPT Slide
Lager Image
Using (25) and (33)-(35), we have for x ct
PPT Slide
Lager Image
PPT Slide
Lager Image
The traveling wave profile u(0.x) for c=2.
PPT Slide
Lager Image
The zoom in profile for the traveling wave u(0.x) for c=2.
On the other hand, using (23), we get
PPT Slide
Lager Image
As a result, (2) is equivalent to the recursive sequence of equation given by
PPT Slide
Lager Image
With the help of (3) and (22), we conclude that (38) holds if
PPT Slide
Lager Image
Now, we apply Lemma 2.1 and 2.3 to show that (39) holds. Let us define
PPT Slide
Lager Image
From the definition of θ ( x ) in Lemma 2.3, we see that Φ( x ) = θ′ ( x ) for | x | < R . By (12), we have θ′ ( x ) = 1 + ( γφ ( x ) − 1) φ′ ( x ) for | x | < R . By Lemma 2.3, we have lim x R φ ( x ) =
PPT Slide
Lager Image
and hence
PPT Slide
Lager Image
Therefore,
PPT Slide
Lager Image
Since αn ≥ 0 for n ≥ 1, one can easily to check that
PPT Slide
Lager Image
This completes the proof of Theorem 3.1. ΋
4. Remarks
In Sections 2-3 we have found explicit expressions for traveling wave solutions to (2) with γ > 1 that travel in the positive x -direction with speed c > 0. In the following, we discuss whether two other equations(the b -family equation and the modified Camassa-Holm equation) have our solution type. We expect our non-smooth traveling wave solutions satisfy to the above both equations. Unfortunately, their equations does not have our traveling wave solution type.
Consider the following two other equations which are the b -family equation [11]
PPT Slide
Lager Image
and the modified Camassa-Holm equation [10]
PPT Slide
Lager Image
Moreover, (40) and (41) are equivalent to the following nonlocal form:
PPT Slide
Lager Image
and
PPT Slide
Lager Image
We just apply our solution to the b -family equation. Let us substitute our solution u ( t, x ) =
PPT Slide
Lager Image
in (22) to (42). In a similar arguments, substituting (23) and (24) into (42) gives us
PPT Slide
Lager Image
Remark 4.1. Since c > 0, (44) implies that βn = 0 for n ≥ 2 and β 1 = c or 0. Therefore, we conclude that (42) have u ( t, x ) = ce −|xct| or u ( t, x ) = 0. In this case, our solution recovers the single peakon solution.
Next we also apply our solution to the modified Camass-Holm equation. In a similar manner, we can plug our solution u ( t, x ) =
PPT Slide
Lager Image
into (43). Then we obtain
PPT Slide
Lager Image
Remark 4.2. Since c > 0, β 1 =
PPT Slide
Lager Image
or 0, β 2 = 0 and βn = 0 for n ≥ 3. Hence (45) have a solution u ( t, x ) =
PPT Slide
Lager Image
e −|xct| or u ( t, x ) = 0. In this case, our solution also recovers single peakon solution which was obtained in [10] .
In the following remark, we explain our solutions fit into the classification previously studied traveling wave solutions of (1).
Remark 4.3. In [13], Lenells study our own equation
PPT Slide
Lager Image
and categorized all weak traveling wave solutions. Our solution corresponds to his category (f) of Theorem 1 in [13] . He proves the existence of a cuspon solution with decay, which corresponds to our solution such that
PPT Slide
Lager Image
where αn and R are defined in (3) and (20).
BIO
Byungsoo Moon received a M.A. from the State University of New York at Buffalo and Ph.D in the University of Texas at Arlington. He is a postdoctoral fellowship of the Pohang Mathematics Institute (PMI) at POSTECH. His current research interests are the study of a certain class of nonlinear partial differential equations referred to as integrable systems.
Pohang Mathematics Institute (PMI), POSTECH, Pohang, Kyungbuk 790-784, Korea
e-mail: bsmoon@postech.ac.kr
References
Benjamin B.T. , Bona L.J. , Mahony J.J. (1972) Model equations for long waves in nonlinear dispersive systems Philos. Trans. Roy. Soc. London Ser. A 272 47 - 48    DOI : 10.1098/rsta.1972.0032
Bougoffa L. , Al-Mazmumy M. (2013) Series solutions to initial-Neumann boundary value problems for parabolic and hyperbolic equations J. Appl. Math. Inform 31 87 - 97    DOI : 10.14317/jami.2013.087
Camassa R. , Holm D. (1993) An integrable shallow water equation with peaked solitons Phys.Rev. Lett. 71 1661 - 1664    DOI : 10.1103/PhysRevLett.71.1661
Cao C. , Holm D. , Titi E. (2004) Traveling wave solutions for a class of one-dimensional nonlinear shallow water wave models J. Dynam. Differential Equations 16 167 - 178    DOI : 10.1023/B:JODY.0000041284.26400.d0
Constantin A. , Strauss W. (2000) Stability of peakons Comm. Pure Appl. Math. 53 603 - 610    DOI : 10.1002/(SICI)1097-0312(200005)53:5<603::AID-CPA3>3.0.CO;2-L
Constantin A. , Strauss W. (2000) Stability of a class of solitary waves in compressible elastic rods Phys. Lett. A 270 140 - 148    DOI : 10.1016/S0375-9601(00)00255-3
Dai H.H. (1998) Model equations for nonlinear dispersive waves in a compressible Mooney-Rivlin rod Acta Mech. 127 193 - 207    DOI : 10.1007/BF01170373
Dai H.H. , Huo Y. (2000) Solitary shock waves and other travelling waves in a general compressible hyperelastic rod R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. Sci. 456 331 - 363    DOI : 10.1098/rspa.2000.0520
Dodd R.K. , Eilbeck J.C. , Gibbon J.D. , Morris H.C. 1982 Solitons and nonlinear wave equations Academic Press London–New York
Gui G. , Liu Y. , Olver P. J. , Qu C. (2013) Wave-breaking and peakons for a modified Camassa-Holm equation Comm. Math. Phys. 319 731 - 759    DOI : 10.1007/s00220-012-1566-0
Holm D.D. , Staley M.F. (2003) Staley, Nonlinear balances and exchange of stablity in dynamics of solitons, peakons, ramp/cliffs and leftons in a 1 + 1 nonlinear evolutionary PDE Phy.Lett. A 308 437 - 444    DOI : 10.1016/S0375-9601(03)00114-2
Johnson R.S. (2003) The Camassa-Holm equation for water waves moving over a shear flow Fluid Dynam. Res. 33 97 - 111    DOI : 10.1016/S0169-5983(03)00036-4
Lenells J. (2006) Traveling waves in compressible elastic rods Discrete Contin. Dyn. Syst. Ser.B 6 151 - 167    DOI : 10.3934/dcdsb.2006.6.151
Wazwaz A.M. (2011) Completely integrable coupled potential KdV equations J. Appl. Math.Inform. 29 847 - 858
Wazwaz A.M. (2012) N-soliton solutions for the sine-Gordon equation of different dimensions J. Appl. Math. Inform. 30 925 - 934