TSWJ The Scientific World Journal 1537-744X 2356-6140 Hindawi Publishing Corporation 716740 10.1155/2014/716740 716740 Research Article Global Existence and Energy Decay Rates for a Kirchhoff-Type Wave Equation with Nonlinear Dissipation http://orcid.org/0000-0002-7929-1343 Kim Daewook 1 Kim Dojin 2 Hong Keum-Shik 3 Jung Il Hyo 1 Baleanu D. Jafari H. 1 Department of Mathematics Pusan National University 30 Jangjeon-dong Geumjeong-gu Busan 609-735 Republic of Korea pusan.ac.kr 2 Department of Mathematics Oregon State University Corvallis, OR 97331 USA oregonstate.edu 3 Department of Cogno-Mechatronics Engineering and School of Mechanical Engineering Pusan National University Busan 609-735 Republic of Korea pusan.ac.kr 2014 742014 2014 23 01 2014 23 02 2014 7 4 2014 2014 Copyright © 2014 Daewook Kim et al. This 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.

The first objective of this paper is to prove the existence and uniqueness of global solutions for a Kirchhoff-type wave equation with nonlinear dissipation of the form Ku+MA1/2u2Au+gu=0 under suitable assumptions on K,A,M(·), and g(·). Next, we derive decay estimates of the energy under some growth conditions on the nonlinear dissipation g. Lastly, numerical simulations in order to verify the analytical results are given.

1. Introduction

A mathematical model for the transverse deflection of an elastic string of length L>0 whose ends are held a fixed distance apart is written in the form of the hyperbolic equation (1)2u(x,t)t2-(α+β0L|u(x,t)x|2dx)2u(x,t)x2=0, which was proposed by Kirchhoff , where u(x,t) is the deflection of the point x of the string at the time t and α>0, β are constants. Kirchhoff first introduced (1) in the study of the oscillations of stretched strings and plates, so that (1) is called the wave equation of Kirchhoff type. The Kirchhoff-type model also appeared in scientific research for beam or plate . Such nonlinear Kirchhoff model gives one way to describe the dynamics of an axially moving string. In recent years, axially moving string-like continua such as wires, belts, chains, and band saws have been the subject of study of researchers .

The mathematical aspects of the natural generalization of the model (1) in Ωn: (2)u′′-M(Ω|u|2dx)Δu+g(u)=0,(3)u(0)=u0,u(0)=u1, under some assumptions on M(·), g(·), have been studied, using different methods, by many authors [6, 8, 1522].

When g(·)=0 and n=1, the problem (2)-(3) was studied by Dickey  and Bernstein  who considered analytic functions as the initial data (see also Yamada  and Ebihara et al. ). In case when g(·)=0 and n1, Pohožaev  obtained the existence and uniqueness of global solutions for the problem (2)-(3). Lions  also formulated Pohožaev's results in an abstract context and obtained better results.

Equation (2) with linear dissipative term, that is, g(u)=δu(δ>0), was investigated by Mizumachi , Nishihara and Yamada , Park et al. , and Jung and Choi . In fact, they studied the existence, uniqueness, and the energy decay rates of solutions for the problem (2)-(3). On the other hand, related works to a Kirchhoff-type equation with Ku′′ instead of u′′ can be found in Levine . Jung and Lee  got the result for a Kirchhoff-type equation with strong dissipative term. But they studied a simple form with the coefficient M(·)1. In case of the equation concerning nonlinear Kirchhoff-type coefficient, recently, Kim et al. , Ghisi and Gobbino , and Aassila and Kaya  have studied existence and energy decay rates of global (or local) solutions for the equation. By giving some suitable smallness conditions on the sizes of the initial data, they assured global existence and energy decay rates for the solutions.

In this paper, we study the existence, uniqueness, and the decay estimates of the energy for a class of Kirchhoff-type wave equations in a Hilbert space H: (4)Ku′′+M(|A1/2u|2)Au+g(u)=0in  H,u(0)=u0,(Ku)(0)=K1/2u1, where K and A are linear operators in H and M(·)C1[0,). For global existence of this problem, we give some suitable smallness conditions. So, the main contribution of these results is to consider a general model which contains the concrete model (2)-(3) and to improve the results of Kouémou-Patcheu  and Jung and Choi . Moreover, as an application, we give some simulation results about solution's shapes and the algebraic decay rate for a Kirchhoff-type wave equation with nonlinear dissipation.

The method applied in this paper is based on the multipliers technique , Galerkin's approximate method, and some integral inequalities due to Haraux .

This paper is organized as follows. In Section 2, we recall the notation, hypotheses, and some necessary preliminaries and prove the existence and uniqueness of global solutions for the system (4) by employing Feado-Galerkin's techniques under suitable smallness condition. In Section 3, we derive the energy decay rates by using the multiplier technique under suitable growth conditions on g. Finally, in Section 4, we give an example and its numerical simulations to illustrate our results.

2. Preliminaries and Existence

Let Ω be a bounded open domain in n having a smooth boundary Γ and H=L2(Ω) with inner product and norm denoted by (·,·) and |·|, respectively. Let K be a linear, positive, and self-adjoint operator on H; that is, there is a constant c>0 such that (5)(Ku,u)c|u|2,uH. Let A be a linear, self-adjoint, and positive operator in H, with domain V:=D(A) dense in H, KA=AK on D(A)D(K), and the graph norm denoted by ·. We assume that the imbedding VH is compact. Identifying H and its dual H, it follows that VHV, where V is the dual of V. Let ·,·V,V denote the duality pairing between V and V and W:=D(A1/2).

Throughout the paper we will make the following assumptions.

M(s) is a C1[0,) real function and M(s)0. Furthermore, there exist some positive constants β and γ0 such that M(s)β>0 for all s0 and |M(s)s|/M(s)γ0.

g: is a nondecreasing continuous function such that g(0)=0 and there is a constant k>0 and q1 such that (6)|g(x)|k(1+|x|q)x. And (g(u),Au)0 for all uD(A)D(A1/2). Note that the last assumption of (G) makes sense. In fact, when A=-Δ and g(u)=|u|αu,α1, we can easily show that (g(u),Au)0 for all uD(A)D(A1/2).

M(s)>M(s)|g(x)|, s[0,), x.

VLq+1(Ω) for some q1.

Let M¯(t) and E(t) be defined as follows: (7)M¯(t)=0tM(s)ds(8)E(t)=12[|K1/2u|2+M¯(|A1/2u|2)].

And also let us consider the functions (9)P(t):=|K1/2u(t)|2M(|A1/2u|2)+|A1/2u(t)|2,Q(t):=|K1/2A1/2u(t)|2M(|A1/2u|2)+|Au(t)|2,G(t):=|K1/2u(t)|M(|A1/2u|2).

Theorem 1.

Let the initial conditions (u0,u1)W×L2q(Ω) satisfy the smallness assumption (10)ML([0,P(0)])B(u0,u1)Q(0)<14, where B(u0,u1)=max{|K1/2u1|/M(|A1/2u0|2), M(|A1/2u0|2)/((M(|A1/2u0|2))  −  g(u1)M(|A1/2u0|2)))Q(0)}. Then there is a unique function uL(0,T;W)W1,(0,T;V)W2,(0,T;H) such that, for any T>0, (11)Ku+M(|A1/2u|2)Au+g(u)=0hhhhhhhhhhhhiinL(q+1)/q(0,T;V),(12)u(0)=u0,(Ku)(0)=K1/2u1.

Proof.

Assume that, for simplicity, V is separable; then there is a sequence (ej)j1 consisting of eigenfunctions of the operator A corresponding to positive real eigenvalues μj tending to + so that Aej=μjej, j1.

Let us denote by Vm the linear hull of e1,e2,,em. Note that (ej)j1 is a basis of H, V, and W and hence it is dense in H, V, and   W.

Approximate Solutions. We search for a function um(t)=j=1mgjm(t)ej such that, for any vVm, um(t) satisfies the approximate equation (13)(Kum′′(t)+M(|A1/2um|2)Aum+g(um),v)=0 and the initial conditions as the projections of u0 and u1 over Vm satisfy (14)um(0)=u0m=j=1m(u0,ej)eju0in  W(15)(Kum)(0)=K1/2u1m=j=1m(u1,ej)ejK1/2u1in  L2q(Ω).

For v=ej, j=1,2,m, the system (13)–(15) of ordinary differential equations of variable t has a solution um(t) in an interval [0,tm).

Now we obtain a priori estimates for the solution um(t) and it can also be extended to [0,T) for all T>0.

A Priori Estimate I. Let us consider v=um in (13). Using (7), we have (16)ddt(|K1/2um(t)|2+M¯(|A1/2um(t)|2))+2(g(um(t)),um(t))=0. Integrating (16) over (0,t), ttm, and using (8), we have (17)2E(0)=[|K1/2um(t)|2+M¯(|A1/2um(t)|2)]+20t(g(um(s)),um(s))ds. Using (5) and (7), we deduce that (18)2E(0)|K1/2um(t)|2+β|A1/2um(t)|22E(0)+20tΩum(s)g(um(s))dxds, where the left-hand side is constant independent of m and t. Thus estimation (18) yields, for any 0<T<, (19)um  bounded  in  L(0,T;H),(20)K1/2umboundedin  L(0,T;H),(21)A1/2um  bounded  in  L(0,T;H),(22)umg(um)  bounded  in  L1([0,T]×Ω).

Now we show that um(t) can be extended to [0,). We need the following smallness assumption: (23)ML([0,P(0)])×max{|K1/2um1|M(|A1/2um0|2),(M(|A1/2um0|2))hhhhhhhh×((M(|A1/2um0|2))(M(|A1/2um0|2))-g(um1)hhhhhhhhhhh×M(|A1/2um0|2)(M(|A1/2um0|2)))-1Q(0)|K1/2um1|M(|A1/2um0|2),}×Q(0)<14, where P(0)=(|K1/2um1|2/M(|A1/2um0|2))+|A1/2um0|2, Q(0)=(|K1/2A1/2um1|2  /  M(|A1/2um0|2))+|Aum0|2.

Let [0,T*) be the maximal interval where the solution exists. Set Z(t):=M(|A1/2um(t)|2) and (24)T:=sup{τ[0,T*)|Z(t)Z(t)|12,Z(t)>0,t[0,τ)}.

With simple computations it follows that (25)P(t)=-1Z(t)(2(g(um(t)),um(t))+Z(t)Z(t)|um(t)|2)P(t)0,(26)Q(t)=-1Z(t)(2(g(um(t)),Aum(t))+Z(t)Z(t)|A1/2um(t)|2)0,(27)(G2)(t)-G(t){2(Z(t)Z(t)-|g(um(t))|)G(t)-2|Aum(t)|}, for all t[0,T).

Next, we show that T=T*. Let us assume by contradiction that T<T*. Since |Z(t)|(1/2)Z(t) in [0,T), we have that (28)0<Z(0)e-T/2Z(T)Z(0)eT/2. Since Z(t) and Z(t) are continuous functions, by the maximality of T we have that necessarily (29)|Z(t)Z(t)|=12. From (88) and (89) it follows that P and Q are nonincreasing functions; hence (30)|A1/2um(t)|2P(t)P(0),|Aum(t)|2Q(t)Q(0). Moreover by Lemma 3.1 in  we have that (31)G(t)max{G(0),Z(0)Z(0)-g(um1)Z(0)Q(0)},  hhhhhhhhhhhhhhhhhhhhhhhhhhhhht[0,T]. By (91)–(31), and the smallness assumption (23), we have that (32)|Z(T)Z(T)|=|2M(|A1/2um(t)|2)(um(T),Aum(T))Z(T)|2max0rP(0)|M(r)||um(T)|Z(T)|Aum(T)|2max0rP(0)|M(r)|×max{G(0),Z(0)Z(0)-g(um1)Z(0)Q(0)}×Q(0)<12.

This contradicts (29). Therefore it follows that um(t) can be extended to [0,T) for any T(0,).

Furthermore, putting v=Aum in (13), we get (33)(Kum′′,Aum)M(|A1/2um|2)+(Aum,Aum)+(g(um),Aum)M(|A1/2um|2)=0. From this we obtain (34)12ddt((Kum,Aum)M(|A1/2um|2)+|Aum|2)+(g(um),Aum)M(|A1/2um|2)=-(Kum,Aum)M(|A1/2um|2)(A1/2um,A1/2um){M(|A1/2um|2)}2. Integrating (34) over (0,t) and taking into account assumptions (M) and (G), and applying Gronwall’s inequality, we obtain (35)Aum  bounded  in  L(0,T;H).

From (6) and (22), it follows that (36)g(um)  bounded  in  L(q+1)/q([0,T]×Ω).

A Priori Estimate II. Taking v=um′′(t) in (13) and choosing t=0, we obtain (37)|K1/2um′′(0)|2+(M(|A1/2u0m|2)Au0m+g(u1m),um′′(0))=0. Thus we have (38)|K1/2um′′(0)|2(|g(u1m)|+|M(|A1/2u0m|2)Au0m|)|um′′(0)|(|g(u1m)|+|M(|A1/2u0|2)Au0|)×|K1/2um′′(0)|. Thanks to the assumption (6), we deduce from (15) that (39)(g(u1m))is  bounded  in  L2(Ω). Therefore we conclude that the right-hand side is bounded; that is, (40)K1/2um′′(0)bounded  in  H.

A Priori Estimate III. For t<T, we apply (13) at points t and t+ζ such that 0<ζ<T-t. By taking the difference v=um(t+ζ)-um(t) in (13) and the assumption (G), we obtain (41)0(Kum′′(t+ζ)-Kum′′(t),um(t+ζ)-um(t))+(M(|A1/2um(t+ζ)|2)Aum(t+ζ)-M(|A1/2um(t)|2)Aum(t),um(t+ζ)-um(t)). Thus we have (42)012ddt[|K1/2(um(t+ζ)-um(t))|2]+M(|A1/2um(t+ζ)|2)×(Aum(t+ζ)-Aum(t),um(t+ζ)-um(t))+[M(|A1/2um(t+ζ)|2)-M(|A1/2um(t)|2)]×(Aum(t),um(t+ζ)-um(t)). Set (43)𝚽ζm(t)=|K1/2(um(t+ζ)-um(t))|2. By using (42), Young's inequality, the assumption (M), and the fact that K is positive self-adjoint operator, we see that 𝚽ζm(t)c𝚽ζm(t). Therefore we deduce (44)𝚽ζm(t)𝚽ζm(0)exp(cT)t[0,T]. Dividing the two sides of (44) by ζ2, letting ζ0, and using (43), we deduce (45)c|um′′|2|K1/2um′′(0)|2. From (40), it follows that |um′′|2C.

Since umC2[0,T], the previous inequality is verified for all t[0,T]. Therefore we conclude that (46)um′′  bounded  in  L(0,T;H). Moreover, using (19) and (46), it follows that (47)um  bounded  in  L2(0,T;H),um′′  bounded  in  L2(0,T;H). Applying a compactness theorem given in , we obtain (48)um  precompact  in  L2(0,T;H).

Passage to the Limit. Applying the Dunford-Pettis theorem, we conclude from (19), (21), (36), and (46)-(48), replacing the sequence um with a subsequence if needed, that (49)umu  weak-star  in  L(0,T;V),(50)um  u  weak-star  in  L(0,T;H),(51)um′′u′′  weak-star  in  L(0,T;H),(52)um  u  a.e  in  Ω×[0,T],(53)g(um  )ψweak-starin  L(q+1)/q(0,T;H),(54)M(|A1/2um|2)Aumχ  weak-star  in  L(0,T;H) for suitable functions uL(0,T;V), χL(0,T;H), and ψL(q+1)/q(Ω×[0,T]).

Now we are going to show that u is a solution of the problem (11)-(12). Indeed, from (49) to (51), we have (55)Ωum(0)ejdxΩu(0)ejdx,  Ωum(0)ejdxΩu(0)ejdx for each fixed j1. So we conclude that, for any j1, (56)Ω(um(0)-u0)ejdx=Ω(u(0)-u1)ejdx  =  0hhhhhhhhhhhhhhhhhhhhhhhhhhhhasm, which shows that (12) holds.

We will prove that, in fact, χ=M(|A1/2u|2)Au; that is, (57)M(|A1/2um|2)AumM(|A1/2u|2)Auhhhhhhhhhhhhiweak-star  in  L(0,;H).

For vL2(0,T;H), we have (58)0T(χ-M(|A1/2u|2)Au,v)dt=0T(χ-M(|A1/2um|2)Aum,v)dt+0TM(|A1/2u|2)(Aum-Au,v)dt+0T(M(|A1/2um|2)-M(|A1/2u|2))×(Aum,v)dt.

We deduce from (49) and (54) that the first and second terms in (58) tend to zero as m. For the last term, using the fact that M is C1 and (21), we can derive (with some positive constants c1,c2) (59)0T(M(|A1/2um|2)-M(|A1/2u|2))(Aum,v)dtc10T|A(um+u),um-u|dtc2(0T|um-u|2dt)1/2. Since um is bounded in L(0,T;V) and the injection of V in H is compact, we have (60)umu  strongly  in  L2(0,T;H). From (58) to (60), we deduce (57). It follows from (49), (51), and (57) that, for each fixed vLq+1(0,T;V), (61)0T(Kum′′+M(|A1/2um|2)Aum,v)dt0T(Ku′′+M(|A1/2u|2)Au,v)dt as m+.

For the nonlinear term, g(u), it remains to show that, for any fixed vLq+1(0,T;V), (62)0TΩvg(um)dxdt0TΩvg(u)dxdt as m.

At this moment we use the following lemma due to Jung and Choi (see [26, page 12]).

Lemma 2.

Suppose that Ω×[0,T] is a bounded open domain of n×; gm and g are in Lq(Ω×[0,T]), 1<q<, such that gmg a.e., in Ω×[0,T]. Then gmg weakly in Lq(Ω×[0,T]).

From (53), g(um)g(u) a.e. in Ω×[0,T]. By (36), we can use the above lemma and so we have ψ=g(u); that is, (63)g(um)g(u)  weak  in  L(q+1)/q(Ω×(0,T)), which implies (62). Therefore we obtain (64)0T(Ku′′+M(|A1/2u|2)Au+g(u),v)dt=0,hhhhhhhhhhhhhhhhhhhhhivLq+1(0,T;V).

The uniqueness is obtained by a standard method, so we omit the proof here.

3. Energy Estimates

In this section we study the energy estimate under suitable growth conditions on g.

Let us assume that there exist a number p1 and positive constants c1, i=1,2, such that (65)c1min{|K1/2x|,|K1/2x|p}|g(x)|c2max{|K1/2x|,|K1/2x|1/p} for all x.

Theorem 3.

Assume that (65) holds. Then one obtains the following energy decay: (66)E(t){c0E(0)e-wtt0,ifp=1,c~0(1+t)-2/(p-1)t0,ifp>1, where c0, w, and c~0 are some positive constants.

Proof.

Let T>0 be arbitrary and fixed and let uL(0,T;V)W2,(0,T;H) be a solution of (11) and (12). Multiplying (11) by u and integrating by parts in Ω×(s,T)(0s<T), we obtain that (67)E(T)-E(s)=-sT(g(u(t)),u(t))dt. By (g(u(t)),u(t))0 and being the primitive of an integrable function, it follows that the energy E is nonincreasing, locally absolutely continuous and (68)E(t)=-(g(u(t)),u(t))a.e.  in[0,). Here and in what follows we will denote by c diverse positive constants. We are going to show that the energy of this solution satisfies (69)sTE(t)(p+1)/2cE(s)0sT<. Once (69) is satisfied, the integral inequalities given in Komornik  and Haraux  will establish (66).

Now, multiplying (11) by E(t)(p-1)/2u and integrating by parts, we have (70)0=sTE(t)(p-1)/2(Ku′′+M(|A1/2u|2)Au+g(u),u)dt=[E(t)(p-1)/2(Ku,u)]sT-p-12sTE(t)(p-3)/2E(t)(Ku,u)dt-sTE(t)(p-1)/2|K1/2u|2dt+sTE(t)(p-1)/2(M(|A1/2u|2)|A1/2u|2,u)dt+sTE(t)(p-1)/2(g(u),u)dt. Note that by the assumption (M) and (21), we can choose some positive number (71)α=maxs[0,|A1/2u|2]{M(s)}< so that 2E(t)|K1/2u|2+α|A1/2u|2. Thus we deduce that (72)2βαsTE(t)(p+1)/2dt-[E(t)(p-1)/2(Ku,u)]sT+p-12sTE(t)(p-3)/2E(t)(Ku,u)dt+sTE(t)(p-1)/2((1+α-1)|K1/2u|2-(g(u),u))dtI1+I2+I3. Using the continuity of the imbedding VH, the Cauchy-Schwarz and the Young inequalities, we obtain (73)|(Ku,u)|c|Ku|ucE(t). Hence, since E(t) is nonincreasing, we obtain (74)I1cE(p-1)/2(0)E(s),I2(p-1)2sTcE(t)(p-1)/2E(t)dtI2cE(p-1)/2(0)E(s).

In order to estimate the last term I3 of (72), we set (75)Ω1={xΩ:|K1/2u(t,x)|1},Ω2={xΩ:|K1/2u(t,x)|>1}. Then we have (76)Ω|K1/2u(t,x)|2dx=Ω1|K1/2u(t,x)|2dx+Ω2|K1/2u(t,x)|2dx. The Hölder inequality yields (77)Ω|K1/2u(t,x)|2dxc(Ω1|K1/2u(t,x)|p+1dx)2/(p+1)+Ω2|K1/2u(t,x)|2dxJ1+J2. Using (65) and (68), we deduce that (78)J1c(Ω1ug(u)dx)2/(p+1)c|E(t)|2/(p+1),J2cΩ2|ug(u)|dxc(-E(t)). Combining these two inequalities with (77), we obtain (79)Ω|K1/2u(t,x)|2dxc(-E(t))2/(p+1)+c(-E(t)). Applying Young's inequality, it follows that, for any ϵ>0, (80)sTE(t)(p-1)/2|K1/2u|2dtϵcsTE(t)(p+1)/2dt+c(ϵ(1-p)/2+E(p-1)/2(0))E(s). It remains to estimate the second term of I3. Using (88) we have (81)|Ω1ug(u)dx|cuL(p+1)/p(Ω1)g(u)Lp+1(Ω1)cuL(p+1)/p(Ω1)(Ω1ug(u)dx)1/(p+1)cE(t)1/2    (-E(t))1/(p+1). Similarly, using (6), we obtain (82)|Ω2ug(u)dx|cuL2(Ω2)g(u)L2(Ω2)cuL2(Ω2)ug(u)L1(Ω2)1/2cE(t)1/2(-E(t))1/2. From (81) and (82), we deduce (83)|Ωug(u)dx|cE(t)1/2(-E(t))1/(p+1)+cE(t)1/2(-E(t))1/2. Using Young's inequality and (84)E(t)p/2(-E(t))1/2=E(t)(p+1)/4(E(t)(p-1)/4(-E(t))1/2), it follows from (82) that, for any ϵ>0, (85)-sTE(t)(p-1)/2(g(u),u)dt=-sTE(t)(p-1)/2Ωug(u)dxdtcsTE(t)p/2(-E(t))1/(p+1)dt+csTE(t)p/2(-E(t))1/2dt2ϵcsTE(t)(p+1)/2dt+c(ϵ-p+ϵ-1E(0)(p-1)/2)×sT(-E(t))dt2ϵcsTE(t)(p+1)/2dt+c(ϵ-p+ϵ-1E(0)(p-1)/2)E(s). Combining (80) with (85) and setting α~=1+α-1, we obtain (86)I3sTE(t)(p-1)/22|K1/2u|2dt+sTE(t)(p-1)/2Ωug(u)dxdt(α~+2)ϵcsTE(t)(p+1)/2dt+c(α~ϵ(1-p)/2+ϵ-p+(α~+ϵ-1)E(0)(p-1)/2)E(s). Therefore we conclude that (87)(2βα-(α~+2)ϵc)sTE(p+1)/2c(α~ϵ(1-p)/2+ϵ-p+(α~+ϵ-1)E(0)(p-1)/2)E(s). Now we choose ϵ as ϵ(0,2β/(3α+1)c); then (69) follows.

4. Numerical Result

In this section, we consider a Kirchhoff-type equation with heterogeneous string as an application: (88)(A(x)ρ)u′′(x,t)-(1+01|u(x,t)|2dx)Δu(x,t)+κ|u(x,t)|2u(x,t)=0,(89)in(x,t)(0,1)×(0,3),(90)u(0,t)=u(1,t)=0on  (0,3),(91)u0=u(x,0)=exp(-64(x-12)2)in  (0,1),(92)u1=ut(x,0)=0in  (0,1), where κ is a positive constant and A(x), ρ are given in Table 1.

Simulation parameters which are satisfied by theoretical conditions.

Symbols Definition Values Reference
A ( x ) Cross-sectional area 0.7853 ( 1 0 - 4 sin ( 2 10 π x ) + 1 )  cm2 
ρ Mass density of the unit length 7.850 g/cm2 

Then, the operators K=A(x)ρI(I:HH;identity  operator), A=-Δ, and the functions M(s)=s+1 and g(x)=κ|x|2x so that we can easily check that the hypotheses (M), (G), (H), and (S) in Preliminaries are satisfied. The smallness condition satisfies (u02+1)Δu020.2131/4. Therefore, by Theorem 1, we can deduce the following results.

Theorem 4.

For any T>0, there is a unique solution uL(0,T;H2(0,1))    W1,(0,T;H01(0,1))W2,(0,T;L2(0,1)) to the system (88)–(92).

The energy for the system (88)–(92) is given by (93)E(t)=12[+12(01|u(x,t)|2dx)201|A(x)ρu(x,t)|2dx+01|u(x,t)|2dxhhhhhhhhi+12(01|u(x,t)|2dx)2]. Next, in order to get the energy decay of (88)–(92), we need the value of the parameter p in (65). We can easily check that p=3 when g(x)=κ|x|2x.

Therefore, by Theorem 3, we get the energy decay rates for the energy E(t) as follows.

Theorem 5.

We obtain the following energy decay: (94)E(t)c1(1+t)-1t0, where c1 is a positive constant.

For the numerical simulation, we use the finite difference methods (FDM) which are the implicit multistep methods in time and second-order central difference methods for the space derivative in space in numerical algorithms (see [8, 9, 11]).

Figures 1(a)1(d) show displacements of solutions to the system (88)–(92) with κ=10 and κ=10-0.3, respectively.

Solution shapes and contour lines with respect to κ=10 and κ=10-0.3.

Temporal and spatial solution shapes in case of κ=10

Temporal and spatial solution shapes in case of κ=10-0.3

Temporal and spatial solution contour line in case of κ=10

Temporal and spatial solution contour line in case of κ=10-0.3

In case of κ=10 and κ=10-0.3, we deduce the algebraic decay rate for the energy as shown in Figure 2, respectively. The blue line and red dotted circled line (or blue circled line) show c1(t+1)1 and E(t) per the two values, respectively, where the parameter value c1=30.2 in (94). This result shows that the energy decay rates for solutions are algebraic in case that the system (88)–(92) with the nonlinear damping term κ|ut|2ut.

Algebraic decay rates of the energy in case of κ=10 and κ=10-0.3.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The first author's research was supported by Basic Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. NRF-2013R1A1A2010704). The corresponding author’s research was supported by Basic Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (Grant no. NRF-2012R1A1B3000599).

Kirchhoff G. Vorlesungen Über Mechanik 1883 Leipzig, Germany Teubner Nguyen Q. C. Hong K. -S. Transverse vibration control of axially moving membranes by regulation of axial velocity IEEE Transactions on Control Systems Technology 2012 20 4 1124 1131 Menzala G. P. Pazoto A. F. Zuazua E. Stabilization of Berger-Timoshenko's equation as limit of the uniform stabilization of the von Kármán system of beams and plates Mathematical Modelling and Numerical Analysis 2002 36 4 657 691 2-s2.0-0036658975 10.1051/m2an:2002029 Eisley J. G. Nonlinear vibration of beams and rectangular plates Zeitschrift für Angewandte Mathematik und Physik 1964 15 2 167 175 2-s2.0-20844445355 10.1007/BF01602658 Woinowsky-Krieger S. The effect of an axial force on the vibration of hinged bars Journal of Applied Mechanics 1950 17 35 36 Chen L. Lim C. W. Ding H. Energetics and conserved quantity of an axially moving string undergoing three-dimensional nonlinear vibration Acta Mechanica Sinica 2008 24 2 215 221 2-s2.0-41249094710 10.1007/s10409-007-0110-5 Hong K.-S. Asymptotic behavior analysis of a coupled time-varying system: application to adaptive systems IEEE Transactions on Automatic Control 1997 42 12 1693 1697 2-s2.0-0031387263 10.1109/9.650018 Kim D. Kim S. Jung I. H. Stabilization for the Kirchhoff type equation from an axially moving heterogeneous string modeling with boundary feedback control Nonlinear Analysis. Theory, Methods and Applications 2012 75 8 3598 3617 2-s2.0-84857918322 10.1016/j.na.2012.01.018 Kim D. Jung I. H. Asymptotic behavior of a nonlinear Kirchhoff type equation with spring boundary conditions Computers and Mathematics with Applications 2011 62 8 3004 3014 2-s2.0-80053368674 10.1016/j.camwa.2011.08.011 Kim D. Kang Y. H. Lee M. J. Jung I. H. Energy decay rate for a quasi-linear wave equation with localized strong dissipation Computers and Mathematics with Applications 2011 62 1 164 172 2-s2.0-79959498645 10.1016/j.camwa.2011.04.064 Kim D. Kang Y. H. Lee J. B. Ko G. R. Jung I. H. Stabilization for a nonlinear Kirchhoff equation by boundary feedback control Journal of Engineering Mathematics 2012 77 197 209 Liu I.-S. Rincon M. A. Effect of moving boundaries on the vibrating elastic string Applied Numerical Mathematics 2003 47 2 159 172 2-s2.0-0141959088 10.1016/S0168-9274(03)00063-1 Pellicano F. Vestroni F. Complex dynamics of high-speed axially moving systems Journal of Sound and Vibration 2002 258 1 31 44 2-s2.0-0344891120 10.1006/jsvi.2002.5070 Shahruz S. M. Parasurama S. A. Suppression of vibration in the axially moving Kirchhoff string by boundary control Journal of Sound and Vibration 1998 214 3 567 575 2-s2.0-0042598844 Bernstein S. Sur une classe d’équations fonctionelles aux dérivées partielles Izvestiya Akademii Nauk SSSR. Seriya Matematicheskaya 1940 4 17 26 Dickey R. W. The initial value problem for a nonlinear semi-infinite string Proceedings of the Royal Society of Edinburgh 1978 82 19 26 Ebihara Y. Medeiros L. A. Milla Miranda M. Local solutions for a nonlinear degenerate Hyperbolic equation Nonlinear Analysis. Theory, Methods and Applications 1986 10 1 27 40 2-s2.0-38249042627 Ha T. G. Kim D. Jung I. H. Grobal existence and uniform decay rates for the semi-linear wave equation with damping and source terms Computers and Mathematics with Applications 2014 67 692 707 Levine H. A. Instability and nonexistence of global solutions to nonlinear equations of the form Putt=Au+(u) Transactions of the American Mathematical Society 1974 192 1 21 Lions J. E. On some questions in boundary value problems of mathematical physics Proceedings of the International Symposium on Continuum Mechanics and Partial Differential Equations 1978 Amsterdam, The Netherlands North Holland Yamada Y. Some nonlinear degenerate wave equations Nonlinear Analysis 1987 11 10 1155 1168 2-s2.0-38249036775 Pohožaev S. I. On a class of quasilinear hyperbolic equations Matematicheskii Sbornik 1975 25 145 158 Mizumachi T. Time decay of solutions to degenerate Kirchhoff type equation Nonlinear Analysis. Theory, Methods and Applications 1998 33 3 235 252 2-s2.0-0032140410 Nishihara K. Yamada Y. On global solutions of some degenerate quasilinear hyperbolic equations with dissipative terms Funkcialaj Ekvacioj 1990 33 151 159 Park J. Y. Jung I. H. Kang Y. H. Generalized quasilinear hyperbolic equations and Yosida approximations Journal of the Australian Mathematical Society 2003 74 1 69 86 2-s2.0-0037302872 Jung I. H. Choi J. Energy decay estimates for a Kirchhoff model with viscosity Bulletin of the Korean Mathematical Society 2006 43 2 245 252 2-s2.0-33745100710 Jung I. H. Lee Y.-H. Exponential decay for the solutions of wave equations of kirchhoff type with strong damping Acta Mathematica Hungarica 2001 92 1-2 163 170 2-s2.0-26444493647 Ghisi M. Gobbino M. Global existence and asymptotic behaviour for a mildly degenerate dissipative hyperbolic equation of Kirchhoff type Asymptotic Analysis 2004 40 1 25 36 2-s2.0-23844506200 Aassila M. Kaya D. On local solutions of a mildly degenerate hyperbolic equation Journal of Mathematical Analysis and Applications 1999 238 2 418 428 2-s2.0-0011241645 Kouémou-Patcheu S. Global existence and exponential decay estimates for a damped quasilinear equation Communications in Partial Differential Equations 1997 22 11-12 2007 2024 2-s2.0-0009329664 Komornik V. Exact Controllability and Stabilizaion—The Multiplier Method 1994 New York, NY, USA John Wiley Haraux A. Oscillations forcées pour certains systèmes dissipatifs nonlinéaires Publications du Laboratoire d’Analyse Numérique, No. 78010 Université Pierre et Marie Curie, Paris, France, 1978 Lions J. L. Quelques Méthodesde Résolution des Prolèmes aux Limites Non Linéarires 1968 Paris, France Dunod Gauthiers-Villars Nguyen Q. C. Hong K. -S. Simultaneous control of longitudinal and transverse vibrations of an axially moving string with velocity tracking Journal of Sound and Vibration 2012 331 13 3006 3019