1. Introduction
The scalar form of the Kuramoto–Sivashinsky equation is
This is taken with initial data
The spatial domain we consider is the n-dimensional torus $\mathbb{T}^{n},$ which is given by
for some given lengths $L_{i} \gt 0,$ $i\in\{1,\ldots,n\}$, with periodic boundary conditions. Equation (1.1) was introduced separately by Kuramoto and Tsuzuki when studying pattern formation in reaction–diffusion equations [Reference Kuramoto and Tsuzuki24] and by Sivashinsky in modelling the evolution of flame fronts [Reference Sivashinsky34]. As a model of flame fronts, the surface $(\vec{x},\phi(\vec{x},t))$ gives the location of the interface between the burnt and unburnt phases of a gas undergoing combustion. As such, the physical cases are n = 1 (representing the interface between 2 two-dimensional gases) and n = 2 (representing the interface between 2 three-dimensional gases).
Demonstrating local well-posedness of the initial value problems (1.1) and (1.2) is straightforward for relatively smooth data. For global well-posedness, the situation is only clear in dimension n = 1. In this case, the nonlinearity has a simpler structure, and there are many results, especially those of Tadmor [Reference Tadmor35], as well as those of Bronski and Gambill [Reference Bronski and Gambill10], Goodman [Reference Goodman19], and Nicolaenko, Scheurer, and Temam [Reference Nicolaenko, Scheurer and Temam30]. All of these articles assume the same regularity on the initial data, which is that $\phi_{0}\in H^{1}.$
In two space dimensions, there are two types of global existence results, both of which have limitations. The earliest global existence result in two dimensions was the thin domain result of Sell and Taboada [Reference Sell and Taboada33]; this was then followed by the other thin-domain results [Reference Benachour, Kukavica, Rusin and Ziane8, Reference Kukavica and Massatt23, Reference Molinet29]. Other than these, D.M.A. and Mazzucato have demonstrated global existence of small solutions for the two-dimensional Kuramoto–Sivashinsky equation for certain domain sizes (i.e., placing certain conditions on L 1 and L 2), but without the anisotropy inherent in the thin-domain results [Reference Ambrose and Mazzucato2, Reference Ambrose and Mazzucato3].
Other global results rely upon modifying either the linear or nonlinear parts of (1.1). For instance, by no longer considering fourth-order linear terms, a maximum principle may be introduced, leading to global existence of solutions [Reference Larios and Yamazaki25, Reference Molinet28]. Changing the power in the nonlinear term leads to global existence or singularity formation, depending on the power, as demonstrated in [Reference Bellout, Benachour and Titi7]. Global existence also follows from the introduction of appropriate transport terms, as shown in [Reference Coti Zelati, Dolce, Feng and Mazzucato14, Reference Feng and Mazzucato16].
Grujić and Kukavica demonstrated existence of solutions for the Kuramoto–Sivashinsky equation in one dimension, with $\partial_x \phi_0 \in L^{\infty},$ and also demonstrated analyticity of the solutions at positive times [Reference Grujić and Kukavica20]. Biswas and Swanson considered the Kuramoto–Sivashinsky equation in general dimension. Their results include improving the assumption made by Grujić and Kukavica, in dimension one, on the regularity of the data (by one); Biswas and Swanson also studied higher regularity through estimates of Gevrey norms [Reference Biswas and Swanson9].
This research is naturally related to work on the Navier–Stokes equations, for which there have been many studies of existence of solutions starting from low-regularity data. The optimal result in critical spaces is due to Koch and Tataru for data in BMO −1 [Reference Koch and Tataru22]. The present work draws more from other studies, such as by Cannone and Karch for data in PM 2 and by Lei and Lin for data in X −1 [Reference Cannone and Karch13, Reference Lei and Lin26] (see § 2 and § 4 in the present work for the definition of these spaces). In [Reference Bae5], Bae proved a version of the Lei–Lin result using a two-norm approach, which also gives analyticity of the solution at positive times and draws upon earlier work [Reference Bae, Biswas and Tadmor6]. The authors of the present article adapted the work of [Reference Bae5] to the spatially periodic case, finding an improved estimate for the radius of analyticity [Reference Ambrose, Lopes Filho and Nussenzveig Lopes1].
The primary contribution of the present work is to weaken the assumed regularity of the initial data as compared to prior works on existence of solutions for the Kuramoto–Sivashinsky equation (1.1). Our primary motivation is to examine how the two-norm approach may be used to improve regularity requirements and analyticity estimates beyond the Navier–Stokes system.
D.M.A. and Mazzucato proved the existence of solutions for the two-dimensional Kuramoto–Sivashinsky equations in the case of small domain sizes with data that have one derivative in the Wiener algebra or one derivative in L 2 [Reference Ambrose and Mazzucato2, Reference Ambrose and Mazzucato3]. Subsequently, Coti Zelati, Dolce, Feng, and Mazzucato treated situations (for an equation with added advection) with data in L 2 [Reference Coti Zelati, Dolce, Feng and Mazzucato14]; Feng and Mazzucato also treated a different class of advective equations, again with L 2 data, in [Reference Feng and Mazzucato16]. Biswas and Swanson treated the whole-space case rather than the spatially periodic case and took data such that the Fourier transform is in an L p space [Reference Biswas and Swanson9], with p ≠ 1 and $p\neq\infty$ (we treat the complementary cases of periodic data with Fourier coefficients in $\ell^{1}$ or $\ell^{\infty}$). Existence of solutions for the Kuramoto–Sivashinsky equation with pseudomeasure data was treated by Miao and Yuan, but only in non-physical spatial dimensions, specificially $n=4,$ $n=5,$ and n = 6 [Reference Miao and Yuan27]. In the present work, we deal with the physically relevant spatial dimensions n = 1 and n = 2. It is notable that our one-dimensional existence theorem allows initial data with Fourier coefficients, which grow as the Fourier variable, $k,$ goes to infinity. The Kuramoto–Sivashinsky equation in dimension n = 3 with pseudomeasure-type data remains unaddressed; while this is not a physical case when regarding the equation as a model for the motion of flame fronts, it may be of physical interest in other areas of application (such as pattern formation in reaction-diffusion equations). We discuss in remark 4.3 how our result unfortunately does not readily extend to the case n = 3.
We prove that our solutions are global in time in the case that the linearized problem has no growing Fourier modes. This amounts to an assumption of smallness of the periodic cell that comprises the spatial domain. In the general case of larger period cells, our results are valid up until a finite time. This is consistent with the lack of general global existence theory for the Kuramoto–Sivashinsky equation in dimension two and higher. In addition to proving that solutions exist, we also prove that they are analytic at positive times, following the approach of Bae [Reference Bae5], which the authors also used previously for the Navier–Stokes equations [Reference Ambrose, Lopes Filho and Nussenzveig Lopes1].
The plan of the article is as follows: we establish some preliminaries in § 2, which includes introducing a number of function spaces and giving an abstract fixed point result; in § 3, we establish existence of solutions with data in a space related to the Wiener algebra; in § 4, we treat existence of solutions with data in pseudomeasure spaces; we establish the associated linear estimates in § 4.1, the nonlinear estimates in one spatial dimension in § 4.2, and the nonlinear estimates in two spatial dimensions in § 4.3; analyticity of all of these solutions at positive times is demonstrated in § 5; the main theorems are the existence theorems, i.e., theorem 3.1 at the beginning of § 3, theorem 4.4 at the beginning of § 4.2, and theorem 4.5 at the beginning of § 4.3, and the analyticity theorems, i.e., theorem 5.3 and theorem 5.4, at the end of § 5; and we close with some concluding remarks in § 6.
2. Preliminaries
We observe that the mean of ϕ does not influence the evolution of ϕ. We thus introduce the projection $\mathbb{P}$, which removes the mean of a periodic function, as follows:
We let $\psi=\mathbb{P}\phi,$ and we note that $\nabla\psi=\nabla\phi;$ we then see that ψ satisfies the equation
Recall the Fourier series of a periodic function, given in terms of its Fourier coefficients:
From this we see directly that the symbol of the partial differential operator $\partial_{x_{i}}$ is
Therefore, the symbols of the Laplacian and bi-Laplacian are
We next introduce some spaces based on the Wiener algebra, which we denote as Y m for $m\in\mathbb{R}.$ We define $\mathbb{Z}^{n}_{*}=\mathbb{Z}^{n}\setminus\{0\}.$ A periodic function, $f,$ is in Y m if the norm given by
is finite. If $m=0,$ then this space is exactly the Wiener algebra. We let T > 0 be given, with T possibly being infinite. On the space–time domain $[0,T]\times\mathbb{T}^{n},$ we also have a related function space, $\mathcal{Y}^{m}.$ The norm for this space is
In practice, we will be dealing with functions with zero mean, so it will be equivalent for us to treat the norms as
We note that the space X −1 as used in [Reference Ambrose, Lopes Filho and Nussenzveig Lopes1, Reference Bae5, Reference Lei and Lin26] is equal to our space $Y^{-1}.$
Given $m\in\mathbb{R},$ we also have a related function space on space–time, $\mathcal{X}^{m}.$ We define the space $\mathcal{X}^{m}$ according to the norm
If f has zero mean for all times, then this becomes simply
In the results to follow, we will typically take m = 2 or m = 4.
We will consider two cases in what follows. We first describe Case A. In Case A, we assume that all $L_{i} \lt 2\pi,$ and we take $T=\infty.$ Because of the size of the $L_{i},$ we have $\sigma(\Delta^{2}+\Delta)(k) \gt 0$ for all $k\in\mathbb{Z}^{n}_{*}.$ Then, we have
In Case B, we let $T\in(0,\infty)$ be given, and we assume there exists at least one $i\in\{1,\ldots,n\}$ such that $L_{i}\geq2\pi.$ Then, there exists $M_{1} \gt 0$ such that
In Case B, we make the decomposition $\mathbb{Z}^{n}_{*}=\Omega_{F}\cup\Omega_{I},$ where for all $k\in\Omega_{F},$ the symbol is non-positive, i.e., $\sigma(\Delta^{2}+\Delta)(k)\leq0.$ Then, on the complement, of course, we have for all $k\in\Omega_{I},$ $\sigma(\Delta^{2}+\Delta)(k) \gt 0.$ Of course, the set $\Omega_{F}$ is finite and $\Omega_{I}$ is infinite.
We may, of course, also consider the decomposition $\mathbb{Z}^{n}_{*}=\Omega_{F}\cup\Omega_{I}$ in Case A as well, and then we simply have $\Omega_{F}=\emptyset.$ In either case, we have $M_{2} \gt 0$ such that
We also introduce M 3 to be the maximum value of $|k|$ for $k\in\Omega_{F},$
We will rely on the following classical abstract result:
Lemma 2.1. Let ($X,$ ${\left\vert\kern-0.25ex\left\vert\kern-0.25ex\left\vert \cdot\right\vert\kern-0.25ex\right\vert\kern-0.25ex\right\vert}_{X}$) be a Banach space. Assume that $\mathcal{B}:X \times X \to X$ is a continuous bilinear operator and let η > 0 satisfy $\eta\geq \|\mathcal{B}\|_{X\times X\rightarrow X}$. Then, for any $x_0 \in X$ such that
there exists one and only one solution to the equation
Moreover, ${\left\vert\kern-0.25ex\left\vert\kern-0.25ex\left\vert x\right\vert\kern-0.25ex\right\vert\kern-0.25ex\right\vert}_{X} \leq 2{\left\vert\kern-0.25ex\left\vert\kern-0.25ex\left\vert x_0\right\vert\kern-0.25ex\right\vert\kern-0.25ex\right\vert}_{X}$.
See [Reference Cannone11, p. 37, lemma 1.2.6] and [Reference Auscher and Tchamitchian4, Reference Cannone12].
We may write the mild formulation of the Kuramoto–Sivashinsky equation (2.1) as
Here, the semigroup operator is
and the bilinear term is
The Fourier coefficients of $B(F,G)$ are
In all of the estimates we will perform, we will use only bounds from above, with respect to the frequency variable k, of $|\widehat{B(F,G)}(t,k)|$. We thus ignore, hereafter, the constants $\displaystyle{\frac{4\pi^2}{L_{i}^2}}$, absorbing them into a positive constant C, which is then normalized to 1.
3. Existence of solutions with data in Y −1
In this section, we will prove the following theorem, giving existence of solutions with initial data taken from the space $Y^{-1}.$
Theorem 3.1 Let T > 0 be given. (If the conditions of Case A hold, then T may be taken to be $T=\infty.$) Let $n\geq1.$ There exists ɛ > 0 such that for any ϕ 0 with $\mathbb{P}\phi_{0}\in Y^{-1},$ if $\|\mathbb{P}\phi_{0}\|_{Y^{-1}} \lt \varepsilon,$ then there exists a unique ϕ with $\mathbb{P}\phi\in\mathcal{Y}^{-1}\cap \mathcal{X}^{3}$ such that ϕ is a mild solution to the initial value problem (1.1) and (1.2).
Proof. To use lemma 2.1, we need to establish the bilinear estimate, and also that $x_{0}=S\mathbb{P}\phi_{0}\in\mathcal{Y}^{-1}\cap\mathcal{X}^{3}.$ Note that the lemma gives both existence and uniqueness of the solution.
For the semigroup, we let $\psi_{0}\in Y^{-1}$ be given, and we must show $S\psi_{0}\in\mathcal{Y}^{-1}\cap\mathcal{X}^{3}.$ We begin by computing the norm in $\mathcal{Y}^{-1}:$
We next compute the norm in $\mathcal{X}^{3}:$
In Case A, the supremum is finite because we may neglect the exponential and $\sigma(\Delta^{2}+\Delta)(k) \gt M_{2}|k|^{4}.$ In Case B, we may take the supremum separately over the sets $\Omega_{F}$ and $\Omega_{I},$ and the reasoning from Case A applies to the supremum over $\Omega_{I}.$ For the supremum over $\Omega_{F},$ we find that it is finite because k is in a bounded set and T is finite. In either case, we have concluded that there exists C > 0 such that $S\psi_{0}\in\mathcal{X}^{3}$ and
This completes the proof of the needed semigroup properties.
We next need to compute $\|B(F,G)\|_{\mathcal{Y}^{-1}}$ and $\|B(F,G)\|_{\mathcal{X}^{3}}.$ We begin to compute the norm in $\mathcal{Y}^{-1}:$
We substitute from (2.10) and make some elementary bounds, arriving at
Since we have
we find
We then multiply and divide both terms by $|k-j||j|,$ arriving at
We note the elementary bounds
This, then, immediately yields the bound
We next will consider the higher norm of $B(F,G),$ attempting to bound $\|B(F,G)\|_{\mathcal{X}^{3}}.$ To begin, we have
We use the triangle inequality and (3.2), finding
We will only include the details for the estimate of $A_{1},$ as the estimate of A 2 is exactly the same. We decompose A 1 further, using the decomposition $\mathbb{Z}^{n}_{*}=\Omega_{F}\cup\Omega_{I}.$ We have
We use the definitions of M 1, (2.4), and M 3, (2.6), to immediately bound A 3 as
We then multiply and divide by $|k-j||j|$ and rearrange, finding
Here we have used the definition of M 3 to make the elementary bound
We next bound $A_{4}.$ We exchange the order of integration and compute the integral with respect to $t,$ finding
Since $\sigma(\Delta^{2}+\Delta)(k)$ is positive for $k\in\Omega_{I},$ we may neglect the exponential in the numerator, and use the definition of $M_{2},$ finding
We then multiply and divide by $|k-j||j|$ and again use (3.3), finding
We have concluded the bound
and we have by symmetry the corresponding estimate for $A_{2},$ namely
These bounds immediately imply the desired conclusion, which is
(It is understood that if the set $\Omega_{F}$ is empty, then we may take $T=\infty,$ and in this case, the combination $M_{1}T$ is understood as $M_{1}T=0.$)
4. Existence of solutions with data in pseudomeasure spaces
We prove existence theorems for the Kuramoto–Sivashinsky equation with pseudomeasure data in dimensions n = 1 and n = 2. We first define the pseudomeasure spaces.
For any $m\in\mathbb{R},$ we define the sets PM m and $\mathcal{PM}^{m}$ by their norms,
As before, we will hereafter assume that all functions considered have zero mean and will therefore only need to use these norms on spaces of functions with zero mean. We then get the simpler expressions
In § 4.1, we give the linear estimates, which are relevant for both dimensions n = 1 and n = 2. In § 4.2, we state the existence theorem in dimension n = 1 and demonstrate the needed bilinear estimates. Then, in § 4.3, we state the existence theorem for dimension n = 2 and again give the needed bilinear estimates.
4.1. Linear estimates
We give the linear estimates for pseudomeasure data in the following two lemmas.
Lemma 4.1. For any $m\in\mathbb{R},$ the semigroup operator S satisfies $S:PM^{m}\rightarrow\mathcal{PM}^{m},$ with the estimate
Proof. We begin to estimate $\|S\psi_{0}\|_{\mathcal{PM}^{m}}.$ We have
Using the definition of $M_{1},$ then, thisestimate becomes the following:
We next estimate $\|S\psi_{0}\|_{\mathcal{X}^{m_{2}}},$ with $\psi_{0}\in PM^{m_{1}}.$
Lemma 4.2. Let m 1 and m 2 be real numbers satisfying $m_{2}-m_{1}-4 \lt -n.$ There exists K > 0 such that for any $\psi_{0}\in PM^{m-1},$ we have $S\psi_{0}\in\mathcal{X}^{m_{2}}$ with the estimate
Proof. We begin by writing out the norm of $S\psi_{0}$ in the space $\mathcal{X}^{m_{2}}:$
We multiply and divide the integrand by $|k|^{m_{1}},$ and we take out the norm of $\psi_{0},$ finding
We then decompose $\mathbb{Z}^{n}_{*}$ into $\Omega_{F}\cup\Omega_{I},$ and we estimate the portion over the finite set $\Omega_{F}.$ This yields
(If we are in Case A, then since $\Omega_{F}=\emptyset,$ then in this case with $T=\infty$, we may take $C(\Omega_{F})=0,$ with the understanding that this would mean that $M_{1}C(\Omega_{F})T=0$.) Then, we evaluate the remaining integral, finding
We may then neglect the exponential and use (2.5) to find
The series on the right-hand side converges, so we have concluded that there exists K > 0 such that
Remark 4.3. The quantity $m_{2}-m_{1}$ describes how many derivatives the solution gains at positive times compared to the data. Of course, one may expect, in L 1-based spaces, to gain four derivatives from a fourth-order parabolic evolution; at the same time, because our nonlinearity only contains first derivatives, we do not require this full gain of four derivatives. Our requirement $m_{2}-m_{1}-4 \lt -n$ implies that in one space dimension, we may take $m_{2}-m_{1} \lt 3,$ and in two space dimensions, we may take $m_{2}-m_{1} \lt 2.$ In these cases, this is sufficient gain of regularity to establish the bilinear estimates. Unfortunately, in three space dimensions, we find no result, as when n = 3, we have $m_{2}-m_{1} \lt 1,$ and this is less than the gain of one full derivative which we need to estimate the nonlinearity by the present method.
4.2. Existence of solutions with pseudomeasure data with n = 1
In one space dimension, we can find the existence of solutions with PM −p data, for any $p\in(0,1/2).$ With the parabolic gain of regularity, we will also have that the solutions are in $\mathcal{X}^{2+p};$ this gain of $2+2p$ derivatives is less than the four full derivatives, which might be possible, but this is sufficient gain to deal with the nonlinearity. Note that these choices satisfy the constraints as discussed in remark 4.3. Specifically, with $m_{1}=-p$ and $m_{2}=2+p,$ and with $p \lt 1/2,$ we have $m_{2}-m_{1}=2+2p \lt 3,$ as desired.
Theorem 4.4 Let $p\in(0,1/2)$ and T > 0 be given. (If the conditions of Case A hold, then T may be taken to be $T=\infty.$) Let n = 1. There exists ɛ > 0 such that for any ϕ 0 and $\mathbb{P}\phi_{0}\in PM^{-p},$ if $\|\mathbb{P}\phi_{0}\|_{PM^{-p}} \lt \varepsilon,$ then there exists a unique ϕ with $\mathbb{P}\phi\in\mathcal{PM}^{-p}\cap \mathcal{X}^{2+p}$ such that ϕ is a mild solution to the initial value problem (1.1) and (1.2).
Proof. We will again use lemma 2.1; recall that this implies both existence and uniqueness of solutions. To apply lemma 2.1, we need to conclude that $x_{0}=S\mathbb{P}\phi_{0}$ is in the space $\mathcal{PM}^{-p}\cap\mathcal{X}^{2+p}.$ This follows from lemmas 4.1 and 4.2 with $m_{1}=-p$ and $m_{2}=2+p,$ since these parameters satisfy the condition $m_{2}-m_{1}-4 \lt -n.$
We begin with the estimate of $B(F,G)$ in $\mathcal{PM}^{-p}.$ Using the definition of $\mathcal{PM}^{-p}$ and the triangle inequality, we have
We then use Young’s inequality on $j(k-j)$ and bound the exponentials by M 1 (recall the definition of M 1 in (2.4)), finding
We multiply and divide by the appropriate powers of $|j|$ and $|k-j|:$
For the first term on the right-hand side of (4.2), we take the supremum with respect to s and k for the factor $\frac{|\hat{G}(s,k-j)|}{|k-j|^{p}}$ in the integrand. For the second term on the right-hand side, we take the supremum with respect to s and j for the factor $\frac{|\hat{F}(s,j)|}{|j|^{p}}$ in the integrand. Also, we again use (3.3). These considerations lead to the following bound:
We may further bound this as
We next must estimate $B(F,G)$ in the space $\mathcal{X}^{2+p}.$ We begin with the definition and use the triangle inequality:
Then, as before, we use Young’s inequality on $j(k-j),$
For the first term on the right-hand side, we find the $\mathcal{PM}^{-p}$-norm of G by multiplying and dividing by $|k-j|^{p}$ and taking a supremum, and we also multiply and divide by $|j|^{p}.$ We treat the second term on the right-hand side similarly, and we arrive at
Using (3.3) with (4.5), we have
In the second term on the right-hand side, we change the variable in the final summation, and we also change the order of integration in both terms on the right-hand side, finding
We will work now with the sum with respect to $k,$ which is the same in both of the terms on the right-hand side. We split it into the sum over $\Omega_{F}$ and the sum over $\Omega_{I}.$ Considering $k\in\Omega_{F},$ we have
(We have said that if $\Omega_{F}=\emptyset$, then we may take $T=\infty,$ and then this product is to be understood as $|\Omega_{F}|T=0.$) Considering $k\in\Omega_{I},$ we evaluate the integral, finding
Since the denominator is positive for $k\in\Omega_{I},$ we may neglect the exponential in the numerator. Then, we use the definition of $M_{2},$ finding
We, of course, have used here that $p \lt 1/2.$
Returning to (4.6), we conclude with the bound
4.3. Existence of solutions with pseudomeasure data with n = 2
We again let $p\in(0,1/2)$ be given. In the case of two space dimensions, we will be taking data in $PM^{1-p}$ and finding solutions in $\mathcal{PM}^{1-p}\cap\mathcal{X}^{2+p}.$ As regards remark 4.3, this means that we have $m_{1}=1-p$ and $m_{2}=2+p,$ so that $m_{2}-m_{1}=1+2p \lt 2,$ as required.
Theorem 4.5 Let $p\in(0,1/2)$ and T > 0 be given. (If the conditions of Case A hold, then T may be taken to be $T=\infty.$) Let n = 2. There exists ɛ > 0 such that for any ϕ 0 with $\mathbb{P}\phi_{0}\in PM^{1-p},$ if $\|\mathbb{P}\phi_{0}\|_{PM^{1-p}} \lt \varepsilon,$ then there exists a unique ϕ with $\mathbb{P}\phi\in\mathcal{PM}^{1-p}\cap \mathcal{X}^{2+p}$ such that ϕ is a mild solution to the initial value problem (1.1) and (1.2).
Proof. To apply lemma 2.1, we need to conclude that $x_{0}=S\mathbb{P}\phi_{0}$ is in the space $\mathcal{PM}^{1-p}\cap\mathcal{X}^{2+p}.$ This follows from lemmas 4.1 and 4.2 with $m_{1}=1-p$ and $m_{2}=2+p,$ since these parameters satisfy the condition $m_{2}-m_{1}-4 \lt -n.$
We estimate $\|B(F,G)\|_{\mathcal{PM}^{1-p}}.$ From the definition of the norm and $B(F,G),$ we have
We use the triangle inequality and the definition of the constant M 1 to find
Bounding $|k|$ as $|k|\leq |k-j|+|j|,$ this becomes
We then proceed as in the proof of theorem 4.4, i.e., we adjust the factors of $|k-j|$ and $|j|$ and use (3.3) as appropriate, until we are able to conclude that
Next we bound $B(F,G)$ in $\mathcal{X}^{2+p}.$ From the definition of the norm and of $B(F,G),$ and again using $|k|\leq |k-j|+|j|$ for just one factor of $|k|,$ we have
We then follow the corresponding steps of the proof of theorem 4.4, including adjusting factors of $|k-j|$ and $|j|$ and using (3.3) as appropriate, until we reach the conclusion
Note that the sum on the right-hand side converges because $p \lt 1/2.$
Remark 4.6. We note that the proofs of theorems 4.4 and 4.5 overall look quite similar, but there is a subtle difference. In the proof of theorem 4.4, in both (4.1) and (4.4), we have bounded $|k-j||j|$ by $\frac{1}{2}(|k-j|^{2}+|j|^{2}),$ while in the proof of theorem 4.5, we have not done so in the corresponding places. In the proof of theorem 4.5, we then take a factor of $|k|$ and bound it as $|k|\leq|k-j|+|j|.$ This first set of manipulations produces a different result from the second set of manipulations. The first set of manipulations requires a greater gain of regularity and allows lower initial regularity than the second set of manipulations. That is, if we were to proceed as in the second way for the n = 1 theorem, then the proof would not work for PM −p data. If we were to proceed in the first way for the n = 2 theorem, then the proof would not work because, for n = 2, we only gain $1+2p \lt 2$ derivatives as compared to the data rather than the $2+2p$ derivatives we gain in the n = 1 case. Recall that this limit on the gain of regularity is determined from our linear theory and is therefore a constraint on how we may conduct the nonlinear estimates.
5. Analyticity
In this section we will show that the solutions produced earlier are analytic within their time of existence, if needed by further restricting the size of the initial data.
Given initial data ψ 0, we recall the mild formulation of the Kuramoto–Sivashinsky equation (2.7):
where the semigroup S was introduced in (2.8) and the bilinear term $B=B(F,G)$ was given in (2.9).
Our approach to establish analyticity follows the one used by Bae in [Reference Bae5], in which one revisits the existence proofs but for an exponentially weighted modification of ψ. More precisely, let $g=g(t)$ be a given function and consider
where $|D|= \sqrt{-\Delta}$. Then, V should satisfy the equation
with $V_0 = \psi_0$. Existence of a solution to this equation for suitable g and sufficiently small V 0 in certain function spaces then implies analyticity of ψ, as will be made precise at the end of this section. The radius of analyticity is bounded from below by g(t).
We rewrite (5.2), separating the linear from the nonlinear term, as
with
and
We will prove existence of a solution V to (5.2) for initial data $V_0 = \psi_0$ in Y −1, in any dimension, and in $PM^{(n-1)-p}$, in dimensions n = 1 and n = 2, with $0 \lt p \lt 1/2$. As before, we use lemma 2.1, so we require the following bounds:
as well as
The Fourier coefficients of $(\mathcal{L} V_0)(t,\cdot)$ are given by
The Fourier coefficients of the nonlinear term are
In what follows, we will consider two kinds of temporal weights:
In order to estimate the linear term $\mathcal{L}$, we will make use of the following technical lemma:
Lemma 5.1. Let $k \in \Omega_I$. Then, if M 2 is as in (2.5), it holds that:
(i) if $g(t) = bt$, with $b \lt \frac{M_2}{2}$, then
\begin{equation*}g(t) |k| - t\sigma(\Delta^2+\Delta)(k) \leq - \frac{M_2 t}{2}|k|^4;\end{equation*}(ii) if $g(t) = a\sqrt[4]{t}$, then there exists $C=C(a) \gt 0$ such that
\begin{equation*}g(t) |k| - t\sigma(\Delta^2+\Delta)(k) \leq C - \frac{M_2 t}{2}|k|^4\end{equation*}for all $t \geq 0$.
Proof. Recall the definition of M 2, from (2.5), such that
Let us first consider the case $g(t) = bt$ with $b \lt \frac{M_2}{2}$. Then, clearly, if k is such that $|k| \geq 1$, it follows that $bt|k| - \frac{M_2}{2} t |k|^4 \leq 0$. It then follows easily that, if $|k| \geq 1$,
as desired. This establishes item (i).
Next, consider the case $g(t) = a\sqrt[4]{t}$. Then, of course, we have
Consider the function $f=f(z)=az - \frac{M_2}{2}z^4$. This function is globally bounded from above. Let
Noting that $a\sqrt[4]{t}|k| - \frac{M_2}{2} t |k|^4 = f(\sqrt[4]{t}|k|)$ we obtain item (ii).
In view of lemma 5.1, all the estimates for $\mathcal{L} V_0 \equiv \mathcal{L}\psi_0$, (5.4), (5.5), (5.6), and (5.7), can be reduced to the corresponding estimates already obtained for $S\psi_0$, namely (3.1) and those obtained in lemmas 4.1 and 4.2.
Next, we prove another technical lemma, which will be used for the nonlinear term $\mathcal{B} (U,W)$.
Lemma 5.2. Let $k \in \Omega_I$. Then, if M 2 is as in (2.5), it holds that:
(i) if $g(t) = bt$, with $b \lt \frac{M_2}{2}$, then it follows that
\begin{equation*}(g(t) - g(s))|k| - (t-s)\sigma(\Delta^2+\Delta)(k) \leq - \frac{M_2 (t-s)}{2}|k|^4\end{equation*}for all t s such that $0\leq s \leq t$ and all $k \in \Omega_I$;
(ii) if $g(t) = a\sqrt[4]{t}$, then there exists $C=C(a) \gt 0$ such that
\begin{equation*}(g(t) - g(s)) |k| - (t-s)\sigma(\Delta^2+\Delta)(k) \leq C - \frac{M_2 (t-s)}{2}|k|^4\end{equation*}for all t s such that $0\leq s \leq t$ and all $k \in \Omega_I$.
Proof. Let us begin with item (i), the case $g(t)=bt$, $b \lt M_2/2$. In view of lemma 5.1 item (i), this is trivial since
and $t-s \geq 0$.
Next consider $g(t) = a\sqrt[4]{t}$. We first note that
using, again, (2.5) and the fact that $t-s\geq 0$. Next, we observe that
where f was introduced in the proof of lemma 5.1. There are two possibilities: either $f(a\sqrt[4]{s})\geq 0$, in which case we may ignore this term and use the boundedness from above of f to obtain (ii), or $f(a\sqrt[4]{s}) \lt 0$. Let us assume the latter and note that f has only two real roots, namely 0 and $\sqrt[3]{\frac{2a}{M_2}} \gt 0$. In addition, f restricted to the positive real axis is only negative on the interval $\left(\sqrt[3]{\frac{2a}{M_2}}, + \infty\right)$, on which it is also decreasing. Therefore, since $s \leq t$, we have $f(a\sqrt[4]{t}) \leq f(a\sqrt[4]{s})$, from which (ii) follows immediately.
Now we rewrite (5.13) in a more convenient form and estimate:
where we used the triangle inequality to estimate $|k|-|k-j|-|j|\leq 0$.
In view of lemma 5.2, it is easy to see that the estimates on the term (5.18), namely (5.8), (5.9), (5.10) and (5.11), can be reduced to the corresponding ones for $B(F,G)$, (3.4), (3.5), (4.3), (4.7), (4.8) and (4.9), established in the previous sections.
We now comment on how these results imply analyticity of solutions. By the periodic analogue of theorem IX.13 of [Reference Reed and Simon32], a function is analytic with radius of analyticity at least ρ if its Fourier series decays like $e^{-\tilde{\rho}|k|}$ for all $\tilde{\rho} \lt \rho.$ With a solution $V\in\mathcal{Y}^{-1},$ then at each time t, we have, for any $\varepsilon \gt 0,$ the existence of c > 0 such that
We see from this that $\hat{\psi}(t,\cdot)$ decays like $e^{-\tilde{\rho}|k|}$ for any $\tilde{\rho} \lt g(t),$ and thus ψ is analytic with radius of analyticity at least $g(t).$ Similarly, for solutions with $V\in\mathcal{PM}^{(n-1)-p},$ the solution of Kuramoto–Sivashinsky, $\psi,$ is again analytic with radius of analyticity at least $g(t).$
By virtue of these considerations, we have established the following results.
Theorem 5.3 Let T > 0 be given. (If the conditions of Case A hold, then T may be taken to be $T=\infty.$) Let $n\geq1.$ There exists ɛ > 0 such that for any ϕ 0 with $\mathbb{P}\phi_{0}\in Y^{-1},$ if $\|\mathbb{P}\phi_{0}\|_{Y^{-1}} \lt \varepsilon,$ then there exists a unique ϕ with $\mathbb{P}\phi\in\mathcal{Y}^{-1}\cap \mathcal{X}^{3}$ such that ϕ is an analytic mild solution to the initial value problem (1.1) and (1.2) with radius of analyticity at least $R(t) = \max\{a\sqrt[4]{t},bt\}$, with $b \lt M_2$ and a > 0.
Theorem 5.4 Let $n \in \{1, \, 2\}$. Let $p\in(0,1/2)$ and T > 0 be given. (If the conditions of Case A hold, then T may be taken to be $T=\infty.$) There exists ɛ > 0 such that for any ϕ 0 with $\mathbb{P}\phi_{0}\in PM^{(n-1)-p},$ if $\|\mathbb{P}\phi_{0}\|_{PM^{(n-1)-p}} \lt \varepsilon,$ then there exists a unique ϕ with $\mathbb{P}\phi\in\mathcal{PM}^{(n-1)-p}\cap \mathcal{X}^{2+p}$ such that ϕ is an analytic mild solution to the initial value problem (1.1) and (1.2) with radius of analyticity at least $R(t) = \max\{a\sqrt[4]{t},bt\}$, with $b \lt M_2$ and a > 0.
6. Concluding remarks
We close now with a few remarks.
First, we comment on our bound for the radius of analyticity. We have shown that our solutions have radius of analyticity that grows at least like $t^{1/4}$ and also at least like $t.$ Of course, the rate $t^{1/4}$ is faster for times near zero, and the rate t is faster for large times. A fractional-power rate has previously been observed for the Navier–Stokes equations (where the rate is $t^{1/2}$) and for the Kuramoto–Sivashinsky equation (where the rate is $t^{1/4}$) for solutions on $\mathbb{R}^{n}$ [Reference Bae5, Reference Foias and Temam18, Reference Grujić and Kukavica20]. For spatially periodic problems, rates like t have been observed previously for the Navier–Stokes equations [Reference Foias and Temam18] or for more general parabolic equations [Reference Ferrari and Titi17]. We have previously observed for the Navier–Stokes equations that in the periodic case, one gets the improvement that both of these rates hold [Reference Ambrose, Lopes Filho and Nussenzveig Lopes1]. The present work shows this improvement in the periodic case of the Kuramoto–Sivashinsky equation. The radius of analyticity of solutions is relevant for the convergence rate of numerical simulations [Reference Doelman and Titi15].
Let $n \in \{1,2\}$. Then, we note that the two function spaces for the initial data we have considered in this work, Y −1 and $PM^{(n-1)-p}$, with $p \lt 1/2$, are not comparable. Consider, for instance, f such that
Then, $f\in PM^{(n-1)-p},$ but $f\notin Y^{-1}.$ On the other hand, let f be such that
Then, $f\in Y^{-1}$, but $f\notin PM^{-1/2},$ and thus $f\notin PM^{(n-1)-p}.$
In the introduction, we mentioned the Navier–Stokes results of Koch and Tataru [Reference Koch and Tataru22], Cannone and Karch [Reference Cannone and Karch13], and Lei and Lin [Reference Lei and Lin26] as works proving existence of solutions for the Navier–Stokes equations with low-regularity data. It should be noted that the function spaces considered in the aforementioned works, $BMO^{-1},$ $PM^{2},$ and $X^{-1},$ respectively, are all critical spaces for the Navier–Stokes equations.
If we discard the unstable Laplacian and consider (1.1) in full space $\mathbb{R}^n$, then it is easy to see that this modified Kuramoto–Sivashinsky equation is invariant under the scaling
Thus, among the hierarchy of spaces considered in this work, the spaces Y −2 and $PM^{n-2}$ are critical spaces, i.e., whose norms are invariant under this scaling. In the present work, we lower the regularity requirements for existence theory for the Kuramoto–Sivashinsky equation as compared to the prior literature and have proved existence of solutions in spaces of negative index, namely Y −1 and $PM^{(n-1)-p}$, for any $0 \lt p \lt 1/2$, but these spaces are not critical. Therefore, there remains work to be done to continue lowering the regularity threshold for the initial data.
We also mentioned in the introduction that solutions of the Kuramoto–Sivashinsky equation have been proved to be global in one spatial dimension, when starting from H 1 data. In the present work, we have shown existence of solutions with rough data, but only until a short time (unless the spatial domain $[0,L]$ satisfies $L \lt 2\pi$). But we have shown that the solutions are analytic at positive times, and thus, the solutions instantaneously become H 1 solutions, which could then be continued for all time. So, our one-dimensional solutions are in fact global. However, the present method would not extend on its own to demonstrate this. The radius of analyticity that we prove grows in time, like both $t^{1/4}$ and $t.$ This growth of the radius for all time is possible in the small-domain case (again, $L \lt 2\pi$), but in the presence of linearly growing modes ($L \gt 2\pi$), one would not expect this. Instead, the solution, in some cases, tends towards coherent structures such as travelling waves or time-periodic waves, and these attracting solutions tend to have finite radius of analyticity. The long-time behaviour of the radius of analyticity for the initial value problem, then, is to tend towards this value of the radius of analyticity rather than to tend towards infinity. This can be seen from computational work such as [Reference Kalogirou, Keaveny and Papageorgiou21, Reference Papageorgiou, Smyrlis and Tomlin31]. Understanding in more detail the time evolution of the radius of analyticity of solutions of the Kuramoto-Sivashinsky problem will be a subject of future work.
Acknowledgements
M.C.L.F. and H.J.N.L. gratefully acknowledge the hospitality of the Department of Mathematics at Drexel University, where part of this research was done.
D.M.A. gratefully acknowledges support from the National Science Foundation through grants DMS-1907684 and DMS-2307638. M.C.L.F. was partially supported by CNPq, through Grant # 304990/2022-1, and FAPERJ, through Grant # E-26/201.209/2021. H.J.N.L. acknowledges the support of CNPq, through Grant # 305309/2022-6, and of FAPERJ, through Grant # E-26/201.027/2022.
Statements and Declarations
Data availability
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.
Competing Interests
The authors declare that they have no competing interests.