Skip to main content

Large time behavior of solutions for the porous medium equation with a nonlinear gradient source

Abstract

This paper deals with the large time behavior of non-negative solutions for the porous medium equation with a nonlinear gradient source u t =Δ u m + | ∇ u l | q , (x,t)∈Ω×(0,∞), where l≥m>1 and 1≤q<2. When lq=m, we prove that the global solution converges to the separate variable solution t − 1 m − 1 f(x). While m<lq≤m+1, we note that the global solution converges to the separate variable solution t − 1 m − 1 f 0 (x). Moreover, when lq>m+1, we show that the global solution also converges to the separate variable solution t − 1 m − 1 f 0 (x) for the small initial data u 0 (x), and we find that the solution u(x,t) blows up in finite time for the large initial data u 0 (x).

MSC:35K55, 35K65, 35B40.

1 Introduction

In this paper, we investigate the large time behavior of non-negative solutions for the following initial-boundary value problem:

{ u t = Δ u m + | ∇ u l | q , ( x , t ) ∈ Ω × ( 0 , ∞ ) , u ( x , t ) = 0 , ( x , t ) ∈ ∂ Ω × ( 0 , ∞ ) , u ( x , 0 ) = u 0 ( x ) , x ∈ Ω ,
(1.1)

where l≥m>1, 1≤q<2, Ω is a bounded domain of R N (N≥2) with smooth boundary ∂ Ω, and the initial function is

u 0 (x)∈ C 0 ( Ω ¯ )= { z ∈ C ( Ω ¯ ) : z = 0  on  ∂ Ω } , u 0 (x)≥0, u 0 (x)≢0.
(1.2)

Equation (1.1) arises in the study of the growth of surfaces and it has been considered as a mathematical model for a variety of physical problems (see [1, 2]). For instance, in ballistic deposition processes, the evolution of the profile of a growing interface is described by the diffusive Hamilton-Jacobi type equation (1.1) with m=l=1 (see [3]).

One of the particular feature of problem (1.1) is that the equation is a slow diffusion equation with nonlinear source term depending on the gradient of a power of the solution. In general, there is no classical solution. Therefore, it turns out that a suitable framework for the well-posedness of the initial-boundary value problem (1.1) is the theory of viscosity solutions (see [4–6]), so we first define the notion of solutions.

Definition 1.1 A non-negative function u(x,t)∈C( Ω ¯ ×(0,∞)) is called a solution of (1.1), if u(x,t) is a viscosity solution to (1.1) in Ω×(0,∞) and satisfies

u(x,t)=0,(x,t)∈∂Ω×(0,∞)andu(x,0)= u 0 (x)∈ C 0 ( Ω ¯ ),x∈ Ω ¯ .
(1.3)

Under some assumptions, the global (local) existence in time, uniqueness and regularity of solutions to reaction-diffusion equations with gradient terms have been extensively investigated by many authors (see [7–12] and the references therein). In particular, in [1], Andreucci proved the existence of solutions for the following degenerate parabolic equation with initial data measures:

{ u t = Δ u m + | ∇ u l | q , ( x , t ) ∈ R N × ( 0 , ∞ ) , u ( x , 0 ) = μ , x ∈ R N ,
(1.4)

where m≥1, ql≥m, 0<q<2, and N≥1.

The main purpose of this paper is to further study the large time behavior of non-negative solutions u(x,t) to (1.4) with homogeneous Dirichlet boundary conditions. In recent years, many authors have investigated the asymptotic behavior of solutions to the viscous Hamilton-Jacobi equations (see [3, 4, 7, 8, 10, 13–21] and the references therein). For example, for the special case m=l=1, Gilding [22] studied the large time behavior of solutions to the following Cauchy problem:

{ u t = Δ u + | ∇ u | q , ( x , t ) ∈ R N × ( 0 , ∞ ) , u ( x , 0 ) = u 0 ( x ) , x ∈ R N ,
(1.5)

and he gave the temporal decay estimates.

In [23], Stinner investigated the asymptotic behavior of solutions for the following one space dimensional viscous Hamilton-Jacobi equation:

{ u t = ( | u x | p − 2 u x ) x + | u x | q , ( x , t ) ∈ ( − R , R ) × ( 0 , ∞ ) , u ( ± R , t ) = 0 , t ∈ ( 0 , ∞ ) , u ( x , 0 ) = u 0 ( x ) , x ∈ ( − R , R ) ,
(1.6)

where R>0, p>2, and 1<q<p−1, and he proved that these solutions converge to the steady states by Lyapunov functional.

In higher dimensional case, Barles et al. [24] studied the large time behavior of solutions for the following initial-boundary value problem:

{ u t = div ( | ∇ u | p − 2 ∇ u ) + | ∇ u | q , ( x , t ) ∈ B ( 0 , 1 ) × ( 0 , ∞ ) , u ( x , t ) = 0 , ( x , t ) ∈ ∂ B ( 0 , 1 ) × ( 0 , ∞ ) , u ( x , 0 ) = u 0 ( x ) , x ∈ B ( 0 , 1 ) ,
(1.7)

where p≥2, 0<q<p−1, and they showed that the non-negative radially symmetric solutions converge to the stationary solution.

Recently, in [25], Laurencot et al. extended the case 0<q<p−1 of the problem (1.7) to the case q≥p−1, and derived that these solutions converge to two different separate variables solutions according to the cases q=p−1 and q>p−1 in the general bounded domain Ω⊂ R N , respectively.

Motivated by the above work, by using the modified comparison argument, the self-similar transformation method, and the half-relaxed limits technique used in [6, 25, 26], we investigate the asymptotic behavior of non-negative solutions to (1.1). Our main results in this paper are stated as follows.

Theorem 1.1 Let l≥m>1, 1≤q<2, and lq=m. Assume that u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2). Then there exists a unique solution u(x,t) to (1.1) in the sense of Definition  1.1 such that

lim t → ∞ ∥ t 1 m − 1 u ( t ) − f ∥ ∞ =0,
(1.8)

where f(x)∈ C 0 ( Ω ¯ ) is the unique positive solution to

−Δ f m − | ∇ f l | q − f m − 1 =0in Î©,f>0in Î©,f=0on âˆ‚Ω.
(1.9)

Moreover, we have |∇u(t)|∈ L ∞ ( Ω ¯ ) for all t≥0 and

ℓ[ u 0 ]:= sup t ≥ 0 { ∥ ∇ u ( t ) ∥ ∞ } <∞.
(1.10)

Theorem 1.2 Let l≥m>1, 1≤q<2, and m<lq≤m+1. Assume that u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2). Then there exists a unique solution u(x,t) to (1.1) in the sense of Definition  1.1 such that

lim t → ∞ ∥ t 1 m − 1 u ( t ) − f 0 ∥ ∞ =0,
(1.11)

where f 0 (x)∈ C 0 ( Ω ¯ ) is the unique positive solution to

−Δ f 0 m − f 0 m − 1 =0in Î©, f 0 >0in Î©, f 0 =0on âˆ‚Ω.
(1.12)

Theorem 1.3 Let l≥m>1, 1≤q<2, and lq>m+1. Assume that u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2) and suppose further that there exists G 0 ∈ W 1 , ∞ (Ω) satisfying (1.2) such that

u 0 (x)≤ G 0 ( x ) ℓ [ G 0 ] ,x∈ Ω ¯ ,
(1.13)

where â„“[ G 0 ] is defined in (1.10). Then there exists a unique solution u(x,t) to (1.1) in the sense of Definition  1.1 such that

lim t → ∞ ∥ t 1 m − 1 u ( t ) − f 0 ∥ ∞ =0,
(1.14)

where f 0 (x)∈ C 0 ( Ω ¯ ) is defined in (1.12).

Theorem 1.4 Let l≥m>1, 1≤q<2, and lq>m+1. Assume that there exists a positive constant K 0 depending only on m, l, q, s, and ε such that ∫ Ω u 0 s (x) ϕ ε (x)dx≥ K 0 , where s∈(0, 1 2 ), ε∈(0,1), and ϕ ε (x)= A ε e − ε | x | with A ε = 1 ∫ Ω e − ε | x | d x . Then the solution u(x,t) of the problem (1.1) blows up in finite time in the sense of weak solution. Moreover, the upper bound of blow-up time is given as follows:

T≤ 2 q l − 1 ( q l q l + s − 1 ) q ( m − 1 ) q l − m ( m ε 2 ) 1 − q l q l − m .
(1.15)

Remark 1.1 Compared to the results in [25], we extend the results of p-Laplacian equation to the porous medium equation (1.1) with a nonlinear gradient source.

Remark 1.2 In Theorem 1.4, we only give an upper estimate of the blow-up time. But the lower estimate of the blow-up time is an open problem.

This paper is organized as follows. In Section 2, we establish the comparison lemmas to prove the uniqueness of the positive solution to (1.9) and the identification of the half-relaxed limits. In Section 3, using the comparison principle, we construct the global solutions to obtain the upper bound and Hölder estimate of solutions to (1.1), and we prove Theorems 1.1 and 1.2 by the half-relaxed limits method. Moreover, we give the large time behavior of solutions to (1.1) with the small initial data u 0 (x) for lq>m+1, and we prove Theorem 1.3 in Section 4. Finally, we obtain the blow-up case, and we prove Theorem 1.4 in Section 5.

2 Comparison lemmas

In this section, we establish the following comparison lemma between positive supersolutions and non-negative subsolutions to the elliptic equation in (1.9): which is an important tool for the uniqueness of the positive solution to (1.9) and the identification of the half-relaxed limits later.

Lemma 2.1 Let l≥m>1, 1≤q<2, and lq=m. Assume that w∈USC( Ω ¯ ) and W∈LSC( Ω ¯ ) are respectively a bounded upper semicontinuous (usc) viscosity subsolution and a bounded lower semicontinuous (lsc) viscosity supersolution to

−Δ ζ m − | ∇ ζ l | q − ζ m − 1 =0in Î©,
(2.1)

such that

w(x)=W(x)=0for x∈∂Ω,
(2.2)

and

W(x)>0for x∈Ω.
(2.3)

Then

w(x)≤W(x)for x∈ Ω ¯ .
(2.4)

Proof The proof is based on the idea as in [25, 27], but with different auxiliary functions.

For i≥ N 0 large enough, it is easy to see that Ω i ={x∈Ω:d(x,∂Ω)> 1 i } is a non-empty open subset of R N .

Since Ω ¯ i is compact and W is lower semicontinuous, the function W has a minimum in  Ω ¯ i . By the positivity (2.3) of W in Ω ¯ i , we have

μ i = min x ∈ Ω ¯ i { W ( x ) } >0.
(2.5)

Similarly, the compactness of Ω ¯ ∖ Ω i , the upper semicontinuity and boundedness of w imply that w has at least one point of maximum x i in Ω ¯ ∖ Ω i and we set

η i =w( x i )= max x ∈ Ω ¯ ∖ Ω i { w ( x ) } .
(2.6)

It follows from ∂Ω⊂ Ω ¯ ∖ Ω i and w vanishes on ∂ Ω by (2.2) that η i ≥0.

Next, we claim that

lim i → ∞ η i =0.
(2.7)

Indeed, owing to the compactness of Ω ¯ and the definition of Ω i , there exist y∈∂Ω and a subsequence of { x i } i ≥ N 0 (not relabeled) such that x i →y as i→∞. Since w(y)=0, we deduce from the upper semicontinuity of w that

lim ε → 0 sup { w ( x ) : x ∈ B ( y , ε ) ∩ Ω ¯ } ≤0.

Given ε>0 small enough, there exists i ε large enough such that

x i ∈B(y,ε)∩ Ω ¯ for i≥ i ε .

Therefore, we have

0≤ η i =w( x i )≤sup { w ( x ) : x ∈ B ( y , ε ) ∩ Ω ¯ } ,i≥ i ε .

Thus

0≤ lim i → ∞ sup η i ≤sup { w ( x ) : x ∈ B ( y , ε ) ∩ Ω ¯ } .

Letting ε→0, we conclude that zero is a cluster point of { η i } i ≥ N 0 as i→∞. The claim (2.7) follows from the monotonicity of { η i } i ≥ N 0 .

Now, fix s∈(0,∞). For δ>0 and i≥ N 0 , we define

z i (t,x)= ( t + s ) − 1 m − 1 w(x)− s − 1 m − 1 η i ,(t,x)∈[0,∞)× Ω ¯

and

Z δ (t,x)= ( t + δ ) − 1 m − 1 W(x),(t,x)∈[0,∞)× Ω ¯ .

It follows from the assumptions on w and W that z i and Z δ are, respectively, a bounded usc viscosity subsolution and a bounded lsc viscosity supersolution to

∂ t ζ−Δ ζ m − | ∇ ζ l | q =0in (0,∞)×Ω,

and satisfy

Z δ (t,x)=0≥− s − 1 m − 1 η i = z i (t,x),(t,x)∈(0,∞)×∂Ω.

Moreover, if

0<δ< [ μ i 1 + ∥ w ∥ ∞ ] m − 1 s,
(2.8)

then it follows from (2.5) and (2.8) that, for x∈ Ω i ,

Z δ (0,x)= δ − 1 m − 1 W(x)≥ δ − 1 m − 1 μ i ≥ s − 1 m − 1 ∥ w ∥ ∞ ≥ z i (0,x).

For x∈ Ω ¯ ∖ Ω i , we deduce from (2.6) that

Z δ (0,x)≥0≥ s − 1 m − 1 ( w ( x ) − η i ) = z i (0,x).

By the comparison principle [6], we have

z i (t,x)≤ Z δ (t,x),(t,x)∈[0,∞)× Ω ¯ ,
(2.9)

for i≥ N 0 and δ>0 satisfying (2.8).

According to (2.8), the parameter δ can be taken to be arbitrarily small in (2.9). Therefore, we deduce that

( t + s ) − 1 m − 1 w(x)− s − 1 m − 1 η i ≤ t − 1 m − 1 W(x),(t,x)∈(0,∞)× Ω ¯ ,

for i≥ N 0 .

Passing to the limit as i→∞, it follows from (2.7) that

( t + s ) − 1 m − 1 w(x)≤ t − 1 m − 1 W(x),(t,x)∈(0,∞)× Ω ¯ .

Finally, let s→0 and take t=1; then we obtain

w(x)≤W(x),(t,x)∈(0,∞)× Ω ¯ .

The proof of Lemma 2.1 is complete. □

A straightforward consequence of Lemma 2.1 is the uniqueness of the solution to (1.9).

Corollary 2.1 There is at most one positive viscosity solution to (1.9).

By the similar argument, we have the following result to (1.12).

Lemma 2.2 Let w∈USC( Ω ¯ ) and W∈LSC( Ω ¯ ) be respectively a bounded upper semicontinuous (usc) viscosity subsolution and a bounded lower semicontinuous (lsc) viscosity supersolution to

−Δ ζ m − ζ m − 1 =0in Î©,
(2.10)

satisfying (2.2) and (2.3). Then

w(x)≤W(x)for x∈ Ω ¯ .
(2.11)

Proof The proof is similar as in Lemma 2.1, so we omit it here. □

3 Proofs of Theorems 1.1 and 1.2

In this section, we obtain the well-posedness and large time behavior of solutions to (1.1) for lq∈[m,m+1], and we prove Theorems 1.1-1.2. To do this, we first obtain the well-posedness to (1.1) by the following proposition.

Proposition 3.1 Assume that l≥m>1, 1≤q<2, lq∈[m,m+1], and u 0 (x)∈ C 0 ( Ω ¯ ) satisfies (1.2). Then there exists a unique solution u(x,t) to (1.1) in the sense of Definition  1.1.

Proof The idea of the proof is same as in [25, 28], so we omit here. □

Next, in order to prove the large time behavior of the solution to (1.1), we shall need several steps: Step 1, we will find that the temporal decay rate of ∥ u ( t ) ∥ ∞ is indeed t − 1 m − 1 . Step 2, we prove the boundary estimates for the large time which guarantee that no loss of boundary condition occurs throughout the time evolution. Step 3, the half-relaxed limits technique is applied to show the expected convergence after introducing self-similar variables. The approach is developed by Laurençot and Stinner in [25, 27]. To do this, we need the following lemmas.

Lemma 3.1 (Upper bound)

Assume that l≥m>1, 1≤q<2, lq=m, and the initial data u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2). Then there exists C 1 >0 depending only on m, l, q, Ω, and ∥ u 0 ∥ ∞ such that

u(t,x)≤ C 1 ( 1 + t ) − 1 m − 1 ,(t,x)∈(0,∞)× Ω ¯ .
(3.1)

Proof Assume that x 0 ∉ Ω ¯ and R 0 >0 such that Ω⊂B( x 0 , R 0 ). For (t,x)∈(0,∞)× Ω ¯ , we define the function

S 1 (t,x)= A 1 ( 1 + t ) − 1 m − 1 φ(r),φ(r)= [ m m + 1 ( e R ( m + 1 ) m − e r ( m + 1 ) m ) ] 1 m with r=|x− x 0 |,

where R> R 0 and A 1 >0 satisfies the following condition:

A 1 ≥max { ( m m − 1 e R ( m + 1 ) m 2 ) 1 m − 1 , ∥ u 0 ∥ ∞ φ ( R 0 ) } .
(3.2)

Since x 0 ∉ Ω ¯ , the function S 1 (t,x) is C ∞ -smooth in [0,∞)× Ω ¯ . Moreover, according the condition l≥m>1, 1≤q<2, and lq=m, we have l=m>1 and q=1. Therefore, for (t,x)∈(0,∞)×Ω and r=|x− x 0 |< R 0 <R, it follows from (3.2) that

( 1 + t ) m m − 1 { ∂ t S 1 − △ S 1 m − | ∇ S 1 l | q } ( t , x ) = − A 1 m − 1 φ ( r ) − A 1 m [ N − 1 r ( φ m ) r + ( φ m ) r r ] − A 1 m | ( φ l ) r | q = − A 1 m − 1 φ ( r ) + A 1 m [ N − 1 r e r ( m + 1 ) m + m + 1 m e r ( m + 1 ) m ] − A 1 m e r ( m + 1 ) m ≥ A 1 [ A 1 m − 1 ( N − 1 r + 1 m ) e r ( m + 1 ) m − 1 m − 1 e R ( m + 1 ) m 2 ] ≥ A 1 ( 1 m A 1 m − 1 − 1 m − 1 e R ( m + 1 ) m 2 ) ≥ 0 .
(3.3)

Hence, the condition (3.2) guarantees that S 1 (t,x) is a supersolution to (1.1) in (0,∞)×Ω. In addition, since r=|x− x 0 |< R 0 <R for x∈Ω, for (t,x)∈(0,∞)×∂Ω, we deduce from (3.2) that

u(t,x)=0≤ A 1 ( 1 + t ) − 1 m − 1 φ( R 0 )≤ S 1 (t,x),

and

u 0 (x)≤ ∥ u 0 ( x ) ∥ ∞ ≤ A 1 φ( R 0 )≤ S 1 (0,x),x∈ Ω ¯ .

By the comparison principle, we have

u(t,x)≤ S 1 (t,x),(t,x)∈[0,∞)× Ω ¯ .

The proof of Lemma 3.1 is complete. □

Lemma 3.2 (Upper bound)

Assume that l≥m>1, 1≤q<2, lq∈(m,m+1], and the initial data u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2). Then there exists C 2 >0 depending only on m, l, q, Ω, and ∥ u 0 ∥ ∞ such that

u(t,x)≤ C 2 ( 1 + t ) − 1 m − 1 ,(t,x)∈(0,∞)× Ω ¯ .
(3.4)

Proof Assume that x 0 ∉ Ω ¯ and R 0 >0 such that Ω⊂B( x 0 , R 0 ). For (t,x)∈(0,∞)× Ω ¯ , we define the function

S 2 (t,x)= A 2 ( 1 + δ t ) − 1 m − 1 σ(r),σ(r)= [ m m + 1 ( R m + 1 m − r m + 1 m ) ] 1 m with r=|x− x 0 |,

where the positive constants A 2 , R, and δ satisfy the following condition:

A 2 = [ m q − 1 ( 1 + m ( N − 1 ) ) 2 l q R ( m + 1 ) ( l − m ) q + m ( q + m − 1 ) m 2 ] 1 l q − m , R = [ R 0 m + 1 m + ( m + 1 ) ∥ u 0 ∥ ∞ m m A 2 m ] m m + 1 > R 0 and δ = ( m − 1 ) ( 1 + m ( N − 1 ) ) A 2 m − 1 2 m R 1 + m 2 m 2 .
(3.5)

Since x 0 ∉ Ω ¯ , the function S 2 (t,x) is C ∞ -smooth in [0,∞)× Ω ¯ . Moreover, for (t,x)∈(0,∞)×Ω and r=|x− x 0 |< R 0 <R, it follows from (3.5) and lq>m that

( 1 + δ t ) m m − 1 { ∂ t S 2 − △ S 2 m − | ∇ S 2 l | q } ( t , x ) = − A 2 δ m − 1 σ ( r ) − A 2 m [ N − 1 r ( σ m ) r + ( σ m ) r r ] − A 2 l q ( 1 + δ t ) − l q − m m − 1 | ( σ l ) r | q = − A 2 δ m − 1 σ ( r ) + A 2 m [ 1 m + N − 1 ] r − m − 1 m − A 2 l q ( 1 + δ t ) − l q − m m − 1 | ( ( σ m ) l m ) r | q = − A 2 δ m − 1 σ ( r ) + A 2 m [ 1 m + N − 1 ] r − m − 1 m − A 2 l q ( l m ) q R ( m + 1 ) ( l − m ) q + m q m 2 ≥ A 2 m [ ( N − 1 + 1 m ) R − m − 1 m − A 2 l q − m ( l m ) q R ( m + 1 ) ( l − m ) q + m q m 2 − δ R m + 1 m 2 ( m − 1 ) A 2 m − 1 ] ≥ A 2 m [ ( 1 2 ( N − 1 + 1 m ) R − m − 1 m − A 2 l q − m ( l m ) q R ( m + 1 ) ( l − m ) q + m q m 2 ) + ( 1 2 ( N − 1 + 1 m ) R − m − 1 m − δ R m + 1 m 2 ( m − 1 ) A 2 m − 1 ) ] ≥ 0 .
(3.6)

Therefore, the condition (3.5) guarantees that S 2 (t,x) is a supersolution to (1.1) in (0,∞)×Ω. In addition, since r=|x− x 0 |< R 0 <R for x∈Ω, for (t,x)∈(0,∞)×∂Ω, we deduce from (3.5) that

u(t,x)=0≤ A 2 ( 1 + t ) − 1 m − 1 σ( R 0 )≤ S 2 (t,x),

and

u 0 (x)≤ ∥ u 0 ( x ) ∥ ∞ ≤ A 2 σ( R 0 )≤ S 2 (0,x),x∈ Ω ¯ .

By the comparison principle, we have

u(t,x)≤ S 2 (t,x),(t,x)∈[0,∞)× Ω ¯ .

The proof of Lemma 3.2 is complete. □

Lemma 3.3 (Hölder estimate)

Assume that l≥m>1, 1≤q<2, lq∈[m,m+1], and the initial data u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2). Then there exists L 1 >0 depending only on m, l, q, Ω, and ∥ u 0 ∥ W 1 , ∞ ( Ω ) such that

| u ( t , x ) | = | u ( t , x ) − u ( t , x 0 ) | ≤ L 1 ( 1 + t ) 1 m − 1 | x − x 0 | 1 m , ( t , x , x 0 ) ∈ [ 0 , ∞ ) × Ω ¯ × ∂ Ω .
(3.7)

Proof Since the boundary ∂ Ω of Ω is smooth, there exists R Ω >0 such that for each x 0 ∈∂Ω, there exists y 0 ∈ R N satisfying | x 0 − y 0 |= R Ω and B( y 0 , R Ω )∩Ω=∅. It follows from the initial data condition u 0 (x)∈ W 1 , ∞ (Ω) that u 0 (x) is Lipschitz continuous, i.e., there exists L 0 >0 such that

| u 0 ( x ) − u 0 ( y ) | ≤ L 0 |x−y|,(x,y)∈ Ω ¯ × Ω ¯ .
(3.8)

Next, we define the open subset U κ , x 0 of R N by

U κ , x 0 = { x ∈ Ω : R Ω < | x − y 0 | < R Ω + κ } ,

where 0<κ<1 satisfies Ω κ ={x∈Ω:d(x,∂Ω)>κ}≠∅, and denote the function

S κ , x 0 (t,x)= A 3 ( 1 + t ) − 1 m − 1 ψ(r),ψ(r)= ( 1 − e − r κ ) 1 m with r=|x− y 0 |− R Ω ,

for (t,x)∈(0,∞)× U ¯ κ , x 0 and A 3 >0.

Moreover, we assume that

A 3 ≥max { e L 0 , e 1 m C ( e − 1 ) 1 m , ( 3 e m − 1 ) 1 m − 1 }
(3.9)

and

0<κ<min { 1 , R Ω 3 ( N − 1 ) , ( m q 3 l q A l q − m ) 1 2 − q } ,
(3.10)

where C=max{ C 1 , C 2 }>0, C 1 , and C 2 are defined in Lemmas 3.1 and 3.2, respectively.

Since y 0 ∉ U ¯ κ , x 0 , the function S κ , x 0 (t,x) is C ∞ -smooth in [0,∞)× U ¯ κ , x 0 . For (t,x)∈(0,∞)× U κ , x 0 , we have r=|x− y 0 |− R Ω ∈(0,κ). By a direct computation, we infer from (3.9)-(3.10) that

( 1 + t ) m m − 1 A 3 m { ∂ t S κ , x 0 − Δ S κ , x 0 m − | ∇ S κ , x 0 l | q } ( t , x ) = − 1 ( m − 1 ) A 3 m − 1 ψ ( r ) − [ N − 1 r + R Ω ( ψ m ) r + ( ψ m ) r r ] − A 3 l q − m ( 1 + t ) − l q − m m − 1 | ( ψ l ) r | q = − 1 ( m − 1 ) A 3 m − 1 ( 1 − e − r κ ) 1 m + 1 κ 2 e − r κ − N − 1 ( r + R Ω ) κ e − r κ − A 3 l q − m ( 1 + t ) − l q − m m − 1 | ( ( ψ m ) l m ) r | q ≥ − 1 ( m − 1 ) A 3 m − 1 + 1 κ 2 e − r κ − N − 1 κ R Ω e − r κ − A 3 l q − m | ( ( ψ m ) l m ) r | q ≥ 1 κ 2 e − r κ [ 1 − N − 1 R Ω κ − A 3 l q − m ( l m ) q κ 2 − q − e ( m − 1 ) A 3 m − 1 ] ≥ 1 κ 2 e − r κ [ ( 1 3 − N − 1 R Ω κ ) + ( 1 3 − A 3 l q − m ( l m ) q κ 2 − q ) + ( 1 3 − e ( m − 1 ) A 3 m − 1 ) ] ≥ 0 .

Therefore, S κ , x 0 (t,x) is a supersolution to (1.1) in (0,∞)× U κ , x 0 . In addition, it follows from (3.8)-(3.10) and the mean value theorem that

S κ , x 0 ( 0 , x ) = A 3 ( 1 − e − | x − y 0 | − R Ω κ ) 1 m ≥ A 3 ( 1 − e − | x − y 0 | − R Ω κ ) ≥ A 3 | x − y 0 | − R Ω e κ = A 3 e κ dist ( x , ∂ B ( y 0 , R Ω ) ) ≥ A 3 e κ dist ( x , ∂ Ω ) ≥ L 0 dist ( x , ∂ Ω ) ≥ u 0 ( x ) , x ∈ U ¯ δ , x 0 .

Moreover, for (t,x)∈(0,∞)×∂ U κ , x 0 , we have either x∈∂Ω or r=|x− y 0 |− R Ω =δ. If x∈∂Ω, then we have

S κ , x 0 (t,x)= A 3 ( 1 + t ) − 1 m − 1 ( 1 − e − r κ ) 1 m ≥0=u(x,t),(t,x)∈(0,∞)×∂Ω.

If r=|x− y 0 |− R Ω =δ, by Lemmas 3.1-3.2 and (3.9), we have

S κ , x 0 (t,x)= A 3 ( 1 + t ) − 1 m − 1 ( 1 − e − 1 ) 1 m ≥C ( 1 + t ) − 1 m − 1 ≥u(t,x).

By the comparison principle [6] we have

u(t,x)≤ S κ , x 0 (t,x),(t,x)∈[0,∞)× U ¯ κ , x 0 .

Consequently,

0≤u(t,x)−u(t, x 0 )=u(t,x)≤ A 3 ( 1 + t ) − 1 m − 1 ( 1 − e − | x − y 0 | − R Ω κ ) 1 m ,

for (t,x)∈[0,∞)× U ¯ κ , x 0 .

Since | x 0 − y 0 |− R Ω =0, we have

0≤ | u ( t , x ) − u ( t , x 0 ) | ≤ A 3 κ 1 m ( 1 + t ) 1 m − 1 | x − x 0 | 1 m ,x∈ U ¯ κ , x 0 .
(3.11)

Finally, we consider x∈Ω and x 0 ∈∂Ω. If |x− x 0 |≥ κ 2 , it follows from Lemmas 3.1-3.2 that

| u ( t , x ) − u ( t , x 0 ) | =u(t,x)≤ 2 1 m C κ 1 m ( 1 + t ) 1 m − 1 | x − x 0 | 1 m ,

where C=max{ C 1 , C 2 }>0, and C 1 and C 2 are defined in Lemmas 3.1 and 3.2, respectively.

If |x− x 0 |< κ 2 , let y 0 ∈ R N satisfy | x 0 − y 0 |= R Ω and B( y 0 , R Ω )∩Ω=∅. Since x∈Ω, we have

|x− y 0 |> R Ω and|x− y 0 |≤|x− x 0 |+| x 0 − y 0 |< R Ω +κ,

which implies x∈ U κ , x 0 . Therefore, we deduce from (3.11) that

| u ( t , x ) − u ( t , x 0 ) | =u(t,x)≤ A 3 κ 1 m ( 1 + t ) 1 m − 1 | x − x 0 | 1 m .

Choosing L 1 =max{ 2 1 m C κ 1 m , A 3 κ 1 m }, we have

| u ( t , x ) | = | u ( t , x ) − u ( t , x 0 ) | ≤ L 1 ( 1 + t ) 1 m − 1 | x − x 0 | 1 m ,(t,x, x 0 )∈[0,∞)× Ω ¯ ×∂Ω.

The proof of Lemma 3.3 is complete. □

We next proceed as in [29] to deduce the Hölder continuity of u(x,t) from Lemma 3.3. Therefore, we obtain the following corollary.

Corollary 3.1 Assume that l≥m>1, 1≤q<2, lq∈[m,m+1], and the initial data u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2). Then there exists L 2 >0 depending only on m, l, q, Ω, and ∥ u 0 ∥ W 1 , ∞ ( Ω ) such that

| u ( t , x ) − u ( t , y ) | ≤ L 2 ( 1 + t ) 1 m − 1 | x − y | 1 m ,(t,x,y)∈[0,∞)× Ω ¯ × Ω ¯ .
(3.12)

Proof The proof is similar to the argument in [25, 29], so we omit here. □

Proofs of Theorems 1.1 and 1.2 The proofs are based on the ideas in [25], but we give the details of the argument for the reader’s convenience. Let U(x,t) be the solution to the porous medium equation with homogeneous Dirichlet boundary conditions

{ ∂ t U − Δ U m = 0 , ( t , x ) ∈ ( 0 , ∞ ) × Ω , U ( t , x ) = 0 , ( t , x ) ∈ ( 0 , ∞ ) × ∂ Ω , U ( 0 , x ) = u 0 ( x ) , x ∈ Ω .
(3.13)

According to the nonnegativity of | ∇ u l | q , it follows from the comparison principle [6] that

0≤U(t,x)≤u(t,x),(t,x)∈(0,∞)× Ω ¯ .
(3.14)

We introduce the scaling variable s=lnt for t>0 and denote the new unknown function v and V by

u(t,x)= t − 1 m − 1 v(lnt,x),(t,x)∈(0,∞)× Ω ¯

and

U(t,x)= t − 1 m − 1 V(lnt,x),(t,x)∈(0,∞)× Ω ¯ .

Then v is a viscosity solution to the following problem:

{ ∂ s v − Δ v m − e − ( l q − m ) s m − 1 | ∇ v l | q − v m − 1 = 0 , ( s , x ) ∈ ( 0 , ∞ ) × Ω , v ( s , x ) = 0 , ( s , x ) ∈ ( 0 , ∞ ) × ∂ Ω , v ( 0 , x ) = u ( 1 , x ) , x ∈ Ω .
(3.15)

In addition, owing to Lemmas 3.1-3.3 and (3.14), we have

V(s,x)≤v(s,x)≤C,(s,x)∈[0,∞)× Ω ¯ ,
(3.16)

and

| v ( s , x ) − v ( s , y ) | ≤ L 1 | x − y | 1 m ,(s,x,y)∈[0,∞)× Ω ¯ ×∂Ω.
(3.17)

Next, for ε∈(0,1), we define

w ε (s,x)=v ( s ε , x ) ,(s,x)∈[0,∞)× Ω ¯ ,

and the half-relaxed limits

w ∗ (x)= lim ( σ , y , ε ) → ( s , x , 0 ) inf w ε (σ,y)and w ∗ (x)= lim ( σ , y , ε ) → ( s , x , 0 ) sup w ε (σ,y),

for (s,x)∈(0,∞)× Ω ¯ .

By (3.16), it is easy to see that w ∗ (x) and w ∗ (x) are well-defined and do not depend on s>0. Moreover, it readily follows from (3.15) and (3.17) that

w ∗ (x)= w ∗ (x)=0,x∈∂Ω.
(3.18)

By a direct computation, w ε (s,x) is a solution to the following initial-boundary problem:

{ ε ∂ s w ε − Δ w ε m − e − ( l q − m ) s ( m − 1 ) ε | ∇ w ε l | q − w ε m − 1 = 0 , ( s , x ) ∈ ( 0 , ∞ ) × Ω , w ε ( s , x ) = 0 , ( s , x ) ∈ ( 0 , ∞ ) × ∂ Ω , w ε ( 0 , x ) = u ( 1 , x ) , x ∈ Ω .
(3.19)

Next, we shall give proofs of Theorems 1.1 and 1.2. To do this, we distinguish the two cases lq=m and lq∈(m,m+1].

Case 1. lq=m. We use the stability of semicontinuous viscosity solutions [6] to deduce from (3.19) that

w ∗ (x) is a supersolution to (2.1) in Î©
(3.20)

and

w ∗ (x) is a subsolution to (2.1) in Î©.
(3.21)

In addition, by [30], V(s)→ f 0 in L ∞ (Ω) as s→∞. Moreover, it follows from (3.16) and the definition of w ∗ (x) and w ∗ (x) that

f 0 (x)≤ w ∗ (x)≤ w ∗ (x)≤C,x∈ Ω ¯ .
(3.22)

Since f 0 >0 in Ω by [30], we deduce from (3.22) that w ∗ (x) and w ∗ (x) are positive and bounded in Ω and vanish on ∂ Ω by (3.18). Owing to (3.20) and (3.21), we infer from Lemma 2.1 that

w ∗ (x)≤ w ∗ (x),x∈ Ω ¯ .

By (3.22), we have

w ∗ (x)= w ∗ (x),x∈ Ω ¯ .

Setting f(x)= w ∗ (x)= w ∗ (x), we deduce from (3.18), (3.20), (3.21), and (3.22) that f(x)∈ C 0 ( Ω ¯ ) is a positive viscosity solution to (2.1) and solves (1.9). Therefore, the existence of a positive solution to (1.9) is proved. Moreover, by Corollary 2.1, we obtain the uniqueness of solution to (1.9).

Furthermore, it follows from the equality w ∗ (x)= w ∗ (x) that

∥ w ε ( 1 ) − f ∥ ∞ →0as Îµâ†’0,

i.e.,

lim s → ∞ ∥ v ( s ) − f ∥ ∞ =0.
(3.23)

Therefore, we infer from the scaling transformation that

lim t → ∞ ∥ t 1 m − 1 u ( t ) − f ∥ ∞ =0.
(3.24)

Finally, Corollary 3.1 gives the last statement of Theorem 1.1. The proof of Theorem 1.1 is complete.

Case 2. lq∈(m,m+1]. We use once more the stability of semicontinuous viscosity solutions [6] to deduce from (3.19) that

w ∗ (x) is a supersolution to (2.10) in Î©
(3.25)

and

w ∗ (x) is a subsolution to (2.10) in Î©.
(3.26)

In addition, by [30], V(s)→ f 0 in L ∞ (Ω) as s→∞. Moreover, it follows from (3.16) and the definition of w ∗ (x) and w ∗ (x) that

f 0 (x)≤ w ∗ (x)≤ w ∗ (x)≤C,x∈ Ω ¯ .
(3.27)

Since f 0 >0 in Ω by [30] and is a solution to (2.10), we can apply Lemma 2.2 to conclude that w ∗ (x)≤ f 0 in Ω ¯ . Owing to (3.27), we have

w ∗ (x)= w ∗ (x)= f 0 (x),x∈ Ω ¯ .

Therefore, we deduce from (3.19) that

∥ w ε ( 1 ) − f 0 ∥ ∞ →0as Îµâ†’0,

i.e.,

lim s → ∞ ∥ v ( s ) − f 0 ∥ ∞ =0.
(3.28)

Thus, we infer from the scaling transformation that

lim t → ∞ ∥ t 1 m − 1 u ( t ) − f 0 ∥ ∞ =0.
(3.29)

The proof of Theorem 1.2 is complete. □

4 Proof of Theorem 1.3

In this section, we shall consider the well-posedness and the large time behavior of solutions to (1.1) with the small initial data u 0 (x) for lq>m+1 by the method used in [25]. To do this, we need the following lemma.

Lemma 4.1 Let l≥m>1 and 1≤q<2. Assume that G is the corresponding solution to (1.1) with the initial data G 0 ∈ W 1 , ∞ (Ω) satisfying (1.2) for the case lq=m, and denote

F(t,x)= 1 ℓ [ G 0 ] G ( t ℓ [ G 0 ] m − 1 , x ) ,(t,x)∈[0,∞)× Ω ¯ ,
(4.1)

where the parameter ℓ[ G 0 ] is defined in (1.10). Then F(t,x) is a solution to (1.1) with the initial data G 0 /ℓ[ G 0 ] and lq=m such that |∇ F l |≤1 for (t,x)∈[0,∞)× Ω ¯ . Moreover, F(t,x) is a supersolution to (1.1) for lq>m.

Proof According to the definition of ℓ[ G 0 ] in (1.10), it is easy to see that |∇ F l |≤1 for (t,x)∈[0,∞)× Ω ¯ .

Next, let φ(t,x)∈ C 2 ((0,∞)×Ω) and F−φ has a local minimum at ( t 0 , x 0 )∈(0,∞)×Ω. Since ℱ is smooth with respect to the time variable and Hölder continuous with respect to the space variable, we obtain

| ∇ φ l ( t 0 , x 0 ) | ≤1.
(4.2)

Moreover, introducing ψ(t,x):=ℓ[ G 0 ]φ(ℓ [ G 0 ] m − 1 t,x) for (t,x)∈[0,∞)× Ω ¯ , the function G−ψ has a local minimum at ( t 0 ℓ [ G 0 ] 1 − m , x 0 ) such that

∂ t ψ ( t 0 ℓ [ G 0 ] 1 − m , x 0 ) −Δ ψ m ( t 0 ℓ [ G 0 ] 1 − m , x 0 ) − | ∇ ψ l ( t 0 ℓ [ G 0 ] 1 − m , x 0 ) | m l ≥0,

i.e.,

∂ t φ( t 0 , x 0 )−Δ φ m ( t 0 , x 0 )− | ∇ φ l ( t 0 , x 0 ) | m l ≥0.
(4.3)

Therefore, F(t,x) is a supersolution to (1.1) with lq=m. In a similar way, it can be shown that F(t,x) is also a subsolution. Hence, F(t,x) is a solution to (1.1) with lq=m.

Furthermore, we deduce from (4.2), (4.3), and lq>m that

∂ t φ ( t 0 , x 0 ) − Δ φ m ( t 0 , x 0 ) − | ∇ φ l ( t 0 , x 0 ) | q ≥ | ∇ φ l ( t 0 , x 0 ) | m l ( 1 − | ∇ φ l ( t 0 , x 0 ) | q − m l ) ≥ 0 .
(4.4)

The proof of Lemma 4.1 is complete. □

Proposition 4.1 Let l≥m>1, 1≤q<2, and lq>m+1. Assume that the initial data u 0 (x)∈ W 1 , ∞ (Ω) satisfies (1.2), moreover, there exists G 0 ∈ W 1 , ∞ (Ω) satisfying (1.2) such that

u 0 (x)≤ G 0 ( x ) ℓ [ G 0 ] ,x∈ Ω ¯ ,
(4.5)

where the parameter â„“[ G 0 ] is defined in (1.10). Then there exists a unique solution u to (1.1) in the sense of Definition  1.1 and it satisfies

u(t,x)≤F(t,x),(t,x)∈[0,∞)× Ω ¯ ,
(4.6)

where F(t,x) is defined in (4.1).

Proof On the one hand, the solution U to the porous medium equation (3.13) is clearly a subsolution to (1.1) in (0,∞)×Ω.

On the other hand, it follows from Lemma 4.1 that the function F(t,x) is a supersolution to (1.1) in (0,∞)×Ω. Therefore, F(t,x) is a supersolution to (3.13).

Since U=F=0 on (0,∞)×∂Ω and U(0,x)= u 0 (x)≤F(0,x) for x∈ Ω ¯ by (4.5), we infer from the comparison principle [6] that

U(t,x)≤F(t,x),(t,x)∈[0,∞)× Ω ¯ .

This property and the simultaneous vanishing of U and ℱ on (0,∞)×∂Ω allow us to use the classical Perron method to establish the existence of a solution u(x,t) to (1.1) in the sense of Definition 1.1 which satisfies (4.6). The uniqueness next follows from the comparison principle [6]. The proof of Proposition 4.1 is complete. □

Proof of Theorem 1.3 We notice that Lemma 3.2 is also valid in that case. It readily follows from Lemma 3.3 and Proposition 4.1 that

0 ≤ u ( t , x ) = u ( t , x ) − u ( t , x 0 ) ≤ F ( t , x ) ≤ L 1 ( 1 + t ) 1 m − 1 | x − x 0 | 1 m , ( t , x , x 0 ) ∈ [ 0 , ∞ ) × Ω ¯ × ∂ Ω .

The convergence proof is similar to that performed in the proof of Theorem 1.2 for lq∈(m,m+1]. The proof of Theorem 1.3 is complete. □

5 Proof of Theorem 1.4

In this section, when lq>m+1, we shall prove that the solution u(x,t) of (1.1) blows up in finite time for the large initial data u 0 (x) in the sense of weak solution by the method used in [26].

In order to obtain a blow-up condition corresponding to (1.1), we have to modify the function e − ε | x | 2 used in [31–33], and introduce a test function ϕ ε (x) as follows:

Ï• ε (x)= A ε e − ε | x | with  A ε = 1 ∫ Ω e − ε | x | d x .

Proof of Theorem 1.4 Suppose that u(x,t) is the solution of the problem (1.1) and T is the blow-up time of the solution. For s∈(0, 1 2 ), we denote

W(t)= 1 s ∫ Ω u s ϕ ε (x)dx,t∈(0,T).
(5.1)

By a direct calculation, we have

W ′ ( t ) = ∫ Ω u s − 1 ϕ ε u t d x = ∫ Ω u s − 1 ϕ ε Δ u m d x + ∫ Ω u s − 1 ϕ ε | ∇ u l | q d x ≥ − ( s − 1 ) ∫ Ω u s − 2 ϕ ε ∇ u m ⋅ ∇ u d x + ε ∫ Ω u s − 1 ϕ ε ∇ u m ⋅ x | x | d x + ∫ Ω u s − 1 ϕ ε | ∇ u l | q d x ≥ m ( 1 − s ) ∫ Ω u m + s − 3 ϕ ε | ∇ u | 2 d x − m ε ∫ Ω u m + s − 2 ϕ ε | ∇ u | d x + ∫ Ω u s − 1 ϕ ε | ∇ u l | q d x .
(5.2)

By Young’s inequality, we obtain

ε ∫ Ω u m + s − 2 ϕ ε |∇u|dx≤ 1 2 ∫ Ω u m + s − 3 ϕ ε | ∇ u | 2 dx+ ε 2 2 ∫ Ω u m + s − 1 ϕ ε dx.
(5.3)

Since 0<s< 1 2 , it follows from (5.2), (5.3), and Poincaré’s inequality that

W ′ ( t ) ≥ ∫ Ω u s − 1 ϕ ε | ∇ u l | q d x − m ε 2 2 ∫ Ω u m + s − 1 ϕ ε d x = ( q l q l + s − 1 ) q ∫ Ω ϕ ε | ∇ u q l + s − 1 q | q d x − m ε 2 2 ∫ Ω u m + s − 1 ϕ ε d x ≥ ( q l q l + s − 1 ) q ∫ Ω ϕ ε u q l + s − 1 d x − m ε 2 2 ∫ Ω u m + s − 1 ϕ ε d x .
(5.4)

According to ql>m+1, 0<s< 1 2 , ∫ Ω ϕ ε (x)dx=1 and Hölder’s inequality, we have

∫ Ω u m + s − 1 ϕ ε dx= ∫ Ω u m + s − 1 ϕ ε m + s − 1 q l + s − 1 ϕ ε q l − m q l + s − 1 dx≤ ( ∫ Ω u q l + s − 1 ϕ ε d x ) m + s − 1 q l + s − 1 .
(5.5)

Thus, by (5.4) and (5.5), we obtain

d W d t ≥ ( ∫ Ω u q l + s − 1 ϕ ε d x ) m + s − 1 q l + s − 1 [ ( q l q l + s − 1 ) q ( ∫ Ω u q l + s − 1 ϕ ε d x ) q l − m q l + s − 1 − m ε 2 2 ] .
(5.6)

Owing to ql>m+1, 0<s< 1 2 , ∫ Ω ϕ ε (x)dx=1, and Jensen’s inequality, we have

( ∫ Ω u q l + s − 1 ϕ ε d x ) q l − m q l + s − 1 ≥ ( ∫ Ω u s ϕ ε d x ) q l − m s .
(5.7)

Therefore, it follows from (5.6) and (5.7) that

d W d t ≥ 1 2 ( q l q l + s − 1 ) q ( ∫ Ω u s ϕ ε d x ) q l + s − 1 s = 1 2 ( q l q l + s − 1 ) q s q l + s − 1 s W q l + s − 1 s (t)
(5.8)

as long as

W(t)≥ 1 s ( q l + s − 1 q l ) q s q l − m ( m ε 2 ) s q l − m for all t∈[0,T).
(5.9)

Taking

K 0 = ( q l + s − 1 q l ) q s q l − m ( m ε 2 ) s q l − m .

Since the initial data u 0 (x) satisfies

∫ Ω u 0 s (x) ϕ ε (x)dx≥ K 0 ,

we have

W(0)≥ 1 s K 0 = 1 s ( q l + s − 1 q l ) q s q l − m ( m ε 2 ) s q l − m .
(5.10)

Therefore, w(t) increases and remains above 1 s K 0 for all t∈[0,T].

By (5.8), integrating over (0,t) yields

W(t)≥ ( W 1 − q l s ( 0 ) − K 1 t ) − s q l − 1 with  K 1 = q l − 1 2 ( q l q l + s − 1 ) q s q l − 1 s .
(5.11)

Hence, it follows from (5.10) and (5.11) that the solution u(x,t) of (1.1) blows up in finite time, T= 1 K 1 W 1 − q l s (0).

Moreover, by (5.10), we obtain the upper estimate on the blow-up time T of the solution u(x,t) as follows:

T≤ 2 q l − 1 ( q l q l + s − 1 ) q ( m − 1 ) q l − m ( m ε 2 ) 1 − q l q l − m .
(5.12)

The proof of Theorem 1.4 is complete. □

References

  1. Andreucci D: Degenerate parabolic equations with initial data measures. Trans. Am. Math. Soc. 1997, 349: 3911-3923. 10.1090/S0002-9947-97-01530-4

    Article  MathSciNet  Google Scholar 

  2. Gilding BH, Guedda M, Kersner R: The Cauchy problem for u t =Δu+ | ∇ u | q . J. Math. Anal. Appl. 2003, 284: 733-755. 10.1016/S0022-247X(03)00395-0

    Article  MathSciNet  Google Scholar 

  3. Li YX, Souplet P: Single-point gradient blow-up on the boundary for diffusive Hamilton-Jacobi equations in planar domains. Commun. Math. Phys. 2010, 293: 499-517. 10.1007/s00220-009-0936-8

    Article  MathSciNet  Google Scholar 

  4. Bardi M, Capuzzo-Dolcetta I: Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations. Birkhäuser, Boston; 1997.

    Book  Google Scholar 

  5. Barles G Mathématiques et Applications 17. In Solutions de Viscosité des Equations d’Hamilton-Jacobi. Springer, Berlin; 1994.

    Google Scholar 

  6. Crandall MG, Ishii H, Lions PL: User’s guide to viscosity solutions of second order partial differential equations. Bull. Am. Math. Soc. (N.S.) 1992, 27: 1-67. 10.1090/S0273-0979-1992-00266-5

    Article  MathSciNet  Google Scholar 

  7. Arrieta JM, Bernal AR, Souplet P: Boundedness of global solutions for nonlinear parabolic equations involving gradient blow-up phenomena. Ann. Sc. Norm. Super. Pisa, Cl. Sci. 2004, 3: 1-15.

    MathSciNet  Google Scholar 

  8. Benachour SD, Dăbuleanu-Hapca S, Laurençot P: Decay estimates for a viscous Hamilton-Jacobi equation with homogeneous Dirichlet boundary conditions. Asymptot. Anal. 2007, 51: 209-229.

    MathSciNet  Google Scholar 

  9. Ding JT, Guo BZ: Global existence and blow-up solutions for quasilinear reaction-diffusion equations with a gradient term. Appl. Math. Lett. 2011, 24: 936-942. 10.1016/j.aml.2010.12.052

    Article  MathSciNet  Google Scholar 

  10. Souplet P, Zhang QS: Global solutions of inhomogeneous Hamilton-Jacobi equations. J. Anal. Math. 2006, 99: 355-396. 10.1007/BF02789452

    Article  MathSciNet  Google Scholar 

  11. Xia L, Yao ZA: Existence, uniqueness and asymptotic behavior of solutions for a singular parabolic equation. J. Math. Anal. Appl. 2009, 358: 182-188. 10.1016/j.jmaa.2009.04.039

    Article  MathSciNet  Google Scholar 

  12. Zhou WS, Lei PD: A one-dimensional nonlinear heat equation with a singular term. J. Math. Anal. Appl. 2010, 368: 711-726. 10.1016/j.jmaa.2010.03.066

    Article  MathSciNet  Google Scholar 

  13. Barles G, Souganidis PE: On the large time behaviour of solutions of Hamilton-Jacobi equations. SIAM J. Math. Anal. 2000, 31: 925-939. 10.1137/S0036141099350869

    Article  MathSciNet  Google Scholar 

  14. Benachour S, Dabuleanu S: Large time behavior for a viscous Hamilton-Jacobi equation with Neumann boundary condition. J. Differ. Equ. 2005, 216: 223-258. 10.1016/j.jde.2005.02.017

    Article  MathSciNet  Google Scholar 

  15. Laurençot P: Convergence to steady states for a one-dimensional viscous Hamilton-Jacobi equation with Dirichlet boundary conditions. Pac. J. Math. 2007, 230: 347-364. 10.2140/pjm.2007.230.347

    Article  Google Scholar 

  16. Namah G, Roquejoffre JM: Remarks on the long time behaviour of the solutions of Hamilton-Jacobi equations. Commun. Partial Differ. Equ. 1999, 24: 883-893. 10.1080/03605309908821451

    Article  MathSciNet  Google Scholar 

  17. Qi YW, Wang MX: The self-similar profiles of generalized KPZ equation. Pac. J. Math. 2001, 201: 223-240. 10.2140/pjm.2001.201.223

    Article  MathSciNet  Google Scholar 

  18. Quittner P, Souplet P: Superlinear Parabolic Problems. Blow-up, Global Existence and Steady States. Birkhäuser, Basel; 2007.

    Google Scholar 

  19. Roquejoffre JM: Convergence to steady states or periodic solutions in a class of Hamilton-Jacobi equations. J. Math. Pures Appl. 2001, 80: 85-104. 10.1016/S0021-7824(00)01183-1

    Article  MathSciNet  Google Scholar 

  20. Shi PH, Wang MX: Global solution of the fast diffusion equation with gradient absorption terms. J. Math. Anal. Appl. 2007, 326: 602-621. 10.1016/j.jmaa.2006.02.094

    Article  MathSciNet  Google Scholar 

  21. Tchamba TT: Large time behavior of solutions of viscous Hamilton-Jacobi equations with superquadratic Hamiltonian. Asymptot. Anal. 2010, 66: 161-186.

    MathSciNet  Google Scholar 

  22. Gilding BH: The Cauchy problem for u t =Δu+ | ∇ u | q : large-time behaviour. J. Math. Pures Appl. 2005, 84: 753-785. 10.1016/j.matpur.2004.11.003

    Article  MathSciNet  Google Scholar 

  23. Stinner C: Convergence to steady states in a viscous Hamilton-Jacobi equation with degenerate diffusion. J. Differ. Equ. 2010, 248: 209-228. 10.1016/j.jde.2009.09.019

    Article  MathSciNet  Google Scholar 

  24. Barles G, Laurençot P, Stinner C: Convergence to steady states for radially symmetric solutions to a quasilinear degenerate diffusive Hamilton-Jacobi equation. Asymptot. Anal. 2010, 67: 229-250.

    MathSciNet  Google Scholar 

  25. Laurençot P, Stinner C: Convergence to separate variable solutions for a degenerate parabolic equation with gradient source. J. Dyn. Differ. Equ. 2012, 24: 29-49. 10.1007/s10884-011-9238-x

    Article  Google Scholar 

  26. Mu CL, Li YH, Wang Y: Life span and a new critical exponent for a quasilinear degenerate parabolic equation with slow decay initial values. Nonlinear Anal., Real World Appl. 2010, 11: 198-206. 10.1016/j.nonrwa.2008.10.048

    Article  MathSciNet  Google Scholar 

  27. Laurençot P, Stinner C: Refined asymptotics for the infinite heat equation with homogeneous Dirichlet boundary conditions. Commun. Partial Differ. Equ. 2011, 36: 532-546.

    Article  Google Scholar 

  28. Da Lio F: Comparison results for quasilinear equations in annular domains and applications. Commun. Partial Differ. Equ. 2002, 27: 283-323. 10.1081/PDE-120002788

    Article  MathSciNet  Google Scholar 

  29. Kawohl B, Kutev N: Comparison principle and Lipschitz regularity for viscosity solutions of some classes of nonlinear partial differential equations. Funkc. Ekvacioj 2000, 43: 241-253.

    MathSciNet  Google Scholar 

  30. Manfredi JJ, Vespri V: Large time behavior of solutions to a class of doubly nonlinear parabolic equations. Electron. J. Differ. Equ. 1994, 2: 1-17.

    MathSciNet  Google Scholar 

  31. Guo JS, Guo YY: On a fast diffusion equation with source. Tohoku Math. J. 2001, 53: 571-579. 10.2748/tmj/1113247801

    Article  MathSciNet  Google Scholar 

  32. Li YH, Mu CL: Life span and a new critical exponent for a degenerate parabolic equation. J. Differ. Equ. 2004, 207: 392-406. 10.1016/j.jde.2004.08.024

    Article  MathSciNet  Google Scholar 

  33. Mukai K, Mochizuki K, Huang Q: Large time behavior and life span for a quasilinear parabolic equation with slowly decaying initial values. Nonlinear Anal. TMA 2000, 39: 33-45. 10.1016/S0362-546X(98)00161-8

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors would like to express sincere gratitude to the referees for their valuable suggestions and comments on the original manuscript. The first author is supported by Scientific Research Fund of Sichuan Provincial Science and Technology Department (2014JY0098); the second author is supported in part by the Fundamental Research Funds for the Central Universities, Project No. CDJXS 12 10 00 14; the third author is supported in part by NSF of China (11371384).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Pan Zheng.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The authors declare that the study was realized in collaboration with the same responsibility. All authors 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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Li, N., Zheng, P., Mu, C. et al. Large time behavior of solutions for the porous medium equation with a nonlinear gradient source. Bound Value Probl 2014, 98 (2014). https://doi.org/10.1186/1687-2770-2014-98

Download citation

  • Received:

  • Accepted:

  • Published:

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

Keywords