Blow-up time estimate for porous-medium problems with gradient terms under Robin boundary conditions

This paper deals with the blow-up phenomena connected to the following porous-medium problem with gradient terms under Robin boundary conditions: {ut=Δum+k1up−k2|∇u|qin Ω×(0,t∗),∂u∂ν+γu=0on ∂Ω×(0,t∗),u(x,0)=u0(x)≥0in Ω‾,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ \textstyle\begin{cases} u_{t}=\Delta u^{m}+k_{1}u^{p}-k_{2} \vert \nabla u \vert ^{q} & \text{in } \Omega\times(0,t^{*}), \\ \frac{\partial u}{\partial\nu}+\gamma u=0 &\text{on } \partial\Omega\times(0,t^{*}), \\ u(x,0)=u_{0}(x)\geq0 &\text{in } \overline{\Omega}, \end{cases} $$\end{document} where Ω⊂Rn\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\Omega\subset\mathbb{R}^{n}$\end{document} (n≥3\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$n\geq3$\end{document}) is a bounded and convex domain with smooth boundary ∂Ω. The constants p, q, m are positive, and p>q>m>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$p>q>m>1$\end{document}, q>2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$q>2$\end{document}. By making use of the Sobolev inequality and the differential inequality technique, we obtain a lower bound for the blow-up time of the solution. In addition, an example is given as an application of the abstract results obtained in this paper. Our results can be regarded as an answer to the open question raised by Li et al. in (Z. Angew. Math. Phys. 70:1–18, 2019).


Introduction
In recent years, blow-up phenomena for parabolic problems have been actively studied by many researchers. Such an issue is quite important, since it is often used to describe many physical and biological phenomena; we refer, for instance, to [2][3][4][5][6][7][8][9]. Most of these papers deal with the existence, blow-up, and other qualitative properties of solutions to related problems. In many situations, the techniques used in the study of blow-up phenomena often lead to the upper bound for the blow-up time when blow-up occurs. However, in applications, the lower bound seems to be more important, due to the explosive nature of the solution (see [10][11][12][13][14]). In this paper, we investigate the blow-up phenomena connected to the following porous-medium problem under Robin boundary conditions: where ⊂ R n (n ≥ 3) is a bounded and convex domain with smooth boundary ∂ , ∂u ∂ν is the normal derivative of u, t * is the maximal existence time of u, and is the closure of . Additionally, we set R + = (0, +∞) and assume that k 1 , k 2 are nonnegative constants, u 0 is a nonnegative C 1 ( ) function satisfying compatibility conditions.
There is a vast literature on the work of the blow-up phenomena for parabolic equations with the gradient terms. We mainly concentrate on the following papers [1,[15][16][17]. Marras, Vernier Piro and Viglialoro in [16] considered the following problem where is a bounded domain in ⊂ R n (n ≥ 2) whose boundary is sufficiently smooth. When ⊂ R 3 , the authors derived a lower bound for the blow-up time under suitable assumptions. Li, Pintus and Viglialoro [1] investigated the following problem where ⊂ R 3 is a bounded and smooth domain. When ⊂ R 3 , the authors obtained a lower bound for the blow-up time by using the differential inequality technique. Furthermore, we note that an open problem in [1] is how to give the lower bound of the blow-up time when the Dirichlet boundary conditions in (3) are replaced by Robin ones. Under this situation and inspired by the aforementioned works, we studied the blow-up phenomena of problem (1). The innovation of this paper is to establish suitable auxiliary functions. Since auxiliary functions defined in problems (2) and (3) are not suitable for our study, we need to construct new auxiliary functions. Moreover, instead of considering the lower bound for the blow-up time of the problem in ⊂ R 3 , we study the more general case when ⊂ R n (n ≥ 3). Combining the Sobolev inequality and a differential inequality technique, we give a lower bound for the blow-up time.
The rest of this paper is organized as follows. In Sect. 2, we derive a lower bound for the blow-up time when blow-up occurs. In Sect. 3, an example is given to illustrate the application of the results obtained in this paper.

Lower bound for blow-up time
In this section, we aim to obtain a lower bound for the blow-up time t * that gives a safe interval of the existence of the solution. To achieve this, we define L 0 = min ∂ (x · ν), d = max |x|, where x is a vector relative to the origin 0 and ν is a unit normal vector directed outward on ∂ . Suppose that k 1 , k 2 are nonnegative constants and m, p, q are positive constants with We now define the following auxiliary function Let λ 1 be the first eigenvalue of the plastic-membrane problem It follows from [18] that Hence, we have where Cs = Cs(n, ) is a Sobolev embedding constant depending on n and . In this section, we need to use the Sobolev inequality (7). Now, we state our main result as follows.
Theorem 2.1 Let u be a nonnegative classical solution of (1). Assume conditions (4)- (6) hold. If the solution u becomes unbounded in the measure (t) at some finite time t * , then Proof Using the divergence theorem and assumption (5), we derive First, we deal with the second term of (12). For this purpose, we consider the following fact that Using the Hölder inequality, we have Applying the general Poincaré inequality [19] to the term u α+q-1 dx in (13), we derive Combining the Lemma in [20], the Hölder inequality and the Young inequality, we obtain We insert (15) into (14) to deduce where M > 0 is given in (11) and assumption (6) guarantees M > 0. From (16) and (13), we derive Hence, multiplying both sides of the inequality (17) by ( q α+q-1 ) q , we obtain the relationship of the second term in (12) We substitute (18) into (12) to obtain Secondly, we pay attention to the third term of (19). It follows from the Hölder inequality and the Young inequality that where 0 < p-m q-m < 1 due to (4). The substitution of (20) into (19) leads to where ε 1 is defined in (9). For the term u α+m-1 dx of (21), using the Hölder inequality and the Sobolev inequality (7), we have It is easy to see that 0 < n(m-1) 2α+n(m-1) < 1. In view of the inequality we can rewrite (22) as . (23) To the first term of the right-hand side of (23), we make use of the Hölder inequality and the Young inequality to derive To the second term of the right-hand side of (23), it follows from the Hölder inequality and the Young inequality that where ε 2 is defined in (10). It follows from (23)-(25) that Thirdly, we insert (26) into (22) to obtain where C 1 is defined in (8). Integrating (27) from 0 to t, we deduce Since u blows up at finite time t * in measure (t), passing the limits as t → t * , we obtain where α+m-1 α > 1 in view of (4) and (5).

Application
In what follows, an example is given to illustrate the application of the abstract results.
Since we assume that u becomes unbounded at t * in measure (t), it follows from the Theorem 2.1 that u blows up in finite time and

Conclusion
In this paper, we make use of the Sobolev inequality (7) to overcome the difficulty mentioned in [1] when finding the lower bound for the blow-up time under Robin boundary conditions. The results of this paper provide an answer to the open question raised in [1].