Next Article in Journal
Molecular Descriptors—Spectral Property Relations for Characterizing Molecular Interactions in Binary and Ternary Solutions, Excited State Dipole Moment Estimation
Previous Article in Journal
Innovative Dual-Stage Blind Noise Reduction in Real-World Images Using Multi-Scale Convolutions and Dual Attention Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Blow-Up Analysis for a Reaction–Diffusion System Coupled via Lα-Norm-Type Sources under Positive Boundary Value Conditions

School of Sciences, Nantong University, Nantong 226019, China
*
Author to whom correspondence should be addressed.
Submission received: 30 August 2023 / Revised: 14 October 2023 / Accepted: 24 October 2023 / Published: 15 November 2023

Abstract

:
This article mainly deals with the blow-up properties of nonnegative solutions for a reaction–diffusion system coupled with norm-type sources under positive boundary value conditions. The local existence of a nonnegative solution and the comparison principle are given. The criteria for the global existence or finite time blow-up of the solutions are obtained by constructing new functions and utilizing the super- and -sub-solution method. The results reveal a correlation between the blow-up profiles of the solutions and the size of the domain, as well as the positive boundary value.

1. Introduction

In this study, we consider the following reaction–diffusion system via L α -norm-type sources, for i = 1 , 2 , , k ,        
u i t = Δ u i m i + a i u i α i p i u i + 1 β i q i ,     u k + 1   = u 1   ,     x Ω , t > 0 ,  
subject to positive boundary value conditions
  u i ( x , t ) = ε 0 > 0     ,         x Ω , t > 0 ,
and initial data
u i ( x , 0 ) = u i , 0 ( x ) ,           x Ω ,      
where Ω R N is a bounded domain with smooth boundary Ω , a i > 0 , m i > 1 ,   α i 1 ,   β i 1 , p i 0 ,     q i > 0 , the initial data u i , 0 ( x ) ε 0   ( i = 1 , 2 , , k ) are bounded functions in Ω , and where · α α = Ω | · | α d x .
The diffusion systems like (1)–(3) emerge in numerous applications across the fields of physics, chemistry, and biology; for instance, the study of the flow of a fluid through a homogeneous isotropic rigid porous medium or the investigation of combustion theory (see [1,2,3,4]), as well as the discussion of population models where communication occurs through chemical means (see [5,6,7,8]).
Over the last few decades, extensive research has been conducted on the blow-up properties of solutions to nonlinear parabolic equations with nonlocal sources under homogeneous Dirichlet boundary conditions. For example, in [9], Deng et al. investigated the following degenerate parabolic equation with the initial and boundary conditions
v t = ( v m ) x x + a l l v q d x ,             x ( l , l ) , t > 0 , v ( ± l , t ) = 0 ,                                                   t > 0 , v ( x , 0 ) = v 0 ( x ) ,                                           x [ l , l ] ,
with l > 0 ,   a > 0 and q > m > 1 . The authors proved that (4) has a global solution if a is small enough, while the solution of (4) blows up in finite time if a l l v 0 q d x 1 and λ < ( 2 a l ) m q 1 a l l φ d x , where λ and φ represent the first eigenvalue and the corresponding eigenfunction of the eigenvalue problem φ x x = λ φ ,   x ( l , l ) , φ ( ± l ) = 0 , respectively. Moreover, the blow-up set is the entire interval [ l , l ] . Later, Duan et al. [4] extended problem (4) to a parabolic system with a nonlocal source and established the uniform blow-up profiles of solutions.
In [10], Deng et al. studied the following problem
u t = Δ u m + a v α p , v t = Δ v n + b u β q     ,             x Ω , t > 0   .
with zero homogeneous Dirichlet boundary data. They proved that if p q < m n , the non-negative solutions are global. For the case p q m n , there are both global solutions and blow-up solutions, which depend on the initial data and the size of domain Ω .
Recently, Liu et al. [11] considered the following degenerate parabolic equations
u t = Δ u m + a v α 1 p 1 u α 2 p 2 ,             v t = Δ v n + b u β 1 q 1 v β 2 q 2 ,         x Ω , t > 0
subject to zero Dirichlet conditions. Based on the studies in [10,12], the authors discussed the criteria for determining whether the solutions of (6) would either exist globally or blow up in finite time. Furthermore, they established various forms of uniform blow-up behavior for simultaneous blow-up solutions.
For other related works on parabolic equations with nonlocal source and more interesting results, we refer the readers to [13,14,15,16,17,18,19,20] and the references therein.
However, there is also a lot of literature concerning the global existence and blow-up properties of solutions for parabolic equations under other boundary conditions, including nonlocal boundary conditions, Neumann conditions, Robin or like-Robin boundary conditions, and so on (see [21,22,23,24,25,26,27,28,29,30] and the references therein). Notably, Ling [21] focused on Equation (5) under positive Dirichlet boundary value conditions and pointed out that small diffusion exponents m ,   n or large coupling exponents p ,   q may lead to the blow-up of solutions. Simultaneously, the author demonstrated that the boundary values also play a significant role in determining the occurrence of blow-up.
Motivated by the works mentioned above, we would like to study the influence of the condition (3) on determining the behavior of both the global and blow-up solutions. Denoting D = i = 1 k ( m i p i ) i = 1 k q i , we have the following results.
Theorem 1.
If the solution of (1)–(3) is global, then  m i p i , for every  i { 1 , 2 , , k } .
Theorem 2.
Assuming that for every i { 1 , 2 , , k } m i > p i D > 0 , then the solution of (1)–(3) exists globally.
Theorem 3.
Assuming that for every i { 1 , 2 , , k } m i > p i D = 0 .
 (1)
If a i   ( i = 1 , 2 , , k ) are sufficiently small, then the solution of (1)–(3) exists globally;
 (2)
If the domain Ω is sufficiently small, then the nonnegative solution of (1)–(3) is global;
 (3)
If the domain Ω contains a sufficiently large ball, then the solution of (1)–(3) blows up in finite time provided that u i , 0 ( x ) are positive and continuous in Ω .
Theorem 4.
Assuming that for every i { 1 , 2 , , k } m i p i D < 0 , then the solution of (1)–(3) blows up in finite time provided that initial data u i , 0 ( x ) are large enough.
The rest of the article is organized as follows. In Section 2, we establish the local existence and the comparison principle which will be used for problems (1)–(3). Then, we give the proof of Theorems 1–4 concerning the global existence and blow-up in finite time in Section 3. Finally, some conclusions are summarized in Section 4.

2. Local Existence and Comparison Principle

Since equations of (1) are degenerate, there are usually no classical solutions. Therefore, we may give a definition of a weak solution for problem (1)–(3). For convenience, we denote Q T = Ω × ( 0 , T ) , S T = Ω × ( 0 , T ) and define a class of test functions as
Φ { ψ C ( Q T ¯ ) ; ψ t , Δ ψ C ( Q T ) L 2 ( Q T ) ; ψ 0 , ψ ( x , t ) | S T = 0 }
Definition 1.
A vector function ( u 1 , u 2 , , u k ) is called a sub- (or super-) solution of (1)–(3) in Q T , if the following conditions hold for every i { 1 , 2 , , k } ,
 (1)
u i L ( Q T ) ;
 (2)
u i ( ) ε 0 ( x , t ) S T , and  u i ( x , 0 ) ( ) u i , 0 ( x ) x Ω ;
 (3)
For every t ( 0 , T )   and  ψ i Φ ,
                                      Ω ( u i ( x , t ) ψ i ( x , t ) u i , 0 ( x ) ψ i ( x , 0 ) ) d x ( )                                                       0 t Ω ( u i ψ i s + u i m i Δ ψ i + a i u i α i p i u i + 1 β i q i ψ i ) d x d s     0 t Ω u i m i ψ i ν d σ d s ,           u k + 1   = u 1   .        
Accordingly, we can say that the vector function ( u 1 , u 2 , , u k ) is a weak solution of problem (1)–(3) if it is both a sub-solution and a super-solution. Moreover, we say ( u 1 , u 2 , , u k ) is global, if ( u 1 , u 2 , , u k ) is a solution of (1)–(3) in Q T for some T < .
Next, we give a maximum principle, which is important for proving the local existence of a solution to (1)–(3).
Lemma 1.
Assuming that w i ( x , t ) C 2 , 1 ( Q T ) C ( Q T ¯ )   ( i = 1 , 2 , , k ) and satisfies
{ w i t d i Δ w i c i 1 w i + c i 2 w i + 1 + c i 3 Ω c i 4 w i + 1 d x + c i 5 Ω c i 6 w i d x ,         w k + 1 = w 1   ,   ( x , t ) Q T w i ( x , t ) ε 0 ,     ( x , t ) S T   w i ( x , 0 ) 0 , x Ω
where  c i j ( i = 1 , 2 , , k ; j = 1 , 2 , 3 , 4 , 5 , 6 )  are bounded functions,
c i j ( x , y ) ,   d i ( x , y ) 0   ( i = 1 , 2 , , k ; j = 1 , 2 , , 6 )
in  Q T .
Then,  w i ( x , t ) 0  on  Q T ¯ .
Proof. 
It can be proven by the similar method in [25,26], so we omit it here. □
The local existence of solutions may be proven by the regularization procedure, so we adopt a similar method (see [15]) and consider the following regularized system, for i = 1 , 2 , , k , and n Z + ,
{ u i , n t = Δ f i , n ( u i , n ) + a i g i , n ( u i , n ) α i p i · g i + 1 , n ( u i + 1 , n ) β i q i ,   u k + 1 , n   = u 1 , n , ( x , t ) Q T , u i , n ( x , t ) = ε 0 + 1 / n ,                                                                                                                                               ( x , t ) S T , u i , n ( x , 0 ) = u i , 0 j + 1 / n ,                                                                                                                                                       x Ω .
where u i , 0 j is a smooth approximation of u i , 0 ( x ) with supp u i , 0 j Ω , and
f i , n ( u i , n ) = { u i , n m i ,                                   u i , n 1 / n , ( 1 / n ) m i ,                   u i , n < 1 / n ,         g i , n ( u i , n ) = { u i , n ,                   u i , n 1 / n , 1 / n ,                 u i , n < 1 / n .
By using a discussion similar to that of Theorems A.1–A.4 in [13], it is known that the problem of (8) has a unique classical solution u i , n j C ( Ω ¯ × [ 0 , T j ( n ) ) ) C 2 , 1 ( Ω   × ( 0 , T j ( n ) ) ) for 0 < T j ( n ) < , where T j ( n ) is the maximal existence time. A direct calculation and the classical maximum principle lead to u i , n j ε 0 + 1 / n . Therefore, u i , n j ( i = 1 , 2 , , k ) satisfies
u i , n j t = Δ ( u i , n j ) m i + a i u i , n j α i p i · u i + 1 , n j β i q i ,   u k + 1 , n j = u 1 , n j ,         ( x , t ) Q T j ( n ) ,
with the corresponding initial and boundary conditions. Obviously, passing to the limit j , it follows that
u i , n lim j u i , n j ,
and ( u 1 , n , u 2 , n , , u k , n ) is a weak solution of
u i , n t = Δ u i , n m i + a i u i , n α i p i · u i + 1 , n β i q i ,     u k + 1 , n = u 1 , n ,             ( x , t ) Q T ( n )
with the corresponding initial and boundary conditions on ( 0 , T ( n ) ) , where T ( n )   lim j T j ( n ) is the maximal existence time. Here, a weak solution of (10) is defined similarly to that for problem (1)–(3), only the equalities for u i , (7) may be replaced with
Ω [ u i , n ( x , t ) ψ i ( x , t ) ( u i , 0 ( x ) + 1 n ) ψ i ( x , 0 ) ] d x   = 0 t Ω [ u i , n ψ i s + u i , n m i Δ ψ i + a i u i , n α i p i u i + 1 , n β i q i ψ i ] d x d s ( ε 0 + 1 n ) m i 0 t Ω ψ i ν d σ d s .
Due to u i , n ε 0 + 1 n and Lemma 1, we have Lemmas 2 and 3 which will be used to prove the local existence of solutions. The proof is standard (see [15]).
Lemma 2.
Assuming that w i ( x , t ) C ( Ω ¯ × [ 0 , T j ( n ) ) ) C 2 , 1 ( Ω × ( 0 , T j ( n ) ) ) for 0 < T j ( n ) < , and ( w 1 , w 2 , , w n ) is a sub- (or super-) solution of 9). Then, ( w 1 , w 2 , , w n ) ( ) ( u 1 , n j , u 2 , n j , , u k , n j ) on Ω ¯ × [ 0 , T j ( n ) ) .
Lemma 3.
If n 1 > n 2 , then u i , n 1 j u i , n 2 j ,   ( i = 1 , 2 , , k ) on Ω ¯ × [ 0 , T j ( n 2 ) ) and T j ( n 1 ) > T j ( n 2 ) .
Hence, the limit T * lim n T ( n ) exists and there are the pointwise limits
u i ( x , t ) lim n u i , n ( x , t ) ,   ( i = 1 , 2 , , k ) ,
for any ( x , t ) Ω ¯ × [ 0 , T * ) . In addition, as the convergence of the sequence u i , n is monotone, passage to the limit   n in identities (11) for   ψ i Φ and   t [ 0 , T * ) , the following theorem can be established by monotone and dominated convergence theorems.
Theorem 5.
(Local existence). Assuming that for every i { 1 , 2 , , k } , nonnegative functions u i , 0 L ( Ω ) , there is some T * = T * ( u i , 0 ) > 0 such that there exists a nonnegative weak solution ( u 1 , u 2 , , u k ) of (1)–(3) for each T < T * . Furthermore, either T * = or lim t T * max x Ω ¯ i = 1 k | u i ( x , t ) | = + .
Proposition 1.
(Comparison principle). Let ( u _ 1 , u _ 2 , , u _ k ) and ( u ¯ 1 , u ¯ 2 , , u ¯ k ) be a nonnegative sub-solution and a super-solution of (1)–(3), respectively. If ( u _ 1 , 0 , u _ 2 , 0 , , u _ k , 0 ) ( u ¯ 1 , 0 , u ¯ 2 , 0 , , u ¯ k , 0 ) and either
Ω u _ i α i d x δ > 0 , Ω u _ i + 1 β i d x δ > 0 ,   u _ i + 1 = u _ 1 ,         i = 1 , 2 , , k ,
or
Ω u ¯ i α i d x δ > 0 , Ω u ¯ i + 1 β i d x δ > 0 ,   u ¯ i + 1 = u ¯ 1 ,         i = 1 , 2 , , k ,
 hold. Then,  ( u _ 1 , u _ 2 , , u _ k ) ( u ¯ 1 , u ¯ 2 , , u ¯ k )  on  Q T .
Although the proof is quite standard and similar to that in [8,15], the comparison principle is very important in proving the existence or blow-up of the solution of (1)–(3). Therefore, we sketch the outline for the reader’s convenience.
Proof. 
Subtracting the first inequalities of (7) for u _ i and u ¯ i yields
Ω [ u _ i ( x , t ) u ¯ i ( x , t ) ] ψ i ( x , t ) d x Ω ( u _ i , 0 u ¯ i , 0 ) ψ i ( x , 0 ) d x                                                                           + 0 t Ω ( u _ i u ¯ i ) [ ψ i s + Φ i ( x , s ) Δ ψ i ] d x d s                                                                           + a i 0 t Ω u _ i + 1 β i q i ( Ω ψ i d x ) D i ( s ) G i ( x , s ) ( u _ i u ¯ i ) d x d s                                                                           + a i 0 t Ω u ¯ i α i p i ( Ω ψ i d x ) E i ( s ) H i ( x , s ) ( u _ i + 1 u ¯ i + 1 ) d x d s ,
where
Φ i ( x , s ) = 0 1 m i [ θ u ¯ i + ( 1 θ ) u _ i ] m i 1 d θ
D i ( s ) = 0 1 p i α i [ θ Ω u ¯ i α i d x + ( 1 θ ) Ω u _ i α i d x ] p i α i 1 d θ
G i ( x , s ) = 0 1 α i [ θ u ¯ i + ( 1 θ ) u _ i ] α i 1 d θ
E i ( s ) = 0 1 q i β i [ θ Ω u ¯ i + 1 β i d x + ( 1 θ ) Ω u _ i + 1 β i d x ] q i β i 1 d θ
H i ( x , s ) = 0 1 β i [ θ u ¯ i + 1 + ( 1 θ ) u _ i + 1 ] β i 1 d θ
Since u _ i and u ¯ i are bounded Q T , it follows from m i > 1 ,   α i ,   β i 1 that Φ i ( x , s ) ,   H i ( x , s ) are bounded. Clearly, if p i / α i ,   q i / β i 1 , then D i ( s ) and E i ( s ) are also bounded. On the other hand, if p i / α i ,   q i / β i < 1 , we have D i ( s ) δ p i / α i 1 ,   E i ( s ) δ q i / β i 1 according to the assumptions (12) or (13). Therefore, we may choose some suitable test function ψ i Φ as in [15] (pp. 118–123) to obtain
Ω [ u _ i ( x , t ) u ¯ i + 1 ( x , t ) ] + d x             ψ i Ω ( u _ i , 0 u ¯ i , 0 ) + d x   + C i 0 t Ω [ ( u _ i u ¯ i ) + + ( u _ i + 1 u ¯ i + 1 ) + ] d x d s ,
where U + max { U ( x , t ) , 0 } and C i are bounded constants. Thus,
Ω i = 1 k [ u _ i ( x , t ) u ¯ i ( x , t ) ] + d x                 max 1 i k { ψ i } Ω i = 1 k ( u _ i , 0 u ¯ i , 0 ) + d x + C 0 t Ω i = 1 k ( u _ i u ¯ i ) + d x d s ,
where C = max 1 i k { C i + C i + 1 } ,   C k + 1 = C 1 . By (15), it follows from the Gronwall lemma that u _ i u ¯ i since u _ i , 0 u ¯ i , 0 . □
Based on the above argument, it can be observed that the solutions of (1)–(3) are unique provided that p i α i ,   q i β i . As a result, we have the following corollary.
Corollary 1.
Assume that p i α i ,   q i β i ( i = 1 , 2 , , k ) . Let ( u _ 1 , u _ 2 , , u _ k ) and ( u ¯ 1 , u ¯ 2 , , u ¯ k ) be a sub-solution and a super-solution of (1)–(3), respectively. Then, ( u _ 1 , u _ 2 , , u _ k ) ( u ¯ 1 , u ¯ 2 , , u ¯ k ) on Q T if u _ i , 0 u ¯ i , 0 .
Denote
A = ( m 1 p 1 q 1 0 0 0 m 2 p 2 0 0 0 0 m k 1 p k 1 q k 1 q k 0 0 m k p k ) , l = ( l 1 l 2 l k 1 l k )
By the theory of linear system, we give the following lemma, which will be used later.
Lemma 4.
Assuming that m i > p i and D = 0 , then there exist positive constants l 1 , l 2 , , l k , such that A l = ( 0 , 0 , , 0 ) T .
Proof. 
Due to D = 0 , then r a n k ( A ) < k , where r a n k ( A ) represents the rank of matrix A . It implies that there must exist nonzero solutions to the linear system A l = ( 0 , 0 , , 0 ) T . Without loss of generality, by taking l 1 = θ > 0 , and l j = i = 1 j 1 ( m i p i ) θ q i > 0 , ( j = 2 , 3 , , k ) , the lemma has been proven directly. □

3. Proof of Global Existence and Blow-Up

In this section, we will employ the super- and sub-solution method to prove the Theorems 1–4. By the comparison principle, we only need to construct appropriate super-solutions or sub-solutions for problem (1)–(3). Firstly, we give the proof of Theorem 1, which is a necessary condition for global solutions.
Proof of Theorem 1.
If the inference is not true, we may assume without loss of generality that p 1 > m 1 and consider the following problem
{ w t = Δ w m 1 + a 1 | Ω | q 1 / β 1 ε 0 q 1 w α 1 p 1 ,                   ( x , t ) Q T ,     w ( x , t ) = ε 0 > 0 ,                                                               ( x , t ) S T ,   w ( x , 0 ) = u 1 , 0 ( x ) ,                                                                   x Ω .
By [17], the solution w of (16) blows up in finite time. Since u i ε 0 ( i = 2 , 3 , , k ) , ε 0 q 1 | Ω | q 1 / β 1 u 1 α 1 p 1 u 1 α 1 p 1 u 2 β 1 q 1 , it implies that ( w , ε 0 , , ε 0 ) is the sub-solution of (1)–(3). Thus, by the comparison principle, the solution ( u 1 , u 2 , , u k ) of (1)–(3) blows up in finite time. □
Proof of Theorem 2.
Now, let φ ( x ) be the unique positive solution of the following elliptic problem:
Δ φ ( x ) = 1 ,   x Ω ;   φ ( x ) = 1 ,                       x Ω .
Taking C 0 = max x Ω φ ( x ) , then 1 φ ( x ) C 0 ,             x Ω       .
Define u ¯ 1 , , u ¯ k as the following:
u ¯ 1 = A l 1 φ 1 / m 1 ( x ) ,   u ¯ 2 = A l 2 φ 1 / m 2 ( x ) , , u ¯ k = A l k φ 1 / m k ( x ) .
where l i ( i = 1 , 2 , , k ) and A > 0 are positive constants to be fixed later. Clearly, u ¯ i A l i . For   i = 1 , 2 , , k , a series of computations yields
Δ u ¯ i m i = A m i l i , a i u ¯ i α i p i u ¯ i + 1 β i q i = a i A p i l i + q i l i + 1 φ 1 / m i α i p i φ 1 / m i + 1 β i q i a i A p i l i + q i l i + 1 C 0 p i m i + q i m i + 1 | Ω | p i α i + q i β i ,
where l k + 1 = l 1 ,     m k + 1 = m 1 . Denote
A i = ( a i C 0 p i m i + q i m i + 1 | Ω | p i α i + q i β i ) 1 / ( m i l i p i l i q i l i + 1 ) ,   i = 1 , 2 , , k .
Due to m i > p i , D > 0 , one may find two positive constants l 1 and l k such that
m k p k q k > l 1 l k > i = 1 k 1 q i i = 1 k 1 ( m i p i )
Furthermore, there exist constants l i > 0   ( i = 2 , 3 , , k 1 ) that satisfy
l 1 l 2 > q 1 m 1 p 1 ,   l 2 l 3 > q 2 m 2 p 2 , , l k 1 l k > q k 1 m k 1 p k 1 .
Thus,
p 1 l 1 + q 1 l 2 < m 1 l 1 ,   p 2 l 2 + q 2 l 3 < m 2 l 2 , , p k l k + q k l 1 < m k l k .
By (20) and (22), we choose A sufficiently large to satisfy
A max 1 i k { A i , ε 0 1 / l i , ( u i , 0 ( x ) ) 1 / l i }
Due to u ¯ i t = 0 ( i = 1 , 2 , , k ), we have
u ¯ i t Δ u ¯ i m i a i u ¯ i α i p i u ¯ i + 1 β i q i 0 ,                       ( x , t ) Q T .
On the other hand, for ( x , t ) S T ,
u ¯ i = A l i φ 1 / m i ( x ) ε 0
and
u ¯ i ( x , 0 ) u i , 0 ( x ) , x Ω .
It follows from (24)–(26) that ( u ¯ 1 , , u ¯ k ) is a positive super-solution of (1)–(3). Hence, ( u 1 , u 2 , , u k ) ( u ¯ 1 , , u ¯ k ) by comparison principle, and it implies that ( u 1 ,   u 2 , , u k ) exists globally. □
Proof of Theorem 3.
Case (1). According to Lemma 4, there exist positive constants l i     ( 1 i k ) satisfying
m l 1 = p 1 l 1 + q 1 l 2 ,   , m k 1 l k 1 = p k 1 l k 1 + q k 1 l k , m k l k = p k l k + q k l 1 .
Clearly, if a i ( i = 1 , 2 , , k ) are sufficiently small such that
a i ( C 0 p i m i + q i m i + 1 | Ω | p i α i + q i β i ) 1   ;   m k + 1 = m 1   .
Choosing A max 1 i k { ε 0 1 / l i , ( u i , 0 ( x ) ) 1 / l i } , it follows from (24)–(26) that ( u ¯ 1 , , u ¯ k ) defined by (18) is a positive super-solution of (1)–(3). This means that the solution ( u 1 , , u k ) of (1)–(3) exists globally.
Case (2). Next, we consider the case where the domain Ω is small enough.
Denote by λ 1 the first eigenvalue of the following eigenvalue problem,
Δ ϕ ( x ) = λ ϕ ( x ) ,     x Ω ;     ϕ ( x ) = 0 ,     x Ω
with the corresponding eigenfunction ϕ 1 ( x ) . Then, λ 1 > 0 and ϕ 1 ( x ) > 0 in Ω . It is well known that λ 1 continuously depends on the domain Ω . Then, for an arbitrary constant λ 2 ( 0 , λ 1 ) , we can find a bounded domain Ω 1 Ω such that λ 2 is the first eigenvalue of the following eigenvalue problem
Δ ϕ ( x ) + λ ϕ ( x ) = 0 ,     x Ω 1 ;     ϕ ( x ) = 0 ,     x Ω 1
Let ϕ 2 ( x ) be the corresponding eigenfunction of (30) with λ 2 . Here, ϕ 2 ( x ) can be normalized as ϕ 2 ( x ) = 1 and ϕ 2 ( x ) > 0 in Ω 1 . So, there exists some constant δ 0 > 0 such that ϕ 2 ( x ) > δ 0 on Ω .
Now, we define the functions u ¯ i ( x , t ) , ( i = 1 , 2 , , k ) as follows:
u ¯ i ( x , t ) = M l i ϕ 2 1 / m i ( x ) ,   x Ω , t > 0 ,
where l i satisfies (27), and M > 0 will be determined later.
Note u ¯ i t = 0 . After a direct calculation, we get
Δ u ¯ i m i + a i u ¯ i α i p i u ¯ i + 1 β i q i = M m i l i ϕ 2 ( λ 2 + c i ϕ 2 1 ) M m i l i ϕ 2 ( λ 2 + c i δ 0 1 )
where
c i = a i ϕ 2 1 / m i ( x ) α i p i ϕ 2 1 / m i + 1 ( x ) β i q i > 0 ,   m k + 1 = m 1
By the relationship between the domain and the corresponding first eigenvalues of problem (30) (see [31]), if the domain Ω is sufficiently small such that λ 1 > λ 0 = max 1 i k { c i δ 0 1 } , then we can choose λ 2 , satisfying λ 0 < λ 2 < λ 1 . Thus, (31) implies that
u ¯ i t Δ u ¯ i m i a i u ¯ i α i p i u ¯ i + 1 β i q i ,   x Ω , t > 0 .
By taking M appropriately large to satisfy M max 1 i k { ( δ 0 1 / m i u i , 0 ) 1 / l i } , we have
u ¯ i ( x , t ) ε 0         ,   x Ω , t > 0 .
and
u ¯ i ( x , 0 ) u i , 0 ( x ) ,
Thus, we prove that ( u ¯ 1 , u ¯ 2 , , u ¯ k ) is a super-solution of (1)–(3). By the comparison principle, the nonnegative solution ( u 1 , , u k ) of (1)–(3) is global.
Case (3). Thanks to i = 1 k ( m i p i ) = i = 1 k q i , there exist positive constants l i ( i = 1 , 2 , ,   k ) that are appropriately large such that
m i p i q i = l i + 1 l i ,   and   ( m i 1 ) l i > 1 ,   l i + 1 = l 1 .
Thus,
l i + 1 < m i l i = p i l i + q i l i + 1
To complete the proof, without loss of generality, we can assume that 0 Ω , and B R Ω is an open ball with center at 0 and radius R .
Denoted by λ B R > 0 the first eigenvalue of the following eigenvalue problem
ϕ ( r ) N 1 r ϕ ( r ) = λ ϕ ( r ) ,   r ( 0 , R ) ;   ϕ ( 0 ) = 0 ,   ( R ) = 0
with the corresponding eigenfunction ϕ R ( r ) . Normalized as ϕ R ( r ) > 0 in B R and ϕ R ( 0 ) = max B R ϕ R ( r ) = 1 . Utilizing the scaling property (let τ = r R ) of eigenvalues and eigenfunctions, we obtain that λ B R = λ B 1 R 2 , ϕ R ( r ) = ϕ 1 ( r R ) = ϕ 1 ( τ ) , where λ B 1 and ϕ 1 ( τ ) are the first eigenvalue and the corresponding normalized eigenfunction of the eigenvalue problem in the unit ball B 1 . Thus, max B 1 ϕ 1 = ϕ 1 ( 0 ) = ϕ R ( 0 ) = max B R ϕ R = 1 .
Constructing functions u _ i ( x , t )   ,     ( i = 1 , 2 , , k ) as follows:
u _ i ( x , t ) = ε 0 ( 1 + d ( 1 c t ) ) l i ϕ R l i ( | x | ) ,             ( x , t ) B R × ( 0 , 1 + d c d )
where c , d > 0 are determined later. Clearly, ( u _ 1 , u _ 2 , , u _ k ) blows up in finite time T 1 + d c d .
For i = 1 , 2 , , k , by directly calculating
u _ i t = l i c d ε 0 ( 1 + d ( 1 c t ) ) l i 1 ϕ R l i , Δ u _ i m i = ( 1 + d ( 1 c t ) ) m i l i [ λ B R m i l i ε 0 m i ϕ R m i l i + m i l i ( m i l i 1 ) ε 0 m i ϕ R m i l i 2 ϕ R 2 ] , u _ i α i p i u _ i + 1 β i q i   = M i , i + 1 ε 0 p i + q i ( 1 + d ( 1 c t ) ) p i l i q i l i + 1 , M i , i + 1 = ϕ R l i α i p i ϕ R l i + 1 β i q i ,   l k + 1 = l 1 .
Then, for ( x , t ) Ω × ( 0 , 1 + d c d ) , we have
u _ i t Δ u _ i m i a i u _ i α i p i u _ i + 1 β i q i   d ε 0 ϕ R l i [ 1 + d ( 1 c t ) ] l i + 1 [ c l i ε 0 m i 1 d [ 1 + d ( 1 c t ) ] m i l i l i 1 ( λ B R m i l i + a i M i , i + 1 ε 0 p i + q i m i ) ] .
Letting
λ i , i + 1 = a i M i , i + 1 ε 0 p i + q i m i m i l i
In view of λ B R = λ B 1 R 2 , we may assume that the ball B R is sufficiently large, i.e., the radius R is large enough such that
λ B R < min 1 i k { λ i , i + 1 }
then
δ i , i + 1 = a i M i , i + 1 ε 0 p i + q i m i λ B R m i l i > 0 , ( i = 1 , 2 , , k ) .
Taking c min 1 i k { ε 0 m i 1 δ i , i + 1 l i d ( 1 + d ) m i l i l i 1 } .
Since u i , 0 are a positive and continuous in Ω , we choose d > 0 large enough to satisfy
u _ i ( x , 0 ) u i , 0 ( x ) , x B R
On the other hand,
u _ i ( x , t ) = 0 < ε 0 ,   x B R × ( 0 , 1 + d c d ) .
By (40)–(43), ( u _ 1 , u _ 2 , , u _ k ) is a positive sub-solution of (1)–(3), which blows up in finite time in the ball B R . It implies that the solution of (1)–(3) blows up in the larger domain Ω . □
Proof of Theorem 4.
Due to m i p i , i = 1 k ( m i p i ) < i = 1 k q i for i { 1 , 2 , , k } , there exist positive constants l i ( i = 1 , 2 , , k ) such that
( m 1 p 1 ) l 1 < q 1 l 2 ,   ( m 2 p 2 ) l 2   < q 2 l 3 , ,   ( m k p k ) l k < q k l 1 ,   and   ( m i 1 ) l i > 1 .
For these positive constants l i   , we construct functions u _ i ( x , t ) as follows:
u _ i ( x , t ) = ε 0 ( 1 + d ( 1 c t ) ) l i e l i ϕ 1 ( x ) ,               ( x , t ) Ω × ( 0 , 1 + d c d ) ,
where c , d > 0 can be determined later, and ϕ 1 ( x ) is the first eigenfunction of the problem (29) with the corresponding first eigenvalue λ 1 > 0 . Here, ϕ 1 ( x ) can be normalized as ϕ 1 = 1 and ϕ 1 ( x ) > 0 in Ω .
Clearly, for ( x , t ) Ω × ( 0 , 1 + d c d )   , and i { 1 , 2 , , k } ,
u _ i ( x , t ) ε 0
Moreover, for ( x , t ) Ω × ( 0 , 1 + d c d ) , we have
  u _ i t = ε 0 l i c d ( 1 + d ( 1 c t ) ) l i 1 e l i ϕ 1 ( x ) , Δ u _ i m i = ( 1 + d ( 1 c t ) ) m i l i [ λ 1 m i l i ε 0 m i e m i l i ϕ 1 ( x ) ϕ 1 + ( m i l i ) 2 ε 0 m i e m i l i ϕ 1 ( x ) | ϕ 1 | 2 ] , u _ i α i p i u _ i + 1 β i q i = b i , i + 1 ε 0 p i + q i ( 1 + d ( 1 c t ) ) p i l i q i l i + 1 ,   b i , i + 1 = e l i ϕ 1 α i p i e l i + 1 ϕ 1 β i q i ,   l k + 1 = l 1
Thus, for every i { 1 , 2 , , k } ,
u _ i t Δ u _ i m i a i u _ i α i p i u _ i + 1 β i q i       l i c d ε 0 [ 1 + d ( 1 c t ) ] l i + 1 e l i ϕ 1 ( x ) + λ 1 m i l i ε 0 m i e m i l i ϕ 1 ( x ) ϕ 1 [ 1 + d ( 1 c t ) ] m i l i a i b i , i + 1 ε 0 p i + q i [ 1 + d ( 1 c t ) ] p i l i + q i l i + 1
= d ε 0 e l i ϕ 1 ( x ) [ 1 + d ( 1 c t ) ] l i + 1 [ c l i + λ 1 m i l i ε 0 m i 1 e ( m i 1 ) l i ϕ 1 ( x ) ϕ 1 d [ 1 + d ( 1 c t ) ] m i l i l i 1 a i b i , i + 1 ε 0 p i + q i 1 e l i ϕ 1 ( x ) d [ 1 + d ( 1 c t ) ] p i l i + q i l i + 1 l i 1 ] d ε 0 e l i ϕ 1 ( x ) [ 1 + d ( 1 c t ) ] l i + 1 [ c l i ε 0 m i 1 d e l i [ 1 + d ( 1 c t ) ] m i l i l i 1 ( λ 1 m i l i e m i l i + a i b i , i + 1 ε 0 p i + q i m i ( 1 + d ) p i l i + q i l i + 1 m i l i ) ] .
Letting
d < min 1 i k { 1 d i , i + 1 d i , i + 1 }
where d i , i + 1 = ( λ 1 m i l i e m i l i a i b i , i + 1 ε 0 p i + q i m i ) 1 / ( p i l i + q i l i + 1 m i l i ) .
Thus,
δ i , i + 1 = λ 1 m i l i e m i l i + a i b i , i + 1 ε 0 p i + q i m i ( 1 + d ) p i l i + q i l i + 1 m i l i > 0 .
Taking
c min 1 i k { ε 0 m i 1 δ i , i + 1 l i d e l i ( 1 + d ) m i l i l i 1 }
By (47)–(49), we have
u _ i t Δ u _ i m i a i u _ i α i p i u _ i + 1 β i q i   0 ,           ( x , t ) Ω × ( 0 , 1 + d c d ) .
Assuming that u i , 0 ( x ) are properly large such that
u _ i ( x , 0 ) u i , 0 ( x ) , x Ω .
Then, it follows from (47), (50), and (51) that ( u _ 1 , u _ 2 , , u _ k ) is a positive sub-solution of (1)–(3). Due to ( u _ 1 , u _ 2 , , u _ k ) , it blows up in time T 1 + d c d , and so does the solution ( u 1 , , u k ) of (1)–(3). □

4. Conclusions

In this article, we discuss the existence and blow-up properties of a diffusion system with positive boundary value conditions. As the above arguments, the positive boundary value conditions guarantee the local existence of solutions of (1)–(3), which also influence the occurrence of the blow-up phenomenon. In the meantime, the blow-up profiles are also affected by the interactions among the multi-nonlinearities m i ,   p i ,   q i . We derive a critical case i = 1 k ( m i p i ) = i = 1 k q i , which belongs to the scenario of global existence (or blow-up) under other assumptions, such as small a i or small (or large) size of the domain Ω. Comparing the results of this article with those under homogeneous Dirichlet boundary conditions, as in [11] (with k = 2 in Equation (1)), we observe differences in the blow-up profiles of the solutions when considering positive value boundary conditions. In this case, the application of specialized processing skills and the construction of new auxiliary functions become necessary to overcome the difficulties caused by boundary values. It is worth noting that the global existence of the two situations appears to be similar, as stated in Theorems 2 and 3(1). The global solutions of (1)–(3) may be independent of the boundary value ε 0 . However, Theorems 3(3) and 4 show that the blow-up situation is slightly different, where ε 0 plays a significant role.
In fact, it follows from the proof of Theorem 3(3) that the global nonexistence of solutions mainly depends on the relation between the boundary value ε 0 and λ B R , where λ B R is the characteristic of the size of B R . That is to say, if we fix the boundary value ε 0 , that λ B R should be properly small (i.e., B R is large enough) such that
λ B R < min 1 i k { a i M i , i + 1 ε 0 p i + q i m i ( m i l i ) 1 }
Correspondingly, for fixed Ω and p i + q i > m i , ε 0 should be suitably large to satisfy
ε 0 max 1 i k { [ λ B R m i l i ( a i M i , i + 1 ) 1 ] 1 / ( p i + q i m i ) }
Similarly, we can observe the relation according to the inequality (49).
As a special case, when the domain Ω is symmetric, particularly when it takes the form of a ball (or an interval), the results of this paper also hold. However, the following problem is still unsolved: how to construct blow-up solutions in an explicit form supporting the theoretical results obtained above for system (1)–(3). Symmetry-based methods (see, e.g., [32,33]) can be useful for its solving.

Author Contributions

Methodology, Q.Z.; Writing—review and editing, G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors express their gratitude to the reviewers and editors for their valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diaz, J.I.; Kersner, R. On a nonlinear degenerate parabolic equation in infiltration or evaporation through a porous medium. J. Differ. Equ. 1987, 69, 368–403. [Google Scholar] [CrossRef]
  2. Vázquez, J.L. The Porous Medium Equation: Mathematical Theory; Oxford University Press: Oxford, UK, 2007. [Google Scholar]
  3. Pao, C.V. Blowing-up of solution for a nonlocal reaction-diffusion problem in combustion theory. J. Math. Anal. App. 1992, 166, 591–600. [Google Scholar] [CrossRef]
  4. Duan, Z.; Deng, W.; Xie, C. Uniform blow-up profile for a degenerate parabolic system with nonlocal source. Comput. Math. Appl. 2004, 47, 977–995. [Google Scholar] [CrossRef]
  5. Meinhardt, H. Models of Biological Pattern Formation; Academic Press: London, UK, 1982. [Google Scholar]
  6. Cantrell, R.S.; Cosner, C. Diffusive logistic equations with indefinite weights: Population models in disrupted environments. Proc. R. Soc. Edinb. Sect. A Math. 1989, 112, 293–318. [Google Scholar] [CrossRef]
  7. Furter, J.; Grinfeld, M. Local vs. non-local interactions in population dynamics. J. Math. Biol. 1989, 27, 65–80. [Google Scholar] [CrossRef]
  8. 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]
  9. Deng, W.; Duan, Z.; Xie, C. The blow-up rate for a degenerate parabolic equation with a non-local source. J. Math. Anal. Appl. 2001, 264, 577–597. [Google Scholar] [CrossRef]
  10. Deng, W.; Li, Y.; Xie, C. Blow-up and global existence for a nonlocal degenerate parabolic system. J. Math. Anal. Appl. 2003, 277, 199–217. [Google Scholar] [CrossRef]
  11. Liu, B.; Zhang, C. Blow-up analysis in degenerate parabolic systems coupled via norm type reactions. Appl. Anal. 2016, 95, 668–689. [Google Scholar] [CrossRef]
  12. Li, M.; Chen, M. Blowup properties for nonlinear degenerate diffusion equations with nonlocal sources. Nonlinear Anal. Real World Appl. 2010, 11, 1122–1130. [Google Scholar] [CrossRef]
  13. Souplet, P. Blow-up in nonlocal reaction-diffusion equations. SIAM J. Math. Anal. 1998, 29, 1301–1334. [Google Scholar] [CrossRef]
  14. Souplet, P. Uniform blow-up profiles and boundary behavior for diffusion equations with nonlocal nonlinear source. J. Differ. Equ. 1999, 153, 374–406. [Google Scholar] [CrossRef]
  15. Anderson, J.R. Local existence and uniqueness of solutions of degenerate parabolic equations. Comm. Part. Differ. Equ. 1991, 16, 105–143. [Google Scholar] [CrossRef]
  16. Du, L. Blow-up for a degenerate reaction–diffusion system with nonlinear nonlocal sources. J. Comput. Appl. Math. 2007, 202, 237–247. [Google Scholar] [CrossRef]
  17. Chen, Y. Blow-up and global existence for a nonlocal degenerate quasilinear parabolic system. Acta Math. Sci. Ser. A 2010, 30, 386–396. [Google Scholar]
  18. Ling, Z.; Wang, Z. Blow-up and global existence for a degenerate parabolic system with nonlocal sources. Discret. Dyn. Nat. Soc. 2012, 2012, 243–248. [Google Scholar]
  19. Zheng, S.; Wang, L. Blow-Up Rate and Profile for a Degenerate Parabolic System Coupled via Nonlocal Sources. Comput. Math. Appl. 2006, 52, 1387–1402. [Google Scholar] [CrossRef]
  20. Huang, S.; Li, F.; Xie, C. Global existence and blow-up of solutions to a nonlocal reaction-diffusion system. Discret. Cont. Dyn. Syst. 2003, 9, 1519–1532. [Google Scholar] [CrossRef]
  21. Ling, Z. Blow-up analysis for a degenerate parabolic System with positive Dirichlet boundary value. Electron. J. Math. Anal. Appl. 2013, 1, 345–354. [Google Scholar]
  22. Lin, Z.; Liu, Y. Uniform blowup profiles for diffusion equations with nonlocal source and nonlocal boundary. Acta Math. Sci. 2004, 24, 443–450. [Google Scholar] [CrossRef]
  23. Cui, Z.; Yang, Z. Roles of weight functions to a nonlinear porous medium equation with nonlocal source and nonlocal boundary condition. J. Math. Anal. Appl. 2008, 342, 559–570. [Google Scholar] [CrossRef]
  24. Li, Y.; Mi, Y.; Mu, C. Properties of positive solutions for a nonlocal nonlinear diffusion equation with nonlocal nonlinear boundary condition. Acta Math. Sci. 2014, 34, 748–758. [Google Scholar] [CrossRef]
  25. Kong, L.; Wang, M. Global existence and blow-up of solutions to a parabolic system with nonlocal sources and boundaries. Sci. China Ser. A Math. 2007, 50, 1251–1266. [Google Scholar] [CrossRef]
  26. Zheng, S.; Kong, L. Roles of weight functions in a nonlinear nonlocal parabolic system. Nonlinear Anal. Theory Methods Appl. 2008, 68, 2406–2416. [Google Scholar] [CrossRef]
  27. Zhong, G.; Tian, L. Blow-up for a parabolic system with nonlocal sources and nonlocal boundary conditions. Bound. Value Probl. 2015, 2015, 61. [Google Scholar] [CrossRef]
  28. Marras, M.; Piro, S.V. Blow-up time estimates in nonlocal reaction-diffusion systems under various boundary conditions. Bound. Value Probl. 2017, 2017, 2. [Google Scholar] [CrossRef]
  29. Liu, B.; Wu, G.; Sun, X.; Li, F. Blow-up estimate in a reaction–diffusion equation with nonlinear nonlocal flux and source. Comptu. Math. Appl. 2019, 78, 1862–1877. [Google Scholar] [CrossRef]
  30. Ma, W.; Yan, B. Global Existence and uniform blow-up to a nonlocal parabolic system with nonlinear boundary conditions arising in a thermal explosion theory. Mathematics 2023, 11, 1993. [Google Scholar] [CrossRef]
  31. Courant, R.; Hilbert, D. Methods of Mathematical Physics; Interscience: New York, NY, USA, 1962. [Google Scholar]
  32. Bluman, G.W.; Cheviakov, A.F.; Anco, S.C. Applications of Symmetry Methods to Partial Differential Equations; Springer: New York, NY, USA, 2010. [Google Scholar]
  33. Cherniha, R.; Serov, M.; Pliukhin, O. Nonlinear Reaction-Diffusion-Convection Equations: Lie and Conditional Symmetry, Exact Solutions and Their Applications; Chapman and Hall/CRC: Boca Raton, FL, USA, 2018. [Google Scholar]
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.

Share and Cite

MDPI and ACS Style

Zhong, G.; Zhang, Q. Blow-Up Analysis for a Reaction–Diffusion System Coupled via Lα-Norm-Type Sources under Positive Boundary Value Conditions. Symmetry 2023, 15, 2074. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15112074

AMA Style

Zhong G, Zhang Q. Blow-Up Analysis for a Reaction–Diffusion System Coupled via Lα-Norm-Type Sources under Positive Boundary Value Conditions. Symmetry. 2023; 15(11):2074. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15112074

Chicago/Turabian Style

Zhong, Guangsheng, and Qinghua Zhang. 2023. "Blow-Up Analysis for a Reaction–Diffusion System Coupled via Lα-Norm-Type Sources under Positive Boundary Value Conditions" Symmetry 15, no. 11: 2074. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15112074

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