Skip to main content

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

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:

  1. a.

    u(x,0)≤(≥) u 0 (x) in Ω;

  2. b.

    u(x,t)≤(≥)0 on ∂Ω×(0,T);

  3. 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 1 Let k, α 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 − θ ,

where C is a constant depending on N, m, r, j, k, and q such 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 1 Suppose 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 dx ∫ Ω udx.

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

∫ Ω u m dx ∫ Ω udx≤ | Ω | 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 dx ∫ Ω u s dx.

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

∫ Ω u m dx ∫ Ω u s dx≤ | Ω | 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 2 If 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 dx ∫ Ω udx.
(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 dx ∫ Ω u s dx.
(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)>2m+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 3 Suppose 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 dx ∫ Ω u s dx≤λ|Ω| ∥ 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.

References

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

    Book  Google Scholar 

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

    Google Scholar 

  3. Furter J, Grinfield M: Local vs. nonlocal interactions in populations dynamics. J. Math. Biol. 1989, 27: 65-80. 10.1007/BF00276081

    Article  MathSciNet  Google Scholar 

  4. Chen YP: Blow-up for a system of heat equations with nonlocal sources and absorptions. Comput. Math. Appl. 2004, 48: 361-372. 10.1016/j.camwa.2004.05.002

    Article  MathSciNet  Google Scholar 

  5. Souplet P: Blow-up in nonlocal reaction-diffusion equations. SIAM J. Math. Anal. 1998, 29: 1301-1334. 10.1137/S0036141097318900

    Article  MathSciNet  Google Scholar 

  6. Souplet P: Uniform blow-up profiles and boundary behavior for diffusion equations with nonlocal nonlinear source. J. Differ. Equ. 1999, 153: 374-406. 10.1006/jdeq.1998.3535

    Article  MathSciNet  Google Scholar 

  7. Wang MX, Wang YM: Properties of positive solutions for non-local reaction-diffusion problems. Math. Methods Appl. Sci. 1996, 19: 1141-1156. 10.1002/(SICI)1099-1476(19960925)19:14<1141::AID-MMA811>3.0.CO;2-9

    Article  MathSciNet  Google Scholar 

  8. Galaktionov VA, Levine HA: A general approach to critical Fujita exponents in nonlinear parabolic problems. Nonlinear Anal. 1998, 34: 1005-1027. 10.1016/S0362-546X(97)00716-5

    Article  MathSciNet  Google Scholar 

  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. 2005, 57: 1687-1711. 10.1007/s11253-006-0024-6

    Article  MathSciNet  Google Scholar 

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

    MathSciNet  Google Scholar 

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

    Article  MathSciNet  Google Scholar 

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

    Google Scholar 

  13. Yin JX, Jin CH: Critical extinction and blow-up exponents for fast diffusive Laplacian with sources. Math. Methods Appl. Sci. 2007, 30: 1147-1167. 10.1002/mma.833

    Article  MathSciNet  Google Scholar 

  14. Zhou J, Mu CL: Critical blow-up and extinction exponents for non-Newton polytropic filtration equation with source. Bull. Korean Math. Soc. 2009, 46: 1159-1173. 10.4134/BKMS.2009.46.6.1159

    Article  MathSciNet  Google Scholar 

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

    Google Scholar 

  16. Fang ZB, Li G: Extinction and decay estimates of solutions for a class of doubly degenerate equations. Appl. Math. Lett. 2012, 25: 1795-1802. 10.1016/j.aml.2012.02.020

    Article  MathSciNet  Google Scholar 

  17. Han YZ, Gao WJ: Extinction for a fast diffusion equation with a nonlinear nonlocal source. Arch. Math. 2011, 97: 353-363. 10.1007/s00013-011-0299-1

    Article  MathSciNet  Google Scholar 

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

    Google Scholar 

  19. Fang ZB, Xu XH: Extinction behavior of solutions for the p -Laplacian equations with nonlocal sources. Nonlinear Anal., Real World Appl. 2012, 13: 1780-1789. 10.1016/j.nonrwa.2011.12.008

    Article  MathSciNet  Google Scholar 

  20. Antontsev S, Diaz JI, Shmarev S Progress in Nonlinear Differential Equations and Their Applications. In Energy Methods for Free Boundary Problems: Applications to Nonlinear PDEs and Fluid Mechanics. Birkhäuser, Boston; 2002.

    Chapter  Google Scholar 

  21. Sacks PE: Continuity of solutions of a singular parabolic equation. Nonlinear Anal. 1983, 7: 387-409. 10.1016/0362-546X(83)90092-5

    Article  MathSciNet  Google Scholar 

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

    Book  Google Scholar 

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

    Google Scholar 

Download references

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.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Zhong Bo Fang.

Additional information

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.

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

Xu, X., Fang, Z.B. & Yi, SC. Extinction and decay estimates of solutions for a porous medium equation with nonlocal source and strong absorption. Bound Value Probl 2013, 24 (2013). https://doi.org/10.1186/1687-2770-2013-24

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-24

Keywords