Preprint
Article

The Second Critical Exponent for a Time Fractional Reaction-Diffusion Equation

Altmetrics

Downloads

91

Views

33

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

27 July 2024

Posted:

29 July 2024

You are already at the latest version

Alerts
Abstract
In this paper, we consider the Cauchy problem of a time fractional nonlinear diffusion equation. According to the Kaplan’s first eigenvalue method, we first prove the blow-up of the solutions in finite time for some sufficient conditions. We next give sufficient conditions for the existence of global solutions by using the result of Zhang and Sun. In conclusions, we find the second critical exponent for the existence of global and non-global solutions via the decay rates of the initial data at spatial infinity.
Keywords: 
Subject: Computer Science and Mathematics  -   Mathematics

MSC:  26A33; 35A01; 35B44; 35K15; 35R11

1. Introduction

We study the Cauchy problem for a time fractional reaction-diffusion equation
t α u = Δ u + | u | p 1 u , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) 0 , x R n ,
where n 1 , 0 < α < 1 , p > 1 , u 0 C 0 ( R n ) : = { f C ( R n ) ; lim | x | f ( x ) = 0 } , and t α denotes the Caputo time fractional derivative of order α defined by
t α u ( t ) = 1 Γ ( 1 α ) 0 t ( t s ) α u s ( x , s ) d s , 0 < α < 1 .
Here, Γ ( · ) is the Gamma function. Moreover, the Caputo time fractional derivative (2) is related to the Riemann-Liouville derivative by
t α u ( t ) = 1 Γ ( 1 α ) d d t ( t α [ u u 0 ] ) ( t ) = t   0 I t 1 α ( u ( x , t ) u 0 ( x ) ) ,
where   0 I t 1 α denotes left Riemann-Liouville fractional integrals of order 1 α and is defined by
  0 I t 1 α u = 1 Γ ( 1 α ) 0 t ( t s ) α u ( s ) d s .
For a given initial data u 0 , let T = T ( u 0 ) be the maximal existence time of the solution of (1). If T = , the solution is global in time. However, if T < , then the solution is not global in time in the sense that it blows up at t = T such as
lim sup t T u ( · , t ) L ( R n ) = .
A lot of significant results on the critical exponents for nonlinear parabolic equations have been obtained during the past decades. Fujita [1] considered the following Cauchy problem:
u t = Δ u + u p , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) 0 , x R n .
In [1], it is shown that (4) possesses the critical Fujita exponent p = 1 + 2 / n such that
  • If 1 < p < p , then the solution blows up in finite time for any nontrivial initial data.
  • If p > p , then there are both global solutions and nonglobal solutions corresponding to small and large initial data, respectively.
According to Hayakawa [2], Kobayashi et al. [3], and Weissler [4], it has been known that p = p = 1 + 2 / n belongs to the blow-up case. In some situations, the size of the initial data required by the global and nonglobal solutions can be determined through the so-called second critical exponent with respect to the decay rates of the initial data as | x | . When p > p = 1 + 2 / n , Lee and Ni [5] established the second critical exponent a = 2 / ( p 1 ) for (4) with initial data u 0 ( x ) = λ ψ ( x ) , where λ > 0 and ψ ( x ) is a bounded continuous function in R n , such that the following conditions hold:
  • If lim inf | x | | x | a ψ ( x ) > 0 for some a ( 0 , a ) and any λ > 0 , then the solution blows up in finite time.
  • If lim sup | x | | x | a ψ ( x ) < for some a ( a , n ) , then there is λ 0 > 0 such that the solution is global in time whenever λ ( 0 , λ 0 ) .
Lee and Ni [5] proved that a = 2 / ( p 1 ) belongs to the global case.
The weighted source case
u t = Δ u + K ( x ) u p , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) 0 , x R n .
with K ( x ) > 0 of the order | x | σ for σ > 1 if n = 1 or for σ > 2 if n 2 was considered with the critical Fujita exponent p = 1 + ( 2 + σ ) / n by Pinsky [6].
The degenerate case
u t = Δ u m + u p , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) 0 , x R n ,
with m > 1 and max ( 0 , 1 2 / n ) < m < 1 was thoroughly studied with the critical Fujita exponent p = m + 2 / n by Galaktionov et al. [7], Qi [8], and Mochizuki and Mukai [9]. Furthermore, Galaktionov [10], Mochizuki and Mukai [9], Kawanago [11], and Mochizuki and Suzuki [12] have shown that p = p = m + 2 / n belongs to the blow-up case. When p > p = m + 2 / n , Mukai et al. [13], and Guo and Guo [14] obtained the second critical exponent a = 2 / ( p m ) for (5).
The extended case
u t = Δ u m + K ( x , t ) u p , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) 0 , x R n .
with K ( x , t ) = t s | x | σ for s 0 , m > max ( 0 , 1 2 / n ) , p > max ( 1 , m ) , σ > 1 if n = 1 , or σ > 2 if n 2 was obtained with the critical Fujita exponent p = m + s ( m 1 ) + ( 2 + 2 s + σ ) / n by Qi [15]. In the case with K ( x , t ) = K ( x ) of the order | x | σ as | x | with σ R in some cone D and K ( x , t ) = 0 if otherwise, Suzuki [16] considered (6) for 1 m < p and obtained the critical Fujita exponent p = m + { 2 + max ( σ , n ) } / n and the second critical exponent a = { 2 + max ( σ , n ) } / ( p m ) .
Winkler [17] considered the nonlinear diffusion equation not in divergence form
u t = u p Δ u + u q , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) > 0 , x R n ,
and obtained the critical Fujita exponent q = p + 1 for p 1 such that
  • Suppose that 1 q < q ( 1 q < 3 / 2 if p = 1 ). If u 0 decreases sufficiently fast in space, all positive solutions of (7) are global and unbounded.
  • Suppose that q = q . Then, all positive solutions of (7) blow-up in finite time.
  • Suppose that q > q . If u 0 is sufficiently large, then any positive solution of (7) blows up in finite time. If u 0 ( x ) f ( | x | ) in R n , then the solutions of (7) are global, where f satisfies for u ( x , t ) = ( 1 + t ) α f ( ( 1 + t ) β | x | )
    f + n 1 r + β r f p f + α f 1 p + f q p = 0 , r ( 0 , ) , f ( 0 ) = f 0 , f ( 0 ) = 0 ,
    with α = 1 / ( q 1 ) , β = ( q p 1 ) / ( 2 q 2 ) , f 0 > 0 , r = ( 1 + t ) β | x | .
Furthermore, Li and Mu [18] also considered (7) and obtained the second critical exponent a = 2 / ( q p 1 ) for p > 1 and q > p + 1 + 2 / n with initial data u 0 ( x ) = λ ψ ( x ) , where λ > 0 and ψ ( x ) is a bounded continuous function in R n , such that
  • Let n 2 . Assume that lim inf | x | | x | a ψ ( x ) > 0 . If 0 < a < a , or a = a and λ is large enough, then the solution u ( x , t ) of (7) blows up in finite time.
  • Assume that lim sup | x | | x | a ψ ( x ) < . If a > a , then there exists λ 0 > 0 such that the solution u ( x , t ) of (7) is global in time whenever λ ( 0 , λ 0 ) .
Yang et al. [19] and the author [20] studied the extended case
u t = u p Δ u + K ( x ) u q , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) > 0 , x R n ,
with the positive weight function K C 0 ( R n ) satisfying
c 1 | x | σ K ( x ) c 2 | x | σ , | x | > R 0 , σ > 2
for some R 0 , c 1 , c 2 > 0 . Then, Yang et al. [19] obtained the critical Fujita exponent q = p + 1 for p 1 such that
  • Suppose that 1 < q < q ( 1 < q < 3 / 2 if p = 1 ). If u 0 decreases sufficiently fast in space, all positive solutions of (8) are global and unbounded.
  • Suppose that q = q . Then, all positive solutions of (8) blow-up in finite time.
  • Suppose that q > q . If u 0 is sufficiently large, then any positive solution of (8) blows up in finite time. If u 0 ( x ) f ( | x | ) in R n , then the solutions of (8) are global, where f satisfies for u ( x , t ) ( 1 + t ) α f ( ( 1 + t ) β | x | )
    f + n 1 r + β r f p f + α f 1 p + C 1 r σ f q p = 0 , r ( 0 , ) , f ( 0 ) = f 0 , f ( 0 ) = 0 ,
    with 2 < σ < 0 , C 1 = sup R n ( K ( x ) | x | σ ) , α = ( σ + 2 ) / ( 2 q 2 + p σ ) , β = ( q p 1 ) / ( 2 q 2 + p σ ) , f 0 > 0 , r = ( 1 + t ) β | x | , and
    f + n 1 r + β r f p f + α f 1 p + C 2 max ( r σ , 1 ) f q p = 0 , r ( 0 , ) , f ( 0 ) = f 0 , f ( 0 ) = 0 ,
    with σ 0 , C 2 = sup R n ( K ( x ) min ( | x | σ , 1 ) ) .
Moreover, Yang et al. [19] and the author [20] obtained the second critical exponent a = ( 2 + σ ) / ( q p 1 ) for p 1 and q > p + 1 with initial data u 0 ( x ) = λ ψ ( x ) , where λ > 0 and ψ ( x ) is a bounded continuous function in R n , such that
  • Assume that lim inf | x | | x | a ψ ( x ) > 0 . If 0 < a < a with σ > 2 , or a = a with σ > 2 and λ is large enough, then the solution u ( x , t ) of (8) blows up in finite time.
  • Assume that lim sup | x | | x | a ψ ( x ) < . If a > a with σ 0 , or a a with 2 < σ < 0 , then there exists λ 0 > 0 such that the solution u ( x , t ) of (8) is global in time whenever λ ( 0 , λ 0 ) .
By reading the literature of time fractional nonlinear diffusion equations, we found that there are no studies on the second critical exponent to the Cauchy problem (1). In this paper, we give the second critical exponent to the Cauchy problem (1) based on the above mentioned literature.
The rest of this paper is organized as follows. In Section 2, we give some preliminaries for the Cauchy problem (1). In Section 3, two sufficient conditions for the blow-up of solutions in finite time are given in Theorem 1. In Section 4, we state the existence of global solutions according to some conditions in Theorem 2. In Section 5, conclusions is presented.

2. Preliminaries

In this section, we give some preliminaries.
We need the following Wright type function
ϕ α ( z ) = k = 0 ( z ) k k ! Γ ( α k + 1 α ) , 0 < α < 1 , z C .
The function ϕ α is an entire function and satisfies the following properties:
(a)
ϕ α ( θ ) 0 for θ 0 and 0 ϕ α ( θ ) d θ = 1 .
(b)
0 ϕ α ( θ ) θ r d θ = Γ ( 1 + r ) Γ ( 1 + α r ) for r > 1 .
The operator A = Δ generates a semigroup { T ( t ) } t 0 on C 0 ( R n ) with domain
D ( A ) = { u C 0 ( R n ) : Δ u C 0 ( R n ) } .
Then T ( t ) is an analytic and contractive semigroup on C 0 ( R n ) , and
[ T ( t ) u 0 ] ( x ) = R n G ( x y , t ) u 0 ( y ) d y ,
where
G ( x , t ) = 1 ( 4 π t ) n / 2 e | x | 2 / ( 4 t ) , x R n , t > 0 .
For t 0 , we define the operators P α ( t ) and S α ( t ) as
[ P α ( t ) u 0 ] ( x ) = 0 ϕ α ( θ ) [ T ( t α θ ) u 0 ] ( x ) d θ
= 0 ϕ α ( θ ) R n G ( x y , t α θ ) u 0 ( y ) d y d θ ,   α ( t ) u 0 ] ( x ) = α 0 θ ϕ α ( θ ) [ T ( t α θ ) u 0 ] ( x ) d θ
= α 0 θ ϕ α ( θ ) R n G ( x y , t α θ ) u 0 ( y ) d y d θ ,
where ϕ α ( θ ) is the function defined by (9). Note that for given x R n { 0 } and t > 0 , G ( x , t α θ ) 0 as θ 0 . Hence, 0 ϕ α ( θ ) G ( x , t α θ ) d θ is well defined. Since 0 ϕ α ( θ ) d θ = 1 and R n G ( x , t ) d x = 1 , we know that
0 ϕ α ( θ ) R n G ( x , t α θ ) d x d θ = 1 f o r t > 0 .
Consider the linear equation
t α u = Δ u + f ( x , t ) , x R n , t > 0 , u ( x , 0 ) = u 0 ( x ) 0 , x R n ,
where u 0 C 0 ( R n ) and f L 1 ( ( 0 , T ) , C 0 ( R n ) ) . If u is a solution of (12), then by [21], it satisfies
u ( x , t ) = [ P α ( t ) u 0 ] ( x ) + 0 t ( t s ) α 1 [ S α ( t s ) f ( · , s ) ] ( x ) d s ,
where P α ( t ) and S α ( t ) are given by (10) and (11), respectively.
Zhang and Sun [21] obtained the following lemmas related to the operators P α ( t ) and S α ( t ) .
Lemma 1 
([21]). If u 0 ( x ) 0 , u 0 ( x ) ¬ 0 , then [ P α ( t ) u 0 ] ( x ) > 0 , [ S α ( t ) u 0 ] ( x ) > 0 and
P α ( t ) u 0 L 1 ( R n ) = u 0 L 1 ( R n ) , S α ( t ) u 0 L 1 ( R n ) = 1 Γ ( α ) u 0 L 1 ( R n ) .
Proof. 
See Lemmas 2.1 (a) and 2.2 (a) in [21]. □
Lemma 2 
([21]). Let 1 p q + and 1 / r = 1 / p 1 / q .
(a)
If 1 / r < 2 / n , then
P α ( t ) u 0 L q ( R n ) ( 4 π t α ) n 2 r Γ ( 1 n 2 r ) Γ ( 1 α n 2 r ) u 0 L p ( R n ) .
(b)
If 1 / r < 4 / n , then
S α ( t ) u 0 L q ( R n ) α ( 4 π t α ) n 2 r Γ ( 2 n 2 r ) Γ ( 1 + α α n 2 r ) u 0 L p ( R n ) .
Proof. 
See Lemmas 2.1 (b) and 2.2 (b) in [21]. □
Next, we give the definition of a mild solution of the Cauchy problem (1).
Definition 1. 
Let u 0 C 0 ( R n ) and T > 0 . We call that u C ( [ 0 , T ] , C 0 ( R n ) ) is a mild solution of the problem (1) if u satisfies the integral equation
u ( x , t ) = [ P α ( t ) u 0 ] ( x ) + 0 t ( t s ) α 1 [ S α ( t s ) | u ( · , s ) | p 1 u ( · , s ) ] ( x ) d s , t [ 0 , T ] .
For the Cauchy problem (1), Zhang and Sun [21] has established the following local existence result.
Proposition 1 (Theorem 3.2 in [21]). Let 0 < α < 1 . For given u 0 C 0 ( R n ) , there exists a maximal existence time T > 0 such that the problem (1) has a unique mild solution u C ( [ 0 , T ] , C 0 ( R n ) ) and either T = + or T < + and u L ( ( 0 , t ) , C 0 ( R n ) ) . In addition, if u 0 ( x ) 0 and u 0 ( x ) ¬ 0 , then u ( x , t ) [ P α ( t ) u 0 ] ( x ) > 0 for t ( 0 , T ) . Moreover, if u 0 L r ( R n ) for some r [ 1 , ) , then u C ( [ 0 , T ) , L r ( R n ) ) .
Furthermore, Zhang and Sun [21] also obtained the following blow-up and global existence results.
Proposition 2 (Theorem 4.3 in [21]). Let 0 < α < 1 , u 0 C 0 ( R n ) and u 0 ( x ) 0 . If
R n u 0 ( x ) χ ( x ) d x > 1 , w h e r e χ ( x ) = R n e n 2 + | x | 2 d x 1 e n 2 + | x | 2 ,
then the mild solutions of (1) blow up in a finite time.
Proposition 3 (Theorem 4.4 in [21]). Let 0 < α < 1 , u 0 C 0 ( R n ) , u 0 ( x ) 0 and u 0 ( x ) ¬ 0 .
(a)
If 1 < p < 1 + 2 / n , then the mild solution of (1) blows up in a finite time.
(b)
If p 1 + 2 / n and u 0 L q c ( R n ) is sufficiently small, where q c = n ( p 1 ) / 2 , then the mild solution of (1) exists globally.

3. Blow-Up of Solution

In this section, we shall prove the following blow-up result.
Theorem 1. 
Let n 1 and 0 < α < 1 . Assume that the initial data u 0 ( x ) = λ ψ ( x ) 0 , where λ > 0 and ψ C 0 ( R n ) . Suppose that one of the following two conditions holds:
(a)
λ > 0 is large enough;
(b)
0 < a < 2 / ( p 1 ) and
lim inf | x | | x | a ψ ( x ) > 0 .
Then, the solution of (1) blows up in finite time.
Proof. 
We take the similar strategy as Theorem 1 in [20] and Theorem 3.7 in [22], using the Kaplan’s first eigenvalue method [23].
Let
B m = x R n ; | x x m | < 1 2 m
for a sequence { x m } m = 1 satisfying | x m | = m for any m N .
Remark 1. 
The method using the sequence of balls B m in (14) was used in [20,24,25,26,27,28].
Let λ m > 0 denote the principal eigenvalue of Δ with Dirichlet problem in B m , and let ϕ m ( x ) denote the corresponding eigenfunction, normalized by
B m ϕ m ( x ) d x = 1 .
Let T ( 0 , T ) be arbitrarily fixed. We define
F m ( t ) = B m u ( x , t ) ϕ m ( x ) d x .
It follows from (1) that
t α B m u ( x , t ) ϕ m ( x ) d x = B m Δ u ( x , t ) ϕ m ( x ) d x + B m u ( x , t ) p ϕ m ( x ) d x f o r t ( 0 , T ) ,
supplemented with the initial condition
F m ( 0 ) = B m u 0 ( x ) ϕ m ( x ) d x .
By integrating by parts, and the fact that ϕ m ( x ) = 0 and ϕ m / ν 0 on B m , where ν denotes the outward unit normal vector to B m at x B m , and applying Green’s formula, we have
t α B m u ( x , t ) ϕ m ( x ) d x B m u ( x , t ) Δ ϕ m ( x ) d x + B m u ( x , t ) p ϕ m ( x ) d x .
Since the principal eigenvalue λ m > 0 and the eigenfunction ϕ m ( x ) satisfy
Δ ϕ m ( x ) = λ m ϕ m ( x ) ,
we obtain
t α B m u ( x , t ) ϕ m ( x ) d x λ m B m u ( x , t ) ϕ m ( x ) d x + B m u ( x , t ) p ϕ m ( x ) d x ,
By (15), (16) and Hölder’ s inequality, we have
F m ( t ) = B m u ( x , t ) ϕ m ( x ) 1 p ϕ m ( x ) 1 1 p d x B m u ( x , t ) p ϕ m ( x ) d x 1 / p .
So, we obtain
B m u ( x , t ) p ϕ m ( x ) d x F m ( t ) p .
Using (16) and (19) in (18), it yields
t α F m ( t ) λ m F m ( t ) + F m ( t ) p f o r t ( 0 , T ) .
Since B m is an n-dimensional ball of radius 1 2 m , it follows that λ m satisfies
λ m c 1 m 2 ,
where c 1 > 0 depends only on the dimension n. Thus, we have
t α F m ( t ) c 1 m 2 F m ( t ) + F m ( t ) p .
Setting H ( ζ ) : = c 1 m 2 ζ + ζ p , then the function H ( ζ ) is convex in ζ > 0 since H C 2 ( 0 , ) and H ( ζ ) 0 . By (3), writing d d t ( k [ F m F m ( 0 ) ] ) ( t ) instead of t α F m ( t ) with k ( t ) = t α Γ ( 1 α ) in (21), we obtain
d d t ( k [ F m F m ( 0 ) ] ) ( t ) H ( F m ( t ) ) f o r t ( 0 , T ) .
It is clear that H ( ζ ) > 0 and H ( ζ ) > 0 for all ζ > ( c 1 m 2 ) 1 p 1 .
Suppose now that
F m ( 0 ) > ( c 1 m 2 ) 1 p 1 .
We claim that (22) implies that F m ( t ) > ( c 1 m 2 ) 1 p 1 for all t ( 0 , T ) . (The fact is stated in the proof of Theorem 3.7 in [22].) Knowing that F m ( t ) F m ( 0 ) > ( c 1 m 2 ) 1 p 1 for all t ( 0 , T ) , it follows from (22) that
t α F m ( t ) = d d t ( k [ F m F m ( 0 ) ] ) ( t ) H ( F m ( t ) ) > 0 f o r a l l t ( 0 , T ) .
Therefore the function F m ( t ) satisfying (24) is an upper solution of the problem
t α ζ = H ( ζ ) = c 1 m 2 ζ + ζ p , ζ ( 0 ) = F m ( 0 ) ,
we have by comparison principle F m ( t ) ζ ( t ) (see Theorem 4.10 in [29]).
On the other hand, since H ( 0 ) = 0 , H ( ζ ) > 0 and H ( ζ ) > 0 for all ζ F m ( 0 ) > ( c 1 m 2 ) 1 p 1 . It then follows from Lemma 3.8 in [22] that v ( t ) = w t α Γ ( α + 1 ) is a lower solution for (25), where v ( t ) satisfies
t α v H ( v ) = c 1 m 2 v + v p , v ( 0 ) F m ( 0 ) ,
and w ( t ) solves the ordinary differential equation
d w d t = H ( w ) = c 1 m 2 w + w p , w ( 0 ) = F m ( 0 ) .
By the comparison principle (see Theorem 4.10 in [29]), we obtain ζ ( t ) v ( t ) . Solving the initial value problem (26), we have the solution
w ( t ) = F m ( 0 ) 1 p 1 exp { ( 1 p ) c 1 m 2 t } c 1 m 2 1 1 p exp c 1 m 2 t .
By the comparison principle (see Theorem 4.10 in [29]), we conclude that
F m ( t ) v ( t ) = w t α Γ ( α + 1 ) = F m ( 0 ) 1 p 1 exp { ( 1 p ) c 2 t α } c 2 Γ ( α + 1 ) 1 1 p exp c 2 t α .
with c 2 = c 1 m 2 Γ ( α + 1 ) . Therefore, from (27), we obtain that v ( t ) as
t log 1 c 2 Γ ( α + 1 ) F m ( 0 ) 1 p ( 1 p ) c 2 1 α ,
and that F m ( t ) . This implies that the solution u ( x , t ) blows up in finite time when (23) holds.
As a result of these arguments, we have the following lemma.
Lemma 3. 
Let F m ( t ) be defined by (16). If F m ( 0 ) as in (17) satisfies (23) for some m N , i.e.,
F m ( 0 ) > A m 2 p 1 f o r s o m e m N ,
where A = c 1 1 p 1 with c 1 as in (20), then u ( x , t ) blows up in finite time.
Here, we shall state the rest of the proof for Theorem 1.
Supposing that u ( x , t ) is a nontrivial global solution, we prove by reductio ad absurdum.
By Lemma 3, then it follows that for any m N
F m ( 0 ) A m 2 p 1 .
Then, by (17) and u 0 ( x ) = λ ψ ( x ) 0 , we obtain
λ B m ψ ( x ) ϕ m ( x ) d x A m 2 p 1 .
Here, if we choose λ to be large enough for any m N , then the left-hand side of (28) is larger than the right-hand side of (28). Thus, we arrive at a contradiction. This completes the proof for the condition (a).
Next, if ψ C 0 ( R n ) satisfies (13), then there is a positive constant L such that ψ ( x ) L | x | a for sufficiently large | x | . Then, we have for sufficiently large m
λ L B m | x | a ϕ m ( x ) d x A m 2 p 1 .
By noting that | x | 3 2 m in B m by (14), we obtain
λ L 3 2 m a B m ϕ m ( x ) d x A m 2 p 1 ,
and then by (15), we have
3 2 a λ L m a A m 2 p 1 .
By multiplying both sides of (29) by m a , we obtain
3 2 a λ L A m a 2 p 1 .
Then, if 0 < a < 2 p 1 and m is sufficiently large, then the left-hand side of (30) is larger than the right-hand side of (30). Thus, we arrive at a contradiction. This completes the proof for the condition (b). □

4. Global Existence

In this section, we state the following global existence result.
Theorem 2. 
Let n 1 and 0 < α < 1 . Assume that the initial data u 0 ( x ) = λ ψ ( x ) 0 , where λ > 0 and ψ C 0 ( R n ) . Suppose that p 1 + 2 / n , and that a > 2 / ( p 1 ) and
lim sup | x | | x | a ψ ( x ) < .
Then, the mild solution of (1) exists globally whenever λ > 0 is small enough.
Proof. 
In what follows, by the letter C we denote generic positive constants, and they may have different values also within the same line.
Since ψ C 0 ( R n ) satisfies (31), there is a constant C > 0 such that
ψ ( x ) C ( 1 + | x | ) a f o r a l l x R n .
Let q c = n ( p 1 ) / 2 . First, if p > 1 + 2 / n and a > 2 / ( p 1 ) , then we know a q c > n . Hence, it follows from (32) that
ψ L q c ( R n ) q c = R n | ψ ( x ) | q c d x C R n ( 1 + | x | ) a q c d x C 0 ( 1 + r ) a q c r n 1 d r C 0 ( 1 + r ) n a q c 1 d r C .
Next, if p = 1 + 2 / n and a > 2 / ( p 1 ) , then we have q c = 1 and a > n = 2 / ( p 1 ) . Hence, it follows from (32) that
ψ L 1 ( R n ) = R n | ψ ( x ) | d x C R n ( 1 + | x | ) a d x C 0 ( 1 + r ) a r n 1 d r C 0 ( 1 + r ) n a 1 d r C .
By (33) and (34), if p 1 + 2 / n and a > 2 / ( p 1 ) , then ψ L q c ( R n ) C . Since u 0 ( x ) = λ ψ ( x ) , u 0 L q c ( R n ) is sufficiently small whenever λ > 0 is small enough. Therefore, the mild solution of (1) exists globally by Proposition 3 (b).
This completes the proof. □

5. Conclusions

In the present paper, we analyze a reaction–diffusion equation with a Caputo fractional derivative in time and with initial conditions. By comparing the conclusions of Zhang and Sun [21] (Proposition 3) and the author (Theorems 1 and 2), we see that the Cauchy problem (1) possesses the critical Fujita exponent
p = 1 + 2 n ,
and the second critical exponent
a = 2 p 1 .
Then we may be summarized in Table 1.
In the case a = a , if 1 < p < p then BU1 by Zhang and Sun [21] (Proposition 3 (a)), but if p p then there are few studies. Therefore, we will consider to study in the case a = a and p p in the future.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This manuscript has no associated data.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fujita, H. On the blowing up of solutions of the Cauchy problem for ut=Δu+u1+α. J. Fac. Sci. Univ. Tokyo Sect. IA Math. 1966, 13, 109–124. [Google Scholar]
  2. Hayakawa, K. On nonexistence of global solution of some semilinear parabolic equation. Proc. Japan Acad. 1973, 49, 503–505. [Google Scholar]
  3. Kobayashi, K.; Sirao, T.; Tanaka, H. On the blowing up problem for semilinear heat equations. J. Math. Soc. Japan 1977, 29, 407–424. [Google Scholar] [CrossRef]
  4. Weissler, F. B. Existence and non-existence of global solutions for semilinear heat equation. Israel J. Math. 1981, 6, 29–40. [Google Scholar] [CrossRef]
  5. Lee, T. Y.; Ni, W. M. Global existence, large time behavior and life span on solution of a semilinear parabolic Cauchy problem. Trans. Amer. Math. Soc. 1992, 333, 365–378. [Google Scholar] [CrossRef]
  6. Pinsky, R. G. Existence and nonexistence of global solutions for ut=Δu+a(x)up in Rn. J. Differ. Eqns. 1997, 133, 152–177. [Google Scholar] [CrossRef]
  7. Galaktionov, V. A.; Kurdyumov, S. P.; Mikhailov, A. P.; Samarskii, A. A. Unbounded solutions of the Cauchy problem for the parabolic equation ut=(uα∇u)+uβ. Soviet Phys. Dokl. 1980, 25, 458–459. [Google Scholar]
  8. Qi, Y. W. On the equation ut=Δuα+uβ. Proc. Roy. Soc. Edinburgh Sect. A. 1993, 123, 373–390. [Google Scholar] [CrossRef]
  9. Mochizuki, K.; Mukai, K. Existence and nonexistence of global solution to fast diffusions with source. Meth. Appl. Anal. 1995, 2, 92–102. [Google Scholar] [CrossRef]
  10. Galaktionov, V. A. Blow-up for quasilinear heat equations with critical Fujita’s exponent. Proc. Roy. Soc. Edinburgh Sect. A 1994, 124, 517–525. [Google Scholar] [CrossRef]
  11. Kawanago, T. Existence and behavior of solutions ut=Δ(um)+ul. Adv. Math. 1997, 7, 367–400. [Google Scholar]
  12. Mochizuki, K.; Suzuki, R. Critical exponent and critical blow-up for quasilinear parabolic equations. Israel J. Math. 1997, 7, 141–156. [Google Scholar] [CrossRef]
  13. Mukai, K.; Mochizuki, K.; Huang, Q. Large time behavior and life span for a quasilinear parabolic equation with slowly decaying initial values. Nonlinear Anal. 2000, 39, 33–45. [Google Scholar] [CrossRef]
  14. Guo, J. S.; Guo, Y. J. On a fast diffusion equation with source. Tohoku Math. J. 2001, 53, 571–579. [Google Scholar] [CrossRef]
  15. Qi, Y. W. The critical exponents of parabolic equations and blow-up in Rn. Proc. Roy. Soc. Edinburgh Sect. A. 1998, 128, 123–136. [Google Scholar] [CrossRef]
  16. Suzuki, R. Existence and nonexistence of global solutions of quasilinear parabolic equations. J. Math. Soc. Japan 2002, 54, 747–792. [Google Scholar] [CrossRef]
  17. Winkler, M. A critical exponent in a degenerate parabolic equation. Math. Meth. Appl. Sci. 2002, 25, 911–925. [Google Scholar] [CrossRef]
  18. Li, Y. H.; Mu, C. L. Life span and a new critical exponent for a degenerate parabolic equation. J. Differ. Eqns. 2004, 207, 392–406. [Google Scholar] [CrossRef]
  19. Yang, C.; Zheng, S.; Zhou, S. Critical exponents in a degenerate parabolic equation with weighted source. Applicable Analysis 2013, 92, 814–830. [Google Scholar] [CrossRef]
  20. Igarashi, T. Blow-up and critical exponents in a degenerate parabolic equation with weighted source. Funkcialaj Ekvacioj 2023, 66, 17–34. [Google Scholar] [CrossRef]
  21. Zhang, Q.G.; Sun, H.R. The blow-up and global existence of solutions of Cauchy problems for a time fractional diffusion equation. Topol. Methods Nonlinear Anal. 2015, 46, 69–92. [Google Scholar] [CrossRef]
  22. Ahmad, B.; Alsaedi, A.; Berbiche, M.; Kirane, M. Existence of global solutions and blow-up of solutions for coupled systems of fractional diffusion equations. Electronic J. Diff. Equ. 2020, 110, 1–28. [Google Scholar] [CrossRef]
  23. Kaplan, S. On the growth of solutions of quasi-linear parabolic equations. Commun. Pure Appl. Math. 1963, 16, 305–330. [Google Scholar] [CrossRef]
  24. Giga, Y.; Umeda, N. Blow-up directions at space infinity for solutions of semilinear heat equations. Bol. Soc. Parana. Mat. 2005, 23, 9–28. [Google Scholar]
  25. Shimoj<i>o</i>¯, M. On blow-up phenomenon at space infinity and its locality for semilinear heat equations (Japanese). Master’s thesis, The University of Tokyo, 2005.
  26. Igarashi, T.; Umeda, N. Existence and nonexistence of global solutions in time for a reaction-diffusion system with inhomogeneous terms. Funkcialaj Ekvacioj 2008, 51, 17–37. [Google Scholar] [CrossRef]
  27. Shimoj<i>o</i>¯, M. The global profile of blow-up at space infinity in semilinear heat equations. J. Math. Kyoto Univ. 2008, 48, 339–361. [Google Scholar]
  28. Shimoj<i>o</i>¯, M.; Umeda, N. Blow-up at space infinity for solutions of cooperative reaction-diffusion systems. Funkcialaj Ekvacioj 2011, 54, 315–334. [Google Scholar]
  29. Li, L.; Liu, J-G. A generalized definition of Caputo derivatives and its application to fractional ODEs. SIAM J. Math. Anal. 2018, 50, 2867–2900. [Google Scholar] [CrossRef]
Table 1. Critical Fujita exponent p and second critical exponent a
Table 1. Critical Fujita exponent p and second critical exponent a
a > a 0 < a < a
1 < p < p BU1 [Proposition 3 (a)] BU1 [Theorem 1 (b)]
BU1 for large λ > 0
p p GE2 for small λ > 0 BU1 [Theorem 1 (b)]
[Theorems 1 (a) and 2]
1BU: Blow-up; 2GE: Global existence
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated