Research

# Extinction and decay estimates of solutions for a porous medium equation with nonlocal source and strong absorption

Xianghui Xu1, Zhong Bo Fang2* and Su-Cheol Yi3

### Author affiliations

1 Department of Mathematics, Pusan National University, Busan, 609-735, Republic of Korea

2 School of Mathematical Sciences, Ocean University of China, Qingdao, 266100, P.R. China

3 Department of Mathematics, Changwon National University, Changwon, 641-773, Republic of Korea

For all author emails, please log on.

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

 Received: 25 September 2012 Accepted: 21 December 2012 Published: 5 March 2013

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

### Abstract

In this paper, we investigate extinction properties of the solutions for the initial Dirichlet boundary value problem of a porous medium equation with nonlocal source and strong absorption terms. We obtain some sufficient conditions for the extinction of nonnegative nontrivial weak solutions and the corresponding decay estimates which depend on the initial data, coefficients, and domains.

### 1 Introduction

We consider the initial Dirichlet boundary value problem for a class of porous medium equations with nonlocal source and strong absorption terms

(1)

(2)

(3)

where 0 < m , k < 1 , d , λ , β , q > 0 , Ω R N ( N 1 ) is a bounded domain with smooth boundary, and u 0 ( x ) L ( Ω ) W 0 1 , 2 ( Ω ) is a nonnegative function. The symbols p and 1 , p , where p 1 , denote L p ( Ω ) - and W 1 , p ( Ω ) -norm, respectively, and | Ω | denotes the measure of Ω.

Equation (1) describes the fast diffusion of concentration of some Newtonian fluids through a porous medium or the density of some biological species in many physical phenomena and biological species theories, while nonlocal source and absorption terms cooperate and interact with each other during the diffusion. It has been known that the nonlocal source term presents a more realistic model for population dynamics; see [1-3]. In the nonlinear diffusion theory, obvious differences exist among the situations of slow ( m > 1 ), fast ( 0 < m < 1 ), and linear ( m = 1 ) diffusions. For example, there is a finite speed propagation in the slow and linear diffusion situations, whereas an infinite speed propagation exists in the fast diffusion situation.

Recently, many scholars have been devoted to the study of blow-up and extinction properties of solutions for nonlinear parabolic equations with nonlocal terms. The blow-up rates and blow-up sets of solutions to equation (1) have been investigated when d = λ = 1 , m 1 , and the linear absorption term is replaced with a nonlinear term with exponent (cf.[4-9]). Extinction is the phenomenon whereby there exists a finite time T > 0 such that the solution is nontrivial for 0 < t < T and then u ( x , t ) 0 for all ( x , t ) Ω × [ T , + ) . In this case, T is called the extinction time. It is also an important property of solutions for nonlinear parabolic equations which have been studied by many researchers. For instance, Evans and Knerr [10] investigated the extinction behaviors of solutions for the Cauchy problem of the semilinear parabolic equation

u t ( x , t ) = Δ u ( x , t ) β ( u ( x , t ) ) , x R n , t > 0 ,

by constructing a suitable comparison function. Ferreira and Vazquez [11] studied the extinction phenomena of solutions for the Cauchy problem of the porous medium equation with an absorption term

u t = ( u m ) x x u p , x R , t > 0 ,

by using the analysis of self-similar solutions. Li and Wu [12] considered the problem of the porous medium equation with a source term

u t = Δ u m + λ u p , x Ω , t > 0 , (4)

subject to (2) and (3). They obtained some conditions for the extinction and non-extinction of solutions to equation (4) and decay estimates by the upper and lower solutions method. On extinctions of solutions to the p-Laplacian equations or the doubly degenerate equations, we refer readers to [13,14] and the references therein.

Replacing the nonlocal term in equation (1) with a local term, Liu [15]et al. considered the initial Dirichlet boundary value problem for a class of porous medium equations

u t = Δ u m + u q β u , x Ω , t > 0 ,

and obtained sufficient conditions for the extinction and non-extinction of solutions to that equation. Thereafter, Fang and Li [16] extended their results to the doubly degenerate equation in the whole dimensional space.

For equation (1) with β = 0 , q > 0 , and N > 2 , Han and Gao [17] showed that q = m is the critical exponent for the occurrence of extinction or non-extinction. When m = 1 , q > 0 , and N > 2 , the conditions for the extinction and non-extinction of solutions and corresponding decay estimates were obtained (cf.[18]). Recently, Fang and Xu [19] considered equation (1) with k = 1 , when the diffusion term was replaced with a p-Laplacian operator in the whole dimensional space, and showed that the extinction of the weak solution is determined by the competition of two nonlinear terms. They also obtained the exponential decay estimates which depend on the initial data, coefficients, and domains. The extinctions of solutions to equation (1) with nonlocal source terms do not depend on the first eigenvalue of the corresponding operator, which is different from the case of local source terms. The extinction and decay estimates for solutions to the nonlocal fast diffusion equations with nonzero coefficients and strong absorption terms, like equation (1), are still being investigated.

Motivated by the above works, we investigate whether the existence of strong absorption can change extinction behaviors for solutions to problem (1)-(3) in the whole dimensional space. The main tools we use are the integral estimate method and the Gagliardo-Nirenberg inequality. This technique has a wide application, especially for equations that do not satisfy the maximum principle (cf.[20]). Our goals are to show that the extinction of nonnegative nontrivial weak solutions to problem (1)-(3) occurs when 0 < k q < 1 and to find the decay estimates depending on the initial data, coefficients, and domains.

Our paper is organized as follows. In Section 2, we give preliminary knowledge including lemmas that are required in the proofs of our results and present the proofs for the results in Section 3.

### 2 Preliminary knowledge

Due to the singularity of equation (1), problem (1)-(3) has no classical solutions in general. To state the definition of the weak solution, we let Q T = Ω × ( 0 , T ) and firstly define the class of nonnegative testing functions

F = { ξ : ξ C ( Q ¯ T ) C 2 , 1 ( Q T ) , ξ t , Δ ξ L 2 ( Q T ) ; ξ 0 , ξ | Ω × ( 0 , T ) = 0 } .

Definition 1 A function u L ( Q T ) is called a weak subsolution (supersolution) of problem (1)-(3) in Q T if the following conditions hold:

a. u ( x , 0 ) ( ) u 0 ( x ) in Ω;

b. u ( x , t ) ( ) 0 on Ω × ( 0 , T ) ;

c. For every t ( 0 , T ) and every ξ F ,

Ω u ( x , t ) ξ ( x , t ) d x Ω u 0 ( x ) ξ ( x , 0 ) d x + 0 T Ω { u ξ s + d u m Δ ξ + λ Ω u q ( y , s ) d y ξ ( x , s ) β u k ( x , s ) ξ ( x , s ) } d x d s .

A function u is called a locally weak solution of problem (1)-(3) if it is both a subsolution and a supersolution for some T > 0 .

Remark 1 The existence and uniqueness of locally nonnegative solutions in time to problem (1)-(3) can be obtained by the standard parabolic regular theory that can be applied to get suitable estimates in the standard limiting process (cf.[2,21,22]). The proof is similar to the ones in the cited references, and so it is omitted here.

Lemma 1Letk, αbe positive constants and k < 1 . If y ( t ) is a nonnegative absolutely continuous function on [ 0 , + ) satisfying the problem

d y d t + α y k 0 , t 0 ; y ( 0 ) 0 ,

then we have the decay estimate

y ( t ) [ y 1 k ( 0 ) α ( 1 k ) t ] 1 1 k , t [ 0 , T ) , y ( t ) = 0 , t [ T , + ) ,

where T = y 1 k ( 0 ) α ( 1 k ) .

Proof

By solving the initial problem

d y d t + α y k = 0 , t > 0 ; y | t = 0 = y ( 0 ) ,

and using the comparison principle, one can easily obtain the result. □

Lemma 2 (The Gagliardo-Nirenberg inequality) [23]

Suppose that u W 0 k , m ( Ω ) , 1 m + , 0 j < k , and 1 1 r 1 m k N . We then have the inequality

D j u q C D k u m θ u r 1 θ ,

whereCis a constant depending onN, m, r, j, k, andqsuch that 1 q = j N + θ ( 1 m k N ) + 1 θ r and θ [ 0 , 1 ) . While if m < N k j , then q [ N r N + r j , N m N ( k j ) m ] , and if m N k j , then q [ N r N + r j , + ] .

### 3 Main results

In this section, we give some extinction properties of nonnegative nontrivial weak solutions of problem (1)-(3) stated in the following theorems. The corresponding decay estimates to the solutions will be presented in the proofs of the theorems for brief expressions instead of in the statements.

Theorem 1Suppose that 0 < k < 1 and 0 < m = q < 1 . Then the nonnegative nontrivial weak solution of problem (1)-(3) vanishes in finite time for any nonnegative initial data provided that either | Ω | orλis sufficiently small.

Proof We first consider the case that N = 1  or  2 . Multiplying both sides of (1) by u and integrating the result over Ω, we have

1 2 d d t u 2 2 + 4 m d ( m + 1 ) 2 u m + 1 2 2 2 + β u k + 1 k + 1 = λ Ω u m d x Ω u d x .

By Hölder’s inequality, we get the inequality

Ω u m d x Ω u d x | Ω | 2 s 1 1 m s 1 u s 1 m + 1 .

In particular, if s 1 = 2 , we get the inequality

1 2 d d t u 2 2 + 4 m d ( m + 1 ) 2 u m + 1 2 2 2 + β u k + 1 k + 1 λ | Ω | 3 m 2 u 2 m + 1 , (5)

from the two expressions above. By using the Sobolev embedding inequality, one can show that there exists an embedding constant γ ( N , Ω ) > 0 such that

u m + 1 2 s 2 γ ( N , Ω ) u m + 1 2 2 ,

i.e.,

u ( m + 1 ) s 2 2 m + 1 γ ( N , Ω ) 2 u m + 1 2 2 2 .

In particular, if s 2 = 4 m + 1 , then the inequality above turns out to be

γ ( N , Ω ) 2 u 2 m + 1 u m + 1 2 2 2 . (6)

Since 0 < m < 1 , we get 4 m + 1 2 , and hence, inequality (5) becomes

1 2 d d t u 2 2 + [ 4 m d ( m + 1 ) 2 λ | Ω | 3 m 2 γ ( N , Ω ) 2 ] u m + 1 2 2 2 + β u k + 1 k + 1 0 . (7)

By Lemma 2, we get the inequality

u 2 C ( N , m , k ) u k + 1 1 θ 1 u m + 1 2 2 2 θ 1 m + 1 , (8)

where θ 1 = m + 1 2 ( 1 k + 1 1 2 ) ( 1 N 1 2 + m + 1 2 1 k + 1 ) 1 = N ( 1 k ) ( m + 1 ) 2 [ 2 ( k + 1 ) + N ( m k ) ] . Since N = 1  or  2 , and 0 < k < 1 , it can be easily seen that 0 < θ 1 < 1 .

It then follows from (8) and Young’s inequality that

u 2 k 1 C ( N , m , k ) k 1 u k + 1 k 1 ( 1 θ 1 ) u m + 1 2 2 2 k 1 θ 1 m + 1 C ( N , m , k ) k 1 ( η 1 u m + 1 2 2 2 + C ( η 1 ) u k + 1 k 1 ( 1 θ 1 ) ( m + 1 ) m + 1 k 1 θ 1 ) , (9)

where k 1 > 1 and η 1 > 0 will be determined later. If we choose

k 1 = ( k + 1 ) ( m + 1 ) ( 1 θ 1 ) ( m + 1 ) + θ 1 ( k + 1 ) = 4 ( k + 1 ) + 2 N ( m k ) 4 + N ( m k ) ,

then 1 < k 1 < 2 and k 1 ( 1 θ 1 ) ( m + 1 ) m + 1 k 1 θ 1 = k + 1 . From (9) we have

C ( N , m , k ) k 1 β C ( η 1 ) u 2 k 1 η 1 β C ( η 1 ) u m + 1 2 2 2 + β u k + 1 k + 1 . (10)

By inequalities (7) and (10), we get the inequality

1 2 d d t u 2 2 + [ 4 m d ( m + 1 ) 2 λ | Ω | 3 m 2 γ ( N , Ω ) 2 η 1 β C ( η 1 ) ] u m + 1 2 2 2 + C ( N , m , k ) k 1 β C ( η 1 ) u 2 k 1 0 .

Here, we can choose η 1 and λ or | Ω | small enough so that

4 m d ( m + 1 ) 2 λ | Ω | 3 m 2 γ ( N , Ω ) 2 η 1 β C ( η 1 ) 0 .

Setting C 1 = C ( N , m , k ) k 1 β C ( η 1 ) , we have d d t u 2 + C 1 u 2 k 1 1 0 .

By Lemma 1, we then obtain

u 2 [ u 0 2 2 k 1 C 1 ( 2 k 1 ) t ] 1 2 k 1 , t [ 0 , T 1 ) , u 2 = 0 , t [ T 1 , + ) ,

where T 1 = u 0 2 2 k 1 C 1 ( 2 k 1 ) , which give the decay estimates in finite time for N = 1  or  2 .

Secondly, we consider the case that N > 2 . If N 2 N + 2 m < 1 , multiplying both sides of (1) by u s ( s > m ) and integrating the result over Ω, we get

1 s + 1 d d t u s + 1 s + 1 + 4 m s d ( m + s ) 2 u m + s 2 2 2 + β u k + s k + s = λ Ω u m d x Ω u s d x .

By Hölder’s inequality, we have the inequality

Ω u m d x Ω u s d x | Ω | 2 s 3 m s s 3 u s 3 m + s .

In particular, if s 3 = s + 1 , we then get the inequality

1 s + 1 d d t u s + 1 s + 1 + 4 m s d ( m + s ) 2 u m + s 2 2 2 + β u k + s k + s λ | Ω | s m + 2 s + 1 u s + 1 m + s (11)

by the two expressions above. By the Sobolev embedding inequality, one can see that there exists an embedding constant C 0 > 0 such that

u m + s 2 2 N N 2 2 C 0 2 u m + s 2 2 2 . (12)

Using Hölder’s inequality again, we have the inequality

u s + 1 m + s | Ω | m + s s + 1 N 2 N u m + s 2 2 N N 2 2 . (13)

From inequalities (11), (12), and (13), we then obtain the inequality

1 s + 1 d d t u s + 1 s + 1 + [ 4 m s d ( m + s ) 2 λ C 0 2 | Ω | 1 + 2 N ] u m + s 2 2 2 + β u k + s k + s 0 . (14)

By Lemma 2, we can also have

u s + 1 C ( N , k , s ) u k + s 1 θ 2 u m + s 2 2 2 θ 2 m + s , (15)

where θ 2 = m + s 2 ( 1 k + s 1 s + 1 ) ( 1 N 1 2 + m + s 2 1 k + s ) 1 = N ( 1 k ) ( m + s ) ( s + 1 ) [ 2 ( k + s ) + N ( m k ) ] . Since N 2 N + 2 m < 1 and 0 < k < 1 , one can easily see that 0 < θ 2 < 1 . Then it follows from (15) and Young’s inequality that

u s + 1 k 2 C ( N , k , s ) k 2 u k + s k 2 ( 1 θ 2 ) u m + s 2 2 2 k 2 θ 2 m + s C ( N , k , s ) k 2 ( η 2 u m + s 2 2 2 + C ( η 2 ) u k + s k 2 ( 1 θ 2 ) ( m + s ) m + s k 2 θ 2 ) , (16)

where k 2 > 0 and η 2 > 0 will be determined later. If we choose

k 2 = ( k + s ) ( m + s ) ( 1 θ 2 ) ( m + s ) + θ 2 ( k + s ) = 2 ( s + 1 ) ( k + s ) + N ( s + 1 ) ( m k ) 2 ( s + 1 ) + N ( m k ) ,

then s < k 2 < s + 1 and k 2 ( 1 θ 2 ) ( m + s ) m + s k 2 θ 2 = k + s . We then have the inequality

C ( N , k , s ) k 2 β C ( η 2 ) u s + 1 k 2 η 2 β C ( η 2 ) u m + s 2 2 2 + β u k + s k + s (17)

by (16). From inequalities (14) and (17), we can also obtain the inequality

1 s + 1 d d t u s + 1 s + 1 + [ 4 m s d ( m + s ) 2 λ C 0 2 | Ω | 1 + 2 N η 2 β C ( η 2 ) ] u m + s 2 2 2 + C ( N , k , s ) k 2 β C ( η 2 ) u s + 1 k 2 0 .

Here, we can choose η 2 and λ or | Ω | small enough so that 4 m s d ( m + s ) 2 λ C 0 2 | Ω | 1 + 2 N η 2 β C ( η 2 ) 0 . Setting C 2 = C ( N , k , s ) k 2 β C ( η 2 ) , we have d d t u s + 1 + C 2 u s + 1 k 2 s 0 from the inequality above. By Lemma 1, we obtain that

u s + 1 [ u 0 s + 1 s + 1 k 2 C 2 ( s + 1 k 2 ) t ] 1 s + 1 k 2 , t [ 0 , T 2 ) , u s + 1 = 0 , t [ T 2 , + ) ,

where T 2 = u 0 s + 1 s + 1 k 2 C 2 ( s + 1 k 2 ) , which give the decay estimates in finite time for N > 2 such that N 2 N + 2 m < 1 .

If 0 < m < N 2 N + 2 , one can show that there exists an embedding constant C 00 > 0 such that

u s + 1 m + s | Ω | 1 ( s + 1 ) ( N + 2 ) N ( m + s ) u N ( m + s ) N 2 m + s C 00 2 | Ω | 1 ( s + 1 ) ( N + 2 ) N ( m + s ) u m + s 2 2 2

by multiplying both sides of (1) by u s ( s > N ( 1 m ) 2 1 ) and integrating the result over Ω, and the Sobolev embedding inequality. By using the inequality above and a similar argument as above, the following decay estimates can be obtained:

u s + 1 [ u 0 s + 1 s + 1 k 2 C ( N , k , s ) k 2 β C ( η 2 ( s + 1 k 2 ) t ] 1 s + 1 k 2 , t ( 0 , T 2 ) , u s + 1 = 0 , t [ T 2 , + ) ,

provided that

4 m s d ( m + s ) 2 λ C 00 2 | Ω | 1 ( s + 1 ) ( N + 2 ) N ( m + s ) + s m + 2 s + 1 η 2 β C ( η 2 ) 0 ,

where T 2 = C ( η 2 ) u 0 s + 1 s + 1 k 2 C ( N , k , s ) k 2 β ( s + 1 k 2 ) . □

Theorem 2If 0 < k < 1 , then the nonnegative nontrivial weak solution of problem (1)-(3) vanishes in finite time provided that u 0 , | Ω | orλis sufficiently small, and q > 2 k ( s + 1 ) + N ( m k ) 2 ( s + 1 ) + N ( m k ) , where if N = 1  or  2 , then s = 1 , and if N > 2 , then s > m .

Proof Assume that q 1 . If N = 1  or  2 , multiplying both sides of (1) by u and integrating the result over Ω, we have the equation

1 2 d d t u 2 2 + 4 m d ( m + 1 ) 2 u m + 1 2 2 2 + β u k + 1 k + 1 = λ Ω u q d x Ω u d x . (18)

By (10) and (18), and using Hölder’s inequality, we get the inequality

Choosing η 1 small enough so that 4 m d ( m + 1 ) 2 η 1 β C ( η 1 ) 0 , we obtain the inequality

d d t u 2 + u 2 k 1 1 [ C ( N , m , k ) k 1 β C ( η 1 ) λ | Ω | 3 q 2 u 2 q k 1 + 1 ] 0 .

Hence, we have the inequality

d d t u 2 + C 3 u 2 k 1 1 0

provided that

u 0 2 < [ C ( N , m , k ) k 1 β C ( η 1 ) λ | Ω | 3 q 2 ] 1 q k 1 + 1 ,

and

q > k 1 1 = 4 k + N ( m k ) 4 + N ( m k ) ,

where C 3 = C ( N , m , k ) k 1 β C ( η 1 ) λ | Ω | 3 q 2 u 0 2 q k 1 + 1 > 0 , from which and a similar argument as the one used in the proof of Theorem 1, the following decay estimates can be obtained:

u 2 [ u 0 2 2 k 1 C 3 ( 2 k 1 ) t ] 1 2 k 1 , t [ 0 , T 3 ) , u 2 = 0 , t [ T 3 , + ) ,

where T 3 = u 0 2 2 k 1 C 3 ( 2 k 1 ) .

If N > 2 and 0 < m < 1 , multiplying both sides of (1) by u s ( s > m ) and integrating the result over Ω, we get the equation

1 s + 1 d d t u s + 1 s + 1 + 4 m s d ( m + s ) 2 u m + s 2 2 2 + β u k + s k + s = λ Ω u q d x Ω u s d x . (19)

By (17) and (19), and using Hölder’s inequality, we obtain the inequality

1 s + 1 d d t u s + 1 s + 1 + [ 4 m s d ( m + s ) 2 η 2 β C ( η 2 ) ] u m + s 2 2 2 + C ( N , k , s ) k 2 β C ( η 2 ) u s + 1 k 2 λ | Ω | s q + 2 s + 1 u s + 1 q + s .

Choosing η 2 small enough so that 4 m s d ( m + s ) 2 η 2 β C ( η 2 ) 0 , we have

d d t u s + 1 + u s + 1 k 2 s [ C ( N , k , s ) k 2 β C ( η 2 ) λ | Ω | s q + 2 s + 1 u s + 1 q k 2 + s ] 0 .

Therefore, we obtain the inequality

d d t u s + 1 + C 4 u s + 1 k 2 s 0

provided that

u 0 s + 1 < [ C ( N , k , s ) k 2 β C ( η 2 ) λ | Ω | s q + 2 s + 1 ] 1 q k 2 + s ,

and

q > k 2 s = 2 k ( s + 1 ) + N ( m k ) 2 ( s + 1 ) + N ( m k ) ,

where C 4 = C ( N , k , s ) k 2 β C ( η 2 ) λ | Ω | s q + 2 s + 1 u 0 s + 1 q k 2 + s > 0 , which yields the following decay estimates:

u s + 1 [ u 0 s + 1 s + 1 k 2 C 4 ( s + 1 k 2 ) t ] 1 s + 1 k 2 , t [ 0 , T 4 ) , u s + 1 = 0 , t [ T 4 , + ) ,

where T 4 = u 0 s + 1 s + 1 k 2 C 4 ( s + 1 k 2 ) .

Since s > m , we have 2 ( s + 1 ) > 2 m + 2 , and hence, if k m , then q > k 2 s m .

Assume that q > 1 . If λ 1 is the first eigenvalue of the boundary problem

Δ ψ ( x ) = λ ψ ( x ) , x Ω ; ψ ( x ) = 0 , x Ω ,

and φ 1 ( x ) 0 , φ 1 = 1 , is an eigenfunction corresponding to the eigenvalue λ 1 , then for sufficiently small a > 0 , it can be easily shown that a φ 1 1 m ( x ) is an upper solution of problem (1)-(3) provided that u 0 ( x ) a φ 1 1 m ( x ) , x Ω . We then have u ( x , t ) a φ 1 ( x ) for t > 0 by the comparison principle. Therefore, from equation (19), we can obtain the inequality

1 s + 1 d d t u s + 1 s + 1 + [ 4 m s d ( m + s ) 2 λ a q m C 0 2 | Ω | 1 + 2 N η 2 β C ( η 2 ) ] u m + s 2 2 2 + C ( N , k , s ) k 2 β C ( η 2 ) u s + 1 k 2 0 ,

from which the following decay estimates can be obtained:

u s + 1 [ u 0 s + 1 s + 1 k 2 C ( N , k , s ) k 2 β C ( η 2 ( s + 1 k 2 ) t ] 1 s + 1 k 2 , t ( 0 , T 4 ) , u s + 1 = 0 , t [ T 4 , + ) ,

provided that

4 m s d ( m + s ) 2 λ a q m C 0 2 | Ω | 1 + 2 N η 2 β C ( η 2 ) 0 ,

where T 4 = C ( η 2 ) u 0 s + 1 s + 1 k 2 C ( N , k , s ) k 2 β ( s + 1 k 2 ) . □

Remark 2 Since the Sobolev embedding inequality cannot be used in the proof of Theorem 2, it is not necessary to consider the cases that N 2 N + 2 m < 1 and 0 < m < N 2 N + 2 , when N > 2 . In addition, if k m , the conditions in Theorem 2 imply that q > m .

Theorem 3Suppose that 0 < k < 1 and m > q k . Then the nonnegative nontrivial weak solution of problem (1)-(3) vanishes in finite time for any nonnegative initial data provided thatβis sufficiently large.

Proof We first consider the case that N = 1  or  2 . Multiplying both sides of (1) by u and integrating the result over Ω, and using Hölder’s inequality, we get

1 2 d d t u 2 2 + 4 m d ( m + 1 ) 2 u m + 1 2 2 2 + β u k + 1 k + 1 = λ Ω u q d x Ω u d x λ | Ω | u q + 1 q + 1 .

By Lemma 2, we have the inequality

u q + 1 C ( N , k , q ) u k + 1 1 θ 3 u m + 1 2 2 2 θ 3 m + 1 , (20)

where θ 3 = m + 1 2 ( 1 k + 1 1 q + 1 ) ( 1 N 1 2 + m + 1 2 1 k + 1 ) 1 = N ( q k ) ( m + 1 ) ( q + 1 ) [ 2 ( k + 1 ) + N ( m k ) ] [ 0 , 1 ) . Since q < m , we have m + 1 ( q + 1 ) θ 3 > 0 . It then follows from (20) and Young’s inequality that

λ | Ω | u q + 1 q + 1 λ | Ω | C ( N , k , q ) q + 1 u k + 1 ( q + 1 ) ( 1 θ 3 ) u m + 1 2 2 2 ( q + 1 ) θ 3 m + 1 λ | Ω | C ( N , k , q ) q + 1 ( η 3 u m + 1 2 2 2 + C ( η 3 ) u k + 1 ( q + 1 ) ( 1 θ 3 ) ( m + 1 ) m + 1 ( q + 1 ) θ 3 ) , (21)

where η 3 will be determined later. From (18) and (21), one can see that

1 2 d d t u 2 2 + [ 4 m d ( m + 1 ) 2 η 3 λ | Ω | C ( N , k , q ) q + 1 ] u m + 1 2 2 2 + β u k + 1 k + 1 C ( η 3 ) λ | Ω | C ( N , k , q ) q + 1 u k + 1 ( q + 1 ) ( 1 θ 3 ) ( m + 1 ) m + 1 ( q + 1 ) θ 3 .

We then obtain the inequality

1 2 d d t u 2 2 + [ 4 m d ( m + 1 ) 2 η 3 λ | Ω | C ( N , k , q ) q + 1 ] γ 2 u 2 m + 1 + u k + 1 k + 1 [ β C ( η 3 ) λ | Ω | C ( N , k , q ) q + 1 u k + 1 α 1 ] 0

by (6) and the inequality above, where α 1 = ( q + 1 ) ( 1 θ 3 ) ( m + 1 ) m + 1 ( q + 1 ) θ 3 ( k + 1 ) = 2 ( q k ) ( k + 1 ) 2 ( k + 1 ) + N ( m q ) 0 . We can choose η 3 small enough so that C 5 = 4 m d ( m + 1 ) 2 η 3 λ | Ω | C ( N , k , q ) q + 1 > 0 . Once η 3 is fixed, we may choose β large enough so that

β C ( η 3 ) λ | Ω | C ( N , k , q ) q + 1 u k + 1 α 1 0 .

Hence, we have the inequality

d d t u 2 + C 5 u 2 m 0 ,

from which the following decay estimates can be obtained by a similar argument as the one used in the proof of Theorem 1:

u 2 [ u 0 2 1 m C 5 ( 1 m ) t ] 1 1 m , t [ 0 , T 5 ) , u 2 = 0 , t [ T 5 , + ) ,

where T 5 = u 0 2 1 m C 5 ( 1 m ) .

Secondly, we consider the case that N > 2 . If N 2 N + 2 m < 1 , multiplying both sides of (1) by u s ( s > m ) and integrating the result over Ω, and then using Hölder’s inequality, we get

1 s + 1 d d t u s + 1 s + 1 + 4 m s d ( m + s ) 2 u m + s 2 2 2 + β u k + s k + s = λ Ω u q d x Ω u s d x λ | Ω | u q + s q + s .

By Lemma 2, it can be shown that

u q + s C ( N , k , q , s ) u k + s 1 θ 4 u m + s 2 2 2 θ 4 m + s , (22)

where θ 4 = m + s 2 ( 1 k + s 1 q + s ) ( 1 N 1 2 + m + s 2 1 k + s ) 1 = N ( q k ) ( m + s ) ( q + s ) [ 2 ( k + s ) + N ( m k ) ] [ 0 , 1 ) . Since q < m , we have m + s ( q + s ) θ 4 > 0 . It then follows from (22) and Young’s inequality that

λ | Ω | u q + s q + s λ | Ω | C ( N , k , q , s ) q + s u k + s ( q + s ) ( 1 θ 4 ) u m + s 2 2 2 ( q + s ) θ 4 m + s λ | Ω | C ( N , k , q , s ) q + s ( η 4 u m + s 2 2 2 + C ( η 4 ) u k + s ( q + s ) ( 1 θ 4 ) ( m + s m + s ( q + s ) θ 4 ) , (23)

where η 4 will be determined later. From (19) and (23), one can see that

1 s + 1 d d t u s + 1 s + 1 + [ 4 m s d ( m + s ) 2 η 4 λ | Ω | C ( N , k , q , s ) q + s ] u m + s 2 2 2 + β u k + s k + s C ( η 4 ) λ | Ω | C ( N , k , q , s ) q + s u k + s ( q + s ) ( 1 θ 4 ) ( m + s ) m + s ( q + s ) θ 4 .

By (12), (13), and the inequality above, we can obtain the inequality

1 s + 1 d d t u s + 1 s + 1 + [ 4 m s d ( m + s ) 2 η 4 λ | Ω | C ( N , k , q , s ) q + s ] C 0 2 | Ω | N 2 N m + s s + 1 u s + 1 m + s + u k + s k + s [ β C ( η 4 ) λ | Ω | C ( N , k , q , s ) q + s u k + s α 2 ] 0 ,

where α 2 = ( q + s ) ( 1 θ 4 ) ( m + s ) m + s ( q + s ) θ 4 ( k + s ) = 2 ( q k ) ( k + s ) 2 ( k + s ) + N ( m q ) 0 . We can choose η 4 small enough so that C 6 = [ 4 m s d ( m + s ) 2 η 4 λ | Ω | C ( N , k , q , s ) q + s ] C 0 2 | Ω | N 2 N m + s s + 1 > 0 . Once η 4 is fixed, we can choose β large enough so that

β C ( η 4 ) λ | Ω | C ( N , k , q , s ) q + s u k + s α 2 0 .

Hence, we can obtain the inequality

d d t u s + 1 + C 6 u s + 1 m 0 ,

from which the following decay estimates can be obtained:

u s + 1 [ u 0 s + 1 1 m C 6 ( 1 m ) t ] 1 1 m , t [ 0 , T 6 ) , u s + 1 = 0 , t [ T 6 , + ) ,

where T 6 = u 0 s + 1 1 m C 6 ( 1 m ) .

Similarly, one can obtain the following decay estimates for 0 < m < N 2 N + 2 :

where T 6 = u 0 s + 1 1 m [ 4 m s d ( m + s ) 2 η 4 λ | Ω | C ( N , k , q , s ) q + s ] C 00 2 | Ω | 1 ( s + 1 ) ( N + 2 ) N ( m + s ) ( 1 m ) . □

Remark 3 One can see from Theorems 1-3 that the extinction of nonnegative nontrivial weak solutions to problem (1)-(3) occurs when 0 < k q < 1 .

Remark 4 Theorems 1-3 all require | Ω | , λ, or u 0 to be sufficiently small or β to be sufficiently large.

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

All authors contributed equally to the manuscript and read and approved the final manuscript.

### Acknowledgements

The second and third authors were supported by the National Science Foundation of Shandong Province of China (ZR2012AM018) and Changwon National University in 2012, respectively. The authors would like to express their sincere gratitude to the anonymous reviewers for their insightful and constructive comments.

### References

1. Bebernes, J, Eberly, D: Mathematical Problems from Combustion Theory, Springer, New York (1989)

2. Pao, CV: Nonlinear Parabolic and Elliptic Equations, Plenum, New York (1992)

3. Furter, J, Grinfield, M: Local vs. nonlocal interactions in populations dynamics. J. Math. Biol.. 27, 65–80 (1989). Publisher Full Text

4. Chen, YP: Blow-up for a system of heat equations with nonlocal sources and absorptions. Comput. Math. Appl.. 48, 361–372 (2004). Publisher Full Text

5. Souplet, P: Blow-up in nonlocal reaction-diffusion equations. SIAM J. Math. Anal.. 29, 1301–1334 (1998). Publisher Full Text

6. Souplet, P: Uniform blow-up profiles and boundary behavior for diffusion equations with nonlocal nonlinear source. J. Differ. Equ.. 153, 374–406 (1999). Publisher Full Text

7. Wang, MX, Wang, YM: Properties of positive solutions for non-local reaction-diffusion problems. Math. Methods Appl. Sci.. 19, 1141–1156 (1996). Publisher Full Text

8. Galaktionov, VA, Levine, HA: A general approach to critical Fujita exponents in nonlinear parabolic problems. Nonlinear Anal.. 34, 1005–1027 (1998). Publisher Full Text

9. Afanas’eva, NV, Tedeev, AF: Theorems on the existence and nonexistence of solutions of the Cauchy problem for degenerate parabolic equations with nonlocal source. Ukr. Math. J.. 57, 1687–1711 (2005). Publisher Full Text

10. Evans, LC, Knerr, BF: Instantaneous shrinking of the support of nonnegative solutions to certain nonlinear parabolic equations and variational inequalities. Ill. J. Math.. 23, 153–166 (1979)

11. Ferreira, R, Vazquez, JL: Extinction behavior for fast diffusion equations with absorption. Nonlinear Anal.. 43, 353–376 (2001)

12. Li, YX, Wu, JC: Extinction for fast diffusion equations with nonlinear sources. Electron. J. Differ. Equ.. 2005, Article ID 23 (2005)

13. Yin, JX, Jin, CH: Critical extinction and blow-up exponents for fast diffusive Laplacian with sources. Math. Methods Appl. Sci.. 30, 1147–1167 (2007). Publisher Full Text

14. Zhou, J, Mu, CL: Critical blow-up and extinction exponents for non-Newton polytropic filtration equation with source. Bull. Korean Math. Soc.. 46, 1159–1173 (2009). Publisher Full Text

15. Liu, WJ, Wang, MX, Wu, B: Extinction and decay estimates of solutions for a class of porous medium equations. J. Inequal. Appl.. 2007, Article ID 087650 (2007)

Article ID 087650

PubMed Abstract | PubMed Central Full Text

16. Fang, ZB, Li, G: Extinction and decay estimates of solutions for a class of doubly degenerate equations. Appl. Math. Lett.. 25, 1795–1802 (2012). Publisher Full Text

17. Han, YZ, Gao, WJ: Extinction for a fast diffusion equation with a nonlinear nonlocal source. Arch. Math.. 97, 353–363 (2011). Publisher Full Text

18. Liu, WJ: Extinction and non-extinction of solutions for a nonlocal reaction-diffusion problem. Electron. J. Qual. Theory Differ. Equ.. 2010, Article ID 15 (2010)

19. Fang, ZB, Xu, XH: Extinction behavior of solutions for the p-Laplacian equations with nonlocal sources. Nonlinear Anal., Real World Appl.. 13, 1780–1789 (2012). Publisher Full Text

20. Antontsev, S, Diaz, JI, Shmarev, S: Energy Methods for Free Boundary Problems: Applications to Nonlinear PDEs and Fluid Mechanics, Birkhäuser, Boston (2002)

21. Sacks, PE: Continuity of solutions of a singular parabolic equation. Nonlinear Anal.. 7, 387–409 (1983). Publisher Full Text

22. Wu, ZQ, Zhao, JN, Yin, JX, Li, HL: Nonlinear Diffusion Equations, World Scientific, River Edge (2001)

23. Ladyzhenskaya, OA, Solonnikov, VA, Ural’tseva, NN: Linear and Quasilinear Equations of Parabolic Type, Am. Math. Soc., Providence (1968)