Skip to main content

A doubly degenerate diffusion equation not in divergence form with gradient term

Abstract

In this paper, we investigate positive solutions to the doubly degenerate parabolic equation not in divergence form with gradient term \(u_{t}=u^{m}\operatorname{div}(|\nabla u|^{p-2}\nabla u)+ \lambda u^{q}+ \gamma u^{r}|\nabla u|^{p}\), subject to the null Dirichlet boundary condition. We first establish the local existence of weak solutions to the problem, and then determine in what way the gradient term affects the behavior of solutions. The conditions for global and non-global solutions are obtained with the critical exponent \(r_{c}= \frac{pm-q}{p-1}\). Here we introduce some precise technique for the ‘concavity method’ to deal with the difficult non-divergence form of the model.

1 Introduction

This paper studies the doubly degenerate parabolic equation not in divergence form with gradient term

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} u_{t}=u^{m}\operatorname{div}(|\nabla u|^{p-2}\nabla u)+\lambda u^{q}+\gamma u^{r}|\nabla u|^{p}, &(x,t)\in\Omega\times(0,T), \\ u=0, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x),& x\in\overline{\Omega}, \end{array}\displaystyle \right . $$
(1.1)

where \(\Omega\subset\mathbb{R}^{N}\) is a bounded domain with smooth boundary Ω, \(m\geq1\), \(p>2\), \(q\ge m\), \(r\ge m-1\), \(\gamma>0\), \(\lambda>0 \), and \(u_{0}(x)\in C(\overline{\Omega})\cap W_{0}^{1,p}(\Omega)\), \(u_{0}(x)> 0\) in Ω.

There has been much work contributed to the degenerate parabolic equations not in divergence form. Friedman-McLeod [1] considered the following problem:

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} u_{t}=u^{p}\Delta u+u^{q}, &(x,t)\in \Omega\times(0,T), \\ u=0, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x),& x\in\overline{\Omega}, \end{array}\displaystyle \right . $$
(1.2)

with \(p=2\) and \(q=p+1=3\), for which it was shown that, for sufficiently large domains, the solutions of (1.2) must blow up in finite time regardless of the size of the initial value. The more general situation with \(p>1\) and \(q>1\) was studied by Wiegner [2, 3]. We refer to [4] for more results on (1.2).

Stinner [5] investigated the non-divergence form parabolic equation with gradient term

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} u_{t}=u^{p}\Delta u+ u^{q}+\kappa u^{r}|\nabla u|^{2}, &(x,t)\in\Omega\times(0,T), \\ u=0, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x),& x\in\overline{\Omega}, \end{array}\displaystyle \right . $$
(1.3)

with \(p>0\), \(q>1\), \(r>-1\). Quite differently from the problem (1.2) without gradient term, it was found that the additional gradient term can enforce blow-up in some cases, with the critical exponent \(r_{c}=2p-q\).

Recently, Jin and Yin [6] studied the doubly degenerate diffusion equation

$$ \left \{ \textstyle\begin{array}{l@{\quad}l}u_{t}=u^{m}\operatorname{div}(|\nabla u|^{p-2}\nabla u)+\lambda u^{q}, &(x,t)\in \Omega\times(0,T), \\ u=0, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x),& x\in\overline{\Omega}, \end{array}\displaystyle \right . $$
(1.4)

where \(m\ge1\), \(p>1\), and they obtained the critical exponent \(q_{c}=p+m-1\), namely, the solutions are global if \(q< q_{c}\), and there exist both global and blow-up solutions if \(q>q_{c}\). In the critical case \(q=q_{c}\), blow-up or not of solutions will be determined by the size of the domain.

As for the doubly degenerate diffusion equation with gradient term

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} u_{t}=u\operatorname{div}(|\nabla u|^{p-2}\nabla u)+\gamma |\nabla u|^{p}, &(x,t)\in \Omega\times(0,T), \\ u=0, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x),& x\in\overline{\Omega}, \end{array}\displaystyle \right . $$
(1.5)

Zhou et al. proved the existence conditions of solutions [710].

A natural question is in what way the additional gradient term in (1.1) affects the behavior of solutions. It will be shown that, depending on the complicated interaction among the multi-nonlinearity parameters m, p, q, and r, the problem (1.1) admits the critical exponent \(r_{c}=\frac{pm-q}{p-1}\), for which there is some substantial difficulty to be overcome due to the doubly degenerate diffusion of non-divergence in (1.1). In particular, to treat the critical case \(r=r_{c}\), we will introduce an auxiliary problem \(w_{t}=f(w)(\operatorname{div}(|\nabla w|^{p-2}\nabla w)+cw^{p-1})\) with \(f(w)=h^{m}(w)h'(w)^{p-2}\) and \(h(s)\) solving an ODE problem. We will explore a ‘concavity method’ where some precise technique is necessary to deal with the difficult non-divergence form with the general \(f(w)\).

Throughout the paper, denote by \(\lambda_{1}\) the first Dirichlet eigenvalue of the problem

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} -\operatorname{div}(|\nabla \varphi|^{p-2}\nabla\varphi)=\lambda\varphi^{p-1}, &x\in\Omega, \\ \varphi>0, &x\in\Omega, \\ \varphi=0,& x\in\partial\Omega, \end{array}\displaystyle \right . $$
(1.6)

with the corresponding eigenfunction \(\varphi\in C(\overline{\Omega})\cap W_{0}^{1,p}(\Omega)\), normalized by \(\varphi>0\) in Ω, \(\|\varphi\|_{\infty}=1\).

2 Local existence of weak solutions

We begin with the local existence of solutions to (1.1). Denote

$$\begin{aligned}& \Omega_{T}=\Omega\times(0,T),\qquad S_{T}=\partial\Omega \times(0,T),\qquad \Sigma_{T}=S_{T}\cup\bigl\{ \Omega\times \{0\}\bigr\} , \\& \mathbb{E}=\bigl\{ u\in L^{\infty}(\Omega_{T}); u_{t} \in L^{1}(\Omega _{T}); \nabla u\in L_{\mathrm{loc}}^{p}( \Omega_{T})\bigr\} , \\& \mathbb{E}_{0}=\{u\in\mathbb{E}; u|_{\partial\Omega}=0\}. \end{aligned}$$

The following comparison principle will play a crucial role in the paper, the proof of which can be found in [11].

Lemma 2.1

Let L be the parabolic differential operator defined by

$$Lu := \partial_{t} u-f(x,t,u)\Delta u+g(x,t,u) $$

with continuous functions f and g, \(f\geq0\). Let \(u_{i}\in C^{0}(\overline{\Omega}_{T})\cap C^{2,1}(\Omega_{T})\), \(i=1,2\), be such that \(Lu_{i}\) is well defined in \(\Omega_{T}\) and

$$Lu_{1}\leq Lu_{2} \quad \textit{in } \Omega_{T}, \qquad u_{1}\leq u_{2} \quad \textit{on } \Sigma_{T}. $$

Assume f and g are Lipschitz with respect to u in a neighborhood of \(u_{i}(\overline{\Omega}_{T})\), \(i=1\textit{ or }2\), and in addition either \(u_{1}< u_{2}\) on \(\Sigma_{T}\) or \(\nabla^{2}u_{i}\in L^{\infty}(\overline{\Omega}_{T})\). Then

$$u_{1}\leq u_{2} \quad \textit{in } \overline{ \Omega}_{T}. $$

Since (1.1) degenerates when \(u=0\) or \(|\nabla u|=0\), the problem does not admit classical solutions in general. Here we deal with nonnegative weak solutions, defined as follows.

Definition 2.1

A nonnegative measurable function \(u\in\mathbb{E}\) is called a weak subsolution of problem (1.1)

$$ \int_{0}^{t} \int_{\Omega}\bigl\{ u_{t}\phi+|\nabla u|^{p-2} \nabla u\cdot \nabla\bigl(u^{m}\phi\bigr)-\gamma u^{r}|\nabla u|^{p}\phi\bigr\} \,dx\,d\tau \leq\lambda \int_{0}^{t} \int_{\Omega}u^{q}\phi \,dx\,d\tau $$
(2.1)

for all bounded test functions \(0\leq\phi\in C_{0}^{1}(\Omega_{T})\). The weak supersolution is defined by the opposite inequality, and u is a weak solution of (1.1) if it is both a subsolution and a supersolution to (1.1).

To show the local solvability of (1.1), consider the following regularized problem:

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} (u_{\epsilon\eta})_{t}=u_{\epsilon\eta}^{m}\operatorname{div}\{ (|\nabla u_{\epsilon\eta}|^{2}+\eta)^{\frac{p-2}{2}}\nabla u_{\epsilon\eta}\}+\lambda u_{\epsilon\eta}^{q} \\ \hphantom{(u_{\epsilon\eta})_{t}={}}{} +\gamma u_{\epsilon\eta}^{r}(|\nabla u_{\epsilon\eta}|^{2}+\eta)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}&\text{in } \Omega_{T}, \\ u_{\epsilon\eta}=\epsilon&\text{on } S_{T}, \\ u_{\epsilon\eta}(x,0)=u_{0}(x)+\epsilon& \text{on } \overline{\Omega}, \end{array}\displaystyle \right . $$
(2.2)

where \(\epsilon\in(0,1)\). Denote the classical solution of problem (2.2) by \(u_{\epsilon\eta} \). It is easy to prove for fixed \(\eta\geq0\) that \(u_{\epsilon\eta}\ge\epsilon\), and \(u_{\epsilon\eta}\) is increasing in ϵ.

Lemma 2.2

For any \(\epsilon\in(0,1)\), there exists a function \(u_{\epsilon}\) with \(u_{\epsilon}-\epsilon\in L^{\infty}(\Omega_{T})\cap L^{p}(0,T;W_{0}^{1,p}(\Omega)) \) for some \(T>0\), such that \(u_{\epsilon}\) is a weak solution of the problem

$$ \left \{ \textstyle\begin{array}{l@{\quad}l}u_{t}=u^{m}\operatorname{div}(|\nabla u|^{p-2}\nabla u)+\lambda u^{q}+\gamma u^{r}|\nabla u|^{p}, &(x,t)\in\Omega\times(0,T), \\ u=\epsilon, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x)+\epsilon,& x\in\overline{\Omega}. \end{array}\displaystyle \right . $$
(2.3)

Proof

Step 1. A priori estimates for \(u_{\epsilon\eta}\).

At first it is easy to show that there exist \(T_{1},M_{1}>0\) such that

$$\|u_{\epsilon\eta}\|_{L^{\infty}(\Omega_{T_{1}})}\leq M_{1}\quad \mbox{for all } \eta\geq0. $$

In fact, let U solve

$$\frac{dU}{dt}=\lambda U^{q},\qquad U(0)=\|u_{0}+1 \|_{L^{\infty}(\Omega)} $$

in \([0,T_{0})\). Set \(T=T_{0}/2\). Then

$$ \|u_{\epsilon\eta}\|_{L^{\infty}(\Omega_{T})}\leq M_{1}=U(T)< \infty $$
(2.4)

by comparison.

Choose s satisfying

$$ s>C_{1}:=\left \{ \textstyle\begin{array}{l@{\quad}l} \gamma M_{1}^{r-m+1}, &r>m-1, \\ \gamma-m, &r=m-1. \end{array}\displaystyle \right . $$

Multiply (2.2) by \(u_{\epsilon\eta}^{s}\), and integrate over \(\Omega_{T}\),

$$\begin{aligned}& \int_{\Omega_{T}} (u_{\epsilon\eta})_{t} u_{\epsilon\eta}^{s}\,dx\,dt+(m+s) \int_{\Omega_{T}}u_{\epsilon\eta }^{m+s-1}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta\bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt \\& \quad \le\lambda \int_{\Omega_{T}} (u_{\epsilon\eta})^{q+s}\,dx\,dt+\gamma \int_{\Omega_{T}}u_{\epsilon\eta}^{r+s}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta\bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt. \end{aligned}$$

We have

$$ \int_{\Omega_{T}}u_{\epsilon\eta}^{m+s-1}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta\bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt\le C_{2}=C_{2}(M_{1},s, \gamma), $$
(2.5)

and hence

$$ \int_{\Omega_{T}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt\le C_{3}=C_{3}(M_{1},s,\gamma,\epsilon), $$
(2.6)

due to \(u_{\epsilon\eta}\ge\epsilon\). Consequently,

$$ \int_{\Omega_{T}}|\nabla u_{\epsilon\eta}|^{p}\,dx\,dt\le \int_{\Omega_{T}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt\le C_{3}. $$
(2.7)

Integrate (2.2) over \(\Omega_{T}\). Noticing \(\frac{\partial u_{\epsilon\eta}}{\partial n} |_{\partial\Omega}\le0\), by (2.4) and (2.7), we have

$$\begin{aligned} \int_{\Omega_{T}}(u_{\epsilon\eta})_{t}\,dx\,dt \le&-m \int_{\Omega_{T}} u_{\epsilon\eta}^{m-1}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta\bigr)^{\frac {p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt+\lambda \int_{\Omega_{T}}u_{\epsilon\eta}^{q}\,dx\,dt \\ &{}+\gamma \int_{\Omega_{T}} u_{\epsilon\eta}^{r}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta\bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt \\ \le& C_{4}=C_{4}(C_{3}, M_{1}). \end{aligned}$$
(2.8)

Step 2. Uniform integrability.

For any \(\varepsilon>0\), choose \(\delta=\frac{\varepsilon}{4C_{5}}\) with \(C_{5}\) to be defined. Then for any measurable set \(E\subset\Omega\), \(\operatorname{meas} (E)<\delta\), there exists such that \(E\subset \widetilde{E}\subset\Omega\) with \(\operatorname{meas}(\widetilde{E})<2\delta\). Take \(\rho(x)\in C_{0}^{1}(\Omega)\) satisfying \(\rho(x)=1\) for \(x\in E\), \(\rho(x)=0\) for \(x\in\Omega\setminus\widetilde{E}\), and \(0\le\rho(x)\le1\), \(\nabla\rho(x)\le C\rho^{\alpha}(x)\) in Ω with \(\frac{p-1}{p}<\alpha<1\). Refer to [12] for such ρ. Multiply (2.2) by \(u_{\epsilon\eta}^{s} \rho\) and integrate over \(\widetilde{E}\times[0,T]\) to get

$$\begin{aligned}& \int_{0}^{T} \int_{\widetilde{E}}(u_{\epsilon\eta})_{t}u_{\epsilon \eta}^{s} \rho \,dx\,dt \\& \qquad {}+ \int_{0}^{T} \int_{\widetilde{E}}\bigl[(m+s)u_{\epsilon\eta }^{m+s-1}-\gamma u_{\epsilon\eta}^{r+s}\bigr]\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\rho \,dx\,dt \\& \quad \le\lambda \int_{0}^{T} \int_{\widetilde{E}}u_{\epsilon\eta }^{q+s}\rho \,dx\,dt- \int_{0}^{T} \int_{\widetilde{E}}u_{\epsilon\eta}^{m+s}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta\bigr)^{\frac{p-2}{2}}\nabla u_{\epsilon\eta}\nabla\rho \,dx\,dt \\& \quad \le2\lambda M_{1}^{q+s}T\delta +CM_{1}^{m+s} \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|\rho^{\alpha}\,dx \,dt. \end{aligned}$$
(2.9)

By Young’s inequality,

$$\begin{aligned}& \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|\rho^{\alpha}\,dx \,dt \\& \quad \le \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-1}{2}}\rho^{\alpha}\,dx\,dt \\& \quad \le\sigma \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)\rho \,dx\,dt+c(\sigma) \int_{0}^{T} \int_{\widetilde{E}}\rho^{(\alpha -1)p+1}\,dx\,dt \\& \quad \le\sigma \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\rho \,dx \,dt+\sigma\eta C_{3}+2c(\sigma)\delta. \end{aligned}$$
(2.10)

Choose \(\sigma<\frac{\varepsilon}{2C_{3}\eta}\) satisfying \(\sigma CM_{1}^{m+s}<\epsilon^{m+s-1}(m+s-\gamma M_{1}^{r-m+1})\). We have by (2.9) and (2.10)

$$ \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\rho \,dx \,dt\le2C_{5} \delta+\sigma\eta C_{3}< \varepsilon, $$
(2.11)

where

$$ C_{5}=\frac{\frac{M_{1}^{s+1}}{s+1}+\lambda M_{1}^{q+s}T+CM_{1}^{m+s}c(\sigma)}{\epsilon^{m+s-1}(m+s-\gamma M_{1}^{r-m+1})-\sigma CM_{1}^{m+s}}>0. $$
(2.12)

It follows from (2.11) that

$$ \int_{0}^{T} \int_{E}|\nabla u_{\epsilon\eta}|^{p}\,dx\,dt\le \int_{0}^{T} \int_{\widetilde{E}}\bigl(|\nabla u_{\epsilon\eta}|^{2}+\eta \bigr)^{\frac{p-2}{2}}|\nabla u_{\epsilon\eta}|^{2}\,dx\,dt< \varepsilon. $$
(2.13)

By using a similarly procedure with (2.13), we can obtain

$$ \int_{0}^{T} \int_{E} \bigl\vert (u_{\epsilon\eta})_{t}\bigr\vert \,dx\,dt < \varepsilon. $$
(2.14)

Furthermore, for any fixed \(\zeta>0\), again by (2.13), we have

$$ \int_{0}^{T} \int_{E}u_{\epsilon\eta}^{\zeta}|\nabla u_{\epsilon\eta}|^{p}\,dx\,dt\le M_{1}^{\zeta}\int_{0}^{T} \int_{E}|\nabla u_{\epsilon\eta}|^{p}\,dx\,dt\le M_{1}^{\zeta}\varepsilon. $$
(2.15)

Step 3. Convergence of \(u_{\epsilon\eta}\).

By Dunford-Pettis theorem, we know from the inequalities (2.4), (2.7), (2.8), (2.13), (2.14), and (2.15) that for any \(\epsilon\in(0,1)\), there exist a subsequence of \(u_{\epsilon\eta}\) (denoted still by \(u_{\epsilon\eta}\)) and a function \(u_{\epsilon}\) with \(u_{\epsilon}-\epsilon\in L^{\infty}(\Omega_{T})\cap L^{p}(0,T;W_{0}^{1,p}(\Omega))\), such that as \(\eta\rightarrow0\),

$$\begin{aligned}& u_{\epsilon\eta}\rightarrow u_{\epsilon}\quad \text{a.e. in } \Omega_{T}, \end{aligned}$$
(2.16)
$$\begin{aligned}& \nabla u_{\epsilon\eta}\rightharpoonup\nabla u_{\epsilon}\quad \text{in } L^{p}(\Omega_{T}), \end{aligned}$$
(2.17)
$$\begin{aligned}& \frac{\partial u_{\epsilon\eta}}{\partial t}\rightharpoonup \frac{\partial u_{\epsilon}}{\partial t} \quad \text{in } L^{1}(\Omega_{T}), \end{aligned}$$
(2.18)
$$\begin{aligned}& |\nabla u_{\epsilon\eta}|^{p}\rightharpoonup|\nabla u_{\epsilon}|^{p} \quad \text{in } L^{1}( \Omega_{T}), \end{aligned}$$
(2.19)
$$\begin{aligned}& u_{\epsilon\eta}^{\zeta}|\nabla u_{\epsilon\eta}|^{p} \rightharpoonup u_{\epsilon}^{\zeta}|\nabla u_{\epsilon}|^{p} \quad \text{in } L^{1}(\Omega_{T}), \end{aligned}$$
(2.20)

where denotes weak convergence, and

$$ \epsilon\le u_{\epsilon}\le M_{1} \quad \text{a.e. in } \Omega_{T}. $$
(2.21)

For any \(\phi(x,t)\in C_{0}^{1}(\Omega_{T})\), \(u_{\epsilon\eta}\) satisfies

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\frac{\partial u_{\epsilon\eta }}{\partial t}\phi \,dx\,dt+ \int_{0}^{T} \int_{\Omega}\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta \bigr)^{(p-2)/2}\nabla u_{\epsilon\eta}\cdot \nabla\bigl(u_{\epsilon\eta}^{m} \phi\bigr)\,dx\,dt \\& \quad =\lambda \int_{0}^{T} \int_{\Omega}u_{\epsilon\eta}^{q}\phi \,dx\,dt+\gamma \int_{0}^{T} \int_{\Omega}u_{\epsilon\eta}^{r}\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{(p-2)/2}|\nabla u_{\epsilon\eta }|^{2}\phi \,dx\,dt. \end{aligned}$$
(2.22)

Furthermore,

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\bigl(u^{m}_{\epsilon\eta}\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{(p-2)/2}\nabla u_{\epsilon\eta}-u^{m}_{\epsilon}|\nabla u_{\epsilon}|^{p-2} \nabla u_{\epsilon}\bigr)\cdot\nabla\phi \,dx\,dt \\& \quad = \int_{0}^{T} \int_{\Omega}u^{m}_{\epsilon\eta} \bigl(\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{(p-2)/2}\nabla u_{\epsilon\eta }-|\nabla u_{\epsilon\eta}|^{p-2}\nabla u_{\epsilon\eta} \bigr)\cdot\nabla \phi \,dx\,dt \\& \qquad {} + \int_{0}^{T} \int_{\Omega}u^{m}_{\epsilon\eta} \bigl(|\nabla u_{\epsilon\eta}|^{p-2}\nabla u_{\epsilon\eta}-|\nabla u_{\epsilon}|^{p-2}\nabla u_{\epsilon}\bigr)\cdot\nabla\phi \,dx\,dt \\& \qquad {} + \int_{0}^{T} \int_{\Omega}\bigl(u^{m}_{\epsilon\eta}-u_{\epsilon}^{m} \bigr)|\nabla u_{\epsilon}|^{p-2}\nabla u_{\epsilon}\cdot \nabla\phi \,dx\,dt \\& \quad =I_{1}+I_{2}+I_{3}. \end{aligned}$$

Obviously, as \(\eta\rightarrow0\), we have \(I_{1}\rightarrow0\) due to \(p>2\), \(I_{2}\rightarrow0\) by (2.20), and \(I_{3}\rightarrow0\) because of (2.16) and the dominated convergence theorem. Consequently,

$$\begin{aligned} \begin{aligned}[b] &\int_{0}^{T} \int_{\Omega}\bigl(u^{m}_{\epsilon\eta}\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{(p-2)/2}\nabla u_{\epsilon\eta}-u^{m}_{\epsilon}|\nabla u_{\epsilon}|^{p-2} \nabla u_{\epsilon}\bigr)\cdot\nabla\phi \,dx\,dt \\ &\quad \rightarrow 0, \quad \mbox{as } \eta\rightarrow0. \end{aligned} \end{aligned}$$
(2.23)

By a similar procedure, we have

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\bigl(u_{\epsilon\eta}^{r}\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{(p-2)/2}|\nabla u_{\epsilon\eta}|^{2} -u_{\epsilon}^{r}|\nabla u_{\epsilon}|^{p} \bigr)\phi \,dx\,dt \\& \quad = \int_{0}^{T} \int_{\Omega}\bigl(u_{\epsilon\eta}^{r}\bigl(|\nabla u_{\epsilon\eta}|^{2}+ \eta\bigr)^{(p-2)/2}|\nabla u_{\epsilon\eta}|^{2}- u_{\epsilon\eta}^{r}|\nabla u_{\epsilon\eta}|^{p} \bigr)\phi \,dx\,dt \\& \qquad {}+ \int_{0}^{T} \int_{\Omega}\bigl(u_{\epsilon\eta}^{r}|\nabla u_{\epsilon\eta}|^{p}-u_{\epsilon\eta}^{r}|\nabla u_{\epsilon}|^{p} \bigr)\phi \,dx\,dt \\& \qquad {} + \int_{0}^{T} \int_{\Omega}\bigl(u_{\epsilon\eta}^{r}|\nabla u_{\epsilon}|^{p}-u_{\epsilon}^{r}|\nabla u_{\epsilon}|^{p} \bigr)\phi \,dx\,dt\rightarrow0, \quad \mbox{as } \eta\rightarrow0. \end{aligned}$$
(2.24)

In summary of (2.16), (2.18), (2.22), (2.23), and (2.24), it is true for any \(\phi\in C_{0}^{1}(\Omega_{T})\) that

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\frac{\partial u_{\epsilon}}{\partial t}\phi \,dx\,dt+ \int_{0}^{T} \int_{\Omega}|\nabla u_{\epsilon}|^{p-2}\nabla u_{\epsilon}\cdot\nabla\bigl(u^{m}_{\epsilon}\phi\bigr)\,dx \,dt \\& \quad =\lambda \int_{0}^{T} \int_{\Omega}u_{\epsilon}^{q}\phi \,dx\,dt+\gamma \int_{0}^{T} \int_{\Omega}u_{\epsilon}^{r}|\nabla u_{\epsilon}|^{p}\phi \,dx\,dt. \end{aligned}$$
(2.25)

Together with the initial and boundary conditions, we conclude that \(u_{\epsilon}\) is a weak solution of problem (2.3). □

Theorem 2.1

Let \(u_{0}\in C(\overline{\Omega})\cap W_{0}^{1,p}(\Omega)\) with \(p> 2\). Then the problem (1.1) admits a weak solution \(u\in \mathbb{E}\).

Proof

From (2.21) and the comparison principle, \(u_{\epsilon}\) is bounded and increasing in ϵ. So,

$$ u_{\epsilon}\rightarrow u\quad \mbox{a.e. in } \Omega_{T}, \text{as } \epsilon\rightarrow0. $$
(2.26)

Moreover, the estimate (2.5) is obviously true for \(\eta=0\), namely,

$$ \int_{\Omega_{T}}u_{\epsilon}^{m+s-1}|\nabla u_{\epsilon}|^{p}\,dx\,dt\le C_{2}=C_{2}(M_{1},s, \gamma). $$
(2.27)

On the other hand, set \(\underline{u}:=c_{1}\mathrm{e}^{-\xi t}\varphi(x)\) for \((x,t)\in\Omega\times[0,\infty)\), with φ defined by (1.6) and \(\xi=c_{1}^{m+p-2}\lambda_{1}\). Then

$$ \underline{u}_{t}-\underline{u}^{m}\operatorname{div}\bigl(| \nabla \underline{u}|^{p-2}\nabla\underline{u}\bigr)-\lambda \underline{u}^{q}-\gamma\underline{u}^{r}|\nabla \underline{u}|^{p}\le\underline{u}_{t}-\underline{u}^{m} \operatorname{div}\bigl(|\nabla \underline{u}|^{p-2}\nabla\underline{u} \bigr)\le0. $$
(2.28)

For any \(K\subset\subset\Omega\), let \(c_{1}\) be small such that \(c_{1}\varphi(x)< u_{0}(x)\) on K. By the comparison principle,

$$u_{\epsilon}\ge\underline{u} \quad \mbox{in } K\times[0,\infty). $$

Since \(\varphi\in C^{1}(\overline{\Omega})\) (see [13]) and \(\varphi>0\) in Ω, there exists \(c_{K}\) such that

$$ u_{\epsilon}\ge c_{K} \quad \mbox{in } K\times(0,\infty). $$
(2.29)

Set \(\Omega_{n}=\{x\in\Omega, \operatorname{dist}(x,\partial\Omega)\ge \frac{1}{n}\}\). Then \(\Omega_{n}\subset\subset\Omega\), and \(\Omega_{n}\rightarrow\Omega\) as \(n\rightarrow\infty\). Together with (2.27) and (2.29), we have

$$ \int_{0}^{T} \int_{\Omega_{1}}|\nabla u_{\epsilon}|^{p}\,dx\,dt\le C_{5}=C_{5}(M_{1},s,\gamma). $$
(2.30)

Thus, there exists a subsequence \(\epsilon=\epsilon_{1k}\rightarrow 0\) such that

$$ \nabla u_{\epsilon_{1k}}\rightharpoonup\nabla u \quad \text{in } L^{p} \bigl(\Omega_{1}\times[0,T]\bigr). $$

It is easy to see the inequality (2.30) is still valid for \(u_{\epsilon_{1k}}\) when \(\Omega_{1}\) is replaced by \(\Omega_{2}\). So, there exists a subsequence of \(\epsilon_{2k}\) such that

$$ \nabla u_{\epsilon_{2k}}\rightharpoonup\nabla u \quad \text{in } L^{p} \bigl(\Omega_{2}\times[0,T]\bigr). $$

By induction, we obtain a subsequence \(\epsilon_{nk}\) such that

$$ \int_{0}^{T} \int_{\Omega_{n}}|\nabla u_{\epsilon_{nk}}|^{p}\,dx\,dt\le C_{5}=C_{5}(M_{1},s,\gamma) $$

and

$$ \nabla u_{\epsilon_{nk}}\rightharpoonup\nabla u \quad \text{in } L^{p} \bigl(\Omega_{n}\times[0,T]\bigr). $$

Similar to the above procedure, we see that the estimates (2.8)-(2.15) hold as well for \(u_{\epsilon_{nk}}\) with \(\Omega_{n}\) instead of Ω. Consequently,

$$\begin{aligned}& \frac{\partial u_{\epsilon_{nk}}}{\partial t}\rightharpoonup \frac{\partial u}{\partial t } \quad \text{in } L^{1}\bigl(\Omega_{n}\times[0,T)\bigr), \end{aligned}$$
(2.31)
$$\begin{aligned}& |\nabla u_{\epsilon_{nk}}|^{p}\rightharpoonup|\nabla u|^{p} \quad \text{in } L^{1}\bigl(\Omega_{n}\times[0,T )\bigr), \end{aligned}$$
(2.32)
$$\begin{aligned}& u_{\epsilon_{nk}}^{\zeta}|\nabla u_{\epsilon_{nk}}|^{p} \rightharpoonup u^{\zeta}|\nabla u|^{p} \quad \text{in } L^{1}\bigl(\Omega_{n}\times[0,T)\bigr). \end{aligned}$$
(2.33)

For any \(\phi(x)\in C_{0}^{1}(\Omega_{T})\), there exists \(K\subset\subset \Omega\) such that \(\phi=0\) in \(\Omega\setminus K\). For such K, there exists \(n\in\mathbb{N}\) such that \(K\subset\subset\Omega_{n}\). By (2.33) and the dominated convergence theorem, we obtain

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\bigl(u_{\epsilon_{nk}}^{r}|\nabla u_{\epsilon}|^{p}-u^{r}|\nabla u|^{p} \bigr) \phi \,dx\,dt \\& \quad = \int_{0}^{T} \int_{\Omega_{n}} \bigl(u_{\epsilon_{nk}}^{r}|\nabla u_{\epsilon_{nk}}|^{p}-u^{r}|\nabla u|^{p} \bigr) \phi \,dx\,dt \\& \quad = \int_{0}^{T} \int_{\Omega_{n}} \bigl(u_{\epsilon_{nk}}^{r}|\nabla u_{\epsilon_{nk}}|^{p}-u_{\epsilon_{nk}}^{r}|\nabla u|^{p} \bigr)\phi \,dx\,dt \\& \qquad {}+ \int_{0}^{T} \int_{\Omega_{n}} \bigl(u_{\epsilon_{nk}}^{r}|\nabla u|^{p}-u^{r}|\nabla u|^{p} \bigr)\phi \,dx\,dt \rightarrow0, \quad \mbox{as } \epsilon_{nk}\rightarrow0. \end{aligned}$$
(2.34)

Similarly,

$$ \int_{0}^{T} \int_{\Omega}\bigl( u_{\epsilon_{nk}}^{m-1}|\nabla u_{\epsilon_{nk}}|^{p}-u^{m-1}|\nabla u|^{p} \bigr) \phi \,dx\,dt\rightarrow 0, \quad \mbox{as } \epsilon_{nk} \rightarrow0 $$
(2.35)

and

$$ \int_{0}^{T} \int_{\Omega}\bigl(u_{\epsilon_{nk}}^{m}|\nabla u_{\epsilon_{nk}}|^{p-2}\nabla u_{\epsilon_{nk}}-u^{m}|\nabla u|^{p-2}\nabla u \bigr)\cdot\nabla\phi \,dx\,dt\rightarrow0, \quad \mbox{as } \epsilon_{nk}\rightarrow 0. $$
(2.36)

Notice that \(u_{\epsilon_{nk}}\) satisfies

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\frac{\partial u_{\epsilon _{nk}}}{\partial t}\phi \,dx\,dt+ \int_{0}^{T} \int_{\Omega}u_{\epsilon_{nk}}^{m}|\nabla u_{\epsilon_{nk}}|^{p-2}\nabla u_{\epsilon_{nk}}\cdot\nabla\phi \,dx\,dt \\& \quad =\lambda \int_{0}^{T} \int_{\Omega}u_{\epsilon_{nk}}^{q}\phi \,dx\,dt+\gamma \int_{0}^{T} \int_{\Omega}u_{\epsilon_{nk}}^{r}|\nabla u_{\epsilon_{nk}}|^{p}\phi \,dx\,dt \\& \qquad {}-m \int_{0}^{T} \int_{\Omega}u_{\epsilon_{nk}}^{m-1}|\nabla u_{\epsilon_{nk}}|^{p}\phi \,dx\,dt. \end{aligned}$$
(2.37)

Letting \(\epsilon_{nk}\rightarrow0\) in (2.37), by (2.26)-(2.36), we conclude for the limit function u that

$$\begin{aligned}& \int_{0}^{T} \int_{\Omega}\frac{\partial u}{\partial t}\phi \,dx\,dt+ \int_{0}^{T} \int_{\Omega}|\nabla u|^{p-2}\nabla u\cdot \nabla \bigl(u^{m}\phi\bigr) \,dx\,dt \\& \quad =\lambda \int_{0}^{T} \int_{\Omega}u^{q}\phi \,dx\,dt+\gamma \int_{0}^{T} \int_{\Omega}u^{r}|\nabla u|^{p}\phi \,dx \,dt. \end{aligned}$$
(2.38)

In addition, u satisfies the initial and boundary conditions of (1.1) (in the sense of trace). This proves that u is a weak solution of (1.1). □

We know by the proof of Theorem 2.1 that if \(\|u(\cdot,t)\|_{L^{\infty}(\Omega)}<\infty\), then \(\int_{0}^{T}\int_{K}|\nabla u|^{p} \,dx\,dt<\infty\) for any \(K\subset\subset\Omega\), namely, \(u\in L^{\infty}(\Omega_{T})\) implies \(\nabla u\in L_{\mathrm{loc}}^{p}(\Omega_{T})\). Let \(T^{*}\) be the maximal existence time of the solution u. We get the following proposition immediately.

Proposition 2.1

If \(T^{*}<\infty\), then \(\lim_{t\rightarrow T^{*}}\|u(\cdot,t)\|_{L^{\infty}(\Omega)}=\infty\).

Remark 1

It is mentioned that the uniqueness of such weak solutions to the problem (1.1) cannot be ensured. In the rest of the paper, the solution of (1.1) always means the maximal solution of (1.1), for which the comparison principle is valid.

3 Global existence and nonexistence of solutions

We discuss the existence and nonexistence of global solutions to the problem (1.1) in this section, via a complete classification on the parameters m, p, q as follows:

$$ (\mathrm{a})\ q< p+m-1, \qquad (\mathrm{b})\ q=p+m-1, \qquad (\mathrm{c})\ q>p+m-1. $$

Correspondingly, we have three theorems for them.

Theorem 3.1

Suppose \(q< p+m-1\).

  1. (i)

    If \(r<\frac{pm-q}{p-1}\), then all solutions are global and bounded.

  2. (ii)

    If \(r>\frac{pm-q}{p-1}\), then the solutions blow up for a large domain or large initial data, and they are global for a small domain with small initial data.

  3. (iii)

    If \(r=\frac{pm-q}{p-1}\), the solutions blow up for a large domain and are global for a small domain.

We will prove Theorem 3.1 in five lemmas.

Lemma 3.1

Suppose \(q< p+m-1\) with \(r<\frac{pm-q}{p-1}\). Then all positive solutions of (1.1) are global and bounded.

Proof

Let \(\psi\in C^{1+\beta}(\overline{\Omega})\) solve \(-\operatorname{div}(|\nabla\psi|^{p-2}\nabla\psi)=1\) in Ω with \(\psi|_{\partial\Omega}=0\), \(\beta\in(0,1)\) [13]. The condition \(q< p+m-1\) with \(r<\frac{pm-q}{p-1}\) implies \(m-r>\frac{q-m}{p-1}\). Choose \(\alpha\in(0,1)\) such that \(\frac{q-m}{p-1}<\alpha<m-r\). Let \(w=M+M^{\alpha}\psi(x)\), with \(M\geq 1\) to be determined. Then

$$\begin{aligned}& w_{t}-w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2} \nabla w\bigr)-\lambda w^{q}-\gamma w^{r}|\nabla w|^{p} \\& \quad =w^{m}\bigl[M^{\alpha(p-1)}-\lambda\bigl(M+M^{\alpha}\psi\bigr)^{q-m}-\gamma \bigl(M+M^{\alpha}\psi \bigr)^{r-m}M^{\alpha p}|\nabla\psi|^{p}\bigr] \\& \quad \geq w^{m}\bigl[M^{\alpha(p-1)}-c\bigl(M^{q-m}+M^{\alpha p+r-m} \bigr)\bigr] \end{aligned}$$

with some \(c>0\). Due to the choice of α, we have

$$q-m< \alpha(p-1), \qquad \alpha p+r-m< \alpha(p-1). $$

Now let \(M>1\) be large enough such that \(u_{0}(x)+\epsilon< M\) in Ω and \(M^{\alpha(p-1)}-c(M^{q-m}+M^{\alpha p+r-m})\geq0\). Therefore, the comparison principle yields \(u_{\epsilon}\leq w\) in \(\Omega\times(0,\infty)\). Letting \(\epsilon\rightarrow0\), we obtain \(u\leq w\) in \(\Omega\times(0,\infty)\). □

Lemma 3.2

Suppose \(q< p+m-1\) with \(r>\frac{pm-q}{p-1}\). If Ω contains a ball with radius R large enough, then all solutions of (1.1) blow up in finite time.

Proof

Suppose for contradiction that for any \(R>0\) such Ω admits a global solution u to (1.1) with suitable initial data \(u_{0}\). Without loss of generality, assume \(\overline{B}_{R}(0)\subset\Omega\). We first show that for any fixed \(M>1\), there exists \(t_{0}>0\) such that \(u\ge M\) in \(B_{1}\times(t_{0},\infty)\).

Let \(R'>R\) be such that \(\overline{B}_{R'}(0)\subset\Omega\). Set \(z:=c\phi_{R'}\), with \(c\in(0,\lambda_{R'}^{-\frac{1}{p+m-1-q}})\) small to be determined, \(\lambda_{R'}\) and \(\phi_{R'}\) the first eigenvalue and eigenfunction in the domain \(B_{R'}\), normalized by \(\phi_{R'}>0\) in \(B_{R'}\), \(\|\phi_{R'}\|_{\infty}=1\). Then

$$ z_{t}-z^{m}\operatorname{div}\bigl(|\nabla z|^{p-2} \nabla z\bigr)-\lambda z^{q}-\gamma z^{r}|\nabla z|^{p}\le\lambda_{R'}c^{p+m-1}\phi_{R'}^{p+m-1}-c^{q} \phi_{R'}^{q}\le0 $$

in \(B_{R'}(0)\times(0,\infty)\). Since \(u_{0}>0\) in Ω, \(u_{0}\in C(\overline{\Omega})\), choose c small enough such that \(u_{0}+\epsilon> c\phi_{R'}\) in \(B_{R'}\). By the comparison principle, \(u\ge z=c\phi_{R'}\) in \(B_{R'}\). Furthermore, there exists \(c_{1}>0\) small such that \(u\ge c_{1}\) in \(B_{R}(0)\times(0,\infty)\), due to \(\phi_{R'}>0\) in \(B_{R'}\) and \(\phi_{R'}\in C(\overline{B}_{R'})\) with \(R< R'\). Define

$$v:=y(t) \bigl(\phi_{R}(x)+\delta\bigr), \quad (x,t)\in B_{R}(0)\times(0,\infty), $$

where \(y(t)\) is a nondecreasing positive function on \([0,\infty)\) to be determined with suitable \(\delta>0\). A direct computation yields

$$\begin{aligned} \begin{aligned} &v_{t}-v^{m}\operatorname{div}\bigl(|\nabla v|^{p-2} \nabla v\bigr)-\lambda v^{q}-\gamma v^{r}|\nabla v|^{p} \\ &\quad \le y'\bigl(\phi_{R}(x)+\delta\bigr)+ \lambda_{R} y^{p+m-1}(\phi_{R}+ \delta)^{m}\phi_{R}^{p-1}-y^{q}( \phi_{R}+\delta)^{q} \\ &\quad = :L_{1}+L_{2}-L_{3} \quad \mbox{in } B_{R}(0)\times(0,\infty), \end{aligned} \end{aligned}$$

where

$$\frac{L_{1}}{L_{3}}=\frac{y'}{y^{q}}(\phi_{R}+\delta)^{1-q} \le\frac{y'}{y^{q}}\delta^{1-q} $$

and

$$\frac{L_{2}}{L_{3}}=\lambda_{R} y^{p+m-1-q}( \phi_{R}+\delta)^{m-q}\phi _{R}^{p-1}\le \lambda_{R} (2y)^{p+m-1-q}. $$

Let \(y_{R}=\frac{1}{2(2\lambda_{R})^{\frac{1}{p+m-1-q}}}\), and choose \(y_{0}\) and δ such that \(y_{0}\le\frac{c_{1}}{2}\), \(\delta\le \min \{\frac{c_{1}}{y_{R}},1 \}\), and \(y'\le \frac{\delta^{q-1}}{2}y^{q}\), \(y(0)=y_{0}\), and \(y(t)\nearrow y_{R}\) as \(t\rightarrow\infty\). We have \(v\le u_{\epsilon}\) on \(\{t=0\}\) and \(\partial B_{R}(0)\). Moreover, \(\frac{L_{1}}{L_{3}}\le\frac{1}{2}\), \(\frac{L_{2}}{L_{3}}\le\frac{1}{2} \). By comparison, \(u_{\epsilon}\ge v\), and hence \(u\ge v\) in \(B_{R}(0)\times(0,\infty)\).

For any fixed \(M>1\), choose \(R>1\) large such that \(\frac{1}{2}y_{R}\phi_{R}\ge M\) in \(B_{1}\), and hence \(y(t_{0})\ge \frac{1}{2}y_{R}\phi_{R}\) for some \(t_{0}>0\). Consequently,

$$ u\ge v\ge\frac{1}{2}y_{R}\phi_{R}\ge M \quad \text{in } B_{1}\times(t_{0},\infty). $$
(3.1)

For \(n\in\mathbb{N}\) with \(n>\max\{p,\frac{pq}{r(p-1)-pm+q}\}\), choose \(0\le\theta\in C_{0}^{1}(\overline{B}_{1})\) such that \(\int_{B_{1}}\theta^{n}\,dx=1\). Multiply (1.1) by \(\theta^{n}\) and integrate over \(B_{1}\),

$$\begin{aligned} \frac{d}{dt} \int_{B_{1}} u\theta^{n} \,dx =&-m \int_{B_{1}} u^{m-1}\theta^{n}|\nabla u|^{p} \,dx-n \int_{B_{1}} u^{m}\theta^{n-1}|\nabla u|^{p-1}\nabla u\cdot\nabla\theta \,dx \\ &{} +\lambda \int_{B_{1}} u^{q}\theta^{n}\,dx+\gamma \int_{B_{1}} u^{r}\theta^{n}|\nabla u|^{p} \,dx \\ =:&-I_{1}-I_{2}+I_{3}+I_{4}. \end{aligned}$$
(3.2)

Since \(r>\frac{pm-q}{p-1}>\frac{pm-p-m+1}{p-1}=m-1\), by (3.1), we have

$$ I_{1}\leq mM^{-(r+1-m)} \int_{B_{1}}u^{r}\theta^{n}|\nabla u|^{p}\,dx\leq \frac{1}{2}I_{4}, \quad t>t_{0}, $$
(3.3)

provided \(M\ge (\frac{2m}{\gamma} )^{\frac{1}{r+1-m}}\). By Young’s inequality,

$$\begin{aligned} I_{2}&\leq n\frac{\gamma}{2n} \int_{B_{1}}u^{r}\theta^{n}|\nabla u|^{p}\,dx+n\frac{2n}{\gamma} \int_{B_{1}}u^{pm-r(p-1)}\theta^{n-p}|\nabla \theta|^{p}\,dx \\ &\leq \frac{1}{2}I_{4}+\frac{2n^{2}}{\gamma} \int_{B_{1}}u^{pm-r(p-1)}\theta ^{n-p}|\nabla \theta|^{p}\,dx. \end{aligned}$$
(3.4)

If \(r\geq\frac{pm}{p-1}\), by (3.1), (3.4), with \(M>1\), we have

$$\begin{aligned} I_{2}&\leq \frac{1}{2}I_{4}+\frac{2n^{2}}{\gamma}M^{pm-r(p-1)} \int_{B_{1}}\theta ^{n-p}|\nabla \theta|^{p} \,dx \\ &\leq \frac{1}{2}I_{4}+\frac{2n^{2}}{\gamma} \int_{B_{1}}\theta^{n-p}|\nabla \theta|^{p} \,dx, \quad t>t_{0}. \end{aligned}$$
(3.5)

In the case \(r<\frac{pm}{p-1}\), due to \(r>\frac{pm-q}{p-1}\) implying \(\beta:= \frac{q}{pm-r(p-1)}>1\), again by Young’s inequality,

$$\begin{aligned}& \int_{B_{1}}u^{pm-r(p-1)}\theta^{n-p}|\nabla \theta|^{p}\,dx \\& \quad \leq \frac{\lambda\gamma}{4n^{2}} \int_{B_{1}} u^{q}\theta^{n}\,dx+ \frac{1}{\beta'} \biggl(\frac{4n^{2}}{\lambda\gamma \beta} \biggr)^{\frac{\beta'}{\beta}} \int_{B_{1}} \theta^{n-p\beta'}|\nabla\theta|^{p\beta'} \,dx \end{aligned}$$
(3.6)

with \(\beta'=\frac{\beta}{\beta-1}=\frac{q}{r(p-1)-pm+q}\).

It follows from (3.4)-(3.6) that

$$\begin{aligned} I_{2} \leq&\frac{1}{2}I_{3}+\frac{1}{2}I_{4}+ \frac{\lambda(r(p-1)-pm+q)}{2(pm-r(p-1))}\cdot \biggl(\frac {4n^{2}(pm-r(p-1))}{q\gamma\lambda} \biggr)^{\frac{q}{r(p-1)-pm+q}} \\ &{}\cdot \int_{B_{1}} \theta^{n-\frac{pq}{r(p-1)-pm+q}}|\nabla\theta|^{\frac {pq}{r(p-1)-pm+q}} \,dx. \end{aligned}$$
(3.7)

Combine (3.2), (3.3), and (3.7) to get

$$ \frac{d}{dt} \int_{B_{1}} u\theta^{n} \,dx\geq\frac{\lambda}{2} \int_{B_{1}} u^{q}\theta^{n}\,dx-C\geq \frac{\lambda}{2} \biggl( \int_{B_{1}} u\theta^{n}\,dx \biggr)^{q} -C $$
(3.8)

with

$$ C=\left \{ \textstyle\begin{array}{l@{\quad}l} \frac{\lambda(r(p-1)-pm+q)}{2(pm-r(p-1))}\cdot (\frac{4n^{2}(pm-r(p-1))}{q\gamma\lambda} )^{\frac{q}{r(p-1)-pm+q}} \\ \quad {}\cdot\int_{B_{1}} \theta^{n-\frac{pq}{r(p-1)-pm+q}}|\nabla\theta|^{\frac{pq}{r(p-1)-pm+q}}\,dx, &r< \frac{pm}{p-1}, \\ \frac{2n^{2}}{\gamma}\int_{B_{1}}\theta^{n-p}|\nabla\theta|^{p}\,dx,& r\geq\frac{pm}{p-1}. \end{array}\displaystyle \right . $$

Choose \(M\ge\max\{ (\frac{2m}{\gamma} )^{\frac{1}{r+1-m}}, (\frac {2C}{\lambda} )^{\frac{1}{q}}\}\). By (3.1), we have

$$ \frac{\lambda}{2} \biggl( \int_{B_{1}} u(\cdot,t_{0})\theta^{n}\,dx \biggr)^{q}\geq\frac{\lambda}{2} \biggl(M \int_{B_{1}} \theta^{n}\,dx \biggr)^{q}= \frac{\lambda}{2}M^{q}>C. $$
(3.9)

We conclude from (3.8) and (3.9) that \(\int_{B_{1}} u\theta^{n} \,dx\) must blow up in finite time. □

Lemma 3.3

Suppose \(q\leq p+m-1\) with \(r>\frac{pm-q}{p-1}\). Let \(u_{0}=b\omega\) with \(\omega\in W_{0}^{1,p}(\Omega)\cap C(\overline{\Omega})\) and \(\omega>0\) in Ω. Then there is \(b_{0}>0\) such that the solution u of (1.1) blows up in finite time provided \(b\geq b_{0}\).

Proof

Suppose for contradiction that u is global with \(u_{0}\ge b_{0} \omega\) for any \(b_{0}>0\). Pick n and M defined in the proof of Lemma 3.2. Let Ω contain a ball with radius 2R, without loss generality, \(B_{2R}(0)\subset\Omega\). It suffices to show there exists \(t_{0}>0\) such that \(u\ge M\) in \(B_{R/2}\times(t_{0},\infty)\).

For \(\sigma>\max \{\frac{p}{p-1},\frac{p(r+1-m)}{(p-1)r+q-pm} \}\), set \(w(x):=\delta\mathrm{e}^{-z}\) for \(x\in B_{R}(0)\) with \(z:=\frac{|x|^{\sigma}}{R^{2}-|x|^{2}}\) and \(\delta>0\) to be determined. For \(x\in B_{R}(0)\), we have

$$\begin{aligned}& w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2}\nabla w\bigr)+ \lambda w^{q}+\gamma w^{r}|\nabla w|^{p} \\& \quad =\delta^{p+m-1}\mathrm{e}^{-(p+m-1)z}\frac{|x|^{\sigma p-\sigma-p}(\sigma R^{2}-(\sigma-2)|x|^{2})^{p-2}}{(R^{2}-|x|^{2})^{2p}} \bigl\{ (p-1)|x|^{\sigma}\bigl(\sigma R^{2}-(\sigma-2)|x|^{2} \bigr)^{2} \\& \qquad {}-n\bigl(R^{2}-|x|^{2}\bigr)^{2}\bigl( \sigma R^{2}-(\sigma-2)|x|^{2}\bigr)-4|x|^{2} \bigl(R^{2}-|x|^{2}\bigr) \bigl(\sigma R^{2}-( \sigma-2)|x|^{2}\bigr) \\& \qquad {}-\sigma(\sigma p-\sigma-p) \bigl(R^{2}-|x|^{2} \bigr)^{3}-(2p-4)|x|^{2}\bigl(R^{2}-|x|^{2} \bigr)^{2} \bigr\} \\& \qquad {}+\lambda\delta^{q}\mathrm{e}^{-qz}+\gamma \delta^{r+p}\mathrm{e}^{-(r+p)z}\frac{|x|^{p(\sigma-1)}(\sigma R^{2}-(\sigma-2)|x|^{2})^{p}}{(R^{2}-|x|^{2})^{2p}}. \end{aligned}$$

Denote

$$\begin{aligned} f(s) :=&(p-1)s^{\sigma}\bigl(\sigma R^{2}-( \sigma-2)s^{2}\bigr)^{2} -n\bigl(R^{2}-s^{2} \bigr)^{2}\bigl(\sigma R^{2}-(\sigma-2)s^{2}\bigr) \\ &{}-4s^{2}\bigl(R^{2}-s^{2}\bigr) \bigl(\sigma R^{2}-(\sigma-2)s^{2}\bigr)2 \\ &{}-\sigma(\sigma p-\sigma-p) \bigl(R^{2}-s^{2} \bigr)^{3}-(2p-4)s^{2}\bigl(R^{2}-s^{2} \bigr)^{2}, \quad s\in[0,R]. \end{aligned}$$

Then there exists \(K>0\) such that \(|f(s)|\leq K\) for \(s\in[0,R]\). Due to \(f(R)=(p-1)2^{p}R^{\sigma+p}>0\), there exists \(c_{2}\in(\frac{1}{2},1)\) such that \(f(s)\geq0\) in \([c_{2}R,R]\), and hence

$$w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2}\nabla w\bigr)+ \lambda w^{q}+\gamma w^{r}|\nabla w|^{p}\geq0,\quad |x|\in[c_{2}R,R]. $$

For the above σ with \(r>\frac{pm-q}{p-1}\), set

$$l\in \biggl(\frac{p+m-1-q}{\sigma p-\sigma-p},\frac{r+1-m}{\sigma} \biggr). $$

Then there exists \(\delta_{0}\geq1\) such that \(c_{1}=\delta^{-l}\in(0,\frac{1}{2})\) for all \(\delta\geq\delta_{0}\). Hence, there exists \(\delta_{1}\geq \delta_{0}\) such that

$$\begin{aligned}& w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2}\nabla w\bigr)+ \lambda w^{q}+\gamma w^{r}|\nabla w|^{p} \\& \quad \geq -\delta^{p+m-1}\frac{c_{1}^{\sigma p-\sigma-p}R^{\sigma p-\sigma-5p}K}{(1-c_{1}^{2})^{2p}}+\lambda\delta^{q} \mathrm{e}^{-\frac{qc_{1}^{\sigma}R^{\sigma-2}}{1-c_{1}^{2}}} \\& \quad \geq-\delta^{p+m-1-l(\sigma p-\sigma-p)}\frac{R^{\sigma p-\sigma-5p}K}{(\frac{3}{4})^{2p}}+\lambda\delta^{q} \mathrm{e}^{-\frac{q}{3}(\frac{R}{2})^{\sigma-2}} \\& \quad \geq 0, \quad |x|\in[0,c_{1}R], \delta\geq\delta_{1}. \end{aligned}$$

Furthermore, there is \(\delta_{2}\geq\delta_{1}\) such that

$$\begin{aligned}& w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2}\nabla w\bigr)+ \lambda w^{q}+\gamma w^{r}|\nabla w|^{p} \\& \quad \geq \frac{|x|^{\sigma p-\sigma-p}(\sigma R^{2}-(\sigma-2)|x|^{2})^{p-2}}{(R^{2}-|x|^{2})^{2p}} \bigl\{ -\delta^{p+m-1}K +4\gamma \delta^{r+p}c_{1}^{\sigma}R^{\sigma+4} \mathrm{e}^{-(r+p)\frac {(c_{2}R)^{\sigma}}{(1-c_{2}^{2})R^{2}}} \bigr\} \\& \quad \geq\frac{|x|^{\sigma p-\sigma-p}(\sigma R^{2}-(\sigma-2)|x|^{2})^{p-2}}{(R^{2}-|x|^{2})^{2p}} \bigl\{ -\delta^{p+m-1}K +4\gamma \delta^{r+p-l\sigma}R^{\sigma+4} \mathrm{e}^{-(r+p)\frac{(c_{2}R)^{\sigma}}{(1-c_{2}^{2})R^{2}}} \bigr\} \\& \quad \geq0, \quad |x|\in[c_{1}R,c_{2}R], \delta\geq \delta_{2}. \end{aligned}$$

Finally, choose \(\delta\ge\delta_{2}\) such that \(w(x)\ge M\) in \(B_{R/2}\). In summary, we obtain

$$w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2}\nabla w\bigr)+ \lambda w^{q}+\gamma w^{r}|\nabla w|^{p}\ge0\quad \text{in } B_{R}(0), $$

and \(w(x)\ge M\) in \(B_{R/2}\).

Let \(b_{0}>0\) be large such that \(b_{0}\omega\ge w\) in \(B_{R}(0)\). By the comparison principle, we have \(u_{\epsilon}\ge w\) in \(B_{R}(0)\times(0,\infty)\). Hence \(u(x,t)\ge w(x)\ge M\) in \(B_{R/2}\times(0,\infty)\).

The rest of the lemma can be proved by the same arguments as those for Lemma 3.2. □

Lemma 3.4

Suppose \(q< p+m-1\) with \(r\ge\frac{pm-q}{p-1}\). Then all positive solutions of (1.1) are bounded, if Ω is contained in a ball with radius R small.

Proof

Without loss of generality, assume \(\Omega\subset \{x\in\mathbb{R}^{N}|R< x_{1}<2R\}\). Set \(z:= Kx_{1}^{\kappa}\) for \((x,t)\in \Omega\times(0,\infty)\), with \(K>0\) and \(0<\kappa<1\) to be determined. A simple calculation yields

$$\begin{aligned}& z_{t}-z^{m}\operatorname{div}\bigl(|\nabla z|^{p-2} \nabla z\bigr)-\lambda z^{q}-\gamma z^{r}|\nabla z|^{p} \\& \quad =(p-1) (1-\kappa)K^{p+m-1}\kappa^{p-1}x_{1}^{(p-1)(\kappa-1)-1+m\kappa }- \lambda K^{q}x_{1}^{q\kappa}-\gamma K^{r+p} \kappa^{p} x_{1}^{p(\kappa-1)+r\kappa } \\& \quad =:L_{1}-L_{2}-L_{3}. \end{aligned}$$

Noticing \(r\ge\frac{pm-q}{p-1}>m-1\), we have

$$\begin{aligned} \frac{2L_{3}}{L_{1}}&=\frac{2\gamma\kappa}{(p-1)(1-\kappa )}K^{r-m+1}x_{1}^{\kappa(r-m+1)} \\ &\le \frac{2\gamma\kappa}{(p-1)(1-\kappa)}K^{r-m+1}(2R)^{\kappa (r-m+1)} \\ &\le1, \end{aligned}$$

provided \(K= (\frac{(p-1)(1-\kappa)}{2\gamma\kappa} )^{\frac {1}{r-m+1}}(2R)^{-\kappa}\).

Choose \(\kappa= \min\{\frac{p-1}{4\gamma M^{r-m+1}+p-1}, \frac{p}{p+m-1-q}\}\), with \(M:=\|u_{0}\|_{L^{\infty}(\Omega)}+1\). Then there exists R small such that

$$K\ge \biggl(\frac{2}{(p-1)(1-\kappa)\kappa^{p-1}} \biggr)^{\frac {1}{p+m-1-q}}(2R)^{\frac{p}{p+m-1-q}-\kappa}, $$

and hence

$$\begin{aligned} \frac{2L_{2}}{L_{1}}&=\frac{2}{(p-1)(1-\kappa)\kappa ^{p-1}}K^{q-(p+m-1)}x_{1}^{q\kappa-m\kappa-(p-1)(\kappa-1)+1} \\ &\le \frac{2}{(p-1)(1-\kappa)\kappa^{p-1}}K^{q-(p+m-1)}(2R)^{p-\kappa (p+m-1-q)} \\ &\le1. \end{aligned}$$

In addition,

$$z=Kx_{1}^{\kappa}\ge \biggl(\frac{(p-1)(1-\kappa)}{2\gamma\kappa} \biggr)^{\frac {1}{r-m+1}}(2R)^{-\kappa}R^{\kappa}\ge M $$

on the parabolic boundary of \(\Omega\times(0,\infty)\) due to the choice of κ. We conclude that z is a time-independent supersolution of (1.1). □

Lemma 3.5

Suppose \(q< p+m-1\) with \(r=\frac{pm-q}{p-1}\). Then all solutions of (1.1) blow up in finite time provided Ω large enough.

Proof

The lemma will be proved in three steps.

Step 1. Set \(u=h(v)\), and substitute into (1.1),

$$\begin{aligned} h'(v)v_{t} =& \bigl[h^{m}(v)h'(v)^{p-2} \bigr]\operatorname{div}\bigl(|\nabla v|^{p-2}\nabla v\bigr)+\lambda h^{q}(v) \\ &{}+ \bigl[(p-1)h(v)^{m} h'(v)^{p-2}h''(v)+ \gamma h^{r}(v)h'(v)^{p} \bigr]|\nabla v|^{p}. \end{aligned}$$
(3.10)

Here \(h\in C^{0}([0,+\infty))\cup C^{2}((0,+\infty))\) satisfies

$$h'(s)=\mathrm{e}^{-\frac{\gamma}{\beta(p-1)}h^{\beta}(s)}, \qquad h(0)=0 $$

with \(\beta=r-m+1=\frac{p+m-1-q}{p-1}\in(0,1]\), and hence \((p-1)h(v)^{m} h''(v)+\gamma h^{r}(v)h'(v)^{2}=0\). Thus,

$$ v_{t}=h^{m}(v)h'(v)^{p-2} \bigl( \operatorname{div}\bigl(|\nabla v|^{p-2}\nabla v\bigr)+\gamma h^{(p-1)(1-\beta)}(v)\mathrm{e}^{\frac{\lambda}{\beta}h^{\beta}(s)} \bigr). $$

Set \(g(v)=h^{(p-1)(1-\beta)}(v)\mathrm{e}^{\frac{\lambda}{\beta}h^{\beta}(s)}\). Similarly to the proof of Theorem 2.5 in [5], we can find constants \(v_{0}\ge2\) and \(c_{0}>0\), only depending on β, λ, and γ, such that

$$ g(v)\ge c_{0} (v+1)^{p-1} \quad \text{for } v\ge v_{0}. $$
(3.11)

So, v is a supersolution of

$$ w_{t}=h^{m}(w)h'(w)^{p-2} \bigl( \operatorname{div}\bigl(|\nabla w|^{p-2}\nabla w\bigr)+c_{0}w^{p-1} \bigr) \quad \text{in } \Omega\times(t_{0},T) $$

whenever \(v\ge v_{0}\).

Step 2. Let w solve

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} w_{t}=f(w) (\operatorname{div}(|\nabla w|^{p-2}\nabla w)+c_{0}w^{p-1} ), &(x,t)\in \Omega\times(0,T), \\ w(x,t)=0, &(x,t)\in\partial\Omega\times(0,T), \\ w(x,0)=w_{0}(x), &x\in\Omega, \end{array}\displaystyle \right . $$
(3.12)

with \(f(w)=h^{m}(w)h'(w)^{p-2}\). We claim that w blows up in finite time for any initial data \(w_{0}>0\) provided Ω large such that \(\lambda_{1}< c_{0}\).

We prove the claim by using the so-called ‘concavity’ method. Define

$$ \mathscr{E}(w)=\frac{1}{p} \int_{\Omega}|\nabla w|^{p} \,dx-\frac{c_{0}}{p} \int_{\Omega}w^{p}\,dx,\qquad \mathscr{H}(w)= \int_{\Omega}\Phi(w)\,dx, $$

where \(\Phi(w)=\int_{0}^{w}\frac{\varrho}{f(\varrho)}\, d\varrho\) with \(f(\varrho)=h^{m}(\varrho)h'(\varrho)^{p-2}\). Such Φ is well defined. In fact, since \(h(0)=0\), there exists \(\varrho_{1}>0\) such that \(h'(\varrho)=\mathrm{e}^{-\frac{\gamma}{\beta(p-1)}h^{\beta}(\varrho)}\ge\frac{1}{2}\) for \(\varrho\in[0,\varrho_{1}]\). Set \(h_{1}(\varrho)= (\alpha\ln (\varrho+1) )^{\frac{1}{\beta}}\), with \(0<\alpha<\frac{\beta(p-1)}{\gamma}\). Then

$$ \frac{h_{1}'(\varrho)}{\mathrm{e}^{-\frac{\gamma}{\beta(p-1)}h_{1}^{\beta}(\varrho)}}=\frac{\frac {1}{\beta} (\alpha\ln(\varrho+1) )^{\frac{1-\beta}{\beta}}\frac {1}{\varrho+1}}{ (\varrho+1)^{-\frac{\alpha\gamma}{\beta(p-1)}}}=\frac{1}{\beta} \bigl(\alpha\ln( \varrho+1) \bigr)^{\frac{1-\beta}{\beta}}(\varrho +1)^{\frac{\alpha\gamma}{\beta(p-1)}-1}. $$

Choose \(\varrho_{2}=\min\{\varrho_{1},\mathrm{e}^{\frac{\beta^{\frac{\beta}{1-\beta}}}{\alpha}-1}\}\). We have \((\alpha\ln(\varrho+1) )^{\frac{1-\beta}{\beta}}\le \beta\) for \(0\le\varrho\le\varrho_{2}\). Since \(0<\alpha<\frac{\beta(p-1)}{\gamma}\) implies \((\varrho+1)^{\frac{\alpha\gamma}{\beta(p-1)}-1}\le1\) for \(0\le \varrho\le\varrho_{2}\), we have

$$ \frac{h_{1}'(\varrho)}{\mathrm{e}^{-\frac{\gamma}{\beta(p-1)}h_{1}^{\beta}(\varrho)}}\le1 \quad \mbox{for } \varrho\in[0,\varrho_{2}]. $$

Hence, by the comparison principle with \(h(0)=h_{1}(0)=0\), we obtain

$$ h(\varrho)\ge h_{1}(\varrho) \quad \mbox{for } \varrho\in[0, \varrho_{2}]. $$

Consequently,

$$ \int_{0}^{\varrho_{2}}\frac{\varrho}{f(\varrho)}\, d\varrho\le \int_{0}^{\varrho_{2}}\frac{\varrho}{2h_{1}^{m}(\varrho)}\, d\varrho\le \int_{0}^{\varrho_{2}}\frac{\varrho}{2\alpha^{\frac{m}{\beta}}\ln ^{\frac{m}{\beta}}(\varrho+1)}\, d\varrho< \infty, $$

and so, \(\Phi(w)=\int_{0}^{w}\frac{\varrho}{f(\varrho)}\, d\varrho\) is well defined.

A simple calculation yields

$$ \mathscr{H}'(w)= \int_{\Omega}\frac{w}{f(w)}w_{t} \,dx=-p \mathscr{E}(w) $$

and

$$ \mathscr{H}''(w)=-p\mathscr{E}'(w)=p \int_{\Omega}\frac{w_{t}^{2}}{f(w)}\,dx. $$

Noticing f is nondecreasing, we have by Hölder’s inequality

$$\begin{aligned} \bigl(\mathscr{H}'(w) \bigr)^{2}&= \biggl( \int_{\Omega}\frac{w}{\sqrt {f(w)}}\frac{w_{t}}{\sqrt{f(w)}} \,dx \biggr)^{2}\le \biggl( \int_{\Omega}\frac{w^{2}}{f(w)}\,dx \biggr) \biggl( \int_{\Omega}\frac {w_{t}^{2}}{f(w)}\,dx \biggr) \\ &\le\frac{2}{p} \biggl( \int_{\Omega}\int_{0}^{w}\frac{\varrho}{f(\varrho)}\,d\varrho \,dx \biggr) \biggl(p \int_{\Omega}\frac{w_{t}^{2}}{f(w)}\,dx \biggr) \\ &=\frac{2}{p}\mathscr{H}(w)\mathscr{H}''(w), \end{aligned}$$

which implies

$$\frac{d^{2}}{dt^{2}}\mathscr{H}^{1-\frac{p}{2}}(t)\le 0. $$

It follows that there exists \(T<\infty\) such that if

$$ \mathscr{E}(w_{0})=\frac{1}{p} \int_{\Omega}|\nabla w_{0}|^{p} \,dx- \frac{c_{0}}{p} \int_{\Omega}w_{0}^{p}\,dx< 0, $$

then

$$ \limsup_{t\rightarrow T}\mathscr{H}(t)=\infty, $$

which implies

$$ \limsup_{t\rightarrow T}\bigl\Vert \Phi\bigl(w(\cdot,t)\bigr)\bigr\Vert _{L^{\infty}(\Omega)}=\infty. $$

Thus, we have

$$ \limsup_{t\rightarrow T}\bigl\Vert w(\cdot,t)\bigr\Vert _{L^{\infty}(\Omega)}=\infty. $$
(3.13)

Let \(\varphi>0\) in Ω be the first eigenfunction of (1.6). We have, for any \(k>0\),

$$ \mathscr{E}(k\psi)=\frac{1}{p} \int_{\Omega}\bigl\vert \nabla (k\psi)\bigr\vert ^{p} \,dx-\frac{c_{0}}{p} \int_{\Omega}(k\psi)^{p}\,dx=\frac{\lambda_{1}-c_{0}}{p} \int_{\Omega}(k\psi)^{p}\,dx< 0, $$

provided Ω large such that the first eigenvalue \(\lambda_{1}< c_{0}\). Choose k small enough such that \(w_{0}\ge k\psi\). Then w blows up in finite time by comparison.

Step 3. By the assumption, Ω contains a closed ball with radius R. Let \(\overline{B}_{R}(0)\subset\Omega\), without loss of generality. Suppose there is initial value \(u_{0}>0\) such that the solution u is global in time. Then we can show in a way similar to the proof of Lemma 3.2 that there is \(t_{0}>0\) such that \(u\ge M>f(v_{0})\) for \((x,t)\in B_{1}\times(t_{0},\infty)\), and hence \(v\ge f^{-1}(u)\ge f^{-1}(M)\ge v_{0}\) in \(B_{1}\times(t_{0},\infty)\). Therefore v is a non-global supersolution of (3.12) in \(B_{1}\times(t_{0},\infty)\) by Step 1. □

Theorem 3.2

Suppose that \(q>p+m-1\). Then the solutions of (1.1) blow up for large initial data and are global for small initial data.

Proof

Let ψ solve \(-\operatorname{div}(|\nabla\psi|^{p-2}\nabla \psi)=1\) in Ω with \(\psi|_{\partial\Omega}=0\). Since \(q>p+m-1\) implies \(\frac{pm-q}{p-1}< m-1\le r\), it follows that \(\frac{p-1}{q-m}<\frac{1}{m-r}\) when \(m>r\). Fix \(\delta>0\) such that

$$ \frac{p-1}{q-m}< \delta< \min \biggl\{ \frac{1}{m-r},1 \biggr\} \quad \mbox{for } m>r, $$
(3.14)

and set

$$w=\mu^{\delta}+\mu \psi(x),\quad (x,t)\in\Omega\times(0,\infty), $$

with \(\mu\in(0,1)\) small to be determined. A simple calculation yields

$$ w_{t}-w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2} \nabla w\bigr)-w^{q}-\gamma w^{r}|\nabla w|^{p}= \mu^{p-1}w^{m}-w^{q}-\gamma\mu^{p} w^{r}|\nabla\psi|^{p}. $$

Since \((q-m)\delta-p+1>0\) by (3.14), we have

$$ \frac{2w^{q}}{\mu^{p-1}w^{m}}=2\bigl(\mu^{\delta}+\mu\psi\bigr)^{q-m} \mu^{1-p} \le 2\mu^{(q-m)\delta-p+1}\bigl(1+\Vert \psi \Vert _{L^{\infty}(\Omega)}\bigr)\le1, $$

provided μ is small enough. If \(r< m\), by (3.14), we have

$$ \frac{2\gamma\mu^{p} w^{r}|\nabla \psi|^{p}}{\mu^{p-1}w^{m}}=2\gamma\mu\bigl(\mu^{\delta}+\mu\psi \bigr)^{r-m}|\nabla \psi|^{p} \le2\gamma\mu^{1-(m-r)\delta}| \nabla\psi|^{p}\le1 $$

for μ small. In the case \(r\ge m\),

$$ \frac{2\gamma\mu^{p} w^{r}|\nabla\psi|^{p}}{\mu^{p-1}w^{m}} \le 2\gamma\mu^{1+(r-m)\delta}\bigl(1+\Vert \psi \Vert _{L^{\infty}(\Omega )}\bigr)^{r-m}|\nabla \psi|^{p}\le1 $$
(3.15)

for μ small. In one word, there exists \(\mu>0\) small enough that w is a time-independent supersolution provided \(u_{0}< \mu^{\delta}+\mu\psi\).

To deal with large initial data, consider the following problem without gradient term:

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} u_{t}=u^{m}\operatorname{div}(|\nabla u|^{p-2}\nabla u)+\lambda u^{q} , &(x,t)\in\Omega\times(0,T), \\ u=0, &(x,t)\in \partial\Omega\times(0,T), \\ u(x,0)=u_{0}(x),& x\in\overline{\Omega}, \end{array}\displaystyle \right . $$

the solutions of which blow up in finite time for large initial data [6]. By the comparison principle, the solutions of (1.1) blow up as well for large initial data. □

Theorem 3.3

Suppose \(q=p+m-1\).

  1. (i)

    If \(r>m-1=\frac{pm-q}{p-1}\), there exist both global and non-global solutions.

  2. (ii)

    Suppose \(r=m-1=\frac{pm-q}{p-1}\). If \(\lambda_{1}< \lambda (\frac{\gamma+p-1}{p-1} )^{\frac{2(p-1)}{p}}\), all solutions of (1.1) blow up in finite time. If \(\lambda_{1}\ge \lambda (\frac{\gamma+p-1}{p-1} )^{\frac{2(p-1)}{p}}\), the solutions are global and bounded.

Proof

(i) It follows from Lemma 3.3 immediately that the solutions of (1.1) blow up in finite time for large initial data.

Next, we will show the solutions are global for small domain and small initial data. Let Ω be small with \(\lambda_{1}(\Omega)>\lambda\), and choose Ω̃ with \(\Omega\subset\subset\widetilde{\Omega}\) such that \(\tilde{\lambda}_{1}:=\lambda_{1}(\widetilde{\Omega})>\lambda\) (see [14]). Normalize φ̃, the eigenfunction corresponding to \(\tilde{\lambda_{1}}\), by \(\tilde{\varphi}>0\) in \(\Omega'\), \(\|\tilde{\varphi}\|_{L^{\infty}(\widetilde{\Omega})}=1\). Then \(\tilde{\varphi}\ge\rho\) in Ω for some \(\rho>0\).

Define \(w=a\tilde{\varphi}^{b}\) with \(a=(\frac{(1-b)(p-1)}{b\gamma})^{\frac{1}{r-m+1}}\), \(b\in (\lambda/\tilde{\lambda_{1}}, 1)\). Then

$$\begin{aligned}& w_{t}-w^{m}\operatorname{div}\bigl(|\nabla w|^{p-2} \nabla w\bigr)-\lambda w^{q}-\gamma u^{r}|\nabla w|^{p} \\& \quad =w^{m} \bigl[(ab)^{p-1}\tilde{\lambda}_{1} \tilde{\varphi }^{(b-1)(p-1)+(p-1)}-(b-1) (p-1) (ab)^{p-1}\tilde{\varphi }^{(b-1)(p-1)-1}|\nabla \tilde{\varphi}|^{p} \\& \qquad {}-\lambda a^{p-1}\tilde{\varphi}^{b(p-1)} -\gamma a^{p+r-m}b^{p}\tilde{\varphi}^{(b-1)p+b(r-m)}|\nabla\tilde{ \varphi }|^{p} \bigr] \\& \quad =w^{m} \bigl[a^{p-1}\tilde{\varphi}^{b(p-1)} \bigl(b^{p-1}\tilde {\lambda}_{1}-\lambda \bigr) \\& \qquad {}-(ab)^{p-1}\tilde{\varphi }^{(b-1)p-b}|\nabla \tilde{ \varphi}|^{p} \bigl((b-1) (p-1)+\gamma b a^{r-m+1}\tilde{ \varphi}^{b(r-m+1)} \bigr) \bigr] \\& \quad \ge0. \end{aligned}$$

If in addition \(u_{0}\) is small such that \(\|u_{0}\|_{L^{\infty}(\Omega)}< a\rho^{b}\), then \(u_{\epsilon}\le w\) in \(\Omega\times(0,\infty)\) by the comparison principle. Thus \(u\le w\) in \(\Omega\times(0,\infty)\).

(ii) By the re-scaling \(v(x,t)=u^{\alpha}(\frac{x}{\beta}, \frac{t}{\alpha})\) with \(\alpha=\frac{\gamma}{p-1}+1\) and \(\beta=\alpha^{\frac{p-1}{p}}\), we transform (1.1) into the Dirichlet problem

$$ \left \{ \textstyle\begin{array}{l@{\quad}l} v_{t}=v^{m_{1}}\operatorname{div}(|\nabla v|^{p-2}\nabla v)+\lambda v^{p+m_{1}-1}, &(x,t)\in U \times(0,\alpha T), \\ v=0, &(x,t)\in \partial U \times(0,\alpha T), \\ v(x,0)=u^{\alpha}_{0}(x),& x\in\overline{U}, \end{array}\displaystyle \right . $$
(3.16)

where \(m_{1}=\frac{\alpha+m-1-(\alpha-1)(p-1)}{\alpha}\), \(U =\{\beta x| x\in\Omega\}\). We know from Theorem 3.2 of [6] that all positive solutions of (3.16) are global and bounded if \(\lambda_{1}(U )\ge\lambda\), and so do the positive solutions of (1.1) whenever

$$ \lambda_{1}:=\lambda_{1}(\Omega)=\beta^{2} \lambda_{1}(U )\ge \lambda\beta^{2}=\lambda \biggl( \frac{\gamma+p-1}{p-1} \biggr)^{\frac {2(p-1)}{p}}. $$

Also by Theorem 3.2 of [6], similarly, if \(\lambda_{1}< \lambda(\frac{\gamma+p-1}{p-1})^{\frac{2(p-1)}{p}}\), all solutions of (1.1) blow up in finite time. □

References

  1. Friedman, A, McLeod, B: Blow-up of solutions of nonlinear parabolic equations. Arch. Ration. Mech. Anal. 79, 62-78 (1986)

    MathSciNet  MATH  Google Scholar 

  2. Wiegner, M: Blow-up for solutions of some degenerate parabolic equations. Differ. Integral Equ. 7, 1641-1647 (1994)

    MathSciNet  MATH  Google Scholar 

  3. Wiegner, M: A degenerate diffusion equation with a nonlinear source team. Nonlinear Anal. 28, 1977-1995 (1997)

    Article  MathSciNet  MATH  Google Scholar 

  4. Winkler, M: Blow-up in a degenerate parabolic equation. Indiana Univ. Math. J. 53, 1414-1442 (2004)

    Article  MathSciNet  MATH  Google Scholar 

  5. Stinner, C, Winkler, M: Boundedness vs. blow-up in a degenerate diffusion equation with gradient nonlinearity. Indiana Univ. Math. J. 56, 2233-2264 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  6. Jin, CH, Yin, JX: Critical exponent of a doubly degenerate parabolic equation in non-divergence form with nonlinear source. Chin. Ann. Math., Ser. A 30, 525-538 (2009)

    MathSciNet  MATH  Google Scholar 

  7. Zhou, WS, Wu, ZQ: Some results on a class of degenerate parabolic equations not in divergence form. Nonlinear Anal. 60, 863-886 (2005)

    Article  MathSciNet  MATH  Google Scholar 

  8. Yao, ZA, Zhou, WS: Some results on a degenerate and singular diffusion equation. Acta Math. Sin. Ser. B 27, 581-601 (2007)

    MathSciNet  MATH  Google Scholar 

  9. Zhou, WS: Some notes on a nonlinear degenerate parabolic equation. Nonlinear Anal. 71, 107-111 (2007)

    Article  MathSciNet  Google Scholar 

  10. Wang, J: Behaviors of solutions to a class of nonlinear degenerate parabolic equations not in divergence form. Appl. Math. Lett. 24, 191-195 (2011)

    Article  MathSciNet  MATH  Google Scholar 

  11. Winkler, M: Lager time behavior of solutions to degenerate parabolic equation with absorption. Nonlinear Differ. Equ. Appl. 8, 343-361 (2001)

    Article  MATH  Google Scholar 

  12. Quittner, P, Souplet, P: Superlinear Parabolic Problems. Birkhäuser, Basel (2007)

    MATH  Google Scholar 

  13. García-Melián, J, Sabina de Lis, J: Maximum and comparison principles for operators involving the p-Laplacian. J. Math. Anal. Appl. 218, 49-65 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  14. Li, YX, Xie, CH: Blow-up for p-Laplacian parabolic equations. Electron. J. Differ. Equ. 2003, 20 (2003)

    MathSciNet  MATH  Google Scholar 

Download references

Acknowledgements

The authors would like to express sincere gratitude to the referees for their valuable suggestions and comments on the original manuscript. The first author is supported by Hunan Provincial Natural Science Foundation of China (14JJ6044); the third author is supported by NSF of China (11501438) and the Natural Science Basic Research Plan in Shaanxi Province of China (2015JQ1014).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shuangshuang Zhou.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

The authors declare that the study was realized in collaboration with the same responsibility. All authors read and approved the final manuscript.

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

Zhou, S., Tang, X. & Yang, C. A doubly degenerate diffusion equation not in divergence form with gradient term. Bound Value Probl 2016, 126 (2016). https://doi.org/10.1186/s13661-016-0633-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-016-0633-4

MSC

Keywords