Open Access
Issue
Wuhan Univ. J. Nat. Sci.
Volume 29, Number 2, April 2024
Page(s) 95 - 105
DOI https://doi.org/10.1051/wujns/2024292095
Published online 14 May 2024

© Wuhan University 2024

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

0 Introduction

In the present article, our attention is focused on the blow-up behavior of the following porous medium equation with superline source and local linear boundary dissipation

{ u t = Δ u m + | u | p - 2 u , x Ω , t > 0 , u ( x , t ) = 0 , x Γ 0 , t > 0 , u m n = - u t , x Γ 1 , t > 0 , u ( x , 0 ) = u 0 ( x ) , x Ω , (1)

where m>1, p>m+1, Ω is a bounded open subset of Rn with C1 boundary Ω and n represents the unit outer normal vector to Ω. {Γ0,Γ1} is a partition of the boundary Ω satisfying Ω=Γ0Γ1,Γ0¯Γ1¯=.

Moreover, Γ0 and Γ1 are measurable over Ω, endowed with (n-1)-dimensional surface measure σ and σ(Γ0)>0.

Problem (1) with m=1 can be regarded as a mathematical model to depict a heat reaction-diffusion process that occurs inside a solid body Ω surrounded by a fluid, with contact Γ1 and having an internal cavity with a contact boundary Γ0. In this physics background, the quantity of heat produced by the reaction is proportional to a superlinear power of the temperature. To avoid an internal explosion in Ω, a refrigeration system is installed in the fluid. The operational mechanism of this refrigeration system lies in the fact that the heat absorbed from the fluid is proportional to the power of the rate of change of the temperature, which can be expressed as: un=-|ut|q-2ut,xΓ1, where un stands for the heat flux from Ω to the fluid.

Evolution equations [1-8] with boundary damping have attracted the attention of mathematicians in the past period. For instance, Fiscella and Vitillaro [9] studied the following problem with local nonlinear boundary dissipation

{ u t - Δ u = | u | p - 2 u , x Ω , t > 0 , u ( x , t ) = 0 x Γ 0 , t > 0 , u n = - | u t | q - 2 u t , x Γ 1 , t > 0 , u ( x , 0 ) = u 0 ( x ) , x Ω ¯ , (2)

where 2p1+2*2 and 2* denotes the Sobolev conjugate of 2. Using the monotonicity method of Lions[9] and a contraction argument, they proved the local well-posedness in the Hadamard sense. Moreover, in the case of a superlinear source, i.e. p>2, under the condition

J ( u 0 ) = 1 2 Ω | u 0 ( x ) | 2 d x - 1 p | u 0 ( x ) | p d x < d : = i n f u H Γ 1 1 ( Ω ) { 0 } s u p λ > 0 J ( λ u ) ,

they gave the global existence and finite time blow-up results. To be precise, if

K ( u 0 ) = Ω | u 0 ( x ) | 2 d x - Ω | u 0 ( x ) | p d x 0 ,

then the weak solution is global, while if

K ( u 0 ) = Ω | u 0 ( x ) | 2 d x - Ω | u 0 ( x ) | p d x 0

and q<q0(p):=2(n+1)p-4(n-1)n(p-2)+4, then the weak solution will blow up in some finite time. In a recent work, the authors [10] considered problem (1) with m=1, and obtained the finite time blow-up result for arbitrary high initial energy. Moreover, under some additional conditions, the authors also gave estimates of the blow-up time. In addition, using some differential inequality techniques, the authors [11, 12] considered the lower bounds for the blow-up time of blow-up solutions to some porous medium equations with null Dirichlet boundary conditions or homogeneous Neumann boundary conditions.

To the best of our knowledge, there is no previous work on the blow-up behavior of the solution to the problem (1). Building on the aforementioned work, we will analyze the effects of the nonlinear diffusion and the local linear boundary dissipation on the blow-up phenomenon of problem (1). In order to deal with the difficulties caused by the nonlinear diffusion term Δum better and more effectively, throughout this paper, we work with the following equivalent formulation of the problem (1) obtained by changing variables um=v,

{ ( v 1 m ) t = Δ v + | v | p - 2 m v 1 m , x Ω , t > 0 , v ( x , t ) = 0 , x Γ 0 , t > 0 , v n = - ( v 1 m ) t , x Γ 1 , t > 0 , v ( x , 0 ) = v 0 ( x ) = u 0 m ( x ) , x Ω . (3)

First, we obtain the finite time blow-up criterion of the solution to the problem (3) by using a modified concavity method (Theorem 1). Second, for any aR,we prove that there exists a v0HΓ01(Ω) with initial energy J(v0)=a that leads to a finite time blow-up solution (Corollary 1). Finally, the lower bounds of the blow-up time are derived by combining the interpolation inequality for Lq-norms, the Sobolev embedding theorem, with some differential inequality techniques (Theorem 2).

The article is organized as follows: In Section 1, we introduce some notations and state some useful lemmas. In Section 2, we give the finite time blow-up criterion and the lower and upper estimates of the blow-up time.

1 Preliminaries

In this section, we first introduce some notations, definitions, and some known results. Throughout this article, we denote q=Lq(Ω), q,Γ1=Lq(Γ1) for some q[1,+], and the Hilbert space

H Γ 0 1 ( Ω ) = { v H 1 ( Ω ) : v = 0   f o r   x Γ 0 } .

( , ) and (,)Γ1 stand for the inner products on the Hilbert spaces L2(Ω) and L2(Γ1) , respectively. From the trace theorem, one knows that there exists a continuous trace mapping HΓ01(Ω)L2(Ω). Moreover, since σ(Γ0)>0, then, Theorem 6.7-5 [13] tells us that a Poincaré-type inequality holds. Therefore, v2 is equivalent to the norm vHΓ01=(v22+v22)12 in the space HΓ01(Ω). On the other hand, since m>1, one can define the following positive optimal constants

S 1 = s u p v H Γ 0 1 v 0 ( Ω ) v m + 1 2 m 2 , Γ 1 2 v 2 2 , S 2 = s u p v H Γ 0 1 v 0 ( Ω ) v m + 1 2 m 2 2 v 2 2 . (4)

The definition of the weak solution to the problem (3) is given as follows.

Definition 1

Suppose that v0HΓ01(Ω), and m+1<p<m(2*-1)+1, where 2* represents the Sobolev conjugate of 2, namely

2 * = { ,   i f   n 2 ; 2 n n - 2 ,   i f   n > 2 .

A function v:=v(x,t)L(0,T;HΓ01(Ω)) with (v1m)tL2((0,T)×Ω) is called a weak solution of problem (3) on [0,T)×Ω if v(x,0)=v0=u0m, and

( ( v 1 m ) t , ϕ ) + ( v ( t ) , ϕ ) + ( ( v 1 m ) t , ϕ ) Γ 1 = Ω | v | p - 2 m v 1 m ϕ d x (5)

holds for a.e. t(0,T) and any ϕHΓ01(Ω). Moreover, the spatial trace of v has a distributional time derivative vt on (0,T)×Ω, belonging to L2((0,T)×Ω).

Definition 2

Suppose that v is a weak solution of problem (3). We say that v blows up in finite time T if

l i m t T - ρ ( t ) = l i m t T - 1 m + 1 ( v m + 1 2 m ( t ) 2 2 + v m + 1 2 m ( t ) 2 , Γ 1 2 ) = + .

In what follows, we introduce the energy-functional J(v)=12v22-mm+p-1vm+p-1mm+p-1m, and Nehari functional

K ( v ) = v 2 2 - v m + p - 1 m m + p - 1 m

in HΓ01(Ω) related to problem (3). Evidently, one knows that J(v) and K(v) are continuous on HΓ01(Ω), and for a.e. t(0,T),

d d t J ( v ( t ) ) = - 4 m ( m + 1 ) 2 ( ( v m + 1 2 m ) t ( t ) 2 2 + ( v m + 1 2 m ) t ( t ) 2 , Γ 1 2 ) 0 (6)

which implies that

J ( v ( t ) ) = J ( v 0 ) - 4 m ( m + 1 ) 2 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ (7)

Now, we introduce the definition of the potential well-depth d=infvNJ(v)=infvHΓ01(Ω)\{0}supλ>0J(λv), where N is the Nehari manifold N={vHΓ01(Ω)|K(v)=0}\{0}. Indeed, d also can be characterized as

d = p - m - 1 2 ( m + p - 1 ) ( s u p v H Γ 0 1 v 0 ( Ω ) v m + p - 1 m v 2 ) - 2 ( m + p - 1 ) p - m - 1 .

We are now able to give some lemmas, which will play a key role in our proof of the main results.

Lemma 1

Suppose that m+1<p<m(2*-1)+1. Then one has vm+p-1mm+p-1m>2(m+p-1)(p-m-1)d for any vN-={vHΓ01(Ω)|K(v)<0}.

Proof

Since v22-vm+p-1mm+p-1m=K(v)<0, one knows that v0. Meanwhile, one can verify that K(λ*v)=0 with

λ * = ( v 2 2 v m + p - 1 m m + p - 1 m ) m p - m - 1 ( 0,1 ) .

Thereupon, one has λ*vN. Furthermore, according to the definition of the potential well depth d, one can arrive at

d = i n f v N J ( v )   J ( λ * v ) = p - m - 1 2 ( p + m - 1 ) λ * v m + p - 1 m m + p - 1 m + 1 2 K ( λ * v )    = p - m - 1 2 ( p + m - 1 ) ( λ * ) m + p - 1 m v m + p - 1 m m + p - 1 m < p - m - 1 2 ( m + p - 1 ) v m + p - 1 m m + p - 1 m (8)

which results in vm+p-1mm+p-1m>2(m+p-1)(p-1-m)d.

The proof of the Lemma 1 is completed.

Lemma 2

Suppose that m+1<p<m(2*-1)+1, and the weak solution v of problem (3) blows up in finite time T. Then there is a t*[0,T) such that v(t*)N-.

Proof

Suppose that K(v(t))0 for any t[0,T). Then, from (6), it follows that

J ( v 0 ) J ( v ( t ) ) = p - m - 1 2 ( m + p - 1 ) v ( t ) 2 2 + m m + p - 1 K ( v ( t ) ) p - m - 1 2 ( m + p - 1 ) v ( t ) 2 2 ,

which contradicts the assumption that v is a finite-time blow-up weak solution. The proof of Lemma 2 is completed.

Lemma 3 [14] Suppose that a positive function F on [0,T] satisfies the following conditions:(i)F is differentiable on [0,T] and F' is absolutely continuous on [0,T] with F'(0)>0; (ii) there exists a positive constant α>0 such that

F ( t ) F   ( t ) - ( 1 + α ) ( F ' ( t ) ) 2 0

holds for a.e. t[0,T]. Then TF(0)αF'(0).

2 The Finite Time Blow-up Results

Recalling that ρ(t)=1m+1(vm+12m(t)22+vm+12m(t)2,Γ12), one can easily check that

d d t ρ ( t ) = 1 m [ ( v 1 m ( t ) , v t ( t ) ) + ( v 1 m ( t ) , v t ( t ) ) Γ 1 ] = - K ( v ( t ) ) (9)

For any T1(0,T), we define an auxiliary function as the form

F ( t ) = 0 t ρ ( τ ) d τ + ( T 1 - t ) ρ ( 0 ) + m β m + 1 ( t + σ ) 2 , t [ 0 , T 1 ] (10)

where β and σ are two positive parameters which will be determined later. It is clear that F is positive on [0,T1]. By a simple calculation, one has

F ' ( t ) = ρ ( t ) - ρ ( 0 ) + 2 m β m + 1 ( t + σ ) = 0 t d d t ρ ( τ ) d τ + 2 m β m + 1 ( t + σ ) = 1 m 0 t [ ( v 1 m , v τ ) + ( v 1 m , v τ ) Γ 1 ] d τ + 2 m β m + 1 ( t + σ ) (11)

and F'(0)=2mβσm+1>0. Moreover, from (6) and (7), it follows that

F   ( t ) = d d t ρ ( t ) + 2 m β m + 1 = - K ( v ( t ) ) + 2 m β m + 1 = - ( m + p - 1 m J ( v ( t ) - p - m - 1 2 m v ( t ) 2 2 ) + 2 m β m + 1 = p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v 0 ) + 2 m β m + 1 + 4 ( m + p - 1 ) ( m + 1 ) 2 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ (12)

On the other hand, employing the Young's inequality and Cauchy-Schwartz inequality, we obtain:

ξ ( t ) = [ 1 m 0 t ( v m + 1 2 m ( τ ) 2 2 + v m + 1 2 m ( τ ) 2 , Γ 1 2 ) d τ + β ( t + σ ) 2 ] [ 4 m ( m + 1 ) 2 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ + 4 m 2 β ( m + 1 ) 2 ] - [ 1 m 0 t [ ( v 1 m , v τ ) + ( v 1 m , v τ ) Γ 1 ] d τ + 2 m β m + 1 ( t + σ ) ] 2

{ 2 m + 1 ( 0 t [ ( v m + 1 2 m , ( v m + 1 2 m ) τ ) + ( v m + 1 2 m , ( v m + 1 2 m ) τ ) Γ 1 ] d τ ) + 2 m β m + 1 ( t + σ ) } 2 - [ 1 m 0 t [ ( v 1 m , v τ ) + ( v 1 m , v τ ) Γ 1 ] d τ + 2 m β m + 1 ( t + σ ) ] 2 = [ 1 m 0 t [ ( v 1 m , v τ ) + ( v 1 m , v τ ) Γ 1 ] d τ + 2 m β m + 1 ( t + σ ) ] 2 - [ 1 m 0 t [ ( v 1 m , v τ ) + ( v 1 m , v τ ) Γ 1 ] d τ + 2 m β m + 1 ( t + σ ) ] 2 = 0 (13)

Now, we are in the position to estimate FF-λ(F')2with λ=m+p-1m+1>1. Combining (10), (11), and (12), one can arrive at:

F F - λ ( F ' ) 2 = F F - λ [ 1 m 0 t [ ( v 1 m , v τ ) + ( v 1 m , v τ ) Γ 1 ] d τ + 2 m β m + 1 ( t + σ ) ] 2 = F F + λ [ ξ ( t ) - m + 1 m ( F ( t ) - ( T - t ) ρ ( 0 ) ) ( 4 m ( m + 1 ) 2 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ + 4 m 2 β ( m + 1 ) 2 ) ] F F - 4 m λ β m + 1 F ( t ) - 4 λ m + 1 F ( t ) 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ = F [ p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v 0 ) + 2 m β m + 1 + 4 ( m + p - 1 ) ( m + 1 ) 2 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ      - 4 λ m + 1 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ - 4 m λ β m + 1 ] = F [ ( 4 ( m + p - 1 ) ( m + 1 ) 2 - 4 λ m + 1 ) 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ (14)

Noticing that, λ=m+p-1m+1, then (14) leads to

F F - m + p - 1 m + 1 ( F ' ) 2 F [ p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 ] (15)

Up to now, from the above discussion, one can summarize the following lemma.

Lemma 4

Suppose that ρ(0)>0 and there is a positive constant such that

p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 0

holds for any t(0,T1]. Then 0<T1(m+1)3ρ(0)m(p-2)2β.

Proof

Since m+p-1m+1=1+p-2m+1>1, then, a direct application of Lemma 3 and (15) yields that

T 1 < F ( 0 ) p - 2 m + 1 F ' ( 0 ) = T 1 ρ ( 0 ) + m β σ 2 m + 1 p - 2 m + 1 2 m β σ m + 1 = ρ ( 0 ) ( m + 1 ) 2 2 m ( p - 2 ) β σ T 1 + m + 1 2 ( p - 2 ) σ (16)

Choosing σ(ρ(0)(m+1)22m(p-2)β,+) to guarantee ρ(0)(m+1)22m(p-2)βσ<1, then (16) tells us that,

T 1 ( 1 - ρ ( 0 ) ( m + 1 ) 2 2 m ( p - 2 ) β σ ) - 1 m + 1 2 ( p - 2 ) σ (17)

By a series of calculations, one can verify that the right side of (17) takes its minimum at the point

σ = σ β = ρ ( 0 ) ( m + 1 ) 2 m ( p - 2 ) β ( ρ ( 0 ) ( m + 1 ) 2 2 m ( p - 2 ) β , + ) .

In other words, one has T1(m+1)3ρ(0)m(p-2)2β.

The proof of Lemma 4 is completed.

Now, we give the finite time blow-up criteria as follows.

Theorem 1

Suppose that m+1<p<m(2*-1)+1 and the initial data v0 belongs to one of the following sets:

1 = { v H Γ 0 1 : J ( v ) < d , K ( v ) < 0 } , 2 = { v H Γ 0 1 : J ( v ) < ( p - m - 1 ) ( v m + 1 2 m 2 2 + v m + 1 2 m 2 , Γ 1 2 ) 2 ( m + p - 1 ) ( S 1 + S 2 ) } .

Then, the weak solution v to the problem (3) blows up in a finite time. Moreover, one has

T 2 m ( 2 p + m - 3 ) ( p - 2 ) 2 ( m + p - 1 ) v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 d - J ( v 0 )

for v01 and

T 2 m ( 2 p + m - 3 ) ( p - 2 ) 2 ( m + p - 1 ) v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 p - m - 1 2 ( m + p - 1 ) ( S 1 + S 2 ) ( v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 ) - J ( v 0 )

for v02.

Proof

If v01. Then, with the help of (7), one has, for any t[0,T),

J ( v ( t ) ) = J ( v 0 ) - 4 m ( m + 1 ) 2 0 t ( ( v m + 1 2 m ) τ ( τ ) 2 2 + ( v m + 1 2 m ) τ ( τ ) 2 , Γ 1 2 ) d τ J ( v 0 ) < d (18)

Suppose that there exists a t1[0,T) such that K(v(t1))=0 and K(v(t))<0 for any t[0,t1). Thereupon, Lemma 1 and the continuity of the mapping tv(t)2 are applicable to produce

v ( t 1 ) 2 2 = l i m t t 1 - v ( t ) 2 2 2 ( m + p - 1 ) p - m - 1 d > 0 ,

which implies that v(t1)N. And hence, one has J(v(t1))infvNJ(v)=d, which contradicts (18). That is to say, for any t[0,T), one can claim that v(t)1 provided that v01.

Selecting

β ( 0 , ( m + p - 1 ) ( m + 1 ) 2 2 m 2 ( 2 p + m - 3 ) ( d - J ( v 0 ) ) ] ,

and keeping Lemma 1 in mind, one has, for any t(0,T),

p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 > p - m - 1 2 m 2 ( m + p - 1 ) p - 1 - m d - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 = m + p - 1 m d - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 0 .

A direct application of Lemma 4 tells us that,

0 < T 2 m ( m + 1 ) ( 2 p + m - 3 ) ρ ( 0 ) ( p - 2 ) 2 ( m + p - 1 ) ( d - J ( v 0 ) ) = 2 m ( 2 p + m - 3 ) ( p - 2 ) 2 ( m + p - 1 ) v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 d - J ( v 0 ) .

If v02. Then from (4) and (9), it follows that,

d d t ρ ( t ) = - K ( v ( t ) ) = p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v ( t ) ) p - m - 1 2 m ( m + 1 ) ρ ( t ) S 1 + S 2 - m + p - 1 m J ( v ( t ) ) = m + p - 1 A m ( ρ ( t ) - A J ( v ( t ) ) ) (19)

where,

A = 2 ( m + p - 1 ) ( S 1 + S 2 ) ( m + 1 ) ( p - m - 1 ) > 0 .

Putting H(t)=ρ(t)-AJ(v(t)) and combining (6) with (19), one has, for a.e. t(0,T),

d d t H ( t ) = d d t ρ ( t ) - A d d t J ( v ( t ) ) d d t ρ ( t ) m + p - 1 A m H ( t ) .

Integrating the above inequality from 0 to t results in:

H ( t ) e m + p - 1 A m H ( 0 ) . (20)

On the other hand, the assumption v02 implies that,

H ( 0 ) = ρ ( 0 ) - A J ( v 0 ) = 1 m + 1 ( v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 ) - 2 ( m + p - 1 ) ( S 1 + S 2 ) ( m + 1 ) ( p - m - 1 ) J ( v 0 ) > 0 (21)

Combining (19), (20) with (21) yields that, for a.e. t(0,T),

d d t ρ ( t ) m + p - 1 A m H ( t ) m + p - 1 A m e m + p - 1 A m H ( 0 ) > 0 ,

which means that ρ(t) is increasing in [0,T). Therefore, one can see that,

p - m - 1 2 m v ( t ) 2 2 - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 p - m - 1 2 m ( m + 1 ) ρ ( t ) S 1 + S 2 - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 p - m - 1 2 m ( m + 1 ) ρ ( 0 ) S 1 + S 2 - m + p - 1 m J ( v 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 = m + p - 1 A m H ( 0 ) - 2 m β ( 2 p + m - 3 ) ( m + 1 ) 2 > 0 (22)

provided that

β ( 0 , ( m + p - 1 ) ( m + 1 ) 2 H ( 0 ) 2 A m 2 ( 2 p + m - 3 ) ] .

According to Lemma 4, one can obtain the following estimate,

0 < T 2 A m ( m + 1 ) ( 2 p + m - 3 ) ρ ( 0 ) ( p - 2 ) 2 ( m + p - 1 ) H ( 0 ) ,

namely,

0 < T 2 m ( 2 p + m - 3 ) ( p - 2 ) 2 ( m + p - 1 ) v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 p - m - 1 2 ( m + p - 1 ) ( S 1 + S 2 ) ( v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 ) - J ( v 0 ) .

The proof of Theorem 1 is completed.

In fact, Theorem 1 tells us the sets 1 and 2 are invariant under the semi-flow associated with problem (3). Namely, if v01, then v(t)1, while v02, then v(t)2. On the other hand, with the help of (4), for any v2, one has

K ( v ) = v 2 2 - v m + p - 1 m m + p - 1 m = m + p - 1 m J ( v ) - p - m - 1 2 m v 2 2 m + p - 1 m ( J ( v ) - ( p - m - 1 ) ( v m + 1 2 m 2 2 + v m + 1 2 m 2 , Γ 1 2 ) 2 ( m + p - 1 ) ( S 1 + S 2 ) ) < 0 .

Therefore, one can claim that v(t)N-={vHΓ01|K(v)<0} for any t[0,T) provided initial data v012. Based on the above arguments, it is natural to ask whether or not the condition v0N- is sufficient enough for a finite-time blow-up. This is a difficult task, and the authors [15] conducted a similar study.

In addition, one can know that both 1 and 2 are non-empty sets. Moreover, Corollary 1 implies that for any aR, there exists a v0HΓ01(Ω) with initial energy J(v0)=a, which leads to finite time blow-up solution.

Corollary 1

For any aR, denote the energy level set by:

J a = { v H Γ 0 1 ( Ω ) | J ( v ) = a } .

Then Ja2.

Proof

Suppose that Ω1 and Ω2 are two disjoint open subdomains of Ω, and

d i s t ( Ω 1 ¯ , Ω ) > 0 , d i s t ( Ω 2 ¯ , Ω ) > 0 , d i s t ( Ω 1 ¯ , Ω 2 ¯ ) > 0 .

From the proof of Theorem 3.7 [16], one knows that there is a sequence {vk}H01(Ω1) such that,

1 2 Ω 1 | v k ( x ) | 2 d x - m m + p - 1 Ω 1 | v k ( x ) | m + p - 1 m d x +   a s   k (23)

On the other hand, choosing an arbitrary nonzero function ωC0(Ω) with supp ωΩ2, then,

a - ( 1 2 Ω 2 | ( r ω ( x ) ) | 2 d x - m m + p - 1 Ω 2 | r ω ( x ) | m + p - 1 m d x ) +   a s   r (24)

and there exists a r0>0 such that

p - m - 1 2 ( p + m - 1 ) ( S 1 + S 2 ) Ω 2 | r m + 1 2 m ω m + 1 2 m ( x ) | 2 d x = r m + 1 m ( p - m - 1 ) 2 ( p + m - 1 ) ( S 1 + S 2 ) Ω 2 | ω ( x ) | m + 1 m d x > a (25)

holds for any r>r0. By (23) and (24), there are k0Z+ and r1>r0 such that

1 2 Ω 1 | v k 0 ( x ) | 2 d x - m m + p - 1 Ω 1 | v k 0 ( x ) | m + p - 1 m d x = a - ( 1 2 Ω 2 | ( r 1 ω ( x ) ) | 2 d x - m m + p - 1 Ω 2 | r 1 ω ( x ) | m + p - 1 m d x ) (26)

Let v0=v˜+r1ω, where,

v ˜ = { 0 , x Ω \ Ω 1 v k 0 ( x ) , x Ω 1 .

It is not difficult to show that v0HΓ01(Ω) and v0(x)=0 in Ω\(Ω1Ω2). From (25) and (26), it follows that,

J ( v 0 ) = 1 2 ( Ω 1 + Ω 2   ) | v 0 ( x ) | 2 d x - m m + p - 1 ( Ω 1 + Ω 2   ) | v 0 ( x ) | m + p - 1 m d x = ( 1 2 Ω 1 | v k 0 ( x ) | 2 d x - m m + p - 1 Ω 1 | v k 0 ( x ) | m + p - 1 m d x ) + ( 1 2 Ω 2 | ( r 1 ω ( x ) ) | 2 d x - m m + p - 1 Ω 2 | r 1 ω ( x ) | m + p - 1 m d x ) = ( 26 ) a < ( 25 ) p - m - 1 2 ( p + m - 1 ) ( S 1 + S 2 ) Ω 2 | r 1 m + 1 2 m ω m + 1 2 m ( x ) | 2 d x = p - m - 1 2 ( p + m - 1 ) ( S 1 + S 2 ) Ω 2 | v 0 m + 1 2 m ( x ) | 2 d x p - m - 1 2 ( m + p - 1 ) ( S 1 + S 2 ) ( v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 ) (27)

which means that v0Ja2. The proof of Corollary 1 is completed.

Theorem 2

Suppose that n>2,m+1<p<(m+1)(1+2n), the weak solution v of problem (3) blows up in finite time T and v(t)N- for any t[0,T). Then

T S ˜ ( v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 ) 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] ,

where σ(0,1) and S˜>0 are two constants given by (28) and (33), respectively.

Proof

Since m>1 and n>2, one can infer that,

( m + 1 ) ( 1 + 2 n ) < 2 m n n - 2 - m + 1 .

Furthermore, from the assumption p<(m+1)(1+2n), it follows that,

1 < m + 1 m < m + p - 1 m < 2 n n - 2 ,

and

σ = ( m m + p - 1 - n - 2 2 n ) ( m m + 1 - n - 2 2 n ) - 1 ( 0,1 ) (28)

Noticing that v(t)N- for any t[0,T), the interpolation inequality for Lq-norms and the embedding H1(Ω)L2nn-2(Ω) can be used to obtain:

v ( t ) 2 2 < v ( t ) m + p - 1 m m + p - 1 m v ( t ) 2 n n - 2 ( m + p - 1 ) ( 1 - σ ) m v m + 1 2 m ( t ) 2 2 ( m + p - 1 ) σ m + 1 C ( m + p - 1 ) ( 1 - σ ) m v ( t ) 2 ( m + p - 1 ) ( 1 - σ ) m v m + 1 2 m ( t ) 2 2 ( m + p - 1 ) σ m + 1 (29)

where C denotes the optimal constant of the embedding H1(Ω)L2nn-2(Ω). Keeping p<(m+1)(1+2n) in mind, one can show that 2-(m+p-1)(1-σ)m>0. And hence, (29) implies that

v ( t ) 2 C ( m + p - 1 ) ( 1 - σ ) 2 m - ( m + p - 1 ) ( 1 - σ ) v m + 1 2 m ( t ) 2 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] . (30)

From (9), (29) and (30), it follows that,

d d t ρ ( t ) = - K ( v ( t ) ) = v ( t ) m + p - 1 m m + p - 1 m - v ( t ) 2 2 < v ( t ) m + p - 1 m m + p - 1 m C ( m + p - 1 ) ( 1 - σ ) m v ( t ) 2 ( m + p - 1 ) ( 1 - σ ) m v m + 1 2 m ( t ) 2 2 ( m + p - 1 ) σ m + 1 C 2 ( m + p - 1 ) ( 1 - σ ) 2 m - ( 1 - σ ) ( m + p - 1 ) v m + 1 2 m ( t ) 2 4 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] < C 2 ( m + p - 1 ) ( 1 - σ ) 2 m - ( 1 - σ ) ( m + p - 1 ) ( v m + 1 2 m ( t ) 2 2 + v m + 1 2 m ( t ) 2 , Γ 1 2 ) 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] = S 3 [ ρ ( t ) ] 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] (31)

where,

S 3 = C 2 ( m + p - 1 ) ( 1 - σ ) 2 m - ( 1 - σ ) ( m + p - 1 ) ( m + 1 ) 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] .

It is not difficult to verify that,

2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] > 1 .

From (31), it follows that,

[ ρ ( t ) ] 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] - [ ρ ( 0 ) ] 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] > S 3 [ 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] ] t (32)

Letting tT, then (32) results in:

- [ ρ ( 0 ) ] 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] > S 3 [ 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] ] T ,

which means that,

T > [ ρ ( 0 ) ] 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] S 3 [ 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] - 1 ] = S ˜ ( v 0 m + 1 2 m 2 2 + v 0 m + 1 2 m 2 , Γ 1 2 ) 1 - 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] ,

where,

S ˜ = S 3 - 1 [ 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] - 1 ] - 1 ( m + 1 ) 2 m σ ( m + p - 1 ) ( m + 1 ) [ 2 m - ( m + p - 1 ) ( 1 - σ ) ] - 1 (33)

The proof of Theorem 2 is completed.

References

  1. Cavalcanti M M, Domingos Cavalcanti V N, Lasiecka I. Well-posedness and optimal decay rates for the wave equation with nonlinear boundary damping-source interaction[J]. Journal of Differential Equations, 2007, 236(2): 407-459. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  2. Chueshov I, Eller M, Lasiecka I. On the attractor for a semilinear wave equation with critical exponent and nonlinear boundary dissipation[J]. Communications in Partial Differential Equations, 2002, 27(9/10): 1901-1951. [CrossRef] [MathSciNet] [Google Scholar]
  3. Hintermann T. Evolution equations with dynamic boundary conditions[J]. Proceedings of the Royal Society of Edinburgh: Section A Mathematics, 1989, 113(1/2): 43-60. [CrossRef] [MathSciNet] [Google Scholar]
  4. Vázquez J L, Vitillaro E. Heat equation with dynamical boundary conditions of reactive-diffusive type[J]. Journal of Differential Equations, 2011, 250(4): 2143-2161. [CrossRef] [MathSciNet] [Google Scholar]
  5. Vitillaro E. On the wave equation with hyperbolic dynamical boundary conditions, interior and boundary damping and source[J]. Archive for Rational Mechanics and Analysis, 2017, 223(3): 1183-1237. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  6. Yang X, Zhou Z F. Blow-up problems for the heat equation with a local nonlinear Neumann boundary condition[J]. Journal of Differential Equations, 2016, 261(5): 2738-2783. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  7. Xue Y Z. Lower bounds of blow up time for a class of slow reaction diffusion equations with inner absorption terms[J].Wuhan University Journal of Natural Sciences, 2023, 28(5): 373-378. [CrossRef] [EDP Sciences] [Google Scholar]
  8. Liu D M, Chen A. Global existence and extinction behaviour for a doubly nonlinear parabolic equation with logarithmic nonlinearity[J].Wuhan University Journal of Natural Sciences, 2023, 28(2): 99-105. [CrossRef] [EDP Sciences] [Google Scholar]
  9. Vitillaro E, Fiscella A. Local Hadamard well: Posedness and blow up for reaction: Diffusion equations with nonlinear dynamical boundary conditions[J]. Discrete and Continuous Dynamical Systems, 2013, 33(11/12): 5015-5047. [CrossRef] [MathSciNet] [Google Scholar]
  10. Sun F L, Wang Y T, Yin H J. Blow-up problems for a parabolic equation coupled with superlinear source and local linear boundary dissipation[J]. Journal of Mathematical Analysis and Applications, 2022, 514(2): 126327. [CrossRef] [Google Scholar]
  11. Huo W T, Fang Z B. Life span bounds for reaction-diffusion equation with a space-time integral source term[J]. Zeitschrift Fiir Angewandte Mathematik Und Physik, 2023, 74(4): 128. [NASA ADS] [CrossRef] [Google Scholar]
  12. Liu D M, Mu C L, Xin Q. Lower bounds estimate for the blow-up time of a nonlinear nonlocal porous medium equation[J]. Acta Mathematica Scientia, 2012, 32(3): 1206-1212. [CrossRef] [MathSciNet] [Google Scholar]
  13. Ciarlet P G. Linear and Nonlinear Functional Analysis with Applications[M]. Philadelphia: Society for Industrial and Applied Mathematics, 2013. [CrossRef] [Google Scholar]
  14. Levine H A. Some nonexistence and instability theorems for solutions of formally parabolic equations of the form Formula [J]. Archive for Rational Mechanics and Analysis, 1973, 51(5): 371-386. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  15. Dickstein F, Mizoguchi N, Souplet P, et al. Transversality of stable and Nehari manifolds for a semilinear heat equation[J]. Calculus of Variations and Partial Differential Equations, 2011, 42(3): 547-562. [CrossRef] [MathSciNet] [Google Scholar]
  16. Willem M. Minimax Theorems, Progress in Nonlinear Differential Equations and Their Applications[M]. Boston: Birkhäuser Boston Inc, 1996. [Google Scholar]

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.