We consider the semilinear heat equation
in the whole space , where and . Unlike the standard case , this equation is not scaling invariant. We construct for this equation a solution which blows up in finite time only at one blowup point , according to the following asymptotic dynamics:
where is the unique positive solution of the ODE
The construction relies on the reduction of the problem to a finite dimensional one and a topological argument based on the index theory to get the conclusion. By the interpretation of the parameters of the finite dimensional problem in terms of the blowup time and the blowup point, we show the stability of the constructed solution with respect to perturbations in initial data. To our knowledge, this is the first successful construction for a genuinely non-scale invariant PDE of a stable blowup solution with the derivation of the blowup profile. From this point of view, we consider our result as a breakthrough.
A non-scaling invariant semilinear heat equation]Construction of a stable blowup solution with a prescribed behavior for a non-scaling invariant semilinear heat equation
G. K. Duong, V. T. Nguyen, H. Zaag ]
\subjclassPrimary: 35K55, 35B40; Secondary: 35L65, 35K57.
G. K. Duong]email@example.com
V. T. Nguyen]Tien.Nguyen@nyu.edu
H. Zaag] Hatem.Zaag@univ-paris13.fr
July 3, 2019
Giao Ky Duong 111 G. K. Duong is fully funded by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreement No 665850., Van Tien Nguyen and Hatem Zaag 222 H. Zaag is supported by the ANR projet ANAÉ ref. ANR -13-BS01-0010-03
Université Paris 13, Sorbonne Paris Cité, LAGA, CNRS (UMR 7539), F-93430, Villetaneuse, France.
New York University in Abu Dhabi, P.O. Box 129188, Abu Dhabi, United Arab Emirates.
We are interested in the semilinear heat equation
where stands for the Laplacian in and
By standard results the model (1.1) is wellposed in thanks to a fixed-point argument. More precisely, there is a unique maximal solution on with . If , then the solution of (1.1) may develop singularities in finite time , in the sense that
In this case, is called the blowup time of . Given , we say that is a blowup point of if and only if there exists as such that as .
In the special case the equation (1.1) becomes the standard semilinear heat equation
This equation is invariant under the following scaling transformation
An extensive literature is devoted to equation (1.3) and no rewiew can be exhaustive. Given our interest in the construction question with a prescribed blowup behavior, we only mention previous work in this direction.
(note that Herrero and Velázquez [Herrero and Velázquez(1992)] proved the same result with a different method; note also that Bressan [Bressan(1992)] made a similar construction in the case of an exponential nonlinearity).
Later, Merle and Zaag [Merle and Zaag(1997)] (see also the note [Merle and Zaag(1996)]) simplified the proof of [Bricmont and Kupiainen(1994)] and proved the stability of the constructed solution verifying the behavior (1.5). Their method relies on the linearization of the similarity variables version around the expected profile. In that setting, the linearized operator has two positive eigenvalues, then a non-negative spectrum. Then, they proceed in two steps:
Reduction of an infinite dimensional problem to finite dimensional one: they show that controlling the similarity variable version around the profile reduces to the control of the components corresponding to the two positive eigenvalues.
Then, they solve the finite dimensional problem thanks to a topological argument based on index theory.
The method of Merle and Zaag [Merle and Zaag(1997)] has been proved to be successful in various situations. This was the case of the complex Ginzgburg-Landau equation by Masmoudi and Zaag [Masmoudi and Zaag(2008)] (see also Zaag [Zaag(1998)] for an ealier work) and also for the case of a complex semilinear heat equation with no variational structure by Nouaili and Zaag [Nouaili and Zaag(2015)]. We also mention the work of Tayachi and Zaag [Tayachi and Zaag(2016)] (see also the note [Tayachi and Zaag(2015)]) and the work of Ghoul, Nguyen and Zaag [Ghoul et al.(2017)Ghoul, Nguyen, and Zaag] dealing with a nonlinear heat equation with a double source depending on the solution and its gradient in a critical way. In [Ghoul et al.(2016)Ghoul, Nguyen, and Zaag], Ghoul, Nguyen and Zaag successfully adapted the method to construct a stable blowup solution for a non variational semilinear parabolic system.
In other words, the method of [Merle and Zaag(1997)] was proved to be efficient even for the case of systems with non variational structure. However, all the previous examples enjoy a common scaling invariant property like (1.4), which seemed at first to be a strong requirement for the method. In fact, this was proved to be untrue.
As matter of fact, Ebde and Zaag [Ebde and Zaag(2011)] were able to adapt the method to construct blowup solutions for the following non scaling invariant equation
These conditions ensure that the perturbation turns out to exponentially small coefficients in the similarity variables. Later, Nguyen and Zaag [Nguyen and Zaag(2016b)] did a more spectacular achievement by addressing the case of stronger perturbation of (1.3), namely
where and . When moving to the similarity variables, the perturbation turns out to have a polynomial decay. Hence, when is small, we are almost in the case of a critical perturbation.
In both cases addressed in [Ebde and Zaag(2011)] and [Nguyen and Zaag(2016b)], the equations are indeed non-scaling invariant, which shows the robustness of the method. However, since both papers proceed by perturbations around the standard case (1.3), it is as if we are still in the scaling invariant case.
In this paper, we aim at trying the approach on a genuinely non-scaling invariant case, namely equation (1.1). This is our main result.
Theorem 1.1 (Blowup solutions for equation (1.1) with a prescribed behavior).
There exists an initial data such that the corresponding solution to equation (1.1) blows up in finite time only at the origin. Moreover, we have
For all , there exists a positive constant such that
where is the unique positive solution of the following ODE
(see Lemma A.1 for the existence and uniqueness of ), and the profile is defined by
There exits such that uniformly on compact sets of , where
From , we see that as , which means that the solution blows up in finite time at . From , we deduce that the solution blows up only at the origin.
Note that the behavior in (1.8) is almost the same as in the standard case treated in [Bricmont and Kupiainen(1994)] and [Merle and Zaag(1997)]. However, the final profile has a difference coming from the extra multiplication of the size , which shows that the nonlinear source in equation (1.1) has a strong effect to the dynamic of the solution in comparing with the standard case .
By the parabolic regularity, one can show that if the initial data , then we have for ,
where is defined by (1.10).
From the technique of Merle [Merle(1992)], we can prove the following result.
As a consequence of our technique, we prove the stability of the solution constructed in Theorem 1.1 under the perturbations of initial data. In particular, we have the following result.
Theorem 1.7 (Stability of the solution constructed in Theorem 1.1).
Consider the solution constructed in Theorem 1.1 and denote by its blowup time. Then there exists a neighborhood of such that for all , equation (1.1) with the initial data has a unique solution blowing up in finite time at a single point . Moreover, the statements and in Theorem 1.1 are satisfied by , and
2 Formulation of the problem.
In this section, we first use the matched asymptotic technique to formally derive the behavior (1.8). Then, we give the formulation of the problem in order to justify the formal result.
2.1 A formal approach.
In this part, we follow the approach of Tayachi and Zaag [Tayachi and Zaag(2016)] to formally explain how to derive the asymptotic behavior (1.8). To do so, we introduce the following self-similarity variables
Note that admits the following asymptotic behavior as ,
(see ii) of Lemma A.5 for the proof of (2.5)). From (2.1), we see that the study of the asymptotic behavior of as is equivalent to the study of the long time behavior of as . In other words, the construction of the solution , which blows up in finite time and verifies the behavior (1.8), reduces to the construction of a global solution for equation (2.2) satisfying
In the following, we will formally explain how to derive the behavior (2.7).
2.1.1 Inner expansion
We remark that are the trivial constant solutions to equation (2.2). Since we are looking for a non zero solution, let us consider the case when as . We now introduce
then from equation (2.2), we see that satisfies
(see Lemma A.6 for the proof of this statement).
Since as and the nonlinear term is quadratic in , we see from equation (2.9) that the linear part will play the main role in the analysis of our solution. Let us recall some properties of . The linear operator is self-adjoint in , where is the weighted space associated with the weight defined by
More precisely, we have
When , all the eigenvalues of are simple and the eigenfunction corresponding to the eigenvalue is the Hermite polynomial defined by
In particular, we have the following orthogonality
When , the eigenspace corresponding to the eigenvalue is defined as follows
Since the eigenfunctions of is a basic of , we can expand in this basic as follows
For simplicity, let us assume that is radially symmetric in . Since with corresponds to negative eigenvalues of , we may consider the solution taking the form
as . we now assume that as , then (2.17) becomes
We consider the following cases:
- Case 1: Either or as , then the second equation in (2.17) becomes
which contradicts with the condition as .
- Case 2: as , then (2.17) becomes
from which we improve the error for as follows
in as . Note that the asymptotic expansion (2.20) also holds for all , K is an arbitrary positive number.
2.1.2 Outer expansion.
The asymptotic behavior of (2.20) suggests that the blowup profile depends on the variable
where is a bounded, smooth function to be determined. From (2.20), we impose the condition
Since is supposed to be bounded, we obtain from Lemma A.7 that
Note also that
from which and the asymptotic behavior (2.20), we find that
In conclusion, we have just derived the following asymptotic profile
2.2 Formulation of the problem.
In this subsection, we set up the problem in order to justify the formal approach presented in the Section 2.1. In particular, we give a formulation to prove item of Theorem 1.1. We aim at constructing for equation (1.1) a solution blowing up in finite time only at the origin and verifying the behavior (1.8). In the similarity variables (2.1), this is equivalent to the construction of a solution for equation (2.2) defined for all and satisfying (2.7). The formal approach given in subsection 2.1 (see (2.25)) suggests to linearize around the profile function defined by (2.26). Let us introduce
where is the linear operator defined by (2.10) and
Since we construct for equation (2.28) a solution verifying as , and the fact that
(see Lemmas A.8, A.9 and A.10 for these estimates), we see that the linear part of equation (2.28) will play an important role in the analysis of the solution. The property of the linear operator has been studied in previous section (see page 2.13), and the potential has the following properties:
Perturbation of effect of inside the blowup region :
For each , there exist and such that
Since is the biggest eigenvalue of , the operator behaves as one with with a fully negative spectrum outside blowup region , which makes the control of the solution in this region easily.
Since the behavior of the potential inside and outside the blowup region is different, we will consider the dynamics of the solution for and for separately for some K to be fixed large. We introduce the following function
and is a positive constant to be fixed large later. We now decompose by
(Note that and ). Since the eigenfunctions of span the whole space , let us write
In particular, we denote and is a symmetric matrix defined explicitly by
Note that and are the components of , and not those of .
3 Proof of the existence assuming some technical results.
In this section, we shall describe the main argument behind the proof of Theorem 1.1. To avoid winding up with details, we shall postpone most of the technicalities involved to the next section. According to the transformations (2.1) and (2.27), proving of Theorem 1.1 is equivalent to showing that there exists an initial data at the time such that the corresponding solution of equation (2.28) satisfies
In particular, we consider the following function
We aim at proving that there exists such that the solution of (2.28) with initial data satisfies
More precisely, we will show that there exists such that the solution belongs to the shrinking set defined as follows:
Definition 3.1 (A shrinking set to zero).
For all we define being the set of all functions such that
where , , , and are defined as in (2.41).
We also denote by being the set
For each , we have the following estimates for all :
We aim at proving the following central proposition which implies Theorem 1.1.
Proposition 3.3 (Existence of a solution trapped in ).
From (3.5), we see that once Proposition 3.3 is proved, item of Theorem 1.1 directly follows. In the following, we shall give all the main arguments for the proof of this proposition assuming some technical results which are left to the next section.
As for the initial data at time defined as in (3.1), we have the following properties.
Proposition 3.4 (Properties of the initial data (3.1)).
For each , there exists such that for all we have the following properties:
There exists such that the mapping