1 Introduction
In this paper, we study the following 4th-order parabolic problem involving the Frobenius norm of a Hessian matrix:
where $\Omega \subset \mathbb {R}^N\left ( 1\leq N\leq 3\right )$ is a general bounded domain with smooth boundary, $\eta $ is the unit outward normal vector on $\partial \Omega $ , the initial datum $u_0\in H_0^2(\Omega )$ , the exponent p is a positive constant, T is the maximal existence time of (1.1), and the Frobenius norm of the Hessian matrix is defined as
By direct computation, problem (1.1) can be rewritten as
Problem (1.1) or (1.2) could be used to describe the growth of thin surfaces when exposed to molecular beam epitaxy (see [Reference Ortiz, Repetto and Si13, Reference Zangwill17, Reference Kim and Das Sarma9, Reference Das Sarma and Ghaisas5]). In particular, u can either represent the absolute thickness of the film or rather the relative surface height – that is, the deviation of the film height at the point x from the mean film thickness at time t (see [Reference Winkler14]); $-\Delta u$ indicates the diffusion due to evaporation-condensation (see [Reference Mullins12]); $\Delta ^{2} u$ indicates capillarity-driven surface diffusion (see [Reference Mullins12]); $\Delta |\nabla u|^{2}-\textrm {div}(2\Delta u\nabla u)$ is related to the equilibration of the inhomogeneous concentration of the diffusing particles on the surface (see [Reference Agélas1]); the source term $u^{p}$ denotes the mean deposition flux of the superlinear growth conditions with respect to u at $\infty $ , which could lead to the singularity of solutions or their derivatives to (1.1) or (1.2).
The parabolic equation in (1.1) or (1.2) is a typical equation of the continuum model of motion for the evolution of the film surface height $u(x, t)$ :
where f is the deposition flux and $\eta $ is the Gaussian random variable which describes the fluctuations in the average deposition flux.
In the case $A_1$ , $A_2> 0$ , $A_3 < 0$ and $A_4 = 0$ , Kohn and Yan in [Reference Kohn and Yan10] considered
where $\Omega \subset \mathbb {R}^2$ is a square domain. They obtained the decay of energy in time.
In the case $A_1$ , $A_3<0$ , $A_2> 0$ , and $A_4 = 0$ , Liu and Li in [Reference Liu and Li11] studied
where $\Omega \subset \mathbb {R}^{N}( N\geq 1)$ is a square domain. They obtain the sufficient conditions on the global existence, asymptotic behavior, and finite time blow-up of weak solutions, but also show exact descriptions of smallness conditions on the initial data.
In the case $A_1$ , $A_2$ , $A_4> 0$ , and $A_3 = 0$ , Winkler in [Reference Winkler14] investigated the following equation by using the computational methods
where $\Omega \subset \mathbb {R}^{N}(1\leq N\leq 3)$ is a bounded convex domain with smooth boundary, $\mu \geq 0$ and $\lambda> 0$ . Under appropriate assumptions on f, the global existence of weak solutions was obtained. Under an additional smallness condition on $\mu $ and the size of f, it was shown that there exists a bounded set which is absorbing for (1.3) in some sense for any solution. Blomker and Gugg in [Reference Blömker and Gugg3] also studied the related problem
where $\Omega \subset \mathbb {R}$ is a bounded interval. The global existence of weak solutions was proved. This result was extended by Blomker et al. in [Reference Blömker, Gugg and Raible4] to the parabolic equation
where $\Omega \subset \mathbb {R}$ is a bounded interval and $\nu> 0$ .
In the case $A_1$ , $A_2$ , $A_4> 0$ and $A_3<0$ , Agélas in [Reference Agélas1] dealt with
where $\Omega = \mathbb {R}^N (N = 1, 2)$ . He proved the existence, uniqueness, and regularity of global weak solutions. Moreover, under the condition $\nu _2\nu _3> \nu _4^{2}$ , the author proved the existence and uniqueness of global strong solutions for sufficiently smooth initial data.
In [Reference Escudero6], Escudero dealt with both the initial and initial-boundary value problems for the partial differential equation $u_t + \Delta ^2u =\textrm {det}(D^{2}u)$ posed either on $\mathbb {R}^2$ or on a bounded subset of the plane, where $\textrm { det}(D^{2}u)$ is the determinant of the Hessian matrix $D^{2}u$ . The author studied the blow-up behavior including the complete blow-up in either finite or infinite time. Moreover, he refined a blow-up criterium that was proved for this evolution equation. The interested authors could find other results in [Reference Xu and Zhou16, Reference Zhou18] and the papers cited therein.
To our knowledge, the 4th-order parabolic problem (1.1) involving a Frobenius type nonlinearity has been rarely considered before. Moreover, the mean deposition flux of the superlinear growth conditions would play an important role in the property of the solutions, including the existence of blow-up, extinction solutions. Inspired by the works [Reference Winkler14, Reference Blömker and Gugg3], we want to study the threshold results on the initial data with respect to the existence of blow-up, global, and extinction solutions of (1.1) or (1.2). Throughout this paper, we denote by $\|\cdot \|_p$ the $L^p(\Omega )$ norm and by $(\cdot ,\cdot )$ the inner product in $L^{2}(\Omega )$ , respectively. For $u\in H_0^{2}(\Omega )$ with norm $\|u\|_{H_0^{2}(\Omega )}=\|\Delta u\|_2$ , we define the energy functional and the Nehari functional, respectively,
which satisfy
We give the weak solutions of problem (1.1) as follows.
Definition 1.1 Let $T>0$ . A function $u(x,t) \in L^{\infty }(0,T;H_0^{2}(\Omega ))$ with $u_t\in L^{2}(0,T;L^{2}(\Omega ))$ is the so-called weak solution to (1.1) or (1.2) in $\Omega \times [0,T)$ , if $u(x,0)=u_0 (x )\in H_0^{2}(\Omega )$ , for any $\varphi (x ) \in H_0^{2}(\Omega )$ ,
Moreover, there is the relationship for the energy of the weak solutions,
Define the Nehari manifold $ \mathcal {N}:=\left \{u\in H_0^{2}(\Omega ) |\ I(u)=0,\ \|\Delta u\|_2\neq 0\right \}$ . The potential well and its corresponding sets are defined by
where $ d:=\operatorname *{inf}\limits _{u \in H_0^{2}(\Omega )\backslash \{0\}}\operatorname *{sup}\limits _{\lambda \geq 0}J(\lambda u)=\operatorname *{inf}\limits _{ u \in \mathcal {N}}J(u)$ is the so-called depth of the potential well $\mathcal {W}$ .
For any $\delta>0$ , we further define the modified functional and the Nehari manifold as
The modified potential wells and their corresponding sets are defined respectively by
Here, $d(\delta ):=\inf _{ u \in \mathcal {N}_\delta }J(u)>0$ is the potential depth of $\mathcal {W}_\delta $ . We also define the open sublevels of J, $ J^{s}:=\left \{u\in H_0^{2}(\Omega )|\ J(u)<s\right \}$ . Furthermore, by the definitions of $J(u)$ , $\mathcal {N}$ and $J^s$ , we see that $ \mathcal {N}^{s}:=\mathcal {N}\cap J^{s}\neq \emptyset $ for $\forall s>d$ . For any $s>d$ , we define
It is clear that $\lambda _s$ is nonincreasing and $\Lambda _s$ is nondecreasing with respect to s, respectively.
We summarize the main results through the following table. The abbreviations “N.E.,” “E.,” “B.U.,” and “G.E.” denote non-extinction, extinction, blow-up, and global existence of weak solutions of (1.1) or (1.2), respectively. It could be checked that if $J(u_0)<d$ , then $I(u_0)\not =0$ .
This paper is arranged as follows. In the next section, we give some important lemmas. Sections 3, 4, and 5 are devoted to the subcritical, the critical, and the supercritical energy cases, respectively. Section 6 gives the upper and the lower bounds of blow-up time of weak solutions. In Section 7, we show some results about non-extinction or extinction of weak solutions.
2 Preliminary Lemmas
In this section, we give ten lemmas which play important roles in the proof of the main results.
Lemma 2.1 Let $p>2$ . For any $u\in H_0^{2}(\Omega )$ with $\|\Delta u\|_2\neq 0$ , we have
-
(i) $\operatorname *{lim}\limits _{\lambda \rightarrow 0^{+}}J(\lambda u)=0$ , $\operatorname *{lim}\limits _{\lambda \rightarrow +\infty }J(\lambda u)=-\infty $ .
-
(ii) There exists an unique constant $\lambda ^{*}=\lambda ^{*}(u)>0$ such that $\frac {\textrm {d}}{\textrm {d}\lambda }J(\lambda u)|_{\lambda =\lambda ^{*}}=0$ . $J(\lambda u)$ is increasing for $0<\lambda <\lambda ^{*}$ , is decreasing for $\lambda ^{*}<\lambda <+\infty $ , and takes its maximum at $\lambda =\lambda ^{*}$ .
-
(iii) $I(\lambda u)>0$ for $0<\lambda < \lambda ^{*}$ , $I(\lambda u)<0$ for $\lambda ^{*}<\lambda <+\infty $ , and $I(\lambda ^{*} u)=0$ .
Proof (i) Define the function $j:\lambda \mapsto J(\lambda u)$ for $\lambda> 0$ . Then
We obtain $\operatorname *{\lim }\limits _{\lambda \rightarrow 0^{+}} J(\lambda u)=0$ and $\operatorname *{\lim }\limits _{\lambda \rightarrow +\infty } J(\lambda u)=-\infty $ .
(ii) Elementary calculations imply that
Let $k(\lambda ):=\lambda ^{-2} j'(\lambda )$ . After direct calculation, we have
Since $\operatorname *{\lim }\limits _{\lambda \rightarrow 0^{+}} k(\lambda )=+\infty $ , $\operatorname *{\lim }\limits _{\lambda \rightarrow +\infty } k(\lambda )=-\infty $ , there exists a unique constant $\lambda ^{*}>0$ such that $k(\lambda )>0$ for $0<\lambda <\lambda ^{*}$ , $k(\lambda )<0$ for $\lambda ^{*}<\lambda <+\infty $ , and $k(\lambda ^{*})=0$ . By $j'(\lambda )=\lambda ^{2}k(\lambda )$ , $I(\lambda u)=\lambda j'(\lambda )$ , cases (ii) and (iii) hold.
Lemma 2.2 Let $p>2$ . The depth d of the potential well $\mathcal {W}$ is positive.
Proof Employing H $\ddot {\textrm {o}}$ lder’s inequality, we obtain
Fix $u\in \mathcal {N}$ . Since $p>2$ and by (1.4),
where $B_1$ is the optimal constant in the embedding $H_0^{2}(\Omega )\hookrightarrow L^{p+1}(\Omega )$ , and $B_2$ is the optimal constant in the embedding $W_0^{1,2}(\Omega )\hookrightarrow L^{4}(\Omega )$ . Let
Then $ J(u) =\frac {1}{6}\|\nabla u\|_2^{2}+\frac {1}{6}\|\Delta u\|_2^{2}+\frac {p-2}{3(p+1)}\| u\|_{p+1}^{p+1}+\frac {1}{3}I(u) \geq \frac {1}{6}\|\Delta u\|_2^{2}\geq \frac {1}{6}B_{0}^{2}>0$ . Therefore, $d=\operatorname *{\inf }\limits _{u\in \mathcal {N}}J(u)>0$ .
Lemma 2.3 Suppose $p>2$ and $u\in H_0^{2}(\Omega )$ . Define $ r(\delta ):=\inf \{x\in (0,+\infty )|\ \delta \leq B_1^{p+1}x^{p-1}+3B_2^{2}x \}$ , where $B_1$ , $B_2$ are defined in (2.1). These are the following results.
-
(i) If $I_{\delta }(u)<0$ , then $\|\Delta u\|_{2}>r(\delta )$ . Specially, if $I(u)<0$ , then $\|\Delta u\|_{2}>r(1)$ .
-
(ii) If $0\leq \|\Delta u\|_{2}\leq r(\delta )$ , then $I_{\delta }(u)\geq 0$ . Specially, if $0\leq \|\Delta u\|_{2}\leq r(1)$ , then $I(u)\geq 0$ .
-
(iii) If $I_{\delta }(u)=0$ , then $\|\Delta u\|_{2}=0$ or $\|\Delta u\|_{2}\geq r(\delta )$ . Specially, if $I(u)=0$ , then $\|\Delta u\|_{2}=0$ or $\|\Delta u\|_{2}\geq r(1)$ .
Proof If $I_{\delta }(u)<0$ and $ \delta \|\Delta u\|_2^{2}<\delta \|\Delta u\|_2^{2}+\|\nabla u\|_{2}^{2}\leq B_1^{p+1}\|\Delta u\|_{2}^{p+1}+3B_2^{2}\|\Delta u\|_{2}^{3}$ , then $\|\Delta u\|_2>r(\delta )$ , and hence, cases (i) and (ii) hold. If $I_{\delta }(u)=0$ , we get $\|\Delta u\|_2\geq r(\delta )$ . If $\|\Delta u\|_2=0$ , $I_{\delta }(u)=0$ .
Lemma 2.4 (Lemma 2.4 in [Reference Han7])
$d(\delta )$ is increasing for $0<\delta \leq 1$ , is decreasing for $\delta \geq 1$ , and takes its maximum $d=d(1)$ at $\delta =1$ .
Lemma 2.5 (Lemma 5 in [Reference Xu and Su15])
Let $p>2$ . Assume $u\in H_0^{2}(\Omega )$ , $0<J(u)<d$ , and ${\delta _1<1<\delta _2}$ is the two roots of the equation $d(\delta )=J(u)$ . Then the sign of $I_{\delta }(u)$ does not change for $\delta _1<\delta <\delta _2$ .
Lemma 2.6 (Lemma 8 in [Reference Xu and Su15])
Let $p>2$ and assume that $u $ is a weak solution of problem (1.1) in $\Omega \times [0,T)$ with $0<J(u_{0})<d$ . Let $\delta _1<1<\delta _2$ be the two roots of the equation $d(\delta )=J(u_{0})$ .
-
(i) If $I(u_{0})>0$ , then $u \in \mathcal {W}_{\delta }$ for $\delta _{1}<\delta <\delta _{2}$ and $0<t<T$ .
-
(ii) If $I(u_{0})<0$ , then $u \in \mathcal {V}_{\delta }$ for $\delta _{1}<\delta <\delta _{2}$ and $0<t<T$ .
Lemma 2.7 Let $p>2$ . $\textrm {dist}(0,\mathcal {N})>0$ and $\textrm {dist}(0,\mathcal {N}_{-})>0$ .
Proof For any $u\in \mathcal {N}$ , by the definition of d, we obtain
which indicates that $ \frac {1}{6}B_{3}^{2}\|\Delta u\|_2^{2}+\frac {1}{6}\|\Delta u\|_2^{2}\geq \frac {d}{2}$ or $\frac {p-2}{3(p+1)}B_{1}^{p+1}\|\Delta u\|_{2}^{p+1}\geq \frac {d}{2}$ . Then $\|\Delta u\|_2\geq \left (\frac {3d}{B_{3}^{2}+1}\right )^{\frac {1}{2}}$ , or $\|\Delta u\|_{2}\geq \left (\frac {3d(p+1)}{2(p-2) B_{1}^{p+1}}\right )^{\frac {1}{p+1}}$ , where $B_{3}$ is the optimal constant in $W_0^{1,2}(\Omega )\hookrightarrow L^{2}(\Omega )$ . Let
Then $\textrm {dist}(0,\mathcal {N})=\operatorname *{\inf }\limits _{u\in \mathcal {N}}\|\Delta u\|_2\geq C_{0}>0$ . For any $u\in \mathcal {N}_{-}$ , we have $ \|\Delta u\|_2^{2}\leq B_1^{p+1}\|\Delta u\|_{2}^{p+1}+3B_2^{2}\|\Delta u\|_{2}^{3}$ , which implies $ \|\Delta u\|_2\geq B_{0}$ . Here, $B_{0}$ is given in (2.1). Then $\textrm {dist}(0,\mathcal {N}_{-})=\operatorname *{\inf }\limits _{u\in \mathcal {N}^{-}}\|\Delta u\|_2\geq B_{0}>0$ .
Lemma 2.8 Let $p>2$ . For any $s>d$ , $u\in J^{s}\cap \mathcal {N}_{+}$ , $ \|\Delta u\|_2<C_{3}:=(6s)^{\frac {1}{2}}$ .
Proof For any $s>d$ and $p>2$ , $u\in J^{s}\cap \mathcal {N}_{+}$ , we have
which yields $ \|\Delta u\|_2<C_{3}$ .
For suitable u and $p>2$ , by using the Gagliardo-Nirenberg inequality, we have
where $a:=\frac {N+4}{8}\in (0,1)$ and $b:=\left (\frac {1}{2}-\frac {1}{p+1}\right )\cdot \frac {N}{2}\in (0,1)$ .
Lemma 2.9 Let $p>2$ . For any $s>d$ , $\lambda _{s}$ and $\Lambda _{s}$ in (1.7) satisfy $0<K_{4}\leq \lambda _{s}\leq \Lambda _{s}\leq K_{1}<+\infty $ , where $C_0$ is defined in (2.2) and the constants a, b, $C_1$ , $C_2$ are defined in (2.3); $K_{1}:=B_{4}C_{3}$ ,
Proof For $u\in \mathcal {N}^{s}$ , we have $ \frac {1}{6}B_{4}^{-2}\|u\|_2^{2}\leq \frac {1}{6}\|\Delta u\|_2^{2}\leq J(u)<s$ , where $B_{4}$ is the optimal constant in the embedding $H_0^{2}(\Omega )\hookrightarrow L^{2}(\Omega )$ . Then $ \|u\|_{2}\leq K_{1}$ . By (2.3), we have
and hence,
By the similar proof of Lemma 2.8, we get $\|\Delta u\|_{2}<C_{3}$ . Combining with (2.4), we have $ \| u\|_{2}\geq K_{2}$ or $\| u\|_{2}\geq K_{3}$ . Then $\|u\|_{2}\geq K_{4}>0$ ; hence, $\lambda _{s}>0$ .
Lemma 2.10 (Lemma 2.1 in [Reference Khelghati and Baghaei8])
Suppose that a positive, twice-differentiable function $\theta (t)$ satisfies the inequality $ \theta "(t)\theta (t)-(1+\beta )\theta '(t)^{2}\geq 0$ , $t>0$ , where $\beta>0$ is a constant. If $\theta (0)>0$ and $\theta '(0)>0$ , then there exists $0<t_{1}<\frac {\theta (0)}{\beta \theta '(0)}$ such that $\theta (t)$ tends to $\infty $ as $t\rightarrow t_{1}$ .
3 The subcritical case
This section is devoted to the property of weak solution of (1.1) or (1.2) under the case $J(u_0)<d$ .
Theorem 3.1 Let $p>2$ . If $J(u_0)<d$ and $I(u_0)>0$ , then problem (1.1) admits a global weak solution $u\in L^{\infty }(0,\infty ;H_{0}^{2}(\Omega ))$ with $u_{t}\in L^{2}(0,\infty ;L^{2}(\Omega ))$ and $u(t)\in \mathcal {W}$ for ${0 \leq t <\infty }$ . Moreover, there exists a constant $\hat C> 0$ such that $\|u\|_{2}^{2}\leq \|u_{0}\|_{2}^{2}\textrm {e}^{-\hat Ct}$ .
Proof The proof is divided into two steps.
Step 1. Global existence. We would use the Galerkin’s approximation with some priori estimates. Let $\left \{\omega _{i}(x)\right \}$ be the orthogonal basis of $H_{0}^{2}(\Omega )$ . Construct the approximate solutions $u_{m}(x,t)$ of (1.1), $ u_{m}(x,t):=\operatorname *{\sum }\limits _{i=1}^{m} a_{mi}(t)\omega _{i}(x)$ , $m=1,2,\cdots $ , $i=1,2,\cdots ,m$ , which satisfy
and $u_{0m}:=\operatorname *{\sum }\limits _{i=1}^{m} b_{mi}(t)\omega _{i}(x)\rightarrow u_{0}(x)$ in $ H_{0}^{2}(\Omega )$ as $m\rightarrow +\infty $ .
By the standard theory of ODEs (e.g., the Peanos theorem), we deduce that the existence of a local solution to (3.1). Multiplying (3.1) by $a_{mi}'(t)$ , summing over i from 1 to m and integrating with respect to t, we have
Due to the convergence of $u_{0m}\rightarrow u_{0}(x)$ in $H_{0}^{2}(\Omega )$ , one has $ J(u_{m}(x,0))\rightarrow J(u_{0}(x))<d$ , $I(u_{m}(x,0))\rightarrow I(u_{0}(x))>0$ . Therefore, for sufficiently large m and any $0\leq t <+\infty $ , we obtain
which implies that $u_{m}(x,0)\in \mathcal {W}$ for sufficiently large m.
By applying the similar discussion of Theorem 8 in [Reference Xu and Su15], one could show from (3.3) that $u_{m}(x,t)\in \mathcal {W}$ for large m and $0\leq t <+\infty $ . Thus, $I(u_{m}(x,t))>0$ , $J(u_{m}(x,t))<d$ for all $t\in [0,T]$ . Then
In addition, by using (3.2–3.4), we get for some positive constant C,
By the uniform estimates (3.5–3.8), it was seen that the local solutions can be extended globally. Thus, by the standard diagonal method and the Aubin-Lions-Simon theorem, we know there exists a function u and a sequence of $\left \{u_{m}\right \}$ (still by $\left \{u_{m}\right \}$ ) such that for each $T>0$ , one could obtain
as $m\rightarrow +\infty $ . Fix $k\in \mathcal {N}$ . In order to show the limit function u in (3.9–3.14) is a weak solution of (1.1), we choose a function $v\in C^1 ([0,T];H_{0}^{2}(\Omega ))$ defined as $ v(x,t):=\operatorname *{\sum }\limits _{i=1}^{k} l_{i}(t)\omega _{i}(x)$ , where $\left \{l_{i}(t)\right \}_{i=1}^{k}$ are arbitrarily given $C^1$ functions. Taking $m\geq k$ in (3.1), multiplying (3.1) by $l_{i}(t)$ , summing for i from 1 to k, and integrating with respect to t, we obtain
Taking $m\rightarrow +\infty $ in (3.15) and recalling the convergence yield that
Since functions of the form in (3.15) are dense in $L^{2}(0,T;H_{0}^{2}(\Omega ))$ , (3.16) also holds for all $v\in L^{2}(0,T;H_{0}^{2}(\Omega ))$ . By the arbitrariness of $T>0$ , we know
Then u is a global weak solution to problem (1.1). To prove (1.6), we first assume that u was smooth enough that $u_{t}\in L^{2}(0,T;H_{0}^{2}(\Omega ))$ . Taking $v=u_{t}$ in (3.16), it is seen that (1.6) is true. By the density of $L^{2}(0,T;H_{0}^{2}(\Omega ))$ in $L^{2}(\Omega \times (0,T))$ , we know (1.6) also holds for weak solutions to (1.1). The existence of global solutions to (1.1) is obtained.
Step 2. Decay rate. Taking $\varphi :=u$ in (1.5), we get $\frac {\textrm {d}}{\textrm {d}t} {\|u\|_2^{2}}=2(u_{t},u)=-2I(u)$ . From Lemma 2.6, we know that $u(x,t)\in \mathcal {W}_{\delta }$ for $\delta _{1} < \delta <\delta _{2}$ and $0 < t < \infty $ under the condition $J(u_0) < d$ and $I(u_0)> 0$ . Thus, $I_{\delta _{1}}(u)\geq 0$ for $0 < t < \infty $ . Therefore,
where $B_{4}>0$ is the best embedding constant from $H_{0}^{2}(\Omega )$ to $L^{2}(\Omega )$ (i.e., $\|u\|_{2}\leq B_{4}\|\Delta u\|_{2}$ for $\forall u \in H_{0}^{2}(\Omega )$ ). Consequently, $ \|u\|_2^{2}\leq \|u_{0}\|_2^{2}\textrm {e}^{-\hat C t}$ with $\hat {C}:=2B_{4}^{-2}(1-\delta _{1})>0$ .
Theorem 3.2 Let $p>2$ and u be a weak solution of (1.2). If $J(u_0)<d$ , $J(u_0)\not =0$ and $I(u_0)<0$ , then u blows up at some finite time T in the sense of $\operatorname *{\lim }\limits _{t\rightarrow T} \int _{0}^{t}\|u\|_{2}^{2}\textrm {d}\tau =+\infty $ .
Proof We employ the concavity method. Assume on the contrary that u was a global weak solution to (1.2) with $J(u_0) < d$ , $J(u_0)\not =0$ , $I(u_0) < 0$ and define $ F(t):= \int _{0}^{t}\|u\|_{2}^{2}\textrm {d}\tau $ , $t\geq 0$ . Then $F'(t) =\|u\|_{2}^{2}$ ,
By (1.6), (1.4), and (3.17), we have
Noticing that $ \frac {1}{4}(F'(t)-F'(0))^{2}= ( \int _{0}^{t} \int _{\Omega }uu'\textrm {d}x\textrm {d}\tau )^{2}\leq \int _{0}^{t}\|u\|_{2}^{2}\textrm {d}\tau \int _{0}^{t}\|u'\|_{2}^{2}\textrm {d}\tau $ , we have
In the forthcoming proof, the cases that $J(u_0)< 0$ and $ 0< J(u_{0}) < d$ will be discussed separately.
(i) If $J(u_0)< 0$ , then (3.18) implies $ F"(t)F(t)-\frac {3}{2}\left (F'(t)\right )^{2} \geq B_{4}^{-2}F(t)F'(t)-3F'(0)F'(t)$ . Now we will prove that $I(u) < 0$ for all $t> 0$ . Otherwise, there must be a constant $t_{0}> 0$ such that $I(u(t_{0})) = 0$ and $I(u) < 0$ for $0 \leq t < t_{0}$ . From Lemma 2.3 (i, iii), $\|\Delta u\|_{2}> r(1)$ for $0 \leq t < t_{0}$ , and $\|\Delta u(t_{0})\|_{2} \geq r(1)$ , where $\|\Delta u(t_{0})\|_{2} \not = 0$ . In fact, using (1.6) and $J(u_0)< 0$ , we have $J(u(t_0))< 0$ . Since $0=F"(t_0)\geq -6J(u(t_0))+\|\Delta u(t_0)\|_{2}^{2}$ , we have $\|\Delta u(t_{0})\|_{2} \not = 0$ . Hence, we have $J(u(t_{0})) \geq d$ , which contradicts to (1.6). By (3.17), we get $F"(t)>0$ for $t\geq 0$ . Since $F'(0) =\|u_{0}\|_{2}^{2}\geq 0$ , $F'(t)>0$ , for large t, $F(t)>3B_{4}^{2}\|u_{0}\|_{2}^{2}$ , we have
(ii) If $0<J(u_{0}) < d$ , then by Lemma 2.6 (ii), we know that $u(x,t)\in \mathcal {V}_{\delta }$ for $t\geq 0$ and $\delta _1 < \delta < \delta _2$ , where $\delta _1<1<\delta _2$ are the two roots of $d(\delta ) = J(u_0)$ . Hence, ${I_{\delta _2}(u)\leq 0}$ and $\|\Delta u\|_{2}\geq r(\delta _2)$ for $t\geq 0$ . It follows from (3.18) that for $t\geq 0$ , $ F"(t) \geq 2(\delta _2-1)r^{2}(\delta _2)$ , which shows for all $t \geq 0$ that $ F'(t)\geq 2(\delta _2-1)r^{2}(\delta _2)t$ and $ F(t)\geq (\delta _2-1)r^{2}(\delta _2)t^{2}$ . Therefore, for sufficiently large t, we have
Then we obtain (3.19). Due to $ (F^{-\frac {1}{2}}(t) )"=-\frac {1}{2}F^{-\frac {5}{2}}(t) [F(t)F"(t)-\frac {3}{2}\left (F'(t)\right )^{2} ]<0$ , $F^{-\frac {1}{2}}(t)$ is concave in $(0,+\infty )$ . So there exists a positive constant $t_{0}$ such that $F(t_0)>0$ , $F'(t_{0})>0$ and $\left (F^{-\frac {1}{2}}(t_{0})\right )'<0$ . Since $F(t)>0,F"(t)>0$ for $t\ge t_0$ and $F'(t_0)>0$ , one can find $(F^{-1/2})'(t)<0$ for $t\ge t_0$ , and hence, there is a constant $T>t_{0}$ such that $\operatorname *{\lim }\limits _{t\rightarrow T}F^{-\frac {1}{2}}(t)=0$ ; that is, $ \operatorname *{\lim }\limits _{t\rightarrow T} \int _{0}^{t}\|u(x,\tau )\|_{2}^{2}\textrm {d}\tau =+\infty $ .
4 The critical case
For $J(u_0)=d$ , the invariance of $\mathcal {W}_{\delta }$ could not be proved in general. By using approximation, we could prove the global existence of weak solutions.
Theorem 4.1 Let $p>2$ . If $J(u_0)=d$ and $I(u_0)\geq 0$ , then problem (1.2) admits a global weak solution $u\in L^{\infty }(0,\infty ;H_{0}^{2}(\Omega ))$ with $u_t\in L^{2}(0,\infty ;L^{2}(\Omega ))$ and $u \in \overline {\mathcal {W}}=\mathcal {W}\cup \partial \mathcal {W}$ for $0 \leq t < \infty $ . Moreover, if $I(u(x, t))> 0$ for all $t>0$ , $u(x,t)$ does not vanish and there exist positive constants $\hat C_{1}$ and $\hat C_{2}$ such that $\|u\|_2^{2}\leq \hat C_{1}\textrm {e}^{-\hat C_{2}t}$ . If not, problem (1.2) admits a solution that vanishes in finite time.
Proof Let $\lambda _{k}=1-\frac {1}{k}$ , $k=2,3,\cdots $ . Consider the approximation problems
Noticing that $I(u_0) \geq 0$ , by Lemma 2.1 (iii), we could deduce that there exists a unique constant $\lambda ^{*}=\lambda ^{*}(u_{0})\geq 1$ such that $I(\lambda ^{*}u_{0})=0$ . By $\lambda _{k}<1\leq \lambda ^{*}$ , we get $I(u_{0}^{k}) = I(\lambda _{k}u_{0})> 0$ and $J(u_{0}^{k})= J(\lambda _{k}u_{0}) <J(u_{0}) = d$ . In view of Theorem 3.1, for each k, problem (4.1) admits a global weak solution $u^{k}(x,t)\in L^{\infty }(0,\infty ;H_{0}^{2}(\Omega ))\cap \mathcal {W}$ with $u_{t}^{k}\in L^{2}(0,\infty ;L^{2}(\Omega ))$ for $0 \leq t < \infty $ satisfying $ \int _{0}^{t}\|u_{\tau }^{k}\|_{2}^{2}\textrm {d}\tau +J(u^{k})=J(u_{0}^{k})<d$ . Applying the similar discussion in Theorem 3.1, there exist a subsequence of $ \{u^{k} \}$ and a function u such that u is a weak solution of (1.2) with $I(u) \geq 0$ and $J(u) \leq d$ for $0 \leq t < \infty $ .
Let us discuss the decay rate of $\|u\|_{2}^{2}$ . First, suppose that $I(u)> 0$ for $0 < t < \infty $ ; then $u $ does not vanish in finite time. Combining with (3.17), we have $u_{t} \not \equiv 0$ . By (1.6), for small $t_0> 0$ , we have $ 0<J(u(t_{0}))<J(u_{0})=d$ . Taking $t = t_0$ as the initial time and by Lemma 2.6, we get $u\in \mathcal {W}_{\delta }$ for $\delta _1 < \delta < \delta _2$ and $t_0<t<\infty $ . Hence, $I_{\delta _1}(u) \geq 0$ for $t_0<t<\infty $ and
which implies $ \|u\|_{2}^{2}\leq \|u(t_{0})\|_{2}^{2}\textrm {e}^{-2B_{4}^{-2}(1-\delta _1)(t-t_{0})}$ . The decay rate holds with $\hat C_{2} := 2B_{4}^{-2}(1-\delta _1)$ and $\hat C_{1}:=\|u(t_{0})\|_2^{2}\textrm {e}^{2 B_{4}^{-2}(1-\delta _{1})t_{0}}$ .
Next, suppose there is a positive constant $t_1$ such that $I(u)> 0$ for $0 < t < t_1$ and $I(u(x, t_1)) = 0$ . Obviously, $u_t \not \equiv 0$ for $0 < t < t_1$ and $ \int _{0}^{t_{1}}\|u_{\tau }\|_{2}^{2}\textrm {d}\tau>0$ . Applying (1.6) again, we have $ J(u(t_{1}))=J(u_{0})- \int _{0}^{t_{1}}\|u_{\tau }\|_{2}^{2}\textrm {d}\tau <J(u_{0})=d$ . By the definition of d, we know $\|u(t_1)\|_{2} = 0$ and $\|u\|_{2} \equiv 0$ for all $t \geq t_1$ . It is seen that such a function $u $ is a weak solution of (1.2) which vanishes in finite time.
Theorem 4.2 Let $p>2$ . If $J(u_0) = d$ and $I(u_0) < 0$ , then there exists a finite time T such that the solution $u $ blows up at that time T in the sense of $\operatorname *{\lim }\limits _{t\rightarrow T} \int _{0}^{t}\|u\|_{2}^{2}\textrm {d}\tau =+\infty $ .
Proof Similarly to the proof of Theorem 4.1, we could get
Since $J(u_0) = d$ , $I(u_0) < 0$ , by the continuity of $J(u)$ and $I(u)$ with respect to t, there exists a constant $t_0> 0$ such that $J(u(x, t))> 0$ and $I(u(x, t)) < 0$ for $0 < t \leq t_0$ . From $(u_t, u) =-I(u)$ , we have $u_t \not \equiv 0$ for $0 < t \leq t_0$ . Furthermore, we have $ J(u(t_{0}))=d- \int _{0}^{t_{0}}\|u_{\tau }\|_{2}^{2}\textrm {d}\tau <d$ . Taking $t = t_0$ as the initial time and by Lemma 2.6, we know that $u(x, t)\in \mathcal {V}_{\delta }$ for $\delta _1 < \delta < \delta _2$ and $t> t_0$ , where $\delta _1 < 1 < \delta _2$ are the two roots of the equation $d(\delta ) = J(u_{0})$ . Thus, $I_{\delta _2}(u)\leq 0$ for $t> t_0$ . The rest of the proof is similar to that of Theorem 3.2.
5 The supercritical case
We obtain some results for arbitrarily high initial energy.
Theorem 5.1 Let $p>2$ and $J(u_0)> d$ .
-
(i) If $u_{0}\in \mathcal {N}_{+}$ and $\|u_{0}\|_{2}\leq \lambda _{J(u_{0})}$ , then the weak solution u of (1.1) exists globally and $u\rightarrow 0$ as $t\rightarrow +\infty $ .
-
(ii) If $u_{0}\in \mathcal {N}_{-}$ and $\|u_{0}\|_{2}\geq \Lambda _{J(u_{0})}$ , then the weak solution u of (1.1) blows up in finite time.
Proof Let $T_{\max }$ be the maximal existence time of (1.1). If $T_{\max } =\infty $ , we denote the $\omega $ -limit set of $u_0$ as $ \omega (u_{0})=\operatorname *{\bigcap }\limits _{t\geq 0}\overline {\left \{u(\tau ):\tau \geq t\right \}}^{H_{0}^{2}(\Omega )}$ .
(i) Assume that $u_0\in \mathcal {N}_{+}$ with $\|u_{0}\|_{2}\leq \lambda _J(u_{0})$ . We first claim that $u(t)\in \mathcal {N}_{+}$ for all $t\in [0, T_{\max })$ . If not, there would exist a constant $t_0\in (0, T_{\max })$ such that $u(t)\in \mathcal {N}_{+}$ for $0 \leq t < t_{0}$ and $I(u(t_0))=0$ . It follows from $I(u(t))=- \int _{\Omega }u_tudx$ that $u_t\not \equiv 0$ for $(x,t)\in \Omega \times [0, t_0)$ . Using (1.6), we have $J(u(t_0)) < J(u_0)$ , which deduces $u(t_0)\in J^{J(u_0)}$ . Therefore, $u(t_0)\in J^{J(u_{0})}\cap \mathcal {N}=\mathcal {N}^{J(u_0)}$ . By the definition of $\lambda _{J(u_0)}$ , we have
Noticing $I(u(t))> 0$ for $t\in [0, t_0)$ , $ \|u(t_{0})\|_{2}<\|u_{0}\|_{2}\leq \lambda _{J(u_{0})}$ , a contradiction to (5.1). So $u(t)\in \mathcal {N}_{+}$ and $u(t)\in J^{J(u_0)}\cap \mathcal {N}_{+}$ for all $t\in [0, T_{\max })$ . Lemma 2.8 shows $ \|\Delta u\|_{2}<C_{3}$ , $t\in [0,T_{\max })$ , so that $T_{\max } =\infty $ . Let $\omega $ be an arbitrary element in $\omega (u_{0})$ . We have $ J(\omega )<J(u_{0})$ , $\|\omega \|_{2}<\lambda _{J(u_{0})}$ , which implies $\omega \notin \mathcal {N}^{J(u_{0})}$ . Recalling the definition of $\lambda _{J(u_0)}$ , $\omega (u_0)\cap \mathcal {N} = \emptyset $ . Hence, $\omega (u_0)=\left \{0\right \}$ . Therefore, the weak solution u of problem (1.1) exists globally and $u\rightarrow 0$ as $t\rightarrow +\infty $ .
(ii) Assume that $u_0\in \mathcal {N}_{-}$ with $\|u_{0}\|_{2}\geq \Lambda _{J(u_{0})}$ . We claim that $u(t)\in \mathcal {N}_{-}$ for all $t\in [0, T_{\max })$ . If not, there would be a constant $t_{0}\in (0, T_{\max })$ such that $u(t)\in \mathcal {N}_{-}$ for $0 \leq t < t_0$ and $I(u(t_{0}))=0$ . Similarly to case (i), one has $J(u(t_0)) < J(u_0)$ , which implies $u(t_0)\in J^{J(u_0)}$ , and $u(t_0)\in \mathcal {N}^{J(u_0)}$ . By the definition of $\Lambda _{J(u_0)}$ , we have
However, $I(u(t)) < 0$ for $t\in [0, t_0)$ , and we get $ \|u(t_{0})\|_{2}>\|u_{0}\|_{2}\geq \Lambda _{J(u_{0})}$ , a contradiction with (5.2). So $u(t)\in \mathcal {N}_{-}$ , $t\in [0,T_{\max })$ . Suppose $T_{\max } = \infty $ ; for every $\omega \in \omega (u_0)$ , we get $\|\omega \|_{2}>\Lambda _{J(u_{0})}$ , $J(\omega )<J(u_{0})$ , and hence, $\omega \in J^{J(u_{0})}$ and $\omega \not \in \mathcal {N}^{J(u_{0})}$ . Recalling the definition of $\Lambda _{J(u_0)}$ , we obtain $\omega (u_0)\cap \mathcal {N}=\emptyset $ . Thus, $\omega (u_0) = \left \{0\right \}$ , a contradiction to Lemma 2.7. Hence, $T_{\max } <+\infty $ .
Remark 5.1 For any $M> d$ , there exists $u_{M}\in \mathcal {N}_{-}$ such that $J(u_M) \geq M$ and the solution of (1.1) blows up in finite time. In fact, by using Theorem 5.1, if the initial data satisfy the following inequality
then the solution of (1.1) blows up in finite time.
6 Blow-up time estimates
Theorem 6.1 For $p>2$ and $J(u_{0})<0$ , the solution $u(x,t)$ blows up at finite time
The upper bound of blow-up rate is given as $ \|u\|_{2}\leq \left [3\cdot 2^{-\frac {3}{2}}\eta \right ]^{-1}(T_{*}-t)^{-1}$ , where $\eta :=|J(u_{0})|L^{-\frac {3}{2}}(0)$ and $L(t):=\frac {1}{2}\|u\|_{2}^{2}$ .
Proof Let $K(t):=-J(u(t))$ . Then $L(0)>0$ and $K(0)>0$ . By (1.6), we get $ K'(t)=\|u'\|_{2}^{2}\geq 0$ , which implies $K(t)\geq K(0)>0$ for all $t\in [0,T)$ . By (1.4), we have
Combining (6.2) with the Cauchy inequality, we obtain
According to (6.3), $ (K(t)L^{-\frac {3}{2}}(t))'=L^{-\frac {5}{2}}(t)(K'(t)L(t)-\frac {3}{2}K(t)L'(t))\geq 0$ . Therefore,
Integrating (6.4) with respect to t, we have $ \eta t\leq -\frac {2}{3} (L^{-\frac {1}{2}}(t)-L^{-\frac {1}{2}}(0) )$ , where ${\eta :=K(0)L^{-\frac {3}{2}}(0)}$ , and hence, we obtain that there exists a constant $T_{*}<+\infty $ such that $\operatorname *{\lim }\limits _{t\rightarrow T_{*}}L(t)=\infty $ (i.e., the weak solution blows up). Hence, (6.1) holds. Similarly, integrating (6.4) from t to $T_{*}$ , we have $ L(t)\leq \left [3\cdot 2^{-1}\eta \right ]^{-2}(T_{*}-t)^{-2}$ .
Theorem 6.2 For $p>2$ and $0\leq J(u_{0})<\frac {1}{6}B_{4}^{-2}\|u_{0}\|_{2}^{2}$ , the solution u blows up at finite time $T_{*} \leq \frac {16\|u_{0}\|_{2}^{2}}{B_{4}^{-2}\|u_{0}\|_{2}^{2}-6J(u_{0})}$ .
Proof We assert that for any $t\in [0,T)$ , $I(u)<0$ . From (1.4), $ I(u_{0}) \leq 3J(u_{0})-\frac {1}{2}\|\Delta u_{0}\|_{2}^{2} <0$ . If the assertion was not true, there would be a constant $t_{0}\in (0,T)$ such that $I(u)<0$ and $I(u(t_{0}))=0$ . From (6.2), we know $\|u\|_{2}^{2}$ is strictly increasing with $t\in [0,t_{0})$ . Therefore,
From the definition of $J(u)$ and $I(u(t_{0}))=0$ , we have $J(u(t_{0}))>\frac {1}{6}\|\Delta u(t_{0})\|_{2}^{2}$ . It contradicts to (6.5). Therefore, we have $I(u)<0$ for any $t\in [0,T)$ , and $\|u\|_{2}^{2}$ is strictly increasing with respect to t. For any $T\in (0,T^{*})$ , $\rho>0$ , $\xi>0$ , we define an auxiliary function
By direct calculation, we have $ G'(t) =\|u\|_{2}^{2}-\|u_0\|_{2}^{2}+2\rho (t+\xi )= \int _{0}^{t} \frac {\textrm {d}}{\textrm {d}\tau }\|u\|_{2}^{2}\textrm {d}\tau + 2\rho (t+\xi )=2 \int _{0}^{t}(u,u_{\tau })\textrm {d}\tau +2\rho (t+\xi )$ , and
It is obvious that $G"(t)>0$ . Then $G'(t)$ is increasing on $[0,T]$ and $ G'(t)\geq G'(0)>0$ . This indicates that $G(t)$ is increasing on $[0,T]$ . Define
By Cauchy-Schwarz inequality, we have $f(t)\geq 0$ . For any $t\in [0,T]$ ,
where $0<\rho \leq -\frac {3}{2}J(u_{0})+\frac {1}{4}B_{4}^{-2}\|u_{0}\|_{2}^{2}:=\hat \rho $ . From Lemma 2.10, $ T\leq \frac {T\|u_{0}\|_{2}^{2}+\rho \xi ^{2}}{\rho \xi }$ . Fixing $\rho _{0}\in (0,\hat \rho ]$ , we have $0<\frac {\|u_{0}\|_{2}^{2}}{\rho \xi }<1$ for any $\xi \in (\frac {\|u_{0}\|_{2}^{2}}{\rho _{0}},+\infty )$ . This indicates $ T\leq \frac {\rho _{0}\xi ^{2}}{\rho _{0}\xi -\|u_{0}\|_{2}^{2}}$ . The right side of the above formula takes the minimum value at $\xi :=\frac {2\|u_{0}\|_{2}^{2}}{\rho _{0}}$ . Then $ T\leq \frac {4\|u_{0}\|_{2}^{2}}{\rho _{0}}$ . For $\rho _{0}\in (0,\hat \rho ]$ , $T<T^{*} \leq \frac {16\|u_{0}\|_{2}^{2}}{B_{4}^{-2}\|u_{0}\|_{2}^{2}-6J(u_{0})}$ .
Theorem 6.3 Let $N=1$ . If u blows up in finite time in its $H_0^2 (\Omega )$ norm, there exists a positive constant such that
where $M(0):=\|\Delta u_{0}\|_{2}^{2}$ , $\tilde \alpha := {\gamma ^{2p}}/{\mu }$ , $\tilde \beta := {1}/{\lambda }$ , and $\gamma $ is the embedding constant from $H_{0}^{2}(\Omega )$ to $L^{2p}(\Omega )$ .
Proof Let $M(t):=\|\Delta u\|_{2}^{2}$ and $M(0)>0$ . When $N=1$ , problem (1.1) can be simplified to
By computation, we have
Let us estimate the integrals in (6.8). We have
with two arbitrary positive constants $\lambda $ , $\mu $ . By replacing (6.9) and (6.10) in (6.8), we obtain
Choosing $\lambda +\mu \leq 2$ , inequality (6.11) reduces to $ M'(t)\leq \tilde \alpha M^{p}(t)+\tilde \beta M(t)$ . Then we get (6.6).
7 Extinction and non-extinction
Theorem 7.1 If $p<1$ and $N=1$ , then the weak solution u of (1.1) becomes extinct in finite time if the initial data satisfies $\|u_{0}\|_{2}^{1-p}\geq D_{2}D_{1}^{-1}$ . An upper bound of extinction rate for u is
where $T_{1}:=\frac {2\|u_{0}\|_{2}^{1-p}}{(1-p)(D_{1}\|u_{0}\|_{2}^{1-p}-D_{2})}$ , $D_{1}:=B_{3}^{-2}+B_{4}^{-2}$ , and $D_{2}:=B_{5}^{p+1}$ .
Proof When $N=1$ , we get (6.7). First, we give the upper bound of extinction rate. Multiplying the equation in (6.7) by u and integrating it over $\Omega \times (t,t+h)$ with $h>0$ and then dividing the result by h yields that
Let $h\rightarrow 0^{+}$ in (7.2). Using the Lebesgue differentiation theorem in [Reference Antontsev and Shmarev2], one could obtain
where $H(t):=\|u\|_{2}^{2}$ . According to the embedding relationship $W_{0}^{1,2}(\Omega )\hookrightarrow L^{2}(\Omega )\hookrightarrow L^{p+1}(\Omega )$ and $H_{0}^{2}(\Omega )\hookrightarrow L^{2}(\Omega )$ , we have $\|u\|_{2}\leq B_{3}\|\nabla u\|_{2}$ , $\|u\|_{2}\leq B_{4}\|\Delta u\|_{2}$ , $\|u\|_{p+1}\leq B_{5}\|u\|_{2}$ . Then (7.3) can be written $ H'(t)+2D_{1}H(t)\leq 2D_{2}H^{\frac {p+1}{2}}(t)$ , where $D_{1}:=B_{3}^{-2}+B_{4}^{-2}$ , $D_{2}:=B_{5}^{p+1}$ . Defining $\varphi (t):=H^{\frac {1-p}{2}}(t)$ , we get the following differential inequality:
It is clear from (7.4) that $\zeta (0) < 0$ . Recalling the continuity of $\zeta (t)$ , there exists a sufficiently small $T_0> 0$ such that $ \zeta (t) <\frac {\zeta (0)}{2}<0$ for $0<t\leq T_{0}$ . Then $\varphi '(t) \leq {\zeta (0)}/{2}$ , which implies that
Obviously, by the definition of $\zeta (t)$ , (7.1) holds.
Corollary 7.1 If $p<1$ , $N=1$ , and $J(u_{0})\leq 0$ , then the weak solution u of (1.1) becomes extinct in finite time with the initial data satisfying $D_{2}D_{1}^{-1}\leq \|u_{0}\|_{2}^{1-p}\leq a D_{2}D_{1}^{-1}$ . The results of Theorem 7.1 hold, and a lower bound of extinction rate for u is
where $T_{2}:=\frac {1}{D_{3}}\log \frac {D_{4}D_{3}^{-1}}{D_{4}D_{3}^{-1}-\|u_{0}\|_{2}^{1-p}}$ , $D_{3}:=\frac {1-p}{2}D_{1}>0$ , $D_{4}:=\frac {(2-p)(1-p)}{p+1}D_{2}>0$ , and $a:=\frac {4-2p}{p+1}$ .
Proof The proof of a upper bound of extinction rate for u is the same as Theorem 7.1. We only give the proof of a lower bound of extinction rate. Let $H(t):=\|u\|_{2}^{2}$ . We have
By the definition of $\varphi (t)$ in Theorem 7.1, we have $ \varphi '(t)\geq D_{3}\varphi (t)-D_{4}$ ; hence, we get (7.5).
Corollary 7.2 Let $p<1$ , $N=1$ , and $\|u_{0}\|_{2}^{1-p}\geq D_{2}D_{1}^{-1}$ . If $J(u_0)> d$ and $I(u_0) = 0$ , then the results of Theorem 7.1 hold.
In the following, we give a result on non-extinction of weak solutions.
Theorem 7.2 Let $J(u_0)\leq 0$ . If one of the following conditions holds: (i) $p<1$ and $\|u_{0}\|_{2}^{1-p}>D_{4}D_{3}^{-1}$ ; (ii) $p=1$ ; (iii) $p\geq 2$ , then the solution u of (1.1) does not vanish in finite time.
Proof Let $H(t):=\|u\|_{2}^{2}$ . In the forthcoming proof, the cases $p<1$ , $p=1$ , and $p\geq 2$ will be discussed separately.
(i) $p<1$ , $\|u_{0}\|_{2}^{1-p}>D_{4}D_{3}^{-1}$ . By the similar proof of Corollary 7.1, we have
(ii) $p=1$ . By (7.6), we obtain $ H'(t)\geq D_{5}H(t)$ , where $D_{5}:=B_{3}^{-2}+B_{4}^{-2}{-B_5^2}$ . Hence, $\|u\|_{2}\geq \|u_{0}\|_{2}^{2}\textrm {e}^{\frac {D_{5}}{2}t}$ .
(iii) $p\geq 2$ . By (7.6), we obtain $ H'(t)\geq D_{6}H(t)$ , where $D_{6}:=B_{3}^{-2}+B_{4}^{-2}$ . Hence, $\|u\|_{2}\geq \|u_{0}\|_{2}^{2}\textrm {e}^{\frac {D_{6}}{2}t}$ .
Acknowledgements
The authors would like to express their sincere thanks to the editor and the reviewers for the constructive comments to improve this paper.
Funding
This paper is supported by Shandong Provincial Natural Science Foundation of China.
Competing interest
The authors declare no competing interests.
Data availability statement
All data generated or analyzed during this study are included in this article.
Use of AI tools declaration
The authors declare they have not used artificial intelligence (AI) tools in the creation of this article.