- Research
- Open access
- Published:
Global and blowup solutions of semilinear heat equation involving the square root of the Laplacian
Boundary Value Problems volume 2015, Article number: 121 (2015)
Abstract
In this paper, we consider nonlinear parabolic equations involving a nonlocal operator: the square root of the Laplacian in a bounded domain with zero Dirichlet boundary condition. We use the method on harmonic extension to study the existence and asymptotic estimates of global solutions, as well as the blowup of the parabolic equation.
1 Introduction
This paper is concerned with the study of global and blowup solutions of semilinear heat equation involving a nonlocal positive operator: the square root of the Laplacian in a bounded domain with zero Dirichlet boundary conditions. We consider the semilinear heat equation of the following form:
where Ω is a smooth boundary domain of \(\mathbb{R}^{n}\) (\(n\geq 2\)) and \((-\Delta)^{1/2}\) stands for the square root of the Laplacian operator −Δ in Ω with zero Dirichlet boundary values on ∂Ω, and \(p=2^{\#}-1=\frac{n+1}{n-1}\), \(2^{\#}=\frac{2n}{n-1}\) is the critical Sobolev exponent.
Nonlinear evolution problems involving fractional Laplacian describing the anomalous diffusion were extensively studied in the mathematical and physical literature (see [1–7] for references). This equation is
where \(0<\alpha< 2\).
However, the parabolic equation in a bounded domain, i.e. the equation
is seldom researched. The main difficulty is that the fractional operator, such as \((-\Delta)^{1/2}\), is nonlocal and nonlinear so that this problem may not possess some geometry structures, for example, the mountain pass structure, the Poincaré inequality structure, etc. In order to overcome this difficulty, we turn to another method, which was introduced by Caffarelli and Silvestre in [8] and since then has been widely used in [9–11]. They used this idea to consider the state equation of (1.3): elliptic problems involving a nonlocal operator. In this paper, we shall employ a similar method to [9] or [8] to study the parabolic equation (1.1) involving the square root of the Laplacian.
Firstly, we define lower (high)-energy initial value.
Associated to problem (1.1), the corresponding energy functional \(I:H_{0}^{1/2}(\Omega )\rightarrow \mathbb{R}\) is defined as follows:
Definition 1.1
We say that a function \(u_{0}(x)\) possesses lower-energy if \(u_{0}(x)\) satisfies \(E(u_{0}(x))< \frac{1}{2n}S^{n}\).
Otherwise, we say that \(u_{0}(x)\) possesses high-energy, where S denotes the best Sobolev constant.
We shall see the existence of global solutions with lower-energy and the difference of asymptotic behavior of blowing-up solutions with lower-energy and those with high-energy in the following results.
Remark 1.1
Lions [12] showed that
Escobar [13] proved that the extremal functions all have the form
where \(x_{0}\in R^{n}\) and \(\varepsilon>0\) are arbitrary. In addition, the best constant is
where \(\sigma_{n}\) denotes the volume of n-dimensional sphere \(R^{n}\subset R^{n+1}\).
In this paper, we consider the following weak solution.
Definition 1.2
We say that a function u is a solution of (1.1) in \(\Omega\times (0,T)\) if
and satisfies (1.1) in the distribution sense.
Remark 1.2
The analogous problem to (1.1) for the Laplacian has been investigated widely in the last decades. This is the problem
see [14] and references therein. In [14], Tan first introduced the energy method to study the existence and asymptotic estimates of global solution of (1.4) in \(R^{n}\) (\(n\geq 3\)), and gave the sufficient conditions of finite time blowup of local solution by the classical concave method. Finally, he considered the asymptotic behavior of any global solution.
For the square root \((-\Delta)^{1/2}\) of the Laplacian, we derive the following results.
Theorem 1.1
Let \(u_{0}(x)\) be a lower-energy initial value and \(E(u_{0}(x))\leq 0\). Then \(u(x,y,t;u_{0})\) blows up in finite time.
Theorem 1.2
Let \(u_{0}(x)\) (≢0) be a lower-energy initial value with \(E(u_{0}(x))> 0\).
- (1∘):
-
If \(\int_{\Omega}|u_{0}|^{2^{\#}}\,dx< S^{n}\), then (1.1) has a global solution \(u(x,t;u_{0})\).
Moreover, there exist \(K_{1}>0\) and \(K_{2}>0\) such that
$$\begin{aligned}& \|u\|_{H_{0}^{1/2}}=O\bigl(e^{-K_{1}t}\bigr),\quad t\rightarrow\infty, \end{aligned}$$(1.5)$$\begin{aligned}& \|u\|_{L^{2}}=O\bigl(e^{-K_{2}t}\bigr),\quad t\rightarrow\infty. \end{aligned}$$(1.6) - (2∘):
-
If \(\int_{\Omega}|u_{0}|^{2^{\#}}\,dx> S^{n}\), then the local solution blows up in finite time.
As \((-\Delta)^{1/2}\) is nonlocal and nonlinear, we realize problem (1.1) through a local problem in one more dimension by a Dirichlet to Neumann map. Here the Sobolev trace embedding comes into play, and its critical exponent \(2^{\#}=\frac{2n}{n-1}\), \(n\geq2\), is the power in Theorem 1.2.
2 Preliminaries
In this section we collect preliminary facts for future reference.
Let \(\{\lambda_{k},\varphi_{k}\}_{k=1}^{\infty}\) be the eigenvalues and corresponding eigenfunctions of the Laplacian operator −Δ in Ω with zero Dirichlet boundary values on ∂Ω, such that \(\|\varphi_{k}\|_{L^{2}(\Omega)}=1\). Let
Denote by \(H_{0}^{-1/2}(\Omega)\) the dual space of \(H_{0}^{1/2}(\Omega)\). \((-\Delta)^{1/2}\) is given by
for \(u= \sum_{k=1}^{\infty}\alpha_{k}\varphi_{k}\in L^{2}(\Omega)\).
Regarding (1.1) as elliptic with respect to x variables, we have, from [9],
That is, we will study the following mixed boundary value problem in a half cylinder:
where n is the unit outer normal to \(\Omega\times\{0\}\times (0,T)\), \(C=\Omega\times(0, \infty)\) and its lateral boundary is denoted by
If v satisfies (2.3), then the trace u on \(\Omega\times\{0\}\times (0,T)\) of the function v will be a solution of problem (1.1). We consider the Sobolev space of a function in \(H^{1}(\mathcal{C})\) whose traces vanish on \(\partial_{L}\mathcal{C}\),
equipped with the norm
Denote by \(\mathcal{V}_{0}(\Omega)\) the space of traces on \(\Omega\times\{0\}\) of functions in \(H_{0,L}^{1}(\mathcal{C})\):
It is easy to see that for every \(\eta\in C^{\infty}(\overline{\mathcal{C}})\cap H^{1}(\mathcal{C})\) and \(\eta\equiv0\) on \(\partial_{L}\mathcal{C}\),
Since the harmonic extension operator is bijective from \(\mathcal {V}_{0}(\Omega)\) to \(H_{0,L}^{1}(\mathcal{C})\), by using the trace theorem we can deduce the following definition.
Definition 2.1
Assume v is the harmonic extension of u (in the weak sense) to \(\mathcal{C}\) and vanishing on \(\partial_{L}\mathcal{C}\). Let us define the operator \((-\Delta)^{1/2}:\mathcal{V}_{0}(\Omega)\rightarrow \mathcal{V}_{0}^{*}(\Omega)\) by
where \(\mathcal{V}_{0}^{*}(\Omega)\) is the dual space of \(\mathcal{V}_{0}(\Omega)\).
Equation (2.3) corresponds formally to the \(L^{2}\) gradient flow associated to the energy functional
One formally sees at once that \(E(v)\) is decreasing in time along the trajectory for
i.e. \(E(v)\) is a Lyapunov functional for this flow.
In the following, we give some properties of the space \(H_{0,L}^{1}(\mathcal{C})\). Denote by \(D^{1,2}(R^{n+1}_{+})\) the closure of the set of smooth functions compactly supported in \(\overline{R}_{+}^{n+1}\) with respect to the norm of \(\|w\|_{D^{1,2}(R^{n+1}_{+})}=(\int_{R^{n+1}_{+}}|\nabla w|^{2}\,dx\,dy)^{1/2}\).
We recall the well-known Sobolev inequality. For \(w\in D^{1,2}(R^{n+1}_{+})\), we have
where C depends only on n.
The Sobolev trace inequality leads directly to the next three lemmas. For \(v\in H_{0,L}^{1}(\mathcal{C})\), its extension by zero in \(R_{+}^{n+1}\backslash\mathcal{C}\) can be approximated by functions compactly supported in \(R_{+}^{n+1}\). Thus the Sobolev trace inequality (2.6) leads to the following.
Lemma 2.1
[9]
Let \(n\geq2\) and \(2^{\#}=\frac{2n}{n-1}\). Then there exists a constant C, depending only n, such that, for all \(v\in H_{0, L}^{1}(\mathcal{C})\),
By Hölder’s inequality, since Ω is bounded, the above lemma leads to the following.
Lemma 2.2
[9]
-
(i)
Let \(1\leq q\leq2^{\#}\) for \(n\geq2\). Then we have, for all \(v\in H_{0,L}^{1}(\mathcal{C})\),
$$ \biggl(\int_{\Omega}\bigl\vert v(x,0)\bigr\vert ^{q}\,dx\biggr)^{1/q}\leq C\biggl(\int_{\mathcal{C}} \bigl\vert \nabla v(x,y)\bigr\vert ^{2}\,dx\,dy\biggr)^{1/2}, $$(2.8)where C depends only on n, q, and the measure of Ω. Moreover, (2.5) also holds for \(1\leq q< \infty\) if \(n=1\).
-
(ii)
Let \(1\leq q< 2^{\#}=\frac{2n}{n-1}\) for \(n\geq2\) and \(1\leq q< \infty\) for \(n=1\). Then \(\operatorname{tr}_{\Omega}(H_{0,L}^{1}(\mathcal{C}))\) is compactly embedded in \(L^{q}(\mathcal{C})\).
Lemma 2.3
[15]
3 Proof of Theorem 1.1
In fact, we can prove a more general result. If there exists some \(t_{0}\) such that \(E(v(t_{0}))\leq0\), then \(v(x,y,t;u_{0})\) blows up in finite time. We shall employ the classical concavity method (see [14, 16, 17]). Suppose that \(t_{\mathrm{max}}=\infty\) and denote \(f(t)=\frac{1}{2}\int_{t_{0}}^{t}\int_{\Omega\times \{0\}}|v|^{2}\,dx\,ds\).
We perform standard manipulations:
From the assumption, \(E(v(t_{0}))\leq0\) such that
for all \(t\geq t_{0}\). If we had \(t_{\mathrm{max}}=\infty\), this inequality would yield
and
and as \(t\rightarrow\infty\) we have, for some \(\alpha> 0\) and \(\forall t\geq t_{0}\),
Hence \(f^{-\alpha}(t)\) is concave on \([t_{0},\infty]\), \(f^{-\alpha}(t)> 0\), and \(\lim_{t\rightarrow\infty}f^{-\alpha}(t)=0\). This contradiction proves that \(t_{\mathrm{max}}< \infty\), which completes the proof of Theorem 1.1.
4 Proof of Theorem 1.2
We divide the proof into several steps.
Step 1: Proof of existence.
(i) As regards a priori estimates and local existence, assume \(v_{n_{0}}\in H_{0,L}^{1}(\mathcal{C})\) such that
and
On the other hand, multiplying (2.3) by \(v_{nt}\) and integrating, we have
For the sake of convenience, define
Now, we show that \(v_{n}(t)\in\Sigma\), for any \(t\geq0\). Suppose that it does not hold and let \(t^{*}\) be the smallest time for which \(v_{n}(t^{*})\notin\Sigma\). Then in virtue of the continuity of \(v_{n}(t)\), we see that \(v_{n}(t^{*})\in\partial\Sigma\). Hence
which contradicts (4.1). Then from (4.1) and noting that if
then
we have
Thus, we obtain
From (4.4), we have
where \(C(T)\) is the constant independent of n.
Using the trace theorem, we see, from the prior estimates (4.3) and (4.4), that there exist a subsequence (not relabeled) and a function v such that
In particular, for every \(\varphi\in H_{0, L}^{1}(\mathcal{C})\), we obtain, from Lemma 2.2, that, as \(n\rightarrow\infty\),
which implies that the function v is a desired local solution of (2.3) and \(v\in H_{0, L}^{1}(\mathcal{C})\).
(ii) As regards global existence, multiplying (2.3) by \(v_{t}\) and integrating, we obtain
Thus,
for any \(t> 0\).
Note if
then
Now we prove that \(v(x,0,t)\in\partial\Sigma\) for any \(t> 0\), then we have
which is a contradiction. Hence
for any \(t>0\). Therefore,
which implies
for any \(T> 0\). Thus, \(v(x,y,t)\) is a global solution of (2.3).
Step 2: Proof of (1.5).
We apply the same argument as in [18], and for the sake of completeness, we give the proof. Let
then by Step 1, we have \(h(v(t))> 0\), for all \(t\geq0\).
As for the Sobolev trace inequality, we have
and the inequality
implies
For simplicity, denote \(\frac{1}{S^{\frac{2^{\#}}{2}}}(2n E(v_{0}))^{\frac{2^{\#}}{2}-1}\) by δ (\(0<\delta< 1\)). Letting \(\gamma=1-\delta\), we have
Let \(T> t_{0}\) be a fixed number, then
Using Lemma 2.2, there exists a positive constant \(C_{1}\) such that
Furthermore, inequality (4.10) implies
on \([t_{0}, \infty)\).
Therefore, by inequalities (4.11) and (4.12), we obtain
on \([t_{0},T]\). On the other hand, inequality (4.10) implies
on \([t_{0}, \infty)\).
By inequalities (4.12) and (4.14), we have
Furthermore, (4.13) and (4.15) give
on \([t_{0},T]\).
Then, from the arbitrariness of \(T>t_{0}\), we have
Let \(T_{0}>t_{0}\) be sufficiently large such that \(C_{1}^{-1}\leq T_{0}\), it follows that
on \([t_{0}, \infty)\).
Setting \(y(t)=\int_{t}^{\infty} E(v(s))\,ds\), it follows from (4.16) that
By (4.12), we have
for some constant \(C_{3}\) for large \(t>T_{0}\). Hence
From Lemma 2.3, we have
Step 3: Proof of (1.6).
Obviously
and
for all \(t>0\). By the same argument as Step 2, we have
We see that the estimate
holds. That is,
Step 4: Proof that if \(\int_{\Omega}|u_{0}|^{2^{\#}}\,dx> S^{n}\), then the local solution blows up in finite time.
We divide the proof into two steps.
(i) First of all, we define a set which consists of the functions that satisfy the following conditions:
We claim that the set is an empty set.
Indeed, let \(v_{0}\) belong to the set. If \(v_{0}\) satisfies
then
and hence
which is contradictory to condition (4.18).
If \(v_{0}\) satisfies \(\int_{\mathcal{C}}|\nabla v_{0}(x,y)|^{2}\,dx\,dy> \int_{\Omega\times \{0\}}|v_{0}(x,0)|^{2^{\#}}\,dx\), then from inequality (4.18), we see that
It implies \(\int_{\Omega\times \{0\}}|v_{0}(x,0)|^{2^{\#}}\,dx< S^{n}\) which is a contradiction because of condition (4.19). Therefore, that set is an empty set.
(ii) Thus, we consider only the following case:
Obviously, in this case we have \(S^{n}< \int_{\mathcal{C}}|\nabla v_{0}(x,y)|^{2}\,dx\,dy< \int_{\Omega\times \{0\}}|v_{0}(x,0)|^{2^{\#}}\,dx\). If \(v(x,y,t)\) is a global solution, then we can deduce that \(v(x,y,t)\) does not converge strongly to 0 in \(H_{0,L}^{1}(\mathcal{C})\). Otherwise, there will exist a \(t^{*}\) (\(t^{*}> 0\)) such that
which contradicts part (i).
To complete the proof of Theorem 1.2(2∘), we first prove the following.
Claim
\(v_{0}\) satisfies (4.20) and \(v(x,y,t;v_{0})\) is a global solution. For \(\forall t\in[0,T]\) the following inequalities hold:
Indeed, if there exists a \(t^{*}\in[0,T_{\mathrm{max}}]\) such that
then we have
But
which is a contradiction. Therefore there exists a constant \(\eta>0\) sufficiently small and independent of t, rely on \(v_{0}\) such that
for any \(t\in[0,\infty]\), which completes the proof of the claim.
Now we can complete the proof of Theorem 1.2(2∘). We shall employ the same argument as the proof of Theorem 1.1. Suppose that \(T_{\mathrm{max}}=\infty\) and denote \(f(t)=\frac{1}{2}\int_{t_{0}}^{t}\int_{\Omega\times \{0\}}|v|^{2}\,dx\,ds\). We obtain (3.1), (3.3), and from (3.3) and (4.22) we have
If we have \(T_{\mathrm{max}}=\infty\), then this inequality would yield
and by (4.20), we have
which implies
and
By the argument of the proof of Theorem 1.1, we obtain a contradiction, which completes the proof of Theorem 1.2(2∘).
References
Biler, P, Karch, G, Woyczynski, WA: Critical nonlinearity exponent and self-similar asymptotics for Lévy conservation laws. Ann. Inst. Henri Poincaré, Anal. Non Linéaire 18, 613-637 (2001)
Fino, A, Karch, G: Decay of mass for nonlinear equation with fractional Laplacian (2006). arXiv:0812.4977v1 [math.AP]
Guo, B, Huo, Z: Global well-posedness for the fractional nonlinear Schrödinger equation. Commun. Partial Differ. Equ. 36, 247-255 (2011)
Guo, B, Han, Y, Xin, J: Existence of the global smooth solution to the period boundary value problem of fractional nonlinear Schrödinger equation. Appl. Math. Comput. 204, 468-477 (2008)
Li, J, Wang, X, Yao, Z: Heat flow for the square root of the negative Laplacian for unit length vectors. Nonlinear Anal. TMA 68, 83-96 (2008)
Silvestre, L: Regularity of the obstacle problem for a fractional power of the Laplace operator. Commun. Pure Appl. Math. 60, 67-112 (2006)
Tan, Z, Xu, Y: Existence and nonexistence of global solutions for a semi-linear heat equation with fractional Laplacian. Acta Math. Sci., Ser. B 28(6), 818-825 (2012)
Caffarelli, L, Silvestre, L: An extension problem related to the fractional Laplacian. Commun. Partial Differ. Equ. 32, 1245-1260 (2007)
Cabre, X, Tan, J: Positive solutions of nonlinear problems involving the square root of the Laplacian. Adv. Math. 224, 2052-2093 (2010)
Tan, JG: The Brezis-Nirenberg type problem involving the square root of the Laplacian. Calc. Var. 42, 21-41 (2011)
Yu, XH: The Nehari manifold for elliptic equation involving the square root of the Laplacian. J. Differ. Equ. 252, 1283-1308 (2012)
Lions, PL: The concentration-compactness principle in the calculus of variations. The limit case II. Rev. Mat. Iberoam. 1, 45-121 (1985)
Escobar, J: Sharp constant in a Sobolev trace inequality. Indiana Univ. Math. J. 37, 687-698 (1988)
Tan, Z: Global solution and blowup of semilinear heat equation with critical Sobolev exponent. Commun. Partial Differ. Equ. 26, 717-741 (2001)
Brändle, C, Colorado, E, De Pablo, A, Sánchez, U: A concave-convex elliptic problem involving the fractional Laplacian. Proc. R. Soc. Edinb., Sect. A 143, 39-71 (2013)
Cazenave, T, Lions, PL: Solutions globales d’équations de la chaleur semi linéaires. Commun. Partial Differ. Equ. 9, 955-978 (1984)
Ishii, H: Asymptotic stability and blowing up of solutions of some nonlinear equations. J. Differ. Equ. 26, 291-319 (1977)
Ikehata, R, Suzuki, T: Semilinear parabolic equations involving critical Sobolev exponent: local and asymptotic behavior of solutions. Preprint
Acknowledgements
The authors would like to thank the anonymous referees for their valuable suggestions, which helped to improve the presentation of the paper. This research is supported by National Natural Science Foundation of China (11271305).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors contributed equally and significantly in writing this article. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Xu, Y., Tan, Z. & Sun, D. Global and blowup solutions of semilinear heat equation involving the square root of the Laplacian. Bound Value Probl 2015, 121 (2015). https://doi.org/10.1186/s13661-015-0384-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-015-0384-7