Skip to main content

Global existence and blow-up of solutions for p-Laplacian evolution equation with nonlinear memory term and nonlocal boundary condition

Abstract

In this paper, we deal with an initial boundary value problem for a p-Laplacian evolution equation with nonlinear memory term and inner absorption term subject to a weighted linear nonlocal boundary condition. We find the effects of a weighted function as regards determining blow-up of nonnegative solutions or not and establish the precise blow-up estimate for the linear diffusion case under some suitable conditions.

1 Introduction

In the past decades, there have been many works dealing with global existence and blow-up properties of solutions for nonlinear parabolic equations, especially the initial boundary value problems with nonlocal terms in equations or boundary conditions, we refer to [110] and references therein. For the study of an initial boundary value problem for local parabolic equations with nonlocal boundary condition, we refer to [14]. For example, Friedman [1] studied the linear parabolic equation

u t Au=0,(x,t)Ω×(0,T),

subject to the weighted linear nonlocal Dirichlet boundary condition

u(x,t)= Ω f(x,y)u(y,t)dy,(x,t)Ω×(0,T),
(1.1)

where A is an elliptic operator,

A= i , j = 1 n a i , j (x) 2 x i x j + i = 1 n b i (x) x i +c(x)andc(x)0,

and the nonnegative continuous function f(x,y) satisfies suitable conditions. He proved that when Ω f(x,y)dyρ<1, the solution approaches to 0 monotonously and exponentially as t. As regards more general discussions on an initial boundary value problem for a linear parabolic equation with a weighted linear nonlocal Neumann boundary condition, one can refer to [2] by Pao, where the following problem:

u t L u = g ( x , u ) , x Ω , t > 0 , B u = Ω f ( x , y ) u ( y , t ) d y , x Ω , t > 0 , u ( x , 0 ) = u 0 ( x ) , x Ω ,

where

Lu= i , j = 1 n a i j (x) u x i x j + j = 1 n b j (x) u x j ,Bu= α 0 u n +u,

was considered. He studied the asymptotic behavior of solutions and found the influence of the weight function on the existence of global and blow-up solutions. Later, Akila [3] adopted the method of an upper-lower solution to consider the semilinear parabolic equation

u t =Δu+g(x,t,u),(x,t)Ω×(0,),

under a similar weighted linear nonlocal boundary condition. Wang et al. [4] studied a porous medium equation with power form source term,

u t = u m + u p ,(x,t)Ω×(0,+),

under the weighted linear nonlocal Dirichlet boundary condition (1.1). By virtue of the method of an upper-lower solution, they obtained global existence, blow-up properties, and blow-up rate of solutions.

For the study of initial boundary value problem with nonlocal parabolic equation, especially the nonlocal problem with time-integral, we refer to [510]. Under a homogeneous Dirichlet boundary condition, Li and Xie [5] studied the nonlinear diffusion equation

u t u= u q 0 t u p ds,(x,t)Ω×(0,T),

where p0, q>0. They obtained the sufficient conditions of global existence and blow-up of solutions under appropriate critical conditions. Furthermore, under the following assumptions:

there exists  t 0 (0,T) such that  u t (x, t 0 )0 for all x Ω ¯ ,
(1.2)

and

Ω = { x R N : | x | < R } , u 0 ( x ) = u 0 ( | x | ) u 0 ( r ) , u 0 ( r ) < 0 and u 0 ( r ) < 0 ,
(1.3)

they derived the following blow-up rate:

C 1 ( T t ) 2 p 1 max x Ω ¯ u(x,t) C 2 ( T t ) 2 p 1 ,tT,
(1.4)

for the special case p>1 and q=0. It is necessary to point out that assumption (1.2) seems to be reasonable, but unfortunately, the authors of [5] did not give a relationship between u 0 (x) and equation (1.2). The characterization of the monotonicity condition (1.2) was given by Souplet in [6], who proved the existence of monotone in time solutions for the above problem and obtained the blow-up rate (1.4) without the assumption of condition (1.3).

Zhou et al. [7] considered the following singular diffusion equation with memory term:

( u k ) t u= u q 0 t u p ds,(x,t)Ω×(0,T),

where p0, q0, 0<k<1. They got similar results by the method of upper-lower solution. We should notice that this kind of equation can be turned into a degenerate porous medium equation by suitable transformation. In addition, for the system of porous medium equations with nonlinear memory terms and a homogeneous Dirichlet boundary condition, one can refer for example to [8, 9].

Recently, Liu and Mu [10] considered the following semilinear parabolic equation with memory term:

u t =u+ u q 0 t u p ds,(x,t)Ω×(0,T),

subject to a weighted nonlinear nonlocal boundary,

u(x,t)= Ω f(x,y) u l (y,t)dy,(x,t)Ω×(0,T),

where p,q0. They gave the conditions of global existence and blow-up of solutions and the blow-up rate of solutions for q=0, l=1 by establishing an auxiliary function.

In view of the works mentioned above, a nonlocal parabolic equation with time-integral term does not seem to be so much investigated as nonlocal equations with space-integral terms. Already at first glance, the problem with a memory term has some difficulties in proving the existence of non-global solutions. First if t is sufficiently small, the nonlinear memory term vanishes, and then it is not clear whether the comparison principle holds in proving the existence of global small solutions. As far as we know, there are a few papers about the blow-up phenomenon for the p-Laplacian evolution equation with nonlinear memory term. Motivated by it, we consider the global existence and blow-up properties of the following p-Laplacian evolution equation with nonlinear memory term and inner absorption term:

u t div ( | u | p 2 u ) = u α 0 t u β dsk u m ,(x,t)Ω×(0,T),
(1.5)

subject to weighted linear nonlocal boundary and initial conditions

u(x,t)= Ω f(x,y)u(y,t)dy,(x,t)Ω×(0,T),
(1.6)
u(x,0)= u 0 (x),xΩ,
(1.7)

where p2, α0, β>0, m>0, k>0, and Ω R N (N1) is a bounded domain with smooth boundary. The weight function f(x,y)0 in the boundary condition is continuous, nonnegative on Ω×Ω, and Ω f(x,y)dy>0 on Ω, while the nonnegative and nontrivial initial data u 0 (x) L (Ω) W 0 1 , p (Ω) satisfies the compatibility conditions u t (x,0)=div( | u 0 ( x ) | p 2 u 0 (x))k u 0 s (x) for xΩ and u 0 (x)= Ω f(x,y) u 0 (y)dy for xΩ, which is the closed relationship for local solvability of our problem (1.5)-(1.7) (see Section 2).

The nonlocal diffusion model like equation (1.5) arises in many natural phenomena. In some sense, this kind of nonlocal problem is closer to the actual model than the local problem, such as the model of non-Newton flux through a porous medium, the model for compressible reactive gases, the model of population dynamics, and the model of biological species with human-controlled distribution (see [2, 1114] and references therein). From a physics point of view, equation (1.5) with p=2, m,α,β=1 and k=0 appears in the theory of nuclear reactor dynamics in which case the nonlocal term with time-integral is called the memory term [15]. In fact, there are some important phenomena formulated as parabolic equations which are coupled with weighted nonlocal boundary conditions in mathematical models, such as thermoelasticity theory. In this case, the solution u(x,t) describes the entropy per volume of the materia1 (see [16, 17]).

Our main goal is to find the effects of weight function on global or non-global existence of solutions for problem (1.5)-(1.7), the suitable range of nonlinear exponent, and to give the blow-up rate estimate under some suitable conditions. In addition, we treat the nonlocal nonlinearity Hölder (non-Lipschitz) cases m or n(0,1), as well as the Lipschitz cases m,n1 in this paper. We get our main results by establishing a modified comparison principle, constructing the suitable upper and lower solutions (including the self-similar lower solutions, the eigenfunction argument and the technique of ordinary differential equation and so on) and the auxiliary function. Moreover, our results extend part of or all results in [810]. The detailed results are stated as follows.

  • For arbitrary f(x,y)>0. If α+β>max{p1,m}, then the solution of problem (1.5)-(1.7) blows up in finite time for sufficiently large initial data.

  • If Ω f(x,y)1, for xΩ. If α+β>m1, then the solution of problem (1.5)-(1.7) blows up in finite time for all strictly positive initial dates with T sufficiently large.

  • If Ω f(x,y)<1, for xΩ. α+β>1, p=2, the initial value u 0 (x) satisfies conditions (H1)-(H2) (see Section 3) and u(x,t) is the blow-up solution of problem (1.2)-(1.4), then the blow-up rate is

    c ( T t ) 1 α + β 1 u(x,t)C ( T t ) 1 α + β 1 ,

where c= [ ( α + β 1 ) T ] 1 α + β 1 , C= [ ( α + β 1 ) δ ] 1 α + β 1 and δ>0 is a constant.

  • If Ω f(x,y)<1, for xΩ,

    1. (i)

      If α+βmax{p1,m}, then the solution of problem (1.5)-(1.7) exists globally for small initial data.

    2. (ii)

      If p1<min{α,β+1,m}, and α+β>m, then the solution of problem (1.5)-(1.7) exists globally for small initial data.

    3. (iii)

      If α+β>m, then the solution of problem (1.5)-(1.7) exists globally for small enough initial data.

The rest of the paper is organized as follows. In Section 2, we give the preliminaries for our research. The proofs of blow-up results and blow-up rate of solutions are given in Section 3. In Section 4, we will deduce the results of global existence.

2 Comparison principle and local existence

Since equation (1.2) is degenerate when p>2, there is no classical solution in general. Hence, it is reasonable to find a weak solution. To this end, we first give the following definition of nonnegative weak solution of problem (1.5)-(1.7).

Definition 1 If the nonnegative function u(x,t) satisfies the following conditions:

( 1 ) u C ( 0 , T ; L ( Ω ) ) L p ( 0 , T ; W 0 1 , p ( Ω ) ) , ( 1 ) u t L 2 ( 0 , T ; L 2 ( Ω ) ) , u ( x , 0 ) = u 0 ( x ) ,
(2.1)
( 2 ) Q T u t ϕ d x d t + Q T | u | p 2 u ϕ d x d t ( 2 ) = Q T ϕ ( u α 0 t u β d s k u m ) d x d t ,
(2.2)

where ϕ L 1 (0,T; W 1 , 2 (Ω)) L 2 ( Q T ) is nonnegative, ϕ(x,t) | Ω =0, and Q T =Ω×(0,T).

(3)u(x,t)= Ω f(x,y)u(y,t)dy,xΩ,
(2.3)

then u(x,t) is called the weak solution of problem (1.5)-(1.7).

If the equalities in equations (2.1)-(2.3) are replaced by ‘≤’ and ‘≥’, we can get u ̲ (x,t) and u ¯ (x,t) which are called the lower solution and upper solution of problem (1.5)-(1.7), respectively.

The following modified comparison principle plays a crucial role in our proofs, which can be obtained by establishing a suitable test function and Gronwall’s inequality.

Proposition 1 (Comparison principle)

Suppose that u ̲ (x,t) and u ¯ (x,t) are the lower and upper solutions of problem (1.2)-(1.4), respectively. If u ̲ (x,0)0, u ¯ (x,0)ε>0 and u ̲ (x,0) u ¯ (x,0), where ε is any positive constant, then u ̲ (x,t) u ¯ (x,t) in Ω ¯ ×(0,T).

Proof For xΩ, since u ̲ (x,t) and u ¯ (x,t) are the lower and upper solutions of problem (1.5)-(1.7), respectively, it follows that

Q T ( u ̲ t u ¯ t ) ϕ d x d t + Q T ( | u ̲ | p 2 u ̲ | u ¯ | p 2 u ¯ ) ϕ d x d t Q T ϕ [ u ̲ α 0 t u ̲ β d s u ¯ α 0 t u ¯ β d s k ( u ̲ m u ¯ m ) ] d x d t = Q T ϕ [ ( u ̲ α u ¯ α ) 0 t u ̲ β d s + u ¯ α 0 t ( u ̲ β u ¯ β ) d s k ( u ̲ m u ¯ m ) ] d x d t .

Choose a test function ϕ= χ [ 0 , t ] ( u ̲ u ¯ ) + =max{ u ̲ u ¯ ,0} for t>0, where χ [ 0 , t ] is a characteristic function defined on [0,t], then we have

Q T ( u ̲ u ¯ ) t ( u ̲ u ¯ ) + d x d t + Q T ( | u ̲ | p 2 u ̲ | u ¯ | p 2 u ¯ ) ( u ̲ u ¯ ) + d x d t Q T ( u ̲ u ¯ ) + [ ( u ̲ α u ¯ α ) 0 t u ̲ β d s + u ¯ α 0 t ( u ̲ β u ¯ β ) d s k ( u ̲ m u ¯ m ) ] d x d t Q T [ Φ 1 ( x , t ) 0 t u ̲ β d s + u ¯ α T Φ 2 ( x , t ) + k Φ 3 ( x , t ) ] ( u ̲ u ¯ ) + 2 d x d t ,

where

Φ 1 ( x , t ) = α 0 1 ( θ u ̲ + ( 1 θ ) u ¯ ) α 1 d θ , Φ 2 ( x , t ) = β 0 1 ( θ u ̲ + ( 1 θ ) u ¯ ) β 1 d θ , Φ 3 ( x , t ) = m 0 1 ( θ u ̲ + ( 1 θ ) u ¯ ) m 1 d θ .

By Lemma 4.10 in [18], we know ( | u ̲ | p 2 u ̲ | u ¯ | p 2 u ¯ ) ( u ̲ u ¯ ) + 0 for p2. Moreover, it follows from m1, α1, β1 that Φ i (i=1,2,3) is bounded, and if 0α<1, or β<1, we have Φ 1 ε α 1 , Φ 2 ε β 1 since u ̲ (x,0)0, u ¯ (x,0)ε>0. Furthermore, because u ̲ (x,t) and u ¯ (x,t) are bounded functions, we can get

Q T ( u ̲ u ¯ ) t ( u ̲ u ¯ ) + dxdtC Q T ( u ̲ u ¯ ) + 2 dxdt.

That is,

Ω ( u ̲ u ¯ ) + 2 dx Ω ( u ̲ ( x , 0 ) u ¯ ( x , 0 ) ) + 2 dx+C Q T ( u ̲ u ¯ ) + 2 dxdt.

Since u ̲ (x,0) u ¯ (x,0), it follows that

Ω ( u ̲ u ¯ ) + 2 dxC Q T ( u ̲ u ¯ ) + 2 dxdt.

By Gronwall’s inequality, we can deduce that ( u ̲ u ¯ ) + =0, and so u ̲ (x,t) u ¯ (x,t) in Ω×(0,T).

For xΩ, yΩ, we have

u ̲ u ¯ Ω f(x,y) ( u ̲ ( y , t ) u ¯ ( y , t ) ) dy0,

in the case of Ω ( u ̲ u ¯ ) + dy=0 in Ω. Therefore, we obtain u ̲ (x,t) u ¯ (x,t) on Ω×(0,T), and u ̲ (x,t) u ¯ (x,t) in Ω ¯ ×(0,T). □

Next, we state the theorem of local existence and uniqueness without proof.

Theorem (Local existence and uniqueness)

Suppose that p2, α0, β>0, m>0 and k>0, the nonnegative initial data u 0 L (Ω) W 0 1 , p (Ω) satisfies the compatibility conditions u t (x,0)=div( | u 0 ( x ) | p 2 u 0 (x))k u 0 s (x) for xΩ and u 0 (x)= Ω f(x,y) u 0 (y)dy for xΩ. Then there exists a constant T >0 such that problem (1.5)-(1.7) admits a nonnegative solution uC(0,T; L (Ω)) L p (0,T; W 0 1 , p (Ω)) for each T< T . Furthermore, either T = or

lim t T sup u ( x , t ) =.

Remark 1 The existence of local nonnegative solutions in time to problem (1.5)-(1.7) can be obtained by combining Theorem 1.2 in [4] with Theorem A4′ in [19]. By the comparison principle above, we can get the uniqueness of the solutions to problem (1.5)-(1.7) with α+β1, m1.

3 Blow-up solutions and blow-up rate

Comparing with the problem under a general homogeneous Dirichlet boundary condition, the existence of weight function in the boundary condition has a great influence on the global and non-global existence of solutions.

Theorem 1 Suppose that α+β>max{p1,m}, then the solution of problem (1.5)-(1.7) blows up in finite time for arbitrary f(x,y)>0 and sufficiently large initial data.

Proof In order to prove the blow-up result, we need to establish a self-similar blow-up solution. Let

u ̲ (x,t)= ω ( r ) v ξ ( r , t ) ,with 0rR,ξ>1,0<t<3T,

where ω(r)= cos 2 ( π r 2 R ), v(x,t)=4 T 2 ω δ ( t T ) 2 = v 1 v 2 , and

v 1 =2T ω δ 2 (tT), v 2 =2T+ ω δ 2 (tT).

It is obvious that u ̲ (r,t) blows up at r=0 as t approaches T. Set

δ2,T 1 2 such that v(r,t)1.

An explicit calculation yields

r x i = x i r , r x i x i = x i x j r 3 , ω x i = ω r x i = π 2 R sin π r R r x i , ω = π 2 R sin π r R , ω x i x j = ( ω r x i ) x j = π 2 2 R 2 cos π r R x i x j r 2 π 2 R sin π r R ( x i x j r 3 ) , ω = π 2 R ( π R cos π r R + n 1 r sin π r R ) , v x i = δ ( t T ) 2 ω δ 1 ω r x i , v = δ ( t T ) 2 ω δ 1 ω , v t = 2 ω δ ( t T ) , u ̲ = ω v ξ ω ξ v ξ 1 v v 2 ξ = 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 v ξ + 1 , u ̲ = ( v + δ ξ ω δ ( t T ) 2 v ξ + 1 ω ) 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 v ξ + 1 ω .

Then there exists r 0 (0,R) satisfying

ω 0 , r 0 r R , 0 ω n π 2 2 R 2 , and ω cos 2 π r 0 2 R , 0 r r 0 ,

and

u ̲ x i = ω r x i v ξ ω ξ v ξ 1 v x i v 2 ξ , 2 u ̲ x i x j = ω v ξ r x i r x j + ξ δ v ξ 1 ( t T ) 2 ω δ 1 ( ω ) 2 r x i r x j + ω v ξ r x i x j + ξ δ 2 v ξ 1 ( t T ) 2 ω δ 1 ( ω ) 2 r x i r x j + ξ δ v ξ 1 ( t T ) 2 ω δ ω r x i r x j + ξ δ 2 ( ξ + 1 ) v ξ 2 ( t T ) 2 ( t T ) 2 ω 2 δ 1 ( ω ) 2 r x i r x j + ξ δ v ξ 1 ( t T ) 2 ω δ ( ω ) r x i x j = v ξ 1 ( 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ω x i x j r 2 ) + v ξ 1 ( 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ω ( x i x j r 3 ) ) + δ ( ξ + 1 ) v ξ 2 ( t T ) 2 ω δ 1 ( ω ) 2 [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] x i x j r 2 + δ ( δ ξ 1 ) v ξ 1 ( t T ) 2 ω δ 1 ( ω ) 2 x i x j r 2 .

Therefore we have

2 u ̲ x i x j v ξ 1 ( t T ) 2 ω δ 1 ( ω ) 2 [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] ω x i x j r 2 ,

and

j = 1 N i = 1 N u ̲ x i 2 u ̲ x i x j u ̲ x j v ξ 1 ( t T ) 2 ω δ 1 ( ω ) 2 [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] ω | ω | 2 .

In view of the above, this gives

u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β d s + k u ̲ m = u ̲ t ( | u ̲ | p 2 u ̲ + ( p 2 ) | u ̲ | p 4 ( u ̲ ) T H x ( u ̲ ) u ̲ ) u ̲ α 0 t u ̲ β d s + k u ̲ m = u ̲ t ( | u ̲ | p 2 u ̲ + ( p 2 ) | u ̲ | p 4 j = 1 N i = 1 N u ̲ x i 2 u ̲ x i x j u ̲ x j ) u ̲ α 0 t u ̲ β d s + k u ̲ m 2 ξ ω δ + 1 ( t T ) v ξ + 1 ( 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 v ξ + 1 ) p 2 { 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 v ξ + 1 + ( p 2 ) 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 v ξ + 1 } ω | ω | p 2 ω α + β v ξ α ( r , t ) 0 t 1 v ξ β ( r , s ) d s + k ω m v m ξ = 2 ξ ω δ + 1 ( t T ) v ξ + 1 ( p 1 ) [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] p 1 v ( ξ + 1 ) ( p 1 ) ω | ω | p 2 ω α + β v ξ α ( r , t ) 0 t 1 v ξ β ( r , s ) d s + k ω m v m ξ .

We will discuss the problem for two cases.

Case 1. 0<t T 2 , 3 T 2 <v4 T 2 . We need to show that for sufficiently small T,

2 ξ ω δ + 1 ( t T ) v ξ + 1 +k ω m v m ξ ( p 1 ) [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] p 1 v ( ξ + 1 ) ( p 1 ) ω | ω | p 2 0.

That is,

ω v ( ξ + 1 ) ( p 2 ) m ξ [ 2 ξ ω δ + 1 ( T t ) v m ξ k ω m v ξ + 1 ] ( p 1 ) [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] p 1 | ω | p 2 2 ξ ω δ + 1 ( T t ) v ( ξ + 1 ) ( p 2 ) ( p 1 ) [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] p 1 | ω | p 2 .

Let δ be sufficiently large, satisfying ξδ>2. By the condition 3 T 2 <v4 T 2 , 0<ω1, and |ω|< π 2 R , we just have to make the following equality hold:

ω 1 p 1 2 ξ ω δ + 1 ( 3 T 2 ) ( ξ + 1 ) ( p 2 ) ( 3 ξ δ T 2 ) p 1 1 ( π 2 R ) p 2 .

It is obvious that it holds for r 0 <rR.

For 0r< r 0 , choose ξ> 1 p 2 , T<min{ 1 2 , [ n π p ( ξ δ ) p 1 2 p 1 3 ξ p 2 ξ 1 R p ( cos π r 0 2 R ) 2 ( δ + 1 ) ] 1 2 ξ p 4 ξ 1 } such that

u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β dsk u ̲ m 0.

Case 2. T 2 <t<3T, choose ξβ>2 and we can get

u ̲ α 0 t u ̲ β d s ω α + β v ξ α v 2 ξ β 0 t 1 v 1 ξ β d s = ω α + β ω δ 2 v ξ α v 2 ξ β 1 β ξ 1 ( 1 v 1 ξ β 1 1 v 1 ξ β 1 ( 0 ) ) = ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 ω α + β δ 2 v ξ α v 2 ξ β v 1 ξ β 1 ( 0 ) 1 β ξ 1 ,

and

u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β d s + k u ̲ m 4 ξ ω δ + 1 T v ξ + 1 ( p 1 ) [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] p 1 v ( ξ + 1 ) ( p 1 ) ω | ω | p 2 ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 + ω α + β δ 2 v ξ α v 2 ξ β v 1 ξ β 1 ( 0 ) 1 β ξ 1 + k ω m v m ξ .
(3.1)
  1. (i)

    r 0 rR; we have ω0, then

    u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β d s + k u ̲ m 4 ξ ω δ + 1 T v ξ + 1 ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 + ω α + β δ 2 v ξ α v 2 ξ β v 1 ξ β 1 ( 0 ) 1 β ξ 1 + k ω m v m ξ .

Since α+β>max{p1,m}, choose ξmax{2 ( α + β 1 ) 1 ,2 β 1 , ( α + β m ) 1 } and δmax{α+β1,2(α+βm),2} to satisfy ξ(α+βp+1m)>p. However,

v 1 v 1 ( 0 ) = 2 T ω δ 2 ( t T ) 2 T + ω δ 2 T 2 T ω δ 2 ( T 2 T ) 2 T + ω δ 2 T = 2 + 1 2 ω δ 2 2 + ω δ 2 <1,

thus we find

v 1 ξ β 1 v 1 ξ β 1 ( 0 ) v 1 v 1 ( 0 ) 2 + 1 2 ω δ 2 2 + ω δ 2 for ξβ1>1.

Then

u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β d s + k u ̲ m ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 [ 4 ξ T ( ξ β 1 ) v 2 v ξ ( α + β 1 ) 2 ω 3 δ 2 ( α + β 1 ) + v 1 ξ β 1 v 1 ξ β 1 ( 0 ) 1 + k ω m α β + δ 2 v 2 v ξ ( m α β ) + 1 ( ξ β 1 ) ] ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 [ 8 ξ T ( ξ β 1 ) ω 3 δ 2 ( α + β 1 ) + 2 + 1 2 ω δ 2 2 + ω δ 2 1 + 4 T k ( ξ β 1 ) ω m α β + δ 2 ] ω α + β δ 2 v ξ ( α + β ) 1 v 2 ( 2 + ω δ 2 ) 1 β ξ 1 [ 8 ξ T ( ξ β 1 ) ω 3 δ 2 ( α + β 1 ) ( 2 + ω δ 2 ) + ( 2 + 1 2 ω δ 2 ) ( 2 + ω δ 2 ) + 4 k T ( 2 + ω δ 2 ) ( ξ β 1 ) ω m α β + δ 2 ] ω α + β δ 2 v ξ ( α + β ) 1 v 2 ( 2 + ω δ 2 ) 1 β ξ 1 [ 24 ξ T ( ξ β 1 ) ω 3 δ 2 ( α + β 1 ) 1 2 ω δ 2 + 12 k T ( ξ β 1 ) ω m α β + δ 2 ]

in the sense of 0<ω1 and v 2 4T2.

In order to get the result, we have to show that

24ξT(ξβ1) ω 3 δ 2 + ( α + β 1 ) +12kT(ξβ1) ω m α β + δ 2 1 2 ω δ 2 .

Note that 0<ω1 and δmax{α+β1,2(α+βm),2}, we choose

Tmin { 1 2 , 1 24 ( 2 ξ + k ) ( ξ β 1 ) } ,

in which case we can get the result.

  1. (ii)

    0r< r 0 ; we have 0ω n π 2 2 R 2 and ω cos 2 π r 0 2 R . Since 0<ω1 and v(r,t)1, we know that 4 ξ ω δ + 1 T v ξ + 1 4 ξ T v ( ξ + 1 ) ( p 1 ) , and

    ( p 1 ) [ 4 T 2 + ( δ ξ 1 ) ω δ ( t T ) 2 ] p 1 v ( ξ + 1 ) ( p 1 ) ( ω ) | ω | p 2 ( p 1 ) ( 4 ξ δ T 2 ) p 1 v ( ξ + 1 ) ( p 1 ) n π 2 2 R 2 ( π 2 R ) p 2 = ( p 1 ) ( 2 ξ δ T 2 ) p 1 v ( ξ + 1 ) ( p 1 ) n π p R p .
    (3.2)

Let A= ( cos 2 π r 0 2 R ) δ , B= 2 + A 2 2 + A , it is obvious that 1 2 <B<1 and ω δ A. Substituting equation (3.2) into equation (3.1) gives

u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β d s + k u ̲ m 4 ξ T v ( ξ + 1 ) ( p 1 ) + ( p 1 ) ( 2 ξ δ T 2 ) p 1 v ( ξ + 1 ) ( p 1 ) n π p R p + ω α + β δ 2 v ξ α v 2 ξ β v 1 ξ β 1 ( 0 ) 1 β ξ 1 ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 + k ω m v m ξ 4 ξ T + ( p 1 ) ( 2 ξ δ T 2 ) p 1 n π p v ( ξ + 1 ) ( p 1 ) R p + ω α + β δ 2 v ξ α v 2 ξ β v 1 ξ β 1 ( 0 ) 1 β ξ 1 ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 + k v m ξ = 4 ξ T + ( p 1 ) ( 2 ξ δ T 2 ) p 1 n π p + k R p v ( ξ + 1 ) ( p 1 ) + m ξ R p + ω α + β δ 2 v ξ α v 2 ξ β 1 β ξ 1 ( 1 v 1 ξ β 1 ( 0 ) B ξ δ 1 v 1 ξ δ 1 ) ( 1 B ξ δ 1 ) ω α + β δ 2 v ξ ( α + β ) 1 v 2 1 β ξ 1 4 ξ T + ( p 1 ) ( 2 ξ δ T 2 ) p 1 n π p + k R p v ( ξ + 1 ) ( p 1 ) + m ξ R p ( 1 B ξ δ 1 ) A α + β δ 1 2 v ξ ( α + β ) 1 1 β ξ 1 + ω α + β δ 2 v ξ α v 2 ξ β 1 β ξ 1 ( 1 v 1 ξ β 1 ( 0 ) B ξ δ 1 v 1 ξ δ 1 ) .

However, since v 1 =2T ω δ 2 (tT)2T+ T 2 ω δ 2 and 1 2 <B<1, it follows that

1 v 1 ξ β 1 ( 0 ) B ξ δ 1 v 1 ξ δ 1 .

Finally, we need to show that

4 ξ T + ( p 1 ) ( 2 ξ δ T 2 ) p 1 n π p + k R p v ( ξ + 1 ) ( p 1 ) + m ξ R p ( 1 B ξ δ 1 ) A α + β δ 1 2 v ξ ( α + β ) 1 1 β ξ 1 .

Because T 1 2 and v4 T 2 1, we have to show

[ 2 ξ + ( p 1 ) ( 2 ξ δ ) p 1 n π p + k R p R p ] ( 4 T 2 ) ξ ( α + β p + 1 m ) p ( 1 B ξ δ 1 ) A α + β δ 1 2 β ξ 1 .

In other words, we just need the following inequality:

4 T 2 [ ( 1 B ξ δ 1 ) A α + β δ 1 2 R p ( ξ β 1 ) ( 4 ξ + ( 2 ξ δ ) p 1 ) n π p ] 1 ξ ( α + β p + 1 m ) p

to hold. So choose T to be small enough, and we can get

u ̲ t div ( | u ̲ | p 2 u ̲ ) u ̲ α 0 t u ̲ β dsk u ̲ m 0.

For xΩ, ω(R)=0,

u ̲ (R,t)= ω ( R ) v ξ ( R , t ) =0 Ω f(x,y) u ̲ dy,

choose u ̲ (r,0)= ω ( r ) v ξ ( r , 0 ) = ω ( r ) 4 T 2 ω δ T 2 u 0 , then u ̲ (r,t) is the lower solution of problem (1.5)-(1.7). This implies that the solution blows up in finite time for large enough initial data. □

Theorem 2 Suppose that Ω f(x,y)dy1 for xΩ. If α+β>m1, then the solution of problem (1.5)-(1.7) blows up in finite time for all strictly positive initial dates with T sufficiently large.

Proof Consider the following problem:

v (t)=T v α + β k v m ,v(0)= v 0 .
(3.3)

As α+β>m, we know v α + β +1> v m , and T v α + β k v m (Tk) v α + β k. Therefore, the solution of equation (3.3) is an upper solution of the following problem:

v (t)=(Tk) v α + β k,v(0)= v 0 .

Since T>k and α+β>1, the solution of this problem blows up in finite time if v 0 > ( k T k ) 1 α + β .

It is obvious that the solution of problem (3.3) is a lower solution of problem (1.5)-(1.7) when Ω φ(x,y)dy1 and u 0 (x)> v 0 . By Proposition 1, u(x,t) is a blow-up solution. □

Suppose that the solution u(x,t) of problem (1.5)-(1.7) with p=2 blows up in finite time, and let U(t)= max x Ω ¯ u(x,t). We suppose that the initial data satisfies the following assumptions:

(H1) u 0 (x)+ u 0 α 0 t u 0 β dsk u 0 m >0.

(H2) There exists a constant δ>0 such that u 0 (x)+ u 0 α 0 t u 0 β dsk u 0 m δ u 0 α + β 0.

Theorem 3 Suppose that Ω f(x,y)dy1 for xΩ. If α+β>1, u 0 (x) satisfies condition (H1)-(H2), and u(x,t) is the blow-up solution of problem (1.5)-(1.7) in finite time T with p=2, then the blow-up rate is

c ( T t ) 1 α + β 1 u(x,t)C ( T t ) 1 α + β 1 ,

where c= [ ( α + β 1 ) T ] 1 α + β 1 and C= [ δ ( α + β 1 ) ] 1 α + β 1 .

Remark 2 Choose Ω=(1,2), α=2, β=1, m=2 and f(x,y)= 1 8 xy, one can easily verify that u 0 (x)= x 2 satisfies (C1)-(C2), conditions in Theorem 3 are thus valid.

Lemma 1 If u 0 (x) satisfies condition (H1)-(H2), p=2, then there exists a positive constant c= [ ( α + β 1 ) T ] 1 α + β 1 such that U(t)c ( T t ) 1 α + β 1 .

Proof It is obvious that U(t) is Lipschitz continuous and differentiable almost everywhere. By equation (1.5) with p=2 and U(t)0, it yields

U (t) U α 0 t U β dsk U m T U α + β ,

and thus

( U 1 α β ( t ) ) T(α+β1).

Integrating the result above over (t,T), we can obtain the conclusion. □

Proof of Theorem 3 Let J= u t δ u α + β , where δ>0 is sufficiently small. Since α+β>1, we have

J t = u t t δ ( α + β ) u α + β 1 u t = u t + α u α 1 u t 0 t u β d s + u α + β k m u m 1 u t δ ( α + β ) u α + β 1 u t , J = u t δ ( α + β ) u α + β 1 u t + δ ( α + β ) u 2 α + β 1 × 0 t u β d s δ ( α + β ) ( α + β 1 ) u α + β 2 | u | 2 , J t J α u α 1 u t 0 t u β d s + u α + β δ ( α + β ) u 2 α + β 1 0 t u β d s k m u m 1 u t = ( α u α 1 0 t u β d s k m u m 1 ) J + u α + β ( 1 δ β u α 1 0 t u β d s + k m u m 1 δ ) ,

so we can choose δ>0 to be small enough and thus obtain

J t Jα u α 1 J 0 t u β ds.

On the other hand, as (x,t)Ω×(0,T), we get

J = u t δ u α + β = Ω f ( x , y ) u t d y δ ( Ω f ( x , y ) u d y ) α + β = Ω f ( x , y ) J d y + δ Ω f ( x , y ) u α + β d y δ ( Ω f ( x , y ) u d y ) α + β .

Let 0<F(x)= Ω f(x,y)dy1, By Jensen’s inequality, this gives

J Ω f ( x , y ) J d y + δ F ( x ) [ Ω f ( x , y ) u d y F ( x ) ] α + β δ ( Ω f ( x , y ) u d y ) α + β Ω f ( x , y ) J d y ,

and

J(x,t)= u t δ u α + β =Δu+ u α 0 t u β dsδ u α + β .

It follows from the assumptions of (H1)-(H2) that J(x,0)0. Therefore, it is easy to deduce that J(x,t)0 for (x,t) Ω ¯ ×[0,T). That is, u t δ u α + β and integrating this over (t,T) yields uC ( T t ) 1 α + β 1 . Combining the results with Lemma 1, we obtain the desired result. □

4 Global existence of solutions

In this section, we give sufficient conditions of the global existence of solutions.

Theorem 4 Suppose that Ω f(x,y)dy<1 for xΩ. If α+βmax{p1,m}, then the solution of problem (1.5)-(1.7) exists globally for small initial data.

Proof Let u ¯ (x,t)=q e n t ψ(x), where q,n>0 are determined later, and ψ(x) solves the following problem:

div | ψ | p 2 ψ = 1 , x Ω , ψ ( x ) = Ω f ( x , y ) d y , x Ω ,

where 0<ψ(x)<1. Let η 1 = max Ω ¯ ψ(x), η 2 = max Ω ¯ ψ(x),

u ¯ t div ( | u ¯ | p 2 u ¯ ) u ¯ α 0 t u ¯ β d s + k u ¯ m = q n e n t ψ + q p 1 e n ( p 1 ) t 1 n β q α + β e n ( α + β ) t ψ α + β + k q m e m n t ψ m q p 1 e n ( p 1 ) t 1 n β q α + β e n ( α + β ) t ψ α + β + k q m e m n t ψ m ;

since α+β<max{p1,m}, choosing

n>max { 1 k β ( q η 1 ) α + β m , 1 β q α + β p + 1 η 2 α + β } ,

we can infer that

u ¯ t div ( | u ¯ | p 2 u ¯ ) u ¯ α 0 t u ¯ β ds+k u ¯ m 0.

For xΩ, 0<ψ(x)<1

u ¯ (x,t)=q e n t Ω f(x,y)dy Ω f(x,y) u ¯ dy.

Selecting u 0 (x)<q η 1 , we can deduce that the result holds. □

Theorem 5 Suppose that Ω f(x,y)dyρ<1 for xΩ. If p1<min{α,β+1,m}, and α+β>m, then the solution of problem (1.5)-(1.7) exists globally for small initial data.

Proof Suppose that η(x) solves the following problem:

div | η | p 2 η = 1 , x Ω , η ( x ) = 0 , x Ω ,

in which 0<η(x)<1. Let C= max Ω ¯ η(x).

Set u ¯ (x,t)= k 1 ( A + t ) r ( ρ 1 ρ + η ( x ) C ), where A, k 1 1 and 0r 1 p 2 , then

u ¯ t = k 1 r ( A + t ) r 1 ( ρ 1 ρ + η ( x ) C ) k 1 r ( A + t ) r 1 1 1 ρ k 1 r ( A + t ) r ( p 1 ) 1 1 ρ ,

and

u ¯ α 0 t u ¯ β d s = k 1 α + β ( A + t ) r α ( η ( x ) C + ρ 1 ρ ) α + β 0 t ( A + t ) r β d s k 1 α + β 1 r β ( A + t ) 1 r α r β ( η ( x ) C + ρ 1 ρ ) α + β .

For p1<min{α,β+1,m}, we choose 0r 1 p 2 such that 1r(α+βm)<0 and rβ1. It follows that

u ¯ α 0 t u ¯ β ds k 1 α + β 1 r β ( A + t ) 1 r α r β ( 1 1 ρ ) α + β ,

and

u ¯ t div ( | u ¯ | p 2 u ¯ ) u ¯ α 0 t u ¯ β d s + k u ¯ m k 1 r ( A + t ) r ( p 1 ) 1 1 ρ + k 1 p 1 ( A + t ) r ( p 1 ) k 1 α + β 1 r β ( A + t ) 1 r α r β ( 1 1 ρ ) α + β + k k 1 m ( A + t ) m r ( η ( x ) C + ρ 1 ρ ) m = k 1 ( A + t ) r ( p 1 ) ( r 1 ρ + k 1 p 2 2 ) + k 1 p 1 ( A + t ) r ( p 1 ) × [ 1 2 k 1 α + β p + 1 1 r β ( A + t ) 1 r ( α + β p + 1 ) ( 1 1 ρ ) α + β + k k 1 m p + 1 ( A + t ) r ( m p + 1 ) ( η ( x ) C + ρ 1 ρ ) m ] .

Since A, k 1 1, choosing

k 1 max { ( 2 r 1 ρ ) p 2 , 1 }

and

Amax { [ k ( 1 r β ) ( 1 ρ k 1 ) α + β m ] 1 1 r ( α + β m ) , 1 } ,

we can get

u ¯ t div ( | u ¯ | p 2 u ¯ ) u ¯ α 0 t u ¯ β ds+k u ¯ m 0.

On the other hand, for xΩ and sufficiently large A, we have

u ¯ ( x , t ) = k ( A + t ) r ρ 1 ρ k ( A + t ) r 1 1 ρ Ω f ( x , y ) d y Ω k ( A + t ) r ( η ( x ) C + ρ 1 ρ ) f ( x , y ) d y = Ω u ¯ f ( x , y ) d y .

Choosing u 0 to be sufficiently small such that u 0 k A r ( η ( x ) C + ρ 1 ρ ), we can conclude that u ¯ is an upper solution of problem (1.5)-(1.7). The proof is completed. □

Theorem 6 Suppose that Ω f(x,y)dy<1 for xΩ. If α+β>m, then the solution of problem (1.5)-(1.7) exists globally for small enough initial data.

Proof Choosing u ¯ (x,t)=A and A>max{ ( k T ) 1 α + β 1 , u 0 (x)}, it is easy to see that the result holds. □

References

  1. Friedman A: Monotonic decay of solutions of parabolic equations with nonlocal boundary conditions. Q. Appl. Math. 1986, 44(3):401-407.

    MATH  Google Scholar 

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

    MATH  Google Scholar 

  3. Akila Y: On a nonlocal parabolic problem. Demonstr. Math. 2009, 42(4):745-755.

    MATH  MathSciNet  Google Scholar 

  4. Wang YL, Mu CL, Xiang ZY: Blow-up of solutions to a porous medium equation with nonlocal boundary condition. Appl. Math. Comput. 2007, 192: 579-585. 10.1016/j.amc.2007.03.036

    Article  MATH  MathSciNet  Google Scholar 

  5. Li YX, Xie CH: Blow-up for semi-linear parabolic equations with nonlinear memory. Z. Angew. Math. Phys. 2004, 55: 15-27. 10.1007/s00033-003-1128-6

    Article  MATH  MathSciNet  Google Scholar 

  6. Souplet P: Monotonicity of solutions and blow-up for semilinear parabolic equations with nonlinear memory. Z. Angew. Math. Phys. 2004, 55: 28-31. 10.1007/s00033-003-1158-0

    Article  MATH  MathSciNet  Google Scholar 

  7. Zhou J, Mu CL, Lu F: Blow-up and global existence to a degenerate reaction-diffusion equation with nonlinear memory. J. Math. Anal. Appl. 2007, 333: 1138-1152. 10.1016/j.jmaa.2006.12.007

    Article  MATH  MathSciNet  Google Scholar 

  8. Du LL, Mu CL: Global existence and blow-up analysis to a degenerate reaction-diffusion system with nonlinear memory. Nonlinear Anal., Real World Appl. 2008, 9: 303-315. 10.1016/j.nonrwa.2006.10.005

    Article  MATH  MathSciNet  Google Scholar 

  9. Zhou J, Mu CL, Fan MS: Global existence and blow-up to a degenerate reaction-diffusion system with nonlinear memory. Nonlinear Anal., Real World Appl. 2008, 9: 1518-1534. 10.1016/j.nonrwa.2007.03.016

    Article  MATH  MathSciNet  Google Scholar 

  10. Liu DM, Mu CL: Blow-up analysis for a semi-linear parabolic equation with nonlinear memory and nonlocal nonlinear boundary condition. Electron. J. Qual. Theory Differ. Equ. 2010, 51: 1-17.

    MathSciNet  Google Scholar 

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

    Book  MATH  Google Scholar 

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

    Article  MATH  MathSciNet  Google Scholar 

  13. Calsina A, Perello C, Saldana J: Non-local reaction-diffusion equations modelling predator-prey co-evolution. Publ. Mat. 1994, 32: 315-325.

    Article  MathSciNet  Google Scholar 

  14. Allegretto W, Fragnelli G, Nistri P, Papin D: Coexistence and optimal control problems for a degenerate predator-prey model. J. Math. Anal. Appl. 2011, 378: 528-540. 10.1016/j.jmaa.2010.12.036

    Article  MATH  MathSciNet  Google Scholar 

  15. Kastenberg WE: Space dependent reactor kinetics with positive feed-back. Nukleonika 1968, 11: 126-130.

    Google Scholar 

  16. Pao CV: Asymptotic behavior of solutions of reaction-diffusion equations with nonlocal boundary conditions. J. Comput. Appl. Math. 1998, 88(1):225-238. 10.1016/S0377-0427(97)00215-X

    Article  MATH  MathSciNet  Google Scholar 

  17. Pao CV: Numerical solutions of reaction-diffusion equations with nonlocal boundary conditions. J. Comput. Appl. Math. 2001, 136(1-2):227-243. 10.1016/S0377-0427(00)00614-2

    Article  MATH  MathSciNet  Google Scholar 

  18. Diaz JI: Nonlinear Partial Differential Equations and Free Boundaries: Elliptic Equations. Vol. I. Pitman, London; 1985.

    Google Scholar 

  19. Soupllet P: Blow-up in nonlocal reaction-diffusion equations. SIAM J. Math. Anal. 1998, 29(6):1301-1334. 10.1137/S0036141097318900

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

This work is supported by the Natural Science Foundation of Shandong Province of China (ZR2012AM018) and the Fundamental Research Funds for the Central Universities (No. 201362032). The authors would like to warmly thank all the 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

Fang, Z.B., Zhang, J. Global existence and blow-up of solutions for p-Laplacian evolution equation with nonlinear memory term and nonlocal boundary condition. Bound Value Probl 2014, 8 (2014). https://doi.org/10.1186/1687-2770-2014-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-8

Keywords