Skip to main content

Global existence and blow-up analysis for parabolic equations with nonlocal source and nonlinear boundary conditions

Abstract

We investigate the following nonlinear parabolic equations with nonlocal source and nonlinear boundary conditions:

$$ \textstyle\begin{cases} (g(u) )_{t} =\sum_{i,j=1}^{N} (a^{ij}(x)u_{x_{i}} ) _{x_{j}}+\gamma _{1}u^{m} (\int _{D} u^{l}{\,\mathrm{d}}x ) ^{p}-\gamma _{2}u^{r}& \mbox{in } D\times (0,t^{*}), \\ \sum_{i,j=1}^{N}a^{ij}(x)u_{x_{i}}\nu _{j}=h(u) & \mbox{on } \partial D\times (0,t^{*}), \\ u(x,0)=u_{0}(x)\geq 0 &\mbox{in } \overline{D}, \end{cases} $$

where p and \(\gamma _{1}\) are some nonnegative constants, m, l, \(\gamma _{2}\), and r are some positive constants, \(D\subset \mathbb{R}^{N}\) (\(N\geq 2\)) is a bounded convex region with smooth boundary ∂D. By making use of differential inequality technique and the embedding theorems in Sobolev spaces and constructing some auxiliary functions, we obtain a criterion to guarantee the global existence of the solution and a criterion to ensure that the solution blows up in finite time. Furthermore, an upper bound and a lower bound for the blow-up time are obtained. Finally, some examples are given to illustrate the results in this paper.

1 Introduction

There has been a vast amount of literature to discuss the global solutions, the blow-up solutions, and the bounds for the blow-up time for nonlinear parabolic equations. We refer the readers to [116]. The global and blow-up solutions for the parabolic problems with nonlinear boundary conditions have been studied in [17]. Recently, some authors have already considered the blow-up phenomena of the parabolic equations with nonlocal source. We recommend the literature [810] and the references therein. In the present work, we study the following nonlinear parabolic equations with nonlocal source and nonlinear boundary conditions:

$$ \textstyle\begin{cases} (g(u) )_{t} =\sum_{i,j=1}^{N} (a^{ij}(x)u_{{x}_{i}} ) _{{x}_{j}}+\gamma _{1}u^{m} (\int _{D} u^{l}{\,\mathrm{d}}x ) ^{p}-\gamma _{2}u^{r} & \mbox{in } D\times (0,t^{*}), \\ \sum_{i,j=1}^{N}a^{ij}(x)u_{{x}_{i}}\nu _{j}=h(u) & \mbox{on } \partial D\times (0,t^{*}), \\ u(x,0)=u_{0}(x)\geq 0 & \mbox{in } \overline{D}, \end{cases} $$
(1.1)

where p and \(\gamma _{1}\) are some nonnegative constants, m, l, \(\gamma _{2}\), and r are some positive constants, \(D\subset \mathbb{R}^{N}\) (\(N\geq 2\)) is a bounded convex region, the boundary ∂D is smooth, \(t^{*}\) is the blow-up time when blow-up occurs, or \(t^{*}=+\infty \), and \(\nu =(\nu _{1},\nu _{2}, \ldots ,\nu _{N})\) is the unit outward normal vector on ∂D. Moreover, \((a^{ij}(x))_{N\times N}\) is a differentiable positive definite matrix; that is, for all \(z=(z_{1},z_{2},\ldots ,z_{N}) \in \mathbb{R}^{N}\), there exists a constant \(\theta >0\) such that

$$ \sum_{i,j=1}^{N}a^{ij}z_{i}z_{j} \geq \theta \vert z \vert ^{2}. $$
(1.2)

Set \(\mathbb{R}_{+}=(0,+\infty )\). We assume that g is a \(C^{2}(\overline{ \mathbb{R}_{+}})\) function which satisfies \(g'(s)>0\) for all \(s \in \overline{\mathbb{R}_{+}}\), h is a nonnegative \(C^{1}(\overline{ \mathbb{R}_{+}})\) function, and \(u_{0}(x)\) is a nonnegative \(C^{1}(\overline{D})\) function which satisfies the compatibility condition. Our study is motivated by the following three papers. Li et al. [1], Baghaei et al. [11], and Zhang et al. [2] studied the following parabolic equations with local source and nonlinear boundary conditions:

$$ \textstyle\begin{cases} u_{t} =\sum_{i,j=1}^{N} (a^{ij}(x)u_{{x}_{i}} )_{{x}_{j}}+k(t)f(u) & \mbox{in } D\times (0,t^{*}), \\ \sum_{i,j=1}^{N}a^{ij}(x)u_{{x}_{i}}\nu _{j}=h(u) &\mbox{on } \partial D\times (0,t^{*}), \\ u(x,0)=u_{0}(x)\geq 0 &\mbox{in } \overline{D}, \end{cases} $$
(1.3)

where D is a bounded convex region in \(\mathbb{R}^{N}\) (\(N\geq 2\)), and the boundary ∂D is smooth. When the function \(k(t)\equiv -1\), Li et al. [1] obtained the conditions to ensure that the global solution exists and the solution blows up in finite time \(t^{*}\). Moreover, they derived upper bounds of the blow-up time in \(D\subset \mathbb{R}^{N}\) (\(N\geq 2\)). Using the three-dimensional Sobolev type inequality that Payne and Schaefer had proven in [7], they got lower bounds in \(D\subset \mathbb{R}^{3}\). When \(k(t)\equiv -1\), in [11], Baghaei et al. used the embedding theorems in Sobolev spaces to get a lower bound of the blow-up time in \(D\subset \mathbb{R}^{N}\) (\(N\geq 3\)). Recently, Zhang et al. [2] considered problem (1.3) when \(k(t)>0\). Their main contribution was to extend the three-dimensional Sobolev type inequality obtained by Payne and Schaefer in [7] to the higher dimensional one. They used this Sobolev inequality in multidimensional space to derive a lower bound of the blow-up time in \(D\subset \mathbb{R}^{N}\) (\(N\geq 3\)). In addition, they obtained the global existence of the solution and an upper bound of the blow-up time when blow-up occurs in \(D\subset \mathbb{R}^{N}\) (\(N\geq 2\)).

In the present paper, we study problem (1.1). In the process of getting lower bounds of the blow-up time, the key of our work is to deal with the nonlocal source. We discuss the blow-up problems of (1.1) by constructing some suitable auxiliary functions and making use of the embedding theorems in Sobolev spaces and differential inequality technique. The present work is organized as follows. In Sect. 2, we establish some conditions on the data g and h to obtain that the solution \(u(x,t)\) exists globally in \(D\subset \mathbb{R}^{N}\) (\(N \geq 2\)). In Sect. 3, we develop conditions on the data of (1.1) to guarantee the blow-up of the solution and derive an upper bound of blow-up time in \(D\subset \mathbb{R}^{N}\) (\(N\geq 2\)). When blow-up does occur, we derive a lower bound of \(t^{*}\) for \(D\subset \mathbb{R} ^{N}\) (\(N\geq 3\)) in Sect. 4 and a lower bound of \(t^{*}\) for \(D\subset \mathbb{R}^{2}\) in Sect. 5. In Sect. 6, some examples are presented to illustrate the results in this paper.

2 The global solution

In this section, we establish some conditions on the data of (1.1) to ensure that the global solution exists. We define the auxiliary functions

$$ \varPhi (t)= \int _{D} G(u)\,\mathrm{d}x,\qquad G(s)=l \int _{0}^{s}y^{l-1}g'(y) \,\mathrm{d}y, \quad s\geq 0. $$
(2.1)

Now we state our main results as follows.

Theorem 2.1

Letube a nonnegative classical solution of problem (1.1). Suppose that the functionsgandhsatisfy

$$ g'(s)\leq \zeta ,\qquad h(s)\leq \gamma s^{q},\quad s\geq 0, $$
(2.2)

whereζ, γ, andqare some positive constants. Assume

$$ l>2,\qquad q>\max \biggl\{ 1, \frac{m+1+lp}{2} \biggr\} ,\qquad r>2q-1. $$
(2.3)

Then\(u(x,t)\)exists for all\(t>0\)in the measure\(\varPhi (t)\).

Proof

It follows from (2.2) that

$$ \varPhi (t)=l \int _{D} \int _{0}^{u}y^{l-1}g'(y) \,\mathrm{d}y\,\mathrm{d}x \leq \zeta \int _{D} u^{l}{\,\mathrm{d}}x; $$

that is,

$$ \int _{D} u^{l}{\,\mathrm{d}}x\geq \frac{1}{\zeta }\varPhi (t). $$
(2.4)

Using (1.2), (2.2), and the divergence theorem, we have

$$\begin{aligned} \varPhi '(t) =&l \int _{D} u^{l-1}g'(u)u_{t} {\,\mathrm{d}}x=l \int _{D} u ^{l-1} \Biggl(\sum _{i,j=1}^{N} \bigl(a^{ij}(x)u_{x_{i}} \bigr)_{x_{j}}+ \gamma _{1} u^{m} \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr)^{p}-\gamma _{2}u^{r} \Biggr) \,\mathrm{d}x \\ =&l \int _{D} \sum_{i,j=1}^{N} \bigl(u^{l-1}a^{ij}(x)u_{x_{i}} \bigr) _{x_{j}}{\,\mathrm{d}}x-l(l-1) \int _{D} u^{l-2}\sum _{i,j=1}^{N}a^{ij}(x)u _{x_{i}}u_{x_{j}}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D} u^{l} {\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x \\ =&l \int _{\partial D}u^{l-1}\sum _{i,j=1}^{N}a^{ij}(x)u_{x_{i}} \nu _{j} {\,\mathrm{d}}S-l(l-1) \int _{D} u^{l-2}\sum _{i,j=1}^{N}a^{ij}(x)u_{x_{i}}u _{x_{j}}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D} u^{l} {\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x \\ \le &l \int _{\partial D} u^{l-1}h(u)\,\mathrm{d}S-\theta l(l-1) \int _{D}u ^{l-2} \vert \nabla u \vert ^{2}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D}u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u ^{l+r-1}{\,\mathrm{d}}x \\ \le &l\gamma \int _{\partial D} u^{l+q-1}{\,\mathrm{d}}S- \frac{4\theta (l-1)}{l} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2} {\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D}u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x. \end{aligned}$$
(2.5)

By the divergence theorem, we get the inequality (see [5])

$$ \int _{\partial D} u^{l+q-1} {\,\mathrm{d}}S \leq \frac{N}{\rho _{0}} \int _{D}u^{l+q-1} {\,\mathrm{d}}x+ \frac{(l+q-1)d}{\rho _{0}} \int _{D}u^{l+q-2} \vert \nabla u \vert { \,\mathrm{d}}x, $$
(2.6)

where

$$ \rho _{0}=\min_{x\in {\partial D}}(x\cdot \nu ),\qquad d=\max _{x\in \overline{D}} \vert x \vert . $$

Applying the Hölder inequality and the Young inequality to the right-hand side of (2.6) yields

$$\begin{aligned} \frac{N}{\rho _{0}} \int _{D}u^{l+q-1}{\,\mathrm{d}}x \le& \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr)^{2}{ \,\mathrm{d}}x \biggr)^{\frac{1}{2}} \biggl( \int _{D} \biggl(\frac{N}{\rho _{0}} u^{\frac{l+2q-2}{2}} \biggr) ^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}} \\ \leq& \frac{1}{2} \int _{D} u ^{l}{\,\mathrm{d}}x+ \frac{N^{2}}{2\rho _{0}^{2}} \int _{D} u^{l+2q-2} {\,\mathrm{d}}x \end{aligned}$$
(2.7)

and

$$\begin{aligned}& \frac{(l+q-1)d}{\rho _{0}} \int _{D} u^{l+q-2} \vert \nabla u \vert { \,\mathrm{d}}x \\& \quad \le \biggl( \int _{D}\frac{1}{\varepsilon } \biggl(\frac{(l+q-1)d}{\rho _{0}} u^{\frac{l+2q-2}{2}} \biggr)^{2}{\,\mathrm{d}}x \biggr)^{ \frac{1}{2}} \biggl( \int _{D} \varepsilon \bigl(u^{\frac{l-2}{2}} \vert \nabla u \vert \bigr)^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}} \\& \quad \le \frac{(l+q-1)^{2}d^{2}}{2\varepsilon \rho _{0}^{2}} \int _{D} u ^{l+2q-2} {\,\mathrm{d}}x+ \frac{\varepsilon }{2} \int _{D}u^{l-2} \vert \nabla u \vert ^{2}{\,\mathrm{d}}x \\& \quad = \frac{(l+q-1)^{2}d^{2}}{2\varepsilon \rho _{0}^{2}} \int _{D} u^{l+2q-2} {\,\mathrm{d}}x+ \frac{2\varepsilon }{l^{2}} \int _{D} \bigl\vert \nabla u^{ \frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \end{aligned}$$
(2.8)

with

$$ \varepsilon =\frac{2\theta (l-1)}{\gamma }>0. $$
(2.9)

Hence, substituting (2.7) and (2.8) into (2.6), we get

$$\begin{aligned} \int _{\partial D} u^{l+q-1} {\,\mathrm{d}}S \le& \frac{1}{2} \int _{D} u ^{l} {\,\mathrm{d}}x+ \biggl( \frac{N^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d ^{2}}{2\varepsilon \rho _{0}^{2}} \biggr) \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \\ &{}+ \frac{2\varepsilon }{l^{2}} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x. \end{aligned}$$
(2.10)

Inserting (2.10) into (2.5) and using (2.9), we obtain

$$\begin{aligned} \varPhi '(t) \le &\frac{1}{2}l\gamma \int _{D} u^{l}{\,\mathrm{d}}x+l \gamma \biggl( \frac{N^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d^{2}}{2 \varepsilon \rho _{0}^{2}} \biggr) \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D}u^{l+m-1} {\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x. \end{aligned}$$
(2.11)

Since (2.3) implies the fact that \(0<\frac{l}{l+2q-2}<1\), \(0< \frac{l+m-1}{l+2q-2}<1\), and \(0<\frac{l+2q-2}{l+s-1}<1\), it follows from the Hölder inequality that

$$\begin{aligned}& \int _{D} u^{l}{\,\mathrm{d}}x\leq \vert D \vert ^{\frac{2(q-1)}{l+2q-2}} \biggl( \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \biggr)^{\frac{l}{l+2q-2}}, \end{aligned}$$
(2.12)
$$\begin{aligned}& \int _{D} u^{l+m-1}{\,\mathrm{d}}x\leq \vert D \vert ^{\frac{2q-m-1}{l+2q-2}} \biggl( \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \biggr)^{\frac{l+m-1}{l+2q-2}}, \end{aligned}$$
(2.13)

and

$$ \int _{D} u^{l+2q-2}{\,\mathrm{d}}x\leq \vert D \vert ^{\frac{r-2q+1}{l+r-1}} \biggl( \int _{D}u^{l+r-1}{\,\mathrm{d}}x \biggr)^{\frac{l+2q-2}{l+r-1}}, $$
(2.14)

where \(|D|\) is the volume of D. By virtue of (2.14), we have

$$ \int _{D}u^{l+r-1}{\,\mathrm{d}}x\geq \vert D \vert ^{\frac{2q-r-1}{l+2q-2}} \biggl( \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \biggr)^{\frac{l+r-1}{l+2q-2}}. $$
(2.15)

We insert (2.12)–(2.13) and (2.15) into (2.11) to derive

$$\begin{aligned} \varPhi '(t) \le & \frac{1}{2}l\gamma \vert D \vert ^{\frac{2(q-1)}{l+2q-2}} \biggl( \int _{D} u ^{l+2q-2}{\,\mathrm{d}}x \biggr)^{\frac{l}{l+2q-2}} +l\gamma \biggl(\frac{N ^{2}}{2\rho _{0}^{2}}+ \frac{(l+q-1)^{2}d^{2}}{2\varepsilon \rho _{0} ^{2}} \biggr) \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \vert D \vert ^{\frac{2q-m-1+2p(q-1)}{l+2q-2}} \biggl( \int _{D}u^{l+2q-2} {\,\mathrm{d}}x \biggr)^{\frac{l+m-1+lp}{l+2q-2}} \\ &{}-l\gamma _{2} \vert D \vert ^{ \frac{2q-r-1}{l+2q-2}} \biggl( \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \biggr) ^{\frac{l+r-1}{l+2q-2}} \\ =& \int _{\varOmega }u^{l+2q-2}{\,\mathrm{d}}x \biggl(l\gamma \biggl(\frac{N ^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d^{2}}{2\varepsilon \rho _{0} ^{2}} \biggr)+\frac{l\gamma }{2} \vert D \vert ^{\frac{2(q-1)}{l+2q-2}} \biggl( \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \biggr)^{\frac{2(1-q)}{l+2q-2}} \\ &{}+l\gamma _{1} \vert D \vert ^{\frac{2q-m-1+2p(q-1)}{l+2q-2}} \biggl( \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \biggr)^{\frac{m-2q+1+lp}{l+2q-2}} \\ &{}-l\gamma _{2} \vert D \vert ^{\frac{2q-r-1}{l+2q-2}} \biggl( \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \biggr) ^{\frac{r-2q+1}{l+2q-2}} \biggr). \end{aligned}$$
(2.16)

It follows from (2.12) and (2.4) that

$$\begin{aligned} \int _{D} u^{l+2q-2}{\,\mathrm{d}}x\geq \vert D \vert ^{\frac{2(1-q)}{l}} \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr)^{\frac{l+2q-2}{l}}\geq \zeta ^{\frac{2-2q-l}{l}} \vert D \vert ^{\frac{2(1-q)}{l}} \varPhi ^{\frac{l+2q-2}{l}}(t). \end{aligned}$$
(2.17)

By (2.3), we know \(\frac{2(1-q)}{l+2q-2}<0\), \(\frac{m-2q+1+lp}{l+2q-2}<0\), and \(\frac{r-2q+1}{l+2q-2}>0\). Substituting (2.17) into (2.16), we deduce

$$\begin{aligned} \varPhi '(t) \le & \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \biggl[l\gamma \biggl(\frac{N ^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d^{2}}{2\varepsilon \rho _{0} ^{2}} \biggr)+ \frac{l\gamma }{2} \zeta ^{\frac{2(q-1)}{l}} \vert D \vert ^{ \frac{2(q-1)}{l}}\varPhi ^{\frac{2(1-q)}{l}}(t) \\ &{}+l\gamma _{1}\zeta ^{\frac{2q-m-1}{l}-p} \vert D \vert ^{\frac{2q-m-1}{l}} \varPhi ^{p+\frac{m-2q+1}{l}}(t) -l\gamma _{2}\zeta ^{\frac{2q-r-1}{l}} \vert D \vert ^{ \frac{2q-r-1}{l}}\varPhi ^{\frac{{r-2q+1}^{\frac{}{}}}{l}}(t) \biggr] \\ =& \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \bigl(C_{1}+C_{2} \varPhi ^{\frac{2(1-q)}{l}}(t)+C_{3} \varPhi ^{p+\frac{m-2q+1}{l}}(t)-C_{4} \varPhi ^{\frac{r-2q+1}{l}}(t) \bigr), \end{aligned}$$
(2.18)

where \(C_{1}=l\gamma (\frac{N^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d ^{2}}{2\varepsilon \rho _{0}^{2}} )\), \(C_{2}=\frac{l\gamma }{2} \zeta ^{\frac{2(q-1)}{l}}|D|^{\frac{2(q-1)}{l}}\), \(C_{3}=l\gamma _{1} \zeta ^{\frac{2q-m-1}{l}-p}|D|^{\frac{2q-m-1}{l}}\), and \(C_{4}=l \gamma _{2}\zeta ^{\frac{2q-r-1}{l}}|D|^{\frac{2q-r-1}{l}}\).

Thus, (2.18) implies that \(u(t)\) cannot blow up for all time \(t>0\) in the measure \(\varPhi (t)\). In fact, if \(u(x,t)\) blows up at finite time \(t^{*}\) in the measure \(\varPhi (t)\), we have

$$ \lim_{t\rightarrow {t^{*-}}}\varPhi (t)=+\infty . $$

It follows from (2.3) that \(\frac{2-2q}{l}<0\), \({p+\frac{m-2q+1}{l}}<0\), and \(\frac{r-2q+1}{l}>0\). From (2.18), we obtain that \(\varPhi '(t)<0\) in some interval \((t_{0}, t^{*})\). So, for any \(t\in [t_{0}, t^{*})\), we get \(\varPhi (t)\leq \varPhi (t_{0})\). We take the limit as \(t\rightarrow {t^{*-}}\) to get

$$ +\infty =\lim_{t\rightarrow {t^{*-}}}\varPhi (t)\leq \varPhi (t_{0}). $$

This is a contradiction. □

3 Upper bounds of the blow-up time \(t^{*}\)

In this section, we set up some conditions on g and k to guarantee that the solution of (1.1) blows up in finite time. An upper bound of the blow-up time \(t^{*}\) is obtained in \(D\subset \mathbb{R}^{N}\) (\(N \geq 2\)). The auxiliary functions of this section are defined as follows:

$$ \varPsi (t)= \int _{D}g(u)\,\mathrm{d}x,\quad t\geq 0. $$
(3.1)

Our main result is the following Theorem 3.1.

Theorem 3.1

Let\(u(x,t)\)be a nonnegative classical solution of problem (1.1). Suppose that the functiongsatisfies

$$ 0\leq g(s)\leq s,\quad s\geq 0. $$
(3.2)

Assume

$$ 1< r< \min \{m, l\} \quad \textit{and}\quad \varPsi (0)> \biggl( \frac{M_{2}}{M_{1}} \biggr)^{\frac{1}{m+lp}} $$
(3.3)

with

$$ \begin{aligned} &M_{1}=\frac{\gamma _{1}}{2} \vert D \vert ^{1+p-m-lp},\qquad M_{2}=\frac{\gamma _{2}(m+lp-r)}{m+lp}\delta ^{\frac{r}{r-m-lp}}, \\ &\delta =\frac{\gamma _{1}(m+lp)}{2r\gamma _{2}} \vert D \vert ^{1+p-\frac{m+lp}{r}}. \end{aligned} $$
(3.4)

Then the solution\(u(x,t)\)must blow up in the measure\(\varPsi (t)\)in finite time\(t^{*}\)and

$$ t^{*}\leq\int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M_{2}}. $$

Proof

Making use of the divergence theorem, we get

$$\begin{aligned} \varPsi '(t) =& \int _{D}g'(u)u_{t}{ \,\mathrm{d}}x = \int _{D} \Biggl(\sum_{i,j=1} ^{N} \bigl(a^{ij}(x)u_{{x}_{i}} \bigr)_{{x}_{j}}+\gamma _{1}u^{m} \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr)^{p}-\gamma _{2}u^{r} \Biggr) \,\mathrm{d}x \\ =& \int _{\partial D}\sum_{i,j=1}^{N}a^{ij}(x)u_{{x}_{i}} \nu _{j}{\,\mathrm{d}}S +\gamma _{1} \int _{D}u^{m}{\,\mathrm{d}}x \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr) ^{p}-\gamma _{2} \int _{D}u^{r}{\,\mathrm{d}}x \\ =& \int _{\partial D}h(u)\,\mathrm{d}S +\gamma _{1} \int _{D}u^{m}{\,\mathrm{d}}x \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr)^{p}-\gamma _{2} \int _{D}u ^{r}{\,\mathrm{d}}x \\ \geq &\gamma _{1} \int _{D}u^{m}{\,\mathrm{d}}x \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr) ^{p}-\gamma _{2} \int _{D}u^{r}{\,\mathrm{d}}x. \end{aligned}$$
(3.5)

By the Hölder inequality and (3.3), we can compute

$$\begin{aligned}& \int _{D}u^{r}{\,\mathrm{d}}x\leq \vert D \vert ^{1-\frac{r}{m}} \biggl( \int _{D}u ^{m}{\,\mathrm{d}}x \biggr)^{\frac{r}{m}},\qquad \int _{D}u^{r}{\,\mathrm{d}}x\leq \vert D \vert ^{1-\frac{r}{l}} \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{\frac{r}{l}}, \\& \int _{D}u\,\mathrm{d}x\leq \vert D \vert ^{1-\frac{1}{r}} \biggl( \int _{D}u^{r}{\,\mathrm{d}}x \biggr) ^{\frac{1}{r}}; \end{aligned}$$

that is,

$$ \begin{aligned} &\int _{D}u^{m}{\,\mathrm{d}}x\geq \vert D \vert ^{1-\frac{m}{r}} \biggl( \int _{D}u ^{r}{\,\mathrm{d}}x \biggr)^{\frac{m}{r}},\qquad \int _{D}u^{l}{\,\mathrm{d}}x\geq \vert D \vert ^{1-\frac{l}{r}} \biggl( \int _{D}u ^{r}{\,\mathrm{d}}x \biggr)^{\frac{l}{r}}, \\ &\int _{D}u^{r}{\,\mathrm{d}}x\geq \vert D \vert ^{1-r} \biggl( \int _{D}u\,\mathrm{d}x \biggr) ^{r}. \end{aligned} $$
(3.6)

Substituting (3.6) into (3.5) and applying the Young inequality and (3.2)–(3.4), we deduce

$$\begin{aligned} \varPsi '(t) \geq &\gamma _{1} \vert D \vert ^{1+p-\frac{m+lp}{r}} \biggl( \int _{D}u ^{r}{\,\mathrm{d}}x \biggr)^{\frac{m+lp}{r}}-\gamma _{2} \int _{D}u^{r} {\,\mathrm{d}}x \\ =&\gamma _{1} \vert D \vert ^{1+p-\frac{m+lp}{r}} \biggl( \int _{D}u^{r}{\,\mathrm{d}}x \biggr) ^{\frac{m+lp}{r}}-\gamma _{2} \biggl(\delta \biggl( \int _{D}u^{r}{\,\mathrm{d}}x \biggr) ^{\frac{m+lp}{r}} \biggr)^{\frac{r}{m+lp}} \bigl(\delta ^{ \frac{r}{r-m-lp}} \bigr)^{\frac{m+lp-r}{m+lp}} \\ \geq & \biggl(\gamma _{1} \vert D \vert ^{1+p-\frac{m+lp}{r}}- \frac{\delta r\gamma _{2}}{m+lp} \biggr) \biggl( \int _{D}u^{r}{\,\mathrm{d}}x \biggr)^{ \frac{m+lp}{r}}-\frac{\gamma _{2}(m+lp-r)}{m+lp}\delta ^{ \frac{r}{r-m-lp}} \\ \geq & \biggl(\gamma _{1} \vert D \vert ^{1+p-\frac{m+lp}{r}}- \frac{\delta r\gamma _{2}}{m+lp} \biggr) \vert D \vert ^{\frac{(1-r)(m+lp)}{r}} \biggl( \int _{D}u \,\mathrm{d}x \biggr)^{m+lp}- \frac{\gamma _{2}(m+lp-r)}{m+lp} \delta ^{\frac{r}{r-m-lp}} \\ =&M_{1} \biggl( \int _{D}u\,\mathrm{d}x \biggr)^{m+lp}-M_{2} \geq M _{1} \biggl( \int _{D}g(u)\,\mathrm{d}x \biggr)^{m+lp}-M_{2} \\ =&M_{1}\varPsi ^{m+lp}(t)-M_{2}. \end{aligned}$$
(3.7)

We know that (3.3) implies

$$ M_{1}\varPsi ^{m+lp}(t)-M_{2}>0,\quad t\geq 0. $$
(3.8)

In fact, if (3.8) does not hold, we let

$$ t'=\min\bigl\{ M_{1}\varPsi ^{m+lp}(t)-M_{2}\leq 0\bigr\} . $$
(3.9)

From (3.9), we have

$$ M_{1}\varPsi ^{m+lp}(t)-M_{2}>0,\quad 0< t< t'. $$
(3.10)

It follows from (3.7) and (3.10) that \(\varPsi '(t)>0\), \(0< t< t'\), and \(\varPsi (t')>\varPsi (0)\). By (3.3), we derive

$$ M_{1}\varPsi ^{m+lp}\bigl(t' \bigr)-M_{2}>0, $$

which contradicts with (3.9).

Integrating (3.7) over \([0,t]\), we obtain

$$ t\leq \int _{\varPsi (0)}^{\varPsi (t)}\frac{d\tau }{M_{1}\tau ^{m+lp}-M_{2}}, $$
(3.11)

which implies that the solution u must blow up at some finite time \(t^{*}\) in the measure \(\varPsi (t)\). In fact, if u does not blow up at \(t^{*}\) in the measure \(\varPsi (t)\), we get

$$ \varPsi (t)< +\infty ,\quad t\geq 0. $$

Therefore,

$$ t\leq \int _{\varPsi (0)}^{\varPsi (t)}\frac{d\tau }{M_{1}\tau ^{m+lp}-M_{2}}< \int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M_{2}},\quad t \geq 0. $$

Letting \(t\rightarrow +\infty \), we have

$$ \int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M_{2}}=+ \infty , $$

which contradicts with

$$ \int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M_{2}}< + \infty . $$

Hence, u blows up at \(t^{*}\) in the measure \(\varPsi (t)\). We pass to the limit as \(t\rightarrow t^{*-}\) in (3.11) to derive

$$ t^{*}\leq \int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M _{2}}. $$

 □

4 Lower bounds of the blow-up time \(t^{*}\) in \(D\subset \mathbb{R}^{N}\) (\(N\geq 3\))

In this section, we impose restriction \(D\subset \mathbb{R}^{N}\)\((N\geq 3)\). Assume that the functions h and g satisfy

$$ h(s)\leq \sigma s^{q},\qquad g'(s)\geq \xi , \quad s\geq 0, $$
(4.1)

where σ, q, and ξ are some positive constants and

$$ q>1. $$
(4.2)

It follows from [17, Corollary 9.14] that

$$ W^{1, 2}\hookrightarrow L^{\frac{2N}{N-2}}, \quad \mbox{where } N\geq 3, $$

from which we have the following Sobolev inequality:

$$ \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr)^{\frac{2N}{N-2}}{ \,\mathrm{d}}x \biggr) ^{\frac{N-2}{2N}} \leq C \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}}, $$
(4.3)

where \(C=C(N,D)\) is a Sobolev embedding constant depending on N (\(N\geq 3\)) and D. In order to derive a lower bound of the blow-up time \(t^{*}\), we use inequality (4.3) and

$$ (a+b)^{\mu }\leq a^{\mu }+b^{\mu },\quad a\geq 0, b\geq 0, 0< \mu \leq 1. $$
(4.4)

We define auxiliary functions as follows:

$$ A(t)= \int _{D} B(u)\,\mathrm{d}x,\quad t\geq 0,\qquad B(s)=l \int _{0}^{s}y^{l-1}g'(y) \,\mathrm{d}y, \quad s\geq 0. $$
(4.5)

Now we present our main results in Theorem 4.1.

Theorem 4.1

Let\(u({x},t)\)be a nonnegative classical solution of (1.1) in\(D\subset \mathbb{R}^{N}\) (\(N\geq 3\)). Suppose that (4.1)(4.2) hold and

$$ r>1, \qquad m>1,\qquad l>\max \biggl\{ 2, \frac{N(m-1)}{2}, N(q-1) \biggr\} . $$
(4.6)

If\(u({x},t)\)becomes unbounded at some finite time\(t^{*}\)in the measure\(A(t)\), then we conclude\(t^{*}\)is bounded from below by

$$ t^{*}\geq \int _{A(0)}^{+\infty }\frac{\mathrm{d}\tau }{J_{1}\tau +J _{2}\tau ^{1+\frac{2(q-1)}{l}}+J_{3}\tau ^{1+\frac{2(q-1)}{l-N(q-1)}}+J _{4}\tau ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} +J_{5} \tau ^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}}}, $$

where

$$\begin{aligned}& J_{1}=\frac{l\sigma }{2\xi },\qquad J_{2}= \frac{l\sigma }{2\rho _{0}^{2}\xi ^{1+\frac{2(q-1)}{l}}}C^{ \frac{2N(q-1)}{l}} \biggl(N^{2}+ \frac{\sigma (l+q-1)^{2}d^{2}}{\theta (l-1)} \biggr), \end{aligned}$$
(4.7)
$$\begin{aligned}& J_{3}=\frac{l-N(q-1)}{\xi ^{1+\frac{2(q-1)}{l-N(q-1)}}} \biggl( \frac{ \theta (l-1)}{lN(q-1)C^{2}} \biggr)^{\frac{N(1-q)}{l-N(q-1)}} \biggl[\frac{ \sigma }{2\rho _{0}^{2}} \biggl(N^{2}+\frac{\sigma (l+q-1)^{2}d^{2}}{ \theta (l-1)} \biggr) \biggr]^{\frac{l}{l-N(q-1)}}, \end{aligned}$$
(4.8)
$$\begin{aligned}& J_{4}=\frac{ \vert D \vert (r-1)}{\xi ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}}}\gamma _{2} ^{\frac{l}{1-r}} \biggl(\frac{l\gamma _{1}}{l+r-1}C^{\frac{N(m-1)}{l}} \biggr) ^{1+\frac{l}{r-1}}, \end{aligned}$$
(4.9)
$$\begin{aligned}& J_{5}=\frac{2l-N(m-1)}{\xi ^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}}} \biggl( \frac{ \gamma _{1}}{2} \biggr)^{\frac{2l}{2l-N(m-1)}} \biggl(\frac{\theta (l-1)}{lN(m-1)C ^{2}} \biggr)^{\frac{N(1-m)}{2l-N(m-1)}}, \end{aligned}$$
(4.10)
$$\begin{aligned}& \rho _{0}=\min_{{x}\in {\partial D}}(x\cdot \nu ), \qquad d=\max_{x\in \overline{D}} \vert x \vert , \end{aligned}$$
(4.11)

andθis defined in (1.2).

Proof

We use (1.2), (4.1), and the divergence theorem to obtain

$$\begin{aligned} A'(t) =&l \int _{D} u^{l-1}g'(u)u_{t} {\,\mathrm{d}}x =l \int _{D} u^{l-1} \Biggl(\sum _{i,j=1}^{N} \bigl(a^{ij}(x)u_{{x}_{i}} \bigr)_{{x}_{j}}+ \gamma _{1}u^{m} \biggl( \int _{D} u^{l}{\,\mathrm{d}}x \biggr)^{p}-\gamma _{2}u^{r} \Biggr) \,\mathrm{d}x \\ =&l \int _{D}\sum_{i,j=1}^{N} \bigl(u^{l-1}a^{ij}({x})u_{{ x}_{i}} \bigr) _{{x}_{j}}{\,\mathrm{d}}x-l(l-1) \int _{D} u^{l-2}\sum _{i,j=1}^{N}a^{ij}(x)u _{{x}_{i}}u_{{x}_{j}}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x \\ \le &l \int _{\partial D}u^{l-1}\sum _{i,j=1}^{N}a^{ij}(x)u_{{ x}_{i}} \nu _{j}{\,\mathrm{d}}S-l\theta (l-1) \int _{D}u^{l-2} \vert \nabla u \vert ^{2}{\,\mathrm{d}}x \\ &{}+l \gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr) ^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x \\ \le &l \int _{\partial D} u^{l-1}h(u)\,\mathrm{d}S- \frac{4\theta (l-1)}{l} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2} {\,\mathrm{d}}x \\ &{}+l \gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr) ^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x \\ \le &l\sigma \int _{\partial D} u^{l+q-1}{\,\mathrm{d}}S- \frac{4\theta (l-1)}{l} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2} {\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr) ^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x. \end{aligned}$$
(4.12)

Repeating the calculation process of (2.6)–(2.8), we have

$$\begin{aligned} \int _{\partial D} u^{l+q-1} {\,\mathrm{d}}S \le& \frac{1}{2} \int _{D} u ^{l} {\,\mathrm{d}}x+ \biggl( \frac{N^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d ^{2}}{2\varepsilon _{1}\rho _{0}^{2}} \biggr) \int _{D} u^{l+2q-2}{\,\mathrm{d}}x \\ &{}+ \frac{2\varepsilon _{1}}{l^{2}} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x, \end{aligned}$$
(4.13)

where \(\rho _{0}\) and d are given in (4.11) and \(\varepsilon _{1}=\frac{ \theta (l-1)}{\sigma }\). Substituting (4.13) into (4.12), we derive

$$\begin{aligned} A'(t) \leq & \frac{l\sigma }{2} \int _{D}u^{l}{\,\mathrm{d}}x+l\sigma \biggl( \frac{N^{2}}{2\rho _{0}^{2}} +\frac{(l+q-1)^{2}d^{2}}{2 \varepsilon _{1}\rho _{0}^{2}} \biggr) \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \\ &{}+ \biggl( \frac{2\varepsilon _{1}\sigma }{l}-\frac{4\theta (l-1)}{l} \biggr) \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x. \end{aligned}$$
(4.14)

First we take into account the second term of (4.14). Assumptions (4.2) and (4.6) imply \(0<\frac{(N-2)(q-1)}{l}<1\) and \(0<\frac{N(q-1)}{l}<1\). Now we apply the Hölder inequality, the Young inequality, and inequalities (4.3)–(4.4) to deduce

$$\begin{aligned}& \int _{D}u^{l+2q-2}{\,\mathrm{d}}x \\& \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{(N-2)(q-1)}{l}} \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr) ^{\frac{2N}{N-2}}{\,\mathrm{d}}x \biggr)^{\frac{(N-2)(q-1)}{l}} \\& \quad = \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1-\frac{(N-2)(q-1)}{l}} \biggl( \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr)^{\frac{2N}{N-2}} {\,\mathrm{d}}x \biggr)^{\frac{N-2}{2N}} \biggr)^{\frac{2N(q-1)}{l}} \\& \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{(N-2)(q-1)}{l}} \biggl(C \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}} \biggr) ^{\frac{2N(q-1)}{l}} \\& \quad = C^{\frac{2N(q-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{(N-2)(q-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{N(q-1)}{l}} \\& \quad \leq C^{\frac{2N(q-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1-\frac{(N-2)(q-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{N(q-1)}{l}}+ \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr) ^{\frac{N(q-1)}{l}} \biggr) \\& \quad = C^{\frac{2N(q-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+\frac{2(q-1)}{l}} + \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{(N-2)(q-1)}{l}} \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr) ^{\frac{N(q-1)}{l}} \biggr) \\& \quad = C^{\frac{2N(q-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+\frac{2(q-1)}{l}} \\& \qquad {}+ \biggl(\varepsilon _{2}^{ \frac{N(1-q)}{l-N(q-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+ \frac{2(q-1)}{l-N(q-1)}} \biggr)^{1-\frac{N(q-1)}{l}} \biggl(\varepsilon _{2} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{ \frac{N(q-1)}{l}} \biggr) \\& \quad \leq C^{\frac{2N(q-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+\frac{2(q-1)}{l}} + \biggl(1-\frac{N(q-1)}{l} \biggr)\varepsilon _{2}^{\frac{N(1-q)}{l-N(q-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+\frac{2(q-1)}{l-N(q-1)}} \\& \qquad {}+\frac{N(q-1)^{\frac{}{}}}{l}\varepsilon _{2} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr), \end{aligned}$$
(4.15)

where \(\varepsilon _{2}=\frac{\theta (l-1)}{I_{1}N(q-1)}\), \(I_{1}=l \sigma C^{\frac{2N(q-1)}{l}} (\frac{N^{2}}{2\rho _{0}^{2}}+\frac{(l+q-1)^{2}d ^{2}}{2\varepsilon _{1}\rho _{0}^{2}} ) \).

Then we deal with the fourth term of (4.14). Due to (4.6), we get \(0<\frac{(m-1)(N-2)}{2l}<1\), \(0<\frac{N(m-1)}{2l}<1\). Using the Hölder inequality, the Young inequality, and (4.3)–(4.4), we derive

$$\begin{aligned}& \int _{D}u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{p} \\ & \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{(m-1)(N-2)}{2l}} \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr) ^{\frac{2N}{N-2}}{\,\mathrm{d}}x \biggr)^{\frac{(m-1)(N-2)}{2l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{p} \\ & \quad = \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p- \frac{(m-1)(N-2)}{2l}} \biggl( \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr) ^{\frac{2N}{N-2}}{\,\mathrm{d}}x \biggr)^{\frac{N-2}{2N}} \biggr)^{ \frac{N(m-1)}{l}} \\ & \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p- \frac{(m-1)(N-2)}{2l}} \biggl(C \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}} \biggr) ^{\frac{N(m-1)}{l}} \\ & \quad = C^{\frac{N(m-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p- \frac{(m-1)(N-2)}{2l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{N(m-1)}{2l}} \\ & \quad \leq C^{\frac{N(m-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+p-\frac{(m-1)(N-2)}{2l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{N(m-1)}{2l}}+ \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr) ^{\frac{N(m-1)}{2l}} \biggr) \\ & \quad = C^{\frac{N(m-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+p+\frac{m-1}{l}}+ \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p- \frac{(m-1)(N-2)}{2l}} \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2} {\,\mathrm{d}}x \biggr)^{\frac{N(m-1)}{2l}} \biggr) \\ & \quad = C^{\frac{N(m-1)}{l}} \biggl( \biggl(\varepsilon _{3} \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{\frac{l+r-1}{l}} \biggr)^{\frac{l}{l+r-1}} \biggl(\varepsilon _{3}^{\frac{l}{1-r}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} \biggr)^{\frac{r-1}{l+r-1}} \\ & \qquad {}+ \biggl(\varepsilon _{4}^{\frac{N(1-m)}{2l-N(m-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}} \biggr) ^{1-\frac{N(m-1)}{2l}} \biggl(\varepsilon _{4} \int _{D} \bigl\vert \nabla u^{ \frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{N(m-1)}{2l}} \biggr) \\ & \quad \leq C^{\frac{N(m-1)}{l}} \biggl(\frac{l}{l+r-1}\varepsilon _{3} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{\frac{l+r-1}{l}}+\frac{r-1}{l+r-1} \varepsilon _{3}^{\frac{l}{1-r}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} \\ & \qquad {}+ \biggl(1-\frac{N(m-1)}{2l} \biggr)\varepsilon _{4}^{ \frac{N(1-m)}{2l-N(m-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}} \\ & \qquad {}+\frac{N(m-1)}{2l}\varepsilon _{4} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr), \end{aligned}$$
(4.16)

where \(\varepsilon _{3}=\frac{\gamma _{2}(l+r-1)}{l\gamma _{1}}C^{ \frac{N(1-m)}{l}}|D|^{\frac{1-r}{l}}\), \(\varepsilon _{4}=\frac{2\theta (l-1)}{l\gamma _{1}N(m-1)}C^{\frac{N(1-m)}{l}}\).

Finally we compute the last term of (4.14). Applying the Hölder inequality, we deduce

$$ \int _{D}u^{l}{\,\mathrm{d}}x\leq \vert D \vert ^{\frac{r-1}{l+r-1}} \biggl( \int _{D}u ^{l+r-1}{\,\mathrm{d}}x \biggr)^{\frac{l}{l+r-1}}; $$

that is,

$$ \int _{D}u^{l+r-1}{\,\mathrm{d}}x\geq \vert D \vert ^{\frac{1-r}{l}} \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{\frac{l+r-1}{l}}. $$
(4.17)

Substituting (4.15)–(4.17) into (4.14), we have

$$\begin{aligned} A'(t) \leq &\frac{l\sigma }{2} \int _{D}u^{l}{\,\mathrm{d}}x+I_{1} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2(q-1)}{l}} \\ &{}+ \frac{I_{1}[l-N(q-1)]}{l}\varepsilon _{2}^{\frac{N(1-q)}{l-N(q-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2(q-1)}{l-N(q-1)}} \\ &{}+\frac{l\gamma _{1}(r-1)}{l+r-1}C^{\frac{N(m-1)}{l}}\varepsilon _{3} ^{\frac{l}{1-r}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}} \\ &{}+\frac{\gamma _{1}[2l-N(m-1)]}{2}C^{\frac{N(m-1)}{l}}\varepsilon _{4} ^{\frac{N(1-m)}{2l-N(m-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}}. \end{aligned}$$
(4.18)

By (4.1), we obtain

$$ A(t)=l \int _{D} \int _{0}^{u}y^{l-1}g'(y) \,\mathrm{d}y\,\mathrm{d}x\geq \xi \int _{D}u^{l}{\,\mathrm{d}}x, $$

which implies

$$ \int _{D}u^{l}{\,\mathrm{d}}x\leq \frac{1}{\xi }A(t). $$
(4.19)

Now, inserting (4.19) into (4.18), we derive

$$\begin{aligned} A'(t) \leq &\frac{l\sigma }{2\xi }A(t) +\frac{I_{1}}{\xi ^{1+ \frac{2(q-1)}{l}}} A^{1+\frac{2(q-1)}{l}}(t) + \frac{I_{1} [l-N(q-1) ]}{l \xi ^{1+\frac{2(q-1)}{l-N(q-1)}}}\varepsilon _{2}^{ \frac{N(1-q)}{l-N(q-1)}} A^{1+\frac{2(q-1)}{l-N(q-1)}}(t) \\ &{}+ \frac{l\gamma _{1}(r-1)}{(l+r-1)\xi ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}}}C ^{\frac{N(m-1)}{l}}\varepsilon _{3}^{\frac{l}{1-r}}A^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}}(t) \\ &{}+ \frac{\gamma _{1}[2l-N(m-1)]}{2\xi ^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}}}C ^{\frac{N(m-1)}{l}}\varepsilon _{4}^{\frac{N(1-m)}{2l-N(m-1)}}A^{1+ \frac{2pl+2(m-1)}{2l-N(m-1)}}(t) \\ =& J_{1}A(t)+J_{2}A^{1+\frac{2(q-1)}{l}}(t)+J_{3}A^{1+ \frac{2(q-1)}{l-N(q-1)}}(t) \\ &{}+J_{4}A^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}}(t)+J _{5}A^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}}(t), \end{aligned}$$
(4.20)

where \(J_{1}\), \(J_{2}\), \(J_{3}\), \(J_{4}\), and \(J_{5}\) are defined in (4.7)–(4.10). Since \(\lim_{t\to t^{*-}}A(t)=+\infty \), we integrate (4.20) from 0 to \(t^{*}\) to get

$$ t^{*}\geq \int _{A(0)}^{+\infty }\frac{\mathrm{d}\tau }{J_{1}\tau +J _{2}\tau ^{1+\frac{2(q-1)}{l}}+J_{3}\tau ^{1+\frac{2(q-1)}{l-N(q-1)}}+J _{4}\tau ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} +J_{5} \tau ^{1+\frac{2pl+2(m-1)}{2l-N(m-1)}}}. $$

 □

5 Lower bounds of the blow-up time \(t^{*}\) in \(D\subset \mathbb{R}^{2}\)

In this section, we search for a lower bound of \(t^{*}\) in \(D\subset \mathbb{R}^{2}\). Suppose that the functions h and g satisfy

$$ h(s)\leq \tilde{\sigma } s^{\tilde{q}},\qquad g'(s)\geq \tilde{\xi },\quad s\geq 0, $$
(5.1)

where σ̃, , and ξ̃ are some positive constants and

$$ \tilde{q}>1. $$
(5.2)

By [17, Corollary 9.14], we have

$$ W^{1, 2}\hookrightarrow L^{4},\quad \mbox{where } N=2, $$

which implies the following Sobolev inequality:

$$ \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr)^{4}{ \,\mathrm{d}}x \biggr) ^{\frac{1}{4}} \leq \tilde{C} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}}, $$
(5.3)

where is a Sobolev embedding constant depending on D. We need to use inequality (5.3) to get a lower bound of \(t^{*}\). Now, we define auxiliary functions as follows:

$$ \tilde{A}(t)= \int _{D} \tilde{B}(u)\,\mathrm{d}x, \quad t\geq 0, \qquad \tilde{B}(s)=l \int _{0}^{s}y^{l-1}g'(y) \,\mathrm{d}y,\quad s\geq 0. $$
(5.4)

The main result of this section is the following Theorem 5.1.

Theorem 5.1

Let\(u({x},t)\)be a nonnegative classical solution of (1.1) in\(D\subset \mathbb{R}^{2}\). Assume that (5.1)(5.2) hold and

$$ r>1,\qquad m>1,\qquad l>\max \bigl\{ 2, 2(m-1), 4( \tilde{q}-1) \bigr\} . $$
(5.5)

In addition, we assume that\(u({x},t)\)blows up at some finite time\(t^{*}\)in the measure\(\tilde{A}(t)\). Then\(t^{*}\)is bounded below by

$$ t^{*}\geq \int _{\tilde{A}(0)}^{+\infty }\frac{\mathrm{d}\tau }{ \tilde{J}_{1}\tau +\tilde{J}_{2}\tau ^{1+\frac{2(\tilde{q}-1)}{l}}+ \tilde{J}_{3}\tau ^{1+\frac{2(\tilde{q}-1)}{l-4(\tilde{q}-1)}}+ \tilde{J}_{4}\tau ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} +\tilde{J}_{5} \tau ^{\frac{l(1+p)-(m-1)}{l-2(m-1)}}}, $$

where

$$\begin{aligned}& \tilde{J}_{1}=\frac{l\tilde{\sigma }}{2\tilde{\xi }},\qquad \tilde{J}_{2}=\frac{l\tilde{\sigma }}{\rho _{0}^{2}\tilde{\xi }^{1+\frac{2( \tilde{q}-1)}{l}}}\tilde{C}^{\frac{8(\tilde{q}-1)}{l}} \biggl(2 +\frac{ \tilde{\sigma }(l+\tilde{q}-1)^{2}\tilde{d}^{2}}{2\theta (l-1)} \biggr), \end{aligned}$$
(5.6)
$$\begin{aligned}& \tilde{J}_{3}=\frac{l-4(\tilde{q}-1)}{\tilde{\xi }^{1+\frac{2( \tilde{q}-1)}{l-4(q-1)}}} \biggl( \frac{\theta (l-1)}{4l(\tilde{q}-1) \tilde{C}^{2}} \biggr)^{\frac{4(1-\tilde{q})}{l-4(\tilde{q}-1)}} \biggl[\frac{ \tilde{\sigma }}{\tilde{\rho }_{0}^{2}} \biggl(2+\frac{\tilde{\sigma } (l+\tilde{q}-1)^{2}\tilde{d}^{2}}{2\theta (l-1)} \biggr) \biggr]^{\frac{l}{l-4( \tilde{q}-1)}}, \end{aligned}$$
(5.7)
$$\begin{aligned}& \tilde{J}_{4}=\frac{ \vert D \vert (r-1)}{\tilde{\xi }^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}}}\gamma _{2}^{\frac{l}{1-r}} \biggl(\frac{l \gamma _{1}}{l+r-1} \tilde{C}^{\frac{4(m-1)}{l}} \biggr)^{1+ \frac{l}{r-1}}, \end{aligned}$$
(5.8)
$$\begin{aligned}& \tilde{J}_{5}=\frac{l-2(m-1)}{\tilde{\xi }^{ \frac{l(1+p)-(m-1)}{l-2(m-1)}}}\gamma _{1}^{\frac{l}{l-2(m-1)}} \biggl(\frac{ \theta (l-1)}{2l(m-1)\tilde{C}^{2}} \biggr)^{\frac{2(1-m)}{l-2(m-1)}}, \end{aligned}$$
(5.9)
$$\begin{aligned}& \tilde{\rho }_{0}=\min_{{x}\in {\partial D}}(x \cdot \nu ), \qquad \tilde{d}=\max_{x\in \overline{D}} \vert x \vert , \end{aligned}$$
(5.10)

andθis defined in (1.2).

Proof

Repeating the calculation process of (4.12)–(4.14) and noting \(N=2\) in this section, we get

$$\begin{aligned} \tilde{A}'(t) \leq & \frac{l\tilde{\sigma }}{2} \int _{D}u^{l}{\,\mathrm{d}}x+l \tilde{\sigma } \biggl(\frac{2}{\tilde{\rho }_{0}^{2}} +\frac{(l+ \tilde{q}-1)^{2}\tilde{d}^{2}}{2\tilde{\varepsilon }_{1}\tilde{\rho } _{0}^{2}} \biggr) \int _{D}u^{l+2\tilde{q}-2}{\,\mathrm{d}}x \\ &{}+ \biggl( \frac{2 \tilde{\varepsilon }_{1}\tilde{\sigma }}{l}-\frac{4\theta (l-1)}{l} \biggr) \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \\ &{}+l\gamma _{1} \int _{D} u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{p}-l\gamma _{2} \int _{D} u^{l+r-1}{\,\mathrm{d}}x, \end{aligned}$$
(5.11)

where \(\tilde{\rho }_{0}\) and are defined in (5.10) and \(\tilde{\varepsilon }_{1}=\frac{\theta (l-1)}{\tilde{\sigma }}\).

First, we focus on the second term of (5.11). By (5.2) and (5.5), we get \(0<\frac{2(\tilde{q}-1)}{l}<1\) and \(0< \frac{4(\tilde{q}-1)}{l}<1\). Now we apply the Hölder inequality, the Young inequality, (4.4), and (5.3) to deduce

$$\begin{aligned}& \int _{D}u^{l+2\tilde{q}-2}{\,\mathrm{d}}x \\& \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{2(\tilde{q}-1)}{l}} \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr) ^{4}{\,\mathrm{d}}x \biggr)^{\frac{2(\tilde{q}-1)}{l}} \\& \quad = \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{2(\tilde{q}-1)}{l}} \biggl( \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr) ^{4}{\,\mathrm{d}}x \biggr)^{\frac{1}{4}} \biggr)^{ \frac{8(\tilde{q}-1)}{l}} \\& \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1- \frac{2(\tilde{q}-1)}{l}} \biggl(\tilde{C} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{ \frac{1}{2}} \biggr)^{\frac{8(\tilde{q}-1)}{l}} \\& \quad = \tilde{C}^{\frac{8(\tilde{q}-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1-\frac{2(\tilde{q}-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{4( \tilde{q}-1)}{l}} \\& \quad \leq \tilde{C}^{\frac{8(\tilde{q}-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1-\frac{2(\tilde{q}-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{4(\tilde{q}-1)}{l}}+ \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2} {\,\mathrm{d}}x \biggr)^{\frac{4(\tilde{q}-1)}{l}} \biggr) \\& \quad = \tilde{C}^{\frac{8(\tilde{q}-1)}{l}} \biggl( \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{1+\frac{2(\tilde{q}-1)}{l}} + \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{1-\frac{2(\tilde{q}-1)}{l}} \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{ \frac{4(\tilde{q}-1)}{l}} \biggr) \\& \quad = \tilde{C}^{\frac{8(\tilde{q}-1)}{l}} \biggl( \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{1+\frac{2(\tilde{q}-1)}{l}} \\& \qquad {}+ \biggl(\tilde{\varepsilon } _{2}^{\frac{4(1-\tilde{q})}{l-4(\tilde{q}-1)}} \biggl( \int _{D}u^{l} {\,\mathrm{d}}x \biggr)^{1+\frac{2(\tilde{q}-1)}{l-4(\tilde{q}-1)}} \biggr) ^{1-\frac{4(\tilde{q}-1)}{l}} \biggl(\tilde{\varepsilon }_{2} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{ \frac{4(\tilde{q}-1)}{l}} \biggr) \\& \quad \leq \tilde{C}^{\frac{8(\tilde{q}-1)}{l}} \biggl( \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2(\tilde{q}-1)}{l}} + \biggl(1-\frac{4( \tilde{q}-1)}{l} \biggr)\tilde{ \varepsilon }_{2}^{\frac{4(1-\tilde{q})}{l-4( \tilde{q}-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2( \tilde{q}-1)}{l-4(\tilde{q}-1)}} \\& \qquad {}+\frac{4(\tilde{q}-1)^{\frac{}{}}}{l}\tilde{\varepsilon } _{2} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr), \end{aligned}$$
(5.12)

where \(\tilde{\varepsilon }_{2}=\frac{\theta (l-1)}{4\tilde{I}_{1}( \tilde{q}-1)}\), \(\tilde{I}_{1}=l\tilde{\sigma }\tilde{C}^{\frac{8( \tilde{q}-1)}{l}} (\frac{2}{\tilde{\rho }_{0}^{2}} +\frac{(l+ \tilde{q}-1)^{2}\tilde{d}^{2}}{2\tilde{\varepsilon }_{1}\tilde{\rho } _{0}^{2}} ) \).

Next, we deal with the fourth term of (5.11). Due to (5.5), we get \(0<\frac{m-1}{l}<1\), \(0<\frac{2(m-1)}{l}<1\). Using the Hölder inequality, the Young inequality, (4.4), and (5.3), we derive

$$\begin{aligned}& \int _{D}u^{l+m-1}{\,\mathrm{d}}x \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{p} \\& \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1-\frac{m-1}{l}} \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr)^{4}{ \,\mathrm{d}}x \biggr) ^{\frac{m-1}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{p} \\& \quad = \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p-\frac{m-1}{l}} \biggl( \biggl( \int _{D} \bigl(u^{\frac{l}{2}} \bigr)^{4}{ \,\mathrm{d}}x \biggr)^{ \frac{1}{4}} \biggr)^{\frac{4(m-1)}{l}} \\& \quad \leq \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p-\frac{m-1}{l}} \biggl(\tilde{C} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u ^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{1}{2}} \biggr)^{ \frac{4(m-1)}{l}} \\& \quad = \tilde{C}^{\frac{4(m-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+p-\frac{m-1}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x+ \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{\frac{2(m-1)}{l}} \\& \quad \leq \tilde{C}^{\frac{4(m-1)}{l}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+p-\frac{m-1}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{2(m-1)}{l}}+ \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr) ^{\frac{2(m-1)}{l}} \biggr) \\& \quad = \tilde{C}^{\frac{4(m-1)}{l}} \biggl( \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{1+p+\frac{m-1}{l}}+ \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+p- \frac{m-1}{l}} \biggl( \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr) ^{\frac{2(m-1)}{l}} \biggr) \\& \quad = \tilde{C}^{\frac{4(m-1)}{l}} \biggl( \biggl(\tilde{\varepsilon }_{3} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{\frac{l+r-1}{l}} \biggr) ^{\frac{l}{l+r-1}} \biggl(\tilde{\varepsilon }_{3}^{\frac{l}{1-r}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}} \biggr)^{\frac{r-1}{l+r-1}} \\& \qquad {}+ \biggl(\tilde{\varepsilon }_{4}^{\frac{2(1-m)}{l-2(m-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{ \frac{l(1+p)-(m-1)}{l-2(m-1)}} \biggr)^{1-\frac{2(m-1)}{l}} \biggl(\tilde{\varepsilon } _{4} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr)^{ \frac{2(m-1)}{l}} \biggr) \\& \quad \leq \tilde{C}^{\frac{4(m-1)}{l}} \biggl(\frac{l}{l+r-1} \tilde{ \varepsilon }_{3} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{ \frac{l+r-1}{l}}+\frac{r-1}{l+r-1}\tilde{\varepsilon }_{3}^{ \frac{l}{1-r}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}} \\& \qquad {}+ \biggl(1-\frac{2(m-1)}{l} \biggr)\tilde{\varepsilon }_{4} ^{\frac{2(1-m)}{l-2(m-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{l(1+p)-(m-1)}{l-2(m-1)}} \\& \qquad {}+\frac{2(m-1)}{l}\tilde{\varepsilon } _{4} \int _{D} \bigl\vert \nabla u^{\frac{l}{2}} \bigr\vert ^{2}{\,\mathrm{d}}x \biggr), \end{aligned}$$
(5.13)

where \(\tilde{\varepsilon }_{3}=\frac{\gamma _{2}(l+r-1)}{l\gamma _{1}} \tilde{C}^{\frac{4(1-m)}{l}}|D|^{\frac{1-r}{l}}\), \(\tilde{\varepsilon }_{4}=\frac{\theta (l-1)}{2l\gamma _{1}(m-1)} \tilde{C}^{\frac{4(1-m)}{l}}\).

Finally, we consider the last term of (5.11). Repeating the calculation process of (4.17), we have

$$ \int _{D}u^{l+r-1}{\,\mathrm{d}}x\geq \vert D \vert ^{\frac{1-r}{l}} \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{\frac{l+r-1}{l}}. $$
(5.14)

Substituting (5.12)–(5.14) into (5.11), we deduce

$$\begin{aligned} \tilde{A}'(t) \leq &\frac{l\tilde{\sigma }}{2} \int _{D}u^{l}{\,\mathrm{d}}x+ \tilde{I}_{1} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2( \tilde{q}-1)}{l}} \\ &{}+\frac{\tilde{I}_{1}[l-4(\tilde{q}-1)]}{l} \tilde{\varepsilon }_{2}^{\frac{4(1-\tilde{q})}{l-4(\tilde{q}-1)}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr)^{1+\frac{2(\tilde{q}-1)}{l-4( \tilde{q}-1)}} \\ &{}+\frac{l\gamma _{1}(r-1)}{l+r-1}\tilde{C}^{\frac{4(m-1)}{l}} \tilde{\varepsilon }_{3}^{\frac{l}{1-r}} \biggl( \int _{D}u^{l}{\,\mathrm{d}}x \biggr) ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} \\ &{}+\gamma _{1}\bigl[l-2(m-1)\bigr]\tilde{C}^{\frac{4(m-1)}{l}} \tilde{\varepsilon }_{4}^{\frac{2(1-m)}{l-2(m-1)}} \biggl( \int _{D}u ^{l}{\,\mathrm{d}}x \biggr)^{\frac{l(1+p)-(m-1)}{l-2(m-1)}}. \end{aligned}$$
(5.15)

We repeat the calculation process of (4.19) to obtain

$$ \int _{D}u^{l}{\,\mathrm{d}}x\leq \frac{1}{\tilde{\xi }}\tilde{A}(t). $$
(5.16)

Now, inserting (5.16) into (5.15), we derive

$$\begin{aligned} \tilde{A}'(t) \leq &\frac{l\tilde{\sigma }}{2\tilde{\xi }}\tilde{A}(t) +\frac{ \tilde{I}_{1}}{\tilde{\xi }^{1+\frac{2(\tilde{q}-1)}{l}}} \tilde{A} ^{1+\frac{2(\tilde{q}-1)}{l}}(t) +\frac{\tilde{I}_{1} [l-4( \tilde{q}-1) ]}{l\tilde{\xi }^{1+\frac{2(\tilde{q}-1)}{l-4( \tilde{q}-1)}}}\tilde{\varepsilon }_{2}^{\frac{4(1-\tilde{q})}{l-4( \tilde{q}-1)}} \tilde{A}^{1+\frac{2(\tilde{q}-1)}{l-4(\tilde{q}-1)}}(t) \\ &{}+\frac{l\gamma _{1}(r-1)}{(l+r-1)\tilde{\xi }^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}}}\tilde{C}^{\frac{4(m-1)}{l}} \tilde{\varepsilon }_{3}^{\frac{l}{1-r}}\tilde{A}^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}}(t) \\ &{}+\frac{\gamma _{1}[l-2(m-1)]}{\tilde{\xi }^{ \frac{l(1+p)-(m-1)}{l-2(m-1)}}}\tilde{C}^{\frac{4(m-1)}{l}} \tilde{\varepsilon }_{4}^{\frac{2(1-m)}{l-2(m-1)}} \tilde{A}^{ \frac{l(1+p)-(m-1)}{l-2(m-1)}}(t) \\ =& \tilde{J}_{1}\tilde{A}(t)+\tilde{J}_{2} \tilde{A}^{1+\frac{2( \tilde{q}-1)}{l}}(t)+\tilde{J}_{3}\tilde{A}^{1+\frac{2(\tilde{q}-1)}{l-4( \tilde{q}-1)}}(t) \\ &{}+ \tilde{J}_{4}\tilde{A}^{ \frac{(lp+m-1)(l+r-1)}{l(r-1)}}(t) +\tilde{J}_{5} \tilde{A}^{ \frac{l(1+p)-(m-1)}{l-2(m-1)}}(t), \end{aligned}$$
(5.17)

where \(\tilde{J}_{1}\), \(\tilde{J}_{2}\), \(\tilde{J}_{3}\), \(\tilde{J} _{4}\), and \(\tilde{J}_{5}\) are defined in (5.6)–(5.9). Since \(\lim_{t\to t^{*-}}\tilde{A}(t)=+\infty \), we integrate (5.17) from 0 to \(t^{*}\) to get

$$ t^{*}\geq \int _{\tilde{A}(0)}^{+\infty }\frac{\mathrm{d}\tau }{ \tilde{J}_{1}\tau +\tilde{J}_{2}\tau ^{1+\frac{2(\tilde{q}-1)}{l}}+ \tilde{J}_{3}\tau ^{1+\frac{2(\tilde{q}-1)}{l-4(\tilde{q}-1)}}+ \tilde{J}_{4}\tau ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} +\tilde{J}_{5} \tau ^{\frac{l(1+p)-(m-1)}{l-2(m-1)}}}. $$

 □

6 Applications

In what follows, three examples are given to demonstrate the results of Theorems 2.15.1 obtained in this paper.

Example 6.1

Let \(u(x,t)\) be a nonnegative classical solution of the following problem:

$$ \textstyle\begin{cases} ( (u+\ln (1+u) ) )_{t} =\sum_{i=1}^{3} ( (1+ \vert x \vert ^{2} )u _{{x}_{i}} )_{{ x}_{i}}+u^{\frac{1}{2}} (\int _{D} u^{3} {\,\mathrm{d}}x )^{\frac{1}{3}}-u^{\frac{5}{2}} &\mbox{in } D \times (0,t^{*}), \\ \sum_{i=1}^{3} (1+ \vert x \vert ^{2} )u_{{x}_{i}}\nu _{i}=\frac{2 \sqrt{5}}{5}u^{\frac{3}{2}} &\mbox{on } \partial D\times (0,t^{*}), \\ u(x,0)=u_{0}(x)=1+ \vert x \vert ^{2} &\mbox{in } \overline{D}, \end{cases} $$

where \(D= \{ x=({x}_{1}, {x}_{2}, {x}_{3}) \mid |x|^{2}=\sum_{i=1}^{3}{x}_{i}^{2} <\frac{1}{4} \} \) is a ball of \(\mathbb{R}^{3}\). Here,

$$\begin{aligned}& g(u)=u+\ln (1+u), \qquad h(u)=\frac{2\sqrt{5}}{5}u^{\frac{3}{2}}, \\& \bigl(a^{ij}(x) \bigr)_{3\times 3}= \begin{pmatrix} 1+ \vert x \vert ^{2} & 0 & 0 \\ 0 & 1+ \vert x \vert ^{2} & 0 \\ 0 & 0& 1+ \vert x \vert ^{2} \end{pmatrix},\qquad u_{0}(x)=1+ \vert x \vert ^{2}. \end{aligned}$$

Now, \(N=3\), \(\gamma _{1}=1\), \(\gamma _{2}=1\), \(m=\frac{1}{2}\), \(l=3\), \(p= \frac{1}{3}\), and \(r=\frac{5}{2}\). Choosing \(\zeta =2\), \(\gamma =\frac{2\sqrt{5}}{5}\), and \(q=\frac{3}{2}\), we can verify that (2.2)–(2.3) are satisfied. Hence, Theorem 2.1 implies that u exists for all time in \(\varPhi (t)\) with

$$ \varPhi (t)= \int _{D} G(u)\,\mathrm{d}x= \int _{D} \biggl(u^{3}+\frac{3}{2}u ^{2}-3u+3\ln (1+u) \biggr)\,\mathrm{d}x. $$

Example 6.2

Let \(u(x,t)\) be a nonnegative classical solution of the following problem:

$$ \textstyle\begin{cases} (\frac{u+\ln (1+u)}{2} )_{t} =\sum_{i=1}^{3} ( (1+ \vert x \vert ^{2} )u _{{x}_{i}} )_{{ x}_{i}}+\frac{9}{4}u^{\frac{3}{2}} (\int _{D} u^{3}{\,\mathrm{d}}x )^{\frac{1}{10}}-\frac{7}{4}u^{ \frac{7}{5}} &\mbox{in } D\times (0,t^{*}), \\ \sum_{i=1}^{3} (1+ \vert x \vert ^{2} )u_{{x}_{i}}\nu _{i}=\frac{3 \sqrt{5}}{5}u^{\frac{5}{4}} &\mbox{on } \partial D\times (0,t^{*}), \\ u(x,0)=u_{0}(x)=1+ \vert x \vert ^{2} &\mbox{in } \overline{D}, \end{cases} $$

where \(D= \{ x=({x}_{1}, {x}_{2}, {x}_{3}) \vert |x|^{2}=\sum_{i=1}^{3}{x}_{i}^{2} . <\frac{9}{16} \} \) is a ball of \(\mathbb{R}^{3}\). Now we have

$$\begin{aligned}& g(u)=\frac{u+\ln (1+u)}{2},\qquad h(u)=\frac{3\sqrt{5}}{5}u^{\frac{5}{4}}, \\& \bigl(a^{ij}(x) \bigr)_{3\times 3}= \begin{pmatrix} 1+ \vert x \vert ^{2} & 0 & 0 \\ 0 & 1+ \vert x \vert ^{2} & 0 \\ 0 & 0& 1+ \vert x \vert ^{2} \end{pmatrix},\qquad u_{0}(x)=1+ \vert x \vert ^{2}. \end{aligned}$$

Here, \(N=3\), \(\gamma _{1}=\frac{9}{4}\), \(\gamma _{2}=\frac{7}{4}\), \(m= \frac{3}{2}\), \(l=3\), \(p=\frac{1}{10}\), \(r=\frac{7}{5}\), and \(|D|=\frac{9 \pi }{16}\). Now, we compute

$$\begin{aligned}& \varPsi (t)= \int _{D} g(u)\,\mathrm{d}x=\frac{1}{2} \int _{D} \bigl(u+\ln (1+u) \bigr) \,\mathrm{d}x, \\& \varPsi (0)=\frac{1}{2} \int _{D} \bigl(u_{0}+\ln (1+u_{0}) \bigr) \,\mathrm{d}x=1.9302. \end{aligned}$$

It follows from (3.4) that \(\delta =0.7436\), \(M_{1}=0.7552\), and \(M_{2}=1.0969\). It is easy to check that (3.2)–(3.3) hold. By Theorem 3.1, we deduce that u must blow up in the measure \(\varPsi (t)\) at some finite time \(t^{*}\) and

$$ t^{*}\leq \int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M _{2}}= \int _{1.9302}^{+\infty } \frac{d\tau }{0.7552\tau ^{\frac{9}{5}}-1.0969}=1.2304. $$
(6.1)

We note that (6.1) gives an upper bound of \(t^{*}\). Next, we use Theorem 4.1 to obtain a lower bound of \(t^{*}\). Now we select \(q=\frac{5}{4}\), \(\sigma =\frac{3\sqrt{5}}{5}\), \(\xi =\frac{1}{2}\), and \(\theta =1\). By (4.11), we derive \(\rho _{0}=d=\frac{3}{4}\). Noting \(N=3\), we can verify that (4.1)–(4.2) and (4.6) hold. Moreover, according to Theorems 2.1 and 3.2 in [18], we obtain the Sobolev embedding constant \(C=7.5931\). Substituting above constants into (4.7)–(4.10), it is easy to get \(J_{1}= 4.0249\), \(J_{2}=287.4\), \(J_{3}=814.14\), \(J_{4}=95979\), and \(J_{5}=68.205\). By (3.6), we obtain

$$\begin{aligned}& B(u)=l \int _{0}^{u}y^{l-1}g'(y) \,\mathrm{d}y=\frac{3}{2} \int _{0}^{u}y ^{2} \biggl(1+ \frac{1}{1+y} \biggr)\,\mathrm{d}y=\frac{1}{2}u^{3}+ \frac{3}{4}u^{2}-\frac{3}{2}u+ \frac{3}{2}\ln (1+u), \\& A(t)= \int _{D} B(u)\,\mathrm{d}x= \int _{D} \biggl(\frac{1}{2}u^{3}+ \frac{3}{4}u^{2}-\frac{3}{2}u+ \frac{3}{2}\ln (1+u) \biggr)\,\mathrm{d}x, \end{aligned}$$

and

$$ A(0)= \int _{D} \biggl(\frac{1}{2}u_{0}^{3}+ \frac{3}{4}u_{0}^{2}- \frac{3}{2}u_{0}+ \frac{3}{2}\ln (1+u_{0}) \biggr)\,\mathrm{d}x =3.2897. $$

Blow-up of u in the measure \(\varPsi (t)\) at \(t^{*}\) means that u blows up at \(t^{*}\). Hence u must also blow up in the measure \(A(t)\) at \(t^{*}\). By Theorem 4.1, we get a lower bound of the blow-up time \(t^{*}\) as follows:

$$\begin{aligned} t^{*} \geq & T_{2}= \int _{A(0)}^{+\infty }\frac{\mathrm{d}\tau }{J_{1} \tau +J_{2}\tau ^{1+\frac{2(q-1)}{l}}+J_{3}\tau ^{1+ \frac{2(q-1)}{l-n(q-1)}}+J_{4}\tau ^{1+p+\frac{m-1}{l}} +J_{5} \tau ^{1+\frac{2pl+2(m-1)}{2l-n(m-1)}}} \\ =& \int _{3.2897}^{+\infty }\frac{\mathrm{d}\tau }{4.0249\tau +287.4 \tau ^{\frac{7}{6}}+814.14\tau ^{\frac{11}{9}}+95979\tau ^{\frac{34}{15}} +68.205\tau ^{\frac{61}{45}}} \\ =&1.8166\times 10^{-6}. \end{aligned}$$
(6.2)

Combining (6.1) and (6.2), we have

$$ 1.8166\times 10^{-6}\leq t^{*}\leq 1.2304. $$

Example 6.3

Let \(u(x,t)\) be a nonnegative classical solution of the following problem:

$$ \textstyle\begin{cases} u_{t} =\sum_{i=1}^{2} ( (1+ \vert x \vert ^{2} )u_{{x}_{i}} ) _{{ x}_{i}}+2u^{\frac{3}{2}} (\int _{D} u^{3}{\,\mathrm{d}}x ) ^{\frac{1}{10}}-2u^{\frac{7}{5}} &\mbox{in } D\times (0,t^{*}), \\ \sum_{i=1}^{2} (1+ \vert x \vert ^{2} )u_{{x}_{i}}\nu _{i}=\frac{3 \sqrt{5}}{5}u^{\frac{5}{4}} &\mbox{on } \partial D\times (0,t^{*}), \\ u(x,0)=u_{0}(x)=1+ \vert x \vert ^{2} &\mbox{in } \overline{D}, \end{cases} $$

where \(D= \{ x=({x}_{1}, {x}_{2}) \mid |x|^{2}=\sum_{i=1}^{2} {x}_{i}^{2} <\frac{9}{16} \} \) is a circular region of \(\mathbb{R}^{2}\). Now

$$\begin{aligned}& g(u)=u,\qquad h(u)=\frac{3\sqrt{5}}{5}u^{\frac{5}{4}}, \\& \bigl(a^{ij}(x) \bigr)_{2\times 2}= \begin{pmatrix} 1+ \vert x \vert ^{2} & 0 \\ 0 & 1+ \vert x \vert ^{2} \end{pmatrix},\qquad u_{0}(x)=1+ \vert x \vert ^{2}. \end{aligned}$$

Here, \(N=2\), \(\gamma _{1}=2\), \(\gamma _{2}=2\), \(m=\frac{3}{2}\), \(l=3\), \(p= \frac{1}{10}\), \(r=\frac{7}{5}\), and \(|D|=\frac{9\pi }{16}\). By (3.1), we have

$$ \varPsi (t)= \int _{D} g(u)\,\mathrm{d}x= \int _{D} u\,\mathrm{d}x, \qquad \varPsi (0)= \int _{D} u_{0} {\,\mathrm{d}}x= \int _{D} \bigl(1+ \vert x \vert ^{2} \bigr) \,\mathrm{d}x=2.2642. $$

By (3.4), we obtain \(\delta =0.6506\), \(M_{1}=0.7552\), and \(M_{2}=2.0004\). After some simple calculations, we know that (3.2)–(3.3) hold. From Theorem 3.1, it follows that u blows up in the measure \(\varPsi (t)\) at \(t^{*}\) and

$$ t^{*}\leq \int _{\varPsi (0)}^{+\infty }\frac{d\tau }{M_{1}\tau ^{m+lp}-M _{2}} = \int _{2.2642}^{+\infty } \frac{d\tau }{0.7552\tau ^{\frac{9}{5}}-2.0004}=1.1291. $$
(6.3)

Next, we apply Theorem 5.1 to get a lower bound of \(t^{*}\). Now we choose \(\tilde{q}=\frac{5}{4}\), \(\tilde{\sigma }=\frac{3\sqrt{5}}{5}\), \(\tilde{\xi }=1\), and \(\theta =1\). It follows from (5.10) that \(\tilde{\rho }_{0}=\tilde{d}=\frac{3}{4}\). Noting \(N=2\), it is easy to verify that (5.1)–(5.2) and (5.5) hold. Moreover, using Theorems 2.1 and 3.2 in [18], we derive the Sobolev embedding constant \(\tilde{C}=10.887\). Inserting the above constants into (5.6)–(5.9), we deduce \(\tilde{J}_{1}= 2.0125\), \(\tilde{J}_{2}=140.34\), \(\tilde{J}_{3}=783.76\), \(\tilde{J}_{4}=997648\), and \(\tilde{J}_{5}=90\). By (5.4), we have

$$ \tilde{B}(u)=l \int _{0}^{u}y^{l-1}g'(y) \,\mathrm{d}y=3 \int _{0}^{u}y^{2} { \,\mathrm{d}}y=u^{3},\qquad \tilde{A}(t)= \int _{D} \tilde{B}(u)\,\mathrm{d}x= \int _{D} u^{3}{\,\mathrm{d}}x, $$

and

$$ \tilde{A}(0)= \int _{D} u_{0}^{3}{ \,\mathrm{d}}x= \int _{D} \bigl(1+ \vert x \vert ^{2} \bigr) ^{3}{\,\mathrm{d}}x =3.8959. $$

Since u blows up in the measure \(\varPsi (t)\) at \(t^{*}\), u must blow up in the measure \(\tilde{A}(t)\) at \(t^{*}\). By Theorem 5.1, we get

$$\begin{aligned} t^{*} \geq & \int _{\tilde{A}(0)}^{+\infty }\frac{\mathrm{d}\tau }{ \tilde{J}_{1}\tau +\tilde{J}_{2}\tau ^{1+\frac{2(\tilde{q}-1)}{l}}+ \tilde{J}_{3}\tau ^{1+\frac{2(\tilde{q}-1)}{l-4(\tilde{q}-1)}}+ \tilde{J}_{4}\tau ^{\frac{(lp+m-1)(l+r-1)}{l(r-1)}} +\tilde{J}_{5} \tau ^{\frac{l(1+p)-(m-1)}{l-2(m-1)}}} \\ =& \int _{3.8959}^{+\infty }\frac{\mathrm{d}\tau }{2.0125\tau +140.34 \tau ^{\frac{7}{6}}+783.76\tau ^{\frac{5}{4}}+997648\tau ^{\frac{34}{15}} +90\tau ^{\frac{7}{5}}} \\ =&1.4132\times 10^{-7}. \end{aligned}$$
(6.4)

From (6.3) and (6.4), it follows that

$$ 1.4132\times 10^{-7}\leq t^{*}\leq 1.1291. $$

7 Conclusion

In this paper, we derive the global existence and bounds for the blow-up time of nonlinear parabolic problem (1.1) with nonlocal source. To deal with nonlocal source, we must establish some new auxiliary functions different from those in [1, 2] and [11]. Furthermore, to obtain the lower bound of the blow-up time in \(D\subset \mathbb{R}^{N}\) (\(N\geq 3\)) and \(D\subset \mathbb{R} ^{2}\), we need to use the embedding theorems in Sobolev spaces \(W^{1, 2}\hookrightarrow L^{\frac{2N}{N-2}}\), \(N\geq 3 \) and \(W^{1, 2}\hookrightarrow L^{4}\), \(N=2\), respectively. Applying these auxiliary functions, the embedding theorems in Sobolev spaces, and the differential inequality technique, we complete our study with the blow-up and global solutions of problem (1.1).

References

  1. Li, F.S., Li, J.L.: Global existence and blow-up phenomena for nonlinear divergence form parabolic equations with inhomogeneous Neumann boundary conditions. J. Math. Anal. Appl. 385, 1005–1014 (2012)

    Article  MathSciNet  Google Scholar 

  2. Zhang, J.Z., Li, F.S.: Global existence and blow-up phenomena for divergence form parabolic equation with time-dependent coefficient in multidimensional space. Z. Angew. Math. Phys. 70, 1–16 (2019)

    Article  MathSciNet  Google Scholar 

  3. Ding, J.T., Kou, W.: Blow-up solutions for reaction diffusion equations with nonlocal boundary conditions. J. Math. Anal. Appl. 470, 1–15 (2019)

    Article  MathSciNet  Google Scholar 

  4. Ding, J.T., Shen, X.H.: Blow-up time estimates in porous medium equations with nonlinear boundary conditions. Z. Angew. Math. Phys. 69, 1–13 (2018)

    Article  MathSciNet  Google Scholar 

  5. Payne, L.E., Philippin, G.A., Vernier Piro, S.: Blow-up phenomena for a semilinear heat equation with nonlinear boundary condition, I. Z. Angew. Math. Phys. 61, 999–1007 (2010)

    Article  MathSciNet  Google Scholar 

  6. Payne, L.E., Philippin, G.A., Vernier Piro, S.: Blow-up phenomena for a semilinear heat equation with nonlinear boundary condition, II. Nonlinear Anal., Real World Appl. 73, 971–978 (2010)

    Article  MathSciNet  Google Scholar 

  7. Payne, L.E., Schaefer, P.W.: Bounds for blow-up time for the heat equation under nonlinear boundary conditions. Proc. R. Soc. Edinb., Sect. A 139, 1289–1296 (2009)

    Article  MathSciNet  Google Scholar 

  8. Ding, J.T., Shen, X.H.: Blow-up analysis in quasilinear reaction–diffusion problems with weighted nonlocal source. Comput. Math. Appl. 75, 1288–1301 (2018)

    Article  MathSciNet  Google Scholar 

  9. Jiang, L.J., Xu, Y.P.: Uniform blow-up rate for parabolic equations with a weighted nonlocal nonlinear source. J. Math. Anal. Appl. 365, 50–59 (2010)

    Article  MathSciNet  Google Scholar 

  10. Ling, Z., Lin, Z.G., Pedersen, M.: Global existence and blowup for a parabolic equation with a non-local source and absorption. Acta Appl. Math. 124, 171–186 (2013)

    Article  MathSciNet  Google Scholar 

  11. Baghaei, K., Hesaaraki, M.: Lower bounds for the blow-up time in the higher-dimensional nonlinear divergence form parabolic equations. C. R. Acad. Sci. Paris, Ser. I 351, 731–735 (2013)

    Article  MathSciNet  Google Scholar 

  12. Baghaei, K., Ghaemi, M.B., Hesaaraki, M.: Lower bounds for the blow-up time in a semilinear parabolic problem involving a variable source. Appl. Math. Lett. 27, 49–52 (2014)

    Article  MathSciNet  Google Scholar 

  13. Enache, C.: Lower bounds for blow-up time in some non-linear parabolic problems under Neumann boundary conditions. Glasg. Math. J. 53, 569–575 (2011)

    Article  MathSciNet  Google Scholar 

  14. Marras, M., Vernier Piro, S.: On global existence and bounds for blow-up time in nonlinear parabolic problems with time dependent coefficients. Discrete Contin. Dyn. Syst. 2013(Suppl.), 535–544 (2013)

    MathSciNet  MATH  Google Scholar 

  15. Marras, M., Vernier Piro, S., Viglialoro, G.: Blow-up phenomena for nonlinear pseudo-parabolic equations with gradient term. Discrete Contin. Dyn. Syst., Ser. B 22, 2291–2300 (2017)

    MathSciNet  MATH  Google Scholar 

  16. Zhang, L.L., Wang, H., Wang, X.Q.: Global and blow-up analysis for a class of nonlinear reaction diffusion model with Dirichlet boundary conditions. Math. Methods Appl. Sci. 17, 7789–7803 (2018)

    Article  MathSciNet  Google Scholar 

  17. Brezis, H.: Functional Analysis, Sobolev Spaces and Partial Differential Equations. Universitext. Springer, New York (2011)

    MATH  Google Scholar 

  18. Mizuguchi, M., Tanaka, K., Sekine, K., Oishi, S.: Estimation of Sobolev embedding constant on a domain dividable into bounded convex domains. J. Inequal. Appl. 2017, 299 (2017)

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

The authors are sincerely grateful to the editors and anonymous referees for careful reading of the original manuscript and useful comments.

Availability of data and materials

Not applicable.

Funding

This work was supported by the National Natural Science Foundation of China (No. 61473180).

Author information

Authors and Affiliations

Authors

Contributions

The authors contributed equally to the manuscript. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Juntang Ding.

Ethics declarations

Competing interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Additional information

Abbreviations

Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kou, W., Ding, J. Global existence and blow-up analysis for parabolic equations with nonlocal source and nonlinear boundary conditions. Bound Value Probl 2020, 37 (2020). https://doi.org/10.1186/s13661-020-01340-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13661-020-01340-5

MSC

Keywords