Skip to main content

Global existence and asymptotic behavior of solutions for the double dispersive-dissipative wave equation with nonlinear damping and source terms

Abstract

In this paper, we consider the initial boundary value problem of the double dispersive-dissipative wave equation with nonlinear damping and source terms. By the combination of the Galerkin method and the monotonicity-compactness method, the existence of global solutions is obtained with the least amount of a priori estimates. Moreover, the asymptotic behavior of global solutions is investigated under some assumptions on the initial data.

1 Introduction

This paper deals with the initial boundary value problem of the double dispersive-dissipative wave equation with nonlinear damping and source terms

$$\begin{aligned} u_{tt}-\triangle{u}-\beta\triangle{u}_{t}+ \gamma\triangle ^{2}u-\delta\triangle{u}_{tt}+g(u_{t})=f(u),\quad x\in\Omega, t>0, \end{aligned}$$
(1.1)

where \(\Omega\subset\mathrm{R}^{n}\) is a suitably smooth bounded domain, β, γ, and δ are some physical parameters. \(g(s)\), \(f(s)\) are given nonlinear functions.

Equation (1.1) includes many important physical models. For example, in the absence of a dissipative term, double dispersive terms (\(\beta =\gamma=\delta=0\)), and a damping term \(g(u_{t})=0\), the model reduces to the common semilinear wave equation

$$\begin{aligned} u_{tt}-\triangle{u}=f(u),\quad x\in\Omega, t>0. \end{aligned}$$
(1.2)

In 1968, Sattinger [1] studied the existence of global weak solutions of (1.2) by using the potential well method. From then on, the potential well theory has become one of the most important methods for studying nonlinear evolution equations. In 1975, Payne and Sattinger [2] made the most important and typical work on the potential well method. They proved the finite time blow up of solutions of (1.2). In 2003, Liu [3] studied the initial boundary value problem of (1.2), where \(f(u)=\vert u\vert ^{p-1}u\). He gave some results on the properties of the family of potential wells. Then, by using them, he obtained some results of the existence and nonexistence of global solutions. In a subsequent article [4], Liu gave a threshold result of the global existence and nonexistence of solutions and proved the global existence of solutions with critical initial conditions \(J(u_{0})=d\). For the semilinear wave equation (1.2) with nonlinear boundary velocity feedback, it is important to cite [5–7] and the references therein. In 1993, Lasiecka and Tataru [5] obtained the uniform decay rates of the energy for (1.2) under the assumption that the boundary velocity feedback is dissipative. In 2002, Vitillaro [6] studied the local existence, blow up, and global existence of the solutions for (1.2) with nonlinear boundary velocity feedback. In 2004, Cavalcanti et al. [7] proved the existence and uniform decay rates of the energy even if the nonlinear boundary velocity feedback has a polynomial growth near the origin.

In the absence of a dissipative term, double dispersive terms (\(\beta =\gamma=\delta=0\)), taking the nonlinear damping term \(g(u_{t})=a\vert u_{t}\vert ^{m-1}u_{t}\) and the nonlinear source term \(f(u)=b\vert u\vert ^{p-1}u\), the model (1.1) reduces to the wave equation with nonlinear damping and source terms

$$\begin{aligned} u_{tt}-\triangle{u}+a\vert u_{t}\vert ^{m-1}u_{t}=b\vert u\vert ^{p-1}u,\quad x\in\Omega, t>0. \end{aligned}$$
(1.3)

In 1994, Georgiev and Todorova [8] investigated the initial boundary value problem of (1.3), where \(a=b=1\), \(p,m>1\). They proved the existence of global solutions under the condition \(1< p\leq {m}\). When \(p\geq{m}>1\), they obtained the finite time blow up of solutions for sufficiently large initial data. In 1996, Ikehata [9] investigated the initial boundary value problem of (1.3), where \(b=1\) and any \(a=\delta>0\). He proved that \(1\leq{m}< p<\infty\) if \(n=1,2\), and \(1\leq{m}< p\leq\frac{n}{n-2}\) if \(n\geq3\), the problem (1.3) had a global solution for sufficiently small initial data. In 1999, Vitillaro [10] proved the global nonexistence of the solutions with positive initial energy. He also gave the application concerning the classical equations of linear elasticity and the damped clamped plate equation. In 2001, Messaoudi [11] changed the nonlinear term \(a\vert u_{t}\vert ^{m-1}u_{t}\), \(b\vert u\vert ^{p-1}u\) into \(a\vert u_{t}\vert ^{m-2}u_{t}\), \(b\vert u\vert ^{p-2}u\). He proved that any strong solution, with negative initial energy, blows up in finite time under the condition \(p>m\). Messaoudi [12, 13] changed the nonlinear term \(a\vert u_{t}\vert ^{m-1}u_{t}\), \(b\vert u\vert ^{p-1}u\) into \(au_{t}(1+\vert u_{t}\vert ^{m-2})\), \(b\vert u\vert ^{p-2}u\). He investigated the global nonexistence and exponential decay of solutions, respectively. For the Cauchy problem associated to (1.3), when Ω is replaced by the entire space \(\mathrm{R}^{n}\), Serrin and Todorova [14] studied the existence of the solutions for the case \(p< m\).

When taking into account no effects of dispersion caused by transverse shearing and nonlinear damping, in one dimension, the model (1.1) reduces to the fourth-order wave equation

$$\begin{aligned} u_{tt}-u_{xx}-u_{xxt}-u_{xxtt}=f(u). \end{aligned}$$
(1.4)

This type of equation was derived by Hayes and Saccomandi [15] in studying the propagation of transverse homogeneous waves in special incompressible viscoelastic solids. In 2000, Shang [16] studied the initial boundary value problem of the fourth-order nonlinear wave equation (1.4). The existence and uniqueness of a global strong solution for the problem were obtained by means of the Galerkin method. The asymptotic behavior and blow up phenomenon of the solution for the problem were discussed under certain conditions.

In the absence of a dispersive term and dissipation (\(\beta=\gamma =0\)), and when the damping term is linear, that is, \(g(u_{t})=u_{t}\), the model (1.1) becomes the 1D fourth-order wave equation

$$\begin{aligned} u_{tt}-u_{xx}+u_{t}-u_{xxtt}=f(u). \end{aligned}$$
(1.5)

This type of equation was derived in the study of a longitudinal vibration of a bar. The weak damping term of (1.5) is introduced to model the contacting of the bar with a rough substrate or a viscous external medium. The dispersive term of the above equation is used to explain the lateral inertia of a bar [17].

In the absence of the dispersive term and the damping term, that is, \(g(u_{t})=0\) and \(\gamma=0\), the model (1.1) reduces to the fourth-order wave equation (\(n\geq1\))

$$\begin{aligned} u_{tt}-\triangle{u}-\triangle{u}_{t}- \triangle{u}_{tt}=f(u),\quad x\in \Omega, t>0. \end{aligned}$$
(1.6)

In 2000, Shang [18] investigated the existence, uniqueness, asymptotic behavior, and the blow up phenomenon of the solutions under some specific assumptions on f and for \(n=1,2,3\). In 2004, Zhang and Hu [19] proved the existence and the stability of global weak solutions. In 2007, Xie and Zhong [20] studied the initial boundary value problem of (1.6). They established the existence of global attractors in \(H_{0}^{1}(\Omega)\times{H}_{0}^{1}(\Omega)\), where the nonlinear term f satisfies a critical exponential growth condition. In 2008, Xu et al. [21] investigated the asymptotic behavior of solutions for (1.6) by using the multiplier method. For more results on the long-time behaviors of global strong solutions of the initial boundary value problem of (1.6), the reader is referred to [22–25].

If \(\gamma=0\), \(\beta=\delta=1\), \(g(u_{t})={u}_{t}\) and \(f(u)=\vert u\vert ^{p-1}u\), the model (1.1) reduces to the fourth-order dispersive-dissipative wave equation

$$\begin{aligned} u_{tt}-\triangle{u}+\triangle{u}_{t}-\triangle {u}_{tt}+{u}_{t}=\vert u\vert ^{p-2}u,\quad x\in \Omega, t>0, \end{aligned}$$
(1.7)

In 2012, Xu and Yang [26] investigated the initial boundary value problem of (1.7), where \(1< p<\infty\) if \(n=1,2\), and \(1< p\leq\frac {n+2}{n-2}\) if \(n\geq3\). A blow up result for certain solutions with arbitrary positive initial energy was established.

As far as we know, there have been no results up till now on the initial boundary value problem for a nonlinear wave equation with double dispersive terms \(\triangle^{2}u\), \(\triangle{u}_{tt}\), the strong dissipation term \(\triangle{u}_{t}\), the nonlinear damping term \(g(u_{t})\), and the nonlinear source term \(f(u)\). So the aim of the present paper is to solve this open problem.

In this work, we investigate the initial boundary value problem for the double dispersive wave equation with a strong dissipation term, a nonlinear damping and source terms:

$$\begin{aligned} \begin{cases} u_{tt}-\triangle{u}-\triangle{u}_{t}+\triangle^{2}u-\triangle {u}_{tt}+a\vert u_{t}\vert ^{m-2}u_{t}=b\vert {u}\vert ^{p-2}{u},\quad x\in\Omega, t>0,\\ u(x,0)=u_{0}(x),\qquad u_{t}(x,0)=u_{1}(x),\quad x\in\Omega,\\ u=0, \qquad \frac{\partial{u}}{\partial\nu}=0, \quad x\in\partial\Omega, t>0, \end{cases} \end{aligned}$$
(1.8)

where \(a,b>0\), \(p,m>2\), ν represents the unit outward normal to ∂Ω, and Ω is a bounded domain of \({R}^{n}\) (\(n\geq1\)) with a smooth boundary ∂Ω. First of all, by the combination of the Galerkin method and the monotonicity-compactness method, the existence of global solutions is obtained with the least amount of a priori estimates. Moreover, the asymptotic behavior of global solutions is investigated under some assumptions on the initial data.

This paper is organized as follows. In Section 2, we introduce some notation, basic definitions, and important lemmas for proving the main theorem. In Section 3, the existence of global weak solutions is proved by the Galerkin method and the monotonicity-compactness method. In Section 4, we consider the asymptotic behavior of global solutions under some assumptions on the initial data.

2 Preliminaries

In this section, we introduce some notation, basic definitions, and important lemmas which will be needed in this paper.

For functions \(u(x,t)\), \(v(x,t)\) defined on Ω, we introduce

$$\begin{aligned}& (u,v)=\int_{\Omega}uv\,dx, \qquad \Vert u\Vert _{{2}}=\biggl( \int_{\Omega }\vert u\vert ^{2}\,dx \biggr)^{\frac{1}{2}}, \\& \Vert u\Vert _{{p}}=\biggl(\int_{\Omega} \vert u \vert ^{p}\,dx\biggr)^{\frac{1}{p}},\qquad \Vert u\Vert _{H^{m}}=\biggl(\sum_{\vert \alpha \vert \leq{m}}\bigl\Vert D^{\alpha}u\bigr\Vert _{{2}}^{2}\biggr)^{\frac{1}{2}}, \\& \Vert u\Vert _{{\infty}}=\operatorname {ess}\sup_{x\in\Omega}\bigl\vert u(x) \bigr\vert . \end{aligned}$$

To obtain the results of this paper, we will introduce the energy function

$$\begin{aligned} E(t)=\frac{1}{2}\Vert {u}_{t}\Vert _{2}^{2}+\frac{1}{2}\Vert \nabla{u} \Vert _{2}^{2}+\frac{1}{2}\Vert \triangle{u} \Vert _{2}^{2}+\frac{1}{2}\Vert \nabla {u}_{t} \Vert _{2}^{2}-\frac{b}{p}\Vert {u}\Vert _{p}^{p}. \end{aligned}$$
(2.1)

For simplicity, we define the weak solutions of (1.8) over the interval \([0,T)\), but it is to be understood throughout that T is either infinity or the limit of the existence interval.

Definition 2.1

We say that \(u(x,t)\) is a weak solution of the problem (1.8) on the interval \(\Omega\times[0,T)\), if \(u\in {L}^{\infty}(0,T;H_{0}^{2}(\Omega)\cap{L}^{p}(\Omega))\), \(u_{t}\in {L}^{\infty}(0,T;H_{0}^{1}(\Omega))\cap{L}^{m}(Q_{T})\) satisfy the following conditions:

  1. (i)

    For any \(v\in{H}_{0}^{2}(\Omega)\cap{L}^{m}(\Omega)\) and a.e. \(0\leq{t}< T\), such that

    $$\begin{aligned} & (u_{tt},v)+(\nabla{u},\nabla{v})+(\nabla{u}_{t}, \nabla {v})+(\triangle{u},\triangle{v})+(\nabla{u}_{tt},\nabla{v}) \\ &\quad {}+\bigl(a\vert u_{t}\vert ^{m-2}u_{t},v \bigr)=\bigl(b\vert {u}\vert ^{p-2}{u},v\bigr); \end{aligned}$$
    (2.2)
  2. (ii)

    \(u(x,0)=u_{0}(x)\) in \({H}_{0}^{2}(\Omega)\cap{L}^{p}(\Omega )\), \(u_{t}(x,0)=u_{1}\) in \({H}_{0}^{1}(\Omega)\cap{L}^{m}(\Omega)\).

Lemma 2.1

Let \(p,m>2\), \(a,b>0\), and u be a solution of (1.8). Then \(E(t)\) is a non-increasing function, that is,

$$\begin{aligned} E'(t)\leq{0}. \end{aligned}$$
(2.3)

Moreover, the following energy inequality holds:

$$\begin{aligned} E(t)+\int_{0}^{t}\bigl\Vert \nabla{u}_{t}(\tau)\bigr\Vert _{2}^{2}\,d\tau \leq{E}(0). \end{aligned}$$
(2.4)

Proof

Multiplying (1.8) by \(u_{t}\), and integrating over Ω, using integration by parts and some manipulation as in [29], we obtain (2.3), (2.4) for any regular solutions. These results remain valid for weak solutions by a simple density argument. We also refer the reader to [5, 7] as regards the method of density arguments.

The following lemma is similar to Lemma 6.1 of Chapter ii of [27] with a slight modification. □

Lemma 2.2

Assume that the function u satisfies \(u\in{L}^{\infty }(0,T;H_{0}^{2}(\Omega))\), \(u_{t}\in{L}^{\infty }(0,T; H_{0}^{1}(\Omega))\cap{L}^{m}(Q_{T})\), \(u(x,0)=u_{0}\), \(u_{t}(x,0)=u_{1}\), and further assume that

$$\begin{aligned} & u_{tt}-\triangle{u}-\triangle{u}_{t}+\triangle^{2}u- \triangle{u}_{tt}=g, \\ &\quad g\in{L}^{2}\bigl(0,T;H^{-1}(\Omega)\bigr)+{L}^{m'}(Q_{T}), \end{aligned}$$
(2.5)

then for any \(t\in(0,T)\), the following inequality holds:

$$\begin{aligned} &\Vert \nabla{u} \Vert _{2}^{2}+\Vert {u}_{t}\Vert _{2}^{2}+\Vert \triangle{u} \Vert _{2}^{2}+\Vert \nabla{u}_{t}\Vert _{2}^{2}+2\int_{0}^{t}\Vert \nabla{u}_{t}\Vert _{2}^{2}\,d\tau \\ &\quad \geq \Vert \nabla{u}_{0}\Vert _{2}^{2}+ \Vert {u}_{1}\Vert _{2}^{2}+\Vert \triangle {u}_{0}\Vert _{2}^{2}+\Vert \nabla{u}_{1}\Vert _{2}^{2}+2\int _{0}^{t}(g,u_{t})\,d\tau. \end{aligned}$$
(2.6)

Proof

By choosing the continuous piecewise linear function \(\theta_{m}(t)\) and the regular sequence \(\rho_{n}(t)\), we construct the smooth function \(v(t)=((\theta_{m}u')\ast\rho_{n}\ast\rho _{n})\theta_{m}\). Then, multiplying (2.5) by \(v(t)\) and integrating over \(Q=\Omega\times(0,T)\), using integration by parts and some manipulation as in Chapter 2 of [27] as \(m,n\rightarrow\infty\), we can obtain the inequality (2.6). This method of smooth approximations is called the method of mollifiers. We also refer reader to [27, 28] about the method of mollifiers.

We construct an approximate weak solution of the problem (1.8) by the Galerkin method. Let \(\{w_{j}\}\) be the system of base functions of \(H_{0}^{2}(\Omega)\cap{L}^{m}(\Omega)\).

Now suppose that the approximate weak solutions of the problem (1.8) can be written

$$\begin{aligned} u_{l}(x,t)=\sum_{j=1}^{l}d_{l}^{j}(t)w_{j}(x),\quad l=1,2,\ldots. \end{aligned}$$
(2.7)

According to the Galerkin method, these coefficients \(d_{l}^{j}(t)\) need to satisfy the following initial value problem of the nonlinear ordinary differential equations:

$$\begin{aligned} \begin{cases} (u_{ltt},w_{j})+(\nabla{u}_{l},\nabla{w}_{j})+(\nabla{u}_{lt},\nabla {w}_{j})+(\triangle{u}_{l},\triangle{w}_{j})+(\nabla{u}_{ltt},\nabla {w}_{j})\\ \quad {}+(a\vert u_{lt}\vert ^{m-2}u_{lt},w_{j})=(b\vert {u}_{l}\vert ^{p-2}{u}_{l},w_{j}), \quad x\in \Omega, t>0,\\ u_{l}(x,0)=\sum_{j=1}^{l}d_{l}^{j}(0)w_{j}(x)\longrightarrow {u}_{0}(x), \quad \mbox{in }H_{0}^{2}(\Omega)\cap{L}^{p}(\Omega),\\ u_{lt}(x,0)=\sum_{j=1}^{l}d_{l}^{j}(0)'w_{j}(x)\longrightarrow {u}_{1}(x), \quad \mbox{in }H_{0}^{1}(\Omega)\cap{L}^{m}(\Omega). \end{cases} \end{aligned}$$
(2.8)

Under some assumption for the nonlinear terms and a priori estimates in Section 3, we prove that the initial value problem (2.8) of the nonlinear ordinary differential equations has global solutions in the interval \([0,T]\). Furthermore, we show that the solutions of the problem (1.8) can be approximated by the functions \(u_{l}(x,t)\). □

3 The existence of global weak solutions

In this section, we establish the existence of global weak solutions for the problem (1.8).

Theorem 3.1

Suppose that \(a,b>0\), \(2< p\leq{m}<\infty\), \({u}_{0}(x)\in {H}_{0}^{2}(\Omega)\cap{L}^{p}(\Omega)\), \({u}_{1}(x)\in {H}_{0}^{1}(\Omega)\cap{L}^{m}(\Omega)\), Then for any \(T>0\), the problem (1.8) admits a global weak solution \(u\in{L}^{\infty }(0,T;H_{0}^{2}(\Omega)\cap{L}^{p}(\Omega))\), \(u_{t}\in{L}^{\infty }(0,T;H_{0}^{1}(\Omega))\cap{L}^{m}(Q_{T})\).

Proof

Multiplying (2.8) by \(d_{l}^{j}(t)'\) and summing for \(j=1,\ldots,l\), we have

$$\begin{aligned} &\frac{d}{dt} \biggl[\frac{1}{2}\Vert u_{lt}\Vert _{2}^{2}+\frac{1}{2}\Vert \nabla{u}_{l}\Vert _{2}^{2}+\frac{1}{2} \Vert \triangle{u}_{l}\Vert _{2}^{2}+ \frac{1}{2}\Vert \nabla{u}_{lt}\Vert _{2}^{2}- \frac{b}{p}\Vert u_{l}\Vert _{p}^{p} \biggr] \\ &\quad {}+\Vert \nabla{u}_{lt}\Vert _{2}^{2}+a \Vert u_{lt}\Vert _{m}^{m}=0. \end{aligned}$$
(3.1)

Integrating (3.1) with respect to t from 0 to t, we find

$$\begin{aligned} &\frac{1}{2}\Vert u_{lt}\Vert _{2}^{2}+\frac{1}{2}\Vert \nabla{u}_{l} \Vert _{2}^{2}+\frac{1}{2}\Vert \triangle{u}_{l}\Vert _{2}^{2}+\frac{1}{2} \Vert \nabla {u}_{lt}\Vert _{2}^{2} \\ &\qquad {} -\frac{b}{p}\Vert u_{l}\Vert _{p}^{p}+ \int_{0}^{t}\Vert \nabla{u}_{lt} \Vert _{2}^{2}\,d\tau+a\int_{0}^{t} \Vert u_{lt}\Vert _{m}^{m}\,d\tau \\ & \quad =\frac{1}{2}\bigl\Vert u_{lt}(0)\bigr\Vert _{2}^{2}+\frac{1}{2}\bigl\Vert \nabla{u}_{l}(0) \bigr\Vert _{2}^{2}+\frac{1}{2}\bigl\Vert \triangle{u}_{l}(0)\bigr\Vert _{2}^{2}+ \frac{1}{2}\bigl\Vert \nabla{u}_{lt}(0)\bigr\Vert _{2}^{2}-\frac{b}{p}\bigl\Vert u_{l}(0) \bigr\Vert _{p}^{p}, \end{aligned}$$

since

$$\begin{aligned} u_{l}(t)=u_{l}(0)+\int_{0}^{t}u_{lt} \,d\tau. \end{aligned}$$

By the Minkowski inequality, we have

$$ \begin{aligned} &\bigl\Vert u_{l}(t)\bigr\Vert _{p}\leq\bigl\Vert u_{l}(0)\bigr\Vert _{p}+\int_{0}^{t} \Vert u_{lt}\Vert _{p}\,d\tau, \\ & \bigl\Vert u_{l}(t)\bigr\Vert _{p}^{p} \leq2^{p} \biggl[\bigl\Vert u_{l}(0)\bigr\Vert _{p}^{p}+\biggl(\int_{0}^{t} \Vert u_{lt}\Vert _{p}\,d\tau\biggr)^{p} \biggr]. \end{aligned} $$
(3.2)

Taking \(q=\frac{p}{p-1}\) and considering the Hölder inequality, we get

$$\begin{aligned} \biggl(\int_{0}^{t}\Vert u_{lt}\Vert _{p}\,d\tau\biggr)^{p}\leq\biggl(\int _{0}^{t}1^{q}\,d\tau \biggr)^{\frac{p}{q}}\int_{0}^{t}\Vert u_{lt}\Vert _{p}^{p}\,d\tau=t^{p-1}\int _{0}^{t}\Vert u_{lt}\Vert _{p}^{p}\,d\tau. \end{aligned}$$
(3.3)

Using the Young inequality, we see that

$$\begin{aligned} t^{p-1}\vert u_{lt}\vert ^{p}\leq \varepsilon\frac {p}{m}\vert u_{lt}\vert ^{m}+c( \varepsilon)\frac{m-p}{m}t^{\beta},\quad m\geq{p}, \end{aligned}$$
(3.4)

where \(\beta=(p-1)\frac{m}{m-p}\), \(c(\varepsilon)=\varepsilon^{\frac{p}{p-m}}\). Hence, we have

$$\begin{aligned} t^{p-1}\int_{0}^{t}\Vert u_{lt}\Vert _{p}^{p}\,d\tau\leq\varepsilon \frac {p}{m}\Vert u_{lt}\Vert ^{m}+c(\varepsilon) \frac{m-p}{m}t^{\beta+1}\vert \Omega \vert . \end{aligned}$$
(3.5)

From (3.2)-(3.5), we get

$$\begin{aligned} \frac{b}{p}\bigl\Vert u_{l}(t)\bigr\Vert _{p}^{p}&\leq\frac{b}{p}2^{p} \biggl[\bigl\Vert u_{l}(0)\bigr\Vert _{p}^{p}+\biggl(\int _{0}^{t}\Vert u_{lt}\Vert _{p}\,d\tau\biggr)^{p} \biggr] \\ & =\frac{b}{p}2^{p}\bigl\Vert u_{l}(0)\bigr\Vert _{p}^{p}+\frac{b}{p}2^{p}\biggl(\int _{0}^{t}\Vert u_{lt}\Vert _{p}\,d\tau\biggr)^{p} \\ & \leq\frac{b}{p}2^{p}\bigl\Vert u_{l}(0)\bigr\Vert _{p}^{p}+\varepsilon\frac {b}{m}2^{p} \Vert u_{lt}\Vert ^{m}+c(\varepsilon) \frac{m-p}{pm}b2^{p}t^{\beta +1}\vert \Omega \vert . \end{aligned}$$
(3.6)

Choosing ε such that \(\varepsilon\frac{b}{m}2^{p}=\frac {a}{2}\), we see that

$$\begin{aligned} &\frac{1}{2}\Vert u_{lt}\Vert _{H^{1}}^{2}+ \frac{1}{2}\Vert \nabla{u}_{l}\Vert _{2}^{2}+ \frac{1}{2}\Vert \triangle{u}_{l}\Vert _{2}^{2}+ \int_{0}^{t}\Vert \nabla{u}_{lt} \Vert _{2}^{2}\,d\tau+\frac{a}{2}\int _{0}^{t}\Vert u_{lt}\Vert _{m}^{m}\,d\tau \\ &\quad \leq\frac{1}{2}\bigl\Vert u_{lt}(0)\bigr\Vert _{H^{1}}^{2}+\frac{1}{2}\bigl\Vert \nabla {u}_{l}(0)\bigr\Vert _{2}^{2}+\frac{1}{2} \bigl\Vert \triangle{u}_{l}(0)\bigr\Vert _{2}^{2} \\ & \qquad {} -\frac{(2^{p}-1)b}{p}\bigl\Vert u_{l}(0)\bigr\Vert _{p}^{p}+c(\varepsilon)\frac {m-p}{pm}b2^{p}t^{\beta+1} \vert \Omega \vert \\ & \quad \leq M(T). \end{aligned}$$
(3.7)

Thus, we find that

$$\Vert u_{lt}\Vert _{H^{1}}^{2},\qquad \Vert \nabla{u}_{l}\Vert _{2}^{2}, \qquad \Vert \triangle {u}_{l}\Vert _{2}^{2},\qquad \int_{0}^{t} \Vert \nabla{u}_{lt}\Vert _{2}^{2}\,d\tau, \qquad \int _{0}^{t}\Vert u_{lt}\Vert _{m}^{m}\,d\tau $$

are finite and can be controlled by a constant depending on T. Together with (3.6), we find that

$$\begin{aligned} \Vert u_{l}\Vert _{p}^{p}\leq{M}(T),\quad 0\leq{t} \leq{T}. \end{aligned}$$

Hence, \(u_{l}\in{L}^{\infty}(0,T;H_{0}^{2}(\Omega)\cap{L}^{p}(\Omega))\), \(u_{lt}\in{L}^{\infty}(0,T;H_{0}^{1}(\Omega))\cap{L}^{m}(Q_{T})\). We see that

$$\begin{aligned} \bigl(a\vert u_{lt}\vert ^{m-2}u_{lt},u_{lt} \bigr)=a\Vert u_{lt}\Vert _{m}^{m},\qquad \bigl(b \vert {u}_{l}\vert ^{p-2}{u}_{l},u_{l} \bigr)=b\Vert u\Vert _{p}^{p}, \end{aligned}$$

so we obtain \(a\vert u_{lt}\vert ^{m-2}u_{lt}\in{L^{m'}}(Q_{T})\), \(b\vert {u}_{l}\vert ^{p-2}{u}_{l}\in{L}^{\infty}(0,T;L^{p'}(\Omega))\).

Thus, we can show that there exist a subsequence \(\{u_{v}\}\) from the sequences \(\{u_{l}\}\) and limit functions u, ξ, η, such that

$$\begin{aligned}& u_{v}\longrightarrow{u} \quad \mbox{in } {L}^{\infty}\bigl(0,T;{H}_{0}^{2}(\Omega)\cap {L}^{p}(\Omega)\bigr)\mbox{ weakly star, }v\longrightarrow\infty, \\& u_{v}\longrightarrow{u}\quad \mbox{a.e. in }Q=\Omega \times[0,T), v\longrightarrow\infty, \\& u_{vt}\longrightarrow{u_{t}} \quad \mbox{in } {L}^{\infty}\bigl(0,T;H_{0}^{1}(\Omega)\bigr)\mbox{ weakly star, }v\longrightarrow\infty, \\& u_{vt}\longrightarrow{u_{t}}\quad \mbox{in } {L}^{m}(Q_{T})\mbox{ weakly, }v\longrightarrow\infty, \\& a\vert u_{vt}\vert ^{m-2}u_{vt} \longrightarrow\xi\quad \mbox{in }{L^{m'}}(Q_{T})\mbox{ weakly, }v \longrightarrow\infty, \\& b\vert {u}_{v}\vert ^{p-2}{u}_{v} \longrightarrow\eta \quad \mbox{in }{L}^{\infty }\bigl(0,T;{L}^{p'}(\Omega) \bigr)\mbox{ weakly star, }v\longrightarrow\infty. \end{aligned}$$

Making use of the Lions lemma in [27], it follows that \(b\vert {u}_{l}\vert ^{p-2}{u}_{l}\longrightarrow{b} \vert {u}\vert ^{p-2}{u}=\eta\).

Integrating (2.8) with respect to t, we get

$$\begin{aligned} &(u_{lt},w_{j})+(\nabla{u}_{lt}, \nabla{w}_{j})+\int_{0}^{t}(\nabla {u}_{l},\nabla{w}_{j})\,d\tau \\ & \qquad {} +(\nabla{u}_{l},\nabla{w}_{j})+\int _{0}^{t}(\triangle {u_{l}}, \triangle{w}_{j})\,d\tau+a\int_{0}^{t} \bigl(\vert u_{lt}\vert ^{m-2}u_{lt},w_{j} \bigr)\,d\tau \\ & \quad =b\int_{0}^{t}\bigl(\vert {u}_{l}\vert ^{p-2}{u}_{l},w_{j}\bigr) \,d\tau +\bigl(u_{lt}(0),w_{j}\bigr)+\bigl( \nabla{u}_{lt}(0),\nabla{w}_{j}\bigr)+\bigl(\nabla {u}_{l}(0),\nabla{w}_{j}\bigr). \end{aligned}$$
(3.8)

Taking \(l=v\rightarrow\infty\) in (3.8), we have

$$\begin{aligned} &(u_{t},w_{j})+(\nabla{u}_{t}, \nabla{w}_{j})+\int_{0}^{t}(\nabla {u}, \nabla{w}_{j})\,d\tau \\ &\qquad {} +(\nabla{u},\nabla{w}_{j})+\int_{0}^{t}( \triangle {u},\triangle{w}_{j})\,d\tau+\int_{0}^{t}( \xi,w_{j})\,d\tau \\ & \quad =b\int_{0}^{t}\bigl(\vert {u}\vert ^{p-2}{u},w_{j}\bigr)\,d\tau+(u_{1},w_{j})+( \nabla {u}_{1},\nabla{w}_{j})+(\nabla{u}_{0}, \nabla{w}_{j}), \end{aligned}$$
(3.9)

and consequently, differentiating (3.9) with respect to t, we deduce

$$\begin{aligned} &(u_{tt},w_{j})+(\nabla{u},\nabla{w}_{j})+( \nabla{u}_{t},\nabla {w}_{j})+(\triangle{u}, \triangle{w}_{j}) \\ &\quad {}+(\nabla{u}_{tt},\nabla{w}_{j})+(\xi,w_{j})= \bigl(b\vert {u}\vert ^{p-2}{u},w_{j}\bigr). \end{aligned}$$
(3.10)

Considering that the basis \(\{w_{j}(x)\}_{j=1}^{\infty}\) is dense in \(H_{0}^{2}(\Omega)\cap{L}^{m}(\Omega)\), we choose a function \(v\in {H_{0}^{2}(\Omega)\cap{L}^{m}(\Omega)}\) having the form \(v=\sum_{j=1}^{\infty}d_{j}w_{j}(x)\), where \(\{d_{j}\}_{1}^{\infty}\) are given functions. Multiplying (3.10) by \(d_{j}\), summing for \(j=1,\ldots \) , then we deduce

$$\begin{aligned} &(u_{tt},v)+(\nabla{u},\nabla{v})+(\nabla{u}_{t},\nabla {v})+(\triangle{u},\triangle{v})+(\nabla{u}_{tt},\nabla{v}) +(\xi,v) \\ &\quad =\bigl(b\vert {u}\vert ^{p-2}{u},v\bigr),\quad \forall v \in{H}_{0}^{2}(\Omega)\cap {L}^{m}(\Omega). \end{aligned}$$

Next, we need to prove that

$$\begin{aligned} \xi=a\vert u_{t}\vert ^{m-2}u_{t}. \end{aligned}$$
(3.11)

In fact, from (2.8) we obtain

$$\begin{aligned} & \frac{1}{2}\Vert u_{lt}\Vert _{H^{1}}^{2}+\frac{1}{2}\Vert \nabla{u}_{l} \Vert _{2}^{2}+\frac{1}{2}\Vert \triangle{u}_{l}\Vert _{2}^{2}+\int _{0}^{t}\Vert \nabla{u}_{lt}\Vert _{2}^{2}\,d\tau+\int_{0}^{t} \bigl(\gamma (u_{lt}),u_{lt}\bigr)\,d\tau \\ & \quad =\frac{1}{2}\bigl\Vert u_{lt}(0)\bigr\Vert _{H^{1}}^{2}+\frac{1}{2}\bigl\Vert \nabla {u}_{l}(0)\bigr\Vert _{2}^{2}+\frac{1}{2} \bigl\Vert \triangle{u}_{l}(0)\bigr\Vert _{2}^{2}+ \int_{0}^{t}\bigl(b\vert {u}_{l} \vert ^{p-2}{u}_{l},u_{lt}\bigr)\,d\tau, \end{aligned}$$

where \(\gamma(u_{lt})=a\vert u_{lt}\vert ^{m-2}u_{lt}\). From what has been discussed above, taking \(l\longrightarrow\infty\), we have

$$\begin{aligned} &\frac{1}{2}\Vert u_{t}\Vert _{H^{1}}^{2}+\frac{1}{2}\Vert \nabla{u} \Vert _{2}^{2}+\frac{1}{2}\Vert \triangle{u} \Vert _{2}^{2}+\int_{0}^{t}\Vert \nabla {u}_{t}\Vert _{2}^{2}\,d\tau+\liminf\int _{0}^{t}\bigl(\gamma (u_{lt}),u_{lt} \bigr)\,d\tau \\ & \quad \leq\frac{1}{2}\Vert u_{1}\Vert _{H^{1}}^{2}+ \frac{1}{2}\Vert \nabla{u}_{0}\Vert _{2}^{2}+ \frac{1}{2}\Vert \triangle{u}_{0}\Vert _{2}^{2}+ \int_{0}^{t}\bigl(b\vert {u}\vert ^{p-2}{u},u_{t}\bigr)\,d\tau. \end{aligned}$$
(3.12)

Combining (3.12) and Lemma 2.2, we deduce

$$\begin{aligned} \int_{0}^{t}(g,u_{t})\,d\tau\leq\int _{0}^{t}\bigl(b\vert {u}\vert ^{p-2}{u},u_{t}\bigr)\,d\tau-\liminf\int_{0}^{t} \bigl(\gamma (u_{lt}),u_{lt}\bigr)\,d\tau. \end{aligned}$$
(3.13)

Noting that \(m'=\frac{m}{m-1}<\frac{p}{p-1}=p'\), it follows that \(b\vert {u}\vert ^{p-2}{u}\in{L}^{\infty}(0,T;{L}^{p'}(\Omega))\subset {L^{m'}}(Q_{T})\). From Lemma 2.2, we know that the function \(g\in {L}^{2}(0,T;H^{-1}(\Omega))+{L}^{m'}(Q_{T})\). Hence, we can choose \(g=b\vert {u}\vert ^{p-2}{u}-\xi\) such that

$$\begin{aligned} \liminf\int_{0}^{t}\bigl(\gamma(u_{lt}),u_{lt} \bigr)\,d\tau\leq\int_{0}^{t}( \xi,u_{t})\,d\tau. \end{aligned}$$
(3.14)

Considering \(u_{vt}\rightarrow{u_{t}}\) weakly in \({L}^{m}(Q_{T})\) and \(a\vert u_{vt}\vert ^{m-2}u_{vt}\rightarrow\xi\) weakly in \({L^{m'}}(Q_{T})\), \(\forall\varphi\in{L}^{m}(Q_{T})\), we have

$$\begin{aligned}& \lim\int_{0}^{t}\bigl( \gamma(u_{lt}),-\varphi\bigr)\,d\tau=\int_{0}^{t}( \xi,-\varphi)\,d\tau, \end{aligned}$$
(3.15)
$$\begin{aligned}& \lim\int_{0}^{t}\bigl(- \gamma(\varphi),u_{lt}-\varphi\bigr)\,d\tau=\int_{0}^{t} \bigl(-\gamma(\varphi),u_{t}-\varphi\bigr)\,d\tau. \end{aligned}$$
(3.16)

By the combination of (3.14), (3.15), and (3.16), it follows that

$$\begin{aligned} \liminf\int_{0}^{t}\bigl(\gamma(u_{lt})- \gamma(\varphi),u_{lt}-\varphi \bigr)\,d\tau\leq\int _{0}^{t}\bigl(\xi-\gamma(\varphi),u_{t}- \varphi\bigr)\,d\tau. \end{aligned}$$
(3.17)

Utilizing the monotonicity of the function \(\gamma(s)=a\vert s\vert ^{m-2}s\), it means that

$$\begin{aligned} \bigl(\gamma(u_{lt})-\gamma(\varphi),u_{lt}-\varphi \bigr)=a(m-1)\int_{\Omega }\vert u_{lt}+\theta\varphi \vert ^{m-2}\vert u_{lt}-\varphi \vert ^{2}\,dx \geq0, \end{aligned}$$
(3.18)

where \(0<\theta<1\). Thus, from (3.17) and (3.18), we have

$$\begin{aligned} \int_{0}^{t}\bigl(\xi-\gamma( \varphi),u_{t}-\varphi\bigr)\,d\tau\geq0. \end{aligned}$$

Taking \(t\longrightarrow{T}\), we obtain

$$\begin{aligned} \int_{0}^{T}\bigl(\xi-\gamma( \varphi),u_{t}-\varphi\bigr)\,d\tau\geq0. \end{aligned}$$
(3.19)

In order to prove (3.11) from (3.19), we use the semi-continuity of the function \(\gamma(s)\) (\(s\in R\)). Let \(\varphi=u_{t}-\lambda{w}_{t}\), \(\lambda>0\), and \(\forall w_{t}\in{L}^{\infty}(0,T;H_{0}^{1}(\Omega))\cap{L}^{m}(Q_{T})\), then

$$\begin{aligned} \lambda\int_{0}^{T}\bigl(\xi- \gamma(u_{t}-\lambda{w}_{t}),w_{t}\bigr)\,d\tau \geq0, \end{aligned}$$

and

$$\begin{aligned} \int_{0}^{T}\bigl(\xi-\gamma(u_{t}- \lambda{w}_{t}),w_{t}\bigr)\,d\tau\geq0. \end{aligned}$$

Passing to the limits as \(\lambda\longrightarrow0\), we obtain

$$\begin{aligned} \int_{0}^{T}\bigl(\xi-\gamma(u_{t}),w_{t} \bigr)\,d\tau\geq0,\quad \forall w_{t}\in {L}^{\infty} \bigl(0,T;H_{0}^{1}(\Omega)\bigr)\cap{L}^{m}(Q_{T}). \end{aligned}$$
(3.20)

In a similar way, let \(\varphi=u_{t}-\lambda{w}_{t}\), \(\lambda<0\), and \(\forall w_{t}\in{L}^{\infty}(0,T;H_{0}^{1}(\Omega))\cap{L}^{m}(Q_{T})\), then

$$\begin{aligned} \int_{0}^{T}\bigl(\xi-\gamma(u_{t}),w_{t} \bigr)\,d\tau\leq0, \quad \forall w_{t}\in {L}^{\infty} \bigl(0,T;H_{0}^{1}(\Omega)\bigr)\cap{L}^{m}(Q_{T}). \end{aligned}$$
(3.21)

From (3.20), (3.21), we see that \(\xi=a\vert u_{t}\vert ^{m-2}u_{t}\). Thus, the theorem is completed. □

4 The asymptotic behavior of global weak solutions

In this section, we consider the asymptotic behavior of global weak solutions for the problem (1.8).

To obtain the results of this section, we now define some functionals as follows:

$$\begin{aligned}& I\bigl(u(t)\bigr)=\Vert \nabla{u} \Vert _{2}^{2}+ \Vert \triangle{u} \Vert _{2}^{2}-{b}\Vert {u}\Vert _{p}^{p}, \\& J\bigl(u(t)\bigr)=\frac{1}{2}\Vert \nabla{u} \Vert _{2}^{2}+\frac{1}{2}\Vert \triangle{u} \Vert _{2}^{2}-\frac{b}{p}\Vert {u}\Vert _{p}^{p}, \end{aligned}$$

and

$$\begin{aligned} E(t) =&E\bigl(u(t)\bigr)=\frac{1}{2}\Vert {u}_{t}\Vert _{2}^{2}+\frac{1}{2}\Vert \nabla {u}_{t}\Vert _{2}^{2}+J\bigl(u(t)\bigr) \\ =&\frac{1}{2}\Vert {u}_{t}\Vert _{2}^{2}+ \frac{1}{2}\Vert \nabla{u}_{t}\Vert _{2}^{2}+ \frac{1}{2}\Vert \nabla{u} \Vert _{2}^{2}+ \frac{1}{2}\Vert \triangle {u}\Vert _{2}^{2}- \frac{b}{p}\Vert {u}\Vert _{p}^{p}. \end{aligned}$$

Next, let us introduce the set

$$\begin{aligned} W= \bigl\{ u\in{H}_{0}^{2}(\Omega)|I\bigl(u(t)\bigr)>0 \bigr\} \cup\{0\}. \end{aligned}$$

Lemma 4.1

Let \((u_{0},u_{1})\in{W}\times{H}_{0}^{1}(\Omega)\) be given. Assume that \(a,b>0\) and

$$\begin{aligned}& 2< p\leq{m}<\infty \quad \textit{if } n=1,2 \quad \textit{and}\quad 2<p\leq{m}\leq\frac{2n}{n-2}\quad \textit{if } n\geq3, \end{aligned}$$
(4.1)
$$\begin{aligned}& \alpha=bC_{*p}^{p} \biggl[ \frac{2p}{p-2}E(0) \biggr]^{\frac{p-2}{2}}< 1, \end{aligned}$$
(4.2)

where \(C_{*p}\) is the Sobolev constant of the space \(H_{0}^{1}(\Omega )\) into \(L^{p}(\Omega)\). Then for any \(t\geq0\), the global weak solutions of the problem (1.8) satisfy \(u(t)\in{W}\).

Proof

Since \(I(u_{0})>0\) and the time continuity of \(I(u)\), then there exists \(t_{1}\) such that \(I(u(t))\geq0\) for all \(t\in [0,t_{1}]\). Thus, we see that

$$ \begin{aligned}[b] J\bigl(u(t)\bigr)&=\frac{1}{2}\Vert \nabla{u} \Vert _{2}^{2}+\frac{1}{2}\Vert \triangle {u}\Vert _{2}^{2}-\frac{b}{p}\Vert {u}\Vert _{p}^{p} \\ & =\frac{1}{2}\bigl[\Vert \nabla{u} \Vert _{2}^{2}+ \Vert \triangle{u} \Vert _{2}^{2}\bigr]-\frac{1}{p} \bigl[\Vert \nabla{u} \Vert _{2}^{2}+\Vert \triangle{u} \Vert _{2}^{2}\bigr]+\frac{1}{p}I\bigl(u(t)\bigr) \\ & \geq\frac{p-2}{2p}\bigl[\Vert \nabla{u} \Vert _{2}^{2}+ \Vert \triangle{u} \Vert _{2}^{2}\bigr], \end{aligned} $$
(4.3)

for any \(t\in[0,t_{1}]\). Hence, we obtain

$$\begin{aligned} \Vert \nabla{u} \Vert _{2}^{2}+\Vert \triangle{u} \Vert _{2}^{2}\leq\frac {2p}{p-2}J(t)\leq \frac{2p}{p-2}E(t)\leq\frac{2p}{p-2}E(0), \end{aligned}$$
(4.4)

for any \(t\in[0,t_{1}]\). If \(\Vert \nabla{u} \Vert _{2}^{2}=0\), we have \(u=0 \) (since \(u\in{H}_{0}^{2}(\Omega)\)). Hence, we have from the definition of the set W: \(u\in{W}\). If \(\Vert \nabla{u} \Vert _{2}^{2}\neq0\), by the Sobolev inequality and (4.2), (4.4), we have

$$\begin{aligned} b\Vert {u}\Vert _{p}^{p} \leq&{b}C_{*p}^{p}\Vert \nabla{u} \Vert _{2}^{p}={b}C_{*p}^{p} \Vert \nabla{u} \Vert _{2}^{p-2}\Vert \nabla{u} \Vert _{2}^{2} \\ \leq&{b}C_{*p}^{p} \biggl[\frac{2p}{p-2}E(0) \biggr]^{\frac {p-2}{2}}\Vert \nabla{u} \Vert _{2}^{2}=\alpha \Vert \nabla{u} \Vert _{2}^{2}< \Vert \nabla {u}\Vert _{2}^{2}, \end{aligned}$$
(4.5)

for any \(t\in[0,t_{1}]\). Hence, we obtain \(\Vert \nabla{u} \Vert _{2}^{2}-{b}\Vert {u}\Vert _{p}^{p}>0\), \(\forall t\in[0,t_{1}]\). This shows that \(u(t)\in {W}\), \(\forall t\in[0,t_{1}]\). We see that

$$\begin{aligned} \lim_{t\longrightarrow{t_{1}}}{b}C_{*p}^{p} \biggl[\frac {2p}{p-2}E(t) \biggr]^{\frac{p-2}{2}}\leq\alpha< 1, \end{aligned}$$
(4.6)

so the above argument may be repeated, and the solution can thus be extended to the time \(t_{1}\leq{t}<{t}_{2}\). Continuing in this way, the assertion of the lemma is proved. □

Lemma 4.2

Assume that (4.1) and (4.2) hold. Then the global weak solutions of the problem (1.8) satisfy

$$\begin{aligned} \bigl\Vert u(t)\bigr\Vert _{m}^{m} \leq{L}E(t), \end{aligned}$$
(4.7)

for some constant L depending on m, p and \(E(0)\) only.

Proof

Using the Sobolev inequality and (4.4), we have

$$\begin{aligned} \bigl\Vert u(t)\bigr\Vert _{m}^{m} \leq&{C}_{*m}^{m}\Vert \nabla{u} \Vert _{2}^{m}={C}_{*m}^{m} \Vert \nabla{u} \Vert _{2}^{m-2}\Vert \nabla{u} \Vert _{2}^{2} \\ \leq&{C}_{*m}^{m} \biggl[\frac{2p}{p-2}E(0) \biggr]^{\frac {m-2}{2}}\frac{2p}{p-2}E(t)\leq{L}E(t), \end{aligned}$$
(4.8)

where \(L={C}_{*m}^{m} [\frac{2p}{p-2}E(0) ]^{\frac {m-2}{2}}\frac{2p}{p-2}\), and \(C_{*m}\) is the Sobolev constant of the space \(H_{0}^{1}(\Omega)\) into \(L^{m}(\Omega)\), and the proof of the lemma is completed. □

Theorem 4.1

Assume that \(a,b>0\), and the conditions (4.1), (4.2) hold. Let \((u_{0},u_{1})\in{W}\times{H}_{0}^{1}(\Omega)\) be given. Then for the global weak solutions of the problem (1.8), there exist positive constants M and k such that

$$\begin{aligned} E(t)\leq{M}e^{-kt},\quad \forall t\geq0. \end{aligned}$$
(4.9)

Proof

From Lemma 4.1, we know that, for any \(t\geq0\), the global weak solutions of the problem (1.8) satisfy \(u(t)\in{W}\). Now defining

$$\begin{aligned} F(t)=E(t)+\varepsilon\int_{\Omega}uu_{t} \,dx+\varepsilon\int_{\Omega }\nabla{u}\nabla{u}_{t} \,dx+\frac{\varepsilon}{2}\int_{\Omega }\vert \nabla{u} \vert ^{2}\,dx, \end{aligned}$$
(4.10)

we can show that for ε small enough, there exist two positive constants \(C_{1}\) and \(C_{2}\) such that

$$\begin{aligned} C_{1}E(t)\leq{F}(t)\leq{C}_{2}E(t). \end{aligned}$$
(4.11)

In fact,

$$\begin{aligned} F(t)&\leq{E}(t)+\frac{\varepsilon}{2}\Vert u_{t}\Vert _{2}^{2}+\frac {\varepsilon}{2}\Vert u\Vert _{2}^{2}+\frac{\varepsilon}{2}\Vert \nabla{u}_{t} \Vert _{2}^{2}+\varepsilon \Vert \nabla{u} \Vert _{2}^{2} \\ & \leq(1+\varepsilon)E(t)+\varepsilon{C}(\Omega)\Vert \nabla{u} \Vert _{2}^{2} \\ & \leq(1+\varepsilon)E(t)+\varepsilon\frac{2p}{p-2}E(t) \\ & \leq {C}_{2}E(t) \end{aligned}$$
(4.12)

and

$$\begin{aligned} F(t)&\geq{E}(t)-\frac{\varepsilon}{4\gamma} \Vert u_{t}\Vert _{2}^{2}-\varepsilon\gamma \Vert u\Vert _{2}^{2}+\frac{\varepsilon}{2}\Vert \nabla{u} \Vert _{2}^{2}-\frac{\varepsilon}{4\gamma} \Vert \nabla{u}_{t} \Vert _{2}^{2}-\varepsilon\gamma \Vert \nabla{u} \Vert _{2}^{2} \\ & \geq E(t)-\frac{\varepsilon}{4\gamma} \bigl[\Vert u_{t}\Vert _{2}^{2}+\Vert \nabla{u}_{t}\Vert _{2}^{2} \bigr]+\varepsilon\biggl(\frac{1}{2}-\gamma - \gamma{C}_{*2}\biggr)\Vert \nabla{u} \Vert _{2}^{2}, \end{aligned}$$
(4.13)

where \(C_{*2}\) is the Sobolev constant of the space \(H_{0}^{1}(\Omega )\) into \(L^{2}(\Omega)\).

By choosing γ small enough, we have

$$\begin{aligned} F(t) \geq&{E}(t)-\frac{\varepsilon}{4\gamma} \bigl[\Vert u_{t} \Vert _{2}^{2}+\Vert \nabla{u}_{t}\Vert _{2}^{2} \bigr] \\ \geq&{J}\bigl(u(t)\bigr)+\biggl(\frac{1}{2}-\frac{\varepsilon}{4\gamma}\biggr) \Vert u_{t}\Vert _{2}^{2}+\biggl( \frac{1}{2}-\frac{\varepsilon}{4\gamma}\biggr)\Vert \nabla{u}_{t}\Vert _{2}^{2}. \end{aligned}$$
(4.14)

Once γ is chosen, we take ε so small that

$$\begin{aligned} F(t)\geq{J}\bigl(u(t)\bigr)+\frac{C_{1}}{2}\Vert u_{t}\Vert _{2}^{2}+\frac{C_{1}}{2}\Vert \nabla{u}_{t}\Vert _{2}^{2}\geq{C}_{1}E(t), \end{aligned}$$
(4.15)

where \(\frac{C_{1}}{2}\leq\frac{1}{2}-\frac{\varepsilon}{4\gamma}\). Now differentiating (4.10) and utilizing (1.8), Lemma 2.1, and the Poincaré inequality, we have

$$\begin{aligned} F'(t) =&{E'}(t)+\varepsilon\int _{\Omega}uu_{tt}\,dx+\varepsilon \Vert u_{t} \Vert _{2}^{2}+\varepsilon\int_{\Omega} \nabla{u}\nabla {u}_{tt}\,dx+\varepsilon \Vert \nabla{u}_{t} \Vert _{2}^{2}+\frac{\varepsilon }{2}\frac{d}{dt}\Vert \nabla{u} \Vert _{2}^{2} \\ =& -\bigl(a\Vert u_{t}\Vert _{m}^{m}+ \Vert \nabla{u}_{t}\Vert _{2}^{2}\bigr)+ \varepsilon\bigl[\Vert u_{t}\Vert _{2}^{2}+ \Vert \nabla{u}_{t}\Vert _{2}^{2}-\Vert \nabla{u} \Vert _{2}^{2}-\Vert \triangle{u} \Vert _{2}^{2}\bigr] \\ &{} -a\varepsilon\int_{\Omega} \vert u_{t}\vert ^{m-2}u_{t}{u}\,dx+\varepsilon {b}\Vert u\Vert _{p}^{p} \\ \leq&-a\Vert u_{t}\Vert _{m}^{m}-[1-C_{*2} \varepsilon-\varepsilon]\Vert \nabla {u}_{t}\Vert _{2}^{2}-\varepsilon \Vert \nabla{u} \Vert _{2}^{2}-\varepsilon \Vert \triangle{u} \Vert _{2}^{2} \\ &{} -a\varepsilon\int_{\Omega} \vert u_{t}\vert ^{m-2}u_{t}{u}\,dx+\varepsilon {b}\Vert u\Vert _{p}^{p}. \end{aligned}$$
(4.16)

Using the energy functionals and the Sobolev inequality, we have

$$ \begin{aligned}[b] b\Vert u\Vert _{p}^{p}={}&\lambda{b} \Vert u\Vert _{p}^{p}+(1-\lambda)b\Vert u\Vert _{p}^{p} \\ \leq{}&\lambda \biggl[\frac{p}{2}\Vert {u}_{t}\Vert _{2}^{2}+\frac{p}{2}\Vert \nabla{u}_{t} \Vert _{2}^{2}+\frac{p}{2}\Vert \nabla{u} \Vert _{2}^{2}+\frac {p}{2}\Vert \triangle{u} \Vert _{2}^{2} \biggr] \\ & {} -p\lambda{E}(t)+(1-\lambda)\alpha \Vert \nabla{u} \Vert _{2}^{2}, \end{aligned} $$
(4.17)

where \(\alpha=bC_{*p}^{p} [\frac{2p}{p-2}E(0) ]^{\frac {p-2}{2}}<1\).

Utilizing Lemma 4.2 and inserting (4.17) into (4.16), we have

$$\begin{aligned} F'(t) \leq& -a\Vert u_{t}\Vert _{m}^{m}-[1-C_{*2}\varepsilon-\varepsilon]\Vert \nabla{u}_{t}\Vert _{2}^{2}-\varepsilon \Vert \nabla{u} \Vert _{2}^{2}-\varepsilon \Vert \triangle{u} \Vert _{2}^{2} \\ &{} +\varepsilon\lambda \biggl[\frac{p}{2}\Vert {u}_{t} \Vert _{2}^{2}+\frac {p}{2}\Vert \nabla{u}_{t}\Vert _{2}^{2}+\frac{p}{2} \Vert \nabla{u} \Vert _{2}^{2}+\frac{p}{2}\Vert \triangle{u} \Vert _{2}^{2} \biggr] \\ &{} -p\varepsilon\lambda{E}(t)+(1-\lambda)\alpha \Vert \nabla{u} \Vert _{2}^{2}\,dx+a\varepsilon\delta\int_{\Omega} \vert u\vert ^{m}\,dx+C(\delta )a\varepsilon\int _{\Omega} \vert u_{t}\vert ^{m}\,dx \\ \leq&-a\Vert u_{t}\Vert _{m}^{m}-[1-C_{*2} \varepsilon-\varepsilon]\Vert \nabla {u}_{t}\Vert _{2}^{2}-\varepsilon \Vert \nabla{u} \Vert _{2}^{2}-\varepsilon \Vert \triangle{u} \Vert _{2}^{2} \\ &{} +\varepsilon\lambda \biggl[\frac{p}{2}\Vert {u}_{t} \Vert _{2}^{2}+\frac {p}{2}\Vert \nabla{u}_{t}\Vert _{2}^{2}+\frac{p}{2} \Vert \nabla{u} \Vert _{2}^{2}+\frac{p}{2}\Vert \triangle{u} \Vert _{2}^{2} \biggr]-p\varepsilon \lambda{E}(t) \\ &{} +(1-\lambda)\alpha \Vert \nabla{u} \Vert _{2}^{2} \,dx+a\varepsilon\delta {L}E(t)+a\varepsilon{C}(\delta) \bigl[\Vert u_{t}\Vert _{m}^{m}+\Vert \nabla {u}_{t}\Vert _{2}^{2} \bigr] \\ \leq&-a \bigl[1-\varepsilon{C}(\delta) \bigr]\Vert u_{t}\Vert _{m}^{m}-\varepsilon [\lambda{p}-a\delta{L} ]E(t) \\ &{} +\varepsilon \biggl[(1-\lambda)\alpha+\frac{p}{2}\lambda -1 \biggr] \Vert \nabla{u} \Vert _{2}^{2}+\varepsilon \biggl[ \frac{p}{2}\lambda -1 \biggr]\Vert \triangle{u} \Vert _{2}^{2} \\ &{} - \biggl\{ 1-\varepsilon \biggl[C_{*2}+1+C(\delta)a+\lambda \frac {p}{2}C_{*2}+\lambda\frac{p}{2} \biggr] \biggr\} \Vert \nabla{u}_{t}\Vert _{2}^{2}, \end{aligned}$$
(4.18)

where L is the constant of Lemma 4.2, δ is any positive constant, and \(C(\delta)\) is a constant depending on δ, m only.

Thus, we see that

$$\begin{aligned} F'(t) \leq&-a \bigl[1-\varepsilon{C}(\delta) \bigr] \Vert u_{t}\Vert _{m}^{m}-\varepsilon [ \lambda{p}-a\delta{L} ]E(t) \\ &{} +\varepsilon \biggl[(1-\lambda)\alpha+\frac{p}{2}\lambda -1 \biggr] \Vert \triangle{u} \Vert _{2}^{2}+\varepsilon \biggl[(1- \lambda )\alpha+\frac{p}{2}\lambda-1 \biggr]\Vert \nabla{u} \Vert _{2}^{2} \\ &{} - \biggl\{ 1-\varepsilon \biggl[C_{*2}+1+C(\delta)a+\lambda \frac {p}{2}C_{*2}+\lambda\frac{p}{2} \biggr] \biggr\} \Vert \nabla{u}_{t}\Vert _{2}^{2} \\ \leq&- \biggl\{ 1-\varepsilon \biggl[C_{*2}+1+C(\delta)a+\lambda \frac{p}{2}C_{*2}+\lambda\frac{p}{2} \biggr] \biggr\} \Vert \nabla {u}_{t}\Vert _{2}^{2}-a \bigl[1- \varepsilon{C}(\delta) \bigr]\Vert u_{t}\Vert _{m}^{m} \\ &{} -\varepsilon [\lambda{p}-a\delta{L} ]E(t)+\varepsilon \biggl[ \frac{p-2}{2}\lambda-(1-\alpha) (1-\lambda) \biggr]\Vert \triangle{u} \Vert _{2}^{2} \\ &{} +\varepsilon \biggl[\frac{p-2}{2}\lambda-(1-\alpha) (1-\lambda ) \biggr]\Vert \nabla{u} \Vert _{2}^{2}. \end{aligned}$$
(4.19)

By choosing λ close to zero such that \(\frac{p-2}{2}\lambda -(1-\alpha)(1-\lambda)\leq0\), then (4.19) becomes

$$\begin{aligned} F'(t) \leq&- \biggl\{ 1-\varepsilon \biggl[C_{*2}+1+C(\delta)a+\lambda \frac{p}{2}C_{*2}+ \lambda\frac{p}{2} \biggr] \biggr\} \Vert \nabla {u}_{t}\Vert _{2}^{2} \\ &{} -\varepsilon(\lambda{p}-a\delta{L})E(t)-a \bigl[1-\varepsilon {C}(\delta) \bigr]\Vert u_{t}\Vert _{m}^{m}. \end{aligned}$$
(4.20)

Once δ is chosen such that \(\lambda{p}-a\delta{L}>0\), we can take ε so small that

$$\begin{aligned} 1-\varepsilon \biggl[C_{*2}+1+C(\delta)a+\lambda\frac {p}{2}C_{*2}+ \lambda\frac{p}{2} \biggr]\geq0,\qquad 1-\varepsilon {C}(\delta)\geq0. \end{aligned}$$

Thus, we see that

$$\begin{aligned} F'(t)\leq-\varepsilon(\lambda{p}-a\delta{L})E(t) \leq-\frac {\varepsilon(\lambda{p}-a\delta{L})}{C_{2}}F(t). \end{aligned}$$
(4.21)

By the Gronwall inequality, we see that

$$\begin{aligned} F(t)\leq{F}(0)e^{-kt}, \quad \forall t\geq0, \end{aligned}$$
(4.22)

where \(k=\frac{\varepsilon(\lambda{p}-a\delta{L})}{C_{2}}\). Combining with (4.11), we obtain

$$\begin{aligned} C_{1}E(t)\leq{F}(t)\leq{F}(0)e^{-kt}, \quad \forall t\geq0, \end{aligned}$$

and

$$\begin{aligned} {E}(t)\leq{M}e^{-kt}, \quad \forall t\geq0, \end{aligned}$$
(4.23)

where \(M=\frac{F(0)}{C_{1}}\). Thus, the proof of the theorem is completed. □

Remark 4.1

From the (2.3), (4.4), (4.5), (4.9), we easily obtain

$$\begin{aligned} \Vert u_{t}\Vert _{2}^{2}+ \Vert \nabla{u} \Vert _{2}^{2}+\Vert \triangle{u} \Vert _{2}^{2}+\Vert \nabla{u}_{t}\Vert _{2}^{2}+\Vert u\Vert _{p}^{p} \leq{C}e^{-{kt}}, \end{aligned}$$
(4.24)

for any \(t\geq0\).

References

  1. Sattinger, DH: On global solution of nonlinear hyperbolic equations. Arch. Ration. Mech. Anal. 30, 148-172 (1968)

    Article  MATH  MathSciNet  Google Scholar 

  2. Payne, LE, Sattinger, DH: Saddle points and instability on nonlinear hyperbolic equations. Isr. J. Math. 22, 273-303 (1975)

    Article  MathSciNet  MATH  Google Scholar 

  3. Liu, YC: On potential wells and vacuum isolating of solutions for semilinear wave equation. J. Differ. Equ. 1992, 155-169 (2003)

    MATH  Google Scholar 

  4. Liu, YC: On potential wells and applications to semilinear hyperbolic equations and parabolic equations. Nonlinear Anal. 64, 2665-2678 (2006)

    Article  MATH  MathSciNet  Google Scholar 

  5. Lasiecka, I, Tataru, D: Uniform boundary stabilization of semilinear wave equations with nonlinear boundary damping. Differ. Integral Equ. 6(3), 507-533 (1993)

    MATH  MathSciNet  Google Scholar 

  6. Vitillaro, E: A potential well method for the wave equation with nonlinear source and boundary damping terms. Glasg. Math. J. 44, 375-395 (2002)

    Article  MATH  MathSciNet  Google Scholar 

  7. Cavalcanti, MM, Cavalcanti, VND, Martinez, P: Existence and decay rate estimates for the wave equation with nonlinear boundary damping and source term. J. Differ. Equ. 203(1), 119-158 (2004)

    Article  MATH  MathSciNet  Google Scholar 

  8. Georgiev, V, Todorova, G: Existence of a solution of the wave equation with nonlinear damping and source terms. J. Differ. Equ. 109, 295-308 (1994)

    Article  MATH  MathSciNet  Google Scholar 

  9. Ikehata, R: Some remarks on the wave equations with nonlinear damping and source terms. Nonlinear Anal. TMA 27, 1165-1175 (1996)

    Article  MATH  MathSciNet  Google Scholar 

  10. Vitillaro, E: Global nonexistence theorems for a class of evolution equations with dissipation. Arch. Ration. Mech. Anal. 149, 155-182 (1999)

    Article  MATH  MathSciNet  Google Scholar 

  11. Messaoudi, SA: Blow up in a nonlinear damped wave equation. Math. Nachr. 231, 105-111 (2001)

    Article  MATH  MathSciNet  Google Scholar 

  12. Messaoudi, SA: Global nonexistence in a nonlinearly damped wave equation. Appl. Anal. 80, 269-277 (2001)

    Article  MATH  MathSciNet  Google Scholar 

  13. Messaoudi, SA: Exponential decay of solutions of a nonlinearly damped wave equation. Nonlinear Differ. Equ. Appl. 12, 391-399 (2005)

    MATH  MathSciNet  Google Scholar 

  14. Serrin, J, Todorova, G, Vitillaro, E: Existence for a nonlinear wave equation with nonlinear damping and source terms. Differ. Integral Equ. 16(1), 13-50 (2003)

    MATH  MathSciNet  Google Scholar 

  15. Hayes, M, Saccomandi, G: Finite amplitude transverse waves in special incompressible viscoelastic solids. J. Elast. 59, 213-225 (2000)

    Article  MATH  Google Scholar 

  16. Shang, YD: Initial-boundary value problem for a class of fourth order nonlinear evolution equations. Math. Appl. 13(1), 7-11 (2000)

    MATH  MathSciNet  Google Scholar 

  17. Villiaggio, P: Mathematical Models for Elastic Structures. Cambridge University Press, Cambridge (1997)

    Book  Google Scholar 

  18. Shang, YD: Initial boundary value problem of equation \(u_{tt}-\triangle {u}-\triangle{u}_{t}-\triangle{u}_{tt}=f(u)\). Acta Math. Appl. Sin. 23, 385-393 (2000)

    MATH  Google Scholar 

  19. Zhang, HW, Hu, QY: Existence of global weak solution and stability of a class nonlinear evolution equation. Acta Math. Sci. 24A, 329-336 (2004)

    MATH  Google Scholar 

  20. Xie, YQ, Zhong, CK: The existence of global attractors for a class nonlinear evolution equation. J. Math. Anal. Appl. 336, 54-69 (2007)

    Article  MATH  MathSciNet  Google Scholar 

  21. Xu, RZ, Zhao, XR, Shen, JH: Asymptotic behavior of solution for fourth order wave equation with dispersive and dissipative terms. Appl. Math. Mech. -Engl. Ed. 29, 259-262 (2008)

    Article  MATH  MathSciNet  Google Scholar 

  22. Xie, YQ, Zhong, CK: Asymptotic behavior of a class of nonlinear evolution equations. Nonlinear Anal. 71(11), 5095-5105 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  23. Carvalho, AN, Cholewa, JW: Local well posedness, asymptotic behavior and asymptotic bootstrapping for a class of semilinear evolution equations of the second order in time. Trans. Am. Math. Soc. 361(5), 2567-2586 (2009)

    Article  MATH  MathSciNet  Google Scholar 

  24. Sun, CY, Yang, L, Duan, JQ: Asymptotic behavior for a semilinear second order evolution equation. Trans. Am. Math. Soc. 363(11), 6085-6109 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  25. Chen, CS, Wang, H, Zhu, SL: Global attractor and decay estimates of solutions to a class of nonlinear evolution equations. Math. Methods Appl. Sci. 34(5), 497-508 (2011)

    Article  MATH  MathSciNet  Google Scholar 

  26. Xu, RZ, Yang, YB: Finite time blow up for the nonlinear fourth-order dispersive-dissipative wave equation at high energy level. Int. J. Math. 23, 1250060 (10 pages) (2012)

    Article  MathSciNet  MATH  Google Scholar 

  27. Lions, JL: Quelques Méthodes de Résolution des Problémes aux Limites Non Linéaires. Dunod Gauthier-Villars, Paris (1969)

    MATH  Google Scholar 

  28. Evance, LC: Partial Differential Equations. Am. Math. Soc., Providence (1988)

    Google Scholar 

  29. Messaoudi, SA: Blow-up of positive-initial-energy solutions of a nonlinear viscoelastic hyperbolic equation. J. Math. Anal. Appl. 320, 902-915 (2006)

    Article  MATH  MathSciNet  Google Scholar 

Download references

Acknowledgements

This work is supported by the NSF of China (11401122, 40890153) and the Scientific Program (2008B080701042) of Guangdong Province.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Yadong Shang.

Additional information

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.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Di, H., Shang, Y. Global existence and asymptotic behavior of solutions for the double dispersive-dissipative wave equation with nonlinear damping and source terms. Bound Value Probl 2015, 29 (2015). https://doi.org/10.1186/s13661-015-0288-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0288-6

MSC

Keywords