- Research
- Open access
- Published:
Uniform blow-up rate for a porous medium equation with a weighted localized source
Boundary Value Problems volume 2011, Article number: 57 (2011)
Abstract
In this article, we investigate the Dirichlet problem for a porous medium equation with a more complicated source term. In some cases, we prove that the solutions have global blow-up and the rate of blow-up is uniform in all compact subsets of the domain. Moreover, in each case, the blow-up rate of |u(t)|∞ is precisely determined.
1 Introduction
Let Ω be a bounded domain in ℝN (N ≥ 1) with smooth boundary ∂Ω. We consider the following parabolic equation with a localized reaction term
where m ≥ 1, q1 ≥ 0, s1 > 0 and x0 ∈ Ω is a fixed point. Throughout this article, we assume the functions a(x) and v0(x) satisfy the following conditions:
(A1) a(x) and v0(x) ∈ C2(Ω); a(x), v0(x) > 0 in Ω and a(x) = v0(x) = 0 on ∂Ω.
When Ω = B = {x ∈ ℝN; |x| < R}, we sometimes assume
(A2) a(x) and v0(x) are radially symmetric; a(r) and v0(r) are non-increasing for r ∈ [0, R].
Problems (1.1)-(1.3) arise in the study of the flow of a fluid through a porous medium with an internal localized source and in the study of population dynamics (see [1–3]). Porous medium equations (m > 1) with or without local sources have been studied by many authors [4–6].
Concerning (1.1)-(1.3), to the best of authors knowledge, a number of articles have studied it from the point of the view of blow-up and global existence [7–10]. Many studies have been devoted to the case m = 1 [10–13]. The case m = 1, a(x) = 1, q1 = 0, s1 ≥ 1 and m = 1, a(x) = 1, q1, s1 > 1 were studied by Souple [10, 11]. Souple [10] demonstrated that the positive solution blows up in finite time if the initial value v0 is large enough. In the case a(x) = 1, q1 = 0, and s1 > 1, Souple [11] showed that the solution v(x, τ) blows up globally and the blow-up rate is precisely determined. The case q1 = 0 and s1 > 0 was studied by Cannon and Yin [12] and Chandam et al. [13]. Cannon and Yin [12] studied its local solvability and Chandam et al. [13] investigated its blow-up properties.
The study of this article is motivated by some recent results of related problems (see [14][15][16]. In the case of a(x)(= constant), the global existence and blow-up behavior have been considered by Chen and Xie [15]. It turns out that if q1 + s1 < m or q1 + s1 = m and a(x)(= constant) is sufficiently small, there exists a global solution of problem (1.1)-(1.3); if q1 + s1 > m, the solution of problem (1.1)-(1.3) blows up for large initial datum while it admits a global solution for small initial datum. Furthermore, Du and Xiang [16] obtained the blow-up rate estimates under some appropriate hypotheses on initial datum. For some related localized models arising in physical phenomena, we refer the readers to [17–19] and the references therein.
For the localized semi-linear parabolic equation of the form
with the Dirichlet boundary condition (1.2) and the initial condition (1.3). In [20], Li and Wang proved that the blow-up set to system (1.2)-(1.4): (a) the system possesses total blow-up when q1 ≤ 1; (b) the system presents single point blow-up patterns when q1 > 1.
We now restrict ourselves to the problem of the form
where q1 ≥ 0, s1 > 0, and q1 + s1 > m > 0. When m = 1, it was proved in [14] that
-
(1)
If 0 ≤ q1 ≤ 1 and q1 + s1 > 1, then the solution of (1.5)-(1.7) blows up in a finite time T.
-
(2)
If q1 > 1, then x = 0 is the only blow-up point for (1.5)-(1.7).
In the meantime, they obtained the blow-up rate estimate but less precise. Namely,
-
(i)
If 0 ≤ q1 < 1, then for any x ∈ B
where , .
-
(ii)
If q1 = 1, then for any x ∈ B
It seems that the results of [14] can be extended to m ≥ 1 and the blow-up rate can be precisely determined. Motivated by this, in this article, we will extend and improve the results of [14].
The purpose of this article is to determine the blow-up rate of solutions for a nonlinear parabolic equation with a weighted localized source, that is, we investigate how the localized source and the local term affect the blow-up properties of the problem (1.5)-(1.7). Indeed, we find that when q1 ≤ 1, the solution of (1.5)-(1.7) blows up at the whole domain with a uniformly blow-up profile.
The rest of this article is organized as follows. The results are stated in Section 2. We then prove these results in Section 3.
2 Preliminaries and Main Results
The following two theorems are our main results.
Theorem 2.1 Assume q1 + s1 > m, (A1) and (A2) hold. Let v(x, t) be the solution of problem (1.5)-(1.7), then v(x,t) blows up provided that the initial value v0(x) is sufficiently large.
The method used in the proof Theorem 2.1 is originally due to [8, 18], and bears much resemblance to that of Theorem 3.2 in [15] and Theorem 1.3 in [16]. Therefore, we omitted them here.
For the case q1 > 1, we do not know how to deal with the uniform blow-up rate of problem (1.5)-(1.7). In the following, we focus only on the case of 0 ≤ q1 ≤ 1.
Theorem 2.2 Assume (A1) and (A2) hold. Let v(x, t) be the blow-up solution of (1.5)-(1.7), which blows up in finite time T and v(x, t) is non-decreasing in time, then the following limits hold uniformly in all compact subsets of B.
-
(i)
If 0 ≤ q1 < 1, then
(2.1)
where .
-
(ii)
If q1 = 1, then
(2.2)
Remark 2.1 The domain we considered here is a ball, it seems that the results of Theorem 2.2 remain valid for the general domain. (It is an open problem in this case.)
To get the blow-up profiles for problem (1.5)-(1.7), we need some transformations. Let u(x, t) = vm(x, τ), t = mτ, then (1.5)-(1.7) becomes
where 0 ≤ p = (m - 1)/m < 1, q = q1/m, and s = s1/m.
Under above transformation, assumptions (A1) and (A2) become
(B1) a(x) and u0(x) ∈ C2(B); a(x), u0(x) > 0 in B and a(x) = u0(x) = 0 on ∂B.
(B2) a(x) and u0(x) are radially symmetric; a(r) and u0(r) are non-increasing for r ∈ [0, R].
In our consideration, a crucial role is played by the Dirichlet eigenvalue problem
Denote λ be the first eigenvalue and by φ the corresponding eigenfunction with φ(x) > 0 in B, normalized by .
3 Proof of Theorem 2.2
For convenience, we denote
Before proving our result, we would like to give a property of the following problem
where 0 ≤ α ≤ 1 and w = u1-q(x, t).
Lemma 3.1 Assume (B1) and (B2) hold. Let w(x, t) be the solution of Equation (3.1), which blows up in a finite time T* and non-decreasing in time t, then the following limits hold uniformly in all compact subsets of B.
-
(i)
If 0 ≤ α < 1, then
-
(ii)
If α = 1, then
Proof. (i) Assumption (B2) implies w r ≤ 0 (r = |x|), it then follows that and Δw(0, t) ≤ 0 for t > 0. From (3.1), we then get
.
Consequently,
which implies
Moreover, it is apparent that limt-T*w(0, t)/g(t) = 0, since s > 1 - q.
Set R1 ∈ (0, R), B1 = {x ∈ ℝN, | x |< R1} and b(x) = 1/a(x), x ∈ B1. Since a'(r) ≤ 0, we obtain that b'(r) ≥ 0, for 0 ≤ r ≤ R1.
We now introduce the function
By a simple calculation, and note that ∇w(x, t)∇b(x) = u r (r, t)b'(r) ≤ 0, then there exist m1, m2 > 0 such that
Setting ε(t) = m2w(0, t)/g(t). From limt→T*w(0, t)/g(t) = 0, we infer that there exists t1 ∈ (0, T*) such that 0 < ε(t) ≤ 1/2 for t1 ≤ t < T*.
Hence, in view of (3.1), we observe
Set g1(t) = (1 - ε(t))g(t), , we then obtain
Obviously, w1(x, t) is a sup-solution of the following equation
By the maximum principle, w1(x, t) ≥ w*(x, t) and . Similar to the proof of (4.15) in [15] that
uniformly in all compact subsets of B1,
Therefore, by the arbitrariness of B1, we obtain that the following inequatlity holds uniformly in all compact subsets of B
In particular,
From (3.2) and (3.4), we deduce
Multiplying both sides of (3.1) by φ and integrating over B × (0, t), 0 < t < T*
Since , so we have
It then follows that
Note that w r ≤ 0, (3.3) and (3.6), it is sufficient to prove
Assume on the contrary that there exists a point x1 ∈ B, x1 ≠ 0 such that
Then there exists a sequence {t n } such that t n → T*
By the continuity of a(x), we deduce that there exists x2 ∈ B such that (1 - α)a(x) < c for B1 = {x ∈ ℝn : |x2| ≤ |x| ≤ |x1|}. Using w r ≤ 0, (3.3) and (3.6), it is easy to check that
which is a contradiction to (3.6). Combining (3.3) and (3.7), Lemma 3.1 (i) is proved. Case (ii) can be treated similarly.
The key step in establishing the result of Theorem 2.2 is the following lemma.
Lemma 3.2 Under the assumption of Lemma 3.1, let u(x, t) be the blow-up solution of (2.3), which blows up in a finite time T* and non-decreasing in time t, then the following statements hold uniformly in all compact subsets of B:
-
(i)
If p + q < 1, then
-
(ii)
If p + q = 1, then
Proof. (i) Since u r ≤ 0 and u t ≥ 0, it then follows that and Δu(0, t) ≤ 0 for t > 0, which imply limt→T*u(0, t) = ∞. Obviously,
which implies
Notice that p + q < 1 and (3.8), hence limt→T*G(t) = ∞ and limt→T*g(t) = ∞.
A simple calculation yields
In view of (2.3), we have, for x ∈ Ω, 0 < t < T*
Multiplying both sides of Equation (3.9) by φ and integrating over B × (0, t), it follows that
for 0 < t < T*. Clearly,
which yields
Setting u1(r, t) = u(1-q)/ 2(r, t)(r = |x|). We may claim that
Indeed, due to limt→T*g(t) = limt→T*us(0, t) = ∞, u r ≤ 0, and s > 1 - q, we then have
Therefore, by Lebesgue's dominated convergence theorem, we infer that
where w n is the surface area of unit ball in ℝN.
Now according to (3.10)-(3.12), we obtain
On the other hand, By (3.9), we find
where γ = p/(1 - q). Consequently, u1-qis a sup-solution of the problem
By the maximum principle, u1-q≥ v in B × (0, T*). Note that 0 ≤ γ < 1, we know from Lemma 3.1 (i) that
uniformly in all compact subsets of B.
Thus,
uniformly in all compact subsets of B.
Next, we prove that
uniformly in all compact subsets of B.
We can verify (3.15) by similar means of (3.7). Therefore, we conclude the proof of case (i).
-
(ii)
Proceeding as (3.8), we have
For any compact subset B1 ∈ B, there exists t1 ∈ (0, T*) such that u(x, t1) ≥ 1 for all , and thus ln u(x, t) ≥ 0 in .
Direct calculation shows
Let λ1 be the first eigenvalue of -Δ in and by φ1 > 0 the corresponding eigenfunction, normalized by . Set . Clearly, limt→T*G(t)/G1(t) = 1.
Multiplying both sides of Equation (3.16) by φ1 and integrating over B1 × (t1, t), we get
The result of case (ii) follows by analogy with the argument used in the proof of case (i).
Proof of Theorem 2.2
-
(i)
By Lemma 3.2 (i), we infer that
hence
Integrating equivalence (3.18) between t and T*, we obtain
Using Lemma 3.2 (i) and substituting p = (m - 1)/m, q = q1/m, s = s1/m, t* = mτ, and u(x, t) = vm(x, τ) into (3.19), we complete the proof of Theorem 2.2 (i).
-
(ii)
To obtain the blow-up rate of the exponent type, we need to be more careful in this case, since exponentiation of equivalents is not permitted. Similar to the proof of Theorem 3 in [14] and Lemma 2.3 in [16], we get
(3.21)
Thanks to Lemma 3.2(ii) and (3.20), we then get the desired result.
4 Discussion
This article deals with the porous medium equation with local and localized source terms, represented by two factors and , respectively. As we all know that, in the absence of weight function, the solutions of model (1.5)-(1.7) have a global blow-up and the rate of blow-up is uniform in all compact subsets of the domain. A natural question is what happens in the model (1.5)-(1.7), where the source term is the product of localized source, local source, and weight function. It is shown by Theorem 2.2 that if 0 ≤ q1 ≤ 1, this equation possesses uniform blow-up profiles. In other words, the localized term plays a leading role in the blow-up profile for this case. Moreover, the blow-up rate estimates in time and space is obtained.
References
Diaz J, Kerser R: On a nonlinear degenerate parabolic equation in infiltration or evaporation through a porous medium. J Diff Equ 1987, 69: 368–403. 10.1016/0022-0396(87)90125-2
Furter J, Grinfeld M: Local vs. nonlocal interactions in population dynamics. J Math Biol 1989, 27: 65–80. 10.1007/BF00276081
Okada A, Fukuda I: Total versus single point blow-up of solution of a semilinear parabolic equation with localized reaction. J Math Anal Appl 2003, 281: 485–500. 10.1016/S0022-247X(03)00133-1
Cantrell R, Cosner C: Diffusive logistic equation with indefinite weights: population models in disrupted environments. II SIAM J Math Anal 1991, 22: 1043–1064. 10.1137/0522068
Levine H: The role of critical exponents in blow-up theorem. SIMA Rev 1990, 32: 268–288.
Anderson J: Local existence and uniqueness of solutions of degenerate parabolic equations. Commun Partial Diff Equ 1991, 16: 105–143. 10.1080/03605309108820753
Chen Y, Liu Q, Gao H: Bounedeness of global solutions of a porous medium equation with a localized source. Nonlinear Anal 2006, 64: 2168–2182. 10.1016/j.na.2005.08.004
Fukuda I, Suzuki R: Blow-up bebavior for a nolinear beat equation with a localized source in a ball. J Diff Equ 2005, 218: 273–291. 10.1016/j.jde.2005.02.013
Chen Y, Liu Q, Gao H: Bounedeness of global positive solutions of a porous medium equation with a moving localized source. J Math Anal Appl 2007, 333: 1008–1023. 10.1016/j.jmaa.2006.11.048
Souple P: Blow-up in non-local reaction-diffusion equations. SIAM J Math Anal 1998, 29(6):1301–1334. 10.1137/S0036141097318900
Souple P: Uniform blow-up profiles and boundary for diffusion equations with nonlocal nonlinear source. J Diff Equ 1999, 153: 374–406. 10.1006/jdeq.1998.3535
Cannon R, Yin M: A class of non-linear non-classical parabolic equations. J Diff Equ 1989, 79: 226–288.
Chandam J, Peirce A, Yin H: The blow-up property of solutions to some diffusion equations with localized nonlinear reactions. J Math Anal Appl 1992, 169: 313–328. 10.1016/0022-247X(92)90081-N
Kong L, Wang L, Zheng S: Asymptotic analysis to a parabolic equation with a weighted localized source. Appl Math Comput 2008, 197: 819–827. 10.1016/j.amc.2007.08.016
Chen Y, Xie C: Blow-up for a porous medium equation with a localized source. Appl Math Comput 2004, 159: 79–93. 10.1016/j.amc.2003.10.032
Du L, Xiang Z: A further blow-up analysis for a localized porous medium equation. Appl Math Comput 2006, 179: 200–208. 10.1016/j.amc.2005.11.100
Wang J, Kong L, Zheng S: Asymptotic analysis for a localized nonlinear diffusion equation. Comput Math Appl 2008, 56: 2294–2304. 10.1016/j.camwa.2008.03.057
Friedman A, Mcleod J: Blow-up of positive solutions of semilinear heat equations. Indiana Univ Math J 1985, 34: 425–447. 10.1512/iumj.1985.34.34025
Rouchon P: Boundeness of global solutions of nonlinear diffusion equations with localized reaction term. Diff Integral Equ 2003, 16(9):1083–1092.
Li H, Wang M: Properties of blow-up solutions to a parabolic system with nonlinear localized terms. Discrete Contin Dyn Syst 2005, 13: 683–700.
Acknowledgements
The authors thank the anonymous referee for their constructive and valuable comments, which helped in improving the presentation of this study. This study was supported by the National Natural Science Foundation of China (60972001), the National Key Technologies R & D Program of China (2009BAG13A06), the Scientific Innovation Research of College Graduate in Jiangsu Province (CXZZ_0163), and the Scientific Research Foundation of Graduate School of Southeast University (YBJJ1140).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
All the authors typed, read, and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Zeng, W., Lu, X. & Liu, Q. Uniform blow-up rate for a porous medium equation with a weighted localized source. Bound Value Probl 2011, 57 (2011). https://doi.org/10.1186/1687-2770-2011-57
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/1687-2770-2011-57