Next Article in Journal
Optimal Purchasing Decisions with Supplier Default in Portfolio Procurement
Next Article in Special Issue
Splitting Methods for Semi-Classical Hamiltonian Dynamics of Charge Transfer in Nonlinear Lattices
Previous Article in Journal
Predicting Perceptual Decision-Making Errors Using EEG and Machine Learning
Previous Article in Special Issue
Sawi Transform Based Homotopy Perturbation Method for Solving Shallow Water Wave Equations in Fuzzy Environment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Blow-Up of the Solution for a Semilinear Parabolic Problem with a Mixed Source

Department of Mathematics, Texas A&M University, Texarkana, TX 75503, USA
Submission received: 2 August 2022 / Revised: 28 August 2022 / Accepted: 30 August 2022 / Published: 2 September 2022
(This article belongs to the Special Issue Computational Methods in Nonlinear Analysis and Their Applications)

Abstract

:
A semilinear parabolic equation with the Dirichlet boundary condition is examined. The reaction source is a mixed nonlinear function. This paper investigates the existence and uniqueness of a solution. A sufficient condition for the solution being blow-up in finite time is found.
MSC:
35K35; 35K55; 35K57; 35K58; 35K61

1. Introduction

Let L be the parabolic operator such that L u = u t u x x , T be a positive real number, c 0 , 1 , and δ x denote the Dirac delta function. In this paper, we study blow-up of the solution of the following semilinear parabolic initial-boundary value problem with a mixed source: a nonlinear function u p combined with a concentrated function δ x c u q
L u = u p + δ x c u q in 0 , 1 × 0 , T ,
u ( x , 0 ) = ψ x for x 0 , 1 , u 0 , t = u 1 , t = 0 for t 0 , T ,
where p 1 and q 1 , and ψ ( x ) is a nontrivial, nonnegative, and differentiable function on 0 , 1 such that ψ 0 = ψ 1 = 0 . Throughout this paper, we assume that ψ ( x ) satisfies the inequality below
ψ x + ψ p ( x ) + δ ( x c ) ψ q x 0 on [ 0 , 1 ] .
We use (3) to show that u is a non-decreasing function in t.
Our study of the problem (1) and (2) is motivated by the papers written by Bimpong-Bota et al. [1] and Ortolvea et al. [2]. In these two articles, the authors investigated instabilities of chemical reactions and transport processes inside a bulk system. Some localized sources, such as electrodes or membranes, are immersed inside this bulk system. They examined the existence of steady state solutions of (1). In (1), u p can represent the rate of bulk reaction and δ x c u q characterizes the rate of localized heterogeneous reaction at the location x = c .
When the source contains only a nonlinear reaction function in the form of f ( u ) , (1) describes the thermal combustion in a bounded domain (cf. [3,4]). Let Ω be a bounded convex domain in R n , and Ω ¯ and Ω be its closure and boundary, respectively. Friedman and McLeod [5] studied (1) and (2) with a nonlinear reaction function in the form of
u t = Δ u + f ( u ) , ( x , t ) Ω × ( 0 , T ) , u ( x , 0 ) = ψ x for x Ω ¯ , u x , t = 0 for x , t Ω × ( 0 , T ) .
They proved that the set of blow-up points is a compact subset of Ω . When n = 1 and f ( u ) = u p with p > 1 , Weissler [6] proved that u blows up at a single point if ψ ( x ) is symmetric and sufficiently large. Caffarrelli and Friedman [7] proved that the blow-up set consists of one or two points if ψ ( x ) has at most two local maxima. When n 1 , Mueller and Weissler [8] proved that u blows up at a single point. Bellout [9] investigated blow-up of u with the Robin boundary condition. If f ( u ) = u p ( x ) , Khelghati and Baghaei [10] showed that u blows up in a finite time when ψ ( x ) is sufficiently large.
When the reaction function contains only a concentrated source, (1) and (2) can describe the physical situation of ignition having a highly localized forcing function (cf. [11,12]). Chan and Tian [11] studied the following degenerate parabolic problem
x r u t u x x = a α r + 2 δ ( x b ) f ( u ( x , t ) ) 0 1 x r | u ( x , t ) | d x α in 0 , 1 × 0 , T ,
where a and α are positive constants, r is a non-negative real number, and b ( 0 , 1 ) . They proved that there exists a unique solution, and u blows up at x = b only.
In Section 2, we prove the existence and uniqueness of the continuous solution u. Then, we prove that u t 0 . Because of the Dirac delta function, it brings challenges to study them. Our approach is through examining the equivalent integral equation. As (1) and (2) has two different heat sources, and each possesses different strength. Without the nonlinear source u p , intuitively, we expect that the blow-up can occur at x = c due to the singularity of heat source. Our main interest is to investigate the location of blow-up point of this problem if u p is present. We want to know whether x = c is still a blow-up point. In Section 3, we prove that the solution u blows up in a finite time under the condition: p + q > 2 . Our method is based on condition (3) alone. Further, an upper bound of blow-up time is obtained. Then, we prove that u blows up at x = c . The proof relies on the integral solution.

2. Existence and Uniqueness of the Solution

Let G x , ξ , t s be the Green’s function corresponding to (1) and (2). For any x and ξ belonging to [ 0 , 1 ] , and t and s are variables within the interval , , G x , ξ , t s satisfies the problem below
L G = δ ( t s ) δ ( x ξ ) for s < t , x ( 0 , 1 ) , G x , ξ , t s = 0 for t < s , G 0 , ξ , t s = G 1 , ξ , t s = 0 .
G x , ξ , t s is a positive function in the set { ( x , ξ , t , s ) : x and ξ are in ( 0 , 1 ) , and t > s 0 } . The integral solution of (1) and (2) is represented by
u ( x , t ) = 0 1 G x , ξ , t ψ ξ d ξ + 0 t 0 1 G ( x , ξ , t s ) u p ( ξ , s ) d ξ d s + 0 t G ( x , c , t s ) u q c , s d s .
From Pinsky [13], G ( x , ξ , t s ) is
G ( x , ξ , t s ) = 2 H ( t s ) j = 1 sin j π x sin j π ξ e j 2 π 2 ( t s ) ,
where H ( t s ) is the Heaviside function. The equivalent expression of G ( x , ξ , t s ) is
G ( x , ξ , t s ) = H ( t s ) 4 π ( t s ) j = e ( x ξ 2 j ) 2 4 ( t s ) e ( x + ξ 2 ( j + 1 ) ) 2 4 ( t s ) .
In R n , because of the singularity ( t s ) n / 2 (corresponding to 1 / t s in the above formula for n = 1 ) of G ( x , c , t s ) in a neighborhood of x = c , the integral solution u ( x , t ) in (4) at x = c does not satisfy (1) and (2) except n = 1 (cf. [14]).
In the following, we want to prove that the integral solution given by (4) is a unique continuous solution of (1) and (2). To achieve it, we construct a sequence of solution u m satisfying the recursive equation below
L u m = u m 1 p + δ x c u m 1 q in ( 0 , 1 ) × 0 , T , u m ( x , 0 ) = ψ ( x ) for x [ 0 , 1 ] , u m 0 , t = u m 1 , t = 0 for t 0 , T ,
for m = 1 , 2 , We assume that u 0 x , t = ψ x on [ 0 , 1 ] × [ 0 , T ) . The integral solution of (6) is
u m ( x , t ) = 0 1 G x , ξ , t ψ ξ d ξ + 0 t 0 1 G ( x , ξ , t s ) u m 1 p ( ξ , s ) d ξ d s + 0 t G ( x , c , t s ) u m 1 q c , s d s .
Since G x , ξ , t s > 0 for t > s , u m > 0 in ( 0 , 1 ) × ( 0 , T ) . The following lemma is going to show that u m is a non-decreasing sequence and u m is a continuous function.
Lemma 1.
The sequence of solution u m of (7) is monotone non-decreasing: ψ x u 1 x , t u m x , t , and u m x , t is continuous on [ 0 , 1 ] × 0 , T .
Proof. 
When m = 1 , from (7) we have
u 1 x , t = 0 1 G ( x , ξ , t ) ψ ξ d ξ + 0 t 0 1 G ( x , ξ , t s ) u 0 p ( ξ , s ) d ξ d s + 0 t G ( x , c , t s ) u 0 q ( c , s ) d s .
By (3), we obtain the inequality below
u 1 x , t ψ ( x ) 0 1 G ( x , ξ , t ) ( ψ ξ ψ ( ξ ) ) d ξ + 0 t 0 1 G ( x , ξ , t s ) ( u 0 p ( ξ , s ) ψ p ( ξ ) ) d ξ d s + 0 t G ( x , c , t s ) ( u 0 q ( c , s ) ψ q ( c ) ) ) d s .
Since G ( x , ξ , t s ) > 0 for t > s , u 1 x , t ψ ( x ) on [ 0 , 1 ] × [ 0 , T ) . To prove that u 1 x , t is continuous on [ 0 , 1 ] × 0 , T , we show that all integrals on the right side of (8) are integrable and lim t 0 u 1 x , t = ψ x on [ 0 , 1 ] . Let us prove that 0 t G ( x , c , t s ) u 0 q ( c , s ) d s is a bounded function. By (5) and switching the summation and integration, the integral becomes
0 t G ( x , c , t s ) u 0 q ( c , s ) d s = 2 j = 1 sin j π x sin j π c 0 t e j 2 π 2 ( t s ) ψ q ( c ) d s .
Let k 1 = max x [ 0 , 1 ] ψ ( x ) . We get
0 t G ( x , c , t s ) u 0 q ( c , s ) d s 2 k 1 q π 2 j = 1 1 j 2 ( 1 e j 2 π 2 t ) .
It is noticed that the right side of above inequality tends to 0 when t approaches 0. From (9),
0 t G ( x , c , t s ) u 0 q ( c , s ) d s 2 π 2 k 1 q j = 1 1 j 2 = 1 3 k 1 q .
Then, we prove 0 t 0 1 G ( x , ξ , t s ) u 0 p ( ξ , s ) d ξ d s being bounded above. By (5),
0 t 0 1 G ( x , ξ , t s ) u 0 p ( ξ , s ) d ξ d s 2 0 t 0 1 j = 1 | sin j π ξ | e j 2 π 2 ( t s ) k 1 p d ξ d s .
After interchanging the summation and integration, we obtain
0 t 0 1 G ( x , ξ , t s ) u 0 p ( ξ , s ) d ξ d s 4 π 3 k 1 p j = 1 1 j 3 ( 1 e j 2 π 2 t ) .
It is noticed that the right side of above inequality tends to 0 when t approaches 0. From (10),
0 t 0 1 G ( x , ξ , t s ) u 0 p ( ξ , s ) d ξ d s 4 π 3 k 1 p j = 1 1 j 3 .
Thus, 0 t 0 1 G ( x , ξ , t s ) u 0 p ( ξ , s ) d ξ d s is bounded above on 0 , 1 × 0 , T . By (5) again, we know
0 1 G ( x , ξ , t ) ψ ξ d ξ 2 k 1 j = 1 1 e j 2 π 2 t .
The right side is bounded when t > 0 . It is noticed that when t = 0 ,
0 1 G ( x , ξ , 0 ) ψ ξ d ξ = 0 1 2 j = 1 sin j π x sin j π ξ ψ ( ξ ) d ξ = ψ x .
Therefore, lim t 0 u 1 x , t = ψ x on [ 0 , 1 ] . Hence, u 1 x , t is continuous on 0 , 1 × 0 , T . Assume that this statement is true for m = i , that is, u i ( x , t ) u i 1 ( x , t ) and u i x , t is continuous on 0 , 1 × 0 , T . When m = i + 1 , we have
u i + 1 x , t u i ( x , t ) = 0 1 G ( x , ξ , t ) ( ψ ξ ψ ( ξ ) ) d ξ + 0 t 0 1 G ( x , ξ , t s ) ( u i p ( ξ , s ) u i 1 p ( ξ , s ) ) d ξ d s + 0 t G ( x , c , t s ) ( u i q ( c , s ) u i 1 q ( c , s ) ) ) d s .
A similar computation from above yields u i + 1 x , t u i x , t and u i + 1 x , t is continuous on [ 0 , 1 ] × 0 , T . By the mathematical induction, we conclude that ψ x u 1 x , t u m x , t and u m x , t is continuous on [ 0 , 1 ] × 0 , T for m = 1 , 2 ,
Let h be a positive real number less than T. Our next result is to prove that the sequence u m t is non-negative.
Lemma 2.
u m t 0 on [ 0 , 1 ] × [ 0 , T ) for m = 1 , 2 ,
Proof. 
To achieve this result, we prove that u m ( x , t + h ) u m ( x , t ) on [ 0 , 1 ] × [ 0 , T h ] . Let m = 1 . From (6), we get
L u 1 ( x , t + h ) = u 0 p x , t + h + δ x c u 0 q x , t + h in ( 0 , 1 ) × ( 0 , T h ] , u 1 ( x , h ) = u 1 ( x , h ) for x [ 0 , 1 ] , u 1 0 , t + h = u 1 1 , t + h = 0 for t ( 0 , T h ] .
If (6) is subtracted from the above problem, it gives
L u 1 ( x , t + h ) u 1 ( x , t ) = u 0 p x , t + h u 0 p ( x , t ) + δ ( x c ) u 0 q x , t + h u 0 q ( x , t ) = 0 .
At t = 0 , u 1 ( x , h ) u 1 ( x , 0 ) = u 1 ( x , h ) ψ ( x ) 0 for x [ 0 , 1 ] . u 1 x , t + h u 1 x , t = 0 at x = 0 and 1 for t ( 0 , T h ] . By the maximum principle, u 1 ( x , t + h ) u 1 ( x , t ) on [ 0 , 1 ] × [ 0 , T h ] . As h 0 , we have u 1 t 0 over [ 0 , 1 ] × [ 0 , T ) . Assume that this statement is true for m = i . When m = i + 1 , a similar computation yields
L u i + 1 ( x , t + h ) u i + 1 ( x , t ) = u i p x , t + h u i p ( x , t ) + δ ( x c ) u i q x , t + h u i q ( x , t ) 0 .
At t = 0 , u i + 1 ( x , h ) u i + 1 ( x , 0 ) = u i + 1 ( x , h ) ψ ( x ) 0 for x [ 0 , 1 ] . u i + 1 x , t + h u i + 1 x , t = 0 at x = 0 and 1 for t ( 0 , T h ] . By the maximum principle, u i + 1 ( x , t + h ) u i + 1 ( x , t ) on [ 0 , 1 ] × [ 0 , T h ] . This implies u i + 1 t 0 on [ 0 , 1 ] × [ 0 , T ) . By the mathematical induction, we have u m t 0 on [ 0 , 1 ] × [ 0 , T ) for m = 1 , 2 ,
Let t ^ > 0 and k 2 be a positive constant such that k 2 > k 1 . The following result shows that { u m } is bounded above by k 2 on [ 0 , 1 ] × [ 0 , t ^ ] .
Lemma 3.
For any k 2 > k 1 , there exists t ^ > 0 such that k 2 > u m on [ 0 , 1 ] × [ 0 , t ^ ] for m = 1 , 2 ,
Proof. 
From Lemma 1, we know that 0 1 G x , ξ , t ψ ξ d ξ , 0 t 0 1 G ( x , ξ , t s ) u m 1 p ( ξ , s ) d ξ d s , and 0 t G ( x , c , t s ) u m 1 q c , s d s are bounded above for m = 1 , 2 , By (9), (10), and 0 1 G ( x , ξ , 0 ) ψ ξ d ξ = ψ x on [ 0 , 1 ] , there exists t ^ > 0 such that
k 2 > 0 1 G x , ξ , t ^ ψ ξ d ξ + 0 t ^ 0 1 G ( x , ξ , t ^ s ) k 2 p d ξ d s + 0 t ^ G ( x , c , t ^ s ) k 2 q d s .
Let m = 1 . Then, by (7)
u 1 ( x , t ^ ) = 0 1 G x , ξ , t ^ ψ ξ d ξ + 0 t ^ 0 1 G ( x , ξ , t ^ s ) u 0 p ( ξ , s ) d ξ d s + 0 t ^ G ( x , c , t ^ s ) u 0 q c , s d s .
As k 2 > k 1 , the above equation is subtracted from (11) to give
k 2 u 1 ( x , t ^ ) > 0 t ^ 0 1 G ( x , ξ , t ^ s ) [ k 2 p u 0 p ( ξ , s ) ] d ξ d s + 0 t ^ G ( x , c , t ^ s ) [ k 2 q u 0 q c , s ] d s > 0 ,
for x [ 0 , 1 ] . By Lemma 2, u 1 is non-decreasing in t for 0 t t ^ . Therefore, k 2 > u 1 ( x , t ) on [ 0 , 1 ] × [ 0 , t ^ ] . Assume that this statement is true for m = i . That is, k 2 > u i ( x , t ) on [ 0 , 1 ] × [ 0 , t ^ ] . When m = i + 1 , a similar computation yields
k 2 u i + 1 ( x , t ^ ) > 0 t ^ 0 1 G ( x , ξ , t ^ s ) [ k 2 p u i p ( ξ , s ) ] d ξ d s + 0 t ^ G ( x , c , t ^ s ) [ k 2 q u i q c , s ] d s > 0 ,
for x [ 0 , 1 ] . Therefore, k 2 > u i + 1 on [ 0 , 1 ] × [ 0 , t ^ ] . By the mathematical induction, we have k 2 > u m on [ 0 , 1 ] × [ 0 , t ^ ] for m = 1 , 2 ,
Let u x , t = lim m u m ( x , t ) . Let us show that (4) has a continuous solution.
Theorem 1.
The integral Equation (4) has a continuous solution u x , t ψ x on [ 0 , 1 ] × [ 0 , t ^ ] .
Proof. 
To obtain this result, we are going to prove u m to converge uniformly. Let S m + 1 = max [ 0 , 1 ] × [ 0 , t ^ ] [ u m + 1 ( x , t ) u m ( x , t ) ] . By (7),
u m + 1 ( x , t ) u m ( x , t ) = 0 t 0 1 G ( x , ξ , t s ) u m p ( ξ , s ) u m 1 p ( ξ , s ) d ξ d s + 0 t G ( x , c , t s ) u m q c , s u m 1 q c , s d s .
By p 1 and q 1 , the mean value theorem, and Lemma 3, we have
u m p ( ξ , s ) u m 1 p ( ξ , s ) p k 2 p 1 [ u m ( ξ , s ) u m 1 ( ξ , s ) ] ,
and
u m q ( c , s ) u m 1 q ( c , s ) q k 2 q 1 [ u m ( c , s ) u m 1 ( c , s ) ] .
Then, we follow a similar computation of (9) and (10) to get
S m + 1 p k 2 p 1 S m 0 t 0 1 G ( x , ξ , t s ) d ξ d s + q k 2 q 1 S m 0 t G ( x , c , t s ) d s 4 p k 2 p 1 π 3 j = 1 1 j 3 ( 1 e j 2 π 2 t ) + 2 q k 2 q 1 π 2 j = 1 1 j 2 ( 1 e j 2 π 2 t ) S m .
By the mean-value theorem, there exists η 1 between 0 and t such that
j = 1 1 j 3 ( 1 e j 2 π 2 t ) = π 2 j = 1 1 j e j 2 π 2 η 1 t .
Similarly, there exists η 2 between 0 and t such that
j = 1 1 j 2 ( 1 e j 2 π 2 t ) = π 2 j = 1 1 e j 2 π 2 η 2 t .
Thus,
S m + 1 4 p k 2 p 1 π j = 1 1 j e j 2 π 2 η 1 + 2 q k 2 q 1 j = 1 1 e j 2 π 2 η 2 t S m .
Since j = 1 1 / ( j e j 2 π 2 η 1 ) and j = 1 1 / e j 2 π 2 η 2 are convergent, we choose t = σ 1 > 0 sufficiently close to 0 such that
4 p k 2 p 1 π j = 1 1 j e j 2 π 2 η 1 + 2 q k 2 q 1 j = 1 1 e j 2 π 2 η 2 σ 1 < 1 .
Then, we have S m + 1 < S m . Therefore, S m m = 1 is a decreasing sequence. By the Dini’s theorem, u m converges uniformly. As u x , t = lim m u m ( x , t ) , u x , t is continuous on [ 0 , 1 ] × [ 0 , σ 1 ] . When t > σ 1 , we repeat the above process through replacing ψ x by u x , σ 1 to obtain that u x , t is continuous on [ 0 , 1 ] × [ 0 , t ^ ] . By the Monotone Convergence Theorem,
u ( x , t ) = 0 1 G ( x , ξ , t ) ψ ξ d ξ + 0 t 0 1 G ( x , ξ , t s ) u p ( ξ , s ) d ξ d s + 0 t G ( x , c , t s ) u q c , s d s ,
is a solution of (1) and (2) on [ 0 , 1 ] × [ 0 , t ^ ] . Since u m ( x , t ) ψ ( x ) on [ 0 , 1 ] × [ 0 , t ^ ] for m = 0 , 1 , 2 , , we have u ( x , t ) ψ ( x ) on [ 0 , 1 ] × [ 0 , t ^ ] . □
Based on a similar proof of Lemma 2, we show that u t is non-negative.
Lemma 4.
u t 0 on [ 0 , 1 ] × [ 0 , t ^ ] .
Proof. 
By Lemma 2, we have u k x , t + h u k ( x , t ) on [ 0 , 1 ] × [ 0 , T ) for k = 0 , 1 , 2 , When k , we have u ( x , t + h ) u ( x , t ) on [ 0 , 1 ] × [ 0 , T ) . Hence, u t 0 on [ 0 , 1 ] × [ 0 , T ) . □
The result below shows the uniqueness of u.
Theorem 2.
The integral solution u in (4) is unique on [ 0 , 1 ] × [ 0 , t ^ ] .
Proof. 
Suppose that u ^ x , t and u ˜ x , t are two different solutions of (4). u ^ x , t u ˜ x , t satisfies
u ^ x , t u ˜ x , t = 0 t 0 1 G x , ξ , t s u ^ p ( ξ , s ) u ˜ p ( ξ , s ) d ξ d s + 0 t G ( x , c , t s ) u ^ q c , s u ˜ q ( c , s ) d s ,
and | u ^ x , t u ˜ x , t | > 0 for some ( x , t ) [ 0 , 1 ] × [ 0 , t ^ ] . By the mean value theorem and (5), we have the inequality below
| u ^ x , t u ˜ x , t | 2 p max 0 , 1 × 0 , t v 1 p 1 max 0 , 1 × 0 , t u ^ u ˜ j = 1 0 t 0 1 sin j π ξ e j 2 π 2 ( t s ) d ξ d s + 2 q max 0 , 1 × 0 , t v 2 q 1 max 0 , 1 × 0 , t u ^ u ˜ j = 1 0 t e j 2 π 2 ( t s ) d s
where v 1 is between u ^ ( ξ , s ) and u ˜ ( ξ , s ) and v 2 is between u ^ c , s and u ˜ ( c , s ) . Let
M = max { max 0 , 1 × 0 , t v 1 p 1 , max 0 , 1 × 0 , t v 2 q 1 } .
We follow a similar computation of Theorem 1 to obtain
max 0 , 1 × 0 , t u ^ u ˜ 4 p π j = 1 1 j e j 2 π 2 η 3 + 2 q j = 1 1 e j 2 π 2 η 4 M t max 0 , 1 × 0 , t u ^ u ˜ ,
where η 3 and η 4 are between 0 and t. We choose t = σ 2 > 0 sufficiently close to 0 such that
4 p π j = 1 1 j e j 2 π 2 η 3 + 2 q j = 1 1 e j 2 π 2 η 4 M σ 2 < 1 .
Then, we have
max 0 , 1 × 0 , σ 2 u ^ x , t u ˜ x , t < max 0 , 1 × 0 , σ 2 u ^ x , t u ˜ x , t .
This contradiction shows that u ^ x , t = u ˜ x , t in 0 , 1 × 0 , σ 2 . We repeat the above process through replacing ψ x by u ^ x , σ 2 to conclude that u ^ x , t = u ˜ x , t on 0 , 1 × [ 0 , t ^ ] . □
Let t b be the supremum of t ^ such that the integral Equation (4) has a unique continuous solution over [ 0 , 1 ] × [ 0 , t ^ ] . The following result shows that u is unbounded if t b is finite.
Theorem 3.
If t b is finite, then u x , t is unbounded somewhere in [ 0 , 1 ] when t t b .
Proof. 
Suppose that u x , t exists at t b for all x [ 0 , 1 ] . We prove this result through showing that u x , t continues to exist for t > t b . Let us consider the following problem
L w = w p + δ x c w q in 0 , 1 × ( t b , t b + h ) , w ( x , t b ) = u x , t b for x 0 , 1 , w 0 , t = w 1 , t = 0 for t t b , t b + h ,
for some h > 0 . We want to prove that the solution w ( x , t ) exists on 0 , 1 × t b , t b + h . Similar to (6) and (7), let us construct a sequence of solution { w m } on [ 0 , 1 ] × [ t b , t b + h ) . The integral solution of w m is given by
w m ( x , t ) = 0 1 G x , ξ , t u ξ , t b d ξ + t b t 0 1 G ( x , ξ , t s ) w m 1 p ( ξ , s ) d ξ d s + t b t G ( x , c , t s ) w m 1 q c , s d s ,
where m = 1 , 2 , and w 0 x , t = u x , t b on [ 0 , 1 ] × [ t b , t b + h ) . We follow a similar proof of Lemma 1 to obtain a continuous monotone non-decreasing sequence solution w m on [ 0 , 1 ] × [ t b , t b + h ) . Let w = lim m w m . We follow from the proof of Theorems 1 and 2 to show that w is a unique continuous solution on [ 0 , 1 ] × [ t b , t b + h ) . Let us define z ( x , t ) as
z ( x , t ) = u ( x , t ) for [ 0 , 1 ] × [ 0 , t b ) , w ( x , t ) for [ 0 , 1 ] × [ t b , t b + h ) .
Then, z ( x , t ) is a unique continuous solution of (1) and (2) on [ 0 , 1 ] × [ 0 , t b + h ) . It leads to a contradiction. Hence, u x , t is unbounded somewhere in [ 0 , 1 ] when t t b . □

3. Blow-Up of the Solution

In this section, we prove that u blows up in a finite time under the condition (3). Our proof does not rely on the first eigenfunction method [15] which requires ψ ( x ) being large enough in the domain. Instead, we modify the method of Lemma 3.4 of [16] to show the finite time blow-up of u.
In the sequel, let us define T = t b , β be a positive real number depending on p and q such that β > 2 , and F be a function in t defined as
F ( t ) = 1 β 0 1 u β d x .
Clearly, F ( t ) > 0 and F ( t ) > 0 for t 0 . Further, F ( t ) has the following property.
Lemma 5.
If p + q > 2 , then
( F ( t ) ) 2 < β 2 ( β 1 ) F ( t ) F ( t ) ,
for t 0 .
Proof. 
The derivative of F ( t ) is
F ( t ) = 0 1 u β 1 u t d x .
By (1), we get
F ( t ) = 0 1 [ u β 1 u x x + u p + β 1 + δ ( x c ) u q + β 1 ] d x .
The derivative of F ( t ) is
F ( t ) = 0 1 [ ( β 1 ) u β 2 u t u x x + u β 1 u x x t ] d x + 0 1 ( p + β 1 ) u p + β 2 u t d x + ( q + β 1 ) u q + β 2 ( c , t ) u t ( c , t ) .
By (1) and a direct computation, we have
F ( t ) = 0 1 ( β 1 ) u β 2 u t [ u t u p δ ( x c ) u q ] d x + 0 1 u β 1 u x x t d x + 0 1 ( p + β 1 ) u p + β 2 u t d x + ( q + β 1 ) u q + β 2 ( c , t ) u t ( c , t ) = ( β 1 ) 0 1 u β 2 ( u t ) 2 d x + 0 1 u β 1 u x x t d x + p 0 1 u p + β 2 u t d x + q u q + β 2 ( c , t ) u t ( c , t ) .
Using integration by parts and u ( 0 , t ) = u ( 1 , t ) = 0 , it yields
F ( t ) = ( β 1 ) 0 1 u β 2 ( u t ) 2 d x + p 0 1 u p + β 2 u t d x + q u q + β 2 ( c , t ) u t ( c , t ) + ( β 1 ) ( β 2 ) 0 1 u β 3 ( u x ) 2 u t d x + ( β 1 ) 0 1 u β 2 u t u x x d x .
By (1) and a direct computation, it gives
F ( t ) = 2 ( β 1 ) 0 1 u β 2 ( u t ) 2 d x + ( β 1 ) ( β 2 ) 0 1 u β 3 ( u x ) 2 u t d x + [ p ( β 1 ) ] 0 1 u p + β 2 u t d x + [ q ( β 1 ) ] u q + β 2 ( c , t ) u t ( c , t ) .
We want to prove that
F ( t ) > 2 ( β 1 ) 0 1 u β 2 ( u t ) 2 d x
through investigating three different cases: (i) q > p 1 , (ii) p > q 1 , (iii) p = q > 1 .
Case (i): If q > p 1 , then p + q > 2 . We let β = p + 1 when p > 1 . Then, by (14)
F ( t ) = 2 p 0 1 u p 1 ( u t ) 2 d x + p ( p 1 ) 0 1 u p 2 ( u x ) 2 u t d x + ( q p ) u q + p 1 ( c , t ) u t ( c , t ) .
By Lemma 4, u t 0 for t 0 when x [ 0 , 1 ] . Then,
F ( t ) > 2 p 0 1 u p 1 ( u t ) 2 d x .
That is, (15). When p = 1 , we set
β = 2 0 1 u β 1 u t d x + ( q + 1 ) u q + β 2 ( c , t ) u t ( c , t ) 0 1 u β 1 u t d x + u q + β 2 ( c , t ) u t ( c , t )
to obtain (15).
Case (ii): If p > q 1 , we let β = q + 1 when q > 1 . We follow the computation of case (i) to obtain (15). When q = 1 , we set
β = ( p + 1 ) 0 1 u p + β 2 u t d x + 2 u β 1 ( c , t ) u t ( c , t ) 0 1 u p + β 2 u t d x + u β 1 ( c , t ) u t ( c , t )
to establish (15).
Case (iii): If p = q > 1 , we let β = p + 1 . By (16), we obtain
F ( t ) = 2 p 0 1 u p 1 ( u t ) 2 d x + p ( p 1 ) 0 1 u p 2 ( u x ) 2 u t d x > 2 ( β 1 ) 0 1 u β 2 ( u t ) 2 d x .
By above three cases, (13), and the Cauchy-Schwartz inequality, we have
( F ( t ) ) 2 0 1 u β d x 0 1 u β 2 ( u t ) 2 d x < β 2 ( β 1 ) F ( t ) F ( t ) ,
for t 0 . The proof is complete. □
The result below provides an estimate of the upper bound of blow-up time of u.
Theorem 4.
If p + q > 2 , then there exists T given by
T = β ( β 2 ) F ( 0 ) F ( 0 )
such that u blows up somewhere in ( 0 , 1 ) for a finite time τ where τ T .
Proof. 
By a direct computation, we obtain
d d t F ( β 2 ) / β = ( β 2 ) β F 2 ( β 1 ) / β F ,
and
d 2 d t 2 F ( β 2 ) / β = 2 ( β 1 ) ( β 2 ) β 2 F ( 3 β + 2 ) / β ( F ) 2 β 2 ( β 1 ) F F .
By (12) and β > 2 , it yields
d 2 d t 2 F ( β 2 ) / β < 0 .
After integrating this expression over [ 0 , t ] , we have
d F ( β 2 ) / β ( t ) d t d F ( β 2 ) / β ( 0 ) d t < 0 .
Equivalently,
d F ( β 2 ) / β ( t ) d t < ( β 2 ) β F 2 ( β 1 ) / β ( 0 ) F ( 0 ) .
Then, we integrate the above expression over [ 0 , t ] to obtain
F ( β 2 ) / β ( t ) < ( β 2 ) β F 2 ( β 1 ) / β ( 0 ) F ( 0 ) t + F ( β 2 ) / β ( 0 ) .
Since F ( 0 ) > 0 and F ( 0 ) > 0 , the right side of above inequality is a decreasing function in t. Thus, there exists a finite time T such that
( β 2 ) β F 2 ( β 1 ) / β ( 0 ) F ( 0 ) T + F ( β 2 ) / β ( 0 ) = 0 .
Then, we solve for T to get
T = β ( β 2 ) F ( 0 ) F ( 0 ) .
Therefore, there exists τ T such that F ( β 2 ) / β ( τ ) = 0 . This implies F ( τ ) = . Hence, u blows up somewhere in ( 0 , 1 ) at τ . The proof is complete. □
The result below is to show that the solution u blows up at x = c .
Theorem 5.
If u blows up in a finite time, then x = c is a blow-up point.
Proof. 
Suppose that u blows up at x * ( 0 , 1 ) where x * c . Let a 1 and a 2 be positive real numbers such that x * and c [ a 1 , a 2 ] ( 0 , 1 ) . By the integral solution (4) at x = x * , it gives
u ( x * , t ) = 0 1 G ( x * , ξ , t ) ψ ( ξ ) d ξ + 0 t 0 1 G ( x * , ξ , t s ) u p ( ξ , s ) d ξ d s + 0 t G ( x * , c , t s ) u q ( c , s ) d s .
By the mean value theorem, there exist η 5 ( 0 , t ) and ξ ^ ( 0 , 1 ) such that
u ( x * , t ) = 0 1 G ( x * , ξ , t ) ψ ( ξ ) d ξ + G ( x * , ξ ^ , t η 5 ) 0 t 0 1 u p ( ξ , s ) d ξ d s + 0 t G ( x * , c , t s ) u q ( c , s ) d s .
By (5), there exists a positive constant k 3 such that G ( x * , ξ ^ , t η 5 ) k 3 (depends on x * and t). Thus,
u ( x * , t ) 0 1 G ( x * , ξ , t ) ψ ξ d ξ + k 3 0 t 0 1 u p ( ξ , s ) d ξ d s + 0 t G ( x * , c , t s ) u q ( c , s ) d s .
Since x * is a blow-up point, u ( x * , t ) when t T but u ( c , t ) remains bounded above when t T . This implies that 0 t 0 1 u p ( ξ , s ) d ξ d s when t T . Since G x , ξ , t > 0 for t > 0 , there exists a positive constant k 4 (depends on x and t) such that G x , ξ , t k 4 for x [ a 1 , a 2 ] . If we consider (4) again at x = c , it yields
u ( c , t ) k 4 [ 0 t 0 1 u p ( ξ , s ) d ξ d s + 0 t u q ( c , s ) d s ] .
Because 0 t 0 1 u p ( ξ , s ) d ξ d s when t T , we have u ( c , t ) when t T . It leads to a contradiction. Thus, c is a blow-up point. □

4. Conclusions

In this paper, the existence and uniqueness of solution of a semilinear parabolic problem with a mixed reaction source: u t u x x = u p + δ ( x c ) u q is studied. We establish this result through investigating the corresponding integral solution. A sufficient condition for finite time blow-up of u is provided. Then, an upper bound of the blow-up time is obtained. Further, we prove that u blows up at c.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The author thanks the anonymous referees for their suggestions.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Bimpong-Bota, K.; Ortoleva, P.; Ross, J. Far-from-equilibrium phenomena at local sites of reaction. J. Chem. Phys. 1974, 60, 3124–3133. [Google Scholar]
  2. Ortoleva, P.; Ross, J. Local structures in chemical reactions with heterogeneous catalysis. J. Chem. Phys. 1972, 56, 4397–4452. [Google Scholar] [CrossRef]
  3. Budd, C.J.; Collins, G.J.; Galaktionov, V.A. An asymptotic and numerical description of self-similar blow-up in quasilinear parabolic equations. J. Comput. Appl. Math. 1998, 97, 51–80. [Google Scholar] [CrossRef]
  4. Chan, C.Y.; Kong, P.C. A thermal explosion model. Appl. Math. Comput. 1995, 71, 201–210. [Google Scholar] [CrossRef]
  5. Friedman, A.; McLeod, B. Blow-up positive solutions of semilinear heat equations. Indiana Univ. Math. J. 1985, 34, 425–447. [Google Scholar] [CrossRef]
  6. Weissler, F.B. Single point blow-up for a semilinear initial value problem. J. Differ. Equ. 1984, 55, 204–224. [Google Scholar] [CrossRef]
  7. Caffarrelli, L.A.; Friedman, A. Blow-up of solutions of nonlinear heat equations. J. Math. Anal. Appl. 1988, 129, 409–419. [Google Scholar] [CrossRef]
  8. Mueller, C.E.; Weissler, F.B. Single point blow-up for a general semilinear heat equation. Indiana Univ. Math. J. 1985, 34, 881–913. [Google Scholar] [CrossRef]
  9. Bellout, H. A criterion for blow-up of solutions to semilinear heat equations. SIAM J. Math. Anal. 1987, 18, 722–727. [Google Scholar] [CrossRef]
  10. Khelghati, A.; Baghaei, K. Blow-up in a semilinear parabolic problem with variable source under positive initial energy. Appl. Anal. 2015, 94, 1888–1896. [Google Scholar] [CrossRef]
  11. Chan, C.Y.; Tian, H.Y. Single-point blow-up for a degenerate parabolic problem with a nonlinear source of local and nonlocal features. Appl. Math. Comput. 2003, 145, 371–390. [Google Scholar] [CrossRef]
  12. Deng, K.; Roberts, C. Quenching for a diffusive equation with a concentrated singularity. Differ. Integral Equ. 1997, 10, 369–379. [Google Scholar]
  13. Pinsky, M.A. Partial Differential Equations and Boundary-Value Problems with Applications; American Mathematical Society: Providence, RI, USA, 2011; p. 454. [Google Scholar]
  14. Chadram, J.M.; Yin, H.M. A diffusion equation with localized chemical reactions. Proc. Edinb. Math. Soc. 1993, 37, 101–118. [Google Scholar] [CrossRef]
  15. Kaplan, S. On the growth of solutions of quasilinear parabolic equations. Comm. Pure. Appl. Math. 1963, 16, 305–333. [Google Scholar] [CrossRef]
  16. Chan, W.Y. Blow-up for degenerate nonlinear parabolic problem. AIMS Math. 2019, 4, 1488–1498. [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

Chan, W.Y. Blow-Up of the Solution for a Semilinear Parabolic Problem with a Mixed Source. Mathematics 2022, 10, 3156. https://0-doi-org.brum.beds.ac.uk/10.3390/math10173156

AMA Style

Chan WY. Blow-Up of the Solution for a Semilinear Parabolic Problem with a Mixed Source. Mathematics. 2022; 10(17):3156. https://0-doi-org.brum.beds.ac.uk/10.3390/math10173156

Chicago/Turabian Style

Chan, Wai Yuen. 2022. "Blow-Up of the Solution for a Semilinear Parabolic Problem with a Mixed Source" Mathematics 10, no. 17: 3156. https://0-doi-org.brum.beds.ac.uk/10.3390/math10173156

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