SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

This article is part of the series Degenerate and Singular Differential Operators with Applications to Boundary Value Problems.

Open Access Research Article

Quenching for a Reaction-Diffusion System with Coupled Inner Singular Absorption Terms

Shouming Zhou* and Chunlai Mu

Author Affiliations

College of Mathematics and Statistics, Chongqing University, Chongqing 400044, China

For all author emails, please log on.

Boundary Value Problems 2010, 2010:797182  doi:10.1155/2010/797182


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2010/1/797182


Received:13 May 2010
Accepted:5 July 2010
Published:20 July 2010

© 2010 The Author(s) Shouming Zhou and Chunlai Mu.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

we devote to investigate the quenching phenomenon for a reaction-diffusion system with coupled singular absorption terms, , . The solutions of the system quenches in finite time for any initial data are obtained, and the blow-up of time derivatives at the quenching point is verified. Moreover, under appropriate hypotheses, the criteria to identify the simultaneous and nonsimultaneous quenching are found, and the four kinds of quenching rates for different nonlinear exponent regions are given. Finally, some numerical experiments are performed, which illustrate our results.

1. Introduction

This paper deals with the following nonlinear parabolic equations with null Neumann boundary conditions:

(11)

where for , is a bounded domain with smooth boundary, the initial data and are positive, smooth, and compatible with the boundary data.

Because of the singular nonlinearity inner absorption terms of (1.1), the so-called finite-time quenching may occur for the model. We say that the solution of the problem (1.1) quenches, if there exists a time ( denotes the quenching time, denotes quenching point), such that

(12)

For a quenching solution of (1.1), the inf norm of one of the components must tend to as tends to the quenching time . The case when quenches and remains bounded from zero is called non-simultaneous quenching. We will call the case, when both components and quench at the same time, as simultaneous quenching. The purpose of this paper is to find a criteria to identify simultaneous and non-simultaneous quenching for (1.1) and then establish quenching rates for the different cases.

In order to motivate the main results for system (1.1), we recall some classical results for the related system. de Pablo et al., firstly distinguished non-simultaneous quenching from simultaneous one in [1]. They considered a heat system coupled via inner absorptions as follows:

(13)

Recently, Zheng and Wang deduced problem (1.3) to -dimensional with positive Dirichlet boundary condition in [2]. Then, Zhou et al. have given a natural continuation for problem (1.3) beyond quenching time for the case of non-simultaneous quenching in [3].

Replacing the coupled inner absorptions in (1.1) by the coupled boundary fluxes, one gets

(14)

Recently, the simultaneous and non-simultaneous quenching for problem (1.4), and what is related to it, was studied by many authors (see [47] and references therein).

In order to investigate the problem (1.1), it is necessary to recall the blow-up problem of the following reaction-diffusion system:

(15)

with positive powers has been extensively studied by many authors for various problems such as global existence and finite time blow-up, Fujita exponents, non-simultaneous and simultaneous blow-up, and blow-up rates, (see [810] and references therein). However, unlike the blow-up problem, there are less papers consider the weakly coupled quenching problem like (1.1), differently from the generally considered, there are two additional singular factors, namely, and for the inner absorptions of and , respectively. In this paper, we will show real contributions of the two additional singular factors to the quenching behavior of solutions. Our main results are stated as follows.

Theorem 1.1.

If and , then the solution of the system (1.1) quenches in finite time for every initial data.

On the other hand, some authors understand quenching as blow-up of time derivatives while the solution itself remains bounded (see [1113]). In present paper, we assume that the initial data satisfy

(16)

Theorem 1.2.

Let and the radial initial function satisfies (1.6), then blows up in finite time.

Next, we characterize the ranges of parameters to distinguish simultaneous and non-simultaneous quenching. In order to simplify our work, we deal with the radial solutions of (1.1) with , and the radial increasing initial data satisfies (1.6). Thus we, see that is the only quenching point (see [2, 14]). Without loss of generality, we only consider the non-simultaneous quenching with remaining strictly positive, and our main results are stated as follows.

Theorem 1.3.

If and , then any quenching in (1.1) must be simultaneous.

Theorem 1.4.

If and , then any quenching in (1.1) is non-simultaneous with being strictly positive.

Theorem 1.5.

If and , then both simultaneous and non-simultaneous quenching may occur in (1.1) depending on the initial data.

Remark 1.6.

In particular, if we choose , then we obtain that the ranges of parameters to distinguish simultaneous and non-simultaneous quenching coincide with the problem (1.3) (see [1, 2]). Moreover, this criteria to identify the simultaneous and non-simultaneous quenching is the same with the problem (1.4) which coupled boundary fluxes (see [6]). This situation also happens for the blow-up problem (see [8, 10, 15]).

Next, we deal with quenching rates. To state our results more conveniently, we introduce the notation which means that there exist two finite positive constants such that , and the two parameters and verifying

(17)

or equivalently,

(18)

In terms of parameters and , the quenching rates of problem (1.1) can be shown as follow.

Theorem 1.7.

If quenching is non-simultaneous and, for instance, is the quenching variable, then .

Theorem 1.8.

If quenching is simultaneous, then for close to , we have

(i) for , or ;

(ii) for and ;

(iii) for and .

The plan of this paper is organized as follows. In Section 2, we distinguish non-simultaneous quenching from simultaneous one. The four kinds of non-simultaneous and simultaneous quenching rates for different nonlinear exponent regions are given in Section 3. In the Section 4, we perform some numerical experiments which illustrate our results.

2. Simultaneous and Non-Simultaneous Quenching

Proof of Theorem 1.1.

Assume that is the classical solution of (1.1) with the maximal existence time . The maximum principle implies and in . Let . Hence, integrating (1.1) in space and using Green's formula, we have

(21)

Consequently,

(22)

Thus, the solution of the problem (1.1) quenches in finite time. The prove of Theorem 1.1 is complete.

In order to prove Theorem 1.2, we need the following Lemma.

Lemma 2.1.

Assume that and the radial nondecreasing initial data satisfy (1.6), then there exists a small such that

(23)

Proof.

Let . Thus,

(24)

Since and are radial and nondecreasing in , we have . A similar computation holds for , and we obtain

(25)

with boundary conditions

(26)

From (1.6), it is easy to deduce in (see [13, 14]). Choosing small enough, we have that the initial data verifying

(27)

Hence, by the comparison result, we derive that

(28)

This proves Lemma 2.1.

Proof of Theorem 1.2.

This theorem is the direct result of Theorem 1.1 and Lemma 2.1.

Next, we characterize the ranges of parameters to distinguish simultaneous and non-simultaneous quenching. By the hypothesis on the initial data, we obtain and for (see [2, 14]). We collect the estimates of the time derivatives obtained before. Clearly, the only quenching point is (see [2]), we only care for the original point,

(29)

(210)

Proof of Theorem 1.3.

We argue by contradiction. Assume that there exists such that on and quenching at the time . Through (2.10), we have , integrating from to we get . Together with (2.9) we have . Integrating in , we obtain

(211)

If , we have the left hand of the above inequality diverged. So, we get a contradiction. The proof of Theorem 1.3 is finished.

Proof of Theorem 1.4.

First, assume that and . Combining (2.9) with (2.10), we get

(212)

Since , integrating the first inequality in the (2.12) from to , we have

(213)

where are positive constants, the above inequality requires that remains positive up to the quenching time. The case can be treated in an analogous way. The proof of Theorem 1.4 is complete.

Proof of Theorem 1.5.

If and the initial data on , thus, it is easy to see that for problem (1.1) simultaneous quenching occurs.

On the other hand, we want to choose small in order that the quenching time (through Theorem 1.1, we get ) be so small that does not have time to vanish.

Let be fixed. From in , we obtain

(214)

Together with the estimate (2.12), we get

(215)

Integrating in , we obtain

(216)

It is easy to see that the last term of the above inequality is strictly positive, if is small enough and , therefore, we prove that under the condition and , for the solution of (1.1) non-simultaneous quenching may occur. The proof of Theorem 1.5 is complete.

3. Quenching Rates

In this section, we deal with the all possible quenching rates in model (1.1).

Proof of Theorem 1.7.

Under the condition of Theorem 1.7, it holds that . By (2.10), we have

(31)

Thus,

(32)

The proof of Theorem 1.7 is complete.

Proof of Theorem 1.8.

(i) Assume that the quenching of problem (1.1) is simultaneous with , integrating (2.12) yields

(33)

where . Since we assume that quench at , we have as .

On the other hand, from and , we get, a positive constant such that

(34)

Similarly, we can show that there exists a positive constant such that

(35)

Consequently,

(36)

Recalling the estimates (2.9) and (2.10), we obtain

(37)

Integrating from to , we get

(38)

If and , we deduce the quenching rate by a bootstrap argument. First, by (2.9), we get , it follows that . Employing (2.10), we get , that is, . Repeating this procedure, we obtain , , where satisfy

(39)

One can check that ( define by (1.8)), and the all positive constants are bounded. Therefore, passing to the limit, we get . The reverse inequalities can be obtained in the same way.

(ii) If and , we have . It is easy to see that as , from (2.9) and (2.10), we obtain

(310)

 (iii) If and , from (2.9), we get

(311)

Recalling the estimate (2.10), we get

(312)

that is,

(313)

Let , we have

(314)

It is known that the incomplete Gamma function satisfies for . With , we obtain

(315)

and hence,

(316)

Next, we deduce the behaviour for . Combining with (2.9) and (3.16), we have

(317)

Integrating from to ,

(318)

Setting , we get

(319)

For the incomplete Gamma function with , we obtain

(320)

The proof of Theorem 1.8 is complete.

4. Numerical Experiments

In this section, we perform some numerical experiments, which illustrate our results. Now we introduce the numerical scheme for the space discretization, we discretize applying linear finite elements with mass lumping in a uniform mesh for the space variable and keeping continuous, it is well known that this discretization in space coincides with the classic central finite difference second-order scheme, (see [16]), Mass lumping is widely used in parabolic problems with blow-up and quenching, (see, e.g., [17, 18]).

Let us consider the uniform partition of size of the interval , (), and its associated standard piecewise linear finite element space . The semidiscrete approximation obtained by the finite element method with mass lumping is defined as

(41)

where the superindex denotes the Lagrange interpolation.

We denote with the values of the numerical approximation at the nodes and the time . Thus,

(42)

where is the standard base of . Then satisfies the following ODE system:

(43)

where is the mass matrix obtained with lumping, is the stiffness matrix, and is the Lagrange interpolation of the initial datum .

We take and . Writing the system (4.3) explicitly, we get the following ODE system:

(44)

where and . In order to show the evolution in time of a numerical solution, we chose , , and which will be choose later.

First, we consider the case , and the initial data , We observe that the solutions of (1.1) quenching only at the origin, if the symmetric initial data with a unique minimum at (see Figure 1), and the quenching is simultaneous (see Figure 2); If we take , and the same initial data (see Figures 3 and 4), then we obtain the results which accords with Theorem 1.3.

thumbnailFigure 1. The value of the solution at the quenching time .

thumbnailFigure 2. Evolution at the point of the solution .

thumbnailFigure 3. The value of the solution at the quenching time .

thumbnailFigure 4. Evolution at the point of the solution .

Next, we take with the same initial data . In this case the quenching in (1.1) is non-simultaneous with being strictly positive (see Figure 5); If we choose with the initial data (see Figure 6), then we can see that our results coincide with Theorem 1.4.

thumbnailFigure 5. Evolution at the point of the solution .

thumbnailFigure 6. Evolution at the point of the solution .

Finally, we choose In Figure 7, we take the initial data , and in Figure 8 we take the different initial data both equal to , we can see that both non-simultaneous quenching and simultaneous quenching may occur in (1.1), depending on the initial data.

thumbnailFigure 7. Evolution at the point of the solution .

thumbnailFigure 8. Evolution at the point of the solution .

Acknowledgments

This work is supported in part by NSF of China (10771226) and in part by Innovative Talent Training Project, the Third Stage of "211 Project", Chongqing University, Project no.: S-09110.

References

  1. de Pablo, A, Quirós, F, Rossi, JD: Nonsimultaneous quenching. Applied Mathematics Letters. 15(3), 265–269 (2002). Publisher Full Text OpenURL

  2. Zheng, S, Wang, W: Non-simultaneous versus simultaneous quenching in a coupled nonlinear parabolic system. Nonlinear Analysis: Theory, Methods & Applications. 69(7), 2274–2285 (2008). PubMed Abstract | Publisher Full Text OpenURL

  3. Zhou, J, He, Y, Mu, C: Incomplete quenching of heat equations with absorption. Applicable Analysis. 87(5), 523–529 (2008). Publisher Full Text OpenURL

  4. Ferreira, R, de Pablo, A, Quirós, F, Rossi, JD: Non-simultaneous quenching in a system of heat equations coupled at the boundary. Zeitschrift für Angewandte Mathematik und Physik. 57(4), 586–594 (2006). PubMed Abstract | Publisher Full Text OpenURL

  5. Ferreira, R, de Pablo, A, Pérez-Llanos, M, Rossi, JD: Incomplete quenching in a system of heat equations coupled at the boundary. Journal of Mathematical Analysis and Applications. 346(1), 145–154 (2008). Publisher Full Text OpenURL

  6. Ji, R, Zheng, S: Quenching behavior of solutions to heat equations with coupled boundary singularities. Applied Mathematics and Computation. 206(1), 403–412 (2008). Publisher Full Text OpenURL

  7. Zheng, S, Song, XF: Quenching rates for heat equations with coupled singular nonlinear boundary flux. Science in China. Series A. 51(9), 1631–1643 (2008). Publisher Full Text OpenURL

  8. Escobedo, M, Levine, HA: Critical blowup and global existence numbers for a weakly coupled system of reaction-diffusion equations. Archive for Rational Mechanics and Analysis. 129(1), 47–100 (1995). Publisher Full Text OpenURL

  9. Guo, J-S, Sasayama, S, Wang, C-J: Blowup rate estimate for a system of semilinear parabolic equations. Communications on Pure and Applied Analysis. 8(2), 711–718 (2009)

  10. Wang, M: Blow-up rate estimates for semilinear parabolic systems. Journal of Differential Equations. 170(2), 317–324 (2001). Publisher Full Text OpenURL

  11. Chan, CY: Recent advances in quenching phenomena. In: Proceedings of Dynamic Systems and Applications, 1996, Atlanta, Ga, USA. 2, 107–113

  12. Kawarada, H: On solutions of initial-boundary problem for . Publications of the Research Institute for Mathematical Sciences. 10(3), 729–736 (1975)

  13. Salin, T: On quenching with logarithmic singularity. Nonlinear Analysis: Theory, Methods & Applications. 52(1), 261–289 (2003). PubMed Abstract | Publisher Full Text OpenURL

  14. Mu, C, Zhou, S, Liu, D: Quenching for a reaction-diffusion system with logarithmic singularity. Nonlinear Analysis: Theory, Methods & Applications. 71(11), 5599–5605 (2009). PubMed Abstract | Publisher Full Text OpenURL

  15. Pinasco, JP, Rossi, JD: Simultaneous versus non-simultaneous blow-up. New Zealand Journal of Mathematics. 29(1), 55–59 (2000)

  16. Ciarlet, PG: The Finite Element Method for Elliptic Problems, Studies in Mathematics and Its Applications,p. xix+530. North-Holland, Amsterdam, The Netherlands (1978)

  17. Ferreira, R: Numerical quenching for the semilinear heat equation with a singular absorption. Journal of Computational and Applied Mathematics. 228(1), 92–103 (2009). Publisher Full Text OpenURL

  18. Groisman, P, Quirós, F, Rossi, JD: Non-simultaneous blow-up in a numerical approximation of a parabolic system. Computational & Applied Mathematics. 21(3), 813–831 (2002). PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL