Next Article in Journal
Distributed Integrated Synthetic Adaptive Multi-Objective Reactive Power Optimization
Next Article in Special Issue
Analysis of Solutions to a Free Boundary Problem with a Nonlinear Gradient Absorption
Previous Article in Journal
A Residual Voltage Data-Driven Prediction Method for Voltage Sag Based on Data Fusion
Previous Article in Special Issue
First-Order Sign-Invariants and Exact Solutions of the Radially Symmetric Nonlinear Diffusion Equations with Gradient-Dependent Diffusivities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Solutions to a Parabolic System with Absorption

1
School of Sciences, Nantong University, Nantong 226019, China
2
College of Mathematics and Statistics, Nanjing University of Information Science and Technology, Nanjing 210044, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2022, 14(6), 1274; https://doi.org/10.3390/sym14061274
Submission received: 30 May 2022 / Revised: 14 June 2022 / Accepted: 17 June 2022 / Published: 20 June 2022
(This article belongs to the Special Issue Symmetry and Biomathematics: Recent Developments and Challenges)

Abstract

:
In this paper, we investigate the blow-up rate and global existence of solutions to a parabolic system with absorption under the homogeneous Dirichlet boundary. By using the comparison principle and super-sub solution method, we obtain some sufficient conditions for the global existence and blow-up in finite time of solutions and establish some estimates of the upper and lower bounds of the blow-up rates. For the special case, if the domain is symmetric, for example, if it is a ball, the results of this paper also hold.

1. Introduction and Main Results

In this paper, we consider the following problem:
u t Δ u m = v p a u r , x Ω , t > 0 , v t Δ v n = u q b v s , x Ω , t > 0 , u ( x , t ) = 0 , v ( x , t ) = 0 , x Ω , t > 0 , u ( x , 0 ) = u 0 ( x ) , v ( x , 0 ) = v 0 ( x ) , x Ω ,
where m , n > 1 , p , q , r , s , a , b > 0 , Ω R N is a bounded and smooth domain. The initial conditions are u 0 , v 0 C 2 + υ ( Ω ¯ ) , with 0 < υ < 1 , u 0 , v 0 , 0 .
Problem (1) can be described as the cooperative interaction of two diffusing biological species [1]. Suppose that each species can find its habitat from the activities of another species (represented by the reaction terms v p and u q ), corresponding, for example, to the actions of invasive species or to overcrowding (represented by the absorption terms u r and v s ). We also refer to [2,3,4,5]. For the case p , q , r , s > 1 , the existence and uniqueness of non-negative solutions of problem (1) can be obtained using the standard contraction mapping theorem [3,4]. If one of the parameters is less than 1, this method is not directly applicable. However, by using a smooth nonlinear approximation, the local existence can still be easily proved [6]. If p , q > 1 , due to the dissipative property of the absorption terms, the solutions are still unique even if r < 1 or s < 1 .
Problem (1) with a = b = 0 was studied in [7,8] (see also [9,10]), and the following conclusions were drawn:
(a)
If p q < m n , then all solutions of (1) with continuous bounded initial values are global;
(b)
If p q > m n , then there exist both nontrivial global solutions and nonglobal solutions of (1);
(c)
If p q = m n and the diameter of the domain is sufficiently small, all solutions of (1) are global.
It is interesting and important to study the properties of solutions in the critical case p q = m n [9,11]. For other porous equations without absorption, we can consult [12,13].
Problem (1) with m = n = 1 was investigated by Bedjaoui and Souplet in [1]. They divided the ( p , q , r , s ) parameter region into three regions: (i) p q > max { r , 1 } max { s , 1 } ; (ii) p q = max { r , 1 } max { s , 1 } ; and (iii) p q < max { r , 1 } max { s , 1 } . They discussed the parameters of the absorption terms and reaction terms for global solutions and blow-up solutions. Xiang et al. [14] also considered the case, and estimates of the blow-up rates were established when the blow-up case occurred.
In 2021, the authors of [15] considered the following fourth-order parabolic problem with absorption
u t + Δ 2 u = u p , in R N × ( 0 , ) , u ( x , 0 ) = u 0 ( x ) , in R N .
They used the relationships to study the asymptotic behavior of solutions and found that complicated asymptotic behavior can occur in the solutions of this equation with absorption.
For a single parabolic equation
u t = Δ u f ( u ) , x Ω , t > 0
with f ( u ) = λ u α , λ > 0 , α ( 0 , 1 ) , the phenomenon called extinction of the solution is observed. This means that the solutions corresponding to bounded non-negative initial data will vanish identically after a finite time. This was observed first by A. S. Kalashnikov in 1974 (see [16,17] and the references therein). Furthermore, for the case f ( u ) = c ( x , t ) u p , the authors of [18] discussed the parabolic equation with nonlinear nonlocal Neumann boundary conditions, in 2017. They proved a comparison principle and the existence of a local solution, and they also studied the solutions for non-uniqueness and uniqueness.
In 2020, Dao [19] studied the instantaneous shrinking of the support for solutions to parabolic problems.
u t = A ( u ) u β χ u , in R N × ( 0 , T ) , u ( x , 0 ) = u 0 ( x ) , in R N ,
where β ( 0 , 1 ) , χ u is the characteristic function and A ( u ) is a diffusion term, expressing Δ u , ( | u x | p 2 u x ) x or ( u m ) x x in one dimension. Some conditions were given for any solution u of (3) to be compactly supported for all time t > 0 under these three cases. In 2021, Palencia [20] discussed porous-medium problems with advection, and he showed the existence, uniqueness and regularity of weak solutions.
The global solutions and blow-up problems for parabolic systems such as problem (1) have been studied extensively (see [21,22,23,24,25,26,27,28] and the references therein). Motivated by such papers, we discuss two aspects of problem (1). Firstly, we study the optimal conditions for the global existence and blow-up solutions to (1). Secondly, we obtain an precise estimate of the blow-up rate when the blow-up phenomenon occurs, by using the scaled variable.
We start with the notions of a subsolution and supersolution. Throughout this paper, we denote Q T = Ω × ( 0 , T ) .
Definition 1.
We say that a non-negative function ( u , v ) [ C 2 , 1 ( Q T ) C ( Q T ) ] 2 is a supersolution of problem (1) in Q T if
u t Δ u m v p a u r , x Ω , t > 0 , v t Δ v n u q b v s , x Ω , t > 0 , u ( x , t ) 0 , v ( x , t ) 0 , x Ω , t > 0 , u ( x , 0 ) u 0 ( x ) , v ( x , 0 ) v 0 ( x ) , x Ω ,
and ( u , v ) [ C 2 , 1 ( Q T ) C ( Q T ) ] 2 is a subsolution if ( u , v ) 0 satisfies (4) in the reverse order. We say that ( u , v ) is a solution of problem (1) in Q T if ( u , v ) is both a subsolution and a supersolution of (1) in Q T . Furthermore, we say that ( u , v ) is a global solution of problem (1) if it is a solution of (1) in Q T for any T > 0 . And if T < , then any solution ( u , v ) blows up in the sense of the L norm:
lim t T ( u ( t , · ) + v ( t , · ) ) = .
Our main results in detail are as follows.
Theorem 1.
If p q < max { r , m } max { s , n } , then all non-negative solutions of (1) are global. Moreover, if r m , s n (hence p q < r s ), then solutions are uniformly bounded.
Theorem 2.
Let p q = max { r , m } max { s , n } . Then:
(i) 
If r m , s n , and a and b are sufficiently large, then all non-negative solutions are global and uniformly bounded.
(ii) 
If r < m or s < n , then all non-negative solutions are global.
(iii) 
If r m , s n , and a and b are sufficiently small, then there exists a non-negative blow-up solution of (1).
Theorem 3.
If p q > max { r , m } max { s , n } , then any non-negative solution of (1) with large initial data blows up in finite time.
Theorem 4.
Under the assumptions of Theorem 3, assume r p q + 1 p + 1 and s q p + 1 q + 1 . Let ( u , v ) be a non-negative blow-up solution of (1) in finite time T. Then, there exists a positive constant c that satisfies
max Ω ¯ × [ 0 , t ] u ( x , τ ) c ( T t ) α , max Ω ¯ × [ 0 , t ] v ( x , τ ) c ( T t ) β ,
where α = p + 1 p q 1 and β = q + 1 p q 1 . Moreover, if r < p < n 1 , s < q < m 1 , then there exists a positive constant C that satisfies
max Ω ¯ × [ 0 , t ] u ( x , τ ) C ( T t ) α , max Ω ¯ × [ 0 , t ] v ( x , τ ) C ( T t ) β .
Corollary 1.
Suppose that all the hypotheses of Theorem 4 hold. Let ( u , v ) be a solution of (1) that blows up in finite time T. Then, there exist two constants c , C > 0 such that
c ( T t ) α max Ω ¯ × [ 0 , t ] u ( x , τ ) C ( T t ) α , c ( T t ) β max Ω ¯ × [ 0 , t ] v ( x , τ ) C ( T t ) β .
This paper is organized as follows. The comparison principle is established in Section 2. Section 3 is devoted to proving Theorems 1–3. In Section 4, we establish estimates of blow-up rates for blow-up solutions in a finite time, which is Theorem 4. In Section 5, some conclusions and observations are discussed.

2. Comparison Principle

Firstly, we give a comparison principle that plays an important role in the study of problem (1).
Lemma 1.
Suppose that ( u ¯ , v ¯ ) and ( u ̲ , v ̲ ) are a supersolution and a subsolution of (1), respectively. Then, ( u ̲ , v ̲ ) ( u ¯ , v ¯ ) in Ω ¯ × [ 0 , T ) , if ( u ̲ ( x , 0 ) , v ̲ ( x , 0 ) ) ( u ¯ ( x , 0 ) , v ¯ ( x , 0 ) ) and ( u ¯ , v ¯ ) has a positive lower bound.
Proof. 
The technique for the proof of the comparison principle with respect to degenerate parabolic equations is standard [29,30,31]. Here, we will give a brief proof. From the definition of supersolutions and subsolutions, we have
( u ̲ u ¯ ) t ( Δ u ̲ m Δ u ¯ m ) ( v ̲ p v ¯ p ) a ( u ̲ r u ¯ r ) , ( v ̲ v ¯ ) t ( Δ v ̲ n Δ v ¯ n ) ( u ̲ q u ¯ q ) b ( v ̲ s v ¯ s ) .
Define a class of test functions,
Ψ = { ψ C ( Q t ¯ ) | ψ t , Δ ψ C ( Q t ) L 2 ( Q t ) ; ψ 0 ; ψ | Ω = 0 } .
Multiplying by φ ( x , t ) Ψ in the first equation of (7) and integrating over Q t for t ( 0 , T ) , we obtain
Ω ( u ̲ ( x , t ) u ¯ ( x , t ) ) φ ( x , t ) d x Ω ( u ̲ ( x , 0 ) u ¯ ( x , 0 ) ) φ ( x , 0 ) d x + Q t ( u ̲ u ¯ ) φ τ + ( u ̲ m u ¯ m ) Δ φ + ( v ̲ p v ¯ p ) φ a ( u ̲ r u ¯ r ) φ d x d τ , Ω ( u ̲ ( x , 0 ) u ¯ ( x , 0 ) ) φ ( x , 0 ) d x + Q t ( u ̲ u ¯ ) φ τ + Φ ( x , τ ) Δ φ + a | ( u ̲ u ¯ ) G ( x , s ) φ | + ( v ̲ v ¯ ) D ( x , τ ) φ d x d τ ,
where
Φ ( x , s ) = 0 1 m ( θ u ¯ + ( 1 θ ) u ̲ ) m 1 d θ , G ( x , s ) = 0 1 r ( θ u ¯ + ( 1 θ ) u ̲ ) r 1 d θ , D ( x , s ) = 0 1 p ( θ v ¯ + ( 1 θ ) v ̲ ) p 1 d θ .
Since ( u ̲ , v ̲ ) and ( u ¯ , v ¯ ) on Q t are bounded, it follows from m , r , p 1 that Φ ( x , s ) , G ( x , s ) and D ( x , s ) are bounded non-negative functions. Now, if r , p < 1 , then we have
G ( x , s ) δ r 1 , D ( x , s ) δ p 1 .
Thus, an appropriate test function may be chosen exactly as in [29] (pp. 118–123) to obtain
Ω ( u ̲ u ¯ ) + d x φ Ω ( u ̲ ( x , 0 ) u ¯ ( x , 0 ) ) + d x + C Q t ( u ̲ u ¯ ) + + ( v ̲ v ¯ ) + d x d τ ,
where w + = max { w , 0 } and C > 0 is a bounded constant. Similarly, we can prove
Ω ( v ̲ v ¯ ) + d x ψ Ω ( v ̲ ( x , 0 ) v ¯ ( x , 0 ) ) + d x + C Q t ( u ̲ u ¯ ) + + ( v ̲ v ¯ ) + d x d τ ,
where ψ ( x , t ) Ψ . Combining (8) with (9), it follows from the Gronwall lemma that u ̲ u ¯ , v ̲ v ¯ , since u ̲ ( x , 0 ) u ¯ ( x , 0 ) , v ̲ ( x , 0 ) v ¯ ( x , 0 ) .  □

3. Global Existence and Blow-Up

In this section, we will give the proof of global existence and blow-up solutions using the comparison principle. Firstly, we will give the proof of the global existence of Theorem 1.
Proof of Theorem 1.
Suppose that p q < max { r , m } max { s , n } .
Case 1:  m r , n s . This implies p q < m n . By the comparison principle, we have u w and v z , where ( w , z ) satisfies
w t Δ w m = z p , x Ω , t > 0 , z t Δ z n = w q , x Ω , t > 0 , w ( x , t ) = 0 , z ( x , t ) = 0 , x Ω , t > 0 , w ( x , 0 ) = u 0 ( x ) , z ( x , 0 ) = v 0 ( x ) , x Ω ,
It follows from [7,8] that ( w , z ) is global, and so is ( u , v ) .
Case 2:  m r , n s . This implies p q < r s . For any C > 0 , our aim is to find a constant supersolution of (1), i.e., ( u ¯ , v ¯ ) = ( A , B ) satisfies
min { A , B } max { C , max Ω ¯ u 0 ( x ) , max Ω ¯ v 0 ( x ) } .
Since ( u ¯ , v ¯ ) is a supersolution,
a A r B p , b B s A q ,
that is,
B p q a q A q r a q b r B r s .
Since r s > p q , then there is a suitable A , B satisfying (11). By the comparison principle, such inequalities imply that the global solutions are uniformly bounded.
Case 3:  m r , n s . Clearly, we have p q < n r . Let ψ ( x ) be the unique solution of
Δ ψ = 1 , in Ω ; ψ ( x ) = 0 , on Ω .
Now, we will seek a global supersolution of (1) in the following form:
( u ¯ , v ¯ ) = ( h ( t ) + c 0 ) n ψ 1 m , C ( h ( t ) + c 0 ) q ψ 1 n ,
where h ( t ) solves h ( t ) = δ ( h ( t ) + c 0 ) m n n + 1 with 0 < h ( 0 ) = h 0 < 1 , c 0 > 0 and 0 < c 0 + h 0 < 1 , where δ and C satisfy
b ψ s 1 n c 0 max { q ( s 1 ) n ( m 1 ) , 0 } δ q C 1 s , ψ q m C n , δ 1 n ψ 1 m , a c 0 n r p q ψ r m p n C p .
The existence of a suitable δ , C is clear if a, b are sufficiently large.
From (13), a direct computation gives
( h + c 0 ) m n 1 n δ ψ 1 m + ( h + c 0 ) p q C p ψ p n + a ( h + c 0 ) n r p q ψ r m 0 , δ C q ( h + c 0 ) q + m n n ψ 1 n + b C s ( h + c 0 ) q s ψ s n + ( h + c 0 ) n q ( C n ψ q m ) 0 .
Noticing that n r > p q , δ 1 n ψ 1 m and a c 0 n r p q ψ r m p n C p , the first equation of (14) holds. The conditions ψ q m C n and b ψ s 1 n c 0 max { q ( s 1 ) n ( m 1 ) , 0 } δ q C 1 s guarantee the second equation of (14). Indeed, if s < 1 , we only need to require b ψ s 1 n δ q C 1 s . If 1 s n and q ( s 1 ) n ( m 1 ) 0 , we only need to require b ψ s 1 n c 0 q ( s 1 ) n ( m 1 ) δ q C 1 s . If 1 s n and q ( s 1 ) n ( m 1 ) < 0 , we need to require b ψ s 1 n δ q C 1 s .
Case 4:  m r , n s . We see that p q < m s . Exchanging the roles of u and v, we can obtain the result.
Combining these four cases, the results of Theorem 1 hold.  □
Proof of Theorem 2.
Theorem 2 (i) can be proved as in case 2 in Theorem 1.
To prove Theorem 2 (ii), we assume that ( u ¯ , v ¯ ) = ( e α t , c e q α t ) , with c = a 1 / p and α = q 1 a 1 / p . It is obvious that ( u ¯ , v ¯ ) is the supersolution of (1). For more details, see p. 204 in [1], since m 1 , n 1 .
Next, we will prove Theorem 2 (iii). It suffices to establish a blow-up subsolution of problem (1). Without loss of generality, assume that Ω contains the origin. Denote
u ̲ ( x , t ) = ( T t ) α 1 U ( | x | ( T t ) β 1 ) , v ̲ ( x , t ) = ( T t ) α 2 V ( | x | ( T t ) β 2 ) ,
where U ( y ) = ( A 2 y 2 ) + 1 m 1 , V ( y ˜ ) = ( A 2 y ˜ 2 ) + 1 n 1 , β 1 = 1 ( m 1 ) α 1 2 > 0 , β 2 = 1 ( n 1 ) α 2 2 > 0 and α 1 , α 2 , A , T > 0 are to be determined later. Note that B A T β ( 0 ) contains the support of u ̲ ( x , t ) and v ̲ ( x , t ) , where β = max { β 1 , β 2 } if T > 1 ; β = min { β 1 , β 2 } if T 1 , which is included in Ω if T is sufficiently small.
Calculating directly,
u ̲ t = ( T t ) α 1 + 1 ( α 1 U + β 1 y U ) , u ̲ m = ( T t ) ( α 1 m + β 1 ) ( U m ) | x | , u ̲ m = ( T t ) ( α 1 m + 2 β 1 ) ( ( U m ) + N 1 y ( U m ) ) ,
where y = | x | ( T t ) β 1 . We need to find suitable parameters such that
( T t ) ( α 1 + 1 ) α 1 U + β 1 y U ( U m ) N 1 y ( U m ) + a ( T t ) α 1 r U r ( T t ) α 2 p V p
and
( T t ) ( α 2 + 1 ) α 2 V + β 2 y ˜ V ( V n ) N 1 y ˜ ( V n ) + b ( T t ) α 2 s V s ( T t ) α 1 q U q ,
where y ˜ = | x | ( T t ) β 2 .
Note that U , V is continuous for C 2 except for y = A , where U , V has a positive jump. Therefore, to obtain a subsolution of (1), we will prove (15) and (16) pointwise for y > 0 , with y A .
It is easy to see that
U ( y ) = 2 m 1 ( A 2 y 2 ) + 2 m m 1 y , ( U m ) ( y ) = 2 m m 1 ( A 2 y 2 ) + 1 m 1 y , ( U m ) ( y ) = 4 m ( m 1 ) 2 ( A 2 y 2 ) + 2 m m 1 y 2 2 m m 1 ( A 2 y 2 ) + 1 m 1 ,
and (15) is trivial for y A . A simple computation shows that (15) is satisfied. We distinguish two steps for 0 < y θ A and θ A < y < A , where
θ = α 1 + 2 m N m 1 α 1 + 2 m N m 1 + 4 m ( m 1 ) 2 1 2 < 1 .
Step 1. For θ A < y < A ,
α 1 U + β 1 y U ( U m ) N 1 y ( U m ) = ( A 2 y 2 ) + 2 m m 1 [ ( α 1 + 2 m N m 1 ) A 2 ( α 1 + 2 m N m 1 + 2 β 1 m 1 + 4 m ( m 1 ) 2 ) y 2 ] ( A 2 y 2 ) + 2 m m 1 [ ( α 1 + 2 m N m 1 ) A 2 ( α 1 + 2 m N m 1 + 2 β 1 m 1 + 4 m ( m 1 ) 2 ) θ 2 A 2 ] 2 β 1 m 1 θ 2 A 2 ( A 2 y 2 ) + 2 m m 1 ,
Step 2. For 0 < y θ A , we have
U ( y ) ( 1 θ 2 ) 1 m 1 A 2 m 1 > 0 , V ( y ) ( 1 θ 2 ) 1 n 1 A 2 n 1 > 0 .
It follows from α 2 p > α 1 + 1 that
( T t ) ( α 1 + 1 ) [ α 1 U + β 1 y U ( U m ) N 1 y ( U m ) ] 1 2 ( T t ) α 2 p V p ,
if T is sufficiently small. If
a ( T t ) α 1 r U r 1 2 ( T t ) α 2 p V p .
then (15) holds.
Similarly, if
b ( T t ) α 2 s V s 1 2 ( T t ) α 1 q U q ,
α 1 q > α 2 + 1 and T is sufficiently small, then (16) holds.
Next, we choose an appropriate α 1 , α 2 to satisfy (18) and (19). It is clear that there exist α 1 > p + 1 p q 1 and α 2 > q + 1 p q 1 , solving the inequalities
α 2 p > α 1 + 1 , α 1 q > α 2 + 1 .
If p q = r s , then we can choose some large α 1 and α 2 to satisfy (20), and p = α 1 α 2 r , hence q = α 2 α 1 s . Therefore, by (18) and (19), we have
a U r 1 2 V p , b V s 1 2 U q ,
that is,
2 q + r a q b r V r s 2 q a q U r q V p q ,
is clearly true, provided 2 q + r a q b r < 1 , since a and b are sufficiently small.
Therefore, (u, v) is a blow-up subsolution of problem (1) with appropriately large ( u 0 , v 0 ) . This completes the proof of Theorem 2 (iii).  □
Proof of Theorem 3.
Firstly, we consider case (i) r m and s n , hence p q > r n . The discussion for the parameters is similar to that in [1], but here we give more details. In order to guarantee that (21) still holds, we require
α 2 p > α 1 + 1 , α 2 p > α 1 r α 1 q > α 2 + 1 > α 2 s .
Set λ = α 1 / α 2 , then by (22),
s q < λ < p r , s 1 < 1 α 2 < min { p λ , λ q 1 } .
If λ p + 1 q + 1 , then min { p λ , λ q 1 } = λ q 1 . Without loss of generality, we assume that
s q < p + 1 q + 1 .
Indeed, since p q > r s , (23) holds or
r p < q + 1 p + 1 .
If (24) holds, we can just exchange the roles of u and v in problem (1). Thus, we only have to guarantee that (23) holds.
To satisfy (22), we need to find a suitable λ satisfying
s q < λ < min { p + 1 q + 1 , p r } ,
which is possible since s q < p r and α 2 > 0 , such that
0 < s 1 < 1 α 2 < λ q 1 .
Therefore, (21) holds. Thus, ( u ̲ , v ̲ ) is a subsolution of (1) for the case r m and s n .
Next, we consider case (ii) r < m or case (iii) s < n . Assume, for instance, that r < m and s n , hence p q > m s . Any solution of (1) is a supersolution of
u t Δ u m v p a u m a , x Ω , t > 0 , v t Δ v n u q b v s , x Ω , t > 0 ,
with homogeneous Dirichlet boundary conditions. Since p q > m s , similarly to the above proof, we can see that ( u ̲ , v ̲ ) is still a subsolution of (25) for suitable u 0 and v 0 , which implies that ( u , v ) blows up. The cases r m , s < n and r < m , s < n can be treated in a similar way.  □

4. Blow-Up Rate

In this section, we always assume p q > max { r , m } max { s , n } . Fix t ( 0 , T ) such that M u ( t ) = sup Ω ¯ × ( 0 , t ] u ( x , τ ) > 1 , M v ( t ) = sup Ω ¯ × ( 0 , t ] v ( x , τ ) > 1 .
Now, we give a lemma which shows the relationship between M u ( t ) and M v ( t ) near the blow-up time T. The lemma also implies that the blow-up phenomena occur at the same time. The proof is similar to [14], and we omit the details.
Lemma 2.
Assume r p q + 1 p + 1 and s q p + 1 q + 1 . Let ( u , v ) , the solution of (1), blow up in finite time T. Then, there exists η ( 0 , 1 ) such that
η M u 1 2 α M v 1 2 β η 1 , t ( T 2 , T ) .
Proof of Theorem 4.
Denote
φ M u ( y , τ ) = 1 M u ( t ) u ( a 1 y , b 1 τ + t ) , ψ M v ( y , τ ) = 1 M v ( t ) v ( a 2 y , b 2 τ + t ) , in B a ˜ R ( 0 ) × ( 0 , S ) ,
where a 1 = M u m 2 M v p 2 , a 2 = M v n 2 M u q 2 , b 1 = M u M v p , b 2 = M v M u q , a ˜ = max { a 1 1 , a 2 1 } and S = min { M u 1 M v p ( T t ) , M v 1 M u q ( T t ) } .
It is clear that ( φ M u , ψ M v ) blows up at τ = S . Moreover, it is a solution of the following problem:
( φ M u ) τ Δ φ M u m = ψ M v p a M v p M u r φ M u r , ( y , τ ) B a ˜ R ( 0 ) × ( 0 , S ) , = φ M u q b M u q M v s ψ M v s , ( y , τ ) B a ˜ R ( 0 ) × ( 0 , S ) , φ M u ( y , τ ) = 0 , ψ M v ( y , τ ) = 0 , ( y , τ ) B a ˜ R ( 0 ) × ( 0 , S ) , φ M u ( y , 0 ) = 1 M u u ( a 1 y , t ) , ψ M v ( y , 0 ) = 1 M v v ( a 2 y , t ) , y B a ˜ R ( 0 ) .
Step 1: We prove (5). Set
W ¯ ( y , τ ) = ( S 1 s ) α f ( ξ 1 ) , Z ¯ ( y , τ ) = ( S 1 s ) β g ( ξ 2 ) ,
where f ( ξ 1 ) = ( L + δ ( L ξ 1 ) + ) 1 m 1 , g ( ξ 2 ) = ( L + δ ( L ξ 2 ) + ) 1 n 1 , ξ 1 = | y | ( S 1 τ ) β 1 , ξ 2 = | y | ( S 1 τ ) β 2 , δ , L , S 1 to be determined later, and α , β , β 1 , β 2 satisfy
α + 1 = α m + 2 β 1 = β p ,
β + 1 = β n + 2 β 2 = α q .
Moreover, (27) and (28) imply
α = p + 1 p q 1 , β = q + 1 p q 1 , β 1 = 1 ( m 1 ) α 2 , β 2 = 1 ( n 1 ) β 2 .
After a direct computation, for 0 < ξ 1 , ξ 2 < L , we have
W ¯ τ ( y , τ ) = α ( S 1 τ ) ( α + 1 ) ( L + δ ( L ξ 1 ) + ) 1 m 1 δ β 1 m 1 ( S 1 τ ) α + 1 ξ 1 ( L + δ ( L ξ 1 ) + ) 2 m m 1 , Z ¯ τ ( y , τ ) = β ( S 1 τ ) ( β + 1 ) ( L + δ ( L ξ 2 ) + ) 1 n 1 δ β 2 n 1 ( S 1 τ ) β + 1 ξ 2 ( L + δ ( L ξ 2 ) + ) 2 n n 1 , Δ W ¯ m ( y , τ ) = ( S 1 τ ) ( α m + 2 β 1 ) [ m δ 2 ( m 1 ) 2 ( L + δ ( L ξ 1 ) + ) 2 m m 1 m δ m 1 ( L + δ ( L ξ 1 ) + ) 1 m 1 N 1 ξ 1 ] , Δ Z ¯ n ( y , τ ) = ( S 1 τ ) ( β n + 2 β 2 ) [ n δ 2 ( n 1 ) 2 ( L + δ ( L ξ 2 ) + ) 2 n n 1 n δ n 1 ( L + δ ( L ξ 2 ) + ) 1 n 1 N 1 ξ 2 ] .
Then, by (27)
W ¯ s Δ W ¯ m Z ¯ p + a M v p M u r W ¯ r ( S 1 τ ) ( α + 1 ) ( L + δ ( L ξ 1 ) + ) 2 m m 1 · [ α ( L + δ ( L ξ 1 ) + ) δ β 1 m 1 L m δ 2 ( m 1 ) 2 ( L + δ ( L ξ 1 ) + ) m 2 m 1 ( L + δ ( L ξ 2 ) + ) p n 1 ] ( S 1 τ ) ( α + 1 ) ( L + δ ( L ξ 1 ) + ) 2 m m 1 · [ α L δ β 1 m 1 L m δ 2 ( m 1 ) 2 L m 2 m 1 ( c 1 L ) p n 1 ] = ( S 1 τ ) ( α + 1 ) ( L + δ ( L ξ 1 ) + ) 2 m m 1 · [ ( α L 4 δ β 1 m 1 L ) + ( α L 4 m δ 2 ( m 1 ) 2 ) + ( α L 2 L m 2 m 1 + p n 1 c 1 p n 1 ) ] 0 ,
where c 1 = 1 if p < n 1 ; c 1 = 2 if p n 1 , and δ , L satisfy
0 < δ < min { 1 , ( m 1 ) α 4 β 1 , m 1 2 ( α L m ) 1 2 } , 0 < 2 L m 2 m 1 + p n 1 1 c 1 p n 1 < α .
Similarly, for 0 < ξ 2 < L , by (28), we have
Z ¯ s Δ Z ¯ n W ¯ q + b M u q M v s Z ¯ s 0 ,
and δ , L also satisfy
0 < δ < min { 1 , ( n 1 ) β 4 β 2 , n 1 2 ( β L n ) 1 2 } , 0 < 2 L n 2 n 1 + q m 1 1 c 2 q m 1 < β .
where c 2 = 1 if q < m 1 ; c 2 = 2 if q m 1 .
Clearly, W ¯ s Δ W ¯ m Z ¯ p + a M v p M u r W ¯ r 0 , Z ¯ s Δ Z ¯ n W ¯ q + b M u q M v s Z ¯ s 0 , for ξ 1 , ξ 2 > L , W ¯ ( y , s ) > 0 , Z ¯ ( y , s ) > 0 on B a ˜ R ( 0 ) × ( 0 , S 1 ) and W ¯ ( y , 0 ) φ M u ( y , 0 ) , Z ¯ ( y , 0 ) ψ M v ( y , 0 ) , provided that S 1 satisfies, in B a ˜ R ( 0 )
min { S 1 α L 1 m 1 , S 1 β L 1 n 1 } > 1 .
Then, we obtain a supersolution independent of M u , M v . Therefore, the blow-up time of φ M u , ψ M v is more than S 1 , That is, min { M u 1 M v p ( T t ) , M v 1 M u q ( T t ) } S 1 . Without loss of generality, we assume
M u 1 M v p ( T t ) = min { M u 1 M v p ( T t ) , M v 1 M u q ( T t ) } S 1 .
From Lemma 2 and (29), there exists some c > 0 , depending only on S 1 , such that
M u 1 α ( T t ) > c , and , M v 1 β ( T t ) > c .
That is
max Ω ¯ × [ 0 , t ] u ( x , τ ) > c ( T t ) α , max Ω ¯ × [ 0 , t ] v ( x , τ ) > c ( T t ) β .
Step 2: We prove (6). Set
W ̲ ( y , τ ) = ( S 2 s ) α U ( ξ 1 ) , Z ̲ ( y , τ ) = ( S 2 s ) β V ( ξ 2 ) ,
where U ( ξ 1 ) = ( A 2 ξ 1 2 ) 1 m 1 , ξ 1 = | y | ( S 2 τ ) β 1 , V ( ξ 2 ) = ( A 2 ξ 2 2 ) 1 n 1 , ξ 2 = | y | ( S 2 τ ) β 2 , α , β , β 1 , β 2 are given by (29) and S 2 is given later. After some computations, we have
W ̲ s = ( S 2 τ ) ( α + 1 ) ( α U + β 1 ξ 1 U ) , Z ̲ s = ( S 2 τ ) ( β + 1 ) ( β V + β 2 ξ 2 V ) , Δ W ̲ m = ( S 2 τ ) ( α m + 2 β 1 ) ( ( U m ) + N 1 ξ 1 ( U m ) ) , Δ Z ̲ n = ( S 2 τ ) ( β n + 2 β 2 ) ( ( v n ) + N 1 ξ 2 ( V n ) ) .
Thus,
W ̲ s Δ W ̲ m Z ̲ p + a M v p M u r W ̲ r = ( S 2 τ ) ( α + 1 ) ( α U + β 1 ξ 1 U ) ( S 2 τ ) ( α m + 2 β 1 ) ( ( U m ) + N 1 ξ 1 ( U m ) ) ( S 2 τ ) p β V p + a M v p M u r ( S 2 τ ) α r U r ,
Z ̲ s Δ Z ̲ n W ̲ q + b M u q M v s Z ̲ s = ( S 2 τ ) ( β + 1 ) ( β V + β 2 ξ 2 V ) ( S 2 τ ) ( β n + 2 β 2 ) ( ( V n ) + N 1 ξ 2 ( V n ) ) ( S 2 τ ) q α V q + b M u q M v s ( S 2 τ ) β s V s ,
Denote M = max { M u ( t ) , M v ( t ) } , and note (27) and (28) and the assumption r < p < n 1 . Similarly to the proof of Theorem 3 for θ A ξ 1 < A ( θ given by (17)), we have
W ̲ s Δ W ̲ m Z ̲ p + a M v p M u r W ̲ r ( S 2 τ ) ( α + 1 ) [ α U + β 1 ξ 1 U ( U m ) N 1 ξ 1 ( U m ) ] + a M r p ( S 2 τ ) α r U r ( S 2 τ ) ( α + 1 ) ( A 2 ξ 1 2 ) + 2 m m 1 [ 2 β 1 m 1 θ 2 A 2 + a M r p ( A 2 ξ 1 2 ) + m + r 2 m 1 ( S 2 τ ) p β r α ] ( S 2 τ ) ( α + 1 ) A 2 ( A 2 ξ 1 2 ) + 2 m m 1 [ 2 β 1 m 1 θ 2 + a M r p S 2 p β r α A 2 ( r 1 ) m 1 ] ( n o t e : r < p ) 0 ,
where M m 1 2 β θ 2 a S 2 p β r α A 2 ( r 1 ) m 1 .
For 0 < ξ 1 < θ A , by the assumption r < p < n 1 , we have
W ̲ s Δ W ̲ m Z ̲ p + a M v p M u r W ̲ r ( A 2 ξ 1 2 ) + 1 m 1 ( S 2 τ ) ( α + 1 ) [ ( α + 2 m N m 1 ) ( A 2 ξ 2 2 ) + p n 1 ( A 2 ξ 1 2 ) + 1 m 1 + a M r p ( A 2 ξ 1 2 ) + r 1 m 1 ( S 2 τ ) p β r α ] ( A 2 ξ 1 2 ) + 1 m 1 ( S 2 τ ) ( α + 1 ) [ ( α + 2 m N m 1 ) A 2 p n 1 ( 1 θ 2 ) 1 m 1 A 2 m 1 + a M r p c 3 r 1 m 1 A 2 ( r 1 ) m 1 S 2 p β r α ] 0 ,
where c 3 = 1 θ 2 if r < m ; c 3 = 1 if r m . In addition, we require that r < p < n 1 and A and M satisfy
A 2 ( p 1 ) n 1 > ( 1 θ 2 ) 1 m 1 ( α + 2 m N m 1 ) , M [ ( A 2 ( p 1 ) n 1 ( 1 θ 2 ) 1 m 1 ( α + 2 m N m 1 ) ) 1 a c 3 r 1 m 1 S 2 p β r α A 2 ( r 1 ) m 1 ] 1 p r .
Therefore, we combine the two cases θ A ξ 1 < A and 0 < ξ 1 < θ A , and we require that M satisfies
M max { ( m 1 2 β θ 2 a S 2 p β r α A 2 ( r 1 ) m 1 ) , [ ( A 2 ( p 1 ) n 1 ( 1 θ 2 ) 1 m 1 ( α + 2 m N m 1 ) ) 1 a c 3 r 1 m 1 S 2 p β r α A 2 ( r 1 ) m 1 ] 1 p r } .
Thus, we have proved
W ̲ s Δ W ̲ m Z ̲ p + a M v p M u r W ̲ r 0 .
Similarly, we can prove
Z ̲ s Δ Z ̲ n W ̲ q + b M u q M v s Z ̲ s 0 .
Finally, we choose some large initial data u 0 ( x ) and v 0 ( x ) , such that u 0 W ̲ ( y , 0 ) and v 0 Z ̲ ( y , 0 ) . Then, we have a subsolution independent of M. Thus, the blow-up time of φ M u , ψ M v is less than S 2 , i.e.,
max { M u 1 M v p ( T t ) , M v 1 M u q ( T t ) } S 2 .
Without loss of generality, we may assume
M u 1 M v p ( T t ) = max { M u 1 M v p ( T t ) , M v 1 M u q ( T t ) } S 2 .
From Lemma 2 and (29), there exists some C > 0 , depending only on S 2 , such that
M u 1 α ( T t ) < C , and , M v 1 β ( T t ) < C .
That is,
max Ω ¯ × [ 0 , t ] u ( x , τ ) C ( T t ) α , max Ω ¯ × [ 0 , t ] v ( x , τ ) C ( T t ) β .
In view of step 1 and step 2, the conclusions of Theorem 4 hold.  □

5. Conclusions and Observations

In this paper, we presented a novel method for deriving the blow-up rate and global existence of solutions to a parabolic system with absorption. Using the classification relationship between parameters, some interesting results were obtained by constructing appropriate auxiliary functions, using the comparison principle and the super and lower solution method. The three cases p q < max { r , m } max { s , n } , p q = max { r , m } max { s , n } and p q > max { r , m } max { s , n } are discussed in Theorems 1–3, respectively. For the case p q < max { r , m } max { s , n } , all non-negative solutions of (1) are global. For the case p q = max { r , m } max { s , n } , the global solution and blow-up solution of (1) are obtained. For the case p q > max { r , m } max { s , n } , any non-negative solution of (1) with a large initial condition blows up in finite time. However, it is unclear that the solution of (1) with a small initial condition will blow up in finite time or that it has global existence.
We also obtained the result that the solutions u and v blow-up simultaneously when the blow-up case occurs. More importantly, we also analyzed the blow-up rates. For the special case, if the domain is symmetric, for example if it is a ball ( N 2 ) or a segment ( N = 1 ), the results of this paper also hold. In this case, more interesting results may be obtained due to the symmetry of the ball and the segment. The relevant proofs are left to the reader.

Author Contributions

Conceptualization, J.W.; formal analysis, W.L.; writing—review and editing, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by National Natural Science Foundation of China (NO. 11501309).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bedjaoui, N.; Souplet, P. Critical blowup exponents for a system of reaction-diffsion equations with absorption. Z. Angew. Math. Phys. 2002, 53, 197–210. [Google Scholar] [CrossRef]
  2. Bebernes, J.; Eberly, D. Mathematical Problems From Combustion Theory. In Applied Mathematical Sciences; Springer: Berlin, Germany, 1989. [Google Scholar]
  3. Pao, C.V. Nonlinear Parabolic and Elliptic Equations; Plenum Press: New York, NY, USA, 1992. [Google Scholar]
  4. Rothe, F. Global Solutions of Reaction-diffusion Systems, Lecture Notes in Mathematics 1072; Springer: Berlin, Germany, 1984. [Google Scholar]
  5. Smoller, J. Shock Waves and Reaction-Diffusion Equations; Springer: Berlin, Germany, 1983. [Google Scholar]
  6. Escobedo, M.; Herrero, M.A. A semilinear parabolic system in a bounded domain. Ann. Mat. Pura Appl. 1993, 165, 315–336. [Google Scholar] [CrossRef]
  7. Galaktionov, V.A.; Kurdyumov, S.P.; Samarskii, A.A. A parabolic system of quasi-linear equations I. Differ. Equ. 1983, 19, 1558–1571. [Google Scholar]
  8. Galaktionov, V.A.; Kurdyumov, S.P.; Samarskii, A.A. A parabolic system of quasi-linear equations II. Differ. Equ. 1985, 21, 1049–1062. [Google Scholar]
  9. Levine, H.A. The role of critical exponents in blowup theorems. SIAM Rev. 1990, 32, 262–288. [Google Scholar] [CrossRef] [Green Version]
  10. Samarskii, A.A.; Galaktionov, V.A.; Kurdyumov, S.P.; Mikhailov, A.P. Blow-Up in Quasilinear Parabolic Equations; Walter de Gruyter: New York, NY, USA, 1995. [Google Scholar]
  11. Wang, M.X. Note on a quasilinear parabolic system. Nonlinear Anal. 1997, 29, 813–821. [Google Scholar] [CrossRef]
  12. Deng, W.B. Global existence andfinite time blow up for a degenerate reaction-diffusion system. Nonlinear Anal. 2005, 60, 977–991. [Google Scholar] [CrossRef]
  13. Du, W.; Li, Z. Asymptotic analysis for reaction-diffusion equations with absorption. Bound. Value Probl. 2012, 2012, 84. Available online: https://link.springer.com/content/pdf/10.1186/1687-2770-2012-84.pdf (accessed on 25 May 2022). [CrossRef] [Green Version]
  14. Xiang, Z.Y.; Chen, Q.; Mu, C.L. Blow-up rate estimates for a system of reaction-diffusion equations with absorption. J. Korean Math. Soc. 2007, 44, 779–786. [Google Scholar] [CrossRef] [Green Version]
  15. Wu, Y.; Yin, J.; Wang, L. Complicated asymptotic behavior of solutions for the fourth-order parabolic equation with absorption. Appl. Math. Lett. 2021, 2, 107278. [Google Scholar] [CrossRef]
  16. Liu, D.; Yang, L. Extinction phenomenon and decay estimate for a quasilinear parabolic equation with a nonlinear source. Adv. Math. Phys. 2021, 2021, 5569043. [Google Scholar] [CrossRef]
  17. Xu, X.; Cheng, T. Extinction for a p-Laplacian equation with gradient source and nonlinear boundary condition. Appl. Anal. 2021, 2021, 1–16. [Google Scholar] [CrossRef]
  18. Gladkov, A. Initial boundary value problem for a semilinear parabolic equation with absorption and nonlinear nonlocal boundary condition. Lith. Math. J. 2017, 57, 468–478. [Google Scholar] [CrossRef] [Green Version]
  19. Dao, N. Instantaneous shrinking of the support of solutions to parabolic equations with a singular absorption. Rev. Real Acad. Cienc. Exactas FÍSicas Nat. Ser. MatemÁTicas 2020, 114, 165. [Google Scholar] [CrossRef]
  20. Palencia, J. Invasive-invaded system of non-Lipschitz porous medium equations with advection. J. Biomath. Engl. Ed. 2021, 14, 29. [Google Scholar] [CrossRef]
  21. Fila, M.; Quittner, P. The blowup rate for a semilinear parabolic system. J. Math. Anal. Appl. 1999, 238, 468–476. [Google Scholar] [CrossRef] [Green Version]
  22. Gmira, A.; Veron, L. Large time behavior of the solutions of a semilinear parabolic equation in Rn. J. Differ. Equ. 1984, 53, 258–276. [Google Scholar] [CrossRef]
  23. Souplet, P. Blow-up in nonlocal reaction-diffusion equations. SIAM J. Math. Anal. 1998, 29, 1301–1334. [Google Scholar] [CrossRef]
  24. Wang, L.W. The blowup for weakly coupled reaction-diffusion systems. Proc. Am. Math. Soc. 2000, 129, 89–95. [Google Scholar] [CrossRef]
  25. Wang, S.; Wang, M.X.; Xie, C.H. A nonlinear degenerate diffusion equation not in divergence form. Z. Angew. Math. Phys. 2000, 51, 149–159. [Google Scholar] [CrossRef]
  26. Xu, X.H.; Fang, Z.B.; Yi, S. Extinction and decay estimates of solutions for a porous medium equation with nonlocal source and strong absorption. Bound. Value Probl. 2013, 2013, 24. [Google Scholar] [CrossRef] [Green Version]
  27. Zhang, H.; Chen, Y. Extinction for the super diffusion equation with a nonlocal absorption and a gradient source. Appl. Math. Lett. 2016, 61, 163–166. [Google Scholar] [CrossRef]
  28. Zheng, S.N. Global existence andglobal nonexistence of solutions to a reaction-diffusion system. Nonlinear Anal. 2000, 39, 327–340. [Google Scholar] [CrossRef]
  29. Anderson, J.R. Local existence and uniqueness of solutions of degenerate parabolic equations. Commun. Partial Differ. Equ. 1991, 16, 105–143. [Google Scholar] [CrossRef]
  30. Aronson, D.G.; Crandall, M.G.; Peletier, L.A. Stabilization of solutions of a degenerate nonlinear diffusion problem. Nonlinear Anal. 1982, 6, 1001–1022. [Google Scholar] [CrossRef]
  31. Anderson, J.R.; Deng, K. Global existence for degenerate parabolic equations with a non-local forcing. Math. Methods Appl. Sci. 1997, 20, 1069–1087. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lu, H.; Wu, J.; Liu, W. Analysis of Solutions to a Parabolic System with Absorption. Symmetry 2022, 14, 1274. https://doi.org/10.3390/sym14061274

AMA Style

Lu H, Wu J, Liu W. Analysis of Solutions to a Parabolic System with Absorption. Symmetry. 2022; 14(6):1274. https://doi.org/10.3390/sym14061274

Chicago/Turabian Style

Lu, Haihua, Jiayuan Wu, and Wenjun Liu. 2022. "Analysis of Solutions to a Parabolic System with Absorption" Symmetry 14, no. 6: 1274. https://doi.org/10.3390/sym14061274

APA Style

Lu, H., Wu, J., & Liu, W. (2022). Analysis of Solutions to a Parabolic System with Absorption. Symmetry, 14(6), 1274. https://doi.org/10.3390/sym14061274

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop