SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Asymptotic profile of solutions to the semilinear beam equation

Yunpeng Zhang1* and Yanling Li2

Author Affiliations

1 College of Electric Power, North China University of Water Resources and Electric Power, Zhengzhou, 450011, China

2 School of Mathematics and Information Sciences, North China University of Water Resources and Electric Power, Zhengzhou, 450011, China

For all author emails, please log on.

Boundary Value Problems 2014, 2014:84  doi:10.1186/1687-2770-2014-84


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2014/1/84


Received:26 January 2014
Accepted:1 April 2014
Published:17 April 2014

© 2014 Zhang and Li; licensee Springer.

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 credited.

Abstract

In this paper, we investigate the initial value problem for the semilinear beam equation. Under a small condition on the initial value, we prove the global existence and optimal decay estimate of solutions. Moreover, we show that as time tends to infinity, the solution is asymptotic to a diffusion wave, which is given explicitly in terms of the solution of parabolic equation.

MSC: 35L30, 35L75.

Keywords:
beam equation; decay estimate; asymptotic profile; diffusion wave

1 Introduction

We investigate the initial value problem for the following semilinear beam equation:

u t t α u x x + u t + u x x x x = f ( u x ) x (1.1)

with the initial value

t = 0 : u = u 0 ( x ) , u t = u 1 ( x ) . (1.2)

Here u = u ( x , t ) is the unknown function of x R and t > 0 , b > 0 , and α > 0 are constants. The nonlinear term f ( v ) is a given smooth function of v R . More precisely,

f ( v ) = f ( 0 ) + f ( 0 ) v + g ( v ) ,

where g ( v ) = O ( v 2 ) for v 0 .

This initial value problem was studied by [1,2] when f satisfies

{ | f ( v ) | C | v | p , | f ( v ) | C | v | p 1 , | f ( v 1 ) f ( v 2 ) | C ( | v 1 | p 1 + | v 2 | p 1 ) | v 1 v 2 | , | f ( v 1 ) f ( v 2 ) | C ( | v 1 | p 2 + | v 2 | p 2 ) | v 1 v 2 | (1.3)

with p 2 , so that f ( 0 ) = 0 . Here C is independent of v, v 1 , and v 2 . [1] proved that there exists a global solution u C ( [ 0 , ) ; H 2 L 1 ) to the problem (1.1), (1.2) under smallness condition on the initial data u 0 W 2 , 1 H 2 and u 1 L 2 L 1 . In particular, they showed the decay estimates:

u L q C ( 1 + t ) 1 2 ( 1 1 q ) , 1 q (1.4)

and

x 2 u L q C ( 1 + t ) 1 2 ( 1 1 q ) 1 , 2 q . (1.5)

In addition to the above assumptions, suppose that the initial data u 0 , u 1 L 1 1 , Takeda and Yoshikawa [2] established the following asymptotic profile of global solution:

u ( t ) M G α ( t ) ( M ˜ + N ) x G α ( t ) L q C ( 1 + t ) 1 2 ( 1 1 q ) 1 2 , 1 q , (1.6)

where

M = ( u 0 + u 1 ) ( x ) d x , G α ( x , t ) = F 1 [ e α ξ 2 t ] , M ˜ = x ( u 0 + u 1 ) ( x ) d x

and

N = 0 f ( x u ( x , t ) ) d x d t .

The main purpose of our present paper is two-fold: first, we try to recover all the results about global existence and decay estimate of solution in Takeda and Yoshikawa [1] under some assumptions on the initial data ( u 0 ( x ) , u 1 ( x ) ) and the nonlinear function f, which is much weaker than those needed in Takeda and Yoshikawa’s arguments. More precisely, the condition u 0 ( x ) W 2 , 1 H 2 , u 1 ( x ) L 2 L 1 and (1.3) have been relaxed to u 0 ( x ) H s + 2 L 1 , u 1 ( x ) H s L 1 ( s 0 ) and the nonlinear function f satisfies f ( 0 ) > α in this paper. Second, we show that the solution is asymptotic to a diffusion wave, given explicitly in terms of the solution of parabolic equation that is different from the one in [2]. For the details, we refer to Theorem 4.1. Moreover, under some additional assumptions on the initial data, we also prove that the convergence rates of our new asymptotic profile are better than that obtained by [2]. For details, we refer to Theorem 4.2.

The study of the global existence and asymptotic behavior of solutions to hyperbolic-type equations has a long history. We refer to [3,4] for hyperbolic equations, [5-7] for the damped wave equation and [8-16] for various aspects of dissipation of the plate equation.

The paper is organized as follows. In Section 2, we study the decay property of the solution to the linear problem. Then, in Section 3, we prove the global existence and decay estimate of the solutions. Finally, we prove that the solution is asymptotic to a diffusion wave, which is given explicitly in terms of the solution of the parabolic equation in Section 4.

Notations We give some notations which are used in this paper. Let F [ u ] denote the Fourier transform of u defined by

u ˆ ( ξ ) = F [ u ] = R n e i ξ x u ( x ) d x ,

and we denote its inverse transform by F 1 .

For 1 p , L p = L p ( R n ) denotes the usual Lebesgue space with the norm L p . The usual Sobolev space of s is defined by W s , p = ( I x 2 ) s 2 L p with the norm f W s , p = ( I x 2 ) s 2 f L p .

Finally, in this paper, we denote every positive constant by the same symbol C or c which will not lead to any confusion. [ ] is the Gauss symbol.

2 Linear problem

2.1 Solution formula

The aim of this subsection is to derive the solution formula for the problem (1.1), (1.2). We first investigate the linearized equation of (1.1):

u t t β u x x + u t + u x x x x = 0 , (2.1)

with the initial data in (1.2), where β = α + f ( 0 ) > 0 . We apply the Fourier transform to (2.1). This yields

u ˆ t t + ξ 2 u ˆ t + ( β ξ 2 + ξ 4 ) u ˆ = 0 . (2.2)

The corresponding initial values are given as

t = 0 : u ˆ = u ˆ 0 ( ξ ) , u ˆ t = u ˆ 1 ( ξ ) . (2.3)

The characteristic equation of (2.2) is

λ 2 + λ + β ξ 2 + ξ 4 = 0 . (2.4)

Let λ = λ ± ( ξ ) be the corresponding eigenvalues of (2.4), we obtain

λ ± ( ξ ) = 1 ± 1 4 ξ 2 ( β + ξ 2 ) 2 . (2.5)

The solution to the problem (2.2), (2.3) in the Fourier space is then given explicitly in the form

u ˆ ( ξ , t ) = G ˆ ( ξ , t ) ( u ˆ 1 ( ξ ) + u ˆ 0 ( ξ ) ) + H ˆ ( ξ , t ) u ˆ 0 ( ξ ) , (2.6)

where

G ˆ ( ξ , t ) = 1 λ + ( ξ ) λ ( ξ ) ( e λ + ( ξ ) t e λ ( ξ ) t ) (2.7)

and

H ˆ ( ξ , t ) = 1 λ + ( ξ ) λ ( ξ ) ( ( λ + ( ξ ) + 1 ) e λ ( ξ ) t ( λ ( ξ ) + 1 ) e λ + ( ξ ) t ) . (2.8)

We define G ( x , t ) and H ( x , t ) by

G ( x , t ) = F 1 [ G ˆ ( ξ , t ) ] ( x ) (2.9)

and

H ( x , t ) = F 1 [ H ˆ ( ξ , t ) ] ( x ) , (2.10)

respectively, where F 1 denotes the inverse Fourier transform. Then, applying F 1 to (2.6), we obtain

u ( t ) = G ( t ) ( u 1 + u 0 ) + H ( t ) u 0 . (2.11)

By the Duhamel principle, we obtain the solution formula to (1.1), (1.2),

u ( t ) = G ( t ) ( u 1 + u 0 ) + H ( t ) u 0 + 0 t G ( t τ ) g ( u x ) x ( τ ) d τ , (2.12)

where g ( v ) is a smooth function; it satisfies g ( v ) = O ( v 2 ) .

2.2 Decay property

The aim of this subsection is to establish decay estimates of the solution operators G ( t ) and H ( t ) appearing in (2.9) and (2.10), respectively.

Lemma 2.1The solution of the problem (2.2), (2.3) satisfies

| u ˆ t ( ξ , t ) | 2 + ( 1 + ξ 2 ) 2 | u ˆ ( ξ , t ) | 2 C e c ρ ( ξ ) t ( | u ˆ 1 ( ξ ) | 2 + ( 1 + ξ 2 ) 2 | u ˆ 0 ( ξ ) | 2 ) (2.13)

for ξ R and t 0 , where ρ ( ξ ) = ξ 2 1 + ξ 2 .

Proof Multiplying (2.2) by u ˆ ¯ t and taking the real part yields

1 2 d d t { | u ˆ t | 2 + β ξ 2 | u ˆ | 2 + ξ 4 | u ˆ | 2 } + | u ˆ t | 2 = 0 . (2.14)

Multiplying (2.2) by u ˆ ¯ and taking the real part, we obtain

1 2 d d t { | u ˆ | 2 + 2 Re ( u ˆ t u ˆ ¯ ) } + β ξ 2 | u ˆ | 2 + ξ 4 | u ˆ | 2 | u ˆ t | 2 = 0 . (2.15)

Multiplying both sides of (2.14) by 2 and summing up the resulting equation and (2.15) yields

d d t E + F = 0 , (2.16)

where

E = | u ˆ t | 2 + ( 1 2 + β ξ 2 + ξ 4 ) | u ˆ | 2 + Re ( u ˆ t u ˆ ¯ )

and

F = | u ˆ t | 2 + β ξ 2 | u ˆ | 2 + ξ 4 | u ˆ | 2 .

A simple computation implies that

C E 0 E C E 0 , (2.17)

where

E 0 = | u ˆ t | 2 + ( 1 + ξ 2 ) 2 | u ˆ | 2 .

Note that

F c ρ ( ξ ) E 0 .

It follows from (2.17) that

F c ρ ( ξ ) E . (2.18)

Using (2.16) and (2.18), we get

d d t E + c ρ ( ξ ) E 0 .

Thus

E ( ξ , t ) e c ρ ( ξ ) t E ( ξ , 0 ) ,

which together with (2.17) proves the desired estimates (2.13). Then we have completed the proof of the lemma. □

Lemma 2.2Let G ˆ ( ξ , t ) and H ˆ ( ξ , t ) be the fundamental solution of (2.1) in the Fourier space, which are given in (2.7) and (2.8), respectively. Then we have the estimates

| G ˆ t ( ξ , t ) | 2 + ( 1 + ξ 2 ) 2 | G ˆ ( ξ , t ) | 2 C e c ρ ( ξ ) t (2.19)

and

| H ˆ t ( ξ , t ) | 2 + ( 1 + ξ 2 ) 2 | H ˆ ( ξ , t ) | 2 C ( 1 + ξ 2 ) 2 e c ρ ( ξ ) t (2.20)

for ξ R and t 0 , where ρ ( ξ ) = ξ 2 1 + ξ 2 .

Proof Firstly, we investigate the problem (2.1), (1.2) with u ˆ 0 ( ξ ) = 0 , from (2.6), we obtain

u ˆ ( ξ , t ) = G ˆ ( ξ , t ) u ˆ 1 ( ξ ) , u ˆ t ( ξ , t ) = G ˆ t ( ξ , t ) u ˆ 1 ( ξ ) .

Substituting the equalities into (3.1) with u ˆ 0 ( ξ ) = 0 , we get (2.19).

In what follows, we consider the problem (2.1), (1.2) with u ˆ 1 ( ξ ) = 0 ; it follows from (2.6) that

u ˆ ( ξ , t ) = ( G ˆ + H ˆ ) ( ξ , t ) u ˆ 0 ( ξ ) , u ˆ t ( ξ , t ) = ( G ˆ t + H ˆ t ) ( ξ , t ) u ˆ 0 ( ξ ) .

Substituting the equalities into (2.13) with u ˆ 1 ( ξ ) = 0 , we get the desired estimate (2.20). The lemma is proved. □

Let

G 0 ( x , t ) = F 1 [ e β ξ 2 t ] ( x ) = 1 4 π β t e e x 2 4 β t

be the fundamental solution to u t β u x x = 0 .

Lemma 2.3Let G ˆ ( ξ , t ) and H ˆ ( ξ , t ) be the fundamental solution of (2.1) in the Fourier space, which are given in (2.6) and (2.7), respectively. Then there is a small positive number R 0 such that if | ξ | R 0 and t 0 , we have the following estimates:

| ( G ˆ G ˆ 0 ) ( ξ , t ) | + | H ˆ ( ξ , t ) | C ξ 2 e c ξ 2 t + C e c t , (2.21)

| G ˆ t ( ξ , t ) | C ξ 2 e c ξ 2 t + C e c t (2.22)

and

| t ( G ˆ G ˆ 0 ) ( ξ , t ) | + | H ˆ t ( ξ , t ) | C ξ 4 e c ξ 2 t + C e c t . (2.23)

Proof For sufficiently small ξ, using the Taylor formula, we get

λ + ( ξ ) = β ξ 2 + O ( ξ 4 ) , λ ( ξ ) = 1 + O ( ξ 2 ) , 1 λ + ( ξ ) λ ( ξ ) = 1 + 2 β ξ 2 + O ( ξ 4 ) . (2.24)

We rewrite G ˆ ( ξ , t ) in (2.6) as

G ˆ ( ξ , t ) = e β ξ 2 t + ( 1 λ + ( ξ ) λ ( ξ ) 1 ) e λ + ( ξ ) t + ( e ( λ + ( ξ ) + β ξ 2 ) t 1 ) e β ξ 2 t e λ ( ξ ) t λ + ( ξ ) λ ( ξ ) . (2.25)

For sufficiently small ξ, from (2.24) and (2.25), we immediately obtain

| ( G ˆ G ˆ 0 ) ( ξ , t ) | C ξ 2 e c ξ 2 t + C e c t .

For sufficiently small ξ, from (2.7) and (2.24), we immediately get

| H ˆ ( ξ , t ) | C ξ 2 e c ξ 2 t + C e c t .

Thus we get (2.21). The other estimates are proved similarly and we omit the details. The proof of Lemma 2.3 is completed. □

Lemma 2.4Let G ( x , t ) and H ( x , t ) be the fundamental solutions of (2.1), which are given in (2.9) and (2.10), respectively. Let 1 p 2 , and letk, j, andlbe nonnegative integers. Then we have

x k G ( t ) ϕ L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k j 2 x j ϕ L p + C e c t x k + l 2 ϕ L 2 , (2.26)

x k H ( t ) ϕ L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 2 j 2 x j ϕ L p + C e c t x k + l ϕ L 2 (2.27)

for 0 j k , where k + l 2 0 in (2.26). Similarly, we have

x k G t ( t ) ϕ L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 2 j 2 x j ϕ L p + C e c t x k + l ϕ L 2 , (2.28)

x k H t ( t ) ϕ L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 4 j 2 x j ϕ L p + C e c t x k + l + 2 ϕ L 2 (2.29)

for 0 j k + 2 in (2.28) and 0 j k + 4 in (2.29).

Proof We only prove (2.26). By the Plancherel theorem and (2.19), the Hausdorff-Young inequality, we obtain

x k G ( t ) ϕ L 2 2 = | ξ | 1 | ξ | 2 k | G ˆ ( ξ , t ) | 2 | ϕ ˆ ( ξ ) | 2 d ξ + | ξ | 1 | ξ | 2 k | G ˆ ( ξ , t ) | 2 | ϕ ˆ ( ξ ) | 2 d ξ C | ξ | 1 | ξ | 2 k e c ξ 2 t | ϕ ˆ ( ξ ) | 2 d ξ + C e c t | ξ | 1 | ξ | 2 k ( 1 + | ξ | 2 ) 2 | ϕ ˆ ( ξ ) | 2 d ξ = : A 1 + A 2 .

For the term A 1 , letting 1 p + 1 p = 1 , we have

A 1 C | ξ | 1 | ξ | 2 k 2 e c ξ 2 t | ϕ ˆ ( ξ ) | 2 d ξ C | ξ | j ϕ ˆ L p 2 ( | ξ | 1 | ξ | 2 ( k j ) p e c q ξ 2 t d ξ ) 1 p C ( 1 + t ) ( 1 p 1 2 ) ( k j ) x j ϕ L ˙ p 2 ,

where we used the Hölder inequality with 2 p + 1 q = 1 and the Hausdorff-Young inequality v ˆ L p C v L p for v = x j ϕ . On the other hand, we can estimate the term A 2 simply as

A 2 C e c t | ξ | 1 | ξ | 2 k 4 | ϕ ˆ ( ξ ) | 2 d ξ C e c t | ξ | 1 | ξ | 2 ( k + l 2 ) | ϕ ˆ ( ξ ) | 2 d ξ C e c t x k + l 2 ϕ L 2 2 ,

where k + l 2 0 .

Combining the above three inequalities yields (2.26). This completes the proof of Lemma 2.4. □

From Lemma 2.4, we immediately have the following corollary.

Corollary 2.1Let G ( x , t ) and H ( x , t ) be the fundamental solution of (2.1), which are given in (2.8) and (2.9), respectively. Let 1 p 2 , and letk, j, andlbe nonnegative integers. Then we have

x k G ( t ) g x L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 1 j 2 x j g L p + C e c t x k + l 1 g L 2 (2.30)

for 0 j k + 1 and k + l 1 0 . Also we have

x k G t ( t ) g x L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 3 j 2 x j g L p + C e c t x k + l + 1 g L 2 (2.31)

for 0 j k + 3 .

Lemma 2.5Let G ( x , t ) be the fundamental solution of (2.1), given in (2.8) and let G 0 ( x , t ) be the fundamental solution of (2.1), given in (2.8). Let 1 p 2 , and letk, j, andlbe nonnegative integers. Then we have

x k ( G G 0 ) ( t ) ϕ L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 2 j 2 x j ϕ L p + C e c t x k + l ϕ L 2 (2.32)

for 0 j k + 2 . Similarly,

x k ( G G 0 ) t ( t ) ϕ L 2 C ( 1 + t ) 1 2 ( 1 p 1 2 ) k + 4 j 2 x j ϕ L p + C e c t x k + l + 2 ϕ L 2 (2.33)

for 0 j k + 4 .

Proof The proof of Lemma 2.5 is similar to the proof of Lemma 2.4. By employing (2.21) and (2.23), we can prove Lemma 2.5. We omit the details. □

3 Global existence and asymptotic behavior of solutions to (1.1), (1.2)

The purpose of this section is to prove global existence and optimal decay estimate of solutions to the initial value problem (1.1), (1.2). We need the following lemma, which comes from [17] (see also [18]).

Lemma 3.1Assume that f = f ( v ) is a smooth function. Suppose that f ( v ) = O ( | v | 1 + θ ) ( θ 1 is an integer) when | v | ν 0 . Then for integer m 0 , if v W m , q ( R n ) L p ( R n ) L ( R n ) and v L ν 0 , the following inequalities hold:

x m f ( v ) L r C v L p x m v L q v L θ 1 , (3.1)

where 1 r = 1 p + 1 q , 1 p , q , r + .

Theorem 3.1Let s 0 . Suppose that f ( v ) is a smooth function and f ( 0 ) > α . Assume that u 0 H s + 2 L 1 , u 1 H s L 1 . Put

E 0 : = u 0 H s + 2 L 1 + u 1 H s L 1 .

Then there exists a positive constant δ 1 such that if E 0 δ 0 , and the initial value problem (1.1), (1.2) has a unique global solution u ( x , t ) satisfying

u C 0 ( [ 0 , ) ; H s + 2 ( R ) ) C 1 ( [ 0 , ) ; H s ( R ) ) .

Moreover, the solution satisfies the decay estimate

x k u ( t ) L 2 C E 0 ( 1 + t ) 1 4 k 2 (3.2)

and

x k u t ( t ) L 2 C E 0 ( 1 + t ) 5 4 k 2 , (3.3)

where 0 k s + 2 in (3.2) and 0 k s in (3.3).

Proof The existence and uniqueness of small solutions can be proved by the contraction mapping principle. Here we only show the decay estimates (3.2) and (3.3) for the solution u of (2.12) satisfying u ( t ) L M 0 with some M 0 . To this end, we introduce the quantity

X ( t ) = k = 0 s + 2 sup 0 τ t ( 1 + τ ) 1 4 + k 2 x k u ( τ ) L 2 . (3.4)

Here we note that

u x ( t ) L C X ( t ) ( 1 + t ) 1 , (3.5)

provided that s 0 . This follows from the Gagliardo-Nirenberg inequality, u x L C u x x L 2 1 2 u x L 2 1 2 and the definition of X ( t ) in (3.4).

Applying x k to (2.12) and taking the L 2 norm, we obtain

x k u ( t ) L 2 x k G ( t ) ( u 0 + u 1 ) L 2 + x k H ( t ) u 0 L 2 + 0 t G ( t τ ) x g ( u x ( τ ) ) L 2 d τ = : I 1 + I 2 + J . (3.6)

Firstly, we estimate I 1 . We apply (2.26) with p = 1 , j = 0 and l = 0 ( l = 2 for k = 0 , l = 1 for k = 1 ). This yields

I 1 C ( 1 + t ) 1 4 k 2 u 0 + u 1 L 1 + C e c t x ( k 2 ) + ( u 0 + u 1 ) L 2 C E 0 ( 1 + t ) 1 4 k 2 , (3.7)

where ( k 2 ) + = max { k 2 , 0 } . Similarly, applying (2.27) with p = 1 , j = 0 , and l = 0 to the term I 2 , we have

I 2 C ( 1 + t ) 1 4 k 2 u 0 L 1 + C e c t x k u 0 L 2 C E 0 ( 1 + t ) 1 4 k 2 . (3.8)

We estimate the nonlinear term J. We divide J into two parts and write J = J 1 + J 2 , where J 1 and J 2 are corresponding to the time intervals [ 0 , t / 2 ] and [ t / 2 , t ] , respectively. For the term J 1 , we apply (2.30) with p = 1 , j = 0 , and l = 0 . This yields

J 1 C 0 t / 2 ( 1 + t τ ) 3 4 k 2 g ( u x ) ( τ ) L 1 d τ + C 0 t / 2 e c ( t τ ) x ( k 1 ) + g ( u x ) ( τ ) L 2 d τ = : J 11 + J 12 . (3.9)

Here we see that g ( u x ) L 1 C u x L 2 2 by Lemma 3.1. Thus we have g ( u x ) ( τ ) L 1 C X ( t ) 2 ( 1 + τ ) 3 2 . Therefore we can estimate the term J 11 as

J 11 C X ( t ) 2 0 t / 2 ( 1 + t τ ) 3 4 k 2 ( 1 + τ ) 3 2 d τ C X ( t ) 2 ( 1 + t ) 3 4 k 2 0 t / 2 ( 1 + τ ) 3 2 d τ C X ( t ) 2 ( 1 + t ) 3 4 k 2 .

On the other hand, we have x ( k 1 ) + g ( u x ) L 2 C u x L x k u L 2 by Lemma 3.1. Therefore, using (3.5), we find that x ( k 1 ) + g ( u x ) L 2 C X ( t ) 2 ( 1 + τ ) 5 4 k 2 . Consequently, we can estimate the term J 12 as

J 12 C X ( t ) 2 0 t / 2 e c ( t τ ) ( 1 + τ ) 5 4 k 2 d τ C X ( t ) 2 e c t .

Finally, we estimate the term J 2 on the time interval [ t / 2 , t ] . Applying (2.30) with p = 2 , j = k , and l = 0 , and using x ( k 1 ) + g ( u x ) ( τ ) L 2 C X ( t ) 2 ( 1 + τ ) 5 4 k 2 , we obtain

J 2 C t / 2 t ( 1 + t τ ) 1 x ( k 1 ) + g ( u x ) ( τ ) L 2 d τ + C t / 2 t e c ( t τ ) x ( k 1 ) + g ( u x ) ( τ ) L 2 d τ C X ( t ) 2 t / 2 t ( 1 + t τ ) 1 ( 1 + τ ) 5 4 k 2 d τ C X ( t ) 2 ( 1 + t ) 5 4 k 2 + ε . (3.10)

Thus we have shown that

J C X ( t ) 2 ( 1 + t ) 3 4 k 2 .

Substituting all these estimates into (3.6), we obtain

( 1 + t ) 1 4 + k 2 x k u ( t ) L 2 C E 0 + C X ( t ) 2

for 0 k s + 2 . Consequently, we have X ( t ) C E 0 + C X ( t ) 2 , from which we can deduce X ( t ) C E 0 , provided that E 0 is suitably small. This proves the decay estimate (3.2).

In what follows, we prove the decay estimate (3.3) for the time derivative u t . For this purpose we differentiate (2.12) with respect to t to obtain

u t ( t ) = G t ( t ) ( u 0 + u 1 ) + H t ( t ) u 0 + 0 t G t ( t τ ) x g ( u x ) ( τ ) d τ . (3.11)

Applying x k to (3.11) and taking the L 2 norm, we have

x k u t ( t ) L 2 x k G t ( t ) ( u 0 + u 1 ) L 2 + x k H t ( t ) u 0 L 2 + 0 t x k G t ( t τ ) x g ( u x ) ( τ ) L 2 d τ = : I 1 + I 2 + J , (3.12)

where 0 k s . For the term I 1 , we apply (2.28) with p = 1 , j = 0 , and l = 0 to get

I 1 C ( 1 + t ) 5 4 k 2 u 0 + u 1 L 1 + e c t x k ( u 1 + u 0 ) L 2 C E 0 ( 1 + t ) 5 4 k 2 .

Also, for the term I 2 , applying (2.29) with p = 1 , j = 0 and l = 0 , we have

I 2 C ( 1 + t ) 9 4 k 2 u 0 L 1 + e c t x k + 2 u 0 L 2 C E 0 ( 1 + t ) 9 4 k 2 .

To estimate the nonlinear term J , we rewrite J = J 1 + J 2 , where J 1 and J 2 correspond to the time intervals [ 0 , t / 2 ] and [ t / 2 , t ] , respectively. For the term J 1 , we apply (2.31) with p = 1 , j = 0 , and l = 0 . This yields

J 1 C 0 t / 2 ( 1 + t τ ) 7 4 k 2 g ( u x ) ( τ ) L 1 d τ + C 0 t / 2 e c ( t τ ) x k + 1 g ( u x ) ( τ ) L 2 d τ = : J 11 + J 12 .

Since g ( u x ) ( τ ) L 1 C X ( t ) 2 ( 1 + τ ) 3 2 as before, we can estimate the term J 11 as

J 11 C X ( t ) 2 0 t / 2 ( 1 + t τ ) 7 4 k 2 ( 1 + τ ) 3 2 d τ C X ( t ) 2 ( 1 + t ) 7 4 k 2 .

Also, the term J 12 is estimated similarly as before and we have J 12 C X ( t ) 2 e c t . Finally, we estimate the term J 2 by applying (2.31) with p = 2 , j = k + 2 , and l = 0 and obtain

J 2 C t / 2 t ( 1 + t τ ) 1 x k + 1 g ( u x ) ( τ ) L 2 d τ + C t / 2 t e c ( t τ ) x k + 1 g ( u x ) ( τ ) L 2 d τ C X ( t ) 2 t / 2 t ( 1 + t τ ) 1 ( 1 + τ ) 9 4 k 2 d τ C X ( t ) 2 ( 1 + t ) 9 4 k 2 + ε ,

where we used the estimate x k + 1 g ( u x ) ( τ ) L 2 C X ( t ) 2 ( 1 + τ ) 9 4 k 2 . Consequently we have shown that

J C X ( t ) 2 ( 1 + t ) 9 4 k 2 .

Substituting all these estimates together with the previous estimate X ( t ) C E 0 into (3.12), we arrive at the desired estimate (3.3) for 0 k s . This completes the proof of Theorem 3.1. □

The above proof of Theorem 3.1 shows that the solution u to the integral equation (2.12) is asymptotic to the linear solution u L given by the formula u L ( t ) : = G ( t ) ( u 0 + u 1 ) + H ( t ) u 0 in (2.11) as t . This result is stated as follows.

Corollary 3.1Assume the same conditions of Theorem 3.1. Then the solutionuof the problem (1.1), (1.2), which is constructed in Theorem 3.1, can be approximated by the solution u L to the linearized problem (2.1), (1.2) as t . More precisely, we have the following asymptotic relations:

x k ( u u L ) ( t ) L 2 C E 0 2 ( 1 + t ) 3 4 k 2 , x k ( u u L ) t ( t ) L 2 C E 0 2 ( 1 + t ) 7 4 k 2

for 0 k s + 2 and 0 k s , respectively, where u L ( t ) : = G ( t ) ( u 0 + u 1 ) + H ( t ) u 0 is the linear solution.

4 Asymptotic profile

In this section, our aim is to establish an asymptotic profile to our global solution that is constructed in Theorem 4.1. In the previous section, we have shown that the solution u to the problem (1.1), (1.2) can be approximated by the linear solution u L . In what follows, we shall derive a simpler asymptotic profile of the linear solution u L .

Let v be the solution to the initial data problem

v t β v x x = 0 , t = 0 : v = ( u 0 + u 1 ) ( x ) . (4.1)

Then

v = G 0 ( x , t ) ( u 0 + u 1 ) ( x ) (4.2)

gives a asymptotic profile of the linear solution u L . In fact we have the following.

Lemma 4.1Let s 0 . Assume that u 0 , u 1 H s + 2 L 1 and put E 0 = ( u 0 , u 1 ) H s + 2 L 1 . Let u L be the linear solution and letvbe defined by (4.2). Then we have

x k ( u L v ) ( t ) L 2 C E 0 ( 1 + t ) 5 4 k 2 (4.3)

for 0 k s + 2 .

Proof Note that ( u L v ) ( t ) = ( G G 0 ) ( t ) ( u 0 + u 1 ) + H ( t ) u 0 , so for the proof of (4.3), it suffices to show the following estimates:

x k ( G G 0 ) ( t ) ϕ L 2 C ( 1 + t ) 5 4 ϕ L 1 + C e c t x k ϕ L 2 , x k H ( t ) ϕ L 2 C ( 1 + t ) 5 4 ϕ L 1 + C e c t x k ϕ L 2 ,

where 0 k s + 2 . These estimates can be obtained by (2.27) and (2.32). Here we omit the details. □

When M 0 = R ( u 0 + u 1 ) ( x ) d x 0 , we call

v ¯ ( x , t ) = M 0 G 0 ( x , t + 1 ) (4.4)

the diffusion wave with the amount M 0 . Obviously, v ¯ ( x , t ) satisfies the following problem:

{ v ¯ t β v ¯ x x = 0 , t = 0 : v ¯ = M 0 G 0 ( x , 1 ) . (4.5)

Therefore, v v ¯ satisfies

{ ( v v ¯ ) t β ( v v ¯ ) x x = 0 , t = 0 : v v ¯ = u 0 + u 1 M 0 G 0 ( x , 1 ) . (4.6)

It is not difficult to prove the following lemma.

Lemma 4.2Assume that x ( u 0 ( y ) + u 1 ( y ) M 0 G 0 ( y , 1 ) ) d y L 1 , then

x k ( v v ¯ ) ( t ) L 2 C ( 1 + t ) 3 4 k 2 . (4.7)

Combining Corollary 3.1, Lemma 4.1, and Lemma 4.2, we immediately have the following.

Theorem 4.1Under the same assumption as Theorem 3.1, and also assuming that M 0 = R ( u 0 + u 1 ) ( x ) d x 0 and x ( u 0 ( y ) + u 1 ( y ) M 0 G 0 ( y , 1 ) ) d y L 1 , we letube the global solution to the problem (1.1), (1.2), which is constructed in Theorem 3.1 and we let v ¯ be the diffusion wave defined by (4.4). Then we have

x k ( u v ¯ ) ( t ) L 2 C ( 1 + t ) 3 4 k 2 . (4.8)

We have M 0 = ( u 0 + u 1 ) ( x ) d x = 0 and M 1 = U 0 ( x ) d x 0 , where U 0 ( x ) : = x ( u 0 + u 1 ) ( y ) d y . We consider the initial value problem

{ ρ t β ρ x x = 0 , t = 0 : ρ = U 0 ( x ) (4.9)

and

{ ρ ¯ t β ρ ¯ x x = 0 , t = 0 : ρ ¯ = M 1 G 0 ( x , 1 ) . (4.10)

When x ( U 0 ( y ) M 0 G 0 ( y , 1 ) ) d y L 1 , by applying Lemma 4.2, we have

x k ( ρ ρ ¯ ) ( t ) L 2 C ( 1 + t ) 3 4 k 2 . (4.11)

We call

v ˜ ( x , t ) = M 1 G 0 x ( x , t + 1 ) (4.12)

the diffusion wave with the amount M 1 . Obviously, v v ˜ satisfies

v v ˜ = x ρ x ρ ¯ . (4.13)

From (4.11) and (4.13), we immediately obtain the following lemma.

Lemma 4.3Assume that x ( U 0 ( y ) M 1 G 0 ( y , 1 ) ) d y L 1 , then

x k ( v v ˜ ) ( t ) L 2 C ( 1 + t ) 5 4 k 2 . (4.14)

Corollary 3.1, Lemma 4.1, and Lemma 4.3 immediately give the following result.

Theorem 4.2Under the same assumption as Theorem 3.1, also assuming that M 0 = R ( u 0 + u 1 ) ( x ) d x = 0 , M 1 = R U 0 ( x ) d x 0 , and x ( U 0 ( y ) M 1 G 0 ( y , 1 ) ) d y L 1 , where U 0 ( x ) : = x ( u 0 + u 1 ) ( y ) d y , we letube the global solution to the problem (1.1), (1.2), which is constructed in Theorem 3.1 and we let v ˜ be the diffusion wave defined by (4.12). Then we have

x k ( u v ˜ ) ( t ) L 2 C ( 1 + t ) 5 4 k 2 . (4.15)

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed to each part of this work equally and read and approved the final manuscript.

References

  1. Takeda, H, Yoshikawa, S: On the initial value problem of the semilinear beam equation with weak damping I: smoothing effect. J. Math. Anal. Appl.. 401, 244–258 (2013). Publisher Full Text OpenURL

  2. Takeda, H, Yoshikawa, S: On the initial data of the semilinear beam equation with weak damping II: asymptotic profiles. J. Differ. Equ.. 253, 3061–3080 (2012). Publisher Full Text OpenURL

  3. Dai, W, Kong, D: Global existence and asymptotic behavior of classical solutions of quasilinear hyperbolic systems with linearly degenerate characteristic fields. J. Differ. Equ.. 235, 127–165 (2007). Publisher Full Text OpenURL

  4. Kong, D, Yang, T: Asymptotic behavior of global classical solutions of quasilinear hyperbolic systems. Commun. Partial Differ. Equ.. 28, 1203–1220 (2003). Publisher Full Text OpenURL

  5. Nakao, M, Ono, K: Existence of global solutions to the Cauchy problem for the semilinear dissipative wave equations. Math. Z.. 214, 325–342 (1993). Publisher Full Text OpenURL

  6. Nishihara, K: Asymptotic behavior of solutions of quasilinear hyperbolic equations with damping. J. Differ. Equ.. 137, 384–395 (1997). Publisher Full Text OpenURL

  7. Ono, K: Global existence and asymptotic behavior of small solutions for semilinear dissipative wave equations. Discrete Contin. Dyn. Syst.. 9, 651–662 (2003)

  8. Takeda, H, Yoshikawa, S: Asymptotic profile of solutions for the limit unstable Cahn-Hilliard equation with inertial term. Differ. Integral Equ.. 25, 341–362 (2012)

  9. Geng, S: Convergence rates to asymptotic profile for solutions of quasilinear hyperbolic equations with linear damping. J. Hyperbolic Differ. Equ.. 8, 115–129 (2011). Publisher Full Text OpenURL

  10. Wang, Y, Wei, Z: Global existence and asymptotic behavior of solutions to Cahn-Hilliard equation with inertial term. Int. J. Math.. 23(9), Article ID 1250087 (2012)

  11. Yang, Z: Longtime behavior of the Kirchhoff type equation with strong damping on R n . J. Differ. Equ.. 242, 269–286 (2007). Publisher Full Text OpenURL

  12. Sugitani, Y, Kawashima, S: Decay estimates of solution to a semi-linear dissipative plate equation. J. Hyperbolic Differ. Equ.. 7, 471–501 (2010). Publisher Full Text OpenURL

  13. Wang, Y, Liu, F, Zhang, Y: Global existence and asymptotic of solutions for a semi-linear wave equation. J. Math. Anal. Appl.. 385, 836–853 (2012). Publisher Full Text OpenURL

  14. Kato, M, Kawashima, S: Asymptotic behavior of solutions to the generalized cubic double dispersion equation in one space dimension. Kinet. Relat. Models. 6, 969–987 (2013)

  15. Wang, S, Xu, H: On the asymptotic behavior of solution for the generalized IBq equation with hydrodynamical damped term. J. Differ. Equ.. 252, 4243–4258 (2012). Publisher Full Text OpenURL

  16. Zhang, Z, Zhang, Y: Global existence and asymptotic behavior of solutions to a class of fourth-order wave equations. Bound. Value Probl.. 2013, Article ID 168 (2013)

  17. Li, T, Chen, Y: Nonlinear Evolution Equations, Scientific Press, Beijing (1989) (in Chinese)

  18. Zheng, S: Nonlinear Evolution Equations, Chapman & Hall/CRC, Boca Raton (2004)