Skip to main content

Non-Newtonian polytropic filtration systems with nonlinear boundary conditions

Abstract

This article deals with the global existence and the blow-up of non-Newtonian polytropic filtration systems with nonlinear boundary conditions. Necessary and sufficient conditions on the global existence of all positive (weak) solutions are obtained by constructing various upper and lower solutions.

Mathematics Subject Classification (2000)

35K50, 35K55, 35K65

Introduction

In this article, we study the global existence and the blow-up of non-Newtonian polytropic filtration systems with nonlinear boundary conditions

(1.1)

where

Ω N is a bounded domain with smooth boundary ∂Ω, ν is the outward normal vector on the boundary ∂Ω, and the constants k i , m i > 0, m ij ≥ 0, i, j = 1,..., n; ui 0(x) (i = 1,..., n) are positive C1 functions, satisfying the compatibility conditions.

The particular feature of the equations in (1.1) is their power- and gradient-dependent diffusibility. Such equations arise in some physical models, such as population dynamics, chemical reactions, heat transfer, and so on. In particular, equations in (1.1) may be used to describe the nonstationary flows in a porous medium of fluids with a power dependence of the tangential stress on the velocity of displacement under polytropic conditions. In this case, the equations in (1.1) are called the non-Newtonian polytropic filtration equations which have been intensively studied (see [14] and the references therein). For the Neuman problem (1.1), the local existence of solutions in time have been established; see the monograph [4].

We note that most previous works deal with special cases of (1.1) (see [513]). For example, Sun and Wang [7] studied system (1.1) with n = 1 (the single-equation case) and showed that all positive (weak) solutions of (1.1) exist globally if and only if m11k1 when k1m1; and exist globally if and only if when k1 > m1. In [13], Wang studied the case n = 2 of (1.1) in one dimension. Recently, Li et al. [5] extended the results of [13] into more general N-dimensional domain.

On the other hand, for systems involving more than two equations when m i = 1(i = 1,..., n), the special case k i = 1(i = 1,..., n) (heat equations) is concerned by Wang and Wang [9], and the case k i ≤ 1(i = 1,..., n) (porous medium equations) is discussed in [12]. In both studies, they obtained the necessary and sufficient conditions to the global existence of solutions. The fast-slow diffusion equations (there exists i(i = 1,..., n) such that k i > 1) is studied by Qi et al. [6], and they obtained the necessary and sufficient blow up conditions for the special case Ω = B R (0) (the ball centered at the origin in N with radius R). However, for the general domain Ω, they only gave some sufficient conditions to the global existence and the blow-up of solutions.

The aim of this article is to study the long-time behavior of solutions to systems (1.1) and provide a simple criterion of the classification of global existence and nonexistence of solutions for general powers k i m i , indices m ij , and number n.

Define

Our main result is

Theorem. All positive solutions of (1.1) exist globally if and only if all of the principal minor determinants of A are non-negative.

Remark. The conclusion of Theorem covers the results of [513]. Moreover, this article provides the necessary and sufficient conditions to the global existence and the blow-up of solutions in the general domain Ω. Therefore, this article improves the results of [6].

The rest of this article is organized as follows. Some preliminaries will be given in next section. The above theorem will be proved in Section 3.

Preliminaries

As is well known that degenerate and singular equations need not possess classical solutions, we give a precise definition of a weak solution to (1.1).

Definition. Let T > 0 and Q T = Ω × (0, T]. A vector function (u1(x, t),.., u n (x, t)) is called a weak upper (or lower) solution to (1.1) in Q T if

(i). ;

(ii). (u1(x, 0),..., u n (x, 0)) ≥ (≤)(u10(x),..., un 0(x));

(iii). for any positive functions ψ i (i = 1,..., n) L1(0, T; W1,2(Ω)) ∩ L2(Q T ), we have

In particular, (u1(x, t),..., u n (x, t)) is called a weak solution of (1.1) if it is both a weak upper and a lower solution. For every T < ∞, if (u1(x, t),..., u n (x, t)) is a solution of (1.1) in Q T , then we say that (u1(x, t),..., u n (x, t)) is global.

Lemma 2.1 (Comparison Principle.) Assume that ui 0(i = 1,..., n) are positivefunctions and (u1,..., u n ) is any weak solution of (1.1). Also assume that (u1,..., u n ) ≥ (δ,..., δ) > 0 andare the lower and upper solutions of (1.1) in Q T , respectively, with nonlinear boundary fluxand, where. Then we havein Q T .

When n = 2, the proof of Lemma 2.1 is given in [5]. When n > 2, the proof is similar.

For convenience, we denote , which are fixed constants, and let .

In the following, we describe three lemmas, which can be obtained directly from Lemmas 2.7-2.9 in [6].

Lemma 2.2 Suppose all the principal minor determinants of A are non-negative. If A is irreducible, then for any positive constant c, there exists α = (α1,..., α n ) T such that A α ≥ 0 and α i > c (i = 1,..., n).

Lemma 2.3 Suppose that all the lower-order principal minor determinants of A are non-negative and A is irreducible. For any positive constant C, there exist large positive constants L i (i = 1,..., n) such that

Lemma 2.4 Suppose that all the lower-order principal minor determinants of A are non-negative and |A| < 0. Then, A is irreducible and, for any positive constant C, there exists α = (α1,..., α n ) T , with α i > 0 (i = 1,..., n) such that

Proof of Theorem

First, we note that if A is reducible, then the full system (1.1) can be reduced to several sub-systems, independent of each other. Therefore, in the following, we assume that A is irreducible. In addition, we suppose that k1 - m1k2 - m2 ≤ · · · k n - m n .

Let be the first eigenfunction of

(3.1)

with the first eigenvalue , normalized by , then , in Ω and and on ∂Ω (see [1416]).

Thus, there exist some positive constants , , , and such that

(3.2)

We also have provided with and some positive constant . For the fixed , there exists a positive constant such that if .

Proof of the sufficiency. We divide this proof into three different cases.

Case 1. (k i < m i (i = 1,..., n)). Let

(3.3)

where Q i satisfies , and constants P i , α i (i = 1,..., n) remain to be determined. Since , by performing direct calculations, we have

in Ω × +. By setting if m i ≥ 1, if m i < 1, we have one the boundary that

we have

if

(3.4)

and

(3.5)

Note that k i < m i (i = 1,..., n). From Lemmas 2.2 and 2.3, we know that inequalities (3.4) and (3.5) hold for suitable choices of P i , α i (i = 1,..., n). Moreover, if we choose P i , α i to be large enough such that

then , . Therefore, we have proved that is a global upper solution of the system (1.1). The global existence of solutions to the problem (1.1) follows from the comparison principle.

Case 2. (k i m i (i = 1,..., n)). Let

(3.6)

where if m i ≥ 1, if m i < 1, , , , are defined in (3.1) and (3.2), α i (i = 1,..., n) are positive constants that remain to be determined, and

Since -ye-y≥ -e-1 for any y > 0, we know that . Thus, for (x, t) Ω × +, a simple computation shows that

In addition, we have

Noting on ∂Ω, we have on the boundary that

Then, we have

if

(3.7)

From Lemma 2.2, we know that inequalities (3.7) hold for suitable choices of α i (i = 1,..., n). Moreover, if we choose ∞ i to be large enough such that

then . Therefore, we have shown that is an upper solution of (1.1) and exists globally. Therefore, , and hence the solution (u1,..., u n ) of (1.1) exists globally.

Case 3. (k i < m i (i = 1,..., s); k i m i (i = s + 1,..., n)). Let be as in (3.3) and

where , and A i are as in case 2. By Lemma 2.3, we choose P i ≥ (log Q i )-1||ui 0|| (i = 1,..., s) and M i ≥ max{1, ||ui 0||} (i = s + 1,..., n) such that

(3.8)

Set

By similar arguments, in cases 1 and 2, we have on the boundary that

Therefore employing (3.8), we see that

if we knew

(3.9)

We deduce from Lemma 2.2 that (3.9) holds for suitable choices of α i (i = 1,..., n). Moreover, we can choose α i large enough to assure that

Then, as in the calculations of cases 1 and 2, we have . We prove that is an upper solution of (1.1), so (u1,..., u n ) exists globally.

Proof of the necessity.

Without loss of generality, we first assume that all the lower-order principal minor determinants of A are non-negative, and |A| < 0, for, if not, there exists some l th-order (1 ≤ l < n) principal minor determinant detAl × lof A = (a ij )n×nwhich is negative. Without loss of generality, we may consider that

and all of the sth-order (1 ≤ sl - 1) principal minor determinants detAs × sof Al × lare non-negative. Then, we consider the following problem:

(3.10)

Note that . If we can prove that the solution (w1,..., w l ) of (3.10) blows up in finite time, then (w1,... w l , δ,..., δ) is a lower solution of (1.1) that blows up in finite time. Therefore, the solution of (1.1) blows up in finite time.

We will complete the proof of the necessity of our theorem in three different cases.

Case 1. (k i < m i (i = 1,..., n)). Let

(3.11)

where , , , , the α i are as given in Lemma 2.4 and satisfy ,

(3.12)

By direct computation for , we have

For , we have

Thus, by (3.12) and Lemma 2.4, we have

We confirm that (u1,..., u n ) is a lower solution of (1.1), which blows up in finite time. We know by the comparison principle that the solution (u1,..., u n ) blows up in finite time.

Case 2. (k i m i (i = 1,..., n)). Let if m i < 1, if m i ≥ 1. for k i m i (i = 1,..., n), set

(3.13)

where α i (i = 1,..., n) are to determined later and

(3.14)
(3.15)
(3.16)

By a direct computation, for x Ω, 0 < t < c/b, we obtain that

(3.17)

If , we have , and thus

(3.18)

On the other hand, since -ye-y≥ -e-1 for any y > 0, we have

(3.19)

We have by (3.16), (3.18), and (3.19) that .

If , then , and then

(3.20)

It follows from (3.16), (3.17), and (3.20) that .

We have on the boundary that

(3.21)

Moreover, by (3.14) and Lemma 2.4, we have that

(3.22)

(3.15), (3.21), and (3.22) imply that . Therefore, (u1,..., u1) is a lower solution of (1.1).

For k i = m i (i = 1,..., n), let

(3.23)

For k i = m i (i = 1,..., s) and k i > m i (i = s + 1,..., n), let as in (3.13) and (3.23). Using similar arguments as above, we can prove that (u1,..., u n ) is a lower solution of (1.1). Therefore, (u1,..., u n ) ≤ (u1,..., u n ). Consequently, (u1,..., u n ) blows up in finite time.

Case 3. (k i < m i (i = 1,..., s); k i m i (i = s + 1,..., n)). Let be as in (3.11) and

where α i 's are to determined later and

Based on arguments in cases 1 and 2, we have for . Furthermore, for , we have

Thus,

holds if

(3.24)

From Lemma 2.4, we know that inequalities (3.24) hold for suitable choices of α i (i = 1,..., n). We show that (u1,..., u n ) is a lower solution of (1.1). Since (u1,..., u n ) blows up in finite time, it follows that the solution of (1.1) blows up in finite time.

References

  1. Kalashnikov AS: Some problems of the qualitative theory of nonlinear degenerate parabolic equations of second order. Russ Math Surv 1987, 42: 169-222.

    Article  MathSciNet  MATH  Google Scholar 

  2. Li ZP, Mu CL: Critical exponents for a fast diffusive polytropic filtration equation with nonlinear boundary condition. J Math Anal Appl 2008, 346: 55-64. 10.1016/j.jmaa.2008.05.035

    Article  MathSciNet  MATH  Google Scholar 

  3. Vazquez JL: The Porous Medium Equations: Mathematical Theory. Oxford Mathematical Monographs. Oxford University Press, Oxford; 2007.

    Google Scholar 

  4. Wu ZQ, Zhao JN, Yin JX, Li HL: Nonlinear Diffusion Equations. Word Scientific Publishing Co. Inc., River Edge, NJ; 2001.

    Google Scholar 

  5. Li ZP, Mu CL, Li YH: Blowup for a non-Newtonian polytropic filtration system coupled via nonlinear boundary flux. Boundary Value Problems 2008., 2008: Article ID 847145

    Google Scholar 

  6. Qi YW, Wang MX, Wang ZJ: Existence and non-existence of global solutions of diffusion systems with nonlinear boundary conditions. Proc R Soc Edinb A 2004, 134: 1199-1217. 10.1017/S030821050000370X

    Article  MathSciNet  MATH  Google Scholar 

  7. Sun WJ, Wang S: Nonlinear degenerate parabolic equation with nonlinear boundary condition. Acta Mathematica Sinica, English Series 2005, 21: 847-854. 10.1007/s10114-004-0512-2

    Article  MathSciNet  MATH  Google Scholar 

  8. Wang MX: Long time behaviors of solutions of quasilinear equations with nonlinear boundary conditions. Acta Math Sinica (in Chinese) 1996, 39: 118-124.

    MATH  Google Scholar 

  9. Wang MX, Wang YM: Reaction diffusion systems with nonlinear boundary conditions. Sci China A 1996, 27: 834-840.

    Google Scholar 

  10. Wang MX, Wang S: Quasilinear reaction-diffusion systems with nonlinear boundary conditions. J Math Anal Appl 1999, 231: 21-33. 10.1006/jmaa.1998.6220

    Article  MathSciNet  MATH  Google Scholar 

  11. Wang MX: Fast-slow diffusion systems with nonlinear boundary conditions. Nonlinear Anal 2001, 46: 893-908. 10.1016/S0362-546X(00)00156-5

    Article  MathSciNet  MATH  Google Scholar 

  12. Wang S, Wang MX, Xie CH: Quasilinear parabolic systems with nonlinear boundary conditions. J Differential Equations 2000, 166: 251-265. 10.1006/jdeq.2000.3784

    Article  MathSciNet  MATH  Google Scholar 

  13. Wang S: Doubly nonlinear degenerate parabolic systems with coupled nonlinear boundary conditions. J Differential Equations 2002, 182: 431-469. 10.1006/jdeq.2001.4101

    Article  MathSciNet  MATH  Google Scholar 

  14. Lindqvist P: On the equation div ( up-2 u ) + λ up-2u = 0. Proc Am Math Soc 1990, 109: 157-164.

    MathSciNet  MATH  Google Scholar 

  15. Lindqvist P: On the equation div ( up-2 u ) + λ up- 2u = 0. Proc Am Math Soc 1992, 116: 583-584.

    MathSciNet  MATH  Google Scholar 

  16. Li YX, Xie CH: Blow-up for p-Laplacian parabolic equations. J Differential Equations 2003, 20: 1-12.

    Google Scholar 

Download references

Acknowledgements

This study was partially supported by the Projects Supported by Scientific Research Fund of SiChuan Provincial Education Department(09ZC011), and partially supported by the Natural Science Foundation Project of China West Normal University (07B046).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Wanjuan Du.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

DW carried out all studies in the paper. LZ participated in the design of the study in the paper.

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.

Reprints and permissions

About this article

Cite this article

Du, W., Li, Z. Non-Newtonian polytropic filtration systems with nonlinear boundary conditions. Bound Value Probl 2011, 2 (2011). https://doi.org/10.1186/1687-2770-2011-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2011-2

Keywords