Skip to main content

Quenching phenomenon for a non-Newtonian filtration equation with singular boundary flux

Abstract

This paper is concerned with the quenching phenomenon for the one-dimensional non-Newtonian filtration equation with both source term and Neumann boundary condition. With two different kinds of initial data, we prove that the solution must quench in a finite time and the time derivative blows up at a quenching point. The corresponding quenching rate and a lower bound for the quenching time are also obtained.

1 Introduction

In this paper, we study the following problem:

$$\begin{aligned} &u_{t}=\bigl(|u_{x}|^{p-2}u_{x} \bigr)_{x}+(1-u)^{-h}, \quad 0< x< 1, t>0, \end{aligned}$$
(1.1)
$$\begin{aligned} &u_{x}(0,t)=0,\qquad u_{x}(1,t)=-u^{-q}(1,t),\quad t>0, \end{aligned}$$
(1.2)
$$\begin{aligned} &u(x,0)=u_{0}(x),\quad 0\le x\le1, \end{aligned}$$
(1.3)

where \(p\geq2\), h, q are positive constants. \(u_{0}(x):[0,1]\rightarrow(0,1)\) and satisfies some compatibility conditions. Equation (1.1) is known as the classical non-Newtonian filtration equation that incorporates the effects of nonlinear reaction source and nonlinear boundary outflux. The quenching behavior describes the phenomenon that there exists a finite time T such that the solution \(u(x,t)\) of the problem (1.1) satisfy

$$\lim_{t\rightarrow T^{-}}\max\bigl\{ u(x,t):0\le x\le1\bigr\} \rightarrow1 \quad\mbox{or}\quad \lim_{t\rightarrow T^{-}}\min\bigl\{ u(x,t):0\le x\le1\bigr\} \rightarrow0. $$

In 1975, Kawarada [1] first studied the quenching phenomenon for the semilinear heat equation \(u_{t}=u_{xx}+(1-u)^{-1}\). He obtained the results that, when the solution reaches level \(u=1\), the reaction term and the time derivative blow up. Since then, quenching phenomena for semilinear parabolic equations have been studied by many researchers; see for examples [1–6] and the references therein. Quenching phenomenon is dependent on the singular term of the model. Different singular terms cause the problem may occur quenching phenomenon at different levels. Recently, more and more researchers have focused on the quenching phenomenon for parabolic problem with two nonlinear heat sources [7–12]. For example, Chan and Yuen [7] investigated the problem with two nonlinear boundary outfluxes:

$$\begin{aligned} &u_{t}=u_{xx},\quad 0< x< a, 0< t< T,\\ &u_{x}(0,t)=\bigl(1-u(0,t)\bigr)^{-p},\qquad u_{x}(a,t)= \bigl(1-u(a,t)\bigr)^{-q},\quad 0< t< T,\\ &u(x,0)=u_{0}(x), \quad 0\le x\le a, \end{aligned}$$

where \(a, p, q > 0\). The authors proved that the solution quenches only at \(x = a\), meanwhile, the time derivative \(u_{t}\) blows up. Moreover, making use of the positive steady states, they have given criteria for nonquenching and quenching. It is worth mentioning the work by Selcuk and Ozalp [8], who considered the problem

$$\begin{aligned} &u_{t}=u_{xx}+(1-u)^{-p},\quad 0< x< 1, 0< t< T,\\ &u_{x}(0,t)=0,\qquad u_{x}(a,t)=-u(1,t)^{-q},\quad 0< t< T,\\ &u(x,0)=u_{0}(x),\quad 0\le x\le1, \end{aligned}$$

where \(p, q > 0\). They showed that quenching occurs only at \(x = 0\) in finite time and they estimated the bounds of the quenching rate and a lower bound for the quenching time.

On the other hand, as is well known, the singular or degenerate parabolic equations have the property of finite speed of propagation, which are more consistent with biological phenomena in the real world. So, it should be more reasonable to discuss some nonlinear diffusion models. There is a natural question if quenching phenomenon may occur for singular or degenerate equations. Few works were concerned with singular or degenerate parabolic equations [13–19], where only models with one nonlinear source are studied. In [13–15], the authors studied a nonlinear equation with homogeneous boundary conditions. In [16–19], a nonlinear model with nonlinear boundary conditions was discussed. It was Nie et al. [14] who considered quenching for a singular and degenerate quasilinear diffusion equation as follows:

$$\begin{aligned} &x^{q}\frac{\partial u}{\partial t}-\frac{\partial^{2} u^{m}}{\partial x^{2}}=f\bigl(u^{m} \bigr), \quad(x,t)\in(0,a)\times(0,T),\\ &u(0,t)=0=u(a,t),\quad t\in(0,T),\\ &u(x,0)=0,\quad x\in(0,a). \end{aligned}$$

Here \(a>0\), \(q\in\mathbb{R}\), \(m\geq1\). They established the existence of a critical length \(a_{*}\) and proved that the solution exists globally if \(0< a< a_{*}\), while the solution quenches if \(a>a_{*}\). They also investigated the set of the quenching points and the blowing up of \(u_{t}\). Deng and Xu [16] studied nonlinear diffusion equation with a singular boundary condition and investigated finite time quenching for the solution. They also gave the quenching set and quenching rate near the quenching time.

However, as far as we know, there are very few papers concerned with the quenching phenomenon for singular or degenerate parabolic problem with two nonlinear heat sources, even if the linear diffusion equation holds. Obviously, in the model (1.1), the internal source \((1-u)^{-h}\) and the boundary flux \(u^{-q}\) exist, both of which may become singular in some finite time, if the solution reach level \(u=1\) or \(u=0\), respectively. In the present paper, we will discuss these two cases by imposing conditions (A1)-(A4) upon the initial datum, which are give below. First of all, motivated by the work of [8], in Section 2, we will study the quenching phenomenon for the solution reaching the level \(u=1\). We will prove that quenching occurs in finite time under condition (A1) and the only quenching point is \(x = 0\) under conditions (A1) and (A2). Furthermore, \(u_{t}\) blows up at quenching time is discussed. Then the bounds for the quenching rate and the lower bound for the quenching time are estimated. Second, in Section 3, we will do research on the quenching phenomenon for the solution reaching the level \(u=0\) under conditions (A2) and (A4). It will be shown that the solution quenches in finite time and \(u_{t}\) blows up at quenching time at the only quenching point \(x=1\). Finally, we will give bounds on the quenching rate.

Furthermore, in this paper, we need the following hypotheses:

(A1):

\((|u_{0}'(x)|^{p-2}u_{0}'(x))'+(1-u_{0}(x))^{-h}\geq0\);

(A2):

\(u_{0}'(x)\le0\);

(A3):

\(u_{0}'(x)\le-xu_{0}^{-q}(x)\);

(A4):

\((|u_{0}'(x)|^{p-2}u_{0}'(x))'+(1-u_{0}(x))^{-h}\le0\).

2 Quenching phenomenon for the solution reaching the level \(u = 1\)

In this section, we study the quenching phenomenon for the problem (1.1) under the conditions (A1) and (A2). Due to the degeneracy of the equation, the classical solutions might not exist and the weak solution should be considered. However, for simplicity of our arguments, we assume that the solution is appropriately smooth, since we may consider some approximate boundary and initial value conditions.

2.1 Quenching on the boundary and blow-up of \(u_{t}\)

In this section, we prove the solution quenches in finite time and blowing up of \(u_{t}\) at the only quenching point \(x=0\).

Remark 2.1

The assumptions (A1) and (A2) on \(u_{0}(x)\) are proper. For example, for \(p=2\), \(h=9\), and \(q=\log_{30/7}3\), we can choose \(u_{0}(x)=0.9-\frac{3}{2}x^{4.5}\), which satisfies (A1), (A2), and compatibility conditions.

In the following, we discuss the properties of the solution to the problem (1.1).

Lemma 2.1

Assume that (A1), (A2) hold and the solution u of the problem (1.1) exists in \((0,T_{0})\) for some \(T_{0}>0\). Then \(u\in C^{2,1}((0,1]\times(0,T_{0}))\) with \(u_{x}(x,t)<0\) and \(u_{t}(x,t)\geq0\) in \((0,1]\times(0,T_{0})\).

Proof

Let \(v=u_{x}\). Then v satisfies

$$ \begin{aligned} &v_{t}=\bigl(|v|^{p-2}v \bigr)_{xx}+h(1-u)^{-h-1}v,\quad 0< x< 1,0< t< T_{0},\\ &v(0,t)=0,\qquad v(1,t)=-u^{-q}(1,t),\quad 0< t< T_{0},\\ &v(x,0)=u_{0}'(x),\quad 0\le x\le1. \end{aligned} $$
(2.1)

Using the maximum principle, we have \(v<0\), that is, \(u_{x}(x,t)<0\) in \((0,1]\times(0,T_{0})\). Then it is easy to see that the problem (2.1) is nondegenerate in \((0,1]\times(0,T_{0})\). So \(u_{x}\) is a classical solution of (2.1).

On the other hand, setting \(w=u_{t}\), then w solves the following:

$$\begin{aligned} &w_{t}=(p-1) \bigl(|u_{x}|^{p-2}w_{x} \bigr)_{x}+h(1-u)^{-h-1}w, \quad 0< x< 1,0< t< T_{0},\\ &w_{x}(0,t)=0,\qquad w_{x}(1,t)=qu^{-q-1}(1,t)w(1,t),\quad 0< t< T_{0},\\ &w(x,0)=(p-1)\bigl|u_{0}'(x)\bigr|^{p-2}u_{0}''(x)+ \bigl(1-u_{0}(x)\bigr)^{-h},\quad 0\le x\le1. \end{aligned}$$

Utilizing the maximum principle, one shows that \(u_{t}(x,t)\geq0\) in \((0,1]\times(0,T_{0})\). Therefore, the solutions of the problem (1.1) \(u\in C^{2,1}((0,1]\times(0,T_{0}))\) and they satisfy \(u_{x}(x,t)<0\) and \(u_{t}(x,t)\geq0\) in \((0,1]\times(0,T_{0})\). □

Now, we are in a position to show the quenching result.

Theorem 2.1

Assume that (A1) and (A2) hold. Then there exists a finite time T, such that every solution of (1.1) quenches in this time, and the only quenching point is \(x=0\).

Proof

The maximum principle leads to \(0< u(\cdot, t)< 1\) for all t in the existence interval. Taking advantage of the assumption (A1), we have

$$\alpha=-u^{-q(p-1)}(1,0)+ \int_{0}^{1}\bigl(1-u(x,0)\bigr)^{-h}\,dx>0. $$

Denote \(A(t)=\int_{0}^{1} (1-u(x,t))\,dx\). By Lemma 2.1, it is easy to see that

$$\begin{aligned} A'(t)=- \int_{0}^{1}u_{t}(x,t)\,dx =u^{-q(p-1)}(1,t)- \int_{0}^{1}\bigl(1-u(x,t)\bigr)^{-h}\,dx\le- \alpha. \end{aligned}$$

Thus \(A(t)\le A(0)-\omega t\), which means that \(A(t_{0})=0\) for some \(t_{0}>0\). In addition, since \(u_{x}<0\) for \(0< x\le1\), we can see that there exists T (\(0< T< t_{0}\)) such that \(\lim_{t\rightarrow T^{-}}u(0,t)=1\). By means of the singular nonlinearity in the source, u must occur quenching on the boundary \(x=0\). Here and below, we use T to denote the quenching time of the solutions u. In the following, we only need to show that the solutions u cannot take place quenching in \((0,1]\times(\eta, T)\) for some η (\(0<\eta<T\)).

Denote

$$B(x,t)=u_{x}+\varepsilon(b_{2}-x) \quad\mbox{in } (b_{1},b_{2})\times [\eta,T), $$

where \(b_{2}\in(0,1]\), \(b_{1}\in(0,b_{2})\), and ε is a positive constant to be specified later. Since \(u_{x}(x,t)<0\) in \((0,1]\times[0,T)\), \(B(x,t)\) satisfies

$$\begin{aligned} &B_{t}-(p-1)|u_{x}|^{p-2}B_{xx}\\ &\quad=-(p-1) (p-2)|u_{x}|^{p-3}(B_{x}+ \varepsilon)^{2}+h(1-u)^{-h-1}u_{x}< 0, \quad\mbox{for } (x,t)\in(b_{1},b_{2})\times[\eta,T). \end{aligned}$$

Hence, according to the maximum principle, \(B(x,t)\) cannot attain a positive interior maximum. Further, on the parabolic boundary, since \(u_{x}(x,t)<0\) in \((0,1]\times[0,T)\), choosing ε small enough, we have

$$\begin{aligned} &B(b_{1},t)=u_{x}(b_{1},t)+\varepsilon(b_{2}-b_{1})< 0,\quad t\in(\eta,T),\\ &B(b_{2},t)=u_{x}(b_{2},t)< 0,\quad t\in(\eta,T),\\ &B(x,\eta)=u_{x}(x,\eta)+\varepsilon(b_{2}-x),\quad x\in(b_{1},b_{2}). \end{aligned}$$

Making use of the maximum principle, we obtain \(B(x,t)<0\), that is,

$$u_{x}< -\varepsilon(b_{2}-x), \quad(x,t)\in[b_{1},b_{2}] \times[\eta,T). $$

Integrating the above inequality with respect to x from \(b_{1}\) to \(b_{2}\) gives

$$u(b_{2},t)< u(b_{1},t)-\frac{\varepsilon(b_{2}-b_{1})^{2}}{2} < 1- \frac{\varepsilon(b_{2}-b_{1})^{2}}{2}< 1, $$

which implies that \(u(x,t)<1\) if \(0< x\le1\). □

Theorem 2.2

Assume that \(h\geq1\). Then \(u_{t}\) blows up at the quenching point \(x=0\).

Proof

We prove the theorem by contradiction. Assume that \(u_{t}\) is bounded on \([0,1]\times[0,T)\). Then there exists a positive constant M such that \(u_{t}< M\). Thus, we have

$$\bigl(|u_{x}|^{p-2}u_{x}\bigr)_{x}+(1-u)^{-h}< M. $$

Multiplying the above inequality by \(u_{x}\), and integrating with respect to x from 0 to x yield

$$\ln\bigl[1-u(0,t)\bigr]>-\frac{p-1}{p}|u_{x}|^{p}+\ln \bigl[1-u(x,t)\bigr]+M\bigl[u(x,t)-u(0,t)\bigr] $$

for \(h=1\) and

$$\frac{(1-u(0,t))^{-h+1}}{-h+1}>-\frac{p-1}{p}|u_{x}|^{p}+ \frac{(1-u(x,t))^{-h+1}}{-h+1}+M\bigl[u(x,t)-u(0,t)\bigr] $$

for \(h\neq1\). It can be seen that the left-hand side tends to negative infinity as \(t\rightarrow T^{-}\), while the right-hand side is finite. This completes the proof of Theorem 2.2. □

2.2 Quenching rate and lower bound for the quenching time

In this section, a bound on the quenching rate is given and a lower bound for the quenching time is obtained. We present the quenching rate in the following:

Theorem 2.3

Assume that (A1), (A2), and (A3) hold. Then there exists a positive constant \(C_{1}\) such that

$$u(0,t)\geq1-C_{1}(T-t)^{1/(h+1)}, $$

for t sufficiently close to T.

Proof

We define a function \(G(x,t)=|u_{x}(x,t)|^{p-2}u_{x}(x,t)+x^{p-1}u^{-q(p-1)}(x,t)\) in \([0,1]\times[0,T)\). Then \(G(x,t)\) solves

$$\begin{aligned} &G_{t}-(p-1)|u_{x}|^{p-2}G_{xx} \\ &\quad=-(p-1)|u_{x}|^{p-1}h(1-u)^{-h-1}-q(p-1)u^{-q(p-1)-1}x^{p-1}(1-u)^{-h} \\ &\qquad{}-(p-1)^{2}(p-2)|u_{x}|^{p-2}x^{p-3}u^{-q(p-1)} -2q(p-1)^{3}|u_{x}|^{p-1}x^{p-2}u^{-q(p-1)-1} \\ &\qquad{}-(p-1)^{2}q\bigl[q(p-1)+1\bigr]x^{p-1}u^{-q(p-1)-2}|u_{x}|^{p}, \end{aligned}$$

since \(u_{x}<0\), \(G(x,t)\) cannot attain a positive interior maximum. On the other hand, it follows from (A3) that

$$G(x,0)=-\bigl(u_{x}(x,0)\bigr)^{p-1}+x^{p-1}u^{-q(p-1)}(x,0) \le0. $$

Also

$$G(0,t)=0,\qquad G(1,t)=0, $$

for \(t\in(0,T)\). The maximum principle yields \(G(x,t)\le0\) for \((x,t)\in[0,1]\times[0,T)\). Therefore

$$G_{x}(0,t)=\lim_{\sigma\rightarrow0^{+}}\frac{G(\sigma,t)-\Phi(0,t)}{\sigma} =\lim _{\sigma\rightarrow0^{+}}\frac{G(\sigma,t)}{\sigma}\le0. $$

Hence, for \(p=2\),

$$\begin{aligned} G_{x}(0,t)=u_{xx}(0,t)+u^{-q}(0,t)=u_{t}(0,t)- \bigl(1-u(0,t)\bigr)^{-h}+u^{-q}(0,t)\le0 \end{aligned}$$

and for \(p>2\)

$$\begin{aligned} G_{x}(0,t)=\bigl(|u_{x}|^{p-2}u_{x} \bigr)_{x}(0,t)=u_{t}(0,t)-\bigl(1-u(0,t) \bigr)^{-h}\le0. \end{aligned}$$

Thus, we get

$$u_{t}(0,t)\le\bigl(1-u(0,t)\bigr)^{-h}. $$

Integrating for t from t to T, we have

$$u(0,t)\geq1-C_{1}(T-t)^{\frac{1}{h+1}}, $$

where \(C_{1}=(h+1)^{1/(h+1)}\). This completes the proof of Theorem 2.3. □

Remark 2.2

According to Theorem 2.3, a lower bound of quenching time T is \((1-u_{0}(0))^{h+1}/(h+1)\). As in Remark 2.1, if \(u_{0}(x)=0.9-\frac{2}{3}x^{4.5}\), then it can be found that \(T=10^{-11}\) for \(h=9\).

3 Quenching phenomenon for the solution reaching the level \(u = 0\)

In this section, we investigate the quenching phenomenon for the problem (1.1) under the conditions (A2) and (A4).

3.1 Quenching on the boundary and blow-up of \(u_{t}\)

In this section, we prove the solution quenches in finite time and blowing up of \(u_{t}\) at the only quenching point \(x = 1\). First of all, we have the following:

Lemma 3.1

Assume that (A2) and (A4) hold and the solution u of the problem (1.1) exists in \((0,\tilde{T_{0}})\) for some \(\tilde{T_{0}}>0\). Then \(u_{x}(x,t)<0\) and \(u_{t}(x,t)<0\) in \((0,1]\times(0,\tilde{T_{0}})\).

The proof is similar to Lemma 2.1, so we omit it.

Theorem 3.1

Assume that (A2) and (A4) hold. Then there exists a finite time T, such that every solution of (1.1) quenches in this time, and the only quenching point is \(x=1\).

Proof

By the maximum principle, we can obtain \(0< u(\cdot, t)< 1\) for all t in the existence interval. Together with (A4), we get

$$\beta=-u^{-q(p-1)}(1,0)+ \int_{0}^{1}\bigl(1-u(x,0)\bigr)^{-h}\,dx< 0. $$

Denote \(I(t)=\int_{0}^{1} u(x,t)\,dx\). By Lemma 3.1, it is easy to see that

$$\begin{aligned} I'(t)= \int_{0}^{1}u_{t}(x,t)\,dx =-u^{-q(p-1)}(1,t)-+ \int_{0}^{1}\bigl(1-u(x,t)\bigr)^{-h}\,dx\le \beta. \end{aligned}$$

Thus \(I(t)\le I(0)+\beta t\), which means that \(I(\tilde{t_{0}})=0\) for some \(\tilde{t_{0}}>0\). In addition, notice that \(u_{x}<0\) for \(0< x\le1\), we can see that there exists T (\(0< T<\tilde{t_{0}}\)) such that \(\lim_{t\rightarrow T^{-}}u(1,t)=0\). Combining with the singular nonlinearity of the boundary flux, u must occur quenching on the boundary \(x=1\). As in Theorem 2.1, in the following, we only need to show that the solutions u cannot take place quenching in \((1/2,1)\times(\gamma, T)\) for some γ (\(0<\gamma<T\)).

Define

$$H(x,t)=u_{x}+\varepsilon \biggl(x-\frac{1}{4} \biggr), \quad(x,t) \in\biggl(\frac{1}{4},1\biggr)\times (\gamma,T), $$

where ε is sufficiently small. Since \(u_{x}(x,t)<0\) in \((0,1]\times[0,T)\), \(H(x,t)\) satisfies

$$\begin{aligned} &H_{t}-(p-1)|u_{x}|^{p-2}H_{xx} =-(p-1) (p-2)|u_{x}|^{p-3}(H_{x}- \varepsilon)^{2}+h(1-u)^{-h-1}u_{x}< 0, \\ &\quad\mbox{for } (x,t)\in\biggl(\frac{1}{4},1\biggr)\times(\gamma,T). \end{aligned}$$

Further, on the parabolic boundary, since \(u_{x}(x, t) < 0\) in \((0,1]\times[0,T)\) and choosing ε small enough, we have

$$\begin{aligned} &H\biggl(\frac{1}{4},t\biggr)=u_{x}\biggl(\frac{1}{4},t \biggr)< 0, \quad\mbox{for } t\in[\gamma,T), \\ &H(1,t)=-u^{-q}(1,t)+\frac{3}{4}\varepsilon\le-1+ \frac{3}{4}\varepsilon< 0, \quad\mbox{for } t\in[\gamma,T), \\ &H(x,\tilde{\eta})\le u_{x}\biggl(\frac{1}{4},\gamma\biggr)+ \frac{3}{4}\varepsilon< 0, \quad\mbox{for } x\in\biggl[\frac{1}{4},1 \biggr]. \end{aligned}$$

Making use of the maximum principle, we obtain \(H(x,t)\le0\) in \((1/4,1)\times(\gamma,T)\), which yields

$$\begin{aligned} -u_{x}\geq\varepsilon\biggl(x-\frac{1}{4} \biggr),\quad (x,t)\in\biggl(\frac{1}{4},1\biggr)\times(\tilde{\eta},T). \end{aligned}$$
(3.1)

Integrating (3.1) with respect to x from x to 1 gives

$$u(x,t)\geq u(1,t)+ \int_{x}^{1}\varepsilon\biggl(x-\frac{1}{4} \biggr)\,dx \geq \int_{x}^{1}\varepsilon\biggl(x-\frac{1}{4} \biggr)\,dx>0, $$

which implies that \(u(x,t)\) if \(x<1\). This completes the proof of Theorem 3.1. □

Theorem 3.2

\(u_{t}\) blows up at the quenching point \(x=1\).

Proof

We prove the theorem by contradiction. Assume that \(u_{t}\) is bounded on \([0,1]\times[0,T)\). Then there exists a positive constant L such that \(u_{t}>-L\). Thus, we have

$$\bigl(|u_{x}|^{p-2}u_{x}\bigr)_{x}+(1-u)^{-h}>-L. $$

Integrating with respect to x from x to 1 yields

$$-u^{-q(p-1)}(1,t)>-u^{-q(p-1)}(x,t)-L-\bigl(1-u(0,t) \bigr)^{-h}. $$

Therefore, it is found that the left-hand side tends to negative infinity as \(t\rightarrow T^{-}\), while the right-hand side is finite. This completes the proof of Theorem 3.2. □

3.2 Quenching rate

Now, we are in a position to investigate the bounds on the quenching rate. First of all, we will show the lower bound of the quenching rate.

Theorem 3.3

Assume that (A2) and (A4) hold. Then there exists a positive constant \(C_{2}\) such that

$$u(1,t)\geq C_{2}(T-t)^{\frac{1}{pq+2}} $$

for t sufficiently close to T.

Proof

Let \(k(u)=-qu^{-q(p-1)(\delta-1)-q-1}\), where \(1-\frac {q+1}{q(p-1)}<\delta<1-\frac{1}{q(p-1)}\). It is easy to see that \(k(u)<0\), \(k'(u)>0\), and \(k''(u)<0\). Letting Ï„ be close to T, we introduce the function

$$Q(x,t)=u_{t}-\varepsilon k(u) (-u_{x})^{(p-1)(2-\delta)}, \quad\mbox{in } (1-T+\tau,1)\times(\tau,T), $$

where ε is a positive constant. Through a fairly complicated calculation, one has

$$\begin{aligned} Q_{t}=(p-1)|u_{x}|^{p-2}Q_{xx}+(p-1) (p-2) (-u_{x})^{p-3}(-u_{xx})Q_{x} +J(x,t)Q+W(x,t). \end{aligned}$$

Here

$$\begin{aligned} J(x,t)={}&h(1-u)^{-h-1} +\varepsilon(2-\delta)\bigl[(p-1) (2-\delta)-1 \bigr](-u_{x})^{(p-1)(1-\delta)-1}k(u)u_{t} \\ &{}+\varepsilon\bigl[(p-1) (5-2\delta)-1\bigr]k'(u) (-u_{x})^{(p-1)(2-\delta)} \\ &{}+\varepsilon^{2}(2-\delta)\bigl[(p-1) (2-\delta )-1 \bigr]k^{2}(u) (-u_{x})^{(p-1)(3-2\delta)-1} \end{aligned}$$

and

$$\begin{aligned} W(x,t)={}&\varepsilon(p-1)k''(u) (-u_{x})^{(p-1)(3-\delta)+1} +\varepsilon^{2}\bigl[(p-1) (5-2\delta)-1\bigr]k'(u)k(u) (-u_{x})^{2(p-1)(2-\delta)} \\ &{}+\varepsilon^{3}(2-\delta)\bigl[(p-1) (2-\delta )-1 \bigr]k^{3}(u) (-u_{x})^{(p-1)(5-3\delta)-1} \\ &{}-2\varepsilon(2-\delta)\bigl[(p-1) (2-\delta)-1\bigr]k(u) (-u_{x})^{(p-1)(1-\delta )-1}u_{t}(1-u)^{-h} \\ &{}+\varepsilon(2-\delta)\bigl[(p-1) (2-\delta)-1\bigr]k(u) (-u_{x})^{(p-1)(1-\delta )}(1-u)^{-2h} \\ &{}-\varepsilon(p-1)\bigl[(5-2\delta)k'(u) (1-u)+(2-\delta )k(u)h \bigr](-u_{x})^{(p-1)(2-\delta)}(1-u)^{-h-1}. \end{aligned}$$

Notice that \(k(u)<0\), \(k'(u)>0\), \(k''(u)<0\), and Ï„ is sufficiently close to T, then \(J(x,t)>0\) and \(W(x,t)<0\). Therefore,

$$\begin{aligned} &Q_{t}< (p-1)|u_{x}|^{p-2}Q_{xx}+(p-1) (p-2) (-u_{x})^{p-3}(-u_{xx})Q_{x} +J(x,t)Q,\\ &\quad (x,t)\in(1-T+\tau, 1)\times(\tau, T). \end{aligned}$$

Further, on the parabolic boundary, in view of the only quenching point \(x=1\) and provided ε sufficient small, both \(Q(1-T+\tau, t)\) and \(Q(x,\tau)\) are negative. On the right boundary \(x=1\), we get

$$\begin{aligned} Q_{x}(1,t)={}&q\bigl[1-\varepsilon(2-\delta)\bigr]u^{-q-1}(1,t)Q(1,t) +\varepsilon q\bigl\{ \bigl[\varepsilon q(2-\delta)+q(p-1) (\delta-1)+1\bigr] \\ &{}\times u^{-qp-1}(1,t)+(2-\delta) \bigl(1-u(1,t)\bigr)^{-h}\bigr\} u^{-q-1}(1,t)\\ \le{}&q\bigl[1-\varepsilon(2-\delta)\bigr]u^{-q-1}(1,t)Q(1,t), \end{aligned}$$

provided ε is sufficiently small and τ is sufficiently close to T. Thus, take advantage of the maximum principle, \(Q(x,t)\le0\) on \([1-T+\tau,1]\times[\tau,T)\). Then we have \(Q(1,t)\le0\), that is,

$$\begin{aligned} u_{t}(1,t)\le\varepsilon k\bigl(u(1,t)\bigr) \bigl(-u_{x}(1,t)\bigr)^{(p-1)(2-\delta)} =-\varepsilon q u^{-qp-1}(1,t). \end{aligned}$$
(3.2)

Integrating (3.2) with respect to t from t to T, it gives

$$\begin{aligned} u(1,t)\geq\bigl[\varepsilon q(qp+2)\bigr]^{\frac{1}{qp+2}}(T-t)^{\frac{1}{qp+2}} =C_{2}(T-t)^{\frac{1}{qp+2}}, \end{aligned}$$

where \(C_{2}=[\varepsilon q(qp+2)]^{\frac{1}{qp+2}}\). This completes the proof of Theorem 3.3. □

To end this section, we present the upper bound on the quenching rate.

Theorem 3.4

Assume that (A2) and (A4) hold. Then there exists a positive constant \(C_{3}\) such that

$$u(1,t)\le C_{3}(T-t)^{\frac{1}{pq+2}} $$

for t sufficiently close to T.

Proof

Denote \(E(x,t)=|u_{x}(x,t)|^{p-2}u_{x}(x,t)+\varrho^{p-1}(x)u^{-q(p-1)}(x,t)\) in \((0,1)\times(0,T)\), where

$$\begin{aligned} \varrho(x)= \left \{ \textstyle\begin{array}{@{}l@{\quad}l} 0,& x\in[0,x_{0}],\\ \frac{(x-x_{0})^{r}}{(1-x_{0})^{r}}, & x\in(x_{0},1], \end{array}\displaystyle \right . \end{aligned}$$

with some \(x_{0}<1\) and choosing \(r\geq3\) large enough so that \(\varrho(x)\le-u_{0}'(x)u^{q}_{0}(x)\) for \(x_{0}< x\le1\). We can easily obtain \(E(0,t)=E(1,t)=0\), and \(E(x,0)\le0\). In addition, in \((0,1)\times(0,T)\), E satisfies

$$\begin{aligned} E_{t}={}&(p-1)|u_{x}|^{p-2}E_{xx} \\ &{}-(p-1)^{2}\bigl[pq^{2}+q\bigr]\varrho^{p-1}(x)|u_{x}|^{p}u^{-q(p-1)-2} \\ &{}-2(p-1)^{3}q\varrho^{p-2}(x)\varrho'(x)|u_{x}|^{p-1}u^{-q(p-1)-1} \\ &{}-(p-1)^{2}|u_{x}|^{p-2}\varrho^{p-3}(x) \bigl[(p-2)\varrho^{\prime2}(x) +\varrho(x) \varrho''(x)\bigr]u^{-q(p-1)} \\ &{}-q(p-1)\varrho^{p-1}(x)u^{-q(p-1)-1}(1-u)^{-h} -(p-1)h|u_{x}|^{p-1}(1-u)^{-h-1}. \end{aligned}$$

According to the definition of \(\varphi(x)\), it is easy to see that \(\varrho(x)\geq0\), \(\varrho'(x)\geq0\), and \(\varrho''(x)\geq0\). Then we have

$$E_{t}\le(p-1)|u_{x}|^{p-2}E_{xx}. $$

Making use of the maximum principle, we get \(E(x,t)\le0\), that is,

$$\varrho(x)u^{-q}(x,t)\le-u_{x}(x,t), \quad\mbox{for } (x,t) \in[0,1]\times[0,T). $$

Furthermore, because of \(E(x,t)\le0\), we have \(E_{x}(1,t)\geq0\). In fact,

$$\begin{aligned} E_{x}(1,t)=\lim_{x\rightarrow1^{-}}\frac{E(x,t)-E(1,t)}{x-1}\geq0, \end{aligned}$$

which implies

$$\begin{aligned} u_{t}(1,t)&\geq -(p-1)\bigl[ \varrho'(1)+qu^{-q-1}(1,t)\bigr]u^{-q(p-1)}(1,t)+ \bigl(1-u(1,t)\bigr)^{-h} \\ &\geq -\tilde{C}q(p-1)u^{-pq-1}(1,t). \end{aligned}$$
(3.3)

Integrating (3.3) with respect to t from t to T gives

$$u(1,t)\le\bigl[\tilde{C}q(p-1) (pq+2)\bigr]^{\frac{1}{pq+2}}(T-t)^{\frac{1}{pq+2}} =C_{3}(T-t)^{\frac{1}{pq+2}}, $$

where \(C_{3}=[\tilde{C}q(p-1)(pq+2)]^{\frac{1}{pq+2}}\), which produces the asserted result. This completes the proof of Theorem 3.4. □

From Theorem 3.3 and Theorem 3.4, we have the following exact quenching rate.

Corollary 3.1

Assume that (A2) and (A4) hold. Then the solution of the problem (1.1) satisfies

$$C_{2}(T-t)^{\frac{1}{pq+2}}\le u(1,t)\le C_{3}(T-t)^{\frac{1}{pq+2}}, $$

that is,

$$u(1,t)\sim(T-t)^{\frac{1}{pq+2}} $$

for t sufficiently close to T. Here \(C_{2}\), \(C_{3}\) are positive constants which are given in Theorem  3.3 and Theorem  3.4.

References

  1. Kawarada, H: On solutions of initial-boundary problem \(u_{t}=u_{xx}+1/(1-u)\). Publ. Res. Inst. Math. Sci. 10, 729-736 (1975)

    Article  MathSciNet  MATH  Google Scholar 

  2. Fila, M, Levine, HA: Quenching on the boundary. Nonlinear Anal. TMA 21, 795-802 (1993)

    Article  MathSciNet  MATH  Google Scholar 

  3. Levine, HA: The quenching of solutions of linear parabolic and hyperbolic equations with nonlinear boundary conditions. SIAM J. Math. Anal. 14, 1139-1153 (1983)

    Article  MathSciNet  MATH  Google Scholar 

  4. Levine, HA: Quenching, nonquenching, and beyond quenching for solutions of some parabolic equations. Ann. Mat. Pura Appl. 155, 243-260 (1989)

    Article  MathSciNet  MATH  Google Scholar 

  5. Ferreira, R, Pablo, AD, Quirós, F, Rossi, JD: Non-simultaneous quenching in a system of heat equations coupled on the boundary. Z. Angew. Math. Phys. 57, 1-9 (2006)

    Article  MathSciNet  Google Scholar 

  6. Jin, CH, Yin, JX, Zhang, XP: Critical quenching exponents for heat equations coupled with nonlinear boundary flux. Northeast. Math. J. 25, 88-96 (2009)

    MathSciNet  MATH  Google Scholar 

  7. Chan, CY, Yuen, SI: Parabolic problems with nonlinear absorptions and releases at the boundaries. Appl. Math. Comput. 121(2-3), 203-209 (2001)

    Article  MathSciNet  MATH  Google Scholar 

  8. Selcuk, B, Ozalp, N: The quenching behavior of a semilinear heat equation with a singular boundary outflux. Q. Appl. Math. 72(4), 747-752 (2014)

    Article  MathSciNet  MATH  Google Scholar 

  9. Zhi, YH, Mu, CL: The quenching behavior of a nonlinear parabolic equation with nonlinear boundary outflux. Appl. Math. Comput. 184(2), 624-630 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  10. Zhao, CL: Blow-up and quenching for solutions of some parabolic equations. Ph.D. thesis, University of Louisiana, Lafayette (2000)

  11. Ozalp, N, Selcuk, B: Blow up and quenching for a problem with nonlinear boundary conditions. Electron. J. Differ. Equ. 2015, 192 (2015)

    Article  MathSciNet  Google Scholar 

  12. Ozalp, N, Selcuk, B: The quenching behavior of a nonlinear parabolic equation with a singular boundary condition. Hacet. J. Math. Stat. 44(3), 615-621 (2015)

    Google Scholar 

  13. Deng, K: Quenching for solutions of a plasma type equation. Nonlinear Anal. TMA 18, 731-742 (1992)

    Article  MATH  Google Scholar 

  14. Nie, YY, Wang, CP, Zhou, Q: Quenching for singular and degenerate quasilinear diffusion equations. Electron. J. Differ. Equ. 2013, 13 (2013)

    Article  MathSciNet  Google Scholar 

  15. Fila, M, Kawohl, B, Levine, HA: Quenching for quasilinear equations. Commun. Partial Differ. Equ. 17, 593-614 (1992)

    MathSciNet  MATH  Google Scholar 

  16. Deng, K, Xu, M: Quenching for a nonlinear diffusion equation with a singular boundary condition. Z. Angew. Math. Phys. 50, 574-584 (1999)

    Article  MathSciNet  MATH  Google Scholar 

  17. Zhou, MJ, Wang, CP, Nie, YY: Quenching of solutions to a class of semilinear parabolic equations with boundary degeneracy. J. Math. Anal. Appl. 421(1), 59-74 (2015)

    Article  MathSciNet  MATH  Google Scholar 

  18. Yang, Y, Yin, JX, Jin, CH: A quenching phenomenon for one-dimensional p-Laplacian with singular boundary flux. Appl. Math. Lett. 23(9), 955-959 (2010)

    Article  MathSciNet  MATH  Google Scholar 

  19. Yang, Y, Yin, JX, Jin, CH: Quenching phenomenon of positive radial solutions for p-Laplacian with singular boundary flux. J. Dyn. Control Syst. (2015). doi:10.1007/s10883-015-9287-2

    Google Scholar 

Download references

Acknowledgements

This work is partially supported by the National Science Foundation of China (11201311, 11301345), Guangdong Natural Science Foundation (2014A030310074) and Natural Science Foundation of SZU (201425, 201545). The author would like to thank the referees for revising the manuscript and for many valuable comments, which helped to clarify this paper.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ying Yang.

Additional information

Competing interests

The author declares to have no competing interests.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yang, Y. Quenching phenomenon for a non-Newtonian filtration equation with singular boundary flux. Bound Value Probl 2015, 233 (2015). https://doi.org/10.1186/s13661-015-0497-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-015-0497-z

MSC

Keywords