Hostname: page-component-586b7cd67f-r5fsc Total loading time: 0 Render date: 2024-11-21T17:26:43.923Z Has data issue: false hasContentIssue false

Existence and analyticity of solutions of the Kuramoto–Sivashinsky equation with singular data

Published online by Cambridge University Press:  20 November 2024

David M. Ambrose
Affiliation:
Department of Mathematics, Drexel University, Philadelphia, PA 19104, USA ([email protected]) (corresponding author)
Milton C. Lopes Filho
Affiliation:
Departamento de Matemática Aplicada, Instituto de Matematica, Universidade Federal do Rio de Janeiro, Caixa Postal 68530, Rio de Janeiro, RJ 21941-909, Brazil ([email protected], [email protected])
Helena J. Nussenzveig Lopes
Affiliation:
Departamento de Matemática, Instituto de Matemática, Universidade Federal do Rio de Janeiro, Caixa Postal 68530, Rio de Janeiro, RJ, Brazil ([email protected])
Rights & Permissions [Opens in a new window]

Abstract

We prove the existence of solutions to the Kuramoto–Sivashinsky equation with low regularity data in function spaces based on the Wiener algebra and in pseudomeasure spaces. In any spatial dimension, we allow the data to have its antiderivative in the Wiener algebra. In one spatial dimension, we also allow data that are in a pseudomeasure space of negative order. In two spatial dimensions, we also allow data that are in a pseudomeasure space one derivative more regular than in the one-dimensional case. In the course of carrying out the existence arguments, we show a parabolic gain of regularity of the solutions as compared to the data. Subsequently, we show that the solutions are in fact analytic at any positive time in the interval of existence.

Type
Research Article
Copyright
© The Author(s), 2024. Published by Cambridge University Press on behalf of The Royal Society of Edinburgh

1. Introduction

The scalar form of the Kuramoto–Sivashinsky equation is

(1.1)\begin{equation} \phi_{t}+\frac{1}{2}|\nabla\phi|^{2}+\Delta^{2}\phi+\Delta\phi=0. \end{equation}

This is taken with initial data

(1.2)\begin{equation} \phi(\cdot,0)=\phi_{0}. \end{equation}

The spatial domain we consider is the n-dimensional torus $\mathbb{T}^{n},$ which is given by

\begin{equation*}\nonumber \mathbb{T}^{n}=\prod_{i=1}^{n}\left[0,L_{i}\right] \end{equation*}

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:

\begin{equation*}\nonumber \mathbb{P}f=f-\frac{1}{L_{1}\cdots L_{n}}\int_{\mathbb{T}^{n}}f(x)\ dx. \end{equation*}

We let $\psi=\mathbb{P}\phi,$ and we note that $\nabla\psi=\nabla\phi;$ we then see that ψ satisfies the equation

(2.1)\begin{equation} \psi_{t}+\frac{1}{2}\mathbb{P}|\nabla\psi|^{2}+\Delta^{2}\psi+\Delta\psi=0. \end{equation}

Recall the Fourier series of a periodic function, given in terms of its Fourier coefficients:

\begin{equation*}\nonumber f(x) \sim \sum_{k\in\mathbb{Z}^{n}} \hat{f}(k)e^{2\pi i k_{1}x_{1}/L_{1}}\cdots e^{2\pi i k_{n}x_{n}/L_{n}}. \end{equation*}

From this we see directly that the symbol of the partial differential operator $\partial_{x_{i}}$ is

\begin{equation*}\nonumber \sigma(\partial_{x_{j}})=\frac{2\pi i}{L_{j}}k_{j}. \end{equation*}

Therefore, the symbols of the Laplacian and bi-Laplacian are

\begin{equation*}\nonumber \sigma(\Delta)=-\sum_{j=1}^{n}\left(\frac{2\pi}{L_{j}}\right)^{2}k_{j}^{2}, \end{equation*}
\begin{equation*}\nonumber \sigma(\Delta^{2})=\left(\sum_{j=1}^{n}\left(\frac{2\pi}{L_{j}}\right)^{2}k_{j}^{2}\right)^{2}. \end{equation*}

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

\begin{equation*}\nonumber \|f\|_{Y^{m}}=|\hat{f}(0)|+\sum_{k\in\mathbb{Z}^{n}_{*}}|k|^{m}|\hat{f}(k)| \end{equation*}

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

\begin{equation*}\nonumber \|f\|_{\mathcal{Y}^{m}}=\sup_{t\in[0,T]}|\hat{f}(t,0)|+\sum_{k\in\mathbb{Z}^{n}}\sup_{t\in[0,T]}|k|^{m}|\hat{f}(t,k)|. \end{equation*}

In practice, we will be dealing with functions with zero mean, so it will be equivalent for us to treat the norms as

\begin{equation*}\nonumber \|f\|_{Y^{m}}=\sum_{k\in\mathbb{Z}^{n}\setminus\{0\}}|k|^{m}|\hat{f}(k)|, \end{equation*}
\begin{equation*}\nonumber \|f\|_{\mathcal{Y}^{m}}=\sum_{k\in\mathbb{Z}^{n}\setminus\{0\}}\sup_{t\in[0,T]}|k|^{m}|\hat{f}(t,k)|. \end{equation*}

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

(2.2)\begin{equation} \|f\|_{\mathcal{X}^{m}}=\int_{0}^{T}|\hat{f}(t,0)|\ dt + \sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}|k|^{m}|\hat{f}(t,k)| \ dt. \end{equation}

If f has zero mean for all times, then this becomes simply

(2.3)\begin{equation} \|f\|_{\mathcal{X}^{m}}= \sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}|k|^{m}|\hat{f}(t,k)| \ dt. \end{equation}

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

\begin{equation*}\nonumber \sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{n}_{*}}e^{-t\sigma(\Delta^{2}+\Delta)(k)} =\sup_{t\in[0,\infty)}\sup_{k\in\mathbb{Z}^{n}_{*}}e^{-t\sigma(\Delta^{2}+\Delta)(k)}=1. \end{equation*}

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

(2.4)\begin{equation} \sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{n}_{*}}e^{-t\sigma(\Delta^{2}+\Delta)(k)}\leq M_{1}. \end{equation}

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

(2.5)\begin{equation} \sigma(\Delta^{2}+\Delta)(k) \gt M_{2}|k|^{4} \qquad \forall k\in\Omega_{I}. \end{equation}

We also introduce M 3 to be the maximum value of $|k|$ for $k\in\Omega_{F},$

(2.6)\begin{equation} |k|\leq M_{3} \qquad \forall k\in\Omega_{F}. \end{equation}

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

\begin{equation*}4\eta {\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} \lt 1,\end{equation*}

there exists one and only one solution to the equation

\begin{equation*}x=x_0+\mathcal{B}(x,x) \qquad \text{with } {\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} \lt \frac{1}{2\eta}.\end{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

(2.7)\begin{equation} \psi=S\psi_{0}-\frac{1}{2}B(\psi,\psi). \end{equation}

Here, the semigroup operator is

(2.8)\begin{equation} S\psi_{0}=e^{-t(\Delta^{2}+\Delta)}\psi_{0}, \end{equation}

and the bilinear term is

(2.9)\begin{equation} B(F,G)=\int_{0}^{t}e^{-(t-s)(\Delta^{2}+\Delta)}\mathbb{P}(\nabla F\cdot\nabla G)\ ds. \end{equation}

The Fourier coefficients of $B(F,G)$ are

(2.10)\begin{align} \widehat{B(F,G)}(t,k)=\int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\mathcal{F}[\mathbb{P}(\nabla F\cdot\nabla G)](s,k)\ ds \end{align}
\begin{align*} = \int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{n}_{*}, j\neq k}\sum_{\ell=1}^{n} \frac{2\pi i}{L_{\ell}}(k_{\ell}-j_{\ell})\hat{F}(s,k-j)\frac{2\pi i}{L_{\ell}} j_{\ell}\hat{G}(s,j)\ ds \end{align*}
\begin{align*} = -\int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{n}_{*}, j\neq k}\sum_{i=1}^{n} \frac{4\pi^2}{L_{i}^2}(k_{i}-j_{i})\hat{F}(s,k-j)j_{i}\hat{G}(s,j)\ ds. \end{align*}

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}:$

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{Y}^{-1}}=\sum_{k\in\mathbb{Z}^{n}_{*}}\sup_{t\in[0,T]}\frac{e^{-t\sigma(\Delta^{2}+\Delta)(k)}}{|k|} |\hat{\psi}_{0}(k)| \leq M_{1}\|\psi_{0}\|_{Y^{-1}}. \end{equation*}

We next compute the norm in $\mathcal{X}^{3}:$

\begin{align*} \nonumber \|S\psi_{0}\|_{\mathcal{X}^{3}}=\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}|k|^{3} e^{-t\sigma(\Delta^{2}+\Delta)(k)}|\hat{\psi}_{0}(k)|\ dt \qquad \qquad \qquad \end{align*}
\begin{align*} \phantom{\|S\psi_{0}\|_{\mathcal{X}^{3}}}\leq\left(\sum_{k\in\mathbb{Z}^{n}_{*}}\frac{|\hat{\psi}_{0}(k)|}{|k|}\right) \sup_{k\in\mathbb{Z}^{n}_{*}}\left(\frac{|k|^{4}(1-e^{-T\sigma(\Delta^{2}+\Delta)(k)})}{\sigma(\Delta^{2}+\Delta)(k)}\right). \end{align*}

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

(3.1)\begin{equation} \|S\psi_{0}\|_{\mathcal{X}^{3} } \leq C \| \psi_0\|_{Y^{-1}}. \end{equation}

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}:$

\begin{equation*}\nonumber \|B(F,G)\|_{\mathcal{Y}^{-1}}=\sum_{k\in\mathbb{Z}^{n}_{*}}\sup_{t\in[0,T]}\frac{|\widehat{B(F,G)}(t,k)|}{|k|}. \end{equation*}

We substitute from (2.10) and make some elementary bounds, arriving at

\begin{equation*}\nonumber \|B(F,G)\|_{\mathcal{Y}^{-1}}\leq nM_{1}\sum_{k\in\mathbb{Z}^{n}_{*}} \int_{0}^{T} \sum_{j\in\mathbb{Z}^{n}_{*}}\frac{1}{|k|}|k-j||\hat{F}(s,k-j)|j||\hat{G}(s,j)|\ ds. \end{equation*}

Since we have

(3.2)\begin{equation} |k-j||j|\leq \frac{1}{2}\left(|k-j|^{2}+|j|^{2}\right), \end{equation}

we find

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{Y}^{-1}}\leq \frac{nM_{1}}{2}\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{n}_{*}} \frac{|k-j|^{2}}{|k|}|\hat{F}(s,k-j)||\hat{G}(s,j)|\ ds \end{align*}
\begin{align*} \qquad \qquad \quad +\frac{nM_{1}}{2}\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{n}_{*}} |\hat{F}(s,k-j)| \frac{|j|^{2}}{|k|}|\hat{G}(s,j)|\ ds. \end{align*}

We then multiply and divide both terms by $|k-j||j|,$ arriving at

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{Y}^{-1}}\leq \frac{nM_{1}}{2}\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{3}|\hat{F}(s,k-j)|\frac{|\hat{G}(s,j)|}{|j|}\left(\frac{|j|}{|k||k-j|}\right)\ ds \end{align*}
\begin{align*} +\frac{nM_{1}}{2}\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{n}_{*}} \frac{|\hat{F}(s,k-j)|}{|k-j|} |j|^{3}|\hat{G}(s,j)|\left(\frac{|k-j|}{|k||j|}\right)\ ds. \end{align*}

We note the elementary bounds

(3.3)\begin{equation} \frac{|j|}{|k||k-j|}\leq\frac{|k|+|k-j|}{|k||k-j|}\leq 2,\qquad \frac{|k-j|}{|k||j|}\leq\frac{|k|+|j|}{|k||j|}\leq 2. \end{equation}

This, then, immediately yields the bound

(3.4)\begin{align} \|B(F,G)\|_{\mathcal{Y}^{-1}}\leq nM_{1}\|F\|_{\mathcal{X}^{3}}\|G\|_{\mathcal{Y}^{-1}} +nM_{1}\|F\|_{\mathcal{Y}^{-1}}\|G\|_{\mathcal{X}^{3}} \end{align}
\begin{align*} \qquad \qquad \qquad \leq 2nM_{1}(\|F\|_{\mathcal{Y}^{-1}}+\|F\|_{\mathcal{X}^{3}}) (\|G\|_{\mathcal{Y}^{-1}}+\|G\|_{\mathcal{X}^{3}}). \end{align*}

We next will consider the higher norm of $B(F,G),$ attempting to bound $\|B(F,G)\|_{\mathcal{X}^{3}}.$ To begin, we have

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{X}^{3}}=\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}|k|^{3}|\widehat{B(F,G)}(k)|\ dt\end{align*}
\begin{align*} = \sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\left|\int_{0}^{t}|k|^{3}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)} \sum_{j\in\mathbb{Z}^{n}_{*}}\sum_{i=1}^{n}(k_{i}-j_{i})\hat{F}(s,k-j)j_{i}\hat{G}(s,j)\ ds\right| dt. \end{align*}

We use the triangle inequality and (3.2), finding

\begin{align*} &\|B(F,G)\|_{\mathcal{X}^{3}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align*}
\begin{align*} \leq \frac{1}{2}\sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\int_{0}^{t}|k|^{3}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{2} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ ds dt \end{align*}
\begin{align*} +\frac{1}{2} \sum_{k\in\mathbb{Z}^{n}_{*}}\int_{0}^{T}\int_{0}^{t}|k|^{3}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{n}_{*}} |\hat{F}(s,k-j)| |j|^{2} |\hat{G}(s,j)|\ ds dt \end{align*}
\begin{align*} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ =A_{1}+A_{2}.\end{align*}

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

\begin{align*} \nonumber A_{1}=\frac{1}{2}\sum_{k\in\Omega_{F}}\int_{0}^{T}\int_{0}^{t}|k|^{3}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{2} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ ds dt \end{align*}
\begin{align*} +\frac{1}{2} \sum_{k\in\Omega_{I}}\int_{0}^{T}\int_{0}^{t}|k|^{3}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{2} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ ds dt \end{align*}
\begin{align*} \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ =A_{3}+A_{4}. \end{align*}

We use the definitions of M 1, (2.4), and M 3, (2.6), to immediately bound A 3 as

\begin{equation*}\nonumber A_{3}\leq \frac{1}{2}M_{1}M_{3}^{3}T\int_{0}^{T}\sum_{k\in\Omega_{F}}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{2} |\hat{F}(s,k-j)| \left[\sup_{\tau\in[0,T]}|\hat{G}(\tau,j)|\right]\ ds. \end{equation*}

We then multiply and divide by $|k-j||j|$ and rearrange, finding

\begin{align*} \nonumber A_{3}\leq \frac{1}{2}M_{1}M_{3}^{3}(M_{3}+1) T\int_{0}^{T}\sum_{k\in\mathbb{Z}^{n}_{*}}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{3} |\hat{F}(s,k-j)| \left[\sup_{\tau\in[0,T]}\frac{|\hat{G}(\tau,j)}{|j|}\right]\ ds. \end{align*}
\begin{align*} \qquad\qquad\qquad\qquad\quad\leq \frac{1}{2}M_{1}M_{3}^{4}(M_{3}+1)T\|F\|_{\mathcal{X}^{3}}\|G\|_{\mathcal{Y}^{-1}}. \end{align*}

Here we have used the definition of M 3 to make the elementary bound

\begin{equation*}\nonumber \sup_{k\in\Omega_{F}}\sup_{j\in\mathbb{Z}^{n},\ j\neq k}\frac{|j|}{|k-j|} = \sup_{k\in\Omega_{F}}\sup_{\ell\in\mathbb{Z}^{n},\ell\neq0} \frac{|\ell+k|}{|\ell|}\leq M_{3}+1. \end{equation*}

We next bound $A_{4}.$ We exchange the order of integration and compute the integral with respect to $t,$ finding

\begin{align*} \nonumber A_{4}=\frac{1}{2}\sum_{k\in\Omega_{I}}\int_{0}^{T}\int_{s}^{T}|k|^{3}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)} \sum_{j\in\mathbb{Z}^{n}_{*}}|k-j|^{2} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ dt ds \end{align*}
\begin{align*} =\frac{1}{2}\sum_{k\in\Omega_{I}}\int_{0}^{T}\frac{|k|^{4}}{|k|} \frac{1-e^{-(T-s)\sigma(\Delta^{2}+\Delta)(k)}}{\sigma(\Delta^{2}+\Delta)(k)} \sum_{j\in\mathbb{Z}^{n}_{*}}|k-j|^{2} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ ds. \end{align*}

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

\begin{equation*}\nonumber A_{4}\leq \frac{1}{2M_{2}}\int_{0}^{T}\sum_{k\in\mathbb{Z}^{n}_{*}}\sum_{j\in\mathbb{Z}^{n}_{*}} \frac{|k-j|^{2}}{|k|} |\hat{F}(s,k-j)| |\hat{G}(s,j)| \ ds. \end{equation*}

We then multiply and divide by $|k-j||j|$ and again use (3.3), finding

\begin{align*} \nonumber A_{4}\leq\frac{1}{2M_{2}}\int_{0}^{T}\sum_{k\in\mathbb{Z}^{n}_{*}}\sum_{j\in\mathbb{Z}^{n}_{*}} |k-j|^{3} |\hat{F}(s,k-j)| \frac{|\hat{G}(s,j)|}{|j|}\left(\frac{|j|}{|k||k-j|}\right) \ ds \end{align*}
\begin{align*} \qquad\qquad\qquad\qquad\qquad\qquad\leq \frac{1}{M_{2}}\|F\|_{\mathcal{X}^{3}}\|G\|_{\mathcal{Y}^{-1}}. \end{align*}

We have concluded the bound

\begin{equation*}\nonumber A_{1}\leq\left(M_{1}M_{3}^{3}(M_{3}+1)T+\frac{1}{M_{2}}\right)\|F\|_{\mathcal{X}^{3}}\|G\|_{\mathcal{Y}^{-1}}, \end{equation*}

and we have by symmetry the corresponding estimate for $A_{2},$ namely

\begin{equation*}\nonumber A_{2}\leq\left(M_{1}M_{3}^{3}(M_{3}+1)T+\frac{1}{M_{2}}\right)\|F\|_{\mathcal{Y}^{-1}}\|G\|_{\mathcal{X}^{3}}. \end{equation*}

These bounds immediately imply the desired conclusion, which is

(3.5)\begin{align} \|B(F,G)\|_{\mathcal{X}^{3}} &\leq \left(M_{1}M_{3}^{3}(M_{3}+1)T +\frac{1}{M_{2}}\right)\nonumber\\ & \quad \times (\|F\|_{\mathcal{Y}^{-1}}+\|F\|_{\mathcal{X}^{3}}) (\|G\|_{\mathcal{Y}^{-1}}+\|G\|_{\mathcal{X}^{3}}). \end{align}

(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,

\begin{equation*}\nonumber \|f\|_{PM^{m}}=|\hat{f}(0)|+\sup_{k\in\mathbb{Z}^{n}}|k|^{m}|\hat{f}(k)|, \end{equation*}
\begin{equation*}\nonumber \|f\|_{\mathcal{PM}^{m}}=\sup_{t\in[0,T]}|\hat{f}(t,0)|+\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{n}}|k|^{m}|\hat{f}(t,k)|. \end{equation*}

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

\begin{equation*}\nonumber \|f\|_{PM^{m}}=\sup_{k\in\mathbb{Z}^{n}_{*}}|k|^{m}|\hat{f}(k)|, \end{equation*}
\begin{equation*}\nonumber \|f\|_{\mathcal{PM}^{m}}=\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{n}_{*}}|k|^{m}|\hat{f}(t,k)|. \end{equation*}

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

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{PM}^{m}}\leq M_{1}\|\psi_{0}\|_{{PM}^{m}}. \end{equation*}

Proof. We begin to estimate $\|S\psi_{0}\|_{\mathcal{PM}^{m}}.$ We have

\begin{align*} \nonumber \|S\psi_{0}\|_{\mathcal{PM}^{m}}=\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{n}_{*}}|k|^{m}e^{-t\sigma(\Delta^{2}+\Delta)(k)} |\hat{\psi}_{0}(k)| \end{align*}
\begin{align*} \qquad\qquad\quad\leq \|\psi_{0}\|_{PM^{m}}\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{n}_{*}}e^{-t\sigma(\Delta^{2}+\Delta)(k)}. \end{align*}

Using the definition of $M_{1},$ then, thisestimate becomes the following:

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{PM}^{m}}\leq M_{1}\|\psi_{0}\|_{PM^{m}}. \end{equation*}

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

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}^{m_{2}}}\leq K\|\psi_{0}\|_{PM^{m_{1}}}. \end{equation*}

Proof. We begin by writing out the norm of $S\psi_{0}$ in the space $\mathcal{X}^{m_{2}}:$

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}^{m_{2}}}=\int_{0}^{T}\sum_{k\in\mathbb{Z}^{n}_{*}}|k|^{m_{2}}e^{-t\sigma(\Delta^{2}+\Delta)(k)} |\hat{\psi}_{0}(k)|\ dt. \end{equation*}

We multiply and divide the integrand by $|k|^{m_{1}},$ and we take out the norm of $\psi_{0},$ finding

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}^{m_{2}}}\leq\|\psi_{0}\|_{PM^{m_{1}}}\int_{0}^{T}\sum_{k\in\mathbb{Z}^{n}_{*}} |k|^{m_{2}-m_{1}}e^{-t(\sigma(\Delta^{2}+\Delta)(k))}\ dt. \end{equation*}

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

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}^{m_{2}}}\leq\|\psi_{0}\|_{PM^{m_{1}}}\left(M_{1}C(\Omega_{F})T+\int_{0}^{T}\sum_{k\in\Omega_{I}} |k|^{m_{2}-m_{1}}e^{-t\sigma(\Delta^{2}+\Delta)(k)}\ dt\right). \end{equation*}

(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

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}_{m_{2}}}\leq\|\psi_{0}\|_{PM^{m_{1}}}\left(M_{1}C(\Omega_{F})T+\sum_{k\in\Omega_{I}} \frac{|k|^{m_{2}-m_{1}}(1-e^{-T(\sigma(\Delta^{2}+\Delta)(k))}}{\sigma(\Delta^{2}+\Delta)(k)}\right). \end{equation*}

We may then neglect the exponential and use (2.5) to find

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}^{m_{2}}}\leq\|\psi_{0}\|_{PM^{m_{1}}}\left(M_{1}C(\Omega_{F})T+\frac{1}{M_{2}} \sum_{k\in\Omega_{I}}|k|^{m_{2}-m_{1}-4}\right). \end{equation*}

The series on the right-hand side converges, so we have concluded that there exists K > 0 such that

\begin{equation*}\nonumber \|S\psi_{0}\|_{\mathcal{X}^{m_{2}}}\leq K\|\psi_{0}\|_{PM^{m_{1}}}. \end{equation*}

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

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{PM}^{-p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align*}
\begin{align*} =\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}_{*}}\left|\int_{0}^{t}\frac{1}{|k|^{p}} e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}_{*}}j\hat{F}(s,j)(k-j)\hat{G}(s,k-j)\ ds\right| \end{align*}
\begin{align*} \leq \sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}_{*}}\int_{0}^{t}\frac{e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}}{|k|^{p}} \sum_{j\in\mathbb{Z}_{*}} |j| |\hat{F}(s,j)| |k-j| |\hat{G}(s,k-j)|\ ds. \end{align*}

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

(4.1)\begin{align} \|B(F,G)\|_{\mathcal{PM}^{-p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align}
\begin{align*} \leq \sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}_{*}}\frac{M_{1}}{2}\int_{0}^{t} \frac{1}{|k|^{p}}\sum_{j\in\mathbb{Z}_{*}} |j|^{2} |\hat{F}(s,j)| |\hat{G}(s,k-j)|\ ds \end{align*}
\begin{align*} +\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}_{*}}\frac{M_{1}}{2}\int_{0}^{t} \frac{1}{|k|^{p}}\sum_{j\in\mathbb{Z}_{*}} |\hat{F}(s,j)| |k-j|^{2}|\hat{G}(s,k-j)|\ ds. \end{align*}

We multiply and divide by the appropriate powers of $|j|$ and $|k-j|:$

(4.2)\begin{align} \|B(F,G)\|_{\mathcal{PM}^{-p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align}
\begin{align*} \leq \sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}_{*}}\frac{M_{1}}{2}\int_{0}^{t} \sum_{j\in\mathbb{Z}_{*}}\left(\frac{|k-j|^{p}}{|k|^{p}|j|^{p}}\right) |j|^{2+p} |\hat{F}(s,j)| \left(\frac{|\hat{G}(s,k-j)|}{|k-j|^{p}}\right)\ ds \end{align*}
\begin{align*} +\sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}_{*}}\frac{M_{1}}{2}\int_{0}^{t} \sum_{j\in\mathbb{Z}_{*}}\left(\frac{|j|^{p}}{|k|^{p}|k-j|^{p}}\right) \left(\frac{|\hat{F}(s,j)|}{|j|^{p}}\right) |k-j|^{2+p}|\hat{G}(s,k-j)|\ ds. \end{align*}

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:

\begin{equation*}\nonumber \|B(F,G)\|_{\mathcal{PM}^{-p}}\\ \leq 2^{p-1}M_{1} \left(\|G\|_{\mathcal{PM}^{-p}}\|F\|_{\mathcal{X}^{2+p}} + \|F\|_{\mathcal{PM}^{-p}}\|G\|_{\mathcal{X}^{2+p}}\right). \end{equation*}

We may further bound this as

(4.3)\begin{equation} \|B(F,G)\|_{\mathcal{PM}^{-p}}\leq 2^{p-1}M_{1}(\|F\|_{\mathcal{PM}^{-p}}+\|F\|_{\mathcal{X}^{2+p}}) (\|G\|_{\mathcal{PM}^{-p}}+\|G\|_{\mathcal{X}^{2+p}}). \end{equation}

We next must estimate $B(F,G)$ in the space $\mathcal{X}^{2+p}.$ We begin with the definition and use the triangle inequality:

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{X}^{2+p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align*}
\begin{align*} \leq \sum_{k\in\mathbb{Z}_{*}}\int_{0}^{T}|k|^{2+p} \int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)} \sum_{j\in\mathbb{Z}_{*}}|j| |\hat{F}(s,j)| |k-j| |\hat{G}(s,k-j)|\ ds dt. \end{align*}

Then, as before, we use Young’s inequality on $j(k-j),$

(4.4)\begin{align} \|B(F,G)\|_{\mathcal{X}^{2+p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align}
\begin{align*} \leq\frac{1}{2}\sum_{k\in\mathbb{Z}_{*}}\int_{0}^{T}|k|^{2+p}\int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)} \sum_{j\in\mathbb{Z}_{*}} |j|^{2}|\hat{F}(s,j)| |\hat{G}(s,k-j)|\ ds dt \end{align*}
\begin{align*} +\sum_{k\in\mathbb{Z}_{*}}\frac{1}{2}\int_{0}^{T}|k|^{2+p}\int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)} \sum_{j\in\mathbb{Z}_{*}} |k-j|^{2}|\hat{F}(s,j)| |\hat{G}(s,k-j)|\ ds dt. \end{align*}

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

(4.5)\begin{align} \|B(F,G)\|_{\mathcal{X}^{2+p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align}
\begin{align*} \leq \frac{\|G\|_{\mathcal{PM}^{-p}}}{2}\sum_{k\in\mathbb{Z}_{*}}\int_{0}^{T}\int_{0}^{t} |k|^{2+p}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}_{*}}\frac{|k-j|^{p}|j|^{2+p}}{|j|^{p}} |\hat{F}(s,j)|\ ds dt \end{align*}
\begin{align*} & \qquad\qquad +\frac{\|F\|_{\mathcal{PM}^{-p}}}{2}\sum_{k\in\mathbb{Z}_{*}}\int_{0}^{T}\int_{0}^{t} |k|^{2+p}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\\ & \qquad\qquad\quad \times \sum_{j\in\mathbb{Z}_{*}}\frac{|k-j|^{2+p}|j|^{p}}{|k-j|^{p}} |\hat{G}(s,k-j)|\ ds dt. \end{align*}

Using (3.3) with (4.5), we have

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{X}^{2+p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align*}
\begin{align*} \leq 2^{p-1}\|G\|_{\mathcal{PM}^{-p}}\sum_{k\in\mathbb{Z}_{*}}\int_{0}^{T}\int_{0}^{t} |k|^{2+2p}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}_{*}}|j|^{2+p} |\hat{F}(s,j)|\ ds dt \end{align*}
\begin{align*} & \qquad\qquad+2^{p-1}\|F\|_{\mathcal{PM}^{-p}}\sum_{k\in\mathbb{Z}_{*}}\int_{0}^{T}\int_{0}^{t} |k|^{2+2p}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\\ & \qquad\qquad\quad \times \sum_{j\in\mathbb{Z}_{*}}|k-j|^{2+p} |\hat{G}(s,k-j)|\ ds dt. \end{align*}

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

(4.6)\begin{align} \|B(F,G)\|_{\mathcal{X}^{2+p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align}
\begin{align*} & \leq 2^{p-1}\|G\|_{\mathcal{PM}^{-p}}\int_{0}^{T} \left(\sum_{j\in\mathbb{Z}_{*}}|j|^{2+p}|\hat{F}(s,j)|\right)\\ & \quad \times \left(\sum_{k\in\mathbb{Z}_{*}}|k|^{2+2p}\int_{s}^{T}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\ dt\right)\ ds \end{align*}
\begin{align*} & \quad +2^{p-1}\|F\|_{\mathcal{PM}^{-p}}\int_{0}^{T} \left(\sum_{j\in\mathbb{Z}_{*}}|j|^{2+p}|\hat{G}(s,j)|\right)\\ & \quad \times \left(\sum_{k\in\mathbb{Z}_{*}}|k|^{2+2p}\int_{s}^{T}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\ dt\right)\ ds. \end{align*}

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

\begin{equation*}\nonumber \sum_{k\in\Omega_{F}}k^{2+2p}\int_{s}^{T}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\ dt\leq |\Omega_{F}|M_{1}M_{3}^{2+2p}T. \end{equation*}

(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

\begin{equation*}\nonumber \sum_{k\in\Omega_{I}}k^{2+2p}\int_{s}^{T}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\ dt =\sum_{k\in\Omega_{I}}k^{2+2p}\left(\frac{1-e^{-(T-s)\sigma(\Delta^{2}+\Delta)(k)}}{\sigma(\Delta^{2}+\Delta)(k)}\right). \end{equation*}

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

\begin{equation*}\nonumber \sum_{k\in\Omega_{I}}k^{2+2p}\int_{s}^{T}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\ dt \leq \frac{1}{M_{2}}\sum_{k\in\mathbb{Z}_{*}} \frac{1}{k^{2-2p}}=\frac{c(p)}{M_{2}} \lt \infty. \end{equation*}

We, of course, have used here that $p \lt 1/2.$

Returning to (4.6), we conclude with the bound

(4.7)\begin{align} & \|B(F,G)\|_{\mathcal{X}^{2+p}}\leq 2^{p-1}\left(|\Omega_{F}|M_{1}M_{3}^{2+2p}T+\frac{c(p)}{M_{2}}\right)\ \ \ \ \ \ \ \ \ \ \ \ \end{align}
\begin{align*} \cdot (\|F\|_{\mathcal{PM}^{-p}}+\|F\|_{\mathcal{X}^{2+p}}) (\|G\|_{\mathcal{PM}^{-p}}+\|G\|_{\mathcal{X}^{2+p}}). \end{align*}

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

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{PM}^{1-p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align*}
\begin{align*} = \sup_{t\in[0,T]}\sup_{k\in\mathbb{Z}^{2}_{*}}|k|^{1-p}\left| \int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{2}_{*}}\sum_{i=1}^{2} (k_{i}-j_{i})\hat{F}(s,k-j)j_{i}\hat{G}(s,j)\ ds\right|. \end{align*}

We use the triangle inequality and the definition of the constant M 1 to find

\begin{equation*}\nonumber \|B(F,G)\|_{\mathcal{PM}^{1-p}} \leq 2M_{1}\sup_{k\in\mathbb{Z}^{2}_{*}}|k|^{1-p}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{2}_{*}} |k-j| |\hat{F}(s,k-j)| |j| |\hat{G}(s,j)|\ ds. \end{equation*}

Bounding $|k|$ as $|k|\leq |k-j|+|j|,$ this becomes

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{PM}^{1-p}} \leq 2M_{1}\sup_{k\in\mathbb{Z}^{2}_{*}}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{2}_{*}} \frac{|k-j|^{2}|j|}{|k|^{p}} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ ds \end{align*}
\begin{align*} &\qquad\qquad\qquad\qquad+2M_{1}\sup_{k\in\mathbb{Z}^{2}_{*}}\int_{0}^{T}\sum_{j\in\mathbb{Z}^{2}_{*}} \frac{|k-j| |j|^{2}}{|k|^{p}} |\hat{F}(s,k-j)| |\hat{G}(s,j)|\ ds. \end{align*}

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

(4.8)\begin{align} \|B(F,G)\|_{\mathcal{PM}^{1-p}}\leq 2^{p+1}M_{1} (\|F\|_{\mathcal{PM}^{1-p}}+\|F\|_{\mathcal{X}^{2+p}}) (\|G\|_{\mathcal{PM}^{1-p}}+ \|G\|_{\mathcal{X}^{2+p}}). \end{align}

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

\begin{align*} \nonumber \|B(F,G)\|_{\mathcal{X}^{2+p}}\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \end{align*}
\begin{align*} \leq 2\int_{0}^{T}\sum_{k\in\mathbb{Z}^{2}_{*}}|k|^{1+p} \int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{2}_{*}} |k-j|^{2} |\hat{F}(s,k-j)| |j| |\hat{G}(s,j)| \ ds dt \end{align*}
\begin{align*} + 2\int_{0}^{T}\sum_{k\in\mathbb{Z}^{2}_{*}}|k|^{1+p} \int_{0}^{t}e^{-(t-s)\sigma(\Delta^{2}+\Delta)(k)}\sum_{j\in\mathbb{Z}^{2}_{*}} |k-j| |\hat{F}(s,k-j)| |j|^{2} |\hat{G}(s,j)| \ ds dt. \end{align*}

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

(4.9)\begin{align} \|B(F,G)\|_{\mathcal{X}^{2+p}}\leq 2^{p+1}\left(|\Omega_{F}|M_{1}M_{3}^{1+2p}T+\frac{1}{M_{2}}\sum_{k\in\mathbb{Z}^{2}_{*}}\frac{1}{|k|^{3-2p}}\right) \end{align}
\begin{align*} \cdot(\|F\|_{\mathcal{PM}^{1-p}}+\|F\|_{\mathcal{X}^{2+p}}) (\|G\|_{\mathcal{PM}^{1-p}}+\|G\|_{\mathcal{X}^{2+p}}). \end{align*}

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):

\begin{equation*}\psi = S\psi_0 - \frac{1}{2} B(\psi,\psi),\end{equation*}

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

(5.1)\begin{equation} V\equiv e^{g(t)|D|}\psi, \end{equation}

where $|D|= \sqrt{-\Delta}$. Then, V should satisfy the equation

(5.2)\begin{equation} V = e^{[g(t)|D|-t(\Delta^2 + \Delta)]}V_0 - \frac{1}{2} \int_0^t e^{[g(t)|D|-(t-s)(\Delta^2 + \Delta)]}[\mathbb{P}(|\nabla e^{-g(s)|D|}V|^2)] \, ds, \end{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

\begin{equation*}\nonumber V = \mathcal{L}V_0 -\frac{1}{2}\mathcal{B}(V,V), \end{equation*}

with

(5.3)\begin{equation} \mathcal{L} V_0 = e^{[g(t)|D|-t(\Delta^2 + \Delta)]}V_0 \end{equation}

and

\begin{equation*}\nonumber \mathcal{B}(U,W)= \int_0^t e^{[g(t)|D|-(t-s)(\Delta^2 + \Delta)]}[\mathbb{P}(\nabla e^{-g(s)|D|}U\cdot \nabla e^{-g(s)|D|}W)] \, ds. \end{equation*}

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:

(5.4)\begin{align} & \|\mathcal{L} V_0\|_{\mathcal{Y}^{-1}}\leq C \|V_0\|_{Y^{-1}} \end{align}
(5.5)\begin{align} & \|\mathcal{L} V_0\|_{\mathcal{X}^{3}}\leq C \|V_0\|_{Y^{-1}} \end{align}
(5.6)\begin{align} & \|\mathcal{L} V_0\|_{\mathcal{PM}^{(n-1)-p}}\leq C \|V_0\|_{PM^{(n-1)-p}} \end{align}
(5.7)\begin{align}& \|\mathcal{L} V_0\|_{\mathcal{X}^{2+p}}\leq C \|V_0\|_{PM^{(n-1)-p}}, \end{align}

as well as

(5.8)\begin{align} & \|\mathcal{B}(U,W)\|_{\mathcal{Y}^{-1}}\leq C (\|U\|_{\mathcal{Y}^{-1}} + \|U\|_{\mathcal{X}^3})(\|W\|_{\mathcal{Y}^{-1}} + \|W\|_{\mathcal{X}^3}) \end{align}
(5.9)\begin{align} & \|\mathcal{B}(U,W)\|_{\mathcal{X}^{3}}\leq C (\|U\|_{\mathcal{Y}^{-1}} + \|U\|_{\mathcal{X}^3})(\|W\|_{\mathcal{Y}^{-1}} + \|W\|_{\mathcal{X}^3}) \end{align}
(5.10)\begin{align} \|\mathcal{B}(U,W)\|_{\mathcal{PM}^{(n-1)-p}} & \leq C (\|U\|_{\mathcal{PM}^{(n-1)-p}} + \|U\|_{\mathcal{X}^{2+p}})\nonumber\\ & \quad \times (\|W\|_{\mathcal{PM}^{(n-1)-p}} + \|W\|_{\mathcal{X}^{2+p}}) \end{align}
(5.11)\begin{align}& \|\mathcal{B}(U,W)\|_{\mathcal{X}^{2+p}}\leq C (\|U\|_{\mathcal{PM}^{(n-1)-p}} + \|U\|_{\mathcal{X}^{2+p}})(\|W\|_{\mathcal{PM}^{(n-1)-p}} + \|W\|_{\mathcal{X}^{2+p}}). \end{align}

The Fourier coefficients of $(\mathcal{L} V_0)(t,\cdot)$ are given by

(5.12)\begin{equation} \mathcal{F}(\mathcal{L} V_0)(t,k) = e^{g(t)|k| -t \sigma(\Delta^2+\Delta)(k)}\hat{V_0}(k). \end{equation}

The Fourier coefficients of the nonlinear term are

(5.13)\begin{align} \mathcal{F}&(\mathcal{B} (U,W))(t,k) = -\int_0^t e^{g(t)|k| -(t-s) \sigma(\Delta^2+\Delta)(k)}\Big[ \end{align}
\begin{align*} & \qquad \sum_{j \in \mathbb{Z}_\ast^n, j \neq k} ((k - j) \cdot j) \, e^{-g(s)|k-j|}\hat{U}(s,k-j) e^{-g(s) |j|} \hat{W}(s,j)\Big] \, ds. \end{align*}

In what follows, we will consider two kinds of temporal weights:

(5.14)\begin{align} & g(t) = a \sqrt[4]{t}, \text{for some constant } a \gt 0; \end{align}
(5.15)\begin{align} & g(t) = bt, \text{for some constant } b \gt 0. \end{align}

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:

  1. (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*}
  2. (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

\begin{equation*}\sigma(\Delta^2+\Delta)(k) \geq M_2 |k|^4.\end{equation*}

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$,

\begin{equation*} g(t) |k| - t\sigma(\Delta^2+\Delta)(k) \leq g(t) |k| - t M_2 |k|^4 \leq - \frac{M_2 t}{2}|k|^4, \end{equation*}

as desired. This establishes item (i).

Next, consider the case $g(t) = a\sqrt[4]{t}$. Then, of course, we have

(5.16)\begin{equation} g(t) |k| - t\sigma(\Delta^2+\Delta)(k) \leq a\sqrt[4]{t}|k| - M_2 t |k|^4. \end{equation}

Consider the function $f=f(z)=az - \frac{M_2}{2}z^4$. This function is globally bounded from above. Let

\begin{equation*}C=C(a)=\max\{\sup f(z), 1\}.\end{equation*}

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:

  1. (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$;

  2. (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

\begin{equation*}(g(t) - g(s))|k| - (t-s)\sigma(\Delta^2+\Delta)(k) = (t-s)b|k| - (t-s)\sigma(\Delta^2+\Delta)(k),\end{equation*}

and $t-s \geq 0$.

Next consider $g(t) = a\sqrt[4]{t}$. We first note that

\begin{equation*}(g(t) - g(s)) |k| - (t-s)\sigma(\Delta^2+\Delta)(k) \leq (g(t) - g(s)) |k| - (t-s)M_2 |k|^4,\end{equation*}

using, again, (2.5) and the fact that $t-s\geq 0$. Next, we observe that

\begin{equation*}(g(t) - g(s)) |k| - (t-s)M_2 |k|^4 = f(a\sqrt[4]{t})-f(a\sqrt[4]{s}) - \frac{M_2 (t-s)}{2}|k|^4,\end{equation*}

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:

(5.18)\begin{align} \nonumber |\mathcal{F}&(\mathcal{B} (U,W))(t,k)| = \left|-\int_0^t e^{(g(t)-g(s))|k| -(t-s) \sigma(\Delta^2+\Delta)(k)}\Big[ \right.\\ & \left. \sum_{j \in \mathbb{Z}_\ast^n, j \neq k} ((k - j)\cdot j) \, e^{g(s)(|k|-|k-j|-|j|)}\hat{U}(s,k-j) \hat{W}(s,j)\Big] \, ds \right| \\ & \leq \int_0^t e^{(g(t)-g(s))|k| -(t-s) \sigma(\Delta^2+\Delta)(k)}\Big[ \sum_{\substack{j \in \mathbb{Z}_\ast^n \\ j \neq k}}|k-j||j||\hat{U}(s,k-j) | |\hat{W}(s,j)|\Big] \, ds, \end{align}

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

\begin{equation*}\nonumber c\sum_{k\in\mathbb{Z}^{n}_{*}}e^{(g(t)-\varepsilon)|k|}|\hat{\psi}(t,k)| \leq \sum_{k\in\mathbb{Z}^{n}_{*}} e^{(g(t)-\varepsilon)|k|}\frac{e^{\varepsilon|k|}}{|k|}|\hat{\psi}(t,k)| \leq \|V\|_{\mathcal{Y}^{-1}}. \end{equation*}

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

\begin{equation*}|\hat{f}(k)|=\frac{1}{|k|^{n-1}} \ \text{for all } k.\end{equation*}

Then, $f\in PM^{(n-1)-p},$ but $f\notin Y^{-1}.$ On the other hand, let f be such that

\begin{equation*} |\hat{f}(j)| = \left\{ \begin{array}{ll} |j|^{3/4} & \text{if } |j|= 2^{4\ell}, \; \ell = 1, 2, \ldots \\ 0 & \text{otherwise. } \end{array} \right. \end{equation*}

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

\begin{equation*}\lambda \mapsto \lambda^2 \psi(\lambda x, \lambda^4 t).\end{equation*}

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.

References

Ambrose, D. M., Lopes Filho, M. C. and Nussenzveig Lopes, H. J.. Existence and analyticity of the Lei-Lin solution of the Navier-Stokes equations on the torus. Proc. Amer. Math. Soc 152 (2023), 781795.Google Scholar
Ambrose, D. M. and Mazzucato, A. L.. Global existence and analyticity for the 2D Kuramoto-Sivashinsky equation. J. Dynam. Differ. Equ. 31 (2019), 15251547 10.1007/s10884-018-9656-0.CrossRefGoogle Scholar
Ambrose, D. M. and Mazzucato, A. L.. Global solutions of the two-dimensional Kuramoto-Sivashinsky equation with a linearly growing mode in each direction. J. Nonlinear Sci. 31 (2021), 10.1007/s00332-021-09748-8.CrossRefGoogle Scholar
Auscher, P. and Tchamitchian, P.. Espaces critiques pour le système des equations de Navier-Stokes incompressibles. No modification to the text. This work was done when the first author was at Université de Picardie. https://hal.archives-ouvertes.fr/hal-00344728, (1999).Google Scholar
Bae, H.. Existence and analyticity of Lei-Lin solution to the Navier-Stokes equations. Proc. Amer. Math. Soc. 143 (2015), 28872892 10.1090/S0002-9939-2015-12266-6.CrossRefGoogle Scholar
Bae, H., Biswas, A. and Tadmor, E.. Analyticity and decay estimates of the Navier-Stokes equations in critical Besov spaces. Arch. Ration. Mech. Anal. 205 (2012), 963991 10.1007/s00205-012-0532-5.CrossRefGoogle Scholar
Bellout, H., Benachour, S. and Titi, E. S.. Finite-time singularity versus global regularity for hyper-viscous Hamilton-Jacobi-like equations. Nonlinearity. 16 (2003), 19671989 10.1088/0951-7715/16/6/305.CrossRefGoogle Scholar
Benachour, S., Kukavica, I., Rusin, W. and Ziane, M.. Anisotropic estimates for the two-dimensional Kuramoto-Sivashinsky equation. J. Dynam. Differ. Equ. 26 (2014), 461476 10.1007/s10884-014-9372-3.CrossRefGoogle Scholar
Biswas, A. and Swanson, D.. Existence and generalized Gevrey regularity of solutions to the Kuramoto-Sivashinsky equation in $\mathbb{R}^n$. J. Differ. Equ. 240 (2007), 145163 10.1016/j.jde.2007.05.022.CrossRefGoogle Scholar
Bronski, J. C. and Gambill, T. N.. Uncertainty estimates and L 2 bounds for the Kuramoto-Sivashinsky equation. Nonlinearity. 19 (2006), 20232039 10.1088/0951-7715/19/9/002.CrossRefGoogle Scholar
Cannone, M.. Ondelettes, Paraproduits et Navier-Stokes. (Diderot Editeur, Paris, 1995). With a preface by Yves Meyer.Google Scholar
Cannone, M.. Harmonic analysis and Navier-Stokes equations with application to the Boltzmann equation. Lectures on the Analysis of Nonlinear Partial Differential equations. Part 2, Morningside Lect. Math. Vol. 2, , (Int. Press, Somerville, MA, 2012).Google Scholar
Cannone, M. and Karch, G.. Smooth or singular solutions to the Navier-Stokes system? J. Differ. Equ. 197 (2004), 247274 10.1016/j.jde.2003.10.003.CrossRefGoogle Scholar
Coti Zelati, M., Dolce, M., Feng, Y. and Mazzucato, A. L.. Global existence for the two-dimensional Kuramoto-Sivashinsky equation with a shear flow. J. Evol. Equ. 21 (2021), 50795099 10.1007/s00028-021-00752-9.CrossRefGoogle Scholar
Doelman, A. and Titi, E. S.. Regularity of solutions and the convergence of the Galerkin method in the Ginzburg-Landau equation. Numer. Funct. Anal. Optim. 14 (1993), 299321 10.1080/01630569308816523.CrossRefGoogle Scholar
Feng, Y. and Mazzucato, A. L.. Global existence for the two-dimensional Kuramoto-Sivashinsky equation with advection. Comm. Partial Differ. Equ. 47 (2022), 279306 10.1080/03605302.2021.1975131.CrossRefGoogle Scholar
Ferrari, A. B. and Titi, E. S.. Gevrey regularity for nonlinear analytic parabolic equations. Comm. Partial Differ. Equ. 23 (1998), 116 10.1080/03605309808821336.CrossRefGoogle Scholar
Foias, C. and Temam, R.. Gevrey class regularity for the solutions of the Navier-Stokes equations. J. Funct. Anal. 87 (1989), 359369 10.1016/0022-1236(89)90015-3.CrossRefGoogle Scholar
Goodman, J.. Stability of the Kuramoto-Sivashinsky and related systems. Comm. Pure Appl. Math. 47 (1994), 293306 10.1002/cpa.3160470304.CrossRefGoogle Scholar
Grujić, Z. and Kukavica, I.. Space analyticity for the Navier-Stokes and related equations with initial data in Lp. J. Funct. Anal. 152 (1998), 447466 10.1006/jfan.1997.3167.CrossRefGoogle Scholar
Kalogirou, A., Keaveny, E. E. and Papageorgiou, D. T.. An in-depth numerical study of the two-dimensional Kuramoto-Sivashinsky equation. Proc. A. 471 (2015), 10.1098/rspa.2014.0932. 20140932.Google ScholarPubMed
Koch, H. and Tataru, D.. Well-posedness for the Navier-Stokes equations. Adv. Math. 157 (2001), 2235 10.1006/aima.2000.1937.CrossRefGoogle Scholar
Kukavica, I. and Massatt, D.. On the global existence for the Kuramoto-Sivashinsky equation. J. Dynam. Differ. Equ. 35 (2023), 6985 10.1007/s10884-021-09985-1.CrossRefGoogle Scholar
Kuramoto, Y. and Tsuzuki, T.. Persistent propagation of concentration waves in dissipative media far from thermal equilibrium. Prog. Theor. Phys. 55 (1976), 356369.CrossRefGoogle Scholar
Larios, A. and Yamazaki, K.. On the well-posedness of an anisotropically-reduced two-dimensional Kuramoto-Sivashinsky equation. Phys. D. 411 (2020), 10.1016/j.physd.2020.132560.CrossRefGoogle Scholar
Lei, Z. and Lin, F.. Global mild solutions of Navier-Stokes equations. Comm. Pure Appl. Math. 64 (2011), 12971304 10.1002/cpa.20361.CrossRefGoogle Scholar
Miao, C. and Yuan, B.. Solutions to some nonlinear parabolic equations in pseudomeasure spaces. Math. Nachr. 280 (2007), 171186 10.1002/mana.200410472.CrossRefGoogle Scholar
Molinet, L.. A bounded global absorbing set for the Burgers-Sivashinsky equation in space dimension two. C. R. Acad. Sci. Paris Sér. I Math. 330 (2000), 635640 10.1016/S0764-4442(00)00224-X.CrossRefGoogle Scholar
Molinet, L.. Local dissipativity in L 2 for the Kuramoto-Sivashinsky equation in spatial dimension 2. J. Dynam. Differ. Equ. 12 (2000), 533556 10.1023/A:1026459527446.CrossRefGoogle Scholar
Nicolaenko, B., Scheurer, B. and Temam, R.. Some global dynamical properties of the Kuramoto-Sivashinsky equations: nonlinear stability and attractors. Phys. D. 16 (1985), 155183 10.1016/0167-2789(85)90056-9.CrossRefGoogle Scholar
Papageorgiou, D. T., Smyrlis, Y. -S. and Tomlin, R. J.. Optimal analyticity estimates for non-linear active-dissipative evolution equations. IMA J. Appl. Math. 87 (2022), 964984 10.1093/imamat/hxac028.CrossRefGoogle Scholar
Reed, M. and Simon, B.. Methods of Modern Mathematical Physics. II. Fourier analysis, self-adjointness. (Academic Press [Harcourt Brace Jovanovich, Publishers], New York-London, 1975).Google Scholar
Sell, G. R. and Taboada, M.. Local dissipativity and attractors for the Kuramoto-Sivashinsky equation in thin $2{{\rm} D}$ domains. Nonlinear Anal. 18 (1992), 671687 10.1016/0362-546X(92)90006-Z.CrossRefGoogle Scholar
Sivashinsky, G. I.. Nonlinear analysis of hydrodynamic instability in laminar flames–I. Derivation of basic equations. Acta Astronaut. 4 (1977), 11771206.CrossRefGoogle Scholar
Tadmor, E.. The well-posedness of the Kuramoto-Sivashinsky equation. SIAM J. Math. Anal. 17 (1986), 884893 10.1137/0517063.CrossRefGoogle Scholar