We study the nonlinear Neumann boundary value problem for semilinear heat equation
where , and is small enough. We investigate the life span of solutions for . Also we study the global in time existence and large time asymptotic behavior of solutions in the case of and .
We consider the nonlinear Neumann boundary value problem for semilinear heat equations
where , , , , for , for , the nonlinearities in the equation and in the Neumann type boundary condition are the following
or
and
We study (1.1) with the critical powers of the nonlinearities .
The integral equation associated with (1.1) is written as follows
(see Appendix A.1), where
with
We also define the Green operator for the case of the Dirichlet type boundary conditions
When , then we can extend to all space as an even function with respect to , then (1.1) converts to the Cauchy problem
which solution can be written by the Duhamel formula
where the Green operator is
When , the large time asymptotics and the life span of solutions to (1.6) was studied in [3, 5, 10–13] in the case of . In [3], existence of blow up of solutions for , and global in time existence of solutions for , were shown. The case was studied in [5] for and in [10, 12] for general n and it was shown that solutions blow up in a finite time. In [11], the global existence, large time behavior or life span of solutions were studied for (1.6) with initial data . In particular, when , , it was shown that the upper bound of the existence time (life span) T is given by
for small and when , , global in time existence of solutions was shown for small data. In [13], global solutions were obtained if and for .
In the case of , , , (1.5) was studied in [6] and it was shown that the large time asymptotic behavior of solutions differs from that of the linear problem under the condition such that . Therefore the power is considered as critical.
We note that the estimates
and
are expected. Therefore if and , then the right hand sides of the above estimates are finite and we find that it is possible to apply the standard contraction mapping principle to show a global existence of small solutions for .
Our purpose in the present paper is to consider and show a lower bound for the existence time of local solutions in the case of and to prove the global in time existence and find the large time asymptotic behavior of solutions in case . Our method can be applied to the case of and for , or and for , .
On the other hand, when , the integral equation (1.4) was studied in [2, 4, 8] and [9].
In [2], a system of heat equations coupled in the boundary conditions such that
was studied, where
We assume ,
Then they showed if , all nonnegative solutions are nonglobal, whereas if , small global solutions exist. When , , then we have a single equation . In the case, their result says if , all nonnegative solutions are nonglobal, while if , small global solutions exist (see [4]). Therefore the power is considered as critical. In [9], necessary and sufficient conditions for the solvability of the problem were obtained, and the strongest singularity of the initial function for which we have existence of solutions was identified. In [8], the result of [9] was applied to get sharp estimates of the life span including the result of Theorem 1.1 when .
In order to state the first result, we define the function space of solutions as
with the norm
where and
Let be a closed ball in with a radius ρ and the center at the origin
We state the local in time existence for the critical case.
We assume that, then there exist anand positive constant c such that the integral equation (
1.4
) with (
1.2
) or (
1.3
) has a unique solutionwhere
By papers [5, 10–12] in the case of (1.2), we can find that the upper bound of the life span T is given by . Indeed, since , when , we have the estimate
Then if , , using the above inequality we can obtain the life span of solutions in the same way as in the proofs of these papers. Therefore there exists a gap between and . However their proofs work well for the boundary problem since we have the integral inequality
which is obtained by the conditions . Therefore we will find (see Appendix A.2 below) that the upper bound for the existence time of solutions is given by . This implies that the result of Theorem 1.1 is sharp (see also [8]).
In order to state the global existence result, we introduce the function space
where
Our basic function space is defined by
with the norm , where
where ,
and is the closed ball in with a radius ρ and the center at the origin. We also introduce the metric in , where and
with .
We now state our main result of this paper. Denote .
We assume that,,. We supposefor the case of (
1.2
) andfor the case of (
1.3
). Then there exists ansuch that the integral equation (
1.4
) with (
1.2
) or (
1.3
) has a unique global solutionMoreover, there exist positive constantssuch that the estimateholds for any, where.
In case and , our proof works well for (1.4) with (1.2) or (1.3) under the conditions of Theorem 1.2 and we have the time decay estimate for solutions such that
for some . This fact implies that dissipative effect from the boundary is stronger than that of the nonlinear forcing term.
In the case of (1.4) with (1.3), we have a-priori estimates for local solutions through the energy method. More precisely, we have
and
Therefore we expect the global in time existence of solutions in the case of large data under some growth conditions on nonlinearities and regularity on the data. However it seems that it is difficult to get time decay estimates presented in Theorem 1.2 by the energy method only.
We organize our paper as follows. In Section 2, we prove some estimates of the evolution operator associated with homogeneous Neumann problem. Section 3 is devoted to the proof of Theorem 1.1. In Section 4, we translate the integral equation of the original problem to another one such that the integral of the nonlinearities over vanishes. In Lemma 4.1, we consider some estimates of nonlinear terms of the transformed integral equation. In order to get an additional time decay of solutions we prepare Lemma 4.2. We need Lemma 4.3 to apply the contraction mapping principle. Section 5 is devoted to the proof of Theorem 1.2. In Appendix A.1, we show that the integral representation of solutions satisfies the boundary and initial conditions. Appendix A.2 is devoted to the estimate of solutions by using the method of [5] which says that the result of Theorem 1.1 is sharp (see also [8]).
We first consider the estimates of . We extend to as an even function with respect to such that
and define by
Then we can write . Also note that
Using the relation and applying Lemma 1.28 from book [7], p. 27, we get the following result.
The following estimates are trueandfor,provided that the right-hand sides are finite.
As a consequence we obtain.
Assume that. Then the estimate is true
By (2.1) and (2.2) from Lemma 2.1, we have the estimates
and
for , . Therefore . This completes the proof of the lemma. □
In order to prove the result for the critical case, we need the estimates of and .
Assume that. Then the estimates are valid
By Lemma 2.1
where . By the Gagliardo–Nirenberg–Sobolev interpolation inequality
and using estimates
and
we find
Therefore
Hence for , we find
if , and for , we have
from which it follows that
Similarly, we find
Note that for , then as above we get
Next by Lemma 2.1 we get
Hence
By the Gagliardo–Nirenberg–Sobolev interpolation inequality
and using estimates
and
we find
Therefore
As above we have for , hence
for and for
Similarly, we find
Then as above we get
This completes the proof of the lemma. □
To prove the local in time existence of solutions, we consider
then by Lemma 2.1 and Lemma 3.1
and by Lemma 3.2 for
We put , then Theorem 1.1 follows via the contraction mapping principle.
In this section, we consider the problem under the conditions such that
Namely, we consider (1.1) with (1.2). We translate the equation (1.1) to another one by using the changing of the dependent variable , , which is initiated in [6]. By a direct calculation, we obtain
Integrating equation (4.1) in , we find
Using the boundary condition , we get
Then it follows that
If we now choose by the condition , i.e.
then we find
and
The corresponding integral equation is written as
where
We denote also
and
then we obtain the following integral equation
where
and
Note that
and
Indeed, integrating (4.2) over , we obtain
since the third term on the right-hand side of (4.2) vanishes . We put
then we can write the above equation as with initial condition , where
Integration in time yields
Therefore we have the conservation law such that
for all . Below we will study equation (4.2) along with the equations for φ
In the case of nonlinearities (1.3), analogously we have
and
and also equations for φ
instead of (4.3).
Next we prove the following result.
Assume that,Also suppose thatfor any. Then
We let (analogously for the case of the nonlinearities (1.3)) and
Then since , which is true by virtue of the condition , we have
Then by Lemma 2.3 we obtain for ,
Hence
By the Gagliardo–Nirenberg–Sobolev interpolation inequality, we find
Then we compute
We extend ϕ to as an even function with respect to to see by the Sobolev embedding inequality (see Theorem 2.44 in [1])
with
Therefore
Using the estimate
for and , we obtain form (4.5) and (4.6)
By (4.4) and (4.7) we find
Then we get from (4.8)
for . We also have for
By (4.9) and (4.10) we obtain
from which it follows that
In the same way as in the proof of (4.11) we have
We take or in (4.11) and (4.12), then we find that
We next consider
or
In what follows, we consider the case (1.2) since the case (1.3) can be treated in the same way. We write
Then by Lemma 2.3 and using
for , , , we get
We take or to find that
In the same way as in the proof of (4.15) we have
By (4.15) and (4.16)
We next consider the estimate of . We again use Lemma 2.3 to obtain with
and
By (4.14)
Therefore by (4.18), (4.19) and the fact that , it follows that
and similarly,
We take or , to find that
from which it follows that
We have the result of the lemma by (4.13), (4.17) and (4.20). □
We consider the estimate of which is defined by
with
or
with
Multiplying both sides of (4.21) by , we have the nonlinear ordinary differential equations
Local in time of existence can be obtained in the standard way. A-priori estimate of local solution of is obtained by the inequality
which implies global existence of solutions. Namely, for any , we have a global solution . Next we prove the two-sided estimate for .
We assume that,,,,andFurthermore, we suppose thatThen there exist anand positive constants,such that the solutionof (
4.21
) satisfies the estimatesfor all.
We start with (4.23). We find that there exists a unique solution . Since , we have
which implies that
Hence
namely,
any , from which it follows that
We apply the estimate to (4.25) to get
for any . We assume that . By a direct calculation we obtain
and
Hence we have
which implies there exists an such that
Integrating in time, we obtain with (4.27)
Therefore we find that there exist an and such that
for . We also have by (4.27)
for . Hence for any , there exists such that
By (4.26), there exists such that
This completes the proof of the lemma. □
We define and
and the sequence as follows:
and
We only treat the case of (1.2), for simplicity. From Lemma 4.2, we find that for any , there exists a unique satisfying the estimate (4.24). Therefore is defined by and so the mapping M is well defined. Since for , by Lemma 2.3 and Lemma 4.1 we have the estimate
for small . We consider the difference between and
where
Assume that. Then
We only consider the case of (1.2), for simplicity. The case of (1.3) is treated in the same way. We have
from which with Lemma 2.1
Since
we have
By the mean value theorem we find that there exists such that
Therefore by Lemma 4.2
which implies
Also
Since for any
we have
for . We also have
Therefore
We consider the estimate of . We have
from which it follows that
Hence
from which it follows that
Therefore
In the same way,
Thus we have the result of the lemma. □
We now consider the linearized version of problem (4.2) which is written as
with the condition such that
We assume that , then by Lemma 4.1, we have
and for we find by Lemma 4.3
Therefore we have the result by the contraction mapping principle.
Footnotes
Acknowledgements
We would like to thank anonymous referee for useful comments and letting us know some references which are related to nonlinear Neumann boundary conditions. The work of N.H. is partially supported by JSPS KAKENHI Grant Numbers JP20K03680, JP19H05597. The work of E.I.K. is partially supported by CONACYT 252053-F. The work of P.I.N. is partially supported by CONACYT project 283698-F and PAPIIT project IN103221. The work of T.O. is partially supported by JSPS KAKENHI Grant Number JP19H05597.
Appendix
References
1.
H.Bahouri, J.-Y.Chemin and R.Danchin, Fourier Analysis and Nonlinear Partial Differential Equations, Springer, Berlin, 2011.
2.
K.Deng, M.Fila and H.A.Levine, On critical exponents for a system of heat equations coupled in the boundary, Acta Math. Univ. Comentian63 (1994), 169–192.
3.
H.Fujita, On the blowing up of solutions of the Cauchy problem for , J. Fac. Sci. Univ. Tokyo (I)13 (1966), 109–124.
4.
V.A.Galaktionov and H.A.Levine, On critical Fujita exponents for heat equations with nonlinear flux conditions on the boundary, Isral J. Math.94 (1996), 125–146. doi:10.1007/BF02762700.
5.
K.Hayakawa, On nonexistence of global solutions of some semilinear parabolic equations, Proc. Japan Acad.49 (1973), 503–505.
6.
N.Hayashi, E.I.Kaikina and P.I.Naumkin, Large time behavior of solutions to the dissipative nonlinear Schrödinger equation, Proc. Royal Soc. Edinburgh130A (2000), 1029–1043.
7.
N.Hayashi, E.I.Kaikina, P.I.Naumkin and I.A.Shishmarev, Asymptotics for Dissipative Nonlinear Equations, Lecture Notes in Mathematics, Vol. 1884, Springer-Verlag, Berlin, 2006, xii+557 pp.
8.
K.Hisa, Sharp estimate of the life span of solutions to the heat equation with a nonlinear boundary condition, in: Geometric Properties for Parabolic and Elliptic PDE’s, Springer INdAM Series, Vol. 47, 2021, pp. 127–149. doi:10.1007/978-3-030-73363-6_7.
9.
K.Hisa and K.Ishige, Solvability of the heat equation with a nonlinear boundary condition, SIAM J. Math. Anal.51(1) (2019), 565–594. doi:10.1137/17M1131416.
10.
K.Kobayashi, T.Sirao and H.Tanaka, On the blowing up problem for semilinear heat equations, J. Math. Soc. Japan29 (1977), 407–424. doi:10.2969/jmsj/02920373.
11.
T.Y.Lee and W.M.Ni, Global existence, large time behavior and life span of solutions of a semilinear parabolic Cauchy problem, Trans. Amer. Math. Soc.333(1) (1992), 365–378. doi:10.1090/S0002-9947-1992-1057781-6.
12.
S.Sugitani, On nonexistence of global solutions for some nonlinear integral equations, Osaka J. Math.12 (1975), 45–51.
13.
F.B.Weissler, Existence and nonexistence of global solutions for a semilinear heat equations, Israel J. Math.38 (1981), 29–40. doi:10.1007/BF02761845.