This paper deals with the properties of positive solutions to a semilinear parabolic system with nonlocal boundary condition. We first give the criteria for finite time blowup or global existence, which shows the important influence of nonlocal boundary. And then we establish the precise blowup rate estimate for small weighted nonlocal boundary.
In this paper, we devote our attention to the singularity analysis of the following semilinear parabolic system:
with nonlocal boundary condition
and initial data
where is a bounded connected domain with smooth boundary and are positive parameters. Most physical settings lead to the default assumption that the functions defined for are nonnegative and continuous, and that the initial data , are nonnegative, which are mathematically convenient and currently followed throughout this paper. We also assume that satisfies the compatibility condition on , and that and for any for the sake of the meaning ofnonlocal boundary.
Over the past few years, a considerable effort has been devoted to studying the blowup properties of solutions to parabolic equations withlocal boundary conditions, say Dirichlet, Neumann, or Robin boundary condition, which can be used to describe heat propagation on the boundary of container (see the survey papers [1, 2]). For example, the system (1.1) and (1.3) with homogeneous Dirichlet boundary condition
(i) All solutions are global if , while there exist both global solutions and finite time blowup solutions depending on the size of initial data when (See ). (ii) The asymptotic behavior near the blowup time is characterized by
On the other hand, there are a number of important phenomena modeled by parabolic equations coupled with nonlocal boundary condition of form (1.2). In this case, the solution could be used to describe the entropy per volume of the material [10–12]. Over the past decades, some basic results such as the global existence and decay property have been obtained for the nonlocal boundary problem (1.1)–(1.3) in the case of scalar equation (see [13–16]). In particular, for the blowup solution of the single equation
under the assumption that , Seo  established the following blowup rate estimate
for any . For the more nonlocal boundary problems, we also mention the recent works [17–22]. In particular, Kong and Wang in , by using some ideas of Souplet , obtained the blowup conditions and blowup profile of the following system:
subject to nonlocal boundary (1.2), and Zheng and Kong in  gave the condition for global existence or nonexistence of solutions to the following similar system:
with nonlocal boundary condition (1.2). The typical characterization of systems (1.8) and (1.9) is the complete couple of the nonlocal sources, which leads to the analysis of simultaneous blowup.
To our surprise, however, it seems that there is no work dealing with singularity analysis of the parabolic system (1.1) with nonlocal boundary condition (1.2) except for the single equation case, although this is a very classical model. Therefore, the basic motivation for the work under consideration was our desire to understand the role of weight function in the blowup properties of that nonlinear system. We first remark by the standard theory [4, 13] that there exist local nonnegative classical solutions to this system.
Our main results read as follows.
Suppose that . All solutions to (1.1)–(1.3) exist globally.
It follows from Theorem 1.2 and Proposition 1.1(i) that any weight perturbation on the boundary has no influence on the global existence when , while the following theorem shows that it plays an important role when . In particular, Theorem 1.3(ii) is completely different from the case of the local boundary (1.4) (by comparing with Proposition 1.1(i)).
Suppose that .
(i) For any nonnegative and , solutions to (1.1)–(1.3) blow up in finite time provided that the initial data are large enough.
(ii) If , for any , then any solutions to (1.1)–(1.3) with positive initial data blow up in finite time.
(iii) If , for any , then solutions to (1.1)–(1.3) with small initial data exist globally in time.
Once we have characterized for which exponents and weights the solution to problem (1.1)–(1.3) can or cannot blow up, we want to study the way the blowing up solutions behave as approaching the blowup time. To this purpose, the first step usually consists in deriving a bound for the blowup rate. For this bound estimate, we will use the classical method initially proposed in Friedman and McLeod . The use of the maximum principle in that process forces us to give the following hypothesis technically.
(H)There exists a constant, such that
However, it seems that such an assumption is necessary to obtain the estimates of type (1.5) or (1.10) unless some additional restrictions on parameters are imposed (for the related problem, we refer to the recent work of Matano and Merle ).
Here to obtain the precise blowup rates, we shall devote to establishing some relationship between the two components and as our problem involves a system, but we encounter the typical difficulties arising from the integral boundary condition. The following theorem shows that we have partially succeeded in this precise blowup characterization.
Suppose that , , , and assumption (H) holds. If the solution of (1.1)–(1.3) with positive initial data blows up in finite time , then
where are both positive constants.
If and , then Theorem 1.4 implies that for the blowup solution of problem (1.6), we have the following precise blowup rate estimate:
which improves the estimate (1.7). Moreover, we relax the restriction on .
By comparing with Proposition 1.1(ii), Theorem 1.4 could be explained as the small perturbation of homogeneous Dirichlet boundary, which leads to the appearance of blowup, does not influence the precise asymptotic behavior of solutions near the blowup time and the blowup rate exponents and are just determined by the corresponding ODE system . Similar phenomena are also noticed in our previous work , where the single porous medium equation is studied.
The rest of this paper is organized as follows. Section 2 is devoted to some preliminaries, which include the comparison principle related to system (1.1)–(1.3). In Section 3, we will study the conditions for the solution to blow up and exist globally and hence prove Theorems 1.2 and 1.3. Proof of Theorem 1.4 is given in Section 4.
In this section, we give some basic preliminaries. For convenience, we denote We begin with the definition of the super- and subsolution of system (1.1)–(1.3).
A pair of functions is called a subsolution of (1.1)–(1.3) if
A supersolution is defined with each inequality reversed.
Suppose that and are nonnegative functions. If satisfy
then on .
Set . Since , by continuity, there exists such that for all . Thus .
We claim that will lead to a contradiction. Indeed, suggests that or for some . Without loss of generality, we suppose that .
If , we first notice that
In addition, it is clear that on boundary and at the initial state . Then it follows from the strong maximum principle that in , which contradicts to .
If , we shall have a contradiction:
In the last inequality, we have used the facts that for any and for any , which is a direct result of the previous case.
Therefore, the claim is true and thus , which implies that on .
If and for any in Lemma 2.2, we can obtain in under the assumption that for . Indeed, for any , we can conclude that in as the proof of Lemma 2.2. Then the desired result follows from the limit procedure .
From the above lemma, we can obtain the following comparison principle by the standard argument.
Let ) and be a subsolution and supersolution of (1.1)–(1.3) in , respectively. If for , then in .
3. Global Existence and Blowup in Finite Time
In this section, we will use the super and subsolution technique to get the global existence or finite time blowup of the solution to (1.1)–(1.3).
Proof of Theorem 1.2.
As , there exist such that
Then we let () be a continuous function satisfying and set
We consider the following auxiliary problem:
Our intention is to show that is a global supersolution of (1.1)–(1.3). Indeed, a direct computation yields
Here we have used the conclusion and inequality (3.1). We still have to consider the boundary and initial conditions. When , in view of Hölder's inequality, we have
Similarly, we have also for that
It is clear that and . Therefore, we get is a global supersolution of (1.1)–(1.3) and hence the solution to (1.1)–(1.3) exists globally by Proposition 2.4.
Proof of Theorem 1.3.
(i) Let be the solution to the homogeneous Dirichlet boundary problem (1.1), (1.4), and (1.3). Then it is well known that for sufficiently large initial data the solution blows up in finite time when (see ). On the other hand, it is obvious that is a subsolution of problem (1.1)–(1.3). Henceforth, the solution of (1.1)–(1.3) with large initial data blows up in finite time provided that .
(ii) We consider the ODE system:
where . Then implies that blows up in finite time (see ). Under the assumption that and for any , is a subsolution of problem (1.1)–(1.3). Therefore, by Proposition 2.4, we see that the solution of problem (1.1)–(1.3) satisfies and then blows up in finite time.
(iii) Let be the positive solution of the linear elliptic problem:
and let be the positive solution of the linear elliptic problem:
where is a positive constant such that . We remark that and ensure the existence of such .
where . We now show that is a supsolution of problem (1.1)–(1.3) for small initial data . Indeed, it follows from that, for ,
Here we used . The above inequalities show that is a supsolution of problem (1.1)–(1.3) whenever . Therefore, system (1.1)–(1.3) has global solutions if and , for any .
4. Blowup Rate Estimate
In this section, we derive the precise blowup rate estimate. To this end, we first establish a partial relationship between the solution components and , which will be very useful in the subsequent analysis. For definiteness, we may assume . If , we can proceed in the same way by changing the role of and and then obtain the corresponding conclusion.
If , and for any , there exists a positive constant such that the solution of problem (1.1)–(1.3) with positive initial data satisfies
Let , where is a positive constant to be chosen. For , a series of calculations show that
If we choose such that , we have
where is a function of and and lies between and .
When , on the other hand, we have
Denote , . Since for any , . It follows from Jensen's inequality, and that
which implies that
For the initial condition, we have
provided that .
Summarily, if we take , then it follows from Theorem in  that , that is,
which is desired.
Using this lemma, we could establish our blowup rate estimate. To derive our conclusion, we shall use some ideas of .
Proof of Theorem 1.4.
For simplicity, we introduce . Let and . A direct computation yields
For , we have from the boundary conditions that
It follows from and Jensen's inequality that the difference in the last brace is nonnegative and thus
By similar arguments, we have
On the other hand, the hypothesis (H) implies that
Hence, from (4.9)–(4.13) and the comparison principle (see Remark 2.3), we get
Let Then and are Lipschitz continuous and thus are differential almost everywhere (see e.g., ). Moreover, we have from equations (1.1) that
We claim that
for some positive constant . Indeed, if we let be the points at which attains its maximum, then relation (4.1) means that
At any point of differentiability of , if ,
From (4.15), (4.18), and (4.19), we can confirm our claim (4.17).
Integrating (4.17) on yields
which gives the upper estimate for . Namely, there exists a constant such that
Then by (4.16) and (4.21), we get
Integrating this equality from to , we obtain
for some positive constant . Thus we have established the upper estimates for .
To obtain the lower estimate for , we notice that (4.16) and (4.18) lead to
for a constant . Integrating above equality on , we see there exists a positive constant such that
Finally, we give the lower estimate for . Indeed, using the relationship (4.16), (4.23) and (4.25), we could prove that is bounded from below; that is, there exists a positive constant such that
To see this, our approach is based on the contradiction arguments. Assume that there would exist two sequences with and with as such that
Then we could choose a corresponding sequence such that , where is a positive constant to be determined later. As , we have
From (4.23) and (4.27), we obtain
Choosing such that , one can get
which would contradict to (4.25) as is large enough since as .
The authors are very grateful to the anonymous referees for their careful reading and useful suggestions, which greatly improved the presentation of the paper. This work is supported in part by Natural Science Foundation Project of CQ CSTC (2007BB2450), China Postdoctoral Science Foundation, the Key Scientific Research Foundation of Xihua University, and Youth Foundation of Science and Technology of UESTC.
Deng, K, Levine, HA: The role of critical exponents in blow-up theorems: the sequel. Journal of Mathematical Analysis and Applications. 243(1), 85–126 (2000). Publisher Full Text
Levine, HA: The role of critical exponents in blowup theorems. SIAM Review. 32(2), 262–288 (1990). Publisher Full Text
Escobedo, M, Herrero, MA: A semilinear parabolic system in a bounded domain. Annali di Matematica Pura ed Applicata. 165, 315–336 (1993). Publisher Full Text
Wang, MX: Blow-up rate estimates for semilinear parabolic systems. Journal of Differential Equations. 170(2), 317–324 (2001). Publisher Full Text
Pao, CV: Asymptotic behavior of solutions of reaction-diffusion equations with nonlocal boundary conditions. Journal of Computational and Applied Mathematics. 88(1), 225–238 (1998). Publisher Full Text
Seo, S: Global existence and decreasing property of boundary values of solutions to parabolic equations with nonlocal boundary conditions. Pacific Journal of Mathematics. 193(1), 219–226 (2000). Publisher Full Text
Kong, L-H, Wang, M-X: Global existence and blow-up of solutions to a parabolic system with nonlocal sources and boundaries. Science in China. Series A. 50(9), 1251–1266 (2007). Publisher Full Text
Wang, YL, Mu, CL, Xiang, ZY: Blowup of solutions to a porous medium equation with nonlocal boundary condition. Applied Mathematics and Computation. 192(2), 579–585 (2007). Publisher Full Text
Yin, H-M: On a class of parabolic equations with nonlocal boundary conditions. Journal of Mathematical Analysis and Applications. 294(2), 712–728 (2004). Publisher Full Text
Yin, YF: On nonlinear parabolic equations with nonlocal boundary condition. Journal of Mathematical Analysis and Applications. 185(1), 161–174 (1994). Publisher Full Text
Zheng, S, Kong, L: Roles of weight functions in a nonlinear nonlocal parabolic system. Nonlinear Analysis: Theory, Methods & Applications. 68(8), 2406–2416 (2008). PubMed Abstract | Publisher Full Text
Souplet, P: Uniform blow-up profiles and boundary behavior for diffusion equations with nonlocal nonlinear source. Journal of Differential Equations. 153(2), 374–406 (1999). Publisher Full Text
Friedman, A, McLeod, B: Blow-up of positive solutions of semilinear heat equations. Indiana University Mathematics Journal. 34(2), 425–447 (1985). Publisher Full Text
Matano, H, Merle, F: Classification of type I and type II behaviors for a supercritical nonlinear heat equation. Journal of Functional Analysis. 256(4), 992–1064 (2009). Publisher Full Text
Xiang, ZY, Hu, XG, Mu, CL: Neumann problem for reaction-diffusion systems with nonlocal nonlinear sources. Nonlinear Analysis: Theory, Methods & Applications. 61(7), 1209–1224 (2005). PubMed Abstract | Publisher Full Text