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.
1. Introduction
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
has been studied extensively (see [3–5] and references therein), and the following proposition was proved.
Proposition 1.1.
(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 [4]). (ii) The asymptotic behavior near the blowup time is characterized by
For the more parabolic problems related to the local boundary, we refer to the recent works [6–9] and references therein.
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 [15] 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 [17], by using some ideas of Souplet [23], obtained the blowup conditions and blowup profile of the following system:
subject to nonlocal boundary (1.2), and Zheng and Kong in [22] 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.
Theorem 1.2.
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)).
Theorem 1.3.
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 [24]. 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 [25]).
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.
Theorem 1.4.
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.
Remark 1.5.
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 .
Remark 1.6.
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 [18], 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.
2. Preliminaries
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).
Definition 2.1.
A pair of functions is called a subsolution of (1.1)–(1.3) if
A supersolution is defined with each inequality reversed.
Lemma 2.2.
Suppose that and are nonnegative functions. If satisfy
then on .
Proof.
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 .
Remark 2.3.
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.
Proposition 2.4.
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:
where is the measure of and . It follows from [13, Theorem ] that exists globally, and indeed (see [13, Theorem ]).
Our intention is to show that is a global supersolution of (1.1)–(1.3). Indeed, a direct computation yields
and thus
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 [4]). 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 [26]). 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 .
Let
where . We now show that is a supsolution of problem (1.1)–(1.3) for small initial data . Indeed, it follows from that, for ,
When ,
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.
Lemma 4.1.
If , and for any , there exists a positive constant such that the solution of problem (1.1)–(1.3) with positive initial data satisfies
Proof.
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 [13] 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 [3].
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
That is,
Let Then and are Lipschitz continuous and thus are differential almost everywhere (see e.g., [24]). 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 .
Acknowledgments
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.
References

Deng, K, Levine, HA: The role of critical exponents in blowup 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

Deng, K: Blowup rates for parabolic systems. Zeitschrift für Angewandte Mathematik und Physik. 47(1), 132–143 (1996). PubMed Abstract  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: Blowup rate estimates for semilinear parabolic systems. Journal of Differential Equations. 170(2), 317–324 (2001). Publisher Full Text

Li, FQ: On initial boundary value problems with equivalued surface for nonlinear parabolic equations. Boundary Value Problems. 2009, (2009)

Xiang, ZY, Mu, CL: Blowup behaviors for degenerate parabolic equations coupled via nonlinear boundary flux. Communications on Pure and Applied Analysis. 6(2), 487–503 (2007)

Xiang, ZY: Global existence and nonexistence for degenerate parabolic equations with nonlinear boundary flux. preprint

Xiang, ZY: Global existence and nonexistence for diffusive polytropic filtration equations with nonlinear boundary conditions. Zeitschrift für Angewandte Mathematik und Physik. In press

Day, WA: A decreasing property of solutions of parabolic equations with applications to thermoelasticity. Quarterly of Applied Mathematics. 40(4), 468–475 (1983)

Day, WA: Heat Conduction within Linear Thermoelasticity, Springer Tracts in Natural Philosophy,p. viii+83. Springer, New York, NY, USA (1985)

Friedman, A: Monotonic decay of solutions of parabolic equations with nonlocal boundary conditions. Quarterly of Applied Mathematics. 44(3), 401–407 (1986)

Deng, K: Comparison principle for some nonlocal problems. Quarterly of Applied Mathematics. 50(3), 517–522 (1992)

Pao, CV: Asymptotic behavior of solutions of reactiondiffusion equations with nonlocal boundary conditions. Journal of Computational and Applied Mathematics. 88(1), 225–238 (1998). Publisher Full Text

Seo, S: Blowup of solutions to heat equations with nonlocal boundary conditions. Kobe Journal of Mathematics. 13(2), 123–132 (1996)

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, LH, Wang, MX: Global existence and blowup 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

Wang, YL, Mu, CL, Xiang, ZY: Properties of positive solution for nonlocal reactiondiffusion equation with nonlocal boundary. Boundary Value Problems. 2007, (2007)

Yin, HM: 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 blowup 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: Blowup 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 reactiondiffusion systems with nonlocal nonlinear sources. Nonlinear Analysis: Theory, Methods & Applications. 61(7), 1209–1224 (2005). PubMed Abstract  Publisher Full Text