Research

Asymptotic behaviour of solution for multidimensional viscoelasticity equation with nonlinear source term

Runzhang Xu1*, Jie Liu1, Yi Niu2 and Shaohua Chen3

Author Affiliations

1 College of Science, Harbin Engineering University, Harbin, 150001, People’s Republic of China

2 College of Automation, Harbin Engineering University, Harbin, 150001, People’s Republic of China

3 Department of Mathematics, Cape Breton University, Sydney, B1P 6L2, Canada

For all author emails, please log on.

Boundary Value Problems 2013, 2013:42  doi:10.1186/1687-2770-2013-42

 Received: 24 April 2012 Accepted: 11 February 2013 Published: 1 March 2013

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

In this paper we study the initial-boundary value problem of the multidimensional viscoelasticity equation with nonlinear source term u t t Δ u t i = 1 N x i σ i ( u x i ) = f ( u ) . By using the potential well method, we first prove the global existence. Then we prove that when time t + , the solution decays to zero exponentially under some assumptions on nonlinear functions and the initial data.

1 Introduction

This paper considers the initial-boundary value problem (IBVP) of the multidimensional viscoelasticity equation with nonlinear source term

(1.1)

(1.2)

(1.3)

where u ( x , t ) is the unknown function with respect to the spacial variable x Ω and the time variable t , Ω R N is a bounded domain.

The viscoelasticity equation

u t t u x x t = σ ( u x ) x (1.4)

was suggested and studied by Greenberg et al.[1,2] from viscoelasticity mechanics in 1968. Under the condition σ ( s ) > 0 and higher smooth conditions on σ ( s ) and the initial data, they obtained the global existence of classical solutions for the initial-boundary value problem of Eq. (1.4).

After that many authors [3-11] studied the global well-posedness of IBVP for Eq. (1.4). In [3-10] the global existence, uniqueness and stability of solution were studied thoroughly. And in [11] the blow up of solution was discussed. Furthermore, in [12-14] the global existence of solution for IBVP of some multidimensional viscoelasticity equation was considered. And in [11] the blow up of solution for IBVP of the multidimensional generalisation of Eq. (1.4) was proved. Recently, in [15] and [16], the IBVP of the multidimensional viscoelasticity equation with nonlinear damping and source terms

(1.5)

(1.6)

(1.7)

was studied, and by using the potential well method, the global existence of weak solution was proved under some assumptions on nonlinear functions σ i ( s ) , f ( s ) , g ( s ) and the initial data. But we do not know how the global solution behaves as the time goes to infinity, namely the asymptotic behaviour of problem (1.1)-(1.3) is still open up to now. In the present paper, we try to study this problem by the multiplier method [17-22].

The main purpose of present paper is to consider the asymptotic behaviour of solution for problem (1.1)-(1.3). Since in the proof of the asymptotic behaviour of solution the global existence theory is required, it is necessary to give the proof of global existence of solution for problem (1.1)-(1.3).

In this paper, suppose that σ ( s ) = ( σ 1 ( s ) , , σ N ( s ) ) and f ( s ) satisfy the following assumptions:

where the constants in (H1) and (H2) are all positive and satisfy

In this paper, we first give some definitions and lemmas (Section 2). Then we prove the global existence of solution (Section 3). Finally, we prove the asymptotic behaviour of solution (Section 4).

In this paper, we denote L p ( Ω ) by p , = L 2 ( Ω ) and ( u , v ) = Ω u v d x .

2 Preliminaries

In this section, we will give some definitions and prove some lemmas for problem (1.1)-(1.3).

For problem (1.1)-(1.3), we define

Remark 2.1 Note that the definitions of J ( u ) and I ( u ) in the present paper are different from those in [11] and [15]. The definitions given in this paper will be shown more natural and rational because they are a part of the total energy E ( t ) .

Lemma 2.2Let (H1) and (H2) hold. Set

σ i ¯ ( s ) = σ i ( s ) a s , G i ¯ ( s ) = 0 s σ i ¯ ( s ) d τ .

Then the following hold:

(i) σ i ¯ ( s ) is increasing and s σ i ¯ ( s ) 0 s R ;

(ii) 0 G i ¯ ( s ) s σ i ¯ ( s ) s R .

Proof This lemma follows from σ i ¯ ( 0 ) = 0 and σ ¯ i ( s ) 0 . □

Lemma 2.3Let (H1) and (H2) hold, u W 0 1 , m + 1 ( Ω ) . Then the following hold:

(i) If 0 < u m + 1 < r ( δ ) , then I δ ( u ) > 0 ;

(ii) If I δ ( u ) < 0 , then u m + 1 > r ( δ ) ;

(iii) If I δ ( u ) = 0 , then u m + 1 r ( δ ) ,

where

r ( δ ) = ( A δ b C q + 1 ) 1 q m , C = sup u W 0 1 , m + 1 ( Ω ) / 0 u q + 1 u m + 1 .

Proof

(i) If 0 < u m + 1 < r ( δ ) , then we have

Ω u f ( u ) d x Ω | u f ( u ) | d x b Ω | u | q + 1 d x = b u q + 1 q + 1 b C q + 1 u m + 1 q + 1 = b C q + 1 A u m + 1 q m A u m + 1 m + 1 < δ i = 1 N Ω u x i σ i ( u x i ) d x ,

which gives I δ ( u ) > 0 .

(ii) If I δ ( u ) < 0 , then we have

δ A u m + 1 m + 1 δ i = 1 N Ω u x i σ i ( u x i ) d x < Ω u f ( u ) d x b C q + 1 u q m m + 1 u m + 1 m + 1 ,

which gives

u m + 1 > r ( δ ) .

(iii) If I δ ( u ) = 0 and u 0 , then by

δ A u m + 1 m + 1 δ i = 1 N Ω u x i σ i ( u x i ) d x < Ω u f ( u ) d x b C q + 1 u m + 1 q m u m + 1 m + 1 ,

we get

u m + 1 r ( δ ) .

□

Lemma 2.4Let (H1) and (H2) hold. Then the following holds:

d d 0 = ( p l ) A ( p + 1 ) ( l + 1 ) ( A b C q + 1 ) m + 1 q m . (2.1)

Proof For any u N , by Lemma 2.3, we have u m + 1 r ( 1 ) and

J ( u ) = i = 1 N Ω G i ( u x i ) d x Ω F ( u ) d x 1 l + 1 i = 1 N Ω u x i σ i ( u x i ) d x 1 p + 1 Ω u f ( u ) d x = ( 1 l + 1 1 p + 1 ) i = 1 N Ω u x i σ i ( u x i ) d x + 1 p + 1 I ( u ) = p l ( p + 1 ) ( l + 1 ) i = 1 N Ω u x i σ i ( u x i ) d x p l ( p + 1 ) ( l + 1 ) A u m + 1 m + 1 p l ( p + 1 ) ( l + 1 ) A r m + 1 ( 1 ) = ( p l ) A ( p + 1 ) ( l + 1 ) ( A b C q + 1 ) m + 1 q m ,

which gives (2.1). □

Now, for problem (1.1)-(1.3), we define

W = { u W 0 1 , m + 1 ( Ω ) I ( u ) > 0 } { 0 } .

3 Global existence of solution

In this section, we prove the global existence of weak solution for problem (1.1)-(1.3).

Definition 3.1 We call u = u ( x , t ) a weak solution of problem (1.1)-(1.3) on Ω × [ 0 , T ) if u L ( 0 , T ; W 0 1 , m + 1 ( Ω ) ) , u t L ( 0 , T ; L 2 ( Ω ) ) L 2 ( 0 , T ; H 0 1 ( Ω ) ) satisfying

(i)

(ii)

u ( x , 0 ) = u 0 ( x ) in  W 0 1 , m + 1 ( Ω ) ; u t ( x , 0 ) = u 1 ( x ) in  L 2 ( Ω ) .

Theorem 3.2Let (H1) and (H2) hold, u 0 ( x ) W 0 1 , m + 1 ( Ω ) , u 1 ( x ) L 2 ( Ω ) . Assume that E ( 0 ) < d , u 0 ( x ) W . Then problem (1.1)-(1.3) admits a global weak solution u L ( 0 , ; W 0 1 , m + 1 ( Ω ) ) and u t L ( 0 , ; L 2 ( Ω ) ) L 2 ( 0 , ; H 0 1 ( Ω ) ) .

Proof Let { w j ( x ) } j = 1 be a system of base functions in W 0 1 , m + 1 ( Ω ) . Construct the approximate solutions of problem (1.1)-(1.3)

u n ( x , t ) = j = 1 n g j n ( t ) w j ( x ) , n = 1 , 2 , ,

satisfying

(3.1)

(3.2)

(3.3)

Multiplying (3.1) by g s n ( t ) and summing for s, we get

d d t E n ( t ) + u n t 2 = 0 , 0 t < (3.4)

and

E n ( t ) + 0 t u n τ 2 d τ = E n ( 0 ) , 0 t < , (3.5)

where

E n ( t ) = 1 2 u n t 2 + J ( u n ) .

From (3.2) and (3.3), we have E n ( 0 ) E ( 0 ) as n . Hence, for sufficiently large n, we have E n ( 0 ) < d and

1 2 u n t 2 + J ( u n ) + 0 t u n τ 2 d τ < d , 0 t < . (3.6)

On the other hand, since W is an open set in W 0 1 , m + 1 ( Ω ) , Eq. (3.2) implies that for sufficiently large n, we have u n ( 0 ) W . Next, we prove that u n ( t ) W for 0 < t < and sufficiently large n. If it is false, then there exists a t 0 > 0 such that u n ( t 0 ) W , i.e. I ( u n ( t 0 ) ) = 0 and u n ( t 0 ) 0 , i.e. u n ( t 0 ) N . So, by the definition of d, we get J ( u n ( t 0 ) ) d , which contradicts (3.6).

From (3.6) we have

1 2 u n t 2 + i = 1 N Ω G i ( u n x i ) d x Ω F ( u n ) d x + 0 t u n τ 2 d τ < d , 0 t < ,

which gives

1 2 u n t 2 + 1 l + 1 i = 1 N Ω u n x i σ i ( u n x i ) d x 1 p + 1 Ω u n f ( u n ) d x + 0 t u n τ 2 d τ < d

and

which together with u n ( t ) W gives

1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) i = 1 N Ω u n x i σ i ( u n x i ) d x + 0 t u n τ 2 d τ < d ,

and

1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) A u n m + 1 m + 1 + 0 t u n τ 2 d τ < d , 0 t < . (3.7)

From (3.7) we can get

(3.8)

(3.9)

(3.10)

(3.11)

(3.12)

Hence there exist u, χ = ( χ 1 , χ 2 , , χ N ) , η and a subsequence { u ν } of { u n } such that as ν , u ν u in u L ( 0 , ; W 0 1 , m + 1 ( Ω ) ) weak-star, and a.e. in Q = Ω × [ 0 , ) , u ν t u t in L ( 0 , ; L 2 ( Ω ) ) weak-star and in L 2 ( 0 , ; H 0 1 ( Ω ) ) weakly, σ i ( u ν x i ) χ i = σ i ( u x i ) in L ( 0 , ; L ( m + 1 ) ( Ω ) ) weak-star, ( m + 1 ) = m + 1 m , f ( u v ) η = f ( u ) in L ( 0 , ; L ( q + 1 ) ( Ω ) ) weak-star, ( q + 1 ) = q + 1 q .

Integrating (3.1) with respect to t, we have

(3.13)

Letting n = ν in (3.13), we get

and

On the other hand, from (3.2) and (3.3), we get u ( x , 0 ) = u 0 ( x ) in W 0 1 , m + 1 ( Ω ) , u t ( x , 0 ) = u 1 ( x ) in L 2 ( Ω ) . Therefore u is a global weak solution of problem (1.1)-(1.3). □

4 Asymptotic behaviour of solution

In this section, we prove the main conclusion of this paper - the asymptotic behaviour of solution for problem (1.1)-(1.3).

Lemma 4.1Let (H1) and (H2) hold, u 0 ( x ) W 0 1 , m + 1 ( Ω ) , u 1 ( x ) L 2 ( Ω ) . Then, for the approximate solutions u n ( x , t ) of problem (1.1)-(1.3) constructed in the proof of Theorem 3.2, the following hold:

(i)

I ( u n ) = u n t 2 d d t ( ( u n t , u n ) + 1 2 u n 2 ) ; (4.1)

(ii) Furthermore, if E ( 0 ) < d 0 and u 0 ( x ) W , then for sufficiently largen, there exists a δ 1 ( 0 , 1 ) such that

I ( u n ) ( 1 δ 1 ) i = 1 N ( σ i ( u n x i ) , u n x i ) . (4.2)

Proof (i) Multiplying (3.1) by g s n ( t ) and summing for s, we get (4.1).

(ii) From

E ( 0 ) < d 0 = p l ( p + 1 ) ( l + 1 ) A ( A b C q + 1 ) m + 1 q m

it follows that there exists a δ 1 ( 0 , 1 ) such that

E ( 0 ) < p l ( p + 1 ) ( l + 1 ) A ( A δ 1 b C q + 1 ) m + 1 q m d ( δ 1 ) . (4.3)

From (3.2), (3.3) and (4.3), it follows that E m ( 0 ) < d ( δ 1 ) for sufficiently large n. Hence from (3.5) we have

and

i = 1 N Ω G i ( u n x i ) d x Ω F ( u n ) d x < d ( δ 1 ) ,

which gives

1 l + 1 i = 1 N Ω u n x i σ i ( u n x i ) d x 1 p + 1 Ω u n f ( u n ) d x < d ( δ 1 )

and

p l ( p + 1 ) ( l + 1 ) i = 1 N Ω u n x i σ i ( u n x i ) d x + 1 p + 1 I ( u n ) < d ( δ 1 ) ,

which together with u n W for sufficiently large n gives

p l ( p + 1 ) ( l + 1 ) A u n m + 1 m + 1 < p l ( p + 1 ) ( l + 1 ) A ( A δ 1 b C q + 1 ) m + 1 q m

and

u n m + 1 < ( A δ 1 b C q + 1 ) 1 q m = r ( δ 1 ) .

Hence, by Lemma 2.3, we have I δ 1 ( u n ) > 0 or u n = 0 . So, we have

I ( u n ) = i = 1 N Ω u n x i σ i ( u n x i ) d x Ω u n f ( u n ) d x = ( 1 δ 1 ) i = 1 N Ω u n x i σ i ( u n x i ) d x + I δ 1 ( u n ) ( 1 δ 1 ) i = 1 N Ω u n x i σ i ( u n x i ) d x .

□

Theorem 4.2Let (H1) and (H2) hold, u 0 ( x ) W 0 1 , m + 1 ( Ω ) , u 1 ( x ) L 2 ( Ω ) . Assume that E ( 0 ) < d 0 , u 0 ( x ) W . Then, for the global weak solutionugiven in Theorem 3.2, there exist positive constantsCandλsuch that

u t 2 + u m + 1 m + 1 C e λ t , 0 t < . (4.4)

Proof Let { u n } be the approximate solutions of problem (1.1)-(1.3) in the proof of Theorem 3.2, then (3.4) holds. Multiplying (3.4) by e α t ( α > 0 ), we get

d d t ( e α t E n ( t ) ) + e α t u n t 2 = α e α t E n ( t )

and

e α t E n ( t ) + 0 t e α τ u n τ 2 d τ = E n ( 0 ) + α 0 t e α τ E n ( τ ) d τ , 0 t < . (4.5)

From (H2), Lemma 2.2 and Lemma 4.1, we get

where

C ( r , δ 1 ) = 1 1 δ 1 r r + 1 + 1 r + 1 .

Hence we have

(4.6)

and

(4.7)

From

1 2 u n t 2 + i = 1 N Ω G i ( u n x i ) d x Ω F ( u n ) d x = E n ( t )

we get

1 2 u n t 2 + 1 l + 1 i = 1 N Ω u n x i σ i ( u n x i ) d x 1 p + 1 Ω u n f ( u n ) d x E n ( t )

and

1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) i = 1 N Ω u n x i σ i ( u n x i ) d x + 1 p + 1 I ( u n ) E n ( t ) ,

which together with u n W for sufficiently large n gives

1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) i = 1 N Ω u n x i σ i ( u n x i ) d x E n ( t ) , 1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) ( a u n 2 + i = 1 N Ω u n x i σ ¯ i ( u n x i ) d x ) E n ( t ) , (4.8)

and

1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) a u n 2 E n ( t ) , 0 t < . (4.9)

From (4.8) and the Poincaré inequality u 2 λ 1 u 2 , it follows that there exists a constant C 0 = C 0 ( p , l , a , λ 1 ) > 0 such that

( u n t 2 + u n 2 + u n 2 ) C 0 E n ( t ) , 0 t < . (4.10)

From (4.5)-(4.10) it follows that there exists a C 0 such that

(4.11)

Choose α such that

0 < α < min { 1 2 C 0 , λ 1 1 2 + C ( r , δ 1 ) } .

Then from (4.11) we get

e α t E n ( t ) 2 ( C 0 α + 1 ) E n ( 0 ) + 2 α 2 C 0 0 t e α τ E n ( τ ) d τ 2 ( C 0 α + 1 ) d 0 + 2 α 2 C 0 0 t e α τ E n ( τ ) d τ 3 d 0 + 2 α C 0 0 t e α τ E n ( τ ) d τ .

From this and the Gronwall inequality, we get

e α t E n ( t ) 3 d 0 e 2 α 2 C 0 t

and

E n ( t ) 3 d 0 e λ t , 0 t < , λ = α ( 1 2 C 0 α ) > 0 . (4.12)

On the other hand, from (4.8) we get

1 2 u n t 2 + p l ( p + 1 ) ( l + 1 ) A u n m + 1 m + 1 E n ( t ) , 0 t < .

Hence, there exists a C 1 = C 1 ( p , l , A ) such that

u n t 2 + u n m + 1 m + 1 C 1 E n ( t ) , 0 t < . (4.13)

Let { u ν } be the subsequence of { u n } in the proof of Theorem 3.2. Then from (4.13) and(4.12), we obtain

which gives (4.4), where C = 3 d 0 C 1 , λ = α ( 1 2 C 0 α ) . □

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The work presented here was carried out in collaboration between all authors. RZ organised this paper and found the topic of this paper. He introduced this problem and suggested the methods and the outline of the proofs. JL finished the proof of the global existence and YN finished the long time behaviour part. SC discussed all the problems arising in the research and provided many good ideas for proving the problems. All authors have contributed to, seen and approved the manuscript.

Acknowledgements

We thank the referees for their valuable suggestions which helped us improve the paper so much. This work was supported by the National Natural Science Foundation of China (11101102), Ph.D. Programs Foundation of the Ministry of Education of China (20102304120022), the Support Plan for the Young College Academic Backbone of Heilongjiang Province (1252G020), the Natural Science Foundation of Heilongjiang Province (A201014), Science and Technology Research Project of the Department of Education of Heilongjiang Province (12521401), Foundational Science Foundation of Harbin Engineering University and Fundamental Research Funds for the Central Universities.

References

1. Greenberg, JM: On the existence, uniqueness and stability of solutions of the equation ρ 0 X x x = E ( X x ) X x x + λ X x x t . J. Math. Anal. Appl.. 25, 575–591 (1969). PubMed Abstract | Publisher Full Text

2. Greenberg, JM, MacCamy, RC, Mizel, JJ: On the existence, uniqueness and stability of solutions of the equation ρ ( u x ) u x x + λ u x x t = ρ 0 u t t . J. Math. Mech.. 17, 707–728 (1968)

3. Andrews, G: On the existence of solutions to the equation u t t u x x t = σ ( u x ) x . J. Differ. Equ.. 35, 200–231 (1980). Publisher Full Text

4. Andrews, G, Ball, JM: Asymptotic behavior and changes in phase in one-dimensional nonlinear viscoelasticity. J. Differ. Equ.. 44, 306–341 (1982). PubMed Abstract | Publisher Full Text

5. Dafermos, CM: The mixed initial-boundary value problem for the equations of nonlinear one-dimensional visco-elasticity. J. Differ. Equ.. 6, 71–86 (1969). Publisher Full Text

6. Davis, P: A quasi-linear hyperbolic and related third-order equation. J. Math. Anal. Appl.. 51, 596–606 (1975). Publisher Full Text

7. Engler, H: Strong solutions for strongly damped quasilinear wave equations. Contemp. Math.. 64, 219–237 (1987)

8. Greenberg, JM, MacCamy, RC: On the exponential stability of solutions of E ( u x ) u x x + λ u x t x = ρ u t t . J. Math. Anal. Appl.. 31, 406–417 (1970). Publisher Full Text

9. Liu, Y, Liu, D: Initial boundary value problem, periodic boundary problem and initial value problem of equation u t t = u x x t + σ ( u x ) x . Chin. Ann. Math., Ser. A. 9, 459–470 (in Chinese) (1988)

10. Yamada, Y: Some remarks on the equation y t t σ ( y x ) y x x y x t x = f . Osaka J. Math.. 17, 303–323 (1980)

11. Yang, Z, Song, C: Blow-up of solutions for a class of quasilinear evolution equations. Nonlinear Anal., Theory Methods Appl.. 28, 2017–2032 (1997). Publisher Full Text

12. Ang, DD, Dinh, PN: Strong solutions of a quasilinear wave, equation with nonlinear damping. SIAM J. Math. Anal.. 102, 337–347 (1988)

13. Clements, J: Existence theorems for a quasilinear evolution equation. SIAM J. Appl. Math.. 226, 745–752 (1974)

14. Clements, J: On the existence and uniqueness of solutions of the equation u t t σ i ( u x i ) / x i D N u t = f . Can. Math. Bull.. 18, 181–187 (1975). Publisher Full Text

15. Yang, Z: Existence and asymptotic behaviour of solutions for a class of quasilinear evolution equations with nonlinear damping and source terms. Math. Methods Appl. Sci.. 25, 795–814 (2002). Publisher Full Text

16. Liu, Y, Zhao, J: Multidimensional viscoelasticity equations with nonlinear damping and source terms. Nonlinear Anal.. 56, 851–865 (2004). Publisher Full Text

17. Xu, R: Asymptotic behavior and blow-up of solutions for the viscous diffusion equation. Appl. Math. Lett.. 20, 255–259 (2007). Publisher Full Text

18. Xu, R, Zhao, X, Shen, J: Asymptotic behaviour of solution for fourth order wave equation with dispersive and dissipative terms. Appl. Math. Mech.. 29, 259–262 (2008). Publisher Full Text

19. Liu, Y, Xu, R, Yu, T: Global existence, nonexistence and asymptotic behavior of solutions for the Cauchy problem of semilinear heat equations. Nonlinear Anal.. 68, 3332–3348 (2008). Publisher Full Text

20. Xu, R, Liu, Y: Asymptotic behavior of solutions for initial-boundary value problems for strongly damped nonlinear wave equations. Nonlinear Anal.. 69, 2492–2495 (2008). Publisher Full Text

21. Xu, R: Global existence, blow up and asymptotic behaviour of solutions for nonlinear Klein-Gordon equation with dissipative term. Math. Methods Appl. Sci.. 33, 831–844 (2010)

22. Jiang, X, Xu, R: Global well-posedness for semilinear hyperbolic equations with dissipative term. J. Appl. Math. Comput.. 38, 467–487 (2012). Publisher Full Text