# General decay for a wave equation of Kirchhoff type with a boundary control of memory type

Shun-Tang Wu

### Author affiliations

General Education Center, National Taipei University of Technology, Taipei 106, Taiwan

Boundary Value Problems 2011, 2011:55  doi:10.1186/1687-2770-2011-55

 Received: 30 July 2011 Accepted: 23 December 2011 Published: 23 December 2011

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

A nonlinear wave equation of Kirchhoff type with memory condition at the boundary in a bounded domain is considered. We establish a general decay result which includes the usual exponential and polynomial decay rates. Furthermore, our results allow certain relaxation functions which are not necessarily of exponential and polynomial decay. This improves earlier results in the literature.

MSC: 35L05; 35L70; 35L75; 74D10.

##### Keywords:
general decay; wave equation; relaxation; memory type; Kirchhoff type; nondissipative

### 1 Introduction

In this article, we study the asymptotic behavior of the energy function related to a nonlinear wave equation of Kirchhoff type subject to memory condition at the boundary as follows:

(1.1)

(1.2)

(1.3)

(1.4)

where Ω is a bounded domain with smooth boundary ∂Ω = Γ0 ∪ Γ1. The partition Γ0 and Γ1 are closed and disjoint, with meas0) > 0, ν represents the unit normal vector directed towards the exterior of Ω, u is the transverse displacement, and g is the relaxation function considered positive and nonincreasing belonging to W1,2 (Ω).

From the physical point of view, we know that the memory effect described in integral equation (1.3) can be caused by the interaction with another viscoelastic element. In fact, the boundary condition (1.3) signifies that Ω is composed of a material which is clamped in a rigid body in the portion Γ0 of its boundary and is clamped in a body with viscoelastic properties in the portion of Γ1.

When Γ1 = ϕ, problem (1.1) has its origin in describing the nonlinear vibrations of an elastic string. More precisely, we have

(1.5)

for 0 < × < L, t ≥ 0; where u is the lateral deflection, x the space coordinate, t the time, E the Young modulus, ρ the mass density, h the cross section area, L the length, p0 the initial axial tension and f the external force. Kirchhoff [1] was the first one who introduced (1.5) to study the oscillations of stretched strings and plates, so that (1.5) is called the wave equation of Kirchhoff type after him. In this direction, problem (1.1) with ∂Ω = Γ0 and l(t) = 0 has been investigated by many authors in recent years, and many results concerning existence, nonexistence and asymptotic behavior have been established, see [2-13].

On the other hand, regarding the viscoelastic wave equations with memory term acting in the boundary or in the domain, there are numerous results related to asymptotic behavior of solutions. For example, in the case where M(s) = 1, Cavalcanti et al. [14] investigated the existence and uniform decay of strong solutions of wave equation (1.1) with a nonlinear boundary damping of memory type and a nonlinear boundary source when l(t) = 0. Cavalcanti and Guesmia [15] considered the following system:

(1.6)

(1.7)

(1.8)

(1.9)

where Ω is a bounded domain with smooth boundary ∂Ω = Γ0 ∪ Γ1. They obtained the general decay result which depends on the relaxation function g. In particular, if the relaxation function g decays exponentially (or polynomially), then the solution also decays exponentially (or polynomially) and with the same decay rate. Moreover, when u0 = 0 on Γ1, they obtained exponential or polynomial decay of solutions, even if the relaxation function g does not converge to 0 at ∞. Later, Messaoudi and Soufyane [16] generalized this result to the case of a system of Timoshenko type. They established general decay rate results, from which the usual exponential and polynomial decay rates are only special cases. Recently, Messaoudi and Soufyane [17] studied the following problem:

in a bounded domain with boundary conditions (1.7)-(1.9). They improved the results of [15] by applying the multiplier techniques. Indeed, they obtained not only a general decay result, but their works also allowed certain relaxation functions which are not necessarily of exponential or polynomial decay. For other related works, we refer the reader to [18-20] and references therein.

Conversely, in the case where M is not a constant function, Santos [21] considered

where μ(t) is a nonincreasing function satisfying μ(t) ≥ μ0 > 0. By denoting k the resolvent kernel of g', he showed that the solution decays exponentially (or polynomially) to zero provided k decays exponentially (or polynomially) to zero. Later on, Santos et al. [22] generalized this result to a nonlinear n-dimensional equation of Kirchhoff type of the form

(1.10)

in a bounded domain with boundary conditions (1.2)-(1.3). In that article, they proved that the energy decays with the same rate of decay of the relaxation function. This latter result improved an earlier one by Park et al. [23], where the authors considered (1.10) in a bounded domain with nonlinear boundary damping and memory term and M(s) = 1 + s and f = 0.

We note that stability of problems with the nonlinear term h(∇u) requires a careful treatment because we do not have any information about the influence of the integral about the sign of the derivative E'(t). Although the subject is important, there are few mathematical results in the presence of the nonlinearity given by h(∇u), see [24-26]. In light of this and previous articles [17,22], it is interesting to investigate whether we still have the similar general decay result as in [17] for nondissipative distributed system (1.1) with the memory-type damping acting on a part of the boundary. Hence, the main purpose of this article is to answer the above question for system (1.1)-(1.4). Consequently, by following the arguments close to the one in [17] with necessary modification required the nature of our problem, we establish a general decay result which includes the usual exponential and polynomial decay rates. Furthermore, our results allow a larger class of relax functions which are not necessarily of exponential and polynomial decay. Therefore, this improves earlier results in the literature [22,27].

In order to obtain our results, we consider system (1.1)-(1.4), under some assumptions on a(x), l(t), M and f. Precisely, we state the general assumptions:

(A1) a(x): Ω → R+ is a function.

(A2) f C1(R) is a function and satisfies

(1.11)

where with

d > 0 and .

(A3) M is a C1 function on [0, ∞) satisfying

(1.12)

Where .

(A4) h : Rn R is a C1 function such that ∇h is bounded and there exists β1 > 0 such that

(1.13)

and l(t) is a positive and nonincreasing function.

The remainder of this article is organized as follows. In Section 2, we introduce some notations, present Lemma 2.1 to describe more general relations between the relaxation function g and the corresponding resolvent kernel k and state the existence result to system (1.1)-(1.4). In Section 3, we give the proof of our main result Theorem 3.5.

### 2 Preliminaries

In this section, we introduce some notations and establish the existence of solutions of the problem (1.1)-(1.4). In what follows, let ||·||p denote the usual Lp (Ω) norm , for 1 ≤ p ≤ ∞. We define the convolution product operator by

(2.1)

and set

(2.2)

(2.3)

Using Hölder's inequality, we observe that

(2.4)

Next, we shall use Equation 1.3 to estimate the boundary term . Differentiating (1.3), we obtain

Assume the function k is the resolvent kernel of the relaxation function g, then

Applying Volterra's inverse operator yields

which implies that

(2.5)

where . Reciprocally, taking u0 = 0 on Γ1, identity (2.5) implies (1.3). As we are interested in relaxation functions of more general decay and the resolvent k appeared in Equation 2.5, we want to know if the resolvent k has the same property with the relaxation function g involved in (1.3). The following lemma answers this question. Let h be a relaxation function and k its resolvent kernel, that is

Lemma 2.1. [15,17,22]If h : [0, ∞) → R+ is continuous, then k is also a positive continuous function. Moreover,

(1) If there exists a positive constant c0 such that

where γ : [0, ∞) → R+, is a nonincreasing function satisfying, for some positive constant ε < 1,

Then, k satisfies

(2) Suppose that

for c0 < p - 1. Then, there exists a positive constant ε < 1 such that

where β > 0 is a constant.

Based on this lemma, we will use (2.5) instead of (1.3), i.e., we can consider system (1.1)-(1.4) as follows:

We notice that, due to the condition (1.2), the solution of system (1.1)-(1.4) must belong to the following space:

which endowed with the norm ||∇·||2 is a Hilbert space. Now, we are ready to give the well-posedness of system (1.1)-(1.4).

Theorem 2.2. Let k W2,1 (R+) ∩ W1,∞ (R+), (u0, u1) ∈ (H2 (Ω) ∩ V)2 and satisfy the compatibility condition

Assume further that (A1)-(A4) hold. Then, there exists a unique solution u of system (1.1)-(1.4) such that

Proof. Using the Galerkin method and procedures similar to that of [22,28], we can obtain the result.    □

### 3 Decay of solutions

In this section, we study the asymptotic behavior of the solutions of system (1.1)-(1.4) when the resolvent kernel k satisfies

(3.1)

where γ : [0, ∞) → R+ is a function satisfying the following condition:

(3.2)

To get our result, we further assume that

(3.3)

Let x0 be a fixed point in Rn. Set

and partition the boundary ∂Ω into two sets

(3.4)

Define the first-order energy function of system (1.1)-(1.4) by

(3.5)

The following lemma is associated with the property of the convolution operator, which is used to estimate the energy identity.

Lemma 3.1. If g, ϕ C1(R+), then

(3.6)

Proof. Our conclusion is followed by differentiating the term g ϕ.   □

Lemma 3.2. Under the assumptions of (A1)-(A4), the energy function E(t) satisfies

(3.7)

Proof. Multiplying Equation 1.1 by ut, and integrating by parts over Ω, we get

Exploiting (2.5), (3.6) and the definition of E(t) by (3.5), we have

Then, using Hölder's inequality and Young's inequality, the inequality (3.7) is obtained.    □

Next, we construct a Lyapunov functional which is equivalent to E(t). To do so, for N > 0 large enough, let

(3.8)

where

(3.9)

for 0 < θ < 1.

For the purpose of achieving our main result, we need the following lemmas.

Lemma 3.3. There exist two positive constants α1 and α2 such that the relation

holds for all t ≥ 0.

Proof. From (3.9) and using Young's inequality, we get

where we have used the fact that by (3.5) and

(3.10)

due to by (A3) and (3.5). Here B1 > 0 is the smallest constant such that

(3.11)

Thus, from (3.8), we deduce that

Hence, selecting

(3.12)

there exist two positive constants α1 and α2 such that the relation

holds.    □

Lemma 3.4. Let (A1)-(A4) and (3.1)-(3.3) hold, with β1 (given by (A4)) small enough and

(3.13)

Then, for some t0 large enough, the functional L(t) verifies, along the solution u of (1.1)-(1.4),

(3.14)

for all t t0, where α = min {2θ, 1 - θ} and ci are positive constants given in the proof, i = 4, 5.

Proof. First, we are going to estimate the derivative of ψ(t). From (3.9) and using Equation 1.1, we have

Performing integration by parts and using Young's inequality, we obtain

(3.15)

where ε > 0, cε and c0 are some positive constants. In the following, we will estimate the last two terms on the right-hand side of (3.15). It follows from (1.13), Hölder's inequality, (3.11), (3.3) and (3.10) that

(3.16)

where . Taking (1.11) and (3.4) into account, we have

(3.17)

A substitution of (3.16)-(3.17) into (3.15), we obtain

(3.18)

Now, we analyze the boundary term on the right-hand side of (3.18). Applying Young's inequality and M(λ) ≥ m0 > 0 by (1.12), we have, for ε1 > 0,

where is a positive constant and B* > 0 is the constant such that

(3.19)

Thus, (3.18) becomes

(3.20)

By rewriting the boundary condition (2.5) as

and, then, combining (3.7) and (3.20), we deduce that

Similarly as in deriving (3.16), we note that

(3.21)

where . This implies that

At this point, we choose

Once ε = ε1 is fixed (hence cε and are also fixed), we pick N large satisfying (3.12) and

(3.22)

at the same time. Then, from the properties of k(t) by (3.1) and noting that by (2.4), we see that

Utilizing the inequality by (1.12) and the definition of E(t) by (3.5), we obtain

which together with (3.19) and (3.10) infers that

where α = min{2θ, 1 - θ}. Besides, we note that there exists t0 large enough satisfying

(3.23)

because of limt→∞ k(t) = 0 by (3.13). Therefore, taking β1 small enough such that

(3.24)

then,

(3.25)

for all t t0, where ci are positive constants, i = 4, 5. This completes the proof.    □

Theorem 3.5. Given that (u0, u1) ∈ (H2 (Ω) ∩ V)2, assume that (A1)-(A4), (3.1)-(3.3) and (3.13)hold, with β1 (given by (A4)) small enough. Assume further that

(3.26)

Then, for some t0 large enough, we have, ∀t t0,

(3.27)

otherwise (if u0 ≠ 0 on Γ1),

(3.28)

where a1 is a fixed positive constant and cis a generic positive constant.

Proof. Multiplying (3.25) by γ(t) and exploiting (3.26), (3.1) and (3.7), we derive that

(3.29)

where c6 = c4γ(0) + c5 and . Employing (3.21) again, (3.29) becomes

where

which is equivalent to E(t) due to Lemma 3.3 and γ(t) is nonincreasing by (3.2). In addition to (3.24), we further require

then, we have

(3.30)

where a1 is a positive constant.

Case I: If u0 = 0 on Γ1, then (3.30) reduces to

Integrating the above inequality over (t0, t) to get

Then, using the fact F1(t) is equivalent to E(t), we obtain, for some positive constant c,

Thus, (3.27) is proved.

Case II: If u0 ≠ 0 on Γ1, then (3.30) gives

where . Direct computations give

An integration over (t0, t) yields

Again using the fact F1(t) is equivalent to E(t), we obtain, for some positive constant c,

This completes the proof of Theorem 3.5.    □

### 4 Conclusion and suggestions

Santos et al. [22] considered problem (1.1)-(1.4) with a = 1 and without a function of the gradient term. They showed the solution decays exponentially (or polynomially) to zero provided the kernel decays exponentially (or polynomially) to zero. Recently, Messaoudi and Soufyane in 2010 [17] considered a semi-linear wave equation, in a bounded domain, where the memory-type damping is acting on the boundary. They established a general decay result, from which the usual exponential and polynomial decay rate are only special cases. Motivated by this, we intended to investigate the decay properties of problem (1.1)-(1.4) using the work of Messaaoudi and Soufyane [17]. Since stability of problems with the nonlinear term h(∇u) requires a careful treatment, it is interesting to investigate whether we still have the similar general decay result as that of [16] in the presence of a function of the gradient term. This is our motivation to consider problem (1.1)-(1.4). And, this problem is not considered before.

By adopting and modifying the method proposed by Messaoudi and Soufyane in 2010 [17], we establish a general decay result, from which the usual exponential and polynomial decay rate are only special cases. Further, our result allows certain kernels which are not necessarily of exponential or polynomial decay. In this way, we improved the results of Santos et al. [22], in which they considered problem (1.1)-(1.4) with a = 1 and in the absence of l(t)h (∇u). Moreover, we note that our result also holds for problem (1.1)-(1.4) with a = 1 and l(t) = 0 and without imposing strong damping term, thus our result improves the one of Bae et al. [27]. More precisely, the estimate (3.27) and (3.28) generalizes the exponential and polynomial decay result given in [22,27]. Indeed, we obtain exponential decay for γ(t) = c and polynomial decay for γ(t) = c(1 + t)-1, where c is a positive constant. Additionally, as in [17], our result allows kernels which satisfy k″(t) ≥ c (-k′)1+q, for 0 < q < 1 instead of the usual assumption . It suffices to take, for example, k(t) = (1 + t)-λ, for λ > 0. Direct computations yield

It is clear that , for λ > 0.

Though we consider the conditions on the term involving the gradient are too restrictive and we combine some known ideas to obtain our result, our findings extend those decay results in [22,27] and these findings are interesting to those with closely concerns. For future work, we will consider not necessarily decreasing kernels and relax the condition of h(∇u).

### Competing interests

The author declare that they have no competing interests.

### Acknowledgements

The author would like to thank the anonymous referees for their valuable and constructive suggestions which improve this work.

### References

1. Kirchhoff, G: Vorlesungen über Mechanik. Leipzig Teubner (1883)

2. Biler, P: Remark on the decay for damped string and beam equations. Nonlinear Anal TMA. 9, 839–842 (1984)

3. Brito, EH: Nonlinear initial boundary value problems. Nonlinear Anal TMA. 11, 125–137 (1987). Publisher Full Text

4. Ikerata, R: On the existence of global solutions for some nonlinear hyperbolic equations with Neumann conditions. TRU Math. 24, 1–17 (1988)

5. Ikerata, R: A note on the global solvability of solutions to some nonlinear wave equations with dissipative terms. Differential Integral Equations. 8, 607–616 (1995)

6. Matos, MP, Pereira, DC: On a hyperbolic equation with strong damping. Funkcial Ekvac. 34, 303–311 (1991)

7. Matsuyama, T, Ikerata, R: On global solutions and energy decay for the wave equations of Kirchhoff type with nonlinear damping terms. J Math Anal Appl. 204, 729–753 (1996). Publisher Full Text

8. Nishihara, K, Yamada, Y: On global solutions of some degenerate quasilinear hyperbolic equations with dissipative terms. Funkcial Ekvac. 33, 151–159 (1990)

9. Nishihara, K: Exponentially decay of solutions of some quasilinear hyperbolic equations with linear damping. Nonlinear Anal TMA. 8, 623–636 (1984). Publisher Full Text

10. Ono, K: On global existence, asymptotic stability and blowing up of solutions for some degenerate nonlinear wave equations of Kirchhoff type with a strong dissipation. Math Methods Appl Sci. 20, 151–177 (1997). Publisher Full Text

11. Ono, K: On global solutions and blow-up solutions of nonlinear Kirchhoff strings with nonlinear dissipation. J Math Anal Appl. 216, 321–342 (1997). Publisher Full Text

12. Wu, ST, Tsai, LY: Blow-up of solutions for some nonlinear wave equations of Kirchhoff type with some dissipation. Nonlinear Anal TMA. 65, 243–264 (2006). Publisher Full Text

13. Yamada, Y: On some quasilinear wave equations with dissipative terms. Nagoya Math J. 87, 17–39 (1982)

14. Cavalcanti, MM, Domingos Cavalcanti, VN, Prates, JS, Soriano, JA: Existence and uniform decay of solutions of a degenerate equation with nonlinear boundary memory source term. Nonlinear Anal TMA. 38, 281–294 (1999). Publisher Full Text

15. Cavalcanti, MM, Guesmia, A: General decay rates of solutions to a nonlinear wave equation with boundary condition of memory type. Differential Integral Equations. 18, 583–600 (2005)

16. Messaoudi, SA, Soufyane, A: Boundary stabilization of solutions of a nonlinear system of Timoshenko system. Nonlinear Anal TMA. 67, 2107–2121 (2007). Publisher Full Text

17. Messaoudi, SA, Soufyane, A: General decay of solutions of a wave equation with a boundary control of memory type. Nonlinear Anal Real World Appl. 11, 2896–2904 (2010). Publisher Full Text

18. Cavalcanti, MM, Domingos Cavalcanti, VN, Prates, JS, Soriano, JA: Existence and uniform decay rates for viscoelastic problems with nonlinear boundary damping. Differential Integral Equations. 14, 85–116 (2001)

19. Cavalcanti, MM, Domingos Cavalcanti, VN, Santos, ML: Existence and uniform decay rates of solutions to a degenerate system with memory conditions at the boundary. Appl Math Comput. 150, 439–465 (2004). Publisher Full Text

20. Jiang, S, Muñoz Rivera, JE: A global existence for the Dirichlet problems in nonlinear n-dimensional viscoelasticity. Differential Integral Equations. 9, 791–810 (1996)

21. Santos, ML: Asymptotic behavior of solutions to wave equations with a memory condition at the boundary. Electronic J Differential Equations. 73, 1–11 (2001)

22. Santos, ML, Ferreira, J, Pereira, DC, Raposo, CA: Global existence and stability for wave equation of Kirchhoff type with memory condition at the boundary. Nonlinear Anal TMA. 54, 959–976 (2003). Publisher Full Text

23. Park, JY, Bae, JJ, Hyo Jung, IL: Uniform decay of solution for wave equation of Kirchhoff type with nonlinear boundary damping and memory term. Nonlinear Anal TMA. 50, 871–884 (2002). Publisher Full Text

24. Aassila, M, Cavalcanti, MM, Domingos Cavalcanti, VN: Existence and uniform decay of the wave equation with nonlinear boundary damping and boundary memory source term. Calc Var. 15, 155–180 (2002). Publisher Full Text

25. Guesmia, A: A new approach of stabilization of nondissipative distributed systems. SIAM J Control Optim. 42, 24–52 (2003). Publisher Full Text

26. Park, JY, Ha, TG: Energy decay for nondissipative distributed systems with boundary damping and source term. Nonlinear Anal TMA. 70, 2416–2434 (2009). Publisher Full Text

27. Bae, JJ, Yoon, SB: On uniform decay of wave equation of carrier model subject to memory condition at the boundary. J Korean Math Soc. 44, 1013–1024 (2007). Publisher Full Text

28. Santos, ML, Junior, F: A boundary condition with memory for Kirchhoff plate equations. Appl Math Comput. 148, 475–496 (2004). Publisher Full Text