Ancient Solutions to Navier–Stokes Equations in Half Space

In the present paper we consider mild bounded ancient (backward) solutions to the Navier–Stokes equations in the half plane. We give two different definitions, prove their equivalence and prove smoothness up to the boundary. Such solutions appear as a result of rescaling around a singular point of the initial boundary value problem for the Navier–Stokes equations in the half-plane.


Introduction
The goal of the paper is to understand properties of the so-called ancient (backward) solutions to the Navier-Stokes equations. The importance of them in the regularity theory for the Navier-Stokes equations, see, for example, papers [1,5,12,15], and [16], and more generally in the theory of PDEs is well understood. They appear as a limit, resulting from rescaling solutions to initial boundary value problems around possible singularities. For the Navier-Stokes equations, this procedure has been described in the above papers.
The weakest version of ancient solutions to the Navier-Stokes equations is as follows. for any q ∈ C ∞ 0 (Q − ). This class of ancient solutions seems to be too wide. Having in mind the problem of regularity for solutions to the Navier-Stokes equations mentioned above, we can put some additional restrictions in the definition of ancient solutions. If an ancient solution u is bounded, we call it a bounded ancient one. We can go further, see [5], and consider an even more narrow class of ancient solutions. We say that a bounded function u is a mild bounded ancient solution if u has the following property: for any A < 0 and for (x, t) ∈ Q A := R 3 ×]A, 0[, This statement seems to be known and we give its prove for completeness. One of the interesting consequences of the above proposition is an alternative proof of smoothness of mild bounded ancient solutions, see [13].
The conjecture that has been made in [5] reads: any mild bounded ancient solution is a constant in Q − . The validity of this conjecture is known in several cases, see details in [5] and [15]. The connection with a possible blowup of a solution to the initial value problem , is as follows. Assume that there is a blowup at t = T , i.e., v(·, t) ∞,R 3 → ∞ as t → T − . Then there exists a mild bounded ancient solution u with |u(0)| = 1. If the aforesaid conjecture is true then u(x, t) = c, where c is a constant vector such that |c| = 1. This would rule out blowups of Type I for which a certain scale-invariant quantity is bounded. Now, let us formulate the main results of the paper about mild bounded ancient solutions to the Navier-Stokes equations in half space, starting with a definition of distributional ancient solutions. From now on we denote Q + − := R 3 + ×] − ∞, 0[. We say that u ∈ L 2 (B + (R)) for any R > 0 is an ancient solution if u satisfies Q + − u · (∂ t ϕ + Δϕ) + u ⊗ u : ∇ϕ dxdt = 0 (1.5) for any ϕ ∈ C ∞ 0,0 (Q − ) with ϕ(x , 0, t) = 0 for any x ∈ R 2 and for any t < 0. Moreover u satisfies for any q ∈ C ∞ 0 (Q − ). We can notice that (1.5) and (1.6) is a weak form for the following ∂ t u + u · ∇u − Δu = −∇p, div u = 0 in Q + − for some distribution p, u(x , 0, t) = 0 for any x ∈ R 2 and any −∞ < t < 0. We shall say that u is a bounded ancient solution in half space if it is ancient and bounded.
Vol. 17 (2015) Ancient Solutions to Navier-Stokes Equations in Half Space 553 From now on define Q + A := R 3 + ×]A, 0[. In order to proceed further, we need to recall how one can construct a solution to the following boundary value for the Stokes equations in half space: It is assumed that f and u 0 are divergence free and f 3 (x , 0, t) = 0. Then a formal solution to the above initial boundary value problem is: The Green's function G has been derived by Solonnikov in [18] and is as follows Let us introduce another potential K = (K mjs ), where the tensor Φ = (Φ ij ) are defined as solutions to the following boundary value problems with ∂Φ mn /∂y 3 (x, y, t) = 0 if n < 3 and with Φ mn (x, y, t) = 0 if n = 3 at y 3 = 0. Now, we are in position to define mild bounded ancient solutions in a half space.

Definition 1.2.
A bounded divergence free function u in Q + − is called a mild bounded ancient solution if, for any A < 0 and any (x, t) ∈ Q + A , To state our main result, we need to introduce the following operator. Given H = (H ij ) ∈ L ∞ (R 3 + ), there exists a unique function p 1 ∈ L 2 (B + (R)) for any R > 0 with [p 1 ] B+ = 0 with the following properties: the even extension of it to R 3 belongs to the space BM O, Barker and G. Seregin JMFM for any ϕ ∈ C ∞ 0 (R 3 + ) with ϕ ,3 (x , 0) = 0 and where A is an absolute constant. We shall use notation p 1 := p 1 H . We notice that if H = (H ij ) is sufficiently smooth and vanishes on the boundary x 3 = 0, the function p 1 H is a solution to the Neumann boundary value problem:

Theorem 1.4. A bounded function u is a mild bounded ancient solution if and only if there exists a pressure p such that
as x 3 → ∞ and for any t < 0; u and p satisfy (1.6) and for any ϕ ∈ C ∞ 0 (Q − ) with ϕ(x , 0, t) = 0 for any x ∈ R 2 and for any t < 0. In [16], there has been conjectured that any mild bounded ancient solution is identically equal to zero in Q + − . At the moment of writing the paper, there are two cases in which the above conjecture is true, see [2] and [14]. Both cases are two-dimensional and additional scale-invariant assumptions have been imposed. In the first paper vorticity preserves its sign, while in the second one kinetic energy is bounded. Now, let us consider the following initial boundary value problem Then there exists a sequence z n = (x (n) , t n ) such that t n > 0, t n → T − , and In conclusion, we notice that the validity of both conjectures allow us to rule out at least Type I blowups of solution to initial boundary value problem for the Navier-Stokes equations in half space.

Proof of Theorem 1.3
Before starting the proof, we remind known facts (due to Solonnikov, see [18,19]), about the Green function and the kernel K.
It is not so difficult to see that the kernel K has the structure is a linear combination of the terms and K ism (x, z, t) is a linear combination of the terms Here, The following estimates for G i and K have been obtained in different papers of Solonnikov, see [18,19]: where α = (α 1 , α 2 ), γ = (γ 1 , γ 2 ), and |γ| = 0 or 1, for l = 0 or 1. Let K 1 and K 2 be generated by G 1 and G 2 , respectively. In particular, we have the estimate In what follows, we are going to use special approximations: and T. Barker and G. Seregin JMFM Here, u , has the additional property that the bounds of D α φ k only depend on |α|. The standard mollifiers are denoted by η and ω, respectively. The properties of the approximation scheme are that u (k) ∈ C ∞ 0 (]A − 2, 1[×R 3 + ) and that (up to subsequence) . It is noticed that u (k)A is a smooth solenodial vector field, with bounded derivatives that all vanish near x 3 = 0. Furthermore, . We let F k := u (k) ⊗ u (k) and then It is not so difficult to infer that (up to subsequence): To treat the second term on the right hand side of representation formula, we are going to use the following statement.
Under assumptions imposed on tensor-valued function F , respectively. In order to construct such a sequence, we proceed as follows. Let ) and ϕ ≡ 1 in B. And finally we can produce F m using (F (h,R) ) with 0 < < h, where (g) is a mollification of the function g.
We also can state that p 1 F m has decay 1 |x| 2 as |x| → ∞. So, we do not need to take care of integrability of functions involved because of Solonnikov estimates and the decay of the pressure. Similarly the decay of the pressure allows one to rigorously justify the integration by parts shown below. This is done by proving slow decay of the kernels, using arguments in [19,20]. Vol. 17 (2015) Ancient Solutions to Navier-Stokes Equations in Half Space 557 Now, letting f m = −divF m − ∇p 1 F m , we have By our assumptions on boundary values of functions Φ and p 1 F m and their derivatives, the integral over the plane x 3 = 0 vanishes. So, we have For the same reason, the surface integral is equal to zero and using the pressure equation, we find So, the formula of the lemma proved for F m , i.e., we have Now, the identity of the lemma can be obtained by passage to the limits in the latter identity as m → ∞.
We then put U (k) = U 1 (k) + U 2 (k) according to splitting of the kernel described in (1.7). Furthermore, . Where for m = 1, 2: T. Barker and G. Seregin JMFM and Notice, that we may integrate by parts in (2.11) to get: This is permissible by the facts that G 2 i3 (x, y, t) = 0, G 1 ij (x, y, t) = 0 (on y 3 = 0), the spatial decay of G 1 and G 2 and that the approximation scheme implies that p 1(k) (y, τ ) has spatial decay of order |y| −2 .
Now we proceed in proving the main body of Theorem 1.3. Most of the proof is split into four main Propositions. The first two Propositions are derived from arguments from [16]. However, we provide some adjustments, simplifications and demonstrate how those arguments interact with the aforementioned approximation scheme.
In what follows, we are going to use additional notation. For p and q between 1 and infinity we say . From now on we use the terms even and odd extensions to mean the following. For for y 3 0 and f even (y) := f (y * ) for y 3 < 0. This is referred to as the even extension of f . The odd extension of f is similarly defined.
. Notice that by classical singular integral theory, we get for the even extension of the pressure: (2.12) By the proerties of the heat kernel and estimates (2.2) it is obtained that From (2.12), along with arguments in [16] (see Lemma 6.1 there), obtain also that (for (x, t) ∈ Q + A ): (2.14) s, A). By the definition of G 1 , the following is + is arbitrary and [p 1(k) ] B(z,2) (t) signifies the average over B(z, 2) for the even extension. Using (2.12) and local regularity theory for heat equation (e.g Appendix of [9]), obtain: .
Thus, using additional properties of the heat kernel and the properties of the approximations u (k) and Hence, we obtain a weak derivative formula for classes of test functions not necessarily zero on C(d, s, A) and that (2.16) holds also for u (note one can Next, fix δ < 0, with |δ| small and let k be sufficiently large such that δ + 1 k < 0. Observing the structure of the approximations u (k) , the analogue of (2.16) (with u) gives for (x, t) ∈ Q + A,δ : Thus the improvement of u gives (for sufficiently large k K(δ)): Now p 1(k) satisfies (for appropriate even and odd extensions of p 1(k) and u i(k) u j(k) ) in Q A . Local regularity theory for Laplace equation gives: To estimate ∇U 2,2 (k) we need the following statement, whose proof is contained in the Appendix.   s, A). The conclusion regarding boundedness of the gradient of u in Q + − is inferred from taking limits and time-shift arguments. The statement regarding p 1

L∞(BMO)
is deduced from (2.12). It remains to prove sup (x,t)∈Q + − |∇p 1 (x, t)| C( u L∞(Q + − ) ). Notice from (1.6) that for (x, t) ∈ Q + A,δ : From the latter, it follows that div u (k) and ∇div u (k) are bounded function in space-time and in k. Local regularity for Laplace equation gives (for the even extension of the pressure): The conclusion is then reached by arguments similar to those previously mentioned. Proposition 2.3 is proven.
Before proceeding the proof of Theorem 1.3, let us adopt the notation: Proposition 2.6. Suppose all the assumptions of Proposition 2.3 hold. Then u also satisfies: Proof of Proposition 2.6. First notice that, due to estimates (2.2), there exists a positive constant C (independent of x ∈ R 3 + and > 0) such that: (2.28) Note that properties of the approximation scheme imply F k jl and p 1(k) are smooth on R 3 + ×]A, 0[ with bounded derivatives (the bounds may be dependent on k). Change variables to get: It follows that differentiation in time is permissible and gives: Notice we have (for > 0): Thus using (2.28) and estimates (2.2)-(2.3) obtain:

c(d, s, A, γ).
Using also (2.25), one obtains higher regularity for (2.15) through local regularity results for the heat equation. A parabolic imbedding theorem then gives (2.32). Proposition 2.6 is proven.
The next Proposition will improve the previously obtained regularity results. But, first let us state a lemma, which is a simplified version of a more general statement proven in the Appendix, see Lemma 5.1.

Lemma 2.7.
Suppose the measurable kernel K : R n \{0} → R satisfies the conditions (see [21]): For suitable f define the singular integral operator: Take a compactly supported function g in L p (R n ), where 1 < p < ∞. Furthermore assume g is in L ∞ (R n 1 ). Then it follows that (almost everywhere) and h 2 L p,unif (R n ) C(n, p, B) g L p,unif (R n ) . Proof. Fix δ < 0, with |δ| small. From Proposition 2.6, we have (for k K(δ) sufficiently large): It is clear that ∂ t p 1(k) = p 1 ∂t(u (k) ⊗u (k) ) and furthermore Lemma 2.7 is applicable to a suitable extension of ∂ t (u (k) ⊗ u (k) ) for n = 3. Hence, it can be written that ∂ t p 1(k) = (∂ t p 1(k) ) 1 + (∂ t p 1(k) ) 2 . Here, Next, it is easy to see that the following is satisfied in D (Q + A,δ ): Here, z ∈ R 3 + is arbitrary. Hence, (2.31), (2.42), (2.44), (2.45 ), and local regularity of heat equation (for sufficiently large s > n + 2) gives: Now, let us examine ∂ t U 2,2 i(k) . One can write: and A. For the second term, decompose as follows: ((x ,0)

Proposition 2.3 implies that the first term is bounded on Q +
for any k= 0, 1. . .
Proof of Proposition 2.9. We give a brief account of the bootstrap arguments. Clearly ∂ t u also satisfies (1.6). By properties of the kernel, if ϕ ∈ C ∞ 0 (Q A ) then the following holds: From Proposition 2.8, we can write: Here, So using Proposition 2.3 along with Green function estimates (2.2) and (2.4) get that: any k, l = 0, 1. . . The third term of (2.53) is dealt with by splitting the integral according to the kernel decomposition (1.7). The arguments are now repeated from Propositions 2.3-2.8. It is possible to repeat this argument indefinitely with higher time derivatives. Proposition 2.9 is proven.
Vol. 17 (2015) Ancient Solutions to Navier-Stokes Equations in Half Space 565 Now, one can recover a pressure p such that in Q + − : for any x ∈ R 2 and any −∞ < t < 0. By considering higher derivatives in time of these equations and Proposition 2.9, one can obtain Theorem 1.3 using the regularity theory of the stationary Stokes system together with bootstrap arguments.  x0,1) )∇ϕ and x0,1) )I)).

Proof of
and v 1 = 0 on ∂ Q + (z 0 , 1). By our assumptions, for any 1 < s < ∞. Therefore we can claim that |∇v 1 | ≤ c on Q + (z 0 , 3/4) with a constant independent of z 0 . Notice that by our assumptions we may write f 2 = div(F 2 ), where So using boundary regularity theory for the heat equation, we can say that ∇v 2 s,Q(z0,3/4) ≤ c. Then we can see that, since with Neumann boundary condition on the flat part of the boundary, we infer that on Q + (z 0 , 1/2) with a constant independent of z 0 . So, we have boundedness near the boundary. To get the interior estimate, we can use the same arguments that used for mild bounded ancient solutions in the whole space. In fact, we have even more ∇u is continuous up the boundary and for any s > 1. Lemma 3.1 is proven.
We wish to show that u has the following properties. For any A < 0, Let us go back to Lemma 2.1 and its proof. Take F and its approximations F m are from that lemma and its proof. Solonnikov showed in [19] that v m given by the formula satisfies the identity: for any divergence free functions ϕ ∈ C ∞ 0 (Q − ) with ϕ(x , 0, t) = 0 for any x ∈ R 2 and for any t < 0. Since we can use boundedness and pass to the limit as m → ∞.
As a result, we have Vol. 17 (2015) Ancient Solutions to Navier-Stokes Equations in Half Space 567 for any q ∈ C ∞ 0 (Q A ), where Q A := R 3 ×]A, 0[. Here, v is defined as u 1 + u 2 . Furthermore, for any divergence free functions ϕ ∈ C ∞ 0 (Q − ) with ϕ(x , 0, t) = 0 for any x ∈ R 2 and for any t < 0: Now assume that F = u ⊗ u. From the the previous pages, it is clear that u is continuous in the completion of Q + (R) := B + (R)×] − R 2 , 0[. Then using cut-off functions in time, we can show that u satisfies the same identity as v. And thus letting w = u − v, we get for any divergence free functions ϕ ∈ C ∞ 0 (Q − ) with ϕ(x , 0, t) = 0 for any x ∈ R 2 and for any t < 0. If we extend w by zero for t ≤ A, we find for any divergence free functions ϕ ∈ C ∞ 0 (Q − ) with ϕ(x , 0, t) = 0 for any x ∈ R 2 and for any t < 0. By the Liouville theorem (see [3,4]), w = w(x 3 , t). We need to show that w ≡ 0 in Q + A . To this end, it is sufficient to show that for x ∈ R 2 and t ∈]A, 0[, one has ∇q 2 (x , x 3 , t) → 0 as x 3 → ∞. Here, q 2 is the pressure for v so that for x ∈ R 2 and A < t < 0, and v 1 (·, 0) = u 0 (·). Thus, On the other hand, T. Barker and G. Seregin JMFM and Using reasoning from Proposition 2.6, obtain Regarding v 2,2 , we have So, we have the required decay for ∇q 2 and thus we have for all A < 0 the following integral representation: for all t > A. Now, our aim is to prove the inverse statement, i.e., we assume that bounded divergence free function satisfied the latter identity for any A < 0. Introducing F = u ⊗ u and tensor H = F + p 1 I and using approximations of u, we can show that u is a distributional solution to the Navier-Stokes equations in Q + − and belongs to the space W 1 ∞ (Q + − ). This can be done in the same way as in [16] (the most difficult part of that paper). We then can introduce the pressure p 2 so that Vol. 17 (2015) Ancient Solutions to Navier-Stokes Equations in Half Space 569 Splitting u = u 1 + u 2 and repeating the aforesaid arguments, we can show that ∇p 2 satisfies all requirements in the definition of bounded mild ancient solutions.

Proof of Proposition 1.1
Assume first that u satisfies the conditions of Proposition 1.1, i.e., there exists a pressure field p ∈ Our aim is to show that, for any A < 0, the function satisfies the integral identity (1.3). First, let us notice that the pressure p (up to a bounded function of time t) is formally represented as follows: We know that mild bounded ancient solutions are infinitely smooth and all partial derivatives, apart from derivatives in time of the pressure, are bounded. The derivatives ∂ k t p, k = 0, 1, . . ., belong to L ∞ (−∞, 0; BM O(R 3 ), see [13]. So, we re-write the Navier-Stokes equations in the following way: We know that f is infinitely smooth and all its derivative are bounded. Then, by Tychonoff's uniqueness theorem, for t > A and for all A < 0. It remains to show that for any A < 0. To this end, we introduce as notation p F which is the BMO-solution to to equation We deduce the required identity from the following lemma.
The proof can be done with the help of suitable approximation of F and the following estimates: The first one is due to Solonnikov, see [17]), and the second one is well known. Inverse statement of Proposition 1.1 can be easily deduced from the above lemma and suitable approximations of u. This completes the proof of Proposition 1.1.
Open Access. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. By (2.22), l > 2 and thus l < 2. So, the last factor is finite. Hence, . Similar arguments to before give, for m 1, . Summing over m and R we then conclude.
Suppose the measurable kernel K : R n × R n \Δ → R is such that there exists M > 0 with: Define the truncation (on L p (R n ), 1 < p < ∞): Suppose, for this kernel, there exists bounded linear operator T : : Furthermore for f ∈ L ∞ (R n ) compactly supported: Here, C(K, n) means that the constant depends on the properties of the Kernel (e.g some smoothness of the kernel) and the dimension of the space. Consider an unbounded domain Ω ⊂ R n that is contained between two n − 1 dimensional parallel hyperplanes (denoted Π 1 and Π 2 respectively) a finite distance 2L apart. Take g to be a compactly supported function in L p (R n ) (for 1 < p < ∞) such that g is non-zero and bounded outside of Ω.
Then it follows that: T g(x) = h 1 (x) + h 2 (x). (5.8) Here, h 1 BMO(R n ) c(K, n) g L∞(R n \Ω) (5.9) and h 2 L p,unif (R n ) c(K, M, n, p, L) g L p,unif (R n ) . (5.10) Proof of Lemma 5.1. For x in R n−1 , denote by Q R n−1 (x , R) a cube in R n−1 with side lengths 2R centred at x . First one shows that, without loss of generality, it is sufficient to reduce to the case where: By rotation invariance of L p (R n ) and BM O(R n ), it is inherited from T and T that there exists a bounded linear operator S : L p (R n ) → L p (R n ) (where 1 < p < ∞) such that (5.5)-(5.7) hold. Since the space L p,unif (R n ) is rotation invariant, rotations of Ω can be considered without loss of generality. Fix z 1 ∈ R n . It can be inferred that: LetK (x, y) := K(x − z 1 , y − z 1 ). (5.17) Using the spaces BM O(R n ), L p (R n ) and L p,unif (R n ) are translation invariant, one can use the aforementioned arguments to show that translations of Ω can be considered without loss of generality. From now on take Π 1 as in (5.11) and Π 2 as in (5.12). Now decompose g: So it is obtained that: From here all conclusions follow immediately.